Search

Abstract To understand the roles of acute-phase viral dynamics and host immune responses in post-acute sequelae of SARS-CoV-2 infection (PASC), we enrolled 136 participants within 5 days of their first positive SARS-CoV-2 real-time PCR test. Participants self-collected up to 21 nasal specimens within the first 28 days post-symptom onset; interviewer-administered questionnaires and blood samples were collected at enrollment, days 9, 14, 21, 28, and month 4 and 8 post-symptom onset. Defining PASC as the presence of any COVID-associated symptom at their 4-month visit, we compared viral markers (quantity and duration of nasal viral RNA load, infectious viral load, and plasma N-antigen level) and host immune markers (IL-6, IL-10, TNF-α, IFN-α, IFN-γ, MCP, IP-10, and Spike IgG) over the acute period. Compared to those who fully recovered, those reporting PASC demonstrated significantly higher maximum levels of SARS-CoV-2 RNA and N-antigen, burden of RNA and infectious viral shedding, and lower Spike-specific IgG levels within 9 days post-illness onset. No significant differences were identified among a panel of host immune markers. Our results suggest early viral dynamics and the associated host immune responses play a role in the pathogenesis of PASC, highlighting the importance of understanding early biological markers in the natural history of PASC.

Abstract A central role of viral capsids is to protect the viral genome from the harsh extracellular environment while facilitating initiation of infection when the virus encounters a target cell. Viruses are thought to have evolved an optimal equilibrium between particle stability and efficiency of cell entry. In this study, we genetically perturb this equilibrium in a non-enveloped virus, enterovirus A71 to determine its structural basis. We isolate a single-point mutation variant with increased particle thermotolerance and decreased efficiency of cell entry. Using cryo-electron microscopy and molecular dynamics simulations, we determine that the thermostable native particles have acquired an expanded conformation that results in a significant increase in protein dynamics. Examining the intermediate states of the thermostable variant reveals a potential pathway for uncoating. We propose a sequential release of the lipid pocket factor, followed by internal VP4 and ultimately the viral RNA.

A better mechanistic understanding of virus-host dependencies can help reveal vulnerabilities and identify opportunities for therapeutic intervention. Of particular interest are essential interactions that enable production of viral proteins, as those could target an early step in the virus lifecycle. Here, we use subcellular proteomics, ribosome profiling analyses and reporter assays to detect changes in protein synthesis dynamics during SARS-CoV-2 (CoV2) infection. We identify specific translation factors and molecular chaperones that are used by CoV2 to promote the synthesis and maturation of its own proteins. These can be targeted to inhibit infection, without major toxicity to the host. We also find that CoV2 non-structural protein 1 (Nsp1) cooperates with initiation factors EIF1 and 1A to selectively enhance translation of viral RNA. When EIF1/1A are depleted, more ribosomes initiate translation from a conserved upstream CUG start codon found in all genomic and subgenomic viral RNAs. This results in higher translation of an upstream open reading frame (uORF1) and lower translation of the main ORF, altering the stoichiometry of viral proteins and attenuating infection. Replacing the upstream CUG with AUG strongly inhibits translation of the main ORF independently of Nsp1, EIF1, or EIF1A. Taken together, our work describes multiple dependencies of CoV2 on host biosynthetic networks and proposes a model for dosage control of viral proteins through Nsp1-mediated control of translation start site selection.

