Long Non-Coding RNA MAGI2-AS3 is a New Player with a Tumor Suppressive Role in High Grade Serous Ovarian Carcinoma

High-Grade Serous Ovarian Carcinoma (HGSC) is the most incidental and lethal subtype of epithelial ovarian cancer (EOC) with a high mortality rate of nearly 65%. Recent findings aimed at understanding the pathogenesis of HGSC have attributed its principal source as the Fallopian Tube (FT). To further comprehend the exact mechanism of carcinogenesis, which is still less known, we performed a transcriptome analysis comparing FT and HGSC. Our study aims at exploring new players involved in the development of HGSC from FT, along with their signaling network, and we chose to focus on non-coding RNAs. Non-coding RNAs (ncRNAs) are increasingly observed to be the major regulators of several cellular processes and could have key functions as biological markers, as well as even a therapeutic approach. The most physiologically relevant and significantly dysregulated non-coding RNAs were identified bioinformatically. After analyzing the trend in HGSC and other cancers, MAGI2-AS3 was observed to be an important player in EOC. We assessed its tumor-suppressive role in EOC by means of various assays. Further, we mapped its signaling pathway using its role as a miRNA sponge to predict the miRNAs binding to MAGI2AS3 and showed it experimentally. We conclude that MAGI2-AS3 acts as a tumor suppressor in EOC, specifically in HGSC by sponging miR-15-5p, miR-374a-5p and miR-374b-5p, and altering downstream signaling of certain mRNAs through a ceRNA network.

Peptides from Extruded Lupin (Lupinus albus L.) Regulate Inflammatory Activity via the p38 MAPK Signal Transduction Pathway in RAW 264.7 Cells

In this study, protein was extracted from extruded lupin and submitted to gastroduodenal digests to obtain lupin peptides, which were characterized using ultraperformance liquid chromatography-tandem mass spectrometry (UPLC-MS/MS). After this, IQDKEGIPPDQQR (IQD), the lupine peptide monomer characterized after UPLC-MS/MS, was screened out by macrophage inflammatory cytokine production assay. RNA-sequencing analysis was performed to explore the mechanisms underlying the anti-inflammatory activity associated with this peptide. The results indicated that lupin peptides effectively inhibited the lipopolysaccharide-induced overproduction of proinflammatory mediators. IQD inhibited the production of tumor necrosis factor-α, interleukin (IL)-6, IL-1β, and monocyte chemoattractant protein-1 by 51.20, 38.52, 44.70, and 40.43%, respectively. RNA-sequencing results showed that IQD inhibited the inflammatory response by regulating the gene expression of the p38 mitogen-activated protein kinase pathway and inhibiting downstream inflammatory cytokines. These bioactive peptides may be used to develop new ingredients for anti-inflammatory nutritional supplements.

Identification of Galectin-3 as Potential Biomarkers for Renal Fibrosis by RNA-Sequencing and Clinicopathologic Findings of Kidney Biopsy

Background: Galectin-3 (Gal-3) is a multifunctional glycan-binding protein shown to be linked to chronic inflammation and fibrogenesis. Plasma Gal-3 is associated with proteinuria and renal dysfunction, but its role has never been confirmed with kidney biopsy results. In our study, we aimed to explore the expression of Gal-3 in biopsy-proven patients, and we tested the hypothesis that chronic kidney disease (CKD) leads to upregulation of plasma Gal-3 expression in corresponding biopsy findings and RNA sequencing analysis. Method: In 249 patients (male/female: 155/94, age: 57.2 ± 16.3 years) who underwent kidney biopsy, plasma levels of Gal-3 were measured to estimate the association of renal fibrosis. Relationships between plasma Gal-3 levels, estimated glomerular filtration rate (eGFR) and renal histology findings were also assessed. We further examined the gene expression of Gal-3 in RNA-sequencing analysis in biopsy-proven patients. Results: Compared to patients without CKD, CKD patients had higher levels of plasma Gal-3 (1,016.3 ± 628.1 pg/mL vs. 811.6 ± 369.6 pg/ml; P = 0.010). Plasma Gal-3 was inversely correlated with eGFR (P = 0.005) but not with proteinuria. Higher Gal-3 levels were associated with interstitial fibrosis, tubular atrophy and vascular intimal fibrosis. RNA-sequencing analysis showed the upregulation of Gal-3 in fibrotic kidney biopsy samples, and the differentially expressed genes were mainly enhanced in immune cell activation and the regulation of cell-cell adhesion. Conclusions: Plasma Gal-3 levels are inverse correlated with eGFR but positively correlated with renal fibrosis, which may be involved in the immune response and associated pathways. These findings support the role of Gal-3 as a predictive marker of renal fibrosis.

