Qtlizer: comprehensive QTL annotation of GWAS results

A Multilayered Post-Genome-Wide Association Study Analysis Pipeline Defines Functional Variants and Target Genes for Systemic Lupus Erythematosus

OBJECTIVE

Systemic lupus erythematosus (SLE), an autoimmune disease with incompletely understood etiology, has a strong genetic component. Although genome-wide association studies (GWASs) have revealed multiple SLE susceptibility loci and associated single-nucleotide polymorphisms (SNPs), the precise causal variants, target genes, cell types, tissues, and mechanisms of action remain largely unknown.

METHODS

Here, we report a comprehensive post-GWAS analysis using extensive bioinformatics, molecular modeling, and integrative functional genomic and epigenomic analyses to optimize fine-mapping. We compile and cross-reference immune cell-specific expression quantitative trait loci (cis- and trans-expression quantitative trait loci) with promoter capture high-throughput capture chromatin conformation (PCHi-C), allele-specific chromatin accessibility, and massively parallel reporter assay data to define predisposing variants and target genes. We experimentally validate a predicted locus using CRISPR/Cas9 genome editing, quantitative polymerase chain reaction, and Western blot.

RESULTS

Anchoring on 452 index SNPs, we selected 9,931 high linkage disequilibrium (r2 > 0.8) SNPs and defined 182 independent non-human leukocyte antigen (HLA) SLE loci. The 3,746 SNPs from 143 loci were identified as regulating 564 unique genes. Target genes are enriched in lupus-related tissues and associated with other autoimmune diseases. Of these, 329 SNPs (106 loci) showed significant allele-specific chromatin accessibility and/or enhancer activity, indicating regulatory potential. Using CRISPR/Cas9, we validated reference SNP identifier 57668933 (rs57668933) as a functional variant regulating multiple targets, including SLE-risk gene ELF1 in B cells.

CONCLUSION

We demonstrate and validate post-GWAS strategies for using multidimensional data to prioritize likely causal variants with cognate gene targets underlying SLE pathogenesis. Our results provide a catalog of significantly SLE-associated SNPs and loci, target genes, and likely biochemical mechanisms to guide experimental characterization.

lncRNA CDKN2B-AS1 regulates collagen expression

