Molecular phylogenetics of blackflies of the Simulium damnosum complex and cytophylogenetic implications

Cytotaxonomic characterization and estimation of migration patterns of onchocerciasis vectors (Simulium damnosum sensu lato) in northwestern Ethiopia based on RADSeq data

BACKGROUND

While much progress has been made in the control and elimination of onchocerciasis across Africa, the extent to which vector migration might confound progress towards elimination or result in re-establishment of endemism in areas where transmission has been eliminated remains unclear. In Northern Ethiopia, Metema and Metekel-two foci located near the Sudan border-exhibit continuing transmission. While progress towards elimination has been faster in Metema, there remains a problematic hotspot of transmission. Whether migration from Metekel contributes to this is currently unknown.

METHODOLOGY/PRINCIPLE FINDINGS

To assess the role of vector migration from Metekel into Metema, we present a population genomics study of 151 adult female vectors using 47,638 RADseq markers and mtDNA CoI sequencing. From additional cytotaxonomy data we identified a new cytoform in Metema, closely related to S. damnosum s.str, here called the Gondar form. RADseq data strongly indicate the existence of two distinctly differentiated clusters within S. damnosum s.l.: one genotypic cluster found only in Metema, and the second found predominantly in Metekel. Because blackflies from both clusters were found in sympatry (in all four collection sites in Metema), but hybrid genotypes were not detected, there may be reproductive barriers preventing interbreeding. The dominant genotype in Metema was not found in Metekel while the dominant genotype in Metekel was found in Metema, indicating that (at the time of sampling) migration is primarily unidirectional, with flies moving from Metekel to Metema. There was strong differentiation between clusters but little genetic differentiation within clusters, suggesting migration and gene flow of flies within the same genetic cluster are sufficient to prevent genetic divergence between sites.

CONCLUSIONS/SIGNIFICANCE

Our results confirm that Metekel and Metema represent different transmission foci, but also indicate a northward movement of vectors between foci that may have epidemiological importance, although its significance requires further study.

Testing the susceptibility of larval stages of Simulium to temephos and Bacillus thuringiensis var israelensis in Germany and Northern Cameroon

Assays to evaluate the susceptibility of Simulium larvae to temephos and Bacillus thuringiensis var israelensis (Bti) were carried out by setting-up an in vitro laboratory test ('bio assay') and a semi-natural test ('système de goutières') to assess the LC50/LC90 values. Larvae of Simulium species in Cameroon (S. damnosum s.l., S. hargreavesi, S. vorax and S. cervicornutum) and (S. (Odagmia) ornatum and S. latipes) in Germany were identified and tested. In the bio-assay, 50 larvae were exposed for 10 min to concentrations from 0.01 to 10 ppm. For the Simulium from Germany, the LC50 (LC90) values after 3 and 6 h were 3.1 (27.9) and 0.14 (1.26) ppm for temephos and for Bti 7.8 (70.2) and 1.7 (15.3) ppm, respectively. For Cameroonian species, the values of LC50 (LC90) were lower, that is, 0.42 (8.04), 0.14 (2.70) and 0.073 (1.38) ppm, respectively, after 3, 6 and 12 h for temephos. In a semi natural condition, the LC50 of 10 min of application of temephos was 0.84 ppm after 3 h and a working solution (2.6 L) of Bti killed 50% after 6 h. To detect an upcoming of any resistance as it happened in Ivory Coast, a study of the occurrence resistance genes should be implemented.

Genomic Epidemiology in Filarial Nematodes: Transforming the Basis for Elimination Program Decisions

