Highly conserved non-coding sequences are associated with vertebrate development

PReMod: a database of genome-wide mammalian cis-regulatory module predictions

We describe PReMod, a new database of genome-wide cis-regulatory module (CRM) predictions for both the human and the mouse genomes. The prediction algorithm, described previously in Blanchette et al. (2006) Genome Res., 16, 656–668, exploits the fact that many known CRMs are made of clusters of phylogenetically conserved and repeated transcription factors (TF) binding sites. Contrary to other existing databases, PReMod is not restricted to modules located proximal to genes, but in fact mostly contains distal predicted CRMs (pCRMs). Through its web interface, PReMod allows users to (i) identify pCRMs around a gene of interest; (ii) identify pCRMs that have binding sites for a given TF (or a set of TFs) or (iii) download the entire dataset for local analyses. Queries can also be refined by filtering for specific chromosomal regions, for specific regions relative to genes or for the presence of CpG islands. The output includes information about the binding sites predicted within the selected pCRMs, and a graphical display of their distribution within the pCRMs. It also provides a visual depiction of the chromosomal context of the selected pCRMs in terms of neighboring pCRMs and genes, all of which are linked to the UCSC Genome Browser and the NCBI. PReMod: .

Experimental validation of predicted mammalian erythroid cis-regulatory modules

Multiple alignments of genome sequences are helpful guides to functional analysis, but predicting cis-regulatory modules (CRMs) accurately from such alignments remains an elusive goal. We predict CRMs for mammalian genes expressed in red blood cells by combining two properties gleaned from aligned, noncoding genome sequences: a positive regulatory potential (RP) score, which detects similarity to patterns in alignments distinctive for regulatory regions, and conservation of a binding site motif for the essential erythroid transcription factor GATA-1. Within eight target loci, we tested 75 noncoding segments by reporter gene assays in transiently transfected human K562 cells and/or after site-directed integration into murine erythroleukemia cells. Segments with a high RP score and a conserved exact match to the binding site consensus are validated at a good rate (50%-100%, with rates increasing at higher RP), whereas segments with lower RP scores or nonconsensus binding motifs tend to be inactive. Active DNA segments were shown to be occupied by GATA-1 protein by chromatin immunoprecipitation, whereas sites predicted to be inactive were not occupied. We verify four previously known erythroid CRMs and identify 28 novel ones. Thus, high RP in combination with another feature of a CRM, such as a conserved transcription factor binding site, is a good predictor of functional CRMs. Genome-wide predictions based on RP and a large set of well-defined transcription factor binding sites are available through servers at http://www.bx.psu.edu/.

ECRbase: database of evolutionary conserved regions, promoters, and transcription factor binding sites in vertebrate genomes

Evolutionary conservation of DNA sequences provides a tool for the identification of functional elements in genomes. We have created a database of evolutionary conserved regions (ECRs) in vertebrate genomes, entitled ECRbase, which is constructed from a collection of whole-genome alignments produced by the ECR Browser. ECRbase features a database of syntenic blocks that recapitulate the evolution of rearrangements in vertebrates and a comprehensive collection of promoters in all vertebrate genomes generated using multiple sources of gene annotation. The database also contains a collection of annotated transcription factor binding sites (TFBSs) in evolutionary conserved and promoter elements. ECRbase currently includes human, rhesus macaque, dog, opossum, rat, mouse, chicken, frog, zebrafish and fugu genomes. It is freely accessible at http://ecrbase.dcode.org.

In vivo enhancer analysis of human conserved non-coding sequences

Identifying the sequences that direct the spatial and temporal expression of genes and defining their function in vivo remains a significant challenge in the annotation of vertebrate genomes. One major obstacle is the lack of experimentally validated training sets. In this study, we made use of extreme evolutionary sequence conservation as a filter to identify putative gene regulatory elements, and characterized the in vivo enhancer activity of a large group of non-coding elements in the human genome that are conserved in human-pufferfish, Takifugu (Fugu) rubripes, or ultraconserved in human-mouse-rat. We tested 167 of these extremely conserved sequences in a transgenic mouse enhancer assay. Here we report that 45% of these sequences functioned reproducibly as tissue-specific enhancers of gene expression at embryonic day 11.5. While directing expression in a broad range of anatomical structures in the embryo, the majority of the 75 enhancers directed expression to various regions of the developing nervous system. We identified sequence signatures enriched in a subset of these elements that targeted forebrain expression, and used these features to rank all approximately 3,100 non-coding elements in the human genome that are conserved between human and Fugu. The testing of the top predictions in transgenic mice resulted in a threefold enrichment for sequences with forebrain enhancer activity. These data dramatically expand the catalogue of human gene enhancers that have been characterized in vivo, and illustrate the utility of such training sets for a variety of biological applications, including decoding the regulatory vocabulary of the human genome.

EPConDB: a web resource for gene expression related to pancreatic development, beta-cell function and diabetes

EPConDB () is a public web site that supports research in diabetes, pancreatic development and beta-cell function by providing information about genes expressed in cells of the pancreas. EPConDB displays expression profiles for individual genes and information about transcripts, promoter elements and transcription factor binding sites. Gene expression results are obtained from studies examining tissue expression, pancreatic development and growth, differentiation of insulin-producing cells, islet or beta-cell injury, and genetic models of impaired beta-cell function. The expression datasets are derived using different microarray platforms, including the BCBC PancChips and Affymetrix gene expression arrays. Other datasets include semi-quantitative RT–PCR and MPSS expression studies. For selected microarray studies, lists of differentially expressed genes, derived from PaGE analysis, are displayed on the site. EPConDB provides database queries and tools to examine the relationship between a gene, its transcriptional regulation, protein function and expression in pancreatic tissues.

Unraveling transcription regulatory networks by protein-DNA and protein-protein interaction mapping

Metazoan genomes contain thousands of protein-coding and noncoding RNA genes, most of which are differentially expressed, i.e., at different locations or at different times during development, function, or pathology of the organism. Differential gene expression is achieved in part by the action of regulatory transcription factors (TFs) that bind to cis-regulatory elements that are often located in or near their target genes. Each TF likely regulates many targets in the context of intricate transcription regulatory networks. Up to 10% of a genome may encode TFs, but only a handful of these have been studied in detail. Here, I will discuss the different steps involved in the mapping and analysis of transcription regulatory networks, including the identification of network nodes (TFs and their target sequences) and edges (TF-TF dimers and TF-DNA target interactions), integration with other data types, and network properties and emerging principles that provide insights into differential gene expression.

Functional variation and evolution of non-coding DNA

The focus of large genomic studies has shifted from only looking at genes and protein-coding sequences to exploring the full set of elements in each genome. The explosion of comparative sequencing data has led to an increase in methodologies, approaches and ideas on how to analyze the unknown fraction of the genome, namely the non-protein-coding fraction. The main issues relate to the discovery, evolutionary analysis and natural variation of non-coding DNA, and the parameters that prevent us from fully understanding the properties of non-coding DNA.

Forces shaping the fastest evolving regions in the human genome

Comparative genomics allow us to search the human genome for segments that were extensively changed in the last approximately 5 million years since divergence from our common ancestor with chimpanzee, but are highly conserved in other species and thus are likely to be functional. We found 202 genomic elements that are highly conserved in vertebrates but show evidence of significantly accelerated substitution rates in human. These are mostly in non-coding DNA, often near genes associated with transcription and DNA binding. Resequencing confirmed that the five most accelerated elements are dramatically changed in human but not in other primates, with seven times more substitutions in human than in chimp. The accelerated elements, and in particular the top five, show a strong bias for adenine and thymine to guanine and cytosine nucleotide changes and are disproportionately located in high recombination and high guanine and cytosine content environments near telomeres, suggesting either biased gene conversion or isochore selection. In addition, there is some evidence of directional selection in the regions containing the two most accelerated regions. A combination of evolutionary forces has contributed to accelerated evolution of the fastest evolving elements in the human genome.

A global genomic transcriptional code associated with CNS-expressed genes

Highly conserved non-coding DNA regions (HCNR) occur frequently in vertebrate genomes, but their functional roles remain unclear. Here, we provide evidence that a large portion of HCNRs are enriched for binding sites for Sox, POU and Homeodomain transcription factors, and such HCNRs can act as cis-regulatory regions active in neural stem cells. Strikingly, these HCNRs are linked to several hundreds of genes expressed in the developing CNS and they may exert locus-wide regulatory effects on multiple genes flanking their genomic location. Moreover, these data imply a unifying transcriptional logic for a large set of CNS-expressed genes in which Sox and POU proteins act as generic promoters of transcription while Homeodomain proteins control the spatial expression of genes through active repression.

Mammalian ultraconserved elements are strongly depleted among segmental duplications and copy number variants

An earlier search in the human, mouse and rat genomes for sequences that are 100% conserved in orthologous segments and > or = 200 bp in length identified 481 distinct sequences. These human-mouse-rat sequences, which represent ultraconserved elements (UCEs), are believed to be important for functions involving DNA binding, RNA processing and the regulation of transcription and development. In vivo and additional computational studies of UCEs and other highly conserved sequences are consistent with these functional associations, with some observations indicating enhancer-like activity for these elements. Here, we show that UCEs are significantly depleted among segmental duplications and copy number variants. Notably, of the UCEs that are found in segmental duplications or copy number variants, the majority overlap exons, indicating, along with other findings presented, that UCEs overlapping exons represent a distinct subset.

Shared long-range regulatory elements coordinate expression of a gene cluster encoding nicotinic receptor heteromeric subtypes

The nicotinic acetylcholine receptor (nAChR) beta4/alpha3/alpha5 gene cluster encodes several heteromeric transmitter receptor subtypes that are essential for cholinergic synaptic transmission in adrenal gland, autonomic ganglia, pineal gland, and several nuclei in the central nervous system. However, the transcriptional mechanisms coordinating expression of these subunit genes in different cell populations are unknown. Here, we used transgenic methods to investigate long-range transcriptional control of the cluster. A 132-kb P1-derived artificial chromosome (PAC) encoding the rat cluster recapitulated the neurally- and endocrine-restricted expression patterns of the endogenous beta4/alpha3/alpha5 genes. Mutation of ETS factor binding sites in an enhancer, beta43', embedded in the beta4 3'-untranslated exon resulted in greatly diminished beta4, alpha3, and alpha5 expression in adrenal gland and to a lesser extent in the superior cervical ganglion (SCG) but not in other tissues. Phylogenetic sequence analyses revealed several conserved noncoding regions (CNRs) upstream of beta4 and alpha5. Deletion of one of them (CNR4) located 20 kb upstream of beta4 resulted in a dramatic decrease in beta4 and alpha3 expression in the pineal gland and SCG. CNR4 was sufficient to direct LacZ transgene expression to SCG neurons, which express the endogenous beta4alpha3alpha5 subunits, and pineal cells, which express the endogenous beta4alpha3 combination. Finally, CNR4 was able to direct transgene expression to major sites of expression of the endogenous cluster in the brain. Together, our findings support a model in which cell type-specific shared long-range regulatory elements are required for coordinate expression of clustered nAChR genes.

Genomic data analysis using DNA structure: an analysis of conserved nongenic sequences and ultraconserved elements

Recent comparative studies of the human and mouse genomes have revealed sets of conserved nongenic sequences (CNGs) and sets of ultraconserved elements (UCEs). Both sets of sequences, which exhibit extremely high levels of conservation, extend over hundreds of bases and have no known function. Since there is no detectable sequence homology between paralogous CNGs or UCEs in either of the species, an alignment-free technique is needed for their analysis. We have previously compiled a database of the structural properties of all 32,896 unique DNA octamers, including information on stability, the minimum energy conformation, and flexibility. We have used Fourier techniques to analyze the UCEs and CNGs in terms of their octamer structural properties, to reveal structural correlations which may indicate possible functions for some of these sequences.

DMRT gene cluster analysis in the platypus: new insights into genomic organization and regulatory regions

We isolated and characterized a cluster of platypus DMRT genes and compared their arrangement, location, and sequence across vertebrates. The DMRT gene cluster on human 9p24.3 harbors, in order, DMRT1, DMRT3, and DMRT2, which share a DM domain. DMRT1 is highly conserved and involved in sexual development in vertebrates, and deletions in this region cause sex reversal in humans. Sequence comparisons of DMRT genes between species have been valuable in identifying exons, control regions, and conserved nongenic regions (CNGs). The addition of platypus sequences is expected to be particularly valuable, since monotremes fill a gap in the vertebrate genome coverage. We therefore isolated and fully sequenced platypus BAC clones containing DMRT3 and DMRT2 as well as DMRT1 and then generated multispecies alignments and ran prediction programs followed by experimental verification to annotate this gene cluster. We found that the three genes have 58-66% identity to their human orthologues, lie in the same order as in other vertebrates, and colocate on 1 of the 10 platypus sex chromosomes, X5. We also predict that optimal annotation of the newly sequenced platypus genome will be challenging. The analysis of platypus sequence revealed differences in structure and sequence of the DMRT gene cluster. Multispecies comparison was particularly effective for detecting CNGs, revealing several novel potential regulatory regions within DMRT3 and DMRT2 as well as DMRT1. RT-PCR indicated that platypus DMRT1 and DMRT3 are expressed specifically in the adult testis (and not ovary), but DMRT2 has a wider expression profile, as it does for other mammals. The platypus DMRT1 expression pattern, and its location on an X chromosome, suggests an involvement in monotreme sexual development.

