The landscape of expression and alternative splicing variation across human traits

Understanding the consequences of individual transcriptome variation is fundamental to deciphering human biology and disease. We implement a statistical framework to quantify the contributions of 21 individual traits as drivers of gene expression and alternative splicing variation across 46 human tissues and 781 individuals from the Genotype-Tissue Expression project. We demonstrate that ancestry, sex, age, and BMI make additive and tissue-specific contributions to expression variability, whereas interactions are rare. Variation in splicing is dominated by ancestry and is under genetic control in most tissues, with ribosomal proteins showing a strong enrichment of tissue-shared splicing events. Our analyses reveal a systemic contribution of types 1 and 2 diabetes to tissue transcriptome variation with the strongest signal in the nerve, where histopathology image analysis identifies novel genes related to diabetic neuropathy. Our multi-tissue and multi-trait approach provides an extensive characterization of the main drivers of human transcriptome variation in health and disease.This study was funded by the HumTranscriptom project with reference PID2019-107937GA-I00. R.G.-P. was supported by a Juan de la Cierva fellowship (FJC2020-044119-I) funded by MCIN/AEI/10.13039/501100011033 and ‘‘European Union NextGenerationEU/PRTR.’’ J.M.R. was supported by a predoctoral fellowship from ‘‘la Caixa’’ Foundation (ID 100010434) with code LCF/BQ/DR22/11950022. A.R.-C. was supported by a Formación Personal Investigador (FPI) fellowship (PRE2019-090193) funded by MCIN/AEI. R.C.-G. was supported by an FPI fellowship (PRE2020-092510) funded by MCIN/AEI. M.M. was supported by a Ramon y Cajal fellowship (RYC-2017-22249).Peer ReviewedPostprint (published version

Efficacy of Anakinra in Refractory Adult-Onset Still's Disease: Multicenter Study of 41 Patients and Literature Review

Adult-onset Still's disease (AOSD) is often refractory to standard therapy. Anakinra (ANK), an interleukin-1 receptor antagonist, has demonstrated efficacy in single cases and small series of AOSD. We assessed the efficacy of ANK in a series of AOSD patients. Multicenter retrospective open-label study. ANK was used due to lack of efficacy to standard synthetic immunosuppressive drugs and in some cases also to at least 1 biologic agent. Forty-one patients (26 women/15 men) were recruited. They had a mean age of 34.4 ± 14 years and a median [interquartile range (IQR)] AOSD duration of 3.5 [2-6] years before ANK onset. At that time the most common clinical features were joint manifestations 87.8%, fever 78%, and cutaneous rash 58.5%. ANK yielded rapid and maintained clinical and laboratory improvement. After 1 year of therapy, the frequency of joint and cutaneous manifestations had decreased to 41.5% and to 7.3% respectively, fever from 78% to 14.6%, anemia from 56.1% to 9.8%, and lymphadenopathy from 26.8% to 4.9%. A dramatic improvement of laboratory parameters was also achieved. The median [IQR] prednisone dose was also reduced from 20 [11.3-47.5] mg/day at ANK onset to 5 [0-10] at 12 months. After a median [IQR] follow-up of 16 [5-50] months, the most important side effects were cutaneous manifestations (n = 8), mild leukopenia (n = 3), myopathy (n = 1), and infections (n = 5). ANK is associated with rapid and maintained clinical and laboratory improvement, even in nonresponders to other biologic agents. However, joint manifestations are more refractory than the systemic manifestations

Association of HLA-B*41:02 with Henoch-Schönlein Purpura (IgA Vasculitis) in Spanish individuals irrespective of the HLA-DRB1 status

INTRODUCTION: A study was conducted to determine whether the human leukocyte antigen (HLA) B alleles are implicated in the susceptibility to Henoch-Schönlein purpura (HSP) in the largest series of Caucasian HSP patients ever assessed for genetic studies. METHODS: The study population was composed of 349 Spanish patients diagnosed with HSP fulfilling the American College of Rheumatology and the Michel et al. classification criteria, and 335 sex and ethnically matched controls. HLA-B phenotypes were determined by sequencing-based typing (SBT) and analyzed by chi-square or Fisher exact test. RESULTS: A statistically significant increase of HLA-B*41:02 allele in HSP patients when compared with controls was found (8.3% versus 1.5% respectively; P = 0.0001; OR (odds ratio) =5.76 [2.15-19.3]). These results remained statistically significant after adjusting for Bonferroni correction (P = 0.0028). An internal validation also confirmed the susceptibility effect on HSP associated with HLA-B*41:02 (OR = 5.70 [1.98-16.44]). Since a former study described an association between HLA-DRB1*01:03 and HSP susceptibility, we also evaluated the implication of HLA-B*41:02 independently of HLA-DRB1*01:03. Interestingly, the association remained statistically significant (P = 0.0004, OR = 4.97 [1.8-16.9]). No HLA-B association with specific HSP clinical features was found. CONCLUSIONS: Our study indicates that HLA-B*41:02 is associated with the susceptibility to HSP in Spanish patients irrespective of HLA-DRB1 status.This study was supported by a grant from ‘Fondo de Investigaciones Sanitarias’ PI12/00193 (Spain). RLM is a recipient of a Sara Borrell postdoctoral fellowship from the Instituto de Salud Carlos III at the Spanish Ministry of Health (Spain) (CD12/00425). FG and BU are supported by funds from the RETICS Program (RIER) (RD12/0009/0013).Ye

HLA-DRB1 association with Henoch-Schonlein purpura

Objective: Henoch-Schönlein purpura (HSP) is the most common vasculitis in children but it is not exceptional in adults. Increased familial occurrence supports a genetic predisposition for HSP. In this context, an association with the human leukocyte antigen-HLA-DRB1*01 phenotype has been suggested in Caucasian individuals with HSP. However, data on the potential association of HSP with HLA-DRB1*01 were based on small case series. To further investigate this issue, we performed HLA-DRB1 genotyping of the largest series of HSP patients ever assessed for genetic studies in Caucasians. Methods: 342 Spanish patients diagnosed with HSP fulfilling the American College of Rheumatology and the Michel et al classification criteria, and 303 sex and ethnically matched controls were assessed. HLA-DRB1 alleles were determined using a PCR-Sequence-Specific-Oligonucleotide Probe (PCR-SSOP) method. Results: A statistically significant increase of HLA-DRB1*01 in HSP patients when compared with controls was found (43% vs 7%, respectively; p<0.001; odds ratio-OR=2.03 [1.43-2.87]). It was due to the increased frequency of HLA-DRB1*0103 phenotype in HSP (14% vs 2%; p<0.001; OR=8.27 [3.46-23.9]). These results remained statistically significant after adjusting for Bonferroni correction. In contrast, a statistically significant decreased frequency of the HLA-DRB1*0301 phenotype was observed in patients compared to controls (5.6% vs 18.1%, respectively; p<0.001, OR=0.26 [0.14-0.47]), even after adjustment for Bonferroni correction. No HLA-DRB1 association with specific features of the disease was found. Conclusion: Our study confirms an association of HSP with HLA-DRB1*01 in Caucasians. Also, a protective effect against the development of HSP appears to exist in Caucasians carrying the HLA-DRB1*03 phenotype

HLA association with the susceptibility to anti-synthetase syndrome

Objective: To investigate the human leukocyte antigen (HLA) association with anti-synthetase syndrome (ASSD). Methods: We conducted the largest immunogenetic HLA-DRB1 and HLA-B study to date in a homogeneous cohort of 168 Caucasian patients with ASSD and 486 ethnically matched healthy controls by sequencing-based-typing. Results: A statistically significant increase of HLA-DRB1*03:01 and HLA-B*08:01 alleles in patients with ASSD compared to healthy controls was disclosed (26.2% versus 12.2%, P=1.56E-09, odds ratio-OR [95% confidence interval-CI]=2.54 [1.84-3.50] and 21.4% versus 5.5%, P=18.95E-18, OR [95% CI]=4.73 [3.18-7.05]; respectively). Additionally, HLA-DRB1*07:01 allele was significantly decreased in patients with ASSD compared to controls (9.2% versus 17.5%, P=0.0003, OR [95% CI]=0.48 [0.31-0.72]). Moreover, a statistically significant increase of HLA-DRB1*03:01 allele in anti-Jo-1 positive compared to anti-Jo-1 negative patients with ASSD was observed (31.8% versus 15.5%, P=0.001, OR [95% CI]=2.54 [1.39-4.81]). Similar findings were observed when HLA carrier frequencies were assessed. The HLA-DRB1*03:01 association with anti-Jo-1 was unrelated to smoking history. No HLA differences in patients with ASSD stratified according to the presence/absence of the most representative non-anti-Jo-1 anti-synthetase autoantibodies (anti-PL-12 and anti-PL-7), arthritis, myositis or interstitial lung disease were observed. Conclusions: Our results support the association of the HLA complex with the susceptibility to ASSD

Influence of MUC5B gene on antisynthetase syndrome

ABSTRACT: MUC5B rs35705950 (G/T) is strongly associated with idiopathic pulmonary fibrosis (IPF) and also contributes to the risk of interstitial lung disease (ILD) in rheumatoid arthritis (RA-ILD) and chronic hypersensitivity pneumonitis (CHP). Due to this, we evaluated the implication of MUC5B rs35705950 in antisynthetase syndrome (ASSD), a pathology characterised by a high ILD incidence. 160 patients with ASSD (142 with ILD associated with ASSD [ASSD-ILD+]), 232 with ILD unrelated to ASSD (comprising 161 IPF, 27 RA-ILD and 44 CHP) and 534 healthy controls were genotyped. MUC5B rs35705950 frequency did not significantly differ between ASSD-ILD+ patients and healthy controls nor when ASSD patients were stratified according to the presence/absence of anti Jo-1 antibodies or ILD. No significant differences in MUC5B rs35705950 were also observed in ASSD-ILD+ patients with a usual interstitial pneumonia (UIP) pattern when compared to those with a non-UIP pattern. However, a statistically significant decrease of MUC5B rs35705950 GT, TT and T frequencies in ASSD-ILD+ patients compared to patients with ILD unrelated to ASSD was observed. In summary, our study does not support a role of MUC5B rs35705950 in ASSD. It also indicates that there are genetic differences between ILD associated with and that unrelated to ASSD.We are indebted to the patients and healthy controls for their essential collaboration to this study. We also thank the National DNA Bank Repository (Salamanca) for supplying part of the control samples. This study was partially supported by grants from the Foundation for Research in Rheumatology (FOREUM). RL-M is a recipient of a Miguel Servet type I programme fellowship from the ‘Instituto de Salud Carlos III’ (ISCIII), co-funded by the European Social Fund (ESF, ‘Investing in your future’) (grant CP16/00033). SR-M is supported by funds of the RETICS Program (RD16/0012/0009), co-funded by the European Regional Development Fund (ERDF). VP-C is supported by a pre-doctoral grant from IDIVAL (PREVAL 18/01). VM is supported by funds of a Miguel Servet type I programme (grant CP16/00033) (ISCIII, co-funded by ESF). LL-G is supported by funds of PI18/00042 (ISCIII, co-funded by ERDF). OG is Staff Personnel of Xunta de Galicia (Servizo Galego de Saude, SERGAS) through a research-staff stabilization contract (ISCIII/SERGAS). OG,is member of RETICS Programme, RD16/0012/0014 (RIER: Red de Investigación en Inflamación y Enfermedades Reumáticas) via Instituto de Salud Carlos III (ISCIII) and FEDER. The work of OG (PI17/00409), was funded by Instituto de Salud Carlos III and FEDER. OG is a beneficiary of a project funded by Research Executive Agency of the European Union in the framework of MSCA-RISE Action of the H2020 Programme (Project number 734899). OG is beneficiary of a grant funded by Xunta de Galicia, Consellería de Educación, Universidade e Formación Profesional and Consellería de Economía, Emprego e Industria (GAIN), GPC IN607B2019/10

Role of PTPN22 and CSK gene polymorphisms as predictors of susceptibility and clinical heterogeneity in patients with Henoch-Schönlein purpura (IgA vasculitis)

INTRODUCTION: To determine whether the PTPN22 (protein tyrosine phosphatase nonreceptor 22)/CSK (c-src tyrosine kinase) pathway is implicated in the susceptibility and clinical heterogeneity of Henoch-Schönlein purpura (HSP) in the largest series of Caucasian HSP patients ever assessed for genetic studies. METHODS: A set of 329 Spanish patients diagnosed with HSP fulfilling the American College of Rheumatology and the Michel et al. classification criteria and 515 sex and ethnically matched controls were recruited in this study. Two well-known CSK (CSK rs34933034 and CSK rs1378942) and two functional PTPN22 (PTPN22 rs2476601 (R620W) and PTPN22 rs33996649 (R263Q)) polymorphisms, previously associated with autoimmunity, were genotyped with TaqMan single nucleotide polymorphism (SNP) genotyping assays. RESULTS: No significant differences in the genotype and allele frequencies between HSP patients and controls were observed when the CSK rs34933034, CSK rs1378942, PTPN22 rs2476601 (R620W) and PTPN22 rs33996649 (R263Q) polymorphisms were analyzed independently. In keeping with this observation, no significant differences were found when we assessed these polymorphisms combined conforming haplotypes. In addition, there were no differences in the allele or genotype frequencies when HSP patients were stratified according the age at disease onset, sex, presence of arthralgia/arthritis, nephritis or gastrointestinal manifestations. CONCLUSIONS: Our results do not support association between PTPN22/CSK and HSP

Genetics of immunoglobulin-A vasculitis (Henoch-Schönlein purpura): An updated review

Immunoglobulin-A vasculitis (IgAV) is classically a childhood small-sized blood vessel vasculitis with predominant involvement of the skin. Gastrointestinal and joint manifestations are common in patients diagnosed with this condition. Nephritis, which is more severe in adults, constitutes the most feared complication of this vasculitis. The molecular bases underlying the origin of IgAV have not been completely elucidated. Nevertheless, several pieces of evidence support the claim that genes play a crucial role in the pathogenesis of this disease. The human leukocyte antigen (HLA) region is, until now, the main genetic factor associated with IgAV pathogenesis. Besides a strong association with HLA class II alleles, specifically HLA-DRB1 alleles, HLA class I alleles also seem to influence on the predisposition of this disease. Other gene polymorphisms located outside the HLA region, including those coding cytokines, chemokines, adhesion molecules as well as those related to T-cells, aberrant glycosylation of IgA1, nitric oxide production, neoangiogenesis, renin-angiotensin system and lipid, Pyrin and homocysteine metabolism, may be implicated not only in the predisposition to IgAV but also in its severity. An update of the current knowledge of the genetic component associated with the pathogenesis of IgAV is detailed in this review.Acknowledgements: RL-Mis supported by the Miguel Servet I programme of the Spanish Ministry of Economy and Competitiveness through the grant CP16/ 00033. FG is recipient of a Sara Borrell postdoctoral fellowship from the “Instituto Carlos III de Salud” at the Spanish Ministry of Health (Spain) (CD15/00095). SR-M is supported by funds from the RETICS Program (RIER) (RD16/0012/0009). FDC is supported by the Ramón y Cajal programme of the Spanish Ministry of Economy and Competitiveness through the grant RYC-2014-16458

Germline variation at 8q24 and prostate cancer risk in men of European ancestry

Chromosome 8q24 is a susceptibility locus for multiple cancers, including prostate cancer. Here we combine genetic data across the 8q24 susceptibility region from 71,535 prostate cancer cases and 52,935 controls of European ancestry to define the overall contribution of germline variation at 8q24 to prostate cancer risk. We identify 12 independent risk signals for prostate cancer (p < 4.28 × 10−15), including three risk variants that have yet to be reported. From a polygenic risk score (PRS) model, derived to assess the cumulative effect of risk variants at 8q24, men in the top 1% of the PRS have a 4-fold (95%CI = 3.62–4.40) greater risk compared to the population average. These 12 variants account for ~25% of what can be currently explained of the familial risk of prostate cancer by known genetic risk factors. These findings highlight the overwhelming contribution of germline variation at 8q24 on prostate cancer risk which has implications for population risk stratification

Fine-mapping of prostate cancer susceptibility loci in a large meta-analysis identifies candidate causal variants

Prostate cancer is a polygenic disease with a large heritable component. A number of common, low-penetrance prostate cancer risk loci have been identified through GWAS. Here we apply the Bayesian multivariate variable selection algorithm JAM to fine-map 84 prostate cancer susceptibility loci, using summary data from a large European ancestry meta-analysis. We observe evidence for multiple independent signals at 12 regions and 99 risk signals overall. Only 15 original GWAS tag SNPs remain among the catalogue of candidate variants identified; the remainder are replaced by more likely candidates. Biological annotation of our credible set of variants indicates significant enrichment within promoter and enhancer elements, and transcription factor-binding sites, including AR, ERG and FOXA1. In 40 regions at least one variant is colocalised with an eQTL in prostate cancer tissue. The refined set of candidate variants substantially increase the proportion of familial relative risk explained by these known susceptibility regions, which highlights the importance of fine-mapping studies and has implications for clinical risk profiling. © 2018 The Author(s).Prostate cancer is a polygenic disease with a large heritable component. A number of common, low-penetrance prostate cancer risk loci have been identified through GWAS. Here we apply the Bayesian multivariate variable selection algorithm JAM to fine-map 84 prostate cancer susceptibility loci, using summary data from a large European ancestry meta-analysis. We observe evidence for multiple independent signals at 12 regions and 99 risk signals overall. Only 15 original GWAS tag SNPs remain among the catalogue of candidate variants identified; the remainder are replaced by more likely candidates. Biological annotation of our credible set of variants indicates significant enrichment within promoter and enhancer elements, and transcription factor-binding sites, including AR, ERG and FOXA1. In 40 regions at least one variant is colocalised with an eQTL in prostate cancer tissue. The refined set of candidate variants substantially increase the proportion of familial relative risk explained by these known susceptibility regions, which highlights the importance of fine-mapping studies and has implications for clinical risk profiling. © 2018 The Author(s).Peer reviewe