Sequence-based species delimitation for the DNA taxonomy of undescribed insects

Cataloging the very large number of undescribed species of insects could be greatly accelerated by automated DNA based approaches, but procedures for large-scale species discovery from sequence data are currently lacking. Here, we use mitochondrial DNA variation to delimit species in a poorly known beetle radiation in the genus Rivacindela from arid Australia. Among 468 individuals sampled from 65 sites and multiple morphologically distinguishable types, sequence variation in three mtDNA genes (cytochrome oxidase subunit 1, cytochrome b, 16S ribosomal RNA) was strongly partitioned between 46 or 47 putative species identified with quantitative methods of species recognition based on fixed unique ("diagnostic") characters. The boundaries between groups were also recognizable from a striking increase in branching rate in clock-constrained calibrated trees. Models of stochastic lineage growth (Yule models) were combined with coalescence theory to develop a new likelihood method that determines the point of transition from species-level (speciation and extinction) to population-level (coalescence) evolutionary processes. Fitting the location of the switches from speciation to coalescent nodes on the ultrametric tree of Rivacindela produced a transition in branching rate occurring at 0.43 Mya, leading to an estimate of 48 putative species (confidence interval for the threshold ranging from 47 to 51 clusters within 2 logL units). Entities delimited in this way exhibited biological properties of traditionally defined species, showing coherence of geographic ranges, broad congruence with morphologically recognized species, and levels of sequence divergence typical for closely related species of insects. The finding of discontinuous evolutionary groupings that are readily apparent in patterns of sequence variation permits largely automated species delineation from DNA surveys of local communities as a scaffold for taxonomy in this poorly known insect group.     

DNA-based taxonomy for associating adults and larvae in multi-species assemblages of chafers (Coleoptera: Scarabaeidae)

DNA sequences provide a universal character system in taxonomy for associating all developmental stages of organisms, but ambiguity arises due to genetic variation within species. The problem is compounded where target groups are less well studied or incompletely represented in DNA databases. Here we investigate the utility of DNA for larval-adult species associations within chafer (Coleoptera: Scarabaeidae) communities from four sites in the tropical lowlands of Nepal. We sequenced ca. 1600 bp of mitochondrial cox1 and rrnL and 700 bp of nuclear 28S rRNA from 250 larval and adult specimens. Individuals were grouped into putative species using statistical parsimony analysis and population aggregation analysis (PAA), whereby specimens from each locality were grouped according to the presence of diagnostic nucleotides. In addition, species membership was determined based on shifts in branching rates on clock-constrained trees to detect the putative transition from speciation to population coalescence patterns. Using these two methods we delineated between 48 and 56 groups, of which 16-20 were composed of larval and adult individuals. Nuclear and mtDNA-based groups were highly congruent; variation of 28S rRNA within groups was very low, while one widespread 28S rRNA genotype was universally found in a paraphyletic group of five mtDNA clusters. Linnean names could be assigned to 19 groups, and hence between 86.1% and 92.7% of larval specimens could be associated to species by their membership in clearly delineated groups that contained fully identified adults. The remaining larvae were delineated as five species, four of which could be assigned to Anomala or Adoretus based on their phylogenetic position. We conclude that the sequence variation was highly structured in this complex assemblage of chafers and that any given individual (larva or adult) can be readily associated with a particular DNA group using the criterion of diagnos ability. The association of different developmental stages therefore becomes a matter of determining the extent of the DNA-based groups, rather than matching of sequences from adult and larval individuals. This indicates the need for a purely sequence-based taxonomic system when associating different life stages via DNA.     

Speciation and DNA barcodes: testing the effects of dispersal on the formation of discrete sequence clusters

