Long-Term m5C Methylome Dynamics Parallel Phenotypic Adaptation in the Cyanobacterium Trichodesmium

A major challenge in modern biology is understanding how the effects of short-term biological responses influence long-term evolutionary adaptation, defined as a genetically determined increase in fitness to novel environments. This is particularly important in globally important microbes experiencing rapid global change, due to their influence on food webs, biogeochemical cycles, and climate. Epigenetic modifications like methylation have been demonstrated to influence short-term plastic responses, which ultimately impact long-term adaptive responses to environmental change. However, there remains a paucity of empirical research examining long-term methylation dynamics during environmental adaptation in nonmodel, ecologically important microbes. Here, we show the first empirical evidence in a marine prokaryote for long-term m5C methylome modifications correlated with phenotypic adaptation to CO₂, using a 7-year evolution experiment (1,000+ generations) with the biogeochemically important marine cyanobacterium Trichodesmium. We identify m5C methylated sites that rapidly changed in response to high (750 µatm) CO₂ exposure and were maintained for at least 4.5 years of CO₂ selection. After 7 years of CO₂ selection, however, m5C methylation levels that initially responded to high-CO₂ returned to ancestral, ambient CO₂ levels. Concurrently, high-CO₂ adapted growth and N₂ fixation rates remained significantly higher than those of ambient CO₂ adapted cell lines irrespective of CO₂ concentration, a trend consistent with genetic assimilation theory. These data demonstrate the maintenance of CO₂-responsive m5C methylation for 4.5 years alongside phenotypic adaptation before returning to ancestral methylation levels. These observations in a globally distributed marine prokaryote provide critical evolutionary insights into biogeochemically important traits under global change.     

Fine‐scale population epigenetic structure in relation to gastrointestinal parasite load in red grouse (Lagopus lagopus scotica)

Epigenetic modification of cytosine methylation states can be elicited by environmental stresses and may be a key process affecting phenotypic plasticity and adaptation. Parasites are potent stressors with profound physiological and ecological effects on their host, but there is little understanding in how parasites may influence host methylation states. Here, we estimate epigenetic diversity and differentiation among 21 populations of red grouse (Lagopus lagopus scotica) in north‐east Scotland and test for association of gastrointestinal parasite load (caecal nematode Trichostrongylus tenuis) with hepatic genome‐wide and locus‐specific methylation states. Following methylation‐sensitive AFLP (MSAP), 129 bands, representing 73 methylation‐susceptible and 56 nonmethylated epiloci, were scored across 234 individuals. The populations differed significantly in genome‐wide methylation levels and were also significantly epigenetically (F_SC = 0.0227; P < 0.001) and genetically (F_SC = 0.0058; P < 0.001) differentiated. Parasite load was not associated with either genome‐wide methylation levels or epigenetic differentiation. Instead, we found eight disproportionately differentiated epilocus‐specific methylation states (F_ST outliers) using bayescan software and significant positive and negative association of 35 methylation states with parasite load from bespoke generalized estimating equations (GEE), simple logistic regression (sam ) and Bayesian environmental analysis (bayenv 2). Following Sanger sequencing, genome mapping and geneontology (go ) annotation, some of these epiloci were linked to genes involved in regulation of cell cycle, signalling, metabolism, immune system and notably rRNA methylation, histone acetylation and small RNAs. These findings demonstrate an epigenetic signature of parasite load in populations of a wild bird and suggest intriguing physiological effects of parasite‐associated cytosine methylation.     

Methylation analysis revealed salicylic acid affects pearl millet defense through external cytosine DNA demethylation

The cytosine DNA methylation and demethylation have a role in regulating plant responses to the environment by affecting the promoter regions of most plant defense-related genes through the CpG islands or the CCGG motifs. Salicylic acid, a defense and signaling plant hormone, is seen playing crucial role in the variation of the methylome. In this study, the effects of salicylic acid and feeding of the millet headminer (Heliocheilus albipunctella de Joannis) on pearl millet DNA methylome changes were evaluated through MSAP epigenotyping during panicle development. The results showed that millet headminer feeding increased the level of genomic methylation while application of salicylic acid caused DNA demethylation occurring mostly at external cytosine and accompanied by a decrease of the number of larvae per panicle. This suggests that hemimethylation (external cytosine methylation) has key role in regulating defense responses and conferring tolerance to pearl millet through salicylic acid application.     

