Positive natural selection in the evolution of human metapneumovirus attachment glycoprotein

Role of human metapneumovirus glycoprotein G in modulation of immune responses

Human metapneumovirus (hMPV) is an important pathogen responsible for acute respiratory tract infections in children, the elderly, and immunocompromised patients, with no effective treatment or vaccine currently available. Knowledge of virus- and host-specific mechanisms contributing to the pathogenesis of hMPV infection is still limited. Studies have shown that hMPV surface glycoprotein G is an important virulence factor, by inhibiting innate immune signaling in airway epithelial cells and immune cells. In this study, we investigated the role of G protein in modulating innate and adaptive immune responses in mice infected with a recombinant virus with deletion of G protein (rhMPV-ΔG). Results show that rhMPV-ΔG was strongly attenuated, as it did not induce significant clinical disease, airway obstruction and airway hyperresponsiveness (AHR), compared to infection with a control strain (rhMPV-WT). By analysis of cells in bronchoalveolar fluid and lung tissue, as well as cytokine production, we found that G protein mediates aspects of both innate and adaptive immune responses, including neutrophils, dendritic cells, natural killer cells and B cells. Lung T cells recruited in response to rhMPV-ΔG had a significantly higher activated phenotype compared to those present after rhMPV-WT infection. Despite highly attenuation characterized by low levels of replication in the lung, rhMPV-ΔG was able to induce neutralizing antibodies and to protect mice from a secondary hMPV challenge. However, challenged mice that had received rhMPV-ΔG as primary infection showed some signs of lung disease at the earliest time points, which were less evident in mice that had received the rhMPV-WT strain as primary infection. These results demonstrate some of the mechanisms by which G protein could contribute to airway disease and modulate immune response to hMPV infection.

Experiments Investigating the Competitive Growth Advantage of Two Different Genotypes of Human Metapneumovirus: Implications for the Alternation of Genotype Prevalence

Human metapneumovirus (hMPV) is an important pathogen that causes upper and lower respiratory tract infections in children worldwide. hMPV has two major genotypes, hMPV-A and hMPV-B. Epidemiological studies have shown that the two hMPV genotypes alternate in predominance worldwide in recent years. Co-circulation of the two genotypes of hMPV was usually observed and there is no study about the interaction between them, such as competitive replication, which maybe the possible mechanisms for alternating prevalence of subtypes. Our present study have used two different genotypes of hMPV (genotype A: NL/1/00; B: NL/1/99) in different proportions in animal model (BALB/c mice) and cell model (Vero-E6) separately. The result showed that the competitive growth does exist in BALB/c mice, genotype B had a strong competitive advantage. However, genotype B did not cause more severe disease than non-predominant (genotype A) or mixed strains in the study, which were evaluated by the body weight, airway hyperresponsiveness and lung pathology of mouse. In cell model, competitive growth and the two genotypes alternately prevalence were observed. In summary, we confirmed that there was a competitive replication between hMPV genotype A and B, and no difference in disease severity caused by the two subtypes. This study shows a new insight to understand the alternation of hMPV genotype prevalence through genotype competition and provide experimental evidence for disease control and vaccine design.

Human metapneumovirus prevalence and patterns of subgroup persistence identified through surveillance of pediatric pneumonia hospital admissions in coastal Kenya, 2007-2016

BACKGROUND

Human metapneumovirus (HMPV) is an important respiratory pathogen that causes seasonal epidemics of acute respiratory illness and contributes significantly to childhood pneumonia. Current knowledge and understanding on its patterns of spread, prevalence and persistence in communities in low resource settings is limited.

METHODS

We present findings of a molecular-epidemiological analysis of nasal samples from children < 5 years of age admitted with syndromic pneumonia between 2007 and 2016 to Kilifi County Hospital, coastal Kenya. HMPV infection was detected using real-time RT-PCR and positives sequenced in the fusion (F) and attachment (G) genes followed by phylogenetic analysis. The association between disease severity and HMPV subgroup was assessed using Fisher's exact test.

RESULTS

Over 10 years, 274/6756 (4.1%) samples screened were HMPV positive. Annual prevalence fluctuated between years ranging 1.2 to 8.7% and lowest in the recent years (2014-2016). HMPV detections were most frequent between October of one year to April of the following year. Genotyping was successful for 205/274 (74.8%) positives revealing clades A2b (41.0%) and A2c (10.7%), and subgroups B1 (23.4%) and B2 (24.9%). The dominance patterns were: clade A2b between 2007 and 11, subgroup B1 between 2012 and 14, and clade A2c in more recent epidemics. Subgroup B2 viruses were present in all the years. Temporal phylogenetic clustering within the subgroups for both local and global sequence data was seen. Subgroups occurring in each epidemic season were comprised of multiple variants. Pneumonia severity did not vary by subgroup (p = 0.264). In both the F and G gene, the sequenced regions were found to be predominantly under purifying selection.

