Comprehensive analysis of the major histocompatibility complex in systemic sclerosis identifies differential HLA associations by clinical and serological subtypes

Objective: The greatest genetic effect reported for systemic sclerosis (SSc) lies in the major histocompatibility complex (MHC) locus. Leveraging the largest SSc genome-wide association study, we aimed to fine-map this region to identify novel human leucocyte antigen (HLA) genetic variants associated with SSc susceptibility and its main clinical and serological subtypes. Methods: 9095 patients with SSc and 17 584 controls genome-wide genotyped were used to impute and test single-nucleotide polymorphisms (SNPs) across the MHC, classical HLA alleles and their composite amino acid residues. Additionally, patients were stratified according to their clinical and serological status, namely, limited cutaneous systemic sclerosis (lcSSc), diffuse cutaneous systemic sclerosis (dcSSc), anticentromere (ACA), antitopoisomerase (ATA) and anti-RNApolIII autoantibodies (ARA). Results: Sequential conditional analyses showed nine SNPs, nine classical alleles and seven amino acids that modelled the observed associations with SSc. This confirmed previously reported associations with HLA-DRB1*11:04 and HLA-DPB1*13:01, and revealed a novel association of HLA-B*08:01. Stratified analyses showed specific associations of HLA-DQA1*02:01 with lcSSc, and an exclusive association of HLA-DQA1*05:01 with dcSSc. Similarly, private associations were detected in HLA-DRB1*08:01 and confirmed the previously reported association of HLA-DRB1*07:01 with ACA-positive patients, as opposed to the HLA-DPA1*02:01 and HLA-DQB1*03:01 alleles associated with ATA presentation. Conclusions: This study confirms the contribution of HLA class II and reveals a novel association of HLA class I with SSc, suggesting novel pathways of disease pathogenesis. Furthermore, we describe specific HLA associations with SSc clinical and serological subtypes that could serve as biomarkers of disease severity and progression.Funding: This work was supported by the Spanish Ministry of Science and Innovation (grant ref. SAF2015-66761-P and RTI20181013 (32-B-100)), Red de Investigación en Inflamación y Enfermedades Reumáticas from Instituto de Salud Carlos III (RD16/0012/0013) and grants from National Institutes of Health (R01AR073284) and DoD (W81XWH-16-1-0296). MAH was funded by the Spanish Ministry of Science and Innovation through the Juan de la Cierva Incorporacion program (ref. IJC2018-035131-I). GO, AB and ALH were supported by the NIHR Manchester Biomedical Research Centre and Versus Arthritis (grant ref 21754)

Genome-Wide Association Study and Gene Expression Analysis Identifies CD84 as a Predictor of Response to Etanercept Therapy in Rheumatoid Arthritis

Anti-tumor necrosis factor alpha (anti-TNF) biologic therapy is a widely used treatment for rheumatoid arthritis (RA). It is unknown why some RA patients fail to respond adequately to anti-TNF therapy, which limits the development of clinical biomarkers to predict response or new drugs to target refractory cases. To understand the biological basis of response to anti-TNF therapy, we conducted a genome-wide association study (GWAS) meta-analysis of more than 2 million common variants in 2,706 RA patients from 13 different collections. Patients were treated with one of three anti-TNF medications: etanercept (n = 733), infliximab (n = 894), or adalimumab (n = 1,071). We identified a SNP (rs6427528) at the 1q23 locus that was associated with change in disease activity score (ΔDAS) in the etanercept subset of patients (P = 8×10-8), but not in the infliximab or adalimumab subsets (P>0.05). The SNP is predicted to disrupt transcription factor binding site motifs in the 3′ UTR of an immune-related gene, CD84, and the allele associated with better response to etanercept was associated with higher CD84 gene expression in peripheral blood mononuclear cells (P = 1×10-11 in 228 non-RA patients and P = 0.004 in 132 RA patients). Consistent with the genetic findings, higher CD84 gene expression correlated with lower cross-sectional DAS (P = 0.02, n = 210) and showed a non-significant trend for better ΔDAS in a subset of RA patients with gene expression data (n = 31, etanercept-treated). A small, multi-ethnic replication showed a non-significant trend towards an association among etanercept-treated RA patients of Portuguese ancestry (n = 139, P = 0.4), but no association among patients of Japanese ancestry (n = 151, P = 0.8). Our study demonstrates that an allele associated with response to etanercept therapy is also associated with CD84 gene expression, and further that CD84 expression correlates with disease activity. These findings support a model in which CD84 genotypes and/or expression may serve as a useful biomarker for response to etanercept treatment in RA patients of European ancestry. © 2013 Cui et al

