Interpreting mammalian evolution using Fugu genome comparisons

The functional verification and analysis of Fugu promoter of cardiac gene tnni1a in zebrafish

Troponin I type 1b (Tnni1b) is thought to be a novel isoform that is expressed only in the zebrafish heart. Knocking down of tnni1b can lead to cardiac defects in zebrafish. Although both the zebrafish tnni1b and human troponin I1 (TNNI1) genes are thought to be closely associated with fatal cardiac development, the regulatory molecular mechanisms of these genes are poorly understood. Analyzing the functionally conserved sequence, especially in the noncoding regulatory region involved in gene expression, clarified these mechanisms. In this study, we isolated a 3 kb fragment upstream of Fugu tnni1a that can regulate green fluorescence protein (GFP) expression in a heart-specific manner, similar to the pattern of zebrafish homologue expression. Three evolutionarily conserved regions (ECRs) in the 5'-flanking sequence of Fugu tnni1a were identified by sequence alignment. Deletion analysis led to the identification of ECR2 as a core sequence that affects the heart-specific expression function of the Fugu tnni1a promoter. Interestingly, both the Fugu tnni1a promoter and ECR2 sequence were functionally conserved in zebrafish, although they shared no sequence similarity. Together, the findings of our study provided further evidence for the important role of tnni1a homologous in cardiac development and demonstrated that two functionally conserved sequences in the zebrafish and Fugu genomes may be ECRs, despite their lack of similarity.

Perfect and imperfect views of ultraconserved sequences

Across the human genome, there are nearly 500 'ultraconserved' elements: regions of at least 200 contiguous nucleotides that are perfectly conserved in both the mouse and rat genomes. Remarkably, the majority of these sequences are non-coding, and many can function as enhancers that activate tissue-specific gene expression during embryonic development. From their first description more than 15 years ago, their extreme conservation has both fascinated and perplexed researchers in genomics and evolutionary biology. The intrigue around ultraconserved elements only grew with the observation that they are dispensable for viability. Here, we review recent progress towards understanding the general importance and the specific functions of ultraconserved sequences in mammalian development and human disease and discuss possible explanations for their extreme conservation.

Genetic Mechanisms Underlying Cortical Evolution in Mammals

The remarkable sensory, motor, and cognitive abilities of mammals mainly depend on the neocortex. Thus, the emergence of the six-layered neocortex in reptilian ancestors of mammals constitutes a fundamental evolutionary landmark. The mammalian cortex is a columnar epithelium of densely packed cells organized in layers where neurons are generated mainly in the subventricular zone in successive waves throughout development. Newborn cells move away from their site of neurogenesis through radial or tangential migration to reach their specific destination closer to the pial surface of the same or different cortical area. Interestingly, the genetic programs underlying neocortical development diversified in different mammalian lineages. In this work, I will review several recent studies that characterized how distinct transcriptional programs relate to the development and functional organization of the neocortex across diverse mammalian lineages. In some primates such as the anthropoids, the neocortex became extremely large, especially in humans where it comprises around 80% of the brain. It has been hypothesized that the massive expansion of the cortical surface and elaboration of its connections in the human lineage, has enabled our unique cognitive capacities including abstract thinking, long-term planning, verbal language and elaborated tool making capabilities. I will also analyze the lineage-specific genetic changes that could have led to the modification of key neurodevelopmental events, including regulation of cell number, neuronal migration, and differentiation into specific phenotypes, in order to shed light on the evolutionary mechanisms underlying the diversity of mammalian brains including the human brain.

Rare and de novo duplications containing SHOX in clubfoot

Introduction

Congenital clubfoot is a common birth defect that affects at least 0.1% of all births. Nearly 25% cases are familial and the remaining are sporadic in inheritance. Copy number variants (CNVs) involving transcriptional regulators of limb development, including PITX1 and TBX4, have previously been shown to cause familial clubfoot, but much of the heritability remains unexplained.

Methods

Exome sequence data from 816 unrelated clubfoot cases and 2645 in-house controls were analysed using coverage data to identify rare CNVs. The precise size and location of duplications were then determined using high-density Affymetrix Cytoscan chromosomal microarray (CMA). Segregation in families and de novo status were determined using qantitative PCR.

Results

Chromosome Xp22.33 duplications involving SHOX were identified in 1.1% of cases (9/816) compared with 0.07% of in-house controls (2/2645) (p=7.98×10⁻⁵, OR=14.57) and 0.27% (38/13592) of Atherosclerosis Risk in Communities/the Wellcome Trust Case Control Consortium 2 controls (p=0.001, OR=3.97). CMA validation confirmed an overlapping 180.28 kb duplicated region that included SHOX exons as well as downstream non-coding regions. In four of six sporadic cases where DNA was available for unaffected parents, the duplication was de novo. The probability of four de novo mutations in SHOX by chance in a cohort of 450 sporadic clubfoot cases is 5.4×10^–10.

Conclusions

Microduplications of the pseudoautosomal chromosome Xp22.33 region (PAR1) containing SHOX and downstream enhancer elements occur in ~1% of patients with clubfoot. SHOX and regulatory regions have previously been implicated in skeletal dysplasia as well as idiopathic short stature, but have not yet been reported in clubfoot. SHOX duplications likely contribute to clubfoot pathogenesis by altering early limb development.

The Conserved Sonic Hedgehog Limb Enhancer Consists of Discrete Functional Elements that Regulate Precise Spatial Expression

Expression of sonic hedgehog (Shh) in the limb bud is regulated by an enhancer called the zone of polarizing activity regulatory sequence (ZRS), which, in evolution, belongs to an ancient group of highly conserved cis regulators found in all classes of vertebrates. Here, we examined the endogenous ZRS in mice, using genome editing to establish the relationship between enhancer composition and embryonic phenotype. We show that enhancer activity is a consolidation of distinct activity domains. Spatial restriction of Shh expression is mediated by a discrete repressor module, whereas levels of gene expression are controlled by large overlapping domains containing varying numbers of HOXD binding sites. The number of HOXD binding sites regulates expression levels incrementally. Substantial portions of conserved sequence are dispensable, indicating the presence of sequence redundancy. We propose a collective model for enhancer activity in which function is an integration of discrete expression activities and redundant components that drive robust expression.

•

The ancient vertebrate enhancer, the ZRS, shows sequence plasticity

•

Discrete regulatory activities are assigned to specific sites in the enhancer

•

The number of HOXD binding sites determines the level of Shh expression

•

Robust expression is a collective of regulatory and redundant information

Identification of Candidate Genes for Reactivity in Guzerat (Bos indicus) Cattle: A Genome-Wide Association Study

Temperament is fundamental to animal production due to its direct influence on the animal-herdsman relationship. When compared to calm animals, the aggressive, anxious or fearful ones exhibit less weight gain, lower reproductive efficiency, decreased milk production and higher herd maintenance costs, all of which contribute to reduced profits. However, temperament is a trait that is complex and difficult to assess. Recently, a new quantitative system, REATEST^®, for assessing reactivity, a phenotype of temperament, was developed. Herein, we describe the results of a Genome-wide association study for reactivity, assessed using REATEST^® with a sample of 754 females from five dual-purpose (milk and meat production) Guzerat (Bos indicus) herds. Genotyping was performed using a 50k SNP chip and a two-step mixed model approach (Grammar-Gamma) with a one-by-one marker regression was used to identify QTLs. QTLs for reactivity were identified on chromosomes BTA1, BTA5, BTA14, and BTA25. Five intronic and two intergenic markers were significantly associated with reactivity. POU1F1, DRD3, VWA3A, ZBTB20, EPHA6, SNRPF and NTN4 were identified as candidate genes. Previous QTL reports for temperament traits, covering areas surrounding the SNPs/genes identified here, further corroborate these associations. The seven genes identified in the present study explain 20.5% of reactivity variance and give a better understanding of temperament biology.

An Ultraconserved Brain-Specific Enhancer Within ADGRL3 (LPHN3) Underpins Attention-Deficit/Hyperactivity Disorder Susceptibility

BACKGROUND

Genetic factors predispose individuals to attention-deficit/hyperactivity disorder (ADHD). Previous studies have reported linkage and association to ADHD of gene variants within ADGRL3. In this study, we functionally analyzed noncoding variants in this gene as likely pathological contributors.

METHODS

In silico, in vitro, and in vivo approaches were used to identify and characterize evolutionary conserved elements within the ADGRL3 linkage region (~207 Kb). Family-based genetic analyses of 838 individuals (372 affected and 466 unaffected patients) identified ADHD-associated single nucleotide polymorphisms harbored in some of these conserved elements. Luciferase assays and zebrafish green fluorescent protein transgenesis tested conserved elements for transcriptional enhancer activity. Electromobility shift assays were used to verify transcription factor-binding disruption by ADHD risk alleles.

RESULTS

An ultraconserved element was discovered (evolutionary conserved region 47) that functions as a transcriptional enhancer. A three-variant ADHD risk haplotype in evolutionary conserved region 47, formed by rs17226398, rs56038622, and rs2271338, reduced enhancer activity by 40% in neuroblastoma and astrocytoma cells (p_Bonferroni < .0001). This enhancer also drove green fluorescent protein expression in the zebrafish brain in a tissue-specific manner, sharing aspects of endogenous ADGRL3 expression. The rs2271338 risk allele disrupts binding of YY1 transcription factor, an important factor in the development and function of the central nervous system. Expression quantitative trait loci analysis of postmortem human brain tissues revealed an association between rs2271338 and reduced ADGRL3 expression in the thalamus.

CONCLUSIONS

These results uncover the first functional evidence of common noncoding variants with potential implications for the pathology of ADHD.

The disruption of a novel limb cis-regulatory element of SHH is associated with autosomal dominant preaxial polydactyly-hypertrichosis

The expression gradient of the morphogen Sonic Hedgehog (SHH) is crucial in establishing the number and the identity of the digits during anteroposterior patterning of the limb. Its anterior ectopic expression is responsible for preaxial polydactyly (PPD). Most of these malformations are due to the gain-of-function of the Zone of Polarizing Activity Regulatory Sequence, the only limb-specific enhancer of SHH known to date. We report a family affected with a novel condition associating PPD and hypertrichosis of the upper back, following an autosomal dominant mode of inheritance. This phenotype is consistent with deregulation of SHH expression during limb and follicle development. In affected members, we identified a 2 kb deletion located ~240 kb upstream from the SHH promoter. The deleted sequence is capable of repressing the transcriptional activity of the SHH promoter in vitro, consistent with a silencer activity. We hypothesize that the deletion of this silencer could be responsible for SHH deregulation during development, leading to a PPD-hypertrichosis phenotype.

Identification and computational analysis of gene regulatory elements

Over the last two decades, advances in experimental and computational technologies have greatly facilitated genomic research. Next-generation sequencing technologies have made de novo sequencing of large genomes affordable, and powerful computational approaches have enabled accurate annotations of genomic DNA sequences. Charting functional regions in genomes must account for not only the coding sequences, but also noncoding RNAs, repetitive elements, chromatin states, epigenetic modifications, and gene regulatory elements. A mix of comparative genomics, high-throughput biological experiments, and machine learning approaches has played a major role in this truly global effort. Here we describe some of these approaches and provide an account of our current understanding of the complex landscape of the human genome. We also present overviews of different publicly available, large-scale experimental data sets and computational tools, which we hope will prove beneficial for researchers working with large and complex genomes.

Characterisation of multiple regulatory domains spanning the major transcriptional start site of the FUS gene, a candidate gene for motor neurone disease

Fused-In-Sarcoma (FUS) is a candidate gene for neurological disorders including motor neurone disease and Parkinson׳s disease in addition to various types of cancer. Recently it has been reported that over expression of FUS causes motor neurone disease in mouse models hence mutations leading to changes in gene expression may contribute to the development of neurodegenerative disease. Genome evolutionary conservation was used to predict important cis-acting DNA regulators of the FUS gene promoter that direct transcription. The putative regulators identified were analysed in reporter gene assays in cells and in chick embryos. Our analysis indicated in addition to regulatory domains 5' of the transcriptional start site an important regulatory domain resides in intron 1 of the gene itself. This intronic domain functioned both in cell lines and in vivo in the neural tube of the chick embryo including developing motor neurones. Our data suggest the interaction of multiple domains including intronic domains are involved in expression of FUS. A better understanding of the regulation of expression of FUS may give insight into how its stimulus inducible expression may be associated with neurological disorders.

