Functions and mechanisms of epigenetic inheritance in animals.

The idea that epigenetic determinants such as DNA methylation, histone modifications or RNA can be passed to the next generation through meiotic products (gametes) is long standing. Such meiotic epigenetic inheritance (MEI) is fairly common in yeast, plants and nematodes, but its extent in mammals has been much debated. Advances in genomics techniques are now driving the profiling of germline and zygotic epigenomes, thereby improving our understanding of MEI in diverse species. Whereas the role of DNA methylation in MEI remains unclear, insights from genome-wide studies suggest that a previously underappreciated fraction of mammalian genomes bypass epigenetic reprogramming during development. Notably, intergenerational inheritance of histone modifications, tRNA fragments and microRNAs can affect gene regulation in the offspring. It is important to note that MEI in mammals rarely constitutes transgenerational epigenetic inheritance (TEI), which spans multiple generations. In this Review, we discuss the examples of MEI in mammals, including mammalian epigenome reprogramming, and the molecular mechanisms of MEI in vertebrates in general. We also discuss the implications of the inheritance of histone modifications and small RNA for embryogenesis in metazoans, with a particular focus on insights gained from genome-wide studies.

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.

Evolutionary conservation of embryonic DNA methylome remodelling in distantly related teleost species.

Methylation of cytosines in the CG context (mCG) is the most abundant DNA modification in vertebrates that plays crucial roles in cellular differentiation and identity. After fertilization, DNA methylation patterns inherited from parental gametes are remodelled into a state compatible with embryogenesis. In mammals, this is achieved through the global erasure and re-establishment of DNA methylation patterns. However, in non-mammalian vertebrates like zebrafish, no global erasure has been observed. To investigate the evolutionary conservation and divergence of DNA methylation remodelling in teleosts, we generated base resolution DNA methylome datasets of developing medaka and medaka-zebrafish hybrid embryos. In contrast to previous reports, we show that medaka display comparable DNA methylome dynamics to zebrafish with high gametic mCG levels (sperm: ∼90%; egg: ∼75%), and adoption of a paternal-like methylome during early embryogenesis, with no signs of prior DNA methylation erasure. We also demonstrate that non-canonical DNA methylation (mCH) reprogramming at TGCT tandem repeats is a conserved feature of teleost embryogenesis. Lastly, we find remarkable evolutionary conservation of DNA methylation remodelling patterns in medaka-zebrafish hybrids, indicative of compatible DNA methylation maintenance machinery in far-related teleost species. Overall, these results suggest strong evolutionary conservation of DNA methylation remodelling pathways in teleosts, which is distinct from the global DNA methylome erasure and reestablishment observed in mammals.

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.

Genome evolution in the allotetraploid frog Xenopus laevis.

To explore the origins and consequences of tetraploidy in the African clawed frog, we sequenced the Xenopus laevis genome and compared it to the related diploid X. tropicalis genome. We characterize the allotetraploid origin of X. laevis by partitioning its genome into two homoeologous subgenomes, marked by distinct families of 'fossil' transposable elements. On the basis of the activity of these elements and the age of hundreds of unitary pseudogenes, we estimate that the two diploid progenitor species diverged around 34 million years ago (Ma) and combined to form an allotetraploid around 17-18 Ma. More than 56% of all genes were retained in two homoeologous copies. Protein function, gene expression, and the amount of conserved flanking sequence all correlate with retention rates. The subgenomes have evolved asymmetrically, with one chromosome set more often preserving the ancestral state and the other experiencing more gene loss, deletion, rearrangement, and reduced gene expression.

Developmental remodelling of non-CG methylation at satellite DNA repeats.

