Skip to main content
Download PDF

Phylogenetic analysis of HA and NA genes of influenza A viruses in immunosuppressed inpatients in Beijing during the 2018–2020 influenza seasons

Virology Journal volume 20, Article number: 101 (2023) Cite this article

Abstract

Background

Influenza A viruses have undergone rapid evolution with virulent; however, complete and comprehensive data on gene evolution and amino acid variation of HA and NA in immunosuppressed patients was few. In this study, we analysed molecular epidemiology and evolution of influenza A viruses in immunosuppressed population, and immunocompetent population were used as controls.

Methods

Full sequences of HA and NA of A(H1N1)pdm09 and A(H3N2) were acquired through reverse transcription-polymerase chain reaction (RT-PCR). HA and NA genes were sequenced using the Sanger method and phylogenetically analysed using ClustalW 2.10 and MEGA software version 11.0.

Results

During the 2018–2020 influenza seasons, 54 immunosuppressed and 46 immunocompetent inpatients screened positive for influenza A viruses by using the quantitative real-time PCR (qRT-PCR) were enrolled. 27 immunosuppressed and 23 immunocompetent nasal swab or bronchoalveolar lavage fluid samples were randomly selected and sequenced using the Sanger method. A(H1N1)pdm09 were detected in 15 samples and the remaining 35 samples were A(H3N2) positive. By analyzing the HA and NA gene sequences of these virus strains, we found that all A(H1N1)pdm09 viruses shared high similarities to each other and the HA and NA genes of these viruses exclusively belonged to subclade 6B.1A.1. Some NA genes of A(H3N2) viruses were not in the same clade as those of A/Singapore/INFIMH-16-0019/2016 and A/Kansas/14/2017, which may have led to A(H3N2) being the dominant strain in the 2019–2020 influenza season. Both A(H1N1)pdm09 and A(H3N2) viruses showed similar evolutionary lineages patterns of HA and NA between immunosuppressed and immunocompetent patients. Compared with the vaccine strains, there were no statistically significant of HA and NA genes and amino acid sequences of influenza A viruses in immunosuppressed and immunocompetent patients. However, the oseltamivir resistance substitution of NA-H275Y and R292K have been observed in immunosuppressed patients.

Conclusions

A(H1N1)pdm09 and A(H3N2) viruses showed similar evolutionary lineages patterns of HA and NA between immunosuppressed and immunocompetent patients. Both immunocompetent and immunosuppressed patients have some key substitutions, which should be of note monitored, especially those with potential to affect the viral antigen.

Background

Influenza viruses include A, B, C and D, and among these four types, type A is the most infectious, which can even cause a life-threatening pandemic [1]. A (H1N1) pdm09 and H3N2 are currently main circulating strains of influenza A virus. The hemagglutinin (HA) gene has the fastest mutation of the eight genes in influenza A viruses, followed by the neuraminidase (NA) gene [2]. The main mechanism of influenza A virus variation are antigenic drift and genetic reassortments [3]. Previous studies have shown that influenza A(H1N1)pdm09 was derived from several swine influenza viruses and transmitted to humans in the months prior to the outbreak [4, 5]. Active monitoring of influenza virus at the molecular level can help to understand the evolution of influenza virus and select better vaccine strains [6, 7].

The previous literatures and surveillances of the Centers for Disease Control and Prevention (CDC) of influenza at the molecular level were mainly from immunocompetent patients [8, 9]. With the progress of medical technology, immunosuppressed patients who have undergone haemopoietic stem cell transplantation (HSCT) or solid-organ transplantation (SOT), patients on chronic haemodialysis, and patients receiving systemic corticosteroids, immunosuppressants and biological reagents have increased each year [10,11,12]. Immunosuppressed patients accounted for 10.3% of the global multi-center cross-sectional study of 35,348 adult influenza patients published in 2020 [13]. Complete and comprehensive data on influenza gene evolution and amino acid variation of HA and NA in immunosuppressed patients are lacking, and it is unclear whether they differ from patients with immunocompetent patients.

Variations in HA and NA may lead to the emergence of new clinical features, such as the HA-D239 site mutation is associated with severity of illness, while NA-H275, E119, R292, Q136, and I223 substitutions can lead to neuraminidase inhibitors (NAIs) resistance [14,15,16]. Compared with immunocompetent patients, immunosuppressed patients have higher morbidity and mortality, a longer duration of viral shedding, more frequent complications, and more antiviral resistance [10, 17,18,19,20,21]. It is necessary to comprehensively monitor gene evolution and amino acid variation of HA and NA in immunosuppressed patients, and explore whether it has an effect on clinical characteristics and outcome of these patients.

Clinical characteristics and antiviral therapy for immunosuppressed influenza patients were analysed in our prior study [22], and in this study, we further investigate the molecular evolution and amino acid variation of HA and NA of influenza A viruses during the 2018–2020 influenza seasons in these population. In addition, we intended to explore vital amino acid variations of HA and NA whether have an effect on patients’ clinical characteristics and outcomes.

Materials and methods

Patients and sample collection

Peking University People’s Hospital (PKUPH) is a national influenza surveillance sentinel unit, receiving at least 100,000 inpatients from all Beijing districts annually. During the 2018–2020 influenza seasons (November to the following March), 54 immunosuppressed and 46 immunocompetent inpatients screened positive for influenza A viruses by using the quantitative real-time PCR (qRT-PCR) method were enrolled in this study. Original nasal swab or bronchoalveolar lavage fluid (BALF) samples from these patients were collected and immediately placed in virus transport media tubes for analysis. Immunosuppressed patients were defined as the presence at least one risk factor as follows: congenital/genetic immunocompromise, HIV infection, SOT, HSCT, malignancies receiving chemotherapy, aplastic anemia, chronic haemodialysis, chronic steroid use, immunosuppressive agents use and biological drug use [10,11,12]. Full sequences of HA and NA of A(H1N1)pdm09 and A(H3N2) were acquired through reverse transcription-polymerase chain reaction (RT–PCR) method [23].

