Molecular epidemiology and evolution of human respiratory syncytial virus and human metapneumovirus

Human respiratory syncytial virus (HRSV) and human metapneumovirus (HMPV) are ubiquitous respiratory pathogens of the Pneumovirinae subfamily of the Paramyxoviridae. Two major surface antigens are expressed by both viruses; the highly conserved fusion (F) protein, and the extremely diverse attachment (G) glycoprotein. Both viruses comprise two genetic groups, A and B. Circulation frequencies of the two genetic groups fluctuate for both viruses, giving rise to frequently observed switching of the predominantly circulating group. Nucleotide sequence data for the F and G gene regions of HRSV and HMPV variants from the UK, The Netherlands, Bangkok and data available from Genbank were used to identify clades of both viruses. Several contemporary circulating clades of HRSV and HMPV were identified by phylogenetic reconstructions. The molecular epidemiology and evolutionary dynamics of clades were modelled in parallel. Times of origin were determined and positively selected sites were identified. Sustained circulation of contemporary clades of both viruses for decades and their global dissemination demonstrated that switching of the predominant genetic group did not arise through the emergence of novel lineages each respiratory season, but through the fluctuating circulation frequencies of pre-existing lineages which undergo proliferative and eclipse phases. An abundance of sites were identified as positively selected within the G protein but not the F protein of both viruses. For HRSV, these were discordant with previously identified residues under selection, suggesting the virus can evade immune responses by generating diversity at multiple sites within linear epitopes. For both viruses, different sites were identified as positively selected between genetic groups.     

Low-pH-induced membrane fusion mediated by human metapneumovirus F protein is a rare, strain-dependent phenomenon

Membrane fusion promoted by human metapneumovirus (HMPV) fusion (F) protein was suggested to require low pH (R. M. Schowalter, S. E. Smith, and R. E. Dutch, J. Virol. 80:10931-10941, 2006). Using prototype F proteins representing the four HMPV genetic lineages, we detected low-pH-dependent fusion only with some lineage A proteins and not with lineage B proteins. A glycine at position 294 was found responsible for the low-pH requirement in lineage A proteins. Only 6% of all HMPV lineage A F sequences have 294G, and none of the lineage B sequences have 294G. Thus, acidic pH is not a general trigger of HMPV F proteins for activity.     

Recombinant human Metapneumovirus lacking the small hydrophobic SH and/or attachment G glycoprotein: deletion of G yields a promising vaccine candidate

Human metapneumovirus (HMPV) has recently been identified as a significant cause of serious respiratory tract disease in humans. In particular, the emerging information on the contribution of HMPV to pediatric respiratory tract disease suggests that it will be important to develop a vaccine against this virus for use in conjunction with those being developed for human respiratory syncytial virus and the human parainfluenza viruses. A recently described reverse genetic system (S. Biacchesi, M. H. Skiadopoulos, K. C. Tran, B. R. Murphy, P. L. Collins, and U. J. Buchholz, Virology 321:247-259, 2004) was used to generate recombinant HMPVs (rHMPVs) that lack the G gene, the SH gene, or both. The DeltaSH, DeltaG, and DeltaSH/G deletion mutants were readily recovered and were found to replicate efficiently during multicycle growth in cell culture. Thus, the SH and G proteins are not essential for growth in cell culture. Apart from the absence of the deleted protein(s), the virions produced by the gene deletion mutants were similar by protein yield and gel electrophoresis protein profile to wild-type HMPV. When administered intranasally to hamsters, the DeltaG and DeltaSH/G mutants replicated in both the upper and lower respiratory tracts, showing that HMPV containing F as the sole viral surface protein is competent for replication in vivo. However, both viruses were at least 40-fold and 600-fold restricted in replication in the lower and upper respiratory tract, respectively, compared to wild-type rHMPV. They also induced high titers of HMPV-neutralizing serum antibodies and conferred complete protection against replication of wild-type HMPV challenge virus in the lungs. Surprisingly, G is dispensable for protection, and the DeltaG and DeltaSH/G viruses represent promising vaccine candidates. In contrast, DeltaSH replicated somewhat more efficiently in hamster lungs compared to wild-type rHMPV (20-fold increase on day 5 postinfection). This indicates that SH is completely dispensable in vivo and that its deletion does not confer an attenuating effect, at least in this rodent model.     

