Molecular systematics of the Philippine malaria vector Anopheles flavirostris

Allozyme and molecular sequence data from the malaria vector Anopheles flavirostris (Ludlow) (Diptera: Culicidae) were analysed from 34 sites throughout the Philippines, including the type locality, to test the hypothesis that this taxon is a single panmictic species. A finer-scaled allozyme study, of mainly Luzon samples, revealed no fixed genetic differences in sympatric sites and only low levels of variation. We obtained data from partial sequences for the internal transcribed spacer 2 (ITS2) (483 bp), the third domain (D3) (330 bp) of the 28S ribosomal DNA subunit and cytochrome c oxidase subunit I (COI) of mitochondrial DNA (261 bp). No sequence variation was observed for ITS2, only a one base pair difference was observed between Philippine and Indonesian D3 sequences and An. flavirostris sequences were unique, confirming their diagnostic value for this taxon. Sixteen COI haplotypes were identified, giving 25 parsimony informative sites. Neighbour-Joining, Maximum Parsimony, Maximum Likelihood and Bayesian phylogenetic analysis of COI sequences for An. flavirostris and outgroup taxa revealed strong branch support for the monophyly of An. flavirostris, thus confirming that Philippine populations of this taxon comprise a single separate species within the Minimus Subgroup of the Funestus Group. Variation in the behaviour of An. flavirostris is likely to be intraspecific rather than interspecific in origin.     

DNA barcoding of Dutch birds

The mitochondrial cytochrome c oxidase subunit I (COI) can serve as a fast and accurate marker for the identification of animal species, and has been applied in a number of studies on birds. We here sequenced the COI gene for 387 individuals of 147 species of birds from the Netherlands, with 83 species being represented by > 2 sequences. The Netherlands occupies a small geographic area and 95% of all samples were collected within a 50 km radius from one another. The intraspecific divergences averaged 0.29% among this assemblage, but most values were lower; the interspecific divergences averaged 9.54%. In all, 95% of species were represented by a unique barcode, with 6 species of gulls and skua (Larus and Stercorarius) having at least one shared barcode. This is best explained by these species representing recent radiations with ongoing hybridization. In contrast, one species, the Lesser Whitethroat Sylvia curruca showed deep divergences, averaging 5.76% and up to 8.68% between individuals. These possibly represent two distinct taxa, S. curruca and S. blythi, both clearly separated in a haplotype network analysis. Our study adds to a growing body of DNA barcodes that have become available for birds, and shows that a DNA barcoding approach enables to identify known Dutch bird species with a very high resolution. In addition some species were flagged up for further detailed taxonomic investigation, illustrating that even in ornithologically well-known areas such as the Netherlands, more is to be learned about the birds that are present.     

Cytogenetic,hybridization and molecular evidence of four cytological forms of Anopheles nigerrimus (Hyrcanus Group) in Thailand and Cambodia

Thirteen isoline colonies of Anopheles nigerrimus were established from individual wild‐caught females collected from cow‐baited traps at locations in Thailand and Cambodia. Three types of X (X₁, X₂, X₃) and 4 types of Y (Y₁, Y₂, Y₃, Y₄) chromosomes were recovered, according to differing amounts of extra heterochromatin. Four karyotypic forms were designed depending upon apparently distinct figures of X and Y chromosomes, i.e., Form A (X₁, X₂, X₃, Y₁), B (X₂, X₃, Y₂), C (X₁, Y₃), and D (X₃, Y₄). Forms C and D were new metaphase karyotypes discovered in this study. Form A appeared to be common in both Thailand and Cambodia. Forms B and D were found to be rather specific to southern and northeastern Thailand, respectively, whereas Form C was confined to Cambodia. Hybridization experiments among the eight isoline colonies, which were representative of four karyotypic forms of An. nigerrimus, demonstrated genetic compatibility in giving viable progenies and synaptic salivary gland polytene chromosomes through F₂‐generations. These results elucidated the conspecific relationship, comprising four cytological forms within this taxon. Supportive evidence was confirmed further by very low intraspecific sequence variations (average genetic distance = 0.002–0.007) of the nucleotide sequences in ribosomal DNA [second internal transcribed spacer (ITS2)] and mitochondrial DNA [cytochrome c oxidase subunit I (COI) and subunit II (COII)].     

