Immune tolerance of vector beetle to its partner plant parasitic nematode modulated by its insect parasitic nematode.

Immune response of insect vectors to transmitted pathogens or insect hosts against parasites are well studied, whereas the mechanism of tripartite interactions remains elusive. In this study, we investigated the immune interactions of the vector beetle Monochamus alternatus ( Ma) to the devastating plant parasitic nematode Bursaphelenchus xylophilus ( Bx) and the insect parasitic nematode Howardula phyllotretae ( Hp). We report the unique immune mechanism by which the vector beetle tolerates many devastating Bx in its trachea, yet that immune tolerance is compromised by the parasitic nematode Hp. Contact with either nematode species triggers epithelial reactive oxygen species (ROS) production in Ma. Only the entry of Bx, not Hp, infection, induces increased expression of antioxidative genes, through which the ROS levels are balanced in the trachea of beetles. Furthermore, we found that up-regulation of antioxidative genes was induced by the interaction of Toll receptors. In contrast, beetles infected by Hp retain high levels of oxidative stress and melanization in trachea, and as a result, decrease Bx loading. This study highlights the role of Toll receptors in mediating the activation of antioxidative genes in immune tolerance to plant parasitic nematodes, and suggests the use of insect parasites as a biologic control.-Zhou, J., Zhao, L.-L., Yu, H.-Y., Wang, Y.-H., Zhang, W., Hu, S.-N., Zou, Z., Sun, J.-H. Immune tolerance of vector beetle to its partner plant parasitic nematode modulated by its insect parasitic nematode.

[Advances of sequencing and assembling technologies for complex genomes].

Complex genomes are noted to be extremely difficult to sequence or assemble by using ordinary methods. Complex genomes are typically characterized as being highly repetitive, highly heterozygous, extremely GC biased, or naturally contaminated, i.e., contaminations which cannot be removed before sequencing. To solve these problems with sequencing and assembling complex genomes, three major techniques include: (1) DNA extraction experiments, (2) Sequencing technologies and platforms, and (3) Algorithms and strategies for assembling. In this review, we summarize these state-of-the-art technologies and strategies used in these directions. We also review the representative projects of complex genome sequencing and address the development of these technologies and strategies for solving the challenges when sequencing or assembling complex genomes.

A Critical Role for CLSP2 in the Modulation of Antifungal Immune Response in Mosquitoes.

Entomopathogenic fungi represent a promising class of bio-insecticides for mosquito control. Thus, detailed knowledge of the molecular mechanisms governing anti-fungal immune response in mosquitoes is essential. In this study, we show that CLSP2 is a modulator of immune responses during anti-fungal infection in the mosquito Aedes aegypti. With a fungal infection, the expression of the CLSP2 gene is elevated. CLSP2 is cleaved upon challenge with Beauveria bassiana conidia, and the liberated CLSP2 CTL-type domain binds to fungal cell components and B. bassiana conidia. Furthermore, CLPS2 RNA interference silencing significantly increases the resistance to the fungal challenge. RNA-sequencing transcriptome analysis showed that the majority of immune genes were highly upregulated in the CLSP2-depleted mosquitoes infected with the fungus. The up-regulated immune gene cohorts belong to melanization and Toll pathways, but not to the IMD or JAK-STAT. A thioester-containing protein (TEP22), a member of α2-macroglobulin family, has been implicated in the CLSP2-modulated mosquito antifungal defense. Our study has contributed to a greater understanding of immune-modulating mechanisms in mosquitoes.

De Novo Sequencing and Transcriptome Analysis of Pleurotus eryngii subsp. tuoliensis (Bailinggu) Mycelia in Response to Cold Stimulation.

