Reciprocal interplays between MicroRNAs and pluripotency transcription factors in dictating stemness features in human cancers.

The interplay between microRNAs (miRNAs) and pluripotency transcription factors (TFs) orchestrates the acquisition of cancer stem cell (CSC) features during the course of malignant transformation, rendering them essential cancer cell dependencies and therapeutic vulnerabilities. In this review, we discuss emerging themes in tumor heterogeneity, including the clonal evolution and the CSC models and their implications in resistance to cancer therapies, and then provide thorough coverage on the roles played by key TFs in maintaining normal and malignant stem cell pluripotency and plasticity. In addition, we discuss the reciprocal interactions between miRNAs and MYC, OCT4, NANOG, SOX2, and KLF4 pluripotency TFs and their contributions to tumorigenesis. We provide our view on the potential to interfere with key miRNA-TF networks through the use of RNA-based therapeutics as single agents or in combination with other therapeutic strategies, to abrogate the CSC state and render tumor cells more responsive to standard and targeted therapies.

Long non-coding RNA and RNA-binding protein interactions in cancer: Experimental and machine learning approaches.

Understanding the complex and specific roles played by non-coding RNAs (ncRNAs), which comprise the bulk of the genome, is important for understanding virtually every hallmark of cancer. This large group of molecules plays pivotal roles in key regulatory mechanisms in various cellular processes. Regulatory mechanisms, mediated by long non-coding RNA (lncRNA) and RNA-binding protein (RBP) interactions, are well documented in several types of cancer. Their effects are enabled through networks affecting lncRNA and RBP stability, RNA metabolism including N6-methyladenosine (m6A) and alternative splicing, subcellular localization, and numerous other mechanisms involved in cancer. In this review, we discuss the reciprocal interplay between lncRNAs and RBPs and their involvement in epigenetic regulation via histone modifications, as well as their key role in resistance to cancer therapy. Other aspects of RBPs including their structural domains, provide a deeper knowledge on how lncRNAs and RBPs interact and exert their biological functions. In addition, current state-of-the-art knowledge, facilitated by machine and deep learning approaches, unravels such interactions in better details to further enhance our understanding of the field, and the potential to harness RNA-based therapeutics as an alternative treatment modality for cancer are discussed.

Therapeutic targeting of the TPX2/TTK network in colorectal cancer.

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.

Identification of a Gene Panel Predictive of Triple-Negative Breast Cancer Response to Neoadjuvant Chemotherapy Employing Transcriptomic and Functional Validation.

Triple-negative breast cancer (TNBC) patients exhibiting pathological complete response (pCR) have better clinical outcomes compared to those with residual disease (RD). Therefore, robust biomarkers that can predict pCR may help with triage and resource prioritization in patients with TNBC. Herein, we identified a gene panel predictive of RD and pCR in TNBC from the discovery (n = 90) treatment-naive tumor transcriptomic data. Eight RD-derived genes were identified as TNBC-essential genes, which were highly predicative of overall survival (OS) and relapse-free survival (RFS) in an additional cohort of basal breast cancer (n = 442). Mechanistically, targeted depletion of the eight genes reduced the proliferation potential of TNBC cell models, while most remarkable effects were for combined SLC39A7, TIMM13, BANF1, and MVD knockdown in conjunction with doxorubicin. Orthogonal partial least squares-discriminant analysis (OPLS-DA) and receiver operating characteristic curve (ROC) analyses revealed significant predictive power for the identified gene panels with an area under the curve (AUC) of 0.75 for the validation cohort (n = 50) to discriminate RD from pCR. Protein-Protein Interaction (PPI) network analysis of the pCR-derived gene signature identified an 87-immune gene signature highly predictive of pCR, which correlated with better OS, RFS, and distant-metastasis-free survival (DMFS) in an independent cohort of basal and, to a lesser extent, HER2+ breast cancer. Our data have identified gene signatures predicative of RD and pCR in TNBC with potential clinical implications.

Noncoding RNAs as potential mediators of resistance to cancer immunotherapy.

Substantial evolution in cancer therapy has been witnessed lately, steering mainly towards immunotherapeutic approaches, replacing or in combination with classical therapies. Whereas the use of various immunotherapy approaches, such as adoptive T cell therapy, genetically-modified T cells, or immune checkpoint inhibitors, has been a triumph for cancer immunotherapy, the great challenge is the ability of the immune system to sustain long lasting anti-tumor response. Additionally, epigenetic changes in a suppressive tumor microenvironment can pertain to T cell exhaustion, limiting their functionality. Noncoding RNAs (ncRNAs) have emerged over the last years as key players in epigenetic regulation. Among those, microRNAs (miRNAs) and long noncoding RNAs (lncRNAs) have been studied extensively for their potential role in regulating tumor immunity through direct regulation of genes involved in immune activation or suppression. In this review, we will provide an overview of contemporary approaches for cancer immunotherapy and will present the current state of knowledge implicating miRNAs and lncRNAs in regulating immune response against human cancer and their potential implications in resistance to cancer immunotherapy, with main emphasis on immune checkpoints regulation.

