the nomenclature system designed by the who / oie / fao h5n1 evolution working group to classify the a / goose / guangdong/1/1996 lineage of eurasian highly pathogenic h5n1 avian influenza viruses has been in place since early 2008 . based on the criteria used to distinguish variant groups of the h5 hemagglutinin ( ha ) gene ,
the system has formally identified 20 distinct clades of the virus . 
 , 
  using phylogenetic characterization and nucleotide sequence divergence , h5n1 virus clades have been defined based on sharing of a common ancestral node and monophyletic evolution with a bootstrap value of 60 at the cladedefining node ( after 1000 neighborjoining bootstrap replicates ) and maintaining the average percentage pairwise nucleotide distances between and within clades of > 15% and < 15% , respectively . 
  since the inception of this approach to h5n1 classification , clade 2.2 viruses ( first detected in the qinghai lake area ) introduced into egypt diverged into a thirdorder clade designated 2.2.1 . 
 
as expected , however , the continued evolution of h5n1 has led to the periodic emergence of new phylogenetic groups of h5n1 in several regions of the world .
therefore , the who / oie / fao h5n1 evolution working group has examined the ha sequence data available as of january 2011 in the context of previous nomenclature assignments to update this classification system .
hemagglutinin sequences of the highly pathogenic avian influenza a ( h5n1 ) viruses within the eurasian a / goose / guangdong/1/1996 lineage were obtained from the epiflu database in gisaid , which contains all of the h5n1 ha sequences available in the influenza virus resource of ncbi , as well as shared data from individual laboratories within the who global influenza surveillance and response system ( gisrs ) and oie / fao global network of expertise on animal influenza ( offlu ) . only sequences with at least 1600 nucleotides of the ha gene were included ( > 90% of the fulllength open reading frame ) . a nucleotide alignment containing 2947 ha sequences
a large neighborjoining tree was initially constructed from the alignment using a gtr+i+g model in paup * v4.0b10 , 
  and clade assignments ( using previously established nomenclature parameters ) were made for new sequences based on their location within a clade ( figure  s1 ) .
cladespecific alignments were then used to generate neighborjoining trees with 1000 bootstraps and rooted to either a / goose / guangdong/1/96 or a common outlier to that clade ( figure  s2a  g ) .
monophyletic groups within each of the clades defined in the 2009 update 
  containing more than four nonredundant sequences with bootstrap values of 60 were selected as potential candidates for new order clade designation . to determine whether a clade had undergone a sufficient divergence to exceed the preestablished limit ( i.e. required splitting ) , the average withingroup nucleotide pdistance for sequences was determined using mega4 
 ; clades having withingroup distances 15% were analyzed further .
only groups containing samples with a collection date of 2007 or later were subjected to this analysis .
new order clade designations were made for those groups that had betweengroup divergences of 15% and withingroup divergences of 15% . as described previously , the new clade number consisted of the name of the originating clade followed by a period and a single numeric digit added on the right ( figure  s2a  g ) . following the identification of new higherorder clades , neighborjoining and bayesian trees of 196 representative sequences
were constructed using paup * v4.0b10 
  and mrbayes , 
  to verify the overall h5n1 ha tree topology and for statistical inferences ( figure  1 ) . 
phylogenetic relationships of recently diverged a / goose / guangdong/1/1996like h5 hemagglutinin ( ha ) genes .
neighborjoining tree of 196 ha gene sequences from h5n1 viruses constructed with 1000 bootstrap replicates ( above branches ) and bayesian posterior probabilities ( below branches ) .
the tree was constructed using paup * v4.0b10 
  with a gtr+i+g model and rooted to a / goose / guangdong/1/1996 .
solid triangles denote ha clades of viruses that have not been in circulation since 2008 or earlier ; the year of the most recent isolation is shown in parenthesis .
the sequence data set used to generate this figure is provided as supporting information in fasta format ( supplementary data s2 ) .
hemagglutinin sequences of the highly pathogenic avian influenza a ( h5n1 ) viruses within the eurasian a / goose / guangdong/1/1996 lineage were obtained from the epiflu database in gisaid , which contains all of the h5n1 ha sequences available in the influenza virus resource of ncbi , as well as shared data from individual laboratories within the who global influenza surveillance and response system ( gisrs ) and oie / fao global network of expertise on animal influenza ( offlu ) . only sequences with at least 1600 nucleotides of the ha gene were included ( > 90% of the fulllength open reading frame ) . a nucleotide alignment containing 2947 ha sequences
a large neighborjoining tree was initially constructed from the alignment using a gtr+i+g model in paup * v4.0b10 , 
  and clade assignments ( using previously established nomenclature parameters ) were made for new sequences based on their location within a clade ( figure  s1 ) .
