Single-Cell Profiling of Ebola Virus Disease In Vivo Reveals Viral and Host Dynamics.

Ebola virus (EBOV) causes epidemics with high mortality yet remains understudied due to the challenge of experimentation in high-containment and outbreak settings. Here, we used single-cell transcriptomics and CyTOF-based single-cell protein quantification to characterize peripheral immune cells during EBOV infection in rhesus monkeys. We obtained 100,000 transcriptomes and 15,000,000 protein profiles, finding that immature, proliferative monocyte-lineage cells with reduced antigen-presentation capacity replace conventional monocyte subsets, while lymphocytes upregulate apoptosis genes and decline in abundance. By quantifying intracellular viral RNA, we identify molecular determinants of tropism among circulating immune cells and examine temporal dynamics in viral and host gene expression. Within infected cells, EBOV downregulates STAT1 mRNA and interferon signaling, and it upregulates putative pro-viral genes (e.g., DYNLL1 and HSPA5), nominating pathways the virus manipulates for its replication. This study sheds light on EBOV tropism, replication dynamics, and elicited immune response and provides a framework for characterizing host-virus interactions under maximum containment.

SARS-CoV-2 Disrupts Splicing, Translation, and Protein Trafficking to Suppress Host Defenses.

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) is a recently identified coronavirus that causes the respiratory disease known as coronavirus disease 2019 (COVID-19). Despite the urgent need, we still do not fully understand the molecular basis of SARS-CoV-2 pathogenesis. Here, we comprehensively define the interactions between SARS-CoV-2 proteins and human RNAs. NSP16 binds to the mRNA recognition domains of the U1 and U2 splicing RNAs and acts to suppress global mRNA splicing upon SARS-CoV-2 infection. NSP1 binds to 18S ribosomal RNA in the mRNA entry channel of the ribosome and leads to global inhibition of mRNA translation upon infection. Finally, NSP8 and NSP9 bind to the 7SL RNA in the signal recognition particle and interfere with protein trafficking to the cell membrane upon infection. Disruption of each of these essential cellular functions acts to suppress the interferon response to viral infection. Our results uncover a multipronged strategy utilized by SARS-CoV-2 to antagonize essential cellular processes to suppress host defenses.

Ebola Virus Glycoprotein with Increased Infectivity Dominated the 2013-2016 Epidemic.

The magnitude of the 2013-2016 Ebola virus disease (EVD) epidemic enabled an unprecedented number of viral mutations to occur over successive human-to-human transmission events, increasing the probability that adaptation to the human host occurred during the outbreak. We investigated one nonsynonymous mutation, Ebola virus (EBOV) glycoprotein (GP) mutant A82V, for its effect on viral infectivity. This mutation, located at the NPC1-binding site on EBOV GP, occurred early in the 2013-2016 outbreak and rose to high frequency. We found that GP-A82V had heightened ability to infect primate cells, including human dendritic cells. The increased infectivity was restricted to cells that have primate-specific NPC1 sequences at the EBOV interface, suggesting that this mutation was indeed an adaptation to the human host. GP-A82V was associated with increased mortality, consistent with the hypothesis that the heightened intrinsic infectivity of GP-A82V contributed to disease severity during the EVD epidemic.

Ebola Virus Epidemiology, Transmission, and Evolution during Seven Months in Sierra Leone.

The 2013-2015 Ebola virus disease (EVD) epidemic is caused by the Makona variant of Ebola virus (EBOV). Early in the epidemic, genome sequencing provided insights into virus evolution and transmission and offered important information for outbreak response. Here, we analyze sequences from 232 patients sampled over 7 months in Sierra Leone, along with 86 previously released genomes from earlier in the epidemic. We confirm sustained human-to-human transmission within Sierra Leone and find no evidence for import or export of EBOV across national borders after its initial introduction. Using high-depth replicate sequencing, we observe both host-to-host transmission and recurrent emergence of intrahost genetic variants. We trace the increasing impact of purifying selection in suppressing the accumulation of nonsynonymous mutations over time. Finally, we note changes in the mucin-like domain of EBOV glycoprotein that merit further investigation. These findings clarify the movement of EBOV within the region and describe viral evolution during prolonged human-to-human transmission.

