Shotgun bisulphite sequencing of the Arabidopsis genome reveals DNA methylation patterning

Expression of genes related to gonadal development and construction of gonadal DNA methylation maps of Trachinotus blochii under hypoxia

Chronic hypoxia can affect the growth and metabolism of fish and potentially impact gonadal development through epigenetic regulation. Trachinotus blochii (Golden pompano) is widely cultured near the coast and is sensitive to low oxygen conditions. We found that hypoxia and reoxygenation processes acted on multiple targets on the HPG axis, leading to endocrine disorders. Changes in the expression of key genes in the brain (gnrh), pituitary (fsh and lh), ovaries (cyp19a1a, foxl2, and er), and testes (dmrt1, ar, sox9, and gsdf) were associated with significant decreases in estrogen and testosterone levels. Hypoxia and reoxygenation lead to changes in DNA methylation levels in the gonads. Hypoxia upregulated the expression of dnmt1, dnmt3a, dnmt3b, tet1, and tet2 in females and dnmt3a and dnmt3b in males, while reoxygenation down-regulated the expression of dnmt1, dnmt3a, dnmt3b, tet1, and tet2 in males. Whole genome methylation sequencing showed that the number of differentially methylated regions was highest on chromosome 10 (5192) and lowest on chromosome 24 (275). Differentially methylated genes in females and males, as well as between males and females, were enriched in the oxytocin signaling pathway, fatty acid metabolism pathway, and HIF-1a pathway. In summary, hypoxia and reoxygenation can induce endocrine disorders, affect the expression of HPG axis genes, change the methylation pattern and modification pattern of gonad DNA, and then have potential effects on gonad development.

Global level of methylation in the sea lamprey (jawless vertebrate) genome is intermediate between invertebrate and jawed vertebrate genomes

In eukaryotes, cytosine methylation is a primary heritable epigenetic modification of the genome that regulates many cellular processes. In invertebrate, methylated cytosine generally located on specific genomic elements (e.g., gene bodies and silenced repetitive elements) to show a "mosaic" pattern. While in jawed vertebrate (teleost and tetrapod), highly methylated cytosine located genome-wide but only absence at regulatory regions (e.g., promoter and enhancer). Many studies imply that the evolution of DNA methylation reprogramming may have helped the transition from invertebrates to jawed vertebrates, but the detail remains largely elusive. In this study, we used the whole-genome bisulfite-sequencing technology to investigate the genome-wide methylation in three tissues (heart, muscle, and sperm) from the sea lamprey, an extant agnathan (jawless) vertebrate. Strikingly, we found that the methylation level of the sea lamprey is very similar to that in sea urchin (a deuterostome) and sea squirt (a chordate) invertebrates. In sum, the global pattern in sea lamprey is intermediate methylation level (around 30%), that is higher than methylation level in the genomes of pre-bilaterians and protostomes (1%-10%), but lower than methylation level appeared in jawed vertebrates (around 70%, teleost and tetrapod). We anticipate that, in addition to genetic dynamics such as genome duplications, epigenetic dynamics such as global methylation reprograming was also orchestrated toward the emergence and evolution of vertebrates.

De novo assembly of transcriptomes and differential gene expression analysis using short-read data from emerging model organisms - a brief guide

Many questions in biology benefit greatly from the use of a variety of model systems. High-throughput sequencing methods have been a triumph in the democratization of diverse model systems. They allow for the economical sequencing of an entire genome or transcriptome of interest, and with technical variations can even provide insight into genome organization and the expression and regulation of genes. The analysis and biological interpretation of such large datasets can present significant challenges that depend on the 'scientific status' of the model system. While high-quality genome and transcriptome references are readily available for well-established model systems, the establishment of such references for an emerging model system often requires extensive resources such as finances, expertise and computation capabilities. The de novo assembly of a transcriptome represents an excellent entry point for genetic and molecular studies in emerging model systems as it can efficiently assess gene content while also serving as a reference for differential gene expression studies. However, the process of de novo transcriptome assembly is non-trivial, and as a rule must be empirically optimized for every dataset. For the researcher working with an emerging model system, and with little to no experience with assembling and quantifying short-read data from the Illumina platform, these processes can be daunting. In this guide we outline the major challenges faced when establishing a reference transcriptome de novo and we provide advice on how to approach such an endeavor. We describe the major experimental and bioinformatic steps, provide some broad recommendations and cautions for the newcomer to de novo transcriptome assembly and differential gene expression analyses. Moreover, we provide an initial selection of tools that can assist in the journey from raw short-read data to assembled transcriptome and lists of differentially expressed genes.

An epigenetic toolbox for conservation biologists

Ongoing climatic shifts and increasing anthropogenic pressures demand an efficient delineation of conservation units and accurate predictions of populations' resilience and adaptive potential. Molecular tools involving DNA sequencing are nowadays routinely used for these purposes. Yet, most of the existing tools focusing on sequence-level information have shortcomings in detecting signals of short-term ecological relevance. Epigenetic modifications carry valuable information to better link individuals, populations, and species to their environment. Here, we discuss a series of epigenetic monitoring tools that can be directly applied to various conservation contexts, complementing already existing molecular monitoring frameworks. Focusing on DNA sequence-based methods (e.g. DNA methylation, for which the applications are readily available), we demonstrate how (a) the identification of epi-biomarkers associated with age or infection can facilitate the determination of an individual's health status in wild populations; (b) whole epigenome analyses can identify signatures of selection linked to environmental conditions and facilitate estimating the adaptive potential of populations; and (c) epi-eDNA (epigenetic environmental DNA), an epigenetic-based conservation tool, presents a non-invasive sampling method to monitor biological information beyond the mere presence of individuals. Overall, our framework refines conservation strategies, ensuring a comprehensive understanding of species' adaptive potential and persistence on ecologically relevant timescales.

Mechanism of heterochromatin remodeling revealed by the DDM1 bound nucleosome structures

The SWI/SNF2 chromatin remodeling factor decreased DNA methylation 1 (DDM1) is essential for the silencing of transposable elements (TEs) in both euchromatic and heterochromatic regions. Here, we determined the cryo-EM structures of DDM1-nucleosomeH2A and DDM1-nucleosomeH2A.W complexes at near-atomic resolution in the presence of the ATP analog ADP-BeFx. The structures show that nucleosomal DNA is unwrapped more on the surface of the histone octamer containing histone H2A than that containing histone H2A.W. DDM1 embraces one DNA gyre of the nucleosome and interacts with the N-terminal tails of histone H4. Although we did not observe DDM1-H2A.W interactions in our structures, the results of the pull-down experiments suggest a direct interaction between DDM1 and the core region of histone H2A.W. Our work provides mechanistic insights into the heterochromatin remodeling process driven by DDM1 in plants.

Improve-RRBS: a novel tool to correct the 3' trimming of reduced representation sequencing reads

Motivation

Reduced Representation Bisulfite Sequencing (RRBS) is a popular approach to determine DNA methylation of the CpG-rich regions of the genome. However, we observed that false positive differentially methylated sites (DMS) are also identified using the standard computational analysis.

Results

During RRBS library preparation the MspI digested DNA undergo end-repair by a cytosine at the 3' end of the fragments. After sequencing, Trim Galore cuts these end-repaired nucleotides. However, Trim Galore fails to detect end-repair when it overlaps with the 3' end of the sequencing reads. We found that these non-trimmed cytosines bias methylation calling, thus, can identify DMS erroneously. To circumvent this problem, we developed improve-RRBS, which efficiently identifies and hides these cytosines from methylation calling with a false positive rate of maximum 0.5%. To test improve-RRBS, we investigated four datasets from four laboratories and two different species. We found non-trimmed 3' cytosines in all datasets analyzed and as much as >50% of false positive DMS under certain conditions. By applying improve-RRBS, these DMS completely disappeared from all comparisons.

