Molecular identification of potential leishmaniasis vector species within the Phlebotomus (Euphlebotomus) argentipes species complex in Sri Lanka
© Gajapathy et al.; licensee BioMed Central Ltd. 2013
Received: 19 July 2013
Accepted: 20 September 2013
Published: 18 October 2013
Leishmaniasis is an emerging vector-borne disease in Sri Lanka. Phlebotomus (Euphlebotomus) argentipes sensu lato Annandale and Brunette 1908 is suspected to be a potential vector. Three sibling species have been reported in the species complex based on analysis of morphological data. A study was carried out in different parts of Sri Lanka including cutaneous leishmaniasis prevailing localities to characterise the sibling species of Phlebotomus (Euphlebotomus) argentipes sensu lato and to establish their possible role in Leishmania transmission.
Sandflies were collected using cattle baited trap nets and mouth aspirator. They were identified based on existing taxonomic keys. Sequences of amplified cytochrome oxidase subunit I (CO I), cytochrome oxidase b (cyt b), internal transcribed spacer 2 (ITS2), 18s and 28s rDNA regions were analysed to confirm the number of sibling species. Vectorial capacity of the sibling species was checked by detecting human and Leishmania DNA.
Sandflies collected using different techniques were processed for identification, parasite detection and molecular characterization. The 18s, 28s rDNA and cytochrome oxidase subunit I (CO I), internal transcribed spacer 2 (ITS2) and cytochrome b oxidase (cytb) sequences confirmed that the species belonged to the Argentipes complex. 18s and 28s sequences did not show any variation among the proposed sibling species. The phylogeny created from mitochondrial CO I and cytochrome b data and from the nuclear ITS2 region supports the existence of only two groups of flies (termed A and B) from Phlebotomus (Euphlebotomus) argentipes complex instead of the previously proposed three. The Leishmania mini-circle kinetoplastid, heat shock protein 70 (hsp70) and internal transcribed spacer I DNA along with human blood were detected from sibling species A only, which has not previously been considered to be a vector.
The taxonomy of the Sri Lankan Argentipes species complex is reassessed based on the molecular data. The existence of two sibling species is proposed; sibling species A has a long sensilla chaetica (> 50% length of the second antennal flagellomere) and sibling species B has a short sensilla cheatica (< 50%). Sibling species A is incriminated as a vector for leishmaniasis in Sri Lanka.
Leishmaniasis is a health threat in as many as 88 countries[1, 2]. The common forms of the disease in the old-world are cutaneous leishmaniasis (CutL), visceral leishmaniasis (VisL) and post Kala-azar dermal leishmaniasis (PKDL). The disease forms are prevalent in different geographic regions. Ninety percent of previously reported VisL patients were from Bangladesh, India, Nepal, Sudan and Brazil whereas around ninety percent of the CutL cases were recorded from areas such as Afghanistan, Algeria, Brazil, Iran, Peru, Saudi Arabia and Syria.
Leishmaniasis was previously considered to be an exotic disease in Sri Lanka. Migrant workers returning from the Middle East were the only diagnosed patients before the early 1990’s. The first case of autochthonous CutL was reported in Mamandala village of Hambantota district in 1989. Since then, the number of cases of CutL has risen and more than 2000 cases have been reported in the last decade from many parts of the country.
The parasite causing CutL in Sri Lanka has been identified as Leishmania donovani zymodeme MON 37. Leishmania donovani in the Donovani species complex is generally associated with VisL and PKDin India and Afrotropical regions[5, 8]. The suspected vector(s) in Sri Lanka are sandflies from the Phlebotomus argentipes sensu lato Annandale & Brunette, 1908 species complex, which are known to be the vector for Leishmania donovani in India.
The occurrence of sibling or cryptic species among insect vectors that have overlapping morphological characters is very well documented. The inability to distinguish these sibling species from one other using standard morphological analysis has led to the use of biochemical tools in other organisms such as Anopheles mosquitoes. These include allozyme analysis and polytene chromosome banding patterns[10, 11]. More recently DNA sequence analysis of conserved regions of Ribosomal DNA (rDNA) or mitochondrial markers (e.g. cytochrome oxidase subunit I or cytochrome b oxidase) have been used. In many cases, only one or two sibling species within a species complex have vectorial capacity. Correct identification of vector sibling species is important as failure to do so may conceal the actual transmission pattern of the disease and will result in inadequate vector/disease control strategy. This has been well demonstrated in cases such as the differential insecticidal resistance observed among sympatric sibling species in Anopheles mosquitoes by[13, 14]. Different micro geographic forms of Phlebotomus papatasi in Sudan and Egypt are suspected to play different role in leishmaniasis epidemiology.
In this study, different methods are used to characterize the relationships between flies within the Phlebotomus (Euphlebotomus) argentipes species complex in Sri Lanka and to establish their potential to be vectors for Leishmania.
Sandfly collection and identification
DNA of the individual flies was extracted by the ethanol precipitation method. Individual flies were homogenized in extraction buffer (Tris base (pH 8.00), NaCl, SDS and EDTA) and proteinase K (4:1). 5M NaCl was added after the overnight incubation at 56°C. The supernatant was separated after centrifuging at 14000 rpm for 5 minutes. 400 μl of cold, 70% ethanol was added to the supernatant. The solution was incubated at -20°C for 1 hour. The solution was centrifuged at 14000 rpm for 30 minutes. Ethanol was removed and the pellet of DNA was washed with cold 70% ethanol. After the removal of ethanol and drying, the pellet was re-suspended with 120 μl of double distilled water.
Specific primers were used to amplify the D3 region of the large subunit ribosomal RNA gene and a section of the 18S small subunit rDNA gene in addition to a variable region of the mitochondrial cytochrome oxidase I (COI) gene. The PCR reaction mixture was prepared as follows: 1X Taq polymerase buffer, 1.5 mM MgCl2, 200 μM dNTPs and 0.625U Taq polymerase (Promega, USA). 0.125 μM (for D3) and 0.6125 μM primers were used in the PCR. The reaction conditions were; 94°C of initial denaturation followed by 94°C (30s), 55°C (45s) and 72°C (1 min.) for 35 cycles and 72°C for 10 minutes for D3. 18S rDNA was amplified using the conditions described by Surendran et al.. The cytochrome oxidase subunit I gene was amplified with an initial denaturation of 94°C (10 min.) and 40 cycles of 94°C (30s), 50°C (30s) and 72°C (1 min.) and a final extension at 72°C for 10 minutes. The cytochrome oxidase b gene was amplified with degenerate primers as previously described. 2.0 mM MgCl2, 1.0U Taq DNA polymerase, 0.1 μM of each primer, 200 μM of dNTPs, 1X Taq buffer and 4 μl of genomic DNA was used in a total volume of 40 μl. The PCR reaction was performed with the initial denaturation of 95°C for 5 minutes, followed by 10 cycles of 94°C for 30s, 42°C for 30s and 72°C for 1.5 minutes. Another 30 cycles of 94°C for 30s, 49°C for 30s and 72°C for 1.5 minutes was followed by the final extension of 72°C for 10 minutes. Universal ITS2 primers were used to amplify and sequence the ITS2 region. The PCR was performed in a total volume of 40 μl with 1x Taq Buffer, 1U Taq DNA polymerase, 200 μM dNTPs, 2.0-2.5 mM MgCl2 (depending on template DNA quality), 100pmoles of each primer and 4 μl of DNA. The amplification was done with an initial denaturation temperature of 95°C for 5 minutes followed by 35 cycles of 95°C for 30s, 55°C for 45s and 72°C for 1 minute. The final extension time was set at 72°C for 10 minutes.
The PCR products were sent to M/s Macrogen, South Korea for sequencing. The sequences were edited in Bioedit (v7.1.3) and aligned and deposited in Genbank (18s- Gen Bank: KC791427, KC791428, KC791429 and KC791430; D3- Genbank: JF312867, JF312868 and JF312869; CO I- Genbank: KC791430, KC791431, KC791432, KC791433, KC791434, KC791435, KC791436, KC791437 cytb- Genbank: KF416345- KF416353; ITS2- Genbank: KF416354- KF416362).
Identification of blood meal sources
A modified precipitation test based on the micro – capillary technique was used to determine the source of the blood meal. 50 μl of diluted sandfly blood meal (in 200 μl saline water) was taken into a micro Hematocrit tubes with 50 μl of human or bovine antiserum (Sigma-Aldrich- separately diluted to 1:1000). Formation of a “ring” at the interface of the two phases (blood and antiserum) indicates that the blood meal contains human or bovine tissue. A PCR based assay was performed using universal cytochrome oxidase b primers as described by Ravasan et al. with human and cattle DNA, which are two likely hosts upon which the sandflies have fed, as positive PCR controls.
Detection of Leishmania in sandflies
Sandflies collected by indoor resting collections were dissected and examined for the presence of Leishmania parasites.The DNA of individual flies was extracted using the method described above. A semi nested PCR was performed in Applied Bio systems 9700 thermal cycler with one common forward primer LINR4 5′(GGG GTT GGT GTA AAA TAG GG-3′) and two reverse primers; LIN17 (reverse) (5′-TTT GAA CGG GAT TTC TG-3′), and LIN19 (reverse) (5′-CAG AAC GCC CCT ACC CG-3′) as described by Aransay et al.. The ITS1 region of Leishmania was amplified with the primers LITSR and L5.8S Schonian et al.. A 600 bp region of the heat shock protein gene was amplified using the primers and method described by Garcia et al.. A negative control was performed using DNA extracted from a male fly collected from a non-Leishmania endemic locality. The PCR product was visualized in 1.5% agarose gel stained with ethidium bromide. Sequencing was done by Macrogen, Europe.
The collection was carried out with the approval from the Department of Wild life Conservation, Sri Lanka and health authorities of relevant sampling sites in accordance with all the local rules and regulations regarding collection of sand flies. The study was approved by the Research committee of the Faculty of Graduate Studies, University of Jaffna.
Results and discussion
Sandfly collection and identification
Collection detail of sandflies from different localities
No. of Phlebotomus
Ph. (Eup.) argentipes
The Argentipes complex was earlier described as a species complex with two morpho species namely A and B. Morphospecies B, which has a small sensilla cheatica in the second antennal flagellomere (the ratio of the length of the sensilla cheatica / the length of antennal flagellomere is less than 0.5) was considered as the vector for Leishmania donovani in South India. The taxonomy of this group was reassessed by Ilango, based on the morphometric and meristic characters, in which the species complex was described as a group of three sibling species. All the three members, namely Phlebotomus (Eup.) glaucus, Ph. (Eup.) argentipes sensu stricto and Phlebotomus (Eup.) annandalei, described by Ilango were recorded in the present collection. Ph. (Eup.) glaucus (females with longer sensilla cheatica in the second antennal flagellomere and males with small genital coxite compared to genital style) was found to be the dominant species.
Sequences of the 18SrDNA and D3 regions of the sibling species of Phlebotomus argentipes sensu lato, (with at least 10 individuals sequenced from each sibling species), did not show any sequence variation within the group.
Cytochrome oxidase subunit I and cytochrome b oxidase sequences were, however, variable within the species complex. The amino acid sequences were checked against the reference sequence, which has been reported from other taxa, and all were found to be coding thus ruling out the presence of any nuclear pseudo genes in the sequences.
Ph. (Eup.) mascomai, Ph. (Phb.) papatasi (Abbasi and Warburg Unpublished) and Ph. (Phb.) alexandri were were used as an out group for CO I phylogeny construction. Ph. (Phb.) papatasi, Ph. (Phb.) alexandri, Ph. (Phb.) tobbi (Absravan et al. unpublished) and Ph. (Phb.) sergenti (Vaziri et al., unpublished) were used as out group in cytochrome b tree construction. Ph. (Phb.) papatasi and Ph. (Phb.) salehi were used as out groups in the creation of ITS2 phylogeny.
Uncorrected p distances were calculated for each gene region sequenced using the software MEGA 5.1 with both transversion and transitions estimated as occurring at a uniform rate. The p distance for the ITS2 was 0.42% within each of the two distinct groups within the phylogeny and ranged from 1.3% - 2.5% between these groups. The inter-group p distance for cytochrome oxidase subunit I was 0.01% while the intra-group p distance is 1.2%-1.4%. Cytochrome b oxidase sequences had no intra-group variation but had 1.29% of inter-group species p distance. These genetic distances along with the morphological differences in the sensilla cheatica support the existence of two distinct groups within the Argentipes complex.
On the basis of our molecular data, we propose that the Sri Lankan Phlebotomus (Euphlebotomus) argentipes complex likely consists of only two, rather than three, sibling species. Males and females in the first of these have a length of sensilla cheatica in second antennal flagellomere/ length of second antennal flagellomere ratio of greater than or equal to 0.50 whereas males and females of the second species have a ratio of less than 0.50. With this cutoff value, the morphometrics used in the classification scheme were tested with the data of 200 flies from each sibling species. Other morphometric features used in the classification, such as genital coxite/ genital style ratio, wing overlap (length of R1 over lap over the R2/ length of R2) and wing index (R2/R2+3) were found to be overlapping. The sensilla cheatica ratio was not over lapping with the distribution pattern of 0.33 (minimum) to 0.52 (maximum, observed in only one fly, which might be a slight error in measurement) for sibling species B and with the minimum of 0.53 to the maximum of 0.75 in sibling species A. Generally males possessed shorter sensilla than females.
Another taxonomic feature which could also be used in classification might be the size of the flies given that sibling species B is larger and wider (across the third abdominal segment) than sibling species A. Sibling species A is also darker compared to the sibling species B.
Identification of blood meal sources
Blood meal analysis results for the Phlebotomus (Euphlebotomus) argentipes sibling species A from Hambantota, Sri Lanka (n = 36)
Positive to human antiserum
Positive to bovine antiserum
Negative to both
Positive to both
Leishmania DNA was present in blood fed as well as unfed females (n = 5). The presence of Leishmania donovani DNA in unfed individuals (as determined by the absence of any blood meal in PCR) confirm that the sibling species A, with long sensilla cheatica, is likely to be a vector for leishmaniasis in Sri-Lanka. None of the tested sibling species B (n = 42) female flies had Leishmania DNA within them.
The combined molecular and morphological data indicate that there are likely two sibling species within the Phlebotomus (Euphlebotomus) argentipes complex in Sri Lanka. The sibling species that is proposed not to be a vector for Leishmania transmission in India seems, in contrast, likely to be a vector for cutaneous leishmaniasis and possibly visceral leishmaniasis in Sri Lanka. It will be interesting to study further evolutionary differences between these sibling species and their disease transmission potential.
The authors wish to thank the University Grants Commission (Sri Lanka), the Commonwealth Scholarship Commission, the National Research Council (Sri Lanka; NRC/06/24) and the University of Nottingham (UK) for financial support. Dr. S.D.P. Warusawithana (WHO Country Office, Sri Lanka), Dr. S.A.H. Liyanage (RDHS, Hambantota), the entomological team members; R.G.W. Dissanayake, H.N. Wickramarachchi, P.P.H. Sumith Danel and J. Migelhewa, Demographic survey staff; C. Weerasooriya, E. L.Lionel, E.T. Sirisena, A.P. Suneth, S. Liyanage, H.U. Nimal, and J.N. Ratnayake, all the entomological team members from Galle and Hambantota and staffs of the regional office of Anti Malaria campaign of Hambantota are greatly acknowledged for monetary (field work) and field support. We acknowledge Prof. D.J. Weligama from the Department of Parasitology of the Faculty of Medicine and Allied Sciences, the University of Rajarata, Sri Lanka. The Department of Wildlife Conservation (Sri Lanka) is acknowledged for the permission to collect sandflies.
- Ashford RW: The Leishmaniasis as emerging and re-emerging globally. Int J Parasitol. 2000, 30: 1269-1281.View ArticlePubMedGoogle Scholar
- Ilango K: Spatial and temporal distribution of Phlebotomine Sand flies in relation to Leishmaniasis in India. Entomology (Ecology & Biodiversity). Edited by: Tyagi BK, Vijay Veer . 2011, India: Scientific Publisher, 99-103. 1Google Scholar
- Lysenko AJ: Distribution of leishmaniasis in the Old World. Bull W H O. 1971, 44: 515-520.PubMed CentralPubMedGoogle Scholar
- de Wijesundara M: Cutaneous leishmaniasis: an emerging health risk in Sri Lanka. Ceylon Med J. 2001, 46: 151-152.Google Scholar
- Thakur CP: Leishmaniasis research- the challenges ahead. Indian J Med Res. 2006, 123: 193-194.PubMedGoogle Scholar
- Athukorale DN, Seneviratne JKK, Ihalamulla RL, Premaratne UN: Locally acquired cutaneous leishmaniasis in Sri Lanka. J Trop Med Hyg. 1992, 95: 432-433.PubMedGoogle Scholar
- Sirwardana HVYD, Noyes HA, Beeching NJ, Chance ML, Karunaweera ND, Bates PA: Leishmania donovani and cutaneous leishmaniasis in Sri Lanka. J Emerg Infect Dis. 2007, 13: 476-478.View ArticleGoogle Scholar
- Ready PD: Leishmaniasis emerging and climate change. Rev Sci Tech Off Epiz. 2008, 27: 399-412.Google Scholar
- Ilango K: Morphological characteristics of the antennal flagellum and its sensilla cheatica with character displacement in the sandfly Phlebotomus argentipes Annandale and Brunetti sensu lato (Diptera: Psychodidae). J Biosci. 2000, 25: 163-172.PubMedGoogle Scholar
- Frizzi G: Salivary gland chromosomes of Anopheles. Nature. 1947, 160: 226-227.View ArticleGoogle Scholar
- Colluzi M, Sabatini A: Cytogenetics observation on species C of the Anopheles gambiae complex. Parasitologia. 1968, 10: 155-156.Google Scholar
- Collins FH, Paskewitz SM: A review of the use of ribosomal DNA to differentiate among cryptic Anopheles species. Insect Mol Biol. 1996, 05: 1-9.View ArticleGoogle Scholar
- Green CA: Malaria epidemiology and anopheline cytogenetics. Cytogen Gen Vect. Edited by: Pal R, Kitzmiller JB, Kanda T. 1981, Amsterdam: Elsevier Biomedical PressGoogle Scholar
- Subbarao SK, Vasantha K, Sharma VP: Responses of Anopheles culicifacies sibling species A and B to DDT and HCH in India: implications in malaria control. Med Vet Entomol. 1988, 2: 219-223.View ArticlePubMedGoogle Scholar
- Khalid NM, Aboud MA, Alrabba FM, Elnaiem D-EA, Tripet F: Evidence for genetic differentiation at the microgeographic scale in Phlebotomus papatasi populations from Sudan. Parasit Vectors. 2012, 5: 249-doi:10.1186/1756-3305-5-249PubMed CentralView ArticlePubMedGoogle Scholar
- Lewis DJ: The phlebotamine sandflies (Diptera, Psychodidae) of the Oriental region. Bull Brit Mus (Nat Hist). 1978, 37: 218-331.Google Scholar
- Ilango K: A taxonomic reassessment of the Phlebotomus argentipes species complex (Diptera:Psychodidae:Phlebotaminae). J Med Entomol. 2010, 47: 1-15.View ArticlePubMedGoogle Scholar
- Litvaitis MK, Nunn G, Thomas WK, Kocher TD: A molecular approach for the identification of Meiofaunal turbellarians (Platyhelminthes, Turbellaria). Mar Biol. 1994, 120: 437-442.View ArticleGoogle Scholar
- Surendran SN, Karunaratne SH, Adams Z, Hemingway J, Hawkes NJ: Molecular and biochemical characterization of a sand fly population from Sri Lanka: evidence for insecticide resistance due to altered esterases and insensitive acetylcholinesterase. Bull Entomol Res. 2005, 95: 371-380.View ArticlePubMedGoogle Scholar
- Simon C, Frati F, Beckenbach A, Crespi B, Liu H, Flook P: Evolution, weighting and phylogenetic utility of mitochondrial gene sequences and a compilation of conserved polymerase chain reaction primers. Ann Entomol Soc Am. 1994, 87: 651-701.View ArticleGoogle Scholar
- Esseghir S, Ready PD, Killick-Kendrick R, Ben-Ismail R: Mitochondrial haplotypes and phylogeography of Phlebotomus vectors of Leishmania major. Insect Mol Biol. 1997, 6: 211-225.View ArticlePubMedGoogle Scholar
- Porter CH, Collins FH: Species diagnostic differences in a ribosomal DNA I nternal transcribed spacer from the sibling species Anopheles freeborni and An. hermsi (Diptera: Culicidae). Am J Trop Med Hyg. 1991, 45: 271-279.PubMedGoogle Scholar
- Tempelis C, Lofy M: A modified Precipitin Method for Identification of Mosquito Blood Meals. Am J Trop Med Hyg. 1963, 12: 825-831.PubMedGoogle Scholar
- Ravasan NM, Oshaghi MA, Javadian E, Rassi Y, Sadraei J, Mohtarami F: Blood meal identification in field-captured sand flies: comparison of PCR-RFLP and ELISA assay. Iran J Arthropod Borne Dis. 2009, 3: 8-18.Google Scholar
- Aransay MA, Scoulica E, Tselentis Y: Detection and Identification of Leishmania DNA within Naturally Infected Sand Flies by Semi nested PCR on Minicircle Kinetoplastic DNA. Appl Environ Microbiol. 2000, 66: 1933-1938.PubMed CentralView ArticlePubMedGoogle Scholar
- Schonian G, Nasereddin A, Dinse N, Schweynoch C, Schallig HD, Presber W, Jaffe CL: PCR diagnosis and characterization of Leishmania in local and imported clinical samples. Diagn Microbiol Infe Dis. 2003, 47: 349-358.View ArticleGoogle Scholar
- Garcia L, Kindt A, Bermudez H, Llanos-Cuentas A, Dedoncker S, Arevalo J, Wilber Quispe Tintaya K, Dujardin JC: Culture-independent species typing of Neotropical Leishmania for clinical validation of a PCR-based assay targeting heat shock protein 70 genes. J Clin Microbiol. 2004, 42: 2294-2297.PubMed CentralView ArticlePubMedGoogle Scholar
- Lunt DH, Zhang DX, Szymura JM, Hewitt GM: The insect cytochrome oxidase I gene: evolutionary patterns and conserved primers for phylogenetic studies. Insect Mol Biol. 1996, 5: 153-165.View ArticlePubMedGoogle Scholar
- Guindon S, Gascuel O: A simple, fast and accurate algorithm to estimate large phylogenies by maximum likelihood. Syst Biol. 2008, 52: 696-704.View ArticleGoogle Scholar
- Depaquit J, Muller F, Leger N: Phlebotomus (Euphlebotomus) barguesae n. sp. from Thailand (Diptera - Psychodidae). Parasit Vectors. 2009, 2: 5-PubMed CentralView ArticlePubMedGoogle Scholar
- Boudabous R, Bounamous A, Jouet D, Depaquit J, Augot D, Ferte H, Berchi S, Couloux A, Veuille M, Babba H: Mitochondrial DNA Differentiation Between Two Closely Related Species, Phlebotomus (Paraphlebotomus) chabaudi and Phlebotomus (Paraphlebotomus) riouxi (Diptera: Psychodidae), Based on Direct Sequencing and Polymerase Chain Reaction-Restriction Fragment Length Polymorphism. Ann Entomol Soc Am. 2009, 102: 347-353.View ArticleGoogle Scholar
- Hamarsheh O, Presber W, Abdeen Z, Sawalha S, Al-Lahem A, Schonian G: Genetic structure of Mediterranean populations of the sandfly Phlebotomus papatasi by mitochondrial cytochrome b haplotype analysis. Med Vet Entomol. 2007, 21: 270-277.View ArticlePubMedGoogle Scholar
- Kruger A, Struven L, Post RJ, Faulde M: The sandflies (Diptera: Psychodidae, Phlebotominae) in military camps in northern Afghanistan (2007–2009), as identified by morphology and DNA 'barcoding’. Ann Trop Med Parasit. 2011, 105: 163-176.PubMed CentralView ArticlePubMedGoogle Scholar
- Hamarsheh O, Presbe P, Abdeen Z, Sawalha S, Al-Lahem A, Schoenian G: Sequence analysis of the ribosomal DNA ITS2 region for Phlebotomus papatasi (Diptera: Psychodidae). J Entomol. 2007, 4: 167-180.View ArticleGoogle Scholar
- Depaquit J, Lienard E, Verzeaux-Griffon A, Ferte H, Bounamous A, Gantier JC, Hanafi HA, Jacobson RL, MaroliM M-VV, Muller F, Ozbel Y, Svobodova M, Volf P, Leger N: Molecular homogeneity in diverse geographical populations of Phlebotomus papatasi (Diptera, Psychodidae) inferred from ND4 mtDNA and ITS2 rDNA Epidemiological consequences. Infect Genet Evol. 2008, 8: 159-170.View ArticlePubMedGoogle Scholar
- Tamura K, Peterson D, Peterson N, Stecher G, Nei M, Kumar S: MEGA5: molecular evolutionary genetics analysis using maximum likelihood, evolutionary distance, and maximum parsimony method. Mol Biol Evol. 2011, 28: 2731-2739.PubMed CentralView ArticlePubMedGoogle Scholar
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.