DNA methylation profiling in mummified human remains from the eighteenth-century

Reconstruction of ancient epigenomes by DNA methylation (DNAm) can shed light into the composition of cell types, disease states, and age at death. However, such analysis is hampered by impaired DNA quality and little is known how decomposition affects DNAm. In this study, we determined if EPIC Illumina BeadChip technology is applicable for specimens from mummies of the eighteenth century CE. Overall, the signal intensity on the microarray was extremely low, but for one of two samples we were able to detect characteristic DNAm signals in a subset of CG dinucleotides (CpGs), which were selected with a stringent processing pipeline. Using only these CpGs we could train epigenetic signatures with reference DNAm profiles of multiple tissues and our predictions matched the fact that the specimen was lung tissue from a 28-year-old woman. Thus, we provide proof of principle that Illumina BeadChips are applicable for DNAm profiling in ancient samples.

CG dinucleotide removal in bioluminescent and fluorescent reporters improves HIV-1 replication and reporter gene expression for dual imaging in humanized mice

Visualizing the transmission and dissemination of human immunodeficiency virus type 1 (HIV-1) in real time in humanized mouse models is a robust tool to investigate viral replication during treatments and in tissue reservoirs. However, the stability and expression of HIV-1 reporter genes are obstacles for long-term serial imaging in vivo . Two replication-competent CCR5-tropic HIV-1 reporter constructs were created that encode either nanoluciferase (nLuc) or a near - infrared fluorescent protein (iRFP) upstream of nef . HIV-1 reporter virus replication and reporter gene expression was measured in cell culture and in humanized mice. While reporter gene expression in vivo correlated initially with plasma viremia, expression decreased after 4-5 weeks despite high plasma viremia. The reporter genes were codon-optimized to remove cytosine/guanine (CG) dinucleotides and new CO-nLuc and CO-iRFP viruses were reconstructed. Removal of CG dinucleotides in HIV-1 reporter viruses improved replication in vitro and reporter expression in vivo and ex vivo . Both codon optimized reporter viruses could be visualized during co-infection and in vivo reporter gene expression during treatment failure preceded detection of plasma viremia. While the dynamic range of CO-iRFP HIV-1 was lower than that of CO-nLuc HIV-1, both viruses could have utility in studying and visualizing HIV-1 infection in humanized mice. Importance Animal models are important for studying HIV-1 pathogenesis and treatments. We developed two viruses each encoding a reporter gene that can be expressed in cells after infection. This study shows that HIV-1 infection can be visualized by noninvasive, whole body imaging in mice with human immune cells over time by reporter expression. We improved reporter expression to reflect HIV-1 replication and showed that two viral variants can be tracked over time in the same animal and can predict failure of antiretroviral therapy to suppress virus.

Investigation of measurable residual disease in acute myeloid leukemia by DNA methylation patterns

Assessment of measurable residual disease (MRD) upon treatment of acute myeloid leukemia (AML) remains challenging. It is usually addressed by highly sensitive PCR- or sequencing-based screening of specific mutations, or by multiparametric flow cytometry. However, not all patients have suitable mutations and heterogeneity of surface markers hampers standardization in clinical routine. In this study, we propose an alternative approach to estimate MRD based on AML-associated DNA methylation (DNAm) patterns. We identified four CG dinucleotides (CpGs) that commonly reveal aberrant DNAm in AML and their combination could reliably discern healthy and AML samples. Interestingly, bisulfite amplicon sequencing demonstrated that aberrant DNAm patterns were symmetric on both alleles, indicating that there is epigenetic crosstalk between homologous chromosomes. We trained shallow-learning and deep-learning algorithms to identify anomalous DNAm patterns. The method was then tested on follow-up samples with and without MRD. Notably, even samples that were classified as MRD negative often revealed higher anomaly ratios than healthy controls, which may reflect clonal hematopoiesis. Our results demonstrate that targeted DNAm analysis facilitates reliable discrimination of malignant and healthy samples. However, since healthy samples also comprise few abnormal-classified DNAm reads the approach does not yet reliably discriminate MRD positive and negative samples.

Multi-species and multi-tissue methylation clocks for age estimation in toothed whales and dolphins

