Parallel evolution of amphioxus and vertebrate small-scale gene duplications.

BACKGROUND: Amphioxus are non-vertebrate chordates characterized by a slow morphological and molecular evolution. They share the basic chordate body-plan and genome organization with vertebrates but lack their 2R whole-genome duplications and their developmental complexity. For these reasons, amphioxus are frequently used as an outgroup to study vertebrate genome evolution and Evo-Devo. Aside from whole-genome duplications, genes continuously duplicate on a smaller scale. Small-scale duplicated genes can be found in both amphioxus and vertebrate genomes, while only the vertebrate genomes have duplicated genes product of their 2R whole-genome duplications. Here, we explore the history of small-scale gene duplications in the amphioxus lineage and compare it to small- and large-scale gene duplication history in vertebrates. RESULTS: We present a study of the European amphioxus (Branchiostoma lanceolatum) gene duplications thanks to a new, high-quality genome reference. We find that, despite its overall slow molecular evolution, the amphioxus lineage has had a history of small-scale duplications similar to the one observed in vertebrates. We find parallel gene duplication profiles between amphioxus and vertebrates and conserved functional constraints in gene duplication. Moreover, amphioxus gene duplicates show levels of expression and patterns of functional specialization similar to the ones observed in vertebrate duplicated genes. We also find strong conservation of gene synteny between two distant amphioxus species, B. lanceolatum and B. floridae, with two major chromosomal rearrangements. CONCLUSIONS: In contrast to their slower molecular and morphological evolution, amphioxus' small-scale gene duplication history resembles that of the vertebrate lineage both in quantitative and in functional terms.

Duplications disrupt chromatin architecture and rewire GPR101-enhancer communication in X-linked acrogigantism.

X-linked acrogigantism (X-LAG) is the most severe form of pituitary gigantism and is characterized by aggressive growth hormone (GH)-secreting pituitary tumors that occur in early childhood. X-LAG is associated with chromosome Xq26.3 duplications (the X-LAG locus typically includes VGLL1, CD40LG, ARHGEF6, RBMX, and GPR101) that lead to massive pituitary tumoral expression of GPR101, a novel regulator of GH secretion. The mechanism by which the duplications lead to marked pituitary misexpression of GPR101 alone was previously unclear. Using Hi-C and 4C-seq, we characterized the normal chromatin structure at the X-LAG locus. We showed that GPR101 is located within a topologically associating domain (TAD) delineated by a tissue-invariant border that separates it from centromeric genes and regulatory sequences. Next, using 4C-seq with GPR101, RBMX, and VGLL1 viewpoints, we showed that the duplications in multiple X-LAG-affected individuals led to ectopic interactions that crossed the invariant TAD border, indicating the existence of a similar and consistent mechanism of neo-TAD formation in X-LAG. We then identified several pituitary active cis-regulatory elements (CREs) within the neo-TAD and demonstrated in vitro that one of them significantly enhanced reporter gene expression. At the same time, we showed that the GPR101 promoter permits the incorporation of new regulatory information. Our results indicate that X-LAG is a TADopathy of the endocrine system in which Xq26.3 duplications disrupt the local chromatin architecture forming a neo-TAD. Rewiring GPR101-enhancer interaction within the new regulatory unit is likely to cause the high levels of aberrant expression of GPR101 in pituitary tumors caused by X-LAG.

CTCF knockout in zebrafish induces alterations in regulatory landscapes and developmental gene expression.

Coordinated chromatin interactions between enhancers and promoters are critical for gene regulation. The architectural protein CTCF mediates chromatin looping and is enriched at the boundaries of topologically associating domains (TADs), which are sub-megabase chromatin structures. In vitro CTCF depletion leads to a loss of TADs but has only limited effects over gene expression, challenging the concept that CTCF-mediated chromatin structures are a fundamental requirement for gene regulation. However, how CTCF and a perturbed chromatin structure impacts gene expression during development remains poorly understood. Here we link the loss of CTCF and gene regulation during patterning and organogenesis in a ctcf knockout zebrafish model. CTCF absence leads to loss of chromatin structure and affects the expression of thousands of genes, including many developmental regulators. Our results demonstrate the essential role of CTCF in providing the structural context for enhancer-promoter interactions, thus regulating developmental genes.

Regulatory landscape fusion in rhabdomyosarcoma through interactions between the PAX3 promoter and FOXO1 regulatory elements.

BACKGROUND: The organisation of vertebrate genomes into topologically associating domains (TADs) is believed to facilitate the regulation of the genes located within them. A remaining question is whether TAD organisation is achieved through the interactions of the regulatory elements within them or if these interactions are favoured by the pre-existence of TADs. If the latter is true, the fusion of two independent TADs should result in the rewiring of the transcriptional landscape and the generation of ectopic contacts. RESULTS: We show that interactions within the PAX3 and FOXO1 domains are restricted to their respective TADs in normal conditions, while in a patient-derived alveolar rhabdomyosarcoma cell line, harbouring the diagnostic t(2;13)(q35;q14) translocation that brings together the PAX3 and FOXO1 genes, the PAX3 promoter interacts ectopically with FOXO1 sequences. Using a combination of 4C-seq datasets, we have modelled the three-dimensional organisation of the fused landscape in alveolar rhabdomyosarcoma. CONCLUSIONS: The chromosomal translocation that leads to alveolar rhabdomyosarcoma development generates a novel TAD that is likely to favour ectopic PAX3:FOXO1 oncogene activation in non-PAX3 territories. Rhabdomyosarcomas may therefore arise from cells which do not normally express PAX3. The borders of this novel TAD correspond to the original 5'- and 3'- borders of the PAX3 and FOXO1 TADs, respectively, suggesting that TAD organisation precedes the formation of regulatory long-range interactions. Our results demonstrate that, upon translocation, novel regulatory landscapes are formed allowing new intra-TAD interactions between the original loci involved.

