Characteristics of the nuclear (18S, 5.8S, 28S and 5S) and mitochondrial (12S and 16S) rRNA genes of Apis mellifera (Insecta: Hymenoptera): structure, organization, and retrotransposable elements

Comparative Mitogenomics Reveals Cryptic Species in Sillago ingenuua McKay, 1985 (Perciformes: Sillaginidae)

It is unreliable to identify marine fishes only by external morphological features. Species misidentification brings great challenges to fishery research, resource monitoring and ecomanagement. Sillago ingenuua is an important part of commercial marine fishes, and in which, the morphological differences between different groups are not obvious. Here, we compared different geographical groups of S. ingenuua which were collected from Xiamen, Dongshan, Keelung, Songkhla and Java. The results showed that all samples of S. ingenuua were similar in external morphological characteristics and the shape of the swim bladder, but there were two distinctive lineages which were flagged as cryptic species based on DNA barcoding. The comparative mitogenomic results showed that S. ingenuua A and S. ingenuua B were identical in structural organization and gene arrangement. Their nucleotide composition and codon usage were also similar. A phylogenetic analysis was performed based on 13 concatenated PCGs from eight Sillago species. The results showed that the genetic distance between S. ingenuua A and S. ingenuua B was large (D = 0.069), and this genetic distance was large enough to reveal that S. ingenuua A and S. ingenuua B might be different species.

Mediating a host cell signaling pathway linked to overwinter mortality offers a promising therapeutic approach for improving bee health

INTRODUCTION

Honey bees provides valuable pollination services for world food crops and wild flowering plants which are habitats of many animal species and remove carbon dioxide from the atmosphere, a powerful tool in the fight against climate change. Nevertheless, the honey bee population has been declining and the majority of colony losses occur during the winter.

OBJECTIVES

The goal of this study was to understand the mechanisms underlying overwinter colony losses and develop novel therapeutic strategies for improving bee health.

METHODS

First, pathogen prevalence in overwintering bees were screened between 2015 and 2018. Second, RNA sequencing (RNA-Seq) for transcriptional profiling of overwintering honey bees was conducted and qRT-PCR was performed to confirm the results of the differential expression of selected genes. Lastly, laboratory bioassays were conducted to measure the effects of cold challenges on bee survivorship and stress responses and to assess the effect of a novel medication for alleviating cold stress in honey bees.

RESULTS

We identified that sirtuin signaling pathway is the most significantly enriched pathway among the down-regulated differentially expressed genes (DEGs) in overwintering diseased bees. Moreover, we showed that the expression of SIRT1 gene, a major sirtuin that regulates energy and immune metabolism, was significantly downregulated in bees merely exposed to cold challenges, linking cold stress with altered gene expression of SIRT1. Furthermore, we demonstrated that activation of SIRT1 gene expression by SRT1720, an activator of SIRT1 expression, could improve the physiology and extend the lifespan of cold-stressed bees.

CONCLUSION

Our study suggests that increased energy consumption of overwintering bees for maintaining hive temperature reduces the allocation of energy toward immune functions, thus making the overwintering bees more susceptible to disease infections and leading to high winter colony losses. The novel information gained from this study provides a promising avenue for the development of therapeutic strategies for mitigating colony losses, both overwinter and annually.

Where Eriophyoidea (Acariformes) Belong in the Tree of Life

Over the past century and a half, the taxonomic placement of Eriophyoidea has been in flux. For much of this period, this group has been treated as a subtaxon within Trombidiformes. However, the vast majority of recent phylogenetic analyses, including almost all phylogenomic analyses, place this group outside Trombidiformes. The few studies that still place Eriophyoidea within Trombidiformes are likely to be biased by incomplete taxon/gene sampling, long branch attraction, the omission of RNA secondary structure in sequence alignment, and the inclusion of hypervariable expansion-contraction rRNA regions. Based on the agreement among a number of independent analyses that use a range of different datasets (morphology; multiple genes; mitochondrial/whole genomes), Eriophyoidea are almost certain to be closely related to Nematalycidae, a family of vermiform mites within Endeostigmata, a basal acariform grade. Much of the morphological evidence in support of this relationship was apparent after the discovery of Nematalycidae in the middle of the 20th century. However, this evidence has largely been disregarded until very recently, perhaps because of overconfidence in the placement of Eriophyoidea within Trombidiformes. Here, we briefly review and identify a number of biases, both molecular- and morphology-based, that can lead to erroneous reconstructions of the position of Eriophyoidea in the tree of life.

Molecular Apomorphies in the Secondary and Tertiary Structures of Length-Variable Regions (LVRs) of 18S rRNA Shed Light on the Systematic Position of the Family Thaumastellidae (Hemiptera: Heteroptera: Pentatomoidea)

The SSU nrDNA, a small subunit of the nuclear ribosomal DNA (coding 18S rRNA), is one of the most frequently sequenced genes in molecular studies in Hexapoda. In insects, including true bugs (Hemiptera: Heteroptera), only its primary structures (i.e., aligned sequences) are predominantly used in phylogenetic reconstructions. It is known that including RNA secondary structures in the alignment procedure is essential for improving accuracy and robustness in phylogenetic tree reconstruction. Moreover, local plasticity in rRNAs might impact their tertiary structures and corresponding functions. To determine the systematic position of Thaumastellidae within the superfamily Pentatomoidea, the secondary and-for the first time among all Hexapoda-tertiary structures of 18S rRNAs in twelve pentatomoid families were compared and analysed. Results indicate that the shapes of the secondary and tertiary structures of the length-variable regions (LVRs) in the 18S rRNA are phylogenetically highly informative. Based on these results, it is suggested that the Thaumastellidae is maintained as an independent family within the superfamily Pentatomoidea, rather than as a part of the family Cydnidae. Moreover, the analyses indicate a close relationship between Sehirinae and Parastrachiidae, expressed in morpho-molecular synapomorphies in the predicted secondary and tertiary structures of the length-variable region L (LVR L).

First detection of Ixodiphagus hookeri (Hymenoptera: Encyrtidae) in Ixodes ricinus ticks (Acari: Ixodidae) from multiple locations in Hungary

The parasitoid wasp, Ixodiphagus hookeri (Hymenoptera: Encyrtidae), is the natural enemy of a wide range of hard and soft tick species. While these encyrtid wasps are supposed to be distributed worldwide, only a few studies report on their actual distribution around the globe. Within a shotgun sequencing-based metagenome analysis, the occurrence of I. hookeri was screened at multiple Ixodes ricinus (Acari: Ixodidae) tick sampling points in Hungary to contribute to the assessment of the distribution patterns of the parasitoid wasps in Central Europe. To our knowledge, the first report of the species in Hungary and the description of the southernmost I. hookeri associated geoposition in Central Europe took place within our study. I. hookeri infested I. ricinus nymphs were detected at five sampling points in Hungary. The results show that the exact distribution range of I. hookeri is still barely studied. At the same time, unprecedented public health issues being brought about by climate change might require steps toward the exploitation of the tick biocontrol potential and as an ecological bioindicator role of the parasitoid wasp in the future.

The complete mitochondrial genome and novel gene arrangement in Nesodiprion zhejiangensis Zhou & Xiao (Hymenoptera: Diprionidae)

