Morphological and molecular characterization of larval trematodes infecting the assassin snail genus Anentome in Thailand.

The assassin snail genus Anentome is widespread in Southeast Asia, and is distributed all over the world via the aquarium trade. One species of genus Anentome, Anentome helena, is known to act as intermediate host of parasitic trematodes. This study investigates the taxonomic diversity of larval trematodes infecting A. helena and Anentome wykoffi in Thailand. Larval trematodes were identified by combining morphological and DNA sequence data (cytochrome c oxidase I and internal transcribed spacer 2). Species delimitation methods were used to explore larval trematode species boundaries. A total of 1107 specimens of Anentome sp. were collected from 25 localities in Thailand. Sixty-two specimens of A. helena (n = 33) and A. wykoffi (n = 29) were infected by zoogonid cercariae, heterophyid metacercariae and echinostome metacercariae, with an overall prevalence of 5.6% (62/1107) and population-level prevalences in the range of 0.0-22.3%. DNA sequence data confirmed that the larval trematodes belong to the families Echinostomatidae, Heterophyidae and Zoogonidae. As such, this is the first report of zoogonid cercariae and heterophyid metacercariae in A. helena, and echinostome metacercariae in A. wykoffi. Moreover, this study provides evidence of tentative species-level differentiation between Thai Echinostoma sp. and Cambodian Echinostoma mekongi, as well as within Echinostoma caproni, Echinostoma trivolvis and Echinostoma revolutum.

DNA identification of species of the Anopheles maculipennis complex and first record of An. daciae in Belgium.

The present study aimed at identifying the members of the Anopheles maculipennis complex (Diptera: Culicidae) occurring in Belgium. Therefore, the second internal transcribed spacer of nuclear ribosomal DNA (ITS2) and the mitochondrial cytochrome oxidase subunit I (COI) loci were sequenced in 175 and 111 specimens, respectively, collected between 2007 and 2019. In parallel, the suitability of two species-diagnostic PCR-RFLP assays was tested. The identified specimens included: An. maculipennis s.s. (N = 105), An. daciae (N = 62), An. atroparvus (N = 6) and An. messeae (N = 2). Each species was characterized by unique ITS2 haplotypes, whereas COI only supported the monophyly of An. atroparvus, a historical malaria vector in Belgium. Species identification results were further supported by unique PCR-RFLP banding patterns. We report for the first time An. daciae in Belgium, where it was found to co-occur with An. maculipennis s.s. The latter was the most prevalent in the collection studied (60%) and appears to have the widest distribution in Belgium. As in other studies, An. daciae and An. messeae appeared the most closely related species, up to the point that their species status remains debatable, while their ecological differences, including vector competences, need further study.

Correction to: Highly polymorphic mitochondrial DNA and deceiving haplotypic differentiation: implications for assessing population genetic differentiation and connectivity.

...

Highly polymorphic mitochondrial DNA and deceiving haplotypic differentiation: implications for assessing population genetic differentiation and connectivity.

BACKGROUND: Hyperdiverse mtDNA with more than 5% of variable synonymous nucleotide sites can lead to erroneous interpretations of population genetic differentiation patterns and parameters (φST, DEST). We illustrate this by using hyperdiverse mtDNA markers to infer population genetic differentiation and connectivity in Melarhaphe neritoides, a NE Atlantic (NEA) gastropod with a high dispersal potential. We also provide a recent literature example of how mtDNA hyperdiversity may have misguided the interpretation of genetic connectivity in the crab Opecarcinus hypostegus. RESULTS: mtDNA variation surveyed throughout the NEA showed that nearly all M. neritoides specimens had haplotypes private to populations, suggesting at first glance a lack of gene flow and thus a strong population genetic differentiation. Yet, the bush-like haplotype network, though visually misleading, showed no signs of phylogeographic or other haplotype structuring. Coalescent-based gene flow estimates were high throughout the NEA, irrespective of whether or not mtDNA hyperdiversity was reduced by removing hypervariable sites. CONCLUSIONS: Melarhaphe neritoides seems to be panmictic over the entire NEA, which is consistent with its long-lived pelagic larval stage. With hyperdiverse mtDNA, the apparent lack of shared haplotypes among populations does not necessarily reflect a lack of gene flow and/or population genetic differentiation by fixation of alternative haplotypes (DEST ≈ 1 does not a fortiori imply φST ≈ 1), but may be due to (1) a too low sampling effort to detect shared haplotypes and/or (2) a very high mutation rate that may conceal the signal of gene flow. Hyperdiverse mtDNA can be used to assess connectivity by coalescent-based methods. Yet, the combined use of φST and DEST can provide a reasonable inference of connectivity patterns from hyperdiverse mtDNA, too.