Onchocerciasis and lymphatic filariasis are targeted for elimination, primarily using mass drug administration at the country and community levels. Elimination of transmission is the onchocerciasis target and global elimination as a public health problem is the end point for lymphatic filariasis. Where program duration, treatment coverage, and compliance are sufficiently high, elimination is achievable for both parasites within defined geographic areas. However, transmission has re-emerged after apparent elimination in some areas, and in others has continued despite years of mass drug treatment. A critical question is whether this re-emergence and/or persistence of transmission is due to persistence of local parasites-i.e., the result of insufficient duration or drug coverage, poor parasite response to the drugs, or inadequate methods of assessment and/or criteria for determining when to stop treatment-or due to re-introduction of parasites via human or vector movement from another endemic area. We review recent genetics-based research exploring these questions in Onchocerca volvulus, the filarial nematode that causes onchocerciasis, and Wuchereria bancrofti, the major pathogen for lymphatic filariasis. We focus in particular on the combination of genomic epidemiology and genome-wide associations to delineate transmission zones and distinguish between local and introduced parasites as the source of resurgence or continuing transmission, and to identify genetic markers associated with parasite response to chemotherapy. Our ultimate goal is to assist elimination efforts by developing easy-to-use tools that incorporate genetic information about transmission and drug response for more effective mass drug distribution, surveillance strategies, and decisions on when to stop interventions to improve sustainability of elimination.

Evidence of multiple colonizations as a driver of black fly diversification in an oceanic island

True oceanic islands typically host reduced species diversity together with high levels of endemism, which make these environmental set-ups ideal for the exploration of species diversification drivers. In the present study, we used black fly species (Diptera: Simuliidae) from Reunion Island as a model to highlight the main drivers of insect species diversification in this young and remote volcanic island located in the Southwestern Indian Ocean. Using local and regional (Comoros and Seychelles archipelagos) samples as well as specimens from continental Africa, we tested the likelihood of two distinct scenarios, i.e. multiple colonizations vs. in-situ diversification. For this, posterior odds were used to test whether species from Reunion did form a monophyletic group and we estimated divergence times between species. Three out of the four previously described Reunion black fly species could be sampled, namely Simulium ruficorne, Simulium borbonense and Simulium triplex. The phylogenies based on nuclear and mitochondrial markers showed that S. ruficorne and S. borbonense are the most closely related species. Interestingly, we report a probable mitochondrial introgression between these two species although they diverged almost six million years ago. Finally, we showed that the three Reunion species did not form a monophyletic group, and, combined with the molecular datation, the results indicated that Reunion black fly diversity resulted from multiple colonization events. Thus, multiple colonizations, rather than in-situ diversification, are likely responsible for an important part of black fly diversity found on this young Darwinian island.

The blackfly vectors and transmission of Onchocerca volvulus in Mahenge, south eastern Tanzania

The Mahenge Mountains onchocerciasis focus in south eastern Tanzania was historically one of the most heavily infected areas in the country. The vectors of Onchocerca volvulus are mainly Simulium damnosum complex blackflies, but a species of the Simulium neavei group may also contribute to transmission in some areas. The only detailed studies of parasite transmission in Mahenge were conducted in the late 1960s. The taxonomy of the S. damnosum complex has since been revised and onchocerciasis control through annual community directed treatment with ivermectin (CDTI) commenced in 1997. This study aimed to provide a cytogenetic and molecular update of the S. damnosum complex cytoforms present in Mahenge, and to evaluate the current status of O. volvulus transmission by blackflies following 19 years of annual CDTI. Rivers were surveyed to identify sites of S. damnosum s.l. breeding among the eastern slopes of the mountains, and human landing collections of adult female blackflies were made close to breeding sites. Identification of S. damnosum complex cytoforms was by cytotaxonomy of late-instar larvae and ITS1 amplicon size polymorphisms of larvae and adults. Adult blackflies were pool screened for O. volvulus infection using a triplex real-time PCR. The cytoforms 'Nkusi', Simulium kilibanum and 'Turiani' were found breeding in perennial rivers. 'Nkusi' and S. kilibanum were collected on human bait at 7/7 catch sites and possessed ITS1 profiles most closely resembling the molecular forms 'Nkusi J' and S. kilibanum 'T'. Whereas 'Turiani' was present in rivers, it was not collected on human bait and appears to be zoophilic. Simulium nyasalandicum was collected in low numbers on human bait at 3/7 catch sites. In total, 12,452 S. damnosum s.l. were pool screened and O. volvulus infection was detected in 97/104 pools of bodies and 51/104 pools of heads. The estimated percentage of S. damnosum s.l. carrying infective L3 stage parasites was 0.57% (95% CI 0.43%-0.74%). Onchocerca volvulus transmission by S. damnosum s.l. is continuing in the Mahenge Mountains after 19 years of annual CDTI. Infection rates appear similar to those reported in the 1960s, but a more detailed study is required to fully understand the epidemiological significance of the ongoing transmission. These results provide further evidence that annual CDTI may be insufficient to eliminate the parasite in formerly hyperendemic foci.

