Whole-genome Bisulfite Sequencing Research Articles

Elucidating the role of DNA methylation in low-salinity adaptation of the swimming crab Portunus trituberculatus

The swimming crab Portunus trituberculatus holds significant economic value in the aquaculture industry in China. Salinity is an important environmental factor that affects its growth and susceptibility to diseases. DNA methylation, a well-studied epigenetic mechanism, plays a pivotal role in the stress response of organisms. Herein we performed whole-genome bisulfite sequencing (WGBS) to assess changes in DNA methylation patterns in P. trituberculatus gill tissues before and after exposure to low-salinity stress. Overall, 33,712 differentially methylated regions (DMRs) were screened and 3546 functional genes were identified, many of which were notably enriched in pathways related to, for example, proteolysis and amino acid metabolism. By integrating WGBS and previous RNA-seq data, we identified 187 overlapping genes between differentially expressed genes and differentially methylated genes. Our findings reveal molecular mechanisms via which DNA methylation participates in low-salinity adaptation by regulating gene expression. Moreover, our data provide preliminary evidence of the substantial role of DNA methylation in facilitating salinity adaptation in P. trituberculatus and contribute to a systematic understanding of the mechanisms underpinning low-salinity adaptation in this species.

PCBS: an R package for fast and accurate analysis of bisulfite sequencing data.

Whole-genome bisulfite sequencing is a powerful tool for analyzing chromatin methylation genome-wide, but analysis of whole-genome bisulfite data is slow, inflexible, and often inaccurate. We developed PCBS (Principal Component BiSulfite), a computationally efficient R package for Whole Genome Bisulfite Sequencing analysis that demonstrates remarkable accuracy and flexibility compared to current tools. PCBS identifies differentially methylated loci, differentially methylated regions, and offers novel functionality that allows for more targeted methylation analyses. PCBS uses minimal computational resources; a complete pipeline in mouse can run on a local RStudio instance in a matter of minutes. PCBS is an R package available under a GNU GPLv3 license on GitHub: https://github.com/katlande/PCBS and CRAN: https://CRAN.R-project.org/package=PCBS. Instructions for use are available at: https://katlande.github.io/PCBS/.

DNA methylation profiling identifies TBKBP1 as potent amplifier of cytotoxic activity in CMV-specific human CD8+ T cells.

Epigenetic mechanisms stabilize gene expression patterns during CD8+ T cell differentiation. Although adoptive transfer of virus-specific T cells is clinically applied to reduce the risk of virus infection or reactivation in immunocompromised individuals, the DNA methylation pattern of virus-specific CD8+ T cells is largely unknown. Hence, we here performed whole-genome bisulfite sequencing of cytomegalovirus-specific human CD8+ T cells and found that they display a unique DNA methylation pattern consisting of 79 differentially methylated regions (DMRs) when compared to memory CD8+ T cells. Among the top demethylated DMRs in cytomegalovirus-specific CD8+ T cells was TBKBP1, coding for TBK-binding protein 1 that can interact with TANK-binding kinase 1 (TBK1) and mediate pro-inflammatory responses in innate immune cells downstream of intracellular virus sensing. Since TBKBP1 has not yet been reported in T cells, we aimed to unravel its role in virus-specific CD8+ T cells. TBKBP1 demethylation in terminal effector CD8+ T cells correlated with higher TBKBP1 expression at both mRNA and protein level, independent of alternative splicing of TBKBP1 transcripts. Notably, the distinct DNA methylation patterns in CD8+ T cell subsets was stable upon long-term in vitro culture. TBKBP1 overexpression resulted in enhanced TBK1 phosphorylation upon stimulation of CD8+ T cells and significantly improved their virus neutralization capacity. Collectively, our data demonstrate that TBKBP1 modulates virus-specific CD8+ T cell responses and could be exploited as therapeutic target to improve adoptive T cell therapies.

Rapid genome-wide profiling of DNA methylation and genetic variation using guide positioning sequencing (GPS).

Whole-genome bisulfite sequencing (WGBS) has been extensively utilized for DNA methylation profiling over the past decade. However, it has shown limitations in terms of high costs and inefficiencies. The productivity and accuracy of DNA methylation detection rely critically on the optimization of methodologies and the continuous refinements of related sequencing platforms. Here, we describe a detailed protocol of guide positioning sequencing (GPS), a bisulfite-based, location-specific sequencing technology designed for comprehensive DNA methylation characterization across the genome. The fundamental principle of GPS lies in the substitution of dCTP with 5-methyl-dCTP (5mC) at the 3'-end of DNA fragments by T4 DNA polymerase, which protects cytosines from bisulfite conversion to preserve the integrity of the base composition. This alteration allows the 3'-end to independently facilitate genetic variation profiling and guides the 5'-end, enriched with methylation information, to align more rapidly to the reference genome. Hence, GPS enables the concurrent detection of both genetic and epigenetic variations. Additionally, we provide an accessible description of the data processing, specifically involving certain software and scripts. Overall, the entire GPS procedure can be completed within a maximum of 15days, starting with a low initial DNA input of 100-500ng, followed by 4-5days for library construction, 8-10days for high-throughput sequencing (HTS) and data analysis, which can greatly facilitate the promotion and application of DNA methylation detection, especially for the rapid clinical diagnosis of diverse disease pathologies associated with concurrent genetic and epigenetic variations.

Whole-Genome Bisulfite Sequencing (WGBS) Analysis of Gossypium hirsutum under High-Temperature Stress Conditions

Background: DNA methylation is an important part of epigenetic regulation and plays an important role in the response of plants to adverse stress. Methods: In this study, whole-genome bisulfite sequencing (WGBS) was performed on the high-temperature-resistant material Xinluzao 36 and the high-temperature-sensitive material Che 61–72 at 0 h and 12 h under high-temperature stress conditions. Results: The results revealed that the Gossypium hirsutum methylation levels of CG and CHG (H = A, C, or T) decreased after the high-temperature stress treatment, and the methylation level of the A subgenome was significantly greater than that of the D subgenome. The methylation level of CHH increased, and the methylation level of CHH in the D subgenome was significantly greater than that in the A subgenome after high-temperature stress treatment. The methylation density of CG is lower than that of CHG and CHH, and the methylation density of the middle region of chromosomes is greater than that of both ends, which is opposite to the distribution density of genes. There were 124 common differentially methylated genes in the CG, CHG, and CHH groups, and 5130 common DEGs and differentially methylated genes were found via joint analysis with RNA-seq; these genes were significantly enriched in the biosynthesis of plant hormones, thiamine metabolism, glutathione metabolism, and tyrosine metabolism pathways. DNA methylation did not affect the expression of many genes (accounting for 85.68% of the differentially methylated genes), DNA methylation-promoted gene expression was located mainly in the downstream region of the gene or gene body, and the expression of inhibitory genes was located mainly in the upstream region of the gene. Conclusions: This study provides a theoretical basis for further exploration of the gene expression and functional regulatory mechanism of G. hirsutum DNA methylation under high-temperature stress conditions.

DNA methylation correlates with transcriptional noise in response to elevated pCO2 in the eastern oyster (Crassostrea virginica).

