Identification of thresholds for dichotomizing DNA methylation data

Long non-coding RNA CCAT1 acts as an oncogene to promote radiation resistance in lung adenocarcinoma: an epigenomics-based investigation

Long non-coding RNAs (lncRNAs) appear to be the crucial modulators in various processes and critically influence the oncogenesis. As one of the LncRNAs, LncRNA CCAT1 has been reported to be closely associated with the progression multiple cancers, but its role in modulating the radioresistance of lung adenocarcinoma (LUAD) remains unclear. In our present study, we screened the potential radioresistance related LncRNAs in LUAD based on the data from The Cancer Genome Atlas (TCGA) database. Data suggested that CCAT1 was abundantly expressed in LUAD and CCAT1 was significantly associated with poor prognosis and radioresistance. Moreover, our in vitro experiments showed that radiation treatment could trigger elevated expression of CCAT1 in the human LUAD cell lines. Further loss/gain-of-function investigations indicated that CCAT1 knockdown significantly inhibited cell proliferation, migration and promoted cell apoptosis in NCI-H1299 cells under irradiation, whereas CCAT1 overexpression in A549 cells yield the opposite effects. In summary, we identified the promoting role of CCAT1 in radioresistance of LUAD, which may provide a theoretical basis for radiotherapy sensitization of LUAD.

Vector genome loss and epigenetic modifications mediate decline in transgene expression of AAV5 vectors produced in mammalian and insect cells

Recombinant adeno-associated virus (rAAV) vectors are often produced in HEK293 or Spodoptera frugiperda (Sf)-based cell lines. We compared expression profiles of "oversized" (∼5,000 bp) and "standard-sized" (4,600 bp) rAAV5-human α1-antitrypsin (rAAV5-hA1AT) vectors manufactured in HEK293 or Sf cells and investigated molecular mechanisms mediating expression decline. C57BL/6 mice received 6 × 10¹³ vg/kg of vector, and blood and liver samples were collected through week 57. For all vectors, peak expression (weeks 12-24) declined by 50% to week 57. For Sf- and HEK293-produced oversized vectors, serum hA1AT was initially comparable, but in weeks 12-57, Sf vectors provided significantly higher expression. For HEK293 oversized vectors, liver genomes decreased continuously through week 57 and significantly correlated with A1AT protein. In RNA-sequencing analysis, HEK293 vector-treated mice had significantly higher inflammatory responses in liver at 12 weeks compared with Sf vector- and vehicle-treated mice. Thus, HEK293 vector genome loss led to decreased transgene protein. For Sf-produced vectors, genomes did not decrease from peak expression. Instead, vector genome accessibility significantly decreased from peak to week 57 and correlated with transgene RNA. Vector DNA interactions with active histone marks (H3K27ac/H3K4me3) were significantly reduced from peak to week 57, suggesting that epigenetic regulation impacts transgene expression of Sf-produced vectors.

Methylation of CpG Sites as Biomarkers Predictive of Drug-Specific Patient Survival in Cancer

Background:

Though the development of targeted cancer drugs continues to accelerate, doctors still lack reliable methods for predicting patient response to standard-of-care therapies for most cancers. DNA methylation has been implicated in tumor drug response and is a promising source of predictive biomarkers of drug efficacy, yet the relationship between drug efficacy and DNA methylation remains largely unexplored.

Method:

In this analysis, we performed log-rank survival analyses on patients grouped by cancer and drug exposure to find CpG sites where binary methylation status is associated with differential survival in patients treated with a specific drug but not in patients with the same cancer who were not exposed to that drug. We also clustered these drug-specific CpG sites based on co-methylation among patients to identify broader methylation patterns that may be related to drug efficacy, which we investigated for transcription factor binding site enrichment using gene set enrichment analysis.

Results:

We identified CpG sites that were drug-specific predictors of survival in 38 cancer-drug patient groups across 15 cancers and 20 drugs. These included 11 CpG sites with similar drug-specific survival effects in multiple cancers. We also identified 76 clusters of CpG sites with stronger associations with patient drug response, many of which contained CpG sites in gene promoters containing transcription factor binding sites.

Conclusion:

These findings are promising biomarkers of drug response for a variety of drugs and contribute to our understanding of drug-methylation interactions in cancer. Investigation and validation of these results could lead to the development of targeted co-therapies aimed at manipulating methylation in order to improve efficacy of commonly used therapies and could improve patient survival and quality of life by furthering the effort toward drug response prediction.

