Characterization of the complete mitochondrial genomes of Maiestas dorsalis and Japananus hyalinus (Hemiptera: Cicadellidae) and comparison with other Membracoidea

Mitogenomic Analysis and Phylogenetic Implications for the Deltocephaline Tribe Chiasmini (Hemiptera: Cicadellidae: Deltocephalinae)

The grassland leafhopper tribe Chiasmini (Cicadellidae: Deltocephalinae) presently comprises 324 described species worldwide, with the highest species diversity occurring in the Nearctic region but a greater diversity of genera occurring in the Old World. In China, this tribe comprises 39 described species in 11 genera, but the fauna remains understudied. The complete mitogenomes of three species of this tribe have been sequenced previously. In order to better understand the phylogenetic position of Chiasmini within the subfamily Deltocephalinae and to investigate relationships among Chiasmini genera and species, we sequenced and analyzed the complete mitogenomes of 13 species belonging to seven genera from China. Comparison of the newly sequenced mitogenomes reveals a closed circular double-stranded structure containing 37 genes with a total length of 14,805 to 16,269 bp and a variable number of non-coding A + T-rich regions. The gene size, gene order, gene arrangement, base composition, codon usage, and secondary structure of tRNAs of the newly sequenced mitogenomes of these 13 species are highly conserved in Chiasmini. The ATN codon is commonly used as the start codon in protein-coding genes (PCGs), except for ND5 in Doratura sp. and ATP6 in Nephotettix nigropictus, which use the rare GTG start codon. Most protein-coding genes have TAA or TAG as the stop codon, but some genes have an incomplete T stop codon. Except for the tRNA for serine (trnS1(AGN)), the secondary structure of the other 21 tRNAs is a typical cloverleaf structure. In addition to the primary type of G-U mismatch, five other types of tRNA mismatches were observed: A-A, A-C, A-G, U-C, and U-U. Chiasmini mitochondrial genomes exhibit gene overlaps with three relatively stable regions: the overlapping sequence between trnW and trnC is AAGTCTTA, the overlapping sequence between ATP8 and ATP6 is generally ATGATTA, and the overlapping sequence between ND4 and ND4L is generally TTATCAT. The largest non-coding region is the control region, which exhibits significant length and compositional variation among species. Some Chiasmini have tandem repeat structures within their control regions. Unlike some other deltocephaline leafhoppers, the sequenced Chiasmini lack mitochondrial gene rearrangements. Phylogenetic analyses of different combinations of protein-coding and ribosomal genes using maximum likelihood and Bayesian methods under different models, using either amino acid or nucleotide sequences, are generally consistent and also agree with results of prior analyses of nuclear and partial mitochondrial gene sequence data, indicating that complete mitochondrial genomes are phylogenetically informative at different levels of divergence within Chiasmini and among leafhoppers in general. Apart from Athysanini and Opsiini, most of the deltocephaline tribes are recovered as monophyletic. The results of ML and BI analyses show that Chiasmini is a monophyletic group with seven monophyletic genera arranged as follows: ((Zahniserius + (Gurawa + (Doratura + Aconurella))) + (Leofa + (Exitianus + Nephotettix))).

First complete mitogenome of Massarineae and its contribution to phylogenetic implications in Pleosporales

Endophytic fungi play an important role in the growth and development of traditional Chinese medicine plants. We isolated a strain of Acrocalymma vagum from the endophytic fungi of the traditional Chinese plants Paris. To accurately identify this endophytic fungal species of interest, we sequenced the mitochondrial genome of A. vagum, which is the first discovered mitochondrial genome in Massarineae. The A. vagum mitochondrial genome consists of a 35,079-bp closed circular DNA molecule containing 36 genes. Then, we compared the general sequence characteristics of A. vagum with those of Pleosporales, and the second structure of the 22 tRNAs was predicted. The phylogenetic relationship of A. vagum was constructed using two different data sets (protein-coding genes and amino acids). The phylogenetic tree shows that A. vagum is located at the root of Pleosporales. The analysis of introns shows that the number of introns increases with the increase in branch length. The results showed that monophyly was confirmed for all families in Pleosporales except for Pleosporaceae. A. vagum is an ancient species in the Pleosporales, and Pleosporaceae may require further revision. In Pleosporales, the number of introns is positively correlated with branch length, providing data for further study on the origin of introns.

Comparative mitochondrial genome analysis of three leafhopper species of the genus Abrus Dai & Zhang (Hemiptera: Cicadellidae: Deltocephalinae) from China with phylogenetic implication

BACKGROUND

The phylogenetic position and classification of Athysanini are poorly defined, as it includes a large group of polyphyletic genera that have historically been assigned to it mainly because they still exhibit the most typical deltocephaline genitalic and external body characters but lack the distinctive characteristics that other tribes possess. The bamboo-feeding leafhopper genus Abrus belong to the tribe Athysanini of subfamily Deltocephalinae, which currently comprises 19 valid described species, and are limited to the Oriental and Palaearctic regions in China. Although the taxonomy of Abrus are well updated, the references on comparative mitogenomic analyses of Abrus species are only known for a single species. In this study, we sequenced and analyzed the complete mitochondrial genomes (mitogenomes) of Abrus daozhenensis Chen, Yang & Li, 2012 (16,391bp) and A. yunshanensis Chen, Yang & Li, 2012 (15,768bp) (Athysanini), and compared with published mitogenome sequence of A. expansivus Xing & Li, 2014 (15,904bp).

RESULTS

These Abrus species shared highly conserved mitogenomes with similar gene order to that of the putative ancestral insect with 37 typical genes and a non-coding A + T-rich region. The nucleotide composition of these genomes is highly biased toward A + T nucleotides (76.2%, 76.3%, and 74.7%), AT-skews (0.091 to 0.095, and 0.095), negative GC-skews (- 0.138, - 0.161, and - 0.138), and codon usage. All 22 tRNA genes had typical cloverleaf secondary structures, except for trnS1 (AGN) which lacks the dihydrouridine arm, and distinctively trnG in the mitogenome of A. expansivus lacks the TψC arm. Phylogenetic analyses based on 13 PCGs, 2 rRNA genes, and 22 tRNA genes consistently recovered the monophyletic Opsiini, Penthimiini, Selenocephalini, Scaphoideini, and Athysanini (except Watanabella graminea, previously sequenced species as Chlorotettix nigromaculatus) based on limited available mitogenome sequence data of 37 species.

CONCLUSION

At present, Abrus belongs to the tribe Athysanini based on both morphological and molecular datasets, which is strongly supported in present phylogenetic analyses in both BI and ML methods using the six concatenated datasets: amino acid sequences and nucleotides from different combinations of protein-coding genes (PCGs), ribosomal RNA (rRNAs), and transfer RNA (tRNAs). Phylogenetic trees reconstructed herein based on the BI and ML analyses consistently recovered monophylitic Athysanini, except Watanabella graminea (Athysanini) in Opsiini with high support values.

Complete Mitochondrial Genome Assembly of an Upland Wild Rice Species, Oryza granulata and Comparative Mitochondrial Genomic Analyses of the Genus Oryza