Functional impact and evolution of a novel human polymorphic inversion that disrupts a gene and creates a fusion transcript

Since the discovery of chromosomal inversions almost 100 years ago, how they are maintained in natural populations has been a highly debated issue. One of the hypotheses is that inversion breakpoints could affect genes and modify gene expression levels, although evidence of this came only from laboratory mutants. In humans, a few inversions have been shown to associate with expression differences, but in all cases the molecular causes have remained elusive. Here, we have carried out a complete characterization of a new human polymorphic inversion and determined that it is specific to East Asian populations. In addition, we demonstrate that it disrupts the ZNF257 gene and, through the translocation of the first exon and regulatory sequences, creates a previously nonexistent fusion transcript, which together are associated to expression changes in several other genes. Finally, we investigate the potential evolutionary and phenotypic consequences of the inversion, and suggest that it is probably deleterious. This is therefore the first example of a natural polymorphic inversion that has position effects and creates a new chimeric gene, contributing to answer an old question in evolutionary biology

Implicating genes, pleiotropy, and sexual dimorphism at blood lipid loci through multi-ancestry meta-analysis

Background: Genetic variants within nearly 1000 loci are known to contribute to modulation of blood lipid levels. However, the biological pathways underlying these associations are frequently unknown, limiting understanding of these findings and hindering downstream translational efforts such as drug target discovery. Results: To expand our understanding of the underlying biological pathways and mechanisms controlling blood lipid levels, we leverage a large multi-ancestry meta-analysis (N = 1,654,960) of blood lipids to prioritize putative causal genes for 2286 lipid associations using six gene prediction approaches. Using phenome-wide association (PheWAS) scans, we identify relationships of genetically predicted lipid levels to other diseases and conditions. We confirm known pleiotropic associations with cardiovascular phenotypes and determine novel associations, notably with cholelithiasis risk. We perform sex-stratified GWAS meta-analysis of lipid levels and show that 3-5% of autosomal lipid-associated loci demonstrate sex-biased effects. Finally, we report 21 novel lipid loci identified on the X chromosome. Many of the sex-biased autosomal and X chromosome lipid loci show pleiotropic associations with sex hormones, emphasizing the role of hormone regulation in lipid metabolism. Conclusions: Taken together, our findings provide insights into the biological mechanisms through which associated variants lead to altered lipid levels and potentially cardiovascular disease risk

Deciphering osteoarthritis genetics across 826,690 individuals from 9 populations

Osteoarthritis affects over 300 million people worldwide. Here, we conduct a genome-wide association study meta-analysis across 826,690 individuals (177,517 with osteoarthritis) and identify 100 independently associated risk variants across 11 osteoarthritis phenotypes, 52 of which have not been associated with the disease before. We report thumb and spine osteoarthritis risk variants and identify differences in genetic effects between weight-bearing and non-weight-bearing joints. We identify sex-specific and early age-at-onset osteoarthritis risk loci. We integrate functional genomics data from primary patient tissues (including articular cartilage, subchondral bone, and osteophytic cartilage) and identify high-confidence effector genes. We provide evidence for genetic correlation with phenotypes related to pain, the main disease symptom, and identify likely causal genes linked to neuronal processes. Our results provide insights into key molecular players in disease processes and highlight attractive drug targets to accelerate translation

Genetics of rheumatoid arthritis contributes to biology and drug discovery

A major challenge in human genetics is to devise a systematic strategy to integrate disease-associated variants with diverse genomic and biological datasets to provide insight into disease pathogenesis and guide drug discovery for complex traits such as rheumatoid arthritis (RA)1. Here, we performed a genome-wide association study (GWAS) meta-analysis in a total of >100,000 subjects of European and Asian ancestries (29,880 RA cases and 73,758 controls), by evaluating ~10 million single nucleotide polymorphisms (SNPs). We discovered 42 novel RA risk loci at a genome-wide level of significance, bringing the total to 1012–4. We devised an in-silico pipeline using established bioinformatics methods based on functional annotation5, cis-acting expression quantitative trait loci (cis-eQTL)6, and pathway analyses7–9 – as well as novel methods based on genetic overlap with human primary immunodeficiency (PID), hematological cancer somatic mutations and knock-out mouse phenotypes – to identify 98 biological candidate genes at these 101 risk loci. We demonstrate that these genes are the targets of approved therapies for RA, and further suggest that drugs approved for other indications may be repurposed for the treatment of RA. Together, this comprehensive genetic study sheds light on fundamental genes, pathways and cell types that contribute to RA pathogenesis, and provides empirical evidence that the genetics of RA can provide important information for drug discovery