Clinical Sequencing Uncovers Origins and Evolution of Lassa Virus.

The 2013-2015 West African epidemic of Ebola virus disease (EVD) reminds us of how little is known about biosafety level 4 viruses. Like Ebola virus, Lassa virus (LASV) can cause hemorrhagic fever with high case fatality rates. We generated a genomic catalog of almost 200 LASV sequences from clinical and rodent reservoir samples. We show that whereas the 2013-2015 EVD epidemic is fueled by human-to-human transmissions, LASV infections mainly result from reservoir-to-human infections. We elucidated the spread of LASV across West Africa and show that this migration was accompanied by changes in LASV genome abundance, fatality rates, codon adaptation, and translational efficiency. By investigating intrahost evolution, we found that mutations accumulate in epitopes of viral surface proteins, suggesting selection for immune escape. This catalog will serve as a foundation for the development of vaccines and diagnostics. VIDEO ABSTRACT.

Programmable Inhibition and Detection of RNA Viruses Using Cas13.

The CRISPR effector Cas13 could be an effective antiviral for single-stranded RNA (ssRNA) viruses because it programmably cleaves RNAs complementary to its CRISPR RNA (crRNA). Here, we computationally identify thousands of potential Cas13 crRNA target sites in hundreds of ssRNA viral species that can potentially infect humans. We experimentally demonstrate Cas13's potent activity against three distinct ssRNA viruses: lymphocytic choriomeningitis virus (LCMV); influenza A virus (IAV); and vesicular stomatitis virus (VSV). Combining this antiviral activity with Cas13-based diagnostics, we develop Cas13-assisted restriction of viral expression and readout (CARVER), an end-to-end platform that uses Cas13 to detect and destroy viral RNA. We further screen hundreds of crRNAs along the LCMV genome to evaluate how conservation and target RNA nucleotide content influence Cas13's antiviral activity. Our results demonstrate that Cas13 can be harnessed to target a wide range of ssRNA viruses and CARVER's potential broad utility for rapid diagnostic and antiviral drug development.

Engineered chimeric receptors for dissecting interferon signaling.

Though interferons (IFNs) were once heralded as panaceas to numerous diseases, how cells decode varying IFN stimuli and subsequently produce (in)appropriate signaling remain unclear. Our labs recently engineered novel erythropoietin receptor-IFN chimeric receptors, and we highlight their utility in two cases uncovering differential genetic determinants of type I (IFN-α/ß) and type III (IFN-λ) IFN signaling. These and other types of synthetic (cytokine) receptors could be expanded to real-time signaling dynamics and in vivo studies.

Zika virus evolution and spread in the Americas.

Although the recent Zika virus (ZIKV) epidemic in the Americas and its link to birth defects have attracted a great deal of attention, much remains unknown about ZIKV disease epidemiology and ZIKV evolution, in part owing to a lack of genomic data. Here we address this gap in knowledge by using multiple sequencing approaches to generate 110 ZIKV genomes from clinical and mosquito samples from 10 countries and territories, greatly expanding the observed viral genetic diversity from this outbreak. We analysed the timing and patterns of introductions into distinct geographic regions; our phylogenetic evidence suggests rapid expansion of the outbreak in Brazil and multiple introductions of outbreak strains into Puerto Rico, Honduras, Colombia, other Caribbean islands, and the continental United States. We find that ZIKV circulated undetected in multiple regions for many months before the first locally transmitted cases were confirmed, highlighting the importance of surveillance of viral infections. We identify mutations with possible functional implications for ZIKV biology and pathogenesis, as well as those that might be relevant to the effectiveness of diagnostic tests.

Evidence of Ebola Virus Replication and High Concentration in Semen of a Patient During Recovery.

In one patient over time, we found that concentration of Ebola virus RNA in semen during recovery is remarkably higher than blood at peak illness. Virus in semen is replication-competent with no change in viral genome over time. Presence of sense RNA suggests replication in cells present in semen.

In vivo Ebola virus infection leads to a strong innate response in circulating immune cells.

