Next Article in Journal
Environmental Surveillance through Next-Generation Sequencing to Unveil the Diversity of Human Enteroviruses beyond the Reported Clinical Cases
Previous Article in Journal
Cellular Receptors Involved in KSHV Infection
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Viruses
Volume 13
Issue 1
10.3390/v13010119
viruses-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Respiratory Syncytial Virus G Protein Sequence Variability among Isolates from St. Petersburg, Russia, during the 2013–2014 Epidemic Season

by
Vera Krivitskaya
1,
Kseniya Komissarova
1,*,
Maria Pisareva
1,
Maria Sverlova
1,
Artem Fadeev
1,
Ekaterina Petrova
1,
Veronika Timonina
2,
Anna Sominina
1 and
Daria Danilenko
1
1
Department of Etiology and Epidemiology, Smorodintsev Research Institute of Influenza, 197376 Saint-Petersburg, Russia
2
Children’s City Hospital of St. Olga, 194017 Saint-Petersburg, Russia
*
Author to whom correspondence should be addressed.
Viruses 2021, 13(1), 119; https://doi.org/10.3390/v13010119
Submission received: 16 November 2020 / Revised: 10 January 2021 / Accepted: 12 January 2021 / Published: 17 January 2021
(This article belongs to the Section Animal Viruses)
Browse Figures
Versions Notes

Abstract

:
Human respiratory syncytial virus (RSV) is the most common cause of upper and lower respiratory tract infections in infants and young children. It is actively evolving under environmental and herd immunity influences. This work presents, for the first time, sequence variability analysis of RSV G gene and G protein using St. Petersburg (Russia) isolates. Viruses were isolated in a cell culture from the clinical samples of 61 children hospitalized (January–April 2014) with laboratory-confirmed RSV infection. Real-time RT-PCR data showed that 56 isolates (91.8%) belonged to RSV-A and 5 isolates (8.2%) belonged to RSV-B. The G genes were sequenced for 27 RSV-A isolates and all of them belonged to genotype ON1/GA2. Of these RSV-A, 77.8% belonged to the ON1(1.1) genetic sub-cluster, and 14.8% belonged to the ON1(1.2) sub-cluster. The ON1(1.3) sub-cluster constituted a minor group (3.7%). Many single-amino acid substitutions were identified in the G proteins of St. Petersburg isolates, compared with the Canadian ON1/GA2 reference virus (ON67-1210A). Most of the amino acid replacements were found in immunodominant B- and T-cell antigenic determinants of G protein. These may affect the antigenic characteristics of RSV and influence the host antiviral immune response to currently circulating viruses.

1. Introduction

Human respiratory syncytial virus (RSV) is the most common cause of respiratory tract infections in infants and young children worldwide and a major public health burden [1,2]. In developed countries, the prevalence of RSV infection reaches 42–63% in children under 3 years of age hospitalized with acute lower respiratory tract illness. In children, RSV is implicated in: 50–90% of bronchiolitis; 5–40% of pneumonia; and 10–30% of tracheobronchitis [3,4,5]. During the 2012–2013 epidemic season in St. Petersburg, the frequency of RSV infection was 13.5% (comparable to that of influenza) among hospitalized children in the first year of life [6].
RSV features a negative-sense, unsegmented RNA genome and belongs to the family Pneumoviridae, genus Orthopneumovirus. Its envelope consists of two transmembrane glycoproteins: fusion protein (F) and attachment protein (G) that interact with the receptors on the cell surface. Human RSV exists as a single serotype, but is divided into two antigenic groups (RSV-A and RSV-B), depending on interactions with monoclonal antibodies [7].
The molecular epidemiology of RSV is based mainly on G gene phylogeny. According to the G gene sequence, human RSV-A is divided into 12 genotypes (GA 1-7, SAA1, NA1, NA2, ON1, ON2); RSV-B is divided into 21 genotypes (GB 1-4, BA 1-11, SAB 1-4, and URU 1-2) [8,9,10]. RSV-A and RSV-B co-circulate, in various proportions, during almost all epidemics. In most cases, RSV-A dominates [8,11,12,13].
The pathogenesis of RSV infection is complex and often ambiguous. Depending on the circumstances, RSV can induce a protective, Th1-dependent immune response or a Th2-mediated immunopathology (bronchial hyperreactivity and obstruction, pulmonary eosinophilia) [14]. Numerous studies have established the roles of certain G protein antigenic determinants in both antiviral protection and immunopathology [15,16]. Changes in certain G protein epitopes can significantly affect viral properties as a whole; these may enable evasion of an immune response, thereby resulting in frequent reinfections [17,18,19].
RSV is actively evolving under the influence of environmental and herd immunity factors. Understanding this pathogen’s evolutionary patterns necessitates the monitoring of the genetic and antigenic variability of viruses in human circulation from different geographic locations. Furthermore, the study of RSV evolution is crucial for vaccine and antiviral drug development. The aim of this study was to identify the types of RSV circulating in St. Petersburg during the 2013–2014 epidemic season and to evaluate the variability of the RSV-A G protein.
The present study confirms previous findings that the RSV-A viruses of the ON1 genotype superseded all the other genotypes. The genetic diversity of RSV isolates from Russia was analyzed with a special focus on the non-synonymous mutations in immunodominant B- and T-cell epitope-coding regions of the G gene.

2. Materials and Methods

2.1. Clinical Samples

Nasopharyngeal swabs were obtained from children (from 2 weeks to 6 years with median age 12.9 mon. (±11.4 mon.)) within the framework of a GIHSN study at the Smorodintsev Research Institute of Influenza (RII) from January to April 2014. All research procedures involving humans were performed in strict accordance with the ethical standards of the committees (local ethical committee and national standards), and with the Helsinki Declaration of 1975, as revised in 2008. The study protocol was approved by the local RII ethics committee (registered under No. 0005131, 2008, International IRBS for Russia). The parents of all patients provided written informed consent.

