Search

Abstract Background Breast cancer is a heterogeneous disease with diverse molecular subtypes, underscoring a better understanding of its molecular features and underlying regulatory mechanisms. Therefore, identifying novel prognostic biomarkers and therapeutic targets is crucial for advancing the current standard of care for breast cancer patients. Methods Ninety-six formalin-fixed paraffin-embedded (FFPE) breast cancer samples underwent miRNAome profiling using QIAseq microRNA library kit and sequencing on Illumina platform. Mature miRNA quantification was conducted using CLC Genomics Workbench v21.0.5, while Relapse-free survival (RFS) analysis was conducted using RStudio 2023.09.1. Gain-of-function studies were conducted using miRNA mimics, while the effects of miRNA exogenous expression on cancer hallmark were assessed using 2-dimentional (2D) proliferation assay, three-dimensional (3D) organotypic culture, and live-dead staining. TargetScan database and Ingenuity Pathway Analysis (IPA) were used for miRNA target identification. Results Hierarchical clustering based on miRNA expression revealed distinct patterns in relation to PAM50 classification and identified miRNAs panels associated with luminal, HER2, and basal subtypes. hsa-miR-5683 emerged as a potential prognostic biomarker, showing a favorable correlation with RFS and suppressing tumorigenicity under 2D and 3D conditions in triple-negative breast cancer (TNBC) models. Findings were further extended to the MCF7 hormone receptor positive (HR+) model. Transcriptomic profiling of hsa-miR-5683 overexpressing TNBC cells revealed its potential role in key oncogenic pathways. Integration of downregulated genes and CRISPR-Cas9 perturbational effects identified ACLY, RACGAP1, AK4, MRPL51, CYB5B, MKRN1, TMEM230, NUP54, ANAPC13, PGAM1, and SOD1 as bona fide gene targets for hsa-miR-5683. Conclusions Our data provides comprehensive miRNA expression atlas in breast cancer subtypes and underscores the prognostic and therapeutic significance of numerous miRNAs, including hsa-miR-5683 in TNBC. The identified gene targets unravel the intricate regulatory network in TNBC progression, suggesting promising avenues for further research and targeted therapeutic interventions.

Background: Stroke refers to an abrupt neurological deficit, caused by an acute focal injury of the central nervous system via infarction or hemorrhage due to impaired vascularity, and remains among the leading causes of disability and death worldwide. Stroke is often preceded by an episode of neuronal deficit termed transient ischemic attack (TIA), which presents an effective opportunity for mitigating the risk of an eminent acute ischemic stroke (AIS). Circulating non-coding RNAs (ncRNAs) have emerged as important biomarkers for stroke, but PIWI-interacting RNAs (piRNAs), a class of small regulatory ncRNAs, have not been previously explored as diagnostic or prognostic biomarkers for stroke. Methods: We conducted comprehensive circulating piRNA profiling of AIS and TIA patients using RNA-seq on serum samples collected within 24 h of clinical diagnosis. The study cohort was divided into discovery and cross-validation datasets to identify replicated piRNAs using stringent analysis cut-offs. The expression levels of the panel of differentially regulated piRNAs between AIS and TIA patients were also compared with healthy controls. Results: We identified a panel of 10 differentially regulated piRNAs between AIS and TIA patients; hsa-piR-28272, -piR-32972, -piR-28247, -piR-24553, -piR-24552, -piR-28275, -piR-28707 and -piR-32882 were upregulated, while hsa-piR-23058 and -piR-23136 were downregulated in AIS patients. Moreover, these 10 piRNAs were also differentially expressed in AIS patients compared to healthy controls. In addition, we investigated the potential gene targets of the dysregulated piRNAs and their plausible involvement in pathophysiological processes affected in stroke. Conclusions: The imbalances in the circulating piRnome of AIS and TIA patients presented herein provide important insights into the roles of piRNAs following ischemic brain injury and potentially provide opportunities to mitigate stroke-induced mortality and morbidity.

Abstract Long noncoding RNAs (lncRNAs) play a critical role in breast cancer pathogenesis, including Triple-Negative Breast Cancer (TNBC) subtype. Identifying the lncRNA expression patterns across different breast cancer subtypes could provide valuable insights into their potential utilization as disease biomarkers and therapeutic targets. In this study, we profiled lncRNA expression in 96 breast cancer cases, revealing significant differences compared to normal breast tissue. Variations across breast cancer subtypes, including Hormone Receptor-positive (HR + ), HER2-positive (HER2 + ), HER2 + HR + , and TNBC, as well as in relation to tumor grade and patients’ age at diagnosis were observed. TNBC and HER2+ subtypes showed distinct clustering, while HER2 + HR+ tumors clustered closer to HR+ tumors based on their lncRNA profiles. Our data identified numerous enriched lncRNAs in TNBC, notably the elevated expression of LINC00960, which was subsequently validated in two additional datasets. Analysis of LINC00960 expression in an independent TNBC cohort (n = 360) revealed elevated expression of LINC00960 to correlate with cell movement, invasion, proliferation, and migration functional categories. Depletion of LINC00960 significantly reduced TNBC cell viability, colony formation, migration, and three-dimensional growth, while increasing cell death. Mechanistically, transcriptomic profiling of LINC00960-depleted cells confirmed its tumor-promoting role, likely through sponging of hsa-miR-34a-5p, hsa-miR-16-5p, and hsa-miR-183-5p, leading to the upregulation of cancer-promoting genes including BMI1, KRAS, and AKT3. Our findings highlight the distinct lncRNA expression patterns in breast cancer subtypes and underscore the crucial role for LINC00960 in promoting TNBC pathogenesis, suggesting its potential utilization as a prognostic marker and therapeutic target.

Abstract Background The mRNA subcellular localization bears substantial impact in the regulation of gene expression, cellular migration, and adaptation. However, the methods employed for experimental determination of this localization are arduous, time-intensive, and come with a high cost. Methods In this research article, we tackle the essential challenge of predicting the subcellular location of messenger RNAs (mRNAs) through Unified mRNA Subcellular Localization Predictor (UMSLP), a machine learning (ML) based approach. We embrace an in silico strategy that incorporate four distinct feature sets: kmer, pseudo k-tuple nucleotide composition, nucleotide physicochemical attributes, and the 3D sequence depiction achieved via Z-curve transformation for predicting subcellular localization in benchmark dataset across five distinct subcellular locales, encompassing nucleus, cytoplasm, extracellular region (ExR), mitochondria, and endoplasmic reticulum (ER). Results The proposed ML model UMSLP attains cutting-edge outcomes in predicting mRNA subcellular localization. On independent testing dataset, UMSLP ahcieved over 87% precision, 94% specificity, and 94% accuracy. Compared to other existing tools, UMSLP outperformed mRNALocator, mRNALoc, and SubLocEP by 11%, 21%, and 32%, respectively on average prediction accuracy for all five locales. SHapley Additive exPlanations analysis highlights the dominance of k-mer features in predicting cytoplasm, nucleus, ER, and ExR localizations, while Z-curve based features play pivotal roles in mitochondria subcellular localization detection. Availability We have shared datasets, code, Docker API for users in GitHub at: https://github.com/smusleh/UMSLP .

Abstract Background Our recent studies have demonstrated the crucial involvement of FOXA2 in the development of human pancreas. Reduction of FOXA2 expression during the differentiation of induced pluripotent stem cells (iPSCs) into pancreatic islets has been found to reduce α-and β-cell masses. However, the extent to which such changes are linked to alterations in the expression profile of long non-coding RNAs (lncRNAs) remains unraveled. Methods Here, we employed our recently established FOXA2-deficient iPSCs (FOXA2 −/− iPSCs) to investigate changes in lncRNA profiles and their correlation with dysregulated mRNAs during the pancreatic progenitor (PP) and pancreatic islet stages. Furthermore, we constructed co-expression networks linking significantly downregulated lncRNAs with differentially expressed pancreatic mRNAs. Results Our results showed that 442 lncRNAs were downregulated, and 114 lncRNAs were upregulated in PPs lacking FOXA2 compared to controls. Similarly, 177 lncRNAs were downregulated, and 59 lncRNAs were upregulated in islet cells lacking FOXA2 compared to controls. At both stages, we observed a strong correlation between lncRNAs and several crucial pancreatic genes and TFs during pancreatic differentiation. Correlation analysis revealed 12 DE-lncRNAs that strongly correlated with key downregulated pancreatic genes in both PPs and islet cell stages. Selected DE-lncRNAs were validated using RT-qPCR. Conclusions Our data indicate that the observed defects in pancreatic islet development due to the FOXA2 loss is associated with significant alterations in the expression profile of lncRNAs. Therefore, our findings provide novel insights into the role of lncRNA and mRNA networks in regulating pancreatic islet development, which warrants further investigations. Video Abstract

Abstract Background While the increased screening, changes in lifestyle, and recent advances in treatment regimen have decreased colorectal cancer (CRC) mortality, metastatic disease and recurrence remains a major clinical challenge. In the era of precision medicine, the identification of actionable novel therapeutic targets could ultimately offer an alternative treatment strategy for CRC. Methods RNA-Seq was conducted using the illumina platform, while bioinformatics analyses were conducted using CLC genomics workbench and iDEP.951. Colony forming unit, flow cytometry, and fluorescent microscopy were used to assess cell proliferation, cell cycle distribution, and cell death, respectively. The growth potential of CRC cells under 3-dimensional (3D) conditions was assessed using Matrigel. STRING database (v11.5) and Ingenuity Pathway Analysis (IPA) tool were used for network and pathway analyses. CRISPR-Cas9 perturbational effects database was used to identify potential therapeutic targets for CRC, through integration with gene-drug interaction database. Structural modeling and molecular docking were used to assess the interaction between candidate drugs and their targets. Results In the current study, we investigated the therapeutic potential of targeting TPX2, TTK, DDX39A, and LRP8, commonly upregulated genes in CRC identified through differential expression analysis in CRC and adjacent non-cancerous tissue. Targeted depletion of TPX2 and TTK impaired CRC proliferation, cell cycle progression, and organoid formation under 3D culture conditions, while suppression of DDX39A and LRP8 had modest effects on CRC colony formation. Differential expression analysis and bioinformatics on TPX2 and TTK-deficient cells identified cell cycle regulation as the hallmark associated with loss of TPX2 and TTK. Elevated expression of TPX2 and TTK correlated with an oncogenic state in tumor tissue from patients with colon adenocarcinoma, thus corroborating an oncogenic role for the TPX2/TTK network in the pathogenesis of CRC. Gene set enrichment and pathway analysis of TPX2high/TTKhigh CRC identified numerous additional gene targets as integral components of the TPX2/TTK network. Integration of TPX2/TTK enriched network with CRISPR-Cas9 functional screen data identified numerous novel dependencies for CRC. Additionally, gene-drug interaction analysis identified several druggable gene targets enriched in the TPX2/TTK network, including AURKA, TOP2A, CDK1, BIRC5, and many others. Conclusions Our data has implicated an essential role for TPX2 and TTK in CRC pathogenesis and identified numerous potential therapeutic targets and their drug interactions, suggesting their potential clinical use as a novel therapeutic strategy for patients with CRC. Video Abstract

IntroductionOsteoporosis is a systemic age-related disease characterized by reduced bone mass and microstructure deterioration, leading to increased risk of bone fragility fractures. Osteoporosis is a worldwide major health care problem and there is a need for preventive approaches.Methods and resultsApigenin and Rutaecarpine are plant-derived antioxidants identified through functional screen of a natural product library (143 compounds) as enhancers of osteoblastic differentiation of human bone marrow stromal stem cells (hBMSCs). Global gene expression profiling and Western blot analysis revealed activation of several intra-cellular signaling pathways including focal adhesion kinase (FAK) and TGFβ. Pharmacological inhibition of FAK using PF-573228 (5 μM) and TGFβ using SB505124 (1μM), diminished Apigenin- and Rutaecarpine-induced osteoblast differentiation. In vitro treatment with Apigenin and Rutaecarpine, of primary hBMSCs obtained from elderly female patients enhanced osteoblast differentiation compared with primary hBMSCs obtained from young female donors. Ex-vivo treatment with Apigenin and Rutaecarpine of organotypic embryonic chick-femur culture significantly increased bone volume and cortical thickness compared to control as estimated by μCT-scanning.DiscussionOur data revealed that Apigenin and Rutaecarpine enhance osteoblastic differentiation, bone formation, and reduce the age-related effects of hBMSCs. Therefore, Apigenin and Rutaecarpine cellular treatment represent a potential strategy for maintaining hBMSCs health during aging and osteoporosis.

Abstract Previous studies have suggested that breast cancer (BC) from the Middle East and North Africa (MENA) is presented at younger age with advanced tumor stage, indicating underlying biological differences. Given the scant transcriptomic data on BC from the MENA region and to better understand the biology of this disease, we performed mRNA and microRNA (miRNA) transcriptomic profiling on a local cohort of BC (n = 96) from Qatar. Our data revealed the differentially expressed genes and miRNAs as function of BC molecular subtypes (HR+, HER2+, HER2+HR+, and TNBC), tumor grade (GIII vs GI-II), patients’ age (young (≤40) vs old (>40)), and ethnicity (MENA vs non-MENA). Our profiling data revealed close similarity between TNBC and HER2+, while the transcriptome of HER2+HR+ tumor was resemblant of that from HR+ tumors. Network analysis identified complex miRNA-mRNA regulatory networks in each BC molecular subtype, in high vs low grade tumors, in tumors from young vs old patients, and in tumors from MENA vs non-MENA, thus implicating miRNA-mediated gene regulation as an essential mechanism in shaping the transcriptome of BC. Integration of our transcriptomic data with CRISPR-Cas9 functional screen data and the OncoKB database identified numerous dependencies and therapeutic vulnerabilities in each BC molecular subtype, while CDC123 was functionally validated as potential therapeutic target for TNBC. Cox regression survival analyses identified mRNA and miRNA-based signatures predicative of worse and better relapse free survival (RFS), which were validated in larger BC cohorts. Our data provides comprehensive transcriptomic profiling and unraveled the miRNA-mRNA regulatory networks in BC patients from the region and identified novel actionable gene targets, employing integrated approach. Findings from the current study have potential implications to improve the current standard-of-care for BC from the MENA as well as patients from other ethnicities.

