Yadav Sapkota

Association between endometriosis and the interleukin 1A (IL1A) locus

STUDY QUESTION Are single-nucleotide polymorphisms (SNPs) at the interleukin 1A (IL1A) gene locus associated with endometriosis risk? SUMMARY ANSWER We found evidence for strong association between IL1A SNPs and endometriosis risk. WHAT IS KNOWN ALREADY Genetic factors contribute substantially to the complex aetiology of endometriosis and the disease has an estimated heritability of ∼51%. We, and others, have conducted genome-wide association (GWA) studies for endometriosis, which identified a total of nine independent risk loci. Recently, two small Japanese studies reported eight SNPs (rs6542095, rs11677416, rs3783550, rs3783525, rs3783553, rs2856836, rs1304037 and rs17561) at the IL1A gene locus as suggestively associated with endometriosis risk. There is also evidence of a link between inflammation and endometriosis. STUDY DESIGN, SIZE, DURATION We sought to further investigate the eight IL1A SNPs for association with endometriosis using an independent sample of 3908 endometriosis cases and 8568 controls of European and Japanese ancestry. The study was conducted between October 2013 and July 2014. PARTICIPANTS/MATERIALS, SETTING, METHODS By leveraging GWA data from our previous multi-ethnic GWA meta-analysis for endometriosis, we imputed variants in the IL1A region, using a recent 1000 Genomes reference panel. After combining summary statistics for the eight SNPs from our European and Japanese imputed data with the published results, a fixed-effect meta-analysis was performed. An additional meta-analysis restricted to endometriosis cases with moderate-to-severe (revised American Fertility Society stage 3 or 4) disease versus controls was also performed. MAIN RESULTS AND THE ROLE OF CHANCE All eight IL1A SNPs successfully replicated at P < 0.014 in the European imputed data with concordant direction and similar size to the effects reported in the original Japanese studies. Of these, three SNPs (rs6542095, rs3783550 and rs3783525) also showed association with endometriosis at a nominal P < 0.05 in our independent Japanese sample. Fixed-effect meta-analysis of the eight SNPs for moderate-to-severe endometriosis produced a genome-wide significant association for rs6542095 (odds ratio = 1.21; 95% confidence interval = 1.13-1.29; P = 3.43 × 10(-8)). LIMITATIONS, REASONS FOR CAUTION The meta-analysis for moderate-to-severe endometriosis included results of moderate-to-severe endometriosis cases from our European data sets and all endometriosis cases from the Japanese data sets, as disease stage information was not available for endometriosis cases in the Japanese data sets. WIDER IMPLICATIONS OF THE FINDINGS SNP rs6542095 is located ∼2.3 kb downstream of the IL1A gene and ∼6.9 kb upstream of cytoskeleton-associated protein 2-like (CKAP2L) gene. The IL1A gene encodes the IL1a protein, a member of the interleukin 1 cytokine family which is involved in various immune responses and inflammatory processes. These results provide important replication in an independent Japanese sample and, for the first time, association of the IL1A locus in endometriosis patients of European ancestry. SNPs within the IL1A locus may regulate other genes, but if IL1A is the target, our results provide supporting evidence for a link between inflammatory responses and the pathogenesis of endometriosis. STUDY FUNDING/COMPETING INTERESTS The research was funded by grants from the Australian National Health and Medical Research Council and Wellcome Trust. None of the authors has competing interests for the study.

Endometriosis risk alleles at 1p36.12 act through inverse regulation of CDC42 and LINC00339

Genome-wide association studies (GWAS) have identified markers within the WNT4 region on chromosome 1p36.12 showing consistent and strong association with increasing endometriosis risk. Fine mapping using sequence and imputed genotype data has revealed strong candidates for the causal SNPs within these critical regions; however, the molecular pathogenesis of these SNPs is currently unknown. We used gene expression data collected from whole blood from 862 individuals and endometrial tissue from 136 individuals from independent populations of European descent to examine the mechanism underlying endometriosis susceptibility. Association mapping results from 7,090 individuals (2,594 cases and 4,496 controls) supported rs3820282 as the SNP with the strongest association for endometriosis risk (P = 1.84 × 10−5, OR = 1.244 (1.126-1.375)). SNP rs3820282 is a significant eQTL in whole blood decreasing expression of LINC00339 (also known as HSPC157) and increasing expression of CDC42 (P = 2.0 ×10−54 and 4.5x10−4 respectively). The largest effects were for two LINC00339 probes (P = 2.0 ×10−54; 1.0 × 10−34). The eQTL for LINC00339 was also observed in endometrial tissue (P = 2.4 ×10−8) with the same direction of effect for both whole blood and endometrial tissue. There was no evidence for eQTL effects for WNT4. Chromatin conformation capture provides evidence for risk SNPs interacting with the promoters of both LINC00339 and CDC4 and luciferase reporter assays suggest the risk SNP rs12038474 is located in a transcriptional silencer for CDC42 and the risk allele increases expression of CDC42. However, no effect of rs3820282 was observed in the LINC00339 expression in Ishikawa cells. Taken together, our results suggest that SNPs increasing endometriosis risk in this region act through CDC42, but further functional studies are required to rule out inverse regulation of both LINC00339 and CDC42.

