Profiling SARS-CoV-2 HLA-I peptidome reveals T cell epitopes from out-of-frame ORFs.

T cell-mediated immunity plays an important role in controlling SARS-CoV-2 infection, but the repertoire of naturally processed and presented viral epitopes on class I human leukocyte antigen (HLA-I) remains uncharacterized. Here, we report the first HLA-I immunopeptidome of SARS-CoV-2 in two cell lines at different times post infection using mass spectrometry. We found HLA-I peptides derived not only from canonical open reading frames (ORFs) but also from internal out-of-frame ORFs in spike and nucleocapsid not captured by current vaccines. Some peptides from out-of-frame ORFs elicited T cell responses in a humanized mouse model and individuals with COVID-19 that exceeded responses to canonical peptides, including some of the strongest epitopes reported to date. Whole-proteome analysis of infected cells revealed that early expressed viral proteins contribute more to HLA-I presentation and immunogenicity. These biological insights, as well as the discovery of out-of-frame ORF epitopes, will facilitate selection of peptides for immune monitoring and vaccine development.

m6A modification controls the innate immune response to infection by targeting type I interferons.

N6-methyladenosine (m6A) is the most common mRNA modification. Recent studies have revealed that depletion of m6A machinery leads to alterations in the propagation of diverse viruses. These effects were proposed to be mediated through dysregulated methylation of viral RNA. Here we show that following viral infection or stimulation of cells with an inactivated virus, deletion of the m6A 'writer' METTL3 or 'reader' YTHDF2 led to an increase in the induction of interferon-stimulated genes. Consequently, propagation of different viruses was suppressed in an interferon-signaling-dependent manner. Significantly, the mRNA of IFNB, the gene encoding the main cytokine that drives the type I interferon response, was m6A modified and was stabilized following repression of METTL3 or YTHDF2. Furthermore, we show that m6A-mediated regulation of interferon genes was conserved in mice. Together, our findings uncover the role m6A serves as a negative regulator of interferon response by dictating the fast turnover of interferon mRNAs and consequently facilitating viral propagation.

Publisher Correction: m6A modification controls the innate immune response to infection by targeting type I interferons.

In the version of this article initially published, the penultimate sentence of the abstract included a typographical error ('cxgenes'). The correct word is 'genes'. The error has been corrected in the HTML and PDF version of the article.

A virally encoded high-resolution screen of cytomegalovirus dependencies.

Genetic screens have transformed our ability to interrogate cellular factor requirements for viral infections1,2, but most current approaches are limited in their sensitivity, biased towards early stages of infection and provide only simplistic phenotypic information that is often based on survival of infected cells2-4. Here, by engineering human cytomegalovirus to express single guide RNA libraries directly from the viral genome, we developed virus-encoded CRISPR-based direct readout screening (VECOS), a sensitive, versatile, viral-centric approach that enables profiling of different stages of viral infection in a pooled format. Using this approach, we identified hundreds of host dependency and restriction factors and quantified their direct effects on viral genome replication, viral particle secretion and infectiousness of secreted particles, providing a multi-dimensional perspective on virus-host interactions. These high-resolution measurements reveal that perturbations altering late stages in the life cycle of human cytomegalovirus (HCMV) mostly regulate viral particle quality rather than quantity, establishing correct virion assembly as a critical stage that is heavily reliant on virus-host interactions. Overall, VECOS facilitates systematic high-resolution dissection of the role of human proteins during the infection cycle, providing a roadmap for in-depth study of host-herpesvirus interactions.

Context-dependent functional compensation between Ythdf m6A reader proteins.