CONCLUSION

Subgroup patterns from this rural African setting temporally map with global strain distribution, suggesting a well-mixed global virus transmission pool of HMPV. Persistence in the local community is characterized by repeated introductions of HMPV variants from the global pool. The factors underlying the declining prevalence of HMPV in this population should be investigated.

Constraints on the Genetic and Antigenic Variability of Measles Virus

Antigenic drift and genetic variation are significantly constrained in measles virus (MeV). Genetic stability of MeV is exceptionally high, both in the lab and in the field, and few regions of the genome allow for rapid genetic change. The regions of the genome that are more tolerant of mutations (i.e., the untranslated regions and certain domains within the N, C, V, P, and M proteins) indicate genetic plasticity or structural flexibility in the encoded proteins. Our analysis reveals that strong constraints in the envelope proteins (F and H) allow for a single serotype despite known antigenic differences among its 24 genotypes. This review describes some of the many variables that limit the evolutionary rate of MeV. The high genomic stability of MeV appears to be a shared property of the Paramyxovirinae, suggesting a common mechanism that biologically restricts the rate of mutation.

Evaluating the adaptive potential of the European eel: is the immunogenetic status recovering?

The recent increased integration of evolutionary theory into conservation programs has greatly improved our ability to protect endangered species. A common application of such theory links population dynamics and indices of genetic diversity, usually estimated from neutrally evolving markers. However, some studies have suggested that highly polymorphic adaptive genes, such as the immune genes of the Major Histocompatibility Complex (MHC), might be more sensitive to fluctuations in population dynamics. As such, the combination of neutrally- and adaptively-evolving genes may be informative in populations where reductions in abundance have been documented. The European eel (Anguilla anguilla) underwent a drastic and well-reported decline in abundance in the late 20th century and still displays low recruitment. Here we compared genetic diversity indices estimated from neutral (mitochondrial DNA and microsatellites) and adaptive markers (MHC) between two distinct generations of European eels. Our results revealed a clear discrepancy between signatures obtained for each class of markers. Although mtDNA and microsatellites showed no changes in diversity between the older and the younger generations, MHC diversity revealed a contemporary drop followed by a recent increase. Our results suggest ongoing gain of MHC genetic diversity resulting from the interplay between drift and selection and ultimately increasing the adaptive potential of the species.

Cell tropism predicts long-term nucleotide substitution rates of mammalian RNA viruses

The high rates of RNA virus evolution are generally attributed to replication with error-prone RNA-dependent RNA polymerases. However, these long-term nucleotide substitution rates span three orders of magnitude and do not correlate well with mutation rates or selection pressures. This substitution rate variation may be explained by differences in virus ecology or intrinsic genomic properties. We generated nucleotide substitution rate estimates for mammalian RNA viruses and compiled comparable published rates, yielding a dataset of 118 substitution rates of structural genes from 51 different species, as well as 40 rates of non-structural genes from 28 species. Through ANCOVA analyses, we evaluated the relationships between these rates and four ecological factors: target cell, transmission route, host range, infection duration; and three genomic properties: genome length, genome sense, genome segmentation. Of these seven factors, we found target cells to be the only significant predictors of viral substitution rates, with tropisms for epithelial cells or neurons (P<0.0001) as the most significant predictors. Further, one-tailed t-tests showed that viruses primarily infecting epithelial cells evolve significantly faster than neurotropic viruses (P<0.0001 and P<0.001 for the structural genes and non-structural genes, respectively). These results provide strong evidence that the fastest evolving mammalian RNA viruses infect cells with the highest turnover rates: the highly proliferative epithelial cells. Estimated viral generation times suggest that epithelial-infecting viruses replicate more quickly than viruses with different cell tropisms. Our results indicate that cell tropism is a key factor in viral evolvability.