Modification of the trypsin-dependent cleavage activation site of the human metapneumovirus fusion protein to be trypsin independent does not increase replication or spread in rodents or nonhuman primates

The contribution of cleavage activation of the fusion F protein of human metapneumovirus (HMPV) to replication and pathogenicity in rodents and nonhuman primates was investigated. Recombinant HMPVs were generated in which the naturally occurring trypsin-dependent cleavage sequence (R-Q-S-R downward arrow) was replaced by each of three sequences whose cleavage in vitro does not depend upon added trypsin. Two of these were multibasic sequences derived from avian metapneumovirus type A (R-R-R-R) or type C (R-K-A-R), with the former containing the consensus furin protease cleavage motif (R-X-R/K-R downward arrow). The third one (R-Q-P-R) was derived from a recently described trypsin independent HMPV isolate (J. H. Schickli, J. Kaur, N. Ulbrandt, R. R. Spaete, and R. S. Tang, J. Virol. 79:10678-10689, 2005). To preclude the possibility of conferring even greater virulence to this significant human pathogen, the modifications were done in an HMPV variant that was attenuated by the deletion of two of the three envelope glycoproteins, SH and G. Each of the introduced cleavage sequences conferred trypsin independent F cleavage and growth to HMPV in vitro. However, they differed in the efficiency of trypsin independent growth and plaque formation in vitro: R-R-R-R > R-K-A-R > R-Q-P-R > R-Q-S-R. The R-R-R-R mutant was the only one whose growth in vitro was not augmented by added trypsin, indicative of highly efficient trypsin independent cleavage. When inoculated intranasally into hamsters, there was essentially no difference in the magnitude of replication in the upper or lower respiratory tract between the mutants, and virus was not detected in organs outside of the respiratory tract. Evaluation of the most cleavage-efficient mutant, R-R-R-R, in African green monkeys showed that there was no detectable change in the magnitude of replication in the upper and lower respiratory tract or in immunogenicity and protective efficacy against HMPV challenge. These results suggest that cleavage activation is not a major determinant of HMPV virulence.     

Differential modes of MHC class IIB gene evolution in cichlid fishes

Cichlid fishes are emblematic models for the study of adaptive radiation, driven by natural and sexual selection. Parasite mediated selection is an important component in these processes, and the evolution of their immune system therefore merits special attention. In this study, light is shed on the phylogeny of the b family of cichlid major histocompatibility complex (MHC) class IIB genes. Full-length coding sequences were used to reconstruct phylogenies using criteria of maximum parsimony, maximum likelihood and Bayesian inference. All analyses suggest monophyly of the b family of cichlid MHC class IIB genes, although sequences of the cichlid sister taxa are currently not available. Two evolutionary lineages of these genes, respectively encompassing the recently defined genomic regions DBB-DEB-DFB and DCB-DDB, show highly contrasting levels of differentiation. To explore putative causes for these differences, exon 2 sequences were screened for variation in recombination rate and strength of selection. The more diversified lineage of cichlid MHC class IIB b genes was found to have higher levels of both recombination and selection. This is consistent with the observation in other taxa that recombination facilitates the horizontal spread of positively selected sites across MHC loci and hence contributes to fast sequence evolution. In contrast, the lineage that showed low diversification might either be under stabilizing selection or is evolutionary constrained by its low recombination rate. We speculate whether this lineage might include MHC genes with non-classical functions.     

Comparative analysis of evolutionary mechanisms of the hemagglutinin and three internal protein genes of influenza B virus: multiple cocirculating lineages and frequent reassortment of the NP, M, and NS genes

