Conserved Noncoding Sequences Research Articles

New cis-regulatory elements in the Rht-D1b locus region of wheat

Fifteen gene-containing BACs with accumulated length of 1.82-Mb from the Rht-D1b locus region were sequenced and compared in detail with the orthologous regions of rice, sorghum, and maize. Our results show that Rht-D1b represents a conserved genomic region as implied by high gene sequence identity, good maintenance of gene colinearity, and the presence of multiple conserved noncoding sequences (CNSs) that are shared by other grass species. Eight cis-regulatory elements in these CNSs around grass DELLA genes were detected.

Evolutionary growth process of highly conserved sequences in vertebrate genomes

Genome sequence comparison between evolutionarily distant species revealed ultraconserved elements (UCEs) among mammals under strong purifying selection. Most of them were also conserved among vertebrates. Because they tend to be located in the flanking regions of developmental genes, they would have fundamental roles in creating vertebrate body plans. However, the evolutionary origin and selection mechanism of these UCEs remain unclear. Here we report that UCEs arose in primitive vertebrates, and gradually grew in vertebrate evolution. We searched for UCEs in two teleost fishes, Tetraodon nigroviridis and Oryzias latipes, and found 554 UCEs with 100% identity over 100bps. Comparison of teleost and mammalian UCEs revealed 43 pairs of common, jawed-vertebrate UCEs (jUCE) with high sequence identities, ranging from 83.1% to 99.2%. Ten of them retain lower similarities to the Petromyzon marinus genome, and the substitution rates of four non-exonic jUCEs were reduced after the teleost-mammal divergence, suggesting that robust conservation had been acquired in the jawed vertebrate lineage. Our results indicate that prototypical UCEs originated before the divergence of jawed and jawless vertebrates and have been frozen as perfect conserved sequences in the jawed vertebrate lineage. In addition, our comparative sequence analyses of UCEs and neighboring regions resulted in a discovery of lineage-specific conserved sequences. They were added progressively to prototypical UCEs, suggesting step-wise acquisition of novel regulatory roles. Our results indicate that conserved non-coding elements (CNEs) consist of blocks with distinct evolutionary history, each having been frozen since different evolutionary era along the vertebrate lineage.

The Medaka zic1/zic4 Mutant Provides Molecular Insights into Teleost Caudal Fin Evolution

Teleosts have an asymmetrical caudal fin skeleton formed by the upward bending of the caudal-most portion of the body axis, the ural region. This homocercal type of caudal fin ensures powerful and complex locomotion and is regarded as one of the most important innovations for teleosts during adaptive radiation in an aquatic environment. However, the mechanisms that create asymmetric caudal fin remain largely unknown. The spontaneous medaka (teleost fish) mutant, Double anal fin (Da), exhibits a unique symmetrical caudal skeleton that resembles the diphycercal type seen in Polypterus and Coelacanth. We performed a detailed analysis of the Da mutant to obtain molecular insight into caudal fin morphogenesis. We first demonstrate that a large transposon, inserted into the enhancer region of the zic1 and zic4 genes (zic1/zic4) in Da, is associated with the mesoderm-specific loss of their transcription. We then show that zic1/zic4 are strongly expressed in the dorsal part of the ural mesenchyme and thereby induce asymmetric caudal fin development in wild-type embryos, whereas their expression is lost in Da. Comparative analysis further indicates that the dorsal mesoderm expression of zic1/zic4 is conserved in teleosts, highlighting the crucial role of zic1/zic4 in caudal fin morphogenesis.

Cis-regulation and conserved non-coding elements in amphioxus

The origin of the vertebrates was a major event in the evolution of morphological diversity and the genetic mechanisms responsible for this diversity, once purely theoretical, can now be approached experimentally in the genome era. With a prototypical chordate genome, vertebrate-like development and simple morphology, amphioxus provides the appropriate model for investigating the origin of the vertebrates. Comparative genomics is revealing that both conservation and divergence of genes and cis-regulatory elements involved in developmental regulatory networks are required to shape different animal body plans. This article reviews the cis-regulatory studies performed in amphioxus, the discovery of conserved non-coding elements (CNEs) across the metazoans and the examination of amphioxus CNEs. Emerging ideas on the evolution of CNEs after large-scale genome duplication events and the state of cephalochordate genomics are also discussed.