BACKGROUND: Ebola virus is the causative agent of a severe syndrome in humans with a fatality rate that can approach 90 %. During infection, the host immune response is thought to become dysregulated, but the mechanisms through which this happens are not entirely understood. In this study, we analyze RNA sequencing data to determine the host response to Ebola virus infection in circulating immune cells. RESULTS: Approximately half of the 100 genes with the strongest early increases in expression were interferon-stimulated genes, such as ISG15, OAS1, IFIT2, HERC5, MX1 and DHX58. Other highly upregulated genes included cytokines CXCL11, CCL7, IL2RA, IL2R1, IL15RA, and CSF2RB, which have not been previously reported to change during Ebola virus infection. Comparing this response in two different models of exposure (intramuscular and aerosol) revealed a similar signature of infection. The strong innate response in the aerosol model was seen not only in circulating cells, but also in primary and secondary target tissues. Conversely, the innate immune response of vaccinated macaques was almost non-existent. This suggests that the innate response is a major aspect of the cellular response to Ebola virus infection in multiple tissues. CONCLUSIONS: Ebola virus causes a severe infection in humans that is associated with high mortality. The host immune response to virus infection is thought to be an important aspect leading to severe pathology, but the components of this overactive response are not well characterized. Here, we analyzed how circulating immune cells respond to the virus and found that there is a strong innate response dependent on active virus replication. This finding is in stark contrast to in vitro evidence showing a suppression of innate immune signaling, and it suggests that the strong innate response we observe in infected animals may be an important contributor to pathogenesis.

DNA methyltransferase DNMT3A associates with viral proteins and impacts HSV-1 infection.

Viral infections can alter the cellular epigenetic landscape, through modulation of either DNA methylation profiles or chromatin remodeling enzymes and histone modifications. These changes can act to promote viral replication or host defense. Herpes simplex virus type 1 (HSV-1) is a prominent human pathogen, which relies on interactions with host factors for efficient replication and spread. Nevertheless, the knowledge regarding its modulation of epigenetic factors remains limited. Here, we used fluorescently-labeled viruses in conjunction with immunoaffinity purification and MS to study virus-virus and virus-host protein interactions during HSV-1 infection in primary human fibroblasts. We identified interactions among viral capsid and tegument proteins, detecting phosphorylation of the capsid protein VP26 at sites within its UL37-binding domain, and an acetylation within the major capsid protein VP5. Interestingly, we found a nuclear association between viral capsid proteins and the de novo DNA methyltransferase DNA (cytosine-5)-methyltransferase 3A (DNMT3A), which we confirmed by reciprocal isolations and microscopy. We show that drug-induced inhibition of DNA methyltransferase activity, as well as siRNA- and shRNA-mediated DNMT3A knockdowns trigger reductions in virus titers. Altogether, our results highlight a functional association of viral proteins with the mammalian DNA methyltransferase machinery, pointing to DNMT3A as a host factor required for effective HSV-1 infection.

A proteomic perspective of inbuilt viral protein regulation: pUL46 tegument protein is targeted for degradation by ICP0 during herpes simplex virus type 1 infection.

Much like the host cells they infect, viruses must also regulate their life cycles. Herpes simples virus type 1 (HSV-1), a prominent human pathogen, uses a promoter-rich genome in conjunction with multiple viral trans-activating factors. Following entry into host cells, the virion-associated outer tegument proteins pUL46 and pUL47 act to increase expression of viral immediate-early (α) genes, thereby helping initiate the infection life cycle. Because pUL46 has gone largely unstudied, we employed a hybrid mass spectrometry-based approach to determine how pUL46 exerts its functions during early stages of infection. For a spatio-temporal characterization of pUL46, time-lapse microscopy was performed in live cells to define its dynamic localization from 2 to 24 h postinfection. Next, pUL46-containing protein complexes were immunoaffinity purified during infection of human fibroblasts and analyzed by mass spectrometry to investigate virus-virus and virus-host interactions, as well as post-translational modifications. We demonstrated that pUL46 is heavily phosphorylated in at least 23 sites. One phosphorylation site matched the consensus 14-3-3 phospho-binding motif, consistent with our identification of 14-3-3 proteins and host and viral kinases as specific pUL46 interactions. Moreover, we determined that pUL46 specifically interacts with the viral E3 ubiquitin ligase ICP0. We demonstrated that pUL46 is partially degraded in a proteasome-mediated manner during infection, and that the catalytic activity of ICP0 is responsible for this degradation. This is the first evidence of a viral protein being targeted for degradation by another viral protein during HSV-1 infection. Together, these data indicate that pUL46 levels are tightly controlled and important for the temporal regulation of viral gene expression throughout the virus life cycle. The concept of a structural virion protein, pUL46, performing nonstructural roles is likely to reflect a theme common to many viruses, and a better understanding of these functions will be important for developing therapeutics.

