Evidence for positive selection in foot-and-mouth disease virus capsid genes from field isolates

Positive selection in the N-terminal extramembrane domain of lung surfactant protein C (SP-C) in marine mammals

Maximum-likelihood models of codon and amino acid substitution were used to analyze the lung-specific surfactant protein C (SP-C) from terrestrial, semi-aquatic, and diving mammals to identify lineages and amino acid sites under positive selection. Site models used the nonsynonymous/synonymous rate ratio (omega) as an indicator of selection pressure. Mechanistic models used physicochemical distances between amino acid substitutions to specify nonsynonymous substitution rates. Site models strongly identified positive selection at different sites in the polar N-terminal extramembrane domain of SP-C in the three diving lineages: site 2 in the cetaceans (whales and dolphins), sites 7, 9, and 10 in the pinnipeds (seals and sea lions), and sites 2, 9, and 10 in the sirenians (dugongs and manatees). The only semi-aquatic contrast to indicate positive selection at site 10 was that including the polar bear, which had the largest body mass of the semi-aquatic species. Analysis of the biophysical properties that were influential in determining the amino acid substitutions showed that isoelectric point, chemical composition of the side chain, polarity, and hydrophobicity were the crucial determinants. Amino acid substitutions at these sites may lead to stronger binding of the N-terminal domain to the surfactant phospholipid film and to increased adsorption of the protein to the air-liquid interface. Both properties are advantageous for the repeated collapse and reinflation of the lung upon diving and resurfacing and may reflect adaptations to the high hydrostatic pressures experienced during diving.

Genome-wide diversity and selective pressure in the human rhinovirus

Background

The human rhinoviruses (HRV) are one of the most common and diverse respiratory pathogens of humans. Over 100 distinct HRV serotypes are known, yet only 6 genomes are available. Due to the paucity of HRV genome sequence, little is known about the genetic diversity within HRV or the forces driving this diversity. Previous comparative genome sequence analyses indicate that recombination drives diversification in multiple genera of the picornavirus family, yet it remains unclear if this holds for HRV.

Results

To resolve this and gain insight into the forces driving diversification in HRV, we generated a representative set of 34 fully sequenced HRVs. Analysis of these genomes shows consistent phylogenies across the genome, conserved non-coding elements, and only limited recombination. However, spikes of genetic diversity at both the nucleotide and amino acid level are detectable within every locus of the genome. Despite this, the HRV genome as a whole is under purifying selective pressure, with islands of diversifying pressure in the VP1, VP2, and VP3 structural genes and two non-structural genes, the 3C protease and 3D polymerase. Mapping diversifying residues in these factors onto available 3-dimensional structures revealed the diversifying capsid residues partition to the external surface of the viral particle in statistically significant proximity to antigenic sites. Diversifying pressure in the pleconaril binding site is confined to a single residue known to confer drug resistance (VP1 191). In contrast, diversifying pressure in the non-structural genes is less clear, mapping both nearby and beyond characterized functional domains of these factors.

Conclusion

This work provides a foundation for understanding HRV genetic diversity and insight into the underlying biology driving evolution in HRV. It expands our knowledge of the genome sequence space that HRV reference serotypes occupy and how the pattern of genetic diversity across HRV genomes differs from other picornaviruses. It also reveals evidence of diversifying selective pressure in both structural genes known to interact with the host immune system and in domains of unassigned function in the non-structural 3C and 3D genes, raising the possibility that diversification of undiscovered functions in these essential factors may influence HRV fitness and evolution.

Genetic variability and molecular evolution of hepatitis A virus

Hepatitis A virus (HAV), the causative agent of type A viral hepatitis, was first identified about three decades ago. Recent findings have shown that HAV possess several characteristics that make it unique among the family Picornaviridae, particularly in terms of its mechanisms of polyprotein processing and virion morphogenesis. HAV circulates in vivo as distributions of closely genetically related variants referred to as quasispecies. HAV exploits all known mechanisms of genetic variation to ensure its survival, including mutation and recombination. Only one serotype and six different genetic groups (three humans and three simian) have been described. HAV mutation rate is significantly lower as compared to other members of the family Picornaviridae. The mode of evolution appears, at least in part, to contribute to the presence of only one known serotype.

Unravelling selection shifts among foot-and-mouth disease virus (FMDV) serotypes

FMDV virus has been increasingly recognised as the most economically severe animal virus with a remarkable degree of antigenic diversity. Using an integrative evolutionary and computational approach we have compelling evidence for heterogeneity in the selection forces shaping the evolution of the seven different FMDV serotypes. Our results show that positive Darwinian selection has governed the evolution of the major antigenic regions of serotypes A, Asia1, O, SAT1 and SAT2, but not C or SAT3. Co-evolution between sites from antigenic regions under positive selection pinpoints their functional communication to generate immune-escape mutants while maintaining their ability to recognise the host-cell receptors. Neural network and functional divergence analyses strongly point to selection shifts between the different serotypes. Our results suggest that, unlike African FMDV serotypes, serotypes with wide geographical distribution have accumulated compensatory mutations as a strategy to ameliorate the effect of slightly deleterious mutations fixed by genetic drift. This strategy may have provided the virus by a flexibility to generate immune-escape mutants and yet recognise host-cell receptors. African serotypes presented no evidence for compensatory mutations. Our results support heterogeneous selective constraints affecting the different serotypes. This points to the possible accelerated rates of evolution diverging serotypes sharing geographical locations as to ameliorate the competition for the host.

Recombination patterns in aphthoviruses mirror those found in other picornaviruses