Sp1 is necessary for gene activation of Adamts17 by estrogen

Adamts17 is a member of a family of secreted metalloproteinases. In this report, we show that knockdown of Adamts17 expression induces apoptosis and inhibits breast cancer cell growth. Adamts17 expression can rapidly be induced by estrogens. siRNA knockdown of Sp1 or Myc demonstrated that Sp1 is required to induce Adamts17 gene expression in response to estrogen. Moreover, reporter assays showed that the proximal promoter and the upstream sequences were not capable of conferring estrogen responsiveness, suggesting that Sp1 elements may be located in the downstream intronic region. We further demonstrated that Sp1 and Myc binding in the proximal promoter region contributed to the Adamts17 basal expression. Furthermore, histone deacetylase (HDAC) and methylase inhibitors also induced Adamts17 expression, indicating that epigenetic alterations, such as aberrant HDAC and/or methylation are associated with dysregulated Adamts17 expression. By meta-analysis using Oncomine microarray data, we found that higher Adamts17 expression is found in several human cancer cell subtypes, especially in breast ductal carcinoma. Moreover, we found that there is an inverse correlation between higher Adamts17 expression and patients' survival. Our study suggests that Adamts17 may support breast cancer cell growth and survival.

Intrinsic properties of Tcf1 and Tcf4 splice variants determine cell-type-specific Wnt/β-catenin target gene expression

T-cell factor (Tcf)/lymphoid-enhancer factor (Lef) proteins are a structurally diverse family of deoxyribonucleic acid-binding proteins that have essential nuclear functions in Wnt/β-catenin signalling. Expression of Wnt/β-catenin target genes is highly dependent on context, but the precise role of Tcf/Lef family members in the generation and maintenance of cell-type-specific Wnt/β-catenin responses is unknown. Herein, we show that induction of a subset of Wnt/β-catenin targets in embryonic stem cells depends on Tcf1 and Tcf4, whereas other co-expressed Tcf/Lef family members cannot induce these targets. The Tcf1/Tcf4-dependent gene responses to Wnt are primarily if not exclusively mediated by C-clamp-containing Tcf1E and Tcf4E splice variants. A combined knockdown of Tcf1/Tcf4 abrogates Wnt-inducible transcription but does not affect the active chromatin conformation of their targets. Thus, the transcriptionally poised state of Wnt/β-catenin targets is maintained independent of Tcf/Lef proteins. Conversely, ectopically overexpressed Tcf1E cannot invade silent chromatin and fails to initiate expression of inactive Wnt/β-catenin targets even if repressive chromatin modifications are abolished. The observed non-redundant functions of Tcf1/Tcf4 isoforms in acute transcriptional activation demonstrated that the cell-type-specific complement of Tcf/Lef proteins is a critical determinant of context-dependent Wnt/β-catenin responses. Moreover, the apparent inability to cope with chromatin uncovers an intrinsic property of Tcf/Lef proteins that prevents false ectopic induction and ensures spatiotemporal stability of Wnt/β-catenin target gene expression.

cJun integrates calcium activity and tlx3 expression to regulate neurotransmitter specification

Neuronal differentiation is accomplished through cascades of intrinsic genetic factors initiated in neuronal progenitors by external gradients of morphogens. Activity was thought to be important only late in development, but recent evidence indicates that activity also regulates early neuronal differentiation. Activity in post-mitotic neurons prior to synapse formation can regulate phenotypic specification, including neurotransmitter choice, but the mechanisms are not clear. Here we identify a mechanism that links endogenous calcium spike activity with an intrinsic genetic pathway to specify neurotransmitter choice in neurons in the dorsal embryonic spinal cord of Xenopus tropicalis. Early activity modulates transcription of the GABAergic/glutamatergic selection gene tlx3 and requires a variant cAMP response element (CRE) in its promoter. The cJun transcription factor binds to this CRE site, modulates transcription, and regulates neurotransmitter phenotype through its transactivation domain. Calcium signals through cJun N-terminal phosphorylation, thus integrating activity-dependent and intrinsic neurotransmitter specification. This mechanism provides a basis for early activity to regulate genetic pathways at critical decision points, switching the phenotype of developing neurons.

Homotypic clusters of transcription factor binding sites are a key component of human promoters and enhancers

