Next Article in Journal
Insecticidal Activity of Plectranthus amboinicus Essential Oil against the Stable Fly Stomoxys calcitrans (Diptera: Muscidae) and the Horse Fly Tabanus megalops (Diptera: Tabanidae)
Previous Article in Journal
Physical Methods for Electrical Trap-and-Kill Fly Traps Using Electrified Insulated Conductors
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

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

1
Institute of Entomology, Guizhou University, Guiyang 550025, China
2
Guizhou Provincial Key Laboratory for Agricultural Pest Management of the Mountainous Region, Guiyang 550025, China
3
College of Tobacco Sciences, Guizhou University, Guiyang 550025, China
*
Author to whom correspondence should be addressed.
Insects 2022, 13(3), 254; https://doi.org/10.3390/insects13030254
Submission received: 14 January 2022 / Revised: 25 February 2022 / Accepted: 1 March 2022 / Published: 3 March 2022
(This article belongs to the Section Insect Molecular Biology and Genomics)

Abstract

:

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.

1. Introduction

Cicadellinae is a relatively large subfamily within the family Cicadellidae (Hemiptera: Auchenorrhyncha: Cicadelloidea), which is the largest family in Hemiptera, with over 23,000 species in 25 subfamilies [1,2]. Cicadellinae is one of the largest and most brightly colored groups of leafhoppers, is distributed throughout all of the zoogeographical realms of the world, and is especially rich in species from tropical and subtropical regions. There are approximately 2400 species in approximately 330 genera distributed worldwide, and 259 species in 23 genera have been recorded in China now that Mileewa, Ujna, and Processina have been placed into the subfamily Mileewinae [3,4,5,6,7]. According to our experience collecting Cicadellinae samples in the wild, they prefer to live in the humid and well-lit forests of the mountains. Some Cicadellinae species are of considerable economic importance because they feed on the sap in the xylem of woody and herbaceous plants and are important vectors of the xylem-limited bacterial pathogen Xylella fastidiosa, which induces diseases of grapevines, citrus, alfalfa, almond, coffee, stone fruits, and landscape ornamentals [8,9,10,11].
Atkinsoniella, a large genus within the subfamily Cicadellinae, is mainly distributed in the Oriental realm and the Palearctic realm, with almost 98 valid species distributed worldwide and 88 species recorded in China [3,7]. Currently, the identification of Cicadellinae insects is mainly based on the external morphological and male genitalia characteristics of adults. However, some adult Cicadellinae species exhibit sexual dimorphism and polymorphism, showing differences in body color, markings, and individual size between male and female individuals, and there are both long-winged and short-winged types in the same sex of one species [3]. Atkinsoniella exhibits extremely similar morphological characteristics among some species, such as the 12 species that are included in this study. Their wings are all yellow with relatively few variable markings. It is very difficult to accurately identify the species without male genitalia, especially for the specimens that have been soaked in ethanol. As important morphological characteristics for identifying Cicadellinae species, the male genitalia of 12 Atkinsoniella species showed little difference, causing some yellow-winged Atkinsoniella species samples to be identified only to the genus level.
Therefore, the identification of yellow-winged species based on male genitalia is also confusing. Females can be identified as Atkinsoniella spp., making the identification of Atkinsoniella at the species level more challenging. Therefore, molecular methods are needed to validate and complement the morphological classification of some difficult-to-identify species. To date, only 41 molecular records of nine Atkinsoniella species (A. alternata Young, 1986; A. dactylia Yang et Li, 2000; A. grahami Young, 1986; A. heiyuana Li, 1992; A. opponens (Walker, 1851); A. sulphurate (Distant, 1908); A. thalia (Distant, 1918); A. trimaculata Li, 1992; A. xanthonota Kuoh, 1994) can be retrieved from the NCBI database, but most of them are molecular fragments (28S, COX1, COX2, H3, 16S, CYTB, SSU), and only three of them are complete mitochondrial genome (mitogenome) sequences: A. grahami, MW533712; A. xanthonota, MW533713, MT596736. Insect mitogenomes are usually double-stranded circular molecules with a typical length of 14.5–17.0 kb, encoding 13 protein-coding genes (PCGs), 22 transfer RNAs (tRNAs), two ribosomal RNAs (rRNAs), and a control region (CR) that contains the initial sites for replication and transcription [12,13,14]. The materially inherited mitogenome is an informative molecular marker for species identification [15,16], population genetics [17,18,19], evolutionary genomics, and phylogenetic analysis of various taxonomic levels [13,20,21,22,23,24,25] due to its highly conserved order and orientation, simple genetic structure, and relatively high evolution rate [12,13,14,26]. Compared with single or multiple genes for phylogenetic analysis, mitogenomes provide relatively rich genetic information, which can better infer phylogenetic relationships and the degrees of relatedness among taxonomic groups [27,28].
Until now, studies of the genus Atkinsoniella have mainly focused on the description of new species, and few have been conducted on mitochondrial genomes and phylogenetic relationships. To better understand the relationships within the genus Atkinsoniella, 12 mitogenomes of Atkinsoniella were sequenced and annotated, and their characteristics were comparatively analyzed. In addition, the comprehensive phylogenetic relationship within the subfamily Cicadellinae was determined based on mitochondrial datasets.

2. Materials and Methods

2.1. Genome Organization

The adult specimens of 12 Atkinsoniella species (A. aurantiaca Cai et Kuoh, 1995; A. curvata Zhang et Kuoh 1993; A. flavipenna Li et Wang, 1992; A. longiuscula Feng et Zhang, 2015; A. thalia (Distant, 1918); A. thaloidea Young, 1986; A. tiani Yang, Meng et Li, 2017; A. uniguttata Li, 1993; A. warpa Yang, Meng et Li, 2017; A. wui Yang, Meng et Li, 2017; A. xanthoabdomena Yang, Meng et Li, 2017; A. yunnanana Yang, Meng et Li, 2017) used in this study were collected in China, and the detailed collection information is listed in Table S1. All of the fresh specimens were preserved in absolute ethanol immediately and stored at −20 °C in the laboratory before DNA extraction. Twelve Atkinsoniella species were identified through morphological characteristics as described by Yang et al. [3]. Total DNA was prepared from the head and thorax muscle tissues of a single adult individual using a DNeasy® Tissue Kit (Qiagen, Hilden, Germany) according to the manufacturer’s protocol. The remainder of the specimens (wings and abdomens) were saved as vouchers and deposited at the Institute of Entomology, Guizhou University, Guiyang, China (GUGC).

2.2. PCR Amplification, Sequencing, and Assembly

The COX1 fragments as the starting references for assembling mitochondrial genomes were amplified using the universal primer pairs LCO1490 (5′-GGTCAACAAATCATAAAGATATTGG-3′) and LCO2198 (5′-TAAACTTCAGGGTGACCAAAAAATCA -3′) [29]. Primer sequences were synthesized by Sangon Biotech Company (Shanghai, China). PCRs were performed in 25 µL reaction volumes in total, which contained 9 µL ddH2O, 13 µL 2X SanTaq PCR Master Mix (Sangon Biotech Company), 1 µL forward primer (0.01 mol/L), 1 µL reverse primer (0.01 mol/L), and 1 µL DNA template. The PCR amplification conditions were as follows: 94 °C for 5 min, followed by 35 cycles of 95 °C for 1 min, 51 °C for 1 min, elongation at 72 °C for 2 min, and a final extension step at 72 °C for 10 min. The PCR products were sequenced by Sangon Biotech Company. Total genomic DNA was used for library preparation and sequenced using the Illumina NovaSeq6000 platform with 150 bp paired-end reads at Berry Genomics (Beijing, China). Approximately 6 Gb clean data of each species were obtained and assembled using Getorganelle 1.7.4.1 [30] or NOVOPlasty 2.7.2 [31] with the corresponding COX1 fragment of each species as a starting reference.

2.3. Mitogenome Annotation and Analysis