Liver-specific Mettl14 deletion induces nuclear heterotypia and dysregulates RNA export machinery.

Modification of RNA with N6-methyladenosine (m6A) has gained attention in recent years as a general mechanism of gene regulation. In the liver, m6A, along with its associated machinery, has been studied as a potential biomarker of disease and cancer, with impacts on metabolism, cell cycle regulation, and pro-cancer state signaling. However these observational data have yet to be causally examined in vivo. For example, neither perturbation of the key m6A writers Mettl3 and Mettl14, nor the m6A readers Ythdf1 and Ythdf2 have been thoroughly mechanistically characterized in vivo as they have been in vitro. To understand the functions of these machineries, we developed mouse models and found that deleting Mettl14 led to progressive liver injury characterized by nuclear heterotypia, with changes in mRNA splicing, processing and export leading to increases in mRNA surveillance and recycling.

Structural basis of specificity and cross-reactivity in T cell receptors specific for cytochrome c-I-E(k).

T cells specific for the cytochrome c Ag are widely used to investigate many aspects of TCR specificity and interactions with peptide-MHC, but structural information has long been elusive. In this study, we present structures for the well-studied 2B4 TCR, as well as a naturally occurring variant of the 5c.c7 TCR, 226, which is cross-reactive with more than half of possible substitutions at all three TCR-sensitive residues on the peptide Ag. These structures alone and in complex with peptide-MHC ligands allow us to reassess many prior mutagenesis results. In addition, the structure of 226 bound to one peptide variant, p5E, shows major changes in the CDR3 contacts compared with wild-type, yet the TCR V-region contacts with MHC are conserved. These and other data illustrate the ability of TCRs to accommodate large variations in CDR3 structure and peptide contacts within the constraints of highly conserved TCR-MHC interactions.

Membrane-proximal motifs encode differences in signaling strength between type I and III interferon receptors.

Interferons (IFNs) play crucial roles in antiviral defenses. Despite using the same Janus-activated kinase (JAK)-signal transducer and activator of transcription (STAT) signaling cascade, type I and III IFN receptors differ in the magnitude and dynamics of their signaling in terms of STAT phosphorylation, gene transcription, and antiviral responses. These differences are not due to ligand-binding affinity and receptor abundance. Here, we investigated the ability of the intracellular domains (ICDs) of IFN receptors to differentiate between type I and III IFN signaling. We engineered synthetic, heterodimeric type I and III IFN receptors that were stably expressed at similar amounts in human cells and responded to a common ligand. We found that our synthetic type I IFN receptors stimulated STAT phosphorylation and gene expression to greater extents than did the corresponding type III IFN receptors. Furthermore, we identified short "box motifs" within ICDs that bind to JAK1 that were sufficient to encode differences between the type I and III IFN receptors. Together, our results indicate that specific regions within the ICDs of IFN receptor subunits encode different downstream signaling strengths that enable type I and III IFN receptors to produce distinct signaling outcomes.

Single-cell massively-parallel multiplexed microbial sequencing (M3-seq) identifies rare bacterial populations and profiles phage infection.

Bacterial populations are highly adaptive. They can respond to stress and survive in shifting environments. How the behaviours of individual bacteria vary during stress, however, is poorly understood. To identify and characterize rare bacterial subpopulations, technologies for single-cell transcriptional profiling have been developed. Existing approaches show some degree of limitation, for example, in terms of number of cells or transcripts that can be profiled. Due in part to these limitations, few conditions have been studied with these tools. Here we develop massively-parallel, multiplexed, microbial sequencing (M3-seq)-a single-cell RNA-sequencing platform for bacteria that pairs combinatorial cell indexing with post hoc rRNA depletion. We show that M3-seq can profile bacterial cells from different species under a range of conditions in single experiments. We then apply M3-seq to hundreds of thousands of cells, revealing rare populations and insights into bet-hedging associated with stress responses and characterizing phage infection.

