Complement genes contribute sex-biased vulnerability in diverse disorders.

Many common illnesses, for reasons that have not been identified, differentially affect men and women. For instance, the autoimmune diseases systemic lupus erythematosus (SLE) and Sjögren's syndrome affect nine times more women than men1, whereas schizophrenia affects men with greater frequency and severity relative to women2. All three illnesses have their strongest common genetic associations in the major histocompatibility complex (MHC) locus, an association that in SLE and Sjögren's syndrome has long been thought to arise from alleles of the human leukocyte antigen (HLA) genes at that locus3-6. Here we show that variation of the complement component 4 (C4) genes C4A and C4B, which are also at the MHC locus and have been linked to increased risk for schizophrenia7, generates 7-fold variation in risk for SLE and 16-fold variation in risk for Sjögren's syndrome among individuals with common C4 genotypes, with C4A protecting more strongly than C4B in both illnesses. The same alleles that increase risk for schizophrenia greatly reduce risk for SLE and Sjögren's syndrome. In all three illnesses, C4 alleles act more strongly in men than in women: common combinations of C4A and C4B generated 14-fold variation in risk for SLE, 31-fold variation in risk for Sjögren's syndrome, and 1.7-fold variation in schizophrenia risk among men (versus 6-fold, 15-fold and 1.26-fold variation in risk among women, respectively). At a protein level, both C4 and its effector C3 were present at higher levels in cerebrospinal fluid and plasma8,9 in men than in women among adults aged between 20 and 50 years, corresponding to the ages of differential disease vulnerability. Sex differences in complement protein levels may help to explain the more potent effects of C4 alleles in men, women's greater risk of SLE and Sjögren's syndrome and men's greater vulnerability to schizophrenia. These results implicate the complement system as a source of sexual dimorphism in vulnerability to diverse illnesses.

Expression of A20 by dendritic cells preserves immune homeostasis and prevents colitis and spondyloarthritis.

Dendritic cells (DCs), which are known to support immune activation during infection, may also regulate immune homeostasis in resting animals. Here we show that mice lacking the ubiquitin-editing molecule A20 specifically in DCs spontaneously showed DC activation and population expansion of activated T cells. Analysis of DC-specific epistasis in compound mice lacking both A20 and the signaling adaptor MyD88 specifically in DCs showed that A20 restricted both MyD88-independent signals, which drive activation of DCs and T cells, and MyD88-dependent signals, which drive population expansion of T cells. In addition, mice lacking A20 specifically in DCs spontaneously developed lymphocyte-dependent colitis, seronegative ankylosing arthritis and enthesitis, conditions stereotypical of human inflammatory bowel disease (IBD). Our findings indicate that DCs need A20 to preserve immune quiescence and suggest that A20-dependent DC functions may underlie IBD and IBD-associated arthritides.

PICS2: next-generation fine mapping via probabilistic identification of causal SNPs.

SUMMARY: The Probabilistic Identification of Causal SNPs (PICS) algorithm and web application was developed as a fine-mapping tool to determine the likelihood that each single nucleotide polymorphism (SNP) in LD with a reported index SNP is a true causal polymorphism. PICS is notable for its ability to identify candidate causal SNPs within a locus using only the index SNP, which are widely available from published GWAS, whereas other methods require full summary statistics or full genotype data. However, the original PICS web application operates on a single SNP at a time, with slow performance, severely limiting its usability. We have developed a next-generation PICS tool, PICS2, which enables performance of PICS analyses of large batches of index SNPs with much faster performance. Additional updates and extensions include use of LD reference data generated from 1000 Genomes phase 3; annotation of variant consequences; annotation of GTEx eQTL genes and downloadable PICS SNPs from GTEx eQTLs; the option of generating PICS probabilities from experimental summary statistics; and generation of PICS SNPs from all SNPs of the GWAS catalog, automatically updated weekly. These free and easy-to-use resources will enable efficient determination of candidate loci for biological studies to investigate the true causal variants underlying disease processes. AVAILABILITY AND IMPLEMENTATION: PICS2 is available at https://pics2.ucsf.edu. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Effect of Immunosuppression on the Immunogenicity of mRNA Vaccines to SARS-CoV-2 : A Prospective Cohort Study.

BACKGROUND: Patients with chronic inflammatory disease (CID) treated with immunosuppressive medications have increased risk for severe COVID-19. Although mRNA-based SARS-CoV-2 vaccination provides protection in immunocompetent persons, immunogenicity in immunosuppressed patients with CID is unclear. OBJECTIVE: To determine the immunogenicity of mRNA-based SARS-CoV-2 vaccines in patients with CID. DESIGN: Prospective observational cohort study. SETTING: Two U.S. CID referral centers. PARTICIPANTS: Volunteer sample of adults with confirmed CID eligible for early COVID-19 vaccination, including hospital employees of any age and patients older than 65 years. Immunocompetent participants were recruited separately from hospital employees. All participants received 2 doses of mRNA vaccine against SARS-CoV-2 between 10 December 2020 and 20 March 2021. Participants were assessed within 2 weeks before vaccination and 20 days after final vaccination. MEASUREMENTS: Anti-SARS-CoV-2 spike (S) IgG+ binding in all participants, and neutralizing antibody titers and circulating S-specific plasmablasts in a subset to assess humoral response after vaccination. RESULTS: Most of the 133 participants with CID (88.7%) and all 53 immunocompetent participants developed antibodies in response to mRNA-based SARS-CoV-2 vaccination, although some with CID developed numerically lower titers of anti-S IgG. Anti-S IgG antibody titers after vaccination were lower in participants with CID receiving glucocorticoids (n = 17) than in those not receiving them; the geometric mean of anti-S IgG antibodies was 357 (95% CI, 96 to 1324) for participants receiving prednisone versus 2190 (CI, 1598 to 3002) for those not receiving it. Anti-S IgG antibody titers were also lower in those receiving B-cell depletion therapy (BCDT) (n = 10). Measures of immunogenicity differed numerically between those who were and those who were not receiving antimetabolites (n = 48), tumor necrosis factor inhibitors (n = 39), and Janus kinase inhibitors (n = 11); however, 95% CIs were wide and overlapped. Neutralization titers seemed generally consistent with anti-S IgG results. Results were not adjusted for differences in baseline clinical factors, including other immunosuppressant therapies. LIMITATIONS: Small sample that lacked demographic diversity, and residual confounding. CONCLUSION: Compared with nonusers, patients with CID treated with glucocorticoids and BCDT seem to have lower SARS-CoV-2 vaccine-induced antibody responses. These preliminary findings require confirmation in a larger study. PRIMARY FUNDING SOURCE: The Leona M. and Harry B. Helmsley Charitable Trust, Marcus Program in Precision Medicine Innovation, National Center for Advancing Translational Sciences, and National Institute of Arthritis and Musculoskeletal and Skin Diseases.

Role of the IL-12/IL-35 balance in patients with Sjögren syndrome.

BACKGROUND: An interferon signature is involved in the pathogenesis of primary Sjögren syndrome (pSS), but whether the signature is type 1 or type 2 remains controversial. Mouse models and genetic studies suggest the involvement of TH1 and type 2 interferon pathways. Likewise, polymorphisms of the IL-12A gene (IL12A), which encodes for IL-12p35, have been associated with pSS. The IL-12p35 subunit is shared by 2 heterodimers: IL-12 and IL-35. OBJECTIVE: We sought to confirm genetic association of the IL12A polymorphism and pSS and elucidate involvement of the IL-12/IL-35 balance in patients with pSS by using functional studies. METHODS: The genetic study involved 673 patients with pSS from 2 French pSS cohorts and 585 healthy French control subjects. Functional studies were performed on sorted monocytes, irrespective of whether they were stimulated. IL12A mRNA expression and IL-12 and IL-35 protein levels were assessed by using quantitative RT-PCR and ELISA and a multiplex kit for IL-35 and IL-12, respectively. RESULTS: We confirmed association of the IL12A rs485497 polymorphism and pSS and found an increased serum protein level of IL-12p70 in patients with pSS carrying the risk allele (P = .016). Serum levels of IL-12p70 were greater in patients than control subjects (P = .0001), especially in patients with more active disease (P = .05); conversely, IL-35 levels were decreased in patients (P = .0001), especially in patients with more active disease (P = .05). In blood cellular subsets both IL12p35 and EBV-induced gene protein 3 (EBI3) mRNAs were detected only in B cells, with a trend toward a lower level among patients with pSS. CONCLUSION: Our findings emphasize involvement of the IL-12/IL-35 balance in the pathogenesis of pSS. Serum IL-35 levels were associated with low disease activity, in contrast with serum IL-12p70 levels, which were associated with more active disease.

