Single-nucleus RNA-seq and ATAC-seq analyses provide molecular insights into cadmium-stress response in alfalfa roots

Abstract

Soil cadmium (Cd) pollution threatens global food security. Elucidating plant cellular responses to Cd stress is critical for the breeding of low-Cd-accumulating crops. Here, we investigated Cd-responsive regulatory mechanisms in alfalfa at single-cell resolution and identified key Cd-associated genes. Eight major root cell types were annotated, with significant remodeling under Cd stress. Under Cd stress, root endodermal and phloem cells adopt distinct adaptive strategies: endodermal cells shift toward a Cd sequestration and detoxification state, whereas phloem cells exhibit a response gradient ranging from basic defense to systemic regulation. Integrated multiomics analyses revealed cell-type-specific genome-wide changes in chromatin accessibility, which positively correlated with gene expression-particularly in promoter regions. Key genes including MsGSH1, MsMT2A, MsHMP47, and MsABCC3 were shown to increase Cd tolerance in yeast. Coexpression network analysis revealed 10 cell-type-specific modules, with the calmodulin-like gene MsCML acting as a highly interconnected hub gene, whose overexpression significantly improved Cd tolerance. These findings provide valuable genetic resources and a theoretical basis for the precise breeding of low-Cd-accumulating forage, with implications for understanding Cd-responsive epigenetic and transcriptional regulation in plants.

References

[1] Rasin P, Ashwathi AV, Basheer SM, et al. Exposure to cadmium and its impacts on human health: a short review. J Hazard Mater Adv. 2025; 17: 100608

[2] Haider FU, Liqun C, Coulter JA, et al. Cadmium toxicity in plants: impacts and remediation strategies. Ecotoxicol Environ Saf. 2021; 211: 111887

[3] Mcgrath SP . Keeping toxic cadmium out of the food chain. Nat Food. 2022; 3: 569-70

[4] Zhang H, Lu L . Transcription factors involved in plant responses to cadmium-induced oxidative stress. Front Plant Sci. 2024; 15: 1397289

[5] Vert G, Grotz N, Dédaldéchamp F, et al. IRT1, an Arabidopsis transporter essential for iron uptake from the soil and for plant growth. Plant Cell. 2002; 14: 1223-33

[6] Li MZ, Hu DW, Liu XQ, et al. The OsZIP2 transporter is involved in root-to-shoot translocation and intervascular transfer of cadmium in rice. Plant Cell Environ. 2024; 47: 3865-81

[7] Cheng Y, Shuai W, Chen J, et al. Wheat ZIP3-2A encodes a metal transporter for Cd influx. Plant Physiol Biochem. 2025; 223: 109893

[8] Tan L, Zhu Y, Fan T, et al. OsZIP7 functions in xylem loading in roots and inter-vascular transfer in nodes to deliver Zn/Cd to grain in rice. Biochem Biophys Res Commun. 2019; 512: 112-8

[9] Chang J, Huang S, Yamaji N, et al. OsNRAMP1 transporter contributes to cadmium and manganese uptake in rice. Plant Cell Environ. 2020; 43: 2476-91

[10] Zhang W, Guan M, Chen M, et al. Mutation of OsNRAMP5 reduces cadmium xylem and phloem transport in rice plants and its physiological mechanism. Environ Pollut. 2024; 341: 122928

[11] Cailliatte R, Lapeyre B, Briat J, et al. The NRAMP6 metal transporter contributes to cadmium toxicity. Biochem J. 2009; 422: 217-28

[12] Cao L, Liu L, Zhang C, et al. The MYC2 and MYB43 transcription factors cooperate to repress HMA2 and HMA4 expression, altering cadmium tolerance in Arabidopsis thaliana. J Hazard Mater. 2024; 479: 135703

[13] Liu H, Zhao H, Wu L, et al. Heavy metal ATPase 3 (HMA3) confers cadmium hypertolerance on the cadmium/zinc hyperaccumulator Sedum plumbizincicola. New Phytol. 2017; 215: 687-98

[14] Fu S, Lu Y, Zhang X, et al. The ABC transporter ABCG36 is required for cadmium tolerance in rice. J Exp Bot. 2019; 70: 5909-18

[15] Kim DY, Bovet L, Maeshima M, et al. The ABC transporter AtPDR8 is a cadmium extrusion pump conferring heavy metal resistance. Plant J. 2007; 50: 207-18

[16] Park J, Song WY, Ko D, et al. The phytochelatin transporters AtABCC1 and AtABCC2 mediate tolerance to cadmium and mercury. Plant J. 2012; 69: 278-88

