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.

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.

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.

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.

The Contribution of Genetics and Epigenetics to Our Understanding of Health Disparities in Rheumatic Diseases.

Socioeconomic determinants of health are associated with worse outcomes in the rheumatic diseases and contribute significantly to health disparities. However, genetic and epigenetic risk factors may affect different populations disproportionally and further exacerbate health disparities. We discuss the role of genetics and epigenetics to the health disparities observed in rheumatic diseases. We review concepts of population genetics and natural selection, current genome-wide genetic and epigenetic studies of several autoimmune diseases, and environmental exposures associated with disease risk in different populations. To understand how genomics influence health disparities in the rheumatic diseases, further studies in different populations worldwide are needed.

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.

Concomitant Guillain-Barre syndrome with COVID-19: a case report.

BACKGROUND: During the Coronavirus disease 2019 (COVID-19) pandemic, different neurological manifestations have been observed. However, only a few cases of Guillain-Barre syndrome (GBS) and COVID-19 have been reported. Therefore, the aim of this study is to investigate a case of concomitant GBS with COVID-19 in Colombia. CASE PRESENTATION: A 39-year-old woman was admitted to a teaching hospital in Barranquilla, Colombia with a history of progressive general weakness with lower limb dominance. Previous symptoms such as ageusia, anosmia and intense headache were reported. Upon admission, facial diplegia, quadriparesis with lower extremity predominance and Medical Research Council muscular strength of 2/5 in the lower limbs and 4/5 in the upper limbs were reported. During clinical evolution, due to general areflexia, hypertensive emergency and progressive diaphragmatic weakness, the patient was admitted to an intensive care unit. The cerebrospinal fluid analysis showed protein-cytological dissociation and the GBS diagnosis was confirmed via a nerve conduction and electromyography test. With regard to the symptoms before hospitalisation, SARS-CoV-2 diagnostic testing was performed with positive results in the second test. The patient was managed with supportive care and was discharged after 20 days of hospitalization with clinical improvement. CONCLUSIONS: Only a few cases of COVID-19 with GBS have been reported. Different subtypes have been previously identified, such as Miller-Fisher syndrome and dysautonomic GBS with SARS-CoV-2 infection. This study investigated the first confirmed case of COVID-19 with concomitant GBS in Colombia. In patients with GBS, several viral and bacterial pathogens have been found in case-control studies but the factors that induce the immune-mediated destruction of the nerve tissues have not been determined. Further studies are needed to determine the possible association between COVID-19 exposure and GBS.

Author Correction: Multiplexed droplet single-cell RNA-sequencing using natural genetic variation.

An amendment to this paper has been published and can be accessed via a link at the top of the paper.

Author Correction: A phenotypic and genomics approach in a multi-ethnic cohort to subtype systemic lupus erythematosus.

An amendment to this paper has been published and can be accessed via a link at the top of the paper.

A phenotypic and genomics approach in a multi-ethnic cohort to subtype systemic lupus erythematosus.

Systemic lupus erythematous (SLE) is a heterogeneous autoimmune disease in which outcomes vary among different racial groups. Here, we aim to identify SLE subgroups within a multiethnic cohort using an unsupervised clustering approach based on the American College of Rheumatology (ACR) classification criteria. We identify three patient clusters that vary according to disease severity. Methylation association analysis identifies a set of 256 differentially methylated CpGs across clusters, including 101 CpGs in genes in the Type I Interferon pathway, and we validate these associations in an external cohort. A cis-methylation quantitative trait loci analysis identifies 744 significant CpG-SNP pairs. The methylation signature is enriched for ethnic-associated CpGs suggesting that genetic and non-genetic factors may drive outcomes and ethnic-associated methylation differences. Our computational approach highlights molecular differences associated with clusters rather than single outcome measures. This work demonstrates the utility of applying integrative methods to address clinical heterogeneity in multifactorial multi-ethnic disease settings.

DNA methylation 101: what is important to know about DNA methylation and its role in SLE risk and disease heterogeneity.

