In this study, 1,833 systemic sclerosis (SSc) cases and 3,466 controls were genotyped with the Immunochip array. Classical alleles, amino acid residues, and SNPs across the human leukocyte antigen (HLA) region were imputed and tested. These analyses resulted in a model composed of six polymorphic amino acid positions and seven SNPs that explained the observed significant associations in the region. In addition, a replication step comprising 4,017 SSc cases and 5,935 controls was carried out for several selected non-HLA variants, reaching a total of 5,850 cases and 9,401 controls of European ancestry. Following this strategy, we identified and validated three SSc risk loci, including DNASE1L3 at 3p14, the SCHIP1-IL12A locus at 3q25, and ATG5 at 6q21, as well as a suggested association of the TREH-DDX6 locus at 11q23. The associations of several previously reported SSc risk loci were validated and further refined, and the observed peak of association in PXK was related to DNASE1L3. Our study has increased the number of known genetic associations with SSc, provided further insight into the pleiotropic effects of shared autoimmune risk factors, and highlighted the power of dense mapping for detecting previously overlooked susceptibility loci.

Systemic lupus erythematosus (SLE) is an autoimmune disease with a strong genetic component and is characterized by chronic inflammation and the production of anti-nuclear auto-antibodies. In the era of genome-wide association studies (GWASs), elucidating the genetic factors present in SLE has been a very successful endeavor; 28 confirmed disease susceptibility loci have been mapped. In this review, we summarize the current understanding of the genetics of lupus and focus on the strongest associated risk loci found to date (P <1.0 × 10−8). Although these loci account for less than 10% of the genetic heritability and therefore do not account for the bulk of the disease heritability, they do implicate important pathways, which contribute to SLE pathogenesis. Consequently, the main focus of the review is to outline the genetic variants in the known associated loci and then to explore the potential functional consequences of the associated variants. We also highlight the genetic overlap of these loci with other autoimmune diseases, which indicates common pathogenic mechanisms. The importance of developing functional assays will be discussed and each of them will be instrumental in furthering our understanding of these associated variants and loci. Finally, we indicate that performing a larger SLE GWAS and applying a more targeted set of methods, such as the ImmunoChip and next generation sequencing methodology, are important for identifying additional loci and enhancing our understanding of the pathogenesis of SLE.

We previously established an 80 kb haplotype upstream of TNFSF4 as a susceptibility locus in the autoimmune disease SLE. SLE-associated alleles at this locus are associated with inflammatory disorders, including atherosclerosis and ischaemic stroke. In Europeans, the TNFSF4 causal variants have remained elusive due to strong linkage disequilibrium exhibited by alleles spanning the region. Using a trans-ancestral approach to fine-map the locus, utilising 17,900 SLE and control subjects including Amerindian/Hispanics (1348 cases, 717 controls), African-Americans (AA) (1529, 2048) and better powered cohorts of Europeans and East Asians, we find strong association of risk alleles in all ethnicities; the AA association replicates in African-American Gullah (152,122). The best evidence of association comes from two adjacent markers: rs2205960-T (P = 1.71×10−34, OR = 1.43[1.26–1.60]) and rs1234317-T (P = 1.16×10−28, OR = 1.38[1.24–1.54]). Inference of fine-scale recombination rates for all populations tested finds the 80 kb risk and non-risk haplotypes in all except African-Americans. In this population the decay of recombination equates to an 11 kb risk haplotype, anchored in the 5′ region proximal to TNFSF4 and tagged by rs2205960-T after 1000 Genomes phase 1 (v3) imputation. Conditional regression analyses delineate the 5′ risk signal to rs2205960-T and the independent non-risk signal to rs1234314-C. Our case-only and SLE-control cohorts demonstrate robust association of rs2205960-T with autoantibody production. The rs2205960-T is predicted to form part of a decameric motif which binds NF-κBp65 with increased affinity compared to rs2205960-G. ChIP-seq data also indicate NF-κB interaction with the DNA sequence at this position in LCL cells. Our research suggests association of rs2205960-T with SLE across multiple groups and an independent non-risk signal at rs1234314-C. rs2205960-T is associated with autoantibody production and lymphopenia. Our data confirm a global signal at TNFSF4 and a role for the expressed product at multiple stages of lymphocyte dysregulation during SLE pathogenesis. We confirm the validity of trans-ancestral mapping in a complex trait.

No abstract

