# Diversity of A(H5N1) clade 2.3.2.1c avian influenza viruses with evidence of reassortment in Cambodia, 2014-2016

Annika Suttie, Songha Tok, Sokhoun Yann, Ponnarath Keo, Srey Viseth Horm, Merryn Roe, Matthew Kaye, San Sorn, Davun Holl, Sothyra Tum, et al.

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

## To cite this version:

Annika Suttie, Songha Tok, Sokhoun Yann, Ponnarath Keo, Srey Viseth Horm, et al. Diversity of A(H5N1) clade 2.3.2.1c avian influenza viruses with evidence of reassortment in Cambodia, 2014-2016. PLoS ONE, 2019, 14 (12), pp.e0226108. 10.1371/journal.pone.0226108. pasteur-03214760

HAL Id: pasteur-03214760

https://riip.hal.science/pasteur-03214760

Submitted on 2 May 2021

HAL is a multi-disciplinary open access archive for the deposit and dissemination of scientific research documents, whether they are published 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_2.png" /></figure>
RESEARCH ARTICLE

# Diversity of A(H5N1) clade 2.3.2.1c avian influenza viruses with evidence of reassortment in Cambodia, 2014-2016

Annika Suttie¹,²,³, Songha Tok¹, Sokhoun Yann¹, Ponnarath Keo¹, Srey Viseth Horm¹, Merryn Roe³, Matthew Kaye³, San Sorn⁴, Davun Holl⁴, Sothyra Tum⁴, Philippe Buchy¹,⁵, Ian Barr³, Aeron Hurt³, Andrew R. Greenhill², Erik A. Karlsson¹, Dhanasekaran Vijaykrishna⁶, Yi-Mo Deng³, Philippe Dussart¹,†*, Paul F. Horwood¹,⁷‡*

1 Virology Unit, Institut Pasteur du Cambodge, Institut Pasteur International Network, Phnom Penh, Cambodia, 2 School of Health and Life Sciences, Federation University, Churchill, Australia, 3 WHO Collaborating Centre for Reference and Research on Influenza, Victorian Infectious Diseases Reference Laboratory, at the Peter Doherty Institute for Infection and Immunity, Melbourne, Victoria, Australia,

4 National Animal Health and Production Research Institute, General Directorate of Animal Health and Production, Cambodian Ministry of Agriculture, Forestry and Fisheries, Phnom Penh, Cambodia,

5 GlaxoSmithKline Vaccines R&D Intercontinental, Singapore, Singapore, 6 Department of Microbiology, Biomedicine Discovery Institute, Monash University, Melbourne, Victoria, Australia, 7 College of Public Health, Medical and Veterinary Sciences, James Cook University, Townsville, Queensland, Australia

‡ These authors contributed equally as joint senior authors on this study

* paul.horwood@jcu.edu.au (PFH); pdussart@pasteur-kh.org (PD)

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

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

## OPEN ACCESS

**Citation:** Suttie A, Tok S, Yann S, Keo P, Horm SV, Roe M, et al. (2019) Diversity of A(H5N1) clade 2.3.2.1c avian influenza viruses with evidence of reassortment in Cambodia, 2014-2016. PLoS ONE 14(12): e0226108. https://doi.org/10.1371/journal.pone.0226108

**Editor:** Charles J. Russell, St. Jude Children's Research Hospital, UNITED STATES

**Received:** August 7, 2019

**Accepted:** November 18, 2019

**Published:** December 9, 2019

**Copyright:** © 2019 Suttie 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 manuscript and its Supporting Information files.

**Funding:** This publication is the result of work conducted under a cooperative agreement with the Office of the Assistant Secretary for Preparedness and Response in the U.S. Department of Health and Human Services (HHS), grant number IDSEP140020-01-00. Its contents and conclusions are solely the responsibility of the authors and do not represent the official views of HHS. The study

## Abstract

In Cambodia, highly pathogenic avian influenza A(H5N1) subtype viruses circulate endemically causing poultry outbreaks and zoonotic human cases. To investigate the genomic diversity and development of endemicity of the predominantly circulating clade 2.3.2.1c A (H5N1) viruses, we characterised 68 AIVs detected in poultry, the environment and from a single human A(H5N1) case from January 2014 to December 2016. Full genomes were generated for 42 A(H5N1) viruses. Phylogenetic analysis shows that five clade 2.3.2.1c genotypes, designated KH1 to KH5, were circulating in Cambodia during this period. The genotypes arose through multiple reassortment events with the neuraminidase (NA) and internal genes belonging to H5N1 clade 2.3.2.1a, clade 2.3.2.1b or A(H9N2) lineages. Phylogenies suggest that the Cambodian AIVs were derived from viruses circulating between Cambodian and Vietnamese poultry. Molecular analyses show that these viruses contained the hemagglutinin (HA) gene substitutions D94N, S133A, S155N, T156A, T188I and K189R known to increase binding to the human-type α2,6-linked sialic acid receptors. Two A (H5N1) viruses displayed the M2 gene S31N or A30T substitutions indicative of adamantane resistance, however, susceptibility testing towards neuraminidase inhibitors (oseltamivir, zanamivir, lananmivir and peramivir) of a subset of thirty clade 2.3.2.1c viruses showed susceptibility to all four drugs. This study shows that A(H5N1) viruses continue to reassort with other A(H5N1) and A(H9N2) viruses that are endemic in the region, highlighting the risk of introduction and emergence of novel A(H5N1) genotypes in Cambodia.
was also funded, in part, by the US Agency for International Development (grant No. AID-442-G-14-00005) and partially funded through the UK Research and Innovation Global Challenges Research Fund to The Consortium of Animal Market Networks to Assess Risk of Emerging Infectious Diseases Through Enhanced Surveillance (CANARIES; grant No. GCRFNGR3 \1497). Annika Suttie is funded by an Australian Government Research Training Program Scholarship and a Faculty of Science and Technology Research Scholarship from Federation University. The Melbourne WHO Collaborating Centre for Reference and Research on Influenza is supported by the Australian Government Department of Health. GlaxoSmithKline Biologicals SA provided support in the form of salary for an author [PB], but did not have any additional role in the study design, data collection and analysis, decision to publish, or preparation of the manuscript. The specific role of this author is articulated in the 'author contributions' section. The authors are solely responsible for final content and interpretation.

**Competing interests:** GlaxoSmithKline Biologicals SA provided support in the form of salary for an author [PB]. This does not alter our adherence to PLOS ONE policies on sharing data and materials.

## Introduction

Avian influenza viruses (AIVs; family Orthomyxoviridae) continuously circulate globally in their natural reservoir, wild waterbirds (Orders Anseriformes and Charadriiformes). Frequent spill over and establishment in domestic poultry species increases risk of zoonotic infection. The majority of AIVs circulate as low pathogenic avian influenza viruses that cause mild or no disease in domestic poultry [1,2]. However, subtype A(H5) and A(H7) AIVs are also capable of mutating to form highly pathogenic avian influenza (HPAI) viruses that can cause high morbidity and mortality in poultry flocks [3,4]. The first outbreak of A(H5N1) HPAI viruses occurred among geese in Guangdong, China during 1996 [5]. Since that time, the A/goose/ Guangdong/1/1996 (Gs/Gd) A(H5N1) lineage has disseminated globally causing intermittent outbreaks in domestic poultry and sporadic human infections [6,7]. As of September 2019, 861 confirmed human cases of A(H5N1) infection, resulting in 455 fatalities (case fatality rate, CFR: 53%) have been reported from 17 countries [8]. The primary risk factor for human infection with AIVs is close contact to poultry [9,10] and in vivo studies have shown that A(H5) viruses (with as few as five amino acid substitutions) can acquire aerosol transmissibility in ferrets [11,12]. Fortunately, sustained transmission of A(H5) AIVs between humans has not been documented, though mutations enabling greater transmissibility among humans greatly increases the pandemic threat [13,14].

Influenza A viruses consist of eight negative sense single-stranded RNA segments, each encoding one or more viral proteins. Influenza A viruses are subtyped based on the hemagglutinin (HA) and neuraminidase (NA) glycoproteins that are present on the surface of the viral envelope. There have been eighteen HA subtypes (H1-H18) and eleven NA subtypes (N1-N11) identified. Subtypes H1-H16 and N1-N9 have mainly been identified in avian species, whereas H17-H18 and N10-N11 have only been detected in bats. The HA protein, which is responsible for initiating viral infection by binding to sialic acid receptors on the surface of host cells, is particularly important in host restriction and pathogenicity. The HAs of AIVs preferentially bind to alpha 2,3 (α2,3)-linked sialic acid receptors found in the intestinal and respiratory tract of avian species [15]. While human seasonal influenza viruses preferentially bind to α2,6-liked sialic acid receptors predominantly found in the human upper respiratory tract. Amino acid substitutions in the receptor binding pocket of the HA gene, such as Q222L and G224S (H5 numbering), have been associated with a switch in receptor binding preference from avian-type α2,3 to human-type α2,6 receptors [16]. The acquisition of α2,6 specificity is of concern as it increases AIV transmissibility in mammalian species, increasing their pandemic potential [11,12].

AIVs are a major concern in Southeast Asia. In Cambodia, AIV outbreaks can have devastating socioeconomic consequences as a large proportion of Cambodians rely on agriculture for their livelihoods [17]. A(H5N1) HPAI viruses were first detected in Cambodia during an AIV outbreak in poultry during January 2004 [6]. Since then, A(H5N1) has become endemic in the country, resulting in 56 human infections with 37 deaths (CFR 66%, June 2019) [18]. AIV surveillance in Cambodian poultry has been conducted since 2006 and consists of active surveillance in prominent live bird markets (LBMs) and passive detection following investigations of disease outbreaks in poultry. Human monitoring for zoonotic AIVs consists of a country-wide influenza-like-illness sentinel system, severe acute respiratory illness and event-based surveillance [19,20].

Based on the A(H5) HA clade nomenclature described by the WHO/OIE/FAO H5 Evolution Working Group [21], all Cambodian A(H5N1) viruses identified prior to 2014 belonged to clade 1, or its associated subclades (1.1, 1.1.1 and 1.1.2) [22]. In 2013, a reassortant clade 1.1.2 A(H5N1) virus emerged in Cambodia. The HA and NA genes were from clade 1.1.2,
whereas the MP and all internal genes clustered with clade 2.3.2.1a [21-23]. This reassortant
caused numerous outbreaks in poultry and was associated with a dramatic increase in human
cases during 2013 (n = 26) [23] and early 2014 (n = 8) [24]. The reassortant virus was subse-
quently replaced by a clade 2.3.2.1c virus after March 2014. Interestingly, only one human case
has been detected since this time. Here, we seek to understand the diversity and molecular evo-
lution of A(H5N1) clade 2.3.2.1c viruses in Cambodia, particularly in regard to human and
avian disease risk. This study reports the genetic diversity and molecular evolution of Cambo-
dian A(H5N1) clade 2.3.2.1c viruses detected from 2014 to 2016 through routine surveillance
in LBMs, from poultry outbreak investigations and human cases.

**Materials and methods**

Ethical approval

