Next Article in Journal
Genetic and Morphological Variation in Hypodontia of Maxillary Lateral Incisors
Next Article in Special Issue
Genetic Diversity in the mtDNA of Physarum polycephalum
Previous Article in Journal
Transcriptome Profiling of Developing Testes and First Wave of Spermatogenesis in the Rat
Previous Article in Special Issue
Complete Mitogenome of the Triplophysa bombifrons: Comparative Analysis and Phylogenetic Relationships among the Members of Triplophysa
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Genes
Volume 14
Issue 1
10.3390/genes14010230
genes-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Communication

The Mitochondrial Genomes of Two Parasitoid Wasps Protapanteles immunis and Parapanteles hyposidrae (Hymenoptera: Braconidae) with Phylogenetic Implications and Novel Gene Rearrangements

by
Dandan Xiao
1,2,3,4,
Ziqi Wang
3,4,
Jiachen Zhu
2,3,4,
Xiaogui Zhou
5,
Pu Tang
1,2,3,4,* and
Xuexin Chen
1,2,3,4
1
Hainan Institute, Zhejiang University, Sanya 572025, China
2
Guangdong Laboratory for Lingnan Modern Agriculture, Guangzhou 510642, China
3
State Key Laboratory of Rice Biology, Ministry of Agriculture Key Lab of Molecular Biology of Crop Pathogens and Insects, and Zhejiang Provincial Key Laboratory of Biology of Crop Pathogens and Insects, Zhejiang University, Hangzhou 310058, China
4
Institute of Insect Sciences, College of Agriculture and Biotechnology, Zhejiang University, Hangzhou 310058, China
5
Ministry of Agriculture Key Laboratory of Tea Quality and Safety Control, Tea Research Institute of Chinese Academy of Agricultural Sciences, Hangzhou 310008, China
*
Author to whom correspondence should be addressed.
Genes 2023, 14(1), 230; https://doi.org/10.3390/genes14010230
Submission received: 19 December 2022 / Revised: 8 January 2023 / Accepted: 10 January 2023 / Published: 16 January 2023
(This article belongs to the Special Issue Advanced Research on Mitochondrial Genome)
Browse Figures
Review Reports Versions Notes

Abstract

:
Parapanteles hypsidrae (Wilkinson, 1928) and Protapanteles immunis (Haliday, 1834) are the most important parasitic wasps of Ectropis grisescens Warren and Ectropis obliqua (Prout). We sequenced and annotated the mitochondrial genomes of Pa. hyposidrae and Pr. immunis, which are 17,063 bp and 16,397 bp in length, respectively, and possess 37 mitochondrial genes. We discovered two novel types of gene rearrangement, the local inversion of nad4L in Pa. hyposidrae and the remote inversion of the block cox3-nad3-nad5-nad4 in Pr. immunis, within the mitogenomes of Braconidae. The phylogenetic analysis supported the subfamily Microgastrinae is a monophyletic group, but the tribes Apantelini and Cotesiini within this subfamily are paraphyletic groups.

1. Introduction

Microgastrinae is one of the most numerous subfamilies of the Braconidae with around 3000 known species worldwide. It is present in all of the major terrestrial ecosystems. The majority of the Microgastrinae are lepidopteran koinobiont endoparasitoids. More than 100 species of microgastrines have been utilized in biological control projects, making them one of the most valuable groups in agricultural and forestry pest biological control [1]. Unfortunately, phylogenetic studies on this subfamily are limited and there is still no robust phylogeny for it. The boundaries of several genera are now unclear and occasionally contradictory, and, in addition, future research is likely to modify the current cognitive way of many groups [2,3,4].
Pa. hypsidrae and Pr. immunis (Hymenoptera: Braconidae: Microgastrinae) are important parasitic wasps of E. grisescens and E. obliqua, which are two of the most destructive chewing pests in China’s tea plantations. Pa. hyposidrae has a black head and mesosoma, a dark brown metasoma with white areas, and a smooth and largely yellowish-green cocoon whereas Pr. immunis has a black body with yellow hind femora, and a white cocoon covered with fluffy cotton-like filaments [5]. When the conditions are suitable, the highest natural parasitism rate of tea geometrid larvae by these two parasitoids can reach more than 90%, playing a vital role in the population control of Ectropis spp. [5]. Under higher temperatures, Pr. immunis outperformed Pa. hyposidrae. It was found that parasitism rates decreased with parasitoid density at different temperatures, resulting in less efficient searching [6].
The mitochondrial genome plays a crucial role in phylogenetic construction, and it is inherited from the maternal lineage and cannot be combined with other mitochondrial lineages of insects [7]. The typical insect mitochondrial genome is a circular molecule with a size of 14–19 kb that contains 37 genes, 13 of which are protein-coding genes (PCGs), two of which are ribosomal RNA (rRNA) genes, and 22 of which are transfer RNA (tRNA) genes [8]. The gene arrangement is relatively conservative in insects [9]. However, rearrangement events involving multiple combinations of protein coding genes, tRNA genes and rRNA genes have been observed in Anoplura, Thysanoptera, Corrodentia, Hymenoptera, and other insects [10,11,12,13]. The evolution of gene rearrangements in insect mitogenomes is a hot topic [14]. The rate of gene rearrangement in hymenopteran mitogenomes is extremely high [15,16,17]. Gene rearrangements are typically restricted to specific lineages, which can aid in phylogenetic reconstruction, such as the Braconidae subfamily level [15]. However, there has been very little progress in sequencing the mitochondrial genomes of Microgastrinae. Currently, only sixteen verified mitochondrial genomes of the subfamily Microgastrinae have been reported in GenBank (https://www.ncbi.nlm.nih.gov/; accessed on 12 November 2022), and only two species have all 37 genes identified.
Two mitogenomes of Microgastrinae were newly sequenced using next generation sequencing in this study. We obtained the mitochondrial genomes of Pa. hyposidrae and Pr. immunis, which provided a detailed description of their genomic characteristics as well as more accurate and extensive information for further studying the gene arrangement and the evolutionary history of Microgastrinae.

2. Materials and Methods

2.1. Sample Identification and DNA Extraction

Before DNA extraction, the obtained specimen was preserved in 100% ethanol at −80 °C. The specimens of two parasitoids were obtained through rearing in the laboratory from the Tea Research Institute, Chinese Academy of Agricultural Sciences. Following the manufacturer’s instructions, genomic DNA was extracted using FastPure Cell/Tissue DNA Isolation Mini Kit (Vazyme Biotech Co., Ltd., Nanjing, China). The voucher specimen was preserved in the Parasitic Hymenoptera Collection (Institute of Insect Sciences, Zhejiang University).

2.2. Next-Generation Sequencing and Assembly

The library was created using the VAHTS™ Universal DNA Library Prep Kit for Illumina® v9.1, and sequenced on an Illumina HiSeq platform with 150 bp pared-end read length by Novogene. FastQC was utilized to examine the raw data, and Trimmomatic was utilized with default settings to trim adaptors and indexes [18,19]. The FastqExtract script was used to filter out the target mitochondrial reads by running BLASTn (E value: 1 × 10−5) against a reference data set of Braconidae mitochondrial genomes [20]. The mitochondrial reads were assembled by SPAdes v3.0 [21] and IDBA v1.1.3 with default values, respectively [22]. Following that, GENEIOUS Prime v2020.0.5 (Biomatters Ltd., San Diego, CA, USA) was used to integrate two assemblies.

2.3. Mitochondrial Genome Annotation and Analysis

The MITOS Web Serve was used to annotate the assembled genomes [23]. Against the reported mitogenomes of Braconidae, the start and stop positions of protein-coding genes were manually modified in Geneious Prime v2020.0.5. The predicted tRNA genes were verified by the tRNAscan-SE search site with their homologs from related species [24]. Two maps of mitochondrial genomes were created by CGView server online V 1.0 (http://cgview.ca/, accessed on 25 February 2022) [25]. The obtained mitogenomes were uploaded to GenBank (OP741148 and OP741149).
Gene rearrangements were investigated by contrasting them with the putative ancestral mitogenome of Drosophila melanogaster. The base composition of all components was estimated using MEGA 11.0 [26]. The following formulas were used to calculate AT-skew and GC-skew: AT-skew = (A% − T%)/(A% + T%) and GC-skew = (G% − C%)/(G% + C%) [27]. Geneious Prime v2020.0.5 was used to calculate the relative synonymous codon usage (RSCU) of all PCGs. DnaSP v6.12.03 was used to compute the Synonymous (Ks) and non-synonymous (Ka) substitution rates of PCGs [28].

2.4. Phylogenetic Analysis

G-INS-i methods implemented in MAFFT v7.464 were used to align the PCGs [29]. PartitionFinder v2 was used to find the optimal partition schemes of substitution models for the matrix, with model selection = BIC and Branch lengths = unlinked across different subsets [30]. Based on nucleotide sequences of all PCGs, Mrbayes v3.2.7a [31] and RAxML-HPC2 v8.2.12 [32] were used to construct the phylogenetic tree, respectively. Four Markov chains were run simultaneously for 10 million generations for Bayesian inference analysis (BI), with tree sampling occurring every 1000 generations and a burn-in of 25% of the trees for phylogenetic analysis. In maximum likelihood (ML) analysis, ML trees with 1000 bootstrap replications are built using the GTRGAMMA model and 200 runs for various individual partitions.

3. Results

3.1. General Features of Mitochondrial Genomes

The mitochondrial genomes of Pa. hyposidrae and Pr. immunis are 17,063 bp and 16,397 bp in length, respectively (Table 1, Figure 1). 37 mitochondrial genes were identified in Pa. hyposidrae and Pr. immunis, containing 13 PCGs, 22 tRNA genes, and 2 ribosomal RNA (rRNA) genes. The complete control region of the two species was unable to be assembled, probably owing to poor similarity among reference sequences and high duplications of A and T, which are typical in hymenopteran mitogenomes [33].
There are six overlapping regions in the mitochondrial genome of Pa. hyposidrae, with a total of 23 bp, as the largest overlap was 7 bp, located at two junctions (trnW-trnC, trnS2-nad1) and the smallest overlap was 1 bp, located at one junction (nad2-trnW). There are 21 intergenic regions with a total length of 2271 bp, as the largest length of the gene spacer was 1129 bp, located at one junction (trnP-trnT) and the smallest length of the gene spacer was 1 bp, located at two junctions (trnQ-nad2, cox2-trnH). The length of rrnL and rrnS is 1127 bp and 751 bp, respectively, and the length of tRNA is 54–74 bp.
There are three overlapping regions in the mitochondrial genome of Pr. immunis with a total of 18 bp, as the largest overlap was 10 bp, located at one junction (atp8-atp6) and the smallest overlap was 1 bp in length, located at one junction (trnS1-trnA). There are 29 intergenic regions with a total length of 1849 bp, as the largest length of the gene spacer was 579 bp, located at one junction (nad5-trnS1) and the smallest length of the gene spacer was 1 bp, located at three junctions (trnQ-nad2, nad4L-trnT, trnP-ad6). The length of rrnL and rrnS is 1137 bp and 661 bp, respectively, and the length of tRNA is 60–73 bp.
All PCGs in Pa. hyposidrae and Pr. immunis mitochondrial genomes started with a typical ATN codon, as in other Hymenoptera insects [34,35,36]. Among the starting codons used by the mitochondrial genome of Pa. hyposidrae, ATT was used the most, seven times, followed by ATG four times, and ATA only twice. All protein coding genes except cox3 use T as a termination codon, and the other twelve protein coding genes use TAA. Incomplete termination codons often exist in the protein coding genes of the mitochondrial genome of arthropods. After being transcribed into mRNA, they are supplemented by 3 ‘end polyadenylation [37,38]. Among the start codons used by the mitochondrial genome of Pr. immunis, ATG was used the most, six times, followed by ATT five times, and ATA only twice. All protein coding genes use TAA as the termination codon.

3.2. Base Composition, and Codon Usage

The sequenced region of the mitogenomes of Pa. hyposidrae and Pr. immunis have 86.15% and 86.56% A + T content, respectively (Table 1). Compared with other orders, the relatively high A + T content in the hymenopteran mitogenomes is not exceptional [39]. The total A + T content of all PCGs of Pa. hyposidrae and Pr. immunis is 85.12% and 85.52%, respectively. Typically, nad6 or atp8 has the largest proportion of A + T content in the hymenopteran mitochondrial genome [40]. In the mitochondrial genomes of Pa. hyposidrae and Pr. immunis, the highest A + T content is atp8 (Pa. hyposidrae: 94.55%, Pr. immunis: 93.10%), and the lowest A + T content is cox1 (Pa. hyposidrae: 77.30%, Pr. immunis: 77.71%) (Figure 2). In Pa. hyposidrae, the A + T content of nad6 is 92.99% while the A + T content of nad6 is 91.67% in Pr. immunis. Most of the AT skew is negative in both Pa. hyposidrae and Pr. immunis. In Pa. hyposidrae, the largest AT skew is nad4 (0.1202), and the lowest is atp6 (−0.2509). In Pr. immunis, the largest AT skew is nad3 (0.1485), and the lowest is atp6 (−0.2076). Similarly, the GC skew also has positive and negative values. Most of the GC skew is positive in both Pa. hyposidrae and Pr. immunis. In Pa. hyposidrae, the largest GC skew is cox2 (0.2137), and the smallest is atp8 (−0.3333). In Pr. immunis, the largest GC skew is nad2 (0.1959), and the smallest is nad3 (−0.2223).
The relative synonymous codon usage (RSCU) in the mitochondrial genomes of Pa. hyposidrae and Pr. immunis displayed a strong preference toward A and T, particularly at the third position. The total number of codons in the Pa. hyposidrae and Pr. immunis mitochondrial genomes is 3751 and 3678, respectively. The four most frequent used amino acids in both Pa. hyposidrae and Pr. immunis, correspond to the following codons: UUA (Leu), UUU (Phe), AUU (Ile), and AUA (Met) (Tables S1 and S2). These findings are similar with reported mitogenomes of other hymenopterans [41,42].

3.3. Gene Rearrangements

Gene rearrangement events are typically divided into four types: translocation, gene shuffling (local translocation), local inversion (inverted in the local position), and remote inversion (translocated and inverted) [11]. Compared with the putative ancestral mitogenome of D. melanogaster, the PCGs and tRNAs have various degrees of rearrangement in the mitogenomes of the two parasitoid wasps (Figure 3).
Frequent gene rearrangements have been reported in Braconidae in previous research. The majority of genes rearranged were tRNA genes, with three hotspots for tRNA gene rearrangement proposed: trnA-trnR-trnN-trnS1-trnE-trnF, trnK-trnD and trnT-trnP [15]. In the two species, Pa. hyposidrae and Pr. Immunis, the three hot spots for tRNA gene rearrangement are also discovered. trnH is remotely inverted to the junction of cox2 and atp8, forming the organization pattern trnH-trnD-trnK in Pa. hyposidrae, and trnH-trnK-trnD in Pr. Immunis. trnA-trnR-trnN-trnS1-trnE-trnF is rearranged into trnN-trnE-trnA-trnS1 and trnS1-trnA-trnE-trnN in Pa. hyposidrae and Pr. Immunis, respectively. trnT-trnP is remotely inverted from downstream to upstream of nad4L in Pa. hyposidrae, forming the arrangement pattern trnP-trnT. Additionally, the clusters trnW-trnC-trnY and trnI-trnQ-trnM are also rearranged. Both species have locally inverted trnY. trnI and trnM were remotely inverted, resulting in the arrangement pattern trnM-trnI-trnQ in Pa. hyposidrae, and trnI-trnM-trnQ in Pr. Immunis.
With the exception of Chelonus sp., Cotesia vestalis, and Stenocorse bruchivorus, all protein coding genes in all Braconidae species were conserved as the putative ancestral arrangement of insects [16]. The rearrangement of protein coding genes was identified in Pa. hyposidrae and Pr. Immunis. Nad4L is locally inverted in Pa. hyposidrae and a large block (cox3-nad3-nad5-nad4) is remotely inverted to nad4-nad5-nad3-cox3 in Pr. Immunis, both of which are unique gene rearrangement patterns in the Braconidae.

3.4. Phylogenetic Analyses

To validate the phylogenetic position of Pa. hyposidrae and Pr. Immunis within Microgastrinae, a phylogenetic analysis based on all PCGs was constructed in this study using Bayesian and maximum likelihood inferences with other 14 organisms, including 12 species from the subfamily Microgastrinae (the tribe concepts follows a classification based on Manson [43], which is widely used) and two from the helconoid subfamily complexes (Macrocentrinae and Agathidinae) as outgroup (Figure 4, Table S3). The topology of BI tree is the same as that of ML tree, though some clades with low support values, which supports that the subfamily Microgastrinae is a monophyletic group, and is a sister group to the helconoid subfamily complexes [15,44]. Apantelini and Cotesiini split into two and three clusters, separately, indicating they are paraphyletic, which is supported by some previous studies [2,4]. The results of this study show that Pa. hyposidrae and Pr. Immunis belong to Apantelini and Cotesiini, respectively. Several genera with unclear boundaries, as Fernandez-Triana et al. pointed out, including Protapanteles and Diolcogaster, are most likely polyphyletic and will need to be split into several genera [3]. Pr. Immunis is sister to C. vestalis and Cotesia flaps rather than Protapanteles sp., indicating that the genera Protapanteles is paraphyletic.

4. Conclusions

In this study, we successfully obtained the two mitogenomes of Pa. hypsidrae and Pr. Immunis, which are the most important parasitic wasps of E. grisescens and E. obliqua. Combined with reported Braconidae mitogenomes, we conducted analyses of base composition, codon usage, gene rearrangement, and phylogeny within Braconidae. Two novel gene rearrangement types were discovered in the two newly acquired mitogenomes, the local inversion of nad4L in Pa. hyposidrae and the remote inversion of the block nad4-nad5-nad3-cox3 in Pr. Immunis. The phylogenetic analysis supported the subfamily Microgastrinae is a monophyletic group, but the tribes Apantelini and Cotesiini within this subfamily are paraphyletic groups. Nevertheless, it is crucial to grasp these results in order to comprehend the evolution of Microgastrinae. For further research into gene organization and the evolutionary history of Microgastrinae, denser taxon sampling yields more accurate and thorough data.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/genes14010230/s1, Table S1: Codon usage in the mitochondrial genome of Pa. hypsidrae; Table S2: Codon usage in the mitochondrial genome of Pr. immunis; Table S3: Information of mitochondrial genomes used in phylogenetic analysis.

Author Contributions

Conceptualization, D.X., X.Z., P.T. and X.C.; methodology, D.X., Z.W. and J.Z.; formal analysis, D.X., Z.W. and J.Z.; data curation, D.X., Z.W. and J.Z.; writing—original draft preparation, D.X.; writing—review and editing, D.X., P.T. and X.C.; project administration, P.T and X.C. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the Key Project of Laboratory of Lingnan Modern Agriculture (NT2021003), the Key International Joint Research Program of National Natural Science Foundation of China (31920103005), the General Program of National Natural Science Foundation of China (32070467), the National Key Research and Development Plan (2019YFD1002100), the Provincial Key Research and Development Plan of Zhejiang (2021C02045), and the Fundamental Research Funds for the Central Universities (2021FZZX001-31).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The data supporting the findings of this study are openly available in the National Center for Biotechnology Information (https://www.ncbi.nlm.nih.gov), accession numbers: OP741148 and OP741149. Raw sequence reads for each specimen-specific library were deposited in the BioProject PRJNA896708. The SRR codes for the two samples are SRR22142645 and SRR22142646.

Acknowledgments

We thank Xiqian Ye for providing access to the cluster computer system used for the phylogenetic analyses.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Whitfield, J.B. Annotated checklist of the Microgastrinae of North America north of Mexico (Hymenoptera: Braconidae). J. Kans. Entomol. Soc. 1995, 68, 245–262. [Google Scholar]
  2. Banks, J.C.; Whitfield, J.B. Dissecting the ancient rapid radiation of microgastrine wasp genera using additional nuclear genes. Mol. Phylogenet. Evol. 2006, 41, 690–703. [Google Scholar] [CrossRef]
  3. Fernandez-Triana, J.; Shaw, M.R.; Boudreault, C.; Beaudin, M.; Broad, G.R. Annotated and illustrated world checklist of Microgastrinae parasitoid wasps (Hymenoptera, Braconidae). ZooKeys 2020, 920, 1. [Google Scholar] [CrossRef] [Green Version]
  4. Whitfield, J.B.; Mardulyn, P.; Austin, A.D.; Dowton, M. Phylogenetic relationships among microgasrinae braconid wasp genera based on data from the 16S, COI and 28S genes and morphology. Syst. Entomol. 2002, 27, 337–359. [Google Scholar] [CrossRef]
  5. Zhou, X.; Tang, P.; Wu, Q.; Guo, H.; Xiao, Q.; Chen, X.X. Identification of Two Common Larval Parasitic Wasps of Ectropis obliqua and Ectropis grisescens (Lepidoptera: Geometridae). Chin. J. Biol. Control 2022, 1–14. [Google Scholar] [CrossRef]
  6. Wang, Z.Q.; Zhou, X.G.; Xiao, Q.; Tang, P.; Chen, X.X. The Potential of Parapanteles hyposidrae and Protapanteles immunis (Hymenoptera: Braconidae) as Biocontrol Agents for the Tea Grey Geometrid Ectropis grisescens (Lepidoptera). Insects 2022, 13, 937. [Google Scholar] [CrossRef]
  7. Cameron, S.L.; Johnson, K.P.; Whiting, M.F. The mitochondrial genome of the screamer louse Bothriometopus (Phthiraptera: Ischnocera): Effects of extensive gene rearrangements on the evolution of the genome. J. Mol. Evol. 2007, 65, 589–604. [Google Scholar] [CrossRef]
  8. Cameron, S.L. Insect Mitochondrial Genomics: Implications for Evolution and Phylogeny. Annu. Rev. Entomol. 2014, 59, 95–117. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  9. Ballard, J.W.O. Comparative genomics of mitochondrial DNA in members of the Drosophila melanogaster subgroup. J. Mol. Evol. 2000, 51, 48–63. [Google Scholar] [CrossRef]
  10. Dowton, M.; Austin, A.D. Evolutionary dynamics of a mitochondrial rearrangement “hot spot” in the hymenoptera. Mol. Biol. Evol. 1999, 16, 298–309. [Google Scholar] [CrossRef] [Green Version]
  11. Dowton, M.; Castro, L.R.; Campbell, S.L.; Bargon, S.D.; Austin, A.D. Frequent mitochondrial gene rearrangements at the hymenopteran nad3-nad5 junction. J. Mol. Evol. 2003, 56, 517–526. [Google Scholar] [CrossRef]
  12. Shao, R.; Campbell, N.J.H.; Barker, S.C. Numerous gene rearrangements in the mitochondrial genome of the wallaby louse, Heterodoxus macropus (Phthiraptera). Mol. Biol. Evol. 2001, 18, 858–865. [Google Scholar] [CrossRef] [Green Version]
  13. Shao, R.F.; Campbell, N.J.H.; Schmidt, E.R.; Barker, S.C. Increased rate of gene rearrangement in the mitochondrial genomes of three orders of hemipteroid insects. Mol. Biol. Evol. 2001, 18, 1828–1832. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  14. Du, C.; Zhang, L.F.; Lu, T.; Ma, J.N.; Zeng, C.J.; Yue, B.S.; Zhang, X.Y. Mitochondrial genomes of blister beetles (Coleoptera, Meloidae) and two large intergenic spacers in Hycleus genera. BMC Genom. 2017, 18, 15. [Google Scholar] [CrossRef] [Green Version]
  15. Li, Q.; Wei, S.J.; Tang, P.; Wu, Q.; Shi, M.; Sharkey, M.J.; Chen, X.X. Multiple Lines of Evidence from Mitochondrial Genomes Resolve Phylogenetic Relationships of Parasitic Wasps in Braconidae. Genome Biol. Evol. 2016, 8, 2651–2662. [Google Scholar] [CrossRef] [Green Version]
  16. Wei, S.J.; Shi, M.; Sharkey, M.J.; van Achterberg, C.; Chen, X.X. Comparative mitogenomics of Braconidae (Insecta: Hymenoptera) and the phylogenetic utility of mitochondrial genomes with special reference to Holometabolous insects. BMC Genom. 2010, 11, 16. [Google Scholar] [CrossRef] [Green Version]
  17. Zheng, X.Y.; Cao, L.J.; Chen, P.Y.; Chen, X.X.; van Achterberg, K.; Hoffmann, A.A.; Liu, J.X.; Wei, S.J. Comparative mitogenomics and phylogenetics of the stinging wasps (Hymenoptera: Aculeata). Mol. Phylogenet. Evol. 2021, 159, 9. [Google Scholar] [CrossRef] [PubMed]
  18. Andrews, S.; Fast, Q.C. Available online: https://qubeshub.org/resources/fastqc (accessed on 30 December 2021).
  19. Bolger, A.M.; Lohse, M.; Usadel, B. Trimmomatic: A flexible trimmer for Illumina sequence data. Bioinformatics 2014, 30, 2114–2120. [Google Scholar] [CrossRef] [Green Version]
  20. Crampton-Platt, A.; Timmermans, M.; Gimmel, M.L.; Kutty, S.N.; Cockerill, T.D.; Khen, C.V.; Vogler, A.P. Soup to Tree: The Phylogeny of Beetles Inferred by Mitochondrial Metagenomics of a Bornean Rainforest Sample. Mol. Biol. Evol. 2015, 32, 2302–2316. [Google Scholar] [CrossRef] [PubMed]
  21. Bankevich, A.; Nurk, S.; Antipov, D.; Gurevich, A.A.; Dvorkin, M.; Kulikov, A.S.; Lesin, V.M.; Nikolenko, S.I.; Pham, S.; Prjibelski, A.D.; et al. SPAdes: A New Genome Assembly Algorithm and Its Applications to Single-Cell Sequencing. J. Comput. Biol. 2012, 19, 455–477. [Google Scholar] [CrossRef] [Green Version]
  22. Peng, Y.; Leung, H.C.M.; Yiu, S.M.; Chin, F.Y.L. IDBA-UD: A de novo assembler for single-cell and metagenomic sequencing data with highly uneven depth. Bioinformatics 2012, 28, 1420–1428. [Google Scholar] [CrossRef]
  23. Bernt, M.; Donath, A.; Juhling, F.; Externbrink, F.; Florentz, C.; Fritzsch, G.; Putz, J.; Middendorf, M.; Stadler, P.F. MITOS: Improved de novo metazoan mitochondrial genome annotation. Mol. Phylogenet. Evol. 2013, 69, 313–319. [Google Scholar] [CrossRef] [PubMed]
  24. Chan, P.P.; Lin, B.Y.; Mak, A.J.; Lowe, T.M. tRNAscan-SE 2.0: Improved detection and functional classification of transfer RNA genes. Nucleic Acids Res. 2021, 49, 9077–9096. [Google Scholar] [CrossRef] [PubMed]
  25. Grant, J.R.; Stothard, P. The CGView Server: A comparative genomics tool for circular genomes. Nucleic Acids Res. 2008, 36, W181–W184. [Google Scholar] [CrossRef]
  26. Tamura, K.; Stecher, G.; Kumar, S. MEGA11 Molecular Evolutionary Genetics Analysis Version 11. Mol. Biol. Evol. 2021, 38, 3022–3027. [Google Scholar] [CrossRef]
  27. Perna, N.T.; Kocher, T.D. Patterns of nucleotide composition at fourfold degenerate sites of animal mitochondrial genomes. J. Mol. Evol. 1995, 41, 353–358. [Google Scholar] [CrossRef]
  28. Rozas, J.; Ferrer-Mata, A.; Sanchez-DelBarrio, J.C.; Guirao-Rico, S.; Librado, P.; Ramos-Onsins, S.E.; Sanchez-Gracia, A. DnaSP 6: DNA Sequence Polymorphism Analysis of Large Data Sets. Mol. Biol. Evol. 2017, 34, 3299–3302. [Google Scholar] [CrossRef]
  29. Katoh, K.; Standley, D.M. MAFFT Multiple Sequence Alignment Software Version 7: Improvements in Performance and Usability. Mol. Biol. Evol. 2013, 30, 772–780. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  30. Lanfear, R.; Calcott, B.; Ho, S.Y.W.; Guindon, S. PartitionFinder: Combined Selection of Partitioning Schemes and Substitution Models for Phylogenetic Analyses. Mol. Biol. Evol. 2012, 29, 1695–1701. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  31. Ronquist, F.; Huelsenbeck, J.P. MrBayes 3: Bayesian phylogenetic inference under mixed models. Bioinformatics 2003, 19, 1572–1574. [Google Scholar] [CrossRef] [Green Version]
  32. Stamatakis, A. RAxML version 8: A tool for phylogenetic analysis and post-analysis of large phylogenies. Bioinformatics 2014, 30, 1312–1313. [Google Scholar] [CrossRef] [Green Version]
  33. Tang, P.; Zhu, J.C.; Zheng, B.Y.; Wei, S.J.; Sharkey, M.; Chen, X.X.; Vogler, A.P. Mitochondrial phylogenomics of the Hymenoptera. Mol. Phylogenet. Evol. 2019, 131, 8–18. [Google Scholar] [CrossRef]
  34. Cameron, S.L.; Dowton, M.; Castro, L.R.; Ruberu, K.; Whiting, M.F.; Austin, A.D.; Diement, K.; Stevens, J. Mitochondrial genome organization and phylogeny of two vespid wasps. Genome 2008, 51, 800–808. [Google Scholar] [CrossRef] [PubMed]
  35. Cha, S.Y.; Yoon, H.J.; Lee, E.M.; Yoon, M.H.; Hwang, J.S.; Jin, B.R.; Han, Y.S.; Kim, I. The complete nucleotide sequence and gene organization of the mitochondrial genome of the bumblebee, Bombus ignitus (Hymenoptera: Apidae). Gene 2007, 392, 206–220. [Google Scholar] [CrossRef] [PubMed]
  36. Wei, S.J.; Shi, M.; He, J.H.; Sharkey, M.; Chen, X.X. The complete mitochondrial genome of Diadegma semiclausum (Hymenoptera: Ichneumonidae) indicates extensive independent evolutionary events. Genome 2009, 52, 308–319. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  37. Ojala, D.; Montoya, J.; Attardi, G. tRNA punctuation model of RNA processing in human mitochondria. Nature 1981, 290, 470–474. [Google Scholar] [CrossRef] [PubMed]
  38. Ruparel, H.; Bi, L.R.; Li, Z.M.; Bai, X.P.; Kim, D.H.; Turro, N.J.; Ju, J.Y. Design and synthesis of a 3 ‘-O-allyl photocleavable fluorescent nucleotide as a reversible terminator for DNA sequencing by synthesis. Proc. Natl. Acad. Sci. USA 2005, 102, 5932–5937. [Google Scholar] [CrossRef] [Green Version]
  39. Oliveira, D.; Raychoudhury, R.; Lavrov, D.V.; Werren, J.H. Rapidly evolving mitochondrial genome and directional selection in mitochondrial genes in the parasitic wasp Nasonia (Hymenoptera: Pteromalidae). Mol. Biol. Evol. 2008, 25, 2167–2180. [Google Scholar] [CrossRef] [Green Version]
  40. Wei, S.J.; Li, Q.; van Achterberg, K.; Chen, X.X. Two mitochondrial genomes from the families Bethylidae and Mutillidae: Independent rearrangement of protein-coding genes and higher-level phylogeny of the Hymenoptera. Mol. Phylogenet. Evol. 2014, 77, 1–10. [Google Scholar] [CrossRef]
  41. Ma, Y.; Zheng, B.Y.; Zhu, J.C.; van Achterberg, C.; Tang, P.; Chen, X.X. The first two mitochondrial genomes of wood wasps (Hymenoptera: Symphyta): Novel gene rearrangements and higher-level phylogeny of the basal hymenopterans. Int. J. Biol. Macromol. 2019, 123, 1189–1196. [Google Scholar] [CrossRef]
  42. Zhu, J.C.; Tang, P.; Zheng, B.Y.; Wu, Q.; Wei, S.J.; Chen, X.X. The first two mitochondrial genomes of the family Aphelinidae with novel gene orders and phylogenetic implications. Int. J. Biol. Macromol. 2018, 118, 386–396. [Google Scholar] [CrossRef] [PubMed]
  43. Mason, W.R.M. The polyphyletic nature of Apanteles Förster (Hymenoptera: Braconidae): A phylogeny and reclassification of Microgastrinae. Mem. Entomol. Soc. Can. 1981, 113, 1–147. [Google Scholar] [CrossRef]
  44. Dowton, M.; Austin, A.D. Phylogenetic Relationships among the Microgastroid Wasps (Hymenoptera: Braconidae): Combined Analysis of 16S and 28S rDNA Genes and Morphological Data. Mol. Phylogenet. Evol. 1998, 10, 354. [Google Scholar] [CrossRef] [PubMed]
Genes 14 00230 g001 550
Figure 1. Maps of the mitochondrial genomes of Pa. hypsidrae (a), and Pr. immunis (b). The circle shows the gene map (PCGs, rRNAs, and tRNAs). The genes shown outside the map are encoded on the majority strand (J strand), while the genes represented within the map are encoded on the minority strand (N strand).
Figure 1. Maps of the mitochondrial genomes of Pa. hypsidrae (a), and Pr. immunis (b). The circle shows the gene map (PCGs, rRNAs, and tRNAs). The genes shown outside the map are encoded on the majority strand (J strand), while the genes represented within the map are encoded on the minority strand (N strand).
Genes 14 00230 g001
Genes 14 00230 g002 550
Figure 2. Base composition of protein-coding genes. (a) Pa. hypsidrae; (b) Pr. immunis.
Figure 2. Base composition of protein-coding genes. (a) Pa. hypsidrae; (b) Pr. immunis.
Genes 14 00230 g002
Genes 14 00230 g003 550
Figure 3. Gene rearrangement in the mitochondrial genome of Pa. hypsidrae and Pr. immunis.
Figure 3. Gene rearrangement in the mitochondrial genome of Pa. hypsidrae and Pr. immunis.
Genes 14 00230 g003
Genes 14 00230 g004 550
Figure 4. Phylogenetic analyses of Microgastrinae and two outgroups based on nucleotide datasets of 13 PCGs. Numbers separated by a slash on the node represent posterior probability (PP) (left) and bootstrap value (BV) (right). The graphic depicts the GenBank accession numbers for all species.
Figure 4. Phylogenetic analyses of Microgastrinae and two outgroups based on nucleotide datasets of 13 PCGs. Numbers separated by a slash on the node represent posterior probability (PP) (left) and bootstrap value (BV) (right). The graphic depicts the GenBank accession numbers for all species.
Genes 14 00230 g004
Table
Table 1. Base composition of Pa. hypsidrae and Pr. immunis.
Table 1. Base composition of Pa. hypsidrae and Pr. immunis.
SpeciesWhole GenomeProtein-Coding Genes
Length (bp)A + T (%)Length (bp)A + T (%)AT-SkewGC-Skew
Pa. hypsidrae17,06386.1511,31185.12−0.06960.1266
Pr. immunis16,39786.5611,08885.52−0.10400.0847
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Xiao, D.; Wang, Z.; Zhu, J.; Zhou, X.; Tang, P.; Chen, X. The Mitochondrial Genomes of Two Parasitoid Wasps Protapanteles immunis and Parapanteles hyposidrae (Hymenoptera: Braconidae) with Phylogenetic Implications and Novel Gene Rearrangements. Genes 2023, 14, 230. https://doi.org/10.3390/genes14010230

AMA Style

Xiao D, Wang Z, Zhu J, Zhou X, Tang P, Chen X. The Mitochondrial Genomes of Two Parasitoid Wasps Protapanteles immunis and Parapanteles hyposidrae (Hymenoptera: Braconidae) with Phylogenetic Implications and Novel Gene Rearrangements. Genes. 2023; 14(1):230. https://doi.org/10.3390/genes14010230

Chicago/Turabian Style

Xiao, Dandan, Ziqi Wang, Jiachen Zhu, Xiaogui Zhou, Pu Tang, and Xuexin Chen. 2023. "The Mitochondrial Genomes of Two Parasitoid Wasps Protapanteles immunis and Parapanteles hyposidrae (Hymenoptera: Braconidae) with Phylogenetic Implications and Novel Gene Rearrangements" Genes 14, no. 1: 230. https://doi.org/10.3390/genes14010230

APA Style

Xiao, D., Wang, Z., Zhu, J., Zhou, X., Tang, P., & Chen, X. (2023). The Mitochondrial Genomes of Two Parasitoid Wasps Protapanteles immunis and Parapanteles hyposidrae (Hymenoptera: Braconidae) with Phylogenetic Implications and Novel Gene Rearrangements. Genes, 14(1), 230. https://doi.org/10.3390/genes14010230

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop