Bulk de novo mitogenome assembly from pooled total DNA elucidates the phylogeny of weevils (Coleoptera: Curculionoidea)

Comparative Analysis of the Mitochondrial Genomes of Chloropidae and Their Implications for the Phylogeny of the Family

Chloropidae, commonly known as grass flies, represent the most taxonomically diverse family of Diptera Carnoidea, comprising over 3000 described species worldwide. Previous phylogenetic studies of this family have predominantly relied on morphological characters, with mitochondrial genomes being reported in a few species. This study presents 11 newly sequenced mitochondrial genomes (10 Chloropidae and 1 Milichiidae) and provides the first comprehensive comparative analysis of mitochondrial genomes for Chloropidae. Apart from 37 standard genes and the control region, three conserved intergenic sequences across Diptera Cyclorrhapha were identified in all available chloropid mitochondrial genomes. Evolutionary rates within Chloropidae exhibit significant variation across subfamilies, with Chloropinae displaying higher rates than the other three subfamilies. Phylogenetic relationships based on mitochondrial genomes were inferred using maximum likelihood and Bayesian methods. The monophyly of Chloropidae and all four subfamilies is consistently strongly supported, while subfamily relationships within Chloropidae remain poorly resolved, possibly due to rapid evolution.

Unusual rearrangements of mitogenomes in Diptera revealed by comparative analysis of 135 tachinid species (Insecta, Diptera, Tachinidae)

The Tachinidae is one of the most speciose families in Diptera, and the exclusively parasitoid species play an important role in regulating populations of many herbivorous insects in ecosystems, including many agricultural pests. To better comprehend the characteristics and evolution of the mitochondrial genome for the Tachinidae, we are adding a massive amount of new molecular data by assembling the mitogenomes for 71 genera and 135 species from all four tachinid subfamilies through next-generation sequencing, and we are presenting the most comprehensive mitogenomic phylogenetic analysis of this family so far. Extensive rearrangements observed in the mitogenome of Admontia podomyia (Exoristinae) are unique for the entire suborder Cyclorrhapha. The rearrangement pattern suggests that the process involved a tandem duplication of the complete mitogenome, followed by both random and nonrandom loss of one copy of each gene. Additionally, five minor mitogenome rearrangements are discovered and described in three subfamilies. We present the largest species-level phylogenetic hypothesis for Tachinidae to date, based on mitogenomes of 152 species of Tachinidae, representing all four subfamilies and with five non-tachinid outgroups. Our analyses support the monophyly of the Tachinidae and most tribes and genera were recovered with good support, but the higher-level phylogenetic relationships within Tachinidae were poorly resolved, indicating that mitogenome data alone are not enough to unambiguously resolve the deeper phylogenetic relationships within Tachinidae.

Mitochondrial genomes provide insights into the Euholognatha (Insecta: Plecoptera)

BACKGROUND

Euholognatha is a monophyletic group within stoneflies comprised by a superfamily Nemouroidea and a family Scopuridae. Based on morphological data, the family-level phylogenetic relationships within Euholognatha are widely accepted, but there is still controversy among different molecular studies. To better understand the phylogeny of all six extant euholognathan families, we sequenced and analyzed seven euholognathan mitogenomes.

RESULTS

The sequence heterogeneity analysis observed a low degree of compositional heterogeneity in euholognathan mitogenomes. Meanwhile, leuctrid mitogenomes were more heterogeneous than other euholognathan families, which may affect the phylogenetic reconstruction. Phylogenetic analyses with various datasets generated three topologies. The Leuctridae was recovered as the earliest branching lineage, and the sister relationship of Capniidae and Taeniopterygidae was supported by most tree topologies and FcLM analyses. When separately excluding sparsely sampled Scopuridae or high heterogeneity leuctrid taxa, phylogenetic analyses under the same methods generated more stable and consistent tree topologies. Finally, based on the results of this study, we reconstructed the relationships within Euholognatha as: Leuctridae + (Scopuridae + ((Taeniopterygidae + Capniidae) + (Nemouridae + Notonemouridae))).

CONCLUSION

Our research shows the potential of data optimizing strategies in reconstructing phylogeny within Euholognatha and provides new insight into the phylogeny of this group.

Relaxed purifying selection pressure drives accelerated and dynamic gene rearrangements in thrips (Insecta: Thysanoptera) mitochondrial genomes

Insect mitochondrial genomes (mitogenome) generally present a typical gene order, which is considered as the ancestral arrangement. All sequenced mitogenomes in the Thysanoptera display high levels of gene rearrangement. Due to limited number of thrips mitogenomes sequenced, how gene rearrangement may be shaped by evolution remain unclear. Here, we analyzed 33 thrips mitogenomes, including 14 newly sequenced. These mitogenomes were diverse in organization, nucleotides substitution and gene arrangements. We found 28 highly rearranged gene orders with the breakpoints of gene rearrangements from 25 to 33. Reconstruction of the ancestors mitochondrial gene arrangements states indicated that Tubulifera have more complex pathways than Terebrantia in the gene order evolution. Molecular calibration estimated that divergence of two suborders occurred in the middle Triassic while the radiation of thrips was associated with the arose and flourish of angiosperm. Our evolutionary hypothesis testing suggests that relaxation of selection pressure enabled the early phase of Thysanoptera evolution, followed by a stronger selective pressure fixed diversification. Our analyses found gene inversion increases the nonsynonymous substitution rates and provide an evolutionary hypothesis driving the diverse gene orders.

In silico prospecting of the mtDNA of Macrobrachium amazonicum from transcriptome data

BACKGROUND

Macrobrachium amazonicum is a freshwater prawn widely distributed in South America that is undergoing speciation, so the denomination "M. amazonicum complex" is used for it. The mitochondrial cytochrome c oxidase subunit I (COI) gene has been used to elucidate this speciation, but heteroplasmies and pseudogenes have been recorded, making separation difficult. Obtaining genes from cDNA (RNA) rather than genomic DNA is an effective tool to mitigate those two types of occurrences. The aim of this study was to assemble in silico the mitochondrial DNA (mtDNA) of the Amazonian coastal population of M. amazonicum inhabiting the state of Pará.

RESULTS

Sequences were obtained from the prawn's transcriptome using the de novo approach. Six libraries of cDNA from the androgen gland, hepatopancreas, and muscle tissue were used. The mtDNA of M. amazonicum was 14,960 bp in length. It contained 13 protein-coding genes, 21 complete transfer RNAs, and the 12S and 16S subunits of ribosomal RNA. All regions were found on the light strand except tRNAGln, which was on the heavy strand. The control region (D-loop) was not recovered, making for a gap of 793 bp. The cladogram showed the formation of the well-defined Macrobrachium clade, with high support value in the established branches (91-100). The three-dimensional spatial conformation of the mtDNA-encoded proteins showed that most of them were mainly composed of major α-helices that typically shows in those proteins inserted in the membrane (mitochondrial).

CONCLUSIONS

It was possible to assemble a large part of the mitochondrial genome of M. amazonicum in silico using data from other genomes deposited in GenBank and to validate it through the similarities between its COI and 16S genes and those from animals of the same region deposited in GenBank. Depositing the M. amazonicum mtDNA sequences in GenBank may help solve the taxonomic problems recorded for the species, in addition to providing complete sequences of candidate coding genes for use as biomarkers in ecological studies.

Novel tRNA Gene Rearrangements in the Mitochondrial Genomes of Poneroid Ants and Phylogenetic Implication of Paraponerinae (Hymenoptera: Formicidae)