Foot-and-mouth disease virus (FMDV) is thought to evolve largely through genetic drift driven by the inherently error-prone nature of its RNA polymerase. There is, however, increasing evidence that recombination is an important mechanism in the evolution of these and other related picornoviruses. Here, we use an extensive set of recombination detection methods to identify 86 unique potential recombination events among 125 publicly available FMDV complete genome sequences. The large number of events detected between members of different serotypes suggests that horizontal flow of sequences among the serotypes is relatively common and does not incur severe fitness costs. Interestingly, the distribution of recombination breakpoints was found to be largely nonrandom. Whereas there are clear breakpoint cold spots within the structural genes, two statistically significant hot spots precisely separate these from the nonstructural genes. Very similar breakpoint distributions were found for other picornovirus species in the genera Enterovirus and Teschovirus. Our results suggest that genome regions encoding the structural proteins of both FMDV and other picornaviruses are functionally interchangeable modules, supporting recent proposals that the structural and nonstructural coding regions of the picornaviruses are evolving largely independently of one another.

Recombination and selection in the evolution of picornaviruses and other Mammalian positive-stranded RNA viruses

Picornaviridae are a large virus family causing widespread, often pathogenic infections in humans and other mammals. Picornaviruses are genetically and antigenically highly diverse, with evidence for complex evolutionary histories in which recombination plays a major part. To investigate the nature of recombination and selection processes underlying the evolution of serotypes within different picornavirus genera, large-scale analysis of recombination frequencies and sites, segregation by serotype within each genus, and sequence selection and composition was performed, and results were compared with those for other nonenveloped positive-stranded viruses (astroviruses and human noroviruses) and with flavivirus and alphavirus control groups. Enteroviruses, aphthoviruses, and teschoviruses showed phylogenetic segregation by serotype only in the structural region; lack of segregation elsewhere was attributable to extensive interserotype recombination. Nonsegregating viruses also showed several characteristic sequence divergence and composition differences between genome regions that were absent from segregating virus control groups, such as much greater amino acid sequence divergence in the structural region, markedly elevated ratios of nonsynonymous-to-synonymous substitutions, and differences in codon usage. These properties were shared with other picornavirus genera, such as the parechoviruses and erboviruses. The nonenveloped astroviruses and noroviruses similarly showed high frequencies of recombination, evidence for positive selection, and differential codon use in the capsid region, implying similar underlying evolutionary mechanisms and pressures driving serotype differentiation. This process was distinct from more-recent sequence evolution generating diversity within picornavirus serotypes, in which neutral or purifying selection was prominent. Overall, this study identifies common themes in the diversification process generating picornavirus serotypes that contribute to understanding of their evolution and pathogenicity.

Population genetics of microbial pathogens estimated from multilocus sequence typing (MLST) data

The inference of population recombination (rho), population mutation (Theta), and adaptive selection is of great interest in microbial population genetics. These parameters can be efficiently estimated using explicit statistical frameworks (evolutionary models) that describe their effect on gene sequences. Within this framework, we estimated rho and Theta using a coalescent approach, and adaptive (or destabilizing) selection under heterogeneous codon-based and amino acid property models in microbial sequences from MLST databases. We analyzed a total of 91 different housekeeping gene regions (loci) corresponding to one fungal and sixteen bacterial pathogens. Our results show that these three population parameters vary extensively across species and loci, but they do not seem to be correlated. For the most part, estimated recombination rates among species agree well with previous studies. Over all taxa, the rho/Theta ratio suggests that each factor contributes similarly to the emergence of variant alleles. Comparisons of Theta estimated under finite- and infinite-site models indicate that recurrent mutation (i.e., multiple mutations at some sites) can increase Theta by up to 39%. Significant evidence of molecular adaptation was detected in 28 loci from 13 pathogens. Three of these loci showed concordant patterns of adaptive selection in two to four different species.

Molecular Signatures of Natural Selection

There is an increasing interest in detecting genes, or genomic regions, that have been targeted by natural selection. The interest stems from a basic desire to learn more about evolutionary processes in humans and other organisms, and from the realization that inferences regarding selection may provide important functional information. This review provides a nonmathematical description of the issues involved in detecting selection from DNA sequences and SNP data and is intended for readers who are not familiar with population genetic theory. Particular attention is placed on issues relating to the analysis of large-scale genomic data sets.

Adaptive Diversification of Vomeronasal Receptor 1 Genes in Rodents

The vomeronasal receptor 1 (V1R) are believed to be pheromone receptors in rodents. Here we used computational methods to identify 95 and 62 new putative V1R genes from the draft rat and mouse genome sequence, respectively. The rat V1R repertoire consists of 11 subfamilies, 10 of which are shared with the mouse, while rat appears to lack the H and I subfamilies found in mouse and possesses one unique subfamily (M). The estimations of the relative divergence times suggest that many subfamilies originated after the split of rodents and primates. The analysis also reveals that these clusters underwent an expansion very close to the split of mouse and rat. In addition, maximum likelihood analysis showed that the nonsynonymous and synonymous rate ratio for most of these clusters was much higher than one, suggesting the role of positive selection in the diversification of these duplicated V1R genes. Because V1R are thought to mediate the process of signal transduction in response to pheromone detection, we speculate that the V1R genes have evolved under positive Darwinian selection to maintain the ability to discriminate between large and complex pheromonal mixtures.

Molecular evolution of the hepatitis delta virus antigen gene: recombination or positive selection?

