The elephant shark methylome reveals conservation of epigenetic regulation across jawed vertebrates

Evolution and Spatiotemporal Expression of ankha and ankhb in Zebrafish

Craniometaphyseal Dysplasia (CMD) is a rare skeletal disorder that can result from mutations in the ANKH gene. This gene encodes progressive anksylosis (ANK), which is responsible for transporting inorganic pyrophosphate (PPi) and ATP from the intracellular to the extracellular environment, where PPi inhibits bone mineralization. When ANK is dysfunctional, as in patients with CMD, the passage of PPi to the extracellular environment is reduced, leading to excess mineralization, particularly in bones of the skull. Zebrafish may serve as a promising model to study the mechanistic basis of CMD. Here, we provide a detailed analysis of the zebrafish Ankh paralogs, Ankha and Ankhb, in terms of their phylogenic relationship with ANK in other vertebrates as well as their spatiotemporal expression patterns during zebrafish development. We found that a closer evolutionary relationship exists between the zebrafish Ankhb protein and its human and other vertebrate counterparts, and stronger promoter activity was predicted for ankhb compared to ankha. Furthermore, we noted distinct temporal expression patterns, with ankha more prominently expressed in early development stages, and both paralogs also being expressed at larval growth stages. Whole-mount in situ hybridization was used to compare the spatial expression patterns of each paralog during bone development, and both showed strong expression in the craniofacial region as well as the notochord and somites. Given the substantial overlap in spatiotemporal expression but only subtle patterning differences, the exact roles of these genes remain speculative. In silico analyses predicted that Ankha and Ankhb have the same function in transporting PPi across the membrane. Nevertheless, this study lays the groundwork for functional analyses of each ankh paralog and highlights the potential of using zebrafish to find possible targeted therapies for CMD.

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.

Pronounced sequence specificity of the TET enzyme catalytic domain guides its cellular function

TET (ten-eleven translocation) enzymes catalyze the oxidation of 5-methylcytosine bases in DNA, thus driving active and passive DNA demethylation. Here, we report that the catalytic domain of mammalian TET enzymes favor CGs embedded within basic helix-loop-helix and basic leucine zipper domain transcription factor-binding sites, with up to 250-fold preference in vitro. Crystal structures and molecular dynamics calculations show that sequence preference is caused by intrasubstrate interactions and CG flanking sequence indirectly affecting enzyme conformation. TET sequence preferences are physiologically relevant as they explain the rates of DNA demethylation in TET-rescue experiments in culture and in vivo within the zygote and germ line. Most and least favorable TET motifs represent DNA sites that are bound by methylation-sensitive immediate-early transcription factors and octamer-binding transcription factor 4 (OCT4), respectively, illuminating TET function in transcriptional responses and pluripotency support.

Paternal hypoxia exposure primes offspring for increased hypoxia resistance

BACKGROUND

In a time of rapid environmental change, understanding how the challenges experienced by one generation can influence the fitness of future generations is critically needed. Using tolerance assays and transcriptomic and methylome approaches, we use zebrafish as a model to investigate cross-generational acclimation to hypoxia.

RESULTS

We show that short-term paternal exposure to hypoxia endows offspring with greater tolerance to acute hypoxia. We detected two hemoglobin genes that are significantly upregulated by more than 6-fold in the offspring of hypoxia exposed males. Moreover, the offspring which maintained equilibrium the longest showed greatest upregulation in hemoglobin expression. We did not detect differential methylation at any of the differentially expressed genes, suggesting that other epigenetic mechanisms are responsible for alterations in gene expression.

CONCLUSIONS

Overall, our findings suggest that an epigenetic memory of past hypoxia exposure is maintained and that this environmentally induced information is transferred to subsequent generations, pre-acclimating progeny to cope with hypoxic conditions.

Age group DNA methylation differences in lemon sharks (Negaprion brevirostris): Implications for future age estimation tools