Ants (Formicidae) are the most diverse eusocial insects in Hymenoptera, distributed across 17 extant subfamilies grouped into 3 major clades, the Formicoid, Leptanilloid, and Poneroid. While the mitogenomes of Formicoid ants have been well studied, there is a lack of published data on the mitogenomes of Poneroid ants, which requires further characterization. In this study, we first present three complete mitogenomes of Poneroid ants: Paraponera clavata, the only extant species from the subfamily Paraponerinae, and two species (Harpegnathos venator and Buniapone amblyops) from the Ponerinae subfamily. Notable novel gene rearrangements were observed in the new mitogenomes, located in the gene blocks CR-trnM-trnI-trnQ-ND2, COX1-trnK-trnD-ATP8, and ND3-trnA-trnR-trnN-trnS1-trnE-trnF-ND5. We reported the duplication of tRNA genes for the first time in Formicidae. An extra trnQ gene was identified in H. venator. These gene rearrangements could be explained by the tandem duplication/random loss (TDRL) model and the slipped-strand mispairing model. Additionally, one large duplicated region containing tandem repeats was identified in the control region of P. clavata. Phylogenetic analyses based on protein-coding genes and rRNA genes via maximum likelihood and Bayes methods supported the monophyly of the Poneroid clade and the sister group relationship between the subfamilies Paraponerinae and Amblyoponinae. However, caution is advised in interpreting the positions of Paraponerinae due to the potential artifact of long-branch attraction.

Phylogenomics illuminates the phylogeny of flower weevils (Curculioninae) and reveals ten independent origins of brood-site pollination mutualism in true weevils

Weevils are an unusually species-rich group of phytophagous insects for which there is increasing evidence of frequent involvement in brood-site pollination. This study examines phylogenetic patterns in the emergence of brood-site pollination mutualism among one of the most speciose beetle groups, the flower weevils (subfamily Curculioninae). We analysed a novel phylogenomic dataset consisting of 214 nuclear loci for 202 weevil species, with a sampling that mainly includes flower weevils as well as representatives of all major lineages of true weevils (Curculionidae). Our phylogenomic analyses establish a uniquely comprehensive phylogenetic framework for Curculioninae and provide new insights into the relationships among lineages of true weevils. Based on this phylogeny, statistical reconstruction of ancestral character states revealed at least 10 independent origins of brood-site pollination in higher weevils through transitions from ancestral associations with reproductive structures in the larval stage. Broadly, our results illuminate the unexpected frequency with which true weevils-typically specialized phytophages and hence antagonists of plants-have evolved mutualistic interactions of ecological significance that are key to both weevil and plant evolutionary fitness and thus a component of their deeply intertwined macroevolutionary success.

Phylogenomics of weevils revisited: data curation and modelling compositional heterogeneity

Weevils represent one of the most prolific radiations of beetles and the most diverse group of herbivores on land. The phylogeny of weevils (Curculionoidea) has received extensive attention, and a largely satisfactory framework for their interfamilial relationships has been established. However, a recent phylogenomic study of Curculionoidea based on anchored hybrid enrichment (AHE) data yielded an abnormal placement for the family Belidae (strongly supported as sister to Nemonychidae + Anthribidae). Here we reanalyse the genome-scale AHE data for Curculionoidea using various models of molecular evolution and data filtering methods to mitigate anticipated systematic errors and reduce compositional heterogeneity. When analysed with the infinite mixture model CAT-GTR or using appropriately filtered datasets, Belidae are always recovered as sister to the clade (Attelabidae, (Caridae, (Brentidae, Curculionidae))), which is congruent with studies based on morphology and other sources of molecular data. Although the relationships of the 'higher Curculionidae' remain challenging to resolve, we provide a consistent and robust backbone phylogeny of weevils. Our extensive analyses emphasize the significance of data curation and modelling across-site compositional heterogeneity in phylogenomic studies.

Climatic oscillation promoted diversification of spinous assassin bugs during Pleistocene glaciation

Insect speciation is among the most fascinating topics in evolutionary biology; however, its underlying mechanisms remain unclear. Allopatric speciation represents one of the major types of speciation and is believed to have frequently occurred during glaciation periods, when climatic oscillation may have caused suitable habitats to be fragmented repeatedly, creating geographical isolation among populations. However, supporting evidence for allopatric speciation of insects in East Asia during the Pleistocene glaciation remains lacking. We aim to investigate the effect of climatic oscillation during the Pleistocene glaciation on the diversification pattern and evolutionary history of hemipteran insects and to test the hypothesis of Pleistocene species stability using spinous assassin bugs Sclomina (Hemiptera: Reduviidae), a small genus widely distributed in southern China but was later found to have cryptic species diversity. Here, using the whole mitochondrial genome (mitogenome) and nuclear ribosomal RNA genes, we investigated both interspecific and intraspecific diversification patterns of spinous assassin bugs. Approximate Bayesian computation, ecological niche modeling, and demographic history analyses were also applied to understand the diversification process and driven factors. Our data suggest that the five species of Sclomina are highly diverged, despite three of them currently being cryptic. Speciation occurred during the Pleistocene when suitable distribution areas were possibly fragmented. Six phylogeographic groups in the type species S. erinacea were identified, among which two groups underwent expansion during the early Last Glacial Period and after Last Glacier Maximum. Our analyses suggest that this genus may have experienced climate-driven habitat fragmentation and postglacial expansion in the Pleistocene, promoting allopatric speciation and intraspecific diversification. Our results reveal underestimated species diversity in a small insect group and illustrate a remarkable example of allopatric speciation of insects in East Asia promoted by Pleistocene climatic oscillations. These findings provide important insights into the speciation processes and aid the conservation of insect species diversity.

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 allow 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.

Characterization of the complete mitochondrial genome of the brazilian cownose ray Rhinoptera brasiliensis (Myliobatiformes, Rhinopteridae) in the western Atlantic and its phylogenetic implications

BACKGROUND

The Brazilian cownose ray, Rhinoptera brasiliensis has undergone a global population reduction and is currently classified by IUCN as Vulnerable. This species is sometimes confused with Rhinoptera bonasus, the only external diagnostic characteristic to distinguish between both species is the number of rows of tooth plates. Both cownose rays overlap geographically from Rio de Janeiro to the western North Atlantic. This calls for a more comprehensive phylogenetic assessment using mitochondria DNA genomes to better understand the relationships and delimitation of these two species.

METHODS AND RESULTS

The mitochondrial genome sequences of R. brasiliensis was obtained by next-generation sequencing. The length of the mitochondrial genome was 17,759 bp containing 13 protein-coding genes (PCGs), two ribosomal RNA (rRNA) genes, 22 transfer RNA (tRNA) genes, and a non-coding control region (D-loop). Each PCG was initiated by an authoritative ATG codon, except for COX1 initiated by a GTG codon. Most of the PCGs were terminated by a complete codon (TAA/TAG), while an incomplete termination codon (TA/T) was found in five out of the 13 PCGs. The phylogenetic analysis showed that R. brasiliensis was closely related to R. steindachneri whereas the reported mitogenome as R. steindachneri (GenBank accession number KM364982), differs from multiple mitocondrial DNA sequences of R. steindachneri and is nearly identical to that of R. javanica.

CONCLUSION

The new mitogenome determined in this study provides new insight into the phylogenetic relationships in Rhinoptera, while providing new molecular data that can be applied to population genetic studies.

Complete Mitochondrial Genome of Scolytoplatypodini Species (Coleoptera: Curculionidae: Scolytinae) and Phylogenetic Implications

The complete mitochondrial genomes (mitogenomes) of beetles in the tribe Scolytoplatypodini (genus Scolytoplatypus) were sequenced and annotated. These included Scolytoplatypus raja (15,324 bp), Scolytoplatypus sinensis (15,394 bp), Scolytoplatypus skyliuae (15,167 bp), and Scolytoplatypus wugongshanensis (15,267 bp). The four mitogenomes contained 37 typical genes, including 13 protein-coding genes (PCGs), 22 transfer RNA genes (tRNAs), and 2 ribosomal RNA genes (rRNAs). The gene orientation and arrangement of the four mitogenomes were similar to other Coleoptera mitogenomes. PCGs mostly started with ATN and terminated with TAA. The Ka/Ks ratio of 13 PCGs in the four species revealed that cox1 had the slowest evolutionary rate and atp8 and nad6 had a higher evolutionary rate. All tRNAs had typical cloverleaf secondary structures, but trnS1 lacked dihydrouridine arm. Partial tRNAs lost the discriminator nucleotide. The trnY did not possess the discriminator nucleotide and also lost three bases, showing a special amino-acyl arm. Bayesian inference (BI) and maximum likelihood (ML) methods were conducted for phylogenetic analyses using 13 PCGs. Scolytoplatypodini was clustered with Hylurgini and Hylastini, and the monophyly of Scolytoplatypodini was supported. The four newly sequenced mitogenomes increase understanding of the evolutionary relationships of Scolytoplatypodini and other Scolytinae species.

Fossil History of Ambrosia Beetles (Coleoptera; Platypodidae) with Description of a New Genus from Dominican Amber

The fossil record of ambrosia beetles is summarized and a new genus and species in the subfamily Mecopelminae, Gongyloceria dominicana gen. et sp. nov. is described. The new genus differs from the extant genus Mecopelmus Blackman, 1944, possessing weakly elongated tarsi, finely faceted eyes, a scape not reaching the posterior margin of the eye, coarsely sculptured pronotum and elytra, and a larger body size. It is the first Miocene record of the Mecopelminae. The genus Xyleborites Wickham, 1913, placem. n. is transferred from the Scolytidae to the Platypodidae. A list of the fossil Platypodidae, including a key to the subfamilies and tribes of ambrosia beetles, is presented. Xyleborites longipennis Wickham, 1913 and Gongyloceria dominicana sp. nov. may have been related to Paullinieae recorded from the late Eocene and the early Miocene of North and Central America. Distribution maps with fossil records for ambrosia beetles are included.

Mitochondrial genomes of Sternochetus species (Coleoptera: Curculionidae) and the phylogenetic implications

The three weevil species, Sternochetus gravis, S. mangiferae, and S. olivieri, have all been reported to be serious pests of mango fruits. Morphology, biology, and various management approaches of these economically important weevils have been well studied. However, no mitochondrial genomes have been reported from the genus Sternochetus. Herein, we assembled mitogenomes of all the three Sternochetus species to reveal their mitogenomic characteristics. A DNA library of 350 bp insert size was constructed and sequenced in Illumina's HiSeq 6000 platform with a pair-end 150 bp sequencing strategy by Novogene. The sequence reads were assembled using GetOrganelle v1.7.1 and the genes were annotated by Geneious Prime 2021.0.3 and MITOS Web Server. Coupled with 61 published mitogenomes from 13 subfamilies of Curculionidae, we reconstructed phylogenetic trees to resolve evolutionary relationships of these closely related species and also examined subfamily-level classification among Curculionidae. All three mitogenomes are double-stranded circular molecules with 22 transfer RNA genes, 13 protein-coding genes (PCGs), 2 ribosomal RNA genes, and 1 noncoding control region as in other insects. Higher interspecific nucleotide divergence (about 10%) of 13 PCGs indicated these three Sternochetus species diverged a long time ago. Phylogenetic analyses using both maximum likelihood and Bayesian inference methods showed that Sternochetus falls into the basal clade of Cryptorhynchini, a tribe in the subfamily Molytinae. The relationship of S. olivieri as a sister species to S. gravis + S. mangiferae was strongly supported. The monophyly of Cryptorhynchini was also well supported whereas Molytinae was suggested to be a polyphyletic group.

Mitogenomics of the Olive Seed Weevil, Anchonocranus oleae Marshall and Implications for Its Phylogenetic Position in Curculionidae

Simple Summary

Anchonocranus oleae is a southern African weevil that feeds on the seeds of the African Wild Olive, a close relative of the European cultivated olive tree. The species is known to occur in the Western Cape of South Africa, the main region of olive production in Southern Africa. We generated reference DNA barcodes and the complete mitogenome of A. oleae as part of our ongoing genetic cataloguing of insects associated with wild and cultivated olives in South Africa. The phylogenetic position of A. oleae in the family Curculionidae was inferred to be in the Curculioninae, Conoderinae, Cossoninae, Molytinae, and Scolytinae (CCCMS) clade but could not be precisely determined due to the paucity of genetic data for adequate taxonomic context, highlighting the need for further coverage of related tribes and genera. Nevertheless, the data generated in this study contribute to the enrichment of baseline information on olive-associated insects, in general, and on the genus Anchonocranus, in particular.

Abstract

Anchonocranus oleae Marshall (Coleoptera: Curculionidae) is a seed-feeding weevil native to southern Africa; its larvae are known to develop in the fruits of the African Wild Olive and, more rarely, cultivated olives. The species has been mainly found in the Western Cape province of South Africa, but it has remained in relative obscurity because it does not seem to represent a current threat to commercial olive production. As part of an ongoing effort to produce baseline genetic data for olive-associated entomofauna in South Africa, we generated reference DNA barcodes for A. oleae collected from wild and cultivated olives and sequenced its mitogenome for assessment of the phylogenetic position of the species in the family Curculionidae. The mitochondrial phylogeny estimate indicated that A. oleae shares a common ancestor with Elaidobius (tribe Derelomini), but a definite and close relationship to this tribe and the precise tribal placement of A. oleae in the subfamily Curculioninae could not be inferred due to the lack of representative mitogenomes of other relevant curculionine tribes and genera. This study will assist future work on the DNA-based species identification, genetic diversity, and phylogenetic position of the genus Anchonocranus and related taxa.

The mitochondrial genome of Gymnopternus bomiensis (Diptera: Dolichopodidae)

The mitochondrial genome of Gymnopternus bomiensis (Yang, 1996) (Diptera: Dolichopodidae) has been reported in this study. This is the first mitogenome representative of Gymnopternus. The sequenced region is determined to be 15,212 bp, including 13 protein-coding genes, 2 ribosomal RNAs, 22 transfer RNAs, and a partial A-T rich region. The nucleotide composition biased toward A and T, and the overall A + T% was up to 74.0％. Additionally, we reconstructed the phylogeny of relative species using 13 PCRs and two rRNAs. Bayesian inference analysis strongly supported the monophyly of Dolichopodinae. It also suggested that Gymnopternus is the sister group of Dolichopus.

Positive Correlation of the Gene Rearrangements and Evolutionary Rates in the Mitochondrial Genomes of Thrips (Insecta: Thysanoptera)

Simple Summary

Aeolothrips, commonly known as banded thrips, is the largest genus of the family Aeolothripidae (predatory thrips). In the current study, we sequenced the mitochondrial genome (mitogenome) of the banded thrip species Aeolothrips xinjiangensis. We found a novel gene arrangement in this mitogenome that has not been reported in Thysanoptera. By comparing the gene order and rearrangement patterns, we found seven identical gene blocks and three identical rearrangement events in two mitogenomes of banded thrips. There was marked variation in the mitochondrial gene order across thrip species, with only two conserved gene blocks shared by all 14 thrips. In addition, we found a positive correlation between the degree of gene rearrangement and evolutionary rate. Our results suggested that the mitogenomes of thrips have tended to be stable since their massive rearrangement.

Abstract

Extensive gene rearrangement is characteristic in the mitogenomes of thrips (Thysanoptera), but the historical process giving rise to the contemporary gene rearrangement pattern remains unclear. To better understand the evolutionary processes of gene rearrangement in the mitogenomes of thrips, we sequenced the mitogenome of the banded thrip species Aeolothrips xinjiangensis. First, we found a novel mitochondrial gene order in this species. This mitogenome is 16,947 bp in length and encodes the typical 37 coding genes (13 protein-coding genes, 22 tRNA genes, and two rRNA genes) of insects. The gene arrangement was dramatically different from the putative ancestral mitogenome, with 26 genes being translocated, eight of which were inverted. Moreover, we found a novel, conserved gene block, trnC-trnY, which has not been previously reported in the mitogenomes of thrips. With this newly assembled mitogenome, we compared mitogenome sequences across Thysanoptera to assess the evolutionary processes giving rise to the current gene rearrangement pattern in thrips. Seven identical gene blocks were shared by two sequenced banded thrip mitogenomes, while the reversal of ND2 combined with TDRL events resulted in the different gene orders of these two species. In phylogenetic analysis, the monophyly of the suborders and families of Thysanoptera was well supported. Across the gene orders of 14 thrips, only two conserved gene blocks, ATP8-ATP6 and ND4-ND4L, could be found. Correlation analysis showed that the degree of gene rearrangement was positively correlated with the non-synonymous substitution rate in thrips. Our study suggests that the mitogenomes of thrips remain stable over long evolutionary timescales after massive rearrangement during early diversification.

A global molecular phylogeny yields insights into the dispersal and invasion history of Junonia, a butterfly genus with remarkable dispersal abilities

The nymphalid butterfly genus Junonia has remarkable dispersal abilities. Occurring on every continent except Europe and Antarctica, Junonia are often among the only butterflies on remote oceanic islands. The biogeography of Junonia has been controversial, plagued by taxonomic disputes, small phylogenetic datasets, incomplete taxon sampling, and shared interspecific mitochondrial haplotypes. Junonia originated in Africa but its route into the New World remains unknown. Presented here is, to our knowledge, the most comprehensive Junonia phylogeny to date, using full mitogenomes and nuclear ribosomal RNA repeats from 40 of 47 described species. Junonia is monophyletic and the genus Salamis is its probable sister clade. Genetic exchange between Indo-Pacific Junonia villida and New World Junonia vestina is evident, suggesting a trans-Pacific route into the New World. However, in both phylogenies, the sister clades to most New World Junonia contain both African and Asian species. Multiple trans-Atlantic or trans-Pacificinvasions could have contributed to New World diversification. Hybridization and lateral transfer of mitogenomes, already well-documented in New World Junonia, also occurs in at least two Old World lineages (Junonia orithya/Junonia hierta and Junonia iphita/Junonia hedonia). Variation associated with reticulate evolution creates challenges for phylogenetic reconstruction, but also may have contributed to patterns of speciation and diversification in this genus.

Metabarcoding versus mapping unassembled shotgun reads for identification of prey consumed by arthropod epigeal predators

Background

A central challenge of DNA gut content analysis is to identify prey in a highly degraded DNA community. In this study, we evaluated prey detection using metabarcoding and a method of mapping unassembled shotgun reads (Lazaro).

Results

In a mock prey community, metabarcoding did not detect any prey, probably owing to primer choice and/or preferential predator DNA amplification, while Lazaro detected prey with accuracy 43–71%. Gut content analysis of field-collected arthropod epigeal predators (3 ants, 1 dermapteran, and 1 carabid) from agricultural habitats in Brazil (27 samples, 46–273 individuals per sample) revealed that 64% of the prey species detections by either method were not confirmed by melting curve analysis and 87% of the true prey were detected in common. We hypothesized that Lazaro would detect fewer true- and false-positive and more false-negative prey with greater taxonomic resolution than metabarcoding but found that the methods were similar in sensitivity, specificity, false discovery rate, false omission rate, and accuracy. There was a positive correlation between the relative prey DNA concentration in the samples and the number of prey reads detected by Lazaro, while this was inconsistent for metabarcoding.

Conclusions

Metabarcoding and Lazaro had similar, but partially complementary, detection of prey in arthropod predator guts. However, while Lazaro was almost 2× more expensive, the number of reads was related to the amount of prey DNA, suggesting that Lazaro may provide quantitative prey information while metabarcoding did not.

Higher-level phylogeny of Chrysomelidae based on expanded sampling of mitogenomes

Chrysomelidae is one of the most diverse lineages of beetles. The classification and phylogeny of Chrysomelidae have been contentious. In this study, we obtained 16 new mitogenome sequences by using next-generation sequencing. Combined with the published mitogenomes, we inferred the phylogenetic relationships of Chrysomelidae. Different data recoding strategies and substitution models were applied to phylogenetic reconstruction. In the Maximum likelihood analyses under the homogeneous model, Dayhoff recoding allowed for the improved phylogenetic resolution due to the decreased level of heterogeneous sequence divergence. Bayesian inference under the heterogeneous model yielded generally well resolved subfamily relationships. The present mitogenome data strongly supported Chrysomelidae as a monophyletic group. Consistent with previous work, we found three major subfamily clades within Chrysomelidae. However, the pattern of the "sagrine" clade plus the "eumolpine" clade being sister to the "chrysomeline" clade contrasted with the prior study. The placement of the genus Syneta with regards to these three clades was ambiguous. Relationships recovered suggested several major chrysomelid lineages, including: (1) Bruchinae in the "sagrine" clade; (2) Donaciinae + Criocerinae; (3) Spilopyrinae + (Cassidinae + (Eumolpinae + (Lamprosomatinae + Cryptocephalinae))); (4) Chrysomelinae + (Alticinae + Galerucinae). Results also suggested the placement of Timarcha outside the major Chrysomelinae.

Development and evaluation of a meat mitochondrial metagenomic (3MG) method for composition determination of meat from fifteen mammalian and avian species

BACKGROUND

Bioassessment and biomonitoring of meat products are aimed at identifying and quantifying adulterants and contaminants, such as meat from unexpected sources and microbes. Several methods for determining the biological composition of mixed samples have been used, including metabarcoding, metagenomics and mitochondrial metagenomics. In this study, we aimed to develop a method based on next-generation DNA sequencing to estimate samples that might contain meat from 15 mammalian and avian species that are commonly related to meat bioassessment and biomonitoring.

RESULTS

In this project, we found the meat composition from 15 species could not be identified with the metabarcoding approach because of the lack of universal primers or insufficient discrimination power. Consequently, we developed and evaluated a meat mitochondrial metagenomics (3MG) method. The 3MG method has four steps: (1) extraction of sequencing reads from mitochondrial genomes (mitogenomes); (2) assembly of mitogenomes; (3) mapping of mitochondrial reads to the assembled mitogenomes; and (4) biomass estimation based on the number of uniquely mapped reads. The method was implemented in a python script called 3MG. The analysis of simulated datasets showed that the method can determine contaminant composition at a proportion of 2% and the relative error was < 5%. To evaluate the performance of 3MG, we constructed and analysed mixed samples derived from 15 animal species in equal mass. Then, we constructed and analysed mixed samples derived from two animal species (pork and chicken) in different ratios. DNAs were extracted and used in constructing 21 libraries for next-generation sequencing. The analysis of the 15 species mix with the method showed the successful identification of 12 of the 15 (80%) animal species tested. The analysis of the mixed samples of the two species revealed correlation coefficients of 0.98 for pork and 0.98 for chicken between the number of uniquely mapped reads and the mass proportion.

CONCLUSION

To the best of our knowledge, this study is the first to demonstrate the potential of the non-targeted 3MG method as a tool for accurately estimating biomass in meat mix samples. The method has potential broad applications in meat product safety.

The mitochondrial genome of Calohilara tibetensis Ding, He, Lin, and Yang, 2020 (Diptera: Empididae)

The dance fly Calohilara tibetensis Ding, He, Lin and Yang, 2020 belongs to the subfamily Empidinae of Empididae. The mitochondrial genome of C. tibetensis was sequenced as the new representative of the subfamily Empidinae. The nearly complete mitogenome was 15,354 bp, consisting of 13 protein-coding genes, two rRNAs, and 22 tRNAs. All genes have the similar locations and strands with that of other published species of Empididae. The nucleotide composition biases toward A and T is 77.4% of the entirety. All PCGs start with ATN codons except COI and NAD1, and end with TAA or incomplete stop codon. Bayesian inference (BI) analysis strongly supported the monophyly of both Empididae and Dolichopodidae and the monophyly of subfamily Empidinae. It suggested that Calohilara is the sister group of Hilara.

Comparative Mitogenomics and Phylogenetic Analyses of Pentatomoidea (Hemiptera: Heteroptera)

Pentatomoidea is the largest superfamily of Pentatomomorpha; however, the phylogenetic relationships among pentatomoid families have been debated for a long time. In the present study, we gathered the mitogenomes of 55 species from eight common families (Acanthosomatidae, Cydnidae, Dinidoridae, Scutelleridae, Tessaratomidae, Plataspidae, Urostylididae and Pentatomidae), including 20 newly sequenced mitogenomes, and conducted comparative mitogenomic studies with an emphasis on the structures of non-coding regions. Heterogeneity in the base composition, and contrasting evolutionary rates were encountered among the mitogenomes in Pentatomoidea, especially in Urostylididae, which may lead to unstable phylogenetic topologies. When the family Urostylididae is excluded in taxa sampling or the third codon positions of protein coding genes are removed, phylogenetic analyses under site-homogenous models could provide more stable tree topologies. However, the relationships between families remained the same in all PhyloBayes analyses under the site-heterogeneous mixture model CAT + GTR with different datasets and were recovered as (Cydnidae + (((Tessaratomidae + Dinidoridae) + (Plataspidae + Scutelleridae)) + ((Acanthosomatidae + Urostylididae) + Pentatomidae)))). Our study showed that data optimizing strategies after heterogeneity assessments based on denser sampling and the use of site-heterogeneous mixture models are essential for further analysis of the phylogenetic relationships of Pentatomoidea.

Middle Jurassic origin in India: a new look at evolution of Vermileonidae and time-scaled relationships of lower brachyceran flies

Vermileonidae (wormlions) comprises 61 described species in 12 genera and is one of the oddest and most rarely collected dipteran families. Larvae of Vermileonidae are famous for their pitfall prey behaviour, whereas the adults are rarely seen. Here we report, for the first time, the complete mitochondrial genome from members of Vermileonidae. Phylogenetic reconstruction based on a representative sampling of the order reveals new insights into relationships between the Vermileonidae and other members of lower brachyceran flies. A sister-group relationship between Vermileonidae and Xylophagidae is supported, and the higher-level clade relationships are Tabanomorpha + (Muscomorpha + (Xylophagomorpha + Stratiomyomorpha)). Combining mitochondrial genome data with a morphological phylogeny, geographical distribution and geological history, we propose that the Vermileonidae originated in India during the Middle Jurassic, spreading to Africa via land bridges during the Late Jurassic, to the Palaearctic after the collision of India with Laurasia in the Late Oligocene, and to the Nearctic in the Early Miocene, via either North Atlantic or Beringian land bridges. Wet forested regions have proved to be barriers to their dispersal, explaining their absence from Central Africa, South America and eastern North America.

Mitogenome of Coprophanaeus ensifer and phylogenetic analysis of the Scarabaeidae family (Coleoptera)

Several studies about the phylogenetic relationships of the Scarabaeinae subfamily (Coleoptera: Scarabaeidae) have been performed, but some phylogenetic uncertainties persist including the relationship and monophyly of different tribes and some genera. The aim of this study was to characterize the mitogenome of Coprophanaeus ensifer in order to establish its position within the Scarabaeidae family and to contribute to the resolution of some phylogenetic uncertainties. The mitogenome was sequenced on the Illumina HiSeq 4000, assembled using the Mitobim software and annotated in MITOS WebServer. The phylogenetic trees were reconstructed by Bayesian inference. The C. ensifer mitogenome is a molecule of 14,964 bp that contains the number and organization of the genes similar to those of most Coleoptera species. Phylogenetic reconstruction suggests monophyly of the tribe Phanaeini and supports the hypothesis that Coprini is a sister group of Phanaeini. The results also revealed the position of the tribe Oniticellini which is grouped with Onthophagini and Onitini. The geographic distribution of these species that form the most ancestral clade suggests with Scarabaeinae originated in Africa.

The mitochondrial genome of Lichtwardtia dentalis Zhang, Masunaga et Yang, 2009 (Diptera: Dolichopodidae)

The long-legged fly Lichtwardtia dentalis Zhang, Masunaga et Yang, 2009 belongs to the subfamily Dolichopodinae of Dolichopodidae. The newly sequenced mitogenome of L. dentalis is a new representative of the subfamily. The nearly complete mitogenome is 15,124 bp in length, consisting of 13 protein-coding genes (PCGs), 2 ribosomal RNA genes (rRNAs), and 22 transfer RNA genes (tRNAs). All genes have similar locations and encoding directions with that of other published mitogenomes of Dolichopodidae. The nucleotide composition biases toward A and T with the overall A + T % is 73.9%. All protein-coding genes initiate with standard start codon ATN except COI and NAD1, and TAA/TAG are conventionally used as stop codons. All tRNAs, ranging from 62 to 71 bp, have a clover-leaf structure. Based on the result of the phylogenetic analysis, Dolichopodidae and Empididae were monophyletic, and the relationships among subfamilies of Dolichopodidae were Diaphorinae + (Peloropeodinae + (Xanthochlorinae + (Medeterinae + Dolichopodinae))). The monophyly of the subfamily Dolichopodinae and the sister relationship between Dolichopus and Lichtwardtia were also strongly supported.

Origin and evolution of fungus farming in wood-boring Coleoptera - a palaeontological perspective

Insect-fungus mutualism is one of the better-studied symbiotic interactions in nature. Ambrosia fungi are an ecological assemblage of unrelated fungi that are cultivated by ambrosia beetles in their galleries as obligate food for larvae. Despite recently increased research interest, it remains unclear which ecological factors facilitated the origin of fungus farming, and how it transformed into a symbiotic relationship with obligate dependency. It is clear from phylogenetic analyses that this symbiosis evolved independently many times in several beetle and fungus lineages. However, there is a mismatch between palaeontological and phylogenetic data. Herein we review, for the first time, the ambrosia system from a palaeontological perspective. Although largely ignored, families such as Lymexylidae and Bostrichidae should be included in the list of ambrosia beetles because some of their species cultivate ambrosia fungi. The estimated origin for some groups of ambrosia fungi during the Cretaceous concurs with a known high diversity of Lymexylidae and Bostrichidae at that time. Although potentially older, the greatest radiation of various ambrosia beetle lineages occurred in the weevil subfamilies Scolytinae and Platypodinae during the Eocene. In this review we explore the evolutionary relationship between ambrosia beetles, fungi and their host trees, which is likely to have persisted for longer than previously supposed.

The mitochondrial genome of Desmomyia sinensis (Diptera: Rhagionidae)

The mitochondrial genome of Desmomyia sinensis Yang et Yang, 1997 was sequenced, the new representative of the family Rhagionidae. The complete mitochondrial genome is 16,430 bp totally, which consists of 13 protein-coding genes, two ribosomal RNAs, 22 transfer RNAs, and one non-coding control region. The nucleotide composition is biased toward A and T, accounting for 77.9% of the total. All genes have the conservational arrangement like other published species of Rhagionidae. Bayesian inference analysis strongly supported the monophyly of Rhagionidae and sister relationship between Coenomyiidae and Xylophagidae. The phylogenetic relationship among families of Tabanomorpha is very clear: (Athericidae + Tabanidae) + (Rhagionidae + (Coenomyiinae + Xylophagidae)).

Unique Duplication of trnN in Odontoptilum angulatum (Lepidoptera: Pyrginae) and Phylogeny within Hesperiidae

