Accelerated evolution of mitochondrial but not nuclear genomes of Hymenoptera: new evidence from crabronid wasps

Reassessing Hybridisation in Australian Tetragonula Stingless Bees Using Multiple Genetic Markers

We re-examined reports of hybridisation in three cryptic stingless bee species in the genus Tetragonula in South East Queensland, Australia (T. carbonaria, T. davenporti and T. hockingsi). Previous studies on this group using microsatellite markers proposed that hybridisation occasionally takes place. In contrast, we find that using 1745 SNPs we could reliably separate the three species, with no evidence of contemporary (or recent) hybridisation. We found identical amplicon sequences of the nuclear gene EF1alpha across most individuals of the three species, but low and moderate species-specific polymorphisms in the nuclear gene Opsin and the mitochondrial 16S rRNA gene, respectively, with no cases of mito-nuclear discordance at these genes. We confirm that nuclear divergence across these species is low, based on 10-26 kb of non-coding sequence flanking EF1alpha and Opsin (0.7%-1% pairwise difference between species). However, we find mitogenomes to be far more diverged than nuclear genomes (21.6%-23.6% pairwise difference between species). Based on these comprehensive analyses of multiple marker types, we conclude there is no ongoing gene flow among the Tetragonula species of South East Queensland, despite their morphological similarity to one another and the low nuclear divergence among them. The higher resolution provided by multiple SNP markers may lead to lower estimates of contemporary hybridisation more generally.

High Microsatellite but No Mitochondrial DNA Variation in an Invasive Japanese Mainland Population of the Parasitoid Wasp Melittobia sosui

Invasive populations are predicted to have reduced genetic diversity due to bottleneck events. The parasitoid wasp Melittobia sosui was previously identified only in the subtropical area of the southern Japanese islands and Taiwan but was recently found in the temperate area of the Japanese mainland. The distribution of this species may have recently expanded northward due to factors such as climatic events and global warming. The population genetics of both the native and invasive regions were investigated using mitochondrial and nuclear microsatellite DNA. As expected, mitochondrial variation was observed in the native region but not in the invasive region, which had only one haplotype. However, the two regions exhibited similar levels of microsatellite variation, and an average of 43% and 38% of alleles were uniquely found in the native and invasive populations, respectively. The difference in genetic variation between mitochondrial and microsatellite DNA in the invasive populations may be explained by the faster mutation rate of microsatellites, as well as the population structure of Melittobia, in which the subdivision into small inbreeding lineages may facilitate the accumulation of mutations. The high proportion of private alleles suggests that the mainland population diverged from the native populations at least 100 years ago, ruling out the possibility that the mainland population was established recently. The present study suggests that M. sosui might have already existed on the mainland but at a low frequency or that the mainland population was derived from a ghost population that diverged from the native populations more than 100 years ago.

Barcode 100K Specimens: In a Single Nanopore Run

It is a global priority to better manage the biosphere, but action must be informed by comprehensive data on the abundance and distribution of species. The acquisition of such information is currently constrained by high costs. DNA barcoding can speed the registration of unknown animal species, the most diverse kingdom of eukaryotes, as the BIN system automates their recognition. However, inexpensive sequencing protocols are critical as the census of all animal species is likely to require the analysis of a billion or more specimens. Barcoding involves DNA extraction followed by PCR and sequencing with the last step dominating costs until 2017. By enabling the sequencing of highly multiplexed samples, the Sequel platforms from Pacific BioSciences slashed costs by 90%, but these instruments are only deployed in core facilities because of their expense. Sequencers from Oxford Nanopore Technologies provide an escape from high capital and service costs, but their low sequence fidelity has, until recently, constrained adoption. However, the improved performance of its latest flow cells (R10.4.1) erases this barrier. This study demonstrates that a MinION flow cell can characterise an amplicon pool derived from 100,000 specimens while a Flongle flow cell can process one derived from several thousand. At $0.01 per specimen, DNA sequencing is now the least expensive step in the barcode workflow.

Rapid evolution of mitochondrion-related genes in haplodiploid arthropods

BACKGROUND

Mitochondrial genes and nuclear genes cooperate closely to maintain the functions of mitochondria, especially in the oxidative phosphorylation (OXPHOS) pathway. However, mitochondrial genes among arthropod lineages have dramatic evolutionary rate differences. Haplodiploid arthropods often show fast-evolving mitochondrial genes. One hypothesis predicts that the small effective population size of haplodiploid species could enhance the effect of genetic drift leading to higher substitution rates in mitochondrial and nuclear genes. Alternatively, positive selection or compensatory changes in nuclear OXPHOS genes could lead to the fast-evolving mitochondrial genes. However, due to the limited number of arthropod genomes, the rates of evolution for nuclear genes in haplodiploid species, besides hymenopterans, are largely unknown. To test these hypotheses, we used data from 76 arthropod genomes, including 5 independently evolved haplodiploid lineages, to estimate the evolutionary rates and patterns of gene family turnover of mitochondrial and nuclear genes.

RESULTS

We show that five haplodiploid lineages tested here have fast-evolving mitochondrial genes and fast-evolving nuclear genes related to mitochondrial functions, while nuclear genes not related to mitochondrion showed no significant evolutionary rate differences. Among hymenopterans, bees and ants show faster rates of molecular evolution in mitochondrial genes and mitochondrion-related nuclear genes than sawflies and wasps. With genome data, we also find gene family expansions and contractions in mitochondrion-related genes of bees and ants.

CONCLUSIONS

Our results reject the small population size hypothesis in haplodiploid species. A combination of positive selection and compensatory changes could lead to the observed patterns in haplodiploid species. The elevated evolutionary rates in OXPHOS complex 2 genes of bees and ants suggest a unique evolutionary history of social hymenopterans.

Dating in the Dark: Elevated Substitution Rates in Cave Cockroaches (Blattodea: Nocticolidae) Have Negative Impacts on Molecular Date Estimates

Rates of nucleotide substitution vary substantially across the Tree of Life, with potentially confounding effects on phylogenetic and evolutionary analyses. A large acceleration in mitochondrial substitution rate occurs in the cockroach family Nocticolidae, which predominantly inhabit subterranean environments. To evaluate the impacts of this among-lineage rate heterogeneity on estimates of phylogenetic relationships and evolutionary timescales, we analyzed nuclear ultraconserved elements (UCEs) and mitochondrial genomes from nocticolids and other cockroaches. Substitution rates were substantially elevated in nocticolid lineages compared with other cockroaches, especially in mitochondrial protein-coding genes. This disparity in evolutionary rates is likely to have led to different evolutionary relationships being supported by phylogenetic analyses of mitochondrial genomes and UCE loci. Furthermore, Bayesian dating analyses using relaxed-clock models inferred much deeper divergence times compared with a flexible local clock. Our phylogenetic analysis of UCEs, which is the first genome-scale study to include all 13 major cockroach families, unites Corydiidae and Nocticolidae and places Anaplectidae as the sister lineage to the rest of Blattoidea. We uncover an extraordinary level of genetic divergence in Nocticolidae, including two highly distinct clades that separated ~115 million years ago despite both containing representatives of the genus Nocticola. The results of our study highlight the potential impacts of high among-lineage rate variation on estimates of phylogenetic relationships and evolutionary timescales.

Adaptive evidence of mitochondrial genes in Pteromalidae and Eulophidae (Hymenoptera: Chalcidoidea)

