Establishment and cryptic transmission of Zika virus in Brazil and the Americas

Transmission of Zika virus (ZIKV) in the Americas was first confirmed in May 2015 in northeast Brazil. Brazil has had the highest number of reported ZIKV cases worldwide (more than 200,000 by 24 December 2016) and the most cases associated with microcephaly and other birth defects (2,366 confirmed by 31 December 2016). Since the initial detection of ZIKV in Brazil, more than 45 countries in the Americas have reported local ZIKV transmission, with 24 of these reporting severe ZIKV-associated disease. However, the origin and epidemic history of ZIKV in Brazil and the Americas remain poorly understood, despite the value of this information for interpreting observed trends in reported microcephaly. Here we address this issue by generating 54 complete or partial ZIKV genomes, mostly from Brazil, and reporting data generated by a mobile genomics laboratory that travelled across northeast Brazil in 2016. One sequence represents the earliest confirmed ZIKV infection in Brazil. Analyses of viral genomes with ecological and epidemiological data yield an estimate that ZIKV was present in northeast Brazil by February 2014 and is likely to have disseminated from there, nationally and internationally, before the first detection of ZIKV in the Americas. Estimated dates for the international spread of ZIKV from Brazil indicate the duration of pre-detection cryptic transmission in recipient regions. The role of northeast Brazil in the establishment of ZIKV in the Americas is further supported by geographic analysis of ZIKV transmission potential and by estimates of the basic reproduction number of the virus.

Transgenic Virus Resistance in Crop-Wild Cucurbita pepo Does Not Prevent Vertical Transmission of Zucchini yellow mosaic virus

Zucchini yellow mosaic virus (ZYMV) is an economically important pathogen of cucurbits that is transmitted both horizontally and vertically. Although ZYMV is seed-transmitted in Cucurbita pepo, the potential for seed transmission in virus-resistant transgenic cultivars is not known. We crossed and backcrossed a transgenic squash cultivar with wild C. pepo, and determined whether seed-to-seedling transmission of ZYMV was possible in seeds harvested from transgenic backcrossed C. pepo. We then compared these transmission rates to those of non-transgenic (backcrossed and wild) C. pepo. The overall seed-to-seedling transmission rate in ZYMV was similar to those found in previous studies (1.37%), with no significant difference between transgenic backcrossed (2.48%) and non-transgenic (1.03%) backcrossed and wild squash. Fewer transgenic backcrossed plants had symptom development (7%) in comparison with all non-transgenic plants (26%) and may be instrumental in preventing yield reduction due to ZYMV. Our study shows that ZYMV is seed transmitted in transgenic backcrossed squash, which may affect the spread of ZYMV via the movement of ZYMV-infected seeds. Deep genome sequencing of the seed-transmitted viral populations revealed that 23% of the variants found in this study were present in other vertically transmitted ZYMV populations, suggesting that these variants may be necessary for seed transmission or are distributed geographically via seeds.

Analysis of viral (zucchini yellow mosaic virus) genetic diversity during systemic movement through a Cucurbita pepo vine

Determining the extent and structure of intra-host genetic diversity and the magnitude and impact of population bottlenecks is central to understanding the mechanisms of viral evolution. To determine the nature of viral evolution following systemic movement through a plant, we performed deep sequencing of 23 leaves that grew sequentially along a single Cucurbita pepo vine that was infected with zucchini yellow mosaic virus (ZYMV), and on a leaf that grew in on a side branch. Strikingly, of 112 genetic (i.e. sub-consensus) variants observed in the data set as a whole, only 22 were found in multiple leaves. Similarly, only three of the 13 variants present in the inoculating population were found in the subsequent leaves on the vine. Hence, it appears that systemic movement is characterized by sequential population bottlenecks, although not sufficient to reduce the population to a single virion as multiple variants were consistently transmitted between leaves. In addition, the number of variants within a leaf increases as a function of distance from the inoculated (source) leaf, suggesting that the circulating sap may serve as a continual source of virus. Notably, multiple mutational variants were observed in the cylindrical inclusion (CI) protein (known to be involved in both cell-to-cell and systemic movement of the virus) that were present in multiple (19/24) leaf samples. These mutations resulted in a conformational change, suggesting that they might confer a selective advantage in systemic movement within the vine. Overall, these data reveal that bottlenecks occur during systemic movement, that variants circulate in the phloem sap throughout the infection process, and that important conformational changes in CI protein may arise during individual infections.

Zucchini yellow mosaic virus (ZYMV, Potyvirus): vertical transmission, seed infection and cryptic infections

The role played by seed transmission in the evolution and epidemiology of viral crop pathogens remains unclear. We determined the seed infection and vertical transmission rates of zucchini yellow mosaic virus (ZYMV), in addition to undertaking Illumina sequencing of nine vertically transmitted ZYMV populations. We previously determined the seed-to-seedling transmission rate of ZYMV in Cucurbita pepo ssp. texana (a wild gourd) to be 1.6%, and herein observed a similar rate (1.8%) in the subsequent generation. We also observed that the seed infection rate is substantially higher (21.9%) than the seed-to-seedling transmission rate, suggesting that a major population bottleneck occurs during seed germination and seedling growth. In contrast, that two thirds of the variants present in the horizontally transmitted inoculant population were also present in the vertically transmitted populations implies that the bottleneck at vertical transmission may not be particularly severe. Strikingly, all of the vertically infected plants were symptomless in contrast to those infected horizontally, suggesting that vertical infection may be cryptic. Although no known virulence determining mutations were observed in the vertically infected samples, the 5' untranslated region was highly variable, with at least 26 different major haplotypes in this region compared to the two major haplotypes observed in the horizontally transmitted population. That the regions necessary for vector transmission are retained in the vertically infected populations, combined with the cryptic nature of vertical infection, suggests that seed transmission may be a significant contributor to the spread of ZYMV.

The evolutionary dynamics of influenza A virus adaptation to mammalian hosts