A new species of Docosia Winnertz from Central Europe,with DNA barcoding based on four gene markers (Diptera,Mycetophilidae)

A new species of Docosia Winnertz, Docosia dentata sp. n., is described and illustrated, based on a single male specimen collected in Muránska planina National Park in Central Slovakia. DNA sequences (COI, COII, CytB, and ITS2) are included and compared for 13 species of Docosia. There was found only little congruence between the molecular results and previous scarce data about interspecific relationships based on morphology. The COI and CytB gene markers showed the highest interspecific gene distances while ITS2 showed the lowest ones. An updated key to the 23 Central European species of Docosia is also presented.     

ITS-2 secondary structures and phylogeny of Anopheles culicifacies species

Background: Second internal transcribed spacer (ITS2) has proven to contain useful biological information at higher taxonomic levels. Objectives: This study was carried out to unravel the biological information in the ITS2 region of An. culicifacies and the internal relationships between the five species of Anopheles culicifacies. Methodology: In achieving these objectives, twenty two ITS2 sequences (~370bp) of An. culicifacies species were retrieved from GenBank and secondary structures were generated. For the refinement of the primary structures, i.e. nucleotide sequence of ITS2 sequences, generated secondary structures were used. The improved ITS2 primary structures sequences were then aligned and used for the construction of phylogenetic trees. Results and discussions: ITS2 secondary structures of culicifacies closely resembled near universal eukaryotes secondary structure and had three helices, and the structures of helix II and distal region of helix III of ITS2 of An. culicifacies were strikingly similar to those regions of other organisms strengthening possible involvement of these regions in rRNA biogenesis. Phylogenetic analysis of improved ITS2 sequences revealed two main clades one representing sibling B, C and E and A and D in the other. Conclusions: Near sequence identity of ITS2 regions of the members in a particular clade indicate that this region is undergoing parallel evolution to perform clade specific RNA biogenesis. The divergence of certain isolates of An. culicifacies from main clades in phylogenetic analyses suggests the possible existence of camouflaged sub-species within the complex of culicifacies. Using the fixed nucleotide differences, we estimate that these two clades have diverged nearly 3.3 million years ago, while the sibling species in clade 2 are under less evolutionary pressure, which may have evolved much later than the members in clade 1.     

DNA barcoding for the identification of eight species members of the Thai Hyrcanus Group and investigation of their stenogamous behavior

Eight species members of the Thai Hyrcanus Group were identified based on the intact morphology and molecular analysis (COI barcoding, 658 bp) of F₁-progenies. Five iso-female lines of each species were pooled in order to establish stock colonies. A stenogamous colony of each species was investigated by making 200 and 300 newly emerged adult females and males co-habit in a 30 cm cubic cage for one week. After ovipositon, the spermathecae of females were examined for sperms. The results revealed that Anopheles argyropus, Anopheles crawfordi, Anopheles nitidus, Anopheles pursati, Anopheles sinensis, Anopheles nigerrimus, Anopheles paraliae and Anopheles peditaeniatus yielded insemination rates of 0%, 0%, 0%, 31%, 33%, 42%, 50% and 77%, respectively. Continuous selection to establish stenogamous colonies indicated that An. sinensis, An. pursati, An. nigerrimus, An. paraliae and An. peditaeniatus provided insemination rates of 33–34%, 27–31%, 42–58%, 43–57% and 61–86% in 1, 2, 5, 6 and 20 generations of passages, respectively.     

Molecular taxonomy of members of the Anopheles hyrcanus group from Thailand and Indonesia

During studies of malaria vectors in Indonesia and Thailand, several specimens identified by field staff as members of the Anopheles barbirostris group (Diptera: Culicidae) were found to belong to the Anopheles hyrcanus group, as shown by marked differences in the size of the nuclear rDNA second internal transcribed spacer (ITS2) between the barbirostris (~1500 bp) and hyrcanus (~600 bp) groups. Identification of the species concerned required a more detailed study of ITS2 sequences and subunit I of the mitochondrial DNA cytochrome oxidase gene (COI). A phylogenetic analysis, based on Bayesian methods, revealed that the hyrcanus group specimens comprised five distinct clades, two of which corresponded with known species, Anopheles peditaeniatus and Anopheles sinensis. The remaining specimens formed three additional clades, for which there are no similar sequences in GenBank and which cannot be linked to previously described species. The misidentification of hyrcanus group species has important implications for malaria vector control; more comprehensive studies employing gene sequences are required to clarify the number of species in the group, their distribution and vector status.     

