Abstract

Background. Colorectal cancer (CRC) is a major health concern worldwide. A series of sequential accumulation of genetic and epigenetic changes are responsible for the initiation and progression of diseases via the normal > adenoma > carcinoma sequence. Genetic variants in crucial cancer-causing genes are known to mediate the risk of cancer. Objective. In this case-control study, we examined single nucleotide polymorphism (SNP) in HER1 (rs763317 and rs3752651) and HER2 (rs1136201 and rs1058808) genes to assess their role in the susceptibility of CRC in a Saudi population. Methods. TaqMan allelic discrimination assay was utilized to identify the genotypes in 163 normal and 143 CRC patients. Results. In the overall analysis, the rs3752651 and rs1136201 were significantly associated with the risk of CRC. Although none of the examined SNPs had any impact on the age at which CRC was diagnosed, interestingly, three SNPs showed a significant association based on gender. The rs3752651 conferred significant protection only in men, whereas rs1136201 diminished the risk and rs1058808 considerably increased the susceptibility of CRC only in women. Conclusions. Our result suggests that these SNPs in HER1 and HER2 after validation in larger cohorts of different ethnicities may be utilized as genetic screening markers for predicting colorectal cancer predisposition.

1. Introduction

Colorectal cancer (CRC) is one of the leading causes of cancer-related mortality and morbidity around the world [1]. As per gender segregation, it is the third commonly diagnosed cancer in men and ranks second in women [2]. In 2018, the estimated total number of new cases was 1,849,518 while the number of deaths was 880,792 [3]. The incidence and mortality of CRC in men are approximately 20% higher than in women. In Saudi Arabia, both the genders combined, CRC is the second most common cancer with an estimated incidence of 3564 (14.6%) cases and mortality of 1603 (15.2%) in the year 2018 [3].

Colorectal tumorigenesis is a multistep process involving accumulation of mutations in oncogenes and tumor suppressor genes leading to the sequential transformation of the colorectal normal epithelium to adenoma followed by invasive and metastatic tumor [4]. Various risk factors including ethnicity, environment, and genetics have been attributed towards the predisposition of CRC [5]. SNPs are considered to be the main source of genetic variation which accounts for about 90 percent of all human genetic variations occurring roughly every 100–300 bases [6]. Most SNPs are usually located in noncoding regions, and their distribution is not under selective pressure [7]. However, SNPs occurring in coding sequences of genes may lead to protein expression changes that directly contribute to a disease phenotype [8]. The location of these SNPs near genes that are involved in a crucial role such as gene expression and regulation, signal transduction, and the control of cell cycle may affect the function of such genes. An important molecular mechanism involved in the development and the progression of CRC is the mitogen-activated protein kinase (MAPK) cascade downstream of the epidermal growth factor receptor (EGFR) [9].

The EGFR also known as human epidermal growth factor receptor 1 (HER1) belongs to the ErbB family of transmembrane protein receptors that function as a receptor tyrosine kinase (RTK). The ErbB family consists of four members including ErbB-1 (HER1/EGFR), ErbB-2 (HER2), ErbB-3 (HER3), and ErbB-4 (HER4) [9]. Previous studies have reported that HER1 and HER2 gene polymorphisms are associated with several malignancies such as colon [10], breast [11, 12], and lung cancers [13]. Additionally, HER2 polymorphism has also been linked with higher expression of the protein in breast and ovarian cancers [14, 15].

Large efforts are being made in identifying trustworthy biomarkers for early detection of human malignant cancers. Genetic variation analysis such as SNPs can be used for predicting novel diagnostic and prognostic genetic markers. In this study, we examine the effect of four SNPs, rs763317 and rs3752651 located in HER1 and rs1058808 and rs1136201 located in HER2, on CRC susceptibility in Saudi Arabian patients by comparing the genotypic and allelic distribution of these SNPs in cancer cases to that of healthy subjects.

2. Materials and Methods

2.1. Sampling and DNA Isolation