Single-cell profiling of lncRNA expression during Ebola virus infection in rhesus macaques.

Long non-coding RNAs (lncRNAs) are involved in numerous biological processes and are pivotal mediators of the immune response, yet little is known about their properties at the single-cell level. Here, we generate a multi-tissue bulk RNAseq dataset from Ebola virus (EBOV) infected and not-infected rhesus macaques and identified 3979 novel lncRNAs. To profile lncRNA expression dynamics in immune circulating single-cells during EBOV infection, we design a metric, Upsilon, to estimate cell-type specificity. Our analysis reveals that lncRNAs are expressed in fewer cells than protein-coding genes, but they are not expressed at lower levels nor are they more cell-type specific when expressed in the same number of cells. In addition, we observe that lncRNAs exhibit similar changes in expression patterns to those of protein-coding genes during EBOV infection, and are often co-expressed with known immune regulators. A few lncRNAs change expression specifically upon EBOV entry in the cell. This study sheds light on the differential features of lncRNAs and protein-coding genes and paves the way for future single-cell lncRNA studies.

Natural history of Ebola virus disease in rhesus monkeys shows viral variant emergence dynamics and tissue-specific host responses.

Ebola virus (EBOV) causes Ebola virus disease (EVD), marked by severe hemorrhagic fever; however, the mechanisms underlying the disease remain unclear. To assess the molecular basis of EVD across time, we performed RNA sequencing on 17 tissues from a natural history study of 21 rhesus monkeys, developing new methods to characterize host-pathogen dynamics. We identified alterations in host gene expression with previously unknown tissue-specific changes, including downregulation of genes related to tissue connectivity. EBOV was widely disseminated throughout the body; using a new, broadly applicable deconvolution method, we found that viral load correlated with increased monocyte presence. Patterns of viral variation between tissues differentiated primary infections from compartmentalized infections, and several variants impacted viral fitness in a EBOV/Kikwit minigenome system, suggesting that functionally significant variants can emerge during early infection. This comprehensive portrait of host-pathogen dynamics in EVD illuminates new features of pathogenesis and establishes resources to study other emerging pathogens.

Phylogenetic analysis of SARS-CoV-2 in Boston highlights the impact of superspreading events.

Analysis of 772 complete severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) genomes from early in the Boston-area epidemic revealed numerous introductions of the virus, a small number of which led to most cases. The data revealed two superspreading events. One, in a skilled nursing facility, led to rapid transmission and significant mortality in this vulnerable population but little broader spread, whereas other introductions into the facility had little effect. The second, at an international business conference, produced sustained community transmission and was exported, resulting in extensive regional, national, and international spread. The two events also differed substantially in the genetic variation they generated, suggesting varying transmission dynamics in superspreading events. Our results show how genomic epidemiology can help to understand the link between individual clusters and wider community spread.

Reporter Assays for Ebola Virus Nucleoprotein Oligomerization, Virion-Like Particle Budding, and Minigenome Activity Reveal the Importance of Nucleoprotein Amino Acid Position 111.

For highly pathogenic viruses, reporter assays that can be rapidly performed are critically needed to identify potentially functional mutations for further study under maximal containment (e.g., biosafety level 4 [BSL-4]). The Ebola virus nucleoprotein (NP) plays multiple essential roles during the viral life cycle, yet few tools exist to study the protein under BSL-2 or equivalent containment. Therefore, we adapted reporter assays to measure NP oligomerization and virion-like particle (VLP) production in live cells and further measured transcription and replication using established minigenome assays. As a proof-of-concept, we examined the NP-R111C substitution, which emerged during the 2013‒2016 Western African Ebola virus disease epidemic and rose to high frequency. NP-R111C slightly increased NP oligomerization and VLP budding but slightly decreased transcription and replication. By contrast, a synthetic charge-reversal mutant, NP-R111E, greatly increased oligomerization but abrogated transcription and replication. These results are intriguing in light of recent structures of NP oligomers, which reveal that the neighboring residue, K110, forms a salt bridge with E349 on adjacent NP molecules. By developing and utilizing multiple reporter assays, we find that the NP-111 position mediates a complex interplay between NP's roles in protein structure, virion budding, and transcription and replication.