Mutation Status and Immunohistochemical Correlation of KRAS, NRAS, and BRAF in 260 Chinese Colorectal and Gastric Cancers

KRAS, NRAS and BRAF are kinases involved in the RAS-RAF-MAPK signaling pathway and also potential tumor-driven genes. Patients with KRAS/NRAS/BRAF mutations are resistant to anti-EGFR monoclonal antibody therapy. The main purpose of this study is to investigate the mutation status and distribution of KRAS/NRAS/BRAF in Chinese colorectal and gastric cancers, and to explore the histopathological changes and related immunohistochemical marker changes caused by these mutations. The mutation status of KRAS (exons 2, codon 12/13), NRAS (exons 2/3/4, codon 12/13/59/61/117/146) and BRAF (exons 15, codon 600) were detected by amplification refractory mutation system polymerase chain reaction (ARMS-PCR) in 86 colon cancer, 140 rectal cancer and 34 gastric cancer tissues. Then, the frequencies and distribution of KRAS/NRAS/BRAF mutations were described in detail. Furthermore, the relationship between KRAS/NRAS/BRAF mutations and the features of histopathological and related immunohistochemical markers were analyzed. The results showed that KRAS/NRAS/BRAF mutation rates in colon cancer were 44.2, 1.2, and 3.5%; in rectal cancer were 37.1, 4.3, and 0.7%; in gastric cancer were none, none and 2.9%. The mutation rate of KRAS in female (48.8%) is significantly higher than that of male (27.8%), and the mutation rate increased with the higher degree of differentiation. Additionally, the mutation rate of BRAF detected by ARMS-PCR (1.77%) was significantly lower than that by immunohistochemistry (4.11%). It also showed that the KRAS/NRAS/BRAF mutation status had a certain relationship with the expression of some immunohistochemical markers. This study provides more data support for clinical research on KRAS/NRAS/BRAF mutation in CRCs or gastric cancers.

Introduction