In vertebrates, DNA methylation predominantly occurs at CG dinucleotides however, widespread non-CG methylation (mCH) has been reported in mammalian embryonic stem cells and in the brain. In mammals, mCH is found at CAC trinucleotides in the nervous system, where it is associated with transcriptional repression, and at CAG trinucleotides in embryonic stem cells, where it positively correlates with transcription. Moreover, CAC methylation appears to be a conserved feature of adult vertebrate brains. Unlike any of those methylation signatures, here we describe a novel form of mCH that occurs in the TGCT context within zebrafish mosaic satellite repeats. TGCT methylation is inherited from both male and female gametes, remodelled during mid-blastula transition, and re-established during gastrulation in all embryonic layers. Moreover, we identify DNA methyltransferase 3ba (Dnmt3ba) as the primary enzyme responsible for the deposition of this mCH mark. Finally, we observe that TGCT-methylated repeats are specifically associated with H3K9me3-marked heterochromatin suggestive of a functional interplay between these two gene-regulatory marks. Altogether, this work provides insight into a novel form of vertebrate mCH and highlights the substrate diversity of vertebrate DNA methyltransferases.

Chromosome-length genome assembly and structural variations of the primal Basenji dog (Canis lupus familiaris) genome.

BACKGROUND: Basenjis are considered an ancient dog breed of central African origins that still live and hunt with tribesmen in the African Congo. Nicknamed the barkless dog, Basenjis possess unique phylogeny, geographical origins and traits, making their genome structure of great interest. The increasing number of available canid reference genomes allows us to examine the impact the choice of reference genome makes with regard to reference genome quality and breed relatedness. RESULTS: Here, we report two high quality de novo Basenji genome assemblies: a female, China (CanFam_Bas), and a male, Wags. We conduct pairwise comparisons and report structural variations between assembled genomes of three dog breeds: Basenji (CanFam_Bas), Boxer (CanFam3.1) and German Shepherd Dog (GSD) (CanFam_GSD). CanFam_Bas is superior to CanFam3.1 in terms of genome contiguity and comparable overall to the high quality CanFam_GSD assembly. By aligning short read data from 58 representative dog breeds to three reference genomes, we demonstrate how the choice of reference genome significantly impacts both read mapping and variant detection. CONCLUSIONS: The growing number of high-quality canid reference genomes means the choice of reference genome is an increasingly critical decision in subsequent canid variant analyses. The basal position of the Basenji makes it suitable for variant analysis for targeted applications of specific dog breeds. However, we believe more comprehensive analyses across the entire family of canids is more suited to a pangenome approach. Collectively this work highlights the importance the choice of reference genome makes in all variation studies.

Sequence determinants, function, and evolution of CpG islands.

In vertebrates, cytosine-guanine (CpG) dinucleotides are predominantly methylated, with ∼80% of all CpG sites containing 5-methylcytosine (5mC), a repressive mark associated with long-term gene silencing. The exceptions to such a globally hypermethylated state are CpG-rich DNA sequences called CpG islands (CGIs), which are mostly hypomethylated relative to the bulk genome. CGIs overlap promoters from the earliest vertebrates to humans, indicating a concerted evolutionary drive compatible with CGI retention. CGIs are characterised by DNA sequence features that include DNA hypomethylation, elevated CpG and GC content and the presence of transcription factor binding sites. These sequence characteristics are congruous with the recruitment of transcription factors and chromatin modifying enzymes, and transcriptional activation in general. CGIs colocalize with sites of transcriptional initiation in hypermethylated vertebrate genomes, however, a growing body of evidence indicates that CGIs might exert their gene regulatory function in other genomic contexts. In this review, we discuss the diverse regulatory features of CGIs, their functional readout, and the evolutionary implications associated with CGI retention in vertebrates and possibly in invertebrates.

Rapid Recovery Gene Downregulation during Excess-Light Stress and Recovery in Arabidopsis.

