The effect of geographical scale of sampling on DNA barcoding

Delimiting species using single-locus data and the Generalized Mixed Yule Coalescent approach: a revised method and evaluation on simulated data sets

DNA barcoding-type studies assemble single-locus data from large samples of individuals and species, and have provided new kinds of data for evolutionary surveys of diversity. An important goal of many such studies is to delimit evolutionarily significant species units, especially in biodiversity surveys from environmental DNA samples. The Generalized Mixed Yule Coalescent (GMYC) method is a likelihood method for delimiting species by fitting within- and between-species branching models to reconstructed gene trees. Although the method has been widely used, it has not previously been described in detail or evaluated fully against simulations of alternative scenarios of true patterns of population variation and divergence between species. Here, we present important reformulations to the GMYC method as originally specified, and demonstrate its robustness to a range of departures from its simplifying assumptions. The main factor affecting the accuracy of delimitation is the mean population size of species relative to divergence times between them. Other departures from the model assumptions, such as varying population sizes among species, alternative scenarios for speciation and extinction, and population growth or subdivision within species, have relatively smaller effects. Our simulations demonstrate that support measures derived from the likelihood function provide a robust indication of when the model performs well and when it leads to inaccurate delimitations. Finally, the so-called single-threshold version of the method outperforms the multiple-threshold version of the method on simulated data: we argue that this might represent a fundamental limit due to the nature of evidence used to delimit species in this approach. Together with other studies comparing its performance relative to other methods, our findings support the robustness of GMYC as a tool for delimiting species when only single-locus information is available.

Phylogenetic relationships and species delimitation in pinus section trifoliae inferrred from plastid DNA

Recent diversification followed by secondary contact and hybridization may explain complex patterns of intra- and interspecific morphological and genetic variation in the North American hard pines (Pinus section Trifoliae), a group of approximately 49 tree species distributed in North and Central America and the Caribbean islands. We concatenated five plastid DNA markers for an average of 3.9 individuals per putative species and assessed the suitability of the five regions as DNA bar codes for species identification, species delimitation, and phylogenetic reconstruction. The ycf1 gene accounted for the greatest proportion of the alignment (46.9%), the greatest proportion of variable sites (74.9%), and the most unique sequences (75 haplotypes). Phylogenetic analysis recovered clades corresponding to subsections Australes, Contortae, and Ponderosae. Sequences for 23 of the 49 species were monophyletic and sequences for another 9 species were paraphyletic. Morphologically similar species within subsections usually grouped together, but there were exceptions consistent with incomplete lineage sorting or introgression. Bayesian relaxed molecular clock analyses indicated that all three subsections diversified relatively recently during the Miocene. The general mixed Yule-coalescent method gave a mixed model estimate of only 22 or 23 evolutionary entities for the plastid sequences, which corresponds to less than half the 49 species recognized based on morphological species assignments. Including more unique haplotypes per species may result in higher estimates, but low mutation rates, recent diversification, and large effective population sizes may limit the effectiveness of this method to detect evolutionary entities.

Delineating species boundaries using an iterative taxonomic approach: the case of soldierless termites (Isoptera, Termitidae, Apicotermitinae)

Species boundaries are traditionally inferred using morphological characters, although morphology sometimes fails to correctly delineate species. To overcome this limitation, researchers have widely taken advantage of alternative methods such as DNA barcoding or analysis of cuticular hydrocarbons (CHs) profiles, but rarely use them simultaneously in an iterative taxonomic approach. Here, we follow such an approach using morphology, DNA barcoding and CHs profiles to precisely discriminate species of soldierless termites, a diversified clade constituting about one-third of the Neotropical termite species richness, but poorly resolved taxonomically due to the paucity of useful characters. We sampled soldierless termites in various forest types of the Nouragues Nature Reserve, French Guiana. Our results show that morphological species determination generally matches DNA barcoding, which only suggests the existence of three cryptic species in the 31 morphological species. Among them, Longustitermes manni is the only species whose splitting is corroborated by ecological data, other widely distributed species being supported by DNA barcoding. On the contrary, although CHs profiles provide a certain taxonomic signal, they often suggest inconsistent groupings which are not supported by other methods. Overall, our data support DNA barcoding and morphology as two efficient methods to distinguish soldierless termite species.

Integrative species delimitation in photosynthetic sea slugs reveals twenty candidate species in three nominal taxa studied for drug discovery, plastid symbiosis or biological control

