Endometrial vezatin and its association with endometriosis risk

Sarah J. Holdsworth-Carson, Jenny N. Fung, Hien T.T. Luong, Yadav Sapkota, Lisa M. Bowdler, Leanne Wallace, Wan Tinn Teh, Joseph E. Powell, Jane E. Girling, Martin Healey, Grant W. Montgomery, Peter A.W. Rogers

Do endometriosis risk-associated single nucleotide polymorphisms (SNPs) found at the 12q22 locus have effects on vezatin (VEZT) expression? The original genome-wide association study (GWAS) SNP (rs10859871), and other newly identified association signals, demonstrate strong evidence for cis-expression quantitative trait loci (eQTL) effects on VEZT expression. GWAS have identified several disease-risk loci (SNPs) associated with endometriosis. The SNP rs10859871 is located within the VEZT gene. VEZT expression is altered in the endometrium of endometriosis patients and is an excellent candidate for having a causal role in endometriosis. Most of the SNPs identified from GWAS are not located within the coding region of the genome. However, they are likely to have an effect on the regulation of gene expression. Genetic variants that affect levels of gene expression are called expression quantitative trait loci (eQTL). Samples for genotyping and VEZT variant screening were drawn from women recruited for genetic studies in Australia/New Zealand and women undergoing surgery in a tertiary care centre. Coding variants for VEZT were screened in blood from 100 unrelated individuals (endometriosis-dense families) from the QIMR Berghofer Medical Research Institute dataset. SNPs at the 12q22 locus were imputed and reanalysed for their association with endometriosis. Reanalysis of endometriosis risk-association was performed on a final combined Australian dataset of 2594 cases and 4496 controls. Gene expression was performed on 136 endometrial samples. eQTL analysis in whole blood was performed on 862 individuals from the Brisbane Systems Genetics Study. Endometrial tissue-specific eQTL analysis was performed on 122 samples (eutopic endometrium) collected following laparoscopic surgery. VEZT protein expression studies employed n = 56 (western blotting) and n = 42 (immunohistochemistry) endometrial samples. The women recruited for this study provided blood and/or endometrial tissue samples in a hospital setting. Genomic DNA was screened for common and coding variants. SNPs of interest in the 12q22 region were genotyped using Agena MassARRAY technology or Taqman SNP genotyping assay. Gene expression profiles from RNA extracted from blood and endometrial tissue samples were generated using Illumina whole-genome expression chips (Human HT-12 v4.0). Whole protein extracted from endometrium was used for VEZT western blots, and paraffin sections of endometrium were employed for VEZT immunohistochemistry semi-quantitative analysis. A total of 11 coding variants of VEZT (including one novel variant) were identified from an endometriosis-dense cohort. Polymorphic coding and imputed SNPs were combined with previous GWAS data to reanalyse the endometriosis risk association of the 12q22 region. The disease association signal at 12q22 was due to coding variants in VEZT or FGD6 (FYVE, RhoGEF and PH domain-containing 6) and SNPs with the strongest signals were either intronic or intergenic. We found strong evidence for VEZT cis-eQTLs with the sentinel SNP (rs10859871) in blood and endometrium, where the endometriosis risk allele (C) was associated with an increase in VEZT expression. We could not demonstrate this genotype-specific effect on VEZT protein expression in endometrium. However, we did observe a menstrual cycle stage specific increase in VEZT protein expression in endometrial glands, specific to the secretory phase (P = 2.0 × 10(-4)). In comparison to the blood sample datasets, the study numbers of endometrial tissues were substantially reduced. Protein studies failed to complement RNA results, also likely a reflection of the low study numbers in these experiments. In silico prediction tools used in this investigation are typically based on cell lines different to our tissues of interest, thus any functional annotations drawn from these approaches should be considered carefully. Therefore, functional studies on VEZT and related pathway components are still warranted to unequivocally implicate a causal role for VEZT in endometriosis pathophysiology. GWAS have proven to be very valuable tools for deciphering complex diseases. Endometriosis is a text-book example of a complex disease, involving genetic, lifestyle and environmental influences. Our focused investigation of the 12q22 region validates an association with increased endometriosis risk. Endometriosis risk SNPs (including rs10859871) located within this locus demonstrated evidence for cis-eQTLs on VEZT expression. By examining women who possess an enhanced genetic risk of developing endometriosis, we have identified an effect on VEZT expression and therefore a potential gene/gene pathway in endometriosis disease establishment and development. Funding for this work was provided by NHMRC Project Grants GNT1012245, GNT1026033, GNT1049472 and GNT1046880. G.W.M. is supported by the NHMRC Fellowship scheme (GNT1078399). S.J.H.-C. is supported by the J.N. Peters Bequest Fellowship. The authors declare no competing interests. N/A.