Molecular epidemiology and phylodynamics of the human respiratory syncytial virus fusion protein in northern Taiwan

Discovery and characterization of BRBV-sheep virus in nasal swabs from domestic sheep in China

Introduction

The escalating occurrence of infectious disease outbreaks in humans and animals necessitates innovative, effective, and integrated research to better comprehend their transmission and dynamics. Viral infection in livestock has led to profound economic losses globally. Pneumonia is the prevalent cause of death in sheep. However, very few studies exist regarding virus-related pathogens in sheep. Metagenomics sequencing technologies in livestock research hold significant potential to elucidate these contingencies and enhance our understanding.

Methods

Therefore, this study aims to characterize respiratory viromes in paired nasal swabs from Inner Mongolian feedlot sheep in China using metaviromic sequencing. Through deep sequencing, de novo assembly, and similarity searches using translated protein sequences, several previously uncharacterized and known viruses were identified in this study.

Results

Among these discoveries, a novel Bovine Rhinitis B Virus (BRBV) (BRBV-sheep) strain was serendipitously detected in the nasal swabs of domestic sheep (Ovis aries). To facilitate further molecular epidemiological studies, the entire genome of BRBV-sheep was also determined. Owing to the unique sequence characteristics and phylogenetic position of BRBV-sheep, genetically distinct lineages of BRBV in sheep may exist. A TaqMan-based qRT-PCR assay targeting the 3D polymerase gene was developed and used to screen 592 clinical sheep specimens. The results showed that 44.59% of the samples (264/592) were positive. These findings suggest that BRBV sheep are widespread among Inner Mongolian herds.

Conclusion

This discovery marks the initial identification of BRBV in sheep within Inner Mongolia, China. These findings contribute to our understanding of the epidemiology and genetic evolution of BRBV. Recognizing the presence of BRBV in sheep informs strategies for disease management and surveillance and the potential development of targeted interventions to control its spread.

Genomic evolution of human respiratory syncytial virus during a decade (2013-2023): bridging the path to monoclonal antibody surveillance

OBJECTIVES

This study investigated the prevalence, genetic diversity, and evolution of human respiratory syncytial virus (HRSV) in Barcelona from 2013 to 2023.

METHODS

Respiratory specimens from patients with RTI suspicion at Hospital Universitari Vall d'Hebron were collected from October 2013 to May 2023 for laboratory-confirmation of respiratory viruses. Next-generation sequencing was performed in randomly-selected samples with Illumina technology. Phylogenetic analyses of whole genome sequences were performed with BEAST v1.10.4. Signals of selection and evolutionary pressures were inferred by population dynamics and evolutionary analyses. Mutations in major surface proteins were genetic and structurally characterised, emphasizing those within antigenic epitopes.

RESULTS

Analyzing 139,625 samples, 5.3% were HRSV-positive (3008 HRSV-A, 3882 HRSV-B, 56 HRSV-A and -B, and 495 unsubtyped HRSV), with a higher prevalence observed in the paediatric population. Pandemic-related shifts in seasonal patterns returned to normal in 2022-2023. A total of 198 whole-genome sequences were obtained for HRSV-A (6.6% of the HRSV-A positive samples) belonging to GA2.3.5 lineage. For HRSV-B, 167 samples were sequenced (4.3% of the HRSV-B positive samples), belonging to GB5.0.2, GB5.0.4a and GB5.0.5a. HRSV-B exhibited a higher evolution rate. Post-SARS-CoV-2 pandemic, both subtypes showed increased evolutionary rates and decreased effective population size initially, followed by a sharp increase. Analyses indicated negative selective pressure on HRSV. Mutations in antigenic epitopes, including S276N and M274I in palivizumab-targeted site II, and I206M, Q209R, and S211N in nirsevimab-targeted site Ø, were identified.

DISCUSSION

Particularly in the context of the large-scale use in 2023-2024 season of nirsevimab, continuous epidemiological and genomic surveillance is crucial.

Detailed Evolutionary Analyses of the F Gene in the Respiratory Syncytial Virus Subgroup A

We performed evolution, phylodynamics, and reinfection-related antigenicity analyses of respiratory syncytial virus subgroup A (RSV-A) fusion (F) gene in globally collected strains (1465 strains) using authentic bioinformatics methods. The time-scaled evolutionary tree using the Bayesian Markov chain Monte Carlo method estimated that a common ancestor of the RSV-A, RSV-B, and bovine-RSV diverged at around 450 years ago, and RSV-A and RSV-B diverged around 250 years ago. Finally, the RSV-A F gene formed eight genotypes (GA1-GA7 and NA1) over the last 80 years. Phylodynamics of RSV-A F gene, including all genotype strains, increased twice in the 1990s and 2010s, while patterns of each RSV-A genotype were different. Phylogenetic distance analysis suggested that the genetic distances of the strains were relatively short (less than 0.05). No positive selection sites were estimated, while many negative selection sites were found. Moreover, the F protein 3D structure mapping and conformational epitope analysis implied that the conformational epitopes did not correspond to the neutralizing antibody binding sites of the F protein. These results suggested that the RSV-A F gene is relatively conserved, and mismatches between conformational epitopes and neutralizing antibody binding sites of the F protein are responsible for the virus reinfection.

