Identification of a new locus and validation of previously reported loci showing differential methylation associated with smoking. The REGICOR study

Smoking increases the risk of many diseases and could act through changes in DNA methylation patterns. The aims of this study were to determine the association between smoking and DNA methylation throughout the genome at cytosine-phosphate-guanine (CpG) site level and genomic regions. A discovery cross-sectional epigenome-wide association study nested in the follow-up of the REGICOR cohort was designed and included 645 individuals. Blood DNA methylation was assessed using the Illumina HumanMethylation450 BeadChip. Smoking status was self-reported using a standardized questionnaire. We identified 66 differentially methylated CpG sites associated with smoking, located in 38 genes. In most of these CpG sites, we observed a trend among those quitting smoking to recover methylation levels typical of never smokers. A CpG site located in a novel smoking-associated gene (cg06394460 in LNX2) was hypomethylated in current smokers. Moreover, we validated two previously reported CpG sites (cg05886626 in THBS1, and cg24838345 in MTSS1) for their potential relation to atherosclerosis and cancer diseases, using several different approaches: CpG site methylation, gene expression, and plasma protein level determinations. Smoking was also associated with higher THBS1 gene expression but with lower levels of thrombospondin-1 in plasma. Finally, we identified differential methylation regions in 13 genes and in four non-coding RNAs. In summary, this study replicated previous findings and identified and validated a new CpG site located in LNX2 associated with smoking.     

Epigenetic loss of putative tumor suppressor SFRP3 correlates with poor prognosis of lung adenocarcinoma patients

Secreted frizzled related protein 3 (SFRP3) contains a cysteine-rich domain (CRD) that shares homology with Frizzled CRD and regulates WNT signaling. Independent studies showed epigenetic silencing of SFRP3 in melanoma and hepatocellular carcinoma. Moreover, a tumor suppressive function of SFRP3 was shown in androgen-independent prostate and gastric cancer cells. The current study is the first to investigate SFRP3 expression and its potential clinical impact on non-small cell lung carcinoma (NSCLC). WNT signaling components present on NSCLC subtypes were preliminary elucidated by expression data of The Cancer Genome Atlas (TCGA). We identified a distinct expression signature of relevant WNT signaling components that differ between adenocarcinoma (LUAD) and squamous cell carcinoma (LUSC). Of interest, canonical WNT signaling is predominant in LUAD samples and non-canonical WNT signaling is predominant in LUSC. In line, high SFRP3 expression resulted in beneficial clinical outcome for LUAD but not for LUSC patients. Furthermore, SFRP3 mRNA expression was significantly decreased in NSCLC tissue compared to normal lung samples. TCGA data verified the reduction of SFRP3 in LUAD and LUSC patients. Moreover, DNA hypermethylation of SFRP3 was evaluated in the TCGA methylation dataset resulting in epigenetic inactivation of SFRP3 expression in LUAD, but not in LUSC, and was validated by pyrosequencing of our NSCLC tissue cohort and in vitro demethylation experiments. Immunohistochemistry confirmed SFRP3 protein downregulation in primary NSCLC and indicated abundant expression in normal lung tissue. Two adenocarcinoma gain-of-function models were used to analyze the functional impact of SFRP3 on cell proliferation and regulation of CyclinD1 expression in vitro. Our results indicate that SFRP3 acts as a novel putative tumor suppressor gene in adenocarcinoma of the lung possibly regulating canonical WNT signaling.     

ABCC9, NKAPL,and TMEM132C are potential diagnostic and prognostic markers in triple‐negative breast cancer

Triple‐negative breast cancer (TNBC) is a highly heterogeneous disease. The aim of this study is to identify the diagnostic and poor prognostic signatures in TNBC by exploring the aberrant DNA methylation and gene expression. Differential expression and methylation analysis of the TNBC and paracancer samples from The Cancer Genome Atlas were performed. Gene set enrichment and protein–protein interaction (PPI) network analysis was used to explore the mechanisms of TNBC. Methylation‐gene expression correlation analysis was performed, and multivariate Cox analysis and receiver operating characteristics analysis were used to further screen the hub genes for TNBC. We identified 1,525 differentially expressed genes and 150 differentially methylated genes between TNBC and paracancer samples. About 96.64% of the methylation sites were located on the CpG island. A total of 17 Gene Ontology biological process terms and 18 signal pathways were significantly enriched. GNG4, GNG11, PENK, MAOA, and AOX1 were identified as the core genes of the PPI network. Methylation‐expression correlations revealed that ABCC9 (cg06951626), NKAPL (cg18675097, cg01031101, and cg17384889), and TMEM132C (cg03530754) showed promise as diagnostic and prognostic markers in TNBC. ABCC9 (cg06951626), NKAPL (cg18675097, cg01031101, and cg17384889), and TMEM132C (cg03530754) were potential diagnostic and prognostic markers in TNBC.     

