Characterization of the complete mitochondrial genomes of four tarantulas (Arachnida: Theraphosidae) with phylogenetic analysis

To better understand the evolution of mitochondrial genomes (mitogenomes) within the family Theraphosidae, we characterized mitogenomes of four tarantulas (Grammostola pulchripes, Phormictopus atrichomatus, Pterinochilus murinus and Pterinopelma sazimai) for the first time. The mitogenomes were all classical circular structures, with lengths ranging from 13,822 bp to 14,011 bp. The constitutive genes and the orientation of the coding strand observed in the four mitogenomes were consistent with those found in other species belonging to the Theraphosidae family. The four mitogenomes were compacted and exhibited a preference for A and T, with the rRNA sequences showing a higher A+T content. Ka/Ks and p-distances analyses showed the ND6 gene had highest evolutionary rate, while the COⅠ gene displayed relatively slower evolution. In contrast to previous phylogenetic studies, our phylogenetic analysis based on mitogenomes provides new phylogenetic relationships among subfamilies. Subfamily Theraphosinae is most closely related to Ornithoctoninae, slightly distant from Harpactirinae, and farthest from Selenocosmiinae. The new data we acquired regarding these mitogenomes will aid in understanding the complex interrelationships among species within the Theraphosidae family.

Phylogenomics resolves the backbone of Poales and identifies signals of hybridization and polyploidy

Poales, as one of the largest orders of angiosperm, holds crucial economic and ecological importance. Nevertheless, achieving a consensus topology has been challenging in previous studies due to limited molecular data and sparse taxon sampling. The uneven distribution of species diversity among families and the factors leading to elevated species richness in certain lineages have also been subjects of ongoing discussion and investigation. In this study, we conducted a comprehensive sampling, including representatives from all 14 families and 85 taxa of Poales, along with five additional outgroups. To reconstruct the phylogeny of Poales, we employed a combination of coalescent and concatenation methods on three nuclear gene sets (1093, 491, 143) and one plastid gene set (53), which were inferenced from genomic data. We also conducted phylogenetic hypothesis analyses to evaluate two major conflicting nodes detected in phylogenetic analyses. As a result, we successfully resolved the backbone of Poales and provided a timeline for its evolutionary history. We recovered the sister relationship between Typhaceae and Bromeliaceae as the earliest diverging families within Poales. The clade consisting of Ecdeiocoleaceae and Joinvilleaceae was recovered as the sister group of Poaceae. Within the xyrid clade, Mayacaceae and Erioaculaceae + Xyridaceae successively diverged along the backbone of Poales. The topology of [Aristidoideae, ((Micrairoideae, Panicoideae), (Arundinoideae, (Chloridoideae, Danthonioideae)))] within the PACMAD clade has received strong support from multiple findings. We also delved into the underlying biological factors that contributed to the conflicting nodes observed in the phylogenetic analysis. Apart from the uncertainty regarding the sister group of Poaceae caused by cytonuclear discordance, frequent hybridization and polyploidy may have contributed to other conflicting nodes. We identified 26 putative whole-genome duplication (WGD) events within Poales. However, apart from the σ-WGD and the ρ-WGD, we did not observe any potential polyploid events that could be directly linked to the species diversification in specific lineages. Furthermore, there was a significant increase in the net diversification rate of Poales following the K-Pg boundary.

Intraspecific diversity of Myrmecophilus acervorum (Orthoptera: Myrmecophilidae) indicating an ongoing cryptic speciation

Myrmecophilus acervorum, previously considered a parthenogenetic species widely-distributed in Europe, has been observed to have both sexes in populations inhabiting the central part of the distribution range. Specimens from those heterosexual populations have been found being infected with Wolbachia. New mitochondrial data (COI and 16S markers) revealed the well-supported differentiation of M. acervorum populations inhabiting western Polesie (Poland) and southern Europe. In turn, analyses of EF1α marker support the hypothesis on the unfinished lineage sorting at the nuclear DNA level. Interestingly, we found that parthenogenetic populations inhabiting western Polesie are infected with Wolbachia belonging to supergroup A, while endosymbionts occurring in sexual populations of M. acervorum observed in Romania belong to supergroup B. Furthermore, new and potentially diagnostic characteristics in the external structures of the eyes of M. acervorum were identified. The surface of ommatidia in specimens occurring in southern Europe was smooth. In contrast, the ommatidia surface of individuals collected in Poland was visibly sculptured. To sum up, the significant genetic variability found in the present case, and the differentiating morphological character, are almost certainly effects of cryptic species being present within M. acervorum. This is indicative of ongoing speciation within the populations of this insect, and of simultaneous unfinished lineage sorting at the nuclear DNA level.

Phylogenetics and species delimitation of the recluse spider, Loxosceles rufescens (Araneae: Sicariidae) populations invading Bangkok, Thailand

The Mediterranean recluse spider, Loxosceles rufescens, has been discovered for the first time inhabiting human dwellings in Bangkok, Thailand. Expeditions across 39 localities revealed five establishments with L. rufescens populations. The highest density was recorded in a storage house on Yaowarat Road, located in the heart of Bangkok's Chinatown, where 315 individuals were found, including adults, juveniles, and spiderlings. This medically significant spider's presence in such a densely populated urban area raises concerns about potential envenomation risks. Thirteen specimens of L. rufescens were extracted for DNA and sequenced for molecular phylogenetic analyses. COI and ITS2 markers were used to investigate relationships within L. rufescens and across available Loxosceles species sequences. Results indicate COI is superior for resolving species-level genetic clusters compared to ITS2. Surprisingly, L. rufescens individuals from the same house were found in significantly distant COI lineages, suggesting mtDNA may not be suitable for studying intra-specific phylogeography in this case. Species delimitation methods ABGD and ASAP demonstrated promising results for both COI and ITS2, while bPTP and GMYC tended to overestimate species numbers. ITS2 exhibited high sequence similarity in L. rufescens, suggesting potential utility as a barcoding marker for identification of this globally distributed species. Genetic distance analyses revealed a potential barcoding gap (K2P) of 8-9 % for COI and <2 % for ITS2 in Loxosceles. This study contributes valuable sequence data for the medically important genus Loxosceles and highlights the need for integrative approaches in understanding its evolution and spread. The findings have important implications for pest management strategies and public health in urban environments.

The mitochondrial genome of Lavandula angustifolia Mill. (Lamiaceae) sheds light on its genome structure and gene transfer between organelles

BACKGROUND

Lavandula angustifolia holds importance as an aromatic plant with extensive applications spanning the fragrance, perfume, cosmetics, aromatherapy, and spa sectors. Beyond its aesthetic and sensory applications, this plant offers medicinal benefits as a natural herbal remedy and finds use in household cleaning products. While extensive genomic data, inclusive of plastid and nuclear genomes, are available for this species, researchers have yet to characterize its mitochondrial genome. This gap in knowledge hampers deeper understanding of the genome organization and its evolutionary significance.

RESULTS

Through the course of this study, we successfully assembled and annotated the mitochondrial genome of L. angustifolia, marking a first in this domain. This assembled genome encompasses 61 genes, which comprise 34 protein-coding genes, 24 transfer RNA genes, and three ribosomal RNA genes. We identified a chloroplast sequence insertion into the mitogenome, which spans a length of 10,645 bp, accounting for 2.94% of the mitogenome size. Within these inserted sequences, there are seven intact tRNA genes (trnH-GUG, trnW-CCA, trnD-GUC, trnS-GGA, trnN-GUU, trnT-GGU, trnP-UGG) and four complete protein-coding genes (psbA, rps15, petL, petG) of chloroplast derivation. Additional discoveries include 88 microsatellites, 15 tandem repeats, 74 palindromic repeats, and 87 forward long repeats. An RNA editing analysis highlighted an elevated count of editing sites in the cytochrome c oxidase genes, notably ccmB with 34 editing sites, ccmFN with 32, and ccmC with 29. All protein-coding genes showed evidence of cytidine-to-uracil conversion. A phylogenetic analysis, utilizing common protein-coding genes from 23 Lamiales species, yielded a tree with consistent topology, supported by high confidence values.

CONCLUSIONS

Analysis of the current mitogenome resource revealed its typical circular genome structure. Notably, sequences originally from the chloroplast genome were found within the mitogenome, pointing to the occurrence of horizontal gene transfer between organelles. This assembled mitogenome stands as a valuable resource for subsequent studies on mitogenome structures, their evolution, and molecular biology.

Notes from the taxonomic disaster zone: Evolutionary drivers of intractable species boundaries in an Australian lizard clade (Scincidae: Ctenotus)

