Single-cell transcriptomic analysis reveals the developmental trajectory and transcriptional regulatory networks of quinoa salt bladders

Abstract

Salt bladders, specialized structures on the surface of quinoa leaves, secrete Na⁺ to mitigate the effects of the plant from abiotic stresses, particularly salt exposure. Understanding the development of these structures is crucial for elucidating quinoa’s salt tolerance mechanisms. In this study, we employed transmission electron microscopy to detail cellular differentiation across the developmental stages of quinoa salt bladders. To further explore the developmental trajectory and underlying molecular mechanisms, we conducted single-cell RNA sequencing on quinoa protoplasts derived from young leaves. This allowed us to construct a cellular atlas, identifying 13 distinct cell clusters. Through pseudotime analysis, we mapped the developmental pathways of salt bladders and identified regulatory factors involved in cell fate decisions. GO and KEGG enrichment analyses, as well as experimental results, revealed the impacts of salt stress and the deprivation of sulfur and nitrogen on the development of quinoa salt bladders. Analysis of the transcription factor interaction network in pre-stalk cells (pre-SC), stalk cells (SC), and epidermal bladder cells (EBCs) indicated that TCP5, YAB5, NAC078, SCL8, GT-3B, and T1P17.40 play crucial roles in EBC development. Based on our findings, we developed an informative model elucidating salt bladder formation. This study provides a vital resource for mapping quinoa leaf cells and contributes to our understanding of its salt tolerance mechanisms.

References

[] AbeM, KatsumataH, KomedaY, TakahashiT. Regulation of shoot epidermal cell differentiation by a pair of homeodomain proteins in Arabidopsis. Development, 2003, 130(4): 635-643

[] AdamsMA, RichterA, HillAK, ColmerTD. Salt tolerance in Eucalyptus spp.: identity and response of putative osmolytes. Plant Cell Environ, 2005, 28(6): 772-787

[] AdolfVI, JacobsenS-E, ShabalaS. Salt tolerance mechanisms in quinoa (Chenopodium quinoa Willd.). Environ Exp Bot, 2013, 92: 43-54

[] AgarieS, ShimodaT, ShimizuY, BaumannK, SunagawaH, KondoA, UenoO, NakaharaT, NoseA, CushmanJC. Salt tolerance, salt accumulation, and ionic homeostasis in an epidermal bladder-cell-less mutant of the common ice plant Mesembryanthemum crystallinum. J Exp Bot, 2007, 58(8): 1957-1967

[] AliA, MaggioA, BressanRA, YunDJ. Role and functional differences of HKT1-Type transporters in plants under salt stress. Int J Mol Sci, 2019, 20(5): 1059

[] Al-MushhinAAM, QariSH, FakhrMA, AlnusairiGSH, AlnusaireTS, AA, A. L., Latef, A., Ali, O. M., Khan, A. A., Soliman, M. H.. Exogenous myo-inositol alleviates salt stress by enhancing antioxidants and membrane stability via the upregulation of stress responsive genes in Chenopodium quinoa L. Plants, 2021, 10(11): 2416

[] AngeliV, Miguel SilvaP, Crispim MassuelaD, KhanMW, HamarA, KhajeheiF, Graeff-HönningerS, PiattiC. Quinoa (Chenopodium quinoa Willd.): an overview of the potentials of the “Golden Grain” and socio-economic and environmental aspects of its cultivation and marketization. Foods, 2020, 9(2): 216

[] AshrafM, FooladMR. Roles of glycine betaine and proline in improving plant abiotic stress resistance. Environ Exp Bot, 2007, 59(2): 206-216

[] Bascuñán-Godoy L, Sanhueza C, Pinto K et al (2018) Nitrogen physiology of contrasting genotypes of Chenopodium quinoa Willd. (Amaranthaceae) Sci Rep 8(1):17524. https://doi.org/10.1038/s41598-018-34656-5

