Cell-Specific Transposable Element and Gene Expression Analysis Across Systemic Lupus Erythematosus Phenotypes.

OBJECTIVE: 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 SLE phenotypes, specifically production of type I IFNs and generation of autoantibodies. METHODS: We profiled cell-sorted RNA-sequencing data (CD4+ T cells, CD14+ monocytes, CD19+ B cells, and natural killer cells) from peripheral blood mononuclear cells of 120 patients with SLE and quantified TE expression identifying 27,135 TEs. We tested for differential TE expression across 10 SLE phenotypes, including autoantibody production and disease activity. RESULTS: We found 731 differentially expressed (DE) TEs across all SLE phenotypes that were mostly cell specific and phenotype specific. DE TEs were enriched for specific families and open reading frames of viral genes encoded in TE sequences. Increased expression of DE TEs was associated with genes involved in antiviral activity, such as LY6E, ISG15, and TRIM22, and pathways such as IFN signaling. CONCLUSION: 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.

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.

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.

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.

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.

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.

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.

Variants in the DDX6-CXCR5 autoimmune disease risk locus influence the regulatory network in immune cells and salivary gland.

Fine mapping and bioinformatic analysis of the DDX6-CXCR5 genetic risk association in Sjögren's Disease (SjD) and Systemic Lupus Erythematosus (SLE) identified five common SNPs with functional evidence in immune cell types: rs4938573, rs57494551, rs4938572, rs4936443, rs7117261. Functional interrogation of nuclear protein binding affinity, enhancer/promoter regulatory activity, and chromatin-chromatin interactions in immune, salivary gland epithelial, and kidney epithelial cells revealed cell type-specific allelic effects for all five SNPs that expanded regulation beyond effects on DDX6 and CXCR5 expression. Mapping the local chromatin regulatory network revealed several additional genes of interest, including lnc-PHLDB1-1. Collectively, functional characterization implicated the risk alleles of these SNPs as modulators of promoter and/or enhancer activities that regulate cell type-specific expression of DDX6, CXCR5, and lnc-PHLDB1-1, among others. Further, these findings emphasize the importance of exploring the functional significance of SNPs in the context of complex chromatin architecture in disease-relevant cell types and tissues.

Functional Testing of Bone Morphogenetic Protein (BMP) Pathway Variants Identified on Whole-Exome Sequencing in a Patient with Delayed-Onset Fibrodysplasia Ossificans Progressiva (FOP) Using ACVR1R206H -Specific Human Cellular and Zebrafish Models.

Bone morphogenetic protein (BMP) signaling is critical in skeletal development. Overactivation can trigger heterotopic ossification (HO) as in fibrodysplasia ossificans progressiva (FOP), a rare, progressive disease of massive HO formation. A small subset of FOP patients harboring the causative ACVR1R206H mutation show strikingly mild or delayed-onset HO, suggesting that genetic variants in the BMP pathway could act as disease modifiers. Whole-exome sequencing of one such patient identified BMPR1AR443C and ACVR2AV173I as candidate modifiers. Molecular modeling predicted significant structural perturbations. Neither variant decreased BMP signaling in ACVR1R206H HEK 293T cells at baseline or after stimulation with BMP4 or activin A (AA), ligands that activate ACVR1R206H signaling. Overexpression of BMPR1AR443C in a Tg(ACVR1-R206Ha) embryonic zebrafish model, in which overactive BMP signaling yields ventralized embryos, did not alter ventralization severity, while ACVR2AV173I exacerbated ventralization. Co-expression of both variants did not affect dorsoventral patterning. In contrast, BMPR1A knockdown in ACVR1R206H HEK cells decreased ligand-stimulated BMP signaling but did not affect dorsoventral patterning in Tg(ACVR1-R206Ha) zebrafish. ACVR2A knockdown decreased only AA-stimulated signaling in ACVR1R206H HEK cells and had no effect in Tg(ACVR1-R206Ha) zebrafish. Co-knockdown in ACVR1R206H HEK cells decreased basal and ligand-stimulated signaling, and co-knockdown/knockout (bmpr1aa/ab; acvr2aa/ab) decreased Tg(ACVR1-R206Ha) zebrafish ventralization phenotypes. Our functional studies showed that knockdown of wild-type BMPR1A and ACVR2A could attenuate ACVR1R206H signaling, particularly in response to AA, and that ACVR2AV173I unexpectedly increased ACVR1R206H -mediated signaling in zebrafish. These studies describe a useful strategy and platform for functionally interrogating potential genes and genetic variants that may impact the BMP signaling pathway. © 2022 The Authors. Journal of Bone and Mineral Research published by Wiley Periodicals LLC on behalf of American Society for Bone and Mineral Research (ASBMR).

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.

