The complete mitochondrial genomes of two erythroneurine leafhoppers (Hemiptera, Cicadellidae, Typhlocybinae, Erythroneurini) with assessment of the phylogenetic status and relationships of tribes of Typhlocybinae

Abstract The number and classification of tribes in the leafhopper subfamily Typhlocybinae are not yet fully clear, and molecular data has recently been used to help resolve the problem. In this study, the mitochondrial genomes of Mitjaevia shibingensis Chen, Song & Webb, 2020 and M. dworakowskae Chen, Song & Webb, 2020 of the tribe Erythroneurini (Cicadellidae, Typhlocybinae) were sequenced. Most protein-coding genes (PCGs) start with ATN and end with TAA or TAG, and the AT content of these three codons were found differ from previous results that show that the first codon has the highest incidence. Two rRNA genes are highly conserved, and the AT content in 16S is higher than that of 12S. The nucleotide diversity and genetic distance among 13 PCGs of the four tribes from Typhlocybinae show that Empoascini nucleotide diversity is significantly less than in the other three tribes, and have the largest distance from the others, while Typhlocybini and Zyginellini have the smallest distance, indicating that the relationship between the two is the closest. The nad2, nad4, nad4L, and nad5 genes have greater nucleotide diversity, showing potential for use as the main markers for species identification. The phylogenetic analysis yielded a well-supported topology with most branches receiving maximum support and a few branches pertaining to relationships within Zyginellini and Typhlocybini receiving lower support. The species of these two tribes are intertwined, and it was impossible to resolve them into separate branches. In addition, the tribes Empoascini and Erythroneurini were recovered as monophyletic, and Alebrini was placed at the base of the tree as the most primitive. These results are broadly in line with other molecular phylogenetical studies which differ from traditional morphological classification.


Introduction
Cicadellidae (leafhoppers) are the largest family of the order Hemiptera. Representatives are important agricultural and forestry pests that feed on a variety of plants such as cereal crops, vegetables, and fruit trees, and they are also vectors of plant pathogens (Morris 1971;Guo 2011;Roddee et al. 2018). Leafhoppers also have many interesting characteristics, such as varied lifestyles and feeding strategies (including utilizing different endosymbionts), covering their bodies with brochosomes, and producing courtship signals through the plant substrate. These characteristics make leafhoppers suitable material for studies on biological evolution and geographical research (Chen et al. 2020a, b). The subfamily of Typhlocybinae is the second largest group of Cicadellidae and is widely distributed in the six major zoogeographic regions of the world. Erythroneurini, the largest tribe of Typhlocybinae, includes ~2,000 species worldwide (Wang et al. 2012).
The traditional classification of leafhopper has attracted much research attention, including the classification of Typhlocybinae. At present, Typhlocybinae contains six tribes (Alebrini, Empoascini, Erythroneurini, Zyginellini, Typhlocybini, Dikraneurini) but this division remains controversial (Hamilton 1998;Balme 2007;Dietrich 2013). Today, the emergence of next-generation sequencing technology is a breakthrough for solving this problem enabling mitochondrial genomic data to verify and reference the existing family-level classification of Typhlocybinae. Many previous attempts have been made to estimate phylogenetic relationships among leafhoppers mostly by using either morphological data or sequence data from a few gene regions (Hamilton 1983;Dietrich et al. 2001;Zahniser and Dietrich 2013;Krishnankutty et al. 2016), but there is very little research on Typhlocybinae. So far, in the National Center for Biotechnology Information (NCBI), Typhlocybinae only has complete mitochondrial genomic data for 19 species (Table 1).
The insect mitochondrial genome (mtDNA) is usually a closed double-stranded DNA molecule with a molecular weight of 14-20 kb. Usually, it contains 37 genes, including 13 protein-coding genes (PCGs), NADH dehydrogenase 1-6 and 4L (nad1-6 and nad4L), cytochrome c oxidase subunits 1-3 (cox1-3), ATPase subunit 6 and 8 (atp6 and atp8), cytochrome b (cytb), two ribosomal RNAs genes (16S and 12S) and 22 transfer RNA (tRNA) genes. A region rich in A + T, the control region, is also present (Boore 1999; Wang et al. 2018). Compared with the nuclear genome, the in-sect mitochondrial genome has the characteristics of a simple structure, low molecular weight, stable composition, conservative arrangement, maternal inheritance, and easy detection. It is very suitable for the study of evolutionary genomics and is widely used to identify the phylogenetic relationships and population structures at different taxonomic levels (Saccone et al. 2002;Cook et al. 2005;Wilson et al. 2008).
To further enrich the mitochondrial genome data of leafhoppers and provide comparative data for closely related species, we sequenced and analyzed the complete mitochondrial genomes of Mitjaevia shibingensis and M. dworakowskae and analyzed their phylogenetic relationship with other Typhlocybinae. The new molecular data obtained will help in the identification of leafhopper species, kinship comparison, and future studies on population genetics and evolution.

