A Mouse-Adapted SARS-CoV-2 Induces Acute Lung Injury and Mortality in Standard Laboratory Mice.

The SARS-CoV-2 pandemic has caused extreme human suffering and economic harm. We generated and characterized a new mouse-adapted SARS-CoV-2 virus that captures multiple aspects of severe COVID-19 disease in standard laboratory mice. This SARS-CoV-2 model exhibits the spectrum of morbidity and mortality of COVID-19 disease as well as aspects of host genetics, age, cellular tropisms, elevated Th1 cytokines, and loss of surfactant expression and pulmonary function linked to pathological features of acute lung injury (ALI) and acute respiratory distress syndrome (ARDS). This model can rapidly access existing mouse resources to elucidate the role of host genetics, underlying molecular mechanisms governing SARS-CoV-2 pathogenesis, and the protective or pathogenic immune responses related to disease severity. The model promises to provide a robust platform for studies of ALI and ARDS to evaluate vaccine and antiviral drug performance, including in the most vulnerable populations (i.e., the aged) using standard laboratory mice.

Protein phosphatase 1 suppresses PKR/EIF2α signaling during human cytomegalovirus infection.

Human cytomegalovirus (HCMV) is a ubiquitous pathogen that infects the majority of the world's population. Lytic HCMV replication in immunocompromised individuals or neonates can lead to severe disease in multiple organ systems and even death. The establishment of lytic replication is driven by the first viral proteins expressed upon infection, the immediate early proteins, which play a key role in creating an intracellular environment conducive to virus replication. Two immediate early proteins, the functional orthologs pTRS1 and pIRS1, stimulate immediate early gene expression by suppressing antiviral PKR/eIF2α signaling and enhance the translation of viral mRNAs independent of PKR antagonism. To better understand the molecular functions of pTRS1, we used proximity labeling proteomics to identify proteins that interact with pTRS1 in infected cells. Multiple novel host and viral interactors were identified, including the catalytic subunits of the protein phosphatase 1 (PP1) holoenzyme. Mutations to a PP1 catalytic subunit known to disrupt binding to PP1 regulatory subunits decreased binding to pTRS1. pTRS1 immune complexes contained phosphatase activity, and inhibition of phosphatase activity in transfected or infected cells reversed the ability of pTRS1 to inhibit the antiviral kinase PKR. Depletion of individual PP1 catalytic subunits decreased virus replication and increased the phosphorylation of the PKR substrate eIF2α. Taken together, our data suggest potential novel functions for pTRS1 and define a novel role for PP1 as an antagonist of the antiviral PKR/eIF2α signaling axis during HCMV infection.IMPORTANCEThe human cytomegalovirus (HCMV) pTRS1 and pIRS1 proteins are critical regulators of HCMV replication, both during primary infection and during reactivation from viral latency. Thus, defining the molecular functions of pTRS1/pIRS1 is important for understanding the molecular events controlling HCMV replication and viral disease. These data provide new insights into potential pTRS1 functional roles, providing a starting point for others to understand new features of infected cell biology. Another important result of this study is the finding that specific protein phosphatase 1 (PP1) regulatory subunits are required to suppress PKR/eIF2α signaling, a critical cellular innate immune defense to viral infection. These data lay the groundwork for future efforts to discover therapeutics that disrupt pTRS1 interaction with PP1 allowing cellular defenses to limit HCMV replication and disease.

Mapping the Human Herpesvirus 6B transcriptome.