A Human MSH6 Germline Variant Associated With Systemic Lupus Erythematosus Induces Lupus-like Disease in Mice.

OBJECTIVE: To determine if single-nucleotide polymorphisms (SNPs) in DNA repair genes are enriched in individuals with systemic lupus erythematosus (SLE) and if they are sufficient to confer a disease phenotype in a mouse model. METHODS: Human exome chip data of 2499 patients with SLE and 1230 healthy controls were analyzed to determine if variants in 10 different mismatch repair genes (MSH4, EXO1, MSH2, MSH6, MLH1, MSH3, POLH, PMS2, ML3, and APEX2) were enriched in individuals with SLE. A mouse model of the MSH6 SNP, which was found to be enriched in individuals with SLE, was created using CRISPR/Cas9 gene targeting. Wildtype mice and mice heterozygous and homozygous for the MSH6 variant were then monitored for 2 years for the development of autoimmune phenotypes, including the presence of high levels of antinuclear antibodies (ANA). Additionally, somatic hypermutation frequencies and spectra of the intronic region downstream of the VH J558-rearranged JH4 immunoglobulin gene was characterized from Peyer's patches. RESULTS: Based on the human exome chip data, the MSH6 variant (rs63750897, p.Ser503Cys) is enriched among patients with SLE versus controls after we corrected for ancestry (odds ratio = 8.39, P = 0.0398). Mice homozygous for the MSH6 variant (Msh6S502C/S502C ) harbor significantly increased levels of ANA. Additionally, the Msh6S502C/S502C mice display a significant increase in the infiltration of CD68+ cells (a marker for monocytes and macrophages) into the lung alveolar space as well as apoptotic cells. Furthermore, characterization of somatic hypermutation in these mice reveals an increase in the DNA polymerase η mutational signature. CONCLUSION: An MSH6 mutation that is enriched in humans diagnosed with lupus was identified. Mice harboring this Msh6 mutation develop increased autoantibodies and an inflammatory lung disease. These results suggest that the human MSH6 variant is linked to the development of SLE.

Dynamics of Methylation of CpG Sites Associated With Systemic Lupus Erythematosus Subtypes in a Longitudinal Cohort.

OBJECTIVE: Findings from cross-sectional studies have revealed associations between DNA methylation and systemic lupus erythematosus (SLE) outcomes. This study was undertaken to investigate the dynamics of DNA methylation by examining participants from an SLE longitudinal cohort using samples collected at 2 time points. METHODS: A total of 101 participants from the California Lupus Epidemiology Study were included in our analysis. DNA was extracted from blood samples collected at the time of enrolment in the cohort and samples collected after 2 years and was analyzed using Illumina EPIC BeadChip kit. Paired t-tests were used to identify genome-wide changes which included 256 CpG sites previously found to be associated with SLE subtypes. Linear mixed models were developed to understand the relationship between DNA methylation and disease activity, medication use, and sample cell-type proportions, adjusted for age, sex, and genetic principal components. RESULTS: The majority of CpGs that were previously determined to be associated with SLE subtypes remained stable over 2 years (185 CpGs [72.3%]; t-test false discovery rate >0.05). Compared to background genome-wide methylation, there was an enrichment of SLE subtype-associated CpGs that changed over time (27.7% versus 0.34%). Changes in cell-type proportions were associated with changes at 67 CpGs (P < 2.70 × 10-5 ), and 15 CpGs had at least 1 significant association with immunosuppressant use. CONCLUSION: In this longitudinal SLE cohort, we identified a subset of SLE subtype-associated CpGs that remained stable over time and may be useful as biomarkers of disease subtypes. Another subset of SLE subtype-associated CpGs changed at a higher proportion compared to the genome-wide methylome. Additional studies are needed to understand the etiology and impact of these changes on methylation of SLE-associated CpGs.

