Bayesian refinement of association signals for 14 loci in 3 common diseases

To further investigate susceptibility loci identified by genome-wide association studies, we genotyped 5,500 SNPs across 14 associated regions in 8,000 samples from a control group and 3 diseases: type 2 diabetes (T2D), coronary artery disease (CAD) and Graves' disease. We defined, using Bayes theorem, credible sets of SNPs that were 95% likely, based on posterior probability, to contain the causal disease-associated SNPs. In 3 of the 14 regions, TCF7L2 (T2D), CTLA4 (Graves' disease) and CDKN2A-CDKN2B (T2D), much of the posterior probability rested on a single SNP, and, in 4 other regions (CDKN2A-CDKN2B (CAD) and CDKAL1, FTO and HHEX (T2D)), the 95% sets were small, thereby excluding most SNPs as potentially causal. Very few SNPs in our credible sets had annotated functions, illustrating the limitations in understanding the mechanisms underlying susceptibility to common diseases. Our results also show the value of more detailed mapping to target sequences for functional studies. © 2012 Nature America, Inc. All rights reserved.

Genome-wide association study using extreme truncate selection identifies novel genes affecting bone mineral density and fracture risk

Osteoporotic fracture is a major cause of morbidity and mortality worldwide. Low bone mineral density (BMD) is a major predisposing factor to fracture and is known to be highly heritable. Site-, gender-, and age-specific genetic effects on BMD are thought to be significant, but have largely not been considered in the design of genome-wide association studies (GWAS) of BMD to date. We report here a GWAS using a novel study design focusing on women of a specific age (postmenopausal women, age 55-85 years), with either extreme high or low hip BMD (age- and gender-adjusted BMD z-scores of +1.5 to +4.0, n = 1055, or -4.0 to -1.5, n = 900), with replication in cohorts of women drawn from the general population (n = 20,898). The study replicates 21 of 26 known BMD-associated genes. Additionally, we report suggestive association of a further six new genetic associations in or around the genes CLCN7, GALNT3, IBSP, LTBP3, RSPO3, and SOX4, with replication in two independent datasets. A novel mouse model with a loss-of-function mutation in GALNT3 is also reported, which has high bone mass, supporting the involvement of this gene in BMD determination. In addition to identifying further genes associated with BMD, this study confirms the efficiency of extreme-truncate selection designs for quantitative trait association studies. © 2011 Duncan et al.

Association of STAT3 and TNFRSF1A with ankylosing spondylitis in Han Chinese

Objectives: Recent association studies by the Australo-Anglo-American Spondyloarthritis Consortium (TASC) in Caucasian European populations from Australia, North America and the UK have identified a number of genes as being associated with ankylosing spondylitis (AS). A candidate gene study in a Han Chinese population was performed based on these findings to identify associated genes in this population. Methods: A case-control study was performed in a Han Chinese population of patients with AS (n=775) and controls (n=1587) from Shanghai and Nanjing. All patients met the modified New York criteria for AS. The cases and controls were genotyped for 115 single nucleotide polymorphisms (SNPs) tagging IL23R, ERAP1, STAT3, JAK2, TNFRSF1A and TRADD, as well as other confirmation SNPs from the TASC study, using the Sequenom iPlex and the ABI OpenArray platforms. Statistical analysis of genotyped SNPs was performed using the Cochran - Armitage test for trend and meta-analysis was performed using METAL. SNPs in AS-associated genes in this study were then imputed using MaCH, and association with AS tested by logistic regression. Results: SNPs in TNFRSF1A (rs4149577, p=8.2×10-4), STAT3 (rs2293152, p=0.0015; rs1053005, p=0.017) and ERAP1 (rs27038, p=0.0091; rs27037, p=0.0092) were significantly associated with AS in Han Chinese. Association was also observed between AS and the intergenic region 2p15 (rs10865331, p=0.023). The lack of association between AS and IL23R in Han Chinese was confirmed (all SNPs p>0.1). Conclusions: The study results demonstrate for the first time that genetic polymorphisms in STAT3, TNFRSF1A and 2p15 are associated with AS in Han Chinese, suggesting common pathogenic mechanisms for the disease in Chinese and Caucasian European populations. Furthermore, previous findings demonstrating that ERAP1, but not IL23R, is associated with AS in Chinese patients were confirmed.

