Long read sequencing characterises a novel structural variant, revealing underactive AKR1C1 with overactive AKR1C2 as a possible cause of severe chronic fatigue

BACKGROUND

Causative genetic variants cannot yet be found for many disorders with a clear heritable component, including chronic fatigue disorders like myalgic encephalomyelitis/chronic fatigue syndrome (ME/CFS). These conditions may involve genes in difficult-to-align genomic regions that are refractory to short read approaches. Structural variants in these regions can be particularly hard to detect or define with short reads, yet may account for a significant number of cases. Long read sequencing can overcome these difficulties but so far little data is available regarding the specific analytical challenges inherent in such regions, which need to be taken into account to ensure that variants are correctly identified. Research into chronic fatigue disorders faces the additional challenge that the heterogeneous patient populations likely encompass multiple aetiologies with overlapping symptoms, rather than a single disease entity, such that each individual abnormality may lack statistical significance within a larger sample. Better delineation of patient subgroups is needed to target research and treatment.

METHODS

We use nanopore sequencing in a case of unexplained severe fatigue to identify and fully characterise a large inversion in a highly homologous region spanning the AKR1C gene locus, which was indicated but could not be resolved by short-read sequencing. We then use GC-MS/MS serum steroid analysis to investigate the functional consequences.

RESULTS

Several commonly used bioinformatics tools are confounded by the homology but a combined approach including visual inspection allows the variant to be accurately resolved. The DNA inversion appears to increase the expression of AKR1C2 while limiting AKR1C1 activity, resulting in a relative increase of inhibitory GABAergic neurosteroids and impaired progesterone metabolism which could suppress neuronal activity and interfere with cellular function in a wide range of tissues.

CONCLUSIONS

This study provides an example of how long read sequencing can improve diagnostic yield in research and clinical care, and highlights some of the analytical challenges presented by regions containing tandem arrays of genes. It also proposes a novel gene associated with a novel disease aetiology that may be an underlying cause of complex chronic fatigue. It reveals biomarkers that could now be assessed in a larger cohort, potentially identifying a subset of patients who might respond to treatments suggested by the aetiology.

Nuclear-embedded mitochondrial DNA sequences in 66,083 human genomes

DNA transfer from cytoplasmic organelles to the cell nucleus is a legacy of the endosymbiotic event-the majority of nuclear-mitochondrial segments (NUMTs) are thought to be ancient, preceding human speciation^1-3. Here we analyse whole-genome sequences from 66,083 people-including 12,509 people with cancer-and demonstrate the ongoing transfer of mitochondrial DNA into the nucleus, contributing to a complex NUMT landscape. More than 99% of individuals had at least one of 1,637 different NUMTs, with 1 in 8 individuals having an ultra-rare NUMT that is present in less than 0.1% of the population. More than 90% of the extant NUMTs that we evaluated inserted into the nuclear genome after humans diverged from apes. Once embedded, the sequences were no longer under the evolutionary constraint seen within the mitochondrion, and NUMT-specific mutations had a different mutational signature to mitochondrial DNA. De novo NUMTs were observed in the germline once in every 10⁴ births and once in every 10³ cancers. NUMTs preferentially involved non-coding mitochondrial DNA, linking transcription and replication to their origin, with nuclear insertion involving multiple mechanisms including double-strand break repair associated with PR domain zinc-finger protein 9 (PRDM9) binding. The frequency of tumour-specific NUMTs differed between cancers, including a probably causal insertion in a myxoid liposarcoma. We found evidence of selection against NUMTs on the basis of size and genomic location, shaping a highly heterogenous and dynamic human NUMT landscape.

Whole-genome sequencing reveals host factors underlying critical COVID-19

Critical COVID-19 is caused by immune-mediated inflammatory lung injury. Host genetic variation influences the development of illness requiring critical care1 or hospitalization2-4 after infection with SARS-CoV-2. The GenOMICC (Genetics of Mortality in Critical Care) study enables the comparison of genomes from individuals who are critically ill with those of population controls to find underlying disease mechanisms. Here we use whole-genome sequencing in 7,491 critically ill individuals compared with 48,400 controls to discover and replicate 23 independent variants that significantly predispose to critical COVID-19. We identify 16 new independent associations, including variants within genes that are involved in interferon signalling (IL10RB and PLSCR1), leucocyte differentiation (BCL11A) and blood-type antigen secretor status (FUT2). Using transcriptome-wide association and colocalization to infer the effect of gene expression on disease severity, we find evidence that implicates multiple genes-including reduced expression of a membrane flippase (ATP11A), and increased expression of a mucin (MUC1)-in critical disease. Mendelian randomization provides evidence in support of causal roles for myeloid cell adhesion molecules (SELE, ICAM5 and CD209) and the coagulation factor F8, all of which are potentially druggable targets. Our results are broadly consistent with a multi-component model of COVID-19 pathophysiology, in which at least two distinct mechanisms can predispose to life-threatening disease: failure to control viral replication; or an enhanced tendency towards pulmonary inflammation and intravascular coagulation. We show that comparison between cases of critical illness and population controls is highly efficient for the detection of therapeutically relevant mechanisms of disease.

Inherited duplications of PPP2R3B predispose to nevi and melanoma via a C21orf91-driven proliferative phenotype

PURPOSE

Much of the heredity of melanoma remains unexplained. We sought predisposing germline copy-number variants using a rare disease approach.

METHODS

Whole-genome copy-number findings in patients with melanoma predisposition syndrome congenital melanocytic nevus were extrapolated to a sporadic melanoma cohort. Functional effects of duplications in PPP2R3B were investigated using immunohistochemistry, transcriptomics, and stable inducible cellular models, themselves characterized using RNAseq, quantitative real-time polymerase chain reaction (qRT-PCR), reverse phase protein arrays, immunoblotting, RNA interference, immunocytochemistry, proliferation, and migration assays.

RESULTS

We identify here a previously unreported genetic susceptibility to melanoma and melanocytic nevi, familial duplications of gene PPP2R3B. This encodes PR70, a regulatory unit of critical phosphatase PP2A. Duplications increase expression of PR70 in human nevus, and increased expression in melanoma tissue correlates with survival via a nonimmunological mechanism. PPP2R3B overexpression induces pigment cell switching toward proliferation and away from migration. Importantly, this is independent of the known microphthalmia-associated transcription factor (MITF)-controlled switch, instead driven by C21orf91. Finally, C21orf91 is demonstrated to be downstream of MITF as well as PR70.

CONCLUSION

This work confirms the power of a rare disease approach, identifying a previously unreported copy-number change predisposing to melanocytic neoplasia, and discovers C21orf91 as a potentially targetable hub in the control of phenotype switching.

An ancestral 10-bp repeat expansion in VWA1 causes recessive hereditary motor neuropathy

