Skip to main content

Applicability of polygenic risk scores in endometriosis clinical presentation

Abstract

Background

Risk prediction is an essential part of preventative medicine and in recent years genomic information has become an interesting factor in risk models. Polygenic risk scores (PRS) combine the effect of many genetic variations into a single score which has been shown to have predictive value for many diseases. This study aimed to investigate the association between PRS for endometriosis and the clinical presentation of the disease.

Methods

Women with endometriosis (N = 172) were identified at the Department of Gynecology. All participants answered questionnaires regarding sociodemographic factors, lifestyle habits and medical history, registered bowel symptoms on the Visual Analog Scale for Irritable Bowel Syndrome and passed blood samples. DNA was extracted and samples were genotyped, and a PRS was calculated based on previous genome-wide association studies of endometriosis. Inflammatory proteins and TSH receptor antibodies (TRAb) in serum were analyzed.

Results

Inverse associations were identified between PRS and spread of endometriosis, involvement of the gastrointestinal tract and hormone treatment. However, significance was lost when calculated as p for trend and the specificity and sensitivity were low. There were no correlations between PRS and TRAb or inflammatory proteins.

Conclusion

The findings indicate that specific PRS should be developed to predict clinical presentations in patient with endometriosis.

Peer Review reports

Background

Endometriosis is a chronic, inflammatory gynecological disease, affecting up to 10% of women of reproductive age [1, 2]. The disease is defined by ectopic endometrial cells and stroma growing outside of the uterus and is associated with dysmenorrhea, pelvic pain, gastrointestinal symptoms and infertility. The most established theory of pathogenesis is Sampson´s theory of retrograde menstruation [3]. However, retrograde menstruation has been seen in up to 90% of women [4], which indicates that other mechanisms, such as genetic, immunological and environmental factors are involved in the development of endometriosis [5]. New theories also suggest that bone marrow-derived stem cells might migrate through the peripheral circulation, causing endometriosis in sites outside of the peritoneal cavity [6].

The hereditary component of endometriosis, explained by genetic variations, has in previous studies been estimated to 47–51% [7, 8]. There are different types of genetic variations contributing to disease. One common form of allelic variation is Single Nucleotide Polymorphism (SNP) and in total, SNPs explain approximately 26% of the polygenic risk of endometriosis [9]. Genome-wide association studies (GWAS) are commonly used to discover genetic variants associated with a disease [10, 11]. Several GWAS of endometriosis have identified genetic variants related to development and regulation of the female reproductive tract [12]. Individually, these genetic variants have low impact on disease risk. Calculating a polygenic risk score (PRS) is a way to assess the risk of a condition based on the combined effect of many genetic variants. Using a PRS has in many diseases been shown to be better for risk prediction than individual SNPs [13, 14]. However, to our knowledge, PRS has never been applied for endometriosis.

Endometriosis is widely known as an inflammatory condition and serum concentrations of the inflammatory proteins Axis inhibition protein 1 (AXIN1), Sulfotransferase Family 1A Member 1 (ST1A1) and C-X-C Motif Chemokine Ligand 9 (CXCL9) have been found to be altered in women with endometriosis [15]. The inflammatory proteins Oncostatin-M (OSM), Monocyte Chemoattractant protein 1 (MCP-1) and TNF Receptor Superfamily Member 9 (TNFRSF9) have been shown to be associated with different symptoms of endometriosis [15]. Several studies have highlighted the association between thyroid dysfunction and endometriosis [16,17,18]. Recently, elevated titers of thyroid-stimulating hormone (TSH) receptor antibodies (TRAb) have been shown to be associated with endometriosis, suggesting a connection between the two pathologies [19].

The primary aim of the present study was to investigate the association between PRS for endometriosis and clinical presentations of the disease. Secondary aims were to investigate associations of PRS for different inflammatory proteins and TRAb.

Materials and methods

Ethical consent for the study was obtained from the Ethics Review Board of Lund University, No 2012/564, 2016/56.

Study participants and design

Women with endometriosis were identified at the Department of Gynecology at Skåne University Hospital, Malmö, according to the ICD-10 classification of endometriosis, N80. The recruitment took place between March 2013 and March 2017 according to previously described inclusion and exclusion criterium [17]. The main inclusion criterium was to have a definite diagnosis of endometriosis confirmed by laparoscopy or laparotomy. General inclusion criteria were to comprehend the Swedish or English language and an age above 18 years. Exclusion criteria were an uncertain diagnosis of endometriosis, multiple or severe somatic or psychiatric comorbidities, a diagnosis of inflammatory bowel disease (IBD), current pregnancy or living too far from the hospital. In total, 172 women were finally included in the study.

All study participants answered a questionnaire and completed the Visual Analog Scale for Irritable Bowel Syndrome (VAS-IBS). They also passed blood samples, and whole blood and sera were kept frozen at − 80 °C. Inflammatory proteins and TRAb in serum were analyzed. DNA was extracted and samples were genotyped, and a PRS was calculated based on previous GWAS of endometriosis.

Study questionnaire

A previously developed questionnaire regarding socioeconomic factors, medical history and pharmacological treatments was completed by all study participants. In the questionnaire they also answered questions about their endometriosis-associated symptoms such as onset of the disease, trigger factors, examinations and treatments (Additional file 1) [19].

The visual analog scale for irritable bowel syndrome

A validated, unlicensed questionnaire called VAS-IBS was used to estimate bowel symptoms (Additional file 2) [20]. The patients used the VAS-IBS questionnaire to estimate the severity of their symptoms over the past 2 weeks regarding abdominal pain; constipation; diarrhea; bloating and flatulence; vomiting and nausea; and intestinal symptom’s influence on daily life, as previously described [17]. They estimated their symptoms on a continuous scale from 0 to 100 mm where 0 represents no symptoms and 100 represents severe symptoms. The scales were inverted from the original format [20]. Reference values are based on controls, recruited among hospital staff at Skåne University Hospital, Malmö, consisting of 52 women who had not undergone prior abdominal surgery, with a median age of 44 (range 22–77) years [21].

Biological sample analyses

A multiplex proximity extension assay, Proseek®Multiplex Inflammation 1 kit (Olink Bioscience, Uppsala, Sweden), was used to analyze 92 inflammatory proteins in serum, as previously described [15]. Measured proteins are expressed in an arbitrary unit on a log2 scale celled Normalized Protein Expression (NPX). Larger NPX values represent higher protein expressions in the sample. The proteins ST1A1, CXL9, OSM, MCP-1 and TNFRSF9 from this analysis, who had previously been shown to be associated with endometriosis and gastrointestinal symptoms, were used in this study [13]. Sandwich ELISA technology (MBS762601, MyBiosource, San Diego, CA, USA, lot.no.H2497C119) was used to analyze AXIN1 in plasma, as previously described [13]. Analysis of TRAb IgG autoantibodies was conducted in serum according to clinical routines, using a competitive Electro Chemi Luminescence Immunoassay (ECLI) detection technique based on Ruthenium deriva [22].

Sample sequencing and quality control

DNA samples of 172 patients were genotyped using Illumina Global Screening Array, on Illumina iScan high-throughput screening system at the Institute of Clinical Molecular Biology (Christian-Albrechts-University, Kiel, Germany). The GenCall algorithm implemented in Illumina GenomeStudio software was used to get the alleles from raw intensity data.

Genotyped data was quality controlled (QC) removing samples and markers using the following pipeline: exclusion of samples with ≥ 15% missing rates; exclusion of markers with non-called alleles; exclusion of markers with missing call rates > 0.05; exclusion of samples with ≥ 5% missing rates; exclusion of related samples (PI-HAT > 0.1875); exclusion of samples whose genotyped sex could not be determined; exclusion of samples with high heterozygosity rate (more than three times standard deviation (SD) from mean); only autosomal SNPs were kept; removal of markers with Hardy–Weinberg equilibrium P-value < 1 × 10–5; removal of markers whose P-value of difference in missingness between cases and control was < 1 × 10–5; and removal of samples which were outliers, identified using principal component analysis (deviation of more than 6 times interquartile range). At the end, 140 samples passed the QC.

Imputation of missing genotypes was done using TOPMed Imputation Server [23]. The reference panel used was TOPMed Version R2 on GRC38; and default pipeline was used. Once imputed, markers with INFO score < 0.80, MAF < 0.01 and non-biallelic markers were removed.

