Abstract

Our study intended to investigate five cytokine gene single nucleotide polymorphisms (SNPs) and their associations with prostate cancer risk. Genotypes of five cytokine gene SNPs were detected by MassARRAY for blood samples from a group of patients with prostate cancer (n = 90) and a control group (n = 140) in central China. The differences in tumor clinical stages, Gleason scores, and PSA values in patients with prostate cancer were also investigated. The frequencies of the five cytokine gene SNPs (L-1β rs16944, IL-4 rs2070874, IL-4rs2227284, IL-16 rs7175701, and IL-16 rs11556218) genotypes were not found to be significantly mutated in prostate cancer patients compared with the control group. In addition, for five cytokine gene SNPs genotypic comparisons, patients with different Gleason scores, clinical stages, and PSA values were grouped into two subgroups. There was also no statistically significant association in all these subgroups. Our study suggests that cytokine gene polymorphisms may not be a risk factor for prostate cancer in a central Chinese population. Nevertheless, more large-scale studies on the Chinese population are necessary to examine our conclusions. The discovery of cytokine gene polymorphisms related to prostate cancer could update our understanding of the etiology and improve our knowledge of the early detection, diagnosis, and treatment of prostate cancer.

1. Introduction

Prostate cancer (PCa) is the most common urinary neoplasm in men, which is estimated to account for about 21% of all cases of newly diagnosed cancers and account for 33,330 cases of deaths in males in America [1]. Although the incidence of prostate cancer in Asian populations is low, the incidence continues to rise. Like other cancers, the pathogenesis and tumorigenesis of prostate cancer are still not clear to the public [2]. Up till now, many factors, including ethnicity, age, or family history, have been identified to affect the incidence of prostate cancer [3]. In addition, lifestyle behaviors [46] and even genetic factors [710] were also regarded as risk factors for prostate cancer. Furthermore, studies have demonstrated that the morbidity of PCa is gene-related at a level of approximately 57% [11], and of those, inflammation-related genes are extensively studied. The diagnosis of prostate cancer is based on a microscopic evaluation of prostate tissue obtained through a needle biopsy. Conventionally, a systematic prostate biopsy is performed using transrectal ultrasound to obtain 10 to 12 gridded tissue samples. Pathologists examine these samples and assign a primary Gleason grade for the predominant histological pattern and a secondary grade for the highest pattern, both on a scale of 1 to 5, based on the microstructure and appearance of the cells. Clinicians have traditionally classified diagnoses as low-risk, intermediate-risk, and high-risk based on the sum of the Gleason pattern, prostate-specific antigen (PSA) level, and clinical stage [12]. Therapeutic methods for prostate cancer, in which surgery and radiation therapy remain treatments for localized disease, but have side effects such as urinary symptoms and sexual dysfunction that negatively impact the quality of life. For metastatic disease, chemotherapy as initial therapy now appears to prolong survival compared with androgen deprivation therapy alone. New vaccines, hormone therapy, and bone-targeting agents have been shown to be effective in men with metastatic prostate cancer resistant to conventional hormone therapy [12, 13].

Infections and inflammatory responses have been reported to be associated with approximately 15% of cancer deaths globally. Chronic inflammation was considered to have an important role in the pathogenetic process of cancers [14, 15] and is now relevant to the early steps of oncogenesis and cancer progression [16, 17]. Inflammation of the prostate was reported to participate in the occurrence and development of prostate cancer [18, 19]. The features of cancer-related inflammation include the infiltration of white blood cells and the production of inflammatory mediators accompanying chronic inflammation, which is often observed in the process of tissue reconstruction and repair [16]. One of the most important components of the inflammatory mediators is the production of cytokines associated with the host immune response. Therefore, genetic variability in altering gene expressions within key cytokines may be related to prostate cancer risk. Cytokine gene single nucleotide polymorphisms (SNPs) of interleukin (IL)-1β, IL-4, and IL-16 genes have been reported to be associated with the risk of cancers, including renal cell carcinoma, gastric carcinoma, and lung cancer [2022]. For PCa, some studies reported that interleukin gene polymorphisms may contribute to the occurrence of prostate cancer [2325]. For example, Luo et al. reported that IL-4 SNPs in PCa were related to cancer risk or mortality [26]. Shailendra Dwivedi et al. concluded that genetic variants of IL-18 and IL-10 promoters influence the circulating levels of these interleukins. Variations at −137 and −607 loci of IL-18 are associated with susceptibility to PCa [27]. Prostate tumor cells produce large amounts of IL-6 and its receptor, IL-6R (gp80), and IL-6 functions as a paracrine growth factor for the human LNCaP androgen-sensitive prostate cancer cells and an autocrine growth factor for the human DU145 and PC3 androgen-insensitive prostate cancer cells [2832].

