Identification of Mechanisms by Which Genetic Susceptibility Loci Influence Systemic Sclerosis Risk Using Functional Genomics in Primary T Cells and Monocytes.

OBJECTIVE: Systemic sclerosis (SSc) is a complex autoimmune disease with a strong genetic component. However, most of the genes associated with the disease are still unknown because associated variants affect mostly noncoding intergenic elements of the genome. We used functional genomics to translate the genetic findings into a better understanding of the disease. METHODS: Promoter capture Hi-C and RNA-sequencing experiments were performed in CD4+ T cells and CD14+ monocytes from 10 SSc patients and 5 healthy controls to link SSc-associated variants with their target genes, followed by differential expression and differential interaction analyses between cell types. RESULTS: We linked SSc-associated loci to 39 new potential target genes and confirmed 7 previously known SSc-associated genes. We highlight novel causal genes, such as CXCR5, as the most probable candidate gene for the DDX6 locus. Some previously known SSc-associated genes, such as IRF8, STAT4, and CD247, showed cell type-specific interactions. We also identified 15 potential drug targets already in use in other similar immune-mediated diseases that could be repurposed for SSc treatment. Furthermore, we observed that interactions were directly correlated with the expression of important genes implicated in cell type-specific pathways and found evidence that chromatin conformation is associated with genotype. CONCLUSION: Our study revealed potential causal genes for SSc-associated loci, some of them acting in a cell type-specific manner, suggesting novel biologic mechanisms that might mediate SSc pathogenesis.

Expression Quantitative Trait Locus Analysis in Systemic Sclerosis Identifies New Candidate Genes Associated With Multiple Aspects of Disease Pathology.

OBJECTIVE: To identify the genetic variants that affect gene expression (expression quantitative trait loci [eQTLs]) in systemic sclerosis (SSc) and to investigate their role in the pathogenesis of the disease. METHODS: We performed an eQTL analysis using whole-blood sequencing data from 333 SSc patients and 524 controls and integrated them with SSc genome-wide association study (GWAS) data. We integrated our findings from expression modeling, differential expression analysis, and transcription factor binding site enrichment with key clinical features of SSc. RESULTS: We detected 49,123 validated cis-eQTLs from 4,539 SSc-associated single-nucleotide polymorphisms (SNPs) (PGWAS < 10-5 ). A total of 1,436 genes were within 1 Mb of the 4,539 SSc-associated SNPs. Of those 1,436 genes, 565 were detected as having ≥1 eQTL with an SSc-associated SNP. We developed a strategy to prioritize disease-associated genes based on their expression variance explained by SSc eQTLs (r2 > 0.05). As a result, 233 candidates were identified, 134 (58%) of them associated with hallmarks of SSc and 105 (45%) of them differentially expressed in the blood cells, skin, or lung tissue of SSc patients. Transcription factor binding site analysis revealed enriched motifs of 24 transcription factors (5%) among SSc eQTLs, 5 of which were found to be differentially regulated in the blood cells (ELF1 and MGA), skin (KLF4 and ID4), and lungs (TBX4) of SSc patients. Ten candidate genes (4%) can be targeted by approved medications for immune-mediated diseases, of which only 3 have been tested in clinical trials in patients with SSc. CONCLUSION: The findings of the present study indicate a new layer to the molecular complexity of SSc, contributing to a better understanding of the pathogenesis of the disease.

Approaching Shared Pathophysiology in Immune-Mediated Diseases through Functional Genomics.

Immune-mediated diseases (IMDs) are complex pathologies that are strongly influenced by environmental and genetic factors. Associations between genetic loci and susceptibility to these diseases have been widely studied, and hundreds of risk variants have emerged during the last two decades, with researchers observing a shared genetic pattern among them. Nevertheless, the pathological mechanism behind these associations remains a challenge that has just started to be understood thanks to functional genomic approaches. Transcriptomics, regulatory elements, chromatin interactome, as well as the experimental characterization of genomic findings, constitute key elements in the emerging understandings of how genetics affects the etiopathogenesis of IMDs. In this review, we will focus on the latest advances in the field of functional genomics, centering our attention on systemic rheumatic IMDs.

Analysis of Anti-RNA Polymerase III Antibody-positive Systemic Sclerosis and Altered GPATCH2L and CTNND2 Expression in Scleroderma Renal Crisis.

OBJECTIVE: Scleroderma renal crisis (SRC) is a life-threatening complication of systemic sclerosis (SSc) strongly associated with anti-RNA polymerase III antibody (ARA) autoantibodies. We investigated genetic susceptibility and altered protein expression in renal biopsy specimens in ARA-positive patients with SRC. METHODS: ARA-positive patients (n = 99) with at least 5 years' follow-up (49% with a history of SRC) were selected from a well characterized SSc cohort (n = 2254). Cases were genotyped using the Illumina Human Omni-express chip. Based on initial regression analysis, 9 single-nucleotide polymorphisms (SNP) were chosen for validation in a separate cohort of 256 ARA-positive patients (40 with SRC). Immunostaining of tissue sections from SRC or control kidney was used to quantify expression of candidate proteins based upon genetic analysis of the discovery cohort. RESULTS: Analysis of 641,489 SNP suggested association of POU2F1 (rs2093658; P = 1.98 × 10-5), CTNND2 (rs1859082; P = 5.58 × 10-5), HECW2 (rs16849716; P = 1.2 × 10-4), and GPATCH2L (rs935332; P = 4.92 × 10-5) with SRC. Further, the validation cohort showed an association between rs935332 within the GPATCH2L region, with SRC (P = 0.025). Immunostaining of renal biopsy sections showed increased tubular expression of GPATCH2L (P = 0.026) and glomerular expression of CTNND2 (P = 0.026) in SRC samples (n = 8) compared with normal human kidney controls (n = 8), despite absence of any genetic replication for the associated SNP. CONCLUSION: Increased expression of 2 candidate proteins, GPATCH2L and CTNND2, in SRC compared with control kidney suggests a potential role in pathogenesis of SRC. For GPATCH2L, this may reflect genetic susceptibility in ARA-positive patients with SSc based upon 2 independent cohorts.

A cross-disease meta-GWAS identifies four new susceptibility loci shared between systemic sclerosis and Crohn's disease.

Genome-wide association studies (GWASs) have identified a number of genetic risk loci associated with systemic sclerosis (SSc) and Crohn's disease (CD), some of which confer susceptibility to both diseases. In order to identify new risk loci shared between these two immune-mediated disorders, we performed a cross-disease meta-analysis including GWAS data from 5,734 SSc patients, 4,588 CD patients and 14,568 controls of European origin. We identified 4 new loci shared between SSc and CD, IL12RB2, IRF1/SLC22A5, STAT3 and an intergenic locus at 6p21.31. Pleiotropic variants within these loci showed opposite allelic effects in the two analysed diseases and all of them showed a significant effect on gene expression. In addition, an enrichment in the IL-12 family and type I interferon signaling pathways was observed among the set of SSc-CD common genetic risk loci. In conclusion, through the first cross-disease meta-analysis of SSc and CD, we identified genetic variants with pleiotropic effects on two clinically distinct immune-mediated disorders. The fact that all these pleiotropic SNPs have opposite allelic effects in SSc and CD reveals the complexity of the molecular mechanisms by which polymorphisms affect diseases.

The Potential Role of Genomic Medicine in the Therapeutic Management of Rheumatoid Arthritis.

During the last decade, important advances have occurred regarding understanding of the pathogenesis and treatment of rheumatoid arthritis (RA). Nevertheless, response to treatment is not universal, and choosing among different therapies is currently based on a trial and error approach. The specific patient's genetic background influences the response to therapy for many drugs: In this sense, genomic studies on RA have produced promising insights that could help us find an effective therapy for each patient. On the other hand, despite the great knowledge generated regarding the genetics of RA, most of the investigations performed to date have focused on identifying common variants associated with RA, which cannot explain the complete heritability of the disease. In this regard, rare variants could also contribute to this missing heritability as well as act as biomarkers that help in choosing the right therapy. In the present article, different aspects of genetics in the pathogenesis and treatment of RA are reviewed, from large-scale genomic studies to specific rare variant analyses. We also discuss the shared genetic architecture existing among autoimmune diseases and its implications for RA therapy, such as drug repositioning.

Analysis of the genetic component of systemic sclerosis in Iranian and Turkish populations through a genome-wide association study.

Objectives: SSc is an autoimmune disease characterized by alteration of the immune response, vasculopathy and fibrosis. Most genetic studies on SSc have been performed in European-ancestry populations. The aim of this study was to analyse the genetic component of SSc in Middle Eastern patients from Iran and Turkey through a genome-wide association study. Methods: This study analysed data from a total of 834 patients diagnosed with SSc and 1455 healthy controls from Iran and Turkey. DNA was genotyped using high-throughput genotyping platforms. The data generated were imputed using the Michigan Imputation Server, and the Haplotype Reference Consortium as a reference panel. A meta-analysis combining both case-control sets was conducted by the inverse variance method. Results: The highest peak of association belonged to the HLA region in both the Iranian and Turkish populations. Strong and independent associations between the classical alleles HLA-DRB1*11: 04 [P = 2.10 × 10-24, odds ratio (OR) = 3.14] and DPB1*13: 01 (P = 5.37 × 10-14, OR = 5.75) and SSc were observed in the Iranian population. HLA-DRB1*11: 04 (P = 4.90 × 10-11, OR = 2.93) was the only independent signal associated in the Turkish cohort. An omnibus test yielded HLA-DRB1 58 and HLA-DPB1 76 as relevant amino acid positions for this disease. Concerning the meta-analysis, we also identified two associations close to the genome-wide significance level outside the HLA region, corresponding to IRF5-TNPO3 rs17424921-C (P = 1.34 × 10-7, OR = 1.68) and NFKB1 rs4648133-C (P = 3.11 × 10-7, OR = 1.47). Conclusion: We identified significant associations in the HLA region and suggestive associations in IRF5-TNPO3 and NFKB1 loci in Iranian and Turkish patients affected by SSc through a genome-wide association study and an extensive HLA analysis.

Genome-wide meta-analysis reveals shared new loci in systemic seropositive rheumatic diseases.

OBJECTIVE: Immune-mediated inflammatory diseases (IMIDs) are heterogeneous and complex conditions with overlapping clinical symptoms and elevated familial aggregation, which suggests the existence of a shared genetic component. In order to identify this genetic background in a systematic fashion, we performed the first cross-disease genome-wide meta-analysis in systemic seropositive rheumatic diseases, namely, systemic sclerosis, systemic lupus erythematosus, rheumatoid arthritis and idiopathic inflammatory myopathies. METHODS: We meta-analysed ~6.5 million single nucleotide polymorphisms in 11 678 cases and 19 704 non-affected controls of European descent populations. The functional roles of the associated variants were interrogated using publicly available databases. RESULTS: Our analysis revealed five shared genome-wide significant independent loci that had not been previously associated with these diseases: NAB1, KPNA4-ARL14, DGQK, LIMK1 and PRR12. All of these loci are related with immune processes such as interferon and epidermal growth factor signalling, response to methotrexate, cytoskeleton dynamics and coagulation cascade. Remarkably, several of the associated loci are known key players in autoimmunity, which supports the validity of our results. All the associated variants showed significant functional enrichment in DNase hypersensitivity sites, chromatin states and histone marks in relevant immune cells, including shared expression quantitative trait loci. Additionally, our results were significantly enriched in drugs that are being tested for the treatment of the diseases under study. CONCLUSIONS: We have identified shared new risk loci with functional value across diseases and pinpoint new potential candidate loci that could be further investigated. Our results highlight the potential of drug repositioning among related systemic seropositive rheumatic IMIDs.

A TNFSF13B functional variant is not involved in systemic sclerosis and giant cell arteritis susceptibility.

BACKGROUND: The TNFSF13B (TNF superfamily member 13b) gene encodes BAFF, a cytokine with a crucial role in the differentiation and activation of B cells. An insertion-deletion variant (GCTGT→A) of this gene, leading to increased levels of BAFF, has been recently implicated in the genetic predisposition to several autoimmune diseases, including multiple sclerosis, systemic lupus erythematosus, and rheumatoid arthritis. Based on the elevated levels of this cytokine found in patients with giant cell arteritis (GCA) and systemic sclerosis (SSc), we aimed to assess whether this functional variant also represents a novel genetic risk factor for these two disorders. METHODS: A total of 1,728 biopsy-proven GCA patients from 4 European cohorts, 4,584 SSc patients from 3 European cohorts and 5,160 ethnically-matched healthy controls were included in the study. The single nucleotide polymorphism (SNP) rs374039502, which colocalizes with the genetic variant previously implicated in autoimmunity, was genotyped using a custom TaqMan assay. First, association analysis was conducted in each independent cohort using χ2 test in Plink (v1.9). Subsequently, different case/control sets were meta-analyzed by the inverse variance method. RESULTS: No statistically significant differences were found when allele distributions were compared between cases and controls for any of the analyzed cohorts. Similarly, combined analysis of the different sets evidenced a lack of association of the rs374039502 variant with GCA (P = 0.421; OR (95% CI) = 0.92 (0.75-1.13)) and SSc (P = 0.500; OR (95% CI) = 1.05 (0.91-1.22)). The stratified analysis considering the main clinical subphenotypes of these diseases yielded similar negative results. CONCLUSION: Our data suggest that the TNFSF13B functional variant does not contribute to the genetic network underlying GCA and SSc.

Association of a rare variant of the TNFSF13B gene with susceptibility to Rheumatoid Arthritis and Systemic Lupus Erythematosus.

A rare variant (BAFF-var) of the tumor necrosis factor superfamily 13b (TNFSF13B) gene has been recently associated with multiple sclerosis (MS) and systemic lupus erythematosus (SLE). The aim of this study was to investigate the association between TNFSF13B BAFF-var and susceptibility to rheumatoid arthritis (RA) and replicate that association in SLE. 6,218 RA patients, 2,575 SLE patients and 4,403 healthy controls from three different countries were included in the study. TNFSF13B BAFF-var was genotyped using TaqMan allelic discrimination assay. PLINK software was used for statistical analyses. TNFSF13B BAFF-var was significantly associated with RA (p = 0.015, OR = 1.21, 95% CI = 1.03-1.41) in the Spanish cohort. A trend of association was observed in the Dutch (p = 0.115) and German (p = 0.228) RA cohorts. A meta-analysis of the three RA cohorts included in this study revealed a statistically significant association (p = 0.002, OR = 1.24, 95% CI = 1.10-1.38). In addition, TNFSF13B BAFF-var was significantly associated with SLE in the Spanish (p = 0.001, OR = 1.41, 95% CI = 1.14-1.74) and the German cohorts (p = 0.030, OR = 1.86, 95% CI = 1.05-3.28), with a statistically significant p-value obtained in the meta-analysis (p = 0.0002, OR = 1.46, 95% CI = 1.09-2.32). The results obtained confirm the known association of TNFSF13B BAFF-var with SLE and, for the first time, demonstrate that this variant contributes to susceptibility to RA.