Buffering role of HSP shapes the molecular evolution of mammalian and human genomes at short and long-term scales

Heat shock proteins in parallel with their main and originally discovered function – maintenance of folded proteins under stressful conditions, can play also background buffering role – by folding proteins with an excess of slightly-deleterious nonsynonymous variants (SDNV). Here we tested several scenarios of this buffering role. On the comparative species scale, we demonstrated that low-Ne species are characterized by a higher expression level of hsp90 which can be explained by the excess of SDNV. On the comparative tissue level, we showed that long-lived tissues have also a higher hsp90 expression level, which can be advantageous to maintain the functionality of proteins. On the comparative gene level, we demonstrated that purifying selection of hsp90 in low-Ne-species did not relax as strongly as it happens for control genes, similar to hsp90. Additionally, we demonstrated that hsp clients versus non-clients are characterised by decreased level of selective constraints; demonstrate stronger relaxation of purifying selection in low-Ne species; have an excess of slightly-deleterious variants associated with complex disease phenotypes in humans; have an excess of pathological variants associated with clinical phenotypes in humans, suggesting that clients, being buffered by hsp90 can degenerate a bit more as compared to non-clients. Altogether, our results show that the secondary role of hsp, buffering of SDNV, is widespread and universal affecting properties of species, tissues and genes. A deep understanding of the buffering role of hsp90 will help to predict the deleterious effect of each variant in the human genome more precisely as well as will extend the application of the effectively-neutral theory of molecular evolution.