Phylogenetic profiles of the genes coding for the hemagglutinin (HA) protein, nucleoprotein (NP), matrix (M) protein, and nonstructural (NS) proteins of influenza B viruses isolated from 1940 to 1998 were analyzed in a parallel manner in order to understand the evolutionary mechanisms of these viruses. Unlike human influenza A (H3N2) viruses, the evolutionary pathways of all four genes of recent influenza B viruses revealed similar patterns of genetic divergence into two major lineages. Although evolutionary rates of the HA, NP, M, and NS genes of influenza B viruses were estimated to be generally lower than those of human influenza A viruses, genes of influenza B viruses demonstrated complex phylogenetic patterns, indicating alternative mechanisms for generation of virus variability. Topologies of the evolutionary trees of each gene were determined to be quite distinct from one another, showing that these genes were evolving in an independent manner. Furthermore, variable topologies were apparently the result of frequent genetic exchange among cocirculating epidemic viruses. Evolutionary analysis done in the present study provided further evidence for cocirculation of multiple lineages as well as sequestering and reemergence of phylogenetic lineages of the internal genes. In addition, comparison of deduced amino acid sequences revealed a novel amino acid deletion in the HA1 domain of the HA protein of recent isolates from 1998 belonging to the B/Yamagata/16/88-like lineage. It thus became apparent that, despite lower evolutionary rates, influenza B viruses were able to generate genetic diversity among circulating viruses through a combination of evolutionary mechanisms involving cocirculating lineages and genetic reassortment by which new variants with distinct gene constellations emerged.     

The two major human metapneumovirus genetic lineages are highly related antigenically, and the fusion (F) protein is a major contributor to this antigenic relatedness

The growth properties and antigenic relatedness of the CAN98-75 (CAN75) and the CAN97-83 (CAN83) human metapneumovirus (HMPV) strains, which represent the two distinct HMPV genetic lineages and exhibit 5 and 63% amino acid divergence in the fusion (F) and attachment (G) proteins, respectively, were investigated in vitro and in rodents and nonhuman primates. Both strains replicated to high titers (> or =6.0 log(10)) in the upper respiratory tract of hamsters and to moderate titers (> or =3.6 log(10)) in the lower respiratory tract. The two lineages exhibited 48% antigenic relatedness based on reciprocal cross-neutralization assay with postinfection hamster sera, and infection with each strain provided a high level of resistance to reinfection with the homologous or heterologous strain. Hamsters immunized with a recombinant human parainfluenza virus type 1 expressing the fusion F protein of the CAN83 strain developed a serum antibody response that efficiently neutralized virus from both lineages and were protected from challenge with either HMPV strain. This result indicates that the HMPV F protein is a major antigenic determinant that mediates extensive cross-lineage neutralization and protection. Both HMPV strains replicated to low titers in the upper and lower respiratory tracts of rhesus macaques but induced high levels of HMPV-neutralizing antibodies in serum effective against both lineages. The level of HMPV replication in chimpanzees was moderately higher, and infected animals developed mild colds. HMPV replicated the most efficiently in the respiratory tracts of African green monkeys, and the infected animals developed a high level of HMPV serum-neutralizing antibodies (1:500 to 1:1,000) effective against both lineages. Reciprocal cross-neutralization assays in which postinfection sera from all three primate species were used indicated that CAN75 and CAN83 are 64 to 99% related antigenically. HMPV-infected chimpanzees and African green monkeys were highly protected from challenge with the heterologous HMPV strain. Taken together, the results from hamsters and nonhuman primates support the conclusion that the two HMPV genetic lineages are highly related antigenically and are not distinct antigenic subtypes or subgroups as defined by reciprocal cross-neutralization in vitro.     

Male and female recombination landscapes of diploid Arabidopsis arenosa