Few questions on infectious disease are more important than understanding how and why avian influenza A viruses successfully emerge in mammalian populations, yet little is known about the rate and nature of the virus' genetic adaptation in new hosts. Here, we measure, for the first time, the genomic rate of adaptive evolution of swine influenza viruses (SwIV) that originated in birds. By using a curated dataset of more than 24 000 human and swine influenza gene sequences, including 41 newly characterized genomes, we reconstructed the adaptive dynamics of three major SwIV lineages (Eurasian, EA; classical swine, CS; triple reassortant, TR). We found that, following the transfer of the EA lineage from birds to swine in the late 1970s, EA virus genes have undergone substantially faster adaptive evolution than those of the CS lineage, which had circulated among swine for decades. Further, the adaptation rates of the EA lineage antigenic haemagglutinin and neuraminidase genes were unexpectedly high and similar to those observed in human influenza A. We show that the successful establishment of avian influenza viruses in swine is associated with raised adaptive evolution across the entire genome for many years after zoonosis, reflecting the contribution of multiple mutations to the coordinated optimization of viral fitness in a new environment. This dynamics is replicated independently in the polymerase genes of the TR lineage, which established in swine following separate transmission from non-swine hosts.

Genes, trees and infections: Molecular evidence in epidemiology

Molecular phylogenies constitute an important new way of tracking the progress of viral epidemics. The phylogenetic analysis of viral sequence data provides information on the origin, spread and maintenance of infections and can be used to reconstruct contact networks of infected individuals. Analysis of the branching structure of phylogenetic trees also allows inferences to be made about the rate of transmission and the distinction between endemic and epidemic infections, and provides estimates of the numbers of infected individuals.

Deep sequencing reveals persistence of intra- and inter-host genetic diversity in natural and greenhouse populations of zucchini yellow mosaic virus

The genetic diversity present in populations of RNA viruses is likely to be strongly modulated by aspects of their life history, including mode of transmission. However, how transmission mode shapes patterns of intra- and inter-host genetic diversity, particularly when acting in combination with de novo mutation, population bottlenecks and the selection of advantageous mutations, is poorly understood. To address these issues, this study performed ultradeep sequencing of zucchini yellow mosaic virus in a wild gourd, Cucurbita pepo ssp. texana, under two infection conditions: aphid vectored and mechanically inoculated, achieving a mean coverage of approximately 10 ,000×. It was shown that mutations persisted during inter-host transmission events in both the aphid vectored and mechanically inoculated populations, suggesting that the vector-imposed transmission bottleneck is not as extreme as previously supposed. Similarly, mutations were found to persist within individual hosts, arguing against strong systemic bottlenecks. Strikingly, mutations were seen to go to fixation in the aphid-vectored plants, suggestive of a major fitness advantage, but remained at low frequency in the mechanically inoculated plants. Overall, this study highlights the utility of ultradeep sequencing in providing high-resolution data capable of revealing the nature of virus evolution, particularly as the full spectrum of genetic diversity within a population may not be uncovered without sequence coverage of at least 2500-fold.

Segmental configuration and putative origin of the reassortant orbivirus, epizootic hemorrhagic disease virus serotype 6, strain Indiana

In 2006, an exotic reassortant orbivirus, epizootic hemorrhagic disease virus serotype 6 (EHDV-6) [strain (Indiana)], was first detected in the United States. To characterize the reassortment configuration of this virus and to conclusively determine the parental virus of each RNA segment, the complete genome of EHDV-6 (Indiana) was sequenced, in addition to the genomes of representative EHDV-6 and EHDV-2 isolates. Based on genomic comparisons to all other EHDV serotypes, we determined that EHDV-6 (Indiana) originated from a reassortment event between the Australian prototype strain of EHDV-6 (CSIRO 753) and the North American topotype of EHDV-2 (Alberta). Additionally, phylogenetic analysis of all EHDV-6 (Indiana) isolates detected in the United States from 2006 to 2010 suggests that the virus may be undergoing continual reassortment with EHDV-2 (Alberta). In 2010, EHDV-6 (CSIRO 753) was detected in Guadeloupe, demonstrating that the parental virus of the reassortment event is circulating in the Caribbean.

Experimental Verification of Seed Transmission of Zucchini yellow mosaic virus

Within two decades of its discovery, Zucchini yellow mosaic virus (ZYMV) achieved a global distribution. However, whether or not seed transmission occurs in this economically significant crop pathogen is controversial, and the relative impact of seed transmission on the epidemiology of ZYMV remains unclear. Using reverse transcription-polymerase chain reaction, we observed a seed transmission rate of 1.6% in Cucurbita pepo subsp. texana and show that seed-infected C. pepo plants are capable of initiating horizontal ZYMV infections, both mechanically and via an aphid vector (Myzus persicae). We also provide evidence that ZYMV-infected seeds may act as effective viral reservoirs, partially accounting for the current geographic distribution of ZYMV. Finally, the observation that ZYMV infection of C. pepo seeds results in virtually symptomless infection, coupled with our finding that an antibody test failed to detect vertically transmitted ZYMV in infected seed, highlights the urgent need to standardize current detection methods for seed infection.

Epidemic dynamics revealed in dengue evolution

Dengue is an emerging tropical disease infecting tens of millions of people annually. A febrile illness with potentially severe hemorrhagic manifestations, dengue is caused by mosquito-borne viruses (DENV-1 to -4) that are maintained in endemic transmission in large urban centers of the tropics with periodic epidemic cycles at 3- to 5-year intervals. Puerto Rico (PR), a major population center in the Caribbean, has experienced increasingly severe epidemics since multiple dengue serotypes were introduced beginning in the late 1970s. We document the phylodynamics of DENV-4 between 1981 and 1998, a period of dramatic ecological expansion during which evolutionary change also occurs. The timescale of viral evolution is sufficiently short that viral transmission dynamics can be elucidated from genetic diversity data. Specifically, by combining virus sequence data with confirmed case counts in PR over these two decades, we show that the pattern of cyclic epidemics is strongly correlated with coalescent estimates of effective population size that have been estimated from sampled virus sequences using Bayesian Markov Chain Monte Carlo methods. Thus, we show that the observed epidemiologic dynamics are correlated with similar fluctuations in diversity, including severe interepidemic reductions in genetic diversity compatible with population bottlenecks that may greatly impact DENV evolutionary dynamics. Mean effective population sizes based on genetic data appear to increase prior to isolation counts, suggesting a potential bias in the latter and justifying more active surveillance of DENV activity. Our analysis explicitly integrates epidemiologic and sequence data in a joint model that could be used to further explore transmission models of infectious disease.

The evolutionary dynamics of canid and mongoose rabies virus in Southern Africa

Two variants of rabies virus (RABV) currently circulate in southern Africa: canid RABV, mainly associated with dogs, jackals, and bat-eared foxes, and mongoose RABV. To investigate the evolutionary dynamics of these variants, we performed coalescent-based analyses of the G-L inter-genic region, allowing for rate variation among viral lineages through the use of a relaxed molecular clock. This revealed that mongoose RABV is evolving more slowly than canid RABV, with mean evolutionary rates of 0.826 and 1.676 x 10(-3) nucleotide substitutions per site, per year, respectively. Additionally, mongoose RABV exhibits older genetic diversity than canid RABV, with common ancestors dating to 73 and 30 years, respectively, and while mongoose RABV has experienced exponential population growth over its evolutionary history in Africa, populations of canid RABV have maintained a constant size. Hence, despite circulating in the same geographic region, these two variants of RABV exhibit striking differences in evolutionary dynamics which are likely to reflect differences in their underlying ecology.

Evolution of the human immunodeficiency virus envelope gene is dominated by purifying selection

The evolution of the human immunodeficiency virus (HIV-1) during chronic infection involves the rapid, continuous turnover of genetic diversity. However, the role of natural selection, relative to random genetic drift, in governing this process is unclear. We tested a stochastic model of genetic drift using partial envelope sequences sampled longitudinally in 28 infected children. In each case the Bayesian posterior (empirical) distribution of coalescent genealogies was estimated using Markov chain Monte Carlo methods. Posterior predictive simulation was then used to generate a null distribution of genealogies assuming neutrality, with the null and empirical distributions compared using four genealogy-based summary statistics sensitive to nonneutral evolution. Because both null and empirical distributions were generated within a coalescent framework, we were able to explicitly account for the confounding influence of demography. From the distribution of corrected P-values across patients, we conclude that empirical genealogies are more asymmetric than expected if evolution is driven by mutation and genetic drift only, with an excess of low-frequency polymorphisms in the population. This indicates that although drift may still play an important role, natural selection has a strong influence on the evolution of HIV-1 envelope. A negative relationship between effective population size and substitution rate indicates that as the efficacy of selection increases, a smaller proportion of mutations approach fixation in the population. This suggests the presence of deleterious mutations. We therefore conclude that intrahost HIV-1 evolution in envelope is dominated by purifying selection against low-frequency deleterious mutations that do not reach fixation.

Cross-protective immunity can account for the alternating epidemic pattern of dengue virus serotypes circulating in Bangkok

Dengue virus, the causative agent of dengue fever and its more serious manifestation dengue hemorrhagic fever, is widespread throughout tropical and subtropical regions. The virus exists as four distinct serotypes, all of which have cocirculated in Bangkok for several decades with epidemic outbreaks occurring every 8-10 years. We analyze time-series data of monthly infection incidence, revealing a distinctive pattern with epidemics of serotypes 1, 2, and 3 occurring at approximately the same time and an isolated epidemic of serotype 4 occurring in the intervening years. Phylogenetic analysis of virus samples collected over the same period shows that clade replacement events are linked to the epidemic cycle and indicates that there is an interserotypic immune reaction. Using an epidemic model with stochastic seasonal forcing showing 8- to 10-year epidemic oscillations, we demonstrate that moderate cross-protective immunity gives rise to persistent out-of-phase oscillations similar to those observed in the data, but that strong or weak cross-protection or cross-enhancement only produces in-phase patterns. This behavior suggests that the epidemic pattern observed in Bangkok is the result of cross-protective immunity and may be significantly altered by changes in the interserotypic immune reaction.

Passive sexual transmission of human immunodeficiency virus type 1 variants and adaptation in new hosts

Human immunodeficiency virus type 1 (HIV-1) genetic diversity is a major obstacle for the design of a successful vaccine. Certain viral polymorphisms encode human leukocyte antigen (HLA)-associated immune escape, potentially overcoming limited vaccine protection. Although transmission of immune escape variants has been reported, the overall extent to which this phenomenon occurs in populations and the degree to which it contributes to HIV-1 viral evolution are unknown. Selection on the HIV-1 env gene at transmission favors neutralization-sensitive variants, but it is not known to what degree selection acts on the internal HIV-1 proteins to restrict or enhance the transmission of immune escape variants. Studies have suggested that HLA class I may determine susceptibility to HIV-1 infection, but a definitive role for HLA at transmission remains unproven. Comparing populations of acute seroconverters and chronically infected patients, we found no evidence of selection acting to restrict transmission of HIV-1 variants. We found that statistical associations previously reported in chronic infection between viral polymorphisms and HLA class I alleles are not present in acute infection, suggesting that the majority of viral polymorphisms in these patients are the result of transmission rather than de novo adaptation. Using four episodes of HIV-1 transmission in which the donors and recipients were both sampled very close to the time of infection we found that, despite a transmission bottleneck, genetic variants of HIV-1 infection are transmitted in a frequency-dependent manner. As HIV-1 infections are seeded by unique donor-adapted viral variants, each episode is a highly individual antigenic challenge. Host-specific, idiosyncratic HIV-1 antigenic diversity will seriously tax the efficacy of immunization based on consensus sequences.

Phylogenetic relationships among rhabdoviruses inferred using the L polymerase gene

