Conditional c-MYC activation in catecholaminergic cells drives distinct neuroendocrine tumors: neuroblastoma vs somatostatinoma.

The MYC proto-oncogenes (c-MYC, MYCN , MYCL ) are among the most deregulated oncogenic drivers in human malignancies including high-risk neuroblastoma, 50% of which are MYCN -amplified. Genetically engineered mouse models (GEMMs) based on the MYCN transgene have greatly expanded the understanding of neuroblastoma biology and are powerful tools for testing new therapies. However, a lack of c-MYC-driven GEMMs has hampered the ability to better understand mechanisms of neuroblastoma oncogenesis and therapy development given that c-MYC is also an important driver of many high-risk neuroblastomas. In this study, we report two transgenic murine neuroendocrine models driven by conditional c-MYC induction in tyrosine hydroxylase (Th) and dopamine ß-hydroxylase (Dbh)-expressing cells. c-MYC induction in Th-expressing cells leads to a preponderance of Pdx1 + somatostatinomas, a type of pancreatic neuroendocrine tumor (PNET), resembling human somatostatinoma with highly expressed gene signatures of δ cells and potassium channels. In contrast, c-MYC induction in Dbh-expressing cells leads to onset of neuroblastomas, showing a better transforming capacity than MYCN in a comparable C57BL/6 genetic background. The c-MYC murine neuroblastoma tumors recapitulate the pathologic and genetic features of human neuroblastoma, express GD2, and respond to anti-GD2 immunotherapy. This model also responds to DFMO, an FDA-approved inhibitor targeting ODC1, which is a known MYC transcriptional target. Thus, establishing c-MYC-overexpressing GEMMs resulted in different but related tumor types depending on the targeted cell and provide useful tools for testing immunotherapies and targeted therapies for these diseases.

Duocarmycin SA Reduces Proliferation and Increases Apoptosis in Acute Myeloid Leukemia Cells In Vitro.

Acute myeloid leukemia (AML) is a hematological malignancy that is characterized by an expansion of immature myeloid precursors. Despite therapeutic advances, the prognosis of AML patients remains poor and there is a need for the evaluation of promising therapeutic candidates to treat the disease. The objective of this study was to evaluate the efficacy of duocarmycin Stable A (DSA) in AML cells in vitro. We hypothesized that DSA would induce DNA damage in the form of DNA double-strand breaks (DSBs) and exert cytotoxic effects on AML cells within the picomolar range. Human AML cell lines Molm-14 and HL-60 were used to perform 3-(4,5-dimethylthiazolyl-2)-2,5-diphenyltetrazolium bromide (MTT), DNA DSBs, cell cycle, 5-ethynyl-2-deoxyuridine (EdU), colony formation unit (CFU), Annexin V, RNA sequencing and other assays described in this study. Our results showed that DSA induced DNA DSBs, induced cell cycle arrest at the G2M phase, reduced proliferation and increased apoptosis in AML cells. Additionally, RNA sequencing results showed that DSA regulates genes that are associated with cellular processes such as DNA repair, G2M checkpoint and apoptosis. These results suggest that DSA is efficacious in AML cells and is therefore a promising potential therapeutic candidate that can be further evaluated for the treatment of AML.

RBM39 degrader invigorates natural killer cells to eradicate neuroblastoma despite cancer cell plasticity.

The cellular plasticity of neuroblastoma is defined by a mixture of two major cell states, adrenergic (ADRN) and mesenchymal (MES), which may contribute to therapy resistance. However, how neuroblastoma cells switch cellular states during therapy remains largely unknown and how to eradicate neuroblastoma regardless of their cell states is a clinical challenge. To better understand the lineage switch of neuroblastoma in chemoresistance, we comprehensively defined the transcriptomic and epigenetic map of ADRN and MES types of neuroblastomas using human and murine models treated with indisulam, a selective RBM39 degrader. We showed that cancer cells not only undergo a bidirectional switch between ADRN and MES states, but also acquire additional cellular states, reminiscent of the developmental pliancy of neural crest cells. The lineage alterations are coupled with epigenetic reprogramming and dependency switch of lineage-specific transcription factors, epigenetic modifiers and targetable kinases. Through targeting RNA splicing, indisulam induces an inflammatory tumor microenvironment and enhances anticancer activity of natural killer cells. The combination of indisulam with anti-GD2 immunotherapy results in a durable, complete response in high-risk transgenic neuroblastoma models, providing an innovative, rational therapeutic approach to eradicate tumor cells regardless of their potential to switch cell states.

Metabolic reprogramming of cancer cells by JMJD6-mediated pre-mRNA splicing associated with therapeutic response to splicing inhibitor.

Dysregulated pre-mRNA splicing and metabolism are two hallmarks of MYC-driven cancers. Pharmacological inhibition of both processes has been extensively investigated as potential therapeutic avenues in preclinical and clinical studies. However, how pre-mRNA splicing and metabolism are orchestrated in response to oncogenic stress and therapies is poorly understood. Here, we demonstrate that jumonji domain containing 6, arginine demethylase, and lysine hydroxylase, JMJD6, acts as a hub connecting splicing and metabolism in MYC-driven human neuroblastoma. JMJD6 cooperates with MYC in cellular transformation of murine neural crest cells by physically interacting with RNA binding proteins involved in pre-mRNA splicing and protein homeostasis. Notably, JMJD6 controls the alternative splicing of two isoforms of glutaminase (GLS), namely kidney-type glutaminase (KGA) and glutaminase C (GAC), which are rate-limiting enzymes of glutaminolysis in the central carbon metabolism in neuroblastoma. Further, we show that JMJD6 is correlated with the anti-cancer activity of indisulam, a 'molecular glue' that degrades splicing factor RBM39, which complexes with JMJD6. The indisulam-mediated cancer cell killing is at least partly dependent on the glutamine-related metabolic pathway mediated by JMJD6. Our findings reveal a cancer-promoting metabolic program is associated with alternative pre-mRNA splicing through JMJD6, providing a rationale to target JMJD6 as a therapeutic avenue for treating MYC-driven cancers.

Histone lysine demethylase 4 family proteins maintain the transcriptional program and adrenergic cellular state of MYCN-amplified neuroblastoma.

Neuroblastoma with MYCN amplification (MNA) is a high-risk disease that has a poor survival rate. Neuroblastoma displays cellular heterogeneity, including more differentiated (adrenergic) and more primitive (mesenchymal) cellular states. Here, we demonstrate that MYCN oncoprotein promotes a cellular state switch in mesenchymal cells to an adrenergic state, accompanied by induction of histone lysine demethylase 4 family members (KDM4A-C) that act in concert to control the expression of MYCN and adrenergic core regulatory circulatory (CRC) transcription factors. Pharmacologic inhibition of KDM4 blocks expression of MYCN and the adrenergic CRC transcriptome with genome-wide induction of transcriptionally repressive H3K9me3, resulting in potent anticancer activity against neuroblastomas with MNA by inducing neuroblastic differentiation and apoptosis. Furthermore, a short-term KDM4 inhibition in combination with conventional, cytotoxic chemotherapy results in complete tumor responses of xenografts with MNA. Thus, KDM4 blockade may serve as a transformative strategy to target the adrenergic CRC dependencies in MNA neuroblastomas.

A spatial cell atlas of neuroblastoma reveals developmental, epigenetic and spatial axis of tumor heterogeneity.