SOX4: Epigenetic regulation and role in tumorigenesis.

Sex-determining region Y-related (SRY) high-mobility group box 4 (SOX4) is a member of the group C subfamily of SOX transcription factors and promotes tumorigenesis by endowing cancer cells with survival, migratory, and invasive capacities. Emerging evidence has highlighted an unequivocal role for this transcription factor in mediating various signaling pathways involved in tumorigenesis, epithelial-to-mesenchymal transition (EMT), and tumor progression. During the last decade, numerous studies have highlighted the epigenetic interplay between SOX4-targeting microRNAs (miRNAs), long noncoding RNAs (lncRNAs) and SOX4 and the subsequent modulation of tumorigenesis, invasion and metastasis. In this review, we summarize the current state of knowledge about the role of SOX4 in cancer development and progression, the epigenetic regulation of SOX4, and the potential utilization of SOX4 as a diagnostic and prognostic biomarker and its depletion as a therapeutic target.

Transcriptional alterations of protein coding and noncoding RNAs in triple negative breast cancer in response to DNA methyltransferases inhibition.

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.

Bone morphogenetic protein 2 (BMP2) induces growth suppression and enhances chemosensitivity of human colon cancer cells.

BACKGROUND: Molecular profiling of colorectal cancer (CRC) based on global gene expression has revealed multiple dysregulated signalling pathways associated with drug resistance and poor prognosis. However, the role of BMP2 signaling in CRC is not fully characterised. METHODS: Bioinformatics data analysis were conducted on the GSE21510 dataset. Leniviral technology was utilized to stably express BMP2 in the HCT116 CRC model. Gene expression profiling was conducted using Agilent microarray platform while data normalization and bioinformatics were conducted using GeneSpring software. Changes in gene expression were assessed using qRT-PCR. AlamarBlue assay was used to assess cell viability in vitro. In vivo experiments were conducted using SCID mice. RESULTS: Our data revealed frequent downregulation of BMP2 in primary CRC tissues. Additionally, interrogation of publically available gene expression datasets revealed significant downregulation of BMP2 in metastatic recurrent compared to non-metastatic cancer (p = 0.02). Global gene expression analysis in CRC cells over-expressing BMP2 revealed multiple dysregulated pathways mostly affecting cell cycle and DNA damage response. Concordantly, lentiviral-mediated re-expression of BMP2 inhibited HCT116 CRC growth, sphere formation, clonogenic potential, cell migration, and sensitized CRC cells to 5-fluorouracil (5-FU) in vitro. Additionally, BMP2 inhibited CRC tumor formation in SCID mice. CONCLUSIONS: Our data revealed an inhibitory role for BMP2 in CRC, suggesting that restoration of BMP2 expression could be a potential therapeutic strategy for CRC.

Apigenin and Rutaecarpine reduce the burden of cellular senescence in bone marrow stromal stem cells.

Introduction: Osteoporosis 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 results: Apigenin 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. Discussion: Our 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.

Long non-coding RNA AC099850.4 correlates with advanced disease state and predicts worse prognosis in triple-negative breast cancer.

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.

Single-Cell Transcriptome Analysis Revealed Heterogeneity and Identified Novel Therapeutic Targets for Breast Cancer Subtypes.

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.

Transcriptome profiling and network enrichment analyses identify subtype-specific therapeutic gene targets for breast cancer and their microRNA regulatory networks.

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.

In vitro differentiation of human skin-derived multipotent stromal cells into putative endothelial-like cells.