Thioredoxin-interacting protein promotes activation and inflammation of monocytes with DNA demethylation in coronary artery disease

Numerous studies have demonstrated that thioredoxin-interacting protein (TXNIP) expression of peripheral blood leucocytes is increased in coronary artery disease (CAD). However, the molecular mechanism of this phenomenon remained unclear. DNA methylation plays important roles in the regulation of gene expression. Therefore, we speculated there might be a close association between the expression of TXNIP and methylation. In this study, we found that compared with controls, DNA methylation at cg19693031 was decreased in CAD, while mRNA expressions of TXNIP and inflammatory factors, NLRP3, IL-1β, IL-18, were increased. Methylation at cg19693031 was negatively associated with TXNIP expression in the cohort, THP-1 and macrophages/foam cells. Furthermore, Transwell assay and co-cultured adhesion assay were performed to investigate functions of TXNIP on the migration of THP-1 or the adhesion of THP-1 on the surface of endothelial cells, respectively. Notably, overexpressed TXNIP promoted the migration and adhesion of THP-1 cells and expressions of NLRP3, IL-18 and IL-1β. Oppositely, knock-down TXNIP inhibited the migration and adhesion of THP-1 and expressions of NLRP3, IL-18. In conclusion, increased TXNIP expression, related to cg19693031 demethylation orientates monocytes towards an inflammatory status through the NLRP3 inflammasome pathway involved in the development of CAD.     

DNA methylation regulates Microtubule-associated tumor suppressor 1 in human non-small cell lung carcinoma

Microtubule associated tumor suppressor 1 (MTUS1) has been recognized as a tumor suppressor gene in multiple cancers. However, the molecular mechanisms underlying the regulation of MTUS1 are yet to be investigated. This study aimed to clarify the significance of DNA methylation in silencing MTUS1 expression. We report that MTUS1 acts as tumor suppressor in non-small cell lung carcinoma (NSCLC). Analysis of in silico database and subsequent knockdown of DNMT1 suggested an inverse correlation between DNMT1 and MTUS1 function. Interestingly, increased methylation at MTUS1 promoter is associated with low expression of MTUS1. Treatment with DNA methyltransferases (DNMTs) inhibitor, 5-aza-2′-deoxycytidine (AZA) leads to both reduced promoter methylation accompanied with enrichment of H3K9Ac and enhanced MTUS1 expression. Remarkably, knockdown of MTUS1 showed increased proliferation and migration of NSCLC cells in contrast to diminished proliferation and migration, upon treatment with AZA. We concluded that low expression of MTUS1 correlates to DNA methylation and histone deacetylation in human NSCLC.     

The overexpression of GRASP might inhibit cell proliferation and invasion in hepatocellular carcinoma

This study aimed to validate the methylation of key genes in hepatocellular carcinoma (HCC) screened by bioinformatics analysis and explore whether they affected HCC cell proliferation, migration, and invasion. Using The Cancer Genome Atlas (TCGA) database, HCC-related differentially methylated positions (DMPs) were screened, genes corresponding to DMPs were selected, and prognosis-related genes were identified. A representative DMP was used to divide the DMPs into hyper- and hypomethylated groups. Expression of key genes in cell lines was detected using quantitative real-time polymerase chain reaction and western blot analysis. After treatment of HepG2 cells with 5-Aza-2′-deoxycytidine (5-Aza-DC), gene expression was observed. Bisulfite sequencing PCR assay was used to detect methylation frequency. Overexpressed GRASP lentiviral vectors were constructed to analyze their influence on cell proliferation, migration, and invasion using cell counting kit-8 and transwell assays. Forty-three HCC prognosis-related genes were screened using the TCGA database. cg00249511 (SCT) was used to divide the DMPs into hyper- and hypomethylated groups, distinguishing between high- and low-risk samples. The prognosis survival model constructed using 12 genes revealed the prognosis type. GRASP messenger RNA was downregulated in HepG2 and upregulated after 5-Aza-DC treatment. In HCC tissues, methylation frequency of GRASP was upregulated. GRASP overexpression inhibited HepG2 cell proliferation, invasion, and G-CSFR expression. Thus, GRASP might be a prognosis-related gene controlled by methylation.     

