Reverse engineering of regulatory networks in human B cells

Cellular phenotypes are determined by the differential activity of networks linking coregulated genes. Available methods for the reverse engineering of such networks from genome-wide expression profiles have been successful only in the analysis of lower eukaryotes with simple genomes. Using a new method called ARACNe (algorithm for the reconstruction of accurate cellular networks), we report the reconstruction of regulatory networks from expression profiles of human B cells. The results are suggestive a hierarchical, scale-free network, where a few highly interconnected genes (hubs) account for most of the interactions. Validation of the network against available data led to the identification of MYC as a major hub, which controls a network comprising known target genes as well as new ones, which were biochemically validated. The newly identified MYC targets include some major hubs. This approach can be generally useful for the analysis of normal and pathologic networks in mammalian cells.     

An epi-allelic series of p53 hypomorphs created by stable RNAi produces distinct tumor phenotypes in vivo

The application of RNA interference (RNAi) to mammalian systems has the potential to revolutionize genetics and produce novel therapies. Here we investigate whether RNAi applied to a well-characterized gene can stably suppress gene expression in hematopoietic stem cells and produce detectable phenotypes in mice. Deletion of the Trp53 tumor suppressor gene greatly accelerates Myc-induced lymphomagenesis, resulting in highly disseminated disease. To determine whether RNAi suppression of Trp53 could produce a similar phenotype, we introduced several Trp53 short hairpin RNAs (shRNAs) into hematopoietic stem cells derived from E(mu)-Myc transgenic mice, and monitored tumor onset and overall pathology in lethally irradiated recipients. Different Trp53 shRNAs produced distinct phenotypes in vivo, ranging from benign lymphoid hyperplasias to highly disseminated lymphomas that paralleled Trp53-/- lymphomagenesis in the E(mu)-Myc mouse. In all cases, the severity and type of disease correlated with the extent to which specific shRNAs inhibited p53 activity. Therefore, RNAi can stably suppress gene expression in stem cells and reconstituted organs derived from those cells. In addition, intrinsic differences between individual shRNA expression vectors targeting the same gene can be used to create an 'epi-allelic series' for dissecting gene function in vivo.     

High frequency of BRAF mutations in nevi

To evaluate the timing of mutations in BRAF (v-raf murine sarcoma viral oncogene homolog B1) during melanocytic neoplasia, we carried out mutation analysis on microdissected melanoma and nevi samples. We observed mutations resulting in the V599E amino-acid substitution in 41 of 60 (68%) melanoma metastases, 4 of 5 (80%) primary melanomas and, unexpectedly, in 63 of 77 (82%) nevi. These data suggest that mutational activation of the RAS/RAF/MAPK pathway in nevi is a critical step in the initiation of melanocytic neoplasia but alone is insufficient for melanoma tumorigenesis.     

Postzygotic HRAS and KRAS mutations cause nevus sebaceous and Schimmelpenning syndrome

Nevus sebaceous is a common congenital cutaneous malformation. Affected individuals may develop benign and malignant secondary tumors in the nevi during life. Schimmelpenning syndrome is characterized by the association of nevus sebaceous with extracutaneous abnormalities. We report that of 65 sebaceous nevi studied, 62 (95%) had mutations in the HRAS gene and 3 (5%) had mutations in the KRAS gene. The HRAS c.37G>C mutation, which results in a p.Gly13Arg substitution, was present in 91% of lesions. Nonlesional tissues from 18 individuals had a wild-type sequence, confirming genetic mosaicism. The HRAS c.37G>C mutation was also found in 8 of 8 associated secondary tumors. Mosaicism for HRAS c.37G>C and KRAS c.35G>A mutations was found in two individuals with Schimmelpenning syndrome. Functional analysis of HRAS c.37G>C mutant cells showed constitutive activation of the MAPK and PI3K-Akt signaling pathways. Our results indicate that nevus sebaceous and Schimmelpenning syndrome are caused by postzygotic HRAS and KRAS mutations. These mutations may predispose individuals to the development of secondary tumors in nevus sebaceous.     