Pteromalidae and Eulophidae are predominant and abundant taxa within Chalcidoidea (Hymenoptera: Apocrita). These taxa are found in diverse ecosystems, ranging from basin deserts (200 m) to alpine grasslands (4500 m). Mitochondria, cellular powerhouses responsible for energy production via oxidative phosphorylation, are sensitive to various environmental factors such as extreme cold, hypoxia, and intense ultraviolet radiation characteristic of alpine regions. Whether the molecular evolution of mitochondrial genes in these parasitoids corresponds to changes in the energy requirements and alpine environmental adaptations remains unknown. In this study, we performed a comparative analysis of mitochondrial protein-coding genes from 11 alpine species of Pteromalidae and Eulophidae, along with 18 lowland relatives, including 16 newly sequenced species. We further examined the codon usage preferences (RSCU, ENC-GC3s, neutrality, and PR2 bias plot) in these mitochondrial protein-coding sequences and conducted positive selection analysis based on their Bayesian phylogenetic relationships, and identified positive selection sites in the ATP6, ATP8, COX1, COX3, and CYTB genes, emphasizing the crucial role of mitochondrial gene adaptive evolution in the adaptation of Pteromalidae and Eulophidae to alpine environments. The phylogenetically independent contrast (PIC) analysis results verified the ω ratio of 13 PCGs from Pteromalidae and Eulophidae increased with elevation, and results from generalized linear model confirm that ATP6, ATP8, COX3, and ND1 are closely correlated with temperature-related environmental factors. This research not only enriched the molecular data of endemic alpine species but also underscores the significance of mitochondrial genes in facilitating the adaptation of these minor parasitoids to plateau habitats.

New species based on the biological species concept within the complex of Lariophagus distinguendus (Hymenoptera, Chalcidoidea, Pteromalidae), a parasitoid of household pests

The pteromalid parasitoid Lariophagus distinguendus (Foerster) belongs to the Hymenoptera, a megadiverse insect order with high cryptic diversity. It attacks stored product pest beetles in human storage facilities. Recently, it has been shown to consist of two separate species. To further study its cryptic diversity, strains were collected to compare their relatedness using barcoding and nuclear genes. Nuclear genes identified two clusters which agree with the known two species, whereas the barcode fragment determined an additional third Clade. Total reproductive isolation (RI) according to the biological species concept (BSC) was investigated in crossing experiments within and between clusters using representative strains. Sexual isolation exists between all studied pairs, increasing from slight to strong with genetic distance. Postzygotic barriers mostly affected hybrid males, pointing to Haldane's rule. Hybrid females were only affected by unidirectional Spiroplasma-induced cytoplasmic incompatibility and behavioural sterility, each in one specific strain combination. RI was virtually absent between strains separated by up to 2.8% COI difference, but strong or complete in three pairs from one Clade each, separated by at least 7.2%. Apparently, each of these clusters represents one separate species according to the BSC, highlighting cryptic diversity in direct vicinity to humans. In addition, these results challenge the recent 'turbo-taxonomy' practice of using 2% COI differences to delimitate species, especially within parasitic Hymenoptera. The gradual increase in number and strength of reproductive barriers between strains with increasing genetic distance also sheds light on the emergence of barriers during the speciation process in L. distinguendus.

Inter- and Intrasexual Variation in Cuticular Hydrocarbons in Trichrysis cyanea (Linnaeus, 1758) (Hymenoptera: Chrysididae)

Simple Summary

The biology of many species of cuckoo wasps (Chrysididae) is largely unknown, except for, if at all, the host species of these brood parasites. Cuticular hydrocarbons (CHCs), which form a waxy layer on the body surface of insects, have been shown to be mimicked by host-specific cuckoo wasp species. We studied the CHC profiles of a rather generalist brood parasitic cuckoo wasp species, Trichrysis cyanea, in detail. We found sex-related differences in CHCs and three different chemotypes among females. Genetic data reject the possibility that these different chemotypes represent different (cryptic) species. The CHC polymorphism could be an adaptation for females to sneak into nesting holes of hosts with the different female chemotypes representing adaptations to a broader host range. However, since information on the CHC profiles of the hosts is missing, it remains unclear whether these different CHC profiles are used for chemical mimicry or are simply age-related, reflecting the females’ reproductive state.

Abstract

Cuckoo wasps (Chrysididae, Hymenoptera) are known for their parasitoid or cleptoparasitic life histories. Indeed, the biology of only a few species has been studied in detail and often only little more is known than the host species. By mimicking their hosts’ cuticular hydrocarbon (CHC) profiles, species that parasitize single (or a few closely related) host species manage to deceive their hosts. However, the variability of the CHC profile in generalist cuckoo-wasp species is still unknown. Here, we used gas chromatography—mass spectrometry (GC-MS) and DNA barcoding to study intraspecific variation in cuticular hydrocarbons of one less host-specific species of cuckoo wasps, Trichrysis cyanea. Cuticular hydrocarbon (CHC) patterns were found to differ between males and females. Additionally, we found chemical polymorphism among females, which formed three distinct chemical subgroups characterized by different alkene patterns. A lack of divergence in the DNA barcoding region suggests that these different chemotypes do not represent cryptic species. Whether this intrasexual CHC-profile variation is an adaptation (mimicry) to different host species, or simply signaling the reproductive status, remains unclear.

Comparative Mitochondrial Genomics of 104 Darwin Wasps (Hymenoptera: Ichneumonidae) and Its Implication for Phylogeny

Simple Summary

Nearly a hundred mitochondrial genomes of ichneumonid wasps are newly reported. Comparative mitogenomics of 104 mitochondrial genomes representing 33 subfamilies of Ichneumonidae, as well as its implications for phylogeny, were studied. We found that the mitochondrial genomes of ichneumonid wasps were highly conserved in their base composition and had low evolutionary rates, but were diverse in gene order. There are 38 types of gene rearrangement events in 104 ichneumonid mitochondrial genomes, of which 30 novel rearrangement types (R3-6, R8-R10, R12-R15, R17-R18, R20-R35 and R38) and a hot spot rearrangement around R1, with a shuffled tRNA cluster trnW-trnY-trnC and trnI-trnQ-trnM, were detected. The relationships among these subfamilies are firstly discussed based on mitochondrial genomes at a large scale. We suggest five subfamily groupings of Ichneumonidae: Brachycyrtiformes, Ichneumoniformes, Ophioniformes, Pimpliformes and Xoridiformes. Two formerly unplaced subfamilies, Eucerotinae and Microleptinae, were placed in Brachycyrtiformes and Ichneumoniformes, respectively.

Abstract

Ichneumonidae is one of the largest families of insects with a mega-diversity of specialized morphological and biological characteristics. We newly sequenced 92 mitochondrial genomes of ichneumonid wasps and found that they have a conserved base composition and a lower evolutionary rate than that of other families of parasitic Hymenoptera. There are 38 types of gene order in the ichneumonid mitochondrial genome, with 30 novel types identified in 104 ichneumonids. We also found that the rearrangement events occur more frequently in Ophioniformes than in Ichneumoniformes and Pimpliformes. Furthermore, the higher Ophioniformes and their relative lineages shared the transposition of trnL2 to trnI-trnQ-trnM tRNA cluster. We confirmed five higher-level groupings of Ichneumonidae: Brachycyrtiformes, Ichneumoniformes, Ophioniformes, Pimpliformes and Xoridiformes. Two formerly unplaced subfamilies, Eucerotinae and Microleptinae, were placed in Brachycyrtiformes and Ichneumoniformes, respectively. The results will improve our understanding of the diversity and evolution of Ichneumonidae.

Assessment of genetic variation in Apis mellifera jemenitica (Hymenoptera: Apidae) using Cytochrome Oxidase I gene sequences