We present the statistical analysis of diversifying selective pressures on the hepatitis D antigen gene (HDAg). Thirty-three distinct HDAg sequences from subtypes I, II, and III were tested for positive selection using maximum likelihood methods based on models of codon substitution that allow variable selective pressures across sites. Such methods have been shown to be sufficiently accurate and successful in detecting positive selection in a variety of viral and nonviral protein-coding genes. About 11% of codon sites in HDAg were estimated to be under diversifying selection. Remarkably, most of the residues predicted to evolve under positive selection were located in the immunogenic domain and the N-terminus region with reported antigenic activity. These sites are potential targets of the host's immune response. Identification of residues mutating to escape immune recognition may help to distinguish the most virulent strains and aid vaccine design. Possible interplay between positive selection and recombination on the gene is discussed but no significant evidence for recombination was found.

Adaptive evolution of cytochrome c oxidase: Infrastructure for a carnivorous plant radiation

Much recent attention in the study of adaptation of organismal form has centered on developmental regulation. As such, the highly conserved respiratory machinery of eukaryotic cells might seem an unlikely target for selection supporting novel morphologies. We demonstrate that a dramatic molecular evolutionary rate increase in subunit I of cytochrome c oxidase (COX) from an active-trapping lineage of carnivorous plants is caused by positive Darwinian selection. Bladderworts (Utricularia) trap plankton when water-immersed, negatively pressured suction bladders are triggered. The resetting of traps involves active ion transport, requiring considerable energy expenditure. As judged from the quaternary structure of bovine COX, the most profound adaptive substitutions are two contiguous cysteines absent in approximately 99.9% of databased COX I sequences from Eukaryota, Archaea, and Bacteria. This motif lies directly at the docking point of COX I helix 3 and cytochrome c, and modeling of bovine COX I suggests the possibility of an unprecedented helix-terminating disulfide bridge that could alter COX/cytochrome c dissociation kinetics. Thus, the key adaptation in Utricularia likely lies in molecular energetic changes that buttressed the mechanisms responsible for the bladderworts' radical morphological evolution. Along with evidence for COX evolution underlying expansion of the anthropoid neocortex, our findings underscore that important morphological and physiological innovations must often be accompanied by specific adaptations in proteins with basic cellular functions.

Darwinian adaptation of proteorhodopsin to different light intensities in the marine environment

Proteorhodopsin, a retinal-binding protein, represents a potentially significant source of light-driven energy production in the world's oceans. The distribution of photochemically divergent proteorhodopsins is stratified according to depth. Here, we present evidence that such photochemical diversity was tuned by Darwinian selection. By using a Bayesian method, we identified sites targeted by Darwinian selection and mapped them to three-dimensional models of proteorhodopsins. We suggest that spectral fine-tuning results from the combined effect of amino acids that directly interact with retinal and those that influence the confirmation of the retinal-binding pocket.

Molecular epidemiology of foot-and-mouth disease viruses in the Adamawa province of Cameroon

Foot-and-mouth disease virus (FMDV) causes a highly contagious viral disease of even-toed ungulates and is one of the most important economic diseases of livestock. Most studies of FMDV are done in countries where control measures are being implemented. In contrast, in areas such as sub-Saharan Africa, where FMDV is endemic and new strains are likely to emerge, there are only sporadic submissions to the World Reference Laboratory, Pirbright, United Kingdom. This paper describes the molecular epidemiology of FMDV in the Adamawa province of Cameroon based on a population sample of cattle herds. Serotypes SAT2 and A were isolated in the cross-sectional study. SAT2 isolates were all similar, with phylogenetic distances of <6%, and were most closely related to published sequences of isolates from Eritrea and Saudi Arabia. Serotype A isolates were more variable, with phylogenetic distances of 0 to 11%, and were most closely related to historic isolates from Cameroon. Use of a population-based sample gives a representative sample of virus diversity and will improve our understanding of the evolution of FMDV and its epidemiology. A supplementary study of pigs passing through the railhead collection yard at Ngaoundere detected a serotype O virus. A third pilot longitudinal study monitored viral persistence in three cattle herds over 12 months, and serotype O and A viruses were recovered from a herd 12 months after it was first recorded as being infected with SAT2 virus. The pig type O isolate was not closely related to that recovered from the cattle, suggesting that the pigs had not introduced the O virus into the cattle herds.

Evidence of unique genotypes of beak and feather disease virus in southern Africa

Psittacine beak and feather disease (PBFD), caused by Beak and feather disease virus (BFDV), is the most significant infectious disease in psittacines. PBFD is thought to have originated in Australia but is now found worldwide; in Africa, it threatens the survival of the indigenous endangered Cape parrot and the vulnerable black-cheeked lovebird. We investigated the genetic diversity of putative BFDVs from southern Africa. Feathers and heparinized blood samples were collected from 27 birds representing 9 psittacine species, all showing clinical signs of PBFD. DNA extracted from these samples was used for PCR amplification of the putative BFDV coat protein (CP) gene. The nucleotide sequences of the CP genes of 19 unique BFDV isolates were determined and compared with the 24 previously described sequences of BFDV isolates from Australasia and America. Phylogenetic analysis revealed eight BFDV lineages, with the southern African isolates representing at least three distinctly unique genotypes; 10 complete genome sequences were determined, representing at least one of every distinct lineage. The nucleotide diversity of the southern African isolates was calculated to be 6.4% and is comparable to that found in Australia and New Zealand. BFDVs in southern Africa have, however, diverged substantially from viruses found in other parts of the world, as the average distance between the southern African isolates and BFDV isolates from Australia ranged from 8.3 to 10.8%. In addition to point mutations, recombination was found to contribute substantially to the level of genetic variation among BFDVs, with evidence of recombination in all but one of the genomes analyzed.