BackgroundDementia is a debilitating neurological disease affecting millions of people worldwide. The exact mechanisms underlying the initiation and progression of the disease remain to be fully defined. There is an increasing body of evidence for the role of immune dysregulation in the pathogenesis of dementia, where blood-borne autoimmune antibodies have been studied as potential markers associated with pathological mechanisms of dementia.MethodsThis study included plasma from 50 cognitively normal individuals, 55 subjects with MCI (mild cognitive impairment), and 22 subjects with dementia. Autoantibody profiling for more than 1,600 antigens was performed using a high throughput microarray platform to identify differentially expressed autoantibodies in MCI and dementia.ResultsThe differential expression analysis identified 33 significantly altered autoantibodies in the plasma of patients with dementia compared to cognitively normal subjects, and 38 significantly altered autoantibodies in the plasma of patients with dementia compared to subjects with MCI. And 20 proteins had significantly altered autoantibody responses in MCI compared to cognitively normal individuals. Five autoantibodies were commonly dysregulated in both dementia and MCI, including anti-CAMK2A, CKS1B, ETS2, MAP4, and NUDT2. Plasma levels of anti-ODF3, E6, S100P, and ARHGDIG correlated negatively with the cognitive performance scores (MoCA) (r2 –0.56 to −0.42, value of p

Abstract Background Subcellular localization of messenger RNA (mRNAs) plays a pivotal role in the regulation of gene expression, cell migration as well as in cellular adaptation. Experiment techniques for pinpointing the subcellular localization of mRNAs are laborious, time-consuming and expensive. Therefore, in silico approaches for this purpose are attaining great attention in the RNA community. Methods In this article, we propose MSLP, a machine learning-based method to predict the subcellular localization of mRNA. We propose a novel combination of four types of features representing k-mer, pseudo k-tuple nucleotide composition (PseKNC), physicochemical properties of nucleotides, and 3D representation of sequences based on Z-curve transformation to feed into machine learning algorithm to predict the subcellular localization of mRNAs. Results Considering the combination of the above-mentioned features, ennsemble-based models achieved state-of-the-art results in mRNA subcellular localization prediction tasks for multiple benchmark datasets. We evaluated the performance of our method in ten subcellular locations, covering cytoplasm, nucleus, endoplasmic reticulum (ER), extracellular region (ExR), mitochondria, cytosol, pseudopodium, posterior, exosome, and the ribosome. Ablation study highlighted k-mer and PseKNC to be more dominant than other features for predicting cytoplasm, nucleus, and ER localizations. On the other hand, physicochemical properties and Z-curve based features contributed the most to ExR and mitochondria detection. SHAP-based analysis revealed the relative importance of features to provide better insights into the proposed approach. Availability We have implemented a Docker container and API for end users to run their sequences on our model. Datasets, the code of API and the Docker are shared for the community in GitHub at: https://github.com/smusleh/MSLP .

IntroductionAutism spectrum disorder (ASD) is a neurodevelopmental condition characterized by defects in two core domains, social/communication skills and restricted/repetitive behaviors or interests. There is no approved biomarker for ASD diagnosis, and the current diagnostic method is based on clinical manifestation, which tends to vary vastly between the affected individuals due to the heterogeneous nature of ASD. There is emerging evidence that supports the implication of the immune system in ASD, specifically autoimmunity; however, the role of autoantibodies in ASD children is not yet fully understood.Materials and methodsIn this study, we screened serum samples from 93 cases with ASD and 28 healthy controls utilizing high-throughput KoRectly Expressed (KREX) i-Ome protein-array technology. Our goal was to identify autoantibodies with differential expressions in ASD and to gain insights into the biological significance of these autoantibodies in the context of ASD pathogenesis.ResultOur autoantibody expression analysis identified 29 differential autoantibodies in ASD, 4 of which were upregulated and 25 downregulated. Subsequently, gene ontology (GO) and network analysis showed that the proteins of these autoantibodies are expressed in the brain and involved in axonal guidance, chromatin binding, and multiple metabolic pathways. Correlation analysis revealed that these autoantibodies negatively correlate with the age of ASD subjects.ConclusionThis study explored autoantibody reactivity against self-antigens in ASD individuals' serum using a high-throughput assay. The identified autoantibodies were reactive against proteins involved in axonal guidance, synaptic function, amino acid metabolism, fatty acid metabolism, and chromatin binding.

Our understanding of the function of long non-coding RNAs (lncRNAs) in health and disease states has evolved over the past decades due to the many advances in genome research. In the current study, we characterized the lncRNA transcriptome enriched in triple-negative breast cancer (TNBC, n = 42) and estrogen receptor (ER+, n = 42) breast cancer compared to normal breast tissue (n = 56). Given the aggressive nature of TNBC, our data revealed selective enrichment of 57 lncRNAs in TNBC. Among those, AC099850.4 lncRNA was chosen for further investigation where it exhibited elevated expression, which was further confirmed in a second TNBC cohort (n = 360) where its expression correlated with a worse prognosis. Network analysis of AC099850.4high TNBC highlighted enrichment in functional categories indicative of cell cycle activation and mitosis. Ingenuity pathway analysis on the differentially expressed genes in AC099850.4high TNBC revealed the activation of the canonical kinetochore metaphase signaling pathway, pyridoxal 5'-phosphate salvage pathway, and salvage pathways of pyrimidine ribonucleotides. Additionally, upstream regulator analysis predicted the activation of several upstream regulator networks including CKAP2L, FOXM1, RABL6, PCLAF, and MITF, while upstream regulator networks of TP53, NUPR1, TRPS1, and CDKN1A were suppressed. Interestingly, elevated expression of AC099850.4 correlated with worse short-term relapse-free survival (log-rank p = 0.01). Taken together, our data are the first to reveal AC099850.4 as an unfavorable prognostic marker in TNBC, associated with more aggressive clinicopathological features, and suggest its potential utilization as a prognostic biomarker and therapeutic target in TNBC.

Triple-negative breast cancer (TNBC) resistance to neoadjuvant chemotherapy (NAC) represents a major clinical challenge; therefore, delineating tumor heterogeneity can provide novel insight into resistance mechanisms and potential therapeutic targets. Herein, we identified the transcriptional landscape associated with TNBC resistance to NAC at the single-cell level by analyzing publicly available transcriptome data from more than 5,000 single cells derived from four extinction (responders) and four persistence (non-responders) patients, revealing remarkable tumor heterogeneity. Employing iterative clustering and guide-gene selection (ICGS) and uniform manifold approximation and projection (UMAP), we classified TNBC single cells into several clusters based on their distinct gene signatures. The presence of clusters indicative of immune cell activation was a hallmark of the extinction group pre-NAC, while post NAC, the extinction tissue consisted mostly of breast, omental fat, and fibroblasts. The persistent gene signatures of pre-NAC resembled the gene signature of lung epithelial, mammary, and salivary glands and acute myeloid leukemia blast cells, which were associated with enhanced cellular movement and activation of FOXM1, NOTCH1, and MYC and suppression of tumor necrosis factor (TNF) and IFNG mechanistic networks. Multivariate survival analysis identified persistence-derived three-gene signature (KIF5BhighHLA-ClowIGHG2low) predictive of relapse-free survival (hazard ratio [HR]: 2.2 [1.6–3.2, p < 0.0001]) in a second cohort of 360 TNBC patients. Mechanistically, loss of function of several upregulated genes in the persistent group (BYSL, FDPS, ENO1, MED20, MRPL9, MRPL37, NDUFB11, PMVK, MYC, and GSTP1) inhibited MDA-MB-231 and BT-549 TNBC models’ colony-forming unit (CFU) potential and enhanced their sensitivity to paclitaxel. Our data unraveled the transcriptional portrait associated with NAC resistance, identified several key genes, and suggested their potential utilization as prognostic markers and therapeutic targets in TNBC.

Abstract Background DNA methylation plays a crucial role in multiple cellular processes such as gene regulation, chromatin stability, and genetic imprinting. In mammals, DNA methylation is achieved by DNA methyltransferases (DNMTs). A number of studies have associated alterations in DNMT activity to tumorigenesis; however, the exact role of DNMTs in shaping the genome in triple negative breast cancer (TNBC) is still being unraveled. Methods In the current study, we employed two DNMT inhibitors (Decitabine and 5-Azacytidine), two TNBC models (MDA-MB-231 and BT-549) and whole transcriptome RNA-Seq and characterized the transcriptional alterations associated with DNMT inhibition. Colony forming unit (CFU), flow cytometry, and fluorescent microscopy were used to assess cell proliferation, cell cycle distribution, and cell death, respectively. Ingenuity pathway analysis (IPA) was used for network and pathway analyses. Results Remarkably, DNMT inhibition induced the expression of genes involved in endoplasmic reticulum response to stress, response to unfolder protein, as well as cobalamin metabolic processes. In contrast, suppression of cellular processes related to cell cycle and mitosis were hallmarks of DNMT inhibition. Concordantly, DNMT inhibition led to significant inhibition of TNBC cell proliferation, G2-M cell cycle arrest and induction of cell death. Mechanistically, DNMT inhibition activated TP53, NUPR1, and NFkB (complex) networks, while RARA, RABL6, ESR1, FOXM1, and ERBB2 networks were suppressed. Our data also identified the long noncoding RNA (lncRNA) transcriptional portrait associated with DNMT inhibition and identified 25 commonly upregulated and 60 commonly downregulated lncRNAs in response to Decitabine and 5-Azacytidinec treatment in both TNBC models. TPT1-AS1 was the most highly induced (6.3 FC), while MALAT1 was the most highly suppressed (− 7.0 FC) lncRNA in response to DNMT inhibition. Conclusions Taken together, our data provides a comprehensive view of transcriptome alterations in the coding and noncoding transcriptome in TNBC in response to DNMT inhibition.

Stroke is the second leading cause of global mortality and continued efforts aim to identify predictive, diagnostic, or prognostic biomarkers to reduce the disease burden. Circulating microRNAs (miRNAs) have emerged as potential biomarkers in stroke. We performed comprehensive circulating miRNA profiling of ischemic stroke patients with or without type 2 diabetes mellitus (T2DM), an important risk factor associated with worse clinical outcomes in stroke. Serum samples were collected within 24 h of acute stroke diagnosis and circulating miRNAs profiled using RNA-Seq were compared between stroke patients with T2DM (SWDM; n = 92) and those without T2DM (SWoDM; n = 98). Our analysis workflow involved random allocation of study cohorts into discovery (n = 96) and validation (n = 94) datasets. Five miRNAs were found to be differentially regulated in SWDM compared to SWoDM patients. Hsa-miR-361-3p and -664a-5p were downregulated, whereas miR-423-3p, -140-5p, and -17-3p were upregulated. We also explored the gene targets of these miRNAs and investigated the downstream pathways associated with them to decipher the potential pathways impacted in stroke with diabetes as comorbidity. Overall, our novel findings provide important insights into the differentially regulated miRNAs, their associated pathways and potential utilization for clinical benefits in ischemic stroke patients with diabetes.

Abstract Colorectal cancer (CRC) remains a global disease burden and a leading cause of cancer related deaths worldwide. The identification of aberrantly expressed messenger RNA (mRNA), long non-coding RNA (lncRNA), and microRNA (miRNA), and the resulting molecular interactions and signaling networks is essential for better understanding of CRC, identification of novel diagnostic biomarkers and potential development of therapeutic interventions. Herein, we performed microRNA (miRNA) sequencing on fifteen CRC and their non-tumor adjacent tissues and whole transcriptome RNA-Seq on six paired samples from the same cohort and identified alterations in miRNA, mRNA, and lncRNA expression. Computational analyses using Ingenuity Pathway Analysis (IPA) identified multiple activated signaling networks in CRC, including ERBB2, RABL6, FOXM1, and NFKB networks, while functional annotation highlighted activation of cell proliferation and migration as the hallmark of CRC. IPA in combination with in silico prediction algorithms and experimentally validated databases gave insight into the complex associations and interactions between downregulated miRNAs and upregulated mRNAs in CRC and vice versa. Additionally, potential interaction between differentially expressed lncRNAs such as H19, SNHG5, and GATA2-AS1 with multiple miRNAs has been revealed. Taken together, our data provides thorough analysis of dysregulated protein-coding and non-coding RNAs in CRC highlighting numerous associations and regulatory networks thus providing better understanding of CRC.

Abstract TGFβ signaling plays crucial role during development and cancer, however the role for TGFβ signaling in regulating the noncoding part of the human genome in triple negative breast cancer (TNBC) is still being unraveled. Herein, we provide the transcriptional landscape of TNBC in response to TGFβ activation and subsequent inhibition employing SB431542, selective TGFβ1 Receptor ALK5 Inhibitor. Our data revealed 72 commonly upregulated [fold change (FC) ≥ 2.0], including PLAU, TPM1, TAGLN, COL1A1, TGFBI, and SNAI1, and 53 downregulated (FC ≤ 2.0) protein coding genes in BT-549 and MDA-MB-231 models in response to TGFβ1 activation. Alignment to the geocode (V33) identified 41 upregulated (FC ≥ 2.0) and 22 downregulated (FC ≤ 2.0) long non-coding RNA (lncRNA) in response to TGFβ1 activation, which were inhibited by concurrent treatment with SB431542. To place our data from the in vitro models into their clinical context, we identified AC015909.1, AC013451.1, CYP1B1-AS1, AC004862.1, LINC01824, AL138828.1, B4GALT1-AS1, AL353751.1, AC090826.3, AC104695.4, ADORA2A-AS1, PTPRG-AS1, LINC01943, AC026954.3, TPM1-AS, ZFPM2-AS1, AC007362.1, AC112721.2, MALAT1, AL513314.2, AC112721.1, AC010343.3, LINC01711, and MAP3K2-DT lncRNA expression to positively correlate with TGFβ1 expression in a cohort of 360 TNBC patients. To provide mechanistic insight into lncRNA regulation by TGFβ signaling, SMAD2/3 ChIp-Seq data from BT-549 TNBC model retrieved from Gene Expression Omnibus (GEO) revealed direct binding of SMAD2/SMAD3 to the promoter of AC112721.1, AC112721.2, MALAT1, HHIP-AS1, LINC00472, and SLC7A11, suggesting their direct regulation by TGFβ1/SMAD2/SMAD3 pathway. Interestingly, AC112721.1, AC112721.2 exhibited higher expression in TNBC compared to normal breast tissue suggesting a possible role for those lncRNA in TNBC biology. Our miRNA analysis in the BT-549 model in response to exogenous TGFB1 revealed several affected miRNAs (2.0 ≤ FC ≤ 2.0), whose expression pattern was reversed in the presence of SB431542, suggesting those miRNA as plausible targets for TGFβ regulation. In particular, we observed hsa-miR-1275 to be downregulated in response to TGFB1 which was highly predicted to regulate PCDH1, FIBCD1, FXYD7, GDNF, STC1, EDN1, ZSWIM4, FGF1, PPP1R9B, NUAK1, PALM2AKAP2, IGFL3, and SPOCK1 whose expression were upregulated in response to TGFβ1 stimulus. On the other hand, hsa-miR-181b-5p was among the top upregulated miRNAs in response to TGFB1, which is also predicted to regulate CDKN1B, TNFRSF11B, SIM1, and ARSJ in the BT-549 model. Taken together, our data is the first to provide such in depth analysis of lncRNA and miRNA epigenetic changes in response to TGFβ signaling in TNBC.

Triple-negative breast cancer (TNBC) patients exhibit variable responses to chemotherapy, suggesting an underlying molecular heterogeneity. In the current study, we analyzed publicly available transcriptome data from 360 TNBC and 88 normal breast tissues, which revealed activation of nucleosome and cell cycle as the hallmarks of TNBC. Mechanistic network analysis identified activation of FOXM1 and ERBB2, and suppression of TP53 and NURP1 networks in TNBC. Employing Iterative Clustering and Guide-gene Selection (ICGS), Uniform Manifold Approximation and Projection (UMAP), and dimensionality reduction analyses, we classified TNBC into seven molecular subtypes, each exhibiting a unique molecular signature, including immune infiltration (CD19, CD8, and macrophages) and mesenchymal signature, which correlated with variable disease outcomes in a larger cohort (1,070) of BC. Mechanistically, depletion of TTK, TPX2, UBE2C, CDCA7, MELK, NFE2L3, DDX39A, and LRP8 led to substantial inhibition of colony formation of TNBC models, which was further enhanced in the presence of paclitaxel. Our data provide novel insights into the molecular heterogeneity of TNBC and identified TTK, TPX2, UBE2C, and LRP8 as main drivers of TNBC tumorigenesis.

Abstract Cumulative evidence suggests added benefit for neoadjuvant chemotherapy (NAC) in a subset of triple-negative breast cancer (TNBC) patients. Herein we identified the long noncoding RNA (lncRNA) transcriptional landscape associated with TNBC resistance to NAC, employing 1758 single cells from three extinction and three persistence TNBC patients. Using Iterative Clustering and Guide-gene Selection (ICGS) and uniform manifold approximation and projection (UMAP) dimensionality reduction analysis, we observed single cells derived from each patient to largely cluster together. Comparing the lncRNA transcriptome from single cells through the course of NAC treatment revealed minimal overlap based on lncRNA transcriptome, suggesting substantial effects of NAC on lncRNA transcription. The differential analysis revealed upregulation of 202 and downregulation of 19 lncRNAs in the persistence group, including upregulation of five different transcripts encoding for the MALAT1 lncRNA. CRISPR/Cas9-mediated MALAT1 promoter deletion in BT-549 TNBC model enhanced sensitivity to paclitaxel and doxorubicin, suggesting a role for MALAT1 in conferring resistance. Mechanistically, whole transcriptome analysis of MALAT1-KO cells revealed multiple affected mechanistic networks as well as oxidative phosphorylation canonical and angiogenesis functional category. Interestingly, lncRNA profiling of MALAT1-depleted TNBC also revealed a number of altered lncRNAs in response to MALAT1 deletion, suggesting a reciprocal relationship between MALAT1 and a number of lncRNAs, including NEAT1, USP3-AS1, and LINC-PINT, in TNBC. Elevated expression of MALAT1, USP3-AS1, and LINC-PINT correlated with worse clinical outcomes in BC patients. Our data revealed the lncRNA transactional portrait and highlighted a complex regulatory network orchestrated by MALAT1 in the context of TNBC resistance to NAC therapy.

Breast cancer (BC) is a heterogeneous disease, which is primarily classified according to hormone receptors and HER2 expression. Despite the many advances in BC diagnosis and management, the identification of novel actionable therapeutic targets expressed by cancerous cells has always been a daunting task due to the large heterogeneity of the disease and the presence of non-cancerous cells (i.e., immune cells and stromal cells) within the tumor microenvironment. In the current study, we employed computational algorithms to decipher the cellular composition of estrogen receptor-positive (ER+), HER2+, ER+HER2+, and triple-negative BC (TNBC) subtypes from a total of 49,899 single cells’ publicly available transcriptomic data derived from 26 BC patients. Restricting the analysis to EPCAM+Lin− tumor epithelial cells, we identified the enriched gene sets in each BC molecular subtype. Integration of single-cell transcriptomic with CRISPR-Cas9 functional screen data identified 13 potential therapeutic targets for ER+, 44 potential therapeutic targets for HER2+, and 29 potential therapeutic targets for TNBC. Interestingly, several of the identified therapeutic targets outperformed the current standard of care for each BC subtype. Given the aggressive nature and lack of targeted therapies for TNBC, elevated expression of ENO1, FDPS, CCT6A, TUBB2A, and PGK1 predicted worse relapse-free survival (RFS) in basal BC (n = 442), while elevated expression of ENO1, FDPS, CCT6A, and PGK1 was observed in the most aggressive BLIS TNBC subtype. Mechanistically, targeted depletion of ENO1 and FDPS halted TNBC cell proliferation, colony formation, and organoid tumor growth under 3-dimensional conditions and increased cell death, suggesting their potential use as novel therapeutic targets for TNBC. Differential expression and gene set enrichment analysis in TNBC revealed enrichment in the cycle and mitosis functional categories in FDPShigh, while ENO1high was associated with numerous functional categories, including cell cycle, glycolysis, and ATP metabolic processes. Taken together, our data are the first to unravel the unique gene signatures and to identify novel dependencies and therapeutic vulnerabilities for each BC molecular subtype, thus setting the foundation for the future development of more effective targeted therapies for BC.

Abstract Targeting regulatory signaling pathways that control human bone marrow stromal (skeletal or mesenchymal) stem cell (hBMSC) differentiation and lineage fate determination is gaining momentum in the regenerative medicine field. Therefore, to identify the central regulatory mechanism of osteoblast differentiation of hBMSCs, the molecular phenotypes of two clonal hBMSC lines exhibiting opposite in vivo phenotypes, namely, bone forming (hBMSC+bone) and non-bone forming (hBMSC−Bone) cells, were studied. Global transcriptome analysis revealed significant downregulation of several TGFβ responsive genes, namely, TAGLN, TMP1, ACTA2, TGFβ2, SMAD6, SMAD9, BMP2, and BMP4 in hBMSC−Bone cells and upregulation on SERPINB2 and NOG. Transcriptomic data was associated with marked reduction in SMAD2 protein phosphorylation, which thereby implies the inactivation of TGFβ and BMP signaling in those cells. Concordantly, activation of TGFβ signaling in hBMSC−Bone cells using either recombinant TGFβ1 protein or knockdown of SERPINB2 TGFβ-responsive gene partially restored their osteoblastic differentiation potential. Similarly, the activation of BMP signaling using exogenous BMP4 or via siRNA-mediated knockdown of NOG partially restored the differentiation phenotype of hBMSC−Bone cells. Concordantly, recombinant NOG impaired ex vivo osteoblastic differentiation of hBMSC+Bone cells, which was associated with SERBINB2 upregulation. Our data suggests the existence of reciprocal relationship between TGFB and BMP signaling that regulates hBMSC lineage commitment and differentiation, whilst provide a plausible strategy for generating osteoblastic committed cells from hBMSCs for clinical applications.

The longevity of COVID-19 as a global pandemic, and the devastating effects it has had on certain subsets of individuals thus far has highlighted the importance of identifying blood-based biomarkers associated with disease severity. We employed computational and transcriptome analyses of publicly available datasets from PBMCs from 126 patients with COVID-19 admitted to ICU (n = 50), COVID-19 not admitted to ICU (n = 50), non-COVID-19 admitted to ICU (n = 16) and non-COVID-19 not admitted to ICU (n = 10), and utilized the Gencode V33 assembly to analyze protein coding mRNA and long noncoding RNA (lncRNA) transcriptomes in the context of disease severity. Our data identified several aberrantly expressed mRNA and lncRNA based biomarkers associated with SARS-CoV-2 severity, which in turn significantly affected canonical, upstream, and disease functions in each group of patients. Immune, interferon, and antiviral responses were severely suppressed in COVID-19 patients admitted to ICU versus those who were not admitted to ICU. Our data suggests a possible therapeutic approach for severe COVID-19 through administration of interferon therapy. Delving further into these biomarkers, roles and their implications on the onset and disease severity of COVID-19 could play a crucial role in patient stratification and identifying varied therapeutic options with diverse clinical implications.

Triple negative breast cancer (TNBC) represents a diverse group of cancers based on their gene expression profiles. While the current mRNA-based classification of TNBC has contributed to our understanding of the heterogeneity of this disease, whether such heterogeneity can be resolved employing a long noncoding RNA (lncRNA) transcriptome has not been established thus far. Herein, we used iterative clustering and guide-gene selection (ICGS) and uniform manifold approximation and projection (UMAP) dimensionality reduction analysis on a large cohort of TNBC transcriptomic data (TNBC = 360, normal = 88) and classified TNBC into four main clusters: LINC00511-enriched, LINC00393-enriched, FIRRE-enriched, and normal tissue-like. Delving into associated gene expression profiles revealed remarkable differences in canonical, casual, upstream, and functional categories among different lncRNA-derived TNBC clusters, suggesting functional consequences for altered lncRNA expression. Correlation and survival analysis comparing mRNA- and lncRNA-based clustering revealed similarities and differences between the two classification approaches. To provide insight into the potential role of the identified lncRNAs in TNBC biology, CRISPR-Cas9 mediated LINC00511 promoter deletion reduced colony formation and enhanced the sensitivity of TNBC cells to paclitaxel, suggesting a role for LINC00511 in conferring tumorigenicity and resistance to therapy. Our data revealed a novel lncRNA-based classification of TNBC and suggested their potential utilization as disease biomarkers and therapeutic targets.

Elevated levels of myeloid-derived suppressor cells (MDSCs), including polymorphonuclear MDSCs (PMN-MDSCs) and immature MDSCs (I-MDSCs), are usually associated with disease progression in cancer patients, including colorectal cancer (CRC). However, biological mechanisms and molecular pathways regulated by MDSC subpopulations in the CRC tumor microenvironment (TME) have not been fully investigated. In this study, we performed transcriptomic analysis of tumor-infiltrating I-MDSCs and PMN-MDSCs isolated from tumor tissues of six CRC patients, compared to antigen-presenting cells (APCs). We also compared the transcriptomic profiles of tumor-infiltrating PMN-MDSCs to I-MDSCs. Our results showed different molecular pathways regulated by each MDSC subset, potentially reflecting their phenotypical/molecular/functional characteristics in the CRC TME. Moreover, we identified gene signatures in PMN-MDSC and I-MDSC of poor overall survival (OS) and disease-free survival (DFS) using the Cancer Genome Atlas (TCGA) dataset from patients with colon adenocarcinoma (COAD). However, functional studies are required to validate these findings.

Myeloid-derived suppressor cells (MDSCs) promote tumor immune evasion and favor tumorigenesis by activating various tumor-promoting downstream signals. MDSC expansion is evident in the circulation and tumor microenvironment of many solid tumors including colorectal cancer (CRC). We have recently reported the transcriptomic profiles of tumor-infiltrating MDSCs in CRC patients and uncovered pathways, which could potentially assist tumor progression. In this study, we sorted different subsets of circulating MDSCs in CRC patients and investigated their transcriptomic profiles in order to disclose pathways, which could potentially contribute to disease progression. The sorted subsets included polymorphonuclear/granulocytic MDSCs (PMN-MDSCs), immature MDSCs (I-MDSCs), and monocytic MDSCs (M-MDSCs). Our functional annotation analyses revealed that multiple pathways including DNA damage-, chemotaxis-, apoptosis-, mitogen-activated protein kinase-, transforming growth factor β-, and myeloid differentiation–related transcripts were higher in PMN-MDSCs, compared with monocytic antigen-presenting cells (APCs) or I-MDSCs. Furthermore, genes related to Janus kinase (JAK)–signal transducer and activator of transcription (STAT) were also elevated in PMN-MDSCs. These data suggest that upregulation of JAK-STAT pathway could trigger multiple downstream targets in PMN-MDSCs, which favor tumor progression. Additionally, we found that pathways including phosphatidyl inositol 3-kinase (PI3K), interleukin 6, and TGF-β in M-MDSCs and cell cycle–related pathways in I-MDSCs were upregulated, compared with monocytic APCs. Moreover, acetylation-related genes were upregulated in both PMN-MDSCs and M-MDSCs. This latter finding implicates that epigenetic modifications could also play a role in the regulation of multiple tumor-promoting genes in PMN-MDSCs and M-MDSCs. Taken together, this study reveals various signaling pathways, which regulate the function of MDSC subsets in the circulation of CRC patients. However, functional studies are warranted to support these findings.

Breast cancer (BC) is the foremost cause of cancer-related deaths in women. BC patients are oftentimes presented with lymph node metastasis (LNM), which increases their risk of recurrence. Compelling data have recently implicated microRNAs in promoting BC metastasis. Therefore, the identification of microRNA (miRNA)-based molecular signature associated with LNM could provide an opportunity for a more personalized treatment for BC patients with high risk of LNM. In current study, we performed comprehensive miRNA profiling in matched primary breast and LNM and identified 40 miRNAs, which were differentially expressed in LNM compared to primary tumors. The expression of 14 miRNAs (Up: hsa-miR-155-5p, hsa-miR-150-5p, hsa-miR-146a-5p, hsa-miR-142-5p and down: hsa-miR-200a-3p, hsa-miR-200b-3p, hsa-miR-200c-3p, hsa-miR-205-5p, hsa-miR-210-3p, hsa-miR-214-3p, hsa-miR-141-3p, hsa-miR-127-3p, hsa-miR-125a-5p, and hsa-let-7c-5p) was subsequently validated in a second cohort of 32 breast and 32 matched LNM tumor tissues. Mechanistically, forced expression of hsa-miR-205-5p, or hsa-miR-214-3p epigenetically inhibited MDA-MB-231 cell proliferation, colony formation, and cell migration. Global gene expression profiling on MDA-MB-231 cells overexpressing hsa-miR-205-5p, or hsa-miR-214-3p in combination with in silico target prediction and ingenuity pathway analyses identified multiple bona fide targets for hsa-miR-205-5p, hsa-miR-214-3p affecting cellular proliferation and migration. Interestingly, interrogation of the expression levels of hsa-miR-205 and hsa-miR-214 in the METABRIC breast cancer dataset revealed significantly poor overall survival in patients with downregulated expression of miR-205 [HR = 0.75 (0.61–0.91)], p = 0.003 and hsa-miR-214 [HR = 0.74 (0.59–0.93) p = 0.008]. Our data unraveled the miRNA-transcriptional landscape associated with LNM and provide novel insight on the role of several miRNAs in promoting BC LNM, and suggest their potential utilization in the clinical management of BC patients.

Gold complex bis(diethyldithiocarbamato-gold(I)) bis(diphenylphosphino) methane (BDG-I) is cytotoxic toward different cancer cell lines. We compared the cytotoxic effect of BDG-I with that of cisplatin in the A549 lung cancer cell line. Additionally, we investigated the molecular mechanism underlying the toxic effect of BDG-I toward the A549 cell line and the identification of cancer-related miRNAs likely to be involved in killing the lung cancer cells. Further, X-ray crystallographic data of the compound were acquired. Using microarray, global miRNA expression profiling in BDG-I-treated A549 cells revealed 64 upregulated and 86 downregulated miRNAs, which targeted 4689 and 2498 genes, respectively. Biological network connectivity of the miRNAs was significantly higher for the upregulated miRNAs than for the downregulated miRNAs. Two of the 10 most upregulated miRNAs (hsa-mir-20a-5p and hsa-mir-15b-5p) were associated with lung cancer. AmiGo2 server and Panther pathway analyses indicated significant enrichment in transcription regulation of miRNA target genes that promote intrinsic kinase-mediated signaling, TGF-β, and GnRH signaling pathways, as well as oxidative stress responses. BDG-I crystal structure X-ray diffraction studies revealed gold–gold intramolecular interaction [Au…Au = 3.1198 (3) Å] for a single independent molecule, reported to be responsible for its activity against cancer. Our present study sheds light on the development of novel gold complex with favorable anti-cancer therapeutic functionality. Keywords: Lung cancer, Chrysotherapeutic agents, Gold, miRNA

Abstract Background Better understanding of the signaling pathways that regulate human bone marrow stromal stem cell (hBMSC) differentiation into bone-forming osteoblasts is crucial for their clinical use in regenerative medicine. Chemical biology approaches using small molecules targeting specific signaling pathways are increasingly employed to manipulate stem cell differentiation fate. Methods We employed alkaline phosphatase activity and staining assays to assess osteoblast differentiation and Alizarin R staining to assess mineralized matrix formation of cultured hBMSCs. Changes in gene expression were assessed using an Agilent microarray platform, and data normalization and bioinformatics were performed using GeneSpring software. For in vivo ectopic bone formation experiments, hMSCs were mixed with hydroxyapatite–tricalcium phosphate granules and implanted subcutaneously into the dorsal surface of 8-week-old female nude mice. Hematoxylin and eosin staining and Sirius Red staining were used to detect bone formation in vivo. Results We identified several compounds which inhibited osteoblastic differentiation of hMSCs. In particular, we identified ruxolitinib (INCB018424) (3 μM), an inhibitor of JAK-STAT signaling that inhibited osteoblastic differentiation and matrix mineralization of hMSCs in vitro and reduced ectopic bone formation in vivo. Global gene expression profiling of ruxolitinib-treated cells identified 847 upregulated and 822 downregulated mRNA transcripts, compared to vehicle-treated control cells. Bioinformatic analysis revealed differential regulation of multiple genetic pathways, including TGFβ and insulin signaling, endochondral ossification, and focal adhesion. Conclusions We identified ruxolitinib as an important regulator of osteoblast differentiation of hMSCs. It is plausible that inhibition of osteoblast differentiation by ruxolitinib may represent a novel therapeutic strategy for the treatment of pathological conditions caused by accelerated osteoblast differentiation and mineralization.

Recent years have witnessed major progress in development of novel therapeutic agents such as chemotherapy, targeted therapy and immune checkpoint inhibitors for breast cancer. However, cancer-related death remains high especially in triple-negative breast cancer (TNBC) due limited therapeutic options. Development of targeted therapies for TNBC requires better understanding of biology and signaling networks that promote disease progression. Fascin, an actin bundling protein, was identified as a key regulator of many signaling pathways that contribute to breast cancer progression. Herein, fascin ShRNA was used to generate stable fascin knockdown (FSCN1KD) in the MDA-MB-231 TNBC cell line and then were subjected to comprehensive mRNA and miRNA transcriptome analysis. We identified 129 upregulated and 114 downregulated mRNA transcripts, while 14 miRNAs were differentially expressed in FSCN1KD. Ingenuity pathway analysis (IPA) was used to predict the impact of differentially expressed transcripts on signaling pathways and functional categories and to construct miRNA-mRNA regulatory networks in the context of FSCN1 knockdown. Compared to FSCN1KD, fascin-positive (FSCN1CON) breast cancer cells showed enrichment in genes promoting cellular proliferation, migration, survival, DNA replication and repair. Expression of FSCN1high (identified in BRCA dataset from TCGA) in conjunction with elevated expression of the top 10 upregulated or decreased expression of the top 10 downregulated genes (identified in our FSCN1CON vs. FSCN1KD) correlates with worst survival outcome. Taken together, these data confirmed fascin’s role in promoting TNBC progression, and identified a novel opportunity for therapeutic interventions via targeting those FSCN1-related transcripts.

Breast cancer (BC) is the leading cause of cancer-related death in women. Therefore, a better understanding of BC biology and signaling pathways might lead to the development of novel biomarkers and targeted therapies. Although a number of transcriptomic studies have been performed on breast cancer patients from various geographic regions, there are almost no such comprehensive studies performed on breast cancer from patients in the gulf region. This study aimed to provide a better understanding of the altered molecular networks in BC from the gulf region. Herein, we compared the transcriptome of BC to adjacent normal tissue from six BC patients and identified 1,108 upregulated and 518 downregulated transcripts. A selected number of genes from the RNA-Seq analysis were subsequently validated using qRT-PCR. Differentially expressed (2.0-fold change, adj. p < 0.05) transcripts were subjected to ingenuity pathway analysis, which revealed a myriad of affected signaling pathways and functional categories. Activation of ERBB2, FOXM1, ESR1, and IGFBP2 mechanistic networks was most prominent in BC tissue. Additionally, BC tissue exhibited marked enrichment in genes promoting cellular proliferation, migration, survival, and DNA replication and repair. The presence of genes indicative of immune cell infiltration and activation was also observed in BC tissue. We observed high concordance [43.5% (upregulated) and 62.1% (downregulated)] between differentially expressed genes in our study group and those reported for the TCGA BC cohort. Our data provide novel insight on BC biology and suggest common altered molecular networks in BC in this geographic region. Our data suggest future development of therapeutic interventions targeting those common signaling pathways.

The epigenetic mechanisms promoting lineage‐specific commitment of human skeletal (mesenchymal or stromal) stem cells (hMSCs) into adipocytes or osteoblasts are still not fully understood. Herein, we performed an epigenetic library functional screen and identified several novel compounds, including abexinostat, which promoted adipocytic and osteoblastic differentiation of hMSCs. Using gene expression microarrays, chromatin immunoprecipitation for H3K9Ac combined with high‐throughput DNA sequencing (ChIP‐seq), and bioinformatics, we identified several key genes involved in regulating stem cell proliferation and differentiation that were targeted by abexinostat. Concordantly, ChIP‐quantitative polymerase chain reaction revealed marked increase in H3K9Ac epigenetic mark on the promoter region of AdipoQ, FABP4, PPARγ, KLF15, CEBPA, SP7, and ALPL in abexinostat‐treated hMSCs. Pharmacological inhibition of focal adhesion kinase (PF‐573228) or insulin‐like growth factor‐1R/insulin receptor (NVP‐AEW51) signaling exhibited significant inhibition of abexinostat‐mediated adipocytic differentiation, whereas inhibition of WNT (XAV939) or transforming growth factor‐β (SB505124) signaling abrogated abexinostat‐mediated osteogenic differentiation of hMSCs. Our findings provide insight into the understanding of the relationship between the epigenetic effect of histone deacetylase inhibitors, transcription factors, and differentiation pathways governing adipocyte and osteoblast differentiation. Manipulating such pathways allows a novel use for epigenetic compounds in hMSC‐based therapies and tissue engineering. Significance This unbiased epigenetic library functional screen identified several novel compounds, including abexinostat, that promoted adipocytic and osteoblastic differentiation of human skeletal (mesenchymal or stromal) stem cells (hMSCs). These data provide new insight into the understanding of the relationship between the epigenetic effect of histone deacetylase inhibitors, transcription factors, and differentiation pathways controlling adipocyte and osteoblast differentiation of hMSCs. Manipulating such pathways allows a novel use for epigenetic compounds in hMSC‐based therapies for tissue engineering, bone disease, obesity, and metabolic‐disorders.

Cumulative data link cytokine storms with coronavirus disease 2019 (COVID-19) severity. The precise identification of immune cell subsets in bronchoalveolar lavage (BAL) and their correlation with COVID-19 disease severity are currently being unraveled. Herein, we employed iterative clustering and guide-gene selection 2 (ICGS2) as well as uniform manifold approximation and projection (UMAP) dimensionality reduction computational algorithms to decipher the complex immune and cellular composition of BAL, using publicly available datasets from a total of 68,873 single cells derived from two healthy subjects, three patients with mild COVID-19, and five patients with severe COVID-19. Our analysis revealed the presence of neutrophils and macrophage cluster-1 as a hallmark of severe COVID-19. Among the identified gene signatures, IFITM2, IFITM1, H3F3B, SAT1, and S100A8 gene signatures were highly associated with neutrophils, while CCL8, CCL3, CCL2, KLF6, and SPP1 were associated with macrophage cluster-1 in severe-COVID-19 patients. Interestingly, although macrophages were also present in healthy subjects and patients with mild COVID-19, they had different gene signatures, indicative of interstitial and cluster-0 macrophage (i.e., FABP4, APOC1, APOE, C1QB, and NURP1). Additionally, MALAT1, NEAT1, and SNGH25 were downregulated in patients with mild and severe COVID-19. Interferon signaling, FCγ receptor-mediated phagocytosis, IL17, and Tec kinase canonical pathways were enriched in patients with severe COVID-19, while PD-1 and PDL-1 pathways were suppressed. A number of upstream regulators (IFNG, PRL, TLR7, PRL, TGM2, TLR9, IL1B, TNF, NFkB, IL1A, STAT3, CCL5, and others) were also enriched in BAL cells from severe COVID-19-affected patients compared to those from patients with mild COVID-19. Further analyses revealed genes associated with the inflammatory response and chemotaxis of myeloid cells, phagocytes, and granulocytes, among the top activated functional categories in BAL from severe COVID-19-affected patients. Transcriptome data from another cohort of COVID-19-derived peripheral blood mononuclear cells (PBMCs) revealed the presence of several genes common to those found in BAL from patients with severe and mild COVID-19 (IFI27, IFITM3, IFI6, IFIT3, MX1, IFIT1, OASL, IFI30, OAS1) or to those seen only in BAL from severe-COVID-19 patients (S100A8, IFI44, IFI44L, CXCL8, CCR1, PLSCR1, EPSTI1, FPR1, OAS2, OAS3, IL1RN, TYMP, BCL2A1). Taken together, our data reveal the presence of neutrophils and macrophage cluster-1 as the main immune cell subsets associated with severe COVID-19 and identify their inflammatory and chemotactic gene signatures, also partially reflected systemically in the circulation, for possible diagnostic and therapeutic interventions.

The continuous and rapid emergence of new viral strains calls for a better understanding of the fundamental changes occurring within the host cell upon viral infection. In this study, we analyzed RNA-seq transcriptome data from Calu-3 human lung epithelial cells infected with severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) compared to five other viruses namely, severe acute respiratory syndrome coronavirus (SARS-CoV), Middle East Respiratory Syndrome (SARS-MERS), influenzavirus A (FLUA), influenzavirus B (FLUB), and rhinovirus (RHINO) compared to mock-infected cells and characterized their coding and noncoding RNA transcriptional portraits. The induction of interferon, inflammatory, and immune response was a hallmark of SARS-CoV-2 infection. Comprehensive bioinformatics revealed the activation of immune response and defense response to the virus as a common feature of viral infection. Interestingly however, the degree of functional categories and signaling pathways activation varied among different viruses. Ingenuity pathways analysis highlighted altered conical and casual pathways related to TNF, IL1A, and TLR7, which are seen more predominantly during SARS-CoV-2 infection. Nonetheless, the activation of chemotaxis and lipid synthesis was prominent in SARS-CoV-2-infected cells. Despite the commonality among all viruses, our data revealed the hyperactivation of chemotaxis and immune cell trafficking as well as the enhanced fatty acid synthesis as plausible mechanisms that could explain the inflammatory cytokine storms associated with severe cases of COVID-19 and the rapid spread of the virus, respectively.

