Genomic Consequences of and Demographic Response to Pervasive Hybridization Over Time in Climate-Sensitive Pikas.

Rare and geographically restricted species may be vulnerable to genetic effects from inbreeding depression in small populations or from genetic swamping through hybridization with common species, but a third possibility is that selective gene flow can restore fitness (genetic rescue). Climate-sensitive pikas (Ochotona spp.) of the Qinghai-Tibetan Plateau (QHTP) and its vicinity have been reduced to residual populations through the movement of climatic zones during the Pleistocene and recent anthropogenic disturbance, whereas the plateau pika (O. curzoniae) remains common. Population-level whole-genome sequencing (n = 142) of six closely related species in the subgenus Ochotona revealed several phases of ancient introgression, lineage replacement, and bidirectional introgression. The strength of gene flow was the greatest from the dominant O. curzoniae to ecologically distinct species in areas peripheral to the QHTP. Genetic analyses were consistent with environmental reconstructions of past population movements. Recurrent periods of introgression throughout the Pleistocene revealed an increase in genetic variation at first but subsequent loss of genetic variation in later phases. Enhanced dispersion of introgressed genomic regions apparently contributed to demographic recovery in three peripheral species that underwent range shifts following climate oscillations on the QHTP, although it failed to drive recovery of northeastern O. dauurica and geographically isolated O. sikimaria. Our findings highlight differences in timescale and environmental background to determine the consequence of hybridization and the unique role of the QHTP in conserving key evolutionary processes of sky island species.

Phylogenetic diversity only weakly mitigates climate-change-driven biodiversity loss in insect communities.

To help address the underrepresentation of arthropods and Asian biodiversity from climate-change assessments, we carried out year-long, weekly sampling campaigns with Malaise traps at different elevations and latitudes in Gaoligongshan National Park in southwestern China. From these 623 samples, we barcoded 10,524 beetles and compared scenarios of climate-change-induced biodiversity loss, by designating seasonal, elevational, and latitudinal subsets of beetles as communities that plausibly could go extinct as a group, which we call "loss sets". The availability of a published mitochondrial-genome-based phylogeny of the Coleoptera allowed us to compare the loss of species diversity with and without accounting for phylogenetic relatedness. We hypothesised that phylogenetic relatedness would mitigate extinction, since the extinction of any loss set would result in the disappearance of all its species but only part of its evolutionary history, which is still extant in the remaining loss sets. We found different patterns of community clustering by season and latitude, depending on whether phylogenetic information was incorporated. However, accounting for phylogeny only slightly mitigated the amount of biodiversity loss under climate change scenarios, against our expectations: there is no phylogenetic "escape clause" for biodiversity conservation. We achieve the same results whether phylogenetic information was derived from the mitogenome phylogeny or from a de novo barcode-gene tree. We encourage interested researchers to use this data set to study lineage-specific community assembly patterns in conjunction with life-history traits and environmental covariates.

Inferring the ecological and evolutionary determinants of community genetic diversity.

Understanding the relative contributions of ecological and evolutionary processes to the structuring of ecological communities is needed to improve our ability to predict how communities may respond to future changes in an increasingly human-modified world. Metabarcoding methods make it possible to gather population genetic data for all species within a community, unlocking a new axis of data to potentially unveil the origins and maintenance of biodiversity at local scales. Here, we present a new eco-evolutionary simulation model for investigating community assembly dynamics using metabarcoding data. The model makes joint predictions of species abundance, genetic variation, trait distributions and phylogenetic relationships under a wide range of parameter settings (e.g. high speciation/low dispersal or vice versa) and across a range of community states, from pristine and unmodified to heavily disturbed. We first demonstrate that parameters governing metacommunity and local community processes leave detectable signatures in simulated biodiversity data axes. Next, using a simulation-based machine learning approach we show that neutral and non-neutral models are distinguishable and that reasonable estimates of several model parameters within the local community can be obtained using only community-scale genetic data, while phylogenetic information is required to estimate those describing metacommunity dynamics. Finally, we apply the model to soil microarthropod metabarcoding data from the Troodos mountains of Cyprus, where we find that communities in widespread forest habitats are structured by neutral processes, while high-elevation and isolated habitats act as an abiotic filter generating non-neutral community structure. We implement our model within the ibiogen R package, a package dedicated to the investigation of island, and more generally community-scale, biodiversity using community-scale genetic data.

Community metabarcoding reveals the relative role of environmental filtering and spatial processes in metacommunity dynamics of soil microarthropods across a mosaic of montane forests.

Disentangling the relative role of environmental filtering and spatial processes in driving metacommunity structure across mountainous regions remains challenging, as the way we quantify spatial connectivity in topographically and environmentally heterogeneous landscapes can influence our perception of which process predominates. More empirical data sets are required to account for taxon- and context-dependency, but relevant research in understudied areas is often compromised by the taxonomic impediment. Here we used haplotype-level community DNA metabarcoding, enabled by stringent filtering of amplicon sequence variants (ASVs), to characterize metacommunity structure of soil microarthropod assemblages across a mosaic of five forest habitats on the Troodos mountain range in Cyprus. We found similar ß diversity patterns at ASV and species (OTU, operational taxonomic unit) levels, which pointed to a primary role of habitat filtering resulting in the existence of largely distinct metacommunities linked to different forest types. Within-habitat turnover was correlated to topoclimatic heterogeneity, again emphasizing the role of environmental filtering. However, when integrating landscape matrix information for the highly fragmented Quercus alnifolia habitat, we also detected a major role of spatial isolation determined by patch connectivity, indicating that stochastic and niche-based processes synergistically govern community assembly. Alpha diversity patterns varied between ASV and OTU levels, with OTU richness decreasing with elevation and ASV richness following a longitudinal gradient, potentially reflecting a decline of genetic diversity eastwards due to historical pressures. Our study demonstrates the utility of haplotype-level community metabarcoding for characterizing metacommunity structure of complex assemblages and improving our understanding of biodiversity dynamics across mountainous landscapes worldwide.

Collective and harmonized high throughput barcoding of insular arthropod biodiversity: Toward a Genomic Observatories Network for islands.

Current understanding of ecological and evolutionary processes underlying island biodiversity is heavily shaped by empirical data from plants and birds, although arthropods comprise the overwhelming majority of known animal species, and as such can provide key insights into processes governing biodiversity. Novel high throughput sequencing (HTS) approaches are now emerging as powerful tools to overcome limitations in the availability of arthropod biodiversity data, and hence provide insights into these processes. Here, we explored how these tools might be most effectively exploited for comprehensive and comparable inventory and monitoring of insular arthropod biodiversity. We first reviewed the strengths, limitations and potential synergies among existing approaches of high throughput barcode sequencing. We considered how this could be complemented with deep learning approaches applied to image analysis to study arthropod biodiversity. We then explored how these approaches could be implemented within the framework of an island Genomic Observatories Network (iGON) for the advancement of fundamental and applied understanding of island biodiversity. To this end, we identified seven island biology themes at the interface of ecology, evolution and conservation biology, within which collective and harmonized efforts in HTS arthropod inventory could yield significant advances in island biodiversity research.

The phylogeny and evolutionary ecology of hoverflies (Diptera: Syrphidae) inferred from mitochondrial genomes.

Hoverflies (Diptera: Syrphidae) are a diverse group of pollinators and a major research focus in ecology, but their phylogenetic relationships remain incompletely known. Using a genome skimming approach we generated mitochondrial genomes for 91 species, capturing a wide taxonomic diversity of the family. To reduce the required amount of input DNA and overall cost of the library construction, sequencing and assembly was conducted on mixtures of specimens, which raises the problem of chimera formation of mitogenomes. We present a novel chimera detection test based on gene tree incongruence, but identified only a single mitogenome of chimeric origin. Together with existing data for a final set of 127 taxa, phylogenetic analysis on nucleotide and amino acid sequences using Maximum Likelihood and Bayesian Inference revealed a basal split of Microdontinae from all other syrphids. The remainder consists of several deep clades assigned to the subfamily Eristalinae in the current classification, including a clade comprising the subfamily Syrphinae (plus Pipizinae). These findings call for a re-definition of subfamilies, but basal nodes had insufficient support to fully justify such action. Molecular-clock dating placed the origin of the Syrphidae crown group in the mid-Cretaceous while the Eristalinae-Syrphinae clade likely originated near the K/Pg boundary. Transformation of larval life history characters on the tree suggests that Syrphidae initially had sap feeding larvae, which diversified greatly in diet and habitat association during the Eocene and Oligocene, coinciding with the diversification of angiosperms and the evolution of various insect groups used as larval host, prey, or mimicry models. Mitogenomes proved to be a powerful phylogenetic marker for studies of Syrphidae at subfamily and tribe levels, allowing dense taxon sampling that provided insight into the great ecological diversity and rapid evolution of larval life history traits of the hoverflies.