Animal sampling was conducted by the National Animal Health and Production Research Institute (NAHPRI) under the direction of the General Directorate for Animal Health and Production, Cambodian Ministry of Agriculture, Forestry and Fisheries as part of routine disease surveillance activities; thus, poultry sampling was not considered as experimental animal research. The analysis of poultry samples for avian influenza testing was approved by the Cambodian National Ethics Committee for Health Research (approval #051NECHR). The Institut Pasteur du Cambodge (IPC) serves as a World Health Organization H5 Reference Laboratory and the Cambodian National Influenza Center, with approvals and infrastructure necessary to work on highly pathogenic avian influenza. No animal experimentation was performed at IPC.

Sample collection

Samples included in this analysis were collected as part of passive and active surveillance sys-
tems used to monitor AIV circulation in Cambodia from 2014 to 2016. Surveillance efforts in
poultry were co-ordinated by the Virology Unit at IPC and NAHPRI under the direction of
the General Directorate for Animal Health and Production, Cambodian Ministry of Agricul-
ture, Forestry and Fisheries. Poultry outbreaks of influenza A(H5N1) were investigated by
NAHPRI and positive samples were sent to IPC for confirmation and viral characterization.

Active surveillance was conducted at two prominent Cambodian LBMs: Phnom Penh, the capital city of Cambodia (Orussey market) and a provincial market in Takeo (Takeo market; Fig 1). Active surveillance was performed from 2015 to 2016, but not in 2014. The LBM surveillance strategies and AIV screening methods have been described previously for the 2015 study [25]. The methods used in the 2016 LBM study were the same as described for 2015, however the sampling strategy varied. In 2016, samples were collected solely from Orussey market and collections were performed around three Cambodian festival periods that are known to have high levels of AIV circulation [22]: Lunar New Year (February), Khmer New Year (April) and Pchum Ben (October).

Surveillance for human zoonotic influenza infections was conducted through: A) a country-wide influenza-like-illness sentinel surveillance system; B) severe acute respiratory illness (SARI) surveillance at two major paediatric hospitals (in Phnom Penh and Siem Reap); and C) event based surveillance focused on SARI cases with a history of contact with sick or dead poultry. All samples were screened as previously described [25] using influenza A virus and A (H5N1) qRT-PCR assays available at the International Reagent Resource (https://www.internationalreagentresource.org/Home.aspx).
<figure><img src="image_5.png" /><figcaption>**Fig 1.** Map of Cambodia showing the locations of live bird market sampling sites, clade 2.3.2.1c AIV human cases and poultry outbreaks from 2014 to 2016. Poultry outbreaks were reported in the provinces: Battambang, Kampong Cham, Kampot, Kandal and Siem Reap. The approximate locations of each outbreak (n = 5) are shown by red circles, with the month and year of the outbreak listed below each site. The LBM surveillance sites were Orussey market in Phnom Penh and Takeo market in Takeo; both markets are indicated by black circles. The location of the single human AIV case caused by a clade 2.3.2.1c virus in Tboung Khmum is indicated by a green star. The map was produced using QGIS version 2.18.4 using public domain data obtained from Natural Earth (http://www.naturalearthdata.com/) [26].</figcaption></figure>

https://doi.org/10.1371/journal.pone.0226108.g001

## Viral isolation

Samples that had a high viral load (cycle threshold (CT) < 30 measured by qRT-PCR) were inoculated into embryonated chicken eggs (ECEs) for viral isolation. High viral load samples were prioritised as it becomes increasingly difficult to isolate viruses from samples with high CT values and the supply of ECEs in Cambodia is limited. Prior to inoculation, original samples were diluted 1:1 with antibiotics (a penicillin-streptomycin solution) and passed through a filter (0.22 µM) to prevent bacterial growth. The filtered solution was injected into the allantoic cavity of 10 to 12 day ECEs and incubated at 35°C for 48 to 72 hours in a humid chamber. After incubation, ECEs were chilled for a minimum of 4 hours to kill the embryo and constrict
blood vessels. A hemagglutination assay (HA) was performed on the collected allantoic fluid using 0.5% chicken red blood cells. To confirm the presence of influenza A virus in HA positive samples a qRT-PCR was performed targeting the influenza A MP gene. Samples with a low or negative HA titre (0 or < 8 HA units) were re-passaged in ECEs. Viral isolation was considered unsuccessful if after three successive passages the allantoic fluid in ECEs had no HA titre and tested negative for influenza A via qRT-PCR.

## Sequencing of viral isolates

Influenza A(H5N1) ECE isolates were shipped to the WHO Collaborating Centre for Reference and Research on Influenza, Melbourne, Australia to obtain whole genome sequences using the Ion Torrent Personal Genome Machine™ (Thermo Fisher Scientific, Waltham, MA, USA). RNA was extracted from isolates using the NucleoMag® VET kit (Macherey-Nagel, Bethlehem, USA) that purifies RNA using a magnetic-bead based system. A single RT-PCR reaction was performed to amplify all eight viral genomic segments [27]. The concentration of the PCR products was assessed using Agilent 2200 Tapestation. PCR products were then normalised and fragmented using the Ion Xpress™ Plus Fragment Library Kit (Thermo Fisher Scientific) to produce a 200 base-read library. The fragmented library was purified using the Agencourt® AMPure® XP Kit (Beckman Coulter, Brea, CA, USA) and barcoded using both the Ion Xpress™ P1 Adapter and a specific Ion Xpress™ Barcode (Thermo Fisher Scientific) for each individual sample. The barcoded library was pooled, purified and quantified by qPCR with the Ion Library TaqMan™ Quantitation Kit (Thermo Fisher Scientific). A 20 pM library was used for the emulsion PCR to obtain template positive, enriched Ion Sphere Particles™ (ISPs). The emulsion PCR was performed using the Ion PGM™ HiQ™ View OT2 Kit and run on the Ion OneTouch™ 2 Instrument (Thermo Fisher Scientific). The enriched ISPs were then loaded onto a 318™ Chip v2 (Life Technologies, Carlsbad, CA, USA) and run on the Ion Torrent PGM (Thermo Fisher Scientific).

Quality control of the resulting NGS reads was performed using CLC Genomic Workbench v10 (Qiagen, Valencia, CA, USA). Reads less than 50 nucleotides in length were removed and the remaining sequences were trimmed to only contain reads with a minimum Phred score of 20. High quality reads were then assembled to a reference genome. Any sequences with low coverage or gaps were completed using Sanger sequencing. A combination of segment specific primers and universal primers were used to amplify gene products for Sanger sequencing. PCR products were sequenced using Big Dye Terminator Reaction Mix (Thermo Fisher Scientific) on an ABI 3500xL Genetic Analyzer. Consensus sequences were generated and collated with available NGS data using Geneious® 9.1.8 (Biomatters Ltd, Newark, NJ, USA). All sequencing accession numbers for the viral sequences generated and/or utilized in this study are listed in S1 Table.

## Phylogenetic analysis

AIV sequences used in the phylogenetic analysis were downloaded from GISAID or GenBank [28], curated and aligned using MAFFT v7.308 [29]. Each gene was tested for evidence of recombination using Genetic Algorithm for Recombination Detection (GARD) [30], using the datamonkey webserver [31]. Maximum Likelihood (ML) trees were estimated for each gene segment using IQ-Tree [32], under the General Time Reversible nucleotide substitution model with a gamma rate of heterogeneity (GTR+ Γ). Phylogenetic support was estimated using 1,000 ultrafast bootstrap replicates [33,34]. ML trees were visualized with the graphical viewer FigTree v1.4.3 [35].
A Bayesian phylogenetic analysis using the Markov chain Monte Carlo (MCMC) frame-
work was performed for the HA and NA genes of Cambodian A(H5N1) clade 2.3.2.1c viruses
using BEAST v1.8.4 [36] run using the CIPRES Science Gateway web portal [37]. Viral collec-
tion dates were parsed with variable precision. The analysis was run using an uncorrelated log-
normal relaxed molecular clock [38] with the SRD06 nucleotide substitution model [39]. This
model separates codons into two partitions, one containing the first and second codon posi-
tion and the second partition containing the third codon position. The GTR + Γ substitution
model was then applied to each partition [34]. Change in relative genetic diversity was esti-
mated through time using a Bayesian skyride analysis utilising a Gaussian Markov Random
Field (GMRF) smoothing prior [40]. For each gene, two independent analyses were performed
for 100 million generations sampled to produce 10,000 states. A summary maximum clade
credibility (MCC) tree was produced using TreeAnnotator v1.8.4 with 10% of the burnin
removed. MCC trees with confidence intervals were visualised using FigTree v1.4.3 [35]. The
GMRF plot was produced using Tracer v1.5.

HA gene lineages were defined by using nomenclature from the WHO/OIE/FAO A(H5N1) Evolution Working Group [21]. MP and internal gene lineages were described in a manner similar to previous reports [41–43]. Lineages were assigned based on gene ML phylogenies. Internal gene lineage designations consisted of the genomic segment number followed by a letter specific for each lineage e.g. gene: PB2, lineage A: 1A (S2a-S2g Figs). For genotypes identified, a two-letter country designation was provided (KH for Cambodia) followed by a systematic number. Genotyping was performed in this manner as it is similar to conventions used in Vietnam and allows us to easily compare A(H5) genetic constellations circulating between the two regions.

Reassortment analysis

Phylogenetic congruence was used to investigate the occurrence of viral reassortment in Cam-
bodian A(H5N1) viruses. ML phylogenetic trees were produced for all eight viral segments, as
mentioned previously, and the topological position of each virus was tracked across the phy-
logenies. Incongruence is evident when there is deviation in the phylogenetic topology of indi-
vidual viruses and can indicate putative reassortment events.

Molecular analysis

The Centers for Disease Control and Prevention (CDC) have created a molecular inventory describing mutations and features in A(H5N1) and their overall effect on viral fitness [44]. This inventory was used in conjunction with an updated inventory of AIV molecular markers produced by Suttie et al., 2019 [45] to screen the Cambodian A(H5N1) clade 2.3.2.1c viruses for mutations of interest. Post translation modifications such as N-linked glycosylation sites were predicted for HA and NA using the NetNGlyc1.0 Server [46]. The numbering systems used throughout this text, unless otherwise specified, are: H5 for HA, N1 for NA. Internal genes are numbered relative to mature proteins from A/Vietnam/1203/2004(H5N1), however deletions in NA and internal gene segments are numbered relative to A/goose/Guangdong/1/96 [47].

Selection pressure

To analyse the site specific selection pressures acting on each gene of the Cambodian A
(H5N1) clade 2.3.2.1c viruses the ratio of non-synonomous substitutions (dN) and synonym-
mous substitutions (dS), the dN/dS ratio denoted omega (ω), was computed using the HyPhy
software package [48] accessed through the datamonkey webserver [31]. Selection is
interpreted based on the value of $ω$, when $ω < 1$ sites are considered to be under negative selection, $ω > 1$ indicates positive selection and $ω = 1$ indicates neutrality. To perform this analysis, a combination of four methods: fixed-effects likelihood (FEL), fast unconstrained Bayesian approximation (FUBAR), mixed effects model of evolution (MEME) and single-like-lihood ancestor (SLAC) were used [49–51]. FEL, FUBAR and SLAC estimate sites experiencing pervasive or diversifying selection pressure. Whereas, MEME detects pervasive or episodic positively selected codons. The analysis was based on the nucleotide alignment, the ML phylogeny and the best-fit nucleotide substitution model for each gene. To limit the number of false-positives, statistically significant sites (FEL, MEME and SLAC p-value <0.1 and FUBAR with a posterior probability ≥0.90) that were identified using more than one method were considered valid.

## Neuraminidase activity assay

A NA activity assay was performed to determine the appropriate dilution of each sample to use for the NAI susceptibility testing as described previously [52]. Briefly, gamma irradiated samples were serially diluted 2-fold and MUNANA substrate was added in a 1:1 volume ratio. After an incubation period of 1 hour at 37°C the reaction was stopped by adding a solution of 0.16M NaOH in absolute ethanol. Fluorescence was measured using a Fluoroskan Ascent™ Microplate Fluorometer (Thermo Scientific) with standard filters for excitation (λ 360 nm) and emission (λ 448nm). The NA enzymatic activity was plotted and a viral dilution in the linear region of the curve was used to perform the NAI inhibition assay.

## Neuraminidase inhibition assay

The susceptibility of each virus to four NAI drugs (oseltamivir, zanamivir, laninamivir and peramivir) was quantified by comparing the fluorescence of the uninhibited virus compared to the virus after it had been incubated with varying concentrations of NAI drugs (final drug concentrations ranged from 0.01 nM to 10,000 nM). Briefly, samples were diluted according to NA activity assay results and incubated at room temperature for 45 minutes with the varying drug concentrations. NA activity was then measured as above. The drug concentration that inhibited NA enzymatic activity by 50% was defined as the IC₅₀ value, and this was calculated using the JASPR™ v1.2 software (CDC, Atlanta, GA, USA).

The WHO Influenza Antiviral Working Group have established criteria for determining AIV susceptibility to NAI drugs based on the fold change in IC₅₀ values [53]. Viruses are described as having either: normal (< 10-fold change), reduced (10 to 100-fold increase) or highly reduced (> 100-fold increase) inhibition. Three control viruses (A/Perth/82/2015, A/Osaka/180/2009 and B/Memphis/20/96-R152K) that have established IC₅₀ ranges were included in the NAI assay in duplicate. If the mean IC₅₀ value of each control virus was within the accepted range, the assay was considered valid.

# Results

## Viral isolation

During 2014–2016, 42 A(H5N1) clade 2.3.2.1c viruses were isolated from poultry (21 chickens and 21 ducks) in Cambodian LBMs (S1 Table). Of these, 38 were obtained from an LBM surveillance study conducted in 2015 [54] and four from spot-testing in Cambodian LBMs during 2016. Additionally, partial A(H5N1) genomes from nine environmental samples, collected from containers used to wash poultry carcasses in 2015 were also included in the analysis [25]. In addition, 16 clade 2.3.2.1c viruses were isolated from five poultry outbreaks during passive
surveillance of poultry illnesses and deaths. In total, 67 clade 2.3.2.1c viruses from Cambodian LBMs and poultry outbreaks were analysed in this study (S1 Table), with 42 full genome sequences and 25 partial sequences.

Since the introduction of A(H5N1) clade 2.3.2.1c into Cambodia only one human case has been detected (February 2014) [18]. This virus (A/Cambodia/Y0219302/2014) could not be isolated due to low viral load in the clinical samples. As such, only limited HA and full length NP genes could be generated for this virus (S1 Table). The HA and NP sequences of this human sample were also included in this study.

## Sequence analysis of A(H5N1) clade 2.3.2.1c genes

Sequence and phylogenetic analyses revealed that the Cambodian clade 2.3.2.1c HA genes formed a monophyletic clade with viruses identified in Vietnam from 2012 to 2017. The HA nucleotide sequence identity of all Cambodian 2.3.2.1c viruses from 2014 to 2016 ranged from 96.2 to 100% and were broadly separated into four main groups (Fig 2). HA group 1 was an outlier to the three other groups and contained all clade 2.3.2.1c viruses identified in 2014 (n = 6; including the human case) and a single virus from 2015. HA groups 2 and 3 contained the majority of isolates detected in the 2015 LBM study (n = 12 and 34, respectively). Additionally, a single isolate from the 2016 LBM study clustered with HA group 3. HA group 4 contained all of the A(H5N1) viruses from 2015 and 2016 poultry outbreaks (n = 7 and 4 respectively), as well as the majority (3/4) of the 2016 LBM viruses. The HA groups are closely related; Group 2 is the parental clade to groups 3 and 4, which are sister clades.

Bayesian estimates of the time to most recent common ancestor (TMRCA) for the HA gene of the Cambodian clade 2.3.2.1c viruses indicated that the six groups diverged from a common ancestor around April, 2012 with a 95% Highest Posterior Density (HPD) interval ranging from February to June, 2012 (S1 Fig). The viruses from HA groups 3 and 4, along with viruses sampled in Vietnam during 2014–2017, shared a more recent common ancestor around April 2013 (95% HPD: November 2012 to August 2013).

Similar to the HA gene, the NA and internal genomic segments all clustered with Vietnamese viruses identified from 2012 to 2017. The defined HA groups were also relatively well conserved for NA and the internal genomic segments (S2a-S2g Figs). The nucleotide sequence identity of the Cambodian A(H5N1) genes were as follows; NA: 96.9–100%, PB2: 96.6–100%, PB1: 88.7–100%, PA: 95.8–100%, NP: 94.1–100%, MP: 85.6–100%, NS: 93.2–100%. The lower sequence identity of the PB1 and MP segments were the result of reassortment between A (H5N1) and A(H9N2). Excluding the A(H9N2) reassortant viral segments the sequence identity of PB1 and MP ranged from 97–100% and 95.4–100% respectively.

The human clade 2.3.2.1c virus (A/Cambodia/Y0219302/2014) HA gene was 99.5% identical to that of two viruses detected in poultry around the same period (A/duck/Cambodia/Y0210303/2014 and A/duck/Cambodia/Y0210304/2014). These three viruses differed by only eight bases resulting in two amino acid changes in the human AIV, I151L and F364Y. Neither mutation has been documented to be associated with a significant increase in viral fitness. The NP gene for A/Cambodia/Y0219302/2014 was 99.7% identical to A/duck/Cambodia/Y0224302/2014 and differed by only five synonymous nucleotide changes.

## Detection of reassortment in A(H5N1) clade 2.3.2.1c viruses

The phylogenetic analyses performed for all eight genomic segments showed no evidence of recombination, however there was evidence of numerous reassortment events. Analysis of the genomic constellations shows five genotypes (KH1-KH5) of clade 2.3.2.1c viruses circulated in Cambodia from 2014 to 2016 (Fig 3). Genotypes KH1 and KH2 were originally identified in
<figure><img src="image_6.png" /><figcaption>**Fig 2.** Maximum likelihood phylogenetic tree for the HA gene of A(H5) AIVs. The tree was generated with IQ-Tree using GTR+ Γ and 1,000 ultrafast bootstrap replicates. Cambodian viruses are coloured based on the year they were collected. Viruses detected prior to 2013 are coloured orange, viruses from 2013 are green, viruses from 2014 are purple, 2015 are blue and 2016 are red. Viruses detected during poultry outbreaks are indicated by an asterisk (*) and the single clade 2.3.2.1c Cambodian human A(H5N1) is indicated by a black circle next to the taxa name. Viral clades, HA groups and genotypes are indicated on the right hand side of the tree. Viral clades not listed on the right hand side of the tree have been appended to the end of taxa names. Amino acid differences relative to the closest related WHO candidate vaccine virus (A/duck/Vietnam/NCVD-1584/2012) are shown next to the phylogeny in grey. Mutations listed at branches on the left side of the tree prevail in descendant viruses. Mutations listed next to viral taxa on the right side of the tree are found in the individual virus. Underlined mutations are those that have been previously reported to affect viral fitness. Bootstrap values of 70 or greater are displayed on nodes.</figcaption></figure>

https://doi.org/10.1371/journal.pone.0226108.g002
<figure><img src="image_7.png" /><figcaption>**Fig 3.** Genotypes of highly pathogenic A(H5N1) viruses detected in circulation in Cambodia from 2014 to 2016.
Genomic segments are colour coded to show their possible ancestry. A(H5) clade 2.3.2.1a is shown in red, clade 2.3.2.1b in blue and clade 2.3.2.1c in black. Segments from A(H9N2) viruses are shown in green. Designated KH genotypes (KH1-KH6), any Vietnamese genotype equivalents (VN52, VN53) and the number of Cambodian isolates detected are shown below the associated genome constellations.</figcaption></figure>

https://doi.org/10.1371/journal.pone.0226108.g003

Vietnamese poultry and originated from reassortment events between clades 2.3.2.1a (A/Hubei/1/2010-like), 2.3.2.1b (A/barn swallow/Hong Kong/1161/2010-like) and 2.3.2.1c viruses [41]. KH1 was the only clade 2.3.2.1c genotype detected in 2014 (n = 2), with a single isolate detected in 2015 and was consistent with the Vietnamese genotype VN52. ML trees showed that all KH1 PB2, PA and NP segments clustered with clade 1.1.2 reassortant viruses detected in Vietnam and Cambodia in 2013 that obtained their NA and internal genomic segments from clade 2.3.2.1a viruses (S2b, S2d and S2e Figs). The KH1 MP genes were also from the A/Hubei/1/2010-like lineage, though only the KH1 viruses detected in 2014 cluster with clade 1.1.2 reassortant viruses (S2f Fig). The phylogenetic congruency map shows the KH1 virus from 2015 acquired its MP gene through a separate reassortment event (Fig 4). The NA, PB1 and NS segments were from the A/barn swallow/Hong Kong/1161/2010-like lineage. The KH1 genotype has not been detected since March, 2015.

KH2 was the most frequently identified genotype in 2015 (n = 31) and the sole genotype detected in 2016 (n = 7; Fig 3). The genomic constellation of KH2 –strictly identical to the
<figure><img src="image_8.png" /><figcaption>**Fig 4. Phylogenetic congruence map of Cambodian A(H5N1) viruses.** ML phylogenetic trees were produced for each genomic segment. The scale bars indicate the number of nucleotide substitution per site. For each Cambodian A(H5N1) virus, links have been drawn to connect phylogenetic position of individual viruses across all eight genomic segments. For visual clarity the links are colour coded based on viral genotypes. Viral designations and genotypes are indicated on the left side of the figure. The KH1 genotype is coloured orange, KH2 is dark grey, KH3 is light grey, KH4 is purple and KH5 in blue. Additionally, the putative ancestry of the HA, NA, MP and internal gene segments are shown by vertical shading of each ML tree. Clade 2.3.2.1a is shaded in green, 2.3.2.1b in brown, 2.3.2.1c in blue and A(H9N2) in mauve.</figcaption></figure>

https://doi.org/10.1371/journal.pone.0226108.g004

Vietnamese genotype VN52 -differed from KH1 only in its NP gene, which was from the A/barn swallow/Hong Kong/1161/2010-like lineage. A single KH2 virus, designated A/duck/ Cambodia/Z934W53M3/2015, had an atypical PB2 gene (relative to this genotype) that clus- tered with clade 1.1.2 reassortant viruses (S2b Fig). The phylogenetic congruency map showed that the KH2 genotype also gave rise to KH5 (n = 1), that acquired a PB1 segment from circu- lating A(H9N2) AIVs (Fig 4).

The second most prevalent genotype in 2015 was KH3 (n = 10), produced from the reas- sortment of clade 2.3.2.1c and clade 1.1.2 reassortant viruses. The NA and PB1 segments clus- ter with clade 2.3.2.1c viruses (A/barn swallow/Hong Kong/1161/2010-like). Whereas the PB2, PA, MP and NS genes all stem from clade 1.1.2 reassortant viruses (S2b, S2d, S2f and S2g Figs). The NP genes of KH3 viruses also stem from clade 2.3.2.1a viruses, forming sister clades with the 1.1.2 reassortant viruses (S2e Fig, group 5A). The KH3 genotype was last detected in March, 2015 though it may have been circulating undetected for a longer period as this geno- type gave rise to KH4 (n = 1) identified in November 2015 (Fig 3). The KH4 genotype differed from KH3 only in the clustering of the MP gene as it was acquired from A(H9N2).
## Spatial phylogenetic patterns

The Cambodian A(H5N1) clade 2.3.2.1c AIVs detected within the LBM system do not show strong spatial segregation. The two main LBM sampling sites in Cambodia were Orussey market and Takeo market, 76 km apart by road (Fig 1). HA groups 2 and 3 that collectively account for 92% (n = 47/51) of the LBM AIVs analysed in this study, all contain viruses from both LBMs with no strong subgrouping based on sampling location. Comparatively, HA groups 1 and 4 contain viruses identified from Orussey market or from poultry outbreaks in multiple Cambodian provinces. HA group 1 contains viruses from two poultry outbreaks in Kampong Cham and Kandal from February 2014. These viruses are closely related, belonging to the same cluster, however they form clearly distinguishable subgroups based on outbreak location. Comparatively, HA group 4 contains viruses from three different poultry outbreaks from 2015 and 2016 with less distinct spatial clusters. This is interesting, considering in some instances these viruses were detected approximately six months and over 400 km apart.

## Phylodynamics

The Bayesian skyride analysis of Cambodian, and closely related Vietnamese, clade 2.3.2.1c A (H5N1) viruses shows fluctuations in genetic diversity from 2012 to 2017. From 2012 to 2013, when clade 2.3.2.1c A(H5N1) was initially detected in Vietnam there was a marked increase in relative genetic diversity (Fig 5) [41]. This is typical when pathogens are newly introduced into a geographical region with a naïve population. From 2013 to 2016, genetic diversity oscillated with peaks occurring each year around December and January. These months have been associated with higher prevalence rates of AIVs in Cambodian LBMs [22,54]. From 2012 to 2015, the effective population size was generally increasing. In 2015, there were more fluctuations in genetic diversity and throughout the year the effective population size was gradually declining. In 2016, after the initial peak around January, genetic diversity became relatively stable, with a slight increase in the effective population size in the latter half of the year.

<figure><img src="image_9.png" /><figcaption>**Fig 5.** Bayesian skyride analysis of Cambodian clade 2.3.2.1c AIV genetic diversity. Genetic diversity was estimated using the Gaussian Markov Random Field (GMRF) model using the HA BEAST dataset. The x-axis measures time in years and the y-axis is an estimate of genetic diversity calculated from $N_{e}$ (effective population size and the generation length in years) shown in log scale. The median estimate of genetic diversity over time is shown as a solid black line and the purple shading represents the 95% HPD intervals.</figcaption></figure>

https://doi.org/10.1371/journal.pone.0226108.g005
## Cambodian A(H5N1) clade 2.3.2.1c viruses display molecular markers of pathogenicity, adaptation and transmission

Markers of increased virulence in mammalian and/or avian models of disease were identified in the PB2, PB1-F2, HA, NP, NA, M1, M2 and NS genes (summarised in Table 1, a full list of markers investigated is available in S2a–S2h Table). All Cambodian A(H5N1) clade 2.3.2.1c viruses analysed in this study (n = 68) possessed HA genes with multibasic cleavage sites characteristic of A(H5) HPAI viruses. HA sequences also contained mutations associated with an increase in sialic acid receptor binding to α2,6 human-type receptors, including: D94N, S133A, S155N, T156A, T188 and K189R. Though, numerous substitutions were also identified

<table><caption>Table 1. Summary of amino acid substitutions in Cambodian A(H5N1) viruses associated with changes in viral fitness.</caption><thead><tr><th>Protein</th><th>Phenotype</th><th>Mutation/Motif</th><th>Cambodian Isolates (%)</th><th>References</th></tr></thead><tbody><tr><td>PB2</td><td>Increased virulence in mice</td><td>I63T</td><td>I (100)</td><td>[57]</td></tr><tr><td rowspan="2">PB1</td><td rowspan="2">Increased polymerase activity, replicative capacity, virulence in mice and contact transmission in guinea pigs</td><td>E627K</td><td>E (99), Q (1)</td><td>[11]</td></tr><tr><td>D701N</td><td>D (100)</td><td>[58]</td></tr><tr><td>PB1-F2</td><td>Increased polymerase activity and virulence in mice</td><td>D622G</td><td>G (100)</td><td>[59]</td></tr><tr><td rowspan="3">PA-X</td><td rowspan="3">Truncations to the 90 aa protein increase AIV pathogenicity in chickens and ducks</td><td rowspan="3">Truncation</td><td>25 aa (4)</td><td rowspan="3">[60,61]</td></tr><tr><td>57 aa (94)</td></tr><tr><td>90 aa (2)</td></tr><tr><td></td><td></td><td></td><td>253 aa (100)</td><td>[62–64]</td></tr><tr><td>HA</td><td>Multibasic cleavage site can increase viral pathogenicity</td><td>Multibasic</td><td>PPRERRRRKR/GLF (1)<br>PQREKRRKR/GLF (1)<br>PQREKRRKR/GLF (1)<br>PQRERRRRKR/GLF (94)<br>PQRERRRRR/GLF (3)</td><td>[65]</td></tr><tr><td rowspan="7">NP</td><td rowspan="7">Increased in specificity for α2,6 human-type receptors</td><td>D94N</td><td>N (99)</td><td>[66]</td></tr><tr><td>S133A</td><td>A (100)</td><td>[67]</td></tr><tr><td>S155N</td><td>N (97)</td><td>[68]</td></tr><tr><td>T188I</td><td>I (1)</td><td>[67]</td></tr><tr><td>K189R</td><td>R (94)</td><td>[68]</td></tr><tr><td>Q222L</td><td>Q (100)</td><td>[55,56]</td></tr><tr><td>G224S</td><td>G (100)</td><td>[68]</td></tr><tr><td rowspan="3">NA</td><td rowspan="3">Increased in specificity for α2,6 human-type receptors, increased transmission in guinea pigs</td><td>T156A</td><td>A (94)</td><td>[68]</td></tr><tr><td>M105V</td><td>V (100)</td><td>[69]</td></tr><tr><td>A184K</td><td>K (100)</td><td>[70]</td></tr><tr><td>M1</td><td>Enhanced virulence in mice</td><td>49–68 deletion</td><td>49–68 deletion (100)</td><td>[71]</td></tr><tr><td rowspan="2">M2</td><td rowspan="2">Enhanced virulence in mice, chickens and ducks</td><td>N30D</td><td>D (100)</td><td rowspan="2">[72]</td></tr><tr><td>T215A</td><td>A (100)</td></tr><tr><td rowspan="4">NS1</td><td rowspan="4">Increased resistance to amantadine and rimantadine</td><td>I43M</td><td>M (100)</td><td>[73]</td></tr><tr><td>A30T</td><td>T (2)</td><td>[74,75]</td></tr><tr><td>S31N</td><td>N (2)</td><td>[76]</td></tr><tr><td></td><td>80–84 aa deletion</td><td>80–84 deletion (100)</td><td>[76]</td></tr><tr><td rowspan="3"></td><td rowspan="3">Decreased antiviral response and increased virulence in mice</td><td>P42S</td><td>S (100)</td><td>[77]</td></tr><tr><td>D87E</td><td>E (77)</td><td>[76]</td></tr><tr><td>L98F</td><td>F (100)</td><td>[78]</td></tr></tbody></table>

https://doi.org/10.1371/journal.pone.0226108.t001
that are associated with a binding preference for α2,3 receptors (S2d Table). The well-known Q222L and G224S substitutions associated with increased binding to α2,6 receptors and decreased α2,3 binding were not observed [55,56].

A number of the mutations identified in the NA, MP and NS genes are conserved in con-
temporary A(H5N1) viruses [44]. For instance, the NA and NS1 stalk deletions, and the NS1
substitutions: P42S, D87E, L98F and I101M residues that have been shown to increase AIV vir-
ulence in mice [76–80]. In the MP gene, all Cambodian A(H5N1) viruses contained the M1
N30D and T215A residues associated with increased mammalian pathogenicity in vivo [72].
No major molecular markers associated with AIV adaptation to mammalian species were
identified in the polymerase proteins, such as PB2 E627K or D701N.

Analysis of selection pressures

The site-specific selection pressures acting on each gene of the Cambodian clade 2.3.2.1c
viruses were analysed using a combination of four methods (FEL, FUBAR, MEME and SLAC)
that calculate omega (ω). The majority of codons in all genes of the Cambodian A(H5N1)
clade 2.3.2.1c viruses were under either neutral or purifying selection pressure (S3 Table). A
total of six codons were detected to be under positive selection, including: HA (169), NA (74),
PB2 (107, 339) and PB1 (54, 739; S3 Table). The HA and NA codons do not lie in major anti-
genic sites. Whereas, the PB1 codons 54 and 739 are part of PB1-PA or PB1-PB2 binding
domains, respectively [81,82]. Substitutions at all the detected codons, excepting PB2 339, have
not been shown to affect viral fitness. The positively selected codon in PB2 at position 339 is
located on the surface of A(H5N1) PB2 at the edge of a putative cap binding site [83]. At this
site Cambodian clade 2.3.2.1c viruses predominantly contained threonine (339T), with a single
isolate from 2015 containing methionine at this site (339M). Reports on the effect of PB2
mutations at position 339 on viral virulence vary and so the effect of mutations at this codon
on viral fitness remain unclear [84–86].

Post translational modifications—N glycosylation

A total of nine N-glycosylation sites were predicted from the HA sequences of the Cambodian clade 2.3.2.1c A(H5N1) viruses (S4 Table). Of these, five glycosylation sites predicted in the HA2 stem domain were identified in all Cambodian viruses, including: ¹¹NNS, ²³NVT, ²⁸⁶NSS, ⁴⁸⁴NGT and ⁵⁴³NGS [87]. In the HA1 subunit four additional glycosylation sites were identified: ¹⁴⁰NSS, ¹⁶⁵NNT, ²³⁶NDT and ²⁷³NCS. The HA N-glycosylation sites have all been identified in A(H5N1) viruses previously [87]. The glycosylation sites in the HA2 stem domain and at position 165 are highly conserved. The ¹⁴⁰NSS glycosylation site emerged in clade 2.3.2.1 viruses and occurs in a HA antigenic site at epitope B [87]. From the Cambodian NA sequences five glycosylation sites were predicted, including: ²⁸NIT, ³⁵NHS, ⁶⁸NSS, ¹²⁶NGT and ²¹⁵NGS (S4 Table). The glycosylation sites ¹²⁶NGT and ²¹⁵NGS were conserved in all Cambodian A(H5N1) NA proteins. The ⁶⁸NSS glycosylation site was predicted for all, barring one, NA segments. These three sites are conserved in A(H5N1) viruses [87]. Whereas, the remaining two glycosylation sites predicted in the Cambodian viruses are less common: ³⁵NHS (n = 13), ²⁸NIT (n = 3) (S4 Table) [87].

Molecular prediction and susceptibility to antiviral drugs

Two Cambodian A(H5N1) viruses belonging to clade 2.3.2.1c had predicted resistance to ada-
mantanes. The first virus, designated A/chicken/Cambodia/Z50W9M1/2015, had an M2 pro-
tein A30T substitution. The second virus, designated A/chicken/Cambodia/Z850W49M1/
2015, was a reassortant with all internal genes from A(H5N1) viruses and an MP from
circulating A(H9N2) viruses, and has been described previously [88]. The A(H9N2) acquired M2 protein contained the S31N substitution known to confer resistance to adamantanes. The A30T and S31N substitutions are well known and many experimental studies have shown these single amino acid substitutions render AIVs highly resistant to adamantanes.

The Cambodian clade 2.3.2.1c viruses did not contain any known markers indicative of resistance to neuraminidase inhibitors (NAI; S2f Table). However, there were two variants detected at position 129: V129I (n = 3) and V129D (n = 1). The mutation V129A was previously reported to mildly decrease AIV susceptibility to zanamivir [89]. The effect of V129I/D substitutions on AIV susceptibility to NAI drugs is unknown. To confirm NAI sensitivity, a subset of thirty isolates were tested against four NAIs (oseltamivir, zanamivir, peramivir and laninamivir) with no resistance detected (S5 Table). All three viruses with V129I substitution were included in the phenotypic testing and no decrease in susceptibility to NAI drugs was observed. The effect of V129D could not be investigated as no isolate was available for this sample. The mean fold change in IC₅₀ (nM) values for oseltamivir, peramivir and zanamivir were 1.0 and for laninamivir 1.18 (S5 Table).

Discussion

AIVs present a threat to agriculture and could potentially cause the next pandemic, therefore it is vital to continually monitor their circulation and evolution. A(H5N1) HPAI viruses were introduced into Cambodia in 2004. They have since become endemic in domestic poultry, causing substantial economic hardship. Continual circulation of A(H5N1) HPAI viruses in Cambodian poultry is also concerning as humans that live and work in close contact to poultry are at risk of zoonotic transmission due to the lack of biosecurity. Therefore, we sought to investigate the phylogenetic and molecular traits of Cambodian A(H5N1) viruses collected from LBM studies, poultry outbreaks and human infections from 2014 to 2016.

Between 2014 to 2016, five A(H5N1) clade 2.3.2.1c poultry outbreaks were reported in Cambodia (Fig 1). In reality the number of AIV outbreaks is likely to be much higher due to reluctance in reporting [90]. Active surveillance was also performed at Cambodian LBMs in 2015 and 2016. Previous studies have demonstrated that the Cambodian LBMs have a high prevalence of AIVs with a diverse range of AIV subtypes detected [22,54,91].

Since March 2014, A(H5N1) viruses detected in Cambodia have exclusively been of clade 2.3.2.1c [25]. Clade 2.3.2.1c viruses are one of three viral subclades (2.3.2.1a-c) that evolved from clade 2.3.2.1 and were initially reported in domestic poultry from Vietnam in 2012 [41,92]. However, reports show that clade 2.3.2.1c began to circulate widely throughout Asia from 2009 onwards [41,92]. Clade 2.3.2.1c viruses have been detected in both poultry and wild waterfowl species. The widespread circulation of clade 2.3.2.1c viruses in wild birds has facilitated their geographical dispersal, with clade 2.3.2.1c viruses reported in the Middle East, Europe and Africa [93–96]. Aside from avian species, this clade has also been detected in mammals, including big cats (a lion and tiger from China) [97,98] and humans (China: A/Hong Kong/6841/2010; Cambodia: A/Cambodia/Y0219302/2014).

In the present study, we investigated the genetic diversity and evolution of Cambodian clade 2.3.2.1c viruses circulating between 2014 and 2016. The close phylogenetic and molecular association of 2.3.2.1c viruses detected in Cambodia and Vietnam suggests that these viruses circulate endemically between the two countries. Similar trends were observed for the introduction and evolution of clade 1 A(H5N1) HPAI viruses that were first detected in Cambodia in 2004 [99]. The close relationship between the Cambodian and Vietnamese viruses is unsurprising considering cross-border trade and farming of poultry between the two countries is known to occur [100]. It is also possible that the transmission of AIVs between Cambodian
and Vietnam is facilitated by the movement of wild birds. However, the role of wild birds in
the introduction of AIVs into Cambodia has not been investigated to date.

The Cambodian clade 2.3.2.1c viruses circulating between 2014 and 2016 were diverse, with five genotypes present over a 3-year period (KH1-KH5). Of these, the KH1 and KH2 genotypes were equivalent to VN52 and VN53, which circulated in Vietnamese poultry in 2013–2015 [42]. Whereas, KH3-KH5 were novel genotypes identified in Cambodia in 2015 that were produced from reassortment between A(H5N1) clades or with A(H9N2) viruses locally. The generation of novel AIV genotypes can have important public health implications. For instance, the acquisition of specific A(H9N2) internal genomic cassettes has been associated with an increase in zoonotic potential of A(H7N9) and A(H10N8) AIVs [101,102]. Reassortment between A(H5) viruses can also produce viruses of concern for human health. Indeed, this was demonstrated in Cambodia in 2013 with the emergence of the novel clade 1.1.2 reassortant viruses that caused numerous human A(H5N1) cases.

The phylogenetic analysis shows limited signs of spatial segregation. Viruses identified at LBMs from 2015 to 2016 cluster closely together in HA groups 2 and 3 (Fig 2). This is unsurprising considering the high degree of poultry movement within the LBM network and that poultry bought at Takeo market are commonly sold to individuals that transport the poultry to the capital city to resell at Orussey market (Fig 1). Spatial segregation was, at times, more evident in viruses detected during poultry outbreaks as seen from the distinct grouping of outbreak samples from Kampong Cham and Kandal, isolated in February, 2014. Comparatively, in outbreak samples detected in 2015 and 2016 spatial segregation was less evident despite the fact that these samples were collected approximately six months apart at sites that were, in some cases, separated by more than 400 km. These findings demonstrate the persistent circulation and spread of highly pathogenic A(H5N1) viruses in Cambodia outside of the LBM network.

As expected, all of the Cambodian A(H5N1) clade 2.3.2.1c viruses had multibasic HA cleavage site motifs characteristic of HPAI viruses. The Cambodian viruses contained six amino acid substitutions that have been shown to increase receptor binding preference for α2,6 human receptors, including: D94N, S133A, S155N, T156A, T188I and K189R [66–68,89]. Typically, five out of six of these mutations were detected. None of the viruses analysed contained the Q222L or G224S substitutions that are associated with mammalian adaptation. Most of the HA molecular sites associated with host specificity that were investigated are indicative of a preference for avian α2,3 receptors (S2d Table). The identified mutations associated with increased specificity for α2,6 are unlikely to switch the receptor specificity of the Cambodian A (H5N1) viruses entirely. No other mutations associated with AIV adaptation to mammals, such as E627K and D701N in PB2 [103,104], were identified in the internal genes of the Cambodian A(H5N1) viruses.

Post translational modifications, such as the N-glycosylation, are important for protein
folding, maturation and biological functionality. In AIVs, the presence of N-glycans on the
HA and NA can affect viral pathogenicity and virulence as N-glycans can shield antigenic sites
enabling viruses to evade detection by the host immune system [105–107]. They have also
been shown to affect HA receptor binding preference and cleavability [108,109]. The T156A
mutation, detected in 64 out of 68 of the Cambodian clade 2.3.2.1c AIVs, removes an impor-
tant glycosylation site in HA. The absence of glycosylation at position 154–156 is common in
viruses that stem from clade 2.3.2.1 [87] and has been shown to increase viral affinity for α2,6
sialic acid receptors [68,110]. Furthermore, multiple independent studies have shown that
mutations at these positions, such as N154D and T156A, are important factors in the transmis-
sibility and pandemic potential of AIVs when combined with Q222L and G224S [11,12].
A small number of codons were determined to be under positive selection pressure, includ-
ing: HA (169), NA (74), PB2 (107, 339) and PB1 (54, 739) (S3 Table). The NA codon 74 is part
of the stalk region, PB1 codon 54 is part of the PB1-PA binding domain, 739 is part of the
PB1-PB2 binding domain [82] and PB2 codon 339 is located at the edge of a PB2 cap binding
site [83]. Substitutions at all the detected codons, excepting PB2 339, have not been shown to
affect viral fitness. At PB2 position 339 Cambodian clade 2.3.2.1c viruses predominantly con-
tain threonine (339T). There are contradictory reports about the effect of 339T on viral replica-
tion and pathogenicity. One report showed that when in combination with 147T and 588T,
339T increased the polymerase activity and replication of H5N1 AIVs in mammalian cells and
increased pathogenicity in mice [84]. Whereas, another report suggested that 339T attenuates
PB2 cap binding capabilities, decreasing viral polymerase activity, replicative efficiency and
decreases H5N1 AIV virulence in mice [85]. A single isolate from 2015 had a methionine sub-
stitution at this position (339M). The T339M mutation has been suggested as an A(H5N1)
adaptation to humans and was shown to enhance the growth in a mammalian cell line, partic-
ularly when in combination with PB2 mutations 249G or 309D [86]. The precise effect of
mutations in PB2 at position 339, particularly 339T, remains unclear.

Vaccination is an effective way to prevent influenza infection [111]. Given that there are no A(H5N1) vaccines commonly used in humans, antiviral prophylaxis and treatment is the best option for limiting A(H5N1) infection and transmission. Therefore, we screened the NA and MP genes from Cambodian A(H5N1) isolates for markers indicative of resistance to the two classes of antiviral drugs available to treat influenza infections: adamantanes (adamantine and rimantadine) and neuraminidase inhibitors (oseltamivir, zanamivir, peramivir and laninamivir). Resistance to adamantanes is widespread and has been reported in seasonal influenza viruses as well as A(H5N1) AIVs [74,112]. In this study, two A(H5N1) clade 2.3.2.1c viruses had M2 amino acid substitutions indicative of resistance to adamantanes, S31N and A30T (Table 1). Whereas none of the viruses had molecular markers indicative of resistance to NAIs. NAI sensitivity was confirmed in a subset of thirty Cambodian A(H5N1) isolates with all viruses tested being highly susceptible to all four NAIs (S5 Table).

The A(H5N1) clade switch that occurred in Cambodia in 2014 suggests that the clade 2.3.2.1c viruses may have better fitness than the clade 1.1.2 reassortant viruses in poultry. However, the decrease in human cases indicates their zoonotic potential may be reduced. The clade 1.1.2 A(H5N1) reassortant viruses caused 64% (n = 34) of the overall Cambodian A(H5N1) human cases during a period of approximately 15 months from 2013 to 2014. Comparatively, only a single human case has been documented in Cambodia as the result of infection with clade 2.3.2.1c, despite nearly 6 years of circulation in the country. A molecular analysis of the clade 1.1.2 A(H5N1) viruses performed previously showed they do not typically contain major markers associated with AIV adaptation to mammals [23]. However, four HA mutations associated with an increase in binding to human-type α2,6 receptors were conserved in the population, including: S123P, S133A, S155N, and K266R [23]. Comparatively, the Cambodian clade 2.3.2.1c A(H5N1) viruses also have the S133A and S155N substitutions in 100% and 97% of viruses, respectively (Table 1). Though they do not have the S123P and K226R substitutions, other molecular markers associated with an increase in α2,6 binding were identified that are not common to clade 1.1.2 AIVs. It is possible differences in the viral genes other than HA contribute to this phenotype. The molecular basis for the decrease in Cambodian clade 2.3.2.1c A(H5N1) transmissibility to humans, compared to the previous circulating clade 1.1.2 viruses remains unclear.

It is important to note that the majority of data available for analysis is from 2015, as limited
LBM sampling was performed in 2014 and 2016. Therefore, it is likely the viral diversity in
these two years was higher than was presented in this study. Analysis of the genetic diversity of
the Cambodian clade 2.3.2.1c A(H5N1) viruses shows fluctuations occurred each year from 2014 to 2016 with peaks typically occurring around December or January (Fig 5). This coincides with the Cambodian dry season that has previously been documented as a period of high AIV prevalence in Cambodia [22,54]. Sampling and analysing the molecular traits of viruses in subsequent years will be crucial to track viral evolution and diversity in Cambodia. Analysing molecular markers associated with changes in viral fitness is useful to rapidly assess the biological characteristics and risk associated with AIVs. However, moving forward it is important to verify the molecular profile of the Cambodian viruses by performing further *in vitro* and *in vivo* virulence testing. Particularly as the effect of substitutions can be context dependent and novel substitutions associated with an increase in viral fitness may emerge. This will provide a more comprehensive assessment of the risk Cambodian A(H5N1) pose to the community.

Overall, in this study we genetically characterised clade 2.3.2.1c A(H5N1) viruses circulating in Cambodia from 2014 to 2016. Analysis shows that reassortment of internal genes between A(H5N1) clades and other AIV subtypes does occur in the region. Continual surveillance and characterisation of AIVs is essential to limit the impact of this disease on animals, the economy and human health.

**Supporting information**

**S1 Figs.** Bayesian maximum clade credibility (MCC) phylogenetic tree of Cambodian clade 2.3.2.1c virus a) HA and b) NA genes detected from 2014 to 2016. Cambodian viruses are coloured based on the year they were detected: viruses from 2014 are purple, 2015 are blue and 2016 are red. The single human Cambodian clade 2.3.2.1c virus is indicated by a black circle. Trees were generated with BEAST v1.84 using GTR+Γ with the SRD06 nucleotide substitution model. The tree branch lengths are time-proportional and the time scale is indicated on the x axis. The proposed tMRCA is displayed at each node.
(PDF)

