Vinita Puri

Zika Virus Antagonizes Type I Interferon Responses during Infection of Human Dendritic Cells

Zika virus (ZIKV) is an emerging mosquito-borne flavivirus that is causally linked to severe neonatal birth defects, including microcephaly, and is associated with Guillain-Barre syndrome in adults. Dendritic cells (DCs) are an important cell type during infection by multiple mosquito-borne flaviviruses, including dengue virus, West Nile virus, Japanese encephalitis virus, and yellow fever virus. Despite this, the interplay between ZIKV and DCs remains poorly defined. Here, we found human DCs supported productive infection by a contemporary Puerto Rican isolate with considerable variability in viral replication, but not viral binding, between DCs from different donors. Historic isolates from Africa and Asia also infected DCs with distinct viral replication kinetics between strains. African lineage viruses displayed more rapid replication kinetics and infection magnitude as compared to Asian lineage viruses, and uniquely induced cell death. Infection of DCs with both contemporary and historic ZIKV isolates led to minimal up-regulation of T cell co-stimulatory and MHC molecules, along with limited secretion of inflammatory cytokines. Inhibition of type I interferon (IFN) protein translation was observed during ZIKV infection, despite strong induction at the RNA transcript level and up-regulation of other host antiviral proteins. Treatment of human DCs with RIG-I agonist potently restricted ZIKV replication, while type I IFN had only modest effects. Mechanistically, we found all strains of ZIKV antagonized type I IFN-mediated phosphorylation of STAT1 and STAT2. Combined, our findings show that ZIKV subverts DC immunogenicity during infection, in part through evasion of type I IFN responses, but that the RLR signaling pathway is still capable of inducing an antiviral state, and therefore may serve as an antiviral therapeutic target.

Comprehensive Genome Scale Phylogenetic Study Provides New Insights on the Global Expansion of Chikungunya Virus

ABSTRACT Since the India and Indian Ocean outbreaks of 2005 and 2006, the global distribution of chikungunya virus (CHIKV) and the locations of epidemics have dramatically shifted. First, the Indian Ocean lineage (IOL) caused sustained epidemics in India and has radiated to many other countries. Second, the Asian lineage has caused frequent outbreaks in the Pacific islands and in 2013 was introduced into the Caribbean, followed by rapid spread to nearly all of the neotropics. Further, CHIKV epidemics, as well as exported cases, have been reported in central Africa after a long period of perceived silence. To understand these changes and to anticipate the future of the virus, the exact distribution, genetic diversity, transmission routes, and future epidemic potential of CHIKV require further assessment. To do so, we conducted the most comprehensive phylogenetic analysis to date, examined CHIKV evolution and transmission, and explored distinct genetic factors associated with the emergence of the East/Central/South African (ECSA) lineage, the IOL, and the Asian lineage. Our results reveal contrasting evolutionary patterns among the lineages, with growing genetic diversities observed in each, and suggest that CHIKV will continue to be a major public health threat with the potential for further emergence and spread. IMPORTANCE Chikungunya fever is a reemerging infectious disease that is transmitted by Aedes mosquitoes and causes severe health and economic burdens in affected populations. Since the unprecedented Indian Ocean and Indian subcontinent outbreaks of 2005 and 2006, CHIKV has further expanded its geographic range, including to the Americas in 2013. Its evolution and transmission during and following these epidemics, as well as the recent evolution and spread of other lineages, require optimal assessment. Using newly obtained genome sequences, we provide a comprehensive update of the global distribution of CHIKV genetic diversity and analyze factors associated with recent outbreaks. These results provide a solid foundation for future evolutionary studies of CHIKV that can elucidate emergence mechanisms and also may help to predict future epidemics.

Respiratory Syncytial Virus whole-genome sequencing identifies convergent evolution of sequence duplication in the C-terminus of the G gene.