Ocean acidification significantly affects marine calcifiers like oysters, warranting the study of molecular mechanisms like DNA methylation that contribute to adaptive plasticity in response to environmental change. However, a consensus has not been reached on the extent to which methylation modules gene expression, and in turn plasticity, in marine invertebrates. In this study, we investigated the impact of pCO2 on gene expression and DNA methylation in the eastern oyster, Crassostrea virginica. After a 30-day exposure to control (572 ppm) or elevated pCO2 (2827 ppm), whole-genome bisulfite sequencing (WGBS) and RNA-seq data were generated from adult female gonad tissue and male sperm samples. Although differentially methylated loci (DMLs) were identified in females (89) and males (2916), there were no differentially expressed genes and only one differentially expressed transcript in females. However, gene body methylation impacted other forms of gene activity in sperm, such as the maximum number of transcripts expressed per gene and changes in the predominant transcript expressed. Elevated pCO2 exposure increased gene expression variability (transcriptional noise) in males but decreased noise in females, suggesting a sex-specific role of methylation in gene expression regulation. Functional annotation of genes with changes in transcript-level expression or containing DMLs revealed several enriched biological processes potentially involved in elevated pCO2 response, including apoptotic pathways and signal transduction, as well as reproductive functions. Taken together, these results suggest that DNA methylation may regulate gene expression variability to maintain homeostasis in elevated pCO2 conditions and could play a key role in environmental resilience in marine invertebrates.

438 DNA methylation patterns and transcriptional regulation during pig fetal skeletal muscle development

Abstract Fetal development is controlled by a complex cascade of highly regulated and coordinated gene expression patterns. Epigenetic mechanisms have important roles in regulating development and differentiation. Among such mechanisms, DNA methylation exhibits context-specific associations with gene expression and has been shown to be highly dynamic during developmental processes. We performed whole-genome bisulfite sequencing (WGBS) to assess DNA methylation in pig longissimus dorsi muscle at 41- and 70-d gestation (dg), as well as RNA- and small RNA-sequencing to identify coordinated changes in methylation and expression between myogenic stages. We identified 45,739 differentially methylated regions (DMRs) between stages, and the majority (n = 34,232) were hypomethylated at 70 vs. 41 dg. Developmental DMRs exhibited feature-specific enrichment in gene regulatory regions, as well as in regions proximal to micro-RNAs (miRNAs) that have known roles in myogenesis. Integration of methylation and transcriptomic data revealed strong associations between differential gene methylation and transcript abundance. We surveyed myogenic regulatory factor (MRF) genes to determine if differential methylation was present in expected genomic regions. Within the MYF5 and MYF6 locus, MYF5 was significantly promoter-hypermethylated at 70 dg, whereas MYF6 was significantly hypomethylated upstream of its transcription start site. MYF5 is the earliest MRF to be expressed and primarily functions in myoblast proliferation and determination, while MYF6 functions in muscle cell differentiation. Thus, these patterns were consistent with expected downregulation of MYF5 and upregulation of MYF6 as muscle development progresses and demonstrate that differential methylation is evident at myogenic transcription factors. Differential miRNA methylation was significantly negatively correlated with abundance, and dynamic expression of assayed miRNAs persisted postnatally. Motif analysis revealed significant enrichment of myogenic regulatory factor motifs among hypomethylated regions, suggesting that DNA hypomethylation may function to increase accessibility of muscle-specific transcription factors. We also show that developmental DMRs are enriched for GWAS SNPs associated with muscle physiology and meat quality traits, demonstrating the potential for epigenetic processes to influence phenotypic diversity. Our results enhance understanding of DNA methylation dynamics in pig fetal skeletal muscle and reveal putative cis-regulatory elements governed by epigenetic processes during porcine myogenesis.

Decoding the Root Lignification Mechanism of Angelica sinensis through Genome-wide Methylation Analysis.

Angelica sinensis is a traditional Chinese herbal medicine with significant economic and medicinal value. However, early bolting and flowering can occur during the second year of the vegetative growth period, rendering the roots unviable for medicinal use and resulting in substantial economic losses. Consequently, the growing interest in studying the molecular mechanisms underlying early bolting or increased root lignification in A. sinensis. Here, we conducted whole-genome bisulfite sequencing and observed an increase in whole-genome DNA methylation levels on chromosomes after bolting in A. sinensis. Comparative analysis revealed methylation patterns in the upstream, gene body, and downstream regions in the context of CG, CHG, and CHH, suggesting a possible association between CHH-type methylation of promoters and phenylpropanoid biosynthesis. Furthermore, joint analysis of transcriptomic and methylomics data revealed a positive correlation between DNA methylation and gene expression. We identified the hyperDMR gene in the CHH background within the promoter region; this gene is also a key gene (AsCOMT1), exhibiting dual catalytic activity and facilitating the synthesis of both ferulic acid and lignin. Enzyme kinetic analysis demonstrated that AsCOMT1 preferentially catalyzes the synthesis of lignin monomer precursors. These findings highlight the important regulatory role of DNA methylation in bolting and the synthesis of secondary metabolites in A. sinensis, providing valuable insights into the underlying molecular mechanisms. Therefore, as DNA methylation plays an important regulatory role in A. sinensis bolting and secondary metabolite synthesis, it has potential significance in the analysis of the underlying molecular mechanism.

Genome-wide analysis of DNA methylation and transcriptional changes associated with overwintering memory in Brassica rapa L. grown in the field

BackgroundWinter rapeseed, the sole overwintering oilseed crop in northern China, emphasizes winter resilience, yet epigenetic regulatory mechanisms governing overwintering memory remain poorly understood.ResultsIn this study, the root collar tissues from the robust cold-resistant variety Longyou-7 were sampled during the pre-winter period (S1), overwintering periods (S2–S5), and re-greening period (S6), to analyze overall genomic DNA methylation levels using high-performance liquid chromatography (HPLC). The result showed that DNA methylation level exceeded 80% in the S1 stage. Throughout the overwintering periods, methylation levels displayed a decreasing trend in S3, followed by an increase in S5, and a pronounced decrease in S6. Consequently, S1, S3, S5, and S6 periods were chosen for whole-genome bisulfite sequencing analyses to elucidate the overwintering memory mechanisms of Longyou-7. The result revealed that DNA methylation primarily occurs in the CG context in Longyou-7. However, methylation of mC sites is most prevalent in the CHH type, gradually decreasing during overwintering periods. Analysis of methylation patterns in specific genomic regions of Longyou-7 showed that the highest methylation levels in the intergenic region. Moreover, mC sites in repeats and transposon elements are distributed differently across the three contexts. Subsequently, differentially methylated regions and promoters of Longyou-7 were identified during various periods compared to the S1 stage, followed by joint analysis with transcriptome sequencing. Functional enrichment analysis highlighted the involvement of most overlapping genes in the MAPK signaling pathway, plant hormone signal transduction, and starch and sucrose metabolism pathways. Changes in candidate gene expression within these three pathways correlated closely with DNA methylation levels.ConclusionsOur findings underscored the critical role of DNA methylation in regulating the expression of overwintering memory genes in winter rapeseed. These results offer a comprehensive insights into the epigenetic regulatory mechanisms governing winter rapeseed's overwintering memory, while identified overwintering memory genes served as crucial genetic resources for multifaceted breeding of winter-resistant varieties.Graphical

Early moderate prenatal alcohol exposure and maternal diet impact offspring DNA methylation across species.

Alcohol consumption in pregnancy can affect genome regulation in the developing offspring but results have been contradictory. We employed a physiologically relevant murine model of short-term moderate prenatal alcohol exposure (PAE) resembling common patterns of alcohol consumption in pregnancy in humans. Early moderate PAE was sufficient to affect site-specific DNA methylation in newborn pups without altering behavioural outcomes in adult littermates. Whole-genome bisulfite sequencing of neonatal brain and liver revealed stochastic influence on DNA methylation that was mostly tissue-specific, with some perturbations likely originating as early as gastrulation. DNA methylation differences were enriched in non-coding genomic regions with regulatory potential indicative of broad effects of alcohol on genome regulation. Replication studies in human cohorts with fetal alcohol spectrum disorder suggested some effects were metastable at genes linked to disease-relevant traits including facial morphology, intelligence, educational attainment, autism, and schizophrenia. In our murine model, a maternal diet high in folate and choline protected against some of the damaging effects of early moderate PAE on DNA methylation. Our studies demonstrate that early moderate exposure is sufficient to affect fetal genome regulation even in the absence of overt phenotypic changes and highlight a role for preventative maternal dietary interventions.

De novo Whole-Genome Assembly of the 10-Gigabase Fokienia Hodginsii Genome to Reveal Differential Epigenetic Events Between Callus and Xylem.

Fokienia hodginsii (F. hodginsii), belonging to the genus Fokienia of the Cupressaceae. F. hodginsii has significant application value due to its wood properties and great research value in evolutionary studies as a gymnosperm. However, the genome of F. hodginsii remains unknown due to the large size of gymnosperms genome. Pacific Bioscience sequencing, Hi-C mapping, whole-genome Bisulfite Sequencing (BS-Seq), long-read isoform sequencing (Iso-Seq), direct RNA sequencing (DRS), quantitative proteomics, and metabonomics analysis are employed to facilitate genome assembly, gene annotation, and investigation into epigenetic mechanisms. In this study, the 10G F. hodginsii genome is assembled into 11 chromosomes. Furthermore, 50521 protein-coding genes are annotated and determined that 65% of F. hodginsii genome comprises repetitive sequences. It is discovered that transposable element (TE)-including introns is associated with higher expression. The DNA methylome of F. hodginsii reveals that xylem has a higher DNA methylation level compared to callus. Moreover, DRS reveals the significant alterations in RNA full-length ratio, which potentially associated with poly(A) length (PAL) and alternative polyadenylation (APA). Finally, the morphology measurement and metabonomics analysis revealed the difference of 14 cultivars. In summary, the genomes and epigenetics datasets provide a molecular basis for callus formation in the gymnosperm family.

DMAP2: A Pipeline for Analysis of Whole-Genome-Scale DNA Methylation Sequencing Data.

DNA methylation is well-established as a major epigenetic mechanism that can control gene expression and is involved in both normal development and disease. Analysis of high-throughput-sequencing-based DNA methylation data is a step toward understanding the relationship between disease and phenotype. Analysis of CpG methylation at single-base resolution is routinely done by bisulfite sequencing, in which methylated Cs remain as C while unmethylated Cs are converted to U, subsequently seen as T nucleotides. Sequence reads are aligned to the reference genome using mapping tools that accept the C-T ambiguity. Then, various statistical packages are used to identify differences in methylation between (groups of) samples. We have previously developed the Differential Methylation Analysis Pipeline (DMAP) as an efficient, fast, and flexible tool for this work, both for whole-genome bisulfite sequencing (WGBS) and reduced-representation bisulfite sequencing (RRBS). The protocol described here includes a series of scripts that simplify the use of DMAP tools and that can accommodate the wider range of input formats now in use to perform analysis of whole-genome-scale DNA methylation sequencing data in various biological and clinical contexts. © 2024 The Author(s). Current Protocols published by Wiley Periodicals LLC. Basic Protocol: DMAP2 workflow for whole-genome bisulfite sequencing (WGBS) and reduced-representation bisulfite sequencing (RRBS).

Role of Epigenetic Factors in Determining the Biological Behavior and Prognosis of Hepatocellular Carcinoma.