The Arabian Honeybee Apis mellifera jemenitica is endemic to the Arabian Peninsula. It is highly adapted to temperature extremes and drought dominating the region. In this study, the mitochondrial Cytochrome Oxidase I (COI) was analyzed in 133 specimens of A. m. jemenitica from eight localities along the Red Sea cost of Saudi Arabia. Results revealed 33 synonymous, and 6 non-synonymous mutations within the COI sequences, resulting in change of 4 amino acids. Phylogenetic analysis based on either type of mutations revealed two main haplogroups accounting for 94% of the samples. In total Eighteen new haplotypes were identified and uploaded in the genebank, Fourteen of them are restricted to one/both haplogroups. All haplotypes identified in this study clustered with reference COI sequences of the sub-lineag Z (African Lineage). However one Haplotype (MW428270) represents high COI variability compared to other haplotypes and may resemble different evolutionary sub-lineage. Tajima's Neutrality Test (Ps = 0.025; D = -1.5) indicated population size expansion that took place after selective sweep and/or purifying selection.

Bivalve molluscs as model systems for studying mitochondrial biology

The class Bivalvia is a highly successful and ancient taxon including ∼25,000 living species. During their long evolutionary history bivalves adapted to a wide range of physicochemical conditions, habitats, biological interactions, and feeding habits. Bivalves can have strikingly different size, and despite their apparently simple body plan, they evolved very different shell shapes, and complex anatomic structures. One of the most striking features of this class of animals is their peculiar mitochondrial biology: some bivalves have facultatively anaerobic mitochondria that allow them to survive prolonged periods of anoxia/hypoxia. Moreover, more than 100 species have now been reported showing the only known evolutionarily stable exception to the strictly maternal inheritance of mitochondria in animals, named doubly uniparental inheritance. Mitochondrial activity is fundamental to eukaryotic life, and thanks to their diversity and uncommon features, bivalves represent a great model system to expand our knowledge about mitochondrial biology, so far limited to a few species. We highlight recent works studying mitochondrial biology in bivalves at either genomic or physiological level. A link between these two approaches is still missing, and we believe that an integrated approach and collaborative relationships are the only possible ways to be successful in such endeavour.

Comparative mitogenomics and phylogenetics of the stinging wasps (Hymenoptera: Aculeata)

The stinging wasps (Hymenoptera: Aculeata) include diverse groups such as vespid wasps, ants and bees. Phylogenetic relationships among major lineages of stinging wasps have been inferred from molecular and morphological data. However, the genomic features of the mitochondrial genomes and their phylogenetic utility remain to be explored. In this study, we determined 23 mitochondrial genomes from the Aculeata. Four Mutillidae species showed relatively low A + T content compared to other species of the Aculeata (69.7%-77.4%). Eleven out of 44 species, mainly from the Chrysididae and the Pompilidae, showed reversals of GC skews. Gene rearrangements occurred across the species. Patterns of tRNA rearrangement were conserved in some groups, including the Chrysididae, Bethylidae, Pompilidae, Scolioidea and Vespoidea. Rearrangement of protein-coding genes were found in 12 out of 44 species of the Aculeata, including all four species from the Chrysididae, both species from the Bethylidae, one species from the Dryinidae, all three Scolioidea species and two Apoidea species. Phylogenetic inference showed a long branch in species with unusual genomic features, such as in the Mutillidae and Bethylidae. By excluding these species, we found paraphyly of the Chrysidoidea and a sister group relationship between the Formicoidea and Vespoidea. These results improve our understanding of the evolution of mitochondrial genomes in the Aculeata and, in general, the evolution across this subclade.

Inferring Evolutionary Timescales without Independent Timing Information: An Assessment of "Universal" Insect Rates to Calibrate a Collembola (Hexapoda) Molecular Clock

Previous estimates of nucleotide substitution rates are routinely applied as secondary or “universal” molecular clock calibrations for estimating evolutionary timescales in groups that lack independent timing information. A major limitation of this approach is that rates can vary considerably among taxonomic groups, but the assumption of rate constancy is rarely evaluated prior to using secondary rate calibrations. Here I evaluate whether an insect mitochondrial DNA clock is appropriate for estimating timescales in Collembola—a group of insect-like arthropods characterized by high levels of cryptic diversity. Relative rates of substitution in cytochrome oxidase subunit 1 (COI) were inferred via Bayesian analysis across a topologically constrained Hexapod phylogeny using a relaxed molecular clock model. Rates for Collembola did not differ significantly from the average rate or from the rates estimated for most other groups (25 of 30), suggesting that (1) their apparent cryptic diversity cannot be explained by accelerated rates of molecular evolution and (2) clocks calibrated using “universal” insect rates may be appropriate for estimating evolutionary timescales in this group. However, of the 31 groups investigated, 10 had rates that deviated significantly from the average (6 higher, 4 lower), underscoring the need for caution and careful consideration when applying secondary insect rate calibrations. Lastly, this study exemplifies a relatively simple approach for evaluating rate constancy within a taxonomic group to determine whether the use of secondary rates are appropriate for molecular clock calibrations.

One Hundred Mitochondrial Genomes of Cicadas

Mitochondrial genomes can provide valuable information on the biology and evolutionary histories of their host organisms. Here, we present and characterize the complete coding regions of 107 mitochondrial genomes (mitogenomes) of cicadas (Insecta: Hemiptera: Auchenorrhyncha: Cicadoidea), representing 31 genera, 61 species, and 83 populations. We show that all cicada mitogenomes retain the organization and gene contents thought to be ancestral in insects, with some variability among cicada clades in the length of a region between the genes nad2 and cox1, which encodes 3 tRNAs. Phylogenetic analyses using these mitogenomes recapitulate a recent 5-gene classification of cicadas into families and subfamilies, but also identify a species that falls outside of the established taxonomic framework. While protein-coding genes are under strong purifying selection, tests of relative evolutionary rates reveal significant variation in evolutionary rates across taxa, highlighting the dynamic nature of mitochondrial genome evolution in cicadas. These data will serve as a useful reference for future research into the systematics, ecology, and evolution of the superfamily Cicadoidea.

Molecular ecology as a tool for understanding pollination and other plant-insect interactions

Advances in molecular ecology offer unprecedented opportunities to understand the ecology and evolution of insects, the complex ways in which they interact and their role in ecosystem functioning. Rapidly developing DNA sequencing technologies are resolving previously intractable questions in taxonomic and functional biodiversity and provide significant potential to determine formerly difficult to observe plant-insect interactions. We provide an overview of the state-of-the-art and critically appraise the range of molecular approaches currently available for the study of insect pollination, host-parasitoid interactions and/or wider food-web studies. Species-interaction data are increasingly being incorporated into ecological network analyses. DNA metabarcoding offers opportunities to scale-up efforts to create large, highly resolved, phylogenetically structured networks within an exciting framework to study pressing questions in ecology and evolution.

Novel gene rearrangement in the mitochondrial genome of Pachyneuron aphidis (Hymenoptera: Pteromalidae)

Species in Hymenoptera usually show exceptionally high rates of mitochondrial molecular evolution and dramatic gene rearrangements, which has been attributed to their parasitic lifestyle. However, mitochondrial genome (mitogenome) of chalcidoid wasps is still poorly understood and the evolution of gene rearrangements is still unclear. In this study, the nearly complete mitogenome of Pachyneuron aphidis, a chalcidoid wasp mainly hyperparasitizes the Aphidius gifuensis, was sequenced using a next-generation sequencing strategy. This genome is 15,137 bp in length, including 13 PCGs, 22 tRNAs, two rRNAs and a partial control region. Alignment with other Chalcidoidea mitogenomes revealed a novel inversion in the srRNA-trnV gene cluster in P. aphidis, which is the first of its kind to be reported in Chalcidoidea. Breakpoint distances analysis showed the high value of chalcidoid wasps compare to the ancestral arrangement pattern, which reflected as extensive gene rearrangements. Despite the high frequency of gene rearrangements in these insects, analyses of gene rearrangement and phylogenetic trees showed that species from the same family and the genus tent to have similar gene orders, and the conserved gene blocks (ND3-trnG, srRNA-trnV and COIII-ATP6-ATP8-trnD-trnK-COII-trnL2-COI) can usually be identified, especially at the family level of chalcidoid wasps.

Cross-contamination and strong mitonuclear discordance in Empria sawflies (Hymenoptera, Tenthredinidae) in the light of phylogenomic data

In several sawfly taxa strong mitonuclear discordance has been observed, with nuclear genes supporting species assignments based on morphology, whereas the barcode region of the mitochondrial COI gene suggests different relationships. As previous studies were based on only a few nuclear genes, the causes and the degree of mitonuclear discordance remain ambiguous. Here, we obtained genomic-scale ddRAD data together with Sanger sequences of mitochondrial COI and two to three nuclear protein coding genes to investigate species limits and mitonuclear discordance in two closely related species groups of the sawfly genus Empria. As found previously based on nuclear ITS and mitochondrial COI sequences, species are in most cases supported as monophyletic based on new nuclear data reported here, but not based on mitochondrial COI. This mitonuclear discordance can be explained by occasional mitochondrial introgression with little or no nuclear gene flow, a pattern that might be common in haplodiploid taxa with slowly evolving mitochondrial genomes. Some species in the E. immersa group are not recovered as monophyletic according to either mitochondrial or nuclear data, but this could partly be because of unresolved taxonomy. Preliminary analyses of ddRAD data did not recover monophyly of E. japonica within the E. longicornis group (three Sanger sequenced nuclear genes strongly supported monophyly), but closer examination of the data and additional Sanger sequencing suggested that both specimens were substantially (possibly 10-20% of recovered loci) cross-contaminated. A reason could be specimen identification tag jumps during sequencing library preparation that in previous studies have been shown to affect up to 2.5% of the sequenced reads. We provide an R script to examine patterns of identical loci among the specimens and estimate that the cross-contamination rate is not unusually high for our ddRAD dataset as a whole (based on counting of identical sequences in the immersa and longicornis groups, which are well separated from each other and probably do not hybridise). The high rate of cross-contamination for both E. japonica specimens might be explained by the small number of recovered loci (~1000) compared to most other specimens (>10 000 in some cases) because of poor sequencing results. We caution against drawing unexpected biological conclusions when closely related specimens are pooled before sequencing and tagged only at one end of the molecule or at both ends using a unique combination of limited number of tags (less than the number of specimens).

Characterization of Hymenopteran Parasitoids of Aphis fabae in an African Smallholder Bean Farming System through Sequencing of COI 'Mini-Barcodes'

Parasitoids are among the most frequently reported natural enemies of insect pests, particularly aphids. The efficacy of parasitoids as biocontrol agents is influenced by biotic and abiotic factors. For example, hyperparasitoids can reduce the abundance of the primary parasitoids as well as modify their behavior. A field study was conducted at three contrasting elevations on Mount Kilimanjaro, Tanzania, to identify the parasitoids of aphids in smallholder bean farming agroecosystems. Sentinel aphids (Aphis fabae) on potted bean plants (Phaseolus vulgaris) were exposed in 15 bean fields at three elevations for 2 days. The sentinel aphids were then kept in cages in a greenhouse until emergence of the parasitoids, which were collected and preserved in 98% ethanol for identification. Of the 214 parasitoids that emerged from sentinel aphids, the greatest abundance (44.86%) were from those placed at intermediate elevations (1000–1500 m a.s.l), compared to 42.52% from the lowest elevations and only 12.62% from the highest elevation farms. Morphological identification of the parasitoids that emerged from parasitized aphids showed that 90% were Aphidius species (Hymenoptera: Braconidae: Aphidiinae). Further characterization by sequencing DNA ‘mini-barcodes’ identified parasitoids with ≥99% sequence similarity to Aphidius colemani, 94–95% sequence similarity to Pachyneuron aphidis and 90% similarity to a Charipinae sp. in the National Center for Biotechnology Information (NCBI) database. These results confidently identified A. colemani as the dominant primary aphid parasitoid of A. fabae in the study area. A Pachyneuron sp., which was most closely related to P. aphidis, and a Charipinae sp. occurred as hyperparasitoids. Thus, interventions to improve landscapes and farming practice should monitor specifically how to augment populations of A. colemani, to ensure any changes enhance the delivery of natural pest regulation. Further studies are needed for continuous monitoring of the hyperparasitism levels and the dynamics of aphids, primary parasitoids, and secondary parasitoids in different cropping seasons and their implications in aphid control.

Ecological and geographical speciation in Lucilia bufonivora: The evolution of amphibian obligate parasitism

Lucilia (Diptera: Calliphoridae) is a genus of blowflies comprised largely of saprophagous and facultative parasites of livestock. Lucilia bufonivora, however, exhibits a unique form of obligate parasitism of amphibians, typically affecting wild hosts. The evolutionary route by which amphibian myiasis arose, however, is not well understood due to the low phylogenetic resolution in existing nuclear DNA phylogenies. Furthermore, the timing of when specificity for amphibian hosts arose in L. bufonivora is also unknown. In addition, this species was recently reported for the first time in North America (Canada) and, to date, no molecular studies have analysed the evolutionary relationships between individuals from Eastern and Western hemispheres. To provide broader insights into the evolution of the amphibian parasitic life history trait and to estimate when the trait first arose, a time-scaled phylogeny was inferred from a concatenated data set comprising mtDNA, nDNA and non-coding rDNA (COX1, per and ITS2 respectively). Specimens from Canada, the UK, Poland, Switzerland, the Netherlands and Germany were analysed, as well as individuals from its sister taxa, the saprophage Lucilia silvarum and a Nearctic species also implicated in amphibian myiasis, Lucilia elongata. Obligate amphibian parasitism appears to have arisen ~4 mya, likely as a result of niche displacement of a saprophagous/facultative parasite ancestor. Consistent paraphyly of L. bufonivora with respect to L. elongata across single-gene phylogenies and high mtDNA genetic distances between Nearctic and Palearctic individuals suggest on-going cryptic speciation facilitated by geographical isolation. These findings suggest that recent reports of L. bufonivora in the Nearctic do not constitute a recent introduction, but instead suggest that it remained unrecorded due to taxonomic confusion and low abundance. This is the first study to confirm the involvement of L. bufonivora in amphibian myiasis in Canada using DNA-based identification methods.

•

Within Lucilia, a small genus of blowflies mostly comprised of carrion-breeding species, obligate parasitism for amphibians evolved around 4 mya.

•

Geographic isolation between Nearctic and Palearctic lineages is facilitating on-going cryptic speciation of Lucilia bufonivora.

•

First positive identification of L. bufonivora from two confirmed cases of amphibian myiasis in North America.

Evolutionary Rate Correlation between Mitochondrial-Encoded and Mitochondria-Associated Nuclear-Encoded Proteins in Insects