Availability and implementation

Improve-RRBS is a freely available python package https://pypi.org/project/iRRBS/ or https://github.com/fothia/improve-RRBS to be implemented in RRBS pipelines.

Genome wide inherited modifications of the tomato epigenome by trans-activated bacterial CG methyltransferase

BACKGROUND

Epigenetic variation is mediated by epigenetic marks such as DNA methylation occurring in all cytosine contexts in plants. CG methylation plays a critical role in silencing transposable elements and regulating gene expression. The establishment of CG methylation occurs via the RNA-directed DNA methylation pathway and CG methylation maintenance relies on METHYLTRANSFERASE1, the homologue of the mammalian DNMT1.

PURPOSE

Here, we examined the capacity to stably alter the tomato genome methylome by a bacterial CG-specific M.SssI methyltransferase expressed through the LhG4/pOP transactivation system.

RESULTS

Methylome analysis of M.SssI expressing plants revealed that their euchromatic genome regions are specifically hypermethylated in the CG context, and so are most of their genes. However, changes in gene expression were observed only with a set of genes exhibiting a greater susceptibility to CG hypermethylation near their transcription start site. Unlike gene rich genomic regions, our analysis revealed that heterochromatic regions are slightly hypomethylated at CGs only. Notably, some M.SssI-induced hypermethylation persisted even without the methylase or transgenes, indicating inheritable epigenetic modification.

CONCLUSION

Collectively our findings suggest that heterologous expression of M.SssI can create new inherited epigenetic variations and changes in the methylation profiles on a genome wide scale. This open avenues for the conception of epigenetic recombinant inbred line populations with the potential to unveil agriculturally valuable tomato epialleles.

Population genomics and epigenomics of Spirodela polyrhiza provide insights into the evolution of facultative asexuality

Many plants are facultatively asexual, balancing short-term benefits with long-term costs of asexuality. During range expansion, natural selection likely influences the genetic controls of asexuality in these organisms. However, evidence of natural selection driving asexuality is limited, and the evolutionary consequences of asexuality on the genomic and epigenomic diversity remain controversial. We analyzed population genomes and epigenomes of Spirodela polyrhiza, (L.) Schleid., a facultatively asexual plant that flowers rarely, revealing remarkably low genomic diversity and DNA methylation levels. Within species, demographic history and the frequency of asexual reproduction jointly determined intra-specific variations of genomic diversity and DNA methylation levels. Genome-wide scans revealed that genes associated with stress adaptations, flowering and embryogenesis were under positive selection. These data are consistent with the hypothesize that natural selection can shape the evolution of asexuality during habitat expansions, which alters genomic and epigenomic diversity levels.

Genome-wide DNA methylation dynamics following recent polyploidy in the allotetraploid Tragopogon miscellus (Asteraceae)

Polyploidy is an important evolutionary force, yet epigenetic mechanisms, such as DNA methylation, that regulate genome-wide expression of duplicated genes remain largely unknown. Here, we use Tragopogon (Asteraceae) as a model system to discover patterns and temporal dynamics of DNA methylation in recently formed polyploids. The naturally occurring allotetraploid Tragopogon miscellus formed in the last 95-100 yr from parental diploids Tragopogon dubius and T. pratensis. We profiled the DNA methylomes of these three species using whole-genome bisulfite sequencing. Genome-wide methylation levels in T. miscellus were intermediate between its diploid parents. However, nonadditive CG and CHG methylation occurred in transposable elements (TEs), with variation among TE types. Most differentially methylated regions (DMRs) showed parental legacy, but some novel DMRs were detected in the polyploid. Differentially methylated genes (DMGs) were also identified and characterized. This study provides the first assessment of both overall and locus-specific patterns of DNA methylation in a recent natural allopolyploid and shows that novel methylation variants can be generated rapidly after polyploid formation. Together, these results demonstrate that mechanisms to regulate duplicate gene expression may arise soon after allopolyploid formation and that these mechanisms vary among genes.

Sulfonamide-induced DNA hypomethylation disturbed sugar metabolism in rice (Oryza sativa L.)

DNA methylation is well-accepted as a bridge to unravel the complex interplay between genome and environmental exposures, and its alteration regulated the cellular metabolic responses towards pollutants. However, the mechanism underlying site-specific aberrant DNA methylation and metabolic disorders under pollutant stresses remained elusive. Herein, the multilevel omics interferences of sulfonamides (i.e., sulfadiazine and sulfamerazine), a group of antibiotics pervasive in farmland soils, towards rice in 14 days of 1 mg/L hydroponic exposure were systematically evaluated. Metabolome and transcriptome analyses showed that 57.1-71.4 % of mono- and disaccharides were accumulated, and the differentially expressed genes were involved in the promotion of sugar hydrolysis, as well as the detoxification of sulfonamides. Most differentially methylated regions (DMRs) were hypomethylated ones (accounting for 87-95 %), and 92 % of which were located in the CHH context (H = A, C, or T base). KEGG enrichment analysis revealed that CHH-DMRs in the promoter regions were enriched in sugar metabolism. To reveal the significant hypomethylation of CHH, multi-spectroscopic and thermodynamic approaches, combined with molecular simulation were conducted to investigate the molecular interaction between sulfonamides and DNA in different sequence contexts, and the result demonstrated that sulfonamides would insert into the minor grooves of DNA, and exhibited a stronger affinity with the CHH contexts of DNA compared to CG or CHG contexts. Computational modeling of DNA 3D structures further confirmed that the binding led to a pitch increase of 0.1 Å and a 3.8° decrease in the twist angle of DNA in the CHH context. This specific interaction and the downregulation of methyltransferase CMT2 (log2FC = -4.04) inhibited the DNA methylation. These results indicated that DNA methylation-based assessment was useful for metabolic toxicity prediction and health risk assessment.

Genome-wide characterization of DNA methyltransferase family genes implies GhDMT6 improving tolerance of salt and drought on cotton

BACKGROUND

DNA methylation is an important epigenetic mode of genomic DNA modification and plays a vital role in maintaining epigenetic content and regulating gene expression. Cytosine-5 DNA methyltransferase (C5-MTase) are the key enzymes in the process of DNA methylation. However, there is no systematic analysis of the C5-MTase in cotton so far, and the function of DNMT2 genes has not been studied.

METHODS

In this study, the whole genome of cotton C5-MTase coding genes was identified and analyzed using a bioinformatics method based on information from the cotton genome, and the function of GhDMT6 was further validated by VIGS experiments and subcellular localization analysis.

RESULTS

33 C5-MTases were identified from three cotton genomes, and were divided into four subfamilies by systematic evolutionary analysis. After the protein domain alignment of C5-MTases in cotton, 6 highly conserved motifs were found in the C-terminus of 33 proteins involved in methylation modification, which indicated that C5-MTases had a basic catalytic methylation function. These proteins were divided into four classes based on the N-terminal difference, of which DNMT2 lacks the N-terminal regulatory domain. The expression of C5-MTases in different parts of cotton was different under different stress treatments, which indicated the functional diversity of cotton C5-MTase gene family. Among the C5-MTases, the GhDMT6 had a obvious up-regulated expression. After silencing GhDMT6 with VIGS, the phenotype of cotton seedlings under different stress treatments showed a significant difference. Compared with cotton seedlings that did not silence GhDMT6, cotton seedlings silencing GhDMT6 showed significant stress resistance.

CONCLUSION

The results show that C5-MTases plays an important role in cotton stress response, which is beneficial to further explore the function of DNMT2 subfamily genes.

Systematic Analysis of DNA Demethylase Gene Families in Foxtail Millet (Setaria italica L.) and Their Expression Variations after Abiotic Stresses

