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

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

The complete mitochondrial genome of Morishitium polonicum (Trematoda, Cyclocoelidae) and its phylogenetic implications

Trematodes can adversely impact the health and survival of wild animals. The trematode family Cyclocoelidae, which includes large digenean bird parasites, lacks molecular analysis, and reclassifications have not been supported. This study produced the first fully assembled and annotated mitochondrial genome sequence for the trematode Morishitium polonicum. The whole length of the M. polonicum (GenBank accession number: OP930879) mitogenome is 14083 bp, containing 22 transfer ribonucleic acids (tRNAs), 2 ribosomal RNAs (rRNAs, rrnL and rrnS), and a noncoding control section (D-loop) 13777 to 13854 bp in length. The 12 PCG areas have 3269 codons and a total length of 10053 bp, which makes up 71.38% of the mitochondrial genome's overall sequence. Most (10/12) of the PCGs that code for proteins begin with ATG, while the nad4L and nad1 genes have a GTG start codon. Phylogenetic analysis using the concatenated nucleotide sequences of 12 PCGs, and the ML tree analysis results showed that M. polonicum is more closely related to with Echinostomatidae and Fasciolidae, which indicates that the family Cyclocoelidae is more closely associated with Echinochasmidae. This study provides mtDNA information, and analysis of mitogenomic structure and evolution. Moreover, we aimed to understand the phylogenetic relationships of this fluke.

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.

De novo mitochondrial genome sequencing of Cladonia subulata and phylogenetic analysis with other dissimilar species

In this study, the complete mitochondrial genome of Cladonia subulata (L.) FH Wigg was sequenced and assembled and then compared with those of other Cladonia species. The mitogenome of Cladonia subulata, the type species of Cladonia, consisted of a circular DNA molecule of 58,895 bp 44 genes (15 protein-coding genes, 2 rRNA genes, and 27 tRNA genes). The base composition had shown an obvious AT preference, and all 27 tRNA genes formed a typical clover structure. Comparison with other 7 Cladonia species indicated that the duplication/loss of tRNAs had occurred during evolution, and introns appeared to explain the variation in cox1 genes in Cladonia, the mitochondrial genome tends to be generally conservative and local dynamic changes. Repeat sequences were mainly located in gene intervals, which were mainly distributed among intergenic spacers and may cause rearrangement of the mitogenome. The phylogenetic results showed that Cladonia subulata and C. polycarpoides were assigned to the Cladonia Subclade. The results add to the available mitochondrial genome sequence information of Cladonia subulata, provide basic data for the systematic development, resource protection, and genetic diversity research in Cladonia subulata, and also provide theoretical support for further genomic research of lichens.

Description of the whole mitochondrial genome of Bhatia longiradiata (Hemiptera, Cicadellidae, Deltocephalinae: Drabescini) and phylogenetic relationship

BACKGROUND

Mitochondrial genomes are extremely conserved in genetic processes and valuable molecular indications for phylogenetic and evolutionary examination, but the mitochondrial genome of Bhatia has not yet been reported.

OBJECTIVE

The target of this writing was to clarify the structural module of the mitochondrial genes of Bhatia longiradiata, verify the monophyletic of Drabescini, and explore the phylogenetic relationship between Drabescini with other leafhoppers.

METHODS

We performed sequencing and explanatory note of the mitochondrion of Bhatia longiradiata. The phylogeny relation was created by ML and Bayesian approaches using three dissimilar datasets (PCG12, PCG12rRNA, and AA), which were constructed to discuss the phylogenetic status of Bhatia longiradiata.

RESULTS

To report the architectural feature of the chondriosome of Bhatia longiradiata is a seal double-stranded annular molecule with 16,122 bp measurement and cover typically 37 genes. Several tandem repetitive units were observed in an AT enrichment area. The analysis showed that the branching relationships among the six trees were generally consistent, and each of the subfamilies was individually clustered into a monophyletic group within Cicadellidae. Bhatia longiradiata and other members of the Drabescini were aggregated into a clade that was situated within the Deltocephalinae.

CONCLUSION

The mitochondrial genome of Bhatia longiradiata covers 37 typical genes and a control region, which covers six tandem repeats. All species of Drabescini procedure a clade within Deltocephalinae. Drabescini and Scaphoideini form a branch and show a sister relationship with strong support. Therefore, we support the relegation of Selenocephalinae to a clan within Deltocephalinae.

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.

First description of the mitogenome and phylogeny:Aedes vexansand Ochlerotatus caspius of the Tribe Aedini (Diptera: Culicidae)

Culicidae, the mosquito family, includes more than 3600 species subdivided into the subfamilies Anophelinae and Culicinae. One-third of mosquitoes belong to the Aedini tribe, which is subordinate to the subfamily Culicinae, which comprises common vectors of viral zoonoses. The tribe of Aedini is extremely diverse in morphology and geographical distribution and has high ecological and medical significance. However, knowledge about the systematics of the Aedini tribe is still limited owing to its large population and the similar morphological characteristics of its species. This study provides the first description of the complete mitochondrial (mt) genome sequence of Aedes vexans and Ochlerotatus caspius belonging to the Aedini tribe. The mt genomes of A. vexans and O. caspius are circular molecules that are 15,861 bp and 15,954 bp in size, with AT contents of 78.54% and 79.36%, respectively. Both the circular mt genomes comprise 37 functional subunits, including 13 protein-coding genes (PCGs), two ribosomal RNA genes, 22 transfer RNA genes (tRNAs), and a control region (also known as the AT-rich region). The most common start codons are ATT/ATG, apart from cox1 (TCG) and nad5 (GTG), while TAA is the termination codon for all PCGs. All tRNAs have a typical clover leaf structure, except tRNA Ser1. Phylogenetic analysis of the concatenated, aligned amino acid sequences of the 13 PCGs showed that A. vexans gathered with Aedes sp. in a sister taxon, and O. caspius gathered with Ochlerotatus sp. in a sister taxon. The findings from the present study support the concept of monophyly of all groups, ratify the current taxonomic classification, and provide vital molecular marker resources for further studies of the taxonomy, population genetics, and systematics of the Aedini tribe.

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.

First Description of the Mitogenome and Phylogeny of Culicinae Species from the Amazon Region

The Culicidae family is distributed worldwide and comprises about 3587 species subdivided into the subfamilies Anophelinae and Culicinae. This is the first description of complete mitochondrial DNA sequences from Aedes fluviatilis, Aedeomyia squamipennis, Coquillettidia nigricans, Psorophora albipes, and Psorophora ferox. The mitogenomes showed an average length of 15,046 pb and 78.02% AT content, comprising 37 functional subunits (13 protein coding genes, 22 tRNAs, and two rRNAs). The most common start codons were ATT/ATG, and TAA was the stop codon for all PCGs. The tRNAs had the typical leaf clover structure, except tRNA^Ser1. Phylogeny was inferred by analyzing the 13 PCGs concatenated nucleotide sequences of 48 mitogenomes. Maximum likelihood and Bayesian inference analysis placed Ps. albipes and Ps. ferox in the Janthinosoma group, like the accepted classification of Psorophora genus. Ae. fluviatilis was placed in the Aedini tribe, but was revealed to be more related to the Haemagogus genus, a result that may have been hampered by the poor sampling of Aedes sequences. Cq. nigricans clustered with Cq. chrysonotum, both related to Mansonia. Ae. squamipennis was placed as the most external lineage of the Culicinae subfamily. The yielded topology supports the concept of monophyly of all groups and ratifies the current taxonomic classification.

Comparative analysis of twelve mitogenomes of Caliscelidae (Hemiptera: Fulgoromorpha) and their phylogenetic implications

Here, the complete mitochondrial genomes (mitogenomes) of 12 Caliscelidae species, Augilina tetraina, Augilina triaina, Symplana brevistrata, Symplana lii, Neosymplana vittatum, Pseudosymplanella nigrifasciata, Symplanella brevicephala, Symplanella unipuncta, Augilodes binghami, Cylindratus longicephalus, Caliscelis shandongensis, and Peltonotellus sp., were determined and comparatively analyzed. The genomes varied from 15,424 to 16,746 bp in size, comprising 37 mitochondrial genes and an A+T-rich region. The typical gene content and arrangement were similar to those of most Fulgoroidea species. The nucleotide compositions of the mitogenomes were biased toward A/T. All protein-coding genes (PCGs) started with a canonical ATN or GTG codon and ended with TAN or an incomplete stop codon, single T. Among 13 PCGs in 16 reported Caliscelidae mitogenomes, cox1 and atp8 showed the lowest and highest nucleotide diversity, respectively. All PCGs evolved under purifying selection, with atp8 considered a comparatively fast-evolving gene. Phylogenetic relationships were reconstructed based on 13 PCGs in 16 Caliscelidae species and five outgroups using maximum likelihood and Bayesian inference analyses. All species of Caliscelidae formed a steadily monophyletic group with high support. Peltonotellini was present at the basal position of the phylogenetic tree. Augilini was the sister group to Caliscelini and Peltonotellini.

A complete sequence of mitochondrial genome of Neolamarckia cadamba and its use for systematic analysis

Neolamarckia cadamba is an important tropical and subtropical tree for timber industry in southern China and is also a medicinal plant because of the secondary product cadambine. N. cadamba belongs to Rubiaceae family and its taxonomic relationships with other species are not fully evaluated based on genome sequences. Here, we report the complete sequences of mitochondrial genome of N. cadamba, which is 414,980 bp in length and successfully assembled in two genome circles (109,836 bp and 305,144 bp). The mtDNA harbors 83 genes in total, including 40 protein-coding genes (PCGs), 31 transfer RNA genes, 6 ribosomal RNA genes, and 6 other genes. The base composition of the whole genome is estimated as 27.26% for base A, 22.63% for C, 22.53% for G, and 27.56% for T, with the A + T content of 54.82% (54.45% in the small circle and 54.79% in the large circle). Repetitive sequences account for ~ 0.14% of the whole genome. A maximum likelihood (ML) tree based on DNA sequences of 24 PCGs supports that N. cadamba belongs to order Gentianales. A ML tree based on rps3 gene of 60 species in family Rubiaceae shows that N. cadamba is more related to Cephalanthus accidentalis and Hymenodictyon parvifolium and belongs to the Cinchonoideae subfamily. The result indicates that N. cadamba is genetically distant from the species and genera of Rubiaceae in systematic position. As the first sequence of mitochondrial genome of N. cadamba, it will provide a useful resource to investigate genetic variation and develop molecular markers for genetic breeding in the future.

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.

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.

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

Characterization of the Complete Mitochondrial Genomes from Two Nitidulid Pests with Phylogenetic Implications

The complete mitochondrial genomes of Xenostrongylusvariegatus and Epuraea sp. were sequenced and analyzed. The total genome lengths are 17,657 and 16,641 bp, with an A+T content of 77.2% and 76.4%, respectively. Each mitochondrial genome consists of 37 coding genes and a non-coding (AT-rich) region. All protein-coding genes (PCGs) start with the standard start codon, ATN, and end with complete stop codons, TAA and TAG, or an incomplete stop codon, T. All tRNAs can be folded into the typical clover-leaf secondary structure, with the exception of trnS1 in both species with a reduced dihydrouridine (DHU) arm. The AT-rich region has tandem repeats differing in both number and length. Genetic distance and Ka/Ks analyses show that nad6 has a higher variability and more rapid evolutionary rate than other PCGs. Both maximum likelihood and Bayesian inference phylogenetic analyses based on 13 PCGs and 2 ribosome DNAs (rDNAs) agree with the previous phylogenies in supporting the Nitidulidae monophyly and the sister-group relationship of Kateretidae + (Monotomidae + Nitidulidae).

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.