Super enhancer loci of EGFR regulate EGFR variant 8 through enhancer RNA and strongly associate with survival in HNSCCs.

Epidermal growth factor receptor (EGFR) has been shown to be overexpressed in human cancers due to mutation, amplification, and epigenetic hyperactivity, which leads to deregulated transcriptional mechanism. Among the eight different EGFR isoforms, the mechanism of regulation of full-length variant 1 is well-known, no studies have examined the function & factors regulating the expression of variant 8. This study aimed to understand the function of EGFR super-enhancer loci and its associated transcription factors regulating the expression of EGFR variant 8. Our study shows that overexpression of variant 8 and its transcription was more prevalent than variant 1 in many cancers and positively correlated with the EGFR-AS1 expression in oral cancer and HNSCC. Notably, individuals overexpressing variant 8 showed shorter overall survival and had a greater connection with other clinical traits than patients with overexpression of variant 1. In this study, TCGA enhancer RNA profiling on the constituent enhancer (CE1 and CE2) region revealed that the multiple enhancer RNAs formed from CE2 by employing CE1 as a promoter. Our bioinformatic analysis further supports the enrichment of enhancer RNA specific chromatin marks H3K27ac, H3K4me1, POL2 and H2AZ on CE2. GeneHancer and 3D chromatin capture analysis showed clustered interactions between CE1, CE2 loci and this interaction may regulates expression of both EGFR-eRNA and variant 8. Moreover, increased expression of SNAI2 and its close relationship to EGFR-AS1 and variant 8 suggest that SNAI2 could regulates variant 8 overexpression by building a MegaTrans complex with both EGFR-eRNA and EGFR-AS1. Our findings show that EGFR variant 8 and its transcriptional regulation & chromatin modification by eRNAs may provide a rationale for targeting RNA splicing in combination with targeted EGFR therapies in cancer.

Linc-ROR genetic variants are associated with the advanced disease in oral squamous cell carcinoma.

OBJECTIVE: The objective of this study is to identify the association between linc-ROR genetic variants and oral squamous cell carcinoma tumorigenesis. DESIGN: Four genetic variants of linc-ROR (rs6420545, rs4801078, rs1942348, and rs9636089) were analyzed in 178 OSCCs and 191 controls of the South Indian population by PCR amplification followed by restriction digestion. In addition, we examined whether these variants alter linc-ROR expression levels and the progression of OSCC. RESULTS: The frequency of linc-ROR rs6420545 and rs4801078 genotypes were significantly associated with advanced tumor grade (>2) (p = 0.002 and p = 0.048), and nodal metastasis (p = 0.001 and p = 0.019), respectively. We observed a significant association of rs6420545 specifically in the over-dominant model [OR 1.77 (95%CI; 1.17-2.68); p = 0.006] and rs9636089 in dominant model [OR 2.17 (95%CI; 1.06 - 4.46); p = 0.03], and allelic model [OR 2.26 (95%CI; 1.13 - 4.53) p = 0.02], respectively. Further, significant upregulation of linc-ROR (p = 0.005) was observed in our cohort, consistent with the HNSCC TCGA dataset (p < 0.0001). CONCLUSIONS: Our findings suggest that the linc-ROR genetic variants could contribute to the metastasis and progression mainly in the late event of tumorigenesis of OSCCs and these variants could be useful in the precision therapeutic management of this cancer particularly in prognosis.

High incidence of PI3K pathway gene mutations in South Indian cervical cancers.