Genetic predisposition to mosaic Y chromosome loss in blood.

Mosaic loss of chromosome Y (LOY) in circulating white blood cells is the most common form of clonal mosaicism1-5, yet our knowledge of the causes and consequences of this is limited. Here, using a computational approach, we estimate that 20% of the male population represented in the UK Biobank study (n = 205,011) has detectable LOY. We identify 156 autosomal genetic determinants of LOY, which we replicate in 757,114 men of European and Japanese ancestry. These loci highlight genes that are involved in cell-cycle regulation and cancer susceptibility, as well as somatic drivers of tumour growth and targets of cancer therapy. We demonstrate that genetic susceptibility to LOY is associated with non-haematological effects on health in both men and women, which supports the hypothesis that clonal haematopoiesis is a biomarker of genomic instability in other tissues. Single-cell RNA sequencing identifies dysregulated expression of autosomal genes in leukocytes with LOY and provides insights into why clonal expansion of these cells may occur. Collectively, these data highlight the value of studying clonal mosaicism to uncover fundamental mechanisms that underlie cancer and other ageing-related diseases.This research has been conducted using the UK Biobank Resource under application 9905 and 19808. This work was supported by the Medical Research Council [Unit Programme number MC_UU_12015/2]. Full study-specific and individual acknowledgements can be found in the supplementary information

Functional annotation of human long noncoding RNAs via molecular phenotyping

Long noncoding RNAs (lncRNAs) constitute the majority of transcripts in the mammalian genomes, and yet, their functions remain largely unknown. As part of the FANTOM6 project, we systematically knocked down the expression of 285 lncRNAs in human dermal fibroblasts and quantified cellular growth, morphological changes, and transcriptomic responses using Capped Analysis of Gene Expression (CAGE). Antisense oligonucleotides targeting the same lncRNAs exhibited global concordance, and the molecular phenotype, measured by CAGE, recapitulated the observed cellular phenotypes while providing additional insights on the affected genes and pathways. Here, we disseminate the largest-todate lncRNA knockdown data set with molecular phenotyping (over 1000 CAGE deep-sequencing libraries) for further exploration and highlight functional roles for ZNF213-AS1 and lnc-KHDC3L-2.Peer reviewe

Contribution of Common Genetic Variants to Risk of Early-Onset Ischemic Stroke

Background and Objectives Current genome-wide association studies of ischemic stroke have focused primarily on late-onset disease. As a complement to these studies, we sought to identify the contribution of common genetic variants to risk of early-onset ischemic stroke. Methods We performed a meta-analysis of genome-wide association studies of early-onset stroke (EOS), ages 18-59 years, using individual-level data or summary statistics in 16,730 cases and 599,237 nonstroke controls obtained across 48 different studies. We further compared effect sizes at associated loci between EOS and late-onset stroke (LOS) and compared polygenic risk scores (PRS) for venous thromboembolism (VTE) between EOS and LOS. Results We observed genome-wide significant associations of EOS with 2 variants in ABO, a known stroke locus. These variants tag blood subgroups O1 and A1, and the effect sizes of both variants were significantly larger in EOS compared with LOS. The odds ratio (OR) for rs529565, tagging O1, was 0.88 (95% confidence interval [CI]: 0.85-0.91) in EOS vs 0.96 (95% CI: 0.92-1.00) in LOS, and the OR for rs635634, tagging A1, was 1.16 (1.11-1.21) for EOS vs 1.05 (0.99-1.11) in LOS; p-values for interaction = 0.001 and 0.005, respectively. Using PRSs, we observed that greater genetic risk for VTE, another prothrombotic condition, was more strongly associated with EOS compared with LOS (p = 0.008). Discussion The ABO locus, genetically predicted blood group A, and higher genetic propensity for venous thrombosis are more strongly associated with EOS than with LOS, supporting a stronger role of prothrombotic factors in EOS.Peer reviewe