Daily rhythms in gene expression of the human parasite Schistosoma mansoni.

BACKGROUND: The consequences of the earth's daily rotation have led to 24-h biological rhythms in most organisms. Even some parasites are known to have daily rhythms, which, when in synchrony with host rhythms, can optimise their fitness. Understanding these rhythms may enable the development of control strategies that take advantage of rhythmic vulnerabilities. Recent work on protozoan parasites has revealed 24-h rhythms in gene expression, drug sensitivity and the presence of an intrinsic circadian clock; however, similar studies on metazoan parasites are lacking. To address this, we investigated if a metazoan parasite has daily molecular oscillations, whether they reveal how these longer-lived organisms can survive host daily cycles over a lifespan of many years and if animal circadian clock genes are present and rhythmic. We addressed these questions using the human blood fluke Schistosoma mansoni that lives in the vasculature for decades and causes the tropical disease schistosomiasis. RESULTS: Using round-the-clock transcriptomics of male and female adult worms collected from experimentally infected mice, we discovered that ~ 2% of its genes followed a daily pattern of expression. Rhythmic processes included a stress response during the host's active phase and a 'peak in metabolic activity' during the host's resting phase. Transcriptional profiles in the female reproductive system were mirrored by daily patterns in egg laying (eggs are the main drivers of the host pathology). Genes cycling with the highest amplitudes include predicted drug targets and a vaccine candidate. These 24-h rhythms may be driven by host rhythms and/or generated by a circadian clock; however, orthologs of core clock genes are missing and secondary clock genes show no 24-h rhythmicity. CONCLUSIONS: There are daily rhythms in the transcriptomes of adult S. mansoni, but they appear less pronounced than in other organisms. The rhythms reveal temporally compartmentalised internal processes and host interactions relevant to within-host survival and between-host transmission. Our findings suggest that if these daily rhythms are generated by an intrinsic circadian clock then the oscillatory mechanism must be distinct from that in other animals. We have shown which transcripts oscillate at this temporal scale and this will benefit the development and delivery of treatments against schistosomiasis.

Duplication and retention biases of essential and non-essential genes revealed by systematic knockdown analyses.

When a duplicate gene has no apparent loss-of-function phenotype, it is commonly considered that the phenotype has been masked as a result of functional redundancy with the remaining paralog. This is supported by indirect evidence showing that multi-copy genes show loss-of-function phenotypes less often than single-copy genes and by direct tests of phenotype masking using select gene sets. Here we take a systematic genome-wide RNA interference approach to assess phenotype masking in paralog pairs in the Caenorhabditis elegans genome. Remarkably, in contrast to expectations, we find that phenotype masking makes only a minor contribution to the low knockdown phenotype rate for duplicate genes. Instead, we find that non-essential genes are highly over-represented among duplicates, leading to a low observed loss-of-function phenotype rate. We further find that duplicate pairs derived from essential and non-essential genes have contrasting evolutionary dynamics: whereas non-essential genes are both more often successfully duplicated (fixed) and lost, essential genes are less often duplicated but upon successful duplication are maintained over longer periods. We expect the fundamental evolutionary duplication dynamics presented here to be broadly applicable.

The genome of the blood fluke Schistosoma mansoni.

Schistosoma mansoni is responsible for the neglected tropical disease schistosomiasis that affects 210 million people in 76 countries. Here we present analysis of the 363 megabase nuclear genome of the blood fluke. It encodes at least 11,809 genes, with an unusual intron size distribution, and new families of micro-exon genes that undergo frequent alternative splicing. As the first sequenced flatworm, and a representative of the Lophotrochozoa, it offers insights into early events in the evolution of the animals, including the development of a body pattern with bilateral symmetry, and the development of tissues into organs. Our analysis has been informed by the need to find new drug targets. The deficits in lipid metabolism that make schistosomes dependent on the host are revealed, and the identification of membrane receptors, ion channels and more than 300 proteases provide new insights into the biology of the life cycle and new targets. Bioinformatics approaches have identified metabolic chokepoints, and a chemogenomic screen has pinpointed schistosome proteins for which existing drugs may be active. The information generated provides an invaluable resource for the research community to develop much needed new control tools for the treatment and eradication of this important and neglected disease.

