The in vivo transcriptome of Schistosoma mansoni in the prominent vector species Biomphalaria pfeifferi with supporting observations from Biomphalaria glabrata

Background The full scope of the genes expressed by schistosomes during intramolluscan development has yet to be characterized. Understanding the gene products deployed by larval schistosomes in their snail hosts will provide insights into their establishment, maintenance, asexual reproduction, ability to castrate their hosts, and their prolific production of human-infective cercariae. Using the Illumina platform, the intramolluscan transcriptome of Schistosoma mansoni was investigated in field-derived specimens of the prominent vector species Biomphalaria pfeifferi at 1 and 3 days post infection (d) and from snails shedding cercariae. These S. mansoni samples were derived from the same snails used in our complementary B. pfeifferi transcriptomic study. We supplemented this view with microarray analyses of S. mansoni from B. glabrata at 2d, 4d, 8d, 16d, and 32d to highlight robust features of S. mansoni transcription, even when a different technique and vector species was used. Principal findings Transcripts representing at least 7,740 (66%) of known S. mansoni genes were expressed during intramolluscan development, with the greatest number expressed in snails shedding cercariae. Many transcripts were constitutively expressed throughout development featuring membrane transporters, and metabolic enzymes involved in protein and nucleic acid synthesis and cell division. Several proteases and protease inhibitors were expressed at all stages, including some proteases usually associated with cercariae. Transcripts associated with G-protein coupled receptors, germ cell perpetuation, and stress responses and defense were well represented. We noted transcripts homologous to planarian anti-bacterial factors, several neural development or neuropeptide transcripts including neuropeptide Y, and receptors that may be associated with schistosome germinal cell maintenance that could also impact host reproduction. In at least one snail the presence of larvae of another digenean species (an amphistome) was associated with repressed S. mansoni transcriptional activity. Conclusions/Significance This in vivo study, emphasizing field-derived snails and schistosomes, but supplemented with observations from a lab model, provides a distinct view from previous studies of development of cultured intramolluscan stages from lab-maintained organisms. We found many highly represented transcripts with suspected or unknown functions, with connection to intramolluscan development yet to be elucidated.


Introduction
The vast majority of the estimated 18,000 species of digenetic trematodes depend on a molluscan host, usually a gastropod, in which to undertake the complex developmental program characterized by extensive asexual reproduction and production of numerous cercariae [1,2]. The extent to which this large lineage of parasites has remained true to its dependence on molluscs, and the evident success achieved by digeneans-including by some species responsible for causing human disease-pose fundamental questions of interest for parasitologists, evolutionary biologists, ecologists, developmental biologists and comparative immunologists.
There is much about digenean-gastropod associations worthy of study: the host specificity often shown; the manner by which digeneans establish intimate infections without provoking destructive host responses; the ability of digeneans to affect and manipulate the energy and resource budgets of their hosts, including to achieve host castration; the intricate developmental program featuring multiple stages, asexual reproduction and the perpetuation of the germinal cell lineage; and finally the tendency for some infections to persist for long periods of time, implying protection of the snail-digenean unit that might involve contributions by the parasite to promote its perpetuation. These common and enduring relationships are also targeted and exploited by other organisms, including competing digenean species. One way forward to gain a deeper understanding of all these processes is to acquire a comprehensive overview of the genes expressed by host snails and larval digeneans during the course of infection. This in turn sets the stage for eventually learning more about how these two sets of gene products interact (the interactome) to influence the outcome of this interaction.
Because S. mansoni causes intestinal schistosomiasis in an estimated 166 million people in the Neotropics, Africa and Southwest Asia, it has long been intensively studied, in part because it is relatively easily maintained in a laboratory setting [3]. Many molecular tools, including a genome are available for S. mansoni [4]. Additionally, B. glabrata, the most important host for S. mansoni in the Neotropics, has become a model gastropod host, including with a recently available genome [5]. In Africa, several Biomphalaria species play an important role in transmission, with the most important being B. pfeifferi. The latter species occurs widely across sub-Saharan Africa, where >90% of the world's cases of schistosomiasis now occur. B. pfeifferi is probably responsible for transmitting more S. mansoni infections than any other snail species [6].
With respect to the intramolluscan development of S. mansoni, following penetration of miracidia, there is a 24 hour period of transformation into a mother sporocyst that can be thought of as one of parasite transition and establishment. The miracidium has a collection of germinal cells that are destined to give rise to the daughter sporocysts. Mitotic division of germinal cells begins as early as 24 hours and germinal cells proliferate notably in an enlarging mother sporocyst [7]. By 6 days after infection, all mother sporocysts have germinal balls (daughter sporocyst embryos) which occupy nearly the entire body cavity [8]. The embryonic daughter sporocysts, of which there are an average of 23 produced per mother sporocyst, continue to grow and elongate [9]. Daughter sporocysts exit at 12-14 days and start their migration to the digestive gland and ovotestis region of the snail [8,10]. Upon release of daughter sporocysts, mother sporocysts collapse and typically do not continue to produce daughter sporocysts, but they nonetheless persist in the head-foot of the snail. By 15-20 days, daughter sporocysts undergo a remarkable transformation, becoming amorphous, wedged between lobules of the digestive gland. Daughter sporocysts occupy a significant proportion of the snail host's body; 65% of the snail's digestive gland can be occupied by daughter sporocysts in a patent, cercariae-producing infection [11]. Within them, cercarial embryos undergo development, culminating in the production of a muscular tail and a body dominated by gland cells that are filled with lytic enzymes [12]. Once again, a separate allotment of germinal cells is sequestered in the cercarial body and these are destined to become the gonads and reproductive cell lineages of adult worms. Around 32 days post-exposure, cercariae exit from daughter sporocysts, migrate through the snail's body and emerge into the water, usually through hemorrhages in the mantle.
The timing of these events is dependent on temperature but cercarial shedding can occur as early as 19 days post-miracidial penetration [13]. In some situations, daughter sporocysts will produce granddaughter sporocysts in lieu of cercariae [13,14]. A typical consequence of infection is that snails are partially or totally castrated, the extent depending on whether they were infected before or after achieving maturity [10,15]. This interaction is remarkable in that some Biomphalaria snails can survive for over a year shedding cercariae daily [16], although there is considerable variability in the duration of survival of infected snails. The productivity of infections within snails has no doubt contributed greatly to the success of all digenetic trematodes, and in the case of human-infecting schistosomes, is a major factor complicating their control.
Despite the significant immunobiological, physiological, and reproductive changes inflicted upon infected Biomphalaria snails [17], we still lack a comprehensive picture of what the parasite is producing to effect such changes. ESTs (expressed sequence tags) from early-developing S. mansoni cercarial germ balls show their transcriptional responses are distinctive from those noted for miracidia, cercariae, schistosomulae, eggs, and adults [18]. Gene products with antioxidant activity, oxidoreductases, and intermolecular binding activity were represented in mother sporocyst-specific genes from a S. mansoni microarray study [19]. Proteomic analyses revealed the significance of proteases, protease inhibitors, heat shock proteins, redox/antioxidant enzymes, ion-binding proteins, and venom allergen-like (SmVAL) proteins released by in vitro cultured miracidia transforming into mother sporocysts [20,21]. Highlights of 6d and 20d cultured mother sporocysts transcript expression, as detected by SAGE, are an up-regulation of HSP 70, HSP 40, egg protein, and trypsinogen 1-like all exclusive to miracidia and sporocyst stages [22]. Differential expression (DE) of transcripts from LongSAGE of in vitro cultured miracidia, 6d and 20d mother sporocysts is dependent on whether or not the sporocysts are conditioned in medium with Bge cell products [22,23], reinforcing the need for S. mansoni in vivo transcript analysis.
Here, our primary focus is on presentation of RNA-Seq results for S. mansoni from the same field-derived Kenyan snails that comprised the B. pfeifferi transcriptomic study of Buddenborg et al. [17]. Field-derived snails were exposed experimentally to S. mansoni miracidia hatched from eggs from fecal samples from local schoolchildren and harvested 1 or 3 days later. Additionally, field snails found to be naturally shedding S. mansoni cercariae were chosen for study. Our goal was to provide in vivo views of establishment of early mother sporocyst development and shedding stages for snails and parasites taken directly from natural transmission sites. We did not investigate longer exposure intervals following experimental exposures because we did not want these snails to lose their field characteristics. We supplemented these observations with previously obtained results from a set of independent microarray experiments of lab-derived S. mansoni in B. glabrata at 2, 4, 8, 16, and 32d. Along with comparable parasite early establishment (2d, 4d) and shedding (32d) time points, additional time points from the array cover supplementary stages in development including production, release and migration of daughter sporocysts. Our approach is distinctive in its focus on in vivo life cycle stages, the inclusion of both snails naturally infected from an endemic area in western Kenya and of laboratory-maintained snails, and the use of two transcriptome technologies particularly as a means to identify robust features of S. mansoni intramolluscan transcription regardless of technique or vector species. We examined specific groups of transcripts to gain distinctive insights on intramolluscan development. The database we provide should also provide helpful information in eventually achieving a deeper understanding of the interactome that is the essence of this dynamic host-parasite interaction. each from pooled fecal samples (5 individuals) for 6 hours in 24-well culture plates and then returned to aquaria. At 1d and 3d post infection, snails were collected and stored in TRIzol until processing. Field-collected, cercariae-producing snails 6-9mm in diameter were used for the shedding sample group. Biological triplicates were sequenced for each sample group using Illumina HiSeq 2000 (Illumina, Carlsbad CA) at the National Center for Genome Resources (NCGR) in Santa Fe, NM. In addition, one naturally shedding B. pfeifferi snail was sequenced on a 454 sequencer (Roche, Basel Switzerland) to improve S. mansoni transcript assembly but these sequences were not used for quantification. See Buddenborg et al. [17] for a thorough explanation of RNA extraction, library preparation, sequencing procedures, and sequencing summaries.