The complete mitochondrial genome (mitogenome) of the sawfly, Nesodiprion zhejiangensis Zhou & Xiao, was sequenced, assembled, and deposited in GenBank (Accession Number: OM501121). The 15,660 bp N. zhejiangensis mitogenome encodes for 2 ribosomal RNAs (rrnL and rrnS), 22 transfer RNAs (tRNAs), 13 protein-coding genes (PCGs), and an AT-rich region of 450 bp in length. The nucleotide composition is biased toward adenine and thymine (A + T = 81.8%). Each PCG is initiated by an ATN codon, except for cox2, which starts with a TTG. Of 13 PCGs, 9 have a TAA termination codon, while the remainder terminate with a TAG or a single T. All tRNAs have the classic cloverleaf structure, except for the dihydrouridine (DHU) arm of tRNA^val, which forms a simple loop. There are 49 helices belonging to 6 domains in rrnL and 30 helices belonging to 4 domains in rrnS. In comparison to the ancestral architecture, N. zhejiangensis has the most rearranged mitogenome in Symphyta, in which rearrangement events of local inversion and transposition are identified in three gene clusters. Specifically, the main hotspot of gene rearrangement occurred between rrnS and trnY, and rearranged from rrnS-(AT-rich region)-I-Q-M-nd2-W-C-Y to rrnS-Q-W-C-nd2-I-M-(AT-rich region)-Y, involving a local inversion event of a large gene cluster and transposition events of some tRNAs. Transposition of trnA and trnR (rearranged from A-R to R-A) was observed at the nd3-nd5 gene junction while shuffling of trnP and trnT (rearranged from T-P to P-T) occurred at the nd4l-nd6 gene junction. While illegitimate inter-mtDNA recombination might explain the opposite orientations of transcription between rrnS and trnY, transposition events of tRNA in some gene blocks can be accounted for by the tandem duplication/random loss (TDRL) model. Our phylogenetic analysis suggests that N. zhejiangensis is closely related to congeneric species N. biremis and N. japonicus, which together form a sister lineage with the European pine sawfly, Neodiprion sertifer.

Eighteen mitochondrial genomes of Syrphidae (Insecta: Diptera: Brachycera) with a phylogenetic analysis of Muscomorpha

In this study, 18 mitochondrial genomes (mitogenomes) of Syrphidae were sequenced. These mitogenomes ranged from 15,648 to 16,405 bp and contained 37 genes that were similar to those from other Syrphidae species. Most protein-coding genes (PCGs) started with a standard ATN codon and ended with TAA/G. All transfer RNAs (tRNAs) could be folded into the cloverleaf secondary structure except tRNA-Ser (AGN), which lacks a dihydrouridine arm. The secondary structures of ribosomal RNAs (rRNAs) were predicted. Six domains (III is absent in arthropods) and 44 helices were included in the 16S rRNA, and three domains and 24 helices were included in the 12S rRNA. We found three conserved fragments in all syrphid mitogenomes. Phylogenetic analyses were performed based on the nucleotide data of 13 PCGs and two rRNAs from 76 Muscomorpha and three outgroup species. In results the paraphyly of Aschiza and Schizophora were supported, the Acalyptratae was also paraphyletic but the relationships of its superfamilies were difficult to determine, the monophyly of Calyptratea was supported with the relationships of Oestroidea and Muscoidea need to be further reconsidered. Within Syrphidae the monophyly of family level was supported, the Syrphinae were clustered into one branch, while the paraphyly of Eristalinae was still well supported.

One-way ticket to the blue: A large-scale, dated phylogeny revealed asymmetric land-to-water transitions in acariform mites (Acari: Acariformes)

Acariform mites are an ancient and megadiverse lineage that may have experienced a complex pattern of invasions into terrestrial and aquatic habitats. These among-realm transitions may relate to periods of turmoil in Earth's history or be simply results of uneven biodiversity patterns across habitats. Here, we inferred a dated, representative acariform phylogeny (five genes, 9,200 bp aligned, 367 terminals belonging to 150 ingroup plus 15 outgroup families, 23 fossil calibration points) which was used to infer transitions between marine/freshwater/terrestrial habitats. We detected four unambiguous transitions from terrestrial to freshwater habitats (Hydrozetes, Naiadacarus, Fusohericia, Afronothrus, Homocaligus); one from freshwater to marine (Pontarachnidae), and four from marine to brackish or freshwater transitions (all among Halacaridae: Acarothrix; Halacarellus petiti; Copidognathus sp.; clade Limnohalacarus + Soldanellonyx + Porohalacarus + Porolohmannella). One transition to the sea was inferred ambiguously with respect to the ancestor being either terrestrial or freshwater (Hyadesiidae), and another must be most carefully examined by adding potential related taxa (Selenoribatidae + Fortuyniidae). Finally, we inferred a single, remarkable transition from aquatic to terrestrial habitats involving early evolution of the large and ecologically diverse lineage: the ancestor of the Halacaridae + Parasitengona clade was probably freshwater given our dataset, thus making terrestrial Parasitengona secondarily terrestrial. Overall, our results suggested a strong asymmetry in environmental transitions: the majority occurred from terrestrial to aquatic habitats. This asymmetry is probably linked to mites' biological properties and uneven biodiversity patterns across habitats rather than Earth's geological history. Since the land holds more acariform diversity than water habitats, a shift from the former is more likely than from the latter. We inferred the following relationships: alicid endeostigmatid + eriophyoid (Alycidae, (Nanorchestidae, (Nematalycidae, Eriophyoidea))) being sister group to the remaining Acariformes: (proteonematalycid Endeostigmata, alicorhagiid Endeostigmata, Trombidiformes, Oribatida (including Astigmata)). Trombidiform relationships had several novel rearrangements: (i) traditional Eupodina lacked support for the inclusion of Bdelloidea; (ii) Teneriffidae, traditionally placed among Anystina, was consistently recovered in a clade including Heterostigmata in Eleutherengona; (iii) several lineages, such as Adamystidae, Paratydeidae, Caeculidae and Erythracaridae, were recovered in a large clade along other Anystina and Eleutherengona, suggesting single origins of several fundamental character states, such as the reduction of the cheliceral fixed digit and development of the palpal thumb-claw complex.

rRNA expansion segment 7 in eukaryotes: from Signature Fold to tentacles

The ribosomal core is universally conserved across the tree of life. However, eukaryotic ribosomes contain diverse rRNA expansion segments (ESs) on their surfaces. Sites of ES insertions are predicted from sites of insertion of micro-ESs in archaea. Expansion segment 7 (ES7) is one of the most diverse regions of the ribosome, emanating from a short stem loop and ranging to over 750 nucleotides in mammals. We present secondary and full-atom 3D structures of ES7 from species spanning eukaryotic diversity. Our results are based on experimental 3D structures, the accretion model of ribosomal evolution, phylogenetic relationships, multiple sequence alignments, RNA folding algorithms and 3D modeling by RNAComposer. ES7 contains a distinct motif, the 'ES7 Signature Fold', which is generally invariant in 2D topology and 3D structure in all eukaryotic ribosomes. We establish a model in which ES7 developed over evolution through a series of elementary and recursive growth events. The data are sufficient to support an atomic-level accretion path for rRNA growth. The non-monophyletic distribution of some ES7 features across the phylogeny suggests acquisition via convergent processes. And finally, illustrating the power of our approach, we constructed the 2D and 3D structure of the entire LSU rRNA of Mus musculus.