The number and placement of meiotic crossover events during meiosis have important implications for the fidelity of chromosome segregation as well as patterns of inheritance. Despite the functional importance of recombination, recombination landscapes vary widely among and within species, and this can have a strong impact on evolutionary processes. A good knowledge of recombination landscapes is important for model systems in evolutionary and ecological genetics, since it can improve interpretation of genomic patterns of differentiation and genome evolution, and provides an important starting point for understanding the causes and consequences of recombination rate variation. Arabidopsis arenosa is a powerful evolutionary genetic model for studying the molecular basis of adaptation and recombination rate evolution. Here, we generate genetic maps for 2 diploid A. arenosa individuals from distinct genetic lineages where we have prior knowledge that meiotic genes show evidence of selection. We complement the genetic maps with cytological approaches to map and quantify recombination rates, and test the idea that these populations might have distinct patterns of recombination. We explore how recombination differs at the level of populations, individuals, sexes and genomic regions. We show that the positioning of crossovers along a chromosome correlates with their number, presumably a consequence of crossover interference, and discuss how this effect can cause differences in recombination landscape among sexes or species. We identify several instances of female segregation distortion. We found that averaged genome-wide recombination rate is lower and sex differences subtler in A. arenosa than in Arabidopsis thaliana.     

Molecular Analysis of a Salmonella Enterica Group E1 Rfb Gene Cluster: O Antigen and the Genetic Basis of the Major Polymorphism

Salmonella enterica is highly polymorphic for the O antigen, a surface polysaccharide that is subject to intense selection by the host immune system. This polymorphism is used for serotyping Salmonella isolates. The genes encoding O antigen biosynthesis are located in the rfb gene cluster. We report here the cloning and sequence of the 19-kb rfb region from strain M32 (serovar anatum, group E1) and compare it with that of strain LT2 (serovar typhimurium, group B). Genes for biosynthetic pathways common to both strains are conserved and have very similar sequences. In contrast, the five genes for CDP-abequose synthesis, present in strain LT2, are absent in strain M32; three open reading frames (ORFs) of strain LT2, thought to include genes for transferases, are not present in strain M32 but are replaced by three different ORFs with little or low level of similarity. Both rfb gene clusters are low in G + C content, indicating that they were transferred from a common ancestral species with low G + C content to S. enterica relatively recently (in the evolutionary sense). We discuss the recombination and lateral transfer events which may have been involved in the evolution of the polymorphism.     

Genome-Wide Analyses of Recombination Suggest That Giardia intestinalis Assemblages Represent Different Species

Giardia intestinalis is a major cause of waterborne enteric disease in humans. The species is divided into eight assemblages suggested to represent separate Giardia species based on host specificities and the genetic divergence of marker genes. We have investigated whether genome-wide recombination occurs between assemblages using the three available G. intestinalis genomes. First, the relative nonsynonymous substitution rates of the homologs were compared for 4,009 positional homologs. The vast majority of these comparisons indicate genetic isolation without interassemblage recombinations. Only a region of 6 kbp suggests genetic exchange between assemblages A and E, followed by gene conversion events. Second, recombination-detecting software fails to identify within-gene recombination between the different assemblages for most of the homologs. Our results indicate very low frequency of recombination between the syntenic core genes, suggesting that G. intestinalis assemblages are genetically isolated lineages and thus should be viewed as separated Giardia species.     

Evolutionary dynamics of Ralstonia solanacearum

We investigated the genetic diversity, extent of recombination, natural selection, and population divergence of Ralstonia solanacearum samples obtained from sources worldwide. This plant pathogen causes bacterial wilt in many crops and constitutes a serious threat to agricultural production due to its very wide host range and aggressiveness. Five housekeeping genes, dispersed around the chromosome, and three virulence-related genes, located on the megaplasmid, were sequenced from 58 strains belonging to the four major phylogenetic clusters (phylotypes). Whereas genetic variation is high and consistent for all housekeeping loci studied, virulence-related gene sequences are more diverse. Phylogenetic and statistical analyses suggest that this organism is a highly diverse bacterial species containing four major, deeply separated evolutionary lineages (phylotypes I to IV) and a weaker subdivision of phylotype II into two subgroups. Analysis of molecular variations showed that the geographic isolation and spatial distance have been the significant determinants of genetic variation between phylotypes. R. solanacearum displays high clonality for housekeeping genes in all phylotypes (except phylotype III) and significant levels of recombination for the virulence-related egl and hrpB genes, which are limited mainly to phylotype strains III and IV. Finally, genes essential for species survival are under purifying selection, and those directly involved in pathogenesis might be under diversifying selection.     