DNA methylation is a highly conserved epigenetic modification involved in many biological processes, including growth and development, stress response, and secondary metabolism. DNA demethylase (DNA-deMTase) genes have been identified in some plant species; however, there are no reports on the identification and analysis of DNA-deMTase genes in Foxtail millet (Setaria italica L.). In this study, seven DNA-deMTases were identified in S. italica. These DNA-deMTase genes were divided into four subfamilies (DML5, DML4, DML3, and ROS1) by phylogenetic and gene structure analysis. Further analysis shows that the physical and chemical properties of these DNA-deMTases proteins are similar, contain the typical conserved domains of ENCO3c and are located in the nucleus. Furthermore, multiple cis-acting elements were observed in DNA-deMTases, including light responsiveness, phytohormone responsiveness, stress responsiveness, and elements related to plant growth and development. The DNA-deMTase genes are expressed in all tissues detected with certain tissue specificity. Then, we investigated the abundance of DNA-deMTase transcripts under abiotic stresses (cold, drought, salt, ABA, and MeJA). The results showed that different genes of DNA-deMTases were involved in the regulation of different abiotic stresses. In total, our findings will provide a basis for the roles of DNA-deMTase in response to abiotic stress.

Ultrasensitive amplification-free quantification of a methyl CpG-rich cancer biomarker by single-molecule kinetic fingerprinting

The most well-studied epigenetic marker in humans is the 5-methyl modification of cytosine in DNA, which has great potential as a disease biomarker in liquid biopsies of cell-free DNA. Currently, quantification of DNA methylation relies heavily on bisulfite conversion followed by PCR amplification and NGS or microarray analysis. PCR is subject to potential bias in differential amplification of bisulfite-converted methylated versus unmethylated sequences. Here, we combine bisulfite conversion with single-molecule kinetic fingerprinting to develop an amplification-free assay for DNA methylation at the branched-chain amino acid transaminase 1 (BCAT1) promoter. Our assay selectively responds to methylated sequences with a limit of detection below 1 fM and a specificity of 99.9999%. Evaluating complex genomic DNA matrices, we reliably distinguish 2-5% DNA methylation at the BCAT1 promoter in whole blood DNA from completely unmethylated whole-genome amplified DNA. Taken together, these results demonstrate the feasibility and sensitivity of our amplification-free, single-molecule quantification approach to improve the early detection of methylated cancer DNA biomarkers.

The unusual predominance of maintenance DNA methylation in Spirodela polyrhiza

Duckweeds are among the fastest reproducing plants, able to clonally divide at exponential rates. However, the genetic and epigenetic impact of clonality on plant genomes is poorly understood. 5-methylcytosine (5mC) is a modified base often described as necessary for the proper regulation of certain genes and transposons and for the maintenance of genome integrity in plants. However, the extent of this dogma is limited by the current phylogenetic sampling of land plant species diversity. Here we analyzed DNA methylomes, small RNAs, mRNA-seq, and H3K9me2 histone modification for Spirodela polyrhiza. S. polyrhiza has lost highly conserved genes involved in de novo methylation of DNA at sites often associated with repetitive DNA, and within genes, however, symmetrical DNA methylation and heterochromatin are maintained during cell division at certain transposons and repeats. Consequently, small RNAs that normally guide methylation to silence repetitive DNA like retrotransposons are diminished. Despite the loss of a highly conserved methylation pathway, and the reduction of small RNAs that normally target repetitive DNA, transposons have not proliferated in the genome, perhaps due in part to the rapid, clonal growth lifestyle of duckweeds.

Regulation and function of transposable elements in cancer genomes

Over half of human genomic DNA is composed of repetitive sequences generated throughout evolution by prolific mobile genetic parasites called transposable elements (TEs). Long disregarded as "junk" or "selfish" DNA, TEs are increasingly recognized as formative elements in genome evolution, wired intimately into the structure and function of the human genome. Advances in sequencing technologies and computational methods have ushered in an era of unprecedented insight into how TE activity impacts human biology in health and disease. Here we discuss the current views on how TEs have shaped the regulatory landscape of the human genome, how TE activity is implicated in human cancers, and how recent findings motivate novel strategies to leverage TE activity for improved cancer therapy. Given the crucial role of methodological advances in TE biology, we pair our conceptual discussions with an in-depth review of the inherent technical challenges in studying repeats, specifically related to structural variation, expression analyses, and chromatin regulation. Lastly, we provide a catalog of existing and emerging assays and bioinformatic software that altogether are enabling the most sophisticated and comprehensive investigations yet into the regulation and function of interspersed repeats in cancer genomes.

HBCR_DMR: A Hybrid Method Based on Beta-Binomial Bayesian Hierarchical Model and Combination of Ranking Method to Detect Differential Methylation Regions in Bisulfite Sequencing Data

DNA methylation is a key epigenetic modification involved in gene regulation, contributing to both physiological and pathological conditions. For a more profound comprehension, it is essential to conduct a precise comparison of DNA methylation patterns between sample groups that represent distinct statuses. Analysis of differentially methylated regions (DMRs) using computational approaches can help uncover the precise relationships between these phenomena. This paper describes a hybrid model that combines the beta-binomial Bayesian hierarchical model with a combination of ranking methods known as HBCR_DMR. During the initial phase, we model the actual methylation proportions of the CpG sites (CpGs) within the replicates. This modeling is achieved through beta-binomial distribution, with parameters set by a group mean and a dispersion parameter. During the second stage, we establish the selection of distinguishing CpG sites based on their methylation status, employing multiple ranking techniques. Finally, we combine the ranking lists of differentially methylated CpG sites through a voting system. Our analyses, encompassing simulations and real data, reveal outstanding performance metrics, including a sensitivity of 0.72, specificity of 0.89, and an F1 score of 0.76, yielding an overall accuracy of 0.82 and an AUC of 0.94. These findings underscore HBCR_DMR's robust capacity to distinguish methylated regions, confirming its utility as a valuable tool for DNA methylation analysis.

Predicting locus-specific DNA methylation levels in cancer and paracancer tissues

Aim: To predict base-resolution DNA methylation in cancerous and paracancerous tissues. Material & methods: We collected six cancer DNA methylation datasets from The Cancer Genome Atlas and five cancer datasets from Gene Expression Omnibus and established machine learning models using paired cancerous and paracancerous tissues. Tenfold cross-validation and independent validation were performed to demonstrate the effectiveness of the proposed method. Results: The developed cross-tissue prediction models can substantially increase the accuracy at more than 68% of CpG sites and contribute to enhancing the statistical power of differential methylation analyses. An XGBoost model leveraging multiple correlating CpGs may elevate the prediction accuracy. Conclusion: This study provides a powerful tool for DNA methylation analysis and has the potential to gain new insights into cancer research from epigenetics.

Lambda3: homology search for protein, nucleotide, and bisulfite-converted sequences

MOTIVATION

Local alignments of query sequences in large databases represent a core part of metagenomic studies and facilitate homology search. Following the development of NCBI Blast, many applications aimed to provide faster and equally sensitive local alignment frameworks. Most applications focus on protein alignments, while only few also facilitate DNA-based searches. None of the established programs allow searching DNA sequences from bisulfite sequencing experiments commonly used for DNA methylation profiling, for which specific alignment strategies need to be implemented.

RESULTS

Here, we introduce Lambda3, a new version of the local alignment application Lambda. Lambda3 is the first solution that enables the search of protein, nucleotide as well as bisulfite-converted nucleotide query sequences. Its protein mode achieves comparable performance to that of the highly optimized protein alignment application Diamond, while the nucleotide mode consistently outperforms established local nucleotide aligners. Combined, Lambda3 presents a universal local alignment framework that enables fast and sensitive homology searches for a wide range of use-cases.

AVAILABILITY AND IMPLEMENTATION

Lambda3 is free and open-source software publicly available at https://github.com/seqan/lambda/.

Epigenetic variation in early and late flowering plants of the rubber-producing Russian dandelion Taraxacum koksaghyz provides insights into the regulation of flowering time

The Russian dandelion (Taraxacum koksaghyz) grows in temperate zones and produces large amounts of poly(cis-1,4-isoprene) in its roots, making it an attractive alternative source of natural rubber. Most T. koksaghyz plants require vernalization to trigger flower development, whereas early flowering varieties that have lost their vernalization dependence are more suitable for breeding and domestication. To provide insight into the regulation of flowering time in T. koksaghyz, we induced epigenetic variation by in vitro cultivation and applied epigenomic and transcriptomic analysis to the resulting early flowering plants and late flowering controls, allowing us to identify differences in methylation patterns and gene expression that correlated with flowering. This led to the identification of candidate genes homologous to vernalization and photoperiodism response genes in other plants, as well as epigenetic modifications that may contribute to the control of flower development. Some of the candidate genes were homologous to known floral regulators, including those that directly or indirectly regulate the major flowering control gene FT. Our atlas of genes can be used as a starting point to investigate mechanisms that control flowering time in T. koksaghyz in greater detail and to develop new breeding varieties that are more suited to domestication.

DNA methylation: an emerging paradigm of gene regulation under drought stress in plants

Drought is an enormous threat to global crop production. In order to ensure food security for the burgeoning population, we must develop drought tolerant crop varieties. This necessitates the identification of drought-responsive genes and understanding the mechanisms involved in their regulation. DNA methylation is a widely studied mechanism of epigenetic regulation of gene expression, which is known to play vital role in conferring tolerance to various biotic and abiotic stress factors. The recent advances in next-generation sequencing (NGS) technologies, has allowed unprecedented access to genome-wide methylation marks, with single base resolution. The most important roles of DNA methylation have been studied in terms of gene body methylation (gbM), which is associated with regulation of both transcript abundance and its stability. The availability of mutants for the various genes encoding enzymes involved in methylation of DNA has allowed ascertainment of the biological significance of methylation. Even though a vast number of reports have emerged in the recent past, where both genome-wide methylation landscape and locus specific changes in DNA methylation have been studied, a conclusive picture with regards to the biological role of DNA methylation is still lacking. Compounding this, is the lack of sufficient evidence supporting the heritability of these epigenetic changes. Amongst the various epigenetic variations, the DNA methylation changes are observed to be the most stable. This review describes the drought-induced changes in DNA methylation identified across different plant species. We also briefly describe the stress memory contributed by these changes. The identification of heritable, drought-induced methylation marks would broaden the scope of crop improvement in the future.

Exploring noncoding variants in genetic diseases: from detection to functional insights

Previous studies on genetic diseases predominantly focused on protein-coding variations, overlooking the vast noncoding regions in the human genome. The development of high-throughput sequencing technologies and functional genomics tools has enabled the systematic identification of functional noncoding variants. These variants can impact gene expression, regulation, and chromatin conformation, thereby contributing to disease pathogenesis. Understanding the mechanisms that underlie the impact of noncoding variants on genetic diseases is indispensable for the development of precisely targeted therapies and the implementation of personalized medicine strategies. The intricacies of noncoding regions introduce a multitude of challenges and research opportunities. In this review, we introduce a spectrum of noncoding variants involved in genetic diseases, along with research strategies and advanced technologies for their precise identification and in-depth understanding of the complexity of the noncoding genome. We will delve into the research challenges and propose potential solutions for unraveling the genetic basis of rare and complex diseases.

scENCORE: leveraging single-cell epigenetic data to predict chromatin conformation using graph embedding

Dynamic compartmentalization of eukaryotic DNA into active and repressed states enables diverse transcriptional programs to arise from a single genetic blueprint, whereas its dysregulation can be strongly linked to a broad spectrum of diseases. While single-cell Hi-C experiments allow for chromosome conformation profiling across many cells, they are still expensive and not widely available for most labs. Here, we propose an alternate approach, scENCORE, to computationally reconstruct chromatin compartments from the more affordable and widely accessible single-cell epigenetic data. First, scENCORE constructs a long-range epigenetic correlation graph to mimic chromatin interaction frequencies, where nodes and edges represent genome bins and their correlations. Then, it learns the node embeddings to cluster genome regions into A/B compartments and aligns different graphs to quantify chromatin conformation changes across conditions. Benchmarking using cell-type-matched Hi-C experiments demonstrates that scENCORE can robustly reconstruct A/B compartments in a cell-type-specific manner. Furthermore, our chromatin confirmation switching studies highlight substantial compartment-switching events that may introduce substantial regulatory and transcriptional changes in psychiatric disease. In summary, scENCORE allows accurate and cost-effective A/B compartment reconstruction to delineate higher-order chromatin structure heterogeneity in complex tissues.

epidecodeR: a functional exploration tool for epigenetic and epitranscriptomic regulation

Recent technological advances in sequencing DNA and RNA modifications using high-throughput platforms have generated vast epigenomic and epitranscriptomic datasets whose power in transforming life science is yet fully unleashed. Currently available in silico methods have facilitated the identification, positioning and quantitative comparisons of individual modification sites. However, the essential challenge to link specific 'epi-marks' to gene expression in the particular context of cellular and biological processes is unmet. To fast-track exploration, we generated epidecodeR implemented in R, which allows biologists to quickly survey whether an epigenomic or epitranscriptomic status of their interest potentially influences gene expression responses. The evaluation is based on the cumulative distribution function and the statistical significance in differential expression of genes grouped by the number of 'epi-marks'. This tool proves useful in predicting the role of H3K9ac and H3K27ac in associated gene expression after knocking down deacetylases FAM60A and SDS3 and N6-methyl-adenosine-associated gene expression after knocking out the reader proteins. We further used epidecodeR to explore the effectiveness of demethylase FTO inhibitors and histone-associated modifications in drug abuse in animals. epidecodeR is available for downloading as an R package at https://bioconductor.riken.jp/packages/3.13/bioc/html/epidecodeR.html.

Tomato plant response to heat stress: a focus on candidate genes for yield-related traits

Climate change and global warming represent the main threats for many agricultural crops. Tomato is one of the most extensively grown and consumed horticultural products and can survive in a wide range of climatic conditions. However, high temperatures negatively affect both vegetative growth and reproductive processes, resulting in losses of yield and fruit quality traits. Researchers have employed different parameters to evaluate the heat stress tolerance, including evaluation of leaf- (stomatal conductance, net photosynthetic rate, Fv/Fm), flower- (inflorescence number, flower number, stigma exertion), pollen-related traits (pollen germination and viability, pollen tube growth) and fruit yield per plant. Moreover, several authors have gone even further, trying to understand the plants molecular response mechanisms to this stress. The present review focused on the tomato molecular response to heat stress during the reproductive stage, since the increase of temperatures above the optimum usually occurs late in the growing tomato season. Reproductive-related traits directly affects the final yield and are regulated by several genes such as transcriptional factors, heat shock proteins, genes related to flower, flowering, pollen and fruit set, and epigenetic mechanisms involving DNA methylation, histone modification, chromatin remodelling and non-coding RNAs. We provided a detailed list of these genes and their function under high temperature conditions in defining the final yield with the aim to summarize the recent findings and pose the attention on candidate genes that could prompt on the selection and constitution of new thermotolerant tomato plant genotypes able to face this abiotic challenge.

Frequency of cytosine methylation in the adjacent regions of soybean retrotransposon SORE-1 depends on chromosomal location

