Comparative genomics of the Hedgehog loci in chordates and the origins of Shh regulatory novelties

Topologically Associating Domains and Regulatory Landscapes in Development, Evolution and Disease

Animal genomes are folded in topologically associating domains (TADs) that have been linked to the regulation of the genes they contain by constraining regulatory interactions between cis-regulatory elements and promoters. Therefore, TADs are proposed as structural scaffolds for the establishment of regulatory landscapes (RLs). In this review, we discuss recent advances in the connection between TADs and gene regulation, their relationship with gene RLs and their dynamics during development and differentiation. Moreover, we describe how restructuring TADs may lead to pathological conditions, which explains their high evolutionary conservation, but at the same time it provides a substrate for the emergence of evolutionary innovations that lay at the origin of vertebrates and other phylogenetic clades.

Ancient Genomic Regulatory Blocks Are a Source for Regulatory Gene Deserts in Vertebrates after Whole-Genome Duplications

We investigated how the two rounds of whole-genome duplication that occurred at the base of the vertebrate lineage have impacted ancient microsyntenic associations involving developmental regulators (known as genomic regulatory blocks, GRBs). We showed that the majority of GRBs identified in the last common ancestor of chordates have been maintained as a single copy in humans. We found evidence that dismantling of the duplicated GRB copies occurred early in vertebrate evolution often through the differential retention of the regulatory gene but loss of the bystander gene’s exonic sequences. Despite the large evolutionary scale, the presence of duplicated highly conserved noncoding regions provided unambiguous proof for this scenario for multiple ancient GRBs. Remarkably, the dismantling of ancient GRB duplicates has contributed to the creation of large gene deserts associated with regulatory genes in vertebrates, providing a potentially widespread mechanism for the origin of these enigmatic genomic traits.

Hedgehog-FGF signaling axis patterns anterior mesoderm during gastrulation

The mechanisms used by embryos to pattern tissues across their axes has fascinated developmental biologists since the founding of embryology. Here, using single-cell technology, we interrogate complex patterning defects and define a Hedgehog (Hh)-fibroblast growth factor (FGF) signaling axis required for anterior mesoderm lineage development during gastrulation. Single-cell transcriptome analysis of Hh-deficient mesoderm revealed selective deficits in anterior mesoderm populations, culminating in defects to anterior embryonic structures, including the pharyngeal arches, heart, and anterior somites. Transcriptional profiling of Hh-deficient mesoderm during gastrulation revealed disruptions to both transcriptional patterning of the mesoderm and FGF signaling for mesoderm migration. Mesoderm-specific Fgf4/Fgf8 double-mutants recapitulated anterior mesoderm defects and Hh-dependent GLI transcription factors modulated enhancers at FGF gene loci. Cellular migration defects during gastrulation induced by Hh pathway antagonism were mitigated by the addition of FGF4 protein. These findings implicate a multicomponent signaling hierarchy activated by Hh ligands from the embryonic node and executed by FGF signals in nascent mesoderm to control anterior mesoderm patterning.

Gene regulatory landscape of the sonic hedgehog locus in embryonic development

The organs of vertebrate species display a wide variety of morphology. A remaining challenge in evolutionary developmental biology is to elucidate how vertebrate lineages acquire distinct morphological features. Developmental programs are driven by spatiotemporal regulation of gene expression controlled by hundreds of thousands of cis-regulatory elements. Changes in the regulatory elements caused by the introduction of genetic variants can confer regulatory innovation that may underlie morphological novelties. Recent advances in sequencing technology have revealed a number of potential regulatory variants that can alter gene expression patterns. However, a limited number of studies demonstrate causal dependence between genetic and morphological changes. Regulation of Shh expression is a good model to understand how multiple regulatory elements organize tissue-specific gene expression patterns. This model also provides insights into how evolution of molecular traits, such as gene regulatory networks, lead to phenotypic novelty.