Neuroblastoma is a pediatric cancer arising from the developing sympathoadrenal lineage with complex inter- and intra-tumoral heterogeneity. To chart this complexity, we generated a comprehensive cell atlas of 55 neuroblastoma patient tumors, collected from two pediatric cancer institutions, spanning a range of clinical, genetic, and histologic features. Our atlas combines single-cell/nucleus RNA-seq (sc/scRNA-seq), bulk RNA-seq, whole exome sequencing, DNA methylation profiling, spatial transcriptomics, and two spatial proteomic methods. Sc/snRNA-seq revealed three malignant cell states with features of sympathoadrenal lineage development. All of the neuroblastomas had malignant cells that resembled sympathoblasts and the more differentiated adrenergic cells. A subset of tumors had malignant cells in a mesenchymal cell state with molecular features of Schwann cell precursors. DNA methylation profiles defined four groupings of patients, which differ in the degree of malignant cell heterogeneity and clinical outcomes. Using spatial proteomics, we found that neuroblastomas are spatially compartmentalized, with malignant tumor cells sequestered away from immune cells. Finally, we identify spatially restricted signaling patterns in immune cells from spatial transcriptomics. To facilitate the visualization and analysis of our atlas as a resource for further research in neuroblastoma, single cell, and spatial-omics, all data are shared through the Human Tumor Atlas Network Data Commons at www.humantumoratlas.org.

PAX3-FOXO1 dictates myogenic reprogramming and rhabdomyosarcoma identity in endothelial progenitors.

Fusion-positive rhabdomyosarcoma (FP-RMS) driven by the expression of the PAX3-FOXO1 (P3F) fusion oncoprotein is an aggressive subtype of pediatric rhabdomyosarcoma. FP-RMS histologically resembles developing muscle yet occurs throughout the body in areas devoid of skeletal muscle highlighting that FP-RMS is not derived from an exclusively myogenic cell of origin. Here we demonstrate that P3F reprograms mouse and human endothelial progenitors to FP-RMS. We show that P3F expression in aP2-Cre expressing cells reprograms endothelial progenitors to functional myogenic stem cells capable of regenerating injured muscle fibers. Further, we describe a FP-RMS mouse model driven by P3F expression and Cdkn2a loss in endothelial cells. Additionally, we show that P3F expression in TP53-null human iPSCs blocks endothelial-directed differentiation and guides cells to become myogenic cells that form FP-RMS tumors in immunocompromised mice. Together these findings demonstrate that FP-RMS can originate from aberrant development of non-myogenic cells driven by P3F.

A Review of Single-Cell RNA-Seq Annotation, Integration, and Cell-Cell Communication.

Single-cell RNA sequencing (scRNA-seq) has emerged as a powerful tool for investigating cellular biology at an unprecedented resolution, enabling the characterization of cellular heterogeneity, identification of rare but significant cell types, and exploration of cell-cell communications and interactions. Its broad applications span both basic and clinical research domains. In this comprehensive review, we survey the current landscape of scRNA-seq analysis methods and tools, focusing on count modeling, cell-type annotation, data integration, including spatial transcriptomics, and the inference of cell-cell communication. We review the challenges encountered in scRNA-seq analysis, including issues of sparsity or low expression, reliability of cell annotation, and assumptions in data integration, and discuss the potential impact of suboptimal clustering and differential expression analysis tools on downstream analyses, particularly in identifying cell subpopulations. Finally, we discuss recent advancements and future directions for enhancing scRNA-seq analysis. Specifically, we highlight the development of novel tools for annotating single-cell data, integrating and interpreting multimodal datasets covering transcriptomics, epigenomics, and proteomics, and inferring cellular communication networks. By elucidating the latest progress and innovation, we provide a comprehensive overview of the rapidly advancing field of scRNA-seq analysis.

Genome-wide mapping of cancer dependency genes and genetic modifiers of chemotherapy in high-risk hepatoblastoma.