RNA viruses of the family Rhabdoviridae include arthropod-borne agents that infect plants, fish and mammals, and also include a variety of non-vector-borne mammalian viruses. Herein is presented a molecular phylogenetic analysis, the largest undertaken to date, of 56 rhabdoviruses, including 20 viruses which are currently unassigned or assigned as tentative species within the Rhabdoviridae. Degenerate primers targeting a region of block III of the L polymerase gene were defined and used for RT-PCR amplification and sequencing. A maximum-likelihood phylogenetic analysis of a 158-residue L polymerase amino acid sequence produced an evolutionary tree containing the six recognized genera of the Rhabdoviridae and also enabled us to identify four more monophyletic groups of currently unclassified rhabdoviruses that we refer to as the 'Hart Park', 'Almpiwar', 'Le Dantec' and 'Tibrogargan' groups. The broad phylogenetic relationships among these groups and genera also indicate that the evolutionary history of rhabdoviruses was strongly influenced by mode of transmission, host species (plant, fish or mammal) and vector (orthopteran, homopteran or dipteran).

A phylogenetic survey of recombination frequency in plant RNA viruses

The severe economic consequences of emerging plant viruses highlights the importance of studies of plant virus evolution. One question of particular relevance is the extent to which the genomes of plant viruses are shaped by recombination. To this end we conducted a phylogenetic survey of recombination frequency in a wide range of positive-sense RNA plant viruses, utilizing 975 capsid gene sequences and 157 complete genome sequences. In total, 12 of the 36 RNA virus species analyzed showed evidence for recombination, comprising 17% of the capsid gene sequence alignments and 44% of the genome sequence alignments. Given the conservative nature of our analysis, we propose that recombination is a relatively common process in some plant RNA viruses, most notably the potyviruses.

A multigene analysis of the phylogenetic relationships among the flaviviruses (Family: Flaviviridae) and the evolution of vector transmission

The genus Flavivirus (family Flaviviridae) presently comprises around 70 single-strand positive-sense RNA viruses. These replicate in a range of vertebrate and invertebrate cells and may be mosquito-borne, tick-borne or have no-known-vector. Since transmission mode correlates strongly with phylogeny, the flaviviruses constitute a valuable model for the evolution of vector-borne disease. Attempts to resolve the higher-level taxonomic relationships of the flaviviruses through molecular phylogenetics have thus far proved inconclusive because of conflicting positions for the three main transmission groups. We conducted the most comprehensive phylogenetic study to date, involving maximum likelihood analyses of the NS3 and NS5 genes and the entire genome sequences available at present. For the first time, we use and test a variety of more robust methods of sequence alignment and appropriate models of amino acid replacement to study these highly divergent sequences, and explicitly test specific hypotheses of tree topology. We show that (i) the NS5 gene contains insufficient phylogenetic signal to choose between competing topological hypotheses, (ii) the NS3 gene and whole genome data indicate that the mosquito-borne flaviviruses represent an outgroup to the remaining flaviviruses, and (iii) that tick-borne transmission is probably a derived trait within the genus.

HIV evolution: CTL escape mutation and reversion after transmission

Within-patient HIV evolution reflects the strong selection pressure driving viral escape from cytotoxic T-lymphocyte (CTL) recognition. Whether this intrapatient accumulation of escape mutations translates into HIV evolution at the population level has not been evaluated. We studied over 300 patients drawn from the B- and C-clade epidemics, focusing on human leukocyte antigen (HLA) alleles HLA-B57 and HLA-B5801, which are associated with long-term HIV control and are therefore likely to exert strong selection pressure on the virus. The CTL response dominating acute infection in HLA-B57/5801-positive subjects drove positive selection of an escape mutation that reverted to wild-type after transmission to HLA-B57/5801-negative individuals. A second escape mutation within the epitope, by contrast, was maintained after transmission. These data show that the process of accumulation of escape mutations within HIV is not inevitable. Complex epitope- and residue-specific selection forces, including CTL-mediated positive selection pressure and virus-mediated purifying selection, operate in tandem to shape HIV evolution at the population level.

Geographic and species association of hepatitis B virus genotypes in non-human primates

Infection with hepatitis B virus (HBV) has been detected in human populations throughout the world, as well as in a number of ape species (Pan troglodytes, Gorilla gorilla, gibbons [Nomascus and Hylobates species] and Pongo pygmaeus). To investigate the distribution of naturally occurring HBV infection in these species and other African Old World monkey species (Cercopithecidae), we screened 137 plasma samples from mainly wild caught animals by polymerase chain reaction (PCR) using several of highly conserved primers from the HB surface (HBs) gene, and for HBs antigen (HBsAg) by ELISA. None of the 93 Cercopithecidae screened (6 species) showed PCR or serology evidence for HBV infection; in contrast 2 from 8 chimpanzees and 5 from 22 gibbons were PCR-positive with each set of primers. Complete genome sequences from each of the positive apes were obtained and compared with all previously published complete and surface gene sequences. This extended phylogenetic analysis indicated that HBV variants from orangutans were interspersed by with HBV variants from southerly distributed gibbon species (H. agilis and H. moloch) occupying overlapping or adjacent habitat ranges with orangutans; in contrast, HBV variants from gibbon species in mainland Asia were phylogenetically distinct. A geographical rather than (sub)species association of HBV would account for the distribution of HBV variants in different subspecies of chimpanzees in Africa, and explain the inlier position of the previously described lowland gorilla sequence in the chimpanzee clade. These new findings have a number of implication for understanding the origins and epidemiology of HBV infection in non-human primates.

Molecular epidemiology of dengue virus type 3 in Venezuela

During the past 40 years, dengue haemorrhagic fever and dengue shock syndrome (DHF/DSS) have emerged in humans, with approximately 3 million cases reported and over 58 000 deaths. Dengue virus serotypes 1, 2 and 4 (DENV-1, -2 and -4) have been co-circulating in Venezuela for at least the past 10 years, causing minor or major outbreaks of dengue fever (DF) and DHF/DSS. The first recorded outbreak due to DENV-3 in Venezuela dates to 1964 and the virus then seems to have disappeared. However, DENV-3 re-appeared recently (in July, 2000) in Venezuela after 32 years of absence and produced a prolonged major outbreak, which, by the end of 2001, involved 83 180 cases of dengue, mostly DF (92 %). Previous phylogenetic studies revealed that the DENV-3 circulating during the 1960s Latin American outbreak was a genotype V virus. To gain a better understanding of the nature of the current epidemic, the complete sequence was determined of the envelope (E) gene of 15 Venezuelan DENV-3 viruses isolated during 2000 and 2001 from patients presenting with different disease severity. Sequence data were used in phylogenetic comparisons with global samples of DENV-3. Analysis revealed that the strain circulating in Venezuela is closely related to isolates that were previously present in Panama and Nicaragua in 1994 and since then have spread through Central American countries and Mexico. This study also confirms previous reports showing that the DENV-3 strain currently circulating in the Americas is related to the strain that caused DHF epidemics in Sri Lanka and India in 1989-1991 (genotype III). Finally, no evidence of the re-emergence of the strain that circulated in Venezuela in the late 1960s and 1970s (genotype V) was found.