Infection of nonhuman primates with recombinant human metapneumovirus lacking the SH, G, or M2-2 protein categorizes each as a nonessential accessory protein and identifies vaccine candidates

Recombinant human metapneumovirus (HMPV) in which the SH, G, or M2 gene or open reading frame was deleted by reverse genetics was evaluated for replication and vaccine efficacy following topical administration to the respiratory tract of African green monkeys, a permissive primate host. Replication of the deltaSH virus was only marginally less efficient than that of wild-type HMPV, whereas the deltaG and deltaM2-2 viruses were reduced sixfold and 160-fold in the upper respiratory tract and 3,200-fold and 4,000-fold in the lower respiratory tract, respectively. Even with the highly attenuated mutants, there was unequivocal HMPV replication at each anatomical site in each animal. Thus, none of these three proteins is essential for HMPV replication in a primate host, although G and M2-2 increased the efficiency of replication. Each gene-deletion virus was highly immunogenic and protective against wild-type HMPV challenge. The deltaG and deltaM2-2 viruses are promising vaccine candidates that are based on independent mechanisms of attenuation and are appropriate for clinical evaluation.     

Contrasting recombination patterns and demographic histories of the plant pathogen Ralstonia solanacearum inferred from MLSA

We used multilocus sequence analysis (MLSA) on a worldwide collection of the plant pathogenic Ralstonia solanacearum (Betaproteobacteria) to retrace its complex evolutionary history. Using genetic imprints left during R. solanacearum evolution, we were able to delineate distinct evolutionary complex displaying contrasting dynamics. Among the phylotypes already described (I, IIA, IIB, III, IV), eight groups of strains with distinct evolutionary patterns, named clades, were identified. From our recombination analysis, we identified 21 recombination events that occurred within and across these lineages. Although appearing the most divergent and ancestral phylotype, phylotype IV was inferred as a gene donor for the majority of the recombination events that we detected. Whereas this phylotype apparently fuelled the species diversity, ongoing diversification was mainly detected within phylotype I, IIA and III. These three groups presented a recent expanding population structure, a high level of homologous recombination and evidences of long-distance migrations. Factors such as adaptation to a specific host or intense trading of infected crops may have promoted this diversification. Whether R. solanacearum lineages will eventually evolve in distinct species remains an open question. The intensification of cropping and increase of geographical dispersion may favour situations of phylotype sympatry and promote higher exchange of key factors for host adaptation from their common genetic pool.     

Sex‐antagonistic genes,XY recombination and feminized Y chromosomes

The canonical model of sex‐chromosome evolution predicts that sex‐antagonistic (SA) genes play an instrumental role in the arrest of XY recombination and ensuing Y chromosome degeneration. Although this model might account for the highly differentiated sex chromosomes of birds and mammals, it does not fit the situation of many lineages of fish, amphibians or nonavian reptiles, where sex chromosomes are maintained homomorphic through occasional XY recombination and/or high turnover rates. Such situations call for alternative explanatory frameworks. A crucial issue at stake is the effect of XY recombination on the dynamics of SA genes and deleterious mutations. Using individual‐based simulations, we show that a complete arrest of XY recombination actually benefits females, not males. Male fitness is maximized at different XY recombination rates depending on SA selection, but never at zero XY recombination. This should consistently favour some level of XY recombination, which in turn generates a recombination load at sex‐linked SA genes. Hill–Robertson interferences with deleterious mutations also impede the differentiation of sex‐linked SA genes, to the point that males may actually fix feminized phenotypes when SA selection and XY recombination are low. We argue that sex chromosomes might not be a good localization for SA genes, and sex conflicts seem better solved through the differential expression of autosomal genes.     