Cold stimulation of Bailinggu's mycelia is the main factor that triggers primordia initiation for successful production of fruiting bodies under commercial cultivation. Yet, the molecular-level mechanisms involved in mycelia response to cold stimulation are still unclear. Here, we performed comparative transcriptomic analysis using RNA-Seq technology to better understand the gene expression regulation during different temporal stages of cold stimulation in Bailinggu. A total of 21,558 Bailinggu mycelia unigenes were de novo assembled and annotated from four libraries (control at 25 °C, plus cold stimulation treatments at -3 °C for a duration of 1-2 days, 5-6 days, and 9-10 days). GO and KEGG pathway analysis indicated that functional groups of differentially expressed unigenes associated with cell wall and membrane stabilization, calcium signaling and mitogen-activated protein kinases (MAPK) pathways, and soluble sugars and protein biosynthesis and metabolism pathways play a vital role in Bailinggu's response to cold stimulation. Six hundred and seven potential EST-based SSRs loci were identified in these unigenes, and 100 EST-SSR primers were randomly selected for validation. The overall polymorphism rate was 92% by using 10 wild strains of Bailinggu. Therefore, these results can serve as a valuable resource for a better understanding of the molecular mechanisms associated with Bailinggu's response to cold stimulation.

Genome-wide mapping of 5-hydroxymethylcytosine in three rice cultivars reveals its preferential localization in transcriptionally silent transposable element genes.

5-Hydroxymethylcytosine (5hmC), a modified form of cytosine that is considered the sixth nucleobase in DNA, has been detected in mammals and is believed to play an important role in gene regulation. In this study, 5hmC modification was detected in rice by employing a dot-blot assay, and its levels was further quantified in DNA from different rice tissues using liquid chromatography-multistage mass spectrometry (LC-MS/MS/MS). The results showed large intertissue variation in 5hmC levels. The genome-wide profiles of 5hmC modification in three different rice cultivars were also obtained using a sensitive chemical labelling followed by a next-generation sequencing method. Thousands of 5hmC peaks were identified, and a comparison of the distributions of 5hmC among different rice cultivars revealed the specificity and conservation of 5hmC modification. The identified 5hmC peaks were significantly enriched in heterochromatin regions, and mainly located in transposable elements (TEs), especially around retrotransposons. The correlation analysis of 5hmC and gene expression data revealed a close association between 5hmC and silent TEs. These findings provide a resource for plant DNA 5hmC epigenetic studies and expand our knowledge of 5hmC modification.

Insights into the genomic evolution and the alkali tolerance mechanisms of Agaricus sinodeliciosus by comparative genomic and transcriptomic analyses.

Agaricus sinodeliciosus is a rare wild edible mushroom from northwest China, and grows naturally in mild saline-alkali soil, which is also unusual in mushrooms. A. sinodeliciosus represents a potential model organism for explaining saline-alkali tolerance mechanisms and revealing related physiological processes in mushrooms. Here, we provide a high-quality genome of A. sinodeliciosus. Comparative genomic analyses reveal A. sinodeliciosus has numerous changes to its genome organization after a solitary evolutionary history under saline-alkali environments, such as gene family contraction, retrotransposon expansion and rapid evolution of adaptative genes. Our saline and alkali tolerance tests show that mycelium growth and fruit body formation of this species are effected by mild alkalinity. Transcriptomic analyses reveal that genes involved in carbon and nitrogen utilization, cell stability and fruit body formation of A. sinodeliciosus could be activated under mildly alkaline conditions. In particular, the 'starch and sucrose metabolism', 'biosynthesis of amino acids' and 'phenylpropanoid biosynthesis' pathways are important for mildly alkaline tolerance of A. sinodeliciosus. Like plants and arbuscular mycorrhizal fungi, in the rot fungus A. sinodeliciosus, the biosynthesis of intracellular small molecules could be enhanced to counter osmotic and oxidative stresses caused by mild alkalinity, and the biosynthesis of monolignol could be suppressed to increase cell wall infiltrates under mildly alkaline conditions. This research provides an understanding of the genomic evolution and mechanisms of A. sinodeliciosus in tolerance to saline-alkali environments. The A. sinodeliciosus genome constitutes a valuable resource for evolutionary and ecological studies of Agaricus.