2.2. Isolation of RSV in Cell Culture

Isolation of RSV from PCR-confirmed clinical samples was carried out in MA-104 cell culture (obtained from the Cell Culture Collection of the Institute of Cytology, Russian Academy of Sciences, St. Petersburg, Russia) using alpha-MEM serum-free support medium (BioloT, Saint Petersburg, Russia).

2.3. RT-PCR of the G Protein Gene

RNA was extracted from MA-104 cell culture supernatants (containing primary passages of RSV isolates) using the QIAampViral RNA Mini Kit (QIAGEN, Dusseldorf, Germany). The One-Step RT-PCR kit (QIAGEN, Dusseldorf, Germany) was used to amplify the G gene open reading frame. An in-house system was used for subtyping RSV-A and RSV-B. M13 sequences were added to products’ 5′ ends during amplification of the RSV-A G gene (Table 1).

2.4. Capillary Sequencing

PCR products were then extracted from RT-PCR reaction mixtures using the QIAquick PCR purification Kit (QIAGEN, Germany). Sequencing was performed using an ABI GA 3130 (Applied Biosystems, Foster City, CA, USA) and the BigDye Terminator v3.1 Cycle Sequencing Kit (ThermoFisher Scientific, Waltham, MA, USA).

2.5. Phylogenetic Analysis of the G Gene (Nucleic Acid) Sequences

For sequencing, primary passages of 26 RSV-A strains isolated in cell culture were used. G protein gene sequences were aligned with available RSV-A sequences (from GenBank) using the MUSCLE algorithm included in MEGA 7. A maximum-likelihood phylogenetic tree was constructed using MEGA 7 (TN93+G, bootstrap 1000 replications). G gene sequences of strains representing each known RSV-A genotype were retrieved from GenBank and included in the tree (Table 2). High-throughput phylogenetic analysis of the GenBank RSV-A dataset was performed using RAxML [20]. A tree was visualized using the R ‘ggtree’ package [21]. Genetic distances were calculated using dist.alignment function in the R ‘seqinr’ package [22].

2.6. RSV Glycosylation Site Analysis

Putative N-glycosylation sites were predicted using the *NetNGlyc 1.0 server (threshold ≥ 0.5).
O-glycosylation sites were evaluated using the **NetOGlyc 4.0 server (G-score ≥ 0.5).

2.7. Statistical Analysis

The Pearson association coefficient (rA) was used to estimate the relationship between the presence of complications and RSV antigenic groups in the analysis of the contingency table.

3. Results

RSV were isolated in cell culture from the clinical samples of 61 children hospitalized (January–April 2014) with laboratory-confirmed RSV infection. Patient ages ranged from 2 weeks to 6 years. According to RT-PCR data, 56 isolates (91.8%) belonged to RSV-A, and only 5 (8.2%)—RSV-B. A complicated course of RSV infection was observed in 14 out of 61 individuals (23.0%) with ages from 6 months to 3 years (median age 9.4 month ± 3.1). Of those 14, the following complications were recorded in patients with RSV-A: nine acute bronchitis cases (64.3%); one bronchiolitis case (7.1%); two otitis media cases (14.3%); and one acute pneumonia case (7.1%). With RSV-B, one case of acute bronchitis (7.1%) was noted.
The RSV G gene was sequenced for 26 RSV-A strains. All of these viruses belonged to the ON1 genotype (ON1/GA2, according to the latest classification revision [23]). Currently, ON1/GA2 viruses are ubiquitously dominant among RSV-A circulating globally [10,13,23,24,25,26,27,28,29].
The genetic diversity of Russian RSV isolates corresponds to the observed genetic diversity of RSV sampled in other regions (Supplementary Figure S1).
All of the viruses analyzed were genetically close to the reference (ON1/GA2 strain ON67-1210A) (Figure 1). According to the Tabatabai classification [11], division of the RSV-A ON1/GA2 genotype into sub-clusters is based on amino acid substitutions at certain positions (274, 298, 303, 304, 310), compared to the G protein of the reference ON67-1210A. Viruses that are as close as possible to the reference virus are classified as sub-cluster ON1(1.1). Sub-cluster ON1(1.2) is characterized by mutations L274P, L298P, and Y304H. The ON1(1.3) viruses contain V303A and L310P substitutions.
Most of the isolates (21/27 strains, 77.8%) belonged to the ON1(1.1) sub-cluster. Four strains (14.8%) could be identified as ON1(1.2). ON1(1.3) constituted a minor group: only two viruses (7.4%), RSVA/St_Petersburg/5002/2014 and RSVA/St_Petersburg/13764/2014, had substitutions in all five critical positions (L274P, L298P, V303A, Y304H, L310P) (Figure S2).
Four strains had substitutions at two positions in the transmembrane region: T42I (RSVA/St_Petersburg/479/2014, RSVA/St_Petersburg/7108/2014, RSVA/St_Petersburg/6021/2014) and S44A (RSVA/St_Petersburg/5002/2014). In two strains, there are several substitutions in the central domain: RSVA/St_Petersburg/4734/2014 (V167G, F168V, S177R) and RSVA/St_Petersburg/6959/2014 (A184S). In 11 strains, we identified substitutions in 6 positions in the cytoplasmic domains. Nineteen St. Petersburg strains had point mutations at 12 positions in the first hypervariable region, and all studied strains have aa substituted in the second hypervariable region. There are also six positions with substitutions in the heparin-binding domain (Figure 2; Supplementary Table S1).
Using the program NetOglyc, 33 and 48 O-bond Ser and/or Thr were predicted as potential O-glycosylation sites in I and II hypervariable regions of the ON67-1210A virus with a high likelihood: the Matthews correlation coefficient (MCC) was 0.9040–0.9920. According to NetNGlyc, the following N-glycosylation sites were identified for the ON67-1210A virus: two potential sites with scores > 0.5 in the first hypervariable region (103NLS105, 135NTT137) and one in the second hypervariable region (237NTT239). Compared with the ON67-1210A strain, O-glycosylation sites were highly variable in St. Petersburg viruses, while N-glycosylation sites were quite conserved.
In 13 of 26 isolates (48.1%), one or two substitutions were found, at six positions in the first hypervariable region, which result in the loss of O-glycosylation sites (S100G, S102F, T112I, S124P, S128F, T138I). In only one strain (RSVA/St_Petersburg/8951/2014), a substitution (N135D) leading to the loss of an N-glycosylation site was observed in addition to those causing loss of two O-glycosylation sites (S100G and S102F). In 15 out of 26 isolates (55.6%), a loss of one or two potential O-glycosylation sites was detected at nine positions in the second hypervariable region (T207I, T259K, S275G, T282A, T288A/N, S291P, S315P, S316P/F, S317C/A). Almost all G protein glycosylation pattern changes occurred at B cell antigenic determinants previously identified in animals (aa 90–110, 129–152, 169–207, 236–263, 263–298) [19] and humans (aa 90–168, 169–207, 236–298, 262–276) [17,30].