Prevalence of Angiostrongylus vasorum, Aelurostrongylus abstrusus and Crenosoma vulpis larvae in native slug populations in Germany.

Metastrongyloid parasites represent sparsely studied parasites of dogs and cats in Germany. Recent European surveys indicate that these parasites are spreading in Europe. Actual data on prevalence of Angiostrongylus vasorum in dogs and foxes reveal several endemic foci in Germany. However, actual data on the prevalence of A. vasorum and other metastrongyloid lungworm larvae in a wide range of slug and snail intermediate hosts, such as Arion lusitanicus, are missing for Germany. To fill this gap, we conducted an epidemiological survey on native German slugs in selected regions of Hesse and Rhineland-Palatinate. The focus was on slugs, because in study areas slugs appear to be more abundant than snails. Slugs were collected throughout different seasons of the year in areas that were previously proven to be hyperendemic for A. vasorum fox infections. Overall, a total of 2701 slugs were collected and examined for lungworm larvae via artificial digestion. The number of A. vasorum larvae per slug varied considerably (1-546 larvae per specimen). Some hotspot areas with high A. vasorum prevalence in slugs (up to 19.4%) were identified. The overall A. vasorum prevalence varied with season with largest number of slugs infected in summer (9.1%) and lowest number in winter (0.8%). The current study revealed a total A. vasorum prevalence of 4.7% in slugs based on microscopic analyses. Confirmation of lungworm species was made by specific duplex-real-time PCRs. Hence, these data demonstrate that final hosts are at a permanent risk for A. vasorum infections during all seasons when living in investigated areas. Besides A. vasorum, other lungworm larvae were also detected, such as Crenosoma vulpis (the fox lungworm, 2.3%) and Aelurostrongylus abstrusus (feline lungworm, 0.2%).

Exploring the shell-based taxonomy of the Sri Lankan land snail Corilla H. and A. Adams, 1855 (Pulmonata: Corillidae) using mitochondrial DNA.

The land-snail genus Corilla is endemic to Sri Lanka and India's Western Ghats. The ten extant Sri Lankan species belong to two distinct shell forms that are associated respectively with lowland and montane rainforest. We here present the first molecular phylogenetic analysis for Corilla. Our dataset includes nine nominal Sri Lankan species and is based on three mitochondrial genes (CO1, ND1 and 16S). Although the deeper nodes in the trees are not fully resolved, the results do suggest speciation in Corilla has involved repeated, ecologically-driven convergence in shell form. The mtDNA data agree with the current shell-based taxonomy for C. adamsi, C. beddomeae, C. carabinata, C. humberti and C. colletti, consistently supporting the first four as monophyletic, and supporting the last also as monophyletic in nearly all analyses. Corilla adamsi, C. beddomeae and C. colletti may each contain at least one additional, previously undescribed species. The relationship between northern and eastern C. odontophora couldn't be reliably resolved, but our results suggest that they are distinct species and that there is further species-level or intraspecific (geographical) differentiation within eastern C. odontophora. The current, morphologically-defined species limits of the three remaining nominal species, C. gudei, C. erronea and C. fryae, are inconsistent with the mtDNA sequence data. Northern and southern C. gudei appear to be distinct species: the sister taxon of southern C. gudei is C. humberti, and most analyses showed that the sister taxon of northern C. gudei is the lowland C. carabinata. Corilla erronea and C. fryae constitute a well supported clade in which both nominal species are paraphyletic. While most intra-clade CO1 p-distances are moderate to relatively large, the phylogenetic structuring within the clade does not seem to correspond to any obvious morphological, elevational or geographical patterns. These results are difficult to interpret, and further detailed study is needed before the taxonomic status of C. erronea and C. fryae can be resolved.

Biomimetic and bio-inspired uses of mollusc shells.

Climate change and ocean acidification are likely to have a profound effect on marine molluscs, which are of great ecological and economic importance. One process particularly sensitive to climate change is the formation of biominerals in mollusc shells. Fundamental research is broadening our understanding of the biomineralization process, as well as providing more informed predictions on the effects of climate change on marine molluscs. Such studies are important in their own right, but their value also extends to applied sciences. Biominerals, organic/inorganic hybrid materials with many remarkable physical and chemical properties, have been studied for decades, and the possibilities for future improved use of such materials for society are widely recognised. This article highlights the potential use of our understanding of the shell biomineralization process in novel bio-inspired and biomimetic applications. It also highlights the potential for the valorisation of shells produced as a by-product of the aquaculture industry. Studying shells and the formation of biominerals will inspire novel functional hybrid materials. It may also provide sustainable, ecologically- and economically-viable solutions to some of the problems created by current human resource exploitation.

Identification of Belgian mosquito species (Diptera: Culicidae) by DNA barcoding.

Since its introduction in 2003, DNA barcoding has proven to be a promising method for the identification of many taxa, including mosquitoes (Diptera: Culicidae). Many mosquito species are potential vectors of pathogens, and correct identification in all life stages is essential for effective mosquito monitoring and control. To use DNA barcoding for species identification, a reliable and comprehensive reference database of verified DNA sequences is required. Hence, DNA sequence diversity of mosquitoes in Belgium was assessed using a 658 bp fragment of the mitochondrial cytochrome oxidase I (COI) gene, and a reference data set was established. Most species appeared as well-supported clusters. Intraspecific Kimura 2-parameter (K2P) distances averaged 0.7%, and the maximum observed K2P distance was 6.2% for Aedes koreicus. A small overlap between intra- and interspecific K2P distances for congeneric sequences was observed. Overall, the identification success using best match and the best close match criteria were high, that is above 98%. No clear genetic division was found between the closely related species Aedes annulipes and Aedes cantans, which can be confused using morphological identification only. The members of the Anopheles maculipennis complex, that is Anopheles maculipennis s.s. and An. messeae, were weakly supported as monophyletic taxa. This study showed that DNA barcoding offers a reliable framework for mosquito species identification in Belgium except for some closely related species.

Population genetic structure of the tree-hole tick Ixodes arboricola (Acari: Ixodidae) at different spatial scales.

The endophilic tick Ixodes arboricola infests cavity-nesting birds, and its dispersal strongly depends on the movements of its host. Population genetic structure of I. arboricola was studied with seven polymorphic microsatellite markers. We collected 268 ticks from 76 nest boxes in four woodlots near Antwerp, Belgium. These nest boxes are mainly used by the principal hosts of I. arboricola, the great tit Parus major and the blue tit Cyanistes caeruleus. As these birds typically return to the same cavity for roosting or breeding, ticks within nest boxes were expected to be highly related, and tick populations were expected to be spatially structured among woodlots and among nest boxes within woodlots. In line with the expectations, genetic population structure was found among woodlots and among nest boxes within woodlots. Surprisingly, there was considerable genetic variation among ticks within nest boxes. This could be explained by continuous gene flow from ticks from nearby tree holes, yet this remains to be tested. A pairwise relatedness analysis conducted for all pairs of ticks within nest boxes showed that relatedness among larvae was much higher than among later instars, which suggests that larvae are the most important instar for tick dispersal. Overall, tick populations at the studied spatial scale are not as differentiated as predicted, which may influence the scale at which host-parasite evolution occurs.

Isolation and characterization of ten polymorphic microsatellite loci in Ixodes arboricola, and cross-amplification in three other Ixodes species.

We characterized ten polymorphic microsatellite loci from the tree-hole tick, Ixodes arboricola. Loci were screened in 11-18 individuals from three Belgian populations and five to ten alleles were found at each locus. Seven loci did not show deviations from Hardy-Weinberg equilibrium conditions and there were no indications for null alleles at these loci. The three other loci showed significant heterozygote deficiencies in at least one population, and a high potential for the occurrence of null alleles. We observed no effect of potential host DNA on the scoring of the microsatellites. Cross-amplification of the microsatellites was tested in eight specimens of three congeneric species: I. ricinus, I. hexagonus and I. frontalis. Depending on the species, six or seven of the loci were amplified in ≥ 4 of the 8 specimens and were polymorphic in each of these species (except for Ixaf 11 in I. frontalis and I. ricinus). These loci thus provide a tool for population genetic analysis of I. arboricola. The suitability of these markers needs to be further investigated in its congeners.

Cryptic diversity and gene flow among three African agricultural pests: Ceratitis rosa, Ceratitis fasciventris and Ceratitis anonae (Diptera, Tephritidae).

