Background: The BRAF protein kinase is widely studied as a cancer driver and therapeutic target. However, the regulation of its expression is not completely understood.
Results: Taking advantage of the RNA-seq data of more than 4800 patients belonging to 9 different cancer types, we show that BRAF mRNA exists as a pool of 3 isoforms (reference BRAF, BRAF-X1, and BRAF-X2) that differ in the last part of their coding sequences, as well as in the length (BRAF-ref: 76 nt; BRAF-X1 and BRAF-X2: up to 7 kb) and in the sequence of their 3'UTRs. The expression levels of BRAF-ref and BRAF-X1/X2 are inversely correlated, while the most prevalent among the three isoforms varies from cancer type to cancer type. In melanoma cells, the X1 isoform is expressed at the highest level in both therapy-naïve cells and cells with acquired resistance to vemurafenib driven by BRAF gene amplification or expression of the Δ[3-10] splicing variant. In addition to the BRAF-ref protein, the BRAF-X1 protein (the full length as well as the Δ[3-10] variant) is also translated. The expression levels of the BRAF-ref and BRAF-X1 proteins are similar, and together they account for BRAF functional activities. In contrast, the endogenous BRAF-X2 protein is hard to detect because the C-terminal domain is selectively recognized by the ubiquitin-proteasome pathway and targeted for degradation.
Conclusions: By shedding light on the repertoire of BRAF mRNA and protein variants, and on the complex regulation of their expression, our work paves the way to a deeper understanding of a crucially important player in human cancer and to a more informed development of new therapeutic strategies.
Keywords: BRAF; Exon-spanning reads; Melanoma; Protein variants; RNA-sequencing; Transcript variants.