Respiratory Syncytial Virus (RSV) is responsible for considerable morbidity and mortality worldwide and is the most important respiratory viral pathogen in infants. Extensive sequence variability within and between RSV group A and B viruses and the ability of multiple clades and sub-clades of RSV to co-circulate are likely mechanisms contributing to the evasion of herd immunity. Surveillance and large-scale whole-genome sequencing of RSV is currently limited but would help identify its evolutionary dynamics and sites of selective immune evasion. In this study, we performed complete-genome next-generation sequencing of 92 RSV isolates from infants in central Tennessee during the 2012–2014 RSV seasons. We identified multiple co-circulating clades of RSV from both the A and B groups. Each clade is defined by signature N- and O-linked glycosylation patterns. Analyses of specific RSV genes revealed high rates of positive selection in the attachment (G) gene. We identified RSV-A viruses in circulation with and without a recently reported 72-nucleotide G gene sequence duplication. Furthermore, we show evidence of convergent evolution of G gene sequence duplication and fixation over time, which suggests a potential fitness advantage of RSV with the G sequence duplication.

A Universal Next-Generation Sequencing Protocol To Generate Noninfectious Barcoded cDNA Libraries from High-Containment RNA Viruses.

This report establishes and validates a standard operating procedure (SOP) for select agents (SAs) and other biosafety level 3 and/or 4 (BSL-3/4) RNA viruses to rapidly generate noninfectious, barcoded cDNA amenable for next-generation sequencing (NGS). This eliminates the burden of testing all processed samples derived from high-consequence pathogens prior to transfer from high-containment laboratories to lower-containment facilities for sequencing. Our established protocol can be scaled up for high-throughput sequencing of hundreds of samples simultaneously, which can dramatically reduce the cost and effort required for NGS library construction. NGS data from this SOP can provide complete genome coverage from viral stocks and can also detect virus-specific reads from limited starting material. Our data suggest that the procedure can be implemented and easily validated by institutional biosafety committees across research laboratories. ABSTRACT Several biosafety level 3 and/or 4 (BSL-3/4) pathogens are high-consequence, single-stranded RNA viruses, and their genomes, when introduced into permissive cells, are infectious. Moreover, many of these viruses are select agents (SAs), and their genomes are also considered SAs. For this reason, cDNAs and/or their derivatives must be tested to ensure the absence of infectious virus and/or viral RNA before transfer out of the BSL-3/4 and/or SA laboratory. This tremendously limits the capacity to conduct viral genomic research, particularly the application of next-generation sequencing (NGS). Here, we present a sequence-independent method to rapidly amplify viral genomic RNA while simultaneously abolishing both viral and genomic RNA infectivity across multiple single-stranded positive-sense RNA (ssRNA+) virus families. The process generates barcoded DNA amplicons that range in length from 300 to 1,000 bp, which cannot be used to rescue a virus and are stable to transport at room temperature. Our barcoding approach allows for up to 288 barcoded samples to be pooled into a single library and run across various NGS platforms without potential reconstitution of the viral genome. Our data demonstrate that this approach provides full-length genomic sequence information not only from high-titer virion preparations but it can also recover specific viral sequence from samples with limited starting material in the background of cellular RNA, and it can be used to identify pathogens from unknown samples. In summary, we describe a rapid, universal standard operating procedure that generates high-quality NGS libraries free of infectious virus and infectious viral RNA. IMPORTANCE This report establishes and validates a standard operating procedure (SOP) for select agents (SAs) and other biosafety level 3 and/or 4 (BSL-3/4) RNA viruses to rapidly generate noninfectious, barcoded cDNA amenable for next-generation sequencing (NGS). This eliminates the burden of testing all processed samples derived from high-consequence pathogens prior to transfer from high-containment laboratories to lower-containment facilities for sequencing. Our established protocol can be scaled up for high-throughput sequencing of hundreds of samples simultaneously, which can dramatically reduce the cost and effort required for NGS library construction. NGS data from this SOP can provide complete genome coverage from viral stocks and can also detect virus-specific reads from limited starting material. Our data suggest that the procedure can be implemented and easily validated by institutional biosafety committees across research laboratories.

Isolation and Characterization of a Novel Gammaherpesvirus from a Microbat Cell Line