Genomic-scale datasets, sophisticated analytical techniques, and conceptual advances have disproportionately failed to resolve species boundaries in some groups relative to others. To understand the processes that underlie taxonomic intractability, we dissect the speciation history of an Australian lizard clade that arguably represents a "worst-case" scenario for species delimitation within vertebrates: the Ctenotus inornatus species group, a clade beset with decoupled genetic and phenotypic breaks, uncertain geographic ranges, and parallelism in purportedly diagnostic morphological characters. We sampled hundreds of localities to generate a genomic perspective on population divergence, structure, and admixture. Our results revealed rampant paraphyly of nominate taxa in the group, with lineages that are either morphologically cryptic or polytypic. Isolation-by-distance patterns reflect spatially continuous differentiation among certain pairs of putative species, yet genetic and geographic distances are decoupled in other pairs. Comparisons of mitochondrial and nuclear gene trees, tests of nuclear introgression, and historical demographic modelling identified gene flow between divergent candidate species. Levels of admixture are decoupled from phylogenetic relatedness; gene flow is often higher between sympatric species than between parapatric populations of the same species. Such idiosyncratic patterns of introgression contribute to species boundaries that are fuzzy while also varying in fuzziness. Our results suggest that "taxonomic disaster zones" like the C. inornatus species group result from spatial variation in the porosity of species boundaries and the resulting patterns of genetic and phenotypic variation. This study raises questions about the origin and persistence of hybridizing species and highlights the unique insights provided by taxa that have long eluded straightforward taxonomic categorization.

Morphological description and genetic analysis of a new black fly species (Diptera: Simuliidae) in the subgenus Asiosimulium from central Thailand

BACKGROUND

Black flies are among the most medically and veterinary important insects, as adult females of certain species are the sole vector of Onchocerca volvulus. Here, a new black fly species belonging to the subgenus Asiosimulium Takaoka & Choochote, 2005, is described and formally named as Simulium (Asiosimulium) kittipati sp. nov.

METHODS

Pupae and larvae of black flies were collected from available substrates in the stream from central Thailand. Pupae were individually separated in plastic tubes and maintained until adult flies emerged. The emerged adult flies associated with their pupal exuviae and cocoon as well as mature larvae preserved in 85% ethanol were used to describe the new species based on an integrated approach of morphological examination and molecular analysis of the COI gene.

RESULTS

The new species is characterized in the female by the medium-long sensory vesicle with a medium-sized opening apically, scutum with three faint longitudinal vittae, and the ellipsoidal spermatheca; in the male by the number of upper-eye (large) facets in 20 vertical columns and 21 horizontal rows, hind basitarsus slender, nearly parallel-sided, and median sclerite much wider and upturned apically; in the pupa by the head and thoracic integument densely covered with tiny tubercles, and the pupal gill of arborescent type with 28-30 filaments; and in the larva by the postgenal cleft deep, nearly reaching the posterior margin of the hypostoma, and dark pigmented sheath of the subesophageal ganglion. The DNA barcode successfully differentiated the new species from its congeners with an interspecific genetic divergence of 1.74-18.72%, confirming the morphological identification that the species is a new member of the subgenus Asiosimulium. Phylogenetic analyses also indicated that the new species is genetically closely related to Simulium phurueaense Tangkawanit, Wongpakam & Pramual, 2018, further supporting its morphological classification.

CONCLUSIONS

This is the ninth species assigned to the subgenus Asiosimulium within the genus Simulium Latreille, 1802. Taxonomic notes and identification keys are given to distinguish this new species from the eight known species members in its same subgenus. Additionally, a distribution map of all species members in this subgenus occurring in Thailand and other countries is provided.

Mitonuclear and phenotypic discordance in an Atlantic Forest frog hybrid zone

Discordance between mitochondrial and nuclear DNA is common among animals and can be the result of a number of evolutionary processes, including incomplete lineage sorting and introgression. Particularly relevant in contact zones, mitonuclear discordance is expected because the mitochondrial genome is haploid and primarily uniparentally inherited, whereas nuclear loci are evolving at slower rates. In addition, when closely related taxa come together in hybrid zones, the distribution of diagnostic phenotypic characters and their concordance with the mitochondrial or nuclear lineages can also inform on historical and ongoing dynamics within hybrid zones. Overall, genetic and phenotypic discordances provide evidence for evolutionary divergence and processes that maintain boundaries among sister species or lineages. In this study, we characterized patterns of genetic and phenotypic variation in a contact zone between Cycloramphus dubius and Cycloramphus boraceiensis, two sister species of frogs endemic to the Atlantic Coastal Forest of Brazil. We examined genomic-scale nuclear diversification across 19 populations, encompassing the two parental forms and a contact zone between them. We compared the distribution of genomic DNA variability with that of a mitochondrial locus (16S) and two morphological traits (dorsal tubercles and body size). Our results reveal multiple divergent lineages with ongoing admixture. We detected discordance in patterns of introgression across the three data types. Cycloramphus dubius males are significantly larger than C. boraceiensis males, and we posit that competition among males in the hybrid zone, coupled with mate choice by females, may be one mechanism leading to patterns of introgression observed between the species.

Discovery of a new tarantula species from the Madrean Sky Islands and the first documented instance of syntopy between two montane endemics (Araneae, Theraphosidae, Aphonopelma): a case of prior mistaken identity

The Chiricahua Mountains in southeastern Arizona are renowned for their exceptional biodiversity and high levels of endemism. Morphological, genomic, behavioral, and distributional data were used to report the discovery of a remarkable new tarantula species from this range. Aphonopelmajacobii sp. nov. inhabits high-elevation mixed conifer forests in these mountains, but also co-occurs and shares its breeding period with A.chiricahua-a related member of the Marxi species group-in mid-elevation Madrean evergreen oak and pine-oak woodlands. This marks the first documented case of syntopy between two montane endemics in the Madrean Archipelago and adds to our knowledge of this threatened region's unmatched tarantula diversity in the United States. An emended diagnosis and redescription for A.chiricahua are also provided based on several newly acquired and accurately identified specimens. Phylogenetic analyses of mitochondrial and genomic-scale data reveal that A.jacobii sp. nov. is more closely related to A.marxi, a species primarily distributed on the Colorado Plateau, than to A.chiricahua or the other Madrean Sky Island taxa. These data provide the evolutionary framework for better understanding the region's complex biogeographic history (e.g., biotic assembly of the Chiricahua Mountains) and conservation of these spiders.

Accelerated mitochondrial evolution and asymmetric fitness of hybrids contribute to the persistence of Helix thessalica in the Helix pomatia range

Interbreeding and introgression between recently diverged species is common. However, the processes that prevent these species from merging where they co-occur are not well understood. We studied the mechanisms that allowed an isolated group of populations of the snail Helix thessalica to persist within the range of the related Helix pomatia despite high gene flow. Using genomic cline analysis, we found that the nuclear gene flow between the two taxa across the mosaic hybrid zone was not different from that expected under neutral admixture, but that the exchange of mtDNA was asymmetric. Tests showed that there is relaxed selection in the mitochondrial genome of H. thessalica and that the substitution rate is elevated compared to that of H. pomatia. A lack of hybrids that combine the mtDNA of H. thessalica with a mainly (>46%) H. pomatia genomic background indicates that the nuclear-encoded mitochondrial proteins of H. pomatia are not well adapted to the more rapidly evolving proteins and RNAs encoded by the mitochondrion of H. thessalica. The presumed reduction of fitness of hybrids with the fast-evolving mtDNA of H. thessalica and a high H. pomatia ancestry, similar to 'Darwin's Corollary to Haldane's rule', resulted in a relative loss of H. pomatia nuclear ancestry compared to H. thessalica ancestry in the hybrid zone. This probably prevents the H. thessalica populations from merging quickly with the surrounding H. pomatia populations and supports the hypothesis that incompatibilities between rapidly evolving mitochondrial genes and nuclear genes contribute to speciation.

Towards a better future for DNA barcoding: Evaluating monophyly- and distance-based species identification using COI gene fragments of Dacini fruit flies

The utility of a universal DNA 'barcode' fragment (658 base pairs of the Cytochrome C Oxidase I [COI] gene) has been established as a useful tool for species identification, and widely criticized as one for understanding the evolutionary history of a group. Large amounts of COI sequence data have been produced that hold promise for rapid species identification, for example, for biosecurity. The fruit fly tribe Dacini holds about a thousand species, of which 80 are pests of economic concern. We generated a COI reference library for 265 species of Dacini containing 5601 sequences that span most of the COI gene using circular consensus sequencing. We compared distance metrics versus monophyly assessments for species identification and although we found a 'soft' barcode gap around 2% pairwise distance, the exceptions to this rule dictate that a monophyly assessment is the only reliable method for species identification. We found that all fragments regularly used for Dacini fruit fly identification >450 base pairs long provide similar resolution. 11.3% of the species in our dataset were non-monophyletic in a COI tree, which is mostly due to species complexes. We conclude with recommendations for the future generation and use of COI libraries. We revise the generic assignment of Dacus transversus stat. rev. Hardy 1982, and Dacus perpusillus stat. rev. Drew 1971 and we establish Dacus maculipterus White 1998 syn. nov. as a junior synonym of Dacus satanas Liang et al. 1993.

Comparative Mitogenomics and Phylogenetic Implications for Nine Species of the Subfamily Meconematinae (Orthoptera: Tettigoniidae)