Large-scale sequencing of short mtDNA fragments for biodiversity inventories ('DNA barcoding') indicates that sequence variation in animal mtDNA is highly structured and partitioned into discrete genetic clusters that correspond broadly to species-level entities. Here we explore how the migration rate, an important demographic parameter that is directly related to population isolation, might affect variation in the strength of mtDNA clustering among taxa. Patterns of mtDNA variation were investigated in two groups of beetles that both contain lineages occupying habitats predicted to select for different dispersal abilities: predacious diving beetles (Dytiscidae) in the genus Bidessus from lotic and lentic habitats across Europe and darkling beetles (Tenebrionidae) in the genus Eutagenia from sand and other soil types in the Aegean Islands. The degree of genetic clustering was determined using the recently developed 'mixed Yule coalescent' (MYC) model that detects the transition from between-species to within-population branching patterns. Lineages from presumed stable habitats, and therefore displaying lower dispersal ability and migration rates, showed greater levels of mtDNA clustering and geographical subdivision than their close relatives inhabiting ephemeral habitats. Simulations of expected patterns of mtDNA variation under island models showed that MYC clusters are only detected when the migration rates are much lower than the value of Nm=1 typically used to define the threshold for neutral genetic divergence. Therefore, discrete mtDNA clusters provide strong evidence for independently evolving populations or species, but their formation is suppressed even under very low levels of dispersal.     

Delineating Species with DNA Barcodes: A Case of Taxon Dependent Method Performance in Moths

The accelerating loss of biodiversity has created a need for more effective ways to discover species. Novel algorithmic approaches for analyzing sequence data combined with rapidly expanding DNA barcode libraries provide a potential solution. While several analytical methods are available for the delineation of operational taxonomic units (OTUs), few studies have compared their performance. This study compares the performance of one morphology-based and four DNA-based (BIN, parsimony networks, ABGD, GMYC) methods on two groups of gelechioid moths. It examines 92 species of Finnish Gelechiinae and 103 species of Australian Elachistinae which were delineated by traditional taxonomy. The results reveal a striking difference in performance between the two taxa with all four DNA-based methods. OTU counts in the Elachistinae showed a wider range and a relatively low (ca. 65%) OTU match with reference species while OTU counts were more congruent and performance was higher (ca. 90%) in the Gelechiinae. Performance rose when only monophyletic species were compared, but the taxon-dependence remained. None of the DNA-based methods produced a correct match with non-monophyletic species, but singletons were handled well. A simulated test of morphospecies-grouping performed very poorly in revealing taxon diversity in these small, dull-colored moths. Despite the strong performance of analyses based on DNA barcodes, species delineated using single-locus mtDNA data are best viewed as OTUs that require validation by subsequent integrative taxonomic work.     

Morphological and molecular variation in tiger beetles of the Cicindela hybrida complex: is an ‘integrative taxonomy’ possible?

Current taxon assignments at the species level are frequently discordant with DNA-based analyses. Recent studies on tiger beetles in the Cicindela hybrida complex identified discordance between mtDNA patterns and the entities currently defined by the taxonomic literature. To test the accuracy of morphologically delimited groups, five named taxa (species) from 24 representative sampling sites across Europe were scored for 41 external morphological characters. Three of the named taxa were 'diagnosable', that is, defined by between one and three characters unique to each group. Newly sequenced ITS1 and existing mitochondrial cox1 markers established 20 and 22 different haplotypes, respectively, but only cox1 produced (four) diagnosable units. Phylogenetic analysis and statistical parsimony networks showed poor congruence of character variation with the taxonomic entities (and each other). Variation in morphological characters was therefore tested directly for association with DNA-based nesting groups at various hierarchical levels using permutational contingency analysis. Significant statistical associations of 11 (of 13 variable) morphological characters were observed with nesting groups from ITS1 and mitochondrial DNA markers, predominantly at the 4-step level. The analysis demonstrates the need for formal tests of congruence with morphological variation at the level of individual characters, a step that is omitted from recent studies of 'integrative taxonomy'. In addition, statistical correlation of particular morphological characters with DNA-based nesting groups can identify the lowest hierarchical level at which various character sets show congruence, as a means to define evolutionarily separated entities supported by diverse data sources.     