Polygenic risk score

For PRS calculation, the results of a genome-wide association study about endometriosis [10] deposited in GWAS catalog with the accession number GCST004549 [24] was used as model. From the available SNPs, the 13 SNPs with p-value < 5 × 10–8 and present in our data were kept for the application of PRS.

Then, PRS was calculated as it is implemented in PLINK software [34] (version 1.9), using unweighted (presence of the risk allele) and weighted (using the beta value of the effect of the risk allele) methods.

Principal components

Four principal components (PC) were calculated for each patient and used as covariates to control for population stratification. The genotyped data was pruned to get SNPs with no linkage disequilibrium using PLINK software [25] and SNPs from high LD regions were excluded. Then, FlashPCA was used to calculate the principal components of SNP data [26].

Data categorization

Data was categorized as previously described [27].

Statistical analyses

Statistical analyses were performed using the software SPSS© statistical computer package version 26 for Windows. Values are presented as median and interquartile range or numbers and percentages. Since the quantitative data was not normally distributed, comparisons of basic characteristics between groups were performed by Mann–Whitney U test and PRS was divided into quartiles or the 10% with highest PRS compared with the rest. Fisher´s exact test was used for dichotomous variables and for comparison of disease distribution among quartiles. Logistic regression was used to calculate odds ratio (OR) with 95% confidence interval (CI) and p for trends with localization, symptom or hormonal treatment as dependent variable and weighted or unweighted PRS as independent variable, adjusted for PC1–PC4, using both lowest and highest quartiles for PRS as reference values. Since age differed between patients with isolated and spread endometriosis and between patients with and without hormone therapy, age was added as a covariate in these calculations. Receiver operating characteristic (ROC) curves were used to show area under the curve (AUC) with 95% CI for PRS. Generalized linear model was used to calculate β value and 95% CI with inflammatory proteins and TRAb as dependent variable and weighted or unweighted PRS and PC1-PC4 together with confounders as independent variables. p < 0.05 was considered statistically significant.

To investigate if combinations of proteins were associated with the PRS, a combination of cross-validation and LASSO-regression was employed. We divided the data into a training and test set using a randomizing split. The shrinkage parameter (λ) was estimated using k-fold cross-validation in the training set. To perform variable selection the estimated shrinkage parameter λCV was used in the test set. The selected variables and the absolute value of the coefficients were saved, and the process was repeated 10 times. Further, variables were ordered by the number of times they were selected and the sum of its estimated coefficients. The lowest ranked variable was removed, and the whole process was repeated until a final model was selected. This method has previously been described in Leandersson et al. (2020) [28]. The final models were estimated with linear regression and compared against an intercept only model with likelihood ratio (LR) tests. Cross-validation and LASSO-regression analyses were carried out using R v 4.0.0 [29].

Results

Subject characteristic

In total, 172 patients were included in the study (Table 1). Patients were divided into subgroups depending on isolated or spread endometriosis, with or without gastrointestinal involvement, gastrointestinal symptoms and hormone therapy (Additional file 3: Table S1–S4).

Table 1 Sociodemographic factors, lifestyle habits and gastrointestinal symptoms in endometriosis patients

Basic characteristics of patients with endometriosis. BMI = Body Mass Index. Values are presented as median (interquartile range) or numbers and percentages (%). Gastrointestinal symptoms were measured on the Visual Analog Scale for Irritable Bowel Syndrome, where 0 mm represents no symptoms and 100 mm represents very severe symptoms [20]. Reference values from healthy controls are shown [21].

