Integrated multi-omics analyses identify anti-viral host factors and pathways controlling SARS-CoV-2 infection.

Host anti-viral factors are essential for controlling SARS-CoV-2 infection but remain largely unknown due to the biases of previous large-scale studies toward pro-viral host factors. To fill in this knowledge gap, we perform a genome-wide CRISPR dropout screen and integrate analyses of the multi-omics data of the CRISPR screen, genome-wide association studies, single-cell RNA-Seq, and host-virus proteins or protein/RNA interactome. This study uncovers many host factors that are currently underappreciated, including the components of V-ATPases, ESCRT, and N-glycosylation pathways that modulate viral entry and/or replication. The cohesin complex is also identified as an anti-viral pathway, suggesting an important role of three-dimensional chromatin organization in mediating host-viral interaction. Furthermore, we discover another anti-viral regulator KLF5, a transcriptional factor involved in sphingolipid metabolism, which is up-regulated, and harbors genetic variations linked to COVID-19 patients with severe symptoms. Anti-viral effects of three identified candidates (DAZAP2/VTA1/KLF5) are confirmed individually. Molecular characterization of DAZAP2/VTA1/KLF5-knockout cells highlights the involvement of genes related to the coagulation system in determining the severity of COVID-19. Together, our results provide further resources for understanding the host anti-viral network during SARS-CoV-2 infection and may help develop new countermeasure strategies.

Faslpr gene dosage tunes the extent of lymphoproliferation and T cell differentiation in lupus.

Sle1 and Faslpr are two lupus susceptibility loci that lead to manifestations of systemic lupus erythematosus. To evaluate the dosage effects of Faslpr in determining cellular and serological phenotypes associated with lupus, we developed a new C57BL/6 (B6) congenic lupus strain, B6.Sle1/Sle1.Faslpr/+ (Sle1homo.lprhet) and compared it with B6.Faslpr/lpr (lprhomo), B6.Sle1/Sle1 (Sle1homo), and B6.Sle1/Sle1.Faslpr/lpr (Sle1homo.lprhomo) strains. Whereas Sle1homo.lprhomo mice exhibited profound lymphoproliferation and early mortality, Sle1homo.lprhet mice had a lifespan comparable to B6 mice, with no evidence of splenomegaly or lymphadenopathy. Compared to B6 monogenic lupus strains, Sle1homo.lprhet mice exhibited significantly elevated serum ANA antibodies and increased proteinuria. Additionally, Sle1homo.lprhet T cells had an increased propensity to differentiate into Th1 cells. Gene dose effects of Faslpr were noted in upregulating serum IL-1⍺, IL-2, and IL-27. Taken together, Sle1homo.lprhet strain is a new C57BL/6-based model of lupus, ideal for genetic studies, autoantibody repertoire investigation, and for exploring Th1 effector cell skewing without early-age lymphoproliferative autoimmunity.

Multiomics analyses reveal DARS1-AS1/YBX1-controlled posttranscriptional circuits promoting glioblastoma tumorigenesis/radioresistance.

The glioblastoma (GBM) stem cell-like cells (GSCs) are critical for tumorigenesis/therapeutic resistance of GBM. Mounting evidence supports tumor-promoting function of long noncoding RNAs (lncRNAs), but their role in GSCs remains poorly understood. By combining CRISPRi screen with orthogonal multiomics approaches, we identified a lncRNA DARS1-AS1-controlled posttranscriptional circuitry that promoted the malignant properties of GBM cells/GSCs. Depleting DARS1-AS1 inhibited the proliferation of GBM cells/GSCs and self-renewal of GSCs, prolonging survival in orthotopic GBM models. DARS1-AS1 depletion also impaired the homologous recombination (HR)-mediated double-strand break (DSB) repair and enhanced the radiosensitivity of GBM cells/GSCs. Mechanistically, DARS1-AS1 interacted with YBX1 to promote target mRNA binding and stabilization, forming a mixed transcriptional/posttranscriptional feed-forward loop to up-regulate expression of the key regulators of G1-S transition, including E2F1 and CCND1. DARS1-AS1/YBX1 also stabilized the mRNA of FOXM1, a master transcription factor regulating GSC self-renewal and DSB repair. Our findings suggest DARS1-AS1/YBX1 axis as a potential therapeutic target for sensitizing GBM to radiation/HR deficiency-targeted therapy.