Background. Hedgehog (Hh) signaling is essential for osteoblast differentiation of mesenchymal progenitors during endochondral bone formation. However, the critical role of Hh signaling during adult bone remodeling remains to be elucidated. Methods. A Smoothened (SMO) antagonist/Hedgehog inhibitor, BMS-833923, identified during a functional screening of a stem cell signaling small molecule library, was investigated for its effects on the osteoblast differentiation of human skeletal (mesenchymal) stem cells (hMSC). Alkaline phosphatase (ALP) activity and Alizarin red staining were employed as markers for osteoblast differentiation and in vitro mineralization capacity, respectively. Global gene expression profiling was performed using the Agilent® microarray platform. Effects on in vivo ectopic bone formation were assessed by implanting hMSC mixed with hydroxyapatite-tricalcium phosphate granules subcutaneously in 8-week-old female nude mice, and the amount of bone formed was assessed using quantitative histology. Results. BMS-833923, a SMO antagonist/Hedgehog inhibitor, exhibited significant inhibitory effects on osteoblast differentiation of hMSCs reflected by decreased ALP activity, in vitro mineralization, and downregulation of osteoblast-related gene expression. Similarly, we observed decreased in vivo ectopic bone formation. Global gene expression profiling of BMS-833923-treated compared to vehicle-treated control cells, identified 348 upregulated and 540 downregulated genes with significant effects on multiple signaling pathways, including GPCR, endochondral ossification, RANK-RANKL, insulin, TNF alpha, IL6, and inflammatory response. Further bioinformatic analysis employing Ingenuity Pathway Analysis revealed significant enrichment in BMS-833923-treated cells for a number of functional categories and networks involved in connective and skeletal tissue development and disorders, e.g., NFκB and STAT signaling. Conclusions. We identified SMO/Hedgehog antagonist (BMS-833923) as a powerful inhibitor of osteoblastic differentiation of hMSC that may be useful as a therapeutic option for treating conditions associated with high heterotopic bone formation and mineralization.

