Genome-wide analysis of DNA methylation and the gene expression change in lung cancer. J Thorac Oncol. 2012;7:20-33. [PMID: 22011669 DOI: 10.1097/jto.0b013e3182307f62]

DNA methylation heterogeneity correlates with field cancerization and prognosis in lung adenocarcinoma patients

BACKGROUND

Lung adenocarcinoma (LUAD) is the most common histological subtype of lung cancer. The distinctive genetic and epigenetic modifications in tumors and paired non-malignant samples, such as adjacent peri-tumor and tumor-distant normal lung tissues, have not been adequately studied.

METHODS

We recruited 57 patients with resectable stage I-III LUAD and collected matched samples of the primary tumor, peri-tumoral tissues, and tumor-distant normal lung tissue. We performed bisulfite sequencing using a custom methylation panel to profile DNA methylation levels and obtained somatic variation landscape through targeted next-generation sequencing (NGS). We attempted to identify differential methylation blocks (DMBs) between the tumor, peri-tumor, and normal tissues.

RESULTS

We analyzed the DNA methylation patterns of matched tumor, peri-tumor, and normal lung tissue samples from 57 LUAD patients. No significantly different methylation blocks were found between peri-tumoral and normal tissues, while they both exhibited distinct methylation profiles compared to tumor tissues. A total of 1329 tumor-specific DMBs, which are potentially associated with aberrant gene expression in LUAD, were identified. Utilizing a consensus clustering algorithm, we classified the tumor samples into two subgroups (C1 and C2) based on distinct methylation profiles, independent of the patient's sex, tumor stage, smoking history, and tumor cell fraction. The C2 subgroup exhibited a higher malignancy density ratio (MD ratio), suggesting a more pronounced degree of field cancerization, while the C1 subgroup was characterized by a higher frequency of EGFR mutations. The DMBs between the two subgroups were enriched in the calcium signaling pathway. Notably, P2RX2 shows significant hypermethylation in the C2 subgroup, and its low expression in the external The Cancer Genome Atlas (TCGA) cohort may correlate with reduced overall survival in LUAD patients.

CONCLUSION

Our findings revealed distinct methylation patterns between tumor and pre-malignant samples, such as peri-tumor and normal tissues. Moreover, our study suggests that distinct clustering based on DNA methylation may indicate different prognoses in LUAD patients.

Breast cancer cells have an increased ferroptosis risk induced by system xc- blockade after deliberately downregulating CYTL1 to mediate malignancy

Cytokine-like protein 1 (CYTL1) expression is deliberately downregulated during the progression of multiple types of cancers, especially breast cancer. However, the metabolic characteristics of cancer progression remain unclear. Here, we uncovered a risk of breast cancer cells harboring low CYTL1 expression, which is metabolically controlled during malignant progression. We performed metabolism comparison and revealed that breast cancer cells with low CYTL1 expression have highly suppressed transsulfuration activity that is driven by cystathionine β-synthase (CBS) and contributes to de novo cysteine synthesis. Mechanistically, CYTL1 activated Nrf2 by promoting autophagic Keap1 degradation, and Nrf2 subsequently transactivated CBS expression. Due to the lack of cellular cysteine synthesis, breast cancer cells with low CYTL1 expression showed hypersensitivity to system xc- blockade-induced ferroptosis in vitro and in vivo. Silencing CBS counteracted CYTL1-mediated ferroptosis resistance. Our results show the importance of exogeneous cysteine in breast cancer cells with low CYTL1 expression and highlight a potential metabolic vulnerability to target.

RNA-Binding Proteins as a Molecular Link between COPD and Lung Cancer

Chronic obstructive pulmonary disease (COPD) represents an independent risk factor for lung cancer development. Accelerated cell senescence, induced by oxidative stress and inflammation, is a common pathogenic determinant of both COPD and lung cancer. The post transcriptional regulation of genes involved in these processes is finely regulated by RNA-binding proteins (RBPs), which regulate mRNA turnover, subcellular localization, splicing and translation. Multiple pro-inflammatory mediators (including cytokines, chemokines, proteins, growth factors and others), responsible of lung microenvironment alteration, are regulated by RBPs. Several mouse models have shown the implication of RBPs in multiple mechanisms that sustain chronic inflammation and neoplastic transformation. However, further studies are required to clarify the role of RBPs in the pathogenic mechanisms shared by lung cancer and COPD, in order to identify novel biomarkers and therapeutic targets. This review will therefore focus on the studies collectively indicating the role of RBPs in oxidative stress and chronic inflammation as common pathogenic mechanisms shared by lung cancer and COPD.

Pre-clinical lung squamous cell carcinoma mouse models to identify novel biomarkers and therapeutic interventions

Primary lung carcinoma or lung cancer (LC) is classified into small-cell or non-small-cell (NSCLC) lung carcinoma. Lung squamous cell carcinoma (LSCC) is the second most common subtype of NSCLC responsible for 30% of all LCs, and its survival remains low with only 24% of patients living for five years or longer post-diagnosis primarily due to the advanced stage of tumors at the time of diagnosis. The pathogenesis of LSCC is still poorly understood and has hampered the development of effective diagnostics and therapies. This review highlights the known risk factors, genetic and epigenetic alterations, miRNA biomarkers linked to the development and diagnosis of LSCC and the lack of therapeutic strategies to target specifically LSCC. We will also discuss existing animal models of LSCC including carcinogen induced, transgenic and xenograft mouse models, and their advantages and limitations along with the chemopreventive studies and molecular studies conducted using them. The importance of developing new and improved mouse models will also be discussed that will provide further insights into the initiation and progression of LSCC, and enable the identification of new biomarkers and therapeutic targets.

Bioinformatics-based analysis reveals elevated CYTL1 as a potential therapeutic target for BRAF-mutated melanoma

Introduction: Despite many recent emerging therapeutic modalities that have prolonged the survival of melanoma patients, the prognosis of melanoma remains discouraging, and further understanding of the mechanisms underlying melanoma progression is needed. Melanoma patients often have multiple genetic mutations, with BRAF mutations being the most common. In this study, public databases were exploited to explore a potential therapeutic target for BRAF-mutated melanoma. Methods: In this study, we analyzed differentially expressed genes (DEGs) in normal tissues and melanomas, Braf wild-type and Braf mutant melanomas using information from TCGA databases and the GEO database. Subsequently, we analyzed the differential expression of CYTL1 in various tumor tissues and its effect on melanoma prognosis, and resolved the mutation status of CYTL1 and its related signalling pathways. By knocking down CYTL1 in melanoma cells, the effects of CYTL1 on melanoma cell proliferation, migration and invasion were further examined by CCK8 assay, Transwell assay and cell migration assay. Results: 24 overlapping genes were identified by analyzing DEGs common to melanoma and normal tissue, BRAF-mutated and BRAF wild-type melanoma. Among them, CYTL1 was highly expressed in melanoma, especially in BRAF-mutated melanoma, and the high expression of CYTL1 was associated with epithelial-mesenchymal transition (EMT), cell cycle, and cellular response to UV. In melanoma patients, especially BRAF-mutated melanoma patients, clinical studies showed a positive correlation between increased CYTL1 expression and shorter overall survival (OS) and disease-free survival (DFS). In vitro experiments further confirmed that the knockdown of CYTL1 significantly inhibited the migration and invasive ability of melanoma cells. Conclusion: CYTL1 is a valuable prognostic biomarker and a potentially effective therapeutic target in melanoma, especially BRAF-mutated melanoma.

Regulatory network and targeted interventions for CCDC family in tumor pathogenesis

CCDC (coiled-coil domain-containing) is a coiled helix domain that exists in natural proteins. There are about 180 CCDC family genes, encoding proteins that are involved in intercellular transmembrane signal transduction and genetic signal transcription, among other functions. Alterations in expression, mutation, and DNA promoter methylation of CCDC family genes have been shown to be associated with the pathogenesis of many diseases, including primary ciliary dyskinesia, infertility, and tumors. In recent studies, CCDC family genes have been found to be involved in regulation of growth, invasion, metastasis, chemosensitivity, and other biological behaviors of malignant tumor cells in various cancer types, including nasopharyngeal carcinoma, lung cancer, colorectal cancer, and thyroid cancer. In this review, we summarize the involvement of CCDC family genes in tumor pathogenesis and the relevant upstream and downstream molecular mechanisms. In addition, we summarize the potential of CCDC family genes as tumor therapy targets. The findings discussed here help us to further understand the role and the therapeutic applications of CCDC family genes in tumors.

Development of a novel autophagy-related gene model for gastric cancer prognostic prediction

