Bipartite viral RNA genome heterodimerization influences genome packaging and virion thermostability.

Multi-segmented viruses often multimerize their genomic segments to ensure efficient and stoichiometric packaging of the correct genetic cargo. In the bipartite Nodaviridae family, genome heterodimerization is also observed and conserved among different species. However, the nucleotide composition and biological function for this heterodimer remain unclear. Using Flock House virus as a model system, we developed a next-generation sequencing approach ("XL-ClickSeq") to probe heterodimer site sequences. We identified an intermolecular base-pairing site which contributed to heterodimerization in both wild-type and defective virus particles. Mutagenic disruption of this heterodimer site exhibited significant deficiencies in genome packaging and encapsidation specificity to viral genomic RNAs. Furthermore, the disruption of this intermolecular interaction directly impacts the thermostability of the mature virions. These results demonstrate that the intermolecular RNA-RNA interactions within the encapsidated genome of an RNA virus have an important role on virus particle integrity and thus may impact its transmission to a new host.IMPORTANCEFlock House virus is a member of Nodaviridae family of viruses, which provides a well-studied model virus for non-enveloped RNA virus assembly, cell entry, and replication. The Flock House virus genome consists of two separate RNA molecules, which can form a heterodimer upon heating of virus particles. Although similar RNA dimerization is utilized by other viruses (such as retroviruses) as a packaging mechanism and is conserved among Nodaviruses, the role of heterodimerization in the Nodavirus replication cycle is unclear. In this research, we identified the RNA sequences contributing to Flock House virus genome heterodimerization and discovered that such RNA-RNA interaction plays an essential role in virus packaging efficiency and particle integrity. This provides significant insight into how the interaction of packaged viral RNA may have a broader impact on the structural and functional properties of virus particles.

RpoS activates Salmonella Typhi biofilms and drives persistence in a Zebrafish model.

The development of strategies for targeting the asymptomatic carriage of Salmonella Typhi in chronic typhoid patients has suffered owing to our basic lack of understanding of the molecular mechanisms that enable the formation of S. Typhi biofilms. Traditionally, studies have relied on cholesterol-attached biofilms formed by a closely related serovar, Typhimurium, to mimic multicellular Typhi communities formed on human gallstones. In long-term infections, S. Typhi adopts the biofilm lifestyle to persist in vivo and survive in the carrier state, ultimately leading to the spread of infections via the fecal-oral route of transmission. In the present work, we studied S. Typhi biofilms directly, applied targeted as well as genome-wide genetic approaches to uncover unique biofilm components that do not conform to the CsgD-dependent pathway as established in S. Typhimurium. We adopted a genome-wide Tn5 mutation screen in S. Typhi in gallstone-mimicking conditions and generated New Generation Sequencing libraries based on the ClickSeq technology to identify the key regulators, IraP and RpoS, and the matrix components as Sth fimbriae, Vi capsule and lipopolysaccharide. We discovered that the starvation sigma factor, RpoS, was required for the transcriptional activation of matrix-encoding genes in vitro, and for S. Typhi colonization in persistent infections in vivo, using a heterologous fish larval model. Overall, our work established a novel RpoS-driven paradigm for the formation of cholesterol-attached Typhi biofilms and emphasized the role(s) of stress signaling pathways for adaptation in chronic infections.

Loss-of-function mutation in Omicron variants reduces spike protein expression and attenuates SARS-CoV-2 infection.

SARS-CoV-2 Omicron variants emerged in 2022 with >30 novel amino acid mutations in the spike protein alone. While most studies focus on receptor binding domain changes, mutations in the C-terminus of S1 (CTS1), adjacent to the furin cleavage site, have largely been ignored. In this study, we examined three Omicron mutations in CTS1: H655Y, N679K, and P681H. Generating a SARS-CoV-2 triple mutant (YKH), we found that the mutant increased spike processing, consistent with prior reports for H655Y and P681H individually. Next, we generated a single N679K mutant, finding reduced viral replication in vitro and less disease in vivo. Mechanistically, the N679K mutant had reduced spike protein in purified virions compared to wild-type; spike protein decreases were further exacerbated in infected cell lysates. Importantly, exogenous spike expression also revealed that N679K reduced overall spike protein yield independent of infection. Although a loss-of-function mutation, transmission competition demonstrated that N679K had a replication advantage in the upper airway over wild-type SARS-CoV-2 in hamsters, potentially impacting transmissibility. Together, the data show that N679K reduces overall spike protein levels during Omicron infection, which has important implications for infection, immunity, and transmission.

