Search

Integration into the nuclear genome of germ line cells can lead to vertical inheritance of retroviral genes as host alleles. For other viruses, germ line integration has only rarely been documented. Nonetheless, we identified endogenous viral elements (EVEs) derived from ten non-retroviral families by systematic in silico screening of animal genomes, including the first endogenous representatives of double-stranded RNA, reverse-transcribing DNA, and segmented RNA viruses, and the first endogenous DNA viruses in mammalian genomes. Phylogenetic and genomic analysis of EVEs across multiple host species revealed novel information about the origin and evolution of diverse virus groups. Furthermore, several of the elements identified here encode intact open reading frames or are expressed as mRNA. For one element in the primate lineage, we provide statistically robust evidence for exaptation. Our findings establish that genetic material derived from all known viral genome types and replication strategies can enter the animal germ line, greatly broadening the scope of paleovirological studies and indicating a more significant evolutionary role for gene flow from virus to animal genomes than has previously been recognized., Author Summary The presence of retrovirus sequences in animal genomes has been recognized since the 1970s, but is readily explained by the fact that these viruses integrate into chromosomal DNA as part of their normal replication cycle. Unexpectedly, however, we identified a large and diverse population of sequences in animal genomes that are derived from non-retroviral viruses. Analysis of these sequences—which represent all known virus genome types and replication strategies—reveals new information about the evolutionary history of viruses, in many cases providing the first and only direct evidence for their ancient origins. Additionally, we provide evidence that the functionality of one of these sequences has been maintained in the host genome over many millions of years, raising the possibility that captured viral sequences may have played a larger than expected role in host evolution.

Despite highly active antiretroviral therapy (HAART), AIDS related lymphoma (ARL) occurs at a significantly higher rate in patients infected with the Human Immunodeficiency Virus (HIV) than in the general population. HIV-infected macrophages are a known viral reservoir and have been shown to have lymphomagenic potential in SCID mice; therefore, there is an interest in determining if a viral component to lymphomagenesis also exists. We sequenced HIV-1 envelope gp120 clones obtained post mortem from several tumor and non-tumor tissues of two patients who died with AIDS-related Non-Hodgkin's lymphoma (ARL-NH). Similar results were found in both patients: 1) high-resolution phylogenetic analysis showed a significant degree of compartmentalization between lymphoma and non-lymphoma viral sub-populations while viral sub-populations from lymph nodes appeared to be intermixed within sequences from tumor and non-tumor tissues, 2) a 100-fold increase in the effective HIV population size in tumor versus non-tumor tissues was associated with the emergence of lymphadenopathy and aggressive metastatic ARL, and 3) HIV gene flow among lymph nodes, normal and metastatic tissues was non-random. The different population dynamics between the viruses found in tumors versus the non-tumor associated viruses suggest that there is a significant relationship between HIV evolution and lymphoma pathogenesis. Moreover, the study indicates that HIV could be used as an effective marker to study the origin and dissemination of lymphomas in vivo.

The effective size of populations (Ne) determines whether selection or genetic drift is the predominant force shaping their genetic structure and evolution. Despite their high mutation rate and rapid evolution, this parameter is poorly documented experimentally in viruses, particularly plant viruses. All available studies, however, have demonstrated the existence of huge within-host demographic fluctuations, drastically reducing Ne upon systemic invasion of different organs and tissues. Notably, extreme bottlenecks have been detected at the stage of systemic leaf colonization in all plant viral species investigated so far, sustaining the general idea that some unknown obstacle(s) imposes a barrier on the development of all plant viruses. This idea has important implications, as it appoints genetic drift as a constant major force in plant virus evolution. By co-inoculating several genetic variants of Cauliflower mosaic virus into a large number of replicate host plants, and by monitoring their relative frequency within the viral population over the course of the host systemic infection, only minute stochastic variations were detected. This allowed the estimation of the CaMV Ne during colonization of successive leaves at several hundreds of viral genomes, a value about 100-fold higher than that reported for any other plant virus investigated so far, and indicated the very limited role played by genetic drift during plant systemic infection by this virus. These results suggest that the barriers that generate bottlenecks in some plant virus species might well not exist, or can be surmounted by other viruses, implying that severe bottlenecks during host colonization do not necessarily apply to all plant-infecting viruses., Author Summary Whether selection or stochastic genetic drift is the major force driving the evolution of a virus depends largely on the size of the viral population, with the former being predominant in large populations and the latter taking over when population sizes are transiently or durably reduced. This question has been intensively debated in both plant and animal viruses, as demographic fluctuations throughout viral life cycles are poorly understood. In plant viruses, an extremely small population size—down to a few founder genome units colonizing each leaf—has been formally estimated in two instances, and all other virus species investigated so far have consistently been shown to undergo extreme demographic bottlenecks during systemic invasion of their host. This situation conveys the general idea that all viruses are confronted with “universal barriers” in plants, imposing repeated transient decreases in their population size, thus making genetic drift a major constant driver of their evolution. Here, using the example of Cauliflower mosaic virus, we mitigate this general idea by showing that at least one virus species can overcome such putative limiting barriers and massively invade leaves with hundreds to thousands of founding genome units.

Research has shown that RNA virus populations are highly variable, most likely due to low fidelity replication of RNA genomes. It is generally assumed that populations of DNA viruses will be less complex and show reduced variability when compared to RNA viruses. Here, we describe the use of high throughput sequencing for a genome wide study of viral populations from urine samples of neonates with congenital human cytomegalovirus (HCMV) infections. We show that HCMV intrahost genomic variability, both at the nucleotide and amino acid level, is comparable to many RNA viruses, including HIV. Within intrahost populations, we find evidence of selective sweeps that may have resulted from immune-mediated mechanisms. Similarly, genome wide, population genetic analyses suggest that positive selection has contributed to the divergence of the HCMV species from its most recent ancestor. These data provide evidence that HCMV, a virus with a large dsDNA genome, exists as a complex mixture of genome types in humans and offer insights into the evolution of the virus., Author Summary Human Cytomegalovirus (HCMV) is a dsDNA virus that is the leading source of birth defects associated with an infectious agent. There is currently no effective HCMV vaccine and few treatment strategies for congenital infections exist. Thus, a better understanding of HCMV infections is warranted. Limited data has shown that HCMV exists as a mixture of a few genotypes in human hosts. Here, we describe our use of high throughput sequencing to study the extent of genome wide variability within HCMV infections sampled from congenital infections. Surprisingly, we find that HCMV populations are as variable as quasispecies RNA viruses; it is commonly believed that DNA viruses are more genetically stable than RNA viruses, and thus produce homogenous populations. Additionally, we find evidence of evolutionary pressures acting on the HCMV genome, both within and among populations. These results provide the first evidence that diversity of DNA virus populations can be comparable to that of RNA virus populations.

The evolutionary benefit of viral genome segmentation is a classical, yet unsolved question in evolutionary biology and RNA genetics. Theoretical studies anticipated that replication of shorter RNA segments could provide a replicative advantage over standard size genomes. However, this question has remained elusive to experimentalists because of the lack of a proper viral model system. Here we present a study with a stable segmented bipartite RNA virus and its ancestor non-segmented counterpart, in an identical genomic nucleotide sequence context. Results of RNA replication, protein expression, competition experiments, and inactivation of infectious particles point to a non-replicative trait, the particle stability, as the main driver of fitness gain of segmented genomes. Accordingly, measurements of the volume occupation of the genome inside viral capsids indicate that packaging shorter genomes involves a relaxation of the packaging density that is energetically favourable. The empirical observations are used to design a computational model that predicts the existence of a critical multiplicity of infection for domination of segmented over standard types. Our experiments suggest that viral segmented genomes may have arisen as a molecular solution for the trade-off between genome length and particle stability. Genome segmentation allows maximizing the genetic content without the detrimental effect in stability derived from incresing genome length. © 2011 Ojosnegros et al., Supported by grants BFU2006-00863 from MEC, BFU2008-02816/BMC from MCI, 36558/06 from FIPSE, and Fundacion R. Areces. CIBERehd is funded by Instituto de Salud Carlos III.

We characterized the evolution of simian immunodeficiency virus (SIV) in the male genital tract by examining blood- and semen-associated virus from experimentally and sham vaccinated rhesus monkeys during primary infection. At the time of peak virus replication, SIV sequences were intermixed between the blood and semen supporting a scenario of high-level virus “spillover” into the male genital tract. However, at the time of virus set point, compartmentalization was apparent in 4 of 7 evaluated monkeys, likely as a consequence of restricted virus gene flow between anatomic compartments after the resolution of primary viremia. These findings suggest that SIV replication in the male genital tract evolves to compartmentalization after peak viremia resolves., Author Summary Methods to reduce the transmission of HIV-1 are hindered by a lack of information regarding early viral dynamics and evolution in the male genital tract. In the present study, we show that SIV in the blood and genital tract are homogeneous during early infection, indicating facile virus gene flow between these compartments. Importantly, the coincidence of the resolution of primary viremia with the decreased virus levels in genital secretions suggest that the dramatic fall in virus replication during early infection underlies the development of viral compartmentalization. Our demonstration of early virus compartmentalization in the male genital tract has important implications for the understanding of early events leading to infection of the male genital tract and the nature of the transmitted virus during primary retrovirus infection.

Perhaps more than any other biological discipline, the study of animal viruses is confined to the present. Virions are simply not the stuff of which robust fossils are made. Phylogenetic analysis can help by revealing deep relationships between extant viral lineages, yet such reconstructions lack detail (telling us nothing about transitional or extinct viral forms, the movement of viruses between species, or the timing of major events in viral evolution), and molecular clock estimates are notoriously imprecise when applied to viruses [1]. Until recently, ancient endogenous retroviruses (ERVs) were the closest thing to a fossil record available to scientists with a proclivity for combining virology and natural history. Happily, a trio of recent studies appearing in PLoS Genetics [2], PLoS Biology [3], and PLoS Pathogens [4] reveal an unexpected wealth of non-retroviral virus sequences embedded in the genome sequence databases, a virtual equivalent of the Burgess Shale, ripe for excavation by eager paleovirologists. Retroviral infection occasionally results in the deposition of a provirus in a host's germline DNA. While germline integration of a provirus may be an exceedingly rare event, across the great expanse of evolutionary time millions of ERV loci have accumulated in animal genomes. Because retroviruses replicate through an integrated DNA intermediate, it is not difficult to imagine how ERVs are generated. For other animal viruses, which do not normally integrate their genomes into host DNA, the formation of germline insertions should be far less likely. Nonetheless, reports of non-retroviral specimens being unearthed from the genomes of animal species are on the rise. Notable examples include functional expression of nudivirus-related structural genes in the genomes of parasitic wasps [5]; Ebolavirus-like sequences, related to modern filoviruses, present in multiple mammalian genomes [6]; and sequences resembling the Bornavirus nucleoprotein gene (N) in the genomes of various mammals including primates, rodents, and elephants [7]. Even some herpesviruses have a propensity for occasional germline insertion and thus, the potential for vertical inheritance [8]. Now, Belyi et al. [4] and Katzourakis and Gifford [2], have unearthed diverse collections of non-retroviral sequences buried in whole genome sequence data from an impressive array of host organisms, including mammals, marsupials, birds, rodents, and insects, using modern viral sequences as bioinformatic probes. A third study from Gilbert and Feschotte specifically reevaluates the macroevolution of hepadnaviruses based on the sequence and distribution of hepadnavirus-like fossils in the genomes of passerine birds [3]. To cope with this newfound abundance, the authors of one of the studies suggest the acronym EVE (for endogenous viral element) as a general term to encompass all virus-derived genomic loci [2]. Two of the studies also took a closer look at a previously described class of EVEs, called EBLNs (for endogenous Bornavirus-like N genes) [2], [4], [7]. While most EVEs were either defective at the time of insertion or rendered functionless by the accumulation of random mutations over the course of millions of years, EBLNs are striking in retaining largely intact protein-coding sequences. In fact, in silico simulations of EBLN evolution estimate that these elements should have accumulated ∼10–20 stop codons since the time of genome insertion. That the EBLN coding sequences appear relatively unscathed suggests that these particular elements provide (or at times provided) a selectively advantageous function, subjecting them to purifying selection. The possibility is not without precedent: for example, at least one human ERV has evolved to provide a cellular function [9], and there are several examples of ERVs that have been subverted by host evolution to serve as inhibitors of retroviral infection [10]–[14]. As a group, viruses are polyphyletic, as evidenced by the variety of unique genome types and distinctive replication strategies they collectively employ. There are double-stranded DNA viruses and single-stranded DNA viruses, double-stranded and single-stranded RNA viruses, and viruses with segmented genomes; among those with single-stranded RNA, there are those with positive polarity (the genome resembles an mRNA) and those with negative sense genomes. Each genome type represents a different starting point for takeover of the host cell, and each requires a different strategy for achieving this fundamental task. For example, replication of some viruses is confined entirely to the cytoplasm, whereas others involve synthesis of DNA or RNA in the nucleus. While the fossil record is still dominated by retroviral sequences, the inventory of known EVE loci now appears to include representatives of all the basic replication strategies exemplified by modern viruses. Non-retroviral EVEs are typically subgenomic, derived from just one or a few viral genes instead of entire viral genomes. Insertion site duplications bracketing some EVEs suggest that retrotransposition in trans, by retrotransposons or possibly retroviruses, may be a predominant mechanism of EVE formation. In fact, for RNA viruses that replicate in the cytoplasm (e.g., filoviruses and rhabdoviruses), retrotransposition is the most plausible mechanism for EVE formation. In such cases, it will be interesting to determine whether the more abundant EVE sequences share some common feature(s) conferring a propensity for retrotransposition. In contrast, hepadnavirus “fossils” lack the hallmarks of retrotransposition (such as flanking insertion-site duplications and poly-A tails), and may instead have resulted from non-homologous end joining and insertion of viral DNA directly into the host genome [3]. When incorporated into phylogenetic trees, many EVEs group as sister taxa to their modern counterparts. Thus, they are not evolutionary intermediates on the path to extant viruses, but rather extinct lineages sharing a common ancestor with modern viruses. From this, one can infer that most of the distinctive replication strategies employed by modern viruses probably originated hundreds of millions of years ago. While virologists intuitively understand this (given the widespread distribution of viruses among living organisms), EVEs constitute direct, physical evidence that modern viral lineages have very ancient roots (Table 1). That modern viruses and ancient EVE sequences are still recognizably related is astonishing, given that they are separated by millions of years of exogenous viral evolution. Table 1 Estimated Minimum Age of Select EVEs. The catalog of EVEs is impressive for what it contains, but even more so for what it does not. Why? Because the known EVEs probably represent a minor and highly skewed sampling of viral prehistory. Minor, because the odds that infection of an individual organism will result in fixation of an EVE are exceedingly small (Figure 1). Skewed, because some viruses may be more prone to germline insertion than others (that retroviral insertions greatly outnumber other EVEs is a particularly striking example of a virus-dependent bias). Thus, as impressive in scope and variety as the EVEs are, they may represent but a drop in the ocean of all the viruses that have buffeted host organisms across the ages. Figure 1 Formation of a hypothetical EVE and relationship to modern viruses. The current EVE record may have other limitations. Just how far back does the EVE fossil record extend? Erosion due to the steady accumulation of mutations must impose an upper limit on how far back the viral fossil record can be deciphered, and theoretical predictions of that limit would be useful. Even in the absence of sequence degradation, some EVEs may be easier to detect than others. For example, the studies described here relied on known viral sequences as queries: if our genomes also harbor ancient viral sequences for which there is no modern counterpart, how would we recognize them for what they once were?

Although overall pollinator populations have declined over the last couple of decades, the honey bee (Apis mellifera) malady, colony collapse disorder (CCD), has caused major concern in the agricultural community. Among honey bee pathogens, RNA viruses are emerging as a serious threat and are suspected as major contributors to CCD. Recent detection of these viral species in bumble bees suggests a possible wider environmental spread of these viruses with potential broader impact. It is therefore vital to study the ecology and epidemiology of these viruses in the hymenopteran pollinator community as a whole. We studied the viral distribution in honey bees, in their pollen loads, and in other non-Apis hymenopteran pollinators collected from flowering plants in Pennsylvania, New York, and Illinois in the United States. Viruses in the samples were detected using reverse transcriptase-PCR and confirmed by sequencing. For the first time, we report the molecular detection of picorna-like RNA viruses (deformed wing virus, sacbrood virus and black queen cell virus) in pollen pellets collected directly from forager bees. Pollen pellets from several uninfected forager bees were detected with virus, indicating that pollen itself may harbor viruses. The viruses in the pollen and honey stored in the hive were demonstrated to be infective, with the queen becoming infected and laying infected eggs after these virus-contaminated foods were given to virus-free colonies. These viruses were detected in eleven other non-Apis hymenopteran species, ranging from many solitary bees to bumble bees and wasps. This finding further expands the viral host range and implies a possible deeper impact on the health of our ecosystem. Phylogenetic analyses support that these viruses are disseminating freely among the pollinators via the flower pollen itself. Notably, in cases where honey bee apiaries affected by CCD harbored honey bees with Israeli Acute Paralysis virus (IAPV), nearby non-Apis hymenopteran pollinators also had IAPV, while those near apiaries without IAPV did not. In containment greenhouse experiments, IAPV moved from infected honey bees to bumble bees and from infected bumble bees to honey bees within a week, demonstrating that the viruses could be transmitted from one species to another. This study adds to our present understanding of virus epidemiology and may help explain bee disease patterns and pollinator population decline in general.

BackgroundHepatitis E virus (HEV) is an enterically transmitted hepatropic virus. It segregates as four genotypes. All genotypes infect humans while only genotypes 3 and 4 also infect several animal species. It has been suggested that hepatitis E is zoonotic, but no study has analyzed the evolutionary history of HEV. We present here an analysis of the evolutionary history of HEV.Methods and findingsThe times to the most recent common ancestors for all four genotypes of HEV were calculated using BEAST to conduct a Bayesian analysis of HEV. The population dynamics for genotypes 1, 3 and 4 were analyzed using skyline plots. Bayesian analysis showed that the most recent common ancestor for modern HEV existed between 536 and 1344 years ago. The progenitor of HEV appears to have given rise to anthropotropic and enzootic forms of HEV, which evolved into genotypes 1 and 2 and genotypes 3 and 4, respectively. Population dynamics suggest that genotypes 1, 3 and 4 experienced a population expansion during the 20(th) century. Genotype 1 has increased in infected population size ∼30-35 years ago. Genotype 3 and 4 have experienced an increase in population size starting late in the 19(th) century until ca.1940-45, with genotype 3 having undergone additional rapid expansion until ca.1960. The effective population size for both genotype 3 and 4 rapidly declined to pre-expansion levels starting in ca.1990. Genotype 4 was further examined as Chinese and Japanese sequences, which exhibited different population dynamics, suggesting that this genotype experienced different evolutionary history in these two countries.ConclusionsHEV appears to have evolved through a series of steps, in which the ancestors of HEV may have adapted to a succession of animal hosts leading to humans. Analysis of the population dynamics of HEV suggests a substantial temporal variation in the rate of transmission among HEV genotypes in different geographic regions late in the 20(th) Century.

Antigenic drift in the influenza A virus hemagglutinin (HA) is responsible for seasonal reformulation of influenza vaccines. Here, we address an important and largely overlooked issue in antigenic drift: how does the number and location of glycosylation sites affect HA evolution in man? We analyzed the glycosylation status of all full-length H1 subtype HA sequences available in the NCBI influenza database. We devised the “flow index” (FI), a simple algorithm that calculates the tendency for viruses to gain or lose consensus glycosylation sites. The FI predicts the predominance of glycosylation states among existing strains. Our analyses show that while the number of glycosylation sites in the HA globular domain does not influence the overall magnitude of variation in defined antigenic regions, variation focuses on those regions unshielded by glycosylation. This supports the conclusion that glycosylation generally shields HA from antibody-mediated neutralization, and implies that fitness costs in accommodating oligosaccharides limit virus escape via HA hyperglycosylation., Author Summary Influenza A virus is highly susceptible to neutralizing antibodies specific for the viral hemagglutinin glycoprotein (HA), and is easily controlled by standard vaccines. Influenza A virus remains an important human pathogen, however, due to its ability to rapidly evade antibody responses. This process, termed antigenic drift, is due to the accumulation of amino acid substitutions that modify HA antigenic sites recognized by neutralizing antibodies. In this study, we perform bioinformatic analysis on thousands of influenza A virus isolates to better understand the influence of N-linked glycosylation on antigenic drift. HA from human IAV isolates can accommodate up to 6 oligosaccharides in its globular domain. We show that for H1, H2, and to a somewhat less extent H3, HAs, the number of glycosylation sites in the globular domain does not greatly modify the total degree of variation in antigenic sites, but rather focuses variation on sites whose access to antibodies is unaffected by glycosylation. Our findings imply that glycosylation protects HA from antibody neutralization, but functional impairment limits the number of oligosaccharides that HA can accommodate.

Viruses are masters of evolution due to high frequency mutations and genetic recombination. In spite of the significance of viral RNA recombination that promotes the emergence of drug-resistant virus strains, the role of host and environmental factors in RNA recombination is poorly understood. Here we report that the host Met22p/Hal2p bisphosphate-3′-nucleotidase regulates the frequency of viral RNA recombination and the efficiency of viral replication. Based on Tomato bushy stunt virus (TBSV) and yeast as a model host, we demonstrate that deletion of MET22 in yeast or knockdown of AHL, SAL1 and FRY1 nucleotidases/phosphatases in plants leads to increased TBSV recombination and replication. Using a cell-free TBSV recombination/replication assay, we show that the substrate of the above nucleotidases, namely 3′-phosphoadenosine-5′-phosphate pAp, inhibits the activity of the Xrn1p 5′-3′ ribonuclease, a known suppressor of TBSV recombination. Inhibition of the activity of the nucleotidases by LiCl and NaCl also leads to increased TBSV recombination, demonstrating that environmental factors could also affect viral RNA recombination. Thus, host factors in combination with environmental factors likely affect virus evolution and adaptation., Author Summary Viral RNA recombination plays a major role in virus evolution. Yet, we know little about the roles of host and environmental factors in viral RNA recombination. In this work, using TBSV and yeast as a model host, we show that MET22 nucleotidase suppresses viral RNA recombination. In vitro experiments with a cell-free extract from yeast revealed that the substrate of Met22p bisphosphate-3′-nucleotidase, namely pAp, could promote TBSV RNA recombination via inhibiting the activity of the Xrn1p 5′-3′ ribonuclease, a known suppressor of viral RNA recombination. Altogether, we demonstrate that the Met22/Xrn1 pathway and environmental factors affecting this pathway, namely salt-stress caused by LiCl and NaCl, play a role in viral RNA recombination. The authors also provide evidence that a similar pathway affects TBSV recombination in plants as well. The most pronounced increase in accumulation of TBSV RNA recombinants was seen in the triple-nucleotidase gene silenced plants treated with LiCl, suggesting that the combined effect of genetic and environmental factors could be critical in regulation of the rate of viral RNA recombination.