Temporospatial sonic hedgehog signalling is essential for neural crest-dependent patterning of the intrinsic tongue musculature

The tongue is a highly specialised muscular organ with a complex anatomy required for normal function. We have utilised multiple genetic approaches to investigate local temporospatial requirements for sonic hedgehog (SHH) signalling during tongue development. Mice lacking a Shh cis-enhancer, MFCS4 (ShhMFCS4/-), with reduced SHH in dorsal tongue epithelium have perturbed lingual septum tendon formation and disrupted intrinsic muscle patterning, with these defects reproduced following global Shh deletion from E10.5 in pCag-CreERTM; Shhflox/flox embryos. SHH responsiveness was diminished in local cranial neural crest cell (CNCC) populations in both mutants, with SHH targeting these cells through the primary cilium. CNCC-specific deletion of orofaciodigital syndrome 1 (Ofd1), which encodes a ciliary protein, in Wnt1-Cre; Ofdfl/Y mice led to a complete loss of normal myotube arrangement and hypoglossia. In contrast, mesoderm-specific deletion of Ofd1 in Mesp1-Cre; Ofdfl/Y embryos resulted in normal intrinsic muscle arrangement. Collectively, these findings suggest key temporospatial requirements for local SHH signalling in tongue development (specifically, lingual tendon differentiation and intrinsic muscle patterning through signalling to CNCCs) and provide further mechanistic insight into the tongue anomalies seen in patients with disrupted hedgehog signalling.

Prevalence of Mutation-Prone Microhomology-Mediated End Joining in a Chordate Lacking the c-NHEJ DNA Repair Pathway

Classical non-homologous end joining (c-NHEJ), a fundamental pathway that repairs double-strand breaks in DNA, is almost universal in eukaryotes and involves multiple proteins highly conserved from yeast to human [1]. The genes encoding these proteins were not detected in the genome of Oikopleura dioica, a new model system of tunicate larvaceans known for its very compact and highly rearranged genome [2-4]. After showing their absence in the genomes of six other larvacean species, the present study examined how O. dioica oocytes and embryos repair double-strand DNA breaks (DSBs), using two approaches: the injection of linearized plasmids, which resulted in their rapid end joining, and a newly established CRISPR Cas9 technique. In both cases, end joining merged short microhomologous sequences surrounding the break (mainly 4 bp long), thus inducing deletions larger than for the tunicate ascidian Ciona intestinalis and human cells. A relatively high frequency of nucleotide insertions was also observed. Finally, a survey of genomic indels supports the involvement of microhomology-mediated repair in natural conditions. Overall, O. dioica repairs DSBs as other organisms do when their c-NHEJ pathway is experimentally rendered deficient, using another mode of end joining with the same effect as alternative NHEJ (a-NHEJ) or microhomology-mediated end joining (MMEJ) [5-7]. We discuss how the exceptional loss of c-NHEJ and its replacement by a more mutation-prone mechanism may have contributed to reshaping this genome and even been advantageous under pressure for genome compaction.

Common genetic causes of holoprosencephaly are limited to a small set of evolutionarily conserved driver genes of midline development coordinated by TGF-β, hedgehog, and FGF signaling

Here, we applied targeted capture to examine 153 genes representative of all the major vertebrate developmental pathways among 333 probands to rank their relative significance as causes for holoprosencephaly (HPE). We now show that comparisons of variant transmission versus nontransmission among 136 HPE Trios indicates some reported genes now lack confirmation, while novel genes are implicated. Furthermore, we demonstrate that variation of modest intrinsic effect can synergize with these driver mutations as gene modifiers.

An evolutionary perspective of regulatory landscape dynamics in development and disease