Untangling the evolution of soldier beetles (Coleoptera: Cantharidae) and the evaluation of the morphological phylogenetic signal in a soft-bodied elateroid lineage.

This study addresses the long-standing uncertainty about the internal classification of soldier beetles (Elateroidea: Cantharidae). Four datasets were compiled and analysed: 66 genes for 14 terminals, 15 mtDNA genes for 79 terminals, one mtDNA and two rRNA genes for 217 terminals, and barcodes for 576 terminals. Based on congruent topologies, Chauliognathinae is proposed as a sister to the remaining Cantharidae, followed by the redefined Malthininae (including Tytthonyxini), the paraphyletic "dysmorphocerine" lineages (Dysmorphocerinae sensu stricto and Heteromastiginae subfam. nov.), and Silinae + Cantharinae as a terminal clade. The present phylogeny supersedes earlier morphology and short-fragment molecular hypotheses that have not converged on a consensus. Few morphological characters corroborate the DNA-based relationships (see the adults and larval keys). However, morphology-based hypotheses have relied on a few informative characters, and no evidence strongly rejects the preferred molecular topology. The interpretation of morphological characters and uncertain polarity resulting from the high phenotypic disparity of Elateroidea are discussed in detail. The dated phylogeny hypothesizes the earliest split within the Cantharidae in the Berriasian stage (Early Cretaceous, ~141 Myr) and the diversification of most extant subfamilies and tribes already in the Late Cretaceous. The most diverse subfamily, Cantharinae, represents a delayed radiation that started during the Eocene climatic optimum, 55.5 Myr. The late origin of Cantharinae questions the classification of Cretaceous Cantharidae as members of Cantharinae. Instead, the results suggest their deeper rooting after separating from dysmorphocerine lineages and before the node between Cantharinae and Silinae.

Demographic History and Genomic Response to Environmental Changes in a Rapid Radiation of Wild Rats.

For organisms to survive and prosper in a harsh environment, particularly under rapid climate change, poses tremendous challenges. Recent studies have highlighted the continued loss of megafauna in terrestrial ecosystems and the subsequent surge of small mammals, such as rodents, bats, lagomorphs, and insectivores. However, the ecological partitioning of these animals will likely lead to large variation in their responses to environmental change. In the present study, we investigated the evolutionary history and genetic adaptations of white-bellied rats (Niviventer Marshall, 1976), which are widespread in the natural terrestrial ecosystems in Asia but also known as important zoonotic pathogen vectors and transmitters. The southeastern Qinghai-Tibet Plateau was inferred as the origin center of this genus, with parallel diversification in temperate and tropical niches. Demographic history analyses from mitochondrial and nuclear sequences of Niviventer demonstrated population size increases and range expansion for species in Southeast Asia, and habitat generalists elsewhere. Unexpectedly, population increases were seen in N. eha, which inhabits the highest elevation among Niviventer species. Genome scans of nuclear exons revealed that among the congeneric species, N. eha has the largest number of positively selected genes. Protein functions of these genes are mainly related to olfaction, taste, and tumor suppression. Extensive genetic modification presents a major strategy in response to global changes in these alpine species.

Unveiling biogeographical patterns in the worldwide distributed Ceratitis capitata (medfly) using population genomics and microbiome composition.