Lens development depends on a pair of highly conserved Sox21 regulatory elements

Highly conserved non-coding elements (CNEs) linked to genes involved in embryonic development have been hypothesised to correspond to cis-regulatory modules due to their ability to induce tissue-specific expression patterns. However, attempts to prove their requirement for normal development or for the correct expression of the genes they are associated with have yielded conflicting results. Here, we show that CNEs at the vertebrate Sox21 locus are crucial for Sox21 expression in the embryonic lens and that loss of Sox21 function interferes with normal lens development. Using different expression assays in zebrafish we find that two CNEs linked to Sox21 in all vertebrates contain lens enhancers and that their removal from a reporter BAC abolishes lens expression. Furthermore inhibition of Sox21 function after the injection of a sox21b morpholino into zebrafish leads to defects in lens development. These findings identify a direct link between sequence conservation and genomic function of regulatory sequences. In addition to this we provide evidence that putative Sox binding sites in one of the CNEs are essential for induction of lens expression as well as enhancer function in the CNS. Our results show that CNEs identified in pufferfish-mammal whole-genome comparisons are crucial developmental enhancers and hence essential components of gene regulatory networks underlying vertebrate embryogenesis.

Diverse Functions of Distal Regulatory Elements at the IFNG Locus

Previous studies have identified multiple conserved noncoding sequences (CNS) at the mouse Ifng locus sufficient for enhancer activity in cell-based assays. These studies do not directly address biology of the human IFNG locus in a genomic setting. IFNG enhancers may be functionally redundant or each may be functionally unique. We test the hypothesis that each IFNG enhancer has a unique necessary function using a bacterial artificial chromosome transgenic model. We find that CNS-30, CNS-4, and CNS+20 are required at distinct stages of Th1 differentiation, whereas CNS-16 has a repressive role in Th1 and Th2 cells. CNS+20 is required for IFN-γ expression by memory Th1 cells and NKT cells. CNS-4 is required for IFN-γ expression by effector Th1 cells. In contrast, CNS-16, CNS-4, and CNS+20 are each partially required for human IFN-γ expression by NK cells. Thus, IFNG CNS enhancers have redundant necessary functions in NK cells but unique necessary functions in Th cells. These results also demonstrate that distinct CNSs are required to transcribe IFNG at each stage of the Th1 differentiation pathway.

The 3′ Enhancer CNS2 Is a Critical Regulator of Interleukin-4-Mediated Humoral Immunity in Follicular Helper T Cells

A main role for interleukin-4 (IL-4) is in humoral immunity, and follicular helper CD4(+) T (Tfh) cells may be an intrinsic IL-4 source. Here we demonstrate that conserved noncoding sequence 2 (CNS2) is an essential enhancer element for IL-4 expression in Tfh cells but not in Th2 cells. Mice with a CNS2 deletion had a reduction in IgG1 and IgE production and in IL-4 expression in Tfh cells. Tracking of CNS2 activity via a GFP reporter mouse demonstrated that CNS2-active cells expressed several markers of Tfh cells: CXCR5, PD-1, and ICOS; the transcriptional master regulator Bcl6; and the cytokines IL-21 and IL-4. These CNS2-active cells were mainly localized in B cell follicles and germinal centers. The GFP(+) Tfh cells were derived from GFP(-) naive T cells after in vivo systemic immunization. These results indicate that CNS2 is an essential enhancer element required for IL-4 expression in Tfh cells controlling humoral immunity.

Transcription of Il17 and Il17f Is Controlled by Conserved Noncoding Sequence 2