Mosquitoes rely on RNA interference (RNAi) as their primary defense against viral infections. To this end, the combination of RNAi and invertebrate cell culture systems has become an invaluable tool in studying virus-vector interactions. Nevertheless, a recent study failed to detect an active RNAi response to West Nile virus (WNV) infection in C6/36 (Aedes albopictus) cells, a mosquito cell line frequently used to study arthropod-borne viruses (arboviruses). Therefore, we sought to determine if WNV actively evades the host's RNAi response or if C6/36 cells have a dysfunctional RNAi pathway. C6/36 and Drosophila melanogaster S2 cells were infected with WNV (Flaviviridae), Sindbis virus (SINV, Togaviridae) and La Crosse virus (LACV, Bunyaviridae) and total RNA recovered from cell lysates. Small RNA (sRNA) libraries were constructed and subjected to high-throughput sequencing. In S2 cells, virus-derived small interfering RNAs (viRNAs) from all three viruses were predominantly 21 nt in length, a hallmark of the RNAi pathway. However, in C6/36 cells, viRNAs were primarily 17 nt in length from WNV infected cells and 26–27 nt in length in SINV and LACV infected cells. Furthermore, the origin (positive or negative viral strand) and distribution (position along viral genome) of S2 cell generated viRNA populations was consistent with previously published studies, but the profile of sRNAs isolated from C6/36 cells was altered. In total, these results suggest that C6/36 cells lack a functional antiviral RNAi response. These findings are analogous to the type-I interferon deficiency described in Vero (African green monkey kidney) cells and suggest that C6/36 cells may fail to accurately model mosquito-arbovirus interactions at the molecular level., Author Summary Cell culture systems are invaluable tools for studying virus-host interactions. These systems are typically easy to maintain and manipulate; however, they can fail to accurately mimic the host environment encountered by viruses. Therefore, defining the limitations of each system is critical to properly interpreting the results. C6/36 Aedes albopictus cells are commonly used to study arthropod-borne viruses (arboviruses), such as West Nile virus (WNV). Recent evidence suggests that the RNA interference (RNAi) pathway, a critical aspect of the cellular innate antiviral immune response in invertebrates, may not actively target WNV in C6/36 cells. However, it is unknown whether this observation is limited to WNV. Therefore, we examined small RNA populations from C6/36 and Drosophila melanogastor S2 cells infected with WNV, Sindbis virus and La Crosse virus by high-throughput sequencing. We demonstrate that the RNAi pathway actively targets each of the three viruses in S2 cells, but does not in C6/36 cells. These findings suggest that C6/36 cells may fail to accurately model mosquito-arbovirus interactions.

Background Predator-prey models for virus-host interactions predict that viruses will cause oscillations of microbial host densities due to an arms race between resistance and virulence. A new form of microbial resistance, CRISPRs (clustered regularly interspaced short palindromic repeats) are a rapidly evolving, sequence-specific immunity mechanism in which a short piece of invading viral DNA is inserted into the host's chromosome, thereby rendering the host resistant to further infection. Few studies have linked this form of resistance to population dynamics in natural microbial populations. Methodology/Principal Findings We examined sequence diversity in 39 strains of the archeaon Sulfolobus islandicus from a single, isolated hot spring from Kamchatka, Russia to determine the effects of CRISPR immunity on microbial population dynamics. First, multiple housekeeping genetic markers identify a large clonal group of identical genotypes coexisting with a diverse set of rare genotypes. Second, the sequence-specific CRISPR spacer arrays split the large group of isolates into two very different groups and reveal extensive diversity and no evidence for dominance of a single clone within the population. Conclusions/Significance The evenness of resistance genotypes found within this population of S. islandicus is indicative of a lack of strain dominance, in contrast to the prediction for a resistant strain in a simple predator-prey interaction. Based on evidence for the independent acquisition of resistant sequences, we hypothesize that CRISPR mediated clonal interference between resistant strains promotes and maintains diversity in this natural population.

The high-risk HPV E6 and E7 proteins cooperate to immortalize primary human cervical cells and the E7 protein can independently transform fibroblasts in vitro, primarily due to its ability to associate with and degrade the retinoblastoma tumor suppressor protein, pRb. The binding of E7 to pRb is mediated by a conserved Leu-X-Cys-X-Glu (LXCXE) motif in the conserved region 2 (CR2) of E7 and this domain is both necessary and sufficient for E7/pRb association. In the current study, we report that the E7 protein of the malignancy-associated canine papillomavirus type 2 encodes an E7 protein that has serine substituted for cysteine in the LXCXE motif. In HPV, this substitution in E7 abrogates pRb binding and degradation. However, despite variation at this critical site, the canine papillomavirus E7 protein still bound and degraded pRb. Even complete deletion of the LXSXE domain of canine E7 failed to interfere with binding to pRb in vitro and in vivo. Rather, the dominant binding site for pRb mapped to the C-terminal domain of canine E7. Finally, while the CR1 and CR2 domains of HPV E7 are sufficient for degradation of pRb, the C-terminal region of canine E7 was also required for pRb degradation. Screening of HPV genome sequences revealed that the LXSXE motif of the canine E7 protein was also present in the gamma HPVs and we demonstrate that the gamma HPV-4 E7 protein also binds pRb in a similar way. It appears, therefore, that the type 2 canine PV and gamma-type HPVs not only share similar properties with respect to tissue specificity and association with immunosuppression, but also the mechanism by which their E7 proteins interact with pRb., Author Summary Human papillomaviruses (HPVs) are estimated to cause the most common sexually transmitted infection in the world, and these infections are recognized as the major cause of cervical cancer. One of the papillomavirus oncoproteins, E7, plays a major role in both the viral life cycle and progression to cancer. In cells E7 associates and inactivates pRb, a tumor suppressor protein. For the vast majority of papillomaviruses, E7 binds to pRb using a small amino acid sequence, LXCXE. However, we have now identified a papillomavirus E7 protein that lacks the LXCXE domain yet still binds and degrades pRb. This E7 protein, derived from a carcinogenic canine virus, uses its C-terminal domain to bind pRb. In addition, we discovered that a family of papillomaviruses, the gamma type HPVs, also lacks the LXCXE domain and binds pRb using a similar mechanism.

