- Research article
- Open Access
Genetic variations of MUC17 are associated with endometriosis development and related infertility
BMC Medical Genetics volume 16, Article number: 60 (2015)
Genetic alterations of mucin genes, such as MUC2 and MUC4, were previously identified to be associated with endometriosis and related infertility. Additionally, gene expression profiling has confirmed MUC17 to be overexpressed in mucinous ovarian carcinoma; however, its associated risk for endometriosis remains unclear. This study was focused on the potential impact of genetic variations in MUC17 on endometriosis development and associated clinical features.
The study subjects included 189 female Taiwanese patients with pathology-proven endometriosis and 191 healthy Taiwanese women as controls. Five single-nucleotide polymorphisms (rs4729645, rs10953316, rs74974199, rs4729655, and rs4729656) within the MUC17 gene were selected and genotyped using the Taqman genotyping assay to examine the allele frequency and genotype distributions of MUC17 polymorphisms.
Genotyping revealed that the A allele at rs10953316 in MUC17 was a protective genetic factor in endometriosis development (p = 0.008; OR = 0.53; 95 % CI: 0.36-0.79). Genetic variation of rs4729655 protected patients against endometriosis-induced infertility, but was associated with a higher cancer antigen 125 (CA125) level. Base-pairing analysis, called MaxExpect, predicted an additional loop in the mRNA structure caused by rs10953316 polymorphism, possibly influencing ribosome sliding and translation efficiency. Such predictions were confirmed by immunohistochemistry that patients with AA genotype at rs10953316 showed low MUC17 levels in their endometrium, patients with GA genotype showed moderate levels, and strong staining could be found in patients with GG genotype.
MUC17 polymorphisms are involved in endometriosis development and the associated infertility in the Taiwanese population.
Endometriosis is a common chronic gynecological disease described as the presence of endometrial glands and stroma located outside the uterine cavity, ovaries, fallopian tubes, and even on the bladder or intestines [1, 2]. It occurs in approximately 10 % of women or teenage girls at reproductive age. Up to 50 % of patients also suffer from infertility [3, 4]. Clinical symptoms of endometriosis include several types of pain, such as excessive menstrual pain, pelvic pain with defecation, and chronic pelvic pain . Epidemiological studies revealed a higher risk for endometriosis patients to develop different types of ovarian cancers [6, 7]. Molecular pathological analyses have provided strong evidence to support the histological transition from benign endometriosis to ovarian malignancy [8, 9].
Mucins are high-molecular-weight o-glycoproteins that function as the protective and lubricative layers on epithelial surfaces, such as the respiratory, gastrointestinal, and reproductive tracts [10–13]. Gene expression analyses have indicated that MUC1 ~ 4 are the major mucins constitutively expressed in endometrial epithelium, and their levels can be controlled by the menstrual cycle [14, 15]. With potent roles in cellular proliferation, MUC1 and MUC4 overexpression has been identified as a tumor marker for endometriotic lesions and certain types of ovarian tumors [16, 17]. Interestingly, genetic polymorphisms in MUC2 and MUC4 were recently reported to be associated with endometriosis development and the related infertility [18, 19]. Because endometriosis shares several similarities with cancer (e.g., highly mobile and proliferative) these data support that alterations in mucins either at the expression level or within their genetic sequences may have a significant impact on endometrium regulation and development.
In this study, we investigated the association between genetic variations of MUC17 (A.K.A. MUC3) and endometriosis development based on the finding of gene expressing profiling in mucinous ovarian carcinoma (MOC), which likely arises from endometriosis [20, 21]. In this study, the mucin genes MUC2, MUC3A, and MUC17 were found to be upregulated in MOC . Similar to MUC4, MUC17 encodes a membrane-bound mucin containing two EGF (epidermal growth factor)-like domains and a SEA (sea urchin sperm protein, enterokinase and agrin) domain in the extracellular portion, which can interact with human EGF receptor 2 (EGFR2) on the adjacent cell surface. Thus, the growth signal can be transduced to stimulate cellular proliferation . Of note, MUC17 overexpression can also be found in colon and pancreatic cancers . Although several lines of evidence support a potent role of MUC17 in tumor development, its role in endometriosis remains unclear. To investigate the possible link with endometriosis, we studied the association of MUC17 genetic variants with the ovarian cancer biomarker CA125 and the related infertility.
One hundred eighty-nine Taiwanese patients who underwent surgery for benign ovarian disease and pathology-proven endometriosis at China Medical University Hospital from 1998 to 2011 were enrolled. These patients were diagnosed with ovarian cysts by ultrasound and confirmed by surgical operation. Symptoms of infertility and CA125 levels were also recorded for each patient during the study. Patients who failed to have pathology-proven endometriosis were excluded from this study. Another 191 healthy women were collected as controls who had regular health checkups at the same hospital. Ultrasound assessment indicated healthy conditions in the abdominal cavities of these controls. Both the controls and patients shared similar age profiles. This study was approved by the Institutional Review Board of China Medical University Hospital. Written informed consent was obtained from patients for scientific research and publication of this study.
Genomic DNA was extracted from peripheral blood leukocytes according to the manufacturer’s protocol (Genomic DNA kit; Qiagen, Valencia, CA, USA). DNA fragments containing the target single-nucleotide polymorphism (SNP) sites were amplified by polymerase chain reaction (PCR) using the Taqman SNP Genotyping Assay System from Applied Biosystems, Inc. (Carlsbad, CA, USA). The probe search and design are available from the ABI SNP genotyping databank and are listed in Additional file 1. The PCR amplification conditions were as follows: 95 °C for 5 minutes, followed by 40 cycles at 95 °C for 10 seconds, 56 °C for 10 seconds and 72 °C for 20 seconds, with one additional cycle at 72 °C for 5 minutes. Genetic variations were detected by reading fluorescence signals of the PCR products. A positive signal indicates a perfect match between the probe and tested DNA, thus identifying the allele type.
The allelic and genotypic frequency distributions for five SNP sites in endometriosis patients and controls were performed by χ2 analysis using SPSS software (version 10.0; SPSS Inc. Chicago, IL, USA). Allelic and genotypic frequencies were expressed as percentages of the total number of alleles and genotypes. Odds ratios (ORs) were calculated for the allelic and genotypic frequencies with 95 % confident intervals (95 % CIs). The major (wild-type) allele was used as the reference group. Lewontin's D' (|D'|) and the linkage disequilibrium (LD) coefficient r2 were determined between selected pairs of biallelic loci. Haploview version 3.2 software (Broad Institute, Cambridge, MA) was used to examine the structure of the LD block.
Prediction of mRNA secondary structures adjacent to the SNP sites
The sequence of rs10953316, as well as the reference sequence, were submitted to the RNAstructure web server (http://rna.urmc.rochester.edu/RNAstructureWeb/index.html) , and structures with maximum expected accuracy were selected. Among the multiple algorithms in the RNAstructure web server, MaxExpect was utilized to generate structures composed of highly probable base pairs with the highest probability of being accurate.
A total of 76 ectopic endometrial tissues were collected from patients with different genotypes at the rs10953316 SNP site. Rabbit anti-MUC17 polyclonal antibody (ab122184) was purchased from Abcam PLC., Cambridge, USA. Immunohistochemistry was performed with an antibody dilution of 1:50 using standard protocols with an EnVision + System peroxidase kit (DAKO Agilent Technologies, Glostrup, Denmark). Immunointensity of MUC17 in endometrial epithelium was independently scored by two pathologists and specific membrane staining was labeled as: weak positive, moderately positive, or strongly positive. For discordant cases, a third investigator was brought in to score and the final intensity score was determined by the majority scores.
MUC17 polymorphisms and endometriosis
To examine whether SNPs in MUC17 play a role in endometriosis, five SNPs were selected for this study, with a frequency greater than 5 % in the Chinese Han in Beijing from the International HapMap Project databank. The SNPs locations in MUC17 are indicated in Fig. 1a. These SNPs were either located in the open reading frame region or the 3′-untranslated region (3′-UTR) of MUC17. The LD values of selected SNPs in patients are shown in Fig. 1B. Genotype analysis and allele distributions of SNPs in MUC17 in the patients and controls are analyzed in Table 1. The major (wild-type) alleles were used as the references. Genetic analysis showed a significant association between endometriosis and genetic variation at the SNP site of rs10953316 with a lower frequency of the A allele in patients than in controls (p = 0.008; OR = 0.53; 95 % CI: 0.36-0.79) (Table 1). Genotype analysis of the rs10953316 site supported the protective potential of the AG genotype in endometriosis development compared with the reference GG genotype (OR = 0.45; 95 % CI: 0.29-0.7). The AA genotype was not found to show a dominant effect due to the limited sample size. Allelic and genotype analyses of the other SNPs showed no statistically significant differences.
MUC17 polymorphisms and infertility
Because endometriosis is considered one of the factors that cause infertility, we next studied the relationship between MUC17 SNPs and infertility. The endometriosis patients were sub-grouped into the fertile and infertile groups for allele distribution analysis. The results shown in Table 2 indicate that the genetic variant rs4729655 is a protective genetic factor for infertility in patients (p = 0.0465; OR = 0.42; 95 % CI: 0.21-0.82). The genetic variant at rs4729656 also showed a protective tendency, although this result did not reach statistical significance after Bonferroni correction.
MUC17 polymorphism and CA125 levels
Cancer antigen 125 (CA125), also known as MUC16, is a commonly used serum biomarker to monitor the progress and recurrence of ovarian malignancy and benign conditions, including endometriosis. It has been shown that serum CA125 levels are markedly elevated in women with cystic ovarian endometriosis and/or deeply infiltrating endometriosis compared with women in the luteal phase with minimal or mild endometriosis . Therefore, we asked whether SNPs in MUC17 can affect the serum CA125 levels in patients. Patients were sub-grouped into two groups for genetic analyses based on the CA125 level: high CA125 (≥35 KU/L) and low CA125 (<35 KU/L) groups [27, 28]. Our data indicated that patients with the rs4729655 SNP showed a significantly high expression level of CA125 (p = 0.038; OR = 0.48; 95 % CI: 0.28-0.83). Genetic variation at rs74974199 also showed a correlation with elevated CA125 expression levels; however, such correlation did not reach statistical significance after Bonferroni correction. By contrast, patients with genetic variation at rs4729656 showed a tendency toward a lower CA125 level (Table 3).
MUC17 rs10953316 polymorphism controls its protein expression level
Based on above findings, the rs10953316 polymorphism was defined as a unique genetic factor regulating endometriosis development. However, this SNP is a silent genetic variation resulting in no amino acid change (Thr2355Thr). Therefore, we asked whether such a SNP can induce an RNA structural conformational change using the MaxExpect algorithm from the RNAstructure web server . MaxExpect analysis revealed that the mRNA with genetic variation at rs10953316 forms an additional loop adjacent to the SNP site (Fig. 2a). Because mRNA structural change could be a hindrance to its stability and translation efficiency , we performed immunohistochemistry on endometriosis tissues to see whether such variation can influence the total amount of MUC17. As shown in Fig. 2B, patients with AA genotype at rs10953316 showed low MUC17 levels in their endometrium, patients with GA genotype showed moderate levels, and strong staining could be found in patients with GG genotype.
In this study, we showed molecular evidence that the A allele of rs10953316 in MUC17 may function as a protective factor in endometriosis development. One of our limitations is the relatively small sample size, so that genotype AA at rs10953316 did not show a stronger effect (OR = 0.77, 95 % CI: 0.05-12.65) than the genotype AG (OR = 0.45; 95 % CI: 0.29-0.7). Although the rs10953316 polymorphism is a silent genetic variation, the total protein amount of MUC17 can be influenced probably due to RNA structural conformational change (Fig. 2). Our finding supports the previous hypothesis that structural change in RNA coding region can lead to ribosome pause and affect translational efficiency . Similar finding was also evidenced by Nackley et al. that haplotypes divergent in synonymous changes in human COMT gene alter the mRNA structure with additional loops and reduced the amount of translated protein, which subsequently affected its function such as pain perception . Because the extracellular domain of MUC17 protein contains several EGF-like functional domains, reduction of MUC17 expression by rs10953316 SNP may down-regulate the erbB signaling leading to inhibition of cell proliferation . In addition, SEA domain at the C-terminal region of MUC17 is critical for its autoproteolysis, which controls cell migration and invasion . These studies provide clues for us to presume that rs10953316 SNP plays a protective role in endometriosis development by reducing total MUC17 expression and the downstream signaling, whereas the detailed mechanism remains to be further investigated.
Another SNP, rs4729655, was linked to endometriosis-related infertility and was associated with abnormally high levels of CA125. The rs4729655 SNP site is located at the 3′-UTR of MUC17, a regulatory region for mRNA stability and miRNA binding. We next analyzed whether there was any possible miRNA targets gained or lost due to genetic variation at the rs4729655 SNP site. Our study found that genetic variation at rs4729655 resulted in the gain of miR-4508 and loss of miR-3158-3p toward the 3′-UTR of MUC17; thus, this SNP site may be involved in multiple miRNA interactions and gene regulation. Our previous study showed that a miR-125b binding site variation at the 3′-UTR of BMPR1B defines endometriosis susceptibility and CA125 levels . Currently, biofunctions and downstream targets of miR-4508 and miR-3158-3p are still largely unknown. However, studies on the miRNAome of different cancer types using a deep sequencing approach provided solid evidence indicating the involvement of these two miRNAs in cancer development [33–36]. Particularly, miR-3158-3p expression was found to be associated with disease states in the female reproductive tract, including ovarian cancer and endometriosis . Whether dysregulation of these two miRNAs could lead to infertility remains unclear. Further studies on these two miRNAs may provide insights to uncover the detailed mechanisms regarding how they regulate the development of female reproductive diseases.
One obstacle of this study was that the SNP associated with endometriosis development did not show an influence on infertility or CA125 levels in patients, and vice versa. Our data may reflect that more than 50 % of patients with endometriosis-related infertility are diagnosed as having minimal or mild clinical symptoms . By contrast, the sensitivity and specificity of the CA125 blood test was also found to be limited, not necessary correlating with clinical stages. Therefore, it was recommended as an indication that women with high CA125 levels should receive further diagnostic screening and treatment [38, 39]. Cytokine/hormone networking, good oocyte quality and the success of embryo implantation may be the key players to determine endometriosis-related infertility. Because mucins have been found to play roles in regulating the immune responses of the reproductive tract and cellular proliferation of endometrial cells, whether MUC17 also plays such a roles is an interesting topic for us to further investigate.
In this study, we analyzed and found the association between MUC17 polymorphisms and endometriosis, endometriosis-related infertility and CA125 in a Taiwanese population. However, the relatively small sample size is one limitation of this study, and the conclusions should be further validated by another study cohort. The detailed regulatory mechanisms regarding how MUC17 functions and collaborates with its binding miRNAs in promoting endometriosis development remain to be elucidated.
cancer antigen 125
mucinous ovarian cancer
Kinkel K, Frei KA, Balleyguier C, Chapron C. Diagnosis of endometriosis with imaging: a review. European radiology. 2006;16(2):285–98.
Jenkins S, Olive DL, Haney AF. Endometriosis: pathogenetic implications of the anatomic distribution. Obstetrics and gynecology. 1986;67(3):335–8.
Ballard K, Lowton K, Wright J. What's the delay? A qualitative study of women's experiences of reaching a diagnosis of endometriosis. Fertility and sterility. 2006;86(5):1296–301.
Eskenazi B, Warner ML. Epidemiology of endometriosis. Obstetrics and gynecology clinics of North America. 1997;24(2):235–58.
Berkley KJ, Rapkin AJ, Papka RE. The pains of endometriosis. Science. 2005;308(5728):1587–9.
Van Gorp T, Amant F, Neven P, Vergote I, Moerman P. Endometriosis and the development of malignant tumours of the pelvis. A review of literature. Best practice & research Clinical obstetrics & gynaecology. 2004;18(2):349–71.
Orezzoli JP, Russell AH, Oliva E, Del Carmen MG, Eichhorn J, Fuller AF. Prognostic implication of endometriosis in clear cell carcinoma of the ovary. Gynecologic oncology. 2008;110(3):336–44.
Kobayashi H, Kajiwara H, Kanayama S, Yamada Y, Furukawa N, Noguchi T, et al. Molecular pathogenesis of endometriosis-associated clear cell carcinoma of the ovary (review). Oncology reports. 2009;22(2):233–40.
Nezhat F, Datta MS, Hanson V, Pejovic T, Nezhat C, Nezhat C. The relationship of endometriosis and ovarian malignancy: a review. Fertility and sterility. 2008;90(5):1559–70.
Yonezawa S, Higashi M, Yamada N, Yokoyama S, Goto M. Significance of mucin expression in pancreatobiliary neoplasms. Journal of hepato-biliary-pancreatic sciences. 2010;17(2):108–24.
Groneberg DA, Peiser C, Dinh QT, Matthias J, Eynott PR, Heppt W, et al. Distribution of respiratory mucin proteins in human nasal mucosa. The Laryngoscope. 2003;113(3):520–4.
Filipe MI. Mucins in the human gastrointestinal epithelium: a review. Investigative & cell pathology. 1979;2(3):195–216.
Gipson IK, Ho SB, Spurr-Michaud SJ, Tisdale AS, Zhan Q, Torlakovic E, et al. Mucin genes expressed by human female reproductive tract epithelia. Biology of reproduction. 1997;56(4):999–1011.
Gollub EG, Goswami S, Kouba D, Marom Z. Regulation of mucin gene expression in secretory epithelial cells. Biochemical and biophysical research communications. 1993;197(2):667–73.
Audie JP, Tetaert D, Pigny P, Buisine MP, Janin A, Aubert JP, et al. Mucin gene expression in the human endocervix. Human reproduction. 1995;10(1):98–102.
Vlad AM, Diaconu I, Gantt KR. MUC1 in endometriosis and ovarian cancer. Immunologic research. 2006;36(1–3):229–36.
Nath S, Mukherjee P. MUC1: a multifaceted oncoprotein with a key role in cancer progression. Trends in molecular medicine. 2014;20(6):332–42.
Chang CY, Chang HW, Chen CM, Lin CY, Chen CP, Lai CH, et al. MUC4 gene polymorphisms associate with endometriosis development and endometriosis-related infertility. BMC medicine. 2011;9:19.
Chang CY, Chen Y, Lin WC, Chen CM, Chen CP, Lee SC, et al. MUC2 polymorphisms are associated with endometriosis development and infertility: a case–control study. BMC medical genetics. 2012;13:15.
Ko ML, Pan HS. Identical twins with ovarian endometriosis and mucinous borderline tumor: an unusual association. Fertility and sterility. 2008;90(4):1200. e1207-1209.
Lee KR, Nucci MR. Ovarian mucinous and mixed epithelial carcinomas of mullerian (endocervical-like) type: a clinicopathologic analysis of four cases of an uncommon variant associated with endometriosis. International journal of gynecological pathology: official journal of the International Society of Gynecological Pathologists. 2003;22(1):42–51.
Heinzelmann-Schwarz VA, Gardiner-Garden M, Henshall SM, Scurry JP, Scolyer RA, Smith AN, et al. A distinct molecular profile associated with mucinous epithelial ovarian cancer. British journal of cancer. 2006;94(6):904–13.
Bafna S, Kaur S, Batra SK. Membrane-bound mucins: the mechanistic basis for alterations in the growth and survival of cancer cells. Oncogene. 2010;29(20):2893–904.
Hattrup CL, Gendler SJ. Structure and function of the cell surface (tethered) mucins. Annual review of physiology. 2008;70:431–57.
Reuter JS, Mathews DH. RNAstructure: software for RNA secondary structure prediction and analysis. BMC bioinformatics. 2010;11:129.
Muyldermans M, Cornillie FJ, Koninckx PR. CA125 and endometriosis. Human reproduction update. 1995;1(2):173–87.
Nossov V, Amneus M, Su F, Lang J, Janco JM, Reddy ST, et al. The early detection of ovarian cancer: from traditional methods to proteomics. Can we really do better than serum CA-125? American journal of obstetrics and gynecology. 2008;199(3):215–23.
Bagan P, Berna P, Assouad J, Hupertan V, Le Pimpec BF, Riquet M. Value of cancer antigen 125 for diagnosis of pleural endometriosis in females with recurrent pneumothorax. The European respiratory journal. 2008;31(1):140–2.
Mortimer SA, Kidwell MA, Doudna JA. Insights into RNA structure and function from genome-wide studies. Nature reviews Genetics. 2014;15(7):469–79.
Nackley AG, Shabalina SA, Tchivileva IE, Satterfield K, Korchynskyi O, Makarov SS, et al. Human catechol-O-methyltransferase haplotypes modulate protein expression by altering mRNA secondary structure. Science. 2006;314(5807):1930–3.
Peng Z, He Y, Yang Y, Zhu R, Bai J, Li Y, et al. Autoproteolysis of the SEA module of rMuc3 C-terminal domain modulates its functional composition. Archives of biochemistry and biophysics. 2010;503(2):238–47.
Chang CY, Chen Y, Lai MT, Chang HW, Cheng J, Chan C, et al. BMPR1B up-regulation via a miRNA binding site variation defines endometriosis susceptibility and CA125 levels. PloS one. 2013;8(12):e80630.
Jima DD, Zhang J, Jacobs C, Richards KL, Dunphy CH, Choi WW, et al. Deep sequencing of the small RNA transcriptome of normal and malignant human B cells identifies hundreds of novel microRNAs. Blood. 2010;116(23):e118–27.
Stark MS, Tyagi S, Nancarrow DJ, Boyle GM, Cook AL, Whiteman DC, et al. Characterization of the Melanoma miRNAome by Deep Sequencing. PloS one. 2010;5(3):e9685.
Creighton CJ, Benham AL, Zhu HF, Khan MF, Reid JG, Nagaraja AK, et al. Discovery of Novel MicroRNAs in Female Reproductive Tract Using Next Generation Sequencing. PloS one. 2010;5(3):e9637.
Persson H, Kvist A, Rego N, Staaf J, Vallon-Christersson J, Luts L, et al. Identification of New MicroRNAs in Paired Normal and Tumor Breast Tissue Suggests a Dual Role for the ERBB2/Her2 Gene. Cancer Res. 2011;71(1):78–86.
Holoch KJ, Lessey BA. Endometriosis and Infertility. Clin Obstet Gynecol. 2010;53(2):429–38.
Baron JA. Screening for cancer with molecular markers: progress comes with potential problems. Nat Rev Cancer. 2012;12(5):368–71.
Fritsche HA, Bast RC. CA 125 in ovarian cancer: Advances and controversy. Clin Chem. 1998;44(7):1379–80.
Authors would like to thank Drs. Lei Wan, Ying-Ju Lin and Shh-Ying Chen at Human Genetic Center, China Medical University Hospital, Taichung, Taiwan for their critical review and suggestions. This work was funded by grant DMR-101-071 from China Medical University Hospital, Taichung, Taiwan.
The authors declare that they have no competing interests.
CWY, YC, PWT and SHY carried out molecular genetic studies. CYYC and MTL collected clinical samples and performed statistical analysis of clinical features. HWC, CMC and SCL performed statistical analysis of genetic studies. MTL and CHL performed IHC staining and scoring. CCL, JJCS and FJT designed the experiments. CWY, CYYC, JJCS and FJT drafted the tables, figures and manuscript. All authors read and approved the final manuscript.
Ching-Wen Yang and Cherry Yin-Yi Chang contributed equally to this work.
About this article
Cite this article
Yang, CW., Chang, C.YY., Lai, MT. et al. Genetic variations of MUC17 are associated with endometriosis development and related infertility. BMC Med Genet 16, 60 (2015). https://doi.org/10.1186/s12881-015-0209-7
- CA125 Level
- Endometriosis Patient
- Ovarian Endometriosis