Mobilization of transposable elements (TEs) is suppressed by epigenetic mechanisms involving cytosine methylation. However, few studies have focused on clarifying relationships between epigenetic influences of TEs on the adjacent DNA regions and time after insertion of TEs into the genome and/or their chromosomal location. Here we addressed these issues using soybean retrotransposon SORE-1. We analyzed SORE-1, inserted in exon 1 of the GmphyA2 gene, one of the newest insertions in this family so far identified. Cytosine methylation was detected in this element but was barely present in the adjacent regions. These results were correlated, respectively, with the presence and absence of the production of short interfering RNAs. Cytosine methylation profiles of 74 SORE-1 elements in the Williams 82 reference genome indicated that methylation frequency in the adjacent regions of SORE-1 was profoundly higher in pericentromeric regions than in euchromatic chromosome arms and was only weakly correlated with the length of time after insertion into the genome. Notably, the higher level of methylation in the 5' adjacent regions of SORE-1 coincided with the presence of repetitive elements in pericentromeric regions. Together, these results suggest that epigenetic influence of SORE-1 on the adjacent regions is influenced by its location on the chromosome.

Knocking down a DNA demethylase gene affects potato plant defense against a specialist insect herbivore

DNA demethylase (DML) is involved in plant development and responses to biotic and abiotic stresses; however, its role in plant-herbivore interaction remains elusive. Here, we found that herbivory by the potato tuber moth, Phthorimaea operculella, rapidly induced the genome-wide DNA methylation and accumulation of DML gene transcripts in potato plants. Herbivory induction of DML transcripts was suppressed in jasmonate-deficient plants, whereas exogenous application of methyl jasmonate (MeJA) improved DML transcripts, indicating that the induction of DML transcripts by herbivory is associated with jasmonate signaling. Moreover, P. operculella larvae grew heavier on DML gene (StDML2) knockdown plants than on wild-type plants, and the decreased biosynthesis of jasmonates in the former may be responsible for this difference, since the larvae feeding on these two genotypes supplemented with MeJA showed similar growth. In addition, P. operculella adult moths preferred to oviposit on StDML2 knockdown plants than on wild-type plants, which was associated with the reduced emission of β-caryophyllene in the former. In addition, supplementing β-caryophyllene to these two genotypes further disrupted moths' oviposit choice preference for them. Interestingly, in StDML2 knockdown plants, hypermethylation was found at the promoter regions for the key genes StAOS and StAOC in the jasmonate biosynthetic pathway, as well as for the key gene StTPS12 in β-caryophyllene production. Our findings suggest that knocking down StDML2 can affect herbivore defense via jasmonate signaling and defense compound production in potato plants.

Phased gap-free genome assembly of octoploid cultivated strawberry illustrates the genetic and epigenetic divergence among subgenomes

The genetic and epigenetic mechanisms underlying the coexistence and coordination of the four diverged subgenomes (ABCD) in octoploid strawberries (Fragaria × ananassa) remains poorly understood. In this study, we have assembled a haplotype-phased gap-free octoploid genome for the strawberry, which allowed us to uncover the sequence, structure, and epigenetic divergences among the subgenomes. The diploid progenitors of the octoploid strawberry, apart from subgenome A (Fragaria vesca), have been a subject of public controversy. Phylogenomic analyses revealed a close relationship between diploid species Fragaria iinumae and subgenomes B, C, and D. Subgenome A, closely related to F. vesca, retains the highest number of genes, exhibits the lowest content of transposable elements (TEs), experiences the strongest purifying selection, shows the lowest DNA methylation levels, and displays the highest expression level compared to the other three subgenomes. Transcriptome and DNA methylome analyses revealed that subgenome A-biased genes were enriched in fruit development biological processes. In contrast, although subgenomes B, C, and D contain equivalent amounts of repetitive sequences, they exhibit diverged methylation levels, particularly for TEs located near genes. Taken together, our findings provide valuable insights into the evolutionary patterns of subgenome structure, divergence and epigenetic dynamics in octoploid strawberries, which could be utilized in strawberry genetics and breeding research.

Ambiguous genes due to aligners and their impact on RNA-seq data analysis

The main scope of the study is ambiguous genes, i.e. genes whose expression is difficult to estimate from the data produced by next-generation sequencing technologies. We focused on the RNA sequencing (RNA-Seq) type of experiment performed on the Illumina platform. It is crucial to identify such genes and understand the cause of their difficulty, as these genes may be involved in some diseases. By giving misleading results, they could contribute to a misunderstanding of the cause of certain diseases, which could lead to inappropriate treatment. We thought that the ambiguous genes would be difficult to map because of their complex structure. So we looked at RNA-seq analysis using different mappers to find genes that would have different measurements from the aligners. We were able to identify such genes using a generalized linear model with two factors: mappers and groups introduced by the experiment. A large proportion of ambiguous genes are pseudogenes. High sequence similarity of pseudogenes to functional genes may indicate problems in alignment procedures. In addition, predictive analysis verified the performance of difficult genes in classification. The effectiveness of classifying samples into specific groups was compared, including the expression of difficult and not difficult genes as covariates. In almost all cases considered, ambiguous genes have less predictive power.

Single-cell bisulfite-free 5mC and 5hmC sequencing with high sensitivity and scalability

Existing single-cell bisulfite-based DNA methylation analysis is limited by low DNA recovery, and the measurement of 5hmC at single-base resolution remains challenging. Here, we present a bisulfite-free single-cell whole-genome 5mC and 5hmC profiling technique, named Cabernet, which can characterize 5mC and 5hmC at single-base resolution with high genomic coverage. Cabernet utilizes Tn5 transposome for DNA fragmentation, which enables the discrimination between different alleles for measuring hemi-methylation status. Using Cabernet, we revealed the 5mC, hemi-5mC and 5hmC dynamics during early mouse embryo development, uncovering genomic regions exclusively governed by active or passive demethylation. We show that hemi-methylation status can be used to distinguish between pre- and post-replication cells, enabling more efficient cell grouping when integrated with 5mC profiles. The property of Tn5 naturally enables Cabernet to achieve high-throughput single-cell methylome profiling, where we probed mouse cortical neurons and embryonic day 7.5 (E7.5) embryos, and constructed the library for thousands of single cells at high efficiency, demonstrating its potential for analyzing complex tissues at substantially low cost. Together, we present a way of high-throughput methylome and hydroxymethylome detection at single-cell resolution, enabling efficient analysis of the epigenetic status of biological systems with complicated nature such as neurons and cancer cells.

Cost-effective DNA methylation profiling by FML-seq

Current methods for profiling DNA methylation require costly reagents, sequencing, and labor time. We introduce fragmentation at methylated loci and sequencing (FML-seq), a sequencing library protocol that greatly reduces all these costs. Relative to other techniques tested on the same human cell lines, FML-seq produces similar measurements of absolute and differential cytosine methylation at a fraction of the price. FML-seq enables inexpensive, high-throughput experimental designs for large-scale epigenetics research projects.

DNA methylation differences between stick insect ecotypes

Epigenetic mechanisms, such as DNA methylation, can influence gene regulation and affect phenotypic variation, raising the possibility that they contribute to ecological adaptation. Beginning to address this issue requires high-resolution sequencing studies of natural populations to pinpoint epigenetic regions of potential ecological and evolutionary significance. However, such studies are still relatively uncommon, especially in insects, and are mainly restricted to a few model organisms. Here, we characterize patterns of DNA methylation for natural populations of Timema cristinae adapted to two host plant species (i.e. ecotypes). By integrating results from sequencing of whole transcriptomes, genomes and methylomes, we investigate whether environmental, host and genetic differences of these stick insects are associated with methylation levels of cytosine nucleotides in the CpG context. We report an overall genome-wide methylation level for T. cristinae of ~14%, with methylation being enriched in gene bodies and impoverished in repetitive elements. Genome-wide DNA methylation variation was strongly positively correlated with genetic distance (relatedness), but also exhibited significant host-plant effects. Using methylome-environment association analysis, we pinpointed specific genomic regions that are differentially methylated between ecotypes, with these regions being enriched for genes with functions in membrane processes. The observed association between methylation variation and genetic relatedness, and with the ecologically important variable of host plant, suggests a potential role for epigenetic modification in T. cristinae adaptation. To substantiate such adaptive significance, future studies could test whether methylation can be transmitted across generations and the extent to which it responds to experimental manipulation in field and laboratory studies.

