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.

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.

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.

Medication Cost Concerns and Disparities in Patient-Reported Outcomes Among a Multiethnic Cohort of Patients With Systemic Lupus Erythematosus.

OBJECTIVE: Concerns about the affordability of medications are common in systemic lupus erythematosus (SLE), but the relationship between medication cost concerns and health outcomes is poorly understood. We assessed the association of self-reported medication cost concerns and patient-reported outcomes (PROs) in a multiethnic SLE cohort. METHODS: The California Lupus Epidemiology Study is a cohort of individuals with physician-confirmed SLE. Medication cost concerns were defined as having difficulties affording SLE medications, skipping doses, delaying refills, requesting lower-cost alternatives, purchasing medications outside the United States, or applying for patient assistance programs. Linear regression and mixed effects models assessed the cross-sectional and longitudinal association of medication cost concerns and PROs, respectively, adjusting for age, sex, race and ethnicity, income, principal insurance, immunomodulatory medications, and organ damage. RESULTS: Of 334 participants, medication cost concerns were reported by 91 (27%). Medication cost concerns were associated with worse Systemic Lupus Activity Questionnaire (SLAQ; beta coefficient [ß] 5.9, 95% CI 4.3-7.6; P < 0.001), 8-item Patient Health Questionnaire depression scale (PHQ-8; ß 2.7, 95% CI 1.4-4.0; P < 0.001), and Patient-Reported Outcomes Measurement Information System (PROMIS; ß for physical function -4.6, 95% CI -6.7 to -2.4; P < 0.001) scores after adjusting for covariates. Medication cost concerns were not associated with significant changes in PROs over 2-year follow-up. CONCLUSION: More than a quarter of participants reported at least 1 medication cost concern, which was associated with worse PROs. Our results reveal a potentially modifiable risk factor for poor outcomes rooted in the unaffordability of SLE care.

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.

Identification of Cell-Specific Differential DNA Methylation Associated With Methotrexate Treatment Response in Rheumatoid Arthritis.

OBJECTIVE: We undertook this study to estimate changes in cell-specific DNA methylation (DNAm) associated with methotrexate (MTX) response using whole blood samples collected from rheumatoid arthritis (RA) patients before and after initiation of MTX treatment. METHODS: Patients included in this study were from the Rheumatoid Arthritis Medication Study (n = 66) and the University of California San Francisco Rheumatoid Arthritis study (n = 11). All patients met the American College of Rheumatology RA classification criteria. Blood samples were collected at baseline and following treatment. Disease Activity Scores in 28 joints using the C-reactive protein level were collected at baseline and after 3-6 months of treatment with MTX. Methylation profiles were generated using the Illumina Infinium HumanMethylation450 and MethylationEPIC v1.0 BeadChip arrays using DNA from whole blood. MTX response was defined using the EULAR response criteria (responders showed good/moderate response; nonresponders showed no response). Differentially methylated positions were identified using the Limma software package and Tensor Composition Analysis, which is a method for identifying cell-specific differential DNAm at the CpG level from tissue-level ("bulk") data. Differentially methylated regions were identified using Comb-p software. RESULTS: We found evidence of differential global methylation between treatment response groups. Further, we found patterns of cell-specific differential global methylation associated with MTX response. After correction for multiple testing, 1 differentially methylated position was associated with differential DNAm between responders and nonresponders at baseline in CD4+ T cells, CD8+ T cells, and natural killer cells. Thirty-nine cell-specific differentially methylated regions associated with MTX treatment response were identified. There were no significant findings in analyses of whole blood samples. CONCLUSION: We identified cell-specific changes in DNAm that were associated with MTX treatment response in RA patients. Future studies of DNAm and MTX treatment response should include measurements of DNAm from sorted cells.

The association of trauma with self-reported flares and disease activity in systemic lupus erythematosus (SLE).

