Molecular evolution of Pr1 proteases depicts ongoing diversification in Metarhizium spp.

The Pr1 family of serine endopeptidases plays an important role in pathogenicity and virulence of entomopathogens such as Metarhizium anisopliae (Ascomycota: Hypocreales). These virulence factors allow for the penetration of the host cuticle, a vital step in the infective process of this fungus, which possesses 11 Pr1 isoforms (Pr1A through Pr1K). The family is divided into two classes with Class II (proteinase K-like) comprising 10 isoforms further split into three subfamilies. It is believed that these isoforms act synergistically and with other virulence factors, allowing pathogenicity to multiple hosts. As virulence coevolves through reciprocal selection with hosts, positive selection may lead to the evolution of new protease families or isoforms of extant ones that can withstand host defenses. This work tests this hypothesis in Class II Pr1 proteins, focusing on M. anisopliae, employing different methods for phylogenetic inference in amino acid and nucleotide datasets in multiple arrangements for Metarhizium spp. and related species. Phylogenies depict groups that match the taxonomy of their respective organisms with high statistical support, with minor discrepancies. Positively selected sites were identified in six out of ten Pr1 isoforms, most of them located in the proteolytic domain and spatially close to the catalytic residues. Moreover, there was evidence of functional divergence in the majority of pairwise comparisons. These results imply the existence of differential selective pressure acting on Pr1 proteins and a potential new isoform, likely affecting host specificities, virulence, or even adapting the organism to different host-independent lifestyles.