DNA conformational dynamics in the context-dependent non-CG CHH methylation by plant methyltransferase DRM2

DNA methylation provides an important epigenetic mechanism that critically regulates gene expression, genome imprinting, and retrotransposon silencing. In plants, DNA methylation is prevalent not only in a CG dinucleotide context but also in non-CG contexts, namely CHG and CHH (H = C, T, or A) methylation. It has been established that plant non-CG DNA methylation is highly context dependent, with the +1- and +2-flanking sequences enriched with A/T nucleotides. How DNA sequence, conformation, and dynamics influence non-CG methylation remains elusive. Here, we report structural and biochemical characterizations of the intrinsic substrate preference of DOMAINS REARRANGED METHYLTRANSFERASE 2 (DRM2), a plant DNA methyltransferase responsible for establishing all cytosine methylation and maintaining CHH methylation. Among nine CHH motifs, the DRM2 methyltransferase (MTase) domain shows marked substrate preference toward CWW (W = A or T) motifs, correlating well with their relative abundance in planta. Furthermore, we report the crystal structure of DRM2 MTase in complex with a DNA duplex containing a flexible TpA base step at the +1/+2-flanking sites of the target nucleotide. Comparative structural analysis of the DRM2-DNA complexes provides a mechanism by which flanking nucleotide composition impacts DRM2-mediated DNA methylation. Furthermore, the flexibility of the TpA step gives rise to two alternative DNA conformations, resulting in different interactions with DRM2 and consequently temperature-dependent shift of the substrate preference of DRM2. Together, this study provides insights into how the interplay between the conformational dynamics of DNA and temperature as an environmental factor contributes to the context-dependent CHH methylation by DRM2.

Estimating genome-wide DNA methylation heterogeneity with methylation patterns

BACKGROUND

In a heterogeneous population of cells, individual cells can behave differently and respond variably to the environment. This cellular diversity can be assessed by measuring DNA methylation patterns. The loci with variable methylation patterns are informative of cellular heterogeneity and may serve as biomarkers of diseases and developmental progression. Cell-to-cell methylation heterogeneity can be evaluated through single-cell methylomes or computational techniques for pooled cells. However, the feasibility and performance of these approaches to precisely estimate methylation heterogeneity require further assessment.

RESULTS

Here, we proposed model-based methods adopted from a mathematical framework originally from biodiversity, to estimate genome-wide DNA methylation heterogeneity. We evaluated the performance of our models and the existing methods with feature comparison, and tested on both synthetic datasets and real data. Overall, our methods have demonstrated advantages over others because of their better correlation with the actual heterogeneity. We also demonstrated that methylation heterogeneity offers an additional layer of biological information distinct from the conventional methylation level. In the case studies, we showed that distinct profiles of methylation heterogeneity in CG and non-CG methylation can predict the regulatory roles between genomic elements in Arabidopsis. This opens up a new direction for plant epigenomics. Finally, we demonstrated that our score might be able to identify loci in human cancer samples as putative biomarkers for early cancer detection.

CONCLUSIONS

We adopted the mathematical framework from biodiversity into three model-based methods for analyzing genome-wide DNA methylation heterogeneity to monitor cellular heterogeneity. Our methods, namely MeH, have been implemented, evaluated with existing methods, and are open to the research community.

Mining the Utricularia gibba genome for insulator-like elements for genetic engineering

Introduction

Gene expression is often controlled via cis-regulatory elements (CREs) that modulate the production of transcripts. For multi-gene genetic engineering and synthetic biology, precise control of transcription is crucial, both to insulate the transgenes from unwanted native regulation and to prevent readthrough or cross-regulation of transgenes within a multi-gene cassette. To prevent this activity, insulator-like elements, more properly referred to as transcriptional blockers, could be inserted to separate the transgenes so that they are independently regulated. However, only a few validated insulator-like elements are available for plants, and they tend to be larger than ideal.

Methods

To identify additional potential insulator-like sequences, we conducted a genome-wide analysis of Utricularia gibba (humped bladderwort), one of the smallest known plant genomes, with genes that are naturally close together. The 10 best insulator-like candidates were evaluated in vivo for insulator-like activity.

Results

We identified a total of 4,656 intergenic regions with expression profiles suggesting insulator-like activity. Comparisons of these regions across 45 other plant species (representing Monocots, Asterids, and Rosids) show low levels of syntenic conservation of these regions. Genome-wide analysis of unmethylated regions (UMRs) indicates ~87% of the targeted regions are unmethylated; however, interpretation of this is complicated because U. gibba has remarkably low levels of methylation across the genome, so that large UMRs frequently extend over multiple genes and intergenic spaces. We also could not identify any conserved motifs among our selected intergenic regions or shared with existing insulator-like elements for plants. Despite this lack of conservation, however, testing of 10 selected intergenic regions for insulator-like activity found two elements on par with a previously published element (EXOB) while being significantly smaller.

Discussion

Given the small number of insulator-like elements currently available for plants, our results make a significant addition to available tools. The high hit rate (2 out of 10) also implies that more useful sequences are likely present in our selected intergenic regions; additional validation work will be required to identify which will be most useful for plant genetic engineering.

Methylome Response to Proteasome Inhibition by Pseudomonas syringae Virulence Factor Syringolin A

DNA methylation is an important epigenetic mark required for proper gene expression and silencing of transposable elements. DNA methylation patterns can be modified by environmental factors such as pathogen infection, in which modification of DNA methylation can be associated with plant resistance. To counter the plant defense pathways, pathogens produce effector molecules, several of which act as proteasome inhibitors. Here, we investigated the effect of proteasome inhibition by the bacterial virulence factor syringolin A (SylA) on genome-wide DNA methylation. We show that SylA treatment results in an increase of DNA methylation at centromeric and pericentromeric regions of Arabidopsis chromosomes. We identify several CHH differentially methylated regions (DMRs) that are enriched in the proximity of transcriptional start sites. SylA treatment does not result in significant changes in small RNA composition. However, significant changes in genome transcriptional activity can be observed, including a strong upregulation of resistance genes that are located on chromosomal arms. We hypothesize that DNA methylation changes could be linked to the upregulation of some atypical members of the de novo DNA methylation pathway, namely AGO3, AGO9, and DRM1. Our data suggests that modification of genome-wide DNA methylation resulting from an inhibition of the proteasome by bacterial effectors could be part of an epi-genomic arms race against pathogens. [Formula: see text] Copyright © 2023 The Author(s). This is an open access article distributed under the CC BY-NC-ND 4.0 International license.

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.

Methylomics and cancer: the current state of methylation profiling and marker development for clinical care

