Highly Conserved Noncoding Elements Research Articles

Accelerated evolution at chaperone promoters among Antarctic notothenioid fishes

BackgroundAntarctic fishes of the Notothenioidei suborder constitutively upregulate multiple inducible chaperones, a highly derived adaptation that preserves proteostasis in extreme cold, and represent a system for studying the evolution of gene frontloading. We screened for Hsf1-binding sites, as Hsf1 is a master transcription factor of the heat shock response, and highly-conserved non-coding elements within proximal promoters of chaperone genes across 10 Antarctic notothens, 2 subpolar notothens, and 17 perciform fishes. We employed phylogenetic models of molecular evolution to determine whether (i) changes in motifs associated with Hsf1-binding and/or (ii) relaxed purifying selection or exaptation at ancestral cis-regulatory elements coincided with the evolution of chaperone frontloading in Antarctic notothens.ResultsAntarctic notothens exhibited significantly fewer Hsf1-binding sites per bp at chaperone promoters than subpolar notothens and Serranoidei, the most closely-related suborder to Notothenioidei included in this study. 90% of chaperone promoters exhibited accelerated substitution rates among Antarctic notothens relative to other perciformes. The proportion of bases undergoing accelerated evolution (i) was significantly greater in Antarctic notothens than in subpolar notothens and Perciformes in 70% of chaperone genes and (ii) increased among bases that were more conserved among perciformes. Lastly, we detected evidence of relaxed purifying selection and exaptation acting on ancestrally conserved cis-regulatory elements in the Antarctic notothen lineage and its major branches.ConclusionA large degree of turnover has occurred in Notothenioidei at chaperone promoter regions that are conserved among perciform fishes following adaptation to the cooling of the Southern Ocean. Additionally, derived reductions in Hsf1-binding site frequency suggest cis-regulatory modifications to the classical heat shock response. Of note, turnover events within chaperone promoters were less frequent in the ancestral node of Antarctic notothens relative to younger Antarctic lineages. This suggests that cis-regulatory divergence at chaperone promoters may be greater between Antarctic notothen lineages than between subpolar and Antarctic clades. These findings demonstrate that strong selective forces have acted upon cis-regulatory elements of chaperone genes among Antarctic notothens.

Chromatin Heterogeneity and Distribution of Regulatory Elements in the Late-Replicating Intercalary Heterochromatin Domains of Drosophila melanogaster Chromosomes.

Late-replicating domains (intercalary heterochromatin) in the Drosophila genome display a number of features suggesting their organization is quite unique. Typically, they are quite large and encompass clusters of functionally unrelated tissue-specific genes. They correspond to the topologically associating domains and conserved microsynteny blocks. Our study aims at exploring further details of molecular organization of intercalary heterochromatin and has uncovered surprising heterogeneity of chromatin composition in these regions. Using the 4HMM model developed in our group earlier, intercalary heterochromatin regions were found to host chromatin fragments with a particular epigenetic profile. Aquamarine chromatin fragments (spanning 0.67% of late-replicating regions) are characterized as a class of sequences that appear heterogeneous in terms of their decompactization. These fragments are enriched with enhancer sequences and binding sites for insulator proteins. They likely mark the chromatin state that is related to the binding of cis-regulatory proteins. Malachite chromatin fragments (11% of late-replicating regions) appear to function as universal transitional regions between two contrasting chromatin states. Namely, they invariably delimit intercalary heterochromatin regions from the adjacent active chromatin of interbands. Malachite fragments also flank aquamarine fragments embedded in the repressed chromatin of late-replicating regions. Significant enrichment of insulator proteins CP190, SU(HW), and MOD2.2 was observed in malachite chromatin. Neither aquamarine nor malachite chromatin types appear to correlate with the positions of highly conserved non-coding elements (HCNE) that are typically replete in intercalary heterochromatin. Malachite chromatin found on the flanks of intercalary heterochromatin regions tends to replicate earlier than the malachite chromatin embedded in intercalary heterochromatin. In other words, there exists a gradient of replication progressing from the flanks of intercalary heterochromatin regions center-wise. The peculiar organization and features of replication in large late-replicating regions are discussed as possible factors shaping the evolutionary stability of intercalary heterochromatin.

