# Genetic and antigenic characterization of influenza A(H3N2) in Cameroon during the 2014-2016 influenza seasons

Gwladys C Monamele, Marie-Astrid Vernet, Mohammed R Njankouo, Kathleen Victoir, Jane Francis Akoachere, Damian Anong, Richard Njouom

► To cite this version:

Gwladys C Monamele, Marie-Astrid Vernet, Mohammed R Njankouo, Kathleen Victoir, Jane Francis Akoachere, et al.. Genetic and antigenic characterization of influenza A(H3N2) in Cameroon during the 2014-2016 influenza seasons. PLoS ONE, 2017, 12 (9), pp.e0184411. 10.1371/journal.pone.0184411.pone-01659289

HAL Id: pasteur-01659289

https://pasteur.hal.science/pasteur-01659289

Submitted on 8 Dec 2017

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
RESEARCH ARTICLE

# Genetic and antigenic characterization of influenza A(H3N2) in Cameroon during the 2014-2016 influenza seasons

Gwladys C. Monamele¹,², Marie-Astrid Vernet¹, Mohammed R. Njankouo¹, Kathleen Victoir³, Jane Francis Akoachere², Damian Anong², Richard Njouom¹*

1 National Influenza Centre, Centre Pasteur du Cameroun, Yaoundé, Cameroon, 2 Department of Microbiology and Parasitology, University of Buea, Buea, Cameroon, 3 International Director, Institut Pasteur de Paris, Paris, France

* njouom@pasteur-yaounde.org

<figure><img src="image_2.png" /></figure>

## Abstract

The first outbreak of influenza A(H3N2) occurred in 1968 and caused the third flu pandemic of the 20th century. It has affected multiple countries over time. The best strategy to reduce the burden of influenza is through vaccination whose efficacy varies with respect to the circulating strains. This study was performed to better understand the molecular evolution of influenza A(H3N2) and assess vaccine efficacy in Cameroon. Complete sequences of three gene segments were obtained from 2014 to 2016 influenza seasons in Cameroon. Hemagglutinin (HA), Neuraminidase (NA) and matrix (M) genes of 35 A(H3N2) virus strains were amplified and sequenced. Predicted vaccine efficacy was measured using the P<sub>epitope</sub> model. Phylogenetic analysis of the HA gene showed that all Cameroonian strains had evolved away from the 3C.1-A/Texas/50/2012-like clade. Globally, 2014 virus strains clustered with the 2015–2016 vaccine strain, 3C.3a-A/Switzerland/9715293/2013, whereas 2015 and 2016 virus strains clustered with the 2016–2017 vaccine strain, 3C.2a-A/HongKong/4801/2014. In order to determine the genotypic drug susceptibility to neuraminidase inhibitors and amantadine, the NA and M2 protein coding sequences were analyzed. There was no strain with characteristic mutation for resistance to neuraminidase inhibitors, per contra; all strains possessed the substitution S31N, peculiar of resistance to adamantanes. There was drift in influenza A(H3N2) dominant epitopes B (2014 and 2015) to epitopes A (2016) with a theoretical efficiency in vaccine ranging from low to moderate. The presence of several antigenic site mutations among H3N2 virus strains between 2014–2016 influenza seasons in Cameroon confirms the progressing evolution of circulating H3N2 strains.

<figure><img src="image_3.png" /></figure>

### OPEN ACCESS

**Citation:** Monamele GC, Vernet M-A, Njankouo MR, Victoir K, Akoachere JF, Anong D, et al. (2017) Genetic and antigenic characterization of influenza A(H3N2) in Cameroon during the 2014-2016 influenza seasons. PLoS ONE 12(9): e0184411. https://doi.org/10.1371/journal.pone.0184411

**Editor:** Juan C. de la Torre, The Scripps Research Institute, UNITED STATES

**Received:** May 16, 2017

**Accepted:** August 23, 2017

**Published:** September 6, 2017

**Copyright:** © 2017 Monamele et al. This is an open access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.

**Data Availability Statement:** All relevant data are within the paper and its Supporting Information files.

**Funding:** This work was supported by the U.S. Department of Health and Human Services (DHHS) grant number 6 DESP060001-01-01 via the International Network of Pasteur Institutes.

**Competing interests:** The authors have declared that no competing interests exist.

## Introduction

Influenza A virus is a major cause of acute respiratory disease in humans and is responsible for approximately 250,000–500,000 deaths annually worldwide [1]. Pandemic influenza A virus infection resulted in significant morbidity and mortality in 1918 (H1N1), 1957 (H2N2), 1968 (H3N2), and 2009 (H1N1) [2, 3]. The first outbreak of influenza A/H3N2 was reported in
1968 and has since affected multiple countries[2]. In several seasons, A(H3N2) strains are more prevalent than the other co-circulating viral sub-types [4, 5] with high morbidity and mortality rates [6, 7].