Epigenetic modifications have long been recognized as an essential level in transcriptional regulation linking behavior and environmental conditions or stimuli with biological processes and disease development. Among them, methylation is the most abundant of these reversible epigenetic marks, predominantly occurring on DNA, RNA, and histones. Methylation modification is intimately involved in regulating gene transcription and cell differentiation, while aberrant methylation status has been linked with cancer development in several malignancies. Early detection and precise restoration of dysregulated methylation form the basis for several epigenetics-based therapeutic strategies. In this review, we summarize the current basic understanding of the regulation and mechanisms responsible for methylation modification and cover several cutting-edge research techniques for detecting methylation across the genome and transcriptome. We then explore recent advances in clinical diagnostic applications of methylation markers of various cancers and address the current state and future prospects of methylation modifications in therapies for different diseases, especially comparing pharmacological methylase/demethylase inhibitors with the CRISPRoff/on methylation editing systems. This review thus provides a resource for understanding the emerging role of epigenetic methylation in cancer, the use of methylation-based biomarkers in cancer detection, and novel methylation-targeted drugs.

Widespread effects of DNA methylation and intra-motif dependencies revealed by novel transcription factor binding models

Several studies suggested that transcription factor (TF) binding to DNA may be impaired or enhanced by DNA methylation. We present MeDeMo, a toolbox for TF motif analysis that combines information about DNA methylation with models capturing intra-motif dependencies. In a large-scale study using ChIP-seq data for 335 TFs, we identify novel TFs that show a binding behaviour associated with DNA methylation. Overall, we find that the presence of CpG methylation decreases the likelihood of binding for the majority of methylation-associated TFs. For a considerable subset of TFs, we show that intra-motif dependencies are pivotal for accurately modelling the impact of DNA methylation on TF binding. We illustrate that the novel methylation-aware TF binding models allow to predict differential ChIP-seq peaks and improve the genome-wide analysis of TF binding. Our work indicates that simplistic models that neglect the effect of DNA methylation on DNA binding may lead to systematic underperformance for methylation-associated TFs.

Omics Views of Mechanisms for Cell Fate Determination in Early Mammalian Development

During mammalian preimplantation development, a totipotent zygote undergoes several cell cleavages and two rounds of cell fate determination, ultimately forming a mature blastocyst. Along with compaction, the establishment of apicobasal cell polarity breaks the symmetry of an embryo and guides subsequent cell fate choice. Although the lineage segregation of the inner cell mass (ICM) and trophectoderm (TE) is the first symbol of cell differentiation, several molecules have been shown to bias the early cell fate through their inter-cellular variations at much earlier stages, including the 2- and 4-cell stages. The underlying mechanisms of early cell fate determination have long been an important research topic. In this review, we summarize the molecular events that occur during early embryogenesis, as well as the current understanding of their regulatory roles in cell fate decisions. Moreover, as powerful tools for early embryogenesis research, single-cell omics techniques have been applied to both mouse and human preimplantation embryos and have contributed to the discovery of cell fate regulators. Here, we summarize their applications in the research of preimplantation embryos, and provide new insights and perspectives on cell fate regulation.

Stochasticity in gene body methylation

Gene body methylation (gbM) is a widely conserved epigenetic feature of plant genomes. Efforts to delineate the mechanisms by which gbM contributes to transcriptional regulation remain largely inconclusive, and its evolutionary significance continues to be debated. Curiously, although steady-state gbM levels are remarkably stable across mitotic and meiotic cell divisions, the methylation status of individual CG dinucleotides in gbM genes is highly stochastic. How can these two seemingly contradictory observations be reconciled? Here, we discuss how stochastic processes relate to gbM maintenance dynamics. We show that a quantitative understanding of these processes can shed deeper insights into the molecular and evolutionary biology of this enigmatic epigenetic trait.

Flavonoid Biosynthesis Pathway May Indirectly Affect Outcrossing Rate of Cytoplasmic Male-Sterile Lines of Soybean

(1) Background: Cytoplasmic male sterility (CMS) is important for exploiting heterosis. Soybean (Glycine max L.) has a low outcrossing rate that is detrimental for breeding sterile lines and producing hybrid seeds. Therefore, the molecular mechanism controlling the outcrossing rate should be elucidated to increase the outcrossing rate of soybean CMS lines; (2) Methods: The male-sterile soybean lines JLCMS313A (with a high outcrossing rate; HL) and JLCMS226A (with a low outcrossing rate; LL) were used for a combined analysis of the transcriptome (RNA-seq) and the targeted phenol metabolome; (3) Results: The comparison between HL and LL detected 5946 differentially expressed genes (DEGs) and 81 phenolic metabolites. The analysis of the DEGs and differentially abundant phenolic metabolites identified only one common KEGG pathway related to flavonoid biosynthesis. The qRT-PCR expression for eight DEGs was almost consistent with the transcriptome data. The comparison of the cloned coding sequence (CDS) regions of the SUS, FLS, UGT, and F3H genes between HL and LL revealed seven single nucleotide polymorphisms (SNPs) only in the F3H CDS. Moreover, five significant differentially abundant phenolic metabolites between HL and LL were associated with flavonoid metabolic pathways. Finally, on the basis of the SNPs in the F3H CDS, one derived cleaved amplified polymorphic sequence (dCAPS) marker was developed to distinguish between HL and LL soybean lines; (4) Conclusions: The flavonoid biosynthesis pathway may indirectly affect the outcrossing rate of CMS sterile lines in soybean.

Herbivory induced methylation changes in the Lombardy poplar: A comparison of results obtained by epiGBS and WGBS

DNA cytosine methylation is an epigenetic mechanism involved in regulation of plant responses to biotic and abiotic stress and its ability to change can vary with the sequence context in which a cytosine appears (CpG, CHG, CHH, where H = Adenine, Thymine, Cytosine). Quantification of DNA methylation in model plant species is frequently addressed by Whole Genome Bisulfite Sequencing (WGBS), which requires a good-quality reference genome. Reduced Representation Bisulfite Sequencing (RRBS) is a cost-effective potential alternative for ecological research with limited genomic resources and large experimental designs. In this study, we provide for the first time a comprehensive comparison between the outputs of RRBS and WGBS to characterize DNA methylation changes in response to a given environmental factor. In particular, we used epiGBS (recently optimized RRBS) and WGBS to assess global and sequence-specific differential methylation after insect and artificial herbivory in clones of Populus nigra cv. 'italica'. We found that, after any of the two herbivory treatments, global methylation percentage increased in CHH, and the shift was detected as statistically significant only by epiGBS. As regards to loci-specific differential methylation induced by herbivory (cytosines in epiGBS and regions in WGBS), both techniques indicated the specificity of the response elicited by insect and artificial herbivory, together with higher frequency of hypo-methylation in CpG and hyper-methylation in CHH. Methylation changes were mainly found in gene bodies and intergenic regions when present at CpG and CHG and in transposable elements and intergenic regions at CHH context. Thus, epiGBS succeeded to characterize global, genome-wide methylation changes in response to herbivory in the Lombardy poplar. Our results support that epiGBS could be particularly useful in large experimental designs aimed to explore epigenetic changes of non-model plant species in response to multiple environmental factors.

Single-molecule targeted accessibility and methylation sequencing of centromeres, telomeres and rDNAs in Arabidopsis

The short read-length of next-generation sequencing makes it challenging to characterize highly repetitive regions (HRRs) such as centromeres, telomeres and ribosomal DNAs. Based on recent strategies that combined long-read sequencing and exogenous enzymatic labelling of open chromatin, we developed single-molecule targeted accessibility and methylation sequencing (STAM-seq) in plants by further integrating nanopore adaptive sampling to investigate the HRRs in wild-type Arabidopsis and DNA methylation mutants that are defective in CG- or non-CG methylation. We found that CEN180 repeats show higher chromatin accessibility and lower DNA methylation on their forward strand, individual rDNA units show a negative correlation between their DNA methylation and accessibility, and both accessibility and CHH methylation levels are lower at telomere compared to adjacent subtelomeric region. Moreover, DNA methylation-deficient mutants showed increased chromatin accessibility at HRRs, consistent with the role of DNA methylation in maintaining heterochromatic status in plants. STAM-seq can be applied to study accessibility and methylation of repetitive sequences across diverse plant species.

