Protocol for the Endometriosis Research Queensland Study (ERQS): an integrated cohort study approach to improve diagnosis and stratify treatment.

INTRODUCTION: Endometriosis is a common gynaecological disease associated with pelvic pain and subfertility. There are no non-invasive diagnostic tests, medical management requires suppression of oestrogens and surgical removal is associated with risk. Endometriosis is a complex genetic disease with variants in at least 27 genetic regions associated with susceptibility. Previous research has implicated a variety of biological mechanisms in multiple cell types. Endometrial and endometriotic epithelial cells acquire somatic mutations at frequency higher than expected in normal tissue. Stromal cells have altered adhesive capacity and immune cells show altered cytotoxicity. Understanding the functional consequences of these genetic variants on each cell type requires the collection of patient symptoms, clinical and genetic data and disease-relevant tissue in an integrated program. METHODS AND ANALYSIS: The aims of this study are to collect tissue associated with endometriosis, chart the genetic architecture related to endometriosis in this tissue, isolate and propagate patient-specific cellular models, understand the functional consequence of these genetic variants and how they interact with environmental factors in pathogenesis and treatment response.We will collect patient information from online questionnaires prior to surgery and at 6 and 12 months postsurgery. Treating physicians will document detailed surgical data. During surgery, we will collect blood, peritoneal fluid, endometrium and endometriotic tissue. Tissue will be used to isolate and propagate in vitro models of individual cells. Genome wide genotyping and gene expression data will be generated. Somatic mutations will be identified via whole genome sequencing. ETHICS AND DISSEMINATION: The study has been approved and will be monitored by the Metro North Human Research Ethics committee (HREC) and research activities at the University of Queensland (UQ) will be overseen by the UQ HREC with annual reports submitted. Research results will be published in peer-reviewed journals and presented at conferences were appropriate. This study involves human participants and was approved by RBWH Human Research Ethics Committee; HREC/2019/QRBW/56763.The University of Queensland; 2017002744. Participants gave informed consent to participate in the study before taking part.

Association between DNA methylation variability and self-reported exposure to heavy metals.

Individuals encounter varying environmental exposures throughout their lifetimes. Some exposures such as smoking are readily observed and have high personal recall; others are more indirect or sporadic and might only be inferred from long occupational histories or lifestyles. We evaluated the utility of using lifetime-long self-reported exposures for identifying differential methylation in an amyotrophic lateral sclerosis cases-control cohort of 855 individuals. Individuals submitted paper-based surveys on exposure and occupational histories as well as whole blood samples. Genome-wide DNA methylation levels were quantified using the Illumina Infinium Human Methylation450 array. We analyzed 15 environmental exposures using the OSCA software linear and MOA models, where we regressed exposures individually by methylation adjusted for batch effects and disease status as well as predicted scores for age, sex, cell count, and smoking status. We also regressed on the first principal components on clustered environmental exposures to detect DNA methylation changes associated with a more generalised definition of environmental exposure. Five DNA methylation probes across three environmental exposures (cadmium, mercury and metalwork) were significantly associated using the MOA models and seven through the linear models, with one additionally across a principal component representing chemical exposures. Methylome-wide significance for four of these markers was driven by extreme hyper/hypo-methylation in small numbers of individuals. The results indicate the potential for using self-reported exposure histories in detecting DNA methylation changes in response to the environment, but also highlight the confounded nature of environmental exposure in cohort studies.

Analysis of DNA methylation associates the cystine-glutamate antiporter SLC7A11 with risk of Parkinson's disease.

An improved understanding of etiological mechanisms in Parkinson's disease (PD) is urgently needed because the number of affected individuals is projected to increase rapidly as populations age. We present results from a blood-based methylome-wide association study of PD involving meta-analysis of 229 K CpG probes in 1,132 cases and 999 controls from two independent cohorts. We identify two previously unreported epigenome-wide significant associations with PD, including cg06690548 on chromosome 4. We demonstrate that cg06690548 hypermethylation in PD is associated with down-regulation of the SLC7A11 gene and show this is consistent with an environmental exposure, as opposed to medications or genetic factors with effects on DNA methylation or gene expression. These findings are notable because SLC7A11 codes for a cysteine-glutamate anti-porter regulating levels of the antioxidant glutathione, and it is a known target of the environmental neurotoxin ß-methylamino-L-alanine (BMAA). Our study identifies the SLC7A11 gene as a plausible biological target in PD.