Crystal Structure and Immunogenicity of the DS-Cav1-Stabilized Fusion Glycoprotein From Respiratory Syncytial Virus Subtype B

Background:

Respiratory syncytial virus (RSV) subtypes, A and B, co-circulate in annual epidemics and alternate in dominance. We have shown that a subtype A RSV fusion (F) glycoprotein, stabilized in its prefusion conformation by DS-Cav1 mutations, is a promising RSV-vaccine immunogen, capable of boosting RSV-neutralizing titers in healthy adults. In both humans and vaccine-tested animals, neutralizing titers elicited by this subtype A DS-Cav1 immunogen were ~ 2- to 3-fold higher against the homologous subtype A virus than against the heterologous subtype B virus.

Methods:

To understand the molecular basis for this subtype difference, we introduced DS-Cav1 mutations into RSV strain B18537 F, determined the trimeric crystal structure, and carried out immunogenicity studies.

Results:

The B18537 DS-Cav1 F structure at 2-Å resolution afforded a precise delineation of prefusion F characteristics, including those of antigenic site Ø, a key trimer-apex site. Structural comparison with the subtype A prefusion F indicated 11% of surface residues to be different, with an alpha-carbon root-mean-square deviation (RMSD) of 1.2 Å; antigenic site Ø, however, differed in 23% of its surface residues and had an alpha-carbon RMSD of 2.2 Å. Immunization of vaccine-tested animals with DS-Cav1-stabilized B18537 F induced neutralizing responses ~100-fold higher than with postfusion B18537 F. Notably, elicited responses neutralized RSV subtypes A and B at similar levels and were directed towards both conserved equatorial and diverse apical regions.

Conclusion:

We propose that structural differences in apical and equatorial sites–coupled to differently focused immune responses–provide a molecular explanation for observed differences in elicited subtype A and B neutralizing responses.

Molecular characterization of the glycoprotein and fusion protein in human respiratory syncytial virus subgroup A: Emergence of ON-1 genotype in Iran

HRSV is a principle cause of infant hospitalization, childhood wheezing and a common pathogen in the elderly. Limited information exists regarding HRSV genotypes in Iran. In order to better understand HRSV strain diversity, we performed an in-depth evaluation of the genetic variability of the HRSV F protein detected in children under two years of age that, presented with acute respiratory symptoms during 2015-2016 in Tehran. A total of 180 nasopharyngeal swabs were evaluated. The HRSV positive samples were genotyped for G and F gene sequences using RT-PCR and sequencing methods. Phylogenetic analysis was performed using the neighbor-joining and maximum likelihood methods. Genetic and antigenic characteristics of the F gene, nucleotide and amino acids in significant positions and immune system binding regions, as well as the p-distance, positive/negative selection site, linear epitopes and glycosylation sites were investigated in all selected sequences. Among the 83 HRSV positive samples, the Fifty-five cases were successfully sequenced. All of them were classified as subgroup A and belonged to the ON-1 genotype, which possessed 72-nt duplication in the G gene. This study is the first report on the emergence of ON-1 in Iran. ON-1 Iranian sequences clustered in three lineages according to virus fusion (F) gene variations. F gene sequence analysis showed that all genetic changes in the isolates from Iran were base substitutions and no deletion/insertions were identified. The low dN/dS ratio and lack of positively selected sites showed that the fusion genes found in the strains from Iran are not under host selective pressure. Continuing and long-term molecular epidemiological surveys for early detection of circulating and newly emerging genotypes are necessary to gain a better understanding of their epidemic potential.

Persistence and continuous evolution of the human respiratory syncytial virus in northern Taiwan for two decades

The study aimed to characterize the molecular epidemiology, phylogenetic relationship, and population dynamics of the G protein gene in clinical respiratory syncytial virus (RSV) strains isolated from northern Taiwan. We analyzed a total of 160 and 116 G protein gene sequences of RSV-A and RSV-B representative strains, respectively, from 804 clinical viral stocks collected between July 2000 and June 2016. Population dynamic patterns of the RSV G protein gene were analyzed using Bayesian inference through the Markov chain Monte Carlo framework. A phylogenetic analysis revealed that RSV-A from Taiwan could be categorized into GA2, GA5, and GA7 lineages. GA2 of RSV-A could be further divided into NA1, NA2, NA4, and ON1 clades. These RSV-A lineages has been replaced over time, whereas RSV-B strains from Taiwan continually evolved from a single lineage with significant time-dependent waves. Four putative positive selection sites were observed in both RSV-A and RSV-B. The Bayesian skyline plot revealed that the local population dynamics of RSV were associated with lineage displacement events. Both circulating subtypes and population dynamics represented a unique local pattern. Our results affirm the necessity of continuing molecular surveillance of RSV to attain a more comprehensive understanding of epidemics.

