Age Prediction Using DNA Methylation Heterogeneity Metrics

Dynamic changes in genomic DNA methylation patterns govern the epigenetic developmental programs and accompany the organism's aging. Epigenetic clock (eAge) algorithms utilize DNA methylation to estimate the age and risk factors for diseases as well as analyze the impact of various interventions. High-throughput bisulfite sequencing methods, such as reduced-representation bisulfite sequencing (RRBS) or whole genome bisulfite sequencing (WGBS), provide an opportunity to identify the genomic regions of disordered or heterogeneous DNA methylation, which might be associated with cell-type heterogeneity, DNA methylation erosion, and allele-specific methylation. We systematically evaluated the applicability of five scores assessing the variability of methylation patterns by evaluating within-sample heterogeneity (WSH) to construct human blood epigenetic clock models using RRBS data. The best performance was demonstrated by the model based on a metric designed to assess DNA methylation erosion with an MAE of 3.686 years. We also trained a prediction model that uses the average methylation level over genomic regions. Although this region-based model was relatively more efficient than the WSH-based model, the latter required the analysis of just a few short genomic regions and, therefore, could be a useful tool to design a reduced epigenetic clock that is analyzed by targeted next-generation sequencing.

Aberrant methylation of placental development genes in chorionic villi of spontaneous abortions with trisomy 16

In humans, aneuploidy is incompatible with the birth of healthy children and mainly leads to the death of embryos in the early stages of development in the first trimester of pregnancy. Trisomy 16 is the most common aneuploidy among spontaneous abortions of the first trimester of pregnancy. However, the mechanisms leading to the death of embryos with trisomy 16 remain insufficiently investigated. One of these potential mechanisms is abnormal placental development, including aberrant remodeling of spiral arteries. Spiral artery remodeling involves the migration of trophoblast cells into the maternal spiral arteries, replacing their endothelium and remodeling to ensure a stable embryonic nutrition and oxygen supply. This is a complex process which depends on many factors from both the embryo and the mother. We analyzed the methylation level of seven genes (ADORA2B, NPR3, PRDM1, PSG2, PHTLH, SV2C, and TICAM2) involved in placental development in the chorionic villi of spontaneous abortions with trisomy 16 (n = 14), compared with spontaneous abortions with a normal karyotype (n = 31) and the control group of induced abortions (n = 10). To obtain sequencing libraries, targeted amplification of individual gene regions using designed oligonucleotide primers for bisulfite-converted DNA was used. The analysis was carried out using targeted bisulfite massive parallel sequencing. In the group of spontaneous abortions with trisomy 16, the level of methylation of the PRDM1 and PSG2 genes was significantly increased compared to induced abortions (p = 0.0004 and p = 0.0015, respectively). In the group of spontaneous abortions, there was no increase in the level of methylation of the PRDM1 and PSG2 genes, but the level of methylation of the ADORA2B gene was significantly increased compared to the induced abortions (p = 0.032). The results obtained indicate the potential mechanisms of the pathogenetic effect of trisomy 16 on the placental development with the participation of the studied genes.

Sperm DNA methylation dynamics after chemotherapy: a longitudinal study of a patient with testicular germ cell tumor treatment

BACKGROUND

An important issue for young men affected by testicular germ cell tumor (TGCT) is how TGCT and its treatment will affect, transiently or permanently, their future reproductive health. Previous studies have reported that xenobiotics can induce changes on human sperm epigenome and have the potential to promote epigenetic alterations in the offspring.

OBJECTIVES

Here, we report the first longitudinal DNA methylation profiling of frozen sperm from a TGCT patient before and up to 2 years after a bleomycin, etoposide, and cisplatin (BEP) chemotherapy.

MATERIALS AND METHODS

A TGCT was diagnosed in a 30-year-old patient. A cryopreservation of spermatozoa was proposed before adjuvant BEP treatment. Semen samples were collected before and after chemotherapy at 6, 9, 12, and 24 months. The DNA methylation status was determined by RRBS to detect DNA differentially methylated regions (DMRs).

RESULTS

The analysis revealed that among the 74 DMRs showing modified methylation status 6 months after therapy, 17 remained altered 24 months after treatment. We next associated DMRs with differentially methylated genes (DMGs), which were subsequently intersected with loci known to be important or expressed during early development.

DISCUSSION AND CONCLUSION

The consequences of the cancer treatment on the sperm epigenome during the recovery periods are topical issues of increasing significance as epigenetic modifications to the paternal genome may have deleterious effects on the offspring. The altered methylated status of these DMGs important for early development might modify their expression pattern and thus affect their function during key stages of embryogenesis, potentially leading to developmental disorders or miscarriages.

Maternal upbringing and selective breeding for voluntary exercise behavior modify patterns of DNA methylation and expression of genes in the mouse brain

Selective breeding has been utilized to study the genetic basis of exercise behavior, but research suggests that epigenetic mechanisms, such as DNA methylation, also contribute to this behavior. In a previous study, we demonstrated that the brains of mice from a genetically selected high runner (HR) line have sex-specific changes in DNA methylation patterns in genes known to be genomically imprinted compared to those from a non-selected control (C) line. Through cross-fostering, we also found that maternal upbringing can modify the DNA methylation patterns of additional genes. Here, we identify an additional set of genes in which DNA methylation patterns and gene expression may be altered by selection for increased wheel-running activity and maternal upbringing. We performed bisulfite sequencing and gene expression assays of 14 genes in the brain and found alterations in DNA methylation and gene expression for Bdnf, Pde4d and Grin2b. Decreases in Bdnf methylation correlated with significant increases in Bdnf gene expression in the hippocampus of HR compared to C mice. Cross-fostering also influenced the DNA methylation patterns for Pde4d in the cortex and Grin2b in the hippocampus, with associated changes in gene expression. We also found that the DNA methylation patterns for Atrx and Oxtr in the cortex and Atrx and Bdnf in the hippocampus were further modified by sex. Together with our previous study, these results suggest that DNA methylation and the resulting change in gene expression may interact with early-life influences to shape adult exercise behavior.

PSTVd infection in Nicotiana benthamiana plants has a minor yet detectable effect on CG methylation

Viroids are small circular RNAs infecting a wide range of plants. They do not code for any protein or peptide and therefore rely on their structure for their biological cycle. Observed phenotypes of viroid infected plants are thought to occur through changes at the transcriptional/translational level of the host. A mechanism involved in such changes is RNA-directed DNA methylation (RdDM). Till today, there are contradictory works about viroids interference of RdDM. In this study, we investigated the epigenetic effect of viroid infection in Nicotiana benthamiana plants. Using potato spindle tuber viroid (PSTVd) as the triggering pathogen and via bioinformatic analyses, we identified endogenous gene promoters and transposable elements targeted by 24 nt host siRNAs that differentially accumulated in PSTVd-infected and healthy plants. The methylation status of these targets was evaluated following digestion with methylation-sensitive restriction enzymes coupled with PCR amplification, and bisulfite sequencing. In addition, we used Methylation Sensitive Amplification Polymorphism (MSAP) followed by sequencing (MSAP-seq) to study genomic DNA methylation of 5-methylcytosine (5mC) in CG sites upon viroid infection. In this study we identified a limited number of target loci differentially methylated upon PSTVd infection. These results enhance our understanding of the epigenetic host changes as a result of pospiviroid infection.

Evaluating Stacked Methylation Markers for Blood-Based Multicancer Detection

The ability to detect several types of cancer using a non-invasive, blood-based test holds the potential to revolutionize oncology screening. We mined tumor methylation array data from the Cancer Genome Atlas (TCGA) covering 14 cancer types and identified two novel, broadly-occurring methylation markers at TLX1 and GALR1. To evaluate their performance as a generalized blood-based screening approach, along with our previously reported methylation biomarker, ZNF154, we rigorously assessed each marker individually or combined. Utilizing TCGA methylation data and applying logistic regression models within each individual cancer type, we found that the three-marker combination significantly increased the average area under the ROC curve (AUC) across the 14 tumor types compared to single markers (p = 1.158 × 10-10; Friedman test). Furthermore, we simulated dilutions of tumor DNA into healthy blood cell DNA and demonstrated increased AUC of combined markers across all dilution levels. Finally, we evaluated assay performance in bisulfite sequenced DNA from patient tumors and plasma, including early-stage samples. When combining all three markers, the assay correctly identified nine out of nine lung cancer plasma samples. In patient plasma from hepatocellular carcinoma, ZNF154 alone yielded the highest combined sensitivity and specificity values averaging 68% and 72%, whereas multiple markers could achieve higher sensitivity or specificity, but not both. Altogether, this study presents a comprehensive pipeline for the identification, testing, and validation of multi-cancer methylation biomarkers with a considerable potential for detecting a broad range of cancer types in patient blood samples.