ABSTRACT Wastewater-based epidemiology (WBE) emerged during the coronavirus disease 2019 (COVID-19) pandemic as a scalable and broadly applicable method for community-level monitoring of infectious disease burden. The lack of high-resolution fecal shedding data for severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) limits our ability to link WBE measurements to disease burden. In this study, we present longitudinal, quantitative fecal shedding data for SARS-CoV-2 RNA, as well as for the commonly used fecal indicators pepper mild mottle virus (PMMoV) RNA and crAss-like phage (crAssphage) DNA. The shedding trajectories from 48 SARS-CoV-2-infected individuals suggest a highly individualized, dynamic course of SARS-CoV-2 RNA fecal shedding. Of the individuals that provided at least three stool samples spanning more than 14 days, 77% had one or more samples that tested positive for SARS-CoV-2 RNA. We detected PMMoV RNA in at least one sample from all individuals and in 96% (352/367) of samples overall. CrAssphage DNA was detected in at least one sample from 80% (38/48) of individuals and was detected in 48% (179/371) of all samples. The geometric mean concentrations of PMMoV and crAssphage in stool across all individuals were 8.7 × 104 and 1.4 × 104 gene copies/milligram-dry weight, respectively, and crAssphage shedding was more consistent for individuals than PMMoV shedding. These results provide us with a missing link needed to connect laboratory WBE results with mechanistic models, and this will aid in more accurate estimates of COVID-19 burden in sewersheds. Additionally, the PMMoV and crAssphage data are critical for evaluating their utility as fecal strength normalizing measures and for source-tracking applications. IMPORTANCE This research represents a critical step in the advancement of wastewater monitoring for public health. To date, mechanistic materials balance modeling of wastewater-based epidemiology has relied on SARS-CoV-2 fecal shedding estimates from small-scale clinical reports or meta-analyses of research using a wide range of analytical methodologies. Additionally, previous SARS-CoV-2 fecal shedding data have not contained sufficient methodological information for building accurate materials balance models. Like SARS-CoV-2, fecal shedding of PMMoV and crAssphage has been understudied to date. The data presented here provide externally valid and longitudinal fecal shedding data for SARS-CoV-2, PMMoV, and crAssphage which can be directly applied to WBE models and ultimately increase the utility of WBE.

Abstract Novel oral poliovirus vaccine type 2 (nOPV2) is being developed to reduce the rare occurrence of disease and outbreaks associated with the genetic instability of the Sabin vaccine strains. Children aged 1 to 5 years were enrolled in two related clinical studies to assess safety, immunogenicity, shedding rates and properties of the shed virus following vaccination with nOPV2 (two candidates) versus traditional Sabin OPV type 2 (mOPV2). The anticipated pattern of reversion and increased virulence was observed for shed Sabin-2 virus, as assessed using a mouse model of poliovirus neurovirulence. In contrast, there were significantly reduced odds of mouse paralysis for shed virus for both nOPV2 candidates when compared to shed Sabin-2 virus. Next-generation sequencing of shed viral genomes was consistent with and further supportive of the observed neurovirulence associated with shed Sabin-2 virus, as well as the reduced reversion to virulence of shed candidate viruses. While shed Sabin-2 showed anticipated A481G reversion in the primary attenuation site in domain V in the 5’ untranslated region to be associated with increased mouse paralysis, the stabilized domain V in the candidate viruses did not show polymorphisms consistent with reversion to neurovirulence. The available data from a key target age group for outbreak response confirm the superior genetic and phenotypic stability of shed nOPV2 strains compared to shed Sabin-2 and suggest that nOPV2 should be associated with less paralytic disease and potentially a lower risk of seeding new outbreaks.

Viruses have evolved mechanisms to modulate cellular pathways to facilitate infection. One such pathway is the formation of stress granules (SG), which are ribonucleoprotein complexes that assemble during translation inhibition following cellular stress. Inhibition of SG assembly has been observed under numerous virus infections across species, suggesting a conserved fundamental viral strategy. However, the significance of SG modulation during virus infection is not fully understood. The 1A protein encoded by the model dicistrovirus, Cricket paralysis virus (CrPV), is a multifunctional protein that can bind to and degrade Ago-2 in an E3 ubiquitin ligase-dependent manner to block the antiviral RNA interference pathway and inhibit SG formation. Moreover, the R146 residue of 1A is necessary for SG inhibition and CrPV infection in both Drosophila S2 cells and adult flies. Here, we uncoupled CrPV-1A's functions and provide insight into its underlying mechanism for SG inhibition. CrPV-1A mediated inhibition of SGs requires the E3 ubiquitin-ligase binding domain and the R146 residue, but not the Ago-2 binding domain. Wild-type but not mutant CrPV-1A R146A localizes to the nuclear membrane which correlates with nuclear enrichment of poly(A)+ RNA. Transcriptome changes in CrPV-infected cells are dependent on the R146 residue. Finally, Nup358/RanBP2 is targeted and degraded in CrPV-infected cells in an R146-dependent manner and the depletion of Nup358 blocks SG formation. We propose that CrPV utilizes a multiprong strategy whereby the CrPV-1A protein interferes with a nuclear event that contributes to SG inhibition in order to promote infection.

Abstract COVID-19’s high virus transmission rates have caused a pandemic that is exacerbated by the high rates of asymptomatic and presymptomatic infections. These factors suggest that face masks and social distance could be paramount in containing the pandemic. We examined the efficacy of each measure and the combination of both measures using an agent-based model within a closed space that approximated real-life interactions. By explicitly considering different fractions of asymptomatic individuals, as well as a realistic hypothesis of face masks protection during inhaling and exhaling, our simulations demonstrate that a synergistic use of face masks and social distancing is the most effective intervention to curb the infection spread. To control the pandemic, our models suggest that high adherence to social distance is necessary to curb the spread of the disease, and that wearing face masks provides optimal protection even if only a small portion of the population comply with social distance. Finally, the face mask effectiveness in curbing the viral spread is not reduced if a large fraction of population is asymptomatic. Our findings have important implications for policies that dictate the reopening of social gatherings.

The impact of vaccination on SARS-CoV-2 infectiousness is not well understood. We compared longitudinal viral shedding dynamics in unvaccinated and fully vaccinated adults. SARS-CoV-2-infected adults were enrolled within 5 days of symptom onset and nasal specimens were self-collected daily for two weeks and intermittently for an additional two weeks. SARS-CoV-2 RNA load and infectious virus were analyzed relative to symptom onset stratified by vaccination status. We tested 1080 nasal specimens from 52 unvaccinated adults enrolled in the pre-Delta period and 32 fully vaccinated adults with predominantly Delta infections. While we observed no differences by vaccination status in maximum RNA levels, maximum infectious titers and the median duration of viral RNA shedding, the rate of decay from the maximum RNA load was faster among vaccinated; maximum infectious titers and maximum RNA levels were highly correlated. Furthermore, amongst participants with infectious virus, median duration of infectious virus detection was reduced from 7.5 days (IQR: 6.0-9.0) in unvaccinated participants to 6 days (IQR: 5.0-8.0) in those vaccinated (P = 0.02). Accordingly, the odds of shedding infectious virus from days 6 to 12 post-onset were lower among vaccinated participants than unvaccinated participants (OR 0.42 95% CI 0.19-0.89). These results indicate that vaccination had reduced the probability of shedding infectious virus after 5 days from symptom onset.

Summary: Understanding viral tropism is an essential step toward reducing severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) transmission, decreasing mortality from coronavirus disease 2019 (COVID-19) and limiting opportunities for mutant strains to arise. Currently, little is known about the extent to which distinct tissue sites in the human head and neck region and proximal respiratory tract selectively permit SARS-CoV-2 infection and replication. In this translational study, we discover key variabilities in expression of angiotensin-converting enzyme 2 (ACE2) and transmembrane serine protease 2 (TMPRSS2), essential SARS-CoV-2 entry factors, among the mucosal tissues of the human proximal airways. We show that SARS-CoV-2 infection is present in all examined head and neck tissues, with a notable tropism for the nasal cavity and tracheal mucosa. Finally, we uncover an association between smoking and higher SARS-CoV-2 viral infection in the human proximal airway, which may explain the increased susceptibility of smokers to developing severe COVID-19. This is at least partially explained by differences in interferon (IFN)-β1 levels between smokers and non-smokers.

The respiratory disease COVID-19 is caused by the coronavirus SARS-CoV-2. Here we report the discovery of ethacridine as a potent drug against SARS-CoV-2 (EC50 ~ 0.08 μM). Ethacridine was identified via high-throughput screening of an FDA-approved drug library in living cells using a fluorescence assay. Plaque assays, RT-PCR and immunofluorescence imaging at various stages of viral infection demonstrate that the main mode of action of ethacridine is through inactivation of viral particles, preventing their binding to the host cells. Consistently, ethacridine is effective in various cell types, including primary human nasal epithelial cells that are cultured in an air-liquid interface. Taken together, our work identifies a promising, potent, and new use of the old drug via a distinct mode of action for inhibiting SARS-CoV-2.

