Differentially Methylated Genes Research Articles

Genome-wide DNA methylation and transcription analysis reveal the potential epigenetic mechanism of heat stress response in the sea cucumber Apostichopus japonicus

DNA methylation is an important epigenetic modification that regulates many biological processes. The sea cucumber Apostichopus japonicus often suffers from heat stress that affects its growth and leads to significant economic losses. In this study, the mRNA expression patterns and DNA methylation characteristics in the body wall of A. japonicus under heat stress were analyzed by whole-genome bisulfite sequencing (WGBS) and transcriptome sequencing (RNA-seq). We found that CpG was the main DNA methylation type, and heat stress caused a significant increase in the overall methylation level and methylation rate, especially in the intergenic region of the A. japonicus genome. In total, 1,409 differentially expressed genes (DEGs) and 17,927 differentially methylated genes (DMGs) were obtained by RNA-seq and WGBS, respectively. Association analysis between DNA methylation and transcription identified 569 negatively correlated genes in both DMGs and DEGs, which indicated that DNA methylation affects on transcriptional regulation in response to heat stress. These negatively correlated genes were significantly enriched in pathways related to energy metabolism and immunoregulation, such as the thyroid hormone signaling pathway, renin secretion, notch signaling pathway and microRNAs in cancer. In addition, potential key genes, including heat shock protein (hsp70), calcium-activated chloride channel regulator 1(clca1), and tenascin R (tnr), were obtained and their expression and methylation were preliminarily verified. The results provide a new perspective for epigenetic and transcriptomic studies of A. japonicus response to heat stress, and provide a reference for breeding sea cucumbers resistant to high temperatures.

Sex-specific transcription and DNA methylation landscapes of the Asian citrus psyllid, a vector of huanglongbing pathogens.

The relationship of DNA methylation and sex-biased gene expression is of high interest, it allows research into mechanisms of sexual dimorphism and the development of potential novel strategies for insect pest control. The Asian citrus psyllid, Diaphorina citri Kuwayama, is a major vector for the causative agents of Huanglongbing (HLB), which presents an unparalleled challenge to citrus production worldwide. Here, we identify the X chromosome of D. citri and investigate differences in the transcription and DNA methylation landscapes between adult virgin males and females. We find a large number of male-biased genes on the autosomes and a depletion of such on the X chromosome. We have also characterised the methylome of D. citri, finding low genome-wide levels, which is unusual for an hemipteran species, as well as evidence for both promoter and TE methylation. Overall, DNA methylation profiles are similar between the sexes but with a small number of differentially methylated genes found to be involved in sex differentiation. There also appears to be no direct relationship between differential DNA methylation and differential gene expression. Our findings lay the groundwork for the development of novel epigenetic-based pest control methods, and given the similarity of theD. citri methylome to some other insect species, these methods could be applicable across agricultural insect pests.

Whole genome methylation combined with RNA-seq reveals the protective effects of Gualou-Xiebai herb pair in foam cells through DNA methylation mediated PI3K-AKT signaling pathway.

DNA methylation, including aberrant hypomethylation and hypermethylation, plays a significant role in atherosclerosis (AS); therefore, targeting the unbalanced methylation in AS is a potential treatment strategy. Gualou-xiebai herb pair (GXHP), a classic herb combination, have been used for the treatment of atherosclerotic-associated diseases in traditional Chinese medicine. However, the effects and underlying mechanism of GXHP on AS remain nebulous. In this study, the CCK-8 method was applied to determine the non-toxic treatment concentrations for GXHP. The formation of foam cells played a critical role in AS, so the foam cells model was established after RAW264.7 cells were treated with ox-LDL. The contents of total cholesterol (TC) and free cholesterol (FC) were determined by Gas chromatography-mass spectrometry (GC-MS). Enzyme-linked immunosorbent assay (ELISA) was used to check the expressions of inflammatory factors including IL-1β, TNF-α, and VCAM-1. Methyl-capture sequencing (MC-seq) and RNA-seq were applied to observe the changes in genome-wide DNA methylation and gene expression, respectively. Kyoto Encyclopedia of Genes and Genomes (KEGG) were performed to analyze differentially methylated genes (DMGs) and differentially expressed genes (DEGs). The targeted signaling pathway was selected and verified using western blotting (WB). The results showed that the lipids and inflammatory factors in foam cells significantly increased. GXHP significantly reduced the expression of TC, FC, and inflammatory factors. MC-seq and RNA-seq showed that GXHP not only corrected the aberrant DNA hypermethylation, but also DNA hypomethylation, thus restored the aberrant DEGs in foam cells induced by ox-LDL. GXHP treatment may target the PI3K-Akt signaling pathway. GXHP reduced the protein levels of phosphorylated(p)-PI3K and p-AKT in foam cells. Our data suggest that treatment with GXHP showed protective effects against AS through the inhibition of DNA methylation mediated PI3K-AKT signaling pathway, suggesting GXHP as a novel methylation-based agent.

