Taxonomy as a hypothesis: testing the status of the Bermuda buckeye butterfly Junonia coenia bergi (Lepidoptera: Nymphalidae)

Species determination and definition in eukaryotes have traditionally been based on morphology, with little focus on genetic differentiation. Molecular methods allow for the independent assessment of morphology‐based taxonomic hypotheses. Three criteria used to define a full species for taxonomic purposes are morphological distinction, formation of a monophyletic lineage, and reproductive isolation. Junonia butterflies (Nymphalidae) are becoming an important experimental model system, but the taxonomy of many New World Junonia species is unclear. One of these species is J. coenia, which contains the subspecies J. coenia coenia, J. coenia grisea and J. coenia bergi. Previous studies suggest that J. coenia grisea may meet the criteria for full species status. Therefore, we evaluated the geographically isolated and rarely studied Bermuda buckeye butterfly J. coenia bergi to determine if it was similarly distinct. Physical examination of specimens and phylogenetic and population genetic analyses of mitochondrial cytochrome c oxidase subunit I, nuclear wingless, and complete mitochondrial genome sequences suggest that while J. coenia bergi is smaller in body size than many Junonia and has distinctive ventral hindwing colouration, it does not form a monophyletic lineage and shows indications of continued gene flow with North American mainland J. coenia coenia populations. Thus, J. coenia bergi does not meet the criteria for full species designation, but geographic isolation, morphological distinctiveness, and cultural importance suggest that it remain recognized as a subspecies of J. coenia. Similar analyses will be useful for addressing further taxonomic questions in Junonia and other taxa, especially where morphology‐based taxonomic determinations are ambiguous.     

Getting western: biogeographical analysis of morphological variation,mitochondrial haplotypes and nuclear markers reveals cryptic species and hybrid zones in the Junonia butterflies of the American southwest and Mexico

The American southwest and northern Mexico has a great degree of endemic diversity compared with the rest of North America. The Pleistocene glaciations and the dispersal of species from glacial refuges in this region have been important engines for the production of biodiversity in the region. The New World Junonia are a recent radiation of butterflies that are thought to have spent time in these refuges during periods of glacial advancement. We have reconstructed the plausible movements and the contemporary geographic distributions of the five taxa (J. coenia, J. grisea, J. litoralis, J. nigrosuffusa and J. zonalis) that occur in the American southwest and northern Mexico using phenotypic and genotypic information primarily from specimens preserved in museum collections, supplemented with additional contemporary collections. Evidence of cryptic species and hybridization events were observed using mitochondrial haplotypes, genotypic variation at the nuclear wingless locus, multilocus DNA fingerprinting, patterns of larval host plant use, variation in life‐history traits and morphological characteristics. Based on these lines of evidence, and in spite of low levels of hybridization between them, we argue that all five Junonia taxa are independent evolutionary lineages. Junonia grisea and J. coenia are morphologically very similar, but differences in morphology, life‐history traits, nuclear genotypes and mitochondrial haplotypes suggest that they are a cryptic species pair, thus elevating J. grisea comb.n. to a full species when it had previously been considered to be a subspecies of J. coenia.     

Phylogenetic relationships and historical biogeography of tribes and genera in the subfamily Nymphalinae (Lepidoptera: Nymphalidae)

We infer for the first time the phylogenetic relationships of genera and tribes in the ecologically and evolutionarily well‐studied subfamily Nymphalinae using DNA sequence data from three genes: 1450 bp of cytochrome oxidase subunit I (COI) (in the mitochondrial genome), 1077 bp of elongation factor 1‐alpha (EF1‐α) and 400–403 bp of wingless (both in the nuclear genome). We explore the influence of each gene region on the support given to each node of the most parsimonious tree derived from a combined analysis of all three genes using Partitioned Bremer Support. We also explore the influence of assuming equal weights for all characters in the combined analysis by investigating the stability of clades to different transition/transversion weighting schemes. We find many strongly supported and stable clades in the Nymphalinae. We are also able to identify ‘rogue’ taxa whose positions are weakly supported (the different gene regions are in conflict with each other) and unstable. Our main conclusions are: (1) the tribe Coeini as currently constituted is untenable, and Smyrna, Colobura and Tigridia are part of Nymphalini; (2) ‘Kallimini’ is paraphyletic with regard to Melitaeini and should be split into three tribes: Kallimini s.s., Junoniini and Victorinini; (3) Junoniini, Victorinini, Melitaeini and the newly circumscribed Nymphalini are strongly supported monophyletic groups, and (4) Precis and Junonia are not synonymous or even sister groups. The species Junonia coenia, a model system in developmental biology, clearly belongs in the genus Junonia. A dispersal‐vicariance analysis suggests that dispersal has had a major effect on the distributions of extant species, and three biotic regions are identified as being centres of diversification of three major clades: the Palaearctic for the Nymphalis‐group, the Afrotropics for Junoniini and the Nearctic for Melitaeini. © 2005 The Linnean Society of London, Biological Journal of the Linnean Society, 2005, 86 , 227–251.     

