Transcriptomic analysis of two Beauveria bassiana strains grown on cuticle extracts of the silkworm uncovers their different metabolic response at early infection stage.

Beauveria bassiana is an important entomopathogenic fungus which not only widely distributes in the environment but also shows phenotypic diversity. However, the mechanism of pathogenic differences among natural B. bassiana strains has not been revealed at transcriptome-wide level. In the present study, in order to explore the mechanism, two B. bassiana strains with different pathogenicity were isolated from silkworms (Bombyx mori L.) and selected to analyze the gene expression of early stage by culturing on cuticle extracts of the silkworm and using RNA-sequencing technique. A total of 2108 up-regulated and 1115 down-regulated genes were identified in B. bassiana strain GXsk1011 (hyper-virulent strain) compared with B. bassiana strain GXtr1009 (hypo-virulent strain), respectively. The function categorization of differential expressed genes (DEGs) showed that most of them involved in metabolic process, biosynthesis of secondary metabolites, catalytic activity, and some involved in nutrition uptake, adhesion and host defense were also noted. Based on our data, distinct pathogenicity among different strains of B. bassiana may largely attribute to unique gene expression pattern which differed at very early infection process. Most of the genes involved in conidia adhesion, cuticle degradation and fungal growth were up-regulated in hyper-virulent B. bassiana strain GXsk1011. Furthermore, in combination with fungal growth analysis, our research provided a clue that fungal growth may also play an important role during early infection process. The results will help to explain why different B. bassiana strains show distinct pathogenicity on the same host even under same condition. Moreover, the transcriptome data were also useful for screening potential virulence factors.