Respiratory syncytial virus (RSV) fusion protein subunit F2, not attachment protein G,determines the specificity of RSV infection

Human respiratory syncytial virus (HRSV) and bovine RSV (BRSV) infect human beings and cattle in a species-specific manner. We have here analyzed the contribution of RSV envelope proteins to species-specific entry into cells. In contrast to permanent cell lines, primary cells of human or bovine origin, including differentiated respiratory epithelia, peripheral blood lymphocytes, and macrophages, showed a pronounced species-specific permissiveness for HRSV and BRSV infection, respectively. Recombinant BRSV deletion mutants lacking either the small hydrophobic (SH) protein gene or both SH and the attachment glycoprotein (G) gene retained their specificity for bovine cells, whereas corresponding mutants carrying the HRSV F gene specifically infected human cells. To further narrow the responsible region of F, two reciprocal chimeric F constructs were assembled from BRSV and HRSV F1 and F2 subunits. The specificity of recombinant RSV carrying only the chimeric F proteins strictly correlated with the origin of the membrane-distal F2 domain. A contribution of G to the specificity of entry could be excluded after reintroduction of BRSV or HRSV G. Virus with F1 and G from BRSV and with only F2 from HRSV specifically infected human cells, whereas virus expressing F1 and G from HRSV and F2 from BRSV specifically infected bovine cells. The introduction of G enhanced the infectiousness of both chimeric viruses to equal degrees. Thus, the role of the nominal attachment protein G is confined to facilitating infection in a non-species-specific manner, most probably by binding to cell surface glycosaminoglycans. The identification of the F2 subunit as the determinant of RSV host cell specificity facilitates identification of virus receptors and should allow for development of reagents specifically interfering with RSV entry.     

Respiratory syncytial virus (RSV) nonstructural (NS) proteins as host range determinants: a chimeric bovine RSV with NS genes from human RSV is attenuated in interferon-competent bovine cells

Bovine respiratory syncytial virus (BRSV) escapes from cellular responses to alpha/beta interferon (IFN-alpha/beta) by a concerted action of the two viral nonstructural proteins, NS1 and NS2. Here we show that the NS proteins of human RSV (HRSV) are also able to counteract IFN responses and that they have the capacity to protect replication of an unrelated rhabdovirus. Even combinations of BRSV and HRSV NS proteins showed a protective activity, suggesting common mechanisms and cellular targets of HRSV and BRSV NS proteins. Although able to cooperate, NS proteins from BRSV and HRSV showed differential protection capacity in cells from different hosts. A chimeric BRSV with HRSV NS genes (BRSV h1/2) was severely attenuated in bovine IFN competent MDBK and Klu cells, whereas it replicated like BRSV in IFN-incompetent Vero cells or in IFN-competent human HEp-2 cells. After challenge with exogenous IFN-alpha, BRSV h1/2 was better protected than wild-type BRSV in human HEp-2 cells. In contrast, in cells of bovine origin, BRSV h1/2 was much less resistant to exogenous IFN than wild-type BRSV. These data demonstrate that RSV NS1 and NS2 proteins are major determinants of host range. The differential IFN escape capacity of RSV NS proteins in cells from different hosts provides a basis for rational development of attenuated live RSV vaccines.     

Molecular evolution and circulation patterns of human respiratory syncytial virus subgroup a: positively selected sites in the attachment g glycoprotein

Human respiratory syncytial virus (HRSV) is the most common etiological agent of acute lower respiratory tract disease in infants and can cause repeated infections throughout life. In this study, we have analyzed nucleotide sequences encompassing 629 bp at the carboxy terminus of the G glycoprotein gene for HRSV subgroup A strains isolated over 47 years, including 112 Belgian strains isolated over 19 consecutive years (1984 to 2002). By using a maximum likelihood method, we have tested the presence of diversifying selection and identified 13 positively selected sites with a posterior probability above 0.5. The sites under positive selection correspond to sites of O glycosylation or to amino acids that were previously described as monoclonal antibody-induced in vitro escape mutants. Our findings suggest that the evolution of subgroup A HRSV G glycoprotein is driven by immune pressure operating in certain codon positions located mainly in the second hypervariable region of the ectodomain. Phylogenetic analysis revealed the prolonged cocirculation of two subgroup A lineages among the Belgian population and the possible extinction of three other lineages. The evolutionary rate of HRSV subgroup A isolates was estimated to be 1.83 x 10(-3) nucleotide substitutions/site/year, projecting the most recent common ancestor back to the early 1940s.     

