An exact nonparametric method for inferring mosaic structure in sequence triplets

Reconstructing the history of maize streak virus strain a dispersal to reveal diversification hot spots and its origin in southern Africa

Maize streak virus strain A (MSV-A), the causal agent of maize streak disease, is today one of the most serious biotic threats to African food security. Determining where MSV-A originated and how it spread transcontinentally could yield valuable insights into its historical emergence as a crop pathogen. Similarly, determining where the major extant MSV-A lineages arose could identify geographical hot spots of MSV evolution. Here, we use model-based phylogeographic analyses of 353 fully sequenced MSV-A isolates to reconstruct a plausible history of MSV-A movements over the past 150 years. We show that since the probable emergence of MSV-A in southern Africa around 1863, the virus spread transcontinentally at an average rate of 32.5 km/year (95% highest probability density interval, 15.6 to 51.6 km/year). Using distinctive patterns of nucleotide variation caused by 20 unique intra-MSV-A recombination events, we tentatively classified the MSV-A isolates into 24 easily discernible lineages. Despite many of these lineages displaying distinct geographical distributions, it is apparent that almost all have emerged within the past 4 decades from either southern or east-central Africa. Collectively, our results suggest that regular analysis of MSV-A genomes within these diversification hot spots could be used to monitor the emergence of future MSV-A lineages that could affect maize cultivation in Africa.

Analysis of the complete DNA sequence of rice tungro bacilliform virus from southern India indicates it to be a product of recombination

The complete nucleotide sequence of an isolate of rice tungro bacilliform virus (RTBV), collected from Kanyakumari, India, where RTBV was reported recently for the first time, has been analyzed. Sequence comparison revealed that the RTBV isolate from Kanyakumari (RTBV-KK) has a high degree of identity to the two previously reported RTBV sequences from India, RTBV-AP and RTBV-WB, which had been collected from field locations about 10 years ago and 1000-2000 km away from the collection site of RTBV-KK. Most of the sequence domains reported previously in other RTBV isolates were found to be conserved in RTBV-KK. Closer inspection revealed RTBV-KK to be a possible recombinant between RTBV-AP and RTBV-WB in the genomic region encompassing the coat protein gene.

Molecular characterization of the complete genomes of two new field isolates of Cowpea chlorotic mottle virus, and their phylogenetic analysis

Cowpea chlorotic mottle virus (CCMV, family Bromoviridae) is found worldwide and has been used as a model virus for a long time, but no data is available about the genetic diversity of field isolates. Recently, two new field isolates (Car1 and Car2) of CCMV obtained from cowpea showed distinct phenotypic symptoms when inoculated to cowpea. CCMV-Car1 induced severe mosaic and interveinal chlorosis, while CCMV-Car2 produced mild mottling and leaf rolling. Both isolates produced asymptomatic infection in Nicotiana benthamiana. The complete genome of both isolates was amplified by reverse transcription-polymerase chain reaction using specific primers against the CCMV sequences available in the GenBank database, cloned and sequenced. Both nucleotide and amino acid sequences were compared between the newly sequenced CCMV isolates and the three previously characterized CCMV strains (T, M1, and R). Phylogenetic analysis of the RNA 1 sequence showed that CCMV-Car1 was in a separate branch from the rest of the CCMV isolates while CCMV-Car2 grouped together with CCMV-R. On the basis of RNA 2 and RNA 3 sequences, two major groupings were obtained. One group included CCMV-Car1 and CCMV-Car2 isolates while the other contained CCMV-T, CCMV-M1, and CCMV-R strains. Recombination programs detected a potential recombination event in the RNA 1 sequence of CCMV-Car2 isolate but not in RNA 2 and RNA 3 sequences. The results showed that both mutations and recombination have played an important role in the genetic diversity of these two new isolates of CCMV.

Conversion events in gene clusters

BACKGROUND

Gene clusters containing multiple similar genomic regions in close proximity are of great interest for biomedical studies because of their associations with inherited diseases. However, such regions are difficult to analyze due to their structural complexity and their complicated evolutionary histories, reflecting a variety of large-scale mutational events. In particular, conversion events can mislead inferences about the relationships among these regions, as traced by traditional methods such as construction of phylogenetic trees or multi-species alignments.

RESULTS

To correct the distorted information generated by such methods, we have developed an automated pipeline called CHAP (Cluster History Analysis Package) for detecting conversion events. We used this pipeline to analyze the conversion events that affected two well-studied gene clusters (α-globin and β-globin) and three gene clusters for which comparative sequence data were generated from seven primate species: CCL (chemokine ligand), IFN (interferon), and CYP2abf (part of cytochrome P450 family 2). CHAP is freely available at http://www.bx.psu.edu/miller_lab.

CONCLUSIONS

These studies reveal the value of characterizing conversion events in the context of studying gene clusters in complex genomes.

Artificial recombination may influence the evolutionary analysis of Newcastle disease virus

The recombination rate in Newcastle disease virus (NDV) was as high as 10% in RDP analysis with full-length NDV genome sequences available in GenBank. We found that two NDV strains, China/Guangxi09/2003 and NDV/03/018, previously reported as recombinants, failed to show any evidence of recombination upon complete genome resequencing. Furthermore, we were able to reproduce artificial recombination by amplification of the M gene in a mixed sample of strains LaSota and ZJ1. It appears that the recombination of NDV is not as common as has been reported. NDV sequences in GenBank should be analyzed with caution during bioinformatic analyses for natural recombination events.

Multi-locus sequence typing of Ehrlichia ruminantium strains from geographically diverse origins and collected in Amblyomma variegatum from Uganda

Background

The rickettsial bacterium Ehrlichia ruminantium is the causative agent of heartwater in ruminants. A better understanding of the population genetics of its different strains is, however, needed for the development of novel diagnostic tools, therapeutics and prevention strategies. Specifically, the development of effective vaccination policies relies on the proper genotyping and characterisation of field isolates. Although multi-locus sequence typing (MLST) has been recently developed, only strains from geographically restricted collections have been analysed so far. The expansion of the MLST database to include global strains with different geographic origins is therefore essential. In this study, we used a panel of reference strains from geographically diverse origins and field samples of E. ruminantium detected from its vector, Amblyomma variegatum, in heartwater-endemic areas in Uganda.

Results