Mitogenome sequencing, assembly, and annotation
For this study, samples of Mitjaevia shibingensis and M. dworakowskae were collected in 100% alcohol and stored at -20 °C in the laboratory. Total DNA was extracted from the entire body without the abdomen and wings. The mitochondrial gene sequences were obtained through second-generation sequencing. Primers were designed to amplify the mtDNA sequence in PCR reactions. The PCR reaction was performed using the LA Taq polymerase. The PCR conditions were as follows: initial denaturation at 94 °C for 2 min, then 35 cycles of denaturation at 94 °C for 30 sec, annealing at 55 °C for 30 sec, and extension at 72 °C for 1 min/kb, followed by the final extension at 72 °C for 10 min. The PCR products were sequenced directly, or, if needed, first cloned into a pMD18-T vector (Takara, JAP) and then sequenced, by the dideoxynucleotide procedure, using an ABI 3730 automatic sequencer (Sanger sequencing) using the same set of primers. After quality-proofing of the obtained fragments, the complete mt genome sequence was assembled manually using DNAStar (Burland 2000), and a homology search was performed by the Blast function in NCBI to verify the amplified sequence as the target sequence (Meng et al. 2013;Yu et al. 2017). The nucleotide base composition, codon usage, and A + T content values were analyzed with MEGA 6.06 (Tamura et al. 2013). The secondary structure of tRNA genes was annotated using online tools tRNAscan-SE 1.21 and ARWEN (Lowe and Eddy 1997;Laslett and Canbäck 2008). The tandem repeat sequence in the control area was determined by the online search tool Tandem Repeats Finder (Benson 1999 (Rozas et al. 2017). And the genetic distance of the four tribes was estimated in MEGA 6.06.

Phylogenetic analysis
The phylogenetic analysis included two sets of data. First, the phylogenetic tree was constructed based on 29 cox1 data among six tribes of Typhlocybiane and two outgroups. Secondly, phylogenetic tree analysis was conducted using a dataset including the complete mitochondrial genomes of the two newly sequenced erythroneurine species, 17 typhlocybiane species, and two outgroups, of which nine sets of data were from team sequencing, while the remaining 10 were obtained from the NCBI database (Table 1). The Gblocks Server online platform was used to eliminate poorly aligned positions and divergent regions of DNA protein alignment, and all alignments were checked and corrected in MEGA 6.06 prior to the phylogenetic analysis (Tamura et al. 2013). Five datasets were generated: (1) cox1 with 573 nucleotides (2) 13 PCGs with 10,452 nucleotides (PCGs); (3) the first and second codon positions of the 13 PCGs with 6968 nucleotides (PCG12); (4) 13 PCGs with 10,452 nucleotides and 2 rRNA with 1615 nucleotides (PCGR); (5) and amino acid sequences of the 13 PCGs with 2666 amino acids (PCGAA).
The trimmed datasets were used to estimate the phylogeny by maximum likelihood (ML) using IQ-TREE and Bayesian inference (BI) using MrBayes 3.2.7 (Zhou et al. 2011;. ML constructed with the IQ-TREE used an ultrafast bootstrap approximation approach with 10,000 replicates and calculated bootstrap scores for each node (BP). BI selected GTR + I + G as the optimal model, running 10 million generations twice, sampling once every 1000 generations, after the average standard deviation of the segmentation frequency drops below 0.01, with the first 25% of the samples are discarded burn-in, and the remaining trees used to generate a consensus tree and calculate the posterior probability (PP) of each branch.

Organization and composition of the genome
The genomic organization and nucleotide composition of the two new mitogenomes sequenced in this study are similar to those of other previously reported Typhlocybina (Tan et al. 2020;Yuan et al. 2020). The complete mitogenomes of M. shibingensis and M. dworakowskae are double-stranded plasmids with 15,788 and 16,399 bp, respectively. Both species contain the usual 13 PCGs, 22 tRNA genes, 2 rRNA genes, and a control region (Fig. 1). Fourteen genes encode in the minority strand (L-strand) while the others encode in the majority strand (H-strand). Mitjaevia shibingensis has a total of 46 bp intergenic space in 12 regions ranging from 1 to 8 bp. Eleven genes were found to overlap by a total of 47 bp. Mitjaevia dworakowskae has a total of 84 bp intergenic space in 12 regions ranging from 2 to 15 bp, and seven genes were found to overlap by a total of 23 bp ( Table 2).
The AT contents and skew statistics are shown in Table 3. The mitochondrial genomes of M. shibingensis and M. dworakowskae show heavy AT nucleotide bias, with A + T% content for the whole sequence was 78.4% and 79.0%, respectively. Similar patterns of nucleotide composition are also found in other leafhopper species Wang et al. 2018;Xian et al. 2020). The control region (CR) has the strongest A + T% bias, while the PCGs shows the lowest A + T% among whole genes. The whole genome has positive AT skews (0.042, 0.051) and negative GC skews (-0.074, -0.104). Analysis of 37 individual genes of the two species show that AT skews are mostly positive, while for GC skews, the genes of both species are mostly negative (Fig. 2). Positive AT skews indicates that the content of base A is higher than that of base T. However, in a few genes, although the AT skews is negative, the difference in absolute value was very small. For negative GC skew, a negative value indicates that the content of base G is lower than that of base C, while a positive value indicates the opposite. In general, the basic composition of these two species is biased towards A and C.    while nad5 and atp8 genes have TTG, a pattern also observed in other leafhopper mitogenomes (Wang et al. 2017;Wang et al. 2020). Conventional stop codons (TAA or TAG) appear in most PCGs, except that cox2 and nad1 use an incomplete codon (a single T--) as the stop codon (see Tables 2, 3).
Research determined the behavior of the PCG codon families and found that codon usage was very similar among Cicadellidae mitogenomes when the results of two species were calculated and summarized (see Table 4, Fig. 3). All 62 available codons (excluding TAA and TAG) are present in M. shibingensis and M. dworakowskae. Synonymous codon usage bias was observed in both mitochondrial genomes, and 22 codons were used more frequently than other codons. The four most abundant codons were AAU (Asn), AAA (Lys), AUU (Ile), and UUA (Leu2). The preferred codons all end with A or U, which contribute to the high A + T bias of the entire mitogenomes.  Transfer RNA and ribosomal RNA genes All 22 typical tRNA genes are present in the M. shibingensis and M. dworakowskae mitochondrial genomes, of which 14 genes were oriented on the major strand (H-strand), whereas the others were transcribed on the minor strand (L-strand). Their nucleotide lengths are almost identical between species, ranging from 60 bp to 71 bp ( Table 2). The average AT content values of tRNAs is 78.5% and 79% in each species, respectively, and the tRNA genes have negligible AT and GC skews (Table 3). Compared to the ancestral insect mitochondrial gene order, no tRNA gene rearrangements were found. All of the tRNA genes can be folded into typical clover-leaf secondary structures except for the trnS1 in both species' mitochondrial genomes, which lacks the dihydrouridine (DHU) stem and forms a simple loop. In the metazoan mt genome, lack of the DHU arm was very common in trnS1 ). Based on the secondary structure, a total of 24 and 22 G-U weak base pairs were found in tRNAs of M. shibingensis and M. dworakowskae, respectively (Figs 4, 5). Most mismatched nucleotides were G-U pairs, which form weak bonds in tRNA and non-classical pairs in tRNA secondary structure, similar to other Cicadellidae (Jia et al. 2010). Leafhopper ribosomal RNA (rRNA) includes 16S RNA and 12S RNA. These two genes are highly conserved and are encoded on the minor strand (L-strand). Similar to other known insects, the content of A + T% in 16S was higher than that of 12S.  The 16s genes of M. shibingensis and M. dworakowskae were 1192 bp and 1852 bp in length, with AT contents of 82.90% and 82.50%, respectively, and located between trnL2 and trnV. The 12S rRNA genes of both were 732 bp and 736 bp in length, with AT contents of 81.80% and 81.90%, respectively, and located after trnV. The rRNA genes showed a positive AT skew and negative GC skew (Table 3).

Control region
Like the typical insect mitochondrial genome, the mt genomes of M. shibingensis and M. dworakowskae have a large non-coding region, which was identified as the control region and located downstream of 12S. Control regions of both species were rich in AT, with lengths of 1429 bp and 2048 bp AT contents of 92.7% and 88.9%, respectively (Table 3). The control regions in the three available Mitjaevia mitogenomes were variable and not highly conserved, and their lengths ranged from 15 and 784 bp with variable numbers of repeat sequences (Fig. 6). Mitjaevia shibingensis included 21 types of repeat unit (R), 17 kinds of repeats (R1, R2) were found in M. dworakowskae with various lengths and copy numbers, and three repeat units were present in M. protuberanta. At present, we were unable to find any correlation in repeating units in the different species, probably because of the limited number of species analyzed in this study. Further comparative studies of additional leafhopper mitogenomes are needed.

Nucleotide diversity and genetic distance analysis
The sliding window analysis shows highly variable nucleotide diversity (Pi values) among 13 PCGs sequences of the four tribes of Typhlocybinae (Fig. 7). Empoascini nucleotide diversity is significantly lower than of the other three tribes. The genes nad2, nad4, nad4L, and nad5 had higher nucleotide diversity, while the genes cox3, cox2, cytb, and cox1 had comparatively low nucleotide diversity when using MEGA 6.06 software, based on Kimura-2 Parameter, and Bootstrap resampling 1000 times to test and analyze the genetic distance of the four tribes of Typhlocybinae. The results show that the genetic distance between Empoascini and the other three tribes is the largest, and between Typhlocybini and Zyginellini are the smallest (Table 5).
Nucleotide diversity analysis, a primary method for identifying the regions with large nucleotide divergence, is especially useful for designing species-specific markers.  These are useful for taxa with highly variable morphological characteristics, especially Typhlocybinae species which belong to groups that are difficult to distinguish by morphology alone (Jia et al. 2010;Ma et al. 2019). Among the four tribes, nad2, nad4, and other highly variable genes have garnered our attention. Whether they can be used as the main marker for species identification or the main related genes that control the appearance of the subfamily is worthy of further study. The genetic distance reflects the distance of the genetic relationships between each tribe. Among the four tribes, Typhlocybini and Zyginellini have the smallest genetic distance, indicating that the relationship between the two is the closest, which is consistent with the results of morphological studies (Zhang 1990;Huang 2013).

Historical review
Typhlocybinae has been divided into tribes based mainly on the characteristics of the wing veins for the past 90 years. Melichar (1903), Distant (1908Distant ( , 1918, and Matsumura (1931), among others, divided Typhlocybinae into Empoascaria and Typhlocybaria according to whether the hindwing apical cell is closed (Fig. 8A). McAtee (1934) further divided Typhlocybinae into four tribes, Alebrini, Dikraneurini, Jorumini, and Eupterygini (Fig. 8B), based on the wing veins. Young (1952), also using the male genitalia, recognised the tribes Alebrini, Dikraneurini, Erythroneurini, and Typhlocybini (Fig. 8C). According to whether the peripheral vein of the hind wings exceeds the end of the R+M vein, it was believed that Erythroneurini evolved from Dikraneurini but Mahmood (1967), when adding the tribe Bakerini, believed this  condition to be an acquired mutation. At the same time, Mahmood (1967) postulated that Erythroneurini might be more closely related to Typhlocybini, but the relationship between Erythroneurini and the other tribes still need to be determined by studying a large number of specimens. Mahmood and Ahmed (1968a), using the characteristic of the peripheral vein of the hind wings extending to the end of the R vein, separated Empoascini from the former Typhlocybini, and Typhlocybinae was divided into six tribes (Fig. 8D). Dworakowska (1970) compared the characteristics of Bakerini, Typhlocybini, and Erythroneurini, and postulated that Bakerini may be a relatively primitive branch of Erythroneurini (Fig. 8E). Dworakowska (1977Dworakowska ( , 1990 discussed Typhlocybini with respect to other tribes, and postulated that the different connection modes of the hindwing peripheral vein and CuA represented different branches, and divided Zyginellini from Typhlocybini resulting in six tribes: Alebrini, Dikraneurini, Empoascini, Erythroneurini, Typhlocybini, and Zyginellin, while Zhang (1990) postulated that Erythroneurini evolved from Empoascini. Since then, Typhlocybinae-related research has followed Dworakowska's six-tribe classification system (Fig. 8F). However, Dietrich (2013) found that the hind wings of leafhoppers have both Zyginellini and Typhlocybini hindwing characteristics when studying the leafhoppers in South America and postulated that the venation characters may not be a stable feature for classification. In terms of overall morphology, Zyginellini and Typhlocybini have similarities present in certain genera and species. Therefore, there is no clear and strong evidence at present to determine whether or not Zyginellini belongs to a natural monophyletic group, and its taxonomic status needs to be further clarified, probably with molecular data.

More recent studies
In recent years, molecular sequencing technology has been widely used in phylogenetic analysis, which can test and verify the results of different levels of more morphology based traditional classifications. Within Typhlocybinae, only a few studies have used the combination of morphological characteristics and molecular data to construct phylogenetic relationships. The amount of data is sparse at present and further data is needed. Dietrich and Dmitriev (2006) used PAUP 4.0b10 to analyze the phylogeny of Typhlocybinae for the first time based on morphological characteristics and concluded that Erythroneurini and Dikraneurini are closely related. However, their analysed samples came mainly from the New World, and whether their results represent the relationship between the tribes of Typhlocybinae remains to be clarified. Balme (2007) combined morphological characteristics with molecular characteristics (16S rDNA, H3) to perform a phylogenetic analysis of Typhlocybinae, and obtained the following topological structure: Alebrini + ((Empoascini + Typhlocybini) + Dikraneurini), but due to the small sample, the results need to be verified. Jiang (2016) used 28S rDNA D2-D3, 16S rDNA sequence and morphological characteristics to make a preliminary exploration of the phylogenetic relationships of Typhlocybinae and obtained the following topological structure: Alebrini + (Empoascini + Erythroneurini), implying that Erythroneurini is a monophyletic group. Song et al. (2020) constructed a phylogenetic tree of Typhlocybinae using 13PCGs of eight species and obtained the following topology: Empoascini+ (Typhlocybini+ (Erythroneurini+Zyginellini)). These four tribes are all monophyletic, and Erythroneurini and Zyginellini are sister groups, differing only slightly from the traditional morphological classification. Jiao (2017) used MP and NJ methods to analyze the phylogenetic relationship between Alebrini and Dikraneurini based only on morphological data. The results showed that the two tribes' monophyly was well supported, and its position in the evolutionary tree was similar to that of Zhang (1990) and showed that Alebrini is more primitive.

Results
This study, based on 29 species of cox1, 19 species of 13 PCGs, and two rRNA mitochondrial genes data of Typhlocybinae produced a slightly different result to the traditional classification with respect to Typhlocybini and Zyginellini. Maximum Likelihood (ML) method was used with IQ-TREE using an ultrafast bootstrap approximation approach with 10,000 replicates. The Bayesian Inference (BI) analysis was performed using MrBayes 3.2.7, with the best fit model GTR+I+ G (Vogler and DeSalle 1993).
Cox1 is one of the mitochondrial protein-coding genes and its bi-terminal sequence is more conservative than cox2. It has a rapid evolution rate and large differences between species, and can provide rich phylogenetic information, hence is an ideal mito-chondrial molecular marker. The gene sequences were obtained in the current study by downloading the cox1 gene sequence of 29 species of Typhlocybinae and two outgroups of Idiocerinae from NCBI to construct a phylogenetic tree. BI and ML analyses generated the same tree topology: (Alebrini + Empoascini) + (Erythroneurini + ((Zyginellini + Typhlocybini) + Dikraneurini))). Most relationships were highly supported, and a few branches pertaining to relationships within Zyginellini and Typhlocybini received lower support (Fig. 9). Also, the tree topology is different from previous research, with Alebrini + Empoascini forming sister groups, and the species of Zyginellini and Typhlocybini are interconnected and cannot be resolved into separate branches. The remaining tribes are monophyletic groups. Alebrini is placed at the base of the tree and is therefore the most primitive. The phylogenetic relationship is generally consistent with the results of previous studies based on morphology and molecules (Balme 2007;Jiang 2016).
At present, the complete mitochondrial genome data of Alebrini and Dikraneurini have not been added to NCBI. Thus, the phylogenetic relationships were analyzed based on the concatenated nucleotide sequences of 13 PCGs and two rRNA from 19 Typhlocybinae (the remaining four tribes) species and two outgroups. Although ML and PB analyses produced inconsistent topologies across the different datasets and models, most relationships were highly supported and consistent in the analyses, and the main difference is the relationship of species between Zyginellini and Typhlocybini. (Figs 10, 11). In this study, Empoascini and Erythroneurini were recovered as monophyletic, always forming a clade with high support values, while Zyginellini and Typhlocybini formed a single branch in every tree, and neither tribe was ever resolved  as monophyletic, which suggests that the hind wing character traditionally used to separate these two tribes is not reliable and that the tribes should probably be treated as synonyms, as was suggested previously by Dietrich (2013). Within the Typhlocybinae, the four species of Empoascini studied constituted one clade and tended to be placed at the basal position of the tree as the sister group to the other tribes. Unlike the previous analyses, the our results support the combination of Zyginellini and Typhlocybini as a tribe. As with other recent studies, our results indicate that sequence data from leafhopper mitogenomes is informative of phylogenetic relationships in the taxonomic hierarchy of this group. Also, the results of the phylogenetic tree and nucleotide diversity are consistent. Empoascini has the lowest nucleotide diversity and is clearly distinguished from the other three tribes. Therefore, we speculate that the richness of nucleotide diversity has an impact on the phylogenetic relationship of Typhlocybinae. However, data are available for only a tiny fraction of species so the addition of more species, and representatives of other major lineages, will be needed to determine the extent to which mitogenome sequence data can resolve leafhopper phylogeny.

Conclusions
This paper describes the complete mitochondrial genomes of M. shibingensis and M. dworakowskae, analyzes the basic composition, location, secondary structure, and other characteristics of PCGs, tRNA genes, rRNA genes, and control regions, and compares them to other Typhlocybinae mitochondrial genomes. The mitogenomes of these two species closely resemble those of most other sequenced leafhoppers in various structural and compositional aspects. The sliding window analysis shows a highly variable nucleotide diversity (Pi values) among 13 PCGs sequences of the four tribes of Typhlocybinae. Empoascini nucleotide diversity is significantly lower than in the other three tribes, and the other three tribes have little difference between them. The genes nad2, nad4, nad4L, and nad5 have higher nucleotide diversity, and whether they can be used as the main markers for species identification or the main related genes that control the appearance of the subfamily is worthy of further study. The genetic distance of the four tribes of Typhlocybinae shows that the Empoascini and the other three tribes are the largest while Typhlocybini and Zyginellini are the smallest and indicates that the relationship between the two is the closest, which is consistent with the results of morphological studies. Phylogenetic analysis of 31 cox1 yielded a well-supported topology with most branches receiving maximum support and a few branches pertaining to relationships within Zyginellini and Typhlocybini receiving lower support; the species of these two tribes are intertwined and cannot be resolved into separate branches, and Alebrini is placed at the base of the tree as the most primitive. Phylogenetic relationships were analyzed based on the concatenated nucleotide sequences of 13 PCGs and two rRNA show that although ML and PB analyses produced inconsistent topologies across the different datasets and models, and most relationships were highly supported and constant in the analyses.
In this study, Empoascini and Erythroneurini were recovered as monophyletic while Zyginellini and Typhlocybini gathered into a single branch and Empoascini tended to be placed at the basal position of the tree as the sister group to the other tribes. This study indicated that mitochondrial genome sequences are informative for leafhopper phylogeny, but unlike the previous analysis (Zhang 1990), the results of this study relocated the taxonomic status and phylogenetic relationship of the six tribes of Typhlocybinae and supported the combination of Zyginellini and Typhlocybini as a single tribe. Also, the results of the phylogenetic tree and nucleotide diversity are consistent. Empoascini has the lowest nucleotide diversity and is clearly distinguished from the other three tribes. Therefore, we speculate that the richness of nucleotide diversity has an impact on the phylogenetic relationship of Typhlocybinae.
Based on the current and previous studies, the classification of the tribes of Typhlocybinae is not yet fully resolved with respect to Typhlocybini and Zyginellini, i.e., one or two tribes. From a molecular perspective, more sequencing data is needed to build a more complete phylogenetic tree to support or modify the traditional morphological classification. To this aim, it is hoped that the new data provided here will facilitate future comparative studies of leafhopper mitogenomes and demonstrate the need for more comparative data.