Complete mitochondrial genome of Neotoxoptera formosana (Takahashi, 1921) (Hemiptera: Aphididae), with the phylogenetic analysis

Abstract In this study, we sequenced and annotated the complete mitochondrial genome (mitogenome) of Neotoxoptera formosana (Takahashi) (Hemiptera: Aphididae). The complete mitogenome of N. formosana is 15,642 bp in length, and includes 13 protein-coding genes (PCGs), 2 ribosome RNA (rRNA) genes, 22 transfer RNA (tRNA) genes, and one control region. The overall base composition was as follows: 45.2% of A, 5.8% of G, 10.5% of C, and 38.4% of T, with a total of A + T content of 83.6%. The phylogenetic tree showed that N. formosana and Myzus persicae were clustered into one branch. This result will enrich the mitogenome of family Aphididae.

Onion aphid, Neotoxoptera formosana (Takahashi) (Hemiptera: Aphididae), originally occurs in Taiwan, China, and is reported in countries all over the world, such as Japan, Brazil, United Kingdom, Italy, Germany and Netherlands (Melo et al. 2005;Pellizzari and Mont a 2005;Piron 2010). Onion aphid is an oligophagous aphid pest of Allium crops which rapidly builds up large colonies and damages the leaves (Hori 2007;Piron 2010). In the last year, we found that N. formosana is a main pest in Allium tuberosum in Guizhou Province, China (Wang et al. 2020). Previous studies mainly focused on ecological characteristic, taxonomic categories and partial sequences of N. formosana. In this study, we determined the mitogenome sequence of N. formosana to accelerate the population genetics and evolution of the genus Neotoxoptera.
The complete mitogenome of N. formosana was a circular double-stranded and 15,642 bp in length (GeneBank accession number MW534268), containing 37 encoding genes (13 PCGs, 22 tRNA genes, and 2 rRNA genes) and one control region. The whole base composition was as follows: 45.2% of A, 5.8% of G, 10.5% of C, and 38.4% of T, with a total of A þ T content of 83.6%. All PCGs began with the typical ATN codons (ATT for COX1, ATP8, ATP6, ND3, ND5, ND6, and ND1; ATG for COX3, ND4, ND4L, and CYTB; ATA for ND2 and COX2). All protein coding genes used TAA as stop codon. The 16 s rRNA gene was 1,263 bp in size and was located between trnL and trnV, while the 12 s rRNA gene was 812 bp in size and was located behind trnV.
The nucleotide sequences of 13 PCGs (delete the third codon position) and 2 rRNA genes from 15 Aphidoidea species and one outgroup taxa from Pyrrhocoroidea species were aligned using MAFFT v7.394 with L-INS-I algorithm (Katoh and Standley 2013). The poorly aligned results were removed by trimAl v1.4.1 (Capella-Guti errez et al. 2009). A phylogenetic tree was constructed for 16 species by maximum likelihood method using IQ-TREE v1.6.3 software under the GTR þ I þ G model (Nguyen et al. 2015). The result confirmed that N. formosana and Myzus persicae were clustered into one branch (Figure 1).

Date availability statement
The data that support the findings of this study are openly available in [NCBI]

Funding
We thank Professor Gexiao Qiao (Chinese Academy of Sciences) for the identification of Neotoxoptera formosana to species. We are also grateful to Chang Han, Likun Zhong and Hongli He (Guizhou University) for technical support. This study was supported by Anshun City Science and Technology Plan Project ( (2020)08