A total of 31 novel alleles (six, four, six, three, two, five, three, and two for gltA, groEL, lepA, lipA, lipB, secY, sodB, and sucA loci, respectively) and 19 novel sequence types (STs) were identified. Both neighbour-joining and minimum spanning tree analyses indicated a high degree of genetic heterogeneity among these strains. No association was observed between genotypes and geographic origins, except for four STs from West African countries. When we performed six different tests for recombination (GeneConv, Bootscan, MaxChi, Chimaera, SiScan, and 3Seq) on concatenated sequences, four possible recombination events were identified in six different STs. All the recombination breakpoints were located near gene borders, indicating the occurrence of intergenic recombination. All four STs that localized to a distinct group in clustering analysis showed evidence of identical recombination events, suggesting that recombination may play a significant role in the diversification of E. ruminantium.

Conclusions

The compilation of MLST data set across the African continent will be particularly valuable for the understanding of the existing genetic diversity of field isolates in African countries. Comprehensive information on the degree of cross-protection between strains and further understanding of possible relationships between genotypes and phenotypes such as vaccine efficacy are expected to lead to the development of region-specific vaccination strategies.

DNA barcoding of lichenized fungi demonstrates high identification success in a floristic context

• Efforts are currently underway to establish a standard DNA barcode region for fungi; we tested the utility of the internal transcribed spacer (ITS) of nuclear ribosomal DNA for DNA barcoding in lichen-forming fungi by sampling diverse species across eight orders. • Amplification of the ITS region (ITS1-5.8S-ITS2) was conducted for 351 samples, encompassing 107, 55 and 28 species, genera and families, respectively, of lichenized fungi. We assessed the ability of the entire ITS vs the ITS2 alone to discriminate between species in a taxonomic dataset (members of the genus Usnea) and a floristic dataset. • In the floristic dataset, 96.3% of sequenced samples could be assigned to the correct species using ITS or ITS2; a barcode gap for ITS is present in 92.1% of species. Although fewer species have a barcode gap in the taxonomic dataset (73.3% with ITS and 68.8% with ITS2), up to 94.1% of samples were assigned to the correct species using BLAST. • While discrimination between the most closely related species will remain challenging, our results demonstrate the potential to identify a high percentage of specimens to the correct species, and the remainder to the correct genus, when using DNA barcoding in a floristic context.

Multilocus sequence typing (MLST) for lineage assignment and high resolution diversity studies in Trypanosoma cruzi

BACKGROUND

Multilocus sequence typing (MLST) is a powerful and highly discriminatory method for analysing pathogen population structure and epidemiology. Trypanosoma cruzi, the protozoan agent of American trypanosomiasis (Chagas disease), has remarkable genetic and ecological diversity. A standardised MLST protocol that is suitable for assignment of T. cruzi isolates to genetic lineage and for higher resolution diversity studies has not been developed.

METHODOLOGY/PRINCIPAL FINDINGS

We have sequenced and diplotyped nine single copy housekeeping genes and assessed their value as part of a systematic MLST scheme for T. cruzi. A minimum panel of four MLST targets (Met-III, RB19, TcGPXII, and DHFR-TS) was shown to provide unambiguous assignment of isolates to the six known T. cruzi lineages (Discrete Typing Units, DTUs TcI-TcVI). In addition, we recommend six MLST targets (Met-II, Met-III, RB19, TcMPX, DHFR-TS, and TR) for more in depth diversity studies on the basis that diploid sequence typing (DST) with this expanded panel distinguished 38 out of 39 reference isolates. Phylogenetic analysis implies a subdivision between North and South American TcIV isolates. Single Nucleotide Polymorphism (SNP) data revealed high levels of heterozygosity among DTUs TcI, TcIII, TcIV and, for three targets, putative corresponding homozygous and heterozygous loci within DTUs TcI and TcIII. Furthermore, individual gene trees gave incongruent topologies at inter- and intra-DTU levels, inconsistent with a model of strict clonality.

CONCLUSIONS/SIGNIFICANCE

We demonstrate the value of systematic MLST diplotyping for describing inter-DTU relationships and for higher resolution diversity studies of T. cruzi, including presence of recombination events. The high levels of heterozygosity will facilitate future population genetics analysis based on MLST haplotypes.

Evolution of endogenous retroviruses in the Suidae: evidence for different viral subpopulations in African and Eurasian host species

Background

Porcine endogenous retroviruses (PERVs) represent remnants of an exogenous form that have become integrated in the domestic pig (Sus scrofa) genome. Although they are usually inactive, the capacity of γ1 ERVs to infect human cells in vitro has raised concerns about xenotransplantation because the viruses could cross the species barrier to humans. Here we have analyzed the evolution of γ1 ERVs in ten species of Suidae (suids, pigs and hogs) from Eurasia and Africa using DNA sequences for their coding domains (gag, pro/pol and env genes). For comparison with γ1 PERVs, we have also analysed γ2 ERVs which in domestic pigs are known to be inactive and do not pose a risk to xenotransplantation.

Results

Phylogenetic analysis using Bayesian inference showed that γ1 and γ2 ERVs have distinctive evolutionary histories. Firstly, two different viral lineages of γ1 ERVs were found and a coevolutionary analysis demonstrated that they correspond broadly to their host phylogeny, one of Eurasian and another of African species, and show no evidence of horizontal transmission. γ2 ERVs, however, show a bush-like evolution, suggesting a rapid viral radiation from a single common ancestor with no correspondence between host and viral evolutionary trees. Furthermore, though γ1 ERV env genes do not possess frequent stop codons, γ2 env genes do. To understand whether γ1 suid ERVs may be still replicating, we have also evaluated their likely mechanism of proliferation by statistically testing internal to terminal branches using nonsynonymous versus synonymous substitution ratios. Our results suggest that γ1 ERVs are increasing in copy number by reinfection, which requires the translocation of the virus from one cell to another.

Conclusions

Evidence of at least two viral subpopulations was observed in γ1 ERVs from Eurasian and African host species. These results should be taken into account in xenotransplantation since γ1 ERVs appear to be codiverging with their host and maintaining ongoing capacity to infect somatic and germ cells.

Analysing recombination in nucleotide sequences