Data collection

Data on demographic factors (sex, age, cause of immunosuppression, diabetes, use of corticosteroids in the preceding 3 months, neuraminidase inhibitor use before admission), clinical presentation and complications (maximum body temperature, headache, muscle soreness, rhinorrhoea, sore throat, cough, dyspnoea, gastrointestinal symptoms, altered mental status, symptom onset, coinfection with other pathogens, complications, antiviral treatment, admission to the intensive care unit (ICU), mechanical ventilation, time to fever clearance, death) and laboratory test results [white blood cell (WBC), lymphocyte, and C-reactive protein] were collected through medical records.

RNA extraction and reverse transcription

We extracted RNA from samples using the QIAamp Viral RNA Mini Kit (Cat. No.52904, Qiagen, Hilden, Germany) and performed the reverse transcription with a commercial kit (Cat. No.18080051, Invitrogen, Carlsbad, CA, USA), following the manufacturer’s instructions. The complementary DNAs (cDNAs) generated from the reverse transcription were stored at − 20 °C until use.

Gene sequencing of HA and NA

For sequencing the HA and NA genes, high-fidelity thermostable DNA polymerase (Cat. No.11304011, Invitrogen) was used. Specific primers for HA and NA were shown in Table 1. The PCR amplification system included the cDNA template (4 μl), Autoclaved, distilled water (12.1 μl), 10X High Fidelity PCR Buffer (2 μl), 50 mM MgSO4 (0.6 μl), 10 mM dNTP Mix (0.4 μl), 10 µM forward primer (0.4 μl), 10 µM reverse primer (0.4 μl), and Platinum® Taq DNA Polymerase High Fidelity (0.1 μl of 5U/µL). The PCR conditions were: 94 °C for 3 min, followed by 40 cycles of 94 °C for 15 s, 60 °C for 30 s and 72 °C for 2 min, with extension at 72 °C for 10 min. PCR products were analyzed by the method of electrophoresis. 27 immunosuppressed and 23 immunocompetent samples were randomly selected and PCR products were directly sequenced using the Sanger method (BioGerm, Shanghai).

Table 1 PCR primers for HA and NA full genome sequencing of A(H1N1) pdm09 and A(H3N2)
Full size table

Phylogenetic analysis of HA and NA genes

The nucleotide sequences generated were aligned with publicly available sequences of influenza A viruses available in the National Center for Biotechnology Information (NCBI) and Global Initiative of Sharing All Influenza Data (GISAID) database. Full-length reference sequences, which were phylogenetically related to the H1N1 and H3N2 viruses, were included with our sequences for phylogenic analysis using MEGA v.11.0. Sequences with (a) evidence of lab errors and (b) 100% similarity were discarded. All HA and NA gene sequences were aligned using ClustalW 2.10. Phylogenetic analysis of HA and NA gene sequences and substitution analysis of HA and NA proteins were performed using MEGA software version 11.0. Sequences of HA and NA genes in our study have been deposited into NCBI with the accession number OQ455958-455972 and OQ456032-456116.

Statistical analysis

Statistical analysis was performed using SPSS statistical software version 22.0 (SPSS Inc., Chicago, IL, USA). Two-group comparisons of normally distributed data were performed with the independent samples t-test. Frequency comparisons were made with the χ2 test. P values < 0.05 were considered statistically significant.

Results

Epidemics characteristics of influenza virus

According to the weekly influenza data released by Chinese National Influenza Center, the epidemiology of influenza virus in northern China from April 2018 to May 2020 is different in the 2018–2019 and 2019–2020 influenza seasons. As shown in Fig. 1, in the 2018–2019 influenza season, the epidemic of A(H1N1)pdm09 started early, peaked at the end of January, and began to decline in February, while A(H3N2) and B gradually increased, and both A(H3N2) and B showed small peaks at the end of March. In the 2019–2020 influenza season, A(H3N2) was dominant and peaked in early January. In the same period, although A(H1N1)pdm09 and B influenza viruses co-existed, they were significantly less than A(H3N2). In February 2020, control measures were taken against the coronavirus disease outbreak, and the circulation of many respiratory diseases was significantly reduced including influenza [24,25,26,27].

Fig. 1
figure 1

Distribution of influenza viruses in northern China. The data were from weekly data of Influenza Laboratory Surveillance Network during April 2018 to May 2020. Red is A(H1N1)pdm09, green is influenza H3N2, and blue is influenza B

Full size image

Distribution of immunosuppressed factors

During the 2018–2020 influenza seasons, 54 immunosuppressed and 46 immunocompetent inpatients were verified using RT–PCR. Our hospital was not designated hospitals of HIV infection, so patients with HIV infection were not involved in this study. The immunosuppressed factors in this study were as follows: 38 patients with malignancies receiving chemotherapy, 8 patients with chronic steroid use (4 patients also have a history of immunosuppressive or biological agents), 6 patients with HSCT (including 1 patient with aplastic anemia), and 2 patients with chronic haemodialysis. There was no congenital/genetic immunocompromise.

Phylogenetic analysis of HA and NA genes

Twenty-seven inpatients with immunosuppression and 23 immunocompetent inpatients of our study were randomly selected and sequenced using the Sanger method. A(H1N1)pdm09 were detected in 15 samples and the remaining 35 samples were A(H3N2) positive. HA and NA gene sequences were compared with other A(H1N1)pdm09 and A(H3N2) sequences on NCBI and GISAID.

