Detection of regulatory variation in mouse genes

Functional polymorphism in genes can be classified as coding variation, altering the amino-acid sequence of the encoded protein, or regulatory variation, affecting the level or pattern of expression of the gene. Coding variation can be recognized directly from DNA sequence, and consequently its frequency and characteristics have been extensively described. By contrast, virtually nothing is known about the extent to which gene regulation varies in populations. Yet it is likely that regulatory variants are important in modulating gene function: alterations in gene regulation have been proposed to influence disease susceptibility and to have been the primary substrate for the evolution of species. Here, we report a systematic study to assess the extent of cis-acting regulatory variation in 69 genes across four inbred mouse strains. We find that at least four of these genes show allelic differences in expression level of 1.5-fold or greater, and that some of these differences are tissue specific. The results show that the impact of regulatory variants can be detected at a significant frequency in a genomic survey and suggest that such variation may have important consequences for organismal phenotype and evolution. The results indicate that larger-scale surveys in both mouse and human could identify a substantial number of genes with common regulatory variation.     

Disruption of an imprinted gene cluster by a targeted chromosomal translocation in mice

Genomic imprinting is an epigenetic process in which the activity of a gene is determined by its parent of origin. Mechanisms governing genomic imprinting are just beginning to be understood. However, the tendency of imprinted genes to exist in chromosomal clusters suggests a sharing of regulatory elements. To better understand imprinted gene clustering, we disrupted a cluster of imprinted genes on mouse distal chromosome 7 using the Cre/loxP recombination system. In mice carrying a site-specific translocation separating Cdkn1c and Kcnq1, imprinting of the genes retained on chromosome 7, including Kcnq1, Kcnq1ot1, Ascl2, H19 and Igf2, is unaffected, demonstrating that these genes are not regulated by elements near or telomeric to Cdkn1c. In contrast, expression and imprinting of the translocated Cdkn1c, Slc22a1l and Tssc3 on chromosome 11 are affected, consistent with the hypothesis that elements regulating both expression and imprinting of these genes lie within or proximal to Kcnq1. These data support the proposal that chromosomal abnormalities, including translocations, within KCNQ1 that are associated with the human disease Beckwith-Wiedemann syndrome (BWS) may disrupt CDKN1C expression. These results underscore the importance of gene clustering for the proper regulation of imprinted genes.     

Interleukin 7 receptor alpha chain (IL7R) shows allelic and functional association with multiple sclerosis

Multiple sclerosis is a demyelinating neurodegenerative disease with a strong genetic component. Previous genetic risk studies have failed to identify consistently linked regions or genes outside of the major histocompatibility complex on chromosome 6p. We describe allelic association of a polymorphism in the gene encoding the interleukin 7 receptor alpha chain (IL7R) as a significant risk factor for multiple sclerosis in four independent family-based or case-control data sets (overall P = 2.9 x 10(-7)). Further, the likely causal SNP, rs6897932, located within the alternatively spliced exon 6 of IL7R, has a functional effect on gene expression. The SNP influences the amount of soluble and membrane-bound isoforms of the protein by putatively disrupting an exonic splicing silencer.     

A 1.5 million-base pair inversion polymorphism in families with Williams-Beuren syndrome.

Williams-Beuren syndrome (WBS) is most often caused by hemizygous deletion of a 1.5-Mb interval encompassing at least 17 genes at 7q11.23 (refs. 1,2). As with many other haploinsufficiency diseases, the mechanism underlying the WBS deletion is thought to be unequal meiotic recombination, probably mediated by the highly homologous DNA that flanks the commonly deleted region. Here, we report the use of interphase fluorescence in situ hybridization (FISH) and pulsed-field gel electrophoresis (PFGE) to identify a genomic polymorphism in families with WBS, consisting of an inversion of the WBS region. We have observed that the inversion is hemizygous in 3 of 11 (27%) atypical affected individuals who show a subset of the WBS phenotypic spectrum but do not carry the typical WBS microdeletion. Two of these individuals also have a parent who carries the inversion. In addition, in 4 of 12 (33%) families with a proband carrying the WBS deletion, we observed the inversion exclusively in the parent transmitting the disease-related chromosome. These results suggest the presence of a newly identified genomic variant within the population that may be associated with the disease. It may result in predisposition to primarily WBS-causing microdeletions, but may also cause translocations and inversions.     

SUMO modification is required for in vivo Hox gene regulation by the Caenorhabditis elegans Polycomb group protein SOP-2

