Phylogenetic, phylogeographic and divergence time analysis of Anopheles subpictus species complex using ITS2 and COI sequences

Abstract

Objective: To address the phylogenetic and phylogeographic relationship between different lineages of Anopheles (An.) subpictus species complex in most parts of the Asian continent by maximum utilization of Internal Transcriber Spacer 2 (ITS2) and cytochrome C oxidase I (COI) sequences deposited at the GenBank.

Methods: Seventy-five ITS2, 210 COI and 26 concatenated sequences available in the NCBI database were used. Phylogenetic analysis was performed using Bayesian likelihood trees, whereas median-joining haplotype networks and time-scale divergence trees were generated for phylogeographic analysis. Genetic diversity indices and genetic differentiation were also calculated.

Results: Two genetically divergent molecular forms of An. subpictus species complex corresponding to sibling species A and B are established. Species A evolved around 37-82 million years ago in Sri Lanka, India, and the Netherlands, and species B evolved around 22-79 million years ago in Sri Lanka, India, and Myanmar. Vietnam, Thailand, and Cambodia have two molecular forms: one is phylogenetically similar to species B. Other forms differ from species A and B and evolved recently in the above mentioned countries, Indonesia and the Philippines. Genetic subdivision among Sri Lanka, India, and the Netherlands is almost absent. A substantial genetic differentiation was obtained for some populations due to isolation by large geographical distances. Genetic diversity indices reveal the presence of a long-established stable mosquito population, at mutation-drift equilibrium, regardless of population fluctuations.

Conclusions: An. subpictus species complex consists of more than two genetically divergent molecular forms. Species A is highly divergent from the rest. Sri Lanka and India contain only species A and B.

Keywords

Molecular systematics / ITS2 / COI / DNA sequences / Phylogeny / Phylogeography

References

[1] Manguin S, Boëte C. Global impact of mosquito biodiversity, human vector-borne diseases and environmental change. In: The importance of biological interactions in the study of biodiversity. London: IntechOpen; 2011, p. 27-50.

[2] Owino EA. Why Kenya should worry about Anopheles stephensi. Asian Pac J Trop Med 2023; 16(3): 99-101.

[3] Karunaweera ND, Galappaththy GN, Wirth DF. On the road to eliminate malaria in Sri Lanka: Lessons from history, challenges, gaps in knowledge and research needs. Malar J 2014; 13(59): 1-10.

[4] Harischandra IN, Dassanayake RS, De Silva BGDNK. Mitochondrial-COII sequence polymorphism reflects spatial genetic clustering of Anopheles culicifacies sibling species E in Sri Lanka. J Vector Borne Dis 2020; 57(2): 139-146.

[5] Venkatesan P. The 2023 WHO world malaria report. Lancet Microbe 2024; 5(3): e214.

[6] Simac JN, Badar S, Farber J, Brako MO, Giudice-Jimenez RL, Raspa S, et al. Malaria elimination in Sri Lanka. J Health Spec 2017; 5(2): 60-65.

[7] Ibáñez-Justicia A, Smitz N, Den Hartog W, Van de Vossenberg B, De Wolf K, Deblauwe I, et al. Detection of exotic mosquito species (Diptera: Culicidae) at international airports in Europe. Int J Environ Res Public Health 2020; 17(10): 3450.

[8] Jayatunga DPW, Harischandra IN, Chandrasekharan NV, De Silva BGDNK. Compensatory base changes reveal sexual incompatibility among members of the Anopheles subpictus Sensu Lato (Diptera: Culicidae) species complex in Sri Lanka. Life 2021; 11(3): 211.

[9] Gunathilaka PADHN. Annotated checklist and review of the mosquito species (Diptera: Culicidae) in Sri Lanka. J Insect Biodivers 2018; 7(3): 38-50.

[10] Rueda LM, Debboun M. Taxonomy, identification, and biology of mosquitoes. In: Mosquitoes, communities, and public health in Texas. New York: Academic Press; 2019, p. 3-7.

[11] Chhilar JS, Chaudhry S. Phylogenetic analysis of Anopheles (Cellia) subpictus Grassi using rDNA-ITS2 sequence. Proc Zool Soc 2012; 65(1): 1-10.

[12] Surendran S, Ramasamy R. The Anopheles culicifacies and Anopheles subpictus species complexes in Sri Lanka and their implications for malaria control in the country. Trop Med Health 2010; 38(1): 1-11.

[13] Sindhania A, Das MK, Sharma G, Surendran SN, Kaushal BR, Lohani HP, et al. Molecular forms of Anopheles subpictus and Anopheles sundaicus in the Indian subcontinent. Malar J 2020; 19(417): 1-17.