DNA taxonomy, phylogeny and Pleistocene diversification of the Cicindela hybrida species group (Coleoptera: Cicindelidae)

Species delimitation is complicated where morphological variation is continuous or poorly subdivided, but for taxonomic convenience it is common practice to separate and name geographical groups to capture this variation. DNA-based approaches may be used to test if these groups in fact represent historically divided, discrete species entities. The Cicindela hybrida complex (Coleoptera: Cicindelidae) is an assemblage of up to seven morphologically recognized species and 15 subspecies with wide distribution in the Palaearctic region. We sequenced a discontinuous segment of 1899 bp of mtDNA including three regions (coxI, rrnL+trnL2+nad1, cob) for a total of 99 specimens from 36 sampling localities across Europe, revealing 48 haplotypes. Four major clades could be identified corresponding to geographical groups from central Iberia, Ukraine, central Europe, and a band from the Atlantic Iberian coast to northern Europe. Taking into account further subdivisions within these clades, four of the six named species included in the analysis were recognizable by applying various procedures for species delimitation. Age estimates from calibrated molecular clocks date the diversification of the hybrida group within the past 2 million years (Myr), and the separation of the northern clade within 0.4 Myr. Nested clade analysis revealed the rapid range expansion of the northern group consistent with postglacial dispersal, but we did not find support for specific source population(s) in the postulated southern refugia. The evolutionary framework based on mtDNA sequences is shown to identify species entities as discrete clusters of closely related sequences and provides an objective system for delineating and recognizing hierarchically structured groups. In the case of the C. hybrida complex, these groups largely coincided with those established from morphology. The study adds further support to the utility of mtDNA-based sequence profiles (the 'DNA taxonomy') as a rapid and objective synthesis of evolutionary diversity and as reference system for communication.     

Structure and evolution of the mitochondrial control region of the pollen beetle Meligethes thalassophilus (Coleoptera: Nitidulidae).

The organization of the mitochondrial DNA (mtDNA) control region (CR) of the pollen beetle Meligethes thalassophilus is described. This mtDNA CR represents the longest sequenced for beetles so far, since the entire nucleotide sequence ranges from approximately 5000 to approximately 5500 bp. The CR of M. thalassophilus is organized in three distinct domains: a conserved domain near the tRNAIle gene, a variable domain flanking the 12S rRNA gene, and a relatively large central tandem array made up of a variable number of approximately 170 bp repeats that is responsible for the intraspecific length variation observed. Like other CRs found in insects, the M. thalassophilus CR contains two long homopolymeric runs that may be involved in mtDNA replication. Furthermore, conserved stem-and-loop structures in the repetitive domain were identified and their possible role in generating length variation is examined. Intraspecific comparison of the tandem repeat elements of M. thalassophilus suggests mechanisms of concerted evolution leading to homogenization of the repetitive region. The utility of such an array of tandem repeats as a genetic marker for assessing population-level variability and evolutionary relationships among populations is discussed. Finally, the technical difficulties found in isolating the mtDNA CR in beetles are remarked upon.     

Genetic relationships among the shads (Alosa) revealed by mitochondrial DNA analysis

We surveyed restriction site differences in mitochondrial DNA. (mtDNA) among five species of shad ( Alosa ) from North America and Europe. Allis shad, Alosa alosa and twaite shad, Alosa fallax shared two divergent genotype groups, suggesting that the two forms are either a single species, or are distinct species that have hybridized. Phenetic and cladistic analyses of the relationships among the mitochondrial genotypes defined two groups of shad, corresponding to the subgenera, Alosa and Pomolobus . The mean estimated sequence divergence between the mtDNAs of these two groups of shad was 6.5%. Taken in conjunction with fossil data, this divergence estimate suggests that the rate of mtDNA divergence between the two subgenera has been almost 10-fold lower than the 'conventional' clock calibration for mtDNA.     

African DNA Lineages in the Mitochondrial Gene Pool of Europeans

Mitochondrial DNA (mtDNA) nucleotide sequences of African origin are found in various European populations at a low frequency (on average, less than 1%). Data on mtDNA variation in Eurasian and African populations have been analyzed, and African mtDNA lineages have been found in Europeans. It has been demonstrated that, despite the high diversity of mtDNA haplotypes of African origin in Europeans, few monophyletic clusters of African lineages are characterized by long-term diversity formed in Europe. Only two such mtDNA clusters (from haplogroups L1b and L3b) have been found, their evolutionary age not exceeding 6500 years. European and African populations have been compared with respect to the frequency distributions of the alleles of autosomal microsatellite loci found in Russian carriers of African mtDNA haplotypes. It has been demonstrated that alleles typical of Europeans are characteristic of the autosomal genotypes of these Russian individuals.     

Falling apart and merging: diversification of slender salamanders (Plethodontidae: Batrachoseps) in the American West

The plethodontid genus Batrachoseps , the slender salamanders, is the most diverse clade of salamanders in western North America, but it has posed taxonomic difficulties because it contains many morphologically cryptic species. A segment of the mitochondrial DNA gene cytochrome b was studied for 278 individuals densely sampled from throughout the range of all 18 described species and several undescribed species. Phylogenetic analyses of the mtDNA data identify six major clades, one corresponding to the subgenus Plethopsis and five within a monophyletic subgenus Batrachoseps. All major clades and most species within these clades display strong phylogeographic structuring. Comparisons of mtDNA and allozyme data show that several allozymically cohesive groups are not monophyletic with respect to mtDNA. We suggest that this phenomenon results from fragmentation of populations, divergence in allopatry, and then recontact and gradual merging of units caused predominantly by male-mediated gene flow. The mtDNA offers evidence that populations were once more isolated than they are now, while the patterns of allozyme variation reflect recent and current interactions among populations. The complex patterns of morphological, allozymic and mtDNA variation associated with the constantly changing geological landscape give insight into the nature of processes responsible for species formation in Batrachoseps .  © 2002 The Linnean Society of London. Biological Journal of the Linnean Society , 2002, 76 , 361–391.     

Sequences, sequence clusters and bacterial species

Whatever else they should share, strains of bacteria assigned to the same species should have house-keeping genes that are similar in sequence. Single gene sequences (or rRNA gene sequences) have very few informative sites to resolve the strains of closely related species, and relationships among similar species may be confounded by interspecies recombination. A more promising approach (multilocus sequence analysis, MLSA) is to concatenate the sequences of multiple house-keeping loci and to observe the patterns of clustering among large populations of strains of closely related named bacterial species. Recent studies have shown that large populations can be resolved into non-overlapping sequence clusters that agree well with species assigned by the standard microbiological methods. The use of clustering patterns to inform the division of closely related populations into species has many advantages for poorly studied bacteria (or to re-evaluate well-studied species), as it provides a way of recognizing natural discontinuities in the distribution of similar genotypes. Clustering patterns can be used by expert groups as the basis of a pragmatic approach to assigning species, taking into account whatever additional data are available (e.g. similarities in ecology, phenotype and gene content). The development of large MLSA Internet databases provides the ability to assign new strains to previously defined species clusters and an electronic taxonomy. The advantages and problems in using sequence clusters as the basis of species assignments are discussed.     

Barcoding Beetles: A Regional Survey of 1872 Species Reveals High Identification Success and Unusually Deep Interspecific Divergences

