Genomic linkage map of the human blood fluke Schistosoma mansoni

Genomic data reveal a north-south split and introgression history of blood fluke (Schistosoma haematobium) populations from across Africa

The human parasitic fluke, Schistosoma haematobium hybridizes with the livestock parasite S. bovis in the laboratory, but the extent of hybridization in nature is unclear. We analyzed 34.6 million single nucleotide variants in 162 samples from 18 African countries, revealing a sharp genetic discontinuity between northern and southern S. haematobium. We found no evidence for recent hybridization. Instead the data reveal admixture events that occurred 257-879 generations ago in northern S. haematobium populations. Fifteen introgressed S. bovis genes are approaching fixation in northern S. haematobium with four genes potentially driving adaptation. We identified 19 regions that were resistant to introgression; these were enriched on the sex chromosomes. These results (i) demonstrate strong barriers to gene flow between these species, (ii) indicate that hybridization may be less common than currently envisaged, but (iii) reveal profound genomic consequences of interspecific hybridization between schistosomes of medical and veterinary importance.

Schistosomiasis

Schistosomiasis is a major cause of morbidity in the world and almost 800 million people worldwide are at risk for schistosomiasis; it is second only to malaria as a major infectious disease. Globally, it is estimated that the disease affects more than 250 million people in 78 countries of the world and is responsible for some 280,000-500,000 deaths each year. The three major schistosomes infecting humans are Schistosoma mansoni, S. japonicum, and S. haematobium. This chapter covers a wide range of aspects of schistosomiasis, including basic biology of the parasites, epidemiology, immunopathology, treatment, control, vaccines, and genomics/proteomics. In this chapter, the reader will understand the significant toll this disease takes in terms of mortality and morbidity. A description of the various life stages of schistosomes is presented, which will be informative for both those unfamiliar with the disease and experienced scientists. Clinical and public health aspects are addressed that cover acute and chronic disease, diagnosis, current treatment regimens and alternative drugs, and schistosomiasis control programs. A brief overview of genomics and proteomics is included that details recent advances in the field that will help scientists investigate the molecular biology of schistosomes. The reader will take away an appreciation for general aspects of schistosomiasis and the current research advances.

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.

Genome-wide analysis of Schistosoma mansoni reveals limited population structure and possible praziquantel drug selection pressure within Ugandan hot-spot communities

Populations within schistosomiasis control areas, especially those in Africa, are recommended to receive regular mass drug administration (MDA) with praziquantel (PZQ) as the main strategy for controlling the disease. The impact of PZQ treatment on schistosome genetics remains poorly understood, and is limited by a lack of high-resolution genetic data on the population structure of parasites within these control areas. We generated whole-genome sequence data from 174 individual miracidia collected from both children and adults from fishing communities on islands in Lake Victoria in Uganda that had received either annual or quarterly MDA with PZQ over four years, including samples collected immediately before and four weeks after treatment. Genome variation within and between samples was characterised and we investigated genomic signatures of natural selection acting on these populations that could be due to PZQ treatment. The parasite population on these islands was more diverse than found in nearby villages on the lake shore. We saw little or no genetic differentiation between villages, or between the groups of villages with different treatment intensity, but slightly higher genetic diversity within the pre-treatment compared to post-treatment parasite populations. We identified classes of genes significantly enriched within regions of the genome with evidence of recent positive selection among post-treatment and intensively treated parasite populations. The differential selection observed in post-treatment and pre-treatment parasite populations could be linked to any reduced susceptibility of parasites to praziquantel treatment.

Genomic analysis of a parasite invasion: Colonization of the Americas by the blood fluke Schistosoma mansoni

Schistosoma mansoni, a snail-borne, blood fluke that infects humans, was introduced into the Americas from Africa during the Trans-Atlantic slave trade. As this parasite shows strong specificity to the snail intermediate host, we expected that adaptation to South American Biomphalaria spp. snails would result in population bottlenecks and strong signatures of selection. We scored 475,081 single nucleotide variants in 143 S. mansoni from the Americas (Brazil, Guadeloupe and Puerto Rico) and Africa (Cameroon, Niger, Senegal, Tanzania, and Uganda), and used these data to ask: (i) Was there a population bottleneck during colonization? (ii) Can we identify signatures of selection associated with colonization? (iii) What were the source populations for colonizing parasites? We found a 2.4- to 2.9-fold reduction in diversity and much slower decay in linkage disequilibrium (LD) in parasites from East to West Africa. However, we observed similar nuclear diversity and LD in West Africa and Brazil, suggesting no strong bottlenecks and limited barriers to colonization. We identified five genome regions showing selection in the Americas, compared with three in West Africa and none in East Africa, which we speculate may reflect adaptation during colonization. Finally, we infer that unsampled populations from central African regions between Benin and Angola, with contributions from Niger, are probably the major source(s) for Brazilian S. mansoni. The absence of a bottleneck suggests that this is a rare case of a serendipitous invasion, where S. mansoni parasites were pre-adapted to the Americas and able to establish with relative ease.

A Schistosoma mansoni tri- and tetramer microsatellite catalog for genetic population diversity and differentiation

All Schistosoma mansoni tri- and tetranucleotide repeat microsatellites published as of December 2018 were identified. All 52 were evaluated for autosomal location, strength of amplification, scorability and behavior as single-copy loci by polyacrylamide and capillary gel electrophoresis. Of these, 27 were unique, autosomal, polymorphic, easily scored and single copy as assessed on pooled adult worm DNA from two different continental origins and adult worm clones. These microsatellites were distributed across all seven autosomal chromosomes. On laboratory strains their heterozygosity ranged from 0.22 to 0.77. Individual markers had 5-13 alleles, allelic richness of 2-10 and an effective allele number of 1.3-8.14. Those infected by Schistosoma mansoni carry many genetically distinct, sexually reproducing parasites, therefore, for an individual infection the complete allele frequency profile of their progeny consists of a pool of DNA from multiple diploid eggs. Using a set of 25 microsatellites, we calculated allele frequency profiles of eggs in fecal samples from people in two Brazilian communities separated by 6 km: Jenipapo (n = 80) and Volta do Rio (n = 38). There were no a priori characteristics that could predict the performance of markers in natural infections based on their performance with laboratory strains. Increasing marker number did not change accuracy for differentiation and diversity but did improve precision. Our data suggest that using a random set of 10-20 microsatellites appears to result in values that exhibit low standard deviations for diversity and differentiation indices. All identified microsatellites as well as PCR conditions, allele size, primer sequences and references for all tri- and tetramer microsatellites markers presented in this work are available at: https://sites.google.com/case.edu/cwru-and-fiocruz-wdrc/home.

Whole-genome sequencing of Schistosoma mansoni reveals extensive diversity with limited selection despite mass drug administration

Control and elimination of the parasitic disease schistosomiasis relies on mass administration of praziquantel. Whilst these programmes reduce infection prevalence and intensity, their impact on parasite transmission and evolution is poorly understood. Here we examine the genomic impact of repeated mass drug administration on Schistosoma mansoni populations with documented reduced praziquantel efficacy. We sequenced whole-genomes of 198 S. mansoni larvae from 34 Ugandan children from regions with contrasting praziquantel exposure. Parasites infecting children from Lake Victoria, a transmission hotspot, form a diverse panmictic population. A single round of treatment did not reduce this diversity with no apparent population contraction caused by long-term praziquantel use. We find evidence of positive selection acting on members of gene families previously implicated in praziquantel action, but detect no high frequency functionally impactful variants. As efforts to eliminate schistosomiasis intensify, our study provides a foundation for genomic surveillance of this major human parasite.