The "omics" revolution of recent years has simplified the study of RNA transcripts produced during viral infection and under specific defined conditions. In the quest to find new and differentially expressed transcripts during the course of human Herpesvirus 6B (HHV-6B) infection, we made use of large-scale RNA sequencing to analyze the HHV-6B transcriptome during productive infection of human Molt-3 T-cells. Analyses were performed at different time points following infection and specific inhibitors were used to classify the kinetic class of each open reading frame (ORF) reported in the annotated genome of HHV-6B Z29 strain. The initial search focussed on HHV-6B-specific reads matching new HHV-6B transcripts. Differential expression of new HHV-6B transcripts were observed in all samples analyzed. The presence of many of these new HHV-6B transcripts were confirmed by RT-PCR and Sanger sequencing. Many of these transcripts represented new splice variants of previously reported ORFs, including some transcripts that have yet to be defined. Overall, our work demonstrates the diversity and the complexity of the HHV-6B transcriptome.IMPORTANCERNA sequencing (RNA-seq) is an important tool for studying RNA transcripts, particularly during active viral infection. We made use of RNA-seq to study human Herpesvirus 6B (HHV-6B) infection. Using six different time points, we were able to identify the presence of differentially spliced genes at 6, 9, 12, 24, 48 and 72 hours post-infection. Determination of the RNA profiles in the presence of cycloheximide (CHX) or phosphonoacetic acid (PAA) also permitted identification of the kinetic class of each ORF described in the annotated GenBank file. We also identified new spliced transcripts for certain genes and evaluated their relative expression over time. These data and next-generation sequencing (NGS) of the viral DNA have led us to propose a new version of the HHV-6B Z29 GenBank annotated file, without changing ORF names in order to facilitate trace back and correlate our work with previous studies on HHV-6B.

Identification of Quinolinones as Antivirals against Venezuelan Equine Encephalitis Virus.

Venezuelan equine encephalitis virus (VEEV) is a reemerging alphavirus that can cause encephalitis resulting in severe human morbidity and mortality. Using a high-throughput cell-based screen, we identified a quinolinone compound that protected against VEEV-induced cytopathic effects. Analysis of viral replication in cells identified several quinolinone compounds with potent inhibitory activity against vaccine and virulent strains of VEEV. These quinolinones also displayed inhibitory activity against additional alphaviruses, such as Mayaro virus and Ross River virus, although the potency was greatly reduced. Time-of-addition studies indicated that these compounds inhibit the early-to-mid stage of viral replication. Deep sequencing and reverse genetics studies identified two unique resistance mutations in the nsP2 gene (Y102S/C; stalk domain) that conferred VEEV resistance on this chemical series. Moreover, introduction of a K102Y mutation into the nsP2 gene enhanced the sensitivity of chikungunya virus (CHIKV) to this chemical series. Computational modeling of CHIKV and VEEV nsP2 identified a highly probable docking alignment for the quinolinone compounds that require a tyrosine residue at position 102 within the helicase stalk domain. These studies identified a class of compounds with antiviral activity against VEEV and other alphaviruses and provide further evidence that therapeutics targeting nsP2 may be useful against alphavirus infection.

Two Genetic Differences between Closely Related Zika Virus Strains Determine Pathogenic Outcome in Mice.

Recent Zika virus (ZIKV) outbreaks and unexpected clinical manifestations of ZIKV infection have prompted an increase in ZIKV-related research. Here, we identify two strain-specific determinants of ZIKV virulence in mice. We found that strain H/PF/2013 caused 100% lethality in Ifnar1-/- mice, whereas PRVABC59 caused no lethality; both strains caused 100% lethality in Ifnar1-/-Ifngr1-/- double-knockout (DKO) mice. Deep sequencing revealed a high-frequency variant in PRVABC59 not present in H/PF/2013: a G-to-T change at nucleotide 1965 producing a Val-to-Leu substitution at position 330 of the viral envelope (E) protein. We show that the V330 variant is lethal on both virus strain backgrounds, whereas the L330 variant is attenuating only on the PRVABC59 background. These results identify a balanced polymorphism in the E protein that is sufficient to attenuate the PRVABC59 strain but not H/PF/2013. The consensus sequences of H/PF/2013 and PRVABC59 differ by 3 amino acids, but these were not responsible for the difference in virulence between the two strains. H/PF/2013 and PRVABC59 differ by an additional 31 noncoding or silent nucleotide changes. We made a panel of chimeric viruses with identical amino acid sequences but nucleotide sequences derived from H/PF/2013 or PRVABC59. We found that 6 nucleotide differences in the 3' quarter of the H/PF/2013 genome were sufficient to confer virulence in Ifnar1-/- mice. Altogether, our work identifies a large and previously unreported difference in virulence between two commonly used ZIKV strains, in two widely used mouse models of ZIKV pathogenesis (Ifnar1-/- and Ifnar1-/- Ifngr1-/- DKO mice).IMPORTANCE Contemporary ZIKV strains are closely related and often used interchangeably in laboratory research. Here, we identify two strain-specific determinants of ZIKV virulence that are evident in only Ifnar1-/- mice but not Ifnar1-/-Ifngr1-/- DKO mice. These results identify a balanced polymorphism in the E protein that is sufficient to attenuate the PRVABC59 strain but not H/PF/2013. We further identify a second virulence determinant in the H/PF/2013 strain, which is driven by the viral nucleotide sequence but not the amino acid sequence. Altogether, our work identifies a large and previously unreported difference in virulence between two commonly used ZIKV strains, in two widely used mouse models of ZIKV pathogenesis. Our results highlight that even very closely related virus strains can produce significantly different pathogenic phenotypes in common laboratory models.