[17] Brunetti P, Zanella L, De Paolis A, et al. Cadmium-inducible expression of the ABC-type transporter AtABCC3 increases phytochelatin-mediated cadmium tolerance in Arabidopsis. J Exp Bot. 2015; 66: 3815-29

[18] Clemens S. Evolution and function of phytochelatin synthases. J Plant Physiol. 2006; 163: 319-32

[19] Wasinger VC, Bustamante S, Najib N, et al. Enzymes drive glutathione shunt to explain oxidative state using an in-parallel multi-omic method. Int J Mol Sci. 2025; 26: 3632

[20] Fu Z, Yao Y, Haq MZU, et al. Glutathione’s role in mitigating cadmium stress in Pogostemon cablin: insights from combined transcriptomic and metabolomic approaches. J Hazard Mater. 2025; 491: 137921

[21] Gupta DK, Palma JM, Corpas FJ . Redox State as a Central Regulator of Plant-Cell Stress Responses. Cham: Springer International Publishing; 2016

[22] Li F, Deng Y, Liu Y, et al. Arabidopsis transcription factor WRKY45 confers cadmium tolerance via activating PCS1 and PCS2 expression. J Hazard Mater. 2023; 460: 132496

[23] Liu Y, Liu L, Qi J, et al. Cadmium activates ZmMPK3-1 and ZmMPK6-1 via induction of reactive oxygen species in maize root. Biochem Biophys Res Commun. 2019; 516: 747-52

[24] Han Y, Chen Y, Yin S, et al. Over-expression of TaEXPB23, a wheat expansin gene, improves oxidative stress tolerance in transgenic tobacco plants. J Plant Physiol. 2015; 173: 62-71

[25] Takahashi R, Ishimaru Y, Shimo H, et al. The OsHMA2 transporter is involved in root-to-shoot translocation of Zn and Cd in rice. Plant Cell Environ. 2012; 35: 1948-57

[26] Uraguchi S, Kamiya T, Clemens S, et al. Characterization of OsLCT1, a cadmium transporter from indica rice (Oryza sativa). Physiol Plant. 2014; 151: 339-47

[27] Luo J, Huang J, Zeng D, et al. A defensin-like protein drives cadmium efflux and allocation in rice. Nat Commun. 2018; 9: 645

[28] Yadav AK, Shankar A, Jha SK, et al. A rice tonoplastic calcium exchanger, OsCCX2 mediates Ca²⁺/cation transport in yeast . Sci Rep. 2015; 5: 17117

[29] Dorrity MW, Alexandre CM, Hamm MO, et al. The regulatory landscape of Arabidopsis thaliana roots at single-cell resolution. Nat Commun. 2021; 12: 3334

[30] Liu Q, Liang Z, Feng D, et al. Transcriptional landscape of rice roots at the single-cell resolution. Mol Plant. 2021; 14: 384-94

[31] Wang T, Wang F, Deng S, et al. Single-cell transcriptomes reveal spatiotemporal heat stress response in maize roots. Nat Commun. 2024; 16: 177

[32] Liu H, Hu D, Du P, et al. Single-cell RNA-seq describes the transcriptome landscape and identifies critical transcription factors in the leaf blade of the allotetraploid peanut (Arachis hypogaea L.). Plant Biotechnol J. 2021; 19: 2261-76

[33] Li R, Wang Z, Wang J, et al. Combining single-cell RNA sequencing with spatial transcriptome analysis reveals dynamic molecular maps of cambium differentiation in the primary and secondary growth of trees. Plant Commun. 2023; 4: 100665

[34] Chen Y, Tong S, Jiang Y, et al. Transcriptional landscape of highly lignified poplar stems at single-cell resolution. Genome Biol. 2021; 22: 319

[35] Pacenza M, Muto A, Chiappetta A, et al. In Arabidopsis thaliana Cd differentially impacts on hormone genetic pathways in the methylation defective ddc mutant compared to wild type. Sci Rep. 2021; 11: 10965

[36] Liu Z, Liu Y, Zhou Y, et al. The crucial role of vascular tissues and COE2-mediated retrograde signaling in mitigating cadmium stress in Arabidopsis thaliana. Stress Biol. 2025; 5: 10

[37] Yu G, Xiang J, Liu J, et al. Single-cell atlases reveal leaf cell-type-specific regulation of metal transporters in the hyperaccumulator Sedum alfredii under cadmium stress. J Hazard Mater. 2024; 480: 136185

[38] Muto Y, Wilson PC, Ledru N, et al. Single cell transcriptional and chromatin accessibility profiling redefine cellular heterogeneity in the adult human kidney. Nat Commun. 2021; 12: 2190