A genetic linkage map of the hermaphrodite teleost fish Sparus aurata L

The gilthead sea bream (Sparus aurata L.) is a marine fish of great importance for fisheries and aquaculture. It has also a peculiar sex-determination system, being a protandrous hermaphrodite. Here we report the construction of a first-generation genetic linkage map for S. aurata, based on 204 microsatellite markers. Twenty-six linkage groups (LG) were found. The total map length was 1241.9 cM. The ratio between sex-specific map lengths was 1:1.2 (male:female). Comparison with a preliminary radiation hybrid (RH) map reveals a good concordance, as all markers located in a single LG are located in a single RH group, except for Ad-25 and CId-31. Comparison with the Tetraodon nigroviridis genome revealed a considerable number of evolutionary conserved regions (ECRs) between the two species. The mean size of ECRs was 182 bp (sequence identity 60-90%). Forty-one ECRs have a known chromosomal location in the pufferfish genome. Despite the limited number of anchoring points, significant syntenic relationships were found. The linkage map presented here provides a robust comparative framework for QTL analysis in S. aurata and is a step toward the identification of genetic loci involved both in the determination of economically important traits and in the individual timing of sex reversal.

Genomic organization, alternative splicing, and multiple regulatory regions of the zebrafish fgf8 gene

Fgf8 is among the members of the fibroblast growth factor (FGF) family that play pivotal roles in vertebrate development. In the present study, the genomic DNA of the zebrafish fgf8 gene was cloned to elucidate the regulatory mechanism behind the temporally and spatially restricted expression of the gene in vertebrate embryos. Structural analysis revealed that the exon-intron organization of fgf8 is highly conserved during vertebrate evolution, from teleosts to mammals. Close inspection of the genomic sequence and reverse transcription-polymerase chain reaction analysis revealed that zebrafish fgf8 encodes two splicing variants, corresponding to Fgf8a and Fgf8b, among the four to seven splicing variants known in mammals. Misexpression of the two variants in zebrafish embryos following mRNA injection showed that both variants have dorsalizing activities on zebrafish embryos, with Fgf8b being more potent. Reporter gene analysis of the transcriptional regulation of zebrafish fgf8 suggested that its complicated expression pattern, which is considered essential for its multiple roles in development, is mediated by combinations of different regulatory regions in the upstream and downstream regions of the gene. Furthermore, comparison of the genomic sequence of fgf8 among different vertebrate species suggests that this regulatory mechanism is conserved during vertebrate evolution.

Conserved distances between vertebrate highly conserved elements

High numbers of sequence element with very high (>95%) sequence conservation between the human and other vertebrate genomes have been reported and ascribed putative cis-regulatory functions. We have investigated the structural relationships between such elements in mammalian genomes and find that not only their sequences, but also the distances between them are significantly (P<2.2x10(-16)) more conserved than corresponding distances between orthologous protein-coding genes or between exons within these genes. Regions of largely conserved distance between consecutive highly conserved elements (HCE) generally overlap previously identified HCE clusters, but may be far longer (up to 20 Mb) and possibly cover close to 25% of the human genome sequence. Similar conservation of distance is found between bird (chicken) and mammalian genomes and is also discernible in comparisons between fish and mammals. The data suggest either that a substantial amount of essential (functionally active) elements with lower sequence conservation occupy the space between the HCEs or that distance itself is an important factor in transcriptional regulation or chromatin modelling.

Scale-invariant structure of strongly conserved sequence in genomic intersections and alignments

A power-law distribution of the length of perfectly conserved sequence from mouse/human whole-genome intersection and alignment is exhibited. Spatial correlations of these elements within the mouse genome are studied. It is argued that these power-law distributions and correlations are comprised in part by functional noncoding sequence and ought to be accounted for in estimating the statistical significance of apparent sequence conservation. These inter-genomic correlations of conservation are placed in the context of previously observed intra-genomic correlations, and their possible origins and consequences are discussed.

Comparative evolutionary genomics of the HADH2 gene encoding Abeta-binding alcohol dehydrogenase/17beta-hydroxysteroid dehydrogenase type 10 (ABAD/HSD10)

Background

The Aβ-binding alcohol dehydrogenase/17β-hydroxysteroid dehydrogenase type 10 (ABAD/HSD10) is an enzyme involved in pivotal metabolic processes and in the mitochondrial dysfunction seen in the Alzheimer's disease. Here we use comparative genomic analyses to study the evolution of the HADH2 gene encoding ABAD/HSD10 across several eukaryotic species.

Results

Both vertebrate and nematode HADH2 genes showed a six-exon/five-intron organization while those of the insects had a reduced and varied number of exons (two to three). Eutherian mammal HADH2 genes revealed some highly conserved noncoding regions, which may indicate the presence of functional elements, namely in the upstream region about 1 kb of the transcription start site and in the first part of intron 1. These regions were also conserved between Tetraodon and Fugu fishes. We identified a conserved alternative splicing event between human and dog, which have a nine amino acid deletion, causing the removal of the strand β_F. This strand is one of the seven strands that compose the core β-sheet of the Rossman fold dinucleotide-binding motif characteristic of the short chain dehydrogenase/reductase (SDR) family members. However, the fact that the substrate binding cleft residues are retained and the existence of a shared variant between human and dog suggest that it might be functional. Molecular adaptation analyses across eutherian mammal orthologues revealed the existence of sites under positive selection, some of which being localized in the substrate-binding cleft and in the insertion 1 region on loop D (an important region for the Aβ-binding to the enzyme). Interestingly, a higher than expected number of nonsynonymous substitutions were observed between human/chimpanzee and orangutan, with six out of the seven amino acid replacements being under molecular adaptation (including three in loop D and one in the substrate binding loop).

Conclusion

Our study revealed that HADH2 genes maintained a reasonable conserved organization across a large evolutionary distance. The conserved noncoding regions identified among mammals and between pufferfishes, the evidence of an alternative splicing variant conserved between human and dog, and the detection of positive selection across eutherian mammals, may be of importance for further research on ABAD/HSD10 function and its implication in the Alzheimer's disease.

Regulation of hypocretin (orexin) expression in embryonic zebrafish

Hypocretins/orexins are neuropeptides involved in the regulation of sleep and energy balance in mammals. Conservation of gene sequence, hypothalamic localization of cell bodies, and projection patterns in adult zebrafish suggest that the architecture and function of the hypocretin system are conserved in fish. We report on the complete genomic structure of the zebrafish and Tetraodon hypocretin genes and the complete predicted hypocretin protein sequences from five teleosts. Using whole mount in situ hybridization, we have traced the development of hypocretin cells in zebrafish from onset of expression at 22 h post-fertilization through the first week of development. Promoter elements of similar size from zebrafish and Tetraodon were capable of driving efficient and specific expression of enhanced green fluorescent protein in developing zebrafish embryos, thus defining a minimal promoter region able to accurately mimic the native hypocretin pattern. This enhanced green fluorescent protein expression also revealed a complex pattern of projections within the hypothalamus, to the midbrain, and to the spinal cord. To further analyze the promoter, a series of deletion and substitution constructs were injected into embryos, and resulting promoter activity was monitored in the first week of development. A critical region of 250 base pairs was identified containing a core 13-base pair element essential for hypocretin expression.