Throughout the living world, genetic recombination and nucleotide substitution are the primary processes that create the genetic variation upon which natural selection acts. Just as analyses of substitution patterns can reveal a great deal about evolution, so too can analyses of recombination. Evidence of genetic recombination within the genomes of apparently asexual species can equate with evidence of cryptic sexuality. In sexually reproducing species, nonrandom patterns of sequence exchange can provide direct evidence of population subdivisions that prevent certain individuals from mating. Although an interesting topic in its own right, an important reason for analysing recombination is to account for its potentially disruptive influences on various phylogenetic-based molecular evolution analyses. Specifically, the evolutionary histories of recombinant sequences cannot be accurately described by standard bifurcating phylogenetic trees. Taking recombination into account can therefore be pivotal to the success of selection, molecular clock and various other analyses that require adequate modelling of shared ancestry and draw increased power from accurately inferred phylogenetic trees. Here, we review various computational approaches to studying recombination and provide guidelines both on how to gain insights into this important evolutionary process and on how it can be properly accounted for during molecular evolution studies.

Molecular data and ploidal levels indicate several putative allopolyploidization events in the genus Potentilla (Rosaceae)

Several naturally occurring hybrids in Potentilla (Rosaceae) have been reported, but no molecular evidence has so far been available to test these hypotheses of hybridization. We have compared a nuclear and a chloroplast gene tree to identify topological incongruences that may indicate hybridization events in the genus. Furthermore, the monophyly and phylogenetic position of the proposed segregated genera Argentina, Ivesia and Horkelia have been tested. The systematic signal from the two morphological characters, style- and anther shape, has also been investigated by ancestral state reconstruction, to elucidate how well these characters concur with the results of the molecular phylogenies. Six major clades, Anserina, Alba, Fragarioides, Reptans, ivesioid and Argentea, have been identified within genus Potentilla. Horkelia, Ivesia and Horkeliella (the ivesioid clade), form a monophyletic group nested within Potentilla. Furthermore, the origin of the proposed segregated genus Argentina (the Anserina clade) is uncertain but not in conflict with a new generic status of the group. We also found style morphology to be an informative character that reflects the phylogenetic relationships within Potentilla. Five well-supported incongruences were found between the nuclear and the chloroplast phylogenies, and three of these involved polyploid taxa. However, further investigations, using low copy molecular markers, are required to infer the phylogeny of these species and to test the hypothesis of hybrid origin.

Evidence of intra-segmental homologous recombination in influenza A virus

The evolution of influenza viruses is remarkably dynamic. Influenza viruses evolve rapidly in sequence and undergo frequent reassortment of different gene segments. Homologous recombination, although commonly seen as an important component of dynamic genome evolution in many other organisms, is believed to be rare in influenza. In this study, 256 gene segments from 32 influenza A genomes were examined for homologous recombination, three recombinant H1N1 strains were detected and they most likely resulted from one recombination event between two closely rated parental sequences. These findings suggest that homologous recombination in influenza viruses tends to take place between strains sharing high sequence similarity. The three recombinant strains were isolated at different time periods and they form a clade, indicating that recombinant strains could circulate. In addition, the simulation results showed that many recombinant sequences might not be detectable by currently existing recombinant detection programs when the parental sequences are of high sequence similarity. Finally, possible ways were discussed to improve the accuracy of the detection for recombinant sequences in influenza.

Molecular epidemiology of Plum pox virus in Japan

For a molecular epidemiological study based on complete genome sequences, 37 Plum pox virus (PPV) isolates were collected from the Kanto region in Japan. Pair-wise analyses revealed that all 37 Japanese isolates belong to the PPV-D strain, with low genetic diversity (less than 0.8%). In phylogenetic analysis of the PPV-D strain based on complete nucleotide sequences, the relationships of the PPV-D strain were reconstructed with high resolution: at the global level, the American, Canadian, and Japanese isolates formed their own distinct monophyletic clusters, suggesting that the routes of viral entry into these countries were independent; at the local level, the actual transmission histories of PPV were precisely reconstructed with high bootstrap support. This is the first description of the molecular epidemiology of PPV based on complete genome sequences.

Genetic diversity of symbiotic Bradyrhizobium elkanii populations recovered from inoculated and non-inoculated Acacia mangium field trials in Brazil

Acacia mangium is a legume tree native to Australasia. Since the eighties, it has been introduced into many tropical countries, especially in a context of industrial plantations. Many field trials have been set up to test the effects of controlled inoculation with selected symbiotic bacteria versus natural colonization with indigenous strains. In the introduction areas, A. mangium trees spontaneously nodulate with local and often ineffective bacteria. When inoculated, the persistence of inoculants and possible genetic recombination with local strains remain to be explored. The aim of this study was to describe the genetic diversity of bacteria spontaneously nodulating A. mangium in Brazil and to evaluate the persistence of selected strains used as inoculants. Three different sites, several hundred kilometers apart, were studied, with inoculated and non-inoculated plots in two of them. Seventy-nine strains were isolated from nodules and sequenced on three housekeeping genes (glnII, dnaK and recA) and one symbiotic gene (nodA). All but one of the strains belonged to the Bradyrhizobium elkanii species. A single case of housekeeping gene transfer was detected among the 79 strains, suggesting an extremely low rate of recombination within B. elkanii, whereas the nodulation gene nodA was found to be frequently transferred. The fate of the inoculant strains varied depending on the site, with a complete disappearance in one case, and persistence in another. We compared our results with the sister species Bradyrhizobium japonicum, both in terms of population genetics and inoculant strain destiny.

Evidence of recombination in the Banana bunchy top virus genome

Viruses serve as good model for evolutionary studies, owing to their short generation times and small genomes. Banana bunchy top virus (BBTV) is a significant subject being multicomponent circular single stranded DNA virus. BBTV belongs to family Nanoviridae and contains DNA-R, -U3, -S, -M, -C, and -N as integral genomic components. Evolutionary studies have shown genetic re-assortment of components among its isolates and revealed a concerted type evolution in non-coding regions of its genome. The DNA U3 having been shown as the most diverse component in our previous studies, was subjected to sequencing from some Pakistani isolates for the first time. Sequence analysis revealed intergenomic recombination in DNA-U3 among the isolates of two sub-groups and a very rare intragenomic recombination in Pakistani BBTV population. This indicates that like other evolutionary processes including intergenomic recombination, intragenomic recombination among the genomic components of the same isolate may also have a significant contribution in the evolution of BBTV genome. Intragenomic recombination therefore appears to be a unique way to generate genetic diversity in the multicomponent ssDNA viruses.