[] BazihizinaN, BöhmJ, MessererM, StigloherC, MüllerHM, CuinTA, MaierhoferT, CabotJ, MayerKFX, FellaC, HuangS, Al-RasheidKAS, AlquraishiS, BreadmoreM, MancusoS, ShabalaS, AcheP, ZhangH, ZhuJK, HedrichR, ScherzerS. Stalk cell polar ion transport provide for bladder-based salinity tolerance in Chenopodium quinoa. New Phytol, 2022, 235(5): 1822-1835

[] BianX, XingTL, YangY, FanJ, MaCM, LiuXF, WangY, HeYY, WangLD, WangB, ZhangN. Effect of soy protein isolate on physical properties of quinoa dough and gluten-free bread quality characteristics. J Sci Food Agric, 2023, 103(1): 118-124

[] Bilkova I, Kjaer A, van der Kooy F et al (2016) Effects of N fertilization on trichome density, leaf size and artemisinin production in artemisia annua leaves. Acta Hortic 1125:369–376. https://doi.org/10.17660/ActaHortic.2016.1125.48

[] BöhmJ, MessererM, MüllerHM, Scholz-StarkeJ, GradognaA, ScherzerS, MaierhoferT, BazihizinaN, ZhangH, StigloherC, AcheP, Al-RasheidKAS, MayerKFX, ShabalaS, CarpanetoA, HabererG, ZhuJ-K, HedrichR. Understanding the Molecular Basis of Salt Sequestration in Epidermal Bladder Cells of Chenopodium quinoa. Curr Biol, 2018, 28(19): 3075-3085.e3077

[] BoseJ, Rodrigo-MorenoA, LaiD, XieY, ShenW, ShabalaS. Rapid regulation of the plasma membrane H+-ATPase activity is essential to salinity tolerance in two halophyte species, Atriplex lentiformis and Chenopodium quinoa. Ann Bot, 2014, 115(3): 481-494

[] ClarkNM, NolanTM, WangP, SongG, MontesC, ValentineCT, GuoH, SozzaniR, YinY, WalleyJW. Integrated omics networks reveal the temporal signaling events of brassinosteroid response in Arabidopsis. Nat Commun, 2021, 12(1): 5858

[] Corteggiani CarpinelliE, TelatinA, VituloN, ForcatoC, D’AngeloM, SchiavonR, VezziA, GiacomettiGM, MorosinottoT, ValleG. Chromosome scale genome assembly and transcriptome profiling of Nannochloropsis gaditana in nitrogen depletion. Mol Plant, 2014, 7(2): 323-335

[] Cui Z, Huang X, Li M et al (2024) Integrated multi-omics analysis reveals genes involved in flavonoid biosynthesis and trichome development of Artemisia argyi. Plant Sci 346

[] DengYN, KashtohH, WangQ, ZhenGX, LiQY, TangLH, GaoHL, ZhangCR, QinL, SuM, LiF, HuangXH, WangYC, XieQ, ClarkeOB, HendricksonWA, ChenYH. Structure and activity of SLAC1 channels for stomatal signaling in leaves. Proc Natl Acad Sci U S A, 2021, 118(18): e2015151118

[] ElhindiKM, El-DinAS, ElgorbanAM. The impact of arbuscular mycorrhizal fungi in mitigating salt-induced adverse effects in sweet basil (Ocimum basilicum L.). Saudi J Biol Sci, 2017, 24(1): 170-179

[] Finken-EigenM, RöhrichtRA, KöhrerK. The VPS4 gene is involved in protein transport out of a yeast pre-vacuolar endosome-like compartment. Curr Genet, 1997, 31(6): 469-480

[] GrafBL, Rojas-SilvaP, RojoLE, Delatorre-HerreraJ, BaldeónME, RaskinI. Innovations in Health Value and Functional Food Development of Quinoa (Chenopodium quinoa Willd.). Compr Rev Food Sci F, 2015, 14(4): 431-445

[] HariadiY, MarandonK, TianY, JacobsenSE, ShabalaS. Ionic and osmotic relations in quinoa (Chenopodium quinoa Willd.) plants grown at various salinity levels. J Exp Bot, 2011, 62(1): 185-193

[] HuS, LiY, ShenJ. A diverse membrane interaction network for plant multivesicular bodies: roles in proteins vacuolar delivery and unconventional secretion. Front Plant Sci, 2020, 11: 425

[] HuangT, IrishVF. Temporal Control of Plant Organ Growth by TCP Transcription Factors. Curr Biol, 2015, 25(13): 1765-1770

[] HuangR, IrishVF. An epigenetic timer regulates the transition from cell division to cell expansion during Arabidopsis petal organogenesis. PLoS Genet, 2024, 20(3): e1011203

[] IchinoT, MaedaK, Hara-NishimuraI, ShimadaT. Arabidopsis ECHIDNA protein is involved in seed coloration, protein trafficking to vacuoles, and vacuolar biogenesis. J Exp Bot, 2020, 71(14): 3999-4009

[] JacobsenSE, MonterosC, CorcueraLJ, BravoLA, ChristiansenJL, MujicaA. Frost resistance mechanisms in quinoa (Chenopodium quinoa Willd.). Eur J Agron, 2007, 26(4): 471-475

[] JacobsenSE, MujicaA, JensenCR. The Resistance of Quinoa (Chenopodium quinoa Willd.) to Adverse Abiotic Factors. Food Rev Int, 2003, 19(1–2): 99-109

[] JedrzejasMJ. Structure, function, and evolution of phosphoglycerate mutases: comparison with fructose-2,6-bisphosphatase, acid phosphatase, and alkaline phosphatase. Prog Biophys Mol Biol, 2000, 73(2–4): 263-287

[] Kiani-PouyaA, RoessnerU, JayasingheNS, LutzA, RupasingheT, BazihizinaN, BohmJ, AlharbiS, HedrichR, ShabalaS. Epidermal bladder cells confer salinity stress tolerance in the halophyte quinoa and Atriplex species. Plant Cell Environ, 2017, 40(9): 1900-1915

[] KimJY, SymeonidiE, PangTY, DenyerT, WeidauerD, BezrutczykM, MirasM, ZöllnerN, HartwigT, WudickMM, LercherM, ChenLQ, TimmermansMCP, FrommerWB. Distinct identities of leaf phloem cells revealed by single cell transcriptomics. Plant Cell, 2021, 33(3): 511-530

[] KumaranMK, BowmanJL, SundaresanV. YABBY polarity genes mediate the repression of KNOX homeobox genes in Arabidopsis. Plant Cell, 2002, 14(11): 2761-2770

[] KuromoriT, FujitaM, UranoK, TanabataT, SugimotoE, ShinozakiK. Overexpression of AtABCG25 enhances the abscisic acid signal in guard cells and improves plant water use efficiency. Plant Sci, 2016, 251: 75-81

[] LiP, LiuQ, WeiY, XingC, XuZ, DingF, LiuY, LuQ, HuN, WangT, ZhuX, ChengS, LiZ, ZhaoZ, LiY, HanJ, CaiX, ZhouZ, WangK, ZhangB, LiuF, JinS, PengR. Transcriptional landscape of cotton roots in response to salt stress at single-cell resolution. Plant Commun, 2024, 5(2): 100740

[] LiuY, JiX, NieX, QuM, ZhengL, TanZ, ZhaoH, HuoL, LiuS, ZhangB, WangY. Arabidopsis AtbHLH112 regulates the expression of genes involved in abiotic stress tolerance by binding to their E-box and GCG-box motifs. New Phytol, 2015, 207(3): 692-709

[] LiuZ, GuoC, WuR, WangJ, ZhouY, YuX, ZhangY, ZhaoZ, LiuH, SunS, HuM, QinA, LiuY, YangJ, BawaG, SunX. Identification of the Regulators of Epidermis Development under Drought- and Salt-Stressed Conditions by Single-Cell RNA-Seq. Int J Mol Sci, 2022, 23(5): 2759

