The little skate genome and the evolutionary emergence of wing-like fins.

Skates are cartilaginous fish whose body plan features enlarged wing-like pectoral fins, enabling them to thrive in benthic environments1,2. However, the molecular underpinnings of this unique trait remain unclear. Here we investigate the origin of this phenotypic innovation by developing the little skate Leucoraja erinacea as a genomically enabled model. Analysis of a high-quality chromosome-scale genome sequence for the little skate shows that it preserves many ancestral jawed vertebrate features compared with other sequenced genomes, including numerous ancient microchromosomes. Combining genome comparisons with extensive regulatory datasets in developing fins-including gene expression, chromatin occupancy and three-dimensional conformation-we find skate-specific genomic rearrangements that alter the three-dimensional regulatory landscape of genes that are involved in the planar cell polarity pathway. Functional inhibition of planar cell polarity signalling resulted in a reduction in anterior fin size, confirming that this pathway is a major contributor to batoid fin morphology. We also identified a fin-specific enhancer that interacts with several hoxa genes, consistent with the redeployment of hox gene expression in anterior pectoral fins, and confirmed its potential to activate transcription in the anterior fin using zebrafish reporter assays. Our findings underscore the central role of genome reorganization and regulatory variation in the evolution of phenotypes, shedding light on the molecular origin of an enigmatic trait.

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.

The emergence of the brain non-CpG methylation system in vertebrates.

Mammalian brains feature exceptionally high levels of non-CpG DNA methylation alongside the canonical form of CpG methylation. Non-CpG methylation plays a critical regulatory role in cognitive function, which is mediated by the binding of MeCP2, the transcriptional regulator that when mutated causes Rett syndrome. However, it is unclear whether the non-CpG neural methylation system is restricted to mammalian species with complex cognitive abilities or has deeper evolutionary origins. To test this, we investigated brain DNA methylation across 12 distantly related animal lineages, revealing that non-CpG methylation is restricted to vertebrates. We discovered that in vertebrates, non-CpG methylation is enriched within a highly conserved set of developmental genes transcriptionally repressed in adult brains, indicating that it demarcates a deeply conserved regulatory program. We also found that the writer of non-CpG methylation, DNMT3A, and the reader, MeCP2, originated at the onset of vertebrates as a result of the ancestral vertebrate whole-genome duplication. Together, we demonstrate how this novel layer of epigenetic information assembled at the root of vertebrates and gained new regulatory roles independent of the ancestral form of the canonical CpG methylation. This suggests that the emergence of non-CpG methylation may have fostered the evolution of sophisticated cognitive abilities found in the vertebrate lineage.

Amphioxus functional genomics and the origins of vertebrate gene regulation.

Vertebrates have greatly elaborated the basic chordate body plan and evolved highly distinctive genomes that have been sculpted by two whole-genome duplications. Here we sequence the genome of the Mediterranean amphioxus (Branchiostoma lanceolatum) and characterize DNA methylation, chromatin accessibility, histone modifications and transcriptomes across multiple developmental stages and adult tissues to investigate the evolution of the regulation of the chordate genome. Comparisons with vertebrates identify an intermediate stage in the evolution of differentially methylated enhancers, and a high conservation of gene expression and its cis-regulatory logic between amphioxus and vertebrates that occurs maximally at an earlier mid-embryonic phylotypic period. We analyse regulatory evolution after whole-genome duplications, and find that-in vertebrates-over 80% of broadly expressed gene families with multiple paralogues derived from whole-genome duplications have members that restricted their ancestral expression, and underwent specialization rather than subfunctionalization. Counter-intuitively, paralogues that restricted their expression increased the complexity of their regulatory landscapes. These data pave the way for a better understanding of the regulatory principles that underlie key vertebrate innovations.

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 (Oryzias latipes) ZRS and enhancer assays identify the existence of ZRS shadow enhancers in both teleost and human genomes. Deletion of both ZRS and shadow ZRS abolishes shh expression and completely truncates 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.

Active DNA demethylation at enhancers during the vertebrate phylotypic period.

The vertebrate body plan and organs are shaped during a conserved embryonic phase called the phylotypic stage. However, the mechanisms that guide the epigenome through this transition and their evolutionary conservation remain elusive. Here we report widespread DNA demethylation of enhancers during the phylotypic period in zebrafish, Xenopus tropicalis and mouse. These enhancers are linked to developmental genes that display coordinated transcriptional and epigenomic changes in the diverse vertebrates during embryogenesis. Binding of Tet proteins to (hydroxy)methylated DNA and enrichment of 5-hydroxymethylcytosine in these regions implicated active DNA demethylation in this process. Furthermore, loss of function of Tet1, Tet2 and Tet3 in zebrafish reduced chromatin accessibility and increased methylation levels specifically at these enhancers, indicative of DNA methylation being an upstream regulator of phylotypic enhancer function. Overall, our study highlights a regulatory module associated with the most conserved phase of vertebrate embryogenesis and suggests an ancient developmental role for Tet dioxygenases.

