Single-cell transcriptomic analysis of hematopoietic progenitor cells from patients with systemic lupus erythematosus reveals interferon-inducible reprogramming in early progenitors.

Introduction: Immune cells that contribute to the pathogenesis of systemic lupus erythematosus (SLE) derive from adult hematopoietic stem and progenitor cells (HSPCs) within the bone marrow (BM). For this reason, we reasoned that fundamental abnormalities in SLE can be traced to a BM-derived HSPC inflammatory signature. Methods: BM samples from four SLE patients, six healthy controls, and two umbilical cord blood (CB) samples were used. CD34+ cells were isolated from BM and CB samples, and single-cell RNA-sequencing was performed. Results: A total of 426 cells and 24,473 genes were used in the analysis. Clustering analysis resulted in seven distinct clusters of cell types. Mutually exclusive markers, which were characteristic of each cell type, were identified. We identified three HSPC subpopulations, one of which consisted of proliferating cells (MKI67 expressing cells), one T-like, one B-like, and two myeloid-like progenitor subpopulations. Differential expression analysis revealed i) cell cycle-associated signatures, in healthy BM of HSPC clusters 3 and 4 when compared with CB, and ii) interferon (IFN) signatures in SLE BM of HSPC clusters 3 and 4 and myeloid-like progenitor cluster 5 when compared with healthy controls. The IFN signature in SLE appeared to be deregulated following TF regulatory network analysis and differential alternative splicing analysis between SLE and healthy controls in HSPC subpopulations. Discussion: This study revealed both quantitative-as evidenced by decreased numbers of non-proliferating early progenitors-and qualitative differences-characterized by an IFN signature in SLE, which is known to drive loss of function and depletion of HSPCs. Chronic IFN exposure affects early hematopoietic progenitors in SLE, which may account for the immune aberrancies and the cytopenias in SLE.

Disentangling the riddle of systemic lupus erythematosus with antiphospholipid syndrome: blood transcriptome analysis reveals a less-pronounced IFN-signature and distinct molecular profiles in venous versus arterial events.

INTRODUCTION: Systemic lupus erythematosus with antiphospholipid syndrome (SLE-APS) represents a challenging SLE endotype whose molecular basis remains unknown. METHODS: We analysed whole-blood RNA-sequencing data from 299 patients with SLE (108 SLE-antiphospholipid antibodies (aPL)-positive, including 67 SLE-APS; 191 SLE-aPL-negative) and 72 matched healthy controls (HC). Pathway enrichment analysis, unsupervised weighted gene coexpression network analysis and machine learning were applied to distinguish disease endotypes. RESULTS: Patients with SLE-APS demonstrated upregulated type I and II interferon (IFN) pathways compared with HC. Using a 100-gene random forests model, we achieved a cross-validated accuracy of 75.6% in distinguishing these two states. Additionally, the comparison between SLE-APS and SLE-aPL-negative revealed 227 differentially expressed genes, indicating downregulation of IFN-α and IFN-γ signatures, coupled with dysregulation of the complement cascade, B-cell activation and neutrophil degranulation. Unsupervised analysis of SLE transcriptome identified 21 gene modules, with SLE-APS strongly linked to upregulation of the 'neutrophilic/myeloid' module. Within SLE-APS, venous thromboses positively correlated with 'neutrophilic/myeloid' and 'B cell' modules, while arterial thromboses were associated with dysregulation of 'DNA damage response (DDR)' and 'metabolism' modules. Anticardiolipin and anti-ß2GPI positivity-irrespective of APS status-were associated with the 'neutrophilic/myeloid' and 'protein-binding' module, respectively. CONCLUSIONS: There is a hierarchical upregulation and-likely-dependence on IFN in SLE with the highest IFN signature observed in SLE-aPL-negative patients. Venous thrombotic events are associated with neutrophils and B cells while arterial events with DDR and impaired metabolism. This may account for their differential requirements for anticoagulation and provide rationale for the potential use of mTOR inhibitors such as sirolimus and the direct fIIa inhibitor dabigatran in SLE-APS.

Transcriptomic and proteomic profiling reveals distinct pathogenic features of peripheral non-classical monocytes in systemic lupus erythematosus.

Peripheral blood monocytes propagate inflammation in systemic lupus erythematosus (SLE). Three major populations of monocytes have been recognized namely classical (CM), intermediate (IM) and non-classical monocytes (NCM). Herein, we performed a comprehensive transcriptomic, proteomic and functional characterization of the three peripheral monocytic subsets from active SLE patients and healthy individuals. Our data demonstrate extensive molecular disruptions in circulating SLE NCM, characterized by enhanced inflammatory features such as deregulated DNA repair, cell cycle and heightened IFN signaling combined with differentiation and developmental cues. Enhanced DNA damage, elevated expression of p53, G0 arrest of cell cycle and increased autophagy stress the differentiation potential of NCM in SLE. This immunogenic profile is associated with an activated macrophage phenotype of NCM exhibiting M1 characteristics in the circulation, fueling the inflammatory response. Together, these findings identify circulating SLE NCM as a pathogenic cell type in the disease that could represent an additional therapeutic target.