Clustering of multiple transcription factor binding sites (TFBSs) for the same transcription factor (TF) is a common feature of cis-regulatory modules in invertebrate animals, but the occurrence of such homotypic clusters of TFBSs (HCTs) in the human genome has remained largely unknown. To explore whether HCTs are also common in human and other vertebrates, we used known binding motifs for vertebrate TFs and a hidden Markov model-based approach to detect HCTs in the human, mouse, chicken, and fugu genomes, and examined their association with cis-regulatory modules. We found that evolutionarily conserved HCTs occupy nearly 2% of the human genome, with experimental evidence for individual TFs supporting their binding to predicted HCTs. More than half of the promoters of human genes contain HCTs, with a distribution around the transcription start site in agreement with the experimental data from the ENCODE project. In addition, almost half of the 487 experimentally validated developmental enhancers contain them as well--a number more than 25-fold larger than expected by chance. We also found evidence of negative selection acting on TFBSs within HCTs, as the conservation of TFBSs is stronger than the conservation of sequences separating them. The important role of HCTs as components of developmental enhancers is additionally supported by a strong correlation between HCTs and the binding of the enhancer-associated coactivator protein Ep300 (also known as p300). Experimental validation of HCT-containing elements in both zebrafish and mouse suggest that HCTs could be used to predict both the presence of enhancers and their tissue specificity, and are thus a feature that can be effectively used in deciphering the gene regulatory code. In conclusion, our results indicate that HCTs are a pervasive feature of human cis-regulatory modules and suggest that they play an important role in gene regulation in the human and other vertebrate genomes.

Human variation in short regions predisposed to deep evolutionary conservation

The landscape of the human genome consists of millions of short islands of conservation that are 100% conserved across multiple vertebrate genomes (termed “bricks”), the majority of which are located in noncoding regions. Several hundred thousand bricks are deeply conserved reaching the genomes of amphibians and fish. Deep phylogenetic conservation of noncoding DNA has been reported to be strongly associated with the presence of gene regulatory elements, introducing bricks as a proxy to the functional noncoding landscape of the human genome. Here, we report a significant overrepresentation of bricks in the promoters of transcription factors and developmental genes, where the high level of phylogenetic conservation correlates with an increase in brick overrepresentation. We also found that the presence of a brick dictates a predisposition to evolutionary constraint, with only 0.7% of the amniota brick central nucleotides being diverged within the primate lineage—an 11-fold reduction in the divergence rate compared with random expectation. Human single-nucleotide polymorphism (SNP) data explains only 3% of primate-specific variation in amniota bricks, thus arguing for a widespread fixation of brick mutations within the primate lineage and prior to human radiation. This variation, in turn, might have been utilized as a driving force for primate- and hominoid-specific adaptation. We also discovered a pronounced deviation from the evolutionary predisposition in the human lineage, with over 20-fold increase in the substitution rate at brick SNP sites over expected values. In addition, contrary to typical brick mutations, brick variation commonly encountered in the human population displays limited, if any, signatures of negative selection as measured by the minor allele frequency and population differentiation (F-statistical measure) measures. These observations argue for the plasticity of gene regulatory mechanisms in vertebrates—with evidence of strong purifying selection acting on the gene regulatory landscape of the human genome, where widespread advantageous mutations in putative regulatory elements are likely utilized in functional diversification and adaptation of species.

Functional and comparative genomics analyses of pmp22 in medaka fish

Background

Pmp22, a member of the junction protein family Claudin/EMP/PMP22, plays an important role in myelin formation. Increase of pmp22 transcription causes peripheral neuropathy, Charcot-Marie-Tooth disease type1A (CMT1A). The pathophysiological phenotype of CMT1A is aberrant axonal myelination which induces a reduction in nerve conduction velocity (NCV). Several CMT1A model rodents have been established by overexpressing pmp22. Thus, it is thought that pmp22 expression must be tightly regulated for correct myelin formation in mammals. Interestingly, the myelin sheath is also present in other jawed vertebrates. The purpose of this study is to analyze the evolutionary conservation of the association between pmp22 transcription level and vertebrate myelin formation, and to find the conserved non-coding sequences for pmp22 regulation by comparative genomics analyses between jawed fishes and mammals.

