Abstract

This study is to investigate the relationship of programmed cell death 1 (PD-1; also known as PDCD1) and programmed death-1-ligand 1 (PD-L1; also known as CD274) single nucleotide polymorphisms (SNPs) with polycystic ovary syndrome (PCOS). This study enrolled 330 PCOS patients and 350 matched controls. ELISA was used to detect the PD-1 and PD-L1 levels in serum. SnaPshot genotyping was performed to analyze the PD-1 and PD-L1 genotyping. Linkage disequilibrium and haplotype of TagSNP loci of PD-1 and PD-L1 genes were also detected. The relationship of genotypes and alleles with PCOS was analyzed. The levels of PD-1 and PD-L1 in the serum of PCOS patients were significantly lower than those in the control group (). The haplotype TT of PD-1 gene at rs10204525 and rs7421861 loci was significantly lower in the PCOS group than in the control group (, , and ). PD-L1 gene SNP loci rs2282055, rs2890658, rs10125854, and rs702275 had linkage disequilibrium. The haplotypes TAAA, GAAC, GAGC, GCAA, and TCGA of PD-L1 gene SNP loci were significantly higher in PCOS patients than in the control group, whereas haplotypes GAAA, TAAC, TCAA, GCGA, GCAC, and TCGC of PD-L1 gene SNP loci were significantly lower in PCOS patients than in the control group. PD-1 and PD-L1 SNPs may be related to the pathogenesis of PCOS. PD-1 gene SNP loci rs10204525 and rs7421861 and PD-L1 gene SNP loci rs2282055, rs2890658, rs10125854, and rs702275 may be new candidate polymorphic loci for PCOS.

1. Introduction

Polycystic ovary syndrome (PCOS) is a metabolic disease of endocrine disorders in women. Its incidence in women of childbearing age is 5%-10% [1–3]. PCOS is mainly manifested as abnormal menstruation, infertility, obesity, hairy, polycystic ovarian changes, hyperandrogenemia, insulin resistance, and induced hyperinsulinemia [4]. With the changes in people’s lifestyle and environment, the incidence of PCOS has been increasing year by year [5]. At present, the pathogenesis of PCOS is still unclear. Epidemiological survey reveals obvious family clustering characteristics of PCOS [6]. Chen et al. have reported 11 susceptible genes related to PCOS disease and extensively tested them in European and American populations [7]. Therefore, genetic factors play a very important role in the etiology of PCOS.

Programmed cell death 1 (PD-1; also known as PDCD1) and programmed death-1-ligand 1 (PD-L1; also known as CD274) play important negative regulatory roles in inflammation, autoimmune diseases, and tumors. They are widely involved in suppressing inflammatory response and tumor immune evasion [8, 9]. Genetic mutations in PD-1 and PD-L1 not only increase the risk of kidney cancer [10] and thyroid cancer [11], but also are closely related to the course, clinical stage, and metastasis of the tumor and even affect the function of kidney transplantation [12]. In a multicenter meta-analysis by Mohammad et al. [13], it is found that PD-1 and PD-L1 gene polymorphisms were related to susceptibility and progression of all cancers. It is reported that PD-1 and PD-L1 gene polymorphisms can affect the susceptibility to ankylosing spondylitis [14, 15]. In addition, PD-L1 gene polymorphisms are reported to be associated with lung cancer susceptibility [16, 17]. Therefore, PD-1 and PD-L1 gene polymorphisms are closely related to the occurrence of diseases such as cancer. However, whether the PD-1 and PD-L1 gene polymorphisms are related to the occurrence and development of PCOS has not been reported.

In this study, we explored the correlation between PD-1 and PD-L1 gene polymorphisms and the pathogenesis of PCOS.

2. Materials and Methods

2.1. Patients