DNA barcoding can highlight taxa in which conventional taxonomy underestimates species richness, identifying mitochondrial lineages that may correspond to unrecognized species. However, key assumptions of barcoding remain untested for many groups of soft-bodied marine invertebrates with poorly resolved taxonomy. Here, we applied an integrative approach for species delimitation to herbivorous sea slugs in clade Sacoglossa, in which unrecognized diversity may complicate studies of drug discovery, plastid endosymbiosis, and biological control. Using the mitochondrial barcoding COI gene and the nuclear histone 3 gene, we tested the hypothesis that three widely distributed "species" each comprised a complex of independently evolving lineages. Morphological and reproductive characters were then used to evaluate whether each lineage was distinguishable as a candidate species. The "circumtropical" Elysia ornata comprised a Caribbean species and four Indo-Pacific candidate species that are potential sources of kahalalides, anti-cancer compounds. The "monotypic" and highly photosynthetic Plakobranchus ocellatus, used for over 60 years to study chloroplast symbiosis, comprised 10 candidate species. Finally, six candidate species were distinguished in the Elysia tomentosa complex, including potential biological control agents for invasive green algae (Caulerpa spp.). We show that a candidate species approach developed for vertebrates effectively categorizes cryptic diversity in marine invertebrates, and that integrating threshold COI distances with non-molecular character data can delimit species even when common assumptions of DNA barcoding are violated.

DNA barcode analysis of butterfly species from Pakistan points towards regional endemism

DNA barcodes were obtained for 81 butterfly species belonging to 52 genera from sites in north-central Pakistan to test the utility of barcoding for their identification and to gain a better understanding of regional barcode variation. These species represent 25% of the butterfly fauna of Pakistan and belong to five families, although the Nymphalidae were dominant, comprising 38% of the total specimens. Barcode analysis showed that maximum conspecific divergence was 1.6%, while there was 1.7–14.3% divergence from the nearest neighbour species. Barcode records for 55 species showed <2% sequence divergence to records in the Barcode of Life Data Systems (BOLD), but only 26 of these cases involved specimens from neighbouring India and Central Asia. Analysis revealed that most species showed little incremental sequence variation when specimens from other regions were considered, but a threefold increase was noted in a few cases. There was a clear gap between maximum intraspecific and minimum nearest neighbour distance for all 81 species. Neighbour-joining cluster analysis showed that members of each species formed a monophyletic cluster with strong bootstrap support. The barcode results revealed two provisional species that could not be clearly linked to known taxa, while 24 other species gained their first coverage. Future work should extend the barcode reference library to include all butterfly species from Pakistan as well as neighbouring countries to gain a better understanding of regional variation in barcode sequences in this topographically and climatically complex region.

Integrative analysis of DNA phylogeography and morphology of the European rose chafer (Cetonia aurata) to infer species taxonomy and patterns of postglacial colonisation in Europe

Integrative taxonomy has been proposed as a framework to unify new conceptual and methodological developments in quantitative assessment of trait variation used in species delimitation, but empirical studies in this young branch of systematics are rare. Here we use standard phylogenetic and parsimony network analyses on nuclear and mitochondrial DNA (Cox1, ITS1) of 230 individuals from 65 European sampling sites in order to deduce population structure of Cetonia beetles from geno- and haplotypes. Statistical measures of population differentiation are inferred on genealogical and geographical scales to test hypotheses about species limits and population history. By combining results of phylogenetic structure with features of morphology, including genital shape morphometrics and discrete external body characters, as well as with measures of population genetics, we attempt to integrate the results as a test of the validity of species limits, in particular of currently recognised subspecies. Despite high Cox1 divergence between some haplotype lineages, even some sympatric lineages (9%, e.g. N2 vs. N4), nDNA and morphology, as well as pattern of geographical and genealogical divergence measured by AMOVA analysis did not support the hypothesis of separate species. Highest divergence in nuclear markers was found among Italian populations of C. aurata pisana and C. a. sicula, and moderately high fixation indices along measurable morphological divergence suggest the correctness of their status as 'subspecies'. Divergence time estimates of the lineages suggest a glacial divergence in different refugia between the major haplogroups, while population differentiation in mtDNA among these was primarily attributable to restricted gene flow caused by geographic isolation.

Bayesian species delimitation reveals generalist and specialist parasitic wasps on Galerucella beetles (Chrysomelidae): sorting by herbivore or plant host

BACKGROUND

To understand the ecological and evolutionary consequences of species interactions in food webs necessitates that interactions are properly identified. Genetic analyses suggest that many supposedly generalist parasitoid species should rather be defined as multiple species with a more narrow diet, reducing the probability that such species may mediate indirect interactions such as apparent competition among hosts. Recent studies showed that the parasitoid Asecodes lucens mediate apparent competition between two hosts, Galerucella tenella and G. calmariensis, affecting both interaction strengths and evolutionary feedbacks. The same parasitoid was also recorded from other species in the genus Galerucella, suggesting that similar indirect effects may also occur for other species pairs.