Progress in the genetics of ankylosing spondylitis

Ankylosing spondylitis (AS) is a common, highly heritable, inflammatory arthropathy. In addition to being strongly associated with HLA-B27, a further 13 genes have been robustly associated with the disease. These genes highlight the involvement of the IL-23 pathway in disease pathogenesis, and indicate overlaps between the pathogenesis of AS, and of inflammatory bowel disease. Genetic associations in B27-positive and -negative disease are similar, with the main exception of association with ERAP1, which is restricted in association to B27-positive cases. This restriction, and the known function of ERAP1 in peptide trimming prior to HLA Class I presentation, indicates that HLA-B27 is likely to operate in AS by a mechanism involving aberrant peptide handling. These advances point to several potential novel therapeutic approaches in AS.

Interaction between ERAP1 and HLA-B27 in ankylosing spondylitis implicates peptide handling in the mechanism for HLA-B27 in disease susceptibility

Ankylosing spondylitis is a common form of inflammatory arthritis predominantly affecting the spine and pelvis that occurs in approximately 5 out of 1,000 adults of European descent. Here we report the identification of three variants in the RUNX3, LTBR-TNFRSF1A and IL12B regions convincingly associated with ankylosing spondylitis (P < 5 × 10-8 in the combined discovery and replication datasets) and a further four loci at PTGER4, TBKBP1, ANTXR2 and CARD9 that show strong association across all our datasets (P < 5 × 10-6 overall, with support in each of the three datasets studied). We also show that polymorphisms of ERAP1, which encodes an endoplasmic reticulum aminopeptidase involved in peptide trimming before HLA class I presentation, only affect ankylosing spondylitis risk in HLA-B27-positive individuals. These findings provide strong evidence that HLA-B27 operates in ankylosing spondylitis through a mechanism involving aberrant processing of antigenic peptides.

Genetic determinants of bone density and fracture risk: State of the art and future directions

Context: Osteoporosis is a common, highly heritable condition that causes substantial morbidity and mortality, the etiopathogenesis of which is poorly understood. Genetic studies are making increasingly rapid progress in identifying the genes involved. Evidence Acquisition and Synthesis: In this review, we will summarize the current understanding of the genetics of osteoporosis based on publications from PubMed from the year 1987 onward. Conclusions: Most genes involved in osteoporosis identified to date encode components of known pathways involved in bone synthesis or resorption, but as the field progresses, new pathways are being identified. Only a small proportion of the total genetic variation involved in osteoporosis has been identified, and new approaches will be required to identify most of the remaining genes.

Suggestive linkage of the parathyroid receptor type 1 to osteoporosis

We have investigated the role of 23 candidate genes in the control of bone mineral density (BMD) by linkage studies in families of probands with osteoporosis (lumbar spine [LS] or femoral neck [FN] BMD T score < -2.5) and low BMD relative to an age- and gender-matched cohort (Z score < -2.0). One hundred and fifteen probands (35 male, 80 female) and 499 of their first- or second-degree relatives (223 males and 276 females) were recruited for the study. BMD was measured at the LS and FN using dual-energy X-ray absorptiometry and expressed as age- and gender-matched Z scores corrected for body mass index. The candidate genes studied were the androgen receptor, type I collagen A1 (COLIA1), COLIA2, COLIIA1, vitamin D receptor (VDR), colony-stimulating factor 1, calcium-sensing receptor, epidermal growth factor (EGF), estrogen receptor 1 (ESR1), fibrillin type 1, insulin-like growth factor 1, interleukin-1 alpha (IL-1α), interleukin-4 (IL-4), interleukin-6 (IL-6), interleukin-11 (IL-11), osteopontin, parathyroid hormone (PTH), PTH-related peptide, PTH receptor type 1 (PTHR1), transforming growth factor-beta 1, and tumor necrosis factors alpha and beta. Sixty-four microsatellites lying close to or within these genes were investigated for linkage with BMD. Using the program MapMaker/Sibs there was suggestive evidence of linkage between BMD and PTHR1 (maximum LOD score obtained [MLS] 2.7-3.5). Moderate evidence of linkage was also observed with EGF (MLS 1.8), COLIA1 (MLS 1.7), COLIIA1/VDR (MLS 1.7), ESR1 (MLS 1.4), IL-1α (MLS 1.4), IL-4 (MLS 1.2), and IL-6 (MLS 1.2). Variance components analysis using the program ACT, correcting for proband-wise ascertainment, also showed evidence of linkage (p ≤0.05) at markers close to or within the candidate genes IL- 1α, PTHR1, IL-6, and COLIIA1/VDR. Further studies will be required to confirm these findings, to refine the location of gene responsible for the observed linkage, and to screen the candidate genes targeted at these loci for mutations.