T helper 17 (Th17) cells specifically transcribe the Il17 and Il17f genes, which are localized in the same chromosome region, but the underlying mechanism is unclear. Here, we report a cis element that we previously named conserved noncoding sequence 2 (CNS2) physically interacted with both Il17 and Il17f gene promoters and was sufficient for regulating their selective transcription in Th17 cells. Targeted deletion of CNS2 resulted in impaired retinoic acid-related orphan receptor gammat (RORγt)-driven IL-17 expression in vitro. CNS2-deficient T cells also produced substantially decreased amounts of IL-17F. These cytokine defects were associated with defective chromatin remodeling in the Il17-Il17f gene locus, possibly because of effects on CNS2-mediated recruitment of histone-modifying enzymes p300 and JmjC domain-containing protein 3 (JMJD3). CNS2-deficient animals were also shown to be resistant to experimental autoimmune encephalomyelitis (EAE). Our results thus suggest that CNS2 is sufficient and necessary for Il17 and optimal Il17f gene transcription in Th17 cells.

Ancient Pbx-Hox signatures define hundreds of vertebrate developmental enhancers

BackgroundGene regulation through cis-regulatory elements plays a crucial role in development and disease. A major aim of the post-genomic era is to be able to read the function of cis-regulatory elements through scrutiny of their DNA sequence. Whilst comparative genomics approaches have identified thousands of putative regulatory elements, our knowledge of their mechanism of action is poor and very little progress has been made in systematically de-coding them.ResultsHere, we identify ancient functional signatures within vertebrate conserved non-coding elements (CNEs) through a combination of phylogenetic footprinting and functional assay, using genomic sequence from the sea lamprey as a reference. We uncover a striking enrichment within vertebrate CNEs for conserved binding-site motifs of the Pbx-Hox hetero-dimer. We further show that these predict reporter gene expression in a segment specific manner in the hindbrain and pharyngeal arches during zebrafish development.ConclusionsThese findings evoke an evolutionary scenario in which many CNEs evolved early in the vertebrate lineage to co-ordinate Hox-dependent gene-regulatory interactions that pattern the vertebrate head. In a broader context, our evolutionary analyses reveal that CNEs are composed of tightly linked transcription-factor binding-sites (TFBSs), which can be systematically identified through phylogenetic footprinting approaches. By placing a large number of ancient vertebrate CNEs into a developmental context, our findings promise to have a significant impact on efforts toward de-coding gene-regulatory elements that underlie vertebrate development, and will facilitate building general models of regulatory element evolution.

Functional Analysis of Conserved Non-Coding Regions Around the Short Stature hox Gene (shox) in Whole Zebrafish Embryos

BackgroundMutations in the SHOX gene are responsible for Leri-Weill Dyschondrosteosis, a disorder characterised by mesomelic limb shortening. Recent investigations into regulatory elements surrounding SHOX have shown that deletions of conserved non-coding elements (CNEs) downstream of the SHOX gene produce a phenotype indistinguishable from Leri-Weill Dyschondrosteosis. As this gene is not found in rodents, we used zebrafish as a model to characterise the expression pattern of the shox gene across the whole embryo and characterise the enhancer domains of different CNEs associated with this gene.Methodology/Principal FindingsExpression of the shox gene in zebrafish was identified using in situ hybridization, with embryos showing expression in the blood, putative heart, hatching gland, brain pharyngeal arch, olfactory epithelium, and fin bud apical ectodermal ridge. By identifying sequences showing 65% identity over at least 40 nucleotides between Fugu, human, dog and opossum we uncovered 35 CNEs around the shox gene. These CNEs were compared with CNEs previously discovered by Sabherwal et al., resulting in the identification of smaller more deeply conserved sub-sequence. Sabherwal et al.'s CNEs were assayed for regulatory function in whole zebrafish embryos resulting in the identification of additional tissues under the regulatory control of these CNEs.Conclusion/SignificanceOur results using whole zebrafish embryos have provided a more comprehensive picture of the expression pattern of the shox gene, and a better understanding of its regulation via deeply conserved noncoding elements. In particular, we identify additional tissues under the regulatory control of previously identified SHOX CNEs. We also demonstrate the importance of these CNEs in evolution by identifying duplicated shox CNEs and more deeply conserved sub-sequences within already identified CNEs.

Epigenetic regulation in murine offspring as a novel mechanism for transmaternal asthma protection induced by microbes

