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

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

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.

Description and complete mitochondrial genome of Atkinsoniella zizhongi sp. nov. (Hemiptera: Cicadellidae: Cicadellinae) from China and its phylogenetic implications

A new species, Atkinsoniella zizhongi sp. nov. of the subfamily Cicadellinae, was described and illustrated from China. The new species is similar to A. nigrominiatula (Jacobi, 1944), A. limba Kuoh, 1991, A. dormana Li, 1992, A. peaka Yang, Meng et Li, 2017, and A. divaricata Yang, Meng et Li, 2017. But the characteristics of aedeagus and pygofer process can be used to distinguish them easily. The complete mitochondrial genome of the paratype was sequenced and assembled. The mitogenome of A. zizhongi sp. nov. was 16,483 bp in length, with an A+T content of 75.9%, containing 37 typical genes and a control region (CR). The gene order was consistent with the inferred insect ancestral mitochondrial genome. All of the PCGs were determined to have the typical stop codon TAA or TAG, while COX2 and ND5 ended with incomplete termination codons T and TA, respectively. In addition, phylogenetic trees were reconstructed based on PCGs and rRNAs using both the maximum likelihood (ML) and Bayesian inference (BI) methods. The results showed that the intergeneric and interspecific relationships within the subfamily Cicadellinae were completely consistent in all of the phylogenetic trees, except that the different interspecific relationships within the genus Bothrogonia were detected in the ML analysis based on the amino acid sequences. This study enriches the species diversity of Cicadellinae and further promotes research on its phylogeny.

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.

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.

Comparative Analysis of Mitochondrial Genomes among Twelve Sibling Species of the Genus Atkinsoniella Distant, 1908 (Hemiptera: Cicadellidae: Cicadellinae) and Phylogenetic Analysis

Simple Summary

Atkinsoniella is a large genus of 98 species across the world and 88 species recorded in China within the globally distributed subfamily Cicadellinae, which is phytophagous, and some of which have been reported as important agricultural pests. Some Atkinsoniella species are very similar in morphological characteristics, making accurate identification at species level confusing. To provide further evidence toward understanding the relationships within the genus Atkinsoniella and subfamily Cicadellinae, mitogenomes of 12 Atkinsoniella sibling species were obtained and annotated. Their characteristics were comparatively analyzed. In addition, the comprehensive phylogenetic relationship within the subfamily Cicadellinae was determined based on three mitochondrial datasets using both the maximum-likelihood (ML) and Bayesian inference (BI) methods. The results suggested that the genus Atkinsoniella was recovered as a monophyletic group. The branches of the 12 newly sequenced species were clearly separated, with most nodes receiving strong support in all analyses, indicating that mitogenomics is an effective method for identifying closely related species and understanding their phylogenetic and evolutionary relationships.

Abstract

The herbivorous leafhopper genus Atkinsoniella Distant, 1908 (Hemiptera: Cicadellidae: Cicadellinae), a large genus of subfamily Cicadellinae, consists of 98 valid species worldwide and 88 species recorded in China. Some species of the genus are very similar in morphological characteristics, so they are difficult to identify accurately. In this study, 12 mitochondrial genomes of Atkinsoniella species with similar morphological characteristics were first obtained through high-throughput sequencing, which featured a typical circular molecule of 15,034–15,988 bp in length. The arrangement and orientation of 37 genes were identical to those of typical Cicadellidae mitogenomes. The phylogenetic relationship within the subfamily Cicadellinae was reconstructed using maximum-likelihood (ML) and Bayesian inference (BI) methods based on three concatenated datasets. The topological structures of the six obtained phylogenetic trees were highly consistent. The results suggested that Atkinsoniella was recovered as a monophyletic group and emerged as a sister group with the monophyletic clade of Bothrogonia, Paracrocampsa (part), and Draeculacephala (part). The branches of the 12 newly sequenced species were clearly separated, with most nodes receiving strong support in all analyses. In addition, the key to the 12 Atkinsoniella species was provided to identify species according to morphological characteristics. This study further promotes research on the classification, genetics, evolution, and phylogeny of the genus Atkinsoniella and subfamily Cicadellinae.

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.

Structural Features and Phylogenetic Implications of 11 New Mitogenomes of Typhlocybinae (Hemiptera: Cicadellidae)

To explore the characteristics of mitogenomes and discuss the phylogenetic relationships and molecular evolution of the six tribes within Typhlocybinae, 11 complete mitogenomes are newly sequenced and comparatively analyzed. In all of these complete mitogenomes, the number and order of the genes are highly conserved in overall organization. The PCGs initiate with ATN/TTG/GTG and terminate with TAA/TAG/T. Almost all tRNAs are folded into the typical clover-leaf secondary structure. The control region is always variable in length and in numbers of multiple tandem repeat units. The atp8 and nad2 exhibits the highest evolution rate among all the PCGs. Phylogenetic analyses based on whole mitogenome sequences, with three different datasets, using both maximum likelihood and Bayesian methods, indicate the monophyly of Typhlocybinae and its inner tribes, respectively, except for Typhlocybini and Zyginellini that are paraphyletic. Finally, we confirm that Erythroneurini is a subtribe of Dikraneurini.

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

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