Age information is often non-existent for most shark populations due to a lack of measurable physiological and morphological traits that can be used to estimate age. Recently, epigenetic clocks have been found to accurately estimate age for mammals, birds, and fish. However, since these clocks rely, among other things, on the availability of reference genomes, their application is hampered in non-traditional model organisms lacking such molecular resources. The technique known as Methyl-Sensitive Amplified Polymorphism (MSAP) has emerged as a valid alternative for studying DNA methylation biomarkers when reference genome information is missing, and large numbers of samples need to be processed. Accordingly, the MSAP technique was used in the present study to characterize global DNA methylation patterns in lemon sharks from three different age groups (juveniles, subadults, and adults). The obtained results reveal that, while MSAP analyses lack enough resolution as a standalone approach to infer age in these organisms, the global DNA methylation patterns observed using this technique displayed significant differences between age groups. Overall, these results confer that DNA methylation does change with age in sharks like what has been seen for other vertebrates and that MSAP could be useful as part of an epigenetics pipeline to infer the broad range of ages found in large samples sizes.

Emergence of distinct syntenic density regimes is associated with early metazoan genomic transitions

Background

Animal genomes are strikingly conserved in terms of local gene order (microsynteny). While some of these microsyntenies have been shown to be coregulated or to form gene regulatory blocks, the diversity of their genomic and regulatory properties across the metazoan tree of life remains largely unknown.

Results

Our comparative analyses of 49 animal genomes reveal that the largest gains of synteny occurred in the last common ancestor of bilaterians and cnidarians and in that of bilaterians. Depending on their node of emergence, we further show that novel syntenic blocks are characterized by distinct functional compositions (Gene Ontology terms enrichment) and gene density properties, such as high, average and low gene density regimes. This is particularly pronounced among bilaterian novel microsyntenies, most of which fall into high gene density regime associated with higher gene coexpression levels. Conversely, a majority of vertebrate novel microsyntenies display a low gene density regime associated with lower gene coexpression levels.

Conclusions