Transgenerational effects of intertidal environment on physiological phenotypes and DNA methylation in Pacific oysters

Climate change and intensifying human activity are posing serious threats to marine organisms. The fluctuating intertidal zone forms a miniature ecosystem of a rapidly changing environment for studying biological adaptation. Transgenerational plasticity (TGP), an evolutionary phenomenon in which parental experience influences offspring phenotypes, provides an avenue for adaptation, but the molecular mechanism was poorly understood in marine molluscs. In this study, wild Pacific oysters (Crassostrea gigas), which were collected from intertidal zones, were used to conduct two-generation breeding in a subtidal area combined with a heat shock experiment in the laboratory to investigate the intertidal environment-induced TGP under temperate subtidal condition and thermally exposed condition, respectively. We showed that TGP could influence the physiological phenotypes related to the status of oxidation and energy in non-stress-exposed subtidal offspring for at least two generations. Genomic DNA methylation exhibited heritable divergence between intertidal and subtidal oysters, and 1655 (or 42.83 %) differentially methylated genes (DMGs) in F0 were continuously reserved to F2, which may mediate physiological TGP by participating in biological processes including macromolecule metabolism, cellular responses to stress, and the positive regulation of molecular function, especially fatty acid metabolism. The intertidal experience also influenced the thermal plasticity of physiological phenotypes within and across generations. Totally, 320 (or 14.74 %) specific thermal response DMGs in the intertidal F0 generation were identified in F1 and F2, participating in pathways including carbohydrate, lipid, and energy metabolism, signal transduction, and the organismal immune system, which suggested transgenerational intertidal effect mediated by these genes could positively contribute to stress adaptation and had potential applications for aquaculture. This study demonstrates an epigenetic mechanism for TGP in stress adaptation in marine molluscs, and provides new avenues to improve the stress adaptation for marine resource conservation and aquaculture.

Transcriptome and methylome dynamics in the gills of large yellow croaker (Larimichthys crocea) during low-salinity adaption

DNA methylation is a critical epigenetic modification that dynamically regulates gene expression in organisms facing abiotic stress. However, few studies have comprehensively examined the role of DNA methylation in marine fish during environmental adaptation. Therefore, this study explored the methylome dynamics and DNA methylation regulation mechanisms in large yellow croaker (Larimichthys crocea) during low-salinity adaption. The methylation level in the gills was notably raised in the S-group (5‰ salinity) compared to C-group (25‰ salinity). A total of 109 differentially methylated promoter target genes and 581 differentially expressed genes were identified via whole-genome bisulfite sequencing (WGBS) and RNA-seq of gills in the two salinity groups, respectively. Moreover, 23 hypo-methylated/up-regulated differentially methylated genes (DMGs) and 28 hyper-methylated/down-regulated DMGs were identified through integrative analysis, which were mainly enriched in signal transduction, ion exchange, energy metabolism, and cytoskeleton system and other biological processes. Collectively, our findings suggested that low-salinity stress can induce adaptive genome-wide DNA methylation changes, which can in turn affect the transcription of genes in large yellow croaker during low-salinity adaptation. Therefore, our findings provide new insights into the regulatory mechanisms of marine fish in response to rapid environmental changes.

Research on the pathological mechanism of rectal adenocarcinoma based on DNA methylation.