Chimeric bovine respiratory syncytial virus with glycoprotein gene substitutions from human respiratory syncytial virus (HRSV): effects on host range and evaluation as a live-attenuated HRSV vaccine

We recently developed a system for the generation of infectious bovine respiratory syncytial virus (BRSV) from cDNA. Here, we report the recovery of fully viable chimeric recombinant BRSVs (rBRSVs) that carry human respiratory syncytial virus (HRSV) glycoproteins in place of their BRSV counterparts, thus combining the replication machinery of BRSV with the major antigenic determinants of HRSV. A cDNA encoding the BRSV antigenome was modified so that the complete G and F genes, including the gene start and gene end signals, were replaced by their HRSV A2 counterparts. Alternatively, the BRSV F gene alone was replaced by that of HRSV Long. Each antigenomic cDNA directed the successful recovery of recombinant virus, yielding rBRSV/A2 and rBRSV/LongF, respectively. The HRSV G and F proteins or the HRSV F in combination with BRSV G were expressed efficiently in cells infected with the appropriate chimeric virus and were efficiently incorporated into recombinant virions. Whereas BRSV and HRSV grew more efficiently in bovine and human cells, respectively, the chimeric rBRSV/A2 exhibited intermediate growth characteristics in a human cell line and grew better than either parent in a bovine line. The cytopathology induced by the chimera more closely resembled that of BRSV. BRSV was confirmed to be highly restricted for replication in the respiratory tract of chimpanzees, a host that is highly permissive for HRSV. Interestingly, the rBRSV/A2 chimeric virus was somewhat more competent than BRSV for replication in chimpanzees but remained highly restricted compared to HRSV. This showed that the substitution of the G and F glycoproteins alone was not sufficient to induce efficient replication in chimpanzees. Thus, the F and G proteins contribute to the host range restriction of BRSV but are not the major determinants of this phenotype. Although rBRSV/A2 expresses the major neutralization and protective antigens of HRSV, chimpanzees infected with this chimeric virus were not significantly protected against subsequent challenge with wild-type HRSV. This suggests that the growth restriction of rBRSV/A2 was too great to provide adequate antigen expression and that the capacity of this chimeric vaccine candidate for replication in primates will need to be increased by the importation of additional HRSV genes.     

The Effect of Branch Length Variation on the Selection of Models of Molecular Evolution

Models of sequence evolution play an important role in molecular evolutionary studies. The use of inappropriate models of evolution may bias the results of the analysis and lead to erroneous conclusions. Several procedures for selecting the best-fit model of evolution for the data at hand have been proposed, like the likelihood ratio test (LRT) and the Akaike (AIC) and Bayesian (BIC) information criteria. The relative performance of these model-selecting algorithms has not yet been studied under a range of different model trees. In this study, the influence of branch length variation upon model selection is characterized. This is done by simulating sequence alignments under a known model of nucleotide substitution, and recording how often this true model is recovered by different model-fitting strategies. Results of this study agree with previous simulations and suggest that model selection is reasonably accurate. However, different model selection methods showed distinct levels of accuracy. Some LRT approaches showed better performance than the AIC or BIC information criteria. Within the LRTs, model selection is affected by the complexity of the initial model selected for the comparisons, and only slightly by the order in which different parameters are added to the model. A specific hierarchy of LRTs, which starts from a simple model of evolution, performed overall better than other possible LRT hierarchies, or than the AIC or BIC. Received: 2 October 2000 / Accepted: 4 January 2001     

Respiratory Syncytial Virus Fusion Protein Mediates Inhibition of Mitogen-Induced T-Cell Proliferation by Contact