Invasive species are among the most important, growing threats to food security and agricultural systems. The Mediterranean medfly, Ceratitis capitata, is one of the most damaging representatives of a group of rapidly expanding species in the family Tephritidae, due to their wide host range and high invasiveness potential. Here, we used restriction site-associated DNA sequencing (RADseq) to investigate the population genomic structure and phylogeographical history of medflies collected from six sampling sites, including Africa (South Africa), the Mediterranean (Spain, Greece), Latin America (Guatemala, Brazil) and Australia. A total of 1907 single nucleotide polymorphisms (SNPs) were used to identify two genetic clusters separating native and introduced ranges, consistent with previous findings. In the introduced range, all individuals were assigned to one genetic cluster except for those in Brazil, which showed introgression of an additional genetic cluster that also appeared in South Africa, and which could not be previously identified using microsatellite markers. Moreover, we assessed the microbial composition variations in medfly populations from selected sampling sites using amplicon sequencing of the 16S ribosomal RNA (V4 region). Microbiome composition and structure were highly similar across geographical regions and host plants, and only the Brazilian specimens showed increased diversity levels and a unique composition of its microbiome compared to other sampling sites. The unique SNP patterns and microbiome features in the Brazilian specimens could point to a direct migration route from Africa with subsequent adaptation of the microbiota to the specific conditions present in Brazil. These findings significantly improve our understanding of the evolutionary history of the global medfly invasions and their adaptation to newly colonized environments.

Community assembly and metaphylogeography of soil biodiversity: Insights from haplotype-level community DNA metabarcoding within an oceanic island.

Most of our understanding of island diversity comes from the study of aboveground systems, while the patterns and processes of diversification and community assembly for belowground biotas remain poorly understood. Here, we take advantage of a relatively young and dynamic oceanic island to advance our understanding of ecoevolutionary processes driving community assembly within soil mesofauna. Using whole organism community DNA (wocDNA) metabarcoding and the recently developed metaMATE pipeline, we have generated spatially explicit and reliable haplotype-level DNA sequence data for soil mesofaunal assemblages sampled across the four main habitats within the island of Tenerife. Community ecological and metaphylogeographic analyses have been performed at multiple levels of genetic similarity, from haplotypes to species and supraspecific groupings. Broadly consistent patterns of local-scale species richness across different insular habitats have been found, whereas local insular richness is lower than in continental settings. Our results reveal an important role for niche conservatism as a driver of insular community assembly of soil mesofauna, with only limited evidence for habitat shifts promoting diversification. Furthermore, support is found for a fundamental role of habitat in the assembly of soil mesofauna, where habitat specialism is mainly due to colonization and the establishment of preadapted species. Hierarchical patterns of distance decay at the community level and metaphylogeographical analyses support a pattern of geographic structuring over limited spatial scales, from the level of haplotypes through to species and lineages, as expected for taxa with strong dispersal limitations. Our results demonstrate the potential for wocDNA metabarcoding to advance our understanding of biodiversity.

Mitochondrial Metagenomics Reveals the Ancient Origin and Phylodiversity of Soil Mites and Provides a Phylogeny of the Acari.

High-throughput DNA methods hold great promise for phylogenetic analysis of lineages that are difficult to study with conventional molecular and morphological approaches. The mites (Acari), and in particular the highly diverse soil-dwelling lineages, are among the least known branches of the metazoan Tree-of-Life. We extracted numerous minute mites from soils in an area of mixed forest and grassland in southern Iberia. Selected specimens representing the full morphological diversity were shotgun sequenced in bulk, followed by genome assembly of short reads from the mixture, which produced >100 mitochondrial genomes representing diverse acarine lineages. Phylogenetic analyses in combination with taxonomically limited mitogenomes available publicly resulted in plausible trees defining basal relationships of the Acari. Several critical nodes were supported by ancestral-state reconstructions of mitochondrial gene rearrangements. Molecular calibration placed the minimum age for the common ancestor of the superorder Acariformes, which includes most soil-dwelling mites, to the Cambrian-Ordovician (likely within 455-552 Ma), whereas the origin of the superorder Parasitiformes was placed later in the Carboniferous-Permian. Most family-level taxa within the Acariformes were dated to the Jurassic and Triassic. The ancient origin of Acariformes and the early diversification of major extant lineages linked to the soil are consistent with a pioneering role for mites in building the earliest terrestrial ecosystems.