Patients with isolated ovarian endometriosis had a higher age compared to patients with spread endometriosis (40.0 (34.0–45.0) years vs. 38.0 (31.0–42.0) years, p = 0.034) and fewer were physically active ≥ 1 h/week (37.3% vs. 57.0%, p = 0.012) (Additional file 3: Table S1). No differences in basal characteristics were observed between patinets with and without endometriosis involving the gastrointestinal tract (Additional file 3: Table S2). Patients without gastrointestinal symptoms had a higher age, (40.0 (34.0–46.5) years versus 38.0 (32.0–43.0) years, p = 0.036) and a lower prevalence of hormone therapy compared to patients with gastrointestinal symptoms (24.0% vs. 50.3%, p = 0.017) (Additional file 3: Table S3). Patients without hormone therapy had a higher age compared to patients with current hormone treatment (39.0 (33.0–44.0) years vs. 37.0 (31.0–42.0) years, p = 0.046) and were to a greater extent married or cohabitants than single or living alone, compared with patients with current hormone therapy (p = 0.002). There were some differences in gastrointestinal symptoms between patients with and without hormone therapy, with more abdominal pain, diarrhea and vomiting and nausea in those with hormone therapy (Additional file 3: Table S4). Smoking habits, alcohol habits, education, occupation, and opioid treatment were similar in the different subgroups (Additional file 3: Table 1–4).

Genetical analyses

In the third quartile of weighted and unweighted PRS, a lesser number of patients had spread endometriosis compared to the lowest quartile (OR: 0.252; 95% CI 0.081–0.782, p = 0.017 and OR: 0.182; 95% CI 0.052–0.0630, p = 0.007, respectively) as well as compared to the highest quartile (OR: 0.409; 95% CI 0.136–1.288, p = 0.111 and OR: 0.245; 95% CI 0.077–0.781, p = 0.017). There was no significant association between the PRS and localization when calculated as p for trend (Table 2 and 3).

Table 2 Association of weighted PRS in different subgroups of endometriosis patients
Table 3 Association of unweighted PRS in different subgroups of endometriosis patients

When calculating the associations between the 10% with highest PRS and localization of endometrios, no associations could be found (OR: 1.853; 95% CI 0.441–7.786; p = 0.400 for unweighted and OR: 1.886; 95% CI 0.428–8.311; p = 0.402 for weighted). The sensitivity and specificity of PRS to predict isolated or spread endometriosis was low (Fig. 1 a–d).

Fig. 1
figure 1

ROC-curves showing the sensitivity and specificity of PRS to predict isolated/spread endometriosis (a), GI tract not involved /GI tract involved (b), no GI symptoms/GI symptoms (c) and no hormone therapy /hormone therapy (d)

The second quartile of weighted PRS was associated with a lower risk of involvement of the gastrointestinal tract (OR: 0.158; 95% CI 0.026–0.949, p = 0.044), and the second quartile of unweighted PRS tended to be associated with a lower risk of gastrointestinal involvement (Table 23). The third quartile of unweighted PRS conferred with a lower risk of having hormone therapy (OR: 0.250; 95% CI 0.075–0.829, p = 0.023) (Table 3). There were no associations between PRS and TRAb or any of the inflammatory proteins investigated in this study (Table 4).

Table 4 Association of weighted and unweighted PRS, and inflammatory proteins and TSH receptor antibodies

Calculations of combinations of proteins did not either show any associations between PRS and TRAb or inflammatory proteins (Figs. 2 and 3).

Fig. 2
figure 2

Residual plot from linear regression model for weighted PRS, comparing a model including MPC1, AXIN1, hormone therapy, CXCL9, TRAb titer, and OSM to an intercept only model. Red and blue horizontal lines are 0 ± 1.96*sd(residual)

Fig. 3
figure 3

Residual plot from linear regression model for unweighted PRS, comparing a model including MPC1, AXIN1, hormone therapy, CXCL9, and TRAb titer to an intercept only model. Red and blue horizontal lines are 0 ± 1.96*sd(residual)

Discussion

The main result in the present study was the inverse associations between the third quartile of weighted and unweighted PRS and spread endometriosis, between the second quartile of weighted PRS and gastrointestinal involvement, and the third quartile of unweighted PRS and hormone therapy. However, there seems to be no clinical use of PRS for endometriosis development to predict the clinical presentation in this entity since sensitivity and specificity for clinical outcome was low. There was no association between PRS and circulating inflammatory proteins or TRAb.

The clinical use of PRS has not yet been defined. For several diseases, PRS seems to be a promising clinical tool to calculate risks and treatment strategies, e.g., coronary artery diseases [30], dementia [31] and generalized epilepsy [32]. On the other hand, PRS is not enough to estimate the risk of coronary artery disease since non-genetic risk factors, such as age and environmental factors, are also involved [33]. When studying the association between PRS of low-density lipoprotein cholesterol (LDL-C) with circulating LDL-C levels and coronary artery disease, there were no associations within a cohort of patients with familial hypercholesterolemia [34]. Accordingly, PRSs are predictive of colorectal cancer in the general population but not in the Lynch syndrome [35]. Further, PRS was not able to have a great impact on discrimination models in some diseases, e.g., breast cancer [36].

