The evolution of Ca2+-ATPases across plants with profiles in Rhododendron and the function of key members in alleviating high calcium stress.

Ca2+-ATPase (CAP) is a key Ca2+ efflux protein in plants. Our previous research suggests that CAPs may play a crucial role in the adaptation of rhododendrons to high calcium environments. However, the evolution, variation, characteristic expression, and subfunctionalization of this gene family in Rhododendron remain unknown. Through the analysis of pan-genomes and pan-transcriptomes, we elucidated the systematic evolution of CAPs in plants, as well as their characteristic expression patterns in Rhododendron. During the evolutionary process from lower to higher plants, CAPs can be divided into six clades and exhibit structural conservation. CAPs have emerged and differentiated in lower plants such as algae, and they have undergone significant amplification in Eudicots plants like rhododendrons. Three Rhododendron species (Rhododendron bailiense, R. delavayi, and R. irroratum) located in the karst province of Guizhou in Southwest China exhibit the highest copies of CAPs, suggesting a strong association between CAP copy number variation and habitat, particularly in high calcium environments. Through multiple transcriptome analyses, we revealed that CAPs are induced under various environmental/developmental conditions (e.g. karst environments, high altitude, early flower development, hormones, etc.). Co-expression network analysis highlighted key members of calcineurin B-like protein (CBL) and CBL-interacting protein kinases (CIPK) that are associated with the high expression of CAPs. Experimental validation demonstrated that CAPb1 and CAPd1 significantly alleviate high calcium stress, and the CAPb1-CIPK1-CBL1 and CAPd1-CIPK2-CBL1 modules can further enhance the alleviation. These findings provide new insights into the evolution, characteristic expression, and function of CAPs, as well as new perspectives on the high calcium adaptability of rhododendrons.