Microarray analysis of gene expression induced by sexual contact in Schistosoma mansoni
BMC Genomics volume 8, Article number: 181 (2007)
The parasitic trematode Schistosoma mansoni is one of the major causative agents of Schistosomiasis, a disease that affects approximately 200 million people, mostly in developing countries. Since much of the pathology is associated with eggs laid by the female worm, understanding the mechanisms involved in oogenesis and sexual maturation is an important step towards the discovery of new targets for effective drug therapy. It is known that the adult female worm only develops fully in the presence of a male worm and that the rates of oviposition and maturation of eggs are significantly increased by mating. In order to study gene transcripts associated with sexual maturation and oviposition, we compared the gene expression profiles of sexually mature and immature parasites using DNA microarrays.
For each experiment, three amplified RNA microarray hybridizations and their dye swaps were analyzed. Our results show that 265 transcripts are differentially expressed in adult females and 53 in adult males when mature and immature worms are compared. Of the genes differentially expressed, 55% are expressed at higher levels in paired females while the remaining 45% are more expressed in unpaired ones and 56.6% are expressed at higher levels in paired male worms while the remaining 43.4% are more expressed in immature parasites. Real-time RT-PCR analysis validated the microarray results. Several new maturation associated transcripts were identified. Genes that were up-regulated in single-sex females were mostly related to energy generation (i.e. carbohydrate and protein metabolism, generation of precursor metabolites and energy, cellular catabolism, and organelle organization and biogenesis) while genes that were down-regulated related to RNA metabolism, reactive oxygen species metabolism, electron transport, organelle organization and biogenesis and protein biosynthesis.
Our results confirm previous observations related to gene expression induced by sexual maturation in female schistosome worms. They also increase the list of S. mansoni maturation associated transcripts considerably, therefore opening new and exciting avenues for the study of the conjugal biology and development of new drugs against schistosomes.
Schistosomiasis is an important public health problem that affects approximately 200 million people, mostly in developing countries, and that poses a risk to another 600 million . Mortality attributed to schistosomiasis is estimated to be around 11,000 deaths/year with an added economic burden of 1.7 million disability adjusted life-years lost per year . The disease is transmitted when parasite eggs in human feces reach fresh water and hatch into miracidia. Upon infection of an appropriate snail species, miracidia transform into sporocysts that asexually reproduce generating hundreds of thousands of cercaria. The cercaria are shed by the snail and swim until they find an appropriate vertebrate host and penetrate through its skin. The worms then undergo differentiation and migrate in the bloodstream until they reach the mesenteric veins where male and female worms pair.
Eggs from paired mature schistosomes are responsible for most of the pathology caused by S. mansoni. Each mature female lays on average 300 eggs per day, some of which are excreted in the feces . The remaining eggs end up in the liver, intestine and other organs where they cause an inflammatory reaction, producing significant scarring which leads to a variety of symptoms depending on the organ where they lodge. For instance, eggs trapped in the bowel wall may cause bloody diarrhea, cramping, and eventually inflammatory colonic polyposis . Eggs that are swept back to the hepatic portal system cause granulomatous reaction in the portal tract which can evolve to hepatoportal fibrosis and portal hypertension. More interestingly, although immature worms may lay eggs which are occasionally shaped normally, those eggs are unable to induce the formation of granulomas .
Male schistosomes are responsible for triggering and maintaining female maturation. In absence of the male, female worms cannot migrate against the blood flow from the portal sites in the liver to the smaller mesenteric circulation where they lay their eggs . Therefore, the survival of S. mansoni couples and the maintenance of their complete life-cycle seem to be dependent on the existence of a permanent association between sexes. Studies by Shaw and Erasmus on praziquantel (the drug of choice for treating Schistosomiasis) have shown some evidence that the drug disrupts the reproductive system of females and that subcurative doses cause a long-lasting regression of both ovary and vitelline gland 24 hours post-treatment . Popiel et al. (1984) have shown that there is a dose-dependent and reversible effect of oxamniquine (another drug which is used for treating S. mansoni infections) on the female reproductive system [8, 9]. Although these effects on the reproductive system are usually reversible, these data might suggest that the reproductive system of female schistosomes could provide good targets for therapeutic agents against these parasites.
In contrast to male schistosomes that undergo normal morphological development regardless of their pairing status, females from single sex infection (single-sex females) show clear differences from those which are paired to males. Virgin females are considerably smaller than paired females . When paired females are separated from males they stop laying eggs and regress to an immature state, however if they are allowed to couple they maturate again [9, 11–13] showing that the maintenance of sexual maturation in the female is dependent upon the presence of the male partner. The stimulus for female maturation is independent of male sperm, species and fertilization [14–17]. In addition to the morphological modifications present in single-sex females, separation of female worms leads to a large variety of changes in the physiology of the worm, culminating in the atrophy of the ovary and vitelline gland . For instance, there is a great increase in DNA synthesis in single-sex females upon pairing with males. Also, paired females utilize more glucose than unpaired ones . In spite of their immaturity, single-sex females reach a limited level of development, allowing them to lay eggs, which are usually malformed and non-viable [5, 14]. The rates of oviposition and maturation of eggs are significantly increased by pairing of the single-sex female worms with males .
Different factors have been proposed to be involved in S. mansoni sexual maturation including physical/tactile contact [11, 14, 15], nutrition [6, 19–21] and chemical stimuli [16, 22, 23]. In support to the importance of tactile contact, Michaels showed that when halves of worms (i.e. one of the two segments from worms that were cut in two halves) are allowed to mate in vitro, like halves (either cranial or caudal halves) always mate normally while unlike halves almost always mate abnormally, postulating that both sexes of worms have linear receptors which are used to determine dating and mating position . In support to the role of nutrition, S. mansoni female worms from single-sex infections have limited pharyngeal musculature, thin intestinal cecal walls, lack of digestive enzymes and reduced intake of red blood cells than females from paired infections . LoVerde and colleagues have suggested that the muscularity of the male worm allows it to assist the weaker females with pumping host blood into their intestine by a massaging effect of the muscular walls of the gynecophoric canal (i.e. a ventral longitudinal groove in the surface of male Schistosomatidae) . Also, Cornford and Huot demonstrated that there is transfer of 14C glucose from male worms to females, implying that male schistosomes indeed help feed female worms . Chemical stimuli such as acetone or ether extracts of male worms have been shown to induce female development in vitro . Also, an observed cholesterol transfer from males to females may suggest a role for hormonal factors .
On the male parasite side, the effects of sexual pairing on male worms are very subtle and therefore have received little attention in comparison to female worms. Nevertheless, some studies show that male S. mansoni and S. haematobium worms from single-sex infections are significantly smaller and have a reduction in the number of testes when compared to worms from homologous paired infections [24, 25]. Using the data published by Basch and Gupta it is possible to observe that pairing also stimulates a higher rate of cell division (as measured by the number of nuclei per section unit area) in males when bisexually paired males and unisexual unpaired males are compared . The stimulatory effect of female schistosomes on male worms has also being studied at the biochemical and molecular level. For instance, worms from single-sex infections have more glutathione than worms from paired-infections and this difference can be reversed by incubation of those worms with females [27, 28]. Also, female worms stimulate tyrosine incorporation, lipid accumulation and lipase utilization in males [23, 29]. Another interesting fact is that paired male schistosomes express a 79-KDa protein in its gynecophoral canal and that this protein detection is severely limited at the surface of single-sex male worms . This protein, which shows wide distribution on the surface of adult female worms has homology to developmentally-regulated homotypic adhesion molecules and is suspected to be essential for continued interaction between sexes and sexual development of schistosomes.
Several expression analyses were recently performed in schistosomes using microarrays to study gender differences in gene expression of S. mansoni and S. japonicum, as well as stage specific, strain specific, maturation specific, and species specific differences [31–38]. Microarrays were also used to investigate the vertebrate hepatic host response to infection with S. mansoni . These studies have provided important information concerning the biology of the parasite and the host response. In this study, we describe the use of DNA microarrays to study the effect of sexual pairing of S. mansoni (Brazilian LE Strain) in the adult gene expression. To accomplish this objective we compared worms from single-sex infections with worms from paired-infections. Our results reveal novel information of genes putatively involved in oogenesis and sexual maturation of schistosomes and allow the identificaton of new possible targets for therapeutic intervention.
Results and discussion
Evaluation of replication quality
In order to investigate the quality of our replicates we tested the correlation between different types of replicates using LIMMA (Linear Models for Microarray Data) . To generate single-sex worms we exposed multiple Biomphalaria glabrata snails each to a single miracidium and harvested the cercariae shed by them after one month. Multiple mice were subsequently infected with either single-sex or mixed-sex cercariae and the worms harvested 51 days post-infection. The experimental design included three types of replicates: biological, technical and in-slide replicates encompassing 6 slides in total for each experiment. In brief, single-sex and paired-infection worms were harvested from multiple mice and pooled separately in pools of 500 worms or more (each pool contained worms from different mice). The RNA from each pool was extracted separately and then hybridized to the slides. Three biological replicates were used for each experiment and each of them was hybridized twice (dye swap). Because each probe was printed in duplicate we computed also the correlation between in-slide replicates.
The approach used was to estimate a common correlation for all the genes within blocks corresponding to each type of replicate. As would be desired, the level of correlation between in slide technical replicates was between 0.90 and 0.92 for male and female experiments respectively. Dye swap replicates displayed a correlation between -0.85 (male experiment) to -0.89 (female experiment), indicating that the fluorescent dyes introduce a very minimal bias in hybridization intensities. The correlation between our biological replicates (using RNA isolated from different pools of worms) was 0.73 for the female experiment and 0.61 for the one investigating sexual pairing on male worms. Our results are in agreement with the data published by Hoffmann and colleagues who found an in-slide replicate correlation of 0.91 and a correlation ranging from 0.81 to 0.90 between slides (biological) [31, 41].
Differentially expressed genes
Statistical analysis of the data using a linear model to assess differential expressed genes revealed a large number of transcripts of gene products putatively related to sexual maturation in both male and female schistosomes (Table 1 and 2). Only genes with at least two-fold difference (log2 ratios ≥ 1 or ≤ -1) and logOdds ≥ 4.6 were considered differentially expressed. Volcano plots, where the magnitude of the gene expression ratio is displayed on the x-axis and the significance of the difference in expression between groups on the y-axis, are shown in (Figure 1). 265 genes were observed to be differentially expressed between mature and immature female parasites and 53 genes between mature and immature male worms.
Our results show that 55% of the genes which are differentially expressed are expressed at higher levels in paired females while the remaining 45% are more abundant in immature female worms. The differentially expressed transcripts (Table 1) in female worms included many previously identified gender related genes including the p14 eggshell protein gene (Contig 1625, SCMESP), other egg shell proteins (Contig 6116, Contig 1557, SMEGGPRO, Contig 867), extracellular superoxide dismutase (Contig 119, SCMSOD2, SEG_SCMOD, SCMSOD3, SCMSODM), tyrosinase (Contig 948, Contig 1617), cathepsin B (SMA312106), dynein light chain (SMU55992), gynecophoral canal protein (SMU47862), aspartate aminotransferase (Contig 1636), diacylglycerol acetyltransferase (Contig1301) and female specific protein fs800 (SCMFS800) [41–43].
Fs800, has already been shown to be up-regulated in the vitelline cells during maturation of female worms and has been linked to the production and maintenance of eggs . The gynecophoral canal protein (SMU47862), also previously shown to be differentially expressed between male and females, has been implicated in male/female communication or interaction . Our data showing down regulation of the small G protein Ras in female is also consistent with published observations . Tyrosinase, extracellular superoxide dismutase, glucose transporter protein I, diacylglycerol acetyltransferase 2-like protein, aspartate aminotransferase and egg shell proteins are known to be differentially expressed between male and females suggesting that they might be related either to sexual determination, maturation, reproduction or other sex specific activities [31, 41]. Our results show that these genes, in addition to being differentially expressed between genders, are also differentially expressed in female worms in different sexual maturation stages. A great number of transcripts not previously associated with sexual maturation were also observed. These included membrane antigens, importin beta-1, S. mansoni glucose transporter, myosin heavy chain, rho2 gtpase, lecithin-cholesterol acytransferase, apoferritin-2, various ribosomal proteins and others (Table 1). Of the 265 genes we identified as differentially expressed in females, 142 did not contain any annotation because they did not show significant similarity with sequences contained in the NCBI at the time the array was printed . These genes represent a multitude of unexplored proteins which might be potential targets for new therapeutic agents.
Of the genes with an associated gene ontology, genes that were up-regulated in single-sex females were mostly related to energy generation (i.e. carbohydrate and protein metabolism, generation of precursor metabolites and energy, cellular catabolism, and organelle organization and biogenesis) while genes that were down-regulated were related to RNA metabolism, reactive oxygen species metabolism, electron transport, organelle organization and biogenesis and protein biosynthesis (Figure 3). The increased amount of antioxidants in mature female is in agreement with their increased need to detoxify hemoglobin byproducts associated with increased blood cell consumption . On the other hand, one would expect that paired-infection females express more carbohydrate and protein metabolism genes than single-sex worms. The fact that we found more of those genes in single-sex females supports Basch's nutritional theory of sexual maturation. If males do indeed help females to feed, this could allow mature worms to direct their metabolism to the more important egg generation (since the male is already taking care of energy generation) and therefore genes related to energy generation should be decreased in mature females in relation to immature. Among the genes that are down-regulated there are six genes coding for eggshell proteins. This is expected since paired females are egg-laying machines capable of converting the equivalent of their own body dry weight into eggs  and consuming about 8 millions erythrocytes per day. In light of this high metabolic activity of mature females, it has been proposed that female-specific gene products associated with the metabolic machinery of egg production and hemoglobin catabolism should be highly represented in molecular biological studies of gender  which is the case in this study. The fact that genes associated with organelle organization were both up and down regulated in immature females could possibly be explained by differences in vitelline cells production among mature and immature females. It has been shown that paired female schistosomes produce several thousands vitteline cells per day . This high mitotic activity must demand high amounts of protein synthesis and the participation of an array of proteins related to organelle biogenesis and organization. The fact that some genes related to this process were up-regulated in immature females is more difficult to explain but it could be due to some compensatory mechanism whereby immature females make up for the lack of nutritional support by males by producing new cells.
For male parasites, our results show that 56.6% of the differentially expressed genes were present in larger quantities in mature males while the remaining 43.4% are more expressed in males from single-sex infections. The fact that more genes were differentially expressed between females than between males is consistent with the fact that immature males present less morphological and biochemical differences than immature females. Nevertheless this data does not agree with a previous study of gene expression induced by sexual maturation which found a similar number of differentially expressed genes in males and females . This difference might have been caused by differences in the experimental design, data analysis parameters, host characteristics or differences in the parasite strain used. Indeed, results from Moertel and colleagues have shown that more than 600 genes are differentially expressed when different strains S. japonicum are compared under identical conditions . When the gene lists of our experiments and Fitzpatrick and Hoffmann  are compared (See Additional file 1) we find less than 10% of genes which are common to the two datasets and in one case (single-sex males) we did not find any genes in common. Our data also shows that mature male worms express more genes related to the metabolism of proteins, nucleotides, macromolecules, ion transport, and metabolite/energy generation than unpaired worms, which suggests that paired worms are more metabolically active. The data agrees with Fitzpatrick and Hoffmann (2006) in that more genes potentially involved in RNA metabolism were differentially expressed between mature and immature male worms , nevertheless we did not find a greater proportion of genes in the structural molecule activity GO category.
Interestingly, the data from Khalil and Mansour show that male worms from single-sex infections are in fact smaller than worms from paired-infections . Although the authors do not explain this phenomenon, our results suggest that the difference in size could be linked to the higher expression of genes involved in energy generation and accumulation (Figure 3). The fact that paired male worms absorb and phosphorylate glucose more rapidly than unpaired worms also reinforces this hypothesis . Nevertheless, we observed a paradoxical higher expression of Cathepsin B Sm31 in unpaired male worms. Gotz and Klinkert demonstrated that this protein is an enzyme which is involved in the degradation of hemoglobin in the digestive tract of schistosomes .
Another interesting observation is that unpaired male worms synthesize more lecithin-cholesterol acyltransferase (LCAT), acyl-CoA binding proteins (ACPB) and fatty acid binding protein (FABP). In humans LCAT is responsible for the formation of cholesteryl esters starting from lecithin and cholesterol and it is the LCAT on the surface of nascent HDL particles that convert cholesterol and phosphatidylcholine of chylomicron and VDLD in cholesteryl esters. The accumulation of cholesteryl esters in the nascent HDL convert it in mature HDL and directs it back to the liver where the cholesterol is usually unloaded. The acyl-CoA binding proteins are small proteins which bind to long and medium chain acyl-CoA esters with high affinity and which can act as intracellular carriers of acyl-CoA esters. FABPs, are a members of a family of cytosolic lipid binding proteins and that have been investigated as possible vaccine candidates . Furthermore, Gobert and colleagues demonstrated that the S. japonicum FAPB is confined to lipid droplets in male worms that are probably involved with the nurture of female worms .
In studies with humans and marmosets (Callithrix jacchus) the plasmatic levels of total cholesterol, esterified cholesterol, triglycerides and phospholipids were significantly reduced after infection with S. mansoni [52, 53]. Similar effect was observed by Doenhoff and colleagues, who reported a reduction in atherosclerosis in rats predisposed to cardiovascular disease (apoE-/-) after infection with S. mansoni . They suggested that this reduction could be due to a modulatory effect of the infection by S. mansoni on lipid metabolism in the host. On the other hand, Ramos and colleagues proposed that the cholesterol reduction could be caused by changes in the LCAT and/or ACAT (acyl-CoA-cholesterol acyl transferase) activity, synthesis or secretion produced by the liver of infected animals . Taking into consideration that schistosomes are not capable of synthesizing cholesterol [55, 56], the fact that sexual pairing has an effect upon so many genes related to fatty acid and cholesterol metabolism in male worms becomes even more relevant.
Validation of microarray results using real time quantitative PCR
In order to validate the microarray results we performed real-time RT PCR. Contigs considered to be differentially expressed were first checked against the S. mansoni genome to evaluate the validity of the probes printed on the arrays and to remove redundant contigs (i.e. different contigs which represented the same gene).
Primer pairs were designed for all genes considered as differentially expressed and then tested against the genome using electronic PCR  and BLAST  to select only those that amplified a unique region of the genome. 41 female and 20 male primer pairs qualified, of which approximately 20 of each were randomly selected for real time quantitative PCR. The results from the quantitative PCR for the male and female experiments were compared against the results of our microarray hybridizations using linear regression analysis  and showed a strong correlation (R = 0.84 and R = 0.67 respectively).
The results presented on this article confirm previous results related to gene expression induced by sexual maturation in schistosome worms. It also contributes a wealth of information about genes (either characterized or unannotated) that may be involved in this process. Since most of the differentially expressed genes do not have homology to other genes, further work is necessary to characterize these transcripts and assign biological functions to them. Once the Schistosoma mansoni genome is fully annotated, it will be possible to extract more interesting information from our data. Nevertheless, we believe that the myriad of up- and down- regulated genes we describe here open new and exciting doors for the study of the conjugal biology and the development of new drugs against schistosomes.
S. mansoni lifecycle was maintained at Instituto de Pesquisa Rene Rachou, Fiocruz (Brazil). S. mansoni (LE Brazilian strain) is maintained in Biomphalaria glabrata as intermediate host. Outbred Swiss mice were used as definite hosts. To obtain single-sex adult worms, B. glabrata snails were exposed to a single miracium, which was generated from eggs obtained from the liver of infected mice . After one month the snails were tested for positive infection by verifying the shedding of cercaria following exposure of the snails to artificial light . 100–200 cercaria from each snail were injected in the peritoneum of female Swiss mice.
Mice were housed conventionally in polypropylene cages with stainless steel screen covers. All animals received lab mouse chow and water ad libitum. The animals were sacrificed at 51 days post-infection according to ethical procedures and adult worms were obtained by portal perfusion . The worms were washed in cold saline solution and checked by microscopy for the presence of possible undesirable mixed-sex infections. We separated the single-sex adult worms in multiple pools (each one containing >500 worms which were originated from multiple mice) which were frozen at -80°C until further processing of the samples.
To obtain worms from paired-infections, B. glabrata snails were exposed to multiple miracidia and the cercaria from these snails used to infect Swiss mice. Each mouse received approximately 100 cercaria and the animals were sacrificed at 51 days post-infection. The perfusion was performed in the same fashion as for single-sex worms. The worms were washed with cold saline solution and carefully separated by their sex with forceps under the microscope. Worms from each sex were pooled separately (>500 worms/pool) and frozen until further processing.
RNA extraction and amplification
Total RNA was extracted using Trizol reagent (Invitrogen Life Technologies, CAT#15596-026) according to the manufacture's instructions. The RNA was quantified using a Nanodrop ND-100 UV/Vis spectrophotometer (NanoDrop Technologies, USA) and the overall RNA quality was assessed using denaturing gel electrophoresis . Two μg of total RNA from each sample were amplified by doing two rounds of linear amplification using the Amino Allyl MessageAmp II kit (Ambion, CAT#1753). The amplification was done according to manufacture's specifications and aaUTP was used on the second round of amplification so that the final product could be labeled using indirect labeling.
RNA labeling and hybridization
Aminoallyl amplified RNA was labeled using Cy3 and Cy5 by indirect labeling according to a modified version of the TIGR's standard operational procedure . In brief, for each hybridization 15 μg of amplified RNA was dried in a speed-vac, resuspended in pH 9.0 carbonate buffer and incubated for 5 minutes at RT. After this incubation, 4.5 μL of one of the two dyes dissolved in DMSO (Amersham, CAT#PA23001 and PA25001) was added to the solution and the samples were incubated at RT in the dark for one hour. Labeled RNA was purified away from unincorporated dye using RNeasy MinElute Cleanup Kit (Quiagen, CAT#74204). The Cy3 and Cy5 labeled samples were then combined, dried again, and resuspended in hybridization buffer (50% formamide, 5× SSC, 0.1% SDS).
Microarray hybridization and experimental design
The pools of female or male worm RNA from single-sex infections were hybridized against RNA from paired-infections from worms of the same sex. Each experiment was performed in duplicate (technical replicate) using a dye swap design in order to account for dye biases. The samples were hybridized using oligonucleotide DNA microarrays obtained from the laboratory of Karl Hoffmann at the University of Cambridge, UK. The arrays contained 7335 oligonucleotides (50 mers) spotted in duplicate. Each oligonucleotide represented either a singleton (a single EST sequence) or a contig (contigous sequences of overlapping EST sequences) generated by assembly with CAP3. The oligonucleotides were designed based on transcriptome information available at the time . Briefly, the slides were pre-hybridized by placing them in coupling jars containing pre-hybridization solution (5× SSC, 0.1% SDS, 1%BSA) for 40 minutes at 42°C. The slides were then washed by dipping 10 times in a beaker containing DI water, the water was changed and the operation repeated once. The slides were spun dry using a table-top high speed microarray centrifuge (TeleChem International Inc., USA). The samples (resuspended in hybridization) were hybridized overnight under cover slips inside Corning® hybridization chambers (Corning, USA) which were kept in a water bath at 42°C in the dark. Slides were washed two times for five minutes each in low stringency wash at 42°C (2× SSC, 0.1% N-lauroysarcosine), followed by two washes for 5 min in medium stringency wash (0.1× SSC, 0.1% N-lauroysarcosine) at RT and 2 washes for 5 minutes each in high stringency wash solution (0.1× SSC). Slides were spun dry and scanned using a microarray dual channel laser scanner (GenePix 4000B, Molecular Devices, USA) at 10 μm resolution, 100% laser power and PMT levels which were adjusted in order to obtain similar distributions of red and green signal intensities.
Real time RT-PCR
A subset of genes predicted to be differentially expressed was selected for validation using real-time RT-PCR. Total amplified RNA (2 μg) from mixed sex infection and single sex infection was used for reverse transcription using TaqMan Reverse Transcription Reagents (Applied Biosystems; Cat# N808-0234). Products were amplified using the Applied SYBR Green Masters Mix kit (Applied Biosystems, Cat#4309155) in an ABI Prism 7900HT Sequence Detection (Applied Biosystems) with the following profile: 50°C for 2 min, 95°C for 10 min; 40 cycles of 95°C for 15 s and 60°C for 1 min. Each reaction was performed using 1 μL of cDNA from the RT reaction using a final volume of 20 μL (PCR Master Mix 1×, 200 nM of each primer). Expression levels of S. mansoni alpha-tubulin (accession number M80214) were used as endogenous control within each sample. Relative levels of gene expression were calculated using the 2-ΔΔCT method . Each sample was analyzed for primer dimer, contamination or mispriming by inspection of their dissociation curves.
Overall, analysis of sexual pairing induced gene expression encompassed 6 slides (for each experiment), incorporating 3 replicates and dye swaps. The data, which are MIAME compliant , were submitted to ArrayExpress at EBI using MIAMExpress . Spots were analyzed using GenePix Pro and flagged according to their quality (i.e. spots were flagged as 'Bad' whenever they were contaminated with particles, were smeared or dilated, had irregular shape or were in areas of high background; otherwise they were flagged as 'Good'). Raw intensity data were analyzed using the R statistical language . The data were inspected for spatial biases on both the red and green channel (background and signal), for print-tip bias, dye bias and intensity dependent bias using the maArray package . The data were then print-tip normalized using a modified version of the robust spline method [40, 70]. The statistical analysis was performed using a linear model incorporating replication information. The in-slide, dye swap and pool's replicate correlations were calculated using the duplicate Correlation function of LIMMA . A list of differentially expressed genes was generated by applying a Bayesian smoothing to the linear model fit. Genes that had log odds > 4.6 (99% probability that the gene is differentially expressed between the conditions being compared) and M values >1 or <-1 (2 fold difference among groups) were considered as significant. Because the slides were designed based on the assembly of ESTs from S. mansoni and because the genome wasn't available at the time the slides were printed, we decided to test the differentially expressed genes against the complete genome. In brief, pairs of primers were designed for all differentially expressed genes (M values <1 or <-1 and logOdds > 4.6) using a web tool (Genscript Corp, USA ) and then tested against the genome and the predicted complete CDSs using mePCR  in order to see if each primer pair would theoretically amplify a single region of the genome. Primer pairs that amplified more than one region were discarded. We further tested the primers to see if two primer pairs would amplify the same gene (predicted CDS) in order to remove redundant genes (i.e. genes that were spotted more than one time at the arrays because of assembly problems).
Further analysis was performed using Blast2GO . This software allows the evaluation of differences in annotation between two groups of data. The analysis of GO terms association for the genes considered differentially expressed in both experimental groups, for the two experiments, was performed by using the combined graphs function of the software. These graphs allow the visualization of the combined annotation of a group of sequences and this can be used to study the biological meaning of a subgroup of sequences. Combined graphs are a good alternative to enrichment analysis because they don't require a reference group and allow the use of small numbers of sequences .
Acyl-CoA-cholesterol acyl transferase
Acyl-CoA Binding Protein
Bovine Serum Albumin
Disability Adjusted Life Years, the sum of years of potential life lost due to premature mortality and the years of productive life lost due to disability.
European Bioinformatics Institute
Expressed Sequence Tags
Fatty acid binding protein
High Density Lipoprotein
Low Density Lipoprotein
Linear Models for Microarray Data
Minimum Information About a Microarray Experiment
Polymerase Chain Reaction
Sodium Chloride Sodium Citrate
Very Low Density Lipoprotein
World Health Organization
TDR: TDR Strategic Direction: Schistosomiasis. 2002, World Health Organization
World Health Organization.: The world health report : Mental Health : New Understanding, New Hope. 2001, Geneva , World Health Organization, v.-
Loverde PT, Chen L: Schistosome female reproductive development. Parasitol Today. 1991, 7 (11): 303-308. 10.1016/0169-4758(91)90263-N.
Behrman AJ: Schistosomiasis. eMedicine. [http://www.emedicine.com/emerg/topic857.htm]
Shaw MK: Schistosoma mansoni: vitelline gland development in females from single sex infections. J Helminthol. 1987, 61 (3): 253-259.
LoVerde PT, Niles EG, Osman A, Wu WJ: Schistosoma mansoni male-female interactions. Can J Zool. 2004, 82 (2): 357-374. 10.1139/z03-217.
Shaw MK, Erasmus DA: Schistosoma mansoni: praziquantel-induced changes to the female reproductive system. Exp Parasitol. 1988, 65 (1): 31-42. 10.1016/0014-4894(88)90104-X.
Popiel I, Erasmus DA: Schistosoma mansoni: ultrastructure of adults from mice treated with oxamniquine. Exp Parasitol. 1984, 58 (3): 254-262. 10.1016/0014-4894(84)90042-0.
Popiel I, Cioli D, Erasmus DA: Reversibility of drug-induced regression in female Schistosoma mansoni upon pairing with male worms in vivo. Z Parasitenkd. 1984, 70 (3): 421-424. 10.1007/BF00927830.
Kunz W: Schistosome male-female interaction: induction of germ-cell differentiation. Trends Parasitol. 2001, 17 (5): 227-231. 10.1016/S1471-4922(01)01893-1.
Popiel I, Basch PF: Reproductive development of female Schistosoma mansoni (Digenea: Schistosomatidae) following bisexual pairing of worms and worm segments. J Exp Zool. 1984, 232 (1): 141-150. 10.1002/jez.1402320117.
Shaw JR: Schistosoma mansoni: pairing in vitro and development of females from single sex infections. Exp Parasitol. 1977, 41 (1): 54-65. 10.1016/0014-4894(77)90129-1.
Clough ER: Morphology and reproductive organs and oogenesis in bisexual and unisexual transplants of mature Schistosoma mansoni females. J Parasitol. 1981, 67 (4): 535-539. 10.2307/3280487.
Armstrong JC: Mating behavior and development of schistosomes in the mouse. J Parasitol. 1965, 51: 605-616. 10.2307/3276242.
Michaels RM: Mating of Schistosoma mansoni in vitro. Exp Parasitol. 1969, 25 (1): 58-71. 10.1016/0014-4894(69)90052-6.
Shaw JR, Marshall I, Erasmus DA: Schistosoma mansoni: in vitro stimulation of vitelline cell development by extracts of male worms. Exp Parasitol. 1977, 42 (1): 14-20. 10.1016/0014-4894(77)90056-X.
Basch PF, Basch N: Intergeneric reproductive stimulation and parthenogenesis in Schistosoma mansoni. Parasitology. 1984, 89 ( Pt 2): 369-376.
Cornford EA, Fitzpatrick AM: Comparative glucose utilization rates in separated and mated schistosomes. Experimental Parasitology. 1987, 64: 448-457. 10.1016/0014-4894(87)90059-2.
Basch PF: Why do schistosomes have separate sexes?. Parasitol Today. 1990, 6 (5): 160-163. 10.1016/0169-4758(90)90339-6.
Gupta BC, Basch PF: The role of Schistosoma mansoni males in feeding and development of female worms. J Parasitol. 1987, 73 (3): 481-486. 10.2307/3282125.
Conford EM, Huot ME: Glucose transfer from male to female schistosomes. Science. 1981, 213 (4513): 1269-1271. 10.1126/science.7268436.
Silveira AM, Friche AA, Rumjanek FD: Transfer of [14C] cholesterol and its metabolites between adult male and female worms of Schistosoma mansoni. Comp Biochem Physiol B. 1986, 85 (4): 851-857. 10.1016/0305-0491(86)90186-0.
Haseeb MA: Schistosoma mansoni: females enhance [C-14]-tyrosine incorporation in males maintained in vitro. Journal of Helminthology. 1998, 72 (2): 123-126.
Khalil SB, Mansour NS: Worm development in hamsters infected with unisex and cross-mated Schistosoma mansoni and Schistosoma haematobium. J Parasitol. 1995, 81 (1): 8-11. 10.2307/3283998.
Zanotti EM, Magalhaes LA, Piedrabuena AE: [Morphology and development of Schistosoma mansoni Sambon, 1907 in unisexual infections experimentally induced in mice]. Rev Saude Publica. 1982, 16 (2): 114-119.
Basch PF, Gupta BC: Homosexual male pairing in Schistosoma mansoni. Int J Parasitol. 1988, 18 (8): 1115-1117. 10.1016/0020-7519(88)90084-7.
Siegel DA, Tracy JW: Effect of pairing in vitro on the glutathione level of male Schistosoma mansoni. J Parasitol. 1988, 74 (4): 524-531. 10.2307/3282167.
Siegel DA, Tracy JW: Schistosoma mansoni: influence of the female parasite on glutathione biosynthesis in the male. Exp Parasitol. 1989, 69 (2): 116-124. 10.1016/0014-4894(89)90179-3.
Haseeb MA, Fried B, Eveland LK: Schistosoma mansoni: female-dependent lipid secretion in males and corresponding changes in lipase activity. Int J Parasitol. 1989, 19 (7): 705-709. 10.1016/0020-7519(89)90054-4.
Bostic JR, Strand M: Molecular cloning of a Schistosoma mansoni protein expressed in the gynecophoral canal of male worms. Mol Biochem Parasitol. 1996, 79 (1): 79-89. 10.1016/0166-6851(96)02640-0.
Hoffmann KF, Johnston DA, Dunne DW: Identification of Schistosoma mansoni gender-associated gene transcripts by cDNA microarray profiling. Genome Biol. 2002, 3 (8): RESEARCH0041-10.1186/gb-2002-3-8-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 (2): 191-209. 10.1016/j.molbiopara.2004.03.014.
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): 1-8. 10.1016/j.ijpara.2005.10.008.
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, 147 (1): 39-47.
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 (3): 160-72. 10.1016/j.exppara.2006.03.003.
Fitzpatrick JM, Hoffmann KF: Dioecius Schistosoma mansoni express divergent gene repertoires regulated by pairing. Int J Parasitol. 2006, 36: 1081-1089. 10.1016/j.ijpara.2006.06.007.
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 (7-8): 919-929. 10.1007/s00018-005-5578-1.
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, 20 (5): 280-289.
Hoffmann KF, McCarty TC, Segal DH, Chiaramonte M, Hesse M, Davis EM, Cheever AW, Meltzer PS, Morse HC, Wynn TA: Disease fingerprinting with cDNA microarrays reveals distinct gene expression profiles in lethal type 1 and type 2 cytokine-mediated inflammatory reactions. FASEB J. 2001, 15 (13): 2545-2547.
Smyth GK: Linear models and empirical Bayes methods for assessing differential expression in microarray experiments. Stat Appl Genet Mol Biol. 2004, 3: Article3-
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): 1-13. 10.1016/j.molbiopara.2005.01.007.
Reis MG, Kuhns J, Blanton R, Davis AH: Localization and pattern of expression of a female specific mRNA in Schistosoma mansoni. Mol Biochem Parasitol. 1989, 32 (2-3): 113-119. 10.1016/0166-6851(89)90062-5.
Chen LL, Rekosh DM, LoVerde PT: Schistosoma mansoni p48 eggshell protein gene: characterization, developmentally regulated expression and comparison to the p14 eggshell protein gene. Mol Biochem Parasitol. 1992, 52 (1): 39-52. 10.1016/0166-6851(92)90034-H.
Aronstein WS, Strand M: A glycoprotein antigen of Schistosoma mansoni expressed on the gynecophoral canal of mature male worms. Am J Trop Med Hyg. 1985, 34 (3): 508-512.
Schussler P, Grevelding CG, Kunz W: Identification of Ras, MAP kinases, and a GAP protein in Schistosoma mansoni by immunoblotting and their putative involvement in male-female interaction. Parasitology. 1997, 115 ( Pt 6): 629-634.
Lawrence JD: The ingestion of red blood cells by Schistosoma mansoni. J Parasitol. 1973, 59 (1): 60-63. 10.2307/3278572.
Hoffmann KF: An historical and genomic view of schistosome conjugal biology with emphasis on sex-specific gene expression. Parasitology. 2004, 128 Suppl: S11-22.
Erasmus DA, Davies TW: Schistosoma mansoni and S. haematobium: calcium metabolism of the vitelline cell. Exp Parasitol. 1979, 47 (1): 91-106. 10.1016/0014-4894(79)90011-0.
Gotz B, Klinkert MQ: Expression and partial characterization of a cathepsin B-like enzyme (Sm31) and a proposed 'haemoglobinase' (Sm32) from Schistosoma mansoni. Biochem J. 1993, 290 ( Pt 3): 801-806.
Kennedy MW, Scott JC, Lo S, Beauchamp J, McManus DP: Sj-FABPc fatty-acid-binding protein of the human blood fluke Schistosoma japonicum: structural and functional characterization and unusual solvent exposure of a portal-proximal tryptophan residue. Biochem J. 2000, 349 (Pt 1): 377-384. 10.1042/0264-6021:3490377.
Gobert GN, Stenzel DJ, Jones MK, McManus DP: Immunolocalization of the fatty acid-binding protein Sj-FABPc within adult Schistosoma japonicum. Parasitology. 1997, 115 ( Pt 1): 33-39. 10.1017/S0031182097008925.
Ramos TM, de Vasconcelos AS, de Carvalho VC, Lima VL: [Alterations in cholesterol, triglyceride and total phospholipid levels in plasma of Callithrix jacchus (sagui) reinfected by Schistosoma mansoni]. Rev Soc Bras Med Trop. 2004, 37 (1): 37-40.
Chechinel YMN, Owen JS: O colesterol plasmático na esquistossomose mansônica. Revista Brasileira de Biologia. 1978, 38: 339-343.
Doenhoff MJ, Stanley RG, Griffiths K, Jackson CL: An anti-atherogenic effect of Schistosoma mansoni infections in mice associated with a parasite-induced lowering of blood total cholesterol. Parasitology. 2002, 125 (Pt 5): 415-421. 10.1017/S0031182002002275.
Brouwers JF, Smeenk IM, van Golde LM, Tielens AG: The incorporation, modification and turnover of fatty acids in adult Schistosoma mansoni. Mol Biochem Parasitol. 1997, 88 (1-2): 175-185. 10.1016/S0166-6851(97)00091-1.
Meyer F, Meyer H, Bueding E: Lipid metabolism in the parasitic and free-living flatworms, Schistosoma mansoni and Dugesia dorotocephala. Biochim Biophys Acta. 1970, 210 (2): 257-266.
Schuler GD: Sequence mapping by electronic PCR. Genome Res. 1997, 7 (5): 541-550.
Altschul SF, Gish W, Miller W, Myers EW, Lipman DJ: Basic local alignment search tool. J Mol Biol. 1990, 215 (3): 403-410.
Steel RGG, Torrie JH, DIckey DA: Principles and Procedures of Statistics. A Biometric Approach. McGraw-Hill Series in Probability and Statistics. Edited by: Rogers M. 1997, WCB McGraw-Hill, Third Edition
Erasmus DA: A comparative study of the reproductive system of mature, immature and "unisexual" female Schistosoma mansoni. Parasitology. 1973, 67 (2): 165-183.
Bogea T, Favre TC, Rotenberg L, Silva HS, Pieri OS: Circadian pattern of cercarial emergence in Schistosoma mansoni (Platyhelminthes:Digenea) from isolated Biomphalaria glabrata. Chronobiol Int. 1996, 13 (2): 93-101.
Smithers SR, Terry RJ: The infection of laboratory hosts with cercariae of Schistosoma mansoni and the recovery of the adult worms. Parasitology. 1965, 55 (4): 695-700.
Sambrook J, Fritsch EF, Maniatis T: Molecular Cloning, A Laboratory Manual. 1989, Cold Spring Harbor Laboratory Press, Second Edition
Hasseman J: Aminoallyl labelling of RNA for microarrays - SOP #: M004. [http://pga.tigr.org/sop/M004_1a.pdf]REVISION LEVEL: 2
Livak KJ, Schmittgen TD: Analysis of relative gene expression data using real-time quantitative PCR and the 2(-Delta Delta C(T)) Method. Methods. 2001, 25 (4): 402-408. 10.1006/meth.2001.1262.
Brazma A, Hingamp P, Quackenbush J, Sherlock G, Spellman P, Stoeckert C, Aach J, Ansorge W, Ball CA, Causton HC, Gaasterland T, Glenisson P, Holstege FC, Kim IF, Markowitz V, Matese JC, Parkinson H, Robinson A, Sarkans U, Schulze-Kremer S, Stewart J, Taylor R, Vilo J, Vingron M: Minimum information about a microarray experiment (MIAME)-toward standards for microarray data. Nat Genet. 2001, 29 (4): 365-371. 10.1038/ng1201-365.
Parkinson H, Sarkans U, Shojatalab M, Abeygunawardena N, Contrino S, Coulson R, Farne A, Lara GG, Holloway E, Kapushesky M, Lilja P, Mukherjee G, Oezcimen A, Rayner T, Rocca-Serra P, Sharma A, Sansone S, Brazma A: ArrayExpress--a public repository for microarray gene expression data at the EBI. Nucleic Acids Res. 2005, 33 (Database issue): D553-5. 10.1093/nar/gki056.
Gentleman R, Ihaka R: The R language. Edited by: Billard L, Fisher N. 1997, The Interface Foundation of North America
Gentleman RC, Carey VJ, Bates DM, Bolstad B, Dettling M, Dudoit S, Ellis B, Gautier L, Ge Y, Gentry J, Hornik K, Hothorn T, Huber W, Iacus S, Irizarry R, Leisch F, Li C, Maechler M, Rossini AJ, Sawitzki G, Smith C, Smyth G, Tierney L, Yang JY, Zhang J: Bioconductor: open software development for computational biology and bioinformatics. Genome Biol. 2004, 5 (10): R80-10.1186/gb-2004-5-10-r80.
Smyth GK, Speed TP: Normalization of cDNA microarray data. Methods. 2003, 31: 265-273. 10.1016/S1046-2023(03)00155-5.
Smyth GK, Michaud J, Scott H: The use of within-array replicate spots for assessing differential expression in microarray experiments. Bioinformatics. 2005, 21 (9): 2067-2075. 10.1093/bioinformatics/bti270.
GenScript C: Real-Time PCR Primer Design. [https://www.genscript.com/ssl-bin/app/primer]
Conesa A, Gotz S, Garcia-Gomez JM, Terol J, Talon M, Robles M: Blast2GO: a universal tool for annotation, visualization and analysis in functional genomics research. Bioinformatics. 2005, 21 (18): 3674-3676. 10.1093/bioinformatics/bti610.
IVIA CG: Blast2Go Tutorial. [http://blast2go.bioinfo.cipf.es/b2gdata/blast2goTutorial.pdf]Version 1.0
The authors would like to acknowledge CAPES, FAPEMIG and NIH for financial support. Work at UFMG was supported by FAPEMIG grant CBB 469/04 to GRF. Work in the El-Sayed laboratory was supported by a NIH grant AI48828 to NES. MW's work at the El-Sayed laboratory was supported by a PDEE scholarship from CAPES and his work at UFMG was supported by a CAPES PhD scholarship. We would like to thank Prof. Karl Hoffman for providing the arrays, Prof. Elida ML Rabelo for thoughtful suggestions, Bryan Frank for technical support and advice concerning the microarray experiments and Delza de Moura Soares Reis and Sueleny Silva Ferreira Texeira for technical assistance with S. mansoni life-cycle maintenance and mouse perfusions.
MW led the design of the study, performed the RNA amplifications, hybridizations, real time PCRs, statistical analysis of the data and drafted the manuscript. LKJP carried out the unisexual and mixed infections of B. glabrata and injection of mice. FPL carried out the mouse perfusions, separation of worms and RNA extractions. GCC participated in the data analysis performing the mePCR, real time PCR oligonucleotides design and subsequent data treatment. GF, NES and OSC participated in the design and coordination of the study and helped to draft the manuscript. All authors read and approved the final manuscript
Electronic supplementary material
About this article
Cite this article
Waisberg, M., Lobo, F.P., Cerqueira, G.C. et al. Microarray analysis of gene expression induced by sexual contact in Schistosoma mansoni. BMC Genomics 8, 181 (2007). https://doi.org/10.1186/1471-2164-8-181