Human respiratory syncytial virus (HRSV) and bovine respiratory syncytial virus (BRSV) are major pathogens in infants and calves, respectively. Experimental BRSV infection of calves and lambs is associated with lymphopenia and a reduction in responsiveness of peripheral blood lymphocytes (PBLs) to mitogens ex vivo. In this report, we show that in vitro mitogen-induced proliferation of PBLs is inhibited after contact with RSV-infected and UV-inactivated cells or with cells expressing RSV envelope proteins on the cell surface. The protein responsible was identified as the RSV fusion protein (F), as cells infected with a recombinant RSV expressing F as the single envelope protein or cells transfected with a plasmid encoding F were able to induce this effect. Thus, direct contact with RSV F is necessary and sufficient to inhibit proliferation of PBLs. Interestingly, F derived from HRSV was more efficient in inhibiting human PBL proliferation, while F from BRSV was more efficient in inhibiting bovine PBLs. Since various T-cell activation markers were upregulated after presenter cell contact, T lymphocytes are viable and may still be activated by mitogen. However, a significant fraction of PBLs were delayed or defective in G0/G1 to S-phase transit.     

Maximum Likelihood Estimation on Large Phylogenies and Analysis of Adaptive Evolution in Human Influenza Virus A

Algorithmic details to obtain maximum likelihood estimates of parameters on a large phylogeny are discussed. On a large tree, an efficient approach is to optimize branch lengths one at a time while updating parameters in the substitution model simultaneously. Codon substitution models that allow for variable nonsynonymous/synonymous rate ratios (ω=d _N/d _S) among sites are used to analyze a data set of human influenza virus type A hemagglutinin (HA) genes. The data set has 349 sequences. Methods for obtaining approximate estimates of branch lengths for codon models are explored, and the estimates are used to test for positive selection and to identify sites under selection. Compared with results obtained from the exact method estimating all parameters by maximum likelihood, the approximate methods produced reliable results. The analysis identified a number of sites in the viral gene under diversifying Darwinian selection and demonstrated the importance of including many sequences in the data in detecting positive selection at individual sites. Received: 25 April 2000 / Accepted: 24 July 2000     

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.     

Antiviral activity and structural characteristics of the nonglycosylated central subdomain of human respiratory syncytial virus attachment (G) glycoprotein

Segments of the cystine noose-containing nonglycosylated central subdomain, residues 149-197, of the attachment (G) glycoprotein of human respiratory syncytial virus (HRSV) have been assessed for impact on the cytopathic effect (CPE) of respiratory syncytial virus (RSV). Nalpha-acetyl residues 149-197-amide (G149-197), G149-189, and G149-177 of the A2 strain of HRSV protected 50% of human epithelial HEp-2 cells from the CPE of the A2 strain at concentrations (IC(50)) between 5 and 80 microm. Cystine noose-containing peptides G171-197 and G173-197 did not inhibit the CPE even at concentrations above 150 microm. Systematic C- and N-terminal truncations from G149-189 and G149-177 and alanine substitutions within G154-177 demonstrated that residues 166-170 (EVFNF), within a sequence that is conserved in HRSV strains, were critical for inhibition. Concordantly, G154-177 of bovine RSV and of an antibody escape mutant of HRSV with residues 166-170 of QTLPY and EVSNP, respectively, were not inhibitory. Surprisingly, a variant of G154-177 with an E166A substitution had an IC(50) of 750 nm. NMR analysis demonstrated that G149-177 adopted a well-defined conformation in solution, clustered around F168 and F170. G154-170, particularly EVFNF, may be important in binding of RSV to host cells. These findings constitute a promising platform for the development of antiviral agents for RSV.     

The Rates of Molecular Evolution in Rodent and Primate Mitochondrial DNA

