The Effects of DNA Methylation on Cytoplasmic Male Sterility in Sugar Beet

DNA methylation is widely found in higher plants and can control gene expression by regulation without changing the DNA sequence. In this study, the whole-genome methylation map of sugar beet was constructed by WGBS (whole-genome bisulfite sequencing) technology, and the results of WGBS were verified by bisulfite transformation, indicating that the results of WGBS technology were reliable. In addition, 12 differential methylation genes (DMGs) were identified, which were related to carbohydrate and energy metabolism, pollen wall development, and endogenous hormone regulation. Quantitative real-time PCR (qRT-PCR) showed that 75% of DMG expression levels showed negative feedback with methylation level, indicating that DNA methylation can affect gene expression to a certain extent. In addition, we found hypermethylation inhibited gene expression, which laid a foundation for further study on the molecular mechanism of DNA methylation at the epigenetic level in sugar beet male sterility.

Single-cell RNA-seq and single-cell bisulfite-sequencing reveal insights into yak preimplantation embryogenesis

Extensive epigenetic reprogramming occurs during preimplantation embryonic development. However, the impact of DNA methylation in plateau yak preimplantation embryos and how epigenetic reprogramming contributes to transcriptional regulatory networks are unclear. In this study, we quantified gene expression and DNA methylation in oocytes and a series of yak embryos at different developmental stages and at single-cell resolution using single-cell bisulfite-sequencing and RNA-seq. We characterized embryonic genome activation and maternal transcript degradation and mapped epigenetic reprogramming events critical for embryonic development. Through cross-species transcriptome analysis, we identified 31 conserved maternal hub genes and 39 conserved zygotic hub genes, including SIN3A, PRC1, HDAC1/2, and HSPD1. Notably, by combining single-cell DNA methylation and transcriptome analysis, we identified 43 candidate methylation driver genes, such as AURKA, NUSAP1, CENPF, and PLK1, that may be associated with embryonic development. Finally, using functional approaches, we further determined that the epigenetic modifications associated with the histone deacetylases HDAC1/2 are essential for embryonic development and that the deubiquitinating enzyme USP7 may affect embryonic development by regulating DNA methylation. Our data represent an extensive resource on the transcriptional dynamics of yak embryonic development and DNA methylation remodeling, and provide new insights into strategies for the conservation of germplasm resources, as well as a better understanding of mammalian early embryonic development that can be applied to investigate the causes of early developmental disorders.

Age-associated DNA methylation changes in Xenopus frogs

Age-associated changes in DNA methylation have been characterized across various animals, but not yet in amphibians, which are of particular interest because they include widely studied model organisms. In this study, we present clear evidence that the aquatic vertebrate species Xenopus tropicalis displays patterns of age-associated changes in DNA methylation. We have generated whole-genome bisulfite sequencing (WGBS) profiles from skin samples of nine frogs representing young, mature, and old adults and characterized the gene- and chromosome-scale DNA methylation changes with age. Many of the methylation features and changes we observe are consistent with what is known in mammalian species, suggesting that the mechanism of age-related changes is conserved. Moreover, we selected a few thousand age-associated CpG sites to build an assay based on targeted DNA methylation analysis (TBSseq) to expand our findings in future studies involving larger cohorts of individuals. Preliminary results of a pilot TBSeq experiment recapitulate the findings obtained with WGBS setting the basis for the development of an epigenetic clock assay. The results of this study will allow us to leverage the unique resources available for Xenopus to study how DNA methylation relates to other hallmarks of ageing.

Genome-wide DNA methylation sequencing reveals epigenetic features and potential biomarkers of Sjögren syndrome

AIM

Sjögren syndrome (SS) is a slowly progressive, inflammatory, autoimmune disease. The aim of this study was to construct the DNA methylation profiles of whole blood of SS patients and healthy controls (HC), and to explore the role of differentially methylated genes in the pathogenesis of the disease.

METHODS

Whole-genome bisulfite sequencing was performed on three SS patients and four HC. The biological function of genes associated with differentially methylated regions (DMRs) was investigated using Gene Ontology functional analysis and Kyoto Encyclopedia of Genes and Genomes pathway analysis, using network-based key driver analysis (KDA) to find KDA genes. In clinical samples of SS patients and controls, the expression levels of KDA genes were validated by quantitative real-time polymerase chain reaction and immunohistochemical analysis. Moreover, the diagnostic value of KDA genes for SS was confirmed using receiver operating characteristic curves.

RESULTS

We identified 322 DMRs, annotated as 162 associated genes. Six genes were selected via the number of networks of KDA genes. Differential expression of genes such as human leukocyte antigen (HLA) class I, ADAR, and OAS2 was observed in patients' peripheral blood mononuclear cells and the minor salivary glands, which can be used as potential diagnostic biomarkers for SS.

CONCLUSION

Clinical sample validation suggested that HLA class I, ADAR, and OAS2 might play a role in the development of SS. Our study shows epigenetic regulatory mechanisms and potential disease markers associated with SS, which in turn will enable us to identify new therapeutic targets.

Effects of In Utero PFOS Exposure on Epigenetics and Metabolism in Mouse Fetal Livers

Prenatal exposure to perfluorooctanesulfonate (PFOS) increases fetus' metabolic risk; however, the investigation of the underlying mechanism is limited. In this study, pregnant mice in the gestational days (GD, 4.5-17.5) were exposed to PFOS (0.3 and 3 μg/g of body weight). At GD 17.5, PFOS perturbed maternal lipid metabolism and upregulated metabolism-regulating hepatokines (Angptl4, Angptl8, and Selenop). Mass-spectrometry imaging and whole-genome bisulfite sequencing revealed, respectively, selective PFOS localization and deregulation of gene methylation in fetal livers, involved in inflammation, glucose, and fatty acid metabolism. PCR and Western blot analysis of lipid-laden fetal livers showed activation of AMPK signaling, accompanied by significant increases in the expression of glucose transporters (Glut2/4), hexose-phosphate sensors (Retsat and ChREBP), and the key glycolytic enzyme, pyruvate kinase (Pk) for glucose catabolism. Additionally, PFOS modulated the expression levels of PPARα and PPARγ downstream target genes, which simultaneously stimulated fatty acid oxidation (Cyp4a14, Acot, and Acox) and lipogenesis (Srebp1c, Acaca, and Fasn). Using human normal hepatocyte (MIHA) cells, the underlying mechanism of PFOS-elicited nuclear translocation of ChREBP, associated with a fatty acid synthesizing pathway, was revealed. Our finding implies that in utero PFOS exposure altered the epigenetic landscape associated with dysregulation of fetal liver metabolism, predisposing postnatal susceptibility to metabolic challenges.

Differentially methylated genomic regions of lettuce seeds relate to divergence across morphologically distinct horticultural types

Heritable cytosine methylation plays a role in shaping plant phenotypes; however, no information is available about DNA methylation in cultivated lettuce (Lactuca sativa), one of the most important leafy vegetables. Whole-genome bisulfite sequencing (WGBS) performed on seeds of 95 accessions from eight morphologically distinct horticultural types (Batavia, butterhead, iceberg, Latin, leaf, oilseed, romaine and stem) revealed a high level of methylation in lettuce genome with an average methylation of 90.6 % in the CG context, 72.9 % in the CHG context and 7.5 % in the CHH context. Although WGBS did not show substantial differences in overall methylation levels across eight horticultural types, 350 differentially methylated regions (DMR) were identified. Majority of the 41 pivotal DMR overlapped with genomic features predicted or confirmed to be involved in plant growth and development. These results provide the first insight into lettuce DNA methylation and indicate a potential role for heritable variation in cytosine methylation in lettuce morphology. The results reveal that differences in methylation profiles of morphologically distinct horticultural types are already detectable in seeds. Identified DMR can be a focus of the future functional studies.

Promoter Methylation-Regulated Differentially Expressed Genes in Breast Cancer

Background

Breast cancer is one of the most common malignancies among women. Recent studies revealed that differentially methylated regions (DMRs) are implicated in regulating gene expression. The goal of this research was to determine which genes and pathways are dysregulated in breast cancer when their promoters are methylated in an abnormal way, leading to differential expression. Whole-genome bisulfite sequencing was applied to analyze DMRs for eight peripheral blood samples collected from five Saudi females diagnosed with stages I and II of breast cancer aligned with three normal females. Three of those patients and three normal samples were used to determine differentially expressed genes (DEG) using Illumina platform NovaSeq PE150.

Results

Based on ontology (GO) and KEGG pathways, the analysis indicated that DMGs and DEG are closely related to associated processes, such as ubiquitin-protein transferase activity, ubiquitin-mediated proteolysis, and oxidative phosphorylation. The findings indicated a potentially significant association between global hypomethylation and breast cancer in Saudi patients. Our results revealed 81 differentially promoter-methylated and expressed genes. The most significant differentially methylated and expressed genes found in gene ontology (GO) are pumilio RNA binding family member 1 (PUM1) and zinc finger AN1-type containing 2B (ZFAND2B) also known as (AIRAPL).

