Unraveling a Tangled Skein: Evolutionary Analysis of the Bacterial Gibberellin Biosynthetic Operon

Gibberellin (GA) phytohormones are ubiquitous regulators of growth and developmental processes in vascular plants. The convergent evolution of GA production by plant-associated bacteria, including both symbiotic nitrogen-fixing rhizobia and phytopathogens, suggests that manipulation of GA signaling is a powerful mechanism for microbes to gain an advantage in these interactions. Although orthologous operons encode GA biosynthetic enzymes in both rhizobia and phytopathogens, notable genetic heterogeneity and scattered operon distribution in these lineages, including loss of the gene for the final biosynthetic step in most rhizobia, suggest varied functions for GA in these distinct plant-microbe interactions. Therefore, deciphering GA operon evolutionary history should provide crucial evidence toward understanding the distinct biological roles for bacterial GA production. To further establish the genetic composition of the GA operon, two operon-associated genes that exhibit limited distribution among rhizobia were biochemically characterized, verifying their roles in GA biosynthesis. This enabled employment of a maximum parsimony ancestral gene block reconstruction algorithm to characterize loss, gain, and horizontal gene transfer (HGT) of GA operon genes within alphaproteobacterial rhizobia, which exhibit the most heterogeneity among the bacteria containing this biosynthetic gene cluster. Collectively, this evolutionary analysis reveals a complex history for HGT of the entire GA operon, as well as the individual genes therein, and ultimately provides a basis for linking genetic content to bacterial GA functions in diverse plant-microbe interactions, including insight into the subtleties of the coevolving molecular interactions between rhizobia and their leguminous host plants. IMPORTANCE While production of phytohormones by plant-associated microbes has long been appreciated, identification of the gibberellin (GA) biosynthetic operon in plant-associated bacteria has revealed surprising genetic heterogeneity. Notably, this heterogeneity seems to be associated with the lifestyle of the microbe; while the GA operon in phytopathogenic bacteria does not seem to vary to any significant degree, thus enabling production of bioactive GA, symbiotic rhizobia exhibit a number of GA operon gene loss and gain events. This suggests that a unique set of selective pressures are exerted on this biosynthetic gene cluster in rhizobia. Through analysis of the evolutionary history of the GA operon in alphaproteobacterial rhizobia, which display substantial diversity in their GA operon structure and gene content, we provide insight into the effect of lifestyle and host interactions on the production of this phytohormone by plant-associated bacteria.