New inferences from tree shape: numbers of missing taxa and population growth rates

The relative positions of branching events in a phylogeny contain information about evolutionary and population dynamic processes. We provide new summary statistics of branching event times and describe how these statistics can be used to infer rates of species diversification from interspecies trees or rates of population growth from intraspecies trees. We also introduce a phylogenetic method for estimating the level of taxon sampling in a clade. Different evolutionary models and different sampling regimes can produce similar patterns of branching events, so it is important to consider explicitly the model assumptions involved when making evolutionary inferences. Results of an analysis of the phylogeny of the mosquito-borne flaviviruses suggest that there could be several thousand currently unidentified viruses in this clade.

Genome sequence analysis of Tamana bat virus and its relationship with the genus Flavivirus

Tamana bat virus (TABV, isolated from the bat Pteronotus parnellii) is currently classified as a tentative species in the genus FLAVIVIRUS: We report here the determination and analysis of its complete coding sequence. Low but significant similarity scores between TABV and member-viruses of the genus Flavivirus were identified in the amino acid sequences of the structural, NS3 and NS5 genes. A series of cysteines located in the envelope protein and the most important enzymatic domains of the virus helicase/NTPase, methyltransferase and RNA-dependent RNA polymerase were found to be highly conserved. In the serine-protease domain, the catalytic sites were conserved, but variations in sequence were found in the putative substrate-binding sites, implying possible differences in the protease specificity. In accordance with this finding, the putative cleavage sites of the TABV polyprotein by the virus protease are substantially different from those of flaviviruses. The phylogenetic position of TABV could not be determined precisely, probably due to the extremely significant genetic divergence from other member-viruses of the family FLAVIVIRIDAE: However, analysis based on both genetic distances and maximum-likelihood confirmed that TABV is more closely related to the flaviviruses than to the other genera. These findings have implications for the evolutionary history and taxonomic classification of the family as a whole: (i) the possibility that flaviviruses were derived from viruses infecting mammals rather than from mosquito viruses cannot be excluded; (ii) using the current criteria for the definition of genera in the family Flaviviridae, TABV should be assigned to a new genus.

Extinction and rapid emergence of strains of dengue 3 virus during an interepidemic period

Strains of dengue 3 (DEN-3) virus circulating in Thailand prior to 1992 appear to have disappeared from that location and to have been replaced by two new lineages which have evolved locally, rather than being introduced. Similar DEN-3 virus extinctions may have occurred previously in Thailand in 1962 and 1973. Although no causal relationship could be shown, this strain replacement event was accompanied by DEN-3 replacing DEN-2 as the serotype recovered most frequently from patients in Thailand. Although this implies a change in selection pressure, we found no evidence for positive natural selection at the level of either the E protein or the E protein gene. Further, the extinction of the pre-1992 strains and the appearance of the new lineages occurred during an interepidemic period, suggesting that a genetic bottleneck, rather than selection, might have been important in the emergence of these two new strains of virus. The pre-1992 DEN-3 virus lineage could still be found in 1998, to the west, in Myanmar. The ratio of nonsynonymous-to-synonymous nucleotide changes within a DEN-3 virus population from a single patient was less than the ratio among the consensus sequences of DEN-3 viruses from different patients, suggesting that many of the nonsynonymous nucleotide changes which occurred naturally in the E protein were deleterious and removed by purifying selection.

Molecular epidemiology of dengue type 2 virus in Venezuela: evidence for in situ virus evolution and recombination

Epidemic outbreaks of dengue fever (DF) were first recorded in Venezuela in 1978 and were followed by the emergence of dengue haemorrhagic fever (DHF) outbreaks in 1989. To gain a better understanding of the nature of these epidemics, the complete envelope (E) gene sequence of 34 Venezuelan dengue type 2 (DEN-2) viruses, isolated between 1997 and 2000 was determined. Of these isolates, 16 were from patients with DF and 17 were from patients diagnosed with DHF. There were no diagnostic sequence differences between them, suggesting that the E gene alone does not determine disease severity. These sequence data were also used in phylogenetic comparisons with a global sample of DEN-2 viruses, including strains collected previously from Venezuela. This analysis revealed that the ancestors of the Venezuelan viruses were Asian in origin, implying that a DEN-2 virus strain from this region was introduced into Venezuela and the wider Caribbean region during the late 1970s or the early 1980s. The phylogenetic trees further indicate that evolution of DEN-2 virus in Venezuela has occurred in situ, with differentiation into a number of distinct but co-circulating lineages, rather than the repeated introduction of new strains from other localities. By incorporating additional sequence data from the virus capsid, premembrane and membrane genes, evidence is provided that a single Venezuelan strain sequenced previously, designated Mara4, is a recombinant virus, incorporating genome sequence from Venezuelan and Asian parental viruses.

Mother-to-child transmission of HIV infection and CTL escape through HLA-A2-SLYNTVATL epitope sequence variation

Cytotoxic T lymphocytes (CTL) play a central role in containment of HIV infection. Evasion of the immune response by CTL escape is associated with progression to disease. It is therefore hypothesised that transmitted viruses encode escape mutations within epitopes that are required for successful control of viraemia. In order to test this hypothesis, escape through the dominant HLA-A2-restricted CTL epitope SLYNTVATL (p17 Gag residues 77-85 SL9) in the setting of mother-to-child-transmission (MTCT) was investigated. Initial data from two families in which the HIV-infected mother expressed HLA-A*0201 and had transmitted the virus to other family members were consistent with this hypothesis. In addition, analysis of the gag sequence phylogeny in one family demonstrated that CTL escape variants can be successfully transmitted both horizontally and vertically. To test the hypothesis further, a larger cohort of transmitting mothers (n=8) and non-transmitters (n=14) were studied. Variation within the SL9 epitope was associated with expression of HLA-A2 (P=0.04) but overall no clear link between variation from the SL9 consensus sequence and MTCT was established. However, the high level of background diversity within p17 Gag served to obscure any possible association between escape and MTCT. In conclusion, these studies highlighted the obstacles to demonstrating CTL escape arising at this particular epitope. Alternative strategies likely to be more definitive are discussed.

Transmission of hepatitis B virus infection in Gambian families revealed by phylogenetic analysis

BACKGROUND/AIMS

Transmission of hepatitis B virus (HBV) in Africa occurs horizontally, with most people becoming infected between the ages of 1 and 5 years. The index cases in such events have been assumed to come from within the family unit or from sources outside the immediate family, such as other families or inhabitants of the same compound or village. Here, we define these routes of transmission by phylogenetic tree analysis of sequences from the entire pre-core/core region of the virus, in Gambian chronic carriers.

METHODS

Amplification by polymerase chain reaction of serum extracted HBV-DNA was followed by direct sequencing of the target region. Following editing and alignment of these sequences, phylogenetic tree analysis was performed using the neighbour-joining and maximum-likelihood methods.

RESULTS

Despite the overall conserved nature of the sequences of the pre-core/core region from 142 chronic carriers, distinct clusters were easily defined at the family and village level, but not on a wider geographical separation.

CONCLUSIONS

Phylogenetic tree analysis of sequences obtained from family members provided strong evidence of intrafamilial transmission of HBV in at least two-thirds of the families studied from Gambia.

The epidemic behavior of the hepatitis C virus

Hepatitis C virus (HCV) is a leading worldwide cause of liver disease. Here, we use a new model of HCV spread to investigate the epidemic behavior of the virus and to estimate its basic reproductive number from gene sequence data. We find significant differences in epidemic behavior among HCV subtypes and suggest that these differences are largely the result of subtype-specific transmission patterns. Our model builds a bridge between the disciplines of population genetics and mathematical epidemiology by using pathogen gene sequences to infer the population dynamic history of an infectious disease.

Evidence for the Non-quasispecies Evolution of RNA Viruses

The quasispecies model of RNA virus evolution differs from those formulated in conventional population genetics in that neutral mutations do not lead to genetic drift of the population, and natural selection acts on the mutant distribution as a whole rather than on individual variants. By computer simulation, we show that this model could be inappropriate for many RNA viruses because the neutral sequence space may be too large to allow the formation of a quasispecies distribution. This view is supported by our analysis of gene sequences from vesicular stomatitis virus, which is considered a prototype RNA virus quasispecies. Our results are relevant to the evolution of RNA systems in general.

Phylogenetic analysis of multiple heterosexual transmission events involving subtype b of HIV type 1

Between 1996 and 1999 thirteen cases of HIV infection were detected in Doncaster, a small town in the north of England (population approximately 250,000). A complex network of shared sexual histories involving local nightclubs linked these cases, with the only known risk factor being heterosexual intercourse. A series of frozen blood samples was collected in 1998-1999 and amplified by PCR to generate full-length gp120 clones. Sequencing demonstrated that all the transmission events in this heterosexual group involved the B subtype of HIV-1. When relationships between the samples were assessed it became clear that these 13 cases represented at least three separate strains of HIV-1, indicating that HIV is well established in this community. Eleven of the 13 cases were related, forming two distinct groups. Further investigation revealed that one group contained five patients whose general health was good and who were not receiving HAART. In contrast, the second group of six patients, including the putative index case, were symptomatic, receiving HAART, and may have been infected with a CXCR-4-utilizing virus. Several of the cases that were linked by genetic criteria were not linked by contact tracing, implying that further undiagnosed cases may exist in this community. To our knowledge, this is the largest outbreak of HIV studied within the heterosexual community in the United Kingdom to date, suggesting that this route of infection is becoming more common within the United Kingdom.

On the origin and evolution of the human immunodeficiency virus (HIV)

The human AIDS viruses--HIV-1 and HIV-2--impose major burdens on the health and economic status of many developing countries. Surveys of other animal species have revealed that related viruses--the SIVs are widespread in a large number of African simian primates where they do not appear to cause disease. Phylogenetic analyses indicate that these SIVs are the reservoirs for the human viruses, with SIVsm from the sooty mangabey monkey the most likely source of HIV-2, and SIVcpz from the common chimpanzee the progenitor population for HIV-1. Although it is clear that AIDS has a zoonotic origin, it is less certain when HIV-1 and HIV-2 first entered human populations and whether cross-species viral transmission is common among primates. Within infected individuals the process of HIV evolution takes the form of an arms race, with the virus continually fixing mutations by natural selection which allow it to escape from host immune responses. The arms race is less intense in SIV-infected monkeys, where a weaker immune response generates less selective pressure on the virus. Such a difference in virus-host interaction, along with a broadening of co-receptor usage such that HIV strains are able to infect cells with both CCR5 and CXCR4 chemokine receptors, may explain the increased virulence of HIV in humans compared to SIV in other primates.

Evolution of base composition and codon usage bias in the genus Flavivirus

The extent to which base composition and codon usage vary among RNA viruses, and the possible causes of this bias, is undetermined in most cases. A maximum-likelihood statistical method was used to test whether base composition and codon usage bias covary with arthropod association in the genus Flavivirus, a major source of disease in humans and animals. Flaviviruses are transmitted by mosquitoes, by ticks, or directly between vertebrate hosts. Those viruses associated with ticks were found to have a significantly lower G+C content than non-vector-borne flaviviruses and this difference was present throughout the genome at all amino acids and codon positions. In contrast, mosquito-borne viruses had an intermediate G+C content which was not significantly different from those of the other two groups. In addition, biases in dinucleotide and codon usage that were independent of base composition were detected in all flaviviruses, but these did not covary with arthropod association. However, the overall effect of these biases was slight, suggesting only weak selection at synonymous sites. A preliminary analysis of base composition, codon usage, and vector specificity in other RNA virus families also revealed a possible association between base composition and vector specificity, although with biases different from those seen in the Flavivirus genus.