cladespecific alignments were then used to generate neighborjoining trees with 1000 bootstraps and rooted to either a / goose / guangdong/1/96 or a common outlier to that clade ( figure  s2a  g ) .
monophyletic groups within each of the clades defined in the 2009 update 
  containing more than four nonredundant sequences with bootstrap values of 60 were selected as potential candidates for new order clade designation .
to determine whether a clade had undergone a sufficient divergence to exceed the preestablished limit ( i.e. required splitting ) , the average withingroup nucleotide pdistance for sequences was determined using mega4 
 ; clades having withingroup distances 15% were analyzed further .
only groups containing samples with a collection date of 2007 or later were subjected to this analysis .
new order clade designations were made for those groups that had betweengroup divergences of 15% and withingroup divergences of 15% . as described previously , the new clade number consisted of the name of the originating clade followed by a period and a single numeric digit added on the right ( figure  s2a  g ) . following the identification of new higherorder clades , neighborjoining and bayesian trees of 196 representative sequences
were constructed using paup * v4.0b10 
  and mrbayes , 
  to verify the overall h5n1 ha tree topology and for statistical inferences ( figure  1 ) . 
phylogenetic relationships of recently diverged a / goose / guangdong/1/1996like h5 hemagglutinin ( ha ) genes .
neighborjoining tree of 196 ha gene sequences from h5n1 viruses constructed with 1000 bootstrap replicates ( above branches ) and bayesian posterior probabilities ( below branches ) .
the tree was constructed using paup * v4.0b10 
  with a gtr+i+g model and rooted to a / goose / guangdong/1/1996 .
solid triangles denote ha clades of viruses that have not been in circulation since 2008 or earlier ; the year of the most recent isolation is shown in parenthesis .
the sequence data set used to generate this figure is provided as supporting information in fasta format ( supplementary data s2 ) .
surveillance for the continuing circulation of highly pathogenic avian influenza a ( h5n1 ) in poultry and wild birds in asia , the middle east , europe , and africa has led to an increase in the number of sequences available for analysis .
the current analysis was applied to 2947 ha gene sequences consisting of at least 1600 nucleotides , of which 1637 sequences ( primarily from 2008 to 2010 ) were added since the previous nomenclature update .
phylogenetic analysis of the neighborjoining tree generated from this data set revealed that the addition of these new sequence data resulted in the generation of one or more monophyletic groups with high bootstrap support within each of the currently circulating clades of h5n1 ( figure  s1 ) .
in addition , many of these groups had long branches separating them from the nearest node in the tree , indicating further nucleotide divergence .
each of the currently circulating clades ( 1 , 2.1.3 , 2.2 , 2.2.1 , 2.3.2 , 2.3.4 , and 7 ) was then examined independently , to measure the average withingroup pairwise nucleotide distances .
all were found to have > 15% withingroup divergence indicating the need to split these groups into new order clades .
after generation of cladespecific trees , monophyletic groups with bootstrap values more than 60 were selected for further testing ( figure  s2a  g ) .
based on the previously defined nomenclature , 12 new second , third , and fourthorder clades were identified ( table  1 ) .  updates to the h5n1 highly pathogenic avian influenza virus clade nomenclature system to summarize these findings and to confirm the consistency of clade topology using a smaller data set , a small tree containing 196 ha sequences was generated using both nj and bayesian methods ( figure  1 ) .
this tree is annotated with the 12 new clades presented here and displays older , noncirculating , and putatively noncirculating clades as collapsed nodes ( figure  1 ) .
both nj bootstrap support and bayesian posterior probabilities were calculated to demonstrate confidence in clade assignments when using smaller data sets .
surveillance for the continuing circulation of highly pathogenic avian influenza a ( h5n1 ) in poultry and wild birds in asia , the middle east , europe , and africa has led to an increase in the number of sequences available for analysis .
the current analysis was applied to 2947 ha gene sequences consisting of at least 1600 nucleotides , of which 1637 sequences ( primarily from 2008 to 2010 ) were added since the previous nomenclature update .
phylogenetic analysis of the neighborjoining tree generated from this data set revealed that the addition of these new sequence data resulted in the generation of one or more monophyletic groups with high bootstrap support within each of the currently circulating clades of h5n1 ( figure  s1 ) .
in addition , many of these groups had long branches separating them from the nearest node in the tree , indicating further nucleotide divergence .
each of the currently circulating clades ( 1 , 2.1.3 , 2.2 , 2.2.1 , 2.3.2 , 2.3.4 , and 7 ) was then examined independently , to measure the average withingroup pairwise nucleotide distances .
all were found to have > 15% withingroup divergence indicating the need to split these groups into new order clades .
after generation of cladespecific trees , monophyletic groups with bootstrap values more than 60 were selected for further testing ( figure  s2a  g ) .
based on the previously defined nomenclature , 12 new second , third , and fourthorder clades were identified ( table  1 ) .  updates to the h5n1 highly pathogenic avian influenza virus clade nomenclature system to summarize these findings and to confirm the consistency of clade topology using a smaller data set , a small tree containing 196 ha sequences was generated using both nj and bayesian methods ( figure  1 ) .
this tree is annotated with the 12 new clades presented here and displays older , noncirculating , and putatively noncirculating clades as collapsed nodes ( figure  1 ) .
both nj bootstrap support and bayesian posterior probabilities were calculated to demonstrate confidence in clade assignments when using smaller data sets .
the highly pathogenic avian influenza virus a ( h5n1 ) lineage was first identified in 1996 and has subsequently spread over much of the eastern hemisphere .
as these viruses have established geographically and ecologically isolated enzootic foci of infection , new clades of viruses continuously evolve .
the who / oie / fao h5n1 evolution working group has previously identified 20 such clades and established specific criteria for naming h5n1 clades . since the most recent recommendations of the h5n1 evolution working group in 2009
, h5n1 viruses have continued to spread and evolve locally and regionally . to address the question of whether this continued evolution has led to the appearance of new clades , we initially constructed a phylogenetic tree composed of nearly 3000 ha sequences from either public databases or ongoing surveillance programs .
as expected , the addition of more than 1500 new sequences to the previous h5n1 nomenclature update tree produced a number of monophyletic groups within each of the currently circulating clades of h5n1 .
after measuring within and betweengroup average nucleotide pdistances , it was observed that all extant clades of h5n1 required splitting into one or more newly defined second , third , and/or fourthorder clades .
the current analysis , therefore , proposes the designation of 12 new clades that meet the criteria of the h5n1 nomenclature system .
the designation of 12 new h5n1 clades ( in contrast with the previous nomenclature update in which only one new clade was identified 
 ) is not surprising , considering the length of time since the last nomenclature update and improved surveillance / reporting of viruses in recent years , increasing the number of sequences available for analysis .
interestingly , the clades that had already diverged into thirdorder groups in previous analyses required the largest number of new clade designations .
previous assignment to thirdorder clades indicates that the viruses had already evolved considerably , presumably because of short generation times and large population sizes in enzootic areas .
the subdivision of these clades suggests the continued evolution and divergence in regions where h5n1 is considered entrenched in poultry . 
 
for example , clade 2.3.4 viruses circulating in china , vietnam , laos , and neighboring regions were split into 2.3.4.1/2.3.4.2/2.3.4.3 , whereas clade 2.1.3 in indonesia was split into 2.1.3.1/2.1.3.2/2.1.3.3 .
the egyptian thirdorder clade 2.2.1 also required splitting into new clade 2.2.1.1 , indicating further divergence of contemporary strains of h5n1 circulating in egyptian poultry . 
 
in addition , clade 2.3.2 enzootic in china was subdivided to accommodate the new clade 2.3.2.1 , which has been identified in poultry populations as well as wild birds in many diverse regions outside of china . 
 , 
 
while these groups have been identified as sufficiently divergent from the previously defined clades , their ancestral populations have not been detected in recent years and are likely to no longer circulate .
one new thirdorder clade , termed 2.2.2 , comprising viruses from india and bangladesh was identified .
this clade diverged from clade 2.2 after its introduction into this region in approximately 2006 . 
 
