Variation rs9929218 and risk of the colorectal Cancer and adenomas: A meta-analysis

Backgrounds Genome-wide association studies (GWAS) have identified multiple common CRC-related (colorectal cancer) SNPs (single nucleotide polymorphisms) including the Cadherin 1(CDH1) rs9929218 may act by increasing the risk of colorectal cancer, colorectal adenoma, or both. These studies, however, reported inconsistent associations. Methods To derive a more accurate approximation of the connection, we carried out a meta-analysis of 12 published pieces of research including 11,590 controls and 8192 cases. We used odds ratios (ORs) and 95% confidence intervals (CIs) to evaluate the associations’ strength. Results Meta-analysis implied considerable association between CRC and rs9929218 (OR = 1.21, 95%CI 1.04–1.42 for GG versus AA; OR = 1.22, 95%CI 1.05–1.42 for GG/AG versus AA). In the subgroup analyses, significantly increased risks were found among Europeans. Conclusions In summary, our meta-analysis studies in different populations confirmed that SNP rs9929218 is significantly associated with CRC risk and that this variant may have a greater impact on Europeans.


Background
More than one million people worldwide are affected by colorectal cancer (CRC) every year [1]. It is currently the third most frequent malignancy and the fourth commonest cause of cancer-related mortality in the world [2], and accounts for approximately 630,000 death from CRC annually. Previous genetic epidemiological studies provided evidence that CRC is a complex disease influenced by environmental and genetic factors and their interactions [3]. The majority of CRC are developed from colorectal adenomas (CRA) [4]. It is well-recognized that the high genetic risk of CRC is due in large part to the susceptibility to adenomas [5]. Single nucleotide polymorphism associated with CRC may increase the risk of colorectal cancer, colorectal adenomas, or both.
Among these SNPs, rs9929218 (16p22.1), located in the intron region of the gene cadherin 1 (CDH1), was identified to be associated with CRC risk [13]. In 2016, Han et al. conducted a meta-analysis of rs9929218 including 16 studies (n = 131,768) which emphasizes a significant association between rs9929218 polymorphism and CRC susceptibility [13]. Nonetheless, its limitation is the absence of raw genotype data that reference dominant and recessive models.
The rs9929218 has been identified as an aroused general interest for CRC susceptibility by recent genomewide association studies and this polymorphism has shown that the G allele is associated with an increased risk of colorectal cancer. However, some of the literature has produced contrary results [14,15].
We, therefore, performed a meta-analysis of the published studies to clarify this antilogy and constitute a comprehensive map of the relationship between 16q22.1 (rs9929218) polymorphism and the CRC susceptibility.

Methods
Our meta-analysis is reported followed the guideline of the PRISMA (Preferred Reporting Items for Systematic Review and Meta-Analysis) statement [16]. There are no ethical issues involved in our study because our data were based on published studies, and no ethical issues were involved in the selection, extraction and analysis of the data.

Identification and eligibility of relevant studies
Embase, PubMed, and ScienceDirect database were hunted invoking the semesters: 'rs9929218', 'Single Nucleotide Polymorphism ','colorectal cancer 'and 'colorectal adenoma 'to gather eligible articles. We searched literature published before February 2020 and the language was restricted to English. This search was implemented by two students independently to determine potential publications meeting our inclusion criteria. If there is a clash, two students consult with the third reviewer and achieve a consensus [3,17]. Studies in our meta-analysis ought to meet the following inclusion standards: (1) assess the association between colorectal cancer risk and the rs9929218 polymorphism; (2) the studies were designed as case-control studies; (3) include the detailed frequency of each genotype.

Data extraction
Two researchers appraised and extracted all data from all qualified publications independently, as reported by the inclusion standards that were listed above. The following data was gathered from each study: the name of the first author, the year of publication, the country of study, the specific number of genotypes in each control and case group; and the source of the control group (population-or hospital-based controls). Varying ethnic ancestries were classified as European, Asians, or Other including more than one ethnicity's subjects [13,17]. There is no duplicate sample was included in these 12 studies.

