BackgroundDespite a high incidence of colorectal carcinoma, data regarding genetic aberrations in colorectal carcinoma (CRC) patients in Pakistan is scarce. This study aimed to determine the frequency of BRAFV600E mutations in colorectal carcinoma tissue in the Pakistani population and to associate BRAFV600E expression with CD133, a marker of colorectal stem cells, and CDX2 marker of differentiation.MethodsSanger Sequencing of exon 15 (426 bp) including the hotspot V600E was performed on formalin-fixed-paraffin-embedded (FFPE) CRC tissue samples of 115 patients. The samples were subjected to immunohistochemistry (IHC) to assess the expression of BRAFV600E, CDX2, and CD133. Additionally, homology modelling and docking were performed to investigate novel deletions revealed in sequencing.ResultsTwenty-four (20.8%) BRAF variants were identified in the coding region, with V600E mutations detected in 14 (12.2% )cases (GenBank: PP003258.1; Pop Set: 2678087296). Moreover, a wide spectrum of novel non-V600E mutations (8.6%) were identified, including deletions and missense variations. In-silico analysis revealed that due to large deletions in the coding region of three samples, the affinity of the anti-BRAF drugs (Encorafenib and Vemurafenib) for the active site decreased in comparison to the wild type. The IHC analysis showed that BRAFV600E expression was significantly associated with CD133 expression (χ2(1, n=115) = 26.351; p = < 0.001) and with CDX2 expression (χ2(1, n=115) = 14.88; p = 0.001). Multivariate analysis using binary logistic regression revealed association of BRAFV600E mutations with age (OR = 1.123; CI = 1.024–1.232; p = 0.014), gender (OR = 0.071; CI = 0.006–0.831; p = 0.035), grade (0.007; CI = 0-0.644) and CD133 expression (OR = 65.649; CI = 2.153-2001.556; p = 0.016).ConclusionThe present study demonstrates a notably high V600E frequency (12.2%) in comparison to global reported data, which ranges from 0.4 to 18%. This finding reflects the importance of upfront BRAF testing of the genetically distinct population of Pakistan. Previously unreported mutations identified in the sample may be of clinical significance and warrant further investigation. The concomitant high expression and significant association between CD133 and BRAFV600E represent vital actionable genes that may be targeted together to improve CRC patient management.