DNA barcodes combined with multi-locus data of representative taxa can generate reliable higher-level phylogenies

Molecular phylogeny of Hesperiidae (Lepidoptera) with an emphasis on Asian and African genera

Despite considerable research efforts in recent years, the deeper phylogenetic relationships among skipper butterflies (Hesperiidae) remain unresolved. This is primarily because of limited sampling, especially within Asian and African lineages. In this study, we consolidated previous data and extensively sampled Asian and African taxa to elucidate the phylogenetic relationships within Hesperiidae. The molecular dataset comprised sequences from two mitochondrial and two nuclear gene regions from 563 species that represented 353 genera. Our analyses revealed seven subfamilies within Hesperiidae: Coeliadinae, Euschemoninae, Eudaminae, Pyrginae, Heteropterinae, Trapezitinae, and Hesperiinae. The systematics of most tribes and genera aligned with those of prior studies. However, notable differences were observed in several tribes and genera. Overall, the position of taxa assigned to incertae sedis in Hesperiinae is largely clarified in this study. Our results strongly support the monophyly of the tribe Tagiadini (Pyrginae), and the systematics of some genera are clarified with comprehensive discussion. We recognize 15 tribes within the subfamily Hesperiinae. Of these, nine tribes are discussed in detail: Aeromachini, Astictopterini, Erionotini, Unkanini (new status), Ancistroidini, Ismini (confirmed status), Plastingini (new status), Gretnini (confirmed status), and Eetionini (confirmed status). We propose four subtribes within Astictopterini: Hypoleucina subtrib.n., Aclerosina, Cupithina, and Astictopterina. Furthermore, we describe a new genus (Hyarotoidesgen.n.) and reinstate two genera (Zeareinst.stat. and Separeinst.stat.) as valid. Additionally, we propose several new combinations: Zea mythecacomb.n.,Sepa bononiacomb.n. & reinst.stat., and Sepa umbrosacomb.n. This study, with extensive sampling of Asian and African taxa, greatly enhances the understanding of the knowledge of the skipper tree of life.

Unique among high passes: Phylogenetic inferences from DNA barcoding of the butterfly fauna of Ladakh Trans-Himalaya, India

The butterfly assemblage of Ladakh Trans-Himalaya demands a thorough analysis of their population genetic structure owing to their typical biogeographic affinity and their adaptability to extreme cold-desert climates. No such effort has been taken till date, and in this backdrop, we created a barcode reference library of 60 specimens representing 23 species. Barcodes were generated from freshly collected leg samples using the Sanger sequencing method, followed by phylogenetic clade analyses and divergence calculation. Our data represents 22% of Ladakh's Rhopaloceran fauna with the novel barcode submission for six species, including one Schedule II species, Paralasa mani . Contrary to the 3% threshold rule, the interspecific divergence between two species pairs of typical mountain genus Hyponephele and Karanasa was found to be 2.3% and 2.2%, respectively. The addition of conspecific global barcodes revealed that most species showed little increase in divergence value, while a two-fold increase was noted in a few species. Bayesian clade clustering outcomes largely aligned with current morphological classifications, forming monophyletic clades of conspecific barcodes, with only minor exceptions observed for the taxonomically complicated genus Polyommatus and misidentified records of Aulocera in the database. We also observed variations within the same phylogenetic clades forming nested lineages, which may be attributed to the taxonomic intricacies present at the subspecies level globally, mostly among Eurasian species. Overall, our effort not only substantiated the effectiveness of DNA Barcoding for the identification and conservation of this climatically vulnerable assemblage but also highlighted the significance of deciphering the unique genetic composition among this geographically isolated population of Ladakh butterflies.

Comprehensive phylogeny of Pieridae butterflies reveals strong correlation between diversification and temperature

Temperature is thought to be a key factor influencing global species richness patterns. We investigate the link between temperature and diversification in the butterfly family Pieridae by combining next generation DNA sequences and published molecular data with fine-grained distribution data. We sampled nearly 600 pierid butterfly species to infer the most comprehensive molecular phylogeny of the family and curated a distribution dataset of more than 800,000 occurrences. We found strong evidence that species in environments with more stable daily temperatures or cooler maximum temperatures in the warm seasons have higher speciation rates. Furthermore, speciation and extinction rates decreased in tandem with global temperatures through geological time, resulting in a constant net diversification.