Influenza A nomenclature is based on the genotypic properties of two surface glycoproteins, hemagglutinin (HA) and neuraminidase (NA). Currently, 18 HA and 11 NA subtypes have been recognized [8]. The primary target of host neutralizing antibodies is the HA glycoprotein where it inhibits binding to the sialic acid receptors present in the respiratory tract [9]. There are five identified antigenic sites for A(H3N2), namely, A, B, C, D and E [10]. Non-recognition of some strains by neutralizing antibodies present in the host is the consequence of mutations at these sites [11]. Continuous build-up of mutations at these antigenic sites is the basis for the evolutionary dynamics observed in the influenza virus resulting in antigenic drift [12]. About seven A(H3N2) clades (designated clades 1 to 7) and many sub-clades have evolved over the past few years, providing a significant challenge for the annual design of effective influenza vaccines [13]. At present, viruses in clade 3 are the dominant group, which has diversified further into three major subclades designated 3A, 3B, and 3C, with subclade 3C containing several distinguishable genetic groups [13, 14].

The best strategy to reduce the burden of influenza is through vaccination whose efficacy varies with respect to the circulating strains[15]. Although influenza vaccination is very limited in most of Africa, particularly Sub-Saharan Africa, there have been previous reports of introduction of viruses that have undergone drift due to antiviral and vaccine pressure [16, 17] from other regions. There is therefore, a crucial need for monitoring the genetic and antigenic characteristics as well as drug susceptibility of influenza viruses as a component of pandemic preparedness [16]. To better understand the molecular evolution of influenza and assess vaccine efficacy in Cameroon, three gene segments of A(H3N2) from 2014 to 2016 were characterized.

# Materials and methods

## Sample collection and preparation

Respiratory samples were collected between January 2014 and June 2016 from outpatients with clinical evidence of influenza-like illness (ILI), defined as person with sudden onset of fever >38°C and cough or sore throat, at 12 sentinel sites involved in influenza surveillance in Cameroon. Naso-pharyngeal and/or oropharyngeal swabs were collected and put in 2 ml cryovials containing virus transport medium and stored at 4°C before transportation to the Pasteur Centre of Cameroon (PCC) for analysis. All samples were stored anonymously. This study was specifically approved by the Cameroon National Ethics committee (N° 2016/08/798/CE/CNERSH/SP).

Viral RNA was extracted then from 140 µL of clinical samples with Qiagen RNA kit according to the manufacturer's instructions into a final volume of 60 µL elution buffer. Extracts were analyzed for the presence of influenza by a real-time reverse transcription-polymerase chain reaction (RT-PCR) in an ABI Prism 7300 or 7500 thermocycler (Applied Biosystems, Foster City, California, USA). Positive samples for influenza H3N2 with Ct below 30 after subtyping were eligible for sequencing. A convenient sample size of 35 was randomly selected from the flu database based on their geographic origin and distribution over time. Samples co-infected with other influenza subtypes were not included.

## Amplification

The first PCR was performed using Superscript III One Step RT-PCR System (Invitrogen, Carlsbad, USA) and a set of primers (S1 Table). The following procedure was used for
amplification of 1192 bp of HA gene, 1459 bp of NA gene and 1027 bp of M gene of influenza A(H3N2). Overall a total reaction volume of 25 µL contained the following reagents: 12.5 µL of 2x PCR buffer, 0.5 µL of 10 µM forward and reverse primer, 0.12 µL RNAsin 40 U/µL, 0.5µL Superscript RT/Platinum Taq enzyme and 5µL RNA. This mixture was run using the following program for HA and M genes: 45°C for 30min, 55°C for 15min, 94°C for 2min, 30 cycles at (94°C for 45sec, 45°C for 45sec, 72°C for 5min) and 72°C for 5min. For NA gene, the following program was used: 45°C for 30min, 55°C for 15min, 94°C for 2min, 15 cycles at (94°C for 30sec, 55°C for 30sec, 72°C for 2min), 15 cycles at (94°C for 30sec, 45°C for 30sec, 72°C for 3min) and 72°C for 5min. The first PCR product was then used as a template for semi-nested PCR. A total reaction volume of 50 µL contained 38.3 µL of water, 5 µL of 10x PCR buffer, 1 µL of 10 mM dNTP mix, 1.45 mM of 50 mM MgCl₂, 1µL of 10µM forward and reverse primer, 0.25 µL of Taq DNA polymerase 5U/µL and 2µL of first PCR product. This mixture was run using the following program for all gene segments: 94°C for 5min, 30 cycles at (94°C for 45sec, 45°C for 45sec, 72°C for 90sec) and 72°C for 5min. PCR products were run in a gel and bands corresponding to the appropriate size were sent for sequencing at GENEWIZ UK, LTD (Hope End, UK) with the corresponding primer sets.

## Phylogenetic analysis and antigenic characterization

Sequences were edited and assembled with CLC Main Workbench version 5.5. The sequences of reference strains of known clades and northern hemisphere vaccine strains recommended by WHO [18] included in phylogenetic analysis were obtained from the GenBank and GISAID databases. The phylogenetic tree of the coding nucleotide sequences was generated by MEGA version 6.0 using a neighbor-joining method. The number of bootstrap replications was set to 1000 and bootstrap values above 50% were labeled on major tree branches for reference. The amino acid residues in epitopes A to E of influenza A(H3N2) was previously identified [19].

