Applying evolutionary genetics to schistosome epidemiology

Population genomics of helminth parasites

Next generation sequencing technologies have facilitated a shift from a few targeted loci in population genetic studies to whole genome approaches. Here, we review the types of questions and inferences regarding the population biology and evolution of parasitic helminths being addressed within the field of population genomics. Topics include parabiome, hybridization, population structure, loci under selection and linkage mapping. We highlight various advances, and note the current trends in the field, particularly a focus on human-related parasites despite the inherent biodiversity of helminth species. We conclude by advocating for a broader application of population genomics to reflect the taxonomic and life history breadth displayed by helminth parasites. As such, our basic knowledge about helminth population biology and evolution would be enhanced while the diversity of helminths in itself would facilitate population genomic comparative studies to address broader ecological and evolutionary concepts.

Prospects for genomic surveillance for selection in schistosome parasites

Schistosomiasis is a neglected tropical disease caused by multiple parasitic Schistosoma species, and which impacts over 200 million people globally, mainly in low- and middle-income countries. Genomic surveillance to detect evidence for natural selection in schistosome populations represents an emerging and promising approach to identify and interpret schistosome responses to ongoing control efforts or other environmental factors. Here we review how genomic variation is used to detect selection, how these approaches have been applied to schistosomes, and how future studies to detect selection may be improved. We discuss the theory of genomic analyses to detect selection, identify experimental designs for such analyses, and review studies that have applied these approaches to schistosomes. We then consider the biological characteristics of schistosomes that are expected to respond to selection, particularly those that may be impacted by control programs. Examples include drug resistance, host specificity, and life history traits, and we review our current understanding of specific genes that underlie them in schistosomes. We also discuss how inherent features of schistosome reproduction and demography pose substantial challenges for effective identification of these traits and their genomic bases. We conclude by discussing how genomic surveillance for selection should be designed to improve understanding of schistosome biology, and how the parasite changes in response to selection.

Integrating genomic and epidemiologic data to accelerate progress toward schistosomiasis elimination

The global community has adopted ambitious goals to eliminate schistosomiasis as a public health problem, and new tools are needed to achieve them. Mass drug administration programs, for example, have reduced the burden of schistosomiasis, but the identification of hotspots of persistent and reemergent transmission threaten progress toward elimination and underscore the need to couple treatment with interventions that reduce transmission. Recent advances in DNA sequencing technologies make whole-genome sequencing a valuable and increasingly feasible option for population-based studies of complex parasites such as schistosomes. Here, we focus on leveraging genomic data to tailor interventions to distinct social and ecological circumstances. We consider two priority questions that can be addressed by integrating epidemiological, ecological, and genomic information: (1) how often do non-human host species contribute to human schistosome infection? and (2) what is the importance of locally acquired versus imported infections in driving transmission at different stages of elimination? These questions address processes that can undermine control programs, especially those that rely heavily on treatment with praziquantel. Until recently, these questions were difficult to answer with sufficient precision to inform public health decision-making. We review the literature related to these questions and discuss how whole-genome approaches can identify the geographic and taxonomic sources of infection, and how such information can inform context-specific efforts that advance schistosomiasis control efforts and minimize the risk of reemergence.

Genotyping of Ascaris spp. infecting humans and pigs in Italy, Slovakia and Colombia

BACKGROUND

The systematics and taxonomy of Ascaris lumbricoides and Ascaris suum, two of the world's most widespread nematodes, still represent a highly debated scientific issue. Two different transmission scenarios have been described according to endemicity: separated host-specific transmission cycles in endemic regions, and a single pool of infection shared by humans and pigs in non-endemic regions. The swine roundworm A. suum is now recognized as an important cause of human ascariasis also in endemic areas such as China, where cross-infections and hybridization have also been reported, as well as in non-endemic regions like Italy. This study aimed to investigate the molecular epidemiology of human and pig ascariasis in three countries representing different epidemiological scenarios: Italy as a non-endemic country, Colombia as an endemic country, and Slovakia as a non-endemic country, but with a poor socio-economic context linked to some focal populations of Roma settlements.

MATERIALS AND METHODS

A total of 237 nematodes were analysed: 46 from Colombia (13 from humans, 33 from pigs), 114 from Slovakia (20 from humans, 94 from pigs) and 77 from Italy (17 from humans and 60 from pigs). Genotyping by PCR-RFLP of nuclear (ITS) and sequencing of mitochondrial (cox1) target regions were performed. ITS genotypes were used to estimate the Hardy-Weinberg (HW) equilibrium according to hosts and country of origin. The partial cox1 sequences were used to analyse genetic polymorphisms according to hosts and country of origin, as well as to infer the network of haplotypes, their evolutionary relationships and geographical distribution.

RESULTS

110 quality cox1 sequences were obtained. Haplotype network revealed three main groups corresponding to clade A, B and C. Clade C included most of the human cases from Italy, while those from Slovakia and Colombia were grouped in clade B. Ascaris from Italian and Colombian pigs showed HW equilibrium at the ITS marker, while disequilibrium was found in A. lumbricoides from Slovak pigs, which suggest a high unexpected amount of roundworms of human origin circulating also in pigs.

CONCLUSIONS

This study updates and extends the current understanding of Ascaris species and genotypes circulating in different epidemiological scenarios, with particular attention to the inclusion of human-derived Ascaris in the phylogenetic cluster C. Despite the evidence of HW equilibrium in the ITS in pig-derived Italian samples, the amount of genetic variation seems to support the existence of two closely related species.

Natural Intra- and Interclade Human Hybrid Schistosomes in Africa with Considerations on Prevention through Vaccination

Causal agents of schistosomiasis are dioecious, digenean schistosomes affecting mankind in 76 countries. Preventive measures are manifold but need to be complemented by vaccination for long-term protection; vaccine candidates in advanced pre-clinical/clinical stages include Sm14, Sm-TSP-2/Sm-TSP-2Al^®, Smp80/SchistoShield^®, and Sh28GST/Bilhvax^®. Natural and anthropogenic changes impact on breaking species isolation barriers favoring introgressive hybridization, i.e., allelic exchange among gene pools of sympatric, interbreeding species leading to instant large genetic diversity. Phylogenetic distance matters, thus the less species differ phylogenetically the more likely they hybridize. PubMed and Embase databases were searched for publications limited to hybridale confirmation by mitochondrial cytochrome c oxidase (COX) and/or nuclear ribosomal internal transcribed spacer (ITS). Human schistosomal hybrids are predominantly reported from West Africa with clustering in the Senegal River Basin, and scattering to Europe, Central and Eastern Africa. Noteworthy is the dominance of Schistosoma haematobium interbreeding with human and veterinary species leading due to hybrid vigor to extinction and homogenization as seen for S. guineensis in Cameroon and S. haematobium in Niger, respectively. Heterosis seems to advantage S. haematobium/S. bovis interbreeds with dominant S. haematobium-ITS/S. bovis-COX1 profile to spread from West to East Africa and reoccur in France. S. haematobium/S. mansoni interactions seen among Senegalese and Côte d’Ivoirian children are unexpected due to their high phylogenetic distance. Detecting pure S. bovis and S. bovis/S. curassoni crosses capable of infecting humans observed in Corsica and Côte d’Ivoire, and Niger, respectively, is worrisome. Taken together, species hybridization urges control and preventive measures targeting human and veterinary sectors in line with the One-Health concept to be complemented by vaccination protecting against transmission, infection, and disease recurrence. Functional and structural diversity of naturally occurring human schistosomal hybrids may impact current vaccine candidates requiring further research including natural history studies in endemic areas targeted for clinical trials.

Investigating the genetic structure of the parasites Anisakis pegreffii and A. berlandi (Nematoda: Anisakidae) in a sympatric area of the southern Pacific Ocean waters using a multilocus genotyping approach: first evidence of their interspecific hybridization

The southern Pacific Ocean, off the New Zealand coast, has been reported as one sympatric area of the two parasite species Anisakis pegreffii and A. berlandi. Here, a multilocus genotyping approach, based on a panel of eleven DNA microsatellite (SSR) loci plus the sequences analysis of the nuclear nas10 nDNA and the mitochondrial mtDNA cox2 gene loci, was applied to a total of N = 344 adults and larvae of Anisakis spp. from cetacean and fish species, respectively. Out of the newly scored SSR loci, Anisl 15 and Anisl 2 showed fixed alternative alleles between A. pegreffii and A. berlandi resulting as 100% diagnostic loci. Out of SSRs Anisl 00314 and Anisl 7 previously disclosed, two additional loci, i.e., Anisl 4 and Anisl 22, were found to be sex-linked. The Bayesian genotypes clustering approach (STRUCTURE) allowed identification, with a 100% of probability value, N = 208 specimens to the "pure parental" A. pegreffii, N = 133 to the "pure parental" A. berlandi, while one adult and two larval stages showed mixed ancestry between the two groups having, in all cases, a Q-value = 0.50. NEWHYBRIDS analysis assigned (100% of probability) those specimens to their F1 hybrid category. This represents the first evidence of contemporary hybridization between the two parasite species in a sympatric area. The pairwise F_ST values estimated at intraspecific and interspecific level, inferred from both SSR loci and mitochondrial mtDNA cox2 sequences, have also demonstrated the existence of two distinct panmictic units in this study area, corresponding respectively to A. pegreffii and A. berlandi. The results obtained support the useful application of a multilocus approach in the identification of sibling species and their hybrid categories in sympatric areas. The possible use of sex-linked SSR loci of the two species of the A. simplex (s. l.), for sex determination of their larval stages, is also suggested.

Population genomic analyses of schistosome parasites highlight critical challenges facing endgame elimination efforts

Schistosomiasis persists in Asian regions despite aggressive elimination measures. To identify factors enabling continued parasite transmission, we performed reduced representation genome sequencing on Schistosoma japonicum miracidia collected across multiple years from transmission hotspots in Sichuan, China. We discovered strong geographic structure, suggesting that local, rather than imported, reservoirs are key sources of persistent infections in the region. At the village level, parasites collected after referral for praziquantel treatment are closely related to local pre-treatment populations. Schistosomes within villages are also highly related, suggesting that only a few parasites from a limited number of hosts drive re-infection. The close familial relationships among miracidia from different human hosts also implicate short transmission routes among humans. At the individual host level, genetic evidence indicates that multiple humans retained infections following referral for treatment. Our findings suggest that end-game schistosomiasis control measures should focus on completely extirpating local parasite reservoirs and confirming successful treatment of infected human hosts.

Schistosomiasis Morbidity Hotspots: Roles of the Human Host, the Parasite and Their Interface in the Development of Severe Morbidity

Schistosomiasis is the second most important human parasitic disease in terms of socioeconomic impact, causing great morbidity and mortality, predominantly across the African continent. For intestinal schistosomiasis, severe morbidity manifests as periportal fibrosis (PPF) in which large tracts of macro-fibrosis of the liver, visible by ultrasound, can occlude the main portal vein leading to portal hypertension (PHT), sequelae such as ascites and collateral vasculature, and ultimately fatalities. For urogenital schistosomiasis, severe morbidity manifests as pathology throughout the urinary system and genitals, and is a definitive cause of squamous cell bladder carcinoma. Preventative chemotherapy (PC) programmes, delivered through mass drug administration (MDA) of praziquantel (PZQ), have been at the forefront of schistosomiasis control programmes in sub-Saharan Africa since their commencement in Uganda in 2003. However, despite many successes, 'biological hotspots' (as distinct from 'operational hotspots') of both persistent high transmission and morbidity remain. In some areas, this failure to gain control of schistosomiasis has devastating consequences, with not only persistently high infection intensities, but both "subtle" and severe morbidity remaining prevalent. These hotspots highlight the requirement to revisit research into severe morbidity and its mechanisms, a topic that has been out of favor during times of PC implementation. Indeed, the focality and spatially-structured epidemiology of schistosomiasis, its transmission persistence and the morbidity induced, has long suggested that gene-environmental-interactions playing out at the host-parasite interface are crucial. Here we review evidence of potential unique parasite factors, host factors, and their gene-environmental interactions in terms of explaining differential morbidity profiles in the human host. We then take the situation of schistosomiasis mansoni within the Albertine region of Uganda as a case study in terms of elucidating the factors behind the severe morbidity observed and the avenues and directions for future research currently underway within a new research and clinical trial programme (FibroScHot).

Prevalence of Schistosoma mansoni Infection and Associated Risk Factors among School Children in Jiga Town, Northwest-Ethiopia: A Cross-Sectional Study

BACKGROUND

Schistosomiasis, one of the 20 Neglected Tropical Diseases, is a human disease which is caused by infection with Schistosoma spp. In Ethiopia, the intestinal form of schistosomiasis (caused by Schistosoma mansoni) is widely distributed. This study is aimed at assessing the prevalence of S. mansoni and associated risk factors among schoolchildren in Jiga town and its surrounding rural "Kebeles", Northwest Ethiopia (Kebele = the smallest administrative unit in Ethiopia).

METHODS

A cross-sectional study was carried out among children from three schools in Jiga town and its surrounding rural "Kebeles". The schools were selected mainly according to their proximity to water bodies that were likely to have been used by the children. A total of 362 randomly selected children were screened for S. mansoni infection. Samples of stool were collected from each pupil and examined by the concentration and microscopic techniques. Semistructured interview questionnaires were administered to the study subjects to identify the possible risk factors of infection with S. mansoni. Data analysis was performed using SPSS version 20 software.

RESULTS

The overall prevalence of S. mansoni was 15.2%. A higher proportion of S. mansoni infection (30.3%) was observed at Chifarit primary school children followed by Tikur Wuha primary school children (17.3%) and Jiga 01 primary school children (10.2%). The habit of contact with river water (AOR, 4.86; 95% CI, 0.991-23.8; P = 0.051) and the irregular wearing of shoes (AOR, 6.06; 95% CI, 2.61-14.13; P = 0.001) were risk factors for S. mansoni infection.

CONCLUSION

S. mansoni infection remains a public health problem in the study area. Snail control, mass deworming, and health promotion are recommended.

