Abstract

Colorectal cancer (CRC) is the 5th leading cancer in China. Alcohol consumption has been reported to be one of the risk factors of CRC. However, it remains unclear whether genetic variants of alcohol metabolic genes are associated with CRC risk. In this study, we tested the coding variants in the alcohol metabolic genes and the risk of CRC, by using 485 cases and 516 controls. A total of 16 germline coding variants in 10 alcohol metabolic genes were genotyped. We identified that rs3741178 in ALDH3B2 was significantly associated with CRC risk with odds ratio being 2.13 (95% CI: 1.24–3.68, ). Further functional annotation suggested that this variant may damage the protein function of ALDH3B2. Our results suggested that ALDH3B2 in the alcohol metabolism pathway contributed to the development of CRC, which may contribute to the prevention of this disease in the future.

1. Introduction

Colorectal cancer (CRC) ranks the fifth leading cause of cancer death among both men and women in China [1]. Epidemiology studies suggest that overnutrition, obesity, tobacco smoking, and heavy alcohol drinking are risk factors of this cancer [2]. Previous genome-wide association study (GWAS) has identified nearly 100 genetic loci that were associated with CRC susceptibility, mostly located in the noncoding region [3]. Recently, several exome-wide association studies identified that functional coding variants also played an important role in the susceptibility of multiple cancers [46]. Therefore, more variants, especially functional coding ones, still need to be further explored.

Alcohol drinking was an important risk factor for colorectal cancer [7]. Gene-drinking interaction was proved to be associated with risk of upper aerodigestive tract cancers, especially for esophageal cancer [8, 9]. Although a few studies found that genetic polymorphisms in the ALDH2 were associated with colorectal cancer risk [1012], genetic variants in the alcohol metabolic genes have not been systematically explored for the susceptibility for this cancer.

In the present study, we searched for variants in the alcohol metabolic genes and performed a case-control study to test whether these variants are associated with CRC susceptibility. We found rs3741178 in ALDH3B2 was significantly associated with CRC risk. Further functional prediction found that this variant may damage the protein function of ALDH3B2.

2. Materials and Methods

2.1. Study Subjects

This study included 485 CRC patients and 516 healthy controls. The CRC cases were enrolled from Beijing Chao-Yang Hospital from January 1, 2016, to December 31, 2018. All cases were histopathologically or cytologically confirmed by at least two local pathologists according to the World Health Organization classification. Health controls were cancer-free individuals who lived in the same residential areas and were seeking for routine physical examination in the same time period from the same hospital where cases were collected. All participants were unrelated Han Chinese descent. The informed consent was obtained from every participant at recruitment, and peripheral blood samples and demographic characteristics such as gender, age, smoking status, drinking status, and ethnicity were collected by interviewers. This study was conducted under the approval of the Institutional Review Board of Beijing Chao-Yang Hospital, Capital Medical University.

2.2. SNP Selection and Genotyping