The CD45 locus in cattle: allelic polymorphism and evidence for exceptional positive natural selection

Cattle in Africa are a genetically diverse population that has resulted from successive introduction of Asian Bos indicus and European B. taurus cattle. However, analysis of mitochondrial genetic diversity in African cattle identified three lineages, one associated with Asian B. indicus, one with European B. taurus, and a third ascribed to an indigenous African sub-species of cattle. Due to their extended coevolution, indigenous African herbivores are generally tolerant to endemic African pathogens. We are interested in identifying alleles derived from the indigenous African cattle that may be associated with tolerance to African pathogens. An analysis of the locus which encodes the abundant plasma membrane-associated tyrosine phosphatase, CD45, identified three highly divergent allelic families in Kenya Boran cattle. Analysis of allelic distribution in a diverse range of cattle populations suggests a European B. taurus, an Asian B. indicus, and an African origin. This demonstrates not only significant allelic polymorphism at the CD45 locus in cattle but also convincing autosomal evidence for a distinct African sub-species of cattle. Furthermore, maximum-likelihood analysis of selection pressures revealed that the CD45 locus is subject to exceptionally strong natural selection which we suggest may be pathogen driven.

Management of acute cardiac failure with mechanical assist: experience with the ABIOMED BVS 5000

BACKGROUND

Mechanical circulatory assist industries have developed ventricular assist devices (VAD) for short-, intermediate-, and long-term use. The purpose of this report is to describe the progress made with the ABIOMED Biventricular System (BVS) 5000 (ABIOMED, Inc, Danvers, MA) short-term VAD.

METHODS

From June 1994 through August 2000, all cardiogenic shock patients who required short-term mechanical assist were supported with the ABIOMED BVS 5000. Insertion criteria included any condition that may potentially result in cardiac recovery. A formal algorithm for timing of insertion was established to standardize implantation criteria.

RESULTS

A total of 45 patients were supported at Hahnemann University Hospital, Philadelphia, PA. There were 26 male and 19 female patients, with a mean age of 57.9 years (range 33 to 80 years). Devices were inserted for postcardiotomy shock in 36 patients (80%) and precardiotomy shock in 9 patients (20%). The average duration of support was 8.3 days (range 1 to 31 days). Overall, there were 22 (49%) patients weaned from support and 14 (31%) discharged from the hospital. For patients in whom the device was implanted in accordance with an established protocol (group A), the wean and discharge rates were 60% and 43%, respectively. The most common morbidities included bleeding and adverse neurologic events.

CONCLUSIONS

The ABIOMED BVS 5000 VAD continues to be a valuable form of short-term mechanical assist for acute cardiogenic shock. The formation of a uniform VAD insertion algorithm has helped to standardize protocols in management.

Clustered mutations in HIV-1 gag are consistently required for escape from HLA-B27-restricted cytotoxic T lymphocyte responses

The immune response to HIV-1 in patients who carry human histocompatibility leukocyte antigen (HLA)-B27 is characterized by an immunodominant response to an epitope in p24 gag (amino acids 263-272, KRWIILGLNK). Substitution of lysine (K) or glycine (G) for arginine (R) at HIV-1 gag residue 264 (R264K and R264G) results in epitopes that bind to HLA-B27 poorly. We have detected a R264K mutation in four patients carrying HLA-B27. In three of these patients the mutation occurred late, coinciding with disease progression. In another it occurred within 1 yr of infection and was associated with a virus of syncytium-inducing phenotype. In each case, R264K was tightly associated with a leucine to methionine change at residue 268. After the loss of the cytotoxic T lymphocyte (CTL) response to this epitope and in the presence of high viral load, reversion to wild-type sequence was observed. In a fifth patient, a R264G mutation was detected when HIV-1 disease progressed. Its occurrence was associated with a glutamic acid to aspartic acid mutation at residue 260. Phylogenetic analyses indicated that these substitutions emerged under natural selection rather than by genetic drift or linkage. Outgrowth of CTL escape viruses required high viral loads and additional, possibly compensatory, mutations in the gag protein.

Homologous recombination in GB virus C/hepatitis G virus

Analysis of 33 GB virus C/hepatitis G virus (GBV-C/HGV) full or nearly full genome sequences revealed several putative inter- and intrasubtype recombinants. The breakpoints of the recombinant regions were mapped using a maximum-likelihood method, and the statistical significance for each region was tested using Monte Carlo simulation. The results were highly significant and provided evidence for the existence of complex mosaic genomes showing as many as nine recombination events, with breakpoints in the 5' UTR and in all of the coding regions except the short NS4b gene. Recombination was confirmed by separate phylogenetic analysis of the various recombinant regions and by Sawyer's runs test. Taken together, these findings demonstrate for the first time that recombination is common in natural populations of GBV-C and that it takes place both within and between subtypes. The wide-ranging implications of such nonclonal history for reconstructing the spread and timescale of GBV-C evolution are discussed.

Variable immune-driven natural selection in the attachment (G) glycoprotein of respiratory syncytial virus (RSV)

A maximum-likelihood analysis of selection pressures acting on the attachment (G) glycoprotein gene of respiratory syncytial virus (RSV) from humans (HRSV) and bovines (BRSV) is presented. Six positively selected sites were identified in both group A and group B of HRSV, although only one site was common between them, while no positively selected sites were detected in BRSV. All positively selected sites were located within the ectodomain of the G protein and showed some association with positions of immunoglobulin (Ig) epitopes and sites of O-glycosylation. These results suggest that immune (antibody)-driven natural selection is an important determinant of RSV evolution and that this selection pressure differs among strains. The passage histories of RSV strains were also shown to affect the distribution of positively selected sites, particularly in HRSV B, and should be considered whenever retrospective analysis of adaptive evolution is undertaken.

Recombination within natural populations of pathogenic bacteria: short-term empirical estimates and long-term phylogenetic consequences