Endometrial vezatin and its association with endometriosis risk

STUDY QUESTION Do endometriosis risk-associated single nucleotide polymorphisms (SNPs) found at the 12q22 locus have effects on vezatin ( ITALIC! VEZT) expression? SUMMARY ANSWER The original genome-wide association study (GWAS) SNP (rs10859871), and other newly identified association signals, demonstrate strong evidence for ITALIC! cis-expression quantitative trait loci (eQTL) effects on ITALIC! VEZT expression. WHAT IS KNOWN ALREADY GWAS have identified several disease-risk loci (SNPs) associated with endometriosis. The SNP rs10859871 is located within the ITALIC! VEZT gene. ITALIC! 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). STUDY DESIGN, SIZE, DURATION Samples for genotyping and ITALIC! 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 ITALIC! 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 ITALIC! n = 56 (western blotting) and ITALIC! n = 42 (immunohistochemistry) endometrial samples. PARTICIPANTS/MATERIALS, SETTING, METHODS 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. MAIN RESULTS AND THE ROLE OF CHANCE A total of 11 coding variants of ITALIC! 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 ITALIC! VEZT or ITALIC! FGD6 (FYVE, RhoGEF and PH domain-containing 6) and SNPs with the strongest signals were either intronic or intergenic. We found strong evidence for ITALIC! VEZT cis-eQTLs with the sentinel SNP (rs10859871) in blood and endometrium, where the endometriosis risk allele (C) was associated with an increase in ITALIC! 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 ( ITALIC! P = 2.0 × 10(-4)). LIMITATIONS, REASONS FOR CAUTION 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. ITALIC! 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. WIDER IMPLICATIONS OF THE FINDINGS 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 ITALIC! cis-eQTLs on ITALIC! VEZT expression. By examining women who possess an enhanced genetic risk of developing endometriosis, we have identified an effect on ITALIC! VEZT expression and therefore a potential gene/gene pathway in endometriosis disease establishment and development. STUDY FUNDING/COMPETING INTERESTS 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. TRIAL REGISTRATION NUMBER N/A.

The genetic regulation of transcription in human endometrial tissue

Study question Do genetic effects regulate gene expression in human endometrium? Summary answer This study demonstrated strong genetic effects on endometrial gene expression and some evidence for genetic regulation of gene expression in a menstrual cycle stage-specific manner. What is known already Genetic effects on expression levels for many genes are tissue specific. Endometrial gene expression varies across menstrual cycle stages and between individuals, but there are limited data on genetic control of expression in endometrium. Study design, size, duration We analysed genome-wide genotype and gene expression data to map cis expression quantitative trait loci (eQTL) in endometrium. Participants/materials, setting, methods We recruited 123 women of European ancestry. DNA samples from blood were genotyped on Illumina HumanCoreExome chips. Total RNA was extracted from endometrial tissues. Whole-transcriptome profiles were characterized using Illumina Human HT-12 v4.0 Expression Beadchips. We performed eQTL mapping with ~8 000 000 genotyped and imputed single nucleotide polymorphisms (SNPs) and 12 329 genes. Main results and the role of chance We identified a total of 18 595 cis SNP-probe associations at a study-wide level of significance (P < 1 × 10-7), which correspond to independent eQTLs for 198 unique genes. The eQTLs with the largest effect in endometrial tissue were rs4902335 for CHURC1 (P = 1.05 × 10-32) and rs147253019 for ZP3 (P = 8.22 × 10-30). We further performed a context-specific eQTL analysis to investigate if genetic effects on gene expression regulation act in a menstrual cycle-specific manner. Interestingly, five cis-eQTLs were identified with a significant stage-by-genotype interaction. The strongest stage interaction was the eQTL for C10ORF33 (PYROXD2) with SNP rs2296438 (P = 2.0 × 10-4), where we observe a 2-fold difference in the average expression levels of heterozygous samples depending on the stage of the menstrual cycle. Large scale data The summary eQTL results are publicly available to browse or download. Limitations, reasons for caution A limitation of the present study was the relatively modest sample size. It was not powered to identify trans-eQTLs and larger sample sizes will also be needed to provide better power to detect cis-eQTLs and cycle stage-specific effects, given the substantial changes in expression across the menstrual cycle for many genes. Wider implications of the findings Identification of endometrial eQTLs provides a platform for better understanding genetic effects on endometriosis risk and other endometrial-related pathologies. Study funding/competing interest(s) Funding for this work was provided by NHMRC Project Grants GNT1026033, GNT1049472, GNT1046880, GNT1050208, GNT1105321 and APP1083405. There are no competing interests.