METHODS

To explore the possibility of such interactions, we sequenced mitochondrial and nuclear genetic markers to resolve the phylogeny of both host and parasitoid and to test the number of parasitoid species involved. We thus collected 139 Galerucella larvae from 8 host plant species and sequenced 31 adult beetle and 108 parasitoid individuals.

RESULTS

The analysis of the Galerucella data, that also included sequences from previous studies, verified the five species previously documented as reciprocally monophyletic, but the Bayesian species delimitation for A. lucens suggested 3-4 cryptic taxa with a more specialised host use than previously suggested. The gene data analyzed under the multispecies coalescent model allowed us to reconstruct the species tree phylogeny for both host and parasitoid and we found a fully congruent coevolutionary pattern suggesting that parasitoid speciation followed upon host speciation.

CONCLUSION

Using multilocus sequence data in a Bayesian species delimitation analysis we propose that hymenopteran parasitoids of the genus Asecodes that infest Galerucella larvae constitute at least three species with narrow diet breath. The evolution of parasitoid Asecodes and host Galerucella show a fully congruent coevolutionary pattern. This finding strengthens the hypothesis that the parasitoid in host search uses cues of the host rather than more general cues of both host and plant.

A multi-locus approach to barcoding in the Anopheles strodei subgroup (Diptera: Culicidae)

Background

The ability to successfully identify and incriminate pathogen vectors is fundamental to effective pathogen control and management. This task is confounded by the existence of cryptic species complexes. Molecular markers can offer a highly effective means of species identification in such complexes and are routinely employed in the study of medical entomology. Here we evaluate a multi-locus system for the identification of potential malaria vectors in the Anopheles strodei subgroup.

Methods

Larvae, pupae and adult mosquitoes (n = 61) from the An. strodei subgroup were collected from 21 localities in nine Brazilian states and sequenced for the COI, ITS2 and white gene. A Bayesian phylogenetic approach was used to describe the relationships in the Strodei Subgroup and the utility of COI and ITS2 barcodes was assessed using the neighbor joining tree and “best close match” approaches.

Results

Bayesian phylogenetic analysis of the COI, ITS2 and white gene found support for seven clades in the An. strodei subgroup. The COI and ITS2 barcodes were individually unsuccessful at resolving and identifying some species in the Subgroup. The COI barcode failed to resolve An. albertoi and An. strodei but successfully identified approximately 92% of all species queries, while the ITS2 barcode failed to resolve An. arthuri and successfully identified approximately 60% of all species queries. A multi-locus COI-ITS2 barcode, however, resolved all species in a neighbor joining tree and successfully identified all species queries using the “best close match” approach.

Conclusions

Our study corroborates the existence of An. albertoi, An. CP Form and An. strodei in the An. strodei subgroup and identifies four species under An. arthuri informally named A-D herein. The use of a multi-locus barcode is proposed for species identification, which has potentially important utility for vector incrimination. Individuals previously found naturally infected with Plasmodium vivax in the southern Amazon basin and reported as An. strodei are likely to have been from An. arthuri C identified in this study.

The evolution of host associations in the parasitic wasp genus Ichneumon (Hymenoptera: Ichneumonidae): convergent adaptations to host pupation sites

Background

The diversification of organisms with a parasitic lifestyle is often tightly linked to the evolution of their host associations. If a tight host association exists, closely related species tend to attack closely related hosts; host associations are less stable if associations are determined by more plastic traits like parasitoid searching and oviposition behaviour. The pupal-parasitoids of the genus Ichneumon attack a variety of macrolepidopteran hosts. They are either monophagous or polyphagous, and therefore offer a promissing system to investigate the evolution of host associations. Ichneumon was previously divided into two groups based on general body shape; however, a stout shape has been suggested as an adaptation to buried host pupation sites, and might thus not represent a reliable phylogenetic character.

Results

We here reconstruct the first molecular phylogeny of the genus Ichneumon using two mitochondrial (CO1 and NADH1) and one nuclear marker (28S). The resulting phylogeny only supports monophyly of Ichneumon when Ichneumon lugens Gravenhorst, 1829 (formerly in Chasmias, stat. rev.) and Ichneumon deliratorius Linnaeus, 1758 (formerly Coelichneumon) are included. Neither parasitoid species that attack hosts belonging to one family nor those attacking butterflies (Rhopalocera) form monophyletic clades. Ancestral state reconstructions suggest multiple transitions between searching for hosts above versus below ground and between a stout versus elongated body shape. A model assuming correlated evolution between the two characters was preferred over independent evolution of host-searching niche and body shape.