Conclusion

The essential outcomes of this study suggested that aberrant hypermethylation at crucial genes that have significant parts in the molecular pathways of breast cancer could be used as a potential prognostic biomarker for breast cancer.

Accelerating the alignment processing speed of the comprehensive end-to-end whole-genome bisulfite sequencing pipeline, wg-blimp

Analyzing whole-genome bisulfite and related sequencing datasets is a time-intensive process due to the complexity and size of the input raw sequencing files and lengthy read alignment step requiring correction for conversion of all unmethylated Cs to Ts genome-wide. The objective of this study was to modify the read alignment algorithm associated with the whole-genome bisulfite sequencing methylation analysis pipeline (wg-blimp) to shorten the time required to complete this phase while retaining overall read alignment accuracy. Here, we report an update to the recently published pipeline wg-blimp achieved by replacing the use of the bwa-meth aligner with the faster gemBS aligner. This improvement to the wg-blimp pipeline has led to a more than ×7 acceleration in the processing speed of samples when scaled to larger publicly available FASTQ datasets containing 80-160 million reads while maintaining nearly identical accuracy of properly mapped reads when compared with data from the previous pipeline. The modifications to the wg-blimp pipeline reported here merge the speed and accuracy of the gemBS aligner with the comprehensive analysis and data visualization assets of the wg-blimp pipeline to provide a significantly accelerated workflow that can produce high-quality data much more rapidly without compromising read accuracy at the expense of increasing RAM requirements up to 48 GB.

Whole-genome DNA methylation profiling reveals epigenetic signatures in developing muscle in Tan and Hu sheep and their offspring

Introduction

The Tan sheep is a popular local breed in China because of its tenderness and flavor. The Hu sheep breed is also famous for its high litter size, and its muscle growth rate is faster than that of Tan sheep. However, the epigenetic mechanism behind these muscle-related phenotypes is unknown.

Methods

In this study, the longissimus dorsi tissue from 18 6 month-old Tan sheep, Hu sheep, and Tan-Hu F2 generation (6 sheep per population) were collected. After genomic DNA extraction, whole-genome bisulfite sequencing (WGBS) and bioinformatics analysis were performed to construct genome-wide DNA methylome maps for the Tan sheep, Hu sheep and their Tan-Hu F2 generation.

Results

Distinct genome-wide DNA methylation patterns were observed between Tan sheep and Hu sheep. Moreover, DNA methylated regions were significantly increased in the skeletal muscle from Tan sheep vs. the F2 generation compared to the Hu sheep vs. F2 generation and the Tan sheep vs. Hu sheep. Compared with Hu sheep, the methylation levels of actin alpha 1 (ACTA1), myosin heavy chain 11 (MYH11), Wiskott-Aldrich syndrome protein (WAS), vav guanine nucleotide exchange factor 1 (VAV1), fibronectin 1 (FN1) and Rho-associated protein kinase 2 (ROCK2) genes were markedly distinct in the Tan sheep. Furthermore, Gene Ontology analysis indicated that these genes were involved in myotube differentiation, myotube cell development, smooth muscle cell differentiation and striated muscle cell differentiation.

Conclusion

The findings from this study, in addition to data from previous research, demonstrated that the ACTA1, MYH11, WAS, VAV1, FN1, and ROCK2 genes may exert regulatory effects on muscle development.

Chromosome-scale de novo genome assembly and annotation of three representative Casuarina species: C. equisetifolia, C. glauca, and C. cunninghamiana

Australian pine (Casuarina spp.) is extensively planted in tropical and subtropical regions for wood production, shelterbelts, environmental protection, and ecological restoration due to their superior biological characteristics, such as rapid growth, wind and salt tolerance, and nitrogen fixation. To analyze the genomic diversity of Casuarina, we sequenced the genomes and constructed de novo genome assemblies of the three most widely planted Casuarina species: C. equisetifolia, C. glauca, and C. cunninghamiana. We generated chromosome-scale genome sequences using both Pacific Biosciences (PacBio) Sequel sequencing and chromosome conformation capture technology (Hi-C). The total genome sizes for C. equisetifolia, C. glauca, and C. cunninghamiana are 268 942 579 bp, 296 631 783 bp, and 293 483 606 bp, respectively, of which 25.91, 27.15, and 27.74% were annotated as repetitive sequences. We annotated 23 162, 24 673, and 24 674 protein-coding genes in C. equisetifolia, C. glauca, and C. cunninghamiana, respectively. We then collected branchlets from male and female individuals for whole-genome bisulfite sequencing (BS-seq) to explore the epigenetic regulation of sex determination in these three species. Transcriptome sequencing (RNA-seq) revealed differential expression of phytohormone-related genes between male and female plants. In summary, we generated three chromosome-level genome assemblies and comprehensive DNA methylation and transcriptome datasets from both male and female material for three Casuarina species, providing a basis for the comprehensive investigation of genomic diversity and functional gene discovery of Casuarina in the future.

Epigenomic signature of major congenital heart defects in newborns with Down syndrome

Background

Congenital heart defects (CHDs) affect approximately half of individuals with Down syndrome (DS) but the molecular reasons for incomplete penetrance are unknown. Previous studies have largely focused on identifying genetic risk factors associated with CHDs in individuals with DS, but comprehensive studies of the contribution of epigenetic marks are lacking. We aimed to identify and characterize DNA methylation differences from newborn dried blood spots (NDBS) of DS individuals with major CHDs compared to DS individuals without CHDs.

Methods

We used the Illumina EPIC array and whole-genome bisulfite sequencing (WGBS) to quantitate DNA methylation for 86 NDBS samples from the California Biobank Program: 1) 45 DS-CHD (27 female, 18 male) and 2) 41 DS non-CHD (27 female, 14 male). We analyzed global CpG methylation and identified differentially methylated regions (DMRs) in DS-CHD vs DS non-CHD comparisons (both sex-combined and sex-stratified) corrected for sex, age of blood collection, and cell type proportions. CHD DMRs were analyzed for enrichment in CpG and genic contexts, chromatin states, and histone modifications by genomic coordinates and for gene ontology enrichment by gene mapping. DMRs were also tested in a replication dataset and compared to methylation levels in DS vs typical development (TD) WGBS NDBS samples.

Results

We found global CpG hypomethylation in DS-CHD males compared to DS non-CHD males, which was attributable to elevated levels of nucleated red blood cells and not seen in females. At a regional level, we identified 58, 341, and 3,938 CHD-associated DMRs in the Sex Combined, Females Only, and Males Only groups, respectively, and used machine learning algorithms to select 19 Males Only loci that could distinguish CHD from non-CHD. DMRs in all comparisons were enriched for gene exons, CpG islands, and bivalent chromatin and mapped to genes enriched for terms related to cardiac and immune functions. Lastly, a greater percentage of CHD-associated DMRs than background regions were differentially methylated in DS vs TD samples.

Conclusions

A sex-specific signature of DNA methylation was detected in NDBS of DS-CHD compared to DS non-CHD individuals. This supports the hypothesis that epigenetics can reflect the variability of phenotypes in DS, particularly CHDs.

Beadchip technology to detect DNA methylation in mouse faithfully recapitulates whole-genome bisulfite sequencing

Aim: To facilitate wide-scale implementation of Illumina Mouse Methylation BeadChip (MMB) technology, array-based measurement of cytosine methylation was compared with the gold-standard assessment of DNA methylation by whole-genome bisulfite sequencing (WGBS). Methods: DNA methylation across two mouse strains (C57B6 and C3H) and both sexes was assessed using the MMB and compared with previously existing deep-coverage WGBS of mice of the same strain and sex. Results & conclusion: The findings demonstrated that 93.3-99.2% of sites had similar measurements of methylation across technologies and that differentially methylated cytosines and regions identified by each technology overlap and enrich for similar biological functions, suggesting that the MMB faithfully recapitulates the findings of WGBS.

The global dynamic of DNA methylation in response to heat stress revealed epigenetic mechanism of heat acclimation in Saccharina japonica

