- Open Access
Whipworms in humans and pigs: origins and demography
Parasites & Vectors volume 9, Article number: 37 (2016)
Trichuris suis and T. trichiura are two different whipworm species that infect pigs and humans, respectively. T. suis is found in pigs worldwide while T. trichiura is responsible for nearly 460 million infections in people, mainly in areas of poor sanitation in tropical and subtropical areas. The evolutionary relationship and the historical factors responsible for this worldwide distribution are poorly understood. In this study, we aimed to reconstruct the demographic history of Trichuris in humans and pigs, the evolutionary origin of Trichuris in these hosts and factors responsible for parasite dispersal globally.
Parts of the mitochondrial nad1 and rrnL genes were sequenced followed by population genetic and phylogenetic analyses. Populations of Trichuris examined were recovered from humans (n = 31), pigs (n = 58) and non-human primates (n = 49) in different countries on different continents, namely Denmark, USA, Uganda, Ecuador, China and St. Kitts (Caribbean). Additional sequences available from GenBank were incorporated into the analyses.
We found no differentiation between human-derived Trichuris in Uganda and the majority of the Trichuris samples from non-human primates suggesting a common African origin of the parasite, which then was transmitted to Asia and further to South America. On the other hand, there was no differentiation between pig-derived Trichuris from Europe and the New World suggesting dispersal relates to human activities by transporting pigs and their parasites through colonisation and trade. Evidence for recent pig transport from China to Ecuador and from Europe to Uganda was also observed from their parasites. In contrast, there was high genetic differentiation between the pig Trichuris in Denmark and China in concordance with the host genetics.
We found evidence for an African origin of T. trichiura which were then transmitted with human ancestors to Asia and further to South America. A host shift to pigs may have occurred in Asia from where T. suis seems to have been transmitted globally by a combination of natural host dispersal and anthropogenic factors.
A range of mammalian hosts is infected with parasitic whipworms belonging to the genus Trichuris. Around 460 million humans are infected with T. trichiura mainly in developing countries in South and South-East Asia, Sub-Saharan Africa and Latin America . The prevalence of whipworms in non-human primates is generally high and despite the fact that the taxonomic status is unsettled, they are historically designated as T. trichuris . T. suis infection in pigs is globally widespread with the highest prevalence in young pigs reared in outdoor production systems [3, 4]. A number of studies have investigated the genetic relationship among pig and primate derived Trichuris but their evolutionary relationship is still poorly understood and in particular the anthropogenic and environmental factors responsible for their global distribution.
Population genetic tools provide a valuable opportunity to investigate the epidemiological history and transmission of parasites and have been adopted to study many parasitic nematodes [5, 6]. For instance, population genetic approaches were applied recently to investigate the pattern of transmission of Ascaris suum and A. lumbricoides between pigs and humans across the globe . Moreover, population genetics can be used to reconstruct the epidemiological and demographic history of micro- and macro- parasites through coalescent analysis coupled with Bayesian approaches . Reconstruction of the epidemiological and demographic history gives us a window into the past to see which factors facilitated spread or the introduction of the parasites to new regions . For instance, the demographic history of Wuchereriabancrofti suggests that this parasite was introduced to India 60,000-70,000 years ago through human migration out of Africa. In contrast, the introduction of W. bancrofti to Papua New Guinea cannot be explained by ancient human migration but must have been introduced through more recent human migration .
During the evolution of genus Homo nearly 4 million years ago, there has been continuous contact with many parasites. Those infecting humans today may have different evolutionary origin and can be categorized into two groups: 1) Parasites transmitted to humans through primate common ancestors and are referred to as “heirloom” and 2) parasites acquired more recently through contact with animals (e.g. during animal domestication in the Neolithic period roughly 10,000 years ago) and are referred to as “souvenirs” . The human whipworm is generally considered heirloom as it is found in the African non-human primates, and parasite eggs were found in human coprolites in archaeological sites before animal domestication and in the New World before the Columbian colonization [10–12]. However, there is no rigorous genetic evidence for this assumption especially with recent studies suggesting that several Trichuris species can be found in non-human primates such as the Trichuris sp. in Francois leaf monkey in China and Trichuris spp. infecting a number of African non-human primates in Uganda [13, 14]. Hence, the genetic and evolutionary relationship between Trichuris from humans and non-human primates is poorly understood.
Different molecular markers have been used to resolve the genetic and evolutionary relationship between Trichuris spp. isolates from different host speices. For instance, several studies have used the nuclear ribosomal internal transcribed spacers (ITS) to investigate the genetic relationship among Trichuris from humans, non-human primates and pigs [13, 15, 16] whereas other studies relied on small subunit ribosomal DNA (18S) and the mitochondrial cox1 gene [17, 18]. Full mitochondrial genome analysis has been used in some studies to investigate the phylogenetic relationship of Trichuris spp. but includes only a limited number of worm samples [14, 19, 20].
Herein, we investigated the genetic and evolutionary relationships between populations of Trichuris parasites derived from humans, non-human primates and pigs from different continents. The mitochondrial nad1 and rrnL genes of 138 worms were sequenced and coalescent simulations applied to infer the demographic history of the different human and pig Trichuris populations and their evolutionary origin.
The human Trichuris from Uganda and Ecuador were recovered from the faeces of children and adults after anthelmintic treatment. Permission was obtained from the Ministry of Health and the National Council of Science and Technology in Uganda and by the Ethical Committees of Liverpool School of Tropical Medicine and Pontificia Universidad Catolica del Ecuador in Ecuador. The Danish Central Medical Ethics Committee approved the study. Informed written consent was obtained from all participants or their guardians (for children) after being informed about the study in both English/Spanish and the local language.
Worms from baboons at Southwest National Primate Research Center, Texas, USA, Copenhagen Zoo, Denmark, Knuthenborg Park, Denmark and African Green Monkeys at Saint Kitts were recovered during post mortem examination which is performed routinely at these places on culled animals.
T. suis was obtained from experimentally infected pigs in Denmark and US. The Animal Experiments Inspectorate, Ministry of Justice, DK approved the animal study protocol and was carried out according to their guidelines (Licence no. 2005/561-1060). In US the animals were handled and managed according to a protocol (Protocol #07-011) specifically approved by the USDA-ARS Beltsville Area Animal Care and Use Committee and followed the Institutional Animal Care and Use Committees (IACUC) guidelines. In Uganda T. suis was obtained from naturally infected pigs raised on private farms, slaughtered and used for local consumption. Permission to recover worms was obtained from the owners. In Ecuador, worms were collected after anthelmintic treatment of a naturally infected pig raised on a private farm. Written consent to recover worms was obtained from the owner. T. suis were obtained from natural infected pigs in China when slaughtered at local abattoirs. Collection was approved by the Animal Ethics Committee of the Lanzhou Veterinary Research Institute, Chinese Academy of Agricultural Sciences.
Parasite isolates, DNA extraction and typing of the worms
A total of 138 worms were collected from humans, pigs and non-human primates from different regions (Table 1). All worms were rinsed with tap water and stored in 70 % ethanol at 5 °C until DNA extraction. The MasterPure DNA Purifications Kit (Epicentre Biotechnologies) was used to extract DNA according to manufacturer’s protocol after homogenization in 300 μl of lysis solution and overnight incubation at 56 °C except for the populations from Ecuador for which the DNA was extracted previously . Worms were initially typed to confirm worm species by Polymerase Chain Reaction-Restriction Fragment Length Polymorphism (PCR-RFLP) on internal transcribed spacer-2 (ITS-2) following the protocol described by . A negative water control was included in all runs. PCR products and RFLP fragments were stained by GelRed (Biotium) and visualized under UV light in 1.5 % agarose gel. All Trichuris from baboons and humans showed banding pattern characteristic of Trichuris from primates, while all T. suis worms showed the banding pattern characteristic of T. suis.
Amplification of genetic markers and sequencing
Partial sequences of the two mitochondrial genes rrnLnad1 were obtained for all samples except for Trichuris from African green monkey, for which only the rrnL gene was sequenced. 562 bp of the nad1 gene was amplified using forward SuiND1_F (5′-CGAGCTTATATAGGTATTTCTCAACG-3′) and reverse SuiND1_R (5′-CGTTGTAGCCTCTTACTAATTCTCTTT-3′) primers while 422 bp of the rrnL gene was amplified using primers, forward TrirrnL_F (5′-TGTAAWTCTCCTGCCCAATGA) and reverse TrirrnL_R (5′-CGGTTTAAACTCAAATCACGTA). The PCR conditions were identical for both markers and were conducted in a total volume of 20 μl using 1 μl DNA as template. PCR ingredients were: 1X PCR buffer, 0.2 mM of each dNTP, 0.4 mM of each primer pair, 2.0 mM MgCl2, and 1 U of Hot Start DNA-polymerase (Ampliqon). PCR conditions were initial denaturation at 95 °C for 15 min followed by 35 cycles consisting of 95 °C for 30 s, 55 °C for 30 s and 72 °C for 1 min and a final extension at 72 °C for 10 min. Agarose gel electrophoresis (1.5 %) was used to verify amplification of a single fragment of the expected size. PCR products were enzymatically cleaned prior to sequencing using 10 μl of PCR product, 1 μl Exonuclease I and 2 μl Fast AP Thermosensitive Alkaline Phosphotase (1 U/μl) (Fermentas). The samples were incubated for 15 min at 37 °C followed by 15 min at 85 °C. Finally, cleaned amplicons were sequenced in both directions using same primers used for PCR by Macrogen Inc. in Seoul, South Korea.
Genetic variation and phylogenetic relationships
Forward and reverse sequences from each sample were checked, edited manually and assembled using vector NTI  and then trimmed using BioEdit . Genetic relatedness and evolutionary relationship were analysed for each of the two markers using 410 bp and 397 bp of the nad1 and rrnL genes, respectively. Two sequences of the nad1 and rrnL from the Chinese human and pig Trichuris mitochondrial genomes (Accession No: GU385218 and GU070737, respectively) were included in the dataset. Also, T. trichiurarrnL gene sequences from humans in China were included in the dataset (Accession no. AM993017-AM993021). Phylogenetic relationship was inferred using NJ and ML phylogenetic trees in MEGA v6.1 . The best-to-fit substitution model was identified using jModelTest0.1.1  under Akaike information criterion (AIC) . Trichinellaspiralis was used as outgroup (Accession No: AF293969). The most parsimonious network was inferred by the Neighbor-net method using SplitsTree v.4.13.1 . Neighbor-net network can reveal ambiguous and incompatible sites which usually appear as a reticulate structure in the network.
Population genetic structure
Analysis of Molecular Variance (AMOVA) was used to estimate the Fixation index, Fst between Trichuris populations using Arlequin v.184.108.40.206 . The Fst for the human and pig derived worms were calculated in separate analysis using the geographical origin of the worms to define the different populations. As the origin of Trichuris collected from baboons was unknown due to previous transport between zoological gardens, we used the phylogenetic analysis to define two populations for which Fst subsequently was estimated for. Lastly, Fst was estimated between the human and baboon worms that were found to cluster together. 10,000 permutations were used to test for differentiation between pairs of populations. The p-distance between the distinct clades identified in the phylogenetic analysis was calculated using MEGA v.6.1. .
Demography, time of divergence and TMRCA
For uniparentally inherited DNA, the time to the most recent common ancestor (TMRCA) in generations is equal to the population size (N). The effective population sizes for the populations of T. trichiura and T. suis were calculated using the formula Θ = 2Neffμ where Θ (theta) is the genetic diversity of a population, μ is the mutation rate per gene and Neff is the effective population size. Theta (Θ) and the ancestral history were estimated from Genetree  using a concatenated dataset of the two markers. First, sequences were aligned and imported to Map modules in the SNAP workbench  to collapse the sequences to haplotypes excluding sites which are indels and infinite site violations. Then, compatibility analysis using CladeEx revealed incompatible sites which were removed . The simulations were repeated 5 times with 10 million runs with different random seeds to ensure convergence of the genealogies. The mutation rate of Caenorhabditiselegans was used which is 1.6 X 10−7 per site per generation  and found not to be significantly different from other free living nematodes . To obtain the mutation rate per gene, the mutation rate per site was multiplied by the number of nucleotides used (807 nt for both markers) giving 1.29 X 10−4 per gene per generation.
BEAST v.1.6.1  was used to infer the phylogeny and the divergence time using the Bayesian statistical framework and the concatenated dataset with Trichuris from pigs and humans as monophyletic groups. Different mutation models were used and the final analysis was done using strict molecular clock with a normal distributed substitution rate of 1.6 X 10−7(±0.3 X 10−7) based on the value of C. elegans . The substitution model used here was Hasegawa-Kishino-Yano (HKY) with gamma distribution as best to fit model based on AIC  in jModelTest0.1.1 . Yule prior, which is suitable for datasets that combine different species, was used as a tree prior with a random starting tree. Markov chain Monte Carlo (MCMC) chains were run for iterations with a burn in value of 1000. Tracer v.1.6 was used to analyse log files and to check whether the MCMC chains were sufficient by recording effective sample size values to be above 200, which was the case for all the parameters. The three log files of the three independent runs were combined using log combiner v1.6.1 . Tree Annotater v1.6.1.  was used to summarize samples from the posterior on maximum credibility tree and the posterior probability limit set to 0.5. Figtree v1.3.1  was used to depict the tree.
Divergence time was also estimated for the human and pig Trichuris populations using IMa2 based on an isolation and migration model  using the concatenated dataset. Priors used for these data sets were: for T. trichiura, t (the upper bound of splitting time) = 295, q (upper bound of population size) = 750, while for T. suis t = 40, q = 100. For both data sets, HKY substitution model was used, no migration between populations after splitting was assumed (m = 0), 20 Markov chains with geometric heating scheme (the first and second heating parameters were 0.96 and 0.90, respectively) and 106 burn-in steps with 105 sampling genealogies were used. Three independent runs were conducted with different seed numbers to assess the convergence.
We used two mitochondrial markers to infer the genetic relationships and the evolutionary history between different Trichuris populations obtained from humans, non-human primates and pigs from various geographical regions (Table 1). Partial sequences of the large ribosomal subunit (rrnL) and NADH dehydrogenase subunit 1 (nad1) were generated for all worms except for Trichuris from African green monkey, for which only the rrnL gene was sequenced (GenBank accession numbers: KU524489-KU524606 and KU524607-KU524730). The phylogenetic relationship was found to be identical when inferred using neighbor joining (NJ) and maximum likelihood (ML) clustering methods and for both genes, and hence only the NJ tree for the rrnL gene is shown in Fig. 1. For the ML tree, Hasegawa-Kishino-Yano with gamma distribution (HKY + G) model was used for the nad1 gene and Tamura-Nei with gamma distribution (TrN + G) model for the rrnL gene as the best-to-fit substitution models. It is noteworthy that two samples (two T. trichiura from humans, one in Uganda and one in China) showed double peaks in the chromatogram of the nad1 gene and gave conflicting signals for the rrnL gene (clustered in the pig clades). This may indicate co-amplification of nuclear mitochondrial pseudogenes (numts) or a heteroplasmy. Hence, these samples were excluded from further analyses.
Phylogeographic distribution among T. suis populations was observed as worms from Uganda and from China were found in separate clades whereas T. suis from Denmark and USA clustered together. The Ecuadorian T. suis were found in two clades, namely the Denmark & USA clade (4 samples) and the China clade (3 samples). In addition, two pig worms from Uganda clustered with T. suis from Denmark & USA. Phylogeographic structure was also observed for worms recovered from humans as T. trichiura from Uganda, China and Ecuador were found in separate clades. The majority of the baboon and green monkey worms clustered with the human T. trichiura from Uganda. Additionally, a few baboon samples (n = 7) and a single green monkey worm grouped together in a different clade (Trichuris non-human primates). The Neighbor-net network identified splits that correspond to the clades in the NJ tree (Fig. 2).
Population genetic structure (Fst)
AMOVA was used to analyse the degree of genetic differentiation (Fst) between Trichuris populations and is summarized in Tables 2 and 3 and is given for each of the two markers, but excluding T. suis from Ecuador as they cluster both with worms from China and Denmark. In general all the populations were highly differentiated as Fst values are above 0.25. Trichuris from baboons and green monkey in the two different clades (Trichuris African non-human primates and Trichuris African Primates) were highly differentiated with Fst values of 0.363 for rrnL (P < 0.01) and 0.471 for nad1 (P < 0.01). In contrast, human, baboon and green monkey Trichuris in the T. trichiura Afican clade represented undifferentiated populations (Fst < 0.05, P > 0.05). The genetic distances (p-distance) between and within each clade are given in Additional file 1: Table S1 for the nad1 and rrnL genes, respectively.
TMRCA and the divergence time of the human and pig Trichuris populations
The estimated Θ was 20.7 ± 6.68 for T. suis and 146.49 ± 50.1 for T. trichiura from humans and the Genetree for their populations is provided in Additional file 2: Figure S1. Hence, the TMRCA in generations for T. suis equals 80,000 generations (upper estimate is 110,000 and lower estimate 60,000) and 560,000 generations for T. trichiura (upper estimate is 760,000 and lower estimate is 390,000). The time of divergence between the different populations as estimated by BEAST is given at each node in number of generations (Fig. 3) and by IMa2 in Additional file 3: Figure S2. BEAST, IMa2 and Genetree gave similar results for the time of divergence of the human Trichuris populations. However, the estimated time of divergence for T. suis populations was nearly three times older for BEAST and IMa2 than Genetree. Genetree relies on importance sampling (IS) algorithm for the coalescent simulations while BEAST and IMa2 each rely on correlated sampling (CS). As the IS algorithm is more suitable for data of low polymorphism , its estimate may be the most reliable in our case. There were two main divergence events in T. suis and T. trichiura populations. For T. suis there was an ancient split between the USA/Denmark populations and the China/Uganda populations (80,000-240,000 generations) and a more recent split between the populations from China and Uganda (32,000-90,000 generations). For T. trichiura, the first split is between the Uganda population and China/Ecuador populations (500,000 generations) and a more recent split between the China and Ecuador populations (120,000 generations).
Herein, we investigated the evolutionary and genetic relationships between populations of Trichuris from primates and pigs from different geographical areas in order to infer the evolutionary and demographic history of the parasites, and to identify possible environmental and anthropological factors driving their spread across the globe. In this study Trichuris from humans and pigs were genetically very distinct with independent demographic histories as summarised in Fig. 4 and discussed below.
The coalescent analysis identified two divergence events in T. suis populations. First an ancient split between the DK/USA and the Chinese/Ugandan populations around 80,000 generations ago with a second more recent divergence between the China and Uganda populations 32,000 generations ago. Interestingly, this latter split between T. suis from China and Uganda is in line with that of their host . Common alleles between the domesticated pigs of Far East and East African origin have been identified suggesting close evolutionary relationships between pigs in these regions [37, 38]. First, this may reflect transport of pigs from the Far East by the European trading routes to Africa a few hundred years ago. Secondly, domesticated pigs may have been introduced to Africa by trading between the ancient civilizations in Africa and Far East or the settlement of Austronesian peoples in East Africa nearly seven thousand years ago . Given the high genetic differentiation between the T. suis populations of China and Uganda and assuming one to three generations for T. suis per year, the divergence between Chinese and Uganda T. suis population happened 32,000-10,700 years ago. Hence, it is unlikely that T. suis was introduced by an European intermediary only a few hundred years ago but may be traced back to the introduction of domesticated pigs from the Far East thousands of years ago. However, this does not exclude other waves of recent introduction of T. suis in domesticated pigs from the Far East. A recent study found that A. lumbricoides in humans on Zanzibar were closely related to worms from Bangladesh  suggesting parasite transportation between the Far East and the Indian subcontinent and East Africa.
Intriguingly, the T. suis population in Uganda was found to be monomorphic for both markers. This may relate to either a founder effect (the establishment of a new population from few individuals derived from a much larger population) or a selective sweep (strong positive natural selection of some few genotypes) or a combination of both factors. Such a selective sweep might have occurred due to new adaptations to host physiology in the new environment or, more likely, due to a recent bottleneck in the pig populations: e.g. African swine fever outbreaks in Uganda  resulted in subsequent reductions in molecular variation among the parasites . However, due to our small sample size further worms should be analysed to confirm this hypothesis.
Several studies have reported high genetic distinctiveness between European and Chinese pigs [41, 42] in line with our observation for T. suis from Denmark and China. However, the divergence time between the pigs in China and Europe has been estimated to be roughly a million years ago  which is much older than the divergence time estimated for the Chinese and Danish T. suis populations in our study (80,000 years assuming one generation/year). In principle, when populations suffer from many bottlenecks, the coalescent will date back to the most recent bottleneck rather than the most recent common ancestor . Considering the bottlenecks that the pig populations have been through during their migration and transport across Eurasia, with the most recent being nearly 20,000 years ago during the last glacial maxima , this may also have resulted in bottlenecks of the associated parasites. In addition, as T. suis can survive in the environment for at least 11 years  the number of generations per year is very hard to estimate and could be as low as 0.09 generations/year.
T. suis from Denmark and USA were found to cluster together (Fig. 1) with their populations being undifferentiated. This is consistent with other parasites of domestic pigs such as Trichinellaspiralis, which was introduced to the Americas by Europeans . However, two T. suis from Uganda clustered with worms from Denmark and USA suggesting recent transport of pigs between these continents. The clustering of T. suis from Ecuador with the populations from China and DK/USA may reflect introgression between pigs from Europe and China during the industrial revolution in Europe in the 20th century as found in a previous study .
As for pig worms two divergence events were also observed for human derived Trichuris, one ancient (~500,000 generations) divergence between Uganda and China and a more recent one between China and Ecuador (~120,000 generations). Assuming the highest and lowest number of generations to be 1 and 3, respectively the divergence times were 500,000-160,000 generations and 120,000-40,000 generations for the Ugandan/Chinese and Chinese/Ecuadorian populations, respectively. The split between the China and Uganda populations preceded the modern human (Homo sapiens) migration out of Africa to South East Asia around 60,000-100,000 years ago [45, 46] and the human settlement in Latin America 14,000-15,000 years ago . There are two possibilities for this discrepancy. Firstly, one of the early human ancestors (e.g. H. erectus) may have transmitted T. trichuris when migrating out of Africa  but this would not explain how the parasite was then introduced into Latin America. Secondly, the mutation rate of the free living nematode (C. elegans) may not be applicable to parasitic nematodes as the mutation rate normally is higher for the latter . T. trichiura eggs in Brazil dating back to 6000–7000 BP  suggests that the parasite was introduced to the New World with the human migration much earlier than the Columbian colonization , which is concordant with our findings.
Green monkeys, olive, yellow and hamadryas baboons are indigenous species in central and western Africa and the introduction of the green monkey into Saint Kitts by the French in the late Seventeenth century  may explain why the majority of the baboon and green monkey worms clustered with human worms from Uganda. However, seven worms from baboons from both Denmark and USA and one from green monkey from St. Kitts were found in a separate clade (Trichuris African non-human primate), suggesting that different populations are circulating among these hosts’ species although they were sampled from the same habitat. Since many of these worms were sampled from unnatural habitats (zoos), the causes for this differentiation could not be investigated. However, in a sympatric natural transmission area in Uganda some Trichuris genotypes were found among all seven investigated non-human primates, whereas other genotypes seem to be more host specific . Hence, in the past, Trichuris in primates in Africa may have been isolated either by geography or by host species leading to population differentiation followed by a more recent secondary sympatry . As the Trichuris population between humans from Uganda and the non-human primates are undifferentiated this suggests continuous or recent gene flow between the host species and suggests Africa as the origin of Trichuris in primates.
Unlike whipworms, it is expected that host shift of the giant roundworm Ascaris in humans and pigs took place during animal domestication in the Neolithic period 10,000 years ago [7, 12]. This is supported by their very close genetic relationship [7, 54] and that Ascaris among non-human primates is observed rarely , suggesting that this parasite represents a ‘souvenir parasite’ . Their different evolutionary history is interesting as Trichuris and Ascaris share several parasitic traits such as mode of transmission and infection and these two parasite species therefore usually co-occur today .
We inferred the possible demographic history of Trichuris from humans and pigs and their potential evolutionary epicenter, namely in primates in Africa. We suggest that Trichuris was dispersed to Asia with human ancestors and that host switching to pigs occurred in China where pigs evolved . T. suis was then spread across the globe mainly by anthropogenic factors. We found that T. trichiura in humans in Africa is genetically similar to Trichuris in non-human primates of African origin suggesting that Trichuris in humans represents a heirloom parasite. Further studies should investigate the genetic relationships of whipworms from different primates and other host species living in natural habitats in order to explore their demographic history and evolutionary origins and to identify the possible switching events between host species. Moreover, since we relied only on mitochondrial markers to resolve the evolutionary history of the parasite which in certain cases might be problematic due to the presence of the mitochondrial pseudogenes (numts), heteroplasmy or incomplete lineage sorting  further studies should also include suitable nuclear markers to confirm our findings.
Pullan RL, Smith JL, Jasrasaria R, Brooker SJ. Global numbers of infection and disease burden of soil transmitted helminth infections in 2010. Parasit Vectors. 2014;7:37. doi:10.1186/1756-3305-7-37.
Abee CR, Mansfield K, Tardif SD, Morris T. Nonhuman primates in biomedical research diseases. 2nd Edition ed. Boston: Academic Press; 2012.
Nansen P, Roepstorff A. Parasitic helminths of the pig: factors influencing transmission and infection levels. Int J Parasitol. 1999;29(6):877–91.
Roepstorff A, Mejer H, Nejsum P, Thamsborg SM. Helminth parasites in pigs: new challenges in pig production and current research highlights. Vet Parasitol. 2011;180(1–2):72–81. doi:10.1016/j.vetpar.2011.05.029.
Archie EA, Luikart G, Ezenwa VO. Infecting epidemiology with genetics: a new frontier in disease ecology. TrendsEcol Evol. 2009;24(1):21–30. doi:10.1016/j.tree.2008.08.008.
Gilabert A, Wasmuth JD. Unravelling parasitic nematode natural history using population genetics. Trends Parasitol. 2013;29(9):438–48. doi:10.1016/j.pt.2013.07.006.
Betson M, Nejsum P, Bendall RP, Deb RM, Stothard JR. Molecular epidemiology of ascariasis: a global perspective on the transmission dynamics of Ascaris in people and pigs. J Infect Dis. 2014;210(6):932–41. doi:10.1093/infdis/jiu193.
Ramesh A, Small ST, Kloos ZA, Kazura JW, Nutman TB, Serre D, et al. The complete mitochondrial genome sequence of the filarial nematode Wuchereria bancrofti from three geographic isolates provides evidence of complex demographic history. Mol Biochem Parasitol. 2012;183(1):32–41. doi:10.1016/j.molbiopara.2012.01.004.
Cox FE. History of human parasitology. Clin Microbiol Rev. 2002;15(4):595–612.
Goncalves ML, Araujo A, Ferreira LF. Human intestinal parasites in the past: new findings and a review. Mem Inst Oswaldo Cruz. 2003;98 Suppl 1:103–18.
Mitchell DP. The origins of human parasites: exploring the evidence forendoparasitism throughout human evolution. Int J Paleopathol. 2013;3:191–8.
Buikstra JE, Roberts CA. A global history of paleopathology: pioneers and prospects. New York: Oxford University Press; 2012.
Ghai RR, Simons ND, Chapman CA, Omeja PA, Davies TJ, Ting N, et al. Hidden population structure and cross-species transmission of whipworms (Trichuris sp.) in humans and non-human primates in Uganda. PLoS Negl Trop Dis. 2014;8(10), e3256. doi:10.1371/journal.pntd.0003256.
Liu GH, Gasser RB, Nejsum P, Wang Y, Chen Q, Song HQ, et al. Mitochondrial and nuclear ribosomal DNA evidence supports the existence of a new Trichuris species in the endangered francois’ leaf-monkey. PLoS One. 2013;8(6), e66249. doi:10.1371/journal.pone.0066249.
Cutillas C, Callejon R, de Rojas M, Tewes B, Ubeda JM, Ariza C, et al. Trichuris suis and Trichuristrichiura are different nematode species. Acta Trop. 2009;111(3):299–307. doi:10.1016/j.actatropica.2009.05.011.
Ravasi DF, O’Riain MJ, Davids F, Illing N. Phylogenetic evidence that two distinct Trichuris genotypes infect both humans and non-human primates. PLoS One. 2012;7(8), e44187. doi:10.1371/journal.pone.0044187.
Callejon R, Cutillas C, Nadler SA. Nuclear and mitochondrial genes for inferring Trichuris phylogeny. Parasitol Res. 2015;114(12):4591–9. doi:10.1007/s00436-015-4705-7.
Dolezalova J, Obornik M, Hajduskova E, Jirku M, Petrzelkova KJ, Bolechova P et al. How many species of whipworms do we share? Whipworms from man and other primates form two phylogenetic lineages. Folia Parasitol. 2015;62. doi:10.14411/fp.2015.063
Liu GH, Gasser RB, Su A, Nejsum P, Peng L, Lin RQ, et al. Clear genetic distinctiveness between human- and pig-derived Trichuris based on analyses of mitochondrial datasets. PLoS Negl Trop Dis. 2012;6(2), e1539. doi:10.1371/journal.pntd.0001539.
Hawash MB, Andersen LO, Gasser RB, Stensvold CR, Nejsum P. Mitochondrial genome analyses suggest multiple Trichuris species in humans, baboons, and pigs from different geographical regions. PLoS Negl Trop Dis. 2015;9(9), e0004059. doi:10.1371/journal.pntd.0004059.
Meekums H, Hawash MB, Sparks AM, Oviedo Y, Sandoval C, Chico ME, et al. A genetic analysis of Trichuris trichiura and Trichuris suis from Ecuador. Parasit Vectors. 2015;8:168. doi:10.1186/s13071-015-0782-9.
Nissen S, Al-Jubury A, Hansen TV, Olsen A, Christensen H, Thamsborg SM, et al. Genetic analysis of Trichuris suis and Trichuris trichiura recovered from humans and pigs in a sympatric setting in Uganda. Vet Parasitol. 2012;188(1–2):68–77. doi:10.1016/j.vetpar.2012.03.004.
Lu G, Moriyama EN. Vector NTI, a balanced all-in-one sequence analysis suite. Brief Bioinform. 2004;5(4):378–88.
Hall TA. BioEdit: a user-friendly biological sequence alignment editor and analysis program for Windows 95/98/NT. Nucl Acids Symp Ser. 1999;41:95–8.
Tamura K, Stecher G, Peterson D, Filipski A, Kumar S. MEGA6: molecular evolutionary genetics analysis version 6.0. Mol Biol Evol. 2013;30(12):2725–9. doi:10.1093/molbev/mst197.
Posada D. jModelTest: phylogenetic model averaging. Mol Biol Evol. 2008;25(7):1253–6. doi:10.1093/molbev/msn083.
Akaike H. A new look at the statistical model identification. IEEE Trans Autom Control. 1974;19(6):716–23. doi:10.1109/TAC.1974.1100705.
Huson DH, Bryant D. Application of phylogenetic networks in evolutionary studies. Mol Biol Evol. 2006;23(2):254–67. doi:10.1093/molbev/msj030.
Excoffier L, Lischer HE. Arlequin suite ver 3.5: a new series of programs to perform population genetics analyses under Linux and Windows. Mol Ecol Res. 2010;10(3):564–7. doi:10.1111/j.1755-0998.2010.02847.x.
Griffiths RC, Tavare S. Sampling theory for neutral alleles in a varying environment. Phil Trans R Soc B. 1994;344(1310):403–10. doi:10.1098/rstb.1994.0079.
Aylor DL, Price EW, Carbone I. SNAP: Combine and Map modules for multilocus population genetic analysis. Bioinformatics. 2006;22(11):1399–401. doi:10.1093/bioinformatics/btl136.
Denver DR, Morris K, Lynch M, Vassilieva LL, Thomas WK. High direct estimate of the mutation rate in the mitochondrial genome of Caenorhabditis elegans. Science. 2000;289(5488):2342–4.
Molnar RI, Bartelmes G, Dinkelacker I, Witte H, Sommer RJ. Mutation rates and intraspecific divergence of the mitochondrial genome of Pristionchus pacificus. Mol Biol Evol. 2011;28(8):2317–26. doi:10.1093/molbev/msr057.
Drummond AJ, Rambaut A. BEAST: Bayesian evolutionary analysis by sampling trees. BMC Evol Biol. 2007;7:214. doi:10.1186/1471-2148-7-214.
Hey J, Nielsen R. Multilocus methods for estimating population sizes, migration rates and divergence time, with applications to the divergence of Drosophila pseudoobscura and D. persimilis. Genetics. 2004;167(2):747–60. doi:10.1534/genetics.103.024182.
Kuhner MK. Coalescent genealogy samplers: windows into population history. Trends Ecol Evol. 2009;24(2):86–93. doi:10.1016/j.tree.2008.09.007.
Amills M, Ramírez O, Galman-Omitogun O, Clop A. Domestic pigs in Africa. AfrArchaeolRev. 2013;30(1):73–82. doi:10.1007/s10437-012-9111-2.
Ramirez O, Ojeda A, Tomas A, Gallardo D, Huang LS, Folch JM, et al. Integrating Y-chromosome, mitochondrial, and autosomal data to analyze the origin of pig breeds. Mol Biol Evol. 2009;26(9):2061–72. doi:10.1093/molbev/msp118.
Atuhaire DK, Afayoa M, Ochwo S, Mwesigwa S, Okuni JB, Olaho-Mukani W, et al. Molecular characterization and phylogenetic study of African swine fever virus isolates from recent outbreaks in Uganda (2010–2013). Virol J. 2013;10:247. doi:10.1186/1743-422X-10-247.
Barrett LG, Thrall PH, Burdon JJ, Linde CC. Life history determines genetic structure and evolutionary potential of host-parasite interactions. Trends Ecol Evol. 2008;23(12):678–85. doi:10.1016/j.tree.2008.06.017.
Giuffra E, Kijas JM, Amarger V, Carlborg O, Jeon JT, Andersson L. The origin of the domestic pig: independent domestication and subsequent introgression. Genetics. 2000;154(4):1785–91.
Groenen MA, Archibald AL, Uenishi H, Tuggle CK, Takeuchi Y, Rothschild MF, et al. Analyses of pig genomes provide insight into porcine demography and evolution. Nature. 2012;491(7424):393–8. doi:10.1038/nature11622.
Burden DJ, Ginnivan MJ. The destruction of pig helminth ova and larvae in a slurry treatment process. Vet Rec. 1978;103(17):373–5.
Rosenthal BM, LaRosa G, Zarlenga D, Dunams D, Chunyu Y, Mingyuan L, et al. Human dispersal of Trichinella spiralis in domesticated pigs. Infect Genet Evol. 2008;8(6):799–805. doi:10.1016/j.meegid.2008.07.008.
Liu W, Jin CZ, Zhang YQ, Cai YJ, Xing S, Wu XJ, et al. Human remains from Zhirendong, South China, and modern human emergence in East Asia. Proc Natl Acad Sci U S A. 2010;107(45):19201–6. doi:10.1073/pnas.1014386107.
Beyin A. Upper Pleistocene human dispersals out of Africa: a review of the current state of the debate. J Evol Biol. 2011;2011:615094. doi:10.4061/2011/615094.
Dillehay TD. The late Pleistocene cultures of South America. Evol Anthropol. 1999;7:206–16.
Finlayson C. Biogeography and evolution of the genus Homo. Trends Ecol Evol. 2005;20(8):457–63. doi:10.1016/j.tree.2005.05.019.
Holterman M, van der Wurff A, van den Elsen S, van Megen H, Bongers T, Holovachov O, et al. Phylum-wide analysis of SSU rDNA reveals deep phylogenetic relationships among nematodes and accelerated evolution toward crown Clades. Mol Biol Evol. 2006;23(9):1792–800. doi:10.1093/molbev/msl044.
Brinkkemper O, van Haaster H. Eggs of intestinal parasites whipworm (Trichuris) and mawworm (Ascaris): Non-pollen palynomorphs in archaeological samples. Rev Palaeobot Palynol. 2012;186:16–21. doi:10.1016/j.revpalbo.2012.07.003.
Araujo A, Reinhard KJ, Ferreira LF, Gardner SL. Parasites as probes for prehistoric human migrations? Trends Parasitol. 2008;24(3):112–5. doi:10.1016/j.pt.2007.11.007.
KA CL. Primate Factsheets: Vervet (Chlorocebus) Taxonomy, Morphology, & Ecology 2006. http://pin.primate.wisc.edu/factsheets/entry/vervet. Accessed June 30 2015.
Ricklefs RE. Host-pathogen coevolution, secondary sympatry and species diversification. Phil Trans R Soc B. 2010;365(1543):1139–47. doi:10.1098/rstb.2009.0279.
Cavallero S, Snabel V, Pacella F, Perrone V, D’Amelio S. Phylogeographical studies of Ascaris spp. based on ribosomal and mitochondrial DNA sequences. PLoS Negl Trop Dis. 2013;7(4):e2170. doi:10.1371/journal.pntd.0002170.
Nejsum P, Grondahl C, Murrell KD. Molecular evidence for the infection of zoo chimpanzees by pig Ascaris. Vet Parasitol. 2006;139(1–3):203–10. doi:10.1016/j.vetpar.2006.02.025.
Leles DRK, Fugassa M, Ferreira LF, Iniguez AM, Araújo A. A parasitological paradox: why is ascarid infection so rare in the prehistoric Americas? J Archaeol Sci. 2010;37:1510–20.
McGuire JA, Linkem CW, Koo MS, Hutchison DW, Lappin AK, Orange DI, et al. Mitochondrial introgression and incomplete lineage sorting through space and time: phylogenetics of crotaphytid lizards. Evolution. 2007;61(12):2879–97. doi:10.1111/j.1558-5646.2007.00239.x.
Hansen TV, Thamsborg SM, Olsen A, Prichard RK, Nejsum P. Genetic variations in the beta-tubulin gene and the internal transcribed spacer 2 region of Trichuris species from man and baboons. Parasit Vectors. 2013;6:236. doi:10.1186/1756-3305-6-236.
We gratefully acknowledge the following for providing worm specimens: J. Urban, United States Department of Agriculture, USA; H. Namwanje, Vector Control Division, Ministry of Health, Kampala, Uganda; S. Nissen, I-H. Poulsen, A. Andreassen, H.H. Petersen University of Copenhagen, Denmark; T. J. C. Anderson and the Southwest National Primate Research Center at Texas Biomedical Research Institute, San Antonio, Texas; Amy Beierschmitt at the Behavioural Science Foundation, St. Kitts. P.J.C. was supported by Wellcome Trust grant 088862/Z/09/Z. This work is supported by a Danish Agency for Science, Technology and Innovation grant awarded to P.N. (FSS 09–066344).
The authors declare that they have no competing interests.
PN and MBFH conceived and designed the study with inputs from MB and DTJL; JK, ALW, MFB, PJC, XQZ provided essential material, MBFH and AAG conducted the molecular work. MBFH and PN analysed the data. MBFH and PN wrote the paper with inputs from MB and DTJL. All authors critically read and approved the final manuscript.
The gene genealogy inferred by Genetree of (A) T. suis populations and (B) T. trichiura populations. Solid circles indicate mutations in the genealogy. (JPG 1237 kb)
Splitting time based on the isolation and migration model between (A) T. suis populations and (B) T. trichiura populations. The horizontal axis represents the number of generations since splitting which were estimated by dividing the splitting times between populations (t0 and t1) and time to most recent common ancestor (tmrca) by the mutation rate per gene per generation (μ) while the vertical axis is the posterior probability density. (JPG 149 kb)