Conformational Changes Triggered by Mg2+ Mediate Transactivator Function

Transactivator protein C of bacteriophage mu is essential for the transition from middle to late gene expression during the phage life cycle. The unusual, multistep activation of mom promoter (Pmom) by C protein involves activator-mediated promoter unwinding to recruit RNA polymerase and subsequent enhanced promoter clearance of the enzyme. To achieve this, C binds its site overlapping the -35 region of the mom promoter with a very high affinity, in Mg2+-dependent fashion. Mg2+-mediated conformational transition in C is necessary for its DNA binding and transactivation. We have determined the residues in C which coordinate Mg2+, to induce allosteric transition in the protein, required for the specific interaction with DNA. Residues E26 and D40 in the putative metal binding motif (E26X10D37X2D40) present toward the N-terminus of the protein are found to be important for Mg2+ ion binding. Mutations in these residues lead to altered Mg2+-induced conformation, compromised DNA binding, and reduced levels of transcription activation. Although Mg2+ is widely used in various DNA transaction reactions, this report provides the first insights on the importance of the metal ion-induced allosteric transitions in regulating transcription factor function.

A linear programming approach for estimating the structure of a sparse linear genetic network from transcript profiling data

Background: A genetic network can be represented as a directed graph in which a node corresponds to a gene and a directed edge specifies the direction of influence of one gene on another. The reconstruction of such networks from transcript profiling data remains an important yet challenging endeavor. A transcript profile specifies the abundances of many genes in a biological sample of interest. Prevailing strategies for learning the structure of a genetic network from high-dimensional transcript profiling data assume sparsity and linearity. Many methods consider relatively small directed graphs, inferring graphs with up to a few hundred nodes. This work examines large undirected graphs representations of genetic networks, graphs with many thousands of nodes where an undirected edge between two nodes does not indicate the direction of influence, and the problem of estimating the structure of such a sparse linear genetic network (SLGN) from transcript profiling data. Results: The structure learning task is cast as a sparse linear regression problem which is then posed as a LASSO (l1-constrained fitting) problem and solved finally by formulating a Linear Program (LP). A bound on the Generalization Error of this approach is given in terms of the Leave-One-Out Error. The accuracy and utility of LP-SLGNs is assessed quantitatively and qualitatively using simulated and real data. The Dialogue for Reverse Engineering Assessments and Methods (DREAM) initiative provides gold standard data sets and evaluation metrics that enable and facilitate the comparison of algorithms for deducing the structure of networks. The structures of LP-SLGNs estimated from the INSILICO1, INSILICO2 and INSILICO3 simulated DREAM2 data sets are comparable to those proposed by the first and/or second ranked teams in the DREAM2 competition. The structures of LP-SLGNs estimated from two published Saccharomyces cerevisae cell cycle transcript profiling data sets capture known regulatory associations. In each S. cerevisiae LP-SLGN, the number of nodes with a particular degree follows an approximate power law suggesting that its degree distributions is similar to that observed in real-world networks. Inspection of these LP-SLGNs suggests biological hypotheses amenable to experimental verification. Conclusion: A statistically robust and computationally efficient LP-based method for estimating the topology of a large sparse undirected graph from high-dimensional data yields representations of genetic networks that are biologically plausible and useful abstractions of the structures of real genetic networks. Analysis of the statistical and topological properties of learned LP-SLGNs may have practical value; for example, genes with high random walk betweenness, a measure of the centrality of a node in a graph, are good candidates for intervention studies and hence integrated computational – experimental investigations designed to infer more realistic and sophisticated probabilistic directed graphical model representations of genetic networks. The LP-based solutions of the sparse linear regression problem described here may provide a method for learning the structure of transcription factor networks from transcript profiling and transcription factor binding motif data.

Growth defect in Grg5 null mice is associated with reduced Ihh signaling in growth plates

Gene-targeted disruption of Grg5, a mouse homologue of Drosophila groucho (gro), results in postnatal growth retardation in mice. The growth defect, most striking in approximately half of the Grg5 null mice, occurs during the first 4-5 weeks of age, but most mice recover retarded growth later. We used the nonlinear mixed-effects model to fit the growth data of wild-type, heterozygous, and Grg5 null mice. On the basis of preliminary evidence suggesting an interaction between Grg5 and the transcription factor Cbfa1/Runx2, critical for skeletal development, we further investigated the skeleton in the mice. A long bone growth plate defect was identified, which included shorter zones of proliferative and hypertrophic chondrocytes and decreased trabecular bone formation. This decreased trabecular bone formation is likely caused by a reduced recruitment of osteoblasts into the growth plate region of Grg5 null mice. Like the growth defect, the growth plate and trabecular bone abnormality improved as the mice grew older. The growth plate defect was associated with reduced Indian hedgehog expression and signaling. We suggest that Grg5, a transcriptional coregulator, modulates the activities of transcription factors, such as Cbfa1/Runx2 in vivo to affect Ihh expression and the function of long bone growth plates.

Cytochrome P450 aromatase in grey mullet: cDNA and promoter isolation; brain, pituitary and ovarian expression during puberty.

In a study towards elucidating the role of aromatases during puberty in female grey mullet, the cDNAs of the brain (muCyp19b) and ovarian (muCyp19a) aromatase were isolated by RT-PCR and their relative expression levels were determined by quantitative real-time RT-PCR. The muCyp19a ORF of 1515 bp encoded 505 predicted amino acid residues, while that of muCyp19b was 1485 bp and encoded 495 predicted amino acid residues. The expression level of muCyp19b significantly increased in the brain as puberty advanced; however, its expression level in the pituitary increased only slightly with pubertal development. In the ovary, the muCyp19a expression level markedly increased as puberty progressed. The promoter regions of the two genes were also isolated and their functionality evaluated in vitro using luciferase as the reporter gene. The muCyp19a promoter sequence (650 bp) contained a consensus TATA box and putative transcription factor binding sites, including two half EREs, an SF-1, an AhR/Arnt, a PR and two GATA-3s. The muCyp19b promoter sequence (2500 bp) showed consensus TATA and CCAAT boxes and putative transcription binding sites, namely: a PR, an ERE, a half ERE, a SP-1, two GATA-binding factor, one half GATA-1, two C/EBPs, a GRE, a NFkappaB, three STATs, a PPAR/RXR, an Ahr/Arnt and a CRE. Basal activity of serially deleted promoter constructs transiently transfected into COS-7, [alpha]T3 and TE671 cells demonstrated the enhancing and silencing roles of the putative transcription factor binding sites. Quinpirole, a dopamine agonist, significantly reduced the promoter activity of muCyp19b in TE671. The results suggest tissue-specific regulation of the muCyp19 genes and a putative alternative promoter for muCyp19b.

Epithelial-mesenchymal transition in oral squamous carcinoma cells

In epithelial-mesenchymal transition (EMT), epithelial cells acquire traits typical for mesenchymal cells, dissociate their cell-cell junctions and gain the ability to migrate. EMT is essential during embryogenesis, but may also mediate cancer progression. Basement membranes are sheets of extracellular matrix that support epithelial cells. They have a major role in maintaining the epithelial phenotype and, in cancer, preventing cell migration, invasion and metastasis. Laminins are the main components of basement membranes and may actively contribute to malignancy. We first evaluated the differences between cell lines obtained from oral squamous cell carcinoma and its recurrence. As the results indicated a change from epithelial to fibroblastoid morphology, E-cadherin to N-cadherin switch, and change in expression of cytokeratins to vimentin intermediate filaments, we concluded that these cells had undergone EMT. We further induced EMT in primary tumour cells to gain knowledge of the effects of transcription factor Snail in this cell model. The E-cadherin repressors responsible for the EMT in these cells were ZEB-1, ZEB-2 and Snail, and ectopic expression of Snail was able to augment the levels of ZEB-1 and ZEB-2. We produced and characterized two monoclonal antibodies that specifically recognized Snail in cell lines and patient samples. By immunohistochemistry, Snail protein was found in mesenchymal tissues during mouse embryonal development, in fibroblastoid cells of healing skin wounds and in fibromatosis and sarcoma specimens. Furthermore, Snail localized to the stroma and borders of tumour cell islands in colon adenocarcinoma, and in laryngeal and cervical squamous cell carcinomas. Immunofluorescence labellings, immunoprecipitations and Northern and Western blots showed that EMT induced a progressive downregulation of laminin-332 and laminin-511 and, on the other hand, an induction of mesenchymal laminin-411. Chromatin immunoprecipitation revealed that Snail could directly bind upstream to the transcription start sites of both laminin α5 and α4 chain genes, thus regulating their expression. The levels of integrin α6β4, a receptor for laminin-332, as well as the hemidesmosomal complex proteins HD1/plectin and BP180 were downregulated in EMT-experienced cells. The expression of Lutheran glycoprotein, a specific receptor for laminin-511, was diminished, whereas the levels of integrins α6β1 and α1β1 and integrin-linked kinase were increased. In quantitative cell adhesion assays, the cells adhered potently to laminin-511 and fibronectin, but only marginally to laminin-411. Western blots and immunoprecipitations indicated that laminin-411 bound to fibronectin and could compromise cell adhesion to fibronectin in a dose-dependent manner. EMT induced a highly migratory and invasive tendency in oral squamous carcinoma cells. Actin-based adhesion and invasion structures, podosomes and invadopodia, were detected in the basal cell membranes of primary tumour and spontaneously transformed cancer cells, respectively. Immunofluorescence labellings showed marked differences in their morphology, as podosomes organized a ring structure with HD1/plectin, αII-spectrin, talin, focal adhesion kinase and pacsin 2 around the core filled with actin, cortactin, vinculin and filamin A. Invadopodia had no division between ring and core and failed to organize the ring proteins, but instead assembled tail-like, narrow actin cables that showed a talin-tensin switch. Time-lapse live-cell imaging indicated that both podosomes and invadopodia were long-lived entities, but the tails of invadopodia vigorously propelled in the cytoplasm and were occasionally released from the cell membrane. Invadopodia could also be externalized outside the cytoplasm, where they still retained the ability to degrade matrix. In 3D confocal imaging combined with in situ gelatin zymography, the podosomes of primary tumour cells were large, cylindrical structures that increased in time, whereas the invadopodia in EMT-driven cells were smaller, but more numerous and degraded the underlying matrix in significantly larger amounts. Fluorescence recovery after photobleaching revealed that the substructures of podosomes were replenished more rapidly with new molecules than those of invadopodia. Overall, our results indicate that EMT has a major effect on the transcription and synthesis of both intra- and extracellular proteins, including laminins and their receptors, and on the structure and dynamics of oral squamous carcinoma cells.

Candidate gene studies on cardiovascular traits