Mapping Whole-Transcriptome Splicing in Mouse Hematopoietic Stem Cells

Hematopoietic stem cells (HSCs) are rare cells that generate all the various types of blood and immune cells. High-quality transcriptome data have enabled the identification of significant genes for HSCs. However, most genes are expressed in various forms by alternative splicing (AS), extending transcriptome complexity. Here, we delineate AS to determine which isoforms are expressed in mouse HSCs. Our analysis of microarray and RNA-sequencing data includes differential expression of splicing factors that may regulate AS, and a complete map of splicing isoforms. Multiple types of isoforms for known HSC genes and unannotated splicing that may alter gene function are presented. Transcriptome-wide identification of genes and their respective isoforms in mouse HSCs will open another dimension for adult stem cells.

Reprogramming of Cancer Cells into Induced Pluripotent Stem Cells Questioned

Background and Objectives

Several recent studies have claimed that cancer cells can be reprogrammed into induced pluripotent stem cells (iPSCs). However, in most cases, cancer cells seem to be resistant to cellular reprogramming. Furthermore, the underlying mechanisms of limited reprogramming in cancer cells are largely unknown. Here, we identified the candidate barrier genes and their target genes at the early stage of reprogramming for investigating cancer reprogramming.

Methods

We tried induction of pluripotency in normal human fibroblasts (BJ) and both human benign (MCF10A) and malignant (MCF7) breast cancer cell lines using a classical retroviral reprogramming method. We conducted RNA-sequencing analysis to compare the transcriptome of the three cell lines at early stage of reprogramming.

Results

We could generate iPSCs from BJ, whereas we were unable to obtain iPSCs from cancer cell lines. To address the underlying mechanism of limited reprogramming in cancer cells, we identified 29 the candidate barrier genes based on RNA-sequencing data. In addition, we found 40 their target genes using Cytoscape software.

Conclusions

Our data suggest that these genes might one of the roadblock for cancer cell reprogramming. Furthermore, we provide new insights into application of iPSCs technology in cancer cell field for therapeutic purposes.

HbxB Is a Key Regulator for Stress Response and β-Glucan Biogenesis in Aspergillus nidulans

Homeobox transcription factors are conserved in eukaryotes and act as multi-functional transcription factors in filamentous fungi. Previously, it was demonstrated that HbxB governs fungal development and spore viability in Aspergillus nidulans. Here, the role of HbxB in A. nidulans was further characterized. RNA-sequencing revealed that HbxB affects the transcriptomic levels of genes associated with trehalose biosynthesis and response to thermal, oxidative, and radiation stresses in asexual spores called conidia. A phenotypic analysis found that hbxB deletion mutant conidia were more sensitive to ultraviolet stress. The loss of hbxB increased the mRNA expression of genes associated with β-glucan degradation and decreased the amount of β-glucan in conidia. In addition, hbxB deletion affected the expression of the sterigmatocystin gene cluster and the amount of sterigmatocystin. Overall, these results indicated that HbxB is a key transcription factor regulating trehalose biosynthesis, stress tolerance, β-glucan degradation, and sterigmatocystin production in A.nidulans conidia.

Iron redistribution induces oxidative burst and resistance in maize against Curvularia lunata

ΔClnps6 induced iron redistribution in maize B73 leaf cells and resulted in reactive oxygen species (ROS) burst to enhance plant resistance against Curvularia lunata. Iron is an indispensable co-factor of various crucial enzymes that are involved in cellular metabolic processes and energy metabolism in eukaryotes. For this reason, plants and pathogens compete for iron to maintain their iron homeostasis, respectively. In our previous study, ΔClnps6, the extracellular siderophore biosynthesis deletion mutant of Curvularia lunata, was sensitive to exogenous hydrogen peroxide and virulence reduction. However, the mechanism was not studied. Here, we report that maize B73 displayed highly resistance to ΔClnps6. The plants recruited more iron at cell wall appositions (CWAs) to cause ROS bursts. Intracellular iron deficiency induced by iron redistribution originated form up-regulated expression of genes involved in intracellular iron consumption in leaves and absorption in roots. The RNA-sequencing data also showed that the expression of respiratory burst oxidase homologue (ZmRBOH4) and NADP-dependent malic enzyme 4 (ZmNADP-ME4) involved in ROS production was up-regulated in maize B73 after ΔClnps6 infection. Simultaneously, jasmonic acid (JA) biosynthesis genes lipoxygenase (ZmLOX), allene oxide synthase (ZmAOS), GA degradation gene gibberellin 2-beta-dioxygenase (ZmGA2OX6) and ABA degradation genes abscisic acid hydroxylase (ZmABH1, ZmABH2) involved in iron homeostasis were up-regulated expression. Ferritin1 (ZmFER1) positive regulated maize resistance against C. lunata via ROS burst under Fe-limiting conditions. Overall, our results showed that iron played vital roles in activating maize resistance in B73-C. lunata interaction.

