Endometriosis is an estrogen dependent chronic inflammatory condition characterized by the presence of endometrial tissue outside the uterus that affects women within reproductive ages worldwide (Giudice & Kao, 2004). Multifactorial processes are involved in the etiopathogenesis of this condition resulting in heterogeneous disease manifestation even among the members within the same family, complicating the establishment of a reliable disease classification strategy and therapy (Matalliotakis et al., 2017). Under these circumstances, the prevalence of endometriosis is estimated to be 10% among the symptomatic cases (Laganà et al., 2019). Moreover, due to varying clinical presentations and lack of markers, the average wait for diagnosis remains 7 years from the onset of symptoms, which leads to an increase of 11% of undiagnosed cases each year (Arruda et al., 2003; Laganà et al., 2019). Women with endometriosis are estimated to have twice as much risk for infertility and that nearly 30% diagnosed with the condition develop chronic pelvic pain resistant to conventional treatments (Prescott et al., 2016). Current theories include the influence of genetic, hormonal, immunological, and environmental factors in the pathogenesis of endometriosis, but the exact causes remain elusive (Lalami et al., 2021; Locci et al., 2013). While genetic studies involving twins and familial cases estimated endometriosis heritability as ~50% (Angioni et al., 2020; Deiana et al., 2019), many sufficiently-powered population-based genome-wide association studies (GWAS) and several whole-genome linkage studies attributed only ~26% risk of endometriosis to common genetic variants (Lee et al., 2013). These highlighted the missing heritability of this complex disorder, implicating the necessity to identify rare genetic variants that are not within the scope of GWAS analyses (Albertsen et al., 2019; Lalami et al., 2021).
In GWAS, thousands to millions of single nucleotide polymorphisms (SNPs) are utilized in capturing the genetic variations between cases and controls to reveal risk-conferring alleles. Apart from limitations to known variants in the genome, GWAS analysis mostly require association of non-coding SNPs through assumptions of their impacts on transcriptional regulation of nearby genes (Deiana et al., 2019). On the other hand, ~85% of mutations related to disease traits have been asserted to reside on protein-coding regions of the genome, emphasizing the power of deleterious rare variants in disease pathogenesis (Choi et al., 2009). In whole exome sequencing (WES), coding regions of the genome are captured and subjected to massive parallel sequencing to reveal many gene sequences that can be filtered for rare pathogenic variants. This approach is commonly utilized in genetic diagnosis and research of both monogenic and multifactorial disorders (Rego & Snyder, 2019; Tetreault et al., 2015; Vinkšel et al., 2021; Xuan et al., 2013). Though different approaches to dissecting genetic basis of diseases, findings from GWAS and familial studies can be complementary and thus are both necessary in unraveling endometriosis pathogenesis.
According to the literature, only limited number of studies focused on identifying rare variants related to endometriosis. In this regard, two familial studies utilized WES to reveal novel hemizygous deletions in two genes in a three-generation family of seven affected women and a novel missense genetic variant in three affected sisters in a two-generation Chinese family (Albertsen et al., 2019; Zhu et al., 2022). Another strategy has been to sequence a linkage region obtained from 32 families with endometriosis that disclosed low frequency coding variants in NPSR1 (Tapmeier et al., 2021). Moreover, a large-scale exome-array genotyping study involving European endometriosis cases and controls highlighted the necessity to involve different populations and high-risk families to provide insights into the molecular pathogenesis of the disease (Sapkota et al., 2017).
Therefore, we recruited a high-risk family for endometriosis that is composed of four affected women in a two-generation family. We performed WES in the three of the affected members to identify rare pathogenic coding variants with an aim to contribute to elucidating genetic infrastructure of endometriosis. Our results, together with additional studies could translate into discovering a biomarker that could serve as a non-invasive diagnostic option for the disease.
MATERIALS AND METHODS Ethical complianceThis study was approved by Istanbul Medical Faculty Clinical Research Ethics Committee.
Patients and clinical assessmentsA family having four affected and two unaffected women were recruited for this study (Figure 1a). Patient recruitment took place at Endometriosis Istanbul Clinic within the study period of 2 years. All cases were evaluated and examined by the same endometriosis specialist. Four affected women received transvaginal ultrasonographic (TVS) examination along with a thorough questioning of their medical history, endometriosis related symptomatology, and surgical history, while two unaffected women were interviewed about their clinical history. Venous blood samples were obtained from all recruits, after obtaining written informed consents in accordance with Istanbul Medical Faculty Clinical Research Ethics Commitee (Protocol no: 2021/178 and Decision no: 81985/04). ‘World Endometriosis Research Foundation Endometriosis Phenome and Biobanking Harmonization Project: III. Fluid biospecimen collection, processing, and storage in endometriosis’ standard operating procedures were followed (Rahmioglu et al., 2014). DNA was extracted from peripheral blood utilizing the PureLink Genomic DNA Mini Kit (Thermo Fisher Scientific, Inc., Waltham, MA, USA).
FIGURE 1. Results on whole exome sequencing and segregation analysis of the family. (a) Pedigree of the family is depicted. Star sign indicates members with available genetic material; plus sign marks members subjected to whole exome sequencing (WES). (b) Integrative Genomics Viewer (IGV) images of variants detected by WES are shown. (c) Allele frequencies of pertinent variants are summarized. (d) Chromatograms of the variants in familial segregation analysis are depicted.
WES was performed on the three affected individuals (case id: 301, 302, and 303) at the service of Izmir Tinaztepe University, School of Medicine, Medical Genetic Diagnostic Center (Izmir, Turkey). Twist Comprehensive Exome kit (Twist Bioscience South San Francisco, CA, USA) was utilized in library preparation that targeted 36.8 Mb of protein coding regions covering >99% of RefSeq, CCDS and GENCODE databases. Libraries were subjected to sequencing on the Illumina NextSeq 550 system to achieve a minimum of 20X reading depth for the targeted bases.
Variant annotations and subsequent filtering were achieved utilizing the SEQ platform version 16.7 (
A parallel approach to variant prioritization was adapted to detect common variants among: (i) all exome sequenced members of the family (case id: 301, 302, 303) and (ii) two sisters (case id: 301 and 302), whose mother (case id: 201) had cancer history. In each approach, initial step was to filter for candidate genes for endometriosis retrieved from GWAS and case–control studies along with two familial studies available in the literature (Table S1). In addition, all variants were filtered with a threshold of MAF ≤ 0.01 to detect rare variants, in which provisions of American College of Medical Genetics (ACMG) guidelines were utilized for pathogenicity classification and interpretation of the variants. In silico prediction tools were utilized to assess novel candidate variants. Shared candidate variants were then subjected to confirmation by Sanger sequencing followed by familial segregation. Sequences were analyzed by Benchling Platform (https://
At the time of our clinical evaluation, cases 301 and 303 were already diagnosed with endometriosis through laparoscopic excision of endometrioma and pathological confirmation of disease at different centers in Istanbul and they both sought consultation due to infertility (Figure 1a). During their TVS examination unilateral endometriomas were observed, indicating a recurrence of the disease. On the other hand, cases 302 and 304 both received their diagnosis for the first time by visualization of unilateral endometriomas in both patients with TVS in our clinic. All of these affected individuals presented with endometriosis related pain symptoms (Table 1). Two unaffected individuals (cases 201 and 202), the mothers of these endometriosis patients were already in their post-menopausal years at the time of recruitment (Figure 1a). Case 202 has no known gynecological disorders, whereas case 201 has received an extensive gynecological procedure involving hysterectomy with bilateral oophorectomy due to endometrial serous adenocarcinoma (ESC). Patients' symptomatology and gynecological history are outlined in Table 1.
TABLE 1 Case descriptive: diagnosis, symptomatology, gynecological history.
Cases | Age (years) | Age at diagnosis (years) | Gravidity/Parity | Infertility | Symptoms | Diagnosis | Endometriosis operation history | Other gynecological conditions | AMH (ng/μL) |
201 | 61 | N/A | 2/2 | No | None | None | No | ESC | N/A |
202 | 51 | N/A | 3/3 | No | None | None | No | None | N/A |
301 | 36 | 32 | 0/0 | Yes |
Dysmenorrhea Dyschezia |
LS | Yes | Leiomyoma | 0.26 |
302 | 31 | 31 | 0/0 | Yes |
Dysmenorrhea Chronic pelvic pain |
TVS | No | None | 0.48 |
303 | 30 | 28 | 0/0 | Yes | Dyspareunia | LS | Yes | None | 0.85 |
304 | 27 | 27 | 0/0 | No | Dysmenorrhea | TVS | No | None | 3.91 |
Abbreviations: AMH, anti-Müllerian Hormone; LS, laparoscopic surgery; TVS, transvaginal ultrasonographic; N/A, not available.
Whole exome sequencing resultsThe quality metrics of WES data and filtering for shared variants are summarized in Table 2. We adapted a parallel filtering approach between two groups composed of all exome sequenced members versus cases 301 and 302, whose mother had ESC history. When MAF ≤ 0.01 was set as a threshold to detect rare variants, 4 and 68 variants across the groups involving three and two cases were detected, respectively. These variants were further evaluated for their known/predicted functions, pathogenicity classifications, and novelties in our population that aided in selecting three rare variants in each group as candidates.
TABLE 2 WES data quality and parallel filtering of WES data.
aPathogenicity determined in accordance with ACMG guidelines.
Sanger sequencing confirmations of the variants followed by segregation analysis in the family resulted in two rare heterozygous novel variants as TNFRSF1B (NM_001066.3: c.1072G>A, p.(Ala358Thr)) and GEN1 (NM_001130009.3: c.1574C>T, p.(Ser525Leu)) in the evaluation of three patients, while a heterozygous novel variant in CRABP1 (NM_004378.3:c.54G>C, p.(Glu18Asp)) was detected only in ESC mother descendant two sisters (Figure 1). All of these variants were not only absent from all available databases as outlined in Figure 1c, but were also not found in SEQ-specific cohorts, provided by the SEQ platform that we utilized as a population control set. These variants were classified as variant of unknown significance (VUS) in accordance with ACMG guidelines and that at least one in silico prediction tool supported the pathogenicity of each of these variants.
DISCUSSIONThe results demonstrated that co-segregation of two novel missense variants (TNFRSF1B p.(Ala358Thr) and GEN1 p.(Ser525Leu)) could have played role in the endometriosis pathogenesis in this family. Moreover, the parallel analysis adapted with respect to ESC history revealed a novel CRABP1 p.(Glu18Asp) only in the descendants of case 201.
The anamnesis of cases 201 and 202 indicated that they did not have fertility problems, since they both had spontaneous pregnancies. Neither of them reported having endometriosis related symptoms during their reproductive ages nor were they diagnosed with endometriosis. Case 201 provided her pathology report that confirmed ESC, so further gynecological evaluations related to endometriosis was not possible. Moreover, since both cases are in their post-menopausal ages, which puts their bodies in a low estrogen state, a gynecological investigation at this stage would have been inconclusive. This is due to the fact that endometriosis is an estrogen-dependent disease and in the majority of cases the disease regresses after menopause (Zondervan et al., 2020). However, both cases are found to carry endometriosis related novel genetic variants that are present in this family. This could be explained by two different scenarios; either the cases were asymptomatic or since they did not exhibit infertility, their symptoms such as pelvic pain was overlooked or confused with other disorders of bowel or bladder. Moreover, it is well acknowledged that patient suffering is not indicative of disease severity, so that if not suspected, these cases would not be evaluated in terms of endometriosis (Maddern et al., 2020). Therefore, since cases 201 and 202 were not subjected to gynecological examination, it is difficult to draw a definite conclusion about their disease status.
Among the novel causative genetic variants of endometriosis, TNFRSF1B (*191191) has been implicated to play a role in the pathophysiology of endometriosis through the biological activities of the cytokine TNF-α. This cytokine has a recognized role in endometriosis pathophysiology through TNF-TNFR system that is involved in inflammation, angiogenesis, programmed cell death and proliferation (Ghezzi & Cerami, 2004). Signaling through this system is achieved by binding of TNF-α to its receptors TNFR1 (TNFRSF1A) or TNFR2 (TNFRSF1B), in which TNFR1 favors apoptosis, while TNFR2 can induce apoptosis or promote survival through proliferation, tissue repair or angiogenesis, depending on the microenvironment (Haider & Knöfler, 2009; Islimye et al., 2011). Due to the facts that the components of this TNF-TNFR system were found to be expressed in the endometrium and the endometrial epithelial cells throughout the whole menstrual cycle (Chegini et al., 1999; Iwabe et al., 2000; Tabibzadeh et al., 1995), and that TNFR2 was defined to act in angiogenesis in endothelial cells (Zhang et al., 2003), this system has been evaluated for plausible roles in endometriosis pathology. Among studies performed on human subjects, the expressions of the soluble forms of TNFR1 and TNFR2 were found to be upregulated in both the peritoneal fluid and the serum of women with endometriosis in contrast to their healthy counterparts (Chae et al., 2008; Koga et al., 2000; Othman et al., 2016; Salmeri et al., 2015). On the other hand, the mRNA and protein expressions of TNFR2 were found to be decreased significantly in the endometrial gland and stromal cells of endometriosis patients in the early stages (I and II) (Kharfi et al., 2003). The authors of these studies have hypothesized that these expression changes during endometriosis progress may disrupt TNF-α mediated apoptosis signaling, which in turn may fail to eliminate implanted endometrial cells and/or to facilitate proliferation of ectopic endometrial tissues, leading to endometriosis onset (Kharfi et al., 2003; Koga et al., 2000; Rojas-Cartagena et al., 2005). In addition, an experimental intestinal endometriosis rat model was developed, which also revealed increased expression of sTNFR1 and sTNFR2 in the peritoneal fluids of affected rats (Rojas-Cartagena et al., 2005). Only a single study has associated the TNFR2 haplotype allele defined by T(676)G, A(1663)G and C(1690)T polymorphisms with endometriosis (Chae et al., 2008), while another study utilized endometrial biopsies to report hypomethylation of TNFRSF1B with no significant changes in gene expression in endometriosis patients (Naqvi et al., 2014). Collectively, these studies investigated the role of TNF-TNFR system mainly through genetic expressions, which are under the influence of genetic variations. Therefore, to the best of our knowledge, the novel missense TNFRSF1B p.(Ala358Thr) variant is the first obtained from whole exome sequencing associated with endometriosis in a familial setting.
On the other hand, GEN1 (*612449), exhibiting the other novel causative variant in this family, has not been associated with endometriosis, yet. It is a member of the Rad2/xeroderma pigmentosum group G nuclease family with a defined role in Holliday junction (HJ) resolution, which is required for chromosome segregation during meiosis and for repairing stalled/collapsed replication forks during mitosis (Rass et al., 2010). Therefore, its function has been evaluated in tumorigenesis with a few studies that focused on breast cancer (Kuligina et al., 2013; Turnbull et al., 2010; Wood et al., 2007). In addition, GEN1 has been depicted as a potential candidate gene for the onset of congenital anomalies of the kidney and urinary tract (CAKUT), and Mullerian anomalies (MA) through mouse models (Wang et al., 2018; Wang, Wang, et al., 2020; Wang, Zhang, et al., 2020). It has been depicted that deletion of Gen1 gene needed to be accompanied by a mutation in Wnt9b to exhibit significantly longer uteri and disordered looser stromal structure that is thought to be responsible for the abnormal development of mouse uterus (Wang, Zhang, et al., 2020). This digenic inheritance model is similar to our findings in which GEN1 variant might act synergistically with TNFRS1B to induce the complex endometriosis phenotype in the affected family members. Moreover, a study utilizing mice depicted that Gen1 mutation might induce CAKUT phenotype through disrupted retinoic acid (RA) signaling, which can partly be rescued by all-trans-retinoic acid (ATRA) administration (Zhang et al., 2019). In addition, disrupted RA signaling in Mullerian duct has been associated both with female reproductive tract abnormalities and with reduced proliferative response to sex hormones, which was reinforced by the absence of Mullerian ducts detected in double knockout mice of RA signaling receptors (Nakajima et al., 2019). In terms of endometriosis, RA signaling was found to be involved in resistance to apoptosis due to its key regulatory role in cell survival and that RA has a therapeutic potential with success shown in suppressing the growth of endometriotic lesions and in inhibiting peritoneal cytokine secretion through human and mouse studies (Pavone et al., 2011; Wieser et al., 2012; Yamagata et al., 2015). In light of these findings, we hypothesize that the novel GEN1 p.(Ser525Leu) variant might play a role in endometriosis pathogenesis, since mutations in this gene has been implicated in urogenital disorders with similar phenotypes observed in deficient RA signaling. However, whether GEN1 variants act upon RA signaling mechanism or in combination with accompanying gene variations to induce endometriosis need to be evaluated through further functional studies.
Among the family members in this study, only case 201 had ESC history, which prompted us to perform a parallel WES analysis that aimed to filter for shared genetic variants in cases 301 and 302 against variants in case 303, which would mark inheritance from case 201. This resulted in the identification of a novel CRABP1 p.(Glu18Asp) variant present only in cases 301, 302 and their mother 201. Cellular retinoic acid–binding protein 1 (CRABP1; *180230) is a tumor suppressor gene with an important role in retinoic acid-mediated differentiation and proliferation, which has been associated with varying cancer types (Banz et al., 2010; Nhieu et al., 2022). It is implicated as a highly conserved gene across mammals, thus variations are hypothesized as potential inducers of tumor formation or progression (Nhieu et al., 2022). However at this stage, it is difficult to draw a definite conclusion about the role of CRABP1 p.(Glu18Asp) variant in ESC onset, since there is only case 201 with a clinical diagnosis. Therefore, being the first-degree relatives, cases 301 and 302 have increased risk for ESC, so these patients are under a long-term follow-up. Since these cases exhibit endometriosis and the literature provides conflicting reports on the association between endometriosis and endometrial cancer (Johnatty et al., 2020; Kvaskoff et al., 2020), their possible diagnosis with ESC will be evaluated for the contribution of CRABP1 p.(Glu18Asp) variant during their regular follow-up.
Our study involves limitations due to availability of only four members in one generation of the family for clinical diagnosis, where only clinical history could be obtained from the mothers (case ids: 201 and 202). Moreover, since all the patients could not be operated due to fertility wish, endometriosis staging could not be provided and that disease tissue could not be obtained for further analysis. However, family members are under clinical follow-up and if surgery is suggested, their surgery materials will be evaluated for the effect of detected novel variants. Familial studies face these limitations, but their power in detecting rare pathogenic variants due to smaller gene pool is of value in identification of enigmatic conditions such as endometriosis.
CONCLUSIONOur findings suggest TNFRSF1B (NM_001066.3: c.1072G>A, p.(Ala358Thr)) and GEN1 (NM_001130009.3: c.1574C>T, p.(Ser525Leu)) as two novel variants that could be responsible for endometriosis diagnosed in this family. Thereby our results not only provided novel genetic markers for endometriosis that merits to be replicated in expanded cohorts or validated by functional analysis but also provided the family with genetic counseling.
AUTHOR CONTRIBUTIONSAmong authors, Busra Gizem Kina took responsibility in WES analysis and confirmation of candidate variants, supervised by Feyza Nur Tuncer. Sevcan Aydin contributed to WES analysis by providing candidate genes list for endometriosis through extensive literature review. Engin Oral provided clinical evaluations of the patients, whereas Nura Fitnat Topbas Selcuki and Pinar Yalcin Bahat helped with clinical follow-up and clinical data collection. Taner Usta helped in writing clinical parts of the manuscript. Busra Gizem Kina, Feyza Nur Tuncer, and Nura Fitnat Topbas wrote the manuscript, whereas Nilufer Rahmioglu provided critical review for the article. Funding was provided through Feyza Nur Tuncer's project. The work was directed equally by Engin Oral and Feyza Nur Tuncer as supervisors for clinical and genetics parts of the study, respectively.
ACKNOWLEDGMENTSThe authors thank the patients for taking part in our study. This work was supported by the Scientific Research Projects Coordination Unit of Istanbul University. Project Number: TYL-2021-37712.
CONFLICT OF INTEREST STATEMENTThe authors declare no conflict of interest.
DISCLOSURE STATEMENTPart of the findings of this work was presented as a poster at European Society of Human Reproduction and Embryology (ESHRE) Meeting during 3–6 July 2022. This work was conducted as a Master of Science Thesis of Busra Gizem Kina at the Genetics Program of Istanbul University Graduate School of Health Sciences, under the supervision of Feyza Nur Tuncer, PhD.
DATA AVAILABILITY STATEMENTThe data that support the findings of this study are available on request from the corresponding author. The data are not publicly available due to privacy or ethical restrictions.
ETHICS STATEMENTThe study was approved by Istanbul Medical Faculty Clinical Research Ethics Committee (Protocol no: 2021/178 and Decision no: 81985/04).
HUMAN RIGHTS STATEMENTS AND INFORMED CONSENTAll methods were performed in accordance with the Declaration of Helsinki and relevant guidelines and regulations. All participants signed written informed consent forms.
You have requested "on-the-fly" machine translation of selected content from our databases. This functionality is provided solely for your convenience and is in no way intended to replace human translation. Show full disclaimer
Neither ProQuest nor its licensors make any representations or warranties with respect to the translations. The translations are automatically generated "AS IS" and "AS AVAILABLE" and are not retained in our systems. PROQUEST AND ITS LICENSORS SPECIFICALLY DISCLAIM ANY AND ALL EXPRESS OR IMPLIED WARRANTIES, INCLUDING WITHOUT LIMITATION, ANY WARRANTIES FOR AVAILABILITY, ACCURACY, TIMELINESS, COMPLETENESS, NON-INFRINGMENT, MERCHANTABILITY OR FITNESS FOR A PARTICULAR PURPOSE. Your use of the translations is subject to all use restrictions contained in your Electronic Products License Agreement and by using the translation functionality you agree to forgo any and all claims against ProQuest or its licensors for your use of the translation functionality and any output derived there from. Hide full disclaimer
© 2024. This work is published under http://creativecommons.org/licenses/by-nc/4.0/ (the “License”). Notwithstanding the ProQuest Terms and Conditions, you may use this content in accordance with the terms of the License.
Abstract
Background
Endometriosis is an estrogen-dependent, chronic inflammatory disease that affects 10% of women during the reproductive ages. Despite the estimated 50% heritability for the condition, only 26% was associated with common genetic variants. Thus, necessity of identifying rare variants for the missing heritability is implicated in the literature. Therefore, our study aimed to identify novel rare genetic variants involved in the pathogenesis of endometriosis utilizing a family of multiple affected members.
Methods
A family composed of four affected women along with their two unaffected mothers were recruited at a single gynecology and infertility clinic specialized in endometriosis. All patients presented with endometriomas, which was visualized by transvaginal ultrasonography. Two affected individuals had received laparoscopic endometrioma excision and therefore were diagnosed with recurrent disease. One mother had a history of endometrial serous adenocarcinoma (ESC) for which she underwent hysterectomy with bilateral oophorectomy. Three endometriosis cases were whole exome sequenced on Illumina NextSeq 550 platform with an average of 90% coverage. Candidate genes were confirmed by Sanger sequencing and followed-up with family segregation.
Results
Novel rare variants were identified in
Conclusion
Novel candidate genetic variants that might contribute to endometriosis were suggested that need replication through independent cohorts or validation by functional studies. The family has also received genetic counseling and that the affected daughters are on clinical follow-up, accordingly.
You have requested "on-the-fly" machine translation of selected content from our databases. This functionality is provided solely for your convenience and is in no way intended to replace human translation. Show full disclaimer
Neither ProQuest nor its licensors make any representations or warranties with respect to the translations. The translations are automatically generated "AS IS" and "AS AVAILABLE" and are not retained in our systems. PROQUEST AND ITS LICENSORS SPECIFICALLY DISCLAIM ANY AND ALL EXPRESS OR IMPLIED WARRANTIES, INCLUDING WITHOUT LIMITATION, ANY WARRANTIES FOR AVAILABILITY, ACCURACY, TIMELINESS, COMPLETENESS, NON-INFRINGMENT, MERCHANTABILITY OR FITNESS FOR A PARTICULAR PURPOSE. Your use of the translations is subject to all use restrictions contained in your Electronic Products License Agreement and by using the translation functionality you agree to forgo any and all claims against ProQuest or its licensors for your use of the translation functionality and any output derived there from. Hide full disclaimer
Details


1 Department of Genetics, Aziz Sancar Institute of Experimental Medicine, Istanbul University, Istanbul, Turkey; Graduate School of Health Sciences, Istanbul University, Istanbul, Turkey
2 Department of Obstetrics and Gynecology, Istanbul Sisli Hamidiye Etfal Training and Research Hospital, University of Health Sciences Turkiye, Istanbul, Turkey
3 Department of Obstetrics and Gynecology, Istanbul Kanuni Sultan Suleyman Training and Research Hospital, University of Health Sciences Turkiye, Istanbul, Turkey
4 Department of Obstetrics and Gynecology, Acibadem Altunizade Hospital, Mehmet Ali Aydinlar University, Istanbul, Turkey
5 Oxford Endometriosis Care Centre, Nuffield Department of Women's and Reproductive Health, University of Oxford, Women's Centre, John Radcliffe Hospital, Oxford, UK; Wellcome Centre for Human Genetics, University of Oxford, Oxford, UK
6 Department of Genetics, Aziz Sancar Institute of Experimental Medicine, Istanbul University, Istanbul, Turkey
7 Department of Obstetrics and Gynecology, Bezmialem Vakif University, Istanbul, Turkey