## Nucleotide sequence accession numbers

Influenza A(H3N2) virus sequences included in the analysis were submitted to Genbank under accession numbers KY653817 to KY653919. S2 Table provides detailed information for the Cameroon influenza gene sequences.

## Prediction of glycosylation sites

The prediction of potential N-liked glycosylation sites was performed with an online server: NetNGlyc 1.0. [20]. This server considers the amino acid alignment Asn-X-Ser/Thr, where X can be any amino acid except Asp or Pro. A threshold value of >0.5 suggests glycosylation.

## Prediction of vaccine efficacy using P<sub>epitope</sub> model

The predicted vaccine efficacy of the influenza A(H3N2) seasonal influenza virus was estimated using the P<sub>epitope</sub> method [21–23]. P<sub>epitope</sub> measures the antigenic distance between the vaccine strain and the dominant circulating strains, and it is calculated as follows: (number of mutations in the dominant epitope)/(number of amino acids in the dominant epitope). The mathematical formula linking vaccine efficacy for influenza A(H3N2) and P<sub>epitope</sub> is given by E = -2.47 × P<sub>epitope</sub> + 0.47 for [23]. When the P<sub>epitope</sub> between the vaccine and circulating strains is null, the influenza A(H3N2) vaccine efficacy is 47% [21]. This model correlates best with vaccine efficacy as compared to phylogenetic analyses or antisera hemagglutination inhibition assay [21].
## Results

A total of 103 gene segments from the 35 virus strains were analyzed (35 HA, 35 NA and 33M). Nucleotide and amino acid sequences for the 3 selected genes were compared within the Cameroon virus strains and with the WHO vaccine strains, WHO reference strains and available sequences of other African virus strains.

## Surface glycoprotein genes and gene segments

**Haemagglutinin.** The HA coding sequences and amino acid sequences of the 35 Cameroonian virus strains were analyzed. Comparison of genes was performed on the strains circulating during the 2014 (N = 6), 2015 (N = 17) and 2016 (N = 12) seasons and sequences from the northern hemisphere vaccine and reference strains. Phylogenetic analysis showed that all Cameroonian virus strains had evolved away from the 3C.1-A/Texas/50/2012-like clade by acquiring several amino acid substitutions (Fig 1, Table 1). Virus strains from 2014 clustered with the 2015–2016 vaccine strain (3C.3a-A/Switzerland/9715293/2013) while virus strains from 2015 and 2016 clustered with the 2016–2017 vaccine strain (3C.2a-A/HongKong/4801/2014) except for one isolate (A/Cameroon/15V-8790/2015) which clustered with the reference strain, 3C.3a-A/Switzerland/9715293/2013.

The 2014 virus strains differed from the 2015 and 2016 virus strains at three antigenic sites: antigenic site A, which contained 2 amino acid substitutions (A138S, R142G), antigenic site B showed 3 substitutions (N128T, F159Y and K160T) and antigenic site C showed a Q311H substitution when compared to the 3C.1-A/Texas/50/2012-like clade (Table 1). Isolate A/Cameroon/15V-8790/2015 differed from all other 2015 virus strains in this study by clustering in the same clade as 2014 strains. Isolate A/Cameroon/14V-4680/2014 differed from all other virus strains in this study in having a G263E substitution, and isolate A/Cameroon/16V-3265/2016 differed from all other 2016 strains by possessing N171K and R142G substitutions.

**Neuraminidase.** The NA sequence analysis showed that all Cameroonian strains evolved away from the 3C.1-A/Texas/50/2012-like clade (Fig 2). Strains from 2014 influenza season clustered with the 3C.3-A/Samara/73/2013 strain while 2015 and 2016 virus strains did not cluster with any strain but were closely related to the 2015–2016 and 2016–2017 vaccine strains (3C.3a-A/Switzerland/9715293/2013 and 3C.2a-A/HongKong/4801/2014, respectively).

The NA genes of the 2014 virus strains differed from those of the 2015 and 2016 virus strains by substitutions of several amino acids: E221D and T267K were found in all 2015 and 2016 virus strains. Whereas, I67V, A110D, V313A and V412I substitutions were found only in 2014 virus strains.

**Matrix gene.** Analysis of the M gene showed that 2014 virus strains clustered with the 2016–2017 vaccine strain, 3C.2a-A/HongKong/4801/2014, while 2015 and 2016 virus strains did not cluster with any reference strain (Fig 3). Globally, there was no significant change observed in the M protein of Cameroonian strains with respect to reference strains. One isolate, A/Cameroon/16V-3265/2016, displayed the amino acid substitution M165L compared to the 2016–2017 vaccine strain, A/HongKong/4801/2014.

## Genotypic antiviral drug susceptibility