Modeling the site-specific variation of selection patterns along lineages

The unambiguous footprint of positive Darwinian selection in protein-coding DNA sequences is revealed by an excess of nonsynonymous substitutions over synonymous substitutions compared with the neutral expectation. Methods for analyzing the patterns of nonsynonymous and synonymous substitutions usually rely on stochastic models in which the selection regime may vary across the sequence but remains constant across lineages for any amino acid position. Despite some work that has relaxed the constraint that selection patterns remain constant over time, no model provides a strong statistical framework to deal with switches between selection processes at individual sites during the course of evolution. This paper describes an approach that allows the site-specific selection process to vary along lineages of a phylogenetic tree. The parameters of the switching model of codon substitution are estimated by using maximum likelihood. The analysis of eight HIV-1 env homologous sequence data sets shows that this model provides a significantly better fit to the data than one that does not take into account switches between selection patterns in the phylogeny at individual sites. We also provide strong evidence that the strength and the frequency of occurrence of selection might not be estimated accurately when the site-specific variation of selection regimes is ignored.

Foot-and-mouth disease

Foot-and-mouth disease (FMD) is a highly contagious disease of cloven-hoofed animals. The disease was initially described in the 16th century and was the first animal pathogen identified as a virus. Recent FMD outbreaks in developed countries and their significant economic impact have increased the concern of governments worldwide. This review describes the reemergence of FMD in developed countries that had been disease free for many years and the effect that this has had on disease control strategies. The etiologic agent, FMD virus (FMDV), a member of the Picornaviridae family, is examined in detail at the genetic, structural, and biochemical levels and in terms of its antigenic diversity. The virus replication cycle, including virus-receptor interactions as well as unique aspects of virus translation and shutoff of host macromolecular synthesis, is discussed. This information has been the basis for the development of improved protocols to rapidly identify disease outbreaks, to differentiate vaccinated from infected animals, and to begin to identify and test novel vaccine candidates. Furthermore, this knowledge, coupled with the ability to manipulate FMDV genomes at the molecular level, has provided the framework for examination of disease pathogenesis and the development of a more complete understanding of the virus and host factors involved.

Widespread adaptive evolution in the human immunodeficiency virus type 1 genome

We investigated variable selective pressures among amino acid sites in HIV-1 genes. Selective pressure at the amino acid level was measured by using the nonsynonymous/synonymous substitution rate ratio (omega = dN/dS). To identify amino acid sites under positive selection with omega > 1, we applied maximum likelihood models that allow variable omega ratios among sites to analyze genomic sequences of 26 HIV-1 lineages including subtypes A, B, and C. Likelihood ratio tests detected sites under positive selection in each of the major genes in the genome: env, gag, pol, vif, and vpr. Positive selection was also detected in nef, tat, and vpu, although those genes are very small. The majority of positive selection sites is located in gp160. Positive selection was not detected if omega was estimated as an average across all sites, indicating the lack of power of the averaging approach. Candidate positive selection sites were mapped onto the available protein tertiary structures and immunogenic epitopes. We measured the physiochemical properties of amino acids and found that those at positive selection sites were more diverse than those at variable sites. Furthermore, amino acid residues at exposed positive selection sites were more physiochemically diverse than at buried positive selection sites. Our results demonstrate genomewide diversifying selection acting on the HIV-1.

Unifying the epidemiological and evolutionary dynamics of pathogens

A key priority for infectious disease research is to clarify how pathogen genetic variation, modulated by host immunity, transmission bottlenecks, and epidemic dynamics, determines the wide variety of pathogen phylogenies observed at scales that range from individual host to population. We call the melding of immunodynamics, epidemiology, and evolutionary biology required to achieve this synthesis pathogen "phylodynamics." We introduce a phylodynamic framework for the dissection of dynamic forces that determine the diversity of epidemiological and phylogenetic patterns observed in RNA viruses of vertebrates. A central pillar of this model is the Evolutionary Infectivity Profile, which captures the relationship between immune selection and pathogen transmission.

Genetic variability of hepatitis A virus

Knowledge of the molecular biology of hepatitis A virus (HAV) has increased exponentially since its identification. HAV exploits all known mechanisms of genetic variation to ensure survival, including mutation and genetic recombination. HAV has been characterized by the emergence of different genotypes, three human antigenic variants and only one major serotype. This paper reviews the genetic variability and molecular epidemiology of HAV. Its evolutionary mechanisms are described with particular emphasis on genetic recombination and HAV mutation rate. Genotypic classification methods are also discussed.

Parasitic exploitation as an engine of diversity

Parasitic exploitation occurs within and between a wide variety of taxa in a plethora of diverse contexts. Theoretical and empirical analyses indicate that parasitic exploitation can generate substantial genetic and phenotypic polymorphism within species. Under some circumstances, parasitic exploitation may also be an important factor causing reproductive isolation and promoting speciation. Here we review research relevant to the relationship between parasitic exploitation, within species-polymorphism, and speciation in some of the major arenas in which such exploitation has been studied. This includes research on the vertebrate major histocompatibility loci, plant-pathogen interactions, the evolution of sexual reproduction, intragenomic conflict, sexual conflict, kin mimicry and social parasitism, tropical forest diversity and the evolution of language. We conclude by discussing some of the issues raised by comparing the effect of parasitic exploitation on polymorphism and speciation in different contexts.

Effect of recombination on the accuracy of the likelihood method for detecting positive selection at amino acid sites