An ecologist's guide for studying DNA methylation variation in wild vertebrates

The field of molecular biology is advancing fast with new powerful technologies, sequencing methods and analysis software being developed constantly. Commonly used tools originally developed for research on humans and model species are now regularly used in ecological and evolutionary research. There is also a growing interest in the causes and consequences of epigenetic variation in natural populations. Studying ecological epigenetics is currently challenging, especially for vertebrate systems, because of the required technical expertise, complications with analyses and interpretation, and limitations in acquiring sufficiently high sample sizes. Importantly, neglecting the limitations of the experimental setup, technology and analyses may affect the reliability and reproducibility, and the extent to which unbiased conclusions can be drawn from these studies. Here, we provide a practical guide for researchers aiming to study DNA methylation variation in wild vertebrates. We review the technical aspects of epigenetic research, concentrating on DNA methylation using bisulfite sequencing, discuss the limitations and possible pitfalls, and how to overcome them through rigid and reproducible data analysis. This review provides a solid foundation for the proper design of epigenetic studies, a clear roadmap on the best practices for correct data analysis and a realistic view on the limitations for studying ecological epigenetics in vertebrates. This review will help researchers studying the ecological and evolutionary implications of epigenetic variation in wild populations.

Intersection clock reveals a rejuvenation event during human embryogenesis

Recent research revealed a rejuvenation event during early development of mice. Here, by examining epigenetic age dynamics of human embryogenesis, we tested whether a similar event exists in humans. For this purpose, we developed an epigenetic clock method, the intersection clock, that utilizes bisulfite sequencing in a way that maximizes the use of informative CpG sites with no missing clock CpG sites in test samples and applied it to human embryo development data. We observed no changes in the predicted epigenetic age between cleavage stage and blastocyst stage embryos; however, a significant decrease was observed between blastocysts and cells representing the epiblast. Additionally, by applying the intersection clock to datasets spanning pre and postimplantation, we found no significant change in the epigenetic age during preimplantation stages; however, the epigenetic age of postimplantation samples was lower compared to the preimplantation stages. We further investigated the epigenetic age of primed (representing early postimplantation) and naïve (representing preimplantation) pluripotent stem cells and observed that in all cases the epigenetic age of primed cells was significantly lower than that of naïve cells. Together, our data suggest that human embryos are rejuvenated during early embryogenesis. Hence, the rejuvenation event is conserved between the mouse and human, and it occurs around the gastrulation stage in both species. Beyond this advance, the intersection clock opens the way for other epigenetic age studies based on human bisulfite sequencing datasets as opposed to methylation arrays.

DNA Methylation in the Fields of Prenatal Diagnosis and Early Detection of Cancers

The central objective of the metamorphosis of discovery science into biomedical applications is to serve the purpose of patients and curtail the global disease burden. The journey from the discovery of DNA methylation (DNAm) as a biological process to its emergence as a diagnostic tool is one of the finest examples of such metamorphosis and has taken nearly a century. Particularly in the last decade, the application of DNA methylation studies in the clinic has been standardized more than ever before, with great potential to diagnose a multitude of diseases that are associated with a burgeoning number of genes with this epigenetic alteration. Fetal DNAm detection is becoming useful for noninvasive prenatal testing, whereas, in very preterm infants, DNAm is also shown to be a potential biological indicator of prenatal risk factors. In the context of cancer, liquid biopsy-based DNA-methylation profiling is offering valuable epigenetic biomarkers for noninvasive early-stage diagnosis. In this review, we focus on the applications of DNA methylation in prenatal diagnosis for delivering timely therapy before or after birth and in detecting early-stage cancers for better clinical outcomes. Furthermore, we also provide an up-to-date commercial landscape of DNAm biomarkers for cancer detection and screening of cancers of unknown origin.

Nucleic acid biomarkers of immune response and cell and tissue damage in children with COVID-19 and MIS-C

Differential host responses in coronavirus disease 2019 (COVID-19) and multisystem inflammatory syndrome in children (MIS-C) remain poorly characterized. Here, we use next-generation sequencing to longitudinally analyze blood samples from pediatric patients with COVID-19 or MIS-C across three hospitals. Profiling of plasma cell-free nucleic acids uncovers distinct signatures of cell injury and death between COVID-19 and MIS-C, with increased multiorgan involvement in MIS-C encompassing diverse cell types, including endothelial and neuronal cells, and an enrichment of pyroptosis-related genes. Whole-blood RNA profiling reveals upregulation of similar pro-inflammatory pathways in COVID-19 and MIS-C but also MIS-C-specific downregulation of T cell-associated pathways. Profiling of plasma cell-free RNA and whole-blood RNA in paired samples yields different but complementary signatures for each disease state. Our work provides a systems-level view of immune responses and tissue damage in COVID-19 and MIS-C and informs future development of new disease biomarkers.

No evidence for epitranscriptomic m5C modification of SARS-CoV-2, HIV and MLV viral RNA

The addition of chemical groups to cellular RNA to modulate RNA fate and/or function is summarized under the term epitranscriptomic modification. More than 170 different modifications have been identified on cellular RNA, such as tRNA, rRNA and, to a lesser extent, on other RNA types. Recently, epitranscriptomic modification of viral RNA has received considerable attention as a possible additional mechanism regulating virus infection and replication. N6-methyladenosine (m6A) and C5-methylcytosine (m5C) have been most broadly studied in different RNA viruses. Various studies, however, reported varying results with regard to number and extent of the modification. Here we investigated the m5C methylome of SARS-CoV-2, and we reexamined reported m5C sites in HIV and MLV. Using a rigorous bisulfite-sequencing protocol and stringent data analysis, we found no evidence for the presence of m5C in these viruses. The data emphasize the necessity for optimizing experimental conditions and bioinformatic data analysis.

Exploring the crop epigenome: a comparison of DNA methylation profiling techniques

Epigenetic modifications play a vital role in the preservation of genome integrity and in the regulation of gene expression. DNA methylation, one of the key mechanisms of epigenetic control, impacts growth, development, stress response and adaptability of all organisms, including plants. The detection of DNA methylation marks is crucial for understanding the mechanisms underlying these processes and for developing strategies to improve productivity and stress resistance of crop plants. There are different methods for detecting plant DNA methylation, such as bisulfite sequencing, methylation-sensitive amplified polymorphism, genome-wide DNA methylation analysis, methylated DNA immunoprecipitation sequencing, reduced representation bisulfite sequencing, MS and immuno-based techniques. These profiling approaches vary in many aspects, including DNA input, resolution, genomic region coverage, and bioinformatics analysis. Selecting an appropriate methylation screening approach requires an understanding of all these techniques. This review provides an overview of DNA methylation profiling methods in crop plants, along with comparisons of the efficacy of these techniques between model and crop plants. The strengths and limitations of each methodological approach are outlined, and the importance of considering both technical and biological factors are highlighted. Additionally, methods for modulating DNA methylation in model and crop species are presented. Overall, this review will assist scientists in making informed decisions when selecting an appropriate DNA methylation profiling method.

Quantitative analysis of MGMT promoter methylation in glioblastoma suggests nonlinear prognostic effect

Background

Epigenetic inhibition of the O6-methylguanine-DNA-methyltransferase (MGMT) gene has emerged as a clinically relevant prognostic marker in glioblastoma (GBM). Methylation of the MGMT promoter has been shown to increase chemotherapy efficacy. While traditionally reported as a binary marker, recent methodological advancements have led to quantitative methods of measuring promoter methylation, providing clearer insight into its functional relationship with survival.

Methods

A CLIA assay and bisulfite sequencing was utilized to develop a quantitative, 17-point, MGMT promoter methylation index. GBMs of 240 newly diagnosed patients were sequenced and risk for mortality was assessed. Nonlinearities were captured by fitting splines to Cox proportional hazard models and plotting smoothed residuals. Covariates included age, Karnofsky performance status, IDH1 mutation, and extent of resection.

Results