A family of conserved noncoding elements derived from an ancient transposable element

The evolutionary origin of the conserved noncoding elements (CNEs) in the human genome remains poorly understood but may hold important clues to their biological functions. Here, we report the discovery of a CNE family with approximately 124 instances in the human genome that demonstrates a clear signature of having been derived from an ancient transposon. The CNE family is also present in the chicken genome, although typically not at orthologous locations. The CNE family is closely related to the active transposon SINE3 in zebrafish and also to a previously uncharacterized transposon in the coelacanth, the so-called "living fossil" belonging to the lobe-finned fish lineage. The mammal, bird, zebrafish, and coelacanth families all share a highly similar core element of approximately 180 bp but have important differences in their 5' and 3' ends. The core element has thus been preserved over 450 million years of evolution, implying an important biological function. In addition, we identify 95 additional CNE families that likely predate the mammalian radiation. The results highlight both the creative role of transposons and the importance of CNE families.

New technologies, new findings, and new concepts in the study of vertebrate cis-regulatory sequences

All vertebrates share a similar early embryonic body plan and use the same regulatory genes for their development. The availability of numerous sequenced vertebrate genomes and significant advances in bioinformatics have resulted in the finding that the genomic regions of many of these developmental regulatory genes also contain highly conserved noncoding sequence. In silico discovery of conserved noncoding regions and of transcription factor binding sites as well as the development of methods for high throughput transgenesis in Xenopus and zebrafish are dramatically increasing the speed with which regulatory elements can be discovered, characterized, and tested in the context of whole live embryos. We review here some of the recent technological developments that will likely lead to a surge in research on how vertebrate genomes encode regulation of transcriptional activity, how regulatory sequences constrain genomic architecture, and ultimately how vertebrate form has evolved.

Control of developmental regulators by Polycomb in human embryonic stem cells

Polycomb group proteins are essential for early development in metazoans, but their contributions to human development are not well understood. We have mapped the Polycomb Repressive Complex 2 (PRC2) subunit SUZ12 across the entire nonrepeat portion of the genome in human embryonic stem (ES) cells. We found that SUZ12 is distributed across large portions of over two hundred genes encoding key developmental regulators. These genes are occupied by nucleosomes trimethylated at histone H3K27, are transcriptionally repressed, and contain some of the most highly conserved noncoding elements in the genome. We found that PRC2 target genes are preferentially activated during ES cell differentiation and that the ES cell regulators OCT4, SOX2, and NANOG cooccupy a significant subset of these genes. These results indicate that PRC2 occupies a special set of developmental genes in ES cells that must be repressed to maintain pluripotency and that are poised for activation during ES cell differentiation.

Searching chromatin for stem cell identity

Stem cells encapsulate the fundamental problem of metazoan biology in miniature: How do cells establish and maintain their fates? Increasing evidence indicates that stem cell chromatin activates proliferation genes and represses differentiation genes. Understanding how these configurations are stabilized by Polycomb group proteins will advance our understanding of embryonic development, tissue homeostasis, regeneration, aging, and oncogenesis.

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.

Close sequence comparisons are sufficient to identify human cis-regulatory elements

Cross-species DNA sequence comparison is the primary method used to identify functional noncoding elements in human and other large genomes. However, little is known about the relative merits of evolutionarily close and distant sequence comparisons. To address this problem, we identified evolutionarily conserved noncoding regions in primate, mammalian, and more distant comparisons using a uniform approach (Gumby) that facilitates unbiased assessment of the impact of evolutionary distance on predictive power. We benchmarked computational predictions against previously identified cis-regulatory elements at diverse genomic loci and also tested numerous extremely conserved human-rodent sequences for transcriptional enhancer activity using an in vivo enhancer assay in transgenic mice. Human regulatory elements were identified with acceptable sensitivity (53%-80%) and true-positive rate (27%-67%) by comparison with one to five other eutherian mammals or six other simian primates. More distant comparisons (marsupial, avian, amphibian, and fish) failed to identify many of the empirically defined functional noncoding elements. Our results highlight the practical utility of close sequence comparisons, and the loss of sensitivity entailed by more distant comparisons. We derived an intuitive relationship between ancient and recent noncoding sequence conservation from whole-genome comparative analysis that explains most of the observations from empirical benchmarking. Lastly, we determined that, in addition to strength of conservation, genomic location and/or density of surrounding conserved elements must also be considered in selecting candidate enhancers for in vivo testing at embryonic time points.

Autism and ultraconserved non-coding sequence on chromosome 7q

OBJECTIVE

Autism has been linked to a broad region on chromosome 7q that contains a large number of genes involved in transcription and development. This region is also enriched for ultraconserved non-coding elements, defined as human-rodent sequences that are 100% aligned over > or =200 base pairs, which have a high likelihood of being functional. Therefore, as only a few rare coding variants have been detected in the autism candidate genes on 7q examined to date, we decided to screen these ultraconserved elements for possible autism susceptibility alleles.

METHODS

We used denaturing high-performance liquid chromatography, and DNA sequencing, to perform variant detection in a total of 146 cases with autism, 96 from the Autism Genetic Resource Exchange and 50 from the Central Valley of Costa Rica, as well as 124 controls from the Polymorphism Discovery Resource Panel. We screened 10 consecutive ultraconserved elements in, or flanking, the genes DLX5/6, AUTS2 and FOXP2 on chromosome 7q.

RESULTS

Although we did find several rare variants in autism cases that were not present in controls, we also observed rare variants present in controls and not cases. The most common variant occurred in controls at a frequency of 3.3%. Interestingly, two ultraconserved elements each harbored three independent variants and one ultraconserved element harbored two independent variants, suggesting that ultraconservation is maintained chiefly by a decreased tendency toward fixation, rather than a significantly lower mutation rate.

CONCLUSIONS

Our results show that these sequences are unlikely to harbor major autism susceptibility alleles.

The Evf-2 noncoding RNA is transcribed from the Dlx-5/6 ultraconserved region and functions as a Dlx-2 transcriptional coactivator