Molecular Variation in the Paragonimus heterotremus Complex in Thailand and Myanmar

Paragonimiasis is an important food-borne parasitic zoonosis caused by infection with lung flukes of the genus Paragonimus. Of the 7 members of the genus known in Thailand until recently, only P. heterotremus has been confirmed as causing human disease. An 8th species, P. pseudoheterotremus, has recently been proposed from Thailand, and has been found in humans. Molecular data place this species as a sister species to P. heterotremus, and it is likely that P. pseudoheterotremus is not specifically distinct from P. heterotremus. In this study, we collected metacercariae of both nominal species (identification based on metacercarial morphology) from freshwater crabs from Phetchabun Province in northern Thailand, Saraburi Province in central Thailand, and Surat Thani Province in southern Thailand. In addition, we purchased freshwater crabs imported from Myanmar at Myawaddy Province, western Thailand, close to the Myanmar-Thailand border. The DNAs extracted from excysted metacercariae were PCR-amplified and sequenced for ITS2 and cox1 genes. The ITS2 sequences were nearly identical among all samples (99-100%). Phylogenies inferred from all available partial cox1 sequences contained several clusters. Sequences from Indian P. heterotremus formed a sister group to sequences from P. pseudoheterotremus-type metacercariae. Sequences of P. heterotremus from Thailand, Vietnam, and China formed a separate distinct clade. One metacercaria from Phitsanulok Province was distinct from all others. There is clearly considerable genetic variation in the P. heterotremus complex in Thailand and the form referred to as P. pseudoheterotremus is widely distributed in Thailand and the Thai-Myanmar border region.     

Late Pleistocene environmental changes lead to unstable demography and population divergence of Anopheles albimanus in the northern Neotropics

We investigated the historical demography of Anopheles albimanus using mosquitoes from five countries and three different DNA regions, the mitochondrial cytochrome oxidase subunit I gene (COI), the single copy nuclear white gene and the ribosomal internal transcribed spacer two (ITS2). All the molecular markers supported the taxonomic status of a single species of An. albimanus. Furthermore, agreement between the COI and the white genes suggested a scenario of Pleistocene geographic fragmentation (i.e., population contraction) and subsequent range expansion across southern Central America.     

Molecular Variation and Distribution of Anopheles fluviatilis (Diptera: Culicidae) Complex in Iran

Anopheles fluviatilis James (Diptera: Culicidae) is one of the known malaria vectors in south and southeastern Iran. Earlier ITS2 sequences analysis of specimens from Iran demonstrated only a single genotype that was identical to species Y in India, which is also the same as species T. We identified 2 haplotypes in the An. fluviatilis populations of Iran based on differences in nucleotide sequences of D3 domain of the 28S locus of ribosomal DNA (rDNA). Comparison of sequence data from 44 Iranian specimens with those publicly available in the Genbank database showed that all of the 28S-D3 sequences from Kazeroun and Khesht regions in Fars Province were identical to the database entry representing species U in India. In other regions, all the individuals showed heterozygosity at the single nucleotide position, which identifies species U and T. It is argued that the 2 species may co-occur in some regions and hybridize; however, the heterozygosity in the 28S-D3 locus was not reflected in ITS2 sequences and this locus for all individuals was identical to species T. This study shows that in a newly diverged species, like members of An. fluviatilis complex, a single molecular marker may not be sufficiently discriminatory to identify all the taxa over a vast geographical area. In addition, other molecular markers may provide more reliable information for species discrimination.     

Mitochondrial DNA variation in the malaria vector Anopheles minimus across China,Thailand and Vietnam: evolutionary hypothesis,population structure and population history