Currently, the subfamily Meconematinae encompasses 1029 species, but whole-mitochondrial-genome assemblies have only been made available for 13. In this study, the whole mitochondrial genomes (mitogenomes) of nine additional species in the subfamily Meconematinae were sequenced. The size ranged from 15,627 bp to 17,461 bp, indicating double-stranded circular structures. The length of the control region was the main cause of the difference in mitochondrial genome length among the nine species. All the mitogenomes including 13 protein-coding genes (PCGs), 22 transfer RNA genes (tRNAs), 2 ribosomal RNA genes (rRNAs) and a control region (CR). The majority strand encoded 23 genes, and the minority strand encoded 14 genes. A phylogenetic analysis reaffirmed the monophyletic status of each subfamily, but the monophysitism of Xizicus, Xiphidiopsis and Phlugiolopsis was not supported.

Integrating mitochondrial and nuclear genomic data to decipher the evolutionary history of Eubranchipus species in Japan

Understanding the genetic diversity and evolutionary history of species is crucial for their conservation and management. In this study, we investigated the genetic diversity and phylogenetic relationships among Eubranchipus species occurring in Japan. Phylogenetic analyses revealed that nuclear and mitochondrial data yield incompatible results. In E. uchidai, nuclear data support the monophyly of the Shimokita area, while mitochondrial data indicate a clustering of Higashidori2 individuals with Hokkaido (Ishikari and Wakkanai) E. uchidai. Similar incongruences were observed in E. hatanakai, where nuclear data favor the monophyly of the Chokai area, while mitochondrial data cluster some Chokai pool 3 individuals with Aizu individuals. These incompatibilities might be caused by mitochondrial gene flow. The findings emphasize the importance of considering both nuclear and mitochondrial data during phylogenetic studies and provide valuable insights into the complex dynamics of migration and genetic exchange in Eubranchipus species.

Wolbachia infection status and molecular diversity in the species of tribe Tagiadini Mabille, 1878 (Lepidoptera: Hesperiidae) collected in China

Wolbachia, one of the most ubiquitous heritable symbionts in lepidopteran insects, can cause mitochondrial introgression in related host species. We recently found mito-nuclear discordance in the Lepidopteran tribe Tagiadini Mabille 1878 from which Wolbachia has not been reported. In this study, we found that 13 of the 46 species of Tagiadini species tested were positive for Wolbachia. Overall, 14% (15/110) of Tagiadini specimens were infected with Wolbachia and nine new STs were found from 15 isolates. A co-phylogenetic comparison, divergence time estimation and Wolbachia recombination analysis revealed that mito-nuclear discordance in Tagiadini species is not mediated by Wolbachia, but Wolbachia acquisition in Tagiadini appears to have occurred mainly through horizontal transmission rather than codivergence.

Metazoa-level USCOs as markers in species delimitation and classification

Metazoa-level universal single-copy orthologs (mzl-USCOs) are universally applicable markers for DNA taxonomy in animals that can replace or supplement single-gene barcodes. Previously, mzl-USCOs from target enrichment data were shown to reliably distinguish species. Here, we tested whether USCOs are an evenly distributed, representative sample of a given metazoan genome and therefore able to cope with past hybridization events and incomplete lineage sorting. This is relevant for coalescent-based species delimitation approaches, which critically depend on the assumption that the investigated loci do not exhibit autocorrelation due to physical linkage. Based on 239 chromosome-level assembled genomes, we confirmed that mzl-USCOs are genetically unlinked for practical purposes and a representative sample of a genome in terms of reciprocal distances between USCOs on a chromosome and of distribution across chromosomes. We tested the suitability of mzl-USCOs extracted from genomes for species delimitation and phylogeny in four case studies: Anopheles mosquitos, Drosophila fruit flies, Heliconius butterflies and Darwin's finches. In almost all instances, USCOs allowed delineating species and yielded phylogenies that corresponded to those generated from whole genome data. Our phylogenetic analyses demonstrate that USCOs may complement single-gene DNA barcodes and provide more accurate taxonomic inferences. Combining USCOs from sources that used different versions of ortholog reference libraries to infer marker orthology may be challenging and, at times, impact taxonomic conclusions. However, we expect this problem to become less severe as the rapidly growing number of reference genomes provides a better representation of the number and diversity of organismal lineages.

Species delimitation, discovery and conservation in a tiger beetle species complex despite discordant genetic data

In an age of species declines, delineating and discovering biodiversity is critical for both taxonomic accuracy and conservation. In recent years, there has been a movement away from using exclusively morphological characters to delineate and describe taxa and an increase in the use of molecular markers to describe diversity or through integrative taxonomy, which employs traditional morphological characters, as well as genetic or other data. Tiger beetles are charismatic, of conservation concern, and much work has been done on the morphological delineation of species and subspecies, but few of these taxa have been tested with genetic analyses. In this study, we tested morphologically based taxonomic hypotheses of polymorphic tiger beetles in the Eunota circumpicta (LaFerté-Sénectère, 1841) species complex using multilocus genomic and mtDNA analyses. We find multiple cryptic species within the previous taxonomic concept of Eunota circumpicta, some of which were historically recognized as subspecies. We found that the mtDNA and genomic datasets did not identify the same taxonomic units and that the mtDNA was most at odds with all other genetic and morphological patterns. Overall, we describe new cryptic diversity, which raises important conservation concerns, and provide a working example for testing species and subspecies validity despite discordant data.

A comprehensive DNA barcoding reference database for Plecoptera of Switzerland

DNA barcoding is an essential tool in modern biodiversity sciences. Despite considerable work to barcode the tree of life, many groups, including insects, remain partially or totally unreferenced, preventing barcoding from reaching its full potential. Aquatic insects, especially the three orders Ephemeroptera, Plecoptera, and Trichoptera (EPT), are key freshwater quality indicators worldwide. Among them, Plecoptera (stoneflies), which are among the most sensitive aquatic insects to habitat modification, play a central role in river monitoring surveys. Here, we present an update of the Plecoptera reference database for (meta)barcoding in Switzerland, now covering all 118 species known from this country. Fresh specimens, mostly from rare or localized species, were collected, and 151 new CO1 barcodes were generated. These were merged with the 422 previously published sequences, resulting in a dataset of 573 barcoded specimens. Our CO1 dataset was delimited in 115 CO1 clusters based on a priori morphological identifications, of which 17% are newly reported for Switzerland, and 4% are newly reported globally. Among the 115 CO1 clusters, 85% showed complete congruence with morphology. Distance-based analysis indicated local barcoding gaps in 97% of the CO1 clusters. This study significantly improves the Swiss reference database for stoneflies, enhancing future species identification accuracy and biodiversity monitoring. Additionally, this work reveals cryptic diversity and incongruence between morphology and barcodes, both presenting valuable opportunities for future integrative taxonomic studies. Voucher specimens, DNA extractions and reference barcodes are available for future developments, including metabarcoding and environmental DNA surveys.

Species delimitation of tea plants (Camellia sect. Thea) based on super-barcodes

BACKGROUND

The era of high throughput sequencing offers new paths to identifying species boundaries that are complementary to traditional morphology-based delimitations. De novo species delimitation using traditional or DNA super-barcodes serve as efficient approaches to recognizing putative species (molecular operational taxonomic units, MOTUs). Tea plants (Camellia sect. Thea) form a group of morphologically similar species with significant economic value, providing the raw material for tea, which is the most popular nonalcoholic caffeine-containing beverage in the world. Taxonomic challenges have arisen from vague species boundaries in this group.

RESULTS

Based on the most comprehensive sampling of C. sect. Thea by far (165 individuals of 39 morphospecies), we applied three de novo species delimitation methods (ASAP, PTP, and mPTP) using plastome data to provide an independent evaluation of morphology-based species boundaries in tea plants. Comparing MOTU partitions with morphospecies, we particularly tested the congruence of MOTUs resulting from different methods. We recognized 28 consensus MOTUs within C. sect. Thea, while tentatively suggesting that 11 morphospecies be discarded. Ten of the 28 consensus MOTUs were uncovered as morphospecies complexes in need of further study integrating other evidence. Our results also showed a strong imbalance among the analyzed MOTUs in terms of the number of molecular diagnostic characters.

CONCLUSION

This study serves as a solid step forward for recognizing the underlying species boundaries of tea plants, providing a needed evidence-based framework for the utilization and conservation of this economically important plant group.

Facilitating taxonomy and phylogenetics: An informative and cost-effective protocol integrating long amplicon PCRs and third-generation sequencing

Phylogenetic inference has become a standard technique in integrative taxonomy and systematics, as well as in biogeography and ecology. DNA barcodes are often used for phylogenetic inference, despite being strongly limited due to their low number of informative sites. Also, because current DNA barcodes are based on a fraction of a single, fast-evolving gene, they are highly unsuitable for resolving deeper phylogenetic relationships due to saturation. In recent years, methods that analyse hundreds and thousands of loci at once have improved the resolution of the Tree of Life, but these methods require resources, experience and molecular laboratories that most taxonomists do not have. This paper introduces a PCR-based protocol that produces long amplicons of both slow- and fast-evolving unlinked mitochondrial and nuclear gene regions, which can be sequenced by the affordable and portable ONT MinION platform with low infrastructure or funding requirements. As a proof of concept, we inferred a phylogeny of a sample of 63 spider species from 20 families using our proposed protocol. The results were overall consistent with the results from approaches based on hundreds and thousands of loci, while requiring just a fraction of the cost and labour of such approaches, making our protocol accessible to taxonomists worldwide.