During replication, RNA viruses accumulate genome alterations, such as mutations and deletions. The interactions between individual variants can determine the fitness of the virus population and, thus, the outcome of infection. To investigate the effects of defective interfering genomes (DI) on wild-type (WT) poliovirus replication, we developed an ordinary differential equation model, which enables exploring the parameter space of the WT and DI competition. We also experimentally examined virus and DI replication kinetics during co-infection, and used these data to infer model parameters. Our model identifies, and our experimental measurements confirm, that the efficiencies of DI genome replication and encapsidation are two most critical parameters determining the outcome of WT replication. However, an equilibrium can be established which enables WT to replicate, albeit to reduced levels.

A novel, genetically-stabilized type 2 oral polio vaccine (nOPV2), developed to assist in the global polio eradication program, was recently the first-ever vaccine granted Emergency Use Listing by the WHO. Lot release tests for this vaccine included—for the first time to our knowledge—the assessment of genetic heterogeneity using next-generation sequencing (NGS). NGS ensures that the genetically-modified regions of the vaccine virus genome remain as designed and that levels of polymorphisms which may impact safety or efficacy are controlled during routine production. The variants present in nOPV2 lots were first assessed for temperature sensitivity and neurovirulence using molecular clones to inform which polymorphisms warranted formal evaluation during lot release. The novel use of NGS as a lot release test required formal validation of the method. Analysis of an nOPV2 lot spiked with the parental Sabin-2 strain enabled performance characteristics of the method to be assessed simultaneously at over 40 positions in the genome. These characteristics included repeatability and intermediate precision of polymorphism measurement, linearity of both spike-induced and nOPV2 lot-specific polymorphisms, and the limit-of-detection of spike-induced polymorphisms. The performance characteristics of the method met pre-defined criteria for 34 spike-induced polymorphic sites and 8 polymorphisms associated with the nOPV2 preparation; these sites collectively spanned most of the viral genome. Finally, the co-location of variants of interest on genomes was evaluated, with implications for the interpretation of NGS discussed.

Dengue virus (DENV) cycles between mosquito and mammalian hosts. To examine how DENV populations adapt to these different host environments, we used serial passage in human and mosquito cell lines and estimated fitness effects for all single-nucleotide variants in these populations using ultra-deep sequencing. This allowed us to determine the contributions of beneficial and deleterious mutations to the collective fitness of the population. Our analysis revealed that the continuous influx of a large burden of deleterious mutations counterbalances the effect of rare, host-specific beneficial mutations to shape the path of adaptation. Beneficial mutations preferentially map to intrinsically disordered domains in the viral proteome and cluster to defined regions in the genome. These phenotypically redundant adaptive alleles may facilitate host-specific DENV adaptation. Importantly, the evolutionary constraints described in our simple system mirror trends observed across DENV and Zika strains, indicating it recapitulates key biophysical and biological constraints shaping long-term viral evolution.

It remains poorly understood whether and how chaperones control protein evolution. Here the authors show how the chaperone Hsp90 shapes the sequence space of its client, poliovirus protein P1, at the polypeptide and RNA level to balance the evolutionary trade-offs between protein stability, aggregation and translation rate.

RNA viruses, such as polioviruses, have a great evolutionary capacity and can adapt quickly during infection. Here, the authors show that poliovirus infection in mice requires adaptation to innate immune microenvironments encountered in different tissues.

Aedes aegypti transmit pathogenic arboviruses while the mosquito itself tolerates the infection. We examine a piRNA-based immunity that relies on the acquisition of viral derived cDNA (vDNA) and how this pathway discriminates between self and non-self. The piRNAs derived from these vDNAs are essential for virus control and Piwi4 has a central role in the pathway. Piwi4 binds preferentially to virus-derived piRNAs but not to transposon-targeting piRNAs. Analysis of episomal vDNA from infected cells reveals that vDNA molecules are acquired through a discriminatory process of reverse-transcription and recombination directed by endogenous retrotransposons. Using a high-resolution Ae. aegypti genomic sequence, we found that vDNAs integrated in the host genome as endogenous viral elements (EVEs), produce antisense piRNAs that are preferentially loaded onto Piwi4. Importantly, EVE-derived piRNAs are specifically loaded onto Piwi4 to inhibit virus replication. Thus, Ae. aegypti employs a sophisticated antiviral mechanism that promotes viral persistence and generates long-lasting adaptive immunity.