Recombination patterns in dicot-infecting mastreviruses mirror those found in monocot-infecting mastreviruses

Recombination has profoundly shaped the evolution of viruses in the family Geminiviridae and has been studied extensively in the two best characterised geminivirus lineages: the dicotyledonous plant infecting begomoviruses and the monocotyledonous plant infecting mastreviruses. Here, we demonstrate that the sizes and distributions of recombination events detectable within the members of a third major geminivirus lineage--the dicotyledonous plant infecting mastreviruses--are very similar to those of the monocot-infecting mastreviruses. This suggests that, despite host range differences, very similar biochemical, ecological and evolutionary factors must underlie recombination patterns in the dicot- and monocot-infecting mastreviruses.

An effective method for detecting gene conversion events in whole genomes

Gene conversion events are often overlooked in analyses of genome evolution. In a conversion event, an interval of DNA sequence (not necessarily containing a gene) overwrites a highly similar sequence. The event creates relationships among genomic intervals that can confound attempts to identify orthologs and to transfer functional annotation between genomes. Here we examine 1,616,329 paralogous pairs of mouse genomic intervals, and detect conversion events in about 7.5% of them. Properties of the putative gene conversions are analyzed, such as the lengths of the paralogous pairs and the spacing between their sources and targets. Our approach is illustrated using conversion events in primate CCL gene clusters. Source code for our program is included in the 3SEQ_2D package, which is freely available at www.bx.psu.edu/miller_lab/ .

Evaluation of methods for detecting conversion events in gene clusters

BACKGROUND

Gene clusters are genetically important, but their analysis poses significant computational challenges. One of the major reasons for these difficulties is gene conversion among the duplicated regions of the cluster, which can obscure their true relationships. Many computational methods for detecting gene conversion events have been released, but their performance has not been assessed for wide deployment in evolutionary history studies due to a lack of accurate evaluation methods.

RESULTS

We designed a new method that simulates gene cluster evolution, including large-scale events of duplication, deletion, and conversion as well as small mutations. We used this simulation data to evaluate several different programs for detecting gene conversion events.

CONCLUSIONS

Our evaluation identifies strengths and weaknesses of several methods for detecting gene conversion, which can contribute to more accurate analysis of gene cluster evolution.

Transmission of single and multiple viral variants in primary HIV-1 subtype C infection

To address whether sequences of viral gag and env quasispecies collected during the early post-acute period can be utilized to determine multiplicity of transmitted HIV's, recently developed approaches for analysis of viral evolution in acute HIV-1 infection [1], [2] were applied. Specifically, phylogenetic reconstruction, inter- and intra-patient distribution of maximum and mean genetic distances, analysis of Poisson fitness, shape of highlighter plots, recombination analysis, and estimation of time to the most recent common ancestor (tMRCA) were utilized for resolving multiplicity of HIV-1 transmission in a set of viral quasispecies collected within 50 days post-seroconversion (p/s) in 25 HIV-infected individuals with estimated time of seroconversion. The decision on multiplicity of HIV infection was made based on the model's fit with, or failure to explain, the observed extent of viral sequence heterogeneity. The initial analysis was based on phylogeny, inter-patient distribution of maximum and mean distances, and Poisson fitness, and was able to resolve multiplicity of HIV transmission in 20 of 25 (80%) cases. Additional analysis involved distribution of individual viral distances, highlighter plots, recombination analysis, and estimation of tMRCA, and resolved 4 of the 5 remaining cases. Overall, transmission of a single viral variant was identified in 16 of 25 (64%) cases, and transmission of multiple variants was evident in 8 of 25 (32%) cases. In one case multiplicity of HIV-1 transmission could not be determined. In primary HIV-1 subtype C infection, samples collected within 50 days p/s and analyzed by a single-genome amplification/sequencing technique can provide reliable identification of transmission multiplicity in 24 of 25 (96%) cases. Observed transmission frequency of a single viral variant and multiple viral variants were within the ranges of 64% to 68%, and 32% to 36%, respectively.

Computational analysis suggests that lyssavirus glycoprotein gene plays a minor role in viral adaptation

The Lyssavirus glycoprotein (G) is a membrane protein responsible for virus entry and protective immune responses. To explore possible roles of the glycoprotein in host shift or adaptation of Lyssavirus, we retrieved 53 full-length glycoprotein gene sequences from NCBI GenBank. The sequences were from different host isolates over a period of 70 years in 21 countries. Computational analyses detected 1 recombinant (AY987478, a dog isolate of CHAND03, genotype 1 in India) with incongruent phylogenetic support. No recombination was detected when AY98748 was excluded in the analyses. We applied different selection models to identify selection pressure on the glycoprotein gene. One codon at amino acid residual 483 was found to be under weak positive selection with marginal probability of 95% by using the maximum likelihood method. We found no significant evidence of positive selection on any site of the glycoprotein gene when the putative recombinant AY987478 was excluded. The computational analyses suggest that the G gene has been under purifying selection and that the evolution of the G gene may not play a significant role in Lyssavirus adaptation.

Influence of molecular resolution on sequence-based discovery of ecological diversity among Synechococcus populations in an alkaline siliceous hot spring microbial mat

Previous research has shown that sequences of 16S rRNA genes and 16S-23S rRNA internal transcribed spacer regions may not have enough genetic resolution to define all ecologically distinct Synechococcus populations (ecotypes) inhabiting alkaline, siliceous hot spring microbial mats. To achieve higher molecular resolution, we studied sequence variation in three protein-encoding loci sampled by PCR from 60°C and 65°C sites in the Mushroom Spring mat (Yellowstone National Park, WY). Sequences were analyzed using the ecotype simulation (ES) and AdaptML algorithms to identify putative ecotypes. Between 4 and 14 times more putative ecotypes were predicted from variation in protein-encoding locus sequences than from variation in 16S rRNA and 16S-23S rRNA internal transcribed spacer sequences. The number of putative ecotypes predicted depended on the number of sequences sampled and the molecular resolution of the locus. Chao estimates of diversity indicated that few rare ecotypes were missed. Many ecotypes hypothesized by sequence analyses were different in their habitat specificities, suggesting different adaptations to temperature or other parameters that vary along the flow channel.