A multi-restriction fragment length polymorphism genotyping approach including the beta-tubulin gene as a new differential nuclear marker for the recognition of the cryptic species Anisakis simplex s.s. and Anisakis pegreffii and their hybridization events

The detection of Anisakis simplex s.s./A. pegreffii putative hybrids has been a controversial issue in spite of the fact that natural hybridization is an extended process across free living and parasitic organisms. Differential traits of biomedical and ecological importance, such as invasive and allergenic potential have been demonstrated in both cryptic species. Therefore, in this work, we discuss about the potential for hybridization between these anisakid species in sympatric zones, implementing a multi-marker Restriction fragment length polymorphism (RFLP) genotyping approach based on the ribosomal DNA internal transcribed spacer 1 (ITS1), the mitochondrial cytochrome C oxidase 2 (Cox-2) and a new nuclear marker, the highly conserved β-tubulin gene (β-TUB). The two cryptic species differed at least in 7 bp in the β-TUB gene and some larvae with heterozygous genotypes at the 7 diagnostic nucleotide positions were found. Taxonomic, population and genealogical analyses served to support the occurrence of hybridization between both species. Predicted restriction endonucleases enzymes were assayed for Cox-2 and β-TUB markers. The implemented multi-marker PCR-RFLP allowed us to detect the two pure parental species, F1 hybrids, hybrid backcrossed progeny and individuals with nuclear-mitochondrial discordance, being a useful, simple and reproducible procedure in any laboratory for epidemiological studies.

Impacts of host gender on Schistosoma mansoni risk in rural Uganda-A mixed-methods approach

Background

The World Health Organization identified Uganda as one of the 10 highly endemic countries for schistosomiasis. Annual mass drug administration (MDA) with praziquantel has led to a decline in intensity of Schistosoma mansoni infections in several areas. However, as hotspots with high (re)infection rates remain, additional research on risk factors and implementing interventions to complement MDA are required to further reduce disease burden in these settings. Through a mixed-methods study we aimed to gain deeper understanding of how gender may impact risk and reinfection in order to inform disease control programmes and ascertain if gender-specific interventions may be beneficial.

Methodology/Principal findings

In Bugoto, Mayuge District, Eastern Uganda we conducted ethnographic observations (n = 16) and examined epidemiology (n = 55) and parasite population genetics (n = 16) in school-aged children (SAC), alongside a community-wide household survey (n = 130). Water contact was frequent at home, school and in the community and was of domestic, personal care, recreational, religious or commercial nature. Qualitative analysis of type of activity, duration, frequency, level of submersion and water contact sites in children showed only few behavioural differences in water contact between genders. However, survey data revealed that adult women carried out the vast majority of household tasks involving water contact. Reinfection rates (96% overall) and genetic diversity were high in boys (pre-He = 0.66; post-He = 0.67) and girls (pre-He = 0.65; post-He = 0.67), but no differences in reinfection rates (p = 0.62) or genetic diversity by gender before (p = 0.54) or after (p = 0.97) treatment were found.

Conclusions/Significance

This mixed methods approach showed complementary findings. Frequent water exposure with few differences between boys and girls was mirrored by high reinfection rates and genetic diversity in both genders. Disease control programmes should consider the high reinfection rates among SAC in remaining hotspots of schistosomiasis and the various purposes and settings in which children and adults are exposed to water.

Globally, over 230 million people are infected with schistosomiasis, an infectious disease caused by parasitic helminths. Humans can get infected when they contact water which contains Schistosoma parasites. Although the disease can be treated with a drug, people get rapidly reinfected in certain high-transmission settings. Drug treatment alone may not be sufficient to eliminate this disease and additional interventions such as health promotion or improvements in water and sanitation need to be scaled up. To provide recommendations to these control programmes we carried out interdisciplinary research in Eastern Uganda to understand the influence of gender on schistosomiasis risk. We found that the water contact behaviour of boys and girls is quite similar, and we did not see differences in reinfection or genetic diversity of the parasite between boys and girls. Differences in water contact between genders is greater in adults, and further research is required for these individuals. In this setting, infection rates are high in school-aged children and there are no differences between genders. These results emphasise improved control efforts for all school-aged children in communities like these. Our interdisciplinary approach provided complementary findings. Such an integrated approach can therefore have more power to meaningfully inform policy on schistosomiasis control.

Two-year longitudinal survey reveals high genetic diversity of Schistosoma mansoni with adult worms surviving praziquantel treatment at the start of mass drug administration in Uganda

BACKGROUND

A key component of schistosomiasis control is mass drug administration with praziquantel. While control interventions have been successful in several endemic regions, mass drug administration has been less effective in others. Here we focus on the impact of repeated praziquantel treatment on the population structure and genetic diversity of Schistosoma mansoni.

METHODS

We examined S. mansoni epidemiology, population genetics, and variation in praziquantel susceptibility in parasites isolated from children across three primary schools in a high endemicity region at the onset of the Ugandan National Control Programme. Children were sampled at 11 timepoints over two years, including one week and four weeks post-praziquantel treatment to evaluate short-term impacts on clearance and evidence of natural variation in susceptibility to praziquantel.

RESULTS

Prevalence of S. mansoni was 85% at baseline. A total of 3576 miracidia larval parasites, isolated from 203 individual children, were genotyped at seven loci. Overall, genetic diversity was high and there was low genetic differentiation, indicating high rates of parasite gene flow. Schistosome siblings were found both pre-treatment and four weeks post-treatment, demonstrating adult worms surviving treatment and natural praziquantel susceptibility variation in these populations at the beginning of mass drug administration. However, we did not find evidence for selection on these parasites. While genetic diversity decreased in the short-term (four weeks post-treatment), diversity did not decrease over the entire period despite four rounds of mass treatment. Furthermore, within-host genetic diversity was affected by host age, host sex, infection intensity and recent praziquantel treatment.

CONCLUSIONS

Our findings suggest that praziquantel treatments have short-term impacts on these parasite populations but impacts were transient and no long-term reduction in genetic diversity was observed. High gene flow reduces the likelihood of local adaptation, so even though parasites surviving treatment were observed, these were likely to be diluted at the beginning of the Ugandan National Control Programme. Together, these results suggest that MDA in isolation may be insufficient to reduce schistosome populations in regions with high genetic diversity and gene flow.

Development of EST-derived microsatellite markers to investigate the population structure of sparganum - the causative agent of zoonotic sparganosis

The plerocercoid (sparganum) of Spirometra erinaceieuropaei is the main aetiological agent of human sparganosis. To improve the current knowledge on S. erinaceieuropaei evolution, we performed multi-locus microsatellite typing of sparganum isolates from China for the first time. All available expressed sequence tag (EST) sequences for the Spirometra were downloaded from the GenBank. The identification and localization of microsatellites in ESTs was accomplished by MISA. Based on the selected microsatellites, the genetic structure of 64 sparganum isolates collected from 11 geographical locations in southwest China were investigated through principal component analysis, STRUCTURE analysis and neighbour-joining clustering. A total of 522 non-redundant ESTs containing 915 simple sequence repeats were identified from 12 481 ESTs screened. Five primer pairs were finally selected. Using these loci, a total of 12 alleles were detected in 64 sparganum isolates. Little variability was observed within each of geographical population, especially among isolates derived from Kunming of Yunnan (YN-KM) province. Both STRUCTURE analysis and the clustering analysis supported that two genotypes existed among the sparganum isolates from southwest China. In conclusion, five microsatellite markers were successfully developed, and sparganum population was observed to harbour low genetic variation, further investigation with deeper sampling was needed to elucidate the population structure.