Landscape Heterogeneity Explains the Genetic Differentiation of a Forest Bird across the Sino-Himalayan Mountains

Mountains are the world's most important centers of biodiversity. The Sino-Himalayan Mountains are global biodiversity hotspot due to their extremely high species richness and endemicity. Ample research investigated the impact of the Qinghai-Tibet Plateau uplift and Quaternary glaciations in driving species diversification in plants and animals across the Sino-Himalayan Mountains. However, little is known about the role of landscape heterogeneity and other environmental features in driving diversification in this region. We utilized whole genomes and phenotypic data in combination with landscape genetic approaches to investigate population structure, demography, and genetic diversity in a forest songbird species native to the Sino-Himalayan Mountains, the red-billed leiothrix (Leiothrix lutea). We identified 5 phylogeographic clades, including 1 in the East of China, 1 in Yunnan, and 3 in Tibet, roughly consistent with differences in song and plumage coloration but incongruent with traditional subspecies boundaries. Isolation-by-resistance model best explained population differentiation within L. lutea, with extensive secondary contact after allopatric isolation leading to admixture among clades. Ecological niche modeling indicated relative stability in the extent of suitable distribution areas of the species across Quaternary glacial cycles. Our results underscore the importance of mountains in the diversification of this species, given that most of the distinct genetic clades are concentrated in a relatively small area in the Sino-Himalayan Mountain region, while a single shallow clade populates vast lower-lying areas to the east. This study highlights the crucial role of landscape heterogeneity in promoting differentiation and provides a deep genomic perspective on the mechanisms through which diversity hotspots form.

Microgeographic population structuring in a genus of California trapdoor spiders and discovery of an enigmatic new species (Euctenizidae: Promyrmekiaphila korematsui sp. nov.)

The recognition and delineation of cryptic species remains a perplexing problem in systematics, evolution, and species delimitation. Once recognized as such, cryptic species complexes provide fertile ground for studying genetic divergence within the context of phenotypic and ecological divergence (or lack thereof). Herein we document the discovery of a new cryptic species of trapdoor spider, Promyrmekiaphila korematsui sp. nov. Using subgenomic data obtained via target enrichment, we document the phylogeography of the California endemic genus Promyrmekiaphila and its constituent species, which also includes P. clathrata and P. winnemem. Based on these data we show a pattern of strong geographic structuring among populations but cannot entirely discount recent gene flow among populations that are parapatric, particularly for deeply diverged lineages within P. clathrata. The genetic data, in addition to revealing a new undescribed species, also allude to a pattern of potential phenotypic differentiation where species likely come into close contact. Alternatively, phenotypic cohesion among genetically divergent P. clathrata lineages suggests that some level of gene flow is ongoing or occurred in the recent past. Despite considerable field collection efforts over many years, additional sampling in potential zones of contact for both species and lineages is needed to completely resolve the dynamics of divergence in Promyrmekiaphila at the population-species interface.

Pseudocryptic diversity and species boundaries in the sea cucumber Stichopus cf. horrens (Echinodermata: Stichopodidae) revealed by mitochondrial and microsatellite markers

Morphologically cryptic and pseudo-cryptic species pose a challenge to taxonomic identification and assessments of species diversity and distributions. Such is the case for the sea cucumber Stichopus horrens, commonly confused with Stichopus monotuberculatus. Here, we used mitochondrial cytochrome oxidase subunit I (COI) and microsatellite markers to examine genetic diversity in Stichopus cf. horrens throughout the Philippine archipelago, to aid species identification and clarify species boundaries. Phylogenetic analysis reveals two recently diverged COI lineages (Clade A and Clade B; c. 1.35-2.54 Mya) corresponding to sequence records for specimens identified as S. monotuberculatus and S. horrens, respectively. Microsatellite markers reveal two significantly differentiated genotype clusters broadly concordant with COI lineages (Cluster 1, Cluster 2). A small proportion of individuals were identified as later-generation hybrids indicating limited contemporary gene flow between genotype clusters, thus confirming species boundaries. Morphological differences in papillae distribution and form are observed for the two species, however tack-like spicules from the dorsal papillae are not a reliable diagnostic character. An additional putative cryptic species was detected within Clade B-Cluster 2 specimens warranting further examination. We propose that these lineages revealed by COI and genotype data be referred to as Stichopus cf. horrens species complex.

Ancient mitochondrial genomes unveil the origins and evolutionary history of New Zealand's enigmatic takahē and moho

Many avian species endemic to Aotearoa New Zealand were driven to extinction or reduced to relict populations following successive waves of human arrival, due to hunting, habitat destruction and the introduction of mammalian predators. Among the affected species were the large flightless South Island takahē (Porphyrio hochstetteri) and the moho (North Island takahē; P. mantelli), with the latter rendered extinct and the former reduced to a single relictual population. Little is known about the evolutionary history of these species prior to their decline and/or extinction. Here we sequenced mitochondrial genomes from takahē and moho subfossils (12 takahē and 4 moho) and retrieved comparable sequence data from takahē museum skins (n = 5) and contemporary individuals (n = 17) to examine the phylogeny and recent evolutionary history of these species. Our analyses suggest that prehistoric takahē populations lacked deep phylogeographic structure, in contrast to moho, which exhibited significant spatial genetic structure, albeit based on limited sample sizes (n = 4). Temporal genetic comparisons show that takahē have lost much of their mitochondrial genetic diversity, likely due to a sudden demographic decline soon after human arrival (~750 years ago). Time-calibrated phylogenetic analyses strongly support a sister species relationship between takahē and moho, suggesting these flightless taxa diverged around 1.5 million years ago, following a single colonisation of New Zealand by a flighted Porphyrio ancestor approximately 4 million years ago. This study highlights the utility of palaeogenetic approaches for informing the conservation and systematic understanding of endangered species whose ranges have been severely restricted by anthropogenic impacts.

Cryptic diversity begets challenges and opportunities in biodiversity research

How many species of life are there on Earth? This is a question that we want to know but cannot yet answer. Some scholars speculate that the number of species may reach 2.2 billion when considering cryptic diversity and that each morphology-based insect species may contain an average of 3.1 cryptic species. With nearly two million described species, such high estimates of cryptic diversity would suggest that cryptic species are widespread. The development of molecular species delimitation has led to the discovery of a large number of cryptic species, and cryptic biodiversity has gradually entered our field of vision and attracted more attention. This paper introduces the concept of cryptic species, how they evolve, and methods by which they may be discovered and confirmed, and provides theoretical and methodological guidance for the study of hidden species. A workflow of how to confirm cryptic species is provided. In addition, the importance and reliability of multi-evidence-based integrated taxonomy are reaffirmed as a way to better standardize decision-making processes. Special focus on cryptic diversity and increased funding for taxonomy is needed to ensure that cryptic species in hyperdiverse groups are discoverable and described. An increased focus on cryptic species in the future will naturally arise as more difficult groups are studied, and thereby, we may finally better understand the rules governing the evolution and maintenance of cryptic biodiversity.

Development of molecular identification methods for Dryophytes suweonensis and D. japonicus, and their hybrids

Background

As hybridization can reduce biodiversity or cause extinction, it is important to identify both purebred parental species and their hybrids prior to conserving them. The Suwon tree frog, Dryophytes suweonensis, is an endangered wildlife species in Korea that shares its habitat and often hybridizes with the Japanese tree frog, D. japonicus. In particular, D. suweonensis, D. japonicus, and their hybrids often have abnormal ovaries and gonads, which are known causes that could threaten their existence.

Methods

We collected 57 individuals from six localities where D. suweonensis is known to be present. High-resolution melting curve (HRM) analysis of the mitochondrial 12S ribosomal RNA gene was performed to determine the maternal species. Thereafter, the DNA sequences of five nuclear genes (SIAH, TYR, POMC, RAG1, and C-MYC) were analyzed to determine their parental species and hybrid status.

Results

The HRM analysis showed that the melting temperature of D. suweonensis was in the range of 79.0-79.3 °C, and that of D. japonicus was 77.7-78.0 °C, which clearly distinguished the two tree frog species. DNA sequencing of the five nuclear genes revealed 37 single-nucleotide polymorphism (SNP) sites, and STRUCTURE analysis showed a two-group structure as the most likely grouping solution. No heterozygous position in the purebred parental sequences with Q values ≥ 0.995 were found, which clearly distinguished the two treefrog species from their hybrids; 11 individuals were found to be D. suweonensis, eight were found to be D. japonicus, and the remaining 38 individuals were found to be hybrids.

Conclusion

Thus, it was possible to unambiguously identify the parental species and their hybrids using HRM analysis and DNA sequencing methods. This study provided fundamental information for D. suweonensis conservation and restoration research.

DNA Barcoding in Species Delimitation: From Genetic Distances to Integrative Taxonomy

Over the past two decades, DNA barcoding has become the most popular exploration approach in molecular taxonomy, whether for identification, discovery, delimitation, or description of species. The present contribution focuses on the utility of DNA barcoding for taxonomic research activities related to species delimitation, emphasizing the following aspects:(1) To what extent DNA barcoding can be a valuable ally for fundamental taxonomic research, (2) its methodological and theoretical limitations, (3) the conceptual background and practical use of pairwise distances between DNA barcode sequences in taxonomy, and (4) the different ways in which DNA barcoding can be combined with complementary means of investigation within a broader integrative framework. In this chapter, we recall and discuss the key conceptual advances that have led to the so-called renaissance of taxonomy, elaborate a detailed glossary for the terms specific to this discipline (see Glossary in Chap. 35 ), and propose a newly designed step-by-step species delimitation protocol starting from DNA barcode data that includes steps from the preliminary elaboration of an optimal sampling strategy to the final decision-making process which potentially leads to nomenclatural changes.

Species Diagnosis and DNA Taxonomy

The use of DNA has helped to improve and speed up species identification and delimitation. However, it also provides new challenges to taxonomists. Incongruence of outcome from various markers and delimitation methods, bias from sampling and skewed species distribution, implemented models, and the choice of methods/priors may mislead results and also may, in conclusion, increase elements of subjectivity in species taxonomy. The lack of direct diagnostic outcome from most contemporary molecular delimitation approaches and the need for a reference to existing and best sampled trait reference systems reveal the need for refining the criteria of species diagnosis and diagnosability in the current framework of nomenclature codes and good practices to avoid nomenclatorial instability, parallel taxonomies, and consequently more and new taxonomic impediment.

iTaxoTools 1.0: Improved DNA Barcode Exploration with TaxI2

The overall availability of user-friendly software tools tailored to the analysis of DNA barcodes is limited. Several obvious functions such as detecting and visualizing the DNA barcode gap, the calculation of matrices of pairwise distances at the level of species, or the filtering and decontaminating of sets of sequences based on comparisons with reference databases can typically be carried out only by complex procedures that involve various programs and/or a substantial manual work of formatting. The iTaxoTools project aims at contributing user-friendly software solutions to improve the speed and quality of the workflow of alpha-taxonomy. In this chapter, we provide detailed protocols for the use of a substantially improved version of the tool TaxI2 for distance-based exploration of DNA barcodes. The program calculates genetic distances from prealigned data sets, or based on pairwise alignments, or with an alignment-free approach. Sequence and metadata input can be formatted as tab-delimited files and TaxI2 then computes tables, matrices and graphs of distances, and distance summary statistics within and between species and genera. TaxI2 also includes modes to compare a set of sequences against one or two reference data sets and output lists of best matches or filter data according to thresholds or reciprocal matches. Here, detailed step-by-step protocols are provided for the use of TaxI2, as well as for the interpretation of the program's output.

Patterns of speciation in a parapatric pair of Saturnia moths as revealed by target capture

The focus of this study has been to understand the evolutionary relationships and taxonomy of a widely distributed parapatric species pair of wild silk moths in Europe: Saturnia pavonia and Saturnia pavoniella (Lepidoptera: Saturniidae). To address species delimitation in these parapatric taxa, target enrichment and mtDNA sequencing was employed alongside phylogenetic, admixture, introgression, and species delimitation analyses. The dataset included individuals from both species close to and farther away from the contact zone as well as two hybrids generated in the lab. Nuclear markers strongly supported both S. pavonia and S. pavoniella as two distinct species, with hybrids forming a sister group to S. pavoniella. However, the Maximum Likelihood (ML) tree generated from mtDNA sequencing data presented a different picture, showing both taxa to be phylogenetically intermixed. This inconsistency is likely attributable to mitonuclear discordance, which can arise from biological factors (e.g., introgressive hybridization and/or incomplete lineage sorting). Our analyses indicate that past introgressions have taken place, but that there is no evidence to suggest an ongoing admixture between the two species, demonstrating that the taxa have reached full postzygotic reproductive isolation and hence represent two distinct biological species. Finally, we discuss our results from an evolutionary point of view taking into consideration the past climatic oscillations that have likely shaped the present dynamics between the two species. Overall, our study demonstrates the effectiveness of the target enrichment approach in resolving shallow phylogenetic relationships under complex evolutionary circumstances and that this approach is useful in establishing robust and well-informed taxonomic delimitations involving parapatric taxa.

How Important Is Budding Speciation for Comparative Studies?

The acknowledgment of evolutionary dependence among species has fundamentally changed how we ask biological questions. Phylogenetic models became the standard approach for studies with 3 or more lineages, in particular those using extant species. Most phylogenetic comparative methods (PCMs) translate relatedness into covariance, meaning that evolutionary changes before lineages split should be interpreted together whereas after the split lineages are expected to change independently. This clever realization has shaped decades of research. Here, we discuss one element of the comparative method often ignored or assumed as unimportant: if nodes of a phylogeny represent the dissolution of the ancestral lineage into two new ones or if the ancestral lineage can survive speciation events (i.e., budding). Budding speciation is often reported in paleontological studies, due to the nature of the evidence for budding in the fossil record, but it is surprisingly absent in comparative methods. Here, we show that many PCMs assume that divergence happens as a symmetric split, even if these methods do not explicitly mention this assumption. We discuss the properties of trait evolution models for continuous and discrete traits and their adequacy under a scenario of budding speciation. We discuss the effects of budding speciation under a series of plausible evolutionary scenarios and show when and how these can influence our estimates. We also propose that long-lived lineages that have survived through a series of budding speciation events and given birth to multiple new lineages can produce evolutionary patterns that challenge our intuition about the most parsimonious history of trait changes in a clade. We hope our discussion can help bridge comparative approaches in paleontology and neontology as well as foster awareness about the assumptions we make when we use phylogenetic trees.

Phylogeography of Two Enigmatic Sulphur Butterflies, Colias mongola Alphéraky, 1897 and Colias tamerlana Staudinger, 1897 (Lepidoptera, Pieridae), with Relations to Wolbachia Infection

The genus Colias Fabricius, 1807 includes numerous taxa and forms with uncertain status and taxonomic position. Among such taxa are Colias mongola Alphéraky, 1897 and Colias tamerlana Staudinger, 1897, interpreted in the literature either as conspecific forms, as subspecies of different but morphologically somewhat similar Colias species or as distinct species-level taxa. Based on mitochondrial and nuclear DNA markers, we reconstructed a phylogeographic pattern of the taxa in question. We recover and include in our analysis DNA barcodes of the century-old type specimens, the lectotype of C. tamerlana deposited in the Natural History Museum (Museum für Naturkunde), Berlin, Germany (ZMHU) and the paralectotype of C. tamerlana and the lectotype of C. mongola deposited in the Zoological Institute, Russian Academy of Sciences, St. Petersburg, Russia (ZISP). Our analysis grouped all specimens within four (HP_I-HP_IV) deeply divergent but geographically poorly structured clades which did not support nonconspecifity of C. mongola-C. tamerlana. We also show that all studied females of the widely distributed haplogroup HP_II were infected with a single Wolbachia strain belonging to the supergroup B, while the males of this haplogroup, as well as all other investigated specimens of both sexes, were not infected. Our data highlight the relevance of large-scale sampling dataset analysis and the need for testing for Wolbachia infection to avoid erroneous phylogenetic reconstructions and species misidentification.

Phylomitogenomics bolsters the high-level classification of Demospongiae (phylum Porifera)

Class Demospongiae is the largest in the phylum Porifera (Sponges) and encompasses nearly 8,000 accepted species in three subclasses: Keratosa, Verongimorpha, and Heteroscleromorpha. Subclass Heteroscleromorpha contains ∼90% of demosponge species and is subdivided into 17 orders. The higher level classification of demosponges underwent major revision as the result of nearly three decades of molecular studies. However, because most of the previous molecular work only utilized partial data from a small number of nuclear and mitochondrial (mt) genes, this classification scheme needs to be tested by larger datasets. Here we compiled a mt dataset for 136 demosponge species-including 64 complete or nearly complete and six partial mt-genome sequences determined or assembled for this study-and used it to test phylogenetic relationships among Demospongiae in general and Heteroscleromorpha in particular. We also investigated the phylogenetic position of Myceliospongia araneosa, a highly unusual demosponge without spicules and spongin fibers, currently classified as Demospongiae incertae sedis, for which molecular data were not available. Our results support the previously inferred sister-group relationship between Heteroscleromorpha and Keratosa + Verongimorpha and suggest five main clades within Heteroscleromorpha: Clade C0 composed of order Haplosclerida; Clade C1 composed of Scopalinida, Sphaerocladina, and Spongillida; Clade C2 composed of Axinellida, Biemnida, Bubarida; Clade C3 composed of Tetractinellida; and Clade C4 composed of Agelasida, Clionaida, Desmacellida, Merliida, Suberitida, Poecilosclerida, Polymastiida, and Tethyida. The inferred relationships among these clades were (C0(C1(C2(C3+C4)))). Analysis of molecular data from M. araneosa placed it in the C3 clade as a sister taxon to the highly skeletonized tetractinellids Microscleroderma sp. and Leiodermatium sp. Molecular clock analysis dated divergences among the major clades in Heteroscleromorpha from the Cambrian to the Early Silurian, the origins of most heteroscleromorph orders in the middle Paleozoic, and the most basal splits within these orders around the Paleozoic to Mesozoic transition. Overall, the results of this study are mostly congruent with the accepted classification of Heteroscleromorpha, but add temporal perspective and new resolution to phylogenetic relationships within this subclass.