A single three-dimensional chromatin compartment in amphioxus indicates a stepwise evolution of vertebrate Hox bimodal regulation.

The HoxA and HoxD gene clusters of jawed vertebrates are organized into bipartite three-dimensional chromatin structures that separate long-range regulatory inputs coming from the anterior and posterior Hox-neighboring regions. This architecture is instrumental in allowing vertebrate Hox genes to pattern disparate parts of the body, including limbs. Almost nothing is known about how these three-dimensional topologies originated. Here we perform extensive 4C-seq profiling of the Hox cluster in embryos of amphioxus, an invertebrate chordate. We find that, in contrast to the architecture in vertebrates, the amphioxus Hox cluster is organized into a single chromatin interaction domain that includes long-range contacts mostly from the anterior side, bringing distant cis-regulatory elements into contact with Hox genes. We infer that the vertebrate Hox bipartite regulatory system is an evolutionary novelty generated by combining ancient long-range regulatory contacts from DNA in the anterior Hox neighborhood with new regulatory inputs from the posterior side.

Genome-wide epigenetic cross-talk between DNA methylation and H3K27me3 in zebrafish embryos.

DNA methylation and histone modifications are epigenetic marks implicated in the complex regulation of vertebrate embryogenesis. The cross-talk between DNA methylation and Polycomb-dependent H3K27me3 histone mark has been reported in a number of organisms [1], [2], [3], [4], [5], [6], [7] and both marks are known to be required for proper developmental progression. Here we provide genome-wide DNA methylation (MethylCap-seq) and H3K27me3 (ChIP-seq) maps for three stages (dome, 24 hpf and 48 hpf) of zebrafish (Danio rerio) embryogenesis, as well as all analytical and methodological details associated with the generation of this dataset. We observe a strong antagonism between the two epigenetic marks present in CpG islands and their compatibility throughout the bulk of the genome, as previously reported in mammalian ESC lines (Brinkman et al., 2012). Next generation sequencing data linked to this project have been deposited in the Gene Expression Omnibus (GEO) database under accession numbers GSE35050 and GSE70847.

Evolutionary comparison reveals that diverging CTCF sites are signatures of ancestral topological associating domains borders.

Increasing evidence in the last years indicates that the vast amount of regulatory information contained in mammalian genomes is organized in precise 3D chromatin structures. However, the impact of this spatial chromatin organization on gene expression and its degree of evolutionary conservation is still poorly understood. The Six homeobox genes are essential developmental regulators organized in gene clusters conserved during evolution. Here, we reveal that the Six clusters share a deeply evolutionarily conserved 3D chromatin organization that predates the Cambrian explosion. This chromatin architecture generates two largely independent regulatory landscapes (RLs) contained in two adjacent topological associating domains (TADs). By disrupting the conserved TAD border in one of the zebrafish Six clusters, we demonstrate that this border is critical for preventing competition between promoters and enhancers located in separated RLs, thereby generating different expression patterns in genes located in close genomic proximity. Moreover, evolutionary comparison of Six-associated TAD borders reveals the presence of CCCTC-binding factor (CTCF) sites with diverging orientations in all studied deuterostomes. Genome-wide examination of mammalian HiC data reveals that this conserved CTCF configuration is a general signature of TAD borders, underscoring that common organizational principles underlie TAD compartmentalization in deuterostome evolution.

Deep conservation of wrist and digit enhancers in fish.

There is no obvious morphological counterpart of the autopod (wrist/ankle and digits) in living fishes. Comparative molecular data may provide insight into understanding both the homology of elements and the evolutionary developmental mechanisms behind the fin to limb transition. In mouse limbs the autopod is built by a "late" phase of Hoxd and Hoxa gene expression, orchestrated by a set of enhancers located at the 5' end of each cluster. Despite a detailed mechanistic understanding of mouse limb development, interpretation of Hox expression patterns and their regulation in fish has spawned multiple hypotheses as to the origin and function of "autopod" enhancers throughout evolution. Using phylogenetic footprinting, epigenetic profiling, and transgenic reporters, we have identified and functionally characterized hoxD and hoxA enhancers in the genomes of zebrafish and the spotted gar, Lepisosteus oculatus, a fish lacking the whole genome duplication of teleosts. Gar and zebrafish "autopod" enhancers drive expression in the distal portion of developing zebrafish pectoral fins, and respond to the same functional cues as their murine orthologs. Moreover, gar enhancers drive reporter gene expression in both the wrist and digits of mouse embryos in patterns that are nearly indistinguishable from their murine counterparts. These functional genomic data support the hypothesis that the distal radials of bony fish are homologous to the wrist and/or digits of tetrapods.