RNA viruses are the fastest evolving human pathogens, making their treatment and control difficult. Compared to DNA viruses, RNA viruses replicate with much lower fidelity, which can explain why RNA viruses evolve significantly faster than most DNA viruses. However, there is tremendous variation among the evolutionary rates of different RNA viruses, which is not explained by variation in mutation rates. Here we present a survey of mammalian RNA virus rates of evolution, and a comprehensive comparison of these rates to different properties of virus genomic architecture and ecology. We found that cell tropism is the most significant predictor of long-term rates of mammalian RNA virus evolution. For instance, viruses targeting epithelial cells evolve significantly faster than viruses that target neurons. Our results provide mechanistic insight into why viruses that infect respiratory and gastrointestinal epithelia have been difficult to control.

Human metapneumovirus G protein is highly conserved within but not between genetic lineages

Human metapneumovirus (HMPV) is an important cause of acute respiratory illnesses in children. HMPV encodes two major surface glycoproteins, fusion (F) and glycoprotein (G). The function of G has not been fully established, though it is dispensable for in vitro and in vivo replication. We analyzed 87 full-length HMPV G sequences from isolates collected over 20 years. The G sequences fell into four subgroups with a mean 63 % amino acid identity (minimum 29 %). The length of G varied from 217 to 241 residues. Structural features such as proline content and N- and O-glycosylation sites were present in all strains but quite variable between subgroups. There was minimal drift within the subgroups over 20 years. The estimated time to the most recent common ancestor was 215 years. HMPV G was conserved within lineages over 20 years, suggesting functional constraints on diversity. However, G was poorly conserved between subgroups, pointing to potentially distinct roles for G among different viral lineages.

Evolutionary dynamics analysis of human metapneumovirus subtype A2: genetic evidence for its dominant epidemic

Human metapneumovirus (hMPV) is a respiratory viral pathogen in children worldwide. hMPV is divided into four subtypes: hMPV_A1, hMPV_A2, hMPV_B1, and hMPV_B2. hMPV_A2 can be further divided into hMPV_A2a and A2b based on phylogenetic analysis. The typical prevalence pattern of hMPV involves a shift of the predominant subtype within one or two years. However, hMPV_A2, in particular hMPV_A2b, has circulated worldwide with a several years long term high epidemic. To study this distinct epidemic behavior of hMPV_A2, we analyzed 294 sequences of partial G genes of the virus from different countries. Molecular evolutionary data indicates that hMPV_A2 evolved toward heterogeneity faster than the other subtypes. Specifically, a Bayesian skyline plot analysis revealed that hMPV_A2 has undergone a generally upward fluctuation since 1997, whereas the other subtypes experienced only one upward fluctuation. Although hMPV_A2 showed a lower value of mean dN/dS than the other subtypes, it had the largest number of positive selection sites. Meanwhile, various styles of mutation were observed in the mutation hotspots of hMPV_A2b. Bayesian phylogeography analysis also revealed two fusions of diffusion routes of hMPV_A2b in India (June 2006) and Beijing, China (June 2008). Sequences of hMPV_A2b retrieved from GenBank boosted simultaneously with the two fusions respectively, indicating that fusion of genetic transmission routes from different regions improved survival of hMPV_A2. Epidemic and evolutionary dynamics of hMPV_A2b were similar to those of hMPV_A2. Overall, our findings provide important molecular insights into hMPV epidemics and viral variation, and explain the occurrence of an atypical epidemic of hMPV_A2, particularly hMPV_A2b.

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.

Soluble recombinant human metapneumovirus G protein is immunogenic but not protective

Human metapneumovirus (HMPV) expresses the major surface glycoproteins F and G. We evaluated the protective efficacy of immunization with G. We generated a recombinant form of G ectodomain (GDeltaTM) that was secreted from mammalian cells and purified by affinity chromatography. We tested the immunogenicity of GDeltaTM in cotton rats. Animals were immunized with PBS, GDeltaTM alone or adjuvanted, or were infected once with HMPV, and challenged with live HMPV at 28 days. Animals vaccinated with adjuvanted and non-adjuvanted GDeltaTM developed high levels of serum antibodies to both recombinant and native G protein; however, vaccinated animals did not develop neutralizing antibodies and were not protected against virus challenge. Unlike the analogous non-fusion glycoproteins of other human paramyxoviruses, HMPV G does not appear to be a protective antigen. This represents an unusual feature of HMPV.

Genetic diversity and evolution of human metapneumovirus fusion protein over twenty years

BACKGROUND

Human metapneumovirus (HMPV) is an important cause of acute respiratory illness in children. We examined the diversity and molecular evolution of HMPV using 85 full-length F (fusion) gene sequences collected over a 20-year period.

RESULTS