The mitochondrion is a pivotal organelle for energy production, and includes components encoded by both the mitochondrial and nuclear genomes. Functional and evolutionary interactions are expected between the nuclear- and mitochondrial-encoded components. The topic is of broad interest in biology, with implications to genetics, evolution, and medicine. Here, we compare the evolutionary rates of mitochondrial proteins and ribosomal RNAs to rates of mitochondria-associated nuclear-encoded proteins, across the major orders of holometabolous insects. There are significant evolutionary rate correlations (ERCs) between mitochondrial-encoded and mitochondria-associated nuclear-encoded proteins, which are likely driven by different rates of mitochondrial sequence evolution and correlated changes in the interacting nuclear-encoded proteins. The pattern holds after correction for phylogenetic relationships and considering protein conservation levels. Correlations are stronger for both nuclear-encoded OXPHOS proteins that are in contact with mitochondrial OXPHOS proteins and for nuclear-encoded mitochondrial ribosomal amino acids directly contacting the mitochondrial rRNAs. We find that ERC between mitochondrial- and nuclear-encoded proteins is a strong predictor of nuclear-encoded proteins known to interact with mitochondria, and ERC shows promise for identifying new candidate proteins with mitochondrial function. Twenty-three additional candidate nuclear-encoded proteins warrant further study for mitochondrial function based on this approach, including proteins in the minichromosome maintenance helicase complex.

The first divergence time estimation of the subfamily Stenogastrinae (Hymenoptera: Vespidae) based on mitochondrial phylogenomics

In this study, the mitochondrial genomes of three Stenogastrinae species, Eustenogaster scitula, Liostenogaster nitidipennis and Parishnogaster mellyi were sequenced and annotated, and a total of 16 vespid mtgenomes are comparatively analyzed. Our results indicate that codon usage bias is mainly dominated by mutational pressure, and affected only slightly by natural selection. Selective pressure analysis of protein-coding genes (PCGs) shows that the highest evolutionary rate is present in NADH complex I, and the lowest in cox1. Compared with the reported mtgenomes of other Vespidae, in Stenogastrinae, trnH is shifted to a new position. Phylogenetic analyses are performed using Bayesian method and Maximum Parsimony. Phylogenetic analysis further confirms that the Stenogastrinae is the sister group of all remaining Vespidae. Divergence time of Stenogastrinae from other Vespidae is estimated at ~ 166 Mya. Our results also support that eusociality evolved twice in the family Vespidae.

Hymenoptera of Canada

A summary of the numbers of species of the 83 families of Hymenoptera recorded in Canada is provided. In total, 8757 described species are recorded compared to approximately 6000 in 1979, which is a 46% increase. Of the families recognized in 1979, three have been newly recorded to Canada since the previous survey: Anaxyelidae (Anaxyleoidea), Liopteridae (Cynipoidea), and Mymarommatidae (Mymarommatoidea). More than 18,400 BINs of Canadian Hymenoptera are available in the Barcode of Life Data Systems (Ratnasingham and Hebert 2007) implying that nearly 9650 undescribed or unrecorded species of Hymenoptera may be present in Canada (and more than 10,300 when taking into account additional species that have not been DNA barcoded). The estimated number of unrecorded species is very similar to that of 1979 (10,637 species), but the percentage of the fauna described/recorded has increased from 36% in 1979 to approximately 45% in 2018. Summaries of the state of knowledge of the major groups of Hymenoptera are presented, including brief comments on numbers of species, biology, changes in classification since 1979, and relevant taxonomic references.

The first mitogenomes of the superfamily Pamphilioidea (Hymenoptera: Symphyta): Mitogenome architecture and phylogenetic inference

The Pamphilioidea represents a small superfamily of the phytophagous suborder Symphyta (Hymenoptera). Here, nearly complete mitochondrial genomes (mitogenomes) of three pamphilioid species: Chinolyda flagellicornis (Pamphiliidae), Megalodontes spiraeae and M. cephalotes (Megalodontesidae) were newly sequenced using next generation sequencing and comparatively analysed with the previously reported symphytan mitogenomes. A positive AT skew (0.013) and a negative GC skew (-0.194) were found in pamphilioid mitogenome, and a deviation from strand asymmetry was also observed in the PCGs encoded on both strands. Several gene rearrangement events were observed in four tRNA gene clusters (WCY, IQM, ARNS1EF and TP clusters), which have not been reported from symphytan mitogenomes to date. As the most parsimonious explanation, compared with the inferred insect ancestral mitogenome architecture, the occurrence of gene rearrangements in pamphilioid mitogenomes requires totally five evolutionary steps, including four transpositions and one inversion. The predicted secondary structures of tRNAs, rrnS and rrnL genes are mostly consistent with reported hymenopteran species. Phylogenetic analyses recovered the monophyly of superfamily Pamphilioidea and indicated the relationship Tenthredinoidea + (Pamphilioidea + (Cephoidea + (Orussoidea + Apocrita))) with strong nodal supports.

The first two mitochondrial genomes of the family Aphelinidae with novel gene orders and phylogenetic implications

Chalcidoidea is one of the most diverse group in Hymenoptera by possessing striking mitochondrial gene arrangement. By using next generation sequencing method, the first two nearly complete mitochondrial genomes in the family Aphelinidae (Insecta, Hymenopetra, Chalcidoidea) were obtained in this study. Almost all previously sequenced mitochondrial genome of Chalcidoidea species have a large inversion including six genes (atp6-atp8-trnD-trnK-cox2-trnL2-cox1) as compared with ancestral mitochondrial genome, but these two Encarsia mitochondrial genomes had a large inversion including nine genes (nad3-trnG-atp6-atp8-trnD-trnK-cox2-trnL2-cox1), which was only congruent with the species in the genus Nasonia. Moreover, we found that one shuffling changes (trnD and trnK) happened in the species E. obtusiclava but not in another species E. formosa within the same genus, of which such shuffling within the same genus at this region was only detected in Polisters within Insecta. Phylogenetic analysis displayed that different data matrix (13PCG+ 2 rRNA or 13 PCG) and inference methods (BI or ML) indicate the identical topology with high nodal supports that Aphelinidae formed a sister group with (Trichogrammatidae + Aganoidae) and the monophyly of Pteramalidae. Our results also indicated the validity of assembling and feasibility of next-generation technology to obtain the mitochondrial genomes of parasitic Hymenoptera.

Looking for variable molecular markers in the chestnut gall wasp Dryocosmus kuriphilus: first comparison across genes

The quick spread of the chestnut gall wasp Dryocosmus kuriphilus in Europe constitutes an outstanding example of recent human-aided biological invasion with dramatic economic losses. We screened for the first time a set of five nuclear and mitochondrial genes from D. kuriphilus collected in the Iberian Peninsula, and compared the sequences with those available from the native and invasive range of the species. We found no genetic variability in Iberia in none of the five genes, moreover, the three genes compared with other European samples showed no variability either. We recorded four cytochrome b haplotypes in Europe; one was genuine mitochondrial DNA and the rest nuclear copies of mitDNA (numts), what stresses the need of careful in silico analyses. The numts formed a separate cluster in the gene tree and at least two of them might be orthologous, what suggests that the invasion might have started with more than one individual. Our results point at a low initial population size in Europe followed by a quick population growth. Future studies assessing the expansion of this pest should include a large number of sampling sites and use powerful nuclear markers (e. g. Single Nucleotide Polymorphisms) to detect genetic variability.

Mitochondrial genome diversity and population structure of two western honey bee subspecies in the Republic of South Africa