We searched for alcohol metabolic genes using Molecular Signature Database (http://www.broadinstitute.org/gsea/msigdb/index.jsp). We included only missense variants in alcohol metabolic genes (ADH and ALDH gene families) with minor allele frequencies (MAF) >0.01 in the Chinese Han population from the 1000 Genome Project. A total of 16 germline coding variants in 10 genes (ADH1C, ADH4, ADH7, ALDH1A1, ALDH1A2, ALDH1A3, ALDH1B1, ALDH3B1, ALDH3B2, and ALDH2) were selected and genotyped. Genotyping was performed using genomic DNA extracted from peripheral blood sample collected from each participant at recruitment. The DNA was extracted using TIANamp® Genomic DNA kit (Tiangen Biotech). SNPs were genotyped by the Sequenom MassARRAY system (San Diego, CA, USA). The case and control samples were mixed in the plates, and persons who performed the genotyping assay were not aware of case or control status.

2.3. Statistical Analysis

The association between SNPs and CRC risk was tested by using logistic regression analysis with adjustment of sex and age. Statistical analyses were performed using SPSS software (18.0). All tests were two-sided.

2.4. Functional Annotation Analysis

Multiple functional annotation tools were used to predict the potential function of rs3741178, including PolyPhen-2 [13], SIFT [14], and AWESOME [15].

3. Results

3.1. Characteristics of Study Subjects

The distributions of selected demographic characteristics including sex and age of the CRC patients and healthy controls are summarized in Table 1. This study consisted of 485 CRC cases and 516 controls. The average age of cases and controls was 61.23 (±12.47) and 61.68 (±11.81), respectively. There are 357 (73.6%) males and 128 (26.4%) females in the cases, while there are 358 (69.4%) males and 158 (30.6%) females in controls. Among cases, 188 (38.8%) and 138 (28.5%) were smokers and drinkers, respectively.

3.2. rs3741178 in ALDH3B2 Was Significantly Associated with CRC Risk

A total of 16 germline coding variants in the alcohol metabolic genes were genotyped, and their association with CRC risk was tested (Table 2). Among these variants, only rs3741178 in ALDH3B2 was significantly associated with increased risk of CRC with FDR <0.25 (Table 2). There were no significant gene-drinking interactions for these variants (Table 2). Compared with CC genotype carriers, CT genotype carriers were associated with CRC risk with odds ratio (OR) being 2.29 (95% CI: 1.26–4.14, ) (Table 3). When combing the CT genotype and TT genotype carriers, the association was also significant with values being 0.0050 (Table 3). There were no significant differences in the stratification analysis for this variant (Table 4).

3.3. rs3741178 in ALDH3B2 Was Predicted to Damage the Protein Function

The rs3741178 variant was located in the 2nd exon of ALDH3B2 with Ala > Thr change. To test the potential function of rs3741178, we used multiple coding variant prediction tools. PolyPhen-2 showed that this variant was predicted to be “PROBABLY DAMAGING” with a score of 1.000, based on a number of features comprising the sequence, phylogenetic, and structural information characterizing the substitution. The SIFT result showed that this variant was predicted to be deleterious with the score being 0.009, indicating that the amino acid (rs3741178 Ala) was highly conserved in the protein family. The AWESOME showed that no significant posttranslational modification was affected by rs3741178. These results suggested that rs3741178 may damage the protein function of ALDH3B2 and thus associated with risk of colorectal cancer.

4. Discussion

In this case-control study, we explored the association between variants in the alcohol metabolic genes and colorectal cancer risk in 485 CRC cases and 516 controls. We found only rs3741178 in ALDH3B2 were significantly associated with an increased risk of CRC. Further functional prediction showed that this variant may damage the protein function of ALDH3B2. These results indicated an important role of ALDH3B2 in the CRC carcinogenesis.

ALDH3B2, also known as ALDH8, was a member of the aldehyde dehydrogenase family. This gene family, including ALDH1A1, ALDH1A2, ALDH1A3, ALDH2, ALDH3A1, ALDH3A2, ALDH3B1, and ALDH3B2, plays a major role in the detoxification of aldehydes generated by alcohol metabolism. ALDH1 and ALDH2 are the most important enzymes for aldehyde oxidation. ALDH1A1 and ALDH1B1 were shown to be markers for colorectal cancer [16, 17]. ALDH1A3 was reported to affect colon cancer proliferation and invasion [18]. The genetic polymorphisms in ALDH2 were reported to be associated with colorectal cancer [1012].

However, little is known about ALDH3B2 in the development of cancer, especially CRC before this study. It was reported that DNA methylation of ALDH3B2 was associated with alcohol dependence [19]. Polymorphism in ALDH3B2 was reported to be associated with esophageal squamous cell carcinoma in a Chinese population [20]. In this study, we found that the polymorphism of ALDH3B2 was associated with CRC risk also in a Chinese population. These results indicated that ALDH3B2 may play an important role in the carcinogenesis of CRC by altering the alcohol metabolism process. The function of this gene in the development of CRC is worth to be investigated. Furthermore, rs3741178 variant may serve as a potential marker for the early detection of CRC.

There are also some limitations for this study. First, our sample size has only 24%–78% power to identify variants with OR being 2.0 with MAF = 0.01 to 0.05, indicating that the statistical power of this study may be insufficient. The results need to be validated in larger samples in the future. Second, the sample size was too small to perform a gene-drinking interaction with sufficient power. The interaction between these variants and drinking in the susceptibility of CRC needs to be investigated in the future. Finally, the function of ALDH3B2 variant was only predicted in this study. In vitro and in vivo experiments need to be conducted to validate the prediction results.

In summary, through a case-control study in a Chinese Han population, we find a significant association between the coding variant in ALDH3B2 and CRC risk. These results expand our insights of CRC carcinogenesis and provide more evidence for the precision medicine of this disease.

Data Availability

The data used to support the findings of this study are available from the corresponding author upon request.

Conflicts of Interest

The authors declare no conflicts of interest.

Authors’ Contributions

Zhi-Gang Gao was the overall principal investigator of this study, who obtained financial support of the project and drafted the manuscript. All authors performed the laboratory analyses and are responsible for patients’ recruitment.

Acknowledgments

This work was supported in part by grants from Beijing Nova Program (grant no. Z171100001117015), Beijing Talents Foundation (grant no. 2017000021223ZK38), and 1351 Talents Program of the Beijing Chao-Yang Hospital (grant no. CYXX-2017-35).