Comparative transcriptome analysis of transporters, phytohormone and lipid metabolism pathways in response to arsenic stress in rice (Oryza sativa).

• Arsenic (As) contamination of rice (Oryza sativa) is a worldwide concern and elucidating the molecular mechanisms of As accumulation in rice may provide promising solutions to the problem. Previous studies using microarray techniques to investigate transcriptional regulation of plant responses to As stress have identified numerous differentially expressed genes. However, little is known about the metabolic and regulatory network remodelings, or their interactions with microRNA (miRNA) in plants upon As(III) exposure. • We used Illumina sequencing to acquire global transcriptome alterations and miRNA regulation in rice under As(III) treatments of varying lengths of time and dosages. • We found that the response of roots was more distinct when the dosage was varied, whereas that of shoots was more distinct when the treatment time was varied. In particular, the genes involved in heavy metal transportation, jasmonate (JA) biosynthesis and signaling, and lipid metabolism were closely related to responses of rice under As(III) stress. Furthermore, we discovered 36 new As(III)-responsive miRNAs, 14 of which were likely involved in regulating gene expression in transportation, signaling, and metabolism. • Our findings highlight the significance of JA signaling and lipid metabolism in response to As(III) stress and their regulation by miRNA, which provides a foundation for subsequent functional research.

Identification of microRNA in the protist Trichomonas vaginalis.

MicroRNAs (miRNAs) are a class of small noncoding RNAs that have important regulatory roles in multicellular organisms. However, miRNA has never been identified experimentally in protist. Direct cloning of 438 expressed miRNA tags by microRNA serial analysis of gene expression from the parasitic protist Trichomonas vaginalis identified nine candidate miRNAs. Bioinformatics analysis of the corresponding genomic region revealed that these miRNA candidates contain a classical stem-loop-stem structure of pre-microRNAs. Analysis of the 20 nt long mature tva-miR-001 showed that it is an intergenic miRNA located at the scaffold DS113596. Tva-miR-001 was differentially expressed in the trophozoite, pseudocyst and amoeboid stages. Based on the experimental results of the present study, we provided solid evidence that protist possesses a miRNA regulating network comparable with multicellular organisms for the first time.

Genetic factors related to the widespread dissemination of ST11 extensively drug-resistant carbapenemase-producing Klebsiella pneumoniae strains within hospital.

BACKGROUND: Carbapenemase-producing Klebsiella pneumoniae (CP-Kp) poses distinct clinical challenges due to extensively drug resistant (XDR) phenotype, and sequence type (ST) 11 is the most dominant blaKPC-2-bearing CP-Kp clone in China. The purpose of this current retrospective study was to explore the genetic factors associated with the success of XDR CP-Kp ST11 strains circulated in the intensive care unit (ICU) of a Chinese tertiary hospital. METHODS: Six ST11 XDR CP-Kp strains were identified between May and December 2014 and validated by minimum inhibitory concentration examination, polymerase chain reaction, and pyrosequencing. The six ST11 XDR CP-Kp, as well as three multi-drug resistant (MDR) and four susceptible strains, were sequenced using single-molecule real-time method. Comprehensively structural and functional analysis based on comparative genomics was performed to identify genomic characteristics of the XDR ST11 CP-Kp strains. RESULTS: We found that ST11 XDR blaKPC-2-bearing CP-Kp strains isolated from inpatients spread in the ICU of the hospital. Functionally, genes associated with information storage and processing of the ST11 XDR CP-Kp strains were more abundant than those of MDR and susceptible strains, especially genes correlative with mobile genetic elements (MGEs) such as transposons and prophages. Structurally, eleven large-scale genetic regions taken for the unique genome in these ST11 XDR CP-Kp strains were identified as MGEs including transposons, integrons, prophages, genomic islands, and integrative and conjugative elements. Three of them were located on plasmids and eight on chromosomes; five of them were with antimicrobial resistance genes and eight with adaptation associated genes. Notably, a new blaKPC-2-bearing ΔΔTn1721-blaKPC-2 transposon, probably transposed and truncated from ΔTn1721-blaKPC-2 by IS903D and ISKpn8, was identified in all six ST11 XDR CP-Kp strains. CONCLUSION: Our findings suggested that together with clonal spread, MGEs identified uniquely in the ST11 XDR CP-Kp strains might contribute to their formidable adaptability, which facilitated their widespread dissemination in hospital.

