Mitochondrial genomic variation and phylogenetic relationships of three groups in the genus Scaphoideus (Hemiptera: Cicadellidae: Deltocephalinae)

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

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 of Cyclopelta obscura (Lepeletier & Serville, 1825) (Hemiptera: Pentatomoidea: Dinidoridae) and phylogenetic analysis of Pentatomoidea species

Cyclopelta obscura is a crop pest, which mainly damages legumes, especially Robinia pseudoacacia and Cercis chinensis. In recent years, many mitochondrial and nuclear DNA sequences of C. obscura have been sequenced and used for phylogenetic inference. However, the complete mitogenome has not been reported yet and studies on the phylogenetic relationships within Dinidoridae are rare. In this study, we sequenced the mitogenome of C. obscura and conducted comparative mitogenomic analyses of seven Dinidoridae species based on several different factors. The length of the mitogenome is 15,426 bp, which includes 37 typical mitochondrial genes (13 protein-coding genes (PCGs), 22 tRNAs, and 2 rRNAs) and a control region (796 bp long), as well as 13 intergenic spacers and 8 overlapping regions. Most PCGs of C. obscura began with the classical start codon ATN, while cox1 and nad4l used TTG, and nad1 used GTG. The Dinidoridae mitogenomes are highly conserved in terms of nucleotide composition, the codon usage of PCGs, and the secondary structure of tRNA. Phylogenetic analysis based on four datasets with two methods recovered the Dinidoridae as a monophyletic group with strong support values. All results indicate that Dinidoridae formed a sister group to Tessaratomidae, and (Tessaratomidae + Dinidoridae) formed a sister group to Cydnidae in most of the phylogenetic trees. Additionally, seven species within the Dinidoridae, we observed the following relationship: (Eumenotes sp. + (Cyclopelta parva + C. obscura)) + ((Megymenum gracilicorne + Megymenum brevicorne) + (Coridius chinensis + Coridius brunneu)).

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.

Comparative Mitogenomics of Flesh Flies: Implications for Phylogeny

Flesh flies (Diptera: Sarcophagidae) represent a rapid radiation belonging to the Calyptratae. With more than 3000 known species, they are extraordinarily diverse in terms of their breeding habits and are therefore of particular importance in human and veterinary medicine, forensics, and ecology. To better comprehend the phylogenetic relationships and evolutionary characteristics of the Sarcophagidae, we sequenced the complete mitochondrial genomes of five species of flesh flies and performed mitogenomic comparisons amongst the three subfamilies. The mitochondrial genomes match the hypothetical condition of the insect ancestor in terms of gene content and gene arrangement. The evolutionary rates of the subfamilies of Sarcophagidae differ significantly, with Miltogramminae exhibiting a higher rate than the other two subfamilies. The monophyly of the Sarcophagidae and each subfamily is strongly supported by phylogenetic analysis, with the subfamily-level relationship inferred as (Sarcophaginae, (Miltogramminae, Paramacronychiinae)). This study suggests that phylogenetic analysis based on mitochondrial genomes may not be appropriate for rapidly evolving groups such as Miltogramminae and that the third-codon positions could play a considerable role in reconstructing the phylogeny of Sarcophagidae. The protein-coding genes ND2 and ND6 have the potential to be employed as DNA markers for species identification and delimitation in flesh flies.

Administration time-dependent effects of poly (I:C) on antioxidant and immune responses along the diurnal time scale in zebrafish

The circadian clock of vertebrates regulates many biological processes, including the immune system. This paper investigated whether responsiveness to poly (I:C), a synthetic analog of double-stranded RNA used as an immunostimulant, exhibits day/night differences in zebrafish. Fish were intraperitoneally (IP) injected with either phosphate-buffered saline (PBS) or poly (I:C) at two different time points: "Zeitgeber Time" (ZT) 4 (day) and ZT16 (night). Then, 6 h later, fish were euthanized, and tissue samples (skin, liver and kidney) were collected. A control group (intact fish) was also sampled at the same time points. The effect of poly (I:C) on the expression of antioxidant and immune genes was time-of-day-dependent, and the response was stronger following poly (I:C) administration in the day than at night. Time-dependent differences were observed for some genes in the PBS and control groups. However, these differences were tissue-specific. In liver, almost all the genes were affected by time of day. In kidney, poly (I:C) affected the expression of all the gene markers regardless of administration time. These findings highlight the importance of considering the time to administer poly (I:C) when evaluating the fish immune response.

Characterization of the complete chloroplast genome of Carya laciniosa (F.Michx.) G.Don (Juglandaceae)

Carya laciniosa (F.Michx.) G.Don is a slow-growing valuable tree in the family Juglandaceae. In this research, we first reported the chloroplast genome sequence of C. laciniosa. The chloroplast genome was a circular form with 160,832 bp in size, comprising four subregions: a large single-copy (LSC) region of 90,065 bp, a small single-copy (SSC) region of 18,791 bp, and two copies of inverted repeats (IRs; IRa and IRb) of 25,988 bp each. A total of 132 genes were annotated in the chloroplast genome of C. laciniosa, including 87 protein-coding genes, 37 tRNA genes, and 8 rRNA genes. The overall GC content of the whole genome was 36.14%. Phylogenetic analysis suggested that C. laciniosa was closely related to C. ovata.

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.

Complete mitochondrial genomes of four species of praying mantises (Dictyoptera, Mantidae) with ribosomal second structure, evolutionary and phylogenetic analyses

Praying mantises are distributed all over the world. Though some Mantodea mitogenomes have been reported, an evolutionary genomic and phylogenetic analysis study lacks the latest taxonomic system. In the present study, four new mitogenomes were sequenced and annotated. Deroplatys truncate, D. lobate, Amorphoscelis chinensis and Macromantis sp. belong to Deroplatyidae, Amorphoscelidae and Photinaidae family, respectively. Our results indicated that the ATP8 gene may be lost in D. truncate and D. lobata mt genome, and four tRNA genes have not been found in D. truncate, D. lobata and Macromantis sp. A dN/dS pair analysis was conducted and it was found that all genes have evolved under purifying selection. Furthermore, we tested the phylogenetic relationships between the eight families of the Mantodea, including 35 species of praying Mantis. Based on the complete mitochondrial genome data, it was also suggested as sister to Deroplatyidae + Mantidae, Metallyticus sp., the only representative of Metallyticidae, is sister to the remaining mantises. Our results support the taxonomic system of Schwarz and Roy and are consistent with previous studies.

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 of the First Complete Mitochondrial Genome of Cyphonocerinae (Coleoptera: Lampyridae) with Implications for Phylogeny and Evolution of Fireflies

Complete mitochondrial genomes are valuable resources for phylogenetics in insects. The Cyphonoceridae represents an important lineage of fireflies. However, no complete mitogenome is available until now. Here, the first complete mitochondrial genome from this subfamily was reported, with Cyphonocerus sanguineus klapperichi as a representative. The mitogenome of C. sanguineus klapperichi was conserved in the structure and comparable to that of others in size and A+T content. Nucleotide composition was A+T-biased, and all genes exhibited a positive AT-skew and negative GC-skew. Two types of tandem repeat sequence units were present in the control region (136 bp × 2; 171 bp × 2 + 9 bp). For reconstruction of Lampyridae's phylogeny, three different datasets were analyzed by both maximum likelihood (ML) and Bayesian inference (BI) methods. As a result, the same topology was produced by both ML analysis of 13 protein-coding genes and 2rRNA and BI analysis of 37 genes. The results indicated that Lampyridae, Lampyrinae, Luciolinae (excluding Emeia) were monophyletic, but Ototretinae was paraphyletic, of which Stenocladius was recovered as the sister taxon to all others, while Drilaster was more closely related to Cyphonocerinae; Phturinae + Emeia were included in a monophyletic clade, which comprised sister groups with Lampyridae. Vesta was deeply rooted in the Luciolinae.

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.

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.

Characterization of the Complete Mitochondrial Genomes of Two Species with Preliminary Investigation on Phylogenetic Status of Zyginellini (Hemiptera: Cicadellidae: Typhlocybinae)