Schistosome W-linked genes inform temporal dynamics of sex chromosome evolution and suggest candidate for sex determination

Schistosomes, the human parasites responsible for snail fever, are female-heterogametic. Different parts of their ZW sex chromosomes have stopped recombining in distinct lineages, creating "evolutionary strata" of various ages. While the Z-chromosome is well characterized at the genomic and molecular level, the W-chromosome has remained largely unstudied from an evolutionary perspective, as only a few W-linked genes have been detected outside of the model species Schistosoma mansoni. Here, we characterize the gene content and evolution of the W-chromosomes of S. mansoni and of the divergent species S. japonicum. We use a combined RNA/DNA k-mer based pipeline to assemble around one hundred candidate W-specific transcripts in each of the species. About half of them map to known protein coding genes, the majority homologous to S. mansoni Z-linked genes. We perform an extended analysis of the evolutionary strata present in the two species (including characterizing a previously undetected young stratum in S. japonicum) to infer patterns of sequence and expression evolution of W-linked genes at different time points after recombination was lost. W-linked genes show evidence of degeneration, including high rates of protein evolution and reduced expression. Most are found in young lineage-specific strata, with only a few high expression ancestral W-genes remaining, consistent with the progressive erosion of non-recombining regions. Among these, the splicing factor U2AF2 stands out as a promising candidate for primary sex determination, opening new avenues for understanding the molecular basis of the reproductive biology of this group.

Rational approach to drug discovery for human schistosomiasis

Human schistosomiasis is a debilitating, life-threatening disease affecting more than 229 million people in as many as 78 countries. There is only one drug of choice effective against all three major species of Schistosoma, praziquantel (PZQ). However, as with many monotherapies, evidence for resistance is emerging in the field and can be selected for in the laboratory. Previously used therapies include oxamniquine (OXA), but shortcomings such as drug resistance and affordability resulted in discontinuation. Employing a genetic, biochemical and molecular approach, a sulfotransferase (SULT-OR) was identified as responsible for OXA drug resistance. By crystallizing SmSULT- OR with OXA, the mode of action of OXA was determined. This information allowed a rational approach to novel drug design. Our team approach with schistosome biologists, medicinal chemists, structural biologists and geneticists has enabled us to develop and test novel drug derivatives of OXA to treat this disease. Using an iterative process for drug development, we have successfully identified derivatives that are effective against all three species of the parasite. One derivative CIDD-0149830 kills 100% of all three human schistosome species within 5 days. The goal is to generate a second therapeutic with a different mode of action that can be used in conjunction with praziquantel to overcome the ever-growing threat of resistance and improve efficacy. The ability and need to design, screen, and develop future, affordable therapeutics to treat human schistosomiasis is critical for successful control program outcomes.

Ancient Hybridization and Adaptive Introgression of an Invadolysin Gene in Schistosome Parasites

Introgression among parasite species has the potential to transfer traits of biomedical importance across species boundaries. The parasitic blood fluke Schistosoma haematobium causes urogenital schistosomiasis in humans across sub-Saharan Africa. Hybridization with other schistosome species is assumed to occur commonly, because genetic crosses between S. haematobium and livestock schistosomes, including S. bovis, can be staged in the laboratory, and sequencing of mtDNA and rDNA amplified from microscopic miracidia larvae frequently reveals markers from different species. However, the frequency, direction, age, and genomic consequences of hybridization are unknown. We hatched miracidia from eggs and sequenced the exomes from 96 individual S. haematobium miracidia from infected patients from Niger and the Zanzibar archipelago. These data revealed no evidence for contemporary hybridization between S. bovis and S. haematobium in our samples. However, all Nigerien S. haematobium genomes sampled show hybrid ancestry, with 3.3-8.2% of their nuclear genomes derived from S. bovis, providing evidence of an ancient introgression event that occurred at least 108-613 generations ago. Some S. bovis-derived alleles have spread to high frequency or reached fixation and show strong signatures of directional selection; the strongest signal spans a single gene in the invadolysin gene family (Chr. 4). Our results suggest that S. bovis/S. haematobium hybridization occurs rarely but demonstrate profound consequences of ancient introgression from a livestock parasite into the genome of S. haematobium, the most prevalent schistosome species infecting humans.

Population Structure and Dynamics of Helminthic Infection: Schistosomiasis