SLE is a complex autoimmune disease that results from the interplay of genetics, epigenetics and environmental exposures. DNA methylation is an epigenetic mechanism that regulates gene expression and tissue differentiation. Among all the epigenetic modifications, DNA methylation perturbations have been the most widely studied in SLE. It mediates processes relevant to SLE, including lymphocyte development, X-chromosome inactivation and the suppression of endogenous retroviruses. The establishment of most DNA methylation marks occurs in utero; however, a small percentage of epigenetic marks are dynamic and can change throughout a person's lifetime and in relation to exposures. In this review, we discuss the current understanding of the biology of DNA methylation and its regulators, the measurement and interpretation of methylation marks, the effects of genetics on DNA methylation and the role of environmental exposures with relevance to SLE. We also summarise research findings associated with SLE disease risk and heterogeneity. The robust finding of hypomethylation of interferon-responsive genes in patients with SLE and new associations beyond interferon-responsive genes such as cell-specific methylation abnormalities are described. We also discuss methylation changes associated with lupus nephritis, autoantibody status and disease activity. Lastly, we explore future research directions, emphasising the need for longitudinal studies, cell tissue and context-specific profiling, as well as integrative approaches. With new technologies, DNA methylation perturbations could be targeted and edited, offering novel therapeutic approaches.

Genetic contributions to lupus nephritis in a multi-ethnic cohort of systemic lupus erythematous patients.

OBJECTIVE: African Americans, East Asians, and Hispanics with systemic lupus erythematous (SLE) are more likely to develop lupus nephritis (LN) than are SLE patients of European descent. The etiology of this difference is not clear, and this study was undertaken to investigate how genetic variants might explain this effect. METHODS: In this cross-sectional study, 1244 SLE patients from multiethnic case collections were genotyped for 817,810 single-nucleotide polymorphisms (SNPs) across the genome. Continental genetic ancestry was estimated utilizing the program ADMIXTURE. Gene-based testing and pathway analysis was performed within each ethnic group and meta-analyzed across ethnicities. We also performed candidate SNP association tests with SNPs previously established as risk alleles for SLE, LN, and chronic kidney disease (CKD). Association testing and logistic regression models were performed with LN as the outcome, adjusted for continental ancestries, sex, disease duration, and age. RESULTS: We studied 255 North European, 263 South European, 238 Hispanic, 224 African American and 264 East Asian SLE patients, of whom 606 had LN (48.7%). In genome-wide gene-based and candidate SNP analyses, we found distinct genes, pathways and established risk SNPs associated with LN for each ethnic group. Gene-based analyses showed significant associations between variation in ZNF546 (p = 1.0E-06), TRIM15 (p = 1.0E-06), and TRIMI0 (p = 1.0E-06) and LN among South Europeans, and TTC34 (p = 8.0E-06) was significantly associated with LN among Hispanics. The SNP rs8091180 in NFATC1 was associated with LN (OR 1.43, p = 3.3E-04) in the candidate SNP meta-analysis with the highest OR among African-Americans (OR 2.17, p = 0.0035). CONCLUSION: Distinct genetic factors are associated with the risk of LN in SLE patients of different ethnicities. CKD risk alleles may play a role in the development of LN in addition to SLE-associated risk variants. These findings may further explain the clinical heterogeneity of LN risk and response to therapy observed between different ethnic groups.

Multiplexed droplet single-cell RNA-sequencing using natural genetic variation.

Droplet single-cell RNA-sequencing (dscRNA-seq) has enabled rapid, massively parallel profiling of transcriptomes. However, assessing differential expression across multiple individuals has been hampered by inefficient sample processing and technical batch effects. Here we describe a computational tool, demuxlet, that harnesses natural genetic variation to determine the sample identity of each droplet containing a single cell (singlet) and detect droplets containing two cells (doublets). These capabilities enable multiplexed dscRNA-seq experiments in which cells from unrelated individuals are pooled and captured at higher throughput than in standard workflows. Using simulated data, we show that 50 single-nucleotide polymorphisms (SNPs) per cell are sufficient to assign 97% of singlets and identify 92% of doublets in pools of up to 64 individuals. Given genotyping data for each of eight pooled samples, demuxlet correctly recovers the sample identity of >99% of singlets and identifies doublets at rates consistent with previous estimates. We apply demuxlet to assess cell-type-specific changes in gene expression in 8 pooled lupus patient samples treated with interferon (IFN)-ß and perform eQTL analysis on 23 pooled samples.