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

Abstract

Ca²⁺-ATPase (CAP) is a key Ca²⁺ 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.

References

[1] Ma Y, Mao X, Wang J, et al. Pervasive hybridization during evolutionary radiation of Rhododendron subgenus Hymenanthes in mountains of Southwest China. Natl Sci Rev. 2022; 9: nwac276

[2] Wen S, Cai X, Yang K, et al. Chromosome-level genome assembly of a rare karst-growing Rhododendron species provides insights into its evolution and environmental adaptation. J Syst Evol. 2025a; 63: 245-67

[3] Wen S, Zhao H, Qiao G, et al. Identification and characterization of genome-wide long terminal repeat retrotransposons provide an insight into elucidating the trait evolution of five Rhododendron species. Plant Biol. 2023a; 25: 813-28

[4] Xia XM, Du HL, Hu XD, et al. Genomic insights into adaptive evolution of the species-rich cosmopolitan plant genus Rhododendron. Cell Rep. 2024; 43: 114745

[5] Xia XM, Yang MQ, Li CL, et al. Spatiotemporal evolution of the global species diversity of Rhododendron. Mol Biol Evol. 2022; 39: msab314

[6] Chandan K, Gupta M, Ahmad A, et al. P-type calcium ATPases play important roles in biotic and abiotic stress signaling. Planta. 2024; 260: 37

[7] Costa A, Resentini F, Buratti S, et al. Plant Ca(2+)-ATPases: from biochemistry to signalling. Biochim Biophys Acta Mol Cell Res. 2023; 1870: 119508

[8] Wang G, Liu X, Liu Y, et al. Efficient enhancement of calcium content in apple by calcium-encapsulated carbon dots. Fruit Res. 2023a; 3: 0

[9] Demidchik V, Shabala S, Isayenkov S, et al. Calcium transport across plant membranes: mechanisms and functions. New Phytol. 2018; 220: 49-69

[10] Li Z, Harper JF, Weigand C, et al. Resting cytosol Ca2+ level maintained by Ca2+ pumps affects environmental responses in Arabidopsis. Plant Physiol. 2023; 191: 2534-50

[11] Rahmati IM, Brown E, Rosenberg A, et al. Arabidopsis Ca2+-ATPases 1, 2, and 7 in the endoplasmic reticulum contribute to growth and pollen fitness. Plant Physiol. 2021; 185: 1966-85

[12] Guan Y, Wu Y, Cao Z, et al. Island biogeography theory and the habitat heterogeneity jointly explain global patterns of Rhododendron diversity. Plant Divers. 2024; 46: 565-74

[13] Ma Y, Chamberlain D, Sun W, et al. A new species of Rhododendron (Ericaceae) from Baili Rhododendron nature reserve, NW Guizhou, China. Phytotaxa. 2015; 195: 197

[14] Liu C, Huang Y, Wu F, et al. Plant adaptability in karst regions. J Plant Res. 2021; 134: 889-906

[15] Wang C, Tang R-J, Kou S, et al. Mechanisms of calcium homeostasis orchestrate plant growth and immunity. Nature. 2024a; 627: 382-8

[16] Kovaka S, Ou S, Jenike KM, et al. Approaching complete genomes, transcriptomes and epi-omes with accurate long-read sequencing. Nat Methods. 2023; 20: 12-6

[17] Liu L, Shu Y, Wang Y, et al. The pan genome analysis of WOX gene family in apple and the two sides of MdWUS-1 in promoting leaf-borne shoot. Hortic Res. 2025; 12: uhaf117

[18] Nie S, Ma H-Y, Shi T-L, et al. Progress in phylogenetics, multi-omics and flower coloration studies in Rhododendron. Ornam Plant Res. 2024; 4: e003

[19] Nie S, Zhao S-W, Shi T-L, et al. Gapless genome assembly of azalea and multi-omics investigation into divergence between two species with distinct flower color. Hortic Res. 2023; 10: uhac241

[20] Shen JS, Lan L, Kan SL, et al. A haplotype-resolved genome for Rhododendron × pulchrum and the expression analysis of heat shock genes. J Syst Evol. 2023; 62: 489-504

