Analysis of influenza epidemic in Fujian Province from 2015 to 2022
The results of PCR identification of original specimens of influenza-like cases from the influenza network surveillance in Fujian Province from 2015 to 2022 are detailed in Table 2 and Fig. 1.Compared to the number of original specimens positive for PCR of influenza-like cases in Fujian Province in 2015–2019 (3523–5551), the number in 2020 (930) is significantly and substantially lower.Compared to 2020, the number in 2021 (2514) is significantly higher, all in the Victoria lineage.The Victoria lineage remained at a high level in January-March 2022 (2006, 99.80%). The proportion of Victoria lineage (40.40%, 46.08%, 100% and 99.80%, respectively) surpassed the other types as the dominant endemic strain in 2016, 2019, 2021 and 2022 (January-March).Taken together, this indicates that the Victoria linage dominated most years in Fujian Province from 2015–2022.
Table 2
Composition ratio of different types of influenza surveillance in Fujian Province from 2015 to 2022
Type | 2015 | 2016 | 2017 | 2018 | 2019 | 2020 | 2021 | 2022 (Jan-Mar) |
B/Victoria | 9.28%(327/3523) | 40.40%(1988/4921) | 9.15%(360/3935) | 8.18%(349/4269) | 46.08%(2558/5551) | 37.42%(348/930) | 100.00%(2514/2514) | 99.80%(2006/2010) |
B/Yamagata | 38.52%(1357/3523) | 9.77%(481/4921) | 26.02%(1024/3935) | 31.32%(1337/4269) | 0.04%(2/5551) | 0.11%(1/930) | 0.00%(0/2514) | 0.00%(0/2010) |
A(H3) | 50.75%(1788/3523) | 27.8%(1368/4921) | 48.26%(1899/3935) | 5.06%(216/4269) | 27.44%(1523/5551) | 37.96(353/930) | 0.00%(0/2514) | 0.2%(4/2010) |
A(H1)pdm09 | 1.45%(51/3523) | 22.03%(1084/4921) | 16.57%(652/3935) | 55.45%(2367/4269) | 26.45%(1468/5551) | 24.52%(228/930) | 0.00%(0/2514) | 0.00%(0/2010) |
Amino Acid Variation Of Ha Gene
The antigenic determinant of HA protein of influenza B virus mainly consists of four regions: 120-loop(116–137), 150-loop (141–150), 160-loop (162–170) and 190-helix (197–205), and the 120-loop of the Victoria lineage includes sites 75 and 77[]. The receptor binding sites (RBS) includes 190-helix (196–205), 240-loop (240–245), 140- loop (136–143), Phe-95, Trp-158, His-194, and Tyr-205[]. Differences in amino acid sites between the 75 Fujian Province strains compared to previous years of northern hemisphere vaccine strains. Only the differences in the antigenic determinant cluster sites and receptor binding sites (RBS) are presented in Table 3. There were also some scattered amino acid differences outside the antigenic determinant cluster.
From 2015 to 2022, there were three circulating strains of HA gene with nucleotide length (1758bp, 1752bp and 1749bp). The prevalent strain 1752bp exhibits a deletion of 2AA at loci 162, 163 (Δ2)[]. The 1749bp prevalent strain exhibited a deletion of 3AA at loci 162, 163, 164 (Δ3).
The high number of loci and strains involved in the difference in strains from 2015–2022 in Fujian Province was 120-loop, followed by 150-loop, 160-loop, 190-helix, and RBS.New variant loci E128K and G133R add in 2019–2020, reverting to E128 and G133 in 2021–2022. New variant loci P241Q add in 2021, reverting to P241 in 2022. New variant loci K136E add in 2019–2022. New variant loci H122N/Q, A127T, E128K, P144L, N150K, S197D, K203R add in 2021–2022. The I117V, N129G/D and S197N variant loci were the same as the vaccine strain. The comprehensive analysis shows that the variation of the HA gene in the Victoria lineage in Fujian Province from 2015 to 2022 is diversified and the variant loci are increasing year by year.
Table 3
Variation of HA gene in influenza B virus Victoria lineage in Fujian Province from 2015 to 2022
| 120-loop | 120-loop RBS | 150- loop RBS | 150-loop | 160-loop | 190-helix RBS | RBS |
Vaccine strain / Year | 117 | 122 | 127 | 128 | 129 | 133 | 136 | 141 | 144 | 146 | 147 | 150 | 162 | 163 | 164 | 169 | 197 | 199 | 202 | 203 | 241 |
B/Brisbane/60/2008 | I | H | A | E | N | G | K | G | P | I | T | N | K | N | D | A | S | T | A | K | P |
B/Colorado/06/2017 | V | . | . | . | G | . | . | . | . | . | . | . | – | – | . | . | N | . | . | . | . |
B/Washington/02/2019 | V | . | . | . | D | R | E | . | . | . | . | . | – | – | – | . | N | . | . | . | . |
2015 (2 strains) | 0 | 0 | 0 | 0 | 0 | 0 | 0 | 0 | 0 | 2(V) | 0 | 0 | 0 | 0 | 0 | 0 | 2(N) | 0 | 0 | 0 | 0 |
2016 (12 strains) | 10(V) | 0 | 0 | 0 | 12(D) | 0 | 0 | 0 | 0 | 0 | 0 | 0 | 0 | 0 | 0 | 0 | 11(N)、1(Y) | 0 | 0 | 0 | 0 |
2017 (12 strains) | 12(V) | 0 | 0 | 0 | 12(D) | 0 | 0 | 0 | 0 | 1(V ) | 0 | 0 | 0 | 0 | 0 | 0 | 12(N) | 1(I) | 0 | 0 | 0 |
2018 (12 strains) | 12(V) | 0 | 0 | 0 | 10(D)、2(G) | 0 | 0 | 0 | 0 | 0 | 1(I) | 0 | 2(-) | 2(-) | 0 | 0 | 12(N) | 0 | 1(T) | 0 | 0 |
2019 (12 strains) | 12(V) | 0 | 0 | 1(K) | 10(D)、1(G) | 11(R) | 11(E) | 0 | 0 | 0 | 0 | 0 | 12(-) | 12(-) | 11(-) | 0 | 12(N) | 0 | 0 | 0 | 0 |
2020 (10 strains) | 10(V) | 0 | 0 | 1(K) | 7(D) | 10(R) | 10(E) | 0 | 0 | 0 | 0 | 0 | 10(-) | 10(-) | 10(-) | 0 | 10(N) | 1(A)、1(I) | 0 | 0 | 0 |
2021 (12 strains) | 11(V) | 2(N)、2(Q) | 2(T) | 0 | 12(D) | 0 | 12(E) | 1(R) | 2(L) | 0 | 0 | 12(K) | 12(-) | 12(-) | 12(-) | 0 | 12(D) | 0 | 0 | 2(R) | 10(Q) |
2022 (3 strains) | 3(V) | 3(Q) | 3(T) | 0 | 3(D) | 0 | 3(E) | 0 | 3(L) | 0 | 0 | 3(K) | 3(-) | 3(-) | 3(-) | 1(T) | 3(D) | 0 | 0 | 3(R) | 0 |
Phylogenetic Tree Analysis
Figure 2 shows that the Clade 1 branch of the evolutionary tree contains three length prevalent strains.The Clade 1 branch HA nucleotide length is mostly 1758bp.The △2 branch has only 4 strains (3 strains in Fujian Province and 1strain in other province; 3 strains in 2018 and 1strain in 2019). The △3 branch contains strains from Fujian Province as well as from other provinces for the years 2019–2022.
In 2015–2020, the Fujian Province and the vaccine strains are located in Clade 1. In 2021–2022, the Fujian Province and other province strains are located in Clade 2.Clade2 branch is different from the branch of vaccine strains B/Washington/02/2019, B/Colorado/06/2017, B/Brisbane/60/2008 (Clade1), which indicates that there is a certain distance of evolutionary differentiation, suggesting that the matching degree of vaccine strains may be reduced.The N150K, G184E/K, and N197D loci began to differ in 2020, and the Clade1 to Clade2 transition suggests a change in prevalence.In addition to the large branches there are still a few scattered small branches in the evolutionary tree, suggesting localised population differences in infection.(See Fig. 2 for details).
Glycosylation Sites
Table 4 shows that the eight glycosylation loci of the HA gene 25, 59, 233, 304, 492, 518, 531 and 563 of the Fujian Province strain are relatively conserved. In contrast, the 145, 166, 171 and 197 glycosylation loci showed differential diversification. For example, there was no glycosylation at locus 197 in 2021–2022, whereas in previous years only individual strains were not glycosylated. 145 and 166 glycosylation loci differed in amino acids of individual strains, and only one strain was glycosylated at locus 171.
Table 4
Glycosylation sites of HA gene of influenza B virus Victoria lineage in Fujian Province from 2015 to 2022
Glycosylation site | Amino acid | Remarks |
25 | NVTG (75 strains) | / |
59 | NCTD (75 strains) | / |
233 | NQTE (75 strains) | / |
304 | NKSK (75 strains) | / |
492 | NQTC (75 strains) | / |
518 | NITA (75 strains) | / |
531 | NHTI (75 strains) | / |
563 | NVSC (75 strains) | / |
145 | NITN (71 strains), NVTN (3 strains), and no glycosylation site (1 strain) | No glycosylation. (B / Fujian-Yanping / 1651 / 2018 C2); NVTN(B/Fujian-Cangshan / 1487 / 2015 C1, B /Fujian-Hanjiang / 31 / 2015 C2, B /Fujian-Xiangcheng / 1228 / 2017 C1). |
166 | NKTA (74 strains), NKTT (1 strain) | NKTT (B / Fujian-Gulou / 1152 / 2022 E1) |
171 | NSST (1 strain), No glycosylation(74 strains) | NSST (B / Fujian-Fengze / 11273 / 2021 C2) |
197 | NETQ (56strains),No glycosylation(19 strains) | 19 strains without glycosylation (1 strain in 2016, 1 in 2017, 2 in 2020 and all 15 strains in 2021–2022). |