The effects of Pleistocene environmental fluctuations on the distribution and diversity of organisms in Southeast Asia are much less well known than in Europe and North America. In these regions, the combination of palaeoenvironmental reconstruction and inferences about population history from genetic data has been very powerful. In Southeast Asia, mosquitoes are good candidates for the genetic approach, with the added benefit that understanding the relative contributions of historical and current processes to population structure can inform management of vector species. Genetic variation among populations of Anopheles minimus was examined using 144 mtDNA COII sequences from 23 sites in China, Thailand and Vietnam. Haplotype diversity was high, with two distinct lineages that have a sequence divergence of over 2% and exhibit different geographical distributions. We compare alternative hypotheses concerning the origin of this pattern. The observed data deviate from the expectations based on a single-panmictic population with or without growth, or a stable but spatially structured population. However, they can be readily accommodated by a model of past fragmentation into eastern and western refugia, followed by growth and range expansion. This is consistent with the palaeoenvironmental reconstructions currently available for the region.     

Description and DNA barcoding of Crematogaster fraxatrix Forel, 1911 and two new closely related species from Cambodia and Indonesia (Hymenoptera,Formicidae)

Crematogaster fraxatrix Forel, 1911 and two new species, C. chhangi sp. n. and C. simboloni sp. n., are described from Cambodia and Indonesia, respectively. DNA sequences were generated for C. fraxarix and the two newly described species using 3 amplications of two regions of the mitochondrial gene COI with a total of 1129 bp. The mean interspecific divergences are 9.4% and 23.5% for C. fraxatrix vs. C. chhangi, C. simboloni, respectively. DNA sequences reveal that C. simboloni is found to be genetically distinct from the other two species, but C. chhangi is not distinct from C. fraxatrix.     

The mosquito Anopheles (Cellia) oreios sp. n., formerly species 6 of the Australasian Anopheles farauti complex,and a critical review of its biology and relation to disease

Species 6 of the Australasian Anopheles farauti sibling species complex (Diptera: Culicidae) is described and formally named Anopheles oreios Bangs & Harbach, sp. n. Adult, pupal and fourth‐instar larval specimens collected in the Baliem Valley, Papua Province, Indonesia, are characterized and compared with those of Anopheles farauti, Anopheles hinesorum, Anopheles irenicus and Anopheles torresiensis (formerly informally denoted as species 1, 2, 7 and 3, respectively). The variable wings of adult females, the male genitalia, the pupa and the fourth‐instar larva of An. oreios are illustrated and DNA sequence data are included for regions coding for sections of the mitochondrial COI and COII genes. The biology of An. oreios and its relation to malaria transmission are discussed in detail and contrasted with the biology and disease relations of some members of the An. farauti and Anopheles punctulatus sibling species complexes.     

Assessing Nuclear and Mitochondrial DNA Sequence Variation Within Steinernema (Rhabditida: Steinernematidae)

DNA sequence analysis was used to characterize the nuclear ribosomal DNA ITS1 region and a portion of the COII and 16S rDNA genes of the mitochondrial genome from Steinernema entomopathogenic nematodes. Nuclear ITS1 nucleotide divergence among seven Steinernema spp. ranged from 6 to 22%, and mtDNA divergence among five species ranged from 12 to 20%. No intraspecific variation was observed among three S. feltiae strains. Phylogenetic analysis of both nuclear and mitochondrial DNA sequences confirms the existing morphological relationships of several Steinernema species. Both the rDNA ITS1 and mtDNA sequences were useful for resolving relationships among Steinernema taxa.     

DNA barcoding and the differentiation between North American and West European Phormia regina (Diptera,?Calliphoridae,Chrysomyinae)

Phormia regina (the black fly) is a common Holarctic blow fly species which serves as a primary indicator taxon to estimate minimal post mortem intervals. It is also a major research model in physiological and neurological studies on insect feeding. Previous studies have shown a sequence divergence of up to 4.3% in the mitochondrial COI gene between W European and N American P. regina populations. Here, we DNA barcoded P. regina specimens from six N American and 17 W European populations and confirmed a mean sequence divergence of ca. 4% between the populations of the two continents, while sequence divergence within each continent was a ten-fold lower. Comparable mean mtDNA sequence divergences were observed for COII (3.7%) and cyt b (5.3%), but mean divergence was lower for 16S (0.4–0.6%). Intercontinental divergence at nuclear DNA was very low (≤ 0.1% for both 28S and ITS2), and we did not detect any morphological differentiation between N American and W European specimens. Therefore, we consider the strong differentiation at COI, COII and cyt b as intraspecific mtDNA sequence divergence that should be taken into account when using P. regina in forensic casework or experimental research.     

Distinct population structure for co-occurring Anopheles goeldii and Anopheles triannulatus in Amazonian Brazil

To evaluate whether environmental heterogeneity contributes to the genetic heterogeneity in Anopheles triannulatus, larval habitat characteristics across the Brazilian states of Roraima and Pará and genetic sequences were examined. A comparison with Anopheles goeldii was utilised to determine whether high genetic diversity was unique to An. triannulatus. Student t test and analysis of variance found no differences in habitat characteristics between the species. Analysis of population structure of An. triannulatus and An. goeldii revealed distinct demographic histories in a largely overlapping geographic range. Cytochrome oxidase I sequence parsimony networks found geographic clustering for both species; however nuclear marker networks depicted An. triannulatus with a more complex history of fragmentation, secondary contact and recent divergence. Evidence of Pleistocene expansions suggests both species are more likely to be genetically structured by geographic and ecological barriers than demography. We hypothesise that niche partitioning is a driving force for diversity, particularly in An. triannulatus.     

Evidence for a conspecific relationship between two morphologically and cytologically different forms of Korean Anopheles pullus mosquito

Despite the medical importance of anopheline mosquitoes as vectors of Korean vivax malaria, differentiation between Korean anopheline mosquitoes by traditional morphological taxonomic criteria is difficult. An. yatsushiroensis is the second most common Anopheles mosquito species in Korea and a possible vector of Korean vivax malaria together with An. sinensis, the predominant anopheline species. Recently, An. yatsushiroensis has been declared a synonym of An. pullus, based on comparisons of egg morphology and adult progeny, although they differ in ecology and morphology. To verify the species status of these two ambiguous forms, we established isofemale lines of Korean An. pullus and An. yatsushiroensis (An. pullus form yatsushiroensis) mosquitoes and investigated their genetic relationship by metaphase karyotype analysis, comparing the DNA sequences of rDNA internal transcribed spacer 2 (ITS2) and mitochondrial cytochrome c oxidase subunits I (COI) and II (COII), and by hybridization experiments. Two isofemale lines had differently shaped X and Y chromosomes. However reciprocal crosses between them yielded viable progeny with completely synaptic salivary gland polytene chromosomes. DNA analyses also strongly supported their conspecificity. The two strains also showed great sequence similarity in the ITS2, COI and COII regions (variation rate = 0.0 to 0.8%). Based on these findings, we conclude that the two forms, though differing distinctly in morphological, cytological and ecological traits, remain interfertile.     

Evidence to support two conspecific cytological races on Anopheles aconitus in Thailand.

Iso-female lines (isolines) of Anopheles aconitus collected from Mae Hong Son, Phet Buri, and Chiang Mai Provinces were successfully identified to karyotypic forms. The results of identification revealed that An. aconitus Form B (X1, X2, Y2) was obtained from four and 48 isolines in Phet Buri and Chiang Mai Provinces, respectively, and Form C (X1, X2, Y3) was recovered from three and 41 isolines in Mae Hong Son and Chiang Mai Provinces, respectively. When comparing band to band on the same arm of ovarian nurse cell polytene chromosomes of An. aconitus Form B (Phet Buri: four isolines) and C (Mae Hong Son: three isolines, Chiang Mai: 20 isolines) to the standard chromosome mapping of An. aconitus Form B (Chiang Mai: 20 isolines), no major chromosomal rearrangements that related to the karyotype variations were demonstrated. The investigations on allelic frequencies of 4th stage larvae and adult females of three (Form C: Mae Hong Son), four (Form B: Phet Buri), 41 (Form C: Chiang Mai) and 48 (Form B: Chiang Mai) isolines suggested that An. aconitus Form B and C of all strains have similar allelic frequencies. This was observed at 10 isoenzymes 16 loci in 4th stage larvae, and 11 isoenzymes 13 loci in adult females. Hybridization tests among the four laboratory-raised isolines of An. aconitus Form B (Chiang Mai and Phet Buri) and C (Chiang Mai and Mae Hong Son) were employed by induced copulation. The results of crosses indicated that they were genetically compatible, providing viable progeny and completely synaptic salivary gland polytene chromosomes. The complete sequences ofrDNA internal-transcribed spacer two (ITS2) and partial sequences of mitochondrial cytochrome c oxidase subunit I and II (COI and COII) from genomic DNA of 12 isolines of An. aconitus Form B and C were identified. Total sequence lengths (ITS2+COI+COII) of An. aconitus isolines varied from 1550bp to 1556bp. Conspecific relationships between the two An. aconitus forms were well supported by low values of intraspecific distances (ranged from 0.1% to 1.0%) and genetic differentiation (d(xy): 0.01322) between the two forms. Based on evidence of no pre- and post-mating isolations, and nearly identical of DNA sequences of ITS2, COI and COII regions between An. aconitus Form B and C, we conclude that they are conspecific cytological races in the Thai population.     

