The complete mitochondrial genome of Aphis gossypii Glover, 1877 (Hemiptera: Aphididae) isolated from Plantago asiatica in Korea

Abstract We have determined mitochondrial genome of A. gossypii isolated from Plantago asiatica in Korea. The circular mitogenome of A. gossypii is 16,045 bp including 13 protein-coding genes, two ribosomal RNA genes, 22 transfer RNAs, and a single control region of 798 bp. Its AT ratio is 83.8%. In comparison this mitogenome to Chinese and Korean A. gossypii mitogenomes, 66 single nucleotide polymorphisms (SNPs) and 176 insertions and deletions (INDELs) and 11 SNPs and 173 INDELs are identified, respectively, presenting similar level to those of Nilaparvata lugens, Laodelphax striatellus, and Spodoptera frugiperda and lower than that of Chilo suppresallis.

Aphis gossypii Glover, 1877 is widely distributed polyphagous and popular pest species to agriculture and horticultural species with influencing mortality of crops and transmitting virus (Ebert and Cartwright 1997). It results that biological control programs for managing Aphid population has been utilized in Europe, Russia, and Korea (Gilkeson and Klein 1981;Vuong et al. 2001). More than 200 species including Plantago asiatica (Inaizumi 1970) have been identified as host plants of A. gossypii (CABI 2014).
A. gossypii mitogenome (GenBank accession is MT430940) is 16,045 bp long, which is the longer than those of Chinese and Korean A. gossypii (Zhang et al. 2016; containing 13 protein-coding genes, two rRNAs, and 22 tRNAs. Its nucleotide composition is AT-biased (A þ T is 83.8%). Control region of 784 bp, which is also longer than two A. gossypii mitogenomes, is found.
Sixty-six single nucleotide polymorphisms (SNPs) and 176 insertions and deletions (INDELs) and 11 SNPs and 173 INDELs are identified against those of Chinese and Korean A. gossypii, respectively. These numbers of sequence variations are similar to those of Nilaparvata lugens Choi et al. 2020), Laodelphax striatellus , and Spodoptera frugiperda (Seo, Lee, et al. 2019). It is smaller than that of Chilo suppresallis . One of the possible reasons why our mitogenome is much different from the remaining two mitogenomes, specifically control region, can be a different host plant like the cases of Acyrthosiphon pisum (Peccoud et al. 2009) and Aphis glycine (Park et al., in preparation).
We inferred the phylogenetic relationship of 23 mitogenomes, including three A. gossypii mitogenomes, with one outgroup species, Bemisia tabaci (Tay et al. 2016). Multiple sequence alignment was conducted by MAFFT 7.450 (Katoh and Standley 2013). Bootstrapped maximum likelihood, neighbor joining, and Bayesian Inference trees were constructed using MEGA X (Kumar et al. 2018) based on multiple alignment of mitogenomes. The result shows that Aphis genus was clearly clustered with the rest genera in Aphidinae, but conspicuously, three phylogenetic trees showed incongruent topology of Aphidinae clade (Figure 1). Taken together, our mitogenome is helpful to understand intraspecific mitogenome variations of A. gossypii as well as phylogenetic relationship among Aphidinae species.

Disclosure statement
No potential conflict of interest was reported by the author(s).

Funding
This study was carried out with the support of InfoBoss Research Grant . This stduy was also supported by Research of Animal and Plant Quarantine Agency of South Korea Fund (PQ20180B008-SP).

Data availability statement
The mitochondrial genomes in this study can be accessed via the NCBI GenBank accession number, MT430940. Phylogenetic tree was drawn based on maximum likelihood tree. The numbers above branches indicate bootstrap support values of neighbor joining, maximum likelihood, and Bayesian inference phylogenetic trees, respectively.