Patients with endometriosis have an impaired psychological well-being compared with controls and are often diagnosed with depression and prescribed antidepressants [17]. Psychiatric comorbidity may exaggerate the experience of pain and gastrointestinal symptoms [37], however the causality in patients with endometriosis is unclear. We do not know whether endometriosis causes depression or whether depressed subjects have more severe symptoms and therefore more often seek health care and are admitted to a tertiary center. In a great community-based cohort, PRS for major depressive disorder, bipolar disease and schizophrenia were associated with all types of mental illness [38]. The comorbidity between depression and endometriosis may be a confounder interacting with SNPs found in GWAS and raise the question whether the associations found are representative for endometriosis or depression.

Our results show that genetic variants involved in developing endometriosis are not enough to explain the clinical presentation of the disease. For example, bowel symptoms, which are common in endometriosis patients, could be caused by other possible conditions such as microbiome dysbiosis [39], small bowel permeability [40] or perineal dysfunction [41]. Bowel symptoms could also be related to opioid treatment which is common in this group of patients [42] or caused by adherent tissue formation following abdominal surgery. Of the patients in our study 165 (95.9%) had undergone prior abdominal surgery. These factors affect the clinical presentation of endometriosis but are not possible to genetically identify. Endometriosis is known to be a chronic inflammatory disease and several studies have reported an increased level of inflammatory mediators in serum and peritoneal fluid [43, 44]. Further, the number of macrophages was higher in ovarian endometriomas tissue than in ovarian cysts [45]. SNPs associated with endometriosis in GWAS are encoding proteins involved in DNA-replication, transcriptional process, embryonic development, cytoskeleton regulation, angiogenic pathways, hormonal pathways and inflammatory pathways [10]. The interleukin 1A gene locus, involved in the regulation of proinflammatory cytokines and chemokine production, has been identified to have a significant association with endometriosis [43]. The associations found in GWAS of endometriosis thus implicate that SNPs involved in inflammatory regulation are associated with endometriosis. Taken together, the associations found in GWAS may represent depression and/or inflammatory conditions and are thus not representative for the clinical or laboratory outcome of endometriosis. These confounders may obscure the genetic picture.

GWAS of endometrios have identified genetic variants to differentiate endometriosis from healthy controls. The available genetic information is thus not intended to interrogate the genetic basis of the clinical presentation. Genetic variants involved in the risk of developing endometriosis seems not to be enough to explain the clinical or laboratory presentations of the disease.

Thus, it is needed to analyze the genetic variants involved in the disease presentation in a large cohort to develop useful PRS for this task. There are several limitations in this study. One of the most important limitations of this study is the small sample size of the cohort, recruited from a tertiary center. Another limitation is the limited characteristics without staging of the endometriosis cohort and no hormonal analyses.

Conclusions

The genetic variants involved in the risk of developing endometriosis are not enough to explain the clinical presentations. Thus, the analysis of genetic variants involved in the clinical manifestations is needed to develop useful PRS for the clinical presentations.

Availability of data and materials

Data available on request due to ethical restrictions.

Abbreviations

AUC:

Area under the curve

AXIN1:

Axis inhibition protein 1

BMI:

Body mass index

CXCL9:

C-X-C motif chemokine ligand 9

ECLI:

Electro chemi luminescence immunoassay

GI:

Gastrointestinal

GnRH:

Gonadotropin-releasing hormone

GWAS:

Genome-wide association study

IBS:

Inflammatory bowel disease

LR:

Likelihood ratio

MCP-1:

Monocyte chemoattractant protein 1

NPX:

Normalized protein expression

OR:

Odds ratio

OSM:

Oncostatin-M

PC:

Principal components

PEA:

Proximity extension assay

PRS:

Polygenic risk score

ROC:

Receiver operating characteristic

SD:

Standard deviation

SNP:

Single nucleotide polymorphism