Frog phylogeny: A time-calibrated, species-level tree based on hundreds of loci and 5,242 species

Large-scale, time-calibrated phylogenies from supermatrix studies have become crucial for evolutionary and ecological studies in many groups of organisms. However, in frogs (anuran amphibians), there is a serious problem with existing supermatrix estimates. Specifically, these trees are based on a limited number of loci (15 or fewer), and the higher-level relationships estimated are discordant with recent phylogenomic estimates based on much larger numbers of loci. Here, we attempted to rectify this problem by generating an expanded supermatrix and combining this with data from phylogenomic studies. To assist in aligning ribosomal sequences for this supermatrix, we developed a new program (TaxonomyAlign) to help perform taxonomy-guided alignments. The new combined matrix contained 5,242 anuran species with data from 307 markers, but with 95% missing data overall. This dataset represented a 71% increase in species sampled relative to the previous largest supermatrix analysis of anurans (adding 2,175 species). Maximum-likelihood analyses generated a tree in which higher-level relationships (and estimated clade ages) were generally concordant with those from phylogenomic analyses but were more discordant with the previous largest supermatrix analysis. We found few obvious problems arising from the extensive missing data in most species. We also generated a set of 100 time-calibrated trees for use in comparative analyses. Overall, we provide an improved estimate of anuran phylogeny based on the largest number of combined taxa and markers to date. More broadly, we demonstrate the potential to combine phylogenomic and supermatrix analyses in other groups of organisms.

Description of Uniphylus gen. nov., a new genus of Carcharodini (Lepidoptera, Hesperiidae, Pyrginae) for Staphylus evemerus Godman & Salvin, 1896

Staphylus evemerus Godman & Salvin, 1896 is a species with a unique set of morphological characters within Carcharodini. Also, mitochondrial cytochrome oxidase subunit 1 (COI) sequences analysis demonstrated a large genetic distance with other related genera of the tribe. Therefore, this paper aims to describe a new genus for this species, which is named as Uniphylus gen. nov. Besides the morphological redescription of the male of Uniphylus evemerus (Godman & Salvin, 1896) new comb., the description of the female is provided for the first time, as well as an updated distributional map with all records known so far for this species.

Molecular phylogeny, divergence time, biogeography and trends in host plant usage in the agriculturally important tortricid tribe Grapholitini (Lepidoptera: Tortricidae: Olethreutinae)

The leaf-roller moth tribe Grapholitini comprises about 1200 described species and contains numerous notorious pests of fruits and seeds. The phylogeny of the tribe has been little studied using contemporary methods, and the monophyly of several genera remains questionable. In order to provide a more robust phylogenetic framework for the group, we conducted a multiple-gene phylogenetic analysis of 104 species representing 27 genera of Grapholitini and 29 outgroup species. Divergence time, ancestral area, and host plant usage were also inferred to explore evolutionary trends in the tribe. Our analyses indicate that Larisa and Corticivora, traditionally assigned to Grapholitini, are best excluded from the tribe. After removal of these two genera, the tribe is found to be monophyletic, represented by two major lineages-a Dichrorampha clade and a Cydia clade, the latter of which can be divided into seven generic groups. The genus Grapholita was found to be polyphyletic, comprising three different clades, and we propose three genera to accommodate these groups: Grapholita (sensu stricto), Aspila (formerly a subgenus of Grapholita) and Ephippiphora (formerly considered a synonym of Grapholita). We summarize each generic group, including related genera not included in our analysis, providing morphological, pheromone and food plant characters that support particular branches within the molecular hypotheses. Biogeographical analyses indicate that Grapholitini probably originated in the Nearctic, Afrotropical and Neotropical regions in the Lutetian of the middle Eocene (ca. 44.3 Ma). Our results also indicate that most groups in Grapholitini originated from Fabaceae-feeding monophagous or oligophagous ancestors, and that host plant shifts probably promoted species diversification within the tribe.

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.

