Abstract

Background. Chronic prostatitis/chronic pelvic pain syndrome (CP/CPPS) refers to a common disorder with unclear etiology and unsatisfactory treatment, which reduces the male’s quality of life. Objective. To examine the effects of genetic polymorphisms of IFNG, IFNGR1, and androgen receptor (AR) on CP/CPPS. Methods. The single nucleotide polymorphisms (SNPs) of IFNG, IFNGR1, and AR were genotyped with the improved multiplex ligation detection reaction. The GTEx, RegulomeDB, HaploReg, and 3DSNP databases were adopted to predict the regulatory functions of the genotyped SNPs. The correlation between SNPs and CP/CPPS was analyzed with the test, logistic regression, and two genetic models (codominant and log-additive models). The nomogram was built to predict the risk of CP/CPPS occurrence. Results. On the whole, 130 CP/CPPS patients and 125 healthy controls were recruited in the study, and 18 SNPs of IFNG, IFNGR1, and AR were genotyped. The results of functional annotation indicated that the 18 genotyped SNPs might have regulatory effects in the whole blood. The rs144488434 was correlated with the elevated CP/CPPS risk (odds ratio (OR): 2.40, 95% confidence interval (CI): 1.12-5.13, , and ) by the test. In the built genetic models, rs10457655 was correlated with the elevated National Institutes of Health Chronic Prostatitis Symptom Index (NIH-CPSI) scores (codominant model: GA/GG: , 95% CI: -1.71-3.67 and AA/GG: , 95% CI: 0.58-17.62, ). In subgroup analysis, rs2069718 was correlated with the elevated CP/CPPS risk (log-additive model: , 95% CI: 1.03-4.64, and ) in years. The nomogram integrating age, rs2069718, rs10457655, and rs144488434 showed good performance to predict the risk of CP/CPPS. Conclusions. Genetic polymorphisms of IFNG, IFNGR1, and AR might act as the genetic factors for CP/CPPS susceptibility, which deserved further explorations.

1. Introduction

According to the National Institutes of Health (NIH) classification, chronic prostatitis/chronic pelvic pain syndrome (CP/CPPS) is defined as category III prostatitis, affecting approximately 8.4% to 25% of males worldwide [1, 2]. The clinical manifestations of CP/CPPS vary significantly, which consist of pain or discomfort in the perineal or pelvic region, sexual dysfunction, and others, thereby affecting males’ quality of life (QoL) [3]. As the etiology of CP/CPPS has been extensively investigated, genetic predisposition towards CP/CPPS was proposed. By analyzing genetic polymorphisms of tumor necrosis factor- (TNF-) α and interleukin- (IL-) 10, Shoskes et al. reported that CPPS patients had a higher proportion of low IL-10 production genotype, and the symptoms of patients with high IL-10 and low TNF-α phenotypes were not mitigated by anti-inflammatory phytotherapy [4]. Additionally, polymorphisms of manganase superoxide dismutase and genes in the X chromosome might be the genetic factors correlated with CP/CPPS development [5, 6]. Hence, genetic polymorphism might contribute to the pathogenesis of CP/CPPS, and a genome-wide association study (GWAS) should be conducted to identify CP/CPPS-related high-risk alleles and genotypes to facilitate its diagnosis and treatment.

Interferon-γ (IFN-γ, encoded by IFNG) and its receptor (encoded by IFNGR1) exert vital effects on the inflammatory response against pathogen infection [7]. According to the experimental autoimmune prostatitis (EAP) model [8], the IFN-γ level increased in the prostate, and IFN-γ-secreting cells specific to the prostate were detected in the spleen, thereby determining the key effects of IFN-γ on the induction of autoimmune prostatitis [9]. For CP/CPPS patients, IFN-γ levels were upregulated in the seminal plasma, and prostate antigen-specific IFN-γ-secreting lymphocytes were identified in the peripheral blood [10]. In our previous studies, the proportion of IFN-γ+-T-helper 1 (Th1) cells was reported to increase in the periphery blood of CP/CPPS patients [11]. Thus, IFN-γ significantly affected the pathogenesis of CP/CPPS. Several single nucleotide polymorphisms (SNPs) of IFNG and IFNGR1 (e.g., rs3799488, rs9376267, and rs9376268) were correlated with the risk of pulmonary tuberculosis and colon and rectal cancer [1214]. Genetic polymorphisms of AR were associated with the risk of testicular cancer and male infertility [15, 16]. Riley and Krieger highlighted the genes located in the X chromosome (e.g., the androgen receptor (AR)) should be studied in depth to examine the relationship between genetic polymorphisms and CP/CPPS susceptibility [6]. However, the study devoted to elucidating the correlation between polymorphisms of IFNG, IFNGR1, and AR and CP/CPPS has been rarely conducted, and the effects of genetic polymorphism in CP/CPPS deserved further exploration to gain insights into CP/CPPS development.