Colorectal cancer (CRC) and gastric cancer (GC) are common gastrointestinal cancers. The latest epidemiological data shows that the incidence of CRC ranks 4th in malignant tumors, and the mortality rate ranks 2nd; the incidence and mortality of GC both ranks the 16th in malignant tumors (1). Symptoms of CRC and GC are occult, most patients are diagnosed until advanced stages. According to statistics from the National Cancer Institute (https://seer.cancer.gov/statfacts/), the 5-year survival rate is 64.5% for CRC and 31.0% for GC under current treatment conditions (2). In recent years, the advent of anti-epidermal growth factor receptor (EGFR) monoclonal antibodies (MoAbs), such as cetuximab and panitumumab, have contributed to improving the 5-year survival of CRC patients. The benefits of individual genetic profiling for the selection of therapy have been proven in clinical use. However, the incidence and mortality of CRC and GC remain high.

The main function of anti-EGFR MoAbs is to compete with endogenous ligands for binding to EGFR, thereby blocking downstream RAS and MAPK signaling pathways, inhibiting proliferation of cancer cells, and prolonging the survival of patients with advanced cancer (3). KRAS, NRAS and BRAF are kinases on the RAS-RAF-MAPK signaling pathway. If the RAS and RAF genes are mutated, the mutated protein will not be regulated by the upstream EGFR signal and remain in the activated state, continuing to activate the downstream MAPK pathway, leading to cell uncontrolled proliferation and canceration (4). What's worse, mutations in the RAS and RAF genes are independent of each other, and mutations in either of them will lead to activation of the RAS-RAF-MAPK signaling pathway. Meanwhile, KRAS, NRAS and BRAF are potential tumor-driven genes themselves (5). Therefore, only patients with wild-type KRAS, NRAS, and BRAF genes can benefit from anti-EGFR targeted therapy (6–8), while patients with KRAS, NRAS, and BRAF mutations are resistant to anti-EGFR MoAbs therapy (9). Detection of KRAS, NRAS, and BRAF gene mutation status in CRC tissue is a direct and effective method for screening patients for using anti-EGFR targeted drugs (10). The 2017 edition of National Comprehensive Cancer Network (NCCN) recommends that KRAS, NRAS, and BRAF gene mutations should be identified in primary or metastatic tumors of patients with metastatic colorectal cancer, as a basis for predicting whether or not the patient should be treated with anti-EGFR MoAbs (11). Therefore, the detection of multiple genes such as KRAS, NRAS, and BRAF can accurately predict the efficacy of anti-EGFR MoAbs, thereby realizing individualized targeted therapy.

98.5% of the KRAS mutation occurs in codons 12 or 13 of exon 2. The common mutation site of NRAS gene is located in exons 2, 3, and 4 (6). About 81.9% BRAF mutations are located at codon 600 with a conversion of valine to glutamic acid (V600E) (4). Several studies indicated that different mutation types of KRAS, NRAS, and BRAF gene in colorectal cancer tissues have different biological characteristics and lead to different biological changes, which may have different effects on patients. For example, a growing number of studies have shown that patients with a mutation in codon 13 of the KRAS gene may have a poorer prognosis but may significantly benefit from an anti-EGFR targeted therapy (12). However, some other studies have denied this conclusion. Apparently, the effects of different mutations on the biological properties of tumors and the real mechanisms that lead to different outcomes need to be further elucidated. Most of the previous studies focused on the frequencies and prognostic values of KRAS, NRAS, and BRAF mutations, however, there is still a lack of understanding of the histopathological changes and other related protein expressions changes caused by these mutations. At the same time, the KRAS, NRAS, and BRAF gene mutation status and the related histopathological changes in GC tissue is still rarely reported.

In the present study, firstly, we detected the common mutation sites of KRAS, NRAS, and BRAF gene in CRC and GC tissues of 260 patients by amplification-refractory mutation system polymerase chain reaction (ARMS-PCR). Then, we investigated the frequencies and biological characteristics of KRAS, NRAS, and BRAF mutations. Subsequently, we analyzed the relationship between KRAS, NRAS, and BRAF mutations and the changes of histopathological features and related protein expressions. In order to better explain the potential effect of KRAS, NRAS, and BRAF mutations on the efficacy of anti-EGFR MoAbs targeted therapy and the prognosis of CRC and GC patients.

Materials and Methods

Samples

Two hundred sixty patients (including 86 cases of colon cancer, 140 cases of rectal cancer and 34 cases of gastric cancer) are consecutively collected at the Second Hospital of Jilin University between November 2016 and June 2018. All cases were diagnosed as CRC or GC by two independent pathologists. For each sample, the histopathological sections were stained by hematoxylin-eosin (HE) and immunohistochemistry (IHC) for clinical pathological diagnosis. No patients had accepted preoperative adjuvant treatment. The patients' information is listed in Table 1.

TABLE 1

Table 1. Clinicopathological characteristics of the patients.

Ethics Statement

The Ethics Committee of the Second Hospital of Jilin University has a detailed understanding of and approved all experimental protocols in this study. This study conforms to the provisions of the Declaration of Helsinki (as revised in Fortaleza, Brazil, October 2013). We informed all participants according to the consent for the use of their specimens, and written consents were obtained from each patient. All involved methods were carried out in accordance with relevant guidelines and regulations of the Ethics Committee of the Second Hospital of Jilin University.

DNA Extraction and Mutation Detection

Genomic DNA was extracted from surgically excised fresh solid tumor tissues. The TIANamp Genomic DNA Kit (Tiangen Biotech, Beijing, China) were used according to the manufacturer's instructions.

For each sample, mutations of KRAS exons 2 (codon 12 and 13) were detected by Human KRAS Gene Seven Mutation Detection Kit (YZY Medical Science & Technology Co., Ltd., Wuhan, China); mutations of NRAS exons 2 (codon 12 and 13), exons 3 (codon 59 and 61), exons 4 (codon 117 and 146) were detected by Human NRAS Gene Mutation Detection Kit (YZY Medical Science & Technology Co., Ltd., Wuhan, China); mutations of BRAF exons 15 (codon 600) was detected by Human BRAF Gene V600E Mutation Detection Kit (YZY Medical Science & Technology Co., Ltd., Wuhan, China). All operations were strictly performed in accordance with the kit manual. Specifically, diluted 30 ng of total DNA sample to 2 μl, then mixed with 0.2 μl polymerase. The mixture was then added to a tube preloaded with a dual fluorescent probe primer. Real-time quantitative polymerase chain reaction was performed by ABI 7500 Fast Dx (Applied Biosystems Co. Ltd., US) as 37°C for 10 min, 95°C for 5 min, then 40 cycles of 95°C for 15 s and 60°C for 60 s.

IHC Staining

All samples were fixed in 4% neutral formaldehyde solution and embedded in paraffin. Tissue block was sliced into 2 μm and dewaxed, hydrated and antigen repaired by PT link (Dako, Agilent Technologies, USA). Specifically, place the slices in the 65°C preheated repair solution, and incubated for 30 min by heating to 90°C, then cooled to 70°C. Subsequently, the slices were washed by PBS. Primary, secondary antibodies and DAB coloring solution were automated incubated by Autostainer Link 48 (Dako, Agilent Technologies, USA). Specifically, incubated with hydrogen peroxide for 10 min, primary antibody for 30 min, and secondary antibody for 20 min in room temperature. Counterstain with hematoxylin, routine dehydration, transparent, and seal.

Statistical Analysis

Statistical analysis was performed using the SPSS version 21 (SPSS Inc., USA). Categorical variables were compared by the Chi-square or Fisher's exact test; quantitative and ordered variables were compared by the Mann-Whitney test. Normally distributed variables were compared by Student's t-test. The correspondence relationship between mutation status and immunohistochemical marker characteristics were analyzed using Canonical Correlation Analysis and Multiple Correspondence Analysis. P < 0.05 indicate the statistically significantly difference. The Kaplan-Meier (KM) method were used to evaluate the time to diagnosis of survival, recurrence and metastases.

Results

Mutation Frequencies and Distributions

General Situation

The distribution of age between KRAS/NRAS/BRAF mutant type (MT) and wild type (WT) was compared by Student's t-test. Additionally, Chi-square test was amplified to analyze the distribution of different age components (divided into two groups by 60 years old) in KRAS/NRAS/BRAF MT and WT. It can be found that KRAS/NRAS/BRAF mutations were not significantly related to patients' age (Figure 1, Table 2). When analyzing the relationship between gender and KRAS/NRAS/BRAF mutations, it can be found that the mutation rate of KRAS gene in female (48.8%) is significantly higher than that of male (27.8%) (p = 0.001). When analyzing the mutation rates of KRAS, NRAS, and BRAF in different locations, it can be found that KRAS gene mutation rate was significantly different in colon cancers (44.2%), rectal cancers (37.1%) and gastric cancers (0%) (p < 0.001). The mutation distribution is shown in Table 2.

FIGURE 1

Figure 1. Age distribution of KRAS/NRAS/BRAF mutant type and wild type patients.

TABLE 2

Table 2. Histopathological characteristics according to KRAS/NRAS/BRAF mutation status.

Colon Cancer

The average age of colon cancer patients was 63.53 ± 11.24 (Table 1). KRAS/NRAS/BRAF mutations were not significantly related to patients' age analyzed by Student's t-test or gender analyzed by Chi-square test in colon cancer. KRAS mutations were detected in 38 out of 86 (44.2%) colon cancer samples (Table 2, Figure 2A), of which 28 (73.7%) samples had mutations in codon 12 and 10 (26.3%) samples had mutations in codon 13 (Table 3). Among mutations in KRAS codon 12, the main mutant type was 12ASP (34.2%), followed by 12VAL (21.1%) (Table 3, Figure 2A). KRAS mutations occurred in all 7 sites included in this study. In contrast, NRAS had a lower mutation rate. NRAS mutations were detected in 1 out of 86 (1.2%) colon cancer samples (Table 2, Figure 2A). This mutation occurred in exon 3 codon 61 and the mutant type was Q61-Mu (Table 3). BRAF exon 15 codon 600 600Glu mutation was detected in 3 out of 86 (3.5%) colon cancer samples (Table 2, Figure 2A).

FIGURE 2

Figure 2. The mutation distribution of KRAS/NRAS/BRAF in (A) colon cancer, (B) rectal cancer, and (C) gastric cancer.

TABLE 3

Table 3. Frequency and distribution of KRAS/NRAS/BRAF mutations.

Rectal Cancer

The average age of rectal cancer patients was 61.53 ± 10.39 (Table 1). KRAS/NRAS/BRAF mutations were not significantly related to patients' age analyzed by Student's t-test in rectal cancer. The mutation rate of KRAS gene in female (53.7%) is significantly higher than that of male (30.3%) (p = 0.012) analyzed by Chi-square test in rectal cancer. KRAS mutations were detected in 52 out of 140 (37.1%) rectal cancer samples (Table 2, Figure 2B), of which 43 (30.7%) samples had mutations in codon 12 and 9 (6.4%) samples had mutations in codon 13 (Table 3). Among mutations in KRAS codon 12, the main mutant type was 12ASP (16.4%), followed by 12VAL (9.3%) (Table 3, Figure 2B). KRAS mutations occurred in all 7 sites included in this study. It is worth noting that, there is one sample harbored both 12ASP and 12SER mutation. NRAS mutations were detected in 6 out of 140 (4.3%) rectal cancer samples (Table 2, Figure 2B), of which 3 (2.1%) samples had mutations in exon 2 codon 12, 2 (1.4%) samples had mutations in exon 2 codon 13, and 1 (0.7%) sample had mutations in exon 3 codon 61 (Table 3). None mutation was detected in NRAS exon 3 codon 59, exon 4 codon 117 and 146 in this study. BRAF exon 15 codon 600 600Glu mutation was detected in 1 out of 140 (0.7%) rectal cancer samples (Table 2, Figure 2B).

Gastric Cancer

The average age of gastric cancer patients was 64.09 ± 8.74 (Table 1). KRAS/NRAS/BRAF mutations were not significantly related to patients' age analyzed by Student's t-test or gender analyzed by Chi-square test in gastric cancer. Compared with colon cancer and rectal cancer, KRAS and NRAS have a lower mutation rate in gastric cancer. In all 34 gastric cancer samples, neither KRAS nor NRAS mutation was detected, and only 1 sample (2.9%) was detected to have a mutation in BRAF exon 15 codon 600 600Glu (Table 2, Table 3, Figure 2C).

Histopathological Characteristics of Mutations

Histological Type

There were 5 histological types (adenocarcinoma, mucinous adenocarcinoma, low adhesion carcinoma, signet-ring cell carcinoma, squamous cell carcinoma) contained in the 260 investigated tumor samples. The mutation rate of KRAS was significantly different among the five histological types (p = 0.033). Among the five histological types, the mutation rates of KRAS were 34.3% in adenocarcinoma, 60.0% in mucinous adenocarcinoma, 0% in low adhesion carcinoma, 50% in signet-ring cell carcinoma, and 50% in squamous cell carcinoma, respectively. NRAS and BRAF mutations were only detected in 3.0 and 2.2% of adenocarcinomas, but not detected in other histological types. The mutation distribution is shown in Table 2.

Differentiation and Staging

KRAS mutation rate was significantly different in different degrees of tissue differentiation (p = 0.036). The mutation rate was 50.0% in well differentiated cancers, 37.9% in moderate differentiated cancers, and 25.0% in poor differentiated cancers. There was no significant correlation between KRAS/NRAS/BRAF mutations and TNM stage, tumor infiltration depth, and lymph node metastasis. The mutation distribution is shown in Table 2.

IHC Characteristics of Mutations

IHC plays an important role in clinical pathology diagnosis. In the diagnostic process of CRC, BRAF (V600E), PMS2, EGFR, CDX2, CD34, Ki67, P53, MLH1, MSH6, and MSH2 are the most commonly used immunohistochemical markers for pathological typing, differential diagnosis of benign and malignant, and prognosis. Since only one sample was detected to have a mutation in BRAF in all 34 gastric cancer samples, we only investigated the IHC characteristics in colon cancers and rectal cancers.

Interestingly, mutations in the BRAF gene (1.77%) were not completely consistent with the IHC results of BRAF (V600E) (4.11%), but their correspondence is significant (p = 0.004). Moreover, there was a significant difference in the expression of EGFR between the NRAS MT and WT (p = 0.049); and there was a significant difference in the expression of MLH1 between the BRAF MT and WT (p = 0.004) (Table 4). When analyzing colon and rectal cancers separately, the results were similar. There was a significant difference in the expression of BRAF between the BRAF (V600E) MT and WT (p = 0.008) in colon cancers; and there was a significant difference in the expression of EGFR between the NRAS MT and WT (p = 0.021) in rectal cancers (Tables S1, S2). The representative IHC images for the markers were presented in Figure S1.

TABLE 4

Table 4. Immunohistochemistry characteristics according to KRAS/NRAS/BRAF mutation status in CRC.

Of particular concern is the sample with double mutation sites on KRAS, and the immunohistochemistry results are as follows: BRAF (V600E) (–), PMS2 (+), EGFR (+), CDX2 (+), CD34 (–), Ki67 (positive rate 70%), P53 (–), MLH1 (+), MSH6 (+), MSH2 (+).

Correspondence Between Mutations and Immunohistochemical Markers

In order to further explore the correlation between KRAS/NRAS/BRAF mutation status and IHC characteristics in CRC, Canonical Correlation Analysis and Multiple Correspondence Analysis were performed. The Canonical Correlation Analysis results showed that there is a strong correlation between mutation status and IHC characteristics (canonical correlation coefficient is 0.544). Among them, BRAF mutation status had a great influence on the mutation status (canonical correlation coefficient is 0.995); BRAF (V600E) expression level had a great influence on the IHC characteristics (canonical correlation coefficient is 0.711), followed by MSH2 (canonical correlation coefficient is −0.547) and MLH1 (canonical correlation coefficient is −0.500). The canonical correlation analysis structural diagram is shown in Figure 3A. The Multiple Correspondence Analysis results showed that there is a strong correspondence between BRAF mutation status and BRAF (V600E) expression level, which is theoretically obvious (Figure 3B). If the factor of BRAF (V600E) expression level is excluded and re-analyzed, the results showed that the KRAS mutation status had a certain relationship with the expression of P53, Ki67, CDX2, and MSH6; the NRAS mutation status had a certain relationship with the expression of EGFR; and the BRAF mutation status had a certain relationship with the expression of CD34 (Figure 3C).

FIGURE 3

Figure 3. The correlation between KRAS/NRAS/BRAF mutation status and IHC characteristics in CRC. (A) The result of Canonical Correlation Analysis. (B) The result of Multiple Correspondence Analysis. (C) The result of Multiple Correspondence Analysis without the factor of BRAF (V600E) expression level.

Correspondence Between Mutations and Prognosis

Since the sampling period of this study is only 19 months, it is difficult to accurately reflect the correspondence between gene mutation and prognosis. Therefore, we only summarized the prognosis information of patients up to the present stage. According to the type of mutation, the correspondences between mutation and survival, recurrence and metastasis were analyzed separately. No significant difference was found between mutant type and wild type patients. The correspondence between mutations and survival was shown in Figure 4A, the correspondence between mutations and recurrence was shown in Figure 4B, the correspondence between mutations and metastasis was shown in Figure 4C.

FIGURE 4

Figure 4. The correspondence between mutations and prognosis. (A) The correspondence between mutations and survival. (B) The correspondence between mutations and recurrence. (C) The correspondence between mutations and metastasis.

Discussion

In this study, we investigated the mutation status of KRAS (exons 2, codon 12/13), NRAS (exons 2/3/4, codon 12/13/59/61/117/146) and BRAF (exons 15, codon 600) in 260 patients, including 86 cases of colon cancer, 140 cases of rectal cancer and 34 cases of gastric cancer. The results showed that KRAS mutations were detected in 44.2% colon cancer, 37.1% rectal cancer and none in gastric cancer; NRAS mutations were detected in 1.2% colon cancer, 4.3% rectal cancer and none in gastric cancer; BRAF mutations were detected in 3.5% colon cancer, 0.7% rectal cancer, and 2.9% in gastric cancer.

The mutation rate of KRAS gene in female (48.8%) is significantly higher than that of male (27.8%). KRAS gene mutation rate was significantly different in colon cancers (44.2%), rectal cancers (37.1%) and gastric cancers (0%), however, when colon cancer and rectal cancer were compared alone, the difference was not significant.

Compared with the previous studies (Table 5), the mutation rate of KRAS in CRC reported in our study was consistent with that reported by Douillard et al. (6) and Gao et al. (25), who also focused on exons 2, codon 12/13. Guo et al. (13) believed that the mutation rate of KRAS was as high as 52.7%. This conclusion may be because they included more genetic loci into the study (exons 2/3/4, codons 12/13/59/61/117/146/147). Compared to the mutation rate of KRAS, there are fewer studies focusing on the mutation rate of NRAS. The mutation rate of NRAS in CRC obtained in this study was in the midstream position compared to other reports (4, 6, 13, 14, 16). Vaughn et al showed that the mutation rate of NRAS in Americans was only 1.2% (24). We measured a low mutation rate of BRAF in CRC. This result is consistent with the research on Chinese carried by Shen et al. (27). In the study of Mao et al. (20), the BRAF mutation rate reached a staggering 25.4%, much higher than other reports. However, the sample size in that study was only 69 cases, and its representativeness was questionable. By comparison, it can be found that the mutation rate of these three genes is not significantly different between Asians and Europeans.

TABLE 5

Table 5. Studies on mutation status of KRAS, NRAS and BRAF in CRC.

There are many reasons for the different mutation rate results. In addition to the influence of the sample size, different mutation sites included in the study will result in different mutation rates. KRAS mutations occur 98.5% in exon 2 at codons 12 and 13. The common mutation site of NRAS gene is located in exons 2, 3, and 4. About 81.9% BRAF mutations are located at codon 600. Therefore, in this study we focused our attention on the mutations at these sites. Another factor that may affect the outcome of the mutation rate is the detection method. Direct sequencing is the most widely used method for mutation detection (15, 17, 18, 22, 28). It is the gold standard for mutation detection, but it is limited by sensitivity, and only mutations with a mutational heterogeneity more than 10% can be detected. Besides, pyrosequencing (24, 26, 27), high-resolution melting (19, 21), ARMS-PCR(14), cycleave PCR(23) and mass spectrometry genotyping (4) are also used for the detection of mutations. In this study we applied two-color fluorescent probe ARMS-PCR. This method is more sensitive than direct sequencing, as little as 1% of heterogeneous mutations in tumor tissue can be detected. The type of sample to be tested may also have a certain effect on the mutation rate results. In this study, fresh tissue without being fixed by paraformaldehyde was used to avoid DNA damage during the fixation process.

The mutation rate of KRAS was significantly different in different histological types, NRAS and BRAF mutations were only detected in adenocarcinomas. Furthermore, we found that KRAS gene mutation rate was significantly different in different degrees of tissue differentiation, but not significantly associated with TNM stage. KRAS mutation rate increased with the higher degree of differentiation. These results were a little different from those reported by Guo et al. (13), who termed KRAS mutations had no significant correlation with clinicopathological characteristics.

None of the previous studies investigated the association between KRAS/NRAS/BRAF mutations and IHC characteristics in CRCs. In our study, we investigated the association between KRAS/NRAS/BRAF mutations and commonly used immunohistochemical markers, including BRAF (V600E), PMS2, EGFR, CDX2, CD34, Ki67, P53, MLH1, MSH6, and MSH2. We were surprised to find that the mutations in BRAF gene were not completely consistent with the IHC results of BRAF (V600E), the mutation rate of BRAF detected by ARMS-PCR (1.77%) was significantly lower than that by IHC (4.11%). Molecular testing is the gold standard for genetic mutation detection. Although many studies have shown that IHC has a good detection effect on BRAF V600E mutations (29–32), but there will still be a part of false positive results (33, 34). Ballester et al. (35) suggested that the highly sensitive molecular assays remain the gold standard for BRAF mutation analysis in paraffin-embedded lesions. Ehsani et al. (36) used IHC to detect BRAF mutations in metastatic malignant melanoma with a false positive rate of 32%. We suggest that if the purpose of detecting BRAF mutations is to guide anti-EGFR targeted therapy, genetic testing will benefit more patients rather than IHC. Moreover, we also found that there was a significant difference in the expression of EGFR between the NRAS MT and WT; and there was a significant difference in the expression of MLH1 between the BRAF MT and WT. Parsons et al. (37) reported that tumor BRAF mutation, and MLH1 promoter “C region” methylation specifically, are strong predictors of negative MMR mutation status in CRCs. Farchoukh et al. (38) also found that although the presence of the BRAF mutation is indicative of a sporadic cancer, up to 30–50% of colorectal carcinomas with MLH1 promoter hypermethylation will lack a BRAF mutation. Some similar studies have also shown that BRAF mutation is closely linked with MLH1 promoter hypermethylation and DNA mismatch repair (MMR) gene mutations, but its specific mechanism needs further study (39).

Furthermore, we employed Canonical Correlation Analysis and Multiple Correspondence Analysis to further explore the correlation between KRAS/NRAS/BRAF mutation status and IHC characteristics in CRC. The results indicated that the KRAS mutation status had a certain relationship with the expression of P53, Ki67, CDX2, and MSH6; the NRAS mutation status had a certain relationship with the expression of EGFR; and the BRAF mutation status had a certain relationship with the expression of CD34. There are few direct reports of correlation studies between these genes. We hypothesize that since KRAS, NRAS, and BRAF are potential tumor-driven genes, their mutation may have some synergistic or inhibitory effects on the expression of genes such as P53, Ki67, CDX2, MSH6, and CD34.

There are few studies on the mutations rate of KRAS/NRAS/BRAF in gastric cancer. And gastric cancer patients benefit little from anti-EGFR MoAbs targeted therapy. Compared with colon cancer rectal cancer, KRAS/NRAS/BRAF have a lower mutation rate in gastric cancer, furthermore, there is no consistent conclusion on the role of KRAS/NRAS/BRAF mutations in gastric cancer (40–43). In this study we found 1 out of 34 gastric cancer cases with BRAF mutation. No KRAS or NRAS mutation was found in gastric cancer in this study. Here we provide these data for further research by peers.

Since the sampling period of this study is between November 2016 and June 2018, it is difficult to accurately reflect the correspondence between gene mutation and prognosis. Therefore, we only summarized the prognosis information of patients up to the present stage. We will continue to follow this group of patients in subsequent studies to delve into the effects of KRAS/NRAS/BRAF mutation on prognosis.

The limitations of our study are its relatively small sample size, and lack of follow-up time which are important for risk assessment of malignant tumor. NRAS and BRAF mutation frequency was too low to analyze its mutation subgroups. The specific mechanism and clinical significance of the relationship between KRAS/NRAS/BRAF mutation and IHC status still need further experiments to confirm.

In this study, we systematically described and statistically analyzed the frequencies and distributions of KRAS/NRAS/BRAF genetic mutation status and their relationship with IHC in 260 cases with colorectal cancer or gastric cancer through retrospective analysis. Based on the analysis results, we draw the following conclusions: (1) KRAS mutations were detected in 44.2% colon cancer, 37.1% rectal cancer and none in gastric cancer; NRAS mutations were detected in 1.2% colon cancer, 4.3% rectal cancer and none in gastric cancer; BRAF mutations were detected in 3.5% colon cancer, 0.7% rectal cancer and 2.9% in gastric cancer; (2) the mutation rate of KRAS in female (48.8%) was significantly higher than that of male (27.8%); (3) the mutation rate increased with the higher degree of differentiation; (4) the mutation rate of BRAF detected by ARMS-PCR (1.77%) was significantly lower than that by IHC (4.11%); (5) the KRAS/NRAS/BRAF mutation status had a certain relationship with the expression of some common immunohistochemical markers. This study provides more data support for clinical research on KRAS/NRAS/BRAF mutation in CRCs or gastric cancers. In the context of precision medicine, more precise classification of genetic profile should be implemented to enhance the clinical experience. Our study suggested that, combining genetic mutations with immunohistochemical phenotypes could help doctors to formulate cancer treatment strategies more accurately with combining chemotherapy, immunotherapy, and targeted therapy. However, the specific mechanism and clinical significance of the relationship between KRAS/NRAS/BRAF mutation and IHC status still need further experiments to confirm.

Author Contributions

QY conducted experimental operations, sample processing, data analysis, and article writing. YS and HZ conducted experimental operations, sample processing, and data analysis. SH, YL, WL and XW performed sample pretreatment. ZD carried out experimental design and sample pretreatment. TL and GZ conceived and designed the experiments.

Funding

This study was supported by National Natural Science Foundation of China (No.81702195); the Department of Science and Technology of Jilin Province, China (No.20180520111JH, 20180520125JH), the Education Department of Jilin Province, China (No.JJKH20170853KJ, JJKH20180103KJ); the Project of Bethune Youth Foundation of Jilin University, China (No.2015409); the Development and Reform Commission of Jilin Province, China (No.2014G073); the Project of Application Demonstration Center of Precision Medicine for Molecular Diagnosis in Jilin Province (2016-2018, NDRC).

Conflict of Interest Statement

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Acknowledgments

We are grateful to Cao Meixia for her concern and support for this research. It was her proposal and encouragement that prompted me to carry out this research. We are grateful to Jia Meng for carefully screening the most representative pathological and immunohistochemical sections, making our research more complete and rigorous.

Supplementary Material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fonc.2018.00487/full#supplementary-material

Figure S1. The representative HE and IHC images for the markers.

Table S1. Immunohistochemistry characteristics according to KRAS/NRAS/BRAF gene mutation status in colon cancer.

Table S2. Immunohistochemistry characteristics according to KRAS/NRAS/BRAF gene mutation status in rectal cancer.

References

1. Siegel RL, Miller KD, Jemal A. Cancer statistics, 2018. CA Cancer J Clin. (2018) 68:7–30. doi: 10.3322/caac.21442

PubMed Abstract | CrossRef Full Text | Google Scholar

2. Siegel RL, Miller KD, Fedewa SA, Ahnen DJ, Meester RGS, Barzi A, et al. Colorectal cancer statistics, 2017. CA Cancer J Clin. (2017) 67:177–93. doi: 10.3322/caac.21395

PubMed Abstract | CrossRef Full Text | Google Scholar

3. Van Cutsem E, Cervantes A, Adam R, Sobrero A, Van Krieken JH, Aderka D, et al. ESMO consensus guidelines for the management of patients with metastatic colorectal cancer. Ann Oncol. (2016) 27:1386–422. doi: 10.1093/annonc/mdw235

PubMed Abstract | CrossRef Full Text | Google Scholar

4. De Roock W, Claes B, Bernasconi D, De Schutter J, Biesmans B, Fountzilas G, et al. Effects of <em>KRAS, BRAF, NRAS</em>, and <em>PIK3CA</em> mutations on the efficacy of cetuximab plus chemotherapy in chemotherapy-refractory metastatic colorectal cancer: a retrospective consortium analysis. Lancet Oncol. (2010) 11:753–62. doi: 10.1016/S1470-2045(10)70130-3

PubMed Abstract | CrossRef Full Text | Google Scholar

5. Li K, Guo Q, Yang J, Chen H, Hu K, Zhao J, et al. FOXD3 is a tumor suppressor of colon cancer by inhibiting EGFR-Ras-Raf-MEK-ERK signal pathway. Oncotarget (2017) 8:5048–56. doi: 10.18632/oncotarget.13790

PubMed Abstract | CrossRef Full Text | Google Scholar

6. Douillard JY, Oliner KS, Siena S, Tabernero J, Burkes R, Barugel M, et al. Panitumumab-FOLFOX4 treatment and RAS mutations in colorectal cancer. N Engl J Med. (2013) 369:1023–34. doi: 10.1056/NEJMoa1305275

PubMed Abstract | CrossRef Full Text | Google Scholar

7. Popovici V, Budinska E, Bosman FT, Tejpar S, Roth AD, Delorenzi M. Context-dependent interpretation of the prognostic value of BRAF and KRAS mutations in colorectal cancer. BMC Cancer (2013) 13:439. doi: 10.1186/1471-2407-13-439

PubMed Abstract | CrossRef Full Text | Google Scholar

8. Dienstmann R, Salazar R, Tabernero J. Personalizing colon cancer adjuvant therapy: selecting optimal treatments for individual patients. J Clin Oncol. (2015) 33:1787–96. doi: 10.1200/jco.2014.60.0213

PubMed Abstract | CrossRef Full Text | Google Scholar

9. Schmoll HJ, Van Cutsem E, Stein A, Valentini V, Glimelius B, Haustermans K, et al. ESMO consensus guidelines for management of patients with colon and rectal cancer. a personalized approach to clinical decision making. Ann Oncol. (2012) 23:2479–516. doi: 10.1093/annonc/mds236

PubMed Abstract | CrossRef Full Text | Google Scholar

10. Baas JM, Krens LL, Guchelaar HJ, Morreau H, Gelderblom H. Concordance of predictive markers for EGFR inhibitors in primary tumors and metastases in colorectal cancer: a review. Oncologist (2011) 16:1239–49. doi: 10.1634/theoncologist.2011-0024

PubMed Abstract | CrossRef Full Text | Google Scholar

11. Benson AB, 3rd, Venook AP, Cederquist L, Chan E, Chen YJ, Cooper HS, et al. Colon Cancer, Version 1.2017, NCCN clinical practice guidelines in oncology. J Natl Compr Canc Netw. (2017) 15:370–98.

PubMed Abstract | Google Scholar

12. Er TK, Chen CC, Bujanda L, Herreros-Villanueva M. Clinical relevance of KRAS mutations in codon 13: where are we? Cancer Lett. (2014) 343:1–5. doi: 10.1016/j.canlet.2013.09.012

PubMed Abstract | CrossRef Full Text | Google Scholar

13. Guo F, Gong H, Zhao H, Chen J, Zhang Y, Zhang L, et al. Mutation status and prognostic values of KRAS, NRAS, BRAF and PIK3CA in 353 Chinese colorectal cancer patients. Sci Rep. (2018) 8:6076. doi: 10.1038/s41598-018-24306-1

PubMed Abstract | CrossRef Full Text | Google Scholar

14. Zhang J, Zheng J, Yang Y, Lu J, Gao J, Lu T, et al. Molecular spectrum of KRAS, NRAS, BRAF and PIK3CA mutations in Chinese colorectal cancer patients: analysis of 1,110 cases. Sci Rep. (2015) 5:18678. doi: 10.1038/srep18678

PubMed Abstract | CrossRef Full Text | Google Scholar

15. Tong JH, Lung RW, Sin FM, Law PP, Kang W, Chan AW, et al. Characterization of rare transforming KRAS mutations in sporadic colorectal cancer. Cancer Biol Ther. (2014) 15:768–76. doi: 10.4161/cbt.28550

PubMed Abstract | CrossRef Full Text | Google Scholar

16. Shen Y, Wang J, Han X, Yang H, Wang S, Lin D, et al. Effectors of epidermal growth factor receptor pathway: the genetic profiling ofKRAS, BRAF, PIK3CA, NRAS mutations in colorectal cancer characteristics and personalized medicine. PLoS ONE (2013) 8:e81628. doi: 10.1371/journal.pone.0081628

PubMed Abstract | CrossRef Full Text | Google Scholar

17. Pu X, Pan Z, Huang Y, Tian Y, Guo H, Wu L, et al. Comparison of KRAS/BRAF mutations between primary tumors and serum in colorectal cancer: Biological and clinical implications. Oncol Lett. (2013) 5:249–54. doi: 10.3892/ol.2012.963

PubMed Abstract | CrossRef Full Text | Google Scholar

18. Wang J, Yang H, Shen Y, Wang S, Lin D, Ma L, et al. Direct sequencing is a reliable assay with good clinical applicability for KRAS mutation testing in colorectal cancer. Cancer Biomark. (2013) 13:89–97. doi: 10.3233/cbm-130334

PubMed Abstract | CrossRef Full Text | Google Scholar

19. Chang YS, Chang SJ, Yeh KT, Lin TH, Chang JG. RAS, BRAF, and TP53 gene mutations in Taiwanese colorectal cancer patients. Onkologie (2013) 36:719–24. doi: 10.1159/000356814

PubMed Abstract | CrossRef Full Text | Google Scholar

20. Mao C, Zhou J, Yang Z, Huang Y, Wu X, Shen H, et al. KRAS, BRAF and PIK3CA mutations and the loss of PTEN expression in Chinese patients with colorectal cancer. PLoS ONE (2012) 7:e36653. doi: 10.1371/journal.pone.0036653

PubMed Abstract | CrossRef Full Text | Google Scholar

21. Hsieh LL, Er TK, Chen CC, Hsieh JS, Chang JG, Liu TC. Characteristics and prevalence of KRAS, BRAF, and PIK3CA mutations in colorectal cancer by high-resolution melting analysis in Taiwanese population. Clin Chim Acta (2012) 413:1605–11. doi: 10.1016/j.cca.2012.04.029

PubMed Abstract | CrossRef Full Text | Google Scholar

22. Li Z, Chen Y, Wang D, Wang G, He L, Suo J. Detection of KRAS mutations and their associations with clinicopathological features and survival in Chinese colorectal cancer patients. J Int Medical Res. (2012) 40:1589–98. doi: 10.1177/147323001204000439

PubMed Abstract | CrossRef Full Text | Google Scholar

23. Yokota T, Ura T, Shibata N, Takahari D, Shitara K, Nomura M, et al. BRAF mutation is a powerful prognostic factor in advanced and recurrent colorectal cancer. Br J Cancer (2011) 104:856–62. doi: 10.1038/bjc.2011.19

PubMed Abstract | CrossRef Full Text | Google Scholar

24. Vaughn CP, Zobell SD, Furtado LV, Baker CL, Samowitz WS. Frequency of KRAS, BRAF, and NRAS mutations in colorectal cancer. Genes Chromosomes Cancer (2011) 50:307–12. doi: 10.1002/gcc.20854

PubMed Abstract | CrossRef Full Text | Google Scholar

25. Gao J, Wang TT, Yu JW, Li YY, Shen L. Wild-Type KRAS and BRAF could predict response to cetuximab in chinese colorectal cancer patients. Chin J Cancer Res. (2011) 23:271–5. doi: 10.1007/s11670-011-0271-4

PubMed Abstract | CrossRef Full Text | Google Scholar

26. Li HT, Lu YY, An YX, Wang X, Zhao QC. KRAS, BRAF and PIK3CA mutations in human colorectal cancer: relationship with metastatic colorectal cancer. Oncol. Rep. (2011) 25:1691–7. doi: 10.3892/or.2011.1217

PubMed Abstract | CrossRef Full Text | Google Scholar

27. Shen H, Yuan Y, Hu HG, Zhong X, Ye XX, Li MD, et al. Clinical significance of K-ras and BRAF mutations in Chinese colorectal cancer patients. World J Gastroenterol. (2011) 17:809–16. doi: 10.3748/wjg.v17.i6.809

PubMed Abstract | CrossRef Full Text | Google Scholar

28. Liou JM, Wu MS, Shun CT, Chiu HM, Chen MJ, Chen CC, et al. Mutations in BRAF correlate with poor survival of colorectal cancers in Chinese population. Int J Colorectal Dis. (2011) 26:1387–95. doi: 10.1007/s00384-011-1229-1

PubMed Abstract | CrossRef Full Text | Google Scholar

29. Kim JK, Seong CY, Bae IE, Yi JW, Yu HW, Kim SJ, et al. Comparison of immunohistochemistry and direct sequencing methods for identification of the BRAF(V600E) mutation in papillary thyroid carcinoma. Ann Surgical Oncol. (2018) 25:1775–81. doi: 10.1245/s10434-018-6460-3

PubMed Abstract | CrossRef Full Text | Google Scholar

30. Fisher KE, Cohen C, Siddiqui MT, Palma JF, Lipford EH, 3rd, Longshore JW. Accurate detection of BRAF p.V600E mutations in challenging melanoma specimens requires stringent immunohistochemistry scoring criteria or sensitive molecular assays. Hum Pathol. (2014) 45:2281–93. doi: 10.1016/j.humpath.2014.07.014

PubMed Abstract | CrossRef Full Text | Google Scholar

31. Ilie MI, Lassalle S, Long-Mira E, Bonnetaud C, Bordone O, Lespinet V, et al. Diagnostic value of immunohistochemistry for the detection of the BRAF(V600E) mutation in papillary thyroid carcinoma: comparative analysis with three DNA-based assays. Thyroid (2014) 24:858–66. doi: 10.1089/thy.2013.0302

PubMed Abstract | CrossRef Full Text | Google Scholar

32. Long GV, Wilmott JS, Capper D, Preusser M, Zhang YE, Thompson JF, et al. Immunohistochemistry is highly sensitive and specific for the detection of V600E BRAF mutation in melanoma. Am J Surg Pathol. (2013) 37:61–5. doi: 10.1097/PAS.0b013e31826485c0

PubMed Abstract | CrossRef Full Text | Google Scholar

33. Jabbar KJ, Luthra R, Patel KP, Singh RR, Goswami R, Aldape KD, et al. Comparison of next-generation sequencing mutation profiling with BRAF and IDH1 mutation-specific immunohistochemistry. Am J Surg Pathol (2015) 39:454–61. doi: 10.1097/pas.0000000000000325

PubMed Abstract | CrossRef Full Text | Google Scholar

34. Zimmermann AK, Camenisch U, Rechsteiner MP, Bode-Lesniewska B, Rossle M. Value of immunohistochemistry in the detection of BRAF(V600E) mutations in fine-needle aspiration biopsies of papillary thyroid carcinoma. Cancer Cytopathol. (2014) 122:48–58. doi: 10.1002/cncy.21352

PubMed Abstract | CrossRef Full Text | Google Scholar

35. Ballester LY, Cantu MD, Lim KPH, Sarabia SF, Ferguson LS, Renee Webb C, et al. The use of BRAF V600E mutation-specific immunohistochemistry in pediatric Langerhans cell histiocytosis. Hematol Oncol. (2018) 36:307–15. doi: 10.1002/hon.2388

PubMed Abstract | CrossRef Full Text | Google Scholar

36. Ehsani L, Cohen C, Fisher KE, Siddiqui MT. BRAF mutations in metastatic malignant melanoma: comparison of molecular analysis and immunohistochemical expression. Appl Immunohistochem Mole Morphol. (2014) 22:648–51. doi: 10.1097/pai.0000000000000013

PubMed Abstract | CrossRef Full Text | Google Scholar

37. Parsons MT, Buchanan DD, Thompson B, Young JP, Spurdle AB. Correlation of tumour BRAF mutations and MLH1 methylation with germline mismatch repair (MMR) gene mutation status: a literature review assessing utility of tumour features for MMR variant classification. J Med Genet. (2012) 49:151–7. doi: 10.1136/jmedgenet-2011-100714

PubMed Abstract | CrossRef Full Text | Google Scholar

38. Farchoukh L, Kuan SF, Dudley B, Brand R, Nikiforova M, Pai RK. MLH1-deficient colorectal carcinoma with wild-type BRAF and MLH1 promoter hypermethylation harbor KRAS mutations and arise from conventional adenomas. Am J Surg Pathol. (2016) 40:1390–9. doi: 10.1097/pas.0000000000000695

PubMed Abstract | CrossRef Full Text | Google Scholar

39. Adar T, Rodgers LH, Shannon KM, Yoshida M, Ma T, Mattia A, et al. A tailored approach to BRAF and MLH1 methylation testing in a universal screening program for Lynch syndrome. Mod Pathol. (2017) 30:440–7. doi: 10.1038/modpathol.2016.211.

PubMed Abstract | CrossRef Full Text | Google Scholar

40. Lee SH, Lee JW, Soung YH, Kim HS, Park WS, Kim SY, et al. BRAF and KRAS mutations in stomach cancer. Oncogene (2003) 22:6942–5. doi: 10.1038/sj.onc.1206749

PubMed Abstract | CrossRef Full Text | Google Scholar

41. Takahashi N, Yamada Y, Taniguchi H, Fukahori M, Sasaki Y, Shoji H, et al. Clinicopathological features and prognostic roles of KRAS, BRAF, PIK3CA and NRAS mutations in advanced gastric cancer. BMC Res. Notes (2014) 7:271. doi: 10.1186/1756-0500-7-271

PubMed Abstract | CrossRef Full Text | Google Scholar

42. Lee SH, Ahn BK, Baek SU, Chang HK. BRAF mutation in multiple primary cancer with colorectal cancer and stomach cancer. Gastroenterol Rep. (2013) 1:70–4. doi: 10.1093/gastro/got004

PubMed Abstract | CrossRef Full Text | Google Scholar

43. Lu W, Wei H, Li M, Wang H, Liu L, Zhang Q, et al. Identification of KRAS and PIK3CA but not BRAF mutations in patients with gastric cancer. Mole Med Rep. (2015) 12:1219–24. doi: 10.3892/mmr.2015.3530

PubMed Abstract | CrossRef Full Text | Google Scholar