The organization of animal genomes into topologically associating domains (TADs) provides a structural scaffold in which cis-regulatory elements (CREs) operate on their target genes. Determining the position of CREs and genes relative to TADs has become instrumental to trace gene expression changes during evolution and in diseases. Here we will review recent studies and discuss TADs as structural units with respect to their conservation and stability during genome reorganization. Furthermore, we describe how TAD restructuring contributed to morphological novelties during evolution but also their deleterious effects associated with disease. Despite considering TADs as structural units, the nested and dynamic scaffold within TADs contributes to tissue-specific gene expression, implying that such changes can also account for gene expression differences during evolution.

Signals from the brain and olfactory epithelium control shaping of the mammalian nasal capsule cartilage

Facial shape is the basis for facial recognition and categorization. Facial features reflect the underlying geometry of the skeletal structures. Here, we reveal that cartilaginous nasal capsule (corresponding to upper jaw and face) is shaped by signals generated by neural structures: brain and olfactory epithelium. Brain-derived Sonic Hedgehog (SHH) enables the induction of nasal septum and posterior nasal capsule, whereas the formation of a capsule roof is controlled by signals from the olfactory epithelium. Unexpectedly, the cartilage of the nasal capsule turned out to be important for shaping membranous facial bones during development. This suggests that conserved neurosensory structures could benefit from protection and have evolved signals inducing cranial cartilages encasing them. Experiments with mutant mice revealed that the genomic regulatory regions controlling production of SHH in the nervous system contribute to facial cartilage morphogenesis, which might be a mechanism responsible for the adaptive evolution of animal faces and snouts.

Holoprosencephaly in the genomics era

Holoprosencephaly (HPE) is the direct consequence of specific genetic and/or environmental insults interrupting the midline specification of the nascent forebrain. Such disturbances can lead to a broad range of phenotypic consequences for the brain and face in humans. This malformation sequence is remarkably common in utero (1 in 250 human fetuses), but 97% typically do not survive to birth. The precise molecular pathogenesis of HPE in these early human embryos remains largely unknown. Here, we outline our current understanding of the principal driving factors leading to HPE pathologies and elaborate our multifactorial integrated genomics approach. Overall, our understanding of the pathogenesis continues to become simpler, rather than more complicated. Genomic technologies now provide unprecedented insight into disease-associated variation, including the overall extent of genetic interactions (coding and noncoding) predicted to explain divergent phenotypes.

A conserved Shh cis-regulatory module highlights a common developmental origin of unpaired and paired fins

Despite their evolutionary, developmental, and functional importance the origin of vertebrate paired appendages remains uncertain. In mice, a single enhancer termed ZRS is solely responsible for Shh expression in limbs. Here, zebrafish and mouse transgenic assays trace the functional equivalence of ZRS across the gnathostome phylogeny. CRISPR/Cas9-mediated deletion of the medaka-ZRS and enhancer assays reveal the existence of ZRS shadow enhancers in both teleost and human genomes. Deletion of both ZRS and shadow ZRS abolish shh expression and completely truncate pectoral fin formation. Strikingly, deletion of ZRS results in an almost complete ablation of the dorsal fin. This finding indicates that a ZRS-Shh regulatory module is shared by paired and median fins, and that paired fins likely emerged by the co‐option of developmental programs established in the median fins of stem gnathostomes. Shh function was later reinforced in pectoral fin development with the recruitment of shadow enhancers, conferring additional robustness.

Topologically associating domains are ancient features that coincide with Metazoan clusters of extreme noncoding conservation

Developmental genes in metazoan genomes are surrounded by dense clusters of conserved noncoding elements (CNEs). CNEs exhibit unexplained extreme levels of sequence conservation, with many acting as developmental long-range enhancers. Clusters of CNEs define the span of regulatory inputs for many important developmental regulators and have been described previously as genomic regulatory blocks (GRBs). Their function and distribution around important regulatory genes raises the question of how they relate to 3D conformation of these loci. Here, we show that clusters of CNEs strongly coincide with topological organisation, predicting the boundaries of hundreds of topologically associating domains (TADs) in human and Drosophila. The set of TADs that are associated with high levels of noncoding conservation exhibit distinct properties compared to TADs devoid of extreme noncoding conservation. The close correspondence between extreme noncoding conservation and TADs suggests that these TADs are ancient, revealing a regulatory architecture conserved over hundreds of millions of years.