In order to determine the genotypic antiviral drug susceptibility of the A(H3N2) virus strains included in this study we analyzed the NA and M2 protein coding sequences for the reported genetic markers of resistance against neuraminidase inhibitors (NAIs; oseltamivir, zanamivir, laninamivir, and peramivir) and adamantanes (amantadine and rimantadine). Mutations E41G, R118K, E119V/I/Q, Q136K (but 1 Q136H), D151A/D/E/V, R152K, I222T, R224K, Q226H, E227D, H274Y, E276D, R292K, N294S, S331R, R371K and E119V+I222L/V in the
<figure><img src="image_4.png" /><figcaption>**Fig 1.** Phylogenetic analysis of HA nucleotide sequences of influenza A(H3N2). ◊ represent reference strains of known clades, ♦ represent the northern hemisphere vaccine strains recommended by WHO and ▲ represent Cameroonian strains.</figcaption></figure>

https://doi.org/10.1371/journal.pone.0184411.g001
<table><caption>Table 1. Comparison of Cameroonian virus strains with respect to A/Texas/50/2012 strain.</caption><thead><tr><th rowspan="2">Antigenic sites</th><th colspan="3">Amino acid substitutions</th></tr><tr><th>2014 strains</th><th>2015 strains</th><th>2016 strains</th></tr></thead><tbody><tr><td>A</td><td>A138S, R142G</td><td>-</td><td>-</td></tr><tr><td>B</td><td>N128A, F159S</td><td>N128T, F159Y, K160T</td><td>N128T, F159Y, K160T</td></tr><tr><td>C</td><td>-</td><td>Q311H</td><td>Q311H</td></tr></tbody></table>

https://doi.org/10.1371/journal.pone.0184411.t001

NA protein have been reported to cause reduced susceptibility to NAIs [24–27]. We laid emphasis on the following substitutions: E119I/V, 245–248 deletion, R292K, N294S and E119V+I222V that have specifically been associated with clinical resistance to neuraminidase inhibitors in A(H3N2) viruses [28–32]. None of the strains had a previously reported mutation peculiar of reduced susceptibility to NAIs however 1 isolate had an amino acid substitution at position 136 (Q136H).

Amino acids at positions 26, 27, 30 and 31 of the transmembrane region of the M2 protein are the residues frequently associated with amantadine- and rimantadine-resistant strains [33]. Amino acid sites were checked as follows for mutations associated with reduced susceptibility or resistance to adamantanes: L26F, V27A/T, A30T/V, S31N/R, G34E, A30V, S193F and R45H [25, 26]. All strains possessed the S31N amino acid substitution.

## Prediction of glycosylation sites

There were 9 predicted glycosylation sites in the HA1 of 2014 virus strains as well as one 2015 isolate (A/Cameroon/15V-8790/2015) at amino acid positions 8(NST), 22(NGT), 38(NAT), 45 (NSS), 63(NCT), 133(NGT), 165(NVT), 246(NST) and 285(NGS), whereas 2015 and 2016 virus strains had 2 additional glycosylation sites at position 126(NWT) and 158(NYT) respectively due to A128T and K160T substitutions.

## Estimation of vaccine efficacy for A(H3N2)

To assess the effect of the accumulated mutations in the HA1 domain on predicted vaccine efficacy in a given year, the $P_{epitope}$ method was used to evaluate how closely the vaccine strain resemble the circulating strain (Table 2). Theoretically, when $P_{epitope}$ in the dominant epitope is higher than 0.19, the vaccine efficacy becomes negative [21, 22]. For 2014, the HA1 sequences mostly had a dominant mutation in epitope B (128, 159, 186, 198) and the $P_{epitope}$ of 0.191 with respect to A/Texas/50/2012 vaccine strain, suggesting that the latter poorly matched the multiple circulating strains present that year consequently leading to a negative vaccine efficacy against these strains of -0.11% (E = -0.05% of 47%, $P_{epitope}$ = 0). For the 2015 and 2016 seasons, HA1 sequences showed antigenic drift mainly in epitopes A. The $P_{epitope}$ of 0.1579 for 2015 strains (dominant epitope = A; substitutions: 138, 142, 144) predicted a vaccine efficacy against these strains of 17.02% (E = 8% of 47%, $P_{epitope}$ = 0) of that of a perfect match with the A/Switzerland/9715293/2013 vaccine strain. For 2016, $P_{epitope}$ of 0.0526 from 9 strains (dominant epitope = A; mutation 144) predicted a minimum vaccine efficacy against these strains of 72.36% (E = 34.01% of 47%, $P_{epitope}$ = 0) of that of a perfect match with the A/HongKong/4801/2014 vaccine strain. Influenza A(H3N2) demonstrated antigenic drift from epitopes B to epitopes A.

The 2015–2016 vaccine strain A/Switzerland/9715293/2013, which belongs to clade 3C.3, was also examined for its ability to match the A(H3N2) strains of the same clade circulating in 2014. The resulting $P_{epitope}$ value of 0 gave an estimated worst-case vaccine efficacy against these strains of 100%. Moreover, clade 3C.2 strains of A(H3N2) in 2015 season and the 2016–
<figure><img src="image_5.png" /><figcaption>**Fig 2.** Phylogenetic analysis of NA nucleotide sequences of influenza A(H3N2). ◊ represent reference strains of known clades, ♦ represent the Northern hemisphere vaccine strains recommended by WHO and ▲ represent Cameroonian strains.</figcaption></figure>