The 'Ceratitis FAR complex' is a species complex of African fruit flies (Diptera, Tephritidae) including the major agricultural pest Ceratitis rosa and the morphologically similar Ceratitis fasciventris and Ceratitis anonae. To resolve their intra- and interspecific genetic relationships and to estimate gene flow within this complex, we surveyed allelic variation at 16 microsatellite loci in 27 African populations of the three morphospecies. Interpopulation genetic distances and individual Bayesian assignments distinguished five genotypic clusters: two involving C. rosa (R1, R2; that may occur in sympatry), two involving C. fasciventris (F1, F2; with parapatric distributions) and one involving C. anonae (A). Intra- and interspecific patterns of genetic differentiation were not hierarchically structured and genetic differentiation between conspecific clusters (F1-F2 and R1-R2) was higher or comparable with differentiation between heterospecific clusters (e.g. F1-A or R2-A). In some cases, gene flow estimates among morphospecies or among heterospecific genotypic clusters were significantly different from zero, showing the lack of reproductive isolation. Genetic differentiation between genotypic clusters was partly supported by morphological differences observed a posteriori in male secondary sexual characters. These results suggest important revisions to current models of ecological niche requirements and invasion risk of the major agricultural pest C. rosa and provide a basis for a taxonomic re-interpretation of the FAR complex.

A quantitative comparison of frugivorous tephritids (Diptera: Tephritidae) in tropical forests and rural areas of the Democratic Republic of Congo.

Most of the current knowledge about African tephritids originates from studies performed in agricultural areas, while information about their distribution in pristine or moderately disturbed environments is extremely scarce. This study aims at (i) describing levels of spatial variability of frugivorous tephritids in tropical forests and small rural villages of the Congo River basin and (ii) verifying if human-mediated activities, such as small-scale agriculture and trade, can affect their distribution patterns. Four locations were sampled along a 250 km stretch of the Congo River. At each location, pristine and disturbed habitats (i.e. tropical forests and small rural villages, respectively) were sampled, with three replicate sites in each combination of habitat and location. Sampling with modified McPhail traps baited with four different attractants yielded 819 tephritid specimens of 29 species from seven genera (Bactrocera, Carpophthoromyia, Ceratitis, Dacus, Celidodacus, Perilampsis, Trirhithrum). The three most abundant species sampled (Dacus bivittatus, D. punctatifrons, Bactrocera invadens) showed significant variations in abundance across locations and sites and accounted for 98.29% of the overall dissimilarity between habitats. Assemblages differed among locations and sites while they showed significant differences between pristine and disturbed habitats in two out of the four locations. This study shows that frugivorous tephritids in central Congo have remarkably patchy distributions with differences among locations and sites representing the main source of variability. Our data show that, in rural villages of central Democratic Republic of Congo, human activities, such as small-scale agriculture and local commerce, are not always sufficient to promote differences between the tephritid assemblages of villages and those of the surrounding tropical forests.

Parallel habitat specialization within the wolf spider genus Hogna from the Galápagos.

Within most island archipelagos, such as the Galápagos, similar ecological gradients are found on geographically isolated islands. Species radiations in response to these ecological gradients may follow different scenarios being (i) a single habitat specialization event followed by secondary colonization of each ecotype on the different islands or (ii) repeated and parallel habitat specialization on each island separately. This latter scenario has been considered less likely as gene flow might hamper such ecotypic differentiation. At least for the Galápagos, the extent to which this process is involved in species radiations remains yet poorly understood. Within the wolf spider genus Hogna, seven species are described that can be divided into three different ecotypes based on general morphology and habitat preference i.e. species that inhabit the pampa vegetation in the highlands, species that occur in coastal dry habitats and one generalist species. Comparison of the species phylogeny based on one mitochondrial (COI) and one nuclear (28S) gene fragment convincingly demonstrates that 'pampa' and 'coastal dry' species evolved in parallel on the islands Santa Cruz and San Cristóbal. Despite the observation that allozymes analysis indicated that each species forms a distinct genetic cluster, phylogenetic divergence within these species complexes was very low and paraphyletic and most likely due to hybridization rather than incomplete lineage sorting, as demonstrated for the Santa Cruz species complex. This suggests that within-island speciation occurred under low levels of gene flow. Species phylogeny in general did not follow the progression of island emergence as a molecular clock analysis suggested that island endemic species may have diverged after as well as before the emergence of the islands. This represents the first clear example of parallel and within-island speciation because of habitat specialization on the Galápagos and that such divergence most likely occurred under historic gene flow.

Macrogeographic population structuring in the cosmopolitan agricultural pest Bactrocera cucurbitae (Diptera: Tephritidae).