Metazoan genomes contain many clusters of conserved noncoding elements. Here, the authors provide evidence that these clusters coincide with distinct topologically associating domains in humans and Drosophila, revealing a conserved regulatory genomic architecture.

Cis-regulatory landscapes in development and evolution

The recent advances in our understanding of the 3D organization of the chromatin together with an almost unlimited ability to detect cis-regulatory elements genome-wide using different biochemical signatures has provided us with an unprecedented power to study gene regulation. It is now possible to profile the complete regulatory apparatus controlling the spatio-temporal expression of any given gene, the so-called gene Regulatory Landscapes (RLs). Here we review several studies over the last two years demonstrating the functional consequences of altering RL structure in development, disease and evolution. These works clearly show that a deep understanding of transcriptional regulation is no longer conceivable without considering the 3D modular organization of animal genomes.

Conserved Noncoding Elements in the Most Distant Genera of Cephalochordates: The Goldilocks Principle

Cephalochordates, the sister group of vertebrates + tunicates, are evolving particularly slowly. Therefore, genome comparisons between two congeners of Branchiostoma revealed so many conserved noncoding elements (CNEs), that it was not clear how many are functional regulatory elements. To more effectively identify CNEs with potential regulatory functions, we compared noncoding sequences of genomes of the most phylogenetically distant cephalochordate genera, Asymmetron and Branchiostoma, which diverged approximately 120-160 million years ago. We found 113,070 noncoding elements conserved between the two species, amounting to 3.3% of the genome. The genomic distribution, target gene ontology, and enriched motifs of these CNEs all suggest that many of them are probably cis-regulatory elements. More than 90% of previously verified amphioxus regulatory elements were re-captured in this study. A search of the cephalochordate CNEs around 50 developmental genes in several vertebrate genomes revealed eight CNEs conserved between cephalochordates and vertebrates, indicating sequence conservation over >500 million years of divergence. The function of five CNEs was tested in reporter assays in zebrafish, and one was also tested in amphioxus. All five CNEs proved to be tissue-specific enhancers. Taken together, these findings indicate that even though Branchiostoma and Asymmetron are distantly related, as they are evolving slowly, comparisons between them are likely optimal for identifying most of their tissue-specific cis-regulatory elements laying the foundation for functional characterizations and a better understanding of the evolution of developmental regulation in cephalochordates.

Cis-regulatory architecture of a brain signaling center predates the origin of chordates

Genomic approaches have predicted hundreds of thousands of tissue-specific cis-regulatory sequences, but the determinants critical to their function and evolutionary history are mostly unknown. Here we systematically decode a set of brain enhancers active in the zona limitans intrathalamica (zli), a signaling center essential for vertebrate forebrain development via the secreted morphogen Sonic hedgehog (Shh). We apply a de novo motif analysis tool to identify six position-independent sequence motifs together with their cognate transcription factors that are essential for zli enhancer activity and Shh expression in the mouse embryo. Using knowledge of this regulatory lexicon, we discover new Shh zli enhancers in mice and a functionally equivalent element in hemichordates, indicating an ancient origin of the Shh zli regulatory network that predates the chordate phylum. These findings support a strategy for delineating functionally conserved enhancers in the absence of overt sequence homologies and over extensive evolutionary distances.

Favorable genomic environments for cis-regulatory evolution: A novel theoretical framework