Colorectal cancer patients (n = 143) included all ages (median age = 58 years) and histopathologically confirmed cases of the disease. Approximately 3 ml of blood samples from colorectal cancer patients attending the endoscopy department of King Khalid University Hospital, King Saud University, was collected in ethylenediaminetetraacetic acid- (EDTA-) containing vacutainers prior to therapy. Additionally, age- and gender-matched normal subjects with no history of cancer were recruited as controls. The study was reviewed and approved by the institutional review board of College of Medicine, King Saud University (Project approval No. E-12-596). Written informed consent was obtained from all the patients. Genomic DNA was isolated from peripheral blood cells using the QIAmp DNA blood mini kit (Qiagen, USA) according to manufacturer’s instructions. Quantitation and purity of DNA were examined on a NanoDrop 8000 spectrophotometer (Thermo Scientific, USA).

2.2. Genotyping

Based on previous studies, four SNPs in HER1, rs763317 and rs3752651, and HER2, rs1058808 and rs1136201, were selected for this study [10, 1619]. TaqMan allelic discrimination assays (Thermo Fisher Scientific, USA) were used for genotyping colorectal cancer and control peripheral blood DNA samples for these SNPs. The TaqMan assay ID and number of samples for which genotype calls were obtained for each SNP are as follows: rs763317 (assay ID = C___2310200_10, cancer = 140, control = 149), rs3752651 (assay ID = C____326500_10, cancer = 143, control = 163), rs1058808 (assay ID = C___1551672_20, cancer = 139, control = 160), and rs1136201 (assay ID = C___7452451_1_, cancer = 141, control = 162). Each reaction contained 20 ng of DNA, 5.0 μL of 2X TaqMan genotyping Master Mix (Applied Biosystems, USA), and 0.25 μL of 40X TaqMan assay in a total volume of 10 μL. Genotyping experiments were conducted by endpoint reading on the QuantStudio 7 Flex Real-Time PCR machine (Applied Biosystems, USA), and genotype calls were made by TaqMan Genotyper Software version 1.6.

2.3. Statistical Analyses

Genotype and allele frequencies of polymorphic variants examined in HER1 and HER2 genes were determined, and tests for association and Hardy–Weinberg equilibrium were performed on a web-based tool at https://ihg.helmholtz-muenchen.de/cgi-bin/hw/hwa1.pl. For each SNP, odds ratios, confidence intervals, chi-square, and significance were calculated to assess the existence of an association between the SNP and colorectal cancer. Statistical significance was kept at a value of ≤0.05. Bonferroni correction was used for multiple comparison with an α = 0.0125 considered as significant.

3. Results

In the present study, the genotypes of four germline SNPs residing in HER1, rs763317 and rs3752651, and HER2, rs1058808 and rs1136201, were evaluated in order to determine their association with colorectal cancer risk in Saudi Arabian population. The rs763317 and rs3752651 are located in the intron of HER1, whereas the rs1058808 (P1170A) and rs1136201 (Ile655Val) are nonsynonymous coding SNPs in HER2.

The distribution of genotypes of these SNPs in both the control and cancer groups was examined for Hardy–Weinberg equilibrium (HWE). All the SNPs followed Hardy–Weinberg equilibrium except rs3752651 () in cancers (Table 1). A high rate of consanguinity in the Saudi Arabian population could be a probable cause for this lack of Hardy–Weinberg equilibrium [20].

3.1. Association of HER1 and HER2 Variants with Colorectal Cancer Risk

The distributions of genotype and allele frequencies of the chosen SNPs within the HER1 and HER2 genes and their association with colorectal cancer were evaluated between both colorectal cancer patients and healthy individuals (Table 2). The rs3752651 within HER1 showed a statistically significant association with colorectal cancer susceptibility. In the control group, the frequency of the CC genotype of rs3752651 was detected at 4%, whereas none of the colorectal cancer patients harbor this genotype. It was estimated that individuals carrying the CC genotype of rs3752651 had about 11-fold lower risk of developing colorectal cancer compared to those with the TT genotype (OR = 0.090, χ2 = 5.03, ). Furthermore, the rs1136201 variant in the HER2 gene revealed a significant protective association with CRC. Individuals with the minor allele G of rs1136201 were 80% less likely to develop CRC compared to those having A allele (OR = 0.555, χ2 = 4.45, ). The other two SNPs, rs763317 (HER1) and rs1058808 (HER2), did not show a significant association with the predisposition of colorectal cancer in the overall analysis (Table 2).