S. mansoni microarray experiments
The M-line strain of B. glabrata infected with S. mansoni PR1 strain was used in the microarray experiments to monitor parasite transcriptional changes that occur during infection. Both snail and trematode were maintained at UNM as previously described [25]. Snails were exposed to 10 miracidia each of S. mansoni for 2d, 4d, 8d, 16d, or 32d (shedding snails), with biological triplicate replicates for each time point. An uninfected B. glabrata group was also used to account for cross-hybridization from mixed snail-trematode samples. Total RNA was extracted as previously described [26] and treated with DNAse I (Ambion UK) to remove gDNA contamination. RNA was quantified on a NanoDrop ND-1000 spectrophotometer and quality-assessed using an Agilent 2100 bioanalyzer. cDNA synthesis, amplification, labeling, and hybridization were performed as previously described [26].
A publicly available S. mansoni microarray (NCBI GEO accession GPL6936) representing 19,244 unique S. mansoni contigs (38,460 total experimental probes) was used with the following modification: all array probes were duplicated to allow for an added level of replicability. The transcript probes contained on the array were designed to profile 15 different developmental stages. Thus, many of the molecules likely important to larval development are present as well. Microarray images were recovered from a GenePix 4100A (Axon Instrument Inc.) dual channel laser scanner.
Raw data was averaged from replicates in each experimental group (2d, 4d, 8d, 16d, 32d), and for replicates in the uninfected snail group (Bg-only). For each experimental group, the mean and standard deviation were calculated, and values falling below one standard deviation from the mean were removed from further analysis. Features that were non-reactive for any of the groups used in this study, amounting to 26,581 probes, were removed as well as those that were cross-reactive with the Bg-only group (787 probes). The average Bg-only value was subtracted from experimental groups for each probe. Calculated expression values less than 1 were removed from the analysis and the remaining values were transformed by log base 2.
An updated annotation of array features was performed using BLASTn with the NCBI nucleotide database (sequence identity >70%, E-value <10 −06 ), BLASTp with the NCBI nonredundant protein database (sequence identity >40%, E-value <10 −06 ). Array features were matched to their homologous S. mansoni transcript by BLASTn against the assembled S. mansoni transcripts. These homologous transcripts were used for analyses comparing across Illumina and array samples.

S. mansoni transcriptome assembly and annotation
An overview of our analysis pipeline is shown in S1 Fig. After pre-processing all Illumina reads, those from Illumina and 454 sequencing that did not map back to the B. pfeifferi transcriptome or identified symbionts were assembled into contigs (assembled, overlapping reads). The separation of host, parasite, and symbiont reads is described in detail in Buddenborg et al. [17]. We employed Trinity v2.2 RNA-Seq de novo assembler [27,28] for de novo and genome-guided transcriptome assembling using paired-end reads only. The S. mansoni de novo assembly consisted of reads that did not map to the B. pfeifferi transcriptome or symbionts after alignment with Bowtie2 v2.2.9 [29]. The genome-guided transcriptome assembly was performed using STAR v.2.5 2-pass alignment [29,30] to the S. mansoni genome (GeneDB: S. mansoni v5.2).
Schistosoma mansoni genome-guided and de novo Trinity assemblies were concatenated and redundancy reduced using CD-Hit-EST at 95% similarity [31]. The resulting sequences were screened against B. glabrata (VectorBase: BglaB1) and S. mansoni genomes, peptides, and mRNAs using BLASTx and BLASTn (sequence identity >70%, E-value < 10 −12 ). Sequences with blast results mapping to B. glabrata were removed and remaining S. mansoni-specific sequences were thereafter referred to as transcripts.
All assembled transcripts were annotated based on their closest homologs and predicted functional domains in the following databases and tools: BLASTp with NCBI non-redundant protein database (sequence identity >40%, E-value <10 −06 ), BLASTn with NCBI nucleotide database (sequence identity >70%, E-value < 10 −06 ),