Characterisation of potato virus Y isolates from Iran

A survey of Potato virus Y (PVY) was conducted in cultivated fields in six Iranian provinces between January 2005 to July 2007. Two hundred samples from potato and tomato were collected and analyzed by enzyme-linked immunosorbent assay (ELISA) for potyviruses. Almost one fourth of the samples were found to be infected by PVY. Analysis of these PVY-positive samples using three monoclonal antibodies (MAbs) facilitating the simultaneous detection of three main strains namely the ordinary (PVY(O)), strain (PVY(N)) and C (PVY(C)) strains. However, the fourth strain (PVY(NTN)) and some others recombinant isolates were also identified by molecular methods. Host range and symptoms analysis using sap inoculation of four different strains of PVY onto a range of plants revealed that the four strains showed biological properties that seemed to be consistent with their molecular grouping. Fourteen isolates of PVY were chosen based on the host and geographical location, primer specificity and serology for further biological and molecular characterisation. The coat protein (CP) and P1 genes and 3'-non-translated region (3'NTR) from 14 representative isolates were sequenced and analysed with the sequences available in GenBank. Composite analysis of the P1, CP and 3'-UTR sequences with all full genome sequences of PVY revealed that there are three potential strains of PVY in Iran, PVY(O), PVY(N)-W and PVY(NTN). Isolate KER.SA(N) was the most divergent of all the 14 isolates reacted with PVY(N) specific MAbs but grouped with PVY(O) strains in maximum likelihood phylogentic analysis. The PVY(NTN) isolates from Iran more closely related to the European than North American PVY(NTN) isolates.

First complete genomic characterization of two tick-borne encephalitis virus isolates obtained from wild rodents in South Korea

We determined for the first time the complete genome sequences of two Korean strains of the tick-borne encephalitis virus (TBEV), designated KrM 93 and KrM 213, isolated from the lung tissues of wild rodents in 2006. The genomes are 11,097 nucleotides (nt) in length and consist of a 132 nt 5'-noncoding region (NCR), a 10,245 nt open reading frame (ORF) containing 10 viral protein-coding regions (3,415 amino acids), and a 720 nt 3'-NCR. Compared with the 31 fully sequenced TBEV strains currently available, KrM 93 and KrM 213 show genomic nucleotide (and deduced amino acid) sequence divergences ranging from 1.8 (0.7) to 19.2 (26.6)% and 1.9 (0.8) to 19.3 (26.7)%, respectively. Phylogenetic and recombination analyses based on the complete genome sequence were performed to identify genetic variations and relationships between the TBEV strains. These showed that the Korean TBEV strains clustered with the Western subtype rather than with Far-Eastern or Siberian subtypes, and phylogenetic trees derived from capsid (C), envelope (E), nonstructural (NS) 4B and NS5 regions represented the same branching pattern shown by the complete genome-based tree. Although no recombination events were identified in these two Korean strains, 11 putative recombination events were identified within the NS5 regions or in the 3'-NCRs of TBEV strains in general. The results provide insight into the genetics of TBEV strains to understand the molecular epidemiology, genetic diversity, and evolution of TBEV.

Genomic analysis of hepatitis B virus reveals antigen state and genotype as sources of evolutionary rate variation

Hepatitis B virus (HBV) genomes are small, semi-double-stranded DNA circular genomes that contain alternating overlapping reading frames and replicate through an RNA intermediary phase. This complex biology has presented a challenge to estimating an evolutionary rate for HBV, leading to difficulties resolving the evolutionary and epidemiological history of the virus. Here, we re-examine rates of HBV evolution using a novel data set of 112 within-host, transmission history (pedigree) and among-host genomes isolated over 20 years from the indigenous peoples of the South Pacific, combined with 313 previously published HBV genomes. We employ Bayesian phylogenetic approaches to examine several potential causes and consequences of evolutionary rate variation in HBV. Our results reveal rate variation both between genotypes and across the genome, as well as strikingly slower rates when genomes are sampled in the Hepatitis B e antigen positive state, compared to the e antigen negative state. This Hepatitis B e antigen rate variation was found to be largely attributable to changes during the course of infection in the preCore and Core genes and their regulatory elements.

Distinctive differences in long terminal repeat sequences between γ1 endogenous retroviruses of African and Eurasian suid species

Diversity of long terminal repeats (LTRs) from γ1 endogenous retroviruses (ERVs) was analysed by DNA sequencing in 10 species of the family Suidae (suids, pigs and hogs). Phylogenetic analysis separated LTR sequences into two groups which correlated approximately with either the previously described cluster I and III, or the clusters II, IV and V. Interestingly, a specific LTR exhibiting a novel molecular rearrangement was identified exclusively within African host species when compared to LTRs previously reported from known ERVs in the domestic pig (Sus scrofa). Furthermore, other sections of LTRs appear to be unique to African suids as suggested by phylogenetic analysis. These differences between African and Eurasian ERV lineages show that these ERVs belong to different viral sub-populations, implying coevolution of endogenous viral sequences with their host species and providing no evidence of transfer of viral sequences between African and Eurasian suids.

Global genetic diversity and geographical and host-species distribution of beak and feather disease virus isolates

Psittacine beak and feather disease (PBFD) has a broad host range and is widespread in wild and captive psittacine populations in Asia, Africa, the Americas, Europe and Australasia. Beak and feather disease circovirus (BFDV) is the causative agent. BFDV has an ∼2 kb single stranded circular DNA genome encoding just two proteins (Rep and CP). In this study we provide support for demarcation of BFDV strains by phylogenetic analysis of 65 complete genomes from databases and 22 new BFDV sequences isolated from infected psittacines in South Africa. We propose 94% genome-wide sequence identity as a strain demarcation threshold, with isolates sharing >94% identity belonging to the same strain, and strain subtypes sharing >98% identity. Currently, BFDV diversity falls within 14 strains, with five highly divergent isolates from budgerigars probably representing a new species of circovirus with three strains (budgerigar circovirus; BCV-A, -B and -C). The geographical distribution of BFDV and BCV strains is strongly linked to the international trade in exotic birds; strains with more than one host are generally located in the same geographical area. Lastly, we examined BFDV and BCV sequences for evidence of recombination, and determined that recombination had occurred in most BFDV and BCV strains. We established that there were two globally significant recombination hotspots in the viral genome: the first is along the entire intergenic region and the second is in the C-terminal portion of the CP ORF. The implications of our results for the taxonomy and classification of circoviruses are discussed.