Systematic Identification of circRNAs in Alzheimer's Disease

Mammalian circRNAs are covalently closed circular RNAs often generated through backsplicing of precursor linear RNAs. Although their functions are largely unknown, they have been found to influence gene expression at different levels and in a wide range of biological processes. Here, we investigated if some circRNAs may be differentially abundant in Alzheimer’s Disease (AD). We identified and analyzed publicly available RNA-sequencing data from the frontal lobe, temporal cortex, hippocampus, and plasma samples reported from persons with AD and persons who were cognitively normal, focusing on circRNAs shared across these datasets. We identified an overlap of significantly changed circRNAs among AD individuals in the various brain datasets, including circRNAs originating from genes strongly linked to AD pathology such as DOCK1, NTRK2, APC (implicated in synaptic plasticity and neuronal survival) and DGL1/SAP97, TRAPPC9, and KIF1B (implicated in vesicular traffic). We further predicted the presence of circRNA isoforms in AD using specialized statistical analysis packages to create approximations of entire circRNAs. We propose that the catalog of differentially abundant circRNAs can guide future investigation on the expression and splicing of the host transcripts, as well as the possible roles of these circRNAs in AD pathogenesis.

Whole transcriptome analysis for T cell receptor-affinity and IRF4-regulated clonal expansion of T cells

Clonal population expansion of T cells during an immune response is dependent on the affinity of the T cell receptor (TCR) for its antigen [1]. However, there is little understanding of how this process is controlled transcriptionally. We found that the transcription factor IRF4 was induced in a manner dependent on TCR-affinity and was critical for the clonal expansion and maintenance of effector function of antigen-specific CD8⁺ T cells. We performed a genome-wide expression profiling experiment using RNA sequencing technology (RNA-seq) to interrogate global expression changes when IRF4 was deleted in CD8⁺ T cells activated with either a low or high affinity peptide ligand. This allowed us not only to determine IRF4-dependent transcriptional changes but also to identify transcripts dependent on TCR-affinity [2]. Here we describe in detail the analyses of the RNA-seq data, including quality control, read mapping, quantification, normalization and assessment of differential gene expression. The RNA-seq data can be accessed from Gene Expression Omnibus database (accession number GSE49929).

d-allose protects brain microvascular endothelial cells from hypoxic/reoxygenated injury by inhibiting endoplasmic reticulum stress

Ischemic stroke is an acute brain disease with a high mortality rate. Currently, the only effective method is to restore the blood supply. But the inflammation and oxidative stress induced by this approach can damage the integrity of the endothelial system, which hampers the patient's outcome. d-allose has the biological activity to protect against ischemia-reperfusion injury, however, the underlying mechanism remains unclear. Here, brain microvascular endothelial cells (RBMECs) were used as the study material to establish an IR-injury model. Cell viability of RBMECs was suppressed after hypoxia/reoxygenation (H/R) treatment and significantly increased after d-allose supplementation. RNAseq results showed 180 differentially expressed genes (DEGs) between the therapy group (H/R + Dal) and the model group (H/R), of which 151 DEGs were restored to control levels by d-allose. Enrichment analysis revealed that DEGs were mainly involved in protein processing in endoplasmic reticulum. 6 DEGs in the unfolded protein response (UPR) pathway were verified by qRT-PCR. All of them were significantly down-regulated by d-allose, indicating that endoplasmic reticulum stress (ERS) was relieved. In addition, d-allose significantly inhibited the phosphorylation level of eIF2α, a marker of ERS. The downstream molecules of Phosphorylation of eIF2α, Gadd45a and Chac1, which trigger cycle arrest and apoptosis, respectively, were also significantly inhibited by d-allose. Thus, we conclude that d-allose inhibits the UPR pathway, attenuates eIF2α phosphorylation and ERS, restores the cell cycle, inhibits apoptosis, and thus enhances endothelial cell tolerance to H/R injury.

