Pfam domain adaptation profiles reflect plant species' evolutionary history

The increase in global temperatures predicted by climate change models presents a serious problem for agriculture because high temperatures reduce crop yields. Protein biochemistry is at the core of plant heat stress response, and understanding the interactions between protein biochemistry and temperature will be key to developing heat-tolerant crop varieties. Current experimental studies of proteome-wide plant thermostability are limited by the complexity of plant proteomes: evaluating function for thousands of proteins across a variety of temperatures is simply not feasible with existing technologies. In this paper, we use homologous prokaryote sequences to predict plant Pfam temperature adaptation and gain insights into how thermostability varies across the proteome for three species: maize, Arabidopsis, and poplar. We find that patterns of Pfam domain adaptation across organelles are consistent and highly significant between species, with cytosolic proteins having the largest range of predicted Pfam stabilities and a long tail of highly-stable ribosomal proteins. Pfam adaptation in leaf and root organs varies between species, and maize root proteins have more low-temperature Pfam domains than do Arabidopsis or poplar root proteins. Both poplar and maize populations have an excess of low-temperature mutations in Pfam domains, but only the mutations identified in poplar accessions have a negative effect on Pfam temperature adaptation overall. These Pfam domain adaptation profiles provide insight into how different plant structures adapt to their surrounding environment and can help inform breeding or protein editing strategies to produce heat-tolerant crops.