Phylogenetic analyses of the HA genes showed the tested A(H1N1)pdm09 viruses felled together with the vaccine strain A/Brisbane/02/2018 rather than A/Michigan/45/2015 (H1N1), and belonged to subclade 6B.1A.1 (Fig. 2). Like HA genes, the NA genes showed the same evolutionary pattern. These results suggest that these viruses shared common evolutionary lineages to the vaccine strain A/Brisbane/02/2018, and HA and NA of A(H1N1)pdm09 viruses showed similar evolutionary lineages patterns between immunosuppressed and immunocompetent patients.

Fig. 2
figure 2

Phylogenetic tree based on HA and NA nucleotide sequences of A(H1N1)pdm09 from 2009 to 2020. Phylogenetic analysis of HA and NA gene sequences was performed with the Hasegawa-Kishino-Yano and Tamura 3-parameter model respectively which were the best fit for our data using MEGA software version 11.0, with gamma-distributed rates. The reliability of the maximum-likelihood tree was run by bootstrap analysis with 1000 replications. Filled triangle represents the vaccine strain. Red represents viral strains from immunosuppressed patients and blue represents viral strains from immunocompetent patients

Full size image

Phylogenetic analyses of the NA genes showed the tested A(H3N2) viruses felled together with the vaccine strain A/Singapore/INFIMH-16-0019/2016 except A/Beijing/32003/2018 which felled into A/Kansas/14/2017 (Fig. 3). While some NA genes of A(H3N2) viruses were not in the same clade as those of A/Singapore/INFIMH-16–0019/2016 and A/Kansas/14/2017, which may have led to A(H3N2) being the dominant strain in the 2019–2020 influenza season. HA and NA of A(H3N2) viruses also showed similar evolutionary lineages patterns between immunosuppressed and immunocompetent patients.

Fig. 3
figure 3

Phylogenetic tree based on HA and NA nucleotide sequences of A(H3N2) from 2009 to 2020. Phylogenetic analysis of HA and NA gene sequences was performed with the General Time Reversible and Tamura 3-parameter model respectively which were the best fit for our data using MEGA software version 11.0, with gamma-distributed rates. The reliability of the maximum-likelihood tree was run by bootstrap analysis with 1000 replications. Filled triangle represents the vaccine strain. Red represents viral strains from immunosuppressed patients and blue represents viral strains from immunocompetent patients

Full size image

Homology analysis of HA and NA genes

The HA and NA genes of the 15 A(H1N1)pdm09 and 35 A(H3N2) viruses were fully sequenced and analysed of homology. HA genes of these A(H1N1)pdm09 viruses shared 98.00% to 100.00% nucleotide similarity and 97.53% to 100.00% amino acid identity to each other. NA genes of these A(H1N1)pdm09 viruses shared 98.47% to 100.00% nucleotide similarity and 96.81% to 100.00% amino acid identity to each other. As shown in Additional files 1 and 2 Supplementary Tables 1, 2, A(H1N1)pdm09 isolates in this study shared 97.41% to 98.94% nucleotide similarity and 97.23% to 99.12% amino acid identity of their HA and NA genes to A/Michigan/45/2015 (H1N1) vaccine strain, respectively. Moreover, A(H1N1)pdm09 isolates shared 98.24% to 99.36% nucleotide similarity and 97.53% to 99.36% amino acid identity of their HA and NA genes to A/Brisbane/02/2018 vaccine strain, respectively. These results suggest that HA and NA of these viruses contained similar gene constellation and possessed high gene identity to the vaccine strain A/Brisbane/02/2018.

Table 2 Key substitutions of A(H1N1)pdm09 virus in this study
Full size table

HA genes of these A(H3N2) viruses shared 94.59% to 100.00% nucleotide similarity and 93.83% to 100.00% amino acid identity to each other. NA genes of these A(H3N2) viruses shared 95.82% to 100.00% nucleotide similarity and 95.11% to 100.00% amino acid identity to each other. As shown in Additional files 3 and 4 Supplementary Tables 3, 4, A(H3N2) isolates in this study shared 95.67% to 98.88% nucleotide similarity and 94.47% to 98.24% amino acid identity of their HA and NA genes to A/Singapore/INFIMH-16-0019/2016 vaccine strain, respectively. Moreover, A(H3N2) isolates shared 94.71% to 99.76% nucleotide similarity and 93.83% to 100.00% amino acid identity of their HA and NA genes to A/Kansas/14/2017 vaccine strain, respectively. These results suggest that HA and NA of these viruses contained similar gene sequences, but some of which were different from above vaccine strains.

Table 3 Key substitutions of A(H3N2) virus in this study
Full size table

Further comparisons between these two groups, we found that there were no statistically significant of HA and NA genes and amino acid sequences of influenza A viruses in immunosuppressed and immunocompetent patients (Additional files 1, 2, 3 and 4 Supplementary Tables 1, 2, 3, 4).

Substitution analysis of HA and NA proteins

Compared with the vaccine strains of A/California/07/2009, A/Michigan/45/2015 and A/Brisbane/02/2018, HA of A(H1N1)pdm09 virus in this study exhibited 11 key substitutions, and NA of A(H1N1)pdm09 virus exhibited 12 key substitutions, as shown in Table 2. Furthermore, the oseltamivir resistance substitution of NA-H275Y have been observed in one immunosuppressed patients.

Compared with the vaccine strains of A/Hong Kong/4801/2014, A/Singapore/INFIMH-16–0019/2016, and A/Kansas/14/2017, HA of A(H3N2) virus in this study exhibited 10 key substitutions, and NA of A(H3N2) virus also exhibited 10 key substitutions, as shown in Table 3. Two virus strains have been observed oseltamivir resistance substitution NA- R292K in immunosuppressed patients.

Clinical characteristics of influenza A infections with key amino acid variations

Both immunocompetent and immunosuppressed patients have some key amino acid variations. The oseltamivir resistance substitution of NA-H275Y and NA-R292K have been observed in immunosuppressed patients. One A(H1N1)pdm09 virus strain showed V190I mutation at the HA receptor binding site 190 helix. HA-A158E mutation of A(H1N1)pdm09 and I156L/K mutation of A(H3N2) were also be found. V190I and I156L mutation came from immunocompetent patients, while A158E and I156K mutation came from immunosuppressed patients. N222K and G298A near enzyme active center of NA in A(H1N1)pdm09 were also derived from immunosuppressed patients.

In this study, we described clinical characteristics of influenza A infections with vital amino acid variations (Additional file 5: Supplementary Table 5). NA-H275Y of A(H1N1)pdm09 came from a HSCT patients with aplastic anemia, who had oseltamivir and hormone exposure before influenza detection, complicated with coinfections, pneumonia and ARDS, and admitted to intensive care unit (ICU). He was treated with double dose of oseltamivir combined with peramivir on the basis of antibiotics and hormone therapy, and body temperature returned to normal 10 days later. NA-R292K of A(H3N2) came from patients receiving chemotherapy for hematopoietic malignancies and breast cancer respectively. None of these two patients had a history of exposure to NAIs before influenza detection. The former had a peak temperature of 40℃, complicated with pneumonia and acute kidney injury, and was treated with standard dose of peramivir for 2 days and oseltamivir for 5 days, and body temperature returned to normal 5 days later. The latter had a peak temperature of 37.7℃, a history of diabetes and no influenza-related complications. She was treated with standard dose of oseltamivir for 5 days, and his fever cleared after 2 days of treatment. Co-occurring mutations of HA and NA occurred in some patients with different clinical characteristics, for example, patient with co-occurring mutations of HA-A158E, NA-I99V and NA-G298A in A(H1N1)pdm09 did not have influenza-related complications such as pneumonia, and fever cleared 3 days after treatment with standard dose of oseltamivir. However, patient with co-occurring mutations of HA-M274I and NA-S340F complicated with pneumonia and ARDS, and were admitted to ICU. Fever cleared 12 days after treatment with oseltamivir and peramivir. NA-N141S mutant strains of A(H3N2) were accompanied by NA-I257V mutation, and NA-R292K mutant strains were accompanied by NA-E277K mutation.

Discussion

Influenza A virus had a significant impact on public health during the 2018–2020 influenza seasons (Fig. 1). The genetic evolution of influenza virus is gradual, and our team conducted active monitoring of influenza at the molecular level since 2012 [28,29,30,31]. In this study, we analysed molecular epidemiology and evolution of HA and NA of influenza A viruses in immunosuppressed population, and immunocompetent population were used as controls. Phylogenetic analyses of HA and NA genes showed the tested A(H1N1)pdm09 viruses fell together with the vaccine strain A/Brisbane/02/2018 rather than A/Michigan/45/2015 (H1N1) during 2018–2020, and belonged to subclade 6B.1A.1, which was basically consistent with foreign reports [32]. Some NA genes of A(H3N2) viruses were not in the same clade as those of A/Singapore/INFIMH-16–0019/2016 and A/Kansas/14/2017, which may have led to A(H3N2) being the dominant strain in the 2019–2020 influenza season. Both A(H1N1)pdm09 and A(H3N2) viruses showed similar evolutionary lineages patterns of HA and NA between immunosuppressed and immunocompetent patients. Compared with the vaccine strains, there were no statistically significant of HA and NA genes and amino acid sequences of influenza A viruses in immunosuppressed and immunocompetent patients. However, the oseltamivir resistance substitution of NA-H275Y and R292K have been observed in immunosuppressed patients.

The epitopes that affect the antigenicity of influenza are mostly located on the surface of HA protein. Antigenic epitopes of A(H1N1)pdm09 include Sa, Sb, Ca, Cb, Pa and Pb, and antigenic epitopes of A(H3N2) include A, B, C, D and E. Both A(H1N1)pdm09 and A(H3N2) have same receptor-binding domain (RBD) on HA, which is composed of 190 helices, 130 rings and 220 rings [33]. In this study, HA-N146D, N173K and K226M substitutions of A(H1N1)pdm09 were at the Ca epitope, A158E was at the Sa epitope, V190I was at Sb epitope, and M274I was at Pa epitope. HA-G94S substitutions of A(H3N2) were located at E epitopes, Q96H, Q213R, S214P and I230V were located at D epitopes, S140N was located at A epitopes, I156L, I156K and V198I were located at B epitopes, and N312S were located at C epitopes. Above mutations of these epitopes may affect the antigenic characteristics of the virus. It has been reported that sites 132, 133, 135, 189, 190, 192, 193, and 197 are aligned to form an antigenic ridge on RBD, which is directly involved in receptor ligand interaction. Other candidate sites (e.g. 155, 156, 158, and 159) are located on a ring surrounding RBD, which can play an important role in immunization [33]. In this study, one A(H1N1)pdm09 virus strain showed V190I mutation at the HA receptor binding site 190 helix. HA-A158E mutation of A(H1N1)pdm09 and I156L/K mutation of A(H3N2) were also be found. In addition to RBS, vestigial esterase domain (ED) also participated in antigenic drift [33], and the variation of G94S derived from a patient with leukemia were located in the ED domain. The 220-ring mutation is often associated with host specificity. For example, leucine at 226 in humans preferentially recognizes alpha-2,6 sialic acid receptors, and glutamine at 226 in birds preferentially recognizes alpha-2,3 sialic acid receptors [34]. In this study, we got H1-K226M mutation from a patient on long-term hormone therapy.

Activity of NA protein is an important component of influenza virus infection. The epitopes of NA protein consist of 83–143, 156–190, 252–303, 330, 332, 340–345, 368, 370, 387–395, 400, 431–435 and 448–468. The enzyme catalytic sites of NA protein include the central site of enzyme activity (R118, D151, R152, R225, E277, R293, R368, Y402) and the auxiliary site (E119, R156, W179, S180, D/N199, I223, E228, H275, E278, N295, E425) [35]. In this study, NA of A(H1N1)pdm09 virus exhibited 12 significant substitutions, and N222K and G298A are close to the enzyme active center. Among of these substitutions, I99V, G298A, S340F and N341D have not been reported in the past. NA-I99V and G298A appeared in one lung cancer chemotherapy patients with HA-A158E, and S340F appeared in one leukemia patients with HA-M274I. Patients with co-occurring mutations of HA-A158E, NA-I99V and NA-G298A in A(H1N1)pdm09 did not have influenza-related complications; however patients with co-occurring mutations of HA-M274I and NA-S340F complicated with pneumonia and ARDS, and were admitted to ICU. HA and NA have complementary and antagonistic effects. Both HA and SA binding affinity and NA enzyme activity may disrupt the balance between HA and NA, but the co-occurring mutations of HA and NA genes may compensate for such imbalance [34, 36,37,38], which requires further study.

One NA-H275Y mutation and two NA-R292K mutations were found, which were the first time that our team found drug-resistant mutations at H275 and R292 sites. In order for oseltamivir to bind correctly, NA must be rearranged to form a pocket, and the key of rearrangements is E277 rotation so as to bind to R225. Vitro modeling and X-ray crystallography have shown that H275Y inhibits the rotation of E277 residues, thereby preventing the formation of pockets [39]. In this study, it was found that the virus strain with R292K mutation was accompanied by E277K mutation. Whether the resistance mechanism is similar to that of H275Y needs further study. In addition to the influence of NAIs, NA subtypes of influenza viruses may also contribute to mutate. H275Y is dominant in N1 subtypes, such as A (H1N1) and A (H5N1), while R292K are molecular markers of oseltamivir resistance of A(H3N2) virus [40, 41]. Immunosuppressed influenza patients have a long virus shedding time, and long-term exposure to neuraminidase inhibitors can promote the selection of drug-resistant variants [42, 43]. In this study, NA-H275Y of A(H1N1)pdm09 came from a HSCT patients with aplastic anemia, who had oseltamivir exposure before influenza detection, complicated with coinfections, pneumonia and ARDS, and admitted to ICU. NA-R292K of A(H3N2) came from patients receiving chemotherapy for hematopoietic malignancies and breast cancer respectively. None of these two patients had a history of exposure to neuraminase inhibitors before influenza detection. Cases of oseltamivir-resistant A(H1N1)pdm09 initially appeared mainly in immunocompromised patients receiving oseltamivir treatment, later in the 2010–2011 influenza season in the UK and other parts of the world, an increasing number of oseltamivir-resistant patients were found to have no history of oseltamivir use [44]. Peramivir binds to sialic acid residues in a manner similar to oseltamivir and is also affected by H275Y mutations, so combination of two NAIs was not recommended [45]. Double dose application of NAIs was recommended for immunosuppressed patients in earlier studies [19], but current studies have found that patients with double dose of NAIs have good tolerance, but no virologic or clinical advantage were observed [46,47,48]. If symptoms persist and tests consistently positive for the virus, especially in immunosuppressed patients, a longer course of treatment may be considered [23]. Due to some differences in the structure of NA, influenza A(H3N2) is less likely to show NAIs resistance than influenza A(H1N1)pdm09, which generally does not lead to significant loss of NA enzymology function and decreased viral fitness [44]. The clinical outcome of patients with R292K mutation in this study was consistent with this.

Conclusions

This study provided complete and comprehensive data on gene evolution and amino acid variation of HA and NA in immunosuppressed patients with influenza A infections, and described clinical characteristics of key amino acid variations of HA and NA. A(H1N1)pdm09 and A(H3N2) viruses showed similar evolutionary lineages patterns of HA and NA between immunosuppressed and immunocompetent patients. Both immunocompetent and immunosuppressed patients have some key substitutions, which should be of note monitored, especially those with potential to affect the viral antigen. However, the number of this study may not be sufficient to draw the firm conclusions, further studies with a larger sample size will be needed to confirm and extend our findings.

Availability of data and materials

The datasets used and/or analysed in the current study are available from the corresponding author upon reasonable request.

Abbreviations

CDC:

Centers for Disease Control and Prevention

ED:

Esterase domain

GISAID:

Global Initiative of Sharing All Influenza Data

HA:

Hemagglutinin

HSCT:

Haemopoietic stem cell transplantation

ICU:

Intensive care unit

NA:

Neuraminidase

NAIs:

Neuraminidase inhibitors

NCBI:

National Center for Biotechnology Information

PKUPH:

Peking University People’s Hospital

qRT-PCR:

Quantitative real-time PCR

RBD:

Receptor-binding domain

RT–PCR:

Reverse transcriptase–polymerase chain reaction

SOT:

Solid-organ transplantation

μl:

Microlitre

WBC:

White blood cell

References

  1. Rambaut A, Pybus OG, Nelson MI, Viboud C, Taubenberger JK, Holmes EC. The genomic and epidemiological dynamics of human influenza A virus. Nature. 2008;453(7195):615–9.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  2. Tewawong N, Prachayangprecha S, Vichiwattana P, Korkong S, Klinfueng S, Vongpunsawad S, et al. Assessing antigenic drift of seasonal Influenza A(H3N2) and A(H1N1)pdm09 viruses. PLoS ONE. 2015;10(10): e0139958.

    Article  PubMed  PubMed Central  Google Scholar 

  3. Al Khatib HA, Al Thani AA, Gallouzi I, Yassine HM. Epidemiological and genetic characterization of pH1N1 and H3N2 influenza viruses circulated in MENA region during 2009–2017. BMC Infect Dis. 2019;19(1):314.

    Article  PubMed  PubMed Central  Google Scholar 

  4. Garten RJ, Davis CT, Russell CA, Shu B, Lindstrom S, Balish A, et al. Antigenic and genetic characteristics of swine-origin 2009 A(H1N1) influenza viruses circulating in humans. Science. 2009;325(5937):197–201.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  5. Mena J, Ariyama N, Navarro C, Quezada M, Brevis C, Rojas D, et al. Ubiquitous influenza A virus in Chilean swine before the H1N1pdm09 introduction. Transbound Emerg Dis. 2021;68(6):3174–9.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  6. Graham M, Liang B, Van Domselaar G, Bastien N, Beaudoin C, Tyler S, et al. Nationwide molecular surveillance of pandemic H1N1 influenza A virus genomes: Canada, 2009. PLoS ONE. 2011;6(1): e16087.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  7. Smith GJ, Vijaykrishna D, Bahl J, Lycett SJ, Worobey M, Pybus OG, et al. Origins and evolutionary genomics of the 2009 swine-origin H1N1 influenza A epidemic. Nature. 2009;459(7250):1122–5.

    Article  CAS  PubMed  Google Scholar 

  8. Kuzmanovska M, Boshevska G, Janchevska E, Buzharova T, Simova M, Peshnacka A, et al. A comprehensive molecular and epidemiological characterization of influenza viruses circulating 2016–2020 in North Macedonia. Front Microbiol. 2021;12: 713408.

    Article  PubMed  PubMed Central  Google Scholar 

  9. Wedde M, Biere B, Wolff T, Schweiger B. Evolution of the hemagglutinin expressed by human influenza A(H1N1)pdm09 and A(H3N2) viruses circulating between 2008–2009 and 2013–2014 in Germany. Int J Med Microbiol. 2015;305(7):762–75.

    Article  CAS  PubMed  Google Scholar 

  10. Kunisaki KM, Janoff EN. Influenza in immunosuppressed populations: a review of infection frequency, morbidity, mortality, and vaccine responses. Lancet Infect Dis. 2009;9(8):493–504.

    Article  PubMed  PubMed Central  Google Scholar 

  11. Di Pasquale MF, Sotgiu G, Gramegna A, Radovanovic D, Terraneo S, Reyes LF, et al. Prevalence and etiology of community-acquired pneumonia in immunocompromised patients. Clin Infect Dis. 2019;68(9):1482–93.

    Article  PubMed  Google Scholar 

  12. Ramirez JA, Musher DM, Evans SE, Dela Cruz C, Crothers KA, Hage CA, et al. Treatment of community-acquired pneumonia in immunocompromised adults: a consensus statement regarding initial strategies. Chest. 2020;158(5):1896–911.

    Article  PubMed  Google Scholar 

  13. Collins JP, Campbell AP, Openo K, Farley MM, Cummings CN, Hill M, et al. Outcomes of immunocompromised adults hospitalized with laboratory-confirmed influenza in the United States, 2011–2015. Clin Infect Dis. 2020;70(10):2121–30.

    Article  PubMed  Google Scholar 

  14. Pascalis H, Temmam S, Wilkinson DA, Dsouli N, Turpin M, de Lamballerie X, et al. Molecular evolutionary analysis of pH1N1 2009 influenza virus in Reunion Island, South West Indian Ocean region: a cohort study. PLoS ONE. 2012;7(8): e43742.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  15. Tramontana AR, George B, Hurt AC, Doyle JS, Langan K, Reid AB, et al. Oseltamivir resistance in adult oncology and hematology patients infected with pandemic (H1N1) 2009 virus, Australia. Emerg Infect Dis. 2010;16(7):1068–75.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  16. L’Huillier AG, Abed Y, Petty TJ, Cordey S, Thomas Y, Bouhy X, et al. E119D neuraminidase mutation conferring pan-resistance to neuraminidase inhibitors in an A(H1N1)pdm09 isolate from a stem-cell transplant recipient. J Infect Dis. 2015;212(11):1726–34.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  17. Memoli MJ, Athota R, Reed S, Czajkowski L, Bristol T, Proudfoot K, et al. The natural history of influenza infection in the severely immunocompromised vs nonimmunocompromised hosts. Clin Infect Dis. 2014;58(2):214–24.

    Article  PubMed  Google Scholar 

  18. Ison MG. Influenza prevention and treatment in transplant recipients and immunocompromised hosts. Influenza Other Respir Viruses. 2013;7(Suppl 3):60–6.

    Article  PubMed  PubMed Central  Google Scholar 

  19. Eshaghi A, Shalhoub S, Rosenfeld P, Li A, Higgins RR, Stogios PJ, et al. Multiple influenza A (H3N2) mutations conferring resistance to neuraminidase inhibitors in a bone marrow transplant recipient. Antimicrob Agents Chemother. 2014;58(12):7188–97.

    Article  PubMed  PubMed Central  Google Scholar 

  20. Weinstock DM, Gubareva LV, Zuccotti G. Prolonged shedding of multidrug-resistant influenza A virus in an immunocompromised patient. N Engl J Med. 2003;348(9):867–8.

    Article  PubMed  Google Scholar 

  21. Nichols WG, Guthrie KA, Corey L, Boeckh M. Influenza infections after hematopoietic stem cell transplantation: risk factors, mortality, and the effect of antiviral therapy. Clin Infect Dis. 2004;39(9):1300–6.

    Article  PubMed  Google Scholar 

  22. Liu Y, Wang Y, Mai H, Chen Y, Liu B, Liu Y, et al. Clinical characteristics, risk factors and antiviral treatments of influenza in immunosuppressed inpatients in Beijing during the 2015–2020 influenza seasons. Virol J. 2022;19(1):11.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  23. Uyeki TM, Bernstein HH, Bradley JS, Englund JA, File TM, Fry AM, et al. Clinical practice guidelines by the Infectious Diseases Society of America: 2018 update on diagnosis, treatment, chemoprophylaxis, and institutional outbreak management of seasonal influenzaa. Clin Infect Dis. 2019;68(6):895–902.

    Article  PubMed  PubMed Central  Google Scholar 

  24. Siegers JY, Dhanasekaran V, Xie R, Deng YM, Patel S, Ieng V, et al. Genetic and antigenic characterization of an Influenza A(H3N2) Outbreak in Cambodia and the Greater Mekong Subregion during the COVID-19 Pandemic, 2020. J Virol. 2021;95(24): e0126721.

    Article  PubMed  Google Scholar 

  25. Liu M, Deng L, Wang D, Jiang T. Influenza activity during the outbreak of coronavirus disease 2019 in Chinese mainland. Biosaf Health. 2020;2(4):206–9.

    Article  PubMed  PubMed Central  Google Scholar 

  26. Fairweather SM, Chang CL, Mansell CJ, Shafuddin E, Hancox RJ. Impact of COVID-19 pandemic restrictions on the cardio-respiratory health of New Zealanders. Respirology. 2021;26(11):1041–8.

    Article  PubMed  PubMed Central  Google Scholar 

  27. Kim MC, Kweon OJ, Lim YK, Choi SH, Chung JW, Lee MK. Impact of social distancing on the spread of common respiratory viruses during the coronavirus disease outbreak. PLoS ONE. 2021;16(6): e0252963.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  28. Fang Q, Gao Y, Chen M, Guo X, Yang X, Yang X, et al. Molecular epidemiology and evolution of A(H1N1)pdm09 and H3N2 virus during winter 2012–2013 in Beijing, China. Infect Genet Evol. 2014;26:228–40.

    Article  CAS  PubMed  Google Scholar 

  29. Fang Q, Gao Y, Chen M, Guo X, Yang X, Wei L, et al. Molecular epidemiology and evolution of influenza A and B viruses during winter 2013–2014 in Beijing, China. Arch Virol. 2015;160(4):1083–95.

    Article  CAS  PubMed  Google Scholar 

  30. Liu Y, Wang Y, Liu B, Cong X, Ji Y, Guo X, et al. Phylogenetic analysis and clinical characteristics of the co-occurring mutations in HA and NA genes of influenza A(H1N1)pdm09 viruses during 2015–2017 in Beijing, China. Virol J. 2020;17(1):182.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  31. Liu B, Wang Y, Liu Y, Chen Y, Liu Y, Cong X, et al. Molecular evolution and characterization of hemagglutinin and neuraminidase of influenza A(H1N1)pdm09 viruses isolated in Beijing, China, during the 2017–2018 and 2018–2019 influenza seasons. Adv Virol. 2021;166(1):179–89.

    CAS  Google Scholar 

  32. Miron VD, Bănică L, Săndulescu O, Paraschiv S, Surleac M, Florea D, et al. Clinical and molecular epidemiology of influenza viruses from Romanian patients hospitalized during the 2019/20 season. PLoS ONE. 2021;16(11): e0258798.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  33. Quan L, Ji C, Ding X, Peng Y, Liu M, Sun J, et al. Cluster-transition determining sites underlying the antigenic evolution of seasonal influenza viruses. Mol Biol Evol. 2019;36(6):1172–86.

    Article  CAS  PubMed  Google Scholar 

  34. Byrd-Leotis L, Cummings RD, Steinhauer DA. The Interplay between the host receptor and influenza virus hemagglutinin and neuraminidase. Int J Mol Sci. 2017;18(7):1541.

    Article  PubMed  PubMed Central  Google Scholar 

  35. Air GM. Influenza neuraminidase. Influenza Other Respir Viruses. 2012;6(4):245–56.

    Article  CAS  PubMed  Google Scholar 

  36. Wagner R, Matrosovich M, Klenk HD. Functional balance between haemagglutinin and neuraminidase in influenza virus infections. Rev Med Virol. 2002;12(3):159–66.

    Article  PubMed  Google Scholar 

  37. Bloom JD, Gong LI, Baltimore D. Permissive secondary mutations enable the evolution of influenza oseltamivir resistance. Science. 2010;328(5983):1272–5.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  38. Mitnaul LJ, Matrosovich MN, Castrucci MR, Tuzikov AB, Bovin NV, Kobasa D, et al. Balanced hemagglutinin and neuraminidase activities are critical for efficient replication of influenza A virus. J Virol. 2000;74(13):6015–20.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  39. Davidson S. Treating influenza infection, from now and into the future. Front Immunol. 2018;9:1946.

    Article  PubMed  PubMed Central  Google Scholar 

  40. Zürcher T, Yates PJ, Daly J, Sahasrabudhe A, Walters M, Dash L, et al. Mutations conferring zanamivir resistance in human influenza virus N2 neuraminidases compromise virus fitness and are not stably maintained in vitro. J Antimicrob Chemother. 2006;58(4):723–32.

    Article  PubMed  Google Scholar 

  41. Baz M, Abed Y, Simon P, Hamelin ME, Boivin G. Effect of the neuraminidase mutation H274Y conferring resistance to oseltamivir on the replicative capacity and virulence of old and recent human influenza A(H1N1) viruses. J Infect Dis. 2010;201(5):740–5.

    Article  CAS  PubMed  Google Scholar 

  42. Gooskens J, Jonges M, Claas EC, Meijer A, van den Broek PJ, Kroes AM. Morbidity and mortality associated with nosocomial transmission of oseltamivir-resistant influenza A(H1N1) virus. JAMA. 2009;301(10):1042–6.

    Article  CAS  PubMed  Google Scholar 

  43. Renaud C, Campbell AP. Changing epidemiology of respiratory viral infections in hematopoietic cell transplant recipients and solid organ transplant recipients. Curr Opin Infect Dis. 2011;24(4):333–43.

    Article  PubMed  PubMed Central  Google Scholar 

  44. Lampejo T. Influenza and antiviral resistance: an overview. Eur J Clin Microbiol Infect Dis. 2020;39(7):1201–8.

    Article  PubMed  PubMed Central  Google Scholar 

  45. de Jong MD, Ison MG, Monto AS, Metev H, Clark C, O’Neil B, et al. Evaluation of intravenous peramivir for treatment of influenza in hospitalized patients. Clin Infect Dis. 2014;59(12):e172–85.

    Article  PubMed  Google Scholar 

  46. Lee N, Hui DS, Zuo Z, Ngai KL, Lui GC, Wo SK, et al. A prospective intervention study on higher-dose oseltamivir treatment in adults hospitalized with influenza a and B infections. Clin Infect Dis. 2013;57(11):1511–9.

    Article  CAS  PubMed  Google Scholar 

  47. South East Asia Infectious Disease Clinical Research Network. Effect of double dose oseltamivir on clinical and virological outcomes in children and adults admitted to hospital with severe influenza: double blind randomised controlled trial. BMJ. 2013;346: f3039.

    Article  PubMed Central  Google Scholar 

  48. Torres A, Loeches IM, Sligl W, Lee N. Severe flu management: a point of view. Intensive Care Med. 2020;46(2):153–62.

    Article  PubMed  PubMed Central  Google Scholar 

Download references

Acknowledgements

We thank all the participants for their contributions to this study.

Funding

This work was supported by National key clinical specialty construction projects, Beijing key clinical specialty construction projects, and the research and development fund of Peking University People’s Hospital (grant number RD 2016-14). The funders had no role in the study design, data collection and analysis, decision to publish or preparation of the manuscript.

Author information

Author notes

  1. Yafen Liu and Yue Wang have contributed equally to this work

Authors and Affiliations

  1. Department of Infectious Diseases, Peking University Hepatology Institute, Peking University People’s Hospital, No. 11, Xizhimen South Street, Xicheng District, Beijing, 100044, People’s Republic of China

    Yafen Liu, Yue Wang, Yanxin Wang, Huan Mai, YuanYuan Chen, Yifan Zhang & Yan Gao

  2. Peking University Hepatology Institute, Peking University People’s Hospital, No. 11, Xizhimen South Street, Xicheng District, Beijing, 100044, People’s Republic of China

    Ying Ji & Xu Cong

Authors
  1. Yafen Liu
    View author publications

    You can also search for this author in PubMed Google Scholar

  2. Yue Wang
    View author publications

    You can also search for this author in PubMed Google Scholar

  3. Yanxin Wang
    View author publications

    You can also search for this author in PubMed Google Scholar

  4. Huan Mai
    View author publications

    You can also search for this author in PubMed Google Scholar

  5. YuanYuan Chen
    View author publications

    You can also search for this author in PubMed Google Scholar

  6. Yifan Zhang
    View author publications

    You can also search for this author in PubMed Google Scholar

  7. Ying Ji
    View author publications

    You can also search for this author in PubMed Google Scholar

  8. Xu Cong
    View author publications

    You can also search for this author in PubMed Google Scholar

  9. Yan Gao
    View author publications

    You can also search for this author in PubMed Google Scholar

Contributions

YFL and YW participated in the experiments and analysed the data. YG conceived and designed the study and helped to modify the manuscript. YXW, HM, and XC helped to perform the experiments. YYC and YFZ helped to collect and analyse the data. YJ collected, transported and saved the samples. The first draft of the manuscript was written by YFL, and all authors commented on previous versions of the manuscript. All authors read and approved the final manuscript.

Corresponding author

Correspondence to Yan Gao.

Ethics declarations

Ethics approval and consent to participate

The study protocol was in accordance with the guidelines of the Declaration of Helsinki and was approved by the ethics committees of the Peking University People’s Hospital (PKUPH, IRB No. 2016PHB100-01). We explained the details of our study to each participant, and written informed consent was obtained from all participants prior to their inclusion in the study. Nasal samples and medical data were collected and analysed anonymously.

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Additional information

Publisher's Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary Information

Additional file 1: Table 1

. Nucleotide similarity of HA and NA genes of Apdm09 compared to vaccine strains.

Additional file 2: Table 2

. Amino acid similarity of HA and NA genes of Apdm09 compared to vaccine strains.

Additional file 3: Table 3

. Nucleotide similarity of HA and NA genes of Acompared to vaccine strains.

Additional file 4: Table 4

. Amino acid similarity of HA and NA genes of Acompared to vaccine strains

Additional file 5: Table 5

. Clinical characteristics of influenza A infections with key amino acid variations.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Liu, Y., Wang, Y., Wang, Y. et al. Phylogenetic analysis of HA and NA genes of influenza A viruses in immunosuppressed inpatients in Beijing during the 2018–2020 influenza seasons. Virol J 20, 101 (2023). https://doi.org/10.1186/s12985-023-02067-2

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12985-023-02067-2

Keywords

Virology Journal

ISSN: 1743-422X