Subfunctionalized expression drives evolutionary retention of ribosomal protein paralogs in vertebrates

The formation of paralogs through gene duplication is a core evolutionary process. For paralogs that encode components of protein complexes such as the ribosome, a central question is whether they encode functionally distinct proteins, or whether they exist to maintain appropriate total expression of equivalent proteins. Here, we systematically tested evolutionary models of paralog function using the mammalian ribosomal protein paralogs eS27 (Rps27) and eS27L (Rps27l) as a case study. We first showed that eS27 and eS27L have inversely correlated mRNA abundance across cell types, with the highest eS27 in lymphocytes and the highest eS27L in mammary alveolar cells and hepatocytes. By endogenously tagging the eS27 and eS27L proteins, we demonstrated that eS27- and eS27L- ribosomes associate preferentially with different transcripts. Furthermore, we generated murine eS27 and eS27L loss-of-function alleles that are homozygous lethal at different developmental stages. However, strikingly, we found that expressing eS27 protein from the endogenous eS27L locus, or vice versa, completely rescues loss-of-function lethality and yields mice with no detectable deficits. Together, these findings suggest that eS27 and eS27L are evolutionarily retained because their subfunctionalized expression patterns render both genes necessary to achieve the requisite total expression of two equivalent proteins across cell types. Our work represents the most in-depth characterization of a mammalian ribosomal protein paralog to date and highlights the importance of considering both protein function and expression when investigating paralogs.