A lack of relevant genetic models and cell lines hampers our understanding of hepatoblastoma pathogenesis and the development of new therapies for this neoplasm. Here, we report an improved MYC-driven hepatoblastoma-like murine model that recapitulates the pathological features of embryonal type of hepatoblastoma, with transcriptomics resembling the high-risk gene signatures of the human disease. Single-cell RNA-sequencing and spatial transcriptomics identify distinct subpopulations of hepatoblastoma cells. After deriving cell lines from the mouse model, we map cancer dependency genes using CRISPR-Cas9 screening and identify druggable targets shared with human hepatoblastoma (e.g., CDK7, CDK9, PRMT1, PRMT5). Our screen also reveals oncogenes and tumor suppressor genes in hepatoblastoma that engage multiple, druggable cancer signaling pathways. Chemotherapy is critical for human hepatoblastoma treatment. A genetic mapping of doxorubicin response by CRISPR-Cas9 screening identifies modifiers whose loss-of-function synergizes with (e.g., PRKDC) or antagonizes (e.g., apoptosis genes) the effect of chemotherapy. The combination of PRKDC inhibition and doxorubicin-based chemotherapy greatly enhances therapeutic efficacy. These studies provide a set of resources including disease models suitable for identifying and validating potential therapeutic targets in human high-risk hepatoblastoma.

Metabolic reprogramming of cancer cells by JMJD6-mediated pre-mRNA splicing is associated with therapeutic response to splicing inhibitor.

Dysregulated pre-mRNA splicing and metabolism are two hallmarks of MYC-driven cancers. Pharmacological inhibition of both processes has been extensively investigated as potential therapeutic avenues in preclinical and clinical studies. However, how pre-mRNA splicing and metabolism are orchestrated in response to oncogenic stress and therapies is poorly understood. Here, we demonstrate that Jumonji Domain Containing 6, Arginine Demethylase and Lysine Hydroxylase, JMJD6, acts as a hub connecting splicing and metabolism in MYC-driven neuroblastoma. JMJD6 cooperates with MYC in cellular transformation by physically interacting with RNA binding proteins involved in pre-mRNA splicing and protein homeostasis. Notably, JMJD6 controls the alternative splicing of two isoforms of glutaminase (GLS), namely kidney-type glutaminase (KGA) and glutaminase C (GAC), which are rate-limiting enzymes of glutaminolysis in the central carbon metabolism in neuroblastoma. Further, we show that JMJD6 is correlated with the anti-cancer activity of indisulam, a "molecular glue" that degrades splicing factor RBM39, which complexes with JMJD6. The indisulam-mediated cancer cell killing is at least partly dependent on the glutamine-related metabolic pathway mediated by JMJD6. Our findings reveal a cancer-promoting metabolic program is associated with alternative pre-mRNA splicing through JMJD6, providing a rationale to target JMJD6 as a therapeutic avenue for treating MYC-driven cancers.

A kinase-independent function of cyclin-dependent kinase 6 promotes outer radial glia expansion and neocortical folding.

The neocortex, the center for higher brain function, first emerged in mammals and has become massively expanded and folded in humans, constituting almost half the volume of the human brain. Primary microcephaly, a developmental disorder in which the brain is smaller than normal at birth, results mainly from there being fewer neurons in the neocortex because of defects in neural progenitor cells (NPCs). Outer radial glia (oRGs), NPCs that are abundant in gyrencephalic species but rare in lissencephalic species, are thought to play key roles in the expansion and folding of the neocortex. However, how oRGs expand, whether they are necessary for neocortical folding, and whether defects in oRGs cause microcephaly remain important questions in the study of brain development, evolution, and disease. Here, we show that oRG expansion in mice, ferrets, and human cerebral organoids requires cyclin-dependent kinase 6 (CDK6), the mutation of which causes primary microcephaly via an unknown mechanism. In a mouse model in which increased Hedgehog signaling expands oRGs and intermediate progenitor cells and induces neocortical folding, CDK6 loss selectively decreased oRGs and abolished neocortical folding. Remarkably, this function of CDK6 in oRG expansion did not require its kinase activity, was not shared by the highly similar CDK4 and CDK2, and was disrupted by the mutation causing microcephaly. Therefore, our results indicate that CDK6 is conserved to promote oRG expansion, that oRGs are necessary for neocortical folding, and that defects in oRG expansion may cause primary microcephaly.