Insights Into Onchocerca volvulus Population Biology Through Multilocus Immunophenotyping

We have developed a serologically based immunophenotyping approach to study Onchocerca volvulus (Ov) population diversity. Using genomic sequence data and polymerase chain reaction-based genotyping, we identified nonsynonymous single-nucleotide polymorphisms (SNPs) in the genes of 16 major immunogenic Ov proteins: Ov-CHI-1/Ov-CHI-2, Ov16, Ov-FAR-1, Ov-CPI-1, Ov-B20, Ov-ASP-1, Ov-TMY-1, OvSOD1, OvGST1, Ov-CAL-1, M3/M4, Ov-RAL-1, Ov-RAL-2, Ov-ALT-1, Ov-FBA-1, and Ov-B8. We assessed the immunoreactivity of onchocerciasis patient sera (n = 152) from the Americas, West Africa, Central Africa, and East Africa against peptides derived from 10 of these proteins containing SNPs. Statistically significant variation in immunoreactivity among the regions was seen in SNP-containing peptides derived from 8 of 10 proteins tested: OVOC1192(1-15), OVOC9988(28-42), OVOC9225(320-334), OVOC7453(22-36), OVOC11517(14-28), OVOC3177(283-297), OVOC7911(594-608), and OVOC12628(174-188). Our data show that differences in immunoreactivity to variant antigenic peptides may be used to characterize Ov populations, thereby elucidating features of Ov population biology previously inaccessible because of the limited availability of parasite material.

Freshwater biodiversity and aquatic insect diversification

Inland waters cover less than 1% of Earth's surface but harbor more than 6% of all insect species: Nearly 100,000 species from 12 orders spend one or more life stages in freshwater. Little is known about how this remarkable diversity arose, although allopatric speciation and ecological adaptation are thought to be primary mechanisms. Freshwater habitats are highly susceptible to environmental change and exhibit marked ecological gradients. Standing waters appear to harbor more dispersive species than running waters, but there is little understanding of how this fundamental ecological difference has affected diversification. In contrast to the lack of evolutionary studies, the ecology and habitat preferences of aquatic insects have been intensively studied, in part because of their widespread use as bioindicators. The combination of phylogenetics with the extensive ecological data provides a promising avenue for future research, making aquatic insects highly suitable models for the study of ecological diversification.

The sandflies (Diptera: Psychodidae, Phlebotominae) in military camps in northern Afghanistan (2007-2009), as identified by morphology and DNA 'barcoding'

As part of a continuous, standardized programme of monitoring the Leishmania vectors in German military camps in northern Afghanistan between 2007 and 2009, a detailed taxonomic analysis of the endemic sandfly fauna, as sampled using light and odour-baited traps, was conducted. Of the 10 sandfly species that were recorded, six may serve as enzootic and/or zooanthroponotic vectors of parasites causing human leishmaniasis. The use of a simple DNA-'barcoding' technique based on the mitochondrial cyt b gene, to identify the collected sandflies to species level, revealed (1) a clear discrimination between the potential vector species, (2) clustering of species within most subgenera, and (3) particularly high heterogeneity within the subgenus Paraphlebotomus (Phlebotomus alexandri being grouped with Ph. papatasi rather than with other Paraphlebotomus species). The data also indicate a high level of genetic heterogeneity within the subgenus Sergentomyia but close similarity between Sergentomyia sintoni and Sergentomyia murgabiensis. The morphological similarity of many medically important sandflies can make species identification difficult, if not impossible. The new DNA-barcoding techniques may provide powerful discriminatory tools in the future.