4. Discussion

An association between the clinical manifestations of infection and genetic traits of currently circulating RSV has been demonstrated [12,13]. However, we did not observe differences in the incidence of complications in patients infected with RSV-A or RSV-B (rA = 0.02, p > 0.5). This may be due to the insufficient size of the RSV-B data.
Despite the relative genetic homogeneity of viruses from St. Petersburg, multiple single-amino acid changes were identified in the G protein, compared with the ON67-1210A strain. It has been shown that changes in RSV G protein antigenic determinants can significantly affect viral properties and help it escape immune response [17,18,19].
In four strains, we find substitutions at two positions in the transmembrane region: T42I and S44A. There is little information available on the variability of the transmembrane domain. However, several changes in this region (S51P, L54P) reduced the affinity of RSV-specific neutralizing antibodies [31].
The central domain of the G protein is involved in interactions with host target cells [18,32], and this may explain its high level of sequence conservation. Only in one strain were substitutions detected in this region: RSVA/St_Petersburg/6959/2014 (A184S). This substitution is of special interest. It is located in the important, regulatory T cell epitope (aa 183–195) that stimulates in vitro production of both Th1-mediated (interferon-ɣ, IL-2) and Th2-mediated (IL-4, IL-5, IL-6, IL-13) cytokines in splenocytes. Depending on experimental conditions, mice pre-immunized with a synthetic peptide corresponding to this epitope develop either Th1-dependent protection or Th2-mediated immunopathology (lung eosinophilia) during subsequent RSV infection [33,34]. In addition, position 184 is part of an RSV-A G protein site (182CWAIC186) that mimics part (aa 202–207) of the fractalkine (the member of the CX3C chemokine family). RSV interacts with the fractalkine receptor CX3CR1 expressed on immunocompetent cells. The synthetic peptide 182–186 acts as a fractalkine antagonist blocking CX3CR1 in RSV-infected mice. This leads to inhibition of the influx of CX3CR1+ inflammatory cells to sites of infection. As a result, a decrease in the RSV-specific protective immune response was observed [35,36].
Surprisingly, the cytoplasmic domains of the investigated viruses had a significant degree of variability. In nine strains (33.3%), five positions showed variation of 1–2 aa residues (T4P/N, K5Q/H, Q7P, R8H, K11R). The sequence variability of this region is also observed in the large GenBank dataset (see Figure 2C) and supported by sequences obtained by using different technologies. Currently, there is no literature data about substitutions in this region, as identified in our isolates. Despite the fact that this region is not exposed on the cell surface, it contains B cell epitopes, mutations in which can affect G protein’s antigenic characteristics. For instance, mutations L33P and L35P led to RSV-A being unable to interact with neutralizing monoclonal antibodies [31].
Few studies are devoted to researching the variations in the first hypervariable region of circulating RSV. In comparison with the ON67-1210A virus, 19 St. Petersburg strains (70.4%) had point mutations at 12 positions in this domain. Altered aa position frequencies were found to be distributed almost evenly at the N- and C-termini of this domain. In 68.4% of mutant strains (13/19), substitutions were identified directly in known B cell epitopes (aa 66–90, 90–110, 129–152) [19] and CD4+ T cell epitopes (aa 104–118) [33]. For instance, the S102F substitution (at the previously evolutionarily conserved position) was first noted as typical for ON1/GA2 viruses circulating in the Netherlands in 2016/2017 [13]. That mutation was detected in five of our viruses, yet was isolated in 2014.
The second hypervariable region is G protein’s most changeable domain. In comparison with the ON67-1210A virus, 100% of the St. Petersburg strains had substitutions in this domain. It was found that 27 positions were variable and asymmetrically distributed. The identified changes were observed two times more often at the C-terminus than at the N-terminus of this region. The same applies to the frequency of strains with substitutions at a different part of the G protein amino acid sequence.
In our isolates, most changes in this domain were found in CD4+ T cell determinants (aa 193–203, 208–222) [37,38] and in B cell immunodominant determinants of both animals (aa 169–207, 236–263, 245–274, 263–298) [19] and humans (aa 229–240, 236–298, 265–273, 283–291) [17,30]. It has been found that some of the substitutions identified in seven of our strains (P206Q, T207I, V225A, L274P, or L310P) can affect the anti-RSV immune response [17,31,39]. It should be noted that the single strain (RSVA/St_Petersburg/11101/2014) that had combined mutations at two positions (V225A, L274P) was isolated from a child, aged 22 days, with severe complicated laboratory-confirmed RSV infection.
In 10 isolates (37.0%), one or two mutations were identified at six positions (K196M, P202L, P206Q, T207I, K216N, V225A/L) in the heparin-binding domain (aa residues 187–231), which can act as a viral receptor. This region is rich in positively charged aa residues and seemed a likely candidate for interaction with negatively charged cell surface glycosaminoglycan-heparan sulfate [40]. Moreover, these positions are part of an epitope (aa 177–220) which can contribute to immunopathological manifestations (lung eosinophilia) in infected mice [41].
Twelve positions, in which substitutions were identified in our strains, have been previously shown to be positively selected in ON1/GA2 viruses [10,12,24,42]. In our case, these are: V225A, L248I, E271K, L274P/V, S275G, H290Q, L298P, V303A, Y304H, L310P, L314I, and S317C/A. This suggests that they are important for the development of antiviral immunity, and their variability may be shaped by host immune pressure.
About 60% of the RSV G protein is glycosylated. The high content of glycans prevents effective antigenic presentation of the G protein to immunocompetent cells and allows RSV to avoid the host immune response [43,44,45].
At some RSV-A G protein positions, the reverse amino acid substitutions have been observed several times during viral evolution [10,31,46]. Such replacements have been termed “flip-flops”. When compared with the ON67-1210A strain, 2 known “flip-flop” substitutions were identified in 6 of 26 (22.2%) St. Petersburg isolates: V225A (V↔A); and/or L274P (L↔P). Those positions (225 and 274) have been reported to be under positive selection in ON1/GA2 viruses; this strongly suggests that their variations reflect changes in the immune status of human populations and are the result of immune pressure [10]. Changes at 225 and 274 have been associated with the decreased affinity of specific neutralizing antibodies [31].
Interestingly, our isolates either did not have substitutions at position 274, or had combined substitutions: L274P and L298P (Figure 2). Thus, these isolates had the haplotype PP at positions 274/298 in contrast with the reference strain (ON67-1210A) with a haplotype LL. For the RSV-A ON1/GA2 genotype, aa residue 298 is part of an insertion (aa residues 284–307) that appeared as a duplication of aa 261–283 [42]. As a result of duplication, aa residue 274 corresponds to aa residue 298. Presumably, substitution at position 298, like at 274, may also be reversible in a “flip-flop” (L↔P) manner. This is supported by RSV (ON1/GA2) G protein structures seen in isolates from 2012 to 2017 in many regions of the world. In addition to viruses with the parent LL haplotype and/or modified PP haplotype at positions 274/298, RSV-A viruses featuring a mixed (LP and/or PL) haplotype have been identified in Kenya, Germany, China, Thailand, Bulgaria, the USA, and the Netherlands [11,13,24,25,26,27,47].

