Molecular evidence of hybridization between pig and human Ascaris indicates an interbred species complex infecting humans
Abstract
Human ascariasis is a major neglected tropical disease caused by the nematode Ascaris lumbricoides. We report a 296 megabase (Mb) reference-quality genome comprised of 17,902 protein-coding genes derived from a single, representative Ascaris worm. An additional 68 worms were collected from 60 human hosts in Kenyan villages where pig husbandry is rare. Notably, the majority of these worms (63/68) possessed mitochondrial genomes that clustered closer to the pig parasite Ascaris suum than to A. lumbricoides. Comparative phylogenomic analyses identified over 11 million nuclear-encoded SNPs but just two distinct genetic types that had recombined across the genomes analyzed. The nuclear genomes had extensive heterozygosity, and all samples existed as genetic mosaics with either A. suum-like or A. lumbricoides-like inheritance patterns supporting a highly interbred Ascaris species genetic complex. As no barriers appear to exist for anthroponotic transmission of these ‘hybrid’ worms, a one-health approach to control the spread of human ascariasis will be necessary.
Introduction
Approximately 447 million people were estimated to be infected with the intestinal nematode Ascaris lumbricoides in 2017, resulting in an estimated 3206 deaths and a loss of over 860,000 Disability-Adjusted Life Years (DALYs, Global Burden of Disease Study, 2017; http://ghdx.healthdata.org/gbd-2017). Many infections go undiagnosed, but like other soil-transmitted helminths (STH), Ascaris spp. infections contribute significantly to global DALYs, perpetuating the cycle of poverty in areas of endemic infection (Brooker, 2010; Hotez et al., 2008; Montresor et al., 2012; Pullan et al., 2014). Despite the large global burden of STH, little is known about A. lumbricoides transmission patterns or the true prevalence of infection with the pig parasite A. suum infection in humans in endemic regions.
Deworming has become more widespread in areas of endemic STH infection (Bundy et al., 2017). Regional health authorities and global health organizations are now looking for strategies to build on these programs by achieving local elimination of STH as a public health problem (Becker et al., 2018). A greater understanding of transmission dynamics (including the frequency of zoonotic transmission) using molecular epidemiological methods in settings where A. lumbricoides prevalence is low but persistent could help move current efforts toward successfully eliminating transmission through more targeted treatment.
Population genetic studies of A. lumbricoides have drawn varying conclusions about whether zoonotic transmission is frequent (Anderson and Jaenike, 1997; Dutto and Petrosillo, 2013; Nejsum et al., 2012; Nejsum et al., 2005a). Some studies have shown that cross-species transmission occurs between pigs and humans living in close proximity (Anderson, 1995; Betson et al., 2014; Miller et al., 2015; Monteiro et al., 2019; Nejsum et al., 2005b; Peng and Criscione, 2012; Sadaow et al., 2018; Takata, 1951; Zhu et al., 1999). This is especially common in non-endemic regions, probably because zoonotic transmission is less likely to be identified in areas where human-to-human transmission is common. The human parasite A. lumbricoides and the pig parasite A. suum have been found to be capable of interbreeding, and 4–7% of worms in Guatemala and China were hybrids (Criscione et al., 2007; Peng and Criscione, 2012). Furthermore, it is unclear whether pigs are an important reservoir of infection in humans worldwide or if A. suum is readily transmitted anthroponotically (Betson et al., 2013; Betson and Stothard, 2016; da Silva Alves et al., 2016; Leles et al., 2012; Nejsum et al., 2012). Studies have generally concluded that the genetic differences between Ascaris worms collected from human populations in different parts of the world (Betson et al., 2014; Peng et al., 1998) are the result of geographic reproductive isolation. Previous studies using Ascaris mitochondrial genomes or genes suggest there are A. lumbricoides-type (human-associated) and A. suum-type (pig-associated) clades (Anderson and Jaenike, 1997; Cavallero et al., 2013; Zhou et al., 2011). Other work suggests multiple clades of worms, only one of which is unique to pigs (Nejsum et al., 2017). Ascaris spp. infections also occur naturally in monkeys and apes, and Ascaris spp. eggs are sometimes found in the feces of dogs but this is likely a result of coprophagy by the dogs, rather than due to infection (https://www.cdc.gov/parasites/ascariasis/biology.html).
In the current study, we constructed a reference-quality Ascaris genome (ALV5) based on sequences from a single female worm collected from a person in Kenya. This person was presumed to be infected with A. lumbricoides as there is a lack of local pig husbandry. Draft A. suum genomes have previously been constructed from worms obtained from pigs in Australia (Jex et al., 2011) and in the United States (Wang et al., 2011; Wang et al., 2012). The Ascaris genome ALV5 was found to be highly similar (99% identity) to the A. suum genome from worms collected from pigs in the United States (Wang et al., 2017). Our mitochondrial and whole-genome analyses from an additional 68 individual worms indicate that A. suum and A. lumbricoides form a genetic complex that is capable of interbreeding. Our data support a model for a recent worldwide, multi-species Ascaris population expansion caused by the movement of humans and/or livestock globally. Ascaris from both pigs and humans may be important in human disease, necessitating a one-health approach to control the spread of human ascariasis.
Results
Human Ascaris reference genome to promote comparative genomic analyses
To generate a human Ascaris spp. germline genome assembly (prior to programmed DNA elimination Wang et al., 2017), ovarian DNA was sequenced from a single female worm collected from a Kenyan study participant who was presumed to be infected with A. lumbricoides using Illumina paired-end and mate-pair libraries of various insert sizes with a total sequence coverage of ~27 fold (Supplementary file 1). Using these data, three different assembly strategies were used. The de novo assembly and semi-de novo strategies produced poor A. lumbricoides germline draft genomes (Table 1). In the semi-de novo assembly, the majority of the >4000 short contigs (making up 15.4 Mb of sequence) that could not be incorporated into the semi-de novo assembly are sequences that aligned to the genome at multiple positions. Comparison of the A. suum gene annotations to this assembly revealed a low A. lumbricoides gene number and high numbers of partial and split genes (Table 1, see footnote 3). These characteristics are typical of highly fragmented genomes or genomes with high levels of mis-assemblies (Wang et al., 2017).
Mapping of the human Ascaris reads to the A. suum reference genome (Wang et al., 2017) revealed an exceptionally high-sequence similarity (>99% identity) between the two species with few human Ascaris reads that could not be mapped to A. suum. Based on this high-sequence similarity, a third reference-based-only assembly strategy was used to generate the human Ascaris germline genome assembly using the A. suum germline genome as a reference (see Materials and methods). This approach led to a reference-quality human Ascaris genome assembly with many fewer gaps (only 0.98 Mb of sequence) and no unplaced contigs. The Ascaris genome assembled into 415 scaffolds with a combined size of 296 Mb. An additional 15.4 Mb of sequence was present in 4072 unscaffolded short contigs. The assembly N50 value was 4.63 Mb, with the largest scaffold measuring 13.2 Mb. The largest 50 scaffolds combined to represent 78% of the genome. The assembly was further polished using additional Illumina reads from the same worm to more accurately reflect single base differences, indels, and any potential local mis-assembled regions.
To evaluate the quality of the assembled genome, we mapped the Ascaris Illumina reads back to the reference-based Ascaris genome assembly and found that >99% of the Illumina reads could be mapped, indicating that the reference-based assembly excluded very few Ascaris reads. We then mapped and transferred the extensive set of A. suum transcripts (Jex et al., 2011; Wang et al., 2017) to the human Ascaris germline assembly to annotate the genome, identifying and classifying 17,902 protein-coding genes (Table 1, Supplementary file 1). As this reference-based assembly exhibits the best assembly attributes, including high continuity with a large N50, low gaps and unplaced sequences, and high-quality protein-coding genes (see Table 1), we suggest that this version should be used as a reference germline genome for a human Ascaris spp. specimen (available in NCBI GenBank with accession number PRJNA515325). The other two assemblies are available online.
Like A. suum embryos, A. lumbricoides embryos undergo programmed DNA elimination during the differentiation of the somatic cells from the germline in early development (Streit et al., 2016; Wang and Davis, 2014). In A. suum, ~30 Mb of 120 bp tandem repeats and ~1000 germline-expressed genes are lost from the germline to form the somatic genome (Wang et al., 2012; Wang et al., 2017). We also sequenced the somatic genome from the intestine of the same female A. lumbricoides worm. Comparison of the germline and somatic genomes revealed that DNA elimination in the human Ascaris sample (including the breaks, sequences, and genes eliminated) was identical to that described for the pig A. suum sample (Wang et al., 2017).
Gene content and Ascaris proteome
Earlier annotations of protein coding genes for A. suum draft genomes were produced by Jex et al., 2011 and Wang et al., 2012 and improved with a recent updated genome (Wang et al., 2017)—although the focus of the recent study was not on protein annotations. Here, we updated, identified, and fully annotated the 17,902 protein-coding genes in the reference-based genome assembly (Supplementary file 2 and Figure 1—figure supplement 1). Our aims were to highlight the phylogenetic relationship with other helminths and between Ascaris spp., to provide potential targets for future diagnostics to differentiate between nematodes and even between pig and human Ascaris, and to detail the potential functions of hypothetical or unknown proteins in the Ascaris genome. Using a custom pipeline (see Materials and methods and Cotton et al., 2017), we classified 48% of the predicted proteome into functional groups (Figure 1A). Although the remaining 52% (9300) of the genes were classified as unknown/uncharacterized, 2515 (27%) of these appear to encode proteins that have signatures indicative of either being secreted or being membrane-bound (some with GPI anchors). To provide a more comprehensive annotation of the transcriptomes of A. suum and A. lumbricoides, we re-mapped the RNA-seq data from A. suum to the current gene models of A. lumbricoides (ALV5) (Supplementary file 2). We performed multivariate analyses of this revised RNA-seq data compilation to generate a comprehensive RNA-seq data set for differential gene expression in diverse stages/tissues (Supplementary file 2).
Phylogenetic trees derived from orthologue analyses of the predicted proteomes of ALV5 with the predicted proteomes of other nematodes across all clades indicated the similarity among the published genomes of A. suum PRJNA62057 and PRJNA80881 in Jex et al., 2011; Wang et al., 2012; Wang et al., 2017 and A. lumbricoides (International Helminth Genomes Consortium, 2019) with ALV5 within the Ascaris branch (Figure 1C). The variation observed within the Ascaris spp. (with relatively weak bootstrap values of 0.3–0.59) is likely due to the differences in protein coding gene annotations and split genes seen in previous assemblies.
Mitochondrial genome assembly
We next took advantage of the abundant reads from the mitochondrial genome in our sequencing data (on average 7690X coverage, see Supplementary file 1) to perform de novo assembly of 68 complete human Ascaris spp. mitochondrial genomes from individual worms (Supplementary file 3). These mitochondrial genomes were then annotated using sequence similarity to well-characterized and annotated mitochondrial genes.
Population structure inferred from mitochondrial cox-1 gene
The mitochondrial cox-1 gene has been frequently used to infer evolutionary distances between species as well as between populations (Cavallero et al., 2013; Amor et al., 2016; Springer et al., 2001; Wiens et al., 2010; Zardoya and Meyer, 1996; Zou et al., 2017) due to its rapid mutation rate, lack of recombination and relatively constant rate of change over time (Brown et al., 1979; Giles et al., 1980; Harrison, 1989). Existing data suggest that mitochondria are inherited maternally in C. elegans (Lim et al., 2019; Zhou et al., 2011; Sato and Sato, 2011; Wang et al., 2017) and Ascaris (Anderson et al., 1995). Previous cox-1 phylogeny studies resolve Ascaris spp. worms into three distinct clades: clade A is predominantly comprised of worms isolated from pigs, clade B is predominantly comprised of worms isolated from humans, and clade C is from worms only isolated from pigs in Europe and Asia (Cavallero et al., 2013). Interestingly, haplotype network analyses revealed that the majority of worms isolated from humans in the Kenyan villages possessed cox-1 haplotypes that were consistent with infection of parasites from clade A (63/68), whereas only six specimens had cox-1 haplotypes consistent with infection by worms from clade B (Figure 2—figure supplement 1 and Figure 2a).
When cox-1 sequences from the present study were compared against those within the Ascaris species complex deposited at NCBI (see Supplementary file 4 and Figure 2B; Cotton et al., 2017; Criscione et al., 2007; Godel et al., 2012; Goldberg et al., 2013) within clade A (which appeared to contain the majority of sequences not only from Kenya but also from other localities), seven unique haplotypes of cox-1 from Kenya were identified. These appeared to be shared not only with other haplotypes from Africa, but also with those from Brazil. In contrast, clade B haplotypes appeared to be even more cosmopolitan, with the three haplotypes from Kenya not only being shared with Zanzibar, but also with haplotypes from Brazil, Denmark, China and Japan. Despite the distinct clustering of haplotypes into the three typical Ascaris clades, there was very little genetic diversity among haplotypes within each of the clades, with the majority of haplotypes being separated by 1–4 nucleotide differences. There were greater levels of genetic divergence between clades; A and B were closer to each other while C was more distinct. Similar findings were seen with nad-4, the most variable gene in the mitochondrial genome (Figure 2—figure supplement 1, Figure 2—figure supplement 2).
Phylogenetic analyses and population structure inferred from complete mitochondrial genomes
Forty-seven SNPs were identified in the human Ascaris mitochondrial genomes. Approximately a quarter of these variants were in non-coding portions of the mitochondrial genome and half were synonymous (Supplementary file 1). As with the cox-1 haplotype analyses, whole mitochondrial genome analysis distinguished two clades (clade A and clade B), but there were no distinct geographically specific sub-clades seen within either clade A or clade B (Figure 2B, Table 2). Clade C was also produced by a single published sequence which was used for comparison. In order to assess the validity of the clades A and B representing two distinct molecular taxonomic units, and thus potentially different species, Birky, 2013 4X ratio was applied to provide a lineage-specific perspective of potential species delimitation. The ratio failed to differentiate clades A and B as distinct species with K/Θ <4 at 2.285 indicating Ascaris is one large population—further supporting the lack of differentiation into separate species (Supplementary file 5). Furthermore, there were no significant associations between mitochondrial sequence variations and other factors (e.g. village, household, time of worm collection, host) based on PERMANOVA (see methods and Table 2) after translating the phylogenetic tree into a distance matrix, suggesting not only a lack of differentiation into distinct species but also a potentially large interbreeding population of worms being transmitted between individuals and across villages.
To account for a potentially large population of interbreeding worms, analyses to detect signatures of population expansion were performed. When the global mitochondrial genome data were compared, the Tajima’s D was negative and significant (Tajima’s D −1.5691; p-value 0.028), indicating an excess of low frequency polymorphisms within the global data set suggesting population size expansion. Despite the Fu’s F not being significant it was positive (Fu’s Fs 8.5673; P-0.975) potentially indicating a deficiency in diversity as would be expected in populations that have recently undergone a bottleneck event. The same pattern was also seen in the Kenyan sequences but neither the Tajima’s D nor the Fu’s were significant (Figure 2—figure supplement 3 and Supplementary file 6). Although there does appear to be a signature of a recent population expansion event in both the global and Kenyan data, the lack of information on the mutation rates of Ascaris and other nematodes prevents the accurate estimate of such an event.
Nuclear genome variation in the Ascaris population
To quantify genetic variation in the Ascaris worms isolated from infected Kenyans, the nuclear genomes of the 68 individual worms were analyzed to assess intraspecific population genetic diversity, heterozygosity, and ploidy. Single-nucleotide polymorphisms (SNPs) and insertion/deletions (indels) across the nuclear genomes were assessed for the first 50 largest scaffolds, which comprised 78% of the genome (see methods). Each Ascaris worm was sequenced to a mean coverage depth of ~27 fold. A total of 11.15 million SNP positions were identified in the first 50 scaffolds among the Ascaris nuclear genomes. Approximately 25% of these variants were intergenic (Supplementary file 1). As an example, SNPs and indels in a single Ascaris chromosome were plotted for two worms collected from humans in Kenya and one worm from a pig in the United States (Figure 2—figure supplement 4). The profiles and the frequency between SNPs and indels are highly consistent within individual worms, with the ratio of indel:SNPs frequency at ~1:7. A comparison of the variations identified between individuals infected with worms that had either A. lumbricoides-like or A. suum-like mitochondrial genomes illustrates that most of the differences appear to be random variations, and there do not appear to be major differences between A. lumbricoides-like and A. suum-like worms. A total of 1.79 million SNPs were unique to individual specimens, presumably representing genetic drift. Of the remaining 9.3 million SNPs, ~32% of these variant positions were present in less than five specimens indicating that the Ascaris genomes sequenced are ~1% polymorphic among the major alleles circulating within the species complex.
Population structure inferred from nuclear genomes
To investigate the evolutionary pressures that account for the high SNP diversity found among the 68 sympatric worms, the ploidy, degree of heterozygosity (He) and allelic diversity were determined. Worms were disomic, with little to no evidence of aneuploidy (Figure 3—figure supplement 1). The vast majority (>98%) of SNP positions were biallelic, and each worm had, on average, 2.3 million variant positions, of which approximately 60% were heterozygous SNPs (Supplementary file 7). SNP density was determined in 10 kb windows for each worm against the reference ALV5 and a patchy, mosaic pattern was resolved. SNP density was structured within the genome, with scaffolds being either SNP poor or SNP dense. For example, Algv5r020 was SNP dense whereas Algv5r019x was SNP poor. In other scaffolds, alternating SNP poor and SNP dense regions were defined within the contig, with distinct transition points, see for example the first half of Algv5b02, the last quarter of Algv5b05, or the middle of Algv5r021x (Figure 3A). In those regions where SNP density was low, the Tajima D statistic was net negative, indicating that allele frequencies within these regions were structured and more limited.
Genome-wide, homozygous SNP regions were found to be unevenly distributed, with some scaffolds possessing long runs of homozygosity, see for example Algv5b02, Algv5r009x, Algv5r013x, Algv5r014x, Algv5r018x, Algv5r019x, Algv5r027x (depicted by solid blue in Figure 3B), and these regions were net negative by the Tajima D test. Conversely, heterozygous SNPs were less structured and appeared randomly distributed throughout the genome (Figure 3B). Overall, three genetic types were resolved by this analysis: in each genome, there existed SNP-poor homozygous regions (colored blue) or SNP dense regions, which either possessed homozygous alternate SNPs (also colored blue) or heterozygous SNPs (colored in ‘red’ or ‘yellow’ blocks depending on the density of heterozygous SNPs resolved in each 10 kb block: one haplotype was similar to ALV5 and the other was different). Only one worm specimen (119_3) was heterozygous genome-wide, and this track is depicted as ‘red’ across all scaffolds in the Circos plot (Figure 3B).
Population genetic structure of Kenyan Ascaris worm specimens
A phylogenetic tree constructed using genome wide SNPs with at least 10x coverage (11.15 million phased SNPs total) from 69 Ascaris worm specimens, including the A. suum reference genome, established that the Kenyan specimens were more similar to each other than they were to the A. suum reference genome, which had many more unique SNPs (Figure 4A). Notably, the nuclear genomes from the worms that possessed A. lumbricoides-like mitochondrial genomes did not clade separately, indicating that the nuclear genomes were incongruent with the mitochondrial genomes, and likely recombinant. A co-ancestry heatmap was generated among the sympatric Ascaris, and this analysis divided the genome into discrete segments and clustered samples along the diagonal based on the greatest number of shared ancestral blocks using the nearest neighbor algorithm from fineSTRUCTURE. The Ascaris genomes resolved as 13 clusters that possessed high frequency nearest-neighbor, or shared ancestry, relationships. In contrast, the A. suum reference genome and specimen 119_3 were anomalous, likely the result of their excess heterozygosity due in part to elevated numbers of unique SNPs. Notably, nine worm specimens did not coalesce into a cluster with shared ancestry. Closer examination of these specimens indicated that their phased genomes possessed limited allelic diversity and were highly recombinant (Figure 4B). This genetic mosaicism was readily resolved by fluctuating intra-scaffold genealogies established using a sliding-window neighbor-joining topology that identified regions with incongruent tree topologies. See for example the trees generated at the scaffolds ALgV5b01, ALgV5b02, and ALgV5r001. Indeed, the pairwise SNP and FST estimates for these specimens identified segments where SNP density was low, but FST was elevated with respect to neighboring segments (see block in ALgV5b02) and the most parsimonious explanation for these results is that recombination of a limited number of distinct alleles had occurred in the regions of increased FST (Figure 4B and C).
To estimate the number of supported ancestries (K) that could be resolved in the Ascaris genomes sequenced, we calculated the Dunn index, which supported 3–6 ancestral populations (Figure 4D). A gradual increase in the Dunn Index after K = 6 was observed for an ancestral population size between 2 and 15 (Figure 4D and Figure 4—figure supplement 1). We next used POPSICLE to calculate the number of clades present within each 10 kb sliding window. Local clades were represented with a different color and painted across the genome to resolve ancestry. The SNP diversity plots across the 68 specimens identified three major ‘parentage blocks’ that were resolved as belonging to ALV5 or were genetically distinct with either both haplotypes sharing the alternate parent (homozygous alternate), or were heterozygous between the two parental haplotypes for the majority of the specimens (Figure 4E, middle Circos plot. Color hues cyan, orange, aqua).
To visualize such shared ancestry across the different Ascaris specimens at chromosome resolution, a color hue representing a local genetic ‘type’ present was assigned and integrated to construct haplotype blocks across each chromosome for the ancestries present. Chromosome painting based on shared ancestry revealed a striking mosaic of large haplotype blocks of different admixed color hues, consistent with limited genetic recombination between a low number of parentage haplotypes. These admixture patterns were readily visualized by shared color blocks between different specimens across entire scaffolds including AlgV5R019X (Figure 5A) and AlgV5R027X (Figure 5B). In low complexity regions such as the left portion of contig ALgV5R019X, only three major haplotypes were resolved (Figure 5A). Strikingly, within each of the six clades resolved, all worm specimens showed a limited, mosaic fingerprint of introgressed sequence blocks indicating that recombination has shaped the population genetic structure among the Ascaris specimens sequenced. Examples of both chromosomal segregation and recombination were seen. For example, specimens 1107E_1 and 2110F_2 shared the same chromosome at ALgV5R019X, but entirely different chromosomes at ALgV5R027X, whereas specimens 107_1, 108_1 and 2110F_2 were identical except at the subtelomeric end of ALgV5R19X. In this region two admixture blocks were resolved; 107_1 and 2110F_2 remained similar to each other but 108_1 now possessed a sequence block that was shared with specimen 119_3. This extensive chimeric pattern in chromosome painting also closely resembled the genome-wide hierarchy tree (Figure 5A). The data support a model in which the specimens are genetic recombinants between A. suum and A. lumbricoides that are predominantly inbreeding.
Geographic and demographic correlates of genetic similarity
To examine genetic clustering of worms in individual human hosts, host households and villages, and study time-points, we statistically compared genetic variation within groups (such as within a village) versus between groups (such as between villages). We found significant genetic separation between worms in different villages (Table 2, Figure 6), although worms from Kenya clustered with worms from around the world based on cox-1, rather than predominantly with each other (Figure 2A). This suggests genetic diversity is present in the population of Ascaris in these Kenyan villages, which is similar to the diversity of populations of Ascaris around the world. It also suggests that a high proportion of Ascaris transmission may occur within villages in this Kenyan setting. There was no evidence from this analysis that the 13 worms collected three months after albendazole treatment were any different than the worms collected prior to albendazole treatment (Table 2).
To expand on our observations that genetically similar worms are found around the world, but that similar worms cluster within a village, based on our nuclear SNPs data, we plotted genetic distances against geographic distances. Surprisingly, we found no significant correlations between genetic and geographic distance, neither across all five studied villages nor within the two most heavily parasitized villages (Figure 6—figure supplement 1).
Discussion
In this study, we generated a high-quality reference genome from a single worm presumed to be human A. lumbricoides. Our comparative phylogenomic analyses of this new Ascaris spp. genome against existing draft genomes of A. lumbricoides and A. suum suggest that A. suum and A. lumbricoides form a genetic complex that is capable of interbreeding, which has apparently undergone a recent worldwide, multi-species Ascaris population expansion.
Our phylogenetic analysis on the complete mitochondrial genomes (from 68 worms collected from human hosts in Kenya and other available sequences) suggests that the worms collected in Kenya mirror the separation into clade A (worms from pigs in non-endemic regions and humans in endemic regions) and clade B (worms from humans and pigs from endemic and non-endemic regions) described elsewhere (Cavallero et al., 2013). It is likely that worms in both these clades are being transmitted from human to human, as pig husbandry is rare in this area of Kenya. Patterns may differ by locality, and it is possible that some of the pig-associated (A. suum-like) worms circulating in this human population in Kenya were acquired, perhaps generations ago, by humans who lived in closer proximity to pigs. It is also possible that these worms were acquired from non-human primates (Nejsum et al., 2010), or some other Ascaris host, rather than from pigs.
However, the SNPs across the whole nuclear Ascaris genome provide significantly greater power in understanding Ascaris speciation. Importantly, our nuclear genome SNP analysis suggests that the 68 Kenyan Ascaris are distributed across multiple clades in a phylogeny based on the nuclear genomes. Overall, data from our study and other studies are consistent with a pattern where hybrid genotypes in Ascaris populations were observed (Betson et al., 2014; Cavallero et al., 2013; Criscione et al., 2007; Jesudoss Chelladurai et al., 2017). Our study represents one of the most detailed accounts of mito-nuclear discordance in nematodes echoing patterns seen in another human nematode: Onchocerca volvulus (Choi et al., 2017). The data in our current study show the occurrence of distinct mitochondrial lineages that could be evidence of early stages of species differentiation. The admixture seen within the nuclear genome, however, appears to disrupt the establishment of defined molecular speciation barriers between the different Ascaris lineages. Such patterns have been recorded in other parasites, including O. volvulus (Choi et al., 2017), the blood fluke Schistosoma (Lawton et al., 2017) and the protist Leishmania (Kato et al., 2019). Each of these studies has implicated definitive hosts in the movement of parasites between otherwise isolated populations, allowing interbreeding to take place. It is most likely the historical movement of humans and their domesticated livestock that has mediated the transport of Ascaris between localities, allowing for extensive interbreeding as shown by the nuclear genomes and resulting in the discordance observed between the mitochondrial and nuclear genomes in our study.
At a more local scale, the insights into the human transmission dynamics of Ascaris showing clustering both within an individual and in villages suggest that villages are appropriate units for interventions and that people are infected with multiple eggs from a single source. These findings are in line with clustering at the village level found in Guatemala (Anderson et al., 1995) and at the sub-village level in Nepal (Criscione et al., 2010), but not in line with the lack of small-scale geographical structuring found in Denmark, Zanzibar and Uganda (Betson et al., 2011; Betson et al., 2012; Nejsum et al., 2005a). Differences could be a result of different patterns in human and livestock movement (Betson et al., 2013).
Although the current genome is, by far, the most continuous assembly for Ascaris, it is not a full chromosome assembly due largely to repetitive sequences, in particular 120 bp tandem repeat clusters and long stretches of subtelomeric repeats. Thus, it is possible that mis-assembly in some scaffolds has increased the frequency of mosaicism detected. It is for this reason that the comparative analyses on the nuclear genome was restricted to the largest 50 scaffolds, most of which are at chromosomal resolution, with only minor localized variation due to the repeat clusters. In these high confidence scaffolds, large haplotype blocks possessing either A. suum, A. lumbricoides or both parental haplotypes (heterozygous) were readily resolved indicating that the genetic mosaicism observed could not be solely attributed to genome mis-assembly. Ultimately, future studies using ultralong PacBio (Rhoads and Au, 2015) or Nanopore (Branton et al., 2009) sequencing combined with chromosome conformation capture (Hi-C) techniques (Belaghzal et al., 2017) will improve the genome to full chromosome assembly to more accurately resolve the true extent to which recombination has impacted the population genetic structure of the Ascaris species genetic complex.
The finding that A. suum and A. lumbricoides form a genetic complex has important public health implications. Reduced treatment efficacy is not currently a common issue in Ascaris infections among humans or pigs (Levecke et al., 2018; Vercruysse et al., 2011; Zuccherato et al., 2018), although low efficacy of benzimidazoles is an issue for Trichuris trichiura in humans (Diawara et al., 2009; Furtado et al., 2016; Olsen et al., 2009) and various intestinal nematodes of veterinary importance (Jaeger and Carvalho-Costa, 2017; Kaplan and Vidyashankar, 2012; Wolstenholme et al., 2004). Extensive albendazole use in either human or pig populations could lead to resistance in both populations, if cross-species infections are common and produce fertile offspring. This study suggests that research and public health interventions targeting A. lumbricoides and A. suum should be more closely integrated, and that extensive work done by the veterinary research community may be highly relevant to mass deworming campaigns that seek to improve human health.
The similarity between Ascaris from different countries and from different vertebrate hosts suggests that Ascaris infection has spread rapidly around the world, leaving little time for it to differentiate. Taken together, these finding have very important implications for parasite control and elimination efforts that only focus on mass deworming of humans for Ascaris. The ability of pig-associated worms to become endemic in human populations indicates that a one-health approach may be necessary for the control of Ascaris. The COVID-19 pandemic has highlighted the importance of one health approaches to zoonotic diseases (Global Burden of Disease, 2020); we must use a one health approach to ensure that pigs do not serve as a reservoir and potential breeding ground for drug resistance in a parasite that can sustain community transmission in humans (Webster et al., 2016).
Materials and methods
Worm collection
Request a detailed protocolWorms were expelled as part of a larger study in rural western Kenya described previously (Easton et al., 2016, Easton et al., 2017). Worms collected from study participants in five villages (Figure 6—figure supplement 2) following treatment with 400 mg albendazole were isolated, washed, labeled and stored frozen (−15ºC). The villages were near the town of Bungoma, located at N 0.57, E 34.56. Temperatures ranged from 15°C to 30°C and rainfall is 1500 mm on average. Chicken, sheep and cattle farming are common, as is subsistence agriculture and growth of sugar cane as a cash crop. The primary spoken language is Bukusu, a dialect of Luhya.
All samples were stored in Kisumu, from which they were subsequently transported to the KEMRI-CDC offices until they were shipped to the NIH (Bethesda, MD, USA) on dry ice.
DNA extraction and sequencing
Request a detailed protocolA modified DNA extraction method was developed based on Phenol/Chloroform and Qiagen methods (available on request) and used on 75 samples (Supplementary file 3). For the five germline samples, DNA was extracted from the uterus, oviduct or ovary of the worms. For the remaining samples, DNA was extracted from somatic tissue: the body wall or the intestine. Our previous work did not reveal any differences between a variety of somatic samples including the intestine and muscle (Wang et al., 2012), thus we do not expect any significant variations in the muscle and intestine genomic DNA used in this study.
Paired-End Genome Libraries – Sixty-eight A. lumbricoides DNA samples were sequenced using Illumina HiSeq 2500 (www.illumina.com) short-read paired-end sequencing. DNA was quantified by UV Spec and Picogreen. A 100 ng of DNA based on picogreen quantification was used as template for NGS library preparation using the TruSeq Nano DNA Sample library prep kit without modification. Primer-dimers in the libraries were removed by additional AMPure beads purification. Sequencing was performed to obtain a minimum genomic depth of 20X coverage for each sample.
Mate-Pair Genome Libraries – Two samples were selected for mate-pair sequencing, based on the quality of the DNA preparation. Three independent DNA isolations (corresponding to what region of the worm or what is the sample for DNA isolation) from specimen ‘119_2.3’ were combined to obtain one μg DNA input. The mate-pair libraries were generated using the Nextera Mate Pair Library Prep Kit, following the gel-free method with the only modification that M-270 Streptavidin binding beads were used instead of M-280 beads. The libraries were amplified for 15 cycles given the low DNA input going into the circularization phase. The mate-pair fragment size averaged 6 kb with a range of 2–10 kb fragments.
Assembly and annotation of A. lumbricoides reference genome
Request a detailed protocolThe A. lumbricoides germline genome assembly was constructed using the A. suum genome as a reference. Briefly, sequencing reads from a single A. lumbricoides worm (libraries #8457, #8458, and #8778) were mapped to the A. suum germline genome assembly (Wang et al., 2017) using BWA (Li and Durbin, 2009) to generate BAM and MPILEUP alignment files. The MPILEUP files were processed with a PERL script that replaced all variation sites in the reference genome with the highest allele frequencies in the A. lumbricoides sample. A. suum genomic regions that represent <5X of A. lumbricoides reads coverage were excluded from the assembly. We further polished the genome with additional Illumina sequencing reads using Pilon and its default parameters (Walker et al., 2014). The A. lumbricoides genome was annotated using the gene models built for A. suum, using the annotation transfer tool RATT (Otto et al., 2011). The protein coding regions were defined using TransDecoder (https://github.com/TransDecoder/TransDecoder/wiki; Haas and Papanicolaou, 2016). To evaluate the gene expression across all stages, we utilized previous RNAseq data from the developmental stages (Wang et al., 2012; Wang et al., 2017), re-mapped the SRA from adult males, females, L3 and L4 stages (Jex et al., 2011) to the current gene models, and quantified the expression using tophat and cufflinks. The re-mapped reads, analyzed by JMP Genomics (SAS) across all the stages and based on the principal component analyses (Figure 1B), were grouped as adult male, adult female, L1, L2, L3 (egg L3, liver L3 and lung L3), L4, carcass, muscle, intestine, embryonic (zygote1, zygote2, zygote3, zygote4, 24 hr, 46 hr, 64 hr, 96 hr, 5d, 7d), ovaries (female mitotic region, female early pachytene, female late pachytene, female diplotene and oocyte) and testis (male mitotic region, spermatogenesis, post meiotic region, seminal vesicles and spermatids). Proteome and comparative genomics analyses were done using an in-house pipeline (Karim et al., 2011). Automated annotation of proteins was done as described earlier (Cotton et al., 2017) and based on a vocabulary of nearly 290 words found in matches to various databases, including Swissprot, Gene Ontology, KOG, Pfam, and SMART, Refseq-invertebrates and a subset of the GenBank sequences containing nematode protein sequences, as well as the presence or absence of signal peptides and transmembrane domains. Signal peptide, SecretomeP, transmembrane domains, furin cleavage sites, and mucin-type glycosylation were determined with software from the Center for Biological Sequence Analysis (Technical University of Denmark, Lyngby, Denmark) (Duckert et al., 2004; Julenius et al., 2005; Sonnhammer et al., 1998). Classification of kinases was done by Kinannote (Goldberg et al., 2013). Interproscan (Jones et al., 2014) analyses were done using the standalone version 5.34. Allergenicity of proteins were predicted by Allerdictor (Dang and Lawrence, 2014), FuzzyApp (Saravanan and Lakshmi, 2014) and AllerTOP (Dimitrov et al., 2014). Genes that had blast scores < 30% of max possible score (self-blast) in other non-Ascaris nematodes with an e-value greater than 1E-05 were considered as ‘unique’. The orthologues of predicted proteome of ALV5 across the publicly available nematode genomes (Ancylostoma caninum [International Helminth Genomes Consortium, 2019], Ancylostoma ceylanicum [International Helminth Genomes Consortium, 2019; Schwarz et al., 2015], Ancylostoma duodenale [International Helminth Genomes Consortium, 2019], Ascaris lumbricoides[International Helminth Genomes Consortium, 2019], Ascaris suum [Jex et al., 2011; Wang et al., 2012; Wang et al., 2017], Brugia malayi [Ghedin et al., 2007], Caenorhabditis elegans C. elegans Sequencing [C. elegans Sequencing Consortium, 1998], Dirofilaria immitis [Godel et al., 2012], Loa loa [Desjardins et al., 2013; Tallon et al., 2014], Necator americanus [Tang et al., 2014], Onchocerca volvulus [Cotton et al., 2017], Strongyloides ratti [Nemetschke et al., 2010], Strongyloides stercoralis [Hunt et al., 2016], Toxocara canis [International Helminth Genomes Consortium, 2019; X.-Q. Zhu et al., 2015], Trichinella spiralis [Korhonen et al., 2016; Mitreva et al., 2011], Trichuris trichiura [Foth et al., 2014], Wuchereria bancrofti International Helminth Genomes Consortium, 2019; Small et al., 2016) were analyzed using OrthoFinder (Emms and Kelly, 2015). The estimated phylogenetic tree generated was graphed using FigTree v1.4.
Further manual annotation was done as required. The data were mapped into a hyperlinked Excel spreadsheet as previously described (Bennuru et al., 2011), available in Supplementary file 2.
Read mapping and SNP analysis for whole genome sequences
Request a detailed protocolThe Illumina paired-end sequence reads of the 68 Ascaris whole genomes were trimmed by removing any adapter sequences with CutAdapt v1.12 (Martin, 2011), then low-quality sequences were filtered and trimmed using the FASTX Toolkit (http://hannonlab.cshl.edu/fastx_toolkit/). Remaining reads were then ref-mapped to the A. lumbricoides genome ALV5 reference genome (described in this paper) using either Bowtie2 v2.2.9 (Langmead and Salzberg, 2012), with very sensitive, no-discordant, and no-mixed settings or using the Burrows-Wheeler Aligner (BWA, v0.7.9) (Li and Durbin, 2009) mem in default parameters and then converted into a bam file for sorted with SAMtools (Li, 2011). Sorted reads were soft-clipped and marked-duplicated using Picard-1.8.4 (http://broadinstitute.github.io/picard; Broad Institute, 2020). Single-nucleotide polymorphisms (SNPs) were obtained using SAMtools (Li, 2011) and BCFtools (Narasimhan et al., 2016) using the mpileup function and –ploidyfile features and taking chromosomal ploidies into account. SNPs were also determined using Genome Analysis Toolkit (GATK) (McKenna et al., 2010). SNPs were called by GATK Haplotype Caller with a read coverage ≥10 x, a Phredscaled SNP quality of ≥30. Mapping statistics were generated in Perl and Awk.
Ploidy determination
Request a detailed protocolThe ploidy of each specimen was calculated using AGELESS software (http://ageless.sourceforge.net/) by dividing the chromosomes into 10 kb sliding windows and averaging the coverage within each window. The windows with zero coverage were not included in any further analyses due to sequencing noise or repeat regions (Inbar et al., 2019).
Genetic diversity
Request a detailed protocolSNPs, pi (Nei and Li, 1979), (Tajima, 1989), and FST (Dunn, 1973) values were calculated using VCFtools (Danecek et al., 2011) in 10 kb sliding windows and plotted using either Circos (Krzywinski et al., 2009) or ggbio (http://bioconductor.org/packages/release/bioc/html/ggbio.html) and VariantAnnotation (http://bioconductor.org/packages/release/bioc/html/VariantAnnotation.html) R packages (v. 3.1.0, URL http://www.R-project.org). The proportions of heterozygous and homozygous SNPs were estimated in 10 kb sliding windows using custom Java scripts to generate histogram plots in Circos (Krzywinski et al., 2009). Red and blue colors indicate the presence of 90% or more heterozygous and homozygous SNPs respectively whereas yellow color was assigned otherwise.
Co-ancestry heatmap
Request a detailed protocolThe SNP data (VCF file) was first phased accurately to estimate the haplotypes using SHAPEIT (Delaneau et al., 2013) after keeping only biallelic SNPs and loci with less than 80% missing data. Co-ancestry heatmaps were generated using the linkage model of ChromoPainter (Lawson et al., 2012) and fineSTRUCTURE (http://www.paintmychromosomes.com) based on the genome-wide phased haplotype data. For fineSTRUCTURE (version 0.02) (Lawson et al., 2012), both the burn-in and Markov Chain Monte Carlo (MCMC) after the burn-in were run for 1000 iterations with default settings. Inference was performed twice at the same parameter values.
Population genetic structure
Request a detailed protocolPopulation genetic structure was constructed using POPSICLE (Shaik et al., 2018) by comparing specimens against the reference sequence ALV5 in 10 kb sliding windows with the number of cluster K = 1 to 15 and then use the Dunn index (Dunn, 1973) to calculate the optimal number of clusters. After calculating the optimal number of clusters, POPSICLE assigned each block to the existing or new clades depending on population structure of specimens and the ancestral state of each block followed by painting in Circos plot (Krzywinski et al., 2009) with color assignment based on number of clusters.
Construction of phylogenetic trees
Request a detailed protocolIn order to determine the phylogenetic relationship between samples, we selected 19005 base positions where variants were detected in a representative sample vs the reference (ALV5), and where each sample had at least 20x coverage for each locus. Using this list, the base calls for each sample were pooled together to generate a single multi-sequence fasta file.
Next, both maximum likelihood (ML) trees and bootstrap (BS) trees were generated with a final ‘best’ tree generated from the best scoring ML and BS trees using RAxML v8.2.10 (Stamatakis, 2014). The tree was visualized in FigTree v1.4.3 (http://tree.bio.ed.ac.uk/software/figtree/).
Permutational multivariate analysis of nuclear phylogeny
Request a detailed protocolSimilarity within and between worms from different villages, households, people and time-points was analyzed based on the distance matrix of the patristic distances from the phylogenetic tree described above, using permutational multivariate analysis of variance (Adonis Vegan in R). The distance matrix underlying the phylogenetic tree was analyzed in order to measure the significance and contribution of different factors to variance between samples. Each factor (village, household, host and time-point) was analyzed both separately and sequentially. The sequence chosen was ordered based on significance of each factor when tested individually. Since multiple groupings were considered using the same dataset, multiple comparison corrections were applied. Sample sizes and descriptions of each group are shown in Table 2. Similar methods were used to analyze the mitochondrial phylogeny along the same groupings.
Mitochondrial genome assembly
Request a detailed protocolWe assembled mitochondrial genomes using a de novo approach from 68 individual Ascaris genomes. For each individual, the Ascaris mitochondrial reads in the total DNA sequencing were identified by mapping the Ascaris reads to the A. suum reference mitochondrial genome (GenBank accession: NC_001327). Adaptor sequences were trimmed prior to de novo assembly. To reduce the complexity of the de novo assembly, we randomly sampled 1000x reads from each individual (the use of higher read coverage often resulted in fragmented scaffolds) and assembled these reads using the SPAdes assembler (Bankevich et al., 2012) with continuous k-mer extension from K = 21 to the maximum k-mer allowed (average extended k-mer size = 91). The assembled scaffolds were corrected with the built-in tool in SPAdes to reduce potential assembly artifacts. Next, the assembled scaffolds were aligned to the A. suum mitochondrial reference genome using BLAST, the order of the scaffolds was adjusted, and they were joined into a single scaffold. Finally, the gaps in the scaffold were filled using GapFiller (Boetzer and Pirovano, 2012) using mitochondrial reads from the same individual to generate a complete mitochondrial genome. Using the same method, we also de novo assembled another five A. suum or A. lumbricoides mitochondrion genomes from previous studies (see Supplementary file 8).
Analysis of mitochondrial genomes
Request a detailed protocolIn order to assess overall evolutionary relationships across the complete mitochondrial genomes, we aligned the genomes using Clustal W and phylogenetic trees constructed using RaxML under the conditions of the general time reversible model (GTR) as described above for the whole genome SNP alignment. Subsequent tree files were formatted in FigTree and MEGA v7. The variation in nucleotide diversity across the mitochondrial genome was measured using sliding window analyses, with a window of 300 bp and a step of 50 bp, using DNAsp v6 (Rozas et al., 2017). In order to assess the validity of potential species groupings in the ML phylogenetic tree the Birky, 2013 X4 ratio was applied to the alignment of the complete mitochondrial genomes including both samples from Kenya and published mitochondrial reference genomes from Tanzania, Uganda, China, USA, Denmark, and the UK. The X4 ratio method of species delimitation compares the ratio of mean pairwise differences between two distinct clades (K) and the mean pairwise differences within each of the clades being compared (Θ). It is considered that if K/Θ >4 this is indicative of the two clades representing two distinct species. Owing to the fact that two clades are being compared there will be two separate values of Θ, as per recommendations of Birky, 2013, the larger Θ value is used to perform the final ratio calculation as this will provide a more conservative result which ultimately will be less likely to provide a false positive result.
Due to the extensive use of mitochondrial genome data in population genetic analyses of Ascaris, several analyses were performed to identify the effect of any population level processes that may be affecting the diversity of the parasites within Kenya. Initially, diversity indices were calculated for each of the genes within the mitochondrial genome across the entire Kenyan data set as well as considering the mitochondrial genome as a whole. In order to account for the diversity within the genic regions, we removed non-coding and tRNA sequences for these analyses. To provide a genealogical perspective of population structure of the Kenya Ascaris worm specimens, we constructed the most parsimonious haplotype network based on the protein coding sequences using the TCS algorithm as implemented in PopArt (Leigh and Bryant, 2015). Further population genetic analyses were also performed to detect the occurrence of selection on the protein coding genes of the mitochondrial genome and if there were any major departures from neutrality. Standard dN/dS ratios were performed to identify the presence of positive selection where both measures equate to 1 = neutral, >1 = positive selection, <1 = purifying selection. Both Tajima’s D and Fu’s Fs were calculated to identify any substantial departure from neutrality which could be indicative of population expansion events (Supplementary file 6). All described analyses were performed using DNAsp6 (Rozas et al., 2017). As both cox-1 and nad-4 have been used in the past for epidemiological studies, single gene phylogenies were also constructed as described previously for comparison against the whole mitochondrial genome phylogeny (Figure 2—figure supplement 2).
Owing to the extensive use of the cox-1 gene for epidemiological studies the gene was extracted from the complete mitochondrial genomes of Kenya and compared to all other available Ascaris lumbricoides and Ascaris suum cox-1 sequences housed by NCBI representing populations from across the globe. Haplotype network analyses was performed to produce the parsimonious network using TCS as implemented through PopArt (Leigh and Bryant, 2015). This provided a genealogical perspective of population structure and allowed genetic connectivity between the Kenyan samples and samples from other locations to be assessed.
Data availability
Data are available under the National Center for Biological Information (NCBI) BioProject numbers; PRJNA511012 for raw sequencing data, and PRJNA515325 for the genomic assembly. Links to all genome assemblies are available at: All (https://s3.amazonaws.com/proj-bip-prod-publicread/his-omics/ALv5/Genome_Assembly/Genome_assemblies.tar.gz), De Novo (https://s3.amazonaws.com/proj-bip-prod-publicread/his-omics/ALv5/Genome_Assembly/AL-version0-genome-assembly.fasta.gz), Semi-De Novo (V1) (https://s3.amazonaws.com/proj-bip-prod-publicread/his-omics/ALv5/Genome_Assembly/AL-version1-genome-assembly.fasta.gz), V2 – (https://s3.amazonaws.com/proj-bip-prod-publicread/his-omics/ALv5/Genome_Assembly/AL-version2-genome-assembly.fasta.gz), V3 – (https://s3.amazonaws.com/proj-bip-prod-publicread/his-omics/ALv5/Genome_Assembly/AL-version3-genome-assembly.fasta.gz), V4 – (https://s3.amazonaws.com/proj-bip-prod-publicread/his-omics/ALv5/Genome_Assembly/AL-version4-genome-assembly.fasta.gz), V5 – (https://s3.amazonaws.com/proj-bip-prod-publicread/his-omics/ALv5/Genome_Assembly/AL-version5-genome-assembly.fasta.gz), Mitochondrial – (https://s3.amazonaws.com/proj-bip-prod-publicread/his-omics/ALv5/Genome_Assembly/mitochondrial_genomes.tar.gz).
-
NCBI BioProjectID PRJNA511012. 68 Ascaris lumbricoides WGS reads from Bungoma County, Kenya.
-
NCBI BioProjectID PRJNA515325. Reference quality genome assembly of Ascaris lumbricoides from the Bungoma region of Kenya.
References
-
SPAdes: a new genome assembly algorithm and its applications to single-cell sequencingJournal of Computational Biology 19:455–477.https://doi.org/10.1089/cmb.2012.0021
-
Toward the 2020 goal of soil-transmitted helminthiasis control and eliminationPLOS Neglected Tropical Diseases 12:e0006606.https://doi.org/10.1371/journal.pntd.0006606
-
A molecular epidemiological investigation of Ascaris on Unguja, Zanzibar using isoenyzme analysis, DNA barcoding and microsatellite DNA profilingTransactions of the Royal Society of Tropical Medicine and Hygiene 105:370–379.https://doi.org/10.1016/j.trstmh.2011.04.009
-
Genetic diversity of Ascaris in southwestern UgandaTransactions of the Royal Society of Tropical Medicine and Hygiene 106:75–83.https://doi.org/10.1016/j.trstmh.2011.10.011
-
BookFrom the twig tips to the deeper branchesIn: Holland C, editors. Ascaris: The Neglected Parasite. Elsevier. pp. 265–285.https://doi.org/10.1016/C2011-0-06705-4
-
Molecular epidemiology of ascariasis: a global perspective on the transmission dynamics of Ascaris in people and pigsThe Journal of Infectious Diseases 210:932–941.https://doi.org/10.1093/infdis/jiu193
-
Ascaris lumbricoides or Ascaris suum : what′s in a Name?Journal of Infectious Diseases 213:1355.2–131356.https://doi.org/10.1093/infdis/jiw037
-
Toward almost closed genomes with GapFillerGenome Biology 13:R56.https://doi.org/10.1186/gb-2012-13-6-r56
-
BookThe potential and challenges of nanopore sequencingIn: Rodgers P, editors. Nanoscience and Technology: A Collection of Reviews From Nature Journals. CO-PUBLISHED WITH World Scientific Publishing Comp. pp. 261–268.
-
Estimating the global distribution and disease burden of intestinal nematode infections: adding up the numbers--a reviewInternational Journal for Parasitology 40:1137–1144.https://doi.org/10.1016/j.ijpara.2010.04.004
-
BookMass deworming programs in middle childhood and adolescenceIn: Bundy D. A. P, de Silva N, Horton S, Jamison D. T, Patton G. C, editors. Child and Adolescent Health and Development (3rd ed.). Washington: The International Bank for Reconstruction and Development / The World Bank. pp. 1–20.
-
Phylogeographical studies of Ascaris spp. based on ribosomal and mitochondrial DNA sequencesPLOS Neglected Tropical Diseases 7:e2170.https://doi.org/10.1371/journal.pntd.0002170
-
The genome of Onchocerca volvulus, agent of river blindnessNature Microbiology 2:16216.https://doi.org/10.1038/nmicrobiol.2016.216
-
Disentangling hybridization and host colonization in parasitic roundworms of humans and pigsProceedings of the Royal Society B: Biological Sciences 274:2669–2677.https://doi.org/10.1098/rspb.2007.0877
-
Landscape genetics reveals focal transmission of a human macroparasitePLOS Neglected Tropical Diseases 4:e665.https://doi.org/10.1371/journal.pntd.0000665
-
The variant call format and VCFtoolsBioinformatics 27:2156–2158.https://doi.org/10.1093/bioinformatics/btr330
-
Genomics of Loa loa, a Wolbachia-free filarial parasite of humansNature Genetics 45:495–500.https://doi.org/10.1038/ng.2585
-
Assays to detect beta-tubulin Codon 200 polymorphism in Trichuris trichiura and Ascaris lumbricoidesPLOS Neglected Tropical Diseases 3:e397.https://doi.org/10.1371/journal.pntd.0000397
-
AllerTOP v.2--a server for in silico prediction of allergensJournal of Molecular Modeling 20:2278.https://doi.org/10.1007/s00894-014-2278-5
-
Prediction of proprotein convertase cleavage sitesProtein Engineering, Design and Selection 17:107–112.https://doi.org/10.1093/protein/gzh013
-
Hybrid Ascaris suum/Lumbricoides (Ascarididae) Infestation in a pig farmer: a rare case of zoonotic ascariasisCentral European Journal of Public Health 21:224–226.https://doi.org/10.21101/cejph.a3798
-
Emerging zoonoses: a one health challengeEClinicalMedicine 19:100300.https://doi.org/10.1016/j.eclinm.2020.100300
-
The genome of the heartworm, Dirofilaria immitis, reveals drug and vaccine targetsThe FASEB Journal 26:4650–4661.https://doi.org/10.1096/fj.12-205096
-
Animal mitochondrial DNA as a genetic marker in population and evolutionary biologyTrends in Ecology & Evolution 4:6–11.https://doi.org/10.1016/0169-5347(89)90006-2
-
Helminth infections: the great neglected tropical diseasesJournal of Clinical Investigation 118:1311–1321.https://doi.org/10.1172/JCI34261
-
The genomic basis of parasitism in the Strongyloides clade of NematodesNature Genetics 48:299–307.https://doi.org/10.1038/ng.3495
-
Comparative genomics of the major parasitic wormsNature Genetics 51:163–174.https://doi.org/10.1038/s41588-018-0262-1
-
Molecular epidemiology of Ascaris infection among pigs in IowaThe Journal of Infectious Diseases 215:131–138.https://doi.org/10.1093/infdis/jiw507
-
InterProScan 5: genome-scale protein function classificationBioinformatics 30:1236–1240.https://doi.org/10.1093/bioinformatics/btu031
-
An inconvenient truth: global worming and anthelmintic resistanceVeterinary Parasitology 186:70–78.https://doi.org/10.1016/j.vetpar.2011.11.048
-
Phylogenomic and biogeographic reconstruction of the Trichinella complexNature Communications 7:10513.https://doi.org/10.1038/ncomms10513
-
Circos: an information aesthetic for comparative genomicsGenome Research 19:1639–1645.https://doi.org/10.1101/gr.092759.109
-
Fast gapped-read alignment with bowtie 2Nature Methods 9:357–359.https://doi.org/10.1038/nmeth.1923
-
Popart : full‐feature software for haplotype network constructionMethods in Ecology and Evolution 6:1110–1116.https://doi.org/10.1111/2041-210X.12410
-
The optimal timing of post-treatment sampling for the assessment of anthelminthic drug efficacy against Ascaris infections in humansInternational Journal for Parasitology. Drugs and Drug Resistance 8:67–69.https://doi.org/10.1016/j.ijpddr.2017.12.004
-
Fndc-1 contributes to paternal mitochondria elimination in C. elegansDevelopmental Biology 454:15–20.https://doi.org/10.1016/j.ydbio.2019.06.016
-
Ascariasis in humans and pigs on small-scale farms, Maine, USA, 2010-2013Emerging Infectious Diseases 21:332–334.https://doi.org/10.3201/eid2102.140048
-
The draft genome of the parasitic nematode Trichinella spiralisNature Genetics 43:228–235.https://doi.org/10.1038/ng.769
-
Preventive chemotherapy and the fight against neglected tropical diseasesExpert Review of Anti-Infective Therapy 10:237–242.https://doi.org/10.1586/eri.11.165
-
Ascariasis is a zoonosis in DenmarkJournal of Clinical Microbiology 43:1142–1148.https://doi.org/10.1128/JCM.43.3.1142-1148.2005
-
Ascaris phylogeny based on multiple whole mtDNA genomesInfection, Genetics and Evolution 48:4–9.https://doi.org/10.1016/j.meegid.2016.12.003
-
A genetic map of the animal-parasitic nematode Strongyloides rattiMolecular and Biochemical Parasitology 169:124–127.https://doi.org/10.1016/j.molbiopara.2009.10.008
-
Albendazole and mebendazole have low efficacy against trichuristrichiura in school-age children in Kabale district, UgandaTransactions of the Royal Society of Tropical Medicine and Hygiene 103:443–446.https://doi.org/10.1016/j.trstmh.2008.12.010
-
RATT: rapid annotation transfer toolNucleic Acids Research 39:e57.https://doi.org/10.1093/nar/gkq1268
-
Ascariasis in people and pigs: new inferences from DNA analysis of worm populationsInfection, Genetics and Evolution 12:227–235.https://doi.org/10.1016/j.meegid.2012.01.012
-
PacBio sequencing and its applicationsGenomics, Proteomics & Bioinformatics 13:278–289.https://doi.org/10.1016/j.gpb.2015.08.002
-
DnaSP 6: dna sequence polymorphism analysis of large data setsMolecular Biology and Evolution 34:3299–3302.https://doi.org/10.1093/molbev/msx248
-
Fuzzy logic for personalized healthcare and diagnostics: fuzzyapp--a fuzzy logic based allergen-protein predictorOMICS: A Journal of Integrative Biology 18:570–581.https://doi.org/10.1089/omi.2014.0021
-
ConferenceA hidden markov model for predicting transmembrane helices in protein sequencesProceedings/International Conference on Intelligent Systems for Molecular Biology ; ISMB. pp. 175–182.
-
Mitochondrial versus nuclear gene sequences in deep-level mammalian phylogeny reconstructionMolecular Biology and Evolution 18:132–143.https://doi.org/10.1093/oxfordjournals.molbev.a003787
-
Gene silencing and sex determination by programmed DNA elimination in parasitic NematodesCurrent Opinion in Microbiology 32:120–127.https://doi.org/10.1016/j.mib.2016.05.012
-
Statistical method for testing the neutral mutation hypothesis by DNA polymorphismGenetics 123:585–595.
-
Experimental infection of man with Ascaris of man and the pigThe Kitasato Archives of Experimental Medicine 23:151–159.
-
Genome of the human hookworm Necator americanusNature Genetics 46:261–269.https://doi.org/10.1038/ng.2875
-
Silencing of germline-expressed genes by DNA elimination in somatic cellsDevelopmental Cell 23:1072–1080.https://doi.org/10.1016/j.devcel.2012.09.020
-
Comparative genome analysis of programmed DNA elimination in NematodesGenome Research 27:2001–2014.https://doi.org/10.1101/gr.225730.117
-
Programmed DNA elimination in multicellular organismsCurrent Opinion in Genetics & Development 27:26–34.https://doi.org/10.1016/j.gde.2014.03.012
-
One health - an ecological and evolutionary framework for tackling neglected zoonotic diseasesEvolutionary Applications 9:313–333.https://doi.org/10.1111/eva.12341
-
Discordant mitochondrial and nuclear gene phylogenies in emydid turtles: implications for speciation and conservationBiological Journal of the Linnean Society 99:445–461.https://doi.org/10.1111/j.1095-8312.2009.01342.x
-
Drug resistance in veterinary helminthsTrends in Parasitology 20:469–476.https://doi.org/10.1016/j.pt.2004.07.010
-
Phylogenetic performance of mitochondrial protein-coding genes in resolving relationships among vertebratesMolecular Biology and Evolution 13:933–942.https://doi.org/10.1093/oxfordjournals.molbev.a025661
-
Phylogeography of Ascaris lumbricoides and A. suum from ChinaParasitology Research 109:329–338.https://doi.org/10.1007/s00436-011-2260-4
-
Characterisation of Ascaris from human and pig hosts by nuclear ribosomal DNA sequencesInternational Journal for Parasitology 29:469–478.https://doi.org/10.1016/s0020-7519(98)00226-4
-
Genetic blueprint of the zoonotic pathogen Toxocara canisNature Communications 6:6145.https://doi.org/10.1038/ncomms7145
Article and author information
Author details
Funding
National Institutes of Health (AI114054)
- Richard E Davis
National Institutes of Health (AI125869)
- Jianbin Wang
Bill and Melinda Gates Foundation
- Roy Anderson
National Institutes of Health
- Thomas B Nutman
Wellcome Trust (KEMRI)
- Roy Anderson
London Centre for Neglected Tropical Disease Research
- Roy Anderson
The funders had no role in study design, data collection and interpretation, or the decision to submit the work for publication.
Acknowledgements
We thank the school children, schoolteachers, and Bungoma administrators for their support. We extend special thanks to all the members of the study team: Bungoma County Hospital, Siangwe, Siaka, Sang’alo, Nasimbo and Ranje village administrators and Community Health Workers. Particular thanks to Dr. Charles S Mwandawiro, Prof. Sammy Njenga, and Dr. Jimmy H Kihara (KEMRI), and Dr Simon J Brooker (BMGF) for making the fieldwork possible in Kenya, and for their invaluable scientific and logistical advice.
Ethics
Human subjects: This study was approved by the Ethics Review Committee of the Kenya Medical Research Institute (Scientific Steering Committee protocol number 2688) and the Imperial College Research Ethics Committee (ICREC_ 13_1_15). Informed written consent was obtained from all adults and parents or guardians of each child. Minor assent was obtained from all children aged 12-17. Anyone found to be infected with any STH was treated with 400 mg ALB during each phase of the study, and all previously-untreated village residents were offered ALB at the end of each study phase.
Copyright
This is an open-access article, free of all copyright, and may be freely reproduced, distributed, transmitted, modified, built upon, or otherwise used by anyone for any lawful purpose. The work is made available under the Creative Commons CC0 public domain dedication.
Metrics
-
- 2,426
- views
-
- 295
- downloads
-
- 45
- citations
Views, downloads and citations are aggregated across all versions of this paper published by eLife.
Download links
Downloads (link to download the article as PDF)
Open citations (links to open the citations from this article in various online reference manager services)
Cite this article (links to download the citations from this article in formats compatible with various reference manager tools)
Further reading
-
- Epidemiology and Global Health
Artificially sweetened beverages containing noncaloric monosaccharides were suggested as healthier alternatives to sugar-sweetened beverages. Nevertheless, the potential detrimental effects of these noncaloric monosaccharides on blood vessel function remain inadequately understood. We have established a zebrafish model that exhibits significant excessive angiogenesis induced by high glucose, resembling the hyperangiogenic characteristics observed in proliferative diabetic retinopathy (PDR). Utilizing this model, we observed that glucose and noncaloric monosaccharides could induce excessive formation of blood vessels, especially intersegmental vessels (ISVs). The excessively branched vessels were observed to be formed by ectopic activation of quiescent endothelial cells (ECs) into tip cells. Single-cell transcriptomic sequencing analysis of the ECs in the embryos exposed to high glucose revealed an augmented ratio of capillary ECs, proliferating ECs, and a series of upregulated proangiogenic genes. Further analysis and experiments validated that reduced foxo1a mediated the excessive angiogenesis induced by monosaccharides via upregulating the expression of marcksl1a. This study has provided new evidence showing the negative effects of noncaloric monosaccharides on the vascular system and the underlying mechanisms.
-
- Epidemiology and Global Health
- Microbiology and Infectious Disease
Influenza viruses continually evolve new antigenic variants, through mutations in epitopes of their major surface proteins, hemagglutinin (HA) and neuraminidase (NA). Antigenic drift potentiates the reinfection of previously infected individuals, but the contribution of this process to variability in annual epidemics is not well understood. Here, we link influenza A(H3N2) virus evolution to regional epidemic dynamics in the United States during 1997—2019. We integrate phenotypic measures of HA antigenic drift and sequence-based measures of HA and NA fitness to infer antigenic and genetic distances between viruses circulating in successive seasons. We estimate the magnitude, severity, timing, transmission rate, age-specific patterns, and subtype dominance of each regional outbreak and find that genetic distance based on broad sets of epitope sites is the strongest evolutionary predictor of A(H3N2) virus epidemiology. Increased HA and NA epitope distance between seasons correlates with larger, more intense epidemics, higher transmission, greater A(H3N2) subtype dominance, and a greater proportion of cases in adults relative to children, consistent with increased population susceptibility. Based on random forest models, A(H1N1) incidence impacts A(H3N2) epidemics to a greater extent than viral evolution, suggesting that subtype interference is a major driver of influenza A virus infection ynamics, presumably via heterosubtypic cross-immunity.