To explore the characteristics of mitogenomes and reveal phylogenetic relationships of the tribes of Zyginellini and Typhlocybini in Typhlocybinae, mitogenomes of two species of the Zyginellini, Parazyginella tiani and Limassolla sp., were sequenced. Mitogenomes of both species contain 13 protein-coding genes (PCGs), 22 transfer RNA genes (tRNAs), two ribosomal RNA genes (rRNAs) and a large non-coding region (A + T-rich region). These characteristics are similar to other Membracoidea mitogenomes. All PCGs initiate with the standard start codon of ATN and terminate with the complete stop codon of TAA/G or with an incomplete T codon. All tRNAs have the typical clover-leaf structure, except trnS1 which has a reduced DHU arm and the acceptor stem of trnR is 5 or 6 bp in some species, an unusual feature here reported for the first time in Typhlocybinae. The A + T-rich region is highly variable in length and in numbers of tandem repeats present. Our analyses indicate that nad6 and atp6 exhibit higher evolutionary rates compared to other PCGs. Phylogenetic analyses by both maximum likelihood and Bayesian methods based on 13 protein-coding genes of 12 species of Typhlocybinae suggest that Zyginellini are paraphyletic with respect to Typhlocybini.

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.

Complete mitochondrial genome sequence of Empoascanara dwalata (Hemiptera: Cicadellidae: Typhlocybinae)

In this study, the complete mitochondrial genome of the leafhopper species Empoascanara dwalata (Hemiptera: Cicadellidae: Typhlocybinae) was sequenced and annotated, and the sequence is available in the GenBank with accession number: MT350235. The circular mitogenome with 15,271 bp in length, including 13 protein-coding genes (PCGs), 22 tRNA genes, two rRNA genes, and one AT-rich region. The nucleotide composition is as follows: A (41.7%), T (34.7%), G (10.6%), and C (13.0%). 11 PCGs have ATN as the start codon, except for atp8 and nad5 genes have TTG. The conventional termination codons (TAA or TAG) occur in 11 PCGs, while cox2 and nad5 use incomplete codon (T––) as termination codon. Phylogenetic analysis using 13 PCGs showed that E. dwalata was clustered with three other Typhlocybine species, which was consistent with the conventional classification.

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.

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.

The complete mitochondrial genome of Empoascanara sipra (Hemiptera:Cicadellidae:Typhlocybinae) with phylogenetic consideration

The complete mitochondrial genome of the leafhopper Empoascanara sipra Dworakowska, 1980 was sequenced and annotated from this study. The circular molecule is 14,827 bp in length, including 13 protein-coding genes (PCGs), 22 tRNA genes, two rRNA genes, and one AT-rich region. The base composition of the genome is as follows: A (42.6%), T (33.9%), G (10.0%), and C (13.4%). 11 PCGs have ATN as the start codon, except for atp8 and nad5 genes have TTG. The conventional termination codons (TAA or TAG) occur in 11 PCGs, whereas cox2 and nad2 use incomplete codon (T--) as termination codon. Phylogenetic analysis using 13 PCGs showed that E. sipra was clustered with two Typhlocybinae species, which was consistent with the conventional classification.

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.

The mitochondrial genomes of three skippers: Insights into the evolution of the family Hesperiidae (Lepidoptera)

We sequenced the mitogenomes of Astictopterus jama, Isoteinon lamprospilus and Notocrypta curvifascia to obtain further insight into the mitogenomic architecture evolution and performed phylogenetic reconstruction using 29 Hesperiidae mitogenome sequences. The complete mitogenome sequences of A. jama, I. lamprospilus and N. curvifascia are 15,430, 15,430 and 15,546 bp in size, respectively. All contain 13 protein-coding genes, 2 ribosomal RNA genes, 22 transfer RNA genes, and an A + T-rich region. Nucleotide composition is A + T biased, and the majority of the protein-coding genes exhibit a negative AT-skew, which is reflected in the nucleotide composition, codon, and amino acid usage. The A + T-rich region is comprised of nonrepetitive sequences, including the motif ATAGA followed by a poly-T stretch, a microsatellite-like element next to the ATTTA motif, and a poly-A adjacent to tRNAs. Although most genes evolve under a strong purifying selection, the entire nad gene family (especially nad6) exhibits somewhat relaxed purifying selection, and atp8, evolving under a highly relaxed selection, is an outlier in the family Hesperiidae. Several different approaches relatively consistently indicated that nad6, atp8 and nad4 are comparatively fast-evolving genes in this family, which may have implications for future phylogenetic, population genetics and species diagnostics studies. For phylogenetic analyses of Hesperiidae, we tested a few datasets, and found that the one comprising all 37 genes produced the highest node support, indicating that the inclusion of RNAs improves the phylogenetic signal. Results indicate that subfamilies Euschemoninae, Heteropterinae, and Coeliadinae are monophyletic with strong nodal support, but Pyrginae and Eudaminae are paraphyletic. Finally, we confirm that A. jama and I. lamprospilus are close relatives.