ViReMa: a virus recombination mapper of next-generation sequencing data characterizes diverse recombinant viral nucleic acids.

BACKGROUND: Genetic recombination is a tremendous source of intrahost diversity in viruses and is critical for their ability to rapidly adapt to new environments or fitness challenges. While viruses are routinely characterized using high-throughput sequencing techniques, characterizing the genetic products of recombination in next-generation sequencing data remains a challenge. Viral recombination events can be highly diverse and variable in nature, including simple duplications and deletions, or more complex events such as copy/snap-back recombination, intervirus or intersegment recombination, and insertions of host nucleic acids. Due to the variable mechanisms driving virus recombination and the different selection pressures acting on the progeny, recombination junctions rarely adhere to simple canonical sites or sequences. Furthermore, numerous different events may be present simultaneously in a viral population, yielding a complex mutational landscape. FINDINGS: We have previously developed an algorithm called ViReMa (Virus Recombination Mapper) that bootstraps the bowtie short-read aligner to capture and annotate a wide range of recombinant species found within virus populations. Here, we have updated ViReMa to provide an "error density" function designed to accurately detect recombination events in the longer reads now routinely generated by the Illumina platforms and provide output reports for multiple types of recombinant species using standardized formats. We demonstrate the utility and flexibility of ViReMa in different settings to report deletion events in simulated data from Flock House virus, copy-back RNA species in Sendai viruses, short duplication events in HIV, and virus-to-host recombination in an archaeal DNA virus.

SARS-CoV-2 Uses Nonstructural Protein 16 To Evade Restriction by IFIT1 and IFIT3.

Understanding the molecular basis of innate immune evasion by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) is an important consideration for designing the next wave of therapeutics. Here, we investigate the role of the nonstructural protein 16 (NSP16) of SARS-CoV-2 in infection and pathogenesis. NSP16, a ribonucleoside 2'-O-methyltransferase (MTase), catalyzes the transfer of a methyl group to mRNA as part of the capping process. Based on observations with other CoVs, we hypothesized that NSP16 2'-O-MTase function protects SARS-CoV-2 from cap-sensing host restriction. Therefore, we engineered SARS-CoV-2 with a mutation that disrupts a conserved residue in the active site of NSP16. We subsequently show that this mutant is attenuated both in vitro and in vivo, using a hamster model of SARS-CoV-2 infection. Mechanistically, we confirm that the NSP16 mutant is more sensitive than wild-type SARS-CoV-2 to type I interferon (IFN-I) in vitro. Furthermore, silencing IFIT1 or IFIT3, IFN-stimulated genes that sense a lack of 2'-O-methylation, partially restores fitness to the NSP16 mutant. Finally, we demonstrate that sinefungin, an MTase inhibitor that binds the catalytic site of NSP16, sensitizes wild-type SARS-CoV-2 to IFN-I treatment and attenuates viral replication. Overall, our findings highlight the importance of SARS-CoV-2 NSP16 in evading host innate immunity and suggest a target for future antiviral therapies. IMPORTANCE Similar to other coronaviruses, disruption of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) NSP16 function attenuates viral replication in a type I interferon-dependent manner. In vivo, our results show reduced disease and viral replication at late times in the hamster lung, but an earlier titer deficit for the NSP16 mutant (dNSP16) in the upper airway. In addition, our results confirm a role for IFIT1 but also demonstrate the necessity of IFIT3 in mediating dNSP16 attenuation. Finally, we show that targeting NSP16 activity with a 2'-O-methyltransferase inhibitor in combination with type I interferon offers a novel avenue for antiviral development.

SARS-CoV-2 Uses Nonstructural Protein 16 to Evade Restriction by IFIT1 and IFIT3.

Understanding the molecular basis of innate immune evasion by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) is an important consideration for designing the next wave of therapeutics. Here, we investigate the role of the nonstructural protein 16 (NSP16) of SARS-CoV-2 in infection and pathogenesis. NSP16, a ribonucleoside 2'- O methyltransferase (MTase), catalyzes the transfer of a methyl group to mRNA as part of the capping process. Based on observations with other CoVs, we hypothesized that NSP16 2'- O MTase function protects SARS-CoV-2 from cap-sensing host restriction. Therefore, we engineered SARS-CoV-2 with a mutation that disrupts a conserved residue in the active site of NSP16. We subsequently show that this mutant is attenuated both in vitro and in vivo , using a hamster model of SARS-CoV-2 infection. Mechanistically, we confirm that the NSP16 mutant is more sensitive to type I interferon (IFN-I) in vitro . Furthermore, silencing IFIT1 or IFIT3, IFN-stimulated genes that sense a lack of 2'- O methylation, partially restores fitness to the NSP16 mutant. Finally, we demonstrate that sinefungin, a methyltransferase inhibitor that binds the catalytic site of NSP16, sensitizes wild-type SARS-CoV-2 to IFN-I treatment. Overall, our findings highlight the importance of SARS-CoV-2 NSP16 in evading host innate immunity and suggest a possible target for future antiviral therapies. Importance: Similar to other coronaviruses, disruption of SARS-CoV-2 NSP16 function attenuates viral replication in a type I interferon-dependent manner. In vivo , our results show reduced disease and viral replication at late times in the hamster lung, but an earlier titer deficit for the NSP16 mutant (dNSP16) in the upper airway. In addition, our results confirm a role for IFIT1, but also demonstrate the necessity of IFIT3 in mediating dNSP16 attenuation. Finally, we show that targeting NSP16 activity with a 2'- O methyltransferase inhibitor in combination with type I interferon offers a novel avenue for antiviral development.

Zika Virus Infection Alters Gene Expression and Poly-Adenylation Patterns in Placental Cells.

Flaviviruses are small RNA viruses that are mainly transmitted via arthropod vectors and are found in tropic and sub-tropical regions. Most infections are asymptomatic (90-95%), but symptoms can be as severe as hemorrhagic fever and encephalitis. One recently emerged flavivirus is Zika virus (ZIKV), which was originally isolated from rhesus monkeys in Uganda roughly 70 years ago but has recently spread east, reaching S. America in 2015-2016. This outbreak was associated with the development of Guillain-Barré syndrome in adults and microcephaly in infants born to expectant mothers infected early in pregnancy. ZIKV must traverse the placenta to impact the development of the fetus, but the mechanisms responsible are unknown. While flaviviruses are known to disrupt splicing patterns in host cells, little is known about how flaviviruses such as ZIKV impact the alternative polyadenylation (APA) of host transcripts. This is important as APA is well-established as a mechanism in the regulation of mRNA metabolism and translation. Thus, we sought to characterize transcriptomic changes including APA in human placental (JEG3) cells in response to ZIKV infection using Poly(A)-ClickSeq (PAC-Seq). We used our differential Poly(A)-cluster (DPAC) analysis pipeline to characterize changes in differential gene expression, alternative poly-adenylation (APA) and the use of alternative terminal exons. We identified 98 upregulated genes and 28 downregulated genes. Pathway enrichment analysis indicated that many RNA processing and immune pathways were upregulated in ZIKV-infected JEG3 cells. We also updated DPAC to provide additional metrics of APA including the percentage-distal usage index (PDUI), which revealed that APA was extensive and the 3' UTRs of 229 genes were lengthened while 269 were shortened. We further found that there were 214 upregulated and 59 downregulated poly(A)-clusters (PACs). We extracted the nucleotide sequences surrounding these PACs and found that the canonical signals for poly-adenylation (binding site for poly-A binding protein (PABP) upstream and a GU-rich region down-stream of the PAC) were only enriched in the downregulated PACs. These results indicate that ZIKV infection makes JEG3 cells more permissive to non-canonical poly-adenylation signals.

QTQTN motif upstream of the furin-cleavage site plays a key role in SARS-CoV-2 infection and pathogenesis.

