Tests of conspecificity for allopatric vectors: Simulium nodosum and Simulium shirakii (Diptera: Simuliidae) in Asia
© Low et al. 2015
Received: 7 April 2015
Accepted: 20 May 2015
Published: 29 May 2015
Allopatric populations present challenges for biologists working with vectors. We suggest that conspecificity can be concluded in these cases when data from four character sets—chromosomal, ecological, molecular, and morphological—express variation no greater between the allopatric populations than between corresponding sympatric populations. We use this approach to test the conspecificity of Simulium nodosum Puri on the mainland of Southeast Asia and Simulium shirakii Kono & Takahasi in Taiwan. The validity of these two putative species has long been disputed given that they are morphologically indistinguishable.
The mitochondria-encoded cytochrome c oxidase subunit I (COI), 12S rRNA, and 16S rRNA genes and the nuclear-encoded 28S rRNA gene support the conspecific status of S. nodosum from Myanmar, Thailand, and Vietnam and S. shirakii from Taiwan; 0 to 0.19 % genetic differences between the two taxa suggest intraspecific polymorphism. The banding patterns of the polytene chromosomes of the insular Taiwanese population of S. shirakii and mainland populations of S. nodosum are congruent. The overlapping ranges of habitat characteristics and hosts of S. nodosum and S. shirakii corroborate the chromosomal, molecular, and morphological data.
Four independent sources of evidence (chromosomes, DNA, ecology, and morphology) support the conspecificity of S. nodosum and S. shirakii. We, therefore, synonymize S. shirakii with S. nodosum. This study provides a guide for applying the procedure of testing conspecificity to other sets of allopatric vectors.
KeywordsSimulium nobile species group Black flies Phylogeny Polytene chromosomes Systematics Multi-locus Conspecificity
Simuliid vectors of human and animal pathogens typically are complexes of evolutionarily and ecologically distinct cryptic species, for which accurate identification is an essential first step in the epidemiological understanding and control of arthropod-borne diseases [1–4]. Compounding the challenge of recognizing cryptic species is the difficulty of evaluating the species status of allopatric populations. The problem is particularly acute when populations are widely disjunct on the mainland, between mainlands, or between an island and the mainland . Assigning the same or different names to allopatric populations has consequences for understanding vector potential and developing control strategies .
A classic example of the challenge of evaluating species status for allopatric populations involves Simulium nodosum Puri, distributed from India across southern China, Myanmar, Thailand, and Vietnam , and the morphologically similar Simulium shirakii Kono & Takahasi from Taiwan, more than 130 km from the Chinese mainland. Although S. shirakii has been treated for 80 years as a separate species, its species status has long been questioned because of morphological similarity with S. nodosum . Simulium nodosum is mammalophilic, attacking humans and bovids [9–11]. It is also a vector of the agents of filariasis to ruminants .
Although allopatry confounds the evaluation of reproductive isolation, we suggest that allopatric populations, including putative species, can be considered conspecific when differences in their molecular sequences, chromosomal profiles, morphology, and ecology are within the range of variation for a geographically cohesive, panmictic set of populations. All four character sources require evaluation, given the existence of homosequential sibling species  and the failure of up to three of the character sets to discriminate reproductively isolated sympatric species .
We applied this test to two members of the Simulium nobile species group, S. nodosum and S. shirakii. The detailed morphological comparisons by Takaoka & Suzuki  between S. nodosum and S. shirakii originally suggested conspecificity and prompted the present study. We used a multi-locus phylogenetic analysis of one nuclear and three mitochondrial genes, a comparative analysis of the polytene chromosomes, and an evaluation of ecological data associated with our collections and in the literature [9, 14–16].
No national permissions were required for this study, which did not involve endangered or protected species. No specific permissions were required to access the study sites; the collections were made on public lands.
Collection data for populations used in chromosomal and molecular analyses of Simulium nodosum and S. shirakii
San Village, Chiang-Tung, Myanmar
Da Chais, Lam Dong, Vietnam
Doi Saket, Chiang Mai, Thailand
Guangfu, Hualien, Taiwan
Ruisui, Hualien, Taiwan
The nucleotide sequences of the mitochondria-encoded COI, 12S rRNA, and 16S rRNA genes and the nuclear-encoded 28S rRNA gene were used. These genetic markers have been used to differentiate other simuliid species [20–23]. In addition to S. nodosum and S. shirakii, we included two nominal members of the S. nobile species group—Simulium nobile De Meijere from Gombak, Selangor, Malaysia (collected 23/07/14) and Simulium timorense Takaoka, Hadi & Sigit from Kupang, Timor Island, Indonesia (collected 27/02/14)—for phylogenetic analysis.
Genomic DNA was extracted from each of five specimens per location, using the i-genomic CTB DNA Extraction Mini Kit (iNtRON Biotechnology Inc., Seongnam, South Korea). Amplifications of the mitochondria-encoded COI, 12S rRNA, and 16S rRNA, and nuclear-encoded 28S rRNA genes were undertaken in a final volume of 50 μL containing 50–100 ng genomic DNA, 25 μL of ExPrime Taq Master Mix (GENETBIO Inc., Daejeon, South Korea), and 10 pmol of each forward and reverse primer. The primers used in this study were adopted from Folmer et al.  for COI, Kocher et al.  and Simon et al.  for 12S rRNA, Xiong & Kocher  for 16S rRNA, and Low et al.  for 28S rRNA.
Ranges of intra- and interspecific genetic distances (uncorrected p, expressed as percentages) based on concatenated COI, 12S rRNA, 16S rRNA, and 28S rRNA sequences
1. S. nodosum (Myanmar)
2. S. nodosum (Vietnam)
3. S. nodosum (Thailand)
4. S. shirakii (Taiwan)
5. S. nobile
6. S. timorense
We used the Feulgen technique and chromosomal slide-mounting procedures outlined by Adler et al.  to prepare the polytene chromosomes of 51 larvae from two sites (19 larvae from Guangfu, 32 larvae from Ruisui) in Taiwan (Table 1). Larval carcasses were deposited in the Clemson University Arthropod Collection. The chromosomes of all 51 larvae (25 females, 26 males) were compared band-for-band with the standard map for the subgenus Simulium [33, 34] and with the chromosomes of larvae analyzed by Tangkawanit et al.  from nine sites in Thailand.
The phylogenetic tree for the concatenated dataset comprised two main clusters (Fig. 1). One cluster, consisting of S. nobile and S. timorense, was supported with full posterior probability/bootstrap values (BI = 1.00, NJ = 100 %, MP = 100 %, ML = 100 %). The other cluster, consisting of S. nodosum and S. shirakii from different geographical regions, was supported with high to full posterior probability/bootstrap values (BI = 0.99, NJ = 100 %, MP = 100 %, ML = 100 %).
No phylogenetic tree from the concatenated dataset nor the single-locus analyses separated S. nodosum from S. shirakii. Both putative species were distributed randomly in the trees (Additional file 1: Figures S1-S4). The pairwise genetic distance analysis revealed intraspecific variation ranging from 0 to 0.23 % for S. nodosum from Vietnam, Myanmar, and Thailand; 0 to 0.08 % for S. shirakii from Taiwan; and 0 to 0.19 % for S. nodosum/S. shirakii (Table 2). The interspecific variation between S. nodosum/S. shirakii and S. nobile, and between S. nodosum/S. shirakii and S. timorense, ranged from 5.10 to 5.45 % and 5.18 to 5.45 %, respectively. Thus, the level of polymorphism (0–0.19 %) between S. nodosum and S. shirakii was less than that between all pairs of mainland populations and within the Myanmar population (0.04–0.23 %) and far below that of S. nobile or S. timorense.
Chromosomal banding patterns of all Taiwanese larvae were four fixed inversions (IS-1, IIL-a, IIIL-b, IIIL-2) removed from the Simulium subgeneric standard sequence, and matched the fixed banding sequence of all 247 Thai larvae studied by Tangkawanit et al. . Taiwanese and Thai populations had the nucleolar organizer in the end of IIS, undifferentiated sex chromosomes, and only two autosomal polymorphisms each, one of which (IIIL-1) was shared. The frequency of IIIL-1 was 0.99—nearly fixed—in Taiwan, and an average of 0.28 in Thai populations analyzed by Tangkawanit et al. . The chromosome map of Tangkawanit et al.  for IIIL, however, shows the homozygous sequence for the IIIL-1 inversion rather than the claimed standard sequence; in addition, the distal breakpoint is shown as extended by one band beyond the actual breakpoint. IIL-1 was an infrequent inversion (average frequency = 0.04) in Thailand  and absent in Taiwan. IL-1 was a new, but rare (0.01) inversion, with breakpoints before the last bands in each of sections 39 and 40, in one male larva in Taiwan. Ectopic pairing of centromeres formed a loose pseudochromocenter in populations in Thailand  but was present in only about 1 % of nuclei per larva in Taiwan.
Habitat characteristics for our collection sites of S. nodosum and S. shirakii, such as stream width, overlapped broadly. Altitude, however, was considerably greater for our mainland populations (733–1439 m) than for our Taiwanese populations (15–70 m) (Table 1).
All four data sources—molecular, chromosomal, morphological, and ecological—independently support the conspecificity of S. nodosum and S. shirakii across a longitudinal range of 2365 km, providing a powerful consensus that S. shirakii and mainland populations in Myanmar, Thailand, and Vietnam are a single species. Our multi-locus analysis demonstrates that genetic differences within mainland populations of S. nodosum are greater than the differences between S. nodosum collectively and S. shirakii.
Simulium nodosum from Thailand and S. shirakii from Taiwan are identical in all details of their fixed chromosomal inversions and sex chromosomes, and share one of three autosomal inversions, the other two being rare. Inversion IIIL-1, which is found in about one-third of all homologues of Thai larvae and is nearly fixed in Taiwanese larvae, possibly expresses clinal variation, with lowest polymorphism in the Taiwan population, a common characteristic of insular and peripheral populations .
Habitat characteristics that typically differ between closely related species, especially stream size , are broadly overlapping among populations from Taiwan and the mainland, in agreement with records from the literature [14, 15]. Known hosts, which can differ between closely related species , are consistent (bovids and humans) across the distribution [9, 11, 14]. Altitude, however, which can be associated with genetic isolation , differs between our populations in Taiwan (15–70 m) and on the mainland (733–1439 m). The altitudinal distinction, however, narrows or disappears when we draw from published information , suggesting that some of the discrepancy is attributable to sampling artifact. Elevation for S. nodosum in Thailand, for instance, ranged from 168 to 800 m . The available ecological data suggest that a broad range of habitat conditions is typical for S. nodosum and S. shirakii. Insular populations, in particular, are more likely to have a broader ecological niche, resulting in part from ecological release .
The molecular separation of S. nodosum and S. shirakii from S. nobile and S. timorense is in concordance with their morphological characters. The simple claws of the females and the number and form of the pupal gill filaments (three inflated horn-like tubes) separate S. nodosum/S. shirakii from other members of the S. nobile species group .
Given the conspecificity of S. shirakii with mainland populations, what taxonomic decision should be applied to populations across the entire range of S. nodosum and S. shirakii, from the type locality in India to Taiwan? Attempts by colleagues from India to collect S. nodosum from its type locality, approximately 2500 km from our nearest site of analysis (Thailand), did not produce S. nodosum. Lacking topotypical material for molecular and chromosomal analyses, we acknowledge two alternative possibilities: (1) recognition of two species—S. nodosum, represented by the type from India, and S. shirakii, represented by the material in our study, from Myanmar to Taiwan, or (2) recognition of a single species, S. nodosum, across the full range from India to Taiwan. We provisionally select the latter option on pragmatic grounds; the same name, when applied across the entire range, holds greater information content, emphasizing the morphological and ecological similarity of all populations. We, therefore, synonymize S. shirakii with S. nodosum.
This study was financially supported by research grants from the University of Malaya (RP003C-13SUS and RP003A-13SUS) and the Fundamental Research Grant Scheme (FP016-2012A).
- Post RJ, Cheke RA, Boakye DA, Wilson MD, Osei-Atweneboana MY, Tetteh-Kumah A, et al. Stability and change in the distribution of cytospecies of the Simulium damnosum complex (Diptera: Simuliidae) in southern Ghana from 1971 to 2011. Parasit Vectors. 2013;13(6):205.View ArticleGoogle Scholar
- Procunier WS. Cytological approaches to simuliid biosystematics in relation to the epidemiology and control of human onchocerciasis. Genome. 1989;32(4):559–69.View ArticlePubMedGoogle Scholar
- Post RJ, Mustapha M, Krueger A. Taxonomy and inventory of the cytospecies and cytotypes of the Simulium damnosum complex (Diptera: Simuliidae) in relation to onchocerciasis. Trop Med Int Health. 2007;12(11):1342–53.View ArticlePubMedGoogle Scholar
- Adler PH, Cheke RA, Post RJ. Evolution, epidemiology, and population genetics of black flies (Diptera: Simuliidae). Infect Genet Evol. 2010;10(7):846–65.View ArticlePubMedGoogle Scholar
- Adler PH, Giberson DJ, Purcell LA. Insular black flies (Diptera: Simuliidae) of North America: tests of colonization hypotheses. J Biogeogr. 2005;32(2):211–20.View ArticleGoogle Scholar
- Adler PH. Biodiversity of biting flies: implications for humanity. In: Foottit RG, Adler PH, editors. Insect biodiversity: science and society. Chichester: Wiley-Blackwell Publishing; 2009. p. 523–45.View ArticleGoogle Scholar
- Adler PH, Crosskey RW. World blackflies (Diptera: Simuliidae): a comprehensive revision of the taxonomic and geographical inventory. 2015, http://www.clemson.edu/cafls/biomia/pdfs/blackflyinventory.pdf. Accessed 26 Mar 2015.
- Takaoka H, Suzuki H. The black flies from Thailand (Diptera: Simuliidae). Jap J Sanit Zool. 1984;35(1):7–45.Google Scholar
- Datta M. An overview of the Simuliidae (Diptera) of West Bengal. India J Bengal Nat Hist Soc. 1992;11:41–62.Google Scholar
- Takaoka H, Choochote W, Aoki C, Fukuda M, Bain O. Black flies (Diptera: Simuliidae) attracted to humans and water buffalos and natural infections with filarial larvae, probably Onchocerca sp., in northern Thailand. Parasite. 2003;10(1):3–8.View ArticlePubMedGoogle Scholar
- Choochote W, Takaoka H, Fukuda M, Otsuka Y, Aoki C, Eshima N. Seasonal abundance and daily flying activity of black flies (Diptera: Simuliidae) attracted to human baits in Doi Inthanon National Park, northern Thailand. Med Entomol Zool. 2005;56(4):335–48.Google Scholar
- Henderson CAP. Homosequential species 2a and 2b within the Prosimulium onychodactylum complex (Diptera): temporal heterogeneity, linkage disequilibrium, and Wahlund effect. Can J Zool. 1986;64(4):859–66.View ArticleGoogle Scholar
- Ilmonen J, Adler PH, Malmqvist B, Cywinska A. The Simulium vernum group (Diptera: Simuliidae) in Europe: multiple character sets for assessing species status. Zool J Linn Soc. 2009;156(4):847–63.View ArticleGoogle Scholar
- Takaoka H. The black flies of Taiwan (Diptera: Simuliidae). Pac Insect. 1979;20(4):365–403.Google Scholar
- Tangkawanit U, Kuvangkadilok C, Trinachartvanit W, Baimai V. Cytotaxonomy, morphology, and ecology of the Simulium nobile species group (Diptera: Simuliidae) in Thailand. Cytogenet Genome Res. 2011;134(4):308–18.View ArticlePubMedGoogle Scholar
- Datta M. Note on the evolution of Simulium asishi Datta from Simulium nodosum Puri (Diptera: Simuliidae) in West Bengal. India Proc Zool Soc (Calcutta). 1992;45:123–7.Google Scholar
- Puri IM. Studies on Indian Simuliidae. Part VII. Descriptions of larva, pupa and female of Simulium nodosum sp. nov., with an appendix dealing with S. novolineatum nov. nom. (= S. lineatum Puri). Ind J Med Res. 1933;20(3):813–7.Google Scholar
- Takaoka H, Davies DM. The black flies (Diptera: Simuliidae) of Java, Indonesia. Bishop Museum Bulletin in Entomology 6. Honolulu: Bishop Museum Press; 1996.Google Scholar
- Takaoka H, Hadi UK, Sigit SH. The black flies (Diptera: Simuliidae) of Flores and Timor. Indonesia Med Entomol Zool. 2006;57(1):1–26.Google Scholar
- Xiong B, Kocher TD. Comparison of mitochondrial DNA sequences of seven morphospecies of black flies (Diptera: Simuliidae). Genome. 1991;34(2):306–11.View ArticlePubMedGoogle Scholar
- Phayuhasena S, Colgan DJ, Kuvangkadilok C, Pramual P, Baimai V. Phylogenetic relationships among the black fly species (Diptera: Simuliidae) of Thailand based on multiple gene sequences. Genetica. 2010;138(6):633–48.View ArticlePubMedGoogle Scholar
- Conflitti IM, Spironello M, Currie DC. Molecular phylogeny of the North American Simulium malyschevi-group species (Diptera: Simuliidae). Syst Entomol. 2012;37(3):571–7.View ArticleGoogle Scholar
- Low VL, Takaoka H, Adler PH, Ya’cob Z, Norma-Rashid Y, Chen CD, et al. A multi-locus approach resolves the phylogenetic relationships of the Simulium asakoae and Simulium ceylonicum species groups (Diptera: Simuliidae) in Malaysia: evidence for distinct evolutionary lineages. Med Vet Entomol. doi: https://doi.org/10.1111/mve.12120.
- Folmer O, Black M, Hoeh W, Lutz R, Vrijenhoek R. DNA primers for amplification of mitochondrial cytochrome c oxidase subunit I from diverse metazoan invertebrates. Mol Mar Biol Biotechnol. 1994;3(5):294–9.PubMedGoogle Scholar
- Kocher TD, Thomas WK, Meyer A, Edwards SV, Pääbo S, Villablanca FX, et al. Dynamics of mitochondrial DNA evolution in animals: amplification and sequencing with conserved primers. Proc Natl Acad Sci U S A. 1989;86(16):6196–200.View ArticlePubMed CentralPubMedGoogle 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(6):651–701.View ArticleGoogle Scholar
- Low VL, Adler PH, Takaoka H, Ya’cob Z, Lim PE, Tan TK, et al. Mitochondrial DNA markers reveal high genetic diversity but low genetic differentiation in the black fly Simulium tani Takaoka & Davies along an elevational gradient in Malaysia. PLoS One. 2014;9(6):e100512.View ArticlePubMed CentralPubMedGoogle Scholar
- Hall TA. BioEdit: a user-friendly biological sequence alignment editor and analysis program for Windows 95/98/NT. Nucleic Acids Symp Ser. 1999;41:95–8.Google Scholar
- Swofford DL. PAUP*. Phylogenetic Analysis Using Parsimony (*and other Methods). USA: Sinauer Associates; 2002.Google Scholar
- Huelsenbeck JP, Ronquist F. MRBAYES: Bayesian inference of phylogenetic trees. Bioinformatics. 2001;17(8):754–5.View ArticlePubMedGoogle Scholar
- Jobb G, von Haeseler A, Strimmer K. Treefinder: a powerful graphical analysis environment for molecular phylogenetics. BMC Evol Biol. 2004;4:18.View ArticlePubMed CentralPubMedGoogle Scholar
- Adler PH, Huang YT, Reeves WK, Kim SK, Otsuka Y, Takaoka H. Macrogenomic evidence for the origin of the black fly Simulium suzukii (Diptera: Simuliidae) on Okinawa Island, Japan. PLoS One. 2013;8(8):1–13.View ArticleGoogle Scholar
- Rothfels K, Feraday R, Kaneps A. A cytological description of sibling species of Simulium venustum and S. verecundum with standard maps for the subgenus Simulium Davies [sic] (Diptera). Can J Zool. 1978;56(5):1110–28.View ArticleGoogle Scholar
- Tangkawanit U, Kuvangkadilok C, Baimai V, Adler PH. Cytosystematics of the Simulium tuberosum group (Diptera: Simuliidae) in Thailand. Zool J Linn Soc. 2009;155(2):289–315.View ArticleGoogle Scholar
- Adler PH. Ecology of black fly sibling species. In: Kim KC, Merritt RW, editors. Black flies: ecology, population management and annotated world list. Pennsylvania: Pennsylvania State University Press; 1988. p. 63–76.Google Scholar
- Krueger A. Zoophily and animal filariae in Simulium damnosum s.l.: on the indicative property of tarsal claws in females. Acta Entomol Serbica. 2006;206(Suppl):39–50.Google Scholar
- Pramual P, Pangjanda S. Effects of habitat specialization on population genetic structure of black fly Simulium weji Takaoka (Diptera: Simuliidae). J Asia-Pacific Entomol. 2015;18(1):33–7.View ArticleGoogle Scholar
- Osborne PL. Tropical ecosystems and ecological concepts. 2nd ed. New York: Cambridge University Press; 2007.Google Scholar