The spread of tomato yellow leaf curl virus from the Middle East to the world

The ongoing global spread of Tomato yellow leaf curl virus (TYLCV; Genus Begomovirus, Family Geminiviridae) represents a serious looming threat to tomato production in all temperate parts of the world. Whereas determining where and when TYLCV movements have occurred could help curtail its spread and prevent future movements of related viruses, determining the consequences of past TYLCV movements could reveal the ecological and economic risks associated with similar viral invasions. Towards this end we applied Bayesian phylogeographic inference and recombination analyses to available TYLCV sequences (including those of 15 new Iranian full TYLCV genomes) and reconstructed a plausible history of TYLCV's diversification and movements throughout the world. In agreement with historical accounts, our results suggest that the first TYLCVs most probably arose somewhere in the Middle East between the 1930s and 1950s (with 95% highest probability density intervals 1905-1972) and that the global spread of TYLCV only began in the 1980s after the evolution of the TYLCV-Mld and -IL strains. Despite the global distribution of TYLCV we found no convincing evidence anywhere other than the Middle East and the Western Mediterranean of epidemiologically relevant TYLCV variants arising through recombination. Although the region around Iran is both the center of present day TYLCV diversity and the site of the most intensive ongoing TYLCV evolution, the evidence indicates that the region is epidemiologically isolated, which suggests that novel TYLCV variants found there are probably not direct global threats. We instead identify the Mediterranean basin as the main launch-pad of global TYLCV movements.

Characterization, genetic diversity, and evolutionary link of Cucumber mosaic virus strain New Delhi from India

The genome of Cucumber mosaic virus New Delhi strain (CMV-ND) from India, obtained from tomato, was completely sequenced and compared with full genome sequences of 14 known CMV strains from subgroups I and II, for their genetic diversity. Sequence analysis suggests CMV-ND shares maximum sequence identity at the nucleotide level with a CMV strain from Taiwan. Among all 15 strains of CMV, the encoded protein 2b is least conserved, whereas the coat protein (CP) is most conserved. Sequence identity values and phylogram results indicate that CMV-ND belongs to subgroup I. Based on the recombination detection program result, it appears that CMV is prone to recombination, and different RNA components of CMV-ND have evolved differently. Recombinational analysis of all 15 CMV strains detected maximum recombination breakpoints in RNA2; CP showed the least recombination sites.

Statistical Analysis on Detecting Recombination Sites in DNA-β Satellites Associated with Old World Geminiviruses

Although exchange of genetic information by recombination plays an important role in the evolution of viruses, it is not clear how it generates diversity. Understanding recombination events helps with the study of the evolution of new virus strains or new viruses. Geminiviruses are plant viruses which have ambisense single-stranded circular DNA genomes and are one of the most economically important plant viruses in agricultural production. Small circular single-stranded DNA satellites, termed DNA-β, have recently been found to be associated with some geminivirus infections. In this paper we analyze several DNA-β sequences of geminiviruses for recombination events using phylogenetic and statistical analysis and we find that one strain from ToLCMaB has a recombination pattern and is a recombinant molecule between two strains from two species, PaLCuB-[IN:Chi:05] (major parent) and ToLCB-[IN:CP:04] (minor parent). We propose that this recombination event contributed to the evolution of the strain of ToLCMaB in South India. The Hidden Markov Chain (HMM) method developed by Webb et al. (2009) estimating phylogenetic tree through out the whole alignment provide us a recombination history of these DNA-β strains. It is the first time that this statistic method has been used on DNA-β recombination study and give a clear recombination history of DNA-β recombination.

The complete genome sequence and proteomics of Yersinia pestis phage Yep-phi

Yep-phi, a lytic phage of Yersinia pestis, was isolated in China and is routinely used as a diagnostic phage for the identification of the plague pathogen. Yep-phi has an isometric hexagonal head containing dsDNA and a short non-contractile conical tail. In this study, we sequenced the Yep-phi genome (GenBank accession no. HQ333270) and performed proteomics analysis. The genome consists of 38 ,616 bp of DNA, including direct terminal repeats of 222 bp, and is predicted to contain 45 ORFs. Most structural proteins were identified by proteomics analysis. Compared with the three available genome sequences of lytic phages for Y. pestis, the phages could be divided into two subgroups. Yep-phi displays marked homology to the bacteriophages Berlin (GenBank accession no. AM183667) and Yepe2 (GenBank accession no. EU734170), and these comprise one subgroup. The other subgroup is represented by bacteriophage ΦA1122 (GenBank accession no. AY247822). Potential recombination was detected among the Yep-phi subgroup.

Mitochondrial genome evolution in fire ants (Hymenoptera: Formicidae)

Background

Complete mitochondrial genome sequences have become important tools for the study of genome architecture, phylogeny, and molecular evolution. Despite the rapid increase in available mitogenomes, the taxonomic sampling often poorly reflects phylogenetic diversity and is often also biased to represent deeper (family-level) evolutionary relationships.

Results

We present the first fully sequenced ant (Hymenoptera: Formicidae) mitochondrial genomes. We sampled four mitogenomes from three species of fire ants, genus Solenopsis, which represent various evolutionary depths. Overall, ant mitogenomes appear to be typical of hymenopteran mitogenomes, displaying a general A+T-bias. The Solenopsis mitogenomes are slightly more compact than other hymentoperan mitogenomes (~15.5 kb), retaining all protein coding genes, ribosomal, and transfer RNAs. We also present evidence of recombination between the mitogenomes of the two conspecific Solenopsis mitogenomes. Finally, we discuss potential ways to improve the estimation of phylogenies using complete mitochondrial genome sequences.

Conclusions

The ant mitogenome presents an important addition to the continued efforts in studying hymenopteran mitogenome architecture, evolution, and phylogenetics. We provide further evidence that the sampling across many taxonomic levels (including conspecifics and congeners) is useful and important to gain detailed insights into mitogenome evolution. We also discuss ways that may help improve the use of mitogenomes in phylogenetic analyses by accounting for non-stationary and non-homogeneous evolution among branches.

RDP3: a flexible and fast computer program for analyzing recombination