Conclusions

Host relations, both in terms of phylogeny and ecology, evolved at a high pace in the genus Ichneumon. Numerous switches between hosts of different lepidopteran families have occurred, a pattern that seems to be the rule among idiobiont parasitoids. A stout body and antennal shape in the parasitoid female is confirmed as an ecological adaptation to host pupation sites below ground and has evolved convergently several times. Morphological characters that might be involved in adaptation to hosts should be avoided as diagnostic characters for phylogeny and classification, as they can be expected to show high levels of homoplasy.

Reexamination of the species assignment of Diacavolinia pteropods using DNA barcoding

Thecosome pteropods (Mollusca, Gastropoda) are an ecologically important, diverse, and ubiquitous group of holoplanktonic animals that are the focus of intense research interest due to their external aragonite shell and vulnerability to ocean acidification. Characterizing the response of these animals to low pH and other environmental stressors has been hampered by continued uncertainty in their taxonomic identification. An example of this confusion in species assignment is found in the genus Diacavolinia. All members of this genus were originally indentified as a single species, Cavolinia longirostris, but over the past fifty years the taxonomy has been revisited multiple times; currently the genus comprises 22 different species. This study examines five species of Diacavolinia, including four sampled in the Northeast Atlantic (78 individuals) and one from the Eastern tropical North Pacific (15 individuals). Diacavolina were identified to species based on morphological characteristics according to the current taxonomy, photographed, and then used to determine the sequence of the "DNA barcoding" region of the cytochrome c oxidase subunit I (COI). Specimens from the Atlantic, despite distinct differences in shell morphology, showed polyphyly and a genetic divergence of <3% (K2P distance) whereas the Pacific and Atlantic samples were more distant (≈ 19%). Comparisons of Diacavolinia spp. with other Cavolinia spp. reveal larger distances (≈ 24%). These results indicate that specimens from the Atlantic comprise a single monophyletic species and suggest possible species-level divergence between Atlantic and Pacific populations. The findings support the maintenance of Diacavolinia as a separate genus, yet emphasize the inadequacy of our current taxonomic understanding of pteropods. They highlight the need for accurate species identifications to support estimates of biodiversity, range extent and natural exposure of these planktonic calcifiers to environmental variability; furthermore, the apparent variation of the pteropods shell may have implications for our understanding of the species' sensitivity to ocean acidification.

Improving environmental and biodiversity monitoring in the Baltic Sea using DNA barcoding of Chironomidae (Diptera)

As for many other regions, environmental and biodiversity monitoring of the brackish Baltic Sea suffers from low species resolution for several taxa. One such case is the benthic larvae of midges Chironomidae (Diptera), which are estimated to constitute about 30% of the macrozoobenthos species of the Baltic Sea and are important indicators of environmental quality. We assessed the usefulness of COI (cytochrome oxidase I) gene barcoding to improve species resolution and its potential for implementation in monitoring programmes. Neighbour-Joining, Maximum parsimony and Bayesian-inference analyses all provided high congruency with morphological analyses of adult males for almost all 42 species studied. Barcoding was helpful to elucidate some cases of taxonomical difficulties, such as synonyms. In contrast to the high identification accuracy when using our local database, there were a number of cases where matching with GenBank and BOLD provided puzzling results. For reliable species identification at least 15-30 specimens from 5-10 well-distributed sites within the geographical range of the species might be needed in a database to adequately cover the intraspecific variability of chironomids. Implementation of DNA barcoding, as applied here, in monitoring would result in an increase from at present less than 10% to more than 90% successful chironomid species identification of Baltic Sea benthic samples, as it also would for many nearby lakes. Routine monitoring of benthic environmental samples based on Next-Generation sequencing techniques would provide a cost effective way to obtain a taxonomically much more complete assessment of environmental quality and biodiversity, as required by EU directives and national legislation.

Endemism and diversification in freshwater insects of Madagascar revealed by coalescent and phylogenetic analysis of museum and field collections