MSAP markers and global cytosine methylation in plants: a literature survey and comparative analysis for a wild‐growing species

Methylation of DNA cytosines affects whether transposons are silenced and genes are expressed, and is a major epigenetic mechanism whereby plants respond to environmental change. Analyses of methylation‐sensitive amplification polymorphism (MS‐AFLP or MSAP) have been often used to assess methyl‐cytosine changes in response to stress treatments and, more recently, in ecological studies of wild plant populations. MSAP technique does not require a sequenced reference genome and provides many anonymous loci randomly distributed over the genome for which the methylation status can be ascertained. Scoring of MSAP data, however, is not straightforward, and efforts are still required to standardize this step to make use of the potential to distinguish between methylation at different nucleotide contexts. Furthermore, it is not known how accurately MSAP infers genome‐wide cytosine methylation levels in plants. Here, we analyse the relationship between MSAP results and the percentage of global cytosine methylation in genomic DNA obtained by HPLC analysis. A screening of literature revealed that methylation of cytosines at cleavage sites assayed by MSAP was greater than genome‐wide estimates obtained by HPLC, and percentages of methylation at different nucleotide contexts varied within and across species. Concurrent HPLC and MSAP analyses of DNA from 200 individuals of the perennial herb Helleborus foetidus confirmed that methyl‐cytosine was more frequent in CCGG contexts than in the genome as a whole. In this species, global methylation was unrelated to methylation at the inner CG site. We suggest that global HPLC and context‐specific MSAP methylation estimates provide complementary information whose combination can improve our current understanding of methylation‐based epigenetic processes in nonmodel plants.     

Dnmt2 is the most evolutionary conserved and enigmatic cytosine DNA methyltransferase in eukaryotes

Dnmt2 is the most strongly conserved cytosine DNA methyltransferase in eukaryotes. It has been found in all organisms possessing methyltransferases of the Dnmt1 and Dnmt3 families, whereas in many others Dnmt2 is the sole cytosine DNA methyltransferase. The Dnmt2 molecule contains all conserved motifs of cytosine DNA methyltransferases. It forms 3D complexes with DNA very similar to those of bacterial DNA methyltransferases and performs cytosine methylation by a catalytic mechanism common to all cytosine DNA methyltransferases. Catalytic activity of the purified Dnmt2 with DNA substrates is very low and could hardly be detected in direct biochemical assays. Dnmt2 is the sole cytosine DNA methyltransferase in Drosophila and other dipteran insects. Its overexpression as a transgene leads to DNA hypermethylation in all sequence contexts and to an extended life span. On the contrary, a null-mutation of the Dnmt2 gene leads to a diminished life span, though no evident anomalies in development are observed. Dnmt2 is also the sole cytosine DNA methyltransferase in several protists. Similar to Drosophila these protists have a very low level of DNA methylation. Some limited genome compartments, such as transposable sequences, are probably the methylation targets in these organisms. Dnmt2 does not participate in genome methylation in mammals, but seems to be an RNA methyltransferase modifying the 38th cytosine residue in anticodon loop of certain tRNAs. This modification enhances stability of tRNAs, especially in stressful conditions. Dnmt2 is the only enzyme known to perform RNA methylation by a catalytic mechanism characteristic of DNA methyltransferases. The Dnmt2 activity has been shown in mice to be necessary for paramutation establishment, though the precise mechanisms of its participation in this form of epigenetic heredity are unknown. It seems likely, that either of the two Dnmt2 activities could become a predominant one during the evolution of different species. The high level of the Dnmt2 evolutionary conservation proves its activity to have a significant adaptive value in natural environment.     

Genome-Wide Methylation Patterns in Salmonella enterica Subsp. enterica Serovars