In the present study, the SNPs of IFNG, IFNGR1, and AR and CP/CPPS risk were analyzed, and a nomogram was built to predict the risk of CP/CPPS based on age and the identified SNPs, which exhibited high performance.

2. Materials and Methods

2.1. Participants and Samples

All patients diagnosed as CP/CPPS originated from the urology clinic of the First Affiliated Hospital of Anhui Medical University between September 2019 and January 2020. The ethnicity and residency-matched healthy subjects came from the healthy check-up center of our hospital. All individuals were Han Chinese and lived in Anhui province. After the history was critically reviewed and the physical examination was performed, CP/CPPS patients were diagnosed by a senior urologist based on the NIH classification [1, 17]. All patients completed the National Institutes of Health Chronic Prostatitis Symptom Index (NIH-CPSI) questionnaire [18]. The inclusion and exclusion criteria are presented below: the inclusion criteria: (i) patients suffering ≥3 months of symptoms in the last 6 months, covering pain or discomfort in the perineal, lower abdominal, or pelvic region, with the and (ii) normal results of expressed prostatic secretion (EPS) and urine analyses. The exclusion criteria: patients with other urological disorders consisting of acute or chronic bacterial prostatitis, urinary tract infection or urethritis in the last year, interstitial cystitis, urolithiasis, history of urogenital cancer, and prostate surgery history, as well as neurological disorders. All patients recruited provided the written informed consents, and this study was approved by the ethics committee of the First Affiliated Hospital of Anhui Medical University (PJ-2021-08-45).

2.2. Genotyping

Genomic DNA was isolated from the 5 ml of collected venous blood and stored at -80°C to be genotyped. On the whole, 18 SNPs (2 SNPs of IFNG, 6 SNPs of IFNGR1, and 10 SNPs of AR) were genotyped with the improved multiplex ligation detection reaction (iMLDR) system by adopting the multiplex polymerase chain reaction- (PCR-) ligase detection method (Shanghai Genesky Biotechnologies Inc., China). The different fluorescent labels (FAM, VIC, NED, and PET) and different extended 3-end lengths were exploited to distinguish alleles and SNPs, respectively. The double-distilled water-replaced templates and the primer-free samples acted as the two negative controls. Duplicate tests were performed to verify the consistency of the results. The GeneMapper 4.0 software (Applied Biosystems, America) was adopted to process the raw data.

2.3. Functional Annotation of the 18 Genotyped SNPs

To assess the effects of genetic variants of IFNG, IFNGR1, and AR in CP/CPPS-related symptoms, the correlation between genotypes and NIH-CPSI scores was analyzed. The GTEx database [19] was adopted to explore the effects of genetic variants on the expressions of their corresponding genes in whole blood. The HaploReg [20], RegulomeDB [21], and 3DSNP [22] were used to examine the potential effects of the 18 genotyped SNPs in our study. The HaploReg v4.1 online tool is aimed at annotating the noncoding genome by predicting the potential target genes, relevant variants in the haplotype block, etc. Moreover, the annotation of the 18 genotyped SNPs in HaploReg was performed by setting the parameters below: East Asian population, , and the 15-state core model source for epigenome analysis. The 3DSNP database was employed to predict the genetically related SNPs by three-dimensional (3D) chromatin looping.

2.4. ELISA Assay