A two-CpG-based prognostic signature for oral squamous cell carcinoma overall survival

Oral squamous cell carcinoma (OSCC) represents one of the most common head and neck cancer that with dire prognosis due partly to the lack of reliable prognostic biomarker. Here, we aimed to develop a CpG site–based prognostic signature through which we could accurately predict overall survival (OS) of patients with OSCC. We obtained OSCC-related DNA methylation and gene expression data sets from the public accessible Gene Expression Omnibus. Correlations between methylation level of CpG sites and OS of patients with OSCC were assessed by univariate Cox regression analysis followed by robust likelihood-based survival analysis on those CpG sites with permutation P < 0.05 for further screening the optimal CpG sites for OSCC OS prediction based on the risk score formula that composed of the methylation level of optimal CpG sites weighted by their regression coefficients. Besides, differential expression genes (DEGs) and differential methylation genes (DMGs) in OSCC samples compared with normal samples were obtained and shared genes were considered as vital genes in OSCC tumorgenesis and progression. As a result, two CpG sites including cg17892178 and cg17378966 that located in NID2 and IDO1, respectively, were identified as the optimal prognostic signatures for OSCC OS. In addition, 12 overlapping genes between DEGs and DMGs that closely associated with inflammation or blood and tissue development–related biological processes were obtained. In conclusions, this study should provide valuable signatures for OSCC diagnosis and treatment.     

Downregulation of miR-16 promotes growth and motility by targeting HDGF in non-small cell lung cancer cells

MicroRNAs play important roles in the development and progression of non-small cell lung cancer (NSCLC). miR-16 functions as a tumor-suppressor and is inhibited in several malignancies. Herein, we validated that miR-16 is downregulated in NSCLC tissue samples and cell lines. Ectopic expression of miR-16 significantly inhibited cell proliferation and colony formation. Moreover, miR-16 suppressed cell migration and invasion in NSCLC cells. Hepatoma-derived growth factor (HDGF) was found to be a direct target of miR-16 in NSCLC cell lines. Rescue experiments showed that the suppressive effect of miR-16 on cell proliferation, colony formation, migration, and invasion is partially mediated by inhibiting HDGF expression. This study indicates that miR-16 might be associated with NSCLC progression, and suggests an essential role for miR-16 in NSCLC.     

An integrated analysis of miRNA and mRNA expressions in non-small cell lung cancers

Using DNA microarrays, we generated both mRNA and miRNA expression data from 6 non-small cell lung cancer (NSCLC) tissues and their matching normal control from adjacent tissues to identify potential miRNA markers for diagnostics. We demonstrated that hsa-miR-96 is significantly and consistently up-regulated in all 6 NSCLCs. We validated this result in an independent set of 35 paired tumors and their adjacent normal tissues, as well as their sera that are collected before surgical resection or chemotherapy, and the results suggested that hsa-miR-96 may play an important role in NSCLC development and has great potential to be used as a noninvasive marker for diagnosing NSCLC. We predicted potential miRNA target mRNAs based on different methods (TargetScan and miRanda). Further classification of miRNA regulated genes based on their relationship with miRNAs revealed that hsa-miR-96 and certain other miRNAs tend to down-regulate their target mRNAs in NSCLC development, which have expression levels permissive to direct interaction between miRNAs and their target mRNAs. In addition, we identified a significant correlation of miRNA regulation with genes coincide with high density of CpG islands, which suggests that miRNA may represent a primary regulatory mechanism governing basic cellular functions and cell differentiations, and such mechanism may be complementary to DNA methylation in repressing or activating gene expression.     

The relationship among Girdin DNA methylation,its high expression,and immune infiltration in hepatocellular carcinoma: Clues from in silico analysis