5. Conclusions

The data collected and presented in this study show that, during the 2013–2014 epidemic season in St. Petersburg, mainly RSV-A viruses were circulating. According to RSV-A G gene sequence analysis, isolates are close to the ON67-1210A reference strain and belong to the ON1/GA2 genotype, without formation of a separate clade. However, multiple point mutations, identified in antigenic immunodominant B- and T-cell epitopes of the G protein, may be the result of immune pressure and were able to affect the antigenic characteristics of currently circulating RSV-A. The wide variability of surface RSV glycoprotein, responsible for the induction of protective antiviral reactions, may contribute to the emergence of RSV variants with a potential to escape immune response.

Supplementary Materials

The following are available online at https://www.mdpi.com/1999-4915/13/1/119/s1. Figure S1: Genetic diversity of RSV sampled in different regions of the world (based on 1491 sequences available in the GenBank database). Genetic distance was calculated from the matrix containing squared root of the pairwase distances using dist.alignment function from “seqinr” package [22]. Black line in the box represents the median value, the length of the box represents the interquartile range (IQR), i.e. the difference between the lower (Q1) and the upper (Q3) quartile. Figure S2: G protein amino acid alignment of St. Petersburg RSV-A isolates, 2013–2014 epidemic season. With grey shadow for 2.1 subcluster, orange—3.1 subcluster. All substitutions discussed are filled with grey, 2.1 subcluster defining substitutions are filled with green, extra substitution defining 3.1 subcluster is filled with brown. Table S1: The occurrence of amino acid substitutions with predicted functional significance in RSV-A G protein dataset (CD—central domain, TMD—transmembrane domain, I HVR—I hypervariable region, II HVR—II hypervariable region). Unique substitutions observed in Russian sequences only are highlighted.

Author Contributions

Conceptualization, V.K. and K.K.; methodology, K.K., V.T., E.P., M.S., M.P.; software, K.K., A.F.; investigation, A.S.; data curation, K.K.; writing—original draft preparation, V.K.; writing—review and editing, K.K.; visualization, K.K.; project administration, D.D. All authors have read and agreed to the published version of the manuscript.

Funding

This study was partially funded by the “Foundation for Influenza Epidemiology” within the framework of the international “Global Influenza Hospital-based Surveillance Network” project.

Institutional Review Board Statement

All research procedures involving humans were performed in strict accordance with the ethical standards of the committees (local ethical committee and national standards), and with the Helsinki Declaration of 1975, as revised in 2008. The study protocol was approved by the local RII ethics committee (registered under No. 0005131, 2008, International IRBS for Russia). The parents of all patients provided written informed consent.

Informed Consent Statement

Informed consent was obtained from all subjects involved in the study.

Data Availability Statement

The data presented in this study are openly available in GenBank, reference number MH142221–MH142238 and MK386440–MK386445.

Conflicts of Interest

The authors declare no conflict of interest. Funding sources had no role in: study design; data collection, analysis, or interpretation; or submission/publication decisions.

