Naïve and primed pluripotent states retain distinct molecular properties, yet limited knowledge exists on how their state transitions are regulated. Here, we identify Mettl3, an N(6)-methyladenosine (m(6)A) transferase, as a regulator for terminating murine naïve pluripotency. Mettl3 knockout preimplantation epiblasts and naïve embryonic stem cells are depleted for m(6)A in mRNAs, yet are viable. However, they fail to adequately terminate their naïve state and, subsequently, undergo aberrant and restricted lineage priming at the postimplantation stage, which leads to early embryonic lethality. m(6)A predominantly and directly reduces mRNA stability, including that of key naïve pluripotency-promoting transcripts. This study highlights a critical role for an mRNA epigenetic modification in vivo and identifies regulatory modules that functionally influence naïve and primed pluripotency in an opposing manner.

Ribosome profiling suggests that ribosomes occupy many regions of the transcriptome thought to be noncoding, including 5' UTRs and long noncoding RNAs (lncRNAs). Apparent ribosome footprints outside of protein-coding regions raise the possibility of artifacts unrelated to translation, particularly when they occupy multiple, overlapping open reading frames (ORFs). Here, we show hallmarks of translation in these footprints: copurification with the large ribosomal subunit, response to drugs targeting elongation, trinucleotide periodicity, and initiation at early AUGs. We develop a metric for distinguishing between 80S footprints and nonribosomal sources using footprint size distributions, which validates the vast majority of footprints outside of coding regions. We present evidence for polypeptide production beyond annotated genes, including the induction of immune responses following human cytomegalovirus (HCMV) infection. Translation is pervasive on cytosolic transcripts outside of conserved reading frames, and direct detection of this expanded universe of translated products enables efforts at understanding how cells manage and exploit its consequences.

Viruses are by definition fully dependent on the cellular translation machinery, and develop diverse mechanisms to co-opt this machinery for their own benefit. Unlike many viruses, human cytomegalovirus (HCMV) does suppress the host translation machinery, and the extent to which translation machinery contributes to the overall pattern of viral replication and pathogenesis remains elusive. Here, we combine RNA sequencing and ribosomal profiling analyses to systematically address this question. By simultaneously examining the changes in transcription and translation along HCMV infection, we uncover extensive transcriptional control that dominates the response to infection, but also diverse and dynamic translational regulation for subsets of host genes. We were also able to show that, at late time points in infection, translation of viral mRNAs is higher than that of cellular mRNAs. Lastly, integration of our translation measurements with recent measurements of protein abundance enabled comprehensive identification of dozens of host proteins that are targeted for degradation during HCMV infection. Since targeted degradation indicates a strong biological importance, this approach should be applicable for discovering central host functions during viral infection. Our work provides a framework for studying the contribution of transcription, translation and degradation during infection with any virus.

Primary infection with human cytomegalovirus (HCMV) results in a lifelong infection due to its ability to establish latent infection, with one characterized viral reservoir being hematopoietic cells. Although reactivation from latency causes serious disease in immunocompromised individuals, our molecular understanding of latency is limited. Here, we delineate viral gene expression during natural HCMV persistent infection by analyzing the massive transcriptome RNA sequencing (RNA-seq) atlas generated by the Genotype-Tissue Expression (GTEx) project. This systematic analysis reveals that HCMV persistence in vivo is prevalent in diverse tissues. Notably, we find only viral transcripts that resemble gene expression during various stages of lytic infection with no evidence of any highly restricted latency-associated viral gene expression program. To further define the transcriptional landscape during HCMV latent infection, we also used single-cell RNA-seq and a tractable experimental latency model. In contrast to some current views on latency, we also find no evidence for any highly restricted latency-associated viral gene expression program. Instead, we reveal that latency-associated gene expression largely mirrors a late lytic viral program, albeit at much lower levels of expression. Overall, our work has the potential to revolutionize our understanding of HCMV persistence and suggests that latency is governed mainly by quantitative changes, with a limited number of qualitative changes, in viral gene expression.IMPORTANCE Human cytomegalovirus is a prevalent pathogen, infecting most of the population worldwide and establishing lifelong latency in its hosts. Although reactivation from latency causes significant morbidity and mortality in immunocompromised hosts, our molecular understanding of the latent state remains limited. Here, we examine the viral gene expression during natural and experimental latent HCMV infection on a transcriptome-wide level. In contrast to the classical views on herpesvirus latency, we find no evidence for a restricted latency-associated viral gene expression program. Instead, we reveal that latency gene expression largely resembles a late lytic viral profile, albeit at much lower levels of expression. Taken together, our data transform the current view of HCMV persistence and suggest that latency is mainly governed by quantitative rather than qualitative changes in viral gene expression.

Herpesviruses undergo life-long latent infection which can be life-threatening in the immunocompromised. Models of latency and reactivation of human cytomegalovirus (HCMV) include primary myeloid cells, cells known to be important for HCMV latent carriage and reactivation in vivo. However, primary cells are limited in availability, and difficult to culture and to genetically modify; all of which have hampered our ability to fully understand virus/host interactions of this persistent human pathogen. We have now used iPSCs to develop a model cell system to study HCMV latency and reactivation in different cell types after their differentiation down the myeloid lineage. Our results show that iPSCs can effectively mimic HCMV latency/reactivation in primary myeloid cells, allowing molecular interrogations of the viral latent/lytic switch. This model may also be suitable for analysis of other viruses, such as HIV and Zika, which also infect cells of the myeloid lineage.

Human herpesvirus-6 (HHV-6) A and B are ubiquitous betaherpesviruses, infecting the majority of the human population. They encompass large genomes and our understanding of their protein coding potential is far from complete. Here, we employ ribosome-profiling and systematic transcript-analysis to experimentally define HHV-6 translation products. We identify hundreds of new open reading frames (ORFs), including upstream ORFs (uORFs) and internal ORFs (iORFs), generating a complete unbiased atlas of HHV-6 proteome. By integrating systematic data from the prototypic betaherpesvirus, human cytomegalovirus, we uncover numerous uORFs and iORFs conserved across betaherpesviruses and we show uORFs are enriched in late viral genes. We identified three highly abundant HHV-6 encoded long non-coding RNAs, one of which generates a non-polyadenylated stable intron appearing to be a conserved feature of betaherpesviruses. Overall, our work reveals the complexity of HHV-6 genomes and highlights novel features conserved between betaherpesviruses, providing a rich resource for future functional studies.

mRNAs carry two layers of information, the genetic code and the information that dictates their post-transcriptional fate. The latter function relies on a complex interplay between cis-elements and trans-regulators, and unbiased identification of these elements is still challenging. To identify cis-elements that control gene expression, we use dimethyl sulfate (DMS) mutational profiling with sequencing and map changes in mRNA secondary structure following viral infection. Our dynamic structural data reveal a major role for ribosomes in unwinding secondary structures, which is further supported by the relationship we uncover between structure and translation efficiency. Moreover, our analysis revealed dozens of regions in viral and cellular mRNAs that exhibit changes in secondary structure. In-depth analysis of these regions reveals cis-elements in 3' UTRs that regulate mRNA stability and elements within coding sequences that control translation. Overall, our study demonstrates how mapping dynamic changes in mRNA structure allows unbiased identification of functional regulatory elements.

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) is the cause of the ongoing coronavirus disease 19 (COVID-19) pandemic. Despite its urgency, we still do not fully understand the molecular basis of SARS-CoV-2 pathogenesis and its ability to antagonize innate immune responses. SARS-CoV-2 leads to shutoff of cellular protein synthesis and over-expression of nsp1, a central shutoff factor in coronaviruses, inhibits cellular gene translation. However, the diverse molecular mechanisms nsp1 employs as well as its functional importance in infection are still unresolved. By overexpressing various nsp1 mutants and generating a SARS-CoV-2 mutant in which nsp1 does not bind ribosomes, we untangle the effects of nsp1. We uncover that nsp1, through inhibition of translation and induction of mRNA degradation, is the main driver of host shutoff during SARS-CoV-2 infection. Furthermore, we find the propagation of nsp1 mutant virus is inhibited specifically in cells with intact interferon (IFN) response as well as in-vivo , in infected hamsters, and this attenuation is associated with stronger induction of type I IFN response. This illustrates that nsp1 shutoff activity has an essential role mainly in counteracting the IFN response. Overall, our results reveal the multifaceted approach nsp1 uses to shut off cellular protein synthesis and uncover the central role it plays in SARS-CoV-2 pathogenesis, explicitly through blockage of the IFN response.

Definition of functional genomic elements is one of the greater challenges of the genomic era. Traditionally, putative short open reading frames (sORFs) coding for less than 100 amino acids were disregarded due to computational and experimental limitations; however, it has become clear over the past several years that translation of sORFs is pervasive and serves diverse functions. The development of ribosome profiling, allowing identification of translated sequences genome wide, revealed wide spread, previously unidentified translation events. New computational methodologies as well as improved mass spectrometry approaches also contributed to the task of annotating translated sORFs in different organisms. Viruses are of special interest due to the selective pressure on their genome size, their rapid and confining evolution, and the potential contribution of novel peptides to the host immune response. Indeed, many functional viral sORFs were characterized to date, and ribosome profiling analyses suggest that this may be the tip of the iceberg. Our computational analyses of sORFs identified by ribosome profiling in DNA viruses demonstrate that they may be enriched in specific features implying that at least some of them are functional. Combination of systematic genome editing strategies with synthetic tagging will take us into the next step-elucidation of the biological relevance and function of this intriguing class of molecules.

The mechanism of action of natural killer (NK) cells in type 1 diabetes is still unknown. Here we show that the activating receptor NKp46 recognizes mouse and human ligands on pancreatic beta cells. NK cells appeared in the pancreas when insulitis progressed to type 1 diabetes, and NKp46 engagement by beta cells led to degranulation of NK cells. NKp46-deficient mice had less development of type 1 diabetes induced by injection of a low dose of streptozotocin. Injection of soluble NKp46 proteins into nonobese diabetic mice during the early phase of insulitis and the prediabetic stage prevented the development of type 1 diabetes. Our findings demonstrate that NKp46 is essential for the development of type 1 diabetes and highlight potential new therapeutic modalities for this disease.

Host shutoff is a common strategy used by viruses to repress cellular mRNA translation and concomitantly allow the efficient translation of viral mRNAs. Here we use RNA-sequencing and ribosome profiling to explore the mechanisms that are being utilized by the Influenza A virus (IAV) to induce host shutoff. We show that viral transcripts are not preferentially translated and instead the decline in cellular protein synthesis is mediated by viral takeover on the mRNA pool. Our measurements also uncover strong variability in the levels of cellular transcripts reduction, revealing that short transcripts are less affected by IAV. Interestingly, these mRNAs that are refractory to IAV infection are enriched in cell maintenance processes such as oxidative phosphorylation. Furthermore, we show that the continuous oxidative phosphorylation activity is important for viral propagation. Our results advance our understanding of IAV-induced shutoff, and suggest a mechanism that facilitates the translation of genes with important housekeeping functions.

Thousands of long noncoding RNA (lncRNA) genes are encoded in the human genome, and hundreds of them are evolutionarily conserved, but their functions and modes of action remain largely obscure. Particularly enigmatic lncRNAs are those that are exported to the cytoplasm, including NORAD-an abundant and highly conserved cytoplasmic lncRNA. Here we show that most of the sequence of NORAD is comprised of repetitive units that together contain at least 17 functional binding sites for the two mammalian Pumilio homologues. Through binding to PUM1 and PUM2, NORAD modulates the mRNA levels of their targets, which are enriched for genes involved in chromosome segregation during cell division. Our results suggest that some cytoplasmic lncRNAs function by modulating the activities of RNA-binding proteins, an activity which positions them at key junctions of cellular signalling pathways.

Death associated protein 5 (DAP5/eIF4G2/NAT1) is a member of the eIF4G translation initiation factors that has been shown to mediate noncanonical and/or cap-independent translation. It is essential for embryonic development and for differentiation of embryonic stem cells (ESCs), specifically its ability to drive translation of specific target mRNAs. In order to expand the repertoire of DAP5 target mRNAs, we compared ribosome profiles in control and DAP5 knockdown (KD) human ESCs (hESCs) to identify mRNAs with decreased ribosomal occupancy upon DAP5 silencing. A cohort of 68 genes showed decreased translation efficiency in DAP5 KD cells. Mass spectrometry confirmed decreased protein abundance of a significant portion of these targets. Among these was KMT2D, a histone methylase previously shown to be essential for ESC differentiation and embryonic development. We found that nearly half of the cohort of DAP5 target mRNAs displaying reduced translation efficiency of their main coding sequences upon DAP5 KD contained upstream open reading frames (uORFs) that are actively translated independently of DAP5. This is consistent with previously suggested mechanisms by which DAP5 mediates leaky scanning through uORFs and/or reinitiation at the main coding sequence. Crosslinking protein-RNA immunoprecipitation experiments indicated that a significant subset of DAP5 mRNA targets bound DAP5, indicating that direct binding between DAP5 protein and its target mRNAs is a frequent but not absolute requirement for DAP5-dependent translation of the main coding sequence. Thus, we have extended DAP5's function in translation of specific mRNAs in hESCs by a mechanism allowing translation of the main coding sequence following upstream translation of short ORFs.