Cis-regulatory changes are arguably the primary evolutionary source of animal morphological diversity. With the recent explosion of genome-wide comparisons of the cis-regulatory content in different animal species is now possible to infer general principles underlying enhancer evolution. However, these studies have also revealed numerous discrepancies and paradoxes, suggesting that the mechanistic causes and modes of cis-regulatory evolution are still not well understood and are probably much more complex than generally appreciated. Here, we argue that the mutational mechanisms and genomic regions generating new regulatory activities must comply with the constraints imposed by the molecular properties of cis-regulatory elements (CREs) and the organizational features of long-range chromatin interactions. Accordingly, we propose a new integrative evolutionary framework for cis-regulatory evolution based on two major premises for the origin of novel enhancer activity: (i) an accessible chromatin environment and (ii) compatibility with the 3D structure and interactions of pre-existing CREs. Mechanisms and DNA sequences not fulfilling these premises, will be less likely to have a measurable impact on gene expression and as such, will have a minor contribution to the evolution of gene regulation. Finally, we discuss current comparative cis-regulatory data under the light of this new evolutionary model, and propose that the two most prominent mechanisms for the evolution of cis-regulatory changes are the overprinting of ancestral CREs and the exaptation of transposable elements.

Gene regulation in amphioxus: An insight from transgenic studies in amphioxus and vertebrates

Cephalochordates, commonly known as amphioxus or lancelets, are the most basal subphylum of chordates. Cephalochordates are thus key to understanding the origin of vertebrates and molecular mechanisms underlying vertebrate evolution. The evolution of developmental control mechanisms during invertebrate-to-vertebrate transition involved not only gene duplication events, but also specific changes in spatial and temporal expression of many genes. To get insight into the spatiotemporal regulation of gene expression during invertebrate-to-vertebrate transition, functional studies of amphioxus gene regulatory elements are highly warranted. Here, we review transgenic studies performed in amphioxus and vertebrates using promoters and enhancers derived from the genome of Branchiostoma floridae. We describe the current methods of transgenesis in amphioxus, provide evidence of Tol2 transposon-generated transgenic embryos of Branchiostoma lanceolatum and discuss possible future directions. We envision that comparative transgenic analysis of gene regulatory sequences in the context of amphioxus and vertebrate embryos will likely provide an important mechanistic insight into the evolution of vertebrate body plan.

Evolutionary genomics and adaptive evolution of the Hedgehog gene family (Shh, Ihh and Dhh) in vertebrates

The Hedgehog (Hh) gene family codes for a class of secreted proteins composed of two active domains that act as signalling molecules during embryo development, namely for the development of the nervous and skeletal systems and the formation of the testis cord. While only one Hh gene is found typically in invertebrate genomes, most vertebrates species have three (Sonic hedgehog – Shh; Indian hedgehog – Ihh; and Desert hedgehog – Dhh), each with different expression patterns and functions, which likely helped promote the increasing complexity of vertebrates and their successful diversification. In this study, we used comparative genomic and adaptive evolutionary analyses to characterize the evolution of the Hh genes in vertebrates following the two major whole genome duplication (WGD) events. To overcome the lack of Hh-coding sequences on avian publicly available databases, we used an extensive dataset of 45 avian and three non-avian reptilian genomes to show that birds have all three Hh paralogs. We find suggestions that following the WGD events, vertebrate Hh paralogous genes evolved independently within similar linkage groups and under different evolutionary rates, especially within the catalytic domain. The structural regions around the ion-binding site were identified to be under positive selection in the signaling domain. These findings contrast with those observed in invertebrates, where different lineages that experienced gene duplication retained similar selective constraints in the Hh orthologs. Our results provide new insights on the evolutionary history of the Hh gene family, the functional roles of these paralogs in vertebrate species, and on the location of mutational hotspots.

Genome-wide analysis of functional and evolutionary features of tele-enhancers

We investigated sequence features of enhancers separated from their target gene by at least one intermediate gene/exon (named tele-enhancers in this study) and enhancers residing inside their target gene locus. In this study, we used whole genome enhancer maps and gene expression profiles to establish a large panel of tele-enhancers. By contrasting tele-enhancers to proximal enhancers targeting heart genes, we observed that heart tele-enhancers use unique regulatory mechanisms based on the cardiac transcription factors SRF, TEAD, and NKX-2.5, whereas proximal heart enhancers rely on GATA4 instead. A functional analysis showed that tele-enhancers preferentially regulate house-keeping genes and genes with a metabolic role during heart development. In addition, tele-enhancers are significantly more conserved than their proximal counterparts. Similar trends have been observed for non-heart tissues and cell types, suggesting that our findings represent general characteristics of tele-enhancers.

Ancient cis-regulatory constraints and the evolution of genome architecture

The order of genes along metazoan chromosomes has generally been thought to be largely random, with few implications for organismal function. However, two recent studies, reporting hundreds of pairs of genes that have remained linked in diverse metazoan species over hundreds of millions of years of evolution, suggest widespread functional implications for gene order. These associations appear to largely reflect cis-regulatory constraints, with either (i) multiple genes sharing transcriptional regulatory elements, or (ii) regulatory elements for a developmental gene being found within a neighboring 'bystander' gene (known as a genomic regulatory block). We discuss implications, questions raised, and new research directions arising from these studies, as well as evidence for similar phenomena in other eukaryotic groups.

Differences in enhancer activity in mouse and zebrafish reporter assays are often associated with changes in gene expression

Background

Phenotypic evolution in animals is thought to be driven in large part by differences in gene expression patterns, which can result from sequence changes in cis-regulatory elements (cis-changes) or from changes in the expression pattern or function of transcription factors (trans-changes). While isolated examples of trans-changes have been identified, the scale of their overall contribution to regulatory and phenotypic evolution remains unclear.

Results

Here, we attempt to examine the prevalence of trans-effects and their potential impact on gene expression patterns in vertebrate evolution by comparing the function of identical human tissue-specific enhancer sequences in two highly divergent vertebrate model systems, mouse and zebrafish. Among 47 human conserved non-coding elements (CNEs) tested in transgenic mouse embryos and in stable zebrafish lines, at least one species-specific expression domain was observed in the majority (83%) of cases, and 36% presented dramatically different expression patterns between the two species. Although some of these discrepancies may be due to the use of different transgenesis systems in mouse and zebrafish, in some instances we found an association between differences in enhancer activity and changes in the endogenous gene expression patterns between mouse and zebrafish, suggesting a potential role for trans-changes in the evolution of gene expression.

Conclusions

In total, our results: (i) serve as a cautionary tale for studies investigating the role of human enhancers in different model organisms, and (ii) suggest that changes in the trans environment may play a significant role in the evolution of gene expression in vertebrates.

Extensive conservation of ancient microsynteny across metazoans due to cis-regulatory constraints

The order of genes in eukaryotic genomes has generally been assumed to be neutral, since gene order is largely scrambled over evolutionary time. Only a handful of exceptional examples are known, typically involving deeply conserved clusters of tandemly duplicated genes (e.g., Hox genes and histones). Here we report the first systematic survey of microsynteny conservation across metazoans, utilizing 17 genome sequences. We identified nearly 600 pairs of unrelated genes that have remained tightly physically linked in diverse lineages across over 600 million years of evolution. Integrating sequence conservation, gene expression data, gene function, epigenetic marks, and other genomic features, we provide extensive evidence that many conserved ancient linkages involve (1) the coordinated transcription of neighboring genes, or (2) genomic regulatory blocks (GRBs) in which transcriptional enhancers controlling developmental genes are contained within nearby bystander genes. In addition, we generated ChIP-seq data for key histone modifications in zebrafish embryos, which provided further evidence of putative GRBs in embryonic development. Finally, using chromosome conformation capture (3C) assays and stable transgenic experiments, we demonstrate that enhancers within bystander genes drive the expression of genes such as Otx and Islet, critical regulators of central nervous system development across bilaterians. These results suggest that ancient genomic functional associations are far more common than previously thought—involving ∼12% of the ancestral bilaterian genome—and that cis-regulatory constraints are crucial in determining metazoan genome architecture.