Using ChEMBL to Complement Schistosome Drug Discovery.

Schistosomiasis is one of the most important neglected tropical diseases. Until an effective vaccine is registered for use, the cornerstone of schistosomiasis control remains chemotherapy with praziquantel. The sustainability of this strategy is at substantial risk due to the possibility of praziquantel insensitive/resistant schistosomes developing. Considerable time and effort could be saved in the schistosome drug discovery pipeline if available functional genomics, bioinformatics, cheminformatics and phenotypic resources are systematically leveraged. Our approach, described here, outlines how schistosome-specific resources/methodologies, coupled to the open-access drug discovery database ChEMBL, can be cooperatively used to accelerate early-stage, schistosome drug discovery efforts. Our process identified seven compounds (fimepinostat, trichostatin A, NVP-BEP800, luminespib, epoxomicin, CGP60474 and staurosporine) with ex vivo anti-schistosomula potencies in the sub-micromolar range. Three of those compounds (epoxomicin, CGP60474 and staurosporine) also demonstrated potent and fast-acting ex vivo effects on adult schistosomes and completely inhibited egg production. ChEMBL toxicity data were also leveraged to provide further support for progressing CGP60474 (as well as luminespib and TAE684) as a novel anti-schistosomal compound. As very few compounds are currently at the advanced stages of the anti-schistosomal pipeline, our approaches highlight a strategy by which new chemical matter can be identified and quickly progressed through preclinical development.

A drug repurposing screen for whipworms informed by comparative genomics.

Hundreds of millions of people worldwide are infected with the whipworm Trichuris trichiura. Novel treatments are urgently needed as current drugs, such as albendazole, have relatively low efficacy. We have investigated whether drugs approved for other human diseases could be repurposed as novel anti-whipworm drugs. In a previous comparative genomics analysis, we identified 409 drugs approved for human use that we predicted to target parasitic worm proteins. Here we tested these ex vivo by assessing motility of adult worms of Trichuris muris, the murine whipworm, an established model for human whipworm research. We identified 14 compounds with EC50 values of ≤50 µM against T. muris ex vivo, and selected nine for testing in vivo. However, the best worm burden reduction seen in mice was just 19%. The high number of ex vivo hits against T. muris shows that we were successful at predicting parasite proteins that could be targeted by approved drugs. In contrast, the low efficacy of these compounds in mice suggest challenges due to their chemical properties (e.g. lipophilicity, polarity, molecular weight) and pharmacokinetics (e.g. absorption, distribution, metabolism, and excretion) that may (i) promote absorption by the host gastrointestinal tract, thereby reducing availability to the worms embedded in the large intestine, and/or (ii) restrict drug uptake by the worms. This indicates that identifying structural analogues that have reduced absorption by the host, and increased uptake by worms, may be necessary for successful drug development against whipworms.

Identification of anti-schistosomal, anthelmintic and anti-parasitic compounds curated and text-mined from the scientific literature.

More than a billion people are infected with parasitic worms, including nematodes, such as hookworms, and flatworms, such as blood flukes. Few drugs are available to treat worm infections, but high-throughput screening approaches hold promise to identify novel drug candidates. One problem for researchers who find an interesting 'hit' from a high-throughput screen is to identify whether that compound, or a similar compound has previously been published as having anthelmintic or anti-parasitic activity. Here, we present (i) data sets of 2,828 anthelmintic compounds, and 1,269 specific anti-schistosomal compounds, manually curated from scientific papers and books, and (ii) a data set of 24,335 potential anthelmintic and anti-parasitic compounds identified by text-mining PubMed abstracts. We provide their structures in simplified molecular-input line-entry system (SMILES) format so that researchers can easily compare 'hits' from their screens to these anthelmintic compounds and anti-parasitic compounds and find previous literature on them to support/halt their progression in drug discovery pipelines.

Genome sequences and comparative genomics of two Lactobacillus ruminis strains from the bovine and human intestinal tracts.

BACKGROUND: The genus Lactobacillus is characterized by an extraordinary degree of phenotypic and genotypic diversity, which recent genomic analyses have further highlighted. However, the choice of species for sequencing has been non-random and unequal in distribution, with only a single representative genome from the L. salivarius clade available to date. Furthermore, there is no data to facilitate a functional genomic analysis of motility in the lactobacilli, a trait that is restricted to the L. salivarius clade. RESULTS: The 2.06 Mb genome of the bovine isolate Lactobacillus ruminis ATCC 27782 comprises a single circular chromosome, and has a G+C content of 44.4%. In silico analysis identified 1901 coding sequences, including genes for a pediocin-like bacteriocin, a single large exopolysaccharide-related cluster, two sortase enzymes, two CRISPR loci and numerous IS elements and pseudogenes. A cluster of genes related to a putative pilin was identified, and shown to be transcribed in vitro. A high quality draft assembly of the genome of a second L. ruminis strain, ATCC 25644 isolated from humans, suggested a slightly larger genome of 2.138 Mb, that exhibited a high degree of synteny with the ATCC 27782 genome. In contrast, comparative analysis of L. ruminis and L. salivarius identified a lack of long-range synteny between these closely related species. Comparison of the L. salivarius clade core proteins with those of nine other Lactobacillus species distributed across 4 major phylogenetic groups identified the set of shared proteins, and proteins unique to each group. CONCLUSIONS: The genome of L. ruminis provides a comparative tool for directing functional analyses of other members of the L. salivarius clade, and it increases understanding of the divergence of this distinct Lactobacillus lineage from other commensal lactobacilli. The genome sequence provides a definitive resource to facilitate investigation of the genetics, biochemistry and host interactions of these motile intestinal lactobacilli.

Caenorhabditis evolution: if they all look alike, you aren't looking hard enough.

Caenorhabditis elegans is widely known as a model organism for cell, molecular, developmental and neural biology, but it is also being used for evolutionary studies. A recent meeting of researchers in Portugal covered topics as diverse as phylogenetics, genetic mapping of quantitative and qualitative intraspecific variation, evolutionary developmental biology and population genetics. Here, we summarize the main findings of the meeting, which marks the formal birth of a research community dedicated to Caenorhabditis species evolution.

TreeFam: 2008 Update.