VF-4 and DR-8 Derived from Salted Egg White Inhibit Inflammatory Activity via NF-κB/PI3K-Akt/MAPK Signal Transduction Pathways in HT-29 Cells Induced by TNF-α

SCOPE

Inflammation is the pathological basis of many chronic diseases, and persistent intestinal inflammation is a key factor in the further development of colon cancer. Egg-derived peptides have been proven to have anti-intestinal inflammation activity. Egg white treated with salt contains a lot of rich protein, whether its peptides have anti-inflammatory activity and how their mechanism of action is still unclear.

METHODS AND RESULTS

In this study, ELISA is used to determine the anti-inflammatory activity of the peptides (VF-4 and DR-8 from salted egg white), and then RNA-seq is used to explore the mechanism of their anti-inflammatory activity, and then verified by western blotting and inhibitors. The results show that VF-4 and DR-8 significantly inhibit TNF-α-induced IL-8 secretion in HT-29 cells in a concentration-dependent manner, and VF-4 show a more significant anti-inflammatory effect than DR-8. The anti-inflammatory mechanism of VF-4 and DR-8 is through inhibiting the activation of Nuclear factor kappa B (NF-κB), mitogen-activated protein kinase (MAPK), and phosphatidylinositol 3' -kinase(PI3K)-Akt pathways, reducing the production of inflammatory mediators.

CONCLUSION

VF-4 and DR-8 have obvious anti-inflammatory activity, which can reduce intestinal inflammation and inhibit its further development into colon cancer.

Reduced Function of the Glutathione S-Transferase S1 Suppresses Behavioral Hyperexcitability in Drosophila Expressing Mutant Voltage-Gated Sodium Channels

Voltage-gated sodium (Na_v) channels play a central role in the generation and propagation of action potentials in excitable cells such as neurons and muscles. To determine how the phenotypes of Na_v-channel mutants are affected by other genes, we performed a forward genetic screen for dominant modifiers of the seizure-prone, gain-of-function Drosophila melanogaster Na_v-channel mutant, para^Shu. Our analyses using chromosome deficiencies, gene-specific RNA interference, and single-gene mutants revealed that a null allele of glutathione S-transferase S1 (GstS1) dominantly suppresses para^Shu phenotypes. Reduced GstS1 function also suppressed phenotypes of other seizure-prone Na_v-channel mutants, para^GEFS+ and para^bss. Notably, para^Shu mutants expressed 50% less GstS1 than wild-type flies, further supporting the notion that para^Shu and GstS1 interact functionally. Introduction of a loss-of-function GstS1 mutation into a para^Shu background led to up- and down-regulation of various genes, with those encoding cytochrome P450 (CYP) enzymes most significantly over-represented in this group. Because GstS1 is a fly ortholog of mammalian hematopoietic prostaglandin D synthase, and in mammals CYPs are involved in the oxygenation of polyunsaturated fatty acids including prostaglandins, our results raise the intriguing possibility that bioactive lipids play a role in GstS1-mediated suppression of para^Shu phenotypes.

Genetics and metabolic responses of Artemisia annua L to the lake of phosphorus under the sparingly soluble phosphorus fertilizer: evidence from transcriptomics analysis

The medicinal herb Artemisia annua L. is prized for its capacity to generate artemisinin, which is used to cure malaria. Potentially influencing the biomass and secondary metabolite synthesis of A. annua is plant nutrition, particularly phosphorus (P). However, most soil P exist as insoluble inorganic and organic phosphates, which results to low P availability limiting plant growth and development. Although plants have developed several adaptation strategies to low P levels, genetics and metabolic responses to P status remain largely unknown. In a controlled greenhouse experiment, the sparingly soluble P form, hydroxyapatite (Ca5OH(PO4)3/CaP) was used to simulate calcareous soils with low P availability. In contrast, the soluble P form KH2PO4/KP was used as a control. A. annua's morphological traits, growth, and artemisinin concentration were determined, and RNA sequencing was used to identify the differentially expressed genes (DEGs) under two different P forms. Total biomass, plant height, leaf number, and stem diameter, as well as leaf area, decreased by 64.83%, 27.49%, 30.47%, 38.70%, and 54.64% in CaP compared to KP; however, LC-MS tests showed an outstanding 37.97% rise in artemisinin content per unit biomass in CaP contrary to KP. Transcriptome analysis showed 2015 DEGs (1084 up-regulated and 931 down-regulated) between two P forms, including 39 transcription factor (TF) families. Further analysis showed that DEGs were mainly enriched in carbohydrate metabolism, secondary metabolites biosynthesis, enzyme catalytic activity, signal transduction, and so on, such as tricarboxylic acid (TCA) cycle, glycolysis, starch and sucrose metabolism, flavonoid biosynthesis, P metabolism, and plant hormone signal transduction. Meanwhile, several artemisinin biosynthesis genes were up-regulated, including DXS, GPPS, GGPS, MVD, and ALDH, potentially increasing artemisinin accumulation. Furthermore, 21 TF families, including WRKY, MYB, bHLH, and ERF, were up-regulated in reaction to CaP, confirming their importance in P absorption, internal P cycling, and artemisinin biosynthesis regulation. Our results will enable us to comprehend how low P availability impacts the parallel transcriptional control of plant development, growth, and artemisinin production in A. annua. This study could lay the groundwork for future research into the molecular mechanisms underlying A. annua's low P adaptation.

AN1284 attenuates steatosis, lipogenesis, and fibrosis in mice with pre-existing non-alcoholic steatohepatitis and directly affects aryl hydrocarbon receptor in a hepatic cell line

Non-alcoholic steatohepatitis (NASH) is an aggressive form of fatty liver disease with hepatic inflammation and fibrosis for which there is currently no drug treatment. This study determined whether an indoline derivative, AN1284, which significantly reduced damage in a model of acute liver disease, can reverse steatosis and fibrosis in mice with pre-existing NASH and explore its mechanism of action. The mouse model of dietary-induced NASH reproduces most of the liver pathology seen in human subjects. This was confirmed by RNA-sequencing analysis. The Western diet, given for 4 months, caused steatosis, inflammation, and liver fibrosis. AN1284 (1 mg or 5 mg/kg/day) was administered for the last 2 months of the diet by micro-osmotic-pumps (mps). Both doses significantly decreased hepatic damage, liver weight, hepatic fat content, triglyceride, serum alanine transaminase, and fibrosis. AN1284 (1 mg/kg/day) given by mps or in the drinking fluid significantly reduced fibrosis produced by carbon tetrachloride injections. In human HUH7 hepatoma cells incubated with palmitic acid, AN1284 (2.1 and 6.3 ng/ml), concentrations compatible with those in the liver of mice treated with AN1284, decreased lipid formation by causing nuclear translocation of the aryl hydrocarbon receptor (AhR). AN1284 downregulated fatty acid synthase (FASN) and sterol regulatory element-binding protein 1c (SREBP-1c) and upregulated Acyl-CoA Oxidase 1 and Cytochrome P450-a1, genes involved in lipid metabolism. In conclusion, chronic treatment with AN1284 (1mg/kg/day) reduced pre-existing steatosis and fibrosis through AhR, which affects several contributors to the development of fatty liver disease. Additional pathways are also influenced by AN1284 treatment.

Decellularized Antler Cancellous Bone Matrix Material Can Serve as Potential Bone Tissue Scaffold

Due to the limited supply of autologous bone grafts, there is a need to develop more bone matrix materials to repair bone defects. Xenograft bone is expected to be used for clinical treatment due to its exact structural similarity to natural bone and its high biocompatibility. In this study, decellularized antler cancellous bone matrix (DACB) was first prepared, and then the extent of decellularization of DACB was verified by histological staining, which demonstrated that it retained the extracellular matrix (ECM). The bioactivity of DACB was assessed using C3H10T1/2 cells, revealing that DACB enhanced cell proliferation and facilitated cell adhesion and osteogenic differentiation. When evaluated by implanting DACB into nude mice, there were no signs of necrosis or inflammation in the epidermal tissues. The bone repair effect of DACB was verified in vivo using sika deer during the antler growth period as an animal model, and the molecular mechanisms of bone repair were further evaluated by transcriptomic analysis of the regenerated tissues. Our findings suggest that the low immunogenicity of DACB enhances the production of bone extracellular matrix components, leading to effective osseointegration between bone and DACB. This study provides a new reference for solving bone defects.

The Transcription Factor ZmNAC89 Gene Is Involved in Salt Tolerance in Maize (Zea mays L.)

The NAC gene family has transcription factors specific to plants, which are involved in development and stress response and adaptation. In this study, ZmNAC89, an NAC gene in maize that plays a role in saline-alkaline tolerance, was isolated and characterized. ZmNAC89 was localized in the nucleus and had transcriptional activation activity during in vitro experiments. The expression of ZmNAC89 was strongly upregulated under saline-alkaline, drought and ABA treatments. Overexpression of the ZmNAC89 gene in transgenic Arabidopsis and maize enhanced salt tolerance at the seedling stage. Differentially expressed genes (DEGs) were then confirmed via RNA-sequencing analysis with the transgenic maize line. GO analyses showed that oxidation-reduction process-regulated genes were involved in ZmNAC89-mediated salt-alkaline stress. ZmNAC89 may regulate maize saline-alkali tolerance through the REDOX pathway and ABA signal transduction pathway. From 140 inbred maize lines, 20 haplotypes and 16 SNPs were found in the coding region of the ZmNAC89 gene, including the excellent haplotype HAP20. These results contribute to a better understanding of the response mechanism of maize to salt-alkali stress and marker-assisted selection during maize breeding.

Key Regulatory Differentially Expressed Genes in the Blood of Atrial Septal Defect Children Treated With Occlusion Devices

Atrial septal defects (ASDs) are the most common types of cardiac septal defects in congenital heart defects. In addition to traditional therapy, interventional closure has become the main treatment method. However, the molecular events and mechanisms underlying the repair progress by occlusion device remain unknown. In this study, we aimed to characterize differentially expressed genes (DEGs) in the blood of patients treated with occlusion devices (metal or poly-L-lactic acid devices) using RNA-sequencing, and further validated them by qRT-PCR analysis to finally determine the expression of key mediating genes after closure of ASD treatment. The result showed that total 1,045 genes and 1,523 genes were expressed differently with significance in metal and poly-L-lactic acid devices treatment, respectively. The 115 overlap genes from the different sub-analyses are illustrated. The similarities and differences in gene expression reflect that the body response process involved after interventional therapy for ASDs has both different parts that do not overlap and the same part that crosses. The same portion of body response regulatory genes are key regulatory genes expressed in the blood of patients with ASDs treated with closure devices. The gene ontology enrichment analysis showed that biological processes affected in metal device therapy are immune response with CXCR4 genes and poly-L-lactic acid device treatment, and the key pathways are nuclear-transcribed mRNA catabolic process and proteins targeting endoplasmic reticulum process with ribosomal proteins (such as RPS26). We confirmed that CXCR4, TOB1, and DDIT4 gene expression are significantly downregulated toward the pre-therapy level after the post-treatment in both therapy groups by qRT-PCR. Our study suggests that the potential role of CXCR4, DDIT4, and TOB1 may be key regulatory genes in the process of endothelialization in the repair progress of ASDs, providing molecular insights into this progress for future studies.

Altered pathways of keratinization, extracellular matrix generation, angiogenesis, and stromal stem cells proliferation in patients with systemic sclerosis

Objective

Systemic sclerosis is characterized by endothelial dysfunction, autoimmunity abnormalities, and fibrosis of the skin and internal organs. The pathogenetic mechanisms underlying systemic sclerosis vasculopathy are still not clarified. A complex cellular and extracellular network of interactions has been studied, but it is currently unclear what drives the activation of fibroblasts/myofibroblasts and the extracellular matrix deposition.

Methods

Using RNA sequencing, the aim of the work was to identify potential functional pathways implied in systemic sclerosis pathogenesis and markers of endothelial dysfunction and fibrosis in systemic sclerosis patients. RNA-sequencing analysis was performed on RNA obtained from biopsies from three systemic sclerosis patients and three healthy controls enrolled in our University Hospital. RNA was used to generate sequencing libraries that were sequenced according to proper transcriptomic analyses. Subsequently, we performed gene set enrichment analysis of differentially expressed genes on the entire list of genes that compose the RNA-sequencing expression matrix.

Results

Gene set enrichment analysis revealed that healthy controls were characterized by gene signatures related to stromal stem cells proliferation, cytokine-cytokine receptor interaction, macrophage-enriched metabolic network, whereas systemic sclerosis tissues were enriched in signatures associated with keratinization, cornification, retinoblastoma 1 and tumor suppressor 53 signaling.

Conclusion

According to our data, RNA-sequencing and pathway analysis revealed that systemic sclerosis subjects display a discrete pattern of gene expression associated with keratinization, extracellular matrix generation, and negative regulation of angiogenesis and stromal stem cells proliferation. Further analysis on larger numbers of patients is needed; however, our findings provide an interesting framework for the development of biomarkers useful to explore potential future therapeutic approaches.