Saccharina japonica is an ecologically and economically important kelp in cold-temperate regions. When it is cultivated on a large scale in the temperate and even subtropical zones, heat stress is a frequent abiotic stress. This study is the first attempt to reveal the regulatory mechanism of the response to heat stress from the perspective of DNA methylation in S. japonica. We firstly obtained the characteristics of variation in the methylome under heat stress, and observed that heat stress caused a slight increase in the overall methylation level and methylation rate, especially in the non-coding regions of the genome. Secondly, we noted that methylation was probably one of factors affecting the expression of genes, and that methylation within the gene body was positively correlated with the gene expression (rho = 0.0784). Moreover, it was found that among the differentially expressed genes regulated by methylation, many genes were related to heat stress response, such as HSP gene family, genes of antioxidant enzymes, genes related to proteasome-ubiquitination pathway, and plant cell signaling pathways. This study demonstrated that DNA methylation is involved in regulating the response to heat stress, laying a foundation for studying the acclimation and adaptation of S. japonica to heat stress from an epigenetic perspective.

Whole-genome resequencing and bisulfite sequencing provide new insights into the feeding habit domestication in mandarin fish (Siniperca chuatsi)

Mandarin fish (Siniperca chuatsi) is one of the most economically important fish in China. However, it has the peculiar feeding habit that it feeds solely on live prey fish since first-feeding, while refuses dead prey fish or artificial diets. After the specific training procedure, partial individuals could accept dead prey fish and artificial diets. The genetic basis of individual difference in artificial diet feeding habit is still unknown. In the present study, the resequencing was performed between 10 individuals which could be domesticated to accept artificial diets and 10 individuals which could not. Through the selective sweep analysis based on heterozygosity (Hp) and population differentiation coefficient (Fst), 57 candidate windows were identified as the putative selected regions for feeding habit domestication of mandarin fish, involved in 149 genes. These genes were related to memory, vision and olfaction function, which could be potential targets of molecular marker assistant breeding of artificial diet feeding trait. Beside of the DNA sequence, we also explored the potential role of DNA methylation in feeding habit domestication in mandarin fish. Whole-genome bisulfite sequencing was performed between the individuals which could be domesticated to accept artificial diets and those could not. 5,976 differentially methylated regions were identified, referring to 3,522 genes, such as the genes involved in cAMP signaling pathway. The DNA methylation changes of these genes might contribute to the adaption of artificial diets in mandarin fish. In conclusion, the putative selected regions and the differentially methylated regions were identified in the whole genome, providing new insights into the feeding habit domestication from live prey fish to artificial diets in mandarin fish. And the involved genes were identified as the candidate genes for molecular breeding of artificial diet utilization in mandarin fish.

Evaluation of cell-free DNA approaches for multi-cancer early detection

In the Circulating Cell-free Genome Atlas (NCT02889978) substudy 1, we evaluate several approaches for a circulating cell-free DNA (cfDNA)-based multi-cancer early detection (MCED) test by defining clinical limit of detection (LOD) based on circulating tumor allele fraction (cTAF), enabling performance comparisons. Among 10 machine-learning classifiers trained on the same samples and independently validated, when evaluated at 98% specificity, those using whole-genome (WG) methylation, single nucleotide variants with paired white blood cell background removal, and combined scores from classifiers evaluated in this study show the highest cancer signal detection sensitivities. Compared with clinical stage and tumor type, cTAF is a more significant predictor of classifier performance and may more closely reflect tumor biology. Clinical LODs mirror relative sensitivities for all approaches. The WG methylation feature best predicts cancer signal origin. WG methylation is the most promising technology for MCED and informs development of a targeted methylation MCED test.

Promoter-Adjacent DNA Hypermethylation Can Downmodulate Gene Expression: TBX15 in the Muscle Lineage

TBX15, which encodes a differentiation-related transcription factor, displays promoter-adjacent DNA hypermethylation in myoblasts and skeletal muscle (psoas) that is absent from non-expressing cells in other lineages. By whole-genome bisulfite sequencing (WGBS) and enzymatic methyl-seq (EM-seq), these hypermethylated regions were found to border both sides of a constitutively unmethylated promoter. To understand the functionality of this DNA hypermethylation, we cloned the differentially methylated sequences (DMRs) in CpG-free reporter vectors and tested them for promoter or enhancer activity upon transient transfection. These cloned regions exhibited strong promoter activity and, when placed upstream of a weak promoter, strong enhancer activity specifically in myoblast host cells. In vitro CpG methylation targeted to the DMR sequences in the plasmids resulted in 86−100% loss of promoter or enhancer activity, depending on the insert sequence. These results as well as chromatin epigenetic and transcription profiles for this gene in various cell types support the hypothesis that DNA hypermethylation immediately upstream and downstream of the unmethylated promoter region suppresses enhancer/extended promoter activity, thereby downmodulating, but not silencing, expression in myoblasts and certain kinds of skeletal muscle. This promoter-border hypermethylation was not found in cell types with a silent TBX15 gene, and these cells, instead, exhibit repressive chromatin in and around the promoter. TBX18, TBX2, TBX3 and TBX1 display TBX15-like hypermethylated DMRs at their promoter borders and preferential expression in myoblasts. Therefore, promoter-adjacent DNA hypermethylation for downmodulating transcription to prevent overexpression may be used more frequently for transcription regulation than currently appreciated.

DNA Methylation Variation Is a Possible Mechanism in the Response of Haemaphysalis longicornis to Low-Temperature Stress

Ticks are notorious ectoparasites and transmit the greatest variety of pathogens than any other arthropods. Cold tolerance is a key determinant of tick abundance and distribution. While studies have shown that DNA methylation is one of the important epigenetic regulations found across many species and plays a significant role in their response to low-temperature stress, its role in the response of ticks to low-temperature stress remains unexplored. Herein, we explored the DNA methylation profile of the tick, Haemaphysalis longicornis, exposed to low-temperature stress (4 °C) using whole-genome bisulfite sequencing (WGBS). We found that approximately 0.95% and 0.94% of the genomic C sites were methylated in the control and low-temperature groups, respectively. Moreover, the methylation level under the CG context was about 3.86% and 3.85% in the control and low-temperature groups, respectively. In addition, a total of 6087 differentially methylated regions (DMRs) were identified between the low-temperature and control groups, including 3288 hypermethylated and 2799 hypomethylated DMRs. Further, Gene ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis of differentially methylated genes revealed that most of the DMGs were significantly enriched in binding and RNA transport pathways. Taken together, this research confirmed, for the first time, the whole genome DNA methylation profile of H. longicornis and provided new insights into the DNA methylation changes relating to low-temperature stress in H. longicornis, as well as provided a foundation for future studies on the epigenetic mechanism underlying the responses of ticks to abiotic stress.

Integrated analyses of the methylome and transcriptome to unravel sex differences in the perirenal fat from suckling lambs

In sheep, differences were observed regarding fat accumulation and fatty acid (FA) composition between males and females, which may impact the quality and organoleptic characteristics of the meat. The integration of different omics technologies is a relevant approach for investigating biological and genetic mechanisms associated with complex traits. Here, the perirenal tissue of six male and six female Assaf suckling lambs was evaluated using RNA sequencing and whole-genome bisulfite sequencing (WGBS). A multiomic discriminant analysis using multiblock (s)PLS-DA allowed the identification of 314 genes and 627 differentially methylated regions (within these genes), which perfectly discriminate between males and females. These candidate genes overlapped with previously reported QTLs for carcass fat volume and percentage of different FAs in milk and meat from sheep. Additionally, differentially coexpressed (DcoExp) modules of genes between males (nine) and females (three) were identified that harbour 22 of these selected genes. Interestingly, these DcoExp were significantly correlated with fat percentage in different deposits (renal, pelvic, subcutaneous and intramuscular) and were associated with relevant biological processes for adipogenesis, adipocyte differentiation, fat volume and FA composition. Consequently, these genes may potentially impact adiposity and meat quality traits in a sex-specific manner, such as juiciness, tenderness and flavour.

Whole-genome bisulfite sequencing analysis of circulating tumour DNA for the detection and molecular classification of cancer

Background

Cancer cell–specific variation and circulating tumour DNA (ctDNA) methylation are promising biomarkers for non‐invasive cancer detection and molecular classification. Nevertheless, the applications of ctDNA to the early detection and screening of cancer remain highly challenging due to the scarcity of cancer cell–specific ctDNA, the low signal‐to‐noise ratio of DNA variation, and the lack of non‐locus‐specific DNA methylation technologies.

Methods

We enrolled three cohorts of breast cancer (BC) patients from two hospitals in China (BC: n = 123; healthy controls: n = 40). We developed a ctDNA whole‐genome bisulfite sequencing technology employing robust trace ctDNA capture from up to 200 μL plasma, mini‐input (1 ng) library preparation, unbiased genome‐wide coverage and comprehensive computational methods.

Results

A diagnostic signature comprising 15 ctDNA methylation markers exhibited high accuracy in the early (area under the curve [AUC] of 0.967) and advanced (AUC of 0.971) BC stages in multicentre patient cohorts. Furthermore, we revealed a ctDNA methylation signature that discriminates estrogen receptor status (Training set: AUC of 0.984 and Test set: AUC of 0.780). Different cancer types, including hepatocellular carcinoma and lung cancer, could also be well distinguished.