Reactogenicity of the Messenger RNA SARS-CoV-2 Vaccines Associated With Immunogenicity in Patients With Autoimmune and Inflammatory Disease.

OBJECTIVE: Little is known regarding the reactogenicity and related SARS-CoV-2 vaccine response in patients with chronic inflammatory disease (CID). Our objective was to characterize the adverse event profile of CID patients following SARS-CoV-2 vaccination and understand the relationship between reactogenicity and immunogenicity of SARS-CoV-2 vaccines. METHODS: CID patients and healthy controls eligible to receive messenger RNA (mRNA) SARS-CoV-2 vaccines participated in 3 study visits (pre-vaccine, after dose 1, and after dose 2) in which blood and clinical data were collected. Assessment of adverse events were solicited within 7 days of receiving each dose. Serum anti-SARS-CoV-2 spike IgG ± antibody titers were quantified following vaccination. Statistical analysis was performed utilizing mixed models and tobit regressions, with adjustment for covariates. RESULTS: The present study included 441 participants (322 CID patients and 119 control subjects). Compared to controls, CID patients reported greater symptom severity after dose 1 (P = 0.0001), including more myalgia and fatigue (P < 0.05). For immunogenicity, a higher symptom severity after dose 1 and a higher number of symptoms after dose 2 was associated with higher antibody titers (P ≤ 0.05). Each increase of 1 symptom was associated with a 15.1% increase in antibody titer. Symptom association was strongest with site pain after dose 1 (105%; P = 0.03) and fatigue after dose 2 (113%; P = 0.004). CONCLUSION: Patients with CID have a distinct reactogenicity profile following SARS-CoV-2 vaccination compared to controls. Furthermore, there is an association between increased reactogenicity and increased vaccine response. This finding may speak to the more variable immunogenicity in CID patients and may be an important indicator of vaccine response to the novel SARS-CoV-2 vaccines.

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.

Correction: Hypomethylation mediates genetic association with the major histocompatibility complex genes in Sjögren's syndrome.

[This corrects the article DOI: 10.1371/journal.pone.0248429.].

Glucocorticoids and B Cell Depleting Agents Substantially Impair Immunogenicity of mRNA Vaccines to SARS-CoV-2.

BACKGROUND: Individuals with chronic inflammatory diseases (CID) are frequently treated with immunosuppressive medications that can increase their risk of severe COVID-19. While novel mRNA-based SARS-CoV-2 vaccination platforms provide robust protection in immunocompetent individuals, the immunogenicity in CID patients on immunosuppression is not well established. Therefore, determining the effectiveness of SARS-CoV-2 vaccines in the setting of immunosuppression is essential to risk-stratify CID patients with impaired protection and provide clinical guidance regarding medication management. METHODS: We conducted a prospective assessment of mRNA-based vaccine immunogenicity in 133 adults with CIDs and 53 immunocompetent controls. Blood from participants over 18 years of age was collected before initial immunization and 1-2 weeks after the second immunization. Serum anti-SARS-CoV-2 spike (S) IgG + binding, neutralizing antibody titers, and circulating S-specific plasmablasts were quantified to assess the magnitude and quality of the humoral response following vaccination. RESULTS: Compared to immunocompetent controls, a three-fold reduction in anti-S IgG titers (P=0.009) and SARS-CoV-2 neutralization (p<0.0001) were observed in CID patients. B cell depletion and glucocorticoids exerted the strongest effect with a 36- and 10-fold reduction in humoral responses, respectively (p<0.0001). Janus kinase inhibitors and antimetabolites, including methotrexate, also blunted antibody titers in multivariate regression analysis (P<0.0001, P=0.0023, respectively). Other targeted therapies, such as TNF inhibitors, IL-12/23 inhibitors, and integrin inhibitors, had only modest impacts on antibody formation and neutralization. CONCLUSIONS: CID patients treated with immunosuppressive therapies exhibit impaired SARS-CoV-2 vaccine-induced immunity, with glucocorticoids and B cell depletion therapy more severely impeding optimal responses.