Gastric cancer (GC) is a major global health issue and one of the leading causes of tumor-associated mortality worldwide. Autophagy is thought to play a critical role in the development and progression of GC, and this process is controlled by a set of conserved regulators termed autophagy-related genes (ATGs). However, the complex contribution of autophagy to cancers is not completely understood. Accordingly, we aimed to develop a prognostic model based on the specific role of ATGs in GC to improve the prediction of GC outcomes. First, we screened 148 differentially expressed ATGs between GC and normal tissues in The Cancer Genome Atlas (TCGA) cohort. Consensus clustering in these ATGs was performed, and based on that, 343 patients were grouped into two clusters. According to Kaplan–Meier survival analysis, cluster C2 had a worse prognosis than cluster C1. Then, a disease risk model incorporating nine differentially expressed ATGs was constructed based on the least absolute shrinkage and selection operator (LASSO) regression analysis, and the ability of this model to stratify patients into high- and low-risk groups was verified. The predictive value of the model was confirmed using both training and validation cohorts. In addition, the results of functional enrichment analysis suggested that GC risk is correlated with immune status. Moreover, autophagy inhibition increased sensitivity to cisplatin and exacerbated reactive oxygen species accumulation in GC cell lines. Collectively, the results indicated that this novel constructed risk model is an effective and reliable tool for predicting GC outcomes and could help with individual treatment through ATG targeting.

Epigenome-wide DNA methylation signature of plasma zinc and their mediation roles in the association of zinc with lung cancer risk

Essential trace element zinc is associated with decreased lung cancer risk, but underlying mechanisms remain unclear. This study aimed to investigate role of DNA methylation in zinc-lung cancer association. We conducted a case-cohort study within prospective Dongfeng-Tongji cohort, including 359 incident lung cancer cases and a randomly selected sub-cohort of 1399 participants. Epigenome-wide association study (EWAS) was used to examine association of plasma zinc with DNA methylation in peripheral blood. For the zinc-related CpGs, their mediation effects on zinc-lung cancer association were assessed; their diagnostic performance for lung cancer was testified in the case-cohort study and further validated in another 126 pairs of lung cancer case-control study. We identified 28 CpGs associated with plasma zinc at P < 1.0 × 10^-5 and seven of them (cg07077080, cg01077808, cg17749033, cg15554270, cg26125625, cg10669424, and cg15409013 annotated to GSR, CALR3, SLC16A3, PHLPP2, SLC12A8, VGLL4, and ADAMTS16, respectively) were associated with incident risk of lung cancer. Moreover, the above seven CpGs were differently methylated between 126 pairs of lung cancer and adjacent normal lung tissues and had the same directions with EWAS of zinc. They could mediate a separate 7.05%∼22.65% and a joint 29.42% of zinc-lung cancer association. Compared to using traditional factors, addition of methylation risk score exerted improved discriminations for lung cancer both in case-cohort study [area under the curve (AUC) = 0.818 vs. 0.738] and in case-control study (AUC = 0.816 vs. 0.646). Our results provide new insights for the biological role of DNA methylation in the inverse association of zinc with incident lung cancer.

Intracellular CYTL1, a novel tumor suppressor, stabilizes NDUFV1 to inhibit metabolic reprogramming in breast cancer

Loss-of-function mutations frequently occur in tumor suppressor genes, i.e., p53, during the malignant progression of various cancers. Whether any intrinsic suppressor carries a rare mutation is largely unknown. Here, we demonstrate that intracellular cytokine-like protein 1 (CYTL1) plays a key role in preventing the robust glycolytic switching characteristic of breast cancer. A low intracellular CYTL1 level, not its mutation, is required for metabolic reprogramming. Breast cancer cells expressing an intracellular form of CYTL1 lacking a 1-22 aa signal peptide, ΔCYTL1, show significantly attenuated glucose uptake and lactate production, which is linked to the inhibition of cell growth and metastasis in vitro and in vivo. Mechanistically, CYTL1 competitively binds the N-terminal sequence of NDUFV1 to block MDM2-mediated degradation by the proteasome, leading to the stability of the NDUFV1 protein. In addition to inducing increased NAD⁺ levels, NDUFV1 interacts with Src to attenuate LDHA phosphorylation at tyrosine 10 and reduce lactate production. Our results reveal, for the first time, that CYTL1 is a novel tumor suppressor. Its function in reversing metabolic reprogramming toward glycolysis may be very important for the development of novel antitumor strategies.

Aberrant Methylation of SLIT2 Gene in Plasma Cell-Free DNA of Non-Small Cell Lung Cancer Patients

Simple Summary

Despite significant advances in the detection, prevention, and treatment of lung cancer, the prognosis of the patients is still very poor due in part to micrometastasis of cancer cells to surrounding tissues at the time of diagnosis. Therefore, identifying biomarkers for early detection of lung cancer is very important for prolonging the lifespan of patients with lung cancer. The methylation statuses of SLIT1, SLIT2, SLIT3 genes were analyzed in bronchial washing, bronchial biopsy, sputum, tumor and matched normal tissues, or plasma samples obtained from a total of 208 non-small cell lung cancer (NSCLC) patients and 121 cancer-free patients to understand the feasibility of the genes as biomarkers for early detection and survival prediction of NSCLC. The present study suggests that aberrant methylation of SLIT2 in plasma cell-free DNA might be a potential biomarker for the early detection and prognosis prediction of NSCLC patient.

Abstract

This study aimed to understand aberrant methylation of SLITs genes as a biomarker for the early detection and prognosis prediction of non-small cell lung cancer (NSCLC). Methylation levels of SLITs were determined using the Infinium HumanMethylation450 BeadChip or pyrosequencing. Five CpGs at the CpG island of SLIT1, SLIT2 or SLIT3 genes were significantly (Bonferroni corrected p < 0.05) hypermethylated in tumor tissues obtained from 42 NSCLC patients than in matched normal tissues. Methylation levels of these CpGs did not differ significantly between bronchial washings obtained from 76 NSCLC patients and 60 cancer-free patients. However, methylation levels of SLIT2 gene were significantly higher in plasma cell-free DNA of 72 NSCLC patients than in that of 61 cancer-free patients (p = 0.001, Wilcoxon rank sum test). Prediction of NSCLC using SLIT2 methylation was achieved with a sensitivity of 73.7% and a specificity of 61.9% in a plasma test dataset (N = 40). A Cox proportional hazards model showed that SLIT2 hypermethylation in plasma cell-free DNA was significantly associated with poor recurrence-free survival (hazards ratio = 2.19, 95% confidence interval = 1.21–4.36, p = 0.01). The present study suggests that aberrant methylation of SLIT2 in plasma cell-free DNA is a valuable biomarker for the early detection of NSCLC and prediction of recurrence-free survival. However, further research is needed with larger sample size to confirm results.

Epigenome-wide DNA methylation signature of benzo[a]pyrene exposure and their mediation roles in benzo[a]pyrene-associated lung cancer development

Benzo[a]pyrene (B[a]P) is a typical carcinogen associated with increased lung cancer risk, but the underlying mechanisms remain unclear. This study aimed to investigate epigenome-wide DNA methylation associated with B[a]P exposure and their mediation effects on B[a]P-lung cancer association in two lung cancer case-control studies of 462 subjects. Their plasma levels of benzo[a]pyrene diol epoxide-albumin (BPDE-Alb) adducts and genome-wide DNA methylations were separately detected in peripheral blood by using enzyme-linked immunosorbent assay (ELISA) and genome-wide methylation arrays. The epigenome-wide meta-analysis was performed to analyze the associations between BPDE-Alb adducts and DNA methylations. Mediation analysis was applied to assess effect of DNA methylation on the B[a]P-lung cancer association. We identified 15 CpGs associated with BPDE-Alb adducts (P_-meta < 1.0 × 10^-5), among which the methylation levels at five loci (cg06245338, cg24256211, cg15107887, cg02211741, and cg04354393 annotated to UBE2O, SAMD4A, ACBD6, DGKZ, and SLFN13, respectively) mediated a separate 38.5%, 29.2%, 41.5%, 47.7%, 56.5%, and a joint 58.2% of the association between BPDE-Alb adducts and lung cancer risk. Compared to the traditional factors [area under the curve (AUC) = 0.788], addition of these CpGs exerted improved discriminations for lung cancer, with AUC ranging 0.828-0.861. Our results highlight DNA methylation alterations as potential mediators in lung tumorigenesis induced by B[a]P exposure.

DNA methylome analysis identifies BMI-related epigenetic changes associated with non-small cell lung cancer susceptibility

Background

Body mass index (BMI) has been reported to be inversely associated with incident risk of non‐small cell lung cancer (NSCLC). However, the underlying mechanism is still unclear. This study aimed to investigate the role of DNA methylation in the relationship between BMI and NSCLC.

Methods

We carried out a genome‐wide DNA methylation study of BMI in peripheral blood among 2266 Chinese participants by using Illumina Methylation arrays. For the BMI‐related DNA methylation changes, their associations with NSCLC risk were further analyzed and their mediation effects on BMI‐NSCLC association were also evaluated.

Results

The methylation levels of four CpGs (cg12593793, cg17061862, cg11024682, and cg06500161, annotated to LMNA, ZNF143, SREBF1, and ABCG1, respectively) were found to be significantly associated with BMI. Methylation levels of cg12593793, cg11024682, and cg06500161 were observed to be inversely associated with NSCLC risk [OR (95%CI) =0.22 (0.16, 0.31), 0.39 (0.30, 0.50), and 0.66 (0.53, 0.82), respectively]. Additionally, cg11024682 in SREBF1 and cg06500161 in ABCG1 mediated 45.3% and 19.5% of the association between BMI and decreased NSCLC risk, respectively.