Using SHAPE-MaP To Model RNA Secondary Structure and Identify 3'UTR Variation in Chikungunya Virus.

Chikungunya virus (CHIKV) is a mosquito-borne alphavirus associated with debilitating arthralgia in humans. RNA secondary structure in the viral genome plays an important role in the lifecycle of alphaviruses; however, the specific role of RNA structure in regulating CHIKV replication is poorly understood. Our previous studies found little conservation in RNA secondary structure between alphaviruses, and this structural divergence creates unique functional structures in specific alphavirus genomes. Therefore, to understand the impact of RNA structure on CHIKV biology, we used SHAPE-MaP to inform the modeling of RNA secondary structure throughout the genome of a CHIKV isolate from the 2013 Caribbean outbreak. We then analyzed regions of the genome with high levels of structural specificity to identify potentially functional RNA secondary structures and identified 23 regions within the CHIKV genome with higher than average structural stability, including four previously identified, functionally important CHIKV RNA structures. We also analyzed the RNA flexibility and secondary structures of multiple 3'UTR variants of CHIKV that are known to affect virus replication in mosquito cells. This analysis found several novel RNA structures within these 3'UTR variants. A duplication in the 3'UTR that enhances viral replication in mosquito cells led to an overall increase in the amount of unstructured RNA in the 3'UTR. This analysis demonstrates that the CHIKV genome contains a number of unique, specific RNA secondary structures and provides a strategy for testing these secondary structures for functional importance in CHIKV replication and pathogenesis.IMPORTANCE Chikungunya virus (CHIKV) is a mosquito-borne RNA virus that causes febrile illness and debilitating arthralgia in humans. CHIKV causes explosive outbreaks but there are no approved therapies to treat or prevent CHIKV infection. The CHIKV genome contains functional RNA secondary structures that are essential for proper virus replication. Since RNA secondary structures have only been defined for a small portion of the CHIKV genome, we used a chemical probing method to define the RNA secondary structures of CHIKV genomic RNA. We identified 23 highly specific structured regions of the genome, and confirmed the functional importance of one structure using mutagenesis. Furthermore, we defined the RNA secondary structure of three CHIKV 3'UTR variants that differ in their ability to replicate in mosquito cells. Our study highlights the complexity of the CHIKV genome and describes new systems for designing compensatory mutations to test the functional relevance of viral RNA secondary structures.

Structural divergence creates new functional features in alphavirus genomes.

Alphaviruses are mosquito-borne pathogens that cause human diseases ranging from debilitating arthritis to lethal encephalitis. Studies with Sindbis virus (SINV), which causes fever, rash, and arthralgia in humans, and Venezuelan equine encephalitis virus (VEEV), which causes encephalitis, have identified RNA structural elements that play key roles in replication and pathogenesis. However, a complete genomic structural profile has not been established for these viruses. We used the structural probing technique SHAPE-MaP to identify structured elements within the SINV and VEEV genomes. Our SHAPE-directed structural models recapitulate known RNA structures, while also identifying novel structural elements, including a new functional element in the nsP1 region of SINV whose disruption causes a defect in infectivity. Although RNA structural elements are important for multiple aspects of alphavirus biology, we found the majority of RNA structures were not conserved between SINV and VEEV. Our data suggest that alphavirus RNA genomes are highly divergent structurally despite similar genomic architecture and sequence conservation; still, RNA structural elements are critical to the viral life cycle. These findings reframe traditional assumptions about RNA structure and evolution: rather than structures being conserved, alphaviruses frequently evolve new structures that may shape interactions with host immune systems or co-evolve with viral proteins.

Functional Analyses of the RsmY and RsmZ Small Noncoding Regulatory RNAs in Pseudomonas aeruginosa.

Pseudomonas aeruginosa is a Gram-negative opportunistic pathogen with distinct acute and chronic virulence phenotypes. Whereas acute virulence is typically associated with expression of a type III secretion system (T3SS), chronic virulence is characterized by biofilm formation. Many of the phenotypes associated with acute and chronic virulence are inversely regulated by RsmA and RsmF. RsmA and RsmF are both members of the CsrA family of RNA-binding proteins and regulate protein synthesis at the posttranscriptional level. RsmA activity is controlled by two small noncoding regulatory RNAs (RsmY and RsmZ). Bioinformatic analyses suggest that RsmY and RsmZ each have 3 or 4 putative RsmA binding sites. Each predicted binding site contains a GGA sequence presented in the loop portion of a stem-loop structure. RsmY and RsmZ regulate RsmA, and possibly RsmF, by sequestering these proteins from target mRNAs. In this study, we used selective 2'-hydroxyl acylation analyzed by primer extension and mutational profiling (SHAPE-MaP) chemistry to determine the secondary structures of RsmY and RsmZ and functional assays to characterize the contribution of each GGA site to RsmY/RsmZ activity. Our data indicate that RsmA has two preferential binding sites on RsmY and RsmZ, while RsmF has one preferential binding site on RsmY and two sites on RsmZ. Despite RsmF and RsmA sharing a common consensus site, RsmF binding properties are more restrictive than those of RsmA.IMPORTANCE CsrA homologs are present in many bacteria. The opportunistic pathogen Pseudomonas aeruginosa uses RsmA and RsmF to inversely regulate factors associated with acute and chronic virulence phenotypes. RsmA has an affinity for RsmY and RsmZ higher than that of RsmF. The goal of this study was to understand the differential binding properties of RsmA and RsmF by using the RsmY and RsmZ regulatory small RNAs (sRNAs) as a model. Mutagenesis of the predicted RsmA/RsmF binding sites on RsmY and RsmZ revealed similarities in the sites required to control RsmA and RsmF activity in vivo Whereas binding by RsmA was relatively tolerant of binding site mutations, RsmF was sensitive to disruption to all but two of the sites, further demonstrating that the requirements for RsmF binding activity in vivo and in vitro are more stringent than those for RsmA.

Mutations in the E2 Glycoprotein and the 3' Untranslated Region Enhance Chikungunya Virus Virulence in Mice.

Chikungunya virus (CHIKV) is a mosquito-transmitted alphavirus that causes debilitating musculoskeletal pain and inflammation and can persist for months to years after acute infection. Although studies of humans and experimentally infected animals suggest that CHIKV infection persists in musculoskeletal tissues, the mechanisms for this remain poorly understood. To evaluate this further, we isolated CHIKV from the serum of persistently infected Rag1 -/- mice at day 28. When inoculated into naive wild-type (WT) mice, this persistently circulating CHIKV strain displayed a capacity for earlier dissemination and greater pathogenicity than the parental virus. Sequence analysis revealed a nonsynonymous mutation in the E2 glycoprotein (E2 K200R) and a deletion within the 3' untranslated region (3'-UTR). The introduction of these changes into the parental virus conferred enhanced virulence in mice, although primary tropism for musculoskeletal tissues was maintained. The E2 K200R mutation was largely responsible for enhanced viral dissemination and pathogenicity, although these effects were augmented by the 3'-UTR deletion. Finally, studies with Irf3/Irf7 -/- and Ifnar1 -/- mice suggest that the E2 K200R mutation enhances viral dissemination from the site of inoculation independently of interferon regulatory factor 3 (IRF3)-, IRF7-, and IFNAR1-mediated responses. As our findings reveal viral determinants of CHIKV dissemination and pathogenicity, their further study should help to elucidate host-virus interactions that determine acute and chronic CHIKV infection.IMPORTANCE CHIKV is a globally spreading, mosquito-transmitted virus that causes debilitating acute and chronic musculoskeletal disease in humans. The viral genetic determinants that dictate the severity of acute and chronic diseases are not understood. To improve our understanding of CHIKV pathogenesis, we evaluated a CHIKV strain isolated from the serum of chronically infected immunocompromised mice. Sequence analysis of this persistent CHIKV strain identified two mutations, an amino acid change in the E2 viral attachment protein and a deletion within the 3'-UTR of the viral genome. We identified roles for these mutations in the enhancement of viral dissemination from the inoculation site and in disease severity. These data improve our understanding of the viral determinants of CHIKV pathogenesis and adaptive changes that occur during viral persistence.