Maximum-likelihood methods based on models of codon substitution accounting for heterogeneous selective pressures across sites have proved to be powerful in detecting positive selection in protein-coding DNA sequences. Those methods are phylogeny based and do not account for the effects of recombination. When recombination occurs, such as in population data, no unique tree topology can describe the evolutionary history of the whole sequence. This violation of assumptions raises serious concerns about the likelihood method for detecting positive selection. Here we use computer simulation to evaluate the reliability of the likelihood-ratio test (LRT) for positive selection in the presence of recombination. We examine three tests based on different models of variable selective pressures among sites. Sequences are simulated using a coalescent model with recombination and analyzed using codon-based likelihood models ignoring recombination. We find that the LRT is robust to low levels of recombination (with fewer than three recombination events in the history of a sample of 10 sequences). However, at higher levels of recombination, the type I error rate can be as high as 90%, especially when the null model in the LRT is unrealistic, and the test often mistakes recombination as evidence for positive selection. The test that compares the more realistic models M7 (beta) against M8 (beta and omega) is more robust to recombination, where the null model M7 allows the positive selection pressure to vary between 0 and 1 (and so does not account for positive selection), and the alternative model M8 allows an additional discrete class with omega = d(N)/d(S) that could be estimated to be >1 (and thus accounts for positive selection). Identification of sites under positive selection by the empirical Bayes method appears to be less affected than the LRT by recombination.

Evidence of recombination in natural populations of hepatitis A virus

Genetic analysis of selected genome regions of hepatitis A virus (HAV) suggested that distinct genotypes of HAV could be found in different geographical regions. At least seven HAV genotypes have been identified all over the world, including four human genotypes (I, II, III, and VII) and three simian strains (IV, V, and VI). Phylogenetic analysis using full-length VP1 sequences revealed that human strain 9F94 has a close genetic relation with strain SLF-88 (sub-genotype VII). Nevertheless, the same analysis using full-length VP2 or VP3 sequences revealed that strain 9F94 has a close genetic relation with strain MBB (sub-genotype IB). To test the possibility of genetic recombination, phylogenetic studies were carried out, revealing that a crossing over had taken place in the VP1 capsid protein. These findings indicate that capsid-recombination can play a significant role in shaping the genetic diversity of HAV and, as such, can have important implications for its evolution, biology, and control.

The implications of virus diversity within the SAT 2 serotype for control of foot-and-mouth disease in sub-Saharan Africa

SAT 2 is the serotype most often associated with outbreaks of foot-and-mouth disease (FMD) in livestock in southern and western Africa and is the only SAT type to have been recorded outside the African continent in the last decade. Its epidemiology is complicated by the presence of African buffalo (Syncerus caffer), which play an important role in virus maintenance and transmission. To assess the level of genetic complexity of this serotype among viruses associated with both domestic livestock and wildlife, complete VP1 gene sequences of 53 viruses from 17 countries and three different host species were analysed. Phylogenetic analysis revealed eleven virus lineages, differing from each other by at least 20 % in pairwise nucleotide comparisons, four of which fall within the southern African region, two in West Africa and the remaining five in central and East Africa. No evidence of recombination between these lineages was detected, and thus we conclude that these are independently evolving virus lineages which occur primarily in discrete geographical localities in accordance with the FMD virus topotype concept. Applied to the whole phylogeny, rates of nucleotide substitution are significantly different between topotypes, but most individual topotypes evolve in accordance with a molecular clock at an average rate of approximately 0.002 substitutions per site per year. This study provides an indication of the intratypic complexity of the SAT 2 serotype at the continental level and emphasizes the value of molecular characterization of diverse FMD field strains for tracing the origin of outbreaks.

Emergence of a new strain of type O foot-and-mouth disease virus: its phylogenetic and evolutionary relationship with the PanAsia pandemic strain

In India, Foot-and-mouth disease virus (FMDV) serotype O has been associated with more than 75% of the outbreaks. Previous studies with this serotype have indicated that the viruses circulating in India belong to a single genotype. Recent (February 2001) FMD epidemics in Europe have focussed global attention on the source of the virus and have been traced to a strain, PanAsia (serotype O), which is present in India since 1990. In this study, to further characterize the isolates belonging to the PanAsian strain, we sequenced the complete VP1-encoding (1 D) gene for 71 FMDV serotype O isolates from India recovered from the field outbreaks during the last 4 decades (1962-2001). All the isolates in the tree were distributed in to three major branches (designated as A, B and C); the branch C is further divided into four groups (I-IV), of which the group IV belongs to the PanAsia strain. Furthermore, we show that the PanAsia strain has been circulating endemically since 1982 (not 1990 as reported earlier) and has been the most dominant outbreak strain in the recent years and distributed at least in 17 states of the country. During the year 2001, another new group (group III) of virus with genetic divergence of 5.4-11.1% at nucleotide level from the PanAsia strain is found to co-circulate endemically, and is slowly replacing it. At amino acid level this strain differed from PanAsia strain at five amino acid positions in the VP1. Although these strains are divergent at nucleotide level, they maintained a good antigenic relationship with one of the vaccine strains (IND R2/75) widely used in the country. Given the ability of the PanAsia virus to persist, spread and to outcompete other strains, the present trend could be of serious concern as the newly emerging virus is replacing it. If this is true, then there is another equally divergent strain as PanAsia that may pose a serious threat to the global dairy and meat industries.

Molecular evolution meets the genomics revolution

