Limited cutaneous systemic sclerosis skin demonstrates distinct molecular subsets separated by a cardiovascular development gene expression signature

Differentially expressed genes in systemic sclerosis: Towards predictive medicine with new molecular tools for clinicians

Systemic sclerosis (SSc) is a rare and chronic autoimmune disease characterized by a pathogenic triad of immune dysregulation, vasculopathy, and progressive fibrosis. Clinical tools commonly used to assess patients, such as the modified Rodnan skin score, difference between limited or diffuse forms of skin involvement, presence of lung, heart or kidney involvement, or of various autoantibodies, are important prognostic factors, but still fail to reflect the large heterogeneity of the disease. SSc treatment options are diverse, ranging from conventional drugs to autologous hematopoietic stem cell transplantation, and predicting response is challenging. Genome-wide technologies, such as high throughput microarray analyses and RNA sequencing, allow accurate, unbiased, and broad assessment of alterations in expression levels of multiple genes. In recent years, many studies have shown robust changes in the gene expression profiles of SSc patients compared to healthy controls, mainly in skin tissues and peripheral blood cells. The objective analysis of molecular patterns in SSc is a powerful tool that can further classify SSc patients with similar clinical phenotypes and help predict response to therapy. In this review, we describe the journey from the first discovery of differentially expressed genes to the identification of enriched pathways and intrinsic subsets identified in SSc, using machine learning algorithms. Finally, we discuss the use of these new tools to predict the efficacy of various treatments, including stem cell transplantation. We suggest that the use of RNA gene expression-based classifications according to molecular subsets may bring us one step closer to precision medicine in Systemic Sclerosis.

Alterations of the Primary Cilia Gene SPAG17 and SOX9 Locus Noncoding RNAs Identified by RNA-Sequencing Analysis in Patients With Systemic Sclerosis

OBJECTIVE

Systemic sclerosis (SSc) is characterized by immune activation, vasculopathy, and unresolving fibrosis in the skin, lungs, and other organs. We performed RNA-sequencing analysis on skin biopsy samples and peripheral blood mononuclear cells (PBMCs) from SSc patients and unaffected controls to better understand the pathogenesis of SSc.

METHODS

We analyzed these data 1) to test for case/control differences and 2) to identify genes whose expression levels correlate with SSc severity as measured by local skin score, modified Rodnan skin thickness score (MRSS), forced vital capacity (FVC), or diffusing capacity for carbon monoxide (DLco).

RESULTS

We found that PBMCs from SSc patients showed a strong type I interferon signature. This signal was found to be replicated in the skin, with additional signals for increased extracellular matrix (ECM) genes, classical complement pathway activation, and the presence of B cells. Notably, we observed a marked decrease in the expression of SPAG17, a cilia component, in SSc skin. We identified genes that correlated with the MRSS, DLco, and FVC in SSc PBMCs and skin using weighted gene coexpression network analysis. These genes were largely distinct from the case/control differentially expressed genes. In PBMCs, type I interferon signatures negatively correlated with the DLco. In SSc skin, ECM gene expression positively correlated with the MRSS. Network analysis of SSc skin genes that correlated with clinical features identified the noncoding RNAs SOX9-AS1 and ROCR, both near the SOX9 locus, as highly connected, "hub-like" genes in the network.

CONCLUSION

These results identify noncoding RNAs and SPAG17 as novel factors potentially implicated in the pathogenesis of SSc.

Current Concepts on the Pathogenesis of Systemic Sclerosis

From the clinical standpoint, systemic sclerosis (SSc) is characterized by skin and internal organ fibrosis, diffuse fibroproliferative vascular modifications, and autoimmunity. Clinical presentation and course are highly heterogenous and life expectancy variably affected mostly dependent on lung and heart involvement. SSc touches more women than men with differences in disease severity and environmental exposure. Pathogenetic events originate from altered homeostasis favored by genetic predisposition, environmental cues and a variety of endogenous and exogenous triggers. Epigenetic modifications modulate SSc pathogenesis which strikingly associate profound immune-inflammatory dysregulation, abnormal endothelial cell behavior, and cell trans-differentiation into myofibroblasts. SSc myofibroblasts show enhanced survival and enhanced extracellular matrix deposition presenting altered structure and altered physicochemical properties. Additional cell types of likely pathogenic importance are pericytes, platelets, and keratinocytes in conjunction with their relationship with vessel wall cells and fibroblasts. In SSc, the profibrotic milieu is favored by cell signaling initiated in the one hand by transforming growth factor-beta and related cytokines and in the other hand by innate and adaptive type 2 immune responses. Radical oxygen species and invariant receptors sensing danger participate to altered cell behavior. Conventional and SSc-specific T cell subsets modulate both fibroblasts as well as endothelial cell dysfunction. Beside autoantibodies directed against ubiquitous antigens important for enhanced clinical classification, antigen-specific agonistic autoantibodies may have a pathogenic role. Recent studies based on single-cell RNAseq and multi-omics approaches are revealing unforeseen heterogeneity in SSc cell differentiation and functional states. Advances in system biology applied to the wealth of data generated by unbiased screening are allowing to subgroup patients based on distinct pathogenic mechanisms. Deciphering heterogeneity in pathogenic mechanisms will pave the way to highly needed personalized therapeutic approaches.