Bronchial asthma is a chronic inflammatory disease resulting from complex gene-environment interactions. Natural microbial exposure has been identified as an important environmental condition that provides asthma protection in a prenatal window of opportunity. Epigenetic regulation is an important mechanism by which environmental factors might interact with genes involved in allergy and asthma development. This study was designed to test whether epigenetic mechanisms might contribute to asthma protection conferred by early microbial exposure. Pregnant maternal mice were exposed to the farm-derived gram-negative bacterium Acinetobacter lwoffii F78. Epigenetic modifications in the offspring were analyzed in T(H)1- and T(H)2-relevant genes of CD4(+) T cells. Prenatal administration of A lwoffii F78 prevented the development of an asthmatic phenotype in the progeny, and this effect was IFN-γ dependent. Furthermore, the IFNG promoter of CD4(+) T cells in the offspring revealed a significant protection against loss of histone 4 (H4) acetylation, which was closely associated with IFN-γ expression. Pharmacologic inhibition of H4 acetylation in the offspring abolished the asthma-protective phenotype. Regarding T(H)2-relevant genes only at the IL4 promoter, a decrease could be detected for H4 acetylation but not at the IL5 promoter or the intergenic T(H)2 regulatory region conserved noncoding sequence 1 (CNS1). These data support the hygiene concept and indicate that microbes operate by means of epigenetic mechanisms. This provides a new mechanism in the understanding of gene-environment interactions in the context of allergy protection.

The Prx1 limb enhancers: Targeted gene expression in developing zebrafish pectoral fins

Limbs represent an excellent model to study the induction, growth, and patterning of several organs. A breakthrough to study gene function in various tissues has been the characterization of regulatory elements that allow tissue-specific interference of gene function. The mouse Prx1 promoter has been used to generate limb-specific mutants and overexpress genes in tetrapod limbs. Although zebrafish possess advantages that favor their use to study limb morphogenesis, there is no driver described suitable for specifically interfering with gene function in developing fins. We report the generation of zebrafish lines that express enhanced green fluorescent protein (EGFP) driven by the mouse Prx1 enhancer in developing pectoral fins. We also describe the expression pattern of the zebrafish prrx1 genes and identify three conserved non-coding elements (CNEs) that we use to generate fin-specific EGFP reporter lines. Finally, we show that the mouse and zebrafish regulatory elements may be used to modify gene function in pectoral fins.

Mouse transgenesis identifies conserved functional enhancers and cis-regulatory motif in the vertebrate LIM homeobox gene Lhx2 locus.

The vertebrate Lhx2 is a member of the LIM homeobox family of transcription factors. It is essential for the normal development of the forebrain, eye, olfactory system and liver as well for the differentiation of lymphoid cells. However, despite the highly restricted spatio-temporal expression pattern of Lhx2, nothing is known about its transcriptional regulation. In mammals and chicken, Crb2, Dennd1a and Lhx2 constitute a conserved linkage block, while the intervening Dennd1a is lost in the fugu Lhx2 locus. To identify functional enhancers of Lhx2, we predicted conserved noncoding elements (CNEs) in the human, mouse and fugu Crb2-Lhx2 loci and assayed their function in transgenic mouse at E11.5. Four of the eight CNE constructs tested functioned as tissue-specific enhancers in specific regions of the central nervous system and the dorsal root ganglia (DRG), recapitulating partial and overlapping expression patterns of Lhx2 and Crb2 genes. There was considerable overlap in the expression domains of the CNEs, which suggests that the CNEs are either redundant enhancers or regulating different genes in the locus. Using a large set of CNEs (810 CNEs) associated with transcription factor-encoding genes that express predominantly in the central nervous system, we predicted four over-represented 8-mer motifs that are likely to be associated with expression in the central nervous system. Mutation of one of them in a CNE that drove reporter expression in the neural tube and DRG abolished expression in both domains indicating that this motif is essential for expression in these domains. The failure of the four functional enhancers to recapitulate the complete expression pattern of Lhx2 at E11.5 indicates that there must be other Lhx2 enhancers that are either located outside the region investigated or divergent in mammals and fishes. Other approaches such as sequence comparison between multiple mammals are required to identify and characterize such enhancers.

Informatic and Functional Approaches to Identifying a Regulatory Region for the Cardiac Sodium Channel

Although multiple lines of evidence suggest that variable expression of the cardiac sodium channel gene SCN5A plays a role in susceptibility to arrhythmia, little is known about its transcriptional regulation. We used in silico and in vitro experiments to identify possible noncoding sequences important for transcriptional regulation of SCN5A. The results were extended to mice in which a putative regulatory region was deleted. We identified 92 noncoding regions highly conserved (>70%) between human and mouse SCN5A orthologs. Three conserved noncoding sequences (CNS) showed significant (>5-fold) activity in luciferase assays. Further in vitro studies indicated one, CNS28 in intron 1, as a potential regulatory region. Using recombinase-mediated cassette exchange (RMCE), we generated mice in which a 435-base pair region encompassing CNS28 was removed. Animals homozygous for the deletion showed significant increases in SCN5A transcripts, Na(V)1.5 protein abundance, and sodium current measured in isolated ventricular myocytes. ECGs revealed a significantly shorter QRS (10.7±0.2 ms in controls versus 9.7±0.2 ms in knockouts), indicating more rapid ventricular conduction. In vitro analysis of CNS28 identified a short 3' segment within this region required for regulatory activity and including an E-box motif. Deletion of this segment reduced reporter activity to 3.6%±0.3% of baseline in CHO cells and 16%±3% in myocytes (both P<0.05), and mutation of individual sites in the E-box restored activity to 62%±4% and 57%±2% of baseline in CHO cells and myocytes, respectively (both P<0.05). These findings establish that regulation of cardiac sodium channel expression modulates channel function in vivo, and identify a noncoding region underlying this regulation.

Mechanisms of tolerance

Mechanisms of tolerance

Positive and Negative Selection in Murine Ultraconserved Noncoding Elements

There are many more selectively constrained noncoding than coding nucleotides in the mammalian genome, but most mammalian noncoding DNA is subject to weak selection, on average. One of the most striking discoveries to have emerged from comparisons among mammalian genomes is the hundreds of noncoding elements of more than 200 bp in length that show absolute conservation among mammalian orders. These elements represent the tip of the iceberg of a much larger class of conserved noncoding elements (CNEs). Much evidence suggests that CNEs are selectively constrained and not mutational cold-spots, and there is evidence that some CNEs play a role in the regulation of development. Here, we quantify negative and positive selection acting in murine CNEs by analyzing within-species nucleotide variation and between-species divergence of CNEs that we identified using a phylogenetically independent comparison. The distribution of fitness effects of new mutations in CNEs, inferred from within-species polymorphism, suggests that CNEs receive a higher number of strongly selected deleterious mutations and many fewer nearly neutral mutations than amino acid sites of protein-coding genes or regulatory elements close to genes. However, we also show that CNEs experience a far higher proportion of adaptive substitutions than any known category of genomic sites in murids. The absolute rate of adaptation of CNEs is similar to that of amino acid sites of proteins. This result suggests that there is widespread adaptation in mammalian conserved noncoding DNA elements, some of which have been implicated in the regulation of crucially important processes, including development.

Evolutionary divergence and limits of conserved non-coding sequence detection in plant genomes.

The discovery of regulatory motifs embedded in upstream regions of plants is a particularly challenging bioinformatics task. Previous studies have shown that motifs in plants are short compared with those found in vertebrates. Furthermore, plant genomes have undergone several diversification mechanisms such as genome duplication events which impact the evolution of regulatory motifs. In this article, a systematic phylogenomic comparison of upstream regions is conducted to further identify features of the plant regulatory genomes, the component of genomes regulating gene expression, to enable future de novo discoveries. The findings highlight differences in upstream region properties between major plant groups and the effects of divergence times and duplication events. First, clear differences in upstream region evolution can be detected between monocots and dicots, thus suggesting that a separation of these groups should be made when searching for novel regulatory motifs, particularly since universal motifs such as the TATA box are rare. Second, investigating the decay rate of significantly aligned regions suggests that a divergence time of ∼100 mya sets a limit for reliable conserved non-coding sequence (CNS) detection. Insights presented here will set a framework to help identify embedded motifs of functional relevance by understanding the limits of bioinformatics detection for CNSs.

Comparative analysis of duplicated sox21 genes in zebrafish

Sox21 is thought to function as a counteracting partner of SoxB1 (Sox1, 2, 3) genes and is involved in cell fate determination. In this study, we comparatively analyzed the expression patterns and conserved cis-regulatory elements of the duplicated sox21 genes in zebrafish. In embryogenesis, sox21b is predominantly expressed in the telencephalon, hypothalamus, mesencephalon and lens, and sox21a is solely expressed in the midbrain-hindbrain boundary, olfactory placode and lateral line, while both genes are expressed in the hindbrain, spinal cord and ear. In adult, sox21a is expressed in the brain, skin, ovary and intestine, while sox21b is expressed in the brain and testis. Interestingly, all 16 pan-vertebrate conserved non-coding elements (CNEs) are asymmetrically preserved in the sox21b locus, whereas two fish-specific elements are kept in the sox21a locus, and this is correlated with increased evolutionary rate of the sox21a protein sequence. Transient transgenic reporter analysis revealed that six sox21b CNEs and two sox21a CNEs drove green fluorescent protein (GFP) expression in tissues correlated with the partitioning of expression in two orthologues. These results indicate that sox21a and sox21b have reciprocally lost expression domains of the ancestral gene reflected by degeneration of certain CNEs in their genomic loci and provide clear evidence for evolution of the duplicated sox21 genes by subfunctionalization. In addition, our data suggest that some CNEs-based regulatory pathways have been predominantly preserved in the sox21b locus.

Identification of two novel mutations in Shh long-range regulator associated with familial pre-axial polydactyly

Pre-axial polydactyly type II (PPDII, MIM #174500), Werner mesomelic syndrome (MIM %188770) and Haas polysyndactyly (MIM #186200) are a group of closely related conditions caused by mutations in a long-range Sonic hedgehog (SHH, MIM *600725) regulator called ZRS. To date, 19 point mutations, 10 duplications and 1 triplication of the ZRS associated with those pre-axial polydactylies have been reported in humans, mice, cats and chickens. Some of these have been shown to cause ectopic expression of Shh in the limb bud in mice, leading to a polydactylous phenotype, but the precise mechanism by which ZRS mutations generate this phenotype remains unknown. We present two PPDII families with fully penetrant point mutations in ultra-conserved predicted binding sites for transcription factors SOX9 and PAX3, two possible candidates for regulating SHH expression. Screening for point mutations or copy-number variation of the ZRS, high-resolution array-CGH, and screening of other conserved non-coding sequences (CNS) surrounding SHH in a third family are negative. This is the sixth PPDII pedigree with possible linkage to 7q36 that presents with no detectable ZRS mutation. We hypothesize that another nearby regulatory sequence, or an undetected position effect between ZRS and SHH, could be responsible for negative familial cases linked to 7q36.

Minor change, major difference: divergent functions of highly conserved cis-regulatory elements subsequent to whole genome duplication events

Within the vertebrate lineage, a high proportion of duplicate genes have been retained after whole genome duplication (WGD) events. It has been proposed that many of these duplicate genes became indispensable because the ancestral gene function was divided between them. In addition, novel functions may have evolved, owing to changes in cis-regulatory elements. Functional analysis of the PAX2/5/8 gene subfamily appears to support at least the first part of this hypothesis. The collective role of these genes has been widely retained, but sub-functions have been differentially partitioned between the genes in different vertebrates. Conserved non-coding elements (CNEs) represent an interesting and readily identifiable class of putative cis-regulatory elements that have been conserved from fish to mammals, an evolutionary distance of 450 million years. Within the PAX2/5/8 gene subfamily, PAX2 is associated with the highest number of CNEs. An additional WGD experienced in the teleost lineage led to two copies of pax2, each of which retained a large proportion of these CNEs. Using a reporter gene assay in zebrafish embryos, we have exploited this rich collection of regulatory elements in order to determine whether duplicate CNEs have evolved different functions. Remarkably, we find that even highly conserved sequences exhibit more functional differences than similarities. We also discover that short flanking sequences can have a profound impact on CNE function. Therefore, if CNEs are to be used as candidate enhancers for transgenic studies or for multi-species comparative analyses, it is paramount that the CNEs are accurately delineated.