https://doi.org/10.1371/journal.pone.0184411.g002
<figure><img src="image_6.png" /><figcaption>**Fig 3.** Phylogenetic analysis of M nucleotide sequences of influenza A(H3N2). ◊ represent reference strains of known clades, ♦ represent the Northern hemisphere vaccine strains recommended by WHO and ▲ represent Cameroonian strains.</figcaption></figure>

https://doi.org/10.1371/journal.pone.0184411.g003
<table><caption>Table 2. Efficacy among the vaccine strains and number of mutations found on the dominant epitopes of influenza A(H3N2) circulating in Cameroon.</caption><thead><tr><th>Vaccine strain</th><th>No. of strains</th><th>Dominant epitope</th><th>No. of mutations (residue differences)</th><th>P<sub>epitope</sub></th><th>Predicted vaccine efficacy (%)</th></tr></thead><tbody><tr><td>3C.1-A/Texas/50/2012 (2014–2015 vaccine strain)</td><td>6</td><td>B</td><td>4 (128, 159, 186, 198)</td><td>0.1905</td><td>-0.11</td></tr><tr><td rowspan="3">3C.3a-A/Switzerland/9715293/2013 (2015–2016 vaccine strain)</td><td>15</td><td>A</td><td>3 (138, 142, 144)</td><td>0.1579</td><td>17.02</td></tr><tr><td>1</td><td>B</td><td>1 (198)</td><td>0.0476</td><td>74.98</td></tr><tr><td>1</td><td>A</td><td>4 (135, 138, 142, 144)</td><td>0.2105</td><td>-10.62</td></tr><tr><td rowspan="4">3C.2a-A/HongKong/4801/2014 (2016–2017 vaccine strain)</td><td>9</td><td>A</td><td>1 (144)</td><td>0.0526</td><td>72.36</td></tr><tr><td>1</td><td>A</td><td>1 (142)</td><td>0.0526</td><td>72.36</td></tr><tr><td>1</td><td>A</td><td>2 (142, 144)</td><td>0.1053</td><td>44.66</td></tr><tr><td>1</td><td>E</td><td>1 (261)</td><td>0.0455</td><td>76.09</td></tr></tbody></table>

https://doi.org/10.1371/journal.pone.0184411.t002

2017 vaccine strain, A/HongKong/4801/2014, showed a value of 0.0526, indicating a worst-case vaccine efficacy against these strains of 72.36%.

Discussion

This is the first report of genome characterization of influenza H3N2 in Cameroon. We analyzed the coding sequences of 3 gene segments of influenza A/H3N2 virus strains from 2014 to 2016 influenza seasons in Cameroon and compared them with those of the WHO reference strains, WHO recommended vaccine strains [18] and with other epidemiologically relevant African strains from both GenBank and GISAID databases. We observed different phylogenetic clusters generated for different gene segments.

The HA gene sequences indicated that all the Cameroon strains had evolved away from the 3C.1-A/Texas/50/2012-like clade by acquiring the genetic markers N145S, Y186G, P198S and F219S which has been identified as being characteristic of clades 3C.2 and 3C.3 [13, 34]; these markers were reported in most other H3N2 influenza viruses isolated in Africa between 2014–2016. There was no cluster between the Cameroonian strains and the recommended 2014–2015 reference vaccine strain for the northern hemispheres, 3C.1-A/Texas/50/2012-like clade. None of the selected African virus strains circulating between 2014 and 2016 clustered in the 3C.1-A/Texas/50/2012 clade either. Cameroonian strains formed two distinct clusters: 2014 virus strains clustered with the 3C.2a-A/Switzerland/9715293/2013-like clade while the 2015 and 2016 virus strains clustered with the 3C.2a-A/Hong Kong/4801/2014-like clade. Interestingly, the Cameroonian NA gene sequences had also evolved away from the 3C.1-A/Texas/50/2012-like clade by acquiring several amino acid substitutions. A unique 2015 isolate, A/Cameroon/15V-8790/2015, was comparable in many regards to the 2014 virus strains; however, its M and NA genes were unlike the remaining 2015 virus strains.

Theoretically, a clear H3N2 vaccine mismatch was observed during the 2014–2015 influenza seasons as reported by Tewawong *et al.* in the Southern hemisphere [22]. Several amino acid mutations on epitope A of the HA gene were responsible for the drift between the circulating strains and the 2014–2015 vaccine strains, 3C.1-A/Texas/50/2012-like strain. This corroborates with the P<sub>epitope</sub> > 0.19 observed, and eventually resulting in a negative vaccine efficacy. The vaccine change to the A/Switzerland/9715293/2013-like strain for the following season resulted in more than 80% of the circulating H3N2 strains not covered by the selected vaccine. Notwithstanding the poor results with the previous influenza seasons, the 2016–2017 selected vaccine strain showed moderate efficacy when compared to the circulating virus strains. It is noteworthy that vaccination against influenza virus, is not a common
phenomenon in the Cameroon population, however a better selection of annual vaccine with the use of $P_{epitope}$ as suggested by other authors, [21] will improve vaccine efficacy and encourage uptake of flu vaccine. This may result in reduced burden caused by influenza due to herd immunity [35].