Cardiovascular diseases (CVD) are a major cause of death and disability in Western countries and a growing health problem in the developing world. The genetic component of both coronary heart disease (CHD) and ischemic stroke events has been established in twin studies, and the traits predisposing to CVD, such as hypertension, dyslipidemias, obesity, diabetes, and smoking behavior, are all partly hereditary. Better understanding of the pathophysiology of CVD-related traits could help to target disease prevention and clinical treatment to individuals at an especially high disease risk and provide novel pharmaceutical interventions. This thesis aimed to clarify the genetic background of CVD at a population level using large Nordic population cohorts and a candidate gene approach. The first study concentrated on the allelic diversity of the thrombomodulin (THBD) gene in two Finnish cohorts, FINRISK-92 and FINRISK-97. The results from this study implied that THBD variants do not substantially contribute to CVD risk. In the second study, three other candidate genes were added to the analyses. The study investigated the epistatic effects of coagulation factor V (F5), intercellular adhesion molecule -1 (ICAM1), protein C (PROC), and THBD in the same FINRISK cohorts. The results were encouraging; we were able to identify several single SNPs and SNP combinations associating with CVD and mortality. Interestingly, THBD variants appeared in the associating SNP combinations despite the negative results from Study I, suggesting that THBD contributes to CVD through gene-gene interactions. In the third study, upstream transcription factor -1 (USF1) was analyzed in a cohort of Swedish men. USF1 was associated with metabolic syndrome, characterized by accumulation of different CVD risk factors. A putative protective and a putative risk variant were identified. A direct association with CVD was not observed. The longitudinal nature of the study also clarified the effect of USF1 variants on CVD risk factors followed in four examinations throughout adulthood. The three studies provided valuable information on the study of complex traits, highlighting the use of large study samples, the importance of replication, and the full coverage of the major allelic variants of the target genes to assure reliable findings. Although the genetic basis of coronary heart disease and ischemic stroke remains unknown, single genetic findings may facilitate the recognition of high-risk subgroups.

Molecular Mechanisms of Androgen Receptor Interactions

The androgen receptor (AR) mediates the effects of the male sex-steroid hormones (androgens), testosterone and 5?-dihydrotestosterone. Androgens are critical in the development and maintenance of male sexual characteristics. AR is a member of the steroid receptor ligand-inducible transcription factor family. The steroid receptor family is a subgroup of the nuclear receptor superfamily that also includes receptors for the active forms of vitamin A, vitamin D3, and thyroid hormones. Like all nuclear receptors, AR has a conserved modular structure consisting of a non-conserved amino-terminal domain (NTD), containing the intrinsic activation function 1, a highly conserved DNA-binding domain, and a conserved ligand-binding domain (LBD) that harbors the activation function 2. Each of these domains plays an important role in receptor function and signaling, either via intra- and inter-receptor interactions, interactions with specific DNA sequences, termed hormone response elements, or via functional interactions with domain-specific proteins, termed coregulators (coactivators and corepressors). Upon binding androgens, AR acquires a new conformational state, translocates to the nucleus, binds to androgen response elements, homodimerizes and recruits sequence-specific coregulatory factors and the basal transcription machinery. This set of events is required to activate gene transcription (expression). Gene transcription is a strictly modulated process that governs cell growth, cell homeostasis, cell function and cell death. Disruptions of AR transcriptional activity caused by receptor mutations and/or altered coregulator interactions are linked to a wide spectrum of androgen insensitivity syndromes, and to the pathogenesis of prostate cancer (CaP). The treatment of CaP usually involves androgen depletion therapy (ADT). ADT achieves significant clinical responses during the early stages of the disease. However, under the selective pressure of androgen withdrawal, androgen-dependent CaP can progress to an androgen-independent CaP. Androgen-independent CaP is invariably a more aggressive and untreatable form of the disease. Advancing our understanding of the molecular mechanisms behind the switch in androgen-dependency would improve our success of treating CaP and other AR related illnesses. This study evaluates how clinically identified AR mutations affect the receptor s transcriptional activity. We reveal that a potential molecular abnormality in androgen insensitivity syndrome and CaP patients is caused by disruptions of the important intra-receptor NTD/LBD interaction. We demonstrate that the same AR LBD mutations can also disrupt the recruitment of the p160 coactivator protein GRIP1. Our investigations reveal that 30% of patients with advanced, untreated local CaP have somatic mutations that may lead to increases in AR activity. We report that somatic mutations that activate AR may lead to early relapse in ADT. Our results demonstrate that the types of ADT a CaP patient receives may cause a clustering of mutations to a particular region of the receptor. Furthermore, the mutations that arise before and during ADT do not always result in a receptor that is more active, indicating that coregulator interactions play a pivotal role in the progression of androgen-independent CaP. To improve CaP therapy, it is necessary to identify critical coregulators of AR. We screened a HeLa cell cDNA library and identified small carboxyl-terminal domain phosphatase 2 (SCP2). SCP2 is a protein phosphatase that directly interacts with the AR NTD and represses AR activity. We demonstrated that reducing the endogenous cellular levels of SCP2 causes more AR to load on to the prostate specific antigen (PSA) gene promoter and enhancer regions. Additionally, under the same conditions, more RNA polymerase II was recruited to the PSA promoter region and overall there was an increase in androgen-dependent transcription of the PSA gene, revealing that SCP2 could play a role in the pathogenesis of CaP.

Functional evaluation of genetic variants associated with endometriosis near GREB1

STUDY QUESTION: Do DNA variants in the growth regulation by estrogen in breast cancer 1 (GREB1) region regulate endometrial GREB1 expression and increase the risk of developing endometriosis in women? SUMMARY ANSWER: We identified new single nucleotide polymorphisms (SNPs) with strong association with endometriosis at the GREB1 locus although we did not detect altered GREB1 expression in endometriosis patients with defined genotypes. WHAT IS ALREADY KNOWN: Genome-wide association studies have identified the GREB1 region on chromosome 2p25.1 for increasing endometriosis risk. The differential expression of GREB1 has also been reported by others in association with endometriosis disease phenotype. STUDY DESIGN, SIZE, DURATION: Fine mapping studies comprehensively evaluated SNPs within the GREB1 region in a large-scale data set (>2500 cases and >4000 controls). Publicly available bioinformatics tools were employed to functionally annotate SNPs showing the strongest association signal with endometriosis risk. Endometrial GREB1 mRNA and protein expression was studied with respect to phases of the menstrual cycle (n = 2-45 per cycle stage) and expression quantitative trait loci (eQTL) analysis for significant SNPs were undertaken for GREB1 [mRNA (n = 94) and protein (n = 44) in endometrium]. PARTICIPANTS/MATERIALS, SETTING, METHODS: Participants in this study are females who provided blood and/or endometrial tissue samples in a hospital setting. The key SNPs were genotyped using Sequenom MassARRAY. The functional roles and regulatory annotations for identified SNPs are predicted by various publicly available bioinformatics tools. Endometrial GREB1 expression work employed qRT-PCR, western blotting and immunohistochemistry studies. MAIN RESULTS AND THE ROLE OF CHANCE: Fine mapping results identified a number of SNPs showing stronger association (0.004 < P < 0.032) with endometriosis risk than the original GWAS SNP (rs13394619) (P = 0.034). Some of these SNPs were predicted to have functional roles, for example, interaction with transcription factor motifs. The haplotype (a combination of alleles) formed by the risk alleles from two common SNPs showed significant association (P = 0.026) with endometriosis and epistasis analysis showed no evidence for interaction between the two SNPs, suggesting an additive effect of SNPs on endometriosis risk. In normal human endometrium, GREB1 protein expression was altered depending on the cycle stage (significantly different in late proliferative versus late secretory, P < 0.05) and cell type (glandular epithelium, not stromal cells). However, GREB1 expression in endometriosis cases versus controls and eQTL analyses did not reveal any significant changes. LIMITATIONS, REASONS FOR CAUTION: In silico prediction tools are generally based on cell lines different to our tissue and disease of interest. Functional annotations drawn from these analyses should be considered with this limitation in mind. We identified cell-specific and hormone-specific changes in GREB1 protein expression. The lack of a significant difference observed following our GREB1 expression studies may be the result of moderate power on mixed cell populations in the endometrial tissue samples. WIDER IMPLICATIONS OF THE FINDINGS: This study further implicates the GREB1 region on chromosome 2p25.1 and the GREB1 gene with involvement in endometriosis risk. More detailed functional studies are required to determine the role of the novel GREB1 transcripts in endometriosis pathophysiology. STUDY FUNDING/COMPETING INTERESTS: Funding for this work was provided by NHMRC Project Grants APP1012245, APP1026033, APP1049472 and APP1046880. There are no competing interests.