Bats are of significant interest as reservoirs for zoonotic viral pathogens; however, tools to dissect bat-virus interactions are limited in availability. This study serendipitously identified, in an established bat cell line, a fully replication-competent gammaherpesvirus; determined the complete genome sequence of the virus; and generated a viral transcript map. This virus can replicate in select human and nonhuman primate cell lines. However, analyses of viral sequences support a bat origin for this virus; we therefore refer to the virus as bat gammaherpesvirus 8 (BGHV8). The viral genome contains unique open reading frames that likely encode modulators of bat innate and adaptive immune signaling pathways and expresses viral miRNAs. The virus and its gene products should provide a unique tool to dissect both bat and gammaherpesvirus biology. ABSTRACT While employing deep sequencing and de novo assembly to characterize the mRNA transcript profile of a cell line derived from the microbat Myotis velifer incautus, we serendipitously identified mRNAs encoding proteins with a high level of identity to herpesviruses. A majority were closely related to proteins of equine herpesvirus 2 (EHV-2), a horse gammaherpesvirus. We demonstrated by electron microscopy the presence of herpesvirus-like particles in the microbat cells. Passage of supernatants from microbat cells to Vero cells resulted in syncytium formation, and expression of viral genes and amplification of viral DNA were demonstrated by quantitative PCR. Susceptibility of human cell lines to productive infection was also demonstrated. Next-generation sequencing and de novo assembly of the viral genome from supernatants from Vero cells yielded a single contig of approximately 130 kb with at least 77 open reading frames (ORFs), predicted microRNAs (miRNAs), and a gammaherpesvirus genomic organization. Phylogenic analysis of the envelope glycoprotein (gB) and DNA polymerase (POLD1) revealed similarity to multiple gammaherpesviruses, including those from as-yet-uncultured viruses of the Rhadinovirus genus that were obtained by deep sequencing of bat tissues. Moreover, the assembled genome revealed ORFs that share little or no homology to known ORFs in EHV-2 but are similar to accessory proteins of other gammaherpesviruses. Some also have striking homology to predicted Myotis bat proteins. Cumulatively, this study provides the first isolation and characterization of a replication-competent bat gammaherpesvirus. IMPORTANCE Bats are of significant interest as reservoirs for zoonotic viral pathogens; however, tools to dissect bat-virus interactions are limited in availability. This study serendipitously identified, in an established bat cell line, a fully replication-competent gammaherpesvirus; determined the complete genome sequence of the virus; and generated a viral transcript map. This virus can replicate in select human and nonhuman primate cell lines. However, analyses of viral sequences support a bat origin for this virus; we therefore refer to the virus as bat gammaherpesvirus 8 (BGHV8). The viral genome contains unique open reading frames that likely encode modulators of bat innate and adaptive immune signaling pathways and expresses viral miRNAs. The virus and its gene products should provide a unique tool to dissect both bat and gammaherpesvirus biology.

The Ebola virus VP35 protein binds viral immunostimulatory and host RNAs identified through deep sequencing

Ebola virus and Marburg virus are members of the Filovirdae family and causative agents of hemorrhagic fever with high fatality rates in humans. Filovirus virulence is partially attributed to the VP35 protein, a well-characterized inhibitor of the RIG-I-like receptor pathway that triggers the antiviral interferon (IFN) response. Prior work demonstrates the ability of VP35 to block potent RIG-I activators, such as Sendai virus (SeV), and this IFN-antagonist activity is directly correlated with its ability to bind RNA. Several structural studies demonstrate that VP35 binds short synthetic dsRNAs; yet, there are no data that identify viral immunostimulatory RNAs (isRNA) or host RNAs bound to VP35 in cells. Utilizing a SeV infection model, we demonstrate that both viral isRNA and host RNAs are bound to Ebola and Marburg VP35s in cells. By deep sequencing the purified VP35-bound RNA, we identified the SeV copy-back defective interfering (DI) RNA, previously identified as a robust RIG-I activator, as the isRNA bound by multiple filovirus VP35 proteins, including the VP35 protein from the West African outbreak strain (Makona EBOV). Moreover, RNAs isolated from a VP35 RNA-binding mutant were not immunostimulatory and did not include the SeV DI RNA. Strikingly, an analysis of host RNAs bound by wild-type, but not mutant, VP35 revealed that select host RNAs are preferentially bound by VP35 in cell culture. Taken together, these data support a model in which VP35 sequesters isRNA in virus-infected cells to avert RIG-I like receptor (RLR) activation.

Analysis of the Aedes albopictus C6/36 genome provides insight into cell line utility for viral propagation