Results

A transgenic pmp22 over-expression medaka fish line was established. The transgenic fish had approximately one fifth the peripheral NCV values of controls, and aberrant myelination of transgenic fish in the peripheral nerve system (PNS) was observed. We successfully confirmed that medaka fish pmp22 has the same exon-intron structure as mammals, and identified some known conserved regulatory motifs. Furthermore, we found novel conserved sequences in the first intron and 3'UTR.

Conclusion

Medaka fish undergo abnormalities in the PNS when pmp22 transcription increases. This result indicates that an adequate pmp22 transcription level is necessary for correct myelination of jawed vertebrates. Comparison of pmp22 orthologs between distantly related species identifies evolutionary conserved sequences that contribute to precise regulation of pmp22 expression.

Variable locus length in the human genome leads to ascertainment bias in functional inference for non-coding elements

MOTIVATION

Several functional gene annotation databases have been developed in the recent years, and are widely used to infer the biological function of gene sets, by scrutinizing the attributes that appear over- and underrepresented. However, this strategy is not directly applicable to the study of non-coding DNA, as the non-coding sequence span varies greatly among different gene loci in the human genome and longer loci have a higher likelihood of being selected purely by chance. Therefore, conclusions involving the function of non-coding elements that are drawn based on the annotation of neighboring genes are often biased. We assessed the systematic bias in several particular Gene Ontology (GO) categories using the standard hypergeometric test, by randomly sampling non-coding elements from the human genome and inferring their function based on the functional annotation of the closest genes. While no category is expected to occur significantly over- or underrepresented for a random selection of elements, categories such as 'cell adhesion', 'nervous system development' and 'transcription factor activities' appeared to be systematically overrepresented, while others such as 'olfactory receptor activity'-underrepresented.

RESULTS

Our results suggest that functional inference for non-coding elements using gene annotation databases requires a special correction. We introduce a set of correction coefficients for the probabilities of the GO categories that accounts for the variability in the length of the non-coding DNA across different loci and effectively eliminates the ascertainment bias from the functional characterization of non-coding elements. Our approach can be easily generalized to any other gene annotation database.

Evolution of conserved secondary structures and their function in transcriptional regulation networks

Background

Many conserved secondary structures have been identified within conserved elements in the human genome, but only a small fraction of them are known to be functional RNAs. The evolutionary variations of these conserved secondary structures in human populations and their biological functions have not been fully studied.

Results

We searched for polymorphisms within conserved secondary structures and identified a number of SNPs within these elements even though they are highly conserved among species. The density of SNPs in conserved secondary structures is about 65% of that of their flanking, non-conserved, sequences. Classification of sites as stems or as loops/bulges revealed that the density of SNPs in stems is about 62% of that found in loops/bulges. Analysis of derived allele frequency data indicates that sites in stems are under stronger evolutionary constraint than sites in loops/bulges. Intergenic conserved secondary structures tend to associate with transcription factor-encoding genes with genetic distance being the measure of regulator-gene associations. A substantial fraction of intergenic conserved secondary structures overlap characterized binding sites for multiple transcription factors.

Conclusion

Strong purifying selection implies that secondary structures are probably important carriers of biological functions for conserved sequences. The overlap between intergenic conserved secondary structures and transcription factor binding sites further suggests that intergenic conserved secondary structures have essential roles in directing gene expression in transcriptional regulation networks.

NcDNAlign: plausible multiple alignments of non-protein-coding genomic sequences

Genome-wide multiple sequence alignments (MSAs) are a necessary prerequisite for an increasingly diverse collection of comparative genomic approaches. Here we present a versatile method that generates high-quality MSAs for non-protein-coding sequences. The NcDNAlign pipeline combines pairwise BLAST alignments to create initial MSAs, which are then locally improved and trimmed. The program is optimized for speed and hence is particulary well-suited to pilot studies. We demonstrate the practical use of NcDNAlign in three case studies: the search for ncRNAs in gammaproteobacteria and the analysis of conserved noncoding DNA in nematodes and teleost fish, in the latter case focusing on the fate of duplicated ultra-conserved regions. Compared to the currently widely used genome-wide alignment program TBA, our program results in a 20- to 30-fold reduction of CPU time necessary to generate gammaproteobacterial alignments. A showcase application of bacterial ncRNA prediction based on alignments of both algorithms results in similar sensitivity, false discovery rates, and up to 100 putatively novel ncRNA structures. Similar findings hold for our application of NcDNAlign to the identification of ultra-conserved regions in nematodes and teleosts. Both approaches yield conserved sequences of unknown function, result in novel evolutionary insights into conservation patterns among these genomes, and manifest the benefits of an efficient and reliable genome-wide alignment package. The software is available under the GNU Public License at http://www.bioinf.uni-leipzig.de/Software/NcDNAlign/.