The identification of ultraconserved noncoding sequences in vertebrates has been associated with developmental regulators and DNA-binding proteins. One of the first of these was identified in the intergenic region between the Dlx-5 and Dlx-6 genes, members of the Dlx/dll homeodomain-containing protein family. In previous experiments, we showed that Sonic hedgehog treatment of forebrain neural explants results in the activation of Dlx-2 and the novel noncoding RNA (ncRNA), Evf-1. In this report, we show that the Dlx-5/6 ultraconserved region is transcribed to generate an alternatively spliced form of Evf-1, the ncRNA Evf-2. Evf-2 specifically cooperates with Dlx-2 to increase the transcriptional activity of the Dlx-5/6 enhancer in a target and homeodomain-specific manner. A stable complex containing the Evf-2 ncRNA and the Dlx-2 protein forms in vivo, suggesting that the Evf-2 ncRNA activates transcriptional activity by directly influencing Dlx-2 activity. These experiments identify a novel mechanism whereby transcription is controlled by the cooperative actions of an ncRNA and a homeodomain protein. The possibility that a subset of vertebrate ultraconserved regions may function at both the DNA and RNA level to control key developmental regulators may explain why ultraconserved sequences exhibit 90% or more conservation even after 450 million years of vertebrate evolution.

Short blocks from the noncoding parts of the human genome have instances within nearly all known genes and relate to biological processes

Using an unsupervised pattern-discovery method, we processed the human intergenic and intronic regions and catalogued all variable-length patterns with identically conserved copies and multiplicities above what is expected by chance. Among the millions of discovered patterns, we found a subset of 127,998 patterns, termed pyknons, which have additional nonoverlapping instances in the untranslated and protein-coding regions of 30,675 transcripts from 20,059 human genes. The pyknons arrange combinatorially in the untranslated and coding regions of numerous human genes where they form mosaics. Consecutive instances of pyknons in these regions show a strong bias in their relative placement, favoring distances of approximately 22 nucleotides. We also found pyknons to be enriched in a statistically significant manner in genes involved in specific processes, e.g., cell communication, transcription, regulation of transcription, signaling, transport, etc. For approximately 1/3 of the pyknons, the intergenic/intronic instances of their reverse complement lie within 380,084 nonoverlapping regions, typically 60-80 nucleotides long, which are predicted to form double-stranded, energetically stable, hairpin-shaped RNA secondary structures; additionally, the pyknons subsume approximately 40% of the known microRNA sequences, thus suggesting a possible link with posttranscriptional gene silencing and RNA interference. Cross-genome comparisons reveal that many of the pyknons have instances in the 3' UTRs of genes from other vertebrates and invertebrates where they are overrepresented in similar biological processes, as in the human genome. These unexpected findings suggest potential unique functional connections between the coding and noncoding parts of the human genome.

Conserved co-regulation and promoter sharing of hoxb3a and hoxb4a in zebrafish

The expression of zebrafish hoxb3a and hoxb4a has been found to be mediated through five transcripts, hoxb3a transcripts I-III and hoxb4a transcripts I-II, driven by four promoters. A "master" promoter, located about 2 kb downstream of hoxb5a, controls transcription of a pre-mRNA comprising exon sequences of both genes. This unique gene structure is proposed to provide a novel mechanism to ensure overlapping, tissue-specific expression of both genes in the posterior hindbrain and spinal cord. Transgenic approaches were used to analyze the functions of zebrafish hoxb3a/hoxb4a promoters and enhancer sequences containing regions of homology that were previously identified by comparative genomics. Two neural enhancers were shown to establish specific anterior expression borders within the hindbrain and mediate expression in defined neuronal populations derived from hindbrain rhombomeres (r) 5 to 8, suggesting a late role of the genes in neuronal cell lineage specification. Species comparison showed that the zebrafish hoxb3a r5 and r6 enhancer corresponded to a sequence within the mouse HoxA cluster controlling activity of Hoxa3 in r5 and r6, whereas a homologous region within the HoxB cluster activated Hoxb3 expression but limited to r5. We conclude that the similarity of hoxb3a/Hoxa3 regulatory mechanisms reflect the shared descent of both genes from a single ancestral paralog group 3 gene.

A distal enhancer and an ultraconserved exon are derived from a novel retroposon

Hundreds of highly conserved distal cis-regulatory elements have been characterized so far in vertebrate genomes. Many thousands more are predicted on the basis of comparative genomics. However, in stark contrast to the genes that they regulate, in invertebrates virtually none of these regions can be traced by using sequence similarity, leaving their evolutionary origins obscure. Here we show that a class of conserved, primarily non-coding regions in tetrapods originated from a previously unknown short interspersed repetitive element (SINE) retroposon family that was active in the Sarcopterygii (lobe-finned fishes and terrestrial vertebrates) in the Silurian period at least 410 million years ago (ref. 4), and seems to be recently active in the 'living fossil' Indonesian coelacanth, Latimeria menadoensis. Using a mouse enhancer assay we show that one copy, 0.5 million bases from the neuro-developmental gene ISL1, is an enhancer that recapitulates multiple aspects of Isl1 expression patterns. Several other copies represent new, possibly regulatory, alternatively spliced exons in the middle of pre-existing Sarcopterygian genes. One of these, a more than 200-base-pair ultraconserved region, 100% identical in mammals, and 80% identical to the coelacanth SINE, contains a 31-amino-acid-residue alternatively spliced exon of the messenger RNA processing gene PCBP2 (ref. 6). These add to a growing list of examples in which relics of transposable elements have acquired a function that serves their host, a process termed 'exaptation', and provide an origin for at least some of the many highly conserved vertebrate-specific genomic sequences.

Identification of cis-regulatory elements for MECP2 expression

Rett syndrome (RTT) is an X-linked dominant disabling neurodevelopmental disorder caused by loss of function mutations in the MECP2 gene, located at Xq28, which encodes a multifunctional protein. MECP2 expression is regulated in a developmental stage and cell-type-specific manner. The need for tightly controlled MeCP2 levels in brain is strongly suggested by neurologically abnormal phenotypes of mouse models with mild overexpression and by mental retardation in human males with MECP2 duplication. We set out to identify long-range cis-regulatory sequences that differentially regulate MECP2 transcription and, when mutated, may contribute to the pathogenesis of RTT, autism or X-linked mental retardation. By inter-species sequence comparisons, we detected 27 highly conserved non-coding DNA sequences within a 210 kb region covering MECP2. We functionally confirmed four enhancer and two silencer elements by performing luciferase reporter assays in four different human cell lines. The transcription factor binding capability of the identified regulatory elements was tested by gel shift assays. To locate the human MECP2 core promoter, we dissected the promoter region by reporter assays with deletion constructs. We then used chromosome conformation capture methods to document long-range interactions of three enhancers and two silencers with the MECP2 promoter. Acting over distances of up to 130 kb, these elements may influence chromatin configurations and regulate MECP2 transcription. Our study has defined the "MECP2 functional expression module" and identified enhancer and silencer elements that are likely to be responsible for the tissue-specific, developmental stage-specific or splice-variant-specific control of MeCP2 protein expression.

Evolutionary genomics of the recently duplicated amphioxus Hairy genes

Amphioxus Hairy genes have gone through a number of lineage-specific duplications, resulting in eight members, some of which are differentially expressed in the embryo. In order to gain insights into the evolution and function of this gene family we have compared their genomic structure and searched for conserved non-coding sequence elements. We have found that introns have been lost independently from these genes at least twice and after the duplication events. By carrying out phylogenetic footprinting between paralogues expressed in the embryo, we have found a differential distribution of conserved elements that could explain the limited overlap in expression patterns of Hairy genes in the amphioxus embryo. Furthermore, clustering of RBP-Jk binding sites in these conserved elements suggests that amphioxus Hairy genes are downstream targets of the Notch signaling pathway, as occurs in vertebrates. All of this evidence suggests that amphioxus Hairy genes have gone through a process of subfunctionalization shortly after their duplication, representing an extreme and rapid case of the duplication-degeneration-complementation model.

Genome-wide computational prediction of transcriptional regulatory modules reveals new insights into human gene expression

The identification of regulatory regions is one of the most important and challenging problems toward the functional annotation of the human genome. In higher eukaryotes, transcription-factor (TF) binding sites are often organized in clusters called cis-regulatory modules (CRM). While the prediction of individual TF-binding sites is a notoriously difficult problem, CRM prediction has proven to be somewhat more reliable. Starting from a set of predicted binding sites for more than 200 TF families documented in Transfac, we describe an algorithm relying on the principle that CRMs generally contain several phylogenetically conserved binding sites for a few different TFs. The method allows the prediction of more than 118,000 CRMs within the human genome. A subset of these is shown to be bound in vivo by TFs using ChIP-chip. Their analysis reveals, among other things, that CRM density varies widely across the genome, with CRM-rich regions often being located near genes encoding transcription factors involved in development. Predicted CRMs show a surprising enrichment near the 3' end of genes and in regions far from genes. We document the tendency for certain TFs to bind modules located in specific regions with respect to their target genes and identify TFs likely to be involved in tissue-specific regulation. The set of predicted CRMs, which is made available as a public database called PReMod (http://genomequebec.mcgill.ca/PReMod), will help analyze regulatory mechanisms in specific biological systems.

A comparative genomics strategy for targeted discovery of single-nucleotide polymorphisms and conserved-noncoding sequences in orphan crops

Completed genome sequences provide templates for the design of genome analysis tools in orphan species lacking sequence information. To demonstrate this principle, we designed 384 PCR primer pairs to conserved exonic regions flanking introns, using Sorghum/Pennisetum expressed sequence tag alignments to the Oryza genome. Conserved-intron scanning primers (CISPs) amplified single-copy loci at 37% to 80% success rates in taxa that sample much of the approximately 50-million years of Poaceae divergence. While the conserved nature of exons fostered cross-taxon amplification, the lesser evolutionary constraints on introns enhanced single-nucleotide polymorphism detection. For example, in eight rice (Oryza sativa) genotypes, polymorphism averaged 12.1 per kb in introns but only 3.6 per kb in exons. Curiously, among 124 CISPs evaluated across Oryza, Sorghum, Pennisetum, Cynodon, Eragrostis, Zea, Triticum, and Hordeum, 23 (18.5%) seemed to be subject to rigid intron size constraints that were independent of per-nucleotide DNA sequence variation. Furthermore, we identified 487 conserved-noncoding sequence motifs in 129 CISP loci. A large CISP set (6,062 primer pairs, amplifying introns from 1,676 genes) designed using an automated pipeline showed generally higher abundance in recombinogenic than in nonrecombinogenic regions of the rice genome, thus providing relatively even distribution along genetic maps. CISPs are an effective means to explore poorly characterized genomes for both DNA polymorphism and noncoding sequence conservation on a genome-wide or candidate gene basis, and also provide anchor points for comparative genomics across a diverse range of species.

Ancient duplicated conserved noncoding elements in vertebrates: a genomic and functional analysis

Fish-mammal genomic comparisons have proved powerful in identifying conserved noncoding elements likely to be cis-regulatory in nature, and the majority of those tested in vivo have been shown to act as tissue-specific enhancers associated with genes involved in transcriptional regulation of development. Although most of these elements share little sequence identity to each other, a small number are remarkably similar and appear to be the product of duplication events. Here, we searched for duplicated conserved noncoding elements in the human genome, using comparisons with Fugu to select putative cis-regulatory sequences. We identified 124 families of duplicated elements, each containing between two and five members, that are highly conserved within and between vertebrate genomes. In 74% of cases, we were able to assign a specific set of paralogous genes with annotation relating to transcriptional regulation and/or development to each family, thus removing much of the ambiguity in identifying associated genes. We find that duplicate elements have the potential to up-regulate reporter gene expression in a tissue-specific manner and that expression domains often overlap, but are not necessarily identical, between family members. Over two thirds of the families are conserved in duplicate in fish and appear to predate the large-scale duplication events thought to have occurred at the origin of vertebrates. We propose a model whereby gene duplication and the evolution of cis-regulatory elements can be considered in the context of increased morphological diversity and the emergence of the modern vertebrate body plan.

Genome-wide prediction of mammalian enhancers based on analysis of transcription-factor binding affinity

Understanding the regulation of human gene expression requires knowledge of the "second genetic code," which consists of the binding specificities of transcription factors (TFs) and the combinatorial code by which TF binding sites are assembled to form tissue-specific enhancer elements. Using a novel high-throughput method, we determined the DNA binding specificities of GLIs 1-3, Tcf4, and c-Ets1, which mediate transcriptional responses to the Hedgehog (Hh), Wnt, and Ras/MAPK signaling pathways. To identify mammalian enhancer elements regulated by these pathways on a genomic scale, we developed a computational tool, enhancer element locator (EEL). We show that EEL can be used to identify Hh and Wnt target genes and to predict activated TFs based on changes in gene expression. Predictions validated in transgenic mouse embryos revealed the presence of multiple tissue-specific enhancers in mouse c-Myc and N-Myc genes, which has implications for organ-specific growth control and tumor-type specificity of oncogenes.

HOXD10 M319K mutation in a family with isolated congenital vertical talus

Congenital vertical talus (CVT) is a primary dislocation of the talonavicular joint that often occurs in neuromusculoskeletal syndromes, but may also be seen as an isolated abnormality. Six families with isolated CVT were ascertained. DNA was isolated from 21 affected individuals and 17 unaffected individuals from these families, as well as from five sporadic patients with CVT. Variable expressivity was noted in three families, manifesting as clubfoot in three individuals. Genome-wide linkage analysis generated a maximum two-point logarithm of odds score on chromosome 2q with D2S1353 (Zmax = 1.43 at theta(max) = 0.1), 17 Mb from the HOXD gene cluster. DNA from one affected individual of each family was subjected to mutational analysis of the HOXD10 gene. A single missense mutation was identified (M319K, 956T > A) in the homeodomain recognition helix of the HOXD10 gene that segregated with disease in one large British family. This mutation was recently described in a family of Italian descent with CVT and Charcot-Marie-Tooth deformity HOXD10 gene mutations were not identified in any of the other families or sporadic patients with CVT, suggesting that genetic heterogeneity underlies this disorder.