Analysis of expressed sequence tags in porcine uterus tissue.

Two non-normalized cDNA libraries of uteri from Danish Landrace and Chinese Erhualian pigs were constructed, and 13,756 expressed sequence tags (ESTs) were randomly sequenced. The ESTs were clustered by Phrap software, and 6,139 distinct tentative consensus sequences were produced, including 2,730 contigs and 3,409 singlets. Using Blast tools, these 6,139 candidate genes were compared to the nr and nt databases; 5,210 of them were assigned putative functions, whereas 929 potentially represent new genes. Highly expressed genes appear to be associated with basic energy metabolism, transferase activity, localization, cellular physiological process, protein binding, and nucleic acid binding. Antileukoproteinase was the most highly expressed gene, corresponding to endometrial differentiation and conceptus or fetal development.

Comparative analysis of ESTs in response to drought stress in chickpea (C. arietinum L.).

Chickpea (Cicer arietinum L.) is an important pulse crop grown mainly in the arid and semi-arid regions. To identify the water-stress-induced genes, two non-normalized cDNA libraries were constructed from the seedling leaves of a drought-tolerant chickpea cultivar under PEG-treated and -nontreated conditions. About 2500 clones from each library were selected randomly for sequencing analysis. Based on IDEG6 online software analysis, 92 genes were differentially expressed, and these genes were involved in diverse biological progresses, such as metabolism, transcription, signal transduction, protein synthesis and others. Most of the up-regulated genes were related to drought tolerance, and the down-regulated genes were mainly involved in photosynthesis. The differential expression patterns of five functional unigenes were confirmed by quantitative real-time PCR (qPCR). The results will help in understanding the molecular basis of drought tolerance in chickpea.

[Analysis of expressed sequence tags (EST) from Panax quinquefolium root].

To investigate the profile of gene expression in American ginseng (Panax quinquefolium L.) and discover its functional genes, for the first time, expressed sequence tags (EST) library of four-year-old American ginseng roots has been established. According to BLAST and Gene Ontology analysis, eleven genes, encoding cytochrome P450, glucosyltransferase, farnesyltransferase and cyclase family protein, are found to be associated with ginsenosides biosynthesis. Six other genes are obtained encoding auxin-regulated protein, auxin response factor 4 and auxin-repressed protein in the roots of American ginseng. In addition, thirteen expressed transcripts are stress-connected proteins and twelve expressed other transcripts are closely related to plant defense in four-year-old American ginseng roots. Furthermore, 62 genes no hit in BLAST and in Interproscan may be new genes. These results indicate EST is an useful tool for research on functional genomics of P. quinquefolium and it can be applied to the molecular modification of the ginsenosides biosynthetic pathway ultimately for improving the quality of American ginseng germplasm.

Comparative analysis of the Monochamus alternatus immune system.

The pine sawyer beetle, Monochamus alternatus, is regarded as a notorious forest pest in Asia, vectoring an invasive pathogenic nematode, Bursaphelenchus xylophilus, which is known to cause pine wilt disease. However, little sequence information is available for this vector beetle. This hampered the research on its immune system. Based on the transcriptome of M. alternatus, we have identified and characterized 194 immunity-related genes in M. alternatus, and compared them with homologues molecules from other species known to exhibit immune responses against invading microbes. The lower number of putative immunity-related genes in M. alternatus were attributed to fewer C-type lectin, serine protease (SP) and anti-microbial peptide (AMP) genes. Phylogenetic analysis revealed that M. alternatus had a unique recognition gene, galectin3, orthologues of which were not identified in Tribolium castaneum, Drosophila melanogastor, Anopheles gambiae and Apis mellifera. This suggested a lineage-specific gene evolution for coleopteran insects. Our study provides the comprehensive sequence resources of the immunity-related genes of M. alternatus, presenting valuable information for better understanding of the molecular mechanism of innate immunity processes in M. alternatus against B. xylophilus.

Oral Microbiota Community Dynamics Associated With Oral Squamous Cell Carcinoma Staging.

Oral squamous cell carcinoma (OSCC) is a highly aggressive cancer and the fourth leading malignancy among males in Taiwan. Some pathogenic bacteria are associated with periodontitis and oral cancer. However, the comprehensive profile of the oral microbiome during the cancer's progression from the early stage to the late stage is still unclear. We profiled the oral microbiota and identified bacteria biomarkers associated with OSCC. The microbiota of an oral rinse from 51 healthy individuals and 197 OSCC patients at different stages were investigated using 16S rRNA V3V4 amplicon sequencing, followed by bioinformatics and statistical analyses. The oral microbiota communities from stage 4 patients showed significantly higher complexity than those from healthy controls. The populations also dynamically changed with the cancer's progression from stage 1 to stage 4. The predominant phyla in the oral samples showed variation in the relative abundance of Fusobacteria, Bacteroidetes, and Actinobacteria. The abundance of Fusobacteria increased significantly with the progression of oral cancer from the healthy controls (2.98%) to OSCC stage 1 (4.35%) through stage 4 (7.92%). At the genus level, the abundance of Fusobacterium increased, while the number of Streptococcus, Haemophilus, Porphyromonas, and Actinomyces decreased with cancer progression. Fusobacterium periodonticum, Parvimonas micra, Streptococcus constellatus, Haemophilus influenza, and Filifactor alocis were associated with OSCC, and they progressively increased in abundance from stage 1 to stage 4. The abundances of Streptococcus mitis, Haemophilus parainfluenzae, and Porphyromonas pasteri were inversely associated with OSCC progression. We selected a bacterial marker panel of three bacteria (upregulated F. periodonticum, down-regulated S. mitis, and P. pasteri), which had an AUC of 0.956 (95% CI = 0.925-0.986) in discriminating OSCC stage 4 from the healthy controls. Furthermore, the functional prediction of oral bacterial communities showed that genes involved in carbohydrate-related metabolism, such as methane metabolism, and energy-metabolism-related parameters, such as oxidative phosphorylation and carbon fixation in photosynthetic organisms, were enriched in late-stage OSCC, while those responsible for amino acid metabolism, such as folate biosynthesis and valine, leucine, and isoleucine biosynthesis, were significantly associated with the healthy controls. In conclusion, our results provided evidence of oral bacteria community changes during oral cancer progression and suggested the possibility of using bacteria as OSCC diagnostic markers.

The calcium-dependent protein kinase (CDPK) and CDPK-related kinase gene families in Hevea brasiliensis-comparison with five other plant species in structure, evolution, and expression.