Phylogenetic Structure Revealed through Combining DNA Barcodes with Multi-Gene Data for Agrodiaetus Blue Butterflies (Lepidoptera, Lycaenidae)

The need for multi-gene analysis in evolutionary and taxonomic studies is generally accepted. However, the sequencing of multiple genes is not always possible. For various reasons, short mitochondrial DNA barcodes are the only source of molecular information for some species in many genera, although multi-locus data are available for other species of the same genera. In particular, such situation exists in the species-rich butterfly subgenus Polyommatus (Agrodiaetus). Here, we analyzed the partitioning of this subgenus into species groups by using three data sets. The first data set was represented by short mitochondrial DNA barcodes for all analyzed samples. The second and third data sets were represented by a combination of short mitochondrial DNA barcodes for part of the taxa with longer mitochondrial sequences COI + tRNA-Leu + COII (data set 2) and with longer mitochondrial COI + tRNA-Leu + COII and nuclear 5.8S rDNA + ITS2 + 28S rDNA sequences (data set 3) for the remaining species. We showed that the DNA barcoding approach (data set 1) failed to reveal the phylogenetic structure, resulting in numerous polytomies in the tree obtained. Combined analysis of the mitochondrial and nuclear sequences (data sets 2 and 3) revealed the species groups and the position within these species groups, even for the taxa for which only short DNA barcodes were available.

Chromosomal conservatism vs chromosomal megaevolution: enigma of karyotypic evolution in Lepidoptera

In the evolution of many organisms, periods of slow genome reorganization (= chromosomal conservatism) are interrupted by bursts of numerous chromosomal changes (= chromosomal megaevolution). Using comparative analysis of chromosome-level genome assemblies, we investigated these processes in blue butterflies (Lycaenidae). We demonstrate that the phase of chromosome number conservatism is characterized by the stability of most autosomes and dynamic evolution of the sex chromosome Z, resulting in multiple variants of NeoZ chromosomes due to autosome-sex chromosome fusions. In contrast during the phase of rapid chromosomal evolution, the explosive increase in chromosome number occurs mainly through simple chromosomal fissions. We show that chromosomal megaevolution is a highly non-random canalized process, and in two phylogenetically independent Lysandra lineages, the drastic parallel increase in number of fragmented chromosomes was achieved, at least partially, through reuse of the same ancestral chromosomal breakpoints. In species showing chromosome number doubling, we found no blocks of duplicated sequences or duplicated chromosomes, thus refuting the hypothesis of polyploidy. In the studied taxa, long blocks of interstitial telomere sequences (ITSs) consist of (TTAGG)_n arrays interspersed with telomere-specific retrotransposons. ITSs are sporadically present in rapidly evolving Lysandra karyotypes, but not in the species with ancestral chromosome number. Therefore, we hypothesize that the transposition of telomeric sequences may be triggers of the rapid chromosome number increase. Finally, we discuss the hypothetical genomic and population mechanisms of chromosomal megaevolution and argue that the disproportionally high evolutionary role of the Z sex chromosome can be additionally reinforced by sex chromosome-autosome fusions and Z-chromosome inversions.

Diversification and evolutionary history of the African laminated-toothed rats (Rodentia, Otomyini)

The African continent was subjected to periodic climatic shifts during the Pliocene and Pleistocene. These habitat changes greatly affected the evolutionary processes and tempo of diversification in numerous, widely distributed mammals. The Otomyini (Family Muridae) comprises three African rodent genera, Parotomys, Otomys and Myotomys, characterized by unique laminated-shaped molars. Species within this tribe generally prefer open-habitat and show low dispersal capabilities, with previous studies suggesting that their diversification was closely associated with climatic oscillations over the last four million years. Our phylogenetic reconstructions, based on three mitochondrial (mtDNA) genes (Cytb, COI and 12S) and four nuclear introns (EF, SPTBN, MGF and THY), identified eight major genetic clades that are distributed across southern, eastern and western Africa. Our data permit the re-examination of the taxonomic status of the three genera as well as the previously proposed mesic-arid dichotomy of the 10 South African species. Moreover, multiple mtDNA species delimitation methods incorporating 168 specimens estimated the number of Otomyini species to be substantially higher than the ∼30 recognized, suggesting that the current taxonomy will necessitate an integrative approach to delimit extant species diversity within the Otomyini. The data suggests that the origin of the tribe can be dated back to ∼5.7 million years ago (Ma) in southern Africa. The distribution and phylogenetic associations among the eight major otomyine evolutionary lineages can best be explained by several waves of northward colonization from southern Africa, complemented by independent reversed dispersals from eastern back to southern Africa at different time periods. There is strong support for the hypothesis that the radiation, dispersion, and diversification of the otomyine rodents is closely linked to recent Plio-Pleistocene climatic oscillations.