Post-translational modification of proteins by the ubiquitin-like molecule SUMO (sumoylation) regulates their subcellular localization and affects their functional properties in vitro, but the physiological function of sumoylation in multicellular organisms is largely unknown. Here, we show that the C. elegans Polycomb group (PcG) protein SOP-2 interacts with the SUMO-conjugating enzyme UBC-9 through its evolutionarily conserved SAM domain. Sumoylation of SOP-2 is required for its localization to nuclear bodies in vivo and for its physiological repression of Hox genes. Global disruption of sumoylation phenocopies a sop-2 mutation by causing ectopic Hox gene expression and homeotic transformations. Chimeric constructs in which the SOP-2 SAM domain is replaced with that derived from fruit fly or mammalian PcG proteins, but not those in which the SOP-2 SAM domain is replaced with the SAM domains of non-PcG proteins, confer appropriate in vivo nuclear localization and Hox gene repression. These observations indicate that sumoylation of PcG proteins, modulated by their evolutionarily conserved SAM domain, is essential to their physiological repression of Hox genes.     

Epigenetic remodeling in colorectal cancer results in coordinate gene suppression across an entire chromosome band

We report a new mechanism in carcinogenesis involving coordinate long-range epigenetic gene silencing. Epigenetic silencing in cancer has always been envisaged as a local event silencing discrete genes. However, in this study of silencing in colorectal cancer, we found common repression of the entire 4-Mb band of chromosome 2q.14.2, associated with global methylation of histone H3 Lys9. DNA hypermethylation within the repressed genomic neighborhood was localized to three separate enriched CpG island 'suburbs', with the largest hypermethylated suburb spanning 1 Mb. These data change our understanding of epigenetic gene silencing in cancer cells: namely, epigenetic silencing can span large regions of the chromosome, and both DNA-methylated and neighboring unmethylated genes can be coordinately suppressed by global changes in histone modification. We propose that loss of gene expression can occur through long-range epigenetic silencing, with similar implications as loss of heterozygosity in cancer.     

MicroRNA-responsive 'sensor' transgenes uncover Hox-like and other developmentally regulated patterns of vertebrate microRNA expression

MicroRNAs (miRNAs) are a class of short ( approximately 22-nt) noncoding RNA molecules that downregulate expression of their mRNA targets. Since their discovery as regulators of developmental timing in Caenorhabditis elegans, hundreds of miRNAs have been identified in both animals and plants. Here, we report a technique for visualizing detailed miRNA expression patterns in mouse embryos. We elucidate the tissue-specific expression of several miRNAs during embryogenesis, including two encoded by genes embedded in homeobox (Hox) clusters, miR-10a and miR-196a. These two miRNAs are expressed in patterns that are markedly reminiscent of those of Hox genes. Furthermore, miR-196a negatively regulates Hoxb8, indicating that its restricted expression pattern probably reflects a role in the patterning function of the Hox complex.     

A negative element in SMN2 exon 7 inhibits splicing in spinal muscular atrophy

Spinal muscular atrophy (SMA) is a relatively common neurodegenerative disease caused by homozygous loss of the survival motor neuron 1 (SMN1) gene. Humans possess a linked, nearly identical gene, SMN2, which produces a functional SMN protein but at levels insufficient to compensate for loss of SMN1 (refs. 1,2). A C/T transition at position +6 in exon 7 is all that differentiates the two genes, but this is sufficient to prevent efficient exon 7 splicing in SMN2 (refs. 2,3). Here we show that the C/T transition functions not to disrupt an exonic splicing enhancer (ESE) in SMN1 (ref. 4), as previously suggested, but rather to create an exonic splicing silencer (ESS) in SMN2. We show that this ESS functions as a binding site for a known repressor protein, hnRNP A1, which binds to SMN2 but not SMN1 exon 7 RNA. We establish the physiological importance of these results by using small interfering RNAs to reduce hnRNP A protein levels in living cells and show that this results in efficient SMN2 exon 7 splicing. Our findings not only define a new mechanism underlying the inefficient splicing of SMN2 exon 7 but also illustrate more generally the remarkable sensitivity and precision that characterizes control of mRNA splicing.     

Regulation of anterior/posterior patterning of the axial skeleton by growth/differentiation factor 11.

The bones that comprise the axial skeleton have distinct morphological features characteristic of their positions along the anterior/posterior axis. We previously described a novel TGF-beta family member, myostatin (encoded by the gene Mstn, formerly Gdf8), that has an essential role in regulating skeletal muscle mass. We also identified a gene related to Mstn by low-stringency screening. While the work described here was being completed, the cloning of this gene, designated Gdf11 (also called Bmp11), was also reported by other groups. Here we show that Gdf11, a new transforming growth factor beta(TGFbeta) superfamily member, has an important role in establishing this skeletal pattern. During early mouse embryogenesis, Gdf11 is expressed in the primitive streak and tail bud regions, which are sites where new mesodermal cells are generated. Homozygous mutant mice carrying a targeted deletion of Gdf11 exhibit anteriorly directed homeotic transformations throughout the axial skeleton and posterior displacement of the hindlimbs. The effect of the mutation is dose dependent, as Gdf11+/- mice have a milder phenotype than Gdf11-/- mice. Mutant embryos show alterations in patterns of Hox gene expression, suggesting that Gdf11 acts upstream of the Hox genes. Our findings suggest that Gdf11 is a secreted signal that acts globally to specify positional identity along the anterior/posterior axis.     

Rnx deficiency results in congenital central hypoventilation

The genes Tlx1 (Hox11), Enx (Hox11L2, Tlx-2) and Rnx (Hox11L2, Tlx-3) constitute a family of orphan homeobox genes. In situ hybridization has revealed considerable overlap in their expression within the nervous system, but Rnx is singularly expressed in the developing dorsal and ventral region of the medulla oblongata. Tlx1-deficient and Enx-deficient mice display phenotypes in tissues where the mutated gene is singularly expressed, resulting in asplenogenesis and hyperganglionic megacolon, respectively. To determine the developmental role of Rnx, we disrupted the locus in mouse embryonic stem (ES) cells. Rnx deficient mice developed to term, but all died within 24 hours after birth from a central respiratory failure. The electromyographic activity of intercostal muscles coupled with the C4 ventral root activity assessed in a medulla-spinal cord preparation revealed a high respiratory rate with short inspiratory duration and frequent apnea. Furthermore, a coordinate pattern existed between the abnormal activity of inspiratory neurons in the ventrolateral medulla and C4 motorneuron output, indicating a central respiratory defect in Rnx mice. Thus, Rnx is critical for the development of the ventral medullary respiratory centre and its deficiency results in a syndrome resembling congenital central hypoventilation.     

Integrated genomic and epigenomic analyses pinpoint biallelic gene inactivation in tumors

Aberrant methylation of CpG islands and genomic deletion are two predominant mechanisms of gene inactivation in tumorigenesis, but the extent to which they interact is largely unknown. The lack of an integrated approach to study these mechanisms has limited the understanding of tumor genomes and cancer genes. Restriction landmark genomic scanning (RLGS; ref. 1) is useful for global analysis of aberrant methylation of CpG islands, but has not been amenable to alignment with deletion maps because the identity of most RLGS fragments is unknown. Here, we determined the nucleotide sequence and exact chromosomal position of RLGS fragments throughout the genome using the whole chromosome of origin of the fragments and in silico restriction digestion of the human genome sequence. To study the interaction of these gene-inactivation mechanisms in primary brain tumors, we integrated RLGS-based methylation analysis with high-resolution deletion maps from microarray-based comparative genomic hybridization (array CGH; ref. 3). Certain subsets of gene-associated CpG islands were preferentially affected by convergent methylation and deletion, including genes that exhibit tumor-suppressor activity, such as CISH1 (encoding SOCS1; ref. 4), as well as genes such as COE3 that have been missed by traditional non-integrated approaches. Our results show that most aberrant methylation events are focal and independent of deletions, and the rare convergence of these mechanisms can pinpoint biallelic gene inactivation without the use of positional cloning.     

Fine mapping of regulatory loci for mammalian gene expression using radiation hybrids

We mapped regulatory loci for nearly all protein-coding genes in mammals using comparative genomic hybridization and expression array measurements from a panel of mouse-hamster radiation hybrid cell lines. The large number of breaks in the mouse chromosomes and the dense genotyping of the panel allowed extremely sharp mapping of loci. As the regulatory loci result from extra gene dosage, we call them copy number expression quantitative trait loci, or ceQTLs. The -2log10P support interval for the ceQTLs was <150 kb, containing an average of <2-3 genes. We identified 29,769 trans ceQTLs with -log10P > 4, including 13 hotspots each regulating >100 genes in trans. Further, this work identifies 2,761 trans ceQTLs harboring no known genes, and provides evidence for a mode of gene expression autoregulation specific to the X chromosome.