Calcium-dependent protein kinases (CDPKs or CPKs) play important roles in various physiological processes of plants, including growth and development, stress responses and hormone signaling. Although the CDPK gene family has been characterized in several model plants, little is known about this gene family in Hevea brasiliensis (the Para rubber tree). Here, we characterize the entire H. brasiliensis CDPK and CDPK-related kinase (CRK) gene families comprising 30 CDPK genes (HbCPK1 to 30) and nine CRK genes (HbCRK1 to 9). Structure and phylogeny analyses of these CDPK and CRK genes demonstrate evolutionary conservation in these gene families across H. brasiliensis and other plant species. The expression of HbCPK and HbCRK genes was investigated via Solexa sequencing in a range of experimental conditions (different tissues, phases of leaf development, ethylene treatment, and various abiotic stresses). The results suggest that HbCPK and HbCRK genes are important components in growth, development, and stress responses of H. brasiliensis. Parallel studies on the CDPK and CRK gene families were also extended to five other plant species (Arabidopsis thaliana, Oryza sativa, Populus trichocarpa, Manihot esculenta, and Ricinus communis). The CDPK and CRK genes from different plant species that exhibit similar expression patterns tend to cluster together, suggesting a coevolution of gene structure and expression behavior in higher plants. The results serve as a foundation to further functional studies of these gene families in H. brasiliensis as well as in the whole plant kingdom.

Analysis of 10,000 ESTs from lymphocytes of the cynomolgus monkey to improve our understanding of its immune system.

BACKGROUND: The cynomolgus monkey (Macaca fascicularis) is one of the most widely used surrogate animal models for an increasing number of human diseases and vaccines, especially immune-system-related ones. Towards a better understanding of the gene expression background upon its immunogenetics, we constructed a cDNA library from Epstein-Barr virus (EBV)-transformed B lymphocytes of a cynomolgus monkey and sequenced 10,000 randomly picked clones. RESULTS: After processing, 8,312 high-quality expressed sequence tags (ESTs) were generated and assembled into 3,728 unigenes. Annotations of these uniquely expressed transcripts demonstrated that out of the 2,524 open reading frame (ORF) positive unigenes (mitochondrial and ribosomal sequences were not included), 98.8% shared significant similarities (E-value less than 1e-10) with the NCBI nucleotide (nt) database, while only 67.7% (E-value less than 1e-5) did so with the NCBI non-redundant protein (nr) database. Further analysis revealed that 90.0% of the unigenes that shared no similarities to the nr database could be assigned to human chromosomes, in which 75 did not match significantly to any cynomolgus monkey and human ESTs. The mapping regions to known human genes on the human genome were described in detail. The protein family and domain analysis revealed that the first, second and fourth of the most abundantly expressed protein families were all assigned to immunoglobulin and major histocompatibility complex (MHC)-related proteins. The expression profiles of these genes were compared with that of homologous genes in human blood, lymph nodes and a RAMOS cell line, which demonstrated expression changes after transformation with EBV. The degree of sequence similarity of the MHC class I and II genes to the human reference sequences was evaluated. The results indicated that class I molecules showed weak amino acid identities (<90%), while class II showed slightly higher ones. CONCLUSION: These results indicated that the genes expressed in the cynomolgus monkey could be used to identify novel protein-coding genes and revise those incomplete or incorrect annotations in the human genome by comparative methods, since the old world monkeys and humans share high similarities at the molecular level, especially within coding regions. The identification of multiple genes involved in the immune response, their sequence variations to the human homologues, and their responses to EBV infection could provide useful information to improve our understanding of the cynomolgus monkey immune system.

[HLA-A, B, DRB1 gene polymorphism of Beijing population was studied by high-resolution polymerase chain reaction sequence-based typing].

OBJECTIVE: To investigate the gene polymorphism of human leukocyte antigen (HLA)-A, B, DRB1 loci in the population of Beijing region, and research on the application feasibility of polymerase chain reaction sequence-based typing (PCR-SBT) method. METHODS: PCR-SBT method was applied to determine HLA- A, B, DRB1 genotypes of 618 unrelated healthy individuals of Beijing region. RESULTS: A total of 84 different alleles and 199 genotypes of HLA-A, 143 alleles and 366 genotypes of HLA-B, 122 alleles and 286 genotypes of HLA-DRB1 were detected. CONCLUSION: The results showed the characteristics of HLA-A, B, DRB1 distributions, and provided more comprehensive and accurate gene data that may serve as normal reference values for all of Beijing people.

[Frequencies distribution of human leukocyte antigen-B27 subtypes in healthy Chinese].

OBJECTIVE: To investigate the frequencies distribution of human leukocyte antigen (HLA)-B27 subtypes in unrelated healthy Chinese. METHOD: Polymerase chain reaction sequence-based typing (PCRSBT) was used to determine HLA high-resolution genotypes of 825 unrelated healthy Chinese. RESULTS: A total of 25 HLA-B27-positive individuals and 8 HLA-B27 subtypes were detected. These subtypes and their corresponding frequencies were B * 2704 (30.77%) , B * 2705 (23.08%), B * 2707 (19.23%), B * 2711 (7.69%), B * 2712 (7.69%), B * 2701 (3.85%), B * 2713 (3.85%) and B * 2721 (3.85%). CONCLUSION: The data obtained through PCR-SBT method may serve as important reference for the research of relationship between HLA-B27 subtypes and some diseases such as ankylosing spondylitis.

[Exclusion of the association of five known mutations with congenital stationary nyctalopia in a large Chinese family].

OBJECTIVE: To detect gene mutations associated with autosomal dominant congenital stationary night blindness(ADCSNB) in a large Chinese family. METHODS: Genomic DNAs were extracted from peripheral blood samples of 16 affected and 14 unaffected family members. According to 5 missense mutations in 3 genes reported previously, 4 pairs of primers were designed and corresponding exons containing the five mutation sites were amplified by polymerase chain reaction. Amplified products were purified and sequenced by MegaBACE1000 capillary array electrophoresis DNA sequencer. Full field electroretinogram (ERG, ISCEV) of patients was recorded and analyzed by Roland Consult System. RESULTS: Dark-adapted ERG showed a-wave was normal, but b-wave of the patients was markedly decreased. None of the five missense mutations were detected in 16 affected and 14 unaffected family members. CONCLUSION: The molecular pathogenesis of ADCSNB in this family does not involve point mutations or deletions of these five sites, which indicates the heterogeneity of ADCSNB.

Transcriptomic insight into the immune defenses in the ghost moth, Hepialus xiaojinensis, during an Ophiocordyceps sinensis fungal infection.

Hepialus xiaojinensis is an economically important species of Lepidopteran insect. The fungus Ophiocordyceps sinensis can infect its larvae, which leads to mummification after 5-12 months, providing a valuable system with which to study interactions between the insect hosts and pathogenic fungi. However, little sequence information is available for this insect. A time-course analysis of the fat body transcriptome was performed to explore the host immune response to O. sinensis infection. In total, 50,164 unigenes were obtained by assembling the reads from two high-throughput approaches: 454 pyrosequencing and Illumina Hiseq2000. Hierarchical clustering and functional examination revealed four major gene clusters. Clusters 1-3 included transcripts markedly induced by the fungal infection within 72 h. Cluster 4, with a lower number of transcripts, was suppressed during the early phase of infection but returned to normal expression levels sometime before 1 year. Based on sequence similarity to orthologs known to participate in immune defenses, 258 candidate immunity-related transcripts were identified, and their functions were hypothesized. The genes were more primitive than those in other Lepidopteran insects. In addition, lineage-specific family expansion of the clip-domain serine proteases and C-type lectins were apparent and likely caused by selection pressures. Global expression profiles of immunity-related genes indicated that H. xiaojinensis was capable of a rapid response to an O. sinensis challenge; however, the larvae developed tolerance to the fungus after prolonged infection, probably due to immune suppression. Specifically, antimicrobial peptide mRNAs could not be detected after chronic infection, because key components of the Toll pathway (MyD88, Pelle and Cactus) were downregulated. Taken together, this study provides insights into the defense system of H. xiaojinensis, and a basis for understanding the molecular aspects of the interaction between the host and the entomopathogen.