**S2 Figs.** Maximum likelihood phylogenetic trees for NA, MP and internal genomic segments of Cambodian A(H5N1) viruses detected from 2014 to 2016. a) NA, b) PB2 c) PB1 d) PA e) NP f) MP and g) NS. Trees were generated with IQ-Tree using GTR+Γ and 1,000 ultra-fast bootstrap replicates. Taxa names show viral subtype, HA clade designation and viral strain name. Cambodian viruses are coloured based on the year they were collected. Viruses detected prior to 2013 are coloured orange, viruses from 2013 are green, viruses from 2014 are purple, 2015 are blue and 2016 are red. Segment lineages are indicated on the right hand side of the tree. For NA amino acid differences relative to the closest related WHO candidate vaccine virus (A/duck/Vietnam/NCVD-1584/2012) are shown next to the phylogeny in grey. Mutations listed at branches on the left hand side of the tree prevail in descendant viruses. Mutations listed next to viral taxa on the right hand side of the tree are found in the individual virus. Underlined mutations are those that have been previously reported to affect viral virulence. Bootstrap values of 70 or greater are displayed on nodes.
(PDF)

**S1 Table.** List of Cambodian A(H5N1) viruses detected between 2014 and 2016 that were included in this analysis with details on sample collection, AIV genotypes and sequencing accession numbers.
(XLSX)

