|Home | About | Journals | Submit | Contact Us | Français|
Previous microarray analyses identified 22 microRNAs (miRNAs) differentially expressed in paired ectopic and eutopic endometrium of women with and without endometriosis. To investigate further the role of these miRNAs in women with endometriosis, we conducted an association study aiming to explore the relationship between endometriosis risk and single-nucleotide polymorphisms (SNPs) in miRNA target sites for these differentially expressed miRNAs. A panel of 102 SNPs in the predicted miRNA binding sites were evaluated for an endometriosis association study and an ingenuity pathway analysis was performed. Fourteen rare variants were identified in this study. We found SNP rs14647 in the Wolf–Hirschhorn syndrome candidate gene1 (WHSC1) 3′UTR (untranslated region) was associated with endometriosis-related infertility presenting an odds ratio of 12.2 (95% confidence interval = 2.4–60.7, P = 9.03 × 10−5). SNP haplotype AGG in the solute carrier family 22, member 23 (SLC22A23) 3′UTR was associated with endometriosis-related infertility and more severe disease. With the individual genotyping data, ingenuity pathways analysis identified the tumour necrosis factor and cyclin-dependant kinase inhibitor as major factors in the molecular pathways. Significant associations between WHSC1 alleles and endometriosis-related infertility and SLC22A23 haplotypes and the disease severe stage were identified. These findings may help focus future research on subphenotypes of this disease. Replication studies in independent large sample sets to confirm and characterize the involvement of the gene variation in the pathogenesis of endometriosis are needed.
Endometriosis is a common disease causing pelvic pain and reduced fertility in 6–10% of reproductive aged women at considerable cost to health systems (Giudice and Kao, 2004; Berkley et al., 2005; Guo and Wang, 2006). Although the aetiology of endometriosis is still not well understood, there is wide support for Sampson's theory that endometrial fragments are carried in retrograde menstrual flow through the fallopian tubes and implant in the pelvic cavity (Sampson, 1927). Other aetiological factors have been proposed to influence establishment and progression of the disease, including hormonal (Nyholt et al., 2009), immune (Weiss et al., 2009) and environmental alterations (Guo, 2009). Genetic inheritance has been shown to be a major risk factor for endometriosis in studies that show the incidence of endometriosis is increased among the relatives of endometriosis patients compared with population controls (Zondervan et al., 2001).
Molecular studies have focused considerable effort on expression, regulation and inheritance of genes in endometriotic disease. In recent years, it has been shown that small non-coding microRNAs (miRNAs), are important components of complex gene regulatory networks (Reinhart and Bartel, 2002; Aravin et al., 2007). miRNAs are endogenous ~19–22 nt RNAs that are a component of the miRNA-induced silencing effector complex (mRISC). The microRNA sequence directs the sequence-specific binding of the mRISC complex to one or more target mRNAs, repressing their translation to protein (Bartel, 2009). In this way, protein production can be inhibited and cellular functions altered without altering the target gene's sequence. However, a genetic aberration in miRNA binding could result in an inherited defect in protein translation of a target mRNA. We hypothesized that an inherited single-nucleotide polymorphism (SNP) and/or a mutation near a miRNA binding site may have potential impact in women with endometriosis and contribute to the development of endometeriotic disease. Recent research (Halvorsen et al., 2010) has shown that SNPs can affect the RNA structure of untranslated regions (UTRs), and these changes in conformation are known to affect the accessibility of miRNA target sites (Doench and Sharp, 2004). Furthermore, research (Wang et al., 2006) swapping miRNA target sites between UTRs-demonstrated miRNA target sites are repressionally active due to contextual signals surrounding target sites. Another study highlighted the importance of the sequence surrounding miRNA target sites for normal miRNA regulation by mutating the regions surrounding the target sites and showing loss of miRNA repression (Mencia et al., 2009).
Recent evidence has shown that miRNAs are involved in many biological processes, including cell differentiation, proliferation and apoptosis that commonly occur in endometriosis (Pan and Chegini, 2008; Glinskii et al., 2009). Compelling studies have implicated miRNAs in various human diseases, such as cancers (Kontorovich et al., 2009; Satzger et al., 2009), diabetes (Yang and Kaye, 2009), schizophrenia (Coyle, 2009), inflammation and viral infection (Hufner et al., 2009; Mattes et al., 2009), as well as in endometriosis (Pan et al., 2007; Burney et al., 2009; Ohlsson Teague et al., 2009b). Ye et al. (2008) found an association between oesophageal cancer risk and a functional SNP located in the pre-mir423 binding site. In addition, Chin et al. and Tian et al. provided evidence for association between lung cancer and the number of miRNA functional SNPs (Chin et al., 2008; Tian et al., 2009).
Previous microarray analyses identified 22 miRNAs differentially expressed in paired ectopic and eutopic endometrium of women with and without endometriosis (Ohlsson Teague et al., 2009b). To investigate the potential impact of these miRNAs in women with endometriosis, we conducted an association study aiming to explore the relationship between SNPs in target sites for these differentially expressed miRNA and endometriosis risk. We designed and evaluated a panel of 102 SNPs in predicted miRNA target sites and conducted an association study in endometriosis cases and controls. Given the key regulatory function of miRNAs in gene expression stability, understanding the mechanisms of how endometrial miRNAs are regulated and identifying their specific target genes might lead to the development of novel therapeutic treatments for endometriosis.
Participants were 958 endometriosis cases and 959 controls from an Australian study of endometriosis. The vast majority of cases (n=926) were drawn from families with affected sister pairs (ASPs) (one per family), and 32 were from non-ASP families. The cases with surgically confirmed endometriosis and with the most severe stage of disease were chosen from ASPs (Zhao et al., 2006). Disease severity was assessed using the revised American Fertility Society (rAFS) classification system (The American Fertility Society, 1985). However, the time from diagnosis to recruitment was variable and there was variable detail in the clinical records available. We therefore collapsed stage of disease into two classes for analysis (rAFS I/II and rAFS III/IV) (Zhao et al., 2007). Cases with rAFS Stages I and II were more common (59%) than the cases with rAFS Stages III and IV (41%). Twenty-five percent of cases (n=236) self-reported a history of ‘problems conceiving’, and were therefore considered to be a subfertility group of women diagnosed from endometriosis. The remainder reported no problems conceiving (44.2%) or had never tried to conceive (29.3%).
A similar number of 959 unrelated female controls were selected from a twin study of gynaecological health (Treloar et al., 1999). These controls had never been diagnosed with the disease on the basis of either self-report or medical records were available. No evidence of endometriosis was reported in the 27% of control women who reported having a hysterectomy and/or laparoscopy, therefore supporting a low risk of endometriosis in our control sample. The mean ages (±SD) of the cases and controls at the time of data collection were 35.82±8.87 and 45.60±11.98 years, respectively. Ethics approval was obtained from the Human Research Ethics Committee of the Queensland Institute of Medical Research and the Australian Twin Registry.
Twenty-two miRNAs were previously identified by microarray analysis (Ohlsson Teague et al., 2009a). Potential miRNA target sites in genomic sequences for the list of miRNAs were assessed by using the Miranda algorithm (Betel et al., 2008). The binding energy of hybridization between the miRNA and target as ≤−30 kcal/mole and a cut-off score of 120 were considered to be efficient for functional binding specificity. The genome locations of these target sites were determined and searched against dbSNP (build 126) to identify genes within target sites. In this case if a gene was considered a target based on the criteria above, all target sites for that miRNA within that gene were searched for SNPs regardless of the binding energies of the individual sites. When searching for SNPs within the target sites, 2 nt were added to each side of all the target sites since SNPs which are close to the target site may also affect miRNA binding. In total, there were 2657 target sites across 145 genes included in the search of dbSNP database. A total of 243 SNPs were identified within target sites using the UC Santa Cruz genome browser (UCSC version 18, 2006).
We selected a subset of these SNPs for genotyping by first selecting any target SNPs located in regions of significant or suggestive linkage to endometriosis on chromosomes 10, 7 and 20 (Treloar et al., 2005; Zondervan et al., 2007). A total of 24 SNPs were identified including 6 SNPs on chromosome 10, 9 SNPs on chromosome 7 and 9 SNPs on chromosome 20. To supplement this list, we then screened the SNPs using ingenuity pathways analysis (IPA) software (http://www.ingenuity.com/) and the UCSC genome browser (http://genome.ucsc.edu/cgi-bin/hgGateway). Of the remaining 219 target SNPs, 78 SNPs were also selected into the study excluding 123 unmapped SNPs or the SNPs located in unknown genes. Eighteen SNPs failed in the assay design stage, as determined by Sequenom Assay Design software (Sequenom, version 3.1).
SNP sequences were downloaded from the National Center for Biotechnology Information (http://www.ncbi.nlm.nih.gov/) and were cross checked in Sequenom databases (https://mysequenom.com) before assay design. Multiplex assays were designed for these 102 miRNA target SNPs using the Sequenom MassARRAY Assay Design software (version 3.1). SNPs were typed using iPLEXTM Gold chemistry and analysed using a Sequenom MassARRAY Compact Mass Spectrometer (Sequenom Inc., San Diego, CA, USA). The PCR and the iplex reactions were performed according to the manufacturer's instructions as described in our previous study (Zhao et al., 2006). Briefly, the 2.5 ml PCRs were performed using 12.5 ng genomic DNA, 1 unit of Taq polymerase (Sequenom Inc.), 500 µM of each dNTP, 2 mM of MgCl2 and 100 nM of each PCR primer (Integrated DNA Technologies Inc., Coralville, IA, USA). PCR thermal cycling program and a two-step 200 short cycles program for extension reaction were used as described in our previous study (Zhao et al., 2006). The post-PCR products were spotted on a SpectroChip (Sequenom Inc.), and data were processed and analysed by MassARRAY TYPER 4.0 software (Sequenom Inc.).
SNP genotypes were tested for departures from Hardy–Weinberg equilibrium (HWE). The PLINK genetic analysis package (http://pngu.mgh.harvard.edu/purcell/plink/) was used to calculate allele frequencies and to access the association of miRNA target SNPs with endometriosis. To obtain nominal P-values, 10 000 permutation tests were performed to obtain a region-wide empirical P-value for each SNP and SNP haplotype. To identify specific relationships between disease phenotype and SNP genotype, we conducted analyses stratifying cases by disease severity and phenotypic characteristics. Each subgroup was compared with the large Australian control sample. We assessed statistical significance at P<0.05 after correction for testing data for multiple SNPs, which for the endometriosis subgroup analyses became 0.01 (0.05/3) after permutation testing. Pairwise linkage disequilibrium (LD), haplotype frequencies and blocks were determined by Haploview version 4.1 using the default method of Gabriel et al. (2002). To test for non-independence of association signals, we analysed multiple SNPs jointly by logistic regression using the PLINK program (Arya et al., 2009).
We searched for SNPs in target sites for published miRNAs differentially expressed in ectopic versus eutopic endometrial tissues and genotyped 102 miRNA target SNPs in an Australian sample comprising 958 endometriosis cases and 959 controls. The miRNA binding sites for each SNP are shown in the Supplementary data, table. Genotypes were obtained for most individuals with a mean completion rate of 99.4%. Of the 102 markers, 47 were not polymorphic (monomorphic) in our case–control sample. There were 41 polymorphic variants with allele minor allele frequencies >0.1%. Genotype frequencies for the 41 polymorphisms were in HWE in the control sample. The minor allele frequencies of the 41 miRNA target SNPs ranged from 0.16 to 46.4% in the cases and 0.05 to 46.8% in the controls (Table I). There was some evidence for allelic association between endometriosis and SNPs rs35091219 and rs1736215 with an asymptotic pointwise P<0.01 (Table I). The minor T allele of SNP rs35091219 was associated with endometriosis (case frequency=0.017; control frequency=0.007), whereas the minor C allele of SNP rs1736215 was more frequent in controls (case frequency=0.421; control frequency=0.469). However, the differences were not significant after correcting for multiple testing (P>0.1).
Fourteen variants were rare in our sample with the variant allele seen in no more than three individuals and these were not included in the association analysis. Of these, six were found only in endometriosis cases, five were found only in unrelated controls and three were found in both cases and controls (Table II).
We further evaluated association results by stratifying cases according to disease severity and phenotypic characteristics. Differences were observed with nominal P-values<0.05 for several SNPs (rs3211066, rs7894262, rs1736215, rs874232 and rs3803947) and higher minor allele frequencies were observed for SNP rs35091219 in all different case groups compared with controls. These differences were not significant after correcting for multiple testing (data not shown). The smallest pointwise P-value was 9×10−5 [χ2=15.3, odd ratio (95% confidence interval)=12.2 (2.4–60.8)] for SNP rs14647, located in the Wolf–Hirschhorn syndrome candidate 1 gene (WHSC1) 3′UTR in 236 endometriosis subfertility cases compared with 959 unrelated controls (Table III). This difference remained significant after correction for multiple testing (P=0.005). It was noted, the minor A allele frequencies for this SNP were very low, at 0.013 in the 236 subfertility cases and 0.001 in the 959 controls (Table III).
To identify a relationship between miR-99b and SNP rs14647 in the WHSC1 gene, a detailed search was conducted using the Miranda package (Enright et al., 2003) and the UCSC genome browser (hg18, www.genome.ucsc.edu). The mature miRNA sequence of has-miR-99b (MIMAT0000689) was obtained from mirBase (http://www.mirbase.org; Fig. 1A) and a miR-99b target site in the WHSC1 gene was identified using the Miranda algorithm. The target site in the WHSC1 3′UTR has a strong binding energy to miR-99b and a 3′ 8 nt seed region (Fig. 1B). SNP rs14647 is located at the extreme end of the WHSC1 3′UTR region and is within 50 nt of the miR-99b binding site (Fig. 1C). We used the UCSC ‘PhastCons'track (http://compgen.bscb.cornell.edu/phast/) to evaluate the conservation of the target binding region (Fig. 1D). Seed conservation was observed in some placental mammals which maintained a good binding energy to the WHSC1 3′UTR, indicating that a mutation in the surrounding regions of the target site could disrupt miRNA regulation of WHSC1 protein production and contribute to the development of the disease. Previous mutation studies have shown mutations near miRNA binding sites can affect mRNA stability/translation, presumably by altering the local RNA secondary structure in the usually highly folded 3′UTR region (Doench and Sharp, 2004; Didiano et al., 2010).
We performed haplotype analysis on genes that have more than two SNPs genotyped in our sample set. SNP haplotypes in solute carrier family 22, member 23 (SLC22A23) were shown to be associated with endometriosis subgroups after multiple testing. Therefore, we present detailed information on SLC22A23 only. We typed five miRNA target SNPs in the 3′UTR of SLC22A23 gene located on chromosome 6p25.2. Figure 2 shows the LD plot for the five miRNA SNPs from SLC22A23. The risk allele C of SNP rs7742745 was detected only in 3/958 cases (~0.16%), whereas the observed minor allele frequency of SNP rs35091219 was very low (~0.7%, see above), so both SNPs were excluded from further haplotype analysis. Therefore, we conducted haplotype analysis for the three SNPs: rs3813486, rs1127473 and rs3211066 (Fig. 2). There are four major haplotypes deduced, according to the criteria of Gabriel et al. (2002).
One major haplotype AGC (haplotype 1 in Table IV) was predicted in our all samples, representing 85.2% of alleles, with three minor haplotypes being GAG (haplotype 2 in Table IV) at 7%, GGG (haplotype 3 in Table IV) at 5.8% and AGG (haplotype 4 in Table IV) at 1.6%. The frequencies of the predicted haplotypes for the three SLC22A23 miRNA SNPs listed from highest to lowest in disease groups are shown in Table IV. The major haplotype AGC represented 80.3% of alleles in subfertility cases and 86.1% of alleles in controls. Evidence for association was detected with a corrected P-value of 0.0064 (Table IV). The strongest associations were found for the fourth haplotype, with P-values of 5.0 × 10−4, 1 × 10−4 and 1.3 × 10−3 in the rAFS Stages III/IV group, subfertility and total endometriosis samples when compared with total controls, respectively. The P-values were significant after correcting for multiple testing (P=0.002, 0.002 and 0.005 for Stages III/IV group, subfertility and total cases, respectively), these measured haplotypes account for 2.2% of risk of endometriosis in our population.
The logistic regression analysis was performed based on four functional SNPs: rs14647, rs3813486, rs1127473 and rs3211066, because there was some evidence for association with endometriosis subfertility. SNP rs14647, located on chromosome 4p16.3 in the WHSC1 gene, is clearly not in LD (as it is on a different chromosome) with any of the chromosome 6 SNPs rs3813486, rs1127473 or rs3211066 (SLC22A23). Both SNPs rs14647 and rs3211066 gave the strongest evidence for association with unconditioned P-value of 0.002 and 0.001, respectively. The effects of SNPs rs14647 and rs3211066 on endometriosis-related infertility were completely independent from conditioned logistic regression analysis.
Among the three SNPs in close proximity to each other on chromosome 6 (rs3813486, rs1127473 and rs3211066), SNP rs3211066 showed the strongest evidence for association with endometriosis-related infertility for both unconditional and conditional models. After conditioning on rs3211066, evidence for association at rs3813486 and rs1127473 was no longer significant (P=0.1186, P=0.1178, respectively), indicating their association is not independent and may be due to them being in LD with rs3211066.
We have assessed 41 of the miRNA target SNPs with individual genotyping data, focusing on the molecular pathways they are potentially involved in using IPA. Several significant common pathways overlapping the data set were identified, mainly with TNF and cyclin-dependent kinase inhibitor (CDKN) centrally located in the molecular network (P=6.9E–04 to 3.6E–02). One of these, the cell cycle, cellular growth and proliferation pathway, was predicted to be involved in endometriosis development (Leyendecker et al., 1998) (Fig. 3). IPA identified 12 SNP-related genes in the pathway with molecular networks. The remaining other pathways were significantly enriched in the pathogenesis of endometriosis, including cell morphology, cellular assembly and organization; gene expression, nervous system development and function and carbohydrate metabolism.
Studies in endometriosis have identified several miRNAs differentially expressed in ectopic versus eutopic endometrial tissues (Ohlsson Teague et al., 2009b), endometrial cells (Pan et al., 2007) and in women with endometriosis versus controls (Burney et al., 2009). In general, few miRNAs were detected in more than one study suggesting complex relationships of miRNAs in gene regulatory networks during development of endometrial lesions. Although the relationship between miRNA and mRNA in endometriosis has been evaluated in two studies (Pan et al., 2007; Ohlsson Teague et al., 2009b), a link between miRNAs and effects of DNA variation at their target sites has not been previously investigated. We hypothesized that mutations in target binding sites for the differentially expressed miRNAs could alter transcriptional regulation resulting in changes in expression of the target gene(s) and modifying endometriosis risk.
We used stringent criteria to identify target genes for 22 miRNAs differentially expressed (Ohlsson Teague et al., 2009b) in ectopic and eutopic endometrium. Once target genes were detected, we searched for SNPs in all corresponding sites for each miRNA within the target gene. This search identified 243 SNPs and we genotyped a subset of 102 SNPs in an Australian sample comprising 958 endometriosis cases and 959 controls. Forty-six percent of SNPs were not polymorphic in our sample. Fourteen SNPs showed evidence for rare minor alleles present in no more than three individuals. These were distributed roughly equally between being present in cases only, controls only or in both cases and controls.
Infertility is a major symptom associated with endometriosis and almost half of women with endometriosis are infertile (American Fertility Society, 1985; Guo and Wang, 2006). There is growing evidence that molecular dysfunction in eutopic endometrium contributes to endometriosis susceptibility and may also contribute to endometriosis-related infertility. A number of genes that code for proteins such as integrin (alpha V, beta 3), nitric oxide synthase, heat shock protein and interleukin-11, are reported to have an important role in endometrial receptivity (Agarwal et al., 2005; Bondza et al., 2009). The WHSC1 gene is located on chromosome 4p16.3 and is associated with Wolf–Hirschhorn syndrome. Patients with this gene defect suffer from craniofacial defects, intellectual disability, growth delay and 30% have genital anomalies including anomalies of the Mullerian system (female reproductive system). In a Whsc1-deficient mouse model reported that surviving Whsc1 heterozygotes were mostly fertile, but homozygous embryos showed growth retardation and died within 10 days of birth (Nimura et al., 2009). SNP rs14647 is located near a miRNA target site, between WHSC1 and WHSC2 and in the 3′UTR region of both genes (the two genes are transcribed in opposite directions). UTRs are the regulatory elements of genes, and play important roles in RNA's translation and decay. The WHSC1 and WHSC2 3′UTR region may be involved in dsRNA degradation when both genes are transcribed with overlapping UTRs. The current study is a genetic study in a large case–control sample set aiming to investigate the potential impact of 102 miRNA SNPs in women with endometriosis. Until now, most SNPs identified in genetic studies are not fully confirmed by functional study. However, it is well established that genetic studies not only attempt to identify functional SNPs but also ‘tag’ the approximate location of disease variants (Wang et al., 2010). Another unknown SNP in high LD with rs14647 might either increase or decrease the accessibility of the miRNA target site, leading to inactivation of the target site by this or other unknown mechanisms. A recent study by Halvorsen et al. showed that multiple specific SNPs in UTRs have a structural consequence and can result in a disease phenotype (Halvorsen et al., 2010). We found the rs14647A allele in the WHSC1 gene was most strongly associated with risk of endometriosis-related infertility with an increased frequency up to 13-fold when compared with our controls, suggesting that the SNP may have local effects on the RNA structural consequence of the associated genes.
We also found significant evidence for association between SLC22A23 haplotypes with SNPs in miR-125a-3p binding sites and predisposition to endometriosis. SLC22A23 belongs to a large family of transmembrane proteins that function as uniporters, symporters and antiporters to transport organic ions across cell membranes (Jacobsson et al., 2007). The SLC22A23 target site harbours five miRNA target SNPs: rs3813486, rs1127473, rs7742745, rs35091219 and rs3211066. Mutation in the functional SNP rs7742745 was only detected in three cases and not in 959 controls, suggesting a possible role in risk for endometriosis in these patients. Haplotype analysis for three SNPs (rs3813486, rs1127473 and rs3211066), excluding the rare variants, supports a role for variation in SLC22A23 in endometriosis risk. Haplotypes in SLC22A23 were strongly associated with infertility and more severe stage of the disease. The three SNPs are in one major haplotype block, with the common haplotype AGC representing 80.3% of alleles in subfertility cases and 86.1% of alleles in controls. The minor haplotype AGG accounted for 2.2% of alleles in all cases and 0.9% of alleles in all controls. The association detected with these haplotypes suggests that these variants may play a role in expression of the SLC22A23 transcript levels. The association results remained significant after correction of multiple testing, but need to be replicated in independent samples. Comprehensive genotyping in large, powerful, samples with appropriate data on disease subphenotypes will be required to ultimately confirm and characterise the involvement of the gene variation in the pathogenesis of endometriosis.
To evaluate the miRNA target SNPs that potentially contribute to molecular pathways in endometriosis, we performed pathway analysis using IPA software. By assessing 41 miRNA SNPs with the genotype data, IPA identified several overlapping common pathways constituted by two major molecular networks: TNF and CDKN, both are centrally located in the cellular growth and proliferation pathways. TNF is a pro-inflammatory cytokine and is known to play important roles in the pathological processes involved in endometriosis (Barcz et al., 2000; Braun et al., 2002). CDKN2A is an important tumour suppressor gene and was found to be less expressed in endometriosis lesions when compared with normal controls from immunohistochemical studies (Cameron et al., 2002; Goumenou et al., 2006). The functional relationship between TNF/NF-κB (nuclear factor of kappa light polypeptide gene enhancer in B-cells) signalling and CDKN2A in tumorigenesis has been explored in a mouse model (Altomare et al., 2009), but their functions and relationships in the pathogenesis of endometriosis are unclear. In women with endometriosis, activated macrophages secrete TNF-α, which then activates NF-κB and other inflammatory mediators to promote epithelial cell proliferation and in turn contributes to the disease progress. Two reports from same group provided evidence that inhibition of TNF reduces endometriosis. Treatment with recombinant human TNFRSF1A (tumour necrosis factor receptor superfamily, member 1A) or anti-TNF antibody in baboons, effectively inhibited the development of endometriosis (D'Hooghe et al., 2006; Falconer et al., 2006). TNF-α has been shown to inhibit asbestos-induced cytotoxicity via a NF-κB-dependent pathway (Yang et al., 2006) and NF-κB inhibition reduces cell proliferation, but increases apoptosis of endometriotic lesions (Gonzalez-Ramos et al., 2008). A previous study provided the information that supports the gene-to-function relationship between binding of promoter fragment from the human CDKN1A gene and human WHSC2 protein occurs in Hct 116 cells (Gomes et al., 2006). We found evidence for association with functional miRNA SNPs/SNP haplotypes in WHSC1 and SLC22A23 loci. Both loci have not previously been reported to be associated with endometriosis. Given that these miRNA SNPs impact molecular networks and their potential targets in the pathogenesis of endometriosis, our findings may provide valuable insights for the future investigations.
Although the present work supports our hypothesis that miRNA target SNPs might be involved in the development of endometriosis, it remains to be determined whether the association between these SNPs and endometriosis risk is due to an effect on miRNA binding. Functional studies in ectopic and eutopic endometrial tissues would be required, but the current results need first to be replicated in independent large sample sets to confirm and characterise the involvement of the gene variation in the pathogenesis of endometriosis.
The principal authors’ roles were as follows: Z.Z.: design, acquisition of data, analysis and interpretation, drafting manuscript and final approval. L.C.: conception and design, interpretation, revising for critical content and final approval. D.R.N.: analysis and interpretation of data, revising for critical content and final approval. B.C.: acquisition of data, revising for critical content and final approval. S.A.T.: conception and design, interpretation, revising for critical content and final approval. M.L.H.: conception and design, interpretation, revising for critical content and final approval. G.W.M.: conception and design, interpretation, revising for critical content and final approval.
This study was supported by grants to G.W.M. from the National Institute of Child Health and Human Development (HD050537) and National Health and Medical Research Council of Australia (339430 and 339446).
We thank Dr Daniel T. O'Connor for confirmation of diagnosis and staging of disease from clinical records of 295 of the 958 cases; Barbara Haddon for co-ordination of family recruitment, blood and phenotype collection; Genetic and Molecular Epidemiology Laboratories staff for sample processing and DNA preparation.