Integrated study of DNA methylation and transcriptome: a new perspective for exploring the pathogenesis of Sjögren's syndrome

Epigenetic modifications like DNA methylation can alter an organism's phenotype without changing its DNA sequence. Exposure to environmental toxicants has the potential to change the resilience of aquatic species. However, little information is available on the dynamics of DNA methylation in fish gonadal tissues in response to organophosphates. In the present work, reduced-representation bisulfite sequencing was performed to identify DNA methylation patterns in the ovarian tissues of Anabas testudienus exposed to organophosphates, specifically monocrotophos (MCP). Through sequencing, an average of 41,087 methylated cytosine sites were identified and distributed in different parts of genes, i.e., in transcription start sites (TSS), promoters, exons, etc. A total of 1058 and 1329 differentially methylated regions (DMRs) were detected as hyper-methylated and hypo-methylated in ovarian tissues, respectively. Utilizing whole-genome data of the climbing perch, the DMRs, and their associated overlapping genes revealed a total of 22 genes within exons, 45 genes at transcription start sites (TSS), and 218 genes in intergenic regions. Through gene ontology analysis, a total of 16 GO terms particularly involved in ovarian follicular development, response to oxidative stress, oocyte maturation, and multicellular organismal response to stress associated with reproductive biology were identified. After functional enrichment analysis, relevant DMGs such as steroid hormone biosynthesis (Cyp19a, 11-beta-HSD, 17-beta-HSD), hormone receptors (ar, esrrga), steroid metabolism (StAR), progesterone-mediated oocyte maturation (igf1ar, pgr), associated with ovarian development in climbing perch showed significant differential methylation patterns. The differentially methylated genes (DMGs) were subjected to analysis using real-time PCR, which demonstrated altered gene expression levels. This study revealed a molecular-level alteration in genes associated with ovarian development in response to chemical exposure. This work provides evidence for understanding the relationship between DNA methylation and gene regulation in response to chemicals that affect the reproductive fitness of aquatic animals.

BackgroundDNA methylation is an important epigenetic modification. It can regulate the expression of many key genes without changing the primary structure of the genomic DNA, and plays a vital role in the growth and development of the organism. The genome-wide DNA methylation profile of the cytoplasmic male sterile (CMS) line in soybean has not been reported so far.ResultsIn this study, genome-wide comparative analysis of DNA methylation between soybean CMS line NJCMS5A and its maintainer NJCMS5B was conducted by whole-genome bisulfite sequencing. The results showed 3527 differentially methylated regions (DMRs) and 485 differentially methylated genes (DMGs), including 353 high-credible methylated genes, 56 methylated genes coding unknown protein and 76 novel methylated genes with no known function were identified. Among them, 25 DMRs were further validated that the genome-wide DNA methylation data were reliable through bisulfite treatment, and 9 DMRs were confirmed the relationship between DNA methylation and gene expression by qRT-PCR. Finally, 8 key DMGs possibly associated with soybean CMS were identified.ConclusionsGenome-wide DNA methylation profile of the soybean CMS line NJCMS5A and its maintainer NJCMS5B was obtained for the first time. Several specific DMGs which participated in pollen and flower development were further identified to be probably associated with soybean CMS. This study will contribute to further understanding of the molecular mechanism behind soybean CMS.

1. Valgus-varus deformity (VVD) is a common leg bone problem in broilers that causes serious economic losses to the breeding industry. The genetic aetiology of VVD is not clear, which restricts the genetic control of VVD. 2. In this study, knee cartilage of 35-day-old VVD and normal broilers was sequenced by whole-genome bisulphite sequencing (WGBS). The unique whole-genome DNA methylation profile of VVD broilers was described, and the methylation data and transcription data were used for joint analysis. 3. The mean methylation level of the VVD group was greater than that in the normal group. A total of 4315 differentially methylated regions (DMRs) were detected from methylation data, with the highest DMR density on chromosomes 25, 27, 31 and 33. DMRs were mainly located in introns, which accounted for more than 60%, followed by promoter and exon regions. 4. A total of 2326 differentially methylated genes (DMGs) were identified from DMRs, including 1159 genes with upregulated DMRs, 936 genes with downregulated DMRs, and 231 genes with two types of DMRs. 5. The ESPL1 gene may be an important epigenetic gene of VVD. The methylation of particular CpG17, CpG18 and CpG19 sites in the promoter region of the ESPL1 gene may hinder the binding of transcription factors and promoters and increase the expression of ESPL1.

Reconstructing Denisovan Anatomy Using DNA Methylation Maps.

Ovarian tissue is critical for goose reproduction. This study aimed to investigate gene regulation by DNA methylation in relation to the reproductive traits of geese. We performed whole-genome bisulfite sequencing (WGBS) on ovarian tissues from Sichuan white geese (high-laying-rate group: HLRG, ♀ = 3; low-laying-rate group: LLRG, ♀ = 3) during the laying period. The results showed a higher level of hypermethylated differentially methylated regions (DMRs) in the HLRG, indicating a higher overall methylation level compared to the LLRG. In total, we identified 2831 DMRs and 733 differentially methylated genes (DMGs), including 363 genes with upregulated methylation. These DMGs were significantly enriched in pathways related to microtubule function (GO:0005874; GO:0000226), GnRH secretion, thyroid hormone signaling, ECM-receptor interaction, and PI3K-Akt signaling. Integration with RNA-seq data identified eight overlapping genes between DMGs and differentially expressed genes (DEGs), with five genes (CUL9, MEGF6, EML6, SYNE2, AK1BA) exhibiting a correlation between hypomethylation and high expression. EML6, in particular, emerged as a promising candidate, potentially regulating follicle growth and development in Sichuan white geese. Future studies should focus on further verifying the role of the EML6 gene. In conclusion, this study provides important insights into the regulatory mechanisms of DNA methylation influencing reproductive traits in geese, offering novel candidate markers for future goose breeding programs.

Deviant DNA methylation play a key role in the pollen abortion of Gossypium barbadense L. cytoplasmic male sterility

A high-fat (HF) diet is associated with progression of liver diseases. To illustrate genome-wide landscape of DNA methylation in liver of rats fed either a control or HF diet, two enrichment-based methods, namely methyl-DNA immunoprecipitation assay with high-throughput sequencing (MeDIP-seq) and methylation-sensitive restriction enzyme sequencing (MRE-seq), were performed in our study. Rats fed with the HF diet exhibited an increased body weight and liver fat accumulation compared with that of the control group when they were 12 wk of age. Genome-wide analysis of differentially methylated regions (DMRs) showed that 12,494 DMRs induced by HF diet were hypomethylated and 6,404 were hypermethylated. DMRs with gene annotations [differentially methylated genes (DMGs)] were further analyzed to show gene-specific methylation profile. There were 88, 2,680, and 95 hypomethylated DMGs identified with changes in DNA methylation in the promoter, intragenic and downstream regions, respectively, compared with fewer hypermethylated DMGs (45, 1,623, and 50 in the respective regions). Some of these genes also contained an ACGT cis-acting motif whose DNA methylation status may affect gene expression. Pathway analysis showed that these DMGs were involved in critical hepatic signaling networks related to hepatic development. Therefore, HF diet had global impacts on DNA methylation profile in the liver of rats, leading to differential expression of genes in hepatic pathways that may involve in functional changes in liver development.

BackgroundBean pyralid is one of the major leaf-feeding insects that affect soybean crops. DNA methylation can control the networks of gene expressions, and it plays an important role in responses to biotic stress. However, at present the genome-wide DNA methylation profile of the soybean resistance to bean pyralid has not been reported so far.ResultsUsing whole-genome bisulfite sequencing (WGBS) and RNA-sequencing (RNA-seq), we analyzed the highly resistant material (Gantai-2-2, HRK) and highly susceptible material (Wan82–178, HSK), under bean pyralid larvae feeding 0 h and 48 h, to clarify the molecular mechanism of the soybean resistance and explore its insect-resistant genes. We identified 2194, 6872, 39,704 and 40,018 differentially methylated regions (DMRs), as well as 497, 1594, 9596 and 9554 differentially methylated genes (DMGs) in the HRK0/HRK48, HSK0/HSK48, HSK0/HRK0 and HSK48/HRK48 comparisons, respectively. Through the analysis of global methylation and transcription, 265 differentially expressed genes (DEGs) were negatively correlated with DMGs, there were 34, 49, 141 and 116 negatively correlated genes in the HRK0/HRK48, HSK0/HSK48, HSK0/HRK0 and HSK48/HRK48, respectively. The MapMan cluster analysis showed that 114 negatively correlated genes were clustered in 24 pathways, such as protein biosynthesis and modification; primary metabolism; secondary metabolism; cell cycle, cell structure and component; RNA biosynthesis and processing, and so on. Moreover, CRK40; CRK62; STK; MAPK9; L-type lectin-domain containing receptor kinase VIII.2; CesA; CSI1; fimbrin-1; KIN-14B; KIN-14 N; KIN-4A; cytochrome P450 81E8; BEE1; ERF; bHLH25; bHLH79; GATA26, were likely regulatory genes involved in the soybean responses to bean pyralid larvae. Finally, 5 DMRs were further validated that the genome-wide DNA data were reliable through PS-PCR and 5 DEGs were confirmed the relationship between DNA methylation and gene expression by qRT-PCR. The results showed an excellent agreement with deep sequencing.ConclusionsGenome-wide DNA methylation profile of soybean response to bean pyralid was obtained for the first time. Several specific DMGs which participated in protein kinase, cell and organelle, flavonoid biosynthesis and transcription factor were further identified to be likely associated with soybean response to bean pyralid. Our data will provide better understanding of DNA methylation alteration and their potential role in soybean insect resistance.

