Integration of multi-omics analysis reveals metabolic alterations of B lymphocytes in systemic lupus erythematosus.

OBJECTIVE: To link changes in the B-cell transcriptome from systemic lupus erythematosus (SLE) patients with those in their macroenvironment, including cellular and fluidic components. METHODS: Analysis was performed on 363 patients and 508 controls, encompassing transcriptomics, metabolomics, and clinical data. B-cell and whole-blood transcriptomes were analysed using DESeq and GSEA. Plasma and urine metabolomics peak changes were quantified and annotated using Ceu Mass Mediator database. Common sources of variation were identified using MOFA integration analysis. RESULTS: Cellular macroenvironment was enriched in cytokines, stress responses, lipidic synthesis/mobility pathways and nucleotide degradation. B cells shared these pathways, except nucleotide degradation diverted to nucleotide salvage pathway, and distinct glycosylation, LPA receptors and Schlafen proteins. CONCLUSIONS: B cells showed metabolic changes shared with their macroenvironment and unique changes directly or indirectly induced by IFN-α signalling. This study underscores the importance of understanding the interplay between B cells and their macroenvironment in SLE pathology.

Scoring personalized molecular portraits identify Systemic Lupus Erythematosus subtypes and predict individualized drug responses, symptomatology and disease progression.

OBJECTIVES: Systemic Lupus Erythematosus is a complex autoimmune disease that leads to significant worsening of quality of life and mortality. Flares appear unpredictably during the disease course and therapies used are often only partially effective. These challenges are mainly due to the molecular heterogeneity of the disease, and in this context, personalized medicine-based approaches offer major promise. With this work we intended to advance in that direction by developing MyPROSLE, an omic-based analytical workflow for measuring the molecular portrait of individual patients to support clinicians in their therapeutic decisions. METHODS: Immunological gene-modules were used to represent the transcriptome of the patients. A dysregulation score for each gene-module was calculated at the patient level based on averaged z-scores. Almost 6100 Lupus and 750 healthy samples were used to analyze the association among dysregulation scores, clinical manifestations, prognosis, flare and remission events and response to Tabalumab. Machine learning-based classification models were built to predict around 100 different clinical parameters based on personalized dysregulation scores. RESULTS: MyPROSLE allows to molecularly summarize patients in 206 gene-modules, clustered into nine main lupus signatures. The combination of these modules revealed highly differentiated pathological mechanisms. We found that the dysregulation of certain gene-modules is strongly associated with specific clinical manifestations, the occurrence of relapses or the presence of long-term remission and drug response. Therefore, MyPROSLE may be used to accurately predict these clinical outcomes. CONCLUSIONS: MyPROSLE (https://myprosle.genyo.es) allows molecular characterization of individual Lupus patients and it extracts key molecular information to support more precise therapeutic decisions.

Molecular characterisation of lupus low disease activity state (LLDAS) and DORIS remission by whole-blood transcriptome-based pathways in a pan-European systemic lupus erythematosus cohort.

OBJECTIVES: To unveil biological milieus underlying low disease activity (LDA) and remission versus active systemic lupus erythematosus (SLE). METHODS: We determined differentially expressed pathways (DEPs) in SLE patients from the PRECISESADS project (NTC02890121) stratified into patients fulfilling and not fulfilling the criteria of (1) Lupus LDA State (LLDAS), (2) Definitions of Remission in SLE remission, and (3) LLDAS exclusive of remission. RESULTS: We analysed data from 321 patients; 40.8% were in LLDAS, and 17.4% in DORIS remission. After exclusion of patients in remission, 28.3% were in LLDAS. Overall, 604 pathways differed significantly in LLDAS versus non-LLDAS patients with an false-discovery rate-corrected p (q)<0.05 and a robust effect size (dr)≥0.36. Accordingly, 288 pathways differed significantly between DORIS remitters and non-remitters (q<0.05 and dr≥0.36). DEPs yielded distinct molecular clusters characterised by differential serological, musculoskeletal, and renal activity. Analysis of partially overlapping samples showed no DEPs between LLDAS and DORIS remission. Drug repurposing potentiality for treating SLE was unveiled, as were important pathways underlying active SLE whose modulation could aid attainment of LLDAS/remission, including toll-like receptor (TLR) cascades, Bruton tyrosine kinase (BTK) activity, the cytotoxic T lymphocyte antigen 4 (CTLA-4)-related inhibitory signalling, and the nucleotide-binding oligomerization domain leucine-rich repeat-containing protein 3 (NLRP3) inflammasome pathway. CONCLUSIONS: We demonstrated for the first time molecular signalling pathways distinguishing LLDAS/remission from active SLE. LLDAS/remission was associated with reversal of biological processes related to SLE pathogenesis and specific clinical manifestations. DEP clustering by remission better grouped patients compared with LLDAS, substantiating remission as the ultimate treatment goal in SLE; however, the lack of substantial pathway differentiation between the two states justifies LLDAS as an acceptable goal from a biological perspective.

Distinct gene dysregulation patterns herald precision medicine potentiality in systemic lupus erythematosus.

OBJECTIVES: We aimed at investigating the whole-blood transcriptome, expression quantitative trait loci (eQTLs), and levels of selected serological markers in patients with SLE versus healthy controls (HC) to gain insight into pathogenesis and identify drug targets. METHODS: We analyzed differentially expressed genes (DEGs) and dysregulated gene modules in a cohort of 350 SLE patients and 497 HC from the European PRECISESADS project (NTC02890121), split into a discovery (60%) and a replication (40%) set. Replicated DEGs qualified for eQTL, pathway enrichment, regulatory network, and druggability analysis. For validation purposes, a separate gene module analysis was performed in an independent cohort (GSE88887). RESULTS: Analysis of 521 replicated DEGs identified multiple enriched interferon signaling pathways through Reactome. Gene module analysis yielded 18 replicated gene modules in SLE patients, including 11 gene modules that were validated in GSE88887. Three distinct gene module clusters were defined i.e., "interferon/plasma cells", "inflammation", and "lymphocyte signaling". Predominant downregulation of the lymphocyte signaling cluster denoted renal activity. By contrast, upregulation of interferon-related genes indicated hematological activity and vasculitis. Druggability analysis revealed several potential drugs interfering with dysregulated genes within the "interferon" and "PLK1 signaling events" modules. STAT1 was identified as the chief regulator in the most enriched signaling molecule network. Drugs annotated to 15 DEGs associated with cis-eQTLs included bortezomib for its ability to modulate CTSL activity. Belimumab was annotated to TNFSF13B (BAFF) and daratumumab was annotated to CD38 among the remaining replicated DEGs. CONCLUSIONS: Modulation of interferon, STAT1, PLK1, B and plasma cell signatures showed promise as viable approaches to treat SLE, pointing to their importance in SLE pathogenesis.

A comprehensive database for integrated analysis of omics data in autoimmune diseases.

BACKGROUND: Autoimmune diseases are heterogeneous pathologies with difficult diagnosis and few therapeutic options. In the last decade, several omics studies have provided significant insights into the molecular mechanisms of these diseases. Nevertheless, data from different cohorts and pathologies are stored independently in public repositories and a unified resource is imperative to assist researchers in this field. RESULTS: Here, we present Autoimmune Diseases Explorer ( https://adex.genyo.es ), a database that integrates 82 curated transcriptomics and methylation studies covering 5609 samples for some of the most common autoimmune diseases. The database provides, in an easy-to-use environment, advanced data analysis and statistical methods for exploring omics datasets, including meta-analysis, differential expression or pathway analysis. CONCLUSIONS: This is the first omics database focused on autoimmune diseases. This resource incorporates homogeneously processed data to facilitate integrative analyses among studies.

Genome-wide whole blood transcriptome profiling in a large European cohort of systemic sclerosis patients.

OBJECTIVES: The analysis of annotated transcripts from genome-wide expression studies may help to understand the pathogenesis of complex diseases, such as systemic sclerosis (SSc). We performed a whole blood (WB) transcriptome analysis on RNA collected in the context of the European PRECISESADS project, aiming at characterising the pathways that differentiate SSc from controls and that are reproducible in geographically diverse populations. METHODS: Samples from 162 patients and 252 controls were collected in RNA stabilisers. Cases and controls were divided into a discovery (n=79+163; Southern Europe) and validation cohort (n=83+89; Central-Western Europe). RNA sequencing was performed by an Illumina assay. Functional annotations of Reactome pathways were performed with the Functional Analysis of Individual Microarray Expression (FAIME) algorithm. In parallel, immunophenotyping of 28 circulating cell populations was performed. We tested the presence of differentially expressed genes/pathways and the correlation between absolute cell counts and RNA transcripts/FAIME scores in regression models. Results significant in both populations were considered as replicated. RESULTS: Overall, 15 224 genes and 1277 functional pathways were available; of these, 99 and 225 were significant in both sets. Among replicated pathways, we found a deregulation in type-I interferon, Toll-like receptor cascade, tumour suppressor p53 protein function, platelet degranulation and activation. RNA transcripts or FAIME scores were jointly correlated with cell subtypes with strong geographical differences; neutrophils were the major determinant of gene expression in SSc-WB samples. CONCLUSIONS: We discovered a set of differentially expressed genes/pathways validated in two independent sets of patients with SSc, highlighting a number of deregulated processes that have relevance for the pathogenesis of autoimmunity and SSc.

NGSmethDB 2017: enhanced methylomes and differential methylation.

The 2017 update of NGSmethDB stores whole genome methylomes generated from short-read data sets obtained by bisulfite sequencing (WGBS) technology. To generate high-quality methylomes, stringent quality controls were integrated with third-part software, adding also a two-step mapping process to exploit the advantages of the new genome assembly models. The samples were all profiled under constant parameter settings, thus enabling comparative downstream analyses. Besides a significant increase in the number of samples, NGSmethDB now includes two additional data-types, which are a valuable resource for the discovery of methylation epigenetic biomarkers: (i) differentially methylated single-cytosines; and (ii) methylation segments (i.e. genome regions of homogeneous methylation). The NGSmethDB back-end is now based on MongoDB, a NoSQL hierarchical database using JSON-formatted documents and dynamic schemas, thus accelerating sample comparative analyses. Besides conventional database dumps, track hubs were implemented, which improved database access, visualization in genome browsers and comparative analyses to third-part annotations. In addition, the database can be also accessed through a RESTful API. Lastly, a Python client and a multiplatform virtual machine allow for program-driven access from user desktop. This way, private methylation data can be compared to NGSmethDB without the need to upload them to public servers. Database website: http://bioinfo2.ugr.es/NGSmethDB.

sRNAtoolbox: an integrated collection of small RNA research tools.

Small RNA research is a rapidly growing field. Apart from microRNAs, which are important regulators of gene expression, other types of functional small RNA molecules have been reported in animals and plants. MicroRNAs are important in host-microbe interactions and parasite microRNAs might modulate the innate immunity of the host. Furthermore, small RNAs can be detected in bodily fluids making them attractive non-invasive biomarker candidates. Given the general broad interest in small RNAs, and in particular microRNAs, a large number of bioinformatics aided analysis types are needed by the scientific community. To facilitate integrated sRNA research, we developed sRNAtoolbox, a set of independent but interconnected tools for expression profiling from high-throughput sequencing data, consensus differential expression, target gene prediction, visual exploration in a genome context as a function of read length, gene list analysis and blast search of unmapped reads. All tools can be used independently or for the exploration and downstream analysis of sRNAbench results. Workflows like the prediction of consensus target genes of parasite microRNAs in the host followed by the detection of enriched pathways can be easily established. The web-interface interconnecting all these tools is available at http://bioinfo5.ugr.es/srnatoolbox.

NGSmethDB: an updated genome resource for high quality, single-cytosine resolution methylomes.

The updated release of 'NGSmethDB' (http://bioinfo2.ugr.es/NGSmethDB) is a repository for single-base whole-genome methylome maps for the best-assembled eukaryotic genomes. Short-read data sets from NGS bisulfite-sequencing projects of cell lines, fresh and pathological tissues are first pre-processed and aligned to the corresponding reference genome, and then the cytosine methylation levels are profiled. One major improvement is the application of a unique bioinformatics protocol to all data sets, thereby assuring the comparability of all values with each other. We implemented stringent quality controls to minimize important error sources, such as sequencing errors, bisulfite failures, clonal reads or single nucleotide variants (SNVs). This leads to reliable and high-quality methylomes, all obtained under uniform settings. Another significant improvement is the detection in parallel of SNVs, which might be crucial for many downstream analyses (e.g. SNVs and differential-methylation relationships). A next-generation methylation browser allows fast and smooth scrolling and zooming, thus speeding data download/upload, at the same time requiring fewer server resources. Several data mining tools allow the comparison/retrieval of methylation levels in different tissues or genome regions. NGSmethDB methylomes are also available as native tracks through a UCSC hub, which allows comparison with a wide range of third-party annotations, in particular phenotype or disease annotations.

Interferon and B-cell Signatures Inform Precision Medicine in Lupus Nephritis.

Introduction: Current therapeutic management of lupus nephritis (LN) fails to induce long-term remission in over 50% of patients, highlighting the urgent need for additional options. Methods: We analyzed differentially expressed genes (DEGs) in peripheral blood from patients with active LN (n = 41) and active nonrenal lupus (n = 62) versus healthy controls (HCs) (n = 497) from the European PRECISESADS project (NTC02890121), and dysregulated gene modules in a discovery (n = 26) and a replication (n = 15) set of active LN cases. Results: Replicated gene modules qualified for correlation analyses with serologic markers, and regulatory network and druggability analysis. Unsupervised coexpression network analysis revealed 20 dysregulated gene modules and stratified the active LN population into 3 distinct subgroups. These subgroups were characterized by low, intermediate, and high interferon (IFN) signatures, with differential dysregulation of the "B cell" and "plasma cells/Ig" modules. Drugs annotated to the IFN network included CC-motif chemokine receptor 1 (CCR1) inhibitors, programmed death-ligand 1 (PD-L1) inhibitors, and irinotecan; whereas the anti-CD38 daratumumab and proteasome inhibitor bortezomib showed potential for counteracting the "plasma cells/Ig" signature. In silico analysis demonstrated the low-IFN subgroup to benefit from calcineurin inhibition and the intermediate-IFN subgroup from B-cell targeted therapies. High-IFN patients exhibited greater anticipated response to anifrolumab whereas daratumumab appeared beneficial to the intermediate-IFN and high-IFN subgroups. Conclusion: IFN upregulation and B and plasma cell gene dysregulation patterns revealed 3 subgroups of LN, which may not necessarily represent distinct disease phenotypes but rather phases of the inflammatory processes during a renal flare, providing a conceptual framework for precision medicine in LN.

Molecular subtypes explain lupus epigenomic heterogeneity unveiling new regulatory genetic risk variants.

The heterogeneity of systemic lupus erythematosus (SLE) can be explained by epigenetic alterations that disrupt transcriptional programs mediating environmental and genetic risk. This study evaluated the epigenetic contribution to SLE heterogeneity considering molecular and serological subtypes, genetics and transcriptional status, followed by drug target discovery. We performed a stratified epigenome-wide association studies of whole blood DNA methylation from 213 SLE patients and 221 controls. Methylation quantitative trait loci analyses, cytokine and transcription factor activity - epigenetic associations and methylation-expression correlations were conducted. New drug targets were searched for based on differentially methylated genes. In a stratified approach, a total of 974 differential methylation CpG sites with dependency on molecular subtypes and autoantibody profiles were found. Mediation analyses suggested that SLE-associated SNPs in the HLA region exert their risk through DNA methylation changes. Novel genetic variants regulating DNAm in disease or in specific molecular contexts were identified. The epigenetic landscapes showed strong association with transcription factor activity and cytokine levels, conditioned by the molecular context. Epigenetic signals were enriched in known and novel drug targets for SLE. This study reveals possible genetic drivers and consequences of epigenetic variability on SLE heterogeneity and disentangles the DNAm mediation role on SLE genetic risk and novel disease-specific meQTLs. Finally, novel targets for drug development were discovered.

NGSmethDB: a database for next-generation sequencing single-cytosine-resolution DNA methylation data.

Next-generation sequencing (NGS) together with bisulphite conversion allows the generation of whole genome methylation maps at single-cytosine resolution. This allows studying the absence of methylation in a particular genome region over a range of tissues, the differential tissue methylation or the changes occurring along pathological conditions. However, no database exists fully addressing such requirements. We propose here NGSmethDB (http://bioinfo2.ugr.es/NGSmethDB/gbrowse/) for the storage and retrieval of methylation data derived from NGS. Two cytosine methylation contexts (CpG and CAG/CTG) are considered. Through a browser interface coupled to a MySQL backend and several data mining tools, the user can search for methylation states in a set of tissues, retrieve methylation values for a set of tissues in a given chromosomal region, or display the methylation of promoters among different tissues. NGSmethDB is currently populated with human, mouse and Arabidopsis data, but other methylomes will be incorporated through an automatic pipeline as soon as new data become available. Dump downloads for three coverage levels (1, 5 or 10 reads) are available. NGSmethDB will be useful for experimental researchers, as well as for bioinformaticians, who might use the data as input for further research.

Clustering of DNA words and biological function: a proof of principle.

Relevant words in literary texts (key words) are known to be clustered, while common words are randomly distributed. Given the clustered distribution of many functional genome elements, we hypothesize that the biological text per excellence, the DNA sequence, might behave in the same way: k-length words (k-mers) with a clear function may be spatially clustered along the one-dimensional chromosome sequence, while less-important, non-functional words may be randomly distributed. To explore this linguistic analogy, we calculate a clustering coefficient for each k-mer (k=2-9bp) in human and mouse chromosome sequences, then checking if clustered words are enriched in the functional part of the genome. First, we found a positive general trend relating clustering level and word enrichment within exons and Transcription Factor Binding Sites (TFBSs), while a much weaker relation exists for repeats, and no relation at all exists for introns. Second, we found that 38.45% of the 200 top-clustered 8-mers, but only 7.70% of the non-clustered words, are represented in known motif databases. Third, enrichment/depletion experiments show that highly clustered words are significantly enriched in exons and TFBSs, while they are depleted in introns and repetitive DNA. Considering exons and TFBSs together, 1417 (or 72.26%) in human and 1385 (or 72.97%) in mouse of the top-clustered 8-mers showed a statistically significant association to either exons or TFBSs, thus strongly supporting the link between word clustering and biological function. Lastly, we identified a subset of clustered, diagnostic words that are enriched in exons but depleted in introns, and therefore might help to discriminate between these two gene regions. The clustering of DNA words thus appears as a novel principle to detect functionality in genome sequences. As evolutionary conservation is not a prerequisite, the proof of principle described here may open new ways to detect species-specific functional DNA sequences and the improvement of gene and promoter predictions, thus contributing to the quest for function in the genome.

Whole blood DNA methylation analysis reveals respiratory environmental traits involved in COVID-19 severity following SARS-CoV-2 infection.

SARS-CoV-2 infection can cause an inflammatory syndrome (COVID-19) leading, in many cases, to bilateral pneumonia, severe dyspnea, and in ~5% of these, death. DNA methylation is known to play an important role in the regulation of the immune processes behind COVID-19 progression, however it has not been studied in depth. In this study, we aim to evaluate the implication of DNA methylation in COVID-19 progression by means of a genome-wide DNA methylation analysis combined with DNA genotyping. The results reveal the existence of epigenomic regulation of functional pathways associated with COVID-19 progression and mediated by genetic loci. We find an environmental trait-related signature that discriminates mild from severe cases and regulates, among other cytokines, IL-6 expression via the transcription factor CEBP. The analyses suggest that an interaction between environmental contribution, genetics, and epigenetics might be playing a role in triggering the cytokine storm described in the most severe cases.

geno5mC: A Database to Explore the Association between Genetic Variation (SNPs) and CpG Methylation in the Human Genome.

Genetic variation, gene expression and DNA methylation influence each other in a complex way. To study the impact of sequence variation and DNA methylation on gene expression, we generated geno5mC, a database that contains statistically significant SNP-CpG associations that are biologically classified either through co-localization with known regulatory regions (promoters and enhancers), or through known correlations with the expression levels of nearby genes. The SNP rs727563 can be used to illustrate the usefulness of this approach. This SNP has been associated with inflammatory bowel disease through GWAS, but it is not located near any gene related to this phenotype. However, geno5mC reveals that rs727563 is associated with the methylation state of several CpGs located in promoter regions of genes reported to be involved in inflammatory processes. This case exemplifies how geno5mC can be used to infer relevant and previously unknown interactions between described disease-associated SNPs and their functional targets.

An elevated polyclonal free light chain level reflects a strong interferon signature in patients with systemic autoimmune diseases.

High amount of polyclonal free light chains (FLC) are reported in systemic autoimmune diseases (SAD) and we took advantage of the PRECISESADS study to better characterize them. Serum FLC levels were explored in 1979 patients with SAD (RA, SLE, SjS, Scl, APS, UCTD, MCTD) and 614 healthy controls. Information regarding clinical parameters, disease activity, medications, autoantibodies (Ab) and the interferon α and/or γ scores were recorded. Among SAD patients, 28.4% had raised total FLC (from 12% in RA to 30% in SLE and APS) with a normal kappa/lambda ratio. Total FLC levels were significantly higher in SAD with inflammation, active disease in SLE and SjS, and an impaired pulmonary functional capacity in SSc, while independent from kidney impairment, infection, cancer and treatment. Total FLC concentrations were positively correlated among the 10/17 (58.8%) autoantibodies (Ab) tested with anti-RNA binding protein Ab (SSB, SSA-52/60 kDa, Sm, U1-RNP), anti-dsDNA/nucleosome Ab, rheumatoid factor and negatively correlated with complement fractions C3/C4. Finally, examination of interferon (IFN) expression as a potential driver of FLC overexpression was tested showing an elevated level of total FLC among patients with a high IFNα and IFNγ Kirou's score, a strong IFN modular score, and the detection in the sera of B-cell IFN dependent factors, such as TNF-R1/TNFRSF1A and CXCL10/IP10. In conclusion, an elevated level of FLC, in association with a strong IFN signature, defines a subgroup of SAD patients, including those without renal affectation, characterized by increased disease activity, autoreactivity, and complement reduction.

Expression Quantitative Trait Locus Analysis in Systemic Sclerosis Identifies New Candidate Genes Associated With Multiple Aspects of Disease Pathology.

OBJECTIVE: To identify the genetic variants that affect gene expression (expression quantitative trait loci [eQTLs]) in systemic sclerosis (SSc) and to investigate their role in the pathogenesis of the disease. METHODS: We performed an eQTL analysis using whole-blood sequencing data from 333 SSc patients and 524 controls and integrated them with SSc genome-wide association study (GWAS) data. We integrated our findings from expression modeling, differential expression analysis, and transcription factor binding site enrichment with key clinical features of SSc. RESULTS: We detected 49,123 validated cis-eQTLs from 4,539 SSc-associated single-nucleotide polymorphisms (SNPs) (PGWAS < 10-5 ). A total of 1,436 genes were within 1 Mb of the 4,539 SSc-associated SNPs. Of those 1,436 genes, 565 were detected as having ≥1 eQTL with an SSc-associated SNP. We developed a strategy to prioritize disease-associated genes based on their expression variance explained by SSc eQTLs (r2 > 0.05). As a result, 233 candidates were identified, 134 (58%) of them associated with hallmarks of SSc and 105 (45%) of them differentially expressed in the blood cells, skin, or lung tissue of SSc patients. Transcription factor binding site analysis revealed enriched motifs of 24 transcription factors (5%) among SSc eQTLs, 5 of which were found to be differentially regulated in the blood cells (ELF1 and MGA), skin (KLF4 and ID4), and lungs (TBX4) of SSc patients. Ten candidate genes (4%) can be targeted by approved medications for immune-mediated diseases, of which only 3 have been tested in clinical trials in patients with SSc. CONCLUSION: The findings of the present study indicate a new layer to the molecular complexity of SSc, contributing to a better understanding of the pathogenesis of the disease.

Integrative epigenomics in Sjögren´s syndrome reveals novel pathways and a strong interaction between the HLA, autoantibodies and the interferon signature.

Primary Sjögren's syndrome (SS) is a systemic autoimmune disease characterized by lymphocytic infiltration and damage of exocrine salivary and lacrimal glands. The etiology of SS is complex with environmental triggers and genetic factors involved. By conducting an integrated multi-omics study, we confirmed a vast coordinated hypomethylation and overexpression effects in IFN-related genes, what is known as the IFN signature. Stratified and conditional analyses suggest a strong interaction between SS-associated HLA genetic variation and the presence of Anti-Ro/SSA autoantibodies in driving the IFN epigenetic signature and determining SS. We report a novel epigenetic signature characterized by increased DNA methylation levels in a large number of genes enriched in pathways such as collagen metabolism and extracellular matrix organization. We identified potential new genetic variants associated with SS that might mediate their risk by altering DNA methylation or gene expression patterns, as well as disease-interacting genetic variants that exhibit regulatory function only in the SS population. Our study sheds new light on the interaction between genetics, autoantibody profiles, DNA methylation and gene expression in SS, and contributes to elucidate the genetic architecture of gene regulation in an autoimmune population.