Multiple conformations are a conserved and regulatory feature of the RB1 5' UTR.

Folding to a well-defined conformation is essential for the function of structured ribonucleic acids (RNAs) like the ribosome and tRNA. Structured elements in the untranslated regions (UTRs) of specific messenger RNAs (mRNAs) are known to control expression. The importance of unstructured regions adopting multiple conformations, however, is still poorly understood. High-resolution SHAPE-directed Boltzmann suboptimal sampling of the Homo sapiens Retinoblastoma 1 (RB1) 5' UTR yields three distinct conformations compatible with the experimental data. Private single nucleotide variants (SNVs) identified in two patients with retinoblastoma each collapse the structural ensemble to a single but distinct well-defined conformation. The RB1 5' UTRs from Bos taurus (cow) and Trichechus manatus latirostris (manatee) are divergent in sequence from H. sapiens (human) yet maintain structural compatibility with high-probability base pairs. SHAPE chemical probing of the cow and manatee RB1 5' UTRs reveals that they also adopt multiple conformations. Luciferase reporter assays reveal that 5' UTR mutations alter RB1 expression. In a traditional model of disease, causative SNVs disrupt a key structural element in the RNA. For the subset of patients with heritable retinoblastoma-associated SNVs in the RB1 5' UTR, the absence of multiple structures is likely causative of the cancer. Our data therefore suggest that selective pressure will favor multiple conformations in eukaryotic UTRs to regulate expression.

Primary and Secondary Sequence Structure Requirements for Recognition and Discrimination of Target RNAs by Pseudomonas aeruginosa RsmA and RsmF.

UNLABELLED: CsrA family RNA-binding proteins are widely distributed in bacteria and regulate gene expression at the posttranscriptional level. Pseudomonas aeruginosa has a canonical member of the CsrA family (RsmA) and a novel, structurally distinct variant (RsmF). To better understand RsmF binding properties, we performed parallel systematic evolution of ligands by exponential enrichment (SELEX) experiments for RsmA and RsmF. The initial target library consisted of 62-nucleotide (nt) RNA transcripts with central cores randomized at 15 sequential positions. Most targets selected by RsmA and RsmF were the expected size and shared a common consensus sequence (CANGGAYG) that was positioned in a hexaloop region of the stem-loop structure. RsmA and RsmF also selected for longer targets (≥96 nt) that were likely generated by rare PCR errors. Most of the long targets contained two consensus-binding sites. Representative short (single consensus site) and long (two consensus sites) targets were tested for RsmA and RsmF binding. Whereas RsmA bound the short targets with high affinity, RsmF was unable to bind the same targets. RsmA and RsmF both bound the long targets. Mutation of either consensus GGA site in the long targets reduced or eliminated RsmF binding, suggesting a requirement for two tandem binding sites. Conversely, RsmA bound long targets containing only a single GGA site with unaltered affinity. The RsmF requirement for two binding sites was confirmed with tssA1, an in vivo regulatory target of RsmA and RsmF. Our findings suggest that RsmF binding requires two GGA-containing sites, while RsmA binding requirements are less stringent. IMPORTANCE: The CsrA family of RNA-binding proteins is widely conserved in bacteria and plays important roles in the posttranscriptional regulation of protein synthesis. P. aeruginosa has two CsrA proteins, RsmA and RsmF. Although RsmA and RsmF share a few RNA targets, RsmF is unable to bind to other targets recognized by RsmA. The goal of the present study was to better understand the basis for differential binding by RsmF. Our data indicate that RsmF binding requires target RNAs with two consensus-binding sites, while RsmA recognizes targets with just a single binding site. This information should prove useful to future efforts to define the RsmF regulon and its contribution to P. aeruginosa physiology and virulence.

