The Transcriptional Landscape of Coding and Noncoding RNAs in Recurrent and Nonrecurrent Colon Cancer.

A considerable number of colon cancer patients with local or local advanced disease suffer from recurrence and there is an urgent need for better prognostic biomarkers in this setting. Here, the transcriptomic landscape of mRNAs, long noncoding RNAs, snRNAs, small nucleolar RNAs (snoRNAs), small Cajal body-specific RNAs, pseudogenes, and circular RNAs, as well as RNAs denoted as miscellaneous RNAs, was profiled by total RNA sequencing. In addition to well-known coding and noncoding RNAs, differential expression analysis also uncovered transcripts that have not been implicated previously in colon cancer, such as RNA5SP149, RNU4-2, and SNORD3A. Moreover, there was a profound global up-regulation of snRNA pseudogenes, snoRNAs, and rRNA pseudogenes in more advanced tumors. A global down-regulation of circular RNAs in tumors relative to normal tissues was observed, although only a few were expressed differentially between tumor stages. Many previously undescribed transcripts, including RNU6-620P, RNU2-20P, VTRNA1-3, and RNA5SP60, indicated strong prognostic biomarker potential in receiver operating characteristics analyses. In summary, this study unveiled numerous differentially expressed RNAs across various classes between recurrent and nonrecurrent colon cancer. Notably, there was a significant global up-regulation of snRNA pseudogenes, snoRNAs, and rRNA pseudogenes in advanced tumors. Many of these newly discovered candidates demonstrated a strong prognostic potential for stage II colon cancer.

Spatial Profiling of Circular RNAs in Cancer Reveals High Expression in Muscle and Stromal Cells.

Circular RNAs (circRNA) are covalently closed molecules that can play important roles in cancer development and progression. Hundreds of differentially expressed circRNAs between tumors and adjacent normal tissues have been identified in studies using RNA sequencing or microarrays, emphasizing a strong translational potential. Most previous studies have been performed using RNA from bulk tissues and lack information on the spatial expression patterns of circRNAs. Here, we showed that the majority of differentially expressed circRNAs from bulk tissue analyses of colon tumors relative to adjacent normal tissues were surprisingly not differentially expressed when comparing cancer cells directly with normal epithelial cells. Manipulating the proliferation rates of cells grown in culture revealed that these discrepancies were explained by circRNAs accumulating to high levels in quiescent muscle cells due to their high stability; on the contrary, circRNAs were diluted to low levels in the fast-proliferating cancer cells due to their slow biogenesis rates. Thus, different subcompartments of colon tumors and adjacent normal tissues exhibited striking differences in circRNA expression patterns. Likewise, the high circRNA content in muscle cells was also a strong confounding factor in bulk analyses of circRNAs in bladder and prostate cancers. Together, these findings emphasize the limitations of using bulk tissues for studying differential circRNA expression in cancer and highlight a particular need for spatial analysis in this field of research. SIGNIFICANCE: The abundance of circRNAs varies systematically between subcompartments of solid tumors and adjacent tissues, implying that differentially expressed circRNAs discovered in bulk tissue analyses may reflect differences in cell type composition between samples.

Impact of U2AF1 mutations on circular RNA expression in myelodysplastic neoplasms.

Mutations in U2AF1 are relatively common in myelodysplastic neoplasms (MDS) and are associated with an inferior prognosis, but the molecular mechanisms underlying this are not fully elucidated. Circular RNAs (circRNAs) have been implicated in cancer, but it is unknown how mutations in splicing factors may impact on circRNA biogenesis. Here, we used RNA-sequencing to investigate the effects of U2AF1 mutations on circRNA expression in K562 cells with a doxycycline-inducible U2AF1S34 mutation, in a mouse model with a doxycycline-inducible U2AF1S34 mutation, and in FACS-sorted CD34+ bone marrow cells from MDS patients with either U2AF1S34 or U2AF1Q157 mutations. In all contexts, we found an increase in global circRNA levels in the U2AF1-mutated setting, which was independent of expression changes in the cognate linear host genes. In patients, the U2AF1S34 and U2AF1Q157 mutations were both associated with an overall increased expression of circRNAs. circRNAs generated by a non-Alu-mediated mechanism generally showed the largest increase in expression levels. Several well-described cancer-associated circRNAs, including circZNF609 and circCSNK1G3, were upregulated in MDS patients with U2AF1 mutations compared to U2AF1-wildtype MDS controls. In conclusion, high circRNA expression is observed in association with U2AF1 mutations in three biological systems, presenting an interesting possibility for biomarker and therapeutic investigation.