The extracellular matrix comprises a network of macromolecules such as collagens, proteoglycans and glycoproteins. VWA1 (von Willebrand factor A domain containing 1) encodes a component of the extracellular matrix that interacts with perlecan/collagen VI, appears to be involved in stabilizing extracellular matrix structures, and demonstrates high expression levels in tibial nerve. Vwa1-deficient mice manifest with abnormal peripheral nerve structure/function; however, VWA1 variants have not previously been associated with human disease. By interrogating the genome sequences of 74 180 individuals from the 100K Genomes Project in combination with international gene-matching efforts and targeted sequencing, we identified 17 individuals from 15 families with an autosomal-recessive, non-length dependent, hereditary motor neuropathy and rare biallelic variants in VWA1. A single disease-associated allele p.(G25Rfs*74), a 10-bp repeat expansion, was observed in 14/15 families and was homozygous in 10/15. Given an allele frequency in European populations approaching 1/1000, the seven unrelated homozygote individuals ascertained from the 100K Genomes Project represents a substantial enrichment above expected. Haplotype analysis identified a shared 220 kb region suggesting that this founder mutation arose >7000 years ago. A wide age-range of patients (6-83 years) helped delineate the clinical phenotype over time. The commonest disease presentation in the cohort was an early-onset (mean 2.0 ± 1.4 years) non-length-dependent axonal hereditary motor neuropathy, confirmed on electrophysiology, which will have to be differentiated from other predominantly or pure motor neuropathies and neuronopathies. Because of slow disease progression, ambulation was largely preserved. Neurophysiology, muscle histopathology, and muscle MRI findings typically revealed clear neurogenic changes with single isolated cases displaying additional myopathic process. We speculate that a few findings of myopathic changes might be secondary to chronic denervation rather than indicating an additional myopathic disease process. Duplex reverse transcription polymerase chain reaction and immunoblotting using patient fibroblasts revealed that the founder allele results in partial nonsense mediated decay and an absence of detectable protein. CRISPR and morpholino vwa1 modelling in zebrafish demonstrated reductions in motor neuron axonal growth, synaptic formation in the skeletal muscles and locomotive behaviour. In summary, we estimate that biallelic variants in VWA1 may be responsible for up to 1% of unexplained hereditary motor neuropathy cases in Europeans. The detailed clinical characterization provided here will facilitate targeted testing on suitable patient cohorts. This novel disease gene may have previously evaded detection because of high GC content, consequential low coverage and computational difficulties associated with robustly detecting repeat-expansions. Reviewing previously unsolved exomes using lower QC filters may generate further diagnoses.

De Novo SOX6 Variants Cause a Neurodevelopmental Syndrome Associated with ADHD, Craniosynostosis, and Osteochondromas

SOX6 belongs to a family of 20 SRY-related HMG-box-containing (SOX) genes that encode transcription factors controlling cell fate and differentiation in many developmental and adult processes. For SOX6, these processes include, but are not limited to, neurogenesis and skeletogenesis. Variants in half of the SOX genes have been shown to cause severe developmental and adult syndromes, referred to as SOXopathies. We here provide evidence that SOX6 variants also cause a SOXopathy. Using clinical and genetic data, we identify 19 individuals harboring various types of SOX6 alterations and exhibiting developmental delay and/or intellectual disability; the individuals are from 17 unrelated families. Additional, inconstant features include attention-deficit/hyperactivity disorder (ADHD), autism, mild facial dysmorphism, craniosynostosis, and multiple osteochondromas. All variants are heterozygous. Fourteen are de novo, one is inherited from a mosaic father, and four offspring from two families have a paternally inherited variant. Intragenic microdeletions, balanced structural rearrangements, frameshifts, and nonsense variants are predicted to inactivate the SOX6 variant allele. Four missense variants occur in residues and protein regions highly conserved evolutionarily. These variants are not detected in the gnomAD control cohort, and the amino acid substitutions are predicted to be damaging. Two of these variants are located in the HMG domain and abolish SOX6 transcriptional activity in vitro. No clear genotype-phenotype correlations are found. Taken together, these findings concur that SOX6 haploinsufficiency leads to a neurodevelopmental SOXopathy that often includes ADHD and abnormal skeletal and other features.

Publisher Correction: c-Maf controls immune responses by regulating disease-specific gene networks and repressing IL-2 in CD4+ T cells

In the version of this article initially published, the Supplementary Data file was an incorrect version. The correct version is now provided. The error has been corrected in the HTML and PDF version of the article.

Low plasma haptoglobin is a risk factor for life-threatening childhood severe malarial anemia and not an exclusive consequence of hemolysis

Severe Malarial Anemia (SMA), a life-threatening childhood Plasmodium falciparum malaria syndrome requiring urgent blood transfusion, exhibits inflammatory and hemolytic pathology. Differentiating between hypo-haptoglobinemia due to hemolysis or that of genetic origin is key to understand SMA pathogenesis. We hypothesized that while malaria-induced hypo-haptoglobinemia should reverse at recovery, that of genetic etiology should not. We carried-out a case-control study of children living under hyper-endemic holoendemic malaria burden in the sub-Saharan metropolis of Ibadan, Nigeria. We show that hypo-haptoglobinemia is a risk factor for childhood SMA and not solely due to intravascular hemolysis from underlying schizogony. In children presenting with SMA, hypo-haptoglobinemia remains through convalescence to recovery suggesting a genetic cause. We identified a haptoglobin gene variant, rs12162087 (g.-1203G > A, frequency = 0.67), to be associated with plasma haptoglobin levels (p = 8.5 × 10^-6). The Homo-Var:(AA) is associated with high plasma haptoglobin while the reference Homo-Ref:(GG) is associated with hypo-haptoglobinemia (p = 2.3 × 10^-6). The variant is associated with SMA, with the most support for a risk effect for Homo-Ref genotype. Our insights on regulatory haptoglobin genotypes and hypo-haptoglobinemia suggest that haptoglobin screening could be part of risk-assessment algorithms to prevent rapid disease progression towards SMA in regions with no-access to urgent blood transfusion where SMA accounts for high childhood mortality rates.

c-Maf controls immune responses by regulating disease-specific gene networks and repressing IL-2 in CD4+ T cells