Sequence Analysis of the Fusion Protein Gene of Human Respiratory Syncytial Virus Circulating in China from 2003 to 2014

The human respiratory syncytial virus (HRSV) fusion (F) protein is important for HRSV infection, but few studies have examined the genetic diversity of the F gene from Chinese samples. In this study, a total of 330 HRSV F sequences collected from different regions of China between 2003 and 2014 were analyzed to understand their genetic characteristics. In addition, these sequences were compared with 1150 HRSV F sequences in Genbank from 18 other countries. In phylogenetic analysis, Chinese HRSV F sequences sorted into a number of clusters containing sequences from China as well as other countries. F sequences from different genotypes (as determined based on the G gene sequences) within a HRSV subgroup could be found in the same clusters in phylogenetic trees generated based on F gene sequences. Amino acid analysis showed that HRSV F sequences from China and other countries were highly conserved. Of interest, F protein sequences from all Chinese samples were completely conserved at the palivizumab binding site, thus predicting the susceptibility of these strains to this neutralizing antibody. In conclusion, HRSV F sequences from China between 2003 and 2014, similar to those from other countries, were highly conserved.

Molecular epidemiology of respiratory syncytial virus

Respiratory syncytial virus (RSV) is a major cause of viral acute respiratory tract infections in young children. The virus is characterised by distinct seasonality that is dependent upon the latitude and its ability to cause reinfection. Respiratory syncytial virus demonstrates a complex molecular epidemiology pattern as multiple strains and/or genotypes cocirculate during a single epidemic. Previous studies have investigated the relationship between RSV genetic diversity, reinfection, and clinical features. Here, we review the evidence behind this relationship together with the impact that the advancement of whole genome sequencing will have upon our understanding and the need for reconsidering the classification of RSV genotypes.

Sequence variability of the respiratory syncytial virus (RSV) fusion gene among contemporary and historical genotypes of RSV/A and RSV/B

Background

The fusion (F) protein of RSV is the major vaccine target. This protein undergoes a conformational change from pre-fusion to post-fusion. Both conformations share antigenic sites II and IV. Pre-fusion F has unique antigenic sites p27, ø, α2α3β3β4, and MPE8; whereas, post-fusion F has unique antigenic site I. Our objective was to determine the antigenic variability for RSV/A and RSV/B isolates from contemporary and historical genotypes compared to a historical RSV/A strain.

Methods

The F sequences of isolates from GenBank, Houston, and Chile (N = 1,090) were used for this analysis. Sequences were compared pair-wise to a reference sequence, a historical RSV/A Long strain. Variability (calculated as %) was defined as changes at each amino acid (aa) position when compared to the reference sequence. Only aa at antigenic sites with variability ≥5% were reported.

Results

A total of 1,090 sequences (822 RSV/A and 268 RSV/B) were analyzed. When compared to the reference F, those domains with the greatest number of non-synonymous changes included the signal peptide, p27, heptad repeat domain 2, antigenic site ø, and the transmembrane domain. RSV/A subgroup had 7 aa changes in the antigenic sites: site I (N = 1), II (N = 1), p27 (N = 4), α2α3β3β4(AM14) (N = 1), ranging in frequency from 7–91%. In comparison, RSV/B had 19 aa changes in antigenic sites: I (N = 3), II (N = 1), p27 (N = 9), ø (N = 4), α2α3β3β4(AM14) (N = 1), and MPE8 (N = 1), ranging in frequency from 79–100%.

Discussion

Although antigenic sites of RSV F are generally well conserved, differences are observed when comparing the two subgroups to the reference RSV/A Long strain. Further, these discrepancies are accented in the antigenic sites in pre-fusion F of RSV/B isolates, often occurring with a frequency of 100%. This could be of importance if a monovalent F protein from the historical GA1 genotype of RSV/A is used for vaccine development.

Molecular evolution of the fusion protein (F) gene in human respiratory syncytial virus subgroup B

In this study, we examined the molecular evolution of the fusion protein (F) gene in human respiratory syncytial virus subgroup B (HRSV-B). First, we performed time-scale evolution analyses using the Bayesian Markov chain Monte Carlo (MCMC) method. Next, we performed genetic distance, linear B-cell epitope prediction, N-glycosylation, positive/negative selection site, and Bayesian skyline plot analyses. We also constructed a structural model of the F protein and mapped the amino acid substitutions and the predicted B-cell epitopes. The MCMC-constructed phylogenetic tree indicated that the HRSV F gene diverged from the bovine respiratory syncytial virus gene approximately 580years ago and had a relatively low evolutionary rate (7.14×10^-4substitutions/site/year). Furthermore, a common ancestor of HRSV-A and -B diverged approximately 290years ago, while HRSV-B diverged into three clusters for approximately 60years. The genetic similarity of the present strains was very high. Although a maximum of 11 amino acid substitutions were observed in the structural model of the F protein, only one strain possessed an amino acid substitution located within the palivizumab epitope. Four epitopes were predicted, although these did not correspond to the neutralization sites of the F protein including the palivizumab epitope. In addition, five N-glycosylation sites of the present HRSV-B strains were inferred. No positive selection sites were identified; however, many sites were found to be under negative selection. The effective population size of the gene has remained almost constant. On the basis of these results, it can be concluded that the HRSV-B F gene is highly conserved, as is the F protein of HRSV-A. Moreover, our prediction of B-cell epitopes does not show that the palivizumab reaction site may be recognized as an epitope during naturally occurring infections.

Molecular evolution of the fusion protein gene in human respiratory syncytial virus subgroup A

We studied the molecular evolution of the fusion protein (F) gene in the human respiratory syncytial virus subgroup A (HRSV-A). We performed time-scaled phylogenetic analyses using the Bayesian Markov chain Monte Carlo (MCMC) method. We also conducted genetic distance (p-distance), positive/negative selection, and Bayesian skyline plot analyses. Furthermore, we mapped the amino acid substitutions of the protein. The MCMC-constructed tree indicated that the HRSV F gene diverged from the bovine RSV (BRSV) gene approximately 550years ago and had a relatively low substitution rate (7.59×10(-4) substitutions/site/year). Moreover, a common ancestor of HRSV-A and -B diverged approximately 280years ago, which has since formed four distinct clusters. The present HRSV-A strains were assigned six genotypes based on F gene sequences and attachment glycoprotein gene sequences. The present strains exhibited high F gene sequence similarity values and low genetic divergence. No positive selection sites were identified; however, 50 negative selection sites were identified. F protein amino acid substitutions at 17 sites were distributed in the F protein. The effective population size of the gene has remained relatively constant, but the population size of the prevalent genotype (GA2) has increased in the last 10years. These results suggest that the HRSV-AF gene has evolved independently and formed some genotypes.

Bioinformatics and Molecular Analysis of the Evolutionary Relationship between Bovine Rhinitis A Viruses and Foot-And-Mouth Disease Virus

Bovine rhinitis viruses (BRVs) cause mild respiratory disease of cattle. In this study, a near full-length genome sequence of a virus named RS3X (formerly classified as bovine rhinovirus type 1), isolated from infected cattle from the UK in the 1960s, was obtained and analyzed. Compared to other closely related Aphthoviruses, major differences were detected in the leader protease (L^pro), P1, 2B, and 3A proteins. Phylogenetic analysis revealed that RS3X was a member of the species bovine rhinitis A virus (BRAV). Using different codon-based and branch-site selection models for Aphthoviruses, including BRAV RS3X and foot-and-mouth disease virus, we observed no clear evidence for genomic regions undergoing positive selection. However, within each of the BRV species, multiple sites under positive selection were detected. The results also suggest that the probability (determined by Recombination Detection Program) for recombination events between BRVs and other Aphthoviruses, including foot-and-mouth disease virus was not significant. In contrast, within BRVs, the probability of recombination increases. The data reported here provide genetic information to assist in the identification of diagnostic signatures and research tools for BRAV.

Epidemiological changes of respiratory syncytial virus (RSV) infections in Israel

RSV is the leading cause of lower respiratory-tract infections in infants and therefore demands in-depth epidemiological characterization. We investigated here the distribution of RSV types in Israel between the years 2005–2012. Clinical samples were collected from 11,018 patients hospitalized due to respiratory illnesses and were evaluated for the presence of various respiratory viruses, including RSV A and RSV B. Until 2008, each year was characterized by the presence of one dominant type of RSV. However, from 2008, both RSV A and B types were detected at significant levels, particularly among infants aged 0–2 years. Furthermore, significant changes in the RSV A and RSV B subtypes circulating in Israel since 2008 were observed. Finally, we demonstrate that, irrespectively of the changes observed in RSV epidemiology, when the pandemic H1N1pdm09 influenza virus appeared in 2009, RSV infections were delayed and were detected when infection with H1N1pdm09 had declined.