Differential DNA Methylation Landscape in Skin Fibroblasts from African Americans with Systemic Sclerosis

The etiology and reasons underlying the ethnic disparities in systemic sclerosis (SSc) remain unknown. African Americans are disproportionally affected by SSc and yet are underrepresented in research. The aim of this study was to comprehensively investigate the association of DNA methylation levels with SSc in dermal fibroblasts from patients of African ancestry. Reduced representation bisulfite sequencing (RRBS) was performed on primary dermal fibroblasts from 15 SSc patients and 15 controls of African ancestry, and over 3.8 million CpG sites were tested for differential methylation patterns between cases and controls. The dermal fibroblasts from African American patients exhibited widespread reduced DNA methylation. Differentially methylated CpG sites were most enriched in introns and intergenic regions while depleted in 5' UTR, promoters, and CpG islands. Seventeen genes and eleven promoters showed significant differential methylation, mostly in non-coding RNA genes and pseudogenes. Gene set enrichment analysis (GSEA) and gene ontology (GO) analyses revealed an enrichment of pathways related to interferon signaling and mesenchymal differentiation. The hypomethylation of DLX5 and TMEM140 was accompanied by these genes' overexpression in patients but underexpression for lncRNA MGC12916. These data show that differential methylation occurs in dermal fibroblasts from African American patients with SSc and identifies novel coding and non-coding genes.

Molecular "omic" signatures in systemic sclerosis

Systemic sclerosis (SSc) is a connective tissue disorder characterized by immunologic, vascular, and extracellular matrix abnormalities. Variation in the proportion and/or timing of activation in the deregulated molecular pathways that underlie SSc may explain the observed clinical heterogeneity in terms of disease phenotype and treatment response. In recent years, SSc research has generated massive amounts of "omics" level data. In this review, we discuss the body of "omics" level work in SSc and how each layer provides unique insight to our understanding of SSc. We posit that effective integration of genomic, transcriptomic, metagenomic, and epigenomic data is an important step toward precision medicine and is vital to the identification of effective therapeutic options for patients with SSc.

Screening and identification of biomarkers for systemic sclerosis via microarray technology

Systemic sclerosis (SSc) is a complex autoimmune disease. The pathogenesis of SSc is currently unclear, although like other rheumatic diseases its pathogenesis is complicated. However, the ongoing development of bioinformatics technology has enabled new approaches to research this disease using microarray technology to screen and identify differentially expressed genes (DEGs) in the skin of patients with SSc compared with individuals with healthy skin. Publicly available data were downloaded from the Gene Expression Omnibus (GEO) database and intra-group data repeatability tests were conducted using Pearson's correlation test and principal component analysis. DEGs were identified using an online tool, GEO2R. Functional annotation of DEGs was performed using Gene Ontology (GO) and the Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis. Finally, the construction and analysis of the protein-protein interaction (PPI) network and identification and analysis of hub genes was carried out. A total of 106 DEGs were detected by the screening of SSc and healthy skin samples. A total of 10 genes [interleukin-6, bone morphogenetic protein 4, calumenin (CALU), clusterin, cysteine rich angiogenic inducer 61, serine protease 23, secretogranin II, suppressor of cytokine signaling 3, Toll-like receptor 4 (TLR4), tenascin C] were identified as hub genes with degrees ≥10, and which could sensitively and specifically predict SSc based on receiver operator characteristic curve analysis. GO and KEGG analysis showed that variations in hub genes were mainly enriched in positive regulation of nitric oxide biosynthetic processes, negative regulation of apoptotic processes, extracellular regions, extracellular spaces, cytokine activity, chemo-attractant activity, and the phosphoinositide 3 kinase-protein kinase B signaling pathway. In summary, bioinformatics techniques proved useful for the screening and identification of biomarkers of disease. A total of 106 DEGs and 10 hub genes were linked to SSc, in particular the TLR4 and CALU genes.

Existing and novel biomarkers for precision medicine in systemic sclerosis

The discovery and validation of biomarkers resulting from technological advances in the analysis of genomic, transcriptomic, lipidomic and metabolomic pathways involved in the pathogenesis of complex human diseases have led to the development of personalized and rationally designed approaches for the clinical management of such disorders. Although some of these approaches have been applied to systemic sclerosis (SSc), an unmet need remains for validated, non-invasive biomarkers to aid in the diagnosis of SSc, as well as in the assessment of disease progression and response to therapeutic interventions. Advances in global transcriptomic technology over the past 15 years have enabled the assessment of microRNAs that circulate in the blood of patients and the analysis of the macromolecular content of a diverse group of lipid bilayer membrane-enclosed extracellular vesicles, such as exosomes and other microvesicles, which are released by all cells into the extracellular space and circulation. Such advances have provided new opportunities for the discovery of biomarkers in SSc that could potentially be used to improve the design and evaluation of clinical trials and that will undoubtedly enable the development of personalized and individualized medicine for patients with SSc.

Identification of Key Genes and Pathways in Female Lung Cancer Patients Who Never Smoked by a Bioinformatics Analysis

Smoking is considered the major risk factor for lung cancer, but only a small portion of female lung adenocarcinoma patients are associated with smoking. Thus, identifying crucial genes and pathways related to nonsmoking female lung cancer patients is of great importance. Gene expression profiles were downloaded from the Gene Expression Omnibus (GEO) and the Cancer Genome Atlas (TCGA) databases. The R software packages were applied to screen the differentially expressed genes (DEGs). GO term enrichment and KEGG pathway analyses were carried out using DAVID tools. The protein-protein interaction (PPI) network was constructed by Cytoscape software. In total, 487 downregulated and 199 upregulated DEGs were identified. The down-regulated DEGs were mainly enriched for behavior and response to wounding, and the upregulated DEGs were significantly enriched for multicellular organismal metabolic process and cell division. The KEGG pathway analysis revealed that the downregulated DEGs were significantly enriched for cell adhesion molecules and neuroactive ligand-receptor interaction, while the upregulated DEGs were mainly enriched for cell cycle and the p53 signaling pathway. The top 10 hub genes and top 3 gene interaction modules were selected from the PPI network. Of the ten hub genes, a high expression of five genes was related to a poor OS in female lung cancer patients who never smoked, including IL6, CXCR2, FPR2, PPBP and HBA1. However, a low expression of GNG11, LRRK2, CDH5, CAV1 and SELE was associated with a worse OS for the female lung cancer patients who never smoked. In conclusion, our study provides novel insight for a better understanding of the pathogenesis of nonsmoking female lung cancer, and these identified DEGs may serve as biomarkers for diagnostics and treatment.

Integrated analysis of hepatic mRNA and miRNA profiles identified molecular networks and potential biomarkers of NAFLD

To enhance our understanding of molecular mechanisms and mine novel biomarkers of non-alcoholic fatty liver disease (NAFLD), RNA sequencing was performed to gain hepatic expression profiles of mRNAs and miRNAs in NAFLD and normal rats. Using DESeq with thresholds of a two-fold change and a false discovery rate (FDR) less than 0.05, 336 mRNAs and 21 miRNAs were identified as differentially expressed. Among those, 17 miRNAs (e.g., miR-144-3p, miR-99a-3p, miR-200b-3p, miR-200b-5p, miR-200c-3p, etc.) might serve as novel biomarkers of NAFLD. MiRNA target genes (13565) were predicted by the miRWalk database. Using DAVID 6.8, the intersection (195 genes) of differentially expressed mRNAs and miRNA-predicted target genes were enriched in 47 gene ontology (GO) terms and 28 Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways. Using Cytoscape, pathway interaction and protein-protein interaction (PPI) networks were constructed, and hub genes (e.g., Abcg8, Cyp1a1, Cyp51, Hmgcr, etc.) associated with NAFLD were obtained. Moreover, 673 miRNA-mRNA negative regulatory pairs were obtained, and networks were constructed. Finally, several representative miRNAs and mRNAs were validated by real-time qPCR. In conclusion, potential molecular mechanisms of NAFLD could be inferred from integrated analysis of mRNA and miRNA profiles, which may indicate novel biomarkers of NAFLD.