Investigations into the isolation of the Tukuyu focus of onchocerciasis (Tanzania) from S. damnosum s.l. vector re-invasion

As part of the feasibility study for an onchocerciasis vector elimination project we investigated the isolation of the Tukuyu focus in Tanzania from possible vector re-invasion. This was achieved by examining the distribution of the Simulium damnosum complex vector cytospecies outside the focus to look for potential sources of re-invasion. Besides cytotaxonomic identifications of the aquatic stages, we applied morphotaxonomic and molecular techniques to identify S. thyolense and confirm it as the anthropophilic species in both the Tukuyu and the neighbouring Ruvuma foci. We detected significant differences in chromosome inversion frequencies between the Tukuyu populations and those breeding to the southwest in the adjacent Songwe river basin and in northern Malawi (where there is no man-biting and no onchocerciasis), suggesting that there is not normally a great deal of migration in either direction. By contrast, populations of S. thyolense from the Tukuyu and Ruvuma foci (150km southeast of Tukuyu) were much more similar in terms of their chromosomal polymorphisms, indicating a higher possibility of re-invasion, although migration is still restricted to some extent, as indicated by some differences in chromosome polymorphisms between the two foci. Future migratory events which might be associated with vector control operations can be monitored by vector cytospecies identification, the frequency of polymorphic inversions which characterise the different vector populations, and the identification of accompanying non-vector cytospecies (e.g. S. plumbeum and cytotype Kasyabone occur exclusively in the two foci, and hence their re-appearance in Tukuyu could have only one outside source). The morphology of the scutal pattern of neonate males may act as a quick test for vector species identification where chromosome squashes are unavailable.

Construction and characterisation of a BAC library made from field specimens of the onchocerciasis vector Simulium squamosum (Diptera: Simuliidae)

A Bacterial Artificial Chromosome (BAC) library was made from wild-caught Simulium squamosum, which is an important vector of human onchocerciasis. The library is composed of 12,288 BACs, with an average insert size of 128 kb, and is expected to contain ~1.54 GB of cloned DNA. Random BAC-end sequencing generated over 95 kb of DNA sequence data from which putative S. squamosum gene sequences and novel repetitive DNA families were identified, including DNA transposons, retrotransposons and simple sequence repeats (SSRs). The sequence survey also provided evidence of DNA of microbial origin, and dissection of sample blackflies indicated that some of those used to prepare the library were likely to be parasitized by the mermithid Isomermis lairdi. Hybridisations with a set of three independent blackfly single-copy genes and two Wolbachia genes suggest that the library provides around 13-fold coverage of the S. squamosum genome and about 12-fold coverage of its Wolbachia endosymbiont.

Suggesting synonymies? Comments on Kiehl et al. (2009) "the European vectors of Bluetongue virus: are there species complexes, single species or races in Culicoides obsoletus and C. pulicaris detectable by sequencing ITS-1, ITS-2 and 18S-rDNA?"

Species recognition and identification are crucial in any biological studies, especially when dealing with insect species involved in pathogen transmission. In recent years, molecular approaches have helped the clarification of systematic schemes and taxonomic status. Kiehl et al. (Parasitol Res 105:331-336, 2009) used molecular data to discuss the taxonomic status of biting midge species in the Palaearctic region. In the present work, the statements that "[Thus] there is no molecular support for the existence of a separate species C. montanus" and "[Therefore] probably C. scoticus should be considered only as a race of C. obsoletus" are discussed.

The species delimitation problem in the Simulium damnosum complex, blackfly vectors of onchocerciasis