Targeting KDM4 for treating PAX3-FOXO1-driven alveolar rhabdomyosarcoma.

Chimeric transcription factors drive lineage-specific oncogenesis but are notoriously difficult to target. Alveolar rhabdomyosarcoma (RMS) is an aggressive childhood soft tissue sarcoma transformed by the pathognomonic Paired Box 3-Forkhead Box O1 (PAX3-FOXO1) fusion protein, which governs a core regulatory circuitry transcription factor network. Here, we show that the histone lysine demethylase 4B (KDM4B) is a therapeutic vulnerability for PAX3-FOXO1+ RMS. Genetic and pharmacologic inhibition of KDM4B substantially delayed tumor growth. Suppression of KDM4 proteins inhibited the expression of core oncogenic transcription factors and caused epigenetic alterations of PAX3-FOXO1-governed superenhancers. Combining KDM4 inhibition with cytotoxic chemotherapy led to tumor regression in preclinical PAX3-FOXO1+ RMS subcutaneous xenograft models. In summary, we identified a targetable mechanism required for maintenance of the PAX3-FOXO1-related transcription factor network, which may translate to a therapeutic approach for fusion-positive RMS.

The molecular characteristics of low-grade and high-grade areas in desmoplastic infantile astrocytoma/ganglioglioma.

AIMS: Desmoplastic infantile astrocytomas and gangliogliomas (DIA/DIGs) are rare brain tumours of infancy. A distinctive feature of their histopathology is a combination of low-grade and high-grade features. Most DIA/DIGs can be surgically resected and have a good prognosis. However, high-grade features often dominate recurrent tumours, some of which have a poor outcome. In this study, we test the hypothesis that low-grade and high-grade areas in DIA/DIGs have distinct molecular characteristics. METHODS: Tissue samples from microdissected low-grade and high-grade areas in 12 DIA/DIGs were analysed by DNA methylation profiling, whole exome sequencing, RNA sequencing and immunohistochemistry to search for potential differences at multiple molecular levels. RESULTS: Copy number variants among tumours and between the two morphologically distinct areas were infrequent. No recurrent genetic alterations were identified across the tumour series, and high-grade areas did not have additional genetic alterations to explain their distinct morphology or biological behaviour. However, high-grade areas showed relative hypomethylation in genes downstream of the transcription factors SOX9 and LEF1 and evidence of a core SOX9 transcription network alongside activation of the BMP, WNT and MAPK signalling pathways. CONCLUSIONS: This study contributes to our knowledge of molecular genetic alterations in DIA/DIGs, uncovers molecular differences between the two distinct cell populations in these tumours and suggests potential therapeutic targets among the more proliferative cell population in DIA/DIGs.

Chromatin architecture at susceptible gene loci in cerebellar Purkinje cells characterizes DNA damage-induced neurodegeneration.

The pathogenesis of inherited genome instability neurodegenerative syndromes remains largely unknown. Here, we report new disease-relevant murine models of genome instability­driven neurodegeneration involving disabled ATM and APTX that develop debilitating ataxia. We show that neurodegeneration and ataxia result from transcriptional interference in the cerebellum via aberrant messenger RNA splicing. Unexpectedly, these splicing defects were restricted to only Purkinje cells, disrupting the expression of critical homeostatic regulators including ITPR1, GRID2, and CA8. Abundant genotoxic R loops were also found at these Purkinje cell gene loci, further exacerbating DNA damage and transcriptional disruption. Using ATAC-seq to profile global chromatin accessibility in the cerebellum, we found a notably unique chromatin conformation specifically in Purkinje chromatin at the affected gene loci, thereby promoting susceptibility to DNA damage. These data reveal the pathogenic basis of DNA damage in the nervous system and suggest chromatin conformation as a feature in directing genome instability­associated neuropathology.