The B7 family member, B7H6, is a ligand for the natural killer cell receptor NKp30. B7H6 is hardly expressed on normal tissues, but undergoes upregulation on different types of tumors, implicating it as an attractive target for cancer immunotherapy. The molecular mechanisms that control B7H6 expression are poorly understood. We report that in contrast to other NK cell ligands, endoplasmic reticulum (ER) stress upregulates B7H6 mRNA levels and surface expression. B7H6 induction by ER stress requires protein kinase R-like ER kinase (PERK), one of the three canonical sensors of the unfolded protein response. PERK phosphorylates eIF2α, which regulates protein synthesis and gene expression. Because eIF2α is phosphorylated by several kinases following different stress conditions, the program downstream to eIF2α phosphorylation is called the integrated stress response (ISR). Several drugs were reported to promote the ISR. Nelfinavir and lopinavir, two clinically approved HIV protease inhibitors, promote eIF2α phosphorylation by different mechanisms. We show that nelfinavir and lopinavir sustainably instigate B7H6 expression at their pharmacologically relevant concentrations. As such, ER stress and ISR conditions sensitize melanoma targets to CAR-T cells directed against B7H6. Our study highlights a novel mechanism to induce B7H6 expression and suggests a pharmacological approach to improve B7H6-directed immunotherapy. KEY MESSAGES: B7H6 is induced by ER stress in a PERK-dependent mechanism. Induction of B7H6 is obtained pharmacologically by HIV protease inhibitors. Exposure of tumor cells to the HIV protease inhibitor nelfinavir improves the recognition by B7H6-directed CAR-T.

Human cytomegalovirus (HCMV) is a widespread pathogen establishing a latent infection in its host. HCMV reactivation is a major health burden in immunocompromised individuals, and is a major cause of morbidity and mortality following hematopoietic stem cell transplantation (HSCT). Here we determined HCMV genomic levels using droplet digital PCR in different peripheral blood mononuclear cell (PBMC) populations in HCMV reactivating HSCT patients. This high sensitivity approach revealed that all PBMC populations harbored extremely low levels of viral DNA at the peak of HCMV DNAemia. Transcriptomic analysis of PBMCs from high-DNAemia samples revealed elevated expression of genes typical of HCMV specific T cells, while regulatory T cell enhancers as well as additional genes related to immune response were downregulated. Viral transcript levels in these samples were extremely low, but remarkably, the detected transcripts were mainly immediate early viral genes. Overall, our data indicate that HCMV DNAemia is associated with distinct signatures of immune response in the blood compartment, however it is not necessarily accompanied by substantial infection of PBMCs and the residual infected PBMCs are not productively infected.

MicroRNAs (miRNAs) exert a broad influence over gene expression by directing effector activities that impinge on translation and stability of mRNAs. We recently discovered that the cap-binding protein 4EHP is a key component of the mammalian miRNA-Induced Silencing Complex (miRISC), which mediates gene silencing. However, little is known about the mRNA repertoire that is controlled by the 4EHP/miRNA mechanism or its biological importance. Here, using ribosome profiling, we identify a subset of mRNAs that are translationally controlled by 4EHP. We show that the Dusp6 mRNA, which encodes an ERK1/2 phosphatase, is translationally repressed by 4EHP and a specific miRNA, miR-145. This promotes ERK1/2 phosphorylation, resulting in augmented cell growth and reduced apoptosis. Our findings thus empirically define the integral role of translational repression in miRNA-induced gene silencing and reveal a critical function for this process in the control of the ERK signaling cascade in mammalian cells.

The urea cycle (UC) is the main pathway by which mammals dispose of waste nitrogen. We find that specific alterations in the expression of most UC enzymes occur in many tumors, leading to a general metabolic hallmark termed "UC dysregulation" (UCD). UCD elicits nitrogen diversion toward carbamoyl-phosphate synthetase2, aspartate transcarbamylase, and dihydrooratase (CAD) activation and enhances pyrimidine synthesis, resulting in detectable changes in nitrogen metabolites in both patient tumors and their bio-fluids. The accompanying excess of pyrimidine versus purine nucleotides results in a genomic signature consisting of transversion mutations at the DNA, RNA, and protein levels. This mutational bias is associated with increased numbers of hydrophobic tumor antigens and a better response to immune checkpoint inhibitors independent of mutational load. Taken together, our findings demonstrate that UCD is a common feature of tumors that profoundly affects carcinogenesis, mutagenesis, and immunotherapy response.

Adenosine deaminase acting on RNA 1 (ADAR1) is the master RNA editor, catalyzing the deamination of adenosine to inosine. RNA editing is vital for preventing abnormal activation of cytosolic nucleic acid sensing pathways by self-double-stranded RNAs. Here we determine, by parallel analysis of RNA secondary structure sequencing (PARS-seq), the global RNA secondary structure changes in ADAR1 deficient cells. Surprisingly, ADAR1 silencing resulted in a lower global double-stranded to single-stranded RNA ratio, suggesting that A-to-I editing can stabilize a large subset of imperfect RNA duplexes. The duplexes destabilized by editing are composed of vastly complementary inverted Alus found in untranslated regions of genes performing vital biological processes, including housekeeping functions and type-I interferon responses. They are predominantly cytoplasmic and generally demonstrate higher ribosomal occupancy. Our findings imply that the editing effect on RNA secondary structure is context dependent and underline the intricate regulatory role of ADAR1 on global RNA secondary structure.

Human cytomegalovirus (HCMV) causes diseases in individuals with immature or compromised immunity. To evade immune control, HCMV evolved numerous antagonists targeting the interferon system at multiple levels. By comparative analysis of naturally arising variants of the most widely studied HCMV strain, AD169, and a panel of targeted mutants, we uncover the UL145 gene as indispensable for STAT2 downregulation. Ribosome profiling confirms the translation of the canonical pUL145 protein (pUL145-Long) and newly identifies a shorter isoform (pUL145-Short). Both isoforms recruit DDB1-containing ubiquitin ligases to induce proteasomal degradation of STAT2. An alanine-scanning mutagenesis discloses the DDB1 interaction motif of pUL145 that resembles the DDB1-binding interface of cellular substrate receptors of DDB1-containing ubiquitin ligases. Thus, pUL145 constitutes a viral DDB1-cullin-associated factor (vDCAF), which mimics cellular DCAFs to exploit the ubiquitin-proteasome system to impede antiviral immunity. Notably, the viral exploitation of the cullins can be targeted to restore the efficacy of the host immune response.

The global spread of SARS-CoV-2 led to major economic and health challenges worldwide. Revealing host genes essential for infection by multiple variants of SARS-CoV-2 can provide insights into the virus pathogenesis, and facilitate the development of novel therapeutics. Here, employing a genome-scale CRISPR screen, we provide a comprehensive data-set of cellular factors that are exploited by wild type SARS-CoV-2 as well as two additional recently emerged variants of concerns (VOCs), Alpha and Beta. We identified several host factors critical for SARS-CoV-2 infection, including various components belonging to the Clathrin-dependent transport pathway, ubiquitination, Heparan sulfate biogenesis and host phosphatidylglycerol biosynthesis. Comparative analysis of the different VOCs revealed the host factors KREMEN2 and SETDB1 as potential unique candidates required only to the Alpha variant. Furthermore, the analysis identified GATA6, a zinc finger transcription factor, as an essential proviral gene for all variants inspected. We show that GATA6 directly regulates ACE2 transcription and accordingly, is critical for SARS-CoV-2 cell entry. Analysis of clinical samples collected from SARS-CoV-2 infected individuals shows elevated levels of GATA6, suggesting a role in COVID-19 pathogenesis. Finally, pharmacological inhibition of GATA6 resulted in down-modulation of ACE2 and inhibition of viral infectivity. Overall, we show GATA6 may represent a target for the development of anti-SARS-CoV-2 therapeutic strategies and reaffirm the value of the CRISPR loss-of-function screens in providing a list of potential new targets for therapeutic interventions.