The macrogeographic population structure of the agricultural pest Bactrocera cucurbitae (Diptera: Tephritidae) was investigated in order to identify the geographic origin of the species and reconstruct its range expansion. Individuals of B. cucurbitae were collected from 25 worldwide-distributed localities (n = 570) and genotyped at 13 microsatellite loci. The Bayesian clustering reveals that B. cucurbitae can be subdivided into five main groups corresponding to populations from (i) the African continent, (ii) La Réunion, (iii) Central Asia, (iv) East Asia and (v) Hawaii. The proportions of inter-regional assignments and the higher values of genetic diversity in populations from Pakistan, India and Bangladesh suggest that B. cucurbitae originated in Central Asia and expanded its range to East Asia and Hawaii on one hand and to Africa and the islands of the Indian Ocean on the other. A number of outliers (10-19 specimens according to different clustering algorithms) show high levels of admixture (Q > 0.70) with populations from different regions and reveal complex patterns of inter-regional gene flow. Anthropogenic transport is the most plausible promoter of this large-scale dispersal. The introduction of individuals from geographically distant sources did not have a relevant role in the most recent African invasions, which originated from the expansion of local populations. These results could provide a useful background to better evaluate invasion risks and establish priorities for the management of this cosmopolitan agricultural pest.

Genetic divergence among sympatric colour morphs of the Dalmatian wall lizard (Podarcis melisellensis).

If alternative phenotypes in polymorphic populations do not mate randomly, they can be used as model systems to study adaptive diversification and possibly the early stages of sympatric speciation. In this case, non random mating is expected to support genetic divergence among the different phenotypes. In the present study, we use population genetic analyses to test putatively neutral genetic divergence (of microsatellite loci) among three colour morphs of the lizard Podarcis melisellensis, which is associated with differences in male morphology, performance and behaviour. We found weak evidence of genetic divergence, indicating that gene flow is somewhat restricted among morphs and suggesting possible adaptive diversification.