OBJECTIVES: Trauma has been linked to incident SLE, but its relationship with SLE disease activity is unknown. This analysis examines associations between trauma exposures and patient-reported SLE disease activity and flares. METHODS: Data were from the California Lupus Epidemiology Study (CLUES). Flares were self-reported as any flare and, of those, flares accompanied by medical care (hospitalization or physician contact). The Systemic Lupus Activity Questionnaire (SLAQ) assessed disease activity. The Brief Trauma Questionnaire (BTQ) assessed all historical trauma exposures. The Adverse Childhood Experiences (ACEs) questionnaire was available for a subset. Multivariable regression analyses (n = 252) examined whether trauma exposure was associated with flares or SLAQ controlling for age, sex, poverty, race/ethnicity, comorbidities, perceived stress, disease duration and self-reported disease damage. RESULTS: Excluding exposure to serious illness, 63.4% reported ≥1 trauma exposure. Any traumatic event, excluding illness, doubled the odds of a flare [OR 2.27 (95% CI 1.24, 4.17)] and was associated with significantly higher SLAQ scores [ß 2.31 (0.86, 3.76)]. Adjusted odds of any flare and flare with medical care were significantly elevated for those with both BTQ and ACE exposures [5.91 (2.21, 15.82) and 4.69 (1.56, 14.07), respectively]. SLAQ scores were also higher for those with both exposures [ß 5.22 (3.00, 7.44)]. CONCLUSION: In this cohort, those with a history of trauma reported more flares and greater disease activity. Identifying mechanisms of associations between trauma and disease activity and flares, as well as interventions to mitigate the effects of trauma exposures is critical, given the high rates of trauma exposures.

Mer-tyrosine kinase: a novel susceptibility gene for SLE related end-stage renal disease.

OBJECTIVE: Lupus nephritis (LN) is a common and severe manifestation of SLE. The genetic risk for nephritis and progression to end-stage renal disease (ESRD) in patients with LN remains unclear. Herein, we aimed to identify novel genetic associations with LN, focusing on subphenotypes and ESRD. METHODS: We analysed genomic data on 958 patients with SLE (discovery cohort: LN=338) with targeted sequencing data from 1832 immunological pathway genes. We used an independent multiethnic cohort comprising 1226 patients with SLE (LN=603) as a replication dataset. Detailed functional annotation and functional epigenomic enrichment analyses were applied to predict functional effects of the candidate variants. RESULTS: A genetic variant (rs56097910) within the MERTK gene was associated with ESRD in both cohorts, meta-analysis OR=5.4 (2.8 to 10.6); p=1.0×10-6. We observed decreased methylation levels in peripheral blood cells from SLE patients with ESRD, compared with patients without renal SLE (p=2.7×10-4), at one CpG site (cg16333401) in close vicinity to the transcription start site of MERTK and located in a DNAse hypersensitivity region in T and B cells. Rs56097910 is linked to altered MERTK expression in kidney tissue in public eQTL databases. Two loci were replicated for association with proliferative LN: PRDM1 (rs6924535, pmeta=1.6×10-5, OR=0.58) and APOA1BP (NAXE) (rs942960, pmeta=1.2×10-5, OR=2.64). CONCLUSION: We identified a novel genetic risk locus, MERTK, associated with SLE-ESRD using the data from two large SLE cohorts. Through DNA methylation analysis and functional annotation, we showed that the risk could be mediated through regulation of gene expression. Our results suggest that variants in the MERTK gene are important for the risk of developing SLE-ESRD and suggest a role for PRDM1 and APOA1BP in proliferative LN.

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.

Single-cell RNA-seq reveals cell type-specific molecular and genetic associations to lupus.

Systemic lupus erythematosus (SLE) is a heterogeneous autoimmune disease. Knowledge of circulating immune cell types and states associated with SLE remains incomplete. We profiled more than 1.2 million peripheral blood mononuclear cells (162 cases, 99 controls) with multiplexed single-cell RNA sequencing (mux-seq). Cases exhibited elevated expression of type 1 interferon-stimulated genes (ISGs) in monocytes, reduction of naïve CD4+ T cells that correlated with monocyte ISG expression, and expansion of repertoire-restricted cytotoxic GZMH+ CD8+ T cells. Cell type-specific expression features predicted case-control status and stratified patients into two molecular subtypes. We integrated dense genotyping data to map cell type-specific cis-expression quantitative trait loci and to link SLE-associated variants to cell type-specific expression. These results demonstrate mux-seq as a systematic approach to characterize cellular composition, identify transcriptional signatures, and annotate genetic variants associated with SLE.

High Disease Severity Among Asian Patients in a US Multiethnic Cohort of Individuals With Systemic Lupus Erythematosus.