Background. Chemical biology approaches using small molecule inhibitors targeting specific signaling pathways are useful tools to dissect the molecular mechanisms governing stem cell differentiation and for their possible use in therapeutic interventions. Methods. Stem cell signaling small molecule library functional screen was performed employing human bone marrow skeletal (mesenchymal) stem cells (hBMSCs). Alkaline phosphatase (ALP) activity and formation of mineralized matrix visualized by Alizarin red staining were employed as markers for osteoblastic differentiation. Global gene expression profiling was conducted using the Agilent microarray platform, and data normalization and bioinformatics were performed using GeneSpring software. Pathway analyses were conducted using the Ingenuity Pathway Analysis (IPA) tool. In vivo ectopic bone formation was performed using hBMSC mixed with hydroxyapatite–tricalcium phosphate granules that were implanted subcutaneously in 8-week-old female nude mice. Hematoxylin and eosin staining and Sirius red staining were performed to identify bone formation in vivo. Results. Among the tested molecules, LY411575, a potent γ-secretase and Notch signaling inhibitor, exhibited significant inhibitory effects on osteoblastic differentiation of hBMSCs manifested by reduced ALP activity, mineralized matrix formation, and decreased osteoblast-specific gene expression as well as in vivo ectopic bone formation. Global gene expression profiling of LY411575-treated cells revealed changes in multiple signaling pathways, including focal adhesion, insulin, TGFβ, IL6, and Notch signaling, and decreased the expression of genes associated with functional categories of tissue development. Among the affected signaling networks were TGFβ1, SPP1, and ERK regulatory networks. Conclusions. We identified γ-secretase inhibitor (LY411575) as a potent regulator of osteoblastic differentiation of hBMSC that may be useful as a therapeutic option for treating conditions associated with ectopic bone formation.

T cell immunoglobulin mucin-3 (TIM-3) is an immune checkpoint identified as one of the key players in regulating T-cell responses. Studies have shown that TIM-3 is upregulated in the tumor microenvironment (TME). However, the precise role of TIM-3 in colorectal cancer (CRC) TME is yet to be elucidated. We performed phenotypic and molecular characterization of TIM-3+ T cells in the TME and circulation of CRC patients by analyzing tumor tissues (TT, TILs), normal tissues (NT, NILs), and peripheral blood mononuclear cells (PBMC). TIM-3 was upregulated on both CD4+ and CD3+CD4− (CD8+) TILs. CD4+TIM-3+ TILs expressed higher levels of T regulatory cell (Tregs)-signature genes, including FoxP3 and Helios, compared with their TIM-3− counterparts. Transcriptomic and ingenuity pathway analyses showed that TIM-3 potentially activates inflammatory and tumor metastatic pathways. Moreover, NF-κB-mediated transcription factors were upregulated in CD4+TIM-3+ TILs, which could favor proliferation/invasion and induce inflammatory and T-cell exhaustion pathways. In addition, we found that CD4+TIM-3+ TILs potentially support tumor invasion and metastasis, compared with conventional CD4+CD25+ Tregs in the CRC TME. However, functional studies are warranted to support these findings. In conclusion, this study discloses some of the functional pathways of TIM-3+ TILs, which could improve their targeting in more specific therapeutic approaches in CRC patients.