Evolutionary Dynamics of Ralstonia solanacearum

We investigated the genetic diversity, extent of recombination, natural selection, and population divergence of Ralstonia solanacearum samples obtained from sources worldwide. This plant pathogen causes bacterial wilt in many crops and constitutes a serious threat to agricultural production due to its very wide host range and aggressiveness. Five housekeeping genes, dispersed around the chromosome, and three virulence-related genes, located on the megaplasmid, were sequenced from 58 strains belonging to the four major phylogenetic clusters (phylotypes). Whereas genetic variation is high and consistent for all housekeeping loci studied, virulence-related gene sequences are more diverse. Phylogenetic and statistical analyses suggest that this organism is a highly diverse bacterial species containing four major, deeply separated evolutionary lineages (phylotypes I to IV) and a weaker subdivision of phylotype II into two subgroups. Analysis of molecular variations showed that the geographic isolation and spatial distance have been the significant determinants of genetic variation between phylotypes. R. solanacearum displays high clonality for housekeeping genes in all phylotypes (except phylotype III) and significant levels of recombination for the virulence-related egl and hrpB genes, which are limited mainly to phylotype strains III and IV. Finally, genes essential for species survival are under purifying selection, and those directly involved in pathogenesis might be under diversifying selection.     

Recombination confounds the early evolutionary history of human immunodeficiency virus type 1: subtype G is a circulating recombinant form

Human immunodeficiency virus type 1 (HIV-1) is classified in nine subtypes (A to D, F, G, H, J, and K), a number of subsubtypes, and several circulating recombinant forms (CRFs). Due to the high level of genetic diversity within HIV-1 and to its worldwide distribution, this classification system is widely used in fields as diverse as vaccine development, evolution, epidemiology, viral fitness, and drug resistance. Here, we demonstrate how the high recombination rates of HIV-1 may confound the study of its evolutionary history and classification. Our data show that subtype G, currently classified as a pure subtype, has in fact a recombinant history, having evolved following recombination between subtypes A and J and a putative subtype G parent. In addition, we find no evidence for recombination within one of the lineages currently classified as a CRF, CRF02_AG. Our analysis indicates that CRF02_AG was the parent of the recombinant subtype G, rather than the two having the opposite evolutionary relationship, as is currently proposed. Our results imply that the current classification of HIV-1 subtypes and CRFs is an artifact of sampling history, rather than reflecting the evolutionary history of the virus. We suggest a reanalysis of all pure subtypes and CRFs in order to better understand how high rates of recombination have influenced HIV-1 evolutionary history.     

The Evolutionary Dynamics of Human Influenza B Virus

Despite their close phylogenetic relationship, type A and B influenza viruses exhibit major epidemiological differences in humans, with the latter both less common and less often associated with severe disease. However, it is unclear what processes determine the evolutionary dynamics of influenza B virus, and how influenza viruses A and B interact at the evolutionary scale. To address these questions we inferred the phylogenetic history of human influenza B virus using complete genome sequences for which the date (day) of isolation was available. By comparing the phylogenetic patterns of all eight viral segments we determined the occurrence of segment reassortment over a 30-year sampling period. An analysis of rates of nucleotide substitution and selection pressures revealed sporadic occurrences of adaptive evolution, most notably in the viral hemagglutinin and compatible with the action of antigenic drift, yet lower rates of overall and nonsynonymous nucleotide substitution compared to influenza A virus. Overall, these results led us to propose a model in which evolutionary changes within and between the antigenically distinct 'Yam88' and 'Vic87' lineages of influenza B virus are the result of changes in herd immunity, with reassortment continuously generating novel genetic variation. Additionally, we suggest that the interaction with influenza A virus may be central in shaping the evolutionary dynamics of influenza B virus, facilitating the shift of dominance between the Vic87 and the Yam88 lineages.