Colorectal cancer is one of the 3 most common cancers worldwide. In this study, a weighted network-based analysis method was proposed to explore the pathological mechanisms and prognostic targets of rectal adenocarcinoma (READ) at the deoxyribonucleic acid (DNA) methylation level. In this study, we downloaded clinical information and DNA methylation data from The Cancer Genome Atlas database. Differentially methylated gene analysis was used to identify the differential methylated genes in READ. Canonical correlation analysis was used to construct the weighted gene regulatory network for READ. Multilevel analysis and association analyses between gene modules and clinical information were used to mine key modules related to tumor metastasis evaluation. Genetic significance analysis was used to identify methylation sites in key modules. Finally, the importance of these methylation sites was confirmed using survival analysis. DNA methylation datasets from 90 cancer tissue samples and 6 paracancerous tissue samples were selected. A weighted gene regulatory network was constructed, and a multilevel algorithm was used to divide the gene co-expression network into 20 modules. From gene ontology enrichment analysis, characteristic M was related to biological processes such as the chemotaxis of fibroblast growth factors and the activation and regulation of immune cells etc and characteristic N was associated with the regulation of cytoskeleton formation, mainly microtubules and flagella, regulation of synapses, and regulation of cell mitosis. Based on the results of survival analysis, 7 key methylation sites were found closely correlated to the survival rate of READ, such as cg04441191 (microtubule-associated protein 4 [MAP4]), cg05658717 (KSR2), cg09622330 (GRIN2A), cg10698404 (YWHAG), cg17047993 (SPAG9), cg24504843 (CEP135), and cg24531267 (CEP250). Mutational and transcriptomic level studies revealed significant differences in DNA methylation, single nucleotide polymorphism, and transcript levels between YWHAG and MAP4 in normal tissues compared to tumor tissues, and differential expression of the 2 proteins in immunohistochemistry. Therefore, potential targeting drugs were screened for these 2 proteins for molecular docking, and artenimol was found to bind to MAP4 protein and 27-hydroxycholesterol to YWHAG. Our study found that key methylation sites played an important role in tumor metastasis and were associated with the prognosis of READ. Mutations and methylation may jointly regulate the transcription and translation of related genes, which in turn affect cancer progression. This may provide some new potential therapeutic targets and thoughts for the prognosis of READ.

Non-Association of Driver Alterations in PTEN with Differential Gene Expression and Gene Methylation in IDH1 Wildtype Glioblastomas

During oncogenesis, alterations in driver genes called driver alterations (DAs) modulate the transcriptome, methylome and proteome through oncogenic signaling pathways. These modulatory effects of any DA may be analyzed by examining differentially expressed mRNAs (DEMs), differentially methylated genes (DMGs) and differentially expressed proteins (DEPs) between tumor samples with and without that DA. We aimed to analyze these modulations with 12 common driver genes in Isocitrate Dehydrogenase 1 wildtype glioblastomas (IDH1-W-GBs). Using Cbioportal, groups of tumor samples with and without DAs in these 12 genes were generated from the IDH1-W-GBs available from "The Cancer Genomics Atlas Firehose Legacy Study Group" (TCGA-FL-SG) on Glioblastomas (GBs). For all 12 genes, samples with and without DAs were compared for DEMs, DMGs and DEPs. We found that DAs in PTEN were unassociated with any DEM or DMG in contrast to DAs in all other drivers, which were associated with several DEMs and DMGs. This contrasting PTEN-related property of being unassociated with differential gene expression or methylation in IDH1-W-GBs was unaffected by concurrent DAs in other common drivers or by the types of DAs affecting PTEN. From the lists of DEMs and DMGs associated with some common drivers other than PTEN, enriched gene ontology terms and insights into the co-regulatory effects of these drivers on the transcriptome were obtained. The findings from this study can improve our understanding of the molecular mechanisms underlying gliomagenesis with potential therapeutic benefits.

Genome-wide DNA methylation analysis of discordant monozygotic twins reveals consistent sites of differential methylation associated with congenital heart disease

BackgroundDespite being essentially genetically identical, monozygotic (MZ) twins can be discordant for congenital heart disease (CHD), thus highlighting the importance of in utero environmental factors for CHD pathogenesis. This study aimed to identify the epigenetic variations between discordant MZ twin pairs that are associated with CHD at birth. MethodsCord blood of CHD-discordant MZ twins from the Chongqing Longitudinal Twin Study Cohort was subjected to whole-genome bisulfite sequencing, then validated by MeDIP-qPCR and qRT-PCR. Results379 DMRs mapped to 175 differentially methylated genes (DMGs) were associated with CHD. Functional enrichment analysis identified these DMGs are involved in histone methylation, actin cytoskeleton organization, the regulation of cell differentiation, and adrenergic signaling in cardiomyocytes. Of note, SPESP1 and NOX5 were hypermethylated in CHD, and associated with lower gene expression levels. ConclusionsSpecific DNA methy (DNAm) variations in cord blood were associated with CHD, thus illustrating new biomarkers and potential interventional targets for CHD. Trial registrationChiCTR-OOC-16008203, registered on 1 April 2016 at the Chinese Clinical Trial Registry.

DNA methylome and transcriptome identified Key genes and pathways involved in Speckled Eggshell formation in aged laying hens