A High-risk Haplotype for Premature Menopause in Childhood Cancer Survivors Exposed to Gonadotoxic Therapy

Background Childhood cancer survivors are at increased risk of therapy-related premature menopause (PM), with a cumulative incidence of 8.0%, but the contribution of genetic factors is unknown. Methods Genome-wide association analyses were conducted to identify single nucleotide polymorphisms (SNPs) associated with clinically diagnosed PM (menopause < 40 years) among 799 female survivors of childhood cancer participating in the St. Jude Lifetime Cohort Study (SJLIFE). Analyses were adjusted for cyclophosphamide equivalent dose of alkylating agents and ovarian radiotherapy (RT) dose (all P values two-sided). Replication was performed using self-reported PM in 1624 survivors participating in the Childhood Cancer Survivor Study (CCSS). Results PM was clinically diagnosed in 30 (3.8%) SJLIFE participants. Thirteen SNPs (70 kb region of chromosome 4q32.1) upstream of the Neuropeptide Receptor 2 gene (NPY2R) were associated with PM prevalence (minimum P = 3.3 × 10-7 for rs9999820, all P < 10-5). Being a homozygous carrier of a haplotype formed by four of the 13 SNPs (seen in one in seven in the general population but more than 50% of SJLIFE clinically diagnosed PM) was associated with markedly elevated PM prevalence among survivors exposed to ovarian RT (odds ratio [OR] = 25.89, 95% confidence interval [CI] = 6.18 to 138.31, P = 8.2 × 10-6); this finding was replicated in an independent second cohort of CCSS in spite of its use of self-reported PM (OR = 3.97, 95% CI = 1.67 to 9.41, P = .002). Evidence from bioinformatics data suggests that the haplotype alters the regulation of NPY2R transcription, possibly affecting PM risk through neuroendocrine pathways. Conclusions The haplotype captures the majority of clinically diagnosed PM cases and, with further validation, may have clinical application in identifying the highest-risk survivors for PM for possible intervention by cryopreservation.

Genome-wide search for higher order epistasis as modifiers of treatment effects on bone mineral density in childhood cancer survivors

Single-nucleotide polymorphisms (SNPs) contributing to interactions between regulatory elements that modulate gene transcription may explain some of the uncharacterized variation for complex traits. We explored this hypothesis among 856 adult survivors of pediatric cancer exposed to curative treatments that adversely affect bone mineral density (BMD). To restrict our search to interactions among SNPs in regulatory elements, our analysis considered 75523 SNPs mapped to putative promoter or enhancer regions. In anticipation that power to detect higher order epistasis would be low using an exhaustive search and a Bonferroni-corrected threshold for genome-wide significance (e.g., P < 5.6 × 10−14), a novel non-exhaustive statistical algorithm was implemented to detect chromosome-wide three-way regulatory interactions. We used a permutation-based evaluation statistic to identify candidate SNP interactions with stronger associations with BMD than expected. Of the six regulatory 3-SNP interactions identified as candidate interactions (P < 3.5 × 10−11) among cancer survivors exposed to treatments, five were replicated in an independent cohort of survivors (N = 1428) as modifiers of treatment effects on BMD (P < 0.05). Analyses with publicly available bioinformatics data revealed that SNPs contributing to replicated interactions were enriched for gene expressions (P = 3.6 × 10−4) and enhancer states (P < 0.05) in cells relevant for bone biology. For each replicated interaction, implicated SNPs were within or directly adjacent to 100-kb windows of genomic regions that plausibly physically interact in lymphoblastoid cells. Our study demonstrates the utility of a hypothesis-driven approach in revealing epistasis associated with complex traits.

Genetic Biomarkers for Endometriosis

