Integrative analyses of N6-methyladenosine-associated single-nucleotide polymorphisms (m6A-SNPs) identify tumor suppressor gene AK9 in lung cancer.

N6 -methyladenosine (m6 A) modification has been identified as one of the most important epigenetic regulation mechanisms in the development of human cancers. However, the association between m6 A-associated single-nucleotide polymorphisms (m6 A-SNPs) and lung cancer risk remains largely unknown. Here, we identified m6 A-SNPs and examined the association of these m6 A-SNPs with lung cancer risk in 13,793 lung cancer cases and 14,027 controls. In silico functional annotation was used to identify causal m6 A-SNPs and target genes. Furthermore, methylated RNA immunoprecipitation and quantitative real-time polymerase chain reaction (MeRIP-qPCR) assay was performed to assess the m6 A modification level of different genotypes of the causal SNP. In vitro assays were performed to validate the potential role of the target gene in lung cancer. A total of 8794 m6 A-SNPs were detected, among which 397 SNPs in nine susceptibility loci were associated with lung cancer risk, including six novel loci. Bioinformatics analyses indicated that rs1321328 in 6q21 was located around the m6 A modification site of AK9 and significantly reduced AK9 expression (ß = -0.15, p = 2.78 × 10-8 ). Moreover, AK9 was significantly downregulated in lung cancer tissues than that in adjacent normal tissues of samples from the Cancer Genome Atlas and Nanjing Lung Cancer Cohort. MeRIP-qPCR assay suggested that C allele of rs1321328 could significantly decrease the m6 A modification level of AK9 compared with G allele. In vitro assays verified the tumor-suppressing role of AK9 in lung cancer. These findings shed light on the pathogenic mechanism of lung cancer susceptibility loci linked with m6 A modification.

The association between METS-IR, an indirect index for insulin resistance, and lung cancer risk.

BACKGROUND: Insulin resistance has been reported to increase the risk of breast, prostate and colorectal cancer. However, the role of insulin resistance and its interaction with genetic risk in the development of lung cancer remains controversial. Therefore, we aimed to explore the association between a novel metabolic score for insulin resistance (METS-IR) and lung cancer risk. METHODS: A total of 395 304 participants without previous cancer at baseline were included. The Cox proportional hazards regression model was performed to investigate the association between METS-IR and lung cancer risk. In addition, a Mendelian randomization analysis was also performed to explore the causal relationship. The joint effects and additive interactions between METS-IR and polygenetic risk score (PRS) of lung cancer were also investigated. RESULTS: During a median follow-up of 11.03 years (Inter-quartile range (IQR): 10.30-11.73), a total of 3161 incident lung cancer cases were diagnosed in 395 304 participants. There was a significant association between METS-IR and lung cancer risk, with an HR of 1.28 (95% CI: 1.17-1.41). Based on the Mendelian randomization analysis, however, no causal associations were observed. We observed a joint effect but no interaction between METS-IR and genetic risk. The lung cancer incidence was estimated to be 100.42 (95% CI: 91.45-109.38) per 100 000 person-year for participants with a high METS-IR and PRS, while only 42.76 (95% CI: 36.94-48.59) with low METS-IR and PRS. CONCLUSIONS: High METS-IR was significantly associated with an increased risk of lung cancer. Keeping a low level of METS-IR might help reduce the long-term incident risk of lung cancer.

Genetic variants associated with expression of TCF19 contribute to the risk of head and neck cancer in Chinese population.

BACKGROUND: Squamous cell carcinoma of the head and neck (SCCHN) is one of the most common cancers worldwide and includes cancers arising from the oral cavity, pharynx and larynx. Genome-wide association studies have found several genetic variants related to the risk of SCCHN; however, they could only explain a small fraction of the heritability. Thus, more susceptibility loci associated with SCCHN need to be identified. METHODS: An association study was conducted by genotyping 555 patients with SCCHN and 1367 controls in a Chinese population. Single-variant association analysis was conducted on 63 373 SNPs, and the promising variants were then confirmed by a two-stage validation with 1875 SCCHN cases and 4637 controls. Bioinformatics analysis and functional assays were applied to uncover the potential pathogenic mechanism of the promising variants and genes associated with SCCHN. RESULTS: We first identified three novel genetic variants significantly associated with the risk of SCCHN (p=7.45×10-7 for rs2517611 at 6p22.1, p=1.76×10-9 for rs2524182 at 6p21.33 and p=2.17×10-10 for rs3131018 at 6p21.33). Further analysis and biochemical assays showed that rs3094187, which was in a region in high linkage disequilibrium with rs3131018, could modify TCF19 expression by regulating the binding affinity of the transcription factor SREBF1 to the promoter of TCF19. In addition, experiments revealed that the inhibition of TCF19 may affect several important pathways involved in tumourigenesis and attenuate the cell proliferation and migration of SCCHN. CONCLUSION: These findings offer important evidence that functional genetic variants could contribute to development of SCCHN and that TCF19 may function as a putative susceptibility gene for SCCHN.

Functional characterization and clinical significance of super-enhancers in lung adenocarcinoma.

Super-enhancers (SEs) are important transcriptional regulators in tumorigenesis; however, the functional characterization and clinical significance of SEs in lung adenocarcinoma (LUAD) remain unclear. By using H3K27ac ChIP-seq data of two LUAD cell lines and eight lung tissues, we detected 1045 cancer-specific and 5032 normal-specific SEs. Compared to normal-specific SEs, cancer-specific SEs have different regulatory mechanisms where associated target genes were enriched in critical tumor-related pathways and tended to be regulated by transcription factors of Fos Proto-Oncogene, AP-1 Transcription Factor Subunit and Jun Proto-Oncogene, AP-1 Transcription Factor Subunit families. By using expression data of 513 LUAD and 57 adjacent samples from The Cancer Genome Atlas and 80 tumor-normal paired LUAD samples from the Nanjing Lung Cancer Cohort study, we performed differential expression analysis of target genes for SEs and defined 243 crucial SEs. Unsupervised clustering of crucial SEs revealed two subtypes with different levels of genomic aberrations (i.e., mutation and copy number alteration) and clinical outcomes (progression-free interval: p = 0.030; disease-free interval: p = 0.047). In addition, patients with adverse clinical outcomes were more sensitive to three small molecule inhibitors (bortezomib, doxorubicin, and etoposide), and their targets (PSMB5 and TOP2A) also have elevated expression levels among these patients. Taken together, our findings provided a comprehensive characterization of SEs in LUAD and emphasized their clinical significance in LUAD therapy.

Identification of A-to-I RNA editing profiles and their clinical relevance in lung adenocarcinoma.

Adenosine-to-inosine (A-to-I) RNA editing is a widespread posttranscriptional modification that has been shown to play an important role in tumorigenesis. Here, we evaluated a total of 19,316 RNA editing sites in the tissues of 80 lung adenocarcinoma (LUAD) patients from our Nanjing Lung Cancer Cohort (NJLCC) and 486 LUAD patients from the TCGA database. The global RNA editing level was significantly increased in tumor tissues and was highly heterogeneous across patients. The high RNA editing level in tumors was attributed to both RNA (ADAR1 expression) and DNA alterations (mutation load). Consensus clustering on RNA editing sites revealed a new molecular subtype (EC3) that was associated with the poorest prognosis of LUAD patients. Importantly, the new classification was independent of classic molecular subtypes based on gene expression or DNA methylation. We further proposed a simplified model including eight RNA editing sites to accurately distinguish the EC3 subtype in our patients. The model was further validated in the TCGA dataset and had an area under the curve (AUC) of the receiver operating characteristic curve of 0.93 (95%CI: 0.91-0.95). In addition, we found that LUAD cell lines with the EC3 subtype were sensitive to four chemotherapy drugs. These findings highlighted the importance of RNA editing events in the tumorigenesis of LUAD and provided insight into the application of RNA editing in the molecular subtyping and clinical treatment of cancer.

Comprehensive analysis of circRNA expression pattern and circRNA-miRNA-mRNA network in oral squamous cell carcinoma.

OBJECTIVE: CircRNAs are critical gene modulators in tumor initiation and progression. However, the expression pattern and molecular pathogenesis of circRNAs in oral squamous cell carcinoma (OSCC) are still poorly characterized. METHODS: RNA sequencing with CIRCexplorer2 pipeline was performed to identify circRNAs in 46 tumor-normal paired tissues from OSCC patients. Another set of 48 head and neck squamous cell carcinoma samples from the MiOncoCirc database were utilized as an independent validation. RESULTS: Of the 1276 identified high-confidence circRNAs, 154 were differentially expressed between tumor and normal tissues (log2|Fold Change|≥1 and false discovery rate < 0.05). CircRNAs expression was globally down-regulated in tumors compared to normal tissues (P = 9.44 × 10-14). Correlation analysis demonstrated that the global expression of circRNAs was positively related to tumor infiltrating lymphocyte (P = 1.10 × 10-4) and stromal signature (P = 2.70 × 10-3) whereas negatively associated with cell proliferation markers (P = 4.32 × 10-2). CircRNAs-miRNAs-mRNAs regulatory network revealed 6574 interactions, and the target genes were enriched in extracellular matrix and immune-related pathways. Survival analysis were performed on target genes in immune-related pathways, and 20 genes were significantly associated with the prognostic status of OSCC in The Cancer Genome Atlas cohort. The risk model constructed with above 20 genes was associated with the prognosis status of OSCC (HR = 3.28, P = 5.06 × 10-11), and the result was validated in an independent study (GSE41613) (HR = 2.06, P = 1.73 × 10-2). CONCLUSION: CircRNAs showed a global down-regulation pattern in OSCC tissues, and genes regulated by circRNAs primarily involved in immune and extracellular matrix pathways, which could also affect the OSCC prognosis, indicating that they may serve as potential prognostic biomarkers.