Statistical analysis
The strength of the association between the CRC risk and rs9929218 polymorphism was measured through odds ratios (ORs) with 95% confidence intervals (CIs). The statistical significance of the pooled OR was determined using the Z-test [17]. First of all, we measured four genetic models: the homozygous model (GG vs. AA), the heterozygous model (GA vs. AA), the dominant model (GG + GA vs. AA), and the recessive model (GG vs. GA + AA) respectively (A: wild allele; G: mutated allele). Stratified analysis based on race, disease type, and source of control.
In return for the probability of heterogeneity throughout the researches, a statistical test for heterogeneity was performed based on the chi-square-based Q test [17,18]. P < 0.05 was considered significant for the heterogeneity. We used two kinds of metaanalysis models including the Mantel-Haenszel (M-H) fixed-effect model and the DerSimonian-Laird (D-L) random-effect model to compute the pooled OR. With the proviso that no significant heterogeneity (I 2 is 50% or less) learnings involved in, the pooled OR is calculated by the fixed-effect model; otherwise (I 2 more than 75%) the OR is calculated by the randomeffect model [17]. We did a sensitivity analysis to evaluate the stability of the results, one single study in the meta-analysis was excluded each time to report the influence of the respective data set to the pooled OR. We offer a potential publication bias with a funnel plot of diagnosis. All analyses were performed using STATA software (Stata Corporation, College Station, TX; version 12.0).