Genome of the marsupial Monodelphis domestica reveals innovation in non-coding sequences

We report a high-quality draft of the genome sequence of the grey, short-tailed opossum (Monodelphis domestica). As the first metatherian ('marsupial') species to be sequenced, the opossum provides a unique perspective on the organization and evolution of mammalian genomes. Distinctive features of the opossum chromosomes provide support for recent theories about genome evolution and function, including a strong influence of biased gene conversion on nucleotide sequence composition, and a relationship between chromosomal characteristics and X chromosome inactivation. Comparison of opossum and eutherian genomes also reveals a sharp difference in evolutionary innovation between protein-coding and non-coding functional elements. True innovation in protein-coding genes seems to be relatively rare, with lineage-specific differences being largely due to diversification and rapid turnover in gene families involved in environmental interactions. In contrast, about 20% of eutherian conserved non-coding elements (CNEs) are recent inventions that postdate the divergence of Eutheria and Metatheria. A substantial proportion of these eutherian-specific CNEs arose from sequence inserted by transposable elements, pointing to transposons as a major creative force in the evolution of mammalian gene regulation.

Predicting tissue-specific enhancers in the human genome

Determining how transcriptional regulatory signals are encoded in vertebrate genomes is essential for understanding the origins of multicellular complexity; yet the genetic code of vertebrate gene regulation remains poorly understood. In an attempt to elucidate this code, we synergistically combined genome-wide gene-expression profiling, vertebrate genome comparisons, and transcription factor binding-site analysis to define sequence signatures characteristic of candidate tissue-specific enhancers in the human genome. We applied this strategy to microarray-based gene expression profiles from 79 human tissues and identified 7187 candidate enhancers that defined their flanking gene expression, the majority of which were located outside of known promoters. We cross-validated this method for its ability to de novo predict tissue-specific gene expression and confirmed its reliability in 57 of the 79 available human tissues, with an average precision in enhancer recognition ranging from 32% to 63% and a sensitivity of 47%. We used the sequence signatures identified by this approach to successfully assign tissue-specific predictions to approximately 328,000 human-mouse conserved noncoding elements in the human genome. By overlapping these genome-wide predictions with a data set of enhancers validated in vivo, in transgenic mice, we were able to confirm our results with a 28% sensitivity and 50% precision. These results indicate the power of combining complementary genomic data sets as an initial computational foray into a global view of tissue-specific gene regulation in vertebrates.

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.

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.

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.

Human-zebrafish non-coding conserved elements act in vivo to regulate transcription

Whole genome comparisons of distantly related species effectively predict biologically important sequences--core genes and cis-acting regulatory elements (REs)--but require experimentation to verify biological activity. To examine the efficacy of comparative genomics in identification of active REs from anonymous, non-coding (NC) sequences, we generated a novel alignment of the human and draft zebrafish genomes, and contrasted this set to existing human and fugu datasets. We tested the transcriptional regulatory potential of candidate sequences using two in vivo assays. Strict selection of non-genic elements which are deeply conserved in vertebrate evolution identifies 1744 core vertebrate REs in human and two fish genomes. We tested 16 elements in vivo for cis-acting gene regulatory properties using zebrafish transient transgenesis and found that 10 (63%) strongly modulate tissue-specific expression of a green fluorescent protein reporter vector. We also report a novel quantitative enhancer assay with potential for increased throughput based on normalized luciferase activity in vivo. This complementary system identified 11 (69%; including 9 of 10 GFP-confirmed elements) with cis-acting function. Together, these data support the utility of comparative genomics of distantly related vertebrate species to identify REs and provide a scaleable, in vivo quantitative assay to define functional activity of candidate REs.

