Skip to main content

The mitochondrial genome of Protostrongylus rufescens – implications for population and systematic studies



Protostrongylus rufescens is a metastrongyloid nematode of small ruminants, such as sheep and goats, causing protostrongylosis. In spite of its importance, the ecology and epidemiology of this parasite are not entirely understood. In addition, genetic data are scant for P. rufescens and related metastrongyloids.


The mt genome was amplified from a single adult worm of P. rufescens (from sheep) by long-PCR, sequenced using 454-technology and annotated using bioinformatic tools. Amino acid sequences inferred from individual genes of the mt genomes were concatenated and subjected to phylogenetic analysis using Bayesian inference.


The circular mitochondrial genome was 13,619 bp in length and contained two ribosomal RNA, 12 protein-coding and 22 transfer RNA genes, consistent with nematodes of the order Strongylida for which mt genomes have been determined. Phylogenetic analysis of the concatenated amino acid sequence data for the 12 mt proteins showed that P. rufescens was closely related to Aelurostrongylus abstrusus, Angiostrongylus vasorum, Angiostrongylus cantonensis and Angiostrongylus costaricensis.


The mt genome determined herein provides a source of markers for future investigations of P. rufescens. Molecular tools, employing such mt markers, are likely to find applicability in studies of the population biology of this parasite and the systematics of lungworms.


Protostrongylus rufescens is a metastrongyloid nematode of small ruminants, including sheep and goats (definitive hosts) in most parts of the world [1]. The dioecious adults of this nematode live in the respiratory system (terminal bronchioles and alveoli) of the definitive host. Here, the females produce eggs, from which first-stage larvae (L1s) hatch within the airways of the lung. L1s then migrate via the bronchial/tracheal escalator to the pharynx, are swallowed and are then excreted in the faeces. L1s infect a molluscan intermediate host (snail) and then develop, under favourable environmental conditions, into third-stage larvae (L3) [1]. L3s within an infected intermediate host are then ingested by the ruminant host, penetrate the gut wall and then migrate via the lymphatic system or blood stream to the lungs, where they develop to adult worms. The prepatent period is reported to be ~ 4–9 weeks [2]. Although P. rufescens infection is widespread, it does not usually cause major clinical disease. Nonetheless, pathological changes, characterized by chronic, eosinophilic, granulomatous pneumonia, can be detected upon post mortem examination. Adult worms reside mainly in the bronchioles and alveoli, and are surrounded by macrophages, giant cells, eosinophils and other inflammatory cells which produce grey or beige plaques (1–2 cm) under the pleura in the dorsal border of the caudal lung lobes [3].

Little is known about fundamental aspects of the epidemiology and ecology of P. rufescens. Molecular tools employing suitable genetic markers can underpin fundamental studies in these areas, with a perspective on investigating transmission patterns linked to particular genotypes of a parasite and on discovering population variants or cryptic species [4, 5]. Advances in nucleic acid sequencing and bioinformatics have provided a foundation for characterizing the mt genomes from parasitic nematodes as a source of genetic markers for such explorations. Here, we used an established, massively parallel sequencing-bioinformatics pipeline [6] for the characterization of the mt genome of P. rufescens, which we compared with those of related metastrongyloid nematodes, for which mt genomic sequence data are available. We also studied the genetic relationships among these lungworms and selected representatives within the order Strongylida, and suggest that selected regions in the genome of P. rufescens should serve well as markers for future studies of the ecology and epidemiology of this nematode around the world.


Parasite and genomic DNA isolation

Adult worms of P. rufescens were collected from the lungs of a fresh sheep cadaver in Victoria, Australia, washed extensively in physiological saline and then stored at −80°C. Upon thawing, genomic DNA was isolated from a single adult male specimen using an established method of sodium dodecyl-sulphate (SDS)/proteinase K digestion and subsequent mini-column purification [7]. The identity of the specimen was verified by PCR-based sequencing (BigDye chemistry v.3.1) of the second internal transcribed spacer (ITS-2) of nuclear ribosomal DNA [7].

Long-PCR, sequencing and mt genome assembly

From the genomic DNA extracted from the single male worm, the complete mt genome was amplified by long-PCR (BD Advantage 2, BD Biosciences) as two overlapping amplicons (~5 kb and ~10 kb), using the protocol described by Hu et al. [8], with appropriate positive (i.e., Haemonchus contortus DNA) and negative (i.e., no template) controls. Amplicons were consistently produced from the positive control samples; in no case was a product detected for the negative controls. Amplicons were then treated with shrimp alkaline phosphatase and exonuclease I [9], and quantified by spectrophotometry. Following agarose electrophoretic analysis, the two amplicons (2.5 μg of each) were pooled and subsequently sequenced using the 454 Genome Sequencer FLX (Roche) [10] according to an established protocol [6]. The mt genome sequence was assembled using the program CAP3 [11] from individual reads (of ~300 bp).