ST1A1:

Sulfotransferase family 1A member 1

TNFRSF9:

TNF receptor superfamily member 9

TSH:

Thyroid stimulating hormone

TRAb:

TSH receptor antibodies

QC:

Quality controlled

VAS-IBS:

Visual analog scale for irritable bowel syndrome

References

  1. Giudice LC, Kao LC. Endometriosis. Lancet. 2004;364(9447):1789–99.

    Article  Google Scholar 

  2. Ahn SH, et al. Pathophysiology and immune dysfunction in endometriosis. Biomed Res Int. 2015;2015: 795976.

    PubMed  PubMed Central  Google Scholar 

  3. Sampson JA. Metastatic or Embolic endometriosis, due to the menstrual dissemination of endometrial tissue into the venous circulation. Am J Pathol. 1927;3(2):93–110.

    CAS  PubMed  PubMed Central  Google Scholar 

  4. Halme J, et al. Retrograde menstruation in healthy women and in patients with endometriosis. Obstet Gynecol. 1984;64(2):151–4.

    CAS  PubMed  Google Scholar 

  5. Laganà AS, et al. The pathogenesis of endometriosis: molecular and cell biology insights. Int J Mol Sci. 2019;20(22):5615. https://doi.org/10.3390/ijms20225615.

    CAS  Article  PubMed Central  Google Scholar 

  6. Laganà AS, et al. Stem cell trafficking during endometriosis: may epigenetics play a pivotal role? Reprod Sci. 2018;25(7):978–9.

    Article  Google Scholar 

  7. Treloar SA, et al. Genetic influences on endometriosis in an Australian twin sample. Fertil Steril. 1999;71(4):701–10. https://doi.org/10.1016/S0015-0282(98)00540-8.

    CAS  Article  PubMed  Google Scholar 

  8. Gao X, et al. Economic burden of endometriosis. Fertil Steril. 2006;86(6):1561–72.

    Article  Google Scholar 

  9. Lee SH, et al. Estimation and partitioning of polygenic variation captured by common SNPs for Alzheimer’s disease, multiple sclerosis and endometriosis. Hum Mol Genet. 2013;22(4):832–41.

    CAS  Article  Google Scholar 

  10. Dehghan A. Genome-wide association studies. Methods Mol Biol. 2018;1793:37–49.

    CAS  Article  Google Scholar 

  11. NHGRI. A Catalog of Published Genome-Wide Association Studies.; Available from: http://www.genome.gov/gwastudies

  12. Sapkota Y, et al. Meta-analysis identifies five novel loci associated with endometriosis highlighting key genes involved in hormone metabolism. Nat Commun. 2017;8:15539.

    CAS  Article  Google Scholar 

  13. Dudbridge F. Power and predictive accuracy of polygenic risk scores. PLoS Genet. 2013;9(3): e1003348.

    CAS  Article  Google Scholar 

  14. Lambert SA, Abraham G, Inouye M. Towards clinical utility of polygenic risk scores. Hum Mol Genet. 2019;28(R2):R133-r142.

    CAS  Article  Google Scholar 

  15. Ek M, et al. AXIN1 in plasma or serum is a potential new biomarker for endometriosis. Int J Mol Sci. 2019;20(1):189. https://doi.org/10.3390/ijms20010189.

    CAS  Article  PubMed Central  Google Scholar 

  16. Sinaii N, et al. High rates of autoimmune and endocrine disorders, fibromyalgia, chronic fatigue syndrome and atopic diseases among women with endometriosis: a survey analysis. Hum Reprod. 2002;17(10):2715–24.

    CAS  Article  Google Scholar 

  17. Poppe K, Velkeniers B. Female infertility and the thyroid. Best Pract Res Clin Endocrinol Metab. 2004;18(2):153–65.

    CAS  Article  Google Scholar 

  18. Yuk JS, et al. Graves disease is associated with endometriosis: a 3-year population-based cross-sectional study. Medicine. 2016;95(10): e2975.

    Article  Google Scholar 

  19. Ek M, et al. Characteristics of endometriosis: a case-cohort study showing elevated IgG titers against the TSH receptor (TRAb) and mental comorbidity. Eur J Obstet Gynecol Reprod Biol. 2018;231:8–14.

    CAS  Article  Google Scholar 

  20. Bengtsson M, Ohlsson B, Ulander K. Development and psychometric testing of the visual analogue scale for irritable bowel syndrome (VAS-IBS). BMC Gastroenterol. 2007;7:16.

    Article  Google Scholar 

  21. Bengtsson M, et al. Evaluation of gastrointestinal symptoms in different patient groups using the visual analogue scale for irritable bowel syndrome (VAS-IBS). BMC Gastroenterol. 2011;11:122.

    Article  Google Scholar 

  22. Laboratoriemedicin Skåne, Analysportalen. [cited 2021 June 29]; Available from: http://www.analysportalen-labmedicin.skane.se.

  23. Das S, et al. Next-generation genotype imputation service and methods. Nat Genet. 2016;48(10):1284–7.

    CAS  Article  Google Scholar 

  24. Buniello A, et al. The NHGRI-EBI GWAS Catalog of published genome-wide association studies, targeted arrays and summary statistics 2019. Nucleic Acids Res. 2019;47(D1):D1005-d1012.

    CAS  Article  Google Scholar 

  25. Chang CC, et al. Second-generation PLINK: rising to the challenge of larger and richer datasets. Gigascience. 2015;4:7.

    Article  Google Scholar 

  26. Abraham G, Qiu Y, Inouye M. FlashPCA2: principal component analysis of Biobank-scale genotype datasets. Bioinformatics. 2017;33(17):2776–8.

    CAS  Article  Google Scholar 

  27. Svensson A, et al. Associations between endometriosis and gut microbiota. Reprod Sci. 2021;28(8):2367–77.

    CAS  Article  Google Scholar 

  28. Leandersson P, et al. A multiplex biomarker assay improves the diagnostic performance of HE4 and CA125 in ovarian tumor patients. PLoS ONE. 2020;15(10): e0240418.

    CAS  Article  Google Scholar 

  29. R: A language and environment for statistical computing. 2018; R Core Team. R Foundation for Statistical Computing]. Available from: http://www.R-project.org/.

  30. Christoffersen M, Tybjærg-Hansen A. Polygenic risk scores: how much do they add? Curr Opin Lipidol. 2021;32(3):157–62.

    CAS  Article  Google Scholar 

  31. Riaz M, et al. Effect of APOE and a polygenic risk score on incident dementia and cognitive decline in a healthy older population. Aging Cell. 2021;20(6):e13384.

    CAS  Article  Google Scholar 

  32. Stevelink R, et al. Shared genetic basis between genetic generalized epilepsy and background electroencephalographic oscillations. Epilepsia. 2021;62(70):1518–27.

    Article  Google Scholar 

  33. Menniti G, et al. multiscale risk factors of cardiovascular disease: CLSA analysis of genetic and psychosocial factors. Front Cardiovasc Med. 2021;8: 599671.

    Article  Google Scholar 

  34. Nomura A, et al. Polygenic risk scores for low-density lipoprotein cholesterol and familial hypercholesterolemia. J Human Genet. 2021;66(11):1079–87. https://doi.org/10.1038/s10038-021-00929-7.

    CAS  Article  Google Scholar 

  35. Jenkins MA, et al. Assessment of a polygenic risk score for colorectal cancer to predict risk of lynch syndrome colorectal cancer. JNCI Cancer Spectr. 2021. https://doi.org/10.1093/jncics/pkab022.

    Article  PubMed  PubMed Central  Google Scholar 

  36. Li SX, et al. Prospective evaluation of the addition of polygenic risk scores to breast cancer risk models. JNCI Cancer Spectr. 2021. https://doi.org/10.1093/jncics/pkab021.

    Article  PubMed  PubMed Central  Google Scholar 

  37. Trivedi MH. The link between depression and physical symptoms. Prim Care Companion J Clin Psychiat. 2004;6(Suppl 1):12–6.

    Google Scholar 

  38. Wainberg M, et al. Symptom dimensions of major depression in a large community-based cohort. Psychol Med. 2021. https://doi.org/10.1017/S0033291721001707.

    Article  PubMed  Google Scholar 

  39. D’Alterio MN, et al. Possible role of microbiome in the pathogenesis of endometriosis. Minerva Obstet Gynecol. 2021;73(2):193–214.

    PubMed  Google Scholar 

  40. Viganó D, et al. How is small bowel permeability in endometriosis patients? a case control pilot study. Gynecol Endocrinol. 2020;36(11):1010–4.

    Article  Google Scholar 

  41. Raimondo D, et al. Pelvic floor dysfunction at transperineal ultrasound and chronic constipation in women with endometriosis. Int J Gynaecol Obstet. 2022. https://doi.org/10.1002/ijgo.14088.

    Article  PubMed  Google Scholar 

  42. Ek M, et al. Gastrointestinal symptoms among endometriosis patients: a case-cohort study. BMC Womens Health. 2015;15:59.

    Article  Google Scholar 

  43. Sapkota Y, et al. Association between endometriosis and the interleukin 1A (IL1A) locus. Hum Reprod. 2015;30(1):239–48.

    CAS  Article  Google Scholar 

  44. Zhang T, et al. The link between immunity, autoimmunity and endometriosis: a literature update. Autoimmun Rev. 2018;17(10):945–55.

    CAS  Article  Google Scholar 

  45. Laganà AS, et al. Evaluation of M1 and M2 macrophages in ovarian endometriomas from women affected by endometriosis at different stages of the disease. Gynecol Endocrinol. 2020;36(5):441–4.

    Article  Google Scholar 