To explore the variation and relationship between gene rearrangement and phylogenetic effectiveness of mitogenomes among lineages of the diversification of the tribe Tagiadini in the subfamily Pyrginae, we sequenced the complete mitogenome of Odontoptilum angulatum. The genome is 15,361 bp with the typical 37 genes, a large AT-rich region and an additional trnN (trnN2), which is completely identical to trnN (sequence similarity: 100%). The gene order differs from the typical Lepidoptera-specific arrangement and is unique to Hesperiidae. The presence of a "pseudo-trnS1" in the non-coding region between trnN1 and trnN2 supports the hypothesis that the presence of an extra trnN can be explained by the tandem duplication-random loss (TDRL) model. Regarding the phylogenetic analyses, we found that the dataset comprising all 37 genes produced the highest node support, as well as a monophyly of Pyrginae, indicating that the inclusion of RNAs improves the phylogenetic signal. Relationships among the subfamilies in Hesperiidae were also in general agreement with the results of previous studies. The monophyly of Tagiadini is strongly supported. Our study provides a new orientation for application of compositional and mutational biases of mitogenomes in phylogenetic analysis of Tagiadini and even all Hesperiidae based on larger taxon sampling in the future.

Global phylogeography and invasion history of the spotted lanternfly revealed by mitochondrial phylogenomics

Biological invasion has been a serious global threat due to increasing international trade and population movements. Tracking the source and route of invasive species and evaluating the genetic differences in their native regions have great significance for the effective monitoring and management, and further resolving the invasive mechanism. The spotted lanternfly Lycorma delicatula is native to China and invaded South Korea, Japan, and the United States during the last decade, causing severe damages to the fruits and timber industries. However, its global phylogeographic pattern and invasion history are not clearly understood. We applied high-throughput sequencing to obtain 392 whole mitochondrial genome sequences from four countries to ascertain the origin, dispersal, and invasion history of the spotted lanternfly. Phylogenomic analyses revealed that the spotted lanternfly originated from southwestern China, diverged into six phylogeographic lineages, and experienced northward expansion across the Yangtze River in the late Pleistocene. South Korea populations were derived from multiple invasions from eastern China and Japan with two different genetic sources of northwestern (Loess Plateau) and eastern (East Plain) lineages in China, whereas the each of Japan and the United States had only one. The United States populations originated through single invasive event from South Korea, which served as a bridgehead of invasion. The environmental conditions, especially the distribution of host Ailanthus trees, and adaptability possibly account for the rapid spread of the spotted lanternfly in the native and introduced regions.

The complete mitochondrial genome of Dicraeus orientalis Becker, 1911 (Diptera: Chloropidae)

Dicraeus orientalis feeds on the seeds of Poaceae. The complete mitochondrial genome of D. orientalis was sequenced and annotated as the first representative of the family Chloropidae. The full length of mitogenome was 16,188 bp, consisting of 13 protein-coding genes (PCGs), 22 transfer RNA genes (tRNAs), and two ribosomal RNA genes (rRNAs). The nucleotide composition was highly A + T biased, accounting for 79.1% of the whole mitogenome. All PCGs start with ATN codons except COI, which end with TAN or incomplete stop codon. ML analysis revealed that Carnoidea was closely related to Ephydroidea and the phylogenetic relationship within Acalyptratae was Tephritoidea + ((Carnoidea + Ephydroidea) + Opomyzoidea).

First phylogenetic analysis of Dryophthorinae (Coleoptera, Curculionidae) based on structural alignment of ribosomal DNA reveals Cenozoic diversification

Dryophthorinae is an economically important, ecologically distinct, and ubiquitous monophyletic group of pantropical weevils with more than 1,200 species in 153 genera. This study provides the first comprehensive phylogeny of the group with the aim to provide insights into the process and timing of diversification of phytophagous insects, inform classification and facilitate predictions. The taxon sampling is the most extensive to date and includes representatives of all five dryophthorine tribes and all but one subtribe. The phylogeny is based on secondary structural alignment of 18S and 28S rRNA totaling 3,764 nucleotides analyzed under Bayesian and maximum likelihood inference. We used a fossil-calibrated relaxed clock model with two approaches, node-dating and fossilized birth-death models, to estimate divergence times for the subfamily. All tribes except the species-rich Rhynchophorini were found to be monophyletic, but higher support is required to ascertain the paraphyly of Rhynchophorini with more confidence. Nephius is closely related to Dryophthorini and Stromboscerini, and there is strong evidence for paraphyly of Sphenophorina. We find a large gap between the divergence of Dryophthorinae from their sister group Platypodinae in the Jurassic-Cretaceous boundary and the diversification of extant species in the Cenozoic, highlighting the role of coevolution with angiosperms in this group.

Higher-level phylogenetic relationships of rove beetles (Coleoptera, Staphylinidae) inferred from mitochondrial genome sequences

Rove beetles (Staphylinidae) and allied families constitute a huge radiation of Coleoptera, but basal relationships in this group remain controversial. In this study, we newly sequenced eight mitogenomes of representatives of Staphylinidae by using next-generation sequencing method. Together with 99 existing mitogenomes of Staphyliniformia, (sub)family relationships were investigated with ML and Bayesian searches under various substitution models and data recoding schemes. The results consistently supported Scydmaenidae and Silphidae to be subordinate groups of Staphylinidae. Within the monophyletic Staphylinidae (including Scydmaenidae and Silphidae), the hypothesis of four major subfamily groups cannot be confirmed. Bayesian inferences under the site-heterogeneous mixture model generally supported the basal position of major clades corresponding to the Omaliine group. At the subfamily level, the monophyly of Pselaphinae, Oxytelinae, Scaphidiinae, Steninae and Staphylininae was supported. However, the subfamilies Omaliinae, Tachyporinae, Aleocharinae and Paederinae were each non-monophyletic.

Tanyptera (Tanyptera) hebeiensis Yang et Yang (Diptera: Tipulidae) newly recorded from Shandong, China: sequencing and phylogenetic analysis of the mitochondrial genome

The genus Tanyptera Latreille, 1804 is recorded from Shandong Province, China for the first time with T. (T.) hebeiensis Yang et Yang, 1988 found in Mount Kunyu, Shandong. In this study, we report the complete mitochondrial genome sequence of T. (T.) hebeiensis, representing the first mitochondrial genome of the subfamily Ctenophorinae (Diptera: Tipulidae), which is a circular molecule of 15,888 bp with an AT content of 77.6%. The mitochondrial genome contains 13 protein-coding genes (PCGs), 22 transfer RNA genes (tRNAs), two ribosomal RNA genes (rRNAs), and a non-coding region. Gene overlaps are found at nine gene junctions, ranging from 1 to 8 bp in length. The canonical mitochondrial start codons for invertebrate mitochondrial genomes are found in 12 PCGs, except for COI which uses the uncanonical start codons TCG. Stop codons of 10 PCGs are invariably complete TAA and TAG, while COII, ND4, and ND5 end with a single thymine stop codon. Phylogenetic analysis reveals that the Pediciidae is a sister group to the remaining Tipuloidea, the Cylindrotomidae has a sister-group relationship with the Tipulidae, and the Limoniidae is not a monophyletic clade.

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.

New Light on Historical Specimens Reveals a New Species of Ladybird (Coleoptera: Coccinellidae): Morphological, Museomic, and Phylogenetic Analyses