The current study's objective is to evaluate the molecular genetic mechanisms influencing the biological behavior of hepatocellular carcinoma (HCC) by analyzing the transcriptomic and epigenetic signatures of the tumors. Transcriptomic data were downloaded from the NCBI GEO database. We investigated the expression differences between the GSE46444 (48 cirrhotic tissues versus 88 HCC tissues) and GSE63898 (168 cirrhotic tissues versus 228 HCC tissues) data sets using GEO2R. Differentially expressed genes were evaluated using GO and KEGG metabolic pathway analysis websites. Whole genome bisulfite sequencing (WGBS) and Methylated DNA Immunoprecipitation Sequencing (MeDIP-Seq) data sets (26 HCC tissues versus 26 adjacent non-tumoral tissues) were also downloaded from the NCBI SRA database. These data sets were analyzed using Bismark and QSEA, respectively. The methylation differences between the groups were assessed using functional enrichment analysis. In the GSE46444 data set, 80 genes were upregulated, and 315 genes were downregulated in the tumor tissue (HCC tissue) compared to the non-tumor cirrhotic tissue. In the GSE63898 data set, 1261 genes were upregulated, and 458 genes were downregulated in the cirrhotic tissue compared to the tumor tissues. WGBS revealed that 20 protein-coding loci were hypermethylated. while the hypomethylated regions were non-protein-coding. The methylated residues of the tumor tissue, non-tumorous cirrhotic tissue, and healthy tissue were comparable. MeDIP-Seq, conducted on tumoral and non-tumoral tissues, identified hypermethylated or hypomethylated areas as protein-coding regions. The functional enrichment analysis indicated that these genes were related to pathways including peroxisome, focal adhesion, mTOR, RAP1, Phospholipase D, Ras, and PI3K/AKT signal transduction. The investigation of transcriptomic and epigenetic mechanisms identified several genes significant in the biological behavior of HCC. These genes present potential targets for the development of targeted therapy.

Genome-wide identification of DNA methylation marks associated with resistance to Streptococcus agalactiae in the GIFT strain of Nile Tilapia (Oreochromis niloticus)

Epigenetic mechanism plays a pivotal role in the immune responses of fish. Streptococcus agalactiae (SA) infection is a significant challenge to tilapia aquaculture. DNA methylation involved in immune responses in tilapia against SA is unclear. Here, by using a QTLseq strategy integrated with whole-genome bisulfite sequencing-based DNA methylation analysis (methyl-QTLseq) and RNA sequencing (RNA-seq), we identified 3608 differentially methylated regions (DMRs) and 4355 differentially expressed genes (DEGs) in the spleen of Nile tilapia Oreochromis niloticus. Compared to the susceptible group, the up-regulated DMRs (2449) were approximately twice as many as the down-regulated DMRs (1159). Compared to the susceptible group, lower methylation levels for CpG, CHG and CHH in TEs and for CHG and CHH in genes and significant decreased expression for 8 of the 11 DNA methyltransferase genes in the resistant group were detected. 46 DEGs with differentially methylated promoter regions (DMPs) were significantly enriched in immune processes such as T cell differentiation and adaptive immunity. We detected 4 methyl-QTLs located on LG8:5332501–5,336,400, LG14:2307301–2,353,500, LG16:16889101–16,903,500, and LG22:25238401–2,536,880 across the tilapia genome. The methyl-QTL region located within the chromosome LG22:25238401–25,368,800 contains the most DMRs (91) and DMGs (7). Interestingly, most of the annotated genes that located within or under the peak of the methyl-QTL interval belong to the HOX protein family including hoxa1, hoxa2, hoxa4a and hoxa5a. Bisulfite sequencing PCR and expression analysis indicated that DNA methylation in the promoter regions of hoxa4a, lcp2a, and ccr9a regulated gene expression in spleen in response to SA infection. Overall, our study provides epigenetics evidences involving in immune responses against SA infection in tilapia.

Maternal Undernutrition Affects Fetal Thymus DNA Methylation, Gene Expression, and, Thereby, Metabolism and Immunopoiesis in Wagyu (Japanese Black) Cattle.