[21] Wang X, Zhou P, Hu X, et al. T2T genome, pan-genome analysis, and heat stress response genes in Rhododendron species. iMeta. 2025b; 4: e70010

[22] Wu X, Zhang L, Wang X, et al. Evolutionary history of two evergreen Rhododendron species as revealed by chromosome-level genome assembly. Front Plant Sci. 2023; 14: 1123707

[23] Zhang L, Xu P, Cai Y, et al. The draft genome assembly of Rhododendron delavayi Franch. var. delavayi. Gigascience. 2017; 6: 1-11

[24] Luo J, Yuan C, Wang H, et al. Study on the genetic diversity characteristics of the endemic plant Rhododendron bailiense in Guizhou, China based on SNP molecular markers. Ecol Evol. 2025; 15: e70966

[25] Wen S, Cai X, Zhou K, et al. Metabolome and comparative genome provide insights into secondary metabolites generation of a rare karst-growing Rhododendron in vitro culture. Plant J. 2025b; 121: e17235

[26] Tong C, Jia Y, Hu H, et al. Pangenome and pantranscriptome as the new reference for gene-family characterization: a case study of basic helix-loop-helix (bHLH) genes in barley. Plant Commun. 2025; 6: 101190

[27] Wang J, Hu Y, Guo D, et al. Evolution and functional divergence of glycosyltransferase genes shaped the quality and cold tolerance of tea plants. Plant Cell. 2024b; 37: koae268

[28] Cao S, Ye Y, Zheng Z, et al. Aux/IAA gene family identification and analysis reveals roles in flower opening and abiotic stress response in Osmanthus fragrans. Ornam Plant Res. 2024; 4: 0

[29] Wang Q, Zhao X, Jiang Y, et al. Isolation, expression and functional analysis of the GbGGPPS1 gene in Ginkgo biloba. Fruit Res. 2023b; 3: 0

[30] Yang FS, Nie S, Liu H, et al. Chromosome-level genome assembly of a parent species of widely cultivated azaleas. Nat Commun. 2020; 11: 5269

[31] Su J, He B, Li P, et al. Overexpression of tonoplast Ca(2+)-ATPase in guard cells synergistically enhances stomatal opening and drought tolerance. J Integr Plant Biol. 2024; 66: 1587-602

[32] Zhang Y, Kobayashi C, Cai X, et al. Multiple sub-state structures of SERCA2b reveal conformational overlap at transition steps during the catalytic cycle. Cell Rep. 2022; 41: 111760

[33] Adimassu Z, Langan S, Johnston R, et al. Impacts of soil and water conservation practices on crop yield, run-off, soil loss and nutrient loss in Ethiopia: review and synthesis. Environ Manag. 2017; 59: 87-101

[34] Verslues PE, Bailey-Serres J, Brodersen C, et al. Burning questions for a warming and changing world: 15 unknowns in plant abiotic stress. Plant Cell. 2022; 35: 67-108

[35] Lyu ZY, Zhou XL, Wang SQ, et al. The first high-altitude autotetraploid haplotype-resolved genome assembled (Rhododendron nivale subsp. boreale) provides new insights into mountaintop adaptation. Gigascience. 2024; 13: giae052

[36] Wang X, Gao Y, Wu X, et al. High-quality evergreen azalea genome reveals tandem duplication-facilitated low-altitude adaptability and floral scent evolution. Plant Biotechnol J. 2021; 19: 2544-60

[37] Dai Y, Feng X, Liu Z, et al. miR1432 negatively regulates cold tolerance by targeting OsACAs. Plant Cell Environ. 2024; 47: 5443-56

[38] Shi T, Gao Z, Chen J, et al. Dosage sensitivity shapes balanced expression and gene longevity of homoeologs after whole-genome duplications in angiosperms. Plant Cell. 2024; 36: 4323-37

[39] Wen S, Zhao H, Zhang M, et al. IRAPs in combination with highly informative ISSRs confer effective potentials for genetic diversity and fidelity assessment in Rhododendron. Int J Mol Sci. 2023b; 24: 6902