The methylation of DNA bases plays an important role in numerous biological processes including development, gene expression, and DNA replication. Salmonella is an important foodborne pathogen, and methylation in Salmonella is implicated in virulence. Using single molecule real-time (SMRT) DNA-sequencing, we sequenced and assembled the complete genomes of eleven Salmonella enterica isolates from nine different serovars, and analysed the whole-genome methylation patterns of each genome. We describe 16 distinct N⁶-methyladenine (m6A) methylated motifs, one N⁴-methylcytosine (m4C) motif, and one combined m6A-m4C motif. Eight of these motifs are novel, i.e., they have not been previously described. We also identified the methyltransferases (MTases) associated with 13 of the motifs. Some motifs are conserved across all Salmonella serovars tested, while others were found only in a subset of serovars. Eight of the nine serovars contained a unique methylated motif that was not found in any other serovar (most of these motifs were part of Type I restriction modification systems), indicating the high diversity of methylation patterns present in Salmonella.     

Chromosome‐level genome assembly of the razor clam Sinonovacula constricta (Lamarck, 1818)

Bivalves, a highly diverse and the most evolutionarily successful class of invertebrates native to aquatic habitats, provide valuable molecular resources for understanding the evolutionary adaptation and aquatic ecology. Here, we reported a high‐quality chromosome‐level genome assembly of the razor clam Sinonovacula constricta using Pacific Bioscience single‐molecule real‐time sequencing, Illumina paired‐end sequencing, 10X Genomics linked‐reads and Hi‐C reads. The genome size was 1,220.85 Mb, containing scaffold N50 of 65.93 Mb and contig N50 of 976.94 Kb. A total of 899 complete (91.92%) and seven partial (0.72%) matches of the 978 metazoa Benchmarking Universal Single‐Copy Orthologs were determined in this genome assembly. And Hi‐C scaffolding of the genome resulted in 19 pseudochromosomes. A total of 28,594 protein‐coding genes were predicted in the S. constricta genome, of which 25,413 genes (88.88%) were functionally annotated. In addition, 39.79% of the assembled genome was composed of repetitive sequences, and 4,372 noncoding RNAs were identified. The enrichment analyses of the significantly expanded and contracted genes suggested an evolutionary adaptation of S. constricta to highly stressful living environments. In summary, the genomic resources generated in this work not only provide a valuable reference genome for investigating the molecular mechanisms of S. constricta biological functions and evolutionary adaptation, but also facilitate its genetic improvement and disease treatment. Meanwhile, the obtained genome greatly improves our understanding of the genetics of molluscs and their comparative evolution.     

Stability of the cytosine methylome during post-testicular sperm maturation in mouse

Beyond the haploid genome, mammalian sperm carry a payload of epigenetic information with the potential to modulate offspring phenotypes. Recent studies show that the small RNA repertoire of sperm is remodeled during post-testicular maturation in the epididymis. Epididymal maturation has also been linked to changes in the sperm methylome, suggesting that the epididymis might play a broader role in shaping the sperm epigenome. Here, we characterize the genome-wide methylation landscape in seven germ cell populations from throughout the male reproductive tract. We find very few changes in the cytosine methylation landscape between testicular germ cell populations and cauda epididymal sperm, demonstrating that the sperm methylome is stable throughout post-testicular maturation. Although our sequencing data suggested that caput epididymal sperm exhibit a highly unusual methylome, follow-up studies revealed that this resulted from contamination of caput sperm by extracellular DNA. Extracellular DNA formed web-like structures that ensnared sperm, and was present only in sperm samples obtained from the caput epididymis and vas deferens of virgin males. Curiously, contaminating extracellular DNA was associated with citrullinated histone H3, potentially resulting from a PAD-driven genome decondensation process. Taken together, our data emphasize the stability of cytosine methylation in mammalian sperm, and identify a surprising, albeit transient, period during which sperm are associated with extracellular DNA.     

Characterization of the level,target sites and inheritance of cytosine methylation in tomato nuclear DNA