Differential DNA Methylation Predicts Response To Combined Treatment Regimens With a DNA Methyltransferase Inhibitor In Acute Myeloid Leukemia (AML)

S-adenosylmethionine Levels Regulate the Schwann Cell DNA Methylome

To study DNA methylation patterns of non-syndromic cleft lip/palate(NSCL/P) using bioinformatic methods, including methylated positions and regions. Whole blood DNA methylation data of NSCL/P samples was download from Gene Expression Omnibus(GEO) database, including 67 NSCL/P cases and 59 controls without birth defects. Data analysis included ①data cleaning, such as probes filtering, quality control and normalization; ②differential methylation analysis, including methylated positions and regions; ③Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis on differential methylated genes. R 3.4.3 software was used for data cleaning, differential methylated positions and regions analysis. DAVID6.8 tool was used for GO and KEGG analysis. 814 differential methylated positions were detected (adjusted P＜0.001，|Δβ|＞0.125）, of which 178 were hypermethylated in NSCL/P patients, and 636 were hypomethylated. In addition, 386 differential methylated regions were identified (P＜0.05), of which 204 were hypermethylation regions and 182 were hypomethylation regions. GO analysis showed that 38 differential methylated genes were involved in 7 kinds of biological processes, 163 differential methylated genes were involved in 3 kinds of molecular functions, and 114 differential methylated genes were involved in 3 kinds of cellular components (P＜0.01). KEGG analysis showed that 59 differential methylated genes were involved in 9 kinds of signal pathways. Abnormal DNA methylation patterns of NSCL/P might be an important epigenetic mechanism affecting the development of NSCL/P. This study might contribute to the identification of identification of biomarkers and targeted interventions of NSCL/P.

Androgenesis is an important chromosome set manipulation technique used in sex control in aquaculture. Haploid embryos exhibit haploid syndrome with body abnormalities and even die during early embryonic development. In this study, we used whole genome bisulfite sequencing (WGBS) to investigate the genome-wide DNA methylation profiles in haploid females (1n-X) and males (1n-Y), and diploid females (2n-XX) and males (2n-XY) of tiger pufferfish (Takifugu rubripes), an economically important fish in China. A total of 96.32 Gb clean data was produced. Differentially methylated regions (DMRs) were found between haploids and diploids, which may be related to abnormal development and early embryonic death in haploids. There were 3,641 hyper-methylated differentially methylated genes (DMGs) and 2,179 hypo-methylated DMGs in haploid vs. diploid comparisons in both females and males. These DMGs were mainly related to genomic stability maintenance and cell cycle regulation. slf1, actr8, gas2, and pbrm1 genes were selected to validate the methylation sequencing. After combining the methylation data with the corresponding transcriptome data, we identified several genes, including guca2a, myoc, fezf2, rprml, telo2, s100a1, and marveld1, which exhibited differential expression levels modulated by DNA methylation. In conclusion, our study revealed different methylation and expression profiles between haploid and diploid T. rubripes for the first time. Several DMGs were identified between different ploidy levels, which may be related to haploid syndrome formation. The results expand the understanding of the effects of ploidy on the early development of teleosts and provide knowledge about target genes and networks to improve the survival rate of haploids.

DNA methylation is a form of epigenetic regulation, having pivotal parts in controlling cellular expansion and expression levels within genes. Although blood DNA methylation has been studied in humans and other species, its prominence in cattle is largely unknown. This study aimed to methodically probe the genomic methylation map of Xinjiang brown (XJB) cattle suffering from bovine respiratory disease (BRD), consequently widening cattle blood methylome ranges. Genome-wide DNA methylation profiling of the XJB blood was investigated through whole-genome bisulfite sequencing (WGBS). Many differentially methylated regions (DMRs) obtained by comparing the cases and controls groups were found within the CG, CHG, and CHH (where H is A, T, or C) sequences (16,765, 7502, and 2656, respectively), encompassing 4334 differentially methylated genes (DMGs). Furthermore, GO/KEGG analyses showed that some DMGs were involved within immune response pathways. Combining WGBS-Seq data and existing RNA-Seq data, we identified 71 significantly differentially methylated (DMGs) and expressed (DEGs) genes (p < 0.05). Next, complementary analyses identified nine DMGs (LTA, STAT3, IKBKG, IRAK1, NOD2, TLR2, TNFRSF1A, and IKBKB) that might be involved in the immune response of XJB cattle infected with respiratory diseases. Although further investigations are needed to confirm their exact implication in the involved immune processes, these genes could potentially be used for a marker-assisted selection of animals resistant to BRD. This study also provides new knowledge regarding epigenetic control for the bovine respiratory immune process.

BackgroundDNA methylation is an epigenetic regulatory form that plays an important role in regulating the gene expression and the tissues development.. However, DNA methylation regulators involved in sheep muscle development remain unclear. To explore the functional importance of genome-scale DNA methylation during sheep muscle growth, this study systematically investigated the genome-wide DNA methylation profiles at key stages of Hu sheep developmental (fetus and adult) using deep whole-genome bisulfite sequencing (WGBS).ResultsOur study found that the expression levels of DNA methyltransferase (DNMT)-related genes were lower in fetal muscle than in the muscle of adults. The methylation levels in the CG context were higher than those in the CHG and CHH contexts, and methylation levels were highest in introns, followed by exons and downstream regions. Subsequently, we identified 48,491, 17, and 135 differentially methylated regions (DMRs) in the CG, CHG, and CHH sequence contexts and 11,522 differentially methylated genes (DMGs). The results of bisulfite sequencing PCR (BSP) correlated well with the WGBS-Seq data. Moreover, Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) functional annotation analysis revealed that some DMGs were involved in regulating skeletal muscle development and fatty acid metabolism. By combining the WGBS-Seq and previous RNA-Seq data, a total of 159 overlap genes were obtained between differentially expressed genes (DEGs) and DMGs (FPKM > 10 and fold change > 4). Finally, we found that 9 DMGs were likely to be involved in muscle growth and metabolism of Hu sheep.ConclusionsWe systemically studied the global DNA methylation patterns of fetal and adult muscle development in Hu sheep, which provided new insights into a better understanding of the epigenetic regulation of sheep muscle development.

DNA methylation has frequently been implicated in sex determination and differentiation in teleost species. In order to detect the DNA methylation patterns established during sexual differentiation in tiger pufferfish T. rubripes, we performed comprehensive whole genome methylation sequencing and analyses of the gonads of male, female, and pseudo male. We obtained a total of 33.12, 32.44, and 31.60 Gb clean data for male, female, and pseudo male, with a sequencing depth of 66.44×, 60.47× and 54.86×, respectively. The methylation level of cytosine (C) residues in the genomic DNA from gonads was 11.016%, 10.428%, and 11.083% in male, female, and pseudo male, respectively. More than 65% of C methylation was at CpG sites, and less than 1% was at CHG and CHH sites. In each regulatory element, there were low methylation levels on both sides of the transcription start site, and higher methylation levels in exons, introns, and downstream of genes. The highest mCpG was on chromosome 8 and the lowest mCpG was on chromosome 5. Comparisons of whole-genome DNA methylation between pairs of samples revealed that there were 3,173 differentially methylated regions (DMRs) between female and male, and 3,037 DMRs between male and pseudo male, corresponding to 0.232% and 0.223% of the length of the genome, respectively. There were only 1,635 DMRs between female and pseudo male, representing 0.127% of the length of the genome. A number of differentially methylated genes (DMGs) related to sex determination and differentiation were selected, such as amhr2 and pfcyp19a. After Bisulfite Sequencing PCR (BSP) verification, amhr2 was exhibited low methylation level in normal males and pseudo male, and high methylation level in normal females but pfcyp19a showed low methylation level in normal females and high methylation level in normal males and pseudo males. These results provide information about the molecular epigenetic mechanisms of DNA methylation during low-temperature induced masculinization of tiger pufferfish, and increase our understanding of the mechanisms of sex determination and differentiation in this important aquaculture fish species.