**S2 Table.** Molecular inventory of the Cambodian A(H5N1) viruses between 2014 and 2016:
a) PB2, b) PB1, c) PA, d) HA, e) NP, f) NA, g) MP, h) NS.
(XLSX)
**S3 Table.** Selection pressure analysis of the Cambodian A(H5N1) genes using FEL, FUBAR, MEME and SLAC.

(XLSX)

**S4 Table.** Predicted HA and NA N-glycosylation sites of Cambodian A(H5N1) viruses between 2014 and 2016.

(XLSX)

**S5 Table.** Sensitivity of Cambodian A(H5N1) viruses to neuraminidase inhibitors (zanamivir, oseltamivir, peramivir and laninamivir).

(XLSX)

**Acknowledgments**

The authors would like to thank the field team from the National Animal Health and Production Research Institute (General Directorate of Animal Health and Production, Cambodian Ministry of Agriculture, Forestry and Fisheries); and the field and laboratory teams from the Virology Unit at the Pasteur Institute in Cambodia.

**Author Contributions**

**Conceptualization:** Philippe Dussart, Paul F. Horwood.

**Data curation:** Annika Suttie.

**Formal analysis:** Annika Suttie, Dhanasekaran Vijaykrishna, Yi-Mo Deng.

**Funding acquisition:** Ian Barr, Philippe Dussart, Paul F. Horwood.