Apis mellifera capensis Eschscholtz and A.m. scutellata Lepeletier are subspecies of western honey bees that are indigenous to the Republic of South Africa (RSA). Both subspecies have invasive potential and are organisms of concern for areas outside their native range, though they are important bees to beekeepers, agriculture, and the environment where they are native. The aim of the present study was to examine genetic differentiation among these subspecies and estimate their phylogenetic relationships using complete mitochondrial genomes sequences. We used 25 individuals that were either assigned to one of the subspecies or designated hybrids using morphometric analyses. Phylogenetic analyses of mitogenome sequences by maximum likelihood (ML) and Bayesian inference identified a monophyletic RSA clade, subdivided into two clades. A haplotype network was consistent with the phylogenetic trees. However, members of both subspecies occurred in both clades, indicating that A.m. capensis and A.m. scutellata are neither reciprocally monophyletic nor do they exhibit paraphyly with one subspecies nested within the other subspecies. Furthermore, no mitogenomic features were diagnostic to either subspecies. All bees analyzed from the RSA expressed a substantial level of haplotype diversity (most samples had unique haplotypes) but limited nucleotide diversity. The number of variable codons across protein-coding genes (PCGs) differed among loci, with CO3 exhibiting the most variation and ATP6 the least.

The molecular evolutionary dynamics of oxidative phosphorylation (OXPHOS) genes in Hymenoptera

BACKGROUND

The primary energy-producing pathway in eukaryotic cells, the oxidative phosphorylation (OXPHOS) system, comprises proteins encoded by both mitochondrial and nuclear genes. To maintain the function of the OXPHOS system, the pattern of substitutions in mitochondrial and nuclear genes may not be completely independent. It has been suggested that slightly deleterious substitutions in mitochondrial genes are compensated by substitutions in the interacting nuclear genes due to positive selection. Among the four largest insect orders, Coleoptera (beetles), Hymenoptera (sawflies, wasps, ants, and bees), Diptera (midges, mosquitoes, and flies) and Lepidoptera (moths and butterflies), the mitochondrial genes of Hymenoptera exhibit an exceptionally high amino acid substitution rate while the evolution of nuclear OXPHOS genes is largely unknown. Therefore, Hymenoptera is an excellent model group for testing the hypothesis of positive selection driving the substitution rate of nuclear OXPHOS genes. In this study, we report the evolutionary rate of OXPHOS genes in Hymenoptera and test for evidence of positive selection in nuclear OXPHOS genes of Hymenoptera.

RESULTS

Our analyses revealed that the amino acid substitution rate of mitochondrial and nuclear OXPHOS genes in Hymenoptera is higher than that in other studied insect orders. In contrast, the amino acid substitution rate of non-OXPHOS genes in Hymenoptera is lower than the rate in other insect orders. Overall, we found the dN/dS ratio of the nuclear OXPHOS genes to be higher in Hymenoptera than in other insect orders. However, nuclear OXPHOS genes with high dN/dS ratio did not always exhibit a high amino acid substitution rate. Using branch-site and site model tests, we identified various codon sites that evolved under positive selection in nuclear OXPHOS genes.

CONCLUSIONS

Our results showed that nuclear OXPHOS genes in Hymenoptera are evolving faster than the genes in other three insect orders. The branch test suggested that while some nuclear OXPHOS genes in Hymenoptera show a signature of positive selection, the pattern is not consistent across all nuclear OXPHOS genes. As only few codon sites were under positive selection, we suggested that positive selection might not be the only factor contributing to the rapid evolution of nuclear OXPHOS genes in Hymenoptera.

Counting animal species with DNA barcodes: Canadian insects

Recent estimates suggest that the global insect fauna includes fewer than six million species, but this projection is very uncertain because taxonomic work has been limited on some highly diverse groups. Validation of current estimates minimally requires the investigation of all lineages that are diverse enough to have a substantial impact on the final species count. This study represents a first step in this direction; it employs DNA barcoding to evaluate patterns of species richness in 27 orders of Canadian insects. The analysis of over one million specimens revealed species counts congruent with earlier results for most orders. However, Diptera and Hymenoptera were unexpectedly diverse, representing two-thirds of the 46 937 barcode index numbers (=species) detected. Correspondence checks between known species and barcoded taxa showed that sampling was incomplete, a result confirmed by extrapolations from the barcode results which suggest the occurrence of at least 94 000 species of insects in Canada, a near doubling from the prior estimate of 54 000 species. One dipteran family, the Cecidomyiidae, was extraordinarily diverse with an estimated 16 000 species, a 10-fold increase from its predicted diversity. If Canada possesses about 1% of the global fauna, as it does for known taxa, the results of this study suggest the presence of 10 million insect species with about 1.8 million of these taxa in the Cecidomyiidae. If so, the global species count for this fly family may exceed the combined total for all 142 beetle families. If extended to more geographical regions and to all hyperdiverse groups, DNA barcoding can rapidly resolve the current uncertainty surrounding a species count for the animal kingdom. A newly detailed understanding of species diversity may illuminate processes important in speciation, as suggested by the discovery that the most diverse insect lineages in Canada employ an unusual mode of reproduction, haplodiploidy.

This article is part of the themed issue ‘From DNA barcodes to biomes’.

Evolutionarily recent, insertional fission of mitochondrial cox2 into complementary genes in bilaterian Metazoa

Background

Mitochondrial genomes (mtDNA) of multicellular animals (Metazoa) with bilateral symmetry (Bilateria) are compact and usually carry 13 protein-coding genes for subunits of three respiratory complexes and ATP synthase. However, occasionally reported exceptions to this typical mtDNA organization prompted speculation that, as in protists and plants, some bilaterian mitogenomes may continue to lose their canonical genes, or may even acquire new genes. To shed more light on this phenomenon, a PCR-based screen was conducted to assess fast-evolving mtDNAs of apocritan Hymenoptera (Arthropoda, Insecta) for genomic rearrangements that might be associated with the modification of mitochondrial gene content.

Results

Sequencing of segmental inversions, identified in the screen, revealed that the cytochrome oxidase subunit II gene (cox2) of Campsomeris (Dielis) (Scoliidae) was split into two genes coding for COXIIA and COXIIB. The COXII-derived complementary polypeptides apparently form a heterodimer, have reduced hydrophobicity compared with the majority of mitogenome-encoded COX subunits, and one of them, COXIIB, features increased content of Cys residues. Analogous cox2 fragmentation is known only in two clades of protists (chlorophycean algae and alveolates), where it has been associated with piecewise relocation of this gene into the nucleus. In Campsomeris mtDNA, cox2a and cox2b loci are separated by a 3-kb large cluster of several antiparallel overlapping ORFs, one of which, qnu, seems to encode a nuclease that may have played a role in cox2 fission.

Conclusions

Although discontinuous mitochondrial protein genes encoding fragmented, complementary polypeptides are known in protists and some plants, split cox2 of Campsomeris is the first case of such a gene arrangement found in animals. The reported data also indicate that bilaterian animal mitogenomes may be carrying lineage-specific genes more often than previously thought, and suggest a homing endonuclease-based mechanism for insertional mitochondrial gene fission.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-017-3626-5) contains supplementary material, which is available to authorized users.

Capturing the Phylogeny of Holometabola with Mitochondrial Genome Data and Bayesian Site-Heterogeneous Mixture Models

After decades of debate, a mostly satisfactory resolution of relationships among the 11 recognized holometabolan orders of insects has been reached based on nuclear genes, resolving one of the most substantial branches of the tree-of-life, but the relationships are still not well established with mitochondrial genome data. The main reasons have been the absence of sufficient data in several orders and lack of appropriate phylogenetic methods that avoid the systematic errors from compositional and mutational biases in insect mitochondrial genomes. In this study, we assembled the richest taxon sampling of Holometabola to date (199 species in 11 orders), and analyzed both nucleotide and amino acid data sets using several methods. We find the standard Bayesian inference and maximum-likelihood analyses were strongly affected by systematic biases, but the site-heterogeneous mixture model implemented in PhyloBayes avoided the false grouping of unrelated taxa exhibiting similar base composition and accelerated evolutionary rate. The inclusion of rRNA genes and removal of fast-evolving sites with the observed variability sorting method for identifying sites deviating from the mean rates improved the phylogenetic inferences under a site-heterogeneous model, correctly recovering most deep branches of the Holometabola phylogeny. We suggest that the use of mitochondrial genome data for resolving deep phylogenetic relationships requires an assessment of the potential impact of substitutional saturation and compositional biases through data deletion strategies and by using site-heterogeneous mixture models. Our study suggests a practical approach for how to use densely sampled mitochondrial genome data in phylogenetic analyses.