The N6-methyladenosine (m6A) modification is the most prevalent post-transcriptional mRNA modification, regulating mRNA decay and splicing. It plays a major role during normal development, differentiation, and disease progression. The modification is regulated by a set of writer, eraser, and reader proteins. The YTH domain family of proteins consists of three homologous m6A-binding proteins, Ythdf1, Ythdf2, and Ythdf3, which were suggested to have different cellular functions. However, their sequence similarity and their tendency to bind the same targets suggest that they may have overlapping roles. We systematically knocked out (KO) the Mettl3 writer, each of the Ythdf readers, and the three readers together (triple-KO). We then estimated the effect in vivo in mouse gametogenesis, postnatal viability, and in vitro in mouse embryonic stem cells (mESCs). In gametogenesis, Mettl3-KO severity is increased as the deletion occurs earlier in the process, and Ythdf2 has a dominant role that cannot be compensated by Ythdf1 or Ythdf3, due to differences in readers' expression pattern across different cell types, both in quantity and in spatial location. Knocking out the three readers together and systematically testing viable offspring genotypes revealed a redundancy in the readers' role during early development that is Ythdf1/2/3 gene dosage-dependent. Finally, in mESCs there is compensation between the three Ythdf reader proteins, since the resistance to differentiate and the significant effect on mRNA decay occur only in the triple-KO cells and not in the single KOs. Thus, we suggest a new model for the Ythdf readers function, in which there is profound dosage-dependent redundancy when all three readers are equivalently coexpressed in the same cell types.

The coding capacity of SARS-CoV-2.

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) is the cause of the ongoing coronavirus disease 2019 (COVID-19) pandemic1. To understand the pathogenicity and antigenic potential of SARS-CoV-2 and to develop therapeutic tools, it is essential to profile the full repertoire of its expressed proteins. The current map of SARS-CoV-2 coding capacity is based on computational predictions and relies on homology with other coronaviruses. As the protein complement varies among coronaviruses, especially in regard to the variety of accessory proteins, it is crucial to characterize the specific range of SARS-CoV-2 proteins in an unbiased and open-ended manner. Here, using a suite of ribosome-profiling techniques2-4, we present a high-resolution map of coding regions in the SARS-CoV-2 genome, which enables us to accurately quantify the expression of canonical viral open reading frames (ORFs) and to identify 23 unannotated viral ORFs. These ORFs include upstream ORFs that are likely to have a regulatory role, several in-frame internal ORFs within existing ORFs, resulting in N-terminally truncated products, as well as internal out-of-frame ORFs, which generate novel polypeptides. We further show that viral mRNAs are not translated more efficiently than host mRNAs; instead, virus translation dominates host translation because of the high levels of viral transcripts. Our work provides a resource that will form the basis of future functional studies.

SARS-CoV-2 uses a multipronged strategy to impede host protein synthesis.

The coronavirus SARS-CoV-2 is the cause of the ongoing pandemic of COVID-191. Coronaviruses have developed a variety of mechanisms to repress host mRNA translation to allow the translation of viral mRNA, and concomitantly block the cellular innate immune response2,3. Although several different proteins of SARS-CoV-2 have previously been implicated in shutting off host expression4-7, a comprehensive picture of the effects of SARS-CoV-2 infection on cellular gene expression is lacking. Here we combine RNA sequencing, ribosome profiling and metabolic labelling of newly synthesized RNA to comprehensively define the mechanisms that are used by SARS-CoV-2 to shut off cellular protein synthesis. We show that infection leads to a global reduction in translation, but that viral transcripts are not preferentially translated. Instead, we find that infection leads to the accelerated degradation of cytosolic cellular mRNAs, which facilitates viral takeover of the mRNA pool in infected cells. We reveal that the translation of transcripts that are induced in response to infection (including innate immune genes) is impaired. We demonstrate this impairment is probably mediated by inhibition of nuclear mRNA export, which prevents newly transcribed cellular mRNA from accessing ribosomes. Overall, our results uncover a multipronged strategy that is used by SARS-CoV-2 to take over the translation machinery and to suppress host defences.

Virus-Induced Changes in mRNA Secondary Structure Uncover cis-Regulatory Elements that Directly Control Gene Expression.

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.

Dynamic changes in tRNA modifications and abundance during T cell activation.

The tRNA pool determines the efficiency, throughput, and accuracy of translation. Previous studies have identified dynamic changes in the tRNA (transfer RNA) supply and mRNA (messenger RNA) demand during cancerous proliferation. Yet dynamic changes may also occur during physiologically normal proliferation, and these are less well characterized. We examined the tRNA and mRNA pools of T cells during their vigorous proliferation and differentiation upon triggering their antigen receptor. We observed a global signature of switch in demand for codons at the early proliferation phase of the response, accompanied by corresponding changes in tRNA expression levels. In the later phase, upon differentiation, the response of the tRNA pool relaxed back to the basal level, potentially restraining excessive proliferation. Sequencing of tRNAs allowed us to evaluate their diverse base-modifications. We found that two types of tRNA modifications, wybutosine and ms2t6A, are reduced dramatically during T cell activation. These modifications occur in the anticodon loops of two tRNAs that decode "slippery codons," which are prone to ribosomal frameshifting. Attenuation of these frameshift-protective modifications is expected to increase the potential for proteome-wide frameshifting during T cell proliferation. Indeed, human cell lines deleted of a wybutosine writer showed increased ribosomal frameshifting, as detected with an HIV gag-pol frameshifting site reporter. These results may explain HIV's specific tropism toward proliferating T cells since it requires ribosomal frameshift exactly on the corresponding codon for infection. The changes in tRNA expression and modifications uncover a layer of translation regulation during T cell proliferation and expose a potential tradeoff between cellular growth and translation fidelity.

Bromodomain proteins regulate human cytomegalovirus latency and reactivation allowing epigenetic therapeutic intervention.

Reactivation of human cytomegalovirus (HCMV) from latency is a major health consideration for recipients of stem-cell and solid organ transplantations. With over 200,000 transplants taking place globally per annum, virus reactivation can occur in more than 50% of cases leading to loss of grafts as well as serious morbidity and even mortality. Here, we present the most extensive screening to date of epigenetic inhibitors on HCMV latently infected cells and find that histone deacetylase inhibitors (HDACis) and bromodomain inhibitors are broadly effective at inducing virus immediate early gene expression. However, while HDACis, such as myeloid-selective CHR-4487, lead to production of infectious virions, inhibitors of bromodomain (BRD) and extraterminal proteins (I-BETs), including GSK726, restrict full reactivation. Mechanistically, we show that BET proteins (BRDs) are pivotally connected to regulation of HCMV latency and reactivation. Through BRD4 interaction, the transcriptional activator complex P-TEFb (CDK9/CycT1) is sequestered by repressive complexes during HCMV latency. Consequently, I-BETs allow release of P-TEFb and subsequent recruitment to promoters via the superelongation complex (SEC), inducing transcription of HCMV lytic genes encoding immunogenic antigens from otherwise latently infected cells. Surprisingly, this occurs without inducing many viral immunoevasins and, importantly, while also restricting viral DNA replication and full HCMV reactivation. Therefore, this pattern of HCMV transcriptional dysregulation allows effective cytotoxic immune targeting and killing of latently infected cells, thus reducing the latent virus genome load. This approach could be safely used to pre-emptively purge the virus latent reservoir prior to transplantation, thereby reducing HCMV reactivation-related morbidity and mortality.

The m1A landscape on cytosolic and mitochondrial mRNA at single-base resolution.

Modifications on mRNA offer the potential of regulating mRNA fate post-transcriptionally. Recent studies suggested the widespread presence of N1-methyladenosine (m1A), which disrupts Watson-Crick base pairing, at internal sites of mRNAs. These studies lacked the resolution of identifying individual modified bases, and did not identify specific sequence motifs undergoing the modification or an enzymatic machinery catalysing them, rendering it challenging to validate and functionally characterize putative sites. Here we develop an approach that allows the transcriptome-wide mapping of m1A at single-nucleotide resolution. Within the cytosol, m1A is present in a low number of mRNAs, typically at low stoichiometries, and almost invariably in tRNA T-loop-like structures, where it is introduced by the TRMT6/TRMT61A complex. We identify a single m1A site in the mitochondrial ND5 mRNA, catalysed by TRMT10C, with methylation levels that are highly tissue specific and tightly developmentally controlled. m1A leads to translational repression, probably through a mechanism involving ribosomal scanning or translation. Our findings suggest that m1A on mRNA, probably because of its disruptive impact on base pairing, leads to translational repression, and is generally avoided by cells, while revealing one case in mitochondria where tight spatiotemporal control over m1A levels was adopted as a potential means of post-transcriptional regulation.