Conclusions

In this study, we identified four DNA methylation sites associated with BMI in the Chinese populations at the genome‐wide significant level. We also found that the BMI‐related methylations of SREBF1 and ABCG1 could mediate about a quintile‐to‐half of the effect of BMI on reduced NSCLC risk, which adds a potential mechanism underlying this association.

Composition and function of ciliary inner-dynein-arm subunits studied in Chlamydomonas reinhardtii

Motile cilia (also interchangeably called "flagella") are conserved organelles extending from the surface of many animal cells and play essential functions in eukaryotes, including cell motility and environmental sensing. Large motor complexes, the ciliary dyneins, are present on ciliary outer-doublet microtubules and drive movement of cilia. Ciliary dyneins are classified into two general types: the outer dynein arms (ODAs) and the inner dynein arms (IDAs). While ODAs are important for generation of force and regulation of ciliary beat frequency, IDAs are essential for control of the size and shape of the bend, features collectively referred to as waveform. Also, recent studies have revealed unexpected links between IDA components and human diseases. In spite of their importance, studies on IDAs have been difficult since they are very complex and composed for several types of IDA motors, each unique in composition and location in the axoneme. Thanks in part to genetic, biochemical, and structural analysis of Chlamydomonas reinhardtii, we are beginning to understand the organization and function of the ciliary IDAs. In this review, we summarize the composition of Chlamydomonas IDAs particularly focusing on each subunit, and discuss the assembly, conservation, and functional role(s) of these IDA subunits. Furthermore, we raise several additional questions/challenges regarding IDAs, and discuss future perspectives of IDA studies.

Upregulation of glucosamine-phosphate N-acetyltransferase 1 is a promising diagnostic and predictive indicator for poor survival in patients with lung adenocarcinoma

Lung adenocarcinoma, a type of non-small cell lung cancer, is the leading cause of cancer death worldwide. Great efforts have been made to identify the underlying mechanism of adenocarcinoma, especially in relation to oncogenes. The present study by integrating computational analysis with western blotting, aimed to understand the role of the upregulation of glucosamine-phosphate N-acetyltransferase 1 (GNPNAT1) in carcinogenesis. In the present study, publicly available gene expression profiles and clinical data were downloaded from The Cancer Genome Atlas to determine the role of GNPNAT1 in lung adenocarcinoma (LUAD). In addition, the association between LUAD susceptibility and GNPNAT1 upregulation were analyzed using Wilcoxon signed-rank test and logistic regression analysis. In LUAD, GNPNAT1 upregulation was significantly associated with disease stage [odds ratio (OR)=2.92, stage III vs. stage I], vital status (dead vs. alive, OR=1.89), cancer status (tumor status vs. tumor-free status, OR=1.85) and N classification (yes vs. no, OR=1.75). Cox regression analysis and the Kaplan-Meier method were utilized to evaluate the association between GNPNAT1 expression and overall survival (OS) time in patients with LUAD. The results demonstrated that patients with increased GNPNAT1 expression levels exhibited a reduced survival rate compared with those with decreased expression levels (P=8.9×10⁻⁵). In addition, Cox regression analysis revealed that GNPNAT1 upregulation was significantly associated with poor OS time [hazard ratio (HR): 1.07; 95% confidence interval (CI): 1.04–1.10; P<0.001]. The gene set enrichment analysis revealed that ‘cell cycle’, ‘oocyte meiosis’, ‘pyrimidine mediated metabolism’, ‘ubiquitin mediated proteolysis’, ‘one carbon pool by folate’, ‘mismatch repair progesterone-mediated oocyte maturation’ and ‘basal transcription factors purine metabolism’ were differentially enriched in the GNPNAT1 high-expression samples compared with GNPNAT1 low-expression samples. The aforementioned pathways are involved in the pathogenesis of LUAD. The findings of the present study suggested that GNPNAT1 upregulation may be considered as a promising diagnostic and prognostic biomarker in patients with LUAD. In addition, the aforementioned pathways may be pivotal pathways perturbed by the abnormal expression of GNPNAT1 in LUAD. The findings of the present study demonstrated the therapeutic value of the regulation of GNPNAT1 in lung adenocarcinoma.

SETDB1 Overexpression Sets an Intertumoral Transcriptomic Divergence in Non-small Cell Lung Carcinoma

An increasing volume of evidence suggests that SETDB1 plays a role in the tumorigenesis of various cancers, classifying SETDB1 as an oncoprotein. However, owing to its numerous protein partners and their global-scale effects, the molecular mechanism underlying SETDB1-involved oncogenesis remains ambiguous. In this study, using public transcriptome data of lung adenocarcinoma (ADC) and squamous-cell carcinoma (SCC), we compared tumors with high-level SETDB1 (SH) and those with low-level SETDB1 (comparable with normal samples; SL). The results of principal component analysis revealed a transcriptomic distinction and divergence between the SH and SL samples in both ADCs and SCCs. The results of gene set enrichment analysis indicated that genes involved in the “epithelial–mesenchymal transition,” “innate immune response,” and “autoimmunity” collections were significantly depleted in SH tumors, whereas those involved in “RNA interference” collections were enriched. Chromatin-modifying genes were highly expressed in SH tumors, and the variance in their expression was incomparably high in SCC-SH, which suggested greater heterogeneity within SCC tumors. DNA methyltransferase genes were also overrepresented in SH samples, and most differentially methylated CpGs (SH/SL) were undermethylated in a highly biased manner in ADCs. We identified interesting molecular signatures associated with the possible roles of SETDB1 in lung cancer. We expect these SETDB1-associated molecular signatures to facilitate the development of biologically relevant targeted therapies for particular types of lung cancer.

CYTL1 Promotes the Activation of Neutrophils in a Sepsis Model

As a novel cytokine, cytokine-like 1 (CYTL1) is a classical secretory protein, and its potential biological function remains to be determined. In this study, we found that expression of CYTL1 was upregulated in neutrophils upon inflammatory stimuli. We demonstrated that CYTL1 enhanced phagocytosis of Escherichia coli by activated neutrophils both in vivo and in vitro through phosphorylation of protein kinase B (Akt). CYTL1-induced chemotactic activity in lipopolysaccharide (LPS) stimulated neutrophils, and the mechanism may be related to CC chemokine receptor 2 (CCR2) mediated action. CYTL1 also increased the release of reactive oxygen species (ROS) in LPS-stimulated neutrophils. These data indicate that upon inflammatory stimulation, neutrophil-derived CYTL1 may play a crucial role in the activation of neutrophils during pathogenic infections.

G9a regulates tumorigenicity and stemness through genome-wide DNA methylation reprogramming in non-small cell lung cancer

Background

Eukaryotic histone methyltransferases 2 (EHMT2 or G9A) has been regarded as a potential target for non-small cell lung cancer (NSCLC) therapy. This study investigated the regulatory roles of G9A in tumorigenesis and stemness in NSCLC. We isolated and enriched tumor-initiating cells (TIC) from surgically resected NSCLC tissues by FACS and sphere formation assays. We then knocked down G9A using shRNA and carried out genome-wide 850K methylation array and RNA sequencing analyses. We carried out in vivo tumorigenecity asssay using mice xenografts and examined G9A interactions with its novel target using chromatin Immunoprecipitation (ChIP).

Results

We identified 67 genes hypomethylated and 143 genes upregulated following G9A knockdown of which 43 genes were both hypomethylated and upregulated. We selected six genes (CDYL2, DPP4, SP5, FOXP1, STAMBPL1, and ROBO1) for validation. In addition, G9A expression was higher in TICs and targeting G9a by shRNA knockdown or by selective inhibitor UNC0642 significantly inhibited the expression of cancer stem cell markers and sphere forming capacity, in vitro proliferation, and in vivo growth. Further, transient overexpression of FOXP1, a protein may promote normal stem cell differentiation, in TICs resulted in downregulation of stem cell markers and sphere forming capacity and cell proliferation in vitro indicating that the genes we identified are directly regulated by G9A through aberrant DNA methylation and subsequent expression. Similarly, ChIP assay has shown that G9a interacts with its target genes through H3K9me2 and downregulation of H3K9me2 following G9a knockdown disrupts its interaction with its target genes.

Conclusions

These data suggest that G9A is involved in lung cancer stemness through epigenetic mechanisms of maintaining DNA methylation of multiple lung cancer stem cell genes and their expression. Further, targeting G9A or its downstream genes could be a novel therapeutic approach in treating NSCLC patients.

Stratifying patients using fast multiple kernel learning framework: case studies of Alzheimer's disease and cancers

BACKGROUND

Predictive patient stratification is greatly emerging, because it allows us to prospectively identify which patients will benefit from what interventions before their condition worsens. In the biomedical research, a number of stratification methods have been successfully applied and have assisted treatment process. Because of heterogeneity and complexity of medical data, it is very challenging to integrate them and make use of them in practical clinic. There are two major challenges of data integration. Firstly, since the biomedical data has a high number of dimensions, combining multiple data leads to the hard problem of vast dimensional space handling. The computation is enormously complex and time-consuming. Secondly, the disparity of different data types causes another critical problem in machine learning for biomedical data. It has a great need to develop an efficient machine learning framework to handle the challenges.