A polymorphic enhancer near GREM1 influences bowel cancer risk through differential CDX2 and TCF7L2 binding.

A rare germline duplication upstream of the bone morphogenetic protein antagonist GREM1 causes a Mendelian-dominant predisposition to colorectal cancer (CRC). The underlying disease mechanism is strong, ectopic GREM1 overexpression in the intestinal epithelium. Here, we confirm that a common GREM1 polymorphism, rs16969681, is also associated with CRC susceptibility, conferring ∼20% differential risk in the general population. We hypothesized the underlying cause to be moderate differences in GREM1 expression. We showed that rs16969681 lies in a region of active chromatin with allele- and tissue-specific enhancer activity. The CRC high-risk allele was associated with stronger gene expression, and higher Grem1 mRNA levels increased the intestinal tumor burden in Apc(Min) mice. The intestine-specific transcription factor CDX2 and Wnt effector TCF7L2 bound near rs16969681, with significantly higher affinity for the risk allele, and CDX2 overexpression in CDX2/GREM1-negative cells caused re-expression of GREM1. rs16969681 influences CRC risk through effects on Wnt-driven GREM1 expression in colorectal tumors.

The genome sequencing of an albino Western lowland gorilla reveals inbreeding in the wild.

BACKGROUND: The only known albino gorilla, named Snowflake, was a male wild born individual from Equatorial Guinea who lived at the Barcelona Zoo for almost 40 years. He was diagnosed with non-syndromic oculocutaneous albinism, i.e. white hair, light eyes, pink skin, photophobia and reduced visual acuity. Despite previous efforts to explain the genetic cause, this is still unknown. Here, we study the genetic cause of his albinism and making use of whole genome sequencing data we find a higher inbreeding coefficient compared to other gorillas. RESULTS: We successfully identified the causal genetic variant for Snowflake's albinism, a non-synonymous single nucleotide variant located in a transmembrane region of SLC45A2. This transporter is known to be involved in oculocutaneous albinism type 4 (OCA4) in humans. We provide experimental evidence that shows that this amino acid replacement alters the membrane spanning capability of this transmembrane region. Finally, we provide a comprehensive study of genome-wide patterns of autozygogosity revealing that Snowflake's parents were related, being this the first report of inbreeding in a wild born Western lowland gorilla. CONCLUSIONS: In this study we demonstrate how the use of whole genome sequencing can be extended to link genotype and phenotype in non-model organisms and it can be a powerful tool in conservation genetics (e.g., inbreeding and genetic diversity) with the expected decrease in sequencing cost.

The developmental epigenomics toolbox: ChIP-seq and MethylCap-seq profiling of early zebrafish embryos.

Genome-wide profiling of DNA methylation and histone modifications answered many questions as to how the genes are regulated on a global scale and what their epigenetic makeup is. Yet, little is known about the function of these marks during early vertebrate embryogenesis. Here we provide detailed protocols for ChIP-seq and MethylCap-seq procedures applied to zebrafish (Danio rerio) embryonic material at four developmental stages. As a proof of principle, we have profiled on a global scale a number of post-translational histone modifications including H3K4me1, H3K4me3 and H3K27ac. We demonstrate that these marks are dynamic during early development and that such developmental transitions can be detected by ChIP-seq. In addition, we applied MethylCap-seq to show that developmentally-regulated DNA methylation remodeling can be detected by such a procedure. Our MethylCap-seq data concur with previous DNA methylation studies of early zebrafish development rendering this method highly suitable for the global assessment of DNA methylation in early vertebrate embryos.

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.

Dynamics of enhancer chromatin signatures mark the transition from pluripotency to cell specification during embryogenesis.

The generation of distinctive cell types that form different tissues and organs requires precise, temporal and spatial control of gene expression. This depends on specific cis-regulatory elements distributed in the noncoding DNA surrounding their target genes. Studies performed on mammalian embryonic stem cells and Drosophila embryos suggest that active enhancers form part of a defined chromatin landscape marked by histone H3 lysine 4 mono-methylation (H3K4me1) and histone H3 lysine 27 acetylation (H3K27ac). Nevertheless, little is known about the dynamics and the potential roles of these marks during vertebrate embryogenesis. Here, we provide genomic maps of H3K4me1/me3 and H3K27ac at four developmental time-points of zebrafish embryogenesis and analyze embryonic enhancer activity. We find that (1) changes in H3K27ac enrichment at enhancers accompany the shift from pluripotency to tissue-specific gene expression, (2) in early embryos, the peaks of H3K27ac enrichment are bound by pluripotent factors such as Nanog, and (3) the degree of evolutionary conservation is higher for enhancers that become marked by H3K27ac at the end of gastrulation, suggesting their implication in the establishment of the most conserved (phylotypic) transcriptome that is known to occur later at the pharyngula stage.

Genome-wide CTCF distribution in vertebrates defines equivalent sites that aid the identification of disease-associated genes.

Many genomic alterations associated with human diseases localize in noncoding regulatory elements located far from the promoters they regulate, making it challenging to link noncoding mutations or risk-associated variants with target genes. The range of action of a given set of enhancers is thought to be defined by insulator elements bound by the 11 zinc-finger nuclear factor CCCTC-binding protein (CTCF). Here we analyzed the genomic distribution of CTCF in various human, mouse and chicken cell types, demonstrating the existence of evolutionarily conserved CTCF-bound sites beyond mammals. These sites preferentially flank transcription factor-encoding genes, often associated with human diseases, and function as enhancer blockers in vivo, suggesting that they act as evolutionarily invariant gene boundaries. We then applied this concept to predict and functionally demonstrate that the polymorphic variants associated with multiple sclerosis located within the EVI5 gene impinge on the adjacent gene GFI1.

An evolutionarily conserved three-dimensional structure in the vertebrate Irx clusters facilitates enhancer sharing and coregulation.

Developmental gene clusters are paradigms for the study of gene regulation; however, the mechanisms that mediate phenomena such as coregulation and enhancer sharing remain largely elusive. Here we address this issue by analysing the vertebrate Irx clusters. We first present a deep enhancer screen of a 2-Mbp span covering the IrxA cluster. Using chromosome conformation capture, we show that enhancer sharing is widespread within the cluster, explaining its evolutionarily conserved organization. We also identify a three-dimensional architecture, probably formed through interactions with CCCTC-binding factor, which is present within both Irx clusters of mouse, Xenopus and zebrafish. This architecture brings the promoters of the first two genes together in the same chromatin landscape. We propose that this unique and evolutionarily conserved genomic architecture of the vertebrate Irx clusters is essential for the coregulation of the first two genes and simultaneously maintains the third gene in a partially independent regulatory landscape.

Genome-wide profiling of p63 DNA-binding sites identifies an element that regulates gene expression during limb development in the 7q21 SHFM1 locus.

Heterozygous mutations in p63 are associated with split hand/foot malformations (SHFM), orofacial clefting, and ectodermal abnormalities. Elucidation of the p63 gene network that includes target genes and regulatory elements may reveal new genes for other malformation disorders. We performed genome-wide DNA-binding profiling by chromatin immunoprecipitation (ChIP), followed by deep sequencing (ChIP-seq) in primary human keratinocytes, and identified potential target genes and regulatory elements controlled by p63. We show that p63 binds to an enhancer element in the SHFM1 locus on chromosome 7q and that this element controls expression of DLX6 and possibly DLX5, both of which are important for limb development. A unique micro-deletion including this enhancer element, but not the DLX5/DLX6 genes, was identified in a patient with SHFM. Our study strongly indicates disruption of a non-coding cis-regulatory element located more than 250 kb from the DLX5/DLX6 genes as a novel disease mechanism in SHFM1. These data provide a proof-of-concept that the catalogue of p63 binding sites identified in this study may be of relevance to the studies of SHFM and other congenital malformations that resemble the p63-associated phenotypes.

Multiple enhancers located in a 1-Mb region upstream of POU3F4 promote expression during inner ear development and may be required for hearing.

POU3F4 encodes a POU-domain transcription factor required for inner ear development. Defects in POU3F4 function are associated with X-linked deafness type 3 (DFN3). Multiple deletions affecting up to ~900-kb upstream of POU3F4 are found in DFN3 patients, suggesting the presence of essential POU3F4 enhancers in this region. Recently, an inner ear enhancer was reported that is absent in most DFN3 patients with upstream deletions. However, two indications suggest that additional enhancers in the POU3F4 upstream region are required for POU3F4 function during inner ear development. First, there is at least one DFN3 deletion that does not eliminate the reported enhancer. Second, the expression pattern driven by this enhancer does not fully recapitulate Pou3f4 expression in the inner ear. Here, we screened a 1-Mb region upstream of the POU3F4 gene for additional cis-regulatory elements and searched for novel DFN3 mutations in the identified POU3F4 enhancers. We found several novel enhancers for otic vesicle expression. Some of these also drive expression in kidney, pancreas and brain, tissues that are known to express Pou3f4. In addition, we report a new and smallest deletion identified so far in a DFN3 family which eliminates 3.9 kb, comprising almost exclusively the previous reported inner ear enhancer. We suggest that multiple enhancers control the expression of Pou3f4 in the inner ear and these may contribute to the phenotype observed in DFN3 patients. In addition, the novel deletion demonstrates that the previous reported enhancer, although not sufficient, is essential for POU3F4 function during inner ear development.