Human cytomegalovirus long noncoding RNA4.9 regulates viral DNA replication.

Viruses are known for their extremely compact genomes composed almost entirely of protein-coding genes. Nonetheless, four long noncoding RNAs (lncRNAs) are encoded by human cytomegalovirus (HCMV). Although these RNAs accumulate to high levels during lytic infection, their functions remain largely unknown. Here, we show that HCMV-encoded lncRNA4.9 localizes to the viral nuclear replication compartment, and that its depletion restricts viral DNA replication and viral growth. RNA4.9 is transcribed from the HCMV origin of replication (oriLyt) and forms an RNA-DNA hybrid (R-loop) through its G+C-rich 5' end, which may be important for the initiation of viral DNA replication. Furthermore, targeting the RNA4.9 promoter with CRISPR-Cas9 or genetic relocalization of oriLyt leads to reduced levels of the viral single-stranded DNA-binding protein (ssDBP), suggesting that the levels of ssDBP are coupled to the oriLyt activity. We further identified a similar, oriLyt-embedded, G+C-rich lncRNA in murine cytomegalovirus (MCMV). These results indicate that HCMV RNA4.9 plays an important role in regulating viral DNA replication, that the levels of ssDBP are coupled to the oriLyt activity, and that these regulatory features may be conserved among betaherpesviruses.

Rho-Associated Coiled-Coil Kinase 1 Translocates to the Nucleus and Inhibits Human Cytomegalovirus Propagation.

Rho-associated coiled-coil kinase (ROCK) protein is a central kinase that regulates numerous cellular functions, including cellular polarity, motility, proliferation, and apoptosis. Here, we demonstrate that ROCK has antiviral properties, and inhibition of its activity results in enhanced propagation of human cytomegalovirus (HCMV). We show that during HCMV infection, ROCK1 translocates to the nucleus and concentrates in the nucleolus, where it colocalizes with the stress-related chaperone heat shock cognate 71-kDa protein (Hsc70). Gene expression measurements show that inhibition of ROCK activity does not seem to affect the cellular stress response. We demonstrate that inhibition of myosin, one of the central targets of ROCK, also increases HCMV propagation, implying that the antiviral activity of ROCK might be mediated by activation of the actomyosin network. Finally, we demonstrate that inhibition of ROCK results in increased levels of the tegument protein UL32 and of viral DNA in the cytoplasm, suggesting ROCK activity hinders the efficient egress of HCMV particles out of the nucleus. Altogether, our findings illustrate ROCK activity restricts HCMV propagation and suggest this inhibitory effect may be mediated by suppression of capsid egress out of the nucleus.IMPORTANCE ROCK is a central kinase in cells that regulates numerous cellular functions, including cellular polarity, motility, proliferation, and apoptosis. Here we reveal a novel antiviral activity of ROCK during infection with HCMV, a prevalent pathogen infecting most of the population worldwide. We reveal ROCK1 is translocated to the nucleus, where it mainly localizes to the nucleolus. Our findings suggest that ROCK's antiviral activity may be related to activation of the actomyosin network and inhibition of capsid egress out of the nucleus.

Identifying genetic susceptibility to Aspergillus fumigatus infection using collaborative cross mice and RNA-Seq approach.

BACKGROUND: Aspergillus fumigatus (Af) is one of the most ubiquitous fungi and its infection potency is suggested to be strongly controlled by the host genetic background. The aim of this study was to search for candidate genes associated with host susceptibility to Aspergillus fumigatus (Af) using an RNAseq approach in CC lines and hepatic gene expression. METHODS: We studied 31 male mice from 25 CC lines at 8 weeks old; the mice were infected with Af. Liver tissues were extracted from these mice 5 days post-infection, and next-generation RNA-sequencing (RNAseq) was performed. The GENE-E analysis platform was used to generate a clustered heat map matrix. RESULTS: Significant variation in body weight changes between CC lines was observed. Hepatic gene expression revealed 12 top prioritized candidate genes differentially expressed in resistant versus susceptible mice based on body weight changes. Interestingly, three candidate genes are located within genomic intervals of the previously mapped quantitative trait loci (QTL), including Gm16270 and Stox1 on chromosome 10 and Gm11033 on chromosome 8. CONCLUSIONS: Our findings emphasize the CC mouse model's power in fine mapping the genetic components underlying susceptibility towards Af. As a next step, eQTL analysis will be performed for our RNA-Seq data. Suggested candidate genes from our study will be further assessed with a human cohort with aspergillosis.

Molecular characterization of human cytomegalovirus infection with single-cell transcriptomics.

Human cytomegalovirus (HCMV) can result in either productive or non-productive infection, with the latter potentially leading to viral latency. The molecular factors dictating these outcomes are poorly understood. Here we used single-cell transcriptomics to analyse HCMV infection progression in monocytes, which are latently infected, and macrophages, considered to be permissive for productive infection. We show that early viral gene expression levels, specifically of those encoding immediate early proteins IE1 and IE2, are a major factor dictating productive infection. We also revealed that intrinsic, not induced, host cell interferon-stimulated gene expression level is a main determinant of infection outcome. Intrinsic interferon-stimulated gene expression is downregulated with monocyte to macrophage differentiation, partially explaining increased macrophage susceptibility to productive HCMV infection. Furthermore, non-productive macrophages could reactivate, making them potential latent virus reservoirs. Overall, we decipher molecular features underlying HCMV infection outcomes and propose macrophages as a potential HCMV reservoir.

Temporal dynamics of HCMV gene expression in lytic and latent infections.

During productive human cytomegalovirus (HCMV) infection, viral genes are expressed in a coordinated cascade that conventionally relies on the dependencies of viral genes on protein synthesis and viral DNA replication. By contrast, the transcriptional landscape of HCMV latency is poorly understood. Here, we examine viral gene expression dynamics during the establishment of both productive and latent HCMV infections. We redefine HCMV gene expression kinetics during productive infection and reveal that viral gene regulation does not represent a simple sequential cascade; many viral genes are regulated by multiple independent modules. Using our improved gene expression classification combined with transcriptome-wide measurements of the effects of a wide array of epigenetic inhibitors on viral gene expression during latency, we show that a defining feature of latency is the unique repression of immediate-early (IE) genes. Altogether, we recharacterize HCMV gene expression kinetics and reveal governing principles of lytic and latent gene expression.

Parsing the role of NSP1 in SARS-CoV-2 infection.

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) leads to shutoff of protein synthesis, and nsp1, a central shutoff factor in coronaviruses, inhibits cellular mRNA translation. However, the diverse molecular mechanisms employed by nsp1 as well as its functional importance are unresolved. By overexpressing various nsp1 mutants and generating a SARS-CoV-2 mutant, we show that nsp1, through inhibition of translation and induction of mRNA degradation, targets translated cellular mRNA and is the main driver of host shutoff during infection. The propagation of nsp1 mutant virus is inhibited exclusively in cells with intact interferon (IFN) pathway as well as in vivo, in hamsters, and this attenuation is associated with stronger induction of type I IFN response. Therefore, although nsp1's shutoff activity is broad, it plays an essential role, specifically in counteracting the IFN response. Overall, our results reveal the multifaceted approach nsp1 uses to shut off cellular protein synthesis and uncover nsp1's explicit role in blocking the IFN response.

Parsing the role of NSP1 in SARS-CoV-2 infection.

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.

The Human Cytomegalovirus pUL145 Isoforms Act as Viral DDB1-Cullin-Associated Factors to Instruct Host Protein Degradation to Impede Innate Immunity.

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.

Comprehensive annotations of human herpesvirus 6A and 6B genomes reveal novel and conserved genomic features.

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.