References

  1. Bont, L.; Checchia, P.A.; Fauroux, B.; Figueras-Aloy, J.; Manzoni, P.; Paes, B.; Simões, E.A.; Carbonell-Estrany, X. Defining the Epidemiology and Burden of Severe Respiratory Syncytial Virus Infection Among Infants and Children in Western Countries. Infect. Dis. Ther. 2016, 5, 271–298. [Google Scholar] [CrossRef] [Green Version]
  2. Manzoni, P.; Figueras-Aloy, J.; Simões, E.A.F.; Checchia, P.A.; Fauroux, B.; Bont, L.; Paes, B.; Carbonell-Estrany, X. Defining the Incidence and Associated Morbidity and Mortality of Severe Respiratory Syncytial Virus Infection Among Children with Chronic Diseases. Infect. Dis. Ther. 2017, 6, 383–411. [Google Scholar] [CrossRef] [Green Version]
  3. Simoes, E.A.; Carbonell-Estrany, X. Impact of severe disease caused by respiratory syncytial virus in children living in developed countries. Pediatr. Infect. Dis. J. 2003, 22, S13–S18, discussion S18–S20. [Google Scholar] [CrossRef] [PubMed]
  4. Nair, H.; Nokes, D.J.; Gessner, B.D.; Dherani, M.; Madhi, S.A.; Singleton, R.J.; O’Brien, K.L.; Roca, A.; Wright, P.F.; Bruce, N.; et al. Global burden of acute lower respiratory infections due to respiratory syncytial virus in young children: A systematic review and meta-analysis. Lancet 2010, 375, 1545–1555. [Google Scholar] [CrossRef] [Green Version]
  5. Karron, R.A. Vaccines, 6th ed.; Plotkin, S.A., Orenstein, W.A., Offit, P.A., Eds.; Elsevier Inc. Saunders: Philadelphia, PA, USA, 2012; pp. 1146–1153. [Google Scholar]
  6. Sominina, A.A.; Pisareva, M.M.; Buzitskaya, Z.V.; Osidak, L.V.; Sukhovetskaya, V.F.; Afanasyeva, O.I.; Voloshchuk, L.V.; Afanas’eva, V.S.; Obraztsova, E.V.; Golovacheva, E.G.; et al. Peculiarities of Etiology of Respiratory Virus Infections in Hospitalized Patients Depending on the Demographic, Socio-Economic Factors and Previous Vaccination. Epidemiol. I Vakcinoprofilaktika 2015, 14, 74–83. (In Russian) [Google Scholar] [CrossRef] [Green Version]
  7. Anderson, L.J.; Hierholzer, J.C.; Tsou, C.; Hendry, R.M.; Fernie, B.F.; Stone, Y.; McIntosh, K. Antigenic Characterization of Respiratory Syncytial Virus Strains with Monoclonal Antibodies. J. Infect. Dis. 1985, 151, 626–633. [Google Scholar] [CrossRef] [PubMed]
  8. Cui, G.; Zhu, R.; Qian, Y.; Deng, J.; Zhao, L.; Sun, Y.; Wang, F. Genetic Variation in Attachment Glycoprotein Genes of Human Respiratory Syncytial Virus Subgroups A and B in Children in Recent Five Consecutive Years. PLoS ONE 2013, 8, e75020. [Google Scholar] [CrossRef]
  9. Nagasawa, K.; Hirano, E.; Kobayashi, M.; Ryo, A.; Oishi, K.; Obuchi, M.; Ishiwada, N.; Noda, M.; Kuroda, M.; Shimojo, N.; et al. Molecular evolution of the hypervariable region of the attachment glycoprotein gene in human respiratory syncytial virus subgroup B genotypes BA9 and BA10. Infect. Genet. Evol. 2015, 36, 217–273. [Google Scholar] [CrossRef]
  10. Duvvuri, V.R.; Granados, A.; Rosenfeld, P.; Bahl, J.; Eshaghi, A.; Gubbay, J.B. Genetic diversity and evolutionary insights of respiratory syncytial virus A ON1genotype: Global and local transmission dynamics. Sci. Rep. 2015, 5, e14268. [Google Scholar] [CrossRef] [Green Version]
  11. Tabatabai, J.; Prifert, C.; Pfeil, J.; Grulich-Henn, J.; Schnitzler, P. Novel Respiratory Syncytial Virus (RSV) Genotype ON1 Predominates in Germany during Winter Season 2012–13. PLoS ONE 2014, 9, e109191. [Google Scholar] [CrossRef] [Green Version]
  12. Esposito, S.; Piralla, A.; Zampiero, A.; Bianchini, S.; Pietro, G.D.; Scala, A.; Pinzani, R.; Fossali, E.; Baldanti, F.; Principi, N. Characteristics and Their Clinical Relevance of Respiratory Syncytial Virus Types and Genotypes Circulating in Northern Italy in Five Consecutive Winter Seasons. PLoS ONE 2015, 10, e0129369. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  13. Vos, L.M.; Oosterheert, J.J.; Kuil, S.D.; Viveen, M.; Bont, L.J.; Hoepelman, A.I.M.; Coenjaerts, F.E.J. High epidemic burden of RSV disease coinciding with genetic alterations causing amino acid substitutions in the RSV G-protein during the 2016/2017 season in The Netherlands. J. Clin. Virol. 2019, 112, 20–26. [Google Scholar] [CrossRef] [PubMed]
  14. Bohmwald, K.; Gálvez, N.M.S.; Canedo-Marroquín, G.; Pizarro-Ortega, M.S.; Andrade-Parra, C.; Gómez-Santander, F.; Kalergis, A.M. Contribution of Cytokines to Tissue Damage During Human Respiratory Syncytial Virus Infection. Front. Immunol. 2019, 10, 452. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  15. Jorquera, P.A.; Anderson, L.; Tripp, R.A. Human Respiratory Syncytial Virus. In Methods and Protocols/Methods in Molecular Biology; Springer Science and Business Media: New York, NY, USA, 2016; Volume 1442, pp. 1–12. [Google Scholar]
  16. Hijano, D.R.; Vu, L.D.; Kauvar, L.M.; Tripp, R.A.; Polack, F.P.; Cormier, S.A. Role of type I interferon (IFN) in the respiratory syncytial virus (RSV) immune response and disease severity. Front. Immunol. 2019, 10, 566. [Google Scholar] [CrossRef] [Green Version]
  17. Cane, P.A. Analysis of linear epitopes recognised by the primary human antibody response to a variable region of the attachment (G) protein of respiratory syncytial virus. J. Med. Virol. 1997, 51, 297–304. [Google Scholar] [CrossRef]
  18. Kim, S.; Joo, D.H.; Lee, J.B.; Shim, B.S.; Cheon, I.S.; Jang, J.E.; Song, H.H.; Kim, K.H.; Song, M.K.; Chang, J. Dual Role of Respiratory Syncytial Virus Glycoprotein Fragment as a Mucosal Immunogen and Chemotactic Adjuvant. PLoS ONE 2012, 7, e32226. [Google Scholar] [CrossRef]
  19. Lee, J.; Klenow, L.; Coyle, E.M.; Golding, H.; Khurana, S. Protective antigenic sites in respiratory syncytial virus G attachment protein outside the central conserved and cysteine noose domains. PLoS Pathog. 2018, 14, e1007262. [Google Scholar] [CrossRef]
  20. Stamatakis, A. RAxML-VI-HPC: Maximum likelihood-based phylogenetic analyses with thousands of taxa and mixed models. Bioinformatics 2006, 22, 2688–2690. [Google Scholar] [CrossRef]
  21. Yu, G. Using ggtree to Visualize Data on Tree-Like Structures. Curr. Protoc. Bioinform. 2020, 69, e96. [Google Scholar] [CrossRef]
  22. Charif, D.; Lobry, J. SeqinR 1.0-2: A contributed package to the R project for statistical computing devoted to biological sequences retrieval and analysis. In Structural Approaches to Sequence Evolution: Molecules, Networks, Populations, Series Biological and Medical Physics, Biomedical Engineering; Bastolla, U., Porto, M., Roman, H., Vendruscolo, M., Eds.; Springer Verlag: New York, NY, USA, 2007; pp. 207–232. ISBN 978-3-540-35305-8. [Google Scholar]
  23. Trento, A.; Ábrego, L.; Rodriguez-Fernandez, R.; González-Sánchez, M.I.; González-Martínez, F.; Delfraro, A.; Pascale, J.M.; Arbiza, J.; Melero, J.A. Conservation of G-Protein Epitopes in Respiratory Syncytial Virus (Group A) Despite Broad Genetic Diversity: Is Antibody Selection Involved in Virus Evolution? J. Virol. 2015, 89, 7776–7785. [Google Scholar] [CrossRef] [Green Version]
  24. Agoti, C.N.; Otieno, J.R.; Gitahi, C.W.; Cane, P.A.; Nokes, D.J. Rapid Spread and Diversification of Respiratory Syncytial Virus Genotype ON1, Kenya. Emerg. Infect. Dis. 2014, 20, 950–959. [Google Scholar] [CrossRef] [PubMed]
  25. Korsun, N.; Angelova, S.; Tzotcheva, I.; Georgieva, I.; Lazova, S.; Parina, S.; Alexiev, I.; Perenovska, P. Prevalence and genetic characterisation of respiratory syncytial viruses circulating in Bulgaria during the 2014/15 and 2015/16 winter seasons. Pathog. Glob. Health 2017, 111, 351–361. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  26. Thongpan, I.; Mauleekoonphairoj, J.; Vichiwattana, P.; Korkong, S.; Wasitthankasem, R.; Vongpun-sawad, S.; Poovorawan, Y. Respiratory syncytial virus genotypes NA1, ON1, and BA9 are prevalent in Thailand, 2012–2015. Peer. J. 2017, 5, e3970. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  27. Lu, B.; Liu, H.; Tabor, D.E.; Tovchigrechko, A.; Qi, Y.; Ruzin, A.; Esser, M.T.; Jin, H. Emergence of new antigenic epitopes in the glycoproteins of human respiratory syncytial virus collected from a US surveillance study, 2015–17. Sci. Rep. 2019, 9, 3898. [Google Scholar]
  28. Ren, L.; Xia, Q.; Xiao, Q.; Zhou, L.; Zang, N.; Long, X.; Xie, X.; Deng, Y.; Wang, L.; Fu, Z.; et al. The genetic variability of glycoproteins among respiratory syncytial virus subtype A in China between 2009 and 2013. Infect. Genet. Evol. 2014, 27, 339–347. [Google Scholar] [CrossRef]
  29. Fall, A.; Dia, N.; Cisseel, H.A.; Kiori, D.E.; Sarr, F.D.; Sy, S.; Goudiaby, D.; Richard, V.; Niang, M.N. Epidemiology and Molecular Characterization of Human Respiratory Syncytial Virus in Senegal after Four Consecutive Years of Surveillance, 2012–2015. PLoS ONE 2016, 11, e0157163. [Google Scholar] [CrossRef] [Green Version]
  30. Fuentes, S.; Coyle, E.M.; Beeler, J.; Golding, H.; Khurana, S. Antigenic Fingerprinting following Primary RSV Infection in Young Children Identifies Novel Antigenic Sites and Reveals Unlinked Evolution of Human Antibody Repertoires to Fusion and Attachment Glycoproteins. PLoS Pathog. 2016, 12, e1005554. [Google Scholar] [CrossRef]
  31. Walsh, E.E.; Hall, C.B.; Schlesinger, J.J.; Brandriss, M.W.; Hildreth, S.; Paradiso, P. Immunization with glycoprotein subunits of respiratory syncytial virus to protect cotton rats against viral infection. J. Gen. Virol. 1989, 70, 2953–2961. [Google Scholar] [CrossRef]
  32. Lee, H.J.; Lee, J.Y.; Park, M.H.; Kim, J.Y.; Chang, J. Monoclonal Antibody against G Glycoprotein Increases Respiratory Syncytial Virus Clearance In Vivo and Prevents Vaccine-Enhanced Diseases. PLoS ONE 2017, 12, e0169139. [Google Scholar] [CrossRef] [Green Version]
  33. Hancock, G.E.; Tebbey, P.W.; Scheuer, C.A.; Pryharski, K.S.; Heers, K.M.; LaPierre, N.A. Immune responses to the nonglycosylated ectodomain of respiratory syncytial virus attachment glycoprotein mediate pulmonary eosinophilia in inbred strains of mice with different MHC haplotypes. J. Med. Virol. 2003, 70, 301–308. [Google Scholar] [CrossRef]
  34. Johnson, T.R.; Varga, S.M.; Braciale, T.J.; Graham, B.S. V14β T Cells Mediate the Vaccine-Enhanced Disease Induced by Immunization with Respiratory Syncytial Virus (RSV) G Glycoprotein but Not with Formalin-Inactivated RSV. J. Virol. 2004, 78, 8753–8760. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  35. Choi, Y.; Mason, C.S.; Jones, L.P.; Crabtree, J.; Jorquera, P.A.; Tripp, R.A. Antibodies to the Central Conserved Region of Respiratory Syncytial Virus (RSV) G Protein Block RSV G Protein CX3C-CX3CR1 Binding and Cross-Neutralize RSV A and B Strains. Viral Immunol. 2012, 25, 193–203. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  36. Ha, B.; Chirkova, T.; Boukhvalova, M.S.; Sun, H.Y.; Walsh, E.E.; Anderson, C.S.; Mariani, T.J.; An-derson, L.J. Mutation of Respiratory Syncytial Virus G Protein’s CX3C Motif Attenuates Infection in Cotton Rats and Primary Human Airway Epithelial Cells. Vaccines 2019, 7, 69. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  37. Sparer, T.E.; Matthews, S.; Hussell, T.; Rae, A.J.; Garcia-Barreno, B.; Malero, J.A.; Openshaw, P. Eliminating a region of respiratory syncytial virus attachment protein allows induction of protective immunity without vaccine-enhanced lung eosinophilia. J. Exp. Med. 1998, 187, 1921–1926. [Google Scholar] [CrossRef] [Green Version]
  38. Tebbey, P.W.; Hagen, M.; Hancock, G.E. Atypical pulmonary eosinophilia is mediated by a specific amino acid sequence of the attachment (G) protein of respiratory syncytial virus. J. Exp. Med. 1998, 188, 1967–1972. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  39. Cane, P.A.; Pringle, C.R. Evolution of subgroup A respiratory syncytial virus: Evidence for progressive accumulation of amino acid changes in the attachment protein. J. Virol. 1995, 69, 2918–2925. [Google Scholar] [CrossRef] [Green Version]
  40. Shields, B.; Mills, J.; Ghildyal, R.; Gooley, P.; Meanger, J. Multiple heparin binding domains of respiratory syncytial virus G mediate binding to mammalian cells. Arch. Virol. 2003, 148, 1987–2003. [Google Scholar] [CrossRef]
  41. Elliott, M.B.; Pryharski, K.S.; Yu, Q.; Boutilier, L.A.; Campeol, N.; Melville, K.; Laughlin, T.S.; Gupta, C.K.; Lerch, R.A.; Randolph, V.B.; et al. Characterization of Recombinant Respiratory Syncytial Viruses with the Region Responsible for Type 2 T-Cell Responses and Pulmonary Eosinophilia Deleted from the Attachment (G) Protein. J. Virol. 2004, 78, 8446–8454. [Google Scholar] [CrossRef] [Green Version]
  42. Eshaghi, A.; Duvvuri, V.R.; Lai, R.; Nadarajah, J.T.; Li, A.; Patel, S.N.; Low, D.E.; Gubbay, J.B. Genetic Variability of Human Respiratory Syncytial Virus A Strains Circulating in Ontario: A Novel Genotype with a 72 Nucleotide G Gene Duplication. PLoS ONE 2012, 7, e32807. [Google Scholar] [CrossRef] [Green Version]
  43. Collins, P.L.; Chanock, R.M.; Murphy, B.R. Fields Virology, 4th ed.; Fields, B.N., Howley, P.M., Eds.; Lippincott Williams and Wilkins: Philadelphia, PA, USA, 2001; Volume 1, pp. 1182–1215. [Google Scholar]
  44. Collins, P.L.; Karron, R.A. Fields Virology, 6th ed.; Fields, B.N., Knipe, D.M., How-Ley, P.M., Eds.; Wolters Kluwer Health; Lippincott Williams, and Wilkins: Philadelphia, PA, USA, 2013; Volume 1, pp. 1086–1123. [Google Scholar]
  45. Melero, J.A.; Mas, V.; McLellan, J.S. Structural, antigenic and immunogenic features of respiratory syncytial virus glycoproteins relevant for vaccine development. Vaccines 2017, 35, 461–468. [Google Scholar] [CrossRef] [Green Version]
  46. Botosso, V.F.; Zanotto, P.M.D.A.; Ueda, M.; Arruda, E.; Gilio, A.E.; Vieira, S.E.; Stewien, K.E.; Peret, T.C.T.; Jamal, L.F.; Pardini, M.I.D.M.C.; et al. Positive Selection Results in Frequent Reversible Amino Acid Replacements in the G Protein Gene of Human Respiratory Syncytial Virus. PLoS Pathog. 2009, 5, e1000254. [Google Scholar] [CrossRef] [PubMed]
  47. Zou, L.; Yi, L.; Wu, J.; Song, Y.; Huang, G.; Zhang, X.; Liang, L.; Ni, H.; Pybus, O.G.; Ke, C.; et al. Evolution and Transmission of Respiratory Syncytial Group A (RSV-A) Viruses in Guangdong, China 2008–2015. Front. Microbiol. 2016, 7, e1263. [Google Scholar] [CrossRef] [PubMed]
