Massive Molecular Parallel Evolution of the HSP90AA1 Gene between High-elevation Anurans

HSP90 AA1 is part of the heat shock protein 90 gene family and has important functions against heat stress. We report a case of molecular level parallel evolution of the HSP90 AA1 gene in high elevation amphibians. HSP90 AA1 gene sequences of four high-elevation anurans, Bufo gargarizans, Nanorana parkeri, Rana kukunoris, and Scutiger boulengeri, were compared along with five of their low-elevation relatives. A total of 16 amino-acid sites were identified as parallel evolution between N. parkeri and R. kukunoris. We generated both model based(Zhang and Kumar’s test) and empirical data based(parallel/divergence plotting) null distributions for non-parallel evolution, and both methods clearly determined that the observed number of parallel substitutions were significantly more than the null expectation. Furthermore, on the HSP90 AA1 gene tree, N. parkeri and R. kukunoris formed a strongly supported clade that was away from their respective relatives. This study provides a clear case of molecular parallel evolution, which may have significant implications in understanding the genetic mechanisms of high-elevation adaptation.