INTRODUCTION: Cervical cancer is the second most common cancer in India. The phosphatidylinositol-3 kinase (PI3K) signaling is one of the most commonly activated pathways in cancer and comprises key molecules commonly targeted in cancer therapy. This study analyzed six PI3K pathway gene mutations. METHODS: We carried out targeted next-generation sequencing of six PI3K pathway genes (PIK3CA, PIK3R1, PTEN, AKT1, TSC2, and mTOR) in a total of 93 South Indian cervical cancer samples and confirmed them by sanger sequencing. RESULTS: The PI3K pathway gene mutations were observed in 54.8% (51/93) of the tumors and PIK3CA was the most mutated (34.4%, 32/93), followed by TSC2 (18.3%, 17/93), and PIK3R1 (14%, 13/93). The PIK3CA hotspot mutations E542K and E545K observed in this study were likely to disrupt the p110α-p85α interaction that could result in the PI3K pathway activation. We also found a few novel mutations in PIK3R1, PTEN, AKT1, TSC2, and mTOR genes while some of the tumors harbored multiple mutations in the genes of the PI3K pathway. The majority of the tumors were positive for high-risk HPV16/18 (60.7%). CONCLUSIONS: The high incidence of the PI3K pathway gene mutations observed in this study could be exploited for the therapeutic management of cervical cancers.

Genetic variant rs10251977 (G>A) in EGFR-AS1 modulates the expression of EGFR isoforms A and D.

Tyrosine kinase inhibitor is an effective chemo-therapeutic drug against tumors with deregulated EGFR pathway. Recently, a genetic variant rs10251977 (G>A) in exon 20 of EGFR reported to act as a prognostic marker for HNSCC. Genotyping of this polymorphism in oral cancer patients showed a similar frequency in cases and controls. EGFR-AS1 expressed significantly high level in tumors and EGFR-A isoform expression showed significant positive correlation (r = 0.6464, p < 0.0001) with reference to EGFR-AS1 expression levels, consistent with larger TCGA HNSCC tumor dataset. Our bioinformatic analysis showed enrichment of alternative splicing marks H3K36me3 and presence of intronic polyA sites spanning around exon 15a and 15b of EGFR facilitates skipping of exon 15b, thereby promoting the splicing of EGFR-A isoform. In addition, high level expression of PTBP1 and its binding site in EGFR and EGFR-AS1 enhances the expression of EGFR-A isoform (r = 0.7404, p < 0.0001) suggesting that EGFR-AS1 expression modulates the EGFR-A and D isoforms through alternative splicing. In addition, this polymorphism creates a binding site for miR-891b in EGFR-AS1 and may negatively regulate the EGFR-A. Collectively, our results suggested the presence of genetic variant in EGFR-AS1 modulates the expression of EGFR-D and A isoforms.

LncRNA OIP5-AS1 is overexpressed in undifferentiated oral tumors and integrated analysis identifies as a downstream effector of stemness-associated transcription factors.

Long non-coding RNAs (lncRNAs) play an important role in the regulation of key cellular processes in early development and cancer. LncRNA Oip5-as1 facilitates stem cell self-renewal in mouse by sponging mmu-miR-7 and modulating NANOG level, yet its role in cancer is less understood. We analyzed OIP5-AS1 expression in oral tumors and in TCGA datasets. We observed overexpression of OIP5-AS1 in oral tumors (P < 0.001) and in tumors of epithelial origin from TCGA. OIP5-AS1 expression was strongly associated with undifferentiated tumors (P = 0.0038). In silico analysis showed miR-7 binding site is conserved in mouse and human OIP5-AS1. However, human NANOG 3'-UTR lost the binding site for hsa-miR-7a-3. Therefore, we screened for other miRNAs that can be sponged by OIP5-AS1 and identified six potential miRNAs and their downstream target genes. Expression analysis showed downregulation of miRNAs and upregulation of downstream target genes, particularly in undifferentiated tumors with high-level of OIP5-AS1 suggesting OIP5-AS1 could post-transcriptionally modulate the downstream target genes. Further, systematic epigenomic analysis of OIP5-AS1 promoter revealed binding motifs for MYC, NANOG and KLF4 suggesting that OIP5-AS1 could be transactivated by stemness-associated transcription factors in cancer. OIP5-AS1 overexpression in undifferentiated oral tumors may be suggestive of enhanced cancer stemness, and consequently, poor clinical outcome.