The long noncoding RNA CDKN2B-AS1 harbors a major coronary artery disease risk haplotype, which is also associated with progressive forms of the oral inflammatory disease periodontitis as well as myocardial infarction (MI). Despite extensive research, there is currently no broad consensus on the function of CDKN2B-AS1 that would explain a common molecular role of this lncRNA in these diseases. Our aim was to investigate the role of CDKN2B-AS1 in gingival cells to better understand the molecular mechanisms underlying the increased risk of progressive periodontitis. We downregulated CDKN2B-AS1 transcript levels in primary gingival fibroblasts with LNA GapmeRs. Following RNA-sequencing, we performed differential expression, gene set enrichment analyses and Western Blotting. Putative causal alleles were searched by analyzing associated DNA sequence variants for changes of predicted transcription factor binding sites. We functionally characterized putative functional alleles using luciferase-reporter and antibody electrophoretic mobility shift assays in gingival fibroblasts and HeLa cells. Of all gene sets analysed, collagen biosynthesis was most significantly upregulated (Padj=9.7 × 10- 5 (AUC > 0.65) with the CAD and MI risk gene COL4A1 showing strongest upregulation of the enriched gene sets (Fold change = 12.13, Padj = 4.9 × 10- 25). The inflammatory "TNFA signaling via NFKB" gene set was downregulated the most (Padj=1 × 10- 5 (AUC = 0.60). On the single gene level, CAPNS2, involved in extracellular matrix organization, was the top upregulated protein coding gene (Fold change = 48.5, P < 9 × 10- 24). The risk variant rs10757278 altered a binding site of the pathogen responsive transcription factor STAT1 (P = 5.8 × 10- 6). rs10757278-G allele reduced STAT1 binding 14.4% and rs10757278-A decreased luciferase activity in gingival fibroblasts 41.2% (P = 0.0056), corresponding with GTEx data. CDKN2B-AS1 represses collagen gene expression in gingival fibroblasts. Dysregulated collagen biosynthesis through allele-specific CDKN2B-AS1 expression in response to inflammatory factors may affect collagen synthesis, and in consequence tissue barrier and atherosclerotic plaque stability.

MRAS in coronary artery disease-Unchartered territory

Genome-wide association studies (GWAS) have identified coronary artery disease (CAD) susceptibility locus on chromosome 3q22.3. This locus contains a cluster of several genes that includes muscle rat sarcoma virus (MRAS). Common MRAS variants are also associated with CAD causing risk factors such as hypertension, dyslipidemia, obesity, and type II diabetes. The MRAS gene is an oncogene that encodes a membrane-bound small GTPase. It is involved in a variety of signaling pathways, regulating cell differentiation and cell survival (mitogen-activated protein kinase [MAPK]/extracellular signal-regulated kinase and phosphatidylinositol 3-kinase) as well as acute phase response signaling (tumor necrosis factor [TNF] and interleukin 6 [IL6] signaling). In this review, we will summarize the role of genetic MRAS variants in the etiology of CAD and its comorbidities with the focus on tissue distribution of MRAS isoforms, cell type/tissue specificity, and mode of action of single nucleotide variants in MRAS associated complex traits. Finally, we postulate that CAD risk variants in the MRAS locus are specific to smooth muscle cells and lead to higher levels of MRAS, particularly in arterial and cardiac tissue, resulting in MAPK-dependent tissue hypertrophy or hyperplasia.

Immunogenetics of lithium response and psychiatric phenotypes in patients with bipolar disorder

The link between bipolar disorder (BP) and immune dysfunction remains controversial. While epidemiological studies have long suggested an association, recent research has found only limited evidence of such a relationship. To clarify this, we performed an exploratory study of the contributions of immune-relevant genetic factors to the response to lithium (Li) treatment and the clinical presentation of BP. First, we assessed the association of a large collection of immune-related genes (4925) with Li response, defined by the Retrospective Assessment of the Lithium Response Phenotype Scale (Alda scale), and clinical characteristics in patients with BP from the International Consortium on Lithium Genetics (ConLi+Gen, N = 2374). Second, we calculated here previously published polygenic scores (PGSs) for immune-related traits and evaluated their associations with Li response and clinical features. Overall, we observed relatively weak associations (p < 1 × 10-4) with BP phenotypes within immune-related genes. Network and functional enrichment analyses of the top findings from the association analyses of Li response variables showed an overrepresentation of pathways participating in cell adhesion and intercellular communication. These appeared to converge on the well-known Li-induced inhibition of GSK-3β. Association analyses of age-at-onset, number of mood episodes, and presence of psychosis, substance abuse and/or suicidal ideation suggested modest contributions of genes such as RTN4, XKR4, NRXN1, NRG1/3 and GRK5 to disease characteristics. PGS analyses returned weak associations (p < 0.05) between inflammation markers and the studied BP phenotypes. Our results suggest a modest relationship between immunity and clinical features in BP. More research is needed to assess the potential therapeutic relevance.

A multilayered post-GWAS analysis pipeline defines functional variants and target genes for systemic lupus erythematosus (SLE)

Objectives

Systemic lupus erythematosus (SLE), an autoimmune disease with incompletely understood etiology, has a strong genetic component. Although genome-wide association studies (GWAS) have revealed multiple SLE susceptibility loci and associated single nucleotide polymorphisms (SNPs), the precise causal variants, target genes, cell types, tissues, and mechanisms of action remain largely unknown.

Methods

Here, we report a comprehensive post-GWAS analysis using extensive bioinformatics, molecular modeling, and integrative functional genomic and epigenomic analyses to optimize fine-mapping. We compile and cross-reference immune cell-specific expression quantitative trait loci ( cis - and trans -eQTLs) with promoter-capture Hi-C, allele-specific chromatin accessibility, and massively parallel reporter assay data to define predisposing variants and target genes. We experimentally validate a predicted locus using CRISPR/Cas9 genome editing, qPCR, and Western blot.

Results

Anchoring on 452 index SNPs, we selected 9,931 high-linkage disequilibrium (r 2 >0.8) SNPs and defined 182 independent non-HLA SLE loci. 3,746 SNPs from 143 loci were identified as regulating 564 unique genes. Target genes are enriched in lupus-related tissues and associated with other autoimmune diseases. Of these, 329 SNPs (106 loci) showed significant allele-specific chromatin accessibility and/or enhancer activity, indicating regulatory potential. Using CRISPR/Cas9, we validated rs57668933 as a functional variant regulating multiple targets, including SLE risk gene ELF1 , in B-cells.

Conclusion

We demonstrate and validate post-GWAS strategies for utilizing multi-dimensional data to prioritize likely causal variants with cognate gene targets underlying SLE pathogenesis. Our results provide a catalog of significantly SLE-associated SNPs and loci, target genes, and likely biochemical mechanisms, to guide experimental characterization.

The potential role of miRNAs and regulation of their expression in the development of mare endometrial fibrosis

Mare endometrial fibrosis (endometrosis), is one of the main causes of equine infertility. Despite the high prevalence, both ethology, pathogenesis and the nature of its progression remain poorly understood. Recent studies have shown that microRNAs (miRNAs) are important regulators in multiple cellular processes and functions under physiological and pathological circumstances. In this article, we reported changes in miRNA expression at different stages of endometrosis and the effect of transforming growth factor (TGF)-β1 on the expression of the most dysregulated miRNAs. We identified 1, 26, and 5 differentially expressed miRNAs (DEmiRs), in categories IIA (mild fibrosis), IIB (moderate fibrosis), and III (severe fibrosis) groups compared to category I (no fibrosis) endometria group, respectively (Padjusted < 0.05, log2FC ≥ 1.0/log2FC ≤  - 1.0). This study indicated the potential involvement of miRNAs in the regulation of the process associated to the development and progression of endometrosis. The functional enrichment analysis revealed, that DEmiRs target genes involved in the mitogen-activated protein kinases, Hippo, and phosphoinositide-3-kinase (PI3K)-Akt signalling pathways, focal adhesion, and extracellular matrix-receptor interaction. Moreover, we demonstrated that the most potent profibrotic cytokine-TGF-β1-downregulated novel-eca-miR-42 (P < 0.05) expression in fibroblasts derived from endometria at early-stage endometrosis (category IIA).

Immunogenetics of lithium response and psychiatric phenotypes in patients with bipolar disorder

The link between bipolar disorder (BP) and immune dysfunction remains controversial. While epidemiological studies have long suggested an association, recent research has found only limited evidence of such a relationship. To clarify this, we investigated the contributions of immune-relevant genetic factors to the response to lithium (Li) treatment and the clinical presentation of BP. First, we assessed the association of a large collection of immune-related genes (4,925) with Li response, defined by the Retrospective Assessment of the Lithium Response Phenotype Scale (Alda scale), and clinical characteristics in patients with BP from the International Consortium on Lithium Genetics (ConLi+Gen, N = 2,374). Second, we calculated here previously published polygenic scores (PGSs) for immune-related traits and evaluated their associations with Li response and clinical features. We found several genes associated with Li response at p < 1×10- 4 values, including HAS3, CNTNAP5 and NFIB. Network and functional enrichment analyses uncovered an overrepresentation of pathways involved in cell adhesion and intercellular communication, which appear to converge on the well-known Li-induced inhibition of GSK-3β. We also found various genes associated with BP's age-at-onset, number of mood episodes, and presence of psychosis, substance abuse and/or suicidal ideation at the exploratory threshold. These included RTN4, XKR4, NRXN1, NRG1/3 and GRK5. Additionally, PGS analyses suggested serum FAS, ECP, TRANCE and cytokine ligands, amongst others, might represent potential circulating biomarkers of Li response and clinical presentation. Taken together, our results support the notion of a relatively weak association between immunity and clinically relevant features of BP at the genetic level.

Multi-trait genome-wide association study identifies a novel endometrial cancer risk locus that associates with testosterone levels

To detect novel endometrial cancer risk variants, we leveraged information from endometrial cancer risk factors in a multi-trait GWAS analysis. We first assessed causal relationships between established and suspected endometrial cancer risk factors, and endometrial cancer using Mendelian randomization. Following multivariable analysis, five independent risk factors (waist circumference, testosterone levels, sex hormone binding globulin levels, age at menarche, and age at natural menopause) were included in a multi-trait Bayesian GWAS analysis. We identified three potentially novel loci that associate with endometrial cancer risk, one of which (7q22.1) replicated in an independent endometrial cancer GWAS dataset and was genome-wide significant in a meta-analysis. This locus may affect endometrial cancer risk through altered testosterone levels. Consistent with this, we observed colocalization between the signals for endometrial cancer risk and expression of CYP3A7, a gene involved in testosterone metabolism. Thus, our findings suggest opportunities for hormone therapy to prevent or treat endometrial cancer.

Genetic Regulation of SMC Gene Expression and Splicing Predict Causal CAD Genes

BACKGROUND

Coronary artery disease (CAD) is the leading cause of death worldwide. Recent meta-analyses of genome-wide association studies have identified over 175 loci associated with CAD. The majority of these loci are in noncoding regions and are predicted to regulate gene expression. Given that vascular smooth muscle cells (SMCs) play critical roles in the development and progression of CAD, we aimed to identify the subset of the CAD loci associated with the regulation of transcription in distinct SMC phenotypes.

METHODS

We measured gene expression in SMCs isolated from the ascending aortas of 151 heart transplant donors of various genetic ancestries in quiescent or proliferative conditions and calculated the association of their expression and splicing with ~6.3 million imputed single-nucleotide polymorphism markers across the genome.

RESULTS

We identified 4910 expression and 4412 splicing quantitative trait loci (sQTLs) representing regions of the genome associated with transcript abundance and splicing. A total of 3660 expression quantitative trait loci (eQTLs) had not been observed in the publicly available Genotype-Tissue Expression dataset. Further, 29 and 880 eQTLs were SMC-specific and sex-biased, respectively. We made these results available for public query on a user-friendly website. To identify the effector transcript(s) regulated by CAD loci, we used 4 distinct colocalization approaches. We identified 84 eQTL and 164 sQTL that colocalized with CAD loci, highlighting the importance of genetic regulation of mRNA splicing as a molecular mechanism for CAD genetic risk. Notably, 20% and 35% of the eQTLs were unique to quiescent or proliferative SMCs, respectively. One CAD locus colocalized with a sex-specific eQTL (TERF2IP), and another locus colocalized with SMC-specific eQTL (ALKBH8). The most significantly associated CAD locus, 9p21, was an sQTL for the long noncoding RNA CDKN2B-AS1, also known as ANRIL, in proliferative SMCs.

CONCLUSIONS

Collectively, our results provide evidence for the molecular mechanisms of genetic susceptibility to CAD in distinct SMC phenotypes.

Genome-wide association study for systemic lupus erythematosus in an egyptian population

Systemic lupus erythematosus (SLE) susceptibility has a strong genetic component. Genome-wide association studies (GWAS) across trans-ancestral populations show both common and distinct genetic variants of susceptibility across European and Asian ancestries, while many other ethnic populations remain underexplored. We conducted the first SLE GWAS on Egyptians-an admixed North African/Middle Eastern population-using 537 patients and 883 controls. To identify novel susceptibility loci and replicate previously known loci, we performed imputation-based association analysis with 6,382,276 SNPs while accounting for individual admixture. We validated the association analysis using adaptive permutation tests (n = 109). We identified a novel genome-wide significant locus near IRS1/miR-5702 (Pcorrected = 1.98 × 10-8) and eight novel suggestive loci (Pcorrected < 1.0 × 10-5). We also replicated (Pperm < 0.01) 97 previously known loci with at least one associated nearby SNP, with ITGAM, DEF6-PPARD and IRF5 the top three replicated loci. SNPs correlated (r 2 > 0.8) with lead SNPs from four suggestive loci (ARMC9, DIAPH3, IFLDT1, and ENTPD3) were associated with differential gene expression (3.5 × 10-95 < p < 1.0 × 10-2) across diverse tissues. These loci are involved in cellular proliferation and invasion-pathways prominent in lupus and nephritis. Our study highlights the utility of GWAS in an admixed Egyptian population for delineating new genetic associations and for understanding SLE pathogenesis.

A minority of somatically mutated genes in pre-existing fatty liver disease have prognostic importance in the development of NAFLD

BACKGROUND

Understanding the genetics of liver disease has the potential to facilitate clinical risk stratification. We recently identified acquired somatic mutations in six genes and one lncRNA in pre-existing fatty liver disease. We hypothesised that germline variation in these genes might be associated with the risk of developing steatosis and contribute to the prediction of disease severity.

METHODS

Genome-wide association study (GWAS) summary statistics were extracted from seven studies (>1.7 million participants) for variants near ACVR2A, ALB, CIDEB, FOXO1, GPAM, NEAT1 and TNRC6B for: aminotransferases, liver fat, HbA1c, diagnosis of NAFLD, ARLD and cirrhosis. Findings were replicated using GWAS data from multiple independent cohorts. A phenome-wide association study was performed to examine for related metabolic traits, using both common and rare variants, including gene-burden testing.

RESULTS

There was no evidence of association between rare germline variants or SNPs near five genes (ACVR2A, ALB, CIDEB, FOXO1 and TNRC6B) and risk or severity of liver disease. Variants in GPAM (proxies for p.Ile43Val) were associated with liver fat (p = 3.6 × 10-13 ), ALT (p = 2.8 × 10-39 ) and serum lipid concentrations. Variants in NEAT1 demonstrated borderline significant associations with ALT (p = 1.9 × 10-11 ) and HbA1c, but not with liver fat, as well as influencing waist-to-hip ratio, adjusted for BMI.

CONCLUSIONS

Despite the acquisition of somatic mutations at these loci during progressive fatty liver disease, we did not find associations between germline variation and markers of liver disease, except in GPAM. In the future, larger sample sizes may identify associations. Currently, germline polygenic risk scores will not capture data from genes affected by somatic mutations.

Microbiome-associated human genetic variants impact phenome-wide disease risk

Human genetic variation associates with the composition of the gut microbiome, yet its influence on clinical traits remains largely unknown. We analyzed the consequences of nearly a thousand gut microbiome-associated variants (MAVs) on phenotypes reported in electronic health records from tens of thousands of individuals. We discovered and replicated associations of MAVs with neurological, metabolic, digestive, and circulatory diseases. Five significant MAVs in these categories correlate with the relative abundance of microbes down to the strain level. We also demonstrate that these relationships are independently observed and concordant with microbe by disease associations reported in case-control studies. Moreover, a selective sweep and population differentiation impacted some disease-linked MAVs. Combined, these findings establish triad relationships among the human genome, microbiome, and disease. Consequently, human genetic influences may offer opportunities for precision diagnostics of microbiome-associated diseases but also highlight the relevance of genetic background for microbiome modulation and therapeutics.

Single-nucleus chromatin accessibility profiling highlights regulatory mechanisms of coronary artery disease risk

Coronary artery disease (CAD) is a complex inflammatory disease involving genetic influences across cell types. Genome-wide association studies (GWAS) have identified over 200 loci associated with CAD, where the majority of risk variants reside in noncoding DNA sequences impacting cis-regulatory elements (CREs). Here, we applied single-nucleus ATAC-seq to profile 28,316 nuclei across coronary artery segments from 41 patients with varying stages of CAD, which revealed 14 distinct cellular clusters. We mapped ~320,000 accessible sites across all cells, identified cell type-specific elements, transcription factors, and prioritized functional CAD risk variants. . We identified elements in smooth muscle cell (SMC) transition states (e.g. fibromyocytes) and functional variants predicted to alter SMC and macrophage-specific regulation of MRAS (3q22) and LIPA (10q23), respectively. We further nominated key driver transcription factors such as PRDM16 and TBX2. Together, this single nucleus atlas provides a critical step towards interpreting regulatory mechanisms across the continuum of CAD risk.

BACH1 Binding Links the Genetic Risk for Severe Periodontitis with ST8SIA1

Genome-wide association studies identified various loci associated with periodontal diseases, but assigning causal alleles remains difficult. Likewise, the generation of biological meaning underlying a statistical association has been challenging. Here, we characterized the genetic association at the gene ST8SIA1 that increases the risk for severe periodontitis in smokers. We used CRISPR/dCas9 activation and RNA-sequencing to identify genetic interaction partners of ST8SIA1 and to determine its function in the cell. We used reporter gene assays to identify regulatory elements at the associated single-nucleotide polymorphisms (SNPs) and to determine effect directions and allele-specific changes of enhancer activity. Antibody electrophoretic mobility shift assays proved allele-specific transcription factor binding at the putative causal SNPs. We found the reported periodontitis risk gene ABCA1 as the top upregulated gene following ST8SIA1 activation. Gene set enrichment analysis showed highest effects on integrin cell surface interactions (area under the curve [AUC] = 0.85; q = 4.9 × 10-6) and cell cycle regulation (AUC = 0.89; q = 1.6 × 10-5). We identified 2 associated repressor elements in the introns of ST8SIA1 that bind the transcriptional repressor BACH1. The putative causative variant rs2012722 decreased BACH1 binding by 40%. We also pinpointed ST8SIA1 as the target gene of the association. ST8SIA1 inhibits cell adhesion with extracellular matrix proteins, integrins, and cell cycle, as well as enhances apoptosis. Likewise, tobacco smoke reportedly results in an inhibition of cell adhesion and a decrease in integrin-positive cells and cell growth. We conclude that impaired ST8SIA1 repression, independently caused by reduced BACH1 binding at the effect T allele, as well as by tobacco smoke, contributes to higher ST8SIA1 levels, and in smokers who carry the effect T allele, both factors would be additive with damaging effects on the gingival barrier integrity. The activity of ST8SIA1 is also linked with the periodontitis risk gene ABCA1.

A compendium of uniformly processed human gene expression and splicing quantitative trait loci

Many gene expression quantitative trait locus (eQTL) studies have published their summary statistics, which can be used to gain insight into complex human traits by downstream analyses, such as fine mapping and co-localization. However, technical differences between these datasets are a barrier to their widespread use. Consequently, target genes for most genome-wide association study (GWAS) signals have still not been identified. In the present study, we present the eQTL Catalogue ( https://www.ebi.ac.uk/eqtl ), a resource of quality-controlled, uniformly re-computed gene expression and splicing QTLs from 21 studies. We find that, for matching cell types and tissues, the eQTL effect sizes are highly reproducible between studies. Although most QTLs were shared between most bulk tissues, we identified a greater diversity of cell-type-specific QTLs from purified cell types, a subset of which also manifested as new disease co-localizations. Our summary statistics are freely available to enable the systematic interpretation of human GWAS associations across many cell types and tissues.

Sex-specific genetic factors affect the risk of early-onset periodontitis in Europeans

AIMS

Various studies have reported that young European women are more likely to develop early-onset periodontitis compared to men. A potential explanation for the observed variations in sex and age of disease onset is the natural genetic variation within the autosomal genomes. We hypothesized that genotype-by-sex (G × S) interactions contribute to the increased prevalence and severity.

MATERIALS AND METHODS

Using the case-only design, we tested for differences in genetic effects between men and women in 896 North-West European early-onset cases, using imputed genotypes from the OmniExpress genotyping array. Population-representative 6823 controls were used to verify that the interacting variables G and S were uncorrelated in the general population.

RESULTS

In total, 20 loci indicated G × S associations (P < 0.0005), 3 of which were previously suggested as risk genes for periodontitis (ABLIM2, CDH13, and NELL1). We also found independent G × S interactions of the related gene paralogs MACROD1/FLRT1 (chr11) and MACROD2/FLRT3 (chr20). G × S-associated SNPs at CPEB4, CDH13, MACROD1, and MECOM were genome-wide-associated with heel bone mineral density (CPEB4, MECOM), waist-to-hip ratio (CPEB4, MACROD1), and blood pressure (CPEB4, CDH13).

CONCLUSIONS

Our results indicate that natural genetic variation affects the different heritability of periodontitis among sexes and suggest genes that contribute to inter-sex phenotypic variation in early-onset periodontitis.

Genetic analyses of gynecological disease identify genetic relationships between uterine fibroids and endometrial cancer, and a novel endometrial cancer genetic risk region at the WNT4 1p36.12 locus

Endometriosis, polycystic ovary syndrome (PCOS) and uterine fibroids have been proposed as endometrial cancer risk factors; however, disentangling their relationships with endometrial cancer is complicated due to shared risk factors and comorbidities. Using genome-wide association study (GWAS) data, we explored the relationships between these non-cancerous gynecological diseases and endometrial cancer risk by assessing genetic correlation, causal relationships and shared risk loci. We found significant genetic correlation between endometrial cancer and PCOS, and uterine fibroids. Adjustment for genetically predicted body mass index (a risk factor for PCOS, uterine fibroids and endometrial cancer) substantially attenuated the genetic correlation between endometrial cancer and PCOS but did not affect the correlation with uterine fibroids. Mendelian randomization analyses suggested a causal relationship between only uterine fibroids and endometrial cancer. Gene-based analyses revealed risk regions shared between endometrial cancer and endometriosis, and uterine fibroids. Multi-trait GWAS analysis of endometrial cancer and the genetically correlated gynecological diseases identified a novel genome-wide significant endometrial cancer risk locus at 1p36.12, which replicated in an independent endometrial cancer dataset. Interrogation of functional genomic data at 1p36.12 revealed biologically relevant genes, including WNT4 which is necessary for the development of the female reproductive system. In summary, our study provides genetic evidence for a causal relationship between uterine fibroids and endometrial cancer. It further provides evidence that the comorbidity of endometrial cancer, PCOS and uterine fibroids may partly be due to shared genetic architecture. Notably, this shared architecture has revealed a novel genome-wide risk locus for endometrial cancer.

Variation in genes implicated in B-cell development and antibody production affects susceptibility to pemphigus

Pemphigus foliaceus (PF) is an autoimmune blistering skin disease characterized by the presence of pathogenic autoantibodies against desmoglein 1, a component of intercellular desmosome junctions. PF occurs sporadically across the globe and is endemic in some Brazilian regions. Because PF is a B-cell-mediated disease, we aimed to study the impact of variants within genes encoding molecules involved in the different steps of B-cell development and antibody production on the susceptibility of endemic PF. We analysed 3,336 single nucleotide polymorphisms (SNPs) from 167 candidate genes genotyped with Illumina microarray in a cohort of 227 PF patients and 193 controls. After quality control and exclusion of non-informative and redundant SNPs, 607 variants in 149 genes remained in the logistic regression analysis, in which sex and ancestry were included as covariates. Our results revealed 10 SNPs within or nearby 11 genes that were associated with susceptibility to endemic PF (OR >1.56; p < 0.005): rs6657275*G (TGFB2); rs1818545*A (RAG1/RAG2/IFTAP);rs10781530*A (PAXX), rs10870140*G and rs10781522*A (TRAF2); rs535068*A (TNFRSF1B); rs324011*A (STAT6);rs6432018*C (YWHAQ); rs17149161*C (YWHAG); and rs2070729*C (IRF1). Interestingly, these SNPs have been previously associated with differential gene expression, mostly in peripheral blood, in publicly available databases. For the first time, we show that polymorphisms in genes involved in B-cell development and antibody production confer differential susceptibility to endemic PF, and therefore are candidates for possible functional studies to understand immunoglobulin gene rearrangement and its impact on diseases.