CRISPR/Cas9 screen uncovers functional translation of cryptic lncRNA-encoded open reading frames in human cancer.

Emerging evidence suggests that cryptic translation within long noncoding RNAs (lncRNAs) may produce novel proteins with important developmental/physiological functions. However, the role of this cryptic translation in complex diseases (e.g., cancer) remains elusive. Here, we applied an integrative strategy combining ribosome profiling and CRISPR/Cas9 screening with large-scale analysis of molecular/clinical data for breast cancer (BC) and identified estrogen receptor α-positive (ER+) BC dependency on the cryptic ORFs encoded by lncRNA genes that were upregulated in luminal tumors. We confirmed the in vivo tumor-promoting function of an unannotated protein, GATA3-interacting cryptic protein (GT3-INCP) encoded by LINC00992, the expression of which was associated with poor prognosis in luminal tumors. GTE-INCP was upregulated by estrogen/ER and regulated estrogen-dependent cell growth. Mechanistically, GT3-INCP interacted with GATA3, a master transcription factor key to mammary gland development/BC cell proliferation, and coregulated a gene expression program that involved many BC susceptibility/risk genes and impacted estrogen response/cell proliferation. GT3-INCP/GATA3 bound to common cis regulatory elements and upregulated the expression of the tumor-promoting and estrogen-regulated BC susceptibility/risk genes MYB and PDZK1. Our study indicates that cryptic lncRNA-encoded proteins can be an important integrated component of the master transcriptional regulatory network driving aberrant transcription in cancer, and suggests that the "hidden" lncRNA-encoded proteome might be a new space for therapeutic target discovery.

Single-cell CRISPR immune screens reveal immunological roles of tumor intrinsic factors.

Genetic screens are widely exploited to develop novel therapeutic approaches for cancer treatment. With recent advances in single-cell technology, single-cell CRISPR screen (scCRISPR) platforms provide opportunities for target validation and mechanistic studies in a high-throughput manner. Here, we aim to establish scCRISPR platforms which are suitable for immune-related screens involving multiple cell types. We integrated two scCRISPR platforms, namely Perturb-seq and CROP-seq, with both in vitro and in vivo immune screens. By leveraging previously generated resources, we optimized experimental conditions and data analysis pipelines to achieve better consistency between results from high-throughput and individual validations. Furthermore, we evaluated the performance of scCRISPR immune screens in determining underlying mechanisms of tumor intrinsic immune regulation. Our results showed that scCRISPR platforms can simultaneously characterize gene expression profiles and perturbation effects present in individual cells in different immune screen conditions. Results from scCRISPR immune screens also predict transcriptional phenotype associated with clinical responses to cancer immunotherapy. More importantly, scCRISPR screen platforms reveal the interactive relationship between targeting tumor intrinsic factors and T cell-mediated antitumor immune response which cannot be easily assessed by bulk RNA-seq. Collectively, scCRISPR immune screens provide scalable and reliable platforms to elucidate molecular determinants of tumor immune resistance.

Integrated multi-omics analyses identify key anti-viral host factors and pathways controlling SARS-CoV-2 infection.

Host anti-viral factors are essential for controlling SARS-CoV-2 infection but remain largely unknown due to the biases of previous large-scale studies toward pro-viral host factors. To fill in this knowledge gap, we performed a genome-wide CRISPR dropout screen and integrated analyses of the multi-omics data of the CRISPR screen, genome-wide association studies, single-cell RNA-seq, and host-virus proteins or protein/RNA interactome. This study has uncovered many host factors that were missed by previous studies, including the components of V-ATPases, ESCRT, and N-glycosylation pathways that modulated viral entry and/or replication. The cohesin complex was also identified as a novel anti-viral pathway, suggesting an important role of three-dimensional chromatin organization in mediating host-viral interaction. Furthermore, we discovered an anti-viral regulator KLF5, a transcriptional factor involved in sphingolipid metabolism, which was up-regulated and harbored genetic variations linked to the COVID-19 patients with severe symptoms. Our results provide a resource for understanding the host anti-viral network during SARS-CoV-2 infection and may help develop new countermeasure strategies.