Assessing the phylogenetic relationship among varieties of Toona ciliata (Meliaceae) in sympatry with chloroplast genomes

Toona ciliata is an endangered species due to over-cutting and low natural regeneration in China. Its genetic conservation is of an increasing concern. However, several varieties are recognized according to the leaf and flower traits, which complicates genetic conservation of T. ciliata. Here, we sequenced the whole chloroplast genome sequences of three samples for each of four varieties (T. ciliata var. ciliata, T. ciliata var. yunnanensis, T. ciliata var. pubescens, and T. ciliata var. henryi) in sympatry and assessed their phylogenetic relationship at a fine spatial scale. The four varieties had genome sizes ranged from 159,546 to 159,617 bp and had small variations in genome structure. Phylogenomic analysis indicated that the four varieties were genetically well-mixed in branch groups. Genetic diversity from the whole chloroplast genome sequences of 12 samples was low among varieties (average π = 0.0003). Besides, we investigated genetic variation of 58 samples of the four varieties in sympatry using two markers (psaA and trnL-trnF) and showed that genetic differentiation was generally insignificant among varieties (Ф st = 0%-5%). Purifying selection occurred in all protein-coding genes except for the ycf2 gene that was under weak positive selection. Most amino acid sites in all protein-coding genes were under purifying selection except for a few sites that were under positive selection. The chloroplast genome-based phylogeny did not support the morphology-based classification. The overall results implicated that a conservation strategy based on the T. ciliata complex rather than on intraspecific taxon was more appropriate.

Inferring phylogenetic structure, taxa hybridization, and divergence times within rock voles of subgenus Aschizomys (Cricetidae: Alticola) using quaddRAD sequencing and a cytb dataset

The subgenus Aschizomys belongs to the genus Alticola (Central Asian mountain vole) and consists of two species: Alticola macrotis and Alticola lemminus. Phylogenetic relationships within the subgenus Aschizomys remain obscure due to limited sampling, an insufficient number of molecular markers used in phylogenetic studies, and paraphyly observed on mitochondrial trees. In this work, to infer reliable phylogenetic relationships and evaluate putative scenarios of ancient hybridization within the subgenus, we applied double-digest restriction site-associated DNA paired-end (quaddRAD) sequencing to 20 DNA samples (20 individuals), including five species of the genus Alticola, and dated the divergence of cytochrome b (cytb) lineages within Aschizomys using a "second calibration" approach. We showed monophyly of the two species on the basis of thousands of nuclear loci and demonstrated traces of introgression also in the nuclear genome. Observed paraphyly in cytb could be explained by an introgression event rather than incomplete lineage sorting. This explanation was confirmed by an analysis of the cytb divergence time. Overall, our results support the hypothesis of extensive migration of the Aschizomys species during the Late Pleistocene, with this migration leading to population divergence and introgression. We expect our article to become a starting point for a series of rigorous studies on the population history of the genus Alticola as a whole.

A genomic test of subspecies in the Eunota togata species group (Coleoptera: Cicindelidae): Morphology masks evolutionary relationships and taxonomy

Most of the world's biodiversity is described primarily or exclusively using morphological traits that may not always reflect the true evolutionary units. Accurate taxonomy is critical for conservation efforts and re-evaluation of traditional taxonomy may often be warranted since species and subspecies are frequently the focus of conservation and faunistic studies. Here, we test comprehensive taxonomic hypotheses of morphologically defined subspecies in the tiger beetle, Eunota togata (LaFerté-Sénectère, 1841). The four recognized subspecies were delineated based mainly on the dorsal coloration and extent of white markings termed maculations. We combine inferences from mtDNA genealogies and genome-wide multilocus data to elucidate the evolutionary relationships within the group and assess the taxonomic implications. Three of the four subspecific taxa delineated by morphology were not supported by the genomic or mtDNA data. In fact, the species-level diversity in this group was underestimated, as E. togata was found to represent three well-supported distinct species in all genetic analyses. Emerging from these analyses, we also document an intriguing example of convergent evolution in lighter colored E. togata adapting to similar white saline backgrounds. Our collective work underscores the importance of using molecular methods to reevaluate morphological based taxonomy for species and subspecies delimitation and conservation.

Baetidae (Insecta: Ephemeroptera) of Aurès Mountains (Algeria): A New Species of the Baetis alpinus Species Group, with Notes on Baetis Laech, 1815 Biogeography within Maghreb

A new species, Baetis (Baetis) dihyaesp. nov., belonging to the Baetis alpinus species group, is described and illustrated based on larval material collected in the Aurès Mountains (northeastern Algeria) in 2020-2021. This new species is closely related to three European species, e.g., Baetis (B.) alpinus (Pictet, 1843); B. (B.) nubecularis Eaton, 1898; and B. (B.) pasquetorum Righetti & Thomas, 2002 by the combination of the following characteristics: (i) more than one short, stout bristle at the tip of segment II of the maxillary palp and (ii) a well-developed paracercus. However, the new species clearly differs from all congeners of the Baetis alpinus species group primarily by the (a) structure of mouthparts-with 14-18 long submarginal setae arranged in a single irregular row on the dorsal surface of the labrum; 2-6 short, stout bristles at the tip of segment II of the maxillary palp; and segment II of the labial palp without a considerably developed apico-internal lobe); (b) setation of abdominal terga, with a few triangular-shaped scales sparsely scattered near the posterior margin only; and (c) a well-developed paracercus, comprised of more than 50 segments. Primary data on the biology and distribution of this new species are provided, and molecular affinities are verified by the analysis of COI (barcode) sequences. Detailed notes on the distribution of mayfly species belonging to the Baetis alpinus species group common in Western Europe and the western part of North Africa are presented. The historical movement of Baetis representatives between Europe, North West Africa, and subsequently Algeria, with the land bridges 'Strait of Gibraltar' and 'Strait of Sicily' as colonization routes, is discussed in detail and identified in the present study as the Western Algeria colonization path and Eastern Algeria colonization path, respectively.

Guidelines and quantitative standards for improved cetacean taxonomy using full mitochondrial genomes

In many organisms, especially those of conservation concern, traditional lines of evidence for taxonomic delineation, such as morphological data, are often difficult to obtain. In these cases, genetic data are often the only source of information available for taxonomic studies. In particular, population surveys of mitochondrial genomes offer increased resolution and precision in support of taxonomic decisions relative to conventional use of the control region or other gene fragments of the mitochondrial genome. To improve quantitative guidelines for taxonomic decisions in cetaceans, we build on a previous effort targeting the control region and evaluate, for whole mitogenome sequences, a suite of divergence and diagnosability estimates for pairs of recognized cetacean populations, subspecies, and species. From this overview, we recommend new guidelines based on complete mitogenomes, combined with other types of evidence for isolation and divergence, which will improve resolution for taxonomic decisions, especially in the face of small sample sizes or low levels of genetic diversity. We further use simulated data to assist interpretations of divergence in the context of varying forms of historical demography, culture, and ecology.

Genomic analysis reveals limited hybridization among three giraffe species in Kenya

BACKGROUND

In the speciation continuum, the strength of reproductive isolation varies, and species boundaries are blurred by gene flow. Interbreeding among giraffe (Giraffa spp.) in captivity is known, and anecdotal reports of natural hybrids exist. In Kenya, Nubian (G. camelopardalis camelopardalis), reticulated (G. reticulata), and Masai giraffe sensu stricto (G. tippelskirchi tippelskirchi) are parapatric, and thus, the country might be a melting pot for these taxa. We analyzed 128 genomes of wild giraffe, 113 newly sequenced, representing these three taxa.

RESULTS

We found varying levels of Nubian ancestry in 13 reticulated giraffe sampled across the Laikipia Plateau most likely reflecting historical gene flow between these two lineages. Although comparatively weaker signs of ancestral gene flow and potential mitochondrial introgression from reticulated into Masai giraffe were also detected, estimated admixture levels between these two lineages are minimal. Importantly, contemporary gene flow between East African giraffe lineages was not statistically significant. Effective population sizes have declined since the Late Pleistocene, more severely for Nubian and reticulated giraffe.

CONCLUSIONS

Despite historically hybridizing, these three giraffe lineages have maintained their overall genomic integrity suggesting effective reproductive isolation, consistent with the previous classification of giraffe into four species.

Deciphering complex reticulate evolution of Asian Buddleja (Scrophulariaceae): insights into the taxonomy and speciation of polyploid taxa in the Sino-Himalayan region

BACKGROUND AND AIMS

Species of the genus Buddleja in Asia are mainly distributed in the Sino-Himalayan region and form a challenging taxonomic group, with extensive hybridization and polyploidization. A phylogenetic approach to unravelling the history of reticulation in this lineage will deepen our understanding of the speciation in biodiversity hotspots.

METHODS

For this study, we obtained 80 accessions representing all the species in the Asian Buddleja clade, and the ploidy level of each taxon was determined by flow cytometry analyses. Whole plastid genomes, nuclear ribosomal DNA, single nucleotide polymorphisms and a large number of low-copy nuclear genes assembled from genome skimming data were used to investigate the reticulate evolutionary history of Asian Buddleja. Complex cytonuclear conflicts were detected through a comparison of plastid and species trees. Gene tree incongruence was also analysed to detect any reticulate events in the history of this lineage.

KEY RESULTS

Six hybridization events were detected, which are able to explain the cytonuclear conflict in Asian Buddleja. Furthermore, PhyloNet analysis combining species ploidy data indicated several allopolyploid speciation events. A strongly supported species tree inferred from a large number of low-copy nuclear genes not only corrected some earlier misinterpretations, but also indicated that there are many Asian Buddleja species that have been lumped mistakenly. Divergent time estimation shows two periods of rapid diversification (8-10 and 0-3 Mya) in the Asian Buddleja clade, which might coincide with the final uplift of the Hengduan Mountains and Quaternary climate fluctuations, respectively.

CONCLUSIONS

This study presents a well-supported phylogenetic backbone for the Asian Buddleja species, elucidates their complex and reticulate evolutionary history and suggests that tectonic activity, climate fluctuations, polyploidization and hybridization together promoted the diversification of this lineage.

DNA barcoding of Japanese earwig species (Insecta, Dermaptera), with sequence diversity analyses of three species of Anisolabididae

Dermaptera is a polyneopteran insect order that includes more than 2,000 described species, commonly known as earwigs, that mainly inhabit tropical, subtropical and warm temperate regions. Although 40 species have been found in Japan, their distribution and habitat preferences have remained ambiguous due to sample misidentification, particularly amongst immature specimens. To overcome this problem, we sequenced and analysed the DNA barcoding region of the mitochondrial cytochrome oxidase I gene (cox1) of dermapteran species recorded from Japan. Including publicly available data, 72.5% of known Japanese dermapteran species were subjected to molecular identification. We extensively sampled three wingless species of subfamily Anisolabidinae (Anisolabididae): Anisolabismaritima, Anisolabellamarginalis and Euborelliapallipes. Although these species exhibit similar habitat preferences as semi-synanthropes, A.maritima, a cosmopolitan species with the highest affinity to seashore, had significantly higher sequence diversity than the latter two species, which are considered endemic to East Asia. A similar trend was observed for (at least partly) winged cosmopolitan species of other families. Introgression with the congener Anisolabisseirokui is also suggested for A.maritima. Possible causes of the varying levels of sequence diversity are discussed.

From bites to barcodes: uncovering the hidden diversity of black flies (Diptera: Simuliidae) in Vietnam

BACKGROUND

Prompt and precise identification of black flies (Simuliidae) is crucial, given their biting behaviour and significant impact on human and animal health. To address the challenges presented by morphology and chromosomes in black fly taxonomy, along with the limited availability of molecular data pertaining to the black fly fauna in Vietnam, this study employed DNA-based approaches. Specifically, we used mitochondrial and nuclear-encoded genes to distinguish nominal species of black flies in Vietnam.

METHODS

In this study, 135 mitochondrial cytochrome c oxidase subunit I (COI) sequences were established for 45 species in the genus Simulium in Vietnam, encompassing three subgenera (Gomphostilbia, Nevermannia, and Simulium), with 64 paratypes of 27 species and 16 topotypes of six species. Of these COI sequences, 71, representing 27 species, are reported for the first time.

RESULTS

Combined with GenBank sequences of specimens from Malaysia, Myanmar, Thailand, and Vietnam, a total of 234 DNA barcodes of 53 nominal species resulted in a 71% success rate for species identification. Species from the non-monophyletic Simulium asakoae, S. feuerborni, S. multistriatum, S. striatum, S. tuberosum, and S. variegatum species groups were associated with ambiguous or incorrect identifications. Pairwise distances, phylogenetics, and species delimitation analyses revealed a high level of cryptic diversity, with discovery of 15 cryptic taxa. The current study also revealed the limited utility of a fast-evolving nuclear gene, big zinc finger (BZF), in discriminating closely related, morphologically similar nominal species of the S. asakoae species group.

CONCLUSION

This study represents the first comprehensive molecular genetic analysis of the black fly fauna in Vietnam to our knowledge, providing a foundation for future research. DNA barcoding exhibits varying levels of differentiating efficiency across species groups but is valuable in the discovery of cryptic diversity.

DNA barcoding of black flies (Diptera: Simuliidae) in Indonesia

BACKGROUND

DNA barcoding is a valuable taxonomic tool for rapid and accurate species identification and cryptic species discovery in black flies. Indonesia has 143 nominal species of black flies, but information on their biological aspects, including vectorial capacity and biting habits, remains underreported, in part because of identification problems. The current study represents the first comprehensive DNA barcoding of Indonesian black flies using mitochondrial cytochrome c oxidase subunit I (COI) gene sequences.

METHODS

Genomic DNA of Indonesian black fly samples were extracted and sequenced, producing 86 COI sequences in total. Two hundred four COI sequences, including 118 GenBank sequences, were analysed. Maximum likelihood (ML) and Bayesian inference (BI) trees were constructed and species delimitation analyses, including ASAP, GMYC and single PTP, were performed to determine whether the species of Indonesian black flies could be delineated. Intra- and interspecific genetic distances were also calculated and the efficacy of COI sequences for species identification was tested.

RESULTS

The DNA barcodes successfully distinguished most morphologically distinct species (> 80% of sampled taxa). Nonetheless, high maximum intraspecific distances (3.32-13.94%) in 11 species suggested cryptic diversity. Notably, populations of the common taxa Simulium (Gomphostilbia) cheongi, S. (Gomphostilbia) sheilae, S. (Nevermannia) feuerborni and S. (Simulium) tani in the islands of Indonesia were genetically distinct from those on the Southeast Asian mainland (Malaysia and Thailand). Integrated morphological, cytogenetic and nuclear DNA studies are warranted to clarify the taxonomic status of these more complex taxa.

CONCLUSIONS

The findings showed that COI barcoding is a promising taxonomic tool for Indonesian black flies. The DNA barcodes will aid in correct identification and genetic study of Indonesian black flies, which will be helpful in the control and management of potential vector species.

Mitochondrial introgression and mito-nuclear discordance obscured the closely related species boundaries in Cletus Stål from China (Heteroptera: Coreidae)

Accurate taxonomy and delimitation are of great importance for pest control strategies and management programs. Here, we focus on Cletus (Insecta: Hemiptera: Coreidae), which includes many crop pests. The species boundaries still conflict and only cytochrome c oxidase subunit I (COI) barcoding has been previously used for molecular studies. We generated new mitochondrial genome and nuclear genome-wide SNPs to explore the species boundaries of 46 Cletus samples from China using multiple species delimitation approaches. All results recovered a monophyly with high support, except for two closely related species in clade I - C. punctiger and C. graminis. Mitochondrial data demonstrated admixture in clade I, while genome-wide SNPs unambiguously identified two separate species, which were confirmed by morphological classification. Inconsistent nuclear and mitochondrial data indicated mito-nuclear discordance. Mitochondrial introgression is the most likely explanation, and more extensive sampling and more comprehensive data are needed to ascertain a pattern. Accurate species delimitation will shed light on species status; thus, an accurate taxonomy is of particular concern, as there is a pressing need to implement precise control of agricultural pests and to perform further research on diversification.

Evolutionary history of an Alpine Archaeognath (Machilis pallida): Insights from different variant

Reconstruction of species histories is a central aspect of evolutionary biology. Patterns of genetic variation within and among populations can be leveraged to elucidate evolutionary processes and demographic histories. However, interpreting genetic signatures and unraveling the contributing processes can be challenging, in particular for non-model organisms with complex reproductive modes and genome organization. One way forward is the combined consideration of patterns revealed by different molecular markers (nuclear vs. mitochondrial) and types of variants (common vs. rare) that differ in their age, mode, and rate of evolution. Here, we applied this approach to RNAseq data generated for Machilis pallida (Archaeognatha), an Alpine jumping bristletail considered parthenogenetic and triploid. We generated de novo transcriptome and mitochondrial assemblies to obtain high-density data to investigate patterns of mitochondrial and common and rare nuclear variation in 17 M. pallida individuals sampled from all known populations. We find that the different variant types capture distinct aspects of the evolutionary history and discuss the observed patterns in the context of parthenogenesis, polyploidy, and survival during glaciation. This study highlights the potential of different variant types to gain insights into evolutionary scenarios even from challenging but often available data and the suitability of M. pallida and the genus Machilis as a study system for the evolution of sexual strategies and polyploidization during environmental change. We also emphasize the need for further research which will be stimulated and facilitated by these newly generated resources and insights.

Delimitation and species discovery in the Profundulidae fish family: Using genetic, environmental and morphologic data to address taxonomic uncertainty

The family Profundulidae includes some one of the most enigmatic freshwater fishes of Mesoamerica: despite many attempts, a robust phylogenetic framework to delimit species is lacking, mainly due to limited morphological variation within the group. The accumulation of molecular data of profundulid fishes has led to advances in the description of new taxa, but relatively less progress has been made estimating evolutionary and phylogenetic relationships for this fish family. Here, we adopt an integrative taxonomy approach including the use of nuclear and mitochondrial DNA sequences, morphometric and ecological data, to test species boundaries in profundulid fishes in the westernmost area of their known distribution range in the states of Guerrero and Oaxaca, Mexico. Using a combination of methods for species discovery and validation based on Bayesian gene tree topologies, our analyses support the delimitation of 15 valid species of profundulid fishes - a combination of previously described species validated by this study, the synonymy of unsupported taxa, and the description of two new species. Using species delimitation methods, examination of phenotypic variation, and ecological niche characterization, we also identify five potentially new lineages which require further evidence to be erected as new species. We demonstrate that the use of an integrative taxonomy approach provides a robust methodology to delimit species in a taxonomically complex group like Profundulidae. Accurate taxonomic and ecological information is crucial for the conservation of these microendemic fishes, as several species are endangered.

Characteristics of plastid genomes in the genus Ceratostigma inhabiting arid habitats in China and their phylogenomic implications

BACKGROUND

Ceratostigma, a genus in the Plumbaginaceae, is an ecologically dominant group of shrubs, subshrub and herb mainly distributed in Qinghai-Tibet Plateau and North China. Ceratostigma has been the focal group in several studies, owing to their importance in economic and ecological value and unique breeding styles. Despite this, the genome information is limited and interspecific relationships within the genus Cerotastigma remains unexplored. Here we sequenced, assembled and characterized the 14 plastomes of five species, and conducted phylogenetic analyses of Cerotastigma using plastomes and nuclear ribosomal DNA (nrDNA) data.

RESULTS

Fourteen Cerotastigma plastomes possess typical quadripartite structures with lengths from 164,076 to 168,355 bp that consist of a large single copy, a small single copy and a pair of inverted repeats, and contain 127-128 genes, including 82-83 protein coding genes, 37 transfer RNAs and eight ribosomal RNAs. All plastomes are highly conservative and similar in gene order, simple sequence repeats (SSRs), long repeat repeats and codon usage patterns, but some structural variations in the border of single copy and inverted repeats. Mutation hotspots in coding (Pi values > 0.01: matK, ycf3, rps11, rps3, rpl22 and ndhF) and non-coding regions (Pi values > 0.02: trnH-psbA, rps16-trnQ, ndhF-rpl32 and rpl32-trnL) were identified among plastid genomes that could be served as potential molecular markers for species delimitation and genetic variation studies in Cerotastigma. Gene selective pressure analysis showed that most protein-coding genes have been under purifying selection except two genes. Phylogenetic analyses based on whole plastomes and nrDNA strongly support that the five species formed a monophyletic clade. Moreover, interspecific delimitation was well resolved except C. minus, individuals of which clustered into two main clades corresponding to their geographic distributions. The topology inferred from the nrDNA dataset was not congruent with the tree derived from the analyses of the plastid dataset.

CONCLUSION

These findings represent the first important step in elucidating plastome evolution in this widespread distribution genus Cerotastigma in the Qinghai-Tibet Plateau. The detailed information could provide a valuable resource for understanding the molecular dynamics and phylogenetic relationship in the family Plumbaginaceae. Lineage genetic divergence within C. minus was perhaps promoted by geographic barriers in the Himalaya and Hengduan Mountains region, but introgression or hybridization could not be completely excluded.

Population genomics indicate three different modes of divergence and speciation with gene flow in the green-winged teal duck complex

The processes leading to divergence and speciation can differ broadly among taxa with different life histories. We examine these processes in a small clade of ducks with historically uncertain relationships and species limits. The green-winged teal (Anas crecca) complex is a Holarctic species of dabbling duck currently categorized as three subspecies (Anas crecca crecca, A. c. nimia, and A. c. carolinensis) with a close relative, the yellow-billed teal (Anas flavirostris) from South America. A. c. crecca and A. c. carolinensis are seasonal migrants, while the other taxa are sedentary. We examined divergence and speciation patterns in this group, determining their phylogenetic relationships and the presence and levels of gene flow among lineages using both mitochondrial and genome-wide nuclear DNA obtained from 1,393 ultraconserved element (UCE) loci. Phylogenetic relationships using nuclear DNA among these taxa showed A. c. crecca, A. c. nimia, and A. c. carolinensis clustering together to form one polytomous clade, with A. flavirostris sister to this clade. This relationship can be summarized as (crecca, nimia, carolinensis)(flavirostris). However, whole mitogenomes revealed a different phylogeny: (crecca, nimia)(carolinensis, flavirostris). The best demographic model for key pairwise comparisons supported divergence with gene flow as the probable speciation mechanism in all three contrasts (crecca-nimia, crecca-carolinensis, and carolinensis-flavirostris). Given prior work, gene flow was expected among the Holarctic taxa, but gene flow between North American carolinensis and South American flavirostris (M ∼0.1-0.4 individuals/generation), albeit low, was not expected. Three geographically oriented modes of divergence are likely involved in the diversification of this complex: heteropatric (crecca-nimia), parapatric (crecca-carolinensis), and (mostly) allopatric (carolinensis-flavirostris). Our study shows that ultraconserved elements are a powerful tool for simultaneously studying systematics and population genomics in systems with historically uncertain relationships and species limits.

Taxonomic insights and evolutionary history in East Asian terrestrial slugs of the genus Meghimatium

East Asia, specifically the Japanese Archipelago, is a biodiversity hotspot of both vertebrates and invertebrates. Mollusks represent a burst of species diversity in this region due to the effects of biotic and abiotic factors on their morphological traits, such as shell shape and size. However, the evolutionary history of terrestrial slugs in East Asia remains unknown. In the present study, we investigated the molecular phylogeny of terrestrial slugs of the genus Meghimatium. This genus includes three described and eight undescribed species, and our study used all except for two. Based on phylogeny and the species delimitation tests, the genus Meghimatium was split into many putative species, suggesting higher species diversity than previously thought based on morphological and anatomical studies and that almost undescribed species may be inappropriate. Therefore, morphological traits, such as body size and colour, conventionally considered for classification may easily vary or be similar across geographic region. Moreover, the divergence time of this genus is almost concordant with the geographical time scale of the formation of the Japanese mainland. Our findings suggest that molecular phylogenetics helps classify Japanese Meghimatium slugs, but comprehensive taxonomic revisions using multi-locus analyses are needed.

Genomic insights into local adaptation in the Asiatic toad Bufo gargarizans, and its genomic offset to climate warming

Genomic signatures of local adaptation have been identified in many species but remain sparsely studied in amphibians. Here, we explored genome-wide divergence within the Asiatic toad, Bufo gargarizans, to study local adaptation and genomic offset (i.e., the mismatch between current and future genotype-environment relationships) under climate warming scenarios. We obtained high-quality SNP data for 94 Asiatic toads from 21 populations in China to study spatial patterns of genomic variation, local adaptation, and genomic offset to warming in this wide-ranging species. Population structure and genetic diversity analysis based on high-quality SNPs revealed three clusters of B. gargarizans in the western, central-eastern, and northeastern portions of the species' range in China. Populations generally dispersed along two migration routes, one from the west to the central-east and one from the central-east to the northeast. Both genetic diversity and pairwise F ST were climatically correlated, and pairwise F ST was also correlated with geographic distance. Spatial genomic patterns in B. gargarizans were determined by the local environment and geographic distance. Global warming will increase the extirpation risk of B. gargarizans.

Geographic-genomic and geographic-phenotypic differentiation of the Aquilegia viridiflora complex

How species diverge into different lineages is a central issue in evolutionary biology. Despite the increasing evidence indicating that such divergences do not need geographic isolation, the correlation between lineage divergence and the adaptive ecological divergence of phenotype corresponding to distribution is still unknown. In addition, gene flow has been widely detected during and through such diverging processes. We used one widely distributed Aquilegia viridiflora complex as a model system to examine genomic differentiation and corresponding phenotypic variations along geographic gradients. Our phenotypic analyses of 20 populations from northwest to northeast China identified two phenotypic groups along the geographic cline. All examined traits are distinct from each other, although a few intermediate individuals occur in their contacting regions. We further sequenced the genomes of representative individuals of each population. However, four distinct genetic lineages were detected based on nuclear genomes. In particular, we recovered numerous genetic hybrids in the contact regions of four lineages. Gene flow is widespread and continuous between four lineages but much higher between contacting lineages than geographically isolated lineages. Gene flow and natural selection might result in inconsistency between heredity and phenotype. Moreover, many genes with fast lineage-specific mutations were identified to be involved in local adaptation. Our results suggest that both geographic isolation and local selection exerted by the environment and pollinators may together create geographic distributions of phenotypic variations as well as the underlying genomic divergences in numerous lineages.