RDP3 is a new version of the RDP program for characterizing recombination events in DNA-sequence alignments. Among other novelties, this version includes four new recombination analysis methods (3SEQ, VISRD, PHYLRO and LDHAT), new tests for recombination hot-spots, a range of matrix methods for visualizing over-all patterns of recombination within datasets and recombination-aware ancestral sequence reconstruction. Complementary to a high degree of analysis flow automation, RDP3 also has a highly interactive and detailed graphical user interface that enables more focused hands-on cross-checking of results with a wide variety of newly implemented phylogenetic tree construction and matrix-based recombination signal visualization methods. The new RDP3 can accommodate large datasets and is capable of analyzing alignments ranging in size from 1000 × 10 kilobase sequences to 20 × 2 megabase sequences within 48 h on a desktop PC.

AVAILABILITY

RDP3 is available for free from its web site http://darwin.uvigo.es/rdp/rdp.html.

Generation of genic diversity among Streptococcus pneumoniae strains via horizontal gene transfer during a chronic polyclonal pediatric infection

Although there is tremendous interest in understanding the evolutionary roles of horizontal gene transfer (HGT) processes that occur during chronic polyclonal infections, to date there have been few studies that directly address this topic. We have characterized multiple HGT events that most likely occurred during polyclonal infection among nasopharyngeal strains of Streptococcus pneumoniae recovered from a child suffering from chronic upper respiratory and middle-ear infections. Whole genome sequencing and comparative genomics were performed on six isolates collected during symptomatic episodes over a period of seven months. From these comparisons we determined that five of the isolates were genetically highly similar and likely represented a dominant lineage. We analyzed all genic and allelic differences among all six isolates and found that all differences tended to occur within contiguous genomic blocks, suggestive of strain evolution by homologous recombination. From these analyses we identified three strains (two of which were recovered on two different occasions) that appear to have been derived sequentially, one from the next, each by multiple recombination events. We also identified a fourth strain that contains many of the genomic segments that differentiate the three highly related strains from one another, and have hypothesized that this fourth strain may have served as a donor multiple times in the evolution of the dominant strain line. The variations among the parent, daughter, and grand-daughter recombinant strains collectively cover greater than seven percent of the genome and are grouped into 23 chromosomal clusters. While capturing in vivo HGT, these data support the distributed genome hypothesis and suggest that a single competence event in pneumococci can result in the replacement of DNA at multiple non-adjacent loci.

Bacterial infections have long been studied using Koch's postulates wherein the paradigm is that a single clone leads to a given infection. Over the past decade, it has become clear that chronic bacterial infections often do not fit this paradigm. Instead these are associated with the presence of multiple strains or species (polyclonal) of bacteria that are organized into highly structured communities, termed biofilms, which can persist in the body and are recalcitrant to antibiotic treatment. In addition, there is extensive evidence that bacteria can incorporate genes from neighboring bacteria into their own genomes. This process can produce new strains and is known as horizontal gene transfer. In this study, we investigated for the first time, the tempo and relevance of gene transfer among bacterial strains of Streptococcus pneumoniae during a naturally occurring chronic childhood infection. We identified extensive gene transfer among multiple infecting strains, by sequencing of isolates recovered sequentially over a seven-month period. This gene transfer may serve as a counterpoint to the host's adaptive immune response and help explain the phenomenon of bacterial persistence, since, as occurs with some chronic viral and parasitic infections, the immune system may become overwhelmed by a set of related strains.

No observed effect of homologous recombination on influenza C virus evolution

The occurrence of homologous recombination in influenza viruses has been under some debate recently. To determine the extent of homologous recombination in influenza C virus, recombination analyses of all available gene sequences of influenza C virus were carried out. No recombination signal was found. With the previous evidence in influenza A and B viruses, it seems that homologous recombination has minimal or no effect on influenza virus evolution.

Genetic and Pathological Diversity Among Xanthomonas Strains Responsible for Bacterial Spot on Tomato and Pepper in the Southwest Indian Ocean Region

Bacterial spot of tomato and pepper, a major problem in tropical climates, can be caused by several Xanthomonas genospecies. We examined the genetic and pathological diversity of a collection of 72 strains from the southwest Indian Ocean region as part of a regional research and development program to update inventories of agricultural pests and pathogens. Xanthomonas euvesicatoria, X. perforans, X. gardneri, and X. vesicatoria were identified in our strain collection. The identification of strains at the species level was consistently achieved by amplified fragment length polymorphism (AFLP) and multilocus sequence analysis (MLSA). Overall, X. euvesicatoria was the species recovered prevalently. MLSA data based on four housekeeping genes identified two to three sequence types per genospecies. It suggested that sequence variations primarily consisted of synonymous mutations, although a recombination event spanning several hundred nucleotides was detected for some strains of X. euvesicatoria on the atpD gene coding for the F1-F0-ATPase β subunit. The pathogenicity of strains was consistent with data found in the literature. Some pathological variations were primarily observed among strains identified as X. euvesicatoria. This study provides the first ever comprehensive description of the status of Xanthomonas species that cause bacterial spot of tomato and pepper in the southwest Indian Ocean region.

A 12-year molecular survey of clinical herpes simplex virus type 2 isolates demonstrates the circulation of clade A and B strains in Germany

BACKGROUND

Recently two different herpes simplex virus type 2 (HSV-2) clades (A and B) were described on DNA sequence data of the glycoprotein E (gE), G (gG) and I (gI) genes.

OBJECTIVE

To type the circulating HSV-2 wild-type strains in Germany by a novel approach and to monitor potential changes in the molecular epidemiology between 1997 and 2008.

STUDY DESIGN

A total of 64 clinical HSV-2 isolates were analyzed by a novel approach using the DNA sequences of the complete open reading frames of glycoprotein B (gB) and gG. Recombination analysis of the gB and gG gene sequences was performed to reveal intragenic recombinants.

RESULTS

Based on the phylogenetic analysis of the gB coding DNA sequence 8 of 64 (12%) isolates were classified as clade A strains and 56 of 64 (88%) isolates were classified as clade B strains. Analysis of the gG coding DNA sequence classified 4 (6%) isolates as clade A strains and 60 (94%) isolates as clade B strains. In comparison, the 8 isolates classified as clade A strains using the gB sequence data were classified as clade B strains when using the gG coding DNA sequence, suggesting intergenic recombination events. Intragenic recombination events were not detected.

