Identification of BRAF Inhibitor Resistance-associated lncRNAs Using Genome-scale CRISPR-Cas9 Transcriptional Activation Screening.

Abstract

Approximately 50% of melanomas harbor the BRAF V600E mutation and targeted therapies using BRAF inhibitors improve patient outcomes. Nonetheless, resistance to BRAF inhibitors develops rapidly and remains a challenge in melanoma treatment. In this study, we attempted to isolate long noncoding RNAs (lncRNAs) involved in BRAF inhibitor resistance using a comprehensive screening method. We used a CRISPR-Cas9 synergistic activation mediator (SAM) protein complex in a genome-scale transcriptional activation assay to screen for candidate lncRNA genes related to BRAF inhibitor resistance. Correlation analysis was performed between expression levels of isolated lncRNA genes and IC50 of dabrafenib in a BRAF-mutated melanoma cell line. Next, online databases were used to construct the lncRNA-miRNA-mRNA regulatory network. Finally, we evaluated the significance of the expression levels of these lncRNAs and mRNAs as biomarkers using clinical specimens. We isolated three BRAF inhibitor resistance-associated lncRNA genes, namely SNHG16, NDUFV2-AS1, and LINC01502. We constructed a lncRNA-miRNA-mRNA network of 13 nodes consisting of three lncRNAs, six miRNAs, and four mRNAs. The lncRNAs and target mRNAs from each regulatory axis significantly and positively correlated with each other. Finally, Kaplan-Meier analysis showed that higher expression levels of MITF, which was up-regulated by LINC01502, were significantly associated with worse prognosis in BRAF V600E-mutated melanoma. The identification of these BRAF inhibitor resistance-associated lncRNA genes at the genomic scale and the establishment of the lncRNA-miRNA-mRNA regulatory network provides new insights into the underlying mechanisms of BRAF inhibitor resistance in melanoma.