Tacedinaline (CI-994), a class I HDAC inhibitor, targets intrinsic tumor growth and leptomeningeal dissemination in MYC-driven medulloblastoma while making them susceptible to anti-CD47-induced macrophage phagocytosis via NF-kB-TGM2 driven tumor inflammation.

BACKGROUND: While major advances have been made in improving the quality of life and survival of children with most forms of medulloblastoma (MB), those with MYC-driven tumors (Grp3-MB) still suffer significant morbidity and mortality. There is an urgent need to explore multimodal therapeutic regimens which are effective and safe for children. Large-scale studies have revealed abnormal cancer epigenomes caused by mutations and structural alterations of chromatin modifiers, aberrant DNA methylation, and histone modification signatures. Therefore, targeting epigenetic modifiers for cancer treatment has gained increasing interest, and inhibitors for various epigenetic modulators have been intensively studied in clinical trials. Here, we report a cross-entity, epigenetic drug screen to evaluate therapeutic vulnerabilities in MYC amplified MB, which sensitizes them to macrophage-mediated phagocytosis by targeting the CD47-signal regulatory protein α (SIRPα) innate checkpoint pathway. METHODS: We performed a primary screen including 78 epigenetic inhibitors and a secondary screen including 20 histone deacetylase inhibitors (HDACi) to compare response profiles in atypical teratoid/rhabdoid tumor (AT/RT, n=11), MB (n=14), and glioblastoma (n=14). This unbiased approach revealed the preferential activity of HDACi in MYC-driven MB. Importantly, the class I selective HDACi, CI-994, showed significant cell viability reduction mediated by induction of apoptosis in MYC-driven MB, with little-to-no activity in non-MYC-driven MB, AT/RT, and glioblastoma in vitro. We tested the combinatorial effect of targeting class I HDACs and the CD47-SIRPa phagocytosis checkpoint pathway using in vitro phagocytosis assays and in vivo orthotopic xenograft models. RESULTS: CI-994 displayed antitumoral effects at the primary site and the metastatic compartment in two orthotopic mouse models of MYC-driven MB. Furthermore, RNA sequencing revealed nuclear factor-kB (NF-κB) pathway induction as a response to CI-994 treatment, followed by transglutaminase 2 (TGM2) expression, which enhanced inflammatory cytokine secretion. We further show interferon-γ release and cell surface expression of engulfment ('eat-me') signals (such as calreticulin). Finally, combining CI-994 treatment with an anti-CD47 mAb targeting the CD47-SIRPα phagocytosis checkpoint enhanced in vitro phagocytosis and survival in tumor-bearing mice. CONCLUSION: Together, these findings suggest a dynamic relationship between MYC amplification and innate immune suppression in MYC amplified MB and support further investigation of phagocytosis modulation as a strategy to enhance cancer immunotherapy responses.

VEGFA-targeting miR-agshRNAs combine efficacy with specificity and safety for retinal gene therapy.

Retinal gene therapy using RNA interference (RNAi) to silence targeted genes requires both efficacy and safety. Short hairpin RNAs (shRNAs) are useful for RNAi, but high expression levels and activity from the co-delivered passenger strand may cause undesirable cellular responses. Ago2-dependent shRNAs (agshRNAs) produce no passenger strand activity. To enhance efficacy and to investigate improvements in safety, we have generated VEGFA-targeting agshRNAs and microRNA (miRNA)-embedded agshRNAs (miR-agshRNAs) and inserted these RNAi effectors in Pol II/III-driven expression cassettes and lentiviral vectors (LVs). Compared with corresponding shRNAs, agshRNAs and miR-agshRNAs increased specificity and safety, while retaining a high knockdown efficacy and abolishing passenger strand activity. The agshRNAs also caused significantly smaller reductions in cell viability and reduced competition with the processing of endogenous miR21 compared with their shRNA counterparts. RNA sequencing (RNA-seq) analysis of LV-transduced ARPE19 cells revealed that expression of shRNAs in general leads to more changes in gene expression levels compared with their agshRNA counterparts and activation of immune-related pathways. In mice, subretinal delivery of LVs encoding tissue-specific miR-agshRNAs resulted in retinal pigment epithelium (RPE)-restricted expression and significant knockdown of Vegfa in transduced RPE cells. Collectively, our data suggest that agshRNAs and miR-agshRNA possess important advantages over shRNAs, thereby posing a clinically relevant approach with respect to efficacy, specificity, and safety.

TLR2 and TLR7 mediate distinct immunopathological and antiviral plasmacytoid dendritic cell responses to SARS-CoV-2 infection.

