- Open Access
First outbreaks and phylogenetic analyses of avian influenza H9N2 viruses isolated from poultry flocks in Morocco
Virology Journal volume 13, Article number: 140 (2016)
H9N2 avian influenza viruses continue to spread in poultry and wild birds worldwide. Morocco just faced its first H9N2 influenza virus outbreaks early 2016 affecting different types of poultry production. After its introduction, the virus spread very rapidly throughout the country.
Samples were collected from 11 chicken flocks with high morbidity and mortality rates. Four viruses were successfully isolated from broiler chickens and one from broiler breeders and fully sequenced.
Phylogenetic and molecular markers analyses showed the Moroccan viruses belonged to the G1 lineage and likely originated from the Middle East. As known for H9N2 viruses, the Moroccanisolates possess several genetic markers that enhance virulence in poultry and transmission to humans.
The present study demonstrated that under field conditions H9N2 could have a devastating effect on egg production and mortalities and highlighted a lack of surveillance data on the pathogen in the region.
The influenza A virus (IAV) H9N2 is of great economical concern for poultry production and public health . IAV H9N2 is now enzootic in many parts of the world. In chickens, the first report dates from 1992 in China [2, 3] and the virus was reported in South Korea in 1996  and Iran and Pakistan in 1998 and 1999 respectively [5, 6]. The virus spread in the early 2000s to several countries of the Middle East, including Israel, Saudi Arabia, United Arab Emirates and Sultanate of Oman [7–9]. Then, the virus appeared in North African countries: Egypt, Libya and Tunisia [10, 11].
Despite the fact that H9N2 is considered as low pathogenic virus, it has been reported that several microbial agents of the respiratory tract and some environmental factors exacerbate IAV H9N2 infections, leading to very severe respiratory disease and causing mortality up to 65 % in broiler chickens and up to 70 % drop in egg production in layers and breeders [12, 13]. It has also been reported that IAV H9N2 affected many species of domesticated and wild birds such as falcon, partridge, quail, houbara, pigeon, sparrow and other species. In Pakistan, sparrows were found to play a very important role in the transmission of the virus between farms . IAV H9N2 was detected in mammalian species such as dogs and cats . In affected countries, high seroprevalences to H9N2 virus were detected in workers of poultry farms and poultry processing plants . The virus itself has been isolated from Human and is monitored continuously by the World Health Organization (WHO).
Molecular analyses of H9N2 viruses isolated during the last two decades revealed that these viruses evolve fast and represent a genetically diverse population . H9N2 influenza viruses circulate in poultry and wild birds worldwide and cluster into two main lineages in Eurasia: A/quail/Hong Kong/G1/97-like viruses (G1-like) and A/duck/Hong Kong/Y280/97-like viruses (Y280-like) .
In Morocco, low pathogenic avian influenza subtype H9N2 virus was first diagnosed in poultry in January 2016. As an initial response, the National Sanitary and Security Food Office of Morocco authorized emergency vaccinations of all types of poultry production. This study describes the first report and rapid spread of IAV H9N2 to all areas of the country within a few weeks and its devastating effects in different types of poultry productions. Genetic analyses of Moroccan H9N2 isolates will be presented.
11 poultry flocks from different regions of Morocco showing very high mortality, decrease in feed consumption and very severe respiratory signs including sneezing, coughing, rales and gasping were analysed by real time RT-PCR for the presence of nucleic acids of infectious bronchitis, Newcastle disease, virulent infectious bursal disease viruses and IAV [10, 19–21]. Details on the flocks sampled in this study are presented in Table 1.
Five samples out of 11 with low Ct values in real time RT-PCR for IAV were selected for inoculation into the allantoic cavities of 9–day–old specific pathogen free chickens embryos, incubated for 48 h at 37 °C then chilled at 4 °C for 4 h before harvesting. The allantoic fluid was harvested, clarified, tested by real time RT-PCR for IAV genome, and then stored at −80 °C until use.
Full genome amplification and sequencing of H9N2 Moroccan isolates
PCR amplification was performed in a thermocycler (GeneAmp PCR System 9700, Applied Biosystems) using segment specific primers described previously by Hoffman et al., (2001)  and the One step RT-PCR kit (Qiagen) following instructions of the manufacturer. Reactions were performed according to the following protocol: 50 °C for 30 min and 95 °C for 5 min, followed by 40 cycles of 95 °C for 30 s, 54 °C to 57 °C for 30 s, 72 °C for 1 min30 s, and a final elongation step of 6 min at 72 °C. PCR products of the expected length were purified with the Nucleospin gel and PCR cleanup kit (Macherey Nagel) according to the manufacturer’s instructions. Sanger sequencing method was carried out on all the segments using the same primers as used for the RT-PCR. Sequencing was performed on a 3130XL Applied Biosystems capillary sequencer at the Plateau de Génomique GeT-Purpan, UDEAR UMR 5165 CNRS/UPS, CHU PURPAN, Toulouse, France.
Sequences and phylogenetic analyses
Complete sequence data of five H9N2 isolates were manually assembled using BioEdit software package version 5.0.9 . The open source BLAST program (National Center for Biotechnology Information, Bethesda MD, http://blast.ncbi.nlm.nih.gov/Blast.cgi) was used for sequence comparison. The best fitting nucleotide substitution model was estimated by means of hierarchical likelihood ratio approach using Mega 6.06  for each sequences alignment. Phylogenetic analysis and tree construction were thus generated using the maximum likelihood method, General Time Reversible model, with 500 bootstrap replicates with Mega 6.06, and bootstrap values above 50 were labelled on major tree branches. Sequences were submitted to the EMBL/GenBank databases.
Epidemiology of H9N2 outbreaks in Morocco
Up to the end of December 2015, no official report on the presence of any subtype of avian influenza in commercial poultry farms or backyard flocks had been recorded in Morocco. During the second week of January 2016, 11cases of chickens from different poultry production areas around the cities of Fez, Meknes, Kenitra and Rabat, were received for diagnosis purpose in our laboratory (Table 1). Most cases showed very high decrease in feed consumption and very severe respiratory signs (sneezing, coughing, rales and gasping). All cases showed a very rapid increase in mortality rates (2 % to 10 % per day between the 3rd and the 6th day after onset of disease). All affected flocks were over three weeks of age and some flocks were at the selling period. In a few cases, diarrhoea was reported. All flocks were vaccinated, as routinely performed, against Newcastle, infectious bronchitis and Gumboro diseases. By mid-February, all parts of the country were affected and the disease spread to layers and breeders farms, which showed very severe respiratory signs, mortality rates ranging from 2 to 15 %and drastic drop in egg production (up to 80 %) with no complete recovery after several weeks. Indeed, up to 10 weeks after the infection, egg production remained between 20 and 25 % lower than normal levels. Turkey flocks were similarly affected at different ages, with similar clinical signs as seen in chickens as well as sinusitis with facial swellings, and mortality rates around 10 %, but with a large inter-flock variability, likely due to other contributing factors.
Post mortem examinations revealed that all carcasses were congested, some showing oculonasal discharge and infraorbital sinusitis. The main consistent lesions were congested or haemorrhagic tracheas with fibrinous exudates in the lumen and presence of fibrinous casts in bronchi (Fig. 1), haemorrhagic lungs, pneumonia with fibrinous exudates and airsaculitis. In some cases congested and swollen kidneys, spleen and liver were observed.
Real time RT-PCR
Eleven samples from respiratory tissues collected from different cases and analysed by real time RT-PCR for the presence of Infectious bronchitis, Newcastle disease and virulent IBD viruses werenegative. However10 out of 11 samples were positive for IAV with cycle threshold (Ct) values ranging from 17.5 to 35.1 (Table 1).
The phylogenetic analyses showed that the Moroccan H9N2 viruses of broiler and breeder broiler origins were of the G1 lineage and closely related to each other (Fig. 2). The nucleotide sequence identities among the five Moroccan strains ranged from 98 to 100 %. Genes segments coding for polymerase basic 2 (PB2), polymerase basic 1 (PB1), acidic polymerase (PA), hemagglutinine (HA), nucleoprotein (NP), neuraminidase (NA), matrix (M) and non structural (NS) proteins all clustered with G1 lineage viruses, closely related to an A/pheasant/United Arab Emirates/D1521/2011-like H9N2 virus, as indicated by the trees topology (Figs. 2, 3) and by the high nucleotide sequence identity: 96.1 to 98.2 %.
Hemagglutinin cleavage site and receptor binding domains
The HA of Moroccan H9N2 isolates analyzed in this study had a RSSR*GLF motif at the HA0 cleavage site, a characteristic of the LPAI viruses in poultry. The HA proteins harboured seven potential glycosylation sites with one additional glycosylation site (residue 82) as compared to recently characterized Egyptian H9N2 viruses (Table 2). All the Moroccan viruses had 158 N, H183, 226 L and 391 K substitutions in their receptor binding site (H3 numbering) known to confer preference to humanlike receptors [25–27].
NA structural features
The NA sequences of the five Moroccan isolates had a nucleotide sequence identity between 99.8 and 100 % when compared to each other. The deduced NA amino acid sequence did not have any previously identified mammalian-adaptation associated substitutions and no stalk deletion (Table 2).
Molecular features of internal genes
Three mammalian adaptation markers were identified in Moroccan PB2 proteins: 318R, 590S and 661 T [28–30]. An additional substitution, 504 V, was observed in all 5 isolates. A unique molecular determinant of host specificity described as avian–human signature 13Pwas observed for PB1 . A truncated 52 amino acids PB1-F2 protein was harbored by the 5 Moroccan isolates. The PA proteins had three amino acid substitutions associated within creased virulence: 127 V, 550 L and 672 L as previously reported by Kandeil et al., 2014, Rolling et al., 2009 and Chen et al., 2006 [27, 31, 32]. All Moroccan isolates had 409 N, 100I and 312R substitutions, which are associated with mammalian adaptation . The deduced NP amino acid sequence had 372D substitution previously identified as important for changing host range from avian to human . On the M2 protein, the 10 L residue, marker of the G1 lineage, was conserved. All Moroccan isolates also harbored as substitution at M2 amino acid position 31:S31N, suggesting a resistance to the amantadane [8, 30]. None of the Moroccan isolates contained deletion at amino acid positions 80–84 in NS1. Moreover, all Moroccan isolates possessed the 42S and 189D virulence markers and the NS1 GSEV motif in NS1 (Table 2). The molecular characterization of NS2 protein revealed the presence of 31 M and 56H markers: two known virulence determinants.
This study illustrates clearly how fast the H9N2 IAV did spread to the whole country within a very short time likely through the movement of vehicles transporting live birds to live bird markets (over 70 % of the poultry trade in Morocco) and feed delivery vehicles as previously reported in other countries in the world [27, 30]. According to Wernery et al., (2013) smuggled H9N2-infected falconry birds may contribute to the spread of these viruses to wild birds, domestic poultry, and humans .
It was important to notice that common biosecurity measures applied routinely in breeder and layer production sites were insufficient to prevent the entry of virus. Although H9N2 viruses are considered as low pathogenic, the recent Moroccan situation showed that it could cause indeed very heavy losses in different types of poultry production systems. In order to limit disease spread to the other parts of the country, several measures were quickly implemented including the importation of inactivated H9N2 vaccine to immunise all types of poultry productions. H9N2 viruses are prevalent in terrestrial poultry throughout Asia and have been isolated from poultry outbreaks worldwide. In addition to the economic losses related to drop in egg production and mortalities in different type of poultry flocks, the zoonotic potential of the pathogen should not be under estimated . In fact, H9N2 viruses infect both avian and mammalian species and might constitute significant donors of genetic material to emerging human type A influenza viruses [33, 34]. It was reported that the natural isolates of H9N2 viruses were transmissible among ferrets via respiratory droplets, emphasizing the pandemic potential of H9N2 viruses  and human cases of H9N2 infections have been reported by WHO .
In the present study our findings clearly demonstrate a relationship between viruses that recently emerged in Morocco and viruses isolated in the Middle East of the G1 lineage. In order to evaluate the genetic similarity of the Moroccan and Middle Eastern viruses, phylogenetic and genetic identities were calculated for each gene segment. The limited influenza virus sequence data available from North African and Middle-Eastern countries makes the identification of the origin of Moroccan H9N2 virus outbreaks very difficult. The similarities between our viruses and those from the United Arab Emirates and Israel confirm that the viruses isolated from Morocco in 2016 are the likely direct progenitors of those circulating in the Middle East [9, 27, 30, 37]. On the HA phylogeny (Fig. 2), Moroccan viruses indeed clustered between 2011 United Arab Emirates and 2011–2014 Israeli viruses. For the other 7 gene segments, the Moroccan isolates sequences formed a distinct cluster with pheasant and white-bellied bustard viruses isolated in the United Arab Emirates (supported by high bootstrap values on all trees, over 99 %, Figs. 2 and 3). However only HA sequences are available for the Israeli viruses and the most recent Emirati viruses sequences are from 2011: the clear lack of surveillance data in the region may considerably bias the analysis.
Molecular characterization of the HA protein of Moroccan isolates revealed one addition potential glycosylation site as compared to recent Egyptian H9N2 HAs. Previous studies showed a link between additional glycosylation sites on HA and the adaptation of H9N2 in poultry . The presence of Q226L mutation in HA is known to enable binding to human-like receptors . H9N2 viruses in Morocco possess several genetic markers shown to enhance virulence in poultry and enable transmission to humans. Previously, several studies demonstrated that G1-like H9N2 viruses were the likely donors of the six internal genes of the H5N1 viruses causing the bird flu outbreak in humans in 1997 in Hong Kong . Moreover, H9N2 influenza viruses continue to reassort with circulating H7N9 viruses, resulting in multiple genotypes of H7N9 viruses [2, 34, 41]. Full genome sequencing of H9N2 influenza viruses is therefore essential to better predict putative zoonotic transmissions. How the disease was introduced to the country is still not understood and several investigations are currently in progress.
The present study investigated the first introduction and isolation of avian H9N2 in Morocco. Although H9N2 viruses are considered as low pathogenic, the recent Moroccan situation showed that it could cause indeed very heavy losses in different types of poultry production. Phylogenetic analyses showed that the Moroccan viruses belonged to the G1 lineage and likely originated from the Middle East. How the disease was introduced to the country is still not understood.
HA, hemagglutinine; IAV, Influenza A virus; M, matrix; NA, neuraminidase; NP, nucleoprotein; NS, non structural; PA, acidic polymerase; PB1, polymerase basic 1; PB2, polymerase basic 2; RT-PCR, reverse transcription polymerase chain reaction; WHO, World Health Organization
Butt KM, Smith GJ, Chen H, et al. Human infection with an avian H9N2 influenza A virus in Hong Kong in 2003. J Clin Microbiol. 2005;43:5760–67.
Sun Y, Liu J. H9N2 influenza virus in China: a cause of concern. Protein Cell. 2015;6:18–25.
Li C, Yu K, Tian G, Yu D, Liu L, Jing B, et al. Evolution of H9N2 influenza viruses from domestic poultry in Mainland China. Virology. 2005;340:70–83.
Lee D-H, Song C-S. H9N2 avian influenza virus in Korea: evolution and vaccination. Clin Exp Vaccine Res. 2013;2:26.
Naeem K, Ullah A, Manvell RJ, Alexander DJ. Avian influenza A subtype H9N2 in poultry in Pakistan. Vet Rec. 1999;145:560.
Nili H, Asasi K. Avian influenza (H9N2) outbreak in Iran. Avian Dis. 2003;47:828–31.
Perk S, Banetnoach C, Shihmanter E, Pokamunski S, Pirak M, Lipkind M, et al. Genetic characterization of the H9N2 influenza viruses circulated in the poultry population in Israel. Comp Immunol Microbiol Infect Dis. 2006;29:207–23.
Aamir UB, Wernery U, Ilyushina N, Webster RG. Characterization of avian H9N2 influenza viruses from United Arab Emirates 2000 to 2003. Virology. 2007;361:45–55.
Body MH, Alrarawahi AH, Alhubsy SS, Saravanan N, Rajmony S, Mansoor MK. Characterization of Low Pathogenic Avian Influenza Virus Subtype H9N2 Isolated from Free-Living Mynah Birds (Acridotheres tristis) in the Sultanate of Oman. Avian Dis. 2015;59:329–34.
Monne I, Ormelli S, Salviato A, De Battisti C, Bettini F, Salomoni A, et al. Development and Validation of a One-Step Real-Time PCR Assay for Simultaneous Detection of Subtype H5, H7, and H9 Avian Influenza Viruses. J Clin Microbiol. 2008;46:1769–73.
Kammon A, Heidari A, Dayhum A, Eldaghayes I, Sharif M, Monne I, et al. Characterization of Avian Influenza and Newcastle Disease Viruses from Poultry in Libya. Avian Dis. 2015;59:422–30.
Seifi S, Asasi K, Mohammadi A. An experimental study on broiler chicken co-infected with the specimens containing avian influenza (H9 subtype) and infectious bronchitis (4/91 strain) viruses. Iran J Vet Res. 2012;13:138–42.
Azizpour A, Goudarzi H, Charkhkar S, Momayez R, Hablolvarid MH. Experimental study on tissue tropism and dissemination of H9N2 avian influenza virus and Ornithobacterium rhinotracheale co-infection in SPF chickens. JAPS J Anim Plant Sci. 2014;24:1655–62.
Iqbal M, Yaqub T, Mukhtar N, Shabbir MZ, McCauley JW. Infectivity and transmissibility of H9N2 avian influenza virus in chickens and wild terrestrial birds. Vet Res. 2013;44:100.
Zhang K, Zhang Z, Yu Z, Li L, Cheng K, Wang T, et al. Domestic cats and dogs are susceptible to H9N2 avian influenza virus. Virus Res. 2013;175:52–7.
Khan SU, Anderson BD, Heil GL, Liang S, Gray GC. A Systematic Review and Meta-Analysis of the Seroprevalence of Influenza A(H9N2) Infection Among Humans. J Infect Dis. 2015;212:562-9.
Group TSHW. Assessing the fitness of distinct clades of influenza A (H9N2) viruses. Emerg Microbes Infect. 2013;2:e75.
Xu KM, Li KS, Smith GJD, Li JW, Tai H, Zhang JX, et al. Evolution and Molecular Epidemiology of H9N2 Influenza A Viruses from Quail in Southern China, 2000 to 2005. J Virol. 2007;81:2635–45.
Wise MG, Suarez DL, Seal BS, Pedersen JC, Senne DA, King DJ, et al. Development of a Real-Time Reverse-Transcription PCR for Detection of Newcastle Disease Virus RNA in Clinical Samples. J Clin Microbiol. 2004;42:329–38.
Meir R, Maharat O, Farnushi Y, Simanov L. Development of a real-time TaqMan® RT-PCR assay for the detection of infectious bronchitis virus in chickens, and comparison of RT-PCR and virus isolation. J Virol Methods. 2010;163:190–4.
Tomás G, Hernández M, Marandino A, Panzera Y, Maya L, Hernández D, et al. Development and validation of a TaqMan-MGB real-time RT-PCR assay for simultaneous detection and characterization of infectious bursal disease virus. J Virol Methods. 2012;185:101–7.
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–89.
Hall TA. BioEdit: a user-friendly biological sequence alignment editor and analysis program for Windows 95/98/NT. Nucleic Acids Symp. Ser. 1999;95–8
Tamura K, Stecher G, Peterson D, Filipski A, Kumar S. MEGA6: Molecular Evolutionary Genetics Analysis Version 6.0. Mol Biol Evol. 2013;30:2725–9.
Wagner R, Wolff T, Herwig A, Pleschka S, Klenk HD. Interdependence of hemagglutinin glycosylation and neuraminidase as regulators of influenza virus growth: a study by reverse genetics. J Virol. 2000;74:6316–23.
Wanzeck K, Boyd KL, McCullers JA. Glycan shielding of the influenza virus hemagglutinin contributes to immunopathology in mice. Am J Respir Crit Care Med. 2011;183:767–73.
Kandeil A, El-Shesheny R, Maatouq AM, Moatasim Y, Shehata MM, Bagato O, et al. Genetic and antigenic evolution of H9N2 avian influenza viruses circulating in Egypt between 2011 and 2013. Arch Virol. 2014;159:2861–76.
Shaw M, Cooper L, Xu X, Thompson W, Krauss S, Guan Y, et al. Molecular changes associated with the transmission of avian influenza a H5N1 and H9N2 viruses to humans. J Med Virol. 2002;66:107–14.
Pan C, Cheung B, Tan S, Li C, Li L, Liu S, et al. Genomic signature and mutation trend analysis of pandemic (H1N1) 2009 influenza A virus. PLoS One. 2010;5:e9549.
Wernery U, Shanmuganatham KK, Krylov PS, Joseph S, Friedman K, Krauss S, et al. H9N2 influenza viruses from birds used in falconry. Influenza Other Respir Viruses. 2013;7:1241–5.
Chen G-W, Chang S-C, Mok C, Lo Y-L, Kung Y-N, Huang J-H, et al. Genomic signatures of human versus avian influenza A viruses. Emerg Infect Dis. 2006;12:1353–60.
Rolling T, Koerner I, Zimmermann P, Holz K, Haller O, Staeheli P, et al. Adaptive mutations resulting in enhanced polymerase activity contribute to high virulence of influenza A virus in mice. J Virol. 2009;83:6673–80.
Kageyama T, Fujisaki S, Takashita E, Xu H, Yamada S, Uchida Y, et al. Genetic analysis of novel avian A (H7N9) influenza viruses isolated from patients in China, February to April 2013. Euro Surveill. 2013;18:7–21.
Wang Y, Yuan X, Qi L, Zhang Y, Xu H, Yang J, et al. H9N2 avian influenza virus-derived natural reassortant H5N2 virus in swan containing the hemagglutinin segment from Eurasian H5 avian influenza virus with an in-frame deletion of four basic residues in the polybasic hemagglutinin cleavage site. Infect Genet Evol. 2016;40:17–20.
Li X, Shi J, Guo J, Deng G, Zhang Q, Wang J, et al. Genetics, receptor binding property, and transmissibility in mammals of naturally isolated H9N2 Avian Influenza viruses. PLoS Pathog. 2014;10:e1004508.
WHO, Influenza at the human-animal interface, Summary and assessment, 5 April to 9 May 2016, 2016. http://www.who.int/influenza/human_animal_interface/Influenza_Summary_IRA_HA_interface_05_09_2016.pdf?ua=1, Accessed 12 July 2016.
Roussan DA, Khawaldeh GY, Al Rifai RH, Totanji WS, Shaheen IA. Avian influenza virus H9 subtype in poultry flocks in Jordan. Prev Vet Med. 2009;88:77–81.
Baigent SJ, McCauley JW. Influenza type A in humans, mammals and birds: Determinants of virus virulence, host-range and interspecies transmission. BioEssays. 2003;25:657–71.
Wan H, Perez DR. Amino Acid 226 in the Hemagglutinin of H9N2 Influenza Viruses Determines Cell Tropism and Replication in Human Airway Epithelial Cells. J Virol. 2007;81:5181–91.
Guan Y, Shortridge KF, Krauss S, Webster RG. Molecular characterization of H9N2 influenza viruses: were they the donors of the “internal” genes of H5N1 viruses in Hong Kong? Proc Natl Acad Sci. 1999;96:9363–7.
Wu D, Zou S, Bai T, Li J, Zhao X, Yang L, et al. Poultry farms as a source of avian influenza A (H7N9) virus reassortment and human infection. Sci Rep. 2015;5:7630.
The project was supported by the French-Moroccan Programme Hubert Curien Toubkal/16/25 - Campus France 34654NLand by the Agronomy and Veterinary Institute Hassan II, Rabat, Morocco.
The Agronomy and Veterinary Institute Hassan II was sponsored the screening and isolation of the different samples. The IHAP, Université de Toulouse, INRA, ENVT, was supported the full genome sequencing of the five isolates.
Availability of data and material
Data supporting the conclusions of this article are presented in the manuscript.
M.E. participated in the design of the study and drafted the manuscript. S.F. participated in the study, carried out the molecular genetic studies and the sequence alignment phylogenetic. S.N. helped to draft the manuscript. J.L.G. corrected the manuscript. M.D. carried out the sequencing of the isolates and molecular analyses, participated in the design of the study, and helped draft the manuscript. All authors read and approved the final manuscript.
The authors declare that they have no competing interests.
Ethical approval and consent to participate
This article does not contain any studies with human participants or animals performed by any of the authors.
About this article
Cite this article
EL Houadfi, M., Fellahi, S., Nassik, S. et al. First outbreaks and phylogenetic analyses of avian influenza H9N2 viruses isolated from poultry flocks in Morocco. Virol J 13, 140 (2016) doi:10.1186/s12985-016-0596-1
- Influenza Virus
- Phylogenetic analyses