However, some other studies reported that there is no association observed between interleukin gene polymorphisms and prostate cancer [33, 34]. Thus, the association between candidate cytokine gene SNPs and prostate cancer remains to be an unsettled issue.

In this study, we intend to investigate the association between the cytokine gene SNPs genotypes and allele distributions and prostate cancer risk. IL-1β and IL-16 gene SNPs were designed for their role in promoting inflammation, and IL-4 gene SNPs were designed for their role in suppressing inflammation. These polymorphisms have been studied before with conflicting results. However, few studies have been conducted on Asian populations, especially in China. Therefore, to figure out the association between the cytokine gene SNPs and prostate cancer risk, we performed this hospital-based case-control study to validate the risk of five cytokine gene SNPs on prostate cancer.

2. Materials and Methods

2.1. Study Subjects

This study was performed according to the Declaration of Helsinki and approved by the Ethics Committee of the first affiliated hospital of Anhui Medical University, Hefei, Anhui, China. All patients or participants signed written informed consent to allow analysis of data for research use. In this study, we adopted a value of Cohen’s d of 0.7, a statistical power level of 0.8, and a p value as 0.05. After calculation, the minimum total sample size was 52 in the evaluation and a minimum sample size of 26 for each group. Totally, between September 2014 and July 2016, ninety patients newly diagnosed with prostate cancer who underwent prostate biopsy were included in this study. Patients with prostate cancer were subdivided into two groups according to Gleason score (low grade ≤ 7 and high grade > 7), clinical stages (based on the international tumor node-metastasis system for PCa, localized and advanced), and PSA values (PSA ≤ 20 ng/ml and PSA > 20 ng/ml). Correspondingly, 140 volunteers were recruited into a healthy control group from individuals who came to the hospital for a routine health examination. All of the healthy people were included according to the following criteria: normal physical examination indicators; no genetic family history; no diagnosis of acute infectious diseases; no drinking or smoking habits; no history of drug abuse; and no impairment of the liver, heart, or kidney. All of the PCa patients were included according to the following criteria: no history of cancer or prostate surgery; no diagnosis of acute infectious diseases; no therapeutic medication for PCa; and no impairment of the liver, heart, or kidney. Venous blood samples were collected from all the participants after fasting for at least 12 hours. We also gathered related clinical and pathological data, including body mass index (BMI), systolic blood pressure (SBP), and serum PSA values from the medical record.

2.2. Selection of SNPs for Evaluation and Genotyping