KDM6B promotes activation of the oncogenic CDK4/6-pRB-E2F pathway by maintaining enhancer activity in MYCN-amplified neuroblastoma.

The H3K27me2/me3 histone demethylase KDM6B is essential to neuroblastoma cell survival. However, the mechanism of KDM6B action remains poorly defined. We demonstrate that inhibition of KDM6B activity 1) reduces the chromatin accessibility of E2F target genes and MYCN, 2) selectively leads to an increase of H3K27me3 but a decrease of the enhancer mark H3K4me1 at the CTCF and BORIS binding sites, which may, consequently, disrupt the long-range chromatin interaction of MYCN and E2F target genes, and 3) phenocopies the transcriptome induced by the specific CDK4/6 inhibitor palbociclib. Overexpression of CDK4/6 or Rb1 knockout confers neuroblastoma cell resistance to both palbociclib and the KDM6 inhibitor GSK-J4. These data indicate that KDM6B promotes an oncogenic CDK4/6-pRB-E2F pathway in neuroblastoma cells via H3K27me3-dependent enhancer-promoter interactions, providing a rationale to target KDM6B for high-risk neuroblastoma.

Targeting the spliceosome through RBM39 degradation results in exceptional responses in high-risk neuroblastoma models.

Aberrant alternative pre-mRNA splicing plays a critical role in MYC-driven cancers and therefore may represent a therapeutic vulnerability. Here, we show that neuroblastoma, a MYC-driven cancer characterized by splicing dysregulation and spliceosomal dependency, requires the splicing factor RBM39 for survival. Indisulam, a "molecular glue" that selectively recruits RBM39 to the CRL4-DCAF15 E3 ubiquitin ligase for proteasomal degradation, is highly efficacious against neuroblastoma, leading to significant responses in multiple high-risk disease models, without overt toxicity. Genetic depletion or indisulam-mediated degradation of RBM39 induces significant genome-wide splicing anomalies and cell death. Mechanistically, the dependency on RBM39 and high-level expression of DCAF15 determine the exquisite sensitivity of neuroblastoma to indisulam. Our data indicate that targeting the dysregulated spliceosome by precisely inhibiting RBM39, a vulnerability in neuroblastoma, is a valid therapeutic strategy.

A genetic mouse model with postnatal Nf1 and p53 loss recapitulates the histology and transcriptome of human malignant peripheral nerve sheath tumor.

BACKGROUND: Malignant peripheral nerve sheath tumors (MPNST) are aggressive sarcomas. Somatic inactivation of NF1 and cooperating tumor suppressors, including CDKN2A/B, PRC2, and p53, is found in most MPNST. Inactivation of LATS1/2 of the Hippo pathway was recently shown to cause tumors resembling MPNST histologically, although Hippo pathway mutations are rarely found in MPNST. Because existing genetically engineered mouse (GEM) models of MPNST do not recapitulate some of the key genetic features of human MPNST, we aimed to establish a GEM-MPNST model that recapitulated the human disease genetically, histologically, and molecularly. METHODS: We combined 2 genetically modified alleles, an Nf1;Trp53 cis-conditional allele and an inducible Plp-CreER allele (NP-Plp), to model the somatic, possibly postnatal, mutational events in human MPNST. We also generated conditional Lats1;Lats2 knockout mice. We performed histopathologic analyses of mouse MPNST models and transcriptomic comparison of mouse models and human nerve sheath tumors. RESULTS: Postnatal Nf1;Trp53 cis-deletion resulted in GEM-MPNST that were histologically more similar to human MPNST than the widely used germline Nf1;Trp53 cis-heterozygous (NPcis) model and showed partial loss of H3K27me3. At the transcriptome level, Nf1;p53-driven GEM-MPNST were distinct from Lats-driven GEM-MPNST and resembled human MPNST more closely than do Lats-driven tumors. CONCLUSIONS: The NP-Plp model recapitulates human MPNST genetically, histologically, and molecularly.

Synthetic essentiality between PTEN and core dependency factor PAX7 dictates rhabdomyosarcoma identity.

PTEN promoter hypermethylation is nearly universal and PTEN copy number loss occurs in ~25% of fusion-negative rhabdomyosarcoma (FN-RMS). Here we show Pten deletion in a mouse model of FN-RMS results in less differentiated tumors more closely resembling human embryonal RMS. PTEN loss activated the PI3K pathway but did not increase mTOR activity. In wild-type tumors, PTEN was expressed in the nucleus suggesting loss of nuclear PTEN functions could account for these phenotypes. Pten deleted tumors had increased expression of transcription factors important in neural and skeletal muscle development including Dbx1 and Pax7. Pax7 deletion completely rescued the effects of Pten loss. Strikingly, these Pten;Pax7 deleted tumors were no longer FN-RMS but displayed smooth muscle differentiation similar to leiomyosarcoma. These data highlight how Pten loss in FN-RMS is connected to a PAX7 lineage-specific transcriptional output that creates a dependency or synthetic essentiality on the transcription factor PAX7 to maintain tumor identity.

Retinoblastoma from human stem cell-derived retinal organoids.

Retinoblastoma is a childhood cancer of the developing retina that initiates with biallelic inactivation of the RB1 gene. Children with germline mutations in RB1 have a high likelihood of developing retinoblastoma and other malignancies later in life. Genetically engineered mouse models of retinoblastoma share some similarities with human retinoblastoma but there are differences in their cellular differentiation. To develop a laboratory model of human retinoblastoma formation, we make induced pluripotent stem cells (iPSCs) from 15 participants with germline RB1 mutations. Each of the stem cell lines is validated, characterized and then differentiated into retina using a 3-dimensional organoid culture system. After 45 days in culture, the retinal organoids are dissociated and injected into the vitreous of eyes of immunocompromised mice to support retinoblastoma tumor growth. Retinoblastomas formed from retinal organoids made from patient-derived iPSCs have molecular, cellular and genomic features indistinguishable from human retinoblastomas. This model of human cancer based on patient-derived iPSCs with germline cancer predisposing mutations provides valuable insights into the cellular origins of this debilitating childhood disease as well as the mechanism of tumorigenesis following RB1 gene inactivation.

CovidExpress: an interactive portal for intuitive investigation on SARS-CoV-2 related transcriptomes.

Infection with severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) in humans could cause coronavirus disease 2019 (COVID-19). Since its first discovery in Dec 2019, SARS-CoV-2 has become a global pandemic and caused 3.3 million direct/indirect deaths (2021 May). Amongst the scientific community's response to COVID-19, data sharing has emerged as an essential aspect of the combat against SARS-CoV-2. Despite the ever-growing studies about SARS-CoV-2 and COVID-19, to date, only a few databases were curated to enable access to gene expression data. Furthermore, these databases curated only a small set of data and do not provide easy access for investigators without computational skills to perform analyses. To fill this gap and advance open-access to the growing gene expression data on this deadly virus, we collected about 1,500 human bulk RNA-seq datasets from publicly available resources, developed a database and visualization tool, named CovidExpress (https://stjudecab.github.io/covidexpress). This open access database will allow research investigators to examine the gene expression in various tissues, cell lines, and their response to SARS-CoV-2 under different experimental conditions, accelerating the understanding of the etiology of this disease to inform the drug and vaccine development. Our integrative analysis of this big dataset highlights a set of commonly regulated genes in SARS-CoV-2 infected lung and Rhinovirus infected nasal tissues, including OASL that were under-studied in COVID-19 related reports. Our results also suggested a potential FURIN positive feedback loop that might explain the evolutional advantage of SARS-CoV-2.