DNA barcoding of Chinese snakes reveals hidden diversity and conservation needs

DNA barcoding has greatly facilitated studies of taxonomy, biodiversity, biological conservation, and ecology. Here, we establish a reliable DNA barcoding library for Chinese snakes, unveiling hidden diversity with implications for taxonomy, and provide a standardized tool for conservation management. Our comprehensive study includes 1638 cytochrome c oxidase subunit I (COI) sequences from Chinese snakes that correspond to 17 families, 65 genera, 228 named species (80.6% of named species) and 36 candidate species. A barcode gap analysis reveals gaps, where all nearest neighbour distances exceed maximum intraspecific distances, in 217 named species and all candidate species. Three species-delimitation methods (ABGD, sGMYC, and sPTP) recover 320 operational taxonomic units (OTUs), of which 192 OTUs correspond to named and candidate species. Twenty-eight other named species share OTUs, such as Azemiops feae and A. kharini, Gloydius halys, G. shedaoensis, and G. intermedius, and Bungarus multicinctus and B. candidus, representing inconsistencies most probably caused by imperfect taxonomy, recent and rapid speciation, weak taxonomic signal, introgressive hybridization, and/or inadequate phylogenetic signal. In contrast, 43 species and candidate species assign to two or more OTUs due to having large intraspecific distances. If most OTUs detected in this study reflect valid species, including the 36 candidate species, then 30% more species would exist than are currently recognized. Several OTU divergences associate with known biogeographic barriers, such as the Taiwan Strait. In addition to facilitating future studies, this reliable and relatively comprehensive reference database will play an important role in the future monitoring, conservation, and management of Chinese snakes.

Whole-Genome Analysis Reveals the Dynamic Evolution of Holocentric Chromosomes in Satyrine Butterflies

Butterfly chromosomes are holocentric, i.e., lacking a localized centromere. Potentially, this can lead to rapid karyotypic evolution through chromosome fissions and fusions, since fragmented chromosomes retain kinetic activity, while fused chromosomes are not dicentric. However, the actual mechanisms of butterfly genome evolution are poorly understood. Here, we analyzed chromosome-scale genome assemblies to identify structural rearrangements between karyotypes of satyrine butterfly species. For the species pair Erebia ligea-Maniola jurtina, sharing the ancestral diploid karyotype 2n = 56 + ZW, we demonstrate a high level of chromosomal macrosynteny and nine inversions separating these species. We show that the formation of a karyotype with a low number of chromosomes (2n = 36 + ZW) in Erebia aethiops was based on ten fusions, including one autosome-sex chromosome fusion, resulting in a neo-Z chromosome. We also detected inversions on the Z sex chromosome that were differentially fixed between the species. We conclude that chromosomal evolution is dynamic in the satyrines, even in the lineage that preserves the ancestral chromosome number. We hypothesize that the exceptional role of Z chromosomes in speciation may be further enhanced by inversions and sex chromosome-autosome fusions. We argue that not only fusions/fissions but also inversions are drivers of the holocentromere-mediated mode of chromosomal speciation.

Molecular Phylogeny and Taxonomy of the Butterfly Subtribe Scolitantidina with Special Focus on the Genera Pseudophilotes, Glaucopsyche and Iolana (Lepidoptera, Lycaenidae)