BackgroundThe quality of poultry eggshells is closely related to the profitability of egg production. Eggshell speckles reflect an important quality trait that influences egg appearance and customer preference. However, the mechanism of speckle formation remains poorly understood. In this study, we systematically compared serum immune and antioxidant indices of hens laying speckled and normal eggs. Transcriptome and methylome analyses were used to elucidate the mechanism of eggshell speckle formation.ResultsThe results showed that seven differentially expressed genes (DEGs) were identified between the normal and speckle groups. Gene set enrichment analysis (GSEA) revealed that the expressed genes were mainly enriched in the calcium signaling pathway, focal adhesion, and MAPK signaling pathway. Additionally, 282 differentially methylated genes (DMGs) were detected, of which 15 genes were associated with aging, including ARNTL, CAV1, and GCLC. Pathway analysis showed that the DMGs were associated with T cell-mediated immunity, response to oxidative stress, and cellular response to DNA damage stimulus. Integrative analysis of transcriptome and DNA methylation data identified BFSP2 as the only overlapping gene, which was expressed at low levels and hypomethylated in the speckle group.ConclusionsOverall, these results indicate that aging- and immune-related genes and pathways play a crucial role in the formation of speckled eggshells, providing useful information for improving eggshell quality.

Impact of FecB Mutation on Ovarian DNA Methylome in Small-Tail Han Sheep.

In this study, adult Han ewes with homozygous FecB(B)/FecB(B) mutations (Han BB group) and ewes with FecB(+)/FecB(+) wildtype (Han ++ group) were selected. Methylated DNA immunoprecipitation and high-throughput sequencing (MeDIP-seq) was used to identify differences in methylated genes in ovary tissue. We examined differences in DNA methylation patterns between HanBB and Han ++ sheep. In both sheep, methylated reads were mainly distributed at the gene body regions, CpG islands and introns. The differentially methylated genes were enriched in neurotrophy in signaling pathway, Gonadotropin Releasing Hormone (GnRH) signaling pathway, Wnt signaling pathway, oocyte meiosis, vascular endothelial growth factor (VEGF) signaling pathway, etc. Differentially-methylated genes were co-analyzed with differentially-expressed mRNAs. Several genes which could be associated with female reproduction were identified, such as FOXP3 (forkhead box P3), TMEFF2 (Transmembrane Protein with EGF Like and Two Follistatin Like Domains 2) and ADAT2 (Adenosine Deaminase TRNA Specific 2). We constructed a MeDIP-seq based methylomic study to investigate the ovarian DNA methylation differences between Small-Tail Han sheep with homozygous FecB mutant and wildtype, and successfully identified FecB gene-associated differentially-methylated genes. This study has provided information with which to understand the mechanisms of FecB gene-induced hyperprolificacy in sheep.

DNA methylation regulatory patterns and underlying pathways behind the co-pathogenesis of allergic rhinitis and chronic spontaneous urticaria.

Allergic rhinitis (AR) and chronic spontaneous urticaria (CSU) are often concurrent in patients. Changes in DNA methylation affect T cell biological processes, which may explain the occurrence and progression of comorbidity. However, downstream regulatory pathways of DNA methylation in two diseases and the underlying mechanisms have not been fully elucidated. The GSE50101, GSE72541, GSE50222 and OEP002482 were mined for the identification of differentially expressed genes (DEGs) or co-expressed genes and differentially methylated genes (DMGs) in AR and CSU patients. We applied GO analysis and consensus clustering to study the potential functions and signal pathways of selected genes in two diseases. GSVA and logistic regression analysis were used to find the regulatory pathway between DNA methylation and activation patterns of CD4+ T cells. Besides, we used the Illumina 850k chip to detect DNA methylation expression profiles and recognize the differentially methylated CpG positions (DMPs) on corresponding genes. Finally, we annotated the biological process of these genes using GO and KEGG pathway analysis. The AR-related DEGs were found closely related to the differentiation and activation of CD4+ T cells. The DEGs or co-expressed genes of CD4+ T cells in AR and CSU patients were also clustered using GO and KEGG analysis and we got 57 co-regulatory pathways. Furthermore, logistic regression analysis showed that the regulation of cellular component size was closely related to the activation of CD4+ T cells regulated by DNA methylation. We got self-tested data using the Illumina 850k chip and identified 98 CpGs that were differentially methylated in patients. Finally, we mapped the DMPs to 15 genes and found that they were mainly enriched in the same CD4+T cell regulating pathway. Our study indicated that DNA methylation affected by pollen participated in the activation patterns of CD4 + T cells, providing a novel direction for the symptomatic treatment of the co-occurrence of AR and CSU.

Association of environmental exposure to chromium with differential DNA methylation: An epigenome-wide study.

Introduction: Previous studies have reported that chromium (Cr)-induced epigenetic alterations and DNA methylation play a vital role in the pathogenesis of diseases induced by chromium exposure. Epigenomic analyses have been limited and mainly focused on occupational chromium exposure; their findings are not generalizable to populations with environmental Cr exposure. Methods: We identified the differential methylation of genes and regions to elucidate the mechanisms of toxicity related to environmental chromium exposure. DNA methylation was measured in blood samples collected from individuals in Cr-contaminated (n = 10) and unexposed areas (n = 10) by using the Illumina Infinium HumanMethylation850K array. To evaluate the relationship between chromium levels in urine and CpG methylation at 850 thousand sites, we investigated differentially methylated positions (DMPs) and differentially methylated regions (DMRs) by using linear models and DMRcate method, respectively. The model was adjusted for biologically relevant variables and estimated cell-type compositions. Results: At the epigenome-wide level, we identified five CpGs [cg20690919 (p FDR =0.006), cg00704664 (p FDR =0.024), cg10809143 (p FDR =0.043), cg27057652 (p FDR =0.047), cg05390480 (p FDR =0.024)] and one DMR (chr17: 19,648,718-19,648,972), annotated to ALDH3A1 genes (p < 0.05) as being significantly associated with log2 transformed urinary chromium levels. Discussion: Environmental chromium exposure is associated with DNA methylation, and the significant DMPs and DMR being annotated to cause DNA damage and genomic instability were found in this work. Research involving larger samples is required to further explore the epigenetic effect of environmental chromium exposure on health outcomes through DNA methylation.

A gene signature driven by abnormally methylated DEGs was developed for TP53 wild-type ovarian cancer samples by integrative omics analysis of DNA methylation and gene expression data.

Integrated omics analysis based on transcriptome and DNA methylation data combined with machine learning methods is very promising for the diagnosis, prognosis, and classification of cancer. In this study, the DNA methylation and gene expression data of ovarian cancer (OC) were analyzed to identify abnormally methylated differentially expressed genes (DEGs), screen potential therapeutic agents for OC, and construct a risk model based on the abnormally methylated DEGs to predict patient prognosis. The gene expression and DNA methylation data of primary OC samples with tumor protein 53 (TP53) wild-type and normal samples were obtained from The Cancer Genome Atlas (TCGA) database. DEGs with aberrant methylation were analyzed by screening the intersection between DEGs and differentially methylated genes (DMGs). We attempted to search for potential drugs targeting DEGs with aberrant methylation by employing a network medicine framework. A gene signature based on the DEGs with aberrant methylation was constructed by regularized least absolute shrinkage and selection operator (LASSO) regression analysis. A total of 440 aberrant methylated DEGs were screened. Based on their gene expression profiles and methylation data from different regions, the results of both discriminative pattern recognition analysis and principal component analysis (PCA) showed a significant separation between tumor tissue and healthy ovarian tissue. In total, 126 potential therapeutic drugs were identified for OC by network-based proximity analysis. Five genes were identified in 440 aberrant methylated DEGs, which formed an aberrant methylated DEGs-driven gene signature. This signature could significantly distinguish the different overall survivals (OS) of OC patients and showed better predictive performance in both the training and validation sets. In this study, the DNA methylation and gene expression data of OC were analyzed to identify abnormally methylated DEGs and potential therapeutic drugs, and a gene signature based on five aberrant methylation DEGs was constructed, which could better predict the risk of death in patients.

Abstract A078: Exploration, analysis and visualization of TCGA pancancer data in the cBioPortal: Prostate adenocarcinoma genomics

Abstract Background Prostate cancer cases are more likely among Black men than their counterparts – and only a few epigenetic studies of prostate cancer have been performed. We hypothesized that differential gene methylation may contribute to racial disparity of prostate adenocarcinoma. MethodsFrom the Prostate Adenocarcinoma in the TCGA PanCancer data, we selected 154 samples from 7 Blacks and 147 Whites. We analyzed the dataset by race/ethnicity, gene methylation and genomic alterations. We test our hypothesis by analyzing genetic alterations and methylation; genes with highest frequencies by race/ethnic group and significant methylation were recorded, respectively. Log ratio of mean methylation in Blacks to mean in whites were performed across the data. Multivariable data analysis, to test for association between race and genes, was performed using the Log Ratio statistic. Kaplan-Meier survival analysis was also performed on the sample data. Results A total of 494 prostate cancer patients were identified from the cBioPortal; a sample cohort of 154 patients comprising 7 Blacks and 147 Whites were used to for this analysis. The median age of the selected men was 63 years (range: 42-78). Using log ratio statistic, the data demonstrated significant (q &amp;lt; 0.05) by race/ethnicity. The association between race/ethnicity and methylation for the 60 genes (14 of which were Black and 46 were white) investigated in prostate cancer data also varied differentially. In Blacks, for example, UBL5 (m: 0.62 vs 0.53, q &amp;lt;0.05), ART4 (m: 0.95 vs 0.90, q &amp;lt;0.05), …, KATNB1 (m: 0.90 vs 0.82, q &amp;lt;0.05) were differentially methylated (Table 1) compared to Whites. Similarly, among Whites, ZNF57 (m: 0.30 vs 0.03, q &amp;lt;0.005), LRIG1 (m: 0.23 vs 0.07, q &amp;lt;0.005), …, EN2 (m: 0.08 vs 0.04, q &amp;lt;0.005) were also differentially methylated in comparison with Blacks. In the genomic alteration analysis: Blacks have more genes with frequent (&amp;gt;30%) alterations than Whites.10 genes had high percentage of gene alterations among Blacks (particularly CSMD1 and KHDRBS3) overlapping with 2 in whites (ERG and TMPRSS2). In this data, age groups 44-56 and 56-61 are the most frequent. On the survival analysis, NUDT12 (LR: 0.05, q &amp;lt;0.005), RPL26 (LR: 0.02 q &amp;lt;0.005), and others were associated with survival. Discussion Comparing the race/ethnicity groups using the methylation and the genetic alterations data, we were able to confirm that genetic alteration and DNA methylation are significantly differentiated, respectively. There exists significant association between race category, age category and genes, causing variation in methylation and genomic alterations. Conclusion In the prostate cancer data, race/ethnicity has an impact on genomic methylation. This factor also has an impact on genetic alterations. Prostate cancer epigenetic studies must take race/ethnicity into consideration. Our genomic samples were sufficient to show a difference in the epigenetic landscape between Blacks and whites. Citation Format: Olumide Arigbede, Sarah G. Buxbaum, Sara Falzarano, Suhn Rhie. Exploration, analysis and visualization of TCGA pancancer data in the cBioPortal: Prostate adenocarcinoma genomics [abstract]. In: Proceedings of the 15th AACR Conference on the Science of Cancer Health Disparities in Racial/Ethnic Minorities and the Medically Underserved; 2022 Sep 16-19; Philadelphia, PA. Philadelphia (PA): AACR; Cancer Epidemiol Biomarkers Prev 2022;31(1 Suppl):Abstract nr A078.

Characterization of lung adenocarcinoma based on immunophenotyping and constructing an immune scoring model to predict prognosis.

Background: Lung cancer poses great threat to human health, and lung adenocarcinoma (LUAD) is the main subtype. Immunotherapy has become first line therapy for LUAD. However, the pathogenic mechanism of LUAD is still unclear. Methods: We scored immune-related pathways in LUAD patients using single sample gene set enrichment analysis (ssGSEA) algorithm, and further identified distinct immune-related subtypes through consistent clustering analysis. Next, immune signatures, Kaplan-Meier survival analysis, copy number variation (CNV) analysis, gene methylation analysis, mutational analysis were used to reveal differences between subtypes. pRRophetic method was used to predict the response to chemotherapeutic drugs (half maximal inhibitory concentration). Then, weighted gene co-expression network analysis (WGCNA) was performed to screen hub genes. Significantly, we built an immune score (IMscore) model to predict prognosis of LUAD. Results: Consensus clustering analysis identified three LUAD subtypes, namely immune-Enrich subtype (Immune-E), stromal-Enrich subtype (Stromal-E) and immune-Deprived subtype (Immune-D). Stromal-E subtype had a better prognosis, as shown by Kaplan-Meier survival analysis. Higher tumor purity and lower immune cell scores were found in the Immune-D subtype. CNV analysis showed that homologous recombination deficiency was lower in Stromal-E and higher in Immune-D. Likewise, mutational analysis found that the Stromal-E subtype had a lower mutation frequency in TP53 mutations. Difference in gene methylation (ZEB2, TWIST1, CDH2, CDH1 and CLDN1) among three subtypes was also observed. Moreover, Immune-E was more sensitive to traditional chemotherapy drugs Cisplatin, Sunitinib, Crizotinib, Dasatinib, Bortezomib, and Midostaurin in both the TCGA and GSE cohorts. Furthermore, a 6-gene signature was constructed to predicting prognosis, which performed better than TIDE score. The performance of IMscore model was successfully validated in three independent datasets and pan-cancer.

Maternal SARS-CoV-2 exposure alters infant DNA methylation

BackgroundInfection during pregnancy can increase the risk of neurodevelopmental disorders in offspring. The impact of maternal SARS-CoV-2 infection on infant neurodevelopment is poorly understood. The maternal immune response to infection may be mimicked in rodent models of maternal immune activation which recapitulate altered neurodevelopment and behavioural disturbances in the offspring. In these models, epigenetic mechanisms, in particular DNA methylation, are one pathway through which this risk is conferred in utero to offspring. We hypothesised that in utero exposure to SARS-CoV-2 in humans may alter infant DNA methylation, particularly in genes associated with neurodevelopment. We aimed to test this hypothesis in a pilot sample of children in Victoria, Australia, who were exposed in utero to SARS-CoV-2. MethodsDNA was extracted from buccal swab specimens from (n = 4) SARS-CoV-2 in utero exposed and (n = 4) non-exposed infants and methylation status assessed across 850,000 methylation sites using an Illumina EPIC BeadChip. We also conducted an exploratory enrichment analysis using Gene Ontology annotations. Results1962 hypermethylated CpG sites were identified with an unadjusted p-value of 0.05, where 1133 CpGs mapped to 959 unique protein coding genes, and 716 hypomethylated CpG sites mapped to 559 unique protein coding genes in SARS-CoV-2 exposed infants compared to non-exposed. One differentially methylated position (cg06758191), located in the gene body of AFAP1 that was hypomethylated in the SARS-CoV-2 exposed cohort was significant after correction for multiple testing (FDR-adjusted p-value <0.00083). Two significant differentially methylated regions were identified; a hypomethylated intergenic region located in chromosome 6p proximal to the genes ZP57 and HLA-F (fwer <0.004), and a hypomethylated region in the promoter and body of the gene GAREM2 (fwer <0.036). Gene network enrichment analysis revealed differential methylation in genes corresponding to pathways relevant to neurodevelopment, including the ERBB pathway. ConclusionThese pilot data suggest that exposure to SARS-CoV-2 in utero differentially alters methylation of genes in pathways that play a role in human neurodevelopment.

Identification of Signature Genes in the PD-1 Relative Gastric Cancer Using a Combined Analysis of Gene Expression and Methylation Data.

The morbidity and mortality rates for gastric cancer (GC) rank second among all cancers, indicating the serious threat it poses to human health, as well as human life. This study aims to identify the pathways and genes as well as investigate the molecular mechanisms of tumor-related genes in gastric cancer (GC). We compared differentially expressed genes (DEGs) and differentially methylated genes (DMGs) in gastric cancer and normal tissue samples using The Cancer Genome Atlas (TCGA) data. The Kyoto Encyclopedia of Gene and Genome (KEGG) and the Gene Ontology (GO) enrichment analysis' pathway annotations were conducted on DMGs and DEGs using a clusterProfiler R package to identify the important functions, as well as the biological processes and pathways involved. The intersection of the two was chosen and defined as differentially methylated and expressed genes (DMEGs). For DMEGs, we used the principal component analysis (PCA) to differentiate gastric cancer from adjacent samples. The linear discriminant analysis method was applied to categorize the samples using DMEGs methylation data and DMEGs expression profiles data and was validated using the leave-one-out cross-validation (LOOCV) method. We plotted the ROC curve for the classification and calculated the AUC (area under the ROC curve) value for a more intuitive view of the classification effect. We also used the NetworkAnalyst 3.0 tool to analyze DMEGs, using DrugBank to acquire information on protein-drug interactions and generate a network map of gene-drug interactions. We identified a total of 971 DMGs in 188 PD-1 negative and 187 PD-1 positive gastric cancer samples obtained from TCGA. The KEGG and GO enrichment analysis showed the involvement of the regulation of ion transmembrane transport, collagen-containing extracellular matrix, cell-cell junction, and peptidase regulator activity. We simultaneously obtained 1,189 DEGs, out of which 986 were downregulated, while 203 were upregulated in tumors. The enriched analysis of the GO's and KEGG's pathways indicated that the most significant pathways included an intestinal immune network for IgA production, Staphylococcus aureus infection, cytokine-cytokine receptor interaction, and viral protein interaction with cytokine and cytokine receptor, which have previously been linked with gastric cancer. The compound DB01830 can bind well to the active site of the LCK protein and shows good stability, thus making it a potential inhibitor of the LCK protein. To observe the relationship between DMEGs' expression and prognosis, we observed 10 genes, among which were TRIM29, TSPAN8, EOMES, PPP1R16B, SELL, PCED1B, IYD, JPH1, CEACAM5, and RP11-44K6.2. Their high expressions were related to high risks. Besides, those genes were validated in different internal and external validation sets. These results may provide potential molecular biological therapy for PD-1 negative gastric cancer.

