Next Article in Journal
A Multimodal Attract-and-Kill Device for the Asian Citrus Psyllid Diaphorina citri (Hemiptera: Liviidae)
Previous Article in Journal
Susceptibility of Various Developmental Stages of the Fall Armyworm, Spodoptera frugiperda, to Entomopathogenic Nematodes
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Insects
Volume 11
Issue 12
10.3390/insects11120869
insects-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

The Complete Mitochondrial Genome of Four Hylicinae (Hemiptera: Cicadellidae): Structural Features and Phylogenetic Implications

by
Jiu Tang
†,
Weijian Huang
† and
Yalin Zhang
*
Key Laboratory of Plant Protection Resources and Pest Management, Ministry of Education, Entomological Museum, College of Plant Protection, Northwest A&F University, Yangling 712100, China
*
Author to whom correspondence should be addressed.
†
These two authors contributed equally in this study.
Insects 2020, 11(12), 869; https://doi.org/10.3390/insects11120869
Submission received: 19 November 2020 / Revised: 4 December 2020 / Accepted: 4 December 2020 / Published: 7 December 2020
Browse Figures
Review Reports Versions Notes

Abstract

:

Simple Summary

Hylicinae, containing 43 described species in 13 genera of two tribes, is one of the most morphologically unique subfamilies of Cicadellidae. Phylogenetic studies on this subfamily were mainly based on morphological characters or several gene fragments and just involved single or two taxa. No mitochondrial genome was reported in Hylicinae before. Therefore, we sequenced and analyzed four complete mtgenomes of Hylicinae (Nacolus tuberculatus, Hylica paradoxa, Balala fujiana, and Kalasha nativa) for the first time to reveal mtgenome characterizations and reconstruct phylogenetic relationships of this group. The comparative analyses showed the mtgenome characterizations of Hylicinae are similar to members of Membracoidea. In phylogenetic results, Hylicinae was recovered as a monophyletic group in Cicadellidae and formed to the sister group of Coelidiinae + Iassinae. These results provide the comprehensive framework and worthy information toward the future researches of this subfamily.

Abstract

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

1. Introduction