Viruses 13 00119 g001 550
Figure 1. Phylogenetic trees by G gene sequence, St. Petersburg RSV-A isolates, 2013–2014 epidemic season. Radial tree (1491 sequences from GenBank) with a focus on geographic distribution. Rectangular tree with reference strains for each genotype colored in blue; original Canadian reference virus (ON67-1210A)—green; progenitors of RSV-A (Long and A2 strains) colored in red.
Figure 1. Phylogenetic trees by G gene sequence, St. Petersburg RSV-A isolates, 2013–2014 epidemic season. Radial tree (1491 sequences from GenBank) with a focus on geographic distribution. Rectangular tree with reference strains for each genotype colored in blue; original Canadian reference virus (ON67-1210A)—green; progenitors of RSV-A (Long and A2 strains) colored in red.
Viruses 13 00119 g001
Viruses 13 00119 g002a 550Viruses 13 00119 g002b 550
Figure 2. Variability of G protein structural regions compared with the reference strain. Comparison uses St. Petersburg RSV-A isolates (n = 27) from the 2013–2014 epidemic season. The rate of strains (A) and position frequencies (B) with aa substitutions identified in specific G protein structural domains. In brackets, in the signature to the X-axis, the numbering and quantity (n) of the amino acid residues of the regions are indicated. Above the diagrams, the number of divergent strains (A) or amino acid positions (B) is indicated. The positions of the amino acid substitutions are shown with red triangles (filled and empty for non-unique and unique mutations, respectively) at the Shannon entropy plot (C) (based on 1491 sequences available in the GenBank database). CD—central domain, TMD—transmembrane domain, I HVR—I hypervariable region, II HVR—II hypervariable region.
Figure 2. Variability of G protein structural regions compared with the reference strain. Comparison uses St. Petersburg RSV-A isolates (n = 27) from the 2013–2014 epidemic season. The rate of strains (A) and position frequencies (B) with aa substitutions identified in specific G protein structural domains. In brackets, in the signature to the X-axis, the numbering and quantity (n) of the amino acid residues of the regions are indicated. Above the diagrams, the number of divergent strains (A) or amino acid positions (B) is indicated. The positions of the amino acid substitutions are shown with red triangles (filled and empty for non-unique and unique mutations, respectively) at the Shannon entropy plot (C) (based on 1491 sequences available in the GenBank database). CD—central domain, TMD—transmembrane domain, I HVR—I hypervariable region, II HVR—II hypervariable region.
Viruses 13 00119 g002aViruses 13 00119 g002b
Table
Table 1. Primers used for amplification of the RSV-A G gene.
Table 1. Primers used for amplification of the RSV-A G gene.
PrimerSequence
M13F-hRSVAB-G-F5′-TGTAAAACGACGGCCAGTGCAAATGCAAACATGTCCAAA-3′
M13R-hRSVA-G-R5′-CAGGAAACAGCTATGCAACYATACGCTTTTTAAATGACTA-3′
Table
Table 2. RSV-A G protein’ gene GenBank sequences.
Table 2. RSV-A G protein’ gene GenBank sequences.
StrainGenBank Accession No.StrainGenBank Accession No.
ON67-1210AJN257693MON-1-89Z33422
ATCC VR-26 (Long)AY911262MON-5-90Z33427
A2KT992094MON-1-90Z33494
RSVA/Homo_sapiens/USA/79E-495-01/1979KP856969MON-4-90Z33426
RSVA/Homo_sapiens/USA/81E-025-01/1981KP856967MON-9-91Z33431
RSVA/Homo_sapiens/USA/84E-004-01/1984KP258733MAD-3-92Z33455
RSVA/Homo_sapiens/USA/86E-007-01/1986KP258723VH220557KM402663
MAD-1-93Z33414HN_6560KT781349
MAD/GM2_14/12KP792373RSV_A_14163_GD-CHN_2014KX009698
MAD/GM2_1/12KP792361RSA_073_GD-CHN_2015KX009700
MAD/GM2_2/12KP792362RSA_029_GD-CHN_2015KX009699
MAD/GM2_5/12KP792365RSA_173_GD-CHN_2015KX009686
RSVA/Homo_sapiens/USA/TH_10603/2013KU950596A/Shiraz/12/2016KU716108
MON-3-88Z33425A/Shiraz/19/2016KU716109
G gene nucleotide sequences for St. Petersburg RSV-A isolates were submitted in GenBank: MH142221–MH142238 and MK386440–MK386445.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Krivitskaya, V.; Komissarova, K.; Pisareva, M.; Sverlova, M.; Fadeev, A.; Petrova, E.; Timonina, V.; Sominina, A.; Danilenko, D. Respiratory Syncytial Virus G Protein Sequence Variability among Isolates from St. Petersburg, Russia, during the 2013–2014 Epidemic Season. Viruses 2021, 13, 119. https://doi.org/10.3390/v13010119

AMA Style

Krivitskaya V, Komissarova K, Pisareva M, Sverlova M, Fadeev A, Petrova E, Timonina V, Sominina A, Danilenko D. Respiratory Syncytial Virus G Protein Sequence Variability among Isolates from St. Petersburg, Russia, during the 2013–2014 Epidemic Season. Viruses. 2021; 13(1):119. https://doi.org/10.3390/v13010119

Chicago/Turabian Style

Krivitskaya, Vera, Kseniya Komissarova, Maria Pisareva, Maria Sverlova, Artem Fadeev, Ekaterina Petrova, Veronika Timonina, Anna Sominina, and Daria Danilenko. 2021. "Respiratory Syncytial Virus G Protein Sequence Variability among Isolates from St. Petersburg, Russia, during the 2013–2014 Epidemic Season" Viruses 13, no. 1: 119. https://doi.org/10.3390/v13010119

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop