Intraspecific concerted evolution of the rDNA ITS1 in Anopheles farauti sensu stricto (Diptera: Culicidae) reveals recent patterns of population structure

Morphological and phylogenetic analysis of Raillietina spp. indigenous chickens (Gallus gallus domesticus) in Bangladesh

Raillietina spp. (Cestoda: Davaineidae), the most common cestodes in indigenous chickens, cause a substantial production loss in poultry industry in Bangladesh. Here, we estimated the prevalence, confirmed the species and determined the genetic pattern of species of Raillietina using molecular tools. We collected and examined 375 chickens randomly from household of different villages of Mymensingh sadar and Gouripur upazila, Mymensingh district and adult parasites were isolated and identified. Genomic DNA was extracted from collected parasites, amplified ITS-2 and ND-1 genes, sequenced and analyzed. Out of 375 samples, 270 (72.0%) were found positive with Raillietina species and mean worm burden was 10.46 ± 0.56. Microscopically, three species of Raillietina, such as R. cesticillus (37.9%), R. echinobothrida (41.1%) and R. tetragona (52.8%) were detected on the basis of their morphological features. The total length, length and width of scolex, sucker and rostellum were also measured. Among different factors, age, farming nature and flock size of chickens were significantly (p < 0.05) influenced Raillietina infections. For further validation, the sequences of ITS-2 gene generated in this study were matched with reference sequences of R. cesticillus, R. echinobothrida and R. tetragona and found 99.63% - 100% similarity. The phylogenetic analyses of ITS-2 and ND-1 sequences were clustered together with the reference sequences of R. cesticillus, R. echinobothrida and R. tetragona confirming microscopic identification. This is the first confirmation of species of Raillietina along with the prevalence of the species, which will be helpful for the formulation of a control strategy and provide basic information for further molecular study.

Phylogenetic analysis of Eimeria tenella isolated from the litter of different chicken farms in Mymensingh, Bangladesh

Background

Eimeria tenella is the most pathogenic intracellular protozoan parasite of seven Eimeria species causing chicken coccidiosis around the world. This species is particularly responsible for caecal coccidiosis leading to serious morbidity–mortality and financial loss in poultry production.

Methods

The present study explored the genetic diversity of E. tenella. Litter slurry was collected from 18 broiler farms located in Mymensingh district, Bangladesh. Litter samples were processed for oocyst isolation–identification using parasitological techniques followed by genomic DNA extraction from sporulated oocysts. For molecular analysis, the internaltranscribedspacer1 gene of E. tenella was amplified using species‐specific primers and sequenced. After editing and alignment, 263 bp sequences were used for analysis.

Results

Genetic analysis showed seven distinct genotypes and detected six single nucleotide polymorphisms among the 18 E. tenella isolates. The nucleotide and genotype diversity were 0.00507 and 0.8235, respectively. A phylogenetic tree was constructed with 66 sequences (seven studied genotypes and 59 reference sequences from GenBank database). The neighbour‐joining tree represented that the studied E. tenella isolates were grouped with reference E. tenella isolates with strong nodal support (100%) and the nucleotide sequences of E. tenella, E. necatrix, E. acervulina, E. brunetti, E. maxima, E. mitis and E. praecox formed separate clusters without any geographical boundaries.

Conclusions

This is the first study on the genetic analysis of E. tenella from Mymensingh district, Bangladesh. These findings will provide baseline data on the species conformation and genetic variations of E. tenella. Further extensive investigation will be needed to reveal the population genetic structure of this parasite and thus will facilitate the planning of effective control strategies.

Comparative Phylogeography and Integrative Taxonomy of Ochlerotatus caspius (Dipera: Culicidae) and Ochlerotatus dorsalis

Given that accurately identifying pathogen vectors is vital for designing efficient mosquito control programs based on the proper surveillance of the epidemiologically important species, it has been suggested the complementary use of independently evolving genes and morphometric traits as a reliable approach for the characterization and delimitation of related species. Hence, we examined the spatial distribution of COI mtDNA and ITS2 rDNA variation from the historical perspective of Ochlerotatus caspius (Pallas, 1771) and O. dorsalis (Meigen, 1830), while simultaneously testing the utility of the two markers in integrative species delimitation when combined with phenotypic character analyses of larvae and adults. Despite the striking difference in haplotype diversity (high in COI mtDNA, low in ITS2 rDNA), no evident phylogeographic structure was apparent in the Palearctic O. caspius. The Holarctic O. dorsalis species was subdivided into two highly distinctive COI mtDNA phylogroups which corresponded to the Nearctic and Palearctic regions. Strong support for the independence of the two allopatric evolutionary lineages suggested that geographical barrier and climatic changes during Pleistocene caused vicariance of the ancestral range. COI mtDNA reliably distinguished O. caspius and O. dorsalis, while ITS2 rDNA yet again lacked the proper resolution for solving this problem. An integrative approach based on the larval and adult morphological traits have varying taxonomic applications due to their differential diagnostic values. Thus, by the implementation of an integrative taxonomic approach, we successfully detected species borders between the two epidemiologically relevant species and uncovered the presence of cryptic diversity within O. dorsalis.