METHODS

In this paper, we propose a fast-multiple kernel learning framework, referred to as fMKL-DR, that optimise equations to calculate matrix chain multiplication and reduce dimensions in data space. We applied our framework to two case studies, Alzheimer's disease (AD) patient stratification and cancer patient stratification. We performed several comparative evaluations on various biomedical datasets.

RESULTS

In the case study of AD patients, we enhanced significantly the multiple-ROIs approach based on MRI image data. The method could successfully classify not only AD patients and non-AD patients but also different phases of AD patients with AUC close to 1. In the case study of cancer patients, the framework was applied to six types of cancers, i.e., glioblastoma multiforme cancer, ovarian cancer, lung cancer, breast cancer, kidney cancer, and liver cancer. We efficiently integrated gene expression, miRNA expression, and DNA methylation. The results showed that the classification model basing on integrated datasets was much more accurate than classification model basing on the single data type.

CONCLUSIONS

The results demonstrated that the fMKL-DR remarkably improves computational cost and accuracy for both AD patient and cancer patient stratification. We optimised the data integration, dimension reduction, and kernel fusion. Our framework has great potential for mining large-scale cohort data and aiding personalised prevention.

The roles of MASPIN expression and subcellular localization in non-small cell lung cancer

Accumulating studies have confirmed that mammary serine protease inhibitor (MASPIN) plays an essential role in non-small cell lung cancer (NSCLC). However, results are still controversial or inconsistent. In the present study, we attempted to identify the clinical significance of MASPIN and its potential molecular roles in NSCLC. The correlation of MASPIN with prognosis and clinicopathological characteristics was assessed by meta-analysis. Additionally, the potential molecular mechanisms of MASPIN in NSCLC was also investigated through several online databases. A total of 2220 NSCLC patients from 12 high quality studies were included and the results indicated that up-regulated MASPIN nucleus and cytoplasm expression was associated with poor overall survival (OS) (hazard ratio (HR) = 1.43, 95% confidence interval (CI) = 1.01-2.04, P<0.05), elevated MASPIN cytoplasm expression was associated with poor OS (HR = 1.45, 95% CI = 1.01-2.07, P<0.05), disease-free survival (DFS) (HR = 1.95, 95% CI = 1.31-2.88, P=0.001), and disease-specific survival (DSS) (HR = 2.17, 95% CI = 1.18-3.99, P=0.013). MASPIN both nucleus and cytoplasm location were associated with clinicopathological characteristics. Bioinformatics analysis validated the above results and suggested that human serpin family B member 5 (SERPINB5) hypomethylated levels were negatively correlated with its mRNA expression. Bioinformatics analysis also revealed the 85 most frequently altered neighboring genes of SERPINB5, and gene ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analysis revealed 20 GO terms and 3 KEGG pathways with statistical significance. MASPIN had a statistically negative correlation with NSCLC prognosis, functioning as an oncoprotein by hypomethylation and influencing specific pathways involving the 85 genes identified herein. MASPIN might be a promising prognostic signature in NSCLC.

Promoter hypermethylation of cysteine dioxygenase type 1 in patients with non-small cell lung cancer

In the present study, promoter hypermethylation of cysteine dioxygenase type 1 (CDO1) was evaluated in non-small cell lung cancer (NSCLC) tissues to assess the value of CDO1 as a novel biomarker to improve the diagnosis of NSCLC. Tumor tissue samples and corresponding normal lung tissue samples from 42 patients with NSCLC were obtained at the Department of Thoracic Surgery, The Second Xiangya Hospital (Changsha, China). Conventional methylation-specific PCR (cMSP) and methylation-on-beads followed by quantitative methylation-specific PCR (MOB-qMSP) were used to analyze the tumor and normal lung tissue samples. Using these two methods, promoter DNA hypermethylation of the CDO1 gene was detected in 59.4 and 71.0% of tumor tissues of patients with NSCLC and in 9.4 and 0% of normal lung tissue, respectively. Compared with the rate of methylation in the well-differentiated NSCLC tissues (15.4 and 55.6%, respectively), the rate of CDO1 gene promoter methylation was higher in the poorly differentiated tissues (89.5 and 92.3%, respectively). Overall, it was demonstrated that the MOB-qMSP method had a higher positive detection rate for CDO1 hypermethylation compared with the cMSP method. In conclusion, CDO1 gene promoter hypermethylation was more frequently observed in NSCLC tissues compared with in normal lung tissues, and a high methylation frequency of the CDO1 gene in biopsy specimens of NSCLC was associated with the degree of differentiation.

A risk score system based on DNA methylation levels and a nomogram survival model for lung squamous cell carcinoma

Lung squamous cell carcinoma (LSCC) is one of the primary types of non-small cell lung carcinoma, and patients with recurrent LSCC usually have a poor prognosis. The present study was conducted to build a risk score (RS) system for LSCC. Methylation data on LSCC (training set) and on head and neck squamous cell carcinoma (validation set 2) were obtained from The Cancer Genome Atlas database, and GSE39279 (validation set 1) was retrieved from the Gene Expression Omnibus database. Differentially methylated protein-coding genes (DMGs)/long non-coding RNAs (DM-lncRNAs) between recurrence-associated samples and nonrecurrence samples were screened out using the limma package, and their correlation analysis was conducted using the cor.test() function. Following identification of the optimal combinations of DMGs or DM-lncRNAs using the penalized package in R, RS systems were built, and the system with optimal performance was selected. Using the rms package, a nomogram survival model was then constructed. For the differentially expressed genes (DEGs) between the high- and low-risk groups, pathway enrichment analysis was performed by Gene Set Enrichment Analysis. There were 335 DMGs and DM-lncRNAs in total. Following screening out of the top 10 genes (aldehyde dehydrogenase 7 family member A1, chromosome 8 open reading frame 48, cytokine-like 1, heat shock protein 90 alpha family class A member 1, isovaleryl-CoA dehydrogenase, phosphodiesterase 3A, PNMA family member 2, SAM domain, SH3 domain and nuclear localization signals 1, thyroid hormone receptor interactor 13 and zinc finger protein 878) and 6 top lncRNAs, RS systems were constructed. According to Kaplan-Meier analysis, the DNA methylation level-based RS system exhibited the best performance. In combination with independent clinical prognostic factors, a nomogram survival model was built and successfully predicted patient survival. Furthermore, 820 DEGs between the high- and low-risk groups were identified, and 3 pathways were identified to be enriched in this gene set. The 10-DMG methylation level-based RS system and the nomogram survival model may be applied for predicting the outcomes of patients with LSCC.

[Application of GLAD-PCR Assay for Study on DNA Methylation in Regulatory Regions of Some Tumor-Suppressor Genes in Lung Cancer]

Hypermethylation of the gene regulatory regions are common for many cancer diseases. In this work we applied GLAD-PCR assay for identificating of the aberrantly methylated RCGY sites in the regulatory regions of some downregulated genes in tissue samples of lung cancer (LC). This list includes EFEMP1, EPHA5, HOXA5, HOXA9, LHX1, MYF6, NID2, OTX1, PAX9, RARB, RASSF1A, RXRG, SIX6, SKOR1 and TERT genes. The results of DNA samples from 40 cancer and 25 normal lung tissues showed a good diagnostic potential of selected RCGY sites in regulatory regions of MYF6, SIX6, RXRG, LHX1, RASSF1A and TERT genes with relatively high sensitivity (80.0 %) and specificity (88.0 %) of LC detection in tumor DNA.

Protein Cytl1: its role in chondrogenesis, cartilage homeostasis, and disease

Cytokine-like protein 1 (Cytl1), also named Protein C17 or C4orf4 is located on human chromosome 4p15-p16 and encodes a polypeptide of 126 amino acid residues that displays characteristics of a secretory protein. Cytl1 is expressed by a sub-population of CD34+ human mononuclear cells from bone marrow and cord blood, and by chondrocytes (cartilage-forming cells). In this review, we explore evidence suggesting that Cytl1 may be involved in the regulation of chondrogenesis, cartilage homeostasis and osteoarthritis progression, accompanied by the modulation of Sox9 and insulin-like growth factor 1 expression. In addition, Cytl1 exhibits chemotactic and pro-angiogenic biological effects. Interestingly, CCR2 (C-C chemokine receptor type 2) has been identified as a likely receptor for Cytl1, which mediates the ERK signalling pathway. Cytl1 also appears to mediate the TGF-beta-Smad signalling pathway, which is hypothetically independent of the CCR2 receptor. More recently, studies have also potentially linked Cytl1 with a variety of conditions including cardiac fibrosis, smoking, alcohol dependence risk, and tumours such as benign prostatic hypertrophy, lung squamous cell carcinoma, neuroblastoma and familial colorectal cancer. Defining the molecular structure of Cytl1 and its role in disease pathogenesis will help us to design therapeutic approaches for Cytl1-associated pathological conditions.