While disease and outbreaks are mainly clonal for bacteria and other asexually reproducing organisms, sexual reproduction in schistosomes and other helminths usually results in unique individuals. For sexually reproducing organisms, the traits conserved in clones will instead be conserved in the group of organisms that tends to breed together, the population. While the same tools are applied to characterize DNA, how results are interpreted can be quite different at times (see another article in this collection, http://www.asmscience.org/content/journal/microbiolspec/10.1128/microbiolspec.AME-0002-2018). It is difficult to know what the real effect any control program has on the parasite population without assessing the health of this population, how they respond to the control measure, and how they recover, if they do. This review, part of the Microbiology Spectrum Curated Collection: Advances in Molecular Epidemiology of Infectious Diseases, concentrates on one approach using pooled samples to study schistosome populations and shows how this and other approaches have contributed to our understanding of this parasite family's biology and epidemiology. *This article is part of a curated collection.

Schistosomiasis

Schistosomiasis is a major cause of morbidity in the world; it is second only to malaria as a major infectious disease. Globally, it is estimated that the disease affects over 250 million people in 78 countries of the world and is responsible for some 280,000 deaths each year. The three major schistosomes infecting humans are Schistosoma mansoni, S. japonicum, and S. haematobium. This chapter covers a wide range of aspects of schistosomiasis, including basic biology of the parasites, epidemiology, immunopathology, treatment, control, vaccines, and genomics/proteomics. In this chapter, the reader will understand the significant toll this disease takes in terms of mortality and morbidity. A description of the various life stages of schistosomes is presented, which will be informative for both those unfamiliar with the disease and experienced scientists. Clinical and public health aspects are addressed that cover acute and chronic disease, diagnosis, current treatment regimens and alternative drugs, and schistosomiasis control programs. A brief overview of genomics and proteomics is included that details recent advances in the field that will help scientists investigate the molecular biology of schistosomes. The reader will take away an appreciation for general aspects of schistosomiasis and research advances.

Parent-of-Origin-Dependent Gene Expression in Male and Female Schistosome Parasites

Schistosomes are the causative agents of schistosomiasis, a neglected tropical disease affecting over 230 million people worldwide. Additionally to their major impact on human health, they are also models of choice in evolutionary biology. These parasitic flatworms are unique among the common hermaphroditic trematodes as they have separate sexes. This so-called “evolutionary scandal” displays a female heterogametic genetic sex-determination system (ZZ males and ZW females), as well as a pronounced adult sexual dimorphism. These phenotypic differences are determined by a shared set of genes in both sexes, potentially leading to intralocus sexual conflicts. To resolve these conflicts in sexually selected traits, molecular mechanisms such as sex-biased gene expression could occur, but parent-of-origin gene expression also provides an alternative. In this work we investigated the latter mechanism, that is, genes expressed preferentially from either the maternal or the paternal allele, in Schistosoma mansoni species. To this end, transcriptomes from male and female hybrid adults obtained by strain crosses were sequenced. Strain-specific single nucleotide polymorphism (SNP) markers allowed us to discriminate the parental origin, while reciprocal crosses helped to differentiate parental expression from strain-specific expression. We identified genes containing SNPs expressed in a parent-of-origin manner consistent with paternal and maternal imprints. Although the majority of the SNPs was identified in mitochondrial and Z-specific loci, the remaining SNPs found in male and female transcriptomes were situated in genes that have the potential to explain sexual differences in schistosome parasites. Furthermore, we identified and validated four new Z-specific scaffolds.

Genetic Crosses and Linkage Mapping in Schistosome Parasites

Linkage mapping - utilizing experimental genetic crosses to examine cosegregation of phenotypic traits with genetic markers - is now 100 years old. Schistosome parasites are exquisitely well suited to linkage mapping approaches because genetic crosses can be conducted in the laboratory, thousands of progeny are produced, and elegant experimental work over the last 75 years has revealed heritable genetic variation in multiple biomedically important traits such as drug resistance, host specificity, and virulence. Application of this approach is timely because the improved genome assembly for Schistosoma mansoni and developing molecular toolkit for schistosomes increase our ability to link phenotype with genotype. We describe current progress and potential future directions of linkage mapping in schistosomes.

Evolution of gene dosage on the Z-chromosome of schistosome parasites

XY systems usually show chromosome-wide compensation of X-linked genes, while in many ZW systems, compensation is restricted to a minority of dosage-sensitive genes. Why such differences arose is still unclear. Here, we combine comparative genomics, transcriptomics and proteomics to obtain a complete overview of the evolution of gene dosage on the Z-chromosome of Schistosoma parasites. We compare the Z-chromosome gene content of African (Schistosoma mansoni and S. haematobium) and Asian (S. japonicum) schistosomes and describe lineage-specific evolutionary strata. We use these to assess gene expression evolution following sex-linkage. The resulting patterns suggest a reduction in expression of Z-linked genes in females, combined with upregulation of the Z in both sexes, in line with the first step of Ohno's classic model of dosage compensation evolution. Quantitative proteomics suggest that post-transcriptional mechanisms do not play a major role in balancing the expression of Z-linked genes.

Clonal amplification of Fasciola hepatica in Galba truncatula: within and between isolate variation of triclabendazole-susceptible and -resistant clones

Background

Fasciola hepatica is of worldwide significance, impacting on the health, welfare and productivity of livestock and regarded by WHO as a re-emerging zoonosis. Triclabendazole (TCBZ), the drug of choice for controlling acute fasciolosis in livestock, is also the drug used to treat human infections. However TCBZ-resistance is now considered a major threat to the effective control of F. hepatica. It has yet to be demonstrated whether F. hepatica undergoes a genetic clonal expansion in the snail intermediate host, Galba truncatula, and to what extent amplification of genotypes within the snail facilitates accumulation of drug resistant parasites. Little is known about genotypic and phenotypic variation within and between F. hepatica isolates.

Results

Six clonal isolates of F. hepatica (3× triclabendazole-resistant, TCBZ-R and 3× triclabendazole-susceptible, TCBZ-S) were generated. Snails infected with one miracidium started to shed cercariae 42–56 days post-infection and shed repeatedly up to a maximum of 11 times. A maximum of 884 cercariae were shed by one clonally-infected snail (FhLivS1) at a single time point, with > 3000 clonal metacercariae shed over its lifetime. Following experimental infection all 12 sheep were FEC positive at the time of TCBZ treatment. Sheep infected with one of three putative TCBZ-S clones and treated with TCBZ had no parasites in the liver at post-mortem, whilst sheep each infected with putative TCBZ-R isolates had 35–165 adult fluke at post-mortem, despite TCBZ treatment. All six untreated control animals had between 15–127 parasites. A single multi-locus genotype was reported for every fluke from each of the six clonal isolates. Adult F. hepatica showed considerable variation in weight, ranging from 20–280 mg, with variation in weight evident within and amongst clonal isolates.

Conclusions

A genetic clonal expansion occurs within G. truncatula, highlighting the potential for amplification of drug resistant genotypes of F. hepatica. Variation in the weight of parasites within and between clonal isolates and when comparing isolates that are either susceptible or resistant to TCBZ represent inherent variation in liver fluke and cannot be attributed to their resistance or susceptibility traits.

Size and Content of the Sex-Determining Region of the Y Chromosome in Dioecious Mercurialis annua, a Plant with Homomorphic Sex Chromosomes

Dioecious plants vary in whether their sex chromosomes are heteromorphic or homomorphic, but even homomorphic sex chromosomes may show divergence between homologues in the non-recombining, sex-determining region (SDR). Very little is known about the SDR of these species, which might represent particularly early stages of sex-chromosome evolution. Here, we assess the size and content of the SDR of the diploid dioecious herb Mercurialis annua, a species with homomorphic sex chromosomes and mild Y-chromosome degeneration. We used RNA sequencing (RNAseq) to identify new Y-linked markers for M. annua. Twelve of 24 transcripts showing male-specific expression in a previous experiment could be amplified by polymerase chain reaction (PCR) only from males, and are thus likely to be Y-linked. Analysis of genome-capture data from multiple populations of M. annua pointed to an additional six male-limited (and thus Y-linked) sequences. We used these markers to identify and sequence 17 sex-linked bacterial artificial chromosomes (BACs), which form 11 groups of non-overlapping sequences, covering a total sequence length of about 1.5 Mb. Content analysis of this region suggests that it is enriched for repeats, has low gene density, and contains few candidate sex-determining genes. The BACs map to a subset of the sex-linked region of the genetic map, which we estimate to be at least 14.5 Mb. This is substantially larger than estimates for other dioecious plants with homomorphic sex chromosomes, both in absolute terms and relative to their genome sizes. Our data provide a rare, high-resolution view of the homomorphic Y chromosome of a dioecious plant.

A Genome Resequencing-Based Genetic Map Reveals the Recombination Landscape of an Outbred Parasitic Nematode in the Presence of Polyploidy and Polyandry

The parasitic nematode Haemonchus contortus is an economically and clinically important pathogen of small ruminants, and a model system for understanding the mechanisms and evolution of traits such as anthelmintic resistance. Anthelmintic resistance is widespread and is a major threat to the sustainability of livestock agriculture globally; however, little is known about the genome architecture and parameters such as recombination that will ultimately influence the rate at which resistance may evolve and spread. Here, we performed a genetic cross between two divergent strains of H. contortus, and subsequently used whole-genome resequencing of a female worm and her brood to identify the distribution of genome-wide variation that characterizes these strains. Using a novel bioinformatic approach to identify variants that segregate as expected in a pseudotestcross, we characterized linkage groups and estimated genetic distances between markers to generate a chromosome-scale F1 genetic map. We exploited this map to reveal the recombination landscape, the first for any helminth species, demonstrating extensive variation in recombination rate within and between chromosomes. Analyses of these data also revealed the extent of polyandry, whereby at least eight males were found to have contributed to the genetic variation of the progeny analyzed. Triploid offspring were also identified, which we hypothesize are the result of nondisjunction during female meiosis or polyspermy. These results expand our knowledge of the genetics of parasitic helminths and the unusual life-history of H. contortus, and enhance ongoing efforts to understand the genetic basis of resistance to the drugs used to control these worms and for related species that infect livestock and humans throughout the world. This study also demonstrates the feasibility of using whole-genome resequencing data to directly construct a genetic map in a single generation cross from a noninbred nonmodel organism with a complex lifecycle.

Clonorchis sinensis and Clonorchiasis: The Relevance of Exploring Genetic Variation

Parasitic trematodes (flukes) cause substantial mortality and morbidity in humans. The Chinese liver fluke, Clonorchis sinensis, is one of the most destructive parasitic worms in humans in China, Vietnam, Korea and the Russian Far East. Although C. sinensis infection can be controlled relatively well using anthelmintics, the worm is carcinogenic, inducing cholangiocarcinoma and causing major suffering in ~15 million people in Asia. This chapter provides an account of C. sinensis and clonorchiasis research-covering aspects of biology, epidemiology, pathogenesis and immunity, diagnosis, treatment and control, genetics and genomics. It also describes progress in the area of molecular biology (genetics, genomics, transcriptomics and proteomics) and highlights challenges associated with comparative genomics and population genetics. It then reviews recent advances in the sequencing and characterisation of the mitochondrial and nuclear genomes for a Korean isolate of C. sinensis and summarises salient comparative genomic work and the implications thereof. The chapter concludes by considering how advances in genomic and informatics will enable research on the genetics of C. sinensis and related parasites, as well as the discovery of new fluke-specific intervention targets.

Genetic variability of Taenia solium cysticerci recovered from experimentally infected pigs and from naturally infected pigs using microsatellite markers

The adult Taenia solium, the pork tapeworm, usually lives as a single worm in the small intestine of humans, its only known definitive host. Mechanisms of genetic variation in T. solium are poorly understood. Using three microsatellite markers previously reported [1], this study explored the genetic variability of T. solium from cysts recovered from experimentally infected pigs. It then explored the genetic epidemiology and transmission in naturally infected pigs and adult tapeworms recovered from human carriers from an endemic rural community in Peru. In an initial study on experimental infection, two groups of three piglets were each infected with proglottids from one of two genetically different tapeworms for each of the microsatellites. After 7 weeks, pigs were slaughtered and necropsy performed. Thirty-six (92.3%) out of 39 cysts originated from one tapeworm, and 27 (100%) out of 27 cysts from the other had exactly the same genotype as the parental tapeworm. This suggests that the microsatellite markers may be a useful tool for studying the transmission of T. solium. In the second study, we analyzed the genetic variation of T. solium in cysts recovered from eight naturally infected pigs, and from adult tapeworms recovered from four human carriers; they showed genetic variability. Four pigs had cysts with only one genotype, and four pigs had cysts with two different genotypes, suggesting that multiple infections of genetically distinct parental tapeworms are possible. Six pigs harbored cysts with a genotype corresponding to one of the identified tapeworms from the human carriers. In the dendrogram, cysts appeared to cluster within the corresponding pigs as well as with the geographical origin, but this association was not statistically significant. We conclude that genotyping of microsatellite size polymorphisms is a potentially important tool to trace the spread of infection and pinpoint sources of infection as pigs spread cysts with a shared parental genotype.

Taenia solium, the pork tapeworm, is a major cause of epilepsy in developing countries. Although it has been deemed a potentially eradicable pathogen, it remains prevalent in rural communities. This two-part study aims to evaluate the utility of three microsatellite markers previously reported, to identify parasites and to establish relationships among them. In the first study, we evaluated the genetic variability of the progeny of two individual tapeworms by infecting two groups of three pigs each. We found variation of 8% and 0% in the two groups with respect to the parental tapeworm, indicating that the cysts source may be identifiable. Next, in the second study we described the genetic relationships among tapeworms obtained from four carriers and cysts obtained from eight naturally infected pigs in a rural community. We demonstrated that pigs can have two types of cysts, suggesting multiple infections. In addition, we found relatedness between 6 pigs and one tapeworm identified in the community. Our results indicate the potential for microsatellite markers to identify genetic relationships between parasites and thereby establish routes of transmission. It is likely that the limited number of microsatellites prevented us from establishing relatedness with more precision. Therefore, further evaluation of additional microsatellites is recommended.

A method for single pair mating in an obligate parasitic nematode

Parasitic nematode species have extremely high levels of genetic diversity, presenting a number of experimental challenges for genomic and genetic work. Consequently, there is a need to develop inbred laboratory strains with reduced levels of polymorphism. The most efficient approach to inbred line development is single pair mating, but this is challenging for obligate parasites where the adult sexual reproductive stages are inside the host, and thus difficult to experimentally manipulate. This paper describes a successful approach to single pair mating of a parasitic nematode, Haemonchus contortus. The method allows for polyandrous mating behaviour and involves the surgical transplantation of a single adult male worm with multiple immature adult females directly into the sheep abomasum. We used a panel of microsatellite markers to monitor and validate the single pair mating crosses and to ensure that the genotypes of progeny and subsequent filial generations were consistent with those expected from a mating between a single female parent of known genotype and a single male parent of unknown genotype. We have established two inbred lines that both show a significant overall reduction in genetic diversity based on microsatellite genotyping and genome-wide single nucleotide polymorphism. There was an approximately 50% reduction in heterozygous SNP sites across the genome in the MHco3.N1 line compared with the MoHco3(ISE) parental strain. The MHco3.N1 inbred line has subsequently been used to provide DNA template for whole genome sequencing of H. contortus. This work provides proof of concept and methodologies for forward genetic analysis of obligate parasitic nematodes.

A Genome Wide Comparison to Identify Markers to Differentiate the Sex of Larval Stages of Schistosoma haematobium, Schistosoma bovis and their Respective Hybrids

For scientists working on gonochoric organisms, determining sex can be crucial for many biological questions and experimental studies, such as crossbreeding, but it can also be a challenging task, particularly when no sexual dimorphism is visible or cannot be directly observed. In metazoan parasites of the genus Schistosoma responsible for schistosomiasis, sex is genetically determined in the zygote with a female heterogametic ZW/ZZ system. Adult flukes have a pronounced sexual dimorphism, whereas the sexes of the larval stages are morphologically indistinguishable but can be distinguished uniquely by using molecular methods. Therefore, reliable methods are needed to identify the sex of larvae individuals. Here, we present an endpoint PCR-based assay using female-specific sequences identified using a genome-wide comparative analysis between males and females. This work allowed us to identify sex-markers for Schistosoma haematobium and Schistosoma bovis but also the hybrid between both species that has recently emerged in Corsica (France). Five molecular sex-markers were identified and are female-specific in S. haematobium and the hybrid parasite, whereas three of them are also female-specific in S. bovis. These molecular markers will be useful to conduct studies, such as experimental crosses on these disease-causing blood flukes, which are still largely neglected but no longer restricted to tropical areas.

Current global changes (environmental and anthropogenic) are expected to promote the spread and transmission of infectious diseases. One of the direct consequences of such changes is the modification of the geographical distribution of species, enabling natural hybridization. Such hybridization is already known to occur in schistosomes, and offspring have been shown to have superior virulence and invasive capacities. The recent outbreak of a hybrid between the human- and animal- infecting schistosomes, S. haematobium x S. bovis, in Europe (Corsica, France) clearly demonstrates this invasive capacity and raises the risk of zoonotic transmission. Therefore, it is important to study such hybrids, and experimental crosses are critical to address this issue. Here, we developed molecular sex markers for S. haematobium and S. bovis in order to distinguish gender and to be able to generate differentially introgressed hybrids, allowing us to investigate parasite fitness.

The NIH-NIAID Schistosomiasis Resource Center at the Biomedical Research Institute: Molecular Redux

Schistosomiasis remains a health burden in many parts of the world. The complex life cycle of Schistosoma parasites and the economic and societal conditions present in endemic areas make the prospect of eradication unlikely in the foreseeable future. Continued and vigorous research efforts must therefore be directed at this disease, particularly since only a single World Health Organization (WHO)-approved drug is available for treatment. The National Institutes of Health (NIH)–National Institute of Allergy and Infectious Diseases (NIAID) Schistosomiasis Resource Center (SRC) at the Biomedical Research Institute provides investigators with the critical raw materials needed to carry out this important research. The SRC makes available, free of charge (including international shipping costs), not only infected host organisms but also a wide array of molecular reagents derived from all life stages of each of the three main human schistosome parasites. As the field of schistosomiasis research rapidly advances, it is likely to become increasingly reliant on omics, transgenics, epigenetics, and microbiome-related research approaches. The SRC has and will continue to monitor and contribute to advances in the field in order to support these research efforts with an expanding array of molecular reagents. In addition to providing investigators with source materials, the SRC has expanded its educational mission by offering a molecular techniques training course and has recently organized an international schistosomiasis-focused meeting. This review provides an overview of the materials and services that are available at the SRC for schistosomiasis researchers, with a focus on updates that have occurred since the original overview in 2008.

Sex-Biased Transcriptome of Schistosoma mansoni: Host-Parasite Interaction, Genetic Determinants and Epigenetic Regulators Are Associated with Sexual Differentiation

Background

Among more than 20,000 species of hermaphroditic trematodes, Schistosomatidae are unusual since they have evolved gonochorism. In schistosomes, sex is determined by a female heterogametic system, but phenotypic sexual dimorphism appears only after infection of the vertebrate definitive host. The completion of gonad maturation occurs even later, after pairing. To date, the molecular mechanisms that trigger the sexual differentiation in these species remain unknown, and in vivo studies on the developing schistosomulum stages are lacking. To study the molecular basis of sex determination and sexual differentiation in schistosomes, we investigated the whole transcriptome of the human parasite Schistosoma mansoni in a stage- and sex-comparative manner.

Methodology/ Principal Findings

We performed a RNA-seq on males and females for five developmental stages: cercariae larvae, three in vivo schistosomulum stages and adults. We detected 7,168 genes differentially expressed between sexes in at least one of the developmental stages, and 4,065 of them were functionally annotated. Transcriptome data were completed with H3K27me3 histone modification analysis using ChIP-Seq before (in cercariae) and after (in adults) the phenotypic sexual dimorphism appearance. In this paper we present (i) candidate determinants of the sexual differentiation, (ii) sex-biased players of the interaction with the vertebrate host, and (iii) different dynamic of the H3K27me3 histone mark between sexes as an illustration of sex-biased epigenetic landscapes.

Conclusions/ Significance

Our work presents evidence that sexual differentiation in S. mansoni is accompanied by distinct male and female transcriptional landscapes of known players of the host-parasite crosstalk, genetic determinants and epigenetic regulators. Our results suggest that such combination could lead to the optimized sexual dimorphism of this parasitic species. As S. mansoni is pathogenic for humans, this study represents a promising source of therapeutic targets, providing not only data on the parasite development in interaction with its vertebrate host, but also new insights on its reproductive function.

Parasitic flatworms include more than 20,000 species that are classically hermaphrodites. Among them, the roughly hundred species of Schistosomatidae are intriguing because they are gonochoric. Schistosomes are responsible of the second most important parasitic disease worldwide, and eggs are the main cause of the inflammatory symptoms. Thus, studying the sexual reproduction mechanisms of schistosomes is of particular interest for drug development. Schistosome’s sex is genetically determined by the presence of sex chromosomes: ZZ in males or ZW in females. There is, however, no phenotypic dimorphism in the larval stages: sexual dimorphism appears only in the vertebrate host. In order to understand the molecular mechanisms underlying phenotypic sexual dimorphism, we performed a transcriptome analysis (RNA-Seq) in five different stages of the parasite lifecycle as well as a chromatin status analysis (ChIP-Seq) in the non-differentiated stage cercariae and in the adult differentiated stage, for males and females separately. Our work presents evidence that sexual differentiation in S. mansoni is accompanied by distinct male and female transcriptional landscapes of known players of the host-parasite crosstalk, developmental pathways and epigenetic regulators. Our sex-comparative approach provides therefore new potential therapeutic targets to affect development and sexual reproduction of parasite.

Real-time PCR for sexing Schistosoma mansoni cercariae

The gender of cercarial larvae can only be determined using molecular methods. End point PCR methods that amplify repetitive markers on the W chromosome of the female (ZW) parasites have been developed, but sometimes results are ambiguous or incorrect. To more effectively distinguish sexes, and to determine why end point PCR can be incorrect, we quantified the W6 repeat sequence and a specific Z chromosome gene using real-time PCR. The ratio between copy number of W6 and a Z chromosome marker unambiguously identifies gender: females have higher ratios (421-4371) than males (0-21). However, some males have low numbers of W6 elements in their genome, and qPCR demonstrated significantly higher W6/Z marker ratios for male genotypes giving ambiguous end point PCR results compared with males giving clear end point results. The quantitative PCR sexing method developed will be particularly useful where reliable sexing of cercariae is critical, for example when staging genetic crosses.

Uncovering multiloci-ordering by algebraic property of Laplacian matrix and its Fiedler vector

MOTIVATION

The loci-ordering, based on two-point recombination fractions for a pair of loci, is the most important step in constructing a reliable and fine genetic map.

RESULTS

Using the concept from complex graph theory, here we propose a Laplacian ordering approach which uncovers the loci-ordering of multiloci simultaneously. The algebraic property for a Fiedler vector of a Laplacian matrix, constructed from the recombination fraction of the loci-ordering for 26 loci of barley chromosome IV, 846 loci of Arabidopsis thaliana and 1903 loci of Malus domestica, together with the variable threshold uncovers their loci-orders. It offers an alternative yet robust approach for ordering multiloci.

AVAILABILITY AND IMPLEMENTATION

Source code program with data set is available as supplementary data and also in a software category of the website (http://biophysics.dgist.ac.kr)

CONTACT

crkim@pusan.ac.kr or iksoochang@dgist.ac.kr

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

The Pseudoautosomal Regions of the U/V Sex Chromosomes of the Brown Alga Ectocarpus Exhibit Unusual Features

The recombining regions of sex chromosomes (pseudoautosomal regions, PARs) are predicted to exhibit unusual features due to their being genetically linked to the nonrecombining, sex-determining region. This phenomenon is expected to occur in both diploid (XY, ZW) and haploid (UV) sexual systems, with slightly different consequences for UV sexual systems because of the absence of masking during the haploid phase (when sex is expressed) and because there is no homozygous sex in these systems. Despite a considerable amount of theoretical work on PAR genetics and evolution, these genomic regions have remained poorly characterized empirically. We show here that although the PARs of the U/V sex chromosomes of the brown alga Ectocarpus recombine at a similar rate to autosomal regions of the genome, they exhibit many genomic features typical of nonrecombining regions. The PARs were enriched in clusters of genes that are preferentially, and often exclusively, expressed during the sporophyte generation of the life cycle, and many of these genes appear to have evolved since the Ectocarpales diverged from other brown algal lineages. A modeling-based approach was used to investigate possible evolutionary mechanisms underlying this enrichment in sporophyte-biased genes. Our results are consistent with the evolution of the PAR in haploid systems being influenced by differential selection pressures in males and females acting on alleles that are advantageous during the sporophyte generation of the life cycle.

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.

Chromosomal differentiation of schistosomes: what is the message?

As the only group of flukes with dioecism, schistosomes are unique organisms; they not only have intriguing biological and evolutionary aspects but also are responsible for major public health problems in the developing world. Schistosomiasis caused by this fluke affects approximately 210 million people in 76 countries. In order to facilitate the discovery of eradication methods for this disease, fundamental biological outcomes must be made available. Whole genome sequence data represent one such resource applicable to discovering eradication methods and measures. Herein, I describe three remarkable chromosomal changes and briefly discuss the differentiation of the Asian and African groups of this parasite taxon. Chromosomal data and evolutionary aspects will enable us to exploit genomic information for advancing schistosome studies.

Deciphering the glycogenome of schistosomes

Schistosoma mansoni and other Schistosoma sp. are multicellular parasitic helminths (worms) that infect humans and mammals worldwide. Infection by these parasites, which results in developmental maturation and sexual differentiation of the worms over a period of 5–6 weeks, induces antibodies to glycan antigens expressed in surface and secreted glycoproteins and glycolipids. There is growing interest in defining these unusual parasite-synthesized glycan antigens and using them to understand immune responses, their roles in immunomodulation, and in using glycan antigens as potential vaccine targets. A key problem in this area, however, has been the lack of information about the enzymes involved in elaborating the complex repertoire of glycans represented by the schistosome glycome. Recent availability of the nuclear genome sequences for Schistosoma sp. has created the opportunity to define the glycogenome, which represents the specific genes and cognate enzymes that generate the glycome. Here we describe the current state of information in regard to the schistosome glycogenome and glycome and highlight the important classes of glycans and glycogenes that may be important in their generation.

Gaining biological perspectives from schistosome genomes

Characterization of the genomic basis underlying schistosome biology is an important strategy for the development of future treatments and interventions. Genomic sequence is now available for the three major clinically relevant schistosome species, Schistosoma mansoni, S. japonicum and S. haematobium, and this information represents an invaluable resource for the future control of human schistosomiasis. The identification of a biologically important, but distinct from the host, schistosome gene product is the ultimate goal for many research groups. While the initial elucidation of the genome of an organism is critical for most biological research, continued improvement or curation of the genome construction should be an ongoing priority. In this review we will discuss prominent recent findings utilizing a systems approach to schistosome biology, as well as the increased use of interference RNA (RNAi). Both of these research strategies are aiming to place parasite genes into a more meaningful biological perspective.

Efficient linkage mapping using exome capture and extreme QTL in schistosome parasites

Background

Identification of parasite genes that underlie traits such as drug resistance and host specificity is challenging using classical linkage mapping approaches. Extreme QTL (X-QTL) methods, originally developed by rodent malaria and yeast researchers, promise to increase the power and simplify logistics of linkage mapping in experimental crosses of schistosomes (or other helminth parasites), because many 1000s of progeny can be analysed, phenotyping is not required, and progeny pools rather than individuals are genotyped. We explored the utility of this method for mapping a drug resistance gene in the human parasitic fluke Schistosoma mansoni.

Results

We staged a genetic cross between oxamniquine sensitive and resistant parasites, then between two F1 progeny, to generate multiple F2 progeny. One group of F2s infecting hamsters was treated with oxamniquine, while a second group was left untreated. We used exome capture to reduce the size of the genome (from 363 Mb to 15 Mb) and exomes from pooled F2 progeny (treated males, untreated males, treated females, untreated females) and the two parent parasites were sequenced to high read depth (mean = 95-366×) and allele frequencies at 14,489 variants compared. We observed dramatic enrichment of alleles from the resistant parent in a small region of chromosome 6 in drug-treated male and female pools (combined analysis: = 11.07, p = 8.74 × 10^-29). This region contains Smp_089320 a gene encoding a sulfotransferase recently implicated in oxamniquine resistance using classical linkage mapping methods.

Conclusions

These results (a) demonstrate the utility of exome capture for generating reduced representation libraries in Schistosoma mansoni, and (b) provide proof-of-principle that X-QTL methods can be successfully applied to an important human helminth. The combination of these methods will simplify linkage analysis of biomedically or biologically important traits in this parasite.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-617) contains supplementary material, which is available to authorized users.

Recombination in the human Pseudoautosomal region PAR1

The pseudoautosomal region (PAR) is a short region of homology between the mammalian X and Y chromosomes, which has undergone rapid evolution. A crossover in the PAR is essential for the proper disjunction of X and Y chromosomes in male meiosis, and PAR deletion results in male sterility. This leads the human PAR with the obligatory crossover, PAR1, to having an exceptionally high male crossover rate, which is 17-fold higher than the genome-wide average. However, the mechanism by which this obligatory crossover occurs remains unknown, as does the fine-scale positioning of crossovers across this region. Recent research in mice has suggested that crossovers in PAR may be mediated independently of the protein PRDM9, which localises virtually all crossovers in the autosomes. To investigate recombination in this region, we construct the most fine-scale genetic map containing directly observed crossovers to date using African-American pedigrees. We leverage recombination rates inferred from the breakdown of linkage disequilibrium in human populations and investigate the signatures of DNA evolution due to recombination. Further, we identify direct PRDM9 binding sites using ChIP-seq in human cells. Using these independent lines of evidence, we show that, in contrast with mouse, PRDM9 does localise peaks of recombination in the human PAR1. We find that recombination is a far more rapid and intense driver of sequence evolution in PAR1 than it is on the autosomes. We also show that PAR1 hotspot activities differ significantly among human populations. Finally, we find evidence that PAR1 hotspot positions have changed between human and chimpanzee, with no evidence of sharing among the hottest hotspots. We anticipate that the genetic maps built and validated in this work will aid research on this vital and fascinating region of the genome.

Praziquantel treatment decreases Schistosoma mansoni genetic diversity in experimental infections

BACKGROUND

Schistosomiasis has a considerable impact on public health in many tropical and subtropical areas. In the new world, schistosomiasis is caused by the digenetic trematode Schistosoma mansoni. Chemotherapy is the main measure for controlling schistosomiasis, and the current drug of choice for treatment is praziquantel (PZQ). Although PZQ is efficient and safe, its repetitive large-scale use in endemic areas may lead to the selection of resistant strains. Isolates less susceptible to PZQ have been found in the field and selected for in the laboratory. The impact of selecting strains with a decreased susceptibility phenotype on disease dynamics and parasite population genetics is not fully understood. This study addresses the impact of PZQ pressure on the genetics of a laboratory population by analyzing frequency variations of polymorphic genetic markers.

METHODOLOGY

Infected mice were treated with increasing PZQ doses until the highest dose of 3 × 300 mg/Kg was reached. The effect of PZQ treatment on the parasite population was assessed using five polymorphic microsatellite markers. Parasitological and genetic data were compared with those of the untreated control. After six parasite generations submitted to treatment, it was possible to obtain a S. mansoni population with decreased susceptibility to PZQ. In our experiments we also observed that female worms were more susceptible to PZQ than male worms.

CONCLUSIONS

The selective pressure exerted by PZQ led to decreased genetic variability in S. mansoni and increased endogamy. The understanding of how S. mansoni populations respond to successive drug pressure has important implications on the appearance and maintenance of a PZQ resistance phenotype in endemic regions.

Genetic and molecular basis of drug resistance and species-specific drug action in schistosome parasites

Oxamniquine resistance evolved in the human blood fluke (Schistosoma mansoni) in Brazil in the 1970s. We crossed parental parasites differing ~500-fold in drug response, determined drug sensitivity and marker segregation in clonally derived second-generation progeny, and identified a single quantitative trait locus (logarithm of odds = 31) on chromosome 6. A sulfotransferase was identified as the causative gene by using RNA interference knockdown and biochemical complementation assays, and we subsequently demonstrated independent origins of loss-of-function mutations in field-derived and laboratory-selected resistant parasites. These results demonstrate the utility of linkage mapping in a human helminth parasite, while crystallographic analyses of protein-drug interactions illuminate the mode of drug action and provide a framework for rational design of oxamniquine derivatives that kill both S. mansoni and S. haematobium, the two species responsible for >99% of schistosomiasis cases worldwide.

Identification of a sex-linked SNP marker in the salmon louse (Lepeophtheirus salmonis) using RAD sequencing

The salmon louse (Lepeophtheirus salmonis (Krøyer, 1837)) is a parasitic copepod that can, if untreated, cause considerable damage to Atlantic salmon (Salmo salar Linnaeus, 1758) and incurs significant costs to the Atlantic salmon mariculture industry. Salmon lice are gonochoristic and normally show sex ratios close to 1:1. While this observation suggests that sex determination in salmon lice is genetic, with only minor environmental influences, the mechanism of sex determination in the salmon louse is unknown. This paper describes the identification of a sex-linked Single Nucleotide Polymorphism (SNP) marker, providing the first evidence for a genetic mechanism of sex determination in the salmon louse. Restriction site-associated DNA sequencing (RAD-seq) was used to isolate SNP markers in a laboratory-maintained salmon louse strain. A total of 85 million raw Illumina 100 base paired-end reads produced 281,838 unique RAD-tags across 24 unrelated individuals. RAD marker Lsa101901 showed complete association with phenotypic sex for all individuals analysed, being heterozygous in females and homozygous in males. Using an allele-specific PCR assay for genotyping, this SNP association pattern was further confirmed for three unrelated salmon louse strains, displaying complete association with phenotypic sex in a total of 96 genotyped individuals. The marker Lsa101901 was located in the coding region of the prohibitin-2 gene, which showed a sex-dependent differential expression, with mRNA levels determined by RT-qPCR about 1.8-fold higher in adult female than adult male salmon lice. This study’s observations of a novel sex-linked SNP marker are consistent with sex determination in the salmon louse being genetic and following a female heterozygous system. Marker Lsa101901 provides a tool to determine the genetic sex of salmon lice, and could be useful in the development of control strategies.

Sex-biased expression of microRNAs in Schistosoma mansoni

Schistosomiasis is an important neglected tropical disease caused by digenean helminth parasites of the genus Schistosoma. Schistosomes are unusual in that they are dioecious and the adult worms live in the blood system. MicroRNAs play crucial roles during gene regulation and are likely to be important in sex differentiation in dioecious species. Here we characterize 112 microRNAs from adult Schistosoma mansoni individuals, including 84 novel microRNA families, and investigate the expression pattern in different sexes. By deep sequencing, we measured the relative expression levels of conserved and newly identified microRNAs between male and female samples. We observed that 13 microRNAs exhibited sex-biased expression, 10 of which are more abundant in females than in males. Sex chromosomes showed a paucity of female-biased genes, as predicted by theoretical evolutionary models. We propose that the recent emergence of separate sexes in Schistosoma had an effect on the chromosomal distribution and evolution of microRNAs, and that microRNAs are likely to participate in the sex differentiation/maintenance process.

Schistosomiasis is the second most common disease caused by a parasite, affecting over 200 million people. The parasites involved are flatworms of the genus Schistosoma. Unlike most non-parasitic flatworms, Schistosoma species have separate sexes, and the emergence of sex has been associated with the development of a parasitic lifestyle. The identification of gene products that are expressed in a sex-biased fashion permits the study of the origin of sexual dimorphism and, in the case of the schistosomes, the evolution of a human parasite. Here we investigated the differential expression of microRNAs in male and female individuals of the species Schistosoma mansoni. MicroRNAs are crucial gene regulators. We observed that many new microRNAs emerged in the evolutionary lineage leading to the schistosomes. However, many sex-biased microRNAs were present in the hermaphrodite ancestor of the flatworms, and therefore acquired sex-biased expression later on. Our results suggest that changes in microRNA expression patterns were associated with the emergence of separate sexes in the schistosomes.

New approaches for understanding mechanisms of drug resistance in schistosomes

Schistosomes are parasitic flatworms that cause schistosomiasis, a neglected tropical disease that affects hundreds of millions worldwide. Treatment and control of schistosomiasis relies almost entirely on the single drug praziquantel (PZQ), making the prospect of emerging drug resistance particularly worrisome. This review will survey reports of PZQ (and other drug) resistance in schistosomes and other platyhelminths, and explore mechanisms by which drug resistance might develop. Newer genomic and post-genomic strategies that offer the promise of better understanding of how drug resistance might arise in these organisms will be discussed. These approaches could also lead to insights into the mode of action of these drugs and potentially provide markers for monitoring the emergence of resistance.

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.

The Schistosoma mansoni phylome: using evolutionary genomics to gain insight into a parasite's biology

BACKGROUND

Schistosoma mansoni is one of the causative agents of schistosomiasis, a neglected tropical disease that affects about 237 million people worldwide. Despite recent efforts, we still lack a general understanding of the relevant host-parasite interactions, and the possible treatments are limited by the emergence of resistant strains and the absence of a vaccine. The S. mansoni genome was completely sequenced and still under continuous annotation. Nevertheless, more than 45% of the encoded proteins remain without experimental characterization or even functional prediction. To improve our knowledge regarding the biology of this parasite, we conducted a proteome-wide evolutionary analysis to provide a broad view of the S. mansoni's proteome evolution and to improve its functional annotation.

RESULTS

Using a phylogenomic approach, we reconstructed the S. mansoni phylome, which comprises the evolutionary histories of all parasite proteins and their homologs across 12 other organisms. The analysis of a total of 7,964 phylogenies allowed a deeper understanding of genomic complexity and evolutionary adaptations to a parasitic lifestyle. In particular, the identification of lineage-specific gene duplications pointed to the diversification of several protein families that are relevant for host-parasite interaction, including proteases, tetraspanins, fucosyltransferases, venom allergen-like proteins, and tegumental-allergen-like proteins. In addition to the evolutionary knowledge, the phylome data enabled us to automatically re-annotate 3,451 proteins through a phylogenetic-based approach rather than solely sequence similarity searches. To allow further exploitation of this valuable data, all information has been made available at PhylomeDB (http://www.phylomedb.org).

CONCLUSIONS

In this study, we used an evolutionary approach to assess S. mansoni parasite biology, improve genome/proteome functional annotation, and provide insights into host-parasite interactions. Taking advantage of a proteome-wide perspective rather than focusing on individual proteins, we identified that this parasite has experienced specific gene duplication events, particularly affecting genes that are potentially related to the parasitic lifestyle. These innovations may be related to the mechanisms that protect S. mansoni against host immune responses being important adaptations for the parasite survival in a potentially hostile environment. Continuing this work, a comparative analysis involving genomic, transcriptomic, and proteomic data from other helminth parasites, other parasites, and vectors will supply more information regarding parasite's biology as well as host-parasite interactions.

Schistosoma tegument proteins in vaccine and diagnosis development: an update

The development of a vaccine against schistosomiasis and also the availability of a more sensitive diagnosis test are important tools to help chemotherapy in controlling disease transmission. Bioinformatics tools, together with the access to parasite genome, published recently, should help generate new knowledge on parasite biology and search for new vaccines or therapeutic targets and antigens to be used in the disease diagnosis. Parasite surface proteins, especially those expressed in schistosomula tegument, represent interesting targets to be used in vaccine formulations and in the diagnosis of early infections, since the tegument represents the interface between host and parasite and its molecules are responsible for essential functions to parasite survival. In this paper we will present the advances in the development of vaccines and diagnosis tests achieved with the use of the information from schistosome genome focused on parasite tegument as a source for antigens.

Differences in genomic architecture between two distinct geographical strains of the blood fluke Schistosoma japonicum reveal potential phenotype basis

The Chinese (SjC) and Philippine (SjP) strains of the blood fluke Schistosoma japonicum have been shown to present clearly different phenotypes in fecundity, pathology, drug sensitivity and immunology. We used microarray based comparative genomic hybridisation (aCGH) to investigate structural differences in the genomes of the two strains and identified seven distinct regions of the S. japonicum genome that present differential aCGH representing either deletion or duplication regions in SjP. Within these regions, genes predicted to be associated with the recognised phenotypic differences were identified and that may provide new insights into the biology and evolution of the two strains, with implications for the epidemiology and control of schistosomiasis japonica in China and the Philippines.

Evolution of sex chromosomes ZW of Schistosoma mansoni inferred from chromosome paint and BAC mapping analyses

Chromosomes of schistosome parasites among digenetic flukes have a unique evolution because they exhibit the sex chromosomes ZW, which are not found in the other groups of flukes that are hermaphrodites. We conducted molecular cytogenetic analyses for investigating the sex chromosome evolution using chromosome paint analysis and BAC clones mapping. To carry this out, we developed a technique for making paint probes of genomic DNA from a single scraped chromosome segment using a chromosome microdissection system, and a FISH mapping technique for BAC clones. Paint probes clearly identified each of the 8 pairs of chromosomes by a different fluorochrome color. Combination analysis of chromosome paint analysis with Z/W probes and chromosome mapping with 93 BAC clones revealed that the W chromosome of Schistosoma mansoni has evolved by at least four inversion events and heterochromatinization. Nine of 93 BAC clones hybridized with both the Z and W chromosomes, but the locations were different between Z and W chromosomes. The homologous regions were estimated to have moved from the original Z chromosome to the differentiated W chromosome by three inversions events that occurred before W heterohcromatinization. An inversion that was observed in the heterochromatic region of the W chromosome likely occurred after W heterochromatinization. These inversions and heterochromatinization are hypothesized to be the key factors that promoted the evolution of the W chromosome of S. mansoni.

Recent advances in Schistosoma genomics

Schistosome research has entered the genomic era with the publications reporting the Schistosoma mansoni and Schistosoma japonicum genomes. Schistosome genomics is motivated by the need for new control tools. However, much can also be learned about the biology of Schistosoma, which is a tractable experimental model. In this article, we review the recent achievements in the field of schistosome research and discuss future perspectives on genomics and how it can be integrated in a usable format, on the genetic mapping and how it has improved the genome assembly and provided new research approaches, on how epigenetics provides interesting insights into the biology of the species and on new functional genomics tools that will contribute to the understanding of the function of genes, many of which are parasite- or taxon specific.

Chromatin structural changes around satellite repeats on the female sex chromosome in Schistosoma mansoni and their possible role in sex chromosome emergence

BACKGROUND

In the leuphotrochozoan parasitic platyhelminth Schistosoma mansoni, male individuals are homogametic (ZZ) whereas females are heterogametic (ZW). To elucidate the mechanisms that led to the emergence of sex chromosomes, we compared the genomic sequence and the chromatin structure of male and female individuals. As for many eukaryotes, the lower estimate for the repeat content is 40%, with an unknown proportion of domesticated repeats. We used massive sequencing to de novo assemble all repeats, and identify unambiguously Z-specific, W-specific and pseudoautosomal regions of the S. mansoni sex chromosomes.

RESULTS

We show that 70 to 90% of S. mansoni W and Z are pseudoautosomal. No female-specific gene could be identified. Instead, the W-specific region is composed almost entirely of 36 satellite repeat families, of which 33 were previously unknown. Transcription and chromatin status of female-specific repeats are stage-specific: for those repeats that are transcribed, transcription is restricted to the larval stages lacking sexual dimorphism. In contrast, in the sexually dimorphic adult stage of the life cycle, no transcription occurs. In addition, the euchromatic character of histone modifications around the W-specific repeats decreases during the life cycle. Recombination repression occurs in this region even if homologous sequences are present on both the Z and W chromosomes.

CONCLUSION

Our study provides for the first time evidence for the hypothesis that, at least in organisms with a ZW type of sex chromosomes, repeat-induced chromatin structure changes could indeed be the initial event in sex chromosome emergence.

About PAR: the distinct evolutionary dynamics of the pseudoautosomal region

Sex chromosomes differ from other chromosomes in the striking divergence they often show in size, structure, and gene content. Not only do they possess genes controlling sex determination that are restricted to either the X or Y (or Z or W) chromosomes, but in many taxa they also include recombining regions. In these 'pseudoautosomal regions' (PARs), sequence homology is maintained by meiotic pairing and exchange in the heterogametic sex. PARs are unique genomic regions, exhibiting some features of autosomes, but they are also influenced by their partial sex linkage. Here we review the distribution and structure of PARs among animals and plants, the theoretical predictions concerning their evolutionary dynamics, the reasons for their persistence, and the diversity and content of genes that reside within them. It is now clear that the evolution of the PAR differs in important ways from that of genes in either the non-recombining regions of sex chromosomes or the autosomes.

Whole-genome sequence of Schistosoma haematobium

Schistosomiasis is a neglected tropical disease caused by blood flukes (genus Schistosoma; schistosomes) and affecting 200 million people worldwide. No vaccines are available, and treatment relies on one drug, praziquantel. Schistosoma haematobium has come into the spotlight as a major cause of urogenital disease, as an agent linked to bladder cancer and as a predisposing factor for HIV/AIDS. The parasite is transmitted to humans from freshwater snails. Worms dwell in blood vessels and release eggs that become embedded in the bladder wall to elicit chronic immune-mediated disease and induce squamous cell carcinoma. Here we sequenced the 385-Mb genome of S. haematobium using Illumina-based technology at 74-fold coverage and compared it to sequences from related parasites. We included genome annotation based on function, gene ontology, networking and pathway mapping. This genome now provides an unprecedented resource for many fundamental research areas and shows great promise for the design of new disease interventions.

A systematically improved high quality genome and transcriptome of the human blood fluke Schistosoma mansoni

Schistosomiasis is one of the most prevalent parasitic diseases, affecting millions of people in developing countries. Amongst the human-infective species, Schistosoma mansoni is also the most commonly used in the laboratory and here we present the systematic improvement of its draft genome. We used Sanger capillary and deep-coverage Illumina sequencing from clonal worms to upgrade the highly fragmented draft 380 Mb genome to one with only 885 scaffolds and more than 81% of the bases organised into chromosomes. We have also used transcriptome sequencing (RNA-seq) from four time points in the parasite's life cycle to refine gene predictions and profile their expression. More than 45% of predicted genes have been extensively modified and the total number has been reduced from 11,807 to 10,852. Using the new version of the genome, we identified trans-splicing events occurring in at least 11% of genes and identified clear cases where it is used to resolve polycistronic transcripts. We have produced a high-resolution map of temporal changes in expression for 9,535 genes, covering an unprecedented dynamic range for this organism. All of these data have been consolidated into a searchable format within the GeneDB (www.genedb.org) and SchistoDB (www.schistodb.net) databases. With further transcriptional profiling and genome sequencing increasingly accessible, the upgraded genome will form a fundamental dataset to underpin further advances in schistosome research.