CONCLUSION

The first molecular survey of clinical HSV-2 isolates from Germany demonstrated the circulation of clade A and B strains and of intergenic recombinants over a period of 12 years.

Guidelines for identifying homologous recombination events in influenza A virus

The rapid evolution of influenza viruses occurs both clonally and non-clonally through a variety of genetic mechanisms and selection pressures. The non-clonal evolution of influenza viruses comprises relatively frequent reassortment among gene segments and a more rarely reported process of non-homologous RNA recombination. Homologous RNA recombination within segments has been proposed as a third such mechanism, but to date the evidence for the existence of this process among influenza viruses has been both weak and controversial. As homologous recombination has not yet been demonstrated in the laboratory, supporting evidence, if it exists, may come primarily from patterns of phylogenetic incongruence observed in gene sequence data. Here, we review the necessary criteria related to laboratory procedures and sample handling, bioinformatic analysis, and the known ecology and evolution of influenza viruses that need to be met in order to confirm that a homologous recombination event occurred in the history of a set of sequences. To determine if these criteria have an effect on recombination analysis, we gathered 8307 publicly available full-length sequences of influenza A segments and divided them into those that were sequenced via the National Institutes of Health Influenza Genome Sequencing Project (IGSP) and those that were not. As sample handling and sequencing are executed to a very high standard in the IGSP, these sequences should be less likely to be exposed to contamination by other samples or by laboratory strains, and thus should not exhibit laboratory-generated signals of homologous recombination. Our analysis shows that the IGSP data set contains only two phylogenetically-supported single recombinant sequences and no recombinant clades. In marked contrast, the non-IGSP data show a very large amount of potential recombination. We conclude that the presence of false positive signals in the non-IGSP data is more likely than false negatives in the IGSP data, and that given the evidence to date, homologous recombination seems to play little or no role in the evolution of influenza A viruses.

Protein phosphatase 2B (PP2B, calcineurin) in Paramecium: partial characterization reveals that two members of the unusually large catalytic subunit family have distinct roles in calcium-dependent processes

We characterized the calcineurin (CaN) gene family, including the subunits CaNA and CaNB, based upon sequence information obtained from the Paramecium genome project. Paramecium tetraurelia has seven subfamilies of the catalytic CaNA subunit and one subfamily of the regulatory CaNB subunit, with each subfamily having two members of considerable identity on the amino acid level (>or=55% between subfamilies, >or=94% within CaNA subfamilies, and full identity in the CaNB subfamily). Within CaNA subfamily members, the catalytic domain and the CaNB binding region are highly conserved and molecular modeling revealed a three-dimensional structure almost identical to a human ortholog. At 14 members, the size of the CaNA family is unprecedented, and we hypothesized that the different CaNA subfamily members were not strictly redundant and that at least some fulfill different roles in the cell. This was tested by selecting two phylogenetically distinct members of this large family for posttranscriptional silencing by RNA interference. The two targets resulted in differing effects in exocytosis, calcium dynamics, and backward swimming behavior that supported our hypothesis that the large, highly conserved CaNA family members are not strictly redundant and that at least two members have evolved diverse but overlapping functions. In sum, the occurrence of CaN in Paramecium spp., although disputed in the past, has been established on a molecular level. Its role in exocytosis and ciliary beat regulation in a protozoan, as well as in more complex organisms, suggests that these roles for CaN were acquired early in the evolution of this protein family.

Ancestral polymorphism in exon 2 of bluethroat (Luscinia svecica) MHC class II B genes

The genes of the major histocompatibility complex (MHC) are important model genes for understanding selective forces in evolution. Here, we document, using a cloning and sequencing approach, high polymorphism at the exon 2 of the MHC class II B (MHCIIB) genes in the bluethroat (Luscinia svecica); a minimum of 61 unique alleles were detected in 20 individuals, and at least 11 functional loci. In addition, several pseudogenes were revealed. The specimens originated from three different bluethroat subspecies (azuricollis, cyanecula and svecica), and we also analysed four specimens of the closely related thrush nightingale (L. luscinia) for comparison. Phylogenetic analyses of the functional alleles revealed 258 equally parsimonious trees with poor statistical support for the majority of nodes. The distribution of the sequences in the trees point to an ancestral origin of the polymorphism in MHC class II B genes, a portion of which predated the phylogenetic split between the bluethroat and the thrush nightingale. Strong signatures of balancing selection were uncovered for the codons coding for the peptide-binding residues of the functional MHCIIB exon 2 alleles. Our results highlight the importance of duplication and recombination events for shaping passerine MHC and give insights in the evolutionary dynamics of MHC variation among closely related taxa.

The effect of vaccination on the evolution and population dynamics of avian paramyxovirus-1

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 19^th 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.

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.

Genome-scale phylogenetic analyses of chikungunya virus reveal independent emergences of recent epidemics and various evolutionary rates

Chikungunya virus (CHIKV), a mosquito-borne alphavirus, has traditionally circulated in Africa and Asia, causing human febrile illness accompanied by severe, chronic joint pain. In Africa, epidemic emergence of CHIKV involves the transition from an enzootic, sylvatic cycle involving arboreal mosquito vectors and nonhuman primates, into an urban cycle where peridomestic mosquitoes transmit among humans. In Asia, however, CHIKV appears to circulate only in the endemic, urban cycle. Recently, CHIKV emerged into the Indian Ocean and the Indian subcontinent to cause major epidemics. To examine patterns of CHIKV evolution and the origins of these outbreaks, as well as to examine whether evolutionary rates that vary between enzootic and epidemic transmission, we sequenced the genomes of 40 CHIKV strains and performed a phylogenetic analysis representing the most comprehensive study of its kind to date. We inferred that extant CHIKV strains evolved from an ancestor that existed within the last 500 years and that some geographic overlap exists between two main enzootic lineages previously thought to be geographically separated within Africa. We estimated that CHIKV was introduced from Africa into Asia 70 to 90 years ago. The recent Indian Ocean and Indian subcontinent epidemics appear to have emerged independently from the mainland of East Africa. This finding underscores the importance of surveillance to rapidly detect and control African outbreaks before exportation can occur. Significantly higher rates of nucleotide substitution appear to occur during urban than during enzootic transmission. These results suggest fundamental differences in transmission modes and/or dynamics in these two transmission cycles.