The total DNA of all patients and controls was extracted from venous blood samples by using commercial DNA extraction kits (Qiagen, # 51106). In the study, we selected five cytokine gene SNPs through the SNP Consortium database (https://snp.cshl.org/) and the DbSNP database (https://www.hapmap.org/index.html.en) for analysis, and each SNP had a minor allele frequency (MAF) of >5% in the Chinese Han population. The five cytokine gene SNPs included L-1β rs16944, IL-4 rs2070874, IL-4rs2227284, IL-16 rs7175701, and IL-16 rs11556218. For each DNA sample, five cytokine gene polymorphisms were detected by MassARRAY (Sequenom, San Diego, USA). Details of the experiment process are referred to in Qin et al.’s [35] study. The data were analyzed by TYPER 4.0 software (Sequenom).

2.3. Statistical Analysis

The Hardy-Weinberg equilibrium (HWE) of cases and controls were calculated by chi-square analysis. The cancer risk, which was assessed by ORs with 95% CI, was calculated by using logistic regression models; otherwise, calculation results were adjusted for age and BMI. was identified to be statistically significant. The differences in tumor clinical stages, Gleason scores, and PSA values in patients with prostate cancer were also investigated by using the statistical methods mentioned above. The software used for statistical analyses in the study was IBM SPSS Statistics 20.0 (IBM Corp., Armonk, NY, USA).

3. Results

3.1. Characteristics of the Study Subjects

Table 1 showed the detailed demographic characteristics of study subjects. The mean age of PCa patients was 73.00 years old, ranging from 48 to 87 years. The mean age of the controls was 67.66 years old, ranging from 60 to 95 years. Between the two groups, we found a significant difference in age (). According to the World Health Organization recommendations for Asian populations, the normal BMI range is 18.5–22.9 kg/m2. There were no significant differences in BMI and systolic blood pressure between the PCa group and the control group. PCa patients were divided into two groups by PSA, Gleason score, and tumor clinical stage, respectively.

3.2. Association of Five Cytokine Gene SNPs with Prostate Cancer Risk

All five cytokine gene SNPs were in HWE among the PCa group and the control group (). The genotypic distributions and allele frequencies of cytokine gene SNPs between the two groups are described in Table 2. Stratification analyses of genotypes on risk of PCa were performed by age, BMI, and SBP groups (Table 3). Then, we performed a comparative analysis of genotypes in PCa patients. These patients were subdivided into two groups according to different Gleason scores (Gleason score ≤ 7, low-grade group; Gleason score > 7, high-grade group), TNM staging (localized group and advanced group), and PSA levels (PSA ≤ 20 group; PSA > 20 group), respectively (Table 4). The details of the results of five cytokine gene SNPs are listed below:

3.3. IL-1β rs16944

The IL-1β rs16944 genotypic distributions of the control group complied with HWE (). There was no difference between the PCa patients and the control group in the frequencies of rs16944 genotypes, including GG, GA, and AA (Table 2). Compared with GG, the variant GA and AA genotypes were not associated with the risk of prostate cancer after an adjustment of age and BMI (GA vs. GG: adjusted OR = 0.689, 95% CI (0.335, 1.418), ; AA vs. GG: adjusted OR = 0.734, 95% CI (0.316, 1.701), ; GA + AA vs. GG: adjusted OR = 0.705, 95% CI (0.361, 1.379), ). When analyzing frequencies of haplotypes, compared to G allele, the variant A allele was also not associated with PCa risk (). In addition, stratification analyses showed that when grouped by age, compared to GG genotypes, the variant AG () and AG + AA () genotypes showed a protective role for prostate cancer in the age > 70 group. There was also no significant difference in the stratified analysis results of different BMI groups and SBP groups (Table 3).

3.4. IL-4 rs2070874

The IL-4 rs2070874 genotypic distributions of the control group complied with HWE ( ). The distributions of rs2070874 genotypes were similar between the PCa group and the control group (Table 2). Compared with TT, the variant TC and CC genotype were not associated with the risk of prostate cancer after adjusting for age and BMI (TC vs. TT: adjusted OR = 1.159, 95% CI (0.637, 2.109), ; CC vs. TT: adjusted OR = 0.314, 95% CI (0.034, 2.896), ; TC + CC vs. TT: adjusted OR = 1.072, 95% CI (0.597, 1.924), ). When analyzing frequencies of haplotypes, compared to T allele, the variant C allele was also not a risk for prostate cancer (). In addition, stratification analysis showed that there was no significant difference in age, BMI, and SBP groups (Table 3).

3.5. IL-4 rs2227284

The IL-4 rs2227284 genotypic distributions of the control group complied with HWE (). The distributions of rs2227284 genotypes were similar between the PCa group and the control group (Table 2). Compared with TT genotypes, the variant TG and GG genotypes do not affect the risk of prostate cancer after an adjustment of age and BMI(TG vs. TT: adjusted OR = 0.669, 95% CI: (0.343, 1.304), ; GG vs. TT: adjusted OR = 2.718, 95% CI: (0.155,47.524), ; TG + GG vs. TT: adjusted OR = 0.703, 95% CI: (0.365,1.355), ). As for the alleles, compared to T allele, the variant G allele does not influence the risk of prostate cancer (). In addition, stratified analyses showed that when grouped by age, compared with TT, the variant TG () genotypes could decrease the risk of prostate cancer in age ≤ 70 group. There was also no significant difference in the stratified analysis results of different BMI groups and SBP groups (Table 3).

3.6. IL-16 rs7175701

The IL-16 rs7175701 genotypic distributions of the control group complied with HWE (). The distributions of rs7175701 genotypes were similar between the PCa group and the control group (Table 2). Compared with TT, the variant TC genotype does not affect the risk of prostate cancer after an adjustment of age and BMI (TC vs. TT: adjusted OR = 0.224, 95% CI (0.044, 1.143), ). As for the alleles, compared to T allele, the variant C allele does not influence the risk of prostate cancer (). Besides, stratification analysis showed that there was no significant difference in age, BMI, and SBP groups (Table 3).

3.7. IL-16 rs11556218

The IL-16 rs11556218 genotypic distributions of the control group complied with HWE (). The distributions of rs11556218 genotypes were similar between the PCa group and the control group (Table 2). Compared to TT, the variant TG and GG genotypes do not affect the risk of prostate cancer after an adjustment of age and BMI (TG vs. TT: adjusted OR = 0.507), 95% CI: (0.264, 0.976), ; GG vs. TT: adjusted OR = 0.598, 95% CI: (0.138, 2.588), ; TG + GG vs. TT: adjusted OR = 0.525, 95% CI: (0.281, 0.978), . As for the alleles, compared to T allele, the variant G allele showed no association with the risk of prostate cancer (). In addition, stratified analyses showed that when grouped by age, compared with TT, the variant TG () genotypes could decrease the risk of prostate cancer in age > 70 group. No significant differences were also found by stratification analysis when performed by BMI and SBP groups (Table 3).

Distributions of the five cytokine gene SNPs genotypes in patients with prostate cancer were also analyzed. PCa patients with different Gleason scores, clinical stages, and PSA values were subdivided into two groups. As a result, there were no significant differences in the five cytokine gene SNPs genotypic frequencies by Gleason scores, clinical stages, and PSA levels (Table 4).

4. Discussion

Like many other tumors, the cause of prostate cancer is unknown. Numerous factors have been reported to increase the incidence of cancer. In recent years, many studies have reported that there may be a link between gene polymorphisms and PCa risks. From then on, subsequent studies have demonstrated SNPs could affect the occurrence of PCa [36, 37]. Interestingly, inflammation has been suggested as a potential trigger for the development of prostate cancer [38]. Epidemiological studies have demonstrated this viewpoint, and the use of anti-inflammatory agents could decrease the risk of prostate cancer [39]. Thus, we hold the opinion that the variants in relevant genes encoding certain proinflammatory or anti-inflammatory cytokines may affect prostate cancer risk since cytokine gene SNPs could lead to aberrant expression of cytokines [40, 41].

IL-1β is a 17-kDa cytokine member of the interleukin-1 family, which is involved in inflammation and immunity. The IL-1β gene is located on chromosome 2q14, and its expression has been reported to be related to different diseases [42, 43]. IL-1β gene polymorphisms are associated with increased secretion of cytokine IL-1β when compared to the wild-type gene [40]. In recent years, it has been reported that IL-1β gene polymorphisms may be related to various cancer types, including gastric carcinoma [21], hepatic carcinoma [44], and lung cancer [45].

IL-4 is a member of the α-helical cytokine family, which is secreted by activated CD4+ T cells, basophils, and mast cells. IL-4 is an important component of our immune system and is involved in functions such as inhibition of the production of proinflammatory cytokines by monocytes and immunoglobulin class switching in activated B lymphocytes. Thus, IL-4 plays a vital role in the inhibition of inflammation and the surveillance and elimination of transformed cells. Many epidemiologic studies have reported the association of IL-4 gene SNPs with cancer risks, such as breast cancer [46], gastric cancer [47], and lung cancer [22].

IL-16 is a proinflammatory cytokine that has a wide array of functions in multiple immunopathobiological processes [48]. IL-16 stimulates the production of a variety of proinflammatory cytokines by monocytes with related biological sequelae [49]. Evidence of increased IL-16 levels in serum or plasma has also been recognized in multiple malignant tumors in both preclinical and clinical trials [5052]. Studies have reported a relationship between IL-16 gene polymorphisms and cancers, including gastric cancer [53], renal cell cancer [20], and colorectal cancer [54].

For PCa, previous studies reported that IL-1β, IL-4, and IL-16 gene polymorphisms could affect the incidence of prostate cancer [24, 25, 55]. The five cytokine gene SNPs we selected in this study were hardly investigated, especially in China. In the present study conducted in central China, we found that all these five cytokine-related loci, including IL-1β rs16944, IL-4 rs2070874, IL-4 rs2227284, IL-16 rs7175701, and IL-16 rs11556218 were all not significantly associated with the occurrence of prostate cancer. The results of the present study were inconsistent with previous studies reported on the Turkish and African American populations. Faruk Yencilek et al. showed that individuals with the IL1B-31 (rs1143627) and IL1B-511 (rs16944) AG genotypes have a reduced risk of prostate cancer. Of all possible combinations analyzed, the combination of IL1B-31 (rs1143627) G and IL1B-511 (rs16944) G prevented the development of prostate cancer. To the best of our knowledge, the data obtained in this study are the first to highlight the association of IL1B gene variants with prostate cancer in the Turkish population [55]. Meanwhile, in prostatectomy tumor specimens of prostate cancer patients, the expression of IL-16 is positively correlated with Gleason score and pathological stage, and IL-16 is also highly expressed in other types of tumor cells. The level of IL-16 in serum is also elevated in patients with other types of cancer, especially in patients with advanced cancer. Ken Batai et al. provide evidence that IL-16 polymorphisms show great significance in Pca susceptibility in African Americans. To our knowledge, this is the first study that demonstrated the association of these variants with cancer in African-descent populations [25]. However, our study suggests that the effect of the five cytokine gene SNPs on PCa may have a racial difference in different populations, though further investigations will be needed.

Additionally, we conducted the analysis of genotypes comparison in PCa patients. These patients were subdivided into two groups according to different Gleason scores (Gleason score ≤ 7, low-grade group; Gleason score > 7, high-grade group), TNM staging (localized group and advanced group), and PSA levels (PSA ≤ 20 group; PSA > 20 group), respectively. As a result, there was no significant association observed between the cytokine gene polymorphisms and all these subgroups of prostate cancer. Our study suggested that the five cytokine gene SNPs may not be related to cancer progression in prostate cancer.

However, there were some limitations to our study. Firstly, although our participants were recruited from a relatively homogenous Chinese population in Anhui province, a central area of China, we could not thoroughly eliminate the influence of potential genetic background. Secondly, for the case-control study, the sample size may be a limitation. Thirdly, the lack of an independent set of samples to examine our conclusions is another limitation. Thus, we plan to carry out further optimization studies and pay attention to the following issues: first of all, when we screen patients, we investigate their genetic backgrounds to ensure that they are as similar as possible. Second, the age difference between the control group and the case group should be minimized. Finally, select a few provinces in central China, not just Anhui, so that it is representative and the sample size is large enough.

Thus, further large-sample studies should be conducted on Chinese populations to evaluate our results. In addition, the age difference between the cases and control group may affect the results of our study.

Bioinformatics is an emerging discipline that emerged with the launch of the Human Genome Project (HGP), which integrates mathematics, computer science, and biology to elucidate the biological significance of various types of data. Currently, bioinformatics plays a pivotal role in the development of medicine. First, based on different omics datasets such as transcriptome, proteome, and epigenome, using clustering, consensus, and other ideas to achieve disease classification is of great significance for understanding the disease mechanism. [5659] In addition, bioinformatics can be used to predict and analyze gene variants, expression, as well as gene and protein structure and function, which are crucial for disease diagnosis and treatment. Therefore, in the future, we should combine immunological therapy with bioinformatics to deeply research the pathology of diseases.

5. Conclusion

In summary, the present study suggests that cytokine gene polymorphisms, including IL-1β rs16944, IL-4 rs2070874, IL-4 rs2227284, IL-16 rs7175701, and IL-16 rs11556218 may not be risk factors for prostate cancer in a population in central China. Nevertheless, more large-scale studies on the Chinese population are necessary to examine our conclusions. The discovery of prostate cancer-related cytokine gene polymorphisms could update our understanding of the etiology and improve our knowledge of the early detection, diagnosis, and treatment of prostate cancer.

Data Availability

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

Ethical Approval

This study was approved by the Ethics Committee of the first affiliated hospital of Anhui Medical University.

All participants had signed the written informed consent.

Conflicts of Interest

The authors declare that they have no conflicts of interest.

Authors’ Contributions

SX is responsible for the study concepts and design and manuscript editing. YTH is responsible for the experimental studies and manuscript preparation. LGZ is responsible for the literature research and manuscript editing. SYY is responsible for the clinical and experimental studies. LZ is responsible for data acquisition. SX is responsible for the statistical analysis and the manuscript preparation. LZ is responsible for the definition of intellectual content and data analysis. CZL is responsible for the integrity of guarantee of the entire study, study concepts and design, and manuscript review. All authors read and approved the final manuscript. Shun Xu and Yongtao Hu contributed equally to this work.

Acknowledgments

This work was supported bythe Project of National Natural Science Foundation of China [NO.82100814],National Natural Science Foundation of China [No.82170787], and Anhui MedicalUniversity Scientific Research Fund Project [No.2021xkj157].