The furin cleavage site (FCS), an unusual feature in the SARS-CoV-2 spike protein, has been spotlighted as a factor key to facilitating infection and pathogenesis by increasing spike processing. Similarly, the QTQTN motif directly upstream of the FCS is also an unusual feature for group 2B coronaviruses (CoVs). The QTQTN deletion has consistently been observed in in vitro cultured virus stocks and some clinical isolates. To determine whether the QTQTN motif is critical to SARS-CoV-2 replication and pathogenesis, we generated a mutant deleting the QTQTN motif (ΔQTQTN). Here, we report that the QTQTN deletion attenuates viral replication in respiratory cells in vitro and attenuates disease in vivo. The deletion results in a shortened, more rigid peptide loop that contains the FCS and is less accessible to host proteases, such as TMPRSS2. Thus, the deletion reduced the efficiency of spike processing and attenuates SARS-CoV-2 infection. Importantly, the QTQTN motif also contains residues that are glycosylated, and disruption of its glycosylation also attenuates virus replication in a TMPRSS2-dependent manner. Together, our results reveal that three aspects of the S1/S2 cleavage site-the FCS, loop length, and glycosylation-are required for efficient SARS-CoV-2 replication and pathogenesis.

ViReMaShiny: an interactive application for analysis of viral recombination data.

MOTIVATION: Recombination is an essential driver of virus evolution and adaption, giving rise to new chimeric viruses, structural variants, sub-genomic RNAs and defective RNAs. Next-generation sequencing (NGS) of virus samples, either from experimental or clinical settings, has revealed a complex distribution of recombination events that contributes to intrahost diversity. We and others have previously developed alignment tools to discover and map these diverse recombination events in NGS data. However, there is no standard for data visualization to contextualize events of interest, and downstream analysis often requires bespoke coding. RESULTS: We present ViReMaShiny, a web-based application built using the R Shiny framework to allow interactive exploration and point-and-click visualization of viral recombination data provided in BED format generated by computational pipelines such as ViReMa (Viral-Recombination-Mapper). AVAILABILITY AND IMPLEMENTATION: The application is hosted at https://routhlab.shinyapps.io/ViReMaShiny/ with associated documentation at https://jayeung12.github.io/. Code is available at https://github.com/routhlab/ViReMaShiny. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Selective ablation of 3' RNA ends and processive RTs facilitate direct cDNA sequencing of full-length host cell and viral transcripts.

Alternative splicing (AS) is necessary for viral proliferation in host cells and a critical regulatory component of viral gene expression. Conventional RNA-seq approaches provide incomplete coverage of AS due to their short read lengths and are susceptible to biases and artifacts introduced in prevailing library preparation methodologies. Moreover, viral splicing studies are often conducted separately from host cell transcriptome analysis, precluding an assessment of the viral manipulation of host splicing machinery. To address current limitations, we developed a quantitative full-length direct cDNA sequencing strategy to simultaneously profile viral and host cell transcripts. This nanopore-based approach couples processive reverse transcriptases with a novel one-step chemical ablation of 3' RNA ends (termed CASPR), which decreases ribosomal RNA reads and enriches polyadenylated coding sequences. We extensively validate our approach using synthetic reference transcripts and show that CASPR doubles the breadth of coverage per transcript and increases detection of long transcripts (>4 kb), while being functionally equivalent to PolyA+ selection for transcript quantification. We used our approach to interrogate host cell and HIV-1 transcript dynamics during viral reactivation and identified novel putative HIV-1 host factors containing exon skipping or novel intron retentions and delineated the HIV-1 transcriptional state associated with these differentially regulated host factors.

Nucleocapsid mutations in SARS-CoV-2 augment replication and pathogenesis.

While SARS-CoV-2 continues to adapt for human infection and transmission, genetic variation outside of the spike gene remains largely unexplored. This study investigates a highly variable region at residues 203-205 in the SARS-CoV-2 nucleocapsid protein. Recreating a mutation found in the alpha and omicron variants in an early pandemic (WA-1) background, we find that the R203K+G204R mutation is sufficient to enhance replication, fitness, and pathogenesis of SARS-CoV-2. The R203K+G204R mutant corresponds with increased viral RNA and protein both in vitro and in vivo. Importantly, the R203K+G204R mutation increases nucleocapsid phosphorylation and confers resistance to inhibition of the GSK-3 kinase, providing a molecular basis for increased virus replication. Notably, analogous alanine substitutions at positions 203+204 also increase SARS-CoV-2 replication and augment phosphorylation, suggesting that infection is enhanced through ablation of the ancestral 'RG' motif. Overall, these results demonstrate that variant mutations outside spike are key components in SARS-CoV-2's continued adaptation to human infection.

Next-generation sequencing: A new avenue to understand viral RNA-protein interactions.

The genomes of RNA viruses present an astonishing source of both sequence and structural diversity. From intracellular viral RNA-host interfaces to interactions between the RNA genome and structural proteins in virus particles themselves, almost the entire viral lifecycle is accompanied by a myriad of RNA-protein interactions that are required to fulfill their replicative potential. It is therefore important to characterize such rich and dynamic collections of viral RNA-protein interactions to understand virus evolution and their adaptation to their hosts and environment. Recent advances in next-generation sequencing technologies have allowed the characterization of viral RNA-protein interactions, including both transient and conserved interactions, where molecular and structural approaches have fallen short. In this review, we will provide a methodological overview of the high-throughput techniques used to study viral RNA-protein interactions, their biochemical mechanisms, and how they evolved from classical methods as well as one another. We will discuss how different techniques have fueled virus research to characterize how viral RNA and proteins interact, both locally and on a global scale. Finally, we will present examples on how these techniques influence the studies of clinically important pathogens such as HIV-1 and SARS-CoV-2.

Covariation of viral recombination with single nucleotide variants during virus evolution revealed by CoVaMa.

Adaptation of viruses to their environments occurs through the acquisition of both novel single-nucleotide variants (SNV) and recombination events including insertions, deletions, and duplications. The co-occurrence of SNVs in individual viral genomes during their evolution has been well-described. However, unlike covariation of SNVs, studying the correlation between recombination events with each other or with SNVs has been hampered by their inherent genetic complexity and a lack of bioinformatic tools. Here, we expanded our previously reported CoVaMa pipeline (v0.1) to measure linkage disequilibrium between recombination events and SNVs within both short-read and long-read sequencing datasets. We demonstrate this approach using long-read nanopore sequencing data acquired from Flock House virus (FHV) serially passaged in vitro. We found SNVs that were either correlated or anti-correlated with large genomic deletions generated by nonhomologous recombination that give rise to Defective-RNAs. We also analyzed NGS data from longitudinal HIV samples derived from a patient undergoing antiretroviral therapy who proceeded to virological failure. We found correlations between insertions in the p6Gag and mutations in Gag cleavage sites. This report confirms previous findings and provides insights on novel associations between SNVs and specific recombination events within the viral genome and their role in viral evolution.

Nucleocapsid mutations in SARS-CoV-2 augment replication and pathogenesis.

While SARS-CoV-2 continues to adapt for human infection and transmission, genetic variation outside of the spike gene remains largely unexplored. This study investigates a highly variable region at residues 203-205 in the SARS-CoV-2 nucleocapsid protein. Recreating a mutation found in the alpha and omicron variants in an early pandemic (WA-1) background, we find that the R203K+G204R mutation is sufficient to enhance replication, fitness, and pathogenesis of SARS-CoV-2. The R203K+G204R mutant corresponds with increased viral RNA and protein both in vitro and in vivo . Importantly, the R203K+G204R mutation increases nucleocapsid phosphorylation and confers resistance to inhibition of the GSK-3 kinase, providing a molecular basis for increased virus replication. Notably, analogous alanine substitutions at positions 203+204 also increase SARS-CoV-2 replication and augment phosphorylation, suggesting that infection is enhanced through ablation of the ancestral 'RG' motif. Overall, these results demonstrate that variant mutations outside spike are key components in SARS-CoV-2's continued adaptation to human infection. AUTHOR SUMMARY: Since its emergence, SARS-CoV-2 has continued to adapt for human infection resulting in the emergence of variants with unique genetic profiles. Most studies of genetic variation have focused on spike, the target of currently available vaccines, leaving the importance of variation elsewhere understudied. Here, we characterize a highly variable motif at residues 203-205 in nucleocapsid. Recreating the prominent nucleocapsid R203K+G204R mutation in an early pandemic background, we show that this mutation is alone sufficient to enhance SARS-CoV-2 replication and pathogenesis. We also link augmentation of SARS-CoV-2 infection by the R203K+G204R mutation to its modulation of nucleocapsid phosphorylation. Finally, we characterize an analogous alanine double substitution at positions 203-204. This mutant was found to mimic R203K+G204R, suggesting augmentation of infection occurs by disrupting the ancestral sequence. Together, our findings illustrate that mutations outside of spike are key components of SARS-CoV-2's adaptation to human infection.

