Fungal co-expression network analyses identify pathogen gene modules associated with host insect invasion

The broad host range fungal insect pathogen, Beauveria bassiana, has been commercialized as an alternative to chemical insecticides for pest control worldwide. B. bassiana represents a unique model system with which to examine host-pathogen interactions, and a wide range of genes and processes have been studied. However, significant aspects of virulence, particularly on the genomic scale, remain poorly studied. Here, we have combined available transcriptomes with three newly generated data sets for a combined total analysis of 76 deep-sequenced samples covering growth, development, stress responses, and infection during the life cycle of B. bassiana. Co-expression network analyses resulted in the identification of gene modules enriched during two critical stages of the infection process, namely (i) cuticle penetration and (ii) in vivo hyphal body (dimorphic transition) growth capable of avoiding innate and humoral immune defenses. These analyses identify unique signatures of metabolism, signaling, secondary metabolite production, host defense suppression, membrane reorganization, effector production, and secretion for each stage, including genetic regulators and epigenetic patterns. These data provide a comprehensive framework for understanding and probing fungal adaptations to its pathogenic life cycle and expand the candidate repertoire for continued dissection of the host-pathogen interaction. IMPORTANCE Insect fungal pathogens have evolved unique strategies for overcoming host structural and immunological defenses that span from the sclerotized cuticle to innate and humoral cellular responses. Two critical stages of the infection process involve (i) cuticle penetration and (ii) immune evasion within the insect hemocoel. A set of 76 global transcriptomic data for B. bassiana that include the cuticle penetration and hemocoel growth stages were analyzed for patterns (gene modules) of expression, yielding unique insights into these different life stages. These analyses integrate gene networks involved in fungal development, stress response and pathogenesis to further the systematic understanding of the global processes integral to the unique adaptation employed by fungal pathogens of insects.