TreeFam (http://www.treefam.org) was developed to provide curated phylogenetic trees for all animal gene families, as well as orthologue and paralogue assignments. Release 4.0 of TreeFam contains curated trees for 1314 families and automatically generated trees for another 14,351 families. We have expanded TreeFam to include 25 fully sequenced animal genomes, as well as four genomes from plant and fungal outgroup species. We have also introduced more accurate approaches for automatically grouping genes into families, for building phylogenetic trees, and for inferring orthologues and paralogues. The user interface for viewing phylogenetic trees and family information has been improved. Furthermore, a new perl API lets users easily extract data from the TreeFam mysql database.

Quorum sensing N-Acyl homoserine lactones are a new class of anti-schistosomal.

BACKGROUND: Schistosomiasis is a prevalent neglected tropical disease that affects approximately 300 million people worldwide. Its treatment is through a single class chemotherapy, praziquantel. Concerns surrounding the emergence of praziquantel insensitivity have led to a need for developing novel anthelmintics. METHODOLOGY/PRINCIPLE FINDINGS: Through evaluating and screening fourteen compounds (initially developed for anti-cancer and anti-viral projects) against Schistosoma mansoni, one of three species responsible for most cases of human schistosomiasis, a racemic N-acyl homoserine (1) demonstrated good efficacy against all intra mammalian lifecycle stages including schistosomula (EC50 = 4.7 µM), juvenile worms (EC50 = 4.3 µM) and adult worms (EC50 = 8.3 µM). To begin exploring structural activity relationships, a further 8 analogues of this compound were generated, including individual (R)- and (S)- enantiomers. Upon anti-schistosomal screening of these analogues, the (R)- enantiomer retained activity, whereas the (S)- lost activity. Furthermore, modification of the lactone ring to a thiolactone ring (3) improved potency against schistosomula (EC50 = 2.1 µM), juvenile worms (EC50 = 0.5 µM) and adult worms (EC50 = 4.8 µM). As the effective racemic parent compound is structurally similar to quorum sensing signaling peptides used by bacteria, further evaluation of its effect (along with its stereoisomers and the thiolactone analogues) against Gram+ (Staphylococcus aureus) and Gram- (Escherichia coli) species was conducted. While some activity was observed against both Gram+ and Gram- bacteria species for the racemic compound 1 (MIC 125 mg/L), the (R) stereoisomer had better activity (125 mg/L) than the (S) (>125mg/L). However, the greatest antimicrobial activity (MIC 31.25 mg/L against S. aureus) was observed for the thiolactone containing analogue (3). CONCLUSION/SIGNIFICANCE: To the best of our knowledge, this is the first demonstration that N-Acyl homoserines exhibit anthelmintic activities. Furthermore, their additional action on Gram+ bacteria opens a new avenue for exploring these molecules more broadly as part of future anti-infective initiatives.

Large CRISPR-Cas-induced deletions in the oxamniquine resistance locus of the human parasite Schistosoma mansoni.

Background. At least 250 million people worldwide suffer from schistosomiasis, caused by Schistosoma worms. Genome sequences for several Schistosoma species are available, including a high-quality annotated reference for Schistosoma mansoni. There is a pressing need to develop a reliable functional toolkit to translate these data into new biological insights and targets for intervention. CRISPR-Cas9 was recently demonstrated for the first time in S. mansoni, to produce somatic mutations in the omega-1 ( ω1) gene. Methods. We employed CRISPR-Cas9 to introduce somatic mutations in a second gene, SULT-OR, a sulfotransferase expressed in the parasitic stages of S. mansoni, in which mutations confer resistance to the drug oxamniquine. A 262-bp PCR product spanning the region targeted by the gRNA against SULT-OR was amplified, and mutations identified in it by high-throughput sequencing. Results. We found that 0.3-2.0% of aligned reads from CRISPR-Cas9-treated adult worms showed deletions spanning the predicted Cas9 cut site, compared to 0.1-0.2% for sporocysts, while deletions were extremely rare in eggs. The most common deletion observed in adults and sporocysts was a 34 bp-deletion directly upstream of the predicted cut site, but rarer deletions reaching as far as 102 bp upstream of the cut site were also detected. The CRISPR-Cas9-induced deletions, if homozygous, are predicted to cause resistance to oxamniquine by producing frameshifts, ablating SULT-OR transcription, or leading to mRNA degradation via the nonsense-mediated mRNA decay pathway. However, no SULT-OR knock down at the mRNA level was observed, presumably because the cells in which CRISPR-Cas9 did induce mutations represented a small fraction of all cells expressing SULT-OR. Conclusions. Further optimisation of CRISPR-Cas protocols for different developmental stages and particular cell types, including germline cells, will contribute to the generation of a homozygous knock-out in any gene of interest, and in particular the SULT-OR gene to derive an oxamniquine-resistant stable transgenic line.

Large-scale RNAi screening uncovers therapeutic targets in the parasite Schistosoma mansoni.

Schistosome parasites kill 250,000 people every year. Treatment of schistosomiasis relies on the drug praziquantel. Unfortunately, a scarcity of molecular tools has hindered the discovery of new drug targets. Here, we describe a large-scale RNA interference (RNAi) screen in adult Schistosoma mansoni that examined the function of 2216 genes. We identified 261 genes with phenotypes affecting neuromuscular function, tissue integrity, stem cell maintenance, and parasite survival. Leveraging these data, we prioritized compounds with activity against the parasites and uncovered a pair of protein kinases (TAO and STK25) that cooperate to maintain muscle-specific messenger RNA transcription. Loss of either of these kinases results in paralysis and worm death in a mammalian host. These studies may help expedite therapeutic development and invigorate studies of these neglected parasites.

Single-cell atlas of the first intra-mammalian developmental stage of the human parasite Schistosoma mansoni.

Over 250 million people suffer from schistosomiasis, a tropical disease caused by parasitic flatworms known as schistosomes. Humans become infected by free-swimming, water-borne larvae, which penetrate the skin. The earliest intra-mammalian stage, called the schistosomulum, undergoes a series of developmental transitions. These changes are critical for the parasite to adapt to its new environment as it navigates through host tissues to reach its niche, where it will grow to reproductive maturity. Unravelling the mechanisms that drive intra-mammalian development requires knowledge of the spatial organisation and transcriptional dynamics of different cell types that comprise the schistomulum body. To fill these important knowledge gaps, we perform single-cell RNA sequencing on two-day old schistosomula of Schistosoma mansoni. We identify likely gene expression profiles for muscle, nervous system, tegument, oesophageal gland, parenchymal/primordial gut cells, and stem cells. In addition, we validate cell markers for all these clusters by in situ hybridisation in schistosomula and adult parasites. Taken together, this study provides a comprehensive cell-type atlas for the early intra-mammalian stage of this devastating metazoan parasite.

Programmed genome editing of the omega-1 ribonuclease of the blood fluke, Schistosoma mansoni.

CRISPR/Cas9-based genome editing has yet to be reported in species of the Platyhelminthes. We tested this approach by targeting omega-1 (ω1) of Schistosoma mansoni as proof of principle. This secreted ribonuclease is crucial for Th2 polarization and granuloma formation. Schistosome eggs were exposed to Cas9 complexed with guide RNA complementary to ω1 by electroporation or by transduction with lentiviral particles. Some eggs were also transfected with a single stranded donor template. Sequences of amplicons from gene-edited parasites exhibited Cas9-catalyzed mutations including homology directed repaired alleles, and other analyses revealed depletion of ω1 transcripts and the ribonuclease. Gene-edited eggs failed to polarize Th2 cytokine responses in macrophage/T-cell co-cultures, while the volume of pulmonary granulomas surrounding ω1-mutated eggs following tail-vein injection into mice was vastly reduced. Knock-out of ω1 and the diminished levels of these cytokines following exposure showcase the novel application of programmed gene editing for functional genomics in schistosomes.

nGASP--the nematode genome annotation assessment project.

BACKGROUND: While the C. elegans genome is extensively annotated, relatively little information is available for other Caenorhabditis species. The nematode genome annotation assessment project (nGASP) was launched to objectively assess the accuracy of protein-coding gene prediction software in C. elegans, and to apply this knowledge to the annotation of the genomes of four additional Caenorhabditis species and other nematodes. Seventeen groups worldwide participated in nGASP, and submitted 47 prediction sets across 10 Mb of the C. elegans genome. Predictions were compared to reference gene sets consisting of confirmed or manually curated gene models from WormBase. RESULTS: The most accurate gene-finders were 'combiner' algorithms, which made use of transcript- and protein-alignments and multi-genome alignments, as well as gene predictions from other gene-finders. Gene-finders that used alignments of ESTs, mRNAs and proteins came in second. There was a tie for third place between gene-finders that used multi-genome alignments and ab initio gene-finders. The median gene level sensitivity of combiners was 78% and their specificity was 42%, which is nearly the same accuracy reported for combiners in the human genome. C. elegans genes with exons of unusual hexamer content, as well as those with unusually many exons, short exons, long introns, a weak translation start signal, weak splice sites, or poorly conserved orthologs posed the greatest difficulty for gene-finders. CONCLUSION: This experiment establishes a baseline of gene prediction accuracy in Caenorhabditis genomes, and has guided the choice of gene-finders for the annotation of newly sequenced genomes of Caenorhabditis and other nematode species. We have created new gene sets for C. briggsae, C. remanei, C. brenneri, C. japonica, and Brugia malayi using some of the best-performing gene-finders.

Chromosome evolution in eukaryotes: a multi-kingdom perspective.

In eukaryotes, chromosomal rearrangements, such as inversions, translocations and duplications, are common and range from part of a gene to hundreds of genes. Lineage-specific patterns are also seen: translocations are rare in dipteran flies, and angiosperm genomes seem prone to polyploidization. In most eukaryotes, there is a strong association between rearrangement breakpoints and repeat sequences. Current data suggest that some repeats promoted rearrangements via non-allelic homologous recombination, for others the association might not be causal but reflects the instability of particular genomic regions. Rearrangement polymorphisms in eukaryotes are correlated with phenotypic differences, so are thought to confer varying fitness in different habitats. Some seem to be under positive selection because they either trap favorable allele combinations together or alter the expression of nearby genes. There is little evidence that chromosomal rearrangements cause speciation, but they probably intensify reproductive isolation between species that have formed by another route.

Genomix: a method for combining gene-finders' predictions, which uses evolutionary conservation of sequence and intron-exon structure.

MOTIVATION: Correct gene predictions are crucial for most analyses of genomes. However, in the absence of transcript data, gene prediction is still challenging. One way to improve gene-finding accuracy in such genomes is to combine the exons predicted by several gene-finders, so that gene-finders that make uncorrelated errors can correct each other. RESULTS: We present a method for combining gene-finders called Genomix. Genomix selects the predicted exons that are best conserved within and/or between species in terms of sequence and intron-exon structure, and combines them into a gene structure. Genomix was used to combine predictions from four gene-finders for Caenorhabditis elegans, by selecting the predicted exons that are best conserved with C.briggsae and C.remanei. On a set of approximately 1500 confirmed C.elegans genes, Genomix increased the exon-level specificity by 10.1% and sensitivity by 2.7% compared to the best input gene-finder. AVAILABILITY: Scripts and Supplementary Material can be found at http://www.sanger.ac.uk/Software/analysis/genomix

TreeFam: a curated database of phylogenetic trees of animal gene families.

TreeFam is a database of phylogenetic trees of gene families found in animals. It aims to develop a curated resource that presents the accurate evolutionary history of all animal gene families, as well as reliable ortholog and paralog assignments. Curated families are being added progressively, based on seed alignments and trees in a similar fashion to Pfam. Release 1.1 of TreeFam contains curated trees for 690 families and automatically generated trees for another 11 646 families. These represent over 128 000 genes from nine fully sequenced animal genomes and over 45 000 other animal proteins from UniProt; approximately 40-85% of proteins encoded in the fully sequenced animal genomes are included in TreeFam. TreeFam is freely available at http://www.treefam.org and http://treefam.genomics.org.cn.

The genome sequence of Caenorhabditis briggsae: a platform for comparative genomics.

The soil nematodes Caenorhabditis briggsae and Caenorhabditis elegans diverged from a common ancestor roughly 100 million years ago and yet are almost indistinguishable by eye. They have the same chromosome number and genome sizes, and they occupy the same ecological niche. To explore the basis for this striking conservation of structure and function, we have sequenced the C. briggsae genome to a high-quality draft stage and compared it to the finished C. elegans sequence. We predict approximately 19,500 protein-coding genes in the C. briggsae genome, roughly the same as in C. elegans. Of these, 12,200 have clear C. elegans orthologs, a further 6,500 have one or more clearly detectable C. elegans homologs, and approximately 800 C. briggsae genes have no detectable matches in C. elegans. Almost all of the noncoding RNAs (ncRNAs) known are shared between the two species. The two genomes exhibit extensive colinearity, and the rate of divergence appears to be higher in the chromosomal arms than in the centers. Operons, a distinctive feature of C. elegans, are highly conserved in C. briggsae, with the arrangement of genes being preserved in 96% of cases. The difference in size between the C. briggsae (estimated at approximately 104 Mbp) and C. elegans (100.3 Mbp) genomes is almost entirely due to repetitive sequence, which accounts for 22.4% of the C. briggsae genome in contrast to 16.5% of the C. elegans genome. Few, if any, repeat families are shared, suggesting that most were acquired after the two species diverged or are undergoing rapid evolution. Coclustering the C. elegans and C. briggsae proteins reveals 2,169 protein families of two or more members. Most of these are shared between the two species, but some appear to be expanding or contracting, and there seem to be as many as several hundred novel C. briggsae gene families. The C. briggsae draft sequence will greatly improve the annotation of the C. elegans genome. Based on similarity to C. briggsae, we found strong evidence for 1,300 new C. elegans genes. In addition, comparisons of the two genomes will help to understand the evolutionary forces that mold nematode genomes.

Comparative genomics in C. elegans, C. briggsae, and other Caenorhabditis species.

The genome of the nematode Caenorhabditis elegans was the first animal genome sequenced. Subsequent sequencing of the Caenorhabditis briggsae genome enabled a comparison of the genomes of two nematode species. In this chapter, we describe the methods that we used to compare the C. elegans genome to that of C. briggsae. We discuss how these methods could be developed to compare the C. elegans and C. briggsae genomes to those of Caenorhabditis remanei, C. n. sp. represented by strains PB2801 and CB5161, among others (1), and Caenorhabditis japonica, which are currently being sequenced.