The enzyme-linked immunosorbent assay (ELISA) assay kit (cat# E-EL-H0108c; Elabscience Biotechnology) was used to detect the serum IFN-γ concentration in CP/CPPS patients and healthy controls at 450 nm, which was performed by complying with the manufacturer’s manual.

2.5. Statistical Analysis

For continuous variables, data were expressed as . The correlation between allele and genotype frequencies and CP/CPPS was assessed by the test and the logistic regression to establish the genetic model with the SPSS (version 22.0, IBM Corp., Armonk, NY, USA) and the SNPStats tool [23]. The additive model performed better in detecting dominant and additive effects in the unknown inheritance pattern of the variant, whereas it remained poor in examining the recessive effects [24, 25]. Lettre et al. found that the codominant model outperformed other models [26]. Hence, we used the additive model and codominant model to examine the correlation between the genotyped SNPs and CP/CPPS. The two genetic models were defined as follows: codominant model: XX versus YX and XX versus YY and log-additive model: for each Y allele increase (wild allele: X and minor allele: Y). The Haploview (version 4.2) was used to perform the Hardy-Weinberg equilibrium (HWE) test and visualize the linkage disequilibrium (LD) among SNPs. The SHEsis software platform was adopted to perform haplotype construction [27], and the lowest frequency threshold was set to 0.03. The nomogram was built with the rms, rmda, and pROC packages in the R 3.6.3 software.

3. Results

3.1. Basic Information

A total of 130 patients diagnosed with CP/CPPS and 125 healthy controls were recruited in this study. The age was years in the patient group and years in the control group; the NIH-CPSI score was in CP/CPPS patients (Table 1). The information of the 18 SNPs, minor allele frequency (MAF), and HWE test are listed in Supplemental Table 1. Because all subjects in the present study were male, HWE was not assessed in SNPs located in the X chromosome. The primers applied in this study are listed in Supplemental Table 2.

3.2. Potential Effects of the 18 Genotyped SNPs

According to Supplemental Figure 1, no significant correlation was identified between genotypes of the 18 genotyped SNPs and NIH-CPSI scores. The results of the GTEx database analysis revealed that rs1861493, rs2069718, rs7749390, rs9376267, and rs9376268 might affect the expression of the corresponding genes (Supplemental Figure 2). According to Supplemental Table 3, the RegulomeDB scores ranged from 2b to 7 of the 18 genotyped SNPs. According to RegulomeDB, the higher score represents a less regulatory function of a variant [21]. For rs10457655, the RegulomeDB score reached 2b, thereby demonstrating its abilities of TF binding, any motif, a DNase Footprint, and a DNase peak. rs3799488, with a score of 3a, had roles of transcription factor (TF) binding, any motif, and a DNase peak. For the rest of SNPs with 4 to 6 scores, only less binding evidence was reported. The HaploReg results revealed that most of the 18 genotyped SNPs were correlated with motif changed and histone markers, and some variants might act as expression quantitative trait loci (eQTL) to regulate the expression of related genes. Additionally, the levels of serum IFN-γ were also detected in CP/CPPS patients and healthy controls, and no significant difference in serum IFN-γ levels was found between CP/CPPS patients and healthy controls, and genotypes of rs1861493 and rs2069718 were not correlated with the serum IFN-γ levels in CP/CPPS patients (Supplemental Figure 3). Taken together, these 18 SNPs were predicted to have potential regulatory functions, while their effects on CP/CPPS warranted further study.

3.3. Distributions of Allele and Genotype Frequencies and Genetic Model

According to Table 2, as revealed from the results, rs144488434 was correlated with the elevated CP/CPPS risk (T/C, odds ratio (OR): 2.40, 95% confidence interval (CI): 1.12-5.13, , and ), and no significant correlation was reported between the remaining SNPs and CP/CPPS. The codominant and log-additive models were built to examine the correlation among SNPs, CP/CPPS, and NIH-CPSI scores. The Akaike’s Information Criterion (AIC) and Bayesian Information Criterion (BIC) values were calculated to determine the appropriate model. As indicated from our results, rs10457655 was correlated with the elevated NIH-CPSI scores in CP/CPPS patients without being adjusted by age (codominant model: GA/GG: , 95% CI: -1.71-3.67 and AA/GG: , 95% CI: 0.58-17.62, ) (Tables 3 and 4). Accordingly, rs10457655 might be correlated with CP/CPPS.

3.4. Subgroup Analysis

Based on the findings above, we conducted a subgroup analysis according to CP/CPPS patients’ age. As shown in Supplemental Table 4 and 5, in years ( and ), rs2069718 was correlated with the elevated CP/CPPS risk without being adjusted by age (log-additive model: , 95% CI: 1.03-4.64, and ), while no significant correlation was found in years ( and ).

3.5. LD and Haplotype Analysis

We analyzed the LDs among the 18 genotyped SNPs, and LD was defined as and . Strong LDs were found between rs9376267 and rs9376268 in IFNGR1 ( and ) and between rs192195540 and rs5965429 in AR ( and ) (Figures 1(a)–1(d)). There was no correlation between rs9376267-rs9376268 and rs192195540-rs5965429 genotype distribution frequencies in CP/CPPS patients and healthy controls (Supplemental Table 6). Additionally, we predicted that the rs1861493, rs2069718, rs9376267, rs9376268, and rs10457655 might exert regulatory effects through interacting with their genetically related SNPs (LD ) in the 3DSNP database (Supplemental Table 7-11).

3.6. Predictive Nomogram

To gain insights into the effects of SNPs on CP/CPPS, the SNP-based nomogram was built to predict the risk of CP/CPPS. Since rs2069718, rs10457655, and rs144488434 were identified to be possibly correlated with CP/CPPS, age, rs2069718, rs10457655, and rs144488434 were adopted to build the nomogram (Figure 2(a)). The calibration curve of the nomogram indicated the agreement of our cohort (Figure 2(b)), and the concordance index (C-index) was 0.62 (95% CI: 0.54-0.68), and the area under the curve (AUC) of receiver operating characteristics (ROC) reached 0.62 (95% CI: 0.55-0.69). Based on the decision curve and the clinical impact curve, the clinical usefulness of the nomogram of this study was visualized, i.e., the model showed good net benefit (Figures 2(c) and 2(d)). The internal validation was performed to assess the goodness-of-fit of the built nomogram by splitting our cohort into training and validation cohorts at different ratios, and the C-indexes and AUC values in training and validation cohorts showed adequate internal goodness-of-fit of this study’s model (Supplemental Table 12).

4. Discussion

In the present study, we investigated the correlation between SNPs and the risk of CP/CPPS, the results of this study were as follows: (1) functional annotation of the 18 genotyped SNPs predicted their regulatory effects in blood, while these SNPs were not correlated with NIH-CPSI score; (2) rs2069718 in IFNG, rs10457655 in IFNGR1, and rs144488434 in AR might be correlated with the increased risk of CP/CPPS; (3) LD was reported between rs9376267 and rs9376268, as well as between rs192195540 and rs5965429; (4) nomogram based on age, rs2069718, rs10457655, and rs144488434 was built to predict the risk of CP/CPPS, which showed good performance. Taken together, we found that genetic polymorphisms of IFNG, IFNGR1, and AR were correlated with CP/CPPS, which might indicate the genetic predisposition to CP/CPPS development and provided novel insight into understanding the etiology of CP/CPPS.

IFN-γ, a proinflammatory and immunomodulatory cytokine secreted by natural killer (NK), NKT, and T cells, acts on activation of the innate immune system against pathogen invasion by binding to its receptor [28]. Previous studies have demonstrated that genetic polymorphisms of IFNG and IFNGR1 were correlated with tuberculosis infection, systemic lupus erythematosus, etc. [12, 29, 30]. Based on the important effects of IFN-γ in EAP and CP/CPPS [911], we determined the effects of genetic polymorphisms of IFNG and IFNGR1 in CP/CPPS. Initially, functional annotation of the 8 genotyped SNPs of IFNG and IFNGR1 was performed by GTEx database, RegulomeDB, and HaploReg tools, which indicated that these 8 SNPs might have regulatory effects in human whole blood. Nevertheless, no significant association among the mentioned variants, NIH-CPSI scores, and serum IFN-γ levels was observed, and the serum IFN-γ levels were not elevated in CP/CPPS patients. Hence, the 8 genotyped SNPs of IFNG and IFNGR1 might not be correlated with CP/CPPS-related symptoms, and the regulatory effects of the 8 genotyped SNPs in CP/CPPS deserved further study. The results of codominant and log-additive models showed that rs2069718 in IFNG was correlated with the elevated CP/CPPS risk, and rs10457655 in IFNGR1 was correlated with increased NIH-CPSI score. Combined with the autoimmune characteristics of CP/CPPS and the roles of rs2069718 and rs10457655 in inflammatory diseases [3032], we speculated that rs2069718 and rs10457655 might be involved in the pathogenesis of CP/CPPS through other mechanisms as predicted by the HaploReg tool [20], including protein binding and regulatory motif alterations. We also found that rs1861493, rs2069718, rs9376267, rs9376268, and rs10457655 might interact with their genetically related SNPs by the 3DSNP database, and combined with the significant functions of the SNP-SNP interaction network in the pathogenesis of prostate cancer [33] and inflammatory diseases [34, 35], we speculated that rs2069718 and rs10457655 might function through the SNP-SNP interaction pattern in CP/CPPS.

AR (encode by AR gene) can act as a transcription factor and exert its role by regulating androgen-targeted genes [36]. Polymorphisms of AR affected hormone level, sperm parameter, and the growth of the prostate [37, 38]. In function annotation of the 10 genotyped SNPs of AR, these SNPs were predicted to exhibit the regulatory potential by RegumeDB and HaploReg. Reciprocally, these SNPs did not regulate blood AR expression in the GTEx database. Additionally, the mentioned variants in AR were found to be not associated with CP/CPPS scores. Thus, the variants of AR might not regulate the phenotype of CP/CPPS directly, but through indirect mechanisms. In this study, rs144488434 in AR was found to be correlated with the elevated CP/CPPS risk, which might be involved in the pathogenesis of CP/CPPS, and further studies should be conducted to examine the underlying mechanisms. Based on the significant findings of rs2069718, rs10457655, and rs144488434 in CP/CPPS, a nomogram was built, integrating age, and genotypes of rs2069718, rs10457655, and rs144488434 to predict the risk of CP/CPPS, which might apply the findings of the study to clinical practice and assisted clinicians to more effectively identify CP/CPPS patients.

As the common urological disorders in males, approximately 8.4% of subjects had prostatitis-like symptoms in Chinese in our previous study [39]. The risk factors of CP/CPPS have been investigated increasingly (e.g., age, physical inactivity, stress, smoking, and alcohol consumption) [40, 41]. Zhang et al. [42] reported that CP/CPPS is a very problem in Chinese males, with a lifetime prevalence of 25.3%, and some patients might have a good prognosis without treatment, whereas some might require symptom-orientated therapies to facilitate the prognosis. Though therapeutic approaches for CP/CPPS (e.g., anti-inflammatory, antimicrobial agents, and α-blockers) had been applied to clinical practice, monotherapy could not effectively mitigate the symptoms, and combination therapy was recommended [43]. The CP/CPPS patients responded differently to the treatments, and the underlying mechanisms remain unclear. The direct (e.g., medications and office visits) and indirect costs (e.g., wage loss) to treat CP/CPPS were substantial to patients, and the unsatisfactory therapeutic efficacy and characteristic of recurrence imposed the heavy burden on patients and society [44, 45]. A cost-efficient treatment to CP/CPPS should be urgently developed to maximize the therapeutic efficacy, and thus, it could reduce the cost [45]. Given the relationship between therapeutic efficacies and IL-10 and TNF-α phenotypes, Shoskes et al. identified a subgroup of CP/CPPS patients with high IL-10 and low TNF-α genotypes responding poorly to anti-inflammatory phytotherapy [4]. For this reason, patients with different genotypes achieved different therapeutic efficacies. Given the findings of genetic predisposition to CP/CPPS, more studies are aimed at examining the effects of genetic factors on the diagnosis and personalized therapy of CP/CPPS should be conducted [5, 6].

The strengths of this study were that we explored the effects of genetic factors on CP/CPPS, which might provide a novel clue for the diagnosis and treatment of CP/CPPS. Moreover, the functions of the 18 SNPs of IFNG, IFNGR1, and AR were annotated by different databases, and their correlation with the symptoms of CP/CPPS was also analyzed. The nomogram predicting the risk of CP/CPPS was built by integrating age and the identified SNPs, and internal validation showed the adequate goodness-of-fit of the model, which might contribute to the diagnosis of CP/CPPS. The limitations of this study were that only 18 SNPs were genotyped, and other significant SNPs of IFNG (e.g., rs2069705 and rs1861494), IFNGR1 (e.g., rs1327475, rs1327474, and rs2234711), and AR (e.g., CAG, GGC, and GGN repeats and G1733A) were not genotyped. Moreover, the sample size was relatively small for a genetic association study, especially in the subgroup analysis based on participants’ age. The subjects in this study were from a single hospital in China, and external validation was not applicable, thereby restricting the wide use of the results. Moreover, a multicenter study should be conducted to further identify CP/CPPS-related SNPs. In the future, the effects of genetic epidemiology in CP/CPPS should be highlighted. Environmental factors (e.g., the low altitude, higher education level, and more time in the sunlight) were correlated with reduction of prostatitis-like symptoms [46, 47], while cold exposure facilitated the initiation of CP/CPPS [48]. Thus, the application of genetic epidemiology to CP/CPPS might contribute to determining the genetic susceptibility and genetic model of CP/CPPS.

5. Conclusions

Three SNPs were identified (i.e., rs2069718 of IFNG, rs10457655 of IFNGR1, and rs144488434 of AR) to be potentially correlated with CP/CPPS. Genetic polymorphisms of IFNG, IFNGR1, and AR might function as the genetic factors for CP/CPPS susceptibility, and the underlying mechanisms should be further explored.

Data Availability

The data in the current study are available from the corresponding author on reasonable request.

Conflicts of Interest

The authors declare that there is no conflict of interest regarding the publication of this paper.

Authors’ Contributions

Lei Chen and Junyi Chen contributed equally to the work.

Acknowledgments

This work was supported by the National Natural Science Foundation of China (Grant numbers 81970597 and 81630019). Some of the data used in this study were obtained from the Genotype-Tissue Expression (GTEx) Project on July 26, 2021.

Supplementary Materials

See Supplemental Figures 1-3 and Supplemental Tables 1-12 in the Supplementary Materials for comprehensive analysis. (Supplementary Materials)