A network-based approach reveals long non-coding RNAs associated with disease activity in lupus nephritis: key pathways for flare and potential biomarkers to be used as liquid biopsies.

Objective: A blood-based biomarker is needed to assess lupus nephritis (LN) disease activity, minimizing the need for invasive kidney biopsies. Long non-coding RNAs (lncRNAs) are known to regulate gene expression, appear to be stable in human plasma, and can serve as non-invasive biomarkers. Methods: Transcriptomic data of whole blood samples from 74 LN patients and 20 healthy subjects (HC) were analyzed to identify differentially expressed (DE) lncRNAs associated with quiescent disease and flares. Weighted gene co-expression network analysis (WGCNA) was performed to uncover lncRNAs with a central role (hub lncRNAs) in regulating key biological processes that drive LN disease activity. The association of hub lncRNAs with disease activity was validated using RT-qPCR on an independent cohort of 15 LN patients and 9 HC. cis- and trans-targets of validated lncRNAs were explored in silico to examine potential mechanisms of their action. Results: There were 444 DE lncRNAs associated with quiescent disease and 6 DE lncRNAs associated with flares (FDR <0.05). WGCNA highlighted IFN signaling and B-cell activity/adaptive immunity as the most significant processes contributing to nephritis activity. Four disease-activity-associated lncRNAs, namely, NRIR, KLHDC7B-DT, MIR600HG, and FAM30A, were detected as hub genes and validated in an independent cohort. NRIR and KLHDC7B-DT emerged as potential key regulators of IFN-mediated processes. Network analysis suggests that FAM30A and MIR600HG are likely to play a central role in the regulation of B-cells in LN through cis-regulation effects and a competing endogenous RNA mechanism affecting immunoglobulin gene expression and the IFN-λ pathway. Conclusions: The expression of lncRNAs NRIR, KLHDC7B-DT, FAM30A, and MIR600HG were associated with disease activity and could be further explored as blood-based biomarkers and potential liquid biopsy on LN.

Restoration of aberrant gene expression of monocytes in systemic lupus erythematosus via a combined transcriptome-reversal and network-based drug repurposing strategy.

BACKGROUND: Monocytes -key regulators of the innate immune response- are actively involved in the pathogenesis of systemic lupus erythematosus (SLE). We sought to identify novel compounds that might serve as monocyte-directed targeted therapies in SLE. RESULTS: We performed mRNA sequencing in monocytes from 15 patients with active SLE and 10 healthy individuals. Disease activity was assessed with the Systemic Lupus Erythematosus Disease Activity Index 2000 (SLEDAI-2 K). Leveraging the drug repurposing platforms iLINCS, CLUE and L1000CDS2, we identified perturbagens capable of reversing the SLE monocyte signature. We identified transcription factors and microRNAs (miRNAs) that regulate the transcriptome of SLE monocytes, using the TRRUST and miRWalk databases, respectively. A gene regulatory network, integrating implicated transcription factors and miRNAs was constructed, and drugs targeting central components of the network were retrieved from the DGIDb database. Inhibitors of the NF-κB pathway, compounds targeting the heat shock protein 90 (HSP90), as well as a small molecule disrupting the Pim-1/NFATc1/NLRP3 signaling axis were predicted to efficiently counteract the aberrant monocyte gene signature in SLE. An additional analysis was conducted, to enhance the specificity of our drug repurposing approach on monocytes, using the iLINCS, CLUE and L1000CDS2 platforms on publicly available datasets from circulating B-lymphocytes, CD4+ and CD8+ T-cells, derived from SLE patients. Through this approach we identified, small molecule compounds, that could potentially affect more selectively the transcriptome of SLE monocytes, such as, certain NF-κB pathway inhibitors, Pim-1 and SYK kinase inhibitors. Furthermore, according to our network-based drug repurposing approach, an IL-12/23 inhibitor and an EGFR inhibitor may represent potential drug candidates in SLE. CONCLUSIONS: Application of two independent - a transcriptome-reversal and a network-based -drug repurposing strategies uncovered novel agents that might remedy transcriptional disturbances of monocytes in SLE.

Combined - whole blood and skin fibroblasts- transcriptomic analysis in Psoriatic Arthritis reveals molecular signatures of activity, resistance and early response to treatment.