Novel gene variants associated with cardiovascular disease in systemic lupus erythematosus and rheumatoid arthritis.

OBJECTIVES: Patients with systemic lupus erythematosus (SLE) and rheumatoid arthritis (RA) have increased risk of cardiovascular disease (CVD). We investigated whether single nucleotide polymorphisms (SNPs) at autoimmunity risk loci were associated with CVD in SLE and RA. METHODS: Patients with SLE (n=1045) were genotyped using the 200K Immunochip SNP array (Illumina). The allele frequency was compared between patients with and without different manifestations of CVD. Results were replicated in a second SLE cohort (n=1043) and in an RA cohort (n=824). We analysed publicly available genetic data from general population, performed electrophoretic mobility shift assays and measured cytokine levels and occurrence of antiphospholipid antibodies (aPLs). RESULTS: We identified two new putative risk loci associated with increased risk for CVD in two SLE populations, which remained after adjustment for traditional CVD risk factors. An IL19 risk allele, rs17581834(T) was associated with stroke/myocardial infarction (MI) in SLE (OR 2.3 (1.5 to 3.4), P=8.5×10-5) and RA (OR 2.8 (1.4 to 5.6), P=3.8×10-3), meta-analysis (OR 2.5 (2.0 to 2.9), P=3.5×10-7), but not in population controls. The IL19 risk allele affected protein binding, and SLE patients with the risk allele had increased levels of plasma-IL10 (P=0.004) and aPL (P=0.01). An SRP54-AS1 risk allele, rs799454(G) was associated with stroke/transient ischaemic attack in SLE (OR 1.7 (1.3 to 2.2), P=2.5×10-5) but not in RA. The SRP54-AS1 risk allele is an expression quantitative trait locus for four genes. CONCLUSIONS: The IL19 risk allele was associated with stroke/MI in SLE and RA, but not in the general population, indicating that shared immune pathways may be involved in the CVD pathogenesis in inflammatory rheumatic diseases.

Multiple polymorphisms in the TNFAIP3 region are independently associated with systemic lupus erythematosus.

The TNFAIP3 (tumor necrosis factor alpha-induced protein 3) gene encodes a ubiquitin editing enzyme, A20, that restricts NF-kappaB-dependent signaling and prevents inflammation. We show that three independent SNPs in the TNFAIP3 region (rs13192841, rs2230926 and rs6922466) are associated with systemic lupus erythematosus (SLE) among individuals of European ancestry. These findings provide critical links between A20 and the etiology of SLE.

Unraveling multiple MHC gene associations with systemic lupus erythematosus: model choice indicates a role for HLA alleles and non-HLA genes in Europeans.

We have performed a meta-analysis of the major-histocompatibility-complex (MHC) region in systemic lupus erythematosus (SLE) to determine the association with both SNPs and classical human-leukocyte-antigen (HLA) alleles. More specifically, we combined results from six studies and well-known out-of-study control data sets, providing us with 3,701 independent SLE cases and 12,110 independent controls of European ancestry. This study used genotypes for 7,199 SNPs within the MHC region and for classical HLA alleles (typed and imputed). Our results from conditional analysis and model choice with the use of the Bayesian information criterion show that the best model for SLE association includes both classical loci (HLA-DRB1(∗)03:01, HLA-DRB1(∗)08:01, and HLA-DQA1(∗)01:02) and two SNPs, rs8192591 (in class III and upstream of NOTCH4) and rs2246618 (MICB in class I). Our approach was to perform a stepwise search from multiple baseline models deduced from a priori evidence on HLA-DRB1 lupus-associated alleles, a stepwise regression on SNPs alone, and a stepwise regression on HLA alleles. With this approach, we were able to identify a model that was an overwhelmingly better fit to the data than one identified by simple stepwise regression either on SNPs alone (Bayes factor [BF] > 50) or on classical HLA alleles alone (BF > 1,000).