3.2. Association of HER1 and HER2 SNPs with Colorectal Cancer Based on the Age of Onset of the Disease

To determine whether HER1 SNPs s763317 and rs3752651 and HER2 SNPs rs1058808 and rs1136201 have any association with the age at colorectal cancer diagnosis, colorectal cancer patients were stratified based on the median age at the time of disease diagnosis as ≤58 and > 58 years, and the genotype and allele frequencies were compared to the age-matched controls. Our data demonstrated that none of the analyzed SNPs had a noticeable association with colorectal cancer suggesting that these variants do not influence early or late onset of the disease (Table 3).

3.3. HER1 and HER2 SNPs’ Association with Colorectal Cancer Risk Based on Gender

We investigated whether gender plays any role in the association of HER1 SNPs, rs3752651 and rs763317, and HER2 SNPs, rs1058808 and rs1136201, with CRC. Our data indicated that the CC genotype of rs3752651 imparts significant protection against CRC only in men (OR = 0.077, χ2 = 5.60, ). Additionally, both SNPs examined in HER2 showed a significant association with CRC but only in women. The rs1058808 increased the risk while rs1136201 was protectively associated with CRC. For the rs1058808 loci, CC homozygotic women showed approximately 3.8-fold higher risk of developing colorectal cancer compared to those having the GG genotype (OR = 3.875, χ2 = 5.65, ). According to the allelic model as well, the minor allele C of rs1058808 conferred a significant association with CRC (OR = 1.665, χ2 = 4.04, ). The minor allele G of rs1136201 SNP in HER2 lowers the susceptibility of developing colorectal cancer in women about 4-fold compared with those having A allele (OR = 0.244, χ2 = 7.26, ). This association persisted even after Bonferroni correction for multiple comparisons. Moreover, the AG heterozygous women for this SNP also had significant protection against this cancer compared to those with the AA genotype (OR = 0.275, χ2 = 5.33, ) (Table 4).

4. Discussion

The HER family of receptors are transmembrane tyrosine kinases; they play a key role in the regulation of various cellular functions in response to extracellular signaling. The basic aspect of signaling in the HER family is the interaction of two receptors [21]. When HER family members dimerize, it leads to transphosphorylation of their intracellular domain generating the initial signal that activates several downstream signaling pathways [21]. HER family signaling is deregulated in several cancer types, mainly in breast and lung cancers [1113].

In the present study, we evaluated the role of SNPs rs763317 and rs3752651 residing in HER1 and rs1058808 and rs1136201 in HER2 genes in colorectal cancer in a Saudi population. Very few studies have examined the association of these SNPs with colorectal cancers. To the best of our knowledge, this is the first study to evaluate the association of HER1 and HER2 gene polymorphism with colorectal cancer risk in a Saudi Arabian population. We have demonstrated that the HER1 SNP rs3752651 and HER2 exonic SNP rs1136201 are protectively associated with colorectal cancer in the overall analysis. Additionally, gender segregation of the population revealed that the association with CRC was statistically significant in men for rs3752651 and women for rs1136201. While the CC genotype of rs3752651 provided significantly increased protection (OR = 0.077, ) relative to the TT genotype in men, such association was not observed in women. Contrary to our finding of protective association in the overall analysis (OR = 0.090, ), Poole et al. demonstrated that the CC genotype of rs3752651 conferred 40% increased risk of colon cancer relative to those having the TT genotype. However, this association was not statistically significant for rectal cancer as well as when the cases were pooled together as colorectal cancers [10]. In glioma patients, the CC genotype of rs3752651 was found to be significantly associated with decreased overall and progression-free survival compared to the TT genotype [22].

