Genome sequence of Isaria javanica and comparative genome analysis insights into family S53 peptidase evolution in fungal entomopathogens

The fungus Isaria javanica is an important entomopathogen that parasitizes various insects and is effective for pest control. In this study, we sequenced and assembled the genomes (IJ1G and IJ2G) of two I. javanica strains isolated from different insects. The genomes were approximately 35 Mb in size with 11,441 and 11,143 protein-coding genes, respectively. Using a phylogenomic approach, we evaluated genome evolution across five entomopathogenic fungi in Cordycipitaceae. By comparative genome analysis, it was found that family S53 serine peptidases were expanded in Cordycipitaceae entomopathogens, particularly in I. javanica. Gene duplication events were identified based on phylogenetic relationships inferred from 82 S53 peptidases within six entomopathogenic fungal genomes. Moreover, we found that carbohydrate-active enzymes and proteinases were the largest secretory protein groups encoded in the I. javanica genome, especially chitinases (GH18), serine and aspartic peptidases (S53, S08, S10, A01). Pathogenesis-related genes and genes for bacterial-like toxins and secondary metabolites were also identified. By comparative transcriptome analysis, differentially expressed genes in response to insect nutrients (in vitro) were identified. Moreover, most S53 peptidases were detected to be significantly upregulated during the initial fungal infection process in insects (in vivo) by RT-qPCR. Our results provide new clues about understanding evolution of pathogenic proteases and may suggest that abundant S53 peptidases in the I. javanica genome may contribute to its effective parasitism on various insects.