Inhibiting Type I Arginine Methyltransferase Activity Promotes T Cell-Mediated Antitumor Immune Responses.

Protein arginine methyltransferases (PRMT) are a widely expressed class of enzymes responsible for catalyzing arginine methylation on numerous protein substrates. Among them, type I PRMTs are responsible for generating asymmetric dimethylarginine. By controlling multiple basic cellular processes, such as DNA damage responses, transcriptional regulation, and mRNA splicing, type I PRMTs contribute to cancer initiation and progression. A type I PRMT inhibitor, GSK3368715, has been developed and has entered clinical trials for solid and hematologic malignancies. Although type I PRMTs have been reported to play roles in modulating immune cell function, the immunologic role of tumor-intrinsic pathways controlled by type I PRMTs remains uncharacterized. Here, our The Cancer Genome Atlas dataset analysis revealed that expression of type I PRMTs associated with poor clinical response and decreased immune infiltration in patients with melanoma. In cancer cell lines, inhibition of type I PRMTs induced an IFN gene signature, amplified responses to IFN and innate immune signaling, and decreased expression of the immunosuppressive cytokine VEGF. In immunocompetent mouse tumor models, including a model of T-cell exclusion that represents a common mechanism of anti-programmed cell death protein 1 (PD-1) resistance in humans, type I PRMT inhibition increased T-cell infiltration, produced durable responses dependent on CD8+ T cells, and enhanced efficacy of anti-PD-1 therapy. These data indicate that type I PRMT inhibition exhibits immunomodulatory properties and synergizes with immune checkpoint blockade (ICB) to induce durable antitumor responses in a T cell-dependent manner, suggesting that type I PRMT inhibition can potentiate an antitumor immunity in refractory settings.

Systematic functional interrogation of human pseudogenes using CRISPRi.

BACKGROUND: The human genome encodes over 14,000 pseudogenes that are evolutionary relics of protein-coding genes and commonly considered as nonfunctional. Emerging evidence suggests that some pseudogenes may exert important functions. However, to what extent human pseudogenes are functionally relevant remains unclear. There has been no large-scale characterization of pseudogene function because of technical challenges, including high sequence similarity between pseudogene and parent genes, and poor annotation of transcription start sites. RESULTS: To overcome these technical obstacles, we develop an integrated computational pipeline to design the first genome-wide library of CRISPR interference (CRISPRi) single-guide RNAs (sgRNAs) that target human pseudogene promoter-proximal regions. We perform the first pseudogene-focused CRISPRi screen in luminal A breast cancer cells and reveal approximately 70 pseudogenes that affect breast cancer cell fitness. Among the top hits, we identify a cancer-testis unitary pseudogene, MGAT4EP, that is predominantly localized in the nucleus and interacts with FOXA1, a key regulator in luminal A breast cancer. By enhancing the promoter binding of FOXA1, MGAT4EP upregulates the expression of oncogenic transcription factor FOXM1. Integrative analyses of multi-omic data from the Cancer Genome Atlas (TCGA) reveal many unitary pseudogenes whose expressions are significantly dysregulated and/or associated with overall/relapse-free survival of patients in diverse cancer types. CONCLUSIONS: Our study represents the first large-scale study characterizing pseudogene function. Our findings suggest the importance of nuclear function of unitary pseudogenes and underscore their underappreciated roles in human diseases. The functional genomic resources developed here will greatly facilitate the study of human pseudogene function.

Integrating genome-wide CRISPR immune screen with multi-omic clinical data reveals distinct classes of tumor intrinsic immune regulators.