Conclusions

Our study provides a toolset to generate unbiased whole‐genome ctDNA methylomes with a minimal amount of plasma to develop highly specific and sensitive biomarkers for the early diagnosis and molecular subtyping of cancer.

Plasticity across levels: relating epigenomic, transcriptomic, and phenotypic responses to osmotic stress in a halotolerant microalga

Phenotypic plasticity, the ability of a given genotype to produce alternative phenotypes in response to its environment of development, is an important mechanism for coping with variable environments. While the mechanisms underlying phenotypic plasticity are diverse, their relative contributions need to be investigated quantitatively to better understand the evolvability of plasticity across biological levels. This requires relating plastic responses of the epigenome, transcriptome, and organismal phenotype, and investigating how they vary with the genotype. Here we carried out this approach for responses to osmotic stress in Dunaliella salina, a green microalga that is a model organism for salinity tolerance. We compared two strains that show markedly different demographic responses to osmotic stress, and showed that these phenotypic responses involve strain‐ and environment‐specific variation in gene expression levels, but a relative low—albeit significant—effect of strain × environment interaction. We also found an important genotype effect on the genome‐wide methylation pattern, but little contribution from environmental conditions to the latter. However, we did detect a significant marginal effect of epigenetic variation on gene expression, beyond the influence of genetic differences on epigenetic state, and we showed that hypomethylated regions are correlated with higher gene expression. Our results indicate that epigenetic mechanisms are either not involved in the rapid plastic response to environmental change in this species, or involve only few changes in trans that are sufficient to trigger concerted changes in the expression of many genes, and phenotypic responses by multiple traits.

Novel deoxyribonucleic acid methylation perturbations in workers exposed to vinyl chloride

To explore the epigenetic mechanism of deoxyribonucleic acid (DNA) damage induced by vinyl chloride (VC), we studied the micronuclei of peripheral blood lymphocytes in 193 subjects (92 in a VC exposure group employed in a chlorine-alkali plant; 101 in a control group employed in a power plant) and selected three pairs from the subjects (exposed and control) for whole-genome bisulfite sequencing (WGBS). The results showed that the rate of micronucleus formation in the VC exposure group was higher than that of control group (6.05 ± 3.28‰ vs. 2.01 ± 1.79‰). A total of 9534 differentially methylated regions (DMRs) were identified by WGBS, of which 4816 were hypomethylated and 4718 were hypermethylated. The Kyoto encyclopedia of genes and genomes (KEGG) pathway and gene ontology (GO) analyses showed the top three KEGG pathways were cancer , neuroactive ligand-receptor interaction, and axon guidance, and the top three GO-BP pathways enriched were multicellular organismal process, developmental process, and anatomical structure development. In the most enriched DMR pathway (pathways in cancer), we found that BCL2, TJP2, TAOK1, PFKFB3, LIPI, and LIPH were hypermethylated, and the methylation levels of BNIP1 and GRPEL2 were decreased. The methylation of differentially methylated genes (DMGs) mentioned above were verified by methylation-specific PCR (MSP) and agarose gel electrophoresis (AGE) in 50 pairs of subjects, where the coincidence rate was 60-100%. In conclusion, the epigenetic perturbations of specific DMGs (BCL2, TJP2, TAOK1, PFKFB3, LIPI, LIPH, BNIP1, and GRPEL2) may be associated with DNA damage from vinyl chloride exposure.

Comethyl: a network-based methylome approach to investigate the multivariate nature of health and disease

Health outcomes are frequently shaped by difficult to dissect inter-relationships between biological, behavioral, social and environmental factors. DNA methylation patterns reflect such multivariate intersections, providing a rich source of novel biomarkers and insight into disease etiologies. Recent advances in whole-genome bisulfite sequencing enable investigation of DNA methylation over all genomic CpGs, but existing bioinformatic approaches lack accessible system-level tools. Here, we develop the R package Comethyl, for weighted gene correlation network analysis of user-defined genomic regions that generates modules of comethylated regions, which are then tested for correlations with multivariate sample traits. First, regions are defined by CpG genomic location or regulatory annotation and filtered based on CpG count, sequencing depth and variability. Next, correlation networks are used to find modules of interconnected nodes using methylation values within the selected regions. Each module containing multiple comethylated regions is reduced in complexity to a single eigennode value, which is then tested for correlations with experimental metadata. Comethyl has the ability to cover the noncoding regulatory regions of the genome with high relevance to interpretation of genome-wide association studies and integration with other types of epigenomic data. We demonstrate the utility of Comethyl on a dataset of male cord blood samples from newborns later diagnosed with autism spectrum disorder (ASD) versus typical development. Comethyl successfully identified an ASD-associated module containing regions mapped to genes enriched for brain glial functions. Comethyl is expected to be useful in uncovering the multivariate nature of health disparities for a variety of common disorders. Comethyl is available at github.com/cemordaunt/comethyl with complete documentation and example analyses.

Placenta and fetal brain share a neurodevelopmental disorder DNA methylation profile in a mouse model of prenatal PCB exposure

Polychlorinated biphenyls (PCBs) are developmental neurotoxicants implicated as environmental risk factors for neurodevelopmental disorders (NDDs). Here, we report the effects of prenatal exposure to a human-relevant mixture of PCBs on the DNA methylation profiles of mouse placenta and fetal brain. Thousands of differentially methylated regions (DMRs) distinguish placenta and fetal brain from PCB-exposed mice from sex-matched vehicle controls. In both placenta and fetal brain, PCB-associated DMRs are enriched for functions related to neurodevelopment and cellular signaling and enriched within regions of bivalent chromatin. The placenta and brain PCB DMRs overlap significantly and map to a shared subset of genes enriched for Wnt signaling, Slit/Robo signaling, and genes differentially expressed in NDD models. The consensus PCB DMRs also significantly overlap with DMRs from human NDD brain and placenta. These results demonstrate that PCB-exposed placenta contains a subset of DMRs that overlap fetal brain DMRs relevant to an NDD.

Genome-Wide DNA Methylation Profile in Jejunum Reveals the Potential Genes Associated With Paratuberculosis in Dairy Cattle

Paratuberculosis in cattle causes substantial economic losses to the dairy industry. Exploring functional genes and corresponding regulatory pathways related to resistance or susceptibility to paratuberculosis is essential to the breeding of disease resistance in cattle. Co-analysis of genome-wide DNA methylation and transcriptome profiles is a critically important approach to understand potential regulatory mechanism underlying the development of diseases. In this study, we characterized the profiles of DNA methylation of jejunum from nine Holstein cows in clinical, subclinical, and healthy groups using whole-genome bisulfite sequencing (WGBS). The average methylation level in functional regions was 29.95% in the promoter, 29.65% in the 5’ untranslated region (UTR), 68.24% in exons, 71.55% in introns, and 72.81% in the 3’ UTR. A total of 3,911, 4,336, and 4,094 differentially methylated genes (DMGs) were detected in clinical vs. subclinical, clinical vs. healthy, and subclinical vs. healthy comparative group, respectively. Gene ontology (GO) and analysis based on the Kyoto Encyclopedia of Genes and Genomes (KEGG) showed that these DMGs were significantly enriched in specific biological processes related to immune response, such as Th1 and Th2 cell differentiation, wnt, TNF, MAPK, ECM-receptor interaction, cellular senescence, calcium, and chemokine signaling pathways (q value <0.05). The integration of information about DMGs, differentially expressed genes (DEGs), and biological functions suggested nine genes CALCRL, TNC, GATA4, CD44, TGM3, CXCL9, CXCL10, PPARG, and NFATC1 as promising candidates related to resistance/susceptibility to Mycobacterium avium subspecies paratuberculosis (MAP). This study reports on the high-resolution DNA methylation landscapes of the jejunum methylome across three conditions (clinical, subclinical, and healthy) in dairy cows. Our investigations integrated different sources of information about DMGs, DEGs, and pathways, enabling us to find nine functional genes that might have potential application in resisting paratuberculosis in dairy cattle.

Genomic Imprinting at the Porcine DIRAS3 Locus

Simple Summary

DNA methylation associated with one of the two alleles from parents is an important mechanism that causes a silencing of that allele, leading to expression of another allele only. There has been a lack of detailed studies on DNA methylation and expression patterns that are related to the DIRAS3 gene in pigs. The objective of this study was to provide a comprehensive overview of DNA methylation and expression associated with the DIRAS3 gene in pigs by generating an embryonic pig model and analyzing next-generation sequencing using pig embryos and adult pigs. Our results clearly showed the presence of DNA methylation near the DIRAS3 gene in pigs and high expression of DIRAS3 in the hypothalamus from adult pigs and expression of only one allele in all the tested tissues including the hypothalamus. In summary, our findings suggested DNA methylation might be related to those unique gene expression patterns during the development of pigs.