We aimed to determine the effects of maternal nutrient restriction (MNR) on the DNA methylation and gene expression patterns associated with metabolism and immunopoiesis in the thymuses of fetal Wagyu cattle. Pregnant cows were allocated to two groups: a low-nutrition (LN; 60% nutritional requirement; n = 5) and a high-nutrition (HN; 120% nutritional requirement, n = 6) group, until 8.5 months of gestation. Whole-genome bisulfite sequencing (WGBS) and RNA sequencing were used to analyze DNA methylation and gene expression, while capillary electrophoresis-Fourier transform mass spectrometry assessed the metabolome. WGBS identified 4566 hypomethylated and 4303 hypermethylated genes in the LN group, with the intergenic regions most frequently being methylated. Pathway analysis linked hypoDMGs to Ras signaling, while hyperDMGs were associated with Hippo signaling. RNA sequencing found 94 differentially expressed genes (66 upregulated, 28 downregulated) in the LN group. The upregulated genes were tied to metabolic pathways and oxidative phosphorylation; the downregulated genes were linked to natural killer cell cytotoxicity. Key overlapping genes (GRIA1, CACNA1D, SCL25A4) were involved in cAMP signaling. The metabolomic analysis indicated an altered amino acid metabolism in the MNR fetuses. These findings suggest that MNR affects DNA methylation, gene expression, and the amino acid metabolism, impacting immune system regulation during fetal thymus development in Wagyu cattle.

Cashmere cyclic growth affected by different photoperiods alters DNA methylation patterns

Whole-genome bisulfite sequencing was used to study epigenetic differences in hair follicles from Inner Mongolia Albas cashmere goats during long and short photoperiods. Results showed earlier anagen entry in short photoperiod. Genome-wide cytosine methylation was 3.23%, with CG, CHG, and CHH methylation rates at 48.00%, 0.19%, and 0.17%, respectively. Methylation was higher for CG than for CHG and CHH, and higher for introns than exons and downstream regions. Both groups had similar levels of CG, CHG, and CHH methylation, but CHH methylation levels differed slightly between photoperiods. A total of 9,700 differentially methylated regions (DMRs) and 4,187 genes (DMGs) were identified. Gene Ontology (GO) and KEGG analyses revealed significant enrichment in 67 GO terms and 41 signaling pathways related to circadian rhythms, bringing new insights in the role of DNA methylation in cashmere growth and development.

Comprehensive Analysis of Methylome and Transcriptome to Identify Potential Genes Regulating Porcine Testis Development.

DNA methylation plays a critical role in regulating gene expression during testicular development. However, few studies report on candidate genes related to the DNA methylation regulation of porcine testicular development. This study examined the differentially expressed genes (DEGs) and their methylation levels in testicular tissues from pigs at 60 days of age (60 d) and 180 days of age (180 d) using RNA-Seq and whole genome bisulfite sequencing (WGBS). It was determined that DNA methylation primarily occurs in the cytosine-guanine (CG) context, and the analysis identified 106,282 differentially methylated regions (DMRs) corresponding to 12,385 differentially methylated genes (DMGs). Further integrated analysis of RNA-Seq and WGBS data revealed 1083 DMGs negatively correlated with the expression of DEGs. GO analysis showed that these genes were significantly enriched in spermatogenesis, germ cell development, and spermatid differentiation. The screening of enriched genes revealed that hyper-methylation repressed ADAM30, ADAM3A, DPY19L2, H2BC1, MAK, RPL10L, SPATA16, and YBX2, while hypo-methylation elevated CACNA1I, CADM1, CTNNB1, JAM2, and PAFAH1B3 expression. Additionally, the methylation status of the key genes ADAM3A, ADAM30, YBX2, JAM2, PAFAH1B3, and CTNNB1 was detected by bisulfite sequencing PCR (BSP). This study offers insights into the epigenetic regulation mechanisms underlying porcine testicular development.

Genomic context determines the effect of DNA methylation on gene expression in the gut epithelium of Atlantic salmon (Salmo salar)