Natural history collections house an important source of genetic data from yet unexplored biological diversity. Molecular data from museum specimens remain underexploited, which is mainly due to the degradation of DNA from specimens over time. However, Next-Generation Sequencing (NGS) technology can now be used to sequence "old" specimens. Indeed, many of these specimens are unique samples of nomenclatural types and can be crucial for resolving systematic or biogeographic scientific questions. Two ladybird beetle specimens from Patagonia corresponding to a new species of the genus Eriopis Mulsant were found in the collections of the Muséum national d'Histoire naturelle (MNHN), Paris. Here, we describe Eriopis patagonia Salazar, sp. nov. Total DNA of one of the two specimens was sequenced by NGS using a paired-end Illumina approach. We reconstruct and characterize the mitochondrial genome of this species (16,194 bp). Then, the protein-coding genes (PCGs) and ribosomal RNAs (rRNAs) were used to infer by maximum likelihood and Bayesian Inference the phylogenetic position of E. patagonia among 27 representatives of Coccinellidae. Phylogenetic analysis confirmed the position of Eriopis as sister group to Cycloneda Crotch. Hence, we highlight the high potential of sequencing technology for extracting molecular information from old specimens, which are used here for the systematic study of a genus, while demonstrating the importance of preserving biological collections.

Type and distribution of sensilla in the antennae of Euplatypus parallelus (F.) (Coleoptera: Curculionidea, Platypodinae)

Euplatypus parallelus (F.) (Coleoptera: Curculionidea) is the most destructive cosmopolitan insect pest of the Platypodinae. Pheromone-based luring agents are used currently in controlling bark beetle. Antennae are the primary insect organs sensing volatiles of host trees and pheromones of pioneer males. We studied the external morphology of antennae and the type, distribution, and the number of the beetle sensilla. Our results show E. parallelus have a geniculate antenna composed of 6 segments, namely the scape, 4-segmented funicle and club. Ninety-seven percent of the antennal sensors were distributed in the club, and 3% were distributed in the scape and funicle. 6 types of sensilla on the antennae were found, including sensilla trichodea (subtypes: STI, STII and STIII), sensilla basiconica (subtypes: SBI, SBII, SBIII and SBIV), sensilla chaetica (subtypes: SChI, SChII and SChIII), as well as sensilla coeloconica, sensilla campaniform and sensilla furcatea. There was no significant difference in the type, distribution and number of sensilla in males and females. No significant difference in the shape and distribution of antennae was found between sexes, but the length of antennae and the number of SChI, SChII, STI, SBI, SBIII and SBIV were significantly larger in females than males. We revealed the external cuticular structure of the antennae in E. parallelus, which can be used to guide future electrophysiological investigations to understand the ability of this beetle to detect semiochemicals.

Novel tRNA gene rearrangements in the mitochondrial genome of Camarochiloides weiweii (Hemiptera: Pachynomidae)

Pachynomidae is a poorly known heteropteran group currently considered as a distinct family closely related to Reduviidae. Taxonomic position of Pachynomidae has been debated for a long time. Knowledge about mitochondrial genome (mitogenome) of this family also remain limited. In the present study, we describe the complete mitogenome of Camarochiloides weiweii, the first mitogenome sequenced for this family. This mitogenome consists of 13 protein-coding genes, two rRNA genes and 21 tRNA genes. Two tRNA gene rearrangements in this genome, including the translocation of trnA and trnR as well as the transposition of trnH, were observed and could be explained by the tandem duplication-random loss (TDRL) model. Bayesian and maximum likelihood analyses based on sequences of three different mitogenomic datasets recovered the monophyly of Reduvioidea, and the sister relationship between Reduvioidea and the remaining Cimicomorpha lineages. This study enlightens the future application of mitogenomic data in the phylogenetic study of Cimicomorpha based on the large-scale sampling.

Mitochondrial genomes of twelve species of hyperdiverse Trigonopterus weevils

Mitochondrial genomes of twelve species of Trigonopterus weevils are presented, ten of them complete. We describe their gene order and molecular features and test their potential for reconstructing the phylogeny of this hyperdiverse genus comprising > 1,000 species. The complete mitochondrial genomes examined herein ranged from 16,501 bp to 21,007 bp in length, with an average AT content of 64.2% to 69.7%. Composition frequencies and skews were generally lower across species for atp6, cox1-3, and cob genes, while atp8 and genes coded on the minus strand showed much higher divergence at both nucleotide and amino acid levels. Most variation within genes was found at the codon level with high variation at third codon sites across species, and with lesser degree at the coding strand level. Two large non-coding regions were found, CR1 (between rrnS and trnI genes) and CR2 (between trnI and trnQ), but both with large variability in length; this peculiar structure of the non-coding region may be a derived character of Curculionoidea. The nad1 and cob genes exhibited an unusually high interspecific length variation of up to 24 bp near the 3' end. This pattern was probably caused by a single evolutionary event since both genes are only separated by trnS2 and length variation is extremely rare in mitochondrial protein coding genes. We inferred phylogenetic trees using protein coding gene sequences implementing both maximum likelihood and Bayesian approaches, each for both nucleotide and amino acid sequences. While some clades could be retrieved from all reconstructions with high confidence, there were also a number of differences and relatively low support for some basal nodes. The best partition scheme of the 13 protein coding sequences obtained by IQTREE suggested that phylogenetic signal is more accurate by splitting sequence variation at the codon site level as well as coding strand, rather than at the gene level. This result corroborated the different patterns found in Trigonopterus regarding to A+T frequencies and AT and GC skews that also greatly diverge at the codon site and coding strand levels.

Mitochondrial genomes of stick insects (Phasmatodea) and phylogenetic considerations

Phasmatodea represents an order of hemimetabolous insects. This group includes species with extreme forms of masquerade crypsis, whereby they imitate twigs, bark, lichen, moss, and leaves. In this study, we sequenced and annotated three mitochondrial genomes (mitogenomes) from Phasmatodea. The lengths of the novel mitogenomes range from 14,162 bp to 15,879 bp. The gene content and organization correspond to those inferred for the ancestral insect. We conducted phylogenetic analyses together with the existing mitogenomes of polyneopterans and mayflies. In most cases, the Phasmatodea was non-monophyletic, with Embioptera and Zoraptera nested inside. The mitogenome sequences from Embioptera and Zoraptera suffered from high substitution rates and displayed very long branches in phylogenetic trees. The monophyletic Phasmatodea was recovered only when the analysis employed the site-heterogeneous CAT-GTR model in PhyloBayes and used the nucleotide dataset PCG_nt. The Euphasmatodea was well established by various data types and inference methods. In addition, the clade Heteropterygidae and the subfamilies Lonchodinae and Necrosciinae were strongly supported. The Australasian clade Lanceocercata was recovered across analyses. However, the Clitumninae was non-monophyletic.

Macroevolution and shifts in the feeding biology of the New World scarab beetle tribe Phanaeini (Coleoptera: Scarabaeidae: Scarabaeinae)

The New World scarab beetle tribe Phanaeini contains coprophagous, necrophagous, mycetophagous and suspected myrmecophilous species. We analyse the largest tribal molecular dataset assembled, incorporating, for the first time, the enigmatic monobasic genus Megatharsis, the thalassinus group of the subgenus Coprophanaeus (Metallophanaeus), and the subgenus Dendropaemon (Eurypodea) (formerly Tetramereia), unveiling their macroevolutionary and biogeographical history in light of Cenozoic abiotic changes and inferring shifts in feeding biology through time. We recover the contentious genus Gromphas outside an otherwise monophyletic Phanaeini. We infer Megatharsis in a clade containing the apparent myrmecophilous genus Dendropaemon, within the Coprophanaeus clade, and demonstrate that the subgenus Coprophanaeus (Metallophanaeus) is polyphyletic, whilst species groups within the subgenus Coprophanaeus (Coprophanaeus) are monophyletic. Our divergence time analyses and ancestral range estimation indicate an eastern South American origin for Phanaeini in the early Eocene, with subsequent colonization of Central America and the Nearctic during the Oligocene, long before a Panamanian land bridge. A shift to necrophagy in Coprophanaeus is possibly linked to increasing Neotropical small vertebrate diversity since the Eocene and, astonishingly, myrmecophily evolved from necrophagy 35 Mya. These drastic shifts in lifestyle are not concordant with variations in diversification rates and appear unlinked to Quaternary extinction of large mammals.