Enterovirus (EV)-D68 has been associated with epidemics in the United Sates in 2014, 2016 and 2018. This study aims to identify potential viral virulence determinants. We found that neonatal type I interferon receptor knockout mice are susceptible to EV-D68 infection via intraperitoneal inoculation and were able to recapitulate the paralysis process observed in human disease. Among the EV-D68 strains tested, strain US/MO-14-18949 caused no observable disease in this mouse model, whereas the other strains caused paralysis and death. Sequence analysis revealed several conserved genetic changes among these virus strains: nucleotide positions 107 and 648 in the 5′-untranslated region (UTR); amino acid position 88 in VP3; 1, 148, 282 and 283 in VP1; 22 in 2A; 47 in 3A. A series of chimeric and point-mutated infectious clones were constructed to identify viral elements responsible for the distinct virulence. A single amino acid change from isoleucine to valine at position 88 in VP3 attenuated neurovirulence by reducing virus replication in the brain and spinal cord of infected mice.

The accumulation of mutations in RNA viruses is thought to facilitate rapid adaptation to changes in the environment. However, most mutations have deleterious effects on fitness, especially for viruses. Thus, tolerance to mutations should determine the nature and extent of genetic diversity that can be maintained in the population. Here, we combine population genetics theory, computer simulation, and experimental evolution to examine the advantages and disadvantages of tolerance to mutations, also known as mutational robustness. We find that mutational robustness increases neutral diversity and, as expected, can facilitate adaptation to a new environment. Surprisingly, under certain conditions, robustness may also be an impediment for viral adaptation, if a highly diverse population contains a large proportion of previously neutral mutations that are deleterious in the new environment. These findings may inform therapeutic strategies that cause extinction of otherwise robust viral populations.

Viruses interact extensively with host proteins, but the mechanisms controlling these interactions are not well understood. We present a comprehensive analysis of eukaryotic linear motifs (ELMs) in 2,208 viral genomes and reveal that viruses exploit molecular mimicry of host-like ELMs to possibly assist in host-virus interactions. Using a statistical genomics approach, we identify a large number of potentially functional ELMs and observe that the occurrence of ELMs is often evolutionarily conserved but not uniform across virus families. Some viral proteins contain multiple types of ELMs, in striking similarity to complex regulatory modules in host proteins, suggesting that ELMs may act combinatorially to assist viral replication. Furthermore, a simple evolutionary model suggests that the inherent structural simplicity of ELMs often enables them to tolerate mutations and evolve quickly. Our findings suggest that ELMs may allow fast rewiring of host-virus interactions, which likely assists rapid viral evolution and adaptation to diverse environments.

Dengue virus NS5 protein plays multiple functions in the cytoplasm of infected cells, enabling viral RNA replication and counteracting host antiviral responses. Here, we demonstrate a novel function of NS5 in the nucleus where it interferes with cellular splicing. Using global proteomic analysis of infected cells together with functional studies, we found that NS5 binds spliceosome complexes and modulates endogenous splicing as well as minigene-derived alternative splicing patterns. In particular, we show that NS5 alone, or in the context of viral infection, interacts with core components of the U5 snRNP particle, CD2BP2 and DDX23, alters the inclusion/exclusion ratio of alternative splicing events, and changes mRNA isoform abundance of known antiviral factors. Interestingly, a genome wide transcriptome analysis, using recently developed bioinformatics tools, revealed an increase of intron retention upon dengue virus infection, and viral replication was improved by silencing specific U5 components. Different mechanistic studies indicate that binding of NS5 to the spliceosome reduces the efficiency of pre-mRNA processing, independently of NS5 enzymatic activities. We propose that NS5 binding to U5 snRNP proteins hijacks the splicing machinery resulting in a less restrictive environment for viral replication.

Life history theory posits that the sequence and timing of events in an organism's lifespan are fine-tuned by evolution to maximize the production of viable offspring. In a virus, a life history strategy is largely manifested in its replication mode. Here, we develop a stochastic mathematical model to infer the replication mode shaping the structure and mutation distribution of a poliovirus population in an intact single infected cell. We measure production of RNA and poliovirus particles through the infection cycle, and use these data to infer the parameters of our model. We find that on average the viral progeny produced from each cell are approximately five generations removed from the infecting virus. Multiple generations within a single cell infection provide opportunities for significant accumulation of mutations per viral genome and for intracellular selection.