Global effects of identity and aging on the human sperm methylome

BACKGROUND

As the average age of fatherhood increases worldwide, so too does the need for understanding effects of aging in male germline cells. Molecular change, including epigenomic alterations, may impact offspring. Age-associated change to DNA cytosine methylation in the cytosine-guanine (CpG) context is a hallmark of aging tissues, including sperm. Prior studies have led to accurate models that predict a man's age based on specific methylation features in the DNA of sperm, but the relationship between aging and global DNA methylation in sperm remains opaque. Further clarification requires a more complete survey of the methylome with assessment of variability within and between individuals.

RESULTS

We collected sperm methylome data in a longitudinal study of ten healthy fertile men. We used whole-genome bisulfite sequencing of samples collected 10 to 18 years apart from each donor. We found that, overall, variability between donors far exceeds age-associated variation. After controlling for donor identity, we see significant age-dependent genome-wide change to the methylome. Notably, trends of change with age depend on genomic location or annotation, with contrasting signatures that correlate with gene density and proximity to centromeres and promoter regions.

CONCLUSIONS

We uncovered epigenetic signatures that reflect a stable process which begins in early adulthood, progressing steadily through most of the male lifespan, and warrants consideration in any future study of the aging sperm epigenome.

Concepts, mechanisms and implications of long-term epigenetic inheritance

Many modes and mechanisms of epigenetic inheritance have been elucidated in eukaryotes. Most of them are relatively short-term, generally not exceeding one or a few organismal generations. However, emerging evidence indicates that one mechanism, cytosine DNA methylation, can mediate epigenetic inheritance over much longer timescales, which are mostly or completely inaccessible in the laboratory. Here we discuss the evidence for, and mechanisms and implications of, such long-term epigenetic inheritance. We argue that compelling evidence supports the long-term epigenetic inheritance of gene body methylation, at least in the model angiosperm Arabidopsis thaliana, and that variation in such methylation can therefore serve as an epigenetic basis for phenotypic variation in natural populations.

The Progression in Developing Genomic Resources for Crop Improvement

Sequencing technologies have rapidly evolved over the past two decades, and new technologies are being continually developed and commercialized. The emerging sequencing technologies target generating more data with fewer inputs and at lower costs. This has also translated to an increase in the number and type of corresponding applications in genomics besides enhanced computational capacities (both hardware and software). Alongside the evolving DNA sequencing landscape, bioinformatics research teams have also evolved to accommodate the increasingly demanding techniques used to combine and interpret data, leading to many researchers moving from the lab to the computer. The rich history of DNA sequencing has paved the way for new insights and the development of new analysis methods. Understanding and learning from past technologies can help with the progress of future applications. This review focuses on the evolution of sequencing technologies, their significant enabling role in generating plant genome assemblies and downstream applications, and the parallel development of bioinformatics tools and skills, filling the gap in data analysis techniques.

Cost-effective DNA methylation profiling by FML-seq

Current methods for profiling DNA methylation require costly reagents, sequencing, or labor time. We introduce FML-seq, a sequencing library protocol that greatly reduces all these costs. Relative to other techniques tested on the same human cell lines, FML-seq produces similar measurements of absolute and differential cytosine methylation at a fraction of the price. FML-seq enables inexpensive, high-throughput experimental designs for large-scale epigenetics research projects.

MBD3 Regulates Male Germ Cell Division and Sperm Fertility in Arabidopsis thaliana

DNA methylation plays important roles through the methyl-CpG-binding domain (MBD) to realize epigenetic modifications. Thirteen AtMBD proteins have been identified from the Arabidopsis thaliana genome, but the functions of some members are unclear. AtMBD3 was found to be highly expressed in pollen and seeds and it preferably binds methylated CG, CHG, and unmethylated DNA sequences. Then, two mutant alleles at the AtMBD3 locus were obtained in order to further explore its function using CRISPR/Cas9. When compared with 92.17% mature pollen production in the wild type, significantly lower percentages of 84.31% and 78.91% were observed in the mbd3-1 and mbd3-2 mutants, respectively. About 16-21% of pollen from the mbd3 mutants suffered a collapse in reproductive transmission, whereas the other pollen was found to be normal. After pollination, about 16% and 24% of mbd3-1 and mbd3-2 mutant seeds underwent early or late abortion, respectively. Among all the late abortion seeds in mbd3-2 plants, 25% of the abnormal seeds were at the globular stage, 31.25% were at the transition stage, and 43.75% were at the heart stage. A transcriptome analysis of the seeds found 950 upregulated genes and 1128 downregulated genes between wild type and mbd3-2 mutants. Some transcriptional factors involved in embryo development were selected to be expressed, and we found significant differences between wild type and mbd3 mutants, such as WOXs, CUC1, AIB4, and RGL3. Furthermore, we found a gene that is specifically expressed in pollen, named PBL6. PBL6 was found to directly interact with AtMBD3. Our results provide insights into the function of AtMBD3 in plants, especially in sperm fertility.

Multiomics Reveals the Regulatory Mechanisms of Arabidopsis Tissues under Heat Stress

Understanding the mechanisms of responses to high temperatures in Arabidopsis will provide insights into how plants may mitigate heat stress under global climate change. And exploring the interconnections of different modification levels in heat stress response could help us to understand the molecular mechanism of heat stress response in Arabidopsis more comprehensively and precisely. In this paper, we combined multiomics analyses to explore the common heat stress-responsive genes and specific heat-responsive metabolic pathways in Arabidopsis leaf, seedling, and seed tissues. We found that genes such as AT1G54050 play a role in promoting proper protein folding in response to HS (Heat stress). In addition, it was revealed that the binding profile of A1B is altered under elevated temperature conditions. Finally, we also show that two microRNAs, ath-mir156h and ath-mir166b-5p, may be core regulatory molecules in HS. Also elucidated that under HS, plants can regulate specific regulatory mechanisms, such as oxygen levels, by altering the degree of CHH methylation.

Characteristics of Salvia miltiorrhiza methylome and the regulatory mechanism of DNA methylation in tanshinone biosynthesis

Salvia miltiorrhiza is a model medicinal plant with significant economic and medicinal value. Its roots produce a group of diterpenoid lipophilic bioactive components, termed tanshinones. Biosynthesis and regulation of tanshinones has attracted widespread interest. However, the methylome of S. miltiorrhiza has not been analysed and the regulatory mechanism of DNA methylation in tanshinone production is largely unknown. Here we report single-base resolution DNA methylomes from roots and leaves. Comparative analysis revealed differential methylation patterns for CG, CHG, and CHH contexts and the association between DNA methylation and the expression of genes and small RNAs. Lowly methylated genes always had higher expression levels and 24-nucleotide sRNAs could be key players in the RdDM pathway in S. miltiorrhiza. DNA methylation variation analysis showed that CHH methylation contributed mostly to the difference. Go enrichment analysis showed that diterpenoid biosynthetic process was significantly enriched for genes with downstream overlapping with hypoCHHDMR in July_root when comparing with those in March_root. Tanshinone biosynthesis-related enzyme genes, such as DXS2, CMK, IDI1, HMGR2, DXR, MDS, CYP76AH1, 2OGD25, and CYP71D373, were less CHH methylated in gene promoters or downstream regions in roots collected in July than those collected in March. Consistently, gene expression was up-regulated in S. miltiorrhiza roots collected in July compared with March and the treatment of DNA methylation inhibitor 5-azacytidine significantly promoted tanshinone production. It suggests that DNA methylation plays a significant regulatory role in tanshinone biosynthesis in S. miltiorrhiza through changing the levels of CHH methylation in promoters or downstreams of key enzyme genes.