This study included 330 patients diagnosed with PCOS at the Reproductive Fertility Center of the First Affiliated Hospital of Xinjiang Medical University. Meanwhile, 350 non-PCOS women of childbearing age who matched the clinical indicators of PCOS patients were recruited as control group. The inclusion criteria for the PCOS patient group included the following: (1) the diagnosis of PCOS met the diagnostic criteria established by the Rotterdam International Conference [18]. Patients with any two items of the following three items can be diagnosed: little or no ovulation, hyperandrogenemia,  mL on ultrasound imaging, and ≥12 follicles with 2-9 mm in diameter. (2) The patients did not take drugs that may affect hormones, blood glucose, or insulin at 3 months before the blood collection. (3) The patients did not implement planned weight loss. Exclusion criteria for patients with PCOS included the following: (1) patients with congenital adrenal hyperplasia, Cushing syndrome, and ovarian or adrenal tumors; (2) patients with hyperprolactinemia, thyroid disease, or pituitary tumor; (3) patients with severe primary diseases such as cardiovascular disease, liver disease and kidney disease, or mental illness; (4) patients with pregnancy or suspicious pregnancy, lactation, suffering from other endocrine diseases, with long-term oral contraceptives, or subcutaneous implanted contraception. The inclusion criteria of the control group included the following: (1) patients with infertility due to tubal obstruction or male reasons. (2) Their age, body mass index, and other clinical indicators were matched with those in the PCOS group. (3) They had regular menstruation and normal serum androgens. (4) They had no clinical hyperinsulinemia, no symptoms of hyperandrogenemia, and no clinical manifestations of PCOS such as hairy, acne, and skin seborrhea. Written informed consent was obtained from every patient, and the study was approved by the Ethics Committee of the First Affiliated Hospital of Xinjiang Medical University.

2.2. Detection of Clinical Indicators

We recorded the clinical information of all patients and normal controls, including age, weight, height, waist circumference, and hip circumference. Body mass index was calculated as . Waist-to-hip ratio (WHR) was calculated as . Peripheral blood (3-5 mL) was collected from patients and controls. The serum levels of follicle-stimulating hormone (FSH), luteinizing hormone (LH), and testosterone (T) were measured using a chemiluminescence immunoassay analyzer (i2000SR, Abbott). Fasting blood glucose (FBG) was measured using hexokinase method (AU5821, Beckman, Germany). Fasting insulin (FNS) was detected by chemiluminescence method (i2000SR, Abbott). The formula for calculating the insulin resistance index (HOMA-IR) is as follows: .

2.3. ELISA

The serum levels of PD-1 and PD-L1 were detected with corresponding ELISA kits (E-EL-H1534c/E-EL-H1547c Elabscience Biotechnology Co. Ltd, Wuhan, China) according to the kit instructions. OD450 was measured, and the concentrations of PD-1 and PD-L1 in the specimens were calculated by the standard curve.

2.4. DNA Extraction and Genotyping

DNA was extracted from peripheral blood using TIANamp Blood DNA Kit (TIANGEN, Beijing, China). Data on the linkage disequilibrium between PD-1 and PD-L1 gene single nucleotide polymorphisms (SNPs) in the Chinese population were obtained by querying the NCBI database. By setting the linkage disequilibrium coefficient of and the % as the standard, we used the Haploview software to select three TagSNP loci representing the four SNP loci of the entire PD-1 gene and 12 TagSNP loci representing 24 SNP loci for the entire PD-L1 gene. Finally, the HapMap map was generated (Figure 1).

2.5. SnaPshot Typing Test

By using the SnaPshot typing technique [19], we detected a total of 15 SNP loci in PD-1 and PD-L1 genes. The sequences of PD-1 and PD-L1 were retrieved through Gene Bank. The primers were designed using the Primer 5.0 software and synthesized by Shanghai Sangon Biotech (China). The primer sequences are shown in Table 1. Briefly, after PCR amplification of the genomic DNA sample, the amplified PCR product (2 μL) was purified and mixed with 1 μL SnaPshot Mix and 0.2 μL primers, and then, the reaction system was supplemented with ddH2O to 6.0 μL to perform PCR reaction. Then, 1 μL of the product was taken, and 9 μL of the loading Hi-Di was added. The samples were denatured at 95°C for 3 min, and an ice-water bath was performed immediately. After that, the sample was added with Hi-Di Formamide 9 μL, denatured at 95°C for 5 min, and then sequenced with an ABI3730XL sequencer. GeneMapper4.1 (AppliedBiosystems Co. Ltd, USA) was used to analyze the data.