BACKGROUND: Despite approval of immunotherapy for a wide range of cancers, the majority of patients fail to respond to immunotherapy or relapse following initial response. These failures may be attributed to immunosuppressive mechanisms co-opted by tumor cells. However, it is challenging to use conventional methods to systematically evaluate the potential of tumor intrinsic factors to act as immune regulators in patients with cancer. METHODS: To identify immunosuppressive mechanisms in non-responders to cancer immunotherapy in an unbiased manner, we performed genome-wide CRISPR immune screens and integrated our results with multi-omics clinical data to evaluate the role of tumor intrinsic factors in regulating two rate-limiting steps of cancer immunotherapy, namely, T cell tumor infiltration and T cell-mediated tumor killing. RESULTS: Our studies revealed two distinct types of immune resistance regulators and demonstrated their potential as therapeutic targets to improve the efficacy of immunotherapy. Among them, PRMT1 and RIPK1 were identified as a dual immune resistance regulator and a cytotoxicity resistance regulator, respectively. Although the magnitude varied between different types of immunotherapy, genetically targeting PRMT1 and RIPK1 sensitized tumors to T-cell killing and anti-PD-1/OX40 treatment. Interestingly, a RIPK1-specific inhibitor enhanced the antitumor activity of T cell-based and anti-OX40 therapy, despite limited impact on T cell tumor infiltration. CONCLUSIONS: Collectively, the data provide a rich resource of novel targets for rational immuno-oncology combinations.

Genome-wide identification and differential analysis of translational initiation.

Translation is principally regulated at the initiation stage. The development of the translation initiation (TI) sequencing (TI-seq) technique has enabled the global mapping of TIs and revealed unanticipated complex translational landscapes in metazoans. Despite the wide adoption of TI-seq, there is no computational tool currently available for analyzing TI-seq data. To fill this gap, we develop a comprehensive toolkit named Ribo-TISH, which allows for detecting and quantitatively comparing TIs across conditions from TI-seq data. Ribo-TISH can also predict novel open reading frames (ORFs) from regular ribosome profiling (rRibo-seq) data and outperform several established methods in both computational efficiency and prediction accuracy. Applied to published TI-seq/rRibo-seq data sets, Ribo-TISH uncovers a novel signature of elevated mitochondrial translation during amino-acid deprivation and predicts novel ORFs in 5'UTRs, long noncoding RNAs, and introns. These successful applications demonstrate the power of Ribo-TISH in extracting biological insights from TI-seq/rRibo-seq data.

Long Noncoding RNA LINC00673 Is Activated by SP1 and Exerts Oncogenic Properties by Interacting with LSD1 and EZH2 in Gastric Cancer.

Long noncoding RNAs (lncRNAs) have emerged as important regulators in a variety of human diseases, including cancers. However, the biological function of these molecules and the mechanisms responsible for their alteration in gastric cancer (GC) are not fully understood. In this study, we found that lncRNA LINC00673 is significantly upregulated in gastric cancer. Knockdown of LINC00673 inhibited cell proliferation and invasion and induced cell apoptosis, whereas LINC00673 overexpression had the opposite effect. Online transcription factor binding site prediction analysis showed that there are SP1 binding sites in the LINC00673 promoter region. Next, luciferase reporter and chromatin immunoprecipitation (ChIP) assays provided evidence that SP1 could bind directly to the LINC00673 promoter region and activate its transcription. Moreover, mechanistic investigation showed that CADM4, KLF2, and LATS2 might be the underlying targets of LINC00673 in GC cells, and RNA immunoprecipitation, RNA pull-down, and ChIP assays showed that LINC00673 can interact with EZH2 and LSD1, thereby repressing KLF2 and LATS2 expression. Taken together, these findings show that SP1-activated LINC00673 exerts an oncogenic function that promotes GC development and progression, at least in part, by functioning as a scaffold for LSD1 and EZH2 and repressing KLF2 and LATS2 expression.