Paromomycin Reduces Vairimorpha (Nosema) ceranae Infection in Honey Bees but Perturbs Microbiome Levels and Midgut Cell Function

Paromomycin is a naturally occurring aminoglycoside antibiotic that has effects on both prokaryotic and eukaryotic microbes. However, previous reports have indicated that it has little effect on microsporidia, including Vairimorpha (Nosema) ceranae, in cell culture models. V. ceranae is one of a number of microsporidia species that cause disease in honey bees and substantial efforts to find new treatment strategies for bees that are infected with these pathogens are ongoing. When testing compounds for potential activity against V. ceranae in whole organisms, we found that paromomycin reduces the infection intensity of this parasite. Critically, the necessary doses of paromomycin have high activity against the bacteria of the honey bee microbiome and cause evident stress in bees. Microsporidia have been shown to lack an essential binding site on the ribosome that is known to allow for maximal inhibition by paromomycin. Thus, it is possible that paromomycin impacts parasite levels through non-cell autonomous effects on microsporidia infection levels via effects on the microbiome or midgut cellular function. As paromomycin treatment could cause widespread honey bee health issues in agricultural settings, it does not represent an appropriate anti-microsporidia agent for use in the field.

Comparative mitochondrial genomes of the Rhus gall aphid Kaburagia rhusicola subspecies with variable gall shapes

Rhus gall aphids (Hemiptera: Aphididae: Eriosomatinae) stimulate the formation of galls on their primary host plants (sumacs: Rhus spp., Anacardiaceae). The shapes of galls are often used as an extended phenotype to identify the aphid species and subspecies. We collected four Rhus galls with conspicuously different shapes formed by Kaburagia rhusicola aphids, whose sequences of the complete mitochondrial genomes (mitogenomes) were obtained by high-throughput sequencing. Each mitogenome was assembled into a circular molecule containing 13 protein-coding genes, two rRNAs, 22 tRNAs and one control region. All the protein-coding genes had a typical ATN initiation codon and TAA termination codon except for cox1 and nad4, which had a single T as stop codon. All the tRNAs could be folded as a typical clover-leaf secondary structure, except for trnS1 lacking a dihydrouridine (DHU) arm. The relative synonymous codon usage and ratio of nonsynonymous to synonymous substitution rates showed that the four K. rhusicola samples were highly similar to the subspecies K. r. ovogallis. The phylogenetic analyses grouped these samples with K. r. ovogallis in a clade sister to K. r. rhusicola. All these molecular analyses demonstrated that our current samples represented one subspecies of Kaburagia rhusicola, i.e., K. r. ovogallis, and the gall shape was variable even at the subspecies level in Kaburagia gall aphids.

Dive into the sea: first molecular phylogenetic evidence of host expansion from terrestrial/freshwater to marine organisms in Mermithidae (Nematoda: Mermithida)

We report the first mermithid nematode found to be parasitic in a marine tanaidacean crustacean. Ten host tanaidaceans were collected from a depth of 52 m in Otsuchi Bay, Iwate, Japan, north-western Pacific, and identified as a species in the tanaidid genus Zeuxo Templeton, 1840. Nematodes occurred in the host's body cavity; in one case, at least two individuals inhabited a single host. We provide a brief description and illustrations of the morphology of the nematode. In a phylogenetic reconstruction based on the 18S ribosomal RNA gene, the nematode nested in a clade otherwise containing mermithids from terrestrial or freshwater hosts, showing an expansion in host utilization in Mermithidae Braun, 1883 from terrestrial/freshwater hosts to a marine organism.

Comparative analysis of four complete mitogenomes from hoverfly genus Eristalinus with phylogenetic implications

The genus Eristalinus is widely distributed globally. Four complete mitochondrial genomes (i.e., mitogenomes) of Eristalinus were sequenced and analyzed in this study: Eristalinus viridis (Coquillett, 1898), E. quinquestriatus (Fabricius, 1781), E. tarsalis (Macquart, 1855), and E. sp. Within these four sequenced mitogenomes, most protein-coding genes (ND2, CO1, COX2, COX3, ND3, ND5, ND4, ND4L, ND6, and Cytb) began with a typical ATN (T/C/G/A) start codon and ended with a stop codon TAA or incomplete T, whereas ND1 began with the start codon TTG. ND3 ended with TAG. The secondary tRNA structure was that of a typical cloverleaf, and only the tRNA-Ser1 lacked a DHU arm. Three and five domains appeared in the 12S and 16S rRNA secondary structures, respectively. The phylogenetic relationships among the four Eristalinus species combined with the published mitogenomes of Syrphidae were reconstructed using the maximum likelihood and Bayesian inference methods, which support the monophyly of the subfamily Syrphinae but do not support that of the subfamily Eristalinae. Of note, Eristalini and Syrphini are monophyletic groups. The mitogenomes of E. viridis, E. quinquestriatus, E. sp., and E. tarsalis are useful for determining the phylogenetic relationships and evolution of Syrphidae.

Sequence-selective purification of biological RNAs using DNA nanoswitches

Nucleic acid purification is a critical aspect of biomedical research and a multibillion-dollar industry. Here we establish sequence-selective RNA capture, release, and isolation using conformationally responsive DNA nanoswitches. We validate purification of specific RNAs ranging in size from 22 to 401 nt with up to 75% recovery and 99.98% purity in a benchtop process with minimal expense and equipment. Our method compared favorably with bead-based extraction of an endogenous microRNA from cellular total RNA, and can be programmed for multiplexed purification of multiple individual RNA targets from one sample. Coupling our approach with downstream LC/MS, we analyzed RNA modifications in 5.8S ribosomal RNA, and found 2'-O-methylguanosine, 2'-O-methyluridine, and pseudouridine in a ratio of ~1:7:22. The simplicity, low cost, and low sample requirements of our method make it suitable for easy adoption, and the versatility of the approach provides opportunities to expand the strategy to other biomolecules.

The mitochondrial genomes of Tortricidae: nucleotide composition, gene variation and phylogenetic performance

BACKGROUND

Mitochondrial genomes (mitogenomes) have greatly improved our understanding of the backbone phylogeny of Lepidoptera, but few studies on comparative mitogenomics below the family level have been conducted. Here, we generated 13 mitogenomes of eight tortricid species, reannotated 27 previously reported mitogenomes, and systematically performed a comparative analysis of nucleotide composition, gene variation and phylogenetic performance.

RESULTS

The lengths of completely sequenced mitogenomes ranged from 15,440 bp to 15,778 bp, and the gene content and organization were conserved in Tortricidae and typical for Lepidoptera. Analyses of AT-skew and GC-skew, the effective number of codons and the codon bias index all show a base bias in Tortricidae, with little heterogeneity among the major tortricid groups. Variations in the divergence rates among 13 protein-coding genes of the same tortricid subgroup and of the same PCG among tortricid subgroups were detected. The secondary structures of 22 transfer RNA genes and two ribosomal RNA genes were predicted and comparatively illustrated, showing evolutionary heterogeneity among different RNAs or different regions of the same RNA. The phylogenetic uncertainty of Enarmoniini in Tortricidae was confirmed. The synonymy of Bactrini and Olethreutini was confirmed for the first time, with the representative Bactrini consistently nesting in the Olethreutini clade. Nad6 exhibits the highest phylogenetic informativeness from the root to the tip of the resulting tree, and the combination of the third coding positions of 13 protein-coding genes shows extremely high phylogenetic informativeness.

CONCLUSIONS

This study presents 13 mitogenomes of eight tortricid species and represents the first detailed comparative mitogenomics study of Tortricidae. The results further our understanding of the evolutionary architectures of tortricid mitogenomes and provide a basis for future studies of population genetics and phylogenetic investigations in this group.

Mitochondrial Genomes of Hestina persimilis and Hestinalis nama (Lepidoptera, Nymphalidae): Genome Description and Phylogenetic Implications

Simple Summary

In this study, the mitogenomes of Hestina persimilis and Hestinalis nama were obtained via sanger sequencing. Compared with other mitogenomes of Apaturinae butterflies, conclusions can be made that the mitogenomes of Hestina persimilis and Hestinalis nama are highly conservative. The phylogenetic trees build upon mitogenomic data showing that the relationships among Nymphalidae are similar to previous studies. Hestinalisnama is apart from Hestina, and closely related to Apatura, forming a monophyletic clade.

Abstract

In this study, the complete mitochondrial genomes (mitogenomes) of Hestina persimilis and Hestinalis nama (Nymphalidae: Apaturinae) were acquired. The mitogenomes of H. persimilis and H. nama are 15,252 bp and 15,208 bp in length, respectively. These two mitogenomes have the typical composition, including 37 genes and a control region. The start codons of the protein-coding genes (PCGs) in the two mitogenomes are the typical codon pattern ATN, except CGA in the cox1 gene. Twenty-one tRNA genes show a typical clover leaf structure, however, trnS1(AGN) lacks the dihydrouridine (DHU) stem. The secondary structures of rrnL and rrnS of two species were predicted, and there are several new stem loops near the 5′ of rrnL secondary structure. Based on comparative genomic analysis, four similar conservative structures can be found in the control regions of these two mitogenomes. The phylogenetic analyses were performed on mitogenomes of Nymphalidae. The phylogenetic trees show that the relationships among Nymphalidae are generally identical to previous studies, as follows: Libytheinae\Danainae + ((Calinaginae + Satyrinae) + Danainae\Libytheinae + ((Heliconiinae + Limenitidinae) + (Nymphalinae + (Apaturinae + Biblidinae)))). Hestinalisnama is apart from Hestina, and closely related to Apatura, forming monophyly.

New genus of gelechiid moths (Lepidoptera: Gelechiidae) from Borneo, and its four new species

Tawaya gen. n. and its four new species (T. flaventia sp. n., T. rutila sp. n., T. luteola sp. n. and T. armeniaca sp. n.) are described from Borneo. A key to the described species is provided. The taxonomic position of a new genus within the subfamily Dichomeridinae (Gelechiidae) is proposed based on a combination of molecular and morphological analyses.

Characterization and Phylogenetic Analysis of the Mitochondrial Genome Sequence of Nisia fuliginosa (Hemiptera: Fulgoroidea: Meenoplidae)

We explored characterization of the mitochondrial genome (mitogenome or mtGenome) and phylogenetic analysis between 32 Fulgoroid species by sequencing and analyzing the mitogenome of Nisia fuliginosa Yang and Hu, 1985 (Hemiptera: Fulgoroidea: Meenoplidae), thereby making it the first determined mitogenome from the family Meenoplidae. The mitogenome was found to be 15,754 bp in length and contained 13 protein-coding genes (PCGs), 22 tRNA genes, two ribosomal RNA genes (rRNAs), and a control region. All PCGs started with typical ATN codons, except for nad1, which used GTG as the start codon. Canonical TAA termination codons were found in 10 PCGs and the remaining three genes (cox2, nad6, and nad1) had incomplete stop codons T. All tRNAs could fold into typical cloverleaf secondary structures, with the exception of trnC, trnV, and trnS1. Additionally, we compared the AT and GC skews of 13 PCGs of 32 Fulgoroidea mitogenomes, on the L-strand, the AT and GC skews were negative and positive, respectively. However, on the H-strand, the AT skew could be positive or negative and the GC skew was always negative. Phylogenetic results showed that the eight families of Fulgoroidea were divided into two large groups. Delphacidae formed a monophyletic group sister to a clade comprising Meenoplidae and other six families (Fulgoridae, Ricaniidae, Flatidae, Issidae, Caliscelidae, and Achilidae). Meenoplidae was located near the clade of Delphacidae, and Fulgoridae was located near the clade of Meenoplidae. Furthermore, Caliscelidae, Issidae, Ricaniidae, and Flatidae are closely related and they collectively formed a sister group to Achilidae.

18S rRNA variability maps reveal three highly divergent, conserved motifs within Rotifera

BACKGROUND

18S rRNA is a major component of the small subunit of the eukaryotic ribosome and an important phylogenetic marker for many groups, often to the point of being the only marker available for some. A core structure across eukaryotes exists for this molecule that can help to inform about its evolution in different groups. Using an alignment of 18S rDNA for Rotifera as traditionally recognized (=Bdelloidea, Monogononta, and Seisonacea, but not Acanthocephala), I fitted sequences for three exemplar species (Adineta vaga, Brachionus plicatilis, and Seison nebaliae, respectively) to the core structure and used these maps to reveal patterns of evolution for the remainder of this diverse group of microscopic animals.

RESULTS

The obtained variability maps of the 18S rRNA molecule revealed a pattern of high diversity among the three major rotifer clades coupled with strong conservation within each of bdelloids and monogononts. A majority of individual sites (ca. 60%) were constant even across rotifers as a whole with variable sites showing only intermediate rates of evolution. Although the three structural maps each showed good agreement with the inferred core structure for eukaryotic 18S rRNA and so were highly similar to one another at the secondary and tertiary levels, the overall pattern is of three highly distinct, but conserved motifs within the group at the primary sequence level. A novel finding was that of a variably expressed deletion at the 3' end of the V3 hypervariable region among some bdelloid species that occasionally extended into and included the pseudoknot structure following this region as well as the central "square" of the 18S rRNA molecule. Compared to other groups, levels of variation and rates of evolution for 18S rRNA in Rotifera roughly matched those for Gastropoda and Acanthocephala, despite increasing evidence for the latter being a clade within Rotifera.

CONCLUSIONS

The lack of comparative data for comparable groups makes interpretation of the results (i.e., very low variation within each of the three major rotifer clades, but high variation between them) and their potential novelty difficult. However, these findings in combination with the high morphological diversity within rotifers potentially help to explain why no clear consensus has been reached to date with regard to the phylogenetic relationships among the major groups.

Mitochondrial Phylogenomics of Tenthredinidae (Hymenoptera: Tenthredinoidea) Supports the Monophyly of Megabelesesinae as a Subfamily

Tenthredinidae represents one of the external feeders of the most diverse superfamily, Tenthredinoidea, with diverse host plant utilization. In this study, four complete mitochondrial genomes (mitogenomes), those of Cladiucha punctata, Cladiucha magnoliae, Megabeleses magnoliae, and Megabeleses liriodendrovorax, are newly sequenced and comparatively analyzed with previously reported tenthredinid mitogenomes. The close investigation of mitogenomes and the phylogeny of Tenthredinidae leads us to the following conclusions: The subfamilial relationships and phylogenetic placements within Tenthredinidae are mostly found to be similar to the previously suggested phylogenies. However, the present phylogeny supports the monophyly of Megabelesesinae as a subfamily, with the sister-group placement of Cladiucha and Megabeleses outside of Allantinae. The occurrence of the same type of tRNA rearrangements (MQI and ANS1ERF) in the mitogenomes of Megabelesesinae species and the presence of apomorphic morphological characters also provide robust evidence for this new subfamily. The divergence and diversification times of the subfamilies appear to be directly related to colonization of the flowering plants following the Early Cretaceous. The origin time and diversification patterns of Megabelesesinae were also well matched with the divergence times of their host plants from Magnoliaceae.

Comparative Mitogenomic Analysis of Two Cuckoo Bees (Apoidea: Anthophila: Megachilidae) with Phylogenetic Implications

Bees (Hymenoptera, Apoidea and Anthophila) are distributed worldwide and considered the primary pollinators of angiosperm. Megachilidae is one of the largest families of Anthophila. In this study, two complete mitogenomes of cuckoo bees in Megachilidae, namely Coelioxys fenestrata and Euaspis polynesia, were amplified and sequenced, with a length of 17,004 bp (C. fenestrata) and 17,682 bp (E. polynesia). The obtained results show that 37 mitogenomic genes and one putative control region were conserved within Hymenoptera. Truncated stop codon T was found in the cox3 gene of E. polynesia. The secondary structure of small (rrnS) and large (rrnL) rRNA subunits contained three domains (28 helices) and five domains (44 helices) conserved within Hymenoptera, respectively. Compared with ancestral gene order, gene rearrangement events included local inversion and gene shuffling. In order to reveal the phylogenetic position of cuckoo bees, we performed phylogenetic analysis. The results supported that all families of Anthophila were monophyletic, the tribe-level relationship of Megachilidae was Osmiini + (Anthidiini + Megachilini) and Coelioxys fenestrata was clustered to the Megachile genus, which was more closely related to Megachile sculpturalis and Megachile strupigera than Euaspis polynesia.

Molecular phylogeny of the phytoparasitic mite family Phytoptidae (Acariformes: Eriophyoidea) identified the female genitalic anatomy as a major macroevolutionary factor and revealed multiple origins of gall induction

Phytoptidae s.str. is a lineage of eriophyoid mites associated with angiosperms. Based on representative taxon sampling and four gene markers (COI, HSP70, 18S, and 28S), we inferred the molecular phylogeny of this group and performed comparative analyses of cuticle-lined female internal genitalia. Although basal relationships were unclear, several well supported clades were recovered. These clades were supported by geography, host associations, and female genital anatomy, but contradicted the current morphology-based systematics. The monophyly of each of five conventional supraspecific groupings (Fragariocoptes, Phytoptus, Phytoptinae, Sierraphytoptinae, and Sierraphytoptini) is rejected based on a series of statistical tests. Additionally, four morphological characters (the absence of tibial solenidion φ and opisthosomal seta c1, presence of telosomal pseudotagma, and 'morphotype') were found to be homoplasies that cannot be used to confidently delimit supraspecific lineages of phytoptids. However, our molecular topology was highly congruent with female genital characters. Eight molecular clades were unambiguously supported by the shapes and topography of the spermathecal apparatus and genital apodemes. This suggests that the female genital anatomy could be an important factor affecting cladogenesis in Phytoptidae, a conclusion contrasting with the general expectation that host characteristics should be a major macroevolutionary force influencing the evolution of host-specific symbionts. Indeed, despite the high host-specificity, there were no apparent cophylogenetic patterns. Furthermore, we show that gall-inducing ability evolved multiple times in phytoptids. Because gall formation creates nearly instantaneous niche partitioning and the potential loss or reduction of gene flow, we hypothesize that it could be an important evolutionary factor affecting speciation within different host-associated clades of phytoptid mites.

Evaluation of the expression stability of reference genes in Apis mellifera under pyrethroid treatment

Honeybees (Apis mellifera L.), which unquestionably play an economically important role in pollination and agricultural production, are at risk of decline. To study changes in gene expression in insects upon exposure to pesticides or other external stimuli, appropriate reference genes are required for data normalization. Since there is no such gene that is absolutely invariable under all experimental conditions, the aim of this study was to identify the most stable targets suitable for subsequent normalization in quantitative experiments based on real-time polymerase chain reaction in honeybee research. Here, we evaluated the expression of fifteen candidate housekeeping genes from three breeding lines of honeybees treated with pyrethroids to identify the most stable genes. The tested insects were exposed to deltamethrin or lambda-cyhalothrin, and then, changes in the accumulation of selected transcripts were assessed, followed by statistical analyses. We concluded that AmRPL32, AmACT and AmRPL13a were the commonly recorded most stable genes in honeybees treated with the selected pyrethroids.

Comparative mitogenomes of six species in the subfamily Iassinae (Hemiptera: Cicadellidae) and phylogenetic analysis

For elucidating phylogenetic relationships among members of the family Cicadellidae, mitogenomes of six members of the subfamily Iassinae, including Batracomorphus lateprocessus, Iassus dorsalis, Krisna concava, Krisna rufimarginata, Gessius rufidorsus, and Trocnadella arisana, were sequenced. The mitogenomes are 14,724-15,356-bp long. Moreover, typical 37 genes in mitogenomes were identified; arrangement of these genes in the studied species was consistent with that in the inferred ancestral insects, except for tRNA genes, with a simple switch between the positions of trnI and trnQ in Trocnadella arisana. Most protein-coding genes in the Iassinae mitogenomes showed typical ATN start codons (ATA/ATT/ATC/ATG) and TAR (TAA/TAG) or an incomplete stop codon T--; ATP8 of all sequenced species showed the start codon TTG. The secondary structures of 16S rRNA and 12S rRNA were predicted. 16S rRNA comprised 6 domains with 43 helices, and 12S rRNA comprised 3 domains with 25 helices. All subfamilies, except Cicadellinae and Evacanthinae, were recovered as monophyletic. As reported previously, treehoppers originated from paraphyletic Cicadellidae. Iassinae and Coelidiinae, Megophthalminae and treehoppers, and Cicadellinae and Evacanthinae were sister groups with high nodal support.

DNA Extraction of Sumatran Striped Rabbit from Tissue Samples