Germline and somatic genetic variations of TNFAIP3 in lymphoma complicating primary Sjogren's syndrome.

Several autoimmune diseases, including primary Sjögren's syndrome (pSS), are associated with an increased risk for lymphoma. Polymorphisms of TNFAIP3, which encodes the A20 protein that plays a key role in controlling nuclear factor κB activation, have been associated with several autoimmune diseases. Somatic mutations of TNFAIP3 have been observed in the mucosa-associated lymphoid tissue lymphoma subtype frequently associated with pSS. We studied germline and somatic abnormalities of TNFAIP3 in 574 patients with pSS, including 25 with lymphoma. Nineteen additional patients with pSS and lymphoma were available for exome sequence analysis. Functional abnormalities of A20 were assessed by gene reporter assays. The rs2230926 exonic variant was associated with an increased risk for pSS complicated by lymphoma (odds ratio, 3.36 [95% confidence interval, 1.34-8.42], and odds ratio, 3.26 [95% confidence interval, 1.31-8.12], vs controls and pSS patients without lymphoma, respectively; P = .011). Twelve (60%) of the 20 patients with paired germline and lymphoma TNFAIP3 sequence data had functional abnormalities of A20: 6 in germline DNA, 5 in lymphoma DNA, and 1 in both. The frequency was even higher (77%) among pSS patients with mucosa-associated lymphoid tissue lymphoma. Some of these variants showed impaired control of nuclear factor κB activation. These results support a key role for germline and somatic variations of A20 in the transformation between autoimmunity and lymphoma.

Brief Report: Single-nucleotide polymorphisms in VKORC1 are risk factors for systemic lupus erythematosus in Asians.

OBJECTIVE: The increased risk of thrombosis in systemic lupus erythematosus (SLE) may be partially explained by interrelated genetic pathways for thrombosis and SLE. The present study was undertaken to investigate whether 33 established and novel single-nucleotide polymorphisms (SNPs) in 20 genes involved in hemostasis pathways that have been associated with deep venous thrombosis (DVT) in the general population are risk factors for SLE among Asian subjects. METHODS: Patients in the discovery cohort were enrolled in 1 of 2 North American SLE cohorts. Patients in the replication cohort were enrolled in 1 of 4 Asian or 2 North American cohorts. We first genotyped 263 Asian patients with SLE and 357 healthy Asian control subjects for 33 SNPs in the discovery phase, and then genotyped 5 SNPs in up to an additional 1,496 patients and 993 controls in the replication phase. Patients were compared to controls for bivariate association with minor alleles. Principal components analysis was used to control for intra-Asian ancestry in the replication cohort. RESULTS: Two genetic variants in the gene VKORC1 were highly significant in both the discovery and replication cohorts: rs9934438 (in the discovery cohort, odds ratio [OR] 2.45, P=2×10(-9); in the replication cohort, OR 1.54, P=4×10(-6)) and rs9923231 (in the discovery cohort, OR 2.40, P=6×10(-9); in the replication cohort, OR 1.53, P=5×10(-6)). These associations were significant in the replication cohort after adjustment for intra-Asian ancestry: for rs9934438, OR 1.34, P=0.0029; for rs9923231, OR 1.34, P=0.0032. CONCLUSION: Genetic variants in VKORC1, which are involved in vitamin K reduction and associated with DVT, correlate with SLE development in Asian subjects. These results suggest that there may be intersecting genetic pathways for the development of SLE and thrombosis.