The biodiversity and endemism of Madagascar are among the most extraordinary and endangered in the world. This includes the island's freshwater biodiversity, although detailed knowledge of the diversity, endemism, and biogeographic origin of freshwater invertebrates is lacking. The aquatic immature stages of mayflies (Ephemeroptera) are widely used as bio-indicators and form an important component of Malagasy freshwater biodiversity. Many species are thought to be microendemics, restricted to single river basins in forested areas, making them particularly sensitive to habitat reduction and degradation. The Heptageniidae are a globally diverse family of mayflies (>500 species) but remain practically unknown in Madagascar except for two species described in 1996. The standard approach to understanding their diversity, endemism, and origin would require extensive field sampling on several continents and years of taxonomic work followed by phylogenetic analysis. Here we circumvent this using museum collections and freshly collected individuals in a combined approach of DNA taxonomy and phylogeny. The coalescent-based GMYC analysis of DNA barcode data (mitochondrial COI) revealed 14 putative species on Madagascar, 70% of which were microendemics. A phylogenetic analysis that included African and Asian species and data from two mitochondrial and four nuclear loci indicated the Malagasy Heptageniidae are monophyletic and sister to African species. The genus Compsoneuria is shown to be paraphyletic and the genus Notonurus is reinstalled for African and Malagasy species previously placed in Compsoneuria. A molecular clock excluded a Gondwanan vicariance origin and instead favoured a more recent overseas colonization of Madagascar. The observed monophyly and high microendemism highlight their conservation importance and suggest the DNA-based approach can rapidly provide information on the diversity, endemism, and origin of freshwater biodiversity. Our results underline the important role that museum collections can play in molecular studies, especially in critically endangered biodiversity hotspots like Madagascar where entire species or populations may go extinct very quickly.

DNA barcoding of Northern Nearctic Muscidae (Diptera) reveals high correspondence between morphological and molecular species limits

BACKGROUND

Various methods have been proposed to assign unknown specimens to known species using their DNA barcodes, while others have focused on using genetic divergence thresholds to estimate "species" diversity for a taxon, without a well-developed taxonomy and/or an extensive reference library of DNA barcodes. The major goals of the present work were to: a) conduct the largest species-level barcoding study of the Muscidae to date and characterize the range of genetic divergence values in the northern Nearctic fauna; b) evaluate the correspondence between morphospecies and barcode groupings defined using both clustering-based and threshold-based approaches; and c) use the reference library produced to address taxonomic issues.

RESULTS

Our data set included 1114 individuals and their COI sequences (951 from Churchill, Manitoba), representing 160 morphologically-determined species from 25 genera, covering 89% of the known fauna of Churchill and 23% of the Nearctic fauna. Following an iterative process through which all specimens belonging to taxa with anomalous divergence values and/or monophyly issues were re-examined, identity was modified for 9 taxa, including the reinstatement of Phaonia luteva (Walker) stat. nov. as a species distinct from Phaonia errans (Meigen). In the post-reassessment data set, no distinct gap was found between maximum pairwise intraspecific distances (range 0.00-3.01%) and minimum interspecific distances (range: 0.77-11.33%). Nevertheless, using a clustering-based approach, all individuals within 98% of species grouped with their conspecifics with high (>95%) bootstrap support; in contrast, a maximum species discrimination rate of 90% was obtained at the optimal threshold of 1.2%. DNA barcoding enabled the determination of females from 5 ambiguous species pairs and confirmed that 16 morphospecies were genetically distinct from named taxa. There were morphological differences among all distinct genetic clusters; thus, no cases of cryptic species were detected.

CONCLUSIONS

Our findings reveal the great utility of building a well-populated, species-level reference barcode database against which to compare unknowns. When such a library is unavailable, it is still possible to obtain a fairly accurate (within ~10%) rapid assessment of species richness based upon a barcode divergence threshold alone, but this approach is most accurate when the threshold is tuned to a particular taxon.

Allopatry as a gordian knot for taxonomists: patterns of DNA barcode divergence in arctic-alpine lepidoptera

Many cold adapted species occur in both montane settings and in the subarctic. Their disjunct distributions create taxonomic complexity because there is no standardized method to establish whether their allopatric populations represent single or different species. This study employs DNA barcoding to gain new perspectives on the levels and patterns of sequence divergence among populations of 122 arctic-alpine species of Lepidoptera from the Alps, Fennoscandia and North America. It reveals intraspecific variability in the barcode region ranging from 0.00-10.08%. Eleven supposedly different species pairs or groups show close genetic similarity, suggesting possible synonymy in many cases. However, a total of 33 species show evidence of cryptic diversity as evidenced by the presence of lineages with over 2% maximum barcode divergence in Europe, in North America or between the two continents. Our study also reveals cases where taxonomic names have been used inconsistently between regions and exposes misidentifications. Overall, DNA barcodes have great potential to both increase taxonomic resolution and to make decisions concerning the taxonomic status of allopatric populations more objective.