Whole genome amplification and exome sequencing of archived schistosome miracidia

Adult schistosomes live in the blood vessels and cannot easily be sampled from humans, so archived miracidia larvae hatched from eggs expelled in feces or urine are commonly used for population genetic studies. Large collections of archived miracidia on FTA cards are now available through the Schistosomiasis Collection at the Natural History Museum (SCAN). Here we describe protocols for whole genome amplification of Schistosoma mansoni and Schistosome haematobium miracidia from these cards, as well as real time PCR quantification of amplified schistosome DNA. We used microgram quantities of DNA obtained for exome capture and sequencing of single miracidia, generating dense polymorphism data across the exome. These methods will facilitate the transition from population genetics, using limited numbers of markers to population genomics using genome-wide marker information, maximising the value of collections such as SCAN.

Assessing S. mansoni prevalence in Biomphalaria snails in the Gombe ecosystem of western Tanzania: the importance of DNA sequence data for clarifying species identification

Background

Snails are essential for the transmission and maintenance of schistosomiasis in endemic areas, as they serve as intermediate hosts for schistosome parasites. A clear understanding of the snail species present, their local distribution and infection status is therefore a prerequisite for effective control of schistosomiasis. The purpose of this study was to establish the infection status and distribution of Schistosoma mansoni in snails in the Gombe area along the shores of Lake Tanganyika in western Tanzania, using both detection of cercarial shedding and molecular approaches.

Methods

Snails were collected from streams located close to human settlements in Gombe National Park, as well as from nearby villages (Kiziba, Mtanga, Mwamgongo and Bugamba) and the largest town in the region (Kigoma). Snails were individually exposed to light to induce shedding of schistosome larvae, which were examined using a compound light microscope. Additionally, the internal transcribed spacer (ITS) region of the ribosomal RNA gene cluster was simultaneously amplified in both snails and their trematodes using a single polymerase chain reaction (PCR) and sequenced to confirm species identification.

Results

Snails morphologically identified as Biomphalaria pfeifferi were present in all streams except at Mtanga but their distribution was patchy in both time and space. Sequencing of PCR products indicated that not all snails were B. pfeifferi. None of the snails from Gombe or Bugamba shed schistosome larvae, while larvae were shed at all other sites. Overall, an infection prevalence of only 12% was observed in snails based on cercarial shedding. While 47% of the snails were PCR-positive for the 500 bp ITS fragment, which was predicted to indicate infection with S. mansoni, sequence data demonstrated that these bands are not species-specific and can be amplified from other trematode infections. In addition, a 1000 bp band was amplified in 14% of samples, which was identified as a trematode in the family Derogenidae.

Conclusions

The results support the previous assumption that B. pfeifferi snails may be involved in transmitting schistosomiasis in the area but suggest that the community structure of both snails and trematodes may be more complicated than previously thought. This emphasises the importance of confirming species identifications using sequencing, rather than relying only on PCR-based diagnostics or cercarial shedding.

Electronic supplementary material

The online version of this article (10.1186/s13071-017-2525-6) contains supplementary material, which is available to authorized users.

No more time to stay 'single' in the detection of Anisakis pegreffii, A. simplex (s. s.) and hybridization events between them: a multi-marker nuclear genotyping approach

A multi-marker nuclear genotyping approach was performed on larval and adult specimens of Anisakis spp. (N = 689) collected from fish and cetaceans in allopatric and sympatric areas of the two species Anisakis pegreffii and Anisakis simplex (s. s.), in order to: (1) identify specimens belonging to the parental taxa by using nuclear markers (allozymes loci) and sequence analysis of a new diagnostic nuclear DNA locus (i.e. partial sequence of the EF1 α-1 nDNA region) and (2) recognize hybrid categories. According to the Bayesian clustering algorithms, based on those markers, most of the individuals (N = 678) were identified as the parental species [i.e. A. pegreffii or A. simplex (s. s.)], whereas a smaller portion (N = 11) were recognized as F1 hybrids. Discordant results were obtained when using the polymerase chain reaction-restriction fragment length polymorphisms (PCR-RFLPs) of the internal transcribed spacer (ITS) ribosomal DNA (rDNA) on the same specimens, which indicated the occurrence of a large number of 'hybrids' both in sympatry and allopatry. These findings raise the question of possible misidentification of specimens belonging to the two parental Anisakis and their hybrid categories derived from the application of that single marker (i.e. PCR-RFLPs analysis of the ITS of rDNA). Finally, Bayesian clustering, using allozymes and EF1 α-1 nDNA markers, has demonstrated that hybridization between A. pegreffii and A. simplex (s. s.) is a contemporary phenomenon in sympatric areas, while no introgressive hybridization takes place between the two species.

Whole Genome Amplification and Reduced-Representation Genome Sequencing of Schistosoma japonicum Miracidia

Background

In areas where schistosomiasis control programs have been implemented, morbidity and prevalence have been greatly reduced. However, to sustain these reductions and move towards interruption of transmission, new tools for disease surveillance are needed. Genomic methods have the potential to help trace the sources of new infections, and allow us to monitor drug resistance. Large-scale genotyping efforts for schistosome species have been hindered by cost, limited numbers of established target loci, and the small amount of DNA obtained from miracidia, the life stage most readily acquired from humans. Here, we present a method using next generation sequencing to provide high-resolution genomic data from S. japonicum for population-based studies.

Methodology/Principal Findings

We applied whole genome amplification followed by double digest restriction site associated DNA sequencing (ddRADseq) to individual S. japonicum miracidia preserved on Whatman FTA cards. We found that we could effectively and consistently survey hundreds of thousands of variants from 10,000 to 30,000 loci from archived miracidia as old as six years. An analysis of variation from eight miracidia obtained from three hosts in two villages in Sichuan showed clear population structuring by village and host even within this limited sample.

Conclusions/Significance

This high-resolution sequencing approach yields three orders of magnitude more information than microsatellite genotyping methods that have been employed over the last decade, creating the potential to answer detailed questions about the sources of human infections and to monitor drug resistance. Costs per sample range from $50-$200, depending on the amount of sequence information desired, and we expect these costs can be reduced further given continued reductions in sequencing costs, improvement of protocols, and parallelization. This approach provides new promise for using modern genome-scale sampling to S. japonicum surveillance, and could be applied to other schistosome species and other parasitic helminthes.