With 400 K described species, beetles (Insecta: Coleoptera) represent the most diverse order in the animal kingdom. Although the study of their diversity currently represents a major challenge, DNA barcodes may provide a functional, standardized tool for their identification. To evaluate this possibility, we performed the first comprehensive test of the effectiveness of DNA barcodes as a tool for beetle identification by sequencing the COI barcode region from 1872 North European species. We examined intraspecific divergences, identification success and the effects of sample size on variation observed within and between species. A high proportion (98.3%) of these species possessed distinctive barcode sequence arrays. Moreover, the sequence divergences between nearest neighbor species were considerably higher than those reported for the only other insect order, Lepidoptera, which has seen intensive analysis (11.99% vs up to 5.80% mean NN divergence). Although maximum intraspecific divergence increased and average divergence between nearest neighbors decreased with increasing sampling effort, these trends rarely hampered identification by DNA barcodes due to deep sequence divergences between most species. The Barcode Index Number system in BOLD coincided strongly with known species boundaries with perfect matches between species and BINs in 92.1% of all cases. In addition, DNA barcode analysis revealed the likely occurrence of about 20 overlooked species. The current results indicate that DNA barcodes distinguish species of beetles remarkably well, establishing their potential to provide an effective identification tool for this order and to accelerate the discovery of new beetle species.     

Phylogenetic relationships of Hemidactylus geckos from the Gulf of Guinea islands: patterns of natural colonizations and anthropogenic introductions estimated from mitochondrial and nuclear DNA sequences

Mitochondrial DNA (12S rRNA, 16S rRNA, and cytochrome b) sequences and nuclear sequences (C-mos and alpha-Enolase) were analyzed within all known Hemidactylus species from all three volcanic islands in the Gulf of Guinea that have never been connected to the continent. These comprise both endemic and widespread species. Our aim was to determine if the widespread species was introduced anthropogenically, to determine the number of distinct genetic lineages within the islands, and to determine if the endemic forms constituted a monophyletic group. Our results suggest that a previously undescribed species on S?o Tomé is the sister taxon to Hemidactylus newtoni, endemic to Annobon. Genetic variation between populations of Hemidactylus greefii from S?o Tomé and Principe is very high based on mtDNA sequences, but the forms cannot be distinguished using the nuclear DNA sequences. Hemidactylus mabouia appears to have been anthropogenically introduced to all three islands. The island endemics do not form a monophyletic group, suggesting multiple independent colonizations of the islands.     

An empirical genetic assessment of the severity of the northern elephant seal population bottleneck

A bottleneck in population size of a species is often correlated with a sharp reduction in genetic variation. The northern elephant seal (Mirounga angustirostris) has undergone at least one extreme bottleneck, having rebounded from 20-100 individuals a century ago to over 175,000 individuals today. The relative lack of molecular-genetic variation in contemporary populations has been documented, but the extent of variation before the late 19th century remains unknown. We have determined the nucleotide sequence of a 179 base-pair segment of the mitochondrial DNA (mtDNA) control region from seals that lived before, during and after a bottleneck low in 1892. A 'primerless' PCR was used to improve the recovery of information from older samples. Only two mtDNA genotypes were present in all 150+ seals from the 1892 bottleneck on, but we discovered four genotypes in five pre-bottleneck seals. This suggests a much greater amount of mtDNA genotypic variation before this bottleneck, and that the persistence of two genotypes today is a consequence of random lineage sampling. We cannot correlate the loss of mtDNA genotypes with a lowered mean fitness of individuals in the species today. However, we show that the species historically possessed additional genotypes to those present now, and that sampling of ancient DNA could elucidate the genetic consequences of severe reductions in population size.     

The genetic legacy of polyploid Bolivian Daphnia: the tropical Andes as a source for the North and South American D. pulicaria complex

We investigated genetic variation in asexual polyploid members of the water flea Daphnia pulex complex from a set of 12 Bolivian high-altitude lakes. We used nuclear microsatellite markers to study genetic relationships among all encountered multilocus genotypes, and combined this with a phylogenetic approach using DNA sequence data of three mitochondrial genes. Analyses of mitochondrial gene sequence divergence showed the presence of three very distinct clades that likely represent cryptic undescribed species. Our phylogenetic results suggest that the Daphnia pulicaria group, a complex of predominantly North American species that has diversified rapidly since the Pleistocene, has its origin in South America, as specific tests of topology indicated that all three South American lineages are ancestral to the North American members of this species group. A comparison between variation of nuclear and mitochondrial markers revealed that closely related polyploid nuclear genotypes sometimes belonged to very divergent mitochondrial lineages, while distantly related nuclear genotypes often belonged to the same mitochondrial lineage. This discrepancy suggests that these South American water fleas originated through reciprocal hybridization between different endemic, sexually reproducing parental lineages. It is also likely that polyploidy of the investigated lineages resulted from this hybridization. Nevertheless, no putative diploid parental lineages were found in the studied region.     