Contradictory results pertaining to the association of rs1136201 in HER2 with various malignancies exist. While in the overall analysis of our population, the genotypes of rs1136201 were not associated with colorectal cancers, we observed a protective association with the G allele conferring 1.8-fold lower risk of developing the disease relative to the A allele. Similar to our findings, the genotypic analysis of rs1136201 in the Indian as well as Scottish Caucasian population did not show a significant correlation with colorectal carcinoma [17, 18]. Interestingly, when our population were segregated based on gender, the rs1136201 demonstrated that women carriers of the AG genotype and the G allele were at 3.6-fold (OR = 0.275, ) and 4.1-fold (OR = 0.244, ) lower risk of colorectal malignancy relative to the AA genotype and A allele, respectively. Similar to our findings, de Almeida et al. reported a negative association in a dominant model of analysis comparing AG + GG versus AA as well as G allele with breast cancer in a Brazilian population [23]. Additionally, Nelson et al. in a United States population from Wisconsin suggested an inverse relationship of the GG genotype of rs1136201 with invasive breast cancer. This inverse association was stronger in women aged more than 55 years, those without a family history of breast cancer, postmenopausal women with higher body mass index, and cases with regional or distant metastasis [24]. Few investigations tried to address the functional consequence of A or G allele of rs1136201 with HER2 protein overexpression. Millikan et al. found that breast tumors from women that possess the GG genotype displayed higher expression of HER2, but the results did not reach statistical significance [14]. In ovarian cancer patients from a South Indian state of Tamil Nadu, Shanmughapriya et al. demonstrated that the GG homozygotes of rs1136201 strongly expressed HER2 protein [15]. A similar association of the GG genotype and HER2 protein overexpression as well as patient survival was not observed in colorectal cancer patients [18].

Besides, HER2 gene SNP rs1058808 which did not show any association in the overall comparison demonstrated that women harboring the CC genotype were at 3.8-fold increased risk of developing colorectal cancer compared to those with the GG genotype. The rs1058808 (c.3418 C > G), also known as P1170A, is a missense variant located in exon 27 of the HER2 gene [12, 19]. Since rs1058808 and rs1136201 variants are in the coding region, they might affect the tyrosine kinase activity or structure of the protein making them unstable. HER2 gene nsSNP rs1058808 has been reported to show a negative change in Gibb’s free energy; ΔΔG value more than 1 which is −1.07 is considered to be less stable by the I-Mutant 2.0 server. The Pro1170Ala change of rs1058808 is a change from polar to nonpolar amino acid that has been reported to be deleterious according to the SIFT program [25].

Similar to our overall analysis, Han et al. did not find any correlation with genotypes of rs1058808 and rs1136201 and the occurrence of breast cancer in Korean women. However, a haplotype comprising these two nonsynonymous tagging SNPs was linked with higher HER2 expression and significantly poorer prognosis [19]. Su et al. also reported a significant association between rs1058808 and HER2 protein expression. They found that patients with CG and GG genotypes of rs1058808 demonstrated higher HER2 protein expression relative to those having the CC genotype () [26]. The mechanism related to the higher HER2 protein expression with genetic variants needs further investigations. Conflicting reports pertaining to the association of rs1058808 and rs1136201 with different cancer types have been reported [11, 13, 2729]. In view of the involvement of HER1 and HER2 in several malignancies such as breast cancer, ovarian cancer, non-small-cell lung cancer, head and neck cancers, gastric cancer, and hepatocellular carcinoma, it would be interesting to examine the association of SNPs tested in this study with these cancers [30, 31].

Despite significant findings, the drawback of our study is that the sample size is relatively small and from a single hospital in the capital city Riyadh which receives patients from all over the country but may not represent the entire Saudi population. Additionally, since consanguinity is common in Saudi Arabia, it may have introduced some bias in our results. Another limitation of our study is that we took a simplistic view and performed an association of individual SNP with CRC; however, haplotype linkage disequilibrium analysis together with additional SNPs on same chromosome might provide broader insights into the role of HER1 and HER2 variants in the causation of the disease.

In conclusion, our preliminary analysis on selected SNPs in HER1 and HER2 genes suggests a statistically significant association of rs3752651, rs1058808, and rs1136201 with colorectal cancer in a Saudi population. These associations were gender specific and not related to the age of onset of the disease. The rs3752651 in men and rs1136201 in women showed protective association, while rs1058808 increases the risk of colorectal cancers in women. Although our results are significant, due to the small sample size, it should be validated in a larger cohort of diverse ethnicities along with haplotype linkage analysis before utilizing these SNPs as screening markers for colorectal cancers.

Data Availability

All data used to support the findings of this study are included within the article.

Conflicts of Interest

The authors declare no conflicts of interest.

Acknowledgments

The authors are grateful to the Deanship of Scientific Research, King Saud University, for funding through the Vice Deanship of Scientific Research Chairs.