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.

Consensus gene modules strategy identifies candidate blood-based biomarkers for primary Sjögren's disease.

Primary Sjögren disease (pSD) is an autoimmune disease characterized by lymphoid infiltration of exocrine glands leading to dryness of the mucosal surfaces and by the production of autoantibodies. The pathophysiology of pSD remains elusive and no treatment with demonstrated efficacy is available yet. To better understand the biology underlying pSD heterogeneity, we aimed at identifying Consensus gene Modules (CMs) that summarize the high-dimensional transcriptomic data of whole blood samples in pSD patients. We performed unsupervised gene classification on four data sets and identified thirteen CMs. We annotated and interpreted each of these CMs as corresponding to cell type abundances or biological functions by using gene set enrichment analyses and transcriptomic profiles of sorted blood cell subsets. Correlation with independently measured cell type abundances by flow cytometry confirmed these annotations. We used these CMs to reconcile previously proposed patient stratifications of pSD. Importantly, we showed that the expression of modules representing lymphocytes and erythrocytes before treatment initiation is associated with response to hydroxychloroquine and leflunomide combination therapy in a clinical trial. These consensus modules will help the identification and translation of blood-based predictive biomarkers for the treatment of pSD.

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.

Stratification according to autoantibody status in systemic sclerosis reveals distinct molecular signatures.

OBJECTIVES: Systemic sclerosis (SSc) is a heterogeneous disease, complicating its management. Its complexity and the insufficiency of clinical manifestations alone to delineate homogeneous patient groups further challenge this task. However, autoantibodies could serve as relevant markers for the pathophysiological mechanisms driving the disease. Identifying specific immunological mechanisms based on patients' serological statuses might facilitate a deeper understanding of the diversity of the disease. METHODS: A cohort of 206 patients with SSc enrolled in the PRECISESADS cross-sectional study was examined. Patients were stratified based on their anti-centromere (ACA) and anti-SCL70 (SCL70) antibody statuses. Comprehensive omics analyses including transcriptomic, flow cytometric, cytokine and metabolomic data were analysed to characterise the differences between these patient groups. RESULTS: Patients with SCL70 antibodies showed severe clinical features such as diffuse cutaneous sclerosis and pulmonary fibrosis and were biologically distinguished by unique transcriptomic profiles. They exhibit a pro-inflammatory and fibrotic signature associated with impaired tissue remodelling and increased carnitine metabolism. Conversely, ACA-positive patients exhibited an immunomodulation and tissue homeostasis signature and increased phospholipid metabolism. CONCLUSIONS: Patients with SSc display varying biological profiles based on their serological status. The findings highlight the potential utility of serological status as a discriminating factor in disease severity and suggest its relevance in tailoring treatment strategies and future research directions.

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.

Continental-scale genomic analysis suggests shared post-admixture adaptation in the Americas.

American populations are one of the most interesting examples of recently admixed groups, where ancestral components from three major continental human groups (Africans, Eurasians and Native Americans) have admixed within the last 15 generations. Recently, several genetic surveys focusing on thousands of individuals shed light on the geography, chronology and relevance of these events. However, even though gene flow could drive adaptive evolution, it is unclear whether and how natural selection acted on the resulting genetic variation in the Americas. In this study, we analysed the patterns of local ancestry of genomic fragments in genome-wide data for ~ 6000 admixed individuals from 10 American countries. In doing so, we identified regions characterized by a divergent ancestry profile (DAP), in which a significant over or under ancestral representation is evident. Our results highlighted a series of genomic regions with DAPs associated with immune system response and relevant medical traits, with the longest DAP region encompassing the human leukocyte antigen locus. Furthermore, we found that DAP regions are enriched in genes linked to cancer-related traits and autoimmune diseases. Then, analysing the biological impact of these regions, we showed that natural selection could have acted preferentially towards variants located in coding and non-coding transcripts and characterized by a high deleteriousness score. Taken together, our analyses suggest that shared patterns of post admixture adaptation occurred at a continental scale in the Americas, affecting more often functional and impactful genomic variants.

Single-cell immune profiling of Meniere Disease patients.

