Generalized structure and evolution of ITS1 and ITS2 rDNA in black flies (Diptera: Simuliidae)

Integrated taxonomy of black flies (Diptera: Simuliidae) reveals unexpected diversity in the most arid ecosystem of Europe

The family Simuliidae includes more than 2000 species of black flies worldwide. Their morphological uniformity creates difficulty for species identification, which limits our knowledge of their ecology and vectorial role. We investigated the systematics of black flies in a semi-arid area of the Iberian Peninsula, an ecologically harsh environment for these organisms. Sampling adult black flies in three different habitats (by means of CDC traps) and in avian nest boxes and collecting immature stages in high-salinity rills provided a representative sample of the component species. A combination of approaches, including morphological, chromosomal, and molecular (based on the mitochondrial cytochrome C oxidase subunit I (COI) and internal transcribed spacer 2 (ITS2) genes) revealed five species: four common species (Simulium intermedium, S. petricolum, S. pseudequinum, and S. rubzovianum) and the first European record for S. mellah. Barcoding gap and phylogenetic analyses revealed that ITS2 is a key marker to identify the species, whereas the COI marker does not provide enough resolution to identify some species or infer their phylogenetic relationships. Morphological and chromosomal features are also provided to identify S. mellah unequivocally. Our study highlights the need for integrated studies of black flies in ecologically extreme habitats to increase our knowledge of their distribution, ecology, and potential risks for public health.

Genetic diversity and geographic distribution of parasitic flies of the Philornis torquans complex in Argentina

Philornis Meinert 1890 (Diptera: Muscidae) is a genus of flies that parasitize birds in the Neotropical region. The characteristics of the host-parasite interactions and its consequences may depend on the Philornis species involved, and thus precise identification of these parasites is crucial for the interpretation of ecological and epidemiological studies. However, morphological identification of Argentine Philornis species is elusive while molecular evidence points towards the existence of a complex of cryptic species or lineages undergoing a speciation process, which were named the 'Philornis torquans complex'. Herein the authors extended the current knowledge on the systematics and biogeography of parasitic Philornis flies from Argentina, analysing samples collected in several ecoregions, including the Atlantic Forest, Iberá Wetlands, Open Fields and Grasslands, Espinal, Pampa, Dry Chaco, Humid Chaco, Delta and Paraná River Islands, Monte of Plains and Plateaus. The results of the present study strengthen the evidence on previously described Philornis genotypes using four genetic markers (ITS2, COI, ND6, 12S rRNA). The authors report new patterns of occurrence and describe the presence of a novel genotype of subcutaneous Philornis. In addition, the present study unveils ecological niche differences among genotypes of the Philornis torquans complex in southern South America.

Characterization of the complete mitochondrial genome of Simulium (Byssodon) maculatum (Diptera: Simuliidae) and its phylogenetic implications

The mitochondrial (mt) genome of the black fly Simulium (Byssodon) maculatum, a pest of great importance for both humans and livestock, is sequenced and annotated for the first time. The genome structure, gene order and codon usage are typical among Diptera mt genomes. The mt genome is circular and 15,799 bp in length with 13 protein coding genes (PCGs), 22 transfer RNA genes (tRNAs), two ribosomal RNA genes (rRNAs) and a control region (CR), and with weakly positive AT-skew (0.02) and negative GC-skew (-0.12). Phylogenetic relationships of 16 species representing five families of Culicomorpha and two outgroups, based on mt genome data, were analyzed using both Maximum Likelihood and Bayesian methods. The monophyly of Culicomorpha is well supported, while Chironomoidea is indicated as a paraphyletic group. The well supported monophyletic Simuliidae is the sister group to Culicidae.

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.

Technical Note: "Mitochondrial and nuclear DNA approaches for reliable identification of Lucilia (Diptera, Calliphoridae) species of forensic interest from Southern Europe"

In forensic entomology, rapid and unambiguous identification of blowfly species is a critical prerequisite for accurately estimating the post-mortem interval (PMI). The conventional diagnosis of cadaveric entomofauna based on external characters is hampered by the morphological similarities between species, especially in immature stages. Genetic analysis has been shown to allow precise and reliable diagnosis and delimitation of insect species. Nevertheless, the taxonomy of some species remains unresolved. This study was focused on improving the effectiveness and accuracy of analysis based on the widely used cytochrome c oxidase subunit I barcode region (COI barcode, 658 bp), complemented by other mitochondrial and nuclear regions, such as cytochrome b (Cyt-b, 307 bp) and the second internal transcribed spacer (ITS2, 310-331 bp), for the identification of Southern European blowflies. We analyzed a total of 209 specimens, collected from 38 human corpses, belonging to three Calliphoridae genera and seven species: Chrysomya (Ch. albiceps), Calliphora (C. vicina and C. vomitoria), and Lucilia (L. sericata, L. ampullacea, L. caesar and L. illustris). These species are the most common PMI indicators in Portugal. The results revealed that unambiguous separation of species of the Lucilia genus requires different loci from the barcode region. Furthermore, we conclude that the ITS2 (310-331 bp) molecular marker is a promising diagnostic tool because its inter-specific discriminatory power enables unequivocal and consistent distinctions to be made, even between closely related species (L. caesar-L. illustris). This work also contributes new genetic data that may be of interest in performing species diagnosis for Southern European blowflies. Notably, to the best of our knowledge, we provide the first records of the Cyt-b (307 bp) locus for L. illustris and the ITS2 (310-331 bp) region for Iberian Peninsula Lucilia species.

