Article Text
Abstract
Objectives Systemic lupus erythematosus (SLE) is an autoimmune disease with extensive heterogeneity in disease presentation between patients, which is likely due to an underlying molecular diversity. Here, we aimed at elucidating the genetic aetiology of SLE from the immunity pathway level to the single variant level, and stratify patients with SLE into distinguishable molecular subgroups, which could inform treatment choices in SLE.
Methods We undertook a pathway-centred approach, using sequencing of immunological pathway genes. Altogether 1832 candidate genes were analysed in 958 Swedish patients with SLE and 1026 healthy individuals. Aggregate and single variant association testing was performed, and we generated pathway polygenic risk scores (PRS).
Results We identified two main independent pathways involved in SLE susceptibility: T lymphocyte differentiation and innate immunity, characterised by HLA and interferon, respectively. Pathway PRS defined pathways in individual patients, who on average were positive for seven pathways. We found that SLE organ damage was more pronounced in patients positive for the T or B cell receptor signalling pathways. Further, pathway PRS-based clustering allowed stratification of patients into four groups with different risk score profiles. Studying sets of genes with priors for involvement in SLE, we observed an aggregate common variant contribution to SLE at genes previously reported for monogenic SLE as well as at interferonopathy genes.
Conclusions Our results show that pathway risk scores have the potential to stratify patients with SLE beyond clinical manifestations into molecular subsets, which may have implications for clinical follow-up and therapy selection.
- lupus erythematosus
- systemic
- polymorphism
- genetic
- autoimmunity
This is an open access article distributed in accordance with the Creative Commons Attribution 4.0 Unported (CC BY 4.0) license, which permits others to copy, redistribute, remix, transform and build upon this work for any purpose, provided the original work is properly cited, a link to the licence is given, and indication of whether changes were made. See: https://creativecommons.org/licenses/by/4.0/.
Statistics from Altmetric.com
Key messages
What is already known about this subject?
The clinical heterogeneity in systemic lupus erythematosus (SLE) is likely due to an underlying molecular diversity that could have implications for therapy.
In recent years, gene expression, autoantibody profiles and cytokine levels have been used to identify groups of patients with SLE with distinct molecular disease mechanisms.
What does this study add?
We have presented a novel strategy to genetically stratify SLE patients according to involved molecular pathways.
Using genetic information to stratify patients would have the advantages of providing stable molecular markers for early classification.
How might this impact on clinical practice or future developments?
Our results show that pathway risk scores have the potential to stratify SLE patients beyond clinical manifestations into molecular subsets, which may have implications for clinical follow-up and therapy selection.
Introduction
Systemic lupus erythematosus (SLE) is characterised by the production of autoantibodies targeting nucleic acids and associated proteins, immune complex formation and inflammation in multiple organs. There is a wide spectrum of clinical manifestations in SLE and extensive heterogeneity in disease presentation between patients; in addition, the treatment response is often unpredictable.1 The pathogenesis of SLE has partially been clarified during the last years, and important features are increased expression of type I interferon (IFN) regulated genes, defects in the apoptotic process and activated autoreactive B cells.1 2 The reasons behind these abnormalities are both environmental and genetic, and today around 100 SLE susceptibility loci have been identified.3 4 Monogenic forms of SLE exist, but for a majority of patients the environment and the cumulative number of susceptibility alleles will influence the risk of developing the disease.4 5
To date, the contribution of rare genetic variants and the impact of regulatory variants have not been widely explored in SLE. DNA sequencing has the potential to discover novel SLE associated variants not captured by genotyping arrays. Due to the high cost, whole genome sequencing studies (WGS) in SLE have so far mainly focused on families or smaller samples, as have exome sequencing studies (WES).6–9 Today it is feasible to perform targeted sequencing in larger cohorts; however, the number of such studies focusing on SLE is still limited.10 Additionally, association analysis for rare variants discovered through sequencing is hampered by low statistical power. Aggregating variants on the gene level or by molecular pathway information is one approach to increase power and gain biological insight from rare variants.11
The clinical heterogeneity in SLE is likely due to an underlying molecular diversity that could have implications for therapy. In recent years this has started to be addressed, mainly using gene expression, autoantibody profiles and cytokines to identify groups of patients with SLE with distinct molecular disease mechanisms.12–14 Using genetic information to stratify patients would have the advantage of providing stable molecular markers for early classification.
Here, we performed targeted sequencing of regulatory and coding regions in a Swedish SLE case–control cohort. We aimed at elucidating the genetic aetiology of SLE from the immunity pathway level to the single variant level, and stratify patients with SLE into molecular subgroups. Altogether around 9% of all genes in the human genome were analysed based on their role in immune-mediated diseases. Gene regions were extended to include promoters and other potentially regulatory elements based on mammalian conservation.15
Methods
For full details on methods see online supplemental methods.
Supplemental material
Subjects and DNA samples
The Swedish SLE cohorts included patients recruited at five rheumatology clinics and the controls were healthy blood donors and population controls. The quality-controlled dataset comprised 958 patients with SLE and 1026 control individuals. Patients with SLE fulfilled at least four of the classification criteria for SLE as defined by the American College of Rheumatology (ACR).16 17 Clinical characteristics of the patients are available in online supplemental tables S1A and B.
Supplemental material
Targeted DNA sequencing analysis
Targeted DNA sequencing was performed in the Swedish SLE case–control cohorts. A SeqCap EZ Choice XL sequence capture panel was designed, libraries were prepared as described elsewhere18 and sequenced on an Illumina HiSeq 2500. An overview of the variant discovery and quality control steps can be found in online supplemental figure S1. Study subjects falling outside of the European subpopulation of the Human Genome Diversity Project (HGDP) reference set were excluded (online supplemental figure S7).19 The quality-controlled dataset contained 287 354 single-nucleotide variants (SNVs) and covered 1832 of the targeted gene regions.
Supplemental material
Genetic association analyses
Several variant sets were generated for aggregate association testing: (1) 1832 individual gene variant sets; (2) 35 pathway variant sets based on the Kyoto Encyclopaedia of Genes and Genomes (KEGG)20 ; (3) five literature review gene sets: the type I interferon pathway,21 interferonopathy genes,22 23 SLE Genome-Wide Association Study (GWAS) genes,3 4 the complement subset of KEGG hsa04610 and genes causing monogenic SLE or lupus-like disease.24 Aggregate association testing was performed using Sequence Kernel Association Optimal Test (SKAT-O) or GenePy.25 26 Single variant association analyses were performed in PLINK. SLE case-only variants were identified by removing all SNVs present in our Swedish control dataset, the SweGen project or the Genome Aggregation Database European non-Finnish controls.27 28
Risk scores and cluster analysis
Cumulative pathway SLE polygenic risk scores (pathway PRSs) were assigned to each individual based on SNVs associated with SLE at nominal significance. For each independent SNV the natural logarithm of the OR for SLE susceptibility was multiplied by the number of minor alleles in each individual. The sum of all products of all genes in each of the 35 KEGG pathways for each patient was defined as the individual pathway PRS. Hierarchical cluster analysis of pathway PRSs was used to identify groups of patients with SLE.
Replication study and meta-analysis
Replication genotyping in individuals from Norway and Denmark was performed using the MassARRAY system. The Swedish SLE case–control study was expanded to include an additional 1000 control individuals.27 The Scandinavian meta-analysis included 1794 patients with SLE and 3241 control individuals.
Results
We performed a DNA sequencing study in SLE to study immunity pathways, an overview of analyses can be found in online supplemental figure S2.
T lymphocyte differentiation and innate immunity pathways are associated with SLE
The sequencing data analysis focused on 1832 genes with relevance for immune-mediated diseases. These genes mainly belong to 35 molecular signalling pathways as defined by the KEGG database (online supplemental table S2).20 Using an aggregate test for all variants in the genes belonging to each pathway, we found that 21 of the tested pathways were associated with SLE (false discovery rate (FDR) <0.05, table 1 and online supplemental table S3). The most significantly associated pathways included T helper cell differentiation pathways, with Th1 and Th2 cell differentiation as the top result (FDRTh1-2=2.2×10-9; FDRTh17=1.5×10-8), followed by antigen processing and presentation (FDR=3.1×10-9).
We next explored a sequential elimination strategy to identify independent pathway associations. First, removing all Th1 and Th2 pathway genes in the pathway aggregate association test resulted in the antigen processing and presentation pathway as the top result (FDR=4.8×10-6). Second, antigen processing and presentation as well as Th1 and Th2 pathway genes were removed, which resulted in Complement and coagulation cascades as the top result (FDR=0.0091). Third, also genes in this pathway were removed, and the janus kinase-signal transducers and activators of transcription (JAK-STAT) pathway became the top result (FDR=0.014). Lastly, when removing genes in all these four pathways no significant pathways remained. Thus, our data point to two main routes with genetic evidence of association to SLE: T cell differentiation and innate immunity.
To identify the genes that underlie the association signals in the T-cell differentiation, antigen processing and presentation, Complement and coagulation and JAK-STAT pathways, gene-based association testing was performed (figure 1). The top association for the JAK-STAT pathway originated from the IFN kappa (IFNK) gene region. SLE-associated genes in the T cell differentiation and antigen processing and presentation pathways were dominated by genes in the HLA region, and for the complement and coagulation cascade pathway, complement genes located in the HLA region were highly significantly associated with SLE.
Pathway PRS define subsets of patients with SLE
Having identified pathways with genetic association with SLE, we hypothesised that different patients with SLE could have distinct pathways affected. We constructed pathway PRS for each individual and each of the pathways, by combining the burden of common SLE associated alleles from our sequencing data. Individuals with a pathway PRS higher than that observed for the 97.5th percentile of control individuals were classified as positive for that pathway (online supplemental figure S3). The largest proportion of positive SLE patients was observed for the Cytokine-cytokine receptor interaction pathway (41%, figure 2A, and online supplemental table S4). For the Th1 and Th2 cell differentiation, antigen processing and presentation, Complement and coagulation cascades and JAK-STAT signalling pathways 18%, 16%, 21% and 29% of patients with SLE were positive, respectively. On average each SLE patient tested positive for the pathway PRS for seven pathways (figure 2B). As we had previously observed that a high SLE genetic risk score was associated with organ damage in SLE, we investigated whether this could be observed for specific pathways.5 We found that the SLE International Collaborating Clinics Damage Index was significantly higher in the SLE patients positive for the T cell or B cell receptor signalling pathways (figure 3A,B). No other pathways were associated with clinical manifestations of SLE or survival.
We then performed a hierarchical cluster analysis on the pathway PRSs in SLE, to identify groups of patients with similar molecular aetiology. Four clusters of patients were identified (figure 4). The pathway with the most significant difference in PRS between clusters was the antigen processing and presentation pathway, followed by Th17 cell differentiation (online supplemental figure S4). Next, we investigated whether the molecular stratification of patients with SLE also mirrored differences in clinical presentation between groups. We found that the presence of autoantibodies against Sjögren’s syndrome-related antigens SSA and/or SSB was more common among patients in clusters 3 and 4 (figure 3C). We did not observe any significant difference in other clinical features, including survival, between the four patient clusters.
Common variants contribute risk at monogenic risk loci in SLE
We then focused our analysis on gene-sets with prior evidence for involvement in SLE, but which were not defined in KEGG, to investigate the impact of both rare and common variants for these groups of genes. We found that interferon system, interferonopathy, SLE GWAS, complement system and monogenic SLE and lupus-like disease genes in aggregate were associated with SLE when analysing variants of all minor allele frequencies (MAF) (table 2 and online supplemental table S5). Only the monogenic SLE and lupus-like disease gene-set was significantly associated with SLE when separately analysing the rarer variant (MAF <0.01) contribution (table 2). There was a clear common variant (MAF >0.05) contribution to associations for the interferonopathy, SLE GWAS, complement system and monogenic SLE and lupus-like disease gene-sets (table 2).
Potentially novel SLE risk loci
Next, we asked whether we could detect novel SLE risk loci, regardless of pathway or gene-set membership. Two potentially novel gene regions passed a Bonferroni corrected threshold in the gene-based SLE case–control association analyses: PABPC4 (p=4.3×10-8) and IFNK (p=1.2×10-5, online supplemental figure 5A, tables S6 and S7). In single variant association analyses, we observed SNV associations at three potentially novel SLE risk loci, CAPN13, MOB3B/IFNK and HAL, at a suggestive significance threshold (p<1×10-4, online supplemental figure 5B–E, table S8). As the association signals at CAPN13, MOB3B/IFNK and HAL had not been reported in SLE GWAS in other ancestries, we attempted to replicate these findings in additional Scandinavian SLE cases and controls (online supplemental table S1A). However, we did not find additional support for a role of SNVs at these novel loci in SLE (online supplemental table S9).
Patients with SLE carry unique coding variants
We next investigated whether there was an increased rare coding mutational burden for patients with SLE at the 1832 genes. We observed that all individuals carried rare non-synonymous variants, with an average number of around 32 variants per individual for both patients with SLE and control individuals (online supplemental figure S6). None of the patients with SLE were homozygous carriers of rare non-synonymous alleles in genes for monogenic SLE and lupus-like diseases (online supplemental table S10). Next, we hypothesised that protein coding variants observed exclusively in patients with SLE could be causal candidates. A total of 1475 case-only nonsynonymous variants were identified in the 958 patients with SLE (online supplemental table S11). These were variants that were observed in at least one patient with SLE, but not in control individuals of similar ancestry.27 28 The most frequent of these SNVs was found in the MUC5B gene which encodes mucin 5B, the major gel-forming mucin in mucus (table 3). Five patients with SLE carried the same deleterious MUC5B missense mutation (rs773068050, p.Thr2724Pro). MUC5B gene variants have previously been associated with interstitial lung disease (ILD), a condition affecting around 3% of Swedish patients with SLE.29–31 However, there was no evidence of ILD in these five patients, but two of them had suffered from pleuritis (online supplemental table S12). In conclusion, we did not find evidence for SLE patients carrying a generally increased burden of rare coding variants at these genes. However, our analysis identified a number of coding variants observed exclusively in patients with SLE. This catalogue of variants could serve as a resource for future studies investigating the role of case-only SNVs in SLE.
Discussion
We here suggest a novel pathway-based approach to stratify patients with SLE beyond clinical manifestations. Further, we characterise genetic pathway associations and investigate rare variant contributions to the pathogenesis of SLE, all using targeted sequencing of immunity genes.
Using case–control association testing for immunological pathways, we identified two main axes of SLE association: T cell differentiation and innate immunity pathways. T cells have a fundamental role in loss of tolerance, autoimmunity and inflammatory reactions. In SLE, a number of different T cell disturbances have been described, which can contribute to the generation of autoreactive T cells, aberrant cytokine production and impaired T regulatory cell function.32 Besides the direct involvement of pathways connected to Th1 and Th2 cells, we noticed association signals from two pathways related to interleukin 17 (IL-17). A proportion of patients with SLE display raised serum levels of IL-17, elevated numbers of circulating IL-17-producing T cells and increased IL-17 production by lymphocytes, suggesting dysregulation of T regulatory cells.33 Our findings strengthen the recent suggestions that IL-17 inhibition could be a therapeutic option in a subset of patients with SLE.34 Conversely, low-dose IL-2 treatment in SLE to stimulate T regulatory cells has recently shown promising results.35
We observed that the T cell differentiation pathway associations were influenced by genetic associations to HLA, which is not surprising given the essential role of HLA in the immune response. This was further demonstrated by the antigen processing and presentation pathway association dominated by HLA genes. Complement pathway associations are also possibly confounded by the HLA SLE association, since early complement component genes are located in the HLA class III locus on chromosome 6.36 The JAK-STAT pathway was associated with SLE, it is the main route to initiate gene expression and protein synthesis for over 50 cytokines, many of which are involved in the SLE disease process.37 38 Variants of a number of genes in the JAK-STAT pathway have been associated with an increased risk for SLE, for instance STAT4-STAT1 and TYK2.3 4
Our study highlights the importance of the interferon system in SLE. Previous studies have shown genetic associations at a number of genes in the IFN signalling pathway in SLE.2 3 Here, we show that, in aggregate, genetic variation at interferonopathy genes also contribute to risk for SLE. In addition to interferonopathy genes, we also observed an aggregate genetic association for monogenic SLE and lupus-like disease genes with both a rare and a common variant contribution. This supports the hypothesis of a shared genetic basis and consequently disease mechanisms between monogenic and complex forms of disease, where also common non-coding variants can affect the regulation of Mendelian disease genes resulting in clinically similar traits.39
We have previously demonstrated that an SLE genetic risk score was associated with disease severity in SLE.5 We here generated a pathway-centred SLE PRS and found that there was a large variation in the number of affected pathways among the patients, which underscores the heterogeneity of SLE. We observed higher SLE damage indexes in patients with SLE positive for the B or T cell receptor signalling pathways, thus, pathways in the adaptive immune system seem important for the long-term severity of the disease. This is in accordance with previous findings that SLE disease activity correlates with abnormal B lymphocyte activity and T cell abnormalities, as well as the connection between disease activity and accumulation of organ damage.40 41
We attempted to cluster patients into subsets with shared genetic pathway profiles, which suggested four subgroups of patients with SLE. Beside the SSA/SSB antibody profile, these clusters were not connected to clinical disease manifestations such as nephritis or survival. This observation may indicate that the PRS reflects part of the central autoimmune process, which is not translated into specific organ manifestations. Whether the PRS in individual patients with SLE, or the different clusters, contribute to treatment response is an interesting possibility, but could not be assessed in this study. This is one limitation of our study, together with the fact that our conclusions apply specifically to this set of candidate genes.
WGS or WES studies will be required to fully elucidate the role of rare variants and pathways in SLE. As previously shown by us and others, WGS and WES in selected patients can provide information on ultrarare and de novo SNVs in SLE.6 7 42 However, larger sample sizes than those reported to date will be required to paint a complete picture of the genetic aetiology of SLE. We did not find support in additional Scandinavian cohorts for a role in SLE for the novel loci identified in the Swedish cohorts. Possible explanations include overestimated effect sizes in the discovery cohort, differences in genetic background within Scandinavia, or differences in clinical manifestations or characterisation of patients. Lastly, our study identified a large number of case-only coding variants. Variants uniquely identified in patients could be causal candidates in SLE, but their statistical significance is difficult to evaluate.
In summary, we have suggested a novel strategy to genetically stratify patients with SLE according to involved molecular pathways. T cell pathways displayed the strongest association, which highlights the importance of the adaptive immune system in the disease. The strong connection to the JAK-STAT pathway, including the IFN system, is perhaps not surprising given the promising clinical trials of JAK and type I interferon receptor inhibition as treatments for SLE.38 43 44 However, not all patients in these studies respond to treatment, and dissecting affected molecular pathways in responders and non-responders could increase the understanding of treatment outcome. This approach has not been tested clinically, but the future of precision medicine for SLE lies in identifying robust methods to perform molecular stratification of patients.
Acknowledgments
DNA sequencing and genotyping was performed at the SNP&SEQ Technology Platform in Uppsala. The facility is part of the National Genomics Infrastructure (NGI) Sweden and Science for Life Laboratory. Computations were performed on resources provided by the Swedish National Infrastructure for Computing (SNIC) through Uppsala Multidisciplinary Centre for Advanced Computational Science (UPPMAX) under projects SNIC SENS 2017142 and 2017107. The SweGen genotype data were generated by Science for Life Laboratory. The authors would like to thank the Genome Aggregation Database (gnomAD) and the groups that provided exome and genome variant data to this resource. A full list of contributing groups can be found at the gnomAD website. The authors wish to thank the Uppsala Bioresource Karolina Tandre and Västerbotten biobank for providing DNA samples on control individuals, and Cane Amcoff for excellent technical assistance.
Supplementary materials
Supplementary Data
This web only file has been produced by the BMJ Publishing Group from an electronic file supplied by the author(s) and has not been edited for content.
Footnotes
Handling editor Josef S Smolen
Collaborators The DISSECT consortium: Johanna K. Sandling (Department of Medical Sciences, Rheumatology, Uppsala University, Uppsala, Sweden), Pascal Pucholt (Department of Medical Sciences, Rheumatology, Uppsala University, Sweden), Lina Hultin Rosenberg (Science for Life Laboratory, Department of Medical Biochemistry and Microbiology, Uppsala University, Uppsala, Sweden), Fabiana H.G. Farias (Science for Life Laboratory, Department of Medical Biochemistry and Microbiology, Uppsala University, Uppsala, Sweden, and Department of Psychiatry, Washington University, St. Louis, MO, USA), Sergey V. Kozyrev (Science for Life Laboratory, Department of Medical Biochemistry and Microbiology, Uppsala University, Uppsala, Sweden), Maija-Leena Eloranta (Department of Medical Sciences, Rheumatology, Uppsala University, Uppsala, Sweden), Andrei Alexsson (Department of Medical Sciences, Rheumatology, Uppsala University, Uppsala, Sweden), Matteo Bianchi (Science for Life Laboratory, Department of Medical Biochemistry and Microbiology, Uppsala University, Uppsala, Sweden), Leonid Padyukov (Division of Rheumatology, Department of Medicine, Karolinska Institutet and Karolinska University Hospital, Stockholm, Sweden), Christine Bengtsson (Department of Public Health and Clinical Medicine/Rheumatology, Umeå University, Umeå, Sweden), Roland Jonsson (Broegelmann Research Laboratory, Department of Clinical Science, University of Bergen, Bergen, Norway), Roald Omdal (Clinical Immunology unit, Department of Internal Medicine, Stavanger University Hospital, Stavanger, Norway and Broegelmann Research Laboratory, Department of Clinical Science, University of Bergen, Bergen, Norway), Øyvind Molberg (Department of Rheumatology, Oslo University Hospital and Institute of Clinical Medicine, University of Oslo, Oslo, Norway), Ann-Christine Syvänen (Department of Medical Sciences, Molecular Medicine and Science for Life Laboratory, Uppsala University, Uppsala, Sweden), Andreas Jönsen (Lund University, Skane University Hospital, Department of Clinical Sciences Lund, Rheumatology, Lund, Sweden), Iva Gunnarsson (Division of Rheumatology, Department of Medicine Solna, Karolinska Institutet, Karolinska University Hospital, Stockholm, Sweden), Elisabet Svenungsson (Division of Rheumatology, Department of Medicine Solna, Karolinska Institutet, Karolinska University Hospital, Stockholm, Sweden), Solbritt Rantapää-Dahlqvist (Department of Public Health and Clinical Medicine/Rheumatology, Umeå University, Umeå, Sweden), Anders A. Bengtsson (Lund University, Skane University Hospital, Department of Clinical Sciences Lund, Rheumatology, Lund, Sweden), Christopher Sjöwall (Department of Biomedical and Clinical Sciences, Division of Inflammation and Infection, Linköping University, Linköping, Sweden), Dag Leonard (Department of Medical Sciences, Rheumatology, Uppsala University, Uppsala, Sweden), Kerstin Lindblad-Toh (Science for Life Laboratory, Department of Medical Biochemistry and Microbiology, Uppsala University, Uppsala, Sweden and Broad Institute of MIT and Harvard, Cambridge, MA, USA), Lars Rönnblom (Department of Medical Sciences, Rheumatology, Uppsala University, Uppsala, Sweden), Jonas Carlsson Almlöf (Department of Medical Sciences, Molecular Medicine and Science for Life Laboratory, Uppsala University, Uppsala, Sweden), Johanna Dahlqvist (Science for Life Laboratory, Department of Medical Sciences and Department of Medical Biochemistry and Microbiology, Uppsala University, Uppsala, Sweden), Daniel Eriksson (Department of Medicine (Solna), Karolinska Institutet, and Department of Endocrinology, Metabolism and Diabetes Karolinska University Hospital, Stockholm, Sweden), Niklas Hagberg (Department of Medical Sciences, Rheumatology, Uppsala University, Uppsala, Sweden), Ingrid E. Lundberg (Division of Rheumatology, Department of Medicine and Center for Molecular Medicine, Karolinska Institutet, Stockholm, Sweden), Argyri Mathioudaki (Science for Life Laboratory, Department of Medical Biochemistry and Microbiology, Uppsala University, Uppsala, Sweden), Jennifer Meadows (Science for Life Laboratory, Department of Medical Biochemistry and Microbiology, Uppsala University, Uppsala, Sweden), Jessika Nordin (Science for Life Laboratory, Department of Medical Biochemistry and Microbiology, Uppsala University, Uppsala, Sweden), Gunnel Nordmark (Department of Medical Sciences, Rheumatology, Uppsala University, Uppsala, Sweden), Marie Wahren-Herlenius (Department of Medicine, Division of Rheumatology, Karolinska Institutet, Karolinska University Hospital, Stockholm, Sweden and Broegelmann Research Laboratory, Department of Clinical Science, University of Bergen, Norway), Sule Yavuz (Department of Medical Sciences, Rheumatology, Uppsala University, Uppsala, Sweden). The ImmunoArray development consortium: Kerstin Lindblad-Toh (Science for Life Laboratory, Department of Medical Biochemistry and Microbiology, Uppsala University, Uppsala, Sweden and Broad Institute of MIT and Harvard, Cambridge, MA, USA), Gerli Rosengren Pielberg (Science for Life Laboratory, Department of Medical Biochemistry and Microbiology, Uppsala University, Uppsala, Sweden), Anna Lobell (Office for Medicine and Pharmacy, Uppsala University, Uppsala, Sweden), Åsa Karlsson (Science for Life Laboratory, Department of Medical Biochemistry and Microbiology, Uppsala University, Uppsala, Sweden), Eva Murén (Science for Life Laboratory, Department of Medical Biochemistry and Microbiology, Uppsala University, Uppsala, Sweden), Göran Andersson (Department of Animal Breeding and Genetics, Swedish University of Agricultural Sciences, Uppsala, Sweden), Kerstin M. Ahlgren (Department of Surgical Sciences, Uppsala University, Uppsala, Sweden), Lars Rönnblom (Department of Medical Sciences, Rheumatology, Uppsala University, Uppsala, Sweden), Maija-Leena Eloranta (Department of Medical Sciences, Rheumatology, Uppsala University, Uppsala, Sweden), Nils Landegren (Department of Medicine (Solna), Center for Molecular Medicine, Karolinska Institutet, Stockholm, Sweden and Science for Life Laboratory, Department of Medical Sciences, Uppsala University, Uppsala, Sweden), Olle Kämpe (Department of Medicine (Solna), Center for Molecular Medicine, Karolinska Institutet, Stockholm, Sweden, Department of Endocrinology, Metabolism and Diabetes Karolinska University Hospital, Stockholm, Sweden, Science for Life Laboratory, Department of Medical Sciences, Uppsala University, Uppsala, Sweden and KG Jebsen Center for autoimmune diseases, University of Bergen, Norway), Peter Söderkvist (Division of Cell Biology, Department of Biomedical and Clinical Sciences, Linköping University, Linköping, Sweden).
Contributors KLT and LR conceived and designed the experiments. CB, KL, AV, AMT, SJ, AJ, IG, ES, SR-D, AAB, CS and DL characterised the patient samples. M-LE, LP, LM, RS, MAJ, RJ, RO, BAL and STL provided samples and data. The ImmunoArray development consortium members designed the targeted sequencing panel. FHGF, SVK, ÅK, and EM performed the experiments. A-CS managed the sequencing platform. JKS, LHR, PP, AA and MB analysed the data and the DISSECT consortium members provided intellectual input and/or developed analysis pipelines. JKS and LR wrote the manuscript. All authors read, provided critical review and accepted the final version of the manuscript.
Funding This study was supported by an AstraZeneca-Science for Life Laboratory Research Collaboration grant (DISSECT). This study was also supported by grants from the Swedish Research Council for Medicine and Health (Dnr 2018–02399, 2018–02535 and a Distinguished professor award to KL-T), the Swedish Rheumatism Association, King Gustav V’s 80-year Foundation, the Swedish-Heart-Lung foundation, a Wallenberg Scholar Award (to KL-T), and the Swedish Society of Medicine and the Ingegerd Johansson donation. The SNP&SEQ Platform is supported by Science for Life Laboratory, the Swedish Research Council (VR-RFI), Uppsala University and the Knut and Alice Wallenberg Foundation.
Disclaimer The funders had no role in the design of the study or collection, analysis, or interpretation of data or in writing the manuscript.
Competing interests An AstraZeneca research collaboration grant to LR supported the study. LHR is now an employee of Olink Proteomics.
Patient consent for publication Not required.
Ethics approval The study was approved by the regional ethics board in Uppsala (Dnr 2015/450 and 2016/155).
Provenance and peer review Not commissioned; externally peer reviewed.
Data availability statement The datasets generated and/or analysed during the current study are not publicly available due to them containing information that could compromise research participant privacy and consent, but are available from the corresponding authors LR (ORCID 0000-0001-9403-6503) and JKS (ORCID 0000-0003-1382-2321) on reasonable request and on a collaborative basis.
Supplemental material This content has been supplied by the author(s). It has not been vetted by BMJ Publishing Group Limited (BMJ) and may not have been peer-reviewed. Any opinions or recommendations discussed are solely those of the author(s) and are not endorsed by BMJ. BMJ disclaims all liability and responsibility arising from any reliance placed on the content. Where the content includes any translated material, BMJ does not warrant the accuracy and reliability of the translations (including but not limited to local regulations, clinical guidelines, terminology, drug names and drug dosages), and is not responsible for any error and/or omissions arising from translation and adaptation or otherwise.