Explore the vast range of titles available.

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.

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

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.

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. Systemic lupus erythematosus (SLE) is an autoimmune disease of substantial phenotypic heterogeneity in different ethnic groups. Here, using data from a multi-ethnic cohort, the authors describe an approach based on clinical and molecular data to subtype SLE patients into three clusters of severity.

Systemic lupus erythematosus (SLE) is characterized by the development of autoantibodies associated with specific clinical manifestations. Previous studies have shown an association between differential DNA methylation and SLE susceptibility, but have not investigated SLE-related autoantibodies. Our goal was to determine whether DNA methylation is associated with production of clinically relevant SLE-related autoantibodies, with an emphasis on the anti-dsDNA autoantibody. In this study, we characterized the methylation status of 467,314 CpG sites in 326 women with SLE. Using a discovery and replication study design, we identified and replicated significant associations between anti-dsDNA autoantibody production and the methylation status of 16 CpG sites (pdiscovery<1.07E-07 and preplication<0.0029) in 11 genes. Associations were further investigated using multivariable regression to adjust for estimated leukocyte cell proportions and population substructure. The adjusted mean DNA methylation difference between anti-dsDNA positive and negative cases ranged from 1.2% to 19%, and the adjusted odds ratio for anti-dsDNA autoantibody production comparing the lowest and highest methylation tertiles ranged from 6.8 to 18.2. Differential methylation for these CpG sites was also associated with anti-SSA, anti-Sm, and anti-RNP autoantibody production. Overall, associated CpG sites were hypomethylated in autoantibody positive compared to autoantibody negative cases. Differential methylation of CpG sites within the major histocompatibility region was not strongly associated with autoantibody production. Genes with differentially methylated CpG sites represent multiple biologic pathways, and have not been associated with autoantibody production in genetic association studies. In conclusion, hypomethylation of CpG sites within genes from different pathways is associated with anti-dsDNA, anti-SSA, anti-Sm, and anti-RNP production in SLE, and these associations are not explained by genetic variation. Thus, studies of epigenetic mechanisms such as DNA methylation represent a complementary method to genetic association studies to identify biologic pathways that may contribute to the clinical heterogeneity of autoimmune diseases.

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.

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

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⁻¹²⁸) 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⁻⁸. 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.

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.

Background Systemic lupus erythematosus (SLE) has numerous symptoms across organs and an unpredictable flare-remittance pattern. This has made it challenging to understand drivers of long-term SLE outcomes. Our objective was to identify whether changes in DNA methylation over time, in an actively flaring SLE cohort, were associated with remission and whether these changes meaningfully subtype SLE patients. Methods Fifty-nine multi-ethnic SLE patients had clinical visits and DNA methylation profiles at a flare and approximately 3 months later. Methylation was measured using the Illumina EPIC array. We identified sites where methylation change between visits was associated with remission at the follow-up visit using limma package and a time x remission interaction term. Models adjusted for batch, age at diagnosis, time between visits, age at flare, sex, medications, and cell-type proportions. Separately, a paired T-test identified Bonferroni significant methylation sites with ≥ 3% change between visits ( n  = 546). Methylation changes at these sites were used for unsupervised consensus hierarchical clustering. Associations between clusters and patient features were assessed. Results Nineteen patients fully remitted at the follow-up visit. For 1,953 CpG sites, methylation changed differently for remitters vs. non-remitters (Bonferroni p  < 0.05). Nearly half were within genes regulated by interferon. The largest effect was at cg22873177; on average, remitters had 23% decreased methylation between visits while non-remitters had no change. Three SLE patient clusters were identified using methylation differences agnostic of clinical outcomes. All Cluster 1 subjects ( n  = 12) experienced complete flare remission, despite similar baseline disease activity scores, medications, and demographics as other clusters. Methylation changes at six CpG sites, including within immune-related CD45 and IFI genes, were particularly distinct for each cluster, suggesting these may be good candidates for stratifying patients in the future. Conclusions Changes in DNA methylation during active SLE were associated with remission status and identified subgroups of SLE patients with several distinct clinical and biological characteristics. DNA methylation patterns might help inform SLE subtypes, leading to targeted therapies based on relevant underlying biological pathways.