- Open Access
Gene expression patterns during adaptation of a helminth parasite to different environmental niches
Genome Biology volume 8, Article number: R65 (2007)
Schistosome bloodflukes are complex trematodes responsible for 200 million cases of schistosomiasis worldwide. Their life cycle is characterized by a series of remarkable morphological and biochemical transitions between an invertebrate host, an aquatic environment, and a mammalian host. We report a global transcriptional analysis of how this parasite alters gene regulation to adapt to three distinct environments.
Utilizing a genomic microarray made of 12,000 45-50-mer oligonucleotides based on expressed sequence tags, three different developmental stages of the schistosome parasite were analyzed by pair-wise comparisons of transcript hybridization signals. This analysis resulted in the identification of 1,154 developmentally enriched transcripts.
This study expands the repertoire of schistosome genes analyzed for stage-specific expression to over 70% of the predicted genome. Among the new associations identified are the roles of robust protein synthesis and programmed cell death in development of cercariae in the sporocyst stages, the relative paucity of cercarial gene expression outside of energy production, and the remarkable diversity of adult gene expression programs that reflect adaptation to the host bloodstream and an average lifespan that may approach 10 years.
Schistosomiasis is a chronic debilitating parasitic disease affecting some 200 million people across 74 countries within Africa, Asia, the Middle East and South America. In terms of public health and socio-economic impact, it ranks second only to malaria among parasitic diseases [1, 2]. The causative agents of schistosomiasis are schistosome bloodflukes, multicellular trematodes whose life cycle is characterized by a series of striking morphological and biochemical transitions between an intermediate host snail in an aquatic environment, two free-swimming aquatic larval forms, and a warm-blooded mammalian host (Figure 1). As such, the schistosome represents an ideal but challenging biological system in which to identify programs of gene regulation that have evolved to facilitate adaptation of metazoa to different biological microenvironments.
In the intermediate host aquatic snail, sac-like clusters of differentiating larvae called daughter sporocysts give rise to large numbers of free-swimming aquatic cercariae. Cercariae navigate their environment by a side to side motion of their tails and find the mammalian host through chemical and environmental signals [3–6]. Upon entering the mammalian host, the tail is detached, and an early vascular stage called the schistosomulum rapidly readjusts to the new host environment. During this time, the schistosomulum undergoes changes in basic metabolism and surface properties . The schistosomulum navigates the vascular system, including transiting through at least two capillary plexuses. Between 14 and 28 days, the organisms undergo rapid growth and differentiation to adult male and female forms in specific regions of the host vasculature. Male and female worms form mating-pairs and the female produces hundreds to thousands of eggs per day, depending on the species of schistosome. The eggs, in turn, pass out of the host via feces or urine, and hatch in fresh water to yet another morphologically distinct form, the miracidium. This larval form moves through fresh water by means of numerous cilia, and invades the soft parts of the intermediate host snail, transforming into the mother sporocyst. The cycle is complete as the mother sporocyst produces daughter sporocysts, which serve as the developmental stratum for production of cercariae.
Efforts to reduce schistosome prevalence have included combinations of sanitation, health education, snail control, better diagnosis and chemotherapy . Control of parasitemia has relied primarily on praziquantel, an effective chemotherapeutic drug that has been in use for over 20 years. Recurring morbidity and evidence of emerging resistance to praziquantel in endemic countries emphasizes the need for both an effective vaccine and new chemotherapy [9–11]. Several potential vaccine candidates have been identified [12, 13], but to date none is sufficiently effective for practical use. Given the complexity of the schistosome lifecycle and its ability to evade the immune system, a more comprehensive approach to understanding the molecular biology and biochemistry of schistosomes is essential to rationally identify and develop potential vaccine candidates or new drugs.
The approximately 270 Mb genome of Schistosoma mansoni  is organized into 8 chromosomes , including the heterologous female pair WZ . It is estimated to have 30% GC content . Currently, the S. mansoni genome is thought to contain some 14,000 predicted genes [18, 19]. Until recently, there have been relatively few attempts to analyze the functional genomics of schistosomes due to limited genomic sequence information. Identification of genes expressed in a sex-specific manner by adult worms was made using either a 576, 7,335, or 4,608 oligonucleotide cDNA based microarray system [20–22] and an analysis of genes associated with pairing of adult male and female worms was addressed . The transition between miracidium and mother sporocyst was also analyzed using the same 7,335 component oligonucleotide array noted above . Comparisons of adult transcripts of two schistosome species, and comparison of lung stage versus experimentally produced schistosomula, were carried out using a larger array [7, 25]. A comparative study between the S. japonicum proteome and transcriptome was also undertaken .
Now that the genome is nearing completion, we present a comprehensive analysis of gene expression during three of the major stages of S. mansoni development: daughter sporocyst, cercaria, and adult. These stages represent key transitions from intermediate host snail to aquatic environment to mammalian host. We designed a 12,000 oligonucleotide microarray chip made of 45-50-mer oligonucleotides and we analyzed expression of a majority of the predicted 14,000 S. mansoni genes under very stringent conditions. While the present annotation of the genome does not distinguish gene insert orientation, we found that over 9,700 of the oligonucleotides printed yielded clear transcriptional signals in at least one parasite stage. We have examined pair-wise differential expression at each stage and identify those genes whose RNA expression profile is conserved or differentially expressed across developmental stages. Noteworthy gene clusters that support previous hypotheses or provide new insights into the responses of the parasite to environmental transitions during the lifecycle are identified. Highly expressed transcripts in sporocysts include those involved in general protein synthesis and quality control, consistent with the function of sporocysts in production of large numbers of cercariae. Cercarial transcripts are dominated by genes involved in mitochondrial function, supporting the energy production necessary for cercarial swimming. Adult worms express a diverse pattern of transcripts necessary for egg production, energy metabolism, immune evasion and physiological maintenance of a relatively long-lived organism.
The 12,000 element schistosome array
The microarray used in this study contained 12,000 individual 45-50-mer oligonucleotides based on 11,998 tentative consensus sequences (TCs), as documented by the Schistosoma mansoni Genome Index maintained at The Institute for Genomic Research (TIGR) . TCs are created by virtual assemblage of full or partial cDNA sequences into transcripts. Spliced forms of genes are listed separately. Of the 12,912 sequences provided by TIGR, 11,998 TCs were chosen based on the maximum size of the TC available. Thus, this microarray contains a majority of known sequences available in the S. mansoni genome. There have been concerns raised about the orientation of inserts in the current annotation of the TIGR assemblage; however, we found that >9,700 of the 12,000 oligonucleotides printed gave some transcript signals in one or another parasite stage by type II analysis (see Materials and methods) and, as detailed below, many key gene programs could be identified and correlated with environmental transitions.
The microarray chip is specific for S. mansonitranscripts
Extraction of RNA from schistosome daughter sporocysts involves excision of the entire snail hepatopancreas, and co-lysis of the snail and schistosome material. To control for any background hybridization of snail-specific material (to our S. mansoni cDNA-based 12,000-oligonucleotide microarray), we hybridized hepatopancreas RNA from infected snails containing daughter sporocysts versus uninfected snail hepatopancreas RNA across four different chips. With the exception of two genes, we found that intermediate host snail RNA was not detected (data not shown). Even for these two genes (TC8129, similar to beta-1 4-galactosyltransferase, and TC6896, similar to Unknown), hybridization signals were not seen on every run and were detected only by type II analysis. In summary, no significant contribution of snail RNA to the hybridization analysis was found. Daughter sporocysts contain cercariae at various stages of differentiation. Nevertheless, the transcriptome profile of this stage in the snail differed significantly from that of mature aquatic cercariae.
We also found that microarray samples and chips were reproducible and correlative. There was an average Pierson correlation coefficient (log ratio of medians) of 0.7 for biological replicates and 0.96 for technical replicates. Dye switch experiments showed that data differences were not due to dye-labeling bias (Additional data file 1).
S. mansonitranscriptome overview
This clustering analysis includes 431 genes analyzed by pair-wise comparison with at Examination of all developmental stages showed a transcript pattern specific to each. Figure 2 (Table 1) is a clustering analysis of the major transcripts identified in which duplicates of sporocyst (Cy5) were compared to cercariae (Cy3). For example, cluster 1 represents transcripts that are upregulated in the daughter sporocyst stage but not significantly in adults or cercariae. Cluster 2 shows transcripts up-regulated in both sporocysts and adults relative to cercariae. Clusters 3 and 4 show transcripts enriched in cercariae relative to adults. Clusters 5 and 6 show transcripts that are upregulated in adult worms but are either relatively low or not enriched in cercariae or sporocysts. least a 3-fold difference in transcript level between developmental stages. Overall, we identified enrichment of 1,154 genes during the schistosome developmental life cycle, the largest proportion of which is expressed in adult worms (type II analysis). Of the 1,154 genes, 406 represent genes of unknown function (Figure 3).
Validation of microarray data by real time PCR
To test the validity of the microarray results, we performed a quantitative analysis study on the expression levels of several parasite genes by real-time PCR. Among the transcripts analyzed were actin, cathepsin B, cytochrome C, eggshell protein precursor, Sm23 and cathepsin L. Gene-specific primers were designed to validate the gene expression profile for sporocyst, cercariae, and adult developmental stages (Additional data file 2). To ascertain gene expression transitions between cercariae and adults, we also analyzed transcript levels in 24-hour schistosomula, the stage of schistosomal development between cercarial penetration and adult worms. We found that the expression of genes analyzed by RT-PCR correlates with our microarray analysis in 10 of 11 cases (Figure 4).
Genes expressed in daughter sporocysts reflect production of cercariae
The daughter sporocyst stage resides in the intermediate host snail. The key biological function of this stage is to support the differentiation and development of large numbers of cercariae, the aquatic larval stage that will initiate invasion of the mammalian host. The function of this stage in generating cercariae is underscored by the fact that many of the most highly expressed transcripts are gene products that function in general protein synthesis (40S and 60S ribosomal subunits, elongation factor) and post-translational protein folding and fidelity (chaperones, ubiquitin). Furthermore, many of the proteins that will be utilized by cercariae in the initial stages of host skin invasion are also produced in this stage, including cercarial elastase (aka cercarial protease), cercarial muscle proteins (actin, dynein light chain isoforms), and calcium binding proteins that are also abundant in the proteome of cercarial secretions .
As would be expected from the large-scale differentiation of cercariae within the daughter sporocyst, transcription factors such as the Y-box binding protein and stathmin, a phosphoryl protein involved in vertebrate growth and regulation of differentiation  are highly expressed. Sporocysts express a cathepsin L-like cysteine protease homolog (TC9217) that is not significantly enriched in adults (Figure 4f). In adults a cathepsin L is gut localized and functions in digestion . The lack of significant enrichment of other gut-specific genes in sporocysts, such as those encoding cathepsins B, B1, and C, and the unique profile of the cathepsin-like protease homolog, suggests that this protease isoform may function differently in the daughter sporocyst, possibly in the generation of cercarial progeny. In Caenorhabditis elegans, a cathepsin L is essential for embryogenesis and development .
Genes expressed in cercariae reflect energy production and motility
Cercariae are a relatively short-lived 'transitional' stage that are released from snails into fresh water and must swim to find, and ultimately invade, a mammalian host. Cercariae will not survive if they fail to enter the mammalian host before energy sources are exhausted. Consistent with this concept of the biological function of cercariae, transcripts that are upregulated are primarily factors necessary to sustain swimming behavior and invasion (Figure 2, clusters 3 and 4). The highest expressed transcripts are genes involved in mitochondrial function or energy metabolism. These include NADH dehydrogenase and its various subunits, cytochrome C and its homologues, and ATP/ADP carrier proteins. Other transcripts highly expressed in cercariae include structural and motility genes like actin and fibrillin and transcripts coding for a protease that plays a role in host invasion, cercarial elastase . In general, fewer transcripts are found in cercariae relative to other stages.
Genes expressed in adult worms
Compared to sporocysts and cercariae, the most abundant transcripts in adult worms (male and female worm pairs) are significantly more diverse. This reflects the fact that adult worms must evade the mammalian host immune system, maintain motility, acquire and metabolize a variety of nutrients, form mating pairs, and produce large numbers of eggs. In parallel with these functions, highly abundant transcripts in adult worms include a group of genes involved in protein degradation in the gut (cathepsin B, cathepsin L), egg production, and oxidative stress responses [33–36]. Genes involved in gonadal differentiation for both male and female worms, as well as genes coding for proteins for egg-associated proteins are readily identified. These include the anti-mullerian hormone receptor for males, and members of the RXR family and Smad family, which are cell signaling pathways previously associated with function of female gonads [37, 38]. Highly expressed egg transcripts include eggshell protein precursor, major egg antigen, and several homologues of these proteins.
Perhaps most striking is the abundance of transcripts in adult worms from genes coding for surface proteins, including nutrient transporters such as the glucose transporter and a number of surface proteins of unknown function previously identified in research projects aimed at subunit vaccine development. These include Sm23 [39–45], integral membrane protein 25 , Sm14 [47, 48], and 26 and 28 glutathione S-transferases [49–52].
The stage-specific transcriptome of S. mansoni provides an informatics foundation for the study of parasite gene regulation and a correlative for proteomic studies. This microarray study validates and extends observations made with a 7,335-oligonucleotide array from previously available expressed sequence tags (ESTs) , and correlates well with comparative analyses between the transcriptome and proteome of adult Schistosoma japonicum and S. mansoni . With 12,000 oligonucleotides, a more complete picture of gene programs that mark transitions between key schistosome stages in distinct environmental niches is now apparent. The entire schistosome genome is estimated to have 14,000 functional genes. While it has been noted that the TIGR annotation to date has not distinguished the orientation of the inserts sequenced, we found that 9,700 of the 12,000 genes analyzed on this array gave clear positive signals for at least one stage. More importantly, many genes could be functionally annotated and gene programs correlating with environmental transitions of the helminth parasite discovered or validated.
One of the most remarkable aspects of the schistosome parasite is the dramatic morphological change that takes place between life cycle stages (Figure 1). These morphological changes parallel transition of the parasite between three markedly different environments. The daughter sporocyst stage is the late developmental stage in the intermediate host snail - a mollusk of fresh water habitats. Motile cercariae emerge from the snail and must navigate an aquatic environment to find the mammalian host. Following entry into the skin of the mammalian host, the parasite must now adapt to a different osmotic environment and a warm-blooded host. Finally, within the mammalian host, two distinct sexual stages develop and eggs are subsequently produced.
Array analysis comparing daughter sporocysts to the cercarial stage highlighted gene programs responsible for supporting larval development. The sporocyst stage is a very efficient protein synthesis factory. This correlates with the need for mass production of cercariae and the availability of nutrients from an intimate relationship with the snail host. In keeping with the complex morphological differentiation of cercariae with distinct body parts and 'organs' from an embryonic cell mass, daughter sporocyst transcripts included genes involved not only in general protein synthesis such as ribosomal and heat shock genes, but also transcripts annotated to function in programmed cell death (TC18024) and ubiquitination, a process key to morphological differentiation.
In Figure 2, cluster 2, daughter sporocysts share with adult worms a previously unrecognized requirement for transcriptional upregulation of antioxidant genes such as glutathione-S-transferase, thioredoxin, and thioredoxin peroxidase. This may reflect oxidative stress generated within the immediate snail host, or storage of these proteins in developing cercariae, so that they are available in the earliest stages of mammalian host invasion. Utilizing a 7,335 oligonucleotide array, Vermeire et al.  documented gene expression patterns between miracidia, the aquatic stage that invades the snail host, and the mother sporocyst stage, which is the initial stage following invasion of the snail host. As is the case for daughter sporocysts analyzed in our study, they also found upregulation of several genes involved in protein synthesis, the redox pathway, and proteolysis. This suggests that these gene programs are initiated in mother sporocysts following entry of miracidia into the snail, and sustained throughout the daughter sporocyst stage until mature cercariae leave the snail.
Transition to the cercarial stage is marked by a reduction in transcript level of a large number of genes as seen in clusters 1 and 2 of Figure 2. We found 116 genes whose abundance is reduced in cercariae by 3-fold relative to sporocysts (type II analysis where sporocyst intensity units are set >3,000). In comparison to other stages, the cercariae are less transcriptionally active. Fewer than 7 genes are upregulated more than 2-fold relative to adults and only 34 are upregulated when compared pair-wise to adult transcripts.
The cercariae express a cluster of genes consistent with the energy required to move rapidly through water in search of a mammalian host. These include genes functioning in ATP production and utilization, presumably for muscle function and swimming behavior. Consistent with previous Northern blot analyses and biochemical studies, much of the repertoire of proteins that cercariae use to invade the skin of their host, as well as structural proteins, have already been produced during cercarial development in the daughter sporocyst stage within the nutrient rich intermediate host snail (Figure 2a-d). Cercariae do express caspases and related cell death programs, which may be required for morphological remodeling during transition into the schistosomulum stage. Cercariae detach their motile tail and surface glycocalyx shortly after entry into the mammalian host. Large gland structures (acetabular glands) producing invasive proteases involute within 48-72 hours of invasion.
Residence of adult schistosomes in the mammalian host bloodstream is supported by expression of gene families that respond to oxidative stress and genes involved in adjustment and adaptation to a new osmotic environment (aquaporin) [54–56]. Adult worms have a major requirement to digest blood-proteins and acquire other nutrients (cathepsin B, glucose transporter) [33, 35, 57]. Schistosomes differentiate into male and female worms and expression of an anti-mullerian hormone receptor [58, 59] is likely related to sexual differentiation. Female parasites produce hundreds to thousands of eggs per day as reflected in the upregulation of egg shell proteins. Recent analysis of gene expression in earlier intravascular stages can now be compared with the adult (versus cercariae) transcriptome presented here. Using a cDNA array, Dillon et al.  equated seven-day cultured schistosomula with lung schistosomula and compared gene expression across life stages. Chai et al.  utilized schistosomula directly obtained from the lungs of infected mice, and compared transcripts expressed to adult worms, cercariae, and newly transformed schistosomula. By and large, the genes we found expressed in adult parasites versus cercariae were similar to those observed in the comparison of lung worms versus cercariae. These included nutrient acquisition genes, such as those encoding the glucose transporter and the proteolytic cathepsins. The genes down-regulated in adult worms relative to lung schistosomula include several we found expressed in cercariae or daughter sporocysts, including those encoding the anti-inflammatory protein Sm16 and paramyosin.
The intravascular stages of schistosome parasites have a complex and highly adapted relationship with the mammalian host [61, 62]. To support this relationship, intravascular stages upregulate surface proteins or receptors, some of which have homology to mammalian receptors and/or factors involved in signaling cascades (Additional data file 6). These include the anti-Mullerian hormone type II receptor  and the thyroid receptor interacting protein .
Developmental regulation of cell number and type by programmed cell death appears to be an important function in all stages of schistosome development analyzed, as it is in the nematode C. elegans . Homologues for genes involved in programmed cell death include TC11294 and the DAP-1 homolog TC18024 (type II analysis).
Since the initial stages of the genome analysis of schistosomes, it has been clear that transposon-like sequences are common. Efforts to exploit these elements for genetic manipulation are ongoing [66, 67]. It is noteworthy that one transposon (TC17720) is more abundantly expressed in the sporocyst than in any other developmental stage. It is also expressed in adults, albeit at low levels, as corroborated in the study by Gobert et al. . There is also another retroelement (TC7011) highly enriched in adult worms.
This expanded array data set, focusing on life cycle stage transitions, should aid in current attempts to develop transfection and gene knockout studies for schistosomes by identifying those genes that are stage-specific versus others shared among different developmental states. For example, the gene encoding cytochrome C is expressed at comparable levels in all stages studied. In contrast, genes such as those encoding the cercarial elastase in the sporocyst stage, or eggshell components in the adult female stage, represent potential models for identifying and characterizing key spatial and temporal promoter elements, and ultimately molecular mechanisms of gene regulation. Fitzpatrick et al. and Moertel et al. also have studied gender-specific gene transcription in S. mansoni using a 7,335 oligonucleotide array  and a 22,575 combined S. japonicum and S. mansoni array . The results reported here for genes that can be attributed to either male or female worms largely validate these studies.
Analysis of the stage-specific transcriptional program of schistosomes also helps to validate and underscore differences in the levels of gene products noted in proteomic studies of schistosome life cycle stages. The levels of many gene products identified by proteomic analysis parallel transcription levels in this microarray analysis. In contrast to protozoan parasites like Leishmania , this may indicate that less post-transcriptional regulation operates in schistosomes. Having transcriptome data for three major stages, combined with proteomic data, should now facilitate a more focused analysis to determine to what extent RNA stability or other post-transcriptional mechanisms play a role in schistosome gene regulation. It was noted from proteome analysis that many genes are expressed across stages and these were, therefore, seen as less attractive vaccine targets [26, 70]. The identification of several stage-specific patterns of expression in this present study should help to redirect efforts aimed at finding the most logical candidates for a subunit vaccine, and also identify new targets to explore for drug therapy.
We show that the daughter sporocyst stage in the intermediate host snail functions primarily to support the development of invasive larvae (cercariae) by up-regulating expression of genes involved in protein synthesis, cellular differentiation, and programmed cell death. Many of the major structural and functional components of cercariae, utilized to later invade the mammalian host, are expressed and packaged prior to larval release from the snail. The aquatic cercariae themselves are less transcriptionally active than other stages, with an emphasis on production of proteins involved in energy metabolism and motility. This is in keeping with the function of these larvae in swimming from snail to mammalian host. Finally, adult parasites, which have adapted to survival in the mammalian host bloodstream, have a complex transcriptional program that supports adaptation to a new host temperature and chemical environment, evasion of the host immune response, acquisition of nutrients, and production of eggs for transmission to a new host.
Materials and methods
All experiments were performed with S. mansoni of the NMRI Puerto Rican strain maintained routinely through Biomphalaria glabrata snails and Syrian golden hamsters as previously detailed .
Schistosome stage collection
Daughter sporocysts were collected by dissection of whole hepatopancreas from six-week old infected B. glabrata snails, during maturation of daughter sporocysts and prior to cercarial release. Uninfected snail hepatopancreas was also collected as a control. Cercariae were collected in distilled water from infected B. glabrata snails using the light induction method as previously described [71, 72]. Following exposure to light for 2 hours, 50-60 snails shed about 200-300 cercariae/snail. Several collections were pooled and used for RNA extraction. Twenty-four hour schistosomula were mechanically transformed as previously described . Adult worms were recovered from the mesenteric veins by portal perfusion as previously described .
Total RNA was extracted from all biological samples for Combimatrix microarray analysis. Multiple uninfected and infected hepatopancreas, cercariae, and adult worms were homogenized in 1 ml Trizol (Invitrogen, Carlsbad, CA, USA) and RNA was extracted using the standard Invitrogen protocol. RNA was eluted in deionized water, quantified, and checked for RNA quality by UV spectrometry and agarose gel analysis. Total RNA (15 μg) for individual stages (coupled to Cy3 or Cy5) was hybridized to Combimatrix custom array microarray chips, as described in detail below (Combimatrix Corporation, Mukilteo, WA, USA).
S. mansoniDNA oligonucleotide probe design
All probes were designed using the Combimatrix software CombiMatrix Automated Probe-design Suite (PDS), based on 11,997 TC groups from sequences at TIGR. TCs represent a virtual assemblage of ESTs. Of the 12,717 TIGR TC groups, 11,998 were used for the 12,000 spot array based on size of the TC source cDNA. Those with the smallest TC source cDNA sequences were excluded. The 12,000 sequences were loaded into the Combimatrix Custom Array content probe array system for probe design, one of which was spotted in triplicate using different oligonucleotide designs constructed by hand. The probe design system took each TC sequence and designed 45-50-mer probes to be unique to each gene with predictable thermodynamic behaviors.
Microarray hybridization and data analysis
Combimatrix Custom Array microarray chip hybridizations were stringently performed in duplicate or triplicate. Hybridizations (3× SSC, 0.025 M HEPES, 0.28 μg/μl polyA, 0.05% SDS) occurred at 63°C for 18 hours to maintain specific binding efficacy. DNA microarray chips were scanned using an Axon4000B scanner at 5 μ resolution and 100% laser power and the images were analyzed with GenPix Pro 4 (Molecular Devices, Sunnyvale, CA, USA). Microarray data were stored in the NOMAD microarray database  and normalized by a global normalization using unflagged features with a regression correlation coefficient ≥0.75 and median intensity value >0 as previously described . Data were clustered using Cluster 3.0  and visualized in Java TreeView version 1.0.11 . Clustering was done as below: all data were log transformed and filtered for the presence of 100% of genes and a standard deviation value of log2 1.5 or 3-fold gene expression, followed by correlation uncentered clustering of genes and arrays by complete linkage analysis. A second and less stringent analysis was repeated using the same methodology, except clustering was done with a standard deviation value of log2 0.6, or 1.5-fold gene expression (Additional data file 3). Each Combimatrix array has 12,500 spots; 148 are empty controls (no oligonucleotides) and 352 represent various small oligonucleotide controls. For a less sophisticated analysis of gene expression, on or off, each channel (635 and 532) was treated independently. The data were reanalyzed and normalized as described above. The average median value of the 148 empty spots for each chip was used as background and subtracted from the initial intensity unit value to produce a 'normalized value'. Spots from the 'normalized value' (IU) were retained as expressed genes. We refer to this type of analysis as type II analysis. To compare all three stages, the average median of all cercariae median values was used as a standard (Additional data file 4). Transcripts with an IU >2,000 are defined as enhanced.
A total of 1 μg parasite RNA from each stage was used to prepare double stranded cDNA using SuperScript™ II reverse transcriptase (Invitrogen) in the presence of oligo dT. Quantitative PCR (qPCR) experiments were repeated in duplicate using cytochrome c as an internal control gene. The relative quantification of our genes of interest was carried out by mixing cDNA from different parasite stages with SYBR Green PCR master mix in 96 well plates, which was then incubated in the Applied Biosystems (Foster City, CA, USA) 7500 Real-Time PCR System. The program included an initial melting phase for 2 minutes at 50°C, denaturation and hot start for 10 minutes at 95°C, followed by 40 amplification cycles (95°C, 15 s; 60°C, 1 minutes). The sequences of the primers and names of the studied genes are given in Additional data file 2.
Additional data files
The following additional data are available with the online version of this paper. Additional data file 1 is an analysis of Cy3 and Cy5 dye-switch experiments comparing sporocyst to cercariae, and adults to cercariae. Additional data file 2 is a table of primer sequences used for real-time PCR. Additional data file 3 is a table listing transcripts enriched at least 1.5-fold relative to cercariae in adults or sporocysts. Additional data file 4 lists all genetic data analyzed by type II analysis. Additional data file 5 is a table listing the 431 genes highly enriched from the clustering analysis shown in Figure 2. Additional data file 6 is a table listing S. mansoni receptor-like proteins and log2 gene expression ratios. Additional data file 7 lists all oligonucleotide sequences used for the schistosome array. Additional data file 8 is a key to compare the schistosome TC sequence to the updated TIGR schistosome gene index (Version 6).
Chitsulo L, Engels D, Montresor A, Savioli L: The global status of schistosomiasis and its control. Acta Trop. 2000, 77: 41-51. 10.1016/S0001-706X(00)00122-4.
Savioli L, Albonico M, Engels D, Montresor A: Progress in the prevention and control of schistosomiasis and soil-transmitted helminthiasis. Parasitol Int. 2004, 53: 103-113. 10.1016/j.parint.2004.01.001.
Haas W, Grabe K, Geis C, Pach T, Stoll K, Fuchs M, Haberl B, Loy C: Recognition and invasion of human skin by Schistosoma mansoni cercariae: the key-role of L-arginine. Parasitology. 2002, 124: 153-167. 10.1017/S0031182001001032.
Grabe K, Haas W: Navigation within host tissues: cercariae orientate towards dark after penetration. Parasitol Res. 2004, 93: 111-113. 10.1007/s00436-003-1045-9.
Shiff CJ, Cmelik SH, Ley HE, Kriel RL: The influence of human skin lipids on the cercarial penetration responses of Schistosoma haematobium and Schistosoma mansoni. J Parasitol. 1972, 58: 476-480. 10.2307/3278192.
Saladin KS: Schistosoma mansoni: cercarial responses to irradiance changes. J Parasitol. 1982, 68: 120-124. 10.2307/3281333.
Chai M, McManus DP, McInnes R, Moertel L, Tran M, Loukas A, Jonesa MK, Gobert GN: Transcriptome profiling of lung schistosomula, in vitro cultured schistosomula and adult Schistosoma japonicum. Cell Mol Life Sci. 2006, 63: 919-929. 10.1007/s00018-005-5578-1.
Sturrock RF: Schistosomiasis epidemiology and control: how did we get here and where should we go?. Mem Inst Oswaldo Cruz. 2001, 96 (Suppl): 17-27.
Wiest PM: The epidemiology of morbidity of schistosomiasis. Parasitol Today. 1996, 12: 215-220. 10.1016/0169-4758(96)10019-3.
Ismail MM, Farghaly AM, Dyab AK, Afify HA, el-Shafei MA: Resistance to praziquantel, effect of drug pressure and stability test. J Egypt Soc Parasitol. 2002, 32: 589-600.
Ismail M, Botros S, Metwally A, William S, Farghally A, Tao LF, Day TA, Bennett JL: Resistance to praziquantel: direct evidence from Schistosoma mansoni isolated from Egyptian villagers. Am J Trop Med Hyg. 1999, 60: 932-935.
Tran MH, Pearson MS, Bethony JM, Smyth DJ, Jones MK, Duke M, Don TA, McManus DP, Correa-Oliveira R, Loukas A: Tetraspanins on the surface of Schistosoma mansoni are protective antigens against schistosomiasis. Nat Med. 2006, 12: 835-840. 10.1038/nm1430.
Pearce EJ: Progress towards a vaccine for schistosomiasis. Acta Trop. 2003, 86: 309-313. 10.1016/S0001-706X(03)00062-7.
Simpson AJ, Sher A, McCutchan TF: The genome of Schistosoma mansoni: isolation of DNA, its size, bases and repetitive sequences. Mol Biochem Parasitol. 1982, 6: 125-137. 10.1016/0166-6851(82)90070-6.
Short RB, Menzel MY, Pathak S: Somatic chromosomes of Schistosoma mansoni. J Parasitol. 1979, 65: 471-473. 10.2307/3280303.
Grossman AI, McKenzie R, Cain GD: Sex heterochromatin in Schistosoma mansoni. J Parasitol. 1980, 66: 368-370. 10.2307/3280847.
Marx KA, Bizzaro JW, Blake RD, Hsien Tsai M, Feng Tao L: Experimental DNA melting behavior of the three major Schistosoma species. Mol Biochem Parasitol. 2000, 107: 303-307. 10.1016/S0166-6851(99)00235-2.
Johnson DH, Tsao J, Luo M, Carson M: SGCEdb: a flexible database and web interface integrating experimental results and analysis for structural genomics focusing on Caenorhabditis elegans. Nucleic Acids Res. 2006, 34: D471-474. 10.1093/nar/gkj036.
Verjovski-Almeida S, DeMarco R, Martins EA, Guimaraes PE, Ojopi EP, Paquola AC, Piazza JP, Nishiyama MY, Kitajima JP, Adamson RE, et al: Transcriptome analysis of the acoelomate human parasite Schistosoma mansoni. Nat Genet. 2003, 35: 148-157. 10.1038/ng1237.
Hoffmann KF, Johnston DA, Dunne DW: Identification of Schistosoma mansoni gender-associated gene transcripts by cDNA microarray profiling. Genome Biol. 2002, 3: RESEARCH0041-
Fitzpatrick JM, Johansen MV, Johnston DA, Dunne DW, Hoffmann KF: Gender-associated gene expression in two related strains of Schistosoma japonicum. Mol Biochem Parasitol. 2004, 136: 191-209. 10.1016/j.molbiopara.2004.03.014.
DeMarco R, Oliveira KC, Venancio TM, Verjovski-Almeida S: Gender biased differential alternative splicing patterns of the transcriptional cofactor CA150 gene in Schistosoma mansoni. Mol Biochem Parasitol. 2006, 150: 123-131. 10.1016/j.molbiopara.2006.07.002.
Fitzpatrick JM, Hoffmann KF: Dioecious Schistosoma mansoni express divergent gene repertoires regulated by pairing. Int J Parasitol. 2006, 36: 1081-1089. 10.1016/j.ijpara.2006.06.007.
Vermeire JJ, Taft AS, Hoffmann KF, Fitzpatrick JM, Yoshino TP: Schistosoma mansoni: DNA microarray gene expression profiling during the miracidium-to-mother sporocyst transformation. Mol Biochem Parasitol. 2006
Gobert GN, McInnes R, Moertel L, Nelson C, Jones MK, Hu W, McManus DP: Transcriptomics tool for the human Schistosoma blood flukes using microarray gene expression profiling. Exp Parasitol. 2006, 114: 160-172. 10.1016/j.exppara.2006.03.003.
Liu F, Lu J, Hu W, Wang SY, Cui SJ, Chi M, Yan Q, Wang XR, Song HD, Xu XN, et al: New perspectives on host-parasite interplay by comparative transcriptomic and proteomic analyses of Schistosoma japonicum. PLoS Pathog. 2006, 2: e29-10.1371/journal.ppat.0020029.
(TIGR), The Institute for Genomic Research. 2003, [http://www.tigr.org]
Knudsen GM, Medzihradszky KF, Lim KC, Hansell E, McKerrow JH: Proteomic analysis of Schistosoma mansoni cercarial secretions. Mol Cell Proteomics. 2005, 4: 1862-1875. 10.1074/mcp.M500097-MCP200.
Valle C, Festucci A, Calogero A, Macri P, Mecozzi B, Liberti P, Cioli D: Stage-specific expression of a Schistosoma mansoni polypeptide similar to the vertebrate regulatory protein stathmin. J Biol Chem. 1999, 274: 33869-33874. 10.1074/jbc.274.48.33869.
Bogitsh BJ, Dalton JP, Brady CP, Brindley PJ: Gut-associated immunolocalization of the Schistosoma mansoni cysteine proteases, SmCL1 and SmCL2. J Parasitol. 2001, 87: 237-241.
Hashmi S, Britton C, Liu J, Guiliano DB, Oksov Y, Lustigman S: Cathepsin L is essential for embryogenesis and development of Caenorhabditis elegans. J Biol Chem. 2002, 277: 3477-3486. 10.1074/jbc.M106117200.
Salter JP, Choe Y, Albrecht H, Franklin C, Lim KC, Craik CS, McKerrow JH: Cercarial elastase is encoded by a functionally conserved gene family across multiple species of schistosomes. J Biol Chem. 2002, 277: 24618-24624. 10.1074/jbc.M202364200.
Sajid M, McKerrow JH: Cysteine proteases of parasitic organisms. Mol Biochem Parasitol. 2002, 120: 1-21. 10.1016/S0166-6851(01)00438-8.
Caffrey CR, McKerrow JH, Salter JP, Sajid M: Blood 'n' guts: an update on schistosome digestive peptidases. Trends Parasitol. 2004, 20: 241-248. 10.1016/j.pt.2004.03.004.
Sajid M, McKerrow JH, Hansell E, Mathieu MA, Lucas KD, Hsieh I, Greenbaum D, Bogyo M, Salter JP, Lim KC, et al: Functional expression and characterization of Schistosoma mansoni cathepsin B and its trans-activation by an endogenous asparaginyl endopeptidase. Mol Biochem Parasitol. 2003, 131: 65-75. 10.1016/S0166-6851(03)00194-4.
Dvorak J, Delcroix M, Rossi A, Vopalensky V, Pospisek M, Sedinova M, Mikes L, Sajid M, Sali A, McKerrow JH, et al: Multiple cathepsin B isoforms in schistosomula of Trichobilharzia regenti: identification, characterisation and putative role in migration and nutrition. Int J Parasitol. 2005, 35: 895-910. 10.1016/j.ijpara.2005.02.018.
Fantappie MR, Freebern WJ, Osman A, LaDuca J, Niles EG, LoVerde PT: Evaluation of Schistosoma mansoni retinoid X receptor (SmRXR1 and SmRXR2) activity and tissue distribution. Mol Biochem Parasitol. 2001, 115: 87-99. 10.1016/S0166-6851(01)00274-2.
Osman A, Niles EG, LoVerde PT: Identification and characterization of a Smad2 homologue from Schistosoma mansoni, a transforming growth factor-beta signal transducer. J Biol Chem. 2001, 276: 10072-10082. 10.1074/jbc.M005933200.
Ganley-Leal LM, Guarner J, Todd CW, Da'Dara AA, Freeman GL, Boyer AE, Harn DA, Secor WE: Comparison of Schistosoma mansoni irradiated cercariae and Sm23 DNA vaccines. Parasite Immunol. 2005, 27: 341-349. 10.1111/j.1365-3024.2005.00785.x.
Simpson AJ: Tegumental proteins of Schistosoma mansoni: complex biomolecules and potent antigens. Mem Inst Oswaldo Cruz. 1992, 87 (Suppl 4): 11-17.
Wright MD, Henkle KJ, Mitchell GF: An immunogenic Mr 23,000 integral membrane protein of Schistosoma mansoni worms that closely resembles a human tumor-associated antigen. J Immunol. 1990, 144: 3195-3200.
Reynolds SR, Shoemaker CB, Harn DA: T and B cell epitope mapping of SM23, an integral membrane protein of Schistosoma mansoni. J Immunol. 1992, 149: 3995-4001.
Koster B, Hall MR, Strand M: Schistosoma mansoni: immunoreactivity of human sera with the surface antigen Sm23. Exp Parasitol. 1993, 77: 282-294. 10.1006/expr.1993.1086.
Inal J, Bickle Q: Sequence and immunogenicity of the 23-kDa transmembrane antigen of Schistosoma haematobium. Mol Biochem Parasitol. 1995, 74: 217-221. 10.1016/0166-6851(95)02490-5.
Koster B, Strand M: Schistosoma mansoni: Sm23 is a transmembrane protein that also contains a glycosylphosphatidylinositol anchor. Arch Biochem Biophys. 1994, 310: 108-117. 10.1006/abbi.1994.1146.
Petzke MM, Suri PK, Bungiro R, Goldberg M, Taylor SF, Ranji S, Taylor H, McCray JW, Knopf PM: Schistosoma mansoni gene GP22 encodes the tegumental antigen sm25: (1) antibodies to a predicted B-cell epitope of Sm25 cross-react with other candidate vaccine worm antigens; (2) characterization of a recombinant product containing tandem-repeats of this peptide as a vaccine. Parasite Immunol. 2000, 22: 381-395. 10.1046/j.1365-3024.2000.00316.x.
Almeida MS, Torloni H, Lee-Ho P, Vilar MM, Thaumaturgo N, Simpson AJ, Tendler M: Vaccination against Fasciola hepatica infection using a Schistosoma mansoni defined recombinant antigen, Sm14. Parasite Immunol. 2003, 25: 135-137. 10.1046/j.1365-3024.2003.00619.x.
Tendler M, Vilar MM, Brito CA, Freire NM, Katz N, Simpson A: Vaccination against schistosomiasis and fascioliasis with the new recombinant antigen Sm14: potential basis of a multi-valent anti-helminth vaccine?. Mem Inst Oswaldo Cruz. 1995, 90: 255-256.
Mitchell GF, Davern KM, Wood SM, Wright MD, Argyropoulos VP, McLeod KS, Tiu WU, Garcia EG: Attempts to induce resistance in mice to Schistosoma japonicum and Schistosoma mansoni by exposure to crude schistosome antigens plus cloned glutathione-S-transferases. Immunol Cell Biol. 1990, 68 (Pt 6): 377-385.
Davern KM, Tiu WU, Samaras N, Gearing DP, Hall BE, Garcia EG, Mitchell GF: Schistosoma japonicum: monoclonal antibodies to the Mr 26,000 schistosome glutathione S-transferase (Sj26) in an assay for circulating antigen in infected individuals. Exp Parasitol. 1990, 70: 293-304. 10.1016/0014-4894(90)90111-O.
Mitchell GF: Glutathione S-transferases - potential components of anti-schistosome vaccines?. Parasitol Today. 1989, 5: 34-37. 10.1016/0169-4758(89)90185-3.
Grzych JM, De Bont J, Liu J, Neyrinck JL, Fontaine J, Vercruysse J, Capron A: Relationship of impairment of schistosome 28-kilodalton glutathione S-transferase (GST) activity to expression of immunity to Schistosoma mattheei in calves vaccinated with recombinant Schistosoma bovis 28-kilodalton GST. Infect Immun. 1998, 66: 1142-1148.
Fitzpatrick JM, Johnston DA, Williams GW, Williams DJ, Freeman TC, Dunne DW, Hoffmann KF: An oligonucleotide microarray for transcriptome analysis of Schistosoma mansoni and its application/use to investigate gender-associated gene expression. Mol Biochem Parasitol. 2005, 141: 1-13. 10.1016/j.molbiopara.2005.01.007.
Kregel KC: Heat shock proteins: modifying factors in physiological stress responses and acquired thermotolerance. J Appl Physiol. 2002, 92: 2177-2186.
Borgnia M, Nielsen S, Engel A, Agre P: Cellular and molecular biology of the aquaporin water channels. Annu Rev Biochem. 1999, 68: 425-458. 10.1146/annurev.biochem.68.1.425.
Abrami L, Gobin R, Berthonaud V, Thanh HL, Chevalier J, Ripoche P, Verbavatz JM: Localization of the FA-CHIP water channel in frog urinary bladder. Eur J Cell Biol. 1997, 73: 215-221.
Elmendorf JS: Signals that regulate GLUT4 translocation. J Membr Biol. 2002, 190: 167-174. 10.1007/s00232-002-1035-3.
Durlinger AL, Visser JA, Themmen AP: Regulation of ovarian function: the role of anti-Mullerian hormone. Reproduction. 2002, 124: 601-609. 10.1530/rep.0.1240601.
Josso N, di Clemente N, Gouedard L: Anti-Mullerian hormone and its receptors. Mol Cell Endocrinol. 2001, 179: 25-32. 10.1016/S0303-7207(01)00467-1.
Dillon GP, Feltwell T, Skelton JP, Ashton PD, Coulson PS, Quail MA, Nikolaidou-Katsaridou N, Wilson RA, Ivens AC: Microarray analysis identifies genes preferentially expressed in the lung schistosomulum of Schistosoma mansoni. Int J Parasitol. 2006, 36: 1-8. 10.1016/j.ijpara.2005.10.008.
Davies SJ, Lim KC, Blank RB, Kim JH, Lucas KD, Hernandez DC, Sedgwick JD, McKerrow JH: Involvement of TNF in limiting liver pathology and promoting parasite survival during schistosome infection. Int J Parasitol. 2004, 34: 27-36. 10.1016/j.ijpara.2003.10.010.
Davies SJ, McKerrow JH: Developmental plasticity in schistosomes and other helminths. Int J Parasitol. 2003, 33: 1277-1284. 10.1016/S0020-7519(03)00161-9.
Imbeaud S, Faure E, Lamarre I, Mattei MG, di Clemente N, Tizard R, Carre-Eusebe D, Belville C, Tragethon L, Tonkin C, et al: Insensitivity to anti-mullerian hormone due to a mutation in the human anti-mullerian hormone receptor. Nat Genet. 1995, 11: 382-388. 10.1038/ng1295-382.
Lee JW, Choi HS, Gyuris J, Brent R, Moore DD: Two classes of proteins dependent on either the presence or absence of thyroid hormone for interaction with the thyroid hormone receptor. Mol Endocrinol. 1995, 9: 243-254. 10.1210/me.9.2.243.
Metzstein MM, Stanfield GM, Horvitz HR: Genetics of programmed cell death in C. elegans: past, present and future. Trends Genet. 1998, 14: 410-416. 10.1016/S0168-9525(98)01573-X.
Brindley PJ, Laha T, McManus DP, Loukas A: Mobile genetic elements colonizing the genomes of metazoan parasites. Trends Parasitol. 2003, 19: 79-87. 10.1016/S1471-4922(02)00061-2.
Mittleider D, Green LC, Mann VH, Michael SF, Didier ES, Brindley PJ: Sequence survey of the genome of the opportunistic microsporidian pathogen, Vittaforma corneae. J Eukaryot Microbiol. 2002, 49: 393-401. 10.1111/j.1550-7408.2002.tb00218.x.
Moertel L, McManus DP, Piva TJ, Young L, McInnes RL, Gobert GN: Oligonucleotide microarray analysis of strain- and gender-associated gene expression in the human blood fluke, Schistosoma japonicum. Mol Cell Probes. 2006
Holzer TR, McMaster WR, Forney JD: Expression profiling by whole-genome interspecies microarray hybridization reveals differential gene expression in procyclic promastigotes, lesion-derived amastigotes, and axenic amastigotes in Leishmania mexicana. Mol Biochem Parasitol. 2006, 146: 198-218. 10.1016/j.molbiopara.2005.12.009.
Curwen RS, Ashton PD, Johnston DA, Wilson RA: The Schistosoma mansoni soluble proteome: a comparison across four life-cycle stages. Mol Biochem Parasitol. 2004, 138: 57-66. 10.1016/j.molbiopara.2004.06.016.
Lim HK, Heyneman D: Intramolluscan inter-trematode antagonism: a review of factors influencing the host-parasite system and its possible role in biological control. Adv Parasitol. 1972, 10: 191-268.
McKerrow JH, Pino-Heiss S, Lindquist R, Werb Z: Purification and characterization of an elastinolytic proteinase secreted by cercariae of Schistosoma mansoni. J Biol Chem. 1985, 260: 3703-3707.
Colley DG, Wikel SK: Schistosoma mansoni: simplified method for the production of schistosomules. Exp Parasitol. 1974, 35: 44-51. 10.1016/0014-4894(74)90005-8.
Smithers SR, Roodyn DB, Wilson RJ: Biochemical and morphological characteristics of subcellular fractions of male Schistosoma Mansoni. Exp Parasitol. 1965, 16: 195-206. 10.1016/0014-4894(65)90044-5.
Llinas M, Bozdech Z, Wong ED, Adai AT, DeRisi JL: Comparative whole genome transcriptome analysis of three Plasmodium falciparum strains. Nucleic Acids Res. 2006, 34: 1166-1173. 10.1093/nar/gkj517.
Eisen MB, Spellman PT, Brown PO, Botstein D: Cluster analysis and display of genome-wide expression patterns. Proc Natl Acad Sci USA. 1998, 95: 14863-14868. 10.1073/pnas.95.25.14863.
Saldanha AJ: Java Treeview - extensible visualization of microarray data. Bioinformatics. 2004, 20: 3246-3248. 10.1093/bioinformatics/bth349.
S. mansoni Gene Index. 2006, [http://compbio.dfci.harvard.edu/tgi]
We thank Dr P'ng Loke, Dr Adam Carroll, Dr Kael Fisher, Dr Anjan Debnath and members of the McKerrow lab for technical support and helpful discussions. This research was supported in part by a UNCF/Merck Postdoctoral Fellowship (ERJ), the Sandler Family Supporting Foundation, NIAID ROI AI53247 and a US-Egypt Joint Science and Technologygrant forJunior Scientist Development (MMB).
Electronic supplementary material
Additional data file 1: To rule out dye specific effects, we performed dye-switch experiments for stage-specific comparisons. We plotted the log-ratio of the median (Cy5 signal versus Cy3 signal) from one chip verses the other dye-switched chip. Here we use a color scale to represent the density of each data point plotted. Most data points are accumulated at the red regions. If gene expression differences are small or subtle between two stages ((a) sporocyst versus cercariae and (b) adults versus cercariae), we expect their signals would center around the origin and be randomly distributed (red regions). For genes that are highly expressed, the signals would be high in one chip but low on the other chip. The scatter plots show the expected strong negative correlation between the dye-switched chips. The genes in the adult and sporocyst stages that are highly expressed are not artifacts due to the choice of dye. (EPS 1 MB)
Additional data file 3: TC names, annotation and relative ratio of sporocyst to cercariae and adult to cercariae are seen in columns A-F in duplicate. The background was removed from all samples manually as described in type II analysis after normalization in NOMAD, and analyzed using Gen Cluster 3.0 where samples were log transformed and screened to exclude samples with a deviation less than 0.6. All ratios are listed as log2 values. (XLS 496 KB)
Additional data file 4: The data presented are not ratios, but are averages of raw expression data of each stage independently, where 65,367 is the highest value (we consider values <500 as background) (XLS 1 MB)
Additional data file 5: All 431 genes highly enriched from the clustering analysis shown in Figure 2 (PDF 51 KB)
Additional data file 6: Transcript comparisons are between sporocysts (Sp) and cercariae (Ce), and adult worms (Ad) and cercariae (PDF 16 KB)
Additional data file 8: Following experiments done within this paper, the TIGR database containing TC sequences was updated and TC names were changed and expanded. Here is a key mapping TC sequences used in this paper to the new database . (XLS 4 MB)
Authors’ original submitted files for images
Below are the links to the authors’ original submitted files for images.
Rights and permissions
About this article
Cite this article
Jolly, E.R., Chin, CS., Miller, S. et al. Gene expression patterns during adaptation of a helminth parasite to different environmental niches. Genome Biol 8, R65 (2007). https://doi.org/10.1186/gb-2007-8-4-r65
- Additional Data File
- Adult Worm
- Mammalian Host
- Female Worm
- Intermediate Host Snail