The Sumatran Striped Rabbit (Nesolagus netscheri ) is  a naturally rare, which is a protected animal based on Government Regulation Number 7 of 1999 and renewal of Permen LHK No. P.106 which are categorized as Data Deficient  by the IUCN (International Union of Conservation of Nature) since 2019. Samples were obtained from Pagar Alam, South Sumatra. Most commonly found at elevations  above 600 m in montane and sub-montane primary habitat. Molecular genetic characterization of sumatran striped rabbit is to know genetic information and genetic identification of N. netscheri. Primers 12S rRNA were used in this study to characterize Sumatran striped rabbit. A 1 band were detected ranged from 900 to 1000 bp. Molecular markers represent reliable tools which may have a great impact in rabbit breeding and genetic improvement of rabbits. Molecular markers on Nesolagus netscheri is expected to tools in the identification of Sumatran striped rabbits from South Sumatra.

The ribosomal intergenic spacer (IGS) in the potato and tobacco cyst nematodes, Globodera pallida, G. rostochiensis and G. tabacum

The potato cyst nematodes Globodera pallida and G. rostochiensis (PCN), and tobacco cyst nematode (TCN), G. tabacum, are the most important parasitic nematodes of potato and tobacco worldwide. Ribosomal DNA provides useful molecular data for diagnostics, the study of polymorphisms and for evolutionary research in eukaryotic organisms including nematodes. Here we present data on the structure and organization of a rarely studied part of the intergenic spacer (IGS) region of the PCN and TCN genome of cyst nematodes. This region has shown potential for diagnostic purposes and population studies in other organisms including nematodes. In nematodes, the ribosomal RNA gene cluster comprises three genes: 5.8S, 18S and 28S rRNA, which are separated by spacer regions: the intergenic spacer (IGS), non-transcribed spacer (NTS), externally transcribed spacer (EST) and the internally transcribed spacer (ITS). The intergenic spacer (IGS) region consists of an external transcribed spacer (ETS) and a non-transcribed spacer (NTS) which is located between the 28S of one repeat and the 18S gene of the next repeat within the rRNA genes cluster. In this study, the first flanking portion of the IGS was amplified, cloned and sequenced from PCN and TCN. Primers were then designed to amplify the whole IGS sequence. PCR amplification of IGS from G. tabacum, G. pallida, and G. rostochiensis yielded respectively: a single amplicon of 3 kb, three amplicons sized 2.5, 2.6 and 2.9 kb, and two amplicons sized 2.8 and 2.9 kb. Results showed that Globodera spp. has more than one variant copy of the IGS, with both long and short repetitive DNA elements. An approximately 400 bp long region without any internal repetitive elements, were identified in a position between the two repetitive regions suggesting that there is a 5S gene in the IGS of these species.

The SSU rRNA secondary structures of the Plagiorchiida species (Digenea), its applications in systematics and evolutionary inferences

The small subunit ribosomal RNA (SSU rRNA) is widely used phylogenetic marker in broad groups of organisms and its secondary structure increasingly attracts the attention of researchers as supplementary tool in sequence alignment and advanced phylogenetic studies. Its comparative analysis provides a great contribution to evolutionary biology, allowing find out how the SSU rRNA secondary structure originated, developed and evolved. Herein, we provide the first data on the putative SSU rRNA secondary structures of the Plagiorchiida species. The structures were found to be quite conserved across broad range of species studied, well compatible with those of others eukaryotic SSU rRNA and possessed some peculiarities: cross-shaped structure of the ES6b, additional shortened ES6c2 helix, and elongated ES6a helix and h39 + ES9 region. The secondary structures of variable regions ES3 and ES7 appeared to be tissue-specific while ES6 and ES9 were specific at a family level allowing considering them as promising markers for digenean systematics. Their uniqueness more depends on the length than on the nucleotide diversity of primary sequences which evolutionary rates well differ. The findings have important implications for understanding rRNA evolution, developing molecular taxonomy and systematics of Plagiorchiida as well as for constructing new anthelmintic drugs.

Complete mitochondrial genome sequence of Labriocimbex sinicus, a new genus and new species of Cimbicidae (Hymenoptera) from China

Labriocimbex sinicus Yan & Wei gen. et sp. nov. of Cimbicidae is described. The new genus is similar to Praia Andre and Trichiosoma Leach. A key to extant Holarctic genera of Cimbicinae is provided. To identify the phylogenetic placement of Cimbicidae, the mitochondrial genome of L. sinicus was annotated and characterized using high-throughput sequencing data. The complete mitochondrial genome of L. sinicus was obtained with a length of 15,405 bp (GenBank: MH136623; SRA: SRR8270383) and a typical set of 37 genes (22 tRNAs, 13 PCGs, and two rRNAs). The results demonstrated that all PCGs were initiated by ATN codon, and ended with TAA or T stop codons. The study reveals that all tRNA genes have a typical clover-leaf secondary structure, except for trnS1. Remarkably, the secondary structures of the rrnS and rrnL of L. sinicus were much different from those of Corynis lateralis. Phylogenetic analyses verified the monophyly and positions of the three Cimbicidae species within the superfamily Tenthredinoidea and demonstrated a relationship as (Tenthredinidae + Cimbicidae) + (Argidae + Pergidae) with strong nodal supports. Furthermore, we found that the generic relationships of Cimbicidae revealed by the phylogenetic analyses based on COI genes agree quite closely with the systematic arrangement of the genera based on the morphological characters. Phylogenetic tree based on two methods shows that L. sinicus is the sister group of Praia with high support values. We suggest that Labriocimbex belongs to the tribe Trichiosomini of Cimbicinae based on adult morphology and molecular data. Besides, we suggest to promote the subgenus Asitrichiosoma to be a valid genus.

Prospects and challenges of implementing DNA metabarcoding for high-throughput insect surveillance

Trap-based surveillance strategies are widely used for monitoring of invasive insect species, aiming to detect newly arrived exotic taxa as well as track the population levels of established or endemic pests. Where these surveillance traps have low specificity and capture non-target endemic species in excess of the target pests, the need for extensive specimen sorting and identification creates a major diagnostic bottleneck. While the recent development of standardized molecular diagnostics has partly alleviated this requirement, the single specimen per reaction nature of these methods does not readily scale to the sheer number of insects trapped in surveillance programmes. Consequently, target lists are often restricted to a few high-priority pests, allowing unanticipated species to avoid detection and potentially establish populations. DNA metabarcoding has recently emerged as a method for conducting simultaneous, multi-species identification of complex mixed communities and may lend itself ideally to rapid diagnostics of bulk insect trap samples. Moreover, the high-throughput nature of recent sequencing platforms could enable the multiplexing of hundreds of diverse trap samples on a single flow cell, thereby providing the means to dramatically scale up insect surveillance in terms of both the quantity of traps that can be processed concurrently and number of pest species that can be targeted. In this review of the metabarcoding literature, we explore how DNA metabarcoding could be tailored to the detection of invasive insects in a surveillance context and highlight the unique technical and regulatory challenges that must be considered when implementing high-throughput sequencing technologies into sensitive diagnostic applications.

The mitochondrial genome of the springtail Bourletiella arvalis (Symphypleona, Collembola)