Dynamic Epigenetic Control of Highly Conserved Noncoding Elements

BackgroundMany noncoding genomic loci have remained constant over long evolutionary periods, suggesting that they are exposed to strong selective pressures. The molecular functions of these elements have been partially elucidated, but the fundamental reason for their extreme conservation is still unknown.ResultsTo gain new insights into the extreme selection of highly conserved noncoding elements (HCNEs), we used a systematic analysis of multi-omic data to study the epigenetic regulation of such elements during the development of Drosophila melanogaster. At the sequence level, HCNEs are GC-rich and have a characteristic oligomeric composition. They have higher levels of stable nucleosome occupancy than their flanking regions, and lower levels of mononucleosomes and H3.3, suggesting that these regions reside in compact chromatin. Furthermore, these regions showed remarkable modulations in histone modification and the expression levels of adjacent genes during development. Although HCNEs are primarily initiated late in replication, about 10% were related to early replication origins. Finally, HCNEs showed strong enrichment within lamina-associated domains.ConclusionHCNEs have distinct and protective sequence properties, undergo dynamic epigenetic regulation, and appear to be associated with the structural components of the chromatin, replication origins, and nuclear matrix. These observations indicate that such elements are likely to have essential cellular functions, and offer insights into their epigenetic properties.

Characterization of Genomic Regulatory Domains Conserved across the Genus Drosophila

In both vertebrates and insects, the conservation of local gene order among distantly related species (microsynteny) is higher than expected in the presence of highly conserved noncoding elements (HCNEs). Dense clusters of HCNEs, or HCNE peaks, have been proposed to mediate the regulation of sometimes distantly located genes, which are central for the developmental program of the organism. Thus, the regions encompassing HCNE peaks and their targets in different species would form genomic regulatory domains (GRDs), which should presumably enjoy an enhanced stability over evolutionary time. By leveraging genome rearrangement information from nine Drosophila species and using gene functional and phenotypic information, we performed a comprehensive characterization of the organization of microsynteny blocks harboring HCNE peaks and provide a functional portrait of the putative HCNE targets that reside therein. We found that Drosophila HCNE peaks tend to colocalize more often than expected and to be evenly distributed across chromosomal elements. Putative HCNE peak targets are characterized by a tight association with particular promoter motifs, higher incidence of severe mutant phenotypes, and evidence of a more precise regulation of gene expression during important developmental transitions. As for their physical organization, ∼65% of these putative targets are separated by a median of two genes from their nearest HCNE peaks. These observations represent the first functional portrait of this euchromatic fraction of the Drosophila genome with distinctive evolutionary dynamics, which will facilitate future experimental studies on the interactions between HCNE peaks and their targets in a genetically tractable system such as Drosophila melanogaster.

Genomic features and computational identification of human microRNAs under long-range developmental regulation

BackgroundRecent functional studies have demonstrated that many microRNAs (miRNAs) are expressed by RNA polymerase II in a specific spatiotemporal manner during the development of organisms and play a key role in cell-lineage decisions and morphogenesis. They are therefore functionally related to a number of key protein coding developmental genes, that form genomic regulatory blocks (GRBs) with arrays of highly conserved non-coding elements (HCNEs) functioning as long-range enhancers that collaboratively regulate the expression of their target genes. Given this functional similarity as well as recent zebrafish transgenesis assays showing that the miR-9 family is indeed regulated by HCNEs with enhancer activity, we hypothesized that this type of miRNA regulation is prevalent. In this paper, we therefore systematically investigate the regulatory landscape around conserved self-transcribed miRNAs (ST miRNAs), with their own known or computationally inferred promoters, by analyzing the hallmarks of GRB target genes. These include not only the density of HCNEs in their vicinity but also the presence of large CpG islands (CGIs) and distinct patterns of histone modification marks associated with developmental genes.ResultsOur results show that a subset of the conserved ST miRNAs we studied shares properties similar to those of protein-coding GRB target genes: they are located in regions of significantly higher HCNE/enhancer binding density and are more likely to be associated with CGIs. Furthermore, their putative promoters have both activating as well as silencing histone modification marks during development and differentiation. Based on these results we used both an elevated HCNE density in the genomic vicinity as well as the presence of a bivalent promoter to identify 29 putative GRB target miRNAs/miRNA clusters, over two-thirds of which are known to play a role during development and differentiation. Furthermore these predictions include miRNAs of the miR-9 family, which are the only experimentally verified GRB target miRNAs.ConclusionsA subset of the conserved miRNA loci we investigated exhibits typical characteristics of GRB target genes, which may partially explain their complex expression profiles during development.

Dissection of the MYCN locus in Feingold syndrome and isolated oesophageal atresia

Feingold syndrome (FS) is a syndromic microcephaly entity for which MYCN is the major disease-causing gene. We studied the expression pattern of MYCN at different stages of human embryonic development and collected a series of 17 FS and 12 isolated oesophageal atresia (IOA) cases. An MYCN gene deletion/mutation was identified in 47% of FS cases exclusively. We hypothesized that mutations or deletions of highly conserved non-coding elements (HCNEs) at the MYCN locus could lead to its misregulation and thereby to FS and/or IOA. We subsequently sequenced five HCNEs at the MYCN locus and designed a high-density tiling path comparative genomic hybridization array of 3.3 Mb at the MYCN locus. We found no mutations or deletions in this region, supporting the hypothesis of genetic heterogeneity in FS.

Transcriptional features of genomic regulatory blocks

CAGE tag mapping of transcription start sites across different human tissues shows that genomic regulatory blocks have unique features that are the likely cause of their ability to respond to regulatory inputs from very long distances.

Hoxa9 and Meis1 Bind Highly Conserved Elements Near Targets Regulated in Leukemia Cells.

The homeobox transcription factors Hoxa9 and Meis1 are critical mediators of transformation by MLL fusion proteins and have also been associated with poor prognosis in other AML subtypes. Despite their pivotal role in leukemia, their mechanism of cooperation and their direct target genes are largely unknown. We have established an experimental model to address these questions by transforming murine hematopoietic progenitors with epitope-tagged forms of Hoxa9 and Meis1. Chromatin immunoprecipitation and massively parallel Illumina sequencing (ChIP-seq) was used to identify binding sites for Hoxa9 and Meis1 across the genome in leukemic cells. We generated 1.7 Gb of sequence data for immunoprecipitated DNA and an untagged control, using two biological replicates, and filtered out (1) regions with nonspecific binding, (2) enriched regions that overlapped certain types of repetitive elements, and (3) nonreproducible enriched regions. We validated a subset of enriched regions by ChIP-chip and ChIP-PCR. De novo motif discovery verified that the resulting set of regions contained a motif that corresponded (E value = 6E-16) to a published Hoxa9 binding motif (Shen et al.; Mol. Cell. Biol. 1999, 19:3051–61). Overlapping or immediately adjacent binding of Hoxa9 and Meis1 was seen for 52% of Hoxa9 and 33% of Meis1 enriched regions, suggesting that these two transcription factors cooperate directly in oncogenesis rather than acting through parallel pathways. Only 7% of enriched regions fell within 2 Kb of known transcriptional start sites, while 46% of regions were within gene boundaries; the remaining peaks were intergenic. Remarkably, given the low overall level of intergenic sequence conservation, over 95% of enriched regions overlapped evolutionarily highly conserved regions (Phastcons for 17 mammalian species >50%; P value<0.0001 [compared to untagged control peaks]). Our experiments suggest that targets of Hoxa9 and Meis1 in both hematopoiesis and leukemic transformation are regulated in part by highly conserved non-coding elements (HCNE), which are recently recognized genomic elements whose functions remain elusive and which are typically absent from widely used promoter microarrays. RNA expression analysis by our group and others confirms that Hoxa9 and Meis1 influence a large number of downstream targets. We used an inducible Hoxa9 expression system to identify genes that were differentially expressed at 72, 96, and 120 h after Hoxa9 withdrawal; 753, 1116, and 2975 transcripts showed >2-fold change, respectively. Phenotypic changes are evident by 96 h; therefore, the 72 h timepoint is most informative of direct Hoxa9 targets. We compared these data to genes containing ChIP-seq peaks and genes with transcription start sites closest to intergenic peaks. ~10% of genes altered at 72 h were associated with a ChIP-seq peak. Another group of investigators identified targets of Meis1 through comparison of Hoxa9- and Hoxa9+Meis1-immortalized cells (Wang et al.; Blood 2005, 106:254–63); ~14% of reported Meis1 targets were associated with a ChIP-seq peak. Hoxa9 and Meis1 binding sites in selected target genes, such the tyrosine kinase Flt3 and the surface marker Cd34, have been confirmed by ChIP-PCR/ChIP-chip. Our current data support a model in which Hoxa9 and Meis1 directly modulate a substantial number of downstream effectors through co-occupancy of regulatory motifs within proximal promoters, intragenic elements, and distally located regions. Studies are underway to determine whether Hoxa9/Meis1 binding influences gene expression and epigenetic modification at promoters and particularly at HCNEs.