Prognostic Implications and Immune Infiltration Characteristics of Chromosomal Instability-Related Dysregulated CeRNA in Lung Adenocarcinoma

An accumulating body of research indicates that long-noncoding RNAs (lncRNAs) regulate the target genes and act as competitive endogenous RNAs (ceRNAs) playing an indispensable role in lung adenocarcinoma (LUAD). LUAD is frequently accompanied by the feature of chromosomal instability (CIN); however, CIN-related ceRNAs have not been investigated yet. We systematically analyzed and integrated CIN-related dysregulated ceRNAs characteristics in LUAD samples for the first time. In TCGA LUAD cohort, CIN in tumor samples was significantly higher than that in those of adjacent, and patients with high CIN risk tended to have worse clinical outcomes. We constructed a double-weighted CIN-related dysregulated ceRNA network, in which edge weight and node weight represented the disorder extent of ceRNA and the correlation of RNA expression level and prognosis, respectively. After module mining and analysis, a potential prognostic biomarker composed of 12 RNAs (8 mRNAs and 4 lncRNAs) named CIN-related dysregulated ceRNAs (CRDC) was obtained. The CRDC risk score had a positive relation with clinical stage and CIN, and patients with high CRDC risk scores exhibited poor prognosis. Moreover, CRDC tended to be an independent risk factor with high robustness to overcome the effect of multicollinearity among other explanatory variables for disease-specific survival (DSS) in TCGA and two GEO cohorts. The result of functional analysis indicated that CRDC was involved in multiple cancer progresses, especially immune-related pathways. The patients with lower CRDC risk had higher B cell, T cell CD4+, T cell CD8+, neutrophil, macrophage, and myeloid dendritic cell infiltration than the patients with higher CRDC risk. Meanwhile, patients with lower CRDC risk could get more benefits from immunological therapy. The results suggested that the CRDC could be a potential prognostic biomarker and an immunotherapy predictor for lung adenocarcinoma.

Distributional changes in gene-specific methylation associated with temperature

Temperature has been related to mean differences in DNA methylation. However, heterogeneity in these associations may exist across the distribution of methylation outcomes. This study examined whether the association between three-week averaged of temperature and methylation differs across quantiles of the methylation distributions in nine candidate genes. We measured gene-specific blood methylation repeatedly in 777 elderly men participating in the Normative Aging Study (1999-2010). We fit quantile regressions for longitudinal data to investigate whether the associations of temperature on methylation (expressed as %5mC) varied across the distribution of the methylation outcomes. We observed heterogeneity in the associations of temperature across percentiles of methylation in F3, TLR-2, CRAT, iNOS, and ICAM-1 genes. For instance, an increase in three-week temperature exposure was associated with a longer left-tail of the F3 methylation distribution. A 5°C increase in temperature was associated with a 0.15%5mC (95% confidence interval (CI): -0.27,-0.04) decrease on the 20th quantile of F3 methylation, but was not significantly related to the 80th quantile of this distribution (Estimate:0.06%5mC, 95%CI: -0.22, 0.35). Individuals with low values of F3, TLR-2, CRAT, and iNOS methylation, as well as a high value of ICAM-1 methylation, may be more susceptible to temperature effects on systemic inflammation.

Genome-scale methylation assessment did not identify prognostic biomarkers in oral tongue carcinomas

Background

DNA methylation profiling of heterogeneous head and neck squamous cell carcinoma (HNSCC) cohorts has been reported to predict patient outcome. We investigated if a prognostic DNA methylation profile could be found in tumour tissue from a single uniform subsite, the oral tongue. The methylation status of 109 comprehensively annotated oral tongue squamous cell carcinoma (OTSCC) formalin-fixed paraffin-embedded (FFPE) samples from a single institution were examined with the Illumina HumanMethylation450K (HM450K) array. Data pre-processing, quality control and analysis were performed using R packages. Probes mapping to SNPs, sex chromosomes and unreliable probes were accounted for prior to downstream analyses. The relationship between methylation and patient survival was examined using both agnostic approaches and feature selection. The cohort was enlarged by incorporation of 331 The Cancer Genome Atlas (TCGA) HNSCC samples, which included 91 TCGA OTSCC samples with HM450K and survival data available.

Results

Given the use of FFPE-derived DNA, we defined different cohorts for separate analyses. Overall, similar results were found between cohorts. With an unsupervised approach, no distinct hypermethylated group of samples was identified and nor was a prognostic methylation profile identified. The use of multiple downstream feature selection approaches, including a linear models for microarray data (LIMMA), centroid feature selection (CFS), and recursive feature elimination (RFE) support vector machines, similarly failed to identify a significant methylation signature informative for patient prognosis or any clinicopathological data available. Furthermore, we were unable to confirm the prognostic methylation profiles or specific prognostic loci reported within the literature for HNSCC.

Conclusions

With genome-scale assessment of DNA methylation using HM450K in one of the largest OTSCC cohorts to date, we were unable to identify a hypermethylated group of tumours or a prognostic methylation signature. This suggests that either DNA methylation in isolation is not likely to be of prognostic value or larger cohorts are required to identify such a biomarker for OTSCC.

Electronic supplementary material

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

Quantitative methodology is critical for assessing DNA methylation and impacts on correlation with patient outcome

BACKGROUND

DNA hypermethylation is reported as a frequent event and prognostic marker in head and neck squamous cell carcinomas (HNSCC). Methylation has been commonly assessed with non-quantitative methodologies, such as methylation-specific PCR (MSP). We investigated previously reported hypermethylated genes with quantitative methodology in oral tongue squamous cell carcinomas (OTSCC).

RESULTS

The methylation status of 12 genes in 115 OTSCC samples was assessed by one or more of three quantitative analyses: methylation sensitive high resolution melting (MS-HRM), sensitive-melting analysis after real time-methylation specific PCR (SMART-MSP), and bisulfite pyrosequencing. In contrast to much of the literature, either no or infrequent locus-specific methylation was identified by MS-HRM for DAPK1, RASSF1A, MGMT, MLH1, APC, CDH1, CDH13, BRCA1, ERCC1, and ATM. The most frequently methylated loci were RUNX3 (18/108 methylated) and ABO (22/107 methylated). Interrogation of the Cancer Genome Atlas (TCGA) HNSCC cohort confirmed the frequency of significant methylation for the loci investigated. Heterogeneous methylation of RUNX3 (18/108) and ABO (22/107) detected by MS-HRM, conferred significantly worse survival (P = 0.01, and P = 0.03). However, following quantification of methylation levels using pyrosequencing, only four tumors had significant quantities (>15%) of RUNX3 methylation which correlated with a worse patient outcome (P <0.001), while the prognostic significance of ABO hypermethylation was lost. RUNX3 methylation was not prognostic for the TCGA cohort (P = 0.76).

CONCLUSIONS

We demonstrated the critical need for quantification of methylation levels and its impact on correlative analyses. In OTSCC, we found little evidence of significant or frequent hypermethylation of many loci reported to be commonly methylated. It is likely that previous reports have overestimated the frequency of significant methylation events as a consequence of the use of non-quantitative methodology.

Connexin 30 expression inhibits growth of human malignant gliomas but protects them against radiation therapy

BACKGROUND

Glioblastomas remain ominous tumors that almost invariably escape treatment. Connexins are a family of transmembrane, gap junction-forming proteins, some members of which were reported to act as tumor suppressors and to modulate cellular metabolism in response to cytotoxic stress.

METHODS

We analyzed the copy number and expression of the connexin (Cx)30 gene gap junction beta-6 (GJB6), as well as of its protein immunoreactivity in several public and proprietary repositories of glioblastomas, and their influence on patient survival. We evaluated the effect of the expression of this gap junction protein on the growth, DNA repair and energy metabolism, and treatment resistance of these tumors.

RESULTS

The GJB6 gene was deleted in 25.8% of 751 analyzed tumors and mutated in 15.8% of 158 tumors. Cx30 immunoreactivity was absent in 28.9% of 145 tumors. Restoration of Cx30 expression in human glioblastoma cells reduced their growth in vitro and as xenografts in the striatum of immunodeficient mice. Cx30 immunoreactivity was, however, found to adversely affect survival in 2 independent retrospective cohorts of glioblastoma patients. Cx30 was found in clonogenic assays to protect glioblastoma cells against radiation-induced mortality and to decrease radiation-induced DNA damage. This radioprotection correlated with a heat shock protein 90-dependent mitochondrial translocation of Cx30 following radiation and an improved ATP production following this genotoxic stress.

CONCLUSION

These results underline the complex relationship between potential tumor suppressors and treatment resistance in glioblastomas and single out GJB6/Cx30 as a potential biomarker and target for therapeutic intervention in these tumors.