Mouse-adapted SARS-CoV-2 protects animals from lethal SARS-CoV challenge.

The emergence of Severe Acute Respiratory Syndrome Coronavirus 2 (SARS-CoV-2) has resulted in a pandemic causing significant damage to public health and the economy. Efforts to understand the mechanisms of Coronavirus Disease 2019 (COVID-19) have been hampered by the lack of robust mouse models. To overcome this barrier, we used a reverse genetic system to generate a mouse-adapted strain of SARS-CoV-2. Incorporating key mutations found in SARS-CoV-2 variants, this model recapitulates critical elements of human infection including viral replication in the lung, immune cell infiltration, and significant in vivo disease. Importantly, mouse adaptation of SARS-CoV-2 does not impair replication in human airway cells and maintains antigenicity similar to human SARS-CoV-2 strains. Coupled with the incorporation of mutations found in variants of concern, CMA3p20 offers several advantages over other mouse-adapted SARS-CoV-2 strains. Using this model, we demonstrate that SARS-CoV-2-infected mice are protected from lethal challenge with the original Severe Acute Respiratory Syndrome Coronavirus (SARS-CoV), suggesting immunity from heterologous Coronavirus (CoV) strains. Together, the results highlight the use of this mouse model for further study of SARS-CoV-2 infection and disease.

RBFOX2 is critical for maintaining alternative polyadenylation patterns and mitochondrial health in rat myoblasts.

RBFOX2, which has a well-established role in alternative splicing, is linked to heart diseases. However, it is unclear whether RBFOX2 has other roles in RNA processing that can influence gene expression in muscle cells, contributing to heart disease. Here, we employ both 3'-end and nanopore cDNA sequencing to reveal a previously unrecognized role for RBFOX2 in maintaining alternative polyadenylation (APA) signatures in myoblasts. RBFOX2-mediated APA modulates mRNA levels and/or isoform expression of a collection of genes, including contractile and mitochondrial genes. Depletion of RBFOX2 adversely affects mitochondrial health in myoblasts, correlating with disrupted APA of mitochondrial gene Slc25a4. Mechanistically, RBFOX2 regulation of Slc25a4 APA is mediated through consensus RBFOX2 binding motifs near the distal polyadenylation site, enforcing the use of the proximal polyadenylation site. In sum, our results unveil a role for RBFOX2 in fine-tuning expression of mitochondrial and contractile genes via APA in myoblasts relevant to heart diseases.

Chikungunya Virus Infects the Heart and Induces Heart-Specific Transcriptional Changes in an Immunodeficient Mouse Model of Infection.

Chikungunya virus (CHIKV) is a mosquito-transmitted pathogen in family Togaviridae, genus Alphavirus. Although CHIKV is well known for its ability to cause debilitating rheumatoid-like arthritis, it has been also been observed to cause cardiovascular symptoms such as arrhythmias. Here, using samples from a previous study, we sequenced RNA from serum, kidney, skeletal muscle, and cardiac muscle from CHIKV- and mock-infected IFN-αR-/- mice using two sequencing techniques to investigate heart-specific changes in virus mutational profiles and host gene expression. Mutation rates were similar across muscle tissues although heart tissue carried heart-specific CHIKV minority variants, one of which had a coding change in the nsP3 gene and another in the 3'UTR. Importantly, heart-specific transcriptional changes included differential expression of genes critical for ion transport and muscle contraction. These results demonstrate that CHIKV replicates in the hearts of immunodeficient mice and induce heart-specific mutations and host responses with implications for cardiac pathologies.

Tiled-ClickSeq for targeted sequencing of complete coronavirus genomes with simultaneous capture of RNA recombination and minority variants.