Stress recovery may prove to be a promising approach to increase plant performance and, theoretically, mRNA instability may facilitate faster recovery. Transcriptome (RNA-seq, qPCR, sRNA-seq, and PARE) and methylome profiling during repeated excess-light stress and recovery was performed at intervals as short as 3 min. We demonstrate that 87% of the stress-upregulated mRNAs analyzed exhibit very rapid recovery. For instance, HSP101 abundance declined 2-fold every 5.1 min. We term this phenomenon rapid recovery gene downregulation (RRGD), whereby mRNA abundance rapidly decreases promoting transcriptome resetting. Decay constants (k) were modeled using two strategies, linear and nonlinear least squares regressions, with the latter accounting for both transcription and degradation. This revealed extremely short half-lives ranging from 2.7 to 60.0 min for 222 genes. Ribosome footprinting using degradome data demonstrated RRGD loci undergo cotranslational decay and identified changes in the ribosome stalling index during stress and recovery. However, small RNAs and 5'-3' RNA decay were not essential for recovery of the transcripts examined, nor were any of the six excess light-associated methylome changes. We observed recovery-specific gene expression networks upon return to favorable conditions and six transcriptional memory types. In summary, rapid transcriptome resetting is reported in the context of active recovery and cellular memory.

TET enzymes, DNA demethylation and pluripotency.

Ten-eleven translocation (TET) methylcytosine dioxygenases (TET1, TET2, TET3) actively cause demethylation of 5-methylcytosine (5mC) and produce and safeguard hypomethylation at key regulatory regions across the genome. This 5mC erasure is particularly important in pluripotent embryonic stem cells (ESCs) as they need to maintain self-renewal capabilities while retaining the potential to generate different cell types with diverse 5mC patterns. In this review, we discuss the multiple roles of TET proteins in mouse ESCs, and other vertebrate model systems, with a particular focus on TET functions in pluripotency, differentiation, and developmental DNA methylome reprogramming. Furthermore, we elaborate on the recently described non-catalytic roles of TET proteins in diverse biological contexts. Overall, TET proteins are multifunctional regulators that through both their catalytic and non-catalytic roles carry out myriad functions linked to early developmental processes.

Meis1 coordinates a network of genes implicated in eye development and microphthalmia.

Microphthalmos is a rare congenital anomaly characterized by reduced eye size and visual deficits of variable degree. Sporadic and hereditary microphthalmos have been associated with heterozygous mutations in genes fundamental for eye development. Yet, many cases are idiopathic or await the identification of molecular causes. Here we show that haploinsufficiency of Meis1, which encodes a transcription factor with evolutionarily conserved expression in the embryonic trunk, brain and sensory organs, including the eye, causes microphthalmic traits and visual impairment in adult mice. By combining analysis of Meis1 loss-of-function and conditional Meis1 functional rescue with ChIP-seq and RNA-seq approaches we show that, in contrast to its preferential association with Hox-Pbx BSs in the trunk, Meis1 binds to Hox/Pbx-independent sites during optic cup development. In the eye primordium, Meis1 coordinates, in a dose-dependent manner, retinal proliferation and differentiation by regulating genes responsible for human microphthalmia and components of the Notch signaling pathway. In addition, Meis1 is required for eye patterning by controlling a set of eye territory-specific transcription factors, so that in Meis1(-/-) embryos boundaries among the different eye territories are shifted or blurred. We propose that Meis1 is at the core of a genetic network implicated in eye patterning/microphthalmia, and represents an additional candidate for syndromic cases of these ocular malformations.

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.

Temporal uncoupling of the DNA methylome and transcriptional repression during embryogenesis.