The tomato nuclear genome was determined to have a G+C content of 37% which is among the lowest reported for any plant species. Non-coding regions have a G+C content even lower (32% average) whereas coding regions are considerably richer in G+C (46%).5-methyl cytosine was the only modified base detected and on average 23% of the cytosine residues are methylated. Immature tissues and protoplasts have significantly lower levels of cytosine methylation (average 20%) than mature tissues (average 25%). Mature pollen has an intermediate level of methylation (22%). Seeds gave the highest value (27%), suggesting de novo methylation after pollination and during seed development.Based on isoschizomer studies we estimate 55% of the CpG target sites (detected by Msp I/Hpa II) and 85% of the CpNpG target sites (detected by Bst NI/Eco RI)are methylated. Unmethylated target sites (both CpG and CpNpG) are not randomly distributed throughout the genome, but frequently occur in clusters. These clusters resemble CpG islands recently reported in maize and tobacco.The low G+C content and high levels of cytosine methylation in tomato may be due to previous transitions of 5mCT. This is supported by the fact that G+C levels are lowest in non-coding portions of the genome in which selection is relaxed and thus transitions are more likely to be tolerated. This hypothesis is also supported by the general deficiency of methylation target sites in the tomato genome, especially in non-coding regions.Using methylation isoschizomers and RFLP analysis we have also determined that polymorphism between plants, for cytosine methylation at allelic sites, is common in tomato. Comparing DNA from two tomato species, 20% of the polymorphisms detected by Bst NI/Eco RII could be attributed to differential methylation at the CpNpG target sites. With Msp I/Hpa II, 50% of the polymorphisms were attributable to methylation (CpG and CpNpG sites). Moreover, these polymorphisms were demonstrated to be inherited in a mendelian fashion and to co-segregate with the methylation target site and thus do not represent variation for transacting factors that might be involved in methylation of DNA. The potential role of heritable methylation polymorphism in evolution of gene regulation and in RFLP studies is discussed.     

Borderline personality disorder and childhood maltreatment: a genome‐wide methylation analysis

Early life adversity plays a critical role in the emergence of borderline personality disorder (BPD) and this could occur through epigenetic programming. In this perspective, we aimed to determine whether childhood maltreatment could durably modify epigenetic processes by the means of a whole‐genome methylation scan of BPD subjects. Using the Illumina Infinium® HumanMethylation450 BeadChip, global methylation status of DNA extracted from peripheral blood leucocytes was correlated to the severity of childhood maltreatment in 96 BPD subjects suffering from a high level of child adversity and 93 subjects suffering from major depressive disorder (MDD) and reporting a low rate of child maltreatment. Several CpGs within or near the following genes (IL17RA, miR124‐3, KCNQ2, EFNB1, OCA2, MFAP2, RPH3AL, WDR60, CST9L, EP400, A2ML1, NT5DC2, FAM163A and SPSB2) were found to be differently methylated, either in BPD compared with MDD or in relation to the severity of childhood maltreatment. A highly relevant biological result was observed for cg04927004 close to miR124‐3 that was significantly associated with BPD and severity of childhood maltreatment. miR124‐3 codes for a microRNA (miRNA) targeting several genes previously found to be associated with BPD such as NR3C1. Our results highlight the potentially important role played by miRNAs in the etiology of neuropsychiatric disorders such as BPD and the usefulness of using methylome‐wide association studies to uncover such candidate genes. Moreover, they offer new understanding of the impact of maltreatments on biological processes leading to diseases and may ultimately result in the identification of relevant biomarkers.     

Genome-Wide Analysis of DNA Methylation Dynamics during Early Human Development

DNA methylation is globally reprogrammed during mammalian preimplantation development, which is critical for normal development. Recent reduced representation bisulfite sequencing (RRBS) studies suggest that the methylome dynamics are essentially conserved between human and mouse early embryos. RRBS is known to cover 5–10% of all genomic CpGs, favoring those contained within CpG-rich regions. To obtain an unbiased and more complete representation of the methylome during early human development, we performed whole genome bisulfite sequencing of human gametes and blastocysts that covered>70% of all genomic CpGs. We found that the maternal genome was demethylated to a much lesser extent in human blastocysts than in mouse blastocysts, which could contribute to an increased number of imprinted differentially methylated regions in the human genome. Global demethylation of the paternal genome was confirmed, but SINE-VNTR-Alu elements and some other tandem repeat-containing regions were found to be specifically protected from this global demethylation. Furthermore, centromeric satellite repeats were hypermethylated in human oocytes but not in mouse oocytes, which might be explained by differential expression of de novo DNA methyltransferases. These data highlight both conserved and species-specific regulation of DNA methylation during early mammalian development. Our work provides further information critical for understanding the epigenetic processes underlying differentiation and pluripotency during early human development.     

Genome-wide antagonism between 5-hydroxymethylcytosine and DNA methylation in the adult mouse brain

Mounting evidence points to critical roles for DNA modifications, including 5-methylcytosine （5mC） and its oxidized forms, in the development, plasticity and disorders of the mammalian nervous system. The novel DNA base 5- hydroxymethylcytosine （5hmC） is known to be capable of initiating passive or active DNA demethylation, but whether and how extensively 5hmC functions in shaping the post-mitotic neuronal DNA methylome is unclear. Here we report the genome-wide distribution of 5hmC in dentate granule neurons from adult mouse hippocampus in vivo. 5hmC in the neuronal genome is highly enriched in gene bodies, especially in exons, and correlates with gene expression. Direct genome-wide comparison of 5hmC distribution between embryonic stem cells and neurons reveals extensive differences, reflecting the functional disparity between these two cell types. Importantly, integrative analysis of 5hmC, overall DNA methylation and gene expression profiles of dentate granule neurons in vivo reveals the genome-wide antagonism between these two states of cytosine modifications, supporting a role for 5hmC in shaping the neuronal DNA methylome by promoting active DNA demethylation.     

Endogenous pararetrovirus sequences associated with 24 nt small RNAs at the centromeres of Fritillaria imperialis L. (Liliaceae), a species with a giant genome

Endogenous pararetroviral sequences are the most commonly found virus sequences integrated into angiosperm genomes. We describe an endogenous pararetrovirus (EPRV) repeat in Fritillaria imperialis, a species that is under study as a result of its exceptionally large genome (1C = 42 096 Mbp, approximately 240 times bigger than Arabidopsis thaliana). The repeat (FriEPRV) was identified from Illumina reads using the RepeatExplorer pipeline, and exists in a complex genomic organization at the centromere of most, or all, chromosomes. The repeat was reconstructed into three consensus sequences that formed three interconnected loops, one of which carries sequence motifs expected of an EPRV (including the gag and pol domains). FriEPRV shows sequence similarity to members of the Caulimoviridae pararetrovirus family, with phylogenetic analysis indicating a close relationship to Petuvirus. It is possible that no complete EPRV sequence exists, although our data suggest an abundance that exceeds the genome size of Arabidopsis. Analysis of single nucleotide polymorphisms revealed elevated levels of C→T and G→A transitions, consistent with deamination of methylated cytosine. Bisulphite sequencing revealed high levels of methylation at CG and CHG motifs (up to 100%), and 15–20% methylation, on average, at CHH motifs. FriEPRV's centromeric location may suggest targeted insertion, perhaps associated with meiotic drive. We observed an abundance of 24 nt small RNAs that specifically target FriEPRV, potentially providing a signature of RNA‐dependent DNA methylation. Such signatures of epigenetic regulation suggest that the huge genome of F. imperialis has not arisen as a consequence of a catastrophic breakdown in the regulation of repeat amplification.     

The CHH motif in sugar beet satellite DNA: a modulator for cytosine methylation

Methylation of DNA is important for the epigenetic silencing of repetitive DNA in plant genomes. Knowledge about the cytosine methylation status of satellite DNAs, a major class of repetitive DNA, is scarce. One reason for this is that arrays of tandemly arranged sequences are usually collapsed in next‐generation sequencing assemblies. We applied strategies to overcome this limitation and quantified the level of cytosine methylation and its pattern in three satellite families of sugar beet (Beta vulgaris) which differ in their abundance, chromosomal localization and monomer size. We visualized methylation levels along pachytene chromosomes with respect to small satellite loci at maximum resolution using chromosome‐wide fluorescent in situ hybridization complemented with immunostaining and super‐resolution microscopy. Only reduced methylation of many satellite arrays was obtained. To investigate methylation at the nucleotide level we performed bisulfite sequencing of 1569 satellite sequences. We found that the level of methylation of cytosine strongly depends on the sequence context: cytosines in the CHH motif show lower methylation (44–52%), while CG and CHG motifs are more strongly methylated. This affects the overall methylation of satellite sequences because CHH occurs frequently while CG and CHG are rare or even absent in the satellite arrays investigated. Evidently, CHH is the major target for modulation of the cytosine methylation level of adjacent monomers within individual arrays and contributes to their epigenetic function. This strongly indicates that asymmetric cytosine methylation plays a role in the epigenetic modification of satellite repeats in plant genomes.