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

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.

The complete mitochondrial genome of leafhopper Atkinsoniella nigrita (Hemiptera: Cicadellidae) with the shortest 12S rRNA and longest tRNA-Lys of the Atkinsoniella genus

The complete mitochondrial genome (mitogenome) of the leafhopper Atkinsoniella nigrita Zhang & Kuoh, 1993 was determined in this study. The entire sequence was 16,011 base pairs (bp) in length. The new mitogenome consists of a typical set of genes: 13 protein-coding genes (PCGs), two ribosomal RNA (rRNA) genes, 22 transfer RNA (tRNA) genes, and one control region of 1720 bp in length. The base composition of the mitogenome was A = 41.7%, T = 38.2%, C = 10.7%, and G = 9.4%. This is the classical structure for most insect mitogenomes without any gene arrangement observed. Compared with other known mitochondrial genomes of 15 Atkinsoniella species, the new mitogenome with three PCGs (ND2, ND5, and ND4L) shared the same gene base length, start codon and stop codon, and the shortest 12S rRNA (729 bp) and the longest tRNA-Lys (73 bp) within the genus Atkinsoniella. A phylogenetic analysis of 31 Cicadellinae and two Ledrinae concatenated sequences of 13 PCGs of their mitogenomes using Bayesian inference (BI) revealed that A. nigrita belongs to the genus Atkinsoniella with strong nodal support (BI posterior probability = 1).

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

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

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

Characterization and phylogenetic implications of newly sequenced mitogenomes of Five Mileewa and Processina species from China (Hemiptera: Cicadellidae: Mileewinae)

To accumulate evidence of the phylogeny of Mileewinae and the relationships among Mileewa, Processina, and Ujna genera, we sequenced the complete mitochondrial genomes of four Mileewa spp., namely, Mileewa mira, Mileewa lamellata, Mileewa sharpa, and Mileewa amplimacula. The first complete mitogenome of the genus Processina (P. sexmaculata), established by Yang, Deitz & Li from China and comprising five species, was also sequenced in this study. Annotation showed that the five mitogenomes were 14787 -15436 bp in length, and all harbored 37 typical genes. The AT content of the five mitogenomes ranged from 78.3% to 80.2%, which was similar to that of other sequenced Mileewinae species. For protein-coding genes (PCGs), ATN was the start codon, while atp8 and nad5 genes were initiated with TTG, and a great majority of them used TAA or TAG as stop codons, whereas cox2 and nad1 ended with an incomplete codon T-. All tRNAs had a typical cloverleaf secondary structure, except for trnS1, which had a reduced dihydrouridine arm. We further used 59 Membracoidea species and two outgroups to reconstruct phylogenetic trees based on 13 PCGs under an independent partition model with Bayesian inference and Maximum-likelihood methods. Among these two trees, each of the subfamilies Cicadellinae, Typhlocybinae, and Mileewinae were recovered as a monophyletic group with high support values, suggesting that Typhlocybinae was more ancient than Mileewinae and Cicadellinae. Within the Mileewinae subfamily, all species maintained the same relationships and topologies according to both the BI and ML analyses (PP > 0.8, BS > 83) as follows: (M. sharpa + (U. puerana + ((M. ponta + (M. mira + M. lamellata)) + ((M. albovittata + (M. margheritae + M. amplimacula)) + (M. rufivena + (P. sexmaculata + M. alara)))))), and the monophyly of the genera Processina, Mileewa and Ujna were not supported. This study further enriches the Mileewinae mitogenome database and will contribute to future research on the systematics, evolution, and classification of this group.

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

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

Enhanced Resolution of Evolution and Phylogeny of the Moths Inferred from Nineteen Mitochondrial Genomes

The vast majority (approximately 90%) of Lepidoptera species belong to moths whose phylogeny has been widely discussed and highly controversial. For the further understanding of phylogenetic relationships of moths, nineteen nearly complete mitochondrial genomes (mitogenomes) of moths involved in six major lineages were sequenced and characterized. These mitogenomes ranged from 15,177 bp (Cyclidia fractifasciata) to 15,749 bp (Ophthalmitis albosignaria) in length, comprising of the core 37 mitochondrial genes (13 protein-coding genes (PCGs) + 22 tRNAs + two rRNAs) and an incomplete control region. The order and orientation of genes showed the same pattern and the gene order of trnM-trnI-trnQ showed a typical rearrangement of Lepidoptera compared with the ancestral order of trnI-trnQ-trnM. Among these 13 PCGs, ATP8 exhibited the fastest evolutionary rate, and Drepanidae showed the highest average evolutionary rate among six families involved in 66 species. The phylogenetic analyses based on the dataset of 13 PCGs suggested the relationship of (Notodontidae + (Noctuidae + Erebidae)) + (Geometridae + (Sphingidae + Drepanidae)), which suggested a slightly different pattern from previous studies. Most groups were well defined in the subfamily level except Erebidae, which was not fully consistent across bayesian and maximum likelihood methods. Several formerly unassigned tribes of Geometridae were suggested based on mitogenome sequences despite a not very strong support in partial nodes. The study of mitogenomes of these moths can provide fundamental information of mitogenome architecture, and the phylogenetic position of moths, and contributes to further phylogeographical studies and the biological control of pests.

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.

Comparative Analysis of Mitogenomes among Five Species of Filchnerella (Orthoptera: Acridoidea: Pamphagidae) and Their Phylogenetic and Taxonomic Implications

Mitogenomes have been widely used for exploring phylogenetic analysis and taxonomic diagnosis. In this study, the complete mitogenomes of five species of Filchnerella were sequenced, annotated and analyzed. Then, combined with other seven mitogenomes of Filchnerella and four of Pamphagidae, the phylogenetic relationships were reconstructed by maximum likelihood (ML) and Bayesian (BI) methods based on PCGs+rRNAs. The sizes of the five complete mitogenomes are Filchnerella sunanensis 15,656 bp, Filchnerella amplivertica 15,657 bp, Filchnerella nigritibia 15,661 bp, Filchnerella pamphagoides 15,661 bp and Filchnerella dingxiensis 15,666 bp. The nucleotide composition of mitogenomes is biased toward A+T. All tRNAs could be folded into the typical clover-leaf structure, except that tRNA Ser (AGN) lacked a dihydrouridine (DHU) arm. The phylogenetic relationships of Filchnerella species based on mitogenome data revealed a general pattern of wing evolution from long wing to increasingly shortened wing.