Wild upland rice species, including Oryza granulata, possess unique characteristics that distinguish them from other Oryza species. For instance, O. granulata characteristically has a GG genome and is accordingly classified as a basal lineage of the genus Oryza. Here, we deployed a versatile hybrid approach by integrating Illumina and PacBio sequencing data to generate a high-quality mitochondrial genome (mitogenome) assembly for O. granulata. The mitogenome of O. granulata was 509,311 base pairs (bp) with sixty-seven genes comprising two circular chromosomes, five ribosomal RNA (rRNA) coding genes, twenty-five transfer RNA (tRNA) coding genes, and thirty-seven genes coding for proteins. We identified a total of 378 simple sequence repeats (SSRs). The genome also contained 643 pairs of dispersed repeats comprising 340 palindromic and 303 forward. In the O. granulata mitogenome, the length of 57 homologous fragments in the chloroplast genome occupied 5.96% of the mitogenome length. Collinearity analysis of three Oryza mitogenomes revealed high structural variability and frequent rearrangements. Phylogenetic analysis showed that, compared to other related genera, O. granulata had the closest genetic relationship with mitogenomes reported for all members of Oryza, and occupies a position at the base of the Oryza phylogeny. Comparative analysis of complete mitochondrial genome assemblies for Oryza species revealed high levels of mitogenomic diversity, providing a foundation for future conservation and utilization of wild rice biodiversity.

Comparative analysis of codon usage patterns and phylogenetic implications of five mitochondrial genomes of the genus Japanagallia Ishihara, 1955 (Hemiptera, Cicadellidae, Megophthalminae)

Japanagallia is a genus of Cicadomorpha in the family of leafhoppers that are plant piercing-sucking insects, and it is difficult to distinguish by morphological characteristics. So far, only one complete mitochondrial genome data has been reported for the genus Japanagallia. Therefore, in order to better understand this group, we assembled and annotated the complete mitochondrial genomes of five Japanagallia species, and analyzed their codon usage patterns. Nucleotide composition analysis showed that AT content was higher than GC content, and the protein-coding sequences preferred to end with A/T at the third codon position. Relative synonymous codon usage analysis revealed most over-represented codon ends with A or T. Parity plot analysis revealed the codon usage bias of mitochondrial genes was influenced by both natural selection and mutation pressure. In the neutrality plot, the slopes of regression lines were < 0.5, suggesting that natural selection was playing a major role while mutation pressure was of minor importance. The effective number of codons showed that the codon usage bias between genes and genomes was low. Correspondence analysis revealed that the codon usage pattern differed among 13 protein-coding genes. Phylogenetic analyses based on three datasets using two methods (maximum likelihood and Bayesian inference), restored the Megophthalminae monophyly with high support values (bootstrap support values (BS) = 100, Bayesian posterior probability (PP) = 1). In the obtained topology, the seven Japanagallia species were clustered into a monophyletic group and formed a sister group with Durgade. In conclusion, our study can provide a reference for the future research on organism evolution, identification and phylogeny relationships of Japanagallia species.

Phylogenetic Analysis of Two New Mitochondrial Genomes of Singapora shinshana and Seriana bacilla from the Karst Region of Southwest China

Leafhoppers have been identified as a serious threat to different plants. To explore the characteristics of mitogenomes and reveal the phylogenetic positions of two species in the Typhlocybinae, complete mitogenomes of Singapora shinshana and Seriana bacilla were sequenced and annotated for the first time with lengths of 15,402 bp and 15,383 bp, respectively. The two mitogenomes contained 13 PCGs, 22 tRNA genes and 2 rRNA genes. The genome content, gene order, nucleotide composition, codon usage and amino acid composition are similar to those of other typical mitogenomes of Typhlocybinae. All 13 PCGs started with ATN codons, except for atp8 (TTA) and nad5 (TTG). All tRNAs were folded into a typical cloverleaf secondary structure, except for tRNA-Ser1 and tRNA-Val. Moreover, phylogenetic trees were constructed and analyzed based on all the PCGs from 42 mitogenomes using maximum likelihood (ML) and Bayesian inference (BI) methods. The results supported that eleven subfamilies are all monophyletic groups, S. shinshana and S. bacilla are members of Erythroneurini, but S. shinshana and the genus Empoascanara have a very close relationship with ((((Empoascanara sipra+ Empoascanara wengangensis) + Empoascanara dwalata) + Empoascanara gracilis) + S. shinshana), and S. bacilla is closely related to the genus Mitjaevia ((Mitjaevia dworakowskae + Mitjaevia shibingensis) + S. bacilla). These results provide valuable information for future study of evolutionary relationships in Typhlocybinae.

Structural Characteristics and Phylogenetic Analysis of the Mitochondrial Genomes of Four Krisna Species (Hemiptera: Cicadellidae: Iassinae)

Krisna species are insects that have piercing-sucking mouthparts and belong to the Krisnini tribe in the Iassinae subfamily of leafhoppers in the Cicadellidae family. In this study, we sequenced and compared the mitochondrial genomes (mitogenomes) of four Krisna species. The results showed that all four mitogenomes were composed of cyclic double-stranded molecules and contained 13 protein-coding genes (PCGs) and 22 and 2 genes coding for tRNAs and rRNAs, respectively. Those mitogenomes exhibited similar base composition, gene size, and codon usage patterns for the protein-coding genes. The analysis of the nonsynonymous substitution rate (Ka)/synonymous substitution rate (Ks) showed that evolution occurred the fastest in ND4 and the slowest in COI. 13 PCGs that underwent purification selection were suitable for studying phylogenetic relationships within Krisna. ND2, ND6, and ATP6 had highly variable nucleotide diversity, whereas COI and ND1 exhibited the lowest diversity. Genes or gene regions with high nucleotide diversity can provide potential marker candidates for population genetics and species delimitation in Krisna. Analyses of parity and neutral plots showed that both natural selection and mutation pressure affected the codon usage bias. In the phylogenetic analysis, all subfamilies were restored to a monophyletic group; the Krisnini tribe is monophyletic, and the Krisna genus is paraphyletic. Our study provides novel insights into the significance of the background nucleotide composition and codon usage patterns in the CDSs of the 13 mitochondrial PCGs of the Krisna genome, which could enable the identification of a different gene organization and may be used for accurate phylogenetic analysis of Krisna species.

Comparison of mitogenomes of three Petalocephala species (Hemiptera: Cicadellidae: Ledrinae) and their phylogenetic analysis

Ledrinae is a unique group of leafhoppers with a distinct appearance. Petalocephala is the largest Ledrinae genus that is difficult to identify except by dissecting the male genitals. To date, research on Ledrinae is relatively less compared with other leafhoppers. Therefore, to better understand this group, we sequenced and analyzed three complete Petalocephala mitochondrial genomes. We comparatively analyzed these general Petalocephala genomic features (including size, AT content, AT/GC skew, 13 protein-coding gene nucleotide compositions, etc.), and predicted 22 transfer RNA secondary structures. We obtained highly consistent phylogenetic results within Cicadellidae based on mitogenomic data using the maximum likelihood and Bayesian methods. Our results showed that all subfamilies were monophyletic and had a high node support rate, and there was a sister group relationship between Ledrinae and all other leafhopper groups. Furthermore, treehoppers were found to originate from leafhoppers and showed sister group relationships with Megophthalminae. Within Ledrinae, all phylogenetic trees supporting phylogenetic relationships were as follows: ([P. dicondylica + P. gongshanensis] + [Tituria pyramidata + [Ledra auditura + P. gongshanensis]]) Based on the complete mitogenome phylogenetic analysis and the comparison of morphological characteristics, we propose that Petalocephala is not monophyletic.

Mitogenomic phylogeny of Typhlocybinae (Hemiptera: Cicadellidae) reveals homoplasy in tribal diagnostic morphological traits

The subfamily Typhlocybinae is a ubiquitous, highly diverse group of mostly tiny, delicate leafhoppers. The tribal classification has long been controversial and phylogenetic methods have only recently begun to test the phylogenetic status and relationships of tribes. To shed light on the evolution of Typhlocybinae, we performed phylogenetic analyses based on 28 newly sequenced and 19 previously sequenced mitochondrial genomes representing all currently recognized tribes. The results support the monophyly of the subfamily and its sister-group relationship to Mileewinae. The tribe Zyginellini is polyphyletic with some included genera derived independently within Typhlocybini. Ancestral character state reconstruction suggests that some morphological characters traditionally considered important for diagnosing tribes (presence/absence of ocelli, development of hind wing submarginal vein) are homoplastic. Divergence time estimates indicate that the subfamily arose during the Middle Cretaceous and that the extant tribes arose during the Late Cretaceous. Phylogenetic results support establishment of a new genus, Subtilissimia Yan & Yang gen. nov., with two new species, Subtilissimia fulva Yan & Yang sp. nov. and Subtilissimia pellicula Yan & Yang sp. nov.; but indicate that two previously recognized species of Farynala distinguished only by the direction of curvature of the processes of the aedeagus are synonyms, that is, Farynala dextra Yan & Yang, 2017 equals Farynala sinistra Yan & Yang, 2017 syn. nov. A key to tribes of Typhlocybinae is provided.

Novel gene rearrangement in the mitochondrial genome of Anastatus fulloi (Hymenoptera Chalcidoidea) and phylogenetic implications for Chalcidoidea

The genus Anastatus comprises a large group of parasitoids, including several biological control agents in agricultural and forest systems. The taxonomy and phylogeny of these species remain controversial. In this study, the mitogenome of A. fulloi Sheng and Wang was sequenced and characterized. The nearly full-length mitogenome of A. fulloi was 15,692 bp, compromising 13 protein-coding genes (PCGs), 2 rRNA genes, 22 tRNA genes and a control region (CR). The total A + T contents were 83.83%, 82.18%, 87.58%, 87.27%, and 82.13% in the whole mitogenome, 13 PCGs, 22 tRNA genes, 2 rRNA genes, and CR, respectively. The mitogenome presented negative AT skews and positive GC skews, except for the CR. Most PCGs were encoded on the heavy strand, started with ATN codons, and ended with TAA codons. Among the 3736 amino acid-encoding codons, TTA (Leu1), CGA (Arg), TCA (Ser2), and TCT (Ser2) were predominant. Most tRNAs had cloverleaf secondary structures, except trnS1, with the absence of a dihydrouridine (DHU) arm. Compared with mitogenomes of the ancestral insect and another parasitoid within Eupelmidae, large-scale rearrangements were found in the mitogenome of A. fulloi, especially inversions and inverse transpositions of tRNA genes. The gene arrangements of parasitoid mitogenomes within Chalcidoidea were variable. A novel gene arrangement was presented in the mitogenome of A. fulloi. Phylogenetic analyses based on the 13 protein-coding genes of 20 parasitoids indicated that the phylogenetic relationship of 6 superfamilies could be presented as Mymaridae + (Eupelmidae + (Encyrtidae + (Trichogrammatidae + (Pteromalidae + Eulophidae)))). This study presents the first mitogenome of the Anastatus genus and offers insights into the identification, taxonomy, and phylogeny of these parasitoids.

A strategy for a high enrichment of insect mitochondrial DNA for mitogenomic analysis

Next-generation sequencing has dramatically fostered insect mitogenomic research in recent years. However, studies on the insect mitochondrial genome (mitogenome) assembly mainly rely on the sequencing data from total DNA, which is not cost-effective as a huge data from nuclear DNA are wasted. Besides, many mitogenomic studies require genomic information from individual organisms, whereas the DNA yield from small individual insects is too low to meet the sequencing requirements. Here, we describe a strategy for a high enrichment of insect mitochondrial DNA (mtDNA) using rolling circle amplification (RCA) technique. This strategy consists of standard DNA extraction, RCA enrichment, next-generation sequencing and mitogenome assembly. We have evaluated the performance of this strategy on nine insect species representing eight families of insecta, three other invertebrates, and even two vertebrate specimens. Results show that our strategy is especially suitable for insects, which allows almost all tested insect mtDNA contents to reach 80% and above. A further examination of enrichment efficiency of our strategy among different taxa shows that it is also applicable to other invertebrates and even some vertebrates such as Rhacophorus and ptyas species, although its enrichment efficiency in these groups is lower than that of insects. After treatment with our strategy, small flux sequencing data can realize the assembly of mitogenome with deep coverage, providing a solid base for subsequent mitogenome-based studies.

Characterization, Comparison of Four New Mitogenomes of Centrotinae (Hemiptera: Membracidae) and Phylogenetic Implications Supports New Synonymy

To explore the phylogenetic relationships of the subfamily Centrotinae from the mitochondrial genome data, four complete mitogenomes (Anchon lineatus, Anchon yunnanensis, Gargara genistae and Tricentrus longivalvulatus) were sequenced and analyzed. All the newly sequenced mitogenomes contain 37 genes. Among the 13 protein-coding genes (PCGs) of the Centrotinae mitogenomes, a sliding window analysis and the ratio of Ka/Ks suggest that atp8 is a relatively fast evolving gene, while cox1 is the slowest. All PCGs start with ATN, except for nad5 (start with TTG), and stop with TAA or the incomplete stop codon T, except for nad2 and cytb (terminate with TAG). All tRNAs can fold into the typical cloverleaf secondary structure, except for trnS1, which lacks the dihydrouridine (DHU) arm. The BI and ML phylogenetic analyses of concatenated alignments of 13 mitochondrial PCGs among the major lineages produce a well-resolved framework. Phylogenetic analyses show that Membracoidea, Smiliinae and Centrotinae, together with tribes Centrotypini and Leptobelini are recovered as well-supported monophyletic groups. The tribe Gargarini (sensu Wallace et al.) and its monophyly are supported.

Characterizing the Complete Mitochondrial Genome of Arma custos and Picromerus lewisi (Hemiptera: Pentatomidae: Asopinae) and Conducting Phylogenetic Analysis

We characterized the mitochondrial genome (mitogenome) and conducted phylogenetic analyses of 48 Hemiptera species by sequencing and analyzing the mitogenome of Arma custos (Fabricius) and Picromerus lewisi (Scott). The complete mitogenomes of the two predators were 16,024 bp and 19,587 bp in length, respectively, and it contained 37 classical genes, including 13 protein-coding genes (PCGs), 2 ribosomal RNA genes (rRNAs), and 22 transfer RNA genes (tRNAs), and a control region. Most PCGs in these predators use ATN as the start codon. This research revealed that the genes of the two natural enemy species have an A + T content of 75.40% and all tRNAs have a typical cloverleaf structure, with the exception of trnS1, which lacks a dihydrouridine arm. This is the first study to compare the mitochondrial genetic structure of two predatory insects; the mitochondrial genetic structure of individual predatory insects has been sequenced in previous studies. Here, phylogenetic analysis on the basis of amino acid and nucleotide sequences of 13 mitochondrial PCGs using Bayesian inference and maximum likelihood methods were conducted to generate similar tree topologies, which suggested that the two predators with close genetic relationships belong to Asopinae subfamily. Furthermore, the monophyly of the Pentatomoidea superfamily is well accepted despite limited taxon and species sampling. Finally, their complete mitogenome provided data to establish a predator-prey food web, which is the foundation of effective pest management. Our results also enhanced the database of natural enemy insects.

The complete mitochondrial genome of Aconurella prolixa (Lethierry 1885) (Hemiptera: Cicadellidae: Deltocephalinae: Chiasmini)

The complete mitochondrial genome of the widespread leafhopper species Aconurella prolixa (Hemiptera: Cicadellidae: Deltocephalinae: Chiasmini) was obtained via next-generation sequencing. This mitochondrial genome is 14,832 bp in length with the 37 classical eukaryotic mitochondrial genes and a control region. All 13 protein-coding genes (PCGs) are initiated with ATN, except ND5 uses TTG as the start codon, and terminate with TAA or TAG with the exception of COX2 and ND4 which use a single T residue as the stop codon. Twenty-one of the 22 transfer RNA (tRNAs) genes have the typical clover-leaf structure except for trnS1. Unlike some other species of deltocephalinae, no tRNA rearrangements were detected. The monophyly of Cicadellidae and Deltocephalinae, as well as the monophyly of Chiasmini, with a sister relationship between A. prolixa and (Exitianus indicus + Nephotettix cincticeps) is supported by Bayesian inference phylogenetic analyses based on 13 PCGs.

Characterization of the morphology and complete mitochondrial genomes of Eupteryx minusula and Eupteryx gracilirama (Hemiptera: Cicadellidae: Typhlocybinae) from Karst area, Southwest China

Background

The hemipteran insect family Cicadellidae (leafhoppers) includes >2,600 valid genera and >22,000 valid species worldwide, including >2,000 species in China. Typhlocybinae, second largest subfamilies of Cicadellidae, is widely distributed in the six major zoogeographic regions of the world, including >4,000 species worldwide and >1,000 species in China. Previously, morphological analysis are often effective to the way of taxonomy, but it did not combine with molecular biology. Therefore, morphology and mitochondrial genomes (mitogenomes) of two leafhopper species, Eupteryx (Eupteryx) minusculaLindberg, 1929 and Eupteryx (Stacla) graciliramaHou, Zhang & Huang, 2016 were studied and analyzed. This study analyzed the morphological and molecular characteristics of the two leafhoppers, and showed whether the results of the two identifications were consistent.

Methods

Based on the method of comparison, mitogenomes and morphology were analyzed to prove the relationship between the two leafhoppers.

Results

Although two focal species are classified in two different subgenera of the same genus, they still share many morphological features, such as the moderately produced crown fore margin; the milky yellow apical part of scutellum; the pronotum, basal triangles of scutellum, and forewing are dark with several colorless patches on the surface; the light yellow face, without any spots or stripes, and so on. The circular mitogenomes are 16,944 bp long in E. minuscula (GenBank: MN910279) and 17,173 bp long in E. gracilirama (GenBank: MT594485). All of the protein-coding genes are starting with ATN, except for some in mitogenome, which has a single T or TAN as a stop codon. All tRNAs have the typical cloverleaf-shaped structure except for trnS1 (AGN) (E. minuscula) which has a reduced DHU arm. Moreover, these two mitogenomes have trnR with an unpaired base in the acceptor stem. The phylogenetic relationships between E. minuscula and E. gracilirama in respect to related lineages were reconstructed using Maximum likelihood and Maximum parsimony analyses.

Discussion

The result showed that the tribe Typhlocybini is a sister to the tribes Erythroneurini and Empoascini, and five genera, Bolanusoides, Typhlocyba, Eupteryx, Zyginella and Limassolla are forming a single clade. E. minuscula and E. gracilirama are clustered together, supporting the monophyly of the genus Eupteryx. The above conclusions are consistent with the traditional classification of the subfamily.

Complete Mitochondrial Genomes of Metcalfa pruinosa and Salurnis marginella (Hemiptera: Flatidae): Genomic Comparison and Phylogenetic Inference in Fulgoroidea

The complete mitochondrial genomes (mitogenomes) of two DNA barcode-defined haplotypes of Metcalfa pruinosa and one of Salurnis marginella (Hemiptera: Flatidae) were sequenced and compared to those of other Fulgoroidea species. Furthermore, the mitogenome sequences were used to reconstruct phylogenetic relationships among fulgoroid families. The three mitogenomes, including that of the available species of Flatidae, commonly possessed distinctive structures in the 1702-1836 bp A+T-rich region, such as two repeat regions at each end and a large centered nonrepeat region. All members of the superfamily Fulgoroidea, including the Flatidae, consistently possessed a motiflike sequence (TAGTA) at the ND1 and trnS₂ junction. The phylogenetic analyses consistently recovered the familial relationships of (((((Ricaniidae + Issidae) + Flatidae) + Fulgoridae) + Achilidae) + Derbidae) in the amino acid-based analysis, with the placement of Cixiidae and Delphacidae as the earliest-derived lineages of fulgoroid families, whereas the monophyly of Delphacidae was not congruent between tree-constructing algorithms.

Structural features of the mitogenome of the leafhopper genus Cladolidia (Hemiptera: Cicadellidae: Coelidiinae) and phylogenetic implications in Cicadellidae

The first two complete mitogenomes of the leafhopper genus Cladolidia (C. biungulata and C. robusta) were sequenced and annotated to further explore the phylogeny of Cladolidia. Both the newly sequenced mitogenomes have a typical circular structure, with lengths of 15,247 and 15,376 bp and A + T contents of 78.2% and 78%, respectively. We identified a highly conserved genome organization in the two Cladolidia spp. through comparative analysis that included the following assessments: genome content, gene order, nucleotide composition, codon usage, amino acid composition, and tRNA secondary structure. Moreover, we detected the base heterogeneity of Cicadellidae mitogenomic data and constructed phylogenetic trees using the nucleotide alignments of 12 subfamilies of 58 leafhopper species. We noted a weak heterogeneity in the base composition among the Cicadellidae mitogenomes. Phylogenetic analyses showed that the monophyly of each subfamily was generally well supported in the family Cicadellidae; the main topology was as follows: (Deltocephalinae + (Treehoppers + ((Megophthalminae + (Macropsinae + (Hylicinae + (Coelidiinae +Iassinae)) + (Idiocerinae + (Cicadellinae + (Typhlocybinae + (Mileewinae + (Evacanthinae +Ledrinae)))))))))). Within Coelidiinae, phylogenetic analyses revealed that C. biungulata and C. robusta belong to Coelidiinae and the monophyly of Cladolidia is well supported. In addition, on the basis of complete mitogenome phylogenetic analysis and the comparison of morphological characteristics, we further confirm the genus Olidiana as a paraphyletic group, suggesting that the genus may need taxonomic revisions.

Two Complete Mitochondrial Genomes of Mileewinae (Hemiptera: Cicadellidae) and a Phylogenetic Analysis

Simple Summary

Mileewinae is a small subfamily of Cicadellidae containing about 160 described species, extensively distributed in the Oriental, Ethiopian and Neotropical regions. Some species are potential pests in agriculture and forestry. The classification of this group has been unstable over the past few decades. Currently, some controversies remain on the monophyly of Mileewinae and phylogenetic relationships of Mileewinae with other subfamilies. To provide further evidence toward answering these questions, two newly completed mitochondrial genomes of Mileewinae species (Mileewa rufivena and Ujna puerana) have been sequenced and analyzed. Results show these two mitochondrial genomes have quite similar structures and features. In phylogenetic analyses, Mileewinae formed a monophyletic group in Cicadellidae in all trees derived from maximum likelihood (ML) and Bayesian inference (BI) methods. In addition, Mileewinae has a closer phylogenetic relationship with Typhlocybinae compared to the Cicadellinae.

Abstract

More studies are using mitochondrial genomes of insects to explore the sequence variability, evolutionary traits, monophyly of groups and phylogenetic relationships. Controversies remain on the classification of the Mileewinae and the phylogenetic relationships between Mileewinae and other subfamilies remain ambiguous. In this study, we present two newly completed mitogenomes of Mileewinae (Mileewa rufivena Cai and Kuoh 1997 and Ujna puerana Yang and Meng 2010) and conduct comparative mitogenomic analyses based on several different factors. These species have quite similar features, including their nucleotide content, codon usage of protein genes and the secondary structure of tRNA. Gene arrangement is identical and conserved, the same as the putative ancestral pattern of insects. All protein-coding genes of U. puerana began with the start codon ATN, while 5 Mileewa species had the abnormal initiation codon TTG in ND5 and ATP8. Moreover, M. rufivena had an intergenic spacer of 17 bp that could not be found in other mileewine species. Phylogenetic analysis based on three datasets (PCG123, PCG12 and AA) with two methods (maximum likelihood and Bayesian inference) recovered the Mileewinae as a monophyletic group with strong support values. All results in our study indicate that Mileewinae has a closer phylogenetic relationship to Typhlocybinae compared to Cicadellinae. Additionally, six species within Mileewini revealed the relationship (U. puerana + (M. ponta + (M. rufivena + M. alara) + (M. albovittata + M. margheritae))) in most of our phylogenetic trees. These results contribute to the study of the taxonomic status and phylogenetic relationships of Mileewinae.

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.

The complete mitochondrial genomes of two erythroneurine leafhoppers (Hemiptera, Cicadellidae, Typhlocybinae, Erythroneurini) with assessment of the phylogenetic status and relationships of tribes of Typhlocybinae

The number and classification of tribes in the leafhopper subfamily Typhlocybinae are not yet fully clear, and molecular data has recently been used to help resolve the problem. In this study, the mitochondrial genomes of Mitjaeviashibingensis Chen, Song & Webb, 2020 and M.dworakowskae Chen, Song & Webb, 2020 of the tribe Erythroneurini (Cicadellidae, Typhlocybinae) were sequenced. Most protein-coding genes (PCGs) start with ATN and end with TAA or TAG, and the AT content of these three codons were found differ from previous results that show that the first codon has the highest incidence. Two rRNA genes are highly conserved, and the AT content in 16S is higher than that of 12S. The nucleotide diversity and genetic distance among 13 PCGs of the four tribes from Typhlocybinae show that Empoascini nucleotide diversity is significantly less than in the other three tribes, and have the largest distance from the others, while Typhlocybini and Zyginellini have the smallest distance, indicating that the relationship between the two is the closest. The nad2, nad4, nad4L, and nad5 genes have greater nucleotide diversity, showing potential for use as the main markers for species identification. The phylogenetic analysis yielded a well-supported topology with most branches receiving maximum support and a few branches pertaining to relationships within Zyginellini and Typhlocybini receiving lower support. The species of these two tribes are intertwined, and it was impossible to resolve them into separate branches. In addition, the tribes Empoascini and Erythroneurini were recovered as monophyletic, and Alebrini was placed at the base of the tree as the most primitive. These results are broadly in line with other molecular phylogenetical studies which differ from traditional morphological classification.

Mitogenomics of five Olidiana leafhoppers (Hemiptera: Cicadellidae: Coelidiinae) and their phylogenetic implications

Similar morphological characteristics and limited molecular data of Olidiana resulted in their unknown phylogenetic statuses and equivocal relationships. To further understand the genus Olidiana, we sequenced and annotated five Olidiana complete mitochondrial genomes (mitogenomes). Our results show that Olidiana mitogenomes range from 15,205 bp to 15,993 bp in length and include 37 typical genes (13 protein-coding genes, 22 tRNAs, and 2 rRNAs) and a control region. Their nucleotide composition, codon usage, features of control region, and tRNA secondary structures are similar to other members of Cicadellidae. We constructed the phylogenetic tree of Cicadellidae using the maximum likelihood (ML) and Bayesian inference (BI) methods based on all valid mitogenome sequences. The most topological structure of the obtained phylogenetic tree is consistent. Our results support the monophyletic relationships among 10 subfamilies within Cicadellidae and confirm Iassinae and Coelidiinae to be sister groups with high approval ratings. Interestingly, Olidiana was inferred as a paraphyletic group with strong support via both ML and BI analyses. These complete mitogenomes of five Olidiana species could be useful in further studies for species diagnosis, evolution, and phylogeny research within Cicadellidae.

Structural features and phylogenetic implications of Cicadellidae subfamily and two new mitogenomes leafhoppers

Complete mitochondrial genome sequences facilitate species identification and analyses of phylogenetic relationships. However, the available data are limited to the diverse and widespread insect family Cicadellidae. This study analyzes and summarizes the complete mitochondrial genome structure characteristics of 11 leafhopper subfamilies and two newly sequenced Typhlocybinae species, Empoascanara wengangensis and E. gracilis. Moreover, using 13PCGs and rRNA data to analyze the nucleotide diversity, evolution rate, and the phylogenetic relationship between the subfamilies of 56 species, verifying the taxonomic status analysis of E. wengangensis and E. gracilis. The analysis results show that the genome structures of the subfamilies and the newly sequenced two species are very similar, and the size of the CR region is significantly related to the repeat unit. However, in the entire AT-skews and CG-skews, the AT-skews of other subfamilies are all positive, and CG-skews are negative, while Empoascini of Typhlocybinae and Ledrinae are the opposite. Furthermore, among 13PCGs, the AT-skews of 13 species are all negative while CG-skews are positive, which from Empoascini in Typhlocybinae, Idiocerinae, Cicadellinae, Ledrinae, and Evacanthinae. Phylogenetic analysis shows that ML and PB analysis produce almost consistent topologies between different data sets and models, and some relationships are highly supported and remain unchanged. Mileewinae is a monophyletic group and is a sister group with Typhlocybinae, and the sister group of Evacanthinae is Ledrinae + Cicadellinae. Phylogenetic analysis grouped the two newly sequenced species with other species of Typhlocybinae, which was separated from other subfamilies, and all Erythroneurini insects gathered together. However, E. gracilis grouped into a single group, not grouped with species of the same genus (Empoascanara). This result does not match the traditional classification, and other nuclear genes or transcriptome genes may be needed to verify the result. Nucleotide diversity analysis shows that nad4 and nad5 may be evaluated as potential DNA markers defining the Cicadellidae insect species.

The complete mitochondrial genome of Trifida elongate and comparative analysis of 43 leafhoppers

In this study, the mitochondrial genome of Trifida elongate was sequenced, and comparative analysis of T. elongate and other 43 leafhoppers was performed based on the mitochondrial genome. The mitochondrial genome sequence length of T. elongate was 14,924 bp. It comprised 13 protein-coding genes, 2 ribosomal RNA (rRNA) genes, 22 transport RNA (tRNA) genes, and 1 non-coding control region. The control region is located between the rrnS and trnI genes, is characterized by two tandem repeats and three simple sequence repeats. Phylogenetic analysis showed that T. elongate is closely related with Bolanusoides shaanxiensis and Limassolla lingchuanensis (bootstrap value = 92% and posterior probabilities = 1). Analysis of synonymous and non-synonymous nucleotide substitutions showed that Ka/Ks value of the 13 protein-coding genes of 8 subfamily leafhoppers were less than 1 ranging from 0.0315 to 0.9928. atp8 had the highest Ka/Ks value whereas cox1 had the lowest Ka/Ks value. This study provides information on the structure and sequence characteristics of the mitochondrial genome of T. elongata. Typhlocybinae is clustered with (Cicadellinae+(Idiocerinae+(Mileewinae+(Nirvaninae+(Evacanthinae+Ledrinae))))).

Characterization of Two Complete Mitochondrial Genomes of Atkinsoniella (Hemiptera: Cicadellidae: Cicadellinae) and the Phylogenetic Implications

Simple Summary

Atkinsoniella is a large genus of almost 99 species across the world within the subfamily Cicadellinae, which is a large subfamily, comprising more than 2400 species of approximately 330 genera. Some of the Cicadellinae distributed worldwide are known as important agricultural pests. To better understand the mitogenomic characteristics of the genus Atkinsoniella and reveal phylogenetic relationships, the complete mitochondrial genomes of Atkinsoniella grahami and Atkinsoniella xanthonota were sequenced and comparatively analyzed in this study. The mitogenomes of these two Atkinsoniella species were found to be highly conserved, similarly to other Cicadellidae, except for the secondary structure of trnaS1, which formed a loop with the dihydrouridine (DHC) arm. This phenomenon has also been observed in other insect mitogenomes. Phylogenetic analyses, based on mitogenomes using both the maximum likelihood (ML) and Bayesian inference (BI) methods of three datasets, supported the monophyly of Cicadellinae, as well as the other subfamilies, and produced a well-resolved framework of Cicadellidae and valuable data for the phylogenetic study of Cicadellinae.

Abstract

The complete mitochondrial genomes of Atkinsoniella grahami and Atkinsoniella xanthonota were sequenced. The results showed that the mitogenomes of these two species are 15,621 and 15,895 bp in length, with A+T contents of 78.6% and 78.4%, respectively. Both mitogenomes contain 13 protein-coding genes (PCGs), 22 transfer RNA genes (tRNAs), 2 ribosomal RNA genes (rRNAs), and a control region (CR). For all PCGs, a standard start ATN codon (ATT, ATG, or ATA) was found at the initiation site, except for ATP8, for which translation is initiated with a TTG codon. All PCGs terminate with a complete TAA or TAG stop codon, except for COX2, which terminates with an incomplete stop codon T. All tRNAs have the typical cloverleaf secondary structure, except for trnS, which has a reduced dihydrouridine arm. Furthermore, these phylogenetic analyses were reconstructed based on 13 PCGs and two rRNA genes of 73 mitochondrial genome sequences, with both the maximum likelihood (ML) and Bayesian inference (BI) methods. The obtained mitogenome sequences in this study will promote research into the classification, population genetics, and evolution of Cicadellinae insects in the future.

Comparative analysis of mitochondrial genomes of Nirvanini and Evacanthini (Hemiptera: Cicadellidae) reveals an explicit evolutionary relationship

Mitogenomes of five leafhopper species, Chudania hellerina and Concaveplana rufolineata in Nirvanini, Carinata rufipenna, Evacanthus danmainus and E. heimianus representing Evacanthini, were sequenced. The lengths of these five mitogenomes range from 15,044 (C. hellerina) to 15,680 bp (E. heimianus). All five mitogenomes exhibit similar base composition, gene size and codon usage of protein-coding genes. All 22 tRNA genes have typical cloverleaf secondary structures, except for trnS1 (AGN) which appears to lack the dihydrouridine arm. The two included Nirvanini species employ the anticodon TCT instead of the commonly used GCT in trnS1 (AGN). Genes nad2, atp8 and nad6 were highly variable while cox1 and cob showed the lowest nucleotide diversity. Phylogenetic analyses of two concatenated nucleotide datasets, incorporating the newly sequenced taxa and other available membracoid mitogenomes, recovered each included leafhopper subfamily as monophyletic with evacanthine tribes Nirvanini and Evacanthini forming monophyletic sister clades. A relationship among Evacanthinae, Cicadellinae and Typhlocybinae received moderate branch support.

The Complete Mitochondrial Genome of Four Hylicinae (Hemiptera: Cicadellidae): Structural Features and Phylogenetic Implications

To reveal mtgenome characterizations and reconstruct phylogenetic relationships of Hylicinae, the complete mtgenomes of four hylicine species, including Nacolus tuberculatus, Hylica paradoxa, Balala fujiana, and Kalasha nativa, were sequenced and comparatively analyzed for the first time. We also carried out the richest (11) subfamily sampling of Cicadellidae to date, and reconstructed phylogenetic relationships of Membracoidea among 61 species based on three datasets using maximum likelihood and Bayesian inference analyses. All new sequenced mtgenomes are molecules ranging from 14,918 to 16,221 bp in length and are double stranded, circular in shape. The gene composition and arrangement of these mtgenomes are consistent with members of Membracoidea. Among 13 protein-coding genes, most show typical ATN start codons and TAR (TAA/TAG) or an incomplete stop codon T-, and several genes start by TTG/GTG. Results of the analysis for sliding window, nucleotide diversity, and nonsynonymous substitution/synonymous substitution indicate cox1 is a comparatively slower-evolving gene while atp8 is the fastest gene. In line with previous researches, phylogenetic results indicate that treehopper families are paraphyletic with respect to family Cicadellidae and also support the monophyly of all involved subfamilies including Hylicinae. Relationships among the four hylicine genera were recovered as (Hylica + (Nacolus + (Balala + Kalasha))).

Comparative mitochondrial genomes of four species of Sinopodisma and phylogenetic implications (Orthoptera, Melanoplinae)

