Association between ABCB1 (3435C>T) polymorphism and susceptibility of colorectal cancer: A meta-analysis.

Studies on the relationship between ABCB1 3435C>T polymorphism (rs1045642) and colorectal cancer (CRC)susceptibility have yielded inconclusive results. To clarify this issue, we undertook a meta-analysis to investigate the relationship between rs1045642 and CRC risk.Three electronic scientific publication databases (Cochrane Library, Pubmed, Embase) were screened using specific search terms. Relevant literature was identified using literature traceability methods. Selected publications were evaluated according to the inclusion and exclusion criteria. Effect size information (odds ratio and the corresponding 95% confidence interval [CI]) was obtained following quality assessment and data extraction from the included publications, and a meta-analysis conducted. Statistical analysis was performed with the Stata sofz (Version 13.0) software.Overall, 17 case-control studies involving 7129 CRC patients and 7710 healthy control subjects satisfied the criteria for inclusion in the meta-analysis. There was no significant association between ABCB1 3435C>T polymorphism and CRC risk in any of the genetic models. In the CC versus CT model (I = 20.9%, Pheterogeneity = .276), CC versus CT + TT model (I = 45.6%, Pheterogeneity = .102) and CT versus CC + TT model (I = 17.8%, Pheterogeneity = .298) analyses, between-study heterogeneities were detected as significant in Asian populations. In the CT versus TT model (I = 24%, Pheterogeneity = .254) and CC + CT versus TT model (I = 0, Pheterogeneity = .55), between-study heterogeneities were found to be significant in groups of different populations.The meta-analysis described here suggests that the ABCB1 3435C>T polymorphism is not related to CRC susceptibility.

Introduction

In a global context, colorectal cancer (CRC) represents a serious threat to human life and health. What is more, it confers an enormous economic burden on society. CRC is the third most commonly diagnosed cancer in the United States. Its estimated morbidity and mortality incidence will be the third among all carcinomas in the United States in 2019.[ The incidence and mortality rates of CRC vary substantially by race/ethnicity.[ Lifestyle difference is also a vital factor leading to the striking variation in CRC morbidity globally.[ Differences in access to prevention, the quality of treatment technology, and the economic level of the patients greatly affect the mortality rate of CRC patients.[ Although the incidence of CRC in the United States has been tapering off in recent years, worldwide, the outlook does not give grounds for optimism.[

The adenosine triphosphate-binding cassette subfamily B member 1 (ABCB1) gene, also known as multidrug resistance gene 1 (MDR1), is located on the chromosomal region 7q21.1 and encodes the P-glycoprotein (P-gp). P-gp is a 1280-amino acid transporter that serves as a genetically polymorphic efflux transporter that removes foreign substances from cells.[ P-gp mediates multiple drug resistance in cancer cells through various signaling pathways, such as the cyclic adenosine monophosphate/protein kinase A pathway,[ the phosphatidylinositol 3-kinase/protein kinase B pathway,[ the Y-box binding protein 1,[ the phosphatase and tensin homolog,[ p53,[ protein kinase C,[ and other protein kinases.[ Previous studies have shown that ABCB1 is overexpressed in a variety of tumors, such as breast cancer, acute myeloid leukemia, hematological malignancies, childhood tumors, and other solid tumors.[

Single nucleotide polymorphisms (SNPs) of ABCB1 affect its expression and function.[ Up to now, numerous SNPs, including some synonymous ones, have been identified in the coding region.[ Three SNPs in the coding sequence (rs1128503, rs1045642, and rs2032582) are the most widely studied in ABCB1; these are relevant to the substrate and inhibitor-dependent functional modifications observed in vitro and reduced expression in tissues.[ The distributions of rs1128503, rs1045642, and rs2032582 differ significantly among races and ethnicities. It is reported that Africans and African-Americans harbor the lowest frequencies of polymorphic alleles and Asians and Caucasians possess the highest.[

The ABCB1 3435C>T polymorphism is a synonymous SNP with no impact on the structural contribution of the amino acid at position 1145 (Ile) in the second ATP binding domain but does affect the expression of P-gp in tissues.[ So far, numerous epidemiological studies have been performed to assess the association between rs1045642 and risk for CRC. However, due to the limitations of individual studies, the results are inconsistent. To help resolve this matter, we performed a meta-analysis based on a total of 17 independent studies, to obtain a more precise estimation of the association between rs1045642 and the risk of CRC. This meta-analysis suggested ABCB1 3435C>T polymorphism is not related to CRC susceptibility.

Materials and methods

Literature search

We queried the Cochrane Library, Pubmed, and Embase databases on August 1, 2019. Keyword combinations for colorectal neoplasms (colorectal, colorectal tumor, colorectal neoplasm, colorectal tumors, tumor colorectal, neoplasms colorectal, neoplasm colorectal, cancers colorectal, cancer colorectal, CRC, carcinomas colorectal, colorectal carcinoma, carcinoma colorectal, colorectal carcinomas) or colonic neoplasms (colonic neoplasm; colon neoplasm; neoplasms colonic, neoplasm colon, colon neoplasms, neoplasms colon, neoplasm colonic, cancer colonic, cancers colon, cancer of the colon, colonic cancer, colon cancer, colon cancers, colonic cancers, cancers colonic, cancer colon, cancer of colon) and gene symbols, and synonyms for the ABCB1 gene (ABCB1, MDR1, CLCS, P-GP, PGY1, ABC20, CD243, and GP170) and polymorphism (polymorphism, SNP, and variant) were used to form a Boolean query formula. Two reviewers (L.H. and Z.Z.), independently and in duplicate, screened titles and abstracts using a standardized data form tested in pilot runs. Inconsistencies regarding inclusion were resolved through consensus. The meta-analysis did not involve data related to patient personal information and therefore does not require ethical approval.

Inclusion and exclusion criteria

The criteria for inclusion in the study were as follows:

manuscripts from peer-reviewed journals;

case-control studies assessing the association between the ABCB1 3435C>T polymorphism (rs1045642) and CRC;

studies focusing on CRC or colonic cancer;

no inconsistencies in genotype data for either cases or controls; and

studies with enough genotype data to estimate the odds ratio (OR) and 95% confidence interval (CI) in at least one genetic comparison model.

The exclusion criteria were:

not case-control studies,

control population including malignant tumor patients, and

duplicate publications.

Two individual authors (L.H. and Z.Z.) performed the literature selection process. Another author (B.Z.) performed an investigation to reach an eventual agreement if the first 2 reviewers came to contradictory conclusions.

Data extraction

Two investigators reviewed and extracted information from all qualified publications based on the inclusion and exclusion criteria listed above. When there was a conflict, the 2 reviewers reached an agreement through discussion. The following information was extracted from each included study: first author's surname, year of publication, ethnicity, total numbers of cases and controls, as well as numbers of cases and controls with CC, CT, and TT genotypes. Individuals of different descent were categorized as Caucasian and Asian. Individuals of different descent were categorized as Caucasian and Asian.

Methodological quality assessment

The quality of the included studies was evaluated by the total score of quality assessment (TSQA).[ Studies were scored according to TSQA standards (Supplementary Table S1). Studies of high quality were given scores of greater than 9.

Statistical analysis

The measure of effect in these studies was the OR with 95% CI. Summary measures were pooled using random-effects models, with the estimate of heterogeneity taken from the Mantel-Haenszel model. All statistical analyses were conducted in the Stata 13 environment. The aggregated estimate of the OR and corresponding 95% CI were calculated for the dominant model (CC + CT vs TT, with C standing for cytosine and T for thymine), the recessive model (CC vs CT + TT), and the overdominant model (CT vs CC + TT). Cochran's Chi-square-based Q test was used to test the heterogeneity assumption. A value of P < .1 in the Q test indicated that the between-study heterogeneity was significant.[ However, when P≥.1, the pooled ORs and 95% CIs should be measured using a fixed-effect model employing the Mantel–Haenszel algorithm.[ To explore the effect of heterogeneity among the studies on the conclusions of this meta-analysis, we performed subgroup analyses by ethnicity. We examined the ABCB1 3435C>T genotypes using dominant (CC + CT vs TT), recessive (CC vs CT + TT), and overdominant (CT vs CC + TT) genetic models, as well as the allelic model (C vs T). The estimated OR and 95% CI were obtained from Forest plots. Publication bias was graphically detected by funnel plots. The symmetry of the funnel plot was further evaluated by Egger linear regression test. The significance of the intercept was determined by the t test suggested by Egger, where P < .05 was considered representative of statistically significant publication bias.

Results

Retrieval of studies and their characteristics

From the searches for studies on colorectal neoplasms or colonic neoplasms and C3435T genotypes, 153 potentially eligible records were identified. Titles and abstracts of these records were screened for inclusion. Seventeen independent studies met the inclusion criteria, consisting of six Asian and eleven Caucasian populations (Fig. 1).[ In total, 7,179 CRC cases and 7,710 controls were included in the meta-analysis. The characteristics of the selected studies are summarized in Table 1.

Figure 1

Flow diagram of literature search and screen.

Table 1

Main characteristics of studies included in this meta-analysis.

Association of rs1045642 C>T and CRC

A total of 17 independent studies consisting of 7129 CRC patients and 7710 healthy controls were included in the analysis of the association of ABCB1 3435C>T polymorphism with susceptibility to CRC. All possible genetic models were analyzed to seek potential differences in genotypic and allelic frequencies regarding ABCB1 3435 C>T polymorphism amongst CRC cases and controls. We used the random effect model for the analysis of all genetic models: we did not find a significant association between the different genotypes of SNP rs1045642 C>T and susceptibility to CRC in any of the models (Table 2; Fig. S1). Moreover, we found no association between ABCB1 3435 C>T polymorphism and CRC when comparing the C and T alleles (Table 2; Fig. S2).

Table 2

Main results of rs1045642 polymorphism and colorectal cancer risk in this meta-analysis.

Subgroup analysis based on ethnicity

In general, between-study heterogeneities were not significant in any of the genetic models for the association between ABCB1 3435C>T polymorphism and CRC, but in the stratified analysis by ethnicity, they were present in some genetic models (Table 2). Eleven studies consisting of 5369 CRC cases and 5653 controls were included in the Caucasian group, while six studies comprising 1760 CRC cases and 2057 controls were enrolled in the Asian group. We used the random-effect model for examining heterogeneity using the genetic models described above. For the Caucasian group, no significant between-study heterogeneity was detected (Figs. 2–4; Fig. S3). For the Asian group, the between-study heterogeneities increased strikingly. In the CC versus CT model (I = 20.9%, Pheterogeneity = 0.276), CC versus CT + TT model (I = 45.6%, Pheterogeneity = 0.102), and CT versus CC + TT model (I = 17.8%, Pheterogeneity = 0.298), between-study heterogeneities were determined to be significant (Table 2; Fig. 2). Moreover, there was no significant association detected for the C allele versus T allele comparison (Table 2; Fig. S4).

Figure 2

Forest plot on association between rs1045642 polymorphism and colorectal cancer risk, stratified by ethnicity. (A) CC versus CT model, random-effect pooled OR = 1.074, 95% CI: 0.884–1.304, P = .473, I = 76.8%, Pheterogeneity = .000. (B) CC versus CT + TT model, random-effect pooled OR = 1.013, 95% CI: 0.866–1.185, P = .872, I = 67.4%, Pheterogeneity = .000. (C) CT versus CC + TT random-effect pooled OR = 0.899, 95% CI: 0.759–1.065, P = .218, I = 80.3%, Pheterogeneity = .000. 95% CI = 95% confidence interval, OR = odds ratio.

Figure 4

Funnel plot analysis to detect publication bias for rs1045642 polymorphism. (A) Recessive model. (B) Dominant model. (C) Overdominant model. Each point represents a separate study for the indicated association.

Subgroup analysis based on the characteristics of patients used as study controls

To further explore other potential sources of heterogeneity, we stratified all the studies according to attributes of study controls. In eleven studies, the participants were 2588 patients with CRC and 2912 people who were in the hospital for unrelated problems (hospital-based (HB) group), while in six studies the participants were 4541 patients with CRC and 4798 healthy individuals who were selected from the general population (population-based (PB) group). All genetic models were evaluated using the random-effect statistical model. For the PB group, the striking between-study heterogeneities were minimized. For all genetic models except for the CT versus TT model (I = 24%, Pheterogeneity = .254) and CC + CT versus TT model (I = 0, Pheterogeneity = .55), between-study heterogeneities were found not to be significant (Table 2; Fig. 3). For the HB group, there was no significant between-study heterogeneity detected (Table 2; Fig. 3; Fig. S5). Moreover, there was no significant association detected for the comparison of the C and T alleles (Table 2; Fig. S6).

Figure 3

Forest plot on association between rs1045642 polymorphism and colorectal cancer risk, stratified by source. (A) CT versus TT model, random-effect pooled OR = 0.877, 95% CI: 0.712–1.08, P = .215, I = 79.7%, Pheterogeneity = .000. (B) CC + CT versus TT model, random-effect pooled OR = 0.914, 95% CI: 0.767–1.09, P = .318, I = 74.1%, Pheterogeneity = .000. 95% CI = 95% confidence interval, OR = odds ratio.

Sensitivity analyses

The leave-one-out sensitivity analysis confirmed the robustness and reliability of our drawn conclusion. The association between SNP rs1045642 C>T and CRC remained insignificant after the removal of any included study (detailed data not shown).

Publication bias

Begg funnel plot and Egger test were performed to assess the publication bias. No apparent asymmetry of funnel plots was detected on visual inspection. Egger test was used to provide statistical evidence for the funnel plot. In the recessive, dominant, and overdominant models, the P values for Begg funnel plot were .201, .57, and .104, respectively, indicating that there was no significant publication bias (Fig. 4). There was also no publication bias in comparisons of any of the other genetic models (Fig. S7).

Discussion

Although the pathogenesis of CRC is multifactorial and mostly unclear, both internal and external factors are considered to contribute to its etiology. The multifactorial nature of the pathology of CRC calls for the quantitation of the independent risk factors. Due to the leading role of genetic factors in the pathogenesis of CRC, an in-depth understanding of the correlation of gene polymorphisms to CRC will help predict the course of the disease and take preventive measures, as well as identify potential targets for specific drug therapy. According to the findings of our present study, there was no statistically significant correlation between ABCB1 3435C>T polymorphism and the risk of CRC, regardless of the ethnicity of the study subjects and the environment of the healthy controls.

In recent years, researchers have carried out a number of studies on the association between ABCB1 gene polymorphisms and the susceptibility to various tumors. As a systematic approach that uses statistical analysis, meta-analysis is an effective way to come to conclusions about different studies that are inconsistent due to the limitations of the individual studies. A meta-analysis performed by Razi et al suggested that ABCB1 3435C>T polymorphism was not associated with the risk of multiple myeloma.[ The meta-analysis carried out by Sharif et al suggested that ABCB1 3435C>T polymorphism might be a genetic risk factor and a potential biomarker for breast cancer,[ but Tazzite et al reported that it was not associated with breast cancer risk in Morocco.[ Research by Wu et al showed that the ABCB1 C3435T polymorphism was not associated with susceptibility to gastric cancer.[ The study executed by He et al indicated that ABCB1 3435C>T polymorphism was associated with CRC risk in Asians, but Zhang et al found there was no significant association between them.[

In recent years, we have been concerned with the study of ABCB1 C3435T polymorphism and the risk of CRC. In the present work, 17 case-control studies were selected from Cochrane Library, PubMed, Embase databases.[ The results of these studies are inconsistent. Patients’ gender, age, race/ethnicity, pathological type, pathological grade, attributes of study controls, and other factors may account for inconsistencies in the results of the studies. Unfortunately, we did not get full information about patients and controls. Despite that, a comprehensive analysis was performed to identify the possible association between ABCB1 3435C>T polymorphism and CRC susceptibility. Several potential limitations can be noted in the present analysis that are inherent to any meta-analysis. First, the total sample size was not large enough, and subgroup analysis could only be performed on the basis of ethnicity and the selection of study controls. Second, only studies written in English were included in this meta-analysis. Third, some studies were of poor quality. Thus, a larger meta-analysis containing more studies, other ethnicities, and gender should be conducted in the future to improve the reliability of the conclusions.

Acknowledgments

This work was supported by the National Natural Science Foundation of China (81903031), Zhoukou Normal University Guiding Project (ZKNUC2018011). We thank LetPub (www.letpub.com) for its linguistic assistance during the preparation of this manuscript.

Author contributions

Guo-yin Li designed the study and drafted the initial manuscript. Li-li Han, Bai-le Zuo and Wei-liang Cai contributed to initial data analysis and interpretation. Zhen-ni Guo, Bing-hua Tong and contributed to the production of the pictures. Guo-yin Li and Zheng zhu supervised all aspects of the study, critically reviewed and revised the manuscript, and approved the final manuscript as submitted.