[14] Suguna SK, Rathinam KG, Rajavel AR, Da VDH. Morphological and chromosomal descriptions of new species in the Anopheles subpictus complex. Med Vet Entomol 1994; 8(1): 88-94.

[15] Subbarao SK. Anopheline species complexes in South-East Asia. New Delhi: WHO Regional Office for South-East Aisa; 1998, p. 1-69.

[16] Jude PJ, Ramasamy R, Surendran SN. Bionomic aspects of the Anopheles subpictus species complex in Sri Lanka. J Insect Sci 2014; 14(1): 97.

[17] Singh SP. Morphotaxonomical studies to identify the member of the Anopheles subpictus Grassi (Diptera: Culicidae) species complex in villages of District Mewat Haryana State, India. Int Interdiscip Res J 2014; 4(1): 239-244.

[18] Abhayawardana TA, Wijesuriya SR, Dilrukshi RK. Anopheles subpictus complex: Distribution of sibling species in Sri Lanka. Indian J Malariol 1996; 33(2): 53-60.

[19] Cooper RD, Edstein MD, Frances SP, Beebe NW. Malaria vectors of Timor-Leste. Malar J 2010; 9(40): 1-11.

[20] Surendran SN, Singh OP, Jude PJ, Ramasamy R. Genetic evidence for malaria vectors of the Anopheles sundaicus complex in Sri Lanka with morphological characteristics attributed to Anopheles subpictus species B. Malar J 2010; 9(343): 1-9.

[21] Sivabalakrishnan K, Kanapathy G, Jeyadas TTP, Surendran SN. Population genetic structure of Anopheles subpictus species B using COII and Cytb markers. J JSA 2020; 2(1): 38-47.

[22] Jayatunga DPW, Harischandra IN, Chandrasekharan NV, De Silva NK. Alterations and interchange of morphometric characters in different life cycle stages with reference to genomic variations of Anopheles subpictus (Diptera; Culicidae) sibling species complex in Sri Lanka. Insects 2018; 9(3): 89.

[23] Hill SM, Crampton JM. DNA-based methods for the identification of insect vectors. Ann Trop Med Parasitol 1994; 88(3): 227-250.

[24] Sallum MAM, Schultz TR, Foster PG, Aronstein K, Wirtz RA, Wilkerson RC. Phylogeny of Anophelinae (Diptera: Culicidae) based on nuclear ribosomal and mitochondrial DNA sequences. Syst Entomol 2002; 27(3): 361-382.

[25] Zomuanpuii R, Ringngheti L, Brindha S, Gurusubramanian G, Kumar NS. ITS2 characterization and Anopheles species identification of the subgenus Cellia. Acta Trop 2013; 125(3): 309-319.

[26] Murugan K, Vadivalagan C, Karthika P, Panneerselvam C, Paulpandi M, Subramaniam J, et al. DNA barcoding and molecular evolution of mosquito vectors of medical and veterinary importance. Parasitol Res 2016; 115: 107-121.

[27] Khoshdel-nezamiha F, Vatandoost H, Oshaghi MA, Azari-hamidian S, Mianroodi RA, Dabiri F, et al. Molecular characterization of mosquitoes (Diptera: Culicidae) in Northwestern Iran by using rDNA-ITS2. Jpn J Infect Dis 2016; 69(4): 319-322.

[28] Weeraratne TC, Surendran SN, Reimer LJ, Wondji CS, Perera MD, Walton C, et al. Molecular characterization of Anopheline (Diptera: Culicidae) mosquitoes from eight geographical locations of Sri Lanka. Malar J 2017; 16(234): 1-14.

[29] Singh D, Rai MK, Vashist D. Molecular studies on Indian mosquitoes (Diptera: Culicidae): A review. Int J Entomol Res 2018; 3(2): 36-103.

[30] Surendran SN, Sarma DK, Jude PJ, Kemppainen P, Kanthakumaran N, Gajapathy K, et al. Molecular characterization and identification of members of the Anopheles subpictus complex in Sri Lanka. Malar J 2013; 12(304): 1-11.

[31] Weeraratne TC, Surendran SN, Walton C, Karunaratne SHPP. Genetic diversity and population structure of malaria vector mosquitoes: Anopheles subpictus, Anopheles peditaeniatus, and Anopheles vagus in five districts of Sri Lanka. Malar J 2018; 17(271): 1-8.

[32] Coyne JA. The evolutionary genetics of speciation. Philos Trans R Soc B Biol Sci 1998; 353(1366): 287-305.

[33] Provine WB. Ernst Mayr: Genetics and speciation. Genetics 2004; 167(3): 1041-1046.

[34] Mayr E. Origin and history of some terms in systematic and evolutionary biology. Syst Zool 1978; 27(1): 83-88.

[35] Wilai P, Ali RS, Saingamsook J, Saeung A, Junkum A, Walton C, et al. Integrated systematics of Anopheles subpictus (Diptera: Culicidae) in the Oriental Region, with emphasis on forms in Thailand and Sulawesi, Indonesia. Acta Trop 2020; 208: 105503.

[36] Karimian F, Oshaghi MA, Sedaghat MM, Waterhouse RM, Vatandoost H, Hanafi-Boid AA, et al. Phylogenetic analysis of the Oriental-Palearctic-Afrotropical members of Anopheles (Culicidae: Diptera) based on nuclear rDNA and mitochondrial DNA characteristics. Jpn J Infect Dis 2014; 67(5): 361-367.

[37] Moreno M, Marinotti O, Krzywinski J, Tadei WP, James AA, Achee NL, et al. Complete mtDNA genomes of Anopheles darlingi and an approach to anopheline divergence time. Malar J 2010; 9(1): 1-13.

[38] Lorenz C, Alves JMP, Foster PG, Suesdek L, Sallum MAM. Phylogeny and temporal diversification of mosquitoes (Diptera: Culicidae) with an emphasis on the Neotropical fauna. Syst Entomol 2021; 46(4): 798-811.

[39] Depaquit J, Lienard E, Verzeaux-Griffon A, Ferte H, Bounamous A, Gantier JC, et al. 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(2): 159-170.

[40] Rangel-Gamboa L, Martinez-Hernandez F, Maravilla P, Arenas-Guzman R, Flisser A. Update of phylogenetic and genetic diversity of Sporothrix schenckii sensu lato. Med Mycol J 2015; 54(3): 248-255.

[41] Martinez-Villegas L, Assis-Geraldo J, Koerich LB, Collier TC, Lee Y, Main BJ, et al. Characterization of the complete mitogenome of Anopheles aquasalis, and phylogenetic divergences among Anopheles from diverse geographic zones. PLoS One 2019; 14(9): 1-22.

[42] Reidenbach KR, Cook S, Bertone MA, Harbach RE, Wiegmann BM, Besansky NJ. Phylogenetic analysis and temporal diversification of mosquitoes (Diptera: Culicidae) based on nuclear genes and morphology. BMC Evol Biol 2009; 9(298): 1-14.

[43] Harbach RE. The classification of genus Anopheles (Diptera: Culicidae): A working hypothesis of phylogenetic relationships. Bull Entomol Res 2004; 94(6): 537-553.

[44] Singh D, Vashist D. DNA barcoding as a tool supporting species identification and molecular phylogeny for Indian Anopheles (Diptera: Culicidae: Anophelinae). J Entomol Res 2020; 44(3): 409-416.

[45] Mayoke A, Muya SM, Bateta R, Mireji PO, Okoth SO, Onyoyo SG, et al. Genetic diversity and phylogenetic relationships of tsetse flies of the palpalis group in Congo Brazzaville based on mitochondrial cox1 gene sequences. Parasit Vectors 2020; 13(253): 1-16.

[46] Dharmarathne HA, Weerasena OV, Perera KL, Galhena G. Genetic characterization of Aedes aegypti (Diptera: Culicidae) in Sri Lanka based on COI gene. J Vector Borne Dis 2020; 57(2): 157.

[47] Dogan I, Dogan N. Statistical tests for neutrality: Review. Turkiye Klin J Biostat 2017; 9(2): 167-174.

[48] Ebrahimi S, Bordbar A, Parvizi P. Genetic dynamics in the sand fly (Diptera: Psychodidae) nuclear and mitochondrial genotypes: Evidence for vector adaptation at the border of Iran with Iraq. Parasit Vectors 2016; 9(319): 1-13.

[49] Jacquot MP. Genomic diversity of pathogenic bacteria of the Borrelia burgdorferi species complex: Evolution and molecular epdemiology. Doctoral dissertation. Universite Blaise Pascal-Clermont-Ferrad II University; 2014.

[50] Jensen JL, Bohonak AJ, Kelley ST. Isolation by distance, web service. BMC Genet 2005; 6(1): 1-6.

[51] Freitas LA, Russo CAM, Voloch CM, Mutaquiha OCF, Marques LP, Schrago CG. Diversification of the genus Anopheles and a neotropical clade from the Late Cretaceous. PLoS One 2015; 10(8): 1-12.