New tools and insights to assist with the molecular identification of Simulium guianense s.l., main Onchocerca volvulus vector within the highland areas of the Amazonia onchocerciasis focus

Following the success of the Onchocerciasis Elimination Programme for the Americas (OEPA), there is now just one Latin American onchocerciasis focus where onchocerciasis transmission is described as 'on-going:' the Amazonia Onchocerciasis focus. In the hyperendemic highland areas of the Amazonia focus, Simulium guianense s.l. Wise are the most important vectors of the disease. Populations of S. guianense s.l. are, however, known to vary in their cytogenetics and in a range of behaviours, including in their biting habits. In the hypoendemic lowland areas of the Amazonia focus, for example, S. guianense s.l. are generally regarded as zoophilic and consequently unimportant to disease transmission. Robust tools, to discriminate among various populations of S. guianense s.l. have, however, not yet been developed. In the work reported here, we have assessed the utility of a ribosomal DNA sequence fragment spanning the nuclear ribosomal ITS-1, ITS-2 and 5.8S sequence regions and a ∼850 nucleotide portion of the mitochondrial cytochrome oxidase gene (CO1) for species-level identification and for resolving the within species substructuring. We report here how we have generated 78 CO1 sequences from a rich set of both zoophilic and anthropophilic populations of S. guianense s.l. that were collected from eight sites that are broadly distributed across Brazil. Consistent with previous findings, our analysis supports the genetic isolation of Simulium litobranchium from S. guianense s.l. In contrast with previous findings, however, our results did not provide support for the divergence of the two species prior to the radiation of S. guianense s.l. In our analysis of the S. guianense s.l. ribosomal DNA sequence trace files we generated, we provide clear evidence of multiple within-specimen single nucleotide polymorphisms and indels suggesting that S. guianense s.l. ribosomal DNA is not a good target for conventional DNA barcoding. This is the first report of S. guianense s.l. within individual ribosomal DNA variation and thus the first evidence that the species is not subject to the normal effects of concerted evolution. Collectively, these data illustrate the need for diverse sampling in the development of robust molecular tools for vector identification and suggest that ribosomal DNA might be able to assist with resolving S. guianense s.l. species substructuring that C01 barcoding has hitherto failed to.

A close-up view on ITS2 evolution and speciation - a case study in the Ulvophyceae (Chlorophyta, Viridiplantae)

BACKGROUND

The second Internal Transcriber Spacer (ITS2) is a fast evolving part of the nuclear-encoded rRNA operon located between the 5.8S and 28S rRNA genes. Based on crossing experiments it has been proposed that even a single Compensatory Base Change (CBC) in helices 2 and 3 of the ITS2 indicates sexual incompatibility and thus separates biological species. Taxa without any CBC in these ITS2 regions were designated as a 'CBC clade'. However, in depth comparative analyses of ITS2 secondary structures, ITS2 phylogeny, the origin of CBCs, and their relationship to biological species have rarely been performed. To gain 'close-up' insights into ITS2 evolution, (1) 86 sequences of ITS2 including secondary structures have been investigated in the green algal order Ulvales (Chlorophyta, Viridiplantae), (2) after recording all existing substitutions, CBCs and hemi-CBCs (hCBCs) were mapped upon the ITS2 phylogeny, rather than merely comparing ITS2 characters among pairs of taxa, and (3) the relation between CBCs, hCBCs, CBC clades, and the taxonomic level of organisms was investigated in detail.

RESULTS

High sequence and length conservation allowed the generation of an ITS2 consensus secondary structure, and introduction of a novel numbering system of ITS2 nucleotides and base pairs. Alignments and analyses were based on this structural information, leading to the following results: (1) in the Ulvales, the presence of a CBC is not linked to any particular taxonomic level, (2) most CBC 'clades' sensu Coleman are paraphyletic, and should rather be termed CBC grades. (3) the phenetic approach of pairwise comparison of sequences can be misleading, and thus, CBCs/hCBCs must be investigated in their evolutionary context, including homoplasy events (4) CBCs and hCBCs in ITS2 helices evolved independently, and we found no evidence for a CBC that originated via a two-fold hCBC substitution.

CONCLUSIONS

Our case study revealed several discrepancies between ITS2 evolution in the Ulvales and generally accepted assumptions underlying ITS2 evolution as e.g. the CBC clade concept. Therefore, we developed a suite of methods providing a critical 'close-up' view into ITS2 evolution by directly tracing the evolutionary history of individual positions, and we caution against a non-critical use of the ITS2 CBC clade concept for species delimitation.

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.