Annotation and analyses of sequence data

Following assembly, the mt genome of P. rufescens was annotated using the bioinformatic annotation pipeline developed by Jex et al. [6]. Briefly, the open reading frame (ORF) of each protein-coding mt gene was identified (six reading frames) by comparison to those of the mt genome of Angiostrongylus vasorum [GenBank: JX268542; [12]]. The large and small subunits of the mt ribosomal RNA genes (rrn S and rrn L, respectively) were identified by local alignment. The transfer RNA (tRNA) genes were predicted (from both strands) based on their structure, using scalable models based on the standard mt tRNAs for nematodes [5]. Predicted tRNA genes were then grouped according to their anti-codon sequence and identified based on the amino acid encoded by the anti-codon. Two separate tRNA gene groups were predicted each for leucine (Leu) (one each for the anticodons CUN and UUR, respectively) and for serine (Ser) (one each for the anticodons AGN and UCN, respectively), as these tRNA genes are duplicated in many invertebrate mt genomes, including those of nematodes [5]. All predicted tRNAs for each amino acid group were ranked according to the “strength” of their structure (inferred based on minimum nucleotide mismatches in each stem); for each group, the 100 best-scoring structures were compared by BLASTn against a database comprising all tRNA genes for each amino acid for all published mt genome sequences of nematodes (available via; [13]). The tRNA genes were then identified and annotated based on their highest sequence identity to known nematode tRNAs. Annotated sequence data were imported using the program SEQUIN (via, the mt genome structure verified and the final sequence submitted as an SQN file to the GenBank database.

Phylogenetic analysis of concatenated amino acid sequence datasets

The amino acid sequences were predicted from individual mt genes of P. rufescens and of other nematodes, including An. cantonensis, An. costaricensis, An. vasorum, Metastrongylus pudendotectus and M. salmi [GenBank: GQ398122, JX268542, GQ398121, GQ888714 and GQ888715, respectively; Metastrongyloidea]; Ancylostoma caninum and Necator americanus [GenBank: FJ483518 and NC_003416, respectively; Anyclostomatoidea]; H. contortus and Trichostrongylus axei [GenBank: NC_010383 and GQ888719, respectively; Trichostrongyloidea]; Oesophagostomum dentatum and Strongylus vulgaris [GenBank: GQ888716 and GQ888717, respectively; Strongyloidea]; and Strongyloides stercoralis [GenBank: AJ558163; Strongyloidoidea] [6, 12, 1420] (Table 1). All amino acid sequences were aligned using the program MUSCLE [21] and then subjected to phylogenetic analysis. For this analysis, best-fit models of evolution were selected using ProtTest 3.0 [22] employing the Akaike information criterion (AIC) [23]. Bayesian inference analysis was conducted using MrBayes 3.1.2 [24], with a fixed mtREV amino acid substitution model [25], using four rate categories approximating a Γ distribution, four chains and 200,000 generations, sampling every 100th generation. The first 200 generations were removed from the analysis as burn-in.

Table 1 Details of the whole mitochondrial genome sequences used in this study as reference sequences

Results and discussion

Features of the mt genome

The circular mt genome sequence of P. rufescens [GenBank: KF481953] is 13,619 bp in length (Figure 1). It contains two ribosomal genes, 12 protein-coding (cox 1-3, nad 1-6, nad 4L, atp 6 and cyt b) and 22 tRNA genes. The gene arrangement (GA2) in the mt genome of P. rufescens was the same as all other strongylid nematodes studied to date [5, 26]. All of the 36 genes are transcribed in the same direction (5′ to 3′) (Figure 1). Overall, the genome is AT-rich, as expected for strongylid nematodes [12, 20, 27, 28], with T being the most favoured nucleotide and C the least favoured. The nucleotide contents were 25.9% (A), 6.8% (C), 18.6% (G) and 48.6% (T) (Table 2). The longest non-coding (AT-rich) region, located between the genes trn A and trn P, was 223 bp in length (see Figure 1); its AT-content was 83.4%, significantly greater than for all other parts of the mt genome (Table 2).

Figure 1

Schematic representation of the circular mt genome of Protostrongylus rufescens. Each transfer RNA gene is identified by a one-letter amino acid code in the map (external), and the AT rich region is also indicated. All genes are transcribed in the clockwise direction.

Table 2 Nucleotide composition (%) for the entire or regions of the mitochondrial genome of Protostrongylus rufescens

Ribosomal RNA genes

The rrn S and rrn L genes of P. rufescens were identified by sequence comparison with An. vasorum. The rrn S gene was located between trn E and trn S (UCN), and rrn L was between trn H and nad 3. The two genes were separated from one another by the protein-encoding genes nad 3, nad 5, nad 6 and nad 4L (Figure 1). The sizes of the rrn S and rrn L genes of P. rufescens were 683 bp and 959 bp, respectively. The lengths of these two genes were similar to those of other metastrongyloids for which mt genomes are known (694–699 bp for rrn S, and 958–961 bp for rrn L [12, 20, 2628] (Figure 1), and amongst the shortest for metazoan organisms [29].

Protein-coding genes and codon usage

The prediction of initiation and termination codons for the protein-coding genes of P. rufescens (Table 3) revealed that the commonest start codon was ATT (for five of 12 proteins), followed by TTG (four genes), ATA (two genes) and ATG (one gene). Ten mt protein genes of P. rufescens were predicted to have a TAA or TAG translation termination codon. The other two protein genes ended in an abbreviated stop codon, such as T or TA (Table 3).

Table 3 Summary of the mitochondrial genome of Protostrongylus rufescens

The codon usage for the 12 protein-encoding genes of P. rufescens was also compared with that of other metastrongyloid nematodes, Aelurostrongylus (Ae.) abstrusus, An. cantonensis, An. costaricensis and An. vasorum[12, 20, 28] (Table 4). All 64 codons were used. The preferred nucleotide usage at the third codon position of mt protein genes of P. rufescens reflects the overall nucleotide composition of the mt genome. At this position, T was the most frequently, and C the least frequently used. For P. rufescens, the codons ending in A had higher frequencies than the codons ending in G, which is similar to, for example, other members of the order Strongylida and Caenorhabditis elegans (Rhabditida), but distinct from Ascaris suum (Ascaridida) and Onchocerca volvulus (Spirurida) [1417, 30]. As the usage of synonymous codons is proposed to be preferred in gene regions of functional importance, codon bias appears to be linked to selection at silent sites and to translation efficiency [31, 32].

Table 4 Codon usages (%) in mitochondrial protein-encoding genes of Protostrongylus rufescens

The AT bias in the genome is also reflected in the amino acid composition of predicted proteins. The AT-rich codons represent the amino acids Phe, Ile, Met, Tyr, Asn the Lys, and GC-rich codons represent Pro, Ala, Arg the Gly. In the mt genome of P. rufescens, the most frequently used codons were TTT (Phe), TTA (Leu), ATT (Ile), TTG (Leu), TAT (Tyr), GGT (Gly), AAT (Asn) and GTT (Val). Six of these codons are AT-rich, and one of them is GC-rich. Seven of the eight codons contained an A or a T at two positions, except for GGT (Gly), which contained a T only in the third position. None of them had a C at any position. The least frequently used codons were CTC, CTG (Leu), GTC (Val), AGC (Ser), CCC (Pro), GCC (Ala), CAC (His), CGA (Arg), TCC (Ser), GGC (Gly) and ACC (Thr). All four GC-rich codons were represented here, and every codon had at least one C. When the frequencies of synonymous codons within the AT-rich group, such as Phe (TTT, 14.2%; TTC, 1.2%), Ile (ATT, 5.6%; ATC, 0.7%), Tyr (TAT, 5.6%; TAC, 0.9%) and Asn (AAT, 3.8%; AAC, 0.7%), were compared, the frequency was always less if the third position was a C.

Transfer RNA genes

Twenty-two tRNA gene sequences were predicted in the mt genome of P. rufescens. These sequences ranged from 52–63 nt in length. The tRNA structures had a 7 bp amino-acyl arm, a 4 bp DHU arm, a 5 bp anticodon stem, a 7 base anticodon loop, a T always preceding an anticodon as well as a purine always following an anticodon. Twenty of the 22 tRNA genes (i.e. excluding the two trn S genes) have a predicted secondary structure with a 4 bp DHU stem and a DHU loop of 4–10 bases, in which the variable TψC arm and loop are replaced by a “TV-replacement loop” of 4–11 bases, in accordance with most nematodes whose mt genomes have been characterised [5]. The mt trn S for P. rufescens has a secondary structure consisting of a DHU replacement loop of 7 bases, 3 bp TψC arm, TψC loop of 4–6 bases and a variable loop of 3 bases, consistent with other members of the Chromadorea [6, 14, 20, 33], but different from the enoplid nematodes Trichinella murrelli and T. spiralis[29, 34, 35]. Overlaps of one to four nucleotides are found between the genes trn H and rrn L, nad 4L and trn W, trn Y and nad 1, trn I and trn R within the mt genome of P. rufescens.

Amino acid sequence comparisons and genetic relationships of P. rufescens with metastrongyloid and other nematodes

The amino acid sequences predicted from individual protein-encoding mt genes of P. rufescens were compared with those of Ae. abstrusus, An. cantonensis, An. costaricensis, An. vasorum, Dictyocaulus viviparus and D. eckerti (Table 5). Pairwise comparisons of the concatenated sequences revealed identities of 37.0-92.4% between these species. Based on identity, COX1 was the most conserved protein, whereas NAD2 and NAD6 were the least conserved. Phylogenetic analysis of the concatenated amino acid sequence data for the 12 mt proteins showed that P. rufescens was more closely related to Ae. abstrusus, An. cantonensis, An. costaricensis and An. vasorum, (pp = 1.00) than to M. pudendotectus and M. salmi (metastrongloids) (pp = 1.00), to the exclusion of H. contortus, T. axei (trichostrongyloids), Anc. caninum, N. americanus (hookworms; ancylostomatoids) and O. dentatum and S. vulgaris (strongyloids) (Figure 2) (pp = 1.00).

Table 5 Pairwise comparison of the amino acid sequences of the 12 protein-encoding mitochondrial genes
Figure 2

Phylogenetic relationship of Protostrongylus rufescens with other nematodes. Concatenated amino acid sequence data for all protein-encoding mitochondrial genes of Protostrongylus rufescens (bold) and other metastrongyloids, including Aelurostrongylus abstrusus, Angiostrongylus cantonensis, An. costaricensis, An. vasorum, Metastrongylus pudendotectus and M. salmi (metastrongyloids), as well as other concatenated sequence data representing different superfamilies, including Ancylostoma caninum and Necator americanus (hookworms; ancylostomatoids); Haemonchus contortus and Trichostrongylus axei (trichostrongyloids); Oesophagostomum dentatum and Strongylus vulgaris (strongyloids); and Strongyloides stercoralis (a rhabditid outgroup) were analyzed using Bayesian inference. The numbers above each tree branch represent the statistical support for each node (based on posterior probability [pp] score). GenBank accession numbers are in round brackets.


The characterisation of the mt genome of P. rufescens provides genetic markers for future population genetic and systematic studies. As sequence variation in ITS-2 nuclear rDNA is usually low within most species of strongylid nematodes [36], mt DNA is better suited for assessing population genetic variation. Therefore, PCR-based analytical approaches, using cox 1, nad 1 and nad 4 (displaying varying levels of within-species divergence), could be used to study haplotypic variation in P. rufescens populations in sheep and goats and also in molluscan hosts. Given that species complexes are commonly encountered in bursate nematodes [1, 4, 36], it would be interesting to prospect for cryptic species, to assess whether distinct genotypes/haplotypes of P. rufescens exist in sheep and goats as well as snails [37], and to establish whether particular sub-populations of P. rufescens occur in particular environments or geographical regions/countries, and have particular patterns of transmission.

It would also be interesting to assess the genetic structure of P. rufescens populations using PCR-coupled mutation scanning and sequencing of selected mt gene regions (such as cox 1 and nad 4), and mt DNA diversity within populations and the gene flow among populations. Findings for this lungworm (with an indirect life cycle via a molluscan intermediate host) could be compared with those for D. viviparus (with a direct life cycle), which has been reported to have surprisingly low mt DNA diversity within populations and limited gene flow among populations [38, 39].

The complete mt genome of P. rufescens provides a basis for extended comparative mt genomic/proteomic analyses of other protostrongyloids of ruminants, including P. brevispiculum, P. davtiani, P. hobmaieri, P. rushi, P. skrjabini, P. stilesi, Cystocaulus ocreatus, Neostrongylus lineatus, Muellerius capillaris (the latter of which is a particularly pathogenic parasite in goats), and those of other animal hosts, such as lagomorphs and pinnipeds. Given the utility of predicted mt proteomic datasets, high phylogenetic signal and consistently high nodal support values in recent systematic analyses [6, 12, 27, 28, 33] provide an opportunity to reassess the evolutionary relationships of lungworms (order Strongylida). For example, the family Protostrongylidae is distinguished from other metastrongyloids by only a couple of morphological characters, i.e., the gubernaculm and telamon in adult male worms [40], and it is proposed that protostrongyloids of lagomorphs originated from their ancestors primarily infecting sheep, goat, antelopes and deer [41]. Analyses of inferred mt proteomic data sets from a range of protostrongyloids should allow relationships within the family Protostrongylidae and also the origin of the protostrongylids of lagomorphs to be assessed. In addition, there has been considerable debate as to the relationships among suborders within the Strongylida, based on the use of phenotypic characters [42]. On one hand, it has been hypothesized that the suborder Metastrongylina (to which species of Protostrongylus, Metastrongylus, Aelurostrongylus and Angiostrongylus belong) originated from ancestors in the Strongylina [43, 44] or Trichostrongylina [45, 46]. On the other hand, it has been proposed that the Metastrongylina gave rise to the Strongylina [47]. To date, molecular phylogenetic analyses of nuclear ribosomal rDNA sequence data [48, 49] have suggested that the Trichostrongylina are basal to the Metastrongylina, which represented a monophyletic assemblage. However, Jex et al. [6], using mitochondrial sequence data, showed that the major suborders within the Strongylida (e.g., the Metastrongylina, Strongylina and Trichostrongylina) were each resolved as distinct, monophyletic clades with maximum statistical and nodal support (posterior probability = 1.00; bootstrap = 100). A detailed analysis using inferred mt proteomic data sets would allow an independent assessment of the systematic relationships of these suborders.


Comparative analyses of proteomic sequence datasets inferred from the mt genomes of P. rufescens and other lungworms indicate that P. rufescens is closely related to Ae. abstrusus, An. cantonensis, An. costaricensis and An. vasorum. The mt genome determined herein should provide a source of markers for future investigations of P. rufescens. Molecular tools, employing such mt markers, are likely to find applicability in studies of the population biology of this parasite and the systematics of lungworms.

Authors’ information

Abdul Jabbar and Namitha Mohandas shared first authorship.



Adenosine triphosphatase subunit 6


Cytochrome c subunit


Cytochrome b subunit


Second internal transcribed spacer




Nicotinamide dehydrogenase subunit


Open reading frame


Nuclear ribosomal DNA


Ribosomal large subunit


Ribosomal small subunit


Transfer RNA.


  1. 1.

    Anderson RC: Nematode Parasites of Vertebrates. Their Development and Transmission. 2000, Wallingford: CAB International, 167-

    Google Scholar 

  2. 2.

    Forrester DJ, Senger CM: Prenatal infection of bighorn sheep with protostrongylid lungworms. Nature. 1964, 201: 1051-

    CAS  Article  PubMed  Google Scholar 

  3. 3.

    Mansfield LS, Gamble HR: Alveolar mastocytosis and eosinophilia in lambs with naturally acquired nematode infections of Protostrongylus rufescens and Haemonchus contortus. Vet Immunol Immunopathol. 1995, 49: 251-262. 10.1016/0165-2427(95)05466-9.

    CAS  Article  PubMed  Google Scholar 

  4. 4.

    Blouin MS: Molecular prospecting for cryptic species of nematodes: mitochondrial DNA versus internal transcribed spacer. Int J Parasitol. 2002, 32: 527-531. 10.1016/S0020-7519(01)00357-5.

    CAS  Article  PubMed  Google Scholar 

  5. 5.

    Hu M, Gasser RB: Mitochondrial genomes of parasitic nematodes - progress and perspectives. Trends Parasitol. 2006, 22: 78-84. 10.1016/

    CAS  Article  PubMed  Google Scholar 

  6. 6.

    Jex AR, Hall RS, Littlewood DT, Gasser RB: An integrated pipeline for next-generation sequencing and annotation of mitochondrial genomes. Nucleic Acids Res. 2010, 38: 522-533. 10.1093/nar/gkp883.

    PubMed Central  CAS  Article  PubMed  Google Scholar 

  7. 7.

    Gasser RB, Hu M, Chilton NB, Campbell BE, Jex AR, Otranto D, Cafarchia C, Beveridge I, Zhu XQ: Single-strand conformation polymorphism (SSCP) for the analysis of genetic variation. Nat Protoc. 2006, 1: 3121-3128.

    CAS  Article  PubMed  Google Scholar 

  8. 8.

    Hu M, Jex AR, Campbell BE, Gasser RB: Long PCR amplification of the entire mitochondrial genome from individual helminths for direct sequencing. Nat Protoc. 2007, 2: 2339-2344. 10.1038/nprot.2007.358.

    CAS  Article  PubMed  Google Scholar 

  9. 9.

    Werle E, Schneider C, Renner M, Völker M, Fiehn W: Convenient single-step, one tube purification of PCR products for direct sequencing. Nucleic Acids Res. 1994, 22: 4354-4355. 10.1093/nar/22.20.4354.

    PubMed Central  CAS  Article  PubMed  Google Scholar 

  10. 10.

    Margulies M, Egholm M, Altman WE, Attiya S, Bader JS, Bemben LA, Berka J, Braverman MS, Chen YJ, Chen Z, Dewell SB, Du L, Fierro JM, Gomes XV, Godwin BC, He W, Helgesen S, Ho CH, Irzyk GP, Jando SC, Alenquer ML, Jarvie TP, Jirage K, Kim JB, Knight JR, Lanza JR, Leamon JH, Lefkowitz SM, Lei M, Li J, Lohman KL: Genome sequencing in microfabricated high-density picolitre reactors. Nature. 2005, 437: 376-380.

    PubMed Central  CAS  PubMed  Google Scholar 

  11. 11.

    Huang X, Madan A: CAP3: A DNA sequence assembly program. Genome Res. 1999, 9: 868-877. 10.1101/gr.9.9.868.

    PubMed Central  CAS  Article  PubMed  Google Scholar 

  12. 12.

    Gasser RB, Jabbar A, Mohandas N, Schnyder M, Deplazes P, Littlewood DTJ, Jex AR: Mitochondrial genome of Angiostrongylus vasorum: comparison with congeners and implications for studying the population genetics and epidemiology of this parasite. Infect Genet Evol. 2012, 12: 1884-1891. 10.1016/j.meegid.2012.07.022.

    CAS  Article  PubMed  Google Scholar 

  13. 13.

    Jameson D, Gibson AP, Hudelot C, Higgs PG: OGRe: a relational database for comparative analysis of mitochondrial genomes. Nucleic Acids Res. 2003, 31: 202-206. 10.1093/nar/gkg077.

    PubMed Central  CAS  Article  PubMed  Google Scholar 

  14. 14.

    Okimoto R, Macfarlane L, Clary D, Wolstenholme DR: The mitochondrial genomes of two nematodes, Caenorhabditis elegans and Ascaris suum. Genetics. 1992, 130: 471-498.

    PubMed Central  CAS  PubMed  Google Scholar 

  15. 15.

    Hu M, Chilton NB, Gasser RB: The mitochondrial genomes of the human hookworms, Ancylostoma duodenale and Necator americanus (Nematoda: Secernentea). Int J Parasitol. 2002, 32: 145-158. 10.1016/S0020-7519(01)00316-2.

    CAS  Article  PubMed  Google Scholar 

  16. 16.

    Hu M, Chilton NB, Gasser RB: The mitochondrial genome of Strongyloides stercoralis (Nematoda) – idiosyncratic gene order and evolutionary implications. Int J Parasitol. 2003, 33: 1393-1408. 10.1016/S0020-7519(03)00130-9.

    CAS  Article  PubMed  Google Scholar 

  17. 17.

    Li MW, Lin RQ, Song HQ, Wu XY, Zhu XQ: The complete mitochondrial genomes for three Toxocara species of human and animal health significance. BMC Genomics. 2008, 16 (9): 224-

    Article  Google Scholar 

  18. 18.

    Jex AR, Hu M, Littlewood DT, Waeschenbach A, Gasser RB: Using 454 technology for long-PCR based sequencing of the complete mitochondrial genome from single Haemonchus contortus (Nematoda). BMC Genomics. 2008, 11: 9-11.

    Google Scholar 

  19. 19.

    Jex AR, Waeschenbach A, Hu M, van Wyk JA, Beveridge I, Littlewood DT, Gasser RB: The mitochondrial genomes of Ancylostoma caninum and Bunostomum phlebotomum - two hookworms of animal health and zoonotic importance. BMC Genomics. 2009, 11: 10-79.

    Google Scholar 

  20. 20.

    Lv S, Zhang Y, Zhang L, Liu Q, Liu HX, Hu L, Wei FR, Steinmann P, Graeff-Teixeira C, Zhou XN, Utzinger J: The complete mitochondrial genome of the rodent intra-arterial nematodes Angiostrongylus cantonensis and Angiostrongylus costaricensis. Parasitol Res. 2012, 111: 115-123. 10.1007/s00436-011-2807-4.

    Article  PubMed  Google Scholar 

  21. 21.

    Edgar RC: MUSCLE: multiple sequence alignment with high accuracy and high throughput. Nucleic Acids Res. 2004, 32: 1792-1797. 10.1093/nar/gkh340.

    PubMed Central  CAS  Article  PubMed  Google Scholar 

  22. 22.

    Darriba D, Taboada GL, Doallo R, Posada D: ProtTest 3: fast selection of best-fit models of protein evolution. Bioinformatics. 2011, 27: 1164-1165. 10.1093/bioinformatics/btr088.

    CAS  Article  PubMed  Google Scholar 

  23. 23.

    Akaike H: A new look at the statistical model identification. IEEE T Automat Contr. 1974, 19: 16-723.

    Google Scholar 

  24. 24.

    Huelsenbeck JP, Ronquist F: MRBAYES: Bayesian inference of phylogeny. Bioinformatics. 2001, 17: 754-755. 10.1093/bioinformatics/17.8.754.

    CAS  Article  PubMed  Google Scholar 

  25. 25.

    Adachi J, Hasegawa M: Model of amino acid substitution in proteins encoded by mitochondrial DNA. J Mol Evol. 1996, 42: 459-468. 10.1007/BF02498640.

    CAS  Article  PubMed  Google Scholar 

  26. 26.

    Jex AR, Littlewood DT, Gasser RB: Towards next-generation sequencing of mitochondrial genomes - focus on parasitic worms of animals and biotechnological implications. Biotechnol Adv. 2010, 28: 151-159. 10.1016/j.biotechadv.2009.11.002.

    CAS  Article  PubMed  Google Scholar 

  27. 27.

    Gasser RB, Jabbar A, Mohandas N, Höglund J, Hall RS, Littlewood DTJ, Jex AR: Assessment of the genetic relationship between Dictyocaulus species from Bos taurus and Cervus elaphus using complete mitochondrial genomic datasets. Parasit Vectors. 2013, 5: 241-

    Article  Google Scholar 

  28. 28.

    Jabbar A, Jex AR, Mohandas N, Hall RS, Littlewood DTJ, Gasser RB: The mitochondrial genome of Aelurostrongylus abstrusus – diagnostic, epidemiological and systematic implications. Gene. 2013, 516: 294-300. 10.1016/j.gene.2012.10.072.

    CAS  Article  PubMed  Google Scholar 

  29. 29.

    Lavrov DV, Brown WM: Trichinella spiralis mtDNA: a nematode mitochondrial genome that encodes a putative ATP8 and normally structured tRNAs and has a gene arrangement relatable to those of coelomate metazoans. Genetics. 2001, 157: 621-637.

    PubMed Central  CAS  PubMed  Google Scholar 

  30. 30.

    Keddie EM, Higazi T, Unnasch TR: The mitochondrial genome of Onchocerca volvulus: sequence, structure and phylogenetic analysis. Mol Biochem Parasitol. 1998, 95: 111-127. 10.1016/S0166-6851(98)00102-9.

    CAS  Article  PubMed  Google Scholar 

  31. 31.

    Sharp PM, Matassi G: Codon usage and genome evolution. Curr Opin Genet Dev. 1994, 4: 851-860. 10.1016/0959-437X(94)90070-1.

    CAS  Article  PubMed  Google Scholar 

  32. 32.

    Durent L, Mouchiroud D: Expression pattern and, surprisingly, gene length shape codon usage in Caenorhabditis, Drosophila and Arabidopsis. Proc Natl Acad Sci U S A. 1999, 96: 4482-4487. 10.1073/pnas.96.8.4482.

    Article  Google Scholar 

  33. 33.

    Park JK, Sultana T, Lee SH, Kang S, Kim HK, Min GS, Eom KS, Nadler SA: Monophyly of clade III nematodes is not supported by phylogenetic analysis of complete mitochondrial genome sequences. BMC Genomics. 2011, 12: 392-10.1186/1471-2164-12-392.

    PubMed Central  CAS  Article  PubMed  Google Scholar 

  34. 34.

    Webb KM, Rosenthal BM: Deep resequencing of Trichinella spiralis reveals previously un-described single nucleotide polymorphisms and intra-isolate variation within the mitochondrial genome. Infect Genet Evol. 2010, 10: 304-310. 10.1016/j.meegid.2010.01.003.

    CAS  Article  PubMed  Google Scholar 

  35. 35.

    Webb KM, Rosenthal BM: Next-generation sequencing of the Trichinella murrelli mitochondrial genome allows comprehensive comparison of its divergence from the principal agent of human trichinellosis, Trichinella spiralis. Infect Genet Evol. 2011, 11: 116-123. 10.1016/j.meegid.2010.10.001.

    CAS  Article  PubMed  Google Scholar 

  36. 36.

    Gasser RB: Molecular tools – advances, opportunities and prospects. Vet Parasitol. 2006, 136: 69-89. 10.1016/j.vetpar.2005.12.002.

    CAS  Article  PubMed  Google Scholar 

  37. 37.

    Bowman DD: Respiratory system parasites of the dog and cat (Part II): trachea and bronchi, and pulmonary vessels. Companion and Exotic Animal Parasitology. Edited by: Bowman DD. 2000, Ithaca NY, USA: International Veterinary Information Service

    Google Scholar 

  38. 38.

    Hu M, Chilton NB, Gasser RB: The mitochondrial genomics of parasitic nematodes of socio-economic importance: recent progress, and implications for population genetics and systematics. Adv Parasitol. 2004, 56: 133-212.

    Article  PubMed  Google Scholar 

  39. 39.

    Höglund J, Morrison DA, Mattsson JG, Engström A: Population genetics of the bovine/cattle lungworm (Dictyocaulus viviparus) based on mtDNA and AFLP marker techniques. Parasitology. 2006, 133: 89-99. 10.1017/S0031182006009991.

    Article  PubMed  Google Scholar 

  40. 40.

    Boev SN: Protostrongylids. Fundamentals of Nematology. 1975, USSR: Academy of Science of the USSR, Translated by Amerind Publishing Co., Pvt. Ltd. New Delhi,1984, 25

    Google Scholar 

  41. 41.

    Anderson TJC, Blouin MS, Beech RN: Populations biology of parasitic nematodes: applications of genetic markers. Adv Parasitol. 1998, 41: 219-283.

    CAS  Article  PubMed  Google Scholar 

  42. 42.

    Skrjabin KI, Shikhobalova NP, Schulz RS, Popova TI, Boev SN, Delyamure SL: Key to Parasitic Nematodes. III. Strongylata. 1952, Moscow: Akad Nauk SSR, [English translation by Israel Program for Scientific Translations, 1961]

    Google Scholar 

  43. 43.

    Dougherty EC: The phylogeny of the nematode family Metastrongylidae Leiper, (1909): a correlation of host and symbiote evolution. Parasitology. 1949, 39: 222-234. 10.1017/S0031182000083785.

    CAS  Article  PubMed  Google Scholar 

  44. 44.

    Dougherty EC: Evolution of zooparasitic groups in the phylum Nematoda, with special reference to host distribution. J Parasitol. 1951, 31: 353-378.

    Article  Google Scholar 

  45. 45.

    Durette-Desset MC, Beveridge I, Spratt DM: The origins and evolutionary expansion of the Strongylida (Nematoda). Int J Parasitol. 1994, 24: 1139-1165. 10.1016/0020-7519(94)90188-0.

    CAS  Article  PubMed  Google Scholar 

  46. 46.

    Schulz RS: Filogenez nematod podotryada strongilyat i perestroika systemy Metastrongyloidea. Dokl Akad Nauk SSSR. 1951, 8: 293-296.

    Google Scholar 

  47. 47.

    Skrjabin KI: Phylogenetic relationships of nematodes of subclass Phasmidia. Zool Zhurnal. 1941, 20: 327-340.

    Google Scholar 

  48. 48.

    Chilton NB, Huby-Chilton F, Gasser RB, Beveridge I: The evolutionary origins of nematodes within the order Strongylida are related to predilection sites within hosts. Mol Phylogenet Evol. 2006, 40: 118-128. 10.1016/j.ympev.2006.01.003.

    CAS  Article  PubMed  Google Scholar 

  49. 49.

    Carreno RA, Nadler SA: Phylogenetic analysis of the Metastrongyloidea (Nematoda: Strongylida) inferred from ribosomal RNA gene sequences. J Parasitol. 2003, 89: 965-973. 10.1645/GE-76R.

    CAS  Article  PubMed  Google Scholar 

Download references


We thank Professor Ian Beveridge for providing the parasite material and commenting on the draft manuscript. The present study was supported by the Australian Research Council (ARC), the National Health and Medical Research Council (NHMRC) and Melbourne Water Corporation. The Early Career Researcher and Collaborative Research grants (AJ) from The University of Melbourne are gratefully acknowledged.

Author information



Corresponding authors

Correspondence to Abdul Jabbar or Robin B Gasser.

Additional information

Competing interests

The authors declare that they have no competing interests.

Authors’ contributions

RBG & ARJ conceived the project and attracted the funding; AJ carried out molecular laboratory work; NM undertook the bioinformatics analysis; RBG, NM & AJ carried out data analysis and interpretation; RBG, AJ, NM & ARJ wrote the draft manuscript. All authors read and approved the final version of the manuscript.

Authors’ original submitted files for images

Below are the links to the authors’ original submitted files for images.

Authors’ original file for figure 1

Authors’ original file for figure 2

Rights and permissions

This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Reprints and Permissions

About this article

Cite this article

Jabbar, A., Mohandas, N., Jex, A.R. et al. The mitochondrial genome of Protostrongylus rufescens – implications for population and systematic studies. Parasites Vectors 6, 263 (2013).

Download citation


  • Protostrongylus (Nematoda: Strongylida)
  • Lungworm
  • Protostrongylosis
  • Sheep
  • Mitochondrial genome
  • Systematics
  • Epidemiology