Changes in technology in the past decade have had such an impact on the way that molecular evolution research is done that it is difficult now to imagine working in a world without genomics or the Internet. In 1992, GenBank was less than a hundredth of its current size and was updated every three months on a huge spool of tape. Homology searches took 30 minutes and rarely found a hit. Now it is difficult to find sequences with only a few homologs to use as examples for teaching bioinformatics. For molecular evolution researchers, the genomics revolution has showered us with raw data and the information revolution has given us the wherewithal to analyze it. In broad terms, the most significant outcome from these changes has been our newfound ability to examine the evolution of genomes as a whole, enabling us to infer genome-wide evolutionary patterns and to identify subsets of genes whose evolution has been in some way atypical.

Maximum likelihood methods for detecting adaptive evolution after gene duplication

The rapid accumulation of genomic sequences in public databases will finally allow large scale studies of gene family evolution, including evaluation of the role of positive Darwinian selection following a duplication event. This will be possible because recent statistical methods of comparing synonymous and nonsynonymous substitution rates permit reliable detection of positive selection at individual amino acid sites and along evolutionary lineages. Here, we summarize maximum-likelihood based methods, and present a framework for their application to analysis of gene families. Using these methods, we investigated the role of positive Darwinian selection in the ECP-EDN gene family of primates and the Troponin C gene family of vertebrates. We also comment on the limitations of these methods and discuss directions for further improvements.

The evolutionary genomics of pathogen recombination

A pressing problem in studying the evolution of microbial pathogens is to determine the extent to which these genomes recombine. This information is essential for locating pathogenicity loci by using association studies or population genetic approaches. Recombination also complicates the use of phylogenetic approaches to estimate evolutionary parameters such as selection pressures. Reliable methods that detect and estimate the rate of recombination are, therefore, vital. This article reviews the approaches that are available for detecting and estimating recombination in microbial pathogens and how they can be used to understand pathogen evolution and to identify medically relevant loci.

Genetic variation and immunohistochemical differences among geographic isolates of Taura syndrome virus of penaeid shrimp

Taura syndrome virus (TSV) is an important virus infecting penaeid shrimp in the western hemisphere. Genetic variation and immunohistochemical differences of 20 TSV isolates collected from the USA, Taiwan, Mexico and Nicaragua were compared. Capsid protein genes CP1 (546 bp) and CP2 (584 bp) were amplified by RT-PCR and the cDNAs were sequenced. Pairwise comparison of nucleotide sequences showed a 0-2.4% difference in CP1 and a 0-3.5% difference in CP2. Phylogenetic analyses clustered the TSV isolates into two groups: one contained USA, Taiwan and some Mexican isolates, the other contained Mexican isolates only. Immunohistochemical analysis using a TSV-specific monoclonal antibody produced positive results for the USA and Taiwan isolates but negative results for the Mexican and Nicaraguan isolates. Molecular and immunohistochemical data suggest the existence of at least two TSV strains, one of which might have evolved following contact with a new penaeid host, Penaeus stylirostris.

Diversity and evolution of the envelope gene of dengue virus type 1

The genetic diversity and phylogenetic relationships of a collection of strains of dengue virus type 1 (DV-1), isolated from different parts of the world, were investigated. Phylogenetic trees derived from the complete sequence of the E gene of 44 strains suggested the existence of five genetic types defined by a maximum nucleotide divergence within each group of 6%. The 22 strains from America were classified into a single genetic type that included strains associated either with classical dengue or hemorrhagic dengue episodes. Using a maximum likelihood procedure based on a single rate with dated tips model and substitution rates calculated at the third codon position, evolution of the five DV-1 genotypes was shown to conform to a molecular clock. The average rate of evolution was estimated to be approximately 16.2 x 10(-4) substitutions/third codon position site/year. Using this estimate, divergence among the DV-1 genotypes was calculated to have occurred approximately 100 years ago. Very low average value of the ratio of nonsynonymous-to-synonymous nucleotide substitutions, relative to the respective sites (0.046), indicated that the evolution of the E gene of the DV-1 is subject mostly to purifying selection.

Molecular evolution of hepatitis A virus: a new classification based on the complete VP1 protein

Hepatitis A virus (HAV) is a positive-stranded RNA virus in the genus Hepatovirus in the family Picornaviridae So far, analysis of the genetic variability of HAV has been based on two discrete regions, the VP1/2A junction and the VP1 N terminus. In this report, we determined the nucleotide and deduced amino acid sequences of the complete VP1 gene of 81 strains from France, Kosovo, Mexico, Argentina, Chile, and Uruguay and compared them with the sequences of seven strains of HAV isolated elsewhere. Overall strain variation in the complete VP1 gene was found to be as high as 23.7% at the nucleotide level and 10.5% at the amino acid level. Different phylogenetic methods revealed that HAV sequences form five distinct and well-supported genetic lineages. Within these lineages, HAV sequences clustered by geographical origin only for European strains. The analysis of the complete VP1 gene allowed insight into the mode of evolution of HAV and revealed the emergence of a novel variant with a 15-amino-acid deletion located on the VP1 region where neutralization escape mutations were found. This could be the first antigenic variant of HAV so far identified.

Quantities of infectious virus and viral RNA recovered from sheep and cattle experimentally infected with foot-and-mouth disease virus O UK 2001

The profiles of virus production and excretion have been established for sheep experimentally infected with the UK 2001 strain of foot-and-mouth disease (FMD) virus by inoculation and by direct and intensive contact. Virus replicated rapidly in the inoculated sheep, from which a peak infectivity of airborne virus of 10(4.3) TCID(50) per sheep per 24 h was recovered. Around 24 h later, contact-infected sheep excreted airborne virus maximally. Similar amounts of airborne virus were recovered from cattle. The excretion of virus by the sheep under these conditions fell into three phases. First, a highly infectious period of around 7-8 days. Second, a period of 1-3 days soon afterwards when trace amounts of viral RNA were recovered in nasal and rectal swabs. Third, at 4 weeks after exposure, the demonstration, by tests on oesophageal-pharyngeal samples, that 50% of the sheep were carriers. The implications of the results and the variable role that sheep may play in the epidemiology of FMD are discussed.

Serial passage of foot-and-mouth disease virus in sheep reveals declining levels of viraemia over time

If an infectious agent is to maintain itself within a closed population by means of an unbroken serial chain of infections, it must maintain the level of infectiousness of individuals through time, or termination of the transmission chain is inevitable. One possible cause of diminution in infectiousness along serial chains of transmission may be that individuals are unable to amplify and transmit comparable levels of the infectious agent. Here, the results are reported of a novel experiment designed specifically to assess the effects of serial passage of foot-and-mouth disease virus (FMDV) in experimental groups of sheep. A virus isolate taken from an epidemic of foot-and-mouth disease (FMD) characterized by rapid fade-out of infection was passed serially through four groups of sheep housed in an isolation unit. Although it was not possible to measure individual infectiousness directly, blood virus load from infected individuals was quantified using a real-time PCR assay and used as an underlying indicator of the level of infection. The results of this assay concurred well with those of the traditional tissue-culture assay and were shown to be highly repeatable. The level of peak viraemia was shown to fall significantly with the time of infection and with passage group, both in terms of the group mean and regression analysis of individual values, suggesting that this isolate of FMDV may, under certain conditions, be unable to maintain itself indefinitely in susceptible sheep populations. The results of these experiments are discussed in terms of the epidemiology of FMD in sheep.

Aspects of the persistence of foot-and-mouth disease virus in animals--the carrier problem

Foot-and-mouth disease virus (FMDV) is a member of the Aphthovirus genus in the Picornaviridae family. Seven distinct serotypes, each including a wide range of variants, have been defined. FMD, affects wild and domesticated ruminants and pigs, is difficult to control and is the major constraint to international trade in livestock and animal products. After the acute stage of infection, FMDV may cause a prolonged, asymptomatic but persistent infection in ruminants. Also, vaccinated or naturally immune animals subsequently exposed to live virus may become persistently infected (the so-called carriers), a situation which can result in export embargoes if vaccination is included in a country's control policy.

Applications of selective neutrality tests to molecular ecology

This paper reviews how statistical tests of neutrality have been used to address questions in molecular ecology are reviewed. The work consists of four major parts: a brief review of the current status of the neutral theory; a review of several particularly interesting examples of how statistical tests of neutrality have led to insight into ecological problems; a brief discussion of the pitfalls of assuming a strictly neutral model if it is false; and a discussion of some of the opportunities and problems that molecular ecologists face when using neutrality tests to study natural selection.

Accuracy and power of bayes prediction of amino acid sites under positive selection

Bayes prediction quantifies uncertainty by assigning posterior probabilities. It was used to identify amino acids in a protein under recurrent diversifying selection indicated by higher nonsynonymous (d(N)) than synonymous (d(S)) substitution rates or by omega = d(N)/d(S) > 1. Parameters were estimated by maximum likelihood under a codon substitution model that assumed several classes of sites with different omega ratios. The Bayes theorem was used to calculate the posterior probabilities of each site falling into these site classes. Here, we evaluate the performance of Bayes prediction of amino acids under positive selection by computer simulation. We measured the accuracy by the proportion of predicted sites that were truly under selection and the power by the proportion of true positively selected sites that were predicted by the method. The accuracy was slightly better for longer sequences, whereas the power was largely unaffected by the increase in sequence length. Both accuracy and power were higher for medium or highly diverged sequences than for similar sequences. We found that accuracy and power were unacceptably low when data contained only a few highly similar sequences. However, sampling a large number of lineages improved the performance substantially. Even for very similar sequences, accuracy and power can be high if over 100 taxa are used in the analysis. We make the following recommendations: (1) prediction of positive selection sites is not feasible for a few closely related sequences; (2) using a large number of lineages is the best way to improve the accuracy and power of the prediction; and (3) multiple models of heterogeneous selective pressures among sites should be applied in real data analysis.

Codon-substitution models for detecting molecular adaptation at individual sites along specific lineages

The nonsynonymous (amino acid-altering) to synonymous (silent) substitution rate ratio (omega = d(N)/d(S)) provides a measure of natural selection at the protein level, with omega = 1, >1, and <1, indicating neutral evolution, purifying selection, and positive selection, respectively. Previous studies that used this measure to detect positive selection have often taken an approach of pairwise comparison, estimating substitution rates by averaging over all sites in the protein. As most amino acids in a functional protein are under structural and functional constraints and adaptive evolution probably affects only a few sites at a few time points, this approach of averaging rates over sites and over time has little power. Previously, we developed codon-based substitution models that allow the omega ratio to vary either among lineages or among sites. In this paper we extend previous models to allow the omega ratio to vary both among sites and among lineages and implement the new models in the likelihood framework. These models may be useful for identifying positive selection along prespecified lineages that affects only a few sites in the protein. We apply those branch-site models as well as previous branch- and site-specific models to three data sets: the lysozyme genes from primates, the tumor suppressor BRCA1 genes from primates, and the phytochrome (PHY) gene family in angiosperms. Positive selection is detected in the lysozyme and BRCA genes by both the new and the old models. However, only the new models detected positive selection acting on lineages after gene duplication in the PHY gene family. Additional tests on several data sets suggest that the new models may be useful in detecting positive selection after gene duplication in gene family evolution.

Codon-substitution models to detect adaptive evolution that account for heterogeneous selective pressures among site classes

The nonsynonymous to synonymous substitution rate ratio (omega = d(N)/d(S)) provides a sensitive measure of selective pressure at the protein level, with omega values <1, =1, and >1 indicating purifying selection, neutral evolution, and diversifying selection, respectively. Maximum likelihood models of codon substitution developed recently account for variable selective pressures among amino acid sites by employing a statistical distribution for the omega ratio among sites. Those models, called random-sites models, are suitable when we do not know a priori which sites are under what kind of selective pressure. Sometimes prior information (such as the tertiary structure of the protein) might be available to partition sites in the protein into different classes, which are expected to be under different selective pressures. It is then sensible to use such information in the model. In this paper, we implement maximum likelihood models for prepartitioned data sets, which account for the heterogeneity among site partitions by using different omega parameters for the partitions. The models, referred to as fixed-sites models, are also useful for combined analysis of multiple genes from the same set of species. We apply the models to data sets of the major histocompatibility complex (MHC) class I alleles from human populations and of the abalone sperm lysin genes. Structural information is used to partition sites in MHC into two classes: those in the antigen recognition site (ARS) and those outside. Positive selection is detected in the ARS by the fixed-sites models. Similarly, sites in lysin are classified into the buried and solvent-exposed classes according to the tertiary structure, and positive selection was detected at the solvent-exposed sites. The random-sites models identified a number of sites under positive selection in each data set, confirming and elaborating the results of the fixed-sites models. The analysis demonstrates the utility of the fixed-sites models, as well as the power of previous random-sites models, which do not use the prior information to partition sites.

Adaptive evolution in LINE-1 retrotransposons

We traced the sequence evolution of the active lineage of LINE-1 (L1) retrotransposons over the last approximately 25 Myr of human evolution. Five major families (L1PA5, L1PA4, L1PA3B, L1PA2, and L1PA1) of elements have succeeded each other as a single lineage. We found that part of the first open-reading frame (ORFI) had a higher rate of nonsynonymous (amino acid replacement) substitution than synonymous substitution during the evolution of the ancestral L1PA5 through the L1PA3B families. This segment encodes the coiled coil region of the protein-protein interaction domain of the ORFI protein (ORFIp). Statistical analysis of these changes indicates that positive selection had been acting on this region. In contrast, the coiled coil segment hardly changed during the evolution of the L1PA3B to the present L1PA1 family. Therefore, selective pressure on the coiled coil segment has changed over time. We suggest that the fast rate of amino acid replacement in the coiled coil segment reflects the adaptation of L1 either to a changing genomic environment or to host repression factors. In contrast, the second open-reading frame and the nucleic acid-binding domain of the first open-reading frame are extremely well conserved, attesting to the strong purifying selection acting on these regions.

Equine rhinitis B virus: a new serotype

Equine rhinovirus serotype 3 isolate P313/75 was assigned, with an unclassified genus status, to the family PICORNAVIRIDAE: The sequence from the 5' poly(C) tract to the 3' poly(A) tract of P313/75 was determined. The sequence is 8821 bases in length and contains a potential open reading frame for a polyprotein of 2583 amino acids. Sequence comparison and phylogenic analysis suggest that P313/75 is most closely related to the prototype equine rhinitis B virus (ERBV) strain P1436/71, formerly named equine rhinovirus type 2. A high degree of sequence similarity was found in the P2 and P3 regions of the two genomes. However, the deduced amino acid sequences of the P1 region of P313/75 and ERBV strain P1436/71 contained significant differences, which presumably account for the serological segregation of the two viruses. It is suggested that P313/75 can be classified as a new serotype of the genus Erbovirus, tentatively named ERBV2. Seroepidemiological data indicate that ERBV2 infection of horses may be common (24%) in Australia.

The generation and persistence of genetic variation in foot-and-mouth disease virus

Genetic variation in foot-and-mouth disease virus (FMDV) is of interest for at least two reasons. First, changes to the genes encoding capsid proteins results in antigenic variation, and affects vaccine efficiency and effectiveness of vaccination programs; second, genetic changes can lead to important insights into the transport of virus between countries, regions, herds, and even possibly individuals. Current estimates of RNA virus mutation rates suggest that an average of about one base mis-incorporation is likely to occur each time a single FMDV genome replicates. This should result in the introduction of every possible 1-step mutation from the progenitor genotype into the viraemia of a single infected animal many times a day. In the absence of purifying selection, a single infected animal should therefore generate a genetically very diverse population of virus.Viral-capsid sequences obtained from infected animals sampled over long-term FMDV epidemics suggest that these genetic changes accrue in a remarkably linear 'clock-like' fashion and at rates of around 1% change per year. While such a rate is generally regarded as quite high, it is actually somewhat lower than one might expect based on the rate at which viral diversity could be generated within a single animal. The difference might be explained in a variety of possible ways: (1) the mutation rate has been overestimated; (2) purifying selection is stronger than predicted; (3) only a restricted subset of excreted virus is actually infectious; (4) infected animals only excrete virus from a small partitioned subset of amplified virus, and that most of the generated viral diversity is unable to exit the animal; or (5) only a small fraction of all infected animals participate in the actual disease-transmission process.