The transcription factor c-Maf induces the anti-inflammatory cytokine IL-10 in CD4+ T cells in vitro. However, the global effects of c-Maf on diverse immune responses in vivo are unknown. Here we found that c-Maf regulated IL-10 production in CD4+ T cells in disease models involving the TH1 subset of helper T cells (malaria), TH2 cells (allergy) and TH17 cells (autoimmunity) in vivo. Although mice with c-Maf deficiency targeted to T cells showed greater pathology in TH1 and TH2 responses, TH17 cell-mediated pathology was reduced in this context, with an accompanying decrease in TH17 cells and increase in Foxp3+ regulatory T cells. Bivariate genomic footprinting elucidated the c-Maf transcription-factor network, including enhanced activity of NFAT; this led to the identification and validation of c-Maf as a negative regulator of IL-2. The decreased expression of the gene encoding the transcription factor RORγt (Rorc) that resulted from c-Maf deficiency was dependent on IL-2, which explained the in vivo observations. Thus, c-Maf is a positive and negative regulator of the expression of cytokine-encoding genes, with context-specific effects that allow each immune response to occur in a controlled yet effective manner.

Deterministic Evolutionary Trajectories Influence Primary Tumor Growth: TRACERx Renal

The evolutionary features of clear-cell renal cell carcinoma (ccRCC) have not been systematically studied to date. We analyzed 1,206 primary tumor regions from 101 patients recruited into the multi-center prospective study, TRACERx Renal. We observe up to 30 driver events per tumor and show that subclonal diversification is associated with known prognostic parameters. By resolving the patterns of driver event ordering, co-occurrence, and mutual exclusivity at clone level, we show the deterministic nature of clonal evolution. ccRCC can be grouped into seven evolutionary subtypes, ranging from tumors characterized by early fixation of multiple mutational and copy number drivers and rapid metastases to highly branched tumors with >10 subclonal drivers and extensive parallel evolution associated with attenuated progression. We identify genetic diversity and chromosomal complexity as determinants of patient outcome. Our insights reconcile the variable clinical behavior of ccRCC and suggest evolutionary potential as a biomarker for both intervention and surveillance.

Corrigendum: Phylogenetic ctDNA analysis depicts early-stage lung cancer evolution

This corrects the article DOI: 10.1038/nature22364.

Aberrant ribonucleotide incorporation and multiple deletions in mitochondrial DNA of the murine MPV17 disease model

All DNA polymerases misincorporate ribonucleotides despite their preference for deoxyribonucleotides, and analysis of cultured cells indicates that mammalian mitochondrial DNA (mtDNA) tolerates such replication errors. However, it is not clear to what extent misincorporation occurs in tissues, or whether this plays a role in human disease. Here, we show that mtDNA of solid tissues contains many more embedded ribonucleotides than that of cultured cells, consistent with the high ratio of ribonucleotide to deoxynucleotide triphosphates in tissues, and that riboadenosines account for three-quarters of them. The pattern of embedded ribonucleotides changes in a mouse model of Mpv17 deficiency, which displays a marked increase in rGMPs in mtDNA. However, while the mitochondrial dGTP is low in the Mpv17^−/− liver, the brain shows no change in the overall dGTP pool, leading us to suggest that Mpv17 determines the local concentration or quality of dGTP. Embedded rGMPs are expected to distort the mtDNA and impede its replication, and elevated rGMP incorporation is associated with early-onset mtDNA depletion in liver and late-onset multiple deletions in brain of Mpv17^−/− mice. These findings suggest aberrant ribonucleotide incorporation is a primary mtDNA abnormality that can result in pathology.

Molecular evolution of the vertebrate blood coagulation network

In mammalian blood coagulation 5 proteases, factor VII (FVII), factor IX (FIX), factor X (FX), protein C (PC) and prothrombin act with two cofactors factor V and factor VIII to control the generation of fibrin. Biochemical evidence and molecular cloning data have previously indicated that blood coagulation involving tissue factor, prothrombin and fibrinogen is present in all vertebrates. Using degenerate RT-PCR we have isolated and characterized novel cDNAs with sequence identity to the blood coagulation serine proteases and cofactors from chicken and the puffer fish (Fugu rubripes). Sequence alignments, phylogenetic and comparative sequence analysis all support the existence of the Gla-EGF1-EGF2-SP domain serine proteases FVII, FIX, FX, PC and the A1-A2-B-A3-C1-C2 domain protein cofactors FV and FVIII in these species. These results strongly suggest that the blood coagulation network is present in all jawed vertebrates and evolved before the divergence of tetrapods and teleosts over 430 million years ago; and that vertebrate blood coagulation may have benefited from two rounds of gene or whole genome duplication. Sequences identified in Fugu coding for additional FVII-like, FIX-like and PC-like sequences support the possibility of further tandem and large-scale duplications in teleosts. Comparative sequence analyses of amino acid residues in the active site region suggest these additional sequences have evolved new and as yet unknown functions.

Supplementary information to this article available at both http://europium.csc.mrc.ac.uk and www.thrombosis-online.com

Phylogenetic ctDNA analysis depicts early-stage lung cancer evolution

The early detection of relapse following primary surgery for non-small-cell lung cancer and the characterization of emerging subclones, which seed metastatic sites, might offer new therapeutic approaches for limiting tumour recurrence. The ability to track the evolutionary dynamics of early-stage lung cancer non-invasively in circulating tumour DNA (ctDNA) has not yet been demonstrated. Here we use a tumour-specific phylogenetic approach to profile the ctDNA of the first 100 TRACERx (Tracking Non-Small-Cell Lung Cancer Evolution Through Therapy (Rx)) study participants, including one patient who was also recruited to the PEACE (Posthumous Evaluation of Advanced Cancer Environment) post-mortem study. We identify independent predictors of ctDNA release and analyse the tumour-volume detection limit. Through blinded profiling of postoperative plasma, we observe evidence of adjuvant chemotherapy resistance and identify patients who are very likely to experience recurrence of their lung cancer. Finally, we show that phylogenetic ctDNA profiling tracks the subclonal nature of lung cancer relapse and metastasis, providing a new approach for ctDNA-driven therapeutic studies.

Evolution of lineage-specific functions in ancient cis-regulatory modules

Morphological evolution is driven both by coding sequence variation and by changes in regulatory sequences. However, how cis-regulatory modules (CRMs) evolve to generate entirely novel expression domains is largely unknown. Here, we reconstruct the evolutionary history of a lens enhancer located within a CRM that not only predates the lens, a vertebrate innovation, but bilaterian animals in general. Alignments of orthologous sequences from different deuterostomes sub-divide the CRM into a deeply conserved core and a more divergent flanking region. We demonstrate that all deuterostome flanking regions, including invertebrate sequences, activate gene expression in the zebrafish lens through the same ancient cluster of activator sites. However, levels of gene expression vary between species due to the presence of repressor motifs in flanking region and core. These repressor motifs are responsible for the relatively weak enhancer activity of tetrapod flanking regions. Ray-finned fish, however, have gained two additional lineage-specific activator motifs which in combination with the ancient cluster of activators and the core constitute a potent lens enhancer. The exploitation and modification of existing regulatory potential in flanking regions but not in the highly conserved core might represent a more general model for the emergence of novel regulatory functions in complex CRMs.

MPV17 Loss Causes Deoxynucleotide Insufficiency and Slow DNA Replication in Mitochondria

MPV17 is a mitochondrial inner membrane protein whose dysfunction causes mitochondrial DNA abnormalities and disease by an unknown mechanism. Perturbations of deoxynucleoside triphosphate (dNTP) pools are a recognized cause of mitochondrial genomic instability; therefore, we determined DNA copy number and dNTP levels in mitochondria of two models of MPV17 deficiency. In Mpv17 ablated mice, liver mitochondria showed substantial decreases in the levels of dGTP and dTTP and severe mitochondrial DNA depletion, whereas the dNTP pool was not significantly altered in kidney and brain mitochondria that had near normal levels of DNA. The shortage of mitochondrial dNTPs in Mpv17^-/- liver slows the DNA replication in the organelle, as evidenced by the elevated level of replication intermediates. Quiescent fibroblasts of MPV17-mutant patients recapitulate key features of the primary affected tissue of the Mpv17^-/- mice, displaying virtual absence of the protein, decreased dNTP levels and mitochondrial DNA depletion. Notably, the mitochondrial DNA loss in the patients’ quiescent fibroblasts was prevented and rescued by deoxynucleoside supplementation. Thus, our study establishes dNTP insufficiency in the mitochondria as the cause of mitochondrial DNA depletion in MPV17 deficiency, and identifies deoxynucleoside supplementation as a potential therapeutic strategy for MPV17-related disease. Moreover, changes in the expression of factors involved in mitochondrial deoxynucleotide homeostasis indicate a remodeling of nucleotide metabolism in MPV17 disease models, which suggests mitochondria lacking functional MPV17 have a restricted purine mitochondrial salvage pathway.

Mitochondrial DNA depletion syndrome (MDS) is a genetically heterogeneous condition characterized by a decrease of mitochondrial DNA (mtDNA) copy number and decreased activities of respiratory chain enzymes. Depletion of mtDNA has been associated with mutations in several genes, which encode either proteins directly involved in mtDNA replication or factors regulating the homeostasis of the mitochondrial deoxynucleotide pool. However, for some genes the mechanism linking mutations and mtDNA depletion is not known. One such gene is MPV17, whose loss-of-function causes mtDNA abnormalities in human, mouse and yeast. Here we show that MPV17 dysfunction leads to a shortage of the precursors for DNA synthesis in the mitochondria, slowing DNA replication in the organelle. Not only does mtDNA copy number correlate with dNTP pool size in both mouse tissues and human cells, deoxynucleoside supplementation of the growth medium prevents depletion and restores mtDNA copy number in quiescent MPV17-deficient cells. Hence, our study links MPV17 deficiency, insufficiency of mitochondrial dNTPs, and slow replication in mitochondria to depletion of mtDNA manifesting in the human disease, and places MPV17-related disease firmly in the category of mtDNA disorders caused by deoxynucleotide perturbation. The prevention and reversal of mtDNA loss in MPV17 patient-derived cells identifies potential therapeutic strategy for a currently untreatable disease.

Complete re-sequencing of a 2Mb topological domain encompassing the FTO/IRXB genes identifies a novel obesity-associated region upstream of IRX5

BACKGROUND

Association studies have identified a number of loci that contribute to an increased body mass index (BMI), the strongest of which is in the first intron of the FTO gene on human chromosome 16q12.2. However, this region is both non-coding and under strong linkage disequilibrium, making it recalcitrant to functional interpretation. Furthermore, the FTO gene is located within a complex cis-regulatory landscape defined by a topologically associated domain that includes the IRXB gene cluster, a trio of developmental regulators. Consequently, at least three genes in this interval have been implicated in the aetiology of obesity.

METHODS

Here, we sequence a 2 Mb region encompassing the FTO, RPGRIP1L and IRXB cluster genes in 284 individuals from a well-characterised study group of Danish men containing extremely overweight young adults and controls. We further replicate our findings both in an expanded male cohort and an independent female study group. Finally, we compare our variant data with a previous study describing IRX3 and FTO interactions in this region.

RESULTS

We obtain deep coverage across the entire region, allowing accurate and unequivocal determination of almost every single nucleotide polymorphism and short insertion/deletion. As well as confirming previous findings across the interval, we identify a further novel age-dependent association upstream of IRX5 that imposes a similar burden on BMI to the FTO locus.

CONCLUSIONS

Our findings are consistent with the hypothesis that chromatin architectures play a role in regulating gene expression levels across topological domains while our targeted sequence approach represents a widely applicable methodology for high-resolution analysis of regional variation across candidate genomic loci.

Identification and functional characterization of novel transcriptional enhancers involved in regulating human GLI3 expression during early development

The zinc-finger transcription factor GLI3 acts as a primary transducer of Sonic hedgehog (Shh) signaling in a context-dependent combinatorial fashion. GLI3 participates in the patterning and growth of many organs, including the central nervous system (CNS) and limbs. Previously, we reported a subset of human intronic cis-regulators controlling many known aspects of endogenous Gli3 expression in mouse and zebrafish. Here we demonstrate in a transgenic zebrafish assay the potential of two novel tetrapod-teleost conserved non-coding elements (CNEs) docking within GLI3 intronic intervals (intron 3 and 4) to induce reporter gene expression at known sites of endogenous Gli3 transcription in embryonic domains such as the central nervous system (CNS) and limbs. Interestingly, the cell culture based assays reveal harmony with the context dependent dual nature of intra-GLI3 conserved elements. Furthermore, a transgenic zebrafish assay of previously reported limb-specific GLI3 transcriptional enhancers (previously tested in mice and chicken limb buds) induced reporter gene expression in zebrafish blood precursor cells and notochord instead of fin. These results demonstrate that the appendage-specific activity of a subset of GLI3-associated enhancers might be a tetrapod innovation. Taken together with our recent data, these results suggest that during the course of vertebrate evolution Gli3 expression control acquired a complex cis-regulatory landscape for spatiotemporal patterning of CNS and limbs. Comparative data from fish and mice suggest that the functional aspects of a subset of these cis-regulators have diverged significantly between these two lineages.

A Simple Predictive Enhancer Syntax for Hindbrain Patterning Is Conserved in Vertebrate Genomes

Background

Determining the function of regulatory elements is fundamental for our understanding of development, disease and evolution. However, the sequence features that mediate these functions are often unclear and the prediction of tissue-specific expression patterns from sequence alone is non-trivial. Previous functional studies have demonstrated a link between PBX-HOX and MEIS/PREP binding interactions and hindbrain enhancer activity, but the defining grammar of these sites, if any exists, has remained elusive.

Results