Multigene Phylogeography of Bactrocera caudata (Insecta: Tephritidae): Distinct Genetic Lineages in Northern and Southern Hemispheres

Bactrocera caudata is a pest of pumpkin flower. Specimens of B. caudata from the northern hemisphere (mainland Asia) and southern hemisphere (Indonesia) were analysed using the partial DNA sequences of the nuclear 28S rRNA and internal transcribed spacer region 2 (ITS-2) genes, and the mitochondrial cytochrome c oxidase subunit I (COI), cytochrome c oxidase subunit II (COII) and 16S rRNA genes. The COI, COII, 16S rDNA and concatenated COI+COII+16S and COI+COII+16S+28S+ITS-2 nucleotide sequences revealed that B. caudata from the northern hemisphere (Peninsular Malaysia, East Malaysia, Thailand) was distinctly different from the southern hemisphere (Indonesia: Java, Bali and Lombok), without common haplotype between them. Phylogenetic analysis revealed two distinct clades (northern and southern hemispheres), indicating distinct genetic lineage. The uncorrected ‘p’ distance for the concatenated COI+COII+16S nucleotide sequences between the taxa from the northern and southern hemispheres (‘p’ = 4.46-4.94%) was several folds higher than the ‘p’ distance for the taxa in the northern hemisphere (‘p’ = 0.00-0.77%) and the southern hemisphere (‘p’ = 0.00%). This distinct difference was also reflected by concatenated COI+COII+16S+28S+ITS-2 nucleotide sequences with an uncorrected ''p'' distance of 2.34-2.69% between the taxa of northern and southern hemispheres. In accordance with the type locality the Indonesian taxa belong to the nominal species. Thus the taxa from the northern hemisphere, if they were to constitute a cryptic species of the B. caudata species complex based on molecular data, need to be formally described as a new species. The Thailand and Malaysian B. caudata populations in the northern hemisphere showed distinct genetic structure and phylogeographic pattern.     

Genetic diversity of genus Vespa including an invaded species of V. velutina (Hymenoptera: Vespidae) in Korea inferred from DNA barcoding data

With about 5000 known species, the Vespidae is a large family belongs to order Hymenoptera. The genus Vespa with 22 species is one of the four genera of the subfamily Vespinae. In Korea, 10 species and subspecies are recognized. Because of their social behavior, their treat to human health and their impact in apiculture, the reliable and sometimes automated identification of these insects to species level are important. To test the efficacy of DNA barcoding method for identification of species of the genus Vespa in Korea, 30 samples of eight Korean species of genus Vespa were collected and mitochondrial DNAs of 658 bp fragment cytochrome oxidase subunit 1 (CO1) region were sequenced. A Bayesian Inference based on COI gene of the Korean Vespa species was constructed. The phylogenetic tree shoed that identification of all specimens is possible based on COI gene and we found strong relation between the sequences of the collected species from different localities in South Korea which clustered together with 100% support with sequences of the same species in GenBank. The results demonstrate that DNA barcoding is a useful technique for rapid and accurate species recognition in Korean Vespa species. The DNA barcode part of COI for V. binghami is provided for the first time that can help for identification of this species through DNA barcoding. Also, the genetic diversity among Korean Vespa velutina was zero suggests that the invasion might have occurred in a single event with small number of founders.