Comparative genomic analysis of thermophilic fungi reveals convergent evolutionary adaptations and gene losses
Andrei S Steindorff,Maria Victoria Aguilar-Pontes,Aaron J Robinson,Bill Andreopoulos,Kurt LaButti,Alan Kuo,Stephen Mondo,Robert Riley,Robert Otillar,Sajeet Haridas,Anna Lipzen,Jane Grimwood,Jeremy Schmutz,Alicia Clum,Ian D Reid,Marie-Claude Moisan,Gregory Butler,Thi Truc Minh Nguyen,Ken Dewar,Gavin Conant,Elodie Drula,Bernard Henrissat,Colleen Hansel,Steven Singer,Miriam I Hutchinson,Ronald P de Vries,Donald O Natvig,Amy J Powell,Adrian Tsang,Igor V Grigoriev
DOI: https://doi.org/10.1038/s42003-024-06681-w
2024-09-12
Abstract:Thermophily is a trait scattered across the fungal tree of life, with its highest prevalence within three fungal families (Chaetomiaceae, Thermoascaceae, and Trichocomaceae), as well as some members of the phylum Mucoromycota. We examined 37 thermophilic and thermotolerant species and 42 mesophilic species for this study and identified thermophily as the ancestral state of all three prominent families of thermophilic fungi. Thermophilic fungal genomes were found to encode various thermostable enzymes, including carbohydrate-active enzymes such as endoxylanases, which are useful for many industrial applications. At the same time, the overall gene counts, especially in gene families responsible for microbial defense such as secondary metabolism, are reduced in thermophiles compared to mesophiles. We also found a reduction in the core genome size of thermophiles in both the Chaetomiaceae family and the Eurotiomycetes class. The Gene Ontology terms lost in thermophilic fungi include primary metabolism, transporters, UV response, and O-methyltransferases. Comparative genomics analysis also revealed higher GC content in the third base of codons (GC3) and a lower effective number of codons in fungal thermophiles than in both thermotolerant and mesophilic fungi. Furthermore, using the Support Vector Machine classifier, we identified several Pfam domains capable of discriminating between genomes of thermophiles and mesophiles with 94% accuracy. Using AlphaFold2 to predict protein structures of endoxylanases (GH10), we built a similarity network based on the structures. We found that the number of disulfide bonds appears important for protein structure, and the network clusters based on protein structures correlate with the optimal activity temperature. Thus, comparative genomics offers new insights into the biology, adaptation, and evolutionary history of thermophilic fungi while providing a parts list for bioengineering applications.