In this study, the whole mitochondrial genomes (mitogenomes) from four species were sequenced. The complete mitochondrial genomes of Sinopodisma pieli, S. houshana, S. qinlingensis, and S. wulingshanensis are 15,857 bp, 15,818 bp, 15,843 bp, and 15,872 bp in size, respectively. The 13 protein-coding genes (PCGs) begin with typical ATN codons, except for COXI in S. qinlingensis, which begins with ACC. The highest A+T content in all the sequenced orthopteran mitogenomes is 76.8% (S. qinlingensis), followed by 76.5% (S. wulingshanensis), 76.4% (S. pieli) and 76.4% (S. houshana) (measured on the major strand). The long polythymine stretches (T-stretch) in the A+T-rich region of the four species are not adjacent to the trnI locus but are inside the stem-loop sequences on the major strand. Moreover, several repeated elements are found in the A+T-rich region of the four species. Phylogenetic analysis based on 53 mitochondrial genomes using Bayesian Inference (BI) and Maximum Likelihood (ML) revealed that Melanoplinae (Podismini) was a monophyletic group; however, the monophyly of Sinopodisma was not supported. These data will provide important information for a better understanding of the phylogenetic relationship of Melanoplinae.

Characterization of Two Complete Mitochondrial Genomes of Ledrinae (Hemiptera: Cicadellidae) and Phylogenetic Analysis

Mitochondrial genomes are widely used for investigations into phylogeny, phylogeography, and population genetics. More than 70 mitogenomes have been sequenced for the diverse hemipteran superfamily Membracoidea, but only one partial and two complete mtgenomes mitochondrial genomes have been sequenced for the included subfamily Ledrinae. Here, the complete mitochondrial genomes (mitogenomes) of two additional Ledrinae species are newly sequenced and comparatively analyzed. Results show both mitogenomes are circular, double-stranded molecules, with lengths of 14,927 bp (Tituria sagittata) and 14,918 bp (Petalocephala chlorophana). The gene order of these two newly sequenced Ledrinae is highly conserved and typical of members of Membracoidea. Similar tandem repeats in the control region were discovered in Ledrinae. Among 13 protein-coding genes (PCGs) of reported Ledrinae mitogenomes, analyses of the sliding window, nucleotide diversity, and nonsynonymous substitution (Ka)/synonymous substitution (Ks) indicate atp8 is a comparatively fast-evolving gene, while cox1 is the slowest. Phylogenetic relationships were also reconstructed for the superfamily Membracoidea based on expanded sampling and gene data from GenBank. This study shows that all subfamilies (sensu lato) are recovered as monophyletic. In agreement with previous studies, these results indicate that leafhoppers (Cicadellidae) are paraphyletic with respect to the two recognized families of treehoppers (Aetalionidae and Membracidae). Relationships within Ledrinae were recovered as (Ledra + (Petalocephala + Tituria)).

Characterization of the Complete Mitochondrial Genome of Drabescus ineffectus and Roxasellana stellata (Hemiptera: Cicadellidae: Deltocephalinae: Drabescini) and Their Phylogenetic Implications

To explore the mitogenome characteristics and shed light on the phylogenetic relationships and molecular evolution of Drabescini species, we sequenced and analyzed the complete mitochondrial genome of two species including Drabescus ineffectus and Roxasellana stellata. The complete mitogenomes of D. ineffectus and R. stellata are circular, closed and double-stranded molecules with a total length of 15744 bp and 15361 bp, respectively. These two newly sequenced mitogenomes contain the typical 37 genes. Most protein-coding genes (PCGs) began with the start codon ATN and terminated with the terminal codon TAA or TAG, with an exception of a special initiation codon of ND5, which started with TTG, and an incomplete stop codon T-- was found in the Cytb, COX2, ND1 and ND4. All tRNAs could be folded into the canonical cloverleaf secondary structure except for the trnS1, which lacks the DHU arm and is replaced by a simple loop. The multiple tandem repeat units were found in A + T-control region. The sliding window, Ka/Ks and genetic distance analyses indicated that the ATP8 presents a high variability and fast evolutionary rate compared to other PCGs. Phylogenetic analyses based on three different datasets (PCG123, PCG12R and AA) using both Bayesian inference (BI) and maximum likelihood (ML) methods showed strong support for the monophyly of Drabescini.

The Phylogenetic Implications of the Mitochondrial Genomes of Macropsis notata and Oncopsis nigrofasciata

Macropsinae are forest pests that feed on woody plants. They can damage the growth of trees and crops, and some species can also spread plant pathogens. Due to their widespread effects, these leafhoppers are of great economic significance, which is why there is a need to study their genomes. To fill the gap in the mitochondrial genomic data of the subfamily Macropsinae, we sequenced the complete mitochondrial genomes of Macropsis notata and Oncopsis nigrofasciata (which were 16,323 and 15,927 bp long, respectively). These two species are representative species of the leafhoppers group (Cicadellidae); the mitochondrial genomes of these species range from a length of 15,131 bp (Trocnadella arisana) to 16,811 bp (Parocerus laurifoliae). Both mitogenomes contained 37 typical insect mitochondrial genes and a control region; there were no long non-coding sequences. The genes within the mitogenome were very compact. The mitogenomes from both species contained two kinds of parallel repeat units in the control region. The whole mitogenomes of Macropsinae showed a heavy AT nucleotide bias (M. notata 76.8% and O. nigrofasciata 79.0%), a positive AT Skew (0.15 and 0.12), and a negative GC Skew (-0.14 and -0.08). Upon comparative ML and BI analysis, some clade relationships were consistent among the six trees. Most subfamilies were reconstructed into monophyletic groups with strong support in all analyses, with the exception of Evacanthinae and Cicadellinae. Unlike the results of previous research, it was shown that although all Deltocephalinae species are grouped into one clade, they were not the sister group to all other leafhoppers. Further, Cicadellinae and Evacanthinae were occasionally reconstructed as a polyphyletic and a paraphyletic group, respectively, possibly due to the limited numbers of samples and sequences. This mitogenome information for M. notata and O. nigrofasciata could facilitate future studies on the mitogenomic diversity and evolution of the related Membracoidea, and eventually help to control their effects on plants for the betterment of society at large.

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.

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.

Complete mitogenome of Olidiana ritcheriina (Hemiptera: Cicadellidae) and phylogeny of Cicadellidae

BACKGROUND

Coelidiinae, a relatively large subfamily within the family Cicadellidae, includes 129 genera and ∼1,300 species distributed worldwide. However, the mitogenomes of only two species (Olidiana sp. and Taharana fasciana) in the subfamily Coelidiinae have been assembled. Here, we report the first complete mitogenome assembly of the genus Olidiana.

METHODS

Specimens were collected from Wenxian County (Gansu Province, China) and identified on the basis of their morphology. Mitogenomes were sequenced by next-generation sequencing, following which an NGS template was generated, and this was confirmed using polymerase chain reaction and Sanger sequencing. Phylogenic trees were constructed using maximum likelihood and Bayesian analyses.

RESULTS

The mitogenome of O. ritcheriina was 15,166 bp long, with an A + T content of 78.0%. Compared with the mitogenome of other Cicadellidae sp., the gene order, gene content, gene size, base composition, and codon usage of protein-coding genes (PCGs) in O. ritcheriina were highly conserved. The standard start codon of all PCGs was ATN and stop codon was TAA or TAG; COII, COIII, and ND4L ended with a single T. All tRNA genes showed the typical cloverleaf secondary structure, except for trnSer, which did not have the dihydrouridine arm. Furthermore, the secondary structures of rRNAs (rrnL and rrnS) in O. ritcheriina were predicted. Overall, five domains and 42 helices were predicted for rrnL (domain III is absent in arthropods), and three structural domains and 27 helices were predicted for rrnS. Maximum likelihood and Bayesian analyses indicated that O. ritcheriina and other Coelidiinae members were clustered into a clade, indicating the relationships among their subfamilies; the main topology was as follows: (Deltocephalinae + ((Coelidiinae + Iassinae) + ((Typhlocybinae + Cicadellinae) + (Idiocerinae + (Treehopper + Megophthalminae))))). The phylogenetic relationships indicated that the molecular taxonomy of O. ritcheriina is consistent with the current morphological classification.

Structural and phylogenetic implications of the complete mitochondrial genome of Ledra auditura

We sequenced and annotated the first complete mitochondrial genome (mitogenome) of Ledra auditura (Hemiptera: Cicadellidae: Ledrinae) and reconstructed phylogenetic relationships among 47 species (including 2 outgroup species) on the basis of 3 datasets using maximum likelihood (ML) and Bayesian inference (BI) analyses. The complete L. auditura mitogenome (length, 16,094 bp) comprises 37 genes [13 protein-coding genes (PCGs), 22 tRNAs, and 2 rRNAs], 1 control region, and 2 long non-coding regions. The first long non-coding region (length, 211 bp) is located between tRNA-I and tRNA-Q and the second region (length, 994 bp) between tRNA-S2 and ND1. All PCGs show ATN (Met/Ile) as their start codon and TAR as their stop codon. Except tRNA-S1 (AGN), which lacks the dihydrouridine arm, all tRNAs can fold into the typical cloverleaf secondary structure. The complete L. auditura mitogenome shows a base composition bias of 76.3% A + T (A = 29.9%, T = 46.4%, G = 13.3%, and C = 10.5%), negative AT skew of -0.22, and positive GC skew of 0.12. In ML and BI analyses, L. auditura was clustered with Evacanthus heimianus (Hemiptera: Cicadellidae: Evacanthinae) with strong branch support.

Characterization of the complete mitochondrial genomes of two species of the genus Aphaena Guérin-Méneville (Hemiptera: Fulgoridae) and its phylogenetic implications

The complete mitochondrial genomes (mitogenomes) of Aphaena (Callidepsa) amabilis and Aphaena (Aphaena) discolor nigrotibiata were sequenced. The mitogenomes of these two species are 16,237 bp and 16,116 bp in length with an A + T content of 77.9% and 77.0%, respectively. Each contains 13 protein-coding genes (PCGs), 22 transfer RNA genes (tRNAs), two ribosomal RNA genes (rRNAs) and a control region (A + T-rich region). All PCGs initiate with the standard start codon of ATN and terminate with the complete stop codon of TAA or TAG except for atp6, where nad1 ends with an incomplete T codon. All tRNAs have the typical clover-leaf structure except for trnS1 and trnV which have a reduced DHU arm. Moreover, these two mitogenomes have trnL2, trnR and trnT with an unpaired base in the acceptor stem. The putative A + T-rich region includes multiple types of tandem repeat regions. These phylogenetic analyses are reconstructed based on 13 protein-coding genes of 25 auchenorrhynchan mitogenomes, with both maximum likelihood and Bayesian analyses yielding robust identical phylogenetic trees. These results support a monophyletic Auchenorrhyncha and the relationship (Pyrops + (Lycorma + Aphaena)) within Fulgoridae.

Complete mitochondrial genomes of Bittacus strigosus and Panorpa debilis and genomic comparisons of Mecoptera

Mitochondrial genomes play a significant role in reconstructing phylogenetic relationships and revealing molecular evolution of insects. However, only four mitochondrial genomes were reported in Mecoptera to date. Here, we obtained two new complete mitochondrial genomes of the hangingfly Bittacus strigosus Hagen, 1861 and the scorpionfly Panorpa debilis Westwood, 1846. The results show that the complete mitogenome sequences of B. strigosus and P. debilis are 15,825 and 17,018 bp, respectively, both containing 37 genes and one control region. The mecopteran mitogenomes are highly similar in A + T bias, AT-skew, and GC-skew. Tandem repeats of the control region were discovered in Mecoptera for the first time. The sliding window, genetic distance, and Ka/Ks ratio analyses indicate the purifying selection of 13 protein-coding genes, the lowest evolutionary rate of cox1, and the highest sequence variability of atp8. Considering the sufficiently large size, fast evolution, and high ratio of Ka/Ks, nad4L and nad6 are regarded as potential markers for future phylogenetic analyses, population genetics, and species delimitations in Mecoptera. The phylogenetic relationships were reconstructed for four families of Mecoptera based on all six available mitogenomes using Bayesian inference and maximum likelihood methods. The phylogeny is presented as Boreidae + (Nannochoristidae + (Bittacidae + Panorpidae)).

Mitochondrial and chloroplast genomes provide insights into the evolutionary origins of quinoa (Chenopodium quinoa Willd.)

Quinoa has recently gained international attention because of its nutritious seeds, prompting the expansion of its cultivation into new areas in which it was not originally selected as a crop. Improving quinoa production in these areas will benefit from the introduction of advantageous traits from free-living relatives that are native to these, or similar, environments. As part of an ongoing effort to characterize the primary and secondary germplasm pools for quinoa, we report the complete mitochondrial and chloroplast genome sequences of quinoa accession PI 614886 and the identification of sequence variants in additional accessions from quinoa and related species. This is the first reported mitochondrial genome assembly in the genus Chenopodium. Inference of phylogenetic relationships among Chenopodium species based on mitochondrial and chloroplast variants supports the hypotheses that 1) the A-genome ancestor was the cytoplasmic donor in the original tetraploidization event, and 2) highland and coastal quinoas were independently domesticated.

Characterization and phylogenetic analysis of the complete mitochondrial genome of the medicinal fungus Laetiporus sulphureus

The medicinal fungus Laetiporus sulphureus is widely distributed worldwide. To screen for molecular markers potentially useful for phylogenetic analyses of this species and related species, the mitochondrial genome of L. sulphureus was sequenced and assembled. The complete circular mitochondrial genome was 101,111 bp long, and contained 38 protein-coding genes (PCGs), 2 rRNA genes, and 25 tRNA genes. Our BLAST search aligned about 6.1 kb between the mitochondrial and nuclear genomes of L. sulphureus, indicative of possible gene transfer events. Both the GC and AT skews in the L. sulphureus mitogenome were negative, in contrast to the other seven Polyporales species tested. Of the 15 PCGs conserved across the seven species of Polyporales, the lengths of 11 were unique in the L. sulphureus mitogenome. The Ka/Ks of these 15 PCGs were all less than 1, indicating that PCGs were subject to purifying selection. Our phylogenetic analysis showed that three single genes (cox1, cob, and rnl) were potentially useful as molecular markers. This study is the first publication of a mitochondrial genome in the family Laetiporaceae, and will facilitate the study of population genetics and evolution in L. sulphureus and other species in this family.

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.