Systematic human/zebrafish comparative identification of cis-regulatory activity around vertebrate developmental transcription factor genes

Pan-vertebrate developmental cis-regulatory elements are discernible as highly conserved noncoding elements (HCNEs) and are often dispersed over large areas around the pleiotropic genes whose expression they control. On the loci of two developmental transcription factor genes, SOX3 and PAX 6, we demonstrate that HCNEs conserved between human and zebrafish can be systematically and reliably tested for their regulatory function in multiple stable transgenes in zebrafish, and their genomic reach estimated with confidence using synteny conservation and HCNE density along these loci. HCNEs of both human and zebrafish function as specific developmental enhancers in zebrafish. We show that human HCNEs result in expression patterns in zebrafish equivalent to those in mouse, establishing zebrafish as a suitable model for large-scale testing of human developmental enhancers. Orthologous human and zebrafish enhancers underwent functional evolution within their sequence and often directed related but non-identical expression patterns. Despite an evolutionary distance of 450 million years, one pax6 HCNE drove expression in identical areas when comparing zebrafish vs. human HCNEs. HCNEs from the same area often drive overlapping patterns, suggesting that multiple regulatory inputs are required to achieve robust and precise complex expression patterns exhibited by developmental genes.

Highly Conserved Non-Coding Sequences and the 18q Critical Region for Short Stature: A Common Mechanism of Disease?

BackgroundIsolated growth hormone deficiency (IGHD) and multiple pituitary hormone deficiency (MPHD) are heterogeneous disorders with several different etiologies and they are responsible for most cases of short stature. Mutations in different genes have been identified but still many patients did not present mutations in any of the known genes. Chromosomal rearrangements may also be involved in short stature and, among others, deletions of 18q23 defined a critical region for the disorder. No gene was yet identified.Methodology/Principal FindingsWe now report a balanced translocation X;18 in a patient presenting a breakpoint in 18q23 that was surprisingly mapped about 500 Kb distal from the short stature critical region. It separated from the flanking SALL3 gene a region enriched in highly conserved non-coding elements (HCNE) that appeared to be regulatory sequences, active as enhancers or silencers during embryonic development.ConclusionWe propose that, during pituitary development, the 18q rearrangement may alter expression of 18q genes or of X chromosome genes that are translocated next to the HCNEs. Alteration of expression of developmentally regulated genes by translocation of HCNEs may represent a common mechanism for disorders associated to isolated chromosomal rearrangements.

Ancora: a web resource for exploring highly conserved noncoding elements and their association with developmental regulatory genes

Ancora is a web resource that provides data and tools for exploring genomic organization of highly conserved noncoding elements for multiple genomes.