Bone marrow mesenchymal stem cells (BMSCs) are adult multipotent stem cells that can differentiate into mesodermal lineage cells, including adipocytes and osteoblasts. However, the epigenetic mechanisms governing the lineage-specific commitment of BMSCs into adipocytes or osteoblasts are under investigation. Herein, we investigated the epigenetic effect of romidepsin, a small molecule dual inhibitor targeting HDAC1 and HDAC2 identified through an epigenetic library functional screen. BMSCs exposed to romidepsin (5 nM) exhibited enhanced adipocytic and osteoblastic differentiation. Global gene expression and signaling pathway analyses of differentially expressed genes revealed a strong enrichment of genes involved in adipogenesis and osteogenesis in romidepsin-treated BMSCs during induction into adipocytes or osteoblasts, respectively. Pharmacological inhibition of FAK signaling during adipogenesis or inhibition of FAK or TGFβ signaling during osteogenesis diminished the biological effects of romidepsin on BMSCs. The results of chromatin immunoprecipitation combined with quantitative polymerase chain reaction indicated a significant increase in H3K9Ac epigenetic markers in the promoter regions of peroxisome proliferator-activated receptor gamma (PPARγ) and KLF15 (related to adipogenesis) or SP7 (Osterix) and alkaline phosphatase (ALP) (related to osteogenesis) in romidepsin-treated BMSCs. Our data indicated that romidepsin is a novel in vitro modulator of adipocytic and osteoblastic differentiation of BMSCs.