Download references

Acknowledgements

We want to acknowledge Per Ekström, Lil Valentin and Johanna Nordengren from the Department of Gynecology, Skåne University Hospital, for identifying and helping with the recruitments of patients. This research was funded by grants from Bengt Ihre Foundation, Dir Albert Påhlsson´s Foundation and Development Foundation of Region Skåne.

Funding

Open access funding provided by Lund University. This research was funded by grants from Bengt Ihre Foundation, Dir Albert Påhlsson´s Foundation and Development Foundation of Region Skåne. The funding sources had no roles in the design of the study, collection, analysis, and interpretation of data or in writing the manuscript.

Author information

Authors and Affiliations

Authors

Contributions

Conceptualization, BO and M.D.; methodology, AS, AÅ, BO, CB, KGE, and MD; software, AS, AÅ, BO, KGE, MD; validation, AS, BO and ME; formal analysis, AS, AÅ, BO, KGE and MD; investigation, BR and ME; resources, BO and MD; data curation, AS, BO and ME; writing—original draft preparation, AS; writing—review and editing, AS, AÅ, BO, BR, CB, KGE, ME and MD; visualization, AS, AÅ and BO; supervision, BO and MD; project administration, BO; funding acquisition, BO, BR and MD. All authors have read and agreed to the published version of the manuscript. All authors read and approved the final manuscript.