BACKGROUND: Multipotent stem cells have been successfully isolated from various tissues and are currently utilized for tissue-engineering and cell-based therapies. Among the many sources, skin has recently emerged as an attractive source for multipotent cells because of its abundance. Recent literature showed that skin stromal cells (SSCs) possess mesoderm lineage differentiation potential; however, the endothelial differentiation and angiogenic potential of SSC remains elusive. In our study, SSCs were isolated from human neonatal foreskin (hNFSSCs) and adult dermal skin (hADSSCs) using explants cultures and were compared with bone marrow (hMSC-TERT) and adipose tissue-derived mesenchymal stem cells (hADMSCs) for their potential differentiation into osteoblasts, adipocytes, and endothelial cells. RESULTS: Concordant with previous studies, both MSCs and SSCs showed similar morphology, surface protein expression, and were able to differentiate into osteoblasts and adipocytes. Using an endothelial induction culture system combined with an in vitro matrigel angiogenesis assay, hNFSSCs and hADSSCs exhibited the highest tube-forming capability, which was similar to those formed by human umbilical vein endothelial cells (HUVEC), with hNFSSCs forming the most tightly packed, longest, and largest diameter tubules among the three cell types. CD146 was highly expressed on hNFSSCs and HUVEC followed by hADSSCs, and hMSC-TERT, while its expression was almost absent on hADMSCs. Similarly, higher vascular density (based on the expression of CD31, CD34, vWF, CD146 and SMA) was observed in neonatal skin, followed by adult dermal skin and adipose tissue. Thus, our preliminary data indicated a plausible relationship between vascular densities, and the expression of CD146 on multipotent cells derived from those tissues. CONCLUSIONS: Our data is the first to demonstrate that human dermal skin stromal cells can be differentiated into endothelial lineage. Hence, SSCs represents a novel source of stem/stromal cells for tissue regeneration and the vascularization of engineered tissues. Moreover, the CD146 investigations suggested that the microenvironmental niche might contribute to direct stromal cells multipotency toward certain lineages, which warrants further investigation.

Comparative investigation of the differentiation capability of bone-marrow- and adipose-derived mesenchymal stem cells by qualitative and quantitative analysis.

Mesenchymal stem cells (MSCs) hold promise for cell-based therapy in regenerative medicine. To date, MSCs have been obtained from conventional bone marrow via a highly invasive procedure. Therefore, MSCs are now also isolated from sources such as adipose tissue, cord blood and cord stroma, a subject of growing interest. As the characterization and differentiation potential of adipose-derived MSCs (AD-MSCs) and bone-marrow-derived MSCs (BM-MSCs) have not been documented, we have evaluated and compared the characteristics of both MSC types by qualitative and quantitative analyses. Both cell types show similar morphology and surface protein expression, being positive for stromal-associated markers and negative for hematopoietic and endothelial markers. The colony-forming potential of AD-MSCs is distinctly higher than that of BM-MSCs. Nonetheless, similar adipogenic and osteogenic differentiation is observed in both groups of MSCs. Cytochemical qualitative analysis and calcium mineralization demonstrate higher levels toward osteogenic differentiation in BM-MSCs than in AD-MSCs. On the contrary, the percentage of Nile red oil staining for differentiated adipocytes is higher in AD-MSCs than in BM-MSCs. Quantitative real-time polymerase chain reaction shows similar patterns of osteogenic- and adipogenic-associated gene expression in both cell types. Each of theMSCs respond in functional analysis by exhibiting unique properties at the differentiation level according to their micro-environmental niche. Thus, quantitative analysis might be a valuable means of describing stem cell multipotency, in addition to qualitative investigation.

Skin-derived multipotent stromal cells--an archrival for mesenchymal stem cells.

Progenitor stem cells have been identified, isolated and characterized in numerous tissues and organs. However, their therapeutic potential and the use of these stem cells remain elusive except for a few progenitor cells from bone marrow, umbilical cord blood, eyes and dental pulp. The use of bone marrow-derived hematopoietic stem cells (HSC) or mesenchymal stem cells (MSCs) is restricted due to their extreme invasive procedures, low differentiation potential with age and rejection. Thus, we need a clinical grade alternative to progenitor stem cells with a high potential to differentiate, naïve and is relatively easy in in vitro propagation. In this review, we summarize cell populations of adherent and floating spheres derived from different origins of skin, or correctly foreskin, by enzymatic digestion compared with established MSCs. The morphology, phenotype, differentiation capability and immunosuppressive property of the adherent cell populations are comparable with MSCs. Serum-free cultured floating spheres have limited mesodermal but higher neurogenic differentation potential, analogous to neural crest stem cells. Both the populations confirmed their plethora potential in in vitro. Together, it may be noted that the skin-derived adherent cell populations and floating cells can be good alternative sources of progenitor cells especially in cosmetic, plastic and sports regenerative medicine.

LncRNA-Based Classification of Triple Negative Breast Cancer Revealed Inherent Tumor Heterogeneity and Vulnerabilities.

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.

Transcriptional landscape associated with TNBC resistance to neoadjuvant chemotherapy revealed by single-cell RNA-seq.

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.

Epigenetic regulation of triple negative breast cancer (TNBC) by TGF-ß signaling.

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.

Molecular subtyping and functional validation of TTK, TPX2, UBE2C, and LRP8 in sensitivity of TNBC to paclitaxel.

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.

Single-cell long noncoding RNA (lncRNA) transcriptome implicates MALAT1 in triple-negative breast cancer (TNBC) resistance to neoadjuvant chemotherapy.

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.