BACKGROUND: Meniere Disease (MD) is an inner ear syndrome, characterized by episodes of vertigo, tinnitus and fluctuating sensorineural hearing loss. The pathological mechanism leading to sporadic MD is still poorly understood, however an allergic inflammatory response seems to be involved in some patients with MD. OBJECTIVE: Decipher an immune signature associated with the syndrome. METHODS: We performed mass cytometry immune profiling on peripheral blood from MD patients and controls. We analyzed differences in state and differences in abundance of the different cellular subsets. IgE levels were quantified through ELISA on supernatant of cultured whole blood. RESULTS: We have identified two clusters of individuals according to the single cell cytokine profile. These clusters presented differences in IgE levels, immune cell population abundance, including a reduction of CD56dim NK-cells, and changes in cytokine expression with a different response to bacterial and fungal antigens. CONCLUSION: Our results support a systemic inflammatory response in some MD patients that show a type 2 response with allergic phenotype, which could benefit from personalized IL-4 blockers.

A survey of gene expression meta-analysis: methods and applications.

The increasing use of high-throughput gene expression quantification technologies over the last two decades and the fact that most of the published studies are stored in public databases has triggered an explosion of studies available through public repositories. All this information offers an invaluable resource for reuse to generate new knowledge and scientific findings. In this context, great interest has been focused on meta-analysis methods to integrate and jointly analyze different gene expression datasets. In this work, we describe the main steps in the gene expression meta-analysis, from data preparation to the state-of-the art statistical methods. We also analyze the main types of applications and problems that can be approached in gene expression meta-analysis studies and provide a comparative overview of the available software and bioinformatics tools. Moreover, a practical guide for choosing the most appropriate method in each case is also provided.

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.

PeacoQC: Peak-based selection of high quality cytometry data.

In cytometry analysis, a large number of markers is measured for thousands or millions of cells, resulting in high-dimensional datasets. During the measurement of these samples, erroneous events can occur such as clogs, speed changes, slow uptake of the sample etc., which can influence the downstream analysis and can even lead to false discoveries. As these issues can be difficult to detect manually, an automated approach is recommended. In order to filter these erroneous events out, we created a novel quality control algorithm, Peak Extraction And Cleaning Oriented Quality Control (PeacoQC), that allows for automated cleaning of cytometry data. The algorithm will determine density peaks per channel on which it will remove low quality events based on their position in the isolation tree and on their mean absolute deviation distance to these density peaks. To evaluate PeacoQC's cleaning capability, it was compared to three other existing quality control algorithms (flowAI, flowClean and flowCut) on a wide variety of datasets. In comparison to the other algorithms, PeacoQC was able to filter out all different types of anomalies in flow, mass and spectral cytometry data, while the other methods struggled with at least one type. In the quantitative comparison, PeacoQC obtained the highest median balanced accuracy and a similar running time compared to the other algorithms while having a better scalability for large files. To ensure that the parameters chosen in the PeacoQC algorithm are robust, the cleaning tool was run on 16 public datasets. After inspection, only one sample was found where the parameters should be further optimized. The other 15 datasets were analyzed correctly indicating a robust parameter choice. Overall, we present a fast and accurate quality control algorithm that outperforms existing tools and ensures high-quality data that can be used for further downstream analysis. An R implementation is available.

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.

Genomic Risk Score impact on susceptibility to systemic sclerosis.

OBJECTIVES: Genomic Risk Scores (GRS) successfully demonstrated the ability of genetics to identify those individuals at high risk for complex traits including immune-mediated inflammatory diseases (IMIDs). We aimed to test the performance of GRS in the prediction of risk for systemic sclerosis (SSc) for the first time. METHODS: Allelic effects were obtained from the largest SSc Genome-Wide Association Study (GWAS) to date (9 095 SSc and 17 584 healthy controls with European ancestry). The best-fitting GRS was identified under the additive model in an independent cohort that comprised 400 patients with SSc and 571 controls. Additionally, GRS for clinical subtypes (limited cutaneous SSc and diffuse cutaneous SSc) and serological subtypes (anti-topoisomerase positive (ATA+) and anti-centromere positive (ACA+)) were generated. We combined the estimated GRS with demographic and immunological parameters in a multivariate generalised linear model. RESULTS: The best-fitting SSc GRS included 33 single nucleotide polymorphisms (SNPs) and discriminated between patients with SSc and controls (area under the receiver operating characteristic (ROC) curve (AUC)=0.673). Moreover, the GRS differentiated between SSc and other IMIDs, such as rheumatoid arthritis and Sjögren's syndrome. Finally, the combination of GRS with age and immune cell counts significantly increased the performance of the model (AUC=0.787). While the SSc GRS was not able to discriminate between ATA+ and ACA+ patients (AUC<0.5), the serological subtype GRS, which was based on the allelic effects observed for the comparison between ACA+ and ATA+ patients, reached an AUC=0.693. CONCLUSIONS: GRS was successfully implemented in SSc. The model discriminated between patients with SSc and controls or other IMIDs, confirming the potential of GRS to support early and differential diagnosis for SSc.

Stabilization of Human Whole Blood Samples for Multicenter and Retrospective Immunophenotyping Studies.

Whole blood is often collected for large-scale immune monitoring studies to track changes in cell frequencies and responses using flow (FC) or mass cytometry (MC). In order to preserve sample composition and phenotype, blood samples should be analyzed within 24 h after bleeding, restricting the recruitment, analysis protocols, as well as biobanking. Herein, we have evaluated two whole blood preservation protocols that allow rapid sample processing and long-term stability. Two fixation buffers were used, Phosphoflow Fix and Lyse (BD) and Proteomic Stabilizer (PROT) to fix and freeze whole blood samples for up to 6 months. After analysis by an 8-plex panel by FC and a 26-plex panel by MC, manual gating of circulating leukocyte populations and cytokines was performed. Additionally, we tested the stability of a single sample over a 13-months period using 45 consecutive aliquots and a 34-plex panel by MC. We observed high correlation and low bias toward any cell population when comparing fresh and 6 months frozen blood with FC and MC. This correlation was confirmed by hierarchical clustering. Low coefficients of variation (CV) across studied time points indicate good sample preservation for up to 6 months. Cytokine detection stability was confirmed by low CVs, with some differences between fresh and fixed conditions. Thirteen months regular follow-up of PROT samples showed remarkable sample stability. Whole blood can be preserved for phenotyping and cytokine-response studies provided the careful selection of a compatible antibody panel. However, possible changes in cell morphology, differences in antibody affinity, and changes in cytokine-positive cell frequencies when compared to fresh blood should be considered. Our setting constitutes a valuable tool for multicentric and retrospective studies. © 2020 International Society for Advancement of Cytometry.

Precision medicine in autoimmune diseases: fact or fiction.

Much is said about precision medicine, but its real significance and potential are far from certain. Several studies in each of the autoimmune diseases have provided important insights into molecular pathways, but the use of molecular studies, particularly those looking into transcriptome pathways, has seldom approached the possibility of using the data for disease stratification and then for prediction, or for diagnosis. Only the type I IFN signature has been considered for therapeutic purposes, particularly in the case of SLE. This review provides an update on precision medicine, on what can be translated into clinical practice and on what single-cell molecular studies contribute to our knowledge of autoimmune diseases, focusing on a few examples. The main message is that we should try to move from precision medicine of established diseases to preventive medicine in order to predict the development of disease.

B Cells and Microbiota in Autoimmunity.

Trillions of microorganisms inhabit the mucosal membranes maintaining a symbiotic relationship with the host's immune system. B cells are key players in this relationship because activated and differentiated B cells produce secretory immunoglobulin A (sIgA), which binds commensals to preserve a healthy microbial ecosystem. Mounting evidence shows that changes in the function and composition of the gut microbiota are associated with several autoimmune diseases suggesting that an imbalanced or dysbiotic microbiota contributes to autoimmune inflammation. Bacteria within the gut mucosa may modulate autoimmune inflammation through different mechanisms from commensals ability to induce B-cell clones that cross-react with host antigens or through regulation of B-cell subsets' capacity to produce cytokines. Commensal signals in the gut instigate the differentiation of IL-10 producing B cells and IL-10 producing IgA+ plasma cells that recirculate and exert regulatory functions. While the origin of the dysbiosis in autoimmunity is unclear, compelling evidence shows that specific species have a remarkable influence in shaping the inflammatory immune response. Further insight is necessary to dissect the complex interaction between microorganisms, genes, and the immune system. In this review, we will discuss the bidirectional interaction between commensals and B-cell responses in the context of autoimmune inflammation.

Metabolic Disturbances in Urinary and Plasma Samples from Seven Different Systemic Autoimmune Diseases Detected by HPLC-ESI-QTOF-MS.