[39] Marand AP, Chen Z, Gallavotti A, et al. A cis-regulatory atlas in maize at single-cell resolution. Cell. 2021; 184: 3041-55.e21

[40] Zhang X, Luo Z, Marand AP, et al. A spatially resolved multi-omic single-cell atlas of soybean development. Cell. 2025; 188: 550-67.e19

[41] Yan H, Mendieta JP, Zhang X, et al. A single-cell rice atlas integrates multi-species data to reveal cis-regulatory evolution. Nat Plants. 2025; 11: 2050-71

[42] Zhang L, He C, Lai Y, et al. Asymmetric gene expression and cell-type-specific regulatory networks in the root of bread wheat revealed by single-cell multiomics analysis. Genome Biol. 2023; 24: 65

[43] Cui Y, Su Y, Bian J, et al. Single-nucleus RNA and ATAC sequencing analyses provide molecular insights into early pod development of peanut fruit. Plant Commun. 2024; 5: 100979

[44] Farmer A, Thibivilliers S, Ryu KH, et al. Single-nucleus RNA and ATAC sequencing reveals the impact of chromatin accessibility on gene expression in Arabidopsis roots at the single cell level. Mol Plant. 2021; 14: 372-83

[45] Liu Q, Kang J, Du L, et al. Single-cell multiome reveals root hair-specific responses to salt stress. New Phytol. 2025; 246: 2634-51

[46] Li Y, Ma H, Wu Y, et al. Single-cell transcriptome atlas and regulatory dynamics in developing cotton anthers. Adv Sci. 2024; 11: 2304017

[47] Zhang C, Shi S, Wang B, et al. Physiological and biochemical changes in different drought-tolerant alfalfa (Medicago sativa L.) varieties under PEG-induced drought stress. Acta Physiol Plant. 2018; 40: 25

[48] Hussan MU, Mahmood U, Zarbakhsh S, et al. Mitigating Cd stress in alfalfa: the role of melatonin and nano-calcium oxide in enhancing photosynthesis and antioxidant defense. BMC Plant Biol. 2025; 25: 1099

[49] Chen H, Zeng Y, Yang Y, et al. Allele-aware chromosome-level genome assembly and efficient transgene-free genome editing for the autotetraploid cultivated alfalfa. Nat Commun. 2020; 11: 2494

[50] Zhang M, Chang M, Li H, et al. MsYSL6, a metal transporter gene of alfalfa, increases iron accumulation and benefits cadmium resistance. Plants. 2023; 12: 3485

[51] Xia H, Jing X, He H, et al. Genome-wide identification of the HIPPs gene family and functional validation of MsHIPP12 in enhancing cadmium tolerance in Medicago sativa. J Hazard Mater. 2025; 491: 137894

[52] Chen L, Li X, Liu H, et al. Comprehensive analysis of epigenetic modifications in alfalfa under cadmium stress. J Hazard Mater. 2025; 482: 136545

[53] Berni R, Luyckxc M, Xud X, et al. Reactive oxygen species and heavy metal stress in plants: impact on the cell wall and secondary metabolism. Environ Exp Bot. 2019; 161: 98-106

[54] Ding J, Adiconis X, Simmons SK, et al. Systematic comparison of single-cell and single-nucleus RNA-sequencing methods. Nat Biotechnol. 2020; 38: 737-46

[55] Wang K, Zhao C, Xiang S, et al. An optimized FACS-free single-nucleus RNA sequencing (snRNA-seq) method for plant science research. Plant Sci. 2023; 326: 111535

[56] Gu Q, Liu Y, Gu L, et al. GSH1 plays a crucial role in melatonin-mediated cadmium tolerance in Arabidopsis. Plant Cell Rep. 2025; 44: 235

[57] Bendall SC, Davis KL, Amir E, et al. Single-cell trajectory detection uncovers progression and regulatory coordination in human B cell development. Cell. 2014; 157: 714-25

[58] Butler A, Hoffman P, Smibert P, et al. Integrating single-cell transcriptomic data across different conditions, technologies, and species. Nat Biotechnol. 2018; 36: 411-20

[59] Chen Z, Ji Z, Ngiow SF, et al. TCF-1-centered transcriptional network drives an effector versus exhausted CD8 T cell-fate decision. Immunity. 2019; 51: 840-55.e5

[60] Yao C, Sun HW, Lacey NE, et al. Single-cell RNA-seq reveals TOX as a key regulator of CD8+ T cell persistence in chronic infection. Nat Immunol. 2019; 20: 890-901

[61] Herring CA, Banerjee A, Mckinley ET, et al. Unsupervised trajectory analysis of single-cell RNA-Seq and imaging data reveals alternative tuft cell origins in the gut. Cell Syst. 2018; 6: 37-51.e9