Objective: The aim of the present study was to explore the relationship among Girdin DNA methylation, its high expression, and immune infiltration in human hepatocellular carcinoma (HCC).Materials and methods: The Cancer Genome Atlas (TCGA), Gene Expression Omnibus (GEO), and International Cancer Genome Consortium (ICGC) databases were used to compare Girdin mRNA expression between HCC tissues and normal tissues, and determine the relationship between Girdin expression and HCC prognosis. TCGA database was also used to analyze the expression of Girdin and its methylation status, as well as the relationship between Girdin DNA methylation and HCC prognosis. The Tumor IMmune Estimation Resource (TIMER) database was used to explore the correlation between Girdin expression and HCC immune infiltration.Results: Girdin expression was elevated in HCC tissues compared with that in normal tissues. The degree of methylation at cg03188526, a CpG site in the Girdin gene body, was positively correlated with Girdin mRNA expression, while high Girdin expression and cg03188526 hypermethylation were both correlated with poor HCC prognosis. Additionally, HCC tissue with high Girdin expression exhibited abundant immune infiltration, and the high Girdin expression was associated with a worse prognosis in macrophage-enriched HCC specimens.Conclusion: Our findings indicated that Girdin likely functions as an oncogene in HCC and that hypermethylation at cg03188526 in the Girdin gene body may explain the high Girdin expression levels in HCC tissue. Furthermore, we report for the first time that the adverse effects of high Girdin expression in HCC patients may be partially mediated by tumor macrophage infiltration.     

Integrative analysis identifies potential DNA methylation biomarkers for pan-cancer diagnosis and prognosis

DNA methylation status is closely associated with diverse diseases, and is generally more stable than gene expression, thus abnormal DNA methylation could be important biomarkers for tumor diagnosis, treatment and prognosis. However, the signatures regarding DNA methylation changes for pan-cancer diagnosis and prognosis are less explored. Here we systematically analyzed the genome-wide DNA methylation patterns in diverse TCGA cancers with machine learning. We identified seven CpG sites that could effectively discriminate tumor samples from adjacent normal tissue samples for 12 main cancers of TCGA (1216 samples, AUC > 0.99). Those seven potential diagnostic biomarkers were further validated in the other 9 different TCGA cancers and 4 independent datasets (AUC > 0.92). Three out of the seven CpG sites were correlated with cell division, DNA replication and cell cycle. We also identified 12 CpG sites that can effectively distinguish 26 different cancers (7605 samples), and the result was repeatable in independent datasets as well as two disparate tumors with metastases (micro-average AUC > 0.89). Furthermore, a series of potential signatures that could significantly predict the prognosis of tumor patients for 7 different cancer were identified via survival analysis (p-value < 1e-4). Collectively, DNA methylation patterns vary greatly between tumor and adjacent normal tissues, as well as among different types of cancers. Our identified signatures may aid the decision of clinical diagnosis and prognosis for pan-cancer and the potential cancer-specific biomarkers could be used to predict the primary site of metastatic breast and prostate cancers.     

Correlation of DNA methylation levels in blood and saliva DNA in young girls of the LEGACY Girls study

Many epidemiologic studies of environmental exposures and disease susceptibility measure DNA methylation in white blood cells (WBC). Some studies are also starting to use saliva DNA as it is usually more readily available in large epidemiologic studies. However, little is known about the correlation of methylation between WBC and saliva DNA. We examined DNA methylation in three repetitive elements, Sat2, Alu, and LINE-1, and in four CpG sites, including AHRR (cg23576855, cg05575921), cg05951221 at 2q37.1, and cg11924019 at CYP1A1, in 57 girls aged 6–15 years with blood and saliva collected on the same day. We measured all DNA methylation markers by bisulfite-pyrosequencing, except for Sat2 and Alu, which were measured by the MethyLight assay. Methylation levels measured in saliva DNA were lower than those in WBC DNA, with differences ranging from 2.8% for Alu to 14.1% for cg05575921. Methylation levels for the three repetitive elements measured in saliva DNA were all positively correlated with those in WBC DNA. However, there was a wide range in the Spearman correlations, with the smallest correlation found for Alu (0.24) and the strongest correlation found for LINE-1 (0.73). Spearman correlations for cg05575921, cg05951221, and cg11924019 were 0.33, 0.42, and 0.79, respectively. If these findings are replicated in larger studies, they suggest that, for selected methylation markers (e.g., LINE-1), methylation levels may be highly correlated between blood and saliva, while for others methylation markers, the levels may be more tissue specific. Thus, in studies that differ by DNA source, each interrogated site should be separately examined in order to evaluate the correlation in DNA methylation levels across DNA sources.     

Non-linear patterns in age-related DNA methylation may reflect CD4+ T cell differentiation