DNA methylation is a tightly regulated epigenetic mark associated with transcriptional repression. Next-generation sequencing of purified methylated DNA obtained from early Xenopus tropicalis embryos demonstrates that this genome is heavily methylated during blastula and gastrula stages. Although DNA methylation is largely absent from transcriptional start sites marked with histone H3 lysine 4 trimethylation (H3K4me3), we find both promoters and gene bodies of active genes robustly methylated. In contrast, DNA methylation is absent in large H3K27me3 domains, indicating that these two repression pathways have different roles. Comparison with chromatin state maps of human ES cells reveals strong conservation of epigenetic makeup and gene regulation between the two systems. Strikingly, genes that are highly expressed in pluripotent cells and in Xenopus embryos but not in differentiated cells exhibit relatively high DNA methylation. Therefore, we tested the repressive potential of DNA methylation using transient and transgenic approaches and show that methylated promoters are robustly transcribed in blastula- and gastrula-stage embryos, but not in oocytes or late embryos. These findings have implications for reprogramming and the epigenetic regulation of pluripotency and differentiation and suggest a relatively open, pliable chromatin state in early embryos followed by reestablished methylation-dependent transcriptional repression during organogenesis and differentiation.

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 DNA methylome rearrangement during early lamprey embryogenesis.

DNA methylation (5mC) is a repressive gene regulatory mark widespread in vertebrate genomes, yet the developmental dynamics in which 5mC patterns are established vary across species. While mammals undergo two rounds of global 5mC erasure, teleosts, for example, exhibit localized maternal-to-paternal 5mC remodeling. Here, we studied 5mC dynamics during the embryonic development of sea lamprey, a jawless vertebrate which occupies a critical phylogenetic position as the sister group of the jawed vertebrates. We employed 5mC quantification in lamprey embryos and tissues, and discovered large-scale maternal-to-paternal epigenome remodeling that affects ~30% of the embryonic genome and is predominantly associated with partially methylated domains. We further demonstrate that sequences eliminated during programmed genome rearrangement (PGR), are hypermethylated in sperm prior to the onset of PGR. Our study thus unveils important insights into the evolutionary origins of vertebrate 5mC reprogramming, and how this process might participate in diverse developmental strategies.

Tumor Biomechanics Alters Metastatic Dissemination of Triple Negative Breast Cancer via Rewiring Fatty Acid Metabolism.

In recent decades, the role of tumor biomechanics on cancer cell behavior at the primary site has been increasingly appreciated. However, the effect of primary tumor biomechanics on the latter stages of the metastatic cascade, such as metastatic seeding of secondary sites and outgrowth remains underappreciated. This work sought to address this in the context of triple negative breast cancer (TNBC), a cancer type known to aggressively disseminate at all stages of disease progression. Using mechanically tuneable model systems, mimicking the range of stiffness's typically found within breast tumors, it is found that, contrary to expectations, cancer cells exposed to softer microenvironments are more able to colonize secondary tissues. It is shown that heightened cell survival is driven by enhanced metabolism of fatty acids within TNBC cells exposed to softer microenvironments. It is demonstrated that uncoupling cellular mechanosensing through integrin ß1 blocking antibody effectively causes stiff primed TNBC cells to behave like their soft counterparts, both in vitro and in vivo. This work is the first to show that softer tumor microenvironments may be contributing to changes in disease outcome by imprinting on TNBC cells a greater metabolic flexibility and conferring discrete cell survival advantages.

The Australasian dingo archetype: de novo chromosome-length genome assembly, DNA methylome, and cranial morphology.