GWAS studies identified seven genomic regions with robust evidence for genome-wide significant association with endometriosis risk. One important question that arises is whether these genetic markers can be used to predict risk of developing endometriosis for individual women. As with most complex diseases, the effect sizes for genetic markers linked to endometriosis risk are small with odds ratios less than 1.3. If we combine information from all seven markers, we explain only 1.85% of the total phenotypic variance on the liability scale (assuming a population prevalence of endometriosis of 8%) with no predictive power for individual risk.

Genetic risk factors for endometriosis

Advances in genetics and genomics are driving progress in understanding genetic risk factors for endometriosis. Genome-wide association scans (GWAS) in endometriosis have identified 11 genomic regions associated with increased risk of disease. Many of the regions contain interesting candidate genes, but the risk alleles may not always act through the obvious candidates. Functional evidence to identify the causal gene(s) will require multiple steps including better mapping precision, genetic studies on gene expression and epigenetic marks, chromatin looping and functional studies. Evidence from gene expression studies in endometrium and chromatin looping experiments implicate CDC42 on chromosome 1, CDKN2B-AS1 on chromosome 9 and VEZT on chromosome 12 as likely causal genes in these regions. Confirming the causal gene(s) in these and other regions will identify the important pathways increasing risk for endometriosis and identify novel targets for interventions to improve diagnosis and treatment.

Genome-wide haplotype association analysis of primary biliary cholangitis risk in Japanese

Primary biliary cholangitis (PBC) susceptibility loci have largely been discovered through single SNP association testing. In this study, we report genic haplotype patterns associated with PBC risk genome-wide in two Japanese cohorts. Among the 74 genic PBC risk haplotype candidates we detected with a novel methodological approach in a discovery cohort of 1,937 Japanese, nearly two-thirds were replicated (49 haplotypes, Bonferroni-corrected P < 6.8 × 10−4) in an independent Japanese cohort (N = 949). Along with corroborating known PBC-associated loci (TNFSF15, HLA-DRA), risk haplotypes may potentially model cis-interactions that regulate gene expression. For example, one replicated haplotype association (9q32–9q33.1, OR = 1.7, P = 3.0 × 10−21) consists of intergenic SNPs outside of the human leukocyte antigen (HLA) region that overlap regulatory histone mark peaks in liver and blood cells, and are significantly associated with TNFSF8 expression in whole blood. We also replicated a novel haplotype association involving non-HLA SNPs mapped to UMAD1 (7p21.3; OR = 15.2, P = 3.9 × 10−9) that overlap enhancer peaks in liver and memory Th cells. Our analysis demonstrates the utility of haplotype association analyses in discovering and characterizing PBC susceptibility loci.

AB028. Identifying the functional role of VEZT gene for endometriosis risk

Endometriosis is a common, chronic gynaecological disease characterised by pelvic pain and sub-fertility. It is a complex genetic disease affecting of 6-10% women of reproductive age and up to 50% of women with infertility. Evidence from genomewide association (GWA) imputed and meta-analysis studies in women with and without endometriosis identified association on chromosome 12q22, located close to the putative candidate gene VEZT. Fine-mapping this region in our 1,029 Australian cases and 958 controls confirmed strong evidence of association for non-coding variants in the 12q22 GWA region with the best imputed SNP rs4762347 in the 3’UTR of VEZT (P=0.036; 95% CI, 1.01-1.38; OR =1.18) and four top imputed SNPs in high linkage disequilibrium (LD) with rs4762347. Bioinformatic analysis for potential functional roles for these SNPs using ENCODE data show rs4762347 is highly conserved and alters a regulatory motif for the transcription factor Nkx3. We conducted an initial experiment to evaluate gene expression of VEZT in 36 samples of endometrial tissue from endometriosis cases and controls using high-throughput gene expression RT-qPCR on a microfluidic dynamic array IFC 48x48 chip (Fluidigm Biomark). The results showed that VEZT was expressed in endometrial samples from cases and controls. There was a trend for increased expression levels of VEZT mRNAs in endometriosis cases compared to controls as well as during the secretory phase compared to proliferative phase of menstrual cycle, although these differences were not statistically significant. The top imputed SNP rs4762347 was chosen to test for its allelic effects on VEZT expression. There was suggestive evidence for differences in VEZT mRNA expression observed in allelic groups of rs4762347. VEZT mRNA expression decreased for transcripts in carriers of minor C allele (CC and CT) of rs4762347 compared to non carriers (TT) (P<0.01), but the differences were no longer significant after corrections for multiple testing, hence further gene expression studies should be conducted in a larger sample size.