Genome-Wide Association Shows that Pigmentation Genes Play a Role in Skin Aging.

Loss of fine skin patterning is a sign of both aging and photoaging. Studies investigating the genetic contribution to skin patterning offer an opportunity to better understand a trait that influences both physical appearance and risk of keratinocyte skin cancer. We undertook a meta-analysis of genome-wide association studies of a measure of skin pattern (microtopography score) damage in 1,671 twin pairs and 1,745 singletons (N = 5,087) drawn from three independent cohorts. We identified that rs185146 near SLC45A2 is associated with a skin aging trait at genome-wide significance (P = 4.1 × 10-9); to our knowledge this is previously unreported. We also confirm previously identified loci, rs12203592 near IRF4 (P = 8.8 × 10-13) and rs4268748 near MC1R (P = 1.2 × 10-15). At all three loci we highlight putative functionally relevant SNPs. There are a number of red hair/low pigmentation alleles of MC1R; we found that together these MC1R alleles explained 4.1% of variance in skin pattern damage. We also show that skin aging and reported experience of sunburns was proportional to the degree of penetrance for red hair of alleles of MC1R. Our work has uncovered genetic contributions to skin aging and confirmed previous findings, showing that pigmentation is a critical determinant of skin aging.

Meta-analysis identifies five novel loci associated with endometriosis highlighting key genes involved in hormone metabolism.

Endometriosis is a heritable hormone-dependent gynecological disorder, associated with severe pelvic pain and reduced fertility; however, its molecular mechanisms remain largely unknown. Here we perform a meta-analysis of 11 genome-wide association case-control data sets, totalling 17,045 endometriosis cases and 191,596 controls. In addition to replicating previously reported loci, we identify five novel loci significantly associated with endometriosis risk (P<5 × 10-8), implicating genes involved in sex steroid hormone pathways (FN1, CCDC170, ESR1, SYNE1 and FSHB). Conditional analysis identified five secondary association signals, including two at the ESR1 locus, resulting in 19 independent single nucleotide polymorphisms (SNPs) robustly associated with endometriosis, which together explain up to 5.19% of variance in endometriosis. These results highlight novel variants in or near specific genes with important roles in sex steroid hormone signalling and function, and offer unique opportunities for more targeted functional research efforts.

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.

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.

Microsatellite stable colorectal cancers stratified by the BRAF V600E mutation show distinct patterns of chromosomal instability.

The BRAF (V600E) mutation in colorectal cancers that are microsatellite stable (MSS) confers a poor patient prognosis, whereas BRAF mutant microsatellite-unstable (MSI) colorectal cancers have an excellent prognosis. BRAF wild type cancers are typically MSS and display chromosomal instability (CIN). CIN has not been extensively studied on a genome-wide basis in relation to BRAF mutational status in colorectal cancer. BRAF mutant/MSS (BRAFmut/MSS) cancers (n = 33) and BRAF mutant/MSI (BRAFmut/MSI) cancers (n = 30) were compared for presence of copy number aberrations (CNAs) indicative of CIN, with BRAF wild type/MSS (BRAFwt/MSS) cancers (n = 18) using Illumina CytoSNP-12 arrays. BRAFmut/MSS and BRAFwt/MSS cancers showed comparable numbers of CNAs/cancer at 32.8 and 29.8 respectively. However, there were differences in patterns of CNA length between MSS cohorts, with BRAFmut/MSS cancers having significantly greater proportions of focal CNAs compared to BRAFwt/MSS cancers (p<0.0001); whereas whole chromosomal arm CNAs were more common in BRAFwt/MSS cancers (p<0.0001). This related to a reduced average CNA length in BRAFmut/MSS compared to BRAFwt/MSS cancers (20.7 Mb vs 33.4 Mb;p<0.0001); and a smaller average percent of CIN affected genomes in BRAFmut/MSS compared to BRAFwt/MSS cancers (23.9% vs 34.9% respectively). BRAFmut/MSI cancers were confirmed to have low CNA rates (5.4/cancer) and minimal CIN-affected genomes (average of 4.5%) compared to MSS cohorts (p<0.0001). BRAFmut/MSS cancers had more frequent deletion CNAs compared to BRAFwt/MSS cancers on 6p and 17q at loci not typically correlated with colorectal cancer, and greater amplification CNAs on 8q and 18q compared to BRAFwt/MSS cancers. These results indicate that comparable rates of CIN occur between MSS subgroups, however significant differences in their patterns of instability exist, with BRAFmut/MSS cancers showing a 'focal pattern' and BRAFwt/MSS cancers having a 'whole arm pattern' of CIN. This and the genomic loci more frequently affected in BRAFmut/MSS cancers provides further evidence of the biological distinctions of this important cancer subgroup.

Association of OPRD1 polymorphisms with heroin dependence in a large case-control series.

Genes encoding the opioid receptors (OPRM1, OPRD1 and OPRK1) are obvious candidates for involvement in risk for heroin dependence. Prior association studies commonly had samples of modest size, included limited single nucleotide polymorphism (SNP) coverage of these genes and yielded inconsistent results. Participants for the current investigation included 1459 heroin-dependent cases ascertained from maintenance clinics in New South Wales, Australia, 1495 unrelated individuals selected from an Australian sample of twins and siblings as not meeting DSM-IV criteria for lifetime alcohol or illicit drug dependence (non-dependent controls) and 531 controls ascertained from economically disadvantaged neighborhoods in proximity to the maintenance clinics. A total of 136 OPRM1, OPRD1 and OPRK1 SNPs were genotyped in this sample. After controlling for admixture with principal components analysis, our comparison of cases to non-dependent controls found four OPRD1 SNPs in fairly high linkage disequilibrium for which adjusted P values remained significant (e.g. rs2236857; OR 1.25; P=2.95×10(-4) ) replicating a previously reported association. A post hoc analysis revealed that the two SNP (rs2236857 and rs581111) GA haplotype in OPRD1 is associated with greater risk (OR 1.68; P=1.41×10(-5) ). No OPRM1 or OPRK1 SNPs reached more than nominal significance. Comparisons of cases to neighborhood controls reached only nominal significance. Our results replicate a prior report providing strong evidence implicating OPRD1 SNPs and, in particular, the two SNP (rs2236857 and rs581111) GA haplotype in liability for heroin dependence. Support was not found for similar association involving either OPRM1 or OPRK1 SNPs.

ANKK1, TTC12, and NCAM1 polymorphisms and heroin dependence: importance of considering drug exposure.

CONTEXT: The genetic contribution to liability for opioid dependence is well established; identification of the responsible genes has proved challenging. OBJECTIVE: To examine association of 1430 candidate gene single-nucleotide polymorphisms (SNPs) with heroin dependence, reporting here only the 71 SNPs in the chromosome 11 gene cluster (NCAM1, TTC12, ANKK1, DRD2) that include the strongest observed associations. DESIGN: Case-control genetic association study that included 2 control groups (lacking an established optimal control group). SETTING: Semistructured psychiatric interviews. PARTICIPANTS: A total of 1459 Australian cases ascertained from opioid replacement therapy clinics, 531 neighborhood controls ascertained from economically disadvantaged areas near opioid replacement therapy clinics, and 1495 unrelated Australian Twin Registry controls not dependent on alcohol or illicit drugs selected from a twin and family sample. MAIN OUTCOME MEASURE: Lifetime heroin dependence. RESULTS: Comparison of cases with Australian Twin Registry controls found minimal evidence of association for all chromosome 11 cluster SNPs (P ≥ .01); a similar comparison with neighborhood controls revealed greater differences (P ≥ 1.8 × 10(-4)). Comparing cases (n = 1459) with the subgroup of neighborhood controls not dependent on illicit drugs (n = 340), 3 SNPs were significantly associated (correcting for multiple testing): ANKK1 SNP rs877138 (most strongly associated; odds ratio = 1.59; 95% CI, 1.32-1.92; P = 9.7 × 10(-7)), ANKK1 SNP rs4938013, and TTC12 SNP rs7130431. A similar pattern of association was observed when comparing illicit drug-dependent (n = 191) and nondependent (n = 340) neighborhood controls, suggesting that liability likely extends to nonopioid illicit drug dependence. Aggregate heroin dependence risk associated with 2 SNPs, rs877138 and rs4492854 (located in NCAM1), varied more than 4-fold (P = 2.7 × 10(-9) for the risk-associated linear trend). CONCLUSIONS: Our results provide further evidence of association for chromosome 11 gene cluster SNPs with substance dependence, including extension of liability to illicit drug dependence. Our findings highlight the necessity of considering drug exposure history when selecting control groups for genetic investigations of illicit drug dependence.

Examining the association of NRXN3 SNPs with borderline personality disorder phenotypes in heroin dependent cases and socio-economically disadvantaged controls.

BACKGROUND: Borderline personality disorder (BPD) and substance use disorders frequently co-occur; their dual presence predicts poor prognosis. The genetic underpinnings of BPD have not been well-characterized and could offer insight into comorbidity. The current report focuses on the association of neurexin 3 (NRXN3) single nucleotide polymorphisms (SNPs) with BPD symptoms in heroin dependent cases and controls. METHODS: The sample of the Comorbidity and Trauma Study, a genetic association study of heroin dependence, consists of Australian heroin dependent cases ascertained from opioid replacement therapy clinics and controls ascertained in nearby economically disadvantaged neighborhoods. The assessment included a screening instrument for BPD, used previously in Australian population surveys. Genotypic and BPD phenotypic data were available for 1439 cases and 507 controls. We examined the association of 1430 candidate gene SNPs with BPD phenotypes. RESULTS: One or more NRXN3 SNPs were nominally associated with all BPD phenotypes; however, none met the conservative significance threshold we employed to correct for multiple testing. The most strongly associated SNPs included rs10144398 with identity disturbance (p=4.9×10(-5)) and rs10151731 with affective instability (p=8.8×10(-5)). The strongest association with screening positive for BPD was found for the NRXN3 SNP, rs10083466 (p=.0013). Neither the correlation of BPD phenotypes nor the linkage disequilibrium relationships of the SNPs account for the number of observed associations involving NRXN3 SNPs. CONCLUSIONS: Our findings provide intriguing preliminary evidence for the association of NRXN3 with BPD phenotypes. The strongest associations were found for traits (i.e., affective instability; identity disturbance) also observed with other disorders.

Sequence variants in three loci influence monocyte counts and erythrocyte volume.

Blood cells participate in vital physiological processes, and their numbers are tightly regulated so that homeostasis is maintained. Disruption of key regulatory mechanisms underlies many blood-related Mendelian diseases but also contributes to more common disorders, including atherosclerosis. We searched for quantitative trait loci (QTL) for hematology traits through a whole-genome association study, because these could provide new insights into both hemopoeitic and disease mechanisms. We tested 1.8 million variants for association with 13 hematology traits measured in 6015 individuals from the Australian and Dutch populations. These traits included hemoglobin composition, platelet counts, and red blood cell and white blood cell indices. We identified three regions of strong association that, to our knowledge, have not been previously reported in the literature. The first was located in an intergenic region of chromosome 9q31 near LPAR1, explaining 1.5% of the variation in monocyte counts (best SNP rs7023923, p=8.9x10(-14)). The second locus was located on chromosome 6p21 and associated with mean cell erythrocyte volume (rs12661667, p=1.2x10(-9), 0.7% variance explained) in a region that spanned five genes, including CCND3, a member of the D-cyclin gene family that is involved in hematopoietic stem cell expansion. The third region was also associated with erythrocyte volume and was located in an intergenic region on chromosome 6q24 (rs592423, p=5.3x10(-9), 0.6% variance explained). All three loci replicated in an independent panel of 1543 individuals (p values=0.001, 9.9x10(-5), and 7x10(-5), respectively). The identification of these QTL provides new opportunities for furthering our understanding of the mechanisms regulating hemopoietic cell fate.

Common variants in TMPRSS6 are associated with iron status and erythrocyte volume.

We report a genome-wide association study to iron status. We identify an association of SNPs in TPMRSS6 to serum iron (rs855791, combined P = 1.5 x 10(-20)), transferrin saturation (combined P = 2.2 x 10(-23)) and erythrocyte mean cell volume (MCV, combined P = 1.1 x 10(-10)). We also find suggestive evidence of association with blood hemoglobin levels (combined P = 5.3 x 10(-7)). These findings demonstrate the involvement of TMPRSS6 in control of iron homeostasis and in normal erythropoiesis.