Hypomethylation mediates genetic association with the major histocompatibility complex genes in Sjögren's syndrome.

Differential methylation of immune genes has been a consistent theme observed in Sjögren's syndrome (SS) in CD4+ T cells, CD19+ B cells, whole blood, and labial salivary glands (LSGs). Multiple studies have found associations supporting genetic control of DNA methylation in SS, which in the absence of reverse causation, has positive implications for the potential of epigenetic therapy. However, a formal study of the causal relationship between genetic variation, DNA methylation, and disease status is lacking. We performed a causal mediation analysis of DNA methylation as a mediator of nearby genetic association with SS using LSGs and genotype data collected from 131 female members of the Sjögren's International Collaborative Clinical Alliance registry, comprising of 64 SS cases and 67 non-cases. Bumphunter was used to first identify differentially-methylated regions (DMRs), then the causal inference test (CIT) was applied to identify DMRs mediating the association of nearby methylation quantitative trait loci (MeQTL) with SS. Bumphunter discovered 215 DMRs, with the majority located in the major histocompatibility complex (MHC) on chromosome 6p21.3. Consistent with previous findings, regions hypomethylated in SS cases were enriched for gene sets associated with immune processes. Using the CIT, we observed a total of 19 DMR-MeQTL pairs that exhibited strong evidence for a causal mediation relationship. Close to half of these DMRs reside in the MHC and their corresponding meQTLs are in the region spanning the HLA-DQA1, HLA-DQB1, and HLA-DQA2 loci. The risk of SS conferred by these corresponding MeQTLs in the MHC was further substantiated by previous genome-wide association study results, with modest evidence for independent effects. By validating the presence of causal mediation, our findings suggest both genetic and epigenetic factors contribute to disease susceptibility, and inform the development of targeted epigenetic modification as a therapeutic approach for SS.

Transcriptomic analysis of immune cells in a multi-ethnic cohort of systemic lupus erythematosus patients identifies ethnicity- and disease-specific expression signatures.

Systemic lupus erythematosus (SLE) is an autoimmune disease in which outcomes vary among different racial groups. We leverage cell-sorted RNA-seq data (CD14+ monocytes, B cells, CD4+ T cells, and NK cells) from 120 SLE patients (63 Asian and 57 White individuals) and apply a four-tier approach including unsupervised clustering, differential expression analyses, gene co-expression analyses, and machine learning to identify SLE subgroups within this multiethnic cohort. K-means clustering on each cell-type resulted in three clusters for CD4 and CD14, and two for B and NK cells. To understand the identified clusters, correlation analysis revealed significant positive associations between the clusters and clinical parameters including disease activity as well as ethnicity. We then explored differentially expressed genes between Asian and White groups for each cell-type. The shared differentially expressed genes across cells were involved in SLE or other autoimmune-related pathways. Co-expression analysis identified similarly regulated genes across samples and grouped these genes into modules. Finally, random forest classification of disease activity in the White and Asian cohorts showed the best classification in CD4+ T cells in White individuals. The results from these analyses will help stratify patients based on their gene expression signatures to enable SLE precision medicine.