**Investigation:** Annika Suttie, Songha Tok, Sokhoun Yann, Ponnarath Keo, Srey Viseth Horm, Merryn Roe, Matthew Kaye, Yi-Mo Deng.

**Methodology:** Aeron Hurt, Yi-Mo Deng.

**Project administration:** Srey Viseth Horm, San Sorn, Davun Holl, Sothyra Tum, Andrew R. Greenhill, Philippe Dussart, Paul F. Horwood.

**Resources:** Philippe Buchy.

**Supervision:** Ian Barr, Aeron Hurt, Andrew R. Greenhill, Erik A. Karlsson, Dhanasekaran Vijaykrishna, Yi-Mo Deng, Philippe Dussart, Paul F. Horwood.

**Writing – original draft:** Annika Suttie.

**Writing – review & editing:** Philippe Buchy, Ian Barr, Andrew R. Greenhill, Erik A. Karlsson, Dhanasekaran Vijaykrishna, Yi-Mo Deng, Philippe Dussart, Paul F. Horwood.

**References**

1. Webster RG, Bean WJ, Gorman OT, Chambers TM, Kawaoka Y. Evolution and ecology of influenza A viruses. Microbiol Rev 1992; 56:152–179. PMID: 1579108

2. Kuiken T. Is low pathogenic avian influenza virus virulent for wild waterbirds? Proc R Soc B Biol Sci 2013; 280:20130990–20130990. https://doi.org/10.1098/rspb.2013.0990.