Several studies have demonstrated the multipotentiality of human neonatal foreskin stromal cells (hNSSCs) as being able to differentiate into adipocytes and osteoblasts and potentially other cell types. Recently, we demonstrated that hNSSCs play a role during in vitro angiogenesis and appear to possess a capacity to differentiate into endothelial-like cells; however, their angiogenic potential within an ex vivo environment remains unclear. Current study shows hNSSCs to display significant migration potential in the undifferentiated state and high responsiveness in the in vitro wound healing scratch assay. When hNSSCs were seeded onto the top of the CAM, human von Willebrand factor (hVWF), CD31, smooth muscle actin (SMA), and factor XIIIa positive cells were observed in the chick endothelium. CAMs transplanted with endothelial-differentiated hNSSCs displayed a higher number of blood vessels containing hNSSCs compared to CAMs transplanted with undifferentiated hNSSCs. Interestingly, undifferentiated hNSSCs showed a propensity to differentiate towards ectoderm with indication of epidermal formation with cells positive for CD1a, CK5/6, CK19, FXIIIa, and S-100 cells, which warrant further investigation. Our findings imply a potential angiogenic role for hNSSCs ex vivo in the differentiated and undifferentiated state, with potential contribution to blood vessel formation and potential application in tissue regeneration and vascularization.