Development of internal COI primers to improve and extend barcoding of fruit flies (Diptera: Tephritidae: Dacini)

Accurate species-level identifications underpin many aspects of basic and applied biology;however,identifications can be hampered by a lack of discriminating morphological characters,taxonomic expertise or time.Molecular approaches,such as DNA"barcoding"of the cytochrome c oxidase(COI)gene,are argued to overcome these issues.However,nuclear encoding of mitochondrial genes(numts)and poor amplification success of suboptimally preserved specimens can lead to erroneous identifications.One insect group for which these molecular and morphological problems are significant are the dacine fruit flies(Diptera:Tephritidae:Dacini).We addressed these issues associated with COI barcoding in the dacines by first assessing several"universal"COI primers against public mitochondrial genome and numt sequences for dacine taxa.We then modified a set of four primers that more closely matched true dacine COI sequence and amplified two overlapping portions of the COI barcode region.Our new primers were tested alongside universal primers on a selection of dacine species,including both fresh preserved and decades-old dry specimens.Additionally,Bactrocera tiyoni mitochondrial and nuclear genomes were compared to identify putative numts.Four numt clades were identified,three of which were amplified using existing universal primers.In contrast,our new primers preferentially amplified the"true"mitochondrial COI barcode in all dacine species tested.The new primers also successfully amplified partial barcodes from dry specimens for which full length barcodes were unobtainable.Thus we recommend these new primers be incorporated into the suites of primers used by diagnosticians and quarantine labs for the accurate identification of dacine species.     

Genetic study reveals close link between Irish and Northern Spanish specimens of the protected Lusitanian slug Geomalacus maculosus

The slug Geomalacus maculosus is a prominent member of the Lusitanian fauna. As its global distribution is restricted to western Ireland and northern Iberia, it is protected under EU legislation. Nothing is known about the genetic variability and population structure of this species, so, with a special view to shedding light on the origin of the Irish G. maculosus, tissue samples from 78 specimens were collected from 13 locations within Ireland and ten locations within Iberia and partial sequences of the mitochondrial 16S rRNA and cytochrome oxidase subunit 1 (COI) and from the nuclear internal transcribed spacer 1 region (ITS‐1) were compared. The genetic diversity of the Irish G. maculosus was found to be greatly reduced compared with the Iberian populations, with only one (16S rRNA) and two (COI) mitochondrial haplotypes identified respectively. No private Irish ITS‐1 haplotype was found. Based on the COI sequences, the Irish specimens clustered closest to Spanish specimens from Northern Asturias and Cantabria, and the bGMYC analysis identified five further Iberian clades that were highly genetically differentiated suggesting long‐term allopatric divergence.     

DNA barcoding commercially important fish species of Turkey

DNA barcoding was used in the identification of 89 commercially important freshwater and marine fish species found in Turkish ichthyofauna. A total of 1765 DNA barcodes using a 654‐bp‐long fragment of the mitochondrial cytochrome c oxidase subunit I gene were generated for 89 commercially important freshwater and marine fish species found in Turkish ichthyofauna. These species belong to 70 genera, 40 families and 19 orders from class Actinopterygii, and all were associated with a distinct DNA barcode. Nine and 12 of the COI barcode clusters represent the first species records submitted to the BOLD and GenBank databases, respectively. All COI barcodes (except sequences of first species records) were matched with reference sequences of expected species, according to morphological identification. Average nucleotide frequencies of the data set were calculated as T = 29.7%, C = 28.2%, A = 23.6% and G = 18.6%. Average pairwise genetic distance among individuals were estimated as 0.32%, 9.62%, 17,90% and 22.40% for conspecific, congeneric, confamilial and within order, respectively. Kimura 2‐parameter genetic distance values were found to increase with taxonomic level. For most of the species analysed in our data set, there is a barcoding gap, and an overlap in the barcoding gap exists for only two genera. Neighbour‐joining trees were drawn based on DNA barcodes and all the specimens clustered in agreement with their taxonomic classification at species level. Results of this study supported DNA barcoding as an efficient molecular tool for a better monitoring, conservation and management of fisheries.     

A tale of two haplotype groups: evaluating the New World Junonia ring species hypothesis using the distribution of divergent COI haplotypes

The New World Junonia butterflies are a possible ring species with a circum‐Caribbean distribution. Previous reports suggest a steady transition between North and South American forms in Mesoamerica, but in Cuba the forms were thought to co‐exist without interbreeding representing the overlapping ends of the ring. Three criteria establish the existence of a ring species: a ring‐shaped geographic distribution, gene flow among intervening forms and genetic isolation in the region of range overlap. We evaluated mitochondrial cytochrome oxidase I haplotypes in Junonia from nine species in the Western Hemisphere to test the Junonia ring species hypothesis. Junonia species are generally not monophyletic with respect to COI haplotypes, which are shared across species. However, two major COI haplotype groups exist. Group A predominates in South America, and Group B predominates in North and Central America. Therefore, COI haplotypes can be used to assess the degree of genetic influence a population receives from each continent. Junonia shows a ring‐shaped distribution around the Caribbean, and evidence is consistent with gene flow among forms of Junonia, including those from Mesoamerica. However, we detected no discontinuity in gene flow in Cuba or elsewhere in the Caribbean consistent with genetic isolation in the region of overlap. Although sampling is still very limited in the critical region, the only remaining possibility for a circum‐Caribbean discontinuity in gene flow is at the Isthmus of Panama, where there may be a transition from 98% Group B haplotypes in Costa Rica to 85–100% Group A haplotypes in South America.     

Speciation in North American Junonia from a genomic perspective

Delineating species boundaries in phylogenetic groups undergoing recent radiation is a daunting challenge akin to discretizing continuity. Here, we propose a general approach exemplified by American butterflies from the genus Junonia Hübner notorious for the variety of similar phenotypes, ease of hybridization, and the lack of consensus about their classification. We obtain whole-genome shotgun sequences of about 200 specimens. We reason that discreteness emerges from continuity by means of a small number of key players, and search for the proteins that diverged markedly between sympatric populations of different species, while keeping low polymorphism within these species. Being 0.25% of the total number, these three dozen ‘speciation’ proteins indeed partition pairs of Junonia populations into two clusters with a prominent break in between, while all proteins taken together fail to reveal this discontinuity. Populations with larger divergence from each other, comparable to that between two sympatric species, form the first cluster and correspond to different species. The other cluster is characterized by smaller divergence, similar to that between allopatric populations of the same species and comprise conspecific pairs. Using this method, we conclude that J. genoveva (Cramer), J. litoralis Brévignon, J. evarete (Cramer), and J. divaricata C. & R. Felder are restricted to South America. We find that six species of Junonia are present in the United States, one of which is new: Junonia stemosa Grishin, sp.n. (i), found in south Texas and phenotypically closest to J. nigrosuffusa W. Barnes & McDunnough (ii) in its dark appearance. In the pale nudum of the antennal club, these two species resemble J. zonalis C. & R. Felder (iii) from Florida and the Caribbean Islands. The pair of sister species, J. grisea Austin & J. Emmel (iv) and J. coenia Hübner (v), represent the classic west/east U.S.A. split. The mangrove feeder (as caterpillar), dark nudum J. neildi Brévignon (vi) enters south Texas as a new subspecies Junonia neildi varia Grishin ssp.n. characterized by more extensive hybridization with and introgression from J. coenia, and, as a consequence, more variable wing patterns compared with the nominal J. n. neildi in Florida. Furthermore, a new mangrove-feeding species from the Pacific Coast of Mexico is described as Junonia pacoma Grishin sp.n. Finally, genomic analysis suggests that J. nigrosuffusa may be a hybrid species formed by the ancestors of J. grisea and J. stemosa sp.n. This published work has been registered on Zoobank, http://zoobank.org/urn:lsid:zoobank.org:act:C6E70C96-463D-4E6A-95CC-B0384B0EEEBA .     

Darker female pigeons transmit more specific antibodies to their eggs than do paler ones

Melanin‐based coloration is widespread among vertebrates, yet the adaptive significance of such pigments remains elusive, particularly with regard to the link between melanin and immune‐mediated maternal effects. The aim of this study was to investigate whether melanin‐based coloration could signal the ability of mothers to mount a humoral response and to transfer maternal antibodies (Ab) to their young. We injected differently coloured (pale and dark) female feral pigeons (Columba livia) with Chlamydiae (a natural antigen) and Keyhole Limpet Haemocyanin (KLH, an artificial antigen), and found no significant difference in humoral response between differently coloured females. However, darker females transferred more Ab against Chlamydiae into their eggs than paler ones, despite similar circulating levels of Ab. In addition to this, melanin‐based coloration showed a high heritability value. This suggests that a genetically based coloured trait might be linked to the ability of females to transfer specific Ab against Chlamydiae (but not against KLH) to their offspring, independent of their ability to produce Ab. This suggests that transmission of maternal Ab is antigen dependent, and that melanin‐based coloration might signal female ability to transmit specific Ab against natural pathogens. © 2013 The Linnean Society of London     

Molecular identification of mosquitoes (Diptera: Culicidae) in southeastern Australia

DNA barcoding is a modern species identification technique that can be used to distinguish morphologically similar species, and is particularly useful when using small amounts of starting material from partial specimens or from immature stages. In order to use DNA barcoding in a surveillance program, a database containing mosquito barcode sequences is required. This study obtained Cytochrome Oxidase I (COI) sequences for 113 morphologically identified specimens, representing 29 species, six tribes and 12 genera; 17 of these species have not been previously barcoded. Three of the 29 species ─ Culex palpalis, Macleaya macmillani, and an unknown species originally identified as Tripteroides atripes ─ were initially misidentified as they are difficult to separate morphologically, highlighting the utility of DNA barcoding. While most species grouped separately (reciprocally monophyletic), the Cx. pipiens subgroup could not be genetically separated using COI. The average conspecific and congeneric p‐distance was 0.8% and 7.6%, respectively. In our study, we also demonstrate the utility of DNA barcoding in distinguishing exotics from endemic mosquitoes by identifying a single intercepted Stegomyia aegypti egg at an international airport. The use of DNA barcoding dramatically reduced the identification time required compared with rearing specimens through to adults, thereby demonstrating the value of this technique in biosecurity surveillance. The DNA barcodes produced by this study have been uploaded to the ‘Mosquitoes of Australia–Victoria’ project on the Barcode of Life Database (BOLD), which will serve as a resource for the Victorian Arbovirus Disease Control Program and other national and international mosquito surveillance programs.     

Integrative taxonomy of the freshwater worm Rhyacodrilus falciformis s.l. (Clitellata: Naididae), with the description of a new species

The genetic and morphological variation within Rhyacodrilus falciformis Bretscher, 1901 (Clitellata: Naididae) in Europe was explored using an integrative approach, with three unlinked genetic markers [the mitochondrial cytochrome c oxidase subunit I (COI), the nuclear histone 3 (H3) and internal transcribed spacer region (ITS)] combined with morphology, to investigate whether this taxon constitutes a single or several species. Using Automatic Barcode Gap Discovery on the COI data set, the specimens were divided into seven clusters, used as hypothetical species that were further tested with the other data sources. Single‐gene trees were estimated for all three markers, using coalescence analysis and they were in many parts incongruent with each other. Only one of the clusters was supported by all trees; it was also morphologically differentiated from the other clusters by the shape of its modified penial chaetae. This group consists of two specimens from the Crotot Cave in south‐eastern France, and morphologically they fit a previously described but invalid variety, ‘pigueti’, which is here described as a new species, Rhyacodrilus pigueti Achurra & Martinsson sp. n. The study highlights the fact that a single data source (e.g. COI barcodes) seldom provides a sufficient basis for taxonomic decisions such as species delimitation.     

Hybridization relics complicate barcode‐based identification of species in earthworms

Introgressive hybridization results in mito‐nuclear discordance which could obscure the delimitation of closely related taxa. Although such events are increasingly reported, they have been poorly studied in earthworms. Here, we propose a method for investigating the degree of introgressive hybridization between three taxa of the Allolobophora chlorotica aggregate within two field populations (N = 67 and N = 105) using a reference data set including published DNA barcoding and microsatellite data of all known A. chlorotica lineages (N = 85). For this, we used both molecular phylogenetic and population genetic approaches. The test of correspondence between mitochondrial cytochrome c oxidase I (COI) lineages and clusters of nuclear microsatellite genotypes allowed individuals to be sorted in three categories (matching, admixed and nonmatching) and additional markers (mitochondrial NADH dehydrogenase subunit 1, nuclear Histone 3 and Internal transcribed Spacer Region 2) were used for phylogenetic reconstructions in order to check assignments. Although 15 admixed individuals were observed, no early‐generation hybrids were detected within the two populations. Interestingly, 14 nonmatching individuals (i.e. with a mtDNA haplotype that did not correspond to their nuclear cluster) were detected, a pattern that would result after multiple generations of unidirectional hybridization of female from one taxon to male of the other taxon. Because earthworms are simultaneous hermaphrodites, these events of unidirectional hybridization suggest sterility of the male function in several crosses and highlight that some individuals can be misidentified if reliance is placed on COI barcodes alone. These findings could improve the use of these barcodes in earthworms for species delineation.     

DNA barcoding herbaceous and woody plant species at a subalpine forest dynamics plot in Southwest China

Although DNA barcoding has been widely used to identify plant species composition in temperate and tropical ecosystems, relatively few studies have used DNA barcodes to document both herbaceous and woody components of forest plot. A total of 201 species (72 woody species and 129 herbaceous species) representing 135 genera distributed across 64 families of seed plants were collected in a 25 ha CForBio subalpine forest dynamics plot. In total, 491 specimens were screened for three DNA regions of the chloroplast genome (rbcL, matK, and trnH‐psbA) as well as the internal transcribed spacers (ITS) of nuclear ribosomal DNA. We quantified species resolution for each barcode separately or in combination using a ML tree‐based method. Amplification and sequencing success were highest for rbcL, followed by trnH‐psbA, which performed better than ITS and matK. The rbcL + ITS barcode had slightly higher species resolution rates (88.60%) compared with rbcL + matK (86.60%) and rbcL + trnH‐psbA (86.01%). The addition of trnH‐psbA or ITS to the rbcL + matK barcode only marginally increased species resolution rates, although in combination the four barcodes had the highest discriminatory power (90.21%). The situations where DNA barcodes did not discriminate among species were typically associated with higher numbers of co‐occurring con‐generic species. In addition, herbaceous species were much better resolved than woody species. Our study represents one of the first applications of DNA barcodes in a subalpine forest dynamics plot and contributes to our understanding of patterns of genetic divergence among woody and herbaceous plant species.     

DNA barcodes successfully delimit morphospecies in a superdiverse insect genus

Polypedilum Kieffer (Diptera: Chironomidae), with 520 currently known species worldwide, can be extremely difficult to identify species level based on the morphology. We used 3,670 cytochrome c oxidase subunit I (COI) barcodes to explore the efficiency of the COI barcodes to differentiate between species in a superdiverse aquatic insect genus. The Barcode of Life Data System (BOLD) presented 286 BIN clusters in Polypedilum, representing 163 morphospecies, of which 93 were contributed from our laboratory. Molecular operational taxonomic units (OTUs) ranged from 158 to 345, based on Automatic Barcode Gap Discovery (ABGD), the Barcode Index Number (BIN), Bayesian Poisson tree processes (bPTP), generalized mixed Yule coalescent (GMYC), jMOTU, multi‐rate Poisson tree processes (mPTP), neighbor‐joining (NJ) tree and prethreshold clustering. In comparison, GMYC, bPTP, mPTP and BIN suggested more species than warranted by morphology, while ABGD, jMOTU, NJ, prethreshold clustering and ABGD yielded a conservative number of species when setting higher thresholds. Nine species complexes with deep intraspecific divergences indicated 18 potentially cryptic species, which require further taxonomic research including complete life histories and nuclear genetic data to be resolved. The discrimination of Polypedilum species by DNA barcodes proved to be successful in 94.4% of all studied morphological species.     

Cryptic diversity of Melanochlamys sea slugs (Gastropoda,Aglajidae) in the North Pacific

North Pacific specimens of Melanochlamys sea slugs were examined morphologically (including the male reproductive morphology, shell and external coloration) and were sequenced for three genes (mitochondrial COI and 16S and nuclear H3). Phylogenetic and species delimitation analyses were used to determine the evolutionary relationships and species diversity among the specimens examined. Both molecular and morphological data revealed an unexpected high level of cryptic diversity. At least four distinct species occur on the Northwestern Pacific. Melanochlamys ezoensis occurs in Russia and temperate and cold areas in Japan. Three additional undescribed species occur in Japan and/or South Korea. One of the undescribed species occurs both in South Korea and in Japan, but only in Tokyo Bay, suggesting that it could be non‐native in Japan. Two distinct species occur on the Northeastern Pacific coast; Melanochlamys diomedea is widespread from Southern California to Alaska, whereas M. ezoensis was found only in San Francisco Bay, suggesting a human‐mediated introduction. This is further supported by the absence of records of M. ezoensis in San Francisco prior to 2001. The species diversity of Melanochlamys in the Northwestern Pacific is much greater than in the Northeastern Pacific; it is hypothesized that differences in geographic and ocean current system complexity might account for different responses to glacial extinction and postglacial expansion.     

