Abstract
Endofungal Mycetohabitans (formerly Burkholderia) spp. rely on a type III secretion system to deliver mostly unidentified effector proteins when colonizing their host fungus, Rhizopus microsporus. The one known secreted effector family from Mycetohabitans consists of homologues of transcription activator-like (TAL) effectors, which are used by plant pathogenic Xanthomonas and Ralstonia spp. to activate host genes that promote disease. These 'Burkholderia TAL-like (Btl)' proteins bind corresponding specific DNA sequences in a predictable manner, but their genomic target(s) and impact on transcription in the fungus are unknown. Recent phenotyping of Btl mutants of two Mycetohabitans strains revealed that the single Btl in one Mycetohabitans endofungorum strain enhances fungal membrane stress tolerance, while others in a Mycetohabitans rhizoxinica strain promote bacterial colonization of the fungus. The phenotypic diversity underscores the need to assess the sequence diversity and, given that sequence diversity translates to DNA targeting specificity, the functional diversity of Btl proteins. Using a dual approach to maximize capture of Btl protein sequences for our analysis, we sequenced and assembled nine Mycetohabitans spp. genomes using long-read PacBio technology and also mined available short-read Illumina fungal-bacterial metagenomes. We show that btl genes are present across diverse Mycetohabitans strains from Mucoromycota fungal hosts yet vary in sequences and predicted DNA binding specificity. Phylogenetic analysis revealed distinct clades of Btl proteins and suggested that Mycetohabitans might contain more species than previously recognized. Within our data set, Btl proteins were more conserved across M. rhizoxinica strains than across M. endofungorum, but there was also evidence of greater overall strain diversity within the latter clade. Overall, the results suggest that Btl proteins contribute to bacterial-fungal symbioses in myriad ways.