Synthesis cost is a hidden driver of convergent amino acid composition in plastid ribosomal proteins

Protein evolution is a walk in the evolutionary space directed by mutation and selection. While functional and structural constraints serve as the main determinant of amino acid substitution in most proteins, synthesis cost and mutational bias can also alter the direction and rate of amino acid evolution, especially in systems experiencing relaxed selection. Here, we focused on the highly expressed plastid ribosomal proteins (PRP), which comprise 58 conserved proteins encoded by both plastid and nuclear genomes. Relaxed selection has been repeatedly identified in three distantly related plant lineages, providing a valuable comparative framework to investigate the significance of synthesis cost and mutation. We first demonstrated that the hemiparasitic tribe Cymbarieae (Orobanchaceae) represented a new case where concerted cyto-nuclear rate elevation occurs in their PRP. Further investigation revealed convergent shifts in amino acid composition in all four plant lineages attributable to arginine-to-lysine and methionine-to-isoleucine/valine/leucine substitutions. The replacement residues were biophysically similar but had lower molecular weight and shorter side chains, which significantly destabilized protein folding as demonstrated by protein structure modeling. We found that the composition shifts ran counter to the expectation of mutational bias but were consistent with the expectation of synthesis cost minimization, which is potentially adaptive for highly expressed PRP. Further, cost minimization significantly influenced all conservative substitutions between biophysically similar amino acids but was absent in non-conservative substitutions. We thus propose cost minimization as a secondary selective drive for protein evolution in PRP, unmasked in lineages and sites with relaxed selection on their function.