Schistosomiasis is a devastating tropical disease that affects more than 200 million people worldwide. Over the past several decades, transmission control strategies implemented in China have reduced the prevalence and morbidity of Schistosoma japonicum in many areas. Infections still persist, however, and it is therefore of great interest to determine the sources of recurring infections. Genetic analysis is a promising means to achieve this. Towards this aim, we conducted a pilot study to assess the feasibility of using high-throughput sequencing to assess the geographic distribution of schistosome genetic variants. Because DNA yields from miracidia, the most easily accessible life stage, are insufficient for high throughput sequencing, we first employed whole genome amplification to obtain sufficient quantities of DNA. We then employed a technique that reproducibly sequences the same fraction of a genome across numerous samples. We successfully sequenced 6-year old S. japonicum samples from Sichuan Province, China, easily and economically identifying tens of thousands of variable loci, a sufficient number to discriminate fine-scale population structure. Further population sampling will help answer important questions concerning the persistence of infections, the sources of new infections, and whether parasite populations have undergone incipient evolution of drug resistance.

Population genomics of the filarial nematode parasite Wuchereria bancrofti from mosquitoes

Wuchereria bancrofti is a parasitic nematode and the primary cause of lymphatic filariasis--a disease specific to humans. W. bancrofti currently infects over 90 million people throughout the tropics and has been acknowledged by the world health organization as a vulnerable parasite. Current research has focused primarily on the clinical manifestations of disease and little is known about the evolutionary history of W. bancrofti. To improve upon knowledge of the evolutionary history of W. bancrofti, we whole genome sequenced 13 W. bancrofti larvae. We circumvent many of the difficulties of multiple infections by sampling larvae directly from mosquitoes that were experimentally inoculated with infected blood. To begin, we used whole genome data to reconstruct the historical population size. Our results support a history of fluctuating population sizes that can be correlated with human migration and fluctuating mosquito abundances. Next, we reconstructed the putative pedigree of W. bancrofti worms within an infection using the kinship coefficient. We deduced that there are full-sib and half-sib relationships residing within the same larval cohort. Through combined analysis of the mitochondrial and nuclear genomes we concluded that this is likely a results of polyandrous mating, the first time reported for W. bancrofti. Lastly, we scanned the genomes for signatures of natural selection. Annotation of putative selected regions identified proteins that may have aided in a parasitic life style or may have evolved to protect against current drug treatments. We discuss our results in the greater context of understanding the biology of an animal with a unique life history and ecology.

Characterisation of a novel panel of polymorphic microsatellite loci for the liver fluke, Fasciola hepatica, using a next generation sequencing approach

•

2448 microsatellite loci were identified within 83 Mb of F. hepatica genome sequence data.

•

A panel of 15 polymorphic loci were developed and validated using genomic DNA from 46 parasites.

•

The panel was developed and optimised as a multiplex PCR protocol.

•

All loci could be amplified from several F. hepatica lifecycle stages with the multiplex approach.

The liver fluke, Fasciola hepatica is an economically important pathogen of sheep and cattle and has been described by the WHO as a re-emerging zoonosis. Control is heavily reliant on the use of drugs, particularly triclabendazole and as a result resistance has now emerged. The population structure of F. hepatica is not well known, yet it can impact on host–parasite interactions and parasite control with drugs, particularly regarding the spread of triclabendazole resistance. We have identified 2448 potential microsatellites from 83 Mb of F. hepatica genome sequence using msatfinder. Thirty-five loci were developed and optimised for microsatellite PCR, resulting in a panel of 15 polymorphic loci, with a range of three to 15 alleles. This panel was validated on genomic DNA from 46 adult F. hepatica; 38 liver flukes sourced from a Northwest abattoir, UK and 8 liver flukes from an established isolate (Shrewsbury; Ridgeway Research). Evidence for null alleles was found at four loci (Fh_1, Fh_8, Fh_13 and Fh_14), which showed markedly higher levels of homozygosity than the remaining 11 loci. Of the 38 liver flukes isolated from cattle livers (n = 10) at the abattoir, 37 genotypes were identified. Using a multiplex approach all 15 loci could be amplified from several life cycle stages that typically yield low amounts of DNA, including metacercariae, the infective life cycle stage present on pasture, highlighting the utility of this multiplex microsatellite panel. This study reports the largest panel of microsatellite markers available to date for population studies of F. hepatica and the first multiplex panel of microsatellite markers that can be used for several life cycle stages.

Modelling the effects of mass drug administration on the molecular epidemiology of schistosomes

As national governments scale up mass drug administration (MDA) programs aimed to combat neglected tropical diseases (NTDs), novel selection pressures on these parasites increase. To understand how parasite populations are affected by MDA and how to maximize the success of control programmes, it is imperative for epidemiological, molecular and mathematical modelling approaches to be combined. Modelling of parasite population genetic and genomic structure, particularly of the NTDs, has been limited through the availability of only a few molecular markers to date. The landscape of infectious disease research is being dramatically reshaped by next-generation sequencing technologies and our understanding of how repeated selective pressures are shaping parasite populations is radically altering. Genomics can provide high-resolution data on parasite population structure, and identify how loci may contribute to key phenotypes such as virulence and/or drug resistance. We discuss the incorporation of genetic and genomic data, focussing on the recently sequenced Schistosoma spp., into novel mathematical transmission models to inform our understanding of the impact of MDA and other control methods. We summarize what is known to date, the models that exist and how population genetics has given us an understanding of the effects of MDA on the parasites. We consider how genetic and genomic data have the potential to shape future research, highlighting key areas where data are lacking, and how future molecular epidemiology knowledge can aid understanding of transmission dynamics and the effects of MDA, ultimately informing public health policy makers of the best interventions for NTDs.

Molecular epidemiology, phylogeny and evolution of the filarial nematode Wuchereria bancrofti

Wuchereria bancrofti (Wb) is the most widely distributed of the three nematodes known to cause lymphatic filariasis (LF), the other two being Brugia malayi and Brugia timori. Current tools available to monitor LF are limited to diagnostic tests targeting DNA repeats, filarial antigens, and anti-filarial antibodies. While these tools are useful for detection and surveillance, elimination programs have yet to take full advantage of molecular typing for inferring infection history, strain fingerprinting, and evolution. To date, molecular typing approaches have included whole mitochondrial genomes, genotyping, targeted sequencing, and random amplified polymorphic DNA (RAPDs). These studies have revealed much about Wb biology. For example, in one study in Papua New Guinea researchers identified 5 major strains that were widespread and many minor strains some of which exhibit geographic stratification. Genome data, while rare, has been utilized to reconstruct evolutionary relationships among taxa of the Onchocercidae (the clade of filarial nematodes) and identify gene synteny. Their phylogeny reveals that speciation from the common ancestor of both B. malayi and Wb occurred around 5-6 millions years ago with shared ancestry to other filarial nematodes as recent as 15 million years ago. These discoveries hold promise for gene discovery and identifying drug targets in species that are more amenable to in vivo experiments. Continued technological developments in whole genome sequencing and data analysis will likely replace many other forms of molecular typing, multiplying the amount of data available on population structure, genetic diversity, and phylogenetics. Once widely available, the addition of population genetic data from genomic studies should hasten the elimination of LF parasites like Wb. Infectious disease control programs have benefited greatly from population genetics data and recently from population genomics data. However, while there is currently a surplus of data for diseases like malaria and HIV, there is a scarcity of this data for filarial nematodes. With the falling cost of genome sequencing, research on filarial nematodes could benefit from the addition of population genetics statistics and phylogenetics especially in dealing with elimination programs. A comprehensive review focusing on population genetics of filarial nematode does not yet exist. Here our goal is to provide a current overview of the molecular epidemiology of W. bancrofti (Wb) the primary causative agent of LF. We begin by reviewing studies utilizing molecular typing techniques with specific focus on genomic and population datasets. Next, we used whole mitochondrial genome data to construct a phylogeny and examine the evolutionary history of the Onchocercidae. Then, we provide a perspective to aid in understanding how population genetic techniques translate to modern epidemiology. Finally, we introduce the concept of genomic epidemiology and provide some examples that will aid in future studies of Wb.