The acquisition of a mutation in influenza A viruses conferring resistance to an antiviral agent may occur as a result of drug selection, spontaneous mutation or through genetic reas sortment with another drug resistant influenza A virus [26]. Neuraminidase inhibitors target the NA glycoprotein which plays a key role in the production of new viral progeny. There have been rare reports of resistance to neuraminidase inhibitors in influenza H3N2 virus strains [27], and this was confirmed Cameroonian virus strains by the absence of previous mutations characteristic of this resistance. There was however a novel amino an amino acid mutation, Q136H, that had not been previously reported. Further investigation on the clinical relevance of this mutation will be essential for a better understanding of the molecular evolution of influ enza A(H3N2).

Amantadine resistance reported in most H3N2 viruses worldwide was confirmed in the Cameroonian virus strains by an S31N substitution in the M2 protein [25, 27]. High prevalence of amantadine resistance in influenza A viruses was observed in other countries, irrespective of its use [25, 26]. Resistance to M2 inhibitors first appeared following extensive drug use in Asia and the US after the SARS (Severe Acute Respiratory Syndrome) epidemic in 2004. Since there is no history of exposure to anti-viral drugs against influenza and eventually to adamantanes, Cameroonian virus strains with S31N mutation indicated that strains already harbouring drug resistant mutation were most likely introduced in the community. These resistant strains probably have a selective advantage over the sensitive strains that favour their evolution and fitness [36].

As reported by Mostafa *et al.*, glycosylation and deglycosylation are important viral mechanisms to (i) mask antigenic epitopes and thus inhibit binding to neutralizing antibodies, (ii) adjust receptor-binding affinity, or (iii) regulate virulence of influenza [14]. The Cameroonian H3N2 virus strains shared at least nine potential N-glycosylation sites with the 2014–2017 vaccine strains, and had similar glycosylation patterns. Group 3C.2a-like strains had acquired 2 additional glycosylation sites: 126NWT128 and 158NYT160. The predicted HA glycosylation patterns may contribute to a potential antigenic escape of influenza H3N2 viruses from antibodies raised against vaccine strain and previous circulating strains or increase viral fitness by yet-to-be-characterized mechanisms.

**Conclusion**

The presence of several antigenic site mutations with high frequency among H3N2 virus strains between 2014–2016 influenza seasons confirms the progressing evolution of circulating H3N2 strains. These results lay more emphasis on the relevance of routine influenza surveillance as well as the genetic and antigenic characterization in facilitating the prompt and efficient selection of the influenza strains which will be included in the flu vaccine. This data will be very valuable in tracking the phylogenetic evolution of influenza viruses within sub-Saharan Africa.

**Supporting information**

S1 Table. Set of primers for PCR and sequencing.
(PDF)

S2 Table. Identification of Cameroon H3N2 virus strains.
(PDF)
**Acknowledgments**

We are grateful to Yong Poovorawan and Nipaporn Tewawong for their helpful assistance in P<sub>epitope</sub> calculation. We are also thankful to Matthieu Schoenhals for editing this manuscript.
This work was supported by the U.S. Department of Health and Human Services (DHHS) grant number 6 DESP060001-01-01 via the International Network of Pasteur Institutes.

**Author Contributions**

**Conceptualization:** Marie-Astrid Vernet, Kathleen Victoir, Jane Francis Akoachere, Damian Anong, Richard Njouom.

**Formal analysis:** Gwladys C. Monamele, Mohammed R. Njankouo.

**Funding acquisition:** Kathleen Victoir, Richard Njouom.

**Investigation:** Gwladys C. Monamele, Mohammed R. Njankouo.

**Methodology:** Gwladys C. Monamele.

**Resources:** Kathleen Victoir, Richard Njouom.

**Supervision:** Marie-Astrid Vernet, Kathleen Victoir, Jane Francis Akoachere, Damian Anong.

**Validation:** Richard Njouom.

**Writing – original draft:** Gwladys C. Monamele.

**Writing – review & editing:** Gwladys C. Monamele, Marie-Astrid Vernet, Mohammed R. Njankouo, Kathleen Victoir, Jane Francis Akoachere, Damian Anong, Richard Njouom.

**References**

1. WHO. Influenza (Seasonal) Mediacentre: WHO, 2016 Contract No.: January 2017.

2. Kilbourne ED. Influenza pandemics of the 20th century. Emerg Infect Dis. 2006; 12(1):9–14. https://doi.org/10.3201/eid1201.051254 PMID: 16494710