The F gene sequences fell into two major groups, each with two subgroups, which exhibited a mean of 96% identity by predicted amino acid sequences. Amino acid identity within and between subgroups was higher than nucleotide identity, suggesting structural or functional constraints on F protein diversity. There was minimal progressive drift over time, and the genetic lineages were stable over the 20-year period. Several canonical amino acid differences discriminated between major subgroups, and polymorphic variations tended to cluster in discrete regions. The estimated rate of mutation was 7.12 x 10(-4) substitutions/site/year and the estimated time to most recent common HMPV ancestor was 97 years (95% likelihood range 66-194 years). Analysis suggested that HMPV diverged from avian metapneumovirus type C (AMPV-C) 269 years ago (95% likelihood range 106-382 years).

CONCLUSION

HMPV F protein remains conserved over decades. HMPV appears to have diverged from AMPV-C fairly recently.

Population dynamics and rates of molecular evolution of a recently emerged paramyxovirus, avian metapneumovirus subtype C

We report the existence of two distinct sublineages of avian metapneumovirus (MPV) subtype C, a virus which has caused serious economic loss in commercial turkey farms in the United States. This subtype is closely related to human MPV, infects multiple avian species, and is globally distributed. The evolutionary rates of this virus are estimated to be 1.3 x 10(-3) to 7 x 10(-3) substitutions per site per year, and coalescent estimates place its emergence between 1991 and 1996. The four genes examined show a concordant demographic pattern which is characterized by a rapid increase in population size followed by stable population grown until the present.

Evolutionary dynamics of human and avian metapneumoviruses

Human (HMPV) and avian (AMPV) metapneumoviruses are closely related viruses that cause respiratory tract illnesses in humans and birds, respectively. Although HMPV was first discovered in 2001, retrospective studies have shown that HMPV has been circulating in humans for at least 50 years. AMPV was first isolated in the 1970s, and can be classified into four subgroups, A-D. AMPV subgroup C is more closely related to HMPV than to any other AMPV subgroup, suggesting that HMPV has emerged from AMPV-C upon zoonosis. Presently, at least four genetic lineages of HMPV circulate in human populations - A1, A2, B1 and B2 - of which lineages A and B are antigenically distinct. We used a Bayesian Markov Chain Monte Carlo (MCMC) framework to determine the evolutionary and epidemiological dynamics of HMPV and AMPV-C. The rates of nucleotide substitution, relative genetic diversity and time to the most recent common ancestor (TMRCA) were estimated using large sets of sequences of the nucleoprotein, the fusion protein and attachment protein genes. The sampled genetic diversity of HMPV was found to have arisen within the past 119-133 years, with consistent results across all three genes, while the TMRCA for HMPV and AMPV-C was estimated to have existed around 200 years ago. The relative genetic diversity observed in the four HMPV lineages was low, most likely reflecting continual population bottlenecks, with only limited evidence for positive selection.

Inhibition of Human Metapneumovirus Replication by Small Interfering RNA

Background

Human metapneumovirus (hMPV) is a major respiratory viral pathogen in young children, elderly individuals and immunocompromised patients. Despite its major effects related to bronchiolitis, pneumonia and its potential role in recurrent wheezing episodes, there is still no commercial treatment or vaccine available against this paramyxovirus.

Methods

We tested a therapeutic strategy for hMPV that was based on RNA interference.

Results

An hMPV genome-wide search for small interfering RNAs (siRNAs) by computational analysis revealed 200 potentially effective 21-mer siRNAs. Initial screening with a luciferase assay identified 57 siRNAs of interest. Further evaluation of their inhibitory potential against the four hMPV subgroups by quantitative real-time reverse transcriptase PCR and plaque immunoassay identified two highly potent siRNAs with 50% inhibitory concentration (IC50) values in the subnanomolar range. siRNA45 targets the nucleoprotein messenger RNA (mRNA) and had IC50 values <0.078 nM against representative strains from the four hMPV subgroups, whereas siRNA60, which targets the phosphoprotein mRNA, had IC50 values between 0.090– <0.078 nM against the same panel of hMPV strains. Longer 25/27-mer siRNAs known as Dicer substrates designed from the top two siRNA candidates were also evaluated and were at least as effective as their corresponding 21-mer siRNAs. Interestingly, the presence of one or two nucleotide mismatches in the target mRNA sequence of some hMPV subgroups did not always affect hMPV inhibition in vitro.

Conclusions

We successfully identified two highly efficient siRNAs against hMPV targeting essential components of the hMPV replication complex.