The limited spatial scale of dispersal in soil arthropods revealed with whole-community haplotype-level metabarcoding.

Soil arthropod communities are highly diverse and critical for ecosystem functioning. However, our knowledge of spatial structure and the underlying processes of community assembly are scarce, hampered by limited empirical data on species diversity and turnover. We implement a high-throughput sequencing approach to generate comparative data for thousands of arthropods at three hierarchical levels: genetic, species and supra-specific lineages. A joint analysis of the spatial arrangement across these levels can reveal the predominant processes driving the variation in biological assemblages at the local scale. This multihierarchical approach was performed using haplotype-level COI metabarcoding of entire communities of mites, springtails and beetles from three Iberian mountain regions. Tens of thousands of specimens were extracted from deep and superficial soil layers and produced comparative phylogeographic data for >1,000 codistributed species and nearly 3,000 haplotypes. Local assemblage composition differed greatly between grasslands and forests and, within each habitat, showed strong spatial structure and high endemicity. Distance decay was high at all levels, even at the scale of a few kilometres or less. The local distance decay patterns were self-similar for the haplotypes and higher hierarchical entities, and this fractal structure was similar in all regions, suggesting that uniform processes of limited dispersal determine local-scale community assembly. Our results from whole-community metabarcoding provide insight into how dispersal limitations constrain mesofauna community structure within local spatial settings over evolutionary timescales. If generalized across wider areas, the high turnover and endemicity in the soil locally may indicate extremely high richness globally, challenging our current estimations of total arthropod diversity on Earth.

Connecting high-throughput biodiversity inventories: Opportunities for a site-based genomic framework for global integration and synthesis.

High-throughput sequencing (HTS) is increasingly being used for the characterization and monitoring of biodiversity. If applied in a structured way, across broad geographical scales, it offers the potential for a much deeper understanding of global biodiversity through the integration of massive quantities of molecular inventory data generated independently at local, regional and global scales. The universality, reliability and efficiency of HTS data can potentially facilitate the seamless linking of data among species assemblages from different sites, at different hierarchical levels of diversity, for any taxonomic group and regardless of prior taxonomic knowledge. However, collective international efforts are required to optimally exploit the potential of site-based HTS data for global integration and synthesis, efforts that at present are limited to the microbial domain. To contribute to the development of an analogous strategy for the nonmicrobial terrestrial domain, an international symposium entitled "Next Generation Biodiversity Monitoring" was held in November 2019 in Nicosia (Cyprus). The symposium brought together evolutionary geneticists, ecologists and biodiversity scientists involved in diverse regional and global initiatives using HTS as a core tool for biodiversity assessment. In this review, we summarize the consensus that emerged from the 3-day symposium. We converged on the opinion that an effective terrestrial Genomic Observatories network for global biodiversity integration and synthesis should be spatially led and strategically united under the umbrella of the metabarcoding approach. Subsequently, we outline an HTS-based strategy to collectively build an integrative framework for site-based biodiversity data generation.

Mimicry diversification in Papilio dardanus via a genomic inversion in the regulatory region of engrailed-invected.

Polymorphic Batesian mimics exhibit multiple protective morphs that each mimic a different noxious model. Here, we study the genomic transitions leading to the evolution of different mimetic wing patterns in the polymorphic Mocker Swallowtail Papilio dardanus. We generated a draft genome (231 Mb over 30 chromosomes) and re-sequenced individuals of three morphs. Genome-wide single nucleotide polymorphism (SNP) analysis revealed elevated linkage disequilibrium and divergence between morphs in the regulatory region of engrailed, a developmental gene previously implicated in the mimicry switch. The diverged region exhibits a discrete chromosomal inversion (of 40 kb) relative to the ancestral orientation that is associated with the cenea morph, but not with the bottom-recessive hippocoonides morph or with non-mimetic allopatric populations. The functional role of this inversion in the expression of the novel phenotype is currently unknown, but by preventing recombination, it allows the stable inheritance of divergent alleles enabling geographic spread and local coexistence of multiple adaptive morphs.

Mitochondrial phylogenomics of the Hymenoptera.