Understanding the molecular pathways driving the acute antiviral and inflammatory response to SARS-CoV-2 infection is critical for developing treatments for severe COVID-19. Here, we find decreasing number of circulating plasmacytoid dendritic cells (pDCs) in COVID-19 patients early after symptom onset, correlating with disease severity. pDC depletion is transient and coincides with decreased expression of antiviral type I IFNα and of systemic inflammatory cytokines CXCL10 and IL-6. Using an in vitro stem cell-based human pDC model, we further demonstrate that pDCs, while not supporting SARS-CoV-2 replication, directly sense the virus and in response produce multiple antiviral (interferons: IFNα and IFNλ1) and inflammatory (IL-6, IL-8, CXCL10) cytokines that protect epithelial cells from de novo SARS-CoV-2 infection. Via targeted deletion of virus-recognition innate immune pathways, we identify TLR7-MyD88 signaling as crucial for production of antiviral interferons (IFNs), whereas Toll-like receptor (TLR)2 is responsible for the inflammatory IL-6 response. We further show that SARS-CoV-2 engages the receptor neuropilin-1 on pDCs to selectively mitigate the antiviral interferon response, but not the IL-6 response, suggesting neuropilin-1 as potential therapeutic target for stimulation of TLR7-mediated antiviral protection.

Transcriptomics Profiling Identifies Cisplatin-Inducible Death Receptor 5 Antisense Long Non-coding RNA as a Modulator of Proliferation and Metastasis in HeLa Cells.

Cisplatin is a well-known cancer chemotherapeutic agent but how extensively long non-coding RNA (lncRNA) expression is modulated by cisplatin is unknown. It is imperative to employ a comprehensive approach to obtain a better account of cisplatin-mediated changes in the expression of lncRNAs. In this study, we used a transcriptomics approach to profile lncRNAs in cisplatin-treated HeLa cells, which resulted in identification of 10,214 differentially expressed lncRNAs, of which 2,500 were antisense lncRNAs. For functional analyses, we knocked down one of the cisplatin inducible lncRNAs, death receptor 5 antisense (DR5-AS) lncRNA, which resulted in a morphological change in HeLa cell shape without inducing any cell death. A second round of transcriptomics-based profiling revealed differential expression of genes associated with immune system, motility and cell cycle in DR5-AS knockdown HeLa cells. Cellular analyses showed that DR5-AS reduced cell proliferation and caused a cell cycle arrest at S and G2/M phases. Moreover, DR5-AS knockdown reduced the invasive capacity of HeLa cells in zebrafish xenograft model. These results suggest that cisplatin-mediated pleiotropic effects, such as reduction in cell proliferation, metastasis and cell cycle arrest, may be mediated by lncRNAs.

The long non-coding RNA HOTAIRM1 promotes tumor aggressiveness and radiotherapy resistance in glioblastoma.

Glioblastoma is the most common malignant primary brain tumor. To date, clinically relevant biomarkers are restricted to isocitrate dehydrogenase (IDH) gene 1 or 2 mutations and O6-methylguanine DNA methyltransferase (MGMT) promoter methylation. Long non-coding RNAs (lncRNAs) have been shown to contribute to glioblastoma pathogenesis and could potentially serve as novel biomarkers. The clinical significance of HOXA Transcript Antisense RNA, Myeloid-Specific 1 (HOTAIRM1) was determined by analyzing HOTAIRM1 in multiple glioblastoma gene expression data sets for associations with prognosis, as well as, IDH mutation and MGMT promoter methylation status. Finally, the role of HOTAIRM1 in glioblastoma biology and radiotherapy resistance was characterized in vitro and in vivo. We identified HOTAIRM1 as a candidate lncRNA whose up-regulation is significantly associated with shorter survival of glioblastoma patients, independent from IDH mutation and MGMT promoter methylation. Glioblastoma cell line models uniformly showed reduced cell viability, decreased invasive growth and diminished colony formation capacity upon HOTAIRM1 down-regulation. Integrated proteogenomic analyses revealed impaired mitochondrial function and determination of reactive oxygen species (ROS) levels confirmed increased ROS levels upon HOTAIRM1 knock-down. HOTAIRM1 knock-down decreased expression of transglutaminase 2 (TGM2), a candidate protein implicated in mitochondrial function, and knock-down of TGM2 mimicked the phenotype of HOTAIRM1 down-regulation in glioblastoma cells. Moreover, HOTAIRM1 modulates radiosensitivity of glioblastoma cells both in vitro and in vivo. Our data support a role for HOTAIRM1 as a driver of biological aggressiveness, radioresistance and poor outcome in glioblastoma. Targeting HOTAIRM1 may be a promising new therapeutic approach.