The Simulium damnosum Theobald complex (Diptera: Simuliidae) comprises 57 cytoforms grouped into six subcomplexes. Previous phylogenetic studies using gene sequences have not completely resolved the evolutionary relationships of the cytoforms. The present study investigated the systematics of the complex using a phylogeographic approach. The differentiation between eastern and western forms observed in the phylogenetic studies is confirmed in the estimated haplotype networks. However, haplotypes tend to group in geographical clades and not according to cytoforms. Spatial analyses of the molecular variance also resulted in optimal groupings of sequences that did not correspond to cytoform boundaries. Moreover, Mantel tests showed significant correlations, although not strong, between genetic and geographical distances. This suggests an isolation-by-distance model of differentiation. Furthermore, there are instances in which genetic differentiation between cytoforms is low and not significant. These results indicate a lack of clear genetic differentiation between the cytoforms, which may be explained either by a separation of the taxa recent enough to allow the accumulation of few genetic differences or by recombination between the genomes of the cytoforms, which may be the result of hybridization with introgression or of non-independent evolutionary lineages. The results also emphasize the need for further sampling and for the use of more variable markers in order to clarify the evolutionary history of the group.

Confirmation of the species status of the blackfly Simulium galeratum in Britain using molecular taxonomy

Since 1920 Simulium reptans (Linnaeus) (Diptera: Simuliidae) has been reported as exhibiting two different larval morphotypes, a typical S. reptans and an atypical S. reptans var. galeratum, which differ in the markings of the larval head capsule. Inconsistent variation in adults and no apparent variation in the pupae have led taxonomists to conclude that these types in Britain are a single species. We investigated populations in Britain where either the typical form or var. galeratum is found, and one population where the two exist sympatrically. A phylogenetic study based upon a region of the mitochondrial cytochrome c oxidase 1 gene (DNA barcoding) produced a tree that delineated the morphotypes into two distinct monophyletic clades. The average Kimura-2-parameter distances within each clade (i.e. within each morphotype) were very low (0.67% and 0.78%), with the distances between morphotypes being 9-10-fold greater (mean 7.06%). This is concordant with differences within and between species in other taxa; based upon the strict correlation between the molecular variation and the morphotypes, we propose the re-instatement of S. galeratum to species status.

Taxonomy and inventory of the cytospecies and cytotypes of the Simulium damnosum complex (Diptera: Simuliidae) in relation to onchocerciasis

We provide an inventory of all named cytoforms of the Simulium damnosum complex (including those which are now considered invalid), along with all inversions that have been recorded (including synonyms and homonyms). There are 55 valid and distinct cytoforms known from the S. damnosum complex making it the largest sibling species complex of any vectors, and probably of any insect or other animal. All cytoforms are listed along with their fixed and diagnostic inversions and country distribution. There are 183 inversions known from the complex as a whole, of which 49% are fixed and/or diagnostic between cytoforms, and the fixed/diagnostic inversions seem to occur disproportionately on chromosome arm 2L.

Cytogenetic and PCR-based identification of S. damnosum'Nkusi J' as the anthropophilic blackfly in the Uluguru onchocerciasis focus in Tanzania

A combined molecular and cytogenetic analysis was used to identify the local blackfly vectors of onchocerciasis in the Uluguru Mts. focus, Morogoro region, central Tanzania. Cytotaxonomic identification revealed the occurrence of three cytospecies of the Simulium damnosum complex, along altitudinal gradients. The cytoform "Nkusi" was found breeding at low to mid-altitudes (100-500 m), while S. kilibanum was found at mid- and high altitudes (300-1260 m), being sympatric with "Sebwe" only in the highlands. The analyses further revealed the cytoform "Nkusi" from the nearby Kilosa focus. Simulium kilibanum and "Nkusi" are potential vectors of onchocerciasis, while the "Sebwe" form is considered non-anthropophilic and therefore not a vector. By means of molecular typing of wild-caught females, "Nkusi" was incriminated as the man-biting fly at all catching sites, and consequently, is the assumed vector of onchocerciasis over the full altitude range. In the high altitude zones, however, S. neavei group species may act as secondary vectors at the edge of the forest reserve, if transmission actually occurs there at all. Our results also suggest that vector control could be feasible because of the vector species breeding in the foothills below 500 m in relatively accessible open country.