The epidemiology of Plasmodium falciparum and Plasmodium vivax in East Sepik Province, Papua New Guinea, pre- and post-implementation of national malaria control efforts

Background

In the past decade, national malaria control efforts in Papua New Guinea (PNG) have received renewed support, facilitating nationwide distribution of free long-lasting insecticidal nets (LLINs), as well as improvements in access to parasite-confirmed diagnosis and effective artemisinin-combination therapy in 2011–2012.

Methods

To study the effects of these intensified control efforts on the epidemiology and transmission of Plasmodium falciparum and Plasmodium vivax infections and investigate risk factors at the individual and household level, two cross-sectional surveys were conducted in the East Sepik Province of PNG; one in 2005, before the scale-up of national campaigns and one in late 2012-early 2013, after 2 rounds of LLIN distribution (2008 and 2011–2012). Differences between studies were investigated using Chi square (χ²), Fischer’s exact tests and Student’s t-test. Multivariable logistic regression models were built to investigate factors associated with infection at the individual and household level.

Results

The prevalence of P. falciparum and P. vivax in surveyed communities decreased from 55% (2005) to 9% (2013) and 36% to 6%, respectively. The mean multiplicity of infection (MOI) decreased from 1.8 to 1.6 for P. falciparum (p = 0.08) and from 2.2 to 1.4 for P. vivax (p < 0.001). Alongside these reductions, a shift towards a more uniform distribution of infections and illness across age groups was observed but there was greater heterogeneity across the study area and within the study villages. Microscopy positive infections and clinical cases in the household were associated with high rate infection households (> 50% of household members with Plasmodium infection).

Conclusion

After the scale-up of malaria control interventions in PNG between 2008 and 2012, there was a substantial reduction in P. falciparum and P. vivax infection rates in the studies villages in East Sepik Province. Understanding the extent of local heterogeneity in malaria transmission and the driving factors is critical to identify and implement targeted control strategies to ensure the ongoing success of malaria control in PNG and inform the development of tools required to achieve elimination. In household-based interventions, diagnostics with a sensitivity similar to (expert) microscopy could be used to identify and target high rate households.

Recent and old duplications in crustaceans "Internal Transcribed Spacer 1″: structural and phylogenetic implications

Internal Transcribed Spacer structures are important in preserving accessibility to specific enzymes for the maturation of rRNAs. ITS1 sequences reported in the literature in Crustaceans range between 182 and 820 bp and are characterized by the absence of repeats or the presence of only a limited number of microsatellites. Here, we sequenced ITS1 for a range of shrimp families (infraorder Caridea) and show that most taxa have much larger ITS1 sequences. We find a high number of microsatellites in Alpheus hebes and Crangon crangon and we report repeat units in Pandalidae, Palaemonidae and mainly in Alpheidae species. Up to four repeats were found in A. vanderbilti (1915 bp), A. rostratus (1635 bp) and A. lottini (1625 bp). In general, four helices were found in ITS1. Repeat units led to extra hairpins and loops. No conserved positions occurred except in helix 4. Three clades were defined in A. lottini for the first time. We estimated the ITS1 divergence rate for the three clades of A. lottini collected in French Polynesia using existing calibrations of substitution rates. Rates of sequence evolution are largely influenced by repeat units, which likely evolve separately. By comparison with COI marker, we estimated the divergence rate of the whole ITS1 sequence to range from 0.5 to 1.4% Pmy and between 0.12 and 0.5% for the 3' end of ITS1 located outside the repeat units. Given the degree of identity between repeats, we suggest that a duplication event recently occurred in A. floridanus (98% identity) whereas an ancient duplication happened in A. sulcatus (50% identity) early at the origination of the group Alpheidae, approximately 50 mya ago. In conclusion, our results highlight an over representation of shorter ITS1 sequences in public repositories, and underlines the importance to further understand patterns of molecular evolution of this functionally important gene.

Molecular Identification and Phylogenetic Analysis of Nuclear rDNA Sequences of Clonorchis sinensis Isolates From Human Fecal Samples in Heilongjiang Province, China

Studying the genetic diversity of parasite is important for understanding their biogeography and molecular epidemiology, as well as for establishing disease prevention and control strategies. Clonorchis sinensis is an important foodborne parasite worldwide. However, despite its epidemiological significance, the genetic diversity of C. sinensis has not been well studied from human in northeastern China. In this study, a total of 342 fecal specimens were collected from residents living in five villages in Heilongjiang Province and analyzed for the presence of C. sinensis by PCR amplification and sequencing of the internal transcribed spacer 1 (ITS1) and ITS2 regions of nuclear ribosomal DNA. 21.64% (74/342) of fecal samples were found to be positive for C. sinensis by PCR. The sequences of the ITS1 region in 34 of the 74 samples (45.95%) matched that of MK179278, Genetic polymorphisms were observed at six nucleotide sites. The ITS2 gene sequence of 37 of the 74 samples (50%) matched that of MK179281. In conclusion, a low degree of genetic diversity between C. sinensis isolates from China and different geographical regions was found at ITS loci. Despite this conservation, sequencing of the rDNA region has provided important data that will be useful for future studies addressing the molecular evolution, biology, medical implications and ecology of C. sinensis.

DNA barcoding mosquitoes: advice for potential prospectors

Mosquitoes’ importance as vectors of pathogens that drive disease underscores the importance of precise and comparable methods of taxa identification among their species. While several molecular targets have been used to study mosquitoes since the initiation of PCR in the 1980s, its application to mosquito identification took off in the early 1990s. This review follows the research's recent journey into the use of mitochondrial DNA (mtDNA) cytochrome oxidase 1 (COI or COX1) as a DNA barcode target for mosquito species identification – a target whose utility for discriminating mosquitoes is now escalating. The pros and cons of using a mitochondrial genome target are discussed with a broad sweep of the mosquito literature suggesting that nuclear introgressions of mtDNA sequences appear to be uncommon and that the COI works well for distantly related taxa and shows encouraging utility in discriminating more closely related species such as cryptic/sibling species groups. However, the utility of COI in discriminating some closely related groups can be problematic and investigators are advised to proceed with caution as problems with incomplete lineage sorting and introgression events can result in indistinguishable COI sequences appearing in reproductively independent populations. In these – if not all – cases, it is advisable to run a nuclear marker alongside the mtDNA and thus the utility of the ribosomal DNA – and in particular the internal transcribed spacer 2 – is also briefly discussed as a useful counterpoint to the COI.

First Records of Culicoides sonorensis (Diptera: Ceratopogonidae), a Known Vector of Bluetongue Virus, in Southern Ontario

Ceratopogonidae (Diptera) were collected on sheep farms in southern Ontario to establish whether Culicoides spp. pose a threat to the livestock industry. Specimens were collected in modified CO2-baited Centers for Disease Control and Prevention light traps, returned to the laboratory, freeze-killed, and identified to species under a microscope. In addition to Culicoides variipennis (Coquillet), we found that Culicoides sonorensis Wirth & Jones occurred on a number of farms over a 2-yr period. These records represent a significant departure from C. sonorensis' previously known geographical distribution. We present spatial and temporal distribution data for both species, with an emphasis on C. sonorensis. DNA sequence information is presented so that researchers lacking the necessary taxonomic skills can determine whether C. sonorensis is present in their collections. To differentiate C. sonorensis from C. variipennis, taxonomically reliable and informative traits were found in EF1α and, to a lesser extent, in ITS1, whereas the universal barcode region of cytochrome oxidase subunit 1 (CO1) was unsuitable.

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.

Anopheles punctulatus group: evolution, distribution, and control

The major malaria vectors of the Southwest Pacific belong to a group of closely related mosquitoes known as the Anopheles punctulatus group. The group comprises 13 co-occurring species that either are isomorphic or carry overlapping morphological features, and today several species remain informally named. The advent of species-diagnostic molecular tools in the 1990s permitted a new raft of studies into the newly differentiated mosquitoes of this group, and these have revealed five species as the region's primary malaria vectors: An. farauti, An. hinesorum, An. farauti 4, An. koliensis, and An. punctulatus. Species' distributions are now well established across Papua New Guinea, northern Australia, and the Solomon Archipelago, but little has been documented thus far in eastern Indonesia. As each species reveals significant differences in distribution and biology, the relative paucity of knowledge of their biology or ecology in relation to malaria transmission is brought into clearer focus. Only three of the species have undergone some form of spatial or population genetics analyses, and this has revealed striking differences in their genetic signatures throughout the region. This review compiles and dissects the key findings for this important mosquito group and points to where future research should focus to maximize the output of field studies in developing relevant knowledge on these malaria vectors.