Distinct circular RNA expression profiles in pediatric ependymomas.

Pediatric ependymomas frequently develop in the cerebellum and are currently treated using non-specific therapies, in part, because few somatically mutated driver genes are present, and the underlying pathobiology is poorly described. Circular RNAs (circRNAs) constitute as a large class of primarily non-coding RNAs with important roles in tumorigenesis, but they have not been described in pediatric ependymomas. To advance our molecular understanding of ependymomas, we performed Next Generation Sequencing of rRNA-depleted total RNA of 10 primary ependymoma and three control samples. CircRNA expression patterns were correlated to disease stage, outcome, age, and gender. We found a profound global downregulation of circRNAs in ependymoma relative to control samples. Many differentially expressed circRNAs were discovered and circSMARCA5 and circ-FBXW7, which are described as tumor suppressors in glioma and glioblastomas in adults, were among the most downregulated. Moreover, patients with a dismal outcome clustered separately from patients with a good prognosis in unsupervised hierarchical cluster analyses. Next, NanoString nCounter experiments were performed, using a custom-designed panel targeting 66 selected circRNAs, on a larger cohort that also included medulloblastomas and pilocytic astrocytomas. These experiments indicated that circRNA expression profiles are different among distinct pediatric brain tumor subtypes. In particular, circRNAs derived from RMST, LRBA, WDR78, DRC1 and BBS9 genes were specifically upregulated in ependymomas. In conclusion, circRNAs have different expression profiles in ependymomas relative to controls and between survivors and patients with a dismal outcome, suggesting that circRNAs could be exerted as diagnostic and prognostic biomarkers in the future if further validated in larger cohorts.

The long noncoding RNA TP73-AS1 promotes tumorigenicity of medulloblastoma cells.

Medulloblastoma is the most common malignant brain cancer in children. Since previous studies have mainly focused on alterations in the coding genome, our understanding of the contribution of long noncoding RNAs (lncRNAs) to medulloblastoma biology is just emerging. Using patient-derived data, we show that the promoter of lncRNA TP73-AS1 is hypomethylated and that the transcript is highly expressed in the SHH subgroup. Furthermore, high expression of TP73-AS1 is correlated with poor outcome in patients with TP53 wild-type SHH tumors. Silencing TP73-AS1 in medulloblastoma tumor cells induced apoptosis, while proliferation and migration were inhibited in culture. In vivo, silencing TP73-AS1 in medulloblastoma tumor cells resulted in reduced tumor growth, reduced proliferation of tumor cells, increased apoptosis and led to prolonged survival of tumor-bearing mice. Together, our study suggests that the lncRNA TP73-AS1 is a prognostic marker and therapeutic target in medulloblastoma tumors and serves as a proof of concept that lncRNAs are important factors in the disease.

The lncRNA TP73-AS1 is linked to aggressiveness in glioblastoma and promotes temozolomide resistance in glioblastoma cancer stem cells.

Glioblastoma multiform (GBM) is the most common brain tumor characterized by a dismal prognosis. GBM cancer stem cells (gCSC) or tumor-initiating cells are the cell population within the tumor-driving therapy resistance and recurrence. While temozolomide (TMZ), an alkylating agent, constitutes the first-line chemotherapeutic significantly improving survival in GBM patients, resistance against this compound commonly leads to GBM recurrence and treatment failure. Although the roles of protein-coding transcripts, proteins and microRNA in gCSC, and therapy resistance have been comprehensively investigated, very little is known about the role of long noncoding RNAs (lncRNAs) in this context. Using nonoverlapping, independent RNA sequencing and gene expression profiling datasets, we reveal that TP73-AS1 constitutes a clinically relevant lncRNA in GBM. Specifically, we demonstrate significant overexpression of TP73-AS1 in primary GBM samples, which is particularly increased in the gCSC. More importantly, we demonstrate that TP73-AS1 comprises a prognostic biomarker in glioma and in GBM with high expression identifying patients with particularly poor prognosis. Using CRISPRi to downregulate our candidate lncRNA in gCSC, we demonstrate that TP73-AS1 promotes TMZ resistance in gCSC and is linked to regulation of the expression of metabolism- related genes and ALDH1A1, a protein known to be expressed in cancer stem cell markers and protects gCSC from TMZ treatment. Taken together, our results reveal that high TP73-AS1 predicts poor prognosis in primary GBM cohorts and that this lncRNA promotes tumor aggressiveness and TMZ resistance in gCSC.