Background Despite high potential for HIV-1 genetic variation, the emergence of some mutations is constrained by fitness costs, and may be associated with compensatory amino acid (AA) co-variation. To characterize the interplay between Cytotoxic T Lymphocyte (CTL)-mediated pressure and HIV-1 evolutionary pathways, we investigated AA co-variation in Gag sequences obtained from 449 South African individuals chronically infected with HIV-1 subtype C. Methodology/Principal Findings Individuals with CTL responses biased toward Gag presented lower viral loads than individuals with under-represented Gag-specific CTL responses. Using methods that account for founder effects and HLA linkage disequilibrium, we identified 35 AA sites under Human Leukocyte Antigen (HLA)-restricted CTL selection pressure and 534 AA-to-AA interactions. Analysis of two-dimensional distances between co-varying residues revealed local stabilization mechanisms since 40% of associations involved neighboring residues. Key features of our co-variation analysis included sites with a high number of co-varying partners, such as HLA-associated sites, which had on average 55% more connections than other co-varying sites. Conclusions/Significance Clusters of co-varying AA around HLA-associated sites (especially at typically conserved sites) suggested that cooperative interactions act to preserve the local structural stability and protein function when CTL escape mutations occur. These results expose HLA-imprinted HIV-1 polymorphisms and their interlinked mutational paths in Gag that are likely due to opposite selective pressures from host CTL-mediated responses and viral fitness constraints., Bill & Melinda Gates Foundation (#43437), American Foundation for AIDS Research (#107005-43-RFNT), Center for AIDS Research (CFAR), United States. Public Health Service (AI057005)

A novel parvovirus, provisionally named Gorilla Bocavirus species 1 (GBoV1), was identified in four stool samples from Western gorillas (Gorilla gorilla) with acute enteritis. The complete genomic sequence of the new parvovirus revealed three open reading frames (ORFs) with an organization similar to that of known bocaviruses. Phylogenetic analysis using complete capsid and non structural (NS) gene sequence suggested that the new parvovirus is most closely related to human bocaviruses (HBoV). However, the NS ORF is more similar in length to the NS ORF found in canine minute virus and bovine parvovirus than in HBoV. Comparative genetic analysis using GBoV and HBoV genomes enabled characterization of unique splice donor and acceptor sites that appear to be highly conserved among all four HBoV species, and provided evidence for expression of two different NS proteins in all primate bocaviruses. GBoV is the first non-human primate bocavirus identified and provides new insights into the genetic diversity and evolution of this highly prevalent and recently discovered group of parvoviruses.

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.

BackgroundThe HIV-1 subtype B epidemic in Brazil is peculiar because of the high frequency of isolates having the GWGR tetramer at V3 loop region. It has been suggested that GWGR is a distinct variant and less pathogenic than other subtype B isolates.Methodology/principal findingsNinety-four percent of the HIV-1 subtype B worldwide sequences (7689/8131) obtained from the Los Alamos HIV database contain proline at the tetramer of the V3 loop of the env gene (GPGR) and only 0.74% (60/8131) have tryptophan (GWGR). By contrast, 48.4% (161/333) of subtype B isolates from Brazil have proline, 30.6% (102/333) contain tryptophan and 10.5% (35/333) have phenylalanine (F) at the second position of the V3 loop tip. The proportion of tryptophan and phenylalanine in Brazilian isolates is much higher than in worldwide subtype B sequences (chi-square test, p = 0.0001). The combined proportion of proline, tryptophan and phenylalanine (GPGR+GWGR+GFGR) of Brazilian isolates corresponds to 89% of all amino acids in the V3 loop. Phylogenetic analysis revealed that almost all subtype B isolates in Brazil have a common origin regardless of their motif (GWGR, GPGR, GGGR, etc.) at the V3 tetramer. This shared ancestral origin was also observed in CRF28_BF and CRF29_BF in a genome region (free of recombination) derived from parental subtype B. These results imply that tryptophan substitution (e.g., GWGR-to-GxGR), which was previously associated with the change in the coreceptor usage within the host, also occurs at the population level.Conclusions/significanceBased on the current findings and previous study showing that tryptophan and phenylalanine in the V3 loop are related with coreceptor usage, we propose that tryptophan and phenylalanine in subtype B isolates in Brazil are kept by selective mechanisms due to the distinct coreceptor preferences in target cells of GWGR, GFGR and GFGR viruses.

Background: Dengue is currently the most important arthropod-borne viral disease of humans. Recent work has shown dengue virus displays limited replication in its primary vector, the mosquito Aedes aegypti, when the insect harbors the endosymbiotic bacterium Wolbachia pipientis. Wolbachia-mediated inhibition of virus replication may lead to novel methods of arboviral control, yet the functional and cellular mechanisms that underpin it are unknown. Methodology/Principal Findings: Using paired Wolbachia-infected and uninfected Aedes-derived cell lines and dengue virus, we confirm the phenomenon of viral inhibition at the cellular level. Although Wolbachia imposes a fitness cost to cells via reduced proliferation, it also provides a significant degree of protection from virus-induced mortality. The extent of viral inhibition is related to the density of Wolbachia per cell, with highly infected cell lines showing almost complete protection from dengue infection and dramatically reduced virus titers compared to lines not infected with the bacteria. Conclusions/Significance: We have shown that cells infected with Wolbachia display inhibition of dengue virus replication, that the extent of inhibition is related to bacterial density and that Wolbachia infection, although costly, will provide a fitness benefit in some circumstances. Our results parallel findings in mosquitoes and flies, indicating that cell line models will provide useful and experimentally tractable models to study the mechanisms underlying Wolbachia-mediated protection from viruses.

Background The human leukocyte antigen (HLA)-restricted cytotoxic T-lymphocyte (CTL) immune response is one of the major factors determining the genetic diversity of human immunodeficiency virus (HIV). There are few population-based analyses of the amino acid variations associated with the host HLA type and their clinical relevance for the Asian population. Here, we identified HLA-associated polymorphisms in the HIV-1 CRF01_AE Gag protein in infected married couples, and examined the consequences of these HLA-selected mutations after transmission to HLA-unmatched recipients. Methodology/Principal Findings One hundred sixteen HIV-1-infected couples were recruited at a government hospital in northern Thailand. The 1.7-kb gag gene was amplified and directly sequenced. We identified 56 associations between amino acid variations in Gag and HLA alleles. Of those amino acid variations, 35 (62.5%) were located within or adjacent to regions reported to be HIV-specific CTL epitopes restricted by the relevant HLA. Interestingly, a significant number of HLA-associated amino acid variations appear to be unique to the CRF01_AE-infected Thai population. Variations in the capsid protein (p24) had the strongest associations with the viral load and CD4 cell count. The mutation and reversion rates after transmission to a host with a different HLA environment varied considerably. The p24 T242N variant escape from B57/58 CTL had a significant impact on the HIV-1 viral load of CRF01_AE-infected patients. Conclusions/Significance HLA-associated amino acid mutations and the CTL selection pressures on the p24 antigen appear to have the most significant impact on HIV replication in a CRF01_AE-infected Asian population. HLA-associated mutations with a low reversion rate accumulated as a footprint in this Thai population. The novel HLA-associated mutations identified in this study encourage us to acquire more extensive information about the viral dynamics of HLA-associated amino acid polymorphisms in a given population as effective CTL vaccine targets.

Background In March 2009, pandemic influenza A(H1N1) (A(H1N1)pdm) emerged in Mexico and the United States. In Japan, since the first outbreak of A(H1N1)pdm in Osaka and Hyogo Prefectures occurred in the middle of May 2009, the virus had spread over 16 of 47 prefectures as of June 4, 2009. Methods/Principal Findings We analyzed all-segment concatenated genome sequences of 75 isolates of A(H1N1)pdm viruses in Japan, and compared them with 163 full-genome sequences in the world. Two analyzing methods, distance-based and Bayesian coalescent MCMC inferences were adopted to elucidate an evolutionary relationship of the viruses in the world and Japan. Regardless of the method, the viruses in the world were classified into four distinct clusters with a few exceptions. Cluster 1 was originated earlier than cluster 2, while cluster 2 was more widely spread around the world. The other two clusters (clusters 1.2 and 1.3) were suggested to be distinct reassortants with different types of segment assortments. The viruses in Japan seemed to be a multiple origin, which were derived from approximately 28 transported cases. Twelve cases were associated with monophyletic groups consisting of Japanese viruses, which were referred to as micro-clade. While most of the micro-clades belonged to the cluster 2, the clade of the first cases of infection in Japan originated from cluster 1.2. Micro-clades of Osaka/Kobe and the Fukuoka cases, both of which were school-wide outbreaks, were eradicated. Time of most recent common ancestor (tMRCA) for each micro-clade demonstrated that some distinct viruses were transmitted in Japan between late May and early June, 2009, and appeared to spread nation-wide throughout summer. Conclusions Our results suggest that many viruses were transmitted from abroad in late May 2009 irrespective of preventive actions against the pandemic influenza, and that the influenza A(H1N1)pdm had become a pandemic stage in June 2009 in Japan.

The interferon-induced BST-2 protein has the unique ability to restrict the egress of HIV-1, Kaposi's sarcoma–associated herpesvirus (KSHV), Ebola virus, and other enveloped viruses. The observation that virions remain attached to the surface of BST-2-expressing cells led to the renaming of BST-2 as “tetherin”. However, viral proteins such as HIV-1 Vpu, simian immunodeficiency virus Nef, and KSHV K5 counteract BST-2, thereby allowing mature virions to readily escape from infected cells. Since the anti-viral function of BST-2 was discovered, there has been an explosion of research into several aspects of this intriguing interplay between host and virus. This review focuses on recent work addressing the molecular mechanisms involved in BST-2 restriction of viral egress and the species-specific countermeasures employed by various viruses.

Elucidating virus-host interactions responsible for HIV-1 transmission is important for advancing HIV-1 prevention strategies. To this end, single genome amplification (SGA) and sequencing of HIV-1 within the context of a model of random virus evolution has made possible for the first time an unambiguous identification of transmitted/founder viruses and a precise estimation of their numbers. Here, we applied this approach to HIV-1 env analyses in a cohort of acutely infected men who have sex with men (MSM) and found that a high proportion (10 of 28; 36%) had been productively infected by more than one virus. In subjects with multivariant transmission, the minimum number of transmitted viruses ranged from 2 to 10 with viral recombination leading to rapid and extensive genetic shuffling among virus lineages. A combined analysis of these results, together with recently published findings based on identical SGA methods in largely heterosexual (HSX) cohorts, revealed a significantly higher frequency of multivariant transmission in MSM than in HSX [19 of 50 subjects (38%) versus 34 of 175 subjects (19%); Fisher's exact p = 0.008]. To further evaluate the SGA strategy for identifying transmitted/founder viruses, we analyzed 239 overlapping 5′ and 3′ half genome or env-only sequences from plasma viral RNA (vRNA) and blood mononuclear cell DNA in an MSM subject who had a particularly well-documented virus exposure history 3–6 days before symptom onset and 14–17 days before peak plasma viremia (47,600,000 vRNA molecules/ml). All 239 sequences coalesced to a single transmitted/founder virus genome in a time frame consistent with the clinical history, and a molecular clone of this genome encoded replication competent virus in accord with model predictions. Higher multiplicity of HIV-1 infection in MSM compared with HSX is consistent with the demonstrably higher epidemiological risk of virus acquisition in MSM and could indicate a greater challenge for HIV-1 vaccines than previously recognized., Author Summary Understanding the biology of sexual transmission of HIV-1 could contribute importantly to the development of effective prevention measures. However, different routes of virus transmission (vaginal, rectal, penile or oral) and inaccessibility of tissues at or near the time of virus transmission make this goal elusive. Here, we apply single genome amplification and sequencing of plasma HIV-1 and a model of random virus evolution to a cohort of acutely infected men who have sex with men (MSM) and find that MSM are twice as likely as heterosexuals to become infected by multiple viruses as opposed to a single virus. Some MSM subjects were infected by as many as 7 to 10 or more genetically distinct viruses as a consequence of a single exposure event. We go on to molecularly clone the first full-length transmitted/founder subtype B HIV-1 virus and show that it is highly replicative in human CD4+ T-cells but not macrophages. Our study provides the first comparative, quantitative analysis of the multiplicity of HIV-1 infection in the two primary risk groups—MSM and heterosexuals—driving the global pandemic, and we discuss the implications of the findings to HIV-1 vaccine development and prevention research.

Many thousands of endoparasitic wasp species are known to inject polydnavirus (PDV) particles into their caterpillar host during oviposition, causing immune and developmental dysfunctions that benefit the wasp larva. PDVs associated with braconid and ichneumonid wasps, bracoviruses and ichnoviruses respectively, both deliver multiple circular dsDNA molecules to the caterpillar. These molecules contain virulence genes but lack core genes typically involved in particle production. This is not completely unexpected given that no PDV replication takes place in the caterpillar. Particle production is confined to the wasp ovary where viral DNAs are generated from proviral copies maintained within the wasp genome. We recently showed that the genes involved in bracovirus particle production reside within the wasp genome and are related to nudiviruses. In the present work we characterized genes involved in ichnovirus particle production by analyzing the components of purified Hyposoter didymator Ichnovirus particles by LC-MS/MS and studying their organization in the wasp genome. Their products are conserved among ichnovirus-associated wasps and constitute a specific set of proteins in the virosphere. Strikingly, these genes are clustered in specialized regions of the wasp genome which are amplified along with proviral DNA during virus particle replication, but are not packaged in the particles. Clearly our results show that ichnoviruses and bracoviruses particles originated from different viral entities, thus providing an example of convergent evolution where two groups of wasps have independently domesticated viruses to deliver genes into their hosts., Author Summary The polydnaviruses (PDVs) are a unique virus type used by an organism (a parasitic wasp) to manipulate the physiology of another organism (a lepidopteran host) in order to ensure successful parasitism. The evolutionary origin of these unusual viruses, found in ∼17,500 braconid wasps (Bracoviruses) and ∼15,000 ichneumonid wasps (Ichnoviruses), has been a major question for the last decade. We thus undertook an exclusive work aiming at investigating this origin via the characterization of genes encoding structural components for both types of PDVs. The present paper constitutes the first report on the identity and genome organisation of the viral machinery producing Ichnovirus virions. Our results strongly suggest that Ichnoviruses originated from a virus belonging to a group as yet uncharacterized that integrated its genome into that of an ichneumonid wasp ancestor. More importantly, our results demonstrate that the ancestor of Ichnoviruses differs from that of Bracoviruses, which originated from a nudivirus. We have now identified, for the two types of PDVs, the non packaged viral genes and their products involved in producing particles injected into the host during oviposition. Together, these data provide an example of convergent evolution where different groups of wasps have independently domesticated viruses to deliver genes into their hosts.

BACKGROUND: In Nigeria, hepatitis B virus (HBV) infection has reached hyperendemic levels and its nature and origin have been described as a puzzle. In this study, we investigated the molecular epidemiology and epidemic history of HBV infection in two semi-isolated rural communities in North/Central Nigeria. It was expected that only a few, if any, HBV strains could have been introduced and effectively transmitted among these residents, reflecting limited contacts of these communities with the general population in the country. METHODS AND FINDINGS: Despite remoteness and isolation, approximately 11% of the entire population in these communities was HBV-DNA seropositive. Analyses of the S-gene sequences obtained from 55 HBV-seropositive individuals showed the circulation of 37 distinct HBV variants. These HBV isolates belong predominantly to genotype E (HBV/E) (n=53, 96.4%), with only 2 classified as sub-genotype A3 (HBV/A3). Phylogenetic analysis showed extensive intermixing between HBV/E variants identified in these communities and different countries in Africa. Quasispecies analysis of 22 HBV/E strains using end-point limiting-dilution real-time PCR, sequencing and median joining networks showed extensive intra-host heterogeneity and inter-host variant sharing. To investigate events that resulted in such remarkable HBV/E diversity, HBV full-size genome sequences were obtained from 47 HBV/E infected persons and P gene was subjected to Bayesian coalescent analysis. The time to the most recent common ancestor (tMRCA) for these HBV/E variants was estimated to be year 1952 (95% highest posterior density (95% HPD): 1927-1970). Using additional HBV/E sequences from other African countries, the tMRCA was estimated to be year 1948 (95% HPD: 1924-1966), indicating that HBV/E in these remote communities has a similar time of origin with multiple HBV/E variants broadly circulating in West/Central Africa. Phylogenetic analysis and statistical neutrality tests suggested rapid HBV/E population expansion. Additionally, skyline plot analysis showed an increase in the size of the HBV/E-infected population over the last approximately 30-40 years. CONCLUSIONS: Our data suggest a massive introduction and relatively recent HBV/E expansion in the human population in Africa. Collectively, these data show a significant shift in the HBV/E epidemic dynamics in Africa over the last century.

Newcastle Disease Virus (NDV) is a pathogenic strain of avian paramyxovirus (aPMV-1) that is among the most serious of disease threats to the poultry industry worldwide. Viral diversity is high in aPMV-1; eight genotypes are recognized based on phylogenetic reconstruction of gene sequences. Modified live vaccines have been developed to decrease the economic losses caused by this virus. Vaccines derived from avirulent genotype II strains were developed in the 1950s and are in use globally, whereas Australian strains belonging to genotype I were developed as vaccines in the 1970s and are used mainly in Asia. In this study, we evaluated the consequences of attenuated live virus vaccination on the evolution of aPMV-1 genotypes. There was phylogenetic incongruence among trees based on individual genes and complete coding region of 54 full length aPMV-1 genomes, suggesting that recombinant sequences were present in the data set. Subsequently, five recombinant genomes were identified, four of which contained sequences from either genotype I or II. The population history of vaccine-related genotype II strains was distinct from other aPMV-1 genotypes; genotype II emerged in the late 19th century and is evolving more slowly than other genotypes, which emerged in the 1960s. Despite vaccination efforts, genotype II viruses have experienced constant population growth to the present. In contrast, other contemporary genotypes showed population declines in the late 1990s. Additionally, genotype I and II viruses, which are circulating in the presence of homotypic vaccine pressure, have unique selection profiles compared to nonvaccine-related strains. Collectively, these data show that vaccination with live attenuated viruses has changed the evolution of aPMV-1 by maintaining a large effective population size of a vaccine-related genotype, allowing for coinfection and recombination of vaccine and wild type strains, and by applying unique selective pressures on viral glycoproteins., Author Summary Modified live virus (MLV) vaccines have been effective in reducing disease burden and economic loss caused by Newcastle Disease (ND) in domestic poultry. Because the vaccine is a live virus, it is transmissible among birds. Thus, vaccination strategies have the potential to impact the evolutionary genetics of wild type strains of aPMV-1 including those that cause ND. In this report, we provided evidence that viruses isolated from wild and domestic birds have recombined with vaccine strains, because vaccinated birds are protected from disease but not infection with other strains of aPMV-1. Despite the use of vaccines since the 1950s, the population size of the strain from which the most widely used vaccine was derived has steadily increased. In contrast, other contemporary genotypes, which emerged in the 1960s, experienced a decline in population size in 1998, which may reflect a change in poultry farming practices or disease. Vaccination imposed a unique selection profile on the genotypes derived from the vaccine-related strains when compared with nonvaccine-related strains. Although modified live viruses are important for controlling Newcastle Disease, the potential of vaccination strategies to change viral diversity and population dynamics should be considered.

Enterovirus 71 (EV71) is one of the major causative agents of hand, foot, and mouth disease, a common febrile disease in children; however, EV71 has been also associated with various neurological diseases including fatal cases in large EV71 outbreaks particularly in the Asia Pacific region. Recently we identified human P-selectin glycoprotein ligand-1 (PSGL-1) as a cellular receptor for entry and replication of EV71 in leukocytes. PSGL-1 is a sialomucin expressed on the surface of leukocytes, serves as a high affinity counterreceptor for selectins, and mediates leukocyte rolling on the endothelium. The PSGL-1–P-selectin interaction requires sulfation of at least one of three clustered tyrosines and an adjacent O-glycan expressing sialyl Lewis x in an N-terminal region of PSGL-1. To elucidate the molecular basis of the PSGL-1–EV71 interaction, we generated a series of PSGL-1 mutants and identified the post-translational modifications that are critical for binding of PSGL-1 to EV71. We expressed the PSGL-1 mutants in 293T cells and the transfected cells were assayed for their abilities to bind to EV71 by flow cytometry. We found that O-glycosylation on T57, which is critical for PSGL-1–selectin interaction, is not necessary for PSGL-1 binding to EV71. On the other hand, site-directed mutagenesis at one or more potential tyrosine sulfation sites in the N-terminal region of PSGL-1 significantly impaired PSGL-1 binding to EV71. Furthermore, an inhibitor of sulfation, sodium chlorate, blocked the PSGL-1–EV71 interaction and inhibited PSGL-1-mediated viral replication of EV71 in Jurkat T cells in a dose-dependent manner. Thus, the results presented in this study reveal that tyrosine sulfation, but not O-glycosylation, in the N-terminal region of PSGL-1 may facilitate virus entry and replication of EV71 in leukocytes., Author Summary Enterovirus 71 (EV71) is a major causative agent of hand, foot, and mouth disease and a diverse array of neurological diseases, including fatal encephalitis, in children. EV71 has increasingly caused large outbreaks of hand, foot, and mouth disease particularly in the Asia-Pacific region. Recently, we identified human P-selectin glycoprotein ligand-1 (PSGL-1) as a functional receptor for EV71. PSGL-1 on immune cells is a key molecule involved in early inflammatory events and the PSGL-1–selectin interaction is regulated by post-translational modifications of PSGL-1. Here, we found that a post-translational modification, tyrosine sulfation, at the N-terminal region of PSGL-1 is critical for its binding to EV71 and subsequent viral replication in lymphocytes. Important roles for tyrosine sulfation in protein-protein interactions have been widely accepted; however, involvement of tyrosine sulfation of the receptor in the virus-receptor interaction has been reported only for HIV-1. Therefore, this is the second and unique example of the involvement of tyrosine sulfation in specific virus-receptor interactions. Our results shed new light on biological roles for tyrosine-sulfated proteins in cell tropism and the pathogenesis of EV71.

Background: New vaccine designs are needed to control diseases associated with antigenically variable RNA viruses. Foot-and-mouth disease (FMD) is a highly contagious disease of livestock that inflicts severe economic losses. Although the current whole-virus chemically inactivated vaccine has proven effective, it has led to new outbreaks of FMD because of incomplete inactivation of the virus or the escape of infectious virus from vaccine production premises. We have previously shown that serial passages of FMD virus (FMDV) C-S8c1 at high multiplicity of infection in cell culture resulted in virus populations consisting of defective genomes that are infectious by complementation (termed C-S8p260). Principal Finding: Here we evaluate the immunogenicity of C-S8p260, first in a mouse model system to establish a proof of principle, and second, in swine, the natural host of FMDV C-S8c1. Mice were completely protected against a lethal challenge with FMDV C-S8c1, after vaccination with a single dose of C-S8p260. Pigs immunized with different C-S8p260 doses and challenged with FMDV C-S8c1 either did not develop any clinical signs or showed delayed and mild disease symptoms. C-S8p260 induced high titers of both FMDV-specific, neutralizing antibodies and activated FMDV-specific T cells in swine, that correlated with solid protection against FMDV. Conclusions: The defective virus-based vaccine did not produce detectable levels of transmissible FMDV. Therefore, a segmented, replication-competent form of a virus, such as FMDV C-S8p260, can provide the basis of a new generation of attenuated antiviral vaccines with two safety barriers. The design can be extended to any viral pathogen that encodes trans-acting gene products, allowing complementation between replication-competent, defective forms. © 2010 Rodríguez-Calvo et al., Ministerio de Ciencia e Innovación, Spain, and European Union, Network of Excellence, EPIZONE (Contract # FOOD-CT-2006-016236). CIBERehd (Centro de Investigación Biomédica en Red de Enfermedades Hepáticas y Digestivas) is funded by Instituto de Salud Carlos III. Work at Centro de Biología Molecular ‘‘Severo Ochoa’’ (CISC-UAM) was supported by an institutional grant from Fundación Ramón Areces

Background Breast milk transmission of HIV-1 remains a major route of pediatric infection. Defining the characteristics of viral variants to which breastfeeding infants are exposed is important for understanding the genetic bottleneck that occurs in the majority of mother-to-child transmissions. The blood-milk epithelial barrier markedly restricts the quantity of HIV-1 in breast milk, even in the absence of antiretroviral drugs. The basis of this restriction and the genetic relationship between breast milk and blood variants are not well established. Methodology/Principal Findings We compared 356 HIV-1 subtype C gp160 envelope (env) gene sequences from the plasma and breast milk of 13 breastfeeding women. A trend towards lower viral population diversity and divergence in breast milk was observed, potentially indicative of clonal expansion within the breast. No differences in potential N-linked glycosylation site numbers or in gp160 variable loop amino acid lengths were identified. Genetic compartmentalization was evident in only one out of six subjects in whom contemporaneously obtained samples were studied. However, in samples that were collected 10 or more days apart, six of seven subjects were classified as having compartmentalized viral populations, highlighting the necessity of contemporaneous sampling for genetic compartmentalization studies. We found evidence of CXCR4 co-receptor using viruses in breast milk and blood in nine out of the thirteen subjects, but no evidence of preferential localization of these variants in either tissue. Conclusions/Significance Despite marked restriction of HIV-1 quantities in milk, our data indicate intermixing of virus between blood and breast milk. Thus, we found no evidence that a restriction in viral genotype diversity in breast milk accounts for the genetic bottleneck observed following transmission. In addition, our results highlight the rapidity of HIV-1 env evolution and the importance of sample timing in analyses of gene flow.

Human rhinoviruses (HRVs) evolve rapidly due in part to their error-prone RNA polymerase. Knowledge of the diversity of HRV populations emerging during the course of a natural infection is essential and represents a basis for the design of future potential vaccines and antiviral drugs. To evaluate HRV evolution in humans, nasal wash samples were collected daily for five days from 15 immunocompetent volunteers experimentally infected with a reference stock of HRV-39. In parallel, HeLa-OH cells were inoculated to compare HRV evolution in vitro. Nasal wash in vivo assessed by real-time PCR showed a viral load that peaked at 48-72 h. Ultra-deep sequencing was used to compare the low-frequency mutation populations present in the HRV-39 inoculum in two human subjects and one HeLa-OH supernatant collected 5 days post-infection. The analysis revealed hypervariable mutation locations in VP2, VP3, VP1, 2C and 3C genes and conserved regions in VP4, 2A, 2B, 3A, 3B and 3D genes. These results were confirmed by classical sequencing of additional samples, both from inoculated volunteers and independent cell infections, and suggest that HRV inter-host transmission is not associated with a strong bottleneck effect. A specific analysis of the VP1 capsid gene of 15 human cases confirmed the high mutation incidence in this capsid region, but not in the antiviral drug-binding pocket. We could also estimate a mutation frequency in vivo of 3.4x10(-4) mutations/nucleotides and 3.1x10(-4) over the entire ORF and VP1 gene, respectively. In vivo, HRV generate new variants rapidly during the course of an acute infection due to mutations that accumulate in hot spot regions located at the capsid level, as well as in 2C and 3C genes.

A novel, swine-origin influenza H1N1 virus (H1N1pdm) caused the first pandemic of the 21st century. This pandemic, although efficient in transmission, is mild in virulence. This atypical mild pandemic season has raised concerns regarding the potential of this virus to acquire additional virulence markers either through further adaptation or possibly by immune pressure in the human host. Using the mouse model we generated, within a single round of infection with A/California/04/09/H1N1 (Ca/04), a virus lethal in mice—herein referred to as mouse-adapted Ca/04 (ma-Ca/04). Five amino acid substitutions were found in the genome of ma-Ca/04: 3 in HA (D131E, S186P and A198E), 1 in PA (E298K) and 1 in NP (D101G). Reverse genetics analyses of these mutations indicate that all five mutations from ma-Ca/04 contributed to the lethal phenotype; however, the D131E and S186P mutations—which are also found in the 1918 and seasonal H1N1 viruses—in HA alone were sufficient to confer virulence of Ca/04 in mice. HI assays against H1N1pdm demonstrate that the D131E and S186P mutations caused minor antigenic changes and, likely, affected receptor binding. The rapid selection of ma-Ca/04 in mice suggests that a virus containing this constellation of amino acids might have already been present in Ca/04, likely as minor quasispecies., Author Summary The 19th century experienced three major influenza pandemics: the Spanish flu of 1918 (H1N1), the Asian flu of 1957 (H2N2) and the Hong Kong flu of 1968 (H3N2). These pandemics were introduced into the human population through the accumulation of avian and human influenza genes (genetic reassortment), creating with each pandemic a novel influenza virus, one that the human population had not been exposed to. These pandemics were associated with high morbidity and mortality; the 1918 pandemic was responsible for an estimated 40 million deaths. In April 2009, a novel, swine-origin influenza H1N1 virus (H1N1pdm) caused the first influenza pandemic of the 21st century. This pandemic was a result of genetic reassortment between not only human and avian influenza genes, but also swine influenza genes. This H1N1 pandemic, although efficient in human-to-human transmission, differed greatly from the previous pandemics in its mild virulence. This atypical mild pandemic season has raised concerns regarding the potential of this virus to become more virulent. Using a mouse model, we were able to demonstrate that this pandemic strain is amenable to mutations that lead to a more virulent virus.

Cell surface carbohydrates play an important role in virus entry and intracellular trafficking. Bovine Adeno-Associated Virus (BAAV) uses plasma membrane gangliosides for transduction and infection. In addition, independent of the infectious pathway, BAAV also has the ability to pass through barrier epithelia and endothelia using a transcytosis pathway dependent upon the presence of cell surface carbohydrates. Thus, in order to better define the carbohydrate interactions that are necessary for BAAV infection or transcytosis, a glycan microarray composed of both natural and synthetic carbohydrates was probed with HA-tagged BAAV particles. This identified chitotriose, a trimer of beta-1-4-linked N-acetyl glucosamine, as having an interaction with BAAV. Competition experiments showed that the BAAV interaction with this carbohydrate is not necessary for infection but is instead important in the transcytosis pathway. The beta-1-4-linked N-acetyl glucosamine modification has been reported on gp96, a glycoprotein involved in the transcytosis of bacteria and toxins. Significantly, immunoprecipitation and competition experiments with an anti-gp96 antibody and a soluble form of gp96, respectively, showed this glycoprotein can also interact with BAAV to serve as a receptor for its transcytosis.

Over the last fifteen years there have been five pandemics of norovirus (NoV) associated gastroenteritis, and the period of stasis between each pandemic has been progressively shortening. NoV is classified into five genogroups, which can be further classified into 25 or more different human NoV genotypes; however, only one, genogroup II genotype 4 (GII.4), is associated with pandemics. Hence, GII.4 viruses have both a higher frequency in the host population and greater epidemiological fitness. The aim of this study was to investigate if the accuracy and rate of replication are contributing to the increased epidemiological fitness of the GII.4 strains. The replication and mutation rates were determined using in vitro RNA dependent RNA polymerase (RdRp) assays, and rates of evolution were determined by bioinformatics. GII.4 strains were compared to the second most reported genotype, recombinant GII.b/GII.3, the rarely detected GII.3 and GII.7 and as a control, hepatitis C virus (HCV). The predominant GII.4 strains had a higher mutation rate and rate of evolution compared to the less frequently detected GII.b, GII.3 and GII.7 strains. Furthermore, the GII.4 lineage had on average a 1.7-fold higher rate of evolution within the capsid sequence and a greater number of non-synonymous changes compared to other NoVs, supporting the theory that it is undergoing antigenic drift at a faster rate. Interestingly, the non-synonymous mutations for all three NoV genotypes were localised to common structural residues in the capsid, indicating that these sites are likely to be under immune selection. This study supports the hypothesis that the ability of the virus to generate genetic diversity is vital for viral fitness., Author Summary Since 1995, norovirus has caused five pandemics of acute gastroenteritis. These pandemics spread across the globe within a few months, causing great economic burden on society due to medical and social expenses. Norovirus, like influenza virus, has over 40 genotypes circulating within the population at the same time. However, it is only a single genotype, known as genogroup II genotype 4 (GII.4), that causes mass outbreaks and pandemics. Very little research has been conducted to determine why GII.4 viruses can cause pandemics. Consequently, we compared the evolution properties of several pandemic GII.4 strains to non-pandemic strains and found that the GII.4 viruses were undergoing evolution at a much higher rate than the non-pandemic norovirus strains. This phenomenon is similar to influenza virus, where an increase in antigenic drift has been associated with increased outbreaks. This discovery has important implications in understanding norovirus incidence and also the development of a vaccine and treatment for norovirus.

Hepatitis A virus (HAV), the prototype of genus Hepatovirus, has several unique biological characteristics that distinguish it from other members of the Picornaviridae family. Among these, the need for an intact eIF4G factor for the initiation of translation results in an inability to shut down host protein synthesis by a mechanism similar to that of other picornaviruses. Consequently, HAV must inefficiently compete for the cellular translational machinery and this may explain its poor growth in cell culture. In this context of virus/cell competition, HAV has strategically adopted a naturally highly deoptimized codon usage with respect to that of its cellular host. With the aim to optimize its codon usage the virus was adapted to propagate in cells with impaired protein synthesis, in order to make tRNA pools more available for the virus. A significant loss of fitness was the immediate response to the adaptation process that was, however, later on recovered and more associated to a re-deoptimization rather than to an optimization of the codon usage specifically in the capsid coding region. These results exclude translation selection and instead suggest fine-tuning translation kinetics selection as the underlying mechanism of the codon usage bias in this specific genome region. Additionally, the results provide clear evidence of the Red Queen dynamics of evolution since the virus has very much evolved to re-adapt its codon usage to the environmental cellular changing conditions in order to recover the original fitness., Author Summary Each organism has a specific codon usage signature. Translational selection i.e., selection for the codon adaptation to the tRNA pools, is one of the driving forces of codon bias. In the virus world, this implies an adjustment of the virus codon usage to that of the host cell. Hepatitis A virus appears as an exception to the rule, with a highly deoptimized codon usage, suggesting that translational selection is not the underlying mechanism of its codon bias. However, since the virus lacks a mechanism of cellular protein synthesis inhibition, the deoptimized codon usage may be envisaged as a hawk (cell) and dove (hepatitis A virus) competition strategy for tRNAs and translational selection as well. To confirm this possibility, we artificially induced cell protein synthesis shut-off, thus increasing the tRNA pool availability for the virus, and we took advantage of the quasispecies dynamics to elucidate changes in its codon usage. Virus adaptation to the drug results in a re-deoptimization of codon usage in the capsid region, suggesting a requirement of a slow translation rate, i.e., a translation kinetic selection, instead of a translational selection associated with an optimization of the codon usage. Translation kinetics control is based on the right combination of codons (common and rare) that allows a regulated ribosome traffic rate ensuring the proper protein folding. Capsid folding is critical for a virus transmitted through the fecal-oral route with long extracorporeal periods.

Background Mother-to-child transmission of HIV-1 remains a significant problem in the resource-constrained settings where anti-retroviral therapy is still not widely available. Understanding the earliest events during HIV-1 transmission and characterizing the newly transmitted or founder virus is central to intervention efforts. In this study, we analyzed the viral env quasispecies of six mother-infant transmission pairs (MIPs) and characterized the genetic features of envelope glycoprotein that could influence HIV-1 subtype C perinatal transmission. Methodology and Findings The V1-V5 region of env was amplified from 6 MIPs baseline samples and 334 DNA sequences in total were analyzed. A comparison of the viral population derived from the mother and infant revealed a severe genetic bottleneck occurring during perinatal transmission, which was characterized by low sequence diversity in the infant. Phylogenetic analysis indicates that most likely in all our infant subjects a single founder virus was responsible for establishing infection. Furthermore, the newly transmitted viruses from the infant had significantly fewer potential N-linked glycosylation sites in Env V1-V5 region and showed a propensity to encode shorter variable loops compared to the nontransmitted viruses. In addition, a similar intensity of selection was seen between mothers and infants with a higher rate of synonymous (dS) compared to nonsynonymous (dN) substitutions evident (dN/dS

[Background] Fitness recovery of HIV-1 “in vitro” was studied using viral clones that had their fitness decreased as a result of plaque-to-plaque passages., [Principal Findings] After ten large population passages, the viral populations showed an average increase of fitness, although with wide variations among clones. While 5 clones showed significant fitness increases, 3 clones showed increases that were only marginally significant (p, [Conclusions] Thus, the initial fitness recovery of debilitated HIV-1 clones was mediated by an increase in quasispecies heterogeneity. This observation, together with the invariance of the consensus sequence despite fitness increases demonstrates the relevance of quasispecies heterogeneity in the evolution of HIV-1 in cell culture., Work in the Centro Nacional de Microbiología was supported by grants, SAF 2005/03833, SAF 2007/61036 by the Plan Nacional del SIDA and in part by the Red Temática Cooperativa de investigación en SIDA (Red de grupos 173) of the Fondo de Investigaciones Sanitarias de la Seguridad Social (FISss) and Fundación para la Investigación y la Prevención del SIDA en España (FIPSE) grant 36558/06, 36641/07, 36779/08, 36641/07 and 360766/09. Work in the Centro de Biología Molecular “Severo Ochoa” was supported by grant BFU 2008-02816/BMC from Ministerio de Ciencia e Innovación, FIPSE 36558/06 and 360766/09 and Fundación Ramón Areces. Centro de Investigaciones Biomédica en Red (CIBER): enfermedades hepáticas y digestivas (ehd) is funded by Instituto de Salud Carlos III. The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.

Wild aquatic birds in the Orders Anseriformes and Charadriiformes are the main reservoir hosts perpetuating the genetic pool of all influenza A viruses, including pandemic viruses. High viral loads in feces of infected birds permit a fecal-oral route of transmission. Numerous studies have reported the isolation of avian influenza viruses (AIVs) from surface water at aquatic bird habitats. These isolations indicate aquatic environments have an important role in the transmission of AIV among wild aquatic birds. However, the progressive dilution of infectious feces in water could decrease the likelihood of virus/host interactions. To evaluate whether alternate mechanisms facilitate AIV transmission in aquatic bird populations, we investigated whether the preen oil gland secretions by which all aquatic birds make their feathers waterproof could support a natural mechanism that concentrates AIVs from water onto birds' bodies, thus, representing a possible source of infection by preening activity. We consistently detected both viral RNA and infectious AIVs on swabs of preened feathers of 345 wild mallards by using reverse transcription-polymerase chain reaction (RT-PCR) and virus-isolation (VI) assays. Additionally, in two laboratory experiments using a quantitative real-time (qR) RT-PCR assay, we demonstrated that feather samples (n = 5) and cotton swabs (n = 24) experimentally impregnated with preen oil, when soaked in AIV-contaminated waters, attracted and concentrated AIVs on their surfaces. The data presented herein provide information that expands our understanding of AIV ecology in the wild bird reservoir system.

Background The Influenza A pandemic H1N1 2009 (H1N1pdm) virus appeared in India in May 2009 and thereafter outbreaks with considerable morbidity and mortality have been reported from many parts of the country. Continuous monitoring of the genetic makeup of the virus is essential to understand its evolution within the country in relation to global diversification and to track the mutations that may affect the behavior of the virus. Methods H1N1pdm viruses were isolated from both recovered and fatal cases representing major cities and sequenced. Phylogenetic analyses of six concatenated whole genomes and the hemagglutinin (HA) gene of seven more isolates from May-September 2009 was performed with reference to 685 whole genomes of global isolates available as of November 24, 2009. Molecular characterization of all the 8 segments was carried out for known pathogenic markers. Results The first isolate of May 2009 belonged to clade 5. Although clade 7 was the dominant H1N1pdm lineage in India, both clades 6 and 7 were found to be co-circulating. The neuraminidase of all the Indian isolates possessed H275, the marker for sensitivity to the neuraminidase inhibitor Oseltamivir. Some of the mutations in HA are at or in the vicinity of antigenic sites and may therefore be of possible antigenic significance. Among these a D222G mutation in the HA receptor binding domain was found in two of the eight Indian isolates obtained from fatal cases. Conclusions The majority of the 13 Indian isolates grouped in the globally most widely circulating H1N1pdm clade 7. Further, correlations of the mutations specific to clade 7 Indian isolates to viral fitness and adaptability in the country remains to be understood. The D222G mutation in HA from isolates of fatal cases needs to be studied for pathogenicity.

Infectious hypodermal and hematopoietic necrosis virus (IHHNV) is a widely distributed single-stranded DNA parvovirus that has been responsible for major losses in wild and farmed penaeid shrimp populations on the northwestern Pacific coast of Mexico since the early 1990's. IHHNV has been considered a slow-evolving, stable virus because shrimp populations in this region have recovered to pre-epizootic levels, and limited nucleotide variation has been found in a small number of IHHNV isolates studied from this region. To gain insight into IHHNV evolutionary and population dynamics, we analyzed IHHNV capsid protein gene sequences from 89 Penaeus shrimp, along with 14 previously published sequences. Using Bayesian coalescent approaches, we calculated a mean rate of nucleotide substitution for IHHNV that was unexpectedly high (1.39x10(-4) substitutions/site/year) and comparable to that reported for RNA viruses. We found more genetic diversity than previously reported for IHHNV isolates and highly significant subdivision among the viral populations in Mexican waters. Past changes in effective number of infections that we infer from Bayesian skyline plots closely correspond to IHHNV epizootiological historical records. Given the high evolutionary rate and the observed regional isolation of IHHNV in shrimp populations in the Gulf of California, we suggest regular monitoring of wild and farmed shrimp and restriction of shrimp movement as preventative measures for future viral outbreaks.

Adeno-associated virus type 2 (AAV) is known to establish latency by preferential integration in human chromosome 19q13.42. The AAV non-structural protein Rep appears to target a site called AAVS1 by simultaneously binding to Rep-binding sites (RBS) present on the AAV genome and within AAVS1. In the absence of Rep, as is the case with AAV vectors, chromosomal integration is rare and random. For a genome-wide survey of wildtype AAV integration a linker-selection-mediated (LSM)-PCR strategy was designed to retrieve AAV-chromosomal junctions. DNA sequence determination revealed wildtype AAV integration sites scattered over the entire human genome. The bioinformatic analysis of these integration sites compared to those of rep-deficient AAV vectors revealed a highly significant overrepresentation of integration events near to consensus RBS. Integration hotspots included AAVS1 with 10% of total events. Novel hotspots near consensus RBS were identified on chromosome 5p13.3 denoted AAVS2 and on chromsome 3p24.3 denoted AAVS3. AAVS2 displayed seven independent junctions clustered within only 14 bp of a consensus RBS which proved to bind Rep in vitro similar to the RBS in AAVS3. Expression of Rep in the presence of rep-deficient AAV vectors shifted targeting preferences from random integration back to the neighbourhood of consensus RBS at hotspots and numerous additional sites in the human genome. In summary, targeted AAV integration is not as specific for AAVS1 as previously assumed. Rather, Rep targets AAV to integrate into open chromatin regions in the reach of various, consensus RBS homologues in the human genome., Author Summary This is the first unbiased genome-wide analysis of wildtype AAV integration combined with a thorough bioinformatic analysis of preferred genomic motifs and patterns in the neighbourhood of the integration sites identified. The preference of Rep-dependent AAV integration near multiple consensus Rep-binding sites was lost in the case of AAV vector integration in the absence of Rep expression. Our findings challenge the commonly accepted notion of site-specific AAV targeting to AAVS1 on chromosome 19q13.42. Although AAVS1 contains a canonical Rep-binding site, numerous additional sites including the newly identified hotspots AAVS2 on chromosome 5p13.3 and AAVS3 on chromosome 3p24.3 harbour functional Rep-binding sites suitable for AAV integration. AAV vectors are quickly moving forward in the clinic and Rep-dependent vector targeting strategies are being actively pursued. Detailed information of AAV wildtype versus recombinant AAV vector integration sites and preferences are needed to evaluate the safety profile of AAV vectors in gene therapy.

Background: Success of a viral infection requires that each infected cell delivers a sufficient number of infectious particles to allow new rounds of infection. In picornaviruses, viral replication is initiated by the viral polymerase and a viral-coded protein, termed VPg, that primes RNA synthesis. Foot-and-mouth disease virus (FMDV) is exceptional among picornaviruses in that its genome encodes 3 copies of VPg. Why FMDV encodes three VPgs is unknown. Methodology and Principal Findings: We have constructed four mutant FMDVs that encode only one VPg: either VPg1, VPg3, or two chimeric versions containing part of VPg1 and VPg3. All mutants, except that encoding only VPg1, were replication-competent. Unexpectedly, despite being replication-competent, the mutants did not form plaques on BHK-21 cell monolayers. The one-VPg mutant FMDVs released lower amounts of encapsidated viral RNA to the extracellular environment than wild type FMDV, suggesting that deficient plaque formation was associated with insufficient release of infectious progeny. Mutant FMDVs subjected to serial passages in BHK-21 cells regained plaque-forming capacity without modification of the number of copies of VPg. Substitutions in non-structural proteins 2C, 3A and VPg were associated with restoration of plaque formation. Specifically, replacement R55W in 2C was repeatedly found in several mutant viruses that had regained competence in plaque development. The effect of R55W in 2C was to mediate an increase in the extracellular viral RNA release without a detectable increase of total viral RNA that correlated with an enhanced capacity to alter and detach BHK-21 cells from the monolayer, the first stage of cell killing. Conclusions: The results link the VPg copies in the FMDV genome with the cytopathology capacity of the virus, and have unveiled yet another function of 2C: modulation of picornavirus cell-to-cell transmission. Implications for picornaviruses pathogenesis are discussed. © 2010 Arias et al., BFU2008-02816/BMC from Ministerio de Ciencia e Innovacion (MICINN), and Fundacion R. Areces. CIBERehd (Centro de Investigacion Biomedica en Red de Enfermedades Hepaticas y Digestivas) is funded by Instituto de Salud Carlos III

A novel swine-origin pandemic influenza A(H1N1) virus (H1N1pdm, also referred to as S-OIV) was identified as the causative agent of the 21(st) century's first influenza pandemic, but molecular features conferring its ability of human-to-human transmission has not been identified. Here we compared the protein sequences of 2009 H1N1pdm strains with those causing other pandemics and the viruses isolated from humans, swines and avians, and then analyzed the mutation trend of the residues at the signature and non-signature positions, which are species- and non-species-associated, respectively, in the proteins of H1N1pdm during the pandemic of 2009. We confirmed that the host-specific genomic signatures of 2009 H1N1pdm, which are mainly swine-like, were highly identical to those of the 1918 H1N1pdm. During the short period of time when the pandemic alert level was raised from phase 4 to phase 6, one signature residue at the position of NP-100 mutated from valine to isoleucine. Four non-signature residues, at positions NA-91, NA-233, HA-206, and NS1-123, also changed during the epidemic in 2009. All these mutant residues, except that at NA-91, are located in the viral functional domains, suggesting that they may play roles in the human adaption and virulence of 2009 H1N1pdm.

The declaration of the human influenza A pandemic (H1N1) 2009 (H1N1/09) raised important questions, including origin and host range [1], [2]. Two of the three pandemics in the last century resulted in the spread of virus to pigs (H1N1, 1918; H3N2, 1968) with subsequent independent establishment and evolution within swine worldwide [3]. A key public and veterinary health consideration in the context of the evolving pandemic is whether the H1N1/09 virus could become established in pig populations [4]. We performed an infection and transmission study in pigs with A/California/07/09. In combination, clinical, pathological, modified influenza A matrix gene real time RT-PCR and viral genomic analyses have shown that infection results in the induction of clinical signs, viral pathogenesis restricted to the respiratory tract, infection dynamics consistent with endemic strains of influenza A in pigs, virus transmissibility between pigs and virus-host adaptation events. Our results demonstrate that extant H1N1/09 is fully capable of becoming established in global pig populations. We also show the roles of viral receptor specificity in both transmission and tissue tropism. Remarkably, following direct inoculation of pigs with virus quasispecies differing by amino acid substitutions in the haemagglutinin receptor-binding site, only virus with aspartic acid at position 225 (225D) was detected in nasal secretions of contact infected pigs. In contrast, in lower respiratory tract samples from directly inoculated pigs, with clearly demonstrable pulmonary pathology, there was apparent selection of a virus variant with glycine (225G). These findings provide potential clues to the existence and biological significance of viral receptor-binding variants with 225D and 225G during the 1918 pandemic [5].

Resistance of viruses to mutagenic agents is an important problem for the development of lethal mutagenesis as an antiviral strategy. Previous studies with RNA viruses have documented that resistance to the mutagenic nucleoside analogue ribavirin (1-b-D-ribofuranosyl-1-H-1,2,4-triazole-3-carboxamide) is mediated by amino acid substitutions in the viral polymerase that either increase the general template copying fidelity of the enzyme or decrease the incorporation of ribavirin into RNA. Here we describe experiments that show that replication of the important picornavirus pathogen foot- and-mouth disease virus (FMDV) in the presence of increasing concentrations of ribavirin results in the sequential incorporation of three amino acid substitutions (M296I, P44S and P169S) in the viral polymerase (3D). The main biological effect of these substitutions is to attenuate the consequences of the mutagenic activity of ribavirin -by avoiding the biased repertoire of transition mutations produced by this purine analogue-and to maintain the replicative fitness of the virus which is able to escape extinction by ribavirin. This is achieved through alteration of the pairing behavior of ribavirin- triphosphate (RTP), as evidenced by in vitro polymerization assays with purified mutant 3Ds. Comparison of the three- dimensional structure of wild type and mutant polymerases suggests that the amino acid substitutions alter the position of the template RNA in the entry channel of the enzyme, thereby affecting nucleotide recognition. The results provide evidence of a new mechanism of resistance to a mutagenic nucleoside analogue which allows the virus to maintain a balance among mutation types introduced into progeny genomes during replication under strong mutagenic pressure. © 2010 Agudo et al.

Ribonucleotide reductases (RRs) are evolutionarily-conserved enzymes that catalyze the rate-limiting step during dNTP synthesis in mammals. RR consists of both large (R1) and small (R2) subunits, which are both required for catalysis by the R12R22 heterotetrameric complex. Poxviruses also encode RR proteins, but while the Orthopoxviruses infecting humans [e.g. vaccinia (VACV), variola, cowpox, and monkeypox viruses] encode both R1 and R2 subunits, the vast majority of Chordopoxviruses encode only R2 subunits. Using plaque morphology, growth curve, and mouse model studies, we investigated the requirement of VACV R1 (I4) and R2 (F4) subunits for replication and pathogenesis using a panel of mutant viruses in which one or more viral RR genes had been inactivated. Surprisingly, VACV F4, but not I4, was required for efficient replication in culture and virulence in mice. The growth defects of VACV strains lacking F4 could be complemented by genes encoding other Chordopoxvirus R2 subunits, suggesting conservation of function between poxvirus R2 proteins. Expression of F4 proteins encoding a point mutation predicted to inactivate RR activity but still allow for interaction with R1 subunits, caused a dominant negative phenotype in growth experiments in the presence or absence of I4. Co-immunoprecipitation studies showed that F4 (as well as other Chordopoxvirus R2 subunits) form hybrid complexes with cellular R1 subunits. Mutant F4 proteins that are unable to interact with host R1 subunits failed to rescue the replication defect of strains lacking F4, suggesting that F4-host R1 complex formation is critical for VACV replication. Our results suggest that poxvirus R2 subunits form functional complexes with host R1 subunits to provide sufficient dNTPs for viral replication. Our results also suggest that R2-deficient poxviruses may be selective oncolytic agents and our bioinformatic analyses provide insights into how poxvirus nucleotide metabolism proteins may have influenced the base composition of these pathogens., Author Summary Efficient genome replication is central to the virulence of all DNA viruses, including poxviruses. To ensure replication efficiency, many of the more virulent poxviruses encode their own nucleotide metabolism machinery, including ribonucleotide reductase (RR) enzymes, which act to provide ample DNA precursors for replication. RR enzymes require both large (R1) and small (R2) subunit proteins for activity. Curiously, some poxviruses only encode R2 subunits. Other poxviruses, such as the smallpox vaccine strain, vaccinia virus (VACV), encode both R1 and R2 subunits. We report here that the R2, but not the R1, subunit of VACV RR is required for efficient replication and virulence. We also provide evidence that several poxvirus R2 proteins form novel complexes with host R1 subunits and this interaction is required for efficient VACV replication in primate cells. Our study explains why some poxviruses only encode R2 subunits and identifies a role for these proteins in poxvirus pathogenesis. Furthermore, we provide evidence that mutant poxviruses unable to generate R2 proteins may become entirely dependent upon host RR activity. This may restrict their replication to cells that over-express RR proteins such as cancer cells, making them potential therapeutics for human malignancies.

Hop stunt was a mysterious disorder that first emerged in the 1940s in commercial hops in Japan. To investigate the origin of this disorder, we infected hops with natural Hop stunt viroid (HpSVd) isolates derived from four host species (hop, grapevine, plum and citrus), which except for hop represent possible sources of the ancestral viroid. These plants were maintained for 15 years, then analyzed the HpSVd variants present. Here we show that the variant originally found in cultivated grapevines gave rise to various combinations of mutations at positions 25, 26, 54, 193, and 281. However, upon prolonged infection, these variants underwent convergent evolution resulting in a limited number of adapted mutants. Some of them showed nucleotide sequences identical to those currently responsible for hop stunt epidemics in commercial hops in Japan, China, and the United States. Therefore, these results indicate that we have successfully reproduced the original process by which a natural HpSVd variant naturally introduced into cultivated hops was able to mutate into the HpSVd variants that are currently present in commercial hops. Furthermore, and importantly, we have identified cultivated grapevines as a symptomless reservoir in which HSVd can evolve and be transmitted to hop crops to cause epidemics.

Many pathogens associated with chronic infections evolve so rapidly that strains found late in an infection have little in common with the initial strain. This raises questions at different levels of analysis because rapid within-host evolution affects the course of an infection, but it can also affect the possibility for natural selection to act at the between-host level. We present a nested approach that incorporates within-host evolutionary dynamics of a rapidly mutating virus (hepatitis C virus) targeted by a cellular cross-reactive immune response, into an epidemiological perspective. The viral trait we follow is the replication rate of the strain initiating the infection. We find that, even for rapidly evolving viruses, the replication rate of the initial strain has a strong effect on the fitness of an infection. Moreover, infections caused by slowly replicating viruses have the highest infection fitness (i.e., lead to more secondary infections), but strains with higher replication rates tend to dominate within a host in the long-term. We also study the effect of cross-reactive immunity and viral mutation rate on infection life history traits. For instance, because of the stochastic nature of our approach, we can identify factors affecting the outcome of the infection (acute or chronic infections). Finally, we show that anti-viral treatments modify the value of the optimal initial replication rate and that the timing of the treatment administration can have public health consequences due to within-host evolution. Our results support the idea that natural selection can act on the replication rate of rapidly evolving viruses at the between-host level. It also provides a mechanistic description of within-host constraints, such as cross-reactive immunity, and shows how these constraints affect the infection fitness. This model raises questions that can be tested experimentally and underlines the necessity to consider the evolution of quantitative traits to understand the outcome and the fitness of an infection., Author Summary Rapidly mutating viruses, such as hepatitis C virus, can escape host immunity by generating new strains that avoid the immune system. Existing data support the idea that such within-host evolution affects the outcome of the infection. Few theoretical models address this question and most follow viral diversity or qualitative traits, such as drug resistance. Here, we study the evolution of two virus quantitative traits—the replication rate and the ability to be recognised by the immune response—during an infection. We develop an epidemiological framework where transmission events are driven by within-host dynamics. We find that the replication rate of the virus that initially infects the host has a strong influence on the epidemiological success of the disease. Furthermore, we show that the cross-reactive immune response is key to determining the outcome of the infection (acute or chronic). Finally, we show that the timing of the start of an anti-viral treatment has a strong effect on viral evolution, which impacts the efficiency of the treatment. Our analysis suggests a new mechanism to explain infection outcomes and proposes testable predictions that can drive future experimental approaches.

St. Louis encephalitis virus (SLEV; Flaviviridae, flavivirus) was the major cause of epidemic flaviviral encephalitis in the U.S. prior to the introduction of West Nile virus (WNV) in 1999. However, outbreaks of SLEV have been significantly more limited then WNV in terms of levels of activity and geographic dispersal. One possible explanation for these variable levels of activity is that differences in the potential for each virus to adapt to its host cycle exist. The need for arboviruses to replicate in disparate hosts is thought to result in constraints on both evolution and host-specific adaptation. If cycling is the cause of genetic stability observed in nature and arboviruses lack host specialization, then sequential passage should result in both the accumulation of mutations and specialized viruses better suited for replication in that host. Previous studies suggest that WNV and SLEV differ in capacity for both genetic change and host specialization, and in the costs each accrues from specializing. In an attempt to clarify how selective pressures contribute to epidemiological patterns of WNV and SLEV, we evaluated mutant spectra size, consensus genetic change, and phenotypic changes for SLEV in vivo following 20 sequential passages via inoculation in either Culex pipiens mosquitoes or chickens. Results demonstrate that the capacity for genetic change is large for SLEV and that the size of the mutant spectrum is host-dependent using our passage methodology. Despite this, a general lack of consensus change resulted from passage in either host, a result that contrasts with the idea that constraints on evolution in nature result from host cycling alone. Results also suggest that a high level of adaptation to both hosts already exists, despite host cycling. A strain significantly more infectious in chickens did emerge from one lineage of chicken passage, yet other lineages and all mosquito passage strains did not display measurable host-specific fitness gains. In addition, increased infectivity in chickens did not decrease infectivity in mosquitoes, which further contrasts the concept of fitness trade-offs for arboviruses.

Widespread infection of highly pathogenic avian influenza A H5N1 was reported from backyard and commercial poultry in West Bengal (WB), an eastern state of India in early 2008. Infection gradually spread to Tripura, Assam and Sikkim, the northeastern states, with 70 outbreaks reported between January 2008 and May 2009. Whole genome sequence analysis of three isolates from WB, one isolate from Tripura along with the analysis of hemagglutinin (HA) and neuraminidase (NA) genes of 17 other isolates was performed during this study. In the HA gene phylogenetic tree, all the 2008-09 Indian isolates belonged to EMA3 sublineage of clade 2.2. The closest phylogenetic relationship was found to be with the 2007-09 isolates from Bangladesh and not with the earlier 2006 and 2007 Indian isolates implying a third introduction into the country. The receptor-binding pocket of HA1 of two isolates from WB showed S221P mutation, one of the markers predicted to be associated with human receptor specificity. Two substitutions E119A (2 isolates of WB) and N294S (2 other isolates of WB) known to confer resistance to NA inhibitors were observed in the active site of neuraminidase. Several additional mutations were observed within the 2008-09 Indian isolates indicating genetic diversification. Overall, the study is indicative of a possible endemicity in the eastern and northeastern parts of the country, demanding active surveillance specifically in view of the critical mutations that have been observed in the influenza A H5N1 viruses.

Tick-borne encephalitis virus (TBEV) causes human epidemics across Eurasia. Clinical manifestations range from inapparent infections and fevers to fatal encephalitis but the factors that determine disease severity are currently undefined. TBEV is characteristically a hemagglutinating (HA) virus; the ability to agglutinate erythrocytes tentatively reflects virion receptor/fusion activity. However, for the past few years many atypical HA-deficient strains have been isolated from patients and also from the natural European host tick, Ixodes persulcatus. By analysing the sequences of HA-deficient strains we have identified 3 unique amino acid substitutions (D67G, E122G or D277A) in the envelope protein, each of which increases the net charge and hydrophobicity of the virion surface. Therefore, we genetically engineered virus mutants each containing one of these 3 substitutions; they all exhibited HA-deficiency. Unexpectedly, each genetically modified non-HA virus demonstrated increased TBEV reproduction in feeding Ixodes ricinus, not the recognised tick host for these strains. Moreover, virus transmission efficiency between infected and uninfected ticks co-feeding on mice was also intensified by each substitution. Retrospectively, the mutation D67G was identified in viruses isolated from patients with encephalitis. We propose that the emergence of atypical Siberian HA-deficient TBEV strains in Europe is linked to their molecular adaptation to local ticks. This process appears to be driven by the selection of single mutations that change the virion surface thus enhancing receptor/fusion function essential for TBEV entry into the unfamiliar tick species. As the consequence of this adaptive mutagenesis, some of these mutations also appear to enhance the ability of TBEV to cross the human blood-brain barrier, a likely explanation for fatal encephalitis. Future research will reveal if these emerging Siberian TBEV strains continue to disperse westwards across Europe by adaptation to the indigenous tick species and if they are associated with severe forms of TBE.