3. Sturm-Ramirez KM, Ellis T, Bousfield B, Bissett L, Dyrting K, Rehg JE, et al. Reemerging H5N1 Influenza Viruses in Hong Kong in 2002 Are Highly Pathogenic to Ducks. J Virol 2004; 78:4892–901. https://doi.org/10.1128/JVI.78.9.4892-4901.2004 PMID: 15078970

4. Swayne DE. Understanding the Complex Pathobiology of High Pathogenicity Avian Influenza Viruses in Birds. Avian Dis 2007; 51:242–9. https://doi.org/10.1637/7763-110706-REGR.1 PMID: 17494560
5. Xu X, Subbarao K, Cox NJ, Guo Y. Genetic characterization of the pathogenic influenza A/Goose/ Guangdong/1/96 (H5N1) virus: similarity of its hemagglutinin gene to those of H5N1 viruses from the 1997 outbreaks in Hong Kong. Virology 1999; 261:15-19. https://doi.org/10.1006/viro.1999.9820 PMID: 10484749

6. World Health Organization, Food and Agriculture Organization of the United Nations. H5N1 highly pathogenic avian influenza: Timeline of major events. 2014.

7. World Organization for Animal Health. Avian Influenza Portal: Update on avian influenza in animals (types H5 and H7) 2019. http://www.oie.int/en/animal-health-in-the-world/update-on-avian-influenza/ 2019/ (accessed August 28, 2019).

8. World Health Organisation. Cumulative number of confirmed human cases for avian influenza A (H5N1) reported to WHO, 2003-2019. 2019.

9. Claas EC, Osterhaus AD, Van Beek R, De Jong JC, Rimmelzwaan GF, Senne DA, et al. Human influ- enza A H5N1 virus related to a highly pathogenic avian influenza virus. The Lancet 1998; 351:472-7.

10. Mounts AW, Kwong H, Izurieta HS, Ho Y, Au T, Lee M, et al. Case-control study of risk factors for avian influenza A (H5N1) disease, Hong Kong, 1997. J Infect Dis 1999; 180:505-508. https://doi.org/ 10.1086/314903 PMID: 10395870

11. Herfst S, Schrauwen EJA, Linster M, Chutinimitkul S, de Wit E, Munster VJ, et al. Airborne Transmis- sion of Influenza A/H5N1 Virus Between Ferrets. Science 2012; 336:1534-41. https://doi.org/10.1126/ science.1213362 PMID: 22723413

12. Imai M, Watanabe T, Hatta M, Das SC, Ozawa M, Shinya K, et al. Experimental adaptation of an influ- enza H5 HA confers respiratory droplet transmission to a reassortant H5 HA/H1N1 virus in ferrets. Nature 2012. https://doi.org/10.1038/nature10831.

13. Ungchusak K, Auewarakul P, Dowell SF, Kitphati R, Auwanit W, Puthavathana P, et al. Probable per- son-to-person transmission of avian influenza A (H5N1). N Engl J Med 2005; 352:333-340. https://doi. org/10.1056/NEJMoa044021 PMID: 15668219

14. Wang H, Feng Z, Shu Y, Yu H, Zhou L, Zu R, et al. Probable limited person-to-person transmission of highly pathogenic avian influenza A (H5N1) virus in China 2008: 8.

15. Costa T, Chaves AJ, Valle R, Darji A, van Riel D, Kuiken T, et al. Distribution patterns of influenza virus receptors and viral attachment patterns in the respiratory and intestinal tracts of seven avian species. Vet Res 2012; 43:28. https://doi.org/10.1186/1297-9716-43-28 PMID: 22489675

16. Matrosovich M, Gambaryan A, Teneberg S, Piskarev V, Yamnikova S, Lvov D, et al. Avian influenza A viruses differ from human viruses by recognition of sialyloligosaccharides and gangliosides and by a higher conservation of the HA receptor-binding site. Virology 1997; 233:224-234. https://doi.org/10. 1006/viro.1997.8580 PMID: 9201232

17. National Institute of Statistics, Ministry of Planning. Cambodia Socio-Economic Survey 2017 2018. https://www.nis.gov.kh/index.php/en/14-cses/12-cambodia-socio-economic-survey-reports (accessed November 9, 2019).

18. World Health Organization. Cumulative number of confirmed human cases of avian influenza A(H5N1) reported to WHO 2019. https://www.who.int/influenza/human_animal_interface/H5N1.cumulative_ table_archives/en/ (accessed September 28, 2019).

19. 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:e110713. https:// doi.org/10.1371/journal.pone.0110713 PMID: 25340711

20. Horwood PF, Karlsson EA, Horm SV, Ly S, Heng S, Chin S, et al. Circulation and characterization of seasonal influenza viruses in Cambodia. Influenza Other Respir Viruses 2019; 13:465-76. https://doi. org/10.1111/irv.12647 PMID: 31251478

21. Smith GJD, Donis RO, World Health Organization/World Organisation for Animal Health/Food and Agriculture Organization (WHO/OIE/FAO) H5 Evolution Working Group. Nomenclature updates result- ing from the evolution of avian influenza A(H5) virus clades 2.1.3.2a, 2.2.1, and 2.3.4 during 2013-2014. Influenza Other Respir Viruses 2015; 9:271-6. https://doi.org/10.1111/irv.12324 PMID: 25966311

22. Horm SV, Tarantola A, Rith S, Ly S, Gambaretti J, Duong V, et al. Intense circulation of A/H5N1 and other avian influenza viruses in Cambodian live-bird markets with serological evidence of sub-clinical human infections. Emerg Microbes Infect 2016; 5:e70. https://doi.org/10.1038/emi.2016.69 PMID: 27436362

23. Rith S, Davis CT, Duong V, Sar B, Horm SV, Chin S, et al. Identification of Molecular Markers Associ- ated with Alteration of Receptor-Binding Specificity in a Novel Genotype of Highly Pathogenic Avian Influenza A(H5N1) Viruses Detected in Cambodia in 2013. J Virol 2014; 88:13897-909. https://doi. org/10.1128/JVI.01887-14 PMID: 25210193
24. Ly S, Horwood P, Chan M, Rith S, Sorn S, Oeung K, et al. Seroprevalence and Transmission of Human Influenza A(H5N1) Virus before and after Virus Reassortment, Cambodia, 2006–2014. Emerg Infect Dis 2017; 23:300–3. https://doi.org/10.3201/eid2302.161232 PMID: 28098551

25. Horwood PF, Horm SV, Suttie A, Thet S, Y P, Rith S, et al. Co-circulation of Influenza A H5, H7, and H9 Viruses and Co-infected Poultry in Live Bird Markets, Cambodia. Emerg Infect Dis 2018; 24:352–5. https://doi.org/10.3201/eid2402.171360 PMID: 29350140

26. QGIS Development Team. QGIS Geographic Information System. Open Source Geospatial Foundation 2018.

27. Zhou B, Donnelly ME, Scholes DT, St. George K, Hatta M, Kawaoka Y, et al. Single-Reaction Genomic Amplification Accelerates Sequencing and Vaccine Production for Classical and Swine Origin Human Influenza A Viruses. J Virol 2009; 83:10309–13. https://doi.org/10.1128/JVI.01109-09 PMID: 19605485

28. Shu Y, McCauley J. GISAID: Global initiative on sharing all influenza data-from vision to reality. Euro-surveillance 2017; 22.

29. Katoh K, Standley DM. MAFFT Multiple Sequence Alignment Software Version 7: Improvements in Performance and Usability. Mol Biol Evol 2013; 30:772–80. https://doi.org/10.1093/molbev/mst010 PMID: 23329690

30. Kosakovsky Pond SL, Posada D, Gravenor MB, Woelk CH, Frost SDW. Automated Phylogenetic Detection of Recombination Using a Genetic Algorithm. Mol Biol Evol 2006; 23:1891–901. https://doi.org/10.1093/molbev/msl051 PMID: 16818476

31. Weaver S, Shank SD, Spielman SJ, Li M, Muse SV, Kosakovsky Pond SL. Datamonkey 2.0: a modern web application for characterizing selective and other evolutionary processes. Mol Biol Evol 2018; 35:773–777. https://doi.org/10.1093/molbev/msx335 PMID: 29301006

32. Nguyen L-T, Schmidt HA, von Haeseler A, Minh BQ. IQ-TREE: a fast and effective stochastic algorithm for estimating maximum-likelihood phylogenies. Mol Biol Evol 2015; 32:268–274. https://doi.org/10.1093/molbev/msu300 PMID: 25371430