Inbreeding within human Schistosoma mansoni: do host-specific factors shape the genetic composition of parasite populations?

The size, structure and distribution of host populations are key determinants of the genetic composition of parasite populations. Despite the evolutionary and epidemiological merits, there has been little consideration of how host heterogeneities affect the evolutionary trajectories of parasite populations. We assessed the genetic composition of natural populations of the parasite Schistosoma mansoni in northern Senegal. A total of 1346 parasites were collected from 14 snail and 57 human hosts within three villages and individually genotyped using nine microsatellite markers. Human host demographic parameters (age, gender and village of residence) and co-infection with Schistosoma haematobium were documented, and S. mansoni infection intensities were quantified. F-statistics and clustering analyses revealed a random distribution (panmixia) of parasite genetic variation among villages and hosts, confirming the concept of human hosts as 'genetic mixing bowls' for schistosomes. Host gender and village of residence did not show any association with parasite genetics. Host age, however, was significantly correlated with parasite inbreeding and heterozygosity, with children being more infected by related parasites than adults. The patterns may be explained by (1) genotype-dependent 'concomitant immunity' that leads to selective recruitment of genetically unrelated worms with host age, and/or (2) the 'genetic mixing bowl' hypothesis, where older hosts have been exposed to a wider variety of parasite strains than children. The present study suggests that host-specific factors may shape the genetic composition of schistosome populations, revealing important insights into host-parasite interactions within a natural system.

Population genetic structure of Schistosoma mansoni and Schistosoma haematobium from across six sub-Saharan African countries: implications for epidemiology, evolution and control

We conducted the first meta-analysis of ten Schistosoma haematobium (one published and nine unpublished) and eight Schistosoma mansoni (two published and six unpublished) microsatellite datasets collected from individual schistosome-infected school-children across six sub-Saharan Africa countries. High levels of genetic diversity were documented in both S. haematobium and S. mansoni. In S. haematobium populations, allelic richness did not differ significantly between the ten schools, despite widely varying prevalences and intensities of infection, but higher levels of heterozygote deficiency were seen in East than in West Africa. In contrast, S. mansoni populations were more diverse in East than West African schools, but heterozygosity levels did not vary significantly with geography. Genetic structure in both S. haematobium and S. mansoni populations was documented, at both a regional and continental scale. Such structuring might be expected to slow the spread to new areas of anti-schistosomal drug resistance should it develop. There was, however, limited evidence of genetic structure at the individual host level, which might be predicted to promote the development or establishment of drug resistance, particularly if it were a recessive trait. Our results are discussed in terms of their potential implications for the epidemiology and evolution of schistosomes as well as their subsequent control across sub-Saharan Africa.

Reductions in genetic diversity of Schistosoma mansoni populations under chemotherapeutic pressure: the effect of sampling approach and parasite population definition

