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

Integrative taxonomy of giant crested Eusirus in the Southern Ocean, including the description of a new species (Crustacea: Amphipoda: Eusiridae)

Among Antarctic amphipods of the genus Eusirus, a highly distinctive clade of giant species is characterized by a dorsal, blade-shaped tooth on pereionites 5–7 and pleonites 1–3. This lineage, herein named ‘crested Eusirus’, includes two potential species complexes, the Eusirus perdentatus and Eusirus giganteus complexes, in addition to the more distinctive Eusirus propeperdentatus. Molecular phylogenies and statistical parsimony networks (COI, CytB and ITS2) of crested Eusirus are herein reconstructed. This study aims to formally revise species diversity within crested Eusirus by applying several species delimitation methods (Bayesian implementation of the Poisson tree processes model, general mixed Yule coalescent, multi-rate Poisson tree processes and automatic barcode gap discovery) on the resulting phylogenies. In addition, results from the DNA-based methods are benchmarked against a detailed morphological analysis of all available specimens of the E. perdentatus complex. Our results indicate that species diversity of crested Eusirus is underestimated. Overall, DNA-based methods suggest that the E. perdentatus complex is composed of three putative species and that the E. giganteus complex includes four or five putative species. The morphological analysis of available specimens from the E. perdentatus complex corroborates molecular results by identifying two differentiable species, the genuine E. perdentatus and a new species, herein described as Eusirus pontomedon sp. nov.

Patterns of morphological diversification in giant Berberomeloe blister beetles (Coleoptera: Meloidae) reveal an unexpected taxonomic diversity concordant with mtDNA phylogenetic structure

Delimiting species boundaries is a complex challenge usually hindered by overlooked morphological diversification or misinterpretation of geographically structured phenotypic variability. Independent molecular data are extremely useful to characterize and understand such morphological diversity. Morphological and molecular variability of the non-phoretic and apterous, widely distributed, giant blister beetles of the genus Berberomeloe, were investigated within and between lineages across most of the distributional range of the genus. We used two mtDNA gene fragments to characterize genetic variability and to produce a time-calibrated phylogeny of the genus. Our results reveal several mitochondrial lineages, allopatrically, parapatrically and sympatrically distributed. Most clades are not distinguishable between each other based on morphometrics. However, no morphometric overlap is observed between two closely related clades, one of them occurring in sympatry with a distantly congeneric species (B. insignis), suggesting that sympatry could trigger morphological diversification. Although most species share a morphometric space, they can be morphologically identified by a combination of easily observed characteristic qualitative features. Based on the concordance between mtDNA clades and morphological units, we describe six new species of Berberomeloe (B. castuo sp. nov., B. comunero sp. nov., B. indalo sp. nov, B. yebli sp. nov., B. payoyo sp. nov. and B. tenebrosus sp. nov.), revalidate two taxa (B. maculifrons comb. nov. and B. laevigatus comb. nov.) and redefine B. majalis.

Delineation of two new, highly similar species of Taiwanese Cylindera tiger beetles (Coleoptera, Carabidae, Cicindelinae) based on morphological and molecular evidence

Tiger beetles have been recognized primarily based on morphological characters. However, the variations of elytral maculation and coloration sometimes lead to misjudgment in species classification and the overlooking of the existence of cryptic species. Recently, specimens of two endemic species of Taiwanese Cylindera, C.sauteri and C.pseudocylindriformis, exhibit morphologically recognizable forms, indicating that some undescribed species may exist. To clarify their taxonomic status, morphological characteristics including male genitalia were examined and two mitochondrial genes, COI and 16S rDNA, and one nuclear 28S rDNA were analyzed. Molecular phylogenetic inferences indicated that both forms in both species are reciprocally monophyletic. Moreover, molecular dating showed the forms diverged approximately 1.3 million years ago. Two new species, Cylinderaooasp. nov. and Cylinderaautumnalissp. nov., are thereby described. The main recognizable characteristics separating C.ooasp. nov. from C.sauteri are the lack of a triangular spot at the middle edge of elytron and the elongated but not rounded subapical spot. For C.autumnalissp. nov., the apical lunula near the elytral suture is thickened but not linear and slender, and its elytra are more metallic brownish than those of C.pseudocylindriformis. Although their aedeagi characteristics are not distinctive, the body size of the proposed two new species is different. Field observation revealed that niche utilization would be relevant for differentiating these closely related species.

Pleistocene phylogeography and cryptic diversity of a tiger beetle, Calomera littoralis, in North-Eastern Mediterranean and Pontic regions inferred from mitochondrial COI gene sequences

Background. Calomera littoralis is a Palearctic species, widely distributed in Europe; inhabiting predominantly its Atlantic, Mediterranean and Black Sea coastlines. Methods. Its phylogeography on the Balkan Peninsula and on the north-western Black Sea coast was inferred using a 697 bp long portion of the mitochondrial COI gene, amplified from 169 individuals collected on 43 localities. Results. The results revealed two genetically divergent groups/lineages, the southern one inhabiting both the Balkan Peninsula and the Pontic Region and the northern one found exclusively in the Pontic Region. Species delimitation based on DNA barcoding gap suggested an interspecific level of divergence between these groups. Multivariate analysis of eight male and female morphometric traits detected no difference between the groups, implying they may represent cryptic species. The Bayesian time-calibrated reconstruction of phylogeny suggested that the lineages diverged ca. 2.3 Ma, in early Pleistocene. Discussion. The presence of the two genetically divergent groups results most likely from contemporary isolation of the Pontic basin from the Mediterranean that broke the continuous strip of coastal habitats inhabited by C. littoralis. Demographic analyses indicated that both lineages have been in demographic and spatial expansion since ca. 0.15 Ma. It coincides with the terminal stage of MIS-6, i.e., Wartanian/Saalian glaciation, and beginning of MIS-5e, i.e., Eemian interglacial, during which, due to eustatic sea level rise, a wide connection between Mediterranean and the Pontic basin was re-established. This, along with re-appearance of coastal habitats could initiate north-east expansion of the southern lineage and its secondary contact with the northern one. The isolation of the Pontic basin from the Mediterranean during the Weichselian glaciation most likely did not have any effect on their phylogeography.

Molecular taxonomic analysis of the plant associations of adult pollen beetles (Nitidulidae: Meligethinae), and the population structure of Brassicogethes aeneus

Pollen beetles (Nitidulidae: Meligethinae) are among the most abundant flower-visiting insects in Europe. While some species damage millions of hectares of crops annually, the biology of many species is little known. We assessed the utility of a 797 base pair fragment of the cytochrome oxidase 1 gene to resolve molecular operational taxonomic units (MOTUs) in 750 adult pollen beetles sampled from flowers of 63 plant species sampled across the UK and continental Europe. We used the same locus to analyse region-scale patterns in population structure and demography in an economically important pest, Brassicogethes aeneus. We identified 44 Meligethinae at ∼2% divergence, 35 of which contained published sequences. A few specimens could not be identified because the MOTUs containing them included published sequences for multiple Linnaean species, suggesting either retention of ancestral haplotype polymorphism or identification errors in published sequences. Over 90% of UK specimens were identifiable as B. aeneus. Plant associations of adult B. aeneus were found to be far wider taxonomically than for their larvae. UK B. aeneus populations showed contrasting affiliations between the north (most similar to Scandinavia and the Baltic) and south (most similar to western continental Europe), with strong signatures of population growth in the south.

The tiger beetles (Coleoptera, Carabidae, Cicindelinae) of Israel and adjacent lands

Based on field studies, museums collections and literature sources, the current knowledge of the tiger beetle fauna of Israel and adjacent lands is presented. In Israel eight species occur, one of them with two subspecies, while in the Sinai Peninsula nine species of tiger beetles are now known. In the combined regions seven genera from two tribes were found. The Rift Valley with six cicindelids species is the most specious region of Israel. Cylinderacontortavaldenbergi and Cicindelajavetiazari have localized distributions and should be considered regional endemics. A similarity analysis of the tiger beetles faunas of different regions of Israel and the Sinai Peninsula reveal two clusters of species. The first includes the Great Rift Valley and most parts of the Sinai Peninsula, and the second incorporates most regions of Israel together with Central Sinai Foothills. Five distinct adult phenological groups of tiger beetles can be distinguished in these two clusters: active all-year (three species), spring-fall (five species), summer (two species), spring-summer (one species) and spring (one species). The likely origins of the tiger beetle fauna of this area are presented. An annotated list and illustrated identification key of the Cicindelinae of Israel and adjacent lands are provided.

Salt lakes of La Mancha (Central Spain): A hot spot for tiger beetle (Carabidae, Cicindelinae) species diversity

The tiger beetle assemblage of the wetlands of La Mancha (central Spain) comprises nine species: Calomera littoralis littoralis, Cephalota maura maura, Cephalota circumdata imperialis, Cephalota dulcinea, Cicindela campestris campestris, Cicindela maroccana, Cylindera paludosa, Lophyra flexuosa flexuosa, and Myriochila melancholica melancholica. This assemblage represents the largest concentration of tiger beetles in a single 1º latitude / longitude square in Europe. General patterns of spatial and temporal segregation among species are discussed based on observations of 1462 specimens registered during an observation period of one year, from April to August. The different species of Cicindelini appear to be distributed over space and time, with little overlapping among them. Three sets of species replace each other phenologically as the season goes on. Most of the species occupy drying or dried salt lakes and salt marshes, with sparse vegetation cover. Spatial segregation is marked in terms of substrate and vegetation use. Calomera littoralis and Myriochila melancholica have been observed mainly on wet soils; Cephalota circumdata on dry open saline flats; Cephalota dulcinea and Cylindera paludosa in granulated substrates with typical halophytic vegetation; Cephalota maura is often present in man-modified areas. Cephalota circumdata and Cephalota dulcinea are included as species of special interest in the list of protected species in Castilla-La Mancha. Conservation problems for the Cicindelini assemblage arise from agricultural activities and inadequate use of sport vehicles. Attempts at restoring the original habitat, supressing old semi-industrial structures, may affect the spatial heterogeneity of the lakes, and have an effect on Cicindelinae diversity.

Exploring the Leaf Beetle Fauna (Coleoptera: Chrysomelidae) of an Ecuadorian Mountain Forest Using DNA Barcoding

Background

Tropical mountain forests are hotspots of biodiversity hosting a huge but little known diversity of insects that is endangered by habitat destruction and climate change. Therefore, rapid assessment approaches of insect diversity are urgently needed to complement slower traditional taxonomic approaches. We empirically compare different DNA-based species delimitation approaches for a rapid biodiversity assessment of hyperdiverse leaf beetle assemblages along an elevational gradient in southern Ecuador and explore their effect on species richness estimates.

Methodology/Principal Findings

Based on a COI barcode data set of 674 leaf beetle specimens (Coleoptera: Chrysomelidae) of 266 morphospecies from three sample sites in the Podocarpus National Park, we employed statistical parsimony analysis, distance-based clustering, GMYC- and PTP-modelling to delimit species-like units and compared them to morphology-based (parataxonomic) species identifications. The four different approaches for DNA-based species delimitation revealed highly similar numbers of molecular operational taxonomic units (MOTUs) (n = 284–289). Estimated total species richness was considerably higher than the sampled amount, 414 for morphospecies (Chao2) and 469–481 for the different MOTU types. Assemblages at different elevational levels (1000 vs. 2000 m) had similar species numbers but a very distinct species composition for all delimitation methods. Most species were found only at one elevation while this turnover pattern was even more pronounced for DNA-based delimitation.

Conclusions/Significance

Given the high congruence of DNA-based delimitation results, probably due to the sampling structure, our study suggests that when applied to species communities on a regionally limited level with high amount of rare species (i.e. ~50% singletons), the choice of species delimitation method can be of minor relevance for assessing species numbers and turnover in tropical insect communities. Therefore, DNA-based species delimitation is confirmed as a valuable tool for evaluating biodiversity of hyperdiverse insect communities, especially when exact taxonomic identifications are missing.

Molecular phylogenetic analysis of the Amiota taurusata species group within the Chinese species, with descriptions of two new species

The relationships among six species of the Amiota taurusata Takada, Beppu, & Toda (Diptera: Drosophilidae) species group were investigated based on DNA sequence data of the mitochondrial NADH dehydrogenase subunit 2 ( ND2) gene, using three species of the genus Amiota as outgroups. A mitochondrial gene, cytochrome c oxidase I ( COI), can be used to discriminate between species of the taurusata group. Two new species are described from South China: A. protuberantis Shao et Chen, sp. nov. and A. shennongi Shao et Chen, sp. nov. A key to all the species of the taurusata group based on morphological characters is provided.

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.

A long-living species of the hydrophiloid beetles: Helophorus sibiricus from the early Miocene deposits of Kartashevo (Siberia, Russia)

The recent hydrophiloid species Helophorus (Gephelophorus) sibiricus (Motschulsky, 1860) is recorded from the early Miocene deposits of Kartashevo assigned to the Ombinsk Formation. A detailed comparison with recent specimens allowed a confident identification of the fossil specimen, which is therefore the oldest record of a recent species for the Hydrophiloidea. The paleodistribution as well as recent distribution of the species is summarized, and the relevance of the fossil is discussed. In addition, the complex geological settings of the Kartashevo area are briefly summarized.

Exploring diversity in cryptorhynchine weevils (Coleoptera) using distance-, character- and tree-based species delineation

Species boundaries are studied in a group of beetles, the western Palaearctic Cryptorhynchinae. We test for congruence of 'traditionally' identified morphospecies with species inferred through parsimony networks, distance-based clustering and the ultrametric tree-based generalized mixed yule-coalescent (GMYC) approach. For that purpose, we sequenced two variable fragments of mitochondrial DNA (CO1 and 16S) for a total of 791 specimens in 217 species of Cryptorhynchinae. Parsimony networks, morphology-calibrated distance clusters and the different tree-based species inferences all achieved low congruence with morphospecies, at best 60%. Although the degree of match with morphospecies was often similar for the different approaches, the composition of clusters partially varied. A barcoding gap was absent in morphospecies-oriented distances as well as for GMYC species clusters. This demonstrates that not only erroneous taxonomic assignments, incomplete lineage sorting, hybridization, or insufficient sampling can compromise distance-based identification, but also differences in speciation rates and uneven tree structure. The initially low match between morphospecies and the different molecular species delineation methods in this case study shows the necessity of combining the output of various methods in an integrative approach. Thereby we obtain an idea about the reliability of the different results and signals, which enables us to fine-tune sampling, delineation technique and data collection, and to identify species that require taxonomic revision.

Use of a mitochondrial COI sequence to identify species of the subtribe Aphidina (Hemiptera, Aphididae)

Aphids of the subtribe Aphidina are found mainly in the North Temperate Zone. The relative lack of diagnostic morphological characteristics has hindered the identification of species in this group. However, DNA-based taxonomic methods can clarify species relationships within this group. Sequence variation in a partial segment of the mitochondrial COI gene was highly effective for identifying species within Aphidina. Thirty-six species of Aphidina were identified in a neighbor-joining tree. Mean intraspecific sequence divergence in Aphidina was 0.52%, with a range of 0.00% to 2.95%, and the divergences of most species were less than 1%. Mean interspecific divergence within previously recognized genera or morphologically similar species groups was 6.80%, with a range of 0.68% to 11.40%, with variation mainly in the range of 3.50% to 8.00%. Possible reasons for anomalous levels of mean nucleotide divergence within or between some taxa are discussed.

Species delimitation and evolution in morphologically and chemically diverse communities of the lichen-forming genus Xanthoparmelia (Parmeliaceae, Ascomycota) in western North America

PREMISE OF THE STUDY

Accurate species delimitation is important for understanding the diversification of biota and has critical implications for ecological and conservation studies. However, a growing body of evidence indicates that morphology-based species circumspection in lichenized fungi misrepresents fungal diversity. The foliose lichen genus Xanthoparmelia includes over 800 species displaying a complex array of morphological and secondary metabolite diversity.

METHODS

We used a multifaceted approach, applying phylogenetic, population genetic, and genealogical analyses to delimit species in a single well-supported monophyletic clade containing 10 morphologically and chemically diverse Xanthoparmelia species in western North America. Sequence data from four ribosomal and two low-copy, protein-coding markers, along with chemical and morphological data were used to assess species diversity.

KEY RESULTS

We found that traditionally circumscribed species are not supported by molecular data. Rather, all sampled taxa were better represented by three polymorphic population clusters. Our results suggest that secondary metabolite variation may have limited utility in diagnosing lineages within this group, while identified populations clusters did not reflect major phylogeographic or ecological patterns.

CONCLUSIONS

In contrast to studies revealing previously undiscovered fungal lineages masked within lichen species circumscribed by traditional morphological and chemical concepts, the present study suggests that species diversity has been overestimated in the species-rich genus Xanthoparmelia.

Island biogeography of Galápagos lava lizards (Tropiduridae: Microlophus): species diversity and colonization of the archipelago

The "lava lizards" (Microlophus) are distributed throughout the Galápagos Archipelago, and consist of radiations derived from two independent colonizations. The "Eastern Radiation" includes M. bivittatus and M. habeli endemic to San Cristobal and Marchena Islands. The "Western Radiation" includes five to seven historically recognized species distributed across almost the entire Archipelago. We combine dense geographic sampling and multilocus sequence data to estimate a phylogenetic hypothesis for the Western Radiation, to delimit species boundaries in this radiation, and to estimate a time frame for colonization events. Our phylogenetic hypothesis rejects two earlier topologies for the Western Radiation and paraphyly of M. albemarlensis, while providing strong support for single colonizations on each island. The colonization history implied by our phylogeny is consistent with general expectations of an east-to-west route predicted by the putative age of island groups, and prevailing ocean currents in the Archipelago. Additionally, combined evidence suggests that M. indefatigabilis from Santa Fe should be recognized as a full species. Finally, molecular divergence estimates suggest that the two colonization events likely occurred on the oldest existing islands, and the Western Radiation represents a recent radiation that, in most cases, has produced species that are considerably younger than the islands they inhabit.

Molecular diversity assessment of arctic and boreal Agaricus taxa

We provide a phylogenetic diversity assessment study in genus Agaricus as part of our ongoing work to saturate ITS and LSU rDNA sequence diversity of soil-dwelling fungi in Alaska. Pairwise sequence similarity-based groupings and statistical parsimony analyses were applied to delimit operational taxonomic unit (OTU) and were compared to results of full phylogenetic analyses. Our results show that the proportion of section Arvenses taxa is particularly high in the boreal forest and hypo-arctic (low arctic) regions, whereas the genus is represented by section Agaricus in high arctic habitats. Furthermore our findings suggest that the commercially important A. bisporus occurs naturally in the boreal region of interior Alaska, substantially expanding the known northern limit of the species. Delimitations of OTU varied greatly with different methods. In general 95% similarity-based grouping proved to be the least sensitive method, often resulting in section- and subsection-level groups. The 95% connection-limit statistical parsimony separated far more groups. The 98% similarity-based groups and the 98% connection limit networks recognized respectively 11 and 13 OTU containing our specimens. The 98% connection limit statistical parsimony was the only method in which all recognized OTU consisted of members grouped by branches with significant (> .95) posterior probabilities, providing an independent support for the groups. Our results also point out that considerable additional efforts will be needed to elucidate the evolution of this diverse genus and to assess its phylogenetic diversity, given that most taxa in our analyses could not be placed convincingly within well characterized species using ITS/LSU data.

Species concepts and species delimitation

The issue of species delimitation has long been confused with that of species conceptualization, leading to a half century of controversy concerning both the definition of the species category and methods for inferring the boundaries and numbers of species. Alternative species concepts agree in treating existence as a separately evolving metapopulation lineage as the primary defining property of the species category, but they disagree in adopting different properties acquired by lineages during the course of divergence (e.g., intrinsic reproductive isolation, diagnosability, monophyly) as secondary defining properties (secondary species criteria). A unified species concept can be achieved by treating existence as a separately evolving metapopulation lineage as the only necessary property of species and the former secondary species criteria as different lines of evidence (operational criteria) relevant to assessing lineage separation. This unified concept of species has several consequences for species delimitation, including the following: First, the issues of species conceptualization and species delimitation are clearly separated; the former secondary species criteria are no longer considered relevant to species conceptualization but only to species delimitation. Second, all of the properties formerly treated as secondary species criteria are relevant to species delimitation to the extent that they provide evidence of lineage separation. Third, the presence of any one of the properties (if appropriately interpreted) is evidence for the existence of a species, though more properties and thus more lines of evidence are associated with a higher degree of corroboration. Fourth, and perhaps most significantly, a unified species concept shifts emphasis away from the traditional species criteria, encouraging biologists to develop new methods of species delimitation that are not tied to those properties.

Things fall apart: biological species form unconnected parsimony networks

The generality of operational species definitions is limited by problematic definitions of between-species divergence. A recent phylogenetic species concept based on a simple objective measure of statistically significant genetic differentiation uses between-species application of statistical parsimony networks that are typically used for population genetic analysis within species. Here we review recent phylogeographic studies and reanalyse several mtDNA barcoding studies using this method. We found that (i) alignments of DNA sequences typically fall apart into a separate subnetwork for each Linnean species (but with a higher rate of true positives for mtDNA data) and (ii) DNA sequences from single species typically stick together in a single haplotype network. Departures from these patterns are usually consistent with hybridization or cryptic species diversity.

Delimiting cohesion species: extreme population structuring and the role of ecological interchangeability

Species exhibiting morphological homogeneity and strong population structuring present challenging taxonomic problems: morphology-based approaches infer few species, whereas genetic approaches often indicate more. Morphologically cryptic, yet genetically divergent species groups require alternative approaches to delimiting species that assess adaptive divergence and ecological interchangeability of lineages. We apply such an approach to Promyrmekiaphila, a small genus (three nominal taxa) of trapdoor spiders endemic to northern California to define cohesion species (lineages that are genetically exchangeable and ecologically interchangeable). Genetic exchangeability is evaluated using standard phylogeographical techniques (e.g. nested clade analysis); ecological interchangeability is assessed using two GIS-based approaches. First, climatic values are extracted from layer data for each locality point and utilized in a principal components analysis followed by MANOVA. Second, niche-based distribution models of genetically divergent lineages are created using a maximum-entropy modelling approach; the amount of overlap among lineages is calculated and evaluated against a probability distribution of null overlap. Lineages that have significant amounts of predicted overlap are considered ecologically interchangeable. Based on a synthetic evaluation of ecological interchangeability, geographical concordance, and morphological differentiation, we conclude that Promyrmekiaphila comprises six cohesion species, five of which are cryptic (i.e. undetectable by conventional means).

Population structure and species boundary delimitation of crypticDioryctriamoths: an integrative approach

Accurate delimitation of species boundaries is especially important in cryptic taxa where one or more character sources are uninformative or are in conflict. Rather than relying on a single marker to delimit species, integrative taxonomy uses multiple lines of evidence such as molecular, morphological, behavioural and geographic characters to test species limits. We examine the effectiveness of this approach by testing the delimitation of two cryptic Nearctic species of Dioryctria (Lepidoptera: Pyralidae) using three independent molecular markers [cytochrome c oxidase I (COI), second internal transcribed spacer unit (ITS2), and elongation factor 1alpha (EF1alpha)], forewing variation and larval host plant association. Although mitochondrial DNA (mtDNA) haplotypes do not form reciprocally monophyletic clades, restricted gene flow between COI haplotype groups, and concordance with ITS2 genotypes, forewing variation and host plant associations support delimitation of two Nearctic species: eastern Dioryctria reniculelloides and western Dioryctria pseudotsugella. Conversely, EF1alpha genotype variation was incongruent with the two previous markers. A case of discordance between COI and ITS2 was detected, suggesting either introgression due to hybridization or retained ancestral polymorphism due to incomplete coalescence. This study is consistent with other similar literature where molecular loci in closely related species progress from shared to fixed haplotypes/alleles, and from polyphyletic to reciprocally monophyletic relationships, although loci may vary in these characteristics despite maintenance of genomic integrity between distinct species. In particular, mtDNA in other studies generally showed a lower rate of fixation of differences than did X-linked or autosomal loci, reinforcing the need to use an integrative approach for delimiting species.

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.