Median follow-up time and progression-free survival were 16 and 9 months, respectively. A total of 176 subjects experienced death. A one-unit increase in promoter CpG methylation resulted in a 4% reduction in hazard (95% CI 0.93-0.99, P < .005). GBM patients with low levels of promoter methylation (1-6 CpG sites) fared markedly worse (HR = 1.62, 95% CI 1.03-2.54, P < .036) than individuals who were unmethylated. Subjects with medium levels of promoter methylation (7-12 sites) had the greatest reduction in hazard (HR = 0.48, 95% CI 0.29-0.80, P < .004), followed by individuals in the highest promoter methylation tertile (HR = 0.62, 95% CI 0.40-0.97, P < .035).

Conclusions

Our findings suggest that the relationship between the extent of MGMT promoter methylation and survival in GBM may be nonlinear. These findings challenge the current understanding of MGMT and underlines the clinical importance of determining its prognostic utility. Potential limitations include censoring, sample size, and extraneous mutations.

Genome assemblies of Vigna reflexo-pilosa (créole bean) and its progenitors, Vigna hirtella and Vigna trinervia, revealed homoeolog expression bias and expression-level dominance in the allotetraploid

Vigna reflexo-pilosa (créole bean) is a wild legume belonging to the subgenus Ceratoropis and is widely distributed in Asia. Créole bean is the only tetraploid species in the genus Vigna, and it has been shown to derive from the hybridization of Vigna hirtella and Vigna trinervia. In this study, we combined the long-read PacBio technology with the chromatin contact mapping (Hi-C) technique to obtain a chromosome-level assembly of V. reflexo-pilosa. The final assembly contained 998,724,903 bases with an N50 length of 42,545,650 bases. Our gene prediction recovered 99.4% of the highly conserved orthologs based on the BUSCO analysis. To investigate homoeolog expression bias and expression level dominance in the tetraploid, we also sequenced and assembled the genomes of its progenitors. Overall, the majority of the homoeolog pairs (72.9%) displayed no expression bias, and among those that exhibited biased expression, 16.3% showed unbalanced homoeolog expression bias toward the V. trinervia subgenome. Moreover, 41.2% and 36.2% of the expressed gene pairs exhibited transgressive expression and expression level dominance, respectively. Interestingly, the genome-wide expression level dominance in the tetraploid was biased toward the V. trinervia subgenome. The analysis of methylation patterns also revealed that the average methylation levels in coding regions were higher in the V. hirtella subgenome than those in the V. trinervia subgenome. The genomic/transcriptomic resources for these three species are useful not only for the development of elite cultivars in Vigna breeding programs but also to researchers studying comparative genomics and investigating genomic/epigenomic changes following polyploid events.

Identification of differential methylome signatures of white pigmented skin patches in Nili Ravi buffalo of India

The DNA methylation events mark a major epigenetic change in the genome, reflecting non-genetic disease developments and varied phenotypes. The water buffalo is a dairy production animal with wide agro-climatic distribution in India. Breed-wise the coat color of water buffalo varies from ash-gray to jet black. A typical pigmentation pattern is found in one of the breeds of North India, Nili Ravi, with variedly distributed white patches. The DNA methylation pattern could potentially reveal the epigenetic factors responsible for the pigmentation patterns. To address this question, the DNA isolated from the skin tissues of Nili Ravi with varied white pigmentation and black Murrah buffaloes was subjected to reduced representation bisulfite sequencing. DNA methylation analysis revealed, 68.44%, 63.39%, and 47.94% of the promoter regions were hypermethylated in Nili Ravi over-white versus Murrah, Nili Ravi under-white versus Murrah, and Nili Ravi under-white versus Nili Ravi over-white, respectively. Major genes identified to be differentially methylated among over-white and under-white skin tissues in Nili Ravi included TBX2, SNAI2, HERC2, and CITED1. Overall the results have indicated differential methylation patterns to be potentially involved in hyper or hypopigmentation in Nili Ravi and Murrah buffaloes.

Allele-specific cis-regulatory methylation of the gene for vasoactive intestinal peptide in white-throated sparrows

White-throated sparrows (Zonotrichia albicollis) offer a unique opportunity to connect genotype with behavioral phenotype. In this species, a rearrangement of the second chromosome is linked with territorial aggression; birds with a copy of this "supergene" rearrangement are more aggressive than those without it. The supergene has captured the gene VIP, which encodes vasoactive intestinal peptide, a neuromodulator that drives aggression in other songbirds. In white-throated sparrows, VIP expression is higher in the anterior hypothalamus of birds with the supergene than those without it, and expression of VIP in this region predicts the level of territorial aggression regardless of genotype. Here, we aimed to identify epigenetic mechanisms that could contribute to differential expression of VIP both in breeding adults, which exhibit morph differences in territorial aggression, and in nestlings, before territorial behavior develops. We extracted and bisulfite-converted DNA from samples of the hypothalamus in wild-caught adults and nestlings and used high-throughput sequencing to measure DNA methylation of a region upstream of the VIP start site. We found that the allele inside the supergene was less methylated than the alternative allele in both adults and nestlings. The differential methylation was attributed primarily to CpG sites that were shared between the alleles, not to polymorphic sites, which suggests that epigenetic regulation is occurring independently of the genetic differentiation within the supergene. This work represents an initial step toward understanding how epigenetic differentiation inside chromosomal inversions leads to the development of alternative behavioral phenotypes.

Mitotically heritable epigenetic modifications of CmMYB6 control anthocyanin biosynthesis in chrysanthemum

Flower color, which is determined by various chemical pigments, is a vital trait for ornamental plants, in which anthocyanin is a major component. However, the epigenetic regulation of anthocyanin biosynthesis remains poorly understood. During chrysanthemum cultivation, we found a heterochromatic chrysanthemum accession (YP) whose progeny generated by asexual reproduction contained both yellow-flowered (YP-Y) and pink-flowered (YP-P) plants. In this study, we aimed to elucidate the epigenetic mechanisms of different flower colors in the YP plant progeny. Metabolome and transcriptome analyses revealed that the difference in flower color between YP-Y and YP-P was caused by expression variation of the anthocyanin biosynthesis gene CmMYB6. Bisulfite sequencing revealed that methylation at the CmMYB6 promoter, especially in the CHH context, was higher in YP-Y than YP-P. After demethylation of the CmMYB6 promoter using the dCas9-TET1cd system, the flower color returned from yellow to pink. Furthermore, the methylation status of the CmMYB6 promoter was higher in YP-Y over three consecutive generations, indicating that this methylation status was heritable mitotically. Finally, investigation of other chrysanthemum cultivars showed that the methylation of CmMYB6 decreased gradually with the increase in anthocyanin content. These results lay an epigenetic foundation for the improvement of flower color in horticultural plants.

Multimodal analysis of ctDNA methylation and fragmentomic profiles enhances detection of nonmetastatic colorectal cancer

Aims: Early detection of colorectal cancer (CRC) provides substantially better survival rates. This study aimed to develop a blood-based screening assay named SPOT-MAS ('screen for the presence of tumor by DNA methylation and size') for early CRC detection with high accuracy. Methods: Plasma cell-free DNA samples from 159 patients with nonmetastatic CRC and 158 healthy controls were simultaneously analyzed for fragment length and methylation profiles. We then employed a deep neural network with fragment length and methylation signatures to build a classification model. Results: The model achieved an area under the curve of 0.989 and a sensitivity of 96.8% at 97% specificity in detecting CRC. External validation of our model showed comparable performance, with an area under the curve of 0.96. Conclusion: SPOT-MAS based on integration of cancer-specific methylation and fragmentomic signatures could provide high accuracy for early-stage CRC detection.

Targeted Long-Read Bisulfite Sequencing Identifies Differences in the TERT Promoter Methylation Profiles between TERT Wild-Type and TERT Mutant Cancer Cells

Background: TERT promoter methylation, located several hundred base pairs upstream of the transcriptional start site, is cancer specific and correlates with increased TERT mRNA expression and poorer patient outcome. Promoter methylation, however, is not mutually exclusive to TERT activating genetic alterations, as predicted for functionally redundant mechanisms. To annotate the altered patterns of TERT promoter methylation and their relationship with gene expression, we applied a Pacific Biosciences-based, long-read, bisulfite-sequencing technology and compared the differences in the methylation marks between wild-type and mutant cancers in an allele-specific manner. Results: We cataloged TERT genetic alterations (i.e., promoter point mutations or structural variations), allele-specific promoter methylation patterns, and allele-specific expression levels in a cohort of 54 cancer cell lines. In heterozygous mutant cell lines, the mutant alleles were significantly less methylated than their silent, mutation-free alleles (p < 0.05). In wild-type cell lines, by contrast, both epialleles were equally methylated to high levels at the TERT distal promoter, but differentially methylated in the proximal regions. ChIP analysis showed that epialleles with the hypomethylated proximal and core promoter were enriched in the active histone mark H3K4me2/3, whereas epialleles that were methylated in those regions were enriched in the repressive histone mark H3K27me3. Decitabine therapy induced biallelic expression in the wild-type cancer cells, whereas the mutant cell lines were unaffected. Conclusions: Long-read bisulfite sequencing analysis revealed differences in the methylation profiles and responses to demethylating agents between TERT wild-type and genetically altered cancer cell lines. The causal relation between TERT promoter methylation and gene expression remains to be established.

Whole-genome DNA methylation patterns of Oryza sativa (L.) and Oryza glumaepatula (Steud) genotypes associated with aluminum response

Epigenetic mechanisms in crops have emerged as a fundamental factor in plant adaptation and acclimation to biotic and abiotic stresses. Among described epigenetic mechanisms, DNA methylation has been defined as the most studied epigenetic modification involved in several developmental processes. It has been shown that contrasting methylation marks are associated with gene expression variations between cultivated and wild crop species. In this study, we analyzed single-base resolution methylome maps for Oryza sativa (a cultivated species) and Oryza glumaepatula (a wild species) genotypes grown under control conditions. Our results showed that overall, genome-wide methylation profiles are mainly conserved between both species, nevertheless, there are several differentially methylated regions with species-specific methylation patterns. In addition, we analyzed the association of identified DNA methylation marks in relation with Aluminum-tolerance levels of studied genotypes. We found several differentially methylated regions (DMRs) and DMR-associated genes (DAGs) that are linked with Al tolerance. Some of these DAGs have been previously reported as differentially expressed under Al exposure in O. sativa. Complementarily a Transposable Elements (TE) analysis revealed that specific aluminum related genes have associated-TEs potentially regulated by DNA methylation. Interestingly, the DMRs and DAGs between Al-tolerant and susceptible genotypes were different between O. sativa and O. glumaepatula, suggesting that methylation patterns related to Al responses are unique for each rice species. Our findings provide novel insights into DNA methylation patterns in wild and cultivated rice genotypes and their possible role in the regulation of plant stress responses.

Structure and Methylation of 35S rDNA in Allopolyploids Anemone multifida (2n = 4x = 32, BBDD) and Anemone baldensis (2n = 6x = 48, AABBDD) and Their Parental Species Show Evidence of Nucleolar Dominance

Transcriptional silencing of 35S rDNA loci inherited from one parental species is occurring relatively frequently in allopolyploids. However, molecular mechanisms by which it is selected for transcriptional silencing remain unclear. We applied NGS, silver staining and bisulfite sequencing to study the structure, expression and methylation landscape of 35S rDNA in two allopolyploids of common origin, allotetraploid Anemone multifida (2n = 4x = 32, genome composition BBDD) and allohexaploid A. baldensis (2n = 6x = 48, AABBDD), and their genome donors, A. sylvestris (2n = 16, AA), A. cylindrica (2n = 16, BB) and A. parviflora (2n = 16, DD). The size of the recovered 35S rDNA units varied from 10,489 bp in A. cylindrica to 12,084 bp in A. sylvestris. Anemone showed an organization typical of most ribosomal 35S rDNA composed of NTS, ETS, rRNA genes, TTS and TIS with structural features of plant IGS sequences and all functional elements needed for rRNA gene activity. The NTS was more variable than the ETS and consisted of SRs which are highly variable among Anemone. Five to six CpG-rich islands were found within the ETS. CpG island located adjacent to the transcription initiation site (TIS) was highly variable regarding the sequence size and methylation level and exhibited in most of the species lower levels of methylation than CpG islands located adjacent to the 18S rRNA gene. Our results uncover hypomethylation of A. sylvestris- and A. parviflora-derived 35S rDNA units in allopolyploids A. multifida and A. baldensis. Hypomethylation of A. parviflora-derived 35S rDNA was more prominent in A. baldensis than in A. multifida. We showed that A. baldensis underwent coupled A. sylvestris-derived 35S rDNA array expansion and A. parviflora-derived 35S rDNA copy number decrease that was accompanied by lower methylation level of A. sylvestris-derived 35S rDNA units in comparison to A. parviflora-derived 35S rDNA units. These observations suggest that in A. baldensis nucleolar dominance is directed toward A. sylvestris-derived chromosomes. This work broadens our current knowledge of the 35S rDNA organization in Anemone and provides evidence of the progenitor-specific 35S rDNA methylation in nucleolar dominance.

Identification of tissue of origin in cancer of unknown primary using a targeted bisulfite sequencing panel

Aim: To construct a targeted bisulfite sequencing panel predicting origin of cancer of unknown primary. Methods: A bisulfite sequencing panel targeting 2793 tissue-specific markers was performed in 100 clinical samples. Results: The authors' prediction model showed 0.85 accuracy for the 'first-ranked' tissue type and 0.93 accuracy for the 'second-ranked' tissue type using 2793 tissue-specific markers and 0.84 accuracy for the 'first-ranked' tissue type and 0.92 accuracy for the 'second-ranked' tissue type when the number of tissue-specific markers was reduced to 514. Conclusion: Targeted bisulfite sequencing is a useful method for predicting the tissue of origin in patients with cancer of unknown primary.

Methylated ccfDNA from plasma biomarkers of Alzheimer's disease using targeted bisulfite sequencing

Aim: Noninvasive biomarkers such as methylated ccfDNA from plasma could help to support the diagnosis of Alzheimer's disease (AD). Methods: A targeted sequencing protocol was developed to identify candidate biomarkers of AD in methylated ccfDNA extracted from plasma. Results: The authors identified differentially methylated CpGs, regions of which were the same as those identified in previous AD studies. Specifically, a differentially methylated CpG of the LHX2 gene previously identified in a plasma study of AD was replicated in the study. The MBP and DUSP22 regions have been identified in other brain studies of AD and in the authors' study. Conclusion: Although these biomarkers must be validated in other cohorts, methylated ccfDNA could be a relevant noninvasive biomarker in AD.

epiGBS2: Improvements and evaluation of highly multiplexed, epiGBS-based reduced representation bisulfite sequencing

Several reduced‐representation bisulfite sequencing methods have been developed in recent years to determine cytosine methylation de novo in nonmodel species. Here, we present epiGBS2, a laboratory protocol based on epiGBS with a revised and user‐friendly bioinformatics pipeline for a wide range of species with or without a reference genome. epiGBS2 is cost‐ and time‐efficient and the computational workflow is designed in a user‐friendly and reproducible manner. The library protocol allows a flexible choice of restriction enzymes and a double digest. The bioinformatics pipeline was integrated in the snakemake workflow management system, which makes the pipeline easy to execute and modular, and parameter settings for important computational steps flexible. We implemented bismark for alignment and methylation analysis and we preprocessed alignment files by double masking to enable single nucleotide polymorphism calling with freebayes (epifreebayes). The performance of several critical steps in epiGBS2 was evaluated against baseline data sets from Arabidopsis thaliana and great tit (Parus major), which confirmed its overall good performance. We provide a detailed description of the laboratory protocol and an extensive manual of the bioinformatics pipeline, which is publicly accessible on github (https://github.com/nioo‐knaw/epiGBS2) and zenodo (https://doi.org/10.5281/zenodo.4764652).

Quantitative proteomics revealed new functions of ALKBH4

ALKBH4 is a versatile demethylase capable of catalyzing the demethylation of monomethylated lysine-84 on actin and N⁶ -methyladenine in DNA. In this study, we conducted a quantitative proteomic experiment to reveal the altered expression of proteins in HEK293T cells upon genetic ablation of ALKBH4. Our results showed markedly diminished levels of GSTP1 and HSPB1 proteins in ALKBH4-depleted cells, which emanate from an augmented expression level of DNA (cytosine-5)-methyltransferase 1 (DNMT1) and the ensuing elevated cytosine methylation in the promoter regions of GSTP1 and HSPB1 genes. Together, our results revealed a role of ALKBH4 in modulating DNA cytosine methylation through regulating the expression level of DNMT1 protein. This article is protected by copyright. All rights reserved.

CUT&Tag-BS for simultaneous profiling of histone modification and DNA methylation with high efficiency and low cost

It remains a challenge to decipher the complex relationship between DNA methylation, histone modification, and the underlying DNA sequence with limited input material. Here, we developed an efficient, low-input, and low-cost method for the simultaneous profiling of genomic localization of histone modification and methylation status of the underlying DNA at single-base resolution from the same cells in a single experiment by integrating cleavage under targets and tagmentation (CUT&Tag) with tagmentation-based bisulfite sequencing (CUT&Tag-BS). We demonstrated the validity of our method using representative histone modifications of euchromatin and constitutive and facultative heterochromatin (H3K4me1, H3K9me3, and H3K27me3, respectively). Similar histone modification enrichment patterns were observed in CUT&Tag-BS compared with non-bisulfite-treated control, and H3K4me1-marked regions were found to mostly be CpG poor, lack methylation concordance, and exhibit prevalent DNA methylation heterogeneity among mouse embryonic stem cells (mESCs). We anticipate that CUT&Tag-BS will be widely applied to directly address the genomic relationship between DNA methylation and histone modification, especially in low-input scenarios with precious biological samples.

Location-Dependent DNA Methylation Signatures in a Clonal Invasive Crayfish

DNA methylation is an important epigenetic modification that has been repeatedly implied in organismal adaptation. However, many previous studies that have linked DNA methylation patterns to environmental parameters have been limited by confounding factors, such as cell-type heterogeneity and genetic variation. In this study, we analyzed DNA methylation variation in marbled crayfish, a clonal and invasive freshwater crayfish that is characterized by a largely tissue-invariant methylome and negligible genetic variation. Using a capture-based subgenome bisulfite sequencing approach that covers a small, variably methylated portion of the marbled crayfish genome, we identified specific and highly localized DNA methylation signatures for specimens from geographically and ecologically distinct wild populations. These results were replicated both biologically and technically by re-sampling at different time points and by using independent methodology. Finally, we show specific methylation signatures for laboratory animals and for laboratory animals that were reared at a lower temperature. Our results thus demonstrate the existence of context-dependent DNA methylation signatures in a clonal animal.

Invertebrate methylomes provide insight into mechanisms of environmental tolerance and reveal methodological biases

There is a growing focus on the role of DNA methylation in the ability of marine invertebrates to rapidly respond to changing environmental factors and anthropogenic impacts. However, genome-wide DNA methylation studies in nonmodel organisms are currently hampered by a limited understanding of methodological biases. Here, we compare three methods for quantifying DNA methylation at single base-pair resolution-whole genome bisulfite sequencing (WGBS), reduced representation bisulfite sequencing (RRBS), and methyl-CpG binding domain bisulfite sequencing (MBDBS)-using multiple individuals from two reef-building coral species with contrasting environmental sensitivity. All methods reveal substantially greater methylation in Montipora capitata (11.4%) than the more sensitive Pocillopora acuta (2.9%). The majority of CpG methylation in both species occurs in gene bodies and flanking regions. In both species, MBDBS has the greatest capacity for detecting CpGs in coding regions at our sequencing depth, but MBDBS may be influenced by intrasample methylation heterogeneity. RRBS yields robust information for specific loci albeit without enrichment of any particular genome feature and with significantly reduced genome coverage. Relative genome size strongly influences the number and location of CpGs detected by each method when sequencing depth is limited, illuminating nuances in cross-species comparisons. As genome-wide methylation differences, supported by data across bisulfite sequencing methods, may contribute to environmental sensitivity phenotypes in critical marine invertebrate taxa, these data provide a genomic resource for investigating the functional role of DNA methylation in environmental tolerance.

DNA Methylation Level Changes in Transgenic Chinese Cabbage (Brassica rapa ssp. pekinensis) Plants and Their Effects on Corresponding Gene Expression Patterns

Plant tissue culture is an in vitro technique used to manipulate cells, tissues, or organs, and plays an important role in genetic transformation. However, plants cultured in vitro often exhibit unintended genetic and epigenetic variations. Since it is important to secure the stability of endogenous and exogenous gene expressions in transgenic plants, it is preferable to avoid the occurrence of such variations. In this study, we focused on epigenetic variations, exclusively on methylation level changes of DNA, in transgenic Chinese cabbage (Brassica rapa ssp. pekinensis) plants. To detect these methylation level changes of DNA, bisulfite sequencing was performed and the obtained sequences were compared with the ‘CT001’ reference genome. Differentially methylated regions (DMRs) of DNA between the non-transgenic and transgenic lines were detected by bisulfite sequencing, and ten DMRs located in exonic regions were identified. The regions with methylation variations that were inherited and consistently maintained in the next generation lines were selected and validated. We also analyzed the relationship between methylation status and expression levels of transformant-conserved DMR (TCD) genes by quantitative reverse transcription-PCR. These results suggested that the changes in methylation levels of these DMRs might have been related to the plant transformation process, affecting subsequent gene expression. Our findings can be used in fundamental research on methylation variations in transgenic plants and suggest that these variations affect the expression of the associated genes.

Epigenetic Forensics for Suspect Identification and Age Prediction

Background: Genetic testing at crime scenes is an instrumental molecular technique to identify or eliminate suspects, as well as to overturn wrongful convictions. Yet, genotyping alone cannot reveal the age of a sample, which could help advance the utility of crime scene samples for suspect identification. The distribution of cytosine methylation within a DNA sample can be leveraged to determine the epigenetic age of someone's blood. Methodology: We sought to demonstrate the ability of DNA methylation markers to accurately discern the age of blood spots from an actual crime scene, a "mock" crime scene, and also from a tube of blood stored in ethylenediaminetetraacetic acid for >20 years. This was achieved by quantifying methylation within known age-associated genetic loci across each DNA sample. We observed a strong linear coefficient (0.91) and high overall correlation (R ² = 0.963) between the known age of a sample and the predicted age. Conclusion: We show that novel methods for targeted methylation and low-input whole-genome bisulfite sequencing can enable a novel and improved forensic profile of a crime scene that discerns not only who was present at the crime, but also their age. Finally, we use this model to discern the age and provenance of a blood sample that was used in a criminal investigation.

Divergence among rice cultivars reveals roles for transposition and epimutation in ongoing evolution of genomic imprinting

Parent-of-origin-dependent gene expression in mammals and flowering plants results from differing chromatin imprints (genomic imprinting) between maternally and paternally inherited alleles. Imprinted gene expression in the endosperm of seeds is associated with localized hypomethylation of maternally but not paternally inherited DNA, with certain small RNAs also displaying parent-of-origin-specific expression. To understand the evolution of imprinting mechanisms in Oryza sativa (rice), we analyzed imprinting divergence among four cultivars that span both japonica and indica subspecies: Nipponbare, Kitaake, 93-11, and IR64. Most imprinted genes are imprinted across cultivars and enriched for functions in chromatin and transcriptional regulation, development, and signaling. However, 4 to 11% of imprinted genes display divergent imprinting. Analyses of DNA methylation and small RNAs revealed that endosperm-specific 24-nt small RNA-producing loci show weak RNA-directed DNA methylation, frequently overlap genes, and are imprinted four times more often than genes. However, imprinting divergence most often correlated with local DNA methylation epimutations (9 of 17 assessable loci), which were largely stable within subspecies. Small insertion/deletion events and transposable element insertions accompanied 4 of the 9 locally epimutated loci and associated with imprinting divergence at another 4 of the remaining 8 loci. Correlating epigenetic and genetic variation occurred at key regulatory regions-the promoter and transcription start site of maternally biased genes, and the promoter and gene body of paternally biased genes. Our results reinforce models for the role of maternal-specific DNA hypomethylation in imprinting of both maternally and paternally biased genes, and highlight the role of transposition and epimutation in rice imprinting evolution.

ZFP57 dictates allelic expression switch of target imprinted genes

ZFP57 is a master regulator of genomic imprinting. It has both maternal and zygotic functions that are partially redundant in maintaining DNA methylation at some imprinting control regions (ICRs). In this study, we found that DNA methylation was lost at most known ICRs in Zfp57 mutant embryos. Furthermore, loss of ZFP57 caused loss of parent-of-origin-dependent monoallelic expression of the target imprinted genes. The allelic expression switch occurred in the ZFP57 target imprinted genes upon loss of differential DNA methylation at the ICRs in Zfp57 mutant embryos. Specifically, upon loss of ZFP57, the alleles of the imprinted genes located on the same chromosome with the originally methylated ICR switched their expression to mimic their counterparts on the other chromosome with unmethylated ICR. Consistent with our previous study, ZFP57 could regulate the NOTCH signaling pathway in mouse embryos by impacting allelic expression of a few regulators in the NOTCH pathway. In addition, the imprinted Dlk1 gene that has been implicated in the NOTCH pathway was significantly down-regulated in Zfp57 mutant embryos. Our allelic expression switch models apply to the examined target imprinted genes controlled by either maternally or paternally methylated ICRs. Our results support the view that ZFP57 controls imprinted expression of its target imprinted genes primarily through maintaining differential DNA methylation at the ICRs.

Viroids as a Tool to Study RNA-Directed DNA Methylation in Plants

Viroids are plant pathogenic, circular, non-coding, single-stranded RNAs (ssRNAs). Members of the Pospiviroidae family replicate in the nucleus of plant cells through double-stranded RNA (dsRNA) intermediates, thus triggering the host’s RNA interference (RNAi) machinery. In plants, the two RNAi pillars are Post-Transcriptional Gene Silencing (PTGS) and RNA-directed DNA Methylation (RdDM), and the latter has the potential to trigger Transcriptional Gene Silencing (TGS). Over the last three decades, the employment of viroid-based systems has immensely contributed to our understanding of both of these RNAi facets. In this review, we highlight the role of Pospiviroidae in the discovery of RdDM, expound the gradual elucidation through the years of the diverse array of RdDM’s mechanistic details and propose a revised RdDM model based on the cumulative amount of evidence from viroid and non-viroid systems.

A New Approach to Identify the Methylation Sites in the Control Region of Mitochondrial DNA

Mitochondrial DNA (mtDNA) methylation has the potential to be used as a biomarker of human development or disease. However, mtDNA methylation procedures are costly and time-consuming. Therefore, we developed a new approach based on an RT-PCR assay for the base site identification of methylated cytosine in the control region of mtDNA through a simple, fast, specific, and low-cost strategy. Total DNA was purified, and methylation was determined by RT-PCR bisulfite sequencing. This procedure included the DNA purification, bisulfite treatment and RT-PCR amplification of the control region divided into three subregions with specific primers. Sequences obtained with and without the bisulfite treatment were compared to identify the methylated cytosine dinucleotides. Furthermore, the efficiency of C to U conversion of cytosines was assessed by including a negative control. Interestingly, mtDNA methylation was observed mainly within non-Cphosphate- G (non-CpG) dinucleotides and mostly in the regions containing regulatory elements, such as OH or CSBI, CSBII, and CSBIII. This new approach will promote the generation of new information regarding mtDNA methylation patterns in samples from patients with different pathologies or that are exposed to a toxic environment in diverse human populations.

Whole-genome characterization of chronological age-associated changes in methylome and circular RNAs in moso bamboo (Phyllostachys edulis) from vegetative to floral growth

In mammals, DNA methylation is associated with aging. However, age-related DNA methylation changes during phase transitions largely remain unstudied in plants. Moso bamboo (Phyllostachys edulis) requires a very long time to transition from the vegetative to the floral phase. To comprehensively investigate the association of DNA methylation with aging, we present here single-base-resolution DNA methylation profiles using both high-throughput bisulfite sequencing and single-molecule nanopore-based DNA sequencing, covering the long period of vegetative growth and transition to flowering in moso bamboo. We discovered that CHH methylation gradually accumulates from vegetative to reproductive growth in a time-dependent fashion. Differentially methylated regions, correlating with chronological aging, occurred preferentially at both transcription start sites and transcription termination sites. Genes with CG methylation changes showed an enrichment of Gene Ontology (GO) categories in 'vegetative to reproductive phase transition of meristem'. Combining methylation data with mRNA sequencing revealed that DNA methylation in promoters, introns and exons may have different roles in regulating gene expression. Finally, circular RNA (circRNA) sequencing revealed that the flanking introns of circRNAs are hypermethylated and enriched in long terminal repeat (LTR) retrotransposons. Together, the observations in this study provide insights into the dynamic DNA methylation and circRNA landscapes, correlating with chronological age, which paves the way to study further the impact of epigenetic factors on flowering in moso bamboo.

Estimating Copy-Number Proportions: The Comeback of Sanger Sequencing

Determination of the relative copy numbers of mixed molecular species in nucleic acid samples is often the objective of biological experiments, including Single-Nucleotide Polymorphism (SNP), indel and gene copy-number characterization, and quantification of CRISPR-Cas9 base editing, cytosine methylation, and RNA editing. Standard dye-terminator chromatograms are a widely accessible, cost-effective information source from which copy-number proportions can be inferred. However, the rate of incorporation of dye terminators is dependent on the dye type, the adjacent sequence string, and the secondary structure of the sequenced strand. These variable rates complicate inferences and have driven scientists to resort to complex and costly quantification methods. Because these complex methods introduce their own biases, researchers are rethinking whether rectifying distortions in sequencing trace files and using direct sequencing for quantification will enable comparable accurate assessment. Indeed, recent developments in software tools (e.g., TIDE, ICE, EditR, BEEP and BEAT) indicate that quantification based on direct Sanger sequencing is gaining in scientific acceptance. This commentary reviews the common obstacles in quantification and the latest insights and developments relevant to estimating copy-number proportions based on direct Sanger sequencing, concluding that bidirectional sequencing and sophisticated base calling are the keys to identifying and avoiding sequence distortions.

Clinical validation of 13-gene DNA methylation analysis in oral brushing samples for detection of oral carcinoma: Italian multicenter study

BACKGROUND

The aim of this Italian multicenter study was to evaluate the diagnostic performance of a minimally invasive method for the detection of oral squamous cell carcinoma (OSCC) based on 13-gene DNA methylation analysis in oral brushing samples.

METHODS

Oral brushing specimens were collected in 11 oral medicine centers across Italy. Twenty brushing specimens were collected by each center, 10 from patients with OSCC, and 10 from healthy volunteers. DNA methylation analysis was performed in blindness, and each sample was determined as positive or negative based on a predefined cutoff value.

RESULTS

DNA amplification failed in 4 of 220 (1.8%) samples. Of the specimens derived from patients with OSCC, 93.6% (103/110) were detected as positive, and 84.9% (90/106) of the samples from healthy volunteers were negative.

CONCLUSION

These data confirmed the diagnostic performance of our novel procedure in a large cohort of brushing specimens collected from 11 different centers and analyzed in blindness.

DNA Methylation Is Responsive to the Environment and Regulates the Expression of Biosynthetic Gene Clusters, Metabolite Production, and Virulence in Fusarium graminearum

Histone modifications play a significant role in the regulation of biosynthetic gene clusters (BGCs) in the phytopathogen Fusarium graminearum, by contrast, epigenetic regulation by DNA methyltransferases (DNMTs) is less documented. In this study, we characterized two DNMTs (FgDIM-2 and FgRID) in F. graminearum, with homologies to "Deficient in methylation" (DIM-2) and "Repeat-induced point (RIP) deficient" (RID) from Neurospora. The loss of DNMTs resulted in not only a decrease in average methylation density in the nutrient-poor, compared to nutrient-rich conditions, but also differences in the genes expressed between the WT and the DNMT mutant strains, implicating the external environment as an important trigger in altering DNA methylation patterns. Consequently, we observed significant changes in the regulation of multiple BGCs and alterations in the pathogenicity of the fungus.

Perinatal DEHP exposure induces sex- and tissue-specific DNA methylation changes in both juvenile and adult mice

Di(2-ethylhexyl) phthalate (DEHP) is a type of phthalate plasticizer found in a variety of consumer products and poses a public health concern due to its metabolic and endocrine disruption activities. Dysregulation of epigenetic modifications, including DNA methylation, has been shown to be an important mechanism for the pathogenic effects of prenatal exposures, including phthalates. In this study, we used an established mouse model to study the effect of perinatal DEHP exposure on the DNA methylation profile in liver (a primary target tissue of DEHP) and blood (a common surrogate tissue) of both juvenile and adult mice. Despite exposure ceasing at 3 weeks of age (PND21), we identified thousands of sex-specific differential DNA methylation events in 5-month old mice, more than identified at PND21, both in blood and liver. Only a small number of these differentially methylated cytosines (DMCs) overlapped between the time points, or between tissues (i.e. liver and blood), indicating blood may not be an appropriate surrogate tissue to estimate the effects of DEHP exposure on liver DNA methylation. We detected sex-specific DMCs common between 3-week and 5-month samples, pointing to specific DNA methylation alterations that are consistent between weanling and adult mice. In summary, this is the first study to assess the genome-wide DNA methylation profiles in liver and blood at two different aged cohorts in response to perinatal DEHP exposure. Our findings cast light on the implications of using surrogate tissue instead of target tissue in human population-based studies and identify epigenetic biomarkers for DEHP exposure.

Effect of 5-azacytidine (5-aza) on UCP2 expression in human liver and colon cancer cells

The function of the uncoupling protein 2 (UCP2) is different for each cancer cell. However, the mechanism of expression is still unclear. DNA methylation affects protein expression and is one factor that transforms normal cells into cancer cells. In this study, the hepatocellular carcinoma Hep3B and HepG2 cells and colorectal cancer HT-29 cells were treated with 5-azacytidine (5-aza), a DNA demethylation agent, to observe the modification of UCP2 expression and the methylation degree in the UCP2 promoter region. Promoter basal activity and degree of UCP2 expression were measured in Hep3B, HepG2, and HT-29 cells. In addition, methylation-specific PCR (MSP) was performed to investigate the degree of methylation in the UCP2 promoter region. The methylation region in the UCP2 promoter was confirmed based on bisulfite sequencing. In Hep3B cells in which UCP2 mRNA was not transcribed, the promoter basal activity was significantly higher than in HT-29 or HepG2 cells in which UCP2 mRNA was transcribed. Treatment with 5-aza increased UCP2 expression in Hep3B and HT-29 cells; however, the expression in HepG2 cells was unchanged. The UCP2 promoter in Hep3B cells has numerous methylated regions compared with HT-29 and HepG2 cells. The results of the present study revealed that inhibition of UCP2 expression in Hep3B cells was due to methylation of the promoter region. Investigating the mechanism that induces UCP2 expression in cancer cells is important to understand the function of UCP2, which could aid in cancer treatment.

LPCAT2 Methylation, a Novel Biomarker for the Severity of Cedar Pollen Allergic Rhinitis in Japan

BACKGROUND

Recently, the role of the epigenome in allergies has been receiving increasing attention. Although several genes that are methylated in relation to serum immunoglobulin E (IgE) concentration have been reported by epigenome-wide association studies, little is known about the DNA methylation sites associated with the symptoms and severity of cedar pollinosis (CP).

OBJECTIVE

Our aim was to analyze the association between DNA methylation and the symptoms and severity of CP in peripheral blood mononuclear cells (PBMCs) and nasal mucosa scraping cells (NMSCs).

METHODS

We recruited 70 participants during the cedar pollen dispersal season. IgE levels were measured by a fluorescence enzyme immunoassay. We analyzed DNA methylation of acyl-CoA thioesterase 7 (ACOT7), mucin 4 (MUC4), schlafen 12 (SLFN12), lysophosphatidylcholine acyltransferase 2 (LPCAT2), and interleukin-4 (IL4) in PBMCs and NMSCs using bisulfite next-generation sequencing; the correlation of DNA methylation with non-specific IgE and cedar pollen-specific IgE levels in peripheral blood samples was also investigated. Symptom severity and DNA methylation were investigated in 15 untreated CP patients.

RESULTS

Non-specific IgE levels showed a significant negative correlation with average IL4 methylation in PBMCs (r = -0.46, P < 0.0001) but not with methylation of ACOT7, MUC4, SLFN12, and LPCAT2. Cedar pollen-specific IgE levels showed a significant negative correlation with average IL4 and MUC4 methylation in PBMCs (r = -0.31, P = 0.01 and r = -0.241, P = 0.046, respectively) but not with methylation of ACOT7, SLFN12, and LPCAT2. The methylation of some genes in NMSCs was not significantly correlated with IgE levels. The mean methylation of LPCAT2 in NMSCs showed a decreasing trend with increasing severity of CP (P = 0.027).

CONCLUSION

LPCAT2 methylation in NMSCs may reflect the severity of CP and could be used as a novel biomarker to identify suitable treatment options for CP.

Inhibition of DNA Methylation in Picochlorum soloecismus Alters Algae Productivity

Eukaryotic organisms regulate the organization, structure, and accessibility of their genomes through chromatin remodeling that can be inherited as epigenetic modifications. These DNA and histone protein modifications are ultimately responsible for an organism's molecular adaptation to the environment, resulting in distinctive phenotypes. Epigenetic manipulation of algae holds yet untapped potential for the optimization of biofuel production and bioproduct formation; however, epigenetic machinery and modes-of-action have not been well characterized in algae. We sought to determine the extent to which the biofuel platform species Picochlorum soloecismus utilizes DNA methylation to regulate its genome. We found candidate genes with domains for DNA methylation in the P. soloecismus genome. Whole-genome bisulfite sequencing revealed DNA methylation in all three cytosine contexts (CpG, CHH, and CHG). While global DNA methylation is low overall (∼1.15%), it occurs in appreciable quantities (12.1%) in CpG dinucleotides in a bimodal distribution in all genomic contexts, though terminators contain the greatest number of CpG sites per kilobase. The P. soloecismus genome becomes hypomethylated during the growth cycle in response to nitrogen starvation. Algae cultures were treated daily across the growth cycle with 20 μM 5-aza-2'-deoxycytidine (5AZA) to inhibit propagation of DNA methylation in daughter cells. 5AZA treatment significantly increased optical density and forward and side scatter of cells across the growth cycle (16 days). This increase in cell size and complexity correlated with a significant increase (∼66%) in lipid accumulation. Site specific CpG DNA methylation was significantly altered with 5AZA treatment over the time course, though nitrogen starvation itself induced significant hypomethylation in CpG contexts. Genes involved in several biological processes, including fatty acid synthesis, had altered methylation ratios in response to 5AZA; we hypothesize that these changes are potentially responsible for the phenotype of early induction of carbon storage as lipids. This is the first report to utilize epigenetic manipulation strategies to alter algal physiology and phenotype. Collectively, these data suggest these strategies can be utilized to fine-tune metabolic responses, alter growth, and enhance environmental adaption of microalgae for desired outcomes.

AmpliconDesign - an interactive web server for the design of high-throughput targeted DNA methylation assays

Targeted analysis of DNA methylation patterns based on bisulfite-treated genomic DNA (BT-DNA) is considered as a gold-standard for epigenetic biomarker development. Existing software tools facilitate primer design, primer quality control or visualization of primer localization. However, high-throughput design of primers for BT-DNA amplification is hampered by limits in throughput and functionality of existing tools, requiring users to repeatedly perform specific tasks manually. Consequently, the design of PCR primers for BT-DNA remains a tedious and time-consuming process. To bridge this gap, we developed AmpliconDesign, a webserver providing a scalable and user-friendly platform for the design and analysis of targeted DNA methylation studies based on BT-DNA, e.g. deep amplicon bisulfite sequencing (ampBS-seq) or EpiTYPER MassArray. Core functionality of the web server includes high-throughput primer design and binding site validation based on in silico bisulfite-converted DNA sequences, prediction of fragmentation patterns for EpiTYPER MassArray, an interactive quality control as well as a streamlined analysis workflow for ampBS-seq.

Downregulation of DUSP9 Promotes Tumor Progression and Contributes to Poor Prognosis in Human Colorectal Cancer

Background

Dual-specificity phosphatase 9 (DUSP9) belongs to the dual-specificity protein phosphatase subfamily. Recently, increasing attention has been paid on the role of DUSP9 in a variety of cancers. However, its functional role in tumor development is still unclear, especially in colorectal cancer (CRC).

Methods

The functional role of DUSP9 in inhibiting the progression of CRC was verified using colony formation assay, wound healing assay, nude mice xenograft model, etc. RNA-seq was performed to assess the gene expression profiling in SW480 cells with DUSP9 stable knockdown and shControl cells. Bisulfite sequencing (BSE) was performed to reveal the methylation status of CpG island in the promoter of DUSP9.

Results

DUSP9 was significantly downregulated in tumor tissues compared with peritumor tissues. Mechanistically, the high methylation status of CpG island in the promoter of DUSP9 may lead to the downregulation of DUSP9 in CRC. Clinically, low DUSP9 expression in CRC was closely associated with depth of invasion, metastasis (TNM) stage, and poor survival, indicating that DUSP9 may be involved in the progression of CRC. Functional study revealed that DUSP9 inhibited proliferation, migration, invasion, and epithelial–mesenchymal transition of CRC cells both in vitro and in vivo. Transcriptome profiling studies revealed that Erk signaling was involved in the tumor progression mediated by DUSP9 silencing, which is confirmed by cell experiments and clinical tissue sample staining analysis.

Conclusion

Our findings demonstrate that DUSP9 plays a critical role in the progression of CRC, and therapeutic intervention to increase the expression or activity of DUSP9 may be a potential target for CRC treatment in the future.

Msuite: A High-Performance and Versatile DNA Methylation Data-Analysis Toolkit

DNA methylation is a pervasive and important epigenetic regulator in mammalian genome. For DNA methylome profiling, emerging bisulfite-free methods have demonstrated desirable superiority over the conventional bisulfite-treatment-based approaches, although current analysis software could not make full use of their advantages. In this work, we present Msuite, an easy-to-use, all-in-one data-analysis toolkit. Msuite implements a unique 4-letter analysis mode specifically optimized for emerging protocols; it also integrates quality controls, methylation call, and data visualizations. Msuite demonstrates substantial performance improvements over current state-of-the-art tools as well as fruitful functionalities, thus holding the potential to serve as an optimal toolkit to facilitate DNA methylome studies. Source codes and testing datasets for Msuite are freely available at https://github.com/hellosunking/Msuite/.

•

Msuite provides a unique 4-letter analysis mode for emerging bisulfite-free protocols

•

Msuite outperforms current tools in terms of higher accuracy and lower resource usage

•

Msuite has integrated quality control and fruitful data-visualization utilities

•

Msuite provides an all-in-one solution for DNA methylation data analysis

DNA methylation is an essential epigenetic modification responsible for many biological regulation pathways. Despite the fact that various high-throughput methods have been developed for base-resolution DNA methylome profiling, DNA methylation data analysis remains a complex and challenging task. Here, we present Msuite, which has integrated quality control, read alignment, methylation call, and fruitful data-visualization functionalities, aiming to offer an all-in-one package for most of the current DNA methylation profiling assays. Msuite also provides dedicated support for emerging bisulfite-free protocols and outperforms the current tools in terms of higher accuracy and lower computational resource requirement. Hence, Msuite could serve as the optimal toolkit for DNA methylation data analysis as well as facilitating the popularization of emerging bisulfite-free protocols.

An Evolutionary Cancer Epigenetic Approach Revealed DNA Hypermethylation of Ultra-Conserved Non-Coding Elements in Squamous Cell Carcinoma of Different Mammalian Species

Background: Ultra-conserved non-coding elements (UCNEs) are genomic sequences that exhibit > 95% sequence identity between humans, mammals, birds, reptiles, and fish. Recent findings reported their functional role in cancer. The aim of this study was to evaluate the DNA methylation modifications of UNCEs in squamous cell carcinoma (SCC) from different mammal species. Methods: Fifty SCCs from 26 humans, 17 cats, 3 dogs, 1 horse, 1 bovine, 1 badger, and 1 porcupine were investigated. Fourteen feline stomatitis and normal samples from 36 healthy human donors, 7 cats, 5 dogs, 5 horses, 2 bovines and 1 badger were collected as normal controls. Bisulfite next generation sequencing evaluated the DNA methylation level from seven UCNEs (uc.160, uc.283, uc.416, uc.339, uc.270, uc.299, and uc.328). Results: 57/59 CpGs were significantly different according to the Kruskal–Wallis test (p < 0.05) comparing normal samples with SCC. A common DNA hypermethylation pattern was observed in SCCs from all the species evaluated in this study, with an increasing trend of hypermethylation starting from normal mucosa, through stomatitis to SCC. Conclusions: Our findings indicate that UCNEs are hypermethylated in human SCC, and this behavior is also conserved among different species of mammals.

Genome-Wide DNA Methylation and Gene Expression Profiles in Cows Subjected to Different Stress Level as Assessed by Cortisol in Milk

Dairy cattle health, wellbeing and productivity are deeply affected by stress. Its influence on metabolism and immune response is well known, but the underlying epigenetic mechanisms require further investigation. In this study, we compared DNA methylation and gene expression signatures between two dairy cattle populations falling in the high- and low-variant tails of the distribution of milk cortisol concentration (MC), a neuroendocrine marker of stress in dairy cows. Reduced Representation Bisulfite Sequencing was used to obtain a methylation map from blood samples of these animals. The high and low groups exhibited similar amounts of methylated CpGs, while we found differences among non-CpG sites. Significant methylation changes were detected in 248 genes. We also identified significant fold differences in the expression of 324 genes. KEGG and Gene Ontology (GO) analysis showed that genes of both groups act together in several pathways, such as nervous system activity, immune regulatory functions and glucocorticoid metabolism. These preliminary results suggest that, in livestock, cortisol secretion could act as a trigger for epigenetic regulation and that peripheral changes in methylation can provide an insight into central nervous system functions.

Corrigendum: Novel Meiotic miRNAs and Indications for a Role of PhasiRNAs in Meiosis

[This corrects the article DOI: 10.3389/fpls.2016.00762.].

DNA methylation analysis in plants: review of computational tools and future perspectives

Genome-wide DNA methylation studies have quickly expanded due to advances in next-generation sequencing techniques along with a wealth of computational tools to analyze the data. Most of our knowledge about DNA methylation profiles, epigenetic heritability and the function of DNA methylation in plants derives from the model species Arabidopsis thaliana. There are increasingly many studies on DNA methylation in plants-uncovering methylation profiles and explaining variations in different plant tissues. Additionally, DNA methylation comparisons of different plant tissue types and dynamics during development processes are only slowly emerging but are crucial for understanding developmental and regulatory decisions. Translating this knowledge from plant model species to commercial crops could allow the establishment of new varieties with increased stress resilience and improved yield. In this review, we provide an overview of the most commonly applied bioinformatics tools for the analysis of DNA methylation data (particularly bisulfite sequencing data). The performances of a selection of the tools are analyzed for computational time and agreement in predicted methylated sites for A. thaliana, which has a smaller genome compared to the hexaploid bread wheat. The performance of the tools was benchmarked on five plant genomes. We give examples of applications of DNA methylation data analysis in crops (with a focus on cereals) and an outlook for future developments for DNA methylation status manipulations and data integration.

A Practical Guide to the Measurement and Analysis of DNA Methylation

DNA methylation represents a fundamental epigenetic mark that is associated with transcriptional repression during development, maintenance of homeostasis, and disease. In addition to methylation-sensitive PCR and targeted deep-amplicon bisulfite sequencing to measure DNA methylation at defined genomic loci, numerous unsupervised techniques exist to quantify DNA methylation on a genome-wide scale, including affinity enrichment strategies and methods involving bisulfite conversion. Both affinity-enriched and bisulfite-converted DNA can serve as input material for array hybridization or sequencing using next-generation technologies. In this practical guide to the measurement and analysis of DNA methylation, the goal is to convey basic concepts in DNA methylation biology and explore genome-scale bisulfite sequencing as the current gold standard for assessment of DNA methylation. Bisulfite conversion chemistry and library preparation are discussed in addition to a bioinformatics approach to quality assessment, trimming, alignment, and methylation calling of individual cytosine residues. Bisulfite-converted DNA presents challenges for standard next-generation sequencing library preparation protocols and data-processing pipelines, but these challenges can be met with elegant solutions that leverage the power of high-performance computing systems. Quantification of DNA methylation, data visualization, statistical approaches to compare DNA methylation between sample groups, and examples of integrating DNA methylation data with other -omics data sets are also discussed. The reader is encouraged to use this article as a foundation to pursue advanced topics in DNA methylation measurement and data analysis, particularly the application of bioinformatics and computational biology principles to generate a deeper understanding of mechanisms linking DNA methylation to cellular function.