A higher rate of molecular evolution in rodents than in primates at synonymous sites and, to a lesser extent, at amino acid replacement sites has been reported previously for most nuclear genes examined. Thus in these genes the average ratio of amino acid replacement to synonymous substitution rates in rodents is lower than in primates, an observation at odds with the neutral model of molecular evolution. Under Ohta's mildly deleterious model of molecular evolution, these observations are seen as the consequence of the combined effects of a shorter generation time (driving a higher mutation rate) and a larger effective population size (resulting in more effective selection against mildly deleterious mutations) in rodents. The present study reports the results of a maximum-likelihood analysis of the ratio of amino acid replacements to synonymous substitutions for genes encoded in mitochondrial DNA (mtDNA) in these two lineages. A similar pattern is observed: in rodents this ratio is significantly lower than in primates, again consistent only with the mildly deleterious model. Interestingly the lineage-specific difference is much more pronounced in mtDNA-encoded than in nuclear-encoded proteins, an observation which is shown to run counter to expectation under Ohta's model. Finally, accepting certain fossil divergence dates, the lineage-specific difference in amino acid replacement-to-synonymous substitution ratio in mtDNA can be partitioned and is found to be entirely the consequence of a higher mutation rate in rodents. This conclusion is consistent with a replication-dependent model of mutation in mtDNA. Received: 24 September 1999 / Accepted: 18 September 2000     

Structure and Origin of a Novel Dimeric Retroposon B1-dID

Here we describe a new short retroposon family of rodents. Like the primate Alu element consisting of two similar monomers, it is dimeric, but the left and right monomers are different and descend from B1 and ID short retroposons, respectively. Such elements (B1-dID) were found in the genomes of Gliridae, Sciuridae, Castoridae, Caviidae, and Hystricidae. Nucleotide sequences of this retroposon can be assigned to several structural variants. Phylogenetic analysis of B1-dID and related sequences suggests a possible scenario of B1-dID evolution in the context of rodent evolution. Received: 30 August 1999 / Accepted: 20 March 2000     

Why are translationally sub-optimal synonymous codons used in Escherichia coli?

Natural selection favors certain synonymous codons which aid translation in Escherichia coli, yet codons not favored by translational selection persist. We use the frequency distributions of synonymous polymorphisms to test three hypotheses for the existence of translationally sub-optimal codons: (1) selection is a relatively weak force, so there is a balance between mutation, selection, and drift; (2) at some sites there is no selection on codon usage, so some synonymous sites are unaffected by translational selection; and (3) translationally sub-optimal codons are favored by alternative selection pressures at certain synonymous sites. We find that when all the data is considered, model 1 is supported and both models 2 and 3 are rejected as sole explanations for the existence of translationally sub-optimal codons. However, we find evidence in favor of both models 2 and 3 when the data is partitioned between groups of amino acids and between regions of the genes. Thus, all three mechanisms appear to contribute to the existence of translationally sub-optimal codons in E. coli. Received: 18 July 2000 / Accepted: 17 April 2001     

MHC (Major Histocompatibility Complex)-DRB Genes and Polymorphisms in Common Marmoset

A New World monkey, the common marmoset (Callithrix jacchus), will be used as a preclinical animal model to study the feasibility of cell and gene therapy targeting immunological and hematological disorders. For elucidating the immunogenetic background of common marmoset to further studies, in the present study, polymorphisms of MHC-DRB genes in this species were examined. Twenty-one Caja-DRB exon 2 alleles, including seven new ones, were detected by means of subcloning and the polymerase chain reaction-single strand conformation polymorphism (PCR-SSCP) methods followed by nucleotide sequencing. Based on the alignment of these allele sequences, we designed two pairs of specific primers and established a PCR-SSCP method for DNA-based histocompatibility typing of the common marmoset. According to the family segregation data and phylogenetic analyses, we presumed that Caja-DRB alleles could be classified into five different loci. Southern blotting analysis also supported the existence of multiple DRB loci. The patterns of nucleotide substitutions suggests that positive selection operates in the antigen-recognition sites of Caja-DRB genes. Received: 18 February 2000 / Accepted: 17 May 2000     

Accumulation of Species-Specific Amino Acid Replacements That Cause Loss of Particular Protein Functions in Buchnera, an Endocellular Bacterial Symbiont

Endosymbiotic bacteria live in animal cells and are transmitted vertically at the time of the host's reproduction. In view of their small and asexual populations with infrequent chances of recombination, these endocellular bacteria are expected to accumulate mildly deleterious mutations. Previous studies showed that the DNA sequences of these bacteria evolved faster than those of free-living bacteria. In this study, we compared all the ORFs of Buchnera, an endocellular bacterial symbiont of aphids, with those of 34 other prokaryotic organisms and estimated the effect of the accelerated evolution of Buchnera on the functions of its proteins. It was revealed that Buchnera proteins contain many mutations at the sites where sequences are conserved in their orthologues in many other organisms. In addition, amino acid replacements at the conserved sites are mostly changes to physicochemically different amino acids. These results suggest that functions and conformations of Buchnera proteins have been seriously impaired or strongly modified. Indeed, extensive loss of functional motifs was observed in some Buchnera proteins. In many Buchnera proteins mutations were not detected evenly throughout each molecule but tended to accumulate in some functional units, possibly leading to loss of specific functions. As Buchnera has an unusual and limited gene repertory, it is conceivable that the manner of interactions among its proteins has been changed, and thus, functional constraints over their amino acid residues have also been changed during evolution. This may account for the loss of some functional units only in the Buchnera proteins. We obtained evidence that amino acid replacements in Buchnera were not always deleterious, but neutral or, in some cases, even positively selected. Received: 14 December 2000 / Accepted: 12 March 2001     

Bovine respiratory syncytial virus protects cotton rats against human respiratory syncytial virus infection.

Human respiratory syncytial virus (HRSV) is the most frequent cause of severe respiratory infections in infancy. No vaccine against this virus has yet been protective, and antiviral drugs have been of limited utility. Using the cotton rat model of HRSV infection, we examined bovine respiratory syncytial virus (BRSV), a cause of acute respiratory disease in young cattle, as a possible vaccine candidate to protect children against HRSV infection. Cotton rats were primed intranasally with graded doses of BRSV/375 or HRSV/Long or were left unprimed. Three weeks later, they were challenged intranasally with either BRSV/375, HRSV/Long (subgroup A), or HRSV/18537 (subgroup B). At intervals postchallenge, animals were sacrificed for virus titration and histologic evaluation. Serum neutralizing antibody titers were determined at the time of viral challenge. BRSV/375 replicated to low titers in nasal tissues and lungs. Priming with 10(5) PFU of BRSV/375 effected a 500- to 1,000-fold reduction in peak nasal HRSV titer and a greater than 1,000-fold reduction in peak pulmonary HRSV titer upon challenge with HRSV/Long or HRSV/18537. In contrast to priming with HRSV, priming with BRSV did not induce substantial levels of neutralizing antibody against HRSV and was associated with a delayed onset of clearance of HRSV upon challenge. Priming with BRSV/375 caused mild nasal and pulmonary pathology and did not cause exacerbation of disease upon challenge with HRSV/Long. Our findings suggest that BRSV may be a potential vaccine against HRSV and a useful tool for studying the mechanisms of immunity to HRSV.     

Molecular Evolution in a Multisite Nearly Neutral Mutation Model

A simple nearly neutral mutation model of protein evolution was studied using computer simulation assuming a constant population size. In this model, a gene consists of a finite number of codons and there is no recombination within a gene. Each codon has two replacement and one silent sites. The fitness of a gene was determined multiplicatively by amino acids specified by codons (the independent multicodon model). Nucleotide diversity at replacement sites decreases as selection becomes stronger. A reduction of nucleotide diversity at silent sites also occurs as selection intensifies but the magnitude of the reduction is not a monotone function of the intensity of selection. The dispersion index is close to one. The average value of Tajima's and Fu and Li's statistics are negative and their absolute values increases as selection intensifies. However, their powers of detecting selection under the present model were not high unless the number of sites is large or mutation rate is high. The MK test was shown to detect intermediate selection fairly well. For comparison, the house-of-cards model was also investigated and its behavior was shown to be more sensitive to changes of population size than that of the independent multicodon model. The relevance of the present model for explaining protein evolution was discussed comparing its prediction and recent DNA data. Received: 24 May 1999 / Accepted: 17 August 1999     

