NONLINEAR RELATIONSHIPS AMONG EVOLUTIONARY RATES IDENTIFY REGIONS OF FUNCTIONAL DIVERGENCE IN HEAT-SHOCK PROTEIN-70 GENES

The neutral theory predicts that, in comparisons among related genes, the number of amino acid replacements per site in a given gene region should be a linear function of that in another region of the same gene, unless the genes have diverged functionally in one region. Therefore, nonlinearity of this relationship can be used to identify regions of possible functional divergence among members of a multigene family. This method of analysis was applied to members of the heat-shock protein 70 (HSP70) gene family, which encode highly conserved ATP-dependent chaperone proteins found in all organisms. A nonlinear relationship was found between the rate of amino acid replacement in the conserved IA domain of the ATPase portion of the molecule and that in other ATPase domains and the peptide-binding domain. These results suggest that genes in the HSP70 subfamily (dnaK of bacteria and SSC1 of yeast) may have diverged functionally from other subfamilies in the ATPase domains, especially IIB, whereas SSB1 of yeast has diverged markedly in the peptide-binding domain. Functional divergence within these regions is consistent with what is known about functional differences between the HSP70 subfamilies in yeast.