Gene regulation and prognostic indicators of lung squamous cell carcinoma: TCGA-derived miRNA/mRNA sequencing and DNA methylation data

Lung squamous cell carcinoma (LSCC) is a common cancer worldwide, and this study aimed to investigate the key regulatory networks and prognostic indicators of LSCC. MicroRNA (miRNA)/messenger RNA (mRNA) sequencing and DNA methylation data were obtained from the Cancer Genome Atlas. Differentially expressed miRNAs (DEmiRNAs) and genes (DEGs) were identified by the limma package. Then, the transcription factors (TFs) of DEmiRNAs/DEGs, as well as the targets of miRNAs, were predicted by the TFmiR online tool. Using the t test, aberrant methylation was detected in TF binding sites (TFBSs) in promoters. Finally, integrated network and survival analyses were conducted using SPSS software. We obtained 104 DEmiRNAs and 4,491 DEGs, and validated 2,113 DEGs (VDEGs). Then, 103 TFs, 295 TFs, and 14 DEmiRNAs were predicted to target 95 DEmiRNAs, 821 DEGs and 283 DEGs, respectively. After TF-DEmiRNA/DEG and TF-DEmiRNA-DEG networks were constructed (e.g., E2F1-CDC25A, miR29a-RAN, miR326-TBL1XR1), five feedforward loops between ZEB1 and miR-141/200a/200b/200c/429 were found. Furthermore, VDEGs CDC25A, RAN, TBL1XR1 as well as miR-130b and miR-590 were negatively correlated with survival rates. E2F1-CDC25A, miR29a-RAN, miR326-TBL1XR1, and the feedforward loops between ZEB1/ZEB2 and miR-141/200a/200b/200c/429 might participate in LSCC development. Compared with BEAS-2B cells, the SK-MES-1 cells presented a higher expression level of miR-141, miR-200a, miR-200b, miR-200c but a lower expression level of ZEB1. Overexpressed miR-200c significantly attenuated the expression of ZEB1 and ZEB2 and inhibited the proliferation and migration of SK-MES-1 cells (all p < 0.05). In addition, CDC25A, miR-200a, miR-200b, miR-200c, miR-130b, and miR-590 are potential prognostic indicators of LSCC.

Comparing progression molecular mechanisms between lung adenocarcinoma and lung squamous cell carcinoma based on genetic and epigenetic networks: big data mining and genome-wide systems identification

Non-small-cell lung cancer (NSCLC) is the predominant type of lung cancer in the world. Lung adenocarcinoma (LADC) and lung squamous cell carcinoma (LSCC) are subtypes of NSCLC. We usually regard them as different disease due to their unique molecular characteristics, distinct cells of origin and dissimilar clinical response. However, the differences of genetic and epigenetic progression mechanism between LADC and LSCC are complicated to analyze. Therefore, we applied systems biology approaches and big databases mining to construct genetic and epigenetic networks (GENs) with next-generation sequencing data of LADC and LSCC. In order to obtain the real GENs, system identification and system order detection are conducted on gene regulatory networks (GRNs) and protein-protein interaction networks (PPINs) for each stage of LADC and LSCC. The core GENs were extracted via principal network projection (PNP). Based on the ranking of projection values, we got the core pathways in respect of KEGG pathway. Compared with the core pathways, we found significant differences between microenvironments, dysregulations of miRNAs, epigenetic modifications on certain signaling transduction proteins and target genes in each stage of LADC and LSCC. Finally, we proposed six genetic and epigenetic multiple-molecule drugs to target essential biomarkers in each progression stage of LADC and LSCC, respectively.

CYTL1 inhibits tumor metastasis with decreasing STAT3 phosphorylation

CYTL1 is a novel cytokine that was first identified in CD34⁺ hematopoietic cells. We previously prepared recombinant CYTL1 and verified that it chemoattracted human monocytes via the CCR2/ERK pathway. It has been reported that CYTL1 plays contradictory roles in neuroblastoma and lung cancer. We found that the expression level of CYTL1 was notably decreased and it was hypermethylated in various tumors, including breast and lung cancer, by bioinformatics analyses. After validating the expression of CYTL1 in lung cancer, we identified that CYTL1 exerted no obvious effect on tumor cell proliferation but inhibited their migration and invasion, and these effects were accompanied by decreasing STAT3 phosphorylation, using recombinant CYTL1 and CYTL1-overexpressing tumor cell lines. Furthermore, we constructed experimental and spontaneous metastasis models of breast cancer in BALB/c mice and found that CYTL1 significantly inhibited tumor metastasis in vivo. In summary, CYTL1 is a cytokine with tumor-suppressing characteristics that inhibits tumor metastasis and STAT3 phosphorylation in multiple types of tumors.

Genetic and epigenetic alterations of the EGFR and mutually independent association with BRCA1, MGMT, and RASSF1A methylations in Vietnamese lung adenocarcinomas

Genetic and epigenetic alterations importantly contribute to the pathogenesis of lung cancer. In the study, we measured the frequency and distribution of molecular abnormalities of EGFR as well as the aberrant promoter methylations of BRCA1, MGMT, MLH1, and RASSF1A in Vietnamese lung adenocarcinomas. We investigated the association between genetic and epigenetic alteration, and between each abnormality with clinicopathologic parameters. Somatic EGFR mutation that was found in 49/139 (35.3%) lung adenocarcinomas showed a significant association with young age, female gender, and non-smokers. EGFR overexpression was identified in 82 tumors (59.0%) and statistical relationships with EGFR or BRCA1 methylation but not EGFR mutation. In addition, EGFR, BRCA1, MGMT, MLH1, and RASSF1A methylations were found in 33 (23.7%), 41 (29.5%), 46 (33.1%), 28 (20.1%), and 41 (29.5%) cases of a total of 139 lung adenocarcinomas, respectively. The RASSF1A methylation was found to be linked to the smoking habit. Methylations in MGMT and RASSF1A were also found to correlate with metastasis status. Furthermore, the distribution of EGFR mutation and that of BRCA1, MGMT or RASSF1A methylation were significantly exclusive in lung adenocarcinomas. The main finding of our study demonstrate that epigenetic abnormalities might play a critical role for the lung tumorigenesis in patients with smoking history and metastasis, and partly affect the predictive value of EGFR mutations through blocking expression due to promoter EGFR hypermethylation. Mutually exclusive distribution of genetic and epigenetic alterations reflects differently biological characteristics in the etiology of lung adenocarcinomas.

Combinatorial Electrophoresis and Mass Spectrometry-Based Proteomics in Breast Milk for Breast Cancer Biomarker Discovery

Innovations in approaches for early detection and individual risk assessment of different cancers, including breast cancer (BC), are needed to reduce cancer morbidity and associated mortality. The assessment of potential cancer biomarkers in accessible bodily fluids provides a novel approach to identify the risk and/or onset of cancer. Biomarkers are biomolecules, such as proteins, that are indicative of an abnormality or a disease. Human milk is vastly underutilized biospecimen that offers the opportunity to investigate potential protein BC-biomarkers in young, reproductively active women. As a first step, we have examined the entire protein pattern in human milk samples from breastfeeding mothers with cancer, who were diagnosed either before or after milk donation, and from women without cancer, using mass spectrometry (MS)-based proteomics.

Proteomics analysis of human breast milk to assess breast cancer risk

Detection of breast cancer (BC) in young women is challenging because mammography, the most common tool for detecting BC, is not effective on the dense breast tissue characteristic of young women. In addition to the limited means for detecting their BC, young women face a transient increased risk of pregnancy-associated BC. As a consequence, reproductively active women could benefit significantly from a tool that provides them with accurate risk assessment and early detection of BC. One potential method for detection of BC is biochemical monitoring of proteins and other molecules in bodily fluids such as serum, nipple aspirate, ductal lavage, tear, urine, saliva and breast milk. Of all these fluids, only breast milk provides access to a large volume of breast tissue, in the form of exfoliated epithelial cells, and to the local breast environment, in the form of molecules in the milk. Thus, analysis of breast milk is a non-invasive method with significant potential for assessing BC risk. Here we analyzed human breast milk by mass spectrometry (MS)-based proteomics to build a biomarker signature for early detection of BC. Ten milk samples from eight women provided five paired-groups (cancer versus control) for analysis of dysregulatedproteins: two within woman comparisons (milk from a diseased breast versus a healthy breast of the same woman) and three across women comparisons (milk from a woman with cancer versus a woman without cancer). Despite a wide range in the time between milk donation and cancer diagnosis (cancer diagnosis occurred from 1 month before to 24 months after milk donation), the levels of some proteins differed significantly between cancer and control in several of the five comparison groups. These pilot data are supportive of the idea that molecular analysis of breast milk will identify proteins informative for early detection and accurate assessment of BC risk, and warrant further research. Data are available via ProteomeXchange with identifier PXD007066.

Aberrant CpG-methylation affects genes expression predicting survival in lung adenocarcinoma

Lung adenocarcinoma (LUAD) is a common diagnosed disease with high-mortality rate, and its prognostic implications are under discovered. DNA methylation aberrations are not only an important event for dysregulation of gene expression during tumorigenesis but also a revolution in epigenetics by identifying key prognostic biomarkers for multiple cancers. In this study, we analyzed methylation status of 485 578 CpG sites and RNA-seq transcriptomes of 20 532 genes for 1095 LUAD samples in TCGA database. The association between DNA methylation and the prognostic value of the corresponding gene expression was identified as well. In total, ten aberrantly methylated and dysregulated genes (AURKA, BLK, CNTN2, HMGA1, PTTG1, TNS4, DAPK2, MFSD2A, THSD1, and WNT7A) were highlighted which were significantly correlated with overall survival of 492 LUAD patients, which were all reported as tumor-associated genes in other various cancers and worthy of further investigated and might be used as therapeutic targets for LUAD. Together, methylation aberrances regulate gene expression level during tumorigenesis and influence prognosis of LUAD patients. Integrating knowledge of epigenetics and expression of genes can be useful for an in-depth understanding of cancer mechanism and for the eventual purpose of precisely prognostic and therapeutic target verification.

Developing DNA methylation-based diagnostic biomarkers

An emerging paradigm shift for disease diagnosis is to rely on molecular characterization beyond traditional clinical and symptom-based examinations. Although genetic alterations and transcription signature were first introduced as potential biomarkers, clinical implementations of these markers are limited due to low reproducibility and accuracy. Instead, epigenetic changes are considered as an alternative approach to disease diagnosis. Complex epigenetic regulation is required for normal biological functions and it has been shown that distinctive epigenetic disruptions could contribute to disease pathogenesis. Disease-specific epigenetic changes, especially DNA methylation, have been observed, suggesting its potential as disease biomarkers for diagnosis. In addition to specificity, the feasibility of detecting disease-associated methylation marks in the biological specimens collected noninvasively, such as blood samples, has driven the clinical studies to validate disease-specific DNA methylation changes as a diagnostic biomarker. Here, we highlight the advantages of DNA methylation signature for diagnosis in different diseases and discuss the statistical and technical challenges to be overcome before clinical implementation.

Waking up dormant tumor suppressor genes with zinc fingers, TALEs and the CRISPR/dCas9 system

The aberrant epigenetic silencing of tumor suppressor genes (TSGs) plays a major role during carcinogenesis and regaining these dormant functions by engineering of sequence-specific epigenome editing tools offers a unique opportunity for targeted therapies. However, effectively normalizing the expression and regaining tumor suppressive functions of silenced TSGs by artificial transcription factors (ATFs) still remains a major challenge. Herein we describe novel combinatorial strategies for the potent reactivation of two class II TSGs, MASPIN and REPRIMO, in cell lines with varying epigenetic states, using the CRISPR/dCas9 associated system linked to a panel of effector domains (VP64, p300, VPR and SAM complex), as well as with protein-based ATFs, Zinc Fingers and TALEs. We found that co-delivery of multiple effector domains using a combination of CRISPR/dCas9 and TALEs or SAM complex maximized activation in highly methylated promoters. In particular, CRISPR/dCas9 VPR with SAM upregulated MASPIN mRNA (22,145-fold change) in H157 lung cancer cells, with accompanying re-expression of MASPIN protein, which led to a concomitant inhibition of cell proliferation and induction of apoptotic cell death. Consistently, CRISPR/dCas9 VP64 with SAM upregulated REPRIMO (680-fold change), which led to phenotypic reprogramming in AGS gastric cancer cells. Altogether, our results outlined novel sequence-specific, combinatorial epigenome editing approaches to reactivate highly methylated TSGs as a promising therapy for cancer and other diseases.

Differential gene expression profiling of matched primary renal cell carcinoma and metastases reveals upregulation of extracellular matrix genes

Background

The majority of renal cell carcinoma (RCC) studies analyze primary tumors, and the corresponding results are extrapolated to metastatic RCC tumors. However, it is unknown if gene expression profiles from primary RCC tumors differs from patient-matched metastatic tumors. Thus, we sought to identify differentially expressed genes between patient-matched primary and metastatic RCC tumors in order to understand the molecular mechanisms underlying the development of RCC metastases.

Patients and methods

We compared gene expression profiles between patient-matched primary and metastatic RCC tumors using a two-stage design. First, we used Affymetrix microarrays on 15 pairs of primary RCC [14 clear cell RCC (ccRCC), 1 papillary] tumors and patient-matched pulmonary metastases. Second, we used a custom NanoString panel to validate seven candidate genes in an independent cohort of 114 ccRCC patients. Differential gene expression was evaluated using a mixed effect linear model; a random effect denoting patient was included to account for the paired data. Third, The Cancer Genome Atlas (TCGA) data were used to evaluate associations with metastasis-free and overall survival in primary ccRCC tumors.

Results

We identified and validated up regulation of seven genes functionally involved in the formation of the extracellular matrix (ECM): DCN, SLIT2, LUM, LAMA2, ADAMTS12, CEACAM6 and LMO3. In primary ccRCC, CEACAM6 and LUM were significantly associated with metastasis-free and overall survival (P < 0.01).

Conclusions

We evaluated gene expression profiles using the largest set to date, to our knowledge, of patient-matched primary and metastatic ccRCC tumors and identified up regulation of ECM genes in metastases. Our study implicates up regulation of ECM genes as a critical molecular event leading to visceral, bone and soft tissue metastases in ccRCC.

QSEA-modelling of genome-wide DNA methylation from sequencing enrichment experiments

Genome-wide enrichment of methylated DNA followed by sequencing (MeDIP-seq) offers a reasonable compromise between experimental costs and genomic coverage. However, the computational analysis of these experiments is complex, and quantification of the enrichment signals in terms of absolute levels of methylation requires specific transformation. In this work, we present QSEA, Quantitative Sequence Enrichment Analysis, a comprehensive workflow for the modelling and subsequent quantification of MeDIP-seq data. As the central part of the workflow we have developed a Bayesian statistical model that transforms the enrichment read counts to absolute levels of methylation and, thus, enhances interpretability and facilitates comparison with other methylation assays. We suggest several calibration strategies for the critical parameters of the model, either using additional data or fairly general assumptions. By comparing the results with bisulfite sequencing (BS) validation data, we show the improvement of QSEA over existing methods. Additionally, we generated a clinically relevant benchmark data set consisting of methylation enrichment experiments (MeDIP-seq), BS-based validation experiments (Methyl-seq) as well as gene expression experiments (RNA-seq) derived from non-small cell lung cancer patients, and show that the workflow retrieves well-known lung tumour methylation markers that are causative for gene expression changes, demonstrating the applicability of QSEA for clinical studies. QSEA is implemented in R and available from the Bioconductor repository 3.4 (www.bioconductor.org/packages/qsea).

Biochemical and biophysical characterization of cytokine-like protein 1 (CYTL1)

Cytokine-like protein 1 (CYTL1) is a small widely expressed secreted protein lacking significant primary sequence homology to any other known protein. CYTL1 expression appears to be highest in the hematopoietic system and in chondrocytes; however, maintenance of cartilage in mouse models of arthritis is its only reported function in vivo. Despite lacking sequence homology to chemokines, CYTL1 is predicted by computational methods to fold like a chemokine, and has been reported to function as a chemotactic agonist at the chemokine receptor CCR2 in mouse monocyte/macrophages. Nevertheless, since chemokines are defined by structure and chemokine receptors are able to bind many non-chemokine ligands, direct determination of the CYTL1 tertiary structure will ultimately be required to know whether it actually folds as a chemokine and therefore is a chemokine. Towards this goal, we have developed a method for producing functional recombinant human CYTL1 in bacteria, and we provide new evidence about the biophysical and biochemical properties of recombinant CYTL1. Circular dichroism analysis showed that, like chemokines, CYTL1has a higher content of beta-sheet than alpha-helix secondary structure. Furthermore, recombinant CYTL1 promoted calcium flux in chondrocytes. Nevertheless, unlike chemokines, CYTL1 had limited affinity to proteoglycans. Together, these properties further support cytokine-like properties for CYTL1 with some overlap with the chemokines.

Identification of Methylation-Driven, Differentially Expressed STXBP6 as a Novel Biomarker in Lung Adenocarcinoma

DNA methylation is an essential epigenetic marker associated with the silencing of gene expression. Although various genome-wide studies revealed aberrantly methylated gene targets as molecular biomarkers for early detection, the survival rate of lung cancer patients is still poor. In order to identify methylation-driven biomarkers, genome-wide changes in DNA methylation and differential expression in 32 pairs of lung adenocarcinoma and adjacent normal lung tissue in non-smoking women were examined. This concurrent analysis identified 21 negatively correlated probes (r ≤ −0.5), corresponding to 17 genes. Examining the endogenous expression in lung cancer cell lines, five of the genes were found to be significantly down-regulated. Furthermore, in tumor cells alone, 5-aza-2′-deoxycytidine treatment increased the expression levels of STXBP6 in a dose dependent manner and pyrosequencing showed higher percentage of methylation in STXBP6 promoter. Functional analysis revealed that overexpressed STXBP6 in A549 and H1299 cells significantly decreased cell proliferation, colony formation, and migration, and increased apoptosis. Finally, significantly lower survival rates (P < 0.05) were observed when expression levels of STXBP6 were low. Our results provide a basis for the genetic etiology of lung adenocarcinoma by demonstrating the possible role of hypermethylation of STXBP6 in poor clinical outcomes in lung cancer patients.

A prognostic predictor panel with DNA methylation biomarkers for early-stage lung adenocarcinoma in Asian and Caucasian populations

Background

The incidence of lung adenocarcinoma (LUAD) is increasing worldwide with different prognosis even in early-stage patients. We aimed to identify a prognostic panel with multiple DNA methylation biomarkers to predict survival in early-stage LUAD patients of different racial groups.

Methods

The methylation array, pyrosequencing methylation assay, Cox regression and Kaplan-Meier analyses were conducted to build the risk score equations of selected probes in a training cohort of 69 Asian LUAD patients. The risk score model was verified in another cohort of 299 Caucasian LUAD patients in The Cancer Genome Atlas (TCGA) database.

Results

We performed a Cox regression analysis, in which the regression coefficients were obtained for eight probes corresponding to eight genes (AGTRL1, ALDH1A3, BDKRB1, CTSE, EFNA2, NFAM1, SEMA4A and TMEM129). The risk score was derived from sum of each methylated probes multiplied by its corresponding coefficient. Patients with the risk score greater than the median value showed poorer overall survival compared with other patients (p = 0.007). Such a risk score significantly predicted patients showing poor survival in TCGA cohort (p = 0.036). A multivariate analysis was further performed to demonstrate that the eight-probe panel association with poor outcome in early-stage LUAD patients remained significant even after adjusting for different clinical variables including staging parameters (hazard ratio, 2.03; p = 0.039).

Conclusions

We established a proof-of-concept prognostic panel consisting of eight-probe signature to predict survival of early-stage LUAD patients of Asian and Caucasian populations.

Electronic supplementary material

The online version of this article (doi:10.1186/s12929-016-0276-x) contains supplementary material, which is available to authorized users.

miR-196b Is Epigenetically Silenced during the Premalignant Stage of Lung Carcinogenesis

miRNA silencing by promoter hypermethylation may represent a mechanism by which lung cancer develops and progresses, but the miRNAs involved during malignant transformation are unknown. We previously established a model of premalignant lung cancer wherein we treated human bronchial epithelial cells (HBEC) with low doses of tobacco carcinogens. Here, we demonstrate that next-generation sequencing of carcinogen-transformed HBECs treated with the demethylating agent 5-aza-2'deoxycytidine revealed miR-196b and miR-34c-5p to be epigenetic targets. Bisulfite sequencing confirmed dense promoter hypermethylation indicative of silencing in multiple malignant cell lines and primary tumors. Chromatin immunoprecipitation studies further demonstrated an enrichment in repressive histone marks on the miR-196b promoter during HBEC transformation. Restoration of miR-196b expression by transfecting transformed HBECs with specific mimics led to cell-cycle arrest mediated in part through transcriptional regulation of the FOS oncogene, and miR-196b reexpression also significantly reduced the growth of tumor xenografts. Luciferase assays demonstrated that forced expression of miR-196b inhibited the FOS promoter and AP-1 reporter activity. Finally, a case-control study revealed that methylation of miR-196b in sputum was strongly associated with lung cancer (OR = 4.7, P < 0.001). Collectively, these studies highlight miR-196b as a tumor suppressor whose silencing early in lung carcinogenesis may provide a selective growth advantage to premalignant cells. Targeted delivery of miR-196b could therefore serve as a preventive or therapeutic strategy for the management of lung cancer. Cancer Res; 76(16); 4741-51. ©2016 AACR.

Epigenetic Repression of CCDC37 and MAP1B Links Chronic Obstructive Pulmonary Disease to Lung Cancer

INTRODUCTION

Lung cancer and chronic obstructive pulmonary disease (COPD) share environmental risk factors. COPD also increases the risk of lung cancer; however, the molecular mechanisms are unclear.

METHODS

An epigenome-wide association study of lung tumors and cancer-free lung tissue (CFLT) pairs from non-small-cell lung cancer cases with (n = 18) or without (n = 17) COPD was conducted using the HumanMethylation450 beadchip (HM450K). COPD-associated methylation of top-ranked genes was confirmed in a larger sample set, independently validated, and their potential as sputum-based biomarkers was investigated.

RESULTS

Methylation of CCDC37 and MAP1B was more prevalent in lung tumors from COPD than non-COPD cases [54 of 71 (76%) versus 20 of 46 (43%), p = 0.0013] and [48 of 71 (68%) versus 17 of 46 (37%), p = 0.0035], respectively, after adjustment for age, sex, smoking status, and tumor histology. HM450K probes across CCDC37 and MAP1B promoters showed higher methylation in tumors than CFLT with the highest methylation seen in tumors from COPD cases (p < 0.05). These results were independently validated using The Cancer Genome Atlas data. CCDC37 methylation was more prevalent in sputum from COPD than non-COPD smokers (p < 0.005) from two cohorts. CCDC37 and MAP1B expression was dramatically repressed in tumors and CFLT from COPD than non-COPD cases, p less than 0.02.

CONCLUSIONS

The reduced expression of CCDC37 and MAP1B associated with COPD likely predisposes these genes to methylation that in turn, may contribute to lung cancer.

Expressions of SH3BP5, LMO3, and SNAP25 in diffuse large B-cell lymphoma cells and their association with clinical features

Diffuse large B‐cell lymphoma (DLBCL) is clinicopathologically and genetically heterogeneous with variable clinical outcomes. We previously identified signature genes overexpressed in CD5‐positive (CD5⁺) DLBCL, which is a poor prognostic subgroup of DLBCL. To elucidate the clinical significance of the protein expression of the signature genes overexpressed in CD5⁺DLBCL with regard to all DLBCL, not otherwise specified (NOS), 10 genes (SH3BP5,LMO3,SNAP25,SYT5,SV2C,CABP1,FGF1,FGFR2,NEUROD1, and SYN2) were selected and examined immunohistochemically with samples from 28 patients with DLBCL, NOS. Only three protein expressions, SH3BP5, LMO3, and SNAP25, were detected in DLBCL cells and then analyzed further with samples from 187 patients with DLBCL, NOS. The SH3BP5, LMO3, and SNAP25 proteins were expressed in 60% (103/173), 34% (59/175), and 46% (77/168) of DLBCL patients, respectively. These protein expressions were associated with CD5 expression, and only SH3BP5 was frequently expressed in activated B‐cell‐like DLBCL (P = 0.046). Compared to the SH3BP5‐negative group, the SH3BP5⁺ group was correlated with elderly onset (>60 years, P = 0.0096) and advanced‐stage disease (stage III/IV, P = 0.037). The LMO3⁺ group showed a worse performance status (>1, P = 0.0004). The SH3BP5⁺ group and the LMO3⁺ group had significantly worse overall survival than the negative groups (P = 0.030, 0.034; respectively) for the entire group. In a subgroup analysis of patients treated with rituximab‐containing chemotherapy, there was no significant difference between groups. To the best of our knowledge, this is the first report showing the protein expressions of SH3BP5, LMO3, and SNAP25 in DLBCL cells and their clinical significance in patients with DLBCL. The SH3BP5 and LMO3 protein expressions are associated with the baseline clinical characteristics of DLBCL.

Enhancer methylation dynamics contribute to cancer plasticity and patient mortality

During development, enhancers play pivotal roles in regulating gene expression programs; however, their involvement in cancer progression has not been fully characterized. We performed an integrative analysis of DNA methylation, RNA-seq, and small RNA-seq profiles from thousands of patients, including 25 diverse primary malignances and seven body sites of metastatic melanoma. We found that enhancers are consistently the most differentially methylated regions (DMR) as cancer progresses from normal to primary tumors and then to metastases, compared to other genomic features. Remarkably, identification of enhancer DMRs (eDMRs) enabled classification of primary tumors according to physiological organ systems, and in metastasis eDMRs are the most correlated with patient outcome. To further understand the eDMR role in cancer progression, we developed a model to predict genes and microRNAs that are regulated by enhancer and not promotor methylation, which shows high accuracy with chromatin architecture methods and was experimentally validated. Interestingly, among all metastatic melanoma eDMRs, the most correlated with patient survival were eDMRs that "switched" their methylation patterns back and forth between normal, primary, and metastases and target cancer drivers, e.g., KIT We further demonstrated that eDMR target genes were modulated in melanoma by the bone metastasis microenvironment, suggesting that eDMRs respond to microenvironmental cues in metastatic niches. Our findings that aberrant methylation in cancer cells mostly affects enhancers, which contribute to tumor progression and cancer cell plasticity, will facilitate development of epigenetic anticancer approaches.

Aspirin and colorectal cancer: the promise of precision chemoprevention

Aspirin (acetylsalicylic acid) has become one of the most commonly used drugs, given its role as an analgesic, antipyretic and agent for cardiovascular prophylaxis. Several decades of research have provided considerable evidence demonstrating its potential for the prevention of cancer, particularly colorectal cancer. Broader clinical recommendations for aspirin-based chemoprevention strategies have recently been established; however, given the known hazards of long-term aspirin use, larger-scale adoption of an aspirin chemoprevention strategy is likely to require improved identification of individuals for whom the protective benefits outweigh the harms. Such a precision medicine approach may emerge through further clarification of aspirin's mechanism of action.

LARVA: an integrative framework for large-scale analysis of recurrent variants in noncoding annotations

In cancer research, background models for mutation rates have been extensively calibrated in coding regions, leading to the identification of many driver genes, recurrently mutated more than expected. Noncoding regions are also associated with disease; however, background models for them have not been investigated in as much detail. This is partially due to limited noncoding functional annotation. Also, great mutation heterogeneity and potential correlations between neighboring sites give rise to substantial overdispersion in mutation count, resulting in problematic background rate estimation. Here, we address these issues with a new computational framework called LARVA. It integrates variants with a comprehensive set of noncoding functional elements, modeling the mutation counts of the elements with a β-binomial distribution to handle overdispersion. LARVA, moreover, uses regional genomic features such as replication timing to better estimate local mutation rates and mutational hotspots. We demonstrate LARVA's effectiveness on 760 whole-genome tumor sequences, showing that it identifies well-known noncoding drivers, such as mutations in the TERT promoter. Furthermore, LARVA highlights several novel highly mutated regulatory sites that could potentially be noncoding drivers. We make LARVA available as a software tool and release our highly mutated annotations as an online resource (larva.gersteinlab.org).

Promoter Methylation of CDO1 Identifies Clear-Cell Renal Cell Cancer Patients with Poor Survival Outcome

PURPOSE

In this era of molecular diagnostics, prediction of clear-cell renal cell cancer (ccRCC) survival requires optimization, as current prognostic markers fail to determine individual patient outcome. Epigenetic events are promising molecular markers. Promoter CpG island methylation of cysteine dioxygenase type 1 (CDO1), which was identified as prognostic marker for breast cancer, is studied as a potential marker for ccRCC survival.

EXPERIMENTAL DESIGN

We collected primary tissues of 365 ccRCC cases identified within the prospective Netherlands Cohort Study (NLCS). In this population-based series, CDO1 promoter methylation was observed in 124 of 324 (38.3%) patients with successful methylation-specific PCR analysis. Kaplan-Meier curves and Wilcoxon tests were used to evaluate 10-year ccRCC-specific survival. Cox regression analysis was used to obtain crude and multivariate HRs and 95% confidence intervals (CI). The relative prognostic value of multivariate models with and without CDO1 promoter methylation was compared using likelihood-ratio tests.

RESULTS

Patients with CDO1 promoter methylation have a significantly poorer survival than those without (Wilcoxon P = 0.006). Differences in survival were independent of other prognostic factors, including age and sex (HR, 1.66; 95% CI, 1.12-2.45) and TNM stage, tumor size, and Fuhrman grade (HR, 1.89; 95% CI, 1.25-2.85). Multivariate models performed better with than without CDO1 promoter methylation status (likelihood-ratio P = 0.003). Survival curves were validated in an independent series of 280 ccRCC cases from The Cancer Genome Atlas (TCGA; Wilcoxon P < 0.001).

CONCLUSIONS

CDO1 promoter methylation may not substitute common prognostic makers to predict ccRCC survival, but offers additional, relevant prognostic information, indicating that it might be a novel molecular marker to determine ccRCC prognosis.

Distinct differences in global gene expression profiles in non-implanted blastocysts and blastocysts resulting in live birth

Results from animal models points towards the existence of a gene expression profile that is distinguishably different in viable embryos compared with non-viable embryos. Knowledge of human embryo transcripts is however limited, in particular with regard to how gene expression is related to clinical outcome. The purpose of the present study was therefore to determine the global gene expression profiles of human blastocysts. Next Generation Sequencing was used to identify genes that were differentially expressed in non-implanted embryos and embryos resulting in live birth. Three trophectoderm biopsies were obtained from morphologically high quality blastocysts resulting in live birth and three biopsies were obtained from non-implanting blastocysts of a comparable morphology. Total RNA was extracted from all samples followed by complete transcriptome sequencing. Using a set of filtering criteria, we obtained a list of 181 genes that were differentially expressed between trophectoderm biopsies from embryos resulting in either live birth or no implantation (negative hCG), respectively. We found that 37 of the 181 genes displayed significantly differential expression (p<0.05), e.g. EFNB1, CYTL1 and TEX26 and TESK1, MSL1 and EVI5 in trophectoderm biopsies associated with live birth and non-implanting, respectively. Out of the 181 genes, almost 80% (145 genes) were up-regulated in biopsies from un-implanted embryos, whereas only 20% (36 genes) showed an up-regulation in the samples from embryos resulting in live birth. Our findings suggest the presence of molecular differences visually undetectable between implanted and non-implanted embryos, and represent a proof of principle study.

Epigenetic regulation of the novel tumor suppressor cysteine dioxygenase 1 in esophageal squamous cell carcinoma

Esophageal squamous cell carcinoma (ESCC), the most common subtype of esophageal cancer in East Asian countries, is still associated with a poor prognosis because of the high frequency of lymph node metastasis and invasion. In our previous study, we identified a novel methylation gene, cysteine dioxygenase 1 (CDO1) that is involved in the conversion of cysteine to cysteine sulfinate, and plays a key role in taurine biosynthesis. Decreased expression of CDO1 was observed in ESCC cell lines and tumors derived from patient tissues, and CDO1 silencing could be reversed by treatment with 5-aza-2'-deoxycytidine in six ESCC cell lines. Forced expression of CDO1 in three different ESCC cell lines, TE-4, TE-6, and TE-14, significantly decreased tumor cell growth, cell migration, invasion, and the ability of colony formation. Although CDO1 expression was not found to significantly correlate with survival in ESCC patients, our results suggest that methylation-regulated CDO1 may represent a functional tumor suppressor and a potentially valuable diagnostic biomarker for ESCC.

Association of aspirin and NSAID use with risk of colorectal cancer according to genetic variants

IMPORTANCE

Use of aspirin and other nonsteroidal anti-inflammatory drugs (NSAIDs) is associated with lower risk of colorectal cancer.

OBJECTIVE

To identify common genetic markers that may confer differential benefit from aspirin or NSAID chemoprevention, we tested gene × environment interactions between regular use of aspirin and/or NSAIDs and single-nucleotide polymorphisms (SNPs) in relation to risk of colorectal cancer.

DESIGN, SETTING, AND PARTICIPANTS

Case-control study using data from 5 case-control and 5 cohort studies initiated between 1976 and 2003 across the United States, Canada, Australia, and Germany and including colorectal cancer cases (n=8634) and matched controls (n=8553) ascertained between 1976 and 2011. Participants were all of European descent.

EXPOSURES

Genome-wide SNP data and information on regular use of aspirin and/or NSAIDs and other risk factors.

MAIN OUTCOMES AND MEASURES

Colorectal cancer.

RESULTS

Regular use of aspirin and/or NSAIDs was associated with lower risk of colorectal cancer (prevalence, 28% vs 38%; odds ratio [OR], 0.69 [95% CI, 0.64-0.74]; P = 6.2 × 10(-28)) compared with nonregular use. In the conventional logistic regression analysis, the SNP rs2965667 at chromosome 12p12.3 near the MGST1 gene showed a genome-wide significant interaction with aspirin and/or NSAID use (P = 4.6 × 10(-9) for interaction). Aspirin and/or NSAID use was associated with a lower risk of colorectal cancer among individuals with rs2965667-TT genotype (prevalence, 28% vs 38%; OR, 0.66 [95% CI, 0.61-0.70]; P = 7.7 × 10(-33)) but with a higher risk among those with rare (4%) TA or AA genotypes (prevalence, 35% vs 29%; OR, 1.89 [95% CI, 1.27-2.81]; P = .002). In case-only interaction analysis, the SNP rs16973225 at chromosome 15q25.2 near the IL16 gene showed a genome-wide significant interaction with use of aspirin and/or NSAIDs (P = 8.2 × 10(-9) for interaction). Regular use was associated with a lower risk of colorectal cancer among individuals with rs16973225-AA genotype (prevalence, 28% vs 38%; OR, 0.66 [95% CI, 0.62-0.71]; P = 1.9 × 10(-30)) but was not associated with risk of colorectal cancer among those with less common (9%) AC or CC genotypes (prevalence, 36% vs 39%; OR, 0.97 [95% CI, 0.78-1.20]; P = .76).

CONCLUSIONS AND RELEVANCE

In this genome-wide investigation of gene × environment interactions, use of aspirin and/or NSAIDs was associated with lower risk of colorectal cancer, and this association differed according to genetic variation at 2 SNPs at chromosomes 12 and 15. Validation of these findings in additional populations may facilitate targeted colorectal cancer prevention strategies.