Molecular characterization of highly pathogenic H5N1 avian influenza viruses isolated in Sweden in 2006
© Kiss et al; licensee BioMed Central Ltd. 2008
Received: 22 August 2008
Accepted: 06 October 2008
Published: 06 October 2008
The analysis of the nonstructural (NS) gene of the highly pathogenic (HP) H5N1 avian influenza viruses (AIV) isolated in Sweden early 2006 indicated the co-circulation of two sub-lineages of these viruses at that time. In order to complete the information on their genetic features and relation to other HP H5N1 AIVs the seven additional genes of twelve Swedish isolates were amplified in full length, sequenced, and characterized.
The presence of two sub-lineages of HP H5N1 AIVs in Sweden in 2006 was further confirmed by the phylogenetic analysis of approximately the 95% of the genome of twelve isolates that were selected on the base of differences in geographic location, timing and animal species of origin. Ten of the analyzed viruses belonged to sub-clade 2.2.2. and grouped together with German and Danish isolates, while two 2.2.1. sub-clade viruses formed a cluster with isolates of Egyptian, Italian, Slovenian, and Nigerian origin. The revealed amino acid differences between the two sub-groups of Swedish viruses affected the predicted antigenicity of the surface glycoproteins, haemagglutinin and neuraminidase, rather than the nucleoprotein, polymerase basic protein 2, and polymerase acidic protein, the main targets of the cellular immune responses. The distinctive characteristics between members of the two subgroups were identified and described.
The comprehensive genetic characterization of HP H5N1 AIVs isolated in Sweden during the spring of 2006 is reported. Our data support previous findings on the coincidental spread of multiple sub-lineage H5N1 HPAIVs via migrating aquatic birds to large distance from their origin. The detection of 2.2.1. sub-clade viruses in Sweden adds further data regarding their spread in the North of Europe in 2006. The close genetic relationship of Swedish isolates sub-clade 2.2.2. to the contemporary German and Danish isolates supports the proposition of the introduction and spread of a single variant of 2.2.2. sub-clade H5N1 avian influenza viruses in the Baltic region. The presented findings underline the importance of whole genome analysis.
The first reports of outbreaks caused by highly pathogenic avian influenza viruses (HPAIV) of H5N1 subtype in 1996 originated from southern China . Systematic influenza surveillances showed that distinct genetic sub-lineages of H5N1 HPAIVs, reflecting on their geographic origin, have been established since then among domestic poultry and have been transmitted to long distances by migratory waterfowl [2, 3]. Europe experienced a peak of outbreaks of H5N1 HPAI in domestic poultry and wild birds in March 2006 – that was supposedly the consequence of an unusual westward movement of waterfowl from the Black Sea area [4–6]. The recent avian influenza virus strains of European-Middle Eastern-African (EMA) origin were assigned to three clades (EMA-1-3) based on the phylogeny of the complete genomes of the isolates , which are referred as sub-clades 2.2.1.-2.2.3. according to the more recent nomenclature . Further, clade 2.2. was classified into three sub-clades: Clade 2.2.1. appeared in Egypt, southern Germany, Italy, Mongolia, and some regions in sub-Sahara Africa. Clade 2.2.2. viruses were detected in northern Germany, Denmark, Sweden, Scotland, and Nigeria, while clade 2.2.3. viruses were demonstrated in India, Afghanistan, Italy, and Iran . Simultaneous transmission of different strains was reported in several European countries such as Sweden , Germany , France and Italy . Characterization of the Swedish H5N1 HPAIV isolates based on the nonstructural (NS) gene nucleotide sequences demonstrated that all belonged to clade 2.2. The majority of them clustered together with clade 2.2.2., viruses belonging to clade 2.2.1. were also introduced into Sweden .
The aim of this study was to further investigate the Swedish H5N1 HPAI viruses by sequencing twelve selected isolates representing four east-coast provinces of the area affected by the epidemic during March-April 2006. The sequence information was used to study the evolution and epidemiology of the outbreak of H5N1 in Europe during 2006. Further, a H5N1 strain isolated from a mink was investigated to reveal any possible adaptation towards mammals.
Results and discussion
According to the Influenza A Virus Genotype Tool  the studied genes of the investigated Swedish isolates belonged to the following lineages: PB2 (K), PB1 (G), PA (D), HA (5J), NP (F), NA (1J), MP (F), NS (1E).
HA amino acid residue 403 was observed to characterize 2.2.1. (isolates mainly from Southern parts of Germany) and 2.2.2. (German isolates from the North) sub-clade German H5N1 viruses because the former group contained mainly D while the latter N at this position. All sub-clade 2.2.2. Swedish H5N1 viruses possessed N at HA 403 position together with A/tufted duck/Sweden/599/06 sub-clade 2.2.1. isolate, and only 2.2.1. isolate A/herring gull/Sweden/1116/06 had D at this site.
As far as the NA gene concerned residues 34 I/V, 44 C/R, 305 S/N appeared to be discriminative of sub-clade 2.2.1./2.2.2. isolates, respectively, consistently in case of Swedish viruses and predominantly in the analyzed additional 100 sequences. Also, at NA amino acid position 305 sub-clade 2.2.1. Swedish isolates uniquely had an S while all other viruses that were analyzed (sub-clade 2.2.2. and 2.2.3. viruses) possessed N at this position due to a AA T→AG T transition. Sub-clade 2.2.2. Swedish viruses, and A/great crested grebe/Denmark/7498/06, A/grey lag goose/Denmark/6692/06, A/buzzard/Denmark/6370/06, A/tufted duck/Denmark/6540/06, and A/swan/Germany/R65/06 isolates possessed D at NA position 316 while all other analyzed viruses had G at this site.
The separation of the Swedish H5N1 HPAIVs into two subgroups was already demonstrated on the basis of NS gene sequences  and this finding was consistent for all eight genes of the isolates (herein summarized in Additional file 1). No reassortant variant was found among the sequenced twelve Swedish isolates.
Characteristic findings regarding the preservation/substitutions at particular amino acid positions along with potential distinctive molecular markers for the Swedish H5N1 viruses are summarized in Additional file 1.
A single amino acid substitution, from glutamic acid (E) to Lysine (K) in position 627 in PB2 is a determinant of mammalian host range [13, 14]. Most avian isolates have E in this position. The substitution to K in this position converts a nonlethal H5N1 influenza A virus isolated from a human to a lethal virus in mice . Among the H5N1 HPAIV sequences we investigated a larger proportion of those originating from 1998–2005 had PB2-E627 than more recent isolates. The 2.2.2.-like Swedish viruses along with the most closely related Danish and German isolates encoded K at this site while the two sub-clade 2.2.1.-like Swedish isolates (A/tufted duck/Sweden/599/06 and A/herring gull/Sweden/1116/06) possessed E at position 627. The mutations D701N and S714R in PB2 contribute to virulence by enhancing polymerase activity . All Swedish isolates had D and S at position 701 and 714, respectively.
PB1-F2 has been identified as a proapoptotic mitochondrial protein expressed from a second open reading frame of the PB1 gene  and it has been shown to contribute to viral pathogenesis in mice . Aspargine in position 66 in PB1-F2 has been demonstrated to play a key role in the pathogenicity of H5N1 viruses  and its presence was determined in all Swedish viruses. Furthermore, Swedish 2.2.1. subclade viruses had a K26Q substitution compared to 2.2.2. subclade viruses. Isolate A/tufted duck/Sweden/599/06 solely contained a T323I and a H562P, while A/herring gull/Sweden/1116/06 a V719M substitution, respectively. The H5N1 viral polymerase activity is enhanced by the presence of PB2 701N and 714R, PB1 13P, PA 615N, further, NP 319K and 678N . Among the Swedish isolates the presence of PB1 13P was determined.
Surface glycoprotein genes
The HA sequences of isolates A/Mute swan/Sweden/827/06, A/Canada goose/Sweden/978/06, and A/peregrine/Denmark/6632/06 proved to be identical. The amino acid sequence flanking the cleavage site of the HA gene was PQGERRRKKRGLF alike all other 2.2. viruses with the exception of some French isolates that had the PQGERKRKKR/G sequence . The identified amino acid markers of H5N1 influenza viruses isolated at Qinghai and Poyang Lakes from migratory birds (HA-I99, HA-N268, and NA-R110) were present in all Swedish isolates as well . No "sub-clade"-specific amino acid changes were identified in the HA among the two subgroups of Swedish isolates. All the Swedish isolates had the 238Q and 240G (numbered from the H5 start codon) which indicates preferred receptor specificity for the avian α-(2,3) linkage to galactose [19, 20]. All HA sequences contained 6 N-linked potential glycosylation sites, as analysed with NetNglyc server (threshold: 0.5) at the following positions: 27, 39, 181, 302, 500, 559; none of them is adjacent to the cleavage site. Furthermore, the substitutions S145L and A172T, which are associated with viral adaptation to poultry  were not determined in association with the Swedish H5N1 viruses.
The amino acid substitutions R178I and I248V in HA that were found in the domestic birds of the Danish isolates  were not present in any of the Swedish viruses, nor the V73I substitution that was found in the Danish swan isolates. However, the D387N substitution found in the German and most of the Danish isolates was also present in the Swedish isolates.
The H5N1 virus isolated from a mink (A/Sweden/mink/2006/V907) was examined in order to reveal any possible adaptation towards mammals. As a result, a unique E513G substitution was found in the HA gene but no substitutions that could be regarded as host-related were found, which is consistent with previous findings, i.e. that a single passage in mammals is not necessarily associated with changes in receptor-binding sites .
As in the other 2.2. viruses, NA-R110 was present in the Swedish isolates, and a 20 amino acid deletion was also found at positions 49–68 similarly to the majority of the recent H5N1 strains . The N228S substitution was present only in A/Herring gull/1116/06 Swedish 2.2.1. virus (alike with several other member of the sub-clade) and not in A/Tufted duck/Sweden/599/06 isolate. These two isolates differed further in amino acid residues 414 and 434 by bearing N/K and S/G corresponding to A/Herring gull/1116/06 and A/Tufted duck/Sweden/599/06 viruses, respectively. Interestingly, while the Danish and German isolates shared unique amino acids in the NA (G336D), PB1 (K531R) and NS2 (G63E) proteins the Swedish isolates were not homogenous in this regard: although NA-G336D was a characteristic of the Swedish viruses too, two isolates retained the PB1-531K, and NS2-63G. Reported substitutions in NA, inducing oseltamivir resistance , were not found in the Swedish isolates.
The NP and M genes
The NP-10Y amino acid residue, which may affect the pathogenicity of AIVs , was present in all of the Swedish isolates. Concerning the M2 gene, all Swedish viruses contained the L26-V27-A30-S31-G34 amino acid pattern, thus, no adamantan drug resistant variant was revealed . Substitutions S64A and E66A that were present in the M2 genes of H5N1 AIV isolates from Hong Kong  did not appear in Swedish viruses.
The complete characterization of the NS genes from these isolates was described by Zohari et al., , and is not further discussed here.
Differences in number of nucleotide and amino acid compositions, synonymous and nonsynonymous nucleotide substitutions, and predicted antigenic epitopes between sub-clade 2.2.1.-2.2.2. Swedish H5N1 avian influenza viruses.
Region of comparison/nucleotide/
Difference between sub-clade 2.2.1.-2.2.2. Swedish viruses
Number of synonymous/nonsynonymous nucleotide changes
Average number of predicted antigenic epitopes
Average number of nucleotide differences
Average number of amino acid differences
The incursion of H5N1 HPAIV strains falling into three sub-clades into Europe throughout late 2005 and 2007 has been demonstrated earlier . Further reports and the analysis of the corresponding published sequences revealed the introduction of multiple variants of H5N1 HPAIV into several European countries, such as sub-clade 2.2.1. and 2.2.2. viruses into Germany, France, and Sweden [6, 9, 11, 25], and subclade 2.2.1. and 2.2.3. viruses into Italy . The Swedish 2.2.1. sub-clade viruses were closely related to A/Cygnus olor/Italy/808/2006 and A/mallard/Italy/835/2006 and shared several common nucleotide and amino acid motifs, among them, importantly, the PB2-627E, suggesting that they derived from an earlier progenitor of Southeast Asian origin. The detection of these H5N1 HPAIV strains in Sweden adds further data regarding the spread of 2.2.1. viruses in the North. The accumulation of particular mutations reflects that presumably these viruses have been circulating in the South before the transmission to the northern parts of Europe . Sub-clade 2.2.2. Swedish H5N1 HPAIV isolates proved to be closely related to the contemporary German and Danish isolates, which supports the proposition of the introduction and spread of a single variant of 2.2.2. sub-clade H5N1 avian influenza viruses in the Baltic region.
The number and composition of the immune reactive peptides predicted by computing indicated that the surface glycoprotein genes were more affected than the nucleoprotein, polymerase basic protein 2, and polymerase acidic protein, the main targets of the cellular immune responses.
The above observations, alike those with similar objectives, highlight and warrant the importance of whole genome sequencing of HPAIV isolates, in order to improve the surveillance and preparedness against highly pathogenic avian influenza.
List of the H5N1 HPAIV isolates used in this study
RT-PCR and nucleotide sequencing
The collection of specimens, RNA extraction, and RT-PCR amplification of the NS1 sequences was described earlier and the same RNA batches were used for this study that served as targets in the previous investigation . In order to obtain possibly the full length nucleotide sequences of the coding regions of the influenza virus isolates several approaches were combined that comprised of either published protocols/primers [22, 26, 27] or those developed and used by the Influenza Genome Sequencing Project [; the primer sequences were kindly provided by David Spiro, The J. Craig Venter Institute, Rockville, Maryland, USA), or designed by ourselves. The primer and PCR protocols for sequencing are available from the authors upon request.
For the phylogenetic analyses, a set of H5N1 AIVs that were isolated in Europe, Asia and Africa in 2005 – 2006 was selected and used for all genes. These were collected from the Influenza Virus Resource at NCBI  and these were included in the phylogenetic analyses.
Sequence assembly, multiple alignment and alignment trimming were performed with the CLC Combined Workbench 3.0.2. bioinformatics software (CLC bio A/S, Aarhus, Denmark). Distance based neighbor joining and character based maximum parsimony phylogenetic trees were generated using the Molecular Evolutionary Genetics Analysis (MEGA) software v.4.0.  with 1000 bootstrap replicates. For the neighbor-joining trees, the Kimura-2 substitution model was used. Other models were also tested which showed similar topologies. The evolutionary divergence between the sub-clades was investigated by pairwise analyses over all sequence pairs between the groups by using the MEGA software also. The occurrence and distribution of synonymous and nonsynonymus substitutions was investigated by the DNA Sequence Polymorphism software (Version 4.50.3) software . Computational analysis of the antigenic sites was carried out by using the Kolaskar-Tongaonkar method .
Nucleotide sequence accession numbers
Nucleotide sequences from Swedish H5N1 virus isolates included in this study have been submitted to GenBank with the following accession numbers: PB2: EU889035–EU889046, PB1: EU889047–EU889058, PA: EU889059–EU889070, HA: EU889071–EU889082, NP: EU889083–EU889094, NA: EU889095–EU889106, M: EU889107–EU889118.
Thanks are due to Elodie Ghedin and David Spiro for their help during the set-up of the amplification protocols and to Béla Lomniczi for his comments on the manuscript. This work was partly supported by the Swedish Emergency Management Agency, the EPIZONE project (Network of Excellence for Epizootic Disease Diagnosis and Control, FP6-2004-Food-3-A), the Swedish Research Council for Environment, Agricultural Sciences and Spatial Planning (Formas 221-2006-2169 and Formas 221-2007-935) projects, and the FLUTEST EU project (Contract No.: 044429). Elisabeth Weiss was supported by the Leonardo da Vinci Mobilität Programme.
- Xu X, Subbatao 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. 10.1006/viro.1999.9820View ArticlePubMed
- Chen H, Smith GJD, Li KS, Wang J, Fan XH, Rayner JM, Vijaykrishna D, Zhang JX, Zhang LJ, Guo CT, Cheung CL, Xu KM, Duan L, Huang K, Qin K, Leung YHC, Wu WL, Lu HR, Chen Y, Xia NS, Naipospos TSP, Yuen KY, Hassan SS, Bahri S, Nguyen TD, Webster RG, Peiris JSM, Guan Y: Establishment of multiple sublineages of H5N1 influenza virus in Asia: Implications for pandemic control. Proc Natl Acad Sci USA 2006, 103: 2845-2850. 10.1073/pnas.0511120103PubMed CentralView ArticlePubMed
- Nguyen TD, Nguyen TV, Vijaykrishna D, Webster RG, Guan Y, Peiris JSM, Smith GJD: Multiple sublineages of influenza A virus (H5N1), Vietnam, 2005–2007. Emerg Infect Dis 2008., 14: [16 April 2008] [http://www.cdc.gov/EID/content/14/4/632.htm]
- World Health Organization[http://www.who.int/csr/disease/avian_influenza/Timeline_08_06_17.pdf]
- Gilbert M, Xiao X, Domenech J, Lubroth J, Martin V, Slingenbergh J: Anatidae migration in the western Palaearctic and spread of highly pathogenic avian influenza H5N1 virus. Emerg Infect Dis 2006, 12: 1650-1656.PubMed CentralView ArticlePubMed
- Weber S, Harder T, Starick E, Beer M, Werner O, Hoffmann B, Mettenleiter TC, Mundt E: Molecular analysis of highly athogenic avian influenza virus of subtype H5N1 isolated from wild birds and mammals in northern Germany. J Gen Virol 2007, 88: 554-558. 10.1099/vir.0.82300-0View ArticlePubMed
- Salzberg SL, Kingsford C, Cattoli G, Spiro DJ, Janies DA, Aly MM, Brown IH, Couacy-Hymann E, De Mia GM, Dung do H, Guercio A, Joannis T, Maken Ali AS, Osmani A, Padalino I, Saad MD, Saviæ V, Sengamalay NA, Yingst S, Zaborsky J, Zorman-Rojs O, Ghedin E, Capua I: Genome analysis linking recent European and African influenza (H5N1) viruses. Emerg Infect Dis 2007, 13: 713-718.PubMed CentralView ArticlePubMed
- Writing Committee of the Second World Health Organization Consultation on Clinical Aspects of Human Infection with Avian Influenza A (H5N1) Virus: Update on Avian Influenza A (H5N1) Virus Infection in Humans. N Engl J Med 2008, 358: 261-273. 10.1056/NEJMra0707279View Article
- Starick E, Beer M, Hoffmann B, Staubach C, Werner O, Globig A, Strebelow G, Grund C, Durban M, Conraths FJ, Mettenleiter T, Harder T: Phylogenetic analyses of highly pathogenic avian influenza virus isolates from Germany in 2006 and 2007 suggest at least three separate introductions of H5N1 virus. Vet Microbiol 2008, 128: 243-252.View ArticlePubMed
- Zohari S, Gyarmati P, Thorén P, Czifra Gy, Bröjer C, Belák S, Berg M: Genetic characterization of the NS gene indicates co-circulation of two sub-lineages of highly pathogenic avian influenza virus of H5N1 subtype in Northern Europe in 2006. Virus Genes 2008, 36: 117-25. 10.1007/s11262-007-0188-7View ArticlePubMed
- Le Gall-Reculé G, Briand FX, Schmitz A, Guionie O, Massin P, Jestin V: Double introduction of highly pathogenic H5N1 avian influenza virus into France in early 2006. Avian Pathology 2008, 37: 15-23. 10.1080/03079450701774835View Article
- Lu G, Rowley T, Garten R, Donis RO: FluGenome: a web tool for genotyping influenza A virus. Nucleic Acids Res 2007, (35 Web Server):W275-279. 10.1093/nar/gkm365
- Shinya K, Hamm S, Hatta M, Ito H, Ito T, Kawaoka Y: PB2 amino acid at position 627 affects replicative efficiency, but not cell tropism, of Hong Kong H5N1 influenza A viruses in mice. Virology 2004, 320: 258-266. 10.1016/j.virol.2003.11.030View ArticlePubMed
- Subbarao EK, London W, Murphy BR: A single amino acid in the PB2 gene of influenza A virus is a determinant of host range. J Virol 1993, 67: 1761-1764.PubMed CentralPubMed
- Gabriel G, Dauber B, Wolff T, Planz O, Klenk H-D, Stech J: The viral polymerase mediates adaptation of an avian influenza virus to a mammalian host. Proc Natl Acad Sci USA 2005, 102: 18590-18595. 10.1073/pnas.0507415102PubMed CentralView ArticlePubMed
- Chen W, Calvo PA, Malide D, Gibbs J, Schubert U, Bacik I, Basta S, O'Neill R, Schickli J, Palese P, Henklein P, Bennink JR, Yewdell JW: A novel influenza A virus mitochondrial protein that induces cell death. Nat Med 2001, 7: 1306-1312. 10.1038/nm1201-1306View ArticlePubMed
- Zamarin D, Ortigoza MB, Palese P: Influenza A virus PB1-PF2 protein contributes to viral pathogenesis in mice. J Virol 2006, 80: 7976-7983. 10.1128/JVI.00415-06PubMed CentralView ArticlePubMed
- Conenello GM, Zamarin D, Perrone LA, Tumpey T, Palese P: A Single Mutation in the PB1-F2 of H5N1 (HK/97) and 1918 Influenza A Viruses Contributes to Increased Virulence. PLoS Pathog 2007,3(10):e141. 10.1371/journal.ppat.0030141PubMed CentralView Article
- Matrosovich MN, Matrosovich TY, Gray T, Roberts NA, Klenk HD: Human and avian influenza viruses target different cell types in cultures of human airway epithelium. Proc Natl Acad Sci USA 2004, 101: 4620-4624. 10.1073/pnas.0308001101PubMed CentralView ArticlePubMed
- Suzuki Y, Ito T, Suzuki T, Holland RE Jr, Chambers TM, Kiso M, Ishida H, Kawaoka Y: Sialic acid species as a determinant of the host range of influenza A viruses. J Virol 2000, 74: 11825-11831. 10.1128/JVI.74.24.11825-11831.2000PubMed CentralView ArticlePubMed
- Perdue ML, Suarez DL: Structural features of the avian influenza virus hemagglutinin that influence virulence. Veterinary Microbiology 2000, 74: 77-86. 10.1016/S0378-1135(00)00168-1View ArticlePubMed
- Bragstad K, Jorgensen PH, Handberg K, Hammer AS, Kabell S, Fomsgaard A: First introduction of highly pathogenic H5N1 avian influenza viruses in wild and domestic birds in Denmark, Northern Europe. Vir J 2007, 4: 43. 10.1186/1743-422X-4-43View Article
- Subbarao K, Joseph T: Scientific barriers to developing vaccines against avian influenza viruses. Nature Reviews – Immunology 2007, 7: 267-278. 10.1038/nri2054View ArticlePubMed
- Kolaskar AS, Tongaonkar PC: A semi-empirical method for prediction of antigenic determinants on protein antigens. FEBS Lett 1990, 276: 172-4. 10.1016/0014-5793(90)80535-QView ArticlePubMed
- Rinder M, Lang V, Fuchs C, Marx-Hafner A, Bogner KH, Neubauer A, Büttner M, Rinder H: Genetic evidence for multi-event imports of avian influenza virus A (H5N1) into Bavaria, Germany. J Vet Diagn Invest 2007, 19: 279-282.View ArticlePubMed
- Hoffmann E, Stech J, Guan Y, Webster RG, Perez DR: Universal primer set for the full-length amplification of all influenza A viruses. Arch Virol 2001, 146: 2275-2289. 10.1007/s007050170002View ArticlePubMed
- Li OT, Barr I, Leung CY, Chen H, Guan Y, Peiris JS, Poon LL: Reliable universal RT-PCR assays for studying influenza polymerase subunit gene sequences from all 16 haemagglutinin subtypes. J Virol Methods 2007, 142: 218-22. 10.1016/j.jviromet.2007.01.015View ArticlePubMed
- Ghedin E, Sengamalay NA, Shumway M, Zaborsky J, Feldblyum T, Subbu V, Spiro DJ, Sitz J, Koo H, Bolotov P, Dernovoy D, Tatusova T, Bao Y, St George K, Taylor J, Lipman DJ, Fraser CM, Taubenberger JK, Salzberg SL: Large-scale sequencing of human influenza reveals the dynamic nature of viral genome evolution. Nature 2005, 437: 1162-1166. 10.1038/nature04239View ArticlePubMed
- Influenza Virus Resource2004. [http://www.ncbi.nlm.nih.gov/genomes/FLU/FLU.html]
- Tamura K, Dudley J, Nei M, Kumar S: Molecular Evolutionary Genetics Analysis (MEGA) software version 4.0. Mol Biol Evol 2007, 24: 1596-1599. 10.1093/molbev/msm092View ArticlePubMed
- Rozas J, Sanchez-DelBarrio JC, Messeguer X, Rozas R: DnaSP, DNA polymorphism analyses by the coalescent and other methods. Bioinformatics 2003, 19: 2496-2497. 10.1093/bioinformatics/btg359View ArticlePubMed
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.