Background: An interplay between immune cells and resident skin and joint stromal cells is implicated in psoriatic arthritis (PsA), yet the mechanisms remain elusive with a paucity of molecular biomarkers for activity and response. Combined transcriptomic and immunophenotypic analysis of whole blood and skin fibroblasts could provide further insights. Methods: Whole blood RNA-seq was performed longitudinally in 30 subjects with PsA at the beginning, one and six months after treatment, with response defined at six months. As control groups, 10 healthy individuals and 10 subjects with rheumatoid arthritis (RA) were recruited combined with public datasets from patients with psoriasis (PsO) and systemic lupus erythematous (SLE). Differential expression analysis and weighted gene co-expression network analysis were performed to identify gene expression signatures, while deconvolution and flow cytometry to characterize the peripheral blood immune cell profile. In a subset of affected and healthy individuals, RNA-seq of skin fibroblasts was performed and subjected to CellChat analysis to identify the blood-skin fibroblast interaction network. Results: PsA demonstrated a distinct "activity" gene signature in the peripheral blood dominated by TNF- and IFN-driven inflammation, deregulated cholesterol and fatty acid metabolism and expansion of pro-inflammatory non-classical monocytes. Comparison with the blood transcriptome of RA, PsO, and SLE revealed a "PsA-specific signature" enriched in extracellular matrix remodeling. This was further supported by the skin fibroblast gene expression profile, displaying an activated, proliferating phenotype, and by skin-blood interactome analysis revealing interactions with circulating immune cells through WNT, PDGF and immune-related semaphorins. Of note, resistance to treatment was associated with upregulation of genes involved in TGFß signaling and angiogenesis and persistent increase of non-classical monocytes. Differentially expressed genes related to platelet activation and hippo signaling discriminated responders and non-responders as early as one month after treatment initiation. Conclusion: Transcriptome analysis of peripheral blood and skin fibroblasts in PsA reveals a distinct disease activity signature and supports the involvement of skin fibroblasts through their activation and interaction with circulating immune cells. Aberrant TGFß signaling and persistently increased non-classical monocytes characterize treatment-resistant PsA, with pro-inflammatory pathways related to platelet activation and Hippo signaling predicting early response to treatment.

Molecular Taxonomy of Systemic Lupus Erythematosus Through Data-Driven Patient Stratification: Molecular Endotypes and Cluster-Tailored Drugs.

Objectives: Treatment of Systemic Lupus Erythematosus (SLE) is characterized by a largely empirical approach and relative paucity of novel compound development. We sought to stratify SLE patients based on their molecular phenotype and identify putative therapeutic compounds for each molecular fingerprint. Methods: By the use of whole blood RNA-seq data from 120 SLE patients, and in a data-driven, clinically unbiased manner, we established modules of commonly regulated genes (molecular endotypes) and re-stratified patients through hierarchical clustering. Disease activity and severity were assessed using SLEDAI-2K and Lupus Severity Index, respectively. Through an in silico drug prediction pipeline, we investigated drugs currently in use, tested in lupus clinical trials, and listed in the iLINCS prediction databases, for their ability to reverse the gene expression signatures in each molecular endotype. Drug repurposing analysis was also performed to identify perturbagens that counteract group-specific SLE signatures. Results: Molecular taxonomy identified five lupus endotypes, each characterized by a unique gene module enrichment pattern. Neutrophilic signature group consisted primarily of patients with active lupus nephritis, while the B-cell expression group included patients with constitutional features. Patients with moderate severity and serologic activity exhibited a signature enriched for metabolic processes. Mild disease was distributed in two groups, exhibiting enhanced basic cellular functions, myelopoiesis, and autophagy. Bortezomib was predicted to reverse disturbances in the "neutrophilic" cluster, azathioprine and ixazomib in the "B-cell" cluster, and fostamatinib in the "metabolic" patient subgroup. Conclusion: The clinical spectrum of SLE encompasses distinct molecular endotypes, each defined by unique pathophysiologic aberrancies potentially reversible by distinct compounds.

Τhe Complete Mitochondrial Genome of Bactrocera carambolae (Diptera: Tephritidae): Genome Description and Phylogenetic Implications.

Bactrocera carambolae is one of the approximately 100 sibling species of the Bactrocera dorsalis complex and considered to be very closely related to B. dorsalis. Due to their high morphological similarity and overlapping distribution, as well as to their economic impact and quarantine status, the development of reliable markers for species delimitation between the two taxa is of great importance. Here we present the complete mitochondrial genome of B. carambolae sourced from its native range in Malaysia and its invaded territory in Suriname. The mitogenome of B. carambolae presents the typical organization of an insect mitochondrion. Comparisons of the analyzed B. carambolae sequences to all available complete mitochondrial sequences of B. dorsalis revealed several species-specific polymorphic sites. Phylogenetic analysis based on Bactrocera mitogenomes supports that B. carambolae is a differentiated taxon though closely related to B. dorsalis. The present complete mitochondrial sequences of B. carambolae could be used, in the frame of Integrative Taxonomy, for species discrimination and resolution of the phylogenetic relationships within this taxonomically challenging complex, which would facilitate the application of species-specific population suppression strategies, such as the sterile insect technique.