DNA methylation (DNAm) is an important epigenetic process involved in the regulation of gene expression. While many studies have identified thousands of loci associated with age, few have differentiated between linear and non-linear DNAm trends with age. Non-linear trends could indicate early- or late-life gene regulatory processes. Using data from the Illumina 450K array on 336 human peripheral blood samples, we identified 21 CpG sites that associated with age (P<1.03E-7) and exhibited changing rates of DNAm change with age (P<1.94E-6). For 2 of these CpG sites (cg07955995 and cg22285878), DNAm increased with age at an increasing rate, indicating that differential DNAm was greatest among elderly individuals. We observed significant replication for both CpG sites (P<5.0E-8) in a second set of peripheral blood samples. In 8 of 9 additional data sets comprising samples of monocytes, T cell subtypes, and brain tissue, we observed a pattern directionally consistent with DNAm increasing with age at an increasing rate, which was nominally significant in the 3 largest data sets (4.3E-15<P<0.039). cg07955995 and cg22285878 reside in the promoter region of KLF14, which encodes a protein involved in immune cell differentiation via the repression of FOXP3. These findings may suggest a possible role for cg07955995 and cg22285878 in immunosenescence.     

Correlation of DNA methylation levels in blood and saliva DNA in young girls of the LEGACY Girls study

Many epidemiologic studies of environmental exposures and disease susceptibility measure DNA methylation in white blood cells (WBC). Some studies are also starting to use saliva DNA as it is usually more readily available in large epidemiologic studies. However, little is known about the correlation of methylation between WBC and saliva DNA. We examined DNA methylation in three repetitive elements, Sat2, Alu, and LINE-1, and in four CpG sites, including AHRR (cg23576855, cg05575921), cg05951221 at 2q37.1, and cg11924019 at CYP1A1, in 57 girls aged 6–15 years with blood and saliva collected on the same day. We measured all DNA methylation markers by bisulfite-pyrosequencing, except for Sat2 and Alu, which were measured by the MethyLight assay. Methylation levels measured in saliva DNA were lower than those in WBC DNA, with differences ranging from 2.8% for Alu to 14.1% for cg05575921. Methylation levels for the three repetitive elements measured in saliva DNA were all positively correlated with those in WBC DNA. However, there was a wide range in the Spearman correlations, with the smallest correlation found for Alu (0.24) and the strongest correlation found for LINE-1 (0.73). Spearman correlations for cg05575921, cg05951221, and cg11924019 were 0.33, 0.42, and 0.79, respectively. If these findings are replicated in larger studies, they suggest that, for selected methylation markers (e.g., LINE-1), methylation levels may be highly correlated between blood and saliva, while for others methylation markers, the levels may be more tissue specific. Thus, in studies that differ by DNA source, each interrogated site should be separately examined in order to evaluate the correlation in DNA methylation levels across DNA sources.     

Aberrant DNA methylation of WNT pathway genes in the development and progression of CIMP-negative colorectal cancer

The WNT signaling pathway has an essential role in colorectal carcinogenesis and progression, which involves a cascade of genetic and epigenetic changes. We aimed to analyze DNA methylation affecting the WNT pathway genes in colorectal carcinogenesis in promoter and gene body regions using whole methylome analysis in 9 colorectal cancer, 15 adenoma, and 6 normal tumor adjacent tissue (NAT) samples by methyl capture sequencing. Functional methylation was confirmed on 5-aza-2′-deoxycytidine-treated colorectal cancer cell line datasets. In parallel with the DNA methylation analysis, mutations of WNT pathway genes (APC, β-catenin/CTNNB1) were analyzed by 454 sequencing on GS Junior platform. Most differentially methylated CpG sites were localized in gene body regions (95% of WNT pathway genes). In the promoter regions, 33 of the 160 analyzed WNT pathway genes were differentially methylated in colorectal cancer vs. normal, including hypermethylated AXIN2, CHP1, PRICKLE1, SFRP1, SFRP2, SOX17, and hypomethylated CACYBP, CTNNB1, MYC; 44 genes in adenoma vs. NAT; and 41 genes in colorectal cancer vs. adenoma comparisons. Hypermethylation of AXIN2, DKK1, VANGL1, and WNT5A gene promoters was higher, while those of SOX17, PRICKLE1, DAAM2, and MYC was lower in colon carcinoma compared to adenoma. Inverse correlation between expression and methylation was confirmed in 23 genes, including APC, CHP1, PRICKLE1, PSEN1, and SFRP1. Differential methylation affected both canonical and noncanonical WNT pathway genes in colorectal normal-adenoma-carcinoma sequence. Aberrant DNA methylation appears already in adenomas as an early event of colorectal carcinogenesis.