The complete mitochondrial genome of the springtail Bourletiella arvalis (Fitch, 1863) is herein described and applied to a Bayesian phylogenetic analysis, inclusive of all the Collembola mitochondrial DNAs sequenced so far. The gene content and order, as well as the nucleotide composition, conform with the well-known features of hexapods’ mitochondrial genomes. The phylogenetic analysis supports the monophyly of Collembola, Poduromorpha, Entomobryomorpha and Symphypleona. However, no mtDNA from Neelipleona is available to date, therefore limiting the application of mitochondrial genomes to further investigate springtail systematics.

Characterization of the mitochondrial genome of Analcellicampa xanthosoma gen. et sp. nov. (Hymenoptera: Tenthredinidae)

A new genus with a new species of the tribe Hoplocampini of Hoplocampinae was described from China: Analcellicampa xanthosoma Wei & Niu, gen. et sp. nov. Hoplocampa danfengensis G. Xiao 1994 was designated as the type species of the new genus. The characters of Analcellicampa danfengensis (G. Xiao) comb. nov. were briefly discussed. A key to the tribes and known genera of Hoplocampinae was provided. The nearly complete mitochondrial genome of A. xanthosoma was characterized as having a length of 15,512 bp and containing 37 genes (22 tRNAs, 13 protein-coding genes (PCGs), and 2 rRNAs). The gene order of this new specimen was the same as that in the inferred insect ancestral mitochondrial genome. All PCGs were initiated by ATN codons and ended with TAA or T stop codons. All tRNAs had a typical cloverleaf secondary structure, except for trnS1. Remarkably, the helices H991 of rrnS and H47 of rrnL were redundant, while helix H563 of rrnL was highly conserved. A phylogeny based on previously reported symphytan mitochondrial genomes showed that A. xanthosoma is a sister group to Monocellicampa pruni, with high support values. We suggest that A. xanthosoma and M. pruni belong to the tribe Hoplocampini of Hoplocampinae.

The first two mitochondrial genomes of wood wasps (Hymenoptera: Symphyta): Novel gene rearrangements and higher-level phylogeny of the basal hymenopterans

The Symphyta has long been recognized as a paraphyletic grade forming the base of the remaining Hymenopteran, and the superfamily relationships within Symphyta remain controversial. Here, the first two representative mitochondrial genomes from the superfamily Siricoidea and Xiphydrioidea (Hymenoptera: Symphyta) are obtained using next-generation sequencing. The complete mitochondrial genome of Xiphydria sp. is 16,482 bp long with an A + T content of 84.18% while the incomplete one of Tremex columba is 16,847 bp long and A + T content is 81.69%. All 37 typical mitochondrial genes are possessed in both species. The secondary structure of tRNAs and rRNAs for both species are successfully predicted. Compared with the ancestral organization, seven and five tRNA genes are rearranged in mitochondrial genomes of Tremex and Xiphydria, respectively, which are the most rearrangement events within Symphyta. The rearrangement patterns in Tremex and Xiphydria present in this study are all novel to the Symphyta. Phylogenetic relationships among the major lineages of Symphyta are reconstructed using mitochondrial genomes. Both maximum likelihood and Bayesian inference analyses highly support Symphyta is a paraphyletic grade, Xyeloidea + (Tenthredinoidea + (Pamphilioidea + (Xiphydrioidea + (Cephoidea + (Orussoidea + Apocrita))))).

Comparative analysis of the mitochondrial genomes of oriental spittlebug trible Cosmoscartini: insights into the relationships among closely related taxa

Background

Cosmoscartini (Hemiptera: Cercopoidea: Cercopidae) is a large and brightly colored Old World tropical tribe, currently containing over 310 phytophagous species (including some economically important pests of eucalyptus in China) in approximately 17 genera. However, very limited information of Cosmoscartini is available except for some scattered taxonomic studies. Even less is known about its phylogenetic relationship, especially among closely related genera or species. In this study, the detailed comparative genomic and phylogenetic analyses were performed on nine newly sequenced mitochondrial genomes (mitogenomes) of Cosmoscartini, with the purpose of exploring the taxonomic status of the previously defined genus Okiscarta and some closely related species within the genus Cosmoscarta.

Results

Mitogenomes of Cosmoscartini display similar genomic characters in terms of gene arrangement, nucleotide composition, codon usage and overlapping regions. However, there are also many differences in intergenic spacers, mismatches of tRNAs, and the control region. Additionally, the secondary structures of rRNAs within Cercopidae are inferred for the first time.

Based on comparative genomic (especially for the substitution pattern of tRNA secondary structure) and phylogenetic analyses, the representative species of Okiscarta uchidae possesses similar structures with other Cosmoscarta species and is placed consistently in Cosmoscarta. Although Cosmoscarta bimacula is difficult to be distinguished from Cosmoscarta bispecularis by traditional morphological methods, evidence from mitogenomes highly support the relationships of (C. bimacula + Cosmoscarta rubroscutellata) + (C. bispecularis + Cosmoscarta sp.).

Conclusions

This study presents mitogenomes of nine Cosmoscartini species and represents the first detailed comparative genomic and phylogenetic analyses within Cercopidae. It is indicated that knowledge of mitogenomes can be effectively used to resolve phylogenetic relationships at low taxonomic levels. Sequencing more mitogenomes at various taxonomic levels will also improve our understanding of mitogenomic evolution and phylogeny in Cercopidae.

Electronic supplementary material

The online version of this article (10.1186/s12864-018-5365-7) contains supplementary material, which is available to authorized users.

Characterization of the mitochondrial genome of Arge bella Wei & Du sp. nov. (Hymenoptera: Argidae)

We describe Arge bella Wei & Du sp. nov., a large and beautiful species of Argidae from south China, and report its mitochondrial genome based on high-throughput sequencing data. We present the gene order, nucleotide composition of protein-coding genes (PCGs), and the secondary structures of RNA genes. The nearly complete mitochondrial genome of A. bella has a length of 15,576 bp and a typical set of 37 genes (22 tRNAs, 13 PCGs, and 2 rRNAs). Three tRNAs are rearranged in the A. bella mitochondrial genome as compared to the ancestral type in insects: trnM and trnQ are shuffled, while trnW is translocated from the trnW-trnC-trnY cluster to a location downstream of trnI. All PCGs are initiated by ATN codons, and terminated with TAA, TA or T as stop codons. All tRNAs have a typical cloverleaf secondary structure, except for trnS1. H821 of rrnS and H976 of rrnL are redundant. A phylogenetic analysis based on mitochondrial genome sequences of A. bella, 21 other symphytan species, two apocritan representatives, and four outgroup taxa supports the placement of Argidae as sister to the Pergidae within the symphytan superfamily Tenthredinoidea.

Assessing the Efficiency of Molecular Markers for the Species Identification of Gregarines Isolated from the Mealworm and Super Worm Midgut