33. Minh BQ, Nguyen MAT, von Haeseler A. Ultrafast Approximation for Phylogenetic Bootstrap. Mol Biol Evol 2013; 30:1188–95. https://doi.org/10.1093/molbev/mst024 PMID: 23418397

34. Tavaré S. Some probabilistic and statistical problems in the analysis of DNA sequences. Lect Math Life Sci 1986; 17:57–86.

35. Andrew Rambaut. Molecular evolution, phylogenetics and epidemiology: FigTree 2016. http://tree.bio.ed.ac.uk/software/figtree/.

36. Drummond AJ, Suchard MA, Xie D, Rambaut A. Bayesian Phylogenetics with BEAUti and the BEAST 1.7. Mol Biol Evol 2012; 29:1969–73. https://doi.org/10.1093/molbev/mss075 PMID: 22367748

37. Miller MA, Pfeiffer W, Schwartz T. Creating the CIPRES Science Gateway for inference of large phylogenetic trees. 2010 Gatew. Comput. Environ. Workshop GCE, Ieee; 2010, p. 1–8.

38. Drummond AJ, Ho SY, Phillips MJ, Rambaut A. Relaxed phylogenetics and dating with confidence. PLoS Biol 2006; 4:e88. https://doi.org/10.1371/journal.pbio.0040088 PMID: 16683862

39. Shapiro B, Rambaut A, Drummond AJ. Choosing Appropriate Substitution Models for the Phylogenetic Analysis of Protein-Coding Sequences. Mol Biol Evol 2006; 23:7–9. https://doi.org/10.1093/molbev/msj021 PMID: 16177232

40. Minin VN, Bloomquist EW, Suchard MA. Smooth Skyride through a Rough Skyline: Bayesian Coalescent-Based Inference of Population Dynamics. Mol Biol Evol 2008; 25:1459–71. https://doi.org/10.1093/molbev/msn090 PMID: 18408232

41. Creanga A, Thi Nguyen D, Gerloff N, Thi Do H, Balish A, Dang Nguyen H, et al. Emergence of multiple clade 2.3.2.1 influenza A (H5N1) virus subgroups in Vietnam and detection of novel reassortants. Virology 2013; 444:12–20. https://doi.org/10.1016/j.virol.2013.06.005 PMID: 23849789

42. Nguyen DT, Jang Y, Nguyen TD, Jones J, Shepard SS, Yang H, et al. Shifting Clade Distribution, Reassortment, and Emergence of New Subtypes of Highly Pathogenic Avian Influenza A(H5) Viruses Collected from Vietnamese Poultry from 2012 to 2015. J Virol 2017; 91:e01708–16. https://doi.org/10.1128/JVI.01708-16.

43. Nguyen T, Rivailler P, Davis CT, Thi Hoa D, Balish A, Hoang Dang N, et al. Evolution of highly pathogenic avian influenza (H5N1) virus populations in Vietnam between 2007 and 2010. Virology 2012; 432:405–16. https://doi.org/10.1016/j.virol.2012.06.021 PMID: 22818871

44. Centers for Disease Control and Prevention. H5N1 Genetic Changes Inventory: A Tool for Influenza Surveillance and Preparedness. 2012.

45. Suttie A, Deng Y-M, Greenhill AR, Dussart P, Horwood PF, Karlsson EA. Inventory of molecular markers affecting biological characteristics of avian influenza A viruses. Virus Genes 2019. https://doi.org/10.1007/s11262-019-01700-z.
46. Gupta R, Jung E, Brunak S. Prediction of N-glycosylation sites in human proteins. 2004. Ref Type Unpubl Work 2016.

47. Burke DF, Smith DJ. A Recommended Numbering Scheme for Influenza A HA Subtypes. PLoS ONE 2014; 9:e112302. https://doi.org/10.1371/journal.pone.0112302 PMID: 25391151

48. Pond SLK, Muse SV. HyPhy: hypothesis testing using phylogenies. Stat. Methods Mol. Evol., Springer; 2005, p. 125–181.

49. Murrell B, Moola S, Mabona A, Weighill T, Sheward D, Kosakovsky Pond SL, et al. FUBAR: a fast, unconstrained bayesian approximation for inferring selection. Mol Biol Evol 2013; 30:1196–1205. https://doi.org/10.1093/molbev/mst030 PMID: 23420840

50. Murrell B, Wertheim JO, Moola S, Weighill T, Scheffler K, Pond SLK. Detecting individual sites subject to episodic diversifying selection. PLoS Genet 2012; 8:e1002764. https://doi.org/10.1371/journal.pgen.1002764 PMID: 22807683

51. Kosakovsky Pond SL, Frost SD. Not so different after all: a comparison of methods for detecting amino acid sites under selection. Mol Biol Evol 2005; 22:1208–1222. https://doi.org/10.1093/molbev/msi105 PMID: 15703242

52. Hurt AC, Barr IG, Hartel G, Hampson AW. Susceptibility of human influenza viruses from Australasia and South East Asia to the neuraminidase inhibitors zanamivir and oseltamivir. Antiviral Res 2004; 62:37–45. https://doi.org/10.1016/j.antiviral.2003.11.008 PMID: 15026200

53. World Health Organisation and others. Meetings of the WHO working group on surveillance of influenza antiviral susceptibility—Geneva, November 2011 and June 2012. Wkly Epidemiol Rec 2012; 87:369–374. PMID: 23061103

54. Horwood PF, Horm SV, Suttie A, Thet S, Phalla Y, Rith S, et al. Co-circulation of Influenza A H5, H7, and H9 Viruses and Co-infected Poultry in Live Bird Markets, Cambodia. Emerg Infect Dis 2018; 24. https://doi.org/10.3201/eid2402.171360.

55. Peng W, Bouwman KM, McBride R, Grant OC, Woods RJ, Verheije MH, et al. Enhanced Human-Type Receptor Binding by Ferret-Transmissible H5N1 with a K193T Mutation. J Virol 2018; 92. https://doi.org/10.1128/JVI.02016-17.

56. Chutinimitkul S, van Riel D, Munster VJ, van den Brand JMA, Rimmelzwaan GF, Kuiken T, et al. In Vitro Assessment of Attachment Pattern and Replication Efficiency of H5N1 Influenza A Viruses with Altered Receptor Specificity. J Virol 2010; 84:6825–33. https://doi.org/10.1128/JVI.02737-09PMID: 20392847

57. Li J, Li Y, Hu Y, Chang G, Sun W, Yang Y, et al. PB1-mediated virulence attenuation of H5N1 influenza virus in mice is associated with PB2. J Gen Virol 2011; 92:1435–44. https://doi.org/10.1099/vir.0.030718-0 PMID: 21367983

58. Taft AS, Ozawa M, Fitch A, Depasse JV, Halfmann PJ, Hill-Batorski L, et al. Identification of mammalian-adapting mutations in the polymerase complex of an avian H5N1 influenza virus. Nat Commun 2015; 6. https://doi.org/10.1038/ncomms8491.

59. Feng X, Wang Z, Shi J, Deng G, Kong H, Tao S, et al. Glycine at Position 622 in PB1 Contributes to the Virulence of H5N1 Avian Influenza Virus in Mice. J Virol 2016; 90:1872–9. https://doi.org/10.1128/JVI.02387-15 PMID: 26656683

60. Schmolke M, Manicassamy B, Pena L, Sutton T, Hai R, Varga ZT, et al. Differential Contribution of PB1-F2 to the Virulence of Highly Pathogenic H5N1 Influenza A Virus in Mammalian and Avian Species. PLoS Pathog 2011; 7:e1002186. https://doi.org/10.1371/journal.ppat.1002186 PMID: 21852950

61. Kamal RP, Kumar A, Davis CT, Tzeng W-P, Nguyen T, Donis RO, et al. Emergence of Highly Pathogenic Avian Influenza A(H5N1) Virus PB1-F2 Variants and Their Virulence in BALB/c Mice. J Virol 2015; 89:5835–46. https://doi.org/10.1128/JVI.03137-14 PMID: 25787281

62. Hu J, Mo Y, Wang X, Gu M, Hu Z, Zhong L, et al. PA-X Decreases the Pathogenicity of Highly Pathogenic H5N1 Influenza A Virus in Avian Species by Inhibiting Virus Replication and Host Response. J Virol 2015; 89:4126–42. https://doi.org/10.1128/JVI.02132-14 PMID: 25631083

63. Gao H, Sun Y, Hu J, Qi L, Wang J, Xiong X, et al. The contribution of PA-X to the virulence of pandemic 2009 H1N1 and highly pathogenic H5N1 avian influenza viruses. Sci Rep 2015; 5. https://doi.org/10.1038/srep08262.

64. Gao H, Sun Y, Liu X, Sun H, Hu J, Wang J, et al. Twenty amino acids at the C-terminus of PA-X are associated with increased influenza A virus replication and pathogenicity. J Gen Virol 2015; 96: 2036–49. https://doi.org/10.1099/vir.0.000143 PMID: 25877935

65. OFFLU. Influenza A Cleavage Sites 2019.

66. Su Y, Yang H-Y, Zhang B-J, Jia H-L, Tien P. Analysis of a point mutation in H5N1 avian influenza virus hemagglutinin in relation to virus entry into live mammalian cells. Arch Virol 2008; 153:2253–61. https://doi.org/10.1007/s00705-008-0255-y PMID: 19020946
67. Yang Z-Y, Wei C-J, Kong W-P, Wu L, Xu L, Smith DF, et al. Immunization by Avian H5 Influenza Hemagglutinin Mutants with Altered Receptor Binding Specificity. Science 2007; 317:825-8. https://doi.org/10.1126/science.1135165 PMID: 17690300

68. Wang W, Lu B, Zhou H, Suguitan AL, Cheng X, Subbarao K, et al. Glycosylation at 158N of the Hemagglutinin Protein and Receptor Binding Specificity Synergistically Affect the Antigenicity and Immuno-ogenicity of a Live Attenuated H5N1 A/Vietnam/1203/2004 Vaccine Virus in Ferrets. J Virol 2010; 84:6570-7. https://doi.org/10.1128/JVI.00221-10 PMID: 20427525

69. Tada T, Suzuki K, Sakurai Y, Kubo M, Okada H, Itoh T, et al. NP body domain and PB2 contribute to increased virulence of H5N1 highly pathogenic avian influenza viruses in chickens. J Virol 2011; 85:1834-46. https://doi.org/10.1128/JVI.01648-10 PMID: 21123376

70. Wasilenko JL, Sarmento L, Pantin-Jackwood MJ. A single substitution in amino acid 184 of the NP protein alters the replication and pathogenicity of H5N1 avian influenza viruses in chickens. Arch Virol 2009; 154:969-79. https://doi.org/10.1007/s00705-009-0399-4 PMID: 19475480

71. Matsuoka Y, Swayne DE, Thomas C, Rameix-Welti M-A, Naffakh N, Warnes C, et al. Neuraminidase Stalk Length and Additional Glycosylation of the Hemagglutinin Influence the Virulence of Influenza H5N1 Viruses for Mice. J Virol 2009; 83:4704-8. https://doi.org/10.1128/JVI.01987-08 PMID: 19225004

72. Fan S, Deng G, Song J, Tian G, Suo Y, Jiang Y, et al. Two amino acid residues in the matrix protein M1 contribute to the virulence difference of H5N1 avian influenza viruses in mice. Virology 2009; 384:28-32. https://doi.org/10.1016/j.virol.2008.11.044 PMID: 19117585

73. Nao N, Kajihara M, Manzoor R, Maruyama J, Yoshida R, Muramatsu M, et al. A Single Amino Acid in the M1 Protein Responsible for the Different Pathogenic Potentials of H5N1 Highly Pathogenic Avian Influenza Virus Strains. PloS One 2015; 10:e0137989. https://doi.org/10.1371/journal.pone.0137989 PMID: 26368015