Differential genetic associations for systemic lupus erythematosus based on anti-dsDNA autoantibody production.

Systemic lupus erythematosus (SLE) is a clinically heterogeneous, systemic autoimmune disease characterized by autoantibody formation. Previously published genome-wide association studies (GWAS) have investigated SLE as a single phenotype. Therefore, we conducted a GWAS to identify genetic factors associated with anti-dsDNA autoantibody production, a SLE-related autoantibody with diagnostic and clinical importance. Using two independent datasets, over 400,000 single nucleotide polymorphisms (SNPs) were studied in a total of 1,717 SLE cases and 4,813 healthy controls. Anti-dsDNA autoantibody positive (anti-dsDNA +, n = 811) and anti-dsDNA autoantibody negative (anti-dsDNA -, n = 906) SLE cases were compared to healthy controls and to each other to identify SNPs associated specifically with these SLE subtypes. SNPs in the previously identified SLE susceptibility loci STAT4, IRF5, ITGAM, and the major histocompatibility complex were strongly associated with anti-dsDNA + SLE. Far fewer and weaker associations were observed for anti-dsDNA - SLE. For example, rs7574865 in STAT4 had an OR for anti-dsDNA + SLE of 1.77 (95% CI 1.57-1.99, p = 2.0E-20) compared to an OR for anti-dsDNA - SLE of 1.26 (95% CI 1.12-1.41, p = 2.4E-04), with p(heterogeneity)<0.0005. SNPs in the SLE susceptibility loci BANK1, KIAA1542, and UBE2L3 showed evidence of association with anti-dsDNA + SLE and were not associated with anti-dsDNA - SLE. In conclusion, we identified differential genetic associations with SLE based on anti-dsDNA autoantibody production. Many previously identified SLE susceptibility loci may confer disease risk through their role in autoantibody production and be more accurately described as autoantibody propensity loci. Lack of strong SNP associations may suggest that other types of genetic variation or non-genetic factors such as environmental exposures have a greater impact on susceptibility to anti-dsDNA - SLE.

Risk alleles for systemic lupus erythematosus in a large case-control collection and associations with clinical subphenotypes.

Systemic lupus erythematosus (SLE) is a genetically complex disease with heterogeneous clinical manifestations. Recent studies have greatly expanded the number of established SLE risk alleles, but the distribution of multiple risk alleles in cases versus controls and their relationship to subphenotypes have not been studied. We studied 22 SLE susceptibility polymorphisms with previous genome-wide evidence of association (p < 5 x 10⁻¹²8) in 1919 SLE cases from 9 independent Caucasian SLE case series and 4813 independent controls. The mean number of risk alleles in cases was 15.1 (SD 3.1) while the mean in controls was 13.1 (SD 2.8), with trend p = 4 x 10⁻8. We defined a genetic risk score (GRS) for SLE as the number of risk alleles with each weighted by the SLE risk odds ratio (OR). The OR for high-low GRS tertiles, adjusted for intra-European ancestry, sex, and parent study, was 4.4 (95% CI 3.8-5.1). We studied associations of individual SNPs and the GRS with clinical manifestations for the cases: age at diagnosis, the 11 American College of Rheumatology classification criteria, and double-stranded DNA antibody (anti-dsDNA) production. Six subphenotypes were significantly associated with the GRS, most notably anti-dsDNA (OR(high-low) = 2.36, p = 9e-9), the immunologic criterion (OR(high-low) = 2.23, p = 3e-7), and age at diagnosis (OR(high-low) = 1.45, p = 0.0060). Finally, we developed a subphenotype-specific GRS (sub-GRS) for each phenotype with more power to detect cumulative genetic associations. The sub-GRS was more strongly associated than any single SNP effect for 5 subphenotypes (the above plus hematologic disorder and oral ulcers), while single loci are more significantly associated with renal disease (HLA-DRB1, OR = 1.37, 95% CI 1.14-1.64) and arthritis (ITGAM, OR = 0.72, 95% CI 0.59-0.88). We did not observe significant associations for other subphenotypes, for individual loci or the sub-GRS. Thus our analysis categorizes SLE subphenotypes into three groups: those having cumulative, single, and no known genetic association with respect to the currently established SLE risk loci.

Screening of Environmental Chemicals to Characterize Exposures in Participants With Systemic Lupus Erythematosus.

OBJECTIVE: There is a need to characterize exposures associated with the pathogenesis of systemic lupus erythematosus (SLE). In this pilot study, we explore a hypothesis-free approach that can measure thousands of exogenous chemicals in blood ("exposome") in patients with SLE and unaffected controls. METHODS: This cross-sectional study analyzed a cohort of patients with prevalent SLE (n = 285) and controls (n = 106). Plasma was analyzed by liquid chromatography-quadrupole time-of-flight mass spectrometry (LC-QTOF/MS). Mass spectrometry features present in at least 25% of all samples were selected for association analysis (n = 2,737). Features were matched to potential chemicals using available databases. Association analysis of abundances of features with SLE status was performed, adjusting for age and sex. We also explored features associated with SLE phenotypes, sociodemographic factors, and current medication use. RESULTS: We found 30 features significantly associated with SLE status (Bonferroni P < 0.05). Of these, seven matched chemical names based on databases. These seven features included phthalate metabolites, a formetanate metabolite, and eugenol. The abundance of acid pesticides differed between patients with SLE and controls (Bonferroni P < 0.05). Two unmatched features were associated with a history of lupus nephritis, and one with anti-double-stranded DNA antibody production (Bonferroni P < 0.05). Seventeen features varied by self-reported race and ethnicity, including a polyfluoroalkyl substance (analysis of variance P < 1.69 × 10-5). Eleven features correlated with antimalarials, 6 with mycophenolate mofetil, and 29 with prednisone use. CONCLUSION: This proof-of-concept study demonstrates that LC-QTOF/MS is a powerful tool that agnostically detects circulating exogenous compounds. These analyses can generate hypotheses of disease-related exposures for future prospective, longitudinal studies.

Local Ancestry at the Major Histocompatibility Complex Region is Not a Major Contributor to Disease Heterogeneity in a Multiethnic Lupus Cohort.

OBJECTIVE: Systemic lupus erythematosus (SLE) is an autoimmune disease resulting in debilitating clinical manifestations that vary in severity by race and ethnicity with a disproportionate burden in African American, Mestizo, and Asian populations compared with populations of European descent. Differences in global and local genetic ancestry may shed light on the underlying mechanisms contributing to these disparities, including increased prevalence of lupus nephritis, younger age of symptom onset, and presence of autoantibodies. METHODS: A total of 1,139 European, African American, and Mestizos patients with SLE were genotyped using the Affymetrix LAT1 World array. Global ancestry proportions were estimated using ADMIXTURE, and local ancestry was estimated using RFMIXv2.0. We investigated associations between lupus nephritis, age at onset, and autoantibody status with both global and local ancestry proportions within the Major Histocompatibility Complex region. RESULTS: Our results showed small effect sizes that did not meet the threshold for statistical significance for global or local ancestry proportions in either African American or Mestizo patients with SLE who presented with the clinical manifestations of interest compared with those who did not. CONCLUSION: These findings suggest that local genetic ancestry within the Major Histocompatibility Complex region is not a major contributor to these SLE manifestations among patients with SLE from admixed populations.

European genetic ancestry is associated with a decreased risk of lupus nephritis.

OBJECTIVE: African Americans, East Asians, and Hispanics with systemic lupus erythematosus (SLE) are more likely to develop renal disease than are SLE patients of European descent. This study was undertaken to investigate whether European genetic ancestry protects against the development of lupus nephritis, with the aim of exploring the genetic and socioeconomic factors that might explain this effect. METHODS: This was a cross-sectional study of SLE patients from a multiethnic case collection. Participants were genotyped for 126 single-nucleotide polymorphisms (SNPs) informative for ancestry. A subset of participants was also genotyped for 80 SNPs in 14 candidate genes for renal disease in SLE. Logistic regression was used to test the association between European ancestry and renal disease. Analyses were adjusted for continental ancestries, socioeconomic status (SES), and candidate genes. RESULTS: Participants (n = 1,906) had, on average, 62.4% European, 15.8% African, 11.5% East Asian, 6.5% Amerindian, and 3.8% South Asian ancestry. Among the participants, 656 (34%) had renal disease. A 10% increase in the proportion of European ancestry estimated in each participant was associated with a 15% reduction in the odds of having renal disease, after adjustment for disease duration and sex (odds ratio 0.85, 95% confidence interval 0.82-0.87; P = 1.9 × 10(-30) ). Adjustment for other genetic ancestries, measures of SES, or SNPs in the genes most associated with renal disease (IRF5 [rs4728142], BLK [rs2736340], STAT4 [rs3024912], and HLA-DRB1*0301 and DRB1*1501) did not substantively alter this relationship. CONCLUSION: European ancestry is protective against the development of renal disease in SLE, an effect that is independent of other genetic ancestries, candidate risk alleles, and socioeconomic factors.

Identification of Sjögren's syndrome patient subgroups by clustering of labial salivary gland DNA methylation profiles.

Heterogeneity in Sjögren's syndrome (SS), increasingly called Sjögren's disease, suggests the presence of disease subtypes, which poses a major challenge for the diagnosis, management, and treatment of this autoimmune disorder. Previous work distinguished patient subgroups based on clinical symptoms, but it is not clear to what extent symptoms reflect underlying pathobiology. The purpose of this study was to discover clinical meaningful subtypes of SS based on genome-wide DNA methylation data. We performed a cluster analysis of genome-wide DNA methylation data from labial salivary gland (LSG) tissue collected from 64 SS cases and 67 non-cases. Specifically, hierarchical clustering was performed on low dimensional embeddings of DNA methylation data extracted from a variational autoencoder to uncover unknown heterogeneity. Clustering revealed clinically severe and mild subgroups of SS. Differential methylation analysis revealed that hypomethylation at the MHC and hypermethylation at other genome regions characterize the epigenetic differences between these SS subgroups. Epigenetic profiling of LSGs in SS yields new insights into mechanisms underlying disease heterogeneity. The methylation patterns at differentially methylated CpGs are different in SS subgroups and support the role of epigenetic contributions to the heterogeneity in SS. Biomarker data derived from epigenetic profiling could be explored in future iterations of the classification criteria for defining SS subgroups.

Cell-Specific Transposable Element Gene Expression Analysis Identifies Associations with Systemic Lupus Erythematosus Phenotypes.

There is an established yet unexplained link between interferon (IFN) and systemic lupus erythematosus (SLE). The expression of sequences derived from transposable elements (TEs) may contribute to production of type I IFNs and generation of autoantibodies. We profiled cell-sorted RNA-seq data (CD4+ T cells, CD14+ monocytes, CD19+ B cells, and NK cells) from PBMCs of 120 SLE patients and quantified TE expression identifying 27,135 TEs. We tested for differential TE expression across 10 SLE phenotypes including autoantibody production and disease activity and discovered 731 differentially expressed (DE) TEs whose effects were mostly cell-specific and phenotype-specific. DE TEs were enriched for specific families and viral genes encoded in TE sequences. Increased expression of DE TEs was associated with genes involved in antiviral activity such as LY6E, ISG15, TRIM22 and pathways such as interferon signaling. These findings suggest that expression of TEs contributes to activation of SLE-related mechanisms in a cell-specific manner, which can impact disease diagnostics and therapeutics.

GCLiPP: global crosslinking and protein purification method for constructing high-resolution occupancy maps for RNA binding proteins.

GCLiPP is a global RNA interactome capture method that detects RNA-binding protein (RBP) occupancy transcriptome-wide. GCLiPP maps RBP-occupied sites at a higher resolution than phase separation-based techniques. GCLiPP sequence tags correspond with known RBP binding sites and are enriched for sites detected by RBP-specific crosslinking immunoprecipitation (CLIP) for abundant cytosolic RBPs. Comparison of human Jurkat T cells and mouse primary T cells uncovers shared peaks of GCLiPP signal across homologous regions of human and mouse 3' UTRs, including a conserved mRNA-destabilizing cis-regulatory element. GCLiPP signal overlapping with immune-related SNPs uncovers stabilizing cis-regulatory regions in CD5, STAT6, and IKZF1.