Abstract

The epigenetic mechanisms underlying genomic imprinting include DNA methylation and monoallelic expression of genes in close proximity. Although genes imprinted in humans and mice have been widely characterized, there is a lack of detailed and comprehensive studies in livestock species including pigs. The purpose of this study was to investigate a detailed methylation status and parent-of-origin-specific gene expression within the genomic region containing an underexamined porcine DIRAS3 locus. Through whole-genome bisulfite sequencing (WGBS) and RNA sequencing (RNA-seq) of porcine parthenogenetic embryos and analyses of public RNA-seq data from adult pigs, DNA methylation and monoallelic expression pattern were investigated. As a result, maternal hypermethylation at the DIRAS3 locus and hypothalamus-specific and monoallelic expression of the DIRAS3 gene were found in pigs. In conclusion, the findings from this study suggest that the presence of maternal hypermethylation, or imprints, might be maintained and related to monoallelic expression of DIRAS3 during pig development.

Evaluating the Consistency of Gene Methylation in Liver Cancer Using Bisulfite Sequencing Data

Bisulfite sequencing is considered as the gold standard approach for measuring DNA methylation, which acts as a pivotal part in regulating a variety of biological processes without changes in DNA sequences. In this study, we introduced the most prevalent methods for processing bisulfite sequencing data and evaluated the consistency of the data acquired from different measurements in liver cancer. Firstly, we introduced three commonly used bisulfite sequencing assays, i.e., reduced-representation bisulfite sequencing (RRBS), whole-genome bisulfite sequencing (WGBS), and targeted bisulfite sequencing (targeted BS). Next, we discussed the principles and compared different methods for alignment, quality assessment, methylation level scoring, and differentially methylated region identification. After that, we screened differential methylated genes in liver cancer through the three bisulfite sequencing assays and evaluated the consistency of their results. Ultimately, we compared bisulfite sequencing to 450 k beadchip and assessed the statistical similarity and functional association of differentially methylated genes (DMGs) among the four assays. Our results demonstrated that the DMGs measured by WGBS, RRBS, targeted BS and 450 k beadchip are consistently hypo-methylated in liver cancer with high functional similarity.

Epigenetic inactivation of ERF reactivates γ-globin expression in β-thalassemia

The fetal-to-adult hemoglobin switch is regulated in a developmental stage-specific manner and reactivation of fetal hemoglobin (HbF) has therapeutic implications for treatment of β-thalassemia and sickle cell anemia, two major global health problems. Although significant progress has been made in our understanding of the molecular mechanism of the fetal-to-adult hemoglobin switch, the mechanism of epigenetic regulation of HbF silencing remains to be fully defined. Here, we performed whole-genome bisulfite sequencing and RNA sequencing analysis of the bone marrow-derived GYPA⁺ erythroid cells from β-thalassemia-affected individuals with widely varying levels of HbF groups (HbF ≥ 95th percentile or HbF ≤ 5th percentile) to screen epigenetic modulators of HbF and phenotypic diversity of β-thalassemia. We identified an ETS2 repressor factor encoded by ERF, whose promoter hypermethylation and mRNA downregulation are associated with high HbF levels in β-thalassemia. We further observed that hypermethylation of the ERF promoter mediated by enrichment of DNMT3A leads to demethylation of γ-globin genes and attenuation of binding of ERF on the HBG promoter and eventually re-activation of HbF in β-thalassemia. We demonstrated that ERF depletion markedly increased HbF production in human CD34⁺ erythroid progenitor cells, HUDEP-2 cell lines, and transplanted NCG-Kit-V831M mice. ERF represses γ-globin expression by directly binding to two consensus motifs regulating γ-globin gene expression. Importantly, ERF depletion did not affect maturation of erythroid cells. Identification of alterations in DNA methylation of ERF as a modulator of HbF synthesis opens up therapeutic targets for β-hemoglobinopathies.

Koala methylomes reveal divergent and conserved DNA methylation signatures of X chromosome regulation

X chromosome inactivation (XCI) mediated by differential DNA methylation between sexes is an iconic example of epigenetic regulation. Although XCI is shared between eutherians and marsupials, the role of DNA methylation in marsupial XCI remains contested. Here, we examine genome-wide signatures of DNA methylation across fives tissues from a male and female koala (Phascolarctos cinereus), and present the first whole-genome, multi-tissue marsupial ‘methylome atlas’. Using these novel data, we elucidate divergent versus common features of representative marsupial and eutherian DNA methylation. First, tissue-specific differential DNA methylation in koalas primarily occurs in gene bodies. Second, females show significant global reduction (hypomethylation) of X chromosome DNA methylation compared to males. We show that this pattern is also observed in eutherians. Third, on average, promoter DNA methylation shows little difference between male and female koala X chromosomes, a pattern distinct from that of eutherians. Fourth, the sex-specific DNA methylation landscape upstream of Rsx, the primary lncRNA associated with marsupial XCI, is consistent with the epigenetic regulation of female-specific (and presumably inactive X chromosome-specific) expression. Finally, we use the prominent female X chromosome hypomethylation and classify 98 previously unplaced scaffolds as X-linked, contributing an additional 14.6 Mb (21.5%) to genomic data annotated as the koala X chromosome. Our work demonstrates evolutionarily divergent pathways leading to functionally conserved patterns of XCI in two deep branches of mammals.

Whole-genome bisulfite sequencing of abdominal adipose reveals DNA methylation pattern variations in broiler lines divergently selected for fatness

The methylation status of pivotal genes involved in fat deposition in chickens has been extensively studied. However, the whole-genome DNA methylation profiles of broiler abdominal adipose tissue remain poorly understood. Using whole-genome bisulfite sequencing, we generated DNA methylation profiles of chicken abdominal adipose tissue from Northeast Agricultural University broiler lines divergently selected for abdominal fat content. We aimed to explore whether DNA methylation was associated with abdominal fat deposition in broilers. The whole-genome DNA methylation profiles of fat- and lean-line broilers abdominal adipose tissue were constructed. The DNA methylation levels of functional genomic regions in the fat broiler were higher than those in the lean broiler, especially in the 3' untranslated regions (UTRs) and exons in the non-CG contexts. Additionally, we identified 29,631 differentially methylated regions and, subsequently, annotated 6,484 and 2,016 differentially methylated genes (DMGs) in the gene body and promoter regions between the two lines, respectively. Functional annotation showed that the DMGs in promoter regions were significantly enriched mainly in the triglyceride catabolic process, lipid metabolism-related pathways, and extracellular matrix signal pathways. When the DMG in promoter regions and differentially expressed genes were integrated, we identified 30 genes with DNA methylation levels that negatively correlated with their messenger RNA (mRNA) expression, of which CMSS1 reached significant levels (false discovery rate < 0.05). These 30 genes were mainly involved in fatty acid metabolism, peroxisome-proliferator-activated receptor signaling, Wnt signaling pathways, transmembrane transport, RNA degradation, and glycosaminoglycan degradation. Comparing the DNA methylation profiles between fat- and lean-line broilers demonstrated that DNA methylation is involved in regulating broiler abdominal fat deposition. Our study offers a basis for further exploring the underlying mechanisms of abdominal adipose deposition in broilers.

The Landscape of Genomic Imprinting at the Porcine SGCE/PEG10 Locus from Methylome and Transcriptome of Parthenogenetic Embryos

In mammals, imprinted genes often exist in the form of clusters in specific chromosome regions. However, in pigs, genomic imprinting of a relatively few genes and clusters has been identified, and genes within or adjacent to putative imprinted clusters need to be investigated including those at the SGCE/PEG10 locus. The objective of this study was to, using porcine parthenogenetic embryos, investigate imprinting status of genes within the genomic region spans between the COL1A2 and ASB4 genes in chromosome 9. Whole-genome bisulfite sequencing (WGBS) and RNA sequencing (RNA-seq) were conducted with normal and parthenogenetic embryos, and methylome and transcriptome were analyzed. As a result, differentially methylated regions (DMRs) between the embryos were identified, and parental allele-specific expressions of the SGCE and PEG10 genes were verified. The pig imprinted interval was limited between SGCE and PEG10, since both the COL1A2 and CASD1 genes at the centromere-proximal region and the genes between PPP1R9A and ASB4 toward the telomere were non-imprinted and biallelically expressed. Consequently, our combining analyses of methylome, transcriptome, and informative polymorphisms revealed the boundary of imprinting cluster at the SGCE/PEG10 locus in pig chromosome 9 and consolidated the landscape of genomic imprinting in pigs.

The Absence of C-5 DNA Methylation in Leishmania donovani Allows DNA Enrichment from Complex Samples

Cytosine C5 methylation is an important epigenetic control mechanism in a wide array of eukaryotic organisms and generally carried out by proteins of the C-5 DNA methyltransferase family (DNMTs). In several protozoans, the status of this mechanism remains elusive, such as in Leishmania, the causative agent of the disease leishmaniasis in humans and a wide array of vertebrate animals. In this work, we showed that the Leishmania donovani genome contains a C-5 DNA methyltransferase (DNMT) from the DNMT6 subfamily, whose function is still unclear, and verified its expression at the RNA level. We created viable overexpressor and knock-out lines of this enzyme and characterized their genome-wide methylation patterns using whole-genome bisulfite sequencing, together with promastigote and amastigote control lines. Interestingly, despite the DNMT6 presence, we found that methylation levels were equal to or lower than 0.0003% at CpG sites, 0.0005% at CHG sites, and 0.0126% at CHH sites at the genomic scale. As none of the methylated sites were retained after manual verification, we conclude that there is no evidence for DNA methylation in this species. We demonstrated that this difference in DNA methylation between the parasite (no detectable DNA methylation) and the vertebrate host (DNA methylation) allowed enrichment of parasite vs. host DNA using methyl-CpG-binding domain columns, readily available in commercial kits. As such, we depleted methylated DNA from mixes of Leishmania promastigote and amastigote DNA with human DNA, resulting in average Leishmania:human enrichments from 62× up to 263×. These results open a promising avenue for unmethylated DNA enrichment as a pre-enrichment step before sequencing Leishmania clinical samples.

Regulation of Leaf Longevity by DML3-Mediated DNA Demethylation

Leaf senescence is driven by the expression of senescence-associated genes (SAGs). Development-specific genes often undergo DNA demethylation in their promoter and other regions, which regulates gene expression. Whether and how DNA demethylation regulates the expression of SAGs and thus leaf senescence remain elusive. Whole-genome bisulfite sequencing (WGBS) analyses of wild-type (WT) and demeter-like 3 (dml3) Arabidopsis leaves at three developmental stages revealed hypermethylation during leaf senescence in dml3 compared with WT, and 20 556 differentially methylated regions (DMRs) were identified by comparing the methylomes of dml3 and WT in the CG, CHG, and CHH contexts. Furthermore, we identified that 335 DMR-associated genes (DMGs), such as NAC016 and SEN1, are upregulated during leaf senescence, and found an inverse correlation between the DNA methylation levels (especially in the promoter regions) and the transcript abundances of the related SAGs in WT. In contrast, in dml3 the promoters of SAGs were hypermethylated and their transcript levels were remarkably reduced, and leaf senescence was significantly delayed. Collectively, our study unraveled a novel epigenetic regulatory mechanism underlying leaf senescence in which DML3 is expressed at the onset of and during senescence to demethylate promoter, gene body or 3' UTR regions to activate a set of SAGs.

Characterizing the dynamic and functional DNA methylation landscape in the developing human cortex

DNA methylation (DNAm) is a key epigenetic regulator of gene expression across development. The developing prenatal brain is a highly dynamic tissue, but our understanding of key drivers of epigenetic variability across development is limited. We, therefore, assessed genomic methylation at over 39 million sites in the prenatal cortex using whole-genome bisulfite sequencing and found loci and regions in which methylation levels are dynamic across development. We saw that DNAm at these loci was associated with nearby gene expression and enriched for enhancer chromatin states in prenatal brain tissue. Additionally, these loci were enriched for genes associated with neuropsychiatric disorders and genes involved with neurogenesis. We also found autosomal differences in DNAm between the sexes during prenatal development, though these have less clear functional consequences. We lastly confirmed that the dynamic methylation at this critical period is specifically CpG methylation, with generally low levels of CpH methylation. Our findings provide detailed insight into prenatal brain development as well as clues to the pathogenesis of psychiatric traits seen later in life.

Genome-Wide DNA Methylation Enhances Stemness in the Mechanical Selection of Tumor-Repopulating Cells

Background

DNA methylation plays essential roles in tumor occurrence and stemness maintenance. Tumor-repopulating cells (TRCs) are cancer stem cell (CSC)-like cells with highly tumorigenic and self-renewing abilities, which were selected from tumor cells in soft three-dimensional (3D) fibrin gels.

Methods

Here, we presented a genome-wide map of methylated cytosines for time-series samples in TRC selection, in a 3D culture using whole-genome bisulfite sequencing (WGBS).

Results

A comparative analysis revealed that the methylation degrees of many differentially methylated genes (DMGs) were increased by the mechanical environment and changed from 2D rigid to 3D soft. DMGs were significantly enriched in stemness-related terms. In 1-day, TRCs had the highest non-CG methylation rate indicating its strong stemness. We found that genes with continuously increasing or decreasing methylation like CREB5/ADAMTS6/LMX1A may also affect the TRC screening process. Furthermore, results showed that stage-specific/common CSCs markers were biased toward changing their methylation in non-CG (CHG and CHH, where H corresponds to A, T, or C) methylation and enriched in gene body region.

Conclusions

WGBS provides DNA methylome in TRC screening. It was confirmed that non-CG DNA methylation plays an important role in TRC selection, which indicates that it is more sensitive to mechanical microenvironments and affects TRCs by regulating the expression of stemness genes in tumor cells.

Cytoplasm Types Affect DNA Methylation among Different Cytoplasmic Male Sterility Lines and Their Maintainer Line in Soybean (Glycine max L.)

Cytoplasmic male sterility (CMS) lines and their maintainer line have the same nucleus but different cytoplasm types. We used three soybean (Glycine max L.) CMS lines, JLCMS9A, JLCMSZ9A, and JLCMSPI9A, and their maintainer line, JLCMS9B, to explore whether methylation levels differed in their nuclei. Whole-genome bisulfite sequencing of these four lines was performed. The results show that the cytosine methylation level in the maintainer line was lower than in the CMS lines. Compared with JLCMS9B, the Gene Ontology (GO) enrichment analysis of DMR (differentially methylated region, DMR)-related genes of JLCMS9A revealed that their different 5-methylcytosine backgrounds were enriched in molecular function, whereas JLCMSZ9A and JLCMSPI9A were enriched in biological process and cellular component. The Kyoto Encyclopedia of Genes and Genome (KEGG) analysis of DMR-related genes and different methylated promoter regions in different cytosine contexts, hypomethylation or hypermethylation, showed that the numbers of DMR-related genes and promoter regions were clearly different. According to the DNA methylation and genetic distances separately, JLCMS9A clustered with JLCMS9B, and JLCMSPI9A with JLCMSZ9A. Thus, the effects of different cytoplasm types on DNA methylation were significantly different. This may be related to their genetic distances revealed by re-sequencing these lines. The detected DMR-related genes and pathways that are probably associated with CMS are also discussed.

Inter- and intratumor DNA methylation heterogeneity associated with lymph node metastasis and prognosis of esophageal squamous cell carcinoma

Background: Esophageal squamous cell carcinoma (ESCC), one of the leading causes of cancer mortality worldwide, is a heterogeneous cancer with diverse clinical manifestations. However, little is known about the epigenetic heterogeneity and its clinical relevance for this prevalent cancer. Methods: We generated 7.56 Tb single-base resolution whole-genome bisulfite sequencing data for 84 ESCC and paired paraneoplastic tissues. The analysis identified inter- and intratumor DNA methylation (DNAm) heterogeneity, epigenome-wide DNAm alterations together with the functional regulators involved in the hyper- or hypomethylated regions, and their association with clinical features. We then validated the correlation between the methylation level of specific regions and clinical outcomes of 96 ESCC patients in an independent cohort. Results: ESCC manifested substantial inter- and intratumor DNAm heterogeneity. The high intratumor DNAm heterogeneity was associated with lymph node metastasis and worse overall survival. Interestingly, hypermethylated regions in ESCC were enriched in promoters of numerous transcription factors, and demethylated noncoding regions related to RXR transcription factor binding appeared to contribute to the development of ESCC. Furthermore, we identified numerous DNAm alterations associated with carcinogenesis and lymph node metastasis of ESCC. We also validated three novel prognostic markers for ESCC, including one each in the promoter of CLK1, the 3' untranslated region of ZEB2, and the intergenic locus surrounded by several lncRNAs. Conclusions: This study presents the first population-level resource for dissecting base-resolution DNAm variation in ESCC and provides novel insights into the ESCC pathogenesis and progression, which might facilitate diagnosis and prognosis for this prevalent malignancy.

Base-Resolution Methylome of Retinal Pigment Epithelial Cells Used in the First Trial of Human Induced Pluripotent Stem Cell-Based Autologous Transplantation

The first-in-human trial of induced pluripotent stem cell (iPSC)-based autologous transplantation was successfully performed on a female patient with age-related macular degeneration. Here we delineated the base-resolution methylome of the iPSC-derived retinal pigment epithelium (iRPE) used in this trial. The methylome of iRPE closely resembled that of native RPE (nRPE), although partially methylated domains (PMDs) emerged in iRPE but not nRPE. Most differentially methylated regions between iRPE and nRPE appeared to originate from (de)methylation errors during differentiation, whereas errors at reprogramming resulted in aberrant genomic imprinting and X chromosome reactivation. Moreover, non-CpG methylation was prominent in nRPE but not iRPE. Intriguingly, xenotransplantation to mouse remodeled the iRPE methylome to demethylate a subset of suppressed genes and accumulate non-CpG methylation, but failed to resolve PMDs and hypermethylated CpG islands. Although the impacts of these alterations remain elusive, our findings should provide a useful guide for methylome analyses of other iPSC-derived cells.

•

The methylome of iPSC-derived RPE closely resembled that of native RPE

•

Most methylomic differences originated from errors during differentiation

•

Errors at reprogramming induced loss of imprinting and X chromosome reactivation

•

Some of the differences were mitigated by xenotransplantation to mouse

Differences in DNA Methylation Between Disease-Resistant and Disease-Susceptible Chinese Tongue Sole (Cynoglossus semilaevis) Families

DNA methylation, the most widely studied and most well-understood epigenetic modification, has been reported to play crucial roles in diverse processes. Although it has been found that DNA methylation can modulate the expression of immune-related genes in teleosts, a systemic analysis of epigenetic regulation on teleost immunity has rarely been performed. In this research, we employed whole-genome bisulfite sequencing to investigate the genome-wide DNA methylation profiles in select disease-resistant Cynoglossus semilaevis (DR-CS, family 14L006) and disease-susceptible C. semilaevis (DS-CS, family 14L104) against Vibrio harveyi infection. The results showed that following selective breeding, DR-CS had higher DNA methylation levels and different DNA methylation patterns, with 3,311 differentially methylated regions and 6,456 differentially methylated genes. Combining these data with the corresponding transcriptome data, we identified several immune-related genes that exhibited differential expression levels that were modulated by DNA methylation. Specifically, DNA methylation of tumor necrosis factor–like and lipopolysaccharide-binding protein-like was significantly correlated with their expression and significantly contributed to the disease resistance of the selected C. semilaevis family. In conclusion, we suggest that artificial selection for disease resistance in Chinese tongue sole causes changes in DNA methylation levels in important immune-related genes and that these epigenetic changes are potentially involved in multiple immune responses in Chinese tongue sole.

Integrating a genome-wide association study with a large-scale transcriptome analysis to predict genetic regions influencing the glycaemic index and texture in rice

Reliably generating rice varieties with low glycaemic index (GI) is an important nutritional intervention given the high rates of Type II diabetes incidences in Asia where rice is staple diet. We integrated a genome-wide association study (GWAS) with a transcriptome-wide association study (TWAS) to determine the genetic basis of the GI in rice. GWAS utilized 305 re-sequenced diverse indica panel comprising ~2.4 million single nucleotide polymorphisms (SNPs) enriched in genic regions. A novel association signal was detected at a synonymous SNP in exon 2 of LOC_Os05g03600 for intermediate-to-high GI phenotypic variation. Another major hotspot region was predicted for contributing intermediate-to-high GI variation, involves 26 genes on chromosome 6 (GI6.1). These set of genes included GBSSI, two hydrolase genes, genes involved in signalling and chromatin modification. The TWAS and methylome sequencing data revealed cis-acting functionally relevant genetic variants with differential methylation patterns in the hot spot GI6.1 region, narrowing the target to 13 genes. Conversely, the promoter region of GBSSI and its alternative splicing allele (G allele of Wxa ) explained the intermediate-to-high GI variation. A SNP (C˃T) at exon-10 was also highlighted in the preceding analyses to influence final viscosity (FV), which is independent of amylose content/GI. The low GI line with GC haplotype confirmed soft texture, while other two low GI lines with GT haplotype were characterized as hard and cohesive. The low GI lines were further confirmed through clinical in vivo studies. Gene regulatory network analysis highlighted the role of the non-starch polysaccharide pathway in lowering GI.

A non-randomized procedure for large-scale heterogeneous multiple discrete testing based on randomized tests

In the analysis of next-generation sequencing technology, massive discrete data are generated from short read counts with varying biological coverage. Conducting conditional hypothesis testing such as Fisher's Exact Test at every genomic region of interest thus leads to a heterogeneous multiple discrete testing problem. However, most existing multiple testing procedures for controlling the false discovery rate (FDR) assume that test statistics are continuous and become conservative for discrete tests. To overcome the conservativeness, in this article, we propose a novel multiple testing procedure for better FDR control on heterogeneous discrete tests. Our procedure makes decisions based on the marginal critical function (MCF) of randomized tests, which enables achieving a powerful and non-randomized multiple testing procedure. We provide upper bounds of the positive FDR (pFDR) and the positive false non-discovery rate (pFNR) corresponding to our procedure. We also prove that the set of detections made by our method contains every detection made by a naive application of the widely-used q-value method. We further demonstrate the improvement of our method over other existing multiple testing procedures by simulations and a real example of differentially methylated region (DMR) detection using whole-genome bisulfite sequencing (WGBS) data.

Genome-Wide DNA Methylation Comparison between Brassica napus Genic Male Sterile Line and Restorer Line

DNA methylation is an essential epigenetic modification that dynamically regulates gene expression during plant development. However, few studies have determined the DNA methylation profiles of male-sterile rapeseed. Here, we conducted a global comparison of DNA methylation patterns between the rapeseed genic male sterile line 7365A and its near-isogenic fertile line 7365B by whole-genome bisulfite sequencing (WGBS). Profiling of the genome-wide DNA methylation showed that the methylation level in floral buds was lower than that in leaves and roots. Besides, a total of 410 differentially methylated region-associated genes (DMGs) were identified in 7365A relative to 7365B. Traditional bisulfite sequencing polymerase chain reaction (PCR) was performed to validate the WGBS data. Eleven DMGs were found to be involved in anther and pollen development, which were analyzed by quantitative PCR. In particular, Bnams4 was hypo-methylated in 7365A, and its expression was up-regulated, which might affect other DMGs and thus control the male sterility. This study provided genome-wide DNA methylation profiles of floral buds and important clues for revealing the molecular mechanism of genic male sterility in rapeseed.

Cultured bovine embryo biopsy conserves methylation marks from original embryo

A major limitation of embryo epigenotyping by chromatin immunoprecipitation analysis is the reduced amount of sample available from an embryo biopsy. We developed an in vitro system to expand trophectoderm cells from an embryo biopsy to overcome this limitation. This work analyzes whether expanded trophectoderm (EX) is representative of the trophectoderm (TE) methylation or adaptation to culture has altered its epigenome. We took a small biopsy from the trophectoderm (30-40 cells) of in vitro produced bovine-hatched blastocysts and cultured it on fibronectin-treated plates until we obtained ∼4 × 104 cells. The rest of the embryo was allowed to recover its spherical shape and, subsequently, TE and inner cell mass were separated. We examined whether there were DNA methylation differences between TE and EX of three bovine embryos using whole-genome bisulfite sequencing. As a consequence of adaptation to culture, global methylation, including transposable elements, was higher in EX, with 5.3% of quantified regions showing significant methylation differences between TE and EX. Analysis of individual embryos indicated that TE methylation is more similar to its EX counterpart than to TE from other embryos. Interestingly, these similarly methylated regions are enriched in CpG islands, promoters and transcription units near genes involved in biological processes important for embryo development. Our results indicate that EX is representative of the embryo in terms of DNA methylation, thus providing an informative proxy for embryo epigenotyping.

Global, integrated analysis of methylomes and transcriptomes from laser capture microdissected bronchial and alveolar cells in human lung

Gene regulatory analysis of highly diverse human tissues in vivo is essentially constrained by the challenge of performing genome-wide, integrated epigenetic and transcriptomic analysis in small selected groups of specific cell types. Here we performed genome-wide bisulfite sequencing and RNA-seq from the same small groups of bronchial and alveolar cells isolated by laser capture microdissection from flash-frozen lung tissue of 12 donors and their peripheral blood T cells. Methylation and transcriptome patterns differed between alveolar and bronchial cells, while each of these epithelia showed more differences from mesodermally-derived T cells. Differentially methylated regions (DMRs) between alveolar and bronchial cells tended to locate at regulatory regions affecting promoters of 4,350 genes. A large number of pathways enriched for these DMRs including GTPase signal transduction, cell death, and skeletal muscle. Similar patterns of transcriptome differences were observed: 4,108 differentially expressed genes (DEGs) enriched in GTPase signal transduction, inflammation, cilium assembly, and others. Prioritizing using DMR-DEG regulatory network, we highlighted genes, e.g., ETS1, PPARG, and RXRG, at prominent alveolar vs. bronchial cell discriminant nodes. Our results show that multi-omic analysis of small, highly specific cells is feasible and yields unique physiologic loci distinguishing human lung cell types in situ.

Insights into Epigenome Evolution from Animal and Plant Methylomes

Evolutionary studies of DNA methylation offer insights into the mechanisms governing the variation of genomic DNA methylation across different species. Comparisons of gross levels of DNA methylation between distantly related species indicate that the size of the genome and the level of genomic DNA methylation are positively correlated. In plant genomes, this can be reliably explained by the genomic contents of repetitive sequences. In animal genomes, the role of repetitive sequences on genomic DNA methylation is less clear. On a shorter timescale, population-level comparisons demonstrate that genetic variation can explain the observed variability of DNA methylation to some degree. The amount of DNA methylation variation that has been attributed to genetic variation in the human population studies so far is substantially lower than that from Arabidopsis population studies, but this disparity might reflect the differences in the computational and experimental techniques used. The effect of genetic variation on DNA methylation has been directly examined in mammalian systems, revealing several causative factors that govern DNA methylation. On the other hand, studies from Arabidopsis have furthered our understanding of spontaneous mutations of DNA methylation, termed “epimutations.” Arabidopsis has an extremely high rate of spontaneous epimutations, which may play a major role in shaping the global DNA methylation landscape in this genome. Key missing information includes the frequencies of spontaneous epimutations in other lineages, in particular animal genomes, and how population-level variation of DNA methylation leads to species-level differences.

A MBD-seq protocol for large-scale methylome-wide studies with (very) low amounts of DNA

We recently showed that, after optimization, our methyl-CpG binding domain sequencing (MBD-seq) application approximates the methylome-wide coverage obtained with whole-genome bisulfite sequencing (WGB-seq), but at a cost that enables adequately powered large-scale association studies. A prior drawback of MBD-seq is the relatively large amount of genomic DNA (ideally >1 µg) required to obtain high-quality data. Biomaterials are typically expensive to collect, provide a finite amount of DNA, and may simply not yield sufficient starting material. The ability to use low amounts of DNA will increase the breadth and number of studies that can be conducted. Therefore, we further optimized the enrichment step. With this low starting material protocol, MBD-seq performed equally well, or better, than the protocol requiring ample starting material (>1 µg). Using only 15 ng of DNA as input, there is minimal loss in data quality, achieving 93% of the coverage of WGB-seq (with standard amounts of input DNA) at similar false/positive rates. Furthermore, across a large number of genomic features, the MBD-seq methylation profiles closely tracked those observed for WGB-seq with even slightly larger effect sizes. This suggests that MBD-seq provides similar information about the methylome and classifies methylation status somewhat more accurately. Performance decreases with <15 ng DNA as starting material but, even with as little as 5 ng, MBD-seq still achieves 90% of the coverage of WGB-seq with comparable genome-wide methylation profiles. Thus, the proposed protocol is an attractive option for adequately powered and cost-effective methylome-wide investigations using (very) low amounts of DNA.

Genome-Wide Analysis of DNA Methylation During Ovule Development of Female-Sterile Rice fsv1

The regulation of female fertility is an important field of rice sexual reproduction research. DNA methylation is an essential epigenetic modification that dynamically regulates gene expression during development processes. However, few reports have described the methylation profiles of female-sterile rice during ovule development. In this study, ovules were continuously acquired from the beginning of megaspore mother cell meiosis until the mature female gametophyte formation period, and global DNA methylation patterns were compared in the ovules of a high-frequency female-sterile line (fsv1) and a wild-type rice line (Gui99) using whole-genome bisulfite sequencing (WGBS). Profiling of the global DNA methylation revealed hypo-methylation, and 3471 significantly differentially methylated regions (DMRs) were observed in fsv1 ovules compared with Gui99. Based on functional annotation and Kyoto encyclopedia of genes and genomes (KEGG) pathway analysis of differentially methylated genes (DMGs), we observed more DMGs enriched in cellular component, reproduction regulation, metabolic pathway, and other pathways. In particular, many ovule development genes and plant hormone-related genes showed significantly different methylation patterns in the two rice lines, and these differences may provide important clues for revealing the mechanism of female gametophyte abortion.

A comparison of existing global DNA methylation assays to low-coverage whole-genome bisulfite sequencing for epidemiological studies

DNA methylation is an epigenetic mark at the interface of genetic and environmental factors relevant to human disease. Quantitative assessments of global DNA methylation levels have therefore become important tools in epidemiology research, particularly for understanding effects of environmental exposures in complex diseases. Among the available methods of quantitative DNA methylation measurements, bisulfite sequencing is considered the gold standard, but whole-genome bisulfite sequencing (WGBS) has previously been considered too costly for epidemiology studies with high sample numbers. Pyrosequencing of repetitive sequences within bisulfite-treated DNA has been routinely used as a surrogate for global DNA methylation, but a comparison of pyrosequencing to WGBS for accuracy and reproducibility of methylation levels has not been performed. This study compared the global methylation levels measured from uniquely mappable (non-repetitive) WGBS sequences to pyrosequencing assays of several repeat sequences and repeat assay-matched WGBS data and determined uniquely mappable WGBS data to be the most reproducible and accurate measurement of global DNA methylation levels. We determined sources of variation in repetitive pyrosequencing assays to be PCR amplification bias, PCR primer selection bias in methylation levels of targeted sequences, and inherent variability in methylation levels of repeat sequences. Low-coverage, uniquely mappable WGBS showed the strongest correlation between replicates of all assays. By using multiplexing by indexed bar codes, the cost of WGBS can be lowered significantly to improve the accuracy of global DNA methylation assessments for human studies.

Distinct Trends of DNA Methylation Patterning in the Innate and Adaptive Immune Systems

DNA methylation and the localization and post-translational modification of nucleosomes are interdependent factors that contribute to the generation of distinct phenotypes from genetically identical cells. With 112 whole-genome bisulfite sequencing datasets from the BLUEPRINT Epigenome Project, we analyzed the global development of DNA methylation patterns during lineage commitment and maturation of a range of immune system effector cells and the cancers that arise from them. We show clear trends in methylation patterns that are distinct in the innate and adaptive arms of the human immune system, both globally and in relation to consistently positioned nucleosomes. Most notable are a progressive loss of methylation in developing lymphocytes and the consistent occurrence of non-CG methylation in specific cell types. Cancer samples from the two lineages are further polarized, suggesting the involvement of distinct lineage-specific epigenetic mechanisms. We anticipate broad utility for this resource as a basis for further comparative epigenetic analyses.

Decoding the DNA Methylome of Mantle Cell Lymphoma in the Light of the Entire B Cell Lineage

We analyzed the in silico purified DNA methylation signatures of 82 mantle cell lymphomas (MCL) in comparison with cell subpopulations spanning the entire B cell lineage. We identified two MCL subgroups, respectively carrying epigenetic imprints of germinal-center-inexperienced and germinal-center-experienced B cells, and we found that DNA methylation profiles during lymphomagenesis are largely influenced by the methylation dynamics in normal B cells. An integrative epigenomic approach revealed 10,504 differentially methylated regions in regulatory elements marked by H3K27ac in MCL primary cases, including a distant enhancer showing de novo looping to the MCL oncogene SOX11. Finally, we observed that the magnitude of DNA methylation changes per case is highly variable and serves as an independent prognostic factor for MCL outcome.

Allele-Specific Methylome and Transcriptome Analysis Reveals Widespread Imprinting in the Human Placenta

DNA methylation is globally reprogrammed after fertilization, and as a result, the parental genomes have similar DNA-methylation profiles after implantation except at the germline differentially methylated regions (gDMRs). We and others have previously shown that human blastocysts might contain thousands of transient maternally methylated gDMRs (transient mDMRs), whose maternal methylation is lost in embryonic tissues after implantation. In this study, we performed genome-wide allelic DNA methylation analyses of purified trophoblast cells from human placentas and, surprisingly, found that more than one-quarter of the transient-in-embryo mDMRs maintained their maternally biased DNA methylation. RNA-sequencing-based allelic expression analyses revealed that some of the placenta-specific mDMRs were associated with expression of imprinted genes (e.g., TIGAR, SLC4A7, PROSER2-AS1, and KLHDC10), and three imprinted gene clusters were identified. This approach also identified some X-linked gDMRs. Comparisons of the data with those from other mammals revealed that genomic imprinting in the placenta is highly variable. These findings highlight the incomplete erasure of germline DNA methylation in the human placenta; understanding this erasure is important for understanding normal placental development and the pathogenesis of developmental disorders with imprinting effects.