Capturing the population structure of microparasites: using ITS-sequence data and a pooled DNA approach

The internal transcribed spacer (ITS) region of nuclear ribosomal DNA is a common marker not only for the molecular identification of different taxa and strains, but also for the analysis of population structure of wild microparasite communities. Importantly, the multicopy nature of this region allows the amplification of low-quantity samples of the target DNA, a common problem in studies on unicellular, unculturable microparasites. We analysed ITS sequences from the protozoan parasite Caullerya mesnili (class Ichthyosporea) infecting waterflea (Daphnia) hosts, across several host population samples. We showed that analysing representative ITS-types [as identified by statistical parsimony network analysis (SPN)] is a suitable method to address relevant polymorphism. The spatial patterns were consistent regardless of whether parasite DNA was extracted from individual hosts or pooled host samples. Remarkably, the efficiency in detecting different sequence types was even higher after sample pooling. As shown by simulations, an easily manageable number of sequences from pooled DNA samples are sufficient to resolve the spatial population structure in this system. In summary, the ITS region analysed from pooled DNA samples can provide valuable insights into the spatial and temporal dynamics of microparasites. Moreover, the application of SPN analysis is a good alternative to the well-established neighbour-joining method (NJ) for the identification of representative ITS-types. SPN can even outperform NJ by joining most of the singleton sequences to representative sequence clusters.

Glossina palpalis palpalis populations from Equatorial Guinea belong to distinct allopatric clades

Background

Luba is one of the four historical foci of Human African Trypanosomiasis (HAT) on Bioko Island, in Equatorial Guinea. Although no human cases have been detected since 1995, T. b. gambiense was recently observed in the vector Glossina palpalis palpalis. The existence of cryptic species within this vector taxon has been previously suggested, although no data are available regarding the evolutionary history of tsetse flies populations in Bioko.

Methods

A phylogenetic analysis of 60 G. p. palpalis from Luba was performed sequencing three mitochondrial (COI, ND2 and 16S) and one nuclear (rDNA-ITS1) DNA markers. Phylogeny reconstruction was performed by Distance Based, Maximum Likelihood and Bayesian Inference methods.

Results

The COI and ND2 mitochondrial genes were concatenated and revealed 10 closely related haplotypes with a dominant one found in 61.1% of the flies. The sequence homology of the other 9 haplotypes compared to the former ranged from 99.6 to 99.9%. Phylogenetic analysis clearly clustered all island samples with flies coming from the Western African Clade (WAC), and separated from the flies belonging to the Central Africa Clade (CAC), including samples from Mbini and Kogo, two foci of mainland Equatorial Guinea. Consistent with mitochondrial data, analysis of the microsatellite motif present in the ITS1 sequence exhibited two closely related genotypes, clearly divergent from the genotypes previously identified in Mbini and Kogo.

Conclusions

We report herein that tsetse flies populations circulating in Equatorial Guinea are composed of two allopatric subspecies, one insular and the other continental. The presence of these two G. p. palpalis cryptic taxa in Equatorial Guinea should be taken into account to accurately manage vector control strategy, in a country where trypanosomiasis transmission is controlled but not definitively eliminated yet.

Novel genetic diversity within Anopheles punctimacula s.l.: phylogenetic discrepancy between the Barcode cytochrome c oxidase I (COI) gene and the rDNA second internal transcribed spacer (ITS2)

Anopheles punctimacula s.l. is a regional malaria vector in parts of Central America, but its role in transmission is controversial due to its unresolved taxonomic status. Two cryptic species, An. malefactor and An. calderoni, have been previously confused with this taxon, and evidence for further genetic differentiation has been proposed. In the present study we collected and morphologically identified adult female mosquitoes of An. punctimacula s.l. from 10 localities across Panama and one in Costa Rica. DNA sequences from three molecular regions, the three prime end of the mitochondrial cytochrome c oxidase I gene (3' COI), the Barcode region in the five prime end of the COI (5' COI), and the rDNA second internal transcribed spacer (ITS2) were used to test the hypothesis of new molecular lineages within An. punctimacula s.l. Phylogenetic analyses using the 3' COI depicted six highly supported molecular lineages (A-F), none of which was An. malefactor. In contrast, phylogenetic inference with the 5' COI demonstrated paraphyly. Tree topologies based on the combined COI regions and ITS2 sequence data supported the same six lineages as the 3' COI alone. As a whole this evidence suggests that An. punctimacula s.l. comprises two geographically isolated lineages, but it is not clear whether these are true species. The phylogenetic structure of the An. punctimacula cluster as well as that of other unknown lineages (C type I vs C type II; D vs E) appears to be driven by geographic partition, because members of these assemblages did not overlap spatially. We report An. malefactor for the first time in Costa Rica, but our data do not support the presence of An. calderoni in Panama.

Phylogeography of the neotropical Anopheles triannulatus complex (Diptera: Culicidae) supports deep structure and complex patterns

BACKGROUND

The molecular phylogenetic relationships and population structure of the species of the Anopheles triannulatus complex: Anopheles triannulatus s.s., Anopheles halophylus and the putative species Anopheles triannulatus C were investigated.

METHODS

The mitochondrial COI gene, the nuclear white gene and rDNA ITS2 of samples that include the known geographic distribution of these taxa were analyzed. Phylogenetic analyses were performed using Bayesian inference, Maximum parsimony and Maximum likelihood approaches.

RESULTS

Each data set analyzed septely yielded a different topology but none provided evidence for the seption of An. halophylus and An. triannulatus C, consistent with the hypothesis that the two are undergoing incipient speciation. The phylogenetic analyses of the white gene found three main clades, whereas the statistical parsimony network detected only a single metapopulation of Anopheles triannulatus s.l. Seven COI lineages were detected by phylogenetic and network analysis. In contrast, the network, but not the phylogenetic analyses, strongly supported three ITS2 groups. Combined data analyses provided the best resolution of the trees, with two major clades, Amazonian (clade I) and trans-Andean + Amazon Delta (clade II). Clade I consists of multiple subclades: An. halophylus + An. triannulatus C; trans-Andean Venezuela; central Amazonia + central Bolivia; Atlantic coastal lowland; and Amazon delta. Clade II includes three subclades: Panama; cis-Andean Colombia; and cis-Venezuela. The Amazon delta specimens are in both clades, likely indicating local sympatry. Spatial and molecular variance analyses detected nine groups, corroborating some of subclades obtained in the combined data analysis.

CONCLUSION

Combination of the three molecular markers provided the best resolution for differentiation within An. triannulatus s.s. and An. halophylus and C. The latest two species seem to be very closely related and the analyses performed were not conclusive regarding species differentiation. Further studies including new molecular markers would be desirable to solve this species status question. Besides, results of the study indicate a trans-Andean origin for An. triannulatus s.l. The potential implications for malaria epidemiology remain to be investigated.

Population structure, mitochondrial polyphyly and the repeated loss of human biting ability in anopheline mosquitoes from the southwest Pacific

Australia and New Guinea contain high levels of endemism and biodiversity, yet there have been few evaluations of population-level genetic diversity in fauna occurring throughout the Australo-Papuan region. Using extensive geographical sampling, we examined and compared the phylogenetic relationships, phylogeography and population structure of Anopheles farauti, An. hinesorum and An. irenicus throughout their ranges in the southwest Pacific using mitochondrial (mtDNA COI) and nuclear (ribosomal protein S9 and ribosomal DNA ITS2) loci. Phylogenetic analyses suggest that the ability to utilize humans as hosts has been lost repeatedly, coincident with independent colonizations of the Solomon Islands. As some of the species under investigation transmit malaria in the region, this is a medically important finding. Maximum likelihood and Bayesian phylogenetic analyses of nuclear loci also showed that the three species are monophyletic. However, putative introgression of An. hinesorum mtDNA onto a nuclear background of An. farauti was evident in populations from Queensland, Torres Strait and southern New Guinea. Haplotype networks and pairwise F(ST) values show that there is significant genetic structure within New Guinea and Australia in both An. farauti and An. hinesorum, consistent with a long-term history of low gene flow among populations.

Comparative evolution of S7 Intron 1 and ribosomal internal transcribed spacer in Coilia nasus (Clupeiformes: Engraulidae)

Coilia nasus is widely distributed in the Yangtze River, the coastal waters of China, Korea and the Ariake Sound of Japan. Several ecotypes exist and this provides a useful model for the study of comparative diversity between molecular markers. Here we analyze and compare the nucleotide sequences between single-copy ribosomal protein S7 gene intron 1 (rpS7) and multiple-copy ribosomal internal transcribed spacer 1 (ITS1) in this species to compare the phylogenetic signal of the two nuclear genes. Nucleotide substitutions among the two gene sequences and partial sequence of mitochondrial cytochrome c oxidase subunit I (COI) gene were also analyzed. A total of 115 clones for rpS7 and 122 clones for ITS1 were obtained from 37 specimens. The nucleotide sequence length is 741 to 743 bp for rpS7 and 334 to 348 bp for ITS1. Intra- and inter-specimen variation in rpS7 results from nucleotide substitution, while such variation in ITS1 is mainly due to different numbers of short base repeats. The content of G + C is lower in rpS7 (43.5%) than in ITS1 (68.2%). Our results indicate that the proportion of the sequence variable sites is higher in rpS7 (61) than in ITS1 (23); the informative parsimony of rpS7 is evidently higher than that of ITS1 (26 vs. 2); the overall ratio between transitions and transversions in ITS1 is slightly lower than in rpS7, but remarkably lower than in COI. These results suggest that rpS7 is more suitable than ITS1 as a marker for genetic divergence of this group. Furthermore, gene flow is observed between the different geographic populations of C. nasus from the phylogeny of this species based on rpS7, showing that rpS7 has more evolutionary characteristics for understanding the processes of genomic evolution at the intraspecific level.

The dominant Anopheles vectors of human malaria in the Asia-Pacific region: occurrence data, distribution maps and bionomic précis

BACKGROUND

The final article in a series of three publications examining the global distribution of 41 dominant vector species (DVS) of malaria is presented here. The first publication examined the DVS from the Americas, with the second covering those species present in Africa, Europe and the Middle East. Here we discuss the 19 DVS of the Asian-Pacific region. This region experiences a high diversity of vector species, many occurring sympatrically, which, combined with the occurrence of a high number of species complexes and suspected species complexes, and behavioural plasticity of many of these major vectors, adds a level of entomological complexity not comparable elsewhere globally. To try and untangle the intricacy of the vectors of this region and to increase the effectiveness of vector control interventions, an understanding of the contemporary distribution of each species, combined with a synthesis of the current knowledge of their behaviour and ecology is needed.

RESULTS

Expert opinion (EO) range maps, created with the most up-to-date expert knowledge of each DVS distribution, were combined with a contemporary database of occurrence data and a suite of open access, environmental and climatic variables. Using the Boosted Regression Tree (BRT) modelling method, distribution maps of each DVS were produced. The occurrence data were abstracted from the formal, published literature, plus other relevant sources, resulting in the collation of DVS occurrence at 10116 locations across 31 countries, of which 8853 were successfully geo-referenced and 7430 were resolved to spatial areas that could be included in the BRT model. A detailed summary of the information on the bionomics of each species and species complex is also presented.

CONCLUSIONS

This article concludes a project aimed to establish the contemporary global distribution of the DVS of malaria. The three articles produced are intended as a detailed reference for scientists continuing research into the aspects of taxonomy, biology and ecology relevant to species-specific vector control. This research is particularly relevant to help unravel the complicated taxonomic status, ecology and epidemiology of the vectors of the Asia-Pacific region. All the occurrence data, predictive maps and EO-shape files generated during the production of these publications will be made available in the public domain. We hope that this will encourage data sharing to improve future iterations of the distribution maps.

Evolution of the courtship display of Melittobia (Hymenoptera: Eulophidae)

Melittobia, a polyphagous pupal parasitoid of many solitary bees and wasps, has a complex courtship display, the evolutionary history of which is poorly understood. The current phylogeny of this genus suggests the complexity of the courtship has increased over time. We test this hypothesis using Bayesian and Maximum likelihood techniques to construct a molecular phylogeny of Melittobia-the first for this genus using modern techniques--using two intergenic regions (Internal Transcribed Spacer regions 1 and 2) and Cytochrome Oxidase 1 (CO1), and reconstruct the ancestral state of each major node in the phylogeny for five behaviors coincidental to major shifts in display complexity. We show that species-groups identified solely on behavioral characters are reliable, although the relationships among groups are different than those previously hypothesized. We also show that behaviors traditionally used in the courtship studies do not support a hypothesis of increasing display complexity over time.

Evaluation of a PCR-RFLP-ITS2 assay for discrimination of Anopheles species in northern and western Colombia

Anopheles mosquitoes are routinely identified using morphological characters of the female that often lead to misidentification due to interspecies similarity and intraspecies variability. The aim of this work was to evaluate the applicability of a previously developed PCR-RFLP-ITS2 assay for accurate discrimination of anophelines in twelve localities spanning three Colombian malaria epidemiological regions: Atlantic Coast, Pacific Coast, and Uraba-Bajo Cauca-Alto Sinu region. The evaluation of the stability of the PCR-RFLP patterns is required since variability of the ITS2 has been documented and may produce discrepancies in the patterns previously reported. The assay was used to evaluate species assignation of 939 mosquitoes identified by morphology. Strong agreement between the morphological and molecular identification was found for species Anopheles albimanus, Anopheles aquasalis, Anopheles darlingi and Anopheles triannulatus s.l. (p≥0.05, kappa=1). However, disagreement was found for species Anopheles nuneztovari s.l., Anopheles neomaculipalpus, Anopheles apicimacula and Anopheles punctimacula (p≤0.05; kappa ranging from 0.33 to 0.80). The ITS2-PCR-RFLP assay proved valuable for discriminating anopheline species of northern and western Colombia, especially those with overlapping morphology in the Oswaldoi Group.

Incomplete concerted evolution and reproductive isolation at the rDNA locus uncovers nine cryptic species within Anopheles longirostris from Papua New Guinea

Background

Nuclear ribosomal DNA (rDNA) genes and transcribed spacers are highly utilized as taxonomic markers in metazoans despite the lack of a cohesive understanding of their evolution. Here we follow the evolution of the rDNA second internal transcribed spacer (ITS2) and the mitochondrial DNA cytochrome oxidase I subunit in the malaria mosquito Anopheles longirostris from Papua New Guinea (PNG). This morphospecies inhabits a variety of ecological environments indicating that it may comprise a complex of morphologically indistinguishable species. Using collections from over 70 sites in PNG, the mtDNA was assessed via direct DNA sequencing while the ITS2 was assessed at three levels - crude sequence variation through restriction digest, intragenomic copy variant organisation (homogenisation) through heteroduplex analysis and DNA sequencing via cloning.

Results

Genetic evaluation of over 300 individuals revealed that A. longirostris comprises eight ITS2 PCR-RFLP genotypes and nine ITS2 heteroduplex genotypes showing distinct copy variant organization profiles after PCR amplification. Seven of these nine genotypes were found to be sympatric with other genotypes. Phylogenetic analysis of cloned ITS2 PCR products and mtDNA COI confirmed all nine clades with evidence of reproductive isolation at the rDNA locus. Compensatory base changes in the ITS2 secondary structure or in pseudoknots were absent when closely related species were assessed. Individuals from each ITS2 genotype showed the same copy variant heteroduplex profile suggesting that the rDNA array is fixed within each genotype.

Conclusion

The centromere-proximal position of the rDNA array in Anopheles mosquitoes has probably reduced interchromosomal recombination leaving intrachromosomal events responsible for the observed pattern of concerted evolution we see in these mosquitoes. The stability of these intragenomic ITS2 copy variants within individuals and interbreeding populations suggests that rDNA is moving as a single evolutionary unit through natural populations to fixation and has provided a complementary diagnostic tool to the restriction digest for studying genetic discontinuities and species boundaries. In this, the utility of the ITS2 as a universal taxonomic marker is probably contingent on several factors pertaining to spacer dimensions and the genomic location of the rDNA array with respect to recombination and proximity to regions potentially under selection.

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.

Evolution of the nocturnal Nearctic Sphaeropthalminae velvet ants (Hymenoptera: Mutillidae) driven by Neogene orogeny and Pleistocene glaciation

The influence of historical climatic and geological changes on patterns of species diversification was investigated in a widely distributed group of North American nocturnal mutillids (Hymenoptera: Mutillidae: Sphaeropthalminae), with particular focus on Pleistocene glacial cycles and earlier patterns of Neogene mountain building. We collected molecular data from two nuclear intergenic regions (internal transcribed spacer regions 1 and 2; approximately 2600 bp in total) to produce the first phylogeny of nocturnal Nearctic mutillids. Bayesian inference of the combined data returned a well-resolved tree with posterior probabilities of over 95% for most nodes. This tree suggested the monophyly of the nocturnal, primarily Nearctic, Sphaeropthalminae genera, but the paraphyly of the three largest genera (Odontophotopsis, Photomorphus and Sphaeropthalma). Dates of species divergences were obtained using r8s (PL and NPRS) and BEAST with the date of Dominican amber set at three different dates (15 Ma, 20 Ma, and 45 Ma) to account for uncertainty in the fossil age. The derived dates ranged from the Pleistocene to the Middle Miocene, suggesting that both recent Pleistocene glaciation cycles and older orogenic events, albeit to a somewhat greater extent, were both causes of major diversification in western North America. Examination of other phylogeographical studies using North American desert taxa indicated that diversification patterns are explained by either mountain building or Pleistocene climate change, depending on the taxa in question.

Implementation of a novel PCR based method for detecting malaria parasites from naturally infected mosquitoes in Papua New Guinea

Background

Detection of Plasmodium species in mosquitoes is important for designing vector control studies. However, most of the PCR-based detection methods show some potential limitations. The objective of this study was to introduce an effective PCR-based method for detecting Plasmodium vivax and Plasmodium falciparum from the field-caught mosquitoes of Papua New Guinea.

Methods

A method has been developed to concurrently detect mitochondrial cytochrome b (Cyt b) of four human Plasmodium species using PCR (Cytb-PCR). To particularly discriminate P. falciparum from P. vivax, Plasmodium ovale and Plasmodium malariae, a polymerase chain reaction-repeated fragment length polymorphism (PCR-RFLP) has further been developed to use with this method. However, due to limited samples number of P. ovale and P. malariae; this study was mainly confined to P. vivax and P. falciparum. The efficiency of Cytb-PCR was evaluated by comparing it with two 'gold standards' enzyme linked immunosorbent assay specific for circumsporozoite protein (CS-ELISA) using artificially infected mosquitoes; and nested PCR specific for small subunit ribosomal RNA (SSUrRNA) using field caught mosquitoes collected from three areas (Kaboibus, Wingei, and Jawia) of the East Sepic Province of Papua New Guinea.

Results

A total of 90 mosquitoes were artificially infected with three strains of Plasmodium: P. vivax-210 (n = 30), P. vivax-247 (n = 30) and P. falciparum (n = 30). These infected mosquitoes along with another 32 unfed mosquitoes were first checked for the presence of Plasmodium infection by CS-ELISA, and later the same samples were compared with the Cytb-PCR. CS-ELISA for P. vivax-210, P. vivax-247 and P. falciparum detected positive infection in 30, 19 and 18 mosquitoes respectively; whereas Cytb-PCR detected 27, 16 and 16 infections, respectively. The comparison revealed a close agreement between the two assays (κ = 0.862, 0.842 and 0.894, respectively for Pv-210, Pv-247 and P. falciparum groups). It was found that the eight CS-ELISA-positive mosquitoes detected negative by Cytb-PCR were false-positive results. The lowest detection limit of this Cytb-PCR was 10 sporozoites. A highly concordance result was also found between nested PCR and Cytb-PCR using 107 field caught mosquitoes, and both tests concordantly detected P. falciparum in an Anopheles punctulatus mosquito collected from Kaboibus. Both tests thus suggested an overall sporozoite rate of 0.9% (1/107) in the study areas. Subsequently, PCR-RFLP efficiently discriminated P. falciparum from P. vivax for all of the Cytb-PCR positive samples.

Conclusion

A single step PCR based method has been introduced here that is highly sensitive, efficient and reliable for identifying P. vivax and P. falciparum from mosquitoes. The reliability of the technique was confirmed by its ability to detect Plasmodium as efficiently as those of CS-ELISA and nested PCR. Application of the assay offers the opportunity to detect vector species of Papua New Guinea and may contribute for designing further vector control programmes.

Malaria vectors of Papua New Guinea

Understanding malaria transmission in Papua New Guinea (PNG) requires exact knowledge of which Anopheles species are transmitting malaria and is complicated by the cryptic species status of many of these mosquitoes. To identify the malaria vectors in PNG we studied Anopheles specimens from 232 collection localities around human habitation throughout PNG (using CO(2) baited light traps and human bait collections). A total of 22,970mosquitoes were individually assessed using a Plasmodium sporozoite enzyme-linked immunosorbent assay to identify Plasmodiumfalciparum, Plasmodiumvivax and Plasmodiummalariae circumsporozoite proteins. All mosquitoes were identified to species by morphology and/or PCR. Based on distribution, abundance and their ability to develop sporozoites, we identified five species as major vectors of malaria in PNG. These included: Anophelesfarauti, Anopheleshinesorum (incriminated here, to our knowledge, for the first time), Anophelesfarauti 4, Anopheleskoliensis and Anophelespunctulatus. Anopheleslongirostris and Anophelesbancroftii were also incriminated in this study. Surprisingly, An. longirostris showed a high incidence of infections in some areas. A newly identified taxon within the Punctulatus Group, tentatively called An. farauti 8, was also found positive for circumsporozoite protein. These latter three species, together with Anopheleskarwari and Anophelessubpictus, incriminated in other studies, appear to be only minor vectors, while Anophelesfarauti 6 appears to be the major vector in the highland river valleys (>1500m above sea level). The nine remaining Anopheles species found in PNG have been little studied and their bionomics are unknown; most appear to be uncommon with limited distribution and their possible role in malaria transmission has yet to be determined.