Association between systemic lupus erythematosus and myasthenia gravis: A population-based National Study.

BACKGROUND: Systemic lupus erythematosus (SLE) and myasthenia gravis (MG) are autoimmune diseases. Previous case reports and case series suggest an association may exist between these diseases, as well as an increased risk of SLE after thymectomy for MG. We undertook this study to determine whether SLE and MG were associated in large cohorts. METHODS: We searched the IBM Watson Health Explorys platform and the Department of Veterans Affairs Million Veteran Program (MVP) database for diagnoses of SLE and MG. In addition, we examined subjects enrolled in the Lupus Family Registry and Repository (LFRR) as well as controls for a diagnosis of MG. RESULTS: Among 59,780,210 individuals captured in Explorys, there were 25,750 with MG and 65,370 with SLE. 370 subjects had both. Those with MG were >10 times more likely to have SLE than those without MG. Those with both diseases were more likely to be women, African American, and at a younger age than MG subjects without SLE. In addition, the MG patients who underwent thymectomy had an increased risk of SLE compared to MG patients who had not undergone thymectomy (OR 3.11, 95% CI: 2.12 to 4.55). Autoimmune diseases such as pernicious anemia and miscellaneous comorbidities such as chronic kidney disease were significantly more common in MG patients who developed SLE. In the MVP, SLE and MG were also significantly associated. Association of SLE and MG in a large SLE cohort with rigorous SLE classification confirmed the association of SLE with MG at a similar level. CONCLUSION: While the number of patients with both MG and SLE is small, SLE and MG are strongly associated together in very large databases and a large SLE cohort.

Lupus enhancer risk variant causes dysregulation of IRF8 through cooperative lncRNA and DNA methylation machinery.

Despite strong evidence that human genetic variants affect the expression of many key transcription factors involved in autoimmune diseases, establishing biological links between non-coding risk variants and the gene targets they regulate remains a considerable challenge. Here, we combine genetic, epigenomic, and CRISPR activation approaches to screen for functional variants that regulate IRF8 expression. We demonstrate that the locus containing rs2280381 is a cell-type-specific enhancer for IRF8 that spatially interacts with the IRF8 promoter. Further, rs2280381 mediates IRF8 expression through enhancer RNA AC092723.1, which recruits TET1 to the IRF8 promoter regulating IRF8 expression by affecting methylation levels. The alleles of rs2280381 modulate PU.1 binding and chromatin state to regulate AC092723.1 and IRF8 expression differentially. Our work illustrates an integrative strategy to define functional genetic variants that regulate the expression of critical genes in autoimmune diseases and decipher the mechanisms underlying the dysregulation of IRF8 expression mediated by lupus risk variants.

A High Prevalence of Anti-EBNA1 Heteroantibodies in Systemic Lupus Erythematosus (SLE) Supports Anti-EBNA1 as an Origin for SLE Autoantibodies.

BACKGROUND: That Epstein-Barr virus (EBV) infection is associated with systemic lupus erythematosus (SLE) is established. The challenge is to explain mechanistic roles EBV has in SLE pathogenesis. Previous studies identify four examples of autoantibody cross-reactions between SLE autoantigens and Epstein-Barr nuclear antigen 1 (EBNA1). For two of these examples, the earliest detected autoantibody specifically cross-reacts with EBNA1; thereby, defined EBNA1 epitopes induce a robust autoantibody response in animals. These results suggest that the autoantibodies initiating the process leading to SLE may emerge from the anti-EBNA1 heteroimmune response. If this hypothesis is true, then anti-EBNA1 responses would be more frequent in EBV-infected SLE patients than in EBV-infected controls. We tested this prediction. METHODS: We evaluated published East Asian data by selecting those with a positive anti-viral capsid antigen (VCA) antibody immunoglobulin G (IgG) test and determining whether anti-EBNA1 was more common among the EBV-infected SLE cases than among matched EBV-infected controls with conditional logistic regression analysis. RESULTS: All the qualifying SLE patients (100%) in this dataset were EBV-infected compared to age- and sex-matched controls (92.2%) [odds ratio (OR) = 28.6, 95% CI 6.4-∞, p = 8.83 × 10-8], confirming the known close association of EBV infection with SLE. Furthermore, virtually all the SLE cases have both anti-VCA IgG and anti-EBNA1 IgG antibodies [124 of 125 (99.2%)], which are more frequently present than in age- and sex-matched EBV-infected controls [232 of 250 (93.2%)] (OR = 9.7, 95% CI 1.5-414, p = 0.0078) for an 89.7% SLE attributable risk from anti-EBNA1, which is in addition to the 100% SLE risk attributable to EBV infection in these data. CONCLUSIONS: The association of EBV infection with SLE is reconfirmed. The prediction that anti-EBNA1 is more frequent in these SLE cases than in EBV-infected controls is true, consistent with the hypothesis that anti-EBNA1 contributes to SLE. This second EBV-dependent risk factor is consistent with a molecular mimicry model for the generation of SLE, starting with EBV infection, progressing to anti-EBNA1 response; then molecular mimicry leads to anti-EBNA1 antibodies cross-reacting with an SLE autoantigen, causing autoantibody epitope spreading, and culminating in clinical SLE. These results support the anti-EBNA1 heteroimmune response being a foundation from which pathogenic SLE autoimmunity emerges.

Antibody Responses to Epstein-Barr Virus in the Preclinical Period of Rheumatoid Arthritis Suggest the Presence of Increased Viral Reactivation Cycles.

OBJECTIVE: Patients with established rheumatoid arthritis (RA) demonstrate altered immune responses to Epstein-Barr virus (EBV), but the presence and roles of EBV have not been fully explored during the pre-clinical disease period. This study was undertaken to determine if EBV infection, as evidenced by an altered anti-EBV antibody response, either plays an important role in driving the development of RA or is a result of expanded RA-related autoimmunity. METHODS: A total of 83 subjects with RA according to the 1987 American College of Rheumatology (ACR) criteria and 83 age-, sex-, and race-matched control subjects without RA were included in our study. We collected sera from RA subjects and matched controls during the pre-RA and post-RA diagnosis periods and tested the sera for the presence of 5 anti-EBV antibodies (anti-EBV nuclear antigen 1 IgG isotype, anti-viral capsid antigen [anti-VCA] isotypes IgG and IgA, and anti-early antigen [EA] isotypes IgG and IgA), 7 RA-related autoantibodies (rheumatoid factor measured by nephelometry [RF-Neph] as well as isotype-specific IgA-RF, IgM-RF, and IgG-RF, and anti-cyclic citrullinated peptide [anti-CCP] antibodies, including anti-CCP2, anti-CCP3, and anti-CCP3.1), 22 cytokines/chemokines, 36 individual anti-citrullinated protein antibodies, and IgG-cytomegalovirus (CMV) antibodies. Random forest classification, mixed modeling, and joint mixed modeling were used to determine differences in anti-EBV antibody levels between RA subjects and controls. RESULTS: Random forest analysis identified the presence of preclinical EBV antibodies in the serum that differentiated RA subjects from controls without RA. Specifically, IgG-EA antibody levels were higher in RA subjects (mean ± SD 0.82 ± 0.72 international standardized ratio [ISR]) compared to controls (mean ± SD 0.49 ± 0.28 ISR). In subjects with RA, elevated serum IgG-EA levels in the preclinical period before seroconversion were significantly correlated with increased serum IgM-RF levels (P = 0.007), whereas this correlation was not seen in control subjects without RA (P = 0.15). IgG-CMV antibody levels did not differ between groups. CONCLUSION: Subjects whose serum IgG-EA antibody levels are elevated in the preclinical period will eventually develop RA, which suggests that EBV reactivation cycles are increased during the preclinical period of RA. A combination of RF and EBV reactivation may play an important role in the development of RA.

Pleiotropy of systemic lupus erythematosus risk alleles and cardiometabolic disorders: A phenome-wide association study and inverse-variance weighted meta-analysis.

