De novo assembly and transcriptome characterization: novel insights into the natural resistance mechanisms of Microtus fortis against Schistosoma japonicum
© Hu et al.; licensee BioMed Central Ltd. 2014
Received: 21 September 2013
Accepted: 19 May 2014
Published: 2 June 2014
Microtus fortis is a non-permissive host of Schistosoma japonicum. It has natural resistance against schistosomes, although the precise resistance mechanisms remain unclear. The paucity of genetic information for M. fortis limits the use of available immunological methods. Thus, studies based on high-throughput sequencing technologies are required to obtain information about resistance mechanisms against S. japonicum.
Using Illumina single-end technology, a de novo assembly of the M. fortis transcriptome produced 67,751 unigenes with an average length of 868 nucleotides. Comparisons were made between M. fortis before and after infection with S. japonicum using RNA-seq quantification analysis. The highest number of differentially expressed genes (DEGs) occurred two weeks after infection, and the highest number of down-regulated DEGs occurred three weeks after infection. Simultaneously, the strongest pathological changes in the liver were observed at week two. Gene ontology terms and pathways related to the DEGs revealed that up-regulated transcripts were involved in metabolism, immunity and inflammatory responses. Quantitative real-time PCR analysis showed that patterns of gene expression were consistent with RNA-seq results.
After infection with S. japonicum, a defensive reaction in M. fortis commenced rapidly, increasing dramatically in the second week, and gradually decreasing three weeks after infection. The obtained M. fortis transcriptome and DEGs profile data demonstrated that natural and adaptive immune responses, play an important role in M. fortis immunity to S. japonicum. These findings provide a better understanding of the natural resistance mechanisms of M. fortis against schistosomes.
Schistosome infections occur in over 70 tropical and subtropical countries, with 200 million individuals infected, resulting in 200,000 human deaths annually . The pathology of chronic infections with S. japonicum includes granuloma formation in the liver, periportal fibrosis, portal hypertension and the formation of vascular shunts . Serious schistosomiasis reduces the ability of humans to work and impedes the development of agriculture and livestock. It is crucial to control schistosomiasis transmission, improve the level of public health knowledge and increase income in agriculture and livestock. At present, praziquantel (PZQ) is the only safe, cheap and effective drug against human schistosome species. However, massive administration of PZQ in endemic areas has raised serious concerns regarding the development of parasite resistance to PZQ . Isolates of S. mansoni from patients in Kisumu, Kenya, were shown to have lower susceptibility to PZQ . Clarifying the immune mechanisms of S. japonicum infection is an important foundation to screen effective vaccine candidates and novelty drugs against schistosomiasis.
The reed vole (Microtus fortis) has been identified as a non-permissive host of S. japonicum that has a natural resistance mechanism against schistosomes [5–8]. Several studies of the resistance of M. fortis to S. japonicum have been conducted. Eosinophils, neutrophils, macrophages and serum antibodies are thought to be involved in resistance  and nitric oxide, albumin, E77, KPNA2, HSP90 were found to be resistance-associated proteins [9–11]. It has been suggested that apoptosis might be employed by M. fortis for eliminating schistosomes . This study aimed to provide new suggestions to control the spread of schistosomiasis and to further elucidate the resistance mechanisms against schistosomes in M. fortis.
Because M. fortis is a new experimental animal model, biological information about it is very limited. To date, only 328 nucleotides, 68 ESTs and 26 genes of M. fortis are accessible from all NCBI databases, and no monoclonal antibodies specific for these markers are commercially available. It is therefore extremely difficult to study resistance mechanisms using immunological techniques, such as flow cytometry, Western blot or microarrays. Due to the high cost of genome sequencing, transcriptome sequencing is a useful tool to obtain biological information and to discover and identify new genes of M. fortis. The next-generation sequencing (NGS) based RNA-Seq technique has been widely used for de novo transcriptome sequencing assembly, discovery of novel genes and the investigation of gene expression in many non-model organisms . It is therefore suitable for the identification of potential resistance-related genes in M. fortis.
In this study, we sequenced the M. fortis transcriptome using Illumina technology and demonstrated the suitability of short-read sequencing for the de novo assembly and annotation of genes expressed in a eukaryote without prior genome information. In addition, we compared the gene expression profiles of M. fortis up to four weeks after infection with schistosomes by using an RNA-Seq quantification expression system. The assembled and annotated transcriptome sequences and gene expression profiles showed that natural and adaptive immune responses, such as Toll-like receptors (TLR), natural killer (NK) cell activation and IgG and IgE, play important roles against S. japonicum infections. These data provide a valuable resource for the analysis of resistance genes of M. fortis against S. japonicum.
Pathological changes in the liver
Illumina sequencing and de novo transcriptome assembly
To obtain an overview of an M. fortis gene expression profile, cDNA from the liver of M. fortis was generated and sequenced using the Illumina sequencing platform. After removing adaptors and low quality reads, 134,985,444 clean reads were obtained from sequencing. The GC content and Q20 were 48.98% and 98.03%, respectively. These reads were assembled with the Trinity software. After clustering using the TGICL program, we obtained 166,501 contigs and 67,751 unigenes (mean length: 868 nt) with an N50 of 1,631 nt (Additional file 1: Table S1). Simultaneously, an abundance of unannotated unigenes was generated in our study, which might represent a specific gene pool of interest for further M. fortis studies.
Function annotation and classification of assembled transcripts
Pathways of de novo transcriptome involved in immunity and inflammation responses
All genes with pathway annotation (23898)
MAPK signaling pathway
Chemokine signaling pathway
Wnt signaling pathway
Fc gamma R-mediated phagocytosis
Cytokine-cytokine receptor interaction
Leukocyte transendothelial migration
JAK-STAT signaling pathway
NF-kappa B signaling pathway
T cell receptor signaling pathway
B cell receptor signaling pathway
Complement and coagulation cascades
NK cell-mediated cytotoxicity
Notch signaling pathway
TGF-beta signaling pathway
TLR signaling pathway
NOD-like receptor signaling pathway
Antigen-processing and presentation
Intestinal immune network for IgA production
DEGs involved in the response to infection with S. japonicum
DGE sequencing statistics
Map to Gene
Total mapped Reads
≤ 2 bp mismatch
Total unmapped reads
M. fortis 0w
6.91E + 08
M. fortis 1w
7.13E + 08
M. fortis 2w
7.06E + 08
M. fortis 3w
6.87E + 08
M. fortis 4w
6.91E + 08
GO and pathway functional enrichment analysis of DEGs
To understand the functions of the DEGs, we performed GO functional enrichment and KEGG pathway analyses using a hypergeometric distribution model. When comparing M. fortis samples before infection and four weeks after infection, the terms “cell periphery”, “extracellular space”, “plasma membrane”, “external side of plasma membrane”, “extracellular region”, “cell surface” and “cytoplasm” showed the most significant differences (Additional files 4, 5, 6 and 7: Datasets S2, S3, S4, S5).
By analyzing the statistically enriched GO function related to the DEGs, it was revealed that the most significant and most frequently enriched terms for up-regulated transcripts after infection were inflammatory- and immunity-correlated genes. The first class included pattern recognition receptors (PRRs), such as mannose-binding lectin (MBL, K03991), C-reactive protein (CRP, K01672), scavenger receptor (SR, GO:0005044), complement receptor (CR, GO:0002430) and TLRs (GO:0002755, GO:0006954, GO:0034162, GO:0034134). The second class consisted of cytokines and complement-related unigenes, such as interferon-gamma (IFN-γ)-related unigenes (GO:0003924, GO:0000975, GO:0005515, GO:0003723), tumor necrosis factor (TNF)-related unigenes (GO:0005164, GO:0033209, GO:0005031, GO:0043120) and complement-related unigenes (GO:0006956, GO:0006957, GO:0006958, GO:0001851). The third class was comprised of lymphocyte activation-related molecules, such as T- and B lymphocyte-related unigenes (GO:0003823, GO:0015643, GO:0005488, GO:0019901, GO:0005068) and macrophage-related unigenes (GO:0008603, GO:0043621, GO:0070891).
Pathway of RNA-seq involved in immunity and inflammation responses after infection for four weeks
DEGs with pathway annotation
All genes with pathway annotation
1st w (700)
2nd w (946)
3rd w (665)
4th w (432)
Toll-like receptor signaling
Natural killer cell mediated cytotoxicity
Fc gamma R-mediated phagocytosis
Complement and coagulation cascades
B cell receptor signaling pathway
T cell receptor signaling pathway
Cytokine-cytokine receptor interaction
RIG-I-like receptor signaling pathway
NOD-like receptor signaling pathway
TGF-beta signaling pathway
Wnt signaling pathway
Ubiquitin mediated proteolysis
Validation of Illumina expression patterns by qRT-PCR analysis
Potential resistance-associated molecular analysis in M. fortis
Metabolism-related unigenes, such as glutathione S-transferase kappa 1, GTPase family M protein 1, 3 beta-hydroxysteroid dehydrogenase, calcium-dependent phospholipase A2, mitochondrial precursor, apolipoprotein L, insulin-like growth factor-binding protein 1, serpin B9-like, cytochrome P450 4A10;
Immunity-related unigenes, such as immunoglobulin IgG heavy chain, h-2 class I histocompatibility antigen, CD74 antigen, MHC class II antigen, interferon-induced protein with tetratricopeptide repeats 1-like, interferon regulatory factor 1, interferon regulatory factor 7, interferon regulatory factor 8, interferon regulatory factor 9, lymphocyte antigen 86-like, T cell surface glycoprotein CD8 alpha chain-like, B cell linker protein, monocyte differentiation antigen CD14-like;
Inflammation-related unigenes, such as fibronectin precursor, C-X-C motif chemokine 10-like, C-C motif chemokine 5, Janus kinase 3;
Apoptosis-related genes, such as caspase-8 and bcl-2 homologous antagonist.
The results showed that unigenes down-regulated after infection were, among others, salivary gland secretion protein 4, cholesterol 7-alpha-monooxygenase, UDP-glucuronosyltransferase 2B1-like isoform 2, UDP-glucuronosyltransferase 1-9-like, hematopoietically expressed homeobox, zinc finger and BTB domain-containing protein 16. These results showed that after infection of M. fortis with S. japonicum, metabolism, immunity, inflammation, and apoptosis-related gene expression levels increased, whereas synthesis was reduced.
M. fortis and BALB/c mice are non-permissive and permissive hosts of S. japonicum, respectively. After infection, the profiles of gene expression, immune response and pathological changes of the two hosts were significantly different. A previous study revealed that the lungs and liver of M. fortis, rats and mice display different characteristics when infected with S. japonicum for 10 days. M. fortis and rats, both of which are resistant to S. japonicum, developed a stronger immune response and more severe pathological lesions in response to schistosomes than mice during the early phase of the infection . After infection with S. japonicum, levels of cytokines, complements and antibodies were significantly higher in M. fortis than in mice [9, 19]. In our study, the changes in cytokine expression were detected by using protein microarrays in the sera of M. fortis voles and C57BL/6 mice. Levels of Th1, Th2, and Th17 cytokines in M. fortis voles increased significantly during the first 3 weeks post-infection, while there was no significant changes in mice (Data not shown). It has been found that after infection with S. japonicum for 1 week, serine protease inhibitor gene expression was up-regulated and immune-associated gene expression did not change in the lung tissue of mice. Conversely, some important immune-associated genes (such as CD74, MHC-I and MHC-II) were up-regulated, and serine protease inhibitor gene expression did not change in M. fortis . Pathological changes in the liver of M. fortis were significant during the first 3 weeks post-infection, while there was no significant changes in the liver of mice (shown in Figure 1). These results help to clarify molecular mechanisms of infection immunity in different hosts of S. japonicum, which might provide new ways to prevent and control schistosomiasis.
To investigate the changes in gene expression and understand the critical genes in the response of M. fortis to S. japonicum infection, RNA-Seq (quantification) analysis was performed and various differently-expressed genes were obtained. The number of up-regulated genes at the second week after infection was 1,354, which was higher than at other time points. Compared to the second week after infection, the number of down-regulated genes at the third and the fourth week after infection was 1,202 and 1,304, respectively, which was higher than at other weeks. It was thus concluded that the in vivo response of M. fortis was the most intense at the second week after infection. Three weeks after infection, the response to schistosomes in M. fortis gradually leveled off. This result was consistent with pathological changes in the liver after infection. GO and pathway functional enrichment analysis of the DEGs showed that “cell periphery”, “extracellular space”, “plasma membrane”, “external side of plasma membrane”, “extracellular region”, “cell surface” and “cytoplasm” were the most significantly changed categories. These structures include metabolic enzymes, cell surface receptors, signaling molecules, and platelet-derived growth factor-binding structures. The most significantly and most frequently enriched terms that were up-regulated transcripts after infection were inflammatory, immunity-correlated genes, and metabolism-related enzymes (Additional file 8: Dataset S6).
NK cells are a type of cytotoxic lymphocytes critical to the innate immune system. They control immune responses by secreting cytokines and chemokines, including tumor necrosis factor-α (TNF-α) and interferon-γ (IFN-γ), and eliminate target cells by polarized exocytosis of cytotoxic granules containing perforin, granzymes and Fas ligand . Binding of IFN-γ (or type II interferon) to its receptors activates the JAK-STAT pathway and has an immune-modulatory function . IFN-γ (originally called macrophage-activating factor) can stimulate macrophages and induce direct antimicrobial and antitumor mechanisms. It can also up-regulate antigen-processing and antigen-presenting pathways . IFN-γ orchestrates leukocyte attraction, thus enhancing NK cell activity and regulating B cell functions, such as immunoglobulin (Ig) production and class switching . Our results show that unigenes involved in IFN-γ, such as interferon regulatory factor 1, interferon regulatory factor 7, interferon regulatory factor 8, interferon regulatory factor 9 and interferon-induced protein with tetratricopeptide repeats 1-like, were significantly up-regulated. Simultaneously, by using genome oligonucleotide microarrays, some studies also reported that after infection with S. japonicum for 10 days, a gene similar to IFN-γ receptor 2 was up-regulated in the lungs of M. fortis . Members of the TNF/TNF-receptor (TNF-R) superfamily coordinate the immune response at multiple levels, such as by lymphoid neogenesis, augmenting immune responses and apoptosis dependent on caspase 8. Among these, acute TNF-R engagement can induce the rapid and robust activation of NF-κB and MAPK pathways . In our study, the TNF-related genes and pathways were significantly up-regulated. These included TNF-receptor binding (ligand) superfamily, TNF-mediated signaling pathway, TNF-activated receptor activity member 1B precursor, TNF binding, caspase 8 (K04398), apoptosis, NF-κB and MAPK pathways. This indicates that IFN-γ and TNF might play important roles in the response against S. japonicum infections.
TLRs are among the best-characterized PRRs and are key regulators of anti-pathogen immune responses. Schistosomal-derived lysophosphatidylcholine participates in the production of cytokines, such as TNF-α and IL-10, and in eosinophil activation through a TLR2-dependent pathway in S. mansoni and S. haematobium infection [26, 27]. Antigen-presenting cells (APCs) recognize lysophosphatidylserine and schistosomal glycolipids of schistosomes through TLRs, leading to NF-κB and MAPK signaling pathway activation and inducing inflammatory cytokine production .
The CD14 antigen is a glycosylphosphatidylinositol (GPI)-anchored receptor known to serve as a co-receptor for several TLRs that were significantly up-regulated in M. fortis after infection with S. japonicum. Earlier studies showed a role for CD14 as a co-receptor working with TLR4 and facilitating cellular responses to low doses of lipopolysaccharide (LPS). Upon LPS stimulation, CD14 enables the activation of the TLR4-TRAM-TRIF pathway and leads to the activation of NFAT transcription factor family members . This suggests that CD14 and TLRs play an important role in initiating inflammation by responding to LPS stimulation in M. fortis.
Adaptive immunity is an important part of the anti-schistosomiasis mechanism in M. fortis. Three weeks after infection, many unigenes were involved in adaptive immunity, including immunoglobulin IgG heavy chain, Fc receptor of IgE and CD8A antigen protein complex, CD74. Antibodies are major components of the immune system. IgG is the main antibody isotype found in blood and extracellular fluid, allowing it to control the infection of body tissues. IgG can mediate a variety of biological functions, such as the classical pathway of complement activation, ADCC, hypersensitivity, and can block pathogens. The Fc receptors of IgE are proteins called Fc-epsilon receptors (FcϵR) found on the surface of certain cells - including mast cells, eosinophils, basophils and Langerhans cells - that contribute to the protective functions of the immune system . The CD8 antigen is a cell surface glycoprotein found on most cytotoxic T lymphocytes that mediates efficient cell-cell interactions within the immune system. CD74, known as the major histocompatibility complex (MHC) class II invariant chain, regulates the trafficking of MHC-II in APCs and acts as a receptor for macrophage migration inhibitory factor (MIF). After the binding of the MIF, NF-κB and Erk1/2 activation occurs, along with the induction of proinflammatory cytokine secretion . It has been reported that the serum of M. fortis was passively transferred in mice, which not only reduced worm burden but also led to worm shortening . Another report showed that purified IgG3 antibodies from laboratory-bred M. fortis and wild M. fortis effectively killed schistosomula, and that IgG3 antibodies of wild M. fortis induced a higher worm-reduction rate. The death rates of schistosomula due to IgG3 antibody purified from the sera of laboratory-bred M. fortis and wild M. fortis were 2.35 and 5.88 times higher than in Km mice, respectively . It has been shown that the sera and/or spleen cells of M. fortis had an important effect on killing schistosomula in vitro. This study could also show a synergism between sera of M. fortis and spleen cells . In our study, after infection with S. japonicum for three weeks, cytokines (such as IL-4, IL-5, IL-6, IL-13, TNF-α and GM-CSF) in the sera of M. fortis were significantly increased, whereas cytokine levels of the BALB/c mice were much lower (data not shown). These findings suggest that these parasite-killing effects were performed in part by CD8+ T and B lymphocytes.
In summary, we first characterized the M. fortis de novo transcriptome and performed RNA-Seq quantification analysis after infection with S. japonicum. A total of 134,985,444 clean reads, 166,501 contigs and 67,751 unigenes were obtained. Based on the assembled de novo transcriptome, the amount of DEGs that play significant roles in the response to schistosome infection were identified. The numbers of DEGs at weeks 1–4 were 1741, 2178, 1683 and 1197, respectively. We concluded that the in vivo response of M. fortis was the most intense at the second week after infection. GO and pathway enrichment analysis revealed that after infection the up-regulated DEGs were predominately involved in metabolism, inflammation response and immunity response. The natural and adaptive immune responses, particularly NK cell activation producing IFN-γ and TGF-β, TLR activation of macrophages, and IgG and IgE antibodies produced by B lymphocytes, play an important role in M. fortis resistance against S. japonicum. Multiple candidate genes involved in innate immune responses were identified further with qRT-PCR, and expression patterns in M. fortis and BALB/c mice were compared. The results reflected significant differences of innate immune responses and different sensitivities in the hosts of schistosomes. This study provides a valuable molecular basis for the analysis of M. fortis resistance mechanisms against S. japonicum.
This study was performed in strict accordance with the recommendations of the Laboratory of Animal Welfare & Ethics Committee (LAWEC) of China. No specific permissions from Hunan provincial authorities were required because the study did not involve endangered or protected species. The protocol was approved by the LAWEC Committee of the National Institute of Parasitic Diseases, Chinese Center for Disease Control and Prevention (approval ID: IPD 2009–4). All surgery was performed under sodium pentobarbital anesthesia, and all efforts were made to minimize suffering.
Animals and cercariae
M. fortis voles were live captured in the Hunan Dongtin Lakes region, China. The animals were transported to Shanghai and reared in independent ventilated cages. Fifteen M. fortis females (aged 7–8 weeks, weighing 60–70 g) were randomly selected from a sub-generation of wild animals to be used in the experiment. BALB/c mice (aged 7–8 weeks, weighing 20–25 g) were purchased from the Experimental Animal Department of the Chinese Academy of Sciences. The animals were infected cutaneously (see below) with a mainland strain of S. japonicum cercariae. Infected Oncomelanias hupensis snails were provided by the National Institute of Parasitic Diseases, Chinese Center for Disease Control and Prevention. The infected snails were placed in dechlorinated water under artificial light to induce cercarial shedding before animal infections.
M. fortis voles and BALB/c mice were infected percutaneously with S. japonicum cercariae. Based on the susceptibility of each rodent species to S. japonicum, M. fortis voles were infected with 1,000 cercariae and BALB/c mice were infected with 40 cercariae per individual. Fifteen M. fortis voles and 15 BALB/c mice were divided into five groups randomly. Animals in each group were sacrificed either before infection or at 1, 2, 3 and 4 weeks post-infection. Before the animals were sacrificed, they were anesthetized with 0.75% sodium pentobarbital by intraperitoneal injection, and the livers of the animals were removed, frozen in liquid nitrogen and stored at -80°C until further use.
Preparation of cDNA library and Illumina sequencing for transcriptome analysis
Total RNA was extracted using Trizol reagent (Invitrogen, USA) following the manufacturer’s protocol. After total RNA extraction and DNase I treatment, magnetic beads with Oligo (dT) were used to isolate mRNA. The mRNA was fragmented into short fragments when it was mixed with the fragmentation buffer. Then, cDNA was synthesized using the mRNA fragments as templates. The short fragments were purified and resolved with EB buffer for end reparation and single nucleotide A (adenine) addition. Afterwards, the short fragments were connected with adapters. The suitable fragments were selected as templates for PCR amplification. During the QC steps, an Agilent 2100 Bioanalyzer and an ABI StepOnePlus Real-Time PCR System were used for quantification and qualification of the sample library. The library was sequenced using Illumina HiSeq™ 2000.
Raw reads were produced from the sequences. Clean reads were obtained after removal of low quality reads, reads with adaptors and reads with unknown nucleotides larger than 5%. A de novo transcriptome assembly was performed using the short reads assembling program Trinity . The resulting sequences obtained from Trinity are called unigenes. Then, a BLASTx alignment (E-value < 0.00001) between unigenes and protein databases, such as NR, Swiss-Prot, KEGG and COG, was performed, and the best alignment results were used to determine the sequence direction of the unigenes. With NR annotation, we used the Blast2GO program  to obtain the GO annotation of the unigenes. After obtaining the GO annotation for every unigene, we used WEGO software  to perform a GO functional classification for all unigenes and to further understand the distribution of gene functions in the species on a macroscopic level. We predicted genes with different expression levels, and performed GO functional analysis and KEGG pathway analysis on them.
RNA-Seq quantification analysis
Four independent cDNA libraries were constructed for the four liver samples according to the RNA-Seq protocol. Raw image files were collected using the Illumina HiSeqTM2000 sequencing platform in BGI Shenzhen, China (http://en.genomics.cn/navigation/index.action). The analyzed data have been deposited in the NCBI Sequence Read Archive under the Accession No. SRX337491.
Screening of DEGs
Screening of DEGs included the screening of genes that are differentially expressed among samples followed by a GO functional enrichment analysis and a KEGG pathway enrichment analysis for these DEGs. The expression level of each gene was measured via RNA-Seq quantification analysis using the RPKM  method (reads per kb per million reads). Referring to “The significance of digital gene expression profiles” , we developed a strict algorithm (the Poisson distribution) to identify DEGs between two samples. We performed a cluster analysis of gene expression patterns with cluster  software and Java Treeview  software.
GO offers a dynamically-updated controlled vocabulary and a strictly defined concept to comprehensively describe properties of genes and their products in any organism. GO enrichment analysis provides all GO terms that are significantly enriched in DEGs compared to the genome background and filters those DEGs that correspond to biological functions. This method first maps all DEGs to GO terms in a database (http://www.geneontology.org/). Genes usually interact with one another to play roles in certain biological functions. The pathway enrichment analysis identified significantly enriched metabolic pathways or signal transduction pathways in DEGs compared with the entire genome background.
RNA isolation for reverse transcription and real-time quantitative PCR (qPCR)
Total RNA was extracted from the livers of M. fortis and BALB/c mice using a QIAgen RNeasy Mini Kit according to the manufacturer’s specifications. The RNA was determined using a NanoDrop 2000 spectrophotometer (Thermo Scientific, USA), and the integrity was evaluated using electrophoresis on an agarose gel stained with ethidium bromide. Reactions were performed in a GeneAmp® PCR System 9700 (Applied Biosystems, USA) for 15 min at 37°C, followed by heat inactivation of RT for 5 s at 85°C. Real-time PCR was performed using a LightCycler® 480 II Real-time PCR Instrument (Roche, Switzerland) with 10 μl of PCR reaction mixture. Reactions were incubated in a 384-well optical plate (Roche, Switzerland) at 95°C for 10 min, followed by 40 cycles of 95°C for 10 s and 60°C for 30 s. Each sample was analyzed in triplicates. At the end of the PCR cycles, a melting curve analysis was performed to validate the specific generation of the expected PCR product. The primer sequences (Additional file 9: Table S3) were designed in the laboratory and synthesized by Generay Biotech (Generay, China) based on the mRNA sequences obtained from the NCBI database and the RNA-Seq. The expression levels of mRNAs were normalized to GAPDH and were calculated using the 2-ΔΔCt method .
Availability of supporting data
Raw sequence read data were submitted to the NCBI Sequence Read Archive under the Accession No. SRX337491, http://www.ncbi.nlm.nih.gov/bioproject/PRJNA215691).
This work was funded by the National Natrral Science Foundation (No. 30801047 to Y.H., No 81371841 to J.C.), and by the National S&T Major Program (No. 2012ZX10004-201 to J.C). We thank the staff of the Beijing Genomics institution for their help in the bioinformatic analysis. We thank Master Ling Wang for kind help with graph drawing.
- Gryseels B, Polman K, Clerinx J, Kestens L: Human schistosomiasis. Lancet. 2006, 368: 1106-1118. 10.1016/S0140-6736(06)69440-3.PubMedView ArticleGoogle Scholar
- Cai P, Piao X, Liu S, Hou N, Wang H, Chen Q: MicroRNA-gene expression network in murine liver during Schistosoma japonicum infection. PLoS One. 2013, 8: e67037-10.1371/journal.pone.0067037.PubMed CentralPubMedView ArticleGoogle Scholar
- Doenhoff MJ, Pica-Mattocia L: Praziquantel for the treatment of schistosomiasis: its use for control in areas with endemic disease and prospects for drug resistance. Exp Rev Anti Infect Ther. 2006, 4: 199-210. 10.1586/14787126.96.36.199.View ArticleGoogle Scholar
- Melman SD, Steinauer ML, Cunningham C, Kubatko LS, Mwangi IN, Wynn NB, Mutuku MW, Karanja DM, Colley DG, Black CL, Secor WE, Mkoji GM, Loker ES: Reduced susceptibility to praziquantel among naturally occurring Kenyan isolates of Schistosoma mansoni. PLoS Negl Trop Dis. 2009, 3: e504-10.1371/journal.pntd.0000504.PubMed CentralPubMedView ArticleGoogle Scholar
- Lu BL: Identification of Important Rodentia Of China. Identification Manual of Important Medical Animals of China. Edited by: Lu BL. 1982, Beijing: Publishing House for People’s Health, 907-916. 1Google Scholar
- He YX, Salafsky B, Ramaswamy K: Host-parasite relationships of Schistosoma japonicum in mammalian hosts. Trends Parasitol. 2001, 17: 320-324. 10.1016/S1471-4922(01)01904-3.PubMedView ArticleGoogle Scholar
- Li FK, Zhu ZL, Jin BR: Uninfectibility to Schistosoma japonicum of Microtus fortis. Acta Parasitol Sinica. 1965, 2: 103-105.Google Scholar
- Li H, He YY, Lin JJ: The observation for the phenomenon of Microtus Fortis against Schistosoma japonicum. Chin J Vet Parasitol. 2000, 8: 12-15.Google Scholar
- Hu Y, Cao JP: Research progress on schistosomiasis resistance in the non-permissive murine host. Int J Med Parasit Dis. 2010, 37: 94-97.Google Scholar
- Cheng G, Gong Q, Gai N, Xiong DH, Yu YJ, Zeng QR, Hu WX: Karyopherin alpha 2 (KPNA2) is associated with the natural resistance to Schistosoma japanicum infection in Microtus fortis. Biomed Pharmacother. 2011, 65: 230-237. 10.1016/j.biopha.2011.02.012.PubMedView ArticleGoogle Scholar
- Gong Q, Cheng G, Qin ZQ, Xiong DH, Yu YJ, Zeng QR, Hu WX: Identification of the resistance of a novel molecule heat shock protein 90alpha (HSP90alpha) in Microtus fortis to Schistosoma japonicum infection. Acta Trop. 2010, 115: 220-226. 10.1016/j.actatropica.2010.03.007.PubMedView ArticleGoogle Scholar
- Peng J, Gobert GN, Hong Y, Jiang W, Han H, McManus DP, Wang X, Liu J, Fu Z, Shi Y: Apoptosis governs the elimination of Schistosoma japonicum from the non-permissive host Microtus fortis. PLoS One. 2011, 6: e21109-10.1371/journal.pone.0021109.PubMed CentralPubMedView ArticleGoogle Scholar
- Jung KH, Ko HJ, Mguyen MX, Kim SR, Ronald P, An G: Genomewide identification and analysis of early heat stress responsive genes in rice. J Plant Biol. 2012, 55: 458-468. 10.1007/s12374-012-0271-z.View ArticleGoogle Scholar
- Shao GY, Feng J, Bai X, Zhao Y, Liu DH, Gao J, Xie JY, Gao C, Hu JH: Pathology and gene expression in Microtus fortis infected with Schistosoma japonicum. Acta labor Anim Sci Sin. 2013, 21: 56-60.Google Scholar
- Jiang W, Hong Y, Peng J, Fu Z, Feng X, Liu J, Shi Y, Lin J: Study on differences in the pathology, T cell subsets and gene expression in susceptible and non-susceptible hosts infected with Schistosoma japonicum. PLoS One. 2010, 5: e13494-10.1371/journal.pone.0013494.PubMed CentralPubMedView ArticleGoogle Scholar
- Yao LX, Lin JJ, Cai YM: Advances in research of Microtus fortis against schistosomiasis. Chin J Schisto Control. 2003, 15: 471-473.Google Scholar
- Roman LT, Michael YG, Darren AN, Eugene VK: The COG database: a tool for genome-sc ale analysis of protein functions and evolution. Nucleic Acids Res. 2000, 28: 33-36. 10.1093/nar/28.1.33.View ArticleGoogle Scholar
- Li R, Yu C, Li Y, Lam TW, Yiu SM, Kristiansen K, Wang J: SOAP2: An improved ultrafast tool for short read alignment. Bioinformatics. 2009, 25: 1966-1967. 10.1093/bioinformatics/btp336.PubMedView ArticleGoogle Scholar
- He Y, Luo X, Yu X, Lin J, Liu S, Zhang X, Li Y: Preliminary study on the levels of natural antibodies against Schistosoma japonicum in Microtus fortis in Dongting lake area. Zhongguo Ji Sheng Chong Xue Yu Ji Sheng Chong Bing Za Zhi. 1999, 17: 132-134.PubMedGoogle Scholar
- Sun J, Lin JJ, Cheng GF, Li H, Shi YJ, Lu K, Cai YM, Jiang HB: Gene difference expression between mice and Microtus fortis infected with Schistosoma japonicum using cDNA microarrays. Acta Sci Nat Univ Pekin. 2004, 40: 532-537.Google Scholar
- Kwon HJ, Kim HS: Signaling for synergistic activation of natural killer cells. Immune Netw. 2012, 12: 240-246. 10.4110/in.2012.12.6.240.PubMed CentralPubMedView ArticleGoogle Scholar
- Schoenborn JR, Wilson CB: Regulation of interferon-gamma during innate and adaptive immune responses. Adv Immunol. 2007, 96: 41-101.PubMedView ArticleGoogle Scholar
- Hu X, Ivashkiv LB: Cross-regulation of signaling pathways by interferon-gamma: implications for immune responses and autoimmune diseases. Immunity. 2009, 31: 539-550. 10.1016/j.immuni.2009.09.002.PubMed CentralPubMedView ArticleGoogle Scholar
- Schroder K, Hertzog PJ, Ravasi T, Hume DA: Interferon-gamma: an overview of signals, mechanisms and functions. J Leukoc Biol. 2004, 75: 163-189.PubMedView ArticleGoogle Scholar
- Clark J, Vagenas P, Panesar M, Cope AP: What does tumor necrosis factor excess do to the immune system long term?. Ann Rheum Dis. 2005, 64 Suppl 4: iv70-iv76.PubMedGoogle Scholar
- Meurs L, Labuda L, Amoah AS, Mbow M, Ngoa UA, Boakye DA, Mboup S, Dièye TN, Mountford AP, Turner JD, Kremsner PG, Polman K, Yazdanbakhsh M, Adegnika AA: Enhanced pro-inflammatory cytokine responses following Toll-like-receptor ligation in Schistosoma haematobium-infected schoolchildren from rural Gabon. PLoS One. 2011, 6: e24393-10.1371/journal.pone.0024393.PubMed CentralPubMedView ArticleGoogle Scholar
- Magalhães K, Almeida PE, Atella G, Maya-Monteiro CM, Castro-Faria-Neto H, Pelajo-Machado M, Lenzi HL, Bozza MT, Bozza PT: Schistosomal-derived lysophosphatidylcholine are involved in eosinophil activation and recruitment through Toll-like receptor-2-dependent mechanisms. J Infect Dis. 2010, 202: 1369-1379. 10.1086/656477.PubMedView ArticleGoogle Scholar
- Venugopal PG, Nutman TB, Semnani RT: Activation and regulation of toll-like receptors (TLRs) by helminth parasites. Immunol Res. 2009, 43: 252-263. 10.1007/s12026-008-8079-0.PubMed CentralPubMedView ArticleGoogle Scholar
- Zanoni I, Granucci F: Role of CD14 in host protection against infections and in metabolism regulation. Front Cell Infect Microbiol. 2013, 3: 32-PubMed CentralPubMedGoogle Scholar
- Prussin C, Metcalfe D: 5. IgE, mast cells, basophils, and eosinophils. J Allergy Clin Immunol. 2006, 117: S450-S456. 10.1016/j.jaci.2005.11.016.PubMedView ArticleGoogle Scholar
- Beswick EJ, Reyes VE: CD74 in antigen presentation, inflammation, and cancers of the gastrointestinal tract. World J Gastroenterol. 2009, 15: 2855-2861. 10.3748/wjg.15.2855.PubMed CentralPubMedView ArticleGoogle Scholar
- Luo XS, He YK, Yu XL, Lin JL, Liu SX, Li Y, Zhang XY: A study on the protective immunity against Schistosoma japonicum infection after passively transferring Microtus fortis sera to mice. Chin J Zoonoses. 1998, 14: 75-76.Google Scholar
- Jiang SF, Wei MX, Lin JJ, Pan CE, Qiu AW, He YY, Li H, Shi YJ: Effect of IgG3 Antibody Purified from Sera of Microtus fortis Against Schistosoma. Chin J Parasitol Parasit Dis. 2008, 26: 34-36.Google Scholar
- Zhang XY, He YK, Li Y, Luo XS, Liu SX, Yu XL, Lin JL, Li YS: A preliminary observation on the killing of schistosomula of Schistosoma japonicum in vitro by sera and spleen cells of Microtus fortis. Chin J Schisto Control. 2001, 13: 206-208.Google Scholar
- Grabherr MG, Haas BJ, Yassour M, Levin JZ, Thompson DA, Amit I, Adiconis X, Fan L, Raychowdhury R, Zeng Q, Chen Z, Mauceli E, Hacohen N, Gnirke A, Rhind N, di Palma F, Birren BW, Nusbaum C, Lindblad-Toh K, Friedman N, Regev A: Full-length transcriptome assembly from RNA-Seq data without a reference genome. Nat Biotechnol. 2011, 29: 644-652. 10.1038/nbt.1883.PubMed CentralPubMedView ArticleGoogle Scholar
- Conesa A, Götz S, García-Gómez JM, Terol J, Talón M, Robles M: Blast2GO: a universal tool for annotation, visualization and analysis in functional genomics research. Bioinformatics. 2005, 21: 3674-3686. 10.1093/bioinformatics/bti610.PubMedView ArticleGoogle Scholar
- Ye J, Fang L, Zheng H, Zhang Y, Chen J, Zhang Z, Wang J, Li S, Li R, Bolund L, Wang J: WEGO: a web tool for plotting GO annotations. Nucleic Acids Res. 2006, 34 (Web Server issue): W293-W297.PubMed CentralPubMedView ArticleGoogle Scholar
- Mortazavi A, Williams BA, McCue K, Schaeffer L, Wold B: Mapping and quantifying mammalian transcriptomes by RNA-Seq. Nat Methods. 2008, 5: 621-628. 10.1038/nmeth.1226.PubMedView ArticleGoogle Scholar
- Audic S, Claverie JM: The significance of digital gene expression profiles. Genome Res. 1997, 7: 986-995.PubMedGoogle Scholar
- de Hoon MJ, Imoto S, Nolan J, Miyano S: Open source clustering software. Bioinformatics. 2004, 20: 1453-1454. 10.1093/bioinformatics/bth078.PubMedView ArticleGoogle Scholar
- Saldanha AJ: Java Treeview-extensible visualization of microarray data. Bioinformatics. 2004, 20: 3246-3258. 10.1093/bioinformatics/bth349.PubMedView ArticleGoogle Scholar
- Livak KJ, Schmittgen TD: Analysis of relative gene expression data using real-time quantitative PCR and the 22DDCT method. Methods. 2001, 25: 402-408. 10.1006/meth.2001.1262.PubMedView ArticleGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly credited. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.