3. Dawood FS, Iuliano AD, Reed C, Meltzer MI, Shay DK, Cheng P-Y, et al. Estimated global mortality associated with the first 12 months of 2009 pandemic influenza A H1N1 virus circulation: a modelling study. The Lancet Infectious Diseases. 2012; 12(9):687–95. https://doi.org/10.1016/S1473-3099(12)70121-4 PMID: 22738893

4. McAnerney JM, Treurnicht F, Walaza S, Cohen AL, Tempia S, Mtshali S, et al. Evaluation of influenza vaccine effectiveness and description of circulating strains in outpatient settings in South Africa, 2014. Influenza and other respiratory viruses. 2015; 9(4):209–15. https://doi.org/10.1111/irv.12314 PMID: 25865249

5. WHO. Review of the 2014–2015 influenza season in the northern hemisphere. 2015 Jun 5. Report No.: 0049–8114 (Print) 0049–8114 Contract No.: 23.

6. Bedford T, Riley S, Barr IG, Broor S, Chadha M, Cox NJ, et al. Global circulation patterns of seasonal influenza viruses vary with antigenic drift. Nature. 2015; 523(7559):217–20. https://doi.org/10.1038/nature14460 PMID: 26053121

7. Caini S, Huang QS, Ciblak MA, Kusznierz G, Owen R, Wangchuk S, et al. Epidemiological and virological characteristics of influenza B: results of the Global Influenza B Study. Influenza and other respiratory viruses. 2015; 9 Suppl 1:3–12.

8. Tong S, Zhu X, Li Y, Shi M, Zhang J, Bourgeois M, et al. New World Bats Harbor Diverse Influenza A Viruses. PLoS Pathogens. 2013; 9(10):e1003657. https://doi.org/10.1371/journal.ppat.1003657 PMID: 24130481

9. Wiley DC, Wilson IA, Skehel JJ. Structural identification of the antibody-binding sites of Hong Kong influenza haemagglutinin and their involvement in antigenic variation. Nature. 1981; 289(5796):373–8. PMID: 6162101

10. Shih AC, Hsiao TC, Ho MS, Li WH. Simultaneous amino acid substitutions at antigenic sites drive influenza A hemagglutinin evolution. Proc Natl Acad Sci U S A. 2007; 104(15):6283–8. https://doi.org/10.1073/pnas.0701396104 PMID: 17395716
11. Wilson IA, Cox NJ. Structural basis of immune recognition of influenza virus hemagglutinin. Annual review of immunology. 1990; 8:737-71. https://doi.org/10.1146/annurev.iy.08.040190.003513 PMID: 2188678

12. Eshaghi A, Duvvuri VR, Li A, Patel SN, Bastien N, Li Y, et al. Genetic characterization of seasonal influenza A (H3N2) viruses in Ontario during 2010–2011 influenza season: high prevalence of mutations at antigenic sites. Influenza and other respiratory viruses. 2014; 8(2):250–7. https://doi.org/10.1111/irv.12219 PMID: 24313991

13. 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. International journal of medical microbiology: IJMM. 2015; 305(7):762–75. https://doi.org/10.1016/j.ijmm.2015.08.030 PMID: 26416089

14. Mostafa A, Abdelwhab el SM, Slanina H, Hussein MA, Kuznetsova I, Schuttler CG, et al. Phylogenetic analysis of human influenza A/H3N2 viruses isolated in 2015 in Germany indicates significant genetic divergence from vaccine strains. Arch Virol. 2016; 161(6):1505–15. https://doi.org/10.1007/s00705-016-2815-x PMID: 26973232

15. Hanon E. Vaccination strategies against influenza. Bulletin et memoires de l'Academie royale de medecine de Belgique. 2009; 164(10):283-7. PMID: 20669618

16. Njouom R, Mba SA, Noah DN, Gregory V, Collins P, Cappy P, et al. Circulation of human influenza viruses and emergence of Oseltamivir-resistant A(H1N1) viruses in Cameroon, Central Africa. BMC infectious diseases. 2010; 10:56. https://doi.org/10.1186/1471-2334-10-56 PMID: 20205961

17. Escorcia M, Vázquez L, Méndez ST, Rodríguez-Ropón A, Lucio E, Nava GM. Avian influenza: genetic evolution under vaccination pressure. Virology journal. 2008; 5:15. https://doi.org/10.1186/1743-422X-5-15 PMID: 18218105

18. WHO. WHO recommendations on the composition of influenza virus vaccines 2017. Available from: http://www.who.int/influenza/vaccines/virus/recommendations/en/.

19. Lee MS, Chen JS. Predicting antigenic variants of influenza A/H3N2 viruses. Emerg Infect Dis. 2004; 10(8):1385–90. https://doi.org/10.3201/eid1008.040107 PMID: 15496238

20. Prediction of N-glycosylation sites in human proteins. [Internet]. 2004. Available from: http://www.cbs.dtu.dk/services/NetNGlyc/.

21. Gupta V, Earl DJ, Deem MW. Quantifying influenza vaccine efficacy and antigenic distance. Vaccine. 2006; 24(18):3881–8. https://doi.org/10.1016/j.vaccine.2006.01.010 PMID: 16460844