Corresponding author

Correspondence to Agnes Svensson.

Ethics declarations

Ethics approval and consent to participate

The study was conducted according to the guidelines of the Declaration of Helsinki and approved by the Ethics Review Board of Lund University No 2012/564 (09/10/2012) and 2016/56 (03/05/2016). Informed, oral and written, consent was obtained from all subjects involved in the study before entrance to the study. The participants gave their written informed consent to use clinical/patient data and publish the results. The Swedish Biobank approved to use blood samples and the Swedish Authority for Privacy Protection approved the genetic analyses. Diarienr 1565–2012.

Consent for publication

Informed, oral and written, consent was obtained from all subjects involved in the study before entrance to the study.

Competing interests

The authors declare no conflict of interest.

Additional information

Publisher's Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary Information

Additional file 1:

 Patient survey.

Additional file 2:

 Visual Analogue Scale for Irritable Bowel Syndrome (VAS-IBS).

Additional file 3:

 Supplementary Tables.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and Permissions

About this article

Verify currency and authenticity via CrossMark

Cite this article

Svensson, A., Garcia-Etxebarria, K., Åkesson, A. et al. Applicability of polygenic risk scores in endometriosis clinical presentation. BMC Women's Health 22, 208 (2022). https://doi.org/10.1186/s12905-022-01788-w

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12905-022-01788-w

Keywords

  • Endometriosis
  • Gastrointestinal symptoms
  • Genetics
  • Inflammatory proteins
  • Polygenic risk score