The development of a precise blood or skin tissue DNA Epigenetic Aging Clock for Odontocete (OEAC) would solve current age estimation inaccuracies for wild odontocetes. Therefore, we determined genome-wide DNA methylation profiles using a custom array (HorvathMammalMethyl40) across skin and blood samples (n = 446) from known age animals representing nine odontocete species within 4 phylogenetic families to identify age associated CG dinucleotides (CpGs). The top CpGs were used to create a cross-validated OEAC clock which was highly correlated for individuals (r = 0.94) and for unique species (median r = 0.93). Finally, we applied the OEAC for estimating the age and sex of 22 wild Norwegian killer whales. DNA methylation patterns of age associated CpGs are highly conserved across odontocetes. These similarities allowed us to develop an odontocete epigenetic aging clock (OEAC) which can be used for species conservation efforts by provide a mechanism for estimating the age of free ranging odontocetes from either blood or skin samples.

DNA methylation changes during long-term in vitro cell culture are caused by epigenetic drift

Culture expansion of primary cells evokes highly reproducible DNA methylation (DNAm) changes. We have identified CG dinucleotides (CpGs) that become continuously hyper- or hypomethylated during long-term culture of mesenchymal stem cells (MSCs) and other cell types. Bisulfite barcoded amplicon sequencing (BBA-seq) demonstrated that DNAm patterns of neighboring CpGs become more complex without evidence of continuous pattern development and without association to oligoclonal subpopulations. Circularized chromatin conformation capture (4C) revealed reproducible changes in nuclear organization between early and late passages, while there was no enriched interaction with other genomic regions that also harbor culture-associated DNAm changes. Chromatin immunoprecipitation of CTCF did not show significant differences during long-term culture of MSCs, however culture-associated hypermethylation was enriched at CTCF binding sites and hypomethylated CpGs were devoid of CTCF. Taken together, our results support the notion that DNAm changes during culture-expansion are not directly regulated by a targeted mechanism but rather resemble epigenetic drift.

Region-level epimutation rates in Arabidopsis thaliana

Failure to maintain DNA methylation patterns during plant development can occasionally give rise to so-called “spontaneous epimutations”. These stochastic methylation changes are sometimes heritable across generations and thus accumulate in plant genomes over time. Recent evidence indicates that spontaneous epimutations have a major role in shaping patterns of methylation diversity in plant populations. Using single CG dinucleotides as units of analysis, previous work has shown that the epimutation rate is several orders of magnitude higher than the genetic mutation rate. While these large rate differences have obvious implications for understanding genome-methylome co-evolution, the functional relevance of single CG methylation changes remains questionable. In contrast to single CG, solid experimental evidence has linked methylation gains and losses in larger genomic regions with transcriptional variation and heritable phenotypic effects. Here we show that such region-level changes arise stochastically at about the same rate as those at individual CG sites, are only marginal dependent on region size and cytosine density, but strongly dependent on chromosomal location. We also find consistent evidence that region-level epimutations are not restricted to CG contexts but also frequently occur in non-CG regions at the genome-wide scale. Taken together, our results support the view that many differentially methylated regions (DMRs) in natural populations originate from epimutation events and may not be effectively tagged by proximal SNPs. This possibility reinforces the need for epigenome-wide association studies (EWAS) in plants as a way to identify the epigenetic basis of complex traits.

Developmental Accumulation of Gene Body and Transposon Non-CpG Methylation in the Zebrafish Brain

DNA methylation predominantly occurs at CG dinucleotides in vertebrate genomes; however, non-CG methylation (mCH) is also detectable in vertebrate tissues, most notably in the nervous system. In mammals it is well established that mCH is targeted to CAC trinucleotides by DNMT3A during nervous system development where it is enriched in gene bodies and associated with transcriptional repression. Nevertheless, the conservation of developmental mCH accumulation and its deposition by DNMT3A is largely unexplored and has yet to be functionally demonstrated in other vertebrates. In this study, by analyzing DNA methylomes and transcriptomes of zebrafish brains, we identified enrichment of mCH at CAC trinucleotides (mCAC) at defined transposon motifs as well as in developmentally downregulated genes associated with developmental and neural functions. We further generated and analyzed DNA methylomes and transcriptomes of developing zebrafish larvae and demonstrated that, like in mammals, mCH accumulates during post-embryonic brain development. Finally, by employing CRISPR/Cas9 technology, we unraveled a conserved role for Dnmt3a enzymes in developmental mCAC deposition. Overall, this work demonstrates the evolutionary conservation of developmental mCH dynamics and highlights the potential of zebrafish as a model to study mCH regulation and function during normal and perturbed development.

The Long-Term Evolutionary History of Gradual Reduction of CpG Dinucleotides in the SARS-CoV-2 Lineage

Recent studies suggested that the fraction of CG dinucleotides (CpG) is severely reduced in the genome of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2). The CpG deficiency was predicted to be the adaptive response of the virus to evade degradation of the viral RNA by the antiviral zinc finger protein that specifically binds to CpG nucleotides. By comparing all representative genomes belonging to the genus Betacoronavirus, this study examined the potential time of origin of CpG depletion. The results of this investigation revealed a highly significant correlation between the proportions of CpG nucleotide (CpG content) of the betacoronavirus species and their times of divergence from SARS-CoV-2. Species that are distantly related to SARS-CoV-2 had much higher CpG contents than that of SARS-CoV-2. Conversely, closely related species had low CpG contents that are similar to or slightly higher than that of SARS-CoV-2. These results suggest a systematic and continuous reduction in the CpG content in the SARS-CoV-2 lineage that might have started since the Sarbecovirus + Hibecovirus clade separated from Nobecovirus, which was estimated to be 1213 years ago. This depletion was not found to be mediated by the GC contents of the genomes. Our results also showed that the depletion of CpG occurred at neutral positions of the genome as well as those under selection. The latter is evident from the progressive reduction in the proportion of arginine amino acid (coded by CpG dinucleotides) in the SARS-CoV-2 lineage over time. The results of this study suggest that shedding CpG nucleotides from their genome is a continuing process in this viral lineage, potentially to escape from their host defense mechanisms.

Long first exons and epigenetic marks distinguish conserved pachytene piRNA clusters from other mammalian genes

In the male germ cells of placental mammals, 26–30-nt-long PIWI-interacting RNAs (piRNAs) emerge when spermatocytes enter the pachytene phase of meiosis. In mice, pachytene piRNAs derive from ~100 discrete autosomal loci that produce canonical RNA polymerase II transcripts. These piRNA clusters bear 5′ caps and 3′ poly(A) tails, and often contain introns that are removed before nuclear export and processing into piRNAs. What marks pachytene piRNA clusters to produce piRNAs, and what confines their expression to the germline? We report that an unusually long first exon (≥ 10 kb) or a long, unspliced transcript correlates with germline-specific transcription and piRNA production. Our integrative analysis of transcriptome, piRNA, and epigenome datasets across multiple species reveals that a long first exon is an evolutionarily conserved feature of pachytene piRNA clusters. Furthermore, a highly methylated promoter, often containing a low or intermediate level of CG dinucleotides, correlates with germline expression and somatic silencing of pachytene piRNA clusters. Pachytene piRNA precursor transcripts bind THOC1 and THOC2, THO complex subunits known to promote transcriptional elongation and mRNA nuclear export. Together, these features may explain why the major sources of pachytene piRNA clusters specifically generate these unique small RNAs in the male germline of placental mammals.

Developmental accumulation of gene body and transposon non-CpG methylation in the zebrafish brain

DNA methylation predominantly occurs at CG dinucleotides in vertebrate genomes; however, non-CG methylation (mCH) is also detectable in vertebrate tissues, most notably in the nervous system. In mammals it is well established that mCH is targeted to CAC trinucleotides by DNMT3A during nervous system development where it is enriched in gene bodies and associated with transcriptional repression. However, the conservation of developmental mCH accumulation and its deposition by DNMT3A is largely unexplored and has yet to be functionally demonstrated in other vertebrates. In this study, by analyzing DNA methylomes and transcriptomes of zebrafish brains, we identified enrichment of mCH at CAC trinucleotides (mCAC) at defined transposon motifs as well as in developmentally downregulated genes associated with developmental and neural functions. We further generated and analyzed DNA methylomes and transcriptomes of developing zebrafish larvae and demonstrated that, like in mammals, mCH accumulates during post-embryonic brain development. Finally, by employing CRISPR/Cas9 technology, we unraveled a conserved role for Dnmt3a enzymes in developmental mCAC deposition. Overall, this work demonstrates the evolutionary conservation of developmental mCH dynamics and highlights the potential of zebrafish as a model to study mCH regulation and function during normal and perturbed development.