DNA barcoding of Bradysia (Diptera: Sciaridae) for detection of the immature stages on agricultural crops

We investigated the usefulness of mitochondrial cytochrome c oxidase (COI) DNA barcoding of the genus Bradysia for the detection of immature stages and cryptic species complex. Although the larvae of some species in this genus are agricultural pests, immature stages are rarely identified due to the lack of key morphological characteristics. We constructed partial sequences of the COI gene for 25 species of Bradysia as a first step towards a DNA barcode. Using these data, Bradysia impatiens, B. procera and B. peraffinis were identified from larval specimens collected, respectively, from paprika, ginseng and oak sawdust beds used for cultivating shiitake. Our findings reveal a complex of three species within the B. tilicola group. These species were all identified as important pest B. ocellaris based on the morphology of male genital structures; however, the interspecific genetic divergence of the COI region was significantly greater (16.1–19.4%) than the intraspecific variation in each species. Therefore, B. ocellaris may consist of at least three species. The results demonstrate that COI DNA barcodes are useful for Bradysia species identification.     

Effects of prey quality on social wasps when given a choice of prey

The effect of prey quality on foraging behavior and colony demographics of the social wasp Polistes fuscatus was examined by providing a choice between non‐toxic prey (Manduca sexta caterpillars) and sublethally toxic prey (Junonia coenia caterpillars), and then comparing the performance of these colonies to others given only the non‐toxic prey. In the choice, one of two types of Manduca were used: those fed an artificial diet or those fed plantain (Plantago lanceolata), which contains iridoid glycosides (IGs) that Junonia coenia store but which Manduca does not. Despite the negative correlation between the number of Junonia prey used and number of adult offspring produced, when a surplus of non‐toxic prey was available, the wasps did not completely avoid the toxic prey. However, they were more discriminating when the choice was between Manduca fed an artificial diet and Junonia fed plantain vs. when both prey species had eaten the plantain. Because the wasps had a choice of prey types and had a surplus of prey on about one‐third of the days, the wasps were able to take enough non‐toxic prey to avoid some of the negative consequences of IGs. For example, the total number of wasp offspring per nest was not affected, but mean weight of female offspring per colony was less for colonies given both prey types eating plantain, compared to that for colonies fed only non‐toxic prey, or those given a choice of non‐toxic prey vs. toxic prey. In addition, compared to the control (only non‐toxic prey), the proportion of males produced was less in the treatment that provided a clear contrast between non‐toxic and toxic prey. Why these wasps did not avoid the toxic prey is discussed.     

Exploring species boundaries with multiple genetic loci using empirical data from non‐biting midges

Over the past decade, molecular approaches to species delimitation have seen rapid development. However, species delimitation based on a single locus, for example, DNA barcodes, can lead to inaccurate results in cases of recent speciation and incomplete lineage sorting. Here, we compare the performance of Automatic Barcode Gap Discovery (ABGD), Bayesian Poisson tree processes (PTP), networks, generalized mixed Yule coalescent (GMYC) and Bayesian phylogenetics and phylogeography (BPP) models to delineate cryptic species previously detected by DNA barcodes within Tanytarsus (Diptera: Chironomidae) non‐biting midges. We compare the results from analyses of one mitochondrial (cytochrome c oxidase subunit I [COI]) and three nuclear (alanyl‐tRNA synthetase 1 [AATS1], carbamoyl phosphate synthetase 1 [CAD1] and 6‐phosphogluconate dehydrogenase [PGD]) protein‐coding genes. Our results show that species delimitation based on multiple nuclear DNA markers is largely concordant with morphological variation and delimitations using a single locus, for example, the COI barcode. However, ABGD, GMYC, PTP and network models led to conflicting results based on a single locus and delineate species differently than morphology. Results from BPP analyses on multiple loci correspond best with current morphological species concept. In total, 10 lineages of the Tanytarsus curticornis species complex were uncovered. Excluding a Norwegian population of Tanytarsus brundini which might have undergone recent hybridization, this suggests six hitherto unrecognized species new to science. Five distinct species are well supported in the Tanytarsus heusdensis species complex, including two species new to science.