[] LiuZ, ZhouY, GuoJ, LiJ, TianZ, ZhuZ, WangJ, WuR, ZhangB, HuY, SunY, ShangguanY, LiW, LiT, HuY, GuoC, RochaixJD, MiaoY, SunX. Global Dynamic Molecular Profiling of Stomatal Lineage Cell Development by Single-Cell RNA Sequencing. Mol Plant, 2020, 13(8): 1178-1193

[] McFarlaneHE, ShinJJ, BirdDA, SamuelsAL. Arabidopsis ABCG transporters, which are required for export of diverse cuticular lipids, dimerize in different combinations. Plant Cell, 2010, 22(9): 3066-3075

[] MorohashiK, GrotewoldE. A systems approach reveals regulatory circuitry for Arabidopsis trichome initiation by the GL3 and GL1 selectors. PLoS Genet, 2009, 5(2): e1000396

[] Ohashi-ItoK, BergmannDC. Arabidopsis FAMA controls the final proliferation/differentiation switch during stomatal development. Plant Cell, 2006, 18(10): 2493-2505

[] OtterbachSL, KhouryH, RupasingheT, MendisH, KwanKH, LuiV, NateraSHA, KlaiberI, AllenNM, JarvisDE, TesterM, RoessnerU, SchmöckelSM. Characterization of epidermal bladder cells in Chenopodium quinoa. Plant Cell Environ, 2021, 44(12): 3606-3622

[] PillitteriLJ, SloanDB, BogenschutzNL, ToriiKU. Termination of asymmetric cell division and differentiation of stomata. Nature, 2007, 445(7127): 501-505

[] QiuQS, GuoY, DietrichMA, SchumakerKS, ZhuJK. Regulation of SOS1, a plasma membrane Na+/H+ exchanger in Arabidopsis thaliana, by SOS2 and SOS3. Proc Natl Acad Sci U S A, 2002, 99(12): 8436-8441

[] QiuX, MaoQ, TangY, WangL, ChawlaR, PlinerHA, TrapnellC. Reversed graph embedding resolves complex single-cell trajectories. Nat Methods, 2017, 14(10): 979-982

[] RasouliF, Kiani-PouyaA, ShabalaL, LiL, TahirA, YuM, HedrichR, ChenZ, WilsonR, ZhangH, ShabalaS. Salinity effects on guard cell proteome in Chenopodium quinoa. Int J Mol Sci, 2021, 22(1): 428

[] RothmanJH, StevensTH. Protein sorting in yeast: Mutants defective in vacuole biogenesis mislocalize vacuolar proteins into the late secretory pathway. Cell, 1986, 47(6): 1041-1051

[] RuffinoAMC, RosaM, HilalM, GonzálezJA, PradoFE. The role of cotyledon metabolism in the establishment of quinoa (Chenopodium quinoa) seedlings growing under salinity. Plant Soil, 2009, 326(1–2): 213-224

[] RyuKH, HuangL, KangHM, SchiefelbeinJ. Single-Cell RNA Sequencing Resolves Molecular Relationships Among Individual Plant Cells. Plant Physiol, 2019, 179(4): 1444-1456

[] SarojamR, SapplPG, GoldshmidtA, EfroniI, FloydSK, EshedY, BowmanJL. Differentiating Arabidopsis shoots from leaves by combined YABBY activities. Plant Cell, 2010, 22(7): 2113-2130

[] ShabalaL, MackayA, TianY, JacobsenS-E, ZhouD, ShabalaS. Oxidative stress protection and stomatal patterning as components of salinity tolerance mechanism in quinoa (Chenopodium quinoa). Physiol Plant, 2012, 146(1): 26-38

[] ShabalaS, BoseJ, HedrichR. Salt bladders: do they matter?. Trends Plant Sci, 2014, 19(11): 687-691

[] ShiH, QuinteroFJ, PardoJM, ZhuJK. The putative plasma membrane Na(+)/H(+) antiporter SOS1 controls long-distance Na(+) transport in plants. Plant Cell, 2002, 14(2): 465-477

[] SunY, HanY, ShengK, YangP, CaoY, LiH, ZhuQH, ChenJ, ZhuS, ZhaoT. Single-cell transcriptomic analysis reveals the developmental trajectory and transcriptional regulatory networks of pigment glands in Gossypium bickii. Mol Plant, 2023, 16(4): 694-708

[] SzabadosL, SavouréA. Proline: a multifunctional amino acid. Trends Plant Sci, 2010, 15(2): 89-97

[] TangY, ZhaoCY, TanST, XueHW. Arabidopsis Type II Phosphatidylinositol 4-Kinase PI4Kγ5 Regulates Auxin Biosynthesis and Leaf Margin Development through Interacting with Membrane-Bound Transcription Factor ANAC078. PLoS Genet, 2016, 12(8): e1006252

[] TrapnellC, CacchiarelliD, GrimsbyJ, PokharelP, LiS, MorseM, LennonNJ, LivakKJ, MikkelsenTS, RinnJL. The dynamics and regulators of cell fate decisions are revealed by pseudotemporal ordering of single cells. Nat Biotechnol, 2014, 32(4): 381-386

[] TsudaK, MaenoA, OtakeA, KatoK, TanakaW, HibaraKI, NonomuraKI. YABBY and diverged KNOX1 genes shape nodes and internodes in the stem. Science, 2024, 384(6701): 1241-1247

[] VogtT, IbdahM, SchmidtJ, WrayV, NimtzM, StrackD. Light-induced betacyanin and flavonol accumulation in bladder cells of Mesembryanthemum crystallinum. Phytochemistry, 1999, 52(4): 583-592

[] WangCQ, WangBS. Ca2+-Calmodulin is Involved in Betacyanin Accumulation Induced by Dark in C3 Halophyte Suaeda salsa. J Integr Plant Biol, 2007, 49(9): 1378-1385

[] WangY, GongX, LiuW, KongL, SiX, GuoS, SunJ. Gibberellin mediates spermidine-induced salt tolerance and the expression of GT-3b in cucumber. Plant Physiol Biochem, 2020, 152: 147-156

[] WatersS, GillihamM, HrmovaM. Plant High-Affinity Potassium (HKT) Transporters Involved in Salinity Tolerance: Structural Insights to Probe Differences in Ion Selectivity. Int J Mol Sci, 2013, 14(4): 7660-7680

[] XieH, WangQ, ZhangP, ZhangX, HuangT, GuoY, LiuJ, LiL, LiH, QinP. Transcriptomic and Metabolomic Analysis of the Response of Quinoa Seedlings to Low Temperatures. Biomolecules, 2022, 12(7): 977

[] YuH, ZhangL, WangW, TianP, WangW, WangK, GaoZ, LiuS, ZhangY, IrishVF, HuangT. TCP5 controls leaf margin development by regulating KNOX and BEL-like transcription factors in Arabidopsis. J Exp Bot, 2021, 72(5): 1809-1821

[] Zhang Y, Mutailifu A, Lan H (2022) Structure, development, and the salt response of salt bladders in Chenopodium album L. Front Plant Sci 13:989946. https://doi.org/10.3389/fpls.2022.989946

[] ZouC, ChenA, XiaoL, MullerHM, AcheP, HabererG, ZhangM, JiaW, DengP, HuangR, LangD, LiF, ZhanD, WuX, ZhangH, BohmJ, LiuR, ShabalaS, HedrichR, ZhuJK, ZhangH. A high-quality genome assembly of quinoa provides insights into the molecular basis of salt bladder-based salinity tolerance and the exceptional nutritional value. Cell Res, 2017, 27(11): 1327-1340