Sleep Disturbance and SARS-CoV-2 Vaccinations in Patients With Chronic Inflammatory Disease.

OBJECTIVE: Immunocompromised patients with chronic inflammatory disease (CID) may have experienced additional psychosocial burden during the COVID-19 pandemic due to their immunocompromised status. This study was undertaken to determine if vaccination would result in improved patient-reported outcomes longitudinally among individuals with CID undergoing SARS-CoV-2 vaccination regardless of baseline anxiety. METHODS: Data are from a cohort of individuals with CID from 2 sites who underwent SARS-CoV-2 vaccination. Participants completed 3 study visits before and after 2 messenger RNA vaccine doses in the initial vaccination series when clinical data were collected. Patient-reported outcomes were measured using the Patient-Reported Outcomes Measurement Information System 29-item Health Profile and expressed as T scores, with 2 groups stratified by high and low baseline anxiety. Mixed-effects models were used to examine longitudinal changes, adjusting for age, sex, and study site. RESULTS: A total of 72% of the cohort was female with a mean ± SD age of 48.1 ± 15.5 years. Overall, sleep disturbance improved following both doses of SARS-CoV-2 vaccinations, and anxiety decreased after the second dose. Physical function scores worsened but did not meet the minimally important difference threshold. When stratifying by baseline anxiety, improvement in anxiety, fatigue, and social participation were greater in the high anxiety group. Physical function worsened slightly in both groups, and sleep disturbance improved significantly in the high anxiety group. CONCLUSION: Sleep disturbance decreased in a significant and meaningful way in patients with CID upon vaccination. In patients with higher baseline anxiety, social participation increased, and anxiety, fatigue, and sleep disturbance decreased. Overall, results suggest that SARS-CoV-2 vaccination may improve mental health and well-being, particularly among those with greater anxiety.

High-density SNP screening of the major histocompatibility complex in systemic lupus erythematosus demonstrates strong evidence for independent susceptibility regions.

A substantial genetic contribution to systemic lupus erythematosus (SLE) risk is conferred by major histocompatibility complex (MHC) gene(s) on chromosome 6p21. Previous studies in SLE have lacked statistical power and genetic resolution to fully define MHC influences. We characterized 1,610 Caucasian SLE cases and 1,470 parents for 1,974 MHC SNPs, the highly polymorphic HLA-DRB1 locus, and a panel of ancestry informative markers. Single-marker analyses revealed strong signals for SNPs within several MHC regions, as well as with HLA-DRB1 (global p = 9.99 x 10(-16)). The most strongly associated DRB1 alleles were: *0301 (odds ratio, OR = 2.21, p = 2.53 x 10(-12)), *1401 (OR = 0.50, p = 0.0002), and *1501 (OR = 1.39, p = 0.0032). The MHC region SNP demonstrating the strongest evidence of association with SLE was rs3117103, with OR = 2.44 and p = 2.80 x 10(-13). Conditional haplotype and stepwise logistic regression analyses identified strong evidence for association between SLE and the extended class I, class I, class III, class II, and the extended class II MHC regions. Sequential removal of SLE-associated DRB1 haplotypes revealed independent effects due to variation within OR2H2 (extended class I, rs362521, p = 0.006), CREBL1 (class III, rs8283, p = 0.01), and DQB2 (class II, rs7769979, p = 0.003, and rs10947345, p = 0.0004). Further, conditional haplotype analyses demonstrated that variation within MICB (class I, rs3828903, p = 0.006) also contributes to SLE risk independent of HLA-DRB1*0301. Our results for the first time delineate with high resolution several MHC regions with independent contributions to SLE risk. We provide a list of candidate variants based on biologic and functional considerations that may be causally related to SLE risk and warrant further investigation.