Evolution of tRNA gene rearrangement in the mitochondrial genome of ichneumonoid wasps (Hymenoptera: Ichneumonoidea)

Gene rearrangements in the mitochondrial genome (mt genome) are common in certain insect groups and can be an informative character for phylogenetic reconstruction. However, knowledge of the mechanism and biases of gene rearrangement in insect mt genomes is still limited. With an accelerated rate of gene rearrangements, Hymenoptera is an important group for mt genome rearrangements diversity and for understanding the gene rearrangement evolution in mt genomes. Here, we sequenced the complete mt genome of Aphidius gifuensis and analyzed the evolution of tRNA gene rearrangements in the mt genomes of ichneumonoid wasps. Two control regions were detected in A. gifuensis and most of the tRNA rearrangement events occurred around these control regions. tRNA gene rearrangements occurred in almost all of the sequenced mt genomes of Ichneumonoidea and the gene block CR-trnI-trnQ-trnM-ND2-trnW-trnC-trnY was the main hot spot of gene rearrangement. Mapped over the backbone phylogeny of Ichneumonoidea, we found that the inversion and translocation of both trnI and trnM is likely a synapomorphic rearrangement in Braconidae. Our study also demonstrated that the gene block CR-trnI-trnQ-trnM-ND2-trnW-trnC-trnY was important for inferring the gene rearrangement dynamics in Ichneumonoidea.

Sequence analysis of mitochondrial genome of the false and phantom crane-fly Ptychoptera qinggouensis Kang, Yao and Yang, 2013 (Diptera, Ptychopteridae)

The genus Ptychoptera Meigen, 1803 is the largest genus of the family Ptychopteridae with 78 known species. In this study, we report a nearly complete mitochondrial (mt) genome of this genus, which is a circular molecule of more than 15,028 bp. The mt genome contains 13 protein-coding genes, 22 tRNA genes, 2 rRNA genes and a non-coding region. The overall base composition is A (38.1%), T (36.7%), C (14.9%), and G (10.4%), with an AT content of 74.8%. The AT content of N-strand genes (75.7%) is higher than that of the J-strand genes (71.7%). Phylogenetic analysis reveals that the monophyly of Ptychopteridae, Bibiomorpha, Tipulomorpha and Brachycera are strongly supported, and the sister group relationship between Tanyderidae and Ptychopteridae is not supported.

The mitochondrial genome of Paralimna (Paralimna) concors (Diptera: Ephydridae)

The mitogenome of Paralimna (Paralimna) concors was sequenced. The mitogenome was 16,155 bp totally, consisting of 13 protein-coding genes (PCGs), two rRNAs, and 22 transfer RNAs. The nucleotide composition biases toward A and T is 78.6％ of the entirety. All PCGs start with ATN codons except COI and ND1, and end with TAA or incomplete stop codon. Phylogenetic analyses based on 11 dipteran species supported the monophyly of Ephydroidea and the relationship of Opomyzoidea + (Ephydroidea + (Lauxanioidea + (Sphaeroceroidea + (Sciomyzoidea + Tephritoidea)))).

Complete mitochondrial genome of Coenomyia ferruginea (Scopoli). (Diptera, Coenomyiidae)

The dipteran family Coenomyiidae was firstly separated from the Xylophagidae by Akrira Nagatomi in 1975. We sequenced and annotated the mitochondrial genome of Coenomyia ferruginea, the first representative of genera Coenomyia with complete mitochondrial data. This mitogenome is 17,283 bp totally, which consists of 22 transfer RNAs, 13 protein-coding genes, 2 ribosomal RNAs, and 1 non-coding control region. All genes have the conservational arrangement like other published species of brachyceran flies. The nucleotide composition biases toward A and T, the overall A + T% was up to 75.4% of the entire mitogenome. Both Bayesian inference and ML analysis strongly supported the sister relationship between Coenomyiidae and Xylophagidae. Our results also suggested that Xylophagomorpha is the sister group to Stratiomyomorpha.

Beyond DNA barcoding: The unrealized potential of genome skim data in sample identification

Genetic tools are increasingly used to identify and discriminate between species. One key transition in this process was the recognition of the potential of the ca 658bp fragment of the organelle cytochrome c oxidase I (COI) as a barcode region, which revolutionized animal bioidentification and lead, among others, to the instigation of the Barcode of Life Database (BOLD), containing currently barcodes from >7.9 million specimens. Following this discovery, suggestions for other organellar regions and markers, and the primers with which to amplify them, have been continuously proposed. Most recently, the field has taken the leap from PCR-based generation of DNA references into shotgun sequencing-based "genome skimming" alternatives, with the ultimate goal of assembling organellar reference genomes. Unfortunately, in genome skimming approaches, much of the nuclear genome (as much as 99% of the sequence data) is discarded, which is not only wasteful, but can also limit the power of discrimination at, or below, the species level. Here, we advocate that the full shotgun sequence data can be used to assign an identity (that we term for convenience its "DNA-mark") for both voucher and query samples, without requiring any computationally intensive pretreatment (e.g. assembly) of reads. We argue that if reference databases are populated with such "DNA-marks," it will enable future DNA-based taxonomic identification to complement, or even replace PCR of barcodes with genome skimming, and we discuss how such methodology ultimately could enable identification to population, or even individual, level.

Complete mitochondrial genome of Niphades castanea (Coleoptera: Curculionidae)

Niphades castanea Chao is an important insect pest on many plants which belong to genus Castanea. The complete mitochondrial genome of N. castanea was sequenced and analyzed. The phylogenetic relationships between N. castanea and other 41 species in the family Curculionidae were reconstructed using maximum likelihood (ML) methods based on the concatenated nucleotide sequences, the phylogenetic analysis showed that N. castanea is closely related to Hylobitelus xiaoi, which is in accordance with the traditional morphological classification.

Salix transect of Europe: records of willow-associated weevils (Coleoptera: Curculionoidea) from Greece to Arctic Norway, with insights from DNA barcoding

Background

Curculionid beetles associated with willow (Salix spp.) were surveyed at 42 sites across Europe, from Greece (lat. 38.8 °N) to arctic Norway (lat. 69.7 °N). DNA sequence data provide additional verification of identifications and geographic clustering.

New information

In all, 73 curculionid species were collected from willows, of which seven were particularly abundant. The most widespread species were: Acalyptuscarpini Fabricius, 1793 at 15 sites; Tachyergesstigma Germar, 1821 at 13 sites; Phyllobiusoblongus (Linnaeus, 1758) at 11 sites; Phyllobiusmaculicornis Germar, 1824 at 10 sites; and Archariussalicivorus (Paykull, 1792), Melanapionminimum (Herbst, 1797), and Phyllobiuscf.pyri (Linnaeus, 1758) all at nine sites. The mean number of curculionid species collected on willow at each site was 5.5 (range 0-14). Compared to chrysomelids, curculionids were richer in species but the species had relatively low average abundance. Widespread curculionid species appear to have scattered and patchy observed distributions with limited geographical structuring in our data. However, deeper sampling (e.g. over multiple seasons and years), would give a better indication of distribution, and may increase apparent geographical structuring. There is some site-to-site variation in colour in a few taxa, but little notable size variation. DNA barcoding, performed on some of the more common species, provides clear species clusters and definitive separation of the taxonomically more challenging species, as well as some interesting geographic insights. Our northernmost sample of Phyllobiusoblongus is unique in clustering with Canadian samples of this species. On the other hand, our samples of Acalyptuscarpini cluster with European samples and are distinct from a separate Canadian cluster of this species. We provide the first available DNA sequences for Phyllobiusthalassinus Gyllenhal, 1834 (Hungary).