Prediction of cis-regulatory elements using binding site matrices--the successes, the failures and the reasons for both

Protein-DNA interactions control many aspects of animal development and cellular responses to the environment. Although profiling of individual transcription factor binding sites is not a reliable guide for predicting the position of cis-regulatory elements in large genomes, modelling the evolution and the organization of regulatory elements has provided enough information to make some successful predictions. For vertebrate genomes, the field is limited by the lack of sufficient experimental data upon which to build reliable models. Nonetheless, a combination of experimental, computational and comparative data is likely to reveal aspects of complex regulatory networks in vertebrates, just as it has already done for simple eukaryotic genomes.

Identifying synonymous regulatory elements in vertebrate genomes

Synonymous gene regulation, defined by regulatory elements driving shared temporal and/or spatial aspects of gene expression, is most probably predicated on genomic elements that contain similar modules of certain transcription factor binding sites (TFBS). We have developed a method to scan vertebrate genomes for evolutionary conserved modules of TFBS in a predefined configuration, and created a tool, named SynoR that identifies synonymous regulatory elements (SREs) in vertebrate genomes. SynoR performs de novo identification of SREs utilizing known patterns of TFBS in active regulatory elements (REs) as seeds for genome scans. Layers of multiple-species conservation allow the use of differential phylogenetic sequence conservation filters in search of SREs and the results are displayed such as to provide an extensive annotation of the genes containing the detected REs. Gene Ontology categories are utilized to further functionally classify the identified genes, and integrated GNF Expression Atlas 2 data allow the cataloging of tissue-specificities of the predicted SREs. SynoR is publicly available at .

Dcode.org anthology of comparative genomic tools

Comparative genomics provides the means to demarcate functional regions in anonymous DNA sequences. The successful application of this method to identifying novel genes is currently shifting to deciphering the non-coding encryption of gene regulation across genomes. To facilitate the practical application of comparative sequence analysis to genetics and genomics, we have developed several analytical and visualization tools for the analysis of arbitrary sequences and whole genomes. These tools include two alignment tools, zPicture and Mulan; a phylogenetic shadowing tool, eShadow for identifying lineage- and species-specific functional elements; two evolutionary conserved transcription factor analysis tools, rVista and multiTF; a tool for extracting cis-regulatory modules governing the expression of co-regulated genes, Creme 2.0; and a dynamic portal to multiple vertebrate and invertebrate genome alignments, the ECR Browser. Here, we briefly describe each one of these tools and provide specific examples on their practical applications. All the tools are publicly available at the http://www.dcode.org/ website.

Evolution and functional classification of vertebrate gene deserts

Large tracts of the human genome, known as gene deserts, are devoid of protein-coding genes. Dichotomy in their level of conservation with chicken separates these regions into two distinct categories, stable and variable. The separation is not caused by differences in rates of neutral evolution but instead appears to be related to different biological functions of stable and variable gene deserts in the human genome. Gene Ontology categories of the adjacent genes are strongly biased toward transcriptional regulation and development for the stable gene deserts, and toward distinctively different functions for the variable gene deserts. Stable gene deserts resist chromosomal rearrangements and appear to harbor multiple distant regulatory elements physically linked to their neighboring genes, with the linearity of conservation invariant throughout vertebrate evolution.

Megabase deletions of gene deserts result in viable mice

The functional importance of the roughly 98% of mammalian genomes not corresponding to protein coding sequences remains largely undetermined. Here we show that some large-scale deletions of the non-coding DNA referred to as gene deserts can be well tolerated by an organism. We deleted two large non-coding intervals, 1,511 kilobases and 845 kilobases in length, from the mouse genome. Viable mice homozygous for the deletions were generated and were indistinguishable from wild-type littermates with regard to morphology, reproductive fitness, growth, longevity and a variety of parameters assaying general homeostasis. Further detailed analysis of the expression of multiple genes bracketing the deletions revealed only minor expression differences in homozygous deletion and wild-type mice. Together, the two deleted segments harbour 1,243 non-coding sequences conserved between humans and rodents (more than 100 base pairs, 70% identity). Some of the deleted sequences might encode for functions unidentified in our screen; nonetheless, these studies further support the existence of potentially 'disposable DNA' in the genomes of mammals.