Initial annotations of these 12 mitogenomes were conducted with MitoZ 2.4-alpha [32] with the invertebrate mitochondrial genetic codes. The MITOS web server (http://mitos.bioinf.uni-leipzig.de/index.py, accessed on 1 July 2021) [33], tRNAscan-SE search server (http://lowelab.ucsc.edu/tRNAscan-SE/, accessed on 4 July 2021) [34] and open reading frames (ORFs) were employed to reconfirm the locations of tRNAs and PCGs with invertebrate mitochondrial genetic codes. Geneious Prime 2021.1.1 (https://www.geneious.com, accessed on 25 July 2021) was used to align the newly sequenced mitogenomes with the available mitogenomes of Cicadellidae in GenBank for further reconfirmation. The secondary structures of tRNA genes were predicted by the MITOS2 web server (http://mitos2.bioinf.uni-leipzig.de/index.py, accessed on 6 July 2021) and drawn with VARNAv3-92 [35], Rnaviz [36] and Adobe Illustrator SC6. Two rRNA genes were determined according to the locations of adjacent tRNA genes and compared with homologous mitogenome sequences of other Cicadellidae species using the MAFFT [37] algorithm in Geneious Prime 2021.1.1. Strand asymmetry was calculated using the following formulas: AT skew = [A − T]/[A + T] and GC skew = [G − C]/[G + C] [38]. The A + T content of nucleotide sequences, codon usage, and genetic distances (using the Kimura-2-parameter method) of PCGs were computed by MEGA 6.0 [39]. A heatmap of the relative synonymous codon usage (RSCU) values of PCGs (excluding stop codons) in the 12 newly sequenced mitogenomes was drawn using TBtools v1.09832 [40]. The rates of nonsynonymous substitutions (Ka) and synonymous substitutions (Ks) of 13 aligned PCGs were determined by DnaSP v6.12.03 [41]. The structure maps of mitogenomes were visualized using Geneious Prime 2021.1.1 and Adobe Illustrator SC6. These 12 newly sequenced mitogenomes were submitted to GenBank with accession numbers OL677863–OL677874.

2.4. Phylogenetic Analysis

In the phylogenetic analyses, the 12 newly obtained mitogenomes and 18 complete or nearly complete mitogenomes of Cicadellinae species available in GenBank (seven unannotated mitogenomes were annotated using the methods described in Section 2.3) were considered ingroup, and Taharana fasciana and Iassus dorsalis from Coelidiinae and Iassinae served as the outgroup (Table 1). Three datasets were concatenated for phylogenetic analysis: (1) cds_faa: amino acid sequences of the PCGs; (2) cds12_fna: first and second codon positions of the PCGs; (3) cds12_rrna: the first and the second codon positions of the PCGs and two rRNA genes. The tools used to obtain these datasets included SeqKit v.0.10.1 [42] (extract each gene from annotated mitogenomes), MAFFT v.7.3 [37] (align each gene), trimAl v1.4.1 [43] (remove gaps and ambiguous sites), and FASconCAT-G v.1.04 [44] (concatenate nucleotide and protein alignments). The partitioned datasets were analyzed using maximum-likelihood (ML) and Bayesian inferences (BI). ML analysis was reconstructed with IQ-TREE v.1.6.8 [45]. The partitioning schemes and the substitution models were estimated by ModelFinder [46] implemented in the IQ-TREE package with the Bayesian information criterion (BIC). Branch support was calculated with 10000 ultrafast bootstrap and 1000 SH-aLRT replicates. BI analysis was performed on MrBayes 3.2.6 [47] with the default settings by simulating four independent runs for 5–10 million generations with sampling every 1000 generations; after the average standard deviation of split frequencies fell below 0.01, the initial 25% of samples were discarded as burn-in. Partitioning schemes and the best-fitting substitution models were determined in PartitionFinder v.2.1.1 with BIC criterion and greedy-search algorithm [48].

3. Results

3.1. Genome Organization

The complete mitochondrial genomes of the 12 newly sequenced Atkinsoniella species were obtained and deposited in the GenBank database (Table S1). All of the mitogenomes are closed circular molecules ranging from 15,034 bp (A. thalia) to 15,988 bp (A. curvata) in length, encoding 37 typical mitochondrial genes, including 13 protein-coding genes (PCGs), 22 transfer RNA genes (tRNAs), two ribosomal RNA genes (rRNAs), and a control region (CR). A total of 23 genes (9 PCGs and 14 tRNAs) were encoded on the majority strand (N-strand), while the remaining 14 genes (4 PCGs, 8 tRNAs, and 2 rRNAs) were encoded on the minority strand (J-strand). All 12 mitogenomes showed identical gene arrangements and orientations (Figure 1, Table S2), which were consistent with other Cicadellidae species. The mitogenomes had a significant bias toward A and T bases, with A + T contents ranging from 77.7 (A. wui) to 79.8 (A. uniguttata). All the AT skews were positive (ranging from 0.066 to 0.104), and all the GC skews were negative (ranging from −0.141 to −0.055) (Table S3), indicating that the occurrence of As and Cs was slightly higher than that of Ts and Gs. Variable intergenic spacers and overlap were commonly found in the 12 mitogenomes; the longest intergenic spacer was 4 bp between trnY and COX1, and the largest overlapping region was 8 bp, where trnC was included in trnW (Table S2).

3.2. Protein-Coding Genes and Codon Usage

Among the 13 PCGs of the 12 Atkinsoniella mitogenomes, nine were encoded on the J-strand, while the other four were encoded on the N-strand. Almost all of the newly sequenced mitogenomes exhibited similar start and stop codons. All the PCGs were initiated with typical ATN codons (ATT, ATG, ATC, or ATA), except for ATP8 and ND5 which were initiated with a TTG codon, and ATP6 of three species (A. aurantiaca, A. curvata, and A. flavipenna) which started with GTG. Among the 13 PCGs of the 12 Atkinsoniella species, most PCGs terminated with complete termination codons TAA or TAG, whereas COX2, COX3, ND5, and ND6 were terminated with a single T or TA residue (Table S4). Such incomplete terminal codons are widely recognized in invertebrate mitogenomes and are assumed to be generated by the polyadenylation process [14]. Meanwhile, the stop codon TAA occurred more frequently than TAG. Among the 13 PCGs, ATP8 was the shortest (153 bp), and ND5 (1675 bp) was the longest in all the 12 Atkinsoniella mitogenomes. The 13 PCGs were all 10,958 bp in length, encoding 3642 amino acids. The scope of A + T content in 13 PCGs ranged from 76.4 (A. thalia, A. wui) to 78.4 (A. uniguttata), and the A + T content of the third codon positions was much higher than that of the first and second codon positions (PCG1: 72.4–74.1, PCG2: 69.1–69.5, PCG3: 87.0–91.7) (Table S3), which is typical for the family Cicadellidae mitogenomes [59,60,61]. All 13 PCGs showed a negative AT skew (from −0.148 to −0.132) and a positive GC skew (from 0.010 to 0.044). The number of codon usage and relative synonymous codon usage (RSCU) of PCGs of the 12 determined Atkinsoniella mitogenomes were calculated and are shown in Figure 2 and Table S5, which show that the six most frequently used codons were AUU(I), UUA(L), UUU(F), AUA(M), UAU(Y), and AAU(N), while CUG(L), AGG(S), ACG(T), CUC(L), GUC(V), and GCG(A) were rarely used. Furthermore, the overall analysis of RSCU indicated that the codons ending up with T or A were more frequently used than those ending up with C or G (Figure 2).
To evaluate the evolutionary patterns among 13 PCGs, the Ka (nonsynonymous substitutions), Ks (synonymous substitutions), and Ka/Ks (ω) values were estimated for each PCG of the 30 Cicadellinae species and two outgroups (Figure 3). The Ka and Ks values ranged from 0.08 (COX1) to 0.26 (ATP8) and 0.40 (ATP8) to 0.54 (COX1), respectively. The Ka/Ks values of all PCGs were less than one and increased from 0.15 for COX1 to 0.66 for ATP8, indicating that the 13 PCGs were under purifying selection. Furthermore, COX1, which has been widely used as a marker for species identification, showed the lowest p-distance (0.215). The result was consistent with the studies of Ledrinae and Typhlocybinae (Cicadellidae) [62,63], as well as Seirinae (Collembola: Entomobryidae) [64], but different from that of Cleridae (Coleoptera: Cleroidea) [65] (Figure 4). While ND2 possessed the greatest genetic distance with the value of 0.415 (Figure 4). The mitochondrial genes with the highest genetic divergence are variable among different insect groups [62,63,64,65] (Figure 4). In addition, the P-distances were measured based on the dataset of 13 PCGs across the 14 Atkinsoniella species to explore their sequence divergence (Table S6). According to previous studies, the intraspecific divergences are rarely greater than 2% and most are less than 1% [66],. This was also verified in studies of several insect groups, such as Lepidoptera [67], Piophilidae [68], Blattidae [69], Liposcelididae [70], etc. The recent study of leafhoppers based on COX1 gene fragments also showed that the intraspecific divergence was less than 1% [71]. As shown in Table S6 of this study, the minimum intergeneric distance was 2.6% between A. tiani and A. warpa, followed by 3.5% between A. flavipenna and A. longiuscula., and the greatest variation was 21.4% between A. wui and A. xanthonota, indicating that the 12 yellow-winged Atkinsoniella species were separate species with close relationships.

3.3. Transfer and Ribosomal RNA Genes

All 12 Atkinsoniella mitogenomes contained the typical 22 transfer RNAs (tRNAs), 14 of which were encoded by the J-strand, and eight were encoded by the N-strand. The size of 22 tRNAs of the 12 mitogenomes ranged from 61 bp (trnA and trnI in A. xanthoabdomena, and trnH in all 12 Atkinsoniella mitogenomes) to 72 bp (trnK in all newly sequenced Atkinsoniella mitogenomes except for that of A. uniguttata) (Table S2). The total length of 22 tRNAs ranged from 1426 to 1443 bp, and the A + T content ranged from 79.9 to 81.4. The AT skew (0.003–0.021) and GC skew (0.110–0.166) of all tRNAs were both positive. All tRNAs showed typical cloverleaf structures, except for trnS1, which lacked a dihydrouridine (DHU) arm and was instead replaced by a simple loop (Figure 5), which is commonly found in Cicadellidae insect mitogenomes [49,58,61,72,73,74,75,76]. All the anticodons of each tRNA were identical to those reported for Atkinsoniella species [49]. In the predicted secondary structure, the nucleotides on the anticodon loop of all tRNAs were highly conserved for 7 bp compared with the size-variable DHU and TΨC loops.
Two rRNA genes (s-rRNA and l-rRNA) were recognized in all the sequenced mitogenomes, and both were encoded on the N-strand. The l-rRNA was located between trnL1 and trnV, with lengths ranging from 1188 bp (A. warpa) to 1199 bp (A. aurantiaca, A. longiuscula). The s-rRNA was located between trnV and the control region, with lengths ranging from 730 bp (A. flavipenna, A. longiuscula, A. uniguttata, A. xanthobdomena) to 734 bp (A. aurantiaca). The A + T content of rRNAs ranged from 80.7% to 81.4%. Additionally, the two rRNAs showed a negative AT skew (−0.227–−0.142) and a positive GC skew (0.219–0.258).

3.4. Control Region

The control region, which is the longest noncoding region, contains controlling elements for replication and transcription. The control regions of the 12 Atkinsoniella mitogenomes were located between s-rRNA and trnI and were variable in length, ranging from 744 bp (A. thalia) to 1688 bp (A. curvata), with A + T contents ranging from 82.1% (A. yunnanana) to 88.6% (A. tiani). The AT skews of A. aurantiaca, A. flavipenna, and A. longiuscula were negative, with values of −0.02, −0.016, and −0.054, respectively. The remaining species shared a positive AT skew, ranging from 0.007 to 0.144. The GC skews were all negative, except for A. aurantiaca, which shared a positive value (0.046).

3.5. Phylogenetic Analyses

Phylogenetic relationships among 30 species of the subfamily Cicadellinae and two outgroup taxa from Coelidiinae and Iassinae were reconstructed based on 13 PCGs and two rRNAs using ML and BI analyses under the best partitioning scheme and models selected by Modelfinder (Table S7) and PartitionFinder (Table S8). Among the six resulting trees, the topology was highly identical within Cicadellinae and showed clearly separated branches of each species, with most nodes being highly supported in both ML and BI analyses, providing new insights into the phylogenetic relationship among Atkinsoniella, Bothrogonia, Paracrocampsa, Draeculacephala, Cicadella, Cuerna, Cofana, Homalodisca, and Kolla (Figure 6 and Figures S1–S5). All the phylogenetic analyses in this study strongly support that the genus Atkinsoniella was consistently recovered as a monophyletic group with strong support values and can be subdivided into two groups. The first group contained the 12 newly sequenced yellow-winged Atkinsoniella species. The second group was composed of sister taxa A. grahami and A. xanthonota [49]. The topological structure of the 12 newly sequenced species with highly similar morphological characteristics was completely congruent in all six phylogenetic trees, except for the differences observed in branches of A. yunnanana, A. uniguttata, and A. xanthoabdomena. The phylogenetic relationships within clade 1 and clade 2 were consistently recovered as (A. wui +((A. tiani + A. warpa) + (A. thalia + A. thaloidea))) and ((A. aurantiaca + A. curvata) + (A. flavipenna + A. longiuscula)) in all analyses, respectively (Figure 6 and Figures S1–S5). In the cds_faa_BI, cds12_fna_ML, and cds12_rrna_ML analyses, A. yunnanana and clade 1 clustered into one branch, forming a sister group with clade 2 (Figure 6 and Figures S1–S5). In the analyses of cds_faa_ML, cds12_fna_BI, and cds12_rrna_BI, A. yunnanana and clade 2 clustered into one branch, emerging as a sister group with clade 1. Additionally, A. xanthoabdomena was placed as the sister group to the remaining 11 yellow-winged species within all phylogenetic analyses, except that in the cds12_fna_BI analysis, the locations of A. xanthoabdomena and A. uniguttata were exchanged in cds12_fna_BI (Figure 6 and Figures S1–S5).
Traditionally, among the 12 yellow-winged Atkinsoniella species, some are difficult to identify due to their morphological similarities. For example, A. thalia and A. thaloidea are very similar in body color, markings, and even male genitalia; they can only be distinguished by the characteristics of their pygofer process. However, in this study, A. thalia and A. thaloidea were completely separated into two branches and strongly supported, demonstrating that the morphological characteristics that we currently use to identify Atkinsoniella species were proven to be reliable at the organelle-genome level. Furthermore, it also indicated that mitogenome sequences could be one of the effective methods to identify species with similar morphological characteristics and may contribute to the male and female correspondence of sexually dimorphic species as a means to shed light on the studies of insect biodiversity.
Meanwhile, the ML and BI analyses from this study also indicated that the genus Bothrogonia was coherently a monophyletic group, and the relationships within the genus were (((B. ferruginea + B. qiongana) +B. tongmaiana) + B. yunana), forming a sister group with the clade of Paracrocampsa sp. (MK251085) and Draeculacephala sp. (MK251110), which was completely concurrent with a previous study [52], except for the analysis based on the cds12_rrna dataset where Draeculacephala sp. (MK251110) was not included due to the lack of rRNA gene sequences. Meanwhile, Cicadella viridis emerged as a sister group with Cuerna septentrionalis and Cofana yasumatsui in all analyses. Unfortunately, the status of some species in the phylogenetic relationships remained questionable. For example, species from the genera Cofana, Draeculacephala, and Homalodisca were not clustered into a clade (e.g., Cofana yasumatsui and C. unimaculata, Draeculacephala sp. and D. crassicornis, Homalodisca sp. and H. vitripennis), showing that the three genera were polyphyletic groups. Furthermore, the Old World taxa (Atkinsoniella, Bothrogonia, Cicadella, Cofana, and Kolla) consistently formed well-supported clades with the New World taxa (Cuerna, Draeculacephala, Homalodisca, and Paracrocampsa) in all the phylogenetic trees. This was also observed in previous studies [51,52,54,55,77,78]. It is worth to notice that the genera included in the present study were all well-defined by previous studies. However, the species of the same genus from the New World in previous studies were not grouped together, strongly suggesting that some or all of these specimens may be misidentified. Thus, molecular data from additional individuals of each species is much needed in order to address these questions. Furthermore, additional specimens and molecular data are also needed to further clarify the phylogeny at the genus and species levels within the subfamily Cicadellinae.

3.6. Morphology

Although the 12 yellow-winged Atkinsoniella species are very similar in morphological characteristics, they can be distinguished by morphological characteristics. They share common morphological characteristics, with colour golden, pale yellow, pale white, or brown forewings without distinct spots or stripes, while there are some physical characteristics shown in Figure 7 that can be used to identify them in some degree. A. aurantiaca, A. flavipenna, A. longiuscula, A. thalia, A. thaloidea, A. tiani, A. warpa, A. wui, and A. yunnanana have a median black spot at the apical portion of the crown, and the median spot of basal portion is not V-shaped; A. curvata has no black spot at the apical portion of crown, and the spot of the basal portion is V-shaped; A. uniguttata has a median black spot at the apical portion of crown, but its basal portion does not have a spot; A. xanthoabdomena has a minimal black spot at the apical portion of crown in some individuals, and the basal portion without a spot; A. flavipenna and A. thalia have basal side spots on the scutellum. Moreover, their paraphysis and aedeagus in lateral view were shown in Figure 8 for further identification. The common characteristics of Atkinsoniella show apically forked, dorsally curved paraphysis, articulating with aedeagus. The major differences among the 12 species are reflected in the aedeagal size (A. flavipenna and A. longiuscula with aedeagus shorter and broader than others), and the different shape of the aedeagus and paraphysis. The other characteristics that can be used to identify them is provided in the following key:
Key to the 12 Atkinsoniella species
1.Basal portion of crown without median black spots ........................................................... 2
Basal portion of crown with a median black spot ................................................................ 3
2.Male pygofer process without branch ............................................................... A. uniguttata
Male pygofer process with a small branch at subapex .......................... A. xanthoabdomena
3.Apical portion of crown without black spots, the spot of basal portion V-shaped ............
..................................................................................................................................... A. curvata
Apical portion of crown with a median black spot, the spot of basal portion not V-shap
ed ............................................................................................................................................... 4
4.Crown with large basal black spot, distinctly larger than the anterior median one .........
............................................................................................................................................ A. wui
Crown with small basal black spot, as big as or smaller than the anterior median one
.................................................................................................................................................... 5
5.Scutellum with a black spot in each basal angle .................................................................. 6
Scutellum without black spots ............................................................................................... 7
6Male pygofer process with apical 1/3 constricted and contorted; aedeagus broad and
short ........................................................................................................................ A. flavipenna
Male pygofer process normal, with toothed process subapically; aedeagus stout and long,
apex curved dorsad ..................................................................................................... A. thalia
7.Male pygofer process with branch or toothed process subapically .................................. 8
Male pygofer process without branch or process ................................................................ 9
8.Male pygofer process with membranous branch; paraphysis curved dorsad from median
portion, and with apical 1/3 curved posteroventrally ........................... A. yunnanana
Male pygofer process with toothed process subapically, paraphysis curved not as above
................................................................................................................................... A. thaloidea
9.Male pygofer process foliated and contorted medially .................................................... 10
Male pygofer process banding and not contorted medially ............................................ 11
10.Forewing ivory, base orange-red; male pygofer rounded apically ............ A. longiuscula
Forewing orange except apical membrane; male pygofer truncated apically ...................
................................................................................................................................ A. aurantiaca
11.Male pygofer process not extending to apex of pygofer; aedeagus with transverse crack
at ventral margin medially .......................................................................................... A. tiani
Male pygofer process extending beyond apex of pygofer; aedeagus without transverse
crack at ventral margin medially .............................................................................. A. warpa

4. Discussion

In this study, complete mitochondrial genomes of 12 yellow-winged Atkinsoniella species with similar morphological characteristics (A. aurantiaca, A. curvata, A. flavipenna, A. longiuscula, A. thalia, A. thaloidea, A. tiani, A. uniguttata, A. warpa, A. wui, A. xanthoabdomena, A. yunnanana) were newly sequenced and analyzed. The sizes of the 12 mitogenomes varied from 15,034 (A. thalia) to 15,988 bp (A. curvata), which was mainly due to the variation in the lengths of control regions. The comparative analysis showed that all 12 mitogenomes had similar structural characteristics and nucleotide compositions, and the gene order was consistent with hypothetical ancestral insects. All the Atkinsoniella mitogenomes contained 37 typical mitochondrial genes (13 PCGs, 22 tRNAs, two rRNAs) and a CR, and had a significant bias toward A and T bases. All PCGs were initiated with ATN codons, except for ATP8 and ND5 of the 12 Atkinsoniella species which were initiated with a TTG codon, and ATP6 of A. aurantiaca, A. curvata, and A. flavipenna which were initiated with GTG. Meanwhile, most PCGs were terminated with TAA or TAG, whereas COX2, COX3, ND5, and ND6 were terminated with a single T or TA residue.
Both ML and BI analyses based on the concatenated datasets of cds_faa, cds12_fna, and cds12_rrna indicated that the genus Atkinsoniella and Bothrogonia were consistently recovered as monophyletic groups. Additionally, Atkinsoniella can be subdivided into two groups. The first group contained the 12 newly sequenced yellow-winged Atkinsoniella species, which were clearly separated, with most nodes receiving high support values among the six phylogenetic trees. The second group was composed of A. grahami and A. xanthonota, whose forewings were black or dark brown with blood-red, orange-red or yellow-brown longitudinal stripes. The analyses also produced a well-resolved framework for the relationships within the subfamily Cicadellinae. The phylogenetic relationships with the genus Atkinsoniella indicated that mitogenomics is an effective method for identifying closely related species and may help to understand their phylogenetic and evolutionary relationships. However, some generic relationships within Cicadellinae remain uncertain, possibly due to the misidentification of the New World species and limited specimens of each genus. Further studies with more taxon samples from various animal geographic regions, richer molecular data, and a combined approach of mitochondrial and nuclear markers are needed to elucidate the phylogenetic and evolutionary relationships within Cicadellinae.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/insects13030254/s1, Figure S1: Phylogenetic tree inferred by maximum likelihood (ML) based on the amino acids (cds_faa); Figure S2: Phylogenetic tree inferred by maximum likelihood (ML) based on the 1st and 2nd code position of protein-coding genes (cds12_fna); Figure S3: Phylogenetic tree inferred by Bayesian inference (BI) based on the 1st and 2nd code position of protein-coding genes (cds12_fna); Figure S4: Phylogenetic tree inferred by maximum likelihood (ML) based on the 1st and 2nd code position of protein coding genes concatenated rRNA (cds12_rrna); Figure S5: Phylogenetic tree inferred by Bayesian inference (BI) based on the 1st and 2nd code position of protein-coding genes concatenated rRNA (cds12_rrna), Table S1: Taxonomic information of newly sequenced species used for phylogenetic analysis in this study; Table S2: Annotations for the 12 newly sequenced Atkinsoniella mitogenomes; Table S3: Nucleotide composition of the 12 newly sequenced Atkinsoniella mitogenomes; Table S4: The start codons and stop codons of each protein-coding gene in the 12 Atkinsoniella mitogenomes; Table S5: The amount of codon usage of PCGs in the 12 newly sequenced Atkinsoniella mitogenomes; Table S6: Genetic distances of the 14 Atkinsoniella species in this study; Table S7: Best models were calculated by Modelfinder of cds_faa, cds12_fna, cds12_rrna datasets used in analysis; Table S8: Best models were calculated by PartitionFinder2 of cds_faa, cds12_fna, cds12_rrna datasets used in analysis.

Author Contributions

Conceptualization, M.-F.Y. and Y.J.; species identification, M.-F.Y. and Y.J.; methodology and experiments, Y.J., H.-X.L. and X.-F.Y.; data analysis, Y.J. and H.-X.L.; writing—original draft preparation, Y.J.; writing—review and editing, Y.J., H.-X.L., X.-F.Y., and M.-F.Y.; visualization, Y.J. and H.-X.L.; supervision, M.-F.Y.; project administration, M.-F.Y.; funding acquisition, M.-F.Y. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the Program of Excellent Innovation Talents, Guizhou Province, China ([2016]-4022); the National Natural Science Foundation of China (31672335); the International Cooperation Base for Insect Evolutionary Biology and Pest Control ([2016]-5802); the Guizhou Provincial Department of Education Youth Science and Technology Talents Growth Project (Qian Jiao He KY word [2018]103); the Introducing Talent Research Program for Guizhou University (Gui Da Ren Ji He word [2016]70); the National Nature Science Foundation of Guizhou University Cultivating program (Qian Ke He Platform Talent [2017]5788-32); the Graduate Innovation Fund, Guizhou Province, China (Qianjiaohe YJSCXJH(2019)023).

Institutional Review Board Statement

Not applicable.

Data Availability Statement

Date available on request.

Acknowledgments

The authors would like to thank Jian-Yue Qiu (Southwest University, Chongqing, China), Zai-Hua Yang (Guizhou Academy of Forestry, Guiyang, China), Jia-Jia Wang (GUGC), Xian-Yi Wang (GUGC), Xiao-Li Xu (GUGC), and Lu Yang (GUGC) for collecting specimens.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Dietrich, C.H. Phylogeny of the leafhopper subfamily Evacanthinae with a review of Neotropical species and notes on related groups (Hemiptera: Membracoidea: Cicadellidae). Syst. Entomol. 2004, 29, 455–487. [Google Scholar] [CrossRef]
  2. Dietrich, C.H. Keys to the families of Cicadomorpha and subfamilies and tribes of Cicadellidae (Hemiptera: Auchenorrhyncha). Fla. Entomol. 2005, 88, 502–517. [Google Scholar] [CrossRef]
  3. Yang, M.F.; Meng, Z.H.; Li, Z.Z. Hemiptera: Cicadellidae (II): Cicadellinae. Fauna Sinica: Insecta. Vol. 67; Han, X., Jiao, T., Eds.; Science Press: Beijing, China, 2017; ISBN 978-7-03-052950-3. (In Chinese) [Google Scholar]
  4. Feng, L.; Zhang, Y.L. Two new species in the genus Kolla Distant (Hemiptera: Cicadellidae: Cicadellinae) from China, with DNA barcoding data. Zootaxa 2017, 4250, 191–197. [Google Scholar] [CrossRef] [PubMed]
  5. Linnavuori, R.; Delong, D.M. Studies of the Neotropical Mileewaninae (Homoptera: Cicadellidae). J. Kansas Entomol. Soc. 1977, 50, 410–421. [Google Scholar]
  6. Young, D.A. Taxonomic Study of the Cicadellinae (Homoptera: Cicadellidae) Pt. 1: Proconiini; United States National Museum Bulletin 261; Smithsonian Institution Press: Washington, DC, USA, 1968; pp. 1–287. [Google Scholar] [CrossRef]
  7. Naveed, H.; Zhang, Y.L. Newly recorded leafhoppers of the subfamily Cicadellinae (Hemiptera: Cicadellidae) with description of a new species from Pakistan. Zootaxa 2018, 4504, 285–295. [Google Scholar] [CrossRef] [PubMed]
  8. Redak, R.A.; Purcell, A.H.; Lopes, R.S.; Blua, M.J.; Iii, R.F.M.; Andersen, P.C. The Biology of Xylem Fluid-Feeding Insect Vectors of Xylella Fastidiosa and Their Relation to Disease Epidemiology. Annu. Rev. Entomol. 2004, 49, 243–270. [Google Scholar] [CrossRef]
  9. Kleina, H.T.; Kudlawiec, K.; Esteves, M.B.; Dalbó, M.A. Settling and feeding behavior of sharpshooter vectors of Xylella fastidiosa on plum genotypes resistant to leaf scald disease. Eur. J. Plant Pathol. 2020, 158, 633–644. [Google Scholar] [CrossRef]
  10. Hopkins, D.L.; Purcell, A.H. Xylella fastidiosa: Cause of Pierce’s Disease of Grapevine and Other Emergent Diseases. Plant Dis. 2002, 86, 1056–1066. [Google Scholar] [CrossRef] [Green Version]
  11. Krugner, R.; Sisterson, M.S.; Backus, E.A.; Burbank, L.P.; Redak, R.A. Sharpshooters: A review of what moves Xylella fastidiosa. Austral Entomol. 2019, 58, 248–267. [Google Scholar] [CrossRef] [Green Version]
  12. Boore, J.L. Animal mitochondrial genomes. Nucleic Acids Res. 1999, 27, 1767–1780. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  13. Cameron, S.L. Insect mitochondrial genomics: Implications for evolution and phylogeny. Annu. Rev. Entomol. 2014, 59, 95–117. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  14. Wolstenholme, D.R. Animal Mitochondrial DNA: Structure and Evolution. Int. Rev. Cytol. 1992, 141, 173–216. [Google Scholar] [CrossRef] [PubMed]
  15. Yeo, D.; Puniamoorthy, J.; Ngiam, R.W.J.; Meier, R. Towards holomorphology in entomology: Rapid and cost-effective adult–larva matching using NGS barcodes. Syst. Entomol. 2018, 43, 678–691. [Google Scholar] [CrossRef]
  16. Sun, C.H.; Liu, H.Y.; Lu, C.H. Five new mitogenomes of Phylloscopus (Passeriformes, Phylloscopidae): Sequence, structure, and phylogenetic analyses. Int. J. Biol. Macromol. 2020, 146, 638–647. [Google Scholar] [CrossRef] [PubMed]
  17. Battaglia, V.; Gabrieli, P.; Brandini, S.; Capodiferro, M.R.; Javier, P.A.; Chen, X.G.; Achilli, A.; Semino, O.; Gomulski, L.M.; Malacrida, A.R.; et al. The worldwide spread of the tiger mosquito as revealed by mitogenome haplogroup diversity. Front. Genet. 2016, 7, 208. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  18. Du, Z.Y.; Hasegawa, H.; Cooley, J.R.; Simon, C.; Yoshimura, J.; Cai, W.; Sota, T.; Li, H. Mitochondrial genomics reveals shared phylogeographic patterns and demographic history among three periodical cicada species groups. Mol. Biol. Evol. 2019, 36, 1187–1200. [Google Scholar] [CrossRef] [PubMed]
  19. Ma, C.; Yang, P.C.; Jiang, F.; Chapuis, M.P.; Shali, Y.; Sword, G.A.; Kang, L. Mitochondrial genomes reveal the global phylogeography and dispersal routes of the migratory locust. Mol. Ecol. 2012, 21, 4344–4358. [Google Scholar] [CrossRef]
  20. Li, H.; Leavengood, J.M.; Chapman, E.G.; Burkhardt, D.; Song, F.; Jiang, P.; Liu, J.P.; Zhou, X.G.; Cai, W.Z. Mitochondrial phylogenomics of Hemiptera reveals adaptive innovations driving the diversification of true bugs. Proc. R. Soc. B Biol. Sci. 2017, 284, 20171223. [Google Scholar] [CrossRef] [PubMed]
  21. Song, F.; Li, H.; Liu, G.H.; Wang, W.; James, P.; Colwell, D.D.; Tran, A.; Gong, S.Y.; Cai, W.Z.; Shao, R.F. Mitochondrial genome fragmentation unites the parasitic lice of Eutherian mammals. Syst. Biol. 2019, 68, 430–440. [Google Scholar] [CrossRef] [Green Version]
  22. Li, H.; Shao, R.F.; Song, N.; Song, F.; Jiang, P.; Li, Z.H.; Cai, W.Z. Higher-level phylogeny of paraneopteran insects inferred from mitochondrial genome sequences. Sci. Rep. 2014, 5, 8527. [Google Scholar] [CrossRef] [Green Version]
  23. Liu, Y.Q.; Li, H.; Song, F.; Zhao, Y.S.; Wilson, J.J.; Cai, W.Z. Higher-level phylogeny and evolutionary history of Pentatomomorpha (Hemiptera: Heteroptera) inferred from mitochondrial genome sequences. Syst. Entomol. 2019, 44, 810–819. [Google Scholar] [CrossRef]
  24. Chen, L.; Wahlberg, N.; Liao, C.Q.; Wang, C.B.; Ma, F.Z.; Huang, G.H. Fourteen complete mitochondrial genomes of butterflies from the genus Lethe (Lepidoptera, Nymphalidae, Satyrinae) with mitogenome-based phylogenetic analysis. Genomics 2020, 112, 4435–4441. [Google Scholar] [CrossRef] [PubMed]
  25. Ren, L.P.; Zhang, X.Y.; Li, Y.; Shang, Y.J.; Chen, S.; Wang, S.W.; Qu, Y.H.; Cai, J.F.; Guo, Y.D. Comparative analysis of mitochondrial genomes among the subfamily Sarcophaginae (Diptera: Sarcophagidae) and phylogenetic implications. Int. J. Biol. Macromol. 2020, 161, 214–222. [Google Scholar] [CrossRef]
  26. Curole, J.P.; Kocher, T.D. Mitogenomics: Digging deeper with complete mitochondrial genomes. Trends Ecol. Evol. 1999, 14, 394–398. [Google Scholar] [CrossRef]
  27. Qin, J.; Zhang, Y.Z.; Zhou, X.; Kong, X.B.; Wei, S.J.; Ward, R.D.; Zhang, A.B. Mitochondrial phylogenomics and genetic relationships of closely related pine moth (Lasiocampidae: Dendrolimus) species in China, using whole mitochondrial genomes. BMC Genom. 2015, 16, 428. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  28. Na, L.; Gao, J.F.; Liu, G.H.; Fu, X.; Su, X.; Yue, D.M.; Gao, Y.; Zhang, Y.; Wang, C.R. The complete mitochondrial genome of Metorchis orientalis (Trematoda: Opisthorchiidae): Comparison with other closely related species and phylogenetic implications. Infect. Genet. Evol. 2016, 39, 45–50. [Google Scholar] [CrossRef] [PubMed]
  29. Folmer, O.; Black, M.; Hoeh, W.; Lutz, R.; Vrijenhoek, R. DNA primers for amplification of mitochondrial cytochrome c oxidase subunit I from diverse metazoan invertebrates. Mol. Mar. Biol. Biotechnol. 1994, 3, 294–299. [Google Scholar] [CrossRef]
  30. Jin, J.J.; Yu, W.B.; Yang, J.B.; Song, Y.; DePamphilis, C.W.; Yi, T.S.; Li, D.Z. GetOrganelle: A fast and versatile toolkit for accurate de novo assembly of organelle genomes. Genome Biol. 2020, 21, 241. [Google Scholar] [CrossRef]
  31. Dierckxsens, N.; Mardulyn, P.; Smits, G. NOVOPlasty: De novo assembly of organelle genomes from whole genome data. Nucleic Acids Res. 2017, 45, e18. [Google Scholar] [CrossRef] [Green Version]
  32. Meng, G.L.; Li, Y.Y.; Yang, C.T.; Liu, S.L. MitoZ: A toolkit for animal mitochondrial genome assembly, annotation and visualization. Nucleic Acids Res. 2019, 47, e63. [Google Scholar] [CrossRef] [PubMed]
  33. Bernt, M.; Donath, A.; Jühling, F.; Externbrink, F.; Florentz, C.; Fritzsch, G.; Pütz, J.; Middendorf, M.; Stadler, P.F. MITOS: Improved de novo metazoan mitochondrial genome annotation. Mol. Phylogenet. Evol. 2013, 69, 313–319. [Google Scholar] [CrossRef] [PubMed]
  34. Lowe, T.M.; Chan, P.P. tRNAscan-SE On-line: Integrating search and context for analysis of transfer RNA genes. Nucleic Acids Res. 2016, 44, W54–W57. [Google Scholar] [CrossRef] [PubMed]
  35. Darty, K.; Denise, A.; Ponty, Y. VARNA: Interactive drawing and editing of the RNA secondary structure. Bioinformatics 2009, 25, 1974–1975. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  36. De Rijk, P.; Wuyts, J.; De Wachter, R. Rna Viz 2: An improved representation of RNA secondary structure. Bioinformatics 2003, 19, 299–300. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  37. Katoh, K.; Standley, D.M. MAFFT multiple sequence alignment software version 7: Improvements in performance and usability. Mol. Biol. Evol. 2013, 30, 772–780. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  38. Perna, N.T.; Kocher, T.D. Patterns of nucleotide composition at fourfold degenerate sites of animal mitochondrial genomes. J. Mol. Evol. 1995, 41, 353–358. [Google Scholar] [CrossRef] [PubMed]
  39. Tamura, K.; Stecher, G.; Peterson, D.; Filipski, A.; Kumar, S. MEGA6 : Molecular Evolutionary Genetics Analysis Version 6.0. Mol. Biol. Evol. 2013, 30, 2725–2729. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  40. Chen, C.J.; Chen, H.; Zhang, Y.; Thomas, H.R.; Frank, M.H.; He, Y.H.; Xia, R. TBtools: An Integrative Toolkit Developed for Interactive Analyses of Big Biological Data. Mol. Plant 2020, 13, 1194–1202. [Google Scholar] [CrossRef] [PubMed]
  41. Librado, P.; Rozas, J. DnaSP v5: A software for comprehensive analysis of DNA polymorphism data. Bioinformatics 2009, 25, 1451–1452. [Google Scholar] [CrossRef] [Green Version]
  42. Shen, W.; Le, S.; Li, Y.; Hu, F.Q. SeqKit: A cross-platform and ultrafast toolkit for FASTA/Q file manipulation. PLoS ONE 2016, 11, e0163962. [Google Scholar] [CrossRef]
  43. Capella-Gutiérrez, S.; Silla-Martínez, J.M.; Gabaldón, T. trimAl: A tool for automated alignment trimming in large-scale phylogenetic analyses. Bioinformatics 2009, 25, 1972–1973. [Google Scholar] [CrossRef] [PubMed]
  44. Kück, P.; Longo, G.C. FASconCAT-G: Extensive functions for multiple sequence alignment preparations concerning phylogenetic studies. Front. Zool. 2014, 11, 81. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  45. Nguyen, L.T.; Schmidt, H.A.; Von Haeseler, A.; Minh, B.Q. IQ-TREE: A fast and effective stochastic algorithm for estimating maximum-likelihood phylogenies. Mol. Biol. Evol. 2015, 32, 268–274. [Google Scholar] [CrossRef] [PubMed]
  46. Kalyaanamoorthy, S.; Minh, B.Q.; Wong, T.K.F.; von Haeseler, A.; Jermiin, L.S. ModelFinder: Fast model selection for accurate phylogenetic estimates. Nat. Methods 2017, 14, 587–589. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  47. Ronquist, F.; Huelsenbeck, J.P. MrBayes 3: Bayesian phylogenetic inference under mixed models. Bioinformatics 2003, 19, 1572–1574. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  48. Lanfear, R.; Frandsen, P.B.; Wright, A.M.; Senfeld, T.; Calcott, B. Partitionfinder 2: New methods for selecting partitioned models of evolution for molecular and morphological phylogenetic analyses. Mol. Biol. Evol. 2017, 34, 772–773. [Google Scholar] [CrossRef] [Green Version]
  49. Jiang, Y.; Li, H.X.; Yu, X.F.; Yang, M.F. Characterization of Two Complete Mitochondrial Genomes of Atkinsoniella (Hemiptera: Cicadellidae: Cicadellinae) and the Phylogenetic Implications. Insects 2021, 12, 338. [Google Scholar] [CrossRef]
  50. Yu, P.F.; Li, Q.; Wang, M.X.; Cheng, P.; Lin, C.; Han, B.Y. Analysis of Complete Mitochondrial Genome and Phylogenetic Relationship of Bothrogonia ferruginea. J. Agric. Biotechnol. 2019, 27, 1246–1258. (In Chinese) [Google Scholar] [CrossRef]
  51. Xu, X.L.; Yan, B.; Yu, X.F.; Yang, M.F. The complete mitochondrial genome of Bothrogonia qiongana (Hemiptera: Cicadellidae) with phylogenetic analyses. Mitochondrial DNA Part B Resour. 2020, 5, 2750–2751. [Google Scholar] [CrossRef]
  52. Xu, X.L.; Yan, B.; Yu, X.F.; Yang, M.F. The complete mitochondrial genomes of Bothrogonia tongmaiana and Bothrogonia yunana (Hemiptera: Cicadellidae) from China. Mitochondrial DNA Part B Resour. 2021, 6, 613–614. [Google Scholar] [CrossRef]
  53. Song, N.; Zhang, H.; Zhao, T. Insights into the phylogeny of Hemiptera from increased mitogenomic taxon sampling. Mol. Phylogenet. Evol. 2019, 137, 236–249. [Google Scholar] [CrossRef] [PubMed]
  54. Zhong, L.K.; Li, H.X.; Yu, X.F.; Yang, M.F. Complete mitochondrial genome sequence of Cicadella viridis (Hemiptera: Cicadellidae: Cicadellinae). Mitochondrial DNA Part B Resour. 2019, 4, 1287–1288. [Google Scholar] [CrossRef] [Green Version]
  55. Zhong, L.K.; Yang, M.F.; Yu, X.F. The mitochondrial genome of Cofana yasumatsui (Hemiptera: Cicadellidae: Cicadellinae). Mitochondrial DNA Part B Resour. 2020, 5, 1075–1076. [Google Scholar] [CrossRef] [Green Version]
  56. Yu, P.F.; Wu, G.H.; Han, B.Y. Sequencing and analysis of the mitochondrial genome of Kolla paulula (Walker) (Hemiptera: Cicadellidae). J. Anhui Agric. Univ. 2017, 44, 874–881. [Google Scholar]
  57. Wang, J.J.; Li, H.; Dai, R.H. Complete mitochondrial genome of Taharana fasciana (Insecta, Hemiptera: Cicadellidae) and comparison with other Cicadellidae insects. Genetica 2017, 145, 593–602. [Google Scholar] [CrossRef]
  58. Wang, J.J.; Wu, Y.F.; Dai, R.H.; Yang, M.F. Comparative mitogenomes of six species in the subfamily Iassinae (Hemiptera: Cicadellidae) and phylogenetic analysis. Int. J. Biol. Macromol. 2020, 149, 1294–1303. [Google Scholar] [CrossRef]
  59. Du, Y.; Dietrich, C.H.; Dai, W. Complete mitochondrial genome of Macrosteles quadrimaculatus (Matsumura) (Hemiptera: Cicadellidae: Deltocephalinae) with a shared tRNA rearrangement and its phylogenetic implications. Int. J. Biol. Macromol. 2019, 122, 1027–1034. [Google Scholar] [CrossRef] [PubMed]
  60. Zhou, N.; Dong, Y.L.; Qiao, P.P.; Yang, Z.F. Complete Mitogenomic Structure and Phylogenetic Implications of the Genus Ostrinia (Lepidoptera: Crambidae). Insects 2020, 11, 232. [Google Scholar] [CrossRef] [Green Version]
  61. Xu, D.L.; Yu, T.H.; Zhang, Y.L. Characterization of the Complete Mitochondrial Genome of Drabescus ineffectus and Roxasellana stellata (Hemiptera: Cicadellidae: Deltocephalinae: Drabescini) and Their Phylogenetic Implications. Insects 2020, 11, 534. [Google Scholar] [CrossRef] [PubMed]
  62. Huang, W.; Zhang, Y. Characterization of two complete mitochondrial genomes of Ledrinae (Hemiptera: Cicadellidae) and phylogenetic analysis. Insects 2020, 11, 609. [Google Scholar] [CrossRef] [PubMed]
  63. Zhou, X.; Dietrich, C.H.; Huang, M. Characterization of the complete mitochondrial genomes of two species with preliminary investigation on phylogenetic status of zyginellini (Hemiptera: Cicadellidae: Typhlocybinae). Insects 2020, 11, 684. [Google Scholar] [CrossRef] [PubMed]
  64. Godeiro, N.N.; Pacheco, G.; Liu, S.; Gioia Cipola, N.; Berbel-Filho, W.M.; Zhang, F.; Gilbert, M.T.P.; Bellini, B.C. Phylogeny of Neotropical Seirinae (Collembola, Entomobryidae) based on mitochondrial genomes. Zool. Scr. 2020, 49, 329–339. [Google Scholar] [CrossRef]
  65. Yuan, L.; Liu, H.; Ge, X.; Yang, G.; Xie, G.; Yang, Y. A Mitochondrial Genome Phylogeny of Cleridae (Coleoptera, Cleroidea). Insects 2022, 13, 118. [Google Scholar] [CrossRef] [PubMed]
  66. Hebert, P.D.N.; Ratnasingham, S.; DeWaard, J.R. Barcoding animal life: Cytochrome coxidase subunit 1 divergences among closely related species. Proc. R. Soc. B Biol. Sci. 2003, 270 (Suppl. 1), S96–S99. [Google Scholar] [CrossRef] [Green Version]
  67. Hajibabaei, M.; Janzen, D.H.; Burns, J.M.; Hallwachs, W.; Hebert, P.D.N. DNA barcodes distinguish species of tropical Lepidoptera. Proc. Natl. Acad. Sci. USA 2006, 103, 968–971. [Google Scholar] [CrossRef] [Green Version]
  68. Taleb, M.; Tail, G.; Açıkgöz, H.N. DNA barcoding of Stearibia nigriceps (Meigen) and Piophila casei (Linnaeus) (Diptera: Piophilidae) from Algeria and the first African report of Stearibia nigriceps. Int. J. Legal Med. 2020, 134, 895–902. [Google Scholar] [CrossRef]
  69. Ma, J.; Liu, J.; Shen, Y.; Fan, Z.; Yue, B.; Zhang, X. Population genetic structure and intraspecific genetic distance of Periplaneta americana (Blattodea: Blattidae) based on mitochondrial and nuclear DNA markers. Ecol. Evol. 2019, 9, 12928–12939. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  70. Zhao, Z.H.; Cui, B.Y.; Li, Z.H.; Jiang, F.; Yang, Q.Q.; Kucerova, Z.; Stejskal, V.; Opit, G.; Cao, Y.; Li, F.J. The establishment of species-specific primers for the molecular identification of ten stored-product psocids based on ITS2 rDNA. Sci. Rep. 2016, 6, 21022. [Google Scholar] [CrossRef]
  71. Haiyan, Z.; Guanghui, C.; Abduklimu, G.; Xiuying, Z.; Yutao, W. Cicadellidae in Some Areas of Xinjiang: Molecular Identification Based on mtDNA COI Gene. Chin. Agric. Sci. Bull. 2021, 37, 119–129. (In Chinese) [Google Scholar]
  72. Wang, X.Y.; Wang, J.J.; Fan, Z.H.; Dai, R.H. Complete mitogenome of Olidiana ritcheriina (Hemiptera: Cicadellidae) and phylogeny of Cicadellidae. PeerJ 2019, 7, e8072. [Google Scholar] [CrossRef] [Green Version]
  73. Wang, J.J.; Yang, M.F.; Dai, R.H.; Li, H.; Wang, X.Y. Characterization and phylogenetic implications of the complete mitochondrial genome of Idiocerinae (Hemiptera: Cicadellidae). Int. J. Biol. Macromol. 2018, 120, 2366–2372. [Google Scholar] [CrossRef] [PubMed]
  74. Wang, J.J.; Li, D.F.; Li, H.; Yang, M.F.; Dai, R.H. Structural and phylogenetic implications of the complete mitochondrial genome of Ledra auditura. Sci. Rep. 2019, 9, 15746. [Google Scholar] [CrossRef] [PubMed]
  75. Lin, S.H.; Huang, M.; Zhang, Y.L. Structural features and phylogenetic implications of 11 new mitogenomes of Typhlocybinae (Hemiptera: Cicadellidae). Insects 2021, 12, 678. [Google Scholar] [CrossRef] [PubMed]
  76. Yu, T.H.; Zhang, Y.L. Two complete mitochondrial genomes of Mileewinae (Hemiptera: Cicadellidae) and a phylogenetic analysis. Insects 2021, 12, 668. [Google Scholar] [CrossRef]
  77. He, H.; Yang, M. Characterization of the leafhopper mitogenome of Mileewa alara (Hemiptera: Cicadellidae: Mileewinae) and its phylogenetic analysis. Mitochondrial DNA Part B Resour. 2021, 6, 1265–1266. [Google Scholar] [CrossRef] [PubMed]
  78. He, H.; Yang, M. Characterization and phylogenetic analysis of the mitochondrial genome of Mileewa ponta (Hemiptera: Cicadellidae: Mileewinae). Mitochondrial DNA Part B Resour. 2020, 5, 2994–2995. [Google Scholar] [CrossRef] [PubMed]
Figure 1. Characteristics of 12 newly sequenced mitochondrial genomes of Atkinsoniella species. Protein coding and ribosomal genes are shown with standard abbreviations. Transfer RNA (tRNA) genes are indicated using the IUPAC-IUB single letter amino acid codes (L1, CUN; L2, UUR; S2). The green circles in front of the species names represent circular mitochondrial genomes. The sharp angle of color blocks on the right and left indicates that the genes were on the J-strand and N-strand, respectively.
Figure 1. Characteristics of 12 newly sequenced mitochondrial genomes of Atkinsoniella species. Protein coding and ribosomal genes are shown with standard abbreviations. Transfer RNA (tRNA) genes are indicated using the IUPAC-IUB single letter amino acid codes (L1, CUN; L2, UUR; S2). The green circles in front of the species names represent circular mitochondrial genomes. The sharp angle of color blocks on the right and left indicates that the genes were on the J-strand and N-strand, respectively.
Insects 13 00254 g001
Figure 2. The relative synonymous codon usage (RSCU) of PCGs in the 12 newly sequenced mitogenomes. The x-axis and y-axis indicate the hierarchical clustering of codon frequencies and Atkinsoniella species, respectively.
Figure 2. The relative synonymous codon usage (RSCU) of PCGs in the 12 newly sequenced mitogenomes. The x-axis and y-axis indicate the hierarchical clustering of codon frequencies and Atkinsoniella species, respectively.
Insects 13 00254 g002
Figure 3. The evolutionary rate of each PCG in the mitogenomes of 32 Cicadellidae species in this study. Ka is the nonsynonymous nucleotide substitutions per nonsynonymous site, Ks is the synonymous nucleotide substitutions per synonymous site, and Ka/Ks is the ratio of the rate of nonsynonymous substitutions to the rate of synonymous substitutions.
Figure 3. The evolutionary rate of each PCG in the mitogenomes of 32 Cicadellidae species in this study. Ka is the nonsynonymous nucleotide substitutions per nonsynonymous site, Ks is the synonymous nucleotide substitutions per synonymous site, and Ka/Ks is the ratio of the rate of nonsynonymous substitutions to the rate of synonymous substitutions.
Insects 13 00254 g003
Figure 4. Boxplot showing the P-distance of each PCG in all 32 species studied. Outlier values are represented by black dots.
Figure 4. Boxplot showing the P-distance of each PCG in all 32 species studied. Outlier values are represented by black dots.
Insects 13 00254 g004
Figure 5. Predicted secondary cloverleaf structures for the tRNAs of the newly sequenced Atkinsoniella species. The tRNA arms are illustrated as for trnV. Dashes (–), solid dots (•), and pluses (+) indicate the Watson–Crick base pairings, G–U bonds, and mismatches, respectively. The conserved and variable sites among the 12 Atkinsoniella species are indicated using black and red, respectively.
Figure 5. Predicted secondary cloverleaf structures for the tRNAs of the newly sequenced Atkinsoniella species. The tRNA arms are illustrated as for trnV. Dashes (–), solid dots (•), and pluses (+) indicate the Watson–Crick base pairings, G–U bonds, and mismatches, respectively. The conserved and variable sites among the 12 Atkinsoniella species are indicated using black and red, respectively.
Insects 13 00254 g005
Figure 6. The phylogenetic relationships of Cicadellinae using the Bayesian inference (BI) method based on the concatenated amino-acid sequences of the PCG (cds_faa) dataset. Numbers on branches are posterior probabilities (PP). The blue clades represent the 12 newly sequenced Atkinsoniella mitogenomes in this study.
Figure 6. The phylogenetic relationships of Cicadellinae using the Bayesian inference (BI) method based on the concatenated amino-acid sequences of the PCG (cds_faa) dataset. Numbers on branches are posterior probabilities (PP). The blue clades represent the 12 newly sequenced Atkinsoniella mitogenomes in this study.
Insects 13 00254 g006
Figure 7. The dorsal views of 12 Atkinsoniella species. (A) Atkinsoniella aurantiaca; (B) A. curvata; (C) A. flavipenna; (D) A. longiuscula; (E) A. thalia; (F) A. thaloidea; (G) A. tiani; (H) A. uniguttata; (I) A. warpa; (J) A. wui; (K) A. xanthoabdomena; (L) A. yunnanana.
Figure 7. The dorsal views of 12 Atkinsoniella species. (A) Atkinsoniella aurantiaca; (B) A. curvata; (C) A. flavipenna; (D) A. longiuscula; (E) A. thalia; (F) A. thaloidea; (G) A. tiani; (H) A. uniguttata; (I) A. warpa; (J) A. wui; (K) A. xanthoabdomena; (L) A. yunnanana.
Insects 13 00254 g007
Figure 8. The paraphysis and aedeagus of 12 Atkinsoniella species. (A) Atkinsoniella aurantiaca; (B) A. curvata; (C) A. flavipenna; (D) A. longiuscula; (E) A. thalia; (F) A. thaloidea; (G) A. tiani; (H) A. uniguttata; (I) A. warpa; (J) A. wui; (K) A. xanthoabdomena; (L) A. yunnanana.
Figure 8. The paraphysis and aedeagus of 12 Atkinsoniella species. (A) Atkinsoniella aurantiaca; (B) A. curvata; (C) A. flavipenna; (D) A. longiuscula; (E) A. thalia; (F) A. thaloidea; (G) A. tiani; (H) A. uniguttata; (I) A. warpa; (J) A. wui; (K) A. xanthoabdomena; (L) A. yunnanana.
Insects 13 00254 g008
Table 1. Mitochondrial genomes used for the phylogenetic analyses in this study.
Table 1. Mitochondrial genomes used for the phylogenetic analyses in this study.
SubfamilySpeciesSize (bp)A + T (%)Accession NumberReference
CicadellinaeAtkinsoniella aurantiaca15,39479.2OL677863This Study
CicadellinaeAtkinsoniella curvata15,98879.3OL677864This Study
CicadellinaeAtkinsoniella flavipenna15,67179.1OL677864This Study
CicadellinaeAtkinsoniella grahami15,62178.6MW533712[49]
CicadellinaeAtkinsoniella longiuscula15,81579.2OL677864This Study
CicadellinaeAtkinsoniella thalia15,03477.8OL677864This Study
CicadellinaeAtkinsoniella thaloidea15,57178.2OL677868This Study
CicadellinaeAtkinsoniella tiani15,84179.2OL677868This Study
CicadellinaeAtkinsoniella uniguttata15,78479.8OL677868This Study
CicadellinaeAtkinsoniella warpa15,85279.0OL677868This Study
CicadellinaeAtkinsoniella wui15,15977.7OL677868This Study
CicadellinaeAtkinsoniella xanthoabdomena15,46378.9OL677868This Study
CicadellinaeAtkinsoniella xanthonota15,89578.4MW533713[49]
CicadellinaeAtkinsoniella yunnanana15,87579.0OL677874This Study
CicadellinaeBothrogonia ferruginea15,26276.5KU167550[50]
CicadellinaeBothrogonia qiongana15,78876.9NC_049894[51]
CicadellinaeBothrogonia tongmaiana15,53976.5MT500857[52]
CicadellinaeBothrogonia yunana15,58576.4NC_049896[52]
CicadellinaeCicadella sp.15,13178.2MK251092[53]
CicadellinaeCicadella viridis13,46178.1KY752061Unpublished
CicadellinaeCicadella viridis15,89178.8MK335936[54]
CicadellinaeCofana unimaculata14,86577.6MK251095[53]
CicadellinaeCofana yasumatsui15,01977.2NC_049087[55]
CicadellinaeCuerna septentrionalis12,42077.1MK251107[53]
CicadellinaeDraeculacephala crassicornis15,21476.9MK251091[53]
CicadellinaeDraeculacephala sp.12,40075.5MK251110[53]
CicadellinaeHomalodisca sp.14,92178.7MK251090[53]
CicadellinaeHomalodisca vitripennis15,30478.4NC_006899Unpublished
CicadellinaeKolla paulula13,57973.3KY680277[56]
CicadellinaeParacrocampsa sp.14,74975.6MK251085[53]
CoelidiinaeTaharana fasciana15,16177.9NC_036015[57]
IassinaeIassus dorsalis15,17680.1NC_046066[58]
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Jiang, Y.; Li, H.-X.; Yu, X.-F.; Yang, M.-F. Comparative Analysis of Mitochondrial Genomes among Twelve Sibling Species of the Genus Atkinsoniella Distant, 1908 (Hemiptera: Cicadellidae: Cicadellinae) and Phylogenetic Analysis. Insects 2022, 13, 254. https://doi.org/10.3390/insects13030254

AMA Style

Jiang Y, Li H-X, Yu X-F, Yang M-F. Comparative Analysis of Mitochondrial Genomes among Twelve Sibling Species of the Genus Atkinsoniella Distant, 1908 (Hemiptera: Cicadellidae: Cicadellinae) and Phylogenetic Analysis. Insects. 2022; 13(3):254. https://doi.org/10.3390/insects13030254

Chicago/Turabian Style

Jiang, Yan, Hao-Xi Li, Xiao-Fei Yu, and Mao-Fa Yang. 2022. "Comparative Analysis of Mitochondrial Genomes among Twelve Sibling Species of the Genus Atkinsoniella Distant, 1908 (Hemiptera: Cicadellidae: Cicadellinae) and Phylogenetic Analysis" Insects 13, no. 3: 254. https://doi.org/10.3390/insects13030254

APA Style

Jiang, Y., Li, H. -X., Yu, X. -F., & Yang, M. -F. (2022). Comparative Analysis of Mitochondrial Genomes among Twelve Sibling Species of the Genus Atkinsoniella Distant, 1908 (Hemiptera: Cicadellidae: Cicadellinae) and Phylogenetic Analysis. Insects, 13(3), 254. https://doi.org/10.3390/insects13030254

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop