Nuclear genes (but not mitochondrialDNAbarcodes) reveal real species: Evidence from theBrenthisfritillary butterflies (Lepidoptera, Nymphalidae)

Extremely divergent COI sequences within an amphipod species complex: A possible role for endosymbionts?

Some heritable endosymbionts can affect host mtDNA evolution in various ways. Amphipods host diverse endosymbionts, but whether their mtDNA has been influenced by these endosymbionts has yet to be considered. Here, we investigated the role of endosymbionts (microsporidians and Rickettsia) in explaining highly divergent COI sequences in Paracalliope fluviatilis species complex, the most common freshwater amphipods in New Zealand. We first contrasted phylogeographic patterns using COI, ITS, and 28S sequences. While molecular species delimitation methods based on 28S sequences supported 3-4 potential species (N, C, SA, and SB) among freshwater lineages, COI sequences supported 17-27 putative species reflecting high inter-population divergence. The deep divergence between NC and S lineages (~20%; 28S) and the substitution saturation on the 3rd codon position of COI detected even within one lineage (SA) indicate a very high level of morphological stasis. Interestingly, individuals infected and uninfected by Rickettsia comprised divergent COI lineages in one of four populations tested, suggesting a potential influence of endosymbionts in mtDNA patterns. We propose several plausible explanations for divergent COI lineages, although they would need further testing with multiple lines of evidence. Lastly, due to common morphological stasis and the presence of endosymbionts, phylogeographic patterns of amphipods based on mtDNA should be interpreted with caution.

Integrative species delimitation reveals fine-scale allopatric speciation in a good-flying insect: a case study on

Alpha taxonomy is fundamental for many biological fields. Delineation of the species boundary, however, can be challenging in a species complex, where different species share a similar morphology and diagnostic characters may not be available. In this context, integrative approaches that incorporate molecular and morphological data sets, and account for speciation history can be helpful to alpha taxonomy. Different approaches to species delimitation based on different assumptions are complementary and by integrating the results from multiple approaches we can generate a more reliable and objective taxonomic decision. In this study, we applied three molecular approaches to species delimitation and inferred the demographic history based on an isolation with migration model to test a morphologically based taxonomic hypothesis for the Cylindera pseudocylindriformis complex. We discuss the association between genetic divergence and microhabitat specialisation, and further corroborate that C. subtilis sp. nov. is a valid new species by integrating the results from model-based species delimitation and the genealogical divergence index. We argue that genetic endemism can occur at a small geographic scale, even in a winged insect like tiger beetles. Our results also indicated that there may still be undocumented species diversity of Taiwanese Cylindera remaining to be discovered. ZooBank LSID: urn:lsid:zoobank.org:pub:9DEC1432-365C-4872-8D06-73B95F30624F

The genome sequence of the lesser marbled fritillary, Brenthis ino, and evidence for a segregating neo-Z chromosome

The lesser marbled fritillary, Brenthis ino (Rottemburg, 1775), is a species of Palearctic butterfly. Male Brenthis ino individuals have been reported to have between 12 and 14 pairs of chromosomes, a much-reduced chromosome number than is typical in butterflies. Here, we present a chromosome-level genome assembly for Brenthis ino, as well as gene and transposable element annotations. The assembly is 411.8 Mb in length with a contig N50 of 9.6 Mb and a scaffold N50 of 29.5 Mb. We also show evidence that the male individual from which we generated HiC data was heterozygous for a neo-Z chromosome, consistent with inheriting 14 chromosomes from one parent and 13 from the other. This genome assembly will be a valuable resource for studying chromosome evolution in Lepidoptera, as well as for comparative and population genomics more generally.

Chromosomal and DNA barcode analysis of the Melitaea ala Staudinger, 1881 species complex (Lepidoptera, Nymphalidae)

The species of the Melitaea ala Staudinger, 1881 complex are distributed in Central Asia. Here we show that this complex is a monophyletic group including the species, M. ala, M. kotshubeji Sheljuzhko, 1929 and M. enarea Fruhstorfer, 1917. The haploid chromosome number n=29 is found in M. ala and M. kotshubeji and is, most likely, a symplesiomorphy of the M. ala complex. We show that M. ala consists of four subspecies: M. ala zaisana Lukhtanov, 1999 (=M. ala irtyshica Lukhtanov, 1999, syn. nov.) (South Altai, Zaisan Lake valley), M. ala ala (Dzhungarian Alatau), M. ala bicolor Seitz, 1908 (North, East, Central and West Tian-Shan) and M. ala determinata Bryk, 1940 (described from "Fu-Shu-Shi", China). We demonstrate that M. kotshubeji kotshubeji (Peter the Great Mts in Tajikistan) and M. kotshubeji bundeli Kolesnichenko, 1999 (Alai Mts in Tajikistan and Kyrgyzstan) are distinct taxa despite their geographic proximity in East Tajikistan. Melitaea enarea is widely distributed in the southern part of Central Asia and is sympatric with M. kotshubeji.

High resolution DNA barcode library for European butterflies reveals continental patterns of mitochondrial genetic diversity

The study of global biodiversity will greatly benefit from access to comprehensive DNA barcode libraries at continental scale, but such datasets are still very rare. Here, we assemble the first high-resolution reference library for European butterflies that provides 97% taxon coverage (459 species) and 22,306 COI sequences. We estimate that we captured 62% of the total haplotype diversity and show that most species possess a few very common haplotypes and many rare ones. Specimens in the dataset have an average 95.3% probability of being correctly identified. Mitochondrial diversity displayed elevated haplotype richness in southern European refugia, establishing the generality of this key biogeographic pattern for an entire taxonomic group. Fifteen percent of the species are involved in barcode sharing, but two thirds of these cases may reflect the need for further taxonomic research. This dataset provides a unique resource for conservation and for studying evolutionary processes, cryptic species, phylogeography, and ecology.

Overlooked cryptic diversity in Muschampia (Lepidoptera: Hesperiidae) adds two species to the European butterfly fauna

Cryptic species represent a challenge for documenting global biodiversity. Even in well-studied groups, such as European butterflies, the application of integrative approaches has allowed the recognition of an unexpected number of cryptic taxa. Here, we combine the analysis of mitochondrial (cytochrome c oxidase I, COI) and nuclear (internal transcribed spacer 2, ITS2) markers with geometric morphometrics of the male genitalia to study diversity within the butterfly Muschampia proto. The nuclear marker reveals three well-supported and deeply diverged lineages, which are also detected based on mitochondrial DNA, although the latter recovers one of them as paraphyletic with poor support. These lineages also present distinct male genital characters, which allow blind assignment of > 97% of specimens when applying a jackknife procedure. We conclude that M. proto comprises three cryptic species that started to differentiate ~2 Mya: M. proto, distributed in northern Africa, the Iberian Peninsula and southern France; Muschampia alta comb. & stat. nov., occurring in southern Italy and the Balkan Peninsula; and Muschampia proteides, present in the easternmost part of Europe, the Near East and Iran. This discovery adds two new species to the European butterfly fauna and highlights the necessity to continue investigating potential cryptic diversity.

Genetic hitchhiking, mitonuclear coadaptation, and the origins of mt DNA barcode gaps

DNA barcoding based on mitochondrial (mt) nucleotide sequences is an enigma. Neutral models of mt evolution predict DNA barcoding cannot work for recently diverged taxa, and yet, mt DNA barcoding accurately delimits species for many bilaterian animals. Meanwhile, mt DNA barcoding often fails for plants and fungi. I propose that because mt gene products must cofunction with nuclear gene products, the evolution of mt genomes is best understood with full consideration of the two environments that impose selective pressure on mt genes: the external environment and the internal genomic environment. Moreover, it is critical to fully consider the potential for adaptive evolution of not just protein products of mt genes but also of mt transfer RNAs and mt ribosomal RNAs. The tight linkage of genes on mt genomes that do not engage in recombination could facilitate selective sweeps whenever there is positive selection on any element in the mt genome, leading to the purging of mt genetic diversity within a population and to the rapid fixation of novel mt DNA sequences. Accordingly, the most important factor determining whether or not mt DNA sequences diagnose species boundaries may be the extent to which the mt chromosomes engage in recombination.

The conundrum of species delimitation: a genomic perspective on a mitogenetically super-variable butterfly

The Palaearctic butterfly Melitaea didyma stands out as one of the most striking cases of intraspecific genetic differentiation detected in Lepidoptera: 11 partially sympatric mitochondrial lineages have been reported, displaying levels of divergence of up to 7.4%. To better understand the evolutionary processes underlying the diversity observed in mtDNA, we compared mtDNA and genome-wide SNP data using double-digest restriction site-associated DNA sequencing (ddRADseq) results from 93 specimens of M. didyma ranging from Morocco to eastern Kazakhstan. We found that, between ddRADseq and mtDNA results, there is a match only in populations that probably remained allopatric for long periods of time. Other mtDNA lineages may have resulted from introgression events and were probably affected by Wolbachia infection. The five main ddRADseq clades supported by STRUCTURE were parapatric or allopatric and showed high pairwise F_ST values, but some were also estimated to display various levels of gene flow. Melitaea didyma represents one of the first cases of deep mtDNA splits among European butterflies assessed by a genome-wide DNA analysis and reveals that the interpretation of patterns remains challenging even when a high amount of genomic data is available. These findings actualize the ongoing debate of species delimitation in allopatry, an issue probably of relevance to a significant proportion of global biodiversity.

Linking karyotypes with DNA barcodes: proposal for a new standard in chromosomal analysis with an example based on the study of Neotropical Nymphalidae (Lepidoptera)

Chromosomal data are important for taxonomists, cytogeneticists and evolutionary biologists; however, the value of these data decreases sharply if they are obtained for individuals with inaccurate species identification or unclear species identity. To avoid this problem, here we suggest linking each karyotyped sample with its DNA barcode, photograph and precise geographic data, providing an opportunity for unambiguous identification of described taxa and for delimitation of undescribed species. Using this approach, we present new data on chromosome number diversity in neotropical butterflies of the subfamily Biblidinae (genus Vila Kirby, 1871) and the tribe Ithomiini (genera Oleria Hübner, 1816, Ithomia Hübner, 1816, Godyris Boisduval, 1870, Hypothyris Hübner, 1821, Napeogenes Bates, 1862, Pseudoscada Godman et Salvin, 1879 and Hyposcada Godman et Salvin, 1879). Combining new and previously published data we show that the species complex Oleria onega (Hewitson, [1852]) includes three discrete chromosomal clusters (with haploid chromosome numbers n = 15, n = 22 and n = 30) and at least four DNA barcode clusters. Then we discuss how the incomplete connection between these chromosomal and molecular data (karyotypes and DNA barcodes were obtained for different sets of individuals) complicates the taxonomic interpretation of the discovered clusters.