IntroductionClinical diversity in systemic sclerosis (SSc) reflects multifaceted pathogenesis and the effect of key growth factors or cytokines operating within a disease-specific microenvironment. Dermal interstitial fluid sampling offers the potential to examine local mechanisms and identify proteins expressed within lesional tissue. We used multiplex cytokine analysis to profile the inflammatory and immune activity in the lesions of SSc patients.MethodsDermal interstitial fluid sample from the involved forearm skin, and synchronous plasma samples were collected from SSc patients (n = 26, diffuse cutaneous SSc (DcSSc) n = 20, limited cutaneous SSc (LcSSc) n = 6), and healthy controls (HC) (n = 10) and profiled by Luminex® array for inflammatory cytokines, chemokines, and growth factors.ResultsLuminex® profiling of the dermal blister fluid showed increased inflammatory cytokines (median interleukin ( IL)-6 in SSc 39.78 pg/ml, HC 5.51 pg/ml, p = 0.01, median IL-15 in SSc 6.27 pg/ml, HC 4.38 pg/ml, p = 0.03), chemokines (monocyte chemotactic protein (MCP)-3 9.81 pg/ml in SSc, 7.18 pg/ml HC, p = 0.04), and profibrotic growth factors (platelet derived growth factor (PDGF)-AA 10.38 pg/ml versus 6.94 pg/ml in HC, p = 0.03). In general dermal fluid and plasma cytokine levels did not correlate, consistent with predominantly local production of these factors within the dermal lesions, rather than leakage from the serum. In hierarchical clustering and network analysis IL-6 emerged as a key central mediator.ConclusionsOur data confirm that an immuno-inflammatory environment and aberrant vascular repair are intimately linked to fibroblast activation in lesional skin in SSc. This non-invasive method could be used to profile disease activity in the clinic, and identifies key inflammatory or pro-fibrotic proteins that might be targeted therapeutically. Distinct subgroups of SSc may be defined that show innate or adaptive immune cytokine signatures.

No abstract

Extended molecular characterization of HLA genes in the IHWG reference B-lymphoblastoid cell lines (B-LCLs) was one of the major goals for the 17th International HLA and Immunogenetics Workshop (IHIW). Although reference B-LCLs have been examined extensively in previous workshops complete high-resolution typing was not completed for all the classical class I and class II HLA genes. To address this, we conducted a single-blind study where select panels of B-LCL genomic DNA samples were distributed to multiple laboratories for HLA genotyping by next-generation sequencing methods. Identical cell panels comprised of 24 and 346 samples were distributed and typed by at least four laboratories in order to derive accurate consensus HLA genotypes. Overall concordance rates calculated at both 2- and 4-field allele-level resolutions ranged from 90.4% to 100%. Concordance for the class I genes ranged from 91.7 to 100%, whereas concordance for class II genes was variable; the lowest observed at HLA-DRB3 (84.2%). At the maximum allele-resolution 78 B-LCLs were defined as homozygous for all 11 loci. We identified 11 novel exon polymorphisms in the entire cell panel. A comparison of the B-LCLs NGS HLA genotypes with the HLA genotypes catalogued in the IPD-IMGT/HLA Database Cell Repository, revealed an overall allele match at 68.4%. Typing discrepancies between the two datasets were mostly due to the lower-resolution historical typing methods resulting in incomplete HLA genotypes for some samples listed in the IPD-IMGT/HLA Database Cell Repository. Our approach of multiple-laboratory NGS HLA typing of the B-LCLs has provided accurate genotyping data. The data generated by the tremendous collaborative efforts of the 17th IHIW participants is useful for updating the current cell and sequence databases and will be a valuable resource for future studies.

Objective NKX2‐5 is a homeobox transcription factor that is required for the formation of the heart and vessels during development, with significant postnatal down‐regulation and reactivation in disease states, characterized by vascular remodeling. The purpose of this study was to investigate mechanisms that activate NKX2‐5 expression in diseased vessels, such as systemic sclerosis (scleroderma; SSc)–associated pulmonary hypertension (PH), and to identify genetic variability that potentially underlies susceptibility to specific vascular complications.MethodsWe explored NKX2‐5 expression in biopsy samples from patients with SSc‐associated PH and in pulmonary artery smooth muscle cells (PASMCs) from patients with scleroderma. Disease‐associated putative functional single‐nucleotide polymorphisms (SNPs) at the NKX2-5 locus were cloned and studied in reporter gene assays. SNP function was further examined through protein–DNA binding assays, chromatin immunoprecipitation assays, and RNA silencing analyses.ResultsIncreased NKX2‐5 expression in biopsy samples from patients with SSc‐associated PH was localized to remodeled vessels and PASMCs. Meta‐analysis of 2 independent scleroderma cohorts revealed an association of rs3131917 with scleroderma (P = 0.029). We demonstrated that disease‐associated SNPs are located in a novel functional enhancer, which increases NKX2-5 transcriptional activity through the binding of GATA‐6, c‐Jun, and myocyte‐specific enhancer factor 2C. We also characterized an activator/coactivator transcription‐enhancer factor domain 1 (TEAD1)/Yes‐associated protein 1 (YAP1) complex, which was bound at rs3095870, another functional SNP, with TEAD1 binding the risk allele and activating the transcription of NKX2-5.Conclusion NKX2-5 is genetically associated with scleroderma, pulmonary hypertension, and fibrosis. Functional evidence revealed a regulatory mechanism that results in NKX2-5 transcriptional activation in PASMCs through the interaction of an upstream promoter and a novel downstream enhancer. This mechanism can act as a model for NKX2‐5 activation in cardiovascular disease characterized by vascular remodeling.