The Palearctic blue butterfly genus Pseudophilotes Beuret, 1958 is not homogenous regarding the morphology of its genital structures. For this reason, some of its species have been considered to be representatives of other genera of the subtribe Scolitantidina (subfamily Polyommatinae). Here, we address these taxonomic problems by analyzing the phylogenetic relationships between the genera, subgenera, and species of this subtribe inferred via the analysis of five nuclear and two mitochondrial DNA sequences. We demonstrate that the enigmatic Asian species P. panope (Eversmann, 1851) belongs to the genus Pseudophilotes but not to Praephilotes Forster, 1938 or Palaeophilotes Forster, 1938 and does not represent the independent genus Inderskia Korshunov, 2000, as hypothesized previously. We synonymize P. svetlana Yakovlev, 2003 (syn. nov.) and P. marina Zhdanko, 2004 (syn. nov.) with P. panope. We demonstrate a deep genetic divergence between lineages that were previously considered as subspecies of the single species Iolana iolas (Ochsenheimer, 1816). As a result, we confirm the multispecies concept of the genus Iolana Bethune-Baker, 1914. We show that the Holarctic genus Glaucopsyche can be divided into four subgenera: Glaucopsyche Scudder, 1872 (=Shijimiaeoides Beuret, 1958), Apelles Hemming, 1931, Bajluana Korshunov and Ivonin, 1990, and Phaedrotes Scudder, 1876.

Predictors of sequence capture in a large-scale anchored phylogenomics project

Next-generation sequencing (NGS) technologies have revolutionized phylogenomics by decreasing the cost and time required to generate sequence data from multiple markers or whole genomes. Further, the fragmented DNA of biological specimens collected decades ago can be sequenced with NGS, reducing the need for collecting fresh specimens. Sequence capture, also known as anchored hybrid enrichment, is a method to produce reduced representation libraries for NGS sequencing. The technique uses single-stranded oligonucleotide probes that hybridize with pre-selected regions of the genome that are sequenced via NGS, culminating in a dataset of numerous orthologous loci from multiple taxa. Phylogenetic analyses using these sequences have the potential to resolve deep and shallow phylogenetic relationships. Identifying the factors that affect sequence capture success could save time, money, and valuable specimens that might be destructively sampled despite low likelihood of sequencing success. We investigated the impacts of specimen age, preservation method, and DNA concentration on sequence capture (number of captured sequences and sequence quality) while accounting for taxonomy and extracted tissue type in a large-scale butterfly phylogenomics project. This project used two probe sets to extract 391 loci or a subset of 13 loci from over 6,000 butterfly specimens. We found that sequence capture is a resilient method capable of amplifying loci in samples of varying age (0–111 years), preservation method (alcohol, papered, pinned), and DNA concentration (0.020 ng/μl - 316 ng/ul). Regression analyses demonstrate that sequence capture is positively correlated with DNA concentration. However, sequence capture and DNA concentration are negatively correlated with sample age and preservation method. Our findings suggest that sequence capture projects should prioritize the use of alcohol-preserved samples younger than 20 years old when available. In the absence of such specimens, dried samples of any age can yield sequence data, albeit with returns that diminish with increasing age.

An open and continuously updated fern tree of life

Ferns, with about 12,000 species, are the second most diverse lineage of vascular plants after angiosperms. They have been the subject of numerous molecular phylogenetic studies, resulting in the publication of trees for every major clade and DNA sequences from nearly half of all species. Global fern phylogenies have been published periodically, but as molecular systematics research continues at a rapid pace, these become quickly outdated. Here, we develop a mostly automated, reproducible, open pipeline to generate a continuously updated fern tree of life (FTOL) from DNA sequence data available in GenBank. Our tailored sampling strategy combines whole plastomes (few taxa, many loci) with commonly sequenced plastid regions (many taxa, few loci) to obtain a global, species-level fern phylogeny with high resolution along the backbone and maximal sampling across the tips. We use a curated reference taxonomy to resolve synonyms in general compliance with the community-driven Pteridophyte Phylogeny Group I classification. The current FTOL includes 5,582 species, an increase of ca. 40% relative to the most recently published global fern phylogeny. Using an updated and expanded list of 51 fern fossil constraints, we find estimated ages for most families and deeper clades to be considerably older than earlier studies. FTOL and its accompanying datasets, including the fossil list and taxonomic database, will be updated on a regular basis and are available via a web portal (https://fernphy.github.io) and R packages, enabling immediate access to the most up-to-date, comprehensively sampled fern phylogeny. FTOL will be useful for anyone studying this important group of plants over a wide range of taxonomic scales, from smaller clades to the entire tree. We anticipate FTOL will be particularly relevant for macroecological studies at regional to global scales and will inform future taxonomic systems with the most recent hypothesis of fern phylogeny.

Climate, immigration and speciation shape terrestrial and aquatic biodiversity in the European Alps

Quaternary climate fluctuations can affect speciation in regional biodiversity assembly in two non-mutually exclusive ways: a glacial species pump, where isolation in glacial refugia accelerates allopatric speciation, and adaptive radiation in underused adaptive zones during ice-free periods. We detected biogeographic and genetic signatures associated with both mechanisms in the assembly of the biota of the European Alps. Age distributions of endemic and widespread species within aquatic and terrestrial taxa (amphipods, fishes, amphibians, butterflies and flowering plants) revealed that endemic fish evolved only in lakes, are highly sympatric, and mainly of Holocene age, consistent with adaptive radiation. Endemic amphipods are ancient, suggesting preglacial radiation with limited range expansion and local Pleistocene survival, perhaps facilitated by a groundwater-dwelling lifestyle. Terrestrial endemics are mostly of Pleistocene age and are thus more consistent with the glacial species pump. The lack of evidence for Holocene adaptive radiation in the terrestrial biome is consistent with faster recolonization through range expansion of these taxa after glacial retreats. More stable and less seasonal ecological conditions in lakes during the Holocene may also have contributed to Holocene speciation in lakes. The high proportion of young, endemic species makes the Alpine biota vulnerable to climate change, but the mechanisms and consequences of species loss will likely differ between biomes because of their distinct evolutionary histories.

Butterfly–parasitoid–hostplant interactions in Western Palaearctic Hesperiidae: a DNA barcoding reference library

The study of ecological interactions between plants, phytophagous insects and their natural enemies is an essential but challenging component for understanding ecosystem dynamics. Molecular methods such as DNA barcoding can help elucidate these interactions. In this study, we employed DNA barcoding to establish hostplant and parasitoid interactions with hesperiid butterflies, using a complete reference library for Hesperiidae of continental Europe and north-western Africa (53 species, 100% of those recorded) based on 2934 sequences from 38 countries. A total of 233 hostplant and parasitoid interactions are presented, some recovered by DNA barcoding larval remains or parasitoid cocoons. Combining DNA barcode results with other lines of evidence allowed 94% species-level identification for Hesperiidae, but success was lower for parasitoids, in part due to unresolved taxonomy. Potential cases of cryptic diversity, both in Hesperiidae and Microgastrinae, are discussed. We briefly analyse the resulting interaction networks. Future DNA barcoding initiatives in this region should focus attention on north-western Africa and on parasitoids, because in these cases barcode reference libraries and taxonomy are less well developed.

Organelle Phylogenomics and Extensive Conflicting Phylogenetic Signals in the Monocot Order Poales

The Poales is one of the largest orders of flowering plants with significant economic and ecological values. Reconstructing the phylogeny of the Poales is important for understanding its evolutionary history that forms the basis for biological studies. However, due to sparse taxon sampling and limited molecular data, previous studies have resulted in a variety of contradictory topologies. In particular, there are three nodes surrounded by incongruence: the phylogenetic ambiguity near the root of the Poales tree, the sister family of Poaceae, and the delimitation of the xyrid clade. We conducted a comprehensive sampling and reconstructed the phylogenetic tree using plastid and mitochondrial genomic data from 91 to 66 taxa, respectively, representing all the 16 families of Poales. Our analyses support the finding of Bromeliaceae and Typhaceae as the earliest diverging groups within the Poales while having phylogenetic relationships with the polytomy. The clade of Ecdeiocoleaceae and Joinvilleaceae is recovered as the sister group of Poaceae. The three families, Mayacaceae, Eriocaulaceae, and Xyridaceae, of the xyrid assembly diverged successively along the backbone of the Poales phylogeny, and thus this assembly is paraphyletic. Surprisingly, we find substantial phylogenetic conflicts within the plastid genomes of the Poales, as well as among the plastid, mitochondrial, and nuclear data. These conflicts suggest that the Poales could have a complicated evolutionary history, such as rapid radiation and polyploidy, particularly allopolyploidy through hybridization. In sum, our study presents a new perspicacity into the complex phylogenetic relationships and the underlying phylogenetic conflicts within the Poales.

Phylogenomic and mitogenomic data can accelerate inventorying of tropical beetles during the current biodiversity crisis

Conservation efforts must be evidence-based, so rapid and economically feasible methods should be used to quantify diversity and distribution patterns. We have attempted to overcome current impediments to the gathering of biodiversity data by using integrative phylogenomic and three mtDNA fragment analyses. As a model, we sequenced the Metriorrhynchini beetle fauna, sampled from ~700 localities in three continents. The species-rich dataset included ~6,500 terminals, ~1,850 putative species delimited at 5% uncorrected pairwise threshold, possibly ~1,000 of them unknown to science. Neither type of data could alone answer our questions on biodiversity and phylogeny. The phylogenomic backbone enabled the integrative delimitation of robustly defined natural genus-group units that will inform future research. Using constrained mtDNA analysis, we identified the spatial structure of species diversity, very high species-level endemism, and a biodiversity hotspot in New Guinea. We suggest that focused field research and subsequent laboratory and bioinformatic workflow steps would substantially accelerate the inventorying of any hyperdiverse tropical group with several thousand species. The outcome would be a scaffold for the incorporation of further data from environmental sequencing and ecological studies. The database of sequences could set a benchmark for the spatiotemporal evaluation of biodiversity, would support evidence-based conservation planning, and would provide a robust framework for systematic, biogeographic, and evolutionary studies.

Phylogeny of gracillariid leaf-mining moths: evolution of larval behaviour inferred from phylogenomic and Sanger data

Gracillariidae is the most taxonomically diverse cosmopolitan leaf-mining moth family, consisting of nearly 2000 named species in 105 described genera, classified into eight extant subfamilies. The majority of gracillariid species are internal plant feeders as larvae, creating mines and galls in plant tissue. Despite their diversity and ecological adaptations, their phylogenetic relationships, especially among subfamilies, remain uncertain. Genomic data (83 taxa, 589 loci) were integrated with Sanger data (130 taxa, 22 loci), to reconstruct a phylogeny of Gracillariidae. Based on analyses of both datasets combined and analyzed separately, monophyly of Gracillariidae and all its subfamilies, monophyly of the clade "LAMPO" (subfamilies: Lithocolletinae, Acrocercopinae, Marmarinae, Phyllocnistinae, and Oecophyllembiinae) and relationships of its subclade "AMO" (subfamilies: Acrocercopinae, Marmarinae, and Oecophyllembiinae) were strongly supported. A sister-group relationship of Ornixolinae to the remainder of the family, and a monophyletic leaf roller lineage (Callicercops Vári + Parornichinae) + Gracillariinae, as sister to the "LAMPO" clade were supported by the most likely tree. Dating analyses indicate a mid-Cretaceous (105.3 Ma) origin of the family, followed by a rapid diversification into the nine subfamilies predating the Cretaceous-Palaeogene extinction. We hypothesize that advanced larval behaviours, such as making keeled or tentiform blotch mines, rolling leaves and galling, allowed gracillariids to better avoid larval parasitoids allowing them to further diversify. Finally, we stabilize the classification by formally re-establishing the subfamily ranks of Marmarinae stat.rev., Oecophyllembiinae stat.rev. and Parornichinae stat.rev., and erect a new subfamily, Callicercopinae Li, Ohshima and Kawahara to accommodate the enigmatic genus Callicercops.