853
Views
8
CrossRef citations to date
0
Altmetric
Mitogenome Announcement

The complete mitochondrial genome of Ricania speculum (Walker, 1851) (Hemiptera: Ricaniidae): investigation of intraspecific variations on mitochondrial genome

, , , , , ORCID Icon & show all
Pages 3796-3798 | Received 31 Aug 2020, Accepted 11 Oct 2020, Published online: 20 Nov 2020

Abstract

We have determined a mitochondrial genome of Ricania speculum (Walker, 1851) collected in Jeollabuk-do, Republic of Korea. The circular mitogenome of R. speculum is 15,530 bp long which is shorter than that of the previous mitogenome of R. speculum by 199 bp. It includes 13 protein-coding genes, two ribosomal RNA genes, and 22 transfer RNAs. Intraspecific variation between two mitogenome of R. speculum was investigated: 171 SNPs and 18 INDELs were identified, presenting a high level of intraspecific variations on mitochondrial genome.

Ricania speculum (Walker, 1851) known as the black planthopper is an agricultural pest species in Asia, which has recently been introduced to Europe (Mazza et al. Citation2014). It is an polyphagous species that can feed on more than 60 different woody plant species belonging to 33 families (Mazza et al. Citation2014; Rossi and Lucchi Citation2015; Rossi et al. Citation2015) and currently distribute in China, Indonesia, Japan, Korea, Philippines, Taiwan, Vietnam, and Italy (Wilson et al. Citation2016). It is known that Ricania japonica contains yeast-like fungal symbiont which recycles uric acid secreted by the host with its own uricase (urate oxidase; EC1.7.3.3; Hongoh and Ishikawa Citation2000). Recent study of fungal mitochondrial genomes rescued from NGS raw reads of the R. speculum sample presented that it belonged to the family Ophicocordycipitaceae (Park, Xi, Park, Lee Citation2020). It indicates that NGS sequencing is an efficient method in detecting symbionts like previous studies (Bae et al. Citation2020; Park, Xi, Park, Nam, et al. Citation2020).

To investigate the intraspecific variations of R. speculum mitochondrial genome, we analyzed a mitogenome of R. speculum collected in Jungtosa Temple in Jeollabuk-do in Korea (35.6585365 N, 126.887853E; the specimen is stored in Gyeongsang National University, Korea, Accession number: 2018-HJ-56). DNA was extracted using DNeasy Blood & Tissue Kit (QIAGEN, Hilden, Germany). Sequencing library was constructed using Illumina TruSeq Nano DNA Library Preparation Kit (Illumina, San Diego, CA) following manufacturer’s recommendations with around 350-bp DNA fragments. Raw sequences obtained using Illumina HiSeqX (Macrogen, Korea) were filtered by Trimmomatic 0.33 (Bolger et al. Citation2014) and de novo assembled by Velvet 1.2.10 (Zerbino and Birney Citation2008), SOAPGapCloser 1.12 (Zhao et al. Citation2011), BWA 0.7.17 (Li Citation2013), and SAMtools 1.9 (Li et al. Citation2009) under the environment of Genome Information System (GeIS; http://geis.infoboss.co.kr/; Park et al., in preparation). Geneious R11 11.1.5 (Biomatters Ltd, Auckland, New Zealand) was used to annotate based on the previous mitogenome of R. speculum (NC_031369; Zhang et al. Citation2016).

The new R. speculum mitogenome (GenBank accession is MT834932) collected in Korea is 15,530 bp long, shorter than the previous R. speculum collected in China by 199 bp. The order and contents of 37 genes, 13 protein-coding genes (PCGs), 2 ribosomal RNAs (rRNAs), 22 transfer RNAs (tRNAs), were well conserved in R. speculum as in other hemipteran mitochondrial genomes.

Based on comparison of two R. speculum mitogenomes, we identified 171 SNPs and 18 INDELs of which total length is 233 bp. A number of these identified intraspecific variations were in the large side for insect mitochondrial genomes (e.g. 11–66 SNPs and 173–176 INDELs in two Aphis gossypii mitogenomes; Bae et al. Citation2020; Park et al. Citation2019), although these numbers may not be suitable for direct comparison as single mutations can cause different bases of variations. Eighty-six SNPs were in PCGs, 6 were in rRNAs, 5 were in tRNAs, and the rest were in intergenic spaces, especially the A–T rich control region. Twenty-five of the 86 SNPs in PCGs were non-synonymous while 52 were synonymous, which is similar to the ratio found in Hipparchia autonoe (Lepidoptera), 5:11, showing the total number of SNPs is irrelevant to the ratio of the SNPs (Lee et al. Citation2020). The remaining 9 SNPs were found in 3′ end of nad2 and were unable to be determined either synonymous or not as preceding INDELs caused frame shifts.

The largest 175-bp INDEL was in the A-T rich control region created due to repeat region reduction against the previous mitogenome (NC_031369). Two and three INDELs were found near the 3′ ends of nad1 and nad2, respectively, causing only minor frame shifts because they were only a few amino acids far from termination. Two INDELs, 2 bp and 1 bp long, found in nad5, however, were located near the 5′ end, though they also did not cause much mutations as frame shift caused by the first INDEL was recovered by the second one.

We constructed a Bayesian inference tree using MrBayes 3.2.6 (Huelsenbeck and Ronquist Citation2001) based on multiple sequence alignments of 16 planthopper mitogenomes, excluding the highly variable control regions, together with those of a cicada, Magicicada tredecula (NC_041656), and a frog hopper, Philaenus spumarius (NC_005944), as an outgroup. These sequences were aligned by MAFFT v7.450 (Katoh and Standley Citation2013). The tree was congruent with previous phylogenetic studies as family Ricaniidae was grouped with Flatidae, and then with Fulgoridae (Song and Liang Citation2013; ). All branches were supported with maximum posterior probability values ().

Figure 1. Bayesian inference (1,000,000 generations) phylogenetic tree of 16 Fulgoromorpha mitochondrial genomes: Ricania speculum (this study, MT834932), Ricania speculum (NC_031369), Ricania marginalis (NC_019597), Ricania sp. (MK251070), Geisha distinctissima (NC_012617), Aphaena amabilis (NC_045075), Lycorma delicatula (NC_012835), Sogatella furcifera (NC_021417), Sogatella vibix (NC_042180), Laodelphax striatellus (NC_013706), Changeondelphax velitchkovskyi (NC_037181), Nilaparvata bakeri (NC_033388), Nilaparvata muiri (NC_024627), Nilaparvata lugens (NC_021748), Peregrinus maidis (NC_037182) and Saccharosydne procerus (NC_042179) and two out group species: Magicicada tredecula (NC_041656) and Philaenus spumarius (NC_005944). Numbers above branches indicate the posterior probability values of Bayesian inference. Family and infraorder names were displayed with gray bars in the right side of phylogenetic tree.

Figure 1. Bayesian inference (1,000,000 generations) phylogenetic tree of 16 Fulgoromorpha mitochondrial genomes: Ricania speculum (this study, MT834932), Ricania speculum (NC_031369), Ricania marginalis (NC_019597), Ricania sp. (MK251070), Geisha distinctissima (NC_012617), Aphaena amabilis (NC_045075), Lycorma delicatula (NC_012835), Sogatella furcifera (NC_021417), Sogatella vibix (NC_042180), Laodelphax striatellus (NC_013706), Changeondelphax velitchkovskyi (NC_037181), Nilaparvata bakeri (NC_033388), Nilaparvata muiri (NC_024627), Nilaparvata lugens (NC_021748), Peregrinus maidis (NC_037182) and Saccharosydne procerus (NC_042179) and two out group species: Magicicada tredecula (NC_041656) and Philaenus spumarius (NC_005944). Numbers above branches indicate the posterior probability values of Bayesian inference. Family and infraorder names were displayed with gray bars in the right side of phylogenetic tree.

Disclosure statement

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

Data availability statement

Mitochondrial genome sequence can be accessed via accession number MT834932 in GenBank of NCBI at https://www.ncbi.nlm.nih.gov. The associated BioProject, SRA, and Bio-Sample numbers are PRJNA665776, SAMN16268313, and SRR12717137, respectively.

Additional information

Funding

This research was supported by “Cooperative Research Program for Agriculture Science & Technology Development (Project no. PJ013389012020)” Rural Development Administration, Republic of Korea.

References

  • Bae Y, Park J, Lee W. 2020. The complete mitochondrial genome of Aphis gossypii Glover, 1877 (Hemiptera: Aphididae) isolated from Plantago asiatica in Korea. Mitochondrial DNA Part B. 5(3):2896–2898.
  • Bolger AM, Lohse M, Usadel B. 2014. Trimmomatic: a flexible trimmer for Illumina sequence data. Bioinformatics. 30(15):2114–2120.
  • Hongoh Y, Ishikawa H. 2000. Evolutionary studies on uricases of fungal endosymbionts of aphids and planthoppers. J Mol Evol. 51(3):265–277.
  • Huelsenbeck JP, Ronquist F. 2001. MRBAYES: Bayesian inference of phylogenetic trees. Bioinformatics. 17(8):754–755.
  • Katoh K, Standley DM. 2013. MAFFT multiple sequence alignment software version 7: improvements in performance and usability. Mol Biol Evol. 30(4):772–780.
  • Lee Y-D, Lee J, Kim D-S, Park J, Xi H, Roh J, Kim D-S, Nam SJ, Kim S-K, Song J-Y, et al. 2020. The complete mitochondrial genome of Hipparchia autonoe (Esper, 1783) (Lepidoptera: Nymphalidae): investigation of intraspecific variations on mitochondrial genome. Mitochondrial DNA Part B. 5(2):1542–1544.
  • Li H. 2013. Aligning sequence reads, clone sequences and assembly contigs with BWA-MEM. arXiv preprint arXiv:13033997.
  • Li H, Handsaker B, Wysoker A, Fennell T, Ruan J, Homer N, Marth G, Abecasis G, Durbin R. 2009. The sequence alignment/map format and SAMtools. Bioinformatics. 25(16):2078–2079.
  • Mazza G, Pennacchio F, Gargani E, Franceschini I, Roversi PF, Cianferoni F. 2014. First report of Ricania speculum (Walker, 1851) in Europe (Hemiptera: Fulgoromorpha: Ricaniidae). Zootaxa. 3861(3):297–300.
  • Park J, Xi H, Kim Y, Park J, Lee W. 2019. The complete mitochondrial genome of Aphis gossypii Glover, 1877 (Hemiptera: Aphididae) collected in Korean peninsula. Mitochondrial DNA Part B. 4(2):3007–3009.
  • Park J, Xi H, Park J, Lee W. 2020. The complete mitochondrial genome of fungal endosymbiont, Ophiocordycipitaceae sp., isolated from Ricania speculum (Hemiptera: Ricaniidae). Mitochondrial DNA Part B. 5(2):1888–1889.
  • Park J, Xi H, Park J, Nam SJ, Lee Y-D. 2020. Complete genome sequence of the Blochmannia endosymbiont of Camponotus nipponensis. Microbiol Resour Announc. 9(29):e00703–e00720.
  • Rossi E, Lucchi A. 2015. The Asian planthopper Ricania speculum (Walker) (Homoptera: Ricaniidae) on several crops in Italy: a potential threat to the EPPO region? EPPO Bull. 45(1):119–122.
  • Rossi E, StroiŃski A, Lucchi A. 2015. Egg morphology, laying behavior and record of the host plants of Ricania speculum (Walker, 1851), a new alien species for Europe (Hemiptera: Ricaniidae). Zootaxa. 4044(1):93–104.
  • Song N, Liang A-P. 2013. A preliminary molecular phylogeny of planthoppers (Hemiptera: Fulgoroidea) based on nuclear and mitochondrial DNA sequences. PLOS One. 8(3):e58400.
  • Wilson SW1, Rossi E, Lucchi A. 2016. Descriptions of the adult genitalia and immatures of the Asian planthopper Ricania speculum (Hemiptera: Fulgoroidea: Ricaniidae) recently introduced to Italy. Ann Entomol Soc Am. 109(6):899–905.
  • Zerbino DR, Birney E. 2008. Velvet: algorithms for de novo short read assembly using de Bruijn graphs. Genome Res. 18(5):821–829.
  • Zhang Q-X, Guan D-L, Niu Y, Sang L-Q, Zhang X-X, Xu S-Q. 2016. Characterization of the complete mitochondrial genome of the Asian planthopper Ricania speculum (Hemiptera: Fulgoroidea: Ricannidae). Conserv Genet Resour. 8(4):463–466.
  • Zhao QY, Wang Y, Kong YM, Luo D, Li X, Hao P. 2011. Optimizing de novo transcriptome assembly from short-read RNA-Seq data: a comparative study. BMC Bioinformatics; 2011: BioMed Central. 12(Suppl 14):S2.