Characteristics of studies
Based on our inclusion criteria, we selected eligible articles from the PubMed, Embase, and Science Direct databases. Twelve eligible tests were included in the metaanalysis, involving data from 13 groups, including 6191 patients and 9314 controls associated with CRC and 2001 patients and 2276 controls associated with CRA [14,15,[19][20][21][22][23][24][25]. Genotype distribution was consistent with Hardy-Weinberg equilibrium in all controls studied (P > 0.05). The main characteristics of the study are shown in Fig. 1. And the raw data is in Table 1.

Quantitative synthesis
We computed the overall OR by the fixed effect model given no significant heterogeneity in all the selected   studies. The main results of the meta-analysis for rs9929218 were listed in Table 2. For the increased risk of CRC and rs9929218 polymorphism, there was no significant evidence shows that the correlations were found between them when all eligible studies were integrated into the meta-analysis. As shown in the Fig. 2. For the rs9929218 polymorphism and CRC risk, our metaanalysis showed the general allele to that OR was 1.14 (95%CI 1.00-1.28, P = 0.042), and the relevant outcomes underneath the recessive and dominant genetic models were 1.13 (95%CI 1.00-1.27, P = 0.044) and 1.03 (95%CI 0.97-1.9, P = 0.334). Significant associations were consequently discovered for the dominant model and overall model. A subgroup meta-analysis was further conducted by us in the type of disease, ethnicity, and the source of controls. As shown in the Fig. 3. In stratified by ethnic analysis, a significant association between rs9929218 is supported by the results of the European population and the CRC. In the subgroup analysis, GG genotypes were found to have a significantly higher risk (OR = 1) than AA genotypes in Europeans (OR = 1.21, 95%CI 1.04-1.42, P = 0.016), AG versus AA genotype (OR = 1.24, 95%CI 1.06-1.45, P = 0.008) and GG/AG versus AA genotype (OR = 1.22, 95%CI 1.05-1.42, P = 0.010). Such an association, however, was not obtained in Asian and the Others group. Considering the type of disease and the source of control, the significant association between risk and mutation in patients with CRC / CRA was not found, and there was no such correlation in each group of control groups.

Sensitivity analyses and publication bias
We then evaluated the influence of each study on the combined OR by sequentially omitting each study and the outcomes suggested that none of the studies Fig. 2 Forest plots for the meta-analysis of rs9929218 polymorphism for overall using a random model. a the GG versus AA. b the GA versus AA. c the dominant model. d the recessive model. OR, odds ratios. CI, confidence interval. CRC, colorectal cancer. CRA, colorectal adenomas meaningfully changed the pooled OR, indicating that the pooled OR of this polymorphism was robust. Perform funnel plots and Egger tests to estimate publication bias. As shown in Fig. 4, no significant asymmetry in the shape of the funnel plots was observed, indicating that no significant publication bias was presented in the study.

Discussion
CRC is one of the gastrointestinal cancers with a high degree of malignancy, aggressiveness, and rapid metastasis. The pathogenic factors of CRC are complicated, as other complex diseases, including both genes and the environment cause [3]. GWAS has different common single nucleotide polymorphism (SNP) associated with colorectal cancer risk. One of SNP with the strong association signal was rs9929218 at 16q22.1 [26,27]. GWAS and replication studies conducted mostly in European populations but less in other ethnic populations.
Colorectal adenomas (CRA) is well-recognized premalignant lesions of CRC and the majority of CRC are formed from adenomas [3]. It has been reported that many genetic risks of CRC are partially mediated by susceptibility to adenomas. CRC-related SNPs may act by increasing the risk of CRC, CRA, or both [28]. CRA is a recognized precursor of CRC based on epidemiologic, histological, and genetic studies demonstrating shared genetic alterations. Several risk factors are related to the risk of developing CRA that is based on epidemiologic researches including cigarette smoking [29], alcohol intake [30], and obesity [31][32][33][34]. Siegert et al. observed significant interactions between alcohol consumption and rs9929218 [35]. In summary, we did a meta-analysis of the relationship between rs9929218 and increased susceptibility to CRC / CRA.
Rs9929218 is situated in the intron zone of CDH1, which encodes a calcium-dependent glycoprotein (Ecadherin). E-cadherin is encoded, a classical calcium adhesion protein, by The CDH1 gene essential for establishing and sustaining cell polarity, tissue constitution, and cell-to-cell adhesion. E-cadherin's low expression in CRC patients is connected with miserable prognoses, like invasive neoplasm development and metastasis in colorectal cancer [26,[36][37][38]. It is considered Ecadherin to be an important factor in epithelialmesenchymal transition (EMT), an important cellular program during tumor cell adhesion, migration, invasion, and metastasis [39]. The EMT is Every point represents a separate study for the indicated association. Funnel plot (a) for publication bias analysis of the selected studies investigating the association between rs9929218 polymorphism and CRC/CRA. The X-axis stands for the ORs and the Y-axis is the standard error for each of the 12 studies. A linear regression based approach proposed by Egger test (b) and Begg plot (c) is used to evaluate the asymmetry of the funnel plot. CI, confidence interval. SND, standard normal deviate suggested to be the first step in the metastatic of cancer cells by emerging evidence. Loss or reduction of E-cadherin expression is considered to be the primary and most important step in the EMT process, and EMT is a critical step in cancer metastasis [39,40].
Our meta-analysis involved 12 studies, 13 groups of data, including 8192 cases and 11,590 controls. The combined results demonstrated that the rs9929218 polymorphism was an increased risk factor for CRC in European populations and there was no significant heterogeneity throughout the study. Moreover, sensitivity analysis and publication bias suggest that our results are robust.
When stratified by ethnicity, significant associations between rs9929218 polymorphism and increased incidence of CRC were observed in Europeans. No clear association, however, was observed between Asia and other races in all genetic models, indicating that the conceivable causes are distinctions in habitat and genetic backgrounds. Besides, apart from genetic factors, multiple lifestyles for varying ethnic groups ought to also be reckoned for CRC risk.
Our meta-analysis's limitations should be regarded in the interpretation of the consequences. Firstly, the number of our case-control study enrolled in our analysis was insufficient. Secondly, we have collected only English literature in this meta-analysis. Thirdly, the involved ethnicities were limited (only Asian, European, and Other population) and further studies of larger sample sizes were needed to explore the influences of different ethnicities. Eventually, this study did not write about some of the potential factors such as age, sex, smoking, alcohol intake, lifestyles, and environmental factors.

Conclusion
In summary, our meta-analysis studies in different populations confirmed that SNP rs9929218 is significantly associated with CRC risk and that this variant may have a greater impact on Europeans than Asians and Others. Due to the limited research available on the current non-European population, further research is needed, including a broader range of population areas and themes. Furthermore, gene-environment interactions and genegene must also be reckoned in succeeding research.