BACKGROUND: One difficulty in testing the hypothesis that the Australasian dingo is a functional intermediate between wild wolves and domesticated breed dogs is that there is no reference specimen. Here we link a high-quality de novo long-read chromosomal assembly with epigenetic footprints and morphology to describe the Alpine dingo female named Cooinda. It was critical to establish an Alpine dingo reference because this ecotype occurs throughout coastal eastern Australia where the first drawings and descriptions were completed. FINDINGS: We generated a high-quality chromosome-level reference genome assembly (Canfam_ADS) using a combination of Pacific Bioscience, Oxford Nanopore, 10X Genomics, Bionano, and Hi-C technologies. Compared to the previously published Desert dingo assembly, there are large structural rearrangements on chromosomes 11, 16, 25, and 26. Phylogenetic analyses of chromosomal data from Cooinda the Alpine dingo and 9 previously published de novo canine assemblies show dingoes are monophyletic and basal to domestic dogs. Network analyses show that the mitochondrial DNA genome clusters within the southeastern lineage, as expected for an Alpine dingo. Comparison of regulatory regions identified 2 differentially methylated regions within glucagon receptor GCGR and histone deacetylase HDAC4 genes that are unmethylated in the Alpine dingo genome but hypermethylated in the Desert dingo. Morphologic data, comprising geometric morphometric assessment of cranial morphology, place dingo Cooinda within population-level variation for Alpine dingoes. Magnetic resonance imaging of brain tissue shows she had a larger cranial capacity than a similar-sized domestic dog. CONCLUSIONS: These combined data support the hypothesis that the dingo Cooinda fits the spectrum of genetic and morphologic characteristics typical of the Alpine ecotype. We propose that she be considered the archetype specimen for future research investigating the evolutionary history, morphology, physiology, and ecology of dingoes. The female has been taxidermically prepared and is now at the Australian Museum, Sydney.

The Australasian dingo archetype: De novo chromosome-length genome assembly, DNA methylome, and cranial morphology.

Background: One difficulty in testing the hypothesis that the Australasian dingo is a functional intermediate between wild wolves and domesticated breed dogs is that there is no reference specimen. Here we link a high-quality de novo long read chromosomal assembly with epigenetic footprints and morphology to describe the Alpine dingo female named Cooinda. It was critical to establish an Alpine dingo reference because this ecotype occurs throughout coastal eastern Australia where the first drawings and descriptions were completed. Findings: We generated a high-quality chromosome-level reference genome assembly (Canfam_ADS) using a combination of Pacific Bioscience, Oxford Nanopore, 10X Genomics, Bionano, and Hi-C technologies. Compared to the previously published Desert dingo assembly, there are large structural rearrangements on Chromosomes 11, 16, 25 and 26. Phylogenetic analyses of chromosomal data from Cooinda the Alpine dingo and nine previously published de novo canine assemblies show dingoes are monophyletic and basal to domestic dogs. Network analyses show that the mtDNA genome clusters within the southeastern lineage, as expected for an Alpine dingo. Comparison of regulatory regions identified two differentially methylated regions within glucagon receptor GCGR and histone deacetylase HDAC4 genes that are unmethylated in the Alpine dingo genome but hypermethylated in the Desert dingo. Morphological data, comprising geometric morphometric assessment of cranial morphology place dingo Cooinda within population-level variation for Alpine dingoes. Magnetic resonance imaging of brain tissue show she had a larger cranial capacity than a similar-sized domestic dog. Conclusions: These combined data support the hypothesis that the dingo Cooinda fits the spectrum of genetic and morphological characteristics typical of the Alpine ecotype. We propose that she be considered the archetype specimen for future research investigating the evolutionary history, morphology, physiology, and ecology of dingoes. The female has been taxidermically prepared and is now at the Australian Museum, Sydney.

Memory of stochastic single-cell apoptotic signaling promotes chemoresistance in neuroblastoma.

Gene expression noise is known to promote stochastic drug resistance through the elevated expression of individual genes in rare cancer cells. However, we now demonstrate that chemoresistant neuroblastoma cells emerge at a much higher frequency when the influence of noise is integrated across multiple components of an apoptotic signaling network. Using a JNK activity biosensor with longitudinal high-content and in vivo intravital imaging, we identify a population of stochastic, JNK-impaired, chemoresistant cells that exist because of noise within this signaling network. Furthermore, we reveal that the memory of this initially random state is retained following chemotherapy treatment across a series of in vitro, in vivo, and patient models. Using matched PDX models established at diagnosis and relapse from individual patients, we show that HDAC inhibitor priming cannot erase the memory of this resistant state within relapsed neuroblastomas but improves response in the first-line setting by restoring drug-induced JNK activity within the chemoresistant population of treatment-naïve tumors.