https://orcid.org/0000-0001-8091-0407
Abstract
The mitochondrial genome has been widely used in the study of phylogeny and species-level evolution. Here, we sequenced and analyzed the full mitogenome of Tetrastichus howardi, an important natural enemy of many lepidopteran pests. The complete mitochondrial genome has 14,791 nucleotides, containing 13 protein-coding genes (PCGs), 2 ribosomal RNA genes (rRNAs), 22 transfer RNA genes (tRNAs), and a partial control region. All the 13 PCGs started with typical ATN (ATA, ATG, and ATT) codon. Among 13 PCGs, nine genes terminated with the stop codon TAA and four genes terminated with T. Our study provides information on comparative mitogenomics of Eulophidae.
Tetrastichus howardi (Hymenoptera: Eulophidae), an exotic pupae parasitoid of a great number of Lepidoptera pests, has a great potential for biological control (Favoreto et al. Citation2021; Rodrigues et al. Citation2021). The parasitoid was introduced into South Africa from the Philippines as a potential biocontrol agent for two important Lepidoptera pests, the Chilo partellus Swinhoe and Busseola fusca Fuller, which were serious pests of corn and sorghum in the country (Skoroszewski and Hamburg Citation1987; Kfir Citation2001). The phylogenetic relationship of the Eulophidae based on molecular and morphology analyses has been conducted, but the position of the genera in the family is not clear (Tian et al. Citation2021; Tang et al. Citation2021). In this paper, we analyzed the mitochondrial genome of the T. howardi for the first time, which provided valuable information for further study on molecular systematics and phylogenetic relationships within the Eulophidae.
In this study, the specimen of T. howardi was collected from the Environment and Plant Protection Institute, Chinese Academy of Tropical Agriculture Sciences, Hainan, China (19°59′21′′N, 110°20′9′′E) and deposited in the Post-Entry Quarantine Station for Tropical Plant, Haikou Customs District P.R. China (http://www.rdzw.net.cn/, [email protected]) under the voucher number IN07040201-0001-00010. Total genomic DNA of T. howardi was extracted from a single sample by CTAB method (Reineke et al. Citation1998). The mitogenome sequence was generated using Illumina HiSeq X TEN Sequencing System with 150 bp paired-end reads. A total of 1.17 Gb clean data was obtained and assembled by the MITObim software (Hahn et al. Citation2013) with the complete mt genome (MN123622) of Tamarixia radiata as a reference. The annotations were mainly compared with the existing mitochondrial genomes of related species, and the annotation results were confirmed by MITOS webserver (Bernt et al. Citation2013). The sequence was submitted to GenBank under the accession number MZ334468.
The nucleotide composition of the mitogenome sequence is 43.8% A, 41.7% T, 8.0% C, 6.5% G, with a high AT bias of 85.5%. All the genes were distributed on two coding chains, of which 26 genes were encoded on the majority strand (J-chains) and 11 genes were encoded on minority chains (N-chains). All 13 protein-coding genes get off by the conventional ATN as start codon, including six ATGs (ATP6, COX1, COX3, CYTB, NAD4, and NAD6), three ATTs (COX2, NAD3, and NAD4L), and four ATAs (ATP8, NAD1, NAD2, and NAD5). Nine (ATP8, COX1, CYTB, NAD2, NAD3, NAD4, NAD4L, NAD5, and NAD6) of 13 protein-coding genes use TAA for the stop codon, while ATP6, COX2, COX3, and NAD1 genes use the incomplete stop codon, T. The length of tRNA genes ranges from 60 to 71 bp and all of them can be folded into the typical cloverleaf secondary structure, except for trnS1, which lack dihydrouridine (DHU) arm. In previous studies, there is a lack of DHU arm in the trnS1 of many insect mitochondrial genomes (Xiong et al. Citation2019; Tang et al. Citation2021). The lrRNA is 1320 bp in length with an A + T content of 88.1%, and the srRNA is 748 bp long with an A + T content of 81.9%.
Based on the concatenated amino acid sequences of 13 PCGs, the phylogenetic relationship of T. howardi with 22 other parasitoid and an outgroup species (Tetranychus urticae) was constructed by the neighbor-joining method using MEGA 7.0 software (Sudhir et al. Citation2016). The result showed that T. howardi was closely related to T. radiata (), which agrees with the conventional classification.
Figure 1. Phylogenetic tree showing the relationship between Tetrastichus howardi and 22 other parasitoids based on neighbor-joining method. GenBank accession numbers of each species were listed in the tree. The parasitoid determined in this study has been marked red. Bootstrap values are shown above the nodes.
Disclosure statement
No potential conflict of interest was reported by the author(s).
Data availability statement
The genome sequence data that support the findings of this study are openly available in GenBank of NCBI at (https://www.ncbi.nlm.nih.gov/) under the accession no. MZ334468. The associated BioProject, SRA, and Bio-Sample numbers are PRJNA737393, SRR14812766, and SAMN19693530, respectively.
Additional information
Funding
References
- Bernt M, Donath A, Jühling F, Externbrink F, Florentz C, Fritzsch G, Pütz J, Middendorf M, Stadler PF. 2013. MITOS: improved de novo metazoan mitochondrial genome annotation. Mol Phylogenet Evol. 69(2):313–319.
- Favoreto AL, Pavani RF, Ribeiro MF, Zanuncio AJV, Soares MA, Zanuncio JC, Wilcken CF, Favoreto AL, Pavani RF, Ribeiro MF. 2021. Tetrastichus howardi (Hymenoptera: Eulophidae): first report of parasitism in Oxydia vesulia (Lepidoptera: Geometridae). Braz J Biol. 81(2):406–410.
- Hahn C, Bachmann L, Chevreux B. 2013. Reconstructing mitochondrial genomes directly from genomic next-generation sequencing reads—a baiting and iterative mapping approach. Nucleic Acids Res. 41(13):e129.
- Kfir R. 2001. Prospects for biological control of Chilo partellus in grain crops in South Africa. Int J Trop Insect Sci. 21(4):275–280.
- Reineke A, Karlovsky P, Zebitz C. 1998. Preparation and purification of DNA from insects for AFLP analysis. Insect Mol Biol. 7(1):95–99.
- Rodrigues A, Pereira FF, Barbosa PRR, Silva-Torres CSA, Torres JB. 2021. Parasitism behavior of Tetrastichus howardi (Hymenoptera: Eulophidae) on larvae and pupae of sugarcane borers. J Insect Behav. 34:71–81.
- Skoroszewski RW, Hamburg HV. 1987. The release of Apanteles flavipes (Cameron) (Hymenoptera: Braconidae) against stalk-borers of maize and grain sorghum in South Africa. J Entomol Soc S Afr. 50:249–255.
- Sudhir K, Glen S, Koichiro T. 2016. MEGA 7: Molecular Evolutionary Genetics Analysis Version 7.0 for bigger datasets. Mol Biol Evol. 33(7):1870.
- Tang X, Lyu BQ, Lu H, Tang JH, Meng R, Cai B. 2021. The mitochondrial genome of a parasitic wasp, Chouioia cunea Yang (Hymenoptera: Chalcidoidea: Eulophidae) and phylogenetic analysis. Mitochondrial DNA B Resour. 6(3):872–874.
- Tian X-C, Xian X-Q, Zhang G-F, Castañé C, Romeis J, Wan F-H, Zhang Y-B. 2021. Complete mitochondrial genome of a predominant parasitoid, Necremnus tutae (Hymenoptera: Eulophidae) of the South American tomato leafminer Tuta absoluta (Lepidoptera: Gelechiidae). Mitochondrial DNA B Resour. 6(2):562–563.
- Xiong M, Zhou QS, Zhang YZ. 2019. The complete mitochondrial genome of Encyrtus infelix (Hymenoptera: Encyrtidae). Mitochondrial DNA Part B-Resour. 4(1):114–115.