CfDNA Methylation Profiles and T-Cell Differentiation in Women with Endometrial Polyps.

DNA methylation is a part of the regulatory mechanisms of gene expression, including chromatin remodeling and the activity of microRNAs, which are involved in the regulation of T-cell differentiation and function. However, the role of cfDNA methylation in T-cell differentiation is entirely unknown. In patients with endometrial polyps (EPs), we have found an imbalance of T-cell differentiation and an aberrant cfDNA methylation profile, respectively. In this study, we investigated the relationship between cfDNA methylation profiles and T-cell differentiation in 14 people with EPs and 27 healthy controls. We found that several differentially methylated genes (DMGs) were associated with T-cell differentiation in people with EPs (ITGA2-Naïve CD4, r = -0.560, p = 0.037; CST9-EMRA CD4, r = -0.626, p = 0.017; and ZIM2-CM CD8, r = 0.576, p = 0.031), but not in healthy controls (all p &gt; 0.05). When we combined the patients' characteristics, we found a significant association between ITGA2 methylation and polyp diameter (r = 0.562, p = 0.036), but this effect was lost when adjusting the level of Naïve CD4 T-cells (r = 0.038, p = 0.903). Moreover, the circulating sex hormone levels were associated with T-cell differentiation (estradiol-Naïve CD4, r = -0.589, p = 0.027), and the cfDNA methylation profile (testosterone-ZIM2, r = -0.656, p = 0.011). In conclusion, this study has established a link between cfDNA methylation profiles and T-cell differentiation among people with EPs, which may contribute to the etiology of EPs. Further functional studies are warranted.

Genome-wide DNA methylation profile of human dental pulp stem cells during odontogenic differentiation

ObjectivesHuman dental pulp stem cells (hDPSCs) is essential for dentin formation and regeneration, emerging evidence revealed that epigenetic regulation plays vital roles in odontogenic differentiation of hDPSCs. The purpose of this study was to explore the genome-wide DNA methylation changes during odontogenic differentiation of hDPSCs. DesignhDPSCs were isolated from young healthy premolars and reduced representation bisulfite sequencing (RRBS) was taken to detect the genome-wide DNA methylation profile of hDPSCs during odontogenic differentiation. Genome-wide DNA methylation map, differentially methylated region (DMR) analysis, Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis were performed. ResultsWe found a totally different DNA methylation patterns during the odontogenic differentiation of hDPSCs. A total of 9309 differentially methylated genes (DMG) were identified. Bio-information analysis revealed that calcium signaling pathway, pathways in cancer, and HTLV–I infection signaling pathways may play potential roles in odontogenic differentiation of hDPSCs. NOTCH1, WNT7B, and AXIN2 proteins were related with calcium signaling pathway. ConclusionsThis study revealed a comprehensive analysis of global DNA methylation profiles in odontogeinc differentiation of hDPSCs and provided several possible underlying signaling pathways and candidate genes that may regulate the odontogenic differentiation of hDPSCs.

Comparison of DNA methylation profiles of human embryos cultured in either uninterrupted or interrupted incubators.

We aimed to compare the DNA methylation profiles of human embryos cultured in uninterrupted or interrupted incubators. This study included 9 women, ≤ 30years old (range: 20-30years), without a history of genetic diseases or smoking, undergoing ICSI treatment, and each woman donated one oocyte. Embryos were randomly assigned to culture in either time-lapse imaging or standard incubators after ICSI. We compared the DNA methylation profiles of human eight-cell embryos cultured in uninterrupted condition using time-lapse imaging (TLI) incubator (EmbryoScope) to those cultured in interrupted culture model using standard incubators (SI, G185 K-System). Nine single-cell whole-genome bisulfite sequencing (WGBS) datasets were analyzed, including four SI-cultured embryos and five TLI-cultured embryos at the eight-cell stage. A total of 581,140,020 and 732,348,182 clean reads were generated from the TLI and SI groups, respectively. TLI-cultured embryos had similar genome-wide methylation patterns to SI-cultured embryos. There were no significant differences in the methylation and transcription levels of transposable elements and imprinted control regions. Although a total of 198 differentially methylated genes (DMGs) were identified, only five DMGs had significantly different transcription levels between the two groups. This is the first study to compare the DNA methylation profiles of embryos cultured in TLI and SI and will provide a foundation for evaluating the safety of TLI application in assisted reproductive technologies. However, further study with a larger cohort of samples was needed for the data validation.