Here, we identify a shared sequence signature (syntax) within a heterogeneous set of conserved vertebrate hindbrain enhancers composed of spatially co-occurring PBX-HOX and MEIS/PREP transcription factor binding motifs. We use this syntax to accurately predict hindbrain enhancers in 89% of cases (67/75 predicted elements) from a set of conserved non-coding elements (CNEs). Furthermore, mutagenesis of the sites abolishes activity or generates ectopic expression, demonstrating their requirement for segmentally restricted enhancer activity in the hindbrain. We refine and use our syntax to predict over 3,000 hindbrain enhancers across the human genome. These sequences tend to be located near developmental transcription factors and are enriched in known hindbrain activating elements, demonstrating the predictive power of this simple model.

Conclusion

Our findings support the theory that hundreds of CNEs, and perhaps thousands of regions across the human genome, function to coordinate gene expression in the developing hindbrain. We speculate that deeply conserved sequences of this kind contributed to the co-option of new genes into the hindbrain gene regulatory network during early vertebrate evolution by linking patterns of hox expression to downstream genes involved in segmentation and patterning, and evolutionarily newer instances may have continued to contribute to lineage-specific elaboration of the hindbrain.

Purifying selection in deeply conserved human enhancers is more consistent than in coding sequences

Comparison of polymorphism at synonymous and non-synonymous sites in protein-coding DNA can provide evidence for selective constraint. Non-coding DNA that forms part of the regulatory landscape presents more of a challenge since there is not such a clear-cut distinction between sites under stronger and weaker selective constraint. Here, we consider putative regulatory elements termed Conserved Non-coding Elements (CNEs) defined by their high level of sequence identity across all vertebrates. Some mutations in these regions have been implicated in developmental disorders; we analyse CNE polymorphism data to investigate whether such deleterious effects are widespread in humans. Single nucleotide variants from the HapMap and 1000 Genomes Projects were mapped across nearly 2000 CNEs. In the 1000 Genomes data we find a significant excess of rare derived alleles in CNEs relative to coding sequences; this pattern is absent in HapMap data, apparently obscured by ascertainment bias. The distribution of polymorphism within CNEs is not uniform; we could identify two categories of sites by exploiting deep vertebrate alignments: stretches that are non-variant, and those that have at least one substitution. The conserved category has fewer polymorphic sites and a greater excess of rare derived alleles, which can be explained by a large proportion of sites under strong purifying selection within humans--higher than that for non-synonymous sites in most protein coding regions, and comparable to that at the strongly conserved trans-dev genes. Conversely, the more evolutionarily labile CNE sites have an allele frequency distribution not significantly different from non-synonymous sites. Future studies should exploit genome-wide re-sequencing to obtain better coverage in selected non-coding regions, given the likelihood that mutations in evolutionarily conserved enhancer sequences are deleterious. Discovery pipelines should validate non-coding variants to aid in identifying causal and risk-enhancing variants in complex disorders, in contrast to the current focus on exome sequencing.

Dissection of a Ciona regulatory element reveals complexity of cross-species enhancer activity

Vertebrate genomes share numerous conserved non-coding elements, many of which function as enhancer elements and are hypothesised to be under evolutionary constraint due to a need to be bound by combinations of sequence-specific transcription factors. In contrast, few such conserved elements can be detected between vertebrates and their closest invertebrate relatives. Despite this lack of sequence identity, cross-species transgenesis has identified some cases where non-coding DNA from invertebrates drives reporter gene expression in transgenic vertebrates in patterns reminiscent of the expression of vertebrate orthologues. Such instances are presumed to reflect the presence of conserved suites of binding sites in the regulatory regions of invertebrate and vertebrate orthologues, such that both regulatory elements can correctly interpret the trans-activating environment. Shuffling of binding sites has been suggested to lie behind loss of sequence conservation; however this has not been experimentally tested. Here we examine the underlying basis of enhancer activity for the Ciona intestinalis βγ-crystallin gene, which drives expression in the lens of transgenic vertebrates despite the Ciona lineage predating the evolution of the lens. We construct an interactive gene regulatory network (GRN) for vertebrate lens development, allowing network interactions to be robustly catalogued and conserved network components and features to be identified. We show that a small number of binding motifs are necessary for Ciona βγ-crystallin expression, and narrow down the likely factors that bind to these motifs. Several of these overlap with the conserved core of the vertebrate lens GRN, implicating these sites in cross species function. However when we test these motifs in a transgenic vertebrate they prove to be dispensable for reporter expression in the lens. These results show that current models depicting cross species enhancer function as dependent on conserved binding sites can be overly simplistic, with sound evolutionary inference requiring detailed dissection of underlying mechanisms.

•

Analysis of binding motifs in a Ciona enhancer that also works in vertebrate lens.

•

Establishment of candidate transcription factors that may regulate this enhancer.

•

Construction of a curated, interactive gene regulatory network of lens development.

•

Public accessibility of this via a dedicated web site.

•

Experimental test of binding motif function in cross species transgenesis.

A reporter assay in lamprey embryos reveals both functional conservation and elaboration of vertebrate enhancers

The sea lamprey is an important model organism for investigating the evolutionary origins of vertebrates. As more vertebrate genome sequences are obtained, evolutionary developmental biologists are becoming increasingly able to identify putative gene regulatory elements across the breadth of the vertebrate taxa. The identification of these regions makes it possible to address how changes at the genomic level have led to changes in developmental gene regulatory networks and ultimately to the evolution of morphological diversity. Comparative genomics approaches using sea lamprey have already predicted a number of such regulatory elements in the lamprey genome. Functional characterisation of these sequences and other similar elements requires efficient reporter assays in lamprey. In this report, we describe the development of a transient transgenesis method for lamprey embryos. Focusing on conserved non-coding elements (CNEs), we use this method to investigate their functional conservation across the vertebrate subphylum. We find instances of both functional conservation and lineage-specific functional evolution of CNEs across vertebrates, emphasising the utility of functionally testing homologous CNEs in their host species.

Sequencing of the sea lamprey (Petromyzon marinus) genome provides insights into vertebrate evolution

Lampreys are representatives of an ancient vertebrate lineage that diverged from our own ∼500 million years ago. By virtue of this deeply shared ancestry, the sea lamprey (P. marinus) genome is uniquely poised to provide insight into the ancestry of vertebrate genomes and the underlying principles of vertebrate biology. Here, we present the first lamprey whole-genome sequence and assembly. We note challenges faced owing to its high content of repetitive elements and GC bases, as well as the absence of broad-scale sequence information from closely related species. Analyses of the assembly indicate that two whole-genome duplications likely occurred before the divergence of ancestral lamprey and gnathostome lineages. Moreover, the results help define key evolutionary events within vertebrate lineages, including the origin of myelin-associated proteins and the development of appendages. The lamprey genome provides an important resource for reconstructing vertebrate origins and the evolutionary events that have shaped the genomes of extant organisms.

Capturing the regulatory interactions of eukaryote genomes

A key finding from early genomics research is the remarkable consistency in the number of protein-coding regions across diverse species. This has led many researchers to look to the cis-regulatory elements of genes as the fundamental influence behind evolving gene function and subsequent species diversification. Historically, since these elements are often located in vast intergenic and intronic regions of the genome, their identification has been recalcitrant. Now, with the deluge of whole-genome data from representatives of numerous eukaryotic lineages, various approaches have enabled us to begin to recognize features that characterize regulatory regions of the genome. Here we endeavour to collate these approaches in order to give an overview of the complexities involved in extrapolating regulatory signatures. The resource provided by the escalating richness of whole-genome datasets enables more sophisticated modelling of these regulatory signatures yet at the same time introduces increasing potential for noise. While we are only at the advent of making these discoveries, the next decade promises to be a very exciting and rewarding time for genome researchers.

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

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

Ancient Pbx-Hox signatures define hundreds of vertebrate developmental enhancers

BACKGROUND

Gene regulation through cis-regulatory elements plays a crucial role in development and disease. A major aim of the post-genomic era is to be able to read the function of cis-regulatory elements through scrutiny of their DNA sequence. Whilst comparative genomics approaches have identified thousands of putative regulatory elements, our knowledge of their mechanism of action is poor and very little progress has been made in systematically de-coding them.

RESULTS

Here, we identify ancient functional signatures within vertebrate conserved non-coding elements (CNEs) through a combination of phylogenetic footprinting and functional assay, using genomic sequence from the sea lamprey as a reference. We uncover a striking enrichment within vertebrate CNEs for conserved binding-site motifs of the Pbx-Hox hetero-dimer. We further show that these predict reporter gene expression in a segment specific manner in the hindbrain and pharyngeal arches during zebrafish development.

CONCLUSIONS

These findings evoke an evolutionary scenario in which many CNEs evolved early in the vertebrate lineage to co-ordinate Hox-dependent gene-regulatory interactions that pattern the vertebrate head. In a broader context, our evolutionary analyses reveal that CNEs are composed of tightly linked transcription-factor binding-sites (TFBSs), which can be systematically identified through phylogenetic footprinting approaches. By placing a large number of ancient vertebrate CNEs into a developmental context, our findings promise to have a significant impact on efforts toward de-coding gene-regulatory elements that underlie vertebrate development, and will facilitate building general models of regulatory element evolution.

STATISTICAL INFORMATION CHARACTERIZATION OF CONSERVED NON-CODING ELEMENTS IN VERTEBRATES

Recently, a set of highly conserved non-coding elements (CNEs) has been derived from a comparison between the genomes of the puffer fish, Takifugu or Fugu rubripes, and man. In order to facilitate the identification of these conserved elements in silico, we characterize them by a number of statistical features.

We found a pronounced information pattern around CNE borders; although the CNEs themselves are AT rich and have high entropy (complexity), they are flanked by GC-rich regions of low entropy (complexity). We also identified the most abundant motifs within and around of CNEs, and identified those that group around their borders. Like in human promoter regions, the TBP, NF-Y and some other binding motifs are clustered around CNE boundaries, which may suggest a possible transcription regulatory function of CNEs.

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

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

SAND, a new protein family: from nucleic acid to protein structure and function prediction

As a result of genome, EST and cDNA sequencing projects, there are huge numbers of predicted and/or partially characterised protein sequences compared with a relatively small number of proteins with experimentally determined function and structure. Thus, there is a considerable attention focused on the accurate prediction of gene function and structure from sequence by using bioinformatics. In the course of our analysis of genomic sequence from Fugu rubripes, we identified a novel gene, SAND, with significant sequence identity to hypothetical proteins predicted in Saccharomyces cerevisiae, Schizosaccharomyces pombe, Caenorhabditis elegans, a Drosophila melanogaster gene, and mouse and human cDNAs. Here we identify a further SAND homologue in human and Arabidopsis thaliana by use of standard computational tools. We describe the genomic organisation of SAND in these evolutionarily divergent species and identify sequence homologues from EST database searches confirming the expression of SAND in over 20 different eukaryotes. We confirm the expression of two different SAND paralogues in mammals and determine expression of one SAND in other vertebrates and eukaryotes. Furthermore, we predict structural properties of SAND, and characterise conserved sequence motifs in this protein family.

The PAX258 gene subfamily: a comparative perspective

Whole genome duplication events are thought to have substantially contributed to organismal complexity, largely via divergent transcriptional regulation. Members of the vertebrate PAX2, PAX5 and PAX8 gene subfamily derived from an ancient class of paired box genes and arose from such whole genome duplication events. These genes are critical in establishing the midbrain-hindbrain boundary, specifying interneuron populations and for eye, ear and kidney development. Also PAX2 has adopted a unique role in pancreas development, whilst PAX5 is essential for early B-cell differentiation. The contribution of PAX258 genes to their collective role has diverged across paralogues and the animal lineages, resulting in a complex wealth of literature. It is now timely to provide a comprehensive comparative overview of these genes and their ancient and divergent roles. We also discuss their fundamental place within gene regulatory networks and the likely influence of cis-regulatory elements over their differential roles during early animal development.

Early evolution of conserved regulatory sequences associated with development in vertebrates

Comparisons between diverse vertebrate genomes have uncovered thousands of highly conserved non-coding sequences, an increasing number of which have been shown to function as enhancers during early development. Despite their extreme conservation over 500 million years from humans to cartilaginous fish, these elements appear to be largely absent in invertebrates, and, to date, there has been little understanding of their mode of action or the evolutionary processes that have modelled them. We have now exploited emerging genomic sequence data for the sea lamprey, Petromyzon marinus, to explore the depth of conservation of this type of element in the earliest diverging extant vertebrate lineage, the jawless fish (agnathans). We searched for conserved non-coding elements (CNEs) at 13 human gene loci and identified lamprey elements associated with all but two of these gene regions. Although markedly shorter and less well conserved than within jawed vertebrates, identified lamprey CNEs are able to drive specific patterns of expression in zebrafish embryos, which are almost identical to those driven by the equivalent human elements. These CNEs are therefore a unique and defining characteristic of all vertebrates. Furthermore, alignment of lamprey and other vertebrate CNEs should permit the identification of persistent sequence signatures that are responsible for common patterns of expression and contribute to the elucidation of the regulatory language in CNEs. Identifying the core regulatory code for development, common to all vertebrates, provides a foundation upon which regulatory networks can be constructed and might also illuminate how large conserved regulatory sequence blocks evolve and become fixed in genomic DNA.

Recent comparative analyses of vertebrate genomes has resulted in the identification of highly conserved non-coding sequences near genes that coordinate early development. Many of these sequences can activate gene expression and are thought to be important regulatory elements. Surprisingly, a large set of these long, near-identical sequences is found in every jawed vertebrate, including sharks, yet almost completely absent in non-vertebrates. This study looks for this set of sequences in the lamprey, a representative of our most distant vertebrate relatives, in order to determine when and how such a large set of important non-coding regulatory sequences became established in the genome. Although the lamprey divergence is only a little older than the divergence of cartilaginous fish (including sharks), relatively few, and considerably shorter, conserved non-coding sequences are identifiable. Nevertheless, these shorter lamprey sequences are capable of driving gene expression in a precise spatial pattern in zebrafish embryos in the same way as the equivalent human elements. This analysis has shed light on the emergence of these regulatory sequences during early vertebrate evolution, at a time of whole-genome duplications and considerable morphological variation, consistent with a role for these sequences in directing gene regulatory networks for vertebrate development.

Organization of conserved elements near key developmental regulators in vertebrate genomes

Sequence conservation has traditionally been used as a means to target functional regions of complex genomes. In addition to its use in identifying coding regions of genes, the recent availability of whole genome data for a number of vertebrates has permitted high-resolution analyses of the noncoding "dark matter" of the genome. This has resulted in the identification of a large number of highly conserved sequence elements that appear to be preserved in all bony vertebrates. Further positional analysis of these conserved noncoding elements (CNEs) in the genome demonstrates that they cluster around genes involved in developmental regulation. This chapter describes the identification and characterization of these elements, with particular reference to their composition and organization.

Persistence of duplicated PAC1 receptors in the teleost, Sparus auratus

Background:

Duplicated genes are common in vertebrate genomes. Their persistence is assumed to be either a consequence of gain of novel function (neofunctionalisation) or partitioning of the function of the ancestral molecule (sub-functionalisation). Surprisingly few studies have evaluated the extent of such modifications despite the numerous duplicated receptor and ligand genes identified in vertebrate genomes to date. In order to study the importance of function in the maintenance of duplicated genes, sea bream (Sparus auratus) PAC₁ receptors, sequence homologues of the mammalian receptor specific for PACAP (Pituitary Adenylate Cyclase-Activating Polypeptide), were studied. These receptors belong to family 2 GPCRs and most of their members are duplicated in teleosts although the reason why both persist in the genome is unknown.

Results:

Duplicate sea bream PACAP receptor genes (sbPAC₁A and sbPAC₁B), members of family 2 GPCRs, were isolated and share 77% amino acid sequence identity. RT-PCR with specific primers for each gene revealed that they have a differential tissue distribution which overlaps with the distribution of the single mammalian receptor. Furthermore, in common with mammals, the teleost genes undergo alternative splicing and a PAC₁Ahop1 isoform has been characterised. Duplicated orthologous receptors have also been identified in other teleost genomes and their distribution profile suggests that function may be species specific. Functional analysis of the paralogue sbPAC₁s in Cos7 cells revealed that they are strongly stimulated in the presence of mammalian PACAP₂₇ and PACAP₃₈ and far less with VIP (Vasoactive Intestinal Peptide). The sbPAC₁ receptors are equally stimulated (LOG_EC50 values for maximal cAMP production) in the presence of PACAP₂₇ (-8.74 ± 0.29 M and -9.15 ± 0.21 M, respectively for sbPAC₁A and sbPAC₁B, P > 0.05) and PACAP₃₈ (-8.54 ± 0.18 M and -8.92 ± 0.24 M, respectively for sbPAC₁A and sbPAC₁B, P > 0.05). Human VIP was found to stimulate sbPAC₁A (-7.23 ± 0.20 M) more strongly than sbPAC₁B (-6.57 ± 0.14 M, P < 0.05) and human secretin (SCT), which has not so far been identified in fish genomes, caused negligible stimulation of both receptors.

Conclusion:

The existence of functionally divergent duplicate sbPAC₁ receptors is in line with previously proposed theories about the origin and maintenance of duplicated genes. Sea bream PAC₁ duplicate receptors resemble the typical mammalian PAC₁, and PACAP peptides were found to be more effective than VIP in stimulating cAMP production, although sbPAC₁A was more responsive for VIP than sbPAC₁B. These results together with the highly divergent pattern of tissue distribution suggest that a process involving neofunctionalisation occurred after receptor duplication within the fish lineage and probably accounts for their persistence in the genome. The characterisation of further duplicated receptors and their ligands should provide insights into the evolution and function of novel protein-protein interactions associated with the vertebrate radiation.

Comparative genomics using Fugu reveals insights into regulatory subfunctionalization

Fish-mammal genomic alignments were used to compare over 800 conserved non-coding elements that associate with genes that have undergone fish-specific duplication and retention, revealing a pattern of element retention and loss between paralogs indicative of subfunctionalization.

Background

A major mechanism for the preservation of gene duplicates in the genome is thought to be mediated via loss or modification of cis-regulatory subfunctions between paralogs following duplication (a process known as regulatory subfunctionalization). Despite a number of gene expression studies that support this mechanism, no comprehensive analysis of regulatory subfunctionalization has been undertaken at the level of the distal cis-regulatory modules involved. We have exploited fish-mammal genomic alignments to identify and compare more than 800 conserved non-coding elements (CNEs) that associate with genes that have undergone fish-specific duplication and retention.

Results

Using the abundance of duplicated genes within the Fugu genome, we selected seven pairs of teleost-specific paralogs involved in early vertebrate development, each containing clusters of CNEs in their vicinity. CNEs present around each Fugu duplicated gene were identified using multiple alignments of orthologous regions between single-copy mammalian orthologs (representing the ancestral locus) and each fish duplicated region in turn. Comparative analysis reveals a pattern of element retention and loss between paralogs indicative of subfunctionalization, the extent of which differs between duplicate pairs. In addition to complete loss of specific regulatory elements, a number of CNEs have been retained in both regions but may be responsible for more subtle levels of subfunctionalization through sequence divergence.

Conclusion

Comparative analysis of conserved elements between duplicated genes provides a powerful approach for studying regulatory subfunctionalization at the level of the regulatory elements involved.

Divergent evolution of the myosin heavy chain gene family in fish and tetrapods: evidence from comparative genomic analysis