finally , firstorder clades 1 and 7 required designation of secondorder clades as a result of their continued divergence in the mekong delta regions of southern vietnam / cambodia and china , respectively .
while clade 1.1 appears to represent a continuum of evolution in the relatively wellsampled mekong delta region , 
 , 
  clade 7.1 and 7.2 viruses remain widely diverse groups ( average betweengroup divergence of 45% ) , suggesting that these groups represent undersampled virus populations with only sporadic detection ( figure  s2a  g ) . 
 , 
  antigenic analysis of viruses from each of these new clades will inform decisions on whether to update the library of prepandemic vaccine candidates ( see http://www.who.int/influenza/resources/documents/characteristics_virus_vaccines/en/index.html ) .
in addition to the emergence of novel clades of h5n1 , a number of previously circulating clades of h5n1 have not been detected for several years ( since at least 2008 or earlier ) .
although we can not rule out the possibility that the lack of detection of viruses in these clades may in some cases be related to gaps in surveillance , it is more likely that many of them have been supplanted by new clades and become inactive . as shown in figure  1 ,
13 clades ( shown as collapsed nodes in the tree with the last year of detection shown in parentheses ) have not been detected since 2008 or earlier ( 0 , 2.1.1 , 2.1.2 , 2.3.1 , 2.3.3 , 2.4 , 2.6 , 3 , 4 , 5 , 6 , 8 , and 9 ) .
therefore , although 12 new clades have been identified and characterized in this study , 13 older clades may now be considered out of circulation . again , this is not unexpected given the time since the last nomenclature update .
as observed previously , the majority of viruses cluster with others of similar geographic and temporal origin , with few exceptions that might indicate viral spread over longer distances ( i.e. , poultry trade and/or bird migration ) . 

, 
  for example , clade 2.3.2.1 viruses have now been detected in 9 countries since 2007 .
in addition , while the majority of the novel ha gene clades have approximately similar levels of within ( < 15% nucleotide differences ) and between ( > 15% ) group diversity , the clade 2.3.2.1 viruses have a higher withingroup mean nucleotide divergence ( 17% ) .
this result suggests that the clade 2.3.2.1 ha genes have diversified to a greater extent than those in other clades , perhaps owing to the adaptation to different bird species and introduction into multiple geographic regions by migratory birds .
although this average withingroup nucleotide pdistance exceeds the < 15% withingroup mean criterion originally proposed as a cutoff for clade determination , 
 , 
  designation of a single fifthorder clade awaits additional data to determine its significance . despite very successful efforts to eliminate or control the spread of h5n1 in poultry , viruses persisting in enzootic regions seed virus epizootics in new areas via poultry trade or wild birds .
periodic population sweeps in enzootic and epizootic areas with sustained viral circulation may lead to the extinction of some clades and extensive divergence of viruses within more successful clades . as a consequence ,
the classification of h5n1 viruses based on ha evolution requires periodic updating , making classification dynamic as the virus has expanded within several disparate ecosystems and along distinct evolutionary trajectories .
although this report focuses strictly on genetic divergence as a measure of h5n1 classification , clearly additional studies on h5n1 antigenicity and crossclade protective immunity are needed to inform prepandemic vaccine strain selection as well as risk management policies for both veterinary and human health . as always , continued surveillance for h5n1 in animals and humans is crucial for this process to be most effective .
data s2 . the nucleotide sequence alignment used to generate the tree shown in figure 1
( a  g ) cladespecific 
 neighborjoining trees with 1000 bootstrap replicates were 
 generated for each clade and rooted to either 
 agooseguangdong196 or a common outlier 
 to that clade .
 
 please cite this paper as : who / oie / fao . 
 ( 2012 ) continued evolution of highly pathogenic avian influenza a(h5n1 ) : updated nomenclature . influenza and other respiratory viruses 6(1 ) , 15 . 
 
 background   continued evolution of highly pathogenic avian influenza a ( h5n1 ) throughout many regions of the eastern hemisphere has led to the emergence of new phylogenetic groups . 
 a total of 1637 new h5n1 hemagglutinin ( ha ) sequences have become available since the previous nomenclature recommendations described in 2009 by the who / oie / fao h5n1 evolution working group . a comprehensive analysis including all the new data is needed to update ha clade nomenclature . 
 
 methods   phylogenetic trees were constructed from data sets of all available h5n1 ha sequences . 
 new clades were designated on the basis of phylogeny and pdistance using the preestablished nomenclature system ( emerg infec dis 2008 ; 14:e1 ) . 
 each circulating h5n1 clade was subjected to further phylogenetic analysis and nucleotide sequence divergence calculations . 
 
 results   all recently circulating clades ( clade 1 in the mekong river delta , 2.1.3 in indonesia , 2.2 in india / bangladesh , 2.2.1 in egypt , 2.3.2 , 2.3.4 and 7 in asia ) required assignment of divergent ha genes to new second , third , and/or fourthorder clades . at the same time , clades 0 , 3 , 4 , 5 , 6 , 8 , 9 , and several second and thirdorder groups from clade 2 have not been detected since 2008 or earlier . 
 
 conclusions   new designations are recommended for 12 ha clades , named according to previously defined criteria . 
 in addition , viruses from 13 clades have not been detected since 2008 or earlier . 
 the periodic updating of this dynamic classification system allows continued use of a unified nomenclature in all h5n1 studies . 
Introduction
Materials and methods
Data sets and phylogenetic analysis
Genetic distance calculations
Results
Further divergence of the H5N1 HA gene
Discussion
Supporting information