Complete mitochondrial genomes of four entomopathogenic nematode species of the genus Steinernema
© The Author(s). 2016
Received: 30 May 2016
Accepted: 26 July 2016
Published: 5 August 2016
Nematodes belonging to the genus Steinernema are insect parasites and are used as effective biological agents against soil-dwelling insect pests. Although the full nuclear genomes of multiple Steinernema species have become available recently, mitochondrial genome information for the genus is limited. In this study, we sequenced the complete mitochondrial genomes of four species of Steinernema and analysed their structure, codon usage and phylogenetic relationships.
Mitochondrial genomes of Steinernema carpocapsae, S. glaseri, S. kushidai and S. litorale comprised 13,924, 13,851, 15,182 and 21,403 bp, respectively, with highly AT-rich nucleotide contents (AT ratio of 71.05–76.76 %). All the expected genes, including 12 protein-coding genes (encoding ATP6, CYTB, COX1-3, ND1-6 and ND4L), two rRNA genes and 22 tRNA genes were identified in the four genomes. Phylogenetic analyses based on the amino acid sequences of the 12 protein-coding genes identified the Steinernema species as monophyletic, representing a sister clade of Rhabditina and Ascaridida. In addition, they were more closely positioned to other Clade 10 nematodes, including Bursaphelenchus xylophilus, Aphelenchoides besseyi and Panagrellus redivivus, than to Strongyloides species. Gene arrangements and codon usage analyses supported this relationship. Mitochondrial genome comparison of two distinct strains of S. carpocapsae detected high intra-specific diversity.
The mitochondrial genomes of four species of Steinernema determined in this study revealed inter- and intra-species divergences/diversities of mitochondrial genomes in this genus. This information provides useful insights into the phylogenetic position of the genus Steinernema within the Nematoda and represents a useful resource for selecting molecular markers for diagnosis and population studies. These data will increase our understanding of the interesting biology of insect parasites.
Nematodes belonging to the genus Steinernema are insect parasites that are used as biological agents to control soil-dwelling insect pests [1, 2]. Steinernema nematodes form symbiotic associations with enterobacteria of the genus Xenorhabdus and the ability to kill insects is formed by a complex association of the nematodes and the bacteria, which makes the nematodes and the bacteria an attractive model for studying animal-microbe symbiosis [3, 4].
Steinernema nematodes belong to infraorder Panagrolaimomorpha, which includes, among others, the vertebrate-parasitic genus Strongyloides and free-living bacteriovores of the genera Panagrolaimus and Panagrellus . Phylogenetic reconstruction based on the small subunit ribosomal RNA gene placed the Steinernema spp. in Clade IV  and, in a later study, in Clade 10  with the Panagrolaimomorpha nematodes and the plant-parasitic/fungivorous order Aphelenchida, which includes Bursaphelenchus xylophilus and Aphelenchoides besseyi.
The number of species in the genus has increased rapidly in recent years because of extensive surveys around the world that aimed to identify more efficient biological control agents (for examples, see [8–10]). Currently, about 80 described species can be found in the genus in the NCBI taxonomy database (http://www.ncbi.nlm.nih.gov/Taxonomy/). Considering the great diversity of hosts (insects), further specific diversity in this nematode group is expected.
Mitochondrial genome sequences have been widely used in population or evolutionary studies of nematodes. Mitochondrial genomes provide a rich source of molecular variation and widespread utility in population genetics and evolutionary biology . Furthermore, mitochondria are deeply involved in many biological processes, from energy production to programmed cell death and ageing . Studying mitochondrial genomes is, therefore, essential for our understanding of fundamental biology of eukaryotes.
Recently, the full nuclear genomes of five Steinernema species were sequenced , which paved the way to investigate the population structures and genetic factors involved in their key biological processes and parasitism using genome-wide analyses. By contrast, mitochondrial genome information was limited to only one species, S. carpocapsae . Complete mitochondrial genome sequences from other Steinernema species will certainly be an important resource to study their ecology and parasitic biology. In this study, we sequenced the complete mitochondrial genomes of four species of Steinernema and analysed their structure, codon usage and phylogenetic relationships.
Steinernema carpocapsae (strain All), Steinernema glaseri (strain Sds102), Steinernema kushidai (strain Hamakita) and Steinernema litorale (strain IbKt142) were maintained using insect larvae (Galleria mellonella or Anomala cuprea) at the laboratory of plant nematology of Meiji University or Kyushu Okinawa Agricultural Research Center, National Agriculture and Food Research Organization. Infective third-stage juveniles of the nematodes were isolated from the infected insect cadavers using a White trap and further purified by the Berman funnel method, as described previously [15, 16].
Library preparation and sequencing
Genomic DNA was extracted from the nematodes using a QIAamp DNA Mini Kit (Qiagen, Tokyo, Japan) according to the manufacturer’s instructions. Illumina libraries were constructed using a Nextera DNA Sample Prep Kit (Illumina, San Diego, USA) with 100 ng of DNA, according to the manufacturer’s instructions. The libraries were sequenced on an Illumina MiSeq using a v3 Reagent kit (600 cycles), according to the manufacturer’s recommended protocol (https://icom.illumina.com/) to produce 300-bp paired-end reads.
Mitochondrial genome assembly
Mitochondrial genomes were reconstructed from the Illumina reads using MITObim ver. 1.6 . Initial assemblies were generated using SGA assembler  and mitochondrial fragments in the assembly were identified by BlastX using Caenorhabditis elegans mitochondrial genes as queries. Those fragments were extended by iterative mapping of Illumina short reads using MITObim . Gap regions in the assemblies were PCR amplified using Tks Gflex DNA Polymerase (Takara, Shiga, Japan) and sequenced on an ABI 3130 sequencer with a BigDye Terminator v3.1 Cycle Sequencing Kit (Applied Biosystems, Foster City, USA) or Illumina MiSeq, as mentioned above, to obtain a complete mitochondrial genome sequence. The assembled mitochondrial genomes were annotated for protein-coding, tRNA and rRNA genes using the MITOS web server  and by manual curation, with support from sequence similarity to other published mitochondrial genomes using Artemis genome annotation tool . Repeats in the assemblies were detected using Tandem Repeat Finder ver.4.09 . General statistical values were calculated using R (ver3.1.1) and in-house Python scripts.
Ribosomal RNA sequencing
Full-length 18S ribosomal RNA genes were amplified using primers 988F (5′-CTC AAA GAT TAA GCC ATG C-3′) and 2646R (5′-GCT ACC TTG TTA CGA CTT TT-3′) . The PCR products were purified using an innuPREP PCRpure kit (Analytik Jena, Jena, Germany) and sequenced on the ABI 3130 sequencer with the BigDye Terminator v3.1 Cycle Sequencing Kit (Applied Biosystems).
Thirty-nine nematode mitochondrial genomes were selected and included in the analysis, with an emphasis on species from Holterman’s Clade 10, to which Steinernema nematodes belong. Amino acid alignments for 12 protein-coding genes were generated separately using MAFFT , with options (−L-INS-i), and trimmed by Gblocks  with less stringent options to remove non-well-aligned sites. The best substitution model for each alignment was estimated using ProtTest . All the protein alignments were concatenated and maximum likelihood trees were constructed using RaxML ver. 7.2.8 , under the MtArt model, for each gene partition. The tree topologies obtained from ML analyses were evaluated with 1,000 bootstrap pseudoreplications. Bayesian inference was performed using MrBayes 3.2.2  from 0.5 million Markov Chain Monte Carlo generations, under a strict clock model using the Mtrev model for each gene partition. Markov chains were sampled at intervals of 100 generations. The first 0.1 million generations were discarded as ‘burn-in’ and three independent Markov Chain Monte Carlo runs converged to the same posterior probability.
Sequences of the small subunit ribosomal RNA gene were aligned using MAFFT version 7 , with option E-INS-I, and a maximum likelihood tree was obtained using RaxML ver. 7.2.8  under the GTR-gamma model with 500 bootstrap pseudoreplications.
Mitochondrial genome assembly
General statistics for four Steinernema mitochondrial genomes
Number of protein coding genes
Length of protein coding genes (bp)
Number of tRNAs
tRNA length (bp)
Number of rRNAs
rRNA length (bp)
Non-coding region length (bp)
AT% protein coding region
AT% non-coding region
Divergence of encoded proteins in the four Steinernema species. Pairwise p-distances were calculated using the Distmat program at the amino acid level
Mean distance %)
Maximum distance (%)
Minimum distance (%)
Codon usages of the four mitochondrial genomes are shown in Additional file 1: Table S1. In S. glaseri, S. litorale and S. kushidai, ATT was the most frequently used start codon (seven, five and eight cases, respectively) followed by TTG, ATG and ATA. Two GTTs in S. kushidai were the only start codons that were exceptions to the translation table (transl_table 5; invertebrate mitochondrial code). By contrast, in S. carpocapsae, exceptional start codons were used frequently (TTT for six genes and TTA for two genes). ATT and ATA start codons were observed in only two genes each in S. carpocapsae. TAA and TAG were used as termination codons in many of the genes. A truncated stop codon T was also observed in a small number of genes in all four species. The most frequently used codon in all four species was TTT (phe) followed by ATT (Ile) and TTA (leu). Principal component analysis (PCA) suggested that the overall codon usage patterns were similar in the four species compared with other Clade10 nematodes (Additional file 2: Figure S1). Within the Steinernema, the pattern of S. carpocapsae was similar to that of S. glaseri, whereas S. litorale showed a similar pattern to S. kushidai (Additional file 2: Figure S1).
The S. glaseri mitochondrial genome had a long non-coding region (1,525 bp). The non-coding region was highly AT rich (AT ratio of 85.71 %) and contained repeat sequences. Tandem repeat finders identified a short tandem repeat (112 copies of “AT”) and a long repeat (two copies of the 458-bp unit) in the sequence (Fig. 1, Additional file 1: Table S2). Two long non-coding sequences were found in the S. litorale mitochondrial genome: one between tRNA-Ser2 and tRNA-Asn (1,249 bp), and the other between tRNA-Asn and tRNA-Tyr (6,260 bp) (Fig. 1). The latter region contained many repeat sequences with variable unit sizes, with the longest unit being 383 bp (Additional file 1: Table S2).
Within the Steinernema clade, S. carpocapsae occupied the basal position and S. glaseri plus S. kushidai was at the most derived position. The relationships of the four Steinernema species inferred from the nuclear 18S rRNA gene were slightly different (Additional file 2: Figure S2). In the 18S rRNA gene tree, S. litorale and S. kushidai were positioned most internally, although S. carpocapsae was placed at the base of the four species.
Distributions of single nucleotide variations and indels in S. carpocapsae mitochondrial gene
Although the full nuclear genomes of five Steinernema species were published recently , information on the mitochondrial genomes of this genus was limited. In this study, we sequenced four complete mitochondrial genomes using next-generation sequencing techniques. To the best of our knowledge, this is the first study employing multi-species comparison of mitochondrial genomes in this genus.
Inconsistencies between mitochondrial and nuclear rRNA gene phylogenies of Nematoda have been reported previously [27, 28]. The mitochondrial gene tree obtained in this study showed a similar topology to those in previous studies. For example, the non-monophyly of Clade III (sensu Blaxter et al. 1998 ) and the non-monophyly of Tylenchomorpha (sensu De Ley & Blaxter, 2002 ) were also observed in our tree. In our mitochondrial gene phylogeny, the four Steinernema species were clustered as a sister clade of Rhabditina plus Ascaridida, with strong support (Fig. 2). Although Steinernema and Strongyloides have been placed in the same superfamily (Strongyloidoidea), the phylogeny did not support closer relationships between the two genera. This finding agrees with recent studies [14, 28–30]. In this study, we included more Steinernema species and Clade 10 species in the phylogenetic analysis, which confirmed the distant relationship of Steinernema spp. to Strongyloides spp.
Comparisons of codon usage or gene arrangements are useful to obtain insights into phylogenetic relationships. The codon usage in Steinernema species was distinguishable from that of other related species in the PCA plot (Additional file 2: Figure S1). The gene arrangements suggested that Steinernema is more closely related to Rhabditina and Ascaridida species than to Strongyloides species (Fig. 3). These results are consistent with the phylogeny based on amino acid data. However, it is not straightforward to interpret the relationships of the four species in the genus. The codon usage analysis placed S. glaseri and S. carpocapsae in a sub-group (Additional file 2: Figure S1); however, the gene arrangements of S. glaseri were more similar to S. litorale than to S. carpocapsae (Fig. 3). Either of these results are not completely consistent with mitochondrial amino acid phylogeny and 18S phylogeny (Fig. 2; Additional file 2: Figure S2).
Comparison of two strains of S. carpocapsae of different origins revealed 1.5 % variations in the mitochondrial genomes at the nucleotide level, indicating that the mitochondrial genome can provide good markers for population studies. We observed a location-biased distribution of SNVs and indels, suggesting that marker selection will be important for population studies. Interestingly, the rRNA gene regions have fewer SNVs than the protein-coding region, while indels were more frequently represented in the rRNA gene and were not found in protein-coding regions. This is probably because of functional differences between protein-coding and rRNAs. Indels in protein-coding region can be more problematic because they could cause large protein changes by insertions of stop or frame-shift mutations, whereas mismatches by SNVs, rather than indels, could be critical in rRNA genes because they could produce more erratic RNA secondary structures. These are only the results from a comparison of a small number of strains. For a deep understanding of their evolution patterns, it will be interesting to re-sequence and compare more mitochondrial genomes from a larger numbers of strains using high-throughput sequencers. The genome information provided in this study will be useful as references in those kinds of studies.
In this study, the complete mitochondrial genome sequences of four Steinernema species were determined. The complete mitochondrial genomes of S. carpocapsae, S. glaseri, S. kushidai and S. litorale were 13,924, 13,851, 15,182 and 21,403 bp in length, respectively, and all contained the 36 expected genes (12 protein-coding genes, two ribosomal RNA genes and 22 transfer RNA genes), all transcribed in the same direction. Phylogenetic analyses using the amino acid sequences of the 12 protein-coding genes showed that the Steinernema species formed a monophylic clade as a sister clade of Rhabditina and Asacridida. They were also more closely positioned to other Clade 10 nematodes, including B. xylophilus, A. besseyi and P. redivivus, than to Strongyloides species, which belong to the same infraorder as the Steinernema species. The data on codon usage and gene arrangement also supported this result. We also compared the mitochondrial genomes of two strains of S. carpocapsae and detected high intraspecific diversity. This mitochondrial genome information will form a useful resource for evolutionary and population studies of Steinernema nematodes.
The authors thank Dr. Nobuo Ogura for kindly providing nematode strains.
This study was partly supported by the Japan Society for the Promotion of Science KAKENHI (Grant Nos. 16H04722, 16 K15267 and 26292178).
Availability of data and material
The datasets supporting the conclusions of this article are available in the Genome sequence DDBJ/EMBL/GenBank international nucleotide sequence repository under accession numbers AP017463–AP017468.
TK designed the experiments, collected the data and performed the analysis. TA performed the 18S analysis. MY prepared the biological materials. TK and TA wrote the manuscript. All authors read and approved the final version of the manuscript.
The authors declare that they have no competing interests.
Consent for publication
Ethics approval and consent to participate
Open AccessThis article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.
- Burnell AM, Stock SP. Heterorhabditis, Steinernema and their bacterial symbionts lethal pathogens of insects. Nematology. 2000;2(1):31–42.View ArticleGoogle Scholar
- Lacey LA, Georgis R. Entomopathogenic nematodes for control of insect pests above and below ground with comments on commercial production. J Nematol. 2012;44(2):218–25.PubMedPubMed CentralGoogle Scholar
- Murfin KE, Dillman AR, Foster JM, Bulgheresi S, Slatko BE, Sternberg PW, et al. Nematode-bacterium symbioses - Cooperation and conflict revealed in the ‘omics’ age. Biol Bull. 2012;223(1):85–102.PubMedPubMed CentralGoogle Scholar
- Martens EC, Heungens K, Goodrich-Blair H. Early colonization events in the mutualistic association between Steinernema carpocapsae nematodes and Xenorhabdus nematophila bacteria. J Bacteriol. 2003;185(10):3147–54.View ArticlePubMedPubMed CentralGoogle Scholar
- De Ley P, Blaxter M. Systematic Position and Phylogeny. In: Lee DL, editor. The Biology of Nematodes. London: Taylor & Francis; 2002. p. 1–30.View ArticleGoogle Scholar
- Blaxter ML, De Ley P, Garey JR, Liu LX, Scheldeman P, Vierstraete A, et al. A molecular evolutionary framework for the phylum Nematoda. Nature. 1998;392(6671):71–5.View ArticlePubMedGoogle Scholar
- Holterman M, van der Wurff A, van den Elsen S, van Megen H, Bongers T, Holovachov O, et al. Phylum-wide analysis of SSU rDNA reveals deep phylogenetic relationships among nematodes and accelerated evolution toward crown Clades. Mol Biol Evol. 2006;23(9):1792–800.View ArticlePubMedGoogle Scholar
- San-Blas E, Portillo E, Nermut J, Puza V, Morales-Montero P. Steinernema papillatum n. sp. (Rhabditida: Steinernematidae), a new entomopathogenic nematode from Venezuela. Nematology. 2015;17(9):1081–97.View ArticleGoogle Scholar
- Ma J, Chen S, De Clercq P, Waeyenberge L, Han R, Moens M. A new entomopathogenic nematode, Steinernema xinbinense n. sp. (Nematoda: Steinernematidae), from north China. Nematology. 2012;14(6):723–39.View ArticleGoogle Scholar
- Cimen H, Lee M-M, Hatting J, Hazir S, Stock SP. Steinernema tophus sp. n. (Nematoda: Steinernematidae), a new entomopathogenic nematode from South Africa. Zootaxa. 2014;3821(3):337–53.View ArticlePubMedGoogle Scholar
- Galtier N, Nabholz B, GlÉMin S, Hurst GDD. Mitochondrial DNA as a marker of molecular diversity: a reappraisal. Mol Ecol. 2009;18(22):4541–50.View ArticlePubMedGoogle Scholar
- Balaban RS, Nemoto S, Finkel T. Mitochondria, oxidants, and aging. Cell. 2005;120(4):483–95.View ArticlePubMedGoogle Scholar
- Dillman AR, Macchietto M, Porter CF, Rogers A, Williams B, Antoshechkin I, et al. Comparative genomics of Steinernema reveals deeply conserved gene regulatory networks. Genome Biol. 2015;16:200.View ArticlePubMedPubMed CentralGoogle Scholar
- Montiel R, Lucena MA, Medeiros J, Simoes N. The complete mitochondrial genome of the entomopathogenic nematode Steinernema carpocapsae: insights into nematode mitochondrial DNA evolution and phylogeny. J Mol Evol. 2006;62(2):211–25.View ArticlePubMedGoogle Scholar
- White GF. A method for obtaining infective nematode larvae from cultures. Science. 1927;66(1709):302–3.View ArticlePubMedGoogle Scholar
- Ogura N, Nakashima T. Cold tolerance and preconditioning of infective juveniles of Steinernema kushidai (Nematoda: Steinernematidae). Nematologica. 1997;43(1):107–15.View ArticleGoogle Scholar
- Hahn C, Bachmann L, Chevreux B. Reconstructing mitochondrial genomes directly from genomic next-generation sequencing reads-baiting and iterative mapping approach. Nucleic Acids Res. 2013;41(13):e129.View ArticlePubMedPubMed CentralGoogle Scholar
- Simpson JT, Durbin R. Efficient de novo assembly of large genomes using compressed data structures. Genome Res. 2012;22:549–56.View ArticlePubMedPubMed CentralGoogle Scholar
- Bernt M, Donath A, Juhling F, Externbrink F, Florentz C, Fritzsch G, et al. MITOS: Improved de novo metazoan mitochondrial genome annotation. Mol Phylogenet Evol. 2013;69(2):313–9.View ArticlePubMedGoogle Scholar
- Carver T, Harris SR, Berriman M, Parkhill J, McQuillan JA. Artemis: an integrated platform for visualization and analysis of high-throughput sequence-based experimental data. Bioinformatics. 2012;28(4):464–9.View ArticlePubMedGoogle Scholar
- Benson G. Tandem repeats finder: a program to analyze DNA sequences. Nucleic Acids Res. 1999;27(2):573–80.View ArticlePubMedPubMed CentralGoogle Scholar
- Katoh K, Misawa K, Kuma K, Miyata T. MAFFT: a novel method for rapid multiple sequence alignment based on fast Fourier transform. Nucleic Acids Res. 2002;30(14):3059–66.View ArticlePubMedPubMed CentralGoogle Scholar
- Talavera G, Castresana J. Improvement of phylogenies after removing divergent and ambiguously aligned blocks from protein sequence alignments. Syst Biol. 2007;56(4):564–77.View ArticlePubMedGoogle Scholar
- Abascal F, Zardoya R, Posada D. ProtTest: selection of best-fit models of protein evolution. Bioinformatics. 2005;21(9):2104–5.View ArticlePubMedGoogle Scholar
- Stamatakis A. RAxML-VI-HPC: maximum likelihood-based phylogenetic analyses with thousands of taxa and mixed models. Bioinformatics. 2006;22(21):2688–90.View ArticlePubMedGoogle Scholar
- Huelsenbeck JP, Ronquist F. MRBAYES: Bayesian inference of phylogenetic trees. Bioinformatics. 2001;17(8):754–5.View ArticlePubMedGoogle Scholar
- Hunt VL, Tsai IJ, Coghlan A, Reid AJ, Holroyd N, Foth BJ, et al. The genomic basis of parasitism in the Strongyloides clade of nematodes. Nat Genet. 2016;48(3):299–307.View ArticlePubMedGoogle Scholar
- Park J-K, Sultana T, Lee S-H, Kang S, Kim HK, Min G-S, Eom KS, Nadler SA. Monophyly of clade III nematodes is not supported by phylogenetic analysis of complete mitochondrial genome sequences. BMC Genomics. 2011;12(1):1.View ArticleGoogle Scholar
- Kim T, Kim J, Nadler SA, Park J-K. The complete mitochondrial genome of Koerneria sudhausi (Diplogasteromorpha: Nematoda) supports monophyly of Diplogasteromorpha within Rhabditomorpha. Curr Genet. 2016;62(2):391–403.View ArticlePubMedGoogle Scholar
- Kim J, Lee S, Gazi M, Kim T, Jung D, Chun J, Kim S, Seo T, Park C, Baldwin JG. Mitochondrial genomes advance phylogenetic hypotheses for Tylenchina (Nematoda: Chromadorea). Zool Scripta. 2015;44(4):446–62.View ArticleGoogle Scholar