Integrated systematics of Anopheles subpictus (Diptera: Culicidae) in the Oriental Region, with emphasis on forms in Thailand and Sulawesi, Indonesia
Graphical abstract
Introduction
Anopheles subpictus Grassi of the Pyretophorus Series of subgenus Cellia Theobald was first described from a single specimen collected in India (Christophers, 1933). It has a very wide distribution, ranging eastward from Saudi Arabia and Iran across southern Asia (Pakistan, India, Sri Lanka, Bangladesh, Myanmar, Thailand), along southern coastal areas of Cambodia and Vietnam, and coastal areas of Malaysia, Indonesia and Papua New Guinea, and as far east as the Solomon Islands (Sinka et al., 2011). The taxon is a complex of species, the Subpictus Complex, consisting of four species in India, which are informally designated as species A, B, C and D (Suguna et al., 1994). However, current evidence is insufficient to determine which of the four species is conspecific with the nominotypical form. Identification of members of the complex has been based primarily on inversions in the X-arm of polytene chromosomes and morphological characters that include the number of float ridges of eggs, the number of branches of seta 4-M of larvae and seta 7-I of pupae, and pale and dark bands on the maxillary palpi of females. According to Suguna et al., species A, C and D generally occur in freshwater habitats of inland areas, whereas species B is restricted to coastal brackish water habitats. However, recent studies in Sri Lanka revealed that identification based on polytene chromosome inversions and morphological characteristics of the sibling species in India do not correspond to molecular identifications based on ITS2 and COI sequences, by which only two linages corresponding to species A and B were detected in the country, and species B is not restricted to coastal areas (Abhayawardana et al., 1996; Surendran et al., 2013). Similarly, morphological assessment and molecular identification in India using ITS2 and COI sequences (Kumar et al., 2007; Kaura et al., 2010; Chhilar and Chaudhry, 2012) suggested the presence of only species A and B. Additionally, sequences for species C and D are not present in GenBank. Therefore, further investigation is needed to determine whether differences in polytene chromosomes and morphological characters described for species C and D may be variations of either species A or B.
Phylogenetic analyses of the Subpictus Complex have revealed that species A and B do not form a monophyletic clade, but share genetic similarity with closely related species of the Pyretophorus Series: species A is closely related to An. vagus Dönitz whereas species B is closely related to An. epiroticus Linton & Harbach and An. sundaicus (Rodenwaldt) of the Sundaicus Complex (Surendran et al., 2013). The distributions of species A and B are well documented in India and Sri Lanka, but little is known about these forms elsewhere. Species A appears to be confined to India and Sri Lanka whereas species B may have a wider distribution (Surendran et al., 2010, 2013). In general, An. subpictus is primarily zoophilic and exophagic, but indices of anthropophilic behavior may be high (up to 20% human feeding) in some areas (reviewed by Chandra et al., 2010). The role the taxon plays in the transmission of human malarial parasites is well documented in coastal areas of India and Sri Lanka (Abhayawardana et al., 1996; Panicker et al., 1981). In the Indonesian Archipelago, An. subpictus is widely distributed from Sumatra to the Maluku Islands, but the role it plays in the transmission of malarial parasites is only well documented in Java, Sulawesi, the Lesser Sundas and Lombok Islands (reviewed by Elyazar et al., 2013). Whether differences in vectorial capacity reflect genetic distinctions is not known and requires further study.
In Thailand, An. subpictus is distributed throughout the country, including coastal areas in southern provinces, but it is not regarded as a vector of malarial plasmodia due to its zoophilic behavior (Rattanarithikul et al., 2006). The specific identity of An. subpictus in Thailand and elsewhere in Southeast Asia is not clear. Baimai et al. (1996) reported the finding of three karyotypic forms of An. subpictus in Thailand and two forms in Indonesia and the Philippines, but whether these forms reflect interspecific difference is not known. In the present study, we conducted a morphological and molecular study of specimens of An. subpictus collected in inland and coastal areas of Thailand, and in South Sulawesi, Indonesia, and compared ITS2 and COI sequences to sequences available in GenBank for specimens from other Asian countries in order to determine the genetic diversity of the taxon. We also conducted crossing experiments to evaluate the genetic compatibility of three COI clades that were found in specimens collected in Thailand.
Section snippets
Mosquito collections and identification
Adult mosquitoes were collected in cattle sheds by aspirator in Doi Saket, Hang Dong, Mae On, Mueang and San Pa Tong Districts of Chiang Mai Province of Thailand (site 1, Fig. 1 and Table 1), and South Sulawesi Province of Indonesia (site 9). Mosquito larvae were collected by dipping mainly in rice fields in Kamphaengsaen District of Nakhon Pathom Province (site 2) and Tha Yang District of Phetchaburi Province (site 3) in Thailand. Specimens previously collected in Ratchaburi (site 4) and Phang
Phylogenetic analyses
A total of 78 wild-caught specimens from Chiang Mai (29), Nakhon Pathom (1) and Phetchaburi (43) Provinces of Thailand, and South Sulawesi of Indonesia (5) were sequenced for fragments of the ITS2 region and COI gene. The length of the ITS2 and COI sequences was 564 and 422 bp, respectively. The best-fit models of nucleotide substitution for ML analysis were HKY+G for ITS2 and GTR+I + G for COI sequences.
Phylogenetic analysis of unambiguous ITS2 sequences revealed four distinct clades (Clades
Discussion
Morphologically, An. subpictus in Thailand exhibits a mixture of the characters found in species A and B, but differs from species C and D described by Suguna et al., 1994. In Thai females, the apical pale band of the maxillary palpi is longer than the subapical dark band, which is characteristic of species A. In species B, the apical pale band is shorter than the subapical dark band, and in species C and D the dark and pale bands are equal in length. Seta 7-I of pupae has 4–7 branches in Thai
CRediT authorship contribution statement
Parinya Wilai: Investigation, Methodology, Formal analysis. Rusdiyah Sudirman Made Ali: Resources. Jassada Saingamsook: Data curation. Atiporn Saeung: Supervision. Anuluck Junkum: Supervision. Catherine Walton: Supervision, Writing - review & editing. Ralph E Harbach: Writing - review & editing. Pradya Somboon: Conceptualization, Formal analysis, Funding acquisition, Project administration, Supervision, Writing - original draft, Writing - review & editing.
Declaration of Competing Interest
The authors declare no conflict of interest.
Acknowledgements
We thank the Faculty of Medicine, Chiang Mai University for providing a scholarship to Parinya Wilai to support her PhD study. This research was partially supported by Faculty of Medicine (Grant No. PAR-2563-07268) and Office of Research Administration, Chiang Mai University and the Bio & Medical Technology Development Program of the National Research Foundation (NRF) funded by the Korean government (MSIT) (No. 2015M3A9B6073666).
References (28)
- et al.
Molecular identification of mosquitoes of the Anopheles maculatus group of subgenus Cellia (Diptera: Culicidae) in the Indonesian Archipelago
Acta Trop.
(2019) - et al.
A review on Anopheles subpictus Grassi – a biological vector
Acta Trop.
(2010) - et al.
The distribution and bionomics of Anopheles malaria vector mosquitoes in Indonesia. Adv
Parasitol.
(2013) - et al.
Molecular identification of mosquito species
Biol. J. Linn. Soc.
(1999) - et al.
Anopheles subpictus complex: distribution of sibling species in Sri Lanka
Indian J. Malariol.
(1996) - et al.
Metaphase karyotypes of Anopheles of Thailand and Southeast Asia. VI. The Pyretophorus and the Neomyzomyia Series
Subgenus Cellia (Diptera: Culicidae). J. Am. Mosq. Control Assoc.
(1996) - et al.
Discrimination of all members of the Anopheles punctulatus complex by polymerase chain reaction-restriction fragment length polymorphism analysis
Am. J. Trop. Med. Hyg.
(1995) - et al.
Phylogenetic analysis of Anopheles (Cellia) subpictus Grassi using rDNA-ITS2 sequence
Proc. Zool. Soc.
(2012) The fauna of British India, including Ceylon and Burma. Diptera
(1933)- et al.
jModelTest 2: more models, new heuristics and parallel computing
Nat. Methods
(2012)