2.6. Statistical Analysis

SPSS17.0 was used to analyze the data. The data was expressed as . The Hardy-Weinberg (H-WE) equilibrium was used to detect the genotype frequency distribution of TagSNP loci of the PD-1 gene and PD-L1 gene in the PCOS group and the control group, respectively. SHEsis (http://analysis.bio-x.cn/myAnalysis.php) was used to analyze the linkage disequilibrium and haplotype of the PD-1 gene SNP loci rs7421861 and rs10204525 and PD-L1 gene SNP loci rs2282055, rs2890658, rs702275, and rs10125854. The independent t test was used to compare the data between the two groups. Chi-square test was used to compare the genotype and allele frequency, and the correlation was analyzed by logistic regression [20]. The OR value and 95% confidence interval were calculated. was considered statistically significant.

3. Results

3.1. Clinical Information of PCOS Group and Control Group

A total of 330 patients with PCOS and 350 controls were included in this study. They were all from Xinjiang Uygur Autonomous Region, China. They were aged 20-39 years. The clinical data of included subjects are shown in Table 2. There was no statistically significant difference in age, BMI, WHR, and FSH between PCOS group and control group (). Statistical results showed that the levels of T, LH, LH/FSH, FBG, FNS, and HOMA-IR in the PCOS group were significantly higher than those in the control group (Table 2, ). These results indicate that the PCOS group has significant clinical differences compared with the normal control group.

3.2. Serum PD-1 and PD-L1 Level in Patients

ELISA was conducted to measure the serum levels of PD-1 and PD-L1. The results showed that the serum levels of PD-1 (Figure 2(a)) and PD-L1 (Figure 2(b)) in patients with PCOS were  ng/mL and  ng/mL, which were significantly lower than those in the control group ( ng/mL and  ng/mL, respectively) (), suggesting that abnormalities of PD-1 and PD-L1 in the serum of PCOS patients may be related to the onset of PCOS.

3.3. PD-1 and PD-L1 Gene SNPs’ H-WE Test

The SnaPshot typing results were shown in supplementary file 1. H-WE test was used to detect the gene frequency distribution of PD-1 gene and PD-L1 gene TagSNP loci in PCOS group and control group, respectively. The results showed that the genotype distribution of the TagSNP loci of the PD-1 gene and PD-L1 gene in the PCOS group and the control group was consistent with the H-WE with good agreement (Table 3, ). Thus, these genotypes can be used for subsequent haplotype analysis.

3.4. SNPs of PD-1 and PD-L1 in Patients

Through HapMap analysis, we obtained 15 TagSNP loci of the PD-1 gene and PD-L1 gene. Among them, the allele frequency of the rs7421861 and rs10204525 in the PD-1 gene SNP loci and the rs2282055, rs2890658, rs702275, and rs10125854 in the PD-L1 SNP loci gene had statistically significant differences between the PCOS group and the control group (Table 4, ). These results indicate that the SNPs of PD-1 gene and PD-L1 gene might be related to the pathogenesis of PCOS disease.

In addition, we found that the 15 TagSNP loci of PD-1 gene and PD-L1 gene in the PCOS group and the control group all conformed to H-WE. At the same time, by using the SHEsis software, we performed linkage disequilibrium and haplotype analyses for the six SNP sites with statistical differences in Table 4, including rs7421861 and rs10204525 of the PD-1 gene SNP loci and the rs2282055, rs2890658, rs702275, and, rs10125854 of the PD-L1 gene SNP loci. The results showed that the haplotype TT of rs7421861 and rs10204525 of the PD-1 SNP loci had a haplotype frequency of 0.45 in the control group and of 0.36 in the PCOS group. The TT haplotype distribution frequency in the PCOS group was significantly lower than that in the control group (Table 5; , , and ). Moreover, in the PCOS group and the control group, there were significant differences in the haplotypes GAAA, TAAA, GAAC, TAAC, GAGC, GCAA, TCAA, GCGA, TCGA, GCAC, and TCGC of the rs2282055, rs2890658, rs702275, and rs10125854 SNP loci of the PD-L1 gene (Table 6, ). The above results suggest that the haplotypes of the SNP loci of PD-1 and PD-L1 genes are related to the occurrence and development of PCOS disease.

4. Discussion

Chronic inflammation has been a uniform clinical feature in patients with PCOS. The PD-1, a member of the B7/CD28 costimulatory molecule family, is a cosuppression signal molecule that inhibits T cell activation [21]. PD-L1, PD-1 ligand, can induce T cell apoptosis and stimulate cells to secrete IL-10, thereby mediating immunosuppressive effects [22]. The physiological function of PD-L1 is to regulate the inflammatory response and limit tissue damage. Thus, the interaction between PD-1 and PD-L1 can prevent excessive inflammation and suppress the immune response.

In this study, we tested the levels of PD-1 and PD-L1 in serum. The results showed that compared with the control group, the levels of PD-1 and PD-L1 were lower in the serum of PCOS patients, which is consistent with the study of Ni et al. [23]. The low levels of PD-1 and PD-L1 in the serum of PCOS patients may be related to their inability to effectively inhibit T and B cell function and proliferation, thus leading to the pathogenesis of PCOS. However, whether the decrease in PD-1 and PD-L1 levels can be used as an immunological diagnosis and treatment index for PCOS disease requires further research.

Mounting evidence suggests that PCOS might be a complex multigenic disorder strongly influenced by epigenetic and environmental factors, including diet and lifestyle factors [24]. Pedigree studies have shown that PCOS has significant family clustering characteristics [25, 26]. At present, the largest twin studies in PCOS have found that the incidence of PCOS was twice as high in monozygotic twins as in dizygotic twins, indicating that genetic factors play a very important role in PCOS pathogenesis [27]. In this study, we used SnaPshot typing to detect PD-1 and PD-L1 gene polymorphisms in PCOS and analyzed the correlation between PD-1 and PD-L1 gene polymorphisms and PCOS. The results showed that the allele frequencies of the PD-1 gene SNP loci (rs7421861 and rs10204525) and PD-L1 gene SNP loci (rs2282055, rs2890658, rs702275, and rs10125854) were significantly different between the two groups, indicating that these SNP loci may be related to the development of PCOS.

In addition, in order to investigate whether the frequency distribution of genes and genotypes of PCOS patients and control groups has reached genetic balance and to understand whether they are representative of the population, we use the H-WE to detect the genotype frequency distribution of PD-1 gene SNP loci rs7421861 and rs10204525 and PD-L1 gene SNP loci rs2282055, rs2890658, rs702275, and rs10125854. H-WE test can reflect the representative of the population [28, 29]. The results confirmed that these SNP loci conformed to H-WE. These results showed that these SNP loci were not only significantly different between the PCOS group and the control group, but also were representative of the population.

In this study, we analyzed the linkage disequilibrium and haplotype between SNP loci of the PD-1 gene. The results showed that there was a linkage disequilibrium between rs10204525 and rs7421861 of the PD-1 gene SNP loci. At the same time, the haplotype TT of PD-1 gene at rs10204525 and rs7421861 loci was less frequently distributed in the PCOS group than in the control group. It is reported that PD-1 gene polymorphism is related to some inflammatory diseases and autoimmune diseases, such as allergic asthma [30]. The PD-1 gene polymorphism rs2227981 is associated with tumor occurrence, including digestive system tumors and female-specific tumors [31]. PD-1 polymorphism is associated with susceptibility to diseases such as gastric cancer, clear cell renal cancer, HBV, and HCV [32–34]. After the PD-1 gene is knocked out, its original negative regulatory effect would be lost, resulting in increased levels of inflammatory factors, which promotes the imbalance of proinflammatory and anti-inflammatory cytokines and thereby promotes the occurrence of related diseases [35]. PD-1 is also related to immune dysfunction [36]. Moreover, inflammation is one of the clinical features of PCOS women [37]. Thus, we believe that the immune response involving PD-1 gene polymorphism is related to the pathogenicity of PCOS. When the T allele and T allele of PD-1 gene are mutated at the same time, haplotype TT is a protective factor for the occurrence of PCOS disease. The haplotype of PD-1 gene is related to the occurrence and development of PCOS disease, which indicates that the PD-1 gene polymorphism is related to the pathogenesis of PCOS.

In the linkage disequilibrium and haplotype analysis of the PD-L1 gene SNP loci rs2282055, rs2890658, rs10125854, and rs702275, the results showed that there was a linkage disequilibrium between the PD-L1 gene SNP loci rs2282055, rs2890658, rs10125854, and rs702275. And the haplotype TAAA, GAAC, GAGC, GCAA, and TCGA frequencies in the PCOS group were significantly higher than the control group. However, the haplotype GAAA, TAAC, TCAA, GCGA, GCAC, and TCGC frequencies were significantly lower in the PCOS group than in the control group. Studies have reported that the PD-L1 gene rs4143815 may be associated with an increased risk of gastric cancer [38, 39] and non-small-cell lung cancer [40, 41], which may be caused by increased PD-L1 expression and suppression of immune tumor surveillance [42]. Women with PCOS have an increased risk of endometrial cancer, ovarian cancer, endocrine adenocarcinoma, pancreatic cancer, and other cancers [43]. We believe that the polymorphism of the PD-L1 gene affects the expression of PD-L1, which in turn promotes the onset of PCOS. These results suggest that haplotypes TAAA, GAAC, GAGC, GCAA, and TCGA might be susceptibility factors for the occurrence of PCOS, while haplotypes GAAA, TAAC, TCAA, GCGA, GCAC, and TCGC might be protective factors for the occurrence of PCOS. The haplotype of PD-L1 gene SNP loci is related to the occurrence and development of PCOS which indicates that the polymorphism of PD-L1 gene is related to the pathogenesis of PCOS.

However, there are still some limitations in the study. First, the number of cases in this study is limited. Second, the specific pathways and mechanisms of PD-1 gene and PD-L1 gene in the pathogenesis of PCOS were not investigated. Further studies with large sample size are still needed to reveal the underlying mechanisms.

5. Conclusions

In summary, in this study, we reported for the first time that PD-1 gene and PD-L1 gene may be involved in the pathogenesis of PCOS. Our results showed that the genotype distribution, allele frequency, and haplotype of the PD-1 gene SNP loci rs10204525 and rs7421861 and PD-L1 gene rs2282055, rs2890658, rs10125854, and rs702275 were all related to the pathogenesis of PCOS. In the follow-up studies, we will further explore the correlation of PD-1 and PD-L1 gene polymorphisms with the clinical characteristics of PCOS. These results will help screen and protect susceptible populations at the genetic level, will also provide theoretical basis for individualized diagnosis and treatment programs for PCOS patients, and will provide a basis for developing novel biological markers for early screening and diagnosis of PCOS.

Data Availability

The data that support the findings of this study are available on request from the corresponding author.

Conflicts of Interest

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

Authors’ Contributions

Rui Han and Xiaoyun Gong contributed equally to this work.

Acknowledgments

This work was supported by the State Key Laboratory of Pathogenesis, Prevention and Treatment of High Incidence Diseases in Central Asia Fund (grant number SKL-HIDCA-2020-JZ10) and the Regional Science Fund Project of National Natural Science Foundation of China (grant number 81960289).

Supplementary Materials

Supplementary file 1. The SnaPshot typing results. (Supplementary Materials)