OBJECTIVES: To test the hypothesis that genetic predisposition to systemic lupus erythematosus (SLE) increases the risk of cardiometabolic disorders. METHODS: Using 41 single nucleotide polymorphisms (SNPs) associated with SLE, we calculated a weighted genetic risk score (wGRS) for SLE. In a large biobank we tested the association between this wGRS and 9 cardiometabolic phenotypes previously associated with SLE: atrial fibrillation, ischemic stroke, coronary artery disease, type 1 and type 2 diabetes, obesity, chronic kidney disease, hypertension, and hypercholesterolemia. Additionally, we performed a phenome-wide association analysis (pheWAS) to discover novel clinical associations with a genetic predisposition to SLE. Findings were replicated in the Electronic Medical Records and Genomics (eMERGE) Network. To further define the association between SLE-related risk alleles and the selected cardiometabolic phenotypes, we performed an inverse variance weighted regression (IVWR) meta-analysis. RESULTS: The wGRS for SLE was calculated in 74,759 individuals of European ancestry. Among the pre-selected phenotypes, the wGRS was significantly associated with type 1 diabetes (OR [95%CI] =1.11 [1.06, 1.17], P-value = 1.05x10-5). In the PheWAS, the wGRS was associated with several autoimmune phenotypes, kidney disorders, and skin neoplasm; but only the associations with autoimmune phenotypes were replicated. In the IVWR meta-analysis, SLE-related risk alleles were nominally associated with type 1 diabetes (P = 0.048) but the associations were heterogeneous and did not meet the adjusted significance threshold. CONCLUSION: A weighted GRS for SLE was associated with an increased risk of several autoimmune-related phenotypes including type I diabetes but not with cardiometabolic disorders.

Global discovery of lupus genetic risk variant allelic enhancer activity.

Genome-wide association studies of Systemic Lupus Erythematosus (SLE) nominate 3073 genetic variants at 91 risk loci. To systematically screen these variants for allelic transcriptional enhancer activity, we construct a massively parallel reporter assay (MPRA) library comprising 12,396 DNA oligonucleotides containing the genomic context around every allele of each SLE variant. Transfection into the Epstein-Barr virus-transformed B cell line GM12878 reveals 482 variants with enhancer activity, with 51 variants showing genotype-dependent (allelic) enhancer activity at 27 risk loci. Comparison of MPRA results in GM12878 and Jurkat T cell lines highlights shared and unique allelic transcriptional regulatory mechanisms at SLE risk loci. In-depth analysis of allelic transcription factor (TF) binding at and around allelic variants identifies one class of TFs whose DNA-binding motif tends to be directly altered by the risk variant and a second class of TFs that bind allelically without direct alteration of their motif by the variant. Collectively, our approach provides a blueprint for the discovery of allelic gene regulation at risk loci for any disease and offers insight into the transcriptional regulatory mechanisms underlying SLE.

Patients with Proliferative Lupus Nephritis Have Autoantibodies That React to Moesin and Demonstrate Increased Glomerular Moesin Expression.

Kidney involvement in systemic lupus erythematosus (SLE)-termed lupus nephritis (LN)-is a severe manifestation of SLE that can lead to end-stage kidney disease (ESKD). LN is characterized by immune complex deposition and inflammation in the glomerulus. We tested the hypothesis that autoantibodies targeting podocyte and glomerular cell proteins contribute to the development of immune complex formation in LN. We used Western blotting with SLE sera from patients with and without LN to identify target antigens in human glomerular and cultured human-derived podocyte membrane proteins. Using liquid chromatography-tandem mass spectrometry (LC-MS/MS), we identified the proteins in the gel regions corresponding to reactive bands observed with sera from LN patients. We identified 102 proteins that were present in both the podocyte and glomerular samples. We identified 10 high-probability candidates, including moesin, using bioinformatic analysis. Confirmation of moesin as a target antigen was conducted using immunohistochemical analysis (IHC) of kidney biopsy tissue and enzyme-linked immunosorbent assay (ELISA) to detect circulating antibodies. By IHC, biopsies from patients with proliferative lupus nephritis (PLN, class III/IV) demonstrated significantly increased glomerular expression of moesin (p < 0.01). By ELISA, patients with proliferative LN demonstrated significantly increased antibodies against moesin (p < 0.01). This suggests that moesin is a target glomerular antigen in lupus nephritis.

SLE non-coding genetic risk variant determines the epigenetic dysfunction of an immune cell specific enhancer that controls disease-critical microRNA expression.

Since most variants that impact polygenic disease phenotypes localize to non-coding genomic regions, understanding the consequences of regulatory element variants will advance understanding of human disease mechanisms. Here, we report that the systemic lupus erythematosus (SLE) risk variant rs2431697 as likely causal for SLE through disruption of a regulatory element, modulating miR-146a expression. Using epigenomic analysis, genome-editing and 3D chromatin structure analysis, we show that rs2431697 tags a cell-type dependent distal enhancer specific for miR-146a that physically interacts with the miR-146a promoter. NF-kB binds the disease protective allele in a sequence-specific manner, increasing expression of this immunoregulatory microRNA. Finally, CRISPR activation-based modulation of this enhancer in the PBMCs of SLE patients attenuates type I interferon pathway activation by increasing miR-146a expression. Our work provides a strategy to define non-coding RNA functional regulatory elements using disease-associated variants and provides mechanistic links between autoimmune disease risk genetic variation and disease etiology.

Genome-wide Modeling of Polygenic Risk Score in Colorectal Cancer Risk.

Accurate colorectal cancer (CRC) risk prediction models are critical for identifying individuals at low and high risk of developing CRC, as they can then be offered targeted screening and interventions to address their risks of developing disease (if they are in a high-risk group) and avoid unnecessary screening and interventions (if they are in a low-risk group). As it is likely that thousands of genetic variants contribute to CRC risk, it is clinically important to investigate whether these genetic variants can be used jointly for CRC risk prediction. In this paper, we derived and compared different approaches to generating predictive polygenic risk scores (PRS) from genome-wide association studies (GWASs) including 55,105 CRC-affected case subjects and 65,079 control subjects of European ancestry. We built the PRS in three ways, using (1) 140 previously identified and validated CRC loci; (2) SNP selection based on linkage disequilibrium (LD) clumping followed by machine-learning approaches; and (3) LDpred, a Bayesian approach for genome-wide risk prediction. We tested the PRS in an independent cohort of 101,987 individuals with 1,699 CRC-affected case subjects. The discriminatory accuracy, calculated by the age- and sex-adjusted area under the receiver operating characteristics curve (AUC), was highest for the LDpred-derived PRS (AUC = 0.654) including nearly 1.2 M genetic variants (the proportion of causal genetic variants for CRC assumed to be 0.003), whereas the PRS of the 140 known variants identified from GWASs had the lowest AUC (AUC = 0.629). Based on the LDpred-derived PRS, we are able to identify 30% of individuals without a family history as having risk for CRC similar to those with a family history of CRC, whereas the PRS based on known GWAS variants identified only top 10% as having a similar relative risk. About 90% of these individuals have no family history and would have been considered average risk under current screening guidelines, but might benefit from earlier screening. The developed PRS offers a way for risk-stratified CRC screening and other targeted interventions.

Pleiotropy in the Genetic Predisposition to Rheumatoid Arthritis: A Phenome-Wide Association Study and Inverse Variance-Weighted Meta-Analysis.

OBJECTIVE: This study was undertaken to investigate the hypothesis that a genetic predisposition toward rheumatoid arthritis (RA) increases the risk of 10 cardiometabolic and autoimmune disorders previously associated with RA in epidemiologic studies, and to define new genetic pleiotropy present in RA. METHODS: Two approaches were used to test our hypothesis. First, we constructed a weighted genetic risk score (wGRS) and then examined its association with 10 prespecified disorders. Additionally, a phenome-wide association study (PheWAS) was carried out to identify potential new associations. Second, inverse variance-weighted regression (IVWR) meta-analysis was used to characterize the association between genetic susceptibility to RA and the prespecified disorders, with the results expressed as odds ratios (ORs) and 95% confidence intervals (95% CIs). RESULTS: The wGRS for RA was significantly associated with type 1 diabetes mellitus (DM) (OR 1.10 [95% CI 1.04-1.16]; P = 9.82 × 10-4 ) and multiple sclerosis (OR 0.82 [95% CI 0.77-0.88]; P = 1.73 × 10-8 ), but not with other cardiometabolic phenotypes. In the PheWAS, wGRS was also associated with an increased risk of several autoimmune phenotypes including RA, thyroiditis, and systemic sclerosis, and with a decreased risk of demyelinating disorders. In the IVWR meta-analyses, RA was significantly associated with an increased risk of type 1 DM (P = 1.15 × 10-14 ), with evidence of horizontal pleiotropy (Mendelian Randomization-Egger intercept estimate P = 0.001) likely driven by rs2476601, a PTPN22 variant. The association between type 1 DM and RA remained significant (P = 9.53 × 10-9 ) after excluding rs2476601, with no evidence of horizontal pleiotropy (intercept estimate P = 0.939). RA was also significantly associated with type 2 DM and C-reactive protein levels. These associations were driven by variation in the major histocompatibility complex region. CONCLUSION: This study presents evidence of pleiotropy between the genetic predisposition to RA and associated phenotypes found in other autoimmune and cardiometabolic disorders, including type 1 DM.

A Polygenic and Phenotypic Risk Prediction for Polycystic Ovary Syndrome Evaluated by Phenome-Wide Association Studies.

CONTEXT: As many as 75% of patients with polycystic ovary syndrome (PCOS) are estimated to be unidentified in clinical practice. OBJECTIVE: Utilizing polygenic risk prediction, we aim to identify the phenome-wide comorbidity patterns characteristic of PCOS to improve accurate diagnosis and preventive treatment. DESIGN, PATIENTS, AND METHODS: Leveraging the electronic health records (EHRs) of 124 852 individuals, we developed a PCOS risk prediction algorithm by combining polygenic risk scores (PRS) with PCOS component phenotypes into a polygenic and phenotypic risk score (PPRS). We evaluated its predictive capability across different ancestries and perform a PRS-based phenome-wide association study (PheWAS) to assess the phenomic expression of the heightened risk of PCOS. RESULTS: The integrated polygenic prediction improved the average performance (pseudo-R2) for PCOS detection by 0.228 (61.5-fold), 0.224 (58.8-fold), 0.211 (57.0-fold) over the null model across European, African, and multi-ancestry participants respectively. The subsequent PRS-powered PheWAS identified a high level of shared biology between PCOS and a range of metabolic and endocrine outcomes, especially with obesity and diabetes: "morbid obesity", "type 2 diabetes", "hypercholesterolemia", "disorders of lipid metabolism", "hypertension", and "sleep apnea" reaching phenome-wide significance. CONCLUSIONS: Our study has expanded the methodological utility of PRS in patient stratification and risk prediction, especially in a multifactorial condition like PCOS, across different genetic origins. By utilizing the individual genome-phenome data available from the EHR, our approach also demonstrates that polygenic prediction by PRS can provide valuable opportunities to discover the pleiotropic phenomic network associated with PCOS pathogenesis.

Epstein Barr virus nuclear antigen 1 (EBNA-1) peptides recognized by adult multiple sclerosis patient sera induce neurologic symptoms in a murine model.

Multiple sclerosis (MS) is an autoimmune demyelinating disease with progressive neurodegeneration and complex etiology likely involving genetic and environmental factors. MS has been associated with Epstein Barr virus (EBV) infection, with patients often showing enhanced responses to EBV antigens. To determine whether abnormal EBV nuclear antigen-1 (EBNA-1) humoral immunity can serve as an initiator of autoimmune responses in MS, we investigated the fine specificities of the humoral immune response against EBNA-1 in MS patients using solid phase epitope mapping. Antibodies from MS patients recognized an EBNA-1 epitope spanning amino acids 411-426, previously unknown to be recognized specifically by untreated MS patients. Antibodies against this epitope cross-reacted to myelin basic protein (MBP). Furthermore, animals immunized with this EBNA-1 polypeptide mounted a response against MBP and developed signs of experimental autoimmune encephalitis (EAE). These data support a link between MS and EBV through antibodies that cross-react between EBV proteins and the MBP autoantigen.

Association of Epstein-Barr virus serological reactivation with transitioning to systemic lupus erythematosus in at-risk individuals.

OBJECTIVE: Systemic lupus erythematosus (SLE) is a systemic autoimmune disease with unknown aetiology. Epstein-Barr virus (EBV) is an environmental factor associated with SLE. EBV maintains latency in B cells with frequent reactivation measured by antibodies against viral capsid antigen (VCA) and early antigen (EA). In this study, we determined whether EBV reactivation and single nucleotide polymorphisms (SNPs) in EBV-associated host genes are associated with SLE transition. METHODS: SLE patient relatives (n=436) who did not have SLE at baseline were recontacted after 6.3 (±3.9) years and evaluated for interim transitioning to SLE (≥4 cumulative American College of Rheumatology criteria); 56 (13%) transitioned to SLE prior to the follow-up visit. At both visits, detailed demographic, environmental, clinical information and blood samples were obtained. Antibodies against viral antigens were measured by ELISA. SNPs in IL10, CR2, TNFAIP3 and CD40 genes were typed by ImmunoChip. Generalised estimating equations were used to test associations between viral antibody levels and transitioning to SLE. RESULTS: Mean baseline VCA IgG (4.879±1.797 vs 3.866±1.795, p=0.0003) and EA IgG (1.192±1.113 vs 0.7774±0.8484, p=0.0236) levels were higher in transitioned compared with autoantibody negative non-transitioned relatives. Increased VCA IgG and EA IgG were associated with transitioning to SLE (OR 1.28 95% CI 1.07 to 1.53, p=0.007, OR 1.43 95% CI 1.06 to 1.93, p=0.02, respectively). Significant interactions were observed between CD40 variant rs48100485 and VCA IgG levels and IL10 variant rs3024493 and VCA IgA levels in transitioning to SLE. CONCLUSION: Heightened serologic reactivation of EBV increases the probability of transitioning to SLE in unaffected SLE relatives.

Trans-Ethnic Mapping of BANK1 Identifies Two Independent SLE-Risk Linkage Groups Enriched for Co-Transcriptional Splicing Marks.

BANK1 is a susceptibility gene for several systemic autoimmune diseases in several populations. Using the genome-wide association study (GWAS) data from Europeans (EUR) and African Americans (AA), we performed an extensive fine mapping of ankyrin repeats 1 (BANK1). To increase the SNP density, we used imputation followed by univariate and conditional analysis, combined with a haplotypic and expression quantitative trait locus (eQTL) analysis. The data from Europeans showed that the associated region was restricted to a minimal and dependent set of SNPs covering introns two and three, and exon two. In AA, the signal found in the Europeans was split into two independent effects. All of the major risk associated SNPs were eQTLs, and the risks were associated with an increased BANK1 gene expression. Functional annotation analysis revealed the enrichment of repressive B cell epigenomic marks (EZH2 and H3K27me3) and a strong enrichment of splice junctions. Furthermore, one eQTL located in intron two, rs13106926, was found within the binding site for RUNX3, a transcriptional activator. These results connect the local genome topography, chromatin structure, and the regulatory landscape of BANK1 with co-transcriptional splicing of exon two. Our data defines a minimal set of risk associated eQTLs predicted to be involved in the expression of BANK1 modulated through epigenetic regulation and splicing. These findings allow us to suggest that the increased expression of BANK1 will have an impact on B-cell mediated disease pathways.

Whole-exome sequencing uncovers oxidoreductases DHTKD1 and OGDHL as linkers between mitochondrial dysfunction and eosinophilic esophagitis.