74. Cheung C-L, Rayner JM, Smith GJ, Wang P, Naipospos TSP, Zhang J, et al. Distribution of amantadine-resistant H5N1 avian influenza variants in Asia. J Infect Dis 2006; 193:1626-1629. https://doi.org/10.1086/504723 PMID: 16703504

75. Bean WJ, Threlkeld SC, Webster RG. Biologic potential of amantadine-resistant influenza A virus in an avian model. J Infect Dis 1989; 159:1050-6. https://doi.org/10.1093/infdis/159.6.1050 PMID: 2723453

76. Long J-X, Peng D-X, Liu Y-L, Wu Y-T, Liu X-F. Virulence of H5N1 avian influenza virus enhanced by a 15-nucleotide deletion in the viral nonstructural gene. Virus Genes 2008; 36:471-8. https://doi.org/10.1007/s11262-007-0187-8 PMID: 18317917

77. Jiao P, Tian G, Li Y, Deng G, Jiang Y, Liu C, et al. A Single-Amino-Acid Substitution in the NS1 Protein Changes the Pathogenicity of H5N1 Avian Influenza Viruses in Mice. J Virol 2008; 82:1146-54. https://doi.org/10.1128/JVI.01698-07 PMID: 18032512

78. Kuo R-L, Krug RM. Influenza A Virus Polymerase Is an Integral Component of the CPSF30-NS1A Protein Complex in Infected Cells. J Virol 2009; 83:1611-6. https://doi.org/10.1128/JVI.01491-08 PMID: 19052083

79. Heui Seo S, Hoffmann E, Webster RG. Lethal H5N1 influenza viruses escape host anti-viral cytokine responses. Nat Med 2002; 8:950-4. https://doi.org/10.1038/nm757 PMID: 12195436

80. Spesock A, Malur M, Hossain MJ, Chen L-M, Njaa BL, Davis CT, et al. The Virulence of 1997 H5N1 Influenza Viruses in the Mouse Model Is Increased by Correcting a Defect in Their NS1 Proteins. J Virol 2011; 85:7048-58. https://doi.org/10.1128/JVI.00417-11 PMID: 21593152

81. Obayashi E, Yoshida H, Kawai F, Shibayama N, Kawaguchi A, Nagata K, et al. The structural basis for an essential subunit interaction in influenza virus RNA polymerase. Nature 2008; 454:1127-31. https://doi.org/10.1038/nature07225 PMID: 18660801

82. Sugiyama K, Obayashi E, Kawaguchi A, Suzuki Y, Tame JRH, Nagata K, et al. Structural insight into the essential PB1-PB2 subunit contact of the influenza virus RNA polymerase. EMBO J 2009; 28:1803-11. https://doi.org/10.1038/emboj.2009.138 PMID: 19461581

83. Guilligay D, Tarendeau F, Resa-Infante P, Coloma R, Crepin T, Sehr P, et al. The structural basis for cap binding by influenza virus polymerase subunit PB2. Nat Struct Mol Biol 2008; 15:500-6. https://doi.org/10.1038/nsmb.1421 PMID: 18454157

84. Fan S, Hatta M, Kim JH, Halfmann P, Imai M, Macken CA, et al. Novel residues in avian influenza virus PB2 protein affect virulence in mammalian hosts. Nat Commun 2014; 5. https://doi.org/10.1038/ncomms6021.

85. Liu Y, Qin K, Meng G, Zhang J, Zhou J, Zhao G, et al. Structural and functional characterization of K339T substitution identified in the PB2 subunit cap-binding pocket of influenza A virus. J Biol Chem 2013; 288:11013-23. https://doi.org/10.1074/jbc.M112.392878 PMID: 23436652
86. Yamaji R, Yamada S, Le MQ, Li C, Chen H, Qurnianingsih E, et al. Identification of PB2 mutations responsible for the efficient replication of H5N1 influenza viruses in human lung epithelial cells. J Virol 2015; 89:3947-56. https://doi.org/10.1128/JVI.03328-14 PMID: 25609813

87. Chen W, Zhong Y, Qin Y, Sun S, Li Z. The Evolutionary Pattern of Glycosylation Sites in Influenza Virus (H5N1) Hemagglutinin and Neuraminidase. PLoS ONE 2012; 7:e49224. https://doi.org/10.1371/journal.pone.0049224 PMID: 23133677

88. Suttie A, Karlsson EA, Deng Y-M, Horm SV, Yann S, Tok S, et al. Influenza A(H5N1) viruses with A (H9N2) single gene (matrix or PB1) reassortment isolated from Cambodian live bird markets. Virology 2018; 523:22-6. https://doi.org/10.1016/j.virol.2018.07.028 PMID: 30075357

89. Naughtin M, Dyason JC, Mardy S, Sorn S, von Itzstein M, Buchy P. Neuraminidase Inhibitor Sensitivity and Receptor-Binding Specificity of Cambodian Clade 1 Highly Pathogenic H5N1 Influenza Virus. Antimicrob Agents Chemother 2011; 55:2004-10. https://doi.org/10.1128/AAC.01773-10 PMID: 21343450

90. FAO. Investigation of duck production and hatcheries and duckling supply in Cambodia. AHBL—Promoting strategies for prevention and control of HPAI. Rome: 2009.

91. Horm SV, Sorn S, Allal L, Buchy P. Influenza A(H5N1) Virus Surveillance at Live Poultry Markets, Cambodia, 2011. Emerg Infect Dis 2013; 19:305-8. https://doi.org/10.3201/eid1902.121201 PMID: 23347451

92. World Health Organization (WHO/OIE/FAO) H5N1 Evolution Working Group. Revised and updated nomenclature for highly pathogenic avian influenza A (H5N1) viruses. Influenza Other Respir Viruses 2014; 8:384-8. https://doi.org/10.1111/irv.12230 PMID: 24483237

93. Naguib MM, Kinne J, Chen H, Chan K-H, Joseph S, Wong P-C, et al. Outbreaks of highly pathogenic avian influenza H5N1 clade 2.3. 2.1 c in hunting falcons and kept wild birds in Dubai implicate intercontinal virus spread. J Gen Virol 2015; 96:3212-3222. https://doi.org/10.1099/jgv.0.000274 PMID: 26350163

94. Marchenko VY, Susloparov IM, Kolosova NP, Goncharova NI, Shipovalov AV, Ilyicheva TN, et al. Highly pathogenic influenza H5N1 virus of clade 2.3.2.1c in Western Siberia. Arch Virol 2016; 161:1645-9. https://doi.org/10.1007/s00705-016-2800-4 PMID: 26935914

95. Monamele CG, Y P, Karlsson EA, Vernet M-A, Wade A, Okomo M-CA, et al. Evidence of exposure and human seroconversion during an outbreak of avian influenza A(H5N1) among poultry in Cameroon. Emerg Microbes Infect 2019; 8:186-96. https://doi.org/10.1080/22221751.2018.1564631 PMID: 30866772

96. Bi Y, Chen J, Zhang Z, Li M, Cai T, Sharshov K, et al. Highly pathogenic avian influenza H5N1 Clade 2.3.2.1c virus in migratory birds, 2014-2015. Virol Sin 2016; 31:300-5. https://doi.org/10.1007/s12250-016-3750-4 PMID: 27405930

97. He S, Shi J, Qi X, Huang G, Chen H, Lu C. Lethal infection by a novel reassortant H5N1 avian influenza A virus in a zoo-housed tiger. Microbes Infect 2015; 17:54-61. https://doi.org/10.1016/j.micinf.2014.10.004 PMID: 25461468

98. Chen Q, Wang H, Zhao L, Ma L, Wang R, Lei Y, et al. First documented case of avian influenza (H5N1) virus infection in a lion. Emerg Microbes Infect 2016; 5:e125-e125. https://doi.org/10.1038/emi.2016.127 PMID: 27999425

99. Buchy P, Fourment M, Mardy S, Sorn S, Holl D, Ly S, et al. Molecular Epidemiology of Clade 1 Influenza A Viruses (H5N1), Southern Indochina Peninsula, 2004-2007. Emerg Infect Dis 2009; 15: 1641-4. https://doi.org/10.3201/eid1510.090115 PMID: 19861062

100. Van Kerkhove MD, Vong S, Guitian J, Holl D, Mangtani P, San S, et al. Poultry movement networks in Cambodia: Implications for surveillance and control of highly pathogenic avian influenza (HPAI/H5N1). Vaccine 2009; 27:6345-52. https://doi.org/10.1016/j.vaccine.2009.05.004 PMID: 19840671

101. Chen H, Yuan H, Gao R, Zhang J, Wang D, Xiong Y, et al. Clinical and epidemiological characteristics of a fatal case of avian influenza A H10N8 virus infection: a descriptive study. The Lancet 2014; 383:714-721.

102. Pu J, Wang S, Yin Y, Zhang G, Carter RA, Wang J, et al. Evolution of the H9N2 influenza genotype that facilitated the genesis of the novel H7N9 virus. Proc Natl Acad Sci 2015; 112:548-53. https://doi.org/10.1073/pnas.1422456112 PMID: 25548189

103. Hatta M. Molecular Basis for High Virulence of Hong Kong H5N1 Influenza A Viruses. Science 2001; 293:1840-2. https://doi.org/10.1126/science.1062882 PMID: 11546875

104. Li Z, Chen H, Jiao P, Deng G, Tian G, Li Y, et al. Molecular Basis of Replication of Duck H5N1 Influenza Viruses in a Mammalian Mouse Model. J Virol 2005; 79:12058-64. https://doi.org/10.1128/JVI.79.18.12058-12064.2005 PMID: 16140781
105. Skehel JJ, Stevens DJ, Daniels RS, Douglas AR, Knossow M, Wilson IA, et al. A carbohydrate side chain on hemagglutinins of Hong Kong influenza viruses inhibits recognition by a monoclonal antibody. Proc Natl Acad Sci 1984; 81:1779–83. https://doi.org/10.1073/pnas.81.6.1779 PMID: 6584912

106. Wang C-C, Chen J-R, Tseng Y-C, Hsu C-H, Hung Y-F, Chen S-W, et al. Glycans on influenza hemagglutinin affect receptor binding and immune response. Proc Natl Acad Sci 2009; 106:18137–18142. https://doi.org/10.1073/pnas.0909696106 PMID: 19822741

107. Hervé P-L, Lorin V, Jouvion G, Da Costa B, Escriou N. Addition of N-glycosylation sites on the globular head of the H5 hemagglutinin induces the escape of highly pathogenic avian influenza A H5N1 viruses from vaccine-induced immunity. Virology 2015; 486:134–45. https://doi.org/10.1016/j.virol.2015.08.033 PMID: 26433051

108. Deshpande KL, Fried VA, Ando M, Webster RG. Glycosylation affects cleavage of an H5N2 influenza virus hemagglutinin and regulates virulence. Proc Natl Acad Sci 1987; 84:36–40. https://doi.org/10.1073/pnas.84.1.36 PMID: 3467357

109. Kawaoka Y, Webster RG. Interplay between Carbohydrate in the Stalk and the Length of the Connecting Peptide Determines the Cleavability of Influenza Virus Hemagglutinin n.d.:5.

110. Gao Y, Zhang Y, Shinya K, Deng G, Jiang Y, Li Z, et al. Identification of Amino Acids in HA and PB2 Critical for the Transmission of H5N1 Avian Influenza Viruses in a Mammalian Host. PLoS Pathog 2009; 5:e1000709. https://doi.org/10.1371/journal.ppat.1000709 PMID: 20041223

111. Demicheli V, Jefferson T, Ferroni E, Rivetti A, Di Pietrantonj C. Vaccines for preventing influenza in healthy adults. Cochrane Database Syst Rev 2018.

112. Bright RA, Medina M, Xu X, Perez-Oronoz G, Wallis TR, Davis XM, et al. Incidence of adamantane resistance among influenza A (H3N2) viruses isolated worldwide from 1994 to 2005: a cause for concern. The Lancet 2005; 366:1175–1181.