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.

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.