Nodamura Virus (NoV) is a nodavirus originally isolated from insects that can replicate in a wide variety of hosts, including mammals. Because of their simplicity and ability to replicate in many diverse hosts, NoV, and the Nodaviridae in general, provide a unique window into the evolution of viruses and host-virus interactions. Here we show that the C-terminus of the viral polymerase exhibits extreme structural and evolutionary flexibility. Indeed, fewer than 10 positively charged residues from the 110 amino acid-long C-terminal region of protein A are required to support RNA1 replication. Strikingly, this region can be replaced by completely unrelated protein sequences, yet still produce a functional replicase. Structure predictions, as well as evolutionary and mutational analyses, indicate that the C-terminal region is structurally disordered and evolves faster than the rest of the viral proteome. Thus, the function of an intrinsically unstructured protein region can be independent of most of its primary sequence, conferring both functional robustness and sequence plasticity on the protein. Our results provide an experimental explanation for rapid evolution of unstructured regions, which enables an effective exploration of the sequence space, and likely function space, available to the virus.

The control and prevention of communicable disease is directly impacted by the genetic mutability of the underlying etiological agents. In the case of RNA viruses, genetic recombination may impact public health by facilitating the generation of new viral strains with altered phenotypes and by compromising the genetic stability of live attenuated vaccines. The landscape of homologous recombination within a given RNA viral genome is thought to be influenced by several factors; however, a complete understanding of the genetic determinants of recombination is lacking. Here, we utilize gene synthesis and deep sequencing to create a detailed recombination map of the poliovirus 1 coding region. We identified over 50 thousand breakpoints throughout the genome, and we show the majority of breakpoints to be concentrated in a small number of specific "hotspots," including those associated with known or predicted RNA secondary structures. Nucleotide base composition was also found to be associated with recombination frequency, suggesting that recombination is modulated across the genome by predictable and alterable motifs. We tested the predictive utility of the nucleotide base composition association by generating an artificial hotspot in the poliovirus genome. Our results imply that modification of these motifs could be extended to whole genome re-designs for the development of recombination-deficient, genetically stable live vaccine strains.

Honey bees are essential pollinators of numerous agricultural crops. Since 2006, honey bee populations have suffered considerable annual losses that are partially attributed to Colony Collapse Disorder (CCD). CCD is an unexplained phenomenon that correlates with elevated incidence of pathogens, including RNA viruses. Honey bees are eusocial insects that live in colonies of genetically related individuals that work in concert to gather and store nutrients. Their social organization provides numerous benefits, but also facilitates pathogen transmission between individuals. To investigate honey bee antiviral defense mechanisms, we developed an RNA virus infection model and discovered that administration of dsRNA, regardless of sequence, reduced virus infection. Our results suggest that dsRNA, a viral pathogen associated molecular pattern (PAMP), triggers an antiviral response that controls virus infection in honey bees.

Respiratory syncytial virus (RSV) is a major cause of respiratory illness in young children, leading to significant morbidity and mortality worldwide. Despite its medical importance, no vaccine or effective therapeutic interventions are currently available. Therefore, there is a pressing need to identify novel antiviral drugs to combat RSV infections. Hsp90, a cellular protein-folding factor, has been shown to play an important role in the replication of numerous viruses. We here demonstrate that RSV requires Hsp90 for replication. Mechanistic studies reveal that inhibition of Hsp90 during RSV infection leads to the degradation of a viral protein similar in size to the RSV L protein, the viral RNA-dependent RNA polymerase, implicating it as an Hsp90 client protein. Accordingly, Hsp90 inhibitors exhibit antiviral activity against laboratory and clinical isolates of RSV in both immortalized as well as primary differentiated airway epithelial cells. Interestingly, we find a high barrier to the emergence of drug resistance to Hsp90 inhibitors, as extensive growth of RSV under conditions of Hsp90 inhibition did not yield mutants with reduced sensitivity to these drugs. Our results suggest that Hsp90 inhibitors may present attractive antiviral therapeutics for treatment of RSV infections and highlight the potential of chaperone inhibitors as antivirals exhibiting high barriers to development of drug resistance.