Comparison between RAPD and SSR molecular markers in detecting genetic variation in kiwifruit (Actinidia deliciosa A. Chev)

Two different DNA-based techniques, random amplified polymorphic DNA (RAPD) and simple sequence repeat (SSR) markers, were used for fingerprinting kiwifruit genotypes and for detecting undesirable genetic variation in micropropagated plants. The fragments were scored as present (1) or absent (0), and those readings were entered in a computer file as a binary matrix (one for each marker). Two cluster analyses were performed to express - in the form of dendrograms - the relationships among the genotypes and the genetic variability detected. Both DNA-based techniques were able to amplify all of the genotypes, but only SSR markers could detect genetic variation induced in micropropagated plants of cv. Tomuri. Two hypotheses were formulated to explain these results, both of them are in agreement with the results obtained using these two types of molecular markers. We conclude that when the tissue culture technique is used, the analysis of somaclonal variability could require more than one DNA-based technique; in fact, the genetic variation present in different sources could interfere or combine with the more or less polymorphic ability, as our results showed for SSR and RAPD markers.     

Molecular Population Genetics of Mtdna Size Variation in Crickets

Nucleotide sequence analysis of a region of cricket (Gryllus firmus) mtDNA showing discrete length variation revealed tandemly repeated sequences 220 base pairs (bp) in length. The repeats consist of 206 bp sequences bounded by the dyad symmetric sequence 5'GGGGGCATGCCCCC3'. The sequence data showed that mtDNA size variation in this species is due to variation in the number of copies of tandem repeats. Southern blot analysis was used to document the frequency of crickets heteroplasmic for two or more different-sized mtDNAs. In New England populations of G. firmus and a close relative Gryllus pennsylvanicus approximately 60% of the former and 45% of the latter were heteroplasmic. From densitometry of autoradiographs the frequencies of mtDNA size classes were determined for the population samples and are shown to very different in the two species. However, in populations where hybridization between the two species has occurred, the frequencies of size classes and cytoplasmic genotypes in each species' distinct mtDNA lineage were shifted in a manner suggesting nuclear-cytoplasmic interactions. The data were applied to reported diversity indices and hierarchical statistics. The hierarchical statistics indicated that the greatest proportion of variation for mtDNA size was due to variation among individuals in their cytoplasmic genotypes (heteroplasmic or homoplasmic state). The diversity indices were used to estimate a per-generation mutation rate for size variants of 10(-4). The data are discussed in light of the relationship between genetic drift and mutation in maintaining variation for mtDNA size.     

Bovine mtDNA D-loop haplotypes exceed mutations in number despite reduced recombination: an effective alternative for identity control

Mitochondrial (mt) DNA D-loop heterogeneity, haplotype distribution and possible sub-population structures within the relevant populations are important for DNA-based traceability. To gain insight into this distribution, we compared 1515 Bos taurus mtDNA D-loop sequences available from GenBank to 219 sequences that we sequenced de novo. A pronounced ambiguous trace typical of C-track length heteroplasmy was encountered in 5% of the samples, which were excluded from the analysis. Previously undescribed mutations and haplotypes were observed in 6% and 63% of the sequences, respectively. B. taurus haplotypes divided into the taurus, indicus and grunniens types and 302 variable sites formed the 858 taurus haplotypes detected. Fifty-five sites displayed a complex level of variation. As each level represents an independent mutation event, a total of 399 mutations were traced, which could potentially explain independent formation of less than half (47%) of the haplotypes encountered: most haplotypes were derived from different combinations of these mutations. We suggest that a mutational hotspot may explain these results and discuss the usefulness of mtDNA for identity and maternity assurance.