Specific RNA structures in the 5' untranslated region of the human cytomegalovirus major immediate early transcript are critical for efficient virus replication.

Human cytomegalovirus (HCMV) requires the robust expression of two immediate early proteins, IE1 and IE2, immediately upon infection to suppress the antiviral response and promote viral gene expression. While transcriptional control of IE1 and IE2 has been extensively studied, the role of post-transcriptional regulation of IE1 and IE2 expression is relatively unexplored. We previously found that the shared major immediate early 5' untranslated region (MIE 5' UTR) of the mature IE1 and IE2 transcripts plays a critical role in facilitating the translation of the IE1 and IE2 mRNAs. As RNA secondary structure in 5' UTRs can regulate mRNA translation efficiency, we used selective 2'-hydroxyl acylation analyzed by primer extension and mutational profiling (SHAPE-MaP) to identify RNA structures in the shared MIE 5' UTR. We found that the MIE 5' UTR contains three stable stem loop structures. Using a series of recombinant viruses to investigate the role of each stem loop in IE1 and IE2 protein synthesis, we found that the stem loop closest to the 5' end of the MIE 5' UTR (SL1) is both necessary and sufficient for efficient IE1 and IE2 mRNA translation and HCMV replication. The positive effect of SL1 on mRNA translation and virus replication was dependent on its location within the 5' UTR. Surprisingly, a synthetic stem loop with the same free energy as SL1 in its native location also supported wild type levels of IE1 and IE2 mRNA translation and virus replication, suggesting that the presence of RNA structure at a specific location in the 5' UTR, rather than the primary sequence of the RNA, is critical for efficient IE1 and IE2 protein synthesis. These data reveal a novel post-transcriptional regulatory mechanism controlling IE1 and IE2 expression and reinforce the critical role of RNA structure in regulating HCMV protein synthesis and replication.IMPORTANCEThese results reveal a new aspect of immediate early gene regulation controlled by non-coding RNA structures in viral mRNAs. Previous studies have largely focused on understanding viral gene expression at the level of transcriptional control. Our results show that a complete understanding of the control of viral gene expression must include an understanding of viral mRNA translation, which is driven in part by RNA structure(s) in the 5' UTR of viral mRNAs. Our results illustrate the importance of these additional layers of regulation by defining specific 5' UTR RNA structures regulating immediate early gene expression in the context of infection and identify important features of RNA structure that govern viral mRNA translation efficiency. These results may therefore broadly impact current thinking on how viral gene expression is regulated for human cytomegalovirus and other DNA viruses.

The Human Cytomegalovirus vGPCR UL33 is Essential for Efficient Lytic Replication in Epithelial Cells.

Human cytomegalovirus (HCMV) is a ß-herpesvirus which is ubiquitous in the human population. HCMV has the largest genome of all known human herpesviruses, and thus encodes a large array of proteins that affect pathogenesis in different cell types. Given the large genome and the ability of HCMV to replicate in a range of cells, investigators have begun to identify viral proteins required for cell type-specific replication. There are four proteins encoded in the HCMV genome that are homologous to human G protein-coupled receptors (GPCRs); these viral-encoded GPCRs (vGPCRs) are UL33, UL78, US27, and US28. In the current study, we find that deletion of all four vGPCR genes from a clinical isolate of HCMV severely attenuates lytic replication in both primary human salivary gland epithelial cells, as well as ARPE-19 retinal epithelial cells as evidenced by significant decreases in immediate early gene expression and virus production. Deletion of UL33 from the HCMV genome also results in a failure to efficiently replicate in epithelial cells, and this defect is manifested by decreased levels of immediate early, early, and late gene expression, as well as reduced viral production. We find that similar to US28, UL33 constitutively activates Gαq-dependent PLC-ß signaling to high levels in these epithelial cells. We also find that UL33 transcription is more complicated than originally believed, and there is the potential for the virus to utilize various 5' UTRs to create novel UL33 proteins that are all capable of constitutive Gαq signaling. Taken together, these studies suggest that UL33 driven signaling is important for lytic HCMV replication in cells of epithelial origin.