Deregulated expression of c-Myc depletes epidermal stem cells

The beta-catenin/TCF signaling pathway is essential for the maintenance of epithelial stem cells in the small intestine. c-Myc a downstream target of beta-catenin/TCF (ref. 2), can induce differentiation of epidermal stem cells in vitro. To determine the role of c-Myc in epidermal stem cells in vivo, we have targeted expression of human MYC2 to the hair follicles and the basal layer of mouse epidermis using a keratin 14 vector (K14.MYC2). Adult K14.MYC2 mice gradually lose their hair and develop spontaneous ulcerated lesions due to a severe impairment in wound healing; their keratinocytes show impaired migration in response to wounding. The expression of beta1 integrin, which is preferentially expressed in epidermal stem cells is unusually low in the epidermis of K14.MYC2 mice. Label-retaining analysis to identify epidermal stem cells reveals a 75% reduction in the number of stem cells in 3-month-old K14.MYC2 mice, compared with wildtype mice. We conclude that deregulated expression of c-Myc in stem cells reduces beta1 integrin expression, which is essential to both keratinocyte migration and stem cell maintenance.     

Senescence bypass screen identifies TBX2, which represses Cdkn2a (p19(ARF)) and is amplified in a subset of human breast cancers

To identify new immortalizing genes with potential roles in tumorigenesis, we performed a genetic screen aimed to bypass the rapid and tight senescence arrest of primary fibroblasts deficient for the oncogene Bmi1. We identified the T-box member TBX2 as a potent immortalizing gene that acts by downregulating Cdkn2a (p19(ARF)). TBX2 represses the Cdkn2a (p19(ARF)) promoter and attenuates E2F1, Myc or HRAS-mediated induction of Cdkn2a (p19(ARF)). We found TBX2 to be amplified in a subset of primary human breast cancers, indicating that it might contribute to breast cancer development.     

Trans-dominant inhibition of RNA viral replication can slow growth of drug-resistant viruses

The high error rates of viral RNA-dependent RNA polymerases create heterogeneous viral populations whose disparate RNA genomes affect each other's survival. We systematically screened the poliovirus genome and identified four sets of dominant mutations. Mutated alleles in capsid- and polymerase-coding regions resulted in dominant negative phenotypes, probably due to the proteins' oligomeric properties. We also identified dominant mutations in an RNA element required for priming RNA synthesis (CRE) and in the protein primer (VPg), suggesting that nonproductive priming intermediates are inhibitory. Mutations that inhibit the activity of viral proteinase 2A were dominant, arguing that inhibition of its known intramolecular activity creates a toxic product. Viral products that, when defective, dominantly interfere with growth of nondefective viruses will probably be excellent drug targets because drug-sensitive viruses should be dominant over drug-resistant variants. Accordingly, a virus sensitive to anticapsid compound WIN51711 dominantly inhibited the intracellular growth of a drug-resistant virus. Therefore, dominant inhibitor screening should validate or predict targets for antiviral therapy with reduced risk for drug resistance.     

Chromatin profiling using targeted DNA adenine methyltransferase

Chromatin is the highly complex structure consisting of DNA and hundreds of associated proteins. Most chromatin proteins exert their regulatory and structural functions by binding to specific chromosomal loci. Knowledge of the identity of these in vivo target loci is essential for the understanding of the functions and mechanisms of action of chromatin proteins. We report here large-scale mapping of in vivo binding sites of chromatin proteins, using a novel approach based on a combination of targeted DNA methylation and microarray technology. We show that three distinct chromatin proteins in Drosophila melanogaster cells each associate with specific sets of genes. HP1 binds predominantly to pericentric genes and transposable elements. GAGA factor associates with euchromatic genes that are enriched in (GA)n motifs. A Drosophila homolog of Saccharomyces cerevisiae Sir2p is associated with several active genes and is excluded from heterochromatin. High-resolution, genome-wide maps of target loci of chromatin proteins ('chromatin profiles') provide new insights into chromatin structure and gene regulation.     

Common genetic variants account for differences in gene expression among ethnic groups