Systemic autoimmune diseases (SADs) are characterized by dysfunctioning of the immune system, which causes damage in several tissues and organs. Among these pathologies are systemic lupus erythematosus (SLE), systemic sclerosis or scleroderma, Sjögren's syndrome, rheumatoid arthritis, primary antiphospholipid syndrome (PAPS), mixed connective tissue disease (MCTD), and undifferentiated connective tissue disease (UCTD). Early diagnosis is difficult due to similarity in symptoms, signs, and clinical test results. Hence, our aim was to search for differentiating metabolites of these diseases in plasma and urine samples. We performed metabolomic profiling by liquid chromatography-mass spectrometry (LC-MS) of samples from 228 SADs patients and 55 healthy volunteers. Multivariate PLS models were applied to investigate classification accuracies and identify metabolites differentiating SADs and healthy controls. Furthermore, we specifically investigated UCTD against the other SADs. PLS models were able to classify most SADs vs healthy controls (area under the roc curve (AUC) > 0.7), with the exception of MCTD and PAPS. Differentiating metabolites consisted predominantly of unsaturated fatty acids, acylglycines, acylcarnitines, and amino acids. In accordance with the difficulties in defining UCTD, the UCTD metabolome did not differentiate well from the other SADs. However, most UCTD cases were classified as SLE, suggesting that metabolomics may provide a tool to reassess UCTD diagnosis into other conditions for more well-informed therapeutic strategies.

Genetic fine mapping of systemic lupus erythematosus MHC associations in Europeans and African Americans.

Genetic variation within the major histocompatibility complex (MHC) contributes substantial risk for systemic lupus erythematosus, but high gene density, extreme polymorphism and extensive linkage disequilibrium (LD) have made fine mapping challenging. To address the problem, we compared two association techniques in two ancestrally diverse populations, African Americans (AAs) and Europeans (EURs). We observed a greater number of Human Leucocyte Antigen (HLA) alleles in AA consistent with the elevated level of recombination in this population. In EUR we observed 50 different A-C-B-DRB1-DQA-DQB multilocus haplotype sequences per hundred individuals; in the AA sample, these multilocus haplotypes were twice as common compared to Europeans. We also observed a strong narrow class II signal in AA as opposed to the long-range LD observed in EUR that includes class I alleles. We performed a Bayesian model choice of the classical HLA alleles and a frequentist analysis that combined both single nucleotide polymorphisms (SNPs) and classical HLA alleles. Both analyses converged on a similar subset of risk HLA alleles: in EUR HLA- B*08:01 + B*18:01 + (DRB1*15:01 frequentist only) + DQA*01:02 + DQB*02:01 + DRB3*02 and in AA HLA-C*17:01 + B*08:01 + DRB1*15:03 + (DQA*01:02 frequentist only) + DQA*02:01 + DQA*05:01+ DQA*05:05 + DQB*03:19 + DQB*02:02. We observed two additional independent SNP associations in both populations: EUR rs146903072 and rs501480; AA rs389883 and rs114118665. The DR2 serotype was best explained by DRB1*15:03 + DQA*01:02 in AA and by DRB1*15:01 + DQA*01:02 in EUR. The DR3 serotype was best explained by DQA*05:01 in AA and by DQB*02:01 in EUR. Despite some differences in underlying HLA allele risk models in EUR and AA, SNP signals across the extended MHC showed remarkable similarity and significant concordance in direction of effect for risk-associated variants.

ImaGEO: integrative gene expression meta-analysis from GEO database.

SUMMARY: The Gene Expression Omnibus (GEO) database provides an invaluable resource of publicly available gene expression data that can be integrated and analyzed to derive new hypothesis and knowledge. In this context, gene expression meta-analysis (geMAs) is increasingly used in several fields to improve study reproducibility and discovering robust biomarkers. Nevertheless, integrating data is not straightforward without bioinformatics expertise. Here, we present ImaGEO, a web tool for geMAs that implements a complete and comprehensive meta-analysis workflow starting from GEO dataset identifiers. The application integrates GEO datasets, applies different meta-analysis techniques and provides functional analysis results in an easy-to-use environment. ImaGEO is a powerful and useful resource that allows researchers to integrate and perform meta-analysis of GEO datasets to lead robust findings for biomarker discovery studies. AVAILABILITY AND IMPLEMENTATION: ImaGEO is accessible at http://bioinfo.genyo.es/imageo/. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

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.