Analysis of opo cis-regulatory landscape uncovers Vsx2 requirement in early eye morphogenesis.

The self-organized morphogenesis of the vertebrate optic cup entails coupling the activation of the retinal gene regulatory network to the constriction-driven infolding of the retinal epithelium. Yet the genetic mechanisms underlying this coordination remain largely unexplored. Through phylogenetic footprinting and transgenesis in zebrafish, here we examine the cis-regulatory landscape of opo, an endocytosis regulator essential for eye morphogenesis. Among the different conserved enhancers identified, we isolate a single retina-specific element (H6_10137) and show that its activity depends on binding sites for the retinal determinant Vsx2. Gain- and loss-of-function experiments and ChIP analyses reveal that Vsx2 regulates opo expression through direct binding to this retinal enhancer. Furthermore, we show that vsx2 knockdown impairs the primary optic cup folding. These data support a model by which vsx2, operating through the effector gene opo, acts as a central transcriptional node that coordinates neural retina patterning and optic cup invagination in zebrafish.

Comparative epigenomics in distantly related teleost species identifies conserved cis-regulatory nodes active during the vertebrate phylotypic period.

The complex relationship between ontogeny and phylogeny has been the subject of attention and controversy since von Baer's formulations in the 19th century. The classic concept that embryogenesis progresses from clade general features to species-specific characters has often been revisited. It has become accepted that embryos from a clade show maximum morphological similarity at the so-called phylotypic period (i.e., during mid-embryogenesis). According to the hourglass model, body plan conservation would depend on constrained molecular mechanisms operating at this period. More recently, comparative transcriptomic analyses have provided conclusive evidence that such molecular constraints exist. Examining cis-regulatory architecture during the phylotypic period is essential to understand the evolutionary source of body plan stability. Here we compare transcriptomes and key epigenetic marks (H3K4me3 and H3K27ac) from medaka (Oryzias latipes) and zebrafish (Danio rerio), two distantly related teleosts separated by an evolutionary distance of 115-200 Myr. We show that comparison of transcriptome profiles correlates with anatomical similarities and heterochronies observed at the phylotypic stage. Through comparative epigenomics, we uncover a pool of conserved regulatory regions (≈700), which are active during the vertebrate phylotypic period in both species. Moreover, we show that their neighboring genes encode mainly transcription factors with fundamental roles in tissue specification. We postulate that these regulatory regions, active in both teleost genomes, represent key constrained nodes of the gene networks that sustain the vertebrate body plan.

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

The origin and evolution of the complex regulatory landscapes of some vertebrate developmental genes, often spanning hundreds of Kbp and including neighboring genes, remain poorly understood. The Sonic Hedgehog (Shh) genomic regulatory block (GRB) is one of the best functionally characterized examples, with several discrete enhancers reported within its introns, vast upstream gene-free region and neighboring genes (Lmbr1 and Rnf32). To investigate the origin and evolution of this GRB, we sequenced and characterized the Hedgehog (Hh) loci from three invertebrate chordate amphioxus species, which share several early expression domains with Shh. Using phylogenetic footprinting within and between chordate lineages, and reporter assays in zebrafish probing >30 Kbp of amphioxus Hh, we report large sequence and functional divergence between both groups. In addition, we show that the linkage of Shh to Lmbr1 and Rnf32, necessary for the unique gnatostomate-specific Shh limb expression, is a vertebrate novelty occurred between the two whole-genome duplications.

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.

Dioxin receptor and SLUG transcription factors regulate the insulator activity of B1 SINE retrotransposons via an RNA polymerase switch.

Complex genomes utilize insulators and boundary elements to help define spatial and temporal gene expression patterns. We report that a genome-wide B1 SINE (Short Interspersed Nuclear Element) retrotransposon (B1-X35S) has potent intrinsic insulator activity in cultured cells and live animals. This insulation is mediated by binding of the transcription factors dioxin receptor (AHR) and SLUG (SNAI2) to consensus elements present in the SINE. Transcription of B1-X35S is required for insulation. While basal insulator activity is maintained by RNA polymerase (Pol) III transcription, AHR-induced insulation involves release of Pol III and engagement of Pol II transcription on the same strand. B1-X35S insulation is also associated with enrichment of heterochromatin marks H3K9me3 and H3K27me3 downstream of B1-X35S, an effect that varies with cell type. B1-X35S binds parylated CTCF and, consistent with a chromatin barrier activity, its positioning between two adjacent genes correlates with their differential expression in mouse tissues. Hence, B1 SINE retrotransposons represent genome-wide insulators activated by transcription factors that respond to developmental, oncogenic, or toxicological stimuli.