Chikungunya Virus RNA Secondary Structures Impact Defective Viral Genome Production.

Chikungunya virus (CHIKV) is a mosquito-borne RNA virus that poses an emerging threat to humans. In a manner similar to other RNA viruses, CHIKV encodes an error-prone RNA polymerase which, in addition to producing full-length genomes, gives rise to truncated, non-functional genomes, which have been coined defective viral genomes (DVGs). DVGs have been intensively studied in the context of therapy, as they can inhibit viral replication and dissemination in their hosts. In this work, we interrogate the influence of viral RNA secondary structures on the production of CHIKV DVGs. We experimentally map RNA secondary structures of the CHIKV genome using selective 2'-hydroxyl acylation analyzed by primer extension and mutational profiling (SHAPE-MaP), which couples chemical labelling with next-generation sequencing. We correlate the inferred secondary structure with preferred deletion sites of CHIKV DVGs. We document an increased probability of DVG generation with truncations at unpaired nucleotides within the secondary structure. We then generated a CHIKV mutant bearing synonymous changes at the nucleotide level to disrupt the existing RNA secondary structure (CHIKV-D2S). We show that CHIKV-D2S presents altered DVG generation compared to wild-type virus, correlating with the change in RNA secondary structure obtained by SHAPE-MaP. Our work thus demonstrates that RNA secondary structure impacts CHIKV DVG production during replication.

U-CAN-seq: A Universal Competition Assay by Nanopore Sequencing.

RNA viruses quickly evolve subtle genotypic changes that can have major impacts on viral fitness and host range, with potential consequences for human health. It is therefore important to understand the evolutionary fitness of novel viral variants relative to well-studied genotypes of epidemic viruses. Competition assays are an effective and rigorous system with which to assess the relative fitness of viral genotypes. However, it is challenging to quickly and cheaply distinguish and quantify fitness differences between very similar viral genotypes. Here, we describe a protocol for using reverse transcription PCR in combination with commercial nanopore sequencing services to perform competition assays on untagged RNA viruses. Our assay, called the Universal Competition Assay by Nanopore Sequencing (U-CAN-seq), is relatively cheap and highly sensitive. We used a well-studied N24A mutation in the chikungunya virus (CHIKV) nsp3 gene to confirm that we could detect a competitive disadvantage using U-CAN-seq. We also used this approach to show that mutations to the CHIKV 5' conserved sequence element that disrupt sequence but not structure did not affect the fitness of CHIKV. However, similar mutations to an adjacent CHIKV stem loop (SL3) did cause a fitness disadvantage compared to wild-type CHIKV, suggesting that structure-independent, primary sequence determinants in this loop play an important role in CHIKV biology. Our novel findings illustrate the utility of the U-CAN-seq competition assay.

A 44-Nucleotide Region in the Chikungunya Virus 3' UTR Dictates Viral Fitness in Disparate Host Cells.

We previously reported that deletion of a 44-nucleotide element in the 3' untranslated region (UTR) of the Chikungunya virus (CHIKV) genome enhances the virulence of CHIKV infection in mice. Here, we find that while this 44-nucleotide deletion enhances CHIKV fitness in murine embryonic fibroblasts in a manner independent of the type I interferon response, the same mutation decreases viral fitness in C6/36 mosquito cells. Further, the fitness advantage conferred by the UTR deletion in mammalian cells is maintained in vivo in a mouse model of CHIKV dissemination. Finally, SHAPE-MaP analysis of the CHIKV 3' UTR revealed this 44-nucleotide element forms a distinctive two-stem-loop structure that is ablated in the mutant 3' UTR without altering additional 3' UTR RNA secondary structures.

Metatranscriptomics analysis reveals a novel transcriptional and translational landscape during Middle East respiratory syndrome coronavirus infection.

Among all RNA viruses, coronavirus RNA transcription is the most complex and involves a process termed "discontinuous transcription" that results in the production of a set of 3'-nested, co-terminal genomic and subgenomic RNAs during infection. While the expression of the classic canonical set of subgenomic RNAs depends on the recognition of a 6- to 7-nt transcription regulatory core sequence (TRS), here, we use deep sequence and metagenomics analysis strategies and show that the coronavirus transcriptome is even more vast and more complex than previously appreciated and involves the production of leader-containing transcripts that have canonical and noncanonical leader-body junctions. Moreover, by ribosome protection and proteomics analyses, we show that both positive- and negative-sense transcripts are translationally active. The data support the hypothesis that the coronavirus proteome is much vaster than previously noted in the literature.