22. 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. https://doi.org/10.1371/journal.pone.0139958 PMID: 26440103

23. Tewawong N, Suntronwong N, Vichiwattana P, Vongpunsawad S, Theamboonlers A, Poovorawan Y. Genetic and antigenic characterization of hemagglutinin of influenza A/H3N2 virus from the 2015 season in Thailand. Virus genes. 2016; 52(5):711–5. https://doi.org/10.1007/s11262-016-1347-5 PMID: 27146171

24. Carr S, Ilyushina NA, Franks J, Adderson EE, Caniza M, Govorkova EA, et al. Oseltamivir-resistant influenza A and B viruses pre- and postantiviral therapy in children and young adults with cancer. Pediatr Infect Dis J. 2011; 30(4):284–8. https://doi.org/10.1097/INF.0b013e3181ff863b PMID: 21048522

25. Lee HK, Tang JW, Loh TP, Hurt AC, Oon LL, Koay ES. Molecular surveillance of antiviral drug resistance of influenza A/H3N2 virus in Singapore, 2009–2013. PLoS One. 2015; 10(1):e0117822. https://doi.org/10.1371/journal.pone.0117822 PMID: 25635767

26. Potdar VA, Dakhave MR, Kulkarni PB, Tikhe SA, Broor S, Gunashekaran P, et al. Antiviral drug profile of human influenza A & B viruses circulating in India: 2004–2011. The Indian journal of medical research. 2014; 140(2):244–51. PMID: 25297358

27. Zaraket H, Kondo H, Hibino A, Yagami R, Odagiri T, Takemae N, et al. Full Genome Characterization of Human Influenza A/H3N2 Isolates from Asian Countries Reveals a Rare Amantadine Resistance-Conferring Mutation and Novel PB1-F2 Polymorphisms. Frontiers in microbiology. 2016; 7:262. https://doi.org/10.3389/fmicb.2016.00262 PMID: 27014195

28. Simon P, Holder BP, Bouhy X, Abed Y, Beauchemin CAA, Boivin G. The I222V Neuraminidase Mutation Has a Compensatory Role in Replication of an Oseltamivir-Resistant Influenza Virus A/H3N2 E119V Mutant. J Clin Microbiol. 2011; 49(2):715–7. https://doi.org/10.1128/JCM.01732-10 PMID: 21106781

29. Okomo-Adhiambo M, Demmler-Harrison GJ, Deyde VM, Sheu TG, Xu X, Klimov AI, et al. Detection of E119V and E119I mutations in influenza A (H3N2) viruses isolated from an immunocompromised patient: challenges in diagnosis of oseltamivir resistance. Antimicrobial agents and chemotherapy. 2010; 54(5):1834–41. https://doi.org/10.1128/AAC.01608-09 PMID: 20194700
30. Baz M, Abed Y, McDonald J, Boivin G. Characterization of multidrug-resistant influenza A/H3N2 viruses shed during 1 year by an immunocompromised child. Clin Infect Dis. 2006; 43(12):1555-61. https://doi.org/10.1086/508777 PMID: 17109288

31. Abed Y, Baz M, Boivin G. Impact of neuraminidase mutations conferring influenza resistance to neuraminidase inhibitors in the N1 and N2 genetic backgrounds. Antiviral therapy. 2006; 11(8):971-6. PMID: 17302366

32. Abed Y, Baz M, Boivin G. A novel neuraminidase deletion mutation conferring resistance to oseltamivir in clinical influenza A/H3N2 virus. The Journal of infectious diseases. 2009; 199(2):180-3. https://doi.org/10.1086/595736 PMID: 19046066

33. Horm SV, Mardy S, Rith S, Ly S, Heng S, Vong S, et al. Epidemiological and Virological Characteristics of Influenza Viruses Circulating in Cambodia from 2009 to 2011. PLoS One. 2014; 9(10).

34. Chambers Benjamin S, Parkhouse K, Ross Ted M, Alby K, Hensley Scott E. Identification of Hemagglutinin Residues Responsible for H3N2 Antigenic Drift during the 2014-2015 Influenza Season. Cell Reports. 2015; 12(1):1-6. https://doi.org/10.1016/j.celrep.2015.06.005 PMID: 26119736

35. Hodgson D, Baguelin M, van Leeuwen E, Panovska-Griffiths J, Ramsay M, Pebody R, et al. Effect of mass paediatric influenza vaccination on existing influenza vaccination programmes in England and Wales: a modelling and cost-effectiveness analysis. The Lancet Public health. 2017; 2(2):e74-e81. https://doi.org/10.1016/S2468-2667(16)30044-5 PMID: 28299371

36. Melidou A, Kyriazopoulou V, Diza E, Alexiou S, Pierroutsakos Y. Antiviral resistance of influenza A (H3N2) strains isolated in northern Greece between 2004 and 2007. Euro surveillance: bulletin Europeen sur les maladies transmissibles = European communicable disease bulletin. 2009; 14(4).