The mitochondrial ribosomal RNA genes of the nematodes Caenorhabditis elegans and Ascaris suum: consensus secondary-structure models and conserved nucleotide sets for phylogenetic analysis

Mitochondrial genome annotation with MFannot: a critical analysis of gene identification and gene model prediction

Compared to nuclear genomes, mitochondrial genomes (mitogenomes) are small and usually code for only a few dozen genes. Still, identifying genes and their structure can be challenging and time-consuming. Even automated tools for mitochondrial genome annotation often require manual analysis and curation by skilled experts. The most difficult steps are (i) the structural modelling of intron-containing genes; (ii) the identification and delineation of Group I and II introns; and (iii) the identification of moderately conserved, non-coding RNA (ncRNA) genes specifying 5S rRNAs, tmRNAs and RNase P RNAs. Additional challenges arise through genetic code evolution which can redefine the translational identity of both start and stop codons, thus obscuring protein-coding genes. Further, RNA editing can render gene identification difficult, if not impossible, without additional RNA sequence data. Current automated mito- and plastid-genome annotators are limited as they are typically tailored to specific eukaryotic groups. The MFannot annotator we developed is unique in its applicability to a broad taxonomic scope, its accuracy in gene model inference, and its capabilities in intron identification and classification. The pipeline leverages curated profile Hidden Markov Models (HMMs), covariance (CMs) and ERPIN models to better capture evolutionarily conserved signatures in the primary sequence (HMMs and CMs) as well as secondary structure (CMs and ERPIN). Here we formally describe MFannot, which has been available as a web-accessible service (https://megasun.bch.umontreal.ca/apps/mfannot/) to the research community for nearly 16 years. Further, we report its performance on particularly intron-rich mitogenomes and describe ongoing and future developments.

Adaptation of the Romanomermis culicivorax CCA-Adding Enzyme to Miniaturized Armless tRNA Substrates

The mitochondrial genome of the nematode Romanomermis culicivorax encodes for miniaturized hairpin-like tRNA molecules that lack D- as well as T-arms, strongly deviating from the consensus cloverleaf. The single tRNA nucleotidyltransferase of this organism is fully active on armless tRNAs, while the human counterpart is not able to add a complete CCA-end. Transplanting single regions of the Romanomermis enzyme into the human counterpart, we identified a beta-turn element of the catalytic core that-when inserted into the human enzyme-confers full CCA-adding activity on armless tRNAs. This region, originally identified to position the 3'-end of the tRNA primer in the catalytic core, dramatically increases the enzyme's substrate affinity. While conventional tRNA substrates bind to the enzyme by interactions with the T-arm, this is not possible in the case of armless tRNAs, and the strong contribution of the beta-turn compensates for an otherwise too weak interaction required for the addition of a complete CCA-terminus. This compensation demonstrates the remarkable evolutionary plasticity of the catalytic core elements of this enzyme to adapt to unconventional tRNA substrates.

The complete mitochondrial genome of Koerneria sudhausi (Diplogasteromorpha: Nematoda) supports monophyly of Diplogasteromorpha within Rhabditomorpha

Testing hypotheses of monophyly for different nematode groups in the context of broad representation of nematode diversity is central to understanding the patterns and processes of nematode evolution. Herein sequence information from mitochondrial genomes is used to test the monophyly of diplogasterids, which includes an important nematode model organism. The complete mitochondrial genome sequence of Koerneria sudhausi, a representative of Diplogasteromorpha, was determined and used for phylogenetic analyses along with 60 other nematode species. The mtDNA of K. sudhausi is comprised of 16,005 bp that includes 36 genes (12 protein-coding genes, 2 ribosomal RNA genes and 22 transfer RNA genes) encoded in the same direction. Phylogenetic trees inferred from amino acid and nucleotide sequence data for the 12 protein-coding genes strongly supported the sister relationship of K. sudhausi with Pristionchus pacificus, supporting Diplogasteromorpha. The gene order of K. sudhausi is identical to that most commonly found in members of the Rhabditomorpha + Ascaridomorpha + Diplogasteromorpha clade, with an exception of some tRNA translocations. Both the gene order pattern and sequence-based phylogenetic analyses support a close relationship between the diplogasterid species and Rhabditomorpha. The nesting of the two diplogasteromorph species within Rhabditomorpha is consistent with most molecular phylogenies for the group, but inconsistent with certain morphology-based hypotheses that asserted phylogenetic affinity between diplogasteromorphs and tylenchomorphs. Phylogenetic analysis of mitochondrial genome sequences strongly supports monophyly of the diplogasteromorpha.

The complete mitochondrial genome of Meloidogyne graminicola (Tylenchina): a unique gene arrangement and its phylogenetic implications

Meloidogyne graminicola is one of the most economically important plant parasitic-nematodes (PPNs). In the present study, we determined the complete mitochondrial (mt) DNA genome sequence of this plant pathogen. Compared with other PPNs genera, this genome (19,589 bp) is only slightly smaller than that of Pratylenchus vulnus (21,656 bp). The nucleotide composition of the whole mtDNA sequence of M. graminicola is significantly biased toward A and T, with T being the most favored nucleotide and C being the least favored. The A+T content of the entire genome is 83.51%. The mt genome of M. graminicola contains 36 genes (lacking atp8) that are transcribed in the same direction. The gene arrangement of the mt genome of M. graminicola is unique. A total of 21 out of 22 tRNAs possess a DHU loop only, while tRNASer(AGN) lacks a DHU loop. The two large noncoding regions (2,031 bp and 5,063 bp) are disrupted by tRNASer(UCN). Phylogenetic analysis based on concatenated amino acid sequences of 12 protein-coding genes support the monophylies of the three orders Rhabditida, Mermithida and Trichinellida, the suborder Rhabditina and the three infraorders Spiruromorpha, Oxyuridomorpha and Ascaridomorpha, but do not support the monophylies of the two suborders Spirurina and Tylenchina, and the three infraorders Rhabditomorpha, Panagrolaimomorpha and Tylenchomorpha. The four Tylenchomorpha species including M. graminicola, P. vulnus, H. glycines and R. similis from the superfamily Tylenchoidea are placed within a well-supported monophyletic clade, but far from the other two Tylenchomorpha species B. xylophilus and B. mucronatus of Aphelenchoidea. In the clade of Tylenchoidea, M. graminicola is sister to P. vulnus, and H. glycines is sister to R. similis, which suggests root-knot nematodes has a closer relationship to Pratylenchidae nematodes than to cyst nematodes.

The complete mitochondrial genome of Aphelenchoides besseyi (Nematoda: Aphelenchoididae), the first sequenced representative of the subfamily Aphelenchoidinae

The complete mitochondrial genome (mitogenome) ofAphelenchoides besseyiis 16 216 bp in size and has the typical organisation of nematode mitogenomes of Chromadorea, including 12 protein-coding genes (PCGs), two rRNA genes, 22 tRNA genes and the AT-rich non-coding region. The nucleotide composition of the mitogenome ofA. besseyiis AT-biased (80.0%) and the AT skew is −0.289. The most common start codon forA. besseyiis ATT. Thenad3andnad4Lgenes have an incomplete stop codon consisting of just a T and the other PCGs stop with the full stop codons. All the tRNA genes display a non-typical cloverleaf structure of mitochondrial tRNA. The AT-rich non-coding region contains ten tandem repeat units with four different regions. Phylogenetic analysis based on concatenated amino acid sequences of 12 protein-coding genes showed that three Tylenchomorpha species, includingA. besseyi,Bursaphelenchus mucronatusandB. xylophilusfrom the superfamily Aphelenchoidea, are placed within a well-supported monophyletic clade, but far from the other six Tylenchomorpha speciesMeloidogyne chitwoodi,M. graminicola,M. incognita,Pratylenchus vulnus,Heterodera glycinesandRadopholus similisof Tylenchoidea. This phylogeny suggests thatAphelenchoideshas a close relative relationship withBursaphelenchusand that the Tylenchomorpha is not monophyletic.

Complete mitochondrial genomes of chimpanzee- and gibbon-derived Ascaris isolated from a zoological garden in southwest China

Roundworms (Ascaridida: Nematoda), one of the most common soil-transmitted helminths (STHs), can cause ascariasis in various hosts worldwide, ranging from wild to domestic animals and humans. Despite the veterinary and health importance of the Ascaridida species, little or no attention has been paid to roundworms infecting wild animals including non-human primates due to the current taxon sampling and survey bias in this order. Importantly, there has been considerable controversy over the years as to whether Ascaris species infecting non-human primates are the same as or distinct from Ascaris lumbricoides infecting humans. Herein, we first characterized the complete mitochondrial genomes of two representative Ascaris isolates derived from two non-human primates, namely, chimpanzees (Pan troglodytes) and gibbons (Hylobates hoolock), in a zoological garden of southwest China and compared them with those of A. lumbricoides and the congeneric Ascaris suum as well as other related species in the same order, and then used comparative mitogenomics, genome-wide nucleotide sequence identity analysis, and phylogeny to determine whether the parasites from chimpanzees and gibbons represent a single species and share genetic similarity with A. lumbricoides. Taken together, our results yielded strong statistical support for the hypothesis that the chimpanzee- and gibbon-derived Ascaris represent a single species that is genetically similar to A. lumbricoides, consistent with the results of previous morphological and molecular studies. Our finding should enhance public alertness to roundworms originating from chimpanzees and gibbons and the mtDNA data presented here also serves to enrich the resource of markers that can be used in molecular diagnostic, systematic, population genetic, and evolutionary biological studies of parasitic nematodes from either wild or domestic hosts.

The complete mitochondrial genome sequence of the filarial nematode Wuchereria bancrofti from three geographic isolates provides evidence of complex demographic history

Mitochondrial (mt) genome sequences have enabled comparison of population genetics and evolution for numerous free-living and parasitic nematodes. Here we define the complete mt genome of Wuchereria bancrofti through analysis of isolates from Papua New Guinea, India and West Africa. Sequences were assembled for each isolate and annotated with reference to the mt genome sequence for Brugia malayi. The length of the W. bancrofti mt genome is approximately 13,637 nucleotides, contains 2 ribosomal RNAs (rrns), 22 transfer RNAs (trns), 12 protein-coding genes, and is characterized by a 74.6% AT content. The W. bancrofti mt gene order is identical to that reported for Onchocerca volvulus, Dirofilaria immitis, Setaria digitata and B. malayi. In addition to using translational start codons identified previously in the mt protein-coding genes of other filarial nematodes, W. bancrofti appears to be unique in using TGT as a translational start codon. Similarly, use of incomplete stop codons in mt protein-coding genes appears to be more common in W. bancrofti than in other human filarial parasites. The complete mt genome sequence reported here provides new genetic markers for investigating phylogenetic and geographic relationships between isolates, and assessing population diversity within endemic regions. The sequence polymorphism enables new strategies to monitor the progress of public health interventions to control and eliminate this important human parasite. We illustrate the utility of this sequence and single nucleotide polymorphisms by inferring the divergence times between the three W. bancrofti isolates, suggesting predictions into their origin and migration.

The mitochondrial genome of Baylisascaris procyonis

Background

Baylisascaris procyonis (Nematoda: Ascaridida), an intestinal nematode of raccoons, is emerging as an important helminthic zoonosis due to serious or fatal larval migrans in animals and humans. Despite its significant veterinary and public health impact, the epidemiology, molecular ecology and population genetics of this parasite remain largely unexplored. Mitochondrial (mt) genomes can provide a foundation for investigations in these areas and assist in the diagnosis and control of B. procyonis. In this study, the first complete mt genome sequence of B. procyonis was determined using a polymerase chain reaction (PCR)-based primer-walking strategy.

Methodology/Principal Findings

The circular mt genome (14781 bp) of B. procyonis contained 12 protein-coding, 22 transfer RNA and 2 ribosomal RNA genes congruent with other chromadorean nematodes. Interestingly, the B. procyonis mtDNA featured an extremely long AT-rich region (1375 bp) and a high number of intergenic spacers (17), making it unique compared with other secernentean nematodes characterized to date. Additionally, the entire genome displayed notable levels of AT skew and GC skew. Based on pairwise comparisons and sliding window analysis of mt genes among the available 11 Ascaridida mtDNAs, new primer pairs were designed to amplify specific short fragments of the genes cytb (548 bp fragment) and rrnL (200 bp fragment) in the B. procyonis mtDNA, and tested as possible alternatives to existing mt molecular beacons for Ascaridida. Finally, phylogenetic analysis of mtDNAs provided novel estimates of the interrelationships of Baylisasaris and Ascaridida.

Conclusions/Significance

The complete mt genome sequence of B. procyonis sequenced here should contribute to molecular diagnostic methods, epidemiological investigations and ecological studies of B. procyonis and other related ascaridoids. The information will be important in refining the phylogenetic relationships within the order Ascaridida and enriching the resource of markers for systematic, population genetic and evolutionary biological studies of parasitic nematodes of socio-economic importance.

A unique genetic code change in the mitochondrial genome of the parasitic nematode Radopholus similis

BACKGROUND

Mitochondria (mt) contain their own autonomously replicating DNA, constituted as a small circular genome encoding essential subunits of the respiratory chain. Mt DNA is characterized by a genetic code which differs from the standard one. Interestingly, the mt genome of nematodes share some peculiar features, such as small transfer RNAs, truncated ribosomal RNAs and - in the class of Chromadorean nematodes - unidirectional transcription.

FINDINGS

We present the complete mt genomic sequence (16,791 bp) of the plant-parasitic nematode Radopholus similis (class Chromadorea). Although it has a gene content similar to most other nematodes, many idiosyncrasies characterize the extremely AT-rich mt genome of R. similis (85.4% AT). The secondary structure of the large (16S) rRNA is further reduced, the gene order is unique, the large non-coding region contains two large repeats, and most interestingly, the UAA codon is reassigned from translation termination to tyrosine. In addition, 7 out of 12 protein-coding genes lack a canonical stop codon and analysis of transcriptional data showed the absence of polyadenylation. Northern blot analysis confirmed that only one strand is transcribed and processed. Furthermore, using nucleotide content bias methods, regions for the origin of replication are suggested.

CONCLUSION

The extraordinary mt genome of R. similis with its unique genetic code appears to contain exceptional features correlated to DNA decoding. Therefore the genome may provide an incentive to further elucidate these barely understood processes in nematodes. This comprehension may eventually lead to parasitic nematode-specific control targets as healthy mitochondria are imperative for organism survival. In addition, the presented genome is an interesting exceptional event in genetic code evolution.

Evolution of the mitochondrial genome of Metazoa as exemplified by comparison of congeneric species

The mitochondrial genome (mtDNA) of Metazoa is a good model system for evolutionary genomic studies and the availability of more than 1000 sequences provides an almost unique opportunity to decode the mechanisms of genome evolution over a large phylogenetic range. In this paper, we review several structural features of the metazoan mtDNA, such as gene content, genome size, genome architecture and the new parameter of gene strand asymmetry in a phylogenetic framework. The data reviewed here show that: (1) the plasticity of Metazoa mtDNA is higher than previously thought and mainly due to variation in number and location of tRNA genes; (2) an exceptional trend towards stabilization of genomic features occurred in deuterostomes and was exacerbated in vertebrates, where gene content, genome architecture and gene strand asymmetry are almost invariant. Only tunicates exhibit a very high degree of genome variability comparable to that found outside deuterostomes. In order to analyse the genomic evolutionary process at short evolutionary distances, we have also compared mtDNAs of species belonging to the same genus: the variability observed in congeneric species significantly recapitulates the evolutionary dynamics observed at higher taxonomic ranks, especially for taxa showing high levels of genome plasticity and/or fast nucleotide substitution rates. Thus, the analysis of congeneric species promises to be a valuable approach for the assessment of the mtDNA evolutionary trend in poorly or not yet sampled metazoan groups.

The mitochondrial genome of Priapulus caudatus Lamarck (Priapulida: Priapulidae)

We sequenced and annotated the complete mitochondrial (mt) genome of the priapulid Priapulus caudatus in order to provide a source of phylogenetic characters including an assessment of gene order arrangement. The genome was 14,919 bp in its entirety with few, short non-coding regions. A number of protein-coding and tRNA genes overlapped, making the genome relatively compact. The gene order was: cox1, cox2, trnK, trnD, atp8, atp6, cox3, trnG, nad3, trnA, trnR, trnN, rrnS, trnV, rrnL, trnL(yaa), trnL(nag), nad1, -trnS(nga), -cob, -nad6, trnP, -trnT, nad4L, nad4, trnH, nad5, trnF, -trnE, -trnS(nct), trnI, -trnQ, trnM, nad2, trnW, -trnC, -trnY; where '-' indicates genes transcribed on the opposite strand. The gene order, although unique amongst Metazoa, shared the greatest number of gene boundaries and the longest contiguous fragments with the chelicerate Limulus polyphemus. The mt genomes of these taxa differed only by a single inversion of 18 contiguous genes bounded by rrnS and trnS(nct). Other arthropods and nematodes shared fewer gene boundaries but considerably more than the most similar non-ecdysozoan.

The complete mitochondrial genome of the entomopathogenic nematode Steinernema carpocapsae: insights into nematode mitochondrial DNA evolution and phylogeny

We determined the complete sequence of the mitochondrial DNA of the entomopathogenic nematode Steinernema carpocapsae and analyzed its structure and composition as well as the secondary structures predicted for its tRNAs and rRNAs. Almost the complete genome has been amplified in one fragment with long PCR and sequenced using a shotgun strategy. The 13,925-bp genome contains genes for 2 rRNAs, 22 tRNAs, and 12 proteins and lacks an ORF encoding ATPase subunit 8. Four initiation codons were inferred, TTT, TTA, ATA, and ATT, most of the genes ended with TAA or TAG, and only two had a T as an incomplete stop codon. All predicted tRNAs showed the nonconventional secondary structure typical of Secernentea. Although we were able to fold the sequences of trnN, trnD, and trnC into more conventional cloverleaf structures after adding adjacent nucleotides, northern blot experiments showed that the nonstandard tRNAs are actually expressed. Phylogenetic and comparative analyses showed that the mitochondrial genome of S. carpocapsae is more closely related to the genomes of A. suum and C. elegans than to that of Strongyloides stercoralis. This finding does not support the phylogeny based on nuclear small subunit ribosomal DNA sequences previously published. This discrepancy may result from differential reproductive strategies and/or differential selective pressure acting on nuclear and mitochondrial genes. The distinctive characteristics observed among mitochondrial genomes of Secernentea may have arisen to counteract the deleterious effects of Muller's ratchet, which is probably enhanced by the reproductive strategies and selective pressures referred to above.

Isolation and Physiochemical Properties of Protein-Rich Nematode Mitochondrial Ribosomes

In the present study, mitochondrial ribosomes of the nematode Ascaris suum were isolated and their physiochemical properties were compared to ribosomes of Escherichia coli. The sedimentation coefficient and buoyant density of A. suum mitochondrial ribosomes were determined. The sedimentation coefficient of the intact monosome was about 55 S. The buoyant density of formaldehyde-fixed ribosomes in cesium chloride was 1.40 g/cm(3), which suggests that the nematode mitoribosomes have a much higher protein composition than other mitoribosomes. The diffusion coefficients obtained from dynamic light scattering measurements were (1.48 +/- 0.04) x 10(-)(7) cm(2) s(-)(1) for 55 S mitoribosomes and (1.74 +/- 0.04) x 10(-)(7) cm(2) s(-)(1) for the 70 S E. coli monosome. The diameter of mitoribosomes was measured by dynamic light-scattering analysis and electron microscopy. Though the nematode mitoribosome has a larger size than the bacterial ribosome, it does not differ significantly in size from mammalian mitoribosomes.

Complete nucleotide sequence of the mitochondrial genome of Doliolum nationalis with implications for evolution of urochordates

The evolutionary history of the diverse lifestyles adopted by urochordates has attracted intense interest because it may effect the evolutionary history of vertebrates. Here, we report the complete mitochondrial (mt) DNA sequence of the pelagic thaliacean doliolid Doliolum nationalis. The doliolid mt genome shares the unusual tRNAs of trnM(uau) and trnG(ucu) with other ascidians, such as Halocynthia and Ciona. On the other hand, the gene order of the doliolid mt genome is significantly different from that of any ascidian species or vertebrate reported to date. Phylogenetic analyses of the amino acid sequences of 12 protein-coding genes strongly support the sister-grouping of doliolids and the Phlebobranch ascidian Ciona, with the Stolidobranch ascidian alocynthia as the outgroup, thereby providing strong support for the paraphyly of ascidians, as has been suggested by 18S rDNA studies. Given the paraphyletic nature of ascidians, it seems likely that the common ancestor of ascidians and thaliaceans was sessile, as are the present-day ascidians, and that the thaliaceans subsequently evolved a pelagic lifestyle.

Phylogenetic analyses among octocorals (Cnidaria): mitochondrial and nuclear DNA sequences (lsu-rRNA, 16S and ssu-rRNA, 18S) support two convergent clades of branching gorgonians

Gorgonian octocorals lack corroborated hypotheses of phylogeny. This study reconstructs genealogical relationships among some octocoral species based on published DNA sequences from the large ribosomal subunit of the mitochondrial RNA (lsu-rRNA, 16S: 524bp and 21 species) and the small subunit of the nuclear RNA (ssu-rRNA, 18S: 1815bp and 13 spp) using information from insertions-deletions (INDELS) and the predicted secondary structure of the lsu-rRNA (16S). There were seven short (3-10bp) INDELS in the 18S with consistent phylogenetic information. The INDELS in the 16S corresponded to informative signature sequences homologous to the G13 helix found in Escherichia coli. We found two main groups of gorgonian octocorals using a maximum parsimony analysis of the two genes. One group corresponds to deep-water taxa including species from the suborders Calcaxonia and Scleraxonia characterized by an enlargement of the G13 helix. The second group has species from Alcyoniina, Holaxonia and again Scleraxonia characterized by insertions in the 18S. Gorgonian corals, branching colonies with a gorgonin-containing flexible multilayered axis (Holaxonia and Calcaxonia), do not form a monophyletic group. These corroborated results from maternally inherited (16S) and biparentally inherited (18S) genes support a hypothesis of independent evolution of branching in the two octocoral clades.

Molecular identification of natural hybrids between Trichinella nativa and Trichinella T6 provides evidence of gene flow and ongoing genetic divergence

To date, there are no data available on the population genetics of Trichinella due to the lack of genetic markers and the difficulty of working with such small parasites. In the Arctic region of North America and along the Rocky Mountains, there exist two genotypes of Trichinella, Trichinella nativa and Trichinella T6, respectively, which are well differentiated by biochemical and molecular characters. However, both are resistant to freezing, show other common biological characters (e.g. low or no infectivity to rodents and swine) and produce fertile F1 offspring upon interbreeding. To data, these two genotypes have been considered allopatric. In this study, we detected both genotypes in wolves of the same wolf packs in Alaska, suggesting sympatry. A single GTT trinucleotide present in the ITS-2 sequence of T. nativa but not in Trichinella T6 was used as a genetic marker to study gene flow for this character in both a murine infection model and in larvae from naturally-infected Alaskan wolves. Only F1 larvae originating from a cross between T. nativa male and Trichinella T6 female were able to produce F2 offspring. Larvae (F1) originating from a cross between Trichinella T6 male and T. nativa female were not reproductively viable. As expected, all F1 larvae showed a heterozygote pattern for the GTT character upon heteroduplex analysis; however, within the F2 population, the number of observed heterozygotes (n=52) was substantially higher than expected (n=39.08), as supported by the F(is) index, and was not in the Hardy-Weinberg equilibrium. Larvae from two of the 16 Trichinella positive Alaskan wolves, showed the Trichinella T6 pattern or the T. nativa/Trichinella T6 hybrid pattern. Our data demonstrate that T. nativa and Trichinella T6 live in sympatry at least in Alaskan wolves, where T. nativa occurs more frequently (69%) than Trichinella T6 (31%). One explanation for this phenomenon is that glacial periods may have caused a geographical relocation, colonisation and independent evolution of T. nativa within the Rocky Mountains, resulting in a bifurcation of the freeze-resistant genotype. Additional studies will be required to test this hypothesis.

Molecular prospecting for cryptic species of nematodes: mitochondrial DNA versus internal transcribed spacer

DNA sequence divergence at internal transcribed spacer regions (ITS-1 and ITS-2) was compared with divergence at mitochondrial cox1 or nad4 loci in pairs of congeneric nematode species. Mitochondrial sequences accumulate substitutions much more quickly than internal transcribed spacer, the difference being most striking in the most closely related species pairs. Thus, mitochondrial DNA may be the best choice for applications in which one is using sequence data on small numbers of individuals to search for potential cryptic species. On the other hand, internal transcribed spacer remains an excellent tool for DNA diagnostics (quickly distinguishing between known species) owing to its lower level of intraspecific polymorphism.

The mitochondrial genomes of the human hookworms, Ancylostoma duodenale and Necator americanus (Nematoda: Secernentea)

The complete mitochondrial genome sequences were determined for two species of human hookworms, Ancylostoma duodenale (13,721 bp) and Necator americanus (13,604 bp). The circular hookworm genomes are amongst the smallest reported to date for any metazoan organism. Their relatively small size relates mainly to a reduced length in the AT-rich region. Both hookworm genomes encode 12 protein, two ribosomal RNA and 22 transfer RNA genes, but lack the ATP synthetase subunit 8 gene, which is consistent with three other species of Secernentea studied to date. All genes are transcribed in the same direction and have a nucleotide composition high in A and T, but low in G and C. The AT bias had a significant effect on both the codon usage pattern and amino acid composition of proteins. For both hookworm species, genes were arranged in the same order as for Caenorhabditis elegans, except for the presence of a non-coding region between genes nad3 and nad5. In A. duodenale, this non-coding region is predicted to form a stem-and-loop structure which is not present in N. americanus. The mitochondrial genome structure for both hookworms differs from Ascaris suum only in the location of the AT-rich region, whereas there are substantial differences when compared with Onchocerca volvulus, including four gene or gene-block translocations and the positions of some transfer RNA genes and the AT-rich region. Based on genome organisation and amino acid sequence identity, A. duodenale and N. americanus were more closely related to C. elegans than to A. suum or O. volvulus (all secernentean nematodes), consistent with a previous phylogenetic study using ribosomal DNA sequence data. Determination of the complete mitochondrial genome sequences for two human hookworms (the first members of the order Strongylida ever sequenced) provides a foundation for studying the systematics, population genetics and ecology of these and other nematodes of socio-economic importance.

Mitochondrial DNA variation in benzimidazole-resistant and -susceptible populations of the small ruminant parasite Teladorsagia circumcincta

The genetic diversity of the mtDNA ND4 gene in 11 Teladorsagia circumcincta populations from France and Morocco was assessed by sequencing. Some of these nematode populations were resistant to benzimidazole (BZ) anthelmintics, while others were susceptible. The nucleotide diversity in all populations studied was very high, probably due to a high mutation rate in nematodes, but there was no significant difference between them. This suggests that no strong, recurrent bottlenecks occur during the acquisition of BZ resistance by a worm population. The conservation of genetic variations during the acquisition of BZ resistance is probably due to the fact that anthelmintic treatments do not kill all the susceptible adult worms and to the presence of numerous free-living larvae that are not submitted to this anthelmintic pressure. There was no genetic subdivision between worm populations on a small geographical scale (less than 200 km), but significant F(ST)s were found on a larger geographical scale. This kind of subdivision cannot be explained by different genetic flows between populations because all these populations were isolated from each other. This subdivision is probably due to the breeding management practices and the large size of these worm populations, which limit genetic drift.

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

The complete mitochondrial DNA (mtDNA) of the nematode Trichinella spiralis has been amplified in four overlapping fragments and 16,656 bp of its sequence has been determined. This sequence contains the 37 genes typical of metazoan mtDNAs, including a putative atp8, which is absent from all other nematode mtDNAs examined. The genes are transcribed from both mtDNA strands and have an arrangement relatable to those of coelomate metazoans, but not to those of secernentean nematodes. All protein genes appear to initiate with ATN codons, typical for metazoans. Neither TTG nor GTT start codons, inferred for several genes of other nematodes, were found. The 22 T. spiralis tRNA genes fall into three categories: (i) those with the potential to form conventional "cloverleaf" secondary structures, (ii) those with TPsiC arm + variable arm replacement loops, and (iii) those with DHU-arm replacement loops. Mt-tRNA(R) has a 5'-UCG-3' anticodon, as in most other metazoans, instead of the very unusual 5'-ACG-3' present in the secernentean nematodes. The sequence also contains a large repeat region that is polymorphic in size at the population and/or individual level.

The filarial genome project: analysis of the nuclear, mitochondrial and endosymbiont genomes of Brugia malayi

The Filarial Genome Project (FGP) was initiated in 1994 under the auspices of the World Health Organisation. Brugia malayi was chosen as the model organism due to the availability of all life cycle stages for the construction of cDNA libraries. To date, over 20000 cDNA clones have been partially sequenced and submitted to the EST database (dbEST). These ESTs define approximately 7000 new Brugia genes. Analysis of the EST dataset provides useful information on the expression pattern of the most abundantly expressed Brugia genes. Some highly expressed genes have been identified that are expressed in all stages of the parasite's life cycle, while other highly expressed genes appear to be stage-specific. To elucidate the structure of the Brugia genome and to provide a basis for comparison to the Caenorhabditis elegans genome, the FGP is also constructing a physical map of the Brugia chromosomes and is sequencing genomic BAC clones. In addition to the nuclear genome, B. malayi possesses two other genomes: the mitochondrial genome and the genome of a bacterial endosymbiont. Eighty percent of the mitochondrial genome of B. malayi has been sequenced and is being compared to mitochondrial sequences of other nematodes. The bacterial endosymbiont genome found in B. malayi is closely related to the Wolbachia group of rickettsia-like bacteria that infects many insect species. A set of overlapping BAC clones is being assembled to cover the entire bacterial genome. Currently, half of the bacterial genome has been assembled into four contigs. A consortium has been established to sequence the entire genome of the Brugia endosymbiont. The sequence and mapping data provided by the FGP is being utilised by the nematode research community to develop a better understanding of the biology of filarial parasites and to identify new vaccine candidates and drug targets to aid the elimination of human filariasis.

Complete DNA sequence of the mitochondrial genome of the ascidian Halocynthia roretzi (Chordata, Urochordata)

The complete nucleotide sequence of the 14,771-bp-long mitochondrial (mt) DNA of a urochordate (Chordata)-the ascidian Halocynthia roretzi-was determined. All the Halocynthia mt-genes were found to be located on a single strand, which is rich in T and G rather than in A and C. Like nematode and Mytilus edulis mtDNAs, that of Halocynthia encodes no ATP synthetase subunit 8 gene. However, it does encode an additional tRNA gene for glycine (anticodon TCT) that enables Halocynthia mitochondria to use AGA and AGG codons for glycine. The mtDNA carries an unusual tRNA(Met) gene with a TAT anticodon instead of the usual tRNA(Met)(CAT) gene. As in other metazoan mtDNAs, there is not any long noncoding region. The gene order of Halocynthia mtDNA is completely different from that of vertebrate mtDNAs except for tRNA(His)-tRNA(Ser)(GCU), suggesting that evolutionary change in the mt-gene structure is much accelerated in the urochordate line compared with that in vertebrates. The amino acid sequences of Halocynthia mt-proteins deduced from their gene sequences are quite different from those in other metazoans, indicating that the substitution rate in Halocynthia mt-protein genes is also accelerated.

Comparison of the complete protein sets of worm and yeast: orthology and divergence

Comparative analysis of predicted protein sequences encoded by the genomes of Caenorhabditis elegans and Saccharomyces cerevisiae suggests that most of the core biological functions are carried out by orthologous proteins (proteins of different species that can be traced back to a common ancestor) that occur in comparable numbers. The specialized processes of signal transduction and regulatory control that are unique to the multicellular worm appear to use novel proteins, many of which re-use conserved domains. Major expansion of the number of some of these domains seen in the worm may have contributed to the advent of multicellularity. The proteins conserved in yeast and worm are likely to have orthologs throughout eukaryotes; in contrast, the proteins unique to the worm may well define metazoans.

The mitochondrial genome of the sea anemone Metridium senile (Cnidaria): introns, a paucity of tRNA genes, and a near-standard genetic code

The circular, 17,443 nucleotide-pair mitochondrial (mt) DNA molecule of the sea anemone, Metridium senile (class Anthozoa, phylum Cnidaria) is presented. This molecule contains genes for 13 energy pathway proteins and two ribosomal (r) RNAs but, relative to other metazoan mtDNAs, has two unique features: only two transfer RNAs (tRNA(f-Met) and tRNA(Trp)) are encoded, and the cytochrome c oxidase subunit I (COI) and NADH dehydrogenase subunit 5 (ND5) genes each include a group I intron. The COI intron encodes a putative homing endonuclease, and the ND5 intron contains the molecule's ND1 and ND3 genes. Most of the unusual characteristics of other metazoan mtDNAs are not found in M. senile mtDNA: unorthodox translation initiation codons and partial translation termination codons are absent, the use of TGA to specify tryptophan is the only genetic code modification, and both encoded tRNAs have primary and secondary structures closely resembling those of standard tRNAs. Also, with regard to size and secondary structure potential, the mt-s-rRNA and mt-1-rRNA have the least deviation from Escherichia coli 16S and 23S rRNAs of all known metazoan mt-rRNAs. These observations indicate that most of the genetic variations previously reported in metazoan mtDNAs developed after Cnidaria diverged from the common ancestral line of all other Metazoa.

Molecular evolution of a portion of the mitochondrial 16S ribosomal gene region in scleractinian corals

Relationships among families and suborders of scleractinian corals are poorly understood because of difficulties 1) in making inferences about the evolution of the morphological characters used in coral taxonomy and 2) in interpreting their 240-million-year fossil record. Here we describe patterns of molecular evolution in a segment of the mitochondrial (mt) 16S ribosomal gene from taxa of 14 families of corals and the use of this gene segment in a phylogenetic analysis of relationships within the order. We show that sequences obtained from scleractinians are homologous to other metazoan 16S ribosomal sequences and fall into-two distinct clades defined by size of the amplified gene product. Comparisons of sequences from the two clades demonstrate that both sets of sequences are evolving under similar evolutionary constraints: they do not differ in nucleotide composition, numbers of transition and transversion substitutions, spatial patterns of substitutions, or in rates of divergence. The characteristics and patterns observed in these sequences as well as the secondary structures, are similar to those observed in mt 16S ribosomal DNA sequences from other taxa. Phylogenetic analysis of these sequences shows that they are useful for evaluating relationships within the order. The hypothesis generated from this analysis differs from traditional hypotheses for evolutionary relationships among the Scleractinia and suggests that a reevaluation of evolutionary affinities in the order is needed.

The sequence and structure of the meadow grasshopper (Chorthippus parallelus) mitochondrial srRNA, ND2, COI, COII ATPase8 and 9 tRNA genes

The nucleotide sequence of the mitochondrial ND2, COI, COII, ATPase8, srRNA and nine tRNA genes have been sequenced from two individual of the meadow grasshopper Chorthippus parallelus. Comparisons are made to other insects for which the same regions are completely sequenced. Percentage A + T is found to be relatively low in C. parallelus though consistent with that of the other Orthopteran, Locusta migratoria. The relative number of substitutions observed in the different protein-coding genes was analysed between pairs of insect species sharing different levels of relatedness. A clear change in this rate was observed between the within-genus and between-genera comparisons. This change is interpreted in terms of the functional constraints acting on these four different genes. The patterns seem to result from an early saturation of COI and COII genes with synonymous substitutions, and a tolerance of ND2 and ATPase8 function to high levels of amino acid replacements. This analysis highlights a need for further sequence studies and comparisons between taxa of different levels of divergence in order to understand the patterns of mtDNA evolution on which many evolutionary investigations are based.