A large family of ancient repeat elements in the human genome is under strong selection

Although conserved noncoding elements (CNEs) constitute the majority of sequences under purifying selection in the human genome, they remain poorly understood. CNEs seem to be largely unique, with no large families of similar elements reported to date. Here, we search for CNEs among the ancestral repeat classes in the human genome and report the discovery of a large CNE family containing >900 members. This family belongs to the MER121 class of repeats. Although the MER121 family members show considerable sequence variation among one another, the individual copies show striking conservation in orthologous locations across the human, dog, mouse, and rat genomes. The element is also present and conserved in orthologous locations in the marsupial, but its genome-wide dispersal postdates the divergence from birds. The comparative genomic data indicate that MER121 does not encode a family of either protein-coding or RNA genes. Although the precise function of these elements remains unknown, the evidence suggests that this unusual family may play a cis-regulatory or structural role in mammalian genomes.

Fish genomics and biology

The last common ancestor between fish and mammals dates back to the very origin of the vertebrate lineage and today, half of modern vertebrates are fish. It is thus not surprising that several fish species have played important roles in recent years to advance our understanding of vertebrate genome evolution, to inform us on the structure of human genes, and, somewhat more unexpectedly, to provide leads to understanding the function of genes involved in human diseases. Genome sequence comparisons between such distantly related organisms are highly informative due to the accumulation of neutral mutations in nonfunctional regions. Yet humans and fishes share many developmental pathways, organ systems, and physiological mechanisms, making conclusions relevant to human biology. The respective advantages of zebrafish, medaka, Tetraodon, or Takifugu have been well exploited so far with bioinformatics analyses and molecular biology techniques. However the full potential of fish genomics is about to be unleashed with the integration of more traditional disciplines such as biochemistry and physiology, with the study of additional species such as carp, trout, or tilapia and a broadening of its applications to environmental genomics or aquaculture.

Long-range upstream and downstream enhancers control distinct subsets of the complex spatiotemporal Sox9 expression pattern

SOX9 is an evolutionary conserved transcription factor that is expressed in a variety of tissues, with essential functions in cartilage, testis, heart, glial cell, inner ear and neural crest development. By comparing human and pufferfish genomic sequences, we previously identified eight highly conserved sequence elements between 290 kb 5' and 450 kb 3' to human SOX9. In this study, we assayed the regulatory potential of elements E1 to E7 in transgenic mice using a lacZ reporter gene driven by a 529 bp minimal mouse Sox9 promoter. We found that three of these elements and the Sox9 promoter control distinct subsets of the tissue-specific expression pattern of Sox9. E3, located 251 kb 5' to SOX9, directs lacZ expression to cranial neural crest cells and to the inner ear. E1 is located 28 kb 5' to SOX9 and controls expression in the node, notochord, gut, bronchial epithelium and pancreas. Transgene expression in the neuroectoderm is mediated by E7, located 95 kb 3' to SOX9, which regulates expression in the telencephalon and midbrain, and by the Sox9 minimal promoter which controls expression in the ventral spinal cord and hindbrain. We show that E3-directed reporter gene expression in neural crest cells of the first but not of the second and third pharyngeal arch is dependent on beta-catenin, revealing a complex regulation of Sox9 in cranial neural crest cells. Moreover, we identify and discuss highly conserved transcription factor binding sites within enhancer E3 that are in good agreement with current models for neural crest and inner ear development. Finally, we identify enhancer E1 as a cis-regulatory element conserved between vertebrates and invertebrates, indicating that some cis-regulatory sequences that control developmental genes in vertebrates might be phylogenetically ancient.

Highly efficient transgenesis in Xenopus tropicalis using I-SceI meganuclease

In this study, we report a highly efficient transgenesis technique for Xenopus tropicalis based on a method described first for Medaka. This simple procedure entails co-injection of meganuclease I-SceI and a transgene construct flanked by two I-SceI sites into fertilized eggs. Approximately 30% of injected embryos express transgenes in a promoter-dependent manner. About 1/3 of such embryos show incorporation of the transgene at the one-cell stage and the remainder are 'half-transgenics' suggesting incorporation at the two-cell stage. Transgenes from both classes of embryos are shown to be transmitted and expressed in offspring. The procedure also works efficiently in Xenopus laevis. Because the needle injection procedure does not significantly damage embryos, a high fraction develop normally and can, as well, be injected with a second reagent, for example an mRNA or antisense morpholino oligonucleotide, thus allowing one to perform several genetic manipulations on embryos at one time. This simple and efficient technique will be a powerful tool for high-throughput transgenesis assays in founder animals, and for facilitating genetic studies in the fast-breeding diploid frog, X. tropicalis.

Evolutionary constraints in conserved nongenic sequences of mammals

Mammalian genomes contain many highly conserved nongenic sequences (CNGs) whose functional significance is poorly understood. Sets of CNGs have previously been identified by selecting the most conserved elements from a chromosome or genome, but in these highly selected samples, conservation may be unrelated to purifying selection. Furthermore, conservation of CNGs may be caused by mutation rate variation rather than selective constraints. To account for the effect of selective sampling, we have examined conservation of CNGs in taxa whose evolution is largely independent of the taxa from which the CNGs were initially identified, and we have controlled for mutation rate variation in the genome. We show that selective constraints in CNGs and their flanks are about one-half as strong in hominids as in murids, implying that hominids have accumulated many slightly deleterious mutations in functionally important nongenic regions. This is likely to be a consequence of the low effective population size of hominids leading to a reduced effectiveness of selection. We estimate that there are one and two times as many conserved nucleotides in CNGs as in known protein-coding genes of hominids and murids, respectively. Polymorphism frequencies in CNGs indicate that purifying selection operates in these sequences. During hominid evolution, we estimate that a total of about three deleterious mutations in CNGs and protein-coding genes have been selectively eliminated per diploid genome each generation, implying that deleterious mutations are eliminated from populations non-independently and that sex is necessary for long-term population persistence.

Organization of the Caenorhabditis elegans small non-coding transcriptome: genomic features, biogenesis, and expression