The Pseudogene DUXAP8 Promotes Non-small-cell Lung Cancer Cell Proliferation and Invasion by Epigenetically Silencing EGR1 and RHOB.

Recently, the non-protein-coding functional elements in the human genome have been identified as key regulators in postgenomic biology, and a large number of pseudogenes as well as long non-coding RNAs (lncRNAs) have been found to be transcribed in multiple human cancers. However, only a small proportion of these pseudogenes has been functionally characterized. In this study, we screened for pseudogenes associated with human non-small-cell lung cancer (NSCLC) by comparative analysis of several independent datasets from the GEO. We identified a transcribed pseudogene named DUXAP8 that is upregulated in tumor tissues. Patients with higher DUXAP8 expression exhibited shorter survival, suggesting DUXAP8 as a new candidate prognostic marker for NSCLC patients. Knockdown of DUXAP8 impairs cell growth, migration, and invasion, and induces apoptosis both in vitro and in vivo. Mechanistically, DUXAP8 represses the tumor suppressors EGR1 and RHOB by recruiting histone demethylase LSD1 and histone methyltransferase EZH2, thereby promoting cell proliferation, migration, and invasion. These findings indicate that the pseudogene DUXAP8 may act as an oncogene in NSCLC by silencing EGR1 and RHOB transcription by binding with EZH2 and LSD1, which may offer a novel therapeutic target for this disease.

LncRNA HOXA11-AS Promotes Proliferation and Invasion of Gastric Cancer by Scaffolding the Chromatin Modification Factors PRC2, LSD1, and DNMT1.

Long noncoding RNAs (lncRNA) have been implicated in human cancer but their mechanisms of action are mainly undocumented. In this study, we investigated lncRNA alterations that contribute to gastric cancer through an analysis of The Cancer Genome Atlas RNA sequencing data and other publicly available microarray data. Here we report the gastric cancer-associated lncRNA HOXA11-AS as a key regulator of gastric cancer development and progression. Patients with high HOXA11-AS expression had a shorter survival and poorer prognosis. In vitro and in vivo assays of HOXA11-AS alterations revealed a complex integrated phenotype affecting cell growth, migration, invasion, and apoptosis. Strikingly, high-throughput sequencing analysis after HOXA11-AS silencing highlighted alterations in cell proliferation and cell-cell adhesion pathways. Mechanistically, EZH2 along with the histone demethylase LSD1 or DNMT1 were recruited by HOXA11-AS, which functioned as a scaffold. HOXA11-AS also functioned as a molecular sponge for miR-1297, antagonizing its ability to repress EZH2 protein translation. In addition, we found that E2F1 was involved in HOXA11-AS activation in gastric cancer cells. Taken together, our findings support a model in which the EZH2/HOXA11-AS/LSD1 complex and HOXA11-AS/miR-1297/EZH2 cross-talk serve as critical effectors in gastric cancer tumorigenesis and progression, suggesting new therapeutic directions in gastric cancer. Cancer Res; 76(21); 6299-310. ©2016 AACR.

The Emerging Function and Mechanism of ceRNAs in Cancer.

Complex diseases, such as cancer, are often associated with aberrant gene expression at both the transcriptional and post-transcriptional level. Over the past several years, competing endogenous RNAs (ceRNAs) have emerged as an important class of post-transcriptional regulators that alter gene expression through a miRNA-mediated mechanism. Recent studies in both solid tumors and hematopoietic malignancies showed that ceRNAs have significant roles in cancer pathogenesis by altering the expression of key tumorigenic or tumor-suppressive genes. Characterizing the identity, function, and mechanism of the ceRNAs will not only further our fundamental understanding of RNA-mediated cancer pathogenesis, but may also shed light on the development of new RNA-based therapeutic strategies for treating cancer.

Microtubule-associated protein 1 light chain 3 interacts with and contributes to growth inhibiting effect of PML.