Variation in DNA sequence contributes to individual differences in quantitative traits, but in humans the specific sequence variants are known for very few traits. We characterized variation in gene expression in cells from individuals belonging to three major population groups. This quantitative phenotype differs significantly between European-derived and Asian-derived populations for 1,097 of 4,197 genes tested. For the phenotypes with the strongest evidence of cis determinants, most of the variation is due to allele frequency differences at cis-linked regulators. The results show that specific genetic variation among populations contributes appreciably to differences in gene expression phenotypes. Populations differ in prevalence of many complex genetic diseases, such as diabetes and cardiovascular disease. As some of these are probably influenced by the level of gene expression, our results suggest that allele frequency differences at regulatory polymorphisms also account for some population differences in prevalence of complex diseases.     

Integrated transcriptional profiling and linkage analysis for identification of genes underlying disease

Integration of genome-wide expression profiling with linkage analysis is a new approach to identifying genes underlying complex traits. We applied this approach to the regulation of gene expression in the BXH/HXB panel of rat recombinant inbred strains, one of the largest available rodent recombinant inbred panels and a leading resource for genetic analysis of the highly prevalent metabolic syndrome. In two tissues important to the pathogenesis of the metabolic syndrome, we mapped cis- and trans-regulatory control elements for expression of thousands of genes across the genome. Many of the most highly linked expression quantitative trait loci are regulated in cis, are inherited essentially as monogenic traits and are good candidate genes for previously mapped physiological quantitative trait loci in the rat. By comparative mapping we generated a data set of 73 candidate genes for hypertension that merit testing in human populations. Mining of this publicly available data set is expected to lead to new insights into the genes and regulatory pathways underlying the extensive range of metabolic and cardiovascular disease phenotypes that segregate in these recombinant inbred strains.     

A conserved microRNA module exerts homeotic control over Petunia hybrida and Antirrhinum majus floral organ identity

It is commonly thought that deep phylogenetic conservation of plant microRNAs (miRNAs) and their targets indicates conserved regulatory functions. We show that the blind (bl) mutant of Petunia hybrida and the fistulata (fis) mutant of Antirrhinum majus, which have similar homeotic phenotypes, are recessive alleles of two homologous miRNA-encoding genes. The BL and FIS genes control the spatial restriction of homeotic class C genes to the inner floral whorls, but their ubiquitous early floral expression patterns are in contradiction with a potential role in patterning C gene expression. We provide genetic evidence for the unexpected function of the MIRFIS and MIRBL genes in the center of the flower and propose a dynamic mechanism underlying their regulatory role. Notably, Arabidopsis thaliana, a more distantly related species, also contains this miRNA module but does not seem to use it to confine early C gene expression to the center of the flower.     

Evidence for gene transfer and expression of factor IX in haemophilia B patients treated with an AAV vector

Pre-clinical studies in mice and haemophilic dogs have shown that introduction of an adeno-associated viral (AAV) vector encoding blood coagulation factor IX (FIX) into skeletal muscle results in sustained expression of F.IX at levels sufficient to correct the haemophilic phenotype. On the basis of these data and additional pre-clinical studies demonstrating an absence of vector-related toxicity, we initiated a clinical study of intramuscular injection of an AAV vector expressing human F.IX in adults with severe haemophilia B. The study has a dose-escalation design, and all patients have now been enrolled in the initial dose cohort (2 x 10(11) vg/kg). Assessment in the first three patients of safety and gene transfer and expression show no evidence of germline transmission of vector sequences or formation of inhibitory antibodies against F.IX. We found that the vector sequences are present in muscle by PCR and Southern-blot analyses of muscle biopsies and we demonstrated expression of F.IX by immunohistochemistry. We observed modest changes in clinical endpoints including circulating levels of F.IX and frequency of FIX protein infusion. The evidence of gene expression at low doses of vector suggests that dose calculations based on animal data may have overestimated the amount of vector required to achieve therapeutic levels in humans, and that the approach offers the possibility of converting severe haemophilia B to a milder form of the disease.