ABSTRACT The canonical view of DNA methylation, a pivotal epigenetic regulation mechanism in eukaryotes, dictates its role as a suppressor of gene activity, particularly within promoter regions. However, this view is being challenged as it is becoming increasingly evident that the connection between DNA methylation and gene expression varies depending on the genomic location and is therefore more complex than initially thought. We examined DNA methylation levels in the gut epithelium of Atlantic salmon (Salmo salar) using whole-genome bisulfite sequencing, which we correlated with gene expression data from RNA sequencing of the same gut tissue sample (RNA-seq). Assuming epigenetic signals might be pronounced between distinctive phenotypes, we compared large and small fish, finding 22 significant associations between 22 differentially methylated regions and 21 genes. We did not detect significant methylation differences between large and small fish. However, we observed a consistent signal of methylation levels around the transcription start sites (TSS), being negatively correlated with the expression levels of those genes. We found both negative and positive associations of methylation levels with gene expression further upstream or downstream of the TSS, revealing a more unpredictable pattern. The 21 genes showing significant methylation-expression correlations were involved in biological processes related to salmon health, such as growth and immune responses. Deciphering how DNA methylation affects the expression of such genes holds great potential for future applications. For instance, our results suggest the importance of genomic context in targeting epigenetic modifications to improve the welfare of aquaculture species like Atlantic salmon.

DNA methylation analysis reveals local changes in resistant and susceptible soybean lines in response to Phytophthora sansomeana.

Phytophthora sansomeana is an emerging oomycete pathogen causing root rot in many agricultural species including soybean. However, as of now, only one potential resistance gene has been identified in soybean, and our understanding of how genetic and epigenetic regulation in soybean contributes to responses against this pathogen remains largely unknown. In this study, we performed whole genome bisulfite sequencing (WGBS) on two soybean lines, Colfax (resistant) and Williams 82 (susceptible), in response to P. sansomeana at two time points: 4 and 16 hours post-inoculation to compare their methylation changes. Our findings revealed that there were no significant changes in genome-wide CG, CHG (H = A, T, or C), and CHH methylation. However, we observed local methylation changes, specially an increase in CHH methylation around genes and transposable elements (TEs) after inoculation, which occurred earlier in the susceptible line and later in the resistant line. After inoculation, we identified differentially methylated regions (DMRs) in both Colfax and Williams 82, with a predominant presence in TEs. Notably, our data also indicated that more TEs exhibited changes in their methylomes in the susceptible line compared to the resistant line. Furthermore, we discovered 837 DMRs within or flanking 772 differentially expressed genes (DEGs) in Colfax and 166 DMRs within or flanking 138 DEGs in Williams 82. These DEGs had diverse functions, with Colfax primarily showing involvement in metabolic process, defense response, plant and pathogen interaction, anion and nucleotide binding, and catalytic activity, while Williams 82 exhibited a significant association with photosynthesis. These findings suggest distinct molecular responses to P. sansomeana infection in the resistant and susceptible soybean lines.

Discovery and Validation of Methylation Signatures in Circulating Cell-Free DNA for the Detection of Colorectal Cancer.

This study was conducted with the primary objective of assessing the performance of cfDNA methylation in the detection of colorectal cancer (CRC). Five tumor tissue, 20 peripheral blood leucocyte, and 169 cfDNA samples were collected for whole-genome bisulfite sequencing (WGBS) analysis. Bioinformatic analysis was conducted to identify differentially methylated regions (DMRs) and their functional characteristics. Quantitative methylation-specific PCR (qMSP) was used to validate the methylation levels of DMRs in the tissues and leucocytes. cfDNA samples from CRC patients and healthy controls were used to evaluate the performance of the DMR analysis. WGBS analysis revealed a decrease in DNA methylation levels in the CpG context in CRC tumor tissues compared with adjacent normal tissues. A total of 132 DMRs in cfDNA were identified as potential markers for diagnosing CRC. In a cohort of 95 CRC patients and 74 healthy controls, a combination of the three DMRs (DAB1, PPP2R5C, and FAM19A5) yielded an AUC of 0.763, achieving 64.21% sensitivity and 78.38% specificity in discriminating CRC patients from healthy controls. This study provides insights into DNA methylation patterns in CRC and identifies a set of DMRs in cfDNA with potential diagnostic value for CRC. These DMRs hold promise as biomarkers for CRC detection, offering promise for non-invasive CRC diagnosis. Further research is warranted to validate these findings in larger cohorts.