Genomic regulatory blocks underlie extensive microsynteny conservation in insects

Insect genomes contain larger blocks of conserved gene order (microsynteny) than would be expected under a random breakage model of chromosome evolution. We present evidence that microsynteny has been retained to keep large arrays of highly conserved noncoding elements (HCNEs) intact. These arrays span key developmental regulatory genes, forming genomic regulatory blocks (GRBs). We recently described GRBs in vertebrates, where most HCNEs function as enhancers and HCNE arrays specify complex expression programs of their target genes. Here we present a comparison of five Drosophila genomes showing that HCNE density peaks centrally in large synteny blocks containing multiple genes. Besides developmental regulators that are likely targets of HCNE enhancers, HCNE arrays often span unrelated neighboring genes. We describe differences in core promoters between the target genes and the unrelated genes that offer an explanation for the differences in their responsiveness to enhancers. We show examples of a striking correspondence between boundaries of synteny blocks, HCNE arrays, and Polycomb binding regions, confirming that the synteny blocks correspond to regulatory domains. Although few noncoding elements are highly conserved between Drosophila and the malaria mosquito Anopheles gambiae, we find that A. gambiae regions orthologous to Drosophila GRBs contain an equivalent distribution of noncoding elements highly conserved in the yellow fever mosquito Aëdes aegypti and coincide with regions of ancient microsynteny between Drosophila and mosquitoes. The structural and functional equivalence between insect and vertebrate GRBs marks them as an ancient feature of metazoan genomes and as a key to future studies of development and gene regulation.

Retroviral enhancer detection insertions in zebrafish combined with comparative genomics reveal genomic regulatory blocks - a fundamental feature of vertebrate genomes

A large-scale enhancer detection screen was performed in the zebrafish using a retroviral vector carrying a basal promoter and a fluorescent protein reporter cassette. Analysis of insertional hotspots uncovered areas around developmental regulatory genes in which an insertion results in the same global expression pattern, irrespective of exact position. These areas coincide with vertebrate chromosomal segments containing identical gene order; a phenomenon known as conserved synteny and thought to be a vestige of evolution. Genomic comparative studies have found large numbers of highly conserved noncoding elements (HCNEs) spanning these and other loci. HCNEs are thought to act as transcriptional enhancers based on the finding that many of those that have been tested direct tissue specific expression in transient or transgenic assays. Although gene order in hox and other gene clusters has long been known to be conserved because of shared regulatory sequences or overlapping transcriptional units, the chromosomal areas found through insertional hotspots contain only one or a few developmental regulatory genes as well as phylogenetically unrelated genes. We have termed these regions genomic regulatory blocks (GRBs), and show that they underlie the phenomenon of conserved synteny through all sequenced vertebrate genomes. After teleost whole genome duplication, a subset of GRBs were retained in two copies, underwent degenerative changes compared with tetrapod loci that exist as single copy, and that therefore can be viewed as representing the ancestral form. We discuss these findings in light of evolution of vertebrate chromosomal architecture and the identification of human disease mutations.

A Bivalent Chromatin Structure Marks Key Developmental Genes in Embryonic Stem Cells

The most highly conserved noncoding elements (HCNEs) in mammalian genomes cluster within regions enriched for genes encoding developmentally important transcription factors (TFs). This suggests that HCNE-rich regions may contain key regulatory controls involved in development. We explored this by examining histone methylation in mouse embryonic stem (ES) cells across 56 large HCNE-rich loci. We identified a specific modification pattern, termed "bivalent domains," consisting of large regions of H3 lysine 27 methylation harboring smaller regions of H3 lysine 4 methylation. Bivalent domains tend to coincide with TF genes expressed at low levels. We propose that bivalent domains silence developmental genes in ES cells while keeping them poised for activation. We also found striking correspondences between genome sequence and histone methylation in ES cells, which become notably weaker in differentiated cells. These results highlight the importance of DNA sequence in defining the initial epigenetic landscape and suggest a novel chromatin-based mechanism for maintaining pluripotency.