Detection of DNA Methylation by Whole-Genome Bisulfite Sequencing

MethMarkerDB: a comprehensive cancer DNA methylation biomarker database

DNA methylation plays a crucial role in tumorigenesis and tumor progression, sparking substantial interest in the clinical applications of cancer DNA methylation biomarkers. Cancer-related whole-genome bisulfite sequencing (WGBS) data offers a promising approach to precisely identify these biomarkers with differentially methylated regions (DMRs). However, currently there is no dedicated resource for cancer DNA methylation biomarkers with WGBS data. Here, we developed a comprehensive cancer DNA methylation biomarker database (MethMarkerDB, https://methmarkerdb.hzau.edu.cn/), which integrated 658 WGBS datasets, incorporating 724 curated DNA methylation biomarker genes from 1425 PubMed published articles. Based on WGBS data, we documented 5.4 million DMRs from 13 common types of cancer as candidate DNA methylation biomarkers. We provided search and annotation functions for these DMRs with different resources, such as enhancers and SNPs, and developed diagnostic and prognostic models for further biomarker evaluation. With the database, we not only identified known DNA methylation biomarkers, but also identified 781 hypermethylated and 5245 hypomethylated pan-cancer DMRs, corresponding to 693 and 2172 genes, respectively. These novel potential pan-cancer DNA methylation biomarkers hold significant clinical translational value. We hope that MethMarkerDB will help identify novel cancer DNA methylation biomarkers and propel the clinical application of these biomarkers.

EpiTyping: analysis of epigenetic aberrations in parental imprinting and X-chromosome inactivation using RNA-seq

Human pluripotent stem cells (hPSCs) hold a central role in studying human development, in disease modeling and in regenerative medicine. These cells not only acquire genetic modifications when kept in culture, but they may also harbor epigenetic aberrations, mainly involving parental imprinting and X-chromosome inactivation. Here we present a detailed bioinformatic protocol for detecting such aberrations using RNA sequencing data. We provide a pipeline designed to process and analyze RNA sequencing data for the identification of abnormal biallelic expression of imprinted genes, and thus detect loss of imprinting. Furthermore, we show how to differentiate among X-chromosome inactivation, full activation and aberrant erosion of X chromosome in female hPSCs. In addition to providing bioinformatic tools, we discuss the impact of such epigenetic variations in hPSCs on their utility for various purposes. This pipeline can be used by any user with basic understanding of the Linux command line. It is available on GitHub as a software container ( https://github.com/Gal-Keshet/EpiTyping ) and produces reliable results in 1-4 d.

Integrative epigenome profiling of 47XXY provides insights into whole genomic DNA hypermethylation and active chromatin accessibility

Klinefelter syndrome (KS, 47XXY) is a disorder characterized by sex chromosomal aneuploidy, which may lead to changes in epigenetic regulations of gene expression. To define epigenetic architectures in 47XXY, we annotated DNA methylation in euploid males (46XY) and females (46XX), and 47XXY individuals using whole genome bisulfite sequencing (WGBS) and integrated chromatin accessbilty, and detected abnormal hypermethylation in 47XXY. Furthermore, we detected altered chromatin accessibility in 47XXY, in particular in chromosome X, using Assay for Transposase-Accessible Chromatin sequencing (ATAC-seq) in cultured amniotic cells. Our results construct the whole genome-wide DNA methylation map in 47XXY, and provide new insights into the early epigenomic dysregulation resulting from an extra chromosome X in 47XXY.

A nanopore counter for highly sensitive evaluation of DNA methylation and its application in in vitro diagnostics

DNA methylation has been considered an essential epigenetic biomarker for diagnosing various diseases, such as cancer. A simple and sensitive way for DNA methylation level detection is necessary. Inspired by the label-free and ultra-high sensitivity of solid-state nanopores to double-stranded DNA (dsDNA), we proposed a nanopore counter for evaluating DNA methylation by integrating a dual-restriction endonuclease digestion strategy coupled with polymerase chain reaction (PCR) amplification. Simultaneous application of BstUI/HhaI endonucleases can ensure the full digestion of the unmethylated target DNA but shows no effect on the methylated ones. Therefore, only the methylated DNA remains intact and can trigger the subsequent PCR reaction, producing a large quantity of fixed-length PCR amplicons, which can be directly detected through glassy nanopores. By simply counting the event rate of the translocation signals, the concentration of methylated DNA can be determined to range from 1 aM to 0.1 nM, with the detection limit as low as 0.61 aM. Moreover, a 0.01% DNA methylation level was successfully distinguished. The strategy of using the nanopore counter for highly sensitive DNA methylation evaluation would be a low-cost but reliable alternative in the analysis of DNA methylation.

Multi-Omics Approaches to Improve Clubroot Resistance in Brassica with a Special Focus on Brassica oleracea L

Brassica oleracea is an agronomically important species of the Brassicaceae family, including several nutrient-rich vegetables grown and consumed across the continents. But its sustainability is heavily constrained by a range of destructive pathogens, among which, clubroot disease, caused by a biotrophic protist Plasmodiophora brassicae, has caused significant yield and economic losses worldwide, thereby threatening global food security. To counter the pathogen attack, it demands a better understanding of the complex phenomenon of Brassica-P. brassicae pathosystem at the physiological, biochemical, molecular, and cellular levels. In recent years, multiple omics technologies with high-throughput techniques have emerged as successful in elucidating the responses to biotic and abiotic stresses. In Brassica spp., omics technologies such as genomics, transcriptomics, ncRNAomics, proteomics, and metabolomics are well documented, allowing us to gain insights into the dynamic changes that transpired during host-pathogen interactions at a deeper level. So, it is critical that we must review the recent advances in omics approaches and discuss how the current knowledge in multi-omics technologies has been able to breed high-quality clubroot-resistant B. oleracea. This review highlights the recent advances made in utilizing various omics approaches to understand the host resistance mechanisms adopted by Brassica crops in response to the P. brassicae attack. Finally, we have discussed the bottlenecks and the way forward to overcome the persisting knowledge gaps in delivering solutions to breed clubroot-resistant Brassica crops in a holistic, targeted, and precise way.

Heterosis of growth trait regulated by DNA methylation and miRNA in allotriploid fish

Background

Heterosis of growth traits in allotriploid fish has benefited the production of aquaculture for many years, yet its genetic and molecular basis has remained obscure. Now, an allotriploid complex, including two triploids and their diploid inbred parents, has provided an excellent model for investigating the potential regulatory mechanisms of heterosis.

Results

Here, we performed a series of analyses on DNA methylation modification and miRNA expression in combination with gene expression in the allotriploid complex. We first established a model of cis- and trans-regulation related to DNA methylation and miRNA in allotriploids. Then, comparative analyses showed that DNA methylation contributed to the emergence of a dosage compensation effect, which reduced gene expression levels in the triploid to the diploid state. We detected 31 genes regulated by DNA methylation in the subgenomes of the allotriploids. Finally, the patterns of coevolution between small RNAs and their homoeologous targets were classified and used to predict the regulation of miRNA expression in the allotriploids.

Conclusions

Our results uncovered the regulatory network between DNA methylation and miRNAs in allotriploids, which not only helps us understand the regulatory mechanisms of heterosis of growth traits but also benefits the study and application of epigenetics in aquaculture.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13072-022-00455-6.

UHRF1 establishes crosstalk between somatic and germ cells in male reproduction

Sertoli cells (SCs) support and nourish germ cells (GCs) through their crosstalk during spermatogenesis. However, the underlying epigenetic mechanism that ensures SCs’ functions in this process remains unclear. Here, we report that UHRF1, a critical epigenetic regulator, is mainly expressed in human and mouse pre-mature SCs, and is essential for establishing Sertoli-Germ cell crosstalk. SC-specific UHRF1 knockout mice exhibit complete sterility with Sertoli cell (SC) proliferation and differentiation aberrance, blood-testis barrier (BTB) disruption, and immature germ cell (GC) sloughing. RNA sequencing and Whole Genome Bisulfite Sequencing (WGBS) revealed that many extracellular matrix (ECM)-related genes (e.g., Timp1, Trf, and Spp1) appeared upregulated with the DNA hypomethylation status in UHRF1-deficient SCs. Strikingly, overexpression of Timp1, Trf, and Spp1 in SCs in vitro and in vivo could phenocopy the SC-specific UHRF1-deficient mice. Our data demonstrated that UHRF1 regulates the transcriptional program of ECM-related genes in SCs and establishes SC-GC crosstalk.

Particulate matter-induced hypomethylation of Alu and LINE1 in normal human bronchial epithelial cells and epidermal keratinocytes

BACKGROUND

Airborne particulate matter (PM), a complex mixture of organic and inorganic compounds, is a major public health concern due to its adverse health effects. Understanding the biological action of PM is of particular importance in the improvement of public health. Differential methylation of repetitive elements (RE) by PM might have severe consequences for the structural integrity of the genome and on transcriptional activity, thereby affecting human health. This study aimed to evaluate the effect of inhaled and non-inhaled PM (PM2.5, PM10, and PM10-PAH) exposure on DNA methylation. We quantitatively measured the methylation content of Alu and LINE1 in PM-treated normal human bronchial epithelial cells (NHBE) and normal human epidermal keratinocytes (NHEK) by using whole-genome bisulfite sequencing and pyrosequencing.

RESULTS

All PMs exposure significantly lowered Alu and LINE1 methylation in both cells than in mock-treated controls. Hypomethylation was more prominent in PM10-PAH exposed-NHBE and PM10 exposed-NHEK. Alu and LINE1 methylation change exhibited different sensitivity according to the subfamily evolutionary ages, with stronger effects on the oldest L1-M and Alu J in NHBE, and oldest L1-M and youngest Alu S in NHEK.

CONCLUSIONS

These results demonstrate that the differential susceptibility of PM-induced hypomethylation of Alu and LINE1 depends upon RE evolutionary age and PM type.

Concerto on Chromatin: Interplays of Different Epigenetic Mechanisms in Plant Development and Environmental Adaptation

Epigenetic mechanisms such as DNA methylation, histone post-translational modifications, chromatin remodeling, and noncoding RNAs, play important roles in regulating plant gene expression, which is involved in various biological processes including plant development and stress responses. Increasing evidence reveals that these different epigenetic mechanisms are highly interconnected, thereby contributing to the complexity of transcriptional reprogramming in plant development processes and responses to environmental stresses. Here, we provide an overview of recent advances in understanding the epigenetic regulation of plant gene expression and highlight the crosstalk among different epigenetic mechanisms in making plant developmental and stress-responsive decisions. Structural, physical, transcriptional and metabolic bases for these epigenetic interplays are discussed.

Rapid Detection of DNA Methylation with a Novel Real-Time Fluorescence Recombinase-Aided Amplification Assay

Researchers have conducted in-depth research on DNA methylation mechanism, which is related to various diseases such as deficiency of imprinted gene and occurrence of tumors. This study provides a novel rapid quantitative detection assay and real-time fluorescence recombinase-aided amplification assay (RAA) for DNA methylation. Firstly, specific sequence of methylation genes was chosen and primers and fluorogenic probe for RAA experiment were designed and synthesized. Lastly, these amplification products were proven by sequencing and analysis. Results showed that the amplification efficiency and template concentration of RAA had linear dependent (R² > 95%) when the concentration range was 4.64×10⁸ copies/μL˜4.64×10⁴ copies/μL. The test assay can also detect positive samples when the template concentration is below 4.64×10⁴ copies/μL. Remarkably, the entire experiment process only takes 15-20 minutes, so it is beneficial for rapid bedside simple screening of some special DNA methylation sites, such as detection of resistance genes. In a word, this method has very great potential for diseases with DNA methylation in clinical settings, especially if methylation analysis needs to be done quickly and easily.

Exploiting Epigenetic Variations for Crop Disease Resistance Improvement

Pathogen infections seriously threaten plant health and global crop production. Epigenetic processes such as DNA methylation, histone post-translational modifications, chromatin assembly and remodeling play important roles in transcriptional regulation of plant defense responses and could provide a new direction to drive breeding strategies for crop disease resistance improvement. Although past decades have seen unprecedented proceedings in understanding the epigenetic mechanism of plant defense response, most of these advances were derived from studies in model plants like Arabidopsis. In this review, we highlighted the recent epigenetic studies on crop-pathogen interactions and discussed the potentials, challenges, and strategies in exploiting epigenetic variations for crop disease resistance improvement.

Chicken cecal DNA methylome alteration in the response to Salmonella enterica serovar Enteritidis inoculation

BACKGROUND

Salmonella enterica serovar Enteritidis (SE) is one of the pathogenic bacteria, which affects poultry production and poses a severe threat to public health. Chicken meat and eggs are the main sources of human salmonellosis. DNA methylation is involved in regulatory processes including gene expression, chromatin structure and genomic imprinting. To understand the methylation regulation in the response to SE inoculation in chicken, the genome-wide DNA methylation profile following SE inoculation was analyzed through whole-genome bisulfite sequencing in the current study.

RESULTS

There were 185,362,463 clean reads and 126,098,724 unique reads in the control group, and 180,530,750 clean reads and 126,782,896 unique reads in the inoculated group. The methylation density in the gene body was higher than that in the upstream and downstream regions of the gene. There were 8946 differentially methylated genes (3639 hypo-methylated genes, 5307 hyper-methylated genes) obtained between inoculated and control groups. Methylated genes were mainly enriched in immune-related Gene Ontology (GO) terms and metabolic process terms. Cytokine-cytokine receptor interaction, TGF-beta signaling pathway, FoxO signaling pathway, Wnt signaling pathway and several metabolism-related pathways were significantly enriched. The density of differentially methylated cytosines in miRNAs was the highest. HOX genes were widely methylated.

CONCLUSIONS

The genome-wide DNA methylation profile in the response to SE inoculation in chicken was analyzed. SE inoculation promoted the DNA methylation in the chicken cecum and caused methylation alteration in immune- and metabolic- related genes. Wnt signal pathway, miRNAs and HOX gene family may play crucial roles in the methylation regulation of SE inoculation in chicken. The findings herein will deepen the understanding of epigenetic regulation in the response to SE inoculation in chicken.

Epigenetics: Recent Advances and Its Role in the Treatment of Alzheimer's Disease

Objective: This review summarizes recent findings on the epigenetics of Alzheimer's disease (AD) and provides therapeutic strategies for AD.

Methods: We searched the following keywords: “genetics,” “epigenetics,” “Alzheimer's disease,” “DNA methylation,” “DNA hydroxymethylation,” “histone modifications,” “non-coding RNAs,” and “therapeutic strategies” in PubMed.

Results: In this review, we summarizes recent studies of epigenetics in AD, including DNA methylation/hydroxymethylation, histone modifications, and non-coding RNAs. There are no consistent results of global DNA methylation/hydroxymethylation in AD. Epigenetic genome-wide association studies show that many differentially methylated sites exist in AD. Several studies investigate the role of histone modifications in AD; for example, histone acetylation decreases, whereas H3 phosphorylation increases significantly in AD. In addition, non-coding RNAs, such as microRNA-16 and BACE1-antisense transcript (BACE1-AS), are associated with the pathology of AD. These epigenetic changes provide us with novel insights into the pathogenesis of AD and may be potential therapeutic strategies for AD.

Conclusion: Epigenetics is associated with the pathogenesis of AD, including DNA methylation/hydroxymethylation, histone modifications, and non-coding RNAs, which provide potential therapeutic strategies for AD.

Efficient and accurate determination of genome-wide DNA methylation patterns in Arabidopsis thaliana with enzymatic methyl sequencing

Background

5′ methylation of cytosines in DNA molecules is an important epigenetic mark in eukaryotes. Bisulfite sequencing is the gold standard of DNA methylation detection, and whole-genome bisulfite sequencing (WGBS) has been widely used to detect methylation at single-nucleotide resolution on a genome-wide scale. However, sodium bisulfite is known to severely degrade DNA, which, in combination with biases introduced during PCR amplification, leads to unbalanced base representation in the final sequencing libraries. Enzymatic conversion of unmethylated cytosines to uracils can achieve the same end product for sequencing as does bisulfite treatment and does not affect the integrity of the DNA; enzymatic methylation sequencing may, thus, provide advantages over bisulfite sequencing.

Results

Using an enzymatic methyl-seq (EM-seq) technique to selectively deaminate unmethylated cytosines to uracils, we generated and sequenced libraries based on different amounts of Arabidopsis input DNA and different numbers of PCR cycles, and compared these data to results from traditional whole-genome bisulfite sequencing. We found that EM-seq libraries were more consistent between replicates and had higher mapping and lower duplication rates, lower background noise, higher average coverage, and higher coverage of total cytosines. Differential methylation region (DMR) analysis showed that WGBS tended to over-estimate methylation levels especially in CHG and CHH contexts, whereas EM-seq detected higher CG methylation levels in certain highly methylated areas. These phenomena can be mostly explained by a correlation of WGBS methylation estimation with GC content and methylated cytosine density. We used EM-seq to compare methylation between leaves and flowers, and found that CHG methylation level is greatly elevated in flowers, especially in pericentromeric regions.

Conclusion

We suggest that EM-seq is a more accurate and reliable approach than WGBS to detect methylation. Compared to WGBS, the results of EM-seq are less affected by differences in library preparation conditions or by the skewed base composition in the converted DNA. It may therefore be more desirable to use EM-seq in methylation studies.

Nick-seq for single-nucleotide resolution genomic maps of DNA modifications and damage

DNA damage and epigenetic marks are well established to have profound influences on genome stability and cell phenotype, yet there are few technologies to obtain high-resolution genomic maps of the many types of chemical modifications of DNA. Here we present Nick-seq for quantitative, sensitive, and accurate mapping of DNA modifications at single-nucleotide resolution across genomes. Pre-existing breaks are first blocked and DNA modifications are then converted enzymatically or chemically to strand-breaks for both 3'-extension by nick-translation to produce nuclease-resistant oligonucleotides and 3'-terminal transferase tailing. Following library preparation and next generation sequencing, the complementary datasets are mined with a custom workflow to increase sensitivity, specificity and accuracy of the map. The utility of Nick-seq is demonstrated with genomic maps of site-specific endonuclease strand-breaks in purified DNA from Eschericia coli, phosphorothioate epigenetics in Salmonella enterica Cerro 87, and oxidation-induced abasic sites in DNA from E. coli treated with a sublethal dose of hydrogen peroxide. Nick-seq applicability is demonstrated with strategies for >25 types of DNA modification and damage.

3D genome architecture coordinates trans and cis regulation of differentially expressed ear and tassel genes in maize

BACKGROUND

Maize ears and tassels are two separate types of inflorescence which are initiated by similar developmental processes but gradually develop distinct architectures. However, coordinated trans and cis regulation of differentially expressed genes determining ear and tassel architecture within the 3D genome context is largely unknown.

RESULTS

We identify 56,055 and 52,633 open chromatin regions (OCRs) in developing maize ear and tassel primordia using ATAC-seq and characterize combinatorial epigenome features around these OCRs using ChIP-seq, Bisulfite-seq, and RNA-seq datasets. Our integrative analysis of coordinated epigenetic modification and transcription factor binding to OCRs highlights the cis and trans regulation of differentially expressed genes in ear and tassel controlling inflorescence architecture. We further systematically map chromatin interactions at high-resolution in corresponding tissues using in situ digestion-ligation-only Hi-C (DLO Hi-C). The extensive chromatin loops connecting OCRs and genes provide a 3D view on cis- and trans-regulatory modules responsible for ear- and tassel-specific gene expression. We find that intergenic SNPs tend to locate in distal OCRs, and our chromatin interaction maps provide a potential mechanism for trait-associated intergenic SNPs that may contribute to phenotypic variation by influencing target gene expression through chromatin loops.

CONCLUSIONS

Our comprehensive epigenome annotations and 3D genome maps serve as valuable resource and provide a deep understanding of the complex regulatory mechanisms of genes underlying developmental and morphological diversities between maize ear and tassel.

Unravelling the complex genetics of common kidney diseases: from variants to mechanisms

Genome-wide association studies (GWAS) have identified hundreds of loci associated with kidney-related traits such as glomerular filtration rate, albuminuria, hypertension, electrolyte and metabolite levels. However, these impressive, large-scale mapping approaches have not always translated into an improved understanding of disease or development of novel therapeutics. GWAS have several important limitations. Nearly all disease-associated risk loci are located in the non-coding region of the genome and therefore, their target genes, affected cell types and regulatory mechanisms remain unknown. Genome-scale approaches can be used to identify associations between DNA sequence variants and changes in gene expression (quantified through bulk and single-cell methods), gene regulation and other molecular quantitative trait studies, such as chromatin accessibility, DNA methylation, protein expression and metabolite levels. Data obtained through these approaches, used in combination with robust computational methods, can deliver robust mechanistic inferences for translational exploitation. Understanding the genetic basis of common kidney diseases means having a comprehensive picture of the genes that have a causal role in disease development and progression, of the cells, tissues and organs in which these genes act to affect the disease, of the cellular pathways and mechanisms that drive disease, and of potential targets for disease prevention, detection and therapy.

Genome-Wide DNA Methylation Analysis of Soybean Curled-Cotyledons Mutant and Functional Evaluation of a Homeodomain-Leucine Zipper (HD-Zip) I Gene GmHDZ20

DNA methylation is a major, conserved epigenetic modification that influences many biological processes. Cotyledons are specialized tissues that provide nutrition for seedlings at the early developmental stage. To investigate the patterns of genomic DNA methylation of germinated cotyledons in soybean (Glycine max) and its effect on cotyledon development, we performed a genome-wide comparative analysis of DNA methylation between the soybean curled-cotyledons (cco) mutant, which has abnormal cotyledons, and its corresponding wild type (WT) by whole-genome bisulfite sequencing. The cco mutant was methylated at more sites but at a slightly lower level overall than the WT on the whole-genome level. A total of 46 CG-, 92 CHG-, and 9723 CHH- (H = A, C, or T) differentially methylated genes (DMGs) were identified in cotyledons. Notably, hypomethylated CHH-DMGs were enriched in the gene ontology term "sequence-specific DNA binding transcription factor activity." We selected a DMG encoding a homeodomain-leucine zipper (HD-Zip) I subgroup transcription factor (GmHDZ20) for further functional characterization. GmHDZ20 localized to the nucleus and was highly expressed in leaf and cotyledon tissues. Constitutive expression of GmHDZ20 in Arabidopsis thaliana led to serrated rosette leaves, shorter siliques, and reduced seed number per silique. A yeast two-hybrid assay revealed that GmHDZ20 physically interacted with three proteins associated with multiple aspects of plant growth. Collectively, our results provide a comprehensive study of soybean DNA methylation in normal and aberrant cotyledons, which will be useful for the identification of specific DMGs that participate in cotyledon development, and also provide a foundation for future in-depth functional study of GmHDZ20 in soybean.

D-GPM: A Deep Learning Method for Gene Promoter Methylation Inference

Whole-genome bisulfite sequencing generates a comprehensive profiling of the gene methylation levels, but is limited by a high cost. Recent studies have partitioned the genes into landmark genes and target genes and suggested that the landmark gene expression levels capture adequate information to reconstruct the target gene expression levels. This inspired us to propose that the methylation level of the promoters in landmark genes might be adequate to reconstruct the promoter methylation level of target genes, which would eventually reduce the cost of promoter methylation profiling. Here, we propose a deep learning model called Deep-Gene Promoter Methylation (D-GPM) to predict the whole-genome promoter methylation level based on the promoter methylation profile of the landmark genes from The Cancer Genome Atlas (TCGA). D-GPM-15%-7000 × 5, the optimal architecture of D-GPM, acquires the least overall mean absolute error (MAE) and the highest overall Pearson correlation coefficient (PCC), with values of 0.0329 and 0.8186, respectively, when testing data. Additionally, the D-GPM outperforms the regression tree (RT), linear regression (LR), and the support vector machine (SVM) in 95.66%, 92.65%, and 85.49% of the target genes by virtue of its relatively lower MAE and in 98.25%, 91.00%, and 81.56% of the target genes based on its relatively higher PCC, respectively. More importantly, the D-GPM predominates in predicting 79.86% and 78.34% of the target genes according to the model distribution of the least MAE and the highest PCC, respectively.

DNA methylation detection methods used in colorectal cancer

Colorectal cancer (CRC) remains a major contributor to the number of cancer-related deaths that occur annually worldwide. With the development of molecular biology methods, an increasing number of molecular biomarkers have been identified and investigated. CRC is believed to result from an accumulation of epigenetic changes, and detecting aberrant DNA methylation patterns is useful for both the early diagnosis and prognosis of CRC. Numerous studies are focusing on the development of DNA methylation detection methods or DNA methylation panels. Thus, this review will discuss the commonly used techniques and technologies to evaluate DNA methylation, their merits and deficiencies as well as the prospects for new methods.

Genome-Wide DNA Methylation Profiling in the Lotus (Nelumbo nucifera) Flower Showing its Contribution to the Stamen Petaloid

DNA methylation is a vital epigenetic modification. Methylation has a significant effect on the gene expression influencing the regulation of different physiological processes. Current studies on DNA methylation have been conducted on model plants. Lotus (Nelumbo nucifera) is a basic eudicot exhibiting variations during development, especially in flower formation. DNA methylation profiling was conducted on different flower tissues of lotuses through whole genome bisulfite sequencing (WGBS) to investigate the effects of DNA methylation on its stamen petaloid. A map of methylated cytosines at the single base pair resolution for the lotus was constructed. When the stamen was compared with the stamen petaloid, the DNA methylation exhibited a global decrease. Genome-wide relationship analysis between DNA methylation and gene expression identified 31 different methylation region (DMR)-associated genes, which might play crucial roles in floral organ formation, especially in the stamen petaloid. One out of 31 DMR-associated genes, NNU_05638 was homolog with Plant U-box 33 (PUB33). The DNA methylation status of NNU_05638 promoter was distinct in three floral organs, which was confirmed by traditional bisulfite sequencing. These results provide further insights about the regulation of stamen petaloids at the epigenetic level in lotus.

DNA Methylation Analysis of the Citrullus lanatus Response to Cucumber Green Mottle Mosaic Virus Infection by Whole-Genome Bisulfite Sequencing

DNA methylation is an important epigenetic mark associated with plant immunity, butlittle is known about its roles in viral infection of watermelon. We carried out whole-genomebisulfite sequencing of watermelon leaves at 0 h (ck), 48 h, and 25 days post-inoculation withCucumber green mottle mosaic virus (CGMMV). The number of differentially methylated regions(DMRs) increased during CGMMV infection and 2788 DMR-associated genes (DMGs) werescreened out among three libraries. Most DMRs and DMGs were obtained under the CHH context.These DMGs were significantly enriched in the Kyoto Encyclopedia of Genes and Genomes (KEGG)pathways of secondary biosynthesis and metabolism, plant-pathogen interactions, Toll-likereceptor signaling, and ABC transporters. Additionally, DMGs encoding PR1a, CaMs, calciumbindingprotein, RIN4, BAK1, WRKYs, RBOHs, STKs, and RLPs/RLKs were involved in thewatermelon-CGMMV interaction and signaling. The association between DNA methylation andgene expression was analyzed by RNA-seq and no clear relationship was detected. Moreover,downregulation of genes in the RdDM pathway suggested the reduced RdDM-directed CHHmethylation plays an important role in antiviral defense in watermelon. Our findings providegenome-wide DNA methylation profiles of watermelon and will aid in revealing the molecularmechanism in response to CGMMV infection at the methylation level.