Positive Selection Results in Frequent Reversible Amino Acid Replacements in the G Protein Gene of Human Respiratory Syncytial Virus

Human respiratory syncytial virus (HRSV) is the major cause of lower respiratory tract infections in children under 5 years of age and the elderly, causing annual disease outbreaks during the fall and winter. Multiple lineages of the HRSVA and HRSVB serotypes co-circulate within a single outbreak and display a strongly temporal pattern of genetic variation, with a replacement of dominant genotypes occurring during consecutive years. In the present study we utilized phylogenetic methods to detect and map sites subject to adaptive evolution in the G protein of HRSVA and HRSVB. A total of 29 and 23 amino acid sites were found to be putatively positively selected in HRSVA and HRSVB, respectively. Several of these sites defined genotypes and lineages within genotypes in both groups, and correlated well with epitopes previously described in group A. Remarkably, 18 of these positively selected tended to revert in time to a previous codon state, producing a “flip-flop” phylogenetic pattern. Such frequent evolutionary reversals in HRSV are indicative of a combination of frequent positive selection, reflecting the changing immune status of the human population, and a limited repertoire of functionally viable amino acids at specific amino acid sites.     

Phylogenomic Analysis of the α Proteasome Gene Family from Early-Diverging Eukaryotes

We employed a phylogenomic approach to study the evolution of α subunits of the proteasome gene family from early diverging eukaryotes. BLAST similarity searches of the Giardia lamblia genome identified all seven α proteasome genes characteristic of eukaryotes from the crown group. In addition, a PCR strategy for the amplification of multiple α subunit sequences generated single α proteasome products for representatives of the Kinetoplastida (Leishmania major), the Parabasalia (Trichomonas vaginalis), and the Microsporidia (Vairimorpha sp., Nosema sp., Endoreticulata sp., and Spraguea lophii). The kinetoplastid Trypanosoma cruzi and the eukaryote crown group Acanthamoeba castellanii yielded two distinct α proteasome genes each. The presence of seven distinct α proteasome genes in G. lamblia, one of the earliest-diverging eukaryotes, indicates that the α proteasome gene family evolved rapidly from a minimum of one gene in Archaea to seven or more in Eukarya. Results from the phylogenomic analysis are consistent with the idea that the Diplomonida (as represented by G. lamblia), the Kinetoplastida, the Parabasalia, and the Microsporidia diverged after the duplication events that originated the α proteasome gene family. A model for the early origin and evolution of the proteasome gene family is presented. Received: 14 February 2000 / Accepted: 14 August 2000     

Origin and evolution of GBV-C/hepatitis G virus and relationships with ancient human migrations

The GB virus C/hepatitis G virus (GBV-C/HGV) is a newly identified human RNA virus, belonging to the Flaviviridae family. Persistent infection by GBV-C/HGV is common in humans, and genetically divergent isolates have been identified in different parts of the world. Due to the absence of a real pathogenic role of GBV-C/HGV in liver disease and its extremely low mutation rate, this virus is a potential marker to trace prehistoric links between human populations. In this study, origin and evolution of GBV-C/HGV were examined using a set of fully sequenced strains of worldwide origin. A first phylogenetic analysis, addressed to the short (255 nucleotides) NS5A overlapping coding region by the neighbor-joining method, suggested an ancient African origin of GBV-C/HGV. This notion was confirmed when the same analysis was applied to the genomic regions showing the lowest rate of synonymous substitutions, covering one-fourth (2184 nucleotides) of the total coding potential of the virus genome. By using a multivariate statistical method and extending the analysis to the complete coding region, fine details of the evolutionary history of GBV-C/HGV were further elucidated. By this approach, isolates from Southeast Asia appeared to be the most closely related to those of African origin, consistent with a major route of ancient human migrations from Africa to southeastern parts of the Asian continent. Received: 26 October 2000 / Accepted: 28 February 2001