Comparative and phylogenetic analysis of the mitochondrial genomes in basal hymenopterans

The Symphyta is traditionally accepted as a paraphyletic group located in a basal position of the order Hymenoptera. Herein, we conducted a comparative analysis of the mitochondrial genomes in the Symphyta by describing two newly sequenced ones, from Trichiosoma anthracinum, representing the first mitochondrial genome in family Cimbicidae, and Asiemphytus rufocephalus, from family Tenthredinidae. The sequenced lengths of these two mitochondrial genomes were 15,392 and 14,864 bp, respectively. Within the sequenced region, trnC and trnY were rearranged to the upstream of trnI-nad2 in T. anthracinum, while in A. rufocephalus all sequenced genes were arranged in the putative insect ancestral gene arrangement. Rearrangement of the tRNA genes is common in the Symphyta. The rearranged genes are mainly from trnL1 and two tRNA clusters of trnI-trnQ-trnM and trnW-trnC-trnY. The mitochondrial genomes of Symphyta show a biased usage of A and T rather than G and C. Protein-coding genes in Symphyta species show a lower evolutionary rate than those of Apocrita. The Ka/Ks ratios were all less than 1, indicating purifying selection of Symphyta species. Phylogenetic analyses supported the paraphyly and basal position of Symphyta in Hymenoptera. The well-supported phylogenetic relationship in the study is Tenthredinoidea + (Cephoidea + (Orussoidea + Apocrita)).

Complete mitochondrial genome of Camponotus atrox (Hymenoptera: Formicidae): a new tRNA arrangement in Hymenoptera

We sequenced the complete mitochondrial (mt) genome of Camponotus atrox (Hymenoptera: Formicidae), which is only distributed in Korea. The genome was 16 540 bp in size and contained typical sets of genes (13 protein-coding genes, 22 tRNAs, and 2 rRNAs). The C. atrox A+T-rich region, at 1402 bp, was the longest of all sequenced ant genomes and was composed of an identical tandem repeat consisting of six 100-bp copies and one 96-bp copy. A total of 315 bp of intergenic spacer sequence was spread over 23 regions. An alignment of the spacer sequences in ants was largely feasible among congeneric species, and there was substantial sequence divergence, indicating their potential use as molecular markers for congeneric species. The A/T contents at the first and second codon positions of protein-coding genes (PCGs) were similar for ant species, including C. atrox (73.9% vs. 72.3%, on average). With increased taxon sampling among hymenopteran superfamilies, differences in the divergence rates (i.e., the non-synonymous substitution rates) between the suborders Symphyta and Apocrita were detected, consistent with previous results. The C. atrox mt genome had a unique gene arrangement, trnI-trnM-trnQ, at the A+T-rich region and ND2 junction (underline indicates inverted gene). This may have originated from a tandem duplication of trnM-trnI, resulting in trnM-trnI-trnM-trnI-trnQ, and the subsequent loss of the first trnM and second trnI, resulting in trnI-trnM-trnQ.

Phylogenetic relationships of Hemiptera inferred from mitochondrial and nuclear genes

Here, we reconstructed the Hemiptera phylogeny based on the expanded mitochondrial protein-coding genes and the nuclear 18S rRNA gene, separately. The differential rates of change across lineages may associate with long-branch attraction (LBA) effect and result in conflicting estimates of phylogeny from different types of data. To reduce the potential effects of systematic biases on inferences of topology, various data coding schemes, site removal method, and different algorithms were utilized in phylogenetic reconstruction. We show that the outgroups Phthiraptera, Thysanoptera, and the ingroup Sternorrhyncha share similar base composition, and exhibit "long branches" relative to other hemipterans. Thus, the long-branch attraction between these groups is suspected to cause the failure of recovering Hemiptera under the homogeneous model. In contrast, a monophyletic Hemiptera is supported when heterogeneous model is utilized in the analysis. Although higher level phylogenetic relationships within Hemiptera remain to be answered, consensus between analyses is beginning to converge on a stable phylogeny.

Higher-level phylogeny of the Hymenoptera inferred from mitochondrial genomes

Higher-level hymenopteran relationships remain unresolved in both morphological and molecular analyses. In this study, we present the most comprehensive analyses of hymenopteran relationships based on 48 mitochondrial (mt) genomes. One complete and two nearly complete mt genomes representing three hymenopteran superfamilies were newly sequenced. We assessed the influence of inclusion/exclusion of 3rd codon positions, alignment approaches, partition schemes and phylogenetic approaches on topology and nodal support within the Hymenoptera. The results showed that the topologies were sensitive to the variation of dataset and analytical approach. However, some robust and highly supported relationships were recovered: the Ichneumonomorpha was monophyletic; the Trigonalyoidea+Megalyroidea and the Diaprioidea+Chalcidoidea were consistently recovered; the Cynipoidea was generally recovered as the sister group to the Diaprioidea+Chalcidoidea. In addition, the monophyletic Aculeata and Proctotrupomorpha were recovered in some analyses. Several gene rearrangements were detected in each of the three newly sequenced mt genomes. Specifically, the Ibalia leucospoides mt genome harbors a large inversion of a gene block from trnE to trnS2. Inverted, duplicated A+T rich regions were detected in the Ibalia leucospoides mt genome, which probably played an important role during the formation of the large gene block inversion via recombination.

Is it an ant or a butterfly? Convergent evolution in the mitochondrial gene order of Hymenoptera and Lepidoptera

Insect mitochondrial genomes (mtDNA) are usually double helical and circular molecules containing 37 genes that are encoded on both strands. The arrangement of the genes is not constant for all species, and produces distinct gene orders (GOs) that have proven to be diagnostic in defining clades at different taxonomic levels. In general, it is believed that distinct taxa have a very low chance of sharing identically arranged GOs. However, examples of identical, homoplastic local rearrangements occurring in distinct taxa do exist. In this study, we sequenced the complete mtDNAs of the ants Formica fusca and Myrmica scabrinodis (Formicidae, Hymenoptera) and compared their GOs with those of other Insecta. The GO of F. fusca was found to be identical to the GO of Dytrisia (the largest clade of Lepidoptera). This finding is the first documented case of an identical GO shared by distinct groups of Insecta, and it is the oldest known event of GO convergent evolution in animals. Both Hymenoptera and Lepidoptera acquired this GO early in their evolution. Using a phylogenetic approach combined with new bioinformatic tools, the chronological order of the evolutionary events that produced the diversity of the hymenopteran GOs was determined. Additionally, new local homoplastic rearrangements shared by distinct groups of insects were identified. Our study showed that local and global homoplasies affecting the insect GOs are more widespread than previously thought. Homoplastic GOs can still be useful for characterizing the various clades, provided that they are appropriately considered in a phylogenetic and taxonomic context.

Complete mitochondrial genomes of Ceratobaeus sp. and Idris sp. (Hymenoptera: Scelionidae): shared gene rearrangements as potential phylogenetic markers at the tribal level