Illumina-derived S. mansoni transcriptomic characteristics
Throughout this discussion, a "transcript" is defined as assembled S. mansoni contigs formed from overlapping reads with the understanding that this includes both full-length transcripts, partial transcripts, and isoforms. For our Illumina-based study, a total of 23,602 transcripts made up our combined genome-guided and de novo assembled S. mansoni intramolluscan transcriptome. Microarray and Illumina expression data can be found in S1 File. Schistosoma mansoni assembly metrics are provided in S1 Table. When all raw reads from each infected snail were mapped to the S. mansoni transcripts, 1d, 3d, and shedding replicates' mapping percentages ranged from 4.01-5.46, 1.48-4.05, and 4.36-9.72, respectively (S2 Fig). The principal component analysis (PCA) plot (S3 Fig) shows that the percentage of S. mansoni reads varies, and that 1d and 3d groups show more variation between replicates than do shedding replicates. It is not surprising that the transcriptional responses among early replicates at 1 and 3d are more variable in this natural system involving both genetically variable snails and schistosomes, especially as compared to shedding snails which have reached a steady state of continued cercarial production. Also, because the parasite stages at 1 and 3d are small relative to their hosts, uniform sampling of their contributions may be harder to achieve.
Overall Illumina and microarray S. mansoni transcript expression. Transcripts with �1 Log 2 TPM in at least one replicate per group in Illumina samples and features with fluorescence � 1 in microarrays were considered for expression analyses. Based on these cutoffs, over fifteen thousand different transcripts representing 7,252 S. mansoni genes were detected in 1d S. mansoni infections. Following a dip in 3d samples, even more S. mansoni transcripts were expressed in shedding snails ( Fig 1A). The decline noted in the 3d Illumina samples may reflect that at least one replicate returned fewer S. mansoni reads in general or may simply reflect a sampling consideration due to the small size of the parasites relative to the snail at this time point. Sustained expression of a large number of transcripts with a general trend towards higher expression in shedding snails was also noted in the microarray data ( Fig  1B). Particularly for the Illumina results, some of the transcripts enumerated represent different portions of the same original full-length mRNAs as well as different isoforms, so the actual number of expressed genes is approximately half as many as the number of recorded transcripts. Nonetheless, the variety produced is impressive and generally supported by our microarray results as well (at least 6,000 features expressed at all time points).
Our datasets generated by Illumina and microarray analyses might be expected to return different results for at least four different reasons: 1) the two methods are totally different in approach; 2) the host snail species and S. mansoni strains differed; 3) the time points sampled differed; and 4) the transcripts represented on the array are more limited than whole transcriptomic sequencing provided by Illumina. However, they also provide independent views of the same basic process and enable identification of shared, strongly expressed features, so comparisons are warranted, especially so for shedding snails when the same developmental stage could be compared between techniques. Particularly noteworthy in both Illumina and array samples was that a large number of transcripts was shared across all time points (Fig 2). Among Illumina groups, >15% of all transcripts are expressed constitutively and among all microarray groups, >34% of all probes were expressed constitutively. This is suggestive of a core transcriptome required of schistosomes living in snails (see below for more details as to what comprises this core transcriptome). When comparing both early time points (1-4d) and shedding time points from both Illumina and array methods, venn diagrams not surprisingly indicate that Illumina RNA-Seq detects more S. mansoni transcripts.

The intra-molluscan metabolic landscape
After successful penetration of the snail host, digeneans alter their metabolism to depend completely on the resources available in the molluscan host and shift their energy budget towards sporocyst and/or rediae development. One of the unique evolutionary innovations of the Neodermata is the syncytial tegument, a vital aspect of digenean biology providing both protection and a highly efficient route to acquire nutrients from the host species [39]. Through the tegument, schistosomes acquire most nutrients and other key molecules via facilitated or active transport using transmembrane transporters [40]. Glucose transporters are expressed in both adult and larval stages of S. mansoni [41,42]. While miracidia in water employ aerobic energy metabolism, after 24 hours of in vitro cultivation, sporocysts shift their metabolism towards lactate production [43]. Expression of glucose transporters is particularly important in initial establishment (1d) and in shedding snails (Fig 3). By the 3d day post infection, the parasite up-regulates metabolic processes that are part of the purine salvage pathway and nucleotide biosynthesis, highlighting its transition to reproduction processes and mitosis. This is concurrent with a down-regulation of phosphorylation and general mitochondrial metabolic activities. This highlights the transition to the less aerobic regime within the host, and the shift to a tightly regulated reproductive program rather than active migration within the host or the environment. It has also been observed that daughter sporocysts have fewer mitochondria [44,45]. This shift to anaerobic mode of energy production is reversed by the presence of fullyformed cercariae developing in the sporocysts of actively shedding snails and is corroborated by the fact that aerobic respiration is especially active in the tails of cercariae [46]. This is expected due to the fact that cercariae, once released from the snail, have an active lifestyle and must generate enough energy from a limited amount of stored glycogen. Oxidative phosphorylation, aided by the greater availability of oxygen in the aquatic environment, helps cercariae fulfill their demanding energy requirements.
While obtaining organic carbon from the host fulfills the energetic requirements of the parasite, any actual growth is nitrogen dependent. Acquiring amino acids and other important building-block molecules is thus paramount to the parasite's fitness. Tegumental amino acid transporters have not been previously reported in S. mansoni sporocysts [40] but here we provide evidence of the expression of several amino acid transporters across all intramolluscan time points, some of which may be tegumental. A glutamate transporter was the highest expressed amino acid transporter across all replicates in the Illumina samples. This is concurrent with an increase of amino acid biosynthesis by 3d which continues, but to a lesser degree in shedding snails (e.g. Alanine transaminase EC: 2.6.1.2, Glutamate Synthase EC 1.4.1.13). Nucleoside transporters were also abundantly expressed, especially in shedding snails, as noted both by Illumina and microarray results.
Components of receptor-mediated endocytosis are present in the transcriptome of free-living and adult stages of S. mansoni [18]. Transcripts necessary for clathrin-mediated endocytosis, including clathrin assembly proteins, low-density lipoproteins, and adapter complex Ap2 were present in the intramolluscan transcriptome. The regulator of endocytosis, dynamin, was also present. These transcripts may be used in endocytosis to bring in lipids needed to make membranes. Expression of transcripts involved in receptor-mediated endocytosis, and possibly also in exocytosis, was high immediately upon miracidial transformation in 1d S. mansoni mother sporocysts.
We identified additional putative transmembrane transporters using the Transporter Classification Database [47]. The most abundantly expressed transporters at 1d were AAA-ATPase and protein kinase superfamilies whereas the nuclear pore complex, H+ or translocating NADH dehydrogenase, and endoplasmic reticular retrotranslocon families were dominant in 3d and shedding groups. ABC transporters were present in all Illumina samples, with 27 ABC transporter transcripts expressed in shedding snails. Three transcripts, identified as an ATPbinding cassette sub-family F member 2-like isoform X1, isoform X2, and ATP-binding cassette sub-family E member 1-like were highly expressed in every Illumina replicate. It has been suggested that ABC transporters serve an excretory function for adult schistosomes, playing a role in the removal of xenobiotics and/or influencing interactions with the definitive host [48]. The high expression of ABC transporters in intramolluscan stages, particularly in shedding snails, suggests they have an important and as yet not fully appreciated role in development. Perhaps they play a role in elimination of wastes associated with production of cercariae or facilitate release of factors that modify the immediate environment of the daughter sporocysts to favor their continued productivity of cercariae.
Protein kinases phosphorylate intracellular proteins in order to alter gene expression and are responsible for many basic cellular functions. In schistosomes, kinases are predicted to play a role in host invasion, sensory behavior, growth, and development [49]. Because of their importance, kinases have been used as potential pharmaceutical targets against S. mansoni [50]. A BLASTx homology search of kinases from Kinase SARfari (https://www.ebi.ac.uk/ chembl/sarfari/kinasesarfari) confirmed the representation of 19 kinases from 4 different superfamilies on the microarray, and 154 kinases belonging to 7 different superfamilies expressed in Illumina 1d, 3d, and shedding samples (S5 Fig). The highest expressed protein kinases are members of the group CMGC which includes MAPK growth and stress response kinases, cell cycle cyclin dependent kinases, and kinases for splicing and metabolic control.

Protease and protease inhibitor transcripts expressed at different stages of parasite development
The protease-encoding genes of parasitic helminths have undergone gene duplication and divergence, and by enabling helminths to process diverse proteinaceous substrates are believed to be critical to establishment and perpetuation of infection [51,52]. Helminth proteases and protease inhibitors have proven useful as markers for diagnostics purposes, or as targets for drugs or vaccines [52][53][54]]. In the snail host, larval schistosomes use proteases for nutrient acquisition, to create the space needed for their expansive growth, and for defense functions, potentially destroying or inhibiting lytic host proteases [55]. Miracidia release proteases to facilitate entry into the snail host, often into dense tissue of the head-foot [56]. In vitro studies of cultured mother sporocysts have revealed secretion of proteases facilitating degradation of snail hemolymph proteins such as hemoglobin [55].
We observed that intramolluscan S. mansoni devotes considerable effort to making proteases and protease inhibitors with 397 protease transcripts and 77 protease inhibitor transcripts represented in at least one time point (S6 Fig). Replicates of each Illumina time point with the lowest percentage of S. mansoni reads (1d-R1, 3d-R2, shedding-R3) also had the least abundant number of transcripts identified as proteases and protease inhibitors. One-day infections (see 1d-R2 and 1d-R3) with higher read counts indicative of robust development show expression of a gamut of S. mansoni proteases that somewhat surprisingly resemble those produced by S. mansoni in shedding snails. Coincidentally, we noted the snail host up-regulated expression of protease inhibitors especially during larval establishment at 1d and 3d [17].
At all time points more S. mansoni proteases were present than protease inhibitors and, in general, protease inhibitors and proteases increased in abundance and expression as infection progressed. For both Illumina and microarray samples, shedding snails had both the greatest number of proteases and protease inhibitors expressed relative to other time points, and the highest expression levels of proteases and protease inhibitors.
As expected, elastases, an expanded family of serine proteases in S. mansoni, were the most highly expressed proteases in S. mansoni from both B. glabrata and B. pfeifferi shedding snails (Fig 4). We identified 9 elastase transcripts including those previously designated as cercarial elastases 1a and 2b and found in daughter sporocysts and cercariae [57]. Although elastases are known to be used in definitive host skin penetration, active translation of SmCE2b into protein sequences is seen prior to exiting the snail and was postulated to be involved in facilitating egress from the snail [57]. The intramolluscan transcriptome of Schistosoma mansoni The intramolluscan transcriptome of Schistosoma mansoni Our data not only corroborate the presence of SmCE2b in shedding snails, but also reveal this and other S. mansoni elastases to be present at all time points we examined, in contrast to microarray results previously reported with early stage sporocysts from in vitro cultures [23]. For example, even at 1d (see 1d-R2) we found high expression of six S. mansoni elastases, some of which are those noted prominently in cercariae [57]. Our microarray samples also show expression of several elastases at all time points. It is not unusual to think that early-stage larval S. mansoni would express protease activity as they too must implement host penetration. Wu et al.
[20] noted a conspicuous absence of elastase proteins in in vitro larval transformation products but other proteases present suggested an obvious degree of overlap between cercarial versus larval protease repertoires.
Leishmanolysin (also called invadolysin), a metalloprotease, is the second most abundant type of secreted protease of cercariae after elastases [58]. Functional studies of leishmanolyin in larval S. mansoni suggested this protease is capable of interfering with the migration of B. glabrata hemocytes and may influence the establishment of infection [59]. Leishmanolysin has also been detected among the proteins accompanying transformation of miracidia to mother sporocysts [20]. We detected leishmanolysin transcripts at all time points, and they were most abundant in shedding snails, likely indicative of their representation in developing cercariae ( Fig 4B).
Cathepsins are papain-like cysteine proteases and have been identified in the S. mansoni miracidia proteome, transforming miracidia, and mother sporocysts [55,60] and are implicated in tissue penetration, digestion and immune evasion in the definitive host [52,[61][62][63][64]. Cathepsins take the place of tissue-invasive elastases in the cercariae of avian schistosomes [65]. Of two cathepsin B transcripts we noted, we found one expressed in all replicates except from 1d-R1 and 3d-R2, the early-stages replicates noted to have lower S. mansoni read counts (Fig 4C). Schistosoma mansoni expresses cathepsin B in the flame cells of cercariae where they are believed to play a role in osmoregulation and/or secretion [66]. Cathepsin C, involved in acquisition of oligopeptides and free amino acids by larval schistosomes [67], was also identified by Illumina at 1d, 3d, and in shedding snails, with the exception of replicate 3d-R2 which had a pre-patent amphistome infection. Cathepsins L1 and L3 were highly expressed by mother sporocyst stages (2d, 3d, 4d, 8d samples) in the microarray samples. At 16d, when daughter sporocysts are migrating through host tissue and hemolymph to the digestive gland, the proteases produced most closely resemble those from 32d shedding infections, including cathepsin C.
In contrast to proteases, there is relatively little information about protease inhibitors and their roles in parasite development and survival (see [53] for a thorough review of schistosome protease inhibitors). One of the better-characterized groups is the serine protease inhibitors (serpins; MEROPS clan ID, family I4) that may play a role in both post-translational regulation of schistosome proteases and defense against host proteases [68]. Serpins were expressed in all the time points sampled but we observed the highest expression of serpins at 1d and in shedding snails. The most abundant protease inhibitors in the Illumina samples (1d, 3d, shedding) were those that belong to the JF clan which is interesting because it is by no means the most abundantly represented clan, comprised of only one family called cytotoxic T-lymphocyte antigen-2-alpha (CTLA-2a), known to induce apoptosis of T-lymphoma cells in schistosomeinfected mice [69]. This gene homolog is not represented on the S. mansoni microarray which accounts for its absence in those samples. The homologous CTLA-2a transcripts expressed in the intramolluscan stages of S. mansoni may play a similar role in apoptosis or immunomodulation in snails to facilitate maintenance of long-term infections.
Transcripts identified as the protease inhibitor aprotinin (IB clan), a trypsin inhibitor, were moderately expressed in Illumina 1d-R2 and R3 replicates and in all replicates of shedding snails. In the plasma of Biomphalaria, the phenoloxidase enzyme laccase, whose activity is enhanced by trypsin, induces a negative impact on late (7-9 week) S. mansoni infections [70]. We noted an up-regulation of snail-produced trypsins in B. pfeifferi shedding S. mansoni cercariae [17] as compared to uninfected controls. By inhibiting snail-produced trypsins, S. mansoni daughter sporocysts and/or developing cercariae within may disable an important snail defense strategy.

The S. mansoni venom allergen-like proteins (SmVALs)
The venom allergen-like proteins (SmVAL2, 3/23, 9, 15, 26/28, and 27) have been identified as secreted larval transformation proteins [20]. SmVAL proteins can be found throughout miracidia and sporocyst parenchymal cell vesicles and in germinal cells with evidence for involvement in larval tissue remodeling and development by regulating snail matrix metalloproteinases [20,71]. One and 3d Illumina samples showed variable expression of SmVALs 1, 11, 14, and 22 (Fig 5). Replicates from shedding snails had more consistent SmVAL profiles, with 14 different SmVAL homologs found among Illumina replicates and 9 SmVAL homologs in the 32d microarray samples. SmVAL1 was ubiquitously expressed across 1d, 3d, and shedding Illumina samples. Chalmers et al. [72] also noted abundant SmVAL transcripts in the infective stages of S. mansoni, namely miracidia and cercariae. SmVALs 4 and 24 transcripts, localized to the preacetabular glands of developing cercariae [73] were also the highest expressed SmVAL transcripts we found in shedding S. mansoni. SmVAL16 was localized close to the neural ganglia of adult male worms [73]; we detected its expression at 1d, 3d, and shedding time points. The repertoire of SmVAL proteins secreted during transformation may differ from the SmVAL transcripts being produced and this may account for the differences in the SmVAL transcripts we report here versus previously published proteomic findings.

S. mansoni intramolluscan G-protein coupled receptors (GPCRs)
G-protein coupled receptors or GPCRs are the largest superfamily of transmembrane proteins in eukaryotes responsible for facilitating signaling affecting various downstream functions like development, reproduction, neuronal control of musculature and more [74,75]. As receptors, GPCRs are involved in mediating a variety of processes critical to schistosome survival including mediating host-parasite interactions, reproduction, and mating [76]. Praziquantel has been identified as a GPCR ligand acting to modulate serotoninergic signaling [76,77]. Several in silico studies identifying and characterizing the S. mansoni "GPCRome" [78,79] culminated in the classification of a broad range of phylogenetically distinct clades/classes of GPCRs [80]. S. mansoni microarray studies have reported diverse expression patterns of individual GPCRs, with the overall highest expression occurring in 3-7 week worms, indicating that they are associated with complex stage-specific roles [23,80].
In intramolluscan stages, we identified 78 GPCR transcripts from our Illumina samples, and 26 probes from microarray samples (Fig 6A). Many (38%) of the Illumina GPCR transcripts were A FLPR-like, a GPCR class containing receptors similar to FMRFamide GPCRs that invoke muscle fiber contractions in schistosomes by increasing calcium transport across voltage-gated calcium channels [81]. Shedding snails had the most diverse representation of GPCRs. One transcript, homologous to an identified S. mansoni GPCR (Smp_193810) with unknown function, was expressed at all time points with markedly high expression in both Illumina and microarray samples from shedding snails. A GPCR sensing the biogenic amine 5HT (Smp_126730) and known to be distributed throughout the adult worm's nervous system [74], was expressed at 1d and shedding Illumina samples (no homologous probe was found on the microarray). Its presence in intramolluscan stages suggests that serotonin-stimulated movement is essential throughout the life cycle of schistosomes. At 1d and shedding, a type 1 serotonin receptor is down-regulated in B. pfeifferi and at 3d, kynurenine 3-monooxygenase (important for its ability to degrade tryptophan and limit concentrations of serotonin) is upregulated [17]. Serotonin is a molecule of relevance to both the snail and parasite, and interference with its levels may be relevant to castration of snails (see concluding comments).

Neuropeptides and neural development
Studies on neuropeptides (peptide hormones) in planarian flatworms and their homologs in S. mansoni have identified their influence in locomotion, feeding, host location, regeneration, and development [82,83]. Lu et al. [84] reported the expression of putative neuropeptides and transcripts suspected to be involved in neural development from paired and unpaired female and male adult worms. Seventeen transcripts were identified as neuropeptide receptors from the Illumina transcriptome, all of which were GPCRs (Fig 6B).
Neuropeptides and their receptors were mostly absent at 1d and 3d but abundant in shedding-R1 and R2. In shedding-R3, the replicate with a muted S. mansoni response, only one neuropeptide receptor (neuropeptide Y receptor) was expressed. Shedding-R3 was curious in that protein 7b2 and NPP-1 (GFVRIamide) prepropeptide were highly expressed. Only one putative neuropeptide precursor (NPP-1 prepropeptide) was identified on the microarray with only~1 Log 2 fluorescence at 32d and was not present in any other sample. In adult worms, GFVRIamide is localized to neurons that run along the cerebral commissure towards the oral sucker [83]. Allatostatin receptor, a GPCR with ovary-specific transcription in adult S. mansoni [80], is important for reproductive development in Schistosoma japonicum adult females [83,85]. Four transcripts homologous to allatostatin receptor were expressed primarily in shedding S. mansoni replicates. Our results indicate expanded roles for neuropeptides and neural development transcripts previously uncharacterized in intramolluscan stages of S. mansoni.
We identified 33 of the 39 genes found to be involved in neural development by Lu et al. [84] in our Illumina S. mansoni transcriptome (Fig 7). Cell polarity proteins were the highest expressed transcripts involved in neural development at 1d, 3d, 16d, and shedding snails. 2d array S. mansoni showed little activity of transcripts related to neural development. In 4d and 8d samples, notch and septate junction transcripts were the most highly expressed neural development transcripts. Notch transcripts are highly expressed in eggs but not in cercariae and are thought to be mainly involved in S. mansoni oogenesis and embryogenesis within the vertebrate host but have been implicated in neurogenesis [18]. Lu et al. [84] found SOX to be The intramolluscan transcriptome of Schistosoma mansoni transcribed in the ovary of paired and unpaired females and its expression in germ balls has also been established [18,86]. Three transcripts homologous to the S. mansoni SOX transcription factor were present predominantly in 1d and shedding time points reinforcing the role of SOX transcription in embryonic and germinal cell development.

Transcripts associated with germinal cells and asexual reproduction of schistosomes in snails
A prominent feature of the complex developmental program of sporocysts in snails is the presence of germinal cells that give rise to embryos that come to contain both the somatic cells that eventually divide to comprise the bodies of either sporocysts or cercariae and more germinal cells. These germinal cells are then poised to give rise to the next generation. None of this asexual polyembryonic process involves the formation of gametes or evidence of fertilization. Germinal cells in S. mansoni sporocysts have been shown to share common molecular features with planarian neoblasts or stem cells, prompting the suggestion that the digenetic nature of the life cycle of schistosomes and other digenetic trematodes may have evolved because of the adaptation of a system of preservation of these stem cell-like germinal cells [7,18,86].
Consistent with Wang et al. [7] we observed expression of fibroblast growth factor receptors (fgfr), vasa, argonaute2 (ago2), and nanos transcripts shown to be associated with long-term maintenance of neoblast stem cells (Fig 8A). Expression of fgfr2, argonaut-2 and especially vasa are expressed in all samples, suggestive of their importance in intramolluscan development. The microarray had an additional fgfr feature (fgfr4) that was not detected in the Illumina transcriptome. Our results are in agreement with Wang et al. [7] that nanos-1 is not expressed in sporocysts, consistent with their suggestion that nanos-1 expression is exclusive in adult S. mansoni [7]. Nanos-2 expression was observed in every replicate of every time point with the exception of the 2d microarray sample. It has been proposed that there are two populations of germinal cells, nanos + and nanos -, with the latter population proliferating much more rapidly [7]. Vasa is needed for proliferation of both nanos + and nanosstem cell populations and ago2 is required for proliferation of only nanoscells. It is hypothesized that the two populations exist for different purposes: one a more undifferentiated stem cell-like population and the other a more differentiated one ready to enter embryogenesis [87]. The intramolluscan transcriptome of Schistosoma mansoni The proliferation of sporocysts and then cercariae by digenetic trematodes in snails is now generally considered to be an asexual process, one that does not involve gamete formation or fertilization [88], and it is frequently assumed that the progeny produced from a single miracidium are genetically the same. However, there are also persistent claims that the process is best considered as apomictic parthenogenesis [89]. Some observations indicate that the S. mansoni cercariae arising from a single miracidium are not genetically identical but exhibit some variation with respect to representation of repetitive elements that has been attributed to mitotic recombination [90,91]. Khalil and Cable [92] examined germinal development in rediae of Philopthalmus megalurus and concluded the process was diploid parthenogenesis. They observed the presence of cells interpreted to be oögonia entering meiotic prophase I up to the stage of diakinesis that was then followed by the cell returning to interphase rather than proceeding through meiosis. Such a process might also allow for some recombination among the progeny produced during intramolluscan development.
Although the preponderance of evidence is surely against the occurrence of meiosis, gamete formation or fertilization during intramolluscan development [88], there may be peculiar remnants of these processes represented, especially considering that most accounts of the evolution of digenetic trematodes favor the interpretation that the ancestral state was likely the sexually reproducing adult worm which was followed at a later time by the addition of asexual proliferative larval development in molluscs [93]. Might there then be peculiar remnant signatures of meiosis in intramolluscan larvae? We identified homologs to two known meiosis prophase-specific transcripts in our Illumina samples (Fig 8B), which were originally characterized in mice: meiosis express protein 1 (MEIG1) known to be involved in chromosome/chromatin binding in meiosis [93,94] and highly expressed during meiosis prophase 1  [93][94][95], and meiosis-specific nuclear structural protein 1 (MNS1). Anderson et al. [96] identified a MEIG transcript expressed in adult male and female S. mansoni with a potential role in gamete production but no possible functional role was suggested to explain its high expression in eggs. MNS1 is specifically expressed in mice during the pachytene stage of prophase 1 of meiosis. Retinoic acid (RA) initiates meiosis and although retinoic acid is not implicated in development of S. mansoni, we did see expression of retinoic acid receptor RXR. Of the putative meiosis stage-specific homologs, only RXR was present as a feature on the S. mansoni microarray, and it showed increasing expression as intramolluscan development progressed. Further study is warranted to learn if the transcripts we observed from intramolluscan stages are perhaps indicative of some tendency for occasional formation of bivalents without associated gamete formation or fertilization, or of a general repurposing of these molecules for use in many kinds of cellular reproduction, including asexual reproduction.
Six recombinase transcripts were expressed in our Illumina samples: three RAD51 homologs, two cassette chromosome recombinase b homologs, and one trad-d4 homolog. Recombinases like RAD51 are up-regulated in the testis and ovary of adult S. mansoni as compared to whole worm controls [96,97] and in female adult S. japonicum when compared to males [98]. Recombinases can repair breaks in DNA as a result of DNA damage or that occur during homologous recombination during meiosis. At least one transcript of another recombinase, topoisomerase II, was expressed in every time point. Among other functions, topoisomerase II interacts with the meiosis-specific RecA-like protein Dmc1 or RAD51 to facilitate pairing of homologous chromosomes during chromosome strand exchange [99].

Glycosyltransferase expression in intramolluscan stages
Molecular mimicry has been an area of interest with respect to schistosome-snail interactions since the early 1960s with the hypothesis that parasites express host-like molecules to evade host immune responses [100]. Several studies have highlighted antigenic similarities between miracidia and mother sporocysts and B. glabrata hemolymph proteins [101,102], and the glycans on S. mansoni glycoproteins and glycolipids have been extensively studied, including for their potential role in mediating host mimicry [20, [103][104][105][106][107]. Yoshino et al. [108] showed that antibodies to S. mansoni glycotopes bound more extensively to cell-free hemolymph (plasma) from snails susceptible to infection than plasma from resistant strains, and suggested hostmimicking glycotopes could be a determining factor in compatibility during early larval stages. Consequently, we were interested in examining S. mansoni glycosyltransferases because of the role they play in generation of glycan moieties on lipids and proteins. Among others, one group of glycosyltranferases we found to be prominent in S. mansoni intramolluscan stages were fucosyltransferases (FTRs). Several of the surface membrane glycoconjugates of S. mansoni that interact with B. glabrata are fucosylated [103,109] and are suspected to be involved in host mimicry [103]. We found 22 unique Illumina-assembled FTRs transcripts and 8 fucosyltransferase-specific probes represented on the microarray (S7 Fig). Fitzpatrick et al.
[23] observed two major clades of FTRs expressed by S. mansoni, those expressed in miracidia and mother sporocyst stages (alpha 1,6 fucosyltransferases D and E) and those expressed primarily in sexually mature adults (alpha 1,3 fucosyltransferases B, L, F). We did not observe an obvious stage-specific demarcation in FTRs expression but rather observed a broad range of FTR transcripts, including those Fitzpatrick et al.
[23] observed primarily in sexually mature adults. They were expressed at all time points with highest diversity being at 1d and in shedding snails. Alpha 1,6 fucosyltransferase H was expressed ubiquitously across all Illumina and microarray samples. We also saw expression of five O-fucosyltransferase transcripts exclusively at 1d and in shedding snails. O-fucosyltransferases add a fucose residue to the oxygen on a side chain of either serine or threonine residues in a glycoprotein. Heavy expression in shedding snails was not surprising because cercariae possess a prominent fucose-rich glycocalyx [110].
A transcriptional regulatory protein, KRAB-A domain-containing protein, and dolichyldiphosphooligosaccharide-protein glycosyltransferase subunit DAD1 homolog that performs post-translational protein glycosylation, were among the most abundantly expressed transcripts across all Illumina samples.

Sporocyst defense and stress responses
It is reasonable to expect that S. mansoni intramolluscan stages are under some duress from host immune responses, and we noted that snail Cu,Zn superoxide dismutases (SOD) were upregulated at both 1d and 3di in the highly compatible snail B. pfeifferi from which the S. mansoni Illumina transcripts discussed here were also obtained [17]. The H 2 O 2 resulting from SOD activity is known to be toxic to S. mansoni sporocysts and is a main factor responsible for killing early larval S. mansoni in some B. glabrata resistant strains [40,110,111]. Organisms can remove harmful intracellular hydrogen peroxide with catalases, glutathione peroxidases, and peroxiredoxins. Schistosomes lack catalases [112,113] and have low levels of glutathione peroxidases with limited antioxidant abilities [114]. It is suggested that peroxiredoxins are the schistosome's main defense against damage from hydrogen peroxide [115]. In our data, thioredoxin peroxidases, peroxiredoxins that scavenge H 2 O 2 using thioredoxin, are consistently (and highly) expressed throughout all time points in array and Illumina samples (Fig 9). Thioredoxin peroxidases reduce hydroperoxides with thioredoxin as a hydrogen donor. S. mansoni also expresses SOD activity and S. mansoni-encoded Mn-and Cu/Zn-type SODs are expressed in every replicate of both Illumina and array-sequenced samples (Fig 9). Because the level of S. mansoni SOD transcripts was consistently modest, we suggest that their function is not to mount an anti-snail counter-offensive but rather to take care of the intracellular anti-oxidative needs of the parasite.
Cytochrome P450 proteins have been associated with stress responses and in detoxification reactions. Schistosoma mansoni has but a single cytochrome P450-encoding gene and the associated protein activity has been shown to be essential for survival in both adult worms and eggs, although its underlying function in schistosomes remains unclear [116]. Cytochrome P450 transcripts showed minimal expression at 1d, were absent at 3d, and had modest (~3 Log 2 TPM) in all shedding replicates.
Heat shock proteins (HSPs) are often produced under conditions of stress, but they are also constitutively expressed in actively synthetic cells to serve as chaperones and to facilitate protein folding [116,117]. Schistosoma mansoni sHSPs 16 and 20 and HSPs 70 and 90 were all found among proteins released during in vitro miracidium to mother sporocyst transformation [20]. We found sHSPs 16, 20, and 40 and HSPs 60 and 70 but not 90, to be expressed throughout intramolluscan stages with the highest expression seen in HSPs 20, 40, and 70 ( Fig  9C). sHSP 20 contributes up to 15% of the soluble proteins of miracidia [118] and is a prominent protein identified in miracidia by LC-MS/MS [20,21]. sHSP 40 has been identified as a soluble egg antigen responsible for eliciting immunopathological reactions in the definitive host that result in granuloma formations [119]. Ishida and Jolly [120] showed that in the absence of HSP 70, cercariae do not orient or penetrate normally, providing additional functional roles for HSP 70 beyond stress responses. We noted that S. mansoni HSP 70 was expressed at high levels in cercariae-producing shedding samples.
In addition to defending themselves from attack by host immune components, in longlived host-parasite associations as represented by S. mansoni in B. pfeifferi, it might also be reasonable to expect digenean sporocysts to contribute to the stability and maintenance of the host-parasite unit as a whole. Rediae of some digenean species do this in the form of actively attacking newly-colonizing trematode infections [121]. In complex, natural transmission foci such as the one from which our samples originated, S. mansoni-infected B. pfeifferi snails are constantly exposed to a variety of viruses, bacteria, and infectious eukaryotes [17]. Therefore, it seems reasonable that while imposing considerable stresses on its hosts, possibly including immunosuppression, that larval schistosomes might also be expected to contribute to the wellbeing of the host-parasite unit by expressing transcripts that contribute to repression or elimination of additional parasites.
One way to gain insight into S. mansoni sporocyst capabilities in this regard was to review what is known for free-living flatworms, such as the planarian Dugesia japonica, that are regularly challenged by pathogenic and non-pathogenic bacteria in their habitats. Planarians are capable of phagocytosing and destroying pathogens like Staphylococcus aureus and Mycobacterium tuberculosis [122,123]. Conserved homologs to human genes, such as MORN2 (membrane occupation and recognition nexus-2 protein) are present and known to play a role in LC3-associated phagocytosis (LAP) elimination of bacterial pathogens in human macrophages and in flatworms. Homologs of 34 transcripts of putative flatworm anti-bacterial factors [123] were expressed in S. mansoni intramolluscan stages (Fig 9D). Homologs to dual specificity phosphatases were the most prominent group of anti-bacterial factors. They were expressed at all time points, with a general increase in expression with development time. Homologs of MORN2 (membrane occupation and recognition nexus-2 protein) were also present throughout intramolluscan development. MORN2 plays a role elimination of bacterial pathogens in human macrophages as well as in flatworms. MORN2 is present in all replicates in 1d, 3d, and shedding S. mansoni samples. Shedding S. mansoni samples in both Illumina and microarray contain the most flatworm bacterial defense homologs. How these putative schistosome defense factors may be deployed in sporocysts that lack a gut and that do not engage in phagocytosis as far as we know, remains to be seen.
Among possible S. mansoni anti-immune factors we noted to be highly expressed in our samples was calreticulin, previously shown to be present in the excretory/secretory products of S. mansoni sporocysts (Fig 10). Because of its calcium-binding capability, Guillou et al. [124] suggested calreticulin may interfere with hemocyte spreading and interfere with their ability to initiate encapsulation responses.

Evidence of amphistome-mediated suppression of S. mansoni sporocyst development
As noted in Buddenborg et al. [17], Illumina replicate 3d-R2 harbored a pre-patent infection of an amphistome species, presumptive Calicophoron sukari, known to be common in B. pfeifferi in Kenya, from the habitat from which these snails were obtained [125]. The presence of an amphistome infection is of interest because previous studies indicate that amphistomes and schistosomes interact in distinctive ways in the intramolluscan environment, with amphistomes having a permissive effect in enabling development of a schistosome that might not have otherwise developed in a particular snail species (e.g. [125,126]). In our context, the effect of the amphistome appears to be the opposite, based both on field results which suggest that amphistome infections supplant S. mansoni infections (Laidemitt, personal communication), and on our transcriptional results. In general, overall S. mansoni transcription in the 3d replicate with the amphistome was dampened relative to 3d replicates lacking the amphistome. This dampening took the form of both fewer numbers of S. mansoni transcripts expressed, and for those that were expressed, representation at lower copy numbers. It did not appear that specific highly expressed S. mansoni transcripts were targeted in any selected way by the presence of the amphistome. In fact, the S. mansoni transcriptome in the amphistome-containing replicate most closely resembles that seen for S. mansoni in 1d infections (1d-R1), suggestive of an amphistome-imposed delay in development. Another amphistome species, Calicophoron The intramolluscan transcriptome of Schistosoma mansoni daubneyi co-occurs with Fasciola hepatica in the intermediate host snail Galba truncatula [127,128] bringing forth interesting parallels and implications for control of helminths in ruminants.

Concluding comments
As a miracidium penetrates a snail, it rapidly enters a radically different milieu from what it has previously experienced, and a number of pre-made proteins are released into its new surroundings to effect transition to the mother sporocyst stage adapted for intramolluscan existence [20,21,40]. Our 1d Illumina samples include many transcripts distinctive from the proteins associated with transformation, indicative of the switch to the needs of existence as sporocysts. As examples, we found representation of different SmVALs, heat shock proteins, protease transcripts (elastases), neural development proteins or neurohormones among 1d Illumina samples than seen as proteins in miracidial transformation products.
From our earliest Illumina samples, it is evident that S. mansoni orchestrates a complex transcriptional program within its snail hosts, with a significant percentage of its genetic repertoire (an estimated 66% of the S. mansoni genome) engaged (see also [18] who reports that 50-60% genes are expressed in each S. mansoni stage). This is particularly evident at the stage of production of cercariae, when large amounts of parasite tissue are present and production and differentiation of the relatively complex cercarial bodies are underway. Also noteworthy is that a core transcriptome required of life in a snail host can be identified which includes transcripts involved in a central glycolysis pathway and the TCA cycle, and for transmembrane transporters for monosaccharides, amino acids, steroids, purines and pyrimidines, indicative of the dependence of sporocysts on their hosts for key molecular building blocks. Several of these like the amino acid transporters are noted for the first time in schistosome sporocysts. Metabolically, early stage schistosomes focus primarily on acquisition of molecular building blocks and nutrients with a distinct switch from storage to expending these components towards cercariae production in patent stage shedding S. mansoni. A role for receptor-mediated endocytosis in sporocyst nutrition should also not be excluded [129]. Transcripts for enzymes required for macromolecular synthesis and cell proliferation, the latter a prominent and perpetual part of intramolluscan development and cercarial production, were also part of the core transcriptome. Oxidoreductase activity and cell redox homeostasis are among the most abundant functions across all larval stages of S. mansoni.
With respect to particular S. mansoni transcripts that may be key to successful intramolluscan development and/or that comprise parts of the "interactome" with snail transcripts, we highlight several findings below. First, by virtue of using both field-derived B. pfeifferi and S. mansoni from infected children, our Illumina study allows for a broader range of outcomes particularly as measured in the early stages of infection, and the variability that results can provide distinctive insights. For instance, cases where early sporocysts seem to be thriving with respect to read count are accompanied by production of larger quantities of factors associated with infectivity like proteases, fucosyltransferases, SmVALs, and GPCRs. Also, poor transcriptomic productivity for S. mansoni sporocysts has been associated with presence of other digenean species, unknown to us to be present at the time of exposure to S. mansoni, but that interfere with S. mansoni development [126,129] (Laidemitt, personal communication).
One of the surprising things about the genome of S. mansoni and of other parasitic helminths is the dearth of genes such as cytochrome P450s involved in degradation of xenobiotics (S. mansoni has but one such gene with unknown function), potentially including harmful snail-produced factors as well. The ABC transporters we have noted to be highly expressed in sporocysts may function to compensate [48].
Of particular interest was the expression of a diverse array of proteases and protease inhibitors at all intramolluscan stages, including some proteases like elastases characteristic of cercariae that were also produced by early sporocyst stages. The up-regulation by the snail host of protease inhibitors during the larval establishment period at 1d and 3d [17] seems a likely response to prevent parasite establishment. Also of note was expression of S. mansoni protease inhibitors that might inhibit the action of snail trypsin-like proteases up-regulated late in infection [17]. These protease inhibitors may prevent the activation of the phenoloxidase enzyme laccase, whose activity induces a negative impact on late (7-9 week) S. mansoni infections [17,70].
GPCRs were also well represented in S. mansoni intramolluscan stages and are likely to play several important roles in schistosome development. One such GPCR is expressed at 1d and in shedding snails and is known to bind serotonin [48,130]. We noted at the same time points that B. pfeifferi down-regulated production of a type 1 serotonin receptor and additionally, at 3d kynurenine 3-monooxygenase, which degrades tryptophan and can limit concentrations of serotonin, is up-regulated [17]. It seems reasonable to continue to suspect serotonin of playing a role in parasitic castration. It stimulates egg production when given to castrated snails [131]. By expressing the appropriate serotonin GPCR, and possibly down-regulating the host receptor, S. mansoni sporocysts may limit availability of serotonin to the snail.
Other factors also worthy of additional consideration with respect to parasitic castration are S. mansoni neuropeptides Y and F and their receptors which were expressed particularly in shedding snails. In snails, neuropeptide Y has been associated with decreased egg production [132] and neuropeptide Y receptor was up-regulated at 1 day and in shedding B. pfeifferi snails [17]. Whether these neurotransmitters produced by S. mansoni might directly affect snail reproduction is not known. Ovipostatins which have a suppressive effect on egg laying in the snail Lymnaea stagnalis were found to be up-regulated in shedding B. pfeifferi so it is possible their expression may be targeted by S. mansoni in some manner as well. We did not see obvious changes in some snail neuroendocrine factors associated with reproduction like calfluxin or schistosomin [17]. Wang et al. [85] in a proteomic study of neuropeptides from B. glabrata, including snails with 12 day infections with S. mansoni, found lower levels of many snail reproductive neuropeptides. The extent to which S. mansoni and other digenetic trematodes might effect snail reproduction through interference with their neuroendocrine systems as proposed by de Jong-Brink [132] remains a topic worthy of more study. As noted by Humphries [15], it is also possible that castration is more a consequence of depletion of nutrients and alterations of metabolism imposed by metabolically demanding larval schistosomes.
Insights provided by the study of planarians were important to our interpretation of our results in two ways. The first was to confirm in intramolluscan S. mansoni samples the common expression of genes associated with maintenance of stem cell-like germinal cells, including fibroblast growth factor receptors (fgfr), vasa, argonaute2 (ago2), and nanos-2 [15]. The second was to examine sporocysts for evidence of homologs of transcripts known to be involved in antibacterial responses in Dugesia japonica [123], for which 34 were found. Whether these factors are actually deployed in anti-bacterial or other defense responses remains to be seen. Their presence is somewhat peculiar because sporocysts lack phagocytic activity, unlike the gut cells of planarians. However, perhaps anti-bacterial proteins are deployed along sporocyst membranes, or sporocysts may engage in limited forms of endocytosis [129] that might result in engulfment of bacteria or their products. Possible anti-snail hemocyte factors like calreticulin [124] were also expressed by sporocysts. The possibility that sporocysts contribute to discouraging or preventing growth of third party symbionts that could compromise the snail-schistosome functional unit, especially in light of the need of schistosome sporocysts to compromise components of host immunity is also a topic worthy of additional study.
One advantage of next gen sequencing is its potential to provide unexpected insights. We were surprised to see two transcripts (MEIG1 and MNS1) associated specifically with prophase I of meiosis and discuss possible interpretations based on previous cytological studies of germinal cell development [92], observations that might help to explain differences among S. mansoni cercariae in genetic content [90]. The expression of recombinases in sporocysts might be consistent with a partial entry into meiosis up to diakinesis, or with mitotic recombination, the latter suggested by Grevelding [90] to account for genetic differences among cercariae derived from the same miracidium.
Finally, we note that many more highly represented transcripts were found, including those encoding both genes with suspected or unknown functions whose connection with intramolluscan development remain to be elucidated. With ever more complete transcriptional profiles becoming available for schistosomes in their snail hosts, the stage is set for further studies employing the best tools available for gene knockout to address the functional roles of these and the many other transcripts we and others have noted. Of particular interest will be to determine if ingenious use of this information can be made to specifically target and prevent the development of sporocysts and their production of human-infective cercariae, thereby opening a much-needed additional front in the effort to control and eliminate human schistosomiasis.
Supporting information S1 Table.