High-throughput genomics of SARS-CoV-2 is essential to characterize virus evolution and to identify adaptations that affect pathogenicity or transmission. While single-nucleotide variations (SNVs) are commonly considered as driving virus adaption, RNA recombination events that delete or insert nucleic acid sequences are also critical. Whole genome targeting sequencing of SARS-CoV-2 is typically achieved using pairs of primers to generate cDNA amplicons suitable for next-generation sequencing (NGS). However, paired-primer approaches impose constraints on where primers can be designed, how many amplicons are synthesized and requires multiple PCR reactions with non-overlapping primer pools. This imparts sensitivity to underlying SNVs and fails to resolve RNA recombination junctions that are not flanked by primer pairs. To address these limitations, we have designed an approach called 'Tiled-ClickSeq', which uses hundreds of tiled-primers spaced evenly along the virus genome in a single reverse-transcription reaction. The other end of the cDNA amplicon is generated by azido-nucleotides that stochastically terminate cDNA synthesis, removing the need for a paired-primer. A sequencing adaptor containing a Unique Molecular Identifier (UMI) is appended to the cDNA fragment using click-chemistry and a PCR reaction generates a final NGS library. Tiled-ClickSeq provides complete genome coverage, including the 5'UTR, at high depth and specificity to the virus on both Illumina and Nanopore NGS platforms. Here, we analyze multiple SARS-CoV-2 isolates and clinical samples to simultaneously characterize minority variants, sub-genomic mRNAs (sgmRNAs), structural variants (SVs) and D-RNAs. Tiled-ClickSeq therefore provides a convenient and robust platform for SARS-CoV-2 genomics that captures the full range of RNA species in a single, simple assay.

Nanopore sequencing reveals full-length Tropomyosin 1 isoforms and their regulation by RNA-binding proteins during rat heart development.

Alternative splicing (AS) contributes to the diversity of the proteome by producing multiple isoforms from a single gene. Although short-read RNA-sequencing methods have been the gold standard for determining AS patterns of genes, they have a difficulty in defining full-length mRNA isoforms assembled using different exon combinations. Tropomyosin 1 (TPM1) is an actin-binding protein required for cytoskeletal functions in non-muscle cells and for contraction in muscle cells. Tpm1 undergoes AS regulation to generate muscle versus non-muscle TPM1 protein isoforms with distinct physiological functions. It is unclear which full-length Tpm1 isoforms are produced via AS and how they are regulated during heart development. To address these, we utilized nanopore long-read cDNA sequencing without gene-specific PCR amplification. In rat hearts, we identified full-length Tpm1 isoforms composed of distinct exons with specific exon linkages. We showed that Tpm1 undergoes AS transitions during embryonic heart development such that muscle-specific exons are connected generating predominantly muscle-specific Tpm1 isoforms in adult hearts. We found that the RNA-binding protein RBFOX2 controls AS of rat Tpm1 exon 6a, which is important for cooperative actin binding. Furthermore, RBFOX2 regulates Tpm1 AS of exon 6a antagonistically to the RNA-binding protein PTBP1. In sum, we defined full-length Tpm1 isoforms with different exon combinations that are tightly regulated during cardiac development and provided insights into the regulation of Tpm1 AS by RNA-binding proteins. Our results demonstrate that nanopore sequencing is an excellent tool to determine full-length AS variants of muscle-enriched genes.

Mouse Adapted SARS-CoV-2 protects animals from lethal SARS-CoV challenge.

The emergence of SARS-CoV-2 has resulted in a worldwide pandemic causing significant damage to public health and the economy. Efforts to understand the mechanisms of COVID-19 disease have been hampered by the lack of robust mouse models. To overcome this barrier, we utilized a reverse genetic system to generate a mouse-adapted strain of SARS-CoV-2. Incorporating key mutations found in SARSCoV-2 variants, this model recapitulates critical elements of human infection including viral replication in the lung, immune cell infiltration, and significant in vivo disease. Importantly, mouse-adaptation of SARS-CoV-2 does not impair replication in human airway cells and maintains antigenicity similar to human SARS-CoV-2 strains. Utilizing this model, we demonstrate that SARS-CoV-2 infected mice are protected from lethal challenge with the original SARS-CoV, suggesting immunity from heterologous CoV strains. Together, the results highlight the utility of this mouse model for further study of SARS-CoV-2 infection and disease.