Figures
Abstract
The whipworm, Trichuris trichiura, causes trichuriasis in ∼600 million people worldwide, mainly in developing countries. Whipworms also infect other animal hosts, including pigs (T. suis), dogs (T. vulpis) and non-human primates, and cause disease in these hosts, which is similar to trichuriasis of humans. Although Trichuris species are considered to be host specific, there has been considerable controversy, over the years, as to whether T. trichiura and T. suis are the same or distinct species. Here, we characterised the entire mitochondrial genomes of human-derived Trichuris and pig-derived Trichuris, compared them and then tested the hypothesis that the parasites from these two host species are genetically distinct in a phylogenetic analysis of the sequence data. Taken together, the findings support the proposal that T. trichiura and T. suis are separate species, consistent with previous data for nuclear ribosomal DNA. Using molecular analytical tools, employing genetic markers defined herein, future work should conduct large-scale studies to establish whether T. trichiura is found in pigs and T. suis in humans in endemic regions.
Author Summary
Trichuriasis is a neglected tropical disease (NTD) caused by parasitic nematodes of the genus Trichuris (Nematoda), causing significant human and animal health problems as well as considerable socio-economic consequences world-wide. Although Trichuris species are considered to be relatively host specific, there has been significant controversy as to whether Trichuris infecting humans (recognized as T. trichiura) is a distinct species from that found in pigs (recognized as T. suis), or not. In the present study, we sequenced, annotated and compared the complete mitochondrial genomes of Trichuris from these two hosts and undertook a phylogenetic analysis of the mitochondrial datasets. This analysis showed clear genetic distinctiveness and strong statistical support for the hypothesis that T. trichiura and T. suis are separate species, consistent with previous studies using nuclear ribosomal DNA sequence data. Future studies could explore, using mitochondrial genetic markers defined in the present study, cross-transmission of Trichuris between pigs and humans in endemic regions, and the population genetics of T. trichiura and T. suis.
Citation: Liu G-H, Gasser RB, Su A, Nejsum P, Peng L, Lin R-Q, et al. (2012) Clear Genetic Distinctiveness between Human- and Pig-Derived Trichuris Based on Analyses of Mitochondrial Datasets. PLoS Negl Trop Dis 6(2): e1539. https://doi.org/10.1371/journal.pntd.0001539
Editor: Aaron R. Jex, University of Melbourne, Australia
Received: July 14, 2010; Accepted: January 11, 2012; Published: February 21, 2012
Copyright: © 2012 Liu et al. 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 author and source are credited.
Funding: This work was supported, in part, by the Program for Outstanding Scientists in Agricultural Research, the Open Funds of the State Key Laboratory of Veterinary Etiological Biology, Lanzhou Veterinary Research Institute, Chinese Academy of Agricultural Sciences (grant nos. SKLVEB2011KFKT004, SKLVEB2011KFKT011, SKLVEB2010KFKT009, SKLVEB2009KFKT008 and SKLVEB2011KFKT010) and the Yunnan Provincial Program for Introducing High-level Scientists (grant no. 2009CI125). Work in the Gasser Lab was supported mainly by the Australian Research Council (ARC). The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.
Competing interests: The authors have declared that no competing interests exist.
Introduction
Soil-transmitted helminths ( = geohelminths), including whipworm, are responsible for neglected tropical diseases (NTDs) of humans in developing countries [1]–[3]. Trichuris trichiura infects ∼600 million people worldwide. This parasite is transmitted directly via a direct, faecal-oral route. The thick-shelled (infective) eggs are ingested and then hatch, following gastric passage, in the small intestine. First-stage larvae (L1s) are released and migrate to the large intestine (caecum and colon), where they develop, following multiple moults, into adults (∼30–50 mm in length). The worms burrow their thin, thread-like anterior end into the mucosal lining of the large intestinal wall, feed on tissue fluids, mature and produce eggs. In the large intestines, large numbers of worms cause disease ( = trichuriasis), which is usually associated with entero-typhlocolitis and clinical signs, such as dysentery, bloody diarrhoea and/or rectal prolapse, in people with a high intensity of infection. Children (∼5–15 years of age) often harbour the largest numbers of worms [2]. Whipworms also infect other animal hosts, including non-human primates, pigs and dogs, and can cause clinical disease similar to trichuriasis of humans [4]–[6].
Based on current knowledge, Trichuris species are considered to specifically infect a particular host species or a group of related hosts. Trichuris species are usually identified based on host origin and the morphological features of the adult worm (spicule and pericloacal papillae) [7], [8]. However, it is not always possible to unequivocally identify and differentiate Trichuris species based on the morphology of adult worms alone. Importantly, T. trichuria cannot be unequivocally differentiated morphologically from T. suis or Trichuris from some other animals, such as non-human primates [7]. Over the years, there has been considerable discussion as to whether T. trichuira and T. suis are the same or distinct species [9]–[13], and whether humans can become infected with T. suis, and pigs with T. trichiura in endemic countries in which both host species live in close association. Although the authors of a recent molecular study suggested that T. suis is a separate species from T. trichiura [14], only a small number of specimens from one country (Spain) was used in this study, and amplicons (from the first and second internal transcribed spacers, ITS-1 and ITS-2, of nuclear ribosomal DNA) were subjected to cloning prior to sequencing, which has significant potential to lead to artefacts [15], [16]. Therefore, the findings from this study [14] need to be interpreted with some caution at this stage. Moreover, internal transcribed spacers (ITS) of nuclear ribosomal DNA might not be suited as specific markers for enoplid nematodes, because of sequence polymorphism (heterogeneity) that occurs within species (or individuals) [14], [17].
Given this heterogeneity in nuclear rDNA, barcoding from whole mitochondrial (mt) genomes (haploid) has major advantages, particularly when concatenated protein sequences derived from all coding genes are used as markers in comparative, phylogenetic-based analyses [18]–[25]. Therefore, in the present study, we (i) characterised the mt genomes of human-derived Trichuris and pig-derived Trichuris, (ii) compared these genomes and (iii) then tested the hypothesis that human-Trichuris and pig-Trichuris are genetically distinct in a phylogenetic analysis of sequence data sets representing both genomes and those from selected nematodes for comparative purposes.
Materials and Methods
Ethics statement
This study was approved by the Animal Ethics Committee of the Lanzhou Veterinary Research Institute, Chinese Academy of Agricultural Sciences. For the collection of Trichuris from humans, the subjects provided informed, written consent. All pigs, from which Trichuris specimens were collected, were handled in accordance with good animal practices required by the Animal Ethics Procedures and Guidelines of the People's Republic of China.
Parasites and isolation of total genomic DNA
Adult specimens of Trichuris were collected from the caecum of a human patient during surgery in Zhanjiang People's Hospital in Zhanjiang, Guangdong Province, China. Adult specimens of Trichuris were also collected from the caecum from a pig slaughtered in an abattoir in Zhanjiang in the same province. Adult worms from each host were washed separately in physiological saline, identified morphologically [9], [26], fixed in 70% (v/v) ethanol and stored at −20°C until use. Total genomic DNA was isolated separately from two individual worms (coded Ttr2 and TsCS1 for human-Trichuris and pig-Trichuris, respectively) using an established method [27]. The region spanning ITS-1, the 5.8S gene and ITS-2 was amplified from each of these individuals by PCR using previously reported primers [14] and sequenced directly. The ITS-1 sequence of the human-Trichuris sample had 99.3% similarity with that of T. trichiura from human in Thailand (GenBank accession no. GQ352554). The ITS-1 and ITS-2 sequences of the pig-Trichuris sample had 98.6% and 98.5% similarity with that of T. suis from pigs in Spain (GenBank accession nos. AJ781762 and AJ249966, respectively) [14].
Long-range PCR-based sequencing of mt DNA
To obtain some mt gene sequence data for primer design, we amplified regions (400–500 bp) of the cox1 and nad1 genes by using (relatively) conserved primers JB3/JB4.5 and JB11/JB12, respectively [28], and of nad4 and rrnL genes using primers designed in this study (Table 1) by PCR. The amplicons were sequenced from both directions, using BigDye terminator v3.1, ABI PRISM 3730. We then designed primers (see Table 2) to regions within cox1, nad1, nad4 and rrnL and amplified from total genomic DNA (from an individual worm) the entire mt genome in three (for human-Trichuris) or four (for pig-Trichuris) overlapping fragments (of ∼2–4 kb each) between nad1 and nad4, nad4 and rrnL, and rrnL and cox1, cox1 and nad1. The cycling conditions used were 92°C for 2 min (initial denaturation), then 92°C for 10 s (denaturation), 50°C for 30 s (annealing), and 60–68°C for 10 min (extension) for 10 cycles, followed by 92°C for 10 s, 50°C for 30 s, and 60–68°C for 10 min for 20 cycles, with a cycle elongation of 10 s for each cycle and a final extension at 60–68°C for 7 min. Each amplicon, which represented a single band in a 0.8% (w/v) agarose gel, following electrophoresis and ethidium-bromide staining, was column-purified and then sequenced using a primer walking strategy [29].
Sequence analyses
Sequences were assembled manually and aligned against the complete mt genome sequences of other nematodes (available publicly) using the computer program Clustal X 1.83 [30] to infer gene boundaries. The open-reading frames (ORFs) and codon usages of protein-coding genes were predicted using the program MacVector v.4.1.4 (Kodak), and subsequently compared with that of Trichinella spiralis [31]. Translation initiation and translation termination codons were identified based on comparison with those reported previously [31]. Codon usages were examined based on the relationships between the nucleotide composition of codon families and amino acid occurrence, for which codons are partitioned into AT rich codons, GC-rich codons and unbiased codons. The secondary structures of 22 tRNA genes were predicted using tRNAscan-SE [32] and/or manual adjustment [33].
Phylogenetic analyses
Amino acid sequences inferred from the 12 protein-coding genes (i.e. not atp-8) common among all of the nematodes included here were concatenated into a single alignment, and then aligned with those of 9 other enoplid nematodes (Trichinella spiralis, GenBank accession number NC_002681; Xiphinema americanum, NC_005928; Hexamermis agrotis, NC_008828; Agamermis sp., NC_008231; Romanomermis culicivorax, NC_008640; Romanomermis iyengari, NC_008693; Romanomermis nielseni, NC_008692; Strelkovimermis spiculatus, NC_008047; Thaumamermis cosgrovei, NC_008046), using the chromadorean nematode, Brugia malayi (NC_004298) as the outgroup. Any regions of ambiguous alignment were excluded using Gblocks (http://molevol.cmima.csic.es/castresana/Gblocks_server.html; Talavera and Castresana 2007) using stringent selection criteria (do not allow many contiguous nonconserved positions). Phylogenetic analyses were conducted using three methods: Bayesian inference (BI) analysis was conducted with four independent Markov chain runs for 1,000,000 metropolis-coupled MCMC generations, sampling a tree every 100th generation in MrBayes 3.1.1 [34]; the first 2,500 trees represented burn-in, and the remaining trees were used to calculate Bayesian posterior probabilities (pp). Maximum likelihood (ML) analyses were performed using PhyML 3.0 [35], and the (mtREV for amino acid sequences and GTR for rrnL nucleotide sequences) models were determined based on the Akaike information criterion (AIC). Bootstrap support was calculated using 100 bootstrap replicates. Maximum parsimony (MP) analysis was conducted using PAUP 4.0 Beta 10 [36], with indels treated as missing character states; 1,000 random additional searches were performed using TBR. Bootstrap support was calculated using 1,000 bootstrap replicates, and 100 random taxon additions in PAUP. Phylograms were drawn using the program Tree View v.1.65 [37].
Sequencing of rrnL and analysis
Two primers, rrnLF (5′-TAAATGGCCGTCGTAACGTGACTGT-3′) and rrnLR (5′- AAAGAGAATCCATTCTATCTCGCAACG-3′), were employed for PCR amplification and subsequent sequencing of a portion (471 bp for human-Trichuris and 482 bp for pig-Trichuris) of the large subunit of mt ribosomal RNA (rrnL) from multiple individuals of human- and pig-derived Trichuris (Table 3). The rrnL sequence from T. spiralis (accession number NC_002681) [31] was used as the outgroup in phylogenetic analyses, because this morphologically distinct species is related to Trichuris [38]. All rrnL sequences were aligned using Clustal X, and the alignment was modified manually, based on the predicted secondary structure of the rrnL for Trichuris [31], and then subjected to phylogenetic analysis using the same methods as described above.
Results
Features of the mt genomes of Trichuris from the human or pig host
The complete mt genome sequences were 14,046 nt (human-Trichuris) and 14,436 nt (pig-Trichuris) in length, respectively (GenBank accession numbers GU385218 and GU070737). Each mt genome contains 13 protein-coding genes (cox1-3, nad1-6, nad4L, cytb, atp6 and atp8), 22 transfer RNA genes and two ribosomal RNA genes (rrnS and rrnL) (Table 4). The genes nad6 and cox3 overlapped by 25 bp (human-Trichuris); nad5 overlaps by 7 and 1 bp with tRNA-His, tRNA-Ser(AGN) by 8 and 3 bp with rrnS, tRNA-Asp by 13 and 4 bp with atp8 for human-Trichuris and pig-Trichuris, respectively. The atp8 gene is encoded (Fig. 1), as is typical for adenophorean nematodes [31]. The protein-coding genes are transcribed in different directions, as described for T. spiralis and X. americanum [31], [39]. Except for four protein-coding genes (nad2, nad5, nad4 and nad4L) and six tRNA genes (tRNA-Arg, tRNA-His, tRNA-Met, tRNA-Phe, tRNA-Pro and tRNA-Thr) encoded on the L-strand, all other genes were encoded on the H-strand. The AT-rich regions are located between nad1 and tRNA-Lys, and nad3 and tRNA-Ser(UCN), which differs from those of secernentean nematodes [19], [33]. The nucleotide composition of the entire mt genome is biased toward A and T, with T being the most favoured nucleotide and G being least favoured, which is consistent with mt genomes of some other nematodes for which mt genomic data are available [18], [19], [21], [31], [33]. The overall A+T content is 68.1% for human-Trichuris (33.6% A, 34.5% T; 15.0% G and 16.9% C) and 71.5% for pig-Trichuris (35.6% A, 35.9% T; 13.4% G and 15.1% C).
Genes follow standard nomenclature [18], except for the 22 tRNA genes, which are designated using one-letter amino acid codes, with numerals differentiating each of the two leucine- and serine-specifying tRNAs (L1 and L2 for codon families CUN and UUR, respectively; S1 and S2 for codon families AGN and UCN, respectively). “NCR-L” refers to a large non-coding region; “NCR-S” refers to a small non-coding region.
Annotation
Protein-coding genes were annotated by aligning sequences and identifying translation initiation and termination codons by comparison with inference sequences for other nematodes. For both human-Trichuris and pig-Trichuris, the lengths of protein-coding genes were in the following order: nad5 (1548–1557 bp) >cox1>nad4>cytb>nad1>nad2>atp6>cox3>cox2>nad6>nad3>nad4L>atp8 (165–171 bp) (Table 4). The longest gene is nad5, and the lengths of the nad1 and nad3 genes are the same for human-Trichuris and pig-Trichuris (Table 5). The inferred nucleotide and amino acid sequences of each of the 13 mt proteins of human-Trichuris were compared with pig-Trichuris. For individual genes, the nucleotide and amino acid sequence differences between human-Trichuris and pig-Trichuris vary from 25.4 to 37.4% and 13.6 to 62.5%, respectively (Table 5).
A total of 3559 and 3562 amino acids are encoded in the mt genome of human-Trichuris and pig-Trichuris, respectively. As the mt genomes of nematodes can contain non-standard initiation codons [18], the identification of initiation codons can sometimes be challenging. For human-Trichuris, five genes (cox1, cox2, cox3, cytb and nad4) start with ATG and eight genes (nad1, nad2, nad3, nad5, nad6, nad4L, atp6 and atp8) use ATA. All genes have complete termination codons, 11 genes (cox1, cox2, cox3, nad1, nad2, nad3, nad4, nad4L, nad5, nad6 and atp6) use TAA and two genes (atp8 and cytb) use TAG as a termination codon, respectively. For pig-Trichuris, six genes (nad2, nad3, nad4 nad4L nad5 and atp6) start with ATA, four genes (cox1, cox2, cox3 and cytb) with ATG, and one gene (atp8) with TTG. All protein-coding genes have complete termination codon; seven genes (cox2, cox3, nad3, nad4, nad6, atp6 and atp8) stop with TAA and six genes (cox1, nad1, nad2, nad5, nad4L and cytb) use TAG. No abbreviated stop codons, such as TA or T, were detected. Such codons are known to occur in the mt genomes of other nematodes, such as Strongyloides stercoralis (cytb, nad4, nad1 and atp6) and T. spiralis (cytb and nad4) and Caenorhabditis elegans (nad1 and nad3) [19], [31], [40]. The atp8 gene was inferred by comparison with the homologous gene of T. spiralis. The ORF of this gene was located between genes tRNA-Asp and nad3 and inferred by the presence of a methionine for human-Trichuris and a leucine for pig-Trichuris.
Twenty-two tRNA genes were predicted from the mt genomes of human-Trichuris and pig-Trichuris and varied from 50 to 67 nt in length. Most of the tRNA genes are smaller than the corresponding genes in the mt genomes of other nematodes due to a reduced TΨC stem-loop region (TV-replacement loop) or DHU stem-loop region [39]. Most of the tRNA gene sequences can be folded into conventional secondary four-arm cloverleaf structures. In these tRNA, there is a strict conservation of the sizes of the amino acid acceptor stem (11–15 bp) and the anticodon loop (7 bp). Their D-loops consist of 5–9 bp. The two tRNA-Ser each contain the TΨC arm and loop, but lack the DHU arm and loop.
The two ribosomal RNA genes (rrnL and rrnS) of human-Trichuris and pig-Trichuris were inferred based on comparisons with sequences from T. spiralis; rrnL is located between tRNA-Val and atp6, and rrnS is located between tRNA-Ser(AGN) and tRNA-Val. The length of rrnL is 1011 bp for both human-Trichuris and pig-Trichuris. The lengths of the rrnS genes are 698 bp for human-Trichuris and 712 bp for pig-Trichuris. The A+T contents of rrnL for human-Trichuris and pig-Trichuris are 72.5% and 76.4%, respectively. The A+T contents of rrnS for human-Trichuris and pig-Trichuris are 69.9% and 75.4%, respectively.
Two AT-rich non-coding regions (NCRs) were inferred in the mt genomes of both human-Trichuris and pig-Trichuris. For these genomes, the long NCR (designated NCR-L; 162 bp and 144 bp in length, respectively) is located between the nad1 and tRNA-Lys (Fig. 1), has an A+T content of 71–72%. This overall A+T content is lower than those reported for nematodes (77.9–93.1%) studied to date [18], [19], [21], [33]. In this NCR, there are also 26 nt (human-Trichuris) and 17 nt (pig-Trichuris) AT dinucleotide repeats. Similar repeats have been detected in this region in C. elegans and A. suum [40]. For both human-Trichuris and pig-Trichuris, the short NCR (NCR-S; 93 bp and 117 bp in length) is located between genes nad3 and tRNA-Ser (UCN) (Fig. 1), with an A+T content of 65.6% and 84.6%, respectively. This region contains dinucleotide [AT]26 repeats and might form a hairpin loop structure (cf. AAAAAAAATTTTTTTTTT). Although nothing is yet known about the replication process in the mt DNA of parasitic nematodes, the high A+T content and the predicted structure of the AT-rich NCRs suggest an involvement in the initiation of replication [41].
Comparative analyses between human-Trichuris and pig-Trichuris.
The full mt genome sequence of human-Trichuris (accession no. GU385218) was 14046 bp in length, 390 bp shorter than that of pig-Trichuris (accession no. GU070737). The arrangement of the mt genes (i.e., 13 protein genes, 2 rrn genes and 22 tRNA genes) and NCRs were the same. A comparison of the nucleotide sequences of each mt gene and NCR, as well as the amino acid sequences, conceptually translated from all protein genes of the two Trichuris, is given in Tables 4 and 5. The sequence lengths of individual genes and NCRs were the same for human-Trichuris and pig-Trichuris, except for variation of three to nine nucleotides in the each of 13 coding gene, 18 nucleotides in the first non-coding region and of one to nine nucleotides in each of 22 tRNA genes (Table 4). The magnitude of sequence variation in each gene and NCR between human-Trichuris and pig-Trichuris ranged from 24.6–47.4%. Sequence difference across the entire mt genome was 32.87%. The greatest variation was in the atp8 gene (47.4%), whereas least differences (24.6% and 25.1%) were detected in the rrnS and rrnL subunits, respectively (Table 5).
Amino acid sequences inferred from individual mt protein genes of human-Trichuris were compared with those of pig-Trichuris. The amino acid sequence differences ranged from 13.6–62.5%, with COX1 being the most conserved protein, and ATP8 the least conserved. There were 1299 amino acid substitutions (for an alignment length of 3559 and 3562 positions) in the 13 proteins, the majority of which were in the proteins NAD4 (n = 235), NAD5 (n = 220) and NAD2 (n = 120). Of the 1299 amino acid substitutions, 407 (31.3%) represented potentially informative characters, the greatest number of them being in the proteins NAD1 (n = 84), NAD5 (n = 55) and ATP6 (n = 37). The phylogenetic analyses of amino acid sequence datasets using B. malayi as the outgroup reflected the clear genetic distinctiveness between human-Trichuris and pig-Trichuris and also the grouping of these two members of Trichuris with T. spiralis (Trichocephalida), with absolute support, to the exclusion of members of the Dorylaimida and Mermithida (Fig. 2).
The phylogenetic relationship was inferred based on the concatenated amino acid sequences of 12 protein-coding genes (with the exception of atp8) by Bayesian (BI), maximum likelihood (ML) and maximum parsimony (MP) analyses, using Brugia malayi (a chromadorean nematode) (NC_004298) as the outgroup. The numbers along branches indicate bootstrap values resulting from different analyses in the order: BI/ML/MP.
Comparison of the mt genomes of human-Trichuris and pig-Trichuris showed that the rrnS and rrnL were the two most conserved genes (Table 5). Sequence variation in part of the rrnL gene was assessed among 16 individuals of Trichuris from humans and pigs (Table 3). Sequences of the 6 human-Trichuris individuals were of the same length (419 bp). Nucleotide variation among the 6 human-Trichuris individuals was detected at 4 sites (sequence positions 206, 228, 233 and 274; GenBank accession numbers AM993017–AM993022). Sequences of the 10 pig-Trichuris individuals were of the same length (430 bp). Nucleotide variation also occurred at 4 sites (sequence positions 184, 233, 318 and 395, GenBank accession numbers AM993023–AM993032).
The alignment of the partial rrnL sequences revealed that all individuals of human-Trichuris differed at 89 nucleotide positions when compared with pig-Trichuris. These differences included 15 indels, 16 purine transitions (A<->G) and 19 transversion (A<->T). Phylogenetic analyses of the rrnL sequence data revealed strong support for the separation of human-Trichuris from pig-Trichuris individuals into two distinct clades, and the trees produced using the three different methods were essentially the same in topology (Fig. 3). Sixteen of the 89 nucleotide differences were considered as derived (i.e., autapomorphic) characters, using T. spiralis as the outgroup.
The phylogenetic relationship was inferred by Bayesian (BI), maximum likelihood (ML) and maximum parsimony (MP) analyses of mitochondrial rrnL sequence data, using Trichinella spirialis as the outgroup. The numbers along branches indicate bootstrap values from different analyses in the order: BI/ML/MP.
Discussion
A substantial level of nucleotide difference (32.9%) was detected in the complete mt genome between an individual of human-Trichuris and pig-Trichuris from China. The sequence variation detected in the 13 protein-coding genes (25.4–47.4%) and in NCRs (36.6%) was consistent with previous findings of variation in the nucleotide sequences of the nuclear ITS rDNA from human and pig [14]. However, for many nematodes [42], [43], there is usually greater within-species variation in mt protein-coding genes than in the ITS. For example, the magnitude of the nucleotide sequence variation in the 12 common mt protein genes (3–7%) [20] was greater than the 15 (1.8%) variable positions in the ITS (over 852 bp) detected among multiple individuals of the human hookworm, N. americanus [44].
Comparison between human- and pig-derived Trichuris from China also revealed variation at 1299 amino acid positions in the 13 predicted mt protein sequences. This level of amino acid variation (36.4%) is very high, given that mt proteins are usually conserved within a species due to structural and functional constraints [45]. In addition, previous studies of other nematodes have detected little to no within-species variation in protein sequences. For example, no within-species variation was detected in a COX1 region of 131 amino acids for N. americanus and for related hookworms, including A. caninum Ancylostoma and A. duodenale [24], [33]. Similarly, amino acid substitutions were recorded at only two of 196 (1%) positions (based on a comparison of conceptually translated sequences originating from GenBank accession nos. AF303135-AF303159) in partial COX1 among 151 N. americanus samples from four locations in China [46]. In the present study, the greatest numbers of amino acid differences between human-Trichuris and pig-Trichuris were in the NAD4 (n = 235; 40.9%), NAD5 (n = 220; 35.9%) and NAD2 (n = 120; 34.1%) sequences; these percentages were significantly higher than that (4.9–10%) between the two hookworms A. caninum and A. duodenale [24], [33]. The nature, extent and significance of the amino acid sequence variation between Trichuris from the human and pig hosts and from different geographical origins needs to be evaluated further, because there is virtually no published data on the magnitude of within-species variation in mt protein sequences for members of the genus Trichuris.
Genetic variation between human- and pig-Trichuris was also detected here in the two mt ribosomal RNA gene subunits (rrnL and rrnS). These subunits are usually more conserved in sequence than the protein genes [45], which is also supported by the present data. Comparison of the complete mt genomic data set between the two Trichuris individuals (Ttr2 and TsCS1) displayed less sequence variation in rrnS and rrnL (24.6% and 25.1%) compared with most protein genes (25.4–47.4%) and the non-coding regions (36.6%) (Table 5). A region (∼430 bp) in the conserved rrnL gene was used to examine the magnitude of genetic variation in Trichuris between the two different host species. A comparison of the partial rrnL sequences among 16 Trichuris individuals revealed 89 (20.7%) variable positions between human-Trichuris and pig-Trichuris, which is comparable with previous findings of a significant genetic difference (17%) in nuclear ITS between the two operational taxonomic units (OTUs) in Spain [14]. Taken together, the molecular evidence presented here supports the hypothesis that the gene pools of human-Trichuris and pig-Trichuris have been isolated for a substantial period of time and that they represent distinct species. In spite of the genetic distinctiveness recorded here between them, host affiliation is not strict [47]. Cross-infection of Trichuris between humans and pigs (both directions) has been described, but infection in the heterologous host is usually abbreviated [47].
In spite of the compelling evidence of genetic distinctiveness between Trichuris specimens from human and pig hosts, interpretation from this study needs to be somewhat guarded until detailed population genetic investigations have been conducted. Future studies could (i) explore, in detail, nucleotide variation in ribosomal and mt DNAs within and among Trichuris populations from humans and pigs from a range of different countries employing, for example, mutation scanning-coupled sequencing [27], (ii) establish, using accurate molecular tools, whether there is a particular affiliation between Trichuris and host in endemic regions and whether cross-host species infection is common or not, and (iii) attempt to establish an experimental infection of Trichuris of human origin in pigs, in order to be able to investigate the genetic and reproductive relationships between human-Trichuris and pig-Trichuris. Moreover, given the advent of high throughput genomic sequencing technologies, and the recent success in sequencing the nuclear genomes of the parasitic nematodes, B. malayi [48] and Ascaris suum [49], it is conceivable that the genomes of human-Trichuris and pig-Trichuris will be characterized in the near future. The transcriptome, and inferred proteome, characterised recently [50] will assist in future efforts to decode these genomes. Such work will pave the way for future fundamental molecular explorations and the design of new methods for the treatment and control of one of the world's socio-economically important nematodes [3]. This focus is important, given the impact of Trichuris and other soil-transmitted helminths (STHs), which affect billions of people and animals world-wide. Although Trichuris species are seriously neglected, genomics and related approaches provide new opportunities for the discovery of novel intervention strategies, with major implications for improving animal and human health and well being globally. In addition, the implications of genomic studies could also be highly relevant in relation to finding new treatments for immune-pathological diseases of humans [50]. Interestingly, various studies [51]–[55] have indicated that iatrogenic infections of human patients suffering from immunological disorders (such as inflammatory bowel disease, IBD) with nematodes, such as pig-Trichuris eggs can significantly suppress clinical symptoms. Although the mechanisms by which Trichuris modulates the human immune system are still unclear [52], [56], [57], studies have proposed that a modified CD4+ T helper 2 (Th2)-immune response and the production of anti-inflammatory cytokines, such the interleukins (IL-) IL-4 and IL-10, contribute to the inhibition of effector mechanisms [56], [58], [59]. Therefore, detailed investigations of pig-Trichuris at the molecular level could provide enormous scope for studying immuno-molecular mechanisms that take place between the parasite and humans affected by autoimmune or other immune diseases. The mt genetic markers defined in the present study should be useful to verify the specific identity of Trichuris employed in such studies.
Author Contributions
Conceived and designed the experiments: RBG XQZ GHL. Performed the experiments: GHL AS. Analyzed the data: GHL PN MWL RBG XQZ. Contributed reagents/materials/analysis tools: LP RQL MJX. Wrote the paper: RBG GHL XQZ.
References
- 1. Hotez PJ, Molyneux DH, Fenwick A, Kumaresan J, Sachs SE, et al. (2007) Control of neglected tropical diseases. N Engl J Med 357: 1018–1027.
- 2. Hotez PJ (2009) The neglected tropical diseases and their devastating health and economic impact on the member nations of the Organisation of the Islamic Conference. PLoS Negl Trop Dis 3: e539.
- 3. Bethony JM, Cole RN, Guo X, Kamhawi S, Lightowlers MW, et al. (2011) Vaccines to combat the neglected tropical diseases. Immunol Rev 239: 237–270.
- 4. Beer RJ, Lean IJ (1973) Clinical trichuriasis produced experimentally in growing pigs. I. Pathology of infection. Vet Rec 93: 189–195.
- 5. Hall GA, Rutter JM, Beer RJ (1976) A comparative study of the histopathology of the large intestine of conventionally reared, specific pathogen free and gnotobiotic pigs infected with Trichuris suis. J Comp Pathol 86: 285–292.
- 6. Stewart TB, Hale OM (1988) Losses to internal parasites in swine production. J Anim Sci 166: 1548–1554.
- 7. Knight RA (1971) Redescription of Trichuris discolor (von Linstow, 1906) and T. skrjabini (Baskakov, 1924) from domestic ruminants in the United States and comparisons with T. ovis (Abildgaard, 1795). J Parasitol 57: 302–310.
- 8. Spakulová M (1994) Discriminant analysis as a method for the numerical evaluation of taxonomic characters in male trichurid nematodes. Syst Parasitol 29: 113–119.
- 9. Ooi HK, Tenora F, Itoh K, Kamiya M (1993) Comparative study of Trichuris trichiura from non-human primates and from man, and their difference with T. suis. J Vet Med Sci 55: 363–366.
- 10. Chandler AC (1930) Specific characters in the genus Trichuris with a description of the new species Trichuris tenuis from a camel. J Parasitol 16: 198–206.
- 11. Barus V, Kotrla B, Tenora F (1977) A scanning electron microscopic study of spicular sheath of some trichurids (Nematoda). Folia Parasitol 24: 107–110.
- 12. Mahmoud LH (2002) Scanning electron microscopy of Trichuris trichiura. J Egypt Soc Parasitol 32: 469–474.
- 13. Zaman V (1984) Scanning electron microscopy of Trichuris trichiura (Nematoda). Acta Trop 41: 287–292.
- 14. Cutillas C, Callejon R, Rojas MD, Tewes B, Ubeda JM (2009) Trichuris suis and Trichuris trichiura are different nematode species. Acta Trop 111: 299–307.
- 15. Gasser RB (1997) Mutation scanning methods for the analysis of parasite genes. Int J Parasitol 27: 1449–1463.
- 16. Gasser RB (2006) Molecular technologies in parasitology, with an emphasis on genomic approaches for investigating parasitic nematodes. Parassitologia 48: 9–11.
- 17. Cutillas C, de Rojas M, Ariza C, Ubeda JM, Guevara D (2007) Molecular identification of Trichuris vulpis and Trichuris suis isolated from different hosts. Parasitol Res 100: 383–389.
- 18. Li MW, Lin RQ, Song HQ, Wu XY, Zhu XQ (2008) The complete mitochondrial genomes for three Toxocara species of human and animal health significance. BMC Genomics 9: e224.
- 19. Hu M, Chilton NB, Gasser RB (2003) The mitochondrial genome of Strongyloides stercoralis (Nematoda)-idiosyncratic gene order and evolutionary implications. Int J Parasitol 33: 1393–1408.
- 20. Hu M, Chilton NB, Abs El-Osta YG, Gasser RB (2003) Comparative analysis of mitochondrial genome data for Necator americanus from two endemic regions reveals substantial genetic variation. Int J Parasitol 33: 955–963.
- 21. Kim KH, Eom KS, Park JK (2006) The complete mitochondrial genome of Anisakis simplex (Ascaridida: Nematoda) and phylogenetic implications. Int J Parasitol 36: 319–328.
- 22. Lin RQ, Qiu LL, Liu GH, Wu XY, Weng YB, et al. (2011) Characterization of the complete mitochondrial genomes of five Eimeria species from domestic chickens. Gene 480: 28–33.
- 23. Liu GH, Lin RQ, Li MW, Liu W, Liu Y, et al. (2011) The complete mitochondrial genomes of three cestode species of Taenia infecting animals and humans. Mol Biol Rep 38: 2249–2256.
- 24. Jex AR, Waeschenbach A, Hu M, van Wyk JA, Beveridge I, et al. (2009) The mitochondrial genomes of Ancylostoma caninum and Bunostomum phlebotomum-two hookworms of animal health and zoonotic importance. BMC Genomics 10: e79.
- 25. Jex AR, Hall RS, Littlewood DT, Gasser RB (2010) An integrated pipeline for next-generation sequencing and annotation of mitochondrial genomes. Nucleic Acids Res 38: 522–533.
- 26. Jenkins T (1970) A morphological and histochemical study of Trichuris suis (Schrank, 1788) with special reference to the host-parasite relationship. Parasitology 61: 357–374.
- 27. Gasser RB, Hu M, Chilton NB, Campbell BE, Jex AJ, et al. (2006) Single-strand conformation polymorphism (SSCP) for the analysis of genetic variation. Nat Protoc 1: 3121–3128.
- 28. Gasser RB, Zhu XQ, McManus DP (1999) NADH dehydrogenase subunit 1 and cytochrome c oxidase subunit I sequences compared for members of the genus Taenia (Cestoda). Int J Parasitol 29: 1965–1970.
- 29. Hu M, Jex AR, Campbell BE, Gasser RB (2007) Long PCR amplification of the entire mitochondrial genome from individual helminths for direct sequencing. Nature Protoc 2: 2339–2344.
- 30. Thompson JD, Gibson TJ, Plewniak F, Jeanmougin F, Higgins DG (1997) The Clustal X windows interface: flexible strategies for multiple sequence alignment aided by quality analysis tools. Nucleic Acids Res 24: 4876–4882.
- 31. Lavrov DV, Brown WM (2001) 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 157: 621–637.
- 32. Lowe TM, Eddy SR (1997) tRNAscan-SE: A program for improved detection of transfer RNA genes in genomic sequence. Nucleic Acids Res 25: 955–964.
- 33. Hu M, Chilton NB, Gasser RB (2002) The mitochondrial genomes of the human hookworms, Ancylostoma duodenale and Necator americanus (Nematoda: Secernentea). Int J Parasitol 32: 145–158.
- 34. Ronquist F, Huelsenbeck JP (2003) MrBayes 3: Bayesian phylogenetic inference under mixed models. Bioinformatics 19: 1572–1574.
- 35. Guindon S, Gascuel O (2003) A simple, fast, and accurate algorithm to estimate large phylogenies by maximum likelihood. Syst Biol 52: 696–704.
- 36.
Swofford DL (2002) Paup*: Phylogenetic Analysis Using Parsimony, version 4.0b10. Sinauer Associates, Sunderland, MA.
- 37. Page RD (1996) TREEVIEW: an application to display phylogenetic trees on personal computers. Comput Appl Biosci 12: 357–358.
- 38. Blaxter ML, De Ley P, Garey JR, Liu LX, Scheldeman P, et al. (1998) A molecular evolutionary framework for the phylum Nematoda. Nature 392: 71–75.
- 39. He Y, Jones J, Armstrong M, Lamberti F, Moens M (2005) The mitochondrial genome of Xiphinema americanum sensu stricto (Nematoda: Enoplea): considerable economization in the length and structural features of encoded genes. J Mol Evol 61: 819–833.
- 40. Okimoto R, Macfarlane L, Clary D, Wolstenholme DR (1992) The mitochondrial genomes of two nematodes, Caenorhabditis elegans and Ascaris suum. Genetics 130: 471–498.
- 41. Boore JL (1999) Animal mitochondrial genomes. Nucleic Acids Res 27: 1767–1780.
- 42. Blouin MS (2002) Molecular prospecting for cryptic species of nematodes: mitochondrial DNA versus internal transcribed spacer. Int J Parasitol 32: 527–531.
- 43. Gasser RB (2006) Molecular tools-advances, opportunities and prospects. Vet Parasitol 136: 69–89.
- 44. Romstad A, Gasser RB, Nansen P, Polderman AM, Chilton NB (1998) Necator americanus (Nematoda: Ancylostomatidae) from Africa and Malaysia have different ITS-2 rDNA sequences. Int J Parasitol 28: 611–615.
- 45. Hu M, Chilton NB, Gasser RB (2004) The mitochondrial genomics of parasitic nematodes of socio-economic importance: recent progress, and implications for population genetics and systematics. Adv Parasitol 56: 133–212.
- 46. Hawdon JM, Li T, Zhan B, Blouin MS (2001) Genetic structure of populations of the human hookworm, Necator americanus, in China. Mol Ecol 10: 1433–1437.
- 47. Beer RJ (1976) Relationship between Trichuris trichiura (Linnaeus 1758) of man and Trichuris suis (Schrank 1788) of pig. Res Vet Sci 20: 47–54.
- 48. Ghedin E, Wang S, Spiro D, Caler E, Zhao Q, et al. (2007) Draft genome of the filarial nematode parasite Brugia malayi. Science 317: 1756–1760.
- 49. Jex AR, Liu S, Li B, Young ND, Hall RS, et al. (2011) Ascaris suum draft genome. Nature 479: 529–533.
- 50. Cantacessi C, Young ND, Nejsum P, Jex AR, Campbell BE, et al. (2011) The Transcriptome of Trichuris suis - First molecular insights into a parasite with curative properties for key immune diseases of humans. PLoS One 6: e23590.
- 51. Summers RW, Elliott DE, Weinstock JV (2005) Is there a role for helminths in the therapy of inflammatory bowel disease? Nat Clin Pract Gastroenterol Hepatol 2: 62–63.
- 52. Summers RW, Elliott DE, Qadir K, Urban JF Jr, Thompson R, et al. (2005) Trichuris suis therapy in Crohn's disease. Gut 54: 87–90.
- 53. Summers RW, Elliott DE, Weinstock JV (2005) Why Trichuris suis should prove safe for use in inflammatory bowel disease. Inflamm Bowel Dis 11: 783–784.
- 54. Bager P, Arnved J, Ronborg S, Wohlfahrt J, Poulsen LK, et al. (2010) Trichuris suis ova therapy for allergic rhinitis: a randomized, double-blind, placebocontrolled clinical trial. J Allergy Clin Immunol 125: 123–130.
- 55. Hepworth MR, Hamelmann E, Lucius R, Hartmann S (2010) Looking into the future of Trichuris suis therapy. J Allergy Clin Immunol 125: 767–768.
- 56. Erb KJ (2008) Can helminths or helminth-derived products be used in humans to prevent or treat allergic diseases? Trends Immunol 30: 75–82.
- 57. Reddy A, Fried B (2007) The use of Trichuris suis and other helminth therapies to treat Crohn's disease. Parasitol Res 100: 921–927.
- 58. Hunter MM, McKay DM (2004) Helminths as therapeutic agents for inflammatory bowel disease. Aliment Pharmacol Ther 19: 167–177.
- 59. Figueiredo CA, Barreto ML, Rodrigues LC, Cooper PJ, Silva NB, et al. (2010) Chronic intestinal helminth infections are associated with immune hyporesponsiveness and induction of a regulatory network. Infect Immun 78: 3160–3167.