African Dacus (Diptera: Tephritidae: molecular data and host plant associations do not corroborate morphology based classifications.

The genus Dacus Fabricius includes economically important pest fruit flies distributed in the Afrotropical and Indo-Australian regions. Two recent revisions based on morphological characters proposed new and partially discordant classifications synonymizing/revalidating several subgeneric names and forming species groups. Regardless these efforts, the phylogenetic relationships among Dacus species remained largely unresolved mainly because of the difficulties in assigning homologous character states. Therefore we investigated the phylogeny of African Dacus by sequencing 71 representatives of 32 species at two mitochondrial (COI, 16S) and one nuclear (period) gene fragments. Phylogenetic relationships were inferred through Bayesian and Maximum Parsimony methods and hypotheses about the monophyly of Dacus subgenera were tested by Shimodaira-Hasegawa tests. The congruence tests and the analyses of the single gene fragments revealed that the nuclear gene supports similar conclusions as the two mitochondrial genes. Levels of intra- and inter-specific differentiation of Dacus species were highly variable and, in some cases, largely overlapping. The analyses of the concatenated dataset resolved two major bootstrap-supported groups as well as a number of well-supported clades and subclades that often comprised representatives of different subgenera. Additionally, specimens of Dacus humeralis from Eastern and Western African localities formed separate clades, suggesting cryptic differentiation within this taxon. The comparisons between the molecular phylogeny and the morphological classification revealed a number of discrepancies and, in the vast majority of cases, the molecular data were not compatible with the monophyly of the currently recognised subgenera. Conversely, the molecular data showed that Apocynaceae feeders are a monophyletic sister group of species feeding on both Cucurbitaceae and Passifloraceae (these latter being also monophyletic). These results show a clear association between the molecular phylogeny of African Dacus and the evolution of host plant choice and provide a basis towards a more congruent taxonomy of this genus.

Phylogeography of the common ragworm Hediste diversicolor (Polychaeta: Nereididae) reveals cryptic diversity and multiple colonization events across its distribution.

Previous studies on the common ragworm Hediste diversicolor (Polychaeta: Nereididae) revealed a marked genetic fragmentation across its distribution and the occurrence of sibling taxa in the Baltic Sea. These results suggested that the phylogeographic patterns of H. diversicolor could reflect interactions between cryptic differentiation and multiple colonization events. This study aims to describe the large-scale genetic structuring of H. diversicolor and to trace the phylogeographic origins of the genetic types described in the Baltic Sea. Samples of H. diversicolor (2 < n < 28) were collected at 16 locations across the NE Atlantic coasts of Europe and Morocco and in the Mediterranean, Black and Caspian Seas and sequenced at two mitochondrial gene fragments (COI and cytb, 345 and 290 bp, respectively). Bayesian analyses revealed deep phylogeographic splits yielding three main clades corresponding to populations (i) from the NE Atlantic coasts (from Germany to Morocco) and from part of the Western Mediterranean, (ii) from the Mediterranean Sea, and (iii) from the Black and Caspian Seas. These clades are further divided in well-supported subclades including populations from different regions of NE Atlantic and Mediterranean (i.e. Portugal/Morocco, Western Mediterranean, Adriatic Sea). The Baltic Sea comprises three sympatric lineages sharing a common evolutionary history with populations from NE Atlantic, Western Mediterranean and Black/Caspian Seas, respectively. Hence, the current patterns of genetic structuring of H. diversicolor appear as the result of allopatric isolation, multiple colonization events and possible adaptation to local environmental conditions.

Characterization of polymorphic microsatellite markers in the Dalmatian wall lizard Podarcis melisellensis (Squamata: Lacertidae).

We describe polymerase chain reaction primers and amplification conditions for 13 highly polymorphic microsatellite DNA loci isolated from the Dalmatian wall lizard, Podarcis melisellensis. The number of alleles per locus ranged from 12 to 41, with levels of observed heterozygosity between 0.62 and 0.94. Most of these loci were successfully cross-amplified in the closely related species P. sicula, but levels of polymorphism were always lower.

Phylogeography of the Rhabditis (Pellioditis) marina species complex: evidence for long-distance dispersal, and for range expansions and restricted gene flow in the northeast Atlantic.

Pinpointing processes that structure the geographical distribution of genetic diversity of marine species and lead to speciation is challenging because of the lack of obvious dispersal barriers and the likelihood of substantial (passive) dispersal in oceans. In addition, cryptic radiations with sympatric distributions abound in marine species, challenging the allopatric speciation mechanism. Here, we present a phylogeographical study of the marine nematode species complex Rhabditis (Pellioditis) marina to investigate processes shaping genetic structure and speciation. Rhabditis (P.) marina lives on decaying macroalgae in the intertidal, and may therefore disperse over considerable distances. Rhabditis (P.) marina consists of several cryptic species sympatrically distributed at a local scale. Genetic variation in the COI gene was screened in 1362 specimens from 45 locations around the world. Two nuclear DNA genes (ITS and D2D3) were sequenced to infer phylogenetic species. We found evidence for ten sympatrically distributed cryptic species, seven of which show a strong genetic structuring. A historical signature showed evidence for restricted gene flow with occasional long-distance dispersal and range expansions pre-dating the last glacial maximum. Our data also point to a genetic break around the British Isles and a contact zone in the Southern Bight of the North Sea. We provide evidence for the transoceanic distribution of at least one cryptic species (PmIII) and discuss the dispersal capacity of marine nematodes. The allopatric distribution of some intraspecific phylogroups and of closely related cryptic species points to the potential for allopatric speciation in R. (P.) marina.

Molecular evaluation of nominal species in the Ceratitis fasciventris, C. anonae, C. rosa complex (Diptera: Tephritidae).

Ceratitis fasciventris, C. anonae and C. rosa form a complex of economically important fruit fly pests infesting a variety of crops in African countries. Hitherto only adult males of these species can be distinguished easily by morphological characters. Other stages cannot, and for some taxa the taxonomic interpretation and species boundaries remain unclear. In order to clarify phylogenetic relationships and taxonomic status of these species, sequences of mitochondrial (16S, COI, ND6) and nuclear markers (period, ITS1) were analysed in specimens of the three morphospecies throughout the distribution of the complex. Maximum likelihood trees did not recover monophyletic groups corresponding to the morphospecies. Conversely, ND6 and COI divided West African C. fasciventris specimens in two consistent and bootstrap supported clades, involving specimens from Benin and from Mali/Ivory Coast, while the nuclear gene fragments per and ITS1 recovered a well-supported clade corresponding to C. fasciventris from Kenya/Uganda. Hence, the phylogenetic relationships and taxonomic interpretation of the complex appear more intricate than previously hypothesised. The current molecular data do not allow to identify C. fasciventris, C. anonae and C. rosa as distinct phylogenetic species but rather suggest that the morphospecies C. fasciventris is itself a complex of cryptic taxa.