Recent evidence points to considerable transcription occurring in non-protein-coding regions of eukaryote genomes. However, their lack of conservation and demonstrated function have created controversy over whether these transcripts are functional. Applying a novel cloning strategy, we have cloned 100 novel and 61 known or predicted Caenorhabditis elegans full-length ncRNAs. Studying the genomic environment and transcriptional characteristics have shown that two-thirds of all ncRNAs, including many intronic snoRNAs, are independently transcribed under the control of ncRNA-specific upstream promoter elements. Furthermore, the transcription levels of at least 60% of the ncRNAs vary with developmental stages. We identified two new classes of ncRNAs, stem-bulge RNAs (sbRNAs) and snRNA-like RNAs (snlRNAs), both featuring distinct internal motifs, secondary structures, upstream elements, and high and developmentally variable expression. Most of the novel ncRNAs are conserved in Caenorhabditis briggsae, but only one homolog was found outside the nematodes. Preliminary estimates indicate that the C. elegans transcriptome contains approximately 2700 small non-coding RNAs, potentially acting as regulatory elements in nematode development.

Wnt9b is the mutated gene involved in multifactorial nonsyndromic cleft lip with or without cleft palate in A/WySn mice, as confirmed by a genetic complementation test

BACKGROUND

Nonsyndromic cleft lip (CL) with or without cleft palate (CLP) is a common human birth defect with complex genetic etiology. One of the unidentified genes maps to chromosome 17q21. A mouse strain, A/WySn, has CLP with complex genetic etiology that models the human defect, and 1 of its causative genes, clf1, maps to a region homologous to human 17q21. Extensive studies of the candidate region pointed to a novel insertion of an IAP transposon 3' from the gene Wnt9b as the clf1 mutation. Independently a recessive knockout mutation of Wnt9b (Wnt9b-) was reported to cause a lethal syndrome that includes some CLP.

METHODS

A standard genetic test of allelism between clf1 and the Wnt9b- mutation was done. A total of 83 F1 embryos at gestation day 14 (GD 14) from Wnt9b-/+ males crossed with A/WySn females, and 79 BC1 GD 14 embryos from F1 Wnt9b-/clf1 males back-crossed to A/WySn females were observed for CL. Embryo genotypes at clf1 and Wnt9b were obtained from DNA markers. Genotypes for a second unlinked modifier locus from A/WySn, clf2, were similarly obtained.

RESULTS

The compound mutant embryos (Wnt9b-/clf1) had high frequencies of CL: 27% in the F1 and 63% in the BC1. The clf2 modifier gene was found to have 3 alleles segregating in this study and to strongly influence the penetrance of CL in the compound mutant.

CONCLUSIONS

The noncomplementation of clf1 and Wnt9b- confirms that clf1 is a mutation of the Wnt9b gene. The homologous human WNT9B gene and 3' conserved noncoding region should be examined for a role in human nonsyndromic CLP.

Genome-wide identification of human functional DNA using a neutral indel model

It has become clear that a large proportion of functional DNA in the human genome does not code for protein. Identification of this non-coding functional sequence using comparative approaches is proving difficult and has previously been thought to require deep sequencing of multiple vertebrates. Here we introduce a new model and comparative method that, instead of nucleotide substitutions, uses the evolutionary imprint of insertions and deletions (indels) to infer the past consequences of selection. The model predicts the distribution of indels under neutrality, and shows an excellent fit to human-mouse ancestral repeat data. Across the genome, many unusually long ungapped regions are detected that are unaccounted for by the neutral model, and which we predict to be highly enriched in functional DNA that has been subject to purifying selection with respect to indels. We use the model to determine the proportion under indel-purifying selection to be between 2.56% and 3.25% of human euchromatin. Since annotated protein-coding genes comprise only 1.2% of euchromatin, these results lend further weight to the proposition that more than half the functional complement of the human genome is non-protein-coding. The method is surprisingly powerful at identifying selected sequence using only two or three mammalian genomes. Applying the method to the human, mouse, and dog genomes, we identify 90 Mb of human sequence under indel-purifying selection, at a predicted 10% false-discovery rate and 75% sensitivity. As expected, most of the identified sequence represents unannotated material, while the recovered proportions of known protein-coding and microRNA genes closely match the predicted sensitivity of the method. The method's high sensitivity to functional sequence such as microRNAs suggest that as yet unannotated microRNA genes are enriched among the sequences identified. Furthermore, its independence of substitutions allowed us to identify sequence that has been subject to heterogeneous selection, that is, sequence subject to both positive selection with respect to substitutions and purifying selection with respect to indels. The ability to identify elements under heterogeneous selection enables, for the first time, the genome-wide investigation of positive selection on functional elements other than protein-coding genes.

A novel approach to identifying regulatory motifs in distantly related genomes

A two-step procedure for identifying regulatory motifs in distantly related organisms is described that combines the advantages of sequence alignment and motif detection approaches.

Although proven successful in the identification of regulatory motifs, phylogenetic footprinting methods still show some shortcomings. To assess these difficulties, most apparent when applying phylogenetic footprinting to distantly related organisms, we developed a two-step procedure that combines the advantages of sequence alignment and motif detection approaches. The results on well-studied benchmark datasets indicate that the presented method outperforms other methods when the sequences become either too long or too heterogeneous in size.

Conserved noncoding sequences are selectively constrained and not mutation cold spots

Noncoding genetic variants are likely to influence human biology and disease, but recognizing functional noncoding variants is difficult. Approximately 3% of noncoding sequence is conserved among distantly related mammals, suggesting that these evolutionarily conserved noncoding regions (CNCs) are selectively constrained and contain functional variation. However, CNCs could also merely represent regions with lower local mutation rates. Here we address this issue and show that CNCs are selectively constrained in humans by analyzing HapMap genotype data. Specifically, new (derived) alleles of SNPs within CNCs are rarer than new alleles in nonconserved regions (P = 3 x 10(-18)), indicating that evolutionary pressure has suppressed CNC-derived allele frequencies. Intronic CNCs and CNCs near genes show greater allele frequency shifts, with magnitudes comparable to those for missense variants. Thus, conserved noncoding variants are more likely to be functional. Allele frequency distributions highlight selectively constrained genomic regions that should be intensively surveyed for functionally important variation.

Trade-offs in detecting evolutionarily constrained sequence by comparative genomics

As whole-genome sequencing efforts extend beyond more traditional model organisms to include a deep diversity of species, comparative genomic analyses will be further empowered to reveal insights into the human genome and its evolution. The discovery and annotation of functional genomic elements is a necessary step toward a detailed understanding of our biology, and sequence comparisons have proven to be an integral tool for that task. This review is structured to broadly reflect the statistical challenges in discriminating these functional elements from the bulk of the genome that has evolved neutrally. Specifically, we review the comparative genomics literature in terms of specificity, sensitivity, and phylogenetic scope, as well as the trade-offs that relate these factors in standard analyses. We consider the impact of an expanding diversity of orthologous sequences on our ability to resolve functional elements. This impact is assessed through both recent comparative analyses of deep alignments and mathematical modeling.