The infraorder Cicadomorpha in Hemiptera includes three superfamilies: Cicadoidea, Cercopoidea, and Membracoidea. The last one, comprising Cicadellidae, Myerslopiidae, and three treehopper families (Membracidae, Melizoderidae, and Aetalionidae), is the most diverse among these groups [1]. The leafhopper family Cicadellidae, containing >2600 genera and about 21,000 species distributed worldwide, is one of the largest families of Hemiptera [2].
One such group, the hylicine leafhoppers, is one of the most morphologically unique cicadellids and was once treated as a tribe of Cicadellinae [3] as well as an independent family with respect to Cicadellidae [4]. However, it is currently recognized as a cicadellid subfamily based on morphological and molecular evidence [1,5,6]. Among subfamilies of Cicadellidae, Hylicinae can be distinguished from the others easily by the following combination of characters: body brown to black covered with scale-like setae, head usually produced, ocelli on the crown close to eyes, forewing appendix relatively wide. This subfamily is very small and is comprised of 43 described species in 13 genera assigned to two tribes [6,7]. Most (11) genera of Hylicinae are restricted to the Oriental and eastern Palearctic regions. The two genera Wolfella Spinola, 1850, and Karasekia Melichar, 1926, are distributed in the Ethiopian region [8]. Hylicines have been found on grasses and trees in moistly forested areas and on chili in farmland [9,10]. Many leafhoppers spread viruses causing plant diseases, but that ability has not been found in this group yet [11].
Research on this group has mainly focused on species taxonomy [9,12,13,14,15,16]. Phylogenetic analyses are relatively rare due to the difficulty of obtaining specimens due to their limited of distribution, scarcity of species, etc. One phylogenetic study of Membracoidea based on 28S rDNA included one hylicine exemplar [5]. Results showed that the position of Hylicinae varied and the relationship of Hylicinae to other subfamilies had not been resolved. Recently, a phylogenetic analysis using DNA sequences including nearly hundreds of thousands aligned nucleotides from hundreds of loci also suggested that Hylicinae, represented by two species of one tribe, is monophyletic with strong branch support. The results were also unable to resolve the relationships of Hylicinae and related subfamilies [17]. The monophyly and the phylogenetic status in Cicadellidae of Hylicinae need further study based on more exemplars and using different methods.
The typical insect mtgenome is a double-stranded circular molecule and generally contains 37 genes: 13 protein coding genes (PCGs), 22 transfer RNA genes (tRNAs), two ribosomal RNA genes (rRNAs), and one non-coding region (CR) [18,19,20]. The genome-level features that mtgenomes can offer include base composition, gene arrangement, genetic codon variation, and secondary structures of tRNA and rRNA [21,22]. More and more insect mitochondrial genomes have been obtained with the development of sequencing technology. To date, the mtgenomes of Membracoidea have been determined and contributed to GenBank (https://www.ncbi.nlm.nih.gov) more than 150. Unfortunately, there is no mtgenome or even a partial one for the subfamily Hylicinae. Simple genetic structure, small size, strict orthologous genes, high genome copy numbers, less recombination, and fast evolutionary rate are characteristics of mtgenomes [18,19,23]. Those characteristics make mtgenome a new tool. In the past decade, mtgenomes have proven to be a useful source for studies on inter- and intra-specific variation as well as the systematics and phylogeny of arthropods at different levels including leafhoppers [24,25,26,27,28,29,30]. Therefore, mtgenomes may be a better method to analyze the taxonomic and phylogenetic status of Hylicinae. Meanwhile, a lack of mtgenome data for hylicine species impedes a more complete phylogenetic analysis of Membracoidea. So it is necessary and urgent to fill in the gaps in the mtgenome of this subfamily.
In this study, we sequenced and annotated mtgenomes of four Hylicinae leafhoppers (Nacolus tuberculatus (Walker, 1858), Hylica paradoxa Stål, 1863, Balala fujiana Tang and Zhang, 2020 and Kalasha nativa Distant, 1908). Using the mtgenomes of membracoids including these four species, the phylogeny of Membracoidea was reconstructed. The aims of this study were to (1) figure out the mitochondrial structure including gene order, nucleotide composition, codon usage, tRNA secondary structure, gene overlap, and the non-coding control region of hylicine species (2) verify the monophyly of subfamily Hylicinae and explore the phylogenetic relationships of Hylicinae to other major lineages.

2. Materials and Methods

2.1. Sample Collection and DNA Extraction

Nacolus tuberculatus was captured in Qingquan Gully (109.4412° E, 36.2472° N), Ganquan County, Yan’an City, Shaanxi Province, China, on 19 June 2019. Hylica paradoxa was collected in Menglun County (101.2552° E, 21.9412° N), Yunnan Province, China, on 4 May 2009. The B. fujiana specimen was captured in Mengla County (101.5682° E, 21.4822° N), Yunnan Province, China, on 27 April 2019, while the K. nativa specimen was collected in Jianfengling (110.3491° E, 19.9770° N), Hainan Province, 975 m, China, on 15 August 2010. Four species of Hylicinae were identified by external morphology and genitalia [7,31,32,33]. All samples were kept in tubes fulling of absolute ethyl alcohol at minus 20 °C. Voucher specimens were stored in the Entomological Museum, Northwest A&F University (NWAFU), Yangling, Shaanxi, China.
Total genomic DNA was extracted using the Easy Pure Genomic DNA Kit following the manufacturer’s instructions from the thoracic muscle.

2.2. Mtgenome Sequences Acquisition and Bioinformatic Analyses

The whole mtgenomes of four species in Hylicinae were generated by NGS. The paired-end clean reads were assembled by employing the mtgenomes of Evacanthus heimianus Kuoh, 1981, (MG813486) and N. tuberculatus (MW218663) as references in Geneious10.0.5. A total of 16,019,948/22,247,400/24,753,222/15,786,834/paired-end clean reads were assembled with trim sequence and under medium-low sensitivity to assemble. Then we generated consensus sequence and used MAFFT (within Geneious 10.0.5., Biomatters, Auckland, New Zealand) to align the consensus sequences and reference sequences for annotation.
The annotation of four mtgenomes was accomplished in Geneious10.0.5. Twenty-two tRNAs were found by MITOS [34] and tRNAscan-SE v1.21 [35] under the invertebrate mitochondrial genetic code. Two rRNAs were predicted by aligning with a reference mtgenomes and its boundaries were defined by the adjacent tRNAs. All 13 PCGs were determined by the ORF Finder employing codon Table 5 (invertebrate mitochondrion) and comparison with the homologous sequence of other Cicadellidae species. The different types of repeat tandem units were identified using the Tandem Repeats Finder server [36]. The mtgenome map was produced using CGView [37].
The analyses of the four newly sequenced mtgenomes of Hylicinae, including the base composition and relative synonymous codon usage (RSCU), were calculated using MEGA v 7.313 [38] incorporated into PhyloSuite v 1.2.2 [39]. The AT and GC base compositional skewness were computed employing the formulas (A – T)/(A + T) and (G – C)/(G + C), respectively [40].
The Pi value (nucleotide diversity), the ratios between non-synonymous (Ka) and synonymous (Ks) substitutions rates for each PCG, and a sliding window analysis (a sliding window of 200 bp and step size of 20 bp) of 13 PCGs among four Hylicinae species were conducted by DnaSP 6.0 [41]. Mean genetic distances among species were analyzed under Kimura-2-parameter using MEGA v 7.313 [38]. The genetic distances, and Ka/Ks ratios were graphically produced by Prism 6.01. The four mtgenome sequences of Hylicinae (N. tuberculatus, H. paradoxa, B. fujiana, and K. nativa) were registered in GenBank as MW218663, MW218660, MW218661, and MW218662, respectively (Table 1).

2.3. Phylogenetic Analyses

Among the 150 available Membracoidea mtgenomes, there are some repeats of same species. We chose the most complete and highest quality one among these repeats. The aim of this study is to explore the monophyly of Hylicinae and the phylogenetic relationships of Hylicinae to other major subfamilies. So we should choose moderate species representing each subfamily/tribe rather than all. In phylogenetic analyses, we chose 50 Membracoidea samples (covering 44 leafhoppers, six treehoppers) as ingroup and eleven previously available samples of the other three suborders in Hemiptera as outgroups. Four mtgenomes of Hylicinae species (H. paradoxa, K. nativa, N. tuberculatus, and B. fujiana) were sequenced in this study while the other sequences were acquired from the GenBank (Table 1).
PhyloSuite v 1.2.2 was used for genes extraction [39]. All PCGs and rRNAs were aligned individually with the G-INS-I algorithm and Q-INS-I algorithm respectively in MAFFT v 7.313 online service [69]. Ambiguous sites and gaps in the alignments were removed using GBlocks v 0.91b [70].
Three optimized datasets were merged in PhyloSuite v 1.2.2: (1) PCG123 matrix, including 3 codon positions of 13 PCGs (10,725 bp); (2) PCG123RNA matrix, including 3 codon positions of the 13 PCGs and two rRNAs (11,650 bp); and (3) AA matrix, including amino acid sequences of 13 PCGs (3277 bp) [39]. Maximum likelihood (ML) and Bayesian inference (BI) analyses were used for phylogenetic reconstruction of all matrices. The best partitioning strategies were inferred in PartitionFinder 2.1.1 [71] incorporated into PhyloSuite v1.2.2 [39] with the greedy algorithm and BIC criterion (Tables S6 and S7). ML analyses were implemented in IQ-TREE v 1.6.8 [72] under an edge-linked partition model. Branch support analyses were conducted under 10,000 ultrafast bootstrap replicates (UFB) [73]. BI analyses were conducted in MrBayes 3.2.6 [74], as implemented in the CIPRES Science Gateway [75]. Each BI analysis involved 5,000,000–20,000,000 generations. The convergence of the independent runs was indicated by a standard deviation of split frequencies < 0.01 and effective sample size (ESS) > 200. A consensus tree was computed from the trees after the initial 25% trees of each MCMC run were discarded as burn-in.

3. Results and Discussion

3.1. Genome Organization and Base Composition

The mtgenomes of N. tuberculatus (15,737 bp), H. paradoxa (14,762 bp), B. fujiana (16,221 bp), and K. nativa (15,716 bp) were single, closed circular double-stranded molecules (Figure 1 and Figure 2). Among the four mtgenomes of Hylicinae, B. fujiana was the largest in 16,221 bp, while H. paradoxa had the smallest mtgenome of 14,762 bp. Length variation of Membracoidea mtgenomes appears causing by variation in the length of the control region [59,61]. Typical 37 animal mitochondrial genes (13 PCGs, 22 tRNAs, and two rRNAs) and one non-coding region were detected in all mtgenomes of Hylicinae (Figure 1 and Figure 2). The J-strand encoded 23 genes (nine PCGs and 14 tRNAs), while the remaining 14 genes (four PCGs, eight tRNAs, and two rRNAs) were transcribed on the J-strand (Tables S1–S4). All the sequenced mtgenomes of Hylicinae possessed the same ancestral gene order as the typical insect: Drosophila yakuba Burla, 1954 (Diptera: Drosophilidae) [76]. All Hylicinae species shared seven conserved overlap regions in trnI-trnQ (TTG), trnW-trnC (AAGTCTT), trnY-cox1 (AT), trnK-trnD (G), atp8-atp6 (TGAAAATGATAA), trnN-trnS (A) and trnS2-nad1 (TTAATAACTT) (Tables S1–S4). Among these conserved overlap regions, trnW-trnC s (AAGTCTT) was also revealed in Ledrinae [30].
As in other published subfamilies, the AT content of the whole mtgenomes of Hylicinae indicated a strong AT bias, ranging from 75.8% in H. paradoxa to 77.4% in B. fujiana (Table S5). Comparing the AT content of the whole mtgenome, control region, PCGs, tRNAs, and rRNAs, the rRNAs was the highest while the PCGs was the lowest for N. tuberculatus and K. nativa. For H. paradoxa, the AT content of the control region was the highest, and the PCGs was the lowest. In addition, the AT content of the rRNAs was the highest, and the control region was the lowest of B. fujiana. The composition skew analyses showed a positive AT skew and a negative GC skew in all four whole mtgenomes (Table S5).

3.2. Protein-Coding Genes and Codon Usage

All 13 PCGs were detected in the newly sequenced mtgenomes and their length ranging from 10,932 bp in H. paradoxa to 10,962 bp in K. nativa (Table S5). In all sequenced mtgenomes, nine PCGs were assigned on the J-strand, the remaining four were located on the N-strand (Tables S1–S4). The AT skews of the PCGs ranged from −0.104 to −0.126 (Table S5).
The majority of PCGs in the four newly sequenced mtgenomes used standard initial codons for the genetic code 5 (invertebrate mitochondrion): ATN, except for nad5 in three of the four newly sequenced Hylicinae which started with TTG (N. tuberculatus, H. paradoxa, and B. fujiana, Tables S1–S3), as seen in previous leafhoppers [29,30,61], while nad5 in K. nativa started with ATT (Table S4). After counting the terminating codon, the majority of PCGs of four mtgenomes stopped with codon TAA or TAG, whereas the remaining PCGs were stopped with a single T (Tables S1–S4). Such incomplete termination codons also occur in other leafhoppers [29,30,61], and could accomplish completion by posttranscriptional polyadenylation [77].
The relative synonymous codon usage (RSCU) of the four Hylicinae mtgenomes were summarized in Figure 3. The four most frequently used amino acids codons were UUU (Phe), UUA (Leu2), AUA (Met), and AUU (Ile), and all of them are composed with only A and/or U.

3.3. Transfer and Ribosomal RNA Genes

All 22 tRNA genes of N. tuberculatus H. paradoxa, B. fujiana, and K. nativa mtgenomes were determined (Tables S1–S4). The tRNAs length of these four mtgenomes was 1429 bp in N. tuberculatus, 1416 bp in H. paradoxa, 1416 bp in B. fujiana, and 1425 bp in K. nativa. The AT content of tRNA genes was slightly higher than that of the PCGs and lesser than that of the rRNAs, ranging from 75.4% to 78.9% (Table S5). Each position of 22 tRNAs was same as relative positions in D. yakuba [76] and previously published cicadellid species (Table 1), except for a few leafhoppers containing tRNA rearrangements [59,60,78,79]. The sizes of these 22 tRNAs ranging from 61 (trnW, trnC, trnD) to 70 bp (trnL) in B. fujiana, 61 (trnD, trnA, trnS, trnV) to 70 bp (trnQ, trnM) in H. paradoxa, 61 (trnA, trnS, trnV) to 71 bp (trnM, trnP) in K. nativa, and from 61 (trnG, trnA, trnE, trnH) to 71 bp (trnQ, trnM) in N. tuberculatus (Table S1–S4). As shown in Figures S1–S4, all of the 22 tRNAs had the same cloverleaf secondary structure in other Cicadellidae species [30,54,59,60,61,78]. The trnS1 lost the DHU arm and this phenomenon probably happened very early during the evolution of the Metazoa [77]. Seven kinds of nonmatched base pairs (GU, UU, AA, AC, AG, CU, and a single C) were found in the arm structures of the tRNAs (Figures S1–S4), while such nonmatched base pairs were also found in other published cicadellid species [29,30,43,61]. In summary, a large number of GU mismatches were the most common, which were also overrepresented as revealed in other leafhoppers [30,43,54,59,61,63].
Both of lrRNA (rrnL) and srRNA (rrnS) genes were encoded from the N-strand in all Hylicinae mtgenomes. Among four mtgenomes of Hylicinae, the length of rrnL was 1180 to 1207bp, resided between trnL1 (CUN) and trnV, and the rrnS was 730 to 740bp and located between trnV and the control region (Tables S1–S4).

3.4. Control Region

In all mtgenomes, a large non-coding region was found (Figure 4) located between rrnS and trnI and ranges from 531 bp to 1976 bp. The A + T contents were 73.5% to 86.6% in the Hylicinae mtgenomes. Each Hylicinae had different repeat units. Three Hylicinae species had two kinds of short repeat tandem units of 11 and 131 bp in size, respectively, while H. paradoxa only had a tandem repeat ranging from 16 to 17 bp (Figure 4). In addition, except in H. paradoxa, the other three Hylicinae had Poly T/A stretches. The results showed the number of absolute tandem repeat units in control region of each Hylicinae mtgenome are different.

3.5. Pi values and Selection Pressures

The sliding window analysis indicated that Pi values (nucleotide diversity) were highly variable among the 13 PCGs of these Hylicinae mtgenomes (Figure 5A). The PCGs atp8, nad2, and atp6 appear to have a comparatively high Pi values of 0.327, 0.306, and 0.261, respectively, while the PCGs cox1, nad4l and nad1 presented relatively low Pi values of 0.185, 0.176 and 0.173 (Figure 5A). Meanwhile, the similar tendencies, that was the PCGs atp8, nad2, and atp6 with high distances of 0.599, 0.519, and 0.406, and the PCGs cox1, nad4l, and nad1 with low distances of 0.249, 0.234, and 0.227, severally, were produced in the analysis of genetic distance (Figure 5B).
To estimate the selection pressures of the 13 PCGs in Hylicinae species, the mean substitution rates of synonymous (Ks), non-synonymous (Ka) of the 13 PCGs among these Hylicinae mtgenomes were calculated and ranging from 0.150 to 0.984 (Figure 5B), indicating that all PCGs are under purifying selection. The Ka/Ks ratios of PCGs atp8, nad2, and nad6 were relatively high of 0.984, 0.609, and 0.518, while the Ka/Ks ratios of cox1, cox2, and cox3 were comparatively low of 0.150, 0.244, and 0.278, severally (Figure 5B).
Among the 13 PCGs, it was apparent that cox1 was the slowest-evolving gene, while atp8 was the relatively fastest-evolving gene in these Hylicinae species (Figure 5B). Moreover, the results also suggested that nad2 and atp8 would be capable markers for sibling hylicine species delimitation.

3.6. Phylogenetic Relationships

Regardless of the method used, the analyses employing different datasets (P123, P123R, and AA) resulted highly identical phylogenetic topologies and most nodes obtained strong supports (BS = 100, PP = 1) (Figure 6, Figure 7 and Figure 8). Fulgoroidea was placed at the base of Auchenorrhyncha, whereas the remaining superfamilies of Auchenorrhyncha were split into two clades: Membracoidea and Cicadoidea + Cercopoidea, which was consistent with recent studies using mtgenomes [30] and transcriptomes [80]. Our analyses corroborated the monophyly of Membracoidea [17,30,50,78,80,81]. Deltocephalinae was placed as sister group to the remainders of Membracoidea with strong support, which is also in line with previous studies [17,29,30,51,78]. For about the remaining Membracoidea, treehoppers were recovered as a monophyletic lineage and derived from Cicadellidae, which also has been supported by previous studies [17,29,50,78,81].
Recent phylogenetic studies have improved our knowledge of the subfamily-level phylogeny of paraphyletic Cicadellidae, the monophyly of nine subfamilies [30], as well as the sister-group relationship between treehoppers and Megophthalminae [29,30,78]. Our results were consistent with these findings and indicated the monophyly of expanded subfamilies (eleven) in all phylogenetic trees (Figure 6, Figure 7 and Figure 8) with strong support. Typhlocybinae was sometimes treated as sister group to all other subfamilies except Deltocephalinae [30]. Our results show the same topology in ML/BI analyses based on the datasets of P123/P123R (Figure 6 and Figure 7). However, a different relationship, Typhlocybinae + Mileewinae and Cicadellinae + (Evacanthinae + Ledrinae) formed a monophyletic group in ML/BI trees based on AA (Figure 8). The relationship, (Coelidiinae + Iassinae) + Hylicinae, was stable in our analyses (Figure 6, Figure 7 and Figure 8) and the supports were moderate to strong (BS > 86; PP = 1). The status of Eurymelinae and Mileewinae was discordant based on three datasets (Figure 6, Figure 7 and Figure 8) and the supports were low to moderate (BS = 26–92, PP < 0.92). Based on more subfamilies sampled and using different datasets, our results are somewhat inconsistent with previous molecular phylogenies [30,54,78,81]. Although we selected more comprehensive subfamily-level samples, the very short deep internal branches still appeared in our study with less than maximum support, as has been shown in other recent phylogenomic analyses of Membracoidea [17,30,81]. This phenomenon may indicate ancient rapid radiations during the Cretaceous period [17].
Within Hylicinae, the four species (N. tuberculatus, H. paradoxa, B. fujiana, and K. nativa) representing four genera of the two tribes (H. paradoxa representing tribe Hylicini, the remaining three representing tribe Sudrini), gathered as a monophyletic group with maximum support values. The monophyly of Hylicinae is consistent with previous studies based on morphological, anchored hybrid enrichment genomics [17]. Hylicini performs as a sister group to Sudrini in BI/ML trees based on PCGR/AA with middle to high support, while Sudrini is paraphyletic based on PCG with low support (Figure 6, Figure 7 and Figure 8). In our analyses, except for BI/ML analysis topology based on PCGs, (H. paradoxa + (N. tuberculatus + (B. fujiana + K. nativa))) was reconstructed with high posterior probabilities support values (PP > 0.98) and middle to high bootstrap support values (BS = 83–100). Unfortunately, our study also included limited taxon sampling despite having more than other related phylogenetic studies, so that we are unable to provide a robust test the monophyly of Hylicinae tribes here.

4. Conclusions

In our work, we report the whole mtgenome sequences of N. tuberculatus, H. paradoxa, B. fujiana, and K. nativa for Hylicinae for the first time, comparatively analyzing the mtgenomes within this subfamily and reconstructing the phylogenetic relationships of the Membracoidea. The mtgenomes of Hylicinae are highly conserved in terms of genome size, gene content and arrangement, base composition, PCG codon usage, and secondary structure of tRNAs. The variations among length of Hylicinae mtgenomes is mainly caused by the varying length of the control region. This study also shows that nad2 and atp8 are potential DNA markers for species delimitation for Hylicinae. Bayesian inference and maximum likelihood analyses provide the most comprehensive subfamily-rich phylogenetic analysis of Membracoidea up to now which presents a well-resolved topology and is largely congruent with previous studies. This study also demonstrates that the mtgenome provides new insight into resolving the phylogenetic relationships at family and subfamily levels. All phylogenetic analyses uniformly support Cicadellidae being paraphyletic with respect to treehoppers and all involved cicadellid subfamilies being monophyletic. The monophyly of Hylicinae is further confirmed in our study. This study also recovered a stable relationship for the subfamilies ((Coelidiinae + Iassinae) + Hylicinae) with moderate to high support. The monophyly of the two tribes, Hylicini and Sudrini, and the topology of four hylicine genera (H. paradoxa + (N. tuberculatus + (B. fujiana + K. nativa))) were recovered in most analyses. There are still some specific problems with Membracoidea (such as the exact placement of Eurymelinae and Mileewinae) which may be caused by the limited taxon sampling of major lineages. We therefore recommend that further research into the proposed relationships within Membracoidea should include more comprehensive taxon sampling to establish a better framework (selecting all subfamilies and tribes).

Supplementary Materials

The following are available online at https://www.mdpi.com/2075-4450/11/12/869/s1, Table S1: Mtgenome organization of N. tuberculatus, Table S2: Mtgenome organization of H. paradoxa, Table S3: Mtgenome organization of B. fujiana, Table S4: Mtgenome organization of K. nativa, Table S5: Base composition and skewness of four Hylicinae mtgenomes, Table S6: Best strategies for Bayesian inference (BI) analysis, Table S7: Best strategies for Maximum likelihood (ML) analysis, Table S8: Detail of software and kit, Figure S1: Secondary structure for the tRNAs of N. tuberculatus, Figure S2: Secondary structure for the tRNAs of H. paradoxa, Figure S3: Secondary structure for the tRNAs of B. fujiana, Figure S4: Secondary structure for the tRNAs of K. nativa.

Author Contributions

Conceptualization, J.T., W.H., and Y.Z.; specimen identification, J.T. and Y.Z.; methodology and experiments, J.T. and W.H.; data analysis, J.T. and W.H.; writing—original draft preparation, J.T. and W.H.; writing—review and editing, J.T., W.H., and Y.Z.; funding acquisition, Y.Z. All authors have read and agreed to the published version of the manuscript.

Funding

This study was supported by the National Natural Science Foundation of China (31420103911, 31672339, 31801995).

Acknowledgments

We sincerely thank John Richard Schrock (Emporia State University, KS, USA) for reviewing the manuscript of this article.

Conflicts of Interest

The authors declare no conflicting interest.

References

  1. 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]
  2. Camisao, B.M.; Cavichioli, R.R.; Takiya, D.M. Eight new species of Oragua Melichar, 1926 (Insecta: Hemiptera: Cicadellidae) from Amazonas State, Brazil, with description of the female terminalia of Oragua jurua Young, 1977, and new records for the genus. Zootaxa 2014, 3841, 501–527. [Google Scholar] [CrossRef] [PubMed]
  3. Hamilton, K.G.A. Classification, morphology and phylogeny of the family Cicadellidae (Rynchota: Homoptera). In Proceedings of the 1st Internal Workshop on Leafhoppers and Planthoppers of Economic Importance; Commonwealth Institute of Entomology: London, UK, 1983; pp. 15–37. [Google Scholar]
  4. Evans, J.W. A natural classification of leaf-hoppers (Homoptera, Jassoidea). Part 2: Aetalionidae, Hylicidae, Eurymelidae. Trans. R. Entomol. Soc. Lond. 1946, 97, 39–54. [Google Scholar] [CrossRef]
  5. Dietrich, C.H.; Rakitov, R.A.; Holmes, J.L.; Black, W.C. Phylogeny of the major lineages of Membracoidea (Insecta: Hemiptera: Cicadomorpha) based on 28S rDNA sequences. Mol. Phylogenet. Evol. 2001, 18, 293–305. [Google Scholar] [CrossRef] [PubMed]
  6. Oman, P.W.; Knight, W.J.; Nielson, M.W. Leafhoppers (Cicadellidae): A Bibliography, Generic Check-list, and Index to the World Literature 1956–1985; CAB International Institute of Entomology: Wallingford, CT, USA, 1990; p. 368. [Google Scholar]
  7. Tang, J.; Zhang, Y. Review of the oar-head leafhopper genus Nacolus Jacobi (Hemiptera: Cicadellidae: Hylicinae). Zootaxa 2019, 4571, 58–72. [Google Scholar] [CrossRef]
  8. Metcalf, Z.P. General Catalogue of the Homoptera; Fascicle VI. Cicadelloidea. Part 2. Hylicidae; United States Department of Agriculture-Agricultural Research Service: Washington, DC, USA, 1962; pp. 1–18.
  9. Kramer, J.P. A review of the African leafhopper genus Wolfella Spinola (Homoptera: Cicadellidae). Trans. Am. Entomol. Soc. 1965, 91, 167–180. [Google Scholar]
  10. Viraktamath, C.A.; Webb, M.D. A review of the Indian hylicine leafhopper genus Traiguma Dist., with description of a new species (Insecta, Homoptera, Auchenorrhyncha: Cicadellidae). Reichenbachia 1991, 28, 123–133. [Google Scholar]
  11. Nielson, M.W. The Leafhopper Vectors of Phytopathogenic Viruses (Homoptera, Cicadellidae): Taxonomy, Biology, and Virus Transmission; No. 1382; US Department of Agriculture: Washington, DC, USA, 1968.
  12. Linnavuori, R. Revisional studies on African leafhoppers (Homoptera Cicadelloidea). Rev. Zool. Bot. Afr. 1972, 86, 195–251. [Google Scholar]
  13. Kramer, J.P. A review of the Oriental leafhopper genus Sudra Distant (Homoptera: Cicadellidae: Hylicinae). Proc. Biol. Soc. Wash. 1964, 77, 47–52. [Google Scholar]
  14. Shen, L.; Zhang, Y. A New species of the genus Assiringia Stål (Homoptera: Cicadellidae: Hylicinae) from Tibet, China. Entomotaxonomia 1995, 17, 106–108. [Google Scholar]
  15. Shen, L.; Zhang, Y. A new species and a new Chinese record of the genus Kalasha Distant (Homoptera: Cicadellidae: Hylicinae). Entomotaxonomia 1995, 17, 185–188. [Google Scholar]
  16. Shen, L.; Zhang, Y. Two new species of the genus Balala Distant (Homoptera: Cicadellidae: Hylicinae) from China. Entomotaxonomia 1995, 17, 271–275. [Google Scholar]
  17. Dietrich, C.H.; Allen, J.M.; Lemmon, A.R.; Lemmon, E.M.; Takiya, D.M.; Evangelista, O.; Walden, K.K.O.; Grady, P.G.S.; Johnson, K.P. Anchored hybrid enrichment-based phylogenomics of leafhoppers and treehoppers (Hemiptera: Cicadomorpha: Membracoidea). Insect Syst. Divers 2017, 1, 57–72. [Google Scholar] [CrossRef]
  18. Boore, L.J. Animal mitochondrial genomes. Nucleic Acids Res. 1999, 27, 1767–1780. [Google Scholar] [CrossRef] [Green Version]
  19. Cameron, S.L. Insect mitochondrial genomics: Implications for evolution and phylogeny. Annu. Rev. Entomol. 2014, 59, 95–117. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  20. Zhang, D.-X.; Hewitt, G.M. Insect mitochondrial control region: A review of its structure, evolution and usefulness in evolutionary studies. Biochem. Syst. Ecol. 1997, 25, 99–120. [Google Scholar] [CrossRef]
  21. Simon, C.; Buckley, T.R.; Frati, F.; Stewart, J.B.; Beckenbach, A.T. Incorporating molecular evolution into phylogenetic analysis, and a new compilation of conserved polymerase chain reaction primers for animal mitochondrial DNA. Annu. Rev. Ecol. Evol. S 2006, 37, 545–579. [Google Scholar] [CrossRef] [Green Version]
  22. Gissi, C.; Iannelli, F.; Pesole, G. Evolution of the mitochondrial genome of Metazoa as exemplified by comparison of congeneric species. Heredity 2008, 101, 301–320. [Google Scholar] [CrossRef] [Green Version]
  23. Shao, R.F.; Dowton, M.; Murrell, A.; Barker, S.C. Rates of gene rearrangement and nucleotide substitution are correlated in the mitochondrial genomes of insects. Mol. Biol. Evol. 2003, 20, 1612–1619. [Google Scholar] [CrossRef] [Green Version]
  24. Du, Z.Y.; Hasegawa, H.; Cooley, J.R.; Simon, C.; Yoshimura, J.; Cai, W.Z.; 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]
  25. 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, 10. [Google Scholar] [CrossRef] [PubMed]
  26. Song, N.; Liang, A.P.; Bu, C.P. A molecular phylogeny of Hemiptera inferred from mitochondrial genome sequences. PLoS ONE 2012, 7, 13. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  27. Cameron, S.L.; Lambkin, C.L.; Barker, S.C.; Whiting, M.F. A mitochondrial genome phylogeny of Diptera: Whole genome sequence data accurately resolve relationships over broad timescales with high precision. Syst. Entomol. 2007, 32, 40–59. [Google Scholar] [CrossRef]
  28. 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]
  29. Xu, D.; Yu, T.; Zhang, Y. 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]
  30. Huang, W.; Zhang, Y. Characterization of two complete mitochondrial genomes of Ledrinae (Hemiptera: Cicadellidae) and phylogenetic analysis. Insects 2020, 11, 609. [Google Scholar] [CrossRef]
  31. Tang, J.; Zhang, Y. Review of the leafhopper genus Hylica Stal (Hemiptera: Cicadellidae: Hylicinae) with description of one new species. Zootaxa 2018, 4388, 526–536. [Google Scholar] [CrossRef]
  32. Tang, J.; Zhang, Y. Review of the leafhopper genus Kalasha Distant (Hemiptera: Cicadellidae: Hylicinae). Zootaxa 2019, 4545, 408–418. [Google Scholar] [CrossRef]
  33. Tang, J.; Zhang, Y. Review of the Oriental leafhopper genus Balala Distant, with new species and new records (Hemiptera: Cicadellidae: Hylicinae). Zootaxa 2020, 4731, 23–42. [Google Scholar] [CrossRef]
  34. Bernt, M.; Donath, A.; Juhling, F.; Externbrink, F.; Florentz, C.; Fritzsch, G.; Putz, J.; Middendorf, M.; Stadler, P.F. MITOS: Improved de novo metazoan mitochondrial genome annotation. Mol. Phylogenet. Evol. 2013, 69, 313–319. [Google Scholar] [CrossRef]
  35. Lowe, T.M.; Eddy, S.R. tRNAscan-SE: A program for improved detection of transfer RNA genes in genomic sequence. Nucleic Acids Res. 1997, 25, 955–964. [Google Scholar] [CrossRef] [PubMed]
  36. Benson, G. Tandem repeats finder: A program to analyze DNA sequences. Nucleic Acids Res. 1999, 27, 573–580. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  37. Grant, J.R.; Stothard, P. The CGView Server: A comparative genomics tool for circular genomes. Nucleic Acids Res. 2008, 36, W181–W184. [Google Scholar] [CrossRef] [PubMed]
  38. Kumar, S.; Stecher, G.; Tamura, K. MEGA7: Molecular Evolutionary Genetics Analysis version 7.0 for bigger datasets. Mol. Biol. Evol. 2016, 33, 1870–1874. [Google Scholar] [CrossRef] [Green Version]
  39. Zhang, D.; Gao, F.L.; Jakovlic, I.; Zou, H.; Zhang, J.; Li, W.X.; Wang, G.T. PhyloSuite: An integrated and scalable desktop platform for streamlined molecular sequence data management and evolutionary phylogenetics studies. Mol. Ecol. Resour. 2020, 20, 348–355. [Google Scholar] [CrossRef]
  40. 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]
  41. Rozas, J.; Ferrer-Mata, A.; Sanchez-DelBarrio, J.C.; Guirao-Rico, S.; Librado, P.; Ramos-Onsins, S.E.; Sanchez-Gracia, A. DnaSP 6: DNA sequence polymorphism analysis of large data sets. Mol. Biol. Evol. 2017, 34, 3299–3302. [Google Scholar] [CrossRef]
  42. Song, N.; Liang, A. The complete mitochondrial genome sequence of Geisha distinctissima (Hemiptera: Flatidae) and comparison with other hemipteran insects. Acta Biochim. Biophys. Sin. 2009, 41, 206–216. [Google Scholar] [CrossRef] [Green Version]
  43. Zhang, K.J.; Zhu, W.C.; Rong, X.; Zhang, Y.K.; Ding, X.L.; Liu, J.; Chen, D.S.; Du, Y.; Hong, X.Y. The complete mitochondrial genomes of two rice planthoppers, Nilaparvata lugens and Laodelphax striatellus: Conserved genome rearrangement in Delphacidae and discovery of new characteristics of atp8 and tRNA genes. BMC Genom. 2013, 14, 12. [Google Scholar] [CrossRef] [Green Version]
  44. Liu, J.; Bu, C.; Wipfler, B.; Liang, A. Comparative analysis of the mitochondrial genomes of Callitettixini spittlebugs (Hemiptera: Cercopidae) confirms the overall high evolutionary speed of the AT-rich region but reveals the presence of short conservative elements at the tribal level. PLoS ONE 2014, 9. [Google Scholar] [CrossRef] [Green Version]
  45. Liu, J.; Liang, A.P. The complete mitochondrial genome of spittlebug Paphnutius ruficeps (Insecta: Hemiptera: Cercopidae) with a fairly short putative control region. Acta Biochim. Biophys. Sin. 2013, 45, 309–319. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  46. Stewart, J.B.; Beckenbach, A.T. Insect mitochondrial genomics: The complete mitochondrial genome sequence of the meadow spittlebug Philaenus spumarius (Hemiptera: Auchenorrhyncha: Cercopoidae). Genome 2005, 48, 46–54. [Google Scholar] [CrossRef] [PubMed]
  47. Lukasik, P.; Chong, R.A.; Nazario, K.; Matsuura, Y.; Bublitz, A.C.; Campbell, M.A.; Meyer, M.C.; Van Leuven, J.T.; Pessacq, P.; Veloso, C.; et al. One hundred mitochondrial genomes of cicadas. J. Hered. 2019, 110, 247–256. [Google Scholar] [CrossRef] [PubMed]
  48. Liang, A.P.; Gao, J.; Zhao, X. Characterization of the complete mitochondrial genome of the treehopper Darthula hardwickii (Hemiptera: Aetalionidae). Mitochondrial DNA Part A 2016, 27, 3291–3292. [Google Scholar] [CrossRef]
  49. Zhao, X.; Liang, A.P. Complete DNA sequence of the mitochondrial genome of the treehopper Leptobelus gazella (Membracoidea: Hemiptera). Mitochondrial DNA Part A 2016, 27, 3318–3319. [Google Scholar] [CrossRef]
  50. Song, N.; Cai, W.Z.; Li, H. Insufficient power of mitogenomic data in resolving the auchenorrhynchan monophyly. Zool. J. Linn. Soc. 2018, 183, 776–790. [Google Scholar] [CrossRef]
  51. Song, N.; Cai, W.; Li, H. Deep-level phylogeny of Cicadomorpha inferred from mitochondrial genomes sequenced by NGS. Sci. Rep. 2017, 7. [Google Scholar] [CrossRef] [Green Version]
  52. Mao, M.; Yang, X.S.; Bennett, G. The complete mitochondrial genome of Entylia carinata (Hemiptera: Membracidae). Mitochondrial DNA Part B Resour. 2016, 1, 662–663. [Google Scholar] [CrossRef] [Green Version]
  53. 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]
  54. 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, 19. [Google Scholar] [CrossRef] [Green Version]
  55. Wu, Y.F.; Dai, R.H.; Zhan, H.P.; Qu, L. Complete mitochondrial genome of Drabescoides nuchalis (Hemiptera: Cicadellidae). Mitochondrial DNA Part A 2016, 27, 3626–3627. [Google Scholar] [CrossRef] [PubMed]
  56. Wang, J.J.; Yang, M.F.; Dai, R.H.; Li, H. Complete mitochondrial genome of Evacanthus heimianus (Hemiptera: Cicadellidae: Evacanthinae) from China. Mitochondrial DNA Part B Resour. 2019, 4, 284–285. [Google Scholar] [CrossRef] [Green Version]
  57. Yuan, Z.W.; Yang, X.; Li, C.; Song, Y.H. The complete mitochondrial genome of the leafhopper Evacanthus acuminatus (Hemiptera: Cicadellidae: Evacanthinae). Mitochondrial DNA Part B Resour. 2019, 4, 3866–3867. [Google Scholar] [CrossRef] [Green Version]
  58. 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]
  59. 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]
  60. Dai, R.H.; Wang, J.J.; Yang, M.F. The complete mitochondrial genome of the leafhopper Idioscopus clypealis (Hemiptera: Cicadellidae: Idiocerinae). Mitochondrial DNA Part B Resour. 2018, 3, 32–33. [Google Scholar] [CrossRef] [Green Version]
  61. 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]
  62. Li, D.-F.; Dai, R.-H. The complete mitochondrial genome of Tituria pyramidata (Hemiptera: Cicadellidae: Ledrinae) from China. Mitochondrial DNA Part B Resour. 2020, 5, 1757–1758. [Google Scholar] [CrossRef] [Green Version]
  63. 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, 11. [Google Scholar] [CrossRef]
  64. He, H.; Yang, M. The mitogenome of Mileewa margheritae (Hemiptera: Cicadellidae: Mileewinae). Mitochondrial DNA Part B Resour. 2020, 5, 3181–3182. [Google Scholar] [CrossRef]
  65. 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]
  66. Wang, J.J.; Dai, R.H.; Li, H.; Zhan, H.P. Characterization of the complete mitochondrial genome of Japanagallia spinosa and Durgades nigropicta (Hemiptera: Cicadellidae: Megophthalminae). Biochem. Syst. Ecol. 2017, 74, 33–41. [Google Scholar] [CrossRef]
  67. Tan, C.; Chen, X.; Li, C.; Song, Y. The complete mitochondrial genome of Empoascanara sipra (Hemiptera: Cicadellidae: Typhlocybinae) with phylogenetic consideration. Mitochondrial DNA Part B Resour. 2019, 5, 260–261. [Google Scholar] [CrossRef] [Green Version]
  68. Song, Y.H.; Yuan, X.W.; Li, C. The mitochondrial genome of Paraahimia luodianensis (Hemiptera: Cicadellidae: Typhlocybinae), a new genus and species from China. Mitochondrial DNA Part B Resour. 2020, 5, 1351–1352. [Google Scholar] [CrossRef]
  69. 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] [Green Version]
  70. Talavera, G.; Castresana, J. Improvement of phylogenies after removing divergent and ambiguously aligned blocks from protein sequence alignments. Syst. Biol. 2007, 56, 564–577. [Google Scholar] [CrossRef] [Green Version]
  71. 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]
  72. 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]
  73. Minh, B.Q.; Nguyen, M.A.T.; von Haeseler, A. Ultrafast approximation for phylogenetic bootstrap. Mol. Biol. Evol. 2013, 30, 1188–1195. [Google Scholar] [CrossRef]
  74. Ronquist, F.; Teslenko, M.; van der Mark, P.; Ayres, D.L.; Darling, A.; Hohna, S.; Larget, B.; Liu, L.; Suchard, M.A.; Huelsenbeck, J.P. MrBayes 3.2: Efficient Bayesian phylogenetic inference and model choice across a large model space. Syst. Biol. 2012, 61, 539–542. [Google Scholar] [CrossRef] [Green Version]
  75. Miller, M.A.; Pfeier, W.; Schwartz, T. Creating the CIPRES science gateway for inference of large phylogenetic trees. In Proceedings of the Gateway Computing Environments Workshop (GCE), New Orleans, LA, USA, 14 November 2010; Institute of Electrical and Electronics Engineers (IEEE): New Orleans, LA, USA, 2010; Volume 14, pp. 1–8. [Google Scholar]
  76. Clary, D.O.; Wolstenholme, D.R. The mitochondrial DNA molecular of Drosophila yakuba: Nucleotide sequence, gene organization, and genetic code. J. Mol. Evol. 1985, 22, 252–271. [Google Scholar] [CrossRef] [PubMed]
  77. Ojala, D.; Montoya, J.; Attardi, G. tRNA punctuation model of RNA processing in human mitochondria. Nature 1981, 290, 470–474. [Google Scholar] [CrossRef] [PubMed]
  78. Du, Y.M.; Zhang, C.N.; Dietrich, C.H.; Zhang, Y.L.; Dai, W. Characterization of the complete mitochondrial genomes of Maiestas dorsalis and Japananus hyalinus (Hemiptera: Cicadellidae) and comparison with other Membracoidea. Sci. Rep. 2017, 7, 10. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  79. Mao, M.; Yang, X.S.; Bennett, G. The complete mitochondrial genome of Macrosteles quadrilineatus (Hemiptera: Cicadellidae). Mitochondrial DNA Part B Resour. 2017, 2, 173–175. [Google Scholar] [CrossRef] [Green Version]
  80. Johnson, K.P.; Dietrich, C.H.; Friedrich, F.; Beutel, R.G.; Wipfler, B.; Peters, R.S.; Allen, J.M.; Petersen, M.; Donath, A.; Walden, K.K.O.; et al. Phylogenomics and the evolution of hemipteroid insects. Proc. Natl. Acad. Sci. USA 2018, 115, 12775–12780. [Google Scholar] [CrossRef] [Green Version]
  81. Skinner, R.K.; Dietrich, C.H.; Walden, K.K.O.; Gordon, E.; Sweet, A.D.; Podsiadlowski, L.; Petersen, M.; Simon, C.; Takiya, D.M.; Johnson, K.P. Phylogenomics of Auchenorrhyncha (Insecta: Hemiptera) using transcriptomes: Examining controversial relationships via degeneracy coding and interrogation of gene conflict. Syst. Entomol. 2020, 45, 85–113. [Google Scholar] [CrossRef]
Insects 11 00869 g001 550
Figure 1. The circular 15.7 kb N. tuberculatus and 14.7 kb H. paradoxa mtgenomes.
Figure 1. The circular 15.7 kb N. tuberculatus and 14.7 kb H. paradoxa mtgenomes.
Insects 11 00869 g001
Insects 11 00869 g002 550
Figure 2. The circular 16.2 kb B. fujiana and 15.7 kb K. nativa mtgenomes.
Figure 2. The circular 16.2 kb B. fujiana and 15.7 kb K. nativa mtgenomes.
Insects 11 00869 g002
Insects 11 00869 g003 550
Figure 3. Relative synonymous codon usage (RSCU) of the mtgenomes of four Hylicinae species.
Figure 3. Relative synonymous codon usage (RSCU) of the mtgenomes of four Hylicinae species.
Insects 11 00869 g003
Insects 11 00869 g004 550
Figure 4. Organization of the control region in Hylicinae mtgenomes. The blue and red ovals indicate the tandem repeats; the purple and red rounds indicate the Poly T/A; the non-repeat regions are shown with green boxes.
Figure 4. Organization of the control region in Hylicinae mtgenomes. The blue and red ovals indicate the tandem repeats; the purple and red rounds indicate the Poly T/A; the non-repeat regions are shown with green boxes.
Insects 11 00869 g004
Insects 11 00869 g005 550
Figure 5. Nucleotide diversities and selection pressures of 13 protein coding genes (PCGs) in Hylicinae. (A) Sliding window analysis of protein-coding genes among four Hylicinae species. The red curve shows the value of Pi (nucleotide diversity). Pi value of each PCG is shown above the arrows. (B) Genetic distances (on average) and ratio of non-synonymous (Ka) to synonymous (Ks) substitution rates of each protein-coding gene among four Hylicinae species.
Figure 5. Nucleotide diversities and selection pressures of 13 protein coding genes (PCGs) in Hylicinae. (A) Sliding window analysis of protein-coding genes among four Hylicinae species. The red curve shows the value of Pi (nucleotide diversity). Pi value of each PCG is shown above the arrows. (B) Genetic distances (on average) and ratio of non-synonymous (Ka) to synonymous (Ks) substitution rates of each protein-coding gene among four Hylicinae species.
Insects 11 00869 g005
Insects 11 00869 g006 550
Figure 6. Phylogenetic tree produced by IQtree/MrBayes based on the dataset of P123R and MrBayes based on the dataset of P123. Numerals at nodes are bootstrap support values (BS) and Bayesian posterior probabilities (PP). “-” indicates the clades or species are different in Bayesian inference (BI).
Figure 6. Phylogenetic tree produced by IQtree/MrBayes based on the dataset of P123R and MrBayes based on the dataset of P123. Numerals at nodes are bootstrap support values (BS) and Bayesian posterior probabilities (PP). “-” indicates the clades or species are different in Bayesian inference (BI).
Insects 11 00869 g006
Insects 11 00869 g007 550
Figure 7. Maximum likelihood tree inferred from the dataset of PCG using IQ-TREE. Node numbers show bootstrap support values.
Figure 7. Maximum likelihood tree inferred from the dataset of PCG using IQ-TREE. Node numbers show bootstrap support values.
Insects 11 00869 g007
Insects 11 00869 g008 550
Figure 8. Phylogenetic tree produced by IQtree and MrBayes based on the dataset of AA. Numerals at nodes are Bayesian posterior probabilities (PP) and bootstrap support values (BS), respectively. “-” indicates the clades or species are different.
Figure 8. Phylogenetic tree produced by IQtree and MrBayes based on the dataset of AA. Numerals at nodes are Bayesian posterior probabilities (PP) and bootstrap support values (BS), respectively. “-” indicates the clades or species are different.
Insects 11 00869 g008
Table
Table 1. Information of samples for phylogenetic analyses.
Table 1. Information of samples for phylogenetic analyses.
SuperfamilyFamily/SubfamilySpeciesAccession NumberReference
Outgroup
FulgoroideaFulgoridae/FulgorinaeFulgora candelariaNC_019576[26]
Flatidae/FlatinaeGeisha distinctissimaNC_012617[42]
Delphacidae/DelphacinaeNilaparvata lugensNC_021748[43]
CercopoideaCercopidae/CercopinaeCosmoscarta dorsimaculaNC_040115Unpublished
Cercopidae/CallitettixinaeCallitettix braconoidesNC_025497[44]
Cercopidae/IschnorhininaePaphnutius ruficepsNC_021100[45]
Aphrophoridae/AphrophorinaePhilaenus spumariusNC_005944[46]
CicadoideaCicadidae/CicadinaeDiceroprocta semicinctaKM000131Unpublished
Cicadidae/CicadettinaeMagicicada tredeculaMH937705[24]
Cicadidae/TibicininaeTettigades auropilosaKM000129Unpublished
Tettigarctidae/TettigarctinaeTettigarcta crinitaMG737758[47]
Ingroup
MembracoideaAetalionidae/AetalioninaeDarthula hardwickiiNC_026699[48]
Membracidae/CentrotinaeLeptobelus gazellaNC_023219[49]
Membracidae/CentrotinaeCentrotus cornutusKX437728[50]
Membracidae/CentrotinaeTricentrus sp.KY039115[51]
Membracidae/CentrotinaeLeptobelus sp.JQ910984[25]
Membracidae/SmiliinaeEntylia carinataNC_033539[52]
Cicadellidae/CicadellinaeCuerna sp.KX437741[50]
Cicadellidae/CicadellinaeGraphocephala sp.KX437740[50]
Cicadellidae/CicadellinaeCicadella viridisKY752061Unpublished
Cicadellidae/CicadellinaeHomalodisca vitripennisNC_006899Unpublished
Cicadellidae/CoelidiinaeTaharana fascianaNC_036015[53]
Cicadellidae/CoelidiinaeOlidiana sp.KY039119Unpublished
Cicadellidae/CoelidiinaeOlidiana_ritcheriinaMK738125[54]
Cicadellidae/DeltocephalinaeHishimonus phycitisKX437727[50]
Cicadellidae/DeltocephalinaeHishimonus sp.KX437735[50]
Cicadellidae/DeltocephalinaePhlogotettix sp. 1KY039135[51]
Cicadellidae/DeltocephalinaePhlogotettix sp. 2KX437721[50]
Cicadellidae/DeltocephalinaeScaphoideus variusKY817245[51]
Cicadellidae/DeltocephalinaeScaphoideus nigrivalveusKY817244[51]
Cicadellidae/DeltocephalinaeAlobaldia tobaeKY039116[51]
Cicadellidae/DeltocephalinaeAthysanopsis sp.KX437726[50]
Cicadellidae/DeltocephalinaeDrabescoides nuchalisNC_028154[55]
Cicadellidae/EvacanthinaeEvacanthus heimianusMG813486[56]
Cicadellidae/EvacanthinaeEvacanthus acuminatusMK948205[57]
Cicadellidae/EvacanthinaeOnukia onukiiMK251119[58]
Cicadellidae/EurymelinaePopulicerus populiNC_039427[59]
Cicadellidae/EurymelinaeIdiocerus laurifoliaeNC_039741[59]
Cicadellidae/EurymelinaeIdioscopus clypealisNC_0396f42[60]
Cicadellidae/EurymelinaeIdioscopus myricaMH492317[59]
Cicadellidae/HylicinaeKalasha nativaMW218662This study
Cicadellidae/HylicinaeNacolus tuberculatusMW218663This study
Cicadellidae/HylicinaeBalala fujianaMW218661This study
Cicadellidae/HylicinaeHylica paradoxaMW218660This study
Cicadellidae/IassinaeTrocnadella arisanaNC036480[61]
Cicadellidae/IassinaeKrisna rufimarginataNC046068[61]
Cicadellidae/IassinaeKrisna concavaMN577635[61]
Cicadellidae/IassinaeBatracomorphus lateprocessusMG813489[61]
Cicadellidae/LedrinaeTituria pyramidataMN920440[62]
Cicadellidae/LedrinaeLedra audituraMK387845[63]
Cicadellidae/LedrinaePetalocephala ochraceaKX437734[50]
Cicadellidae/LedrinaePetalocephala chlorophanaMT610899[30]
Cicadellidae/LedrinaeTituria sagittataMT610900[30]
Cicadellidae/MileewinaeMileewa margheritaeMT483998[64]
Cicadellidae/MileewinaeMileewa pontaMT497465[65]
Cicadellidae/MegophthalminaeJapanagallia spinosaNC_035685[66]
Cicadellidae/MegophthalminaeDurgades nigropictaNC_035684[66]
Cicadellidae/TyphlocybinaeTyphlocyba sp.KY039138[51]
Cicadellidae/TyphlocybinaeEmpoascanara dwalataMT350235Unpublished
Cicadellidae/TyphlocybinaeEmpoascanara sipraMN604278[67]
Cicadellidae/TyphlocybinaeParaahimia luodianensisNC047464[68]
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Tang, J.; Huang, W.; Zhang, Y. The Complete Mitochondrial Genome of Four Hylicinae (Hemiptera: Cicadellidae): Structural Features and Phylogenetic Implications. Insects 2020, 11, 869. https://doi.org/10.3390/insects11120869

AMA Style

Tang J, Huang W, Zhang Y. The Complete Mitochondrial Genome of Four Hylicinae (Hemiptera: Cicadellidae): Structural Features and Phylogenetic Implications. Insects. 2020; 11(12):869. https://doi.org/10.3390/insects11120869

Chicago/Turabian Style

Tang, Jiu, Weijian Huang, and Yalin Zhang. 2020. "The Complete Mitochondrial Genome of Four Hylicinae (Hemiptera: Cicadellidae): Structural Features and Phylogenetic Implications" Insects 11, no. 12: 869. https://doi.org/10.3390/insects11120869

APA Style

Tang, J., Huang, W., & Zhang, Y. (2020). The Complete Mitochondrial Genome of Four Hylicinae (Hemiptera: Cicadellidae): Structural Features and Phylogenetic Implications. Insects, 11(12), 869. https://doi.org/10.3390/insects11120869

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