Outline
Comptes Rendus

Molecular biology and genetics/Biologie et génétique moléculaires
Fast assembly of the mitochondrial genome of a plant parasitic nematode (Meloidogyne graminicola) using next generation sequencing
Comptes Rendus. Biologies, Volume 337 (2014) no. 5, pp. 295-301.

Abstracts

Little is known about the variations of nematode mitogenomes (mtDNA). Sequencing a complete mtDNA using a PCR approach remains a challenge due to frequent genome reorganizations and low sequence similarities between divergent nematode lineages. Here, a genome skimming approach based on HiSeq sequencing (shotgun) was used to assemble de novo the first complete mtDNA sequence of a root-knot nematode (Meloidogyne graminicola). An AT-rich genome (84.3%) of 20,030 bp was obtained with a mean sequencing depth superior to 300. Thirty-six genes were identified with a semi-automated approach. A comparison with a gene map of the M. javanica mitochondrial genome indicates that the gene order is conserved within this nematode lineage. However, deep genome rearrangements were observed when comparing with other species of the superfamily Hoplolaimoidea. Repeat elements of 111 bp and 94 bp were found in a long non-coding region of 7.5 kb, as similarly reported in Mjavanica and Mhapla. This study points out the power of next generation sequencing to produce complete mitochondrial genomes, even without a reference sequence, and possibly opening new avenues for species/race identification, phylogenetics and population genetics of nematodes.

Supplementary Materials:
Supplementary material for this article is supplied as a separate file:

Il existe peu de connaissances sur la variation des génomes mitochondriaux de nématodes, car leur séquençage avec une approche PCR reste difficile du fait de fréquentes réorganisations génomiques et de la faible similitude de séquence entre des lignées divergentes. Nous avons utilisé une approche d’écrémage de génome pour assembler la première séquence mitochondriale complète d’un nématode à galles (Meloidogyne graminicola). Un génome riche en AT (84,3 %) de 20 030 paires de bases (pb) a été obtenu avec une profondeur de séquençage supérieure à 300. Trente-six gènes ont été identifiés selon une approche semi-automatisée. La comparaison de cette séquence avec la carte physique du génome mitochondrial de M. javanica indique que l’ordre des gènes est conservé chez Meloidogyne. Toutefois, des réarrangements génomiques importants ont été observés entre genres de la superfamille Hoplolaimoidea. Des éléments répétés de 111 et 94 pb ont été détectés dans une longue région non codante d’environ 7500 pb, comme rapporté chez Mjavanica et Mhapla. Cette étude montre le potentiel du séquençage de nouvelle génération pour produire rapidement des génomes mitochondriaux complets, même en l’absence de séquence de référence. Ceci ouvre des perspectives dans les domaines de l’identification d’espèces/races, de la phylogénie et de la génétique des populations de nématodes.

Compléments :
Des compléments sont fournis pour cet article dans le fichier séparé :

Metadata
Received:
Accepted:
Published online:
DOI: 10.1016/j.crvi.2014.03.003
Keywords: Organellar genome assembly, Genome skimming, Meloidogyne, Mitogenomics, Nematode, Next generation sequencing
Keywords: Assemblage de génome d’organelles, Écrémage de génome, Meloidogyne, Mitogénome, Nématode, Séquençage de nouvelle génération

Guillaume Besnard 1; Frank Jühling 2; Élodie Chapuis 3; Loubab Zedane 1; Émeline Lhuillier 4; Thierry Mateille 5; Stéphane Bellafiore 3

1 CNRS–UPS–ENFA, UMR5174, EDB (Laboratoire Évolution & Diversité biologique), 118, route de Narbonne, 31062 Toulouse, France
2 University of Leipzig, Interdisciplinary Center for Bioinformatics, Härtelstraße 16–18, 04107 Leipzig, Germany
3 IRD, UMR186 “Résistance des plantes aux bioagresseurs”, 911, avenue Agropolis, BP 64501, 34394 Montpellier cedex 5, France
4 INRA, GeT–PlaGe, UAR 1209 Département de génétique animale, INRA Auzeville, 31326 Castanet-Tolosan, France
5 IRD, UMR 022 Centre de biologie pour la gestion des populations, campus de Baillarguet, CS30016, 34988 Montferrier-sur-Lez cedex, France
@article{CRBIOL_2014__337_5_295_0,
     author = {Guillaume Besnard and Frank J\"uhling and \'Elodie Chapuis and Loubab Zedane and \'Emeline Lhuillier and Thierry Mateille and St\'ephane Bellafiore},
     title = {Fast assembly of the mitochondrial genome of a plant parasitic nematode {(\protect\emph{Meloidogyne} graminicola}) using next generation sequencing},
     journal = {Comptes Rendus. Biologies},
     pages = {295--301},
     publisher = {Elsevier},
     volume = {337},
     number = {5},
     year = {2014},
     doi = {10.1016/j.crvi.2014.03.003},
     language = {en},
}
TY  - JOUR
AU  - Guillaume Besnard
AU  - Frank Jühling
AU  - Élodie Chapuis
AU  - Loubab Zedane
AU  - Émeline Lhuillier
AU  - Thierry Mateille
AU  - Stéphane Bellafiore
TI  - Fast assembly of the mitochondrial genome of a plant parasitic nematode (Meloidogyne graminicola) using next generation sequencing
JO  - Comptes Rendus. Biologies
PY  - 2014
SP  - 295
EP  - 301
VL  - 337
IS  - 5
PB  - Elsevier
DO  - 10.1016/j.crvi.2014.03.003
LA  - en
ID  - CRBIOL_2014__337_5_295_0
ER  - 
%0 Journal Article
%A Guillaume Besnard
%A Frank Jühling
%A Élodie Chapuis
%A Loubab Zedane
%A Émeline Lhuillier
%A Thierry Mateille
%A Stéphane Bellafiore
%T Fast assembly of the mitochondrial genome of a plant parasitic nematode (Meloidogyne graminicola) using next generation sequencing
%J Comptes Rendus. Biologies
%D 2014
%P 295-301
%V 337
%N 5
%I Elsevier
%R 10.1016/j.crvi.2014.03.003
%G en
%F CRBIOL_2014__337_5_295_0
Guillaume Besnard; Frank Jühling; Élodie Chapuis; Loubab Zedane; Émeline Lhuillier; Thierry Mateille; Stéphane Bellafiore. Fast assembly of the mitochondrial genome of a plant parasitic nematode (Meloidogyne graminicola) using next generation sequencing. Comptes Rendus. Biologies, Volume 337 (2014) no. 5, pp. 295-301. doi : 10.1016/j.crvi.2014.03.003. https://comptes-rendus.academie-sciences.fr/biologies/articles/10.1016/j.crvi.2014.03.003/

Version originale du texte intégral

1 Introduction

Around one billion years old, Nematoda (round worms) constitute the most ancient and diverse animal phylum [1], with estimates of the number of extant species ranging from 25,000 to certainly more than 100 million [2]. As free living organisms as well as parasites of animals and plants, these species represent an ecologically diverse assemblage of roundworms distributed in nearly every environment, including soil and sediments, fresh and sea water [3]. Five main lineages were identified using molecular data, and all diverged from a shared ancestor several hundred million years ago [4]. Today parasitic nematodes are recognized to cause significant impacts on agriculture, domestic animals and human health [5].

Endoparasitic root-knot nematodes (Meloidogyne spp.) are foremost economically damaging genera of plant parasitic nematodes [6]. Meloidogyne spp. belongs to the class Secernentea, order Tylenchida, suborder Hoplolaimina, superfamily Hoplolaimoidea, family Meloidogynidae, subfamily Meloidogyninae. Recent phylogenies indicated that Meloidogyne is a monophyletic group but these analyses also pointed to poorly resolved paraphyletic relationships with the sister lineage that includes genera Pratylenchus, Hirschmanniella and Zygotylenchus [7,8]. The Meloidogyne genus includes more than 60 species, while some of the members even have several races. Meloidogyne species are major pests worldwide with very large plant host ranges (e.g. MjavanicaM. arenariaMincognita and Mhapla) and cause important crop losses locally, as it has been reported for M. graminicola in rice fields [9].

Meloidogyne species use various reproduction strategies (i.e. amphimixis, facultative meiotic parthenogenesis or obligatory mitotic parthenogenesis [10,11]). For instance, Mincognita, Mjavanica and Marenaria have a mitotic partenogenetic reproduction mode whereas Mhapla and Mgraminicola reproduce usually by a facultative meiotic parthenogenetic mode and less so frequently by cross fertilization [10]. Although these nematodes are particularly well studied because of their economic importance and their ecological widespread adaptation, their characterization faces constraints. Species and race identification of Meloidogyne based on phenotypic characterization is time-consuming because adult forms are hardly accessible since males are usually rare and females are distorted (pear-shaped due to their settlement inside roots). Therefore, taxonomists use routinely morphometrics of second stage juveniles (e.g. tapering tail, body and stylet size) and female perineal patterns [12,13]. Consequently, species identification is long and requires trained nematologists. For these reasons, biochemical and molecular tools have been developed to assist the identification of taxa [13–15] as well as to investigate the phylogeography and population dynamics of nematodes (for a recent review on this topic, see [16]). However, the molecular investigation of nematodes still has to face other challenges caused by the difficulty to isolate individuals or colonies for extracting enough protein or DNA. Despite of this, isozymes (esterase) as well as ribosomal and mitochondrial DNA sequences are generally used for routine diagnostics of nematodes [12,15,17]. Polymorphic DNA markers (microsatellites) have been also developed to study the genetic structure of the potato cyst nematode Globodera pallida in a whole country [18], or to detect and locate rare sexual events in predominantly asexual Xiphinema spp. [19]. Lastly, it is worth to mention that the nuclear genome of Meloidogyne spp. is usually small (e.g., M. graminicola exhibits the smallest animal genome referenced to date, C-value = ca. 0.03 pg; [20]) and recent advances about the complete genome sequence of these organisms should provide new insights in nematode diagnostics and ecology [21].

Mitochondrial DNA (mtDNA) is usually maternally inherited in nematodes [22,23], and polymorphism of this genome is thus useful to follow maternal lineages. For instance, the cox1 mitochondrial gene has been used to study the genetic structure of the plant parasitic nematode Xiphinema spp. [24]. No reference mitochondrial genome (or mitogenome) is presently available for root-knot nematodes, although the structure of the Mjavanica mtDNA has been described more than 20 years ago, and was partially sequenced (X57625/X57626; ca. 20% of the genome; [25]). Today, in the superfamily Hoplolaimoidea, only three complete mitogenomes have been deposited in the public DNA databases [26–28] for the distantly related species Heterodera glycines (Heteroderidae), Radopholus similis, and Pratylenchus vulnus (Pratylenchidae). Thus, little is known about the variations in mitogenomes of this important phytoparasitic lineage. This also represents a limitation to develop genetic tools based on mtDNA sequence polymorphism for various purposes from species or race identification to phylogeography or population genetics. In addition, frequent mitogenome reorganizations have been reported in Hoplolaimoidea [28], rendering fastidious the development of PCR-based methods to sequence long mtDNA regions. Lastly, the gene annotation of nematode mitochondrial genomes remains a challenge, especially the typical truncated tRNA (transfer RNA) genes are hard to detect with current tools.

Recent advances in sequencing technologies, so called next generation sequencing or NGS, are bringing about a revolution for the sequencing of organellar genomes [29]. NGS has been used to generate mitochondrial genomes of various metazoan organisms, but the method is still usually based on the sequencing of long-range PCR products. The method however offers the possibility of fast assembly of the most represented DNA regions in a genome with high sequencing depths [29] but this has never been tested on nematodes. In particular, the Illumina technology (e.g. HiSeq 2000; [30]) allows generating tens of millions of short DNA segments (i.e. 100-bp reads) that can be used for this purpose.

Here, the Illumina technology was tested on M. graminicola in order to generate a reference mitogenome for root-knot nematodes. A complete mtDNA genome sequence was released and compared to the mtDNA gene map of M. javanica, and also to other species of the superfamily Hoplolaimoidea. Lastly, the utility of our approach for generating full mtDNA resources in nematodes is briefly discussed.

2 Material and methods

2.1 DNA extraction and DNA sequencing

An isolate of M. graminicola from a rice field in Laurel, Batangas, Philippines [31], was propagated in a hydroponic solution on rice cultivar IR64 during four weeks before egg extraction. Eggs were recovered from infested roots using the hypochlorite extraction method with minor modifications [32] and were collected in 0.8% hypochlorite solution with 0.1% Tween 20 then left for 30 min at room temperature. Eggs were rinsed several times with sterile ddH2O and placed on a strainer covered by three damp Kimwipe tissues on a 50-mL beaker filled with sterile ddH2O and streptomycin for a final concentration of 10 μg/mL. After 24 h in the dark at room temperature, water containing nematodes was discarded and replaced with new sterile ddH2O. Twenty-four hours later, the freshly hatched second stage juveniles (J2s) were collected by centrifugation and conserved at –80 °C. For 100 μL of J2, an equal volume of 2 × lysis buffer containing 20 mM Tris–HCl pH 8.4, 100 mM KCl, 3 mM MgCl2 and 0.9% Tween 20 was added. Fresh DTT (2 mM) and 20 μg of fresh proteinase K were added and the tube was placed in a water bath at 50 °C. After a 2-h treatment, 0.1 μg of RNase A (Qiagen Inc., GmbH, Hilden, Germany) was added and the tube was placed at 37 °C for an additional 30 min. Then, genomic DNA was recovered by a phenol–chloroform method [33]. DNA quality was inspected with a NanoDrop® spectrophotometer and 1 μL of genomic DNA (gDNA) was loaded on 1% agarose gel to control DNA integrity.

Three hundred nanograms of DNA were used for shotgun sequencing with the Illumina technology (e.g., HiSeq 2000; [30]). The library was constructed using the Illumina TruSeq DNA Sample Prep v2 kit and following the instructions of the supplier. Briefly, genomic DNA was sheared by sonication and the fragments were end-repaired. Purified fragments were A-tailed and ligated to sequencing indexed adapters. Fragments with an insert size of approximately 300 to 500 bp were gel-extracted and enriched with 10 cycles of PCR before library quantification and validation. This library was multiplexed with 23 other libraries (generated by other colleagues on various plant species). The pool of libraries was then hybridized to the Hiseq 2000 flow cell using the Illumina TruSeq PE Cluster Kit v3. Bridge amplification was performed to generate clusters, and paired reads of 100 nucleotides were collected on the HiSeq 2000 sequencer using the Illumina TruSeq SBS Kit v3 (200 cycles).

2.2 Read assembly, sequencing depth and gene annotation

The read assembly method is similar to the approach reported by [34] to generate a complete plastid genome in plants. We used the cox2 sequence of Mgraminicola as a seed because it was already available in public DNA databases (JN241929; 530 bp). This sequence allowed initiating the global assembly process that involved two main steps. First, we used the program extractread2 (included in the OBITools package, http://metabarcoing.orf/obitools) to select sequence reads sharing at least 80 bases with the seed. The newly selected reads were subsequently used as probes. This process was iteratively repeated until no new read could be identified. Second, the set of selected reads were assembled into two large contigs using the Velvet program [35]. Long repeat elements in two regions (see below) hampered us to generate one contig of the full mitogenome. All computations were done on the computer cluster of the Genopole Bioinformatic Platform (Toulouse, France).

Two regions, namely 111R and 94R, are composed of 111-bp and 94-bp repeats and could not be assembled with our approach (see below). While region 111R includes about 30 repeats and is difficult to sequence with current methodologies, 94R is relatively short (inferior to 1000 bp). This later region was generated by PCR with the following primers (forward: GAATGAAAACGATCAGAGACATG; reverse: CAACTAATATTTTTGATAGCACTCC). The 94R region was amplified by PCR in a final volume of 50 μL containing 10 ng of gDNA, 500 nM of each PCR primer, 200 μM of each dNTP, 10 μL of 5 × Phusion buffer and one unit of the proof-reading enzyme “Phusion High-Fidelity DNA Polymerase” (Thermo Fisher Scientific Inc., France). The following PCR program was used: 98 °C for 30 s; 35× (98 °C, 10 s; 55 °C, 30 s; 72 °C, 70 s); 72 °C for 10 min. PCR products were gel-purified and eluted in 30 μL ddH2O using a gel purification kit (Qiagen Inc., Germany). Approximately 200 ng of purified amplicon were sent for sequencing with a classical Sanger approach (Cogenics, UK). The generated sequence allowed merging both contigs. Finally, reads were mapped onto the obtained mtDNA sequence using Geneious v6.0.5 (Biomatters Ltd.) to assess the depth of sequencing and the number of repeats in the region 111R.

We annotated the consensus mtDNA sequence using the protein prediction pipeline of Mitos with the invertebrate mitochondrial code [36], but could not detect all common mitochondrial genes, especially rRNA (ribosomal RNA) and tRNA genes. We thus used a semi-automated procedure: all ORFs of the mitochondrial genome of Mgraminicola were extracted and translated with the help of Getorf from a local installation of the Emboss package v6.4.0.0 [37]. Candidate genes were aligned to known nematode mitochondrial proteins extracted from RefSeq 58 using a local installation of Mafft v7.058b [38], and identified manually by sequence similarity. All tRNA genes could be detected with MiTFi [39], using Infernal v.1.0 [40] and covariance models as described by [41]. The applied models additionally contained truncated tRNAs of nematodes of the class Enoplea, which are partially experimentally validated [42]. Using the nematode models without the Enoplea tRNAs, some tRNAs could not be found in M. graminicola. All tRNAs have been checked manually, the start/stop positions of trnG, trnH and trnS1 have been improved. Both rRNA genes were cut using alignments of the Tylenchida (Bursaphelenchus mucronatus, Pratylenchus vulnus, Radopholus similis) extracted from RefSeq 58. For these sequences the 5′ ends of the adjacent tRNAs were used as start position, and the 3′ end of the rrnL gene was trimmed manually to fit with the nad3 gene. The structures have been cleaned manually and all wrongly predicted base pairs, resulting from the Infernal output, were removed.

The newly generated sequence was deposited within GenBank (HG529223). A physical map was finally drawn with Mtviz (http://pacosy.informatik.uni-leipzig.de/mtviz/) and colors were modified. Gene order in the M. graminicola mtDNA was compared with the three available complete mitochondrial genomes of the chromadorean nematode clade (i.e. Heterodera glycines, Radopholus similis and Pratylenchus vulnus; [28]). In order to check the phylogenetic position of M. graminicola in the superfamily Hoplolaimoidea, the 12 protein encoding genes were extracted and aligned with the complete mitochondrial gene sequences of 41 other nematodes as used by [28]. A Maximum Likelihood analysis was conducted with PhyML 3.0 [43]. Statistical support for the different internal branches was assessed through bootstrap resampling (1000 replicates).

3 Results and discussion

3.1 De novo assembly of a root-knot nematode mitochondrial genome

We generated 6,012,732 paired-end 100-bp reads. Using extractread2, we selected 35,597 paired-end reads from the mitochondrial genome (ca. 0.6% of initial sequences). The first assembly generated two contigs of 15,109 and 1813 bp. These contigs were separated by regions with long repeats (namely 111R and 94R) that were not possible to properly assemble. These two repeats are located in the same regions than 102R and 63R in M. javanica ([25]; see below). Sequence comparisons did not show homology between repeats from the different species (see also [44] for 102R in M. hapla). The sequencing depth was estimated on the non-repeated regions at 318.5× ± 121.7 (min. 68×, max. 651×). Based on this value, we estimated that 94R is composed of three repeats of 94 bp that was confirmed by a classical Sanger sequencing. Similarly, sequencing depth of region 111R (9,077× ± 1830) suggested that about 29 tandem repeats of 111 bp composed this region. A consensus sequence of an 111R repeat was generated considering all mismatches: TATTTTTAATTATTRYYATTTTTGATGATTTTTCAGAATTTTGYTTAATTTAAYTTTATTTTTAARCCAATTTTTTGTAAAATATRAATYATAWAGGTAGTTACCTTTTYA. Considering 29 repeats in 111R, the mitochondrial genome size of M. graminicola was estimated at 20,030 bp, a size similar to the mtDNA size reported for M. javanica and M. incognita (ca. 20.5 kb and 19.5 kb, respectively; [25]). Lastly, the mtDNA of M. graminicola was especially AT-rich (84.3%) as previously reported in M. javanica and M. hapla [25,44].

Our complete sequence was also compared to mitochondrial DNA sequences of M. graminicola available in GenBank. Polymorphism was found in both genes cob and cox2: 16 substitutions/531 bp (3%) were found in cob compared to sequences of American isolates (AY757884 and AY757885) and three to five substitutions/530 bp (ca. 1%) were found in cox2 compared to six sequences of Asian isolates (from India, China and Bangladesh; JN241926 to JN241929, JN241939, GU187309, GQ865513). This result indicates that intraspecific variation is present in coding regions and potentially useful for population genetic analyses or race identification, but non-coding regions – such as 94R – may potentially be more variable (due to relaxed selective pressures compared to protein encoding genes) but their informativeness still needs to be assessed.

3.2 Mitochondrial genome organisation

Thirty-six genes were detected in the M. graminicola mtDNA sequence (Table S1, Supplementary data) as in other available mitogenomes of the suborder Hoplolaimina [28]: 12 protein encoding genes, 22 transfer RNA genes (trn) and two ribosomal RNA genes (rrn). Among the typical mitochondrial genes of metazoans, only atp8 is missing as already reported in other nematode mitogenomes [28,45]. An annotated physical map of the mtDNA genome was reconstructed and is shown in Fig. 1. We found the same gene order than in the incomplete gene map of M. Javanica [25], suggesting no important gene reorganisation between these two Meloidogyne species. In contrast, the comparison of the M. graminicola mtDNA with those of the chromadorean nematode clade reveals a deep reorganisation (Fig. 2). The gene order between rrnS and cox2 in the M. graminicola mtDNA is, however, similar to the mitogenome of Pratylenchus vulnus. In this genomic region of about 5000 bp (that includes 16 genes; Table S1, Supplementary data), only one gene is not located at the same position in these two species: trnF, which is located between trnN and trnG in M. graminicola, while it is located between trnC and nad6 in P. vulnus. Additionally, we detected in the mitogenomes of M. graminicola and P. vulnus the same three pairs of genes nad4-trnD, atp6-nad5 and trnQ-trnA plus another region with three genes (rrnL-nad3-cob). In contrast, a comparison of the mtDNA gene orders between M. graminicola and Heterodera glycines or Radopholus similis – two other members of the superfamily Hoplolaimoidea – show deep genome rearrangements. Only small regions with the same gene arrangements were detected between M. graminicola and R. similis [from cox2 to nad3 (four genes), from nad2 to trnN (four genes), from atp6 to trnA (four genes), from trnY to nad1 (three genes) and nad6-nad4L; Fig. 2], and between M. graminicola and H. glycines (from atp6 to trnA (four genes), plus three pairs of genes: trnE-trnS1, trnI-cox3 and rrnL-nad3; Fig. 2). Only four pairs of genes are conserved between the four mitogenomes of Hoplolaimoidea (i.e. atp6-nad5, trnQ-trnA, trnI-cox3 and rrnL-nad3). Note that Pratylenchus is the most related species to Meloidogyne for which a complete mitogenome is available ([7,8]; Supplementary data, Fig. S1) and the genome arrangement may reflect the phylogenetic distance between taxa [45]. The timing of mitogenome rearrangements needs to be investigated on a larger taxon sample in a phylogenetic framework. M. graminicola and P. vulnus also share a large non-coding region (that harbours only one tRNA gene). The length of this region is ca. 7.5 and 7.8 kb in M. graminicola and P. vulnus, respectively. To date, such a long non-coding region has not been reported in other mitogenomes of Hoplolaimoidea. This region is reported to contain the replication origin (the control region) that varies from a few hundred bp to tens of kb in metazoans [28].

Fig. 1

(Color online). Circular gene map of the complete mitochondrial genome of Meloidogyne graminicola. Genes and non-coding regions are represented as boxes. The direction of transcription of all genes is the same and indicated by the arrow. Abbreviations of protein coding and rRNA genes are: NADi = subunit i of NADH dehydrogenase; COXi = subunit i of cytochrome c oxydase; COB = cytochrome b; ATP6 = subunit 6 of ATP-synthase, 12S  =  12S ribosomal RNA, 16S = 16S ribosomal RNA. tRNA genes are named with single-letter amino acid abbreviations except for those coding for leucine and serine, which are named as L1 (anticodon uag), L2 (uaa), S1 (ucu), and S2 (uga). Two regions, namely 111R and 94R, are composed of 111-bp and 94-bp repeats. These elements are located in a large AT-rich non-coding region (that harbours only one tRNA gene, trnV) of 7.5 kb (in grey).

Fig. 2

(Color online). Linearized representation of the mitochondrial gene arrangement in the four available mitogenomes of the superfamily Hoplolaimoidea. Gene and genome size are not to scale. All genes are transcribed in the same direction (from left to right). The tRNAs are designated by single-letter abbreviation (in circles). The AT-rich non-coding regions are not indicated.

3.3 Structures of tRNA and rRNA genes

The short sizes of tRNAs, ranging between 45 and 56 bp, and the fact that a part of them could only be detected by including truncated Enoplea sequences to the covariance models, indicates how minimal the mitochondrial tRNAs of M. graminicola are. Truncated mitochondrial tRNA sequences are found throughout nearly all metazoan phyla [41]. While nearly all nematodes share common truncated tRNA sequences missing either the D- or the T-stem, e.g., Caenorhabditis elegans and Ascaris suum [46], the most truncated sequences were found in Enoplea [41]. Here, tRNAs were found that lack both, the D- and T-stems. All tRNAs in M. graminicola share the typical nematode structure missing the T-stem, but forming a D-stem (data not shown). The only exceptions from this rule are serine tRNAs, which instead miss the D-stem in nematodes. This is also the case for trnS1 and trnS2 in M. graminicola. Both rRNA structures show characteristics that are common for nematodes, with the typical domains. Their short lengths (602/830 bp) again support the trend to minimal structures as related organisms encode larger rRNAs, e.g., Radopholus similis (692/840 bp). Both rRNA structures are shown on Supplementary data, Fig. S2.

The automated annotation of mitochondrial genomes is still a challenging problem. All highly conserved genes could be predicted by the Mitos pipeline [36]. However, the programs failed to predict the truncated non-coding RNAs correctly. We expect that our findings will help to improve these pipelines as most of predictive programs are based on former annotations, which is even more important for truncated sequences as it is the case in nematodes, and especially in the M. graminicola genome. Nice examples are the mitochondrial tRNAs of the American house dust mites. Here, an improvement in the detection methods led to the discovery of extremely short tRNAs [47].

4 Concluding remarks

Complete mitogenomes of nematodes that were recently released in DNA databases were generated with PCR-based methods [28], and this approach can be challenging due to numerous genome rearrangements. Our study indicates that NGS can overcome this difficulty and allows building assemblies of complete organellar genomes even when no sequence reference is available [29], and with a low amount of DNA (300 ng). In addition, our approach is fast and relatively cheap already today, and consequently can be applied to generate complete mitogenomes on a large panel of nematode lineages. In turn, such an accumulation of data will permit an improved automated gene annotation in nematodes. This may open new avenues to identify informative regions for developing genetic tools for race identification or phylogeographic investigations [48], to study the rearrangements of mitogenomes in nematodes [28], and may also allow the development of phylogenomics to depict the long evolutionary history of this still poorly known animal phylum [4].

Disclosure of interest

The authors declare that they have no conflicts of interest concerning this article.

Acknowledgments

We are grateful to Alexander Donath and Matthias Bernt for helpful analysis of protein coding genes and general discussions on gene annotations, and to Malyna Suong for nematode culture and DNA extraction. We also gladly acknowledge the Plateforme Génomique and Bioinformatics platform of Genopole Toulouse Midi-Pyrénées for providing support to analyse NGS data and a referee for his comments. Professor Joong-Ki Park kindly provided the alignment of mitochondrial protein coding genes used in this study. This work was funded by the fellowship ANR-12-AGRI-0002 (Arimnet 2011–PestOlive). This work has been conducted in the “Laboratoire Évolution & Diversité biologique” (EDB), part of the LABEX entitled TULIP managed by “Agence nationale de la recherché” (ANR-10-LABX-0041). Authors also wish to thank a Global Rice Science Partnership (GRiSP) Grant (Menergep Project) for financial support to SB.


References

[1] D.Y.C. Wang; S. Kumar; B.S. Hedges Divergence time estimates for the early history of animal phyla and the origin of plants, animals and fungi, Proc. Roy. Soc. Lond. B, Volume 266 (1999), pp. 163-171

[2] M.L. Blaxter Nematoda: genes, genomes and the evolution of parasitism, Adv. Parasitol., Volume 54 (2003), pp. 101-195

[3] G.O. Poinar Nematoda and nematomorpha (J.H. Thorp; A.P. Covich, eds.), Ecology and Classification of North American Freshwater Invertebrates, Academic Press Inc, New York, 1991, pp. 249-283

[4] M.L. Blaxter; P. de Ley; J.R. Garey; L.X. Liu; P. Scheldeman; A. Vierstraete; J.R. Vanfleteren; L.Y. Mackey; M. Dorris; L.M. Frisse; J.T. Vida; W.K. Thomas A molecular evolutionary framework for the phylum Nematoda, Nature, Volume 392 (1998), pp. 71-75

[5] V.V. Malakhov Nematodes: Structure, Development, Classification, and Phylogeny (W.D. Hope, ed.), Smithsonian Institution Press, Washington/London, 1994

[6] J.T. Jones; A. Haegeman; E. Danchin; H.S. Gaur; J. Helder; M.G.K. Jones; T. Kikuchi; R. Manzanilla-Lopez; J.E. Palomares-Ruis; W.M.L. Wesemael; R.N. Perry Top 10 plant-parasitic nematodes in molecular plant pathology, Mol. Plant. Pathol., Volume 14 (2013), pp. 946-961

[7] W. Bert; F. Leliaert; A.R. Vierstraete; J.R. Vanfleteren; G. Borgonie Molecular phylogeny of the Tylenchina and evolution of the female gonoduct (Nematoda: Rhabditida), Mol. Phylogenet. Evol., Volume 48 (2008), pp. 728-744

[8] M. Holterman; G. Karssen; S. van den Elsen; H. van Megen; J. Bakker; J. Helder Small subunit rDNA-based phylogeny of the Tylenchida sheds light on relationships among some high-impact plant-parasitic nematodes and the evolution of plant feeding, Phytopathology, Volume 99 (2009), pp. 227-235

[9] J.D. Eisenback; H.H. Triantaphyllou Root-Knot Nematodes: Meloidogyne species and races (W.R. Nickle, ed.), Manual of Agricultural Nematology, Marcel Dekker, New York, 1991, pp. 281-286

[10] A.C. Triantaphyllou Cytogenetics, cytotaxonomy and phylogeny of root-knot nematodes (J.N. Sasser; C.C. Carter, eds.), An Advanced Treatise on Meloidogyne, Vol. 1, North Carolina State University Graphics, Raleigh, 1985, pp. 113-126

[11] P. Castagnone-Sereno; E. Danchin; L. Perfus-Babeoch; P. Abad Diversity and evolution of root-knot nematodes, genus Meloidogyne: new insights from the genomic era, Annu. Rev. Phytopathol., Volume 51 (2013), pp. 203-220

[12] S.B. Jepson, CAB International, Wallingford (1987), p. 265

[13] V.C. Blok; T.O. Powers Biochemical and molecular identification (R.N. Perry; M. Moens; J.L. Starr, eds.), Root-Knot Nematodes, CAB International, 2009, pp. 98-118

[14] J.N. Sasser; C.C. Carter, North Carolina University Graphics, Raleigh (1985), p. 422

[15] C. Zijlstra; D.T.H.M. Donkers-Venne; M. Fargette Identification of Meloidogyne incognita, M. javanica and M. arenaria using sequence characterised amplified region (SCAR) based PCR assays, Nematology, Volume 2 (2000), pp. 847-853

[16] A. Gilabert; J.D. Wasmuth Unravelling parasitic nematode natural history using population genetics, Trends Parasitol., Volume 29 (2013), pp. 438-448

[17] M.S. Tigano; R.M.D.G. Carneiro; A. Jeyaprakash; D.W. Dickson; B.J. Adams Phylogeny of Meloidogyne spp. based on 18S rDNA and the intergenic region of mitochondrial DNA sequences, Nematology, Volume 7 (2005), pp. 851-862

[18] D. Picard; O. Plantard What constitutes a population for the plant parasitic nematode Globodera pallida in its native area (Peru)?, Int. J. Parasitol., Volume 36 (2005), pp. 115-122

[19] L. Villate; D. Esmenjaud; M. Van Helden; S. Stoeckel; O. Plantard Genetic signature of amphimixis allows for the detection and fine scale localization of sexual reproduction events in a mainly parthenogenetic nematode, Mol. Ecol., Volume 19 (2010), pp. 856-873

[20] T.R. Gregory; J.A. Nicol; H. Tamm; B. Kullman; K. Kullman; I.J. Leitch; B.G. Murray; D.F. Kapraun; J. Greilhuber; M.D. Bennett Eukaryotic genome size databases, Nucleic Acids Res., Volume 35 (2007), pp. 332-338

[21] P. Abad; J. Gouzy; J.M.M. Aury et al. Genome sequence of the metazoan plant-parasitic nematode Meloidogyne incognita, Nat. Biotechnol., Volume 26 (2008), pp. 909-915

[22] T.J.C. Anderson; M.S. Blouin; R.N. Beech Population biology of parasitic nematodes: applications of genetic markers, Adv. Parasitol., Volume 41 (1998), pp. 219-283

[23] A. Hoolahan; V.C. Blok; T. Gibson; M. Dowton Paternal leakage of mitochondrial DNA in experimental crosses of populations of the potato cyst nematode Globodera pallida, Genetica, Volume 139 (2011), pp. 1509-1519

[24] C. Gutiérrez-Gutiérrez; P. Castillo; C. Cantalapiedra-Navarrete; B.B. Landa; S. Derycke; J.E. Palomares-Rius Genetic structure of Xiphinema pachtaicum and X. index populations based on mitochondrial DNA variation, Phytopathology, Volume 101 (2011), pp. 1168-1175

[25] R. Okimoto; H.M. Chamberlin; J.L. Macfarlane; D.R. Wolstenholme Repeated sequence sets in mitochondrial DNA molecules of root-knot nematodes (Meloidogyne): nucleotide sequences, genome location and potential for host-race identification, Nucleic Acids Res., Volume 19 (1991), pp. 1619-1626

[26] J.E. Jacob; B. Vanholme; T. van Leeuwen; G. Gheysen A unique genetic code change in the mitochondrial genome of the parasitic nematode Radopholus similis, BMC Res. Notes, Volume 2 (2009), p. 192

[27] T. Gibson; D. Farrugia; J. Barrett; D.J. Chitwood; J. Rowe; S. Subbotin; M. Dowton The mitochondrial genome of the soybean cyst nematode, Heterodera glycines, Genome, Volume 54 (2011), pp. 565-574

[28] T. Sultana; J. Kim; S.H. Lee; H. Han; S. Kim; G.S. Min; S.A. Nadler; J.K. Park Comparative analysis of complete mitochondrial genome sequences confirms independent origins of plant-parasitic nematodes, BMC Evol. Biol., Volume 13 (2013), p. 12

[29] C. Hahn; L. Bachmann; B. Chevreux Reconstructing mitochondrial genomes directly from genomic next-generation sequencing reads – a baiting and iterative mapping approach, Nucleic Acids Res., Volume 41 (2013), p. e129

[30] Inc. Illumina HiSeq™ 2000 sequencing system, 2010 ( http://supportres.illumina.com/documents/documentation/system_documentation/hiseq2000/hiseq2000_ug_15011190_r.pdf )

[31] I.R. Soriano; V. Schmit; D.S. Brar; J.C. Prot; G. Reversat Resistance to rice root-knot nematode Meloidogyne graminicola identified in Oryza longistaminata and O. glaberrima, Nematology, Volume 1 (1999), pp. 395-398

[32] S. Bellafiore; Z.X. Shen; M.N. Rosso; P. Abad; P. Shih; S.P. Briggs Direct identification of the Meloidogyne incognita secretome reveals proteins with host cell reprogramming potential, PLoS Pathog., Volume 4 (2008), p. e1000192

[33] J. Sambrook; E.F. Fritsch; T. Maniatis Molecular Cloning: a Laboratory Manual, Cold Spring Harbor Laboratory Press, Cold Spring Harbor, NY, USA, 1989

[34] G. Besnard; P.A. Christin; P.J.G. Malé; E. Coissac; H. Ralimanana; M. Vorontsova Phylogenomics and taxonomy of Lecomtelleae (Poaceae), an isolated panicoid lineage from Madagascar, Ann. Bot., Volume 112 (2013), pp. 1057-1066

[35] D.R. Zerbino; E. Birney Velvet: algorithms for de novo short read assembly using de Bruijn graphs, Genome Res., Volume 18 (2008), pp. 821-829

[36] M. Bernt; A. Donath; F. Jühling; F. Externbrink; C. Florentz; G. Fritzsch; J. Putz; M. Middendorf; P.F. Stadler Mitos: improved de novo metazoan mitochondrial genome annotation, Mol. Phylogenet. Evol., Volume 69 (2013), pp. 313-319

[37] P. Rice; I. Longden; A. Bleasby EMBOSS: the European Molecular Biology Open Software Suite, Trends Genet., Volume 16 (2000), pp. 276-277

[38] K. Katoh; D.M. Standley Mafft multiple sequence alignment software version 7: improvements in performance and usability, Mol. Biol. Evol., Volume 30 (2013), pp. 772-780

[39] F. Jühling; J. Pütz; C. Florentz; P.F. Stadler Armless mitochondrial tRNAs in Enoplea (Nematoda), RNA Biol., Volume 9 (2012), pp. 1161-1166

[40] E.P. Nawrocki; D.L. Kolbe; S.R. Eddy Infernal 1.0: inference of RNA alignments, Bioinformatics, Volume 25 (2009), pp. 1335-1337

[41] F. Jühling; J. Pütz; M. Bernt; A. Donath; M. Middendorf; C. Florentz; P.F. Stadler Improved systematic tRNA gene annotation allows new insights into the evolution of mitochondrial tRNA structures and into the mechanisms of mitochondrial genome rearrangements, Nucleic Acids Res., Volume 40 (2012), pp. 2833-2845

[42] S. Wende; E.G. Platzer; F. Jühling; J. Pütz; C. Florentz; P.F. Stadler; M. Mörl Biological evidence for the world's smallest tRNAs, Biochimie, Volume 100 (2014), pp. 151-158

[43] S. Guindon; J....F. Dufayard; V. Lefort; M. Anisimova; W. Hordijk; O. Gascuel New algorithms and methods to estimate maximum-likelihood phylogenies: assessing the performance of PhyML 3.0, Syst. Biol., Volume 59 (2010), pp. 307-321

[44] A. Hugall; J.... Stanton; C. Moritz Evolution of the AT-rich mitochondrial DNA of the root-knot nematode, Meloidogyne hapla, Mol. Biol. Evol., Volume 14 (1997), pp. 40-48

[45] M. Bernt; A. Braband; B. Schierwater; P.F. Stadler Genetic aspects of mitochondrial genome evolution, Mol. Phylogenet. Evol., Volume 69 (2013), pp. 328-338

[46] D.R. Wolstenholme; J....L. Macfarlane; R. Okimoto; D.O. Clary; J....A. Wahleithner Bizarre tRNAs inferred from DNA sequences of mitochondrial genomes of nematode worms, Proc. Natl. Acad. Sci. U S A, Volume 84 (1987), pp. 1324-1328

[47] P.B. Klimov; B.M. Oconnor Improved tRNA prediction in the American house dust mite reveals widespread occurrence of extremely short minimal tRNAs in acariform mites, BMC Genomics, Volume 10 (2009), p. 598

[48] D.A. Humphreys-Pereira; A.A. Elling Intraspecific variability and genetic structure in Meloidogyne chitwoodi from the USA, Nematology, Volume 15 (2013), pp. 315-327


Comments - Policy