Previously we reported that the expression of promyelocytic leukemia (PML)-retinoic acid receptor alpha (RARα) fusion gene, which is caused by specific translocation (15;17) in acute promyelocytic leukemia, can enhance constitutive autophagic activity in leukemic and nonleukemic cells, and PML overexpression can sequestrate part of microtubule-associated protein light chain 3 (LC3) protein in PML nuclear bodies, suggesting that LC3 protein also distributes into nuclei although it is currently thought to function primarily in the cytoplasm, the site of autophagosomal formation. However, its potential significance of nucleoplasmic localizations remains greatly elusive. Here we demonstrate that PML interacts with LC3 in a cell type-independent manner as assessed by Co-IP assay and co-localization observation. Overexpressed PML significantly coprecipitates with endogenous and nuclear LC3 protein. Furthermore, a fraction of endogenous PML protein is found to be co-localized with LC3 protein under steady state condition, which is further enhanced by IFNα induction, indicating that PML up-regulation potentiates this interaction. Additionally, DsRed-PML associates with EGFP-LC3 during telophase and G1 phase but not in metaphase and anaphase. Two potential LC3-interacting region (LIR) motifs in PML are required for interaction of PML with LC3 while this association is independent of autophagic activity. Finally, we show that interaction between PML and LC3 contributes to cell growth inhibition function of PML. Considering that PML is an important tumor suppressor, we propose that nuclear portion of LC3 protein may associate with PML to control cell growth for prevention and inhibition of cancer occurrence and development.

Chloroquine enhances cobalt chloride-induced leukemic cell differentiation via the suppression of autophagy at the late phase.

We previously reported that moderate hypoxia and hypoxia-mimetic agents including cobalt chloride (CoCl(2)) induce differentiation of human acute myeloid leukemia (AML) cells through hypoxia-inducible factor-1 α (HIF-1 α), which interacts with and enhances transcriptional activity of CCAAT-enhancer binding factor alpha and Runx1/AML1, two important transcriptional factors for hematopoietic cell differentiation. Here, we show that autophagy inhibitor chloroquine (CQ) increases HIF-1 α accumulation, thus potentiating CoCl(2)-induced growth arrest and differentiation of leukemic cells. Furthermore, the increased effect of CQ on differentiation induction is dependent of the inhibition of autophagosome maturation and degradation, since this sensitization could be mimicked by the suppression of expression of both lysosome-associated membrane proteins 1 and 2 (LAMP1 and LAMP2). These findings not only provide the evidence that CQ is a sensitizer for CoCl(2)-induced differentiation of leukemic cells but also possibly propose the new therapeutic strategy for differentiation induction of AML.

PML-RARα enhances constitutive autophagic activity through inhibiting the Akt/mTOR pathway.

Autophagy is a highly conserved, closely regulated homeostatic cellular activity that allows for the bulk degradation of long-lived proteins and cytoplasmic organelles. Its roles in cancer initiation and progression and in determining the response of tumor cells to anticancer therapy are complicated, and only limited investigation has been conducted on the potential significance of autophagy in the pathogenesis and therapeutic response of acute myeloid leukemia. Here we demonstrate that the inducible or transfected expression of the acute promyelocytic leukemia (APL)-specific PML-RARα, but not PLZF-RARα or NPM-RARα, fusion protein upregulates constitutive autophagy activation in leukemic and nonleukemic cells, as evaluated by hallmarks for autophagy including transmission electron microscopy. The significant increase in autophagic activity is also found in the leukemic cells-infiltrated bone marrow and spleen from PML-RARα-transplanted leukemic mice. The autophagy inhibitor 3-methyladenine significantly abrogates the autophagic events upregulated by PML-RARα, while the autophagic flux assay reveals that the fusion protein induces autophagy by increasing the on-rate of autophagic sequestration. Furthermore, this modulation of autophagy by PML-RARα is possibly mediated by a decreased activation of the Akt/mTOR pathway. Finally, we also show that autophagy contributes to the anti-apoptotic function of the PML-RARα protein. Given the critical role of the PML-RARα oncoprotein in APL pathogenesis, this study suggests an important role of autophagy in the development and treatment of this disease.