Detecting potential changes in genetic diversity in schistosome populations following chemotherapy with praziquantel (PZQ) is crucial if we are to fully understand the impact of such chemotherapy with respect to the potential emergence of resistance and/or other evolutionary outcomes of interventions. Doing so by implementing effective, and cost-efficient sampling protocols will help to optimise time and financial resources, particularly relevant to a disease such as schistosomiasis currently reliant on a single available drug. Here we explore the effect on measures of parasite genetic diversity of applying various field sampling approaches, both in terms of the number of (human) hosts sampled and the number of transmission stages (miracidia) sampled per host for a Schistosoma mansoni population in Tanzania pre- and post-treatment with PZQ. In addition, we explore population structuring within and between hosts by comparing the estimates of genetic diversity obtained assuming a 'component population' approach with those using an 'infrapopulation' approach. We found that increasing the number of hosts sampled, rather than the number of miracidia per host, gives more robust estimates of genetic diversity. We also found statistically significant population structuring (using Wright's F-statistics) and significant differences in the measures of genetic diversity depending on the parasite population definition. The relative advantages, disadvantages and, hence, subsequent reliability of these metrics for parasites with complex life-cycles are discussed, both for the specific epidemiological and ecological scenario under study here and for their future application to other areas and schistosome species.

Significant variance in genetic diversity among populations of Schistosoma haematobium detected using microsatellite DNA loci from a genome-wide database

BACKGROUND

Urogenital schistosomiasis caused by Schistosoma haematobium is widely distributed across Africa and is increasingly being targeted for control. Genome sequences and population genetic parameters can give insight into the potential for population- or species-level drug resistance. Microsatellite DNA loci are genetic markers in wide use by Schistosoma researchers, but there are few primers available for S. haematobium.

METHODS

We sequenced 1,058,114 random DNA fragments from clonal cercariae collected from a snail infected with a single Schistosoma haematobium miracidium. We assembled and aligned the S. haematobium sequences to the genomes of S. mansoni and S. japonicum, identifying microsatellite DNA loci across all three species and designing primers to amplify the loci in S. haematobium. To validate our primers, we screened 32 randomly selected primer pairs with population samples of S. haematobium.

RESULTS

We designed >13,790 primer pairs to amplify unique microsatellite loci in S. haematobium, (available at http://www.cebio.org/projetos/schistosoma-haematobium-genome). The three Schistosoma genomes contained similar overall frequencies of microsatellites, but the frequency and length distributions of specific motifs differed among species. We identified 15 primer pairs that amplified consistently and were easily scored. We genotyped these 15 loci in S. haematobium individuals from six locations: Zanzibar had the highest levels of diversity; Malawi, Mauritius, Nigeria, and Senegal were nearly as diverse; but the sample from South Africa was much less diverse.

CONCLUSIONS

About half of the primers in the database of Schistosoma haematobium microsatellite DNA loci should yield amplifiable and easily scored polymorphic markers, thus providing thousands of potential markers. Sequence conservation among S. haematobium, S. japonicum, and S. mansoni is relatively high, thus it should now be possible to identify markers that are universal among Schistosoma species (i.e., using DNA sequences conserved among species), as well as other markers that are specific to species or species-groups (i.e., using DNA sequences that differ among species). Full genome-sequencing of additional species and specimens of S. haematobium, S. japonicum, and S. mansoni is desirable to better characterize differences within and among these species, to develop additional genetic markers, and to examine genes as well as conserved non-coding elements associated with drug resistance.

Non-invasive sampling of schistosomes from humans requires correcting for family structure

For ethical and logistical reasons, population-genetic studies of parasites often rely on the non-invasive sampling of offspring shed from their definitive hosts. However, if the sampled offspring are naturally derived from a small number of parents, then the strong family structure can result in biased population-level estimates of genetic parameters, particularly if reproductive output is skewed. Here, we document and correct for the strong family structure present within schistosome offspring (miracidia) that were collected non-invasively from humans in western Kenya. By genotyping 2,424 miracidia from 12 patients at 12 microsatellite loci and using a sibship clustering program, we found that the samples contained large numbers of siblings. Furthermore, reproductive success of the breeding schistosomes was skewed, creating differential representation of each family in the offspring pool. After removing the family structure with an iterative jacknifing procedure, we demonstrated that the presence of relatives led to inflated estimates of genetic differentiation and linkage disequilibrium, and downwardly-biased estimates of inbreeding coefficients (F_IS). For example, correcting for family structure yielded estimates of F_ST among patients that were 27 times lower than estimates from the uncorrected samples. These biased estimates would cause one to draw false conclusions regarding these parameters in the adult population. We also found from our analyses that estimates of the number of full sibling families and other genetic parameters of samples of miracidia were highly intercorrelated but are not correlated with estimates of worm burden obtained via egg counting (Kato-Katz). Whether genetic methods or the traditional Kato-Katz estimator provide a better estimate of actual number of adult worms remains to be seen. This study illustrates that family structure must be explicitly accounted for when using offspring samples to estimate the genetic parameters of adult parasite populations.

Genetic epidemiology uses genetic data to uncover patterns of disease processes. To acquire data for these analyses, individual pathogens are collected and scored at genetic markers, and the resultant data are analyzed to infer biological patterns about the pathogen populations. In lieu of invasive sampling of adult pathogens in humans, researchers have relied on non-invasive sampling of parasite offspring (often shed in fecal samples). One potential problem with this approach is that analyses using the offspring data will be biased because many of the offspring are related and family sizes are likely to be unequal. We show that this sampling issue is relevant in a natural transmission zone in western Kenya and that it yields biases in three important parameters: genetic differentiation, inbreeding coefficients, and estimates of the amount of non-random association between loci (linkage disequilibrium). We also develop a method to remove these biases by removing the sibling structure present in the dataset. Finally, we suggest that our measure of family number, as well as other genetic measures, may be useful measures of the worm burdens in patients.

Approaches to genotyping individual miracidia of Schistosoma japonicum

Molecular genetic tools are needed to address questions as to the source and dynamics of transmission of the human blood fluke Schistosoma japonicum in regions where human infections have reemerged, and to characterize infrapopulations in individual hosts. The life stage that interests us as a target for collecting genotypic data is the miracidium, a very small larval stage that consequently yields very little DNA for analysis. Here, we report the successful development of a multiplex format permitting genotyping of 17 microsatellite loci in four sequential multiplex reactions using a single miracidium held on a Whatman Classic FTA indicating card. This approach was successful after short storage periods, but after long storage (>4 years), considerable difficulty was encountered in multiplex genotyping, necessitating the use of whole genome amplification (WGA) methods. WGA applied to cards stored for long periods of time resulted in sufficient DNA for accurate and repeatable genotyping. Trials and tests of these methods, as well as application to some field-collected samples, are reported, along with the discussion of the potential insights to be gained from such techniques. These include recognition of sibships among miracidia from a single host, and inference of the minimum number of worm pairs that might be present in a host.

Unravelling parasitic nematode natural history using population genetics

The health and economic importance of parasitic nematodes cannot be overstated. Moreover, they offer a complex and diverse array of life strategies, raising a multitude of evolutionary questions. Researchers are applying population genetics to parasitic nematodes in order to disentangle some aspects of their life strategies, improve our knowledge about disease epidemiology, and design control strategies. However, population genetics studies of nematodes have been constrained due to the difficulty in sampling nematodes and developing molecular markers. In this context, new computational and sequencing technologies represent promising tools to investigate population genomics of parasitic, non-model, nematode species in an epidemiological context.

Elucidating the temporal and spatial dynamics of Biomphalaria glabrata genetic diversity in three Brazilian villages

OBJECTIVE

The freshwater snail Biomphalaria glabrata is the principal intermediate host for the parasite Schistosoma mansoni within Brazil. We assessed the potential effects of snail population dynamics on parasite transmission dynamics via population genetics.

METHODS

We sampled snail populations located within the confines of three schistosome-endemic villages in the state of Minas Gerais, Brazil. Snails were collected from individual microhabitats following seasonal periods of flood and drought over the span of 1 year. Snail spatio-temporal genetic diversity and population differentiation of 598 snails from 12 sites were assessed at seven microsatellite loci.

RESULTS

Average genetic diversity was relatively low, ranging from 4.29 to 9.43 alleles per locus, and overall, subpopulations tended to exhibit heterozygote deficits. Genetic diversity was highly spatially partitioned among subpopulations, while virtually, no partitioning was observed across temporal sampling. Comparison with previously published parasite genetic diversity data indicated that S. mansoni populations are significantly more variable and less subdivided than those of the B. glabrata intermediate hosts.

DISCUSSION

Within individual Brazilian villages, observed distributions of snail genetic diversity indicate temporal stability and very restricted gene flow. This is contrary to observations of schistosome genetic diversity over the same spatial scale, corroborating the expectation that parasite gene flow at the level of individual villages is likely driven by vertebrate host movement.

Introgressive hybridization of Schistosoma haematobium group species in Senegal: species barrier break down between ruminant and human schistosomes

Background

Schistosomes are dioecious parasitic flatworms, which live in the vasculature of their mammalian definitive hosts. They are the causative agent of schistosomiasis, a disease of considerable medical and veterinary importance in tropical and subtropical regions. Schistosomes undergo a sexual reproductive stage within their mammalian host enabling interactions between different species, which may result in hybridization if the species involved are phylogenetically close. In Senegal, three closely related species in the Schistosoma haematobium group are endemic: S. haematobium, which causes urogenital schistosomiasis in humans, and S. bovis and S. curassoni, which cause intestinal schistosomiasis in cows, sheep and goats.

Methodology/Principal Findings

Large-scale multi-loci molecular analysis of parasite samples collected from children and domestic livestock across Senegal revealed that interactions and hybridization were taking place between all three species. Evidence of hybridization between S. haematobium/S. curassoni and S. haematobium/S. bovis was commonly found in children from across Senegal, with 88% of the children surveyed in areas of suspected species overlap excreting hybrid miracidia. No S. haematobium worms or hybrids thereof were found in ruminants, although S. bovis and S. curassoni hybrid worms were found in cows. Complementary experimental mixed species infections in laboratory rodents confirmed that males and females of each species readily pair and produce viable hybrid offspring.

Conclusions/Significance

These data provide indisputable evidence for: the high occurrence of bidirectional hybridization between these Schistosoma species; the first conclusive evidence for the natural hybridisation between S. haematobium and S. curassoni; and demonstrate that the transmission of the different species and their hybrids appears focal. Hybridization between schistosomes has been known to influence the disease epidemiology and enhance phenotypic characteristics affecting transmission, morbidity and drug sensitivity. Therefore, understanding and monitoring such inter-species interactions will be essential for optimizing and evaluating control strategies across such potential hybrid zones.

Schistosome blood flukes are transmitted through water contact and cause a severe debilitating disease in humans and their livestock. Understanding the biology and epidemiology of these parasites is essential to enable the development of better control strategies in endemic areas. Several species of schistosome exist and species barriers are normally maintained by differences in ecology, host specificity, and evolutionary history. However, hybridization between closely related species can occur if parasites infect the same definitive host. Here we report on the introgression between a human (S. haematobium) and two ruminant schistosomes (S. bovis and S. curassoni), the prevalence and distribution of their hybrids and the novel evidence for the presence of S. haematobium/S. curassoni hybrids in Senegalese children. Modern sampling and genotyping techniques have brought to light the extent of these hybridization events which could have been facilitated by the natural progression in farming, population (both human and livestock) movements and expansion, as well as changes in snail distribution, creating areas of sympatric transmission. Hybridization can lead to phenotypic characteristics that can influence disease epidemiology and control success, highlighting the importance of monitoring these evolving hybrid zones.

Genetics: A New Landscape for Medical Geography

The emergence and re-emergence of human pathogens resistant to medical treatment will present a challenge to the international public health community in the coming decades. Geography is uniquely positioned to examine the progressive evolution of pathogens across space and through time, and to link molecular change to interactions between population and environmental drivers. Landscape as an organizing principle for the integration of natural and cultural forces has a long history in geography, and, more specifically, in medical geography. Here, we explore the role of landscape in medical geography, the emergent field of landscape genetics, and the great potential that exists in the combination of these two disciplines. We argue that landscape genetics can enhance medical geographic studies of local-level disease environments with quantitative tests of how human-environment interactions influence pathogenic characteristics. In turn, such analyses can expand theories of disease diffusion to the molecular scale and distinguish the important factors in ecologies of disease that drive genetic change of pathogens.

New frontiers in schistosoma genomics and transcriptomics

Schistosomes are digenean blood flukes of aves and mammals comprising 23 species. Some species are causative agents of human schistosomiasis, the second major neglected disease affecting over 230 million people worldwide. Modern technologies including the sequencing and characterization of nucleic acids and proteins have allowed large-scale analyses of parasites and hosts, opening new frontiers in biological research with potential biomedical and biotechnological applications. Nuclear genomes of the three most socioeconomically important species (S. haematobium, S. japonicum, and S. mansoni) have been sequenced and are under intense investigation. Mitochondrial genomes of six Schistosoma species have also been completely sequenced and analysed from an evolutionary perspective. Furthermore, DNA barcoding of mitochondrial sequences is used for biodiversity assessment of schistosomes. Despite the efforts in the characterization of Schistosoma genomes and transcriptomes, many questions regarding the biology and evolution of this important taxon remain unanswered. This paper aims to discuss some advances in the schistosome research with emphasis on genomics and transcriptomics. It also aims to discuss the main challenges of the current research and to point out some future directions in schistosome studies.

Intestinal schistosomiasis in chimpanzees on Ngamba Island, Uganda: observations on liver fibrosis, schistosome genetic diversity and praziquantel treatment

Despite treatment with praziquantel (PZQ) at 40 mg/kg in food, several chimpanzees on Ngamba Island Chimpanzee Sanctuary (NICS) continue to excrete eggs of Schistosoma mansoni. To monitor disease, 8 animals were closely examined under anaesthesia in March 2011 with portable ultrasonography and by rectal snip biopsy. Schistosome genetic diversity had been previously assayed within 4 of these chimpanzees, finding extensive diversity with 27 DNA barcodes encountered, although none was common to all animals. Calcified schistosome eggs were found in the rectal snips from 5 chimpanzees and liver fibrosis was clearly documented, indicative of progressive disease in 6 animals, the latter being surprisingly advanced in a younger chimpanzee. All 8 animals were treated under anaesthesia by oral gavage with PZQ at 60 mg/kg dosing that was well tolerated. These animals were again re-examined in June 2012 using stool and urine sampling. Only 1 chimpanzee appeared to be free from infection and active egg excretion was confirmed in 6 animals. If intestinal schistosomiasis is to be controlled within this setting, a long-term disease management plan is required which should combine active case-detection with an insistent treatment regime with praziquantel for these chimpanzees, exploring perhaps the performance of even higher dosing.

Testing local-scale panmixia provides insights into the cryptic ecology, evolution, and epidemiology of metazoan animal parasites

When every individual has an equal chance of mating with other individuals, the population is classified as panmictic. Amongst metazoan parasites of animals, local-scale panmixia can be disrupted due to not only non-random mating, but also non-random transmission among individual hosts of a single host population or non-random transmission among sympatric host species. Population genetics theory and analyses can be used to test the null hypothesis of panmixia and thus, allow one to draw inferences about parasite population dynamics that are difficult to observe directly. We provide an outline that addresses 3 tiered questions when testing parasite panmixia on local scales: is there greater than 1 parasite population/species, is there genetic subdivision amongst infrapopulations within a host population, and is there asexual reproduction or a non-random mating system? In this review, we highlight the evolutionary significance of non-panmixia on local scales and the genetic patterns that have been used to identify the different factors that may cause or explain deviations from panmixia on a local scale. We also discuss how tests of local-scale panmixia can provide a means to infer parasite population dynamics and epidemiology of medically relevant parasites.

Genetic diversity of Schistosoma japonicum miracidia from individual rodent hosts

Schistosoma japonicum is an important parasite in terms of clinical, veterinary and socio-economic impacts, and rodents, a long neglected reservoir for the parasite, have recently been found to act as reservoir hosts in some endemic areas of China. Any difference in the host's biological characteristics and/or associated living habitats among rodents may result in different environments for parasites, possibly resulting in a specific population structure of parasites within hosts. Therefore knowledge of the genetic structure of parasites within individual rodents could improve our understanding of transmission dynamics and hence our ability to develop effective control strategies. In this study, we aimed to describe a host-specific structure for S. japonicum and its potential influencing factors. The results showed a significant genetic differentiation among hosts. Two factors, including sampling seasons and the number of miracidia genotyped per host, showed an effect on the genetic diversity of an infrapopulation through a univariable analysis but not a multivariable analysis. A possible scenario of clustered infection foci and the fact of multiple definitive host species, the latter of which is unique to S. japonicum compared with other schistosomes, were proposed to explain the observed results and practical implications for control strategies are recommended.

An infectious topic in reticulate evolution: introgression and hybridization in animal parasites

Little attention has been given to the role that introgression and hybridization have played in the evolution of parasites. Most studies are host-centric and ask if the hybrid of a free-living species is more or less susceptible to parasite infection. Here we focus on what is known about how introgression and hybridization have influenced the evolution of protozoan and helminth parasites of animals. There are reports of genome or gene introgression from distantly related taxa into apicomplexans and filarial nematodes. Most common are genetic based reports of potential hybridization among congeneric taxa, but in several cases, more work is needed to definitively conclude current hybridization. In the medically important Trypanosoma it is clear that some clonal lineages are the product of past hybridization events. Similarly, strong evidence exists for current hybridization in human helminths such as Schistosoma and Ascaris. There remain topics that warrant further examination such as the potential hybrid origin of polyploid platyhelminths. Furthermore, little work has investigated the phenotype or fitness, and even less the epidemiological significance of hybrid parasites.