We sequenced the complete mitochondrial genomes of two sceliond taxa (Ceratobaeus sp. and Idris sp.). An atypical tRNA-Arg which lacks a D-stem was identified in both taxa, and represents a potentially derived character of sceliond wasps. A number of tRNA genes have rearranged in the two mitochondrial genomes compared with the ancestral organization. Some of these derived genome organizations are shared, and thus have much potential as phylogenetic markers at the tribal level in the subfamily Scelioninae. We test the influence of third codon inclusion/exclusion, alignment methods and partition schemes on the reconstruction of phylogenetic relationships. The results show that inclusion of third codon positions does not appear to be problematic when investigating the phylogeny of closely related taxa. Muscle and PartitionFinder schemes significantly improve the likelihood scores.

Evolutionary dynamics of the mitochondrial genome in the evaniomorpha (hymenoptera)—a group with an intermediate rate of gene rearrangement

We determined the complete mitochondrial (mt) genomes of three evaniomorph species, Ceraphron sp. (Ceraphronoidea), Gasteruption sp. (Evanioidea), and Orthogonalys pulchella (Trigonalyoidea) as well as the nearly complete mt genome from another evaniomorph species, Megalyra sp. (Megalyroidea). Each of them possesses dramatic gene rearrangements, including protein-coding or rRNA genes. Gene inversions were identified in all of these mt genomes; for example, the two rRNA genes have inverted and moved into the nad2-cox1 junction in the Megalyra sp. mt genome. In addition, we found two copies of a 10-bp complementary repeat at the beginning of rrnS and at the end of trnL(2) in the Gasteruption sp. mt genome, consistent with recombination as the possible mechanism for gene inversion and long-range movement. Although each of the genomes contains a number of repeats of varying size, there was no consistent association of the size or number of repeats with the extent or type of gene rearrangement. The breakpoint distance analysis showed the Evaniomorpha has an intermediate rate of gene rearrangement. Sequence-based phylogenetic analyses of 13 protein-coding and 2 rRNA genes in 22 hymenopteran taxa recovered a paraphyletic Evaniomorpha with the Aculeata nested within it. Within the Evaniomorpha, our analyses confirmed the Trigonalyoidea + Megalyroidea as the sister group to the Aculeata and recovered a novel clade, Ceraphronoidea + Evanioidea. In contrast to previous hymenopteran phylogenetic studies, the internal relationships of the Evaniomorpha were highly supported and robust to the variation of alignment approach and phylogenetic inference approach.

The mitochondrial genome of the Vespa bicolor Fabricius (Hymenoptera: Vespidae: Vespinae)

Here, we report the first representative mitochondrial genome of the subfamily Vespinae (Hymenopterea: Vespidae) from the Vespa bicolor Fabricius (GenBank accession No. KJ735511). Nearly complete mitochondrial genome was sequenced with a length of 16,937 bp, including 13 protein-coding genes, 2 rRNA and 20 tRNA genes, as well as a portion of A+T-rich region. Two tRNA genes, i.e. trnI and trnY, were failed to sequence, which were presumed to be located within a region between A+T-rich region and trnM-trnQ-nad2. In the V. bicolor mitochondrial genome, at least three tRNA genes were rearranged. trnY was rearranged to the unsequenced region between A+T-rich region and trnM-trnQ-nad2. trnL1 was rearranged from a location between nad1 and rrnL to the upstream of nad1 gene. trnS1 and trnE were shuffled in the tRNA cluster of trnA-trnR-trnN-trnS1-trnE-trnF. Our study showed that the mitochondrial genomes between Vespinae and Polistinae shared more arrangement pattern than that between Vespinae and Eumeninae.

The mitochondrial genome of the leaf-cutter ant Atta laevigata: a mitogenome with a large number of intergenic spacers

In this paper we describe the nearly complete mitochondrial genome of the leaf-cutter ant Atta laevigata, assembled using transcriptomic libraries from Sanger and Illumina next generation sequencing (NGS), and PCR products. This mitogenome was found to be very large (18,729 bp), given the presence of 30 non-coding intergenic spacers (IGS) spanning 3,808 bp. A portion of the putative control region remained unsequenced. The gene content and organization correspond to that inferred for the ancestral pancrustacea, except for two tRNA gene rearrangements that have been described previously in other ants. The IGS were highly variable in length and dispersed through the mitogenome. This pattern was also found for the other hymenopterans in particular for the monophyletic Apocrita. These spacers with unknown function may be valuable for characterizing genome evolution and distinguishing closely related species and individuals. NGS provided better coverage than Sanger sequencing, especially for tRNA and ribosomal subunit genes, thus facilitating efforts to fill in sequence gaps. The results obtained showed that data from transcriptomic libraries contain valuable information for assembling mitogenomes. The present data also provide a source of molecular markers that will be very important for improving our understanding of genomic evolutionary processes and phylogenetic relationships among hymenopterans.

Mitochondrial-nuclear interactions: compensatory evolution or variable functional constraint among vertebrate oxidative phosphorylation genes?

Oxidative phosphorylation (OXPHOS), the major energy-producing pathway in aerobic organisms, includes protein subunits encoded by both mitochondrial (mt) and nuclear (nu) genomes. How these independent genomes have coevolved is a long-standing question in evolutionary biology. Although mt genes evolve faster than most nu genes, maintenance of OXPHOS structural stability and functional efficiency may involve correlated evolution of mt and nu OXPHOS genes. The nu OXPHOS genes might be predicted to exhibit accelerated evolutionary rates to accommodate the elevated substitution rates of mt OXPHOS subunits with which they interact. Evolutionary rates of nu OXPHOS genes should, therefore, be higher than that of nu genes that are not involved in OXPHOS (nu non-OXPHOS). We tested the compensatory evolution hypothesis by comparing the evolutionary rates (synonymous substitution rate dS and nonsynonymous substitution rate dN) among 13 mt OXPHOS genes, 60 nu OXPHOS genes, and 77 nu non-OXPHOS genes in vertebrates (7 fish and 40 mammal species). The results from a combined analysis of all OXPHOS subunits fit the predictions of the hypothesis. However, results from two OXPHOS complexes did not fit this pattern when analyzed separately. We found that the d(N) of nu OXPHOS genes for "core" subunits (those involved in the major catalytic activity) was lower than that of "noncore" subunits, whereas there was no significant difference in d(N) between genes for nu non-OXPHOS and core subunits. This latter finding suggests that compensatory changes play a minor role in the evolution of OXPHOS genes and that the observed accelerated nu substitution rates are due largely to reduced functional constraint on noncore subunits.

Rearrangement of the nad1 gene in Pristaulacus compressus (Spinola) (Hymenoptera: Evanioidea: Aulacidae) mitochondrial genome

The mitochondrial genome of the Pristaulacus compressus (Spinola, 1808) (Hymenoptera: Aulacidae) (GenBank accession No. KF500406) is reported in this study. This is the first sequenced mitochondrial genome from the family Aulacidae of the order Hymenoptera. The length of this mitochondrial genome is 15,563 bp with an A + T content of 84%, including 13 protein-coding, 2 rRNA and 22 tRNA gene, and an A + T-rich region (Table 1). Three tRNA and one protein-coding genes were rearranged in the P. compressus mitochondrial genome, in which, the trnY was inverted, while the trnQ was shuffled to the downstream of tRNA cluster trnI-trnQ-trnM. The trnS1 was translocated to the downstream of the A + T-rich region together with the protein-coding gene nad1. The gene arrangement pattern of this mitochondrial genome is new to the Hymenoptera. All protein-coding genes start with ATN start codon. Ten protein-coding genes stop with termination codon TAA, whereas one protein-coding gene uses incomplete stop codon TA and two use T. The A + T-region is located between rrnS and trnS1 with a length of 780 bp.