Honey bees (Apis mellifera) play a critical role in global food production as pollinators of numerous crops. Recently, honey bee populations in the United States, Canada, and Europe have suffered an unexplained increase in annual losses due to a phenomenon known as Colony Collapse Disorder (CCD). Epidemiological analysis of CCD is confounded by a relative dearth of bee pathogen field studies. To identify what constitutes an abnormal pathophysiological condition in a honey bee colony, it is critical to have characterized the spectrum of exogenous infectious agents in healthy hives over time. We conducted a prospective study of a large scale migratory bee keeping operation using high-frequency sampling paired with comprehensive molecular detection methods, including a custom microarray, qPCR, and ultra deep sequencing. We established seasonal incidence and abundance of known viruses, Nosema sp., Crithidia mellificae, and bacteria. Ultra deep sequence analysis further identified four novel RNA viruses, two of which were the most abundant observed components of the honey bee microbiome (∼10(11) viruses per honey bee). Our results demonstrate episodic viral incidence and distinct pathogen patterns between summer and winter time-points. Peak infection of common honey bee viruses and Nosema occurred in the summer, whereas levels of the trypanosomatid Crithidia mellificae and Lake Sinai virus 2, a novel virus, peaked in January.

A large number of medically important viruses, including HIV, hepatitis C virus, and influenza, have RNA genomes. These viruses replicate with extremely high mutation rates and exhibit significant genetic diversity. This diversity allows a viral population to rapidly adapt to dynamic environments and evolve resistance to vaccines and antiviral drugs. For the last 30 years, quasispecies theory has provided a population-based framework for understanding RNA viral evolution. A quasispecies is a cloud of diverse variants that are genetically linked through mutation, interact cooperatively on a functional level, and collectively contribute to the characteristics of the population. Many predictions of quasispecies theory run counter to traditional views of microbial behavior and evolution and have profound implications for our understanding of viral disease. Here, we discuss basic principles of quasispecies theory and describe its relevance for our understanding of viral fitness, virulence, and antiviral therapeutic strategy.

RNA structures present throughout RNA virus genomes serve as scaffolds to organize multiple factors involved in the initiation of RNA synthesis. Several of these RNA elements play multiple roles in the RNA replication pathway. An RNA structure formed around the 5'- end of the poliovirus genomic RNA has been implicated in the initiation of both negative- and positive-strand RNA synthesis. Dissecting the roles of these multifunctional elements is usually hindered by the interdependent nature of the viral replication processes and often pleiotropic effects of mutations. Here, we describe a novel approach to examine RNA elements with multiple roles. Our approach relies on the duplication of the RNA structure so that one copy is dedicated to the initiation of negative-strand RNA synthesis, while the other mediates positive-strand synthesis. This allows us to study the function of the element in promoting positive-strand RNA synthesis, independently of its function in negative-strand initiation. Using this approach, we demonstrate that the entire 5'-end RNA structure that forms on the positive-strand is required for initiation of new positive-strand RNAs. Also required to initiate positive-strand RNA synthesis are the binding sites for the viral polymerase precursor, 3CD, and the host factor, PCBP. Furthermore, we identify specific nucleotide sequences within "stem a" that are essential for the initiation of positive-strand RNA synthesis. These findings provide direct evidence for a trans-initiation model, in which binding of proteins to internal sequences of a pre-existing positive-strand RNA affects the synthesis of subsequent copies of that RNA, most likely by organizing replication factors around the initiation site.

Many techniques for the study of complex populations provide either specific information on a small number of variants or general information on the entire population. Here we describe a powerful new technique for elucidating mutation frequencies at each genomic position in a complex population. This single base extension (SBE) based microarray platform was designed and optimized using poliovirus as the target genotype, but can be easily adapted to assay populations derived from any organism. The sensitivity of the method was demonstrated by accurate and consistent readouts from a controlled population of mutant genotypes. We subsequently deployed the technique to investigate the effects of the nucleotide analog ribavirin on a typical poliovirus population through two rounds of passage. Our results show that this economical platform can be used to investigate dynamic changes occurring at frequencies below 1% within a complex nucleic acid population. Given that many key aspects of the study and treatment of disease are intimately linked to population-level genomic diversity, our SBE-based technique provides a scalable and cost-effective complement to both traditional and next generation sequencing methodologies.