[40] Long F, Wu H, Li H, et al. Genome-wide analysis of MYB transcription factors and screening of MYBs involved in the red color formation in Rhododendron delavayi. Int J Mol Sci. 2023; 24: 4641

[41] Zhou X-L, Wang S-Q, Dai N-F, et al. Interspecific variance of molecular response to external pH conditions in two Rhododendron species. Plant Soil. 2024; 504: 111-27

[42] Liu B, Wang X-Y, Cao Y, et al. Factors affecting freezing tolerance: a comparative transcriptomics study between field and artificial cold acclimations in overwintering evergreens. Plant J. 2020; 103: 2279-300

[43] Zhang X-M, Li J-T, Xia Y, et al. Early and late transcriptomic and metabolomic responses of Rhododendron ‘Xiaotaohong’ petals to infection with Alternaria sp. Int J Mol Sci. 2023; 24: 12695

[44] Shen J, Cheng H, Li X, et al. Beneficial effect of exogenously applied calcium chloride on the anatomy and fast chlorophyll fluorescence in Rhododendron × pulchrum leaves following short-term heat stress treatment. Agronomy. 2022; 12: 3226

[45] Liu Z, Zhang C, He J, et al. plantGIR: a genomic database of plants. Hortic Res. 2024; 11: uhae342

[46] Wang X, Wei Y, Liu Z, et al. TEGR: A comprehensive Ericaceae genome resource database. J Integr Agric. 2025a; 24: 1140-51

[47] Cantalapiedra CP, Hernandez-Plaza A, Letunic I, et al. eggNOG-mapper v2: functional annotation, orthology assignments, and domain prediction at the metagenomic scale. Mol Biol Evol. 2021; 38: 5825-9

[48] Qiao X, Li QH, Yin H, et al. Gene duplication and evolution in recurring polyploidization-diploidization cycles in plants. Genome Biol. 2019; 20: 38

[49] Wang D, Zhang Y, Zhang Z, et al. KaKs_Calculator 2.0: a toolkit incorporating gamma-series methods and sliding window strategies. Genomics Proteomics Bioinformatics. 2010; 8: 77-80

[50] Wang Y, Tang H, Debarry JD, et al. MCScanX: a toolkit for detection and evolutionary analysis of gene synteny and collinearity. Nucleic Acids Res. 2012; 40: e49

[51] Li H, Durbin R . Fast and accurate long-read alignment with burrows-wheeler transform. Bioinformatics. 2010; 26: 589-95

[52] Talevich E, Shain AH, Botton T, et al. CNVkit: genome-wide copy number detection and visualization from targeted DNA sequencing. PLoS Comput Biol. 2016; 12: e1004873

[53] Bolger AM, Lohse M, Usadel B . Trimmomatic: a flexible trimmer for Illumina sequence data. Bioinformatics. 2014; 30: 2114-20

[54] Kim D, Paggi JM, Park C, et al. Graph-based genome alignment and genotyping with HISAT2 and HISAT-genotype. Nat Biotechnol. 2019; 37: 907-15

[55] Ye Q, Zhang L, Li Q, et al. Genome and GWAS analysis identified genes significantly related to phenotypic state of Rhododendron bark. Hortic Res. 2024; 11: uhae008

[56] Pertea M, Pertea GM, Antonescu CM, et al. StringTie enables improved reconstruction of a transcriptome from RNA-seq reads. Nat Biotechnol. 2015; 33: 290-5

[57] Chen C, Wu Y, Li J, et al. TBtools-II: a “one for all, all for one” bioinformatics platform for biological big-data mining. Mol Plant. 2023; 16: 1733-42

[58] Wang X, Liang S, Yang W, et al. MetMiner: a user-friendly pipeline for large-scale plant metabolomics data analysis. J Integr Plant Biol. 2024c; 66: 2329-45

[59] Abramson J, Adler J, Dunger J, et al. Accurate structure prediction of biomolecular interactions with AlphaFold 3. Nature. 2024; 630: 493-500

[60] Schrodinger L. The PyMOL molecular graphics system. Version. 2015; 1: 8