Abstract Background The 50-year-old Aedes albopictus C6/36 cell line is a resource for the detection, amplification, and analysis of mosquito-borne viruses including Zika, dengue, and chikungunya. The cell line is derived from an unknown number of larvae from an unspecified strain of Aedes albopictus mosquitoes. Toward improved utility of the cell line for research in virus transmission, we present an annotated assembly of the C6/36 genome. Results The C6/36 genome assembly has the largest contig N50 (3.3 Mbp) of any mosquito assembly, presents the sequences of both haplotypes for most of the diploid genome, reveals independent null mutations in both alleles of the Dicer locus, and indicates a male-specific genome. Gene annotation was computed with publicly available mosquito transcript sequences. Gene expression data from cell line RNA sequence identified enrichment of growth-related pathways and conspicuous deficiency in aquaporins and inward rectifier K+ channels. As a test of utility, RNA sequence data from Zika-infected cells were mapped to the C6/36 genome and transcriptome assemblies. Host subtraction reduced the data set by 89%, enabling faster characterization of nonhost reads. Conclusions The C6/36 genome sequence and annotation should enable additional uses of the cell line to study arbovirus vector interactions and interventions aimed at restricting the spread of human disease.

Whole genome sequencing, variant analysis, phylogenetics, and deep sequencing of Zika virus strains

The recent emergence of Zika virus (ZIKV) has been concentrated in the Caribbean, Southeastern United States, and South- and Central America; resulting in travel-based cases being reported around the globe. As multi-disciplinary collaborations are combatting the ZIKV outbreak, the need to validate the sequence of existing strains has become apparent. Here, we report high-quality sequence data for multiple ZIKV strains made publicly available through the National Institutes of Health- (NIH) funded biorepository, BEI Resources (www.beiresources.org). Next-generation sequencing, 3′ rapid amplification of cDNA ends (RACE), and viral genome annotation pipelines generated GenBank sequence records for 16 BEI Resources strains. Minor variants, consensus mutations, and consensus insertions/deletions were identified within the viral stocks using next-generation sequencing (NGS) and consensus changes were confirmed with Sanger sequencing. Bioinformatics analyses of the sequencing results confirm that the virus stocks available to the scientific research community through BEI Resources adequately represent the viral population diversity of ZIKV.

Differing epidemiological dynamics of Chikungunya virus in the Americas during the 2014-2015 epidemic

Chikungunya virus (CHIKV) has been detected sporadically since the 1950s and includes three distinct co-circulating genotypes. In late 2013, the Asian genotype of CHIKV was responsible for the Caribbean outbreak (CO) that rapidly became an epidemic throughout the Americas. There is a limited understanding of the molecular evolution of CHIKV in the Americas during this epidemic. We sequenced 185 complete CHIKV genomes collected mainly from Nicaragua in Central America and Florida in the United States during the 2014–2015 Caribbean/Americas epidemic. Our comprehensive phylogenetic analyses estimated the epidemic history of the Asian genotype and the recent Caribbean outbreak (CO) clade, revealed considerable genetic diversity within the CO clade, and described different epidemiological dynamics of CHIKV in the Americas. Specifically, we identified multiple introductions in both Nicaragua and Florida, with rapid local spread of viruses in Nicaragua but limited autochthonous transmission in Florida in the US. Our phylogenetic analysis also showed phylogeographic clustering of the CO clade. In addition, we identified the significant amino acid substitutions that were observed across the entire Asian genotype during its evolution and examined amino acid changes that were specific to the CO clade. Deep sequencing analysis identified specific minor variants present in clinical specimens below-consensus levels. Finally, we investigated the association between viral phylogeny and geographic/clinical metadata in Nicaragua. To date, this study represents the largest single collection of CHIKV complete genomes during the Caribbean/Americas epidemic and significantly expands our understanding of the emergence and evolution of CHIKV CO clade in the Americas.

Whole-Genome Sequences of Zika Virus FLR Strains after Passage in Vero or C6/36 Cells

ABSTRACT We report 26 complete genomes of Zika virus (ZIKV) isolated after passaging the Zika virus strain FLR in mosquito (C6/36) and mammalian (Vero) cell lines. The consensus ZIKV genomes we recovered show greater than 99% nucleotide identify with each other and with the FLR strain used as input.