Myosin heavy chain genes (MYHs) are the most important functional domains of myosins, which are highly conserved throughout evolution. The human genome contains 15 MYHs, whereas the corresponding number in teleost appears to be much higher. Although teleosts comprise more than one-half of all vertebrate species, our knowledge of MYHs in teleosts is rather limited. A comprehensive analysis of the torafugu (Takifugu rubripes) genome database enabled us to detect at least 28 MYHs, almost twice as many as in humans. RT-PCR revealed that at least 16 torafugu MYH representatives (5 fast skeletal, 3 cardiac, 2 slow skeletal, 1 superfast, 2 smooth, and 3 nonmuscle types) are actually transcribed. Among these, MYH(M743-2) and MYH(M5) of fast and slow skeletal types, respectively, are expressed during development of torafugu embryos. Syntenic analysis reveals that torafugu fast skeletal MYHs are distributed across five genomic regions, three of which form clusters. Interestingly, while human fast skeletal MYHs form one cluster, its syntenic region in torafugu is duplicated, although each locus contains just a single MYH in torafugu. The results of the syntenic analysis were further confirmed by corresponding analysis of MYHs based on databases from Tetraodon, zebrafish, and medaka genomes. Phylogenetic analysis suggests that fast skeletal MYHs evolved independently in teleosts and tetrapods after fast skeletal MYHs had diverged from four ancestral MYHs.

Ultraconserved non-coding sequence element controls a subset of spatiotemporal GLI3 expression

The zinc-finger transcription factor GLI3 acts during vertebrate development in a combinatorial, context-dependent fashion as a primary transducer of sonic hedgehog (SHH) signaling. In humans, mutations affecting this key regulator of development are associated with GLI3-morphopathies, a group of congenital malformations in which forebrain and limb development are preferentially affected. We show that a non-coding element from intron two of GLI3, ultraconserved in mammals and highly conserved in the pufferfish Fugu, is a transcriptional enhancer. In transient transfection assays, it activates reporter gene transcription in human cell cultures expressing endogenous GLI3 but not in GLI3 negative cells. The identified enhancer element is predicted to contain conserved binding sites for transcription factors crucial for developmental steps in which GLI3 is involved. The regulatory potential of this element is conserved and was used to direct tissue-specific expression of a green fluorescent protein reporter gene in zebrafish embryos and of a beta-galactosidase reporter in transgenic mouse embryos. Time, location, and quantity of reporter gene expression are congruent with part of the pattern previously reported for endogenous GLI3 transcription.

CONDOR: a database resource of developmentally associated conserved non-coding elements

Background

Comparative genomics is currently one of the most popular approaches to study the regulatory architecture of vertebrate genomes. Fish-mammal genomic comparisons have proved powerful in identifying conserved non-coding elements likely to be distal cis-regulatory modules such as enhancers, silencers or insulators that control the expression of genes involved in the regulation of early development. The scientific community is showing increasing interest in characterizing the function, evolution and language of these sequences. Despite this, there remains little in the way of user-friendly access to a large dataset of such elements in conjunction with the analysis and the visualization tools needed to study them.

Description

Here we present CONDOR (COnserved Non-coDing Orthologous Regions) available at: . In an interactive and intuitive way the website displays data on > 6800 non-coding elements associated with over 120 early developmental genes and conserved across vertebrates. The database regularly incorporates results of ongoing in vivo zebrafish enhancer assays of the CNEs carried out in-house, which currently number ~100. Included and highlighted within this set are elements derived from duplication events both at the origin of vertebrates and more recently in the teleost lineage, thus providing valuable data for studying the divergence of regulatory roles between paralogs. CONDOR therefore provides a number of tools and facilities to allow scientists to progress in their own studies on the function and evolution of developmental cis-regulation.

Conclusion

By providing access to data with an approachable graphics interface, the CONDOR database presents a rich resource for further studies into the regulation and evolution of genes involved in early development.

Human GLI3 intragenic conserved non-coding sequences are tissue-specific enhancers

The zinc-finger transcription factor GLI3 is a key regulator of development, acting as a primary transducer of Sonic hedgehog (SHH) signaling in a combinatorial context dependent fashion controlling multiple patterning steps in different tissues/organs. A tight temporal and spatial control of gene expression is indispensable, however, cis-acting sequence elements regulating GLI3 expression have not yet been reported. We show that 11 ancient genomic DNA signatures, conserved from the pufferfish Takifugu (Fugu) rubripes to man, are distributed throughout the introns of human GLI3. They map within larger conserved non-coding elements (CNEs) that are found in the tetrapod lineage. Full length CNEs transiently transfected into human cell cultures acted as cell type specific enhancers of gene transcription. The regulatory potential of these elements is conserved and was exploited to direct tissue specific expression of a reporter gene in zebrafish embryos. Assays of deletion constructs revealed that the human-Fugu conserved sequences within the GLI3 intronic CNEs were essential but not sufficient for full-scale transcriptional activation. The enhancer activity of the CNEs is determined by a combinatorial effect of a core sequence conserved between human and teleosts (Fugu) and flanking tetrapod-specific sequences, suggesting that successive clustering of sequences with regulatory potential around an ancient, highly conserved nucleus might be a possible mechanism for the evolution of cis-acting regulatory elements.

Parallel evolution of conserved non-coding elements that target a common set of developmental regulatory genes from worms to humans

BACKGROUND

The human genome contains thousands of non-coding sequences that are often more conserved between vertebrate species than protein-coding exons. These highly conserved non-coding elements (CNEs) are associated with genes that coordinate development, and have been proposed to act as transcriptional enhancers. Despite their extreme sequence conservation in vertebrates, sequences homologous to CNEs have not been identified in invertebrates.

RESULTS

Here we report that nematode genomes contain an alternative set of CNEs that share sequence characteristics, but not identity, with their vertebrate counterparts. CNEs thus represent a very unusual class of sequences that are extremely conserved within specific animal lineages yet are highly divergent between lineages. Nematode CNEs are also associated with developmental regulatory genes, and include well-characterized enhancers and transcription factor binding sites, supporting the proposed function of CNEs as cis-regulatory elements. Most remarkably, 40 of 156 human CNE-associated genes with invertebrate orthologs are also associated with CNEs in both worms and flies.

CONCLUSION

A core set of genes that regulate development is associated with CNEs across three animal groups (worms, flies and vertebrates). We propose that these CNEs reflect the parallel evolution of alternative enhancers for a common set of developmental regulatory genes in different animal groups. This 're-wiring' of gene regulatory networks containing key developmental coordinators was probably a driving force during the evolution of animal body plans. CNEs may, therefore, represent the genomic traces of these 'hard-wired' core gene regulatory networks that specify the development of each alternative animal body plan.

Different words, same meaning: understanding the languages of the genome

There is undeniable value in using sequence comparison to identify putative regulatory sequences. However, a recent report has demonstrated that not all regulatory sequences are evolutionarily conserved. Cis-acting sequences around the RET gene, conserved in mammals but not in fish, are able to reproduce patterns of RET expression in zebrafish embryos. It is as yet unclear whether these sequences are 'below the radar' of current sequence alignment tools or whether their functional homology is not sequence based.