Core mitochondrial genes are down-regulated during SARS-CoV-2 infection of rodent and human hosts.

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) viral proteins bind to host mitochondrial proteins, likely inhibiting oxidative phosphorylation (OXPHOS) and stimulating glycolysis. We analyzed mitochondrial gene expression in nasopharyngeal and autopsy tissues from patients with coronavirus disease 2019 (COVID-19). In nasopharyngeal samples with declining viral titers, the virus blocked the transcription of a subset of nuclear DNA (nDNA)-encoded mitochondrial OXPHOS genes, induced the expression of microRNA 2392, activated HIF-1α to induce glycolysis, and activated host immune defenses including the integrated stress response. In autopsy tissues from patients with COVID-19, SARS-CoV-2 was no longer present, and mitochondrial gene transcription had recovered in the lungs. However, nDNA mitochondrial gene expression remained suppressed in autopsy tissue from the heart and, to a lesser extent, kidney, and liver, whereas mitochondrial DNA transcription was induced and host-immune defense pathways were activated. During early SARS-CoV-2 infection of hamsters with peak lung viral load, mitochondrial gene expression in the lung was minimally perturbed but was down-regulated in the cerebellum and up-regulated in the striatum even though no SARS-CoV-2 was detected in the brain. During the mid-phase SARS-CoV-2 infection of mice, mitochondrial gene expression was starting to recover in mouse lungs. These data suggest that when the viral titer first peaks, there is a systemic host response followed by viral suppression of mitochondrial gene transcription and induction of glycolysis leading to the deployment of antiviral immune defenses. Even when the virus was cleared and lung mitochondrial function had recovered, mitochondrial function in the heart, kidney, liver, and lymph nodes remained impaired, potentially leading to severe COVID-19 pathology.

TARGETED DOWN REGULATION OF CORE MITOCHONDRIAL GENES DURING SARS-COV-2 INFECTION.

Defects in mitochondrial oxidative phosphorylation (OXPHOS) have been reported in COVID-19 patients, but the timing and organs affected vary among reports. Here, we reveal the dynamics of COVID-19 through transcription profiles in nasopharyngeal and autopsy samples from patients and infected rodent models. While mitochondrial bioenergetics is repressed in the viral nasopharyngeal portal of entry, it is up regulated in autopsy lung tissues from deceased patients. In most disease stages and organs, discrete OXPHOS functions are blocked by the virus, and this is countered by the host broadly up regulating unblocked OXPHOS functions. No such rebound is seen in autopsy heart, results in severe repression of genes across all OXPHOS modules. Hence, targeted enhancement of mitochondrial gene expression may mitigate the pathogenesis of COVID-19.

Comparative analysis of coronavirus genomic RNA structure reveals conservation in SARS-like coronaviruses.

Coronaviruses, including SARS-CoV-2 the etiological agent of COVID-19 disease, have caused multiple epidemic and pandemic outbreaks in the past 20 years1-3. With no vaccines, and only recently developed antiviral therapeutics, we are ill equipped to handle coronavirus outbreaks4. A better understanding of the molecular mechanisms that regulate coronavirus replication and pathogenesis is needed to guide the development of new antiviral therapeutics and vaccines. RNA secondary structures play critical roles in multiple aspects of coronavirus replication, but the extent and conservation of RNA secondary structure across coronavirus genomes is unknown5. Here, we define highly structured RNA regions throughout the MERS-CoV, SARS-CoV, and SARS-CoV-2 genomes. We find that highly stable RNA structures are pervasive throughout coronavirus genomes, and are conserved between the SARS-like CoV. Our data suggests that selective pressure helps preserve RNA secondary structure in coronavirus genomes, suggesting that these structures may play important roles in virus replication and pathogenesis. Thus, disruption of conserved RNA secondary structures could be a novel strategy for the generation of attenuated SARS-CoV-2 vaccines for use against the current COVID-19 pandemic.