Our study provides first evidence for evolutionary transitions between different modes of microsyntenic block regulation that coincide with key events of metazoan evolution. Moreover, the microsyntenic profiling strategy and interactive online application (Syntenic Density Browser, available at: http://synteny.csb.univie.ac.at/) we present here can be used to explore regulatory properties of microsyntenic blocks and predict their coexpression in a wide-range of animal genomes.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-022-08304-2.

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.

A chicken DNA methylation clock for the prediction of broiler health

The domestic chicken (Gallus gallus domesticus) is the globally most important source of commercially produced meat. While genetic approaches have played an important role in the development of chicken stocks, little is known about chicken epigenetics. We have systematically analyzed the chicken DNA methylation machinery and DNA methylation landscape. While overall DNA methylation distribution was similar to mammals, sperm DNA appeared hypomethylated, which correlates with the absence of the DNMT3L cofactor in the chicken genome. Additional analysis revealed the presence of low-methylated regions, which are conserved gene regulatory elements that show tissue-specific methylation patterns. We also used whole-genome bisulfite sequencing to generate 56 single-base resolution methylomes from various tissues and developmental time points to establish an LMR-based DNA methylation clock for broiler chicken age prediction. This clock was used to demonstrate epigenetic age acceleration in animals with experimentally induced inflammation. Our study provides detailed insights into the chicken methylome and suggests a novel application of the DNA methylation clock as a marker for livestock health.

Raddatz, Lyko and colleagues use whole-genome bisulfite sequencing data to generate a methylation clock for chicken. This clock was able to detect age acceleration in broiler chickens under experimentally induced inflammation.

Hairpin-bisulfite sequencing of cells exposed to decitabine documents the process of DNA demethylation

Although the mechanism of DNA demethylating drugs has been understood for many years, the direct effect of these drugs on methylation of the complementary strands of DNA has not been formally demonstrated. By using hairpin-bisulphite sequencing, we describe the kinetics and pattern of DNA methylation following treatment of cells by the DNA methyltransferase 1 (DNMT1) inhibitor, decitabine. As expected, we demonstrate complete loss of methylation on the daughter strand following S-phase in selected densely methylated genes in synchronized Jurkat cells. Thereafter, cells showed a heterogeneous pattern of methylation reflecting replication of the unmethylated strand and restoration of methylation.

Cytosine methylation patterns suggest a role of methylation in plastic and adaptive responses to temperature in European grayling (Thymallus thymallus) populations

Temperature is a key environmental parameter affecting both the phenotypes and distributions of organisms, particularly ectotherms. Rapid organismal responses to thermal environmental changes have been described for several ectotherms; however, the underlying molecular mechanisms often remain unclear. Here, we studied whole genome cytosine methylation patterns of European grayling (Thymallus thymallus) embryos from five populations with contemporary adaptations of early life history traits at either 'colder' or 'warmer' spawning grounds. We reared fish embryos in a common garden experiment using two temperatures that resembled the 'colder' and 'warmer' conditions of the natal natural environments. Genome-wide methylation patterns were similar in populations originating from colder thermal origin subpopulations, whereas single nucleotide polymorphisms uncovered from the same data identified strong population structure among isolated populations, but limited structure among interconnected populations. This was surprising because the previously studied gene expression response among populations was mostly plastic, and mainly influenced by the developmental temperature. These findings support the hypothesis of the magnified role of epigenetic mechanisms in modulating plasticity. The abundance of consistently changing methylation loci between two warmer-to-colder thermal origin population pairs suggests that local adaptation has shaped the observed methylation patterns. The dynamic nature of the methylomes was further highlighted by genome-wide and site-specific plastic responses. Our findings support both the presence of a plastic response in a subset of CpG loci, and the evolutionary role of methylation divergence between populations adapting to contrasting thermal environments.

DNA methylation in the vertebrate germline: balancing memory and erasure

Cytosine methylation is a DNA modification that is critical for vertebrate development and provides a plastic yet stable information module in addition to the DNA code. DNA methylation memory establishment, maintenance and erasure is carefully balanced by molecular machinery highly conserved among vertebrates. In mammals, extensive erasure of epigenetic marks, including 5-methylcytosine (5mC), is a hallmark of early embryo and germline development. Conversely, global cytosine methylation patterns are preserved in at least some non-mammalian vertebrates over comparable developmental windows. The evolutionary mechanisms which drove this divergence are unknown, nevertheless a direct consequence of retaining epigenetic memory in the form of 5mC is the enhanced potential for transgenerational epigenetic inheritance (TEI). Given that DNA methylation dynamics remains underexplored in most vertebrate lineages, the extent of information transferred to offspring by epigenetic modification might be underestimated.

Zebrafish preserve global germline DNA methylation while sex-linked rDNA is amplified and demethylated during feminisation

The germline is the only cellular lineage capable of transferring genetic information from one generation to the next. Intergenerational transmission of epigenetic memory through the germline, in the form of DNA methylation, has been proposed; however, in mammals this is largely prevented by extensive epigenetic erasure during germline definition. Here we report that, unlike mammals, the continuously-defined ‘preformed’ germline of zebrafish does not undergo genome-wide erasure of DNA methylation during development. Our analysis also uncovers oocyte-specific germline amplification and demethylation of an 11.5-kb repeat region encoding 45S ribosomal RNA (fem-rDNA). The peak of fem-rDNA amplification coincides with the initial expansion of stage IB oocytes, the poly-nucleolar cell type responsible for zebrafish feminisation. Given that fem-rDNA overlaps with the only zebrafish locus identified thus far as sex-linked, we hypothesise fem-rDNA expansion could be intrinsic to sex determination in this species.

Germline cells transfer genetic information to offspring, and in zebrafish, drive sex determination. Here the authors report that, unlike mammals, the germline of zebrafish does not undergo genome-wide DNA methylation erasure, while amplifying and demethylating sex-linked rDNA during feminisation.

Stress, novel sex genes, and epigenetic reprogramming orchestrate socially controlled sex change

Bluehead wrasses undergo dramatic, socially cued female-to-male sex change. We apply transcriptomic and methylome approaches in this wild coral reef fish to identify the primary trigger and subsequent molecular cascade of gonadal metamorphosis. Our data suggest that the environmental stimulus is exerted via the stress axis and that repression of the aromatase gene (encoding the enzyme converting androgens to estrogens) triggers a cascaded collapse of feminizing gene expression and identifies notable sex-specific gene neofunctionalization. Furthermore, sex change involves distinct epigenetic reprogramming and an intermediate state with altered epigenetic machinery expression akin to the early developmental cells of mammals. These findings reveal at a molecular level how a normally committed developmental process remains plastic and is reversed to completely alter organ structures.

Bridging disciplines to advance elasmobranch conservation: applications of physiological ecology

A strength of physiological ecology is its incorporation of aspects of both species' ecology and physiology; this holistic approach is needed to address current and future anthropogenic stressors affecting elasmobranch fishes that range from overexploitation to the effects of climate change. For example, physiology is one of several key determinants of an organism's ecological niche (along with evolutionary constraints and ecological interactions). The fundamental role of physiology in niche determination led to the development of the field of physiological ecology. This approach considers physiological mechanisms in the context of the environment to understand mechanistic variations that beget ecological trends. Physiological ecology, as an integrative discipline, has recently experienced a resurgence with respect to conservation applications, largely in conjunction with technological advances that extended physiological work from the lab into the natural world. This is of critical importance for species such as elasmobranchs (sharks, skates and rays), which are an especially understudied and threatened group of vertebrates. In 2017, at the American Elasmobranch Society meeting in Austin, Texas, the symposium entitled `Applications of Physiological Ecology in Elasmobranch Research' provided a platform for researchers to showcase work in which ecological questions were examined through a physiological lens. Here, we highlight the research presented at this symposium, which emphasized the strength of linking physiological tools with ecological questions. We also demonstrate the applicability of using physiological ecology research as a method to approach conservation issues, and advocate for a more available framework whereby results are more easily accessible for their implementation into management practices.

Molecular evolution of DNMT1 in vertebrates: Duplications in marsupials followed by positive selection

DNA methylation is mediated by a conserved family of DNA methyltransferases (Dnmts). The human genome encodes three active Dnmts (Dnmt1, Dnmt3a and Dnmt3b), the tRNA methyltransferase Dnmt2, and the regulatory protein Dnmt3L. Despite their high degree of conservation among different species, genes encoding Dnmts have been duplicated and/or lost in multiple lineages throughout evolution, indicating that the DNA methylation machinery has some potential to undergo evolutionary change. However, little is known about the extent to which this machinery, or the methylome, varies among vertebrates. Here, we study the molecular evolution of Dnmt1, the enzyme responsible for maintenance of DNA methylation patterns after replication, in 79 vertebrate species. Our analyses show that all studied species exhibit a single copy of the DNMT1 gene, with the exception of tilapia and marsupials (tammar wallaby, koala, Tasmanian devil and opossum), each of which displays two apparently functional DNMT1 copies. Our phylogenetic analyses indicate that DNMT1 duplicated before the radiation of major marsupial groups (i.e., at least ~75 million years ago), thus giving rise to two DNMT1 copies in marsupials (copy 1 and copy 2). In the opossum lineage, copy 2 was lost, and copy 1 recently duplicated again, generating three DNMT1 copies: two putatively functional genes (copy 1a and 1b) and one pseudogene (copy 1ψ). Both marsupial copies (DNMT1 copies 1 and 2) are under purifying selection, and copy 2 exhibits elevated rates of evolution and signatures of positive selection, suggesting a scenario of neofunctionalization. This gene duplication might have resulted in modifications in marsupial methylomes and their dynamics.

Evolutionary conservation of DNA methylation in CpG sites within ultraconserved noncoding elements

Ultraconserved noncoding elements (UCNEs) constitute less than 1 Mb of vertebrate genomes and are impervious to accumulating mutations. About 4000 UCNEs exist in vertebrate genomes, each at least 200 nucleotides in length, sharing greater than 95% sequence identity between human and chicken. Despite extreme sequence conservation over 400 million years of vertebrate evolution, we show both ordered interspecies and within-species interindividual variation in DNA methylation in these regions. Here, we surveyed UCNEs with high CpG density in 56 species finding half to be intermediately methylated and the remaining near 0% or 100%. Intermediately methylated UCNEs displayed a greater range of methylation between mouse tissues. In a human population, most UCNEs showed greater variation than the LINE1 transposon, a frequently used epigenetic biomarker. Global methylation was found to be inversely correlated to hydroxymethylation across 60 vertebrates. Within UCNEs, DNA methylation is flexible, conserved between related species, and relaxed from the underlying sequence selection pressure, while remaining heritable through speciation.