The insect order Hymenoptera presents marvelous morphological and ecological diversity. Higher-level hymenopteran relationships remain controversial, even after recent phylogenomic analyses, as their taxon sampling was limited. To shed light on the origin and diversification of Hymenoptera, in particular the poorly studied Parasitica, we undertook phylogenetic analyses of 40 newly and 43 previously sequenced mitochondrial genomes representing all major clades of Hymenoptera. Various Bayesian inferences using different data partitions and phylogenetic methods recovered similar phylogenetic trees with strong statistical support for almost all nodes. Novel findings of the mitogenomic phylogeny mainly affected the three infraorders Ichneumonomorpha, Proctotrupomorpha and Evaniomorpha, the latter of which was split into three clades. Basal relationships of Parasitica recovered Stephanoidea + (Gasteruptiidae + Aulacidae) as the sister group to Ichneumonomorpha + (Trigonalyoidea + Megalyroidea). This entire clade is sister to Proctotrupomorpha, and Ceraphronoidea + Evaniidae is sister to Aculeata (stinging wasps). Our divergence time analysis indicates that major hymenopteran lineages originated in the Mesozoic. The radiation of early apocritans may have been triggered by the Triassic-Jurassic mass extinction; all extant families were present by the Cretaceous.

Why the COI barcode should be the community DNA metabarcode for the metazoa.

Metabarcoding of complex metazoan communities is increasingly being used to measure biodiversity in terrestrial, freshwater and marine ecosystems, revolutionizing our ability to observe patterns and infer processes regarding the origin and conservation of biodiversity. A fundamentally important question is which genetic marker to amplify, and although the mitochondrial cytochrome oxidase subunit I (COI) gene is one of the more widely used markers in metabarcoding for the Metazoa, doubts have recently been raised about its suitability. We argue that (a) the extensive coverage of reference sequence databases for COI; (b) the variation it presents; (c) the comparative advantages for denoising protein-coding genes; and (d) recent advances in DNA sequencing protocols argue in favour of standardizing for the use of COI for metazoan community samples. We also highlight where research efforts should focus to maximize the utility of metabarcoding.

Metabarcoding of freshwater invertebrates to detect the effects of a pesticide spill.

Biomonitoring underpins the environmental assessment of freshwater ecosystems and guides management and conservation. Current methodology for surveys of (macro)invertebrates uses coarse taxonomic identification where species-level resolution is difficult to obtain. Next-generation sequencing of entire assemblages (metabarcoding) provides a new approach for species detection, but requires further validation. We used metabarcoding of invertebrate assemblages with two fragments of the cox1 "barcode" and partial nuclear ribosomal (SSU) genes, to assess the effects of a pesticide spill in the River Kennet (southern England). Operational taxonomic unit (OTU) recovery was tested under 72 parameters (read denoising, filtering, pair merging and clustering). Similar taxonomic profiles were obtained under a broad range of parameters. The SSU marker recovered Platyhelminthes and Nematoda, missed by cox1, while Rotifera were only amplified with cox1. A reference set was created from all available barcode entries for Arthropoda in the BOLD database and clustered into OTUs. The River Kennet metabarcoding produced matches to 207 of these reference OTUs, five times the number of species recognized with morphological monitoring. The increase was due to the following: greater taxonomic resolution (e.g., splitting a single morphotaxon "Chironomidae" into 55 named OTUs); splitting of Linnaean binomials into multiple molecular OTUs; and the use of a filtration-flotation protocol for extraction of minute specimens (meiofauna). Community analyses revealed strong differences between "impacted" vs. "control" samples, detectable with each gene marker, for each major taxonomic group, and for meio- and macrofaunal samples separately. Thus, highly resolved taxonomic data can be extracted at a fraction of the time and cost of traditional nonmolecular methods, opening new avenues for freshwater invertebrate biodiversity monitoring and molecular ecology.

The contribution of mitochondrial metagenomics to large-scale data mining and phylogenetic analysis of Coleoptera.

A phylogenetic tree at the species level is still far off for highly diverse insect orders, including the Coleoptera, but the taxonomic breadth of public sequence databases is growing. In addition, new types of data may contribute to increasing taxon coverage, such as metagenomic shotgun sequencing for assembly of mitogenomes from bulk specimen samples. The current study explores the application of these techniques for large-scale efforts to build the tree of Coleoptera. We used shotgun data from 17 different ecological and taxonomic datasets (5 unpublished) to assemble a total of 1942 mitogenome contigs of >3000 bp. These sequences were combined into a single dataset together with all mitochondrial data available at GenBank, in addition to nuclear markers widely used in molecular phylogenetics. The resulting matrix of nearly 16,000 species with two or more loci produced trees (RAxML) showing overall congruence with the Linnaean taxonomy at hierarchical levels from suborders to genera. We tested the role of full-length mitogenomes in stabilizing the tree from GenBank data, as mitogenomes might link terminals with non-overlapping gene representation. However, the mitogenome data were only partly useful in this respect, presumably because of the purely automated approach to assembly and gene delimitation, but improvements in future may be possible by using multiple assemblers and manual curation. In conclusion, the combination of data mining and metagenomic sequencing of bulk samples provided the largest phylogenetic tree of Coleoptera to date, which represents a summary of existing phylogenetic knowledge and a defensible tree of great utility, in particular for studies at the intra-familial level, despite some shortcomings for resolving basal nodes.

Genome sequencing of Rhinorhipus Lawrence exposes an early branch of the Coleoptera.

BACKGROUND: Rhinorhipidae Lawrence, 1988 is an enigmatic beetle family represented by a single species, Rhinorhipus tamborinensis Lawrence, 1988, from Australia, with poorly established affinities near the superfamily Elateroidea (click beetles, soldier beetles and fireflies) or the more inclusive series (infraorder) Elateriformia. Its evolutionary position may inform the basal relationships of the suborder Polyphaga, the largest clade of Coleoptera. RESULTS: We analyzed four densely sampled DNA datasets of major coleopteran lineages for mitogenomes, rRNA genes and single copy nuclear genes. Additionally, genome sequencing was used for incorporation of R. tamborinensis into a set of 4220 orthologs for 24 terminals representing 12 polyphagan superfamilies. Topologies differed to various degrees, but all consistently refute the proposed placement of Rhinorhipidae in Elateroidea and instead indicate either sister relationships with other Elateriformia, frequently together with Nosodendridae, another divergent small family hitherto placed in Derodontoidea, or in an isolated position among the deepest lineages of Polyphaga. The phylogenomic analyses recovered Rhinorhipus in a sister position to all other Elateriformia composed of five superfamilies. Therefore, we erect the new superfamily Rhinorhipoidea Lawrence, 1988, stat. Nov., with the type-family Rhinorhipidae. The origins of the Rhinorhipidae were dated to the Upper Triassic/Lower Jurassic at the very early phase of polyphagan diversification. CONCLUSIONS: Thus, Rhinorhipidae adds another example to several recently recognized ancient relict lineages which are interspersed within contemporaneous hugely species-rich lineages of Coleoptera.

The phylogeny of Galerucinae (Coleoptera: Chrysomelidae) and the performance of mitochondrial genomes in phylogenetic inference compared to nuclear rRNA genes.

With efficient sequencing techniques, full mitochondrial genomes are rapidly replacing other widely used markers, such as the nuclear rRNA genes, for phylogenetic analysis but their power to resolve deep levels of the tree remains controversial. We studied phylogenetic relationships of leaf beetles (Chrysomelidae) in the tribes Galerucini and Alticini (root worms and flea beetles) based on full mitochondrial genomes (103 newly sequenced), and compared their performance to the widely sequenced nuclear rRNA genes (full 18S, partial 28S). Our results show that: (i) the mitogenome is phylogenetically informative from subtribe to family level, and the per-nucleotide contribution to nodal support is higher than that of rRNA genes, (ii) the Galerucini and Alticini are reciprocally monophyletic sister groups, if the classification is adjusted to accommodate several 'problematic genera' that do not fit the dichotomy of lineages based on the presence (Alticini) or absence (Galerucini) of the jumping apparatus, and (iii) the phylogenetic results suggest a new classification system of Galerucini with eight subtribes: Oidina, Galerucina, Hylaspina, Metacyclina, Luperina, Aulacophorina, Diabroticina and Monoleptina.