OBJECTIVE: Knowledge about systemic lupus erythematosus (SLE) outcomes among US Asian patients is lacking. The present study was undertaken to examine SLE disease activity, severity, and damage among Asian patients of primarily Chinese and Filipino descent in a multiethnic cohort. METHODS: California Lupus Epidemiology Study (n = 328) data were analyzed. Data were collected in English, Cantonese, Mandarin, or Spanish using validated instruments for disease activity (Systemic Lupus Erythematosus Disease Activity Index), disease severity (Lupus Severity Index [LSI]), and disease damage (Systemic Lupus International Collaborating Clinics/American College of Rheumatology Damage Index). We assessed differences in SLE outcomes among racial/ethnic groups using multivariable linear regression including interaction terms for age at diagnosis and race/ethnicity. RESULTS: Asian was the largest racial/ethnic group (38% [Chinese = 22%; Filipino = 9%; Other = 7%]). Average age at diagnosis was younger among Asian patients (27.9 years), particularly Filipino patients (22.2 years), compared with White (29.4 years) and Black patients (34.0 years). After adjustment, disease activity and damage were not significantly different across groups. Disease severity among Asian patients was significantly higher than among White patients (LSI score 7.1 versus 6.5; P < 0.05) but similar among Black and Hispanic patients. Early age at diagnosis was associated with greater organ damage among Asian, Black, and Hispanic patients, but not White patients. CONCLUSION: SLE was more severe among US Asian patients compared to White patients. Filipinos were affected at strikingly young ages. Asian patients and non-White groups with younger age at diagnosis had greater organ damage than White patients. Such racial/ethnic distinctions suggest the need for heightened clinical awareness to improve health outcomes among Asian patients with SLE. Further study of SLE outcomes across a range of US Asian subgroups is important.

Associations Between Smoking and Primary Sjögren Syndrome Classification Using the Sjögren's International Collaborative Clinical Alliance Cohort.

OBJECTIVE: The objective of this study was to examine the association of smoking with Primary Sjögren syndrome (pSS) classification and pSS diagnostic test results. We hypothesized that past and current smokers would have lower odds of being classified as having Sjögren syndrome (SS) and lower odds of having abnormal individual SS diagnostic test results compared with nonsmokers. METHODS: Participants with suspected or established pSS were enrolled into the Sjögren's International Collaborative Clinical Alliance (SICCA) registry and had oral, ocular, and rheumatologic examinations performed; blood and saliva samples collected; and labial salivary gland biopsy examinations performed; they also completed questionnaires at baseline. Logistic regression was used to determine whether smoking status was associated with pSS classification and individual pSS diagnostic test results. RESULTS: A total of 3514 participants were enrolled in SICCA. A total of 1541 (52.9%) met classification criteria for pSS. Compared with never smokers, current smokers had reduced odds of being classified as having pSS, reduced odds of having a focus score ≥ 1 and serologic positivity for anti-SSA/anti-SSB antibodies, and lower odds of having abnormal signs or test results of dry eye disease. Compared with never smokers, past smokers did not have a statistically significant reduction in odds of being classified as having pSS and of having abnormal individual pSS diagnostic test results. CONCLUSION: Compared with never smokers, current smokers in the SICCA cohort had lower odds of being classified as having pSS, lower odds of exhibiting abnormal signs and test results for dry eye disease, and lower odds of having a labial salivary gland biopsy supportive of pSS classification. Such negative associations, however, do not suggest that current smoking is of any benefit with respect to pSS.

P2RY8 variants in lupus patients uncover a role for the receptor in immunological tolerance.

B cell self-tolerance is maintained through multiple checkpoints, including restraints on intracellular signaling and cell trafficking. P2RY8 is a receptor with established roles in germinal center (GC) B cell migration inhibition and growth regulation. Somatic P2RY8 variants are common in GC-derived B cell lymphomas. Here, we identify germline novel or rare P2RY8 missense variants in lupus kindreds or the related antiphospholipid syndrome, including a "de novo" variant in a child with severe nephritis. All variants decreased protein expression, F-actin abundance, and GPCR-RhoA signaling, and those with stronger effects increased AKT and ERK activity and cell migration. Remarkably, P2RY8 was reduced in B cell subsets from some SLE patients lacking P2RY8 gene variants. Low P2RY8 correlated with lupus nephritis and increased age-associated B cells and plasma cells. By contrast, P2RY8 overexpression in cells and mice restrained plasma cell development and reinforced negative selection of DNA-reactive developing B cells. These findings uncover a role of P2RY8 in immunological tolerance and lupus pathogenesis.

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.].