The identification of clones within bacterial populations is often taken as evidence for a low rate of recombination, but the validity of this inference is rarely examined. We have used statistical tests of congruence between gene trees to examine the extent and significance of recombination in six bacterial pathogens. For Neisseria meningitidis, Streptococcus pneumoniae, Streptococcus pyogenes, and Staphylococcus aureus, the congruence between the maximum likelihood trees reconstructed using seven house-keeping genes was in most cases no better than that between each tree and trees of random topology. The lack of congruence between gene trees in these four species, which include both naturally transformable and nontransformable species, is in three cases supported by high ratios of recombination to point mutation during clonal diversification (estimates of this parameter were not possible for Strep. pyogenes). In contrast, gene trees constructed for Hemophilus influenzae and pathogenic isolates of Escherichia coli showed a higher degree of congruence, suggesting lower rates of recombination. The impact of recombination therefore varies between bacterial species but in many species is sufficient to obliterate the phylogenetic signal in gene trees.

Endoscopic ultrasound does not accurately assess pathologic stage of esophageal cancer after neoadjuvant chemoradiotherapy

The accuracy of endoscopic ultrasound (EUS) for initial staging of esophageal cancer is widely accepted. There is, however, considerable variability in the reported accuracy of EUS for restaging of esophageal neoplasms after neoadjuvant therapy. From June 1995 through December 1999, we prospectively studied a series of 26 patients who underwent neoadjuvant treatment for esophageal cancer and were subsequently restaged by EUS before resection. Twenty-four patients had adenocarcinoma (92%), and two patients had squamous cell carcinoma (8%). EUS correctly predicted tumor stage in seven of 26 patients for an overall accuracy of 27 per cent. EUS overestimated the depth of tumor penetration in 18 patients (69%) and underestimated depth of penetration in one patient (4%). Lymph nodes were correctly staged in 15 of 26 patients for an overall accuracy of 58 per cent. Levels of sensitivity for detecting N0 and N1 disease were 44 per cent and 80 per cent respectively. Patients with a complete pathologic response were staged as T4N1 (one patient), T3N1 (three patients), T3N0 (one patient), and T2N1 (two patients). EUS cannot distinguish tumor involvement of the esophageal wall and lymph nodes from the postinflammatory changes that characterize effective neoadjuvant treatment. EUS is of limited utility in guiding clinical decision making after neoadjuvant therapy.

The evolution of base composition and phylogenetic inference

Base composition varies at all levels of the phylogenetic hierarchy and throughout the genome, and can be caused by active selection or passive mutation pressure. This variation can make reconstructing trees difficult. However, recent tree-based analyses have shed light on the forces responsible for the evolution of base composition, forces that might be very general. More explicit tree-based work is encouraged.

Genetic diversity of Borrelia burgdorferi sensu lato in ticks from mainland Portugal

To date Borrelia lusitaniae is the only genospecies of Borrelia burgdorferi sensu lato isolated from Ixodes ricinus ticks collected in Portugal and Tunisia. This suggests that the genospecies diversity of B. burgdorferi sensu lato decreases toward the southwestern margin of its Old World subtropical range. In order to further explore the genetic diversity of B. burgdorferi sensu lato from this region, 55 I. ricinus and 27 Hyalomma marginatum questing adults, collected during the spring of 1998 from a sylvatic habitat south of Lisbon, Portugal, were analyzed. Infection prevalences of 75% in I. ricinus ticks and 7% in H. marginatum ticks were detected by a nested PCR that targets the rrf (5S)-rrl (23S) spacer of B. burgdorferi sensu lato. Restriction fragment length polymorphism (RFLP) analysis of the I. ricinus-derived amplicons showed that the sequences in the majority of samples were similar to those of B. lusitaniae type strains (76% for strain PotiB1, 5% for strain PotiB3). Two novel RFLP patterns were obtained from 12% of the samples. The remaining 7% of samples gave mixed RFLP patterns. Phylogenetic analysis of rrf-rrl spacer sequences revealed a diverse population of B. lusitaniae in questing adult I. ricinus ticks (the sequences did not cluster with those of any other genospecies). This population consisted of 10 distinct sequence types, suggesting that multiple strains of B. lusitaniae were present in the local I. ricinus population. We hypothesize that B. lusitaniae has a narrow ecological niche that involves host species restricted to the Mediterranean Basin.

Detection of hepatitis B virus infection in wild-born chimpanzees (Pan troglodytes verus): phylogenetic relationships with human and other primate genotypes

Infection with hepatitis B virus (HBV) was detected by serological testing for HBV surface antigen and by PCR assay for HBV DNA in serum samples from two common chimpanzees (Pan troglodytes subsp. verus) born in West Africa. The complete genome sequences obtained by nucleotide sequencing of overlapping DNA fragments amplified by PCR were compared with HBV variants recovered from other primates and with human genotypes A to F. Both chimpanzee sequences were 3, 182 nucleotides in length, and the surface gene sequence predicted the existence of a, d, and w serological determinants. Neither sequence contained stop codons in the precore region. On phylogenetic analysis, the HBV variants infecting the chimpanzees clustered together with a third chimpanzee HBV isolate independently obtained from an infected captive animal (A. J. Zuckerman, A. Thornton, C. R. Howard, K. N. Tsiquaye, D. M. Jones, and M. R. Brambell, Lancet ii:652-654, 1978), with an overall sequence similarity of >94%. This provides strong evidence for a chimpanzee-specific genotype of HBV which circulates in nature. These findings add to the recent evidence for infection in the wild of other Old and New World primates (gibbon, orangutan, and woolly monkey) with species-specific variants of HBV. There is no evidence for close phylogenetic clustering of variants found so far in primates with any of the established HBV genotypes from humans. With the new evidence for the widespread distribution of HBV in primates, hypotheses for the origins of human infection are reviewed.

The causes and consequences of genetic variation in dengue virus

Despite the fact that dengue is one of the most prevalent viral infections of humans, the mechanisms responsible for its pathogenesis remain uncertain. Evolutionary studies of dengue virus have revealed that its genetic diversity is increasing. This, coupled with evidence that viral strains could naturally differ in virulence, suggests that in the future we might be exposed to viruses with an expanded range of pathogenic properties.