Mitochondrial genes of Schistosoma mansoni

Mitochondrial genomes of two eucotylids as the first representatives from the superfamily Microphalloidea (Trematoda) and phylogenetic implications

Background

The Eucotylidae Cohn, 1904 (Superfamily: Microphalloidea), is a family of digeneans parasitic in kidneys of birds as adults. The group is characterized by the high level of morphological similarities among genera and unclear systematic value of morphological characters traditionally used for their differentiation. In the present study, we sequenced the complete or nearly complete mitogenomes (mt genome) of two eucotylids representing the genera Tamerlania (T. zarudnyi) and Tanaisia (Tanaisia sp.). They represent the first sequenced mt genomes of any member of the superfamily Microphalloidea.

Methods

A comparative mitogenomic analysis of the two newly sequenced eucotylids was conducted for the investigation of mitochondrial gene arrangement, contents and genetic distance. Phylogenetic position of the family Eucotylidae within the order Plagiorchiida was examined using nucleotide sequences of mitochondrial protein-coding genes (PCGs) plus RNAs using maximum likelihood (ML) and Bayesian inference (BI) methods. BI phylogeny based on concatenated amino acids sequences of PCGs was also conducted to determine possible effects of silent mutations.

Results

The complete mt genome of T. zarudnyi was 16,188 bp and the nearly complete mt genome of Tanaisia sp. was 13,953 bp in length. A long string of additional amino acids (about 123 aa) at the 5′ end of the cox1 gene in both studied eucotylid mt genomes has resulted in the cox1 gene of eucotylids being longer than in all previously sequenced digeneans. The rrnL gene was also longer than previously reported in any digenean mitogenome sequenced so far. The TΨC and DHU loops of the tRNAs varied greatly between the two eucotylids while the anticodon loop was highly conserved. Phylogenetic analyses based on mtDNA nucleotide and amino acids sequences (as a separate set) positioned eucotylids as a sister group to all remaining members of the order Plagiorchiida. Both ML and BI phylogenies revealed the paraphyletic nature of the superfamily Gorgoderoidea and the suborder Xiphidiata.

Conclusions

The average sequence identity, combined nucleotide diversity and Kimura-2 parameter distances between the two eucotylid mitogenomes demonstrated that atp6, nad5, nad4L and nad6 genes are better markers than the traditionally used cox1 or nad1 for the species differentiation and population-level studies of eucotylids because of their higher variability. The position of the Dicrocoeliidae and Eucotylidae outside the clade uniting other xiphidiatan trematodes strengthened the argument for the need for re-evaluation of the taxonomic content of the Xiphidiata.

The Complete Mitochondrial Genome of the Caecal Fluke of Poultry, Postharmostomum commutatum, as the First Representative from the Superfamily Brachylaimoidea

Postharmostomum commutatum (Platyhelminthes: Brachylaimoidea), a parasite of the caeca of poultry, has been frequently reported from many countries and regions, including China. However, the molecular epidemiology, population genetics and phylogenetics of this parasite are poorly understood. In the present study, we determined and characterized the complete mitochondrial (mt) genome of P. commutatum, as the first representative from the superfamily Brachylaimoidea. The mt genome of P. commutatum is a circular DNA molecule of 13,799 bp in size and encodes the complete set of 36 genes (12 protein-coding genes, 22 transfer RNA genes, two ribosomal RNA genes) as well as a typical control region. The mt genome of P. commutatum presents a clear bias in nucleotide composition with a negative AT-skew on average (-0.306) and a positive GC-skew on average (0.466). Phylogenetic analyses showed that P. commutatum (superfamily Brachylaimoidea) and other ten members of the order Diplostomida were recovered as sister groups of the order Plagiorchiida, indicating that the order Diplostomida is paraphyletic. This is the first mt genome of any member of the superfamily Brachylaimoidea and should represent a rich source of genetic markers for molecular epidemiological, population genetic and phylogenetic studies of parasitic flukes of socio-economic importance in poultry.

DNA 'barcoding' of Schistosoma mansoni across sub-Saharan Africa supports substantial within locality diversity and geographical separation of genotypes

Schistosoma mansoni is a widespread human helminth and causes intestinal schistosomiasis in 54 countries, mainly across Africa but also in Madagascar, the Arabian Peninsula and the neotropics. The geographical range of this parasite relies on the distribution of certain species of freshwater pulmonate snails of the genus Biomphalaria. Whilst S. mansoni is known to exhibit high population diversity the true extent of this diversity is still to be fully elucidated as sampling of this taxon progressively accrues. Here a DNA 'barcoding' approach is taken using sequence analysis of a 450bp region within the mitochondrial cox1 gene to assess the genetic diversity within a large number of S. mansoni larval stages collected from their natural human hosts across sub-Saharan Africa. Five hundred and sixty one individual parasite samples were examined from 22 localities and 14 countries. Considerable within-species diversity was found with 120 unique haplotypes splitting geographically into five discrete lineages. The highest diversity was found in East Africa with samples forming three of the five lineages. Less diversity was found in the Far and Central Western regions of Africa with haplotypes from the New World showing a close affinity to the Far Western African S. mansoni populations supporting the hypothesis of a colonisation of South America via the West African slave trade. The data are discussed in relation to parasite diversity and disease epidemiology.

The mitochondrial genome of Paramphistomum cervi (Digenea), the first representative for the family Paramphistomidae

We determined the complete mitochondrial DNA (mtDNA) sequence of a fluke, Paramphistomum cervi (Digenea: Paramphistomidae). This genome (14,014 bp) is slightly larger than that of Clonorchis sinensis (13,875 bp), but smaller than those of other digenean species. The mt genome of P. cervi contains 12 protein-coding genes, 22 transfer RNA genes, 2 ribosomal RNA genes and 2 non-coding regions (NCRs), a complement consistent with those of other digeneans. The arrangement of protein-coding and ribosomal RNA genes in the P. cervi mitochondrial genome is identical to that of other digeneans except for a group of Schistosoma species that exhibit a derived arrangement. The positions of some transfer RNA genes differ. Bayesian phylogenetic analyses, based on concatenated nucleotide sequences and amino-acid sequences of the 12 protein-coding genes, placed P. cervi within the Order Plagiorchiida, but relationships depicted within that order were not quite as expected from previous studies. The complete mtDNA sequence of P. cervi provides important genetic markers for diagnostics, ecological and evolutionary studies of digeneans.

Genetic diversity within Schistosoma haematobium: DNA barcoding reveals two distinct groups

Background

Schistosomiasis in one of the most prevalent parasitic diseases, affecting millions of people and animals in developing countries. Amongst the human-infective species S. haematobium is one of the most widespread causing urogenital schistosomiasis, a major human health problem across Africa, however in terms of research this human pathogen has been severely neglected.

Methodology/Principal Findings

To elucidate the genetic diversity of Schistosoma haematobium, a DNA ‘barcoding’ study was performed on parasite material collected from 41 localities representing 18 countries across Africa and the Indian Ocean Islands. Surprisingly low sequence variation was found within the mitochondrial cytochrome oxidase subunit I (cox1) and the NADH-dehydrogenase subunit 1 snad1). The 61 haplotypes found within 1978 individual samples split into two distinct groups; one (Group 1) that is predominately made up of parasites from the African mainland and the other (Group 2) that is made up of samples exclusively from the Indian Ocean Islands and the neighbouring African coastal regions. Within Group 1 there was a dominance of one particular haplotype (H1) representing 1574 (80%) of the samples analyzed. Population genetic diversity increased in samples collected from the East African coastal regions and the data suggest that there has been movement of parasites between these areas and the Indian Ocean Islands.

Conclusions/Significance

The high occurrence of the haplotype (H1) suggests that at some point in the recent evolutionary history of S. haematobium in Africa the population may have passed through a genetic ‘bottleneck’ followed by a population expansion. This study provides novel and extremely interesting insights into the population genetics of S. haematobium on a large geographic scale, which may have consequence for control and monitoring of urogenital schistosomiasis.

Schistosomiasis is a disease caused by parasitic blood flukes of the genus Schistosoma. Species that infect humans are prevalent in developing countries, having a major impact on public health and well-being as well as an impediment to socioeconomic development. More people are infected with Schistosoma haematobium than with all the other schistosome species combined, however mainly due to the inability to maintain S. haematobium in the laboratory system empirical studies on this parasite are minimal. The genetic variation of this Schistosoma species on a wide geographical scale has never been investigated. In this study, we have used a DNA ?barcoding? approach to document the genetic variation and population structure of S. haematobium sampled from 18 countries across Africa and the Indian ocean Islands. The study revealed a distinct genetic separation of S. haematobium from the Indian Ocean Islands and the closely neighbouring coastal regions from S. haematobium found throughout the African mainland, the latter of which exhibited extremely low levels of mitochondrial diversity within and between populations of parasites sampled. The data from this study provides a novel insight into the population genetics of S. haematobium and will have an impact on future research strategies.

Molecular epidemiology ofSchistosoma mansoniin Uganda: DNA barcoding reveals substantial genetic diversity within Lake Albert and Lake Victoria populations

Representative samples of UgandanSchistosoma mansonifrom Lake Albert and Lake Victoria were examined using DNA barcoding, sequence analysis of two partially overlapping regions – ASMIT (396 bp) & MORGAN (617 bp) – of the mitochondrial cytochrome oxidase subunit I (cox1). The Victorian sample exhibited greater nucleotide diversity, 1·4%vs. 1·0%, and a significant population partition appeared as barcodes did not cross-over between lakes. With one exception, Lake Albert populations were more mixed by sampled location, while those from Lake Victoria appeared more secluded. Using statistical parsimony, barcode ASMIT 1 was putatively ancestral to all others and analysis of MORGAN cox1 confirmed population diversity. All samples fell into two of five well-resolved lineages; sub-lineages therein broadly partitioning by lake. It seems that barcode ASMIT 1 (and close variants) was likely widely dispersed throughout the Nilotic environment but later diversifiedin situ, and in parallel, within Lake Albert and Lake Victoria. The genetic uniformity of UgandanS. mansonican no longer be assumed, which might better explain known epidemiological heterogeneities. While it appears plausible that locally evolved heritable traits could spread through most of the Lake Albert populations, it seems unlikely they could quickly homogenise into Lake Victoria or amongst populations therein.

Brazilian studies on the genetics of Schistosoma mansoni

The parasite Schistosoma is known to exhibit variations among species, strains and genera, such as, the levels of infectivity, pathogenicity and immunogenicity. These factors may differ among parasite populations according to the local epidemiological conditions. Diversity observed in Schistosoma mansoni from different geographical regions or within individuals of the same region can be determined by differences in the genotype of each parasite strain. However, until recently, finding adequate genetic markers to investigate infectivity or other epidemiological characteristics of a transmission area proved difficult. Several studies have been conducted to evaluate the genetic variability of S. mansoni, using different techniques. Intraspecific variability was observed in morphological characters, isoenzyme studies, mtDNA, ribosomal gene probes, random amplification of polymorphic DNA (RAPD) and microsatellites. The sequencing of the S. mansoni genome was the most important achievement concerning genetic approaches to the study of this parasite and may improve the development of drugs, vaccines and diagnostics of schistosomiasis. The knowledge of the genetic structure of schistosome populations in relation to epidemiological data and host variability is essential for the understanding of the epidemiology of the disease and the design of control strategies.

Codon usage and bias in mitochondrial genomes of parasitic platyhelminthes

Sequences of the complete protein-coding portions of the mitochondrial (mt) genome were analysed for 6 species of cestodes (including hydatid tapeworms and the pork tapeworm) and 5 species of trematodes (blood flukes and liver- and lung-flukes). A near-complete sequence was also available for an additional trematode (the blood fluke Schistosoma malayensis). All of these parasites belong to a large flatworm taxon named the Neodermata. Considerable variation was found in the base composition of the protein-coding genes among these neodermatans. This variation was reflected in statistically-significant differences in numbers of each inferred amino acid between many pairs of species. Both convergence and divergence in nucleotide, and hence amino acid, composition was noted among groups within the Neodermata. Considerable variation in skew (unequal representation of complementary bases on the same strand) was found among the species studied. A pattern is thus emerging of diversity in the mt genome in neodermatans that may cast light on evolution of mt genomes generally.

Advances and Trends in the Molecular Systematics of the Parasitic Platyhelminthes

The application of molecular systematics to the parasitic Platyhelminthes (Cestoda, Digenea and Monogenea) over the last decade has advanced our understanding of their interrelationships and evolution substantially. Here we review the current state of play and the early works that led to the molecular-based hypotheses that now predominate in the field; advances in their systematics, taxonomy, classification and phylogeny, as well as trends in species circumscription, molecular targets and analytical methods are discussed for each of the three major parasitic groups. A by-product of this effort has been an ever increasing number of parasitic flatworms characterized genetically, and the useful application of these data to the diagnosis of animal and human pathogens, and to the elucidation of life histories are presented. The final section considers future directions in the field, including taxon sampling, molecular targets of choice, and the current and future utility of mitochondrial and nuclear genomics in systematic study.

Molecular taxonomic position of the elephant schistosome, Bivitellobilharzia nairi, newly discovered in Sri Lanka

Bivitellobilharzia nairi (Mudaliar and Ramanujachar, 1945) Dutt and Srivastava, 1955 was first recorded in India. A number of adult worm specimens of this schistosome species were recovered from a domestic elephant, which died in 1999 in Sri Lanka. This is the first report of this schistosome from Sri Lanka. In the present study, in order to clarify the phylogenetic relationship with other species of schistosomes, sequences from the second internal transcribed spacer (ITS2) of the ribosomal gene repeat, part of the 28S ribosomal RNA gene (28S), and part of the mitochondrial cytochrome c oxidase subunit 1 (CO1) gene from B. nairi were analyzed. Two intraspecific variations were seen within 13 individuals in the ITS2 region. In the CO1 region of the mitochondrial DNA, there were four haplotypes in the nucleotide sequences and two haplotypes in the amino acid sequences. Phylogenetic analysis using the nuclear DNA showed that B. nairi was basal to all of species of the genus Schistosoma. The 28S tree also showed that the mammalian lineage was monophyletic. However, phylogenetic analysis using the mitochondrial DNA showed that B. nairi was nested within the genus Schistosoma. The taxonomical position for this species as well as the contradiction between the results from the nuclear and mitochondrial genes were discussed.

Identities of two Paragonimus species from Sri Lanka inferred from molecular sequences

Metacercariae of Paragonimus spp. were obtained from field-collected freshwater crabs in Sri Lanka. Genomic DNA was extracted from single metacercariae. Two gene regions (partial mitochondrial cytochrome c oxidase subunit 1 (CO1) and the second internal transcribed spacer of the nuclear ribosomal gene repeat (ITS2)) were amplified using the polymerase chain reaction. Two differing sequences were obtained for each of these gene regions. Phylogenetic analyses placed the type 1 sequences as sister to a clade containing P. westermani and P. siamensis whereas the type 2 sequences were close to published sequences of P. siamensis from Thailand. The possible taxonomic status of these two types are discussed. This is the first report of molecular data about Paragonimus from Sri Lanka.

A molecular phylogeographic study based on DNA sequences from individual metacercariae of Paragonimus mexicanus from Guatemala and Ecuador

A molecular phylogeographic study of Paragonimus mexicanus collected from Guatemala and Ecuador was performed. Genomic DNA was extracted from individual metacercariae, and two gene regions (partial mitochondrial cytochrome c oxidase subunit 1 (CO1) and the second internal transcribed spacer of the nuclear ribosomal gene repeat (ITS2)) were amplified by the polymerase chain reaction (PCR). Sequences segregated in a phylogenetic tree according to their geographic origins. ITS2 sequences from Ecuador and Guatemala differed at only one site. Pairwise distances among CO1 sequences within a country were always lower than between countries. Nevertheless, genetic distances between countries were less than between geographical forms of P. westermani that have been suggested to be distinct species. This result suggests that populations from Guatemala and Ecuador are genetically differentiated perhaps at the level of subspecies.

Mitochondrial genomes of parasitic flatworms

Complete or near-complete mitochondrial genomes are now available for 11 species or strains of parasitic flatworms belonging to the Trematoda and the Cestoda. The organization of these genomes is not strikingly different from those of other eumetazoans, although one gene (atp8) commonly found in other phyla is absent from flatworms. The gene order in most flatworms has similarities to those seen in higher protostomes such as annelids. However, the gene order has been drastically altered in Schistosoma mansoni, which obscures this possible relationship. Among the sequenced taxa, base composition varies considerably, creating potential difficulties for phylogeny reconstruction. Long non-coding regions are present in all taxa, but these vary in length from only a few hundred to approximately 10000 nucleotides. Among Schistosoma spp., the long non-coding regions are rich in repeats and length variation among individuals is known. Data from mitochondrial genomes are valuable for studies on species identification, phylogenies and biogeography.

Affinities between Asian non-human Schistosoma species, the S. indicum group, and the African human schistosomes

Schistosoma species have traditionally been arranged in groups based on egg morphology, geographical origins, and the genus or family of snail intermediate host. One of these groups is the 'S. indicum group' comprising species from Asia that use pulmonate snails as intermediate hosts. DNA sequences were obtained from the four members of this group (S. indicum, S. spindale, S. nasale and S. incognitum) to provide information concerning their phylogenetic relationships with other Asian and African species and species groups. The sequences came from the second internal transcribed spacer (ITS2) of the ribosomal gene repeat, part of the 28S ribosomal RNA gene (28S), and part of the mitochondrial cytochrome c oxidase subunit 1 (CO1) gene. Tree analyses using both distance and parsimony methods showed the S. indicum group not to be monophyletic. Schistosoma indicum, S. spindale and S. nasale were clustered among African schistosomes, while S. incognitum was placed as sister to the African species (using ITS2 and 28S nucleotide sequences and CO1 amino acid sequences), or as sister to all other species of Schistosoma (CO1 nucleotide sequences). Based on the present molecular data, a scenario for the evolution of the S. indicum group is discussed.

Mitochondrial gene content, arrangement and composition compared in African and Asian schistosomes

Complete sequences were obtained for the coding portions of the mitochondrial (mt) genomes of Schistosoma mansoni (NMRI strain, Puerto Rico; 14 415 bp), S. japonicum (Anhui strain, China; 14 085 bp) and S. mekongi (Khong Island, Laos; 14 072 bp). Each comprises 36 genes: 12 protein-encoding genes (cox1-3, nad1-6, nad4L, atp6 and cob); two ribosomal RNAs, rrnL (large subunit rRNA or 16S) and rrnS (small subunit rRNA or 12S); as well as 22 transfer RNA (tRNA) genes. The atp8 gene is absent. A large segment (9.6 kb) of the coding region (comprising 14 tRNAs, eight complete and two incomplete protein-encoding genes) for S. malayensis (Baling, Malaysian Peninsula) was also obtained. Each genome also possesses a long non-coding region that is divided into two parts (a small and a large non-coding region, the latter not fully sequenced in any species) by one or more tRNAs. The protein-encoding genes are similar in size, composition and codon usage in all species except for cox1 in S. mansoni (609 aa) and cox2 in S. mekongi (219 aa), both of which are longer than homologues in other species. An unexpected finding in all the Schistosoma species was the presence of a leucine zipper motif in the nad4L gene. The gene order in S. mansoni is strikingly different from that seen in the S. japonicum group and other flatworms. There is a high level of identity (87-94% at both the nucleotide and amino acid levels) for all protein-encoding genes of S. mekongi and S. malayensis. The identity between genes of these two species and those of S. japonicum is less (56-83% for amino acids and 73-79% for nucleotides). The identity between the genes of S. mansoni and the Asian schistosomes is far less (33-66% for amino acids and 54-68% for nucleotides), an observation consistent with the known phylogenetic distance between S. mansoni and the other species.

A molecular phylogeny of the genus Ichthyocotylurus (Digenea, Strigeidae)

Three nucleotide data sets, two nuclear (ribosomal internal transcribed spacer regions 1 and 2, ITS1 and ITS2) and one mitochondrial (cytochrome c oxidase subunit 1, CO1), were analysed using distance matrix and maximum likelihood methods to determine the inter-relationships amongst the four species attributed to the genus Ichthyocotylurus Odening, 1969. Sequence data obtained from all gene loci investigated supported the position of Ichthyocotylurus variegatus as a species discrete from Ichthyocotylurus platycephalus. Phylogenetic analyses yielded congruent trees, with I. variegatus isolates comprising a common clade to which I. platycephalus constitutes a sister taxon. Ichthyocotylurus erraticus and Ichthyocotylurus pileatus were found to demonstrate a similarly close inter-specific relationship. The greatest intra-generic divergence occurred in the CO1 region (16% variability), with resultant disparities in three to eight encoded amino acids. PCR amplification yielded multiple ITS1 products for all Ichthyocotylurus spp. Analyses of equivalent-sized amplicons showed 5.4% intra-generic variation and several point mutations between I. variegatus isolates from different geographical localities and from different piscine hosts. The ITS2 locus was extremely conserved, with less than 1% variation between species. No intra-specific variation was recorded for any CO1 or ITS2 sequences.

Complete sequence of the mitochondrial genome of the tapeworm Hymenolepis diminuta: gene arrangements indicate that Platyhelminths are Eutrochozoans

Using "long-PCR," we amplified in overlapping fragments the complete mitochondrial genome of the tapeworm Hymenolepis diminuta (Platyhelminthes: Cestoda) and determined its 13,900-nt sequence. The gene content is the same as that typically found for animal mitochondrial DNA (mtDNA) except that atp8 appears to be lacking, a condition found previously for several other animals. Despite the small size of this mtDNA, there are two large noncoding regions, one of which contains 13 repeats of a 31-nt sequence and a potential stem-loop structure of 25 bp with an 11-member loop. Large potential secondary structures were identified also for the noncoding regions of two other cestode mtDNAS: Comparison of the mitochondrial gene arrangement of H. diminuta with those previously published supports a phylogenetic position of flatworms as members of the Eutrochozoa, rather than placing them basal to either a clade of protostomes or a clade of coelomates.

Mitochondrial genomes of human helminths and their use as markers in population genetics and phylogeny

To date, over 100 complete metazoan mitochondrial (mt) genomes of different phyla have been reported. Here, we briefly summarise mt gene organisation in the Metazoa and review what is known of the mt genomes of nematodes and flatworms parasitic in humans. The availability of complete or almost complete mtDNA sequences for several parasitic helminths provides a rich source of genetic markers for phylogenetic analysis and study of genetic variability in helminth groups. Examples of the application of mtDNA in studies on Ascaris, Onchocerca, Schistosoma, Fasciola, Paragonimus, Echinostoma, Echinococcus and Taenia are described.

Changes in mitochondrial genetic codes as phylogenetic characters: two examples from the flatworms

Shared molecular genetic characteristics other than DNA and protein sequences can provide excellent sources of phylogenetic information, particularly if they are complex and rare and are consequently unlikely to have arisen by chance convergence. We have used two such characters, arising from changes in mitochondrial genetic code, to define a clade within the Platyhelminthes (flatworms), the Rhabditophora. We have sampled 10 distinct classes within the Rhabditophora and find that all have the codon AAA coding for the amino acid Asn rather than the usual Lys and AUA for Ile rather than the usual Met. We find no evidence to support claims that the codon UAA codes for Tyr in the Platyhelminthes rather than the standard stop codon. The Rhabditophora are a very diverse group comprising the majority of the free-living turbellarian taxa and the parasitic Neodermata. In contrast, three other classes of turbellarian flatworm, the Acoela, Nemertodermatida, and Catenulida, have the standard invertebrate assignments for these codons and so are convincingly excluded from the rhabditophoran clade. We have developed a rapid computerized method for analyzing genetic codes and demonstrate the wide phylogenetic distribution of the standard invertebrate code as well as confirming already known metazoan deviations from it (ascidian, vertebrate, echinoderm/hemichordate).

Mitochondrial genome diversity in parasites

Mitochondrial genomes have been sequenced from a wide variety of organisms, including an increasing number of parasites. They maintain some characteristics in common across the spectrum of life-a common core of genes related to mitochondrial respiration being most prominent-but have also developed a great diversity of gene content, organisation, and expression machineries. The characteristics of mitochondrial genomes vary widely among the different groups of protozoan parasites, from the minute genomes of the apicomplexans to amoebae with 20 times as many genes. Kinetoplastid protozoa have a similar number of genes to metazoans, but the details of gene organisation and expression in kinetoplastids require extraordinary mechanisms. Mitochondrial genes in nematodes and trematodes appear quite sedate in comparison, but a closer look shows a strong tendency to unusual tRNA structure and alternative initiation codons among these groups. Mitochondrial genes are increasingly coming into play as aids to phylogenetic and epidemiologic analyses, and mitochondrial functions are being recognised as potential drug targets. In addition, examination of mitochondrial genomes is producing further insights into the diversity of the wide-ranging group of organisms comprising the general category of parasites.

Mitochondrial DNA sequences of human schistosomes: the current status

Sequences generated from the mitochondrial genome provide useful molecular markers for defining population groups, for tracing the genetic history of an individual or a particular group of related individuals, and for constructing deep-branch taxonomic phylogenies. There is every reason to believe that the mitochondrial genome will be as valuable in studies on flatworms, such as the human schistosomes, as it has been for other taxa. To date, however, our knowledge of mitochondrial genomes of flatworms remains limited, and this review summarises the currently available information. In particular, details of the recent sequence obtained for cloned Schistosoma mansoni mitochondrial DNA fragments spanning over half of the mitochondrial genome of this species are emphasised. This and other information, available as a result of the Schistosome Genome Project, provide the basis for obtaining the complete mitochondrial DNA sequence and gene order of S. mansoni and the other human schistosomes. The availability of complete mitochondrial DNA sequences from the different species will facilitate much more in-depth study of genetic diversity and host specificity in schistosomes and the interrelationships between the various forms infecting humans and between these and other flatworms.