Dysregulated long non-coding RNA in Sjögren's disease impacts both interferon and adaptive immune responses

OBJECTIVE

Sjögren's disease (SjD) is an autoimmune disease characterised by inflammatory destruction of exocrine glands. Patients with autoantibodies to Ro/SSA (SjDRo+) exhibit more severe disease. Long non-coding RNAs (lncRNAs) are a functionally diverse class of non-protein-coding RNAs whose role in autoimmune disease pathology has not been well characterised.

METHODS

Whole blood RNA-sequencing (RNA-seq) was performed on SjD cases (n=23 Ro/SSA negative (SjDRo-); n=27 Ro/SSA positive (SjDRo+) and healthy controls (HCs; n=27). Bioinformatics and pathway analyses of differentially expressed (DE) transcripts (log2 fold change ≥2 or ≤0.5; padj<0.05) were used to predict lncRNA function. LINC01871 was characterised by RNA-seq analyses of HSB-2 cells with CRISPR-targeted LINC01871 deletion (LINC01871-/ -) and in vitro stimulation assays.

RESULTS

Whole blood RNA-seq revealed autoantibody-specific transcription profiles and disproportionate downregulation of DE transcripts in SjD cases relative to HCs. Sixteen DE lncRNAs exhibited correlated expression with the interferon (IFN)-regulated gene, RSAD2, in SjDRo+ (r≥0.65 or ≤-0.6); four antisense lncRNAs exhibited IFN-regulated expression in immune cell lines. LINC01871 was upregulated in all SjD cases. RNA-seq and pathway analyses of LINC01871-/ - cells implicated roles in cytotoxic function, differentiation and IFNγ induction. LINC01871 was induced by IFNγ in a myeloid cell line and regulated by calcineurin/NFAT pathway and T cell receptor (TCR) signalling in primary human T cells.

CONCLUSION

LINC01871 influences expression of many immune cell genes and growth factors, is IFNγ inducible, and regulated by calcineurin signalling and TCR ligand engagement. Altered LINC01871 expression may influence the dysregulated T cell inflammatory pathways implicated in SjD.

Genome-wide association study identifies Sjögren's risk loci with functional implications in immune and glandular cells

Sjögren's disease is a complex autoimmune disease with twelve established susceptibility loci. This genome-wide association study (GWAS) identifies ten novel genome-wide significant (GWS) regions in Sjögren's cases of European ancestry: CD247, NAB1, PTTG1-MIR146A, PRDM1-ATG5, TNFAIP3, XKR6, MAPT-CRHR1, RPTOR-CHMP6-BAIAP6, TYK2, SYNGR1. Polygenic risk scores yield predictability (AUROC = 0.71) and relative risk of 12.08. Interrogation of bioinformatics databases refine the associations, define local regulatory networks of GWS SNPs from the 95% credible set, and expand the implicated gene list to >40. Many GWS SNPs are eQTLs for genes within topologically associated domains in immune cells and/or eQTLs in the main target tissue, salivary glands.

Defective Efferocytosis in a Murine Model of Sjögren's Syndrome Is Mediated by Dysfunctional Mer Tyrosine Kinase Receptor

Sjögren's syndrome (SjS) is a chronic autoimmune disease primarily involving the exocrine glands in which the involvement of the innate immune system is largely uncharacterized. Mer signaling has been found to be protective in several autoimmune diseases but remains unstudied in SjS. Here, we investigated the role of Mer signaling in SjS. Mer knockout (MerKO) mice were examined for SjS disease criteria. SjS-susceptible (SjS^S) C57BL/6.NOD-Aec1Aec2 mice were assessed for defective Mer signaling outcomes, soluble Mer (sMer) levels, A disintegrin and metalloprotease 17 (ADAM17) activity, and Rac1 activation. In addition, SjS patient plasma samples were evaluated for sMer levels via ELISA, and sMer levels were correlated to disease manifestations. MerKO mice developed submandibular gland (SMG) lymphocytic infiltrates, SMG apoptotic cells, anti-nuclear autoantibodies (ANA), and reduced saliva flow. Mer signaling outcomes were observed to be diminished in SjS^S mice, as evidenced by reduced Rac1 activation in SjS^S mice macrophages in response to apoptotic cells and impaired efferocytosis. Increased sMer was also detected in SjS^S mouse sera, coinciding with higher ADAM17 activity, the enzyme responsible for cleavage and inactivation of Mer. sMer levels were elevated in patient plasma and positively correlated with focus scores, ocular staining scores, rheumatoid factors, and anti-Ro60 levels. Our data indicate that Mer plays a protective role in SjS, similar to other autoimmune diseases. Furthermore, we suggest a series of events where enhanced ADAM17 activity increases Mer inactivation and depresses Mer signaling, thus removing protection against the loss of self-tolerance and the onset of autoimmune disease in SjS^S mice.

Transcriptomic and Network Analysis of Minor Salivary Glands of Patients With Primary Sjögren's Syndrome

Primary Sjögren’s syndrome (pSS) is a systemic autoimmune disease characterized primarily by immune-mediated destruction of exocrine tissues, such as those of the salivary and lacrimal glands, resulting in the loss of saliva and tear production, respectively. This disease predominantly affects middle-aged women, often in an insidious manner with the accumulation of subtle changes in glandular function occurring over many years. Patients commonly suffer from pSS symptoms for years before receiving a diagnosis. Currently, there is no effective cure for pSS and treatment options and targeted therapy approaches are limited due to a lack of our overall understanding of the disease etiology and its underlying pathology. To better elucidate the underlying molecular nature of this disease, we have performed RNA-sequencing to generate a comprehensive global gene expression profile of minor salivary glands from an ethnically diverse cohort of patients with pSS. Gene expression analysis has identified a number of pathways and networks that are relevant in pSS pathogenesis. Moreover, our detailed integrative analysis has revealed a primary Sjögren’s syndrome molecular signature that may represent important players acting as potential drivers of this disease. Finally, we have established that the global transcriptomic changes in pSS are likely to be attributed not only to various immune cell types within the salivary gland but also epithelial cells which are likely playing a contributing role. Overall, our comprehensive studies provide a database-enriched framework and resource for the identification and examination of key pathways, mediators, and new biomarkers important in the pathogenesis of this disease with the long-term goals of facilitating earlier diagnosis of pSS and to mitigate or abrogate the progression of this debilitating disease.

Unique Sjögren's syndrome patient subsets defined by molecular features

OBJECTIVE

To address heterogeneity complicating primary SS (pSS) clinical trials, research and care by characterizing and clustering patients by their molecular phenotypes.

METHODS

pSS patients met American-European Consensus Group classification criteria and had at least one systemic manifestation and stimulated salivary flow of ⩾0.1 ml/min. Correlated transcriptional modules were derived from gene expression microarray data from blood (n = 47 with appropriate samples). Patients were clustered based on this molecular information using an unbiased random forest modelling approach. In addition, multiplex, bead-based assays and ELISAs were used to assess 30 serum cytokines, chemokines and soluble receptors. Eleven autoantibodies, including anti-Ro/SSA and anti-La/SSB, were measured by Bio-Rad Bioplex 2200.

RESULTS

Transcriptional modules distinguished three clusters of pSS patients. Cluster 1 showed no significant elevation of IFN or inflammation modules. Cluster 2 showed strong IFN and inflammation modular network signatures, as well as high plasma protein levels of IP-10/CXCL10, MIG/CXCL9, BLyS (BAFF) and LIGHT. Cluster 3 samples exhibited moderately elevated IFN modules, but with suppressed inflammatory modules, increased IP-10/CXCL10 and B cell-attracting chemokine 1/CXCL13 and trends toward increased MIG/CXCL9, IL-1α, and IL-21. Anti-Ro/SSA and anti-La/SSB were present in all three clusters.

CONCLUSION

Molecular profiles encompassing IFN, inflammation and other signatures can be used to separate patients with pSS into distinct clusters. In the future, such profiles may inform patient selection for clinical trials and guide treatment decisions.

Sjögren's Syndrome Minor Salivary Gland CD4+ Memory T Cells Associate with Glandular Disease Features and have a Germinal Center T Follicular Helper Transcriptional Profile

To assess the types of salivary gland (SG) T cells contributing to Sjögren's syndrome (SS), we evaluated SG T cell subtypes for association with disease features and compared the SG CD4⁺ memory T cell transcriptomes of subjects with either primary SS (pSS) or non-SS sicca (nSS). SG biopsies were evaluated for proportions and absolute numbers of CD4⁺ and CD8⁺ T cells. SG memory CD4⁺ T cells were evaluated for gene expression by microarray. Differentially-expressed genes were identified, and gene set enrichment and pathways analyses were performed. CD4⁺CD45RA^- T cells were increased in pSS compared to nSS subjects (33.2% vs. 22.2%, p < 0.0001), while CD8⁺CD45RA^- T cells were decreased (38.5% vs. 46.0%, p = 0.0014). SG fibrosis positively correlated with numbers of memory T cells. Proportions of SG CD4⁺CD45RA^- T cells correlated with focus score (r = 0.43, p < 0.0001), corneal damage (r = 0.43, p < 0.0001), and serum Ro antibodies (r = 0.40, p < 0.0001). Differentially-expressed genes in CD4⁺CD45RA^- cells indicated a T follicular helper (Tfh) profile, increased homing and increased cellular interactions. Predicted upstream drivers of the Tfh signature included TCR, TNF, TGF-β1, IL-4, and IL-21. In conclusion, the proportions and numbers of SG memory CD4⁺ T cells associate with key SS features, consistent with a central role in disease pathogenesis.

American Indians Have a Higher Risk of Sjögren's Syndrome and More Disease Activity Than European Americans and African Americans

OBJECTIVE

To describe the clinical and serologic manifestations of Sjögren's syndrome (SS) in ethnic groups of the US.

METHODS

This was a cross-sectional study of 648 patients with primary SS: 20 African American (AA), 164 American Indian (AI), 426 European American (EA), and 38 patients of other races evaluated in a multidisciplinary Sjögren's International Collaborative Clinical Alliance research clinic.

RESULTS

AA subjects comprised 3.1% of the SS cohort, much lower than the percentage of AA in the state of Oklahoma (P = 3.01 × E - 05), the US (P = 2.24E - 13), or a systemic lupus erythematosus (SLE) cohort at the same institution (P = 4.26 × 10E - 27). In contrast, the percentage of AI in the SS cohort (25.3%) was much higher than expected (P = 3.17E - 09 versus SLE cohort, P = 6.36 - 26 versus Oklahoma, and P = 8.14E - 96 versus US population). The SS classification criteria were similar between AA and EA, but subjects of AI ancestry had lower rates of abnormal tear and salivary flow, as well as anti-Ro/SSA and anti-La/SSB antibodies. Paradoxically, AI had higher levels of disease activity (mean ± SD European League Against Rheumatism Sjögren's Syndrome Disease Activity Index score 3.77 ± 4.78) in comparison to EA (2.90 ± 4.12; P = 0.011) and more extraglandular manifestations affecting mainly the articular and glandular domains. Meanwhile, AA patients were characterized by higher rates of hypergammaglobulinemia (odds ratio [OR] 1.39 [95% confidence interval (95% CI) 1.39-8.65]; P = 0.01), elevated erythrocyte sedimentation rate (ESR) (OR 3.95 [95% CI 1.46-9.95]; P = 0.009), and parotid enlargement (OR 4.40 [95% CI 1.49-13.07]; P = 0.02).

CONCLUSION

AI are affected at high rates with SS but present with few classical features, potentially preventing timely diagnosis. In contrast to SLE, SS is infrequent and not more severe among AA, but the triad of hypergammaglobulinemia, increased ESR, and parotid enlargement warrants extra vigilance for lymphomagenesis.

Novel shared antibody specificities in anti-Ro/ La antibody negative Sjögren’s Syndrome

INTRODUCTION

Sjögren’s syndrome (SS) is a rheumatic autoimmune disease characterized by focal lymphocytic infiltrates in the lacrimal and salivary glands, severe dry mouth and eyes, pain and debilitation. Diagnosis requires autoantibodies to ubiquitous Ro antigens or a lip biopsy positive for focal lymphocytic infiltrates. Here we used human proteome arrays to identify novel antibodies in plasma from Ro positive and Ro/La antibody negative SS patients compared with healthy controls.

METHODS

Anti-Ro positive (n=15) and anti-Ro negative (n=15) cases meeting 2016 ACR/EULAR classification criteria for SS were age, race, and sex matched with each other and healthy controls (n=15). Plasma IgG binding to human proteome arrays containing &gt;19,500 recombinant human proteins representing &gt;80% of the human proteome (HuProt v3.2 arrays, CDI Laboratories) was assessed. Data were normalized by the Robust Linear Model using the PAA Bioconductor Package in R and log intensity values for each protein generated. Thresholds of mean + 4SD were established using the controls. Antigens bound by IgG more frequently in cases compared to controls (p&lt;0.05, one-tailed Fisher’s exact test) were considered significant.

RESULTS

IgG from Ro positive SS cases significantly bound 18 proteins, including the canonical SS antigens Ro60 and Ro52. IgG from Ro negative SS cases significantly bound 4 proteins compared to controls, 3 of which were shared with the Ro positive group. Binding to any one of 4 novel proteins identified 73% of the Ro negative SS cases.

CONCLUSION

A total of 17 novel antigen specificities were identified in SS, with 4 antigens being bound by plasma IgG from Ro/La negative SS cases. These antigens may be useful for diagnosing SS without a lip biopsy.

Sjögren’s syndrome minor salivary gland CD4+ T cells associate with glandular disease features and have a germinal center T follicular helper transcriptional profile

Objective

To assess the types of salivary gland (SG) T cells contributing to Sjögren’s syndrome (SS), we evaluated SG T cell subtypes for association with disease features and compared the SG CD4+ memory T cell transcriptomes of primary SS (pSS) and sicca subjects not meeting SS criteria (nSS).

Methods

SG biopsies of pSS and nSS subjects were evaluated for proportions (n=51 pSS, n=69 nSS) and absolute numbers (n=35 pSS, n=57 nSS) of CD4+ and CD8+ T cells. SG memory CD4+ T cells from focus score positive pSS (n=17) and focus score negative nSS subjects (n=15) were evaluated for gene expression by microarray. Differentially-expressed genes were identified and gene set enrichment and pathways analyses were performed.

Results

CD4+CD45RA− T cell frequencies were increased in pSS compared to nSS subjects (33.2 vs. 22.2%, p&lt;0.0001), while CD8+CD45RA− T cells were decreased (38.5 vs. 46.0%, p=0.0014). SG fibrosis positively correlated with numbers of memory T cells, regardless of phenotype. Proportions of SG CD4+CD45RA− T cells correlated with focus score (r=0.43, p&lt;0.0001), corneal damage (r=0.43, p&lt;0.0001), and serum Ro antibodies (r=0.40, p&lt;0.0001). Differentially expressed genes in CD4+CD45RA− cells of pSS versus nSS subjects indicated a T follicular helper (Tfh) profile, increased homing and increased cellular interactions. Predicted upstream drivers of the Tfh signature included TCR, TNF, TGF-β1, IL-4 and IL-21.

Conclusions

The proportions and numbers of SG memory CD4+ T cells associate with key SS features, consistent with a central role in disease pathogenesis. SG memory CD4+ T cells express a germinal center Tfh profile.

Sjögren Syndrome without Focal Lymphocytic Infiltration of the Salivary Glands

OBJECTIVE

Primary Sjögren syndrome (SS) is characterized by a focal lymphocytic infiltrate in exocrine glands. We describe patients who lacked this key feature.

METHODS

We evaluated patients with sicca in a comprehensive clinic at which medical, dental, and ophthalmological examinations were performed. All subjects underwent a minor salivary gland biopsy with focus score calculation. Extraglandular manifestations were also determined. We categorized subjects as high, intermediate, or low in terms of expression of interferon (IFN)-regulated genes.

RESULTS

About 20% (51 of 229, 22%) of those classified as having primary SS had a focus score of zero. Compared to those with anti-Ro positivity and a focus score > 1.0, the patients with focus score of zero (who by classification criteria must be anti-Ro-positive) were statistically less likely to have anti-La (or SSB) and elevated immunoglobulin, as well as less severe corneal staining. The focus score zero patients were less likely to have elevated expression of IFN-regulated genes in peripheral blood mononuclear cells than anti-Ro-positive SS patients with a focal salivary infiltrate.

CONCLUSION

There are only a few clinical differences between patients with primary SS with focus score zero and those with both anti-Ro and a focus score > 1.0. The small subset of focus score zero patients tested did not have elevated expression of IFN-regulated genes, but did have systemic disease. Thus, extraglandular manifestations are perhaps more related to the presence of anti-Ro than increased IFN. This may have relevance to pathogenesis of SS.

Reproducibility of Ocular Surface Staining in the Assessment of Sjögren Syndrome-Related Keratoconjunctivitis Sicca: Implications on Disease Classification

Objective

The objective of this study was to assess the performance and reproducibility of the two currently used ocular surface staining scores in the assessment of keratoconjunctivitis sicca in Sjögren syndrome (SS) research classification.

Methods

In a multidisciplinary clinic for the evaluation of sicca, we performed all tests for the American European Consensus Group (AECG) and the American College of Rheumatology (ACR)/European League Against Rheumatism (EULAR) classification criteria, including the van Bijsterveld score (vBS) and the Ocular Staining Score (OSS), in 994 participants with SS or with non-SS sicca. We analyzed the concordance between the scores, the diagnostic accuracy and correlation with clinical variables, and interrater and intrasubject reproducibility.

Results

A total of 308 (31.1%) participants had a discordant vBS and OSS that was due to extra corneal staining points in the OSS. The presence of one or more of the additional points was highly predictive of SS classification (odds ratio = 3.66; P = 1.65 × 10e-20) and was associated with abnormal results of all measures of autoimmunity and glandular dysfunction. Receiver operating characteristic curves showed optimal cutoff values of four for the vBS (sensitivity = 0.62; specificity = 0.71; Youden's J = 0.33) and five for the OSS (sensitivity = 0.56; specificity = 0.75; Youden's J = 0.31). Notably, there was very poor consistency in interobserver mean scores and distributions (P < 0.0001) and in intrasubject scores after a median of 5.5 years (35% changed status of the ocular criterion).

Conclusion

Ocular surface staining scores are useful for SS research classification; however, they are subject to significant interrater and intrasubject variability, which could result in changes in classification in 5%-10% of all subjects. These results highlight the need for objective and reproducible markers of disease that have thus far remained elusive for SS.

Genetics in Sjögren's syndrome: where we are and where we go

Sjögren's syndrome is a complex autoimmune disease that involves dysregulation of immune responses that preferentially target exocrine glands. Systemic manifestations vary and may involve nearly every organ system. Genetic studies to date are in their infancy relative to other autoimmune diseases such as systemic lupus erythematosus and rheumatoid arthritis, each with more than 100 genetic associations now established. However, recent work in SS has successfully established associations that shed light on pathophysiology and implicate aberrant innate and adaptive immune responses. In this review, we provide an overview of genetic approaches used to identify risk variants in SS, discuss major findings and their relevance to SS, and describe the future directions that are likely to lead to understanding fundamental causes of this disease and new opportunities for improving clinical care.

Genetic fine mapping of systemic lupus erythematosus MHC associations in Europeans and African Americans

Genetic variation within the major histocompatibility complex (MHC) contributes substantial risk for systemic lupus erythematosus, but high gene density, extreme polymorphism and extensive linkage disequilibrium (LD) have made fine mapping challenging. To address the problem, we compared two association techniques in two ancestrally diverse populations, African Americans (AAs) and Europeans (EURs). We observed a greater number of Human Leucocyte Antigen (HLA) alleles in AA consistent with the elevated level of recombination in this population. In EUR we observed 50 different A-C-B-DRB1-DQA-DQB multilocus haplotype sequences per hundred individuals; in the AA sample, these multilocus haplotypes were twice as common compared to Europeans. We also observed a strong narrow class II signal in AA as opposed to the long-range LD observed in EUR that includes class I alleles. We performed a Bayesian model choice of the classical HLA alleles and a frequentist analysis that combined both single nucleotide polymorphisms (SNPs) and classical HLA alleles. Both analyses converged on a similar subset of risk HLA alleles: in EUR HLA- B*08:01 + B*18:01 + (DRB1*15:01 frequentist only) + DQA*01:02 + DQB*02:01 + DRB3*02 and in AA HLA-C*17:01 + B*08:01 + DRB1*15:03 + (DQA*01:02 frequentist only) + DQA*02:01 + DQA*05:01+ DQA*05:05 + DQB*03:19 + DQB*02:02. We observed two additional independent SNP associations in both populations: EUR rs146903072 and rs501480; AA rs389883 and rs114118665. The DR2 serotype was best explained by DRB1*15:03 + DQA*01:02 in AA and by DRB1*15:01 + DQA*01:02 in EUR. The DR3 serotype was best explained by DQA*05:01 in AA and by DQB*02:01 in EUR. Despite some differences in underlying HLA allele risk models in EUR and AA, SNP signals across the extended MHC showed remarkable similarity and significant concordance in direction of effect for risk-associated variants.

Novel genetic associations with interferon in systemic lupus erythematosus identified by replication and fine-mapping of trait-stratified genome-wide screen

BACKGROUND/PURPOSE

High serum interferon alpha (IFN-α) is an important heritable phenotype in systemic lupus erythematosus (SLE) which is involved in primary disease pathogenesis. High vs. low levels of IFN-α are associated with disease severity and account for some of the biological heterogeneity between SLE patients. The aim of the study was to replicate and fine-map previously detected genetic associations with serum IFN-α in SLE.

METHODS

We previously undertook a case-case genome-wide association study of SLE patients stratified by ancestry and extremes of phenotype in serum IFN-α. Single nucleotide polymorphisms (SNPs) in seven loci identified in this screen were selected for follow up in a large independent cohort of 1370 SLE patients (703 European-ancestry, 432 African ancestry, and 235 Amerindian ancestry). Each ancestral background was analyzed separately, and ancestry-informative markers were used to control for ancestry and admixture.

RESULTS

We find a rare haplotype spanning the promoter region of EFNA5 that is strongly associated with serum IFN-α in both African-American and European-American SLE patients (OR = 3.0, p = 3.7 × 10^-6). We also find SNPs in the PPM1H, PTPRM, and NRGN regions associated with IFN-α levels in European-American, Amerindian, and African-American SLE patients respectively. Many of these associations are within regulatory regions of the gene, suggesting an impact on transcription.

CONCLUSION

This study demonstrates the power of molecular sub-phenotypes to reveal genetic factors involved in complex autoimmune disease. The distinct associations observed in different ancestral backgrounds emphasize the heterogeneity of molecular pathogenesis in SLE.

A plausibly causal functional lupus-associated risk variant in the STAT1-STAT4 locus

Systemic lupus erythematosus (SLE or lupus) (OMIM: 152700) is a chronic autoimmune disease with debilitating inflammation that affects multiple organ systems. The STAT1-STAT4 locus is one of the first and most highly replicated genetic loci associated with lupus risk. We performed a fine-mapping study to identify plausible causal variants within the STAT1-STAT4 locus associated with increased lupus disease risk. Using complementary frequentist and Bayesian approaches in trans-ancestral Discovery and Replication cohorts, we found one variant whose association with lupus risk is supported across ancestries in both the Discovery and Replication cohorts: rs11889341. In B cell lines from patients with lupus and healthy controls, the lupus risk allele of rs11889341 was associated with increased STAT1 expression. We demonstrated that the transcription factor HMGA1, a member of the HMG transcription factor family with an AT-hook DNA-binding domain, has enriched binding to the risk allele compared with the non-risk allele of rs11889341. We identified a genotype-dependent repressive element in the DNA within the intron of STAT4 surrounding rs11889341. Consistent with expression quantitative trait locus (eQTL) analysis, the lupus risk allele of rs11889341 decreased the activity of this putative repressor. Altogether, we present a plausible molecular mechanism for increased lupus risk at the STAT1-STAT4 locus in which the risk allele of rs11889341, the most probable causal variant, leads to elevated STAT1 expression in B cells due to decreased repressor activity mediated by increased binding of HMGA1.

Evidence of Alternative Modes of B Cell Activation Involving Acquired Fab Regions of N-Glycosylation in Antibody-Secreting Cells Infiltrating the Labial Salivary Glands of Patients With Sjögren's Syndrome

OBJECTIVE

To better understand the role of B cells, the potential mechanisms responsible for their aberrant activation, and the production of autoantibodies in the pathogenesis of Sjögren's syndrome (SS), this study explored patterns of selection pressure and sites of N-glycosylation acquired by somatic mutation (acN-glyc) in the IgG variable (V) regions of antibody-secreting cells (ASCs) isolated from the minor salivary glands of patients with SS and non-SS control patients with sicca symptoms.

METHODS

A novel method to produce and characterize recombinant monoclonal antibodies (mAb) from single cell-sorted ASC infiltrates was applied to concurrently probe expressed genes (all heavy- and light-chain isotypes as well as any other gene of interest not related to immunoglobulin) in the labial salivary glands of patients with SS and non-SS controls. V regions were amplified by reverse transcription-polymerase chain reaction, sequenced, and analyzed for the incidence of N-glycosylation and selection pressure. For specificity testing, the amplified regions were expressed as either the native mAb or mutant mAb lacking the acN-glyc motif. Protein modeling was used to demonstrate how even an acN-glyc site outside of the complementarity-determining region could participate in, or inhibit, antigen binding.

RESULTS

V-region sequence analyses revealed clonal expansions and evidence of secondary light-chain editing and allelic inclusion, of which neither of the latter two have previously been reported in patients with SS. Increased frequencies of acN-glyc were found in the sequences from patients with SS, and these acN-glyc regions were associated with an increased number of replacement mutations and lowered selection pressure. A clonal set of polyreactive mAb with differential framework region 1 acN-glyc motifs was also identified, and removal of the acN-glyc could nearly abolish binding to autoantigens.

CONCLUSION

These findings support the notion of an alternative mechanism for the selection and proliferation of some autoreactive B cells, involving V-region N-glycosylation, in patients with SS.

Minor salivary gland fibrosis in Sjögren's syndrome is elevated, associated with focus score and not solely a consequence of aging

OBJECTIVES

Evaluate the presence of minor salivary gland (SG) fibrosis in primary Sjögren's syndrome (pSS) as a function of disease pathology or a consequence of ageing.

METHODS

Subjects with sicca symptoms attending a Sjögren's research clinic were classified by American European Consensus Group (AECG) criteria as either pSS or non-SS (nSS). Discovery (n=34 pSS, n=28 nSS) and replication (n=35 pSS, n=31 nSS) datasets were evaluated. Minor SG cross-sections from haematoxylin and eosin stained slides were imaged, digitally reconstructed and analysed for percent area fibrosis. Relationships between SG fibrosis, age, and clinical measures were evaluated using Spearman correlations. Association with SS was assessed by: ROC curve, Variable Selection Using Random Forests (VSURF) and uni- and bi-variate regression analyses.

RESULTS

SS subjects had significantly more fibrotic tissue in their minor labial salivary glands (median 24.39%, range 5.12-51.67%) than nSS participants (median 16.7%, range 5.97-38.65%, p<0.0001); age did not differ between groups (average ± SD pSS 50.2 ±13.9 years, nSS 53.8±12.4 years). In both the discovery and replication data sets, multiple regression models showed that the area of minor salivary gland fibrosis predicted pSS significantly better than age alone. Age-corrected linear regression revealed that the area of minor salivary gland fibrosis positively associated with vanBijsterveld score (p=0.042) and biopsy focus score (p=0.002). ROC curve and VSURF analyses ranked fibrosis as a significantly more important variable for subject discrimination than age.

CONCLUSIONS

SG fibrosis is an element of pSS pathology that is related to focus score and is not solely attributable to age.

How immunological profile drives clinical phenotype of primary Sjögren's syndrome at diagnosis: analysis of 10,500 patients (Sjögren Big Data Project)

OBJECTIVES

To evaluate the influence of the main immunological markers on the disease phenotype at diagnosis in a large international cohort of patients with primary Sjögren's syndrome (SjS).

METHODS

The Big Data Sjögren Project Consortium is an international, multicentre registry created in 2014. As a first step, baseline clinical information from leading centres on clinical research in SjS of the 5 continents was collected. The centres shared a harmonised data architecture and conducted cooperative online efforts in order to refine collected data under the coordination of a big data statistical team. Inclusion criteria were the fulfillment of the 2002 classification criteria. Immunological tests were carried out using standard commercial assays.

RESULTS

By January 2018, the participant centres had included 10,500 valid patients from 22 countries. The cohort included 9,806 (93%) women and 694 (7%) men, with a mean age at diagnosis of primary SjS of 53 years, mainly White (78%) and included from European countries (71%). The frequency of positive immunological markers at diagnosis was 79.3% for ANA, 73.2% for anti-Ro, 48.6% for RF, 45.1% for anti- La, 13.4% for low C3 levels, 14.5% for low C4 levels and 7.3% for cryoglobulins. Positive autoantibodies (ANA, Ro, La) correlated with a positive result in salivary gland biopsy, while hypocomplementaemia and especially cryoglo-bulinaemia correlated with systemic activity (mean ESSDAI score of 17.7 for cryoglobulins, 11.3 for low C3 and 9.2 for low C4, in comparison with 3.8 for negative markers). The immunological markers with a great number of statistically-significant associations (p<0.001) in the organ-by-organ ESS- DAI evaluation were cryoglobulins (9 domains), low C3 (8 domains), anti-La (7 domains) and low C4 (6 domains).

CONCLUSIONS

We confirm the strong influence of immunological markers on the phenotype of primary SjS at diagnosis in the largest multi-ethnic international cohort ever analysed, with a greater influence for cryoglobulinaemic-related markers in comparison with Ro/La autoantibodies and ANA. Immunological patterns play a central role in the phenotypic expression of the disease already at the time of diagnosis, and may guide physicians to design a specific personalised management during the follow-up of patients with primary SjS.

Unique clinical characteristics, autoantibodies and medication use in Native American patients with systemic lupus erythematosus

Objective

Systemic lupus erythematosus (SLE) is a systemic autoimmune disease with varied morbidity and mortality. We assessed clinical presentations, autoantibody specificities and therapeutic interventions in Native American (NA) patients with SLE.

Methods

Patients with SLE meeting 1997 American College of Rheumatology classification criteria (n=3148) were enrolled between 1992 and 2010 in the multiethnic, cross-sectional Lupus Family Registry and Repository. Clinical, demographic and therapeutic information were extracted from medical records using a standardised form and formalised training. Autoantibodies were assessed by indirect immunofluorescence (antinuclear antibodies (ANA) and antidouble-stranded DNA), precipitin (ENA) and ELISA (IgG and IgM anticardiolipins).

Results

NA patients met SLE classification at a younger age (29.89±12.3 years) than European Americans (EA; 32.02±12.87, P=0.0157) and a similar age to African-Americans (AAs) and Hispanics (HIS). More NA patients had concurrent rheumatic diseases or symptoms, such as Raynaud’s phenomenon, interstitial lung disease, Sjӧgren’s syndrome and systemic sclerosis. Compared with EAs, NAs were more likely to have high-titre ANA (≥1:3240; P<0.0001) and had more SLE-associated autoantibodies. Autoantibodies with unknown specificities were more common in NAs (41%) compared with other racial/ethnic groups in this collection (AA: 24%, P=0.0006; EA: 17%, P<0.0001; HIS: 23%, P=0.0050). Fewer NA patients used hydroxychloroquine (68%) compared with others (AA: 74%, P=0.0308; EA: 79%, P=0.0001, HIS: 77%, P=0.0173); this was influenced by lower hydroxychloroquine use in NA patients from Latin America (32%). NA patients had higher rates of methotrexate use (28%) compared with AA (18%, P=0.0006) and HIS patients (14%, P=0.0003), higher azathioprine use (38%) compared with EA patients (30%, P=0.0105) and higher mycophenolate mofetil use (26%) compared with EA (17%, P=0.0012) and HIS patients (11%, P<0.0001).

Conclusions

NA patients are diagnosed with SLE earlier in life and present worse concurrent rheumatic disease symptoms than EA patients. NA patients also are more likely to have expanded autoantibody profiles and precipitins of unknown specificities.

Accumulation of Antigen-Driven Lymphoproliferations in Complement Receptor 2/CD21-/low B Cells From Patients With Sjögren's Syndrome

OBJECTIVE

Patients with Sjögren's syndrome (SS) are prone to develop malignant lymphomas, and a correlation has been established between the lymphoproliferations occurring in these disorders and the presence in patients' blood of an unusual B cell population that down-regulates complement receptor 2/CD21. This study was undertaken to identify the B cell compartment from which these lymphoproliferations emerge and determine the mechanisms that promote clonal B cell expansion in patients with SS.

METHODS

The reactivity of antibodies expressed by CD19+CD10-CD27-IgM+CD21^-/low cells isolated from the blood of patients with SS was tested using a polymerase chain reaction-based approach that allows us to clone and express, in vitro, recombinant antibodies produced by single B cells.

RESULTS

Clonal expansions were identified in CD21^-/low B cells isolated from the peripheral blood of 3 patients with SS. These lymphoproliferations expressed B cell receptors (BCRs) that displayed somatic hypermutation lineage trees characteristic of a strong selection by antigens; one of these antigens was identified as a ribosomal self antigen. When the mutated BCR sequences expressed by the expanded CD21^-/low B cell clones from patients with SS were reverted in vitro to their germline counterparts, one clone remained autoreactive.

CONCLUSION

Clonal lymphoproliferations in patients with SS preferentially accumulate in the autoreactive CD21^-/low B cell compartment often expanded in these subjects, and recognition of self antigens may drive the clonal B cell expansion while further refining BCR self-reactivity.

Characterization and classification of lupus patients based on plasma thermograms

Objective

Plasma thermograms (thermal stability profiles of blood plasma) are being utilized as a new diagnostic approach for clinical assessment. In this study, we investigated the ability of plasma thermograms to classify systemic lupus erythematosus (SLE) patients versus non SLE controls using a sample of 300 SLE and 300 control subjects from the Lupus Family Registry and Repository. Additionally, we evaluated the heterogeneity of thermograms along age, sex, ethnicity, concurrent health conditions and SLE diagnostic criteria.

Methods

Thermograms were visualized graphically for important differences between covariates and summarized using various measures. A modified linear discriminant analysis was used to segregate SLE versus control subjects on the basis of the thermograms. Classification accuracy was measured based on multiple training/test splits of the data and compared to classification based on SLE serological markers.

Results

Median sensitivity, specificity, and overall accuracy based on classification using plasma thermograms was 86%, 83%, and 84% compared to 78%, 95%, and 86% based on a combination of five antibody tests. Combining thermogram and serology information together improved sensitivity from 78% to 86% and overall accuracy from 86% to 89% relative to serology alone. Predictive accuracy of thermograms for distinguishing SLE and osteoarthritis / rheumatoid arthritis patients was comparable. Both gender and anemia significantly interacted with disease status for plasma thermograms (p<0.001), with greater separation between SLE and control thermograms for females relative to males and for patients with anemia relative to patients without anemia.

Conclusion

Plasma thermograms constitute an additional biomarker which may help improve diagnosis of SLE patients, particularly when coupled with standard diagnostic testing. Differences in thermograms according to patient sex, ethnicity, clinical and environmental factors are important considerations for application of thermograms in a clinical setting.

Clinical and Serologic Features in Patients With Incomplete Lupus Classification Versus Systemic Lupus Erythematosus Patients and Controls

OBJECTIVE

Incomplete lupus erythematosus (ILE) involves clinical and/or serologic manifestations consistent with but insufficient for systemic lupus erythematosus (SLE) classification. Because the nature of ILE is poorly understood and no treatment recommendations exist, we examined the clinical manifestations, medication history, and immunologic features in a diverse collection of ILE and SLE patients.

METHODS

Medical records of subjects enrolled in the Lupus Family Registry and Repository were reviewed for medication history and American College of Rheumatology (ACR) classification criteria to identify ILE patients (3 ACR criteria; n = 440) and SLE patients (≥4 ACR criteria; n = 3,397). Participants completed the Connective Tissue Disease Screening Questionnaire. Anticardiolipin and plasma B lymphocyte stimulator (BLyS) were measured by enzyme-linked immunosorbent assay, antinuclear antibodies (ANAs) by indirect immunofluorescence, and 13 autoantibodies by bead-based assays.

RESULTS

On average, ILE patients were older than SLE patients (46.2 years versus 42.0 years; P < 0.0001), and fewer ILE patients were African American (23.9% versus 32.2%; P < 0.001). ILE patients exhibited fewer autoantibody specificities than SLE patients (1.3 versus 2.6; P < 0.0001) and were less likely to have ANA titers ≥1:1,080 (10.5% versus 19.5%; P < 0.0001). BLyS levels were intermediate in ILE patients (controls < ILE; P = 0.016; ILE < SLE; P = 0.008). Pericarditis, renal, or neurologic manifestations occurred in 12.5% of ILE patients and were associated with non-European American race/ethnicity (P = 0.012). Hydroxychloroquine use increased over time, but was less frequent in ILE than SLE patients (65.2% versus 83.1%; P < 0.0001).

CONCLUSION

Although usually characterized by milder symptoms, ILE manifestations may require immunomodulatory treatments. Longitudinal studies are necessary to understand how ILE affects organ damage and future SLE risk, and to delineate molecular pathways unique to ILE.

A new MHC-linked susceptibility locus for primary Sjögren's syndrome: MICA

The association of primary Sjögren's syndrome (pSS) with Major Histocompatibility Complex (MHC) alleles is quintessential of MHC-disease associations. Indeed, although disease associations with classical HLA class I and II alleles/haplotypes are amply documented, further dissection is often prevented by the strong linkage disequilibrium across the entire MHC complex. Here we study the association of pSS, not with HLA genes, but with the non-conventional MHC encoded class I gene, MICA (MHC class I chain-related gene A). MICA is selectively expressed within epithelia, and is the major ligand for the activatory receptor, NKG2D, both attributes relevant to pSS' etiology. MICA-pSS association was studied in two independent (French and UK) cohorts representing a total of 959 cases and 1,043 controls. MICA*008 allele was shown to be significantly associated with pSS (pcor=2.61 × 10-35). A multivariate logistic regression showed that this association was independent of all major known MHC-linked risk loci/alleles, as well as other relevant candidate loci that are in linkage disequilibrium with MICA*008 i.e. HLA-B*08:01, rs3131619 (T), MICB*008, TNF308A, HLA-DRB1*03:01 and HLA-DRB1*15:01 (P = 1.84 × 10-04). Furthermore, independently of the MICA*008 allele, higher levels of soluble MICA proteins were detected in sera of pSS patients compared to healthy controls. This study hence defines MICA as a new, MHC-linked, yet HLA-independent, pSS risk locus and opens a new front in our understanding of the still enigmatic pathophysiology of this disease. The fact that the soluble MICA protein is further amplified in MICA*008 carrying individuals, might also be relevant in other auto-immune diseases and cancer.

Rare X Chromosome Abnormalities in Systemic Lupus Erythematosus and Sjögren's Syndrome

OBJECTIVE

Sjögren's syndrome (SS) and systemic lupus erythematosus (SLE) are related by clinical and serologic manifestations as well as genetic risks. Both diseases are more commonly found in women than in men, at a ratio of ~10 to 1. Common X chromosome aneuploidies, 47,XXY and 47,XXX, are enriched among men and women, respectively, in either disease, suggesting a dose effect on the X chromosome.

METHODS

We examined cohorts of SS and SLE patients by constructing intensity plots of X chromosome single-nucleotide polymorphism alleles, along with determining the karyotype of selected patients.

RESULTS

Among ~2,500 women with SLE, we found 3 patients with a triple mosaic, consisting of 45,X/46,XX/47,XXX. Among ~2,100 women with SS, 1 patient had 45,X/46,XX/47,XXX, with a triplication of the distal p arm of the X chromosome in the 47,XXX cells. Neither the triple mosaic nor the partial triplication was found among the controls. In another SS cohort, we found a mother/daughter pair with partial triplication of this same region of the X chromosome. The triple mosaic occurs in ~1 in 25,000-50,000 live female births, while partial triplications are even rarer.

CONCLUSION

Very rare X chromosome abnormalities are present among patients with either SS or SLE and may inform the location of a gene(s) that mediates an X dose effect, as well as critical cell types in which such an effect is operative.

Fatty infiltration of the minor salivary glands is a selective feature of aging but not Sjögren's syndrome

OBJECTIVE

Determine the presence and assess the extent of fatty infiltration of the minor salivary glands (SG) of primary SS patients (pSS) as compared to those with non-SS sicca (nSS).

METHODS

Minor SG biopsy samples from 134 subjects with pSS (n = 72) or nSS (n = 62) were imaged. Total area and fatty replacement area for each glandular cross-section (n = 4-6 cross-sections per subject) were measured using Image J (National Institutes of Health, Bethesda, MD). The observer was blinded to subject classification status. The average area of fatty infiltration calculated per subject was evaluated by logistic regression and general linearized models (GLM) to assess relationships between fatty infiltration and clinical exam results, extent of fibrosis and age.

RESULTS

The average area of fatty infiltration for subjects with pSS (median% (range) 4.97 (0.05-30.2)) was not significantly different from that of those with nSS (3.75 (0.087-41.9). Infiltration severity varied widely, and subjects with fatty replacement greater than 6% were equivalently distributed between pSS and nSS participants (χ² p = .50). Age accounted for all apparent relationships between fatty infiltration and fibrosis or reduced saliva flow. The all-inclusive GLM for prediction of pSS versus non-SS classification including fibrosis, age, fatty replacement, and focus score was not significantly different from any desaturated model. In no iteration of the model did fatty replacement exert a significant effect on the capacity to predict pSS classification.

CONCLUSIONS

Fatty infiltration is an age-associated phenomenon and not a selective feature of Sjögren's syndrome. Sicca patients who do not fulfil pSS criteria have similar rates of fatty infiltration of the minor SG.

X Chromosome Dose and Sex Bias in Autoimmune Diseases: Increased Prevalence of 47,XXX in Systemic Lupus Erythematosus and Sjögren's Syndrome

OBJECTIVE

More than 80% of autoimmune disease predominantly affects females, but the mechanism for this female bias is poorly understood. We suspected that an X chromosome dose effect accounts for this, and we undertook this study to test our hypothesis that trisomy X (47,XXX; occurring in ∼1 in 1,000 live female births) would be increased in patients with female-predominant diseases (systemic lupus erythematosus [SLE], primary Sjögren's syndrome [SS], primary biliary cirrhosis, and rheumatoid arthritis [RA]) compared to patients with diseases without female predominance (sarcoidosis) and compared to controls.

METHODS

All subjects in this study were female. We identified subjects with 47,XXX using aggregate data from single-nucleotide polymorphism arrays, and, when possible, we confirmed the presence of 47,XXX using fluorescence in situ hybridization or quantitative polymerase chain reaction.

RESULTS

We found 47,XXX in 7 of 2,826 SLE patients and in 3 of 1,033 SS patients, but in only 2 of 7,074 controls (odds ratio in the SLE and primary SS groups 8.78 [95% confidence interval 1.67-86.79], P = 0.003 and odds ratio 10.29 [95% confidence interval 1.18-123.47], P = 0.02, respectively). One in 404 women with SLE and 1 in 344 women with SS had 47,XXX. There was an excess of 47,XXX among SLE and SS patients.

CONCLUSION

The estimated prevalence of SLE and SS in women with 47,XXX was ∼2.5 and ∼2.9 times higher, respectively, than that in women with 46,XX and ∼25 and ∼41 times higher, respectively, than that in men with 46,XY. No statistically significant increase of 47,XXX was observed in other female-biased diseases (primary biliary cirrhosis or RA), supporting the idea of multiple pathways to sex bias in autoimmunity.

Identification of a Sjögren's syndrome susceptibility locus at OAS1 that influences isoform switching, protein expression, and responsiveness to type I interferons

Sjögren's syndrome (SS) is a common, autoimmune exocrinopathy distinguished by keratoconjunctivitis sicca and xerostomia. Patients frequently develop serious complications including lymphoma, pulmonary dysfunction, neuropathy, vasculitis, and debilitating fatigue. Dysregulation of type I interferon (IFN) pathway is a prominent feature of SS and is correlated with increased autoantibody titers and disease severity. To identify genetic determinants of IFN pathway dysregulation in SS, we performed cis-expression quantitative trait locus (eQTL) analyses focusing on differentially expressed type I IFN-inducible transcripts identified through a transcriptome profiling study. Multiple cis-eQTLs were associated with transcript levels of 2'-5'-oligoadenylate synthetase 1 (OAS1) peaking at rs10774671 (PeQTL = 6.05 × 10-14). Association of rs10774671 with SS susceptibility was identified and confirmed through meta-analysis of two independent cohorts (Pmeta = 2.59 × 10-9; odds ratio = 0.75; 95% confidence interval = 0.66-0.86). The risk allele of rs10774671 shifts splicing of OAS1 from production of the p46 isoform to multiple alternative transcripts, including p42, p48, and p44. We found that the isoforms were differentially expressed within each genotype in controls and patients with and without autoantibodies. Furthermore, our results showed that the three alternatively spliced isoforms lacked translational response to type I IFN stimulation. The p48 and p44 isoforms also had impaired protein expression governed by the 3' end of the transcripts. The SS risk allele of rs10774671 has been shown by others to be associated with reduced OAS1 enzymatic activity and ability to clear viral infections, as well as reduced responsiveness to IFN treatment. Our results establish OAS1 as a risk locus for SS and support a potential role for defective viral clearance due to altered IFN response as a genetic pathophysiological basis of this complex autoimmune disease.

Anti-La positive, anti-Ro negative subset of primary Sjögren's syndrome: anti-La is a reality but is the disease?

OBJECTIVES

To characterise the serological and clinical findings in primary Sjögren's syndrome in which anti-La was found without anti-Ro. We hypothesised that a significant portion of these are falsely negative for anti-Ro60.

METHODS

Twenty-nine sera from primary Sjögren's syndrome patients were tested for antibodies directed against La and Ro. Anti-La was detected using bovine La treated with or without DNAase and RNAase to identify potential false positivity. Anti-Ro60 antibodies were detected using HEp-2000 substrate (in which cells are transfected with human Ro60) and HEp-2 substrate. Anti-Ro60 and Ro-52 were also tested by in vitro transcription/translation followed by immunoprecipitation assay.

RESULTS

All 29 sera bound La, even after treatment with DNAase and RNAase. Of the 29 sera, 25 were unequivocally negative on HEp-2000 (1:40 dilution). Four samples were anti-Ro60 positive with a speckled pattern, three of the four at 1:320 dilution. Thus, false negative anti-Ro60 exists in a small fraction (14%) of the Ro-negative/La-positive primary Sjögren's patients. However, all the samples were negative for Ro60 and Ro52 by in vitro immunoprecipitation assay. Clinically these patients tended not to have salivary gland pathology characteristic of Sjögren's syndrome.

CONCLUSIONS

We found only a small fraction of Ro negative/La positive sera to show positive HEp-2000 pattern. These subjects did not have characteristic findings on pathological examination of minor salivary glands, suggesting these subjects have a process distinct from Sjögren's syndrome.

Use of SLICC criteria in a large, diverse lupus registry enables SLE classification of a subset of ACR-designated subjects with incomplete lupus

Objective

SLE is traditionally classified using the American College of Rheumatology (ACR) criteria. The Systemic Lupus International Collaborating Clinics (SLICC) recently validated an alternative system. This study examined large cohorts of subjects with SLE and incomplete lupus erythematosus (ILE) to compare the impact of ACR and SLICC criteria.

Methods

Medical records of subjects in the Lupus Family Registry and Repository were reviewed for documentation of 1997 ACR classification criteria, SLICC classification criteria and medication usage. Autoantibodies were assessed by indirect immunofluorescence (ANA, antidouble-stranded DNA), precipitin (Sm) and ELISA (anticardiolipin). Other relevant autoantibodies were detected by precipitin and with a bead-based multiplex assay.

Results

Of 3575 subjects classified with SLE under at least one system, 3312 (92.6%) were classified as SLE by both systems (SLE^both), 85 only by ACR criteria (SLE^ACR-only) and 178 only by SLICC criteria (SLE^SLICC-only). Of 440 subjects meeting 3 ACR criteria, 33.9% (149/440) were SLE^SLICC-only, while 66.1% (n=291, designated ILE) did not meet the SLICC classification criteria. Under the SLICC system, the complement criterion and the individual autoantibody criteria enabled SLE classification of SLE^SLICC-only subjects, while SLE^ACR-only subjects failed to meet SLICC classification due to the combined acute/subacute cutaneous criterion. The SLICC criteria classified more African-American subjects by the leucopenia/lymphopenia criterion than did ACR criteria. Compared with SLE^ACR-only subjects, SLE^SLICC-only subjects exhibited similar numbers of affected organ systems, rates of major organ system involvement (∼30%: pulmonary, cardiovascular, renal, neurological) and medication history.

Conclusions

The SLICC criteria classify more subjects with SLE than ACR criteria; however, individuals with incomplete lupus still exist under SLICC criteria. Subjects who gain SLE classification through SLICC criteria exhibit heterogeneous disease, including potential major organ involvement. These results provide supportive evidence that SLICC criteria may be more inclusive of SLE subjects for clinical studies.

A missense variant in NCF1 is associated with susceptibility to multiple autoimmune diseases

Systemic lupus erythematosus (SLE) is a heterogeneous autoimmune disease with a strong genetic component characterized by autoantibody production and a type I interferon signature. Here we report a missense variant (g.74779296G>A; p.Arg90His) in NCF1, encoding the p47^phox subunit of the phagocyte NADPH oxidase (NOX2), as the putative underlying causal variant that drives a strong SLE-associated signal detected by the Immunochip in the GTF2IRD1-GTF2I region at 7q11.23 with a complex genomic structure. We show that the p.Arg90His substitution, which is reported to cause reduced reactive oxygen species (ROS) production, predisposes to SLE (odds ratio (OR) = 3.47 in Asians (P_meta = 3.1 × 10^-104), OR = 2.61 in European Americans, OR = 2.02 in African Americans) and other autoimmune diseases, including primary Sjögren's syndrome (OR = 2.45 in Chinese, OR = 2.35 in European Americans) and rheumatoid arthritis (OR = 1.65 in Koreans). Additionally, decreased and increased copy numbers of NCF1 predispose to and protect against SLE, respectively. Our data highlight the pathogenic role of reduced NOX2-derived ROS levels in autoimmune diseases.

Genome-Wide Association Study in an Amerindian Ancestry Population Reveals Novel Systemic Lupus Erythematosus Risk Loci and the Role of European Admixture

OBJECTIVE

Systemic lupus erythematosus (SLE) is a chronic autoimmune disease with a strong genetic component. We undertook the present work to perform the first genome-wide association study on individuals from the Americas who are enriched for Native American heritage.

METHODS

We analyzed 3,710 individuals from the US and 4 countries of Latin America who were diagnosed as having SLE, and healthy controls. Samples were genotyped with HumanOmni1 BeadChip. Data on out-of-study controls genotyped with HumanOmni2.5 were also included. Statistical analyses were performed using SNPtest and SNPGWA. Data were adjusted for genomic control and false discovery rate. Imputation was performed using Impute2 and, for classic HLA alleles, HiBag. Odds ratios (ORs) and 95% confidence intervals (95% CIs) were calculated.

RESULTS

The IRF5-TNPO3 region showed the strongest association and largest OR for SLE (rs10488631: genomic control-adjusted P [Pgcadj ] = 2.61 × 10(-29), OR 2.12 [95% CI 1.88-2.39]), followed by HLA class II on the DQA2-DQB1 loci (rs9275572: Pgcadj  = 1.11 × 10(-16), OR 1.62 [95% CI 1.46-1.80] and rs9271366: Pgcadj  = 6.46 × 10(-12), OR 2.06 [95% CI 1.71-2.50]). Other known SLE loci found to be associated in this population were ITGAM, STAT4, TNIP1, NCF2, and IRAK1. We identified a novel locus on 10q24.33 (rs4917385: Pgcadj  = 1.39 × 10(-8)) with an expression quantitative trait locus (eQTL) effect (Peqtl  = 8.0 × 10(-37) at USMG5/miR1307), and several new suggestive loci. SLE risk loci previously identified in Europeans and Asians were corroborated. Local ancestry estimation showed that the HLA allele risk contribution is of European ancestral origin. Imputation of HLA alleles suggested that autochthonous Native American haplotypes provide protection against development of SLE.

CONCLUSION

Our results demonstrate that studying admixed populations provides new insights in the delineation of the genetic architecture that underlies autoimmune and complex diseases.

Brief Report: Patients With Primary Sjögren's Syndrome Who Are Positive for Autoantibodies to Tripartite Motif-Containing Protein 38 Show Greater Disease Severity

OBJECTIVE

Autoantibodies reactive with Ro52 (tripartite motif-containing protein 21 [TRIM21]) are detected in 70% of patients with primary Sjögren's syndrome (SS). TRIM21 belongs to a 34-member C-IV family of TRIM proteins. Although autoantibodies against other TRIM proteins within the C-IV family have been detected in the sera of patients with primary SS, their clinical relevance remains unclear. This study was undertaken to investigate the frequency of anti-TRIM38 in patients with primary SS and evaluate its association with various clinical measures of the disease.

METHODS

Serum samples from patients with primary SS (n = 235) and controls (n = 50) were analyzed for reactivity with in vitro-transcribed and -translated (35) S-methionine-labeled TRIM38 protein. The associations of anti-TRIM38 with various laboratory and clinical measures of primary SS were evaluated. Reactivity of anti-TRIM38 with different structural domains of TRIM38 was analyzed. Affinity-purified anti-TRIM38 antibodies were used to immunoprecipitate TRIM21.

RESULTS

TRIM38-reactive autoantibodies were detected in the sera of 24 of the 235 patients with primary SS and 2 of the 50 controls. Anti-TRIM38 positivity was significantly associated with the presence of anti-Ro60, anti-Ro52, anti-La, rheumatoid factor, and hypergammaglobulinemia. Clinically, anti-TRIM38 was associated with significantly higher ocular surface staining scores, lower Schirmer's test scores, and minor labial salivary gland biopsy focus scores of ≥3.0. Anti-TRIM38 antibodies mainly recognized the cortactin-binding protein 2 (CortBP-2; amino acids 128-238) and the B30.2/SPRY (amino acids 268-465) domains on TRIM38. Affinity-purified antibodies to TRIM38-CortBP-2 and TRIM38-B30.2/SPRY domains reacted with TRIM21.

CONCLUSION

Our data demonstrate that anti-TRIM38 specificity arising in a subset of patients with primary SS is associated with increased severity of the disease.

Genetics in Sjögren Syndrome

The genes associated with Sjögren syndrome (SS) can be assigned to the NF-kB pathway, the IFN signaling pathway, lymphocyte signaling, and antigen presentation. The frequencies of risk variants show they are common with modest genetic effects. The strongest genetic association outside the human leukocyte antigen region is in IRF5, a gene relevant in the IFN signaling pathway and for B cell differentiation. Although no association has been found with the NF-kB gene itself, associations in TNFAIP3 and TNIP1 (both genome-wide significant), VCAM1 and IRAK1BP (both suggestive), point to genetic explanations for dysregulation of the NF-kB pathway in SS.

Single-cell analysis of glandular T cell receptors in Sjögren's syndrome

CD4⁺ T cells predominate in salivary gland (SG) inflammatory lesions in Sjögren's syndrome (SS). However, their antigen specificity, degree of clonal expansion, and relationship to clinical disease features remain unknown. We used multiplex reverse-transcriptase PCR to amplify paired T cell receptor α (TCRα) and β transcripts of single CD4⁺CD45RA^- T cells from SG and peripheral blood (PB) of 10 individuals with primary SS, 9 of whom shared the HLA DR3/DQ2 risk haplotype. TCRα and β sequences were obtained from a median of 91 SG and 107 PB cells per subject. The degree of clonal expansion and frequency of cells expressing two productively rearranged α genes were increased in SG versus PB. Expanded clones from SG exhibited complementary-determining region 3 (CDR3) sequence similarity both within and among subjects, suggesting antigenic selection and shared antigen recognition. CDR3 similarities were shared among expanded clones from individuals discordant for canonical Ro and La autoantibodies, suggesting recognition of alternative SG antigen(s). The extent of SG clonal expansion correlated with reduced saliva production and increased SG fibrosis, linking expanded SG T cells with glandular dysfunction. Knowledge of paired TCRα and β sequences enables further work toward identification of target antigens and development of novel therapies.

Identification of a Systemic Lupus Erythematosus Risk Locus Spanning ATG16L2, FCHSD2, and P2RY2 in Koreans

OBJECTIVE

Systemic lupus erythematosus (SLE) is a chronic autoimmune disorder whose etiology is incompletely understood, but likely involves environmental triggers in genetically susceptible individuals. Using an unbiased genome-wide association (GWA) scan and replication analysis, we sought to identify the genetic loci associated with SLE in a Korean population.

METHODS

A total of 1,174 SLE cases and 4,246 population controls from Korea were genotyped and analyzed with a GWA scan to identify single-nucleotide polymorphisms (SNPs) significantly associated with SLE, after strict quality control measures were applied. For select variants, replication of SLE risk loci was tested in an independent data set of 1,416 SLE cases and 1,145 population controls from Korea and China.

RESULTS

Eleven regions outside the HLA exceeded the genome-wide significance level (P = 5 × 10(-8) ). A novel SNP-SLE association was identified between FCHSD2 and P2RY2, peaking at rs11235667 (P = 1.03 × 10(-8) , odds ratio [OR] 0.59) on a 33-kb haplotype upstream of ATG16L2. In the independent replication data set, the SNP rs11235667 continued to show a significant association with SLE (replication meta-analysis P = 0.001, overall meta-analysis P = 6.67 × 10(-11) ; OR 0.63). Within the HLA region, the SNP-SLE association peaked in the class II region at rs116727542, with multiple independent effects observed in this region. Classic HLA allele imputation analysis identified HLA-DRB1*1501 and HLA-DQB1*0602, each highly correlated with one another, as most strongly associated with SLE. Ten previously established SLE risk loci were replicated: STAT1-STAT4, TNFSF4, TNFAIP3, IKZF1, HIP1, IRF5, BLK, WDFY4, ETS1, and IRAK1-MECP2. Of these loci, previously unreported, independent second risk effects of SNPs in TNFAIP3 and TNFSF4, as well as differences in the association with a putative causal variant in the WDFY4 region, were identified.

CONCLUSION

Further studies are needed to identify true SLE risk effects in other loci suggestive of a significant association, and to identify the causal variants in the regions of ATG16L2, FCHSD2, and P2RY2.

Blood relatives with incomplete lupus are distinguished from unaffected relatives and relatives with classified disease by altered autoantibodies, soluble inflammatory and regulatory mediators, and clinical features

Blood relatives (Rel) of lupus patients have increased risk of systemic lupus erythematosus (SLE). Some have autoantibodies and clinical features, but do not meet ≥ 4 ACR criteria needed to reach SLE classification (incomplete lupus, ILE). This study examined individuals enrolled in the Lupus Family Registry and Repository who only met 3 ACR classification criteria during medical record review. ILE patients (n=77) were matched to unaffected Rel and unrelated controls (Ctls), with a subset of ILE patients (n=55) matched to medical record-confirmed SLE patients. ILE patients had significantly higher SLE-specific screening questionnaire (CSQ) scores than unaffected Rel and Ctls (p≤0.0001). A number of soluble mediators that positively correlated with CSQ scores were highest in ILE patients, including SCF (p=0.0001), BLyS (p=0.0018), MCP-3 (p=0.0167), and TNFRI (p=0.0196), as well as ANA titer (p&lt;0.01) and number of lupus-associated autoantibodies (p&lt;0.01). Comparing ILE vs. SLE patients, SLE patients had higher rates of arthritis, serositis, and renal disease (p≤0.004), number of autoantibody specificities (p&lt;0.004), and levels of BLyS (p=0.0138), IL-2Rα (p=0.0201), IP-10 (p=0.0269), and TNFRII (p=0.0309). Yet, ILE patients had higher levels of the regulatory mediator TGF-β (p=0.0454). Identification of factors which discern relatives at increased risk of transitioning to classified SLE may be beneficial to curtail inflammatory damage.

Klinefelter's syndrome (47,XXY) is in excess among men with Sjögren's syndrome

Primary Sjögren's syndrome (pSS) has a strong female bias. We evaluated an X chromosome dose effect by analyzing 47,XXY (Klinefelter's syndrome, 1 in 500 live male births) among subjects with pSS. 47,XXY was determined by examination of fluorescence intensity of single nucleotide polymorphisms from the X and Y chromosomes. Among 136 pSS men there were 4 with 47,XXY. This was significantly different from healthy controls (1 of 1254 had 47,XXY, p=0.0012 by Fisher's exact test) as well men with rheumatoid arthritis (0 of 363 with 47,XXY), but not different compared to men with systemic lupus erythematosus (SLE) (4 of 136 versus 8 of 306, Fisher's exact test p=NS). These results are consistent with the hypothesis that the number of X chromosomes is critical for the female bias of pSS, a property that may be shared with SLE but not RA.

Previous diagnosis of Sjögren's Syndrome as rheumatoid arthritis or systemic lupus erythematosus

OBJECTIVE

The diagnosis of SS is often difficult and many patients are symptomatic for years with other diagnoses before confirmation of SS. Our aim was to determine whether overlapping clinical and serologic features with RA and SLE may in part drive the misdiagnoses.

METHODS

A total of 1175 sicca patients were evaluated in a multidisciplinary clinic and classified as having SS based on the American-European Consensus Group Criteria. They were interrogated for a past history of suspicion or diagnosis of RA, SLE or SSc. These diseases were confirmed or ruled out by applying the corresponding classification criteria if the patients responded affirmatively.

RESULTS

Of these, 524 (44.6%) subjects reported previous diagnosis or suspicion of RA, SLE or SSc, which was confirmed in 130 (24.8%) but excluded in 394 (75.2%) subjects. Of those previously diagnosed with another illness, 183 (34.9%) met the criteria for primary SS. RF was present in 70/191 patients with previous diagnosis of RA compared with 445/845 without a prior RA diagnosis (P = 3.38E-05), while 128/146 with a diagnosis of SLE had positive ANA compared with 622/881 without the diagnosis (P = 8.77E-06). Age also influenced former diagnoses: people with suspected RA were older than those without the diagnosis (P = 5.89E-06), while patients with SLE suspicion were younger (P = 0.0003). Interestingly, the previous diagnoses did not significantly delay a final classification of SS.

CONCLUSION

Among subjects classified as SS, the presence of a positive ANA or RF was associated with a previous, apparently erroneous diagnosis of SLE or RA, respectively.

Identification of a New Susceptibility Locus for Systemic Lupus Erythematosus on Chromosome 12 in Individuals of European Ancestry

OBJECTIVE

Genome-wide association studies (GWAS) in individuals of European ancestry identified a number of systemic lupus erythematosus (SLE) susceptibility loci using earlier versions of high-density genotyping platforms. Followup studies on suggestive GWAS regions using larger samples and more markers identified additional SLE loci in subjects of European descent. This multistage study was undertaken to identify novel SLE loci.

METHODS

In stage 1, we conducted a new GWAS of SLE in a North American case-control sample of subjects of European ancestry (n = 1,166) genotyped on Affymetrix Genome-Wide Human SNP Array 6.0. In stage 2, we further investigated top new suggestive GWAS hits by in silico evaluation and meta-analysis using an additional data set of subjects of European descent (>2,500 individuals), followed by replication of top meta-analysis findings in another data set of subjects of European descent (>10,000 individuals) in stage 3.

RESULTS

As expected, our GWAS revealed the most significant associations at the major histocompatibility complex locus (6p21), which easily surpassed the genome-wide significance threshold (P < 5 × 10(-8)). Several other SLE signals/loci previously implicated in Caucasians and/or Asians were also confirmed in the stage 1 discovery sample, and the strongest signals were observed at 2q32/STAT4 (P = 3.6 × 10(-7)) and at 8p23/BLK (P = 8.1 × 10(-6)). Stage 2 meta-analyses identified a new genome-wide significant SLE locus at 12q12 (meta P = 3.1 × 10(-8)), which was replicated in stage 3.

CONCLUSION

Our multistage study identified and replicated a new SLE locus that warrants further followup in additional studies. Publicly available databases suggest that this newly identified SLE signal falls within a functionally relevant genomic region and near biologically important genes.

Detrimental effects of duplicate reads and low complexity regions on RNA- and ChIP-seq data

Background

Adapter trimming and removal of duplicate reads are common practices in next-generation sequencing pipelines. Sequencing reads ambiguously mapped to repetitive and low complexity regions can also be problematic for accurate assessment of the biological signal, yet their impact on sequencing data has not received much attention. We investigate how trimming the adapters, removing duplicates, and filtering out reads overlapping low complexity regions influence the significance of biological signal in RNA- and ChIP-seq experiments.

Methods

We assessed the effect of data processing steps on the alignment statistics and the functional enrichment analysis results of RNA- and ChIP-seq data. We compared differentially processed RNA-seq data with matching microarray data on the same patient samples to determine whether changes in pre-processing improved correlation between the two. We have developed a simple tool to remove low complexity regions, RepeatSoaker, available at https://github.com/mdozmorov/RepeatSoaker, and tested its effect on the alignment statistics and the results of the enrichment analyses.

Results

Both adapter trimming and duplicate removal moderately improved the strength of biological signals in RNA-seq and ChIP-seq data. Aggressive filtering of reads overlapping with low complexity regions, as defined by RepeatMasker, further improved the strength of biological signals, and the correlation between RNA-seq and microarray gene expression data.

Conclusions

Adapter trimming and duplicates removal, coupled with filtering out reads overlapping low complexity regions, is shown to increase the quality and reliability of detecting biological signals in RNA-seq and ChIP-seq data.

Germline variation of TNFAIP3 in primary Sjögren's syndrome-associated lymphoma

Background and objective

A germline and coding polymorphism (rs2230926) of TNFAIP3 (A20), a central gatekeeper of nuclear factor-kappa B (NF-kB) activation, was recently found associated with primary Sjögren's syndrome (pSS)-associated lymphoma in a French cohort. We aimed to replicate this association.

Patients and methods

The rs2230926 polymorphism was genotyped in cases and controls of European ancestry from two independent cohorts from UK and France. Case control association tests were performed (Fisher's test) in the two cohorts, followed by a meta-analysis of the two cohorts.

Results

The UK cohort included 308 controls and 590 patients with pSS including 31 with a history of lymphoma. The French cohort consisted of 448 controls and 589 patients with pSS including 47 with lymphoma. In both cohorts, the rs2230926 missense polymorphism was not associated with pSS. However, in the UK cohort, the rs2230926G variant was significantly associated with pSS-associated lymphoma (OR=2.74, 95% CI (1.07 to 7.03), p=0.0423, compared with patients with pSS without lymphoma, and OR=3.12, 95% CI (1.16 to 8.41), p=0.0314, compared with healthy controls) as observed in the French cohort. The meta-analysis of the two cohorts confirmed these results (OR=2.48, 95% CI (1.87 to 3.28) p=0.0037 and OR=2.60, 95% CI (1.91 to 3.53) p=0.0031, respectively).

Conclusions

This study confirms the role of A20 impairment in pSS-associated lymphoma. Subtle germline abnormalities of genes leading to impaired control of NF-kB activation in B cells continuously stimulated by autoimmunity enhance the risk of lymphoma.

EULAR Sjögren's syndrome disease activity index (ESSDAI): a user guide

The EULAR Sjögren's syndrome (SS) disease activity index (ESSDAI) is a systemic disease activity index that was designed to measure disease activity in patients with primary SS. With the growing use of the ESSDAI, some domains appear to be more challenging to rate than others. The ESSDAI is now in use as a gold standard to measure disease activity in clinical studies, and as an outcome measure, even a primary outcome measure, in current randomised clinical trials. Therefore, ensuring an accurate and reproducible rating of each domain, by providing a more detailed definition of each domain, has emerged as an urgent need. The purpose of the present article is to provide a user guide for the ESSDAI. This guide provides definitions and precisions on the rating of each domain. It also includes some minor improvement of the score to integrate advance in knowledge of disease manifestations. This user guide may help clinicians to use the ESSDAI, and increase the reliability of rating and consequently of the ability to detect true changes over time. This better appraisal of ESSDAI items, along with the recent definition of disease activity levels and minimal clinically important change, will improve the assessment of patients with primary SS and facilitate the demonstration of effectiveness of treatment for patients with primary SS.

Antibody-secreting cell specificity in labial salivary glands reflects the clinical presentation and serology in patients with Sjögren's syndrome

OBJECTIVE

The serologic hallmark of primary Sjögren's syndrome (SS) is the presence of IgG antibodies specific for Ro (SSA) and La (SSB). The molecular characteristics of gland-derived B cells at the site of primary SS inflammation have been described previously; however, parallels between glandular antibody-secreting cells (ASCs) and serologic antibody specificities have not been evaluated. We used recombinant monoclonal antibody (mAb) technology to study the specificities of salivary gland (SG)-derived ASCs, evaluate their molecular characteristics, and identify IgG antibody specificity.

METHODS

Human antibodies were generated from glandular IgG ASCs. Heavy chain and light chain use and immunoglobulin subclass were analyzed by sequencing. Enzyme-linked immunosorbent assay, indirect immunofluorescence, enzyme immunoassay, and (35) S-labeled protein immunoprecipitation analysis were used to determine antibody specificity.

RESULTS

Evaluation of single ASCs in SG biopsy specimens from a patient with primary SS and a patient with SS and overlapping systemic lupus erythematosus revealed significant concordance between serum autoantibody and glandular ASC specificities. Gland-derived ASC heavy chains and light chains were extensively somatically hypermutated, which is indicative of antigen-driven responses. Specifically, we produced the first fully human mAb derived from SGs.

CONCLUSION

In patients with SS, the SGs are a site for the production of antibodies that extend beyond the canonical Ro and/or La SS specificities. Glandular antibody production strongly reflected the serologic humoral response in the 2 patients whom we studied.

Interaction between innate immunity and Ro52-induced antibody causes Sjögren's syndrome-like disorder in mice

OBJECTIVES

Autoantibodies reactive with Ro52 are often found in sera of patients with Sjögren's syndrome (SS). This study was undertaken to investigate the role of Ro52-induced immune responses in pathogenesis of SS.

METHODS

New Zealand Mixed (NZM) 2758 mice were immunised with Ro52 in alum adjuvant. Control mice were immunised either with maltose-binding protein or injected with alum alone. Mice were monitored for anti-Ro52 antibody, sialoadenitis and pilocarpine-induced salivation. Antibody binding to salivary gland (SG) cells was analysed in vivo and in vitro by immunofluorescence. Sera from immunised mice were passively transferred into untreated or alum injected NZM2758 mice.

RESULTS

By day 30 post-immunisation, Ro52 immunised mice generated immunoprecipitating anti-Ro52 antibodies and they had the maximum drop in saliva production. Both Ro52 immunised and control mice showed evidence of mild sialoadenitis. However, only Ro52 immunised mice had antibody deposition in their SG. Passive transfer of Ro52-immune sera induced SG dysfunction in recipient mice, only if the recipients were primed with alum. In vitro, antibodies from Ro52-immune sera were internalised by a SG cell line and this uptake was inhibited by cytochalasin D treatment.

CONCLUSIONS

Our data show for the first time that antibodies induced by Ro52 are capable of inducing SG dysfunction, and that this phenomenon is dependent on the activation of innate immunity. The mouse model described in this study implies that autoantibody deposition in the SG might be an important step in the induction of xerostomia and pathogenesis of SS.

Genetic analysis of the pathogenic molecular sub-phenotype interferon-alpha identifies multiple novel loci involved in systemic lupus erythematosus

Systemic Lupus Erythematosus (SLE) is a chronic autoimmune disorder characterized by inflammation of multiple organ systems and dysregulated interferon responses. SLE is both genetically and phenotypically heterogeneous, greatly reducing the power of case-control studies in SLE. Elevated circulating interferon alpha (IFN-α) is a stable, heritable trait in SLE, which has been implicated in primary disease pathogenesis. 40–50% of patients have high IFN-α, and high levels correspond with clinical differences. To study genetic heterogeneity in SLE, we performed a case-case study comparing patients with high vs. low IFN-α in over 1550 SLE cases, including GWAS and replication cohorts. In meta-analysis, the top associations in European ancestry were PRKG1 rs7897633 (P_Meta=2.75 × 10⁻⁸) and PNP rs1049564 (P_Meta=1.24 × 10⁻⁷). We also found evidence for cross-ancestral background associations with the ANKRD44 and PLEKHF2 loci. These loci have not been previously identified in case-control SLE genetic studies. Bioinformatic analyses implicated these loci functionally in dendritic cells and natural killer cells, both of which are involved in IFN-α production in SLE. As case-control studies of heterogeneous diseases reach a limit of feasibility with respect to subject number and detectable effect size, the study of informative pathogenic subphenotypes becomes an attractive strategy for genetic discovery in complex disease.