Eosinophilic esophagitis (EoE) is an allergic inflammatory esophageal disorder with a complex underlying genetic etiology often associated with other comorbidities. Using whole-exome sequencing (WES) of 63 patients with EoE and 60 unaffected family members and family-based trio analysis, we sought to uncover rare coding variants. WES analysis identified 5 rare, damaging variants in dehydrogenase E1 and transketolase domain-containing 1 (DHTKD1). Rare variant burden analysis revealed an overabundance of putative, potentially damaging DHTKD1 mutations in EoE (P = 0.01). Interestingly, we also identified 7 variants in the DHTKD1 homolog oxoglutarate dehydrogenase-like (OGDHL). Using shRNA-transduced esophageal epithelial cells and/or patient fibroblasts, we further showed that disruption of normal DHTKD1 or OGDHL expression blunts mitochondrial function. Finally, we demonstrated that the loss of DHTKD1 expression increased ROS production and induced the expression of viperin, a gene previously shown to be involved in production of Th2 cytokines in T cells. Viperin had increased expression in esophageal biopsies of EoE patients compared with control individuals and was upregulated by IL-13 in esophageal epithelial cells. These data identify a series of rare genetic variants implicating DHTKD1 and OGDHL in the genetic etiology of EoE and underscore a potential pathogenic role for mitochondrial dysfunction in EoE.

Transcription factors operate across disease loci, with EBNA2 implicated in autoimmunity.

Explaining the genetics of many diseases is challenging because most associations localize to incompletely characterized regulatory regions. Using new computational methods, we show that transcription factors (TFs) occupy multiple loci associated with individual complex genetic disorders. Application to 213 phenotypes and 1,544 TF binding datasets identified 2,264 relationships between hundreds of TFs and 94 phenotypes, including androgen receptor in prostate cancer and GATA3 in breast cancer. Strikingly, nearly half of systemic lupus erythematosus risk loci are occupied by the Epstein-Barr virus EBNA2 protein and many coclustering human TFs, showing gene-environment interaction. Similar EBNA2-anchored associations exist in multiple sclerosis, rheumatoid arthritis, inflammatory bowel disease, type 1 diabetes, juvenile idiopathic arthritis and celiac disease. Instances of allele-dependent DNA binding and downstream effects on gene expression at plausibly causal variants support genetic mechanisms dependent on EBNA2. Our results nominate mechanisms that operate across risk loci within disease phenotypes, suggesting new models for disease origins.

T-bet+CD11c+ B cells are critical for antichromatin immunoglobulin G production in the development of lupus.

BACKGROUND: A hallmark of systemic lupus erythematosus is high titers of circulating autoantibodies. Recently, a novel CD11c+ B-cell subset has been identified that is critical for the development of autoimmunity. However, the role of CD11c+ B cells in the development of lupus is unclear. Chronic graft-versus-host disease (cGVHD) is a lupus-like syndrome with high autoantibody production. The purpose of this study was to explore the role of CD11c+ B cells in the pathogenesis of lupus in cGVHD mice. METHODS: cGVHD was induced by an intraperitoneal injection of 5 × 107 Bm12 splenocytes into B6 mice. Flow cytometry was used to analyze mice splenocytes and human samples. Magnetic beads were used to isolate mice B cells. Gene expression was determined by real-time quantitative polymerase chain reaction (RT-qPCR). Enzyme-linked immunosorbent assay (ELISA) was used to detect antibodies in serum and supernatants. RESULTS: The percentage and absolute number of CD11c+ B cells was increased in cGVHD-induced lupus, with elevated levels of antichromatin immunoglobulin (Ig)G and IgG2a in sera. CD11c+ plasma cells from cGVHD mice produced large amounts of antichromatin IgG2a upon stimulation. Depletion of CD11c+ B cells reduced antichromatin IgG and IgG2a production. T-bet was upregulated in CD11c+ B cells. Knockout of T-bet in B cells alleviated cGVHD-induced lupus. Importantly, the percentage of T-bet+CD11c+ B cells increased in lupus patients and positively correlated with serum antichromatin levels. CONCLUSION: T-bet+CD11c+ B cells promoted high antichromatin IgG production in the lupus-like disease model cGVHD. In lupus patients, the percentage of T-bet+CD11c+ B cells was elevated and positively correlated with antichromatin antibodies. The findings provide potential therapeutic insight into lupus disease treatment.

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.

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.