Protozoa, of the taxon Gregarinasina, are a heterogeneous group of Apicomplexa that includes ~1600 species. They are parasites of a large variety of both marine and terrestrial invertebrates, mainly annelids, arthropods and mollusks. Unlike coccidians and heamosporidians, gregarines have not proven to have a negative effect on human welfare; thus, they have been poorly investigated. This study focuses on the molecular identification and phylogeny of the gregarine species found in the midgut of two insect species that are considered as an alternative source of animal proteins for the human diet: the mealworm Tenebrio molitor, and the super-worm Zophobas atratus (Coleoptera: Tenebrionidae). Gregarine specimens were isolated from the gut of both larval and adult stages of T. molitor specimens, as well as from Z. atratus larvae. The morphological analyses were restricted to the identification of the different parasite morpho-types, likely corresponding either to different life-cycle stages or to alternative gregarine species. The samples were also used for the DNA extraction necessary for their genetic characterization. Finally, the efficiency of different molecular markers (i.e., 18S rDNA gene alone or combined with the Internal Transcribed Spacer 1) was assessed when applied either to gregarine species identification and to phylogenetic inference.

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 Mitochondrial Genomes of Phytophagous Scarab Beetles and Systematic Implications

In this study, we newly sequenced five mitogenomes of representatives of phytophagous scarab beetles (Coleoptera: Scarabaeidae) by using next-generation sequencing technology. Two species have complete (or nearly complete) mitogenome sequences, namely Popillia mutans Newman (Coleoptera: Scarabaeidae) and Holotrichia oblita Faldermann (Coleoptera: Scarabaeidae). The remaining three species have the partial mitogenomes, and the missing genes are mainly located adjacent to the control region. The complete (or nearly complete) mitogenomes have the same genome structure as most of the existing Scarabaeidae mitogenomes. We conducted phylogenetic analyses together with 24 published mitogenomes of Scarabaeoidea. The results supported a basal split of coprophagous and phytophagous Scarabaeidae. The subfamily Sericinae was recovered as sister to all other phytophagous scarab beetles. All analyses supported a non-monophyletic Melolonthinae, which included two different non-sister clades. The Cetoniinae was recovered as sister to a clade including Rutelinae and Dynastinae. Although the Rutelinae was rendered paraphyletic by Dynastinae in the Bayesian trees inferred under the site-heterogeneous CAT-GTR or CAT-MTART model, discordant patterns were given in some of ML trees estimated using the homogeneous GTR model.

A Combined Molecular and Morphological Approach to Explore the Higher Phylogeny of Entimine Weevils (Coleoptera: Curculionidae), with Special Reference to South American Taxa

The Entiminae are broad-nosed weevils constituting the most diverse subfamily of Curculionidae, with over 50 tribes. We performed Bayesian and Maximum Parsimony combined phylogenetic analyses with the main objective of testing higher-level relationships and the naturalness of the major Neotropical and Southern South American (Patagonia and Andes) tribes, including some members from other regions. We compiled a data matrix of 67 terminal units with 63 Entiminae species, as well as four outgroup taxa from Cyclominae, by 3522 molecular (from nuclear 18S rDNA and 28S rDNA, and mitochondrial 16S rDNA and COI gene sequences) and 70 morphological characters. The resulting trees recover a clade Entiminae with a monophyletic Cylydrorhinini and Premnotrypes branching off early. The tree resulting from parsimony analysis shows a clade of Leptopiini from the Australian region and another clade including taxa mainly distributed in the Palaearctic and Neotropical regions, but in the Bayesian tree the South American and Australian Leptopiini are grouped together. The mainly Palaearctic Entiminae (e.g., Brachyderini, Laparocerini, Otiorhynchini, Peritelini, Polydrusini, Phyllobiini and Sciaphylini) form a subclade separated from Southern Hemisphere taxa. Among the latter, the well-supported Naupactini are the sister group of the South American Tanymecini, excluding Platyaspistes, herein transferred to Leptopiini (new placement). Another well-justified clade is Eustylini–Geonemini, which also includes the enigmatic Galapagonotus, and the genus Artipus, thus corroborating its recent exclusion from Naupactini.

Morphological and Molecular Perspectives on the Phylogeny, Evolution, and Classification of Weevils (Coleoptera: Curculionoidea): Proceedings from the 2016 International Weevil Meeting

The 2016 International Weevil Meeting was held immediately after the International Congress of Entomology (ICE). It built on the topics and content of the 2016 ICE weevil symposium Phylogeny and Evolution of Weevils (Coleoptera: Curculionoidea): A Symposium in Honor of Dr. Guillermo "Willy” Kuschel. Beyond catalyzing research and collaboration, the meeting was intended to serve as a forum for identifying priorities and goals for those who study weevils. The meeting consisted of 46 invited and contributed lectures, discussion sessions and introductory remarks presented by 23 speakers along with eight contributed research posters. These were organized into three convened sessions, each lasting one day: (1) weevil morphology; (2) weevil fossils, biogeography and host/habitat associations; and (3) molecular phylogenetics and classification of weevils. Some of the topics covered included the 1K Weevils Project, major morphological character systems of adult and larval weevils, weevil morphological terminology, prospects for future morphological character discovery, phylogenetic analysis of morphological character data, the current status of weevil molecular phylogenetics and evolution, resources available for phylogenetic and comparative genomic studies of weevils, the weevil fossil record, weevil biogeography and evolution, weevil host plants, evolutionary development of the weevil rostrum, resources available for weevil identification and the current status of and challenges in weevil classification.

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.

Cryptic species in a well-known habitat: applying taxonomics to the amphipod genus Epimeria (Crustacea, Peracarida)

Taxonomy plays a central role in biological sciences. It provides a communication system for scientists as it aims to enable correct identification of the studied organisms. As a consequence, species descriptions should seek to include as much available information as possible at species level to follow an integrative concept of 'taxonomics'. Here, we describe the cryptic species Epimeria frankei sp. nov. from the North Sea, and also redescribe its sister species, Epimeria cornigera. The morphological information obtained is substantiated by DNA barcodes and complete nuclear 18S rRNA gene sequences. In addition, we provide, for the first time, full mitochondrial genome data as part of a metazoan species description for a holotype, as well as the neotype. This study represents the first successful implementation of the recently proposed concept of taxonomics, using data from high-throughput technologies for integrative taxonomic studies, allowing the highest level of confidence for both biodiversity and ecological research.

The Complete Mitochondrial Genome of Ugyops sp. (Hemiptera: Delphacidae)

The complete mitochondrial genome (mitogenome) of Ugyops sp. (Hemiptera: Delphacidae) was sequenced, making it the first determined mitogenome from the subfamily Asiracinae, the basal clade of the family Delphacidae. The mitogenome was 15,259 bp in length with A + T content of 77.65% and contained 13 protein-coding genes (PCGs), 22 transfer RNA genes (tRNAs), two ribosomal RNA genes (rRNAs), and a control region. The gene order was identical with that of the ancestral insect. The nucleotide composition analysis indicated that the whole mitogenome was strongly A-skewed (0.288) and highly C-skewed (-0.270). For PCGs on the J-strand, the AT skew was positive, and the GC skew was negative. All PCGs started with canonical ATN codons, except for cox1 and nad5, which used CTG and GTG as start codon, respectively. All tRNAs could fold into typical cloverleaf secondary structures, with the exception of trnS1 (AGN), in which the dihydrouridine arm was reduced to a simple loop. The control region included a poly-T stretch downstream of the small rRNA gene (rrnS), a subregion of higher A + T content and tandemly repeated sequence near trnI. The mitogenome of Ugyops sp. could be very helpful in exploring the diversity and evolution of mitogenomes in Delphacidae.