[62] van den Heuvel A, Mahfouz A, Kloet SL, et al. Single-cell RNA sequencing in facioscapulohumeral muscular dystrophy disease etiology and development. Hum Mol Genet. 2019; 28: 1064-75

[63] Guo WJ, Bundithya W, Goldsbrough PB . Characterization of the Arabidopsis metallothionein gene family: tissue-specific expression and induction during senescence and in response to copper. New Phytol. 2003; 159: 369-81

[64] Li J, Zhang M, Sun J, et al. Heavy metal stress-associated proteins in rice and Arabidopsis: genome-wide identification, phylogenetics, duplication, and expression profiles analysis. Front Genet. 2020; 11: 477

[65] Klemm SL, Shipony Z, Greenleaf WJ . Chromatin accessibility and the regulatory epigenome. Nat Rev Genet. 2019; 20: 207-20

[66] Pastor WA, Liu W, Chen D, et al. TFAP2C regulates transcription in human naive pluripotency by opening enhancers. Nat Cell Biol. 2018; 20: 553-64

[67] Cao S, Yu S, Li D, et al. Chromatin accessibility dynamics during chemical induction of pluripotency. Cell Stem Cell. 2018; 22: 529-42.e5

[68] Kusnadi EP, Timpone C, Topisirovic I, et al. Regulation of gene expression via translational buffering. Biochim Biophys Acta Mol Cell Res. 2022; 1869: 119140

[69] Detomaso D, Yosef N . Hotspot identifies informative gene modules across modalities of single-cell genomics. Cell Syst. 2021; 12: 446-56.e9

[70] Halema AA, El-Beltagi HS, Al-Dossary O, et al. Omics technology draws a comprehensive heavy metal resistance strategy in bacteria. World J Microbiol Biotechnol. 2024; 40: 193

[71] Wen B, Mei Z, Zeng C, et al. metaX: a flexible and comprehensive software for processing metabolomics data. BMC Bioinformatics. 2017; 18: 183

[72] Han L, Wei X, Liu C, et al. Cell transcriptomic atlas of the non-human primate Macaca fascicularis. Nature. 2022; 604: 723-31

[73] Lei Y, Cheng M, Li Z, et al. Spatially resolved gene regulatory and disease-related vulnerability map of the adult macaque cortex. Nat Commun. 2022; 13: 6747

[74] Liu C, Wu T, Fan F, et al. A portable and cost-effective microfluidic system for massively parallel single-cell transcriptome profiling. bioRxiv. 2019; 25: 1-7

[75] Yuhan Hao SH, Andersen-Nissen E, Gottardo R, et al. Integrated analysis of multimodal single-cell data. Cell. 2021; 184: 3573-87.e29

[76] Young MD, Behjati S . SoupX removes ambient RNA contamination from droplet-based single-cell RNA sequencing data. Gigascience. 2020; 9: 1-10

[77] Mcginnis CS, Murrow LM, Gartner ZJ . DoubletFinder: doublet detection in single-cell RNA sequencing data using artificial nearest neighbors. Cell Syst. 2019; 8: 329-37.e4

[78] Guo X, Wang Y, Zhao C, et al. An Arabidopsis single-nucleus atlas decodes leaf senescence and nutrient allocation. Cell. 2025; 188: 2856-71.e16

[79] Granja JM, Corces MR, Pierce SE, et al. ArchR is a scalable software package for integrative single-cell chromatin accessibility analysis. Nat Genet. 2021; 53: 403-11

[80] Korsunsky I, Millard N, Fan J, et al. Fast, sensitive and accurate integration of single-cell data with harmony. Nat Methods. 2019; 16: 1289-96

[81] Lex A, Gehlenborg N, Strobelt H, et al. UpSet: visualization of intersecting sets. IEEE Trans Vis Comput Graph. 2014; 20: 1983-92

[82] Yu G, Wang L, Han Y, et al. clusterProfiler: an R package for comparing biological themes among gene clusters. OMICS. 2012; 16: 284-7

[83] Qiu X, Mao Q, Tang Y, et al. Reversed graph embedding resolves complex single-cell trajectories. Nat Methods. 2017; 14: 979-82

[84] Gulati GS, Sikandar SS, Wesche DJ, et al. Single-cell transcriptional diversity is a hallmark of developmental potential. Science. 2020; 367: 405-11

[85] Modeling transcriptional regulation. Nat Genet. 2015; 47: 1

[86] Zhang X, Henriques R, Lin S, et al. Agrobacterium-mediated transformation of Arabidopsis thaliana using the floral dip method. Nat Protoc. 2006; 1: 641-6