Front. Plant Sci.Frontiers in Plant ScienceFront. Plant Sci.1664-462XFrontiers Media S.A.10.3389/fpls.2022.844430Plant ScienceOriginal ResearchThe Developmental Delay of Seedlings With Cotyledons Only Confers Stress Tolerance to Suaeda aralocaspica (Chenopodiaceae) by Unique Performance on Morphology, Physiology, and Gene ExpressionCaoJing†LiXiaorong†ChenLingHeMeixiangLanHaiyan*Xinjiang Key Laboratory of Biological Resources and Genetic Engineering, College of Life Science and Technology, Xinjiang University, Urumqi, China
Edited by: Mario A. Pagnotta, University of Tuscia, Italy
Reviewed by: Baris Uzilday, Ege University, Turkey; Manish Kumar Patel, Agricultural Research Organization (ARO), Israel
*Correspondence: Haiyan Lan lanhaiyan@xju.edu.cn
This article was submitted to Plant Abiotic Stress, a section of the journal Frontiers in Plant Science
†These authors have contributed equally to this work
This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
Cotyledons play an important role in seedling establishment, although they may just exist for a short time and become senescent upon the emergence of euphylla. So far, the detailed function of cotyledons has not been well understood. Suaeda aralocaspica is an annual halophyte distributed in cold deserts; its cotyledons could exist for a longer time, even last until maturity, and they must exert a unique function in seedling development. Therefore, in this study, we conducted a series of experiments to investigate the morphological and physiological performances of cotyledons under salt stress at different developmental stages. The results showed that the cotyledons kept growing slowly to maintain the normal physiological activities of seedlings by balancing phytohormone levels, accumulating osmoprotectants and antioxidants, and scavenging reactive oxygen species (ROS). Salt stress activated the expression of osmoprotectant-related genes and enhanced the accumulation of related primary metabolites. Furthermore, differentially expressed transcriptional profiles of the cotyledons were also analyzed by cDNA-AFLP to gain an understanding of cotyledons in response to development and salt stress, and the results revealed a progressive increase in the expression level of development-related genes, which accounted for a majority of the total tested TDFs. Meanwhile, key photosynthetic and important salt stress-related genes also actively responded. All these performances suggest that “big cotyledons” are experiencing a delayed but active developmental process, by which S. aralocaspica may survive the harsh condition of the seedling stage.
cotyledondevelopmental delayseedling growthstress toleranceSuaeda aralocaspica3126003732160057National Natural Science Foundation of China10.13039/501100001809Introduction
Seed germination and seedling development are the two most critical stages in the life cycle of a plant, can determine the establishment, persistence, and succession of plant population, and can even be significant in the maintenance of biodiversity (Kitajima and Fenner, 2000; Grime, 2001; Hanley et al., 2004). The plant seedling stage is defined as the beginning of time when a plant is able to perform photosynthesis and live individually with radicle growth and cotyledon development after seed germination (Fenner, 1987; Kitajima and Fenner, 2000), but the end of the seedling stage has not been conclusive yet. Cotyledons, as the most significant organ in the seedling stage, play a critical role in the establishment of a stronger plant (Hanley, 1998). Cotyledons provide various nutrients (fatty acids, carbohydrates, mineral nutrients, etc.) for seed germination and early seedling development (Lovell and Moore, 1971; Milberg and Lamont, 1997; Lamont and Groom, 2002; Shi et al., 2021). In a short period after seedling emergence, cotyledons can carry out photosynthesis immediately to provide energy for seedling growth (Xu et al., 2021), which is the principal source of nutrients before the euphylla appears. Since true leaves appear, with continuous exhaustion of minerals and reserves, as well as decline in photosynthetic capacity, most cotyledons gradually stay in wilting and become senescent (Keighley, 1984; Marek and Stewart, 1992). Generally, the cotyledons of most dicots can exist for 15–30 days (Hanley et al., 2004). However, although cotyledons remain healthy, most of the energy for seedling growth comes from the photosynthesis of true leaves (Ampofo et al., 1976a,b). So far, very limited reports have focused on cotyledon function and the effect of its damage or excision on seedling growth (Meyer, 1998; Hanley and May, 2006; Zhang et al., 2011; Zheng et al., 2011), and documentation on desert plant species is even fewer. Ruffino et al. (2010) preliminarily demonstrated that cotyledons of Chenopodium quinoa seedlings significantly contribute to the absorption of excess ions and accumulation of osmotic protectants under salt stress, which may enhance the adaptability of seedlings to heterogeneous habitats.
Suaeda aralocaspica (Bunge) Freitag and Schütze (Amaranthaceae) is an annual halophyte with a single-cell (SC) C4 photosynthesis pathway (Voznesenskaya et al., 2001) and distributed in the inland cold desert of Junggar Basin in China (Commissione Redactorum Florae Xinjiangensis, 1994). It can produce heteromorphic seeds, which have different dormancy properties and present different responses to salt stress (Wang et al., 2008), and has evolved strong adaptability in arid, saline-alkali, and barren environments (Cao et al., 2015). After many years of field investigation, we found that seedlings of S. aralocaspica developed very slowly in early spring, i.e., they grew for nearly 2 months with only two cotyledons in an extremely saline habitat, and the cotyledons could remain on the seedlings and become much larger without visible true leaves and even preserved on an adult plant for the rest of its life. However, the detailed function of a “big” cotyledon in normal growth and response to adverse stress remains largely unknown. Apparently, this unique developmental mode of S. aralocaspica in early stages must have its own internal mechanisms; further understanding of which should be important for us to elucidate how early seedlings adapt to the heterogeneous habitats in the desert.
Transcriptional profiling is a useful technique to determine key genes related to development and stress responses (Parmar et al., 2017; Fatemi et al., 2019). cDNA-amplified fragment length polymorphism (cDNA-AFLP) is one of the convenient and effective strategies for transcriptome analysis, which costs much lower with observable gel polymorphism in analyzing differential gene expression, even without any prerequisite knowledge of sequences at the molecular level (Bachem et al., 1996; Zhou et al., 2013; Xue et al., 2015). It is also efficient in increasing the resolution of detection of expression levels using a smaller amount of mRNA (Reijans et al., 2003; Amini et al., 2017), carrying out comprehensive and systematic analyses of transcriptome data (Gupta et al., 2012; Xue et al., 2014), construction of transcriptional profile (Cnudde et al., 2003; Sarosh and Meijer, 2007), and separation of differentially expressed genes (Wu et al., 2013; Lang et al., 2014) with high reproducibility, accuracy, reliability, and efficiency (Jayaraman et al., 2008; Li et al., 2015). The sensitivity and specificity of this method allows for detecting even the poorly expressed genes and distinguishing between homologous sequences (Breyne et al., 2003). Recently, an increasing number of reports on molecular mechanisms in response to stress was illustrated based on the cDNA-AFLP technique, and large amounts of stress-inducible genes have been identified (Urbanavičiute et al., 2021), e.g., genes encoding proteins in plant cell protection such as LEA (late embryogenesis abundant protein) (Jin et al., 2019), SOS1 (Na+/H+ antiporter protein) (Foster and Miklavcic, 2019), P5CS (pyrroline-5-carboxylate synthetase) (Székely et al., 2008), and Cu/Zn-SOD (Cu/Zn-superoxide dismutase) (Dvorák et al., 2021), etc.; genes encoding regulatory proteins in response to stress such as bHLH (basic helix-loop-helix transcription factor) (Wang et al., 2017), CDPK (calcium-dependent protein kinase) (Li et al., 2020), PP2C (protein phosphatase) (Shazadee et al., 2019), CaM (calmodulin) (Aldon et al., 2018), and NCED (9-cis-epoxycarotenoid dioxygenase) (Wang et al., 2021b), etc., which largely promote the understanding of the molecular biology of plants in development and response to stress. Nevertheless, the detailed knowledge of temporal alteration in the gene expression of halophyte seedlings is limited.
In this study, to clarify the function of cotyledons in seedling growth and response to stress in S. aralocaspica, we investigated morphological, physiological, and phytohormonal changes, and further analyzed transcriptome information (by cDNA-AFLP) and gene expression patterns during seedling (cotyledon) development, with which we expect to demonstrate the internal mechanism of the delayed development of seedlings of S. aralocaspica in adaptation to heterogeneous habitats.
Materials and MethodsSeedling Growth and Treatments
Mature seeds of S. aralocaspica were harvested from natural plants growing at Gurbantunggut Desert in Wujiaqu 103 regiment (44°29′N, 87°31′E; 430 mH) in October, in Xinjiang Uygur Autonomous Region, China. The seeds were air-dried and cleaned, and then stored at 4°C in sealed brown paper bags. Dimorphic (brown and black) seeds were separately sown in pots containing perlite: vermiculite (1:3, v/v) in a growth chamber and under conditions of a temperature regime of 20–25°C, 15–30% relative humidity, illumination of 150–700 μmol·m−2·s−1, and 14- to 16-h light/8- to 10-h dark photoperiod. The seedlings were cultivated with salt solutions containing 100, 300, or 500 mM NaCl [prepared with half-strength Hoagland solution (Arnon and Hoagland, 1940)] and applied at a 1-week interval using half-strength Hoagland solution only as the control. Cotyledons were harvested on days 3, 18, and 33 after emergence. Four biological replicates (in absence of salt stress) or three biological replicates (after salt treatment) of samples for physiological and gene expression analyses were collected from seedlings at different developmental stages, only 3- and 10-day-old seedlings were sampled under salt treatment for no enough seedlings left at 0 and 100 mM NaCl treatment thereafter. All the samples were immediately frozen in liquid nitrogen upon harvest and then stored at −80°C until use.
Observation of Micro-Structure of Apical Growth Point in Development
Paraffin sections were prepared following Li et al. (2021) to visualize the development of the apical growth point of seedlings in S. aralocaspica. For fixation, an FAA solution was used to fix the tissue of the joint part between two cotyledons and the hypocotyl. For dehydration, tissues were treated in 50% ethanol (30 min, two changes), followed by 1% Safranin O solution (prepared in 70% ethanol) staining for 8–10 h at room temperature, dehydrated in a series of ethanol from 80% to absolute ethanol, and then cleared at different concentrations of xylene and absolute ethanol. First-grade paraffin powder and a small amount of xylene were added to the tissue container in an oven at 38°C overnight. For paraffin inclusion and embedding, the above tissues were incubated at 56°C until the wax was completely melted, transferred into melted second-grade paraffin for 1 h, and then submerged in third-grade paraffin for 3 h. Finally, the melted paraffin containing tissues was poured into a paper tank. For sectioning, paraffin blocks were sliced into 6- to 12-μm sections with a Leica RM2126 microtome (Leica Biosystems, Germany). For deparaffinization and staining, glass slides with sections were successively incubated in xylene, xylene + ethanol, and ethanol, followed by 1% Safranin O staining overnight. Then, the slides were quickly and gently dipped into 1% Fast Green solution for about 10 s and immediately transferred into a series of ethanol, ethanol + xylene, and xylene. After mounting, the sections were inspected under a Leica DM3000 light microscope (Leica Microsystems, Germany) and photographed using the LAS V4.0 program.
Determination of Phytohormones
Seedlings were subjected to different concentrations of NaCl (0, 100, and 500 mM) treatment, and cotyledons were sampled at 0 (dry seeds), 7, 14, 21, 28, and 35 days after emergence. Four replicates of each treatment were homogenized in liquid nitrogen and transferred into 10 ml deionized distilled H2O (ddH2O) at 30°C for 3–4 h to extract phytohormones. High performance liquid chromatography (HPLC) (SCL-10AVP, Shimadzu, Japan) was conducted to determine four endogenous hormones: zeatin (ZT), gibberellin acid 3 (GA3), indole-3-acetic acid (IAA), and abscisic acid (ABA) (Pan et al., 2010). The standard curves, regression equation, and relationship coefficient were generated from standard solutions of ZT, GA3, IAA, ABA, and a series of dilutions (Supplementary Table 1). Conditions for chromatography were: (A) mobile phase: methanol: glacial acetic acid (0.75% in H2O) = 45:55; (B) flow rate: 0.7 ml·min−1; (C) column temperature: 35°C; (D) detection wavelength: 245 nm.
Measurement of Physiological ParametersDetermination of Tissue Water Content, and Na<sup>+</sup> and K<sup>+</sup> Concentrations
Fresh seedlings were weighed immediately to record the fresh weight (FW) and deactivated at 150°C for 20 min; then, they were transferred to an oven at 60°C until a constant dry weight (DW) was achieved. Tissue water content (TWC) was calculated as follows: TWC (%) = (FW–DW) × 100/FW (Tang et al., 2019). Fresh seedlings (1.5 g) were treated as above until DW was obtained. A dry plant material (0.3 g) was dissolved in HNO3, and a total volume of 25 ml was prepared by the addition of ddH2O. Na+ and K+ concentrations were measured by inductively coupled plasma atomic emission spectrometry (Agilent AA240, United States). HNO3 was used as the control (Al-Ashkar et al., 2021).
Determination of Concentration of Osmotic Protectants
For proline, seedlings (0.15 g) were homogenized in aqueous sulfosalicylic acid [5 ml, 3% (w/v)] and boiled for 10 min. The extract (200 μl), combined with a mixture containing acidic ninhydrin [300 μl, 2.5% (w/v)] and glacial acetic acid (200 μl), was incubated in boiling water bath for 40 min and then cooled on ice. After adding toluene (500 μl), the solution was shaken vigorously and the absorbance of the supernatant (200 μl) was measured at 520 nm. Toluene was used as the blank. L-Proline (0–20 μg·ml−1) was used to prepare the standard curve (Khan et al., 2017). For soluble sugar (SS), seedlings (0.12 g) were sheared into pieces and boiled in distilled H2O (dH2O, 10 ml) for 30 min. The extract (50 μl) was added to dH2O (150 μl), reacted with anthrone ethyl acetate [50 μl, 2% (w/v)] and concentrated sulfuric acid (500 μl), and then quickly boiled for 1 min. The absorbance at 630 nm was detected, and a standard curve (0–100 μg·ml−1 of glucose) was prepared to calculate SS content (Cui et al., 2018). For glycine betaine (GB), seedlings (0.15 g) were homogenized in liquid nitrogen and transferred into dH2O (1 ml) in a microtube, and then centrifuged at 10,000 g for 15 min; the supernatant (300 μl) mixed with saturated Reinecke's salt (500 μl) was incubated on ice for 1 h, followed by centrifuging at 10,000 g for 15 min. Finally, the precipitate was dissolved in acetone (70%). A calibration curve (0–200 μg·ml−1 of GB) was used to determine the concentration, and the absorbance was read at 525 nm (Zhao et al., 2016).
Measurement of <inline-formula><mml:math id="M1"><mml:mrow><mml:msubsup><mml:mtext>O</mml:mtext><mml:mn>2</mml:mn><mml:mo>−</mml:mo></mml:msubsup></mml:mrow></mml:math></inline-formula> Production, and H<sub>2</sub>O<sub>2</sub> and Malondialdehyde Concentrations
For O2− production, seedlings (0.15 g) were homogenized with phosphate-buffered saline (PBS, 1.3 ml, 50 mM, pH 7.8), and the slurry was centrifuged at 10,000 g and 4°C for 20 min. The mixture containing the supernatant (1 ml) and hydroxylamine (1 ml, 1 mM) was incubated at 25°C for 1 h, and then was added α-naphthalenamine (1 ml, 7 mM) and 4-aminobenzenesulfonic acid (1 ml, 17 mM). The reaction mixture (500 μl) was incubated at 25°C, and then was eliminated pigments with ether. A standard curve (0–50 μM of NaNO2) was prepared, and the absorbance was recorded at 530 nm (Imran et al., 2021). For H2O2 concentration, seedlings (0.15 g) were homogenized with acetone (2 ml) and then centrifuged at 3,000 g and 4°C for 10 min. The supernatant (1 ml) was added into a mixture of titanium sulfate [100 μl, 5% (w/v)] and ammonia (200 μl), and the precipitate that formed was washed to remove pigments and then dissolved in sulfuric acid (5 ml, 2 M). The absorbance was measured at 415 nm, and a calibration curve (0–100 μM of H2O2) was prepared to determine the concentration (Wang et al., 2021a). For malondialdehyde (MDA) concentration, seedlings (0.12 g) were ground into homogenate with trichloroacetic acid [TCA, 1 ml, 10% (w/v)] and centrifuged at 4,000 g for 10 min. The supernatant (300 μl) was mixed with thiobarbituric acid [TBA, 0.5% (w/v)] and then centrifuged. The absorbance of the supernatant was monitored separately at 450, 532, and 600 nm. The concentration of MDA was calculated using the equation (C) = 6.45 × (D532 – D600) – 0.56 × D450, where D450, D532, and D600 stand for absorbance values at 450, 532, and 600 nm, respectively (Motmainna et al., 2021).
Determination of the Content of Non-Enzymatic Antioxidants
Seedlings (0.15 g) were homogenized with TCA [1 ml, 5% (w/v)] and then centrifuged. For AsA content, the supernatant was mixed with TCA (200 μl, 612 mM), NaH2PO4 (100 μl, 150 mM), FeCl3 (100 μl, 184.9 mM), H3PO4 (200 μl, 4.49 M), and 2,2-biphenyl (200 μl, 214.8 mM) and then incubated. The absorbance value was recorded at 525 nm, and a standard curve (0–70 mM AsA) was plotted to calculate the AsA content (Rabelo et al., 2020). For GSH content, the supernatant was mixed with 5,5-dithio-bis (2-nitrobenzoic acid) (DTNB, 90 μl) and NaH2PO4 (130 μl, 150 mM), and then incubated. The absorbance value was determined at 412 nm, and the GSH content was calculated from a standard curve of GSH (0–0.12 mM) (Castro et al., 2021).
Measurement of the Activity of Antioxidant Enzymes
Seedlings (0.2 g) were homogenized and added into a PBS buffer (1.3 ml, 50 mM, pH 7.8). After centrifugation, the supernatant was used as crude enzyme for the analysis of different enzymatic activities immediately. For superoxide dismutase (SOD), a reaction mixture (3 ml) consisting of crude enzyme extract (50 μl), PBS (50 mM, pH 7.8), EDTA-Na2 (100 μM), methionine (130 mM), nitroblue tetrazolium (NBT, 750 μM), and riboflavin (2 μM) was incubated at 25°C for 25 min under a light source of 50 μmol·m−2·s−1, and then was measured the absorbance at 560 nm. One unit (U) of SOD activity was defined as the amount of enzyme for 50% inhibition of NBT photoreduction rate (Chrysargyris et al., 2018). For peroxidase (POD), a reaction mixture (3 ml) containing the crude enzyme extract (50 μl), PBS (1 ml, pH 7.0), H2O2 [1 ml, 0.3% (v/v)], and guaiacol [0.95 ml, 0.2% (v/v)] was prepared, the absorbance was recorded at 470 nm, and the activity was determined from the oxidation of guaiacol (Abid et al., 2018). For catalase (CAT), the reaction mixture (3 ml) contained the crude enzyme extract (100 μl), H2O2 [1 ml, 0.3% (v/v)], and H2O (1.9 ml). The activity was determined by measuring the decrease in H2O2 absorbance at 240 nm (Gao et al., 2019). The absorbance values of POD and CAT mixtures were recorded with a spectrophotometer (UV-3010; Shimadzu, Japan) within 1 min against a blank without the crude enzyme. The concentration of total proteins was assessed at 595 nm with the method of Bradford (1976). All the enzyme activities were expressed as unit per mg total protein.
Determination of Chlorophyll and Carotenoid Concentrations
Seedlings (0.15 g) were treated with 96% ethanol, and the resulting homogenate was centrifuged at 4,000 g for 10 min. The absorbance of the supernatant (300 μl) was measured at 470, 649, and 665 nm. The concentrations of chlorophylls (Chla and Chlb) and carotenoid (Car) were calculated as Chla = 13.95 D665 − 6.88D649, Chlb = 24.96 D649 − 7.32 D665, and Car = (1,000 D470 − 2.05 Chla − 114.8Chlb)/245, respectively. Total chlorophyll concentration was calculated by the addition of Chla and Chlb (Gao et al., 2021).
Estimation of the Activity of Key Photosynthetic Enzymes
Seedlings (0.4 g) were homogenized in liquid nitrogen with an extraction buffer (1.5 ml) containing Tris-H2SO4 (100 mM, pH 8.2), EDTA (1 mM), β-mercaptoethanol (7 mM), and glycerol [5%(v/v)], and centrifuged at 4,500 g and 4°C for 20 min. The supernatant was used as crude enzyme for assay of the activity of phosphoenolpyruvate carboxylase (PEPC) immediately. A reaction mixture (1 ml) containing NaHCO3 (143 μl, 70 mM), MgSO4 (143 μl, 70 mM), phosphoenolpyruvate (PEP) (286 μl, 14 mM), and NADH (429 μl, 5 mM) was incubated at 30°C for 30 min. After the addition of malate dehydrogenase (4 μl) and the crude enzyme (40 μl) to initiate the reaction, the mixture was recorded with an absorbance value for 3 min at 340 nm. One unit (U) of PEPC activity was defined as 1 nmol of NADH oxidation per minute per mg protein (You et al., 2020).
RNA Extraction
Total RNA was isolated from 0.1 g cotyledons of S. aralocaspica according to the manufacturer's instructions in the RNAprep Pure Plant Kit (DP432; Tiangen, Beijing, China). Cotyledons were sampled from seedlings derived from dimorphic seeds (brown and black seeds) at different developmental stages (3, 18, and 33 days after emergence) or at different NaCl concentrations (0, 100, 300, and 500 mM). The quantity and quality of total RNA were determined with a Nanodrop ND-1000 UV spectrophotometer (Thermo Fisher Scientific, Waltham, MA, United States).
Analysis of cDNA-AFLP
About 1 μg of mRNA [Poly(A) mRNA Purification Kit, SK8810; Sangon, Shanghai, China] was reversely transcribed into double-stranded cDNA using a PrimeScript™ Double Strand cDNA Synthesis kit (6111A; TaKaRa, Dalian, China). The cDNA (500 ng) was digested with EcoRI and MseI (P0525S and R0101V; NEB, Ipswich, MA, United States), and was then precipitated, resuspended, and ready for ligation with adaptors using a T4 DNA ligase (NEB, Ipswich, MA, United States).
Non-selective PCR pre-amplification was performed by mixing the above ligated products and non-selective primers EcoRI or MseI (Supplementary Table 2) under conditions of 94°C and 4 min for denaturation, then 30 cycles at 94°C for 30 s, 56°C for 30 s, and 72°C for 1 min, with a final 7 min extension at 72°C. PCR products were used as templates for subsequent selective amplification (SA) under the same conditions. The SA involved touchdown PCR, with the annealing temperature starting from 65 to 56°C, followed by 25 cycles of annealing temperature at 56°C. The AFLP protocol was the same as described by Banikamali et al. (2020), and a total of 25 selective primer combinations (i.e., 5 primers with EcoRI site and 5 primers with MseI site) were used for transcript profiling (Supplementary Table 2).
SA products (i.e., transcript-derived fragments, TDFs) were mixed with 0.5 times the volume of a denaturing buffer as described by Creste et al. (2001), and were separated on 6% polyacrylamide sequencing-type gel running at 850 V for 3–4 h after a 30-min pre-run. Gels were visualized by silver staining (with AgNO3) (Amini et al., 2017) and dried for scanning of images.
Isolation and Sequencing of TDFs
Differentially displayed polymorphic TDFs in silver-stained gels were identified and eluted according to the procedure suggested by Baisakh et al. (2006), and then re-amplified with appropriately selected primer combinations. Specific bands were sequenced after being recovered from the gels and purified using TIANgen Midi Gel Purification Kit (DP209; Tiangen, Beijing, China). Sequences of TDFs were analyzed for their homology against the publicly available genomic database of S. aralocaspica (Wang et al., 2019a) and non-redundant genes/ESTs/transcripts in the NCBI database (http://www.ncbi.nlm.nih.gov/BLAST) using the BLASTN and BLASTX algorithms. Functional classification of selected genes was performed using the MIPS database (http://mips.gsf.de/proj/funcatDB/search_main_frame.html).
Quantitative Real-Time PCR
qRT-PCR was performed with target TDFs and some primary metabolite-related genes to investigate the expression pattern in seedlings of S. aralocaspica at different developmental stages or under NaCl treatment. Each reversely transcriptional reaction was performed with total RNA (1 μg) in a final volume of 20 μl using the M-MLV RTase cDNA Synthesis Kit (D6130l TaKaRa, Dalian, China). Gene-specific primers of target TDFs and primary metabolite-related genes were designed using the Primer-Blast tool (Supplementary Table 3), and the β-actin of S. aralocaspica was used as internal reference (Cao et al., 2016). The relative expression level of target genes was quantified according to the mathematical model R=2−ΔΔCT (Shi and Chiang, 2005), where ΔΔCT = ΔCTtarget sample − ΔCTcontrol sample, and ΔCTsample= CTtest gene − CTreference gene. The final value of the relative quantification was described as a fold change of gene expression in the test sample compared to the control. Data are presented as mean ± standard error of four biological replicates and two technical replicates (n = 8).
Statistical Analysis
All data were plotted with the software of GraphPad Prism Version 7.0 for Windows (GraphPad Software, San Diego, United States) and analyzed using SPSS Version 26.0 for Windows (SPSS Inc., Chicago, United States). One way ANOVA and two-way ANOVA were conducted to test the significance of different treatments, and Tukey's HSD test was performed for multiple comparisons to determine significant differences between samples at 0.05, 0.01, 0.001, and 0.0001 significance levels. When the homogeneity of variance assumption was not met, differences were analyzed by Welch's ANOVA and Games Howell post-hoc test (McDonald, 2014).
ResultsThe Unique Morphology of the Seedlings in Development
In the natural habitat, S. aralocaspica seeds germinate at the beginning of March, and seedlings are gradually established from March to May. For extremely saline conditions, the seedlings may stay with only two cotyledons for nearly 2 months, presenting significant developmental retardation (Figure 1). Before true leaves start to visibly emerge, the two cotyledons were able to enlarge gradually to a much larger size compared with that of the seedling at emergence. Observation of the morphological change in seedling apical growth point (AGP) under a microscope revealed that the first leaf bud was visible on the 21st day, and appeared to grow faster at a higher salt concentration, i.e., 500 mM NaCl treatment promoted the generation of larger leaf primordium compared with the control (Figure 2). Such a phenomenon suggests that “large” cotyledons may be an adaptation strategy for early seedlings to get through harsh conditions.
Unique phenotype of “big cotyledons” of Suaeda aralocaspica in the natural habitat. C, cotyledon; GP, growth point; Mar, March; Apr, April. (A) Early-stage germinating seedling, (B) one-month-old seedling, (C) 2-month-old seedling, and (D) 3-month-old seedling.
Anatomic structures of the apical growth point (AGP), cotyledon (C), and euphylla (E) of S. aralocaspica after seedling emergence under different salt stress conditions. (A–D) without NaCl treatment; (E–H) 100 mM NaCl; (I–L) 500 mM NaCl. 7 d, 14 d, 21 d, 30 d: represent different developmental days. Scale bar is 50 μm.
Physiological Responses of the Seedlings at Different Developmental StagesChanges in Phytohormones
Four endogenous phytohormones, GA3, IAA, ZT and ABA, showed a general trend in developmental seedlings (0–35 days), in which relatively higher levels of the four hormones were observed in 7-day-old seedlings, followed by significant decreases in 14-day-old seedlings, and then were maintained relatively stable thereafter from 14 to 35 days, except for ABA on the 28th day (Figure 3). Among them, GA3 presented the highest level, IAA was moderate, ZT and ABA were relatively lower, and ABA content fluctuated slightly. Under salt treatment, lower concentration (100 mM) of NaCl stimulated the increase, while higher concentration (500 mM) resulted in decrease in the content of the four hormones, especially in seedlings younger than 14 days. Our results suggest that there are no significant fluctuations in GA3, IAA, and ZT contents during the later stage (14–35 days) of seedling development, at least in our experiment.
Changes in the content of four endogenous phytohormones during seedling development under salt treatment. (A) GA3; (B) IAA; (C) ZT; (D) ABA. Different letters in lowercase indicate a significant difference (P < 0.05) of the same salt concentration at different developmental stages. Values are means ± SE of three replicates.
Changes in Inorganic and Organic Osmolytes
The seedlings showed a significant accumulation of Na+ in cotyledons with the development progression (F2, 18 = 21.42, P < 0.0001), which differed between seedlings derived from two seed morphs (F1, 18 = 4.806, P = 0.0417) (Figure 4A); the level of K+ decreased in the brown seedlings from days 3 to 33 after emergence but increased in the black seedlings on day 18, and no significant change was present between days 3 and 33 in the black seedlings (Figure 4B). As a consequence, K+/Na+ ratios in both types of seedlings were significantly reduced (F2, 18 = 55.85, P < 0.0001) (Figure 4C). The TWC was also decreased by the developmental stage extension (F2, 18 = 9.049, P = 0.0019) (Figure 4D). A significant decrease in proline (F2, 18 = 231.9, P < 0.0001), GB (F2, 18 = 26.52, P < 0.0001), and SS (F2, 18 = 66.15, P < 0.0001) contents was observed in the two types of seedlings, especially proline content, which was reduced by about 9 folds in the brown and 18 folds in the black seedlings on day 33 (Figures 5A–C); whereas the total soluble protein was increased significantly from days 3 to 33 in both types of seedlings (F2, 18 = 12.3, P < 0.001) (Figure 5D). Between the different seedlings, the contents of GB (F1, 18 = 19.84, P < 0.001) and total soluble protein (F1, 18 = 4.497, P = 0.0481) were significantly higher in the black seedlings than that of the brown seedlings.
Changes in Na+, K+ and tissue water content (TWC) during seedling development. (A) Na+; (B) K+; (C) K+/Na+ ratio; (D) TWC. 3 d, 18 d, and 33 d: represent 3, 18, 33 days after seedling emergence, respectively. F-values are given when significance levels are reached (the subscript “D” represents “development”; “S”, “seed type”; *P < 0.05, **P < 0.01; ****P < 0.0001). Bars with different uppercase letters indicate significant differences (P < 0.05) according to Tukey's test. Values are means ± SE of four replicates. The Y-axis scale of TWC was adjusted to 150% to make visible of the difference comparison.
Changes in the concentration of osmotic protectants during seedling development. (A) Proline, (B) glycinbetaine (GB), (C) soluble sugar (SS), and (D) total soluble protein. 3 d, 18 d, and 33 d: represent 3, 18, and 33 days after seedling emergence, respectively. F-values are given when significance levels are reached (the subscript “D” represents “development”; “S”, “seed type”; *P < 0.05, ***P < 0.001; ****P < 0.0001). Bars with different uppercase letters indicate significant differences (P < 0.05) according to Tukey's test. Values are means ± SE of four replicates.
Levels of the Reactive Oxygen Species and Membrane Lipid Peroxidation
With the development progression of cotyledons, the production of O2− (F2, 18 = 4.795, P = 0.0214) and H2O2 (F2, 18 = 4.84, P = 0.0208) was gradually decreased (Figures 6A,B). MDA concentration was distinctly accumulated on day 18 (F2, 18 = 5.069, P = 0.0179) and significantly higher in the brown seedlings than that of the black seedlings (F1, 18 = 5.763, P = 0.0274) (Figure 6C). H2O2 accumulation had no significant difference between the two types of seedlings (F1, 18 = 0.1768, P = 0.6791); whereas O2− level (F2, 18 = 4.632, P = 0.0238) and MDA concentration (F2, 18 = 4.065, P = 0.0349) were significantly affected by the interaction of seedling type and developmental stage. Our results suggest that no stress significantly disturbs the normal growth of the two types of seedlings.
Changes in the ROS level and MDA concentration during seedling development. (A)O2−, (B) H2O2, and (C) MDA. 3 d, 18 d, and 33 d: represent 3, 18, and 33 days after seedling emergence, respectively. F-values are given when significance levels are reached (the subscript “D” represents “development”; “S”, “seed type”; *P < 0.05). Bars with different uppercase letters indicate significant differences (P < 0.05) according to Tukey's test. Values are means ± SE of four replicates.
Responses of Non-Enzymatic Antioxidants and Antioxidant Enzymes
The concentration of ASA (F2, 18 = 6.294, P = 0.0085) showed an increasing tendency with the development progression, and was significantly higher in the black seedlings than that of the brown seedlings (F1, 18 = 7.681, P = 0.0126) (Figure 7A). On the contrary, the contents of GSH and Car (F2, 18 = 18.54, P < 0.001) presented a decreasing trend (Figures 7B,C) and showed no significant difference between the two types of seedlings (F1, 18 = 0.3018, P = 0.5895 for GSH; F1, 18 = 0.1415, P = 0.7112 for Car) from days 3 to 33. The activities of SOD (F2, 18 = 5.487, P = 0.0138), POD (F2, 18 = 6.002, P = 0.0101) and CAT (F2, 18 = 253.7, P < 0.0001) were significantly reduced with the seedling development from days 3 to 33, especially that of CAT, but no significant difference was detected between the two types of seedlings (F1, 18 = 0.3135, P = 0.5825 for SOD; F1, 18 = 1.468, P = 0.2413 for POD; F1, 18 = 1.878, P = 0.1874 for CAT) (Figures 7D–F). Our data suggest that normal physiological activities are proceeding in seedling (cotyledon) development from days 3 to 33.
Changes in the concentration of non-enzymatic antioxidants and activity of antioxidant enzymes during seedling development. (A) AsA, (B) GSH, (C) Car, (D) SOD, (E) POD, and (F) CAT. 3 d, 18 d, and 33 d represent 3, 18, and 33 days after seedling emergence, respectively. F-values are given when significance levels are reached (the subscript “D” represents “development”; “S”, “seed type”; *P < 0.05; **P < 0.01; ****P < 0.0001). Bars with different uppercase letters indicate significant differences (P < 0.05) according to Tukey's test. Values are means ± SE of four replicates.
Changes in Chl Content and PEPC Activity
With the development progression from days 3 to 33, the contents of Chl a (F2, 18 = 7.125, P = 0.0053), Chl b (F2, 18 = 9.042, P = 0.0019), and total Chl (F2, 18 = 9.18, P = 0.0018) of both types of seedlings were significantly increased (Figures 8A–C). Corresponding to the increase in chlorophyll content, the activity of PEPC (F2, 18 = 122.1, P < 0.0001) was also significantly enhanced by seedling growth and achieved the highest value (about 24- and 27-fold higher than that of the control for the brown and black seedlings, respectively) on day 33 (Figure 8D). These results suggest that the photosynthetic function of the cotyledons is strengthened during seedling development.
Changes in the concentration of chlorophyll and activity of PEPC during seedling development. (A) Chl a, (B) Chl b, (C) total Chl, and (D) PEPC. 3 d, 18 d, 33 d represent 3, 18, and 33 days after seedling emergence, respectively. F-values are given when significance levels are reached (the subscript “D” represents “development”; “S”, “seed type”; *P < 0.05; **P < 0.01; ****P < 0.0001). Bars with different uppercase letters indicate significant differences (P < 0.05) according to Tukey's test. Values are means ± SE of four replicates.
Responses of the Seedlings Under Salt StressChanges in the Primary Metabolites
The accumulation of Na+ was significantly increased in cotyledons with elevated NaCl concentration (F3, 8 = 299.1, P < 0.0001 on day 3; F3, 8 = 808.3, P < 0.0001 on day 10), and was approximately 1.6 times higher in the 10-day-old seedlings than that of the 3-day-old seedlings (Figure 9A). On the contrary, salt stress induced a significant decrease in K+ concentration, especially in the 10-day-old seedlings under higher salinity (t4 = 27.67, P < 0.0001 at 300 mM) (Figure 9B), and K+/Na+ ratios were significantly reduced with the development progression and NaCl increasing (F3, 8 = 390.1, P < 0.0001 on day 3; F3, 8 = 1475.7, P < 0.0001 on day 10) (Figure 9C). The contents of proline (F3, 8 = 44.68, P < 0.0001 on day 3; F3, 8 = 15.52, P = 0.0011 on day 10), GB (F3, 8 = 24.08, P < 0.0001 on day 3; F3, 8 = 141.2, P < 0.0001 on day 10), and SS (F3, 8 = 118.39, P < 0.0001 on 10 d) were significantly enhanced by salt stress; whereas they were reduced by development extension (for proline, t4 = 6.239, P = 0.0034 at 100 mM; t4 = 3.143, P = 0.0348 at 300 mM; for SS, t4 = 9.041, P = 0.0008 at 0 mM; t4 = 6.197, P = 0.0034 at 500 mM). Total soluble protein was decreased with increasing salt concentration (F3, 8 = 13.18, P = 0.0018 on day 3; F3, 8 = 16.16, P = 0.0009 on day 10), while it was increased with the development progression at lower NaCl concentration (t4 = 5.998, P = 0.0039 at 100 mM) (Figure 10).
Changes in Na+ and K+ during seedling development under salt stress. (A) Na+, (B) K+, and (C) K+/Na+ ratio. 3 d and 10 d: represent 3 and 10 days after seedling emergence, respectively. 0, 100, 300, and 500: different NaCl concentrations (mM). Bars with different lowercase letters of the same developmental stage indicate significant differences (P < 0.05) at different concentrations; *** and **** indicate a significant difference between 3 d and 10 d at the same salt concentration at 0.001 and 0.0001 levels, respectively. Values are means ± SE of three replicates.
Changes in the concentration of osmotic protectants during seedling development under salt stress. (A) Proline, (B) glycinbetaine (GB), (C) soluble sugar (SS), and (D) total soluble protein. 3 d and 10 d represent 3 and 10 days after seedling emergence, respectively. 0, 100, 300, and 500: different NaCl concentrations (mM). Bars with different lowercase letters of the same developmental stage indicate significant differences (P < 0.05) at different salt concentrations; *, **, and *** indicate a significant difference between 3 d and 10 d at the same salt concentration at the 0.05, 0.01, and 0.001 levels, respectively. Values are means ± SE of three replicates.
Expression Patterns of Primary Metabolite-Related Genes
The expression of four tested genes was all upregulated with increasing NaCl concentration (except for V-PPase at 500 mM) (Figure 11), i.e., the vacuolar pyrophosphatase gene (V-PPase) (Welch's F3, 9.673 = 7.911, P = 0.0058 on day 3; Welch's F3, 9.62 = 38.386, P < 0.0001 on day 10), vacuolar ATP synthase gene (V-ATPase) (Welch's F3, 9.164 = 66.13, P < 0.0001 on day 3; F3, 20 = 7.98, P < 0.0001 on day 10), pyrroline-5-carboxylate synthase gene (P5CS) (Welch's F3, 20 = 191.4, P < 0.0001 on day 3; Welch's F3, 10.294 = 252.3, P < 0.0001 on day 10), and betaine aldehyde dehydrogenase gene (BADH) (F3, 20 = 66.08, P < 0.0001 on day 3; F3, 20 = 19.16, P < 0.0001 on day 10); especially with 300 mM NaCl (except for BADH), the transcript abundance of V-PPase (t10 = 5.771, P = 0.0002) and V-ATPase (t10 = 19.49, P < 0.0001) was significantly increased in the 10-day-old seedlings, and was 1.69- and 2.26-fold higher than that in the 3-day-old seedlings, respectively. Meanwhile, a significant reduction in P5CS expression level in the 10-day-old seedlings was observed under 0 mM (t10 = 17.15, P < 0.0001) and 500 mM NaCl treatments (t10 = 7.611, P < 0.0001) compared with that of the 3-day-old seedlings. The expression pattern of the above genes agreed with the accumulation tendency of primary metabolites involved.
Expression patterns of primary metabolite-related genes during seedling development under salt stress. (A) Vacuolar pyrophosphatase gene (V-PPase), (B) vacuolar ATP synthase gene (V-ATPase), (C) pyrroline-5-carboxylate synthase gene (P5CS), and (D) betaine aldehyde dehydrogenase gene (BADH). 3 d and 10 d represent 3 and 10 days after seedling emergence, respectively. 0, 100, 300, and 500: different NaCl concentrations (mM). Bars with different lowercase letters of the same developmental stage indicate significant differences (P < 0.05) at different salt concentrations; *, **, ***, and **** indicate a significant difference between 3 d and 10 d at the same salt concentration at the 0.05, 0.01, 0.001, and 0.0001 levels, respectively. Values are means ± SE of six replicates.
Identification and Analysis of Differentially Expressed Transcripts During Seedling Development
Differentially expressed TDFs in the two types of seedlings were analyzed with the cDNA-AFLP technique under salt treatment. A representative silver-staining PAGE gel of cDNA expression profiles is shown in Supplementary Figure S1. With 25 pairs of selective primer combinations, a total of 432 TDFs with an average of 18 distinct bands per primer pair were generated, of which more than 90% TDFs were identified between different developmental stages, and only 2% TDFs were differentially expressed between seedlings derived from the two seed morphs.
Among the 432 TDFs, 337 produced reliable sequences with a length ranging from 50 to 700 bp. Each sequence was identified by similarity search using the BLAST program against the S. aralocaspica genomic database and the GenBank non-redundant public sequence database. Sequence comparison of 328 reliable TDFs (9 had no significant matches) revealed that the majority was unique and shared homology with genes of known function in the database; out of these TDFs, 8% belonged to neither classified nor known function proteins. Based on putative function, the rest of the TDFs could be grouped into nine functional categories, transport (14%), transcription (13%), metabolism (12%), protein fate (11%), protein synthesis (11%), cell rescue, defense and virulence (9%), cell fate (9%), signal transduction (7%), and energy (6%) (Figure 12A; Supplementary Table S4), suggesting that substance transportation, gene transcription, protein synthesis, metabolism, etc. are actively in progress in seedling (cotyledon) development. A functional enrichment analysis based on databases of the GO (Gene Ontology) and KEGG (Kyoto Encyclopedia of Genes and Genomes) indicated that the TDFs are highly associated with activities proceeding during development, e.g., replication, recombination and repair, transcription, translation, carbohydrate transport, etc. (Figures 12B,C).
Distribution of differentially expressed TDFs in the cDNA-AFLP analysis during seedling development under salt stress in S. aralocaspica. (A) BLAST functional classification, (B) KEGG functional cluster, and (C) GO functional classification. “Signal trans.” represents signal transduction, “Cell rescue etc. “represents cell rescue, defense and virulence; “Protein syn.” represents protein synthesis.
Validation of Expression Patterns of TDFs During Seedling Development
Based on putative function, we selected 32 differentially expressed TDFs to validate the expression pattern by qRT-PCR analysis. Among them, 15 TDFs encode genes involved in seedling development, 12 encode genes of osmolytes, transporters, or kinases involved in balancing or regulating stress tolerance, and 5 encode photosynthetic enzyme genes including Rubisco (key in C3 cycle) and PEPC (key in C4 cycle).
Our results verified that most development-related genes were upregulated with seedling growth (Figure 13A), especially UVH1 (Welch's F3, 14.62 = 12.17, P = 0.0003), SYNC1 (Welch's F3, 14.64 = 22.32, P < 0.0001), and LST8 (Welch's F3, 14.98 = 15.14, P < 0.0001), the transcript abundance of which was significantly increased in cotyledons of the brown seedlings on day 33, and was 1.2-, 0.82- and 0.74-fold higher than that on day 3, respectively (Figure 14A). The photosynthetic genes had similar expression patterns with significantly higher expression levels observed in mature cotyledons (Figure 13B), especially the PEPC (Welch's F3, 12.36 = 98, P < 0.0001), LUT2 (Welch's F3, 14.22 = 39.27, P < 0.0001), and ACLA1 (Welch's F3, 13.84 = 21.98, P < 0.0001) with increases of 8.42-, 1.59- and 1.29-fold, respectively, on day 33 compared with day 3 in the brown seedlings (Figure 14B). The majority of salt response genes were up-regulated under salt stress in the two types of seedlings with the development progression (Figure 13C). Under 500 mM NaCl treatment, the transcripts of BADH4 (Welch's t7.448 = 9.114, P < 0.0001), SKD1 (Welch's t8.177 = 8.119, P = 0.0008), and MKK4 (Welch's t7.645 = 7.393, P < 0.0001) were increased by ~1.84–2.96-folds compared with the control, especially in cotyledons of the brown seedlings on day 3 (Figure 14C). In general, genes involved in development, photosynthesis, and salt tolerance actively participated in seedling (cotyledon) development from days 3 to 33; among them, photosynthetic or salt-stress response genes presented much greater involvement than development-related genes.
Heatmap of fold expression of 34 TDFs expressed in seedlings of S. aralocaspica at different developmental stages under salt stress. (A) Development-related genes, (B) photosynthesis-related genes, and (C) salt stress response genes. 3 d, 18 d, and 33 d represent 3, 18, and 33 days after seedling emergence, respectively. 0, 100, 300, 500: different NaCl concentrations (mM). Br, brown seedling; Bl, black seedling.
Relative expression of 21 interested TDFs in seedlings of S. aralocaspica at different developmental stages or under different salt stress. (A) Development-related genes, (B) photosynthesis-related genes, and (C) salt stress response genes. 3 d and 33 d represent 3 and 33 days after seedling emergence. 0, 500 mM: different NaCl concentrations. Br, brown seedling; Bl, black seedling. Bars with different lowercase letters indicate a significant differences (P < 0.05) of the same gene among different treatments according to Tukey's test. Values are means ± SE of eight replicates.
Discussion
Suaeda aralocaspica has developed various strategies to adapt to heterogeneous conditions at different developmental stages, e.g., seed heteromorphism in the germination stage, SC C4 assimilation pathway in the adult stage, growth enhancement by salt in the whole life cycle, etc. (Voznesenskaya et al., 2001; Wang et al., 2008; Cao et al., 2015, 2020). However, the specific mechanism for early seedling development in an extreme environment has not been well understood so far. It is known that the seedling stage is usually the most fragile period, and that well-established seedlings should facilitate the survival and successful propagation of the population. Therefore, it is necessary to understand the strategy for early seedling establishment. In the present study, we reported that seedlings of S. aralocaspica could tolerate the harsh habitat of early spring through a developmental retardation phenomenon. Its two cotyledons grew slowly but continuously for about 2 months, eventually forming a much bigger size without any visible euphylla. These two large cotyledons could even be existing on mature plants without wilting. To reveal the living situation of such big cotyledons, we investigated the cotyledons in terms of response to different developmental stages and salt stress based on morphological and physiological performances and gene expression. The results indicate that the cotyledons are able to maintain normal growth, physiological and biochemical activities, and active gene expression for a much longer time instead of senescence, which may be a unique strategy for the S. aralocaspica seedlings to survive in the harsh habitat during early developmental stages.
According to meteorological data from 2008 to 2012 of the natural habitat of S. aralocaspica, the average annual precipitation is about 165 mm but is usually much less during the first 4 months (Cao et al., 2020). The higher moisture of soil in spring along with a much higher evaporation rate promotes the seasonal accumulation of large amounts of salinity on soil surface, which creates a hostile environment for seedling growth (Ma et al., 2018). In the natural habitat, we found that the S. aralocaspica seedlings stayed with only two cotyledons without any visible true leaf for about 2 months, and that the cotyledons continuously grew to have a much larger size instead of becoming withered or senescent (Figure 1). Generally, cotyledons of most plants would gradually become exhausted with minerals and storage substances soon after the emergence of true leaves, accompanying with a sharp decrease in the biomass and photosynthetic capacity from the 6th day to the lowest in 20 days, and eventually enter the withered and senescent state (Ampofo et al., 1976a; Keighley, 1984; Jucknischke and Kutschera, 1998; Hanley et al., 2004). In this study, we found that the first leaf bud at the apical growth point (AGP) experienced much slower development within the first 2–3 weeks after seedling emergence, i.e., two cotyledons grew vertically without opening and extending, in which the lower end of the cotyledons (at the junction of two cotyledons and the hypocotyl) was tightly pressed on the top surface of the AGP, which may largely restrain leaf growth from the shoot tip (Figure 2). Moreover, the appearance of the flat AGP within the first 3 weeks also indicated inactive cell division nearby, because the shoot tip usually presents as a spherical structure and the active apical meristem layer leaves a growing shoot behind, pushing the former forward (Lyndon, 1998). On the contrary, in the cotyledons, levels of endogenous hormones GA3, IAA, and ZT had no significant reduction and ABA was not consistently increased from the 14th day after emergence in this study, corresponding to the fact that normal physiological activities were still performed in the cotyledons. The existing literature indicates that the energy reserved in the cotyledons may largely act as a potential nutrition pool and play an important role in supporting seedling development (Yi et al., 2015). Taken together, our results suggest that slow but sustainable growth of two cotyledons is in progress, and that an inhibitory effect of the cotyledons on the differentiation of AGP is probably applied; such a unique strategy may keep the seedlings alive in the harsh habitat of early spring while waiting for the favorable coming season (Ma et al., 2018).
The physiological and biochemical activities of plants seem to be associated with environmental factors that vary in different developmental stages (Liu et al., 2010; Kong et al., 2017). Primary metabolites such as organic solutes (e.g., proline, GB, and SS) are not only essential for plant growth and development, but also responsible for stress mitigation during abiotic stress (Chaudhary et al., 2020). To maintain the osmotic potential under saline conditions, halophytes typically accumulate organic solutes and inorganic ions (e.g., Na+, K+, and Cl−) in leaves (Guo et al., 2019). These osmoprotectants may provide protection for cells by raising the osmotic potential, stabilizing proteins and membranes, and maintaining the relatively higher water content under stressful conditions (Sadak and Bakry, 2020). In this study, cotyledons of non-salt-grown S. aralocaspica seedlings were able to maintain Na+ and K+ homeostasis for at least 4 weeks after emergence, whereas the level of organic solutes was significantly reduced in the cotyledons with the development progression, implying that seedlings (cotyledons) grow normally after emergence without requiring excess osmolytes in various biological processes. However, when treated with higher salinity, the seedlings of S. aralocaspica significantly accumulated inorganic (Na+) and organic (proline, GB, SS, and soluble proteins) osmolytes with the development progression, along with the upregulation of related genes (V-H+-ATPase, V-H+-PPase, P5CS, and BADH), suggesting that cotyledons can positively respond to environmental stimuli by gene activation and primary metabolite production; these results were in agreement with those found in 60-day-old adult plants of S. aralocaspica and other halophytes (Cao et al., 2015). Interestingly, the proline content in seedlings developed in spring (Figure 5A) was about 5 times lower than that in the winter seedlings (Figure 10A), which may be the effect of seasonal temperature fluctuation on the seedlings, and we also detected such a change in other experiments (data not shown). In Brassica oleracea, low temperature can also cause increase in the proline content of sprout (Šamec et al., 2022). Upon stress, plant cellular homeostasis is disrupted, accompanying with the generation of ROS (e.g., H2O2 and O2−), which may cause peroxidation of membrane lipids and increase MDA concentration (Liu et al., 2021). Excess ROS can be scavenged by plant defense systems involving non-enzymatic antioxidants (e.g., carotenoids, AsA, and GSH) and antioxidant enzymes (e.g., SOD, POD, and CAT) (ElSayed et al., 2020). In this study, the concentrations of O2−, H2O2, and MDA were higher in early stages and decreased gradually with seedling development, which is consistent with the decreasing trend of non-enzymatic antioxidant contents (except for AsA) and antioxidant enzyme activities. These results suggest that it normal for developing cotyledons to carry out various biological activities without needing to cope with the excess ROS produced in the aging process. A similar strategy has been found in Suaeda physophora, in which cotyledons play an important role in seedling establishment mainly by generating oxygen and compartmentalizing Na+ under salt stress (Zhou et al., 2014). Cotyledon greening is critical for the establishment of photosynthetic activity in the post-germination growth of seedlings (Ojeda et al., 2017). The natural senescent process of leaves has been characterized by destruction of chlorophyll (Lohman et al., 1994) and down-regulation of photosynthetic genes (Brown and Hudson, 2015). In this study, the content of total chlorophyll and activity of PEPC in cotyledons were significantly increased with seedling growth, indicating that the photosynthetic process is active in cotyledons with delayed development. Taken together, our results suggest that developmental retardation may assist seedlings in balancing between internal physiological homeostasis and external stressful pressure.
As a unique SC C4 plant species, S. aralocaspica carries out C4 and C3 cycles in an elongated and polarized chlorenchyma cell and has apparently higher photosynthetic efficiency than Kranz C4 species in the same family (Edwards et al., 2004; Smith et al., 2009; Liu et al., 2020). Moreover, S. aralocaspica operates the same SC C4 pathway in either cotyledon or euphylla (Voznesenskaya et al., 2004), which may confer an advantage on early seedlings to tolerate harsh environmental conditions in persistence with two “big cotyledons”. It has been known that C4 structure and biochemistry can be enhanced toward a full C4 photosynthetic system with the development progression (Voznesenskaya et al., 2003; Koteyeva et al., 2014); S. aralocaspica also experiences a gradual developmental transition of the C4 structure and photosynthetic enzymes from a young leaf to a mature leaf, or from the bottom to the top in a single leaf (Koteyeva et al., 2016; Liu et al., 2020). Cotyledons are a key photosynthetic organ in early seedlings of S. aralocaspica, in which photosynthesis enhancement may play an important role in assisting with development and response to stress. Our data on physiological and biochemical parameters, and gene expression profiling, could support this assumption, e.g., the transcriptional level and enzymatic activity of key photosynthetic genes (e.g., SaPEPC) were increased continuously and significantly during cotyledon development from days 3 to 33 (Figures 8, 11). Moreover, suppression of photorespiration is one of the most important events in C4 photosynthesis (Uzilday et al., 2012), which can effectively reduce the oxidative load of a cell (Uzilday et al., 2014). In this study, the accumulation of O2− and H2O2 was reduced with the progression of development (Figures 6, 7); correspondingly, the activity of the antioxidant enzymes (SOD, POD, and CAT) was decreased in the cotyledons during seedling development, which is matched with the process of establishment for a complete SCC4 photosynthetic system. This is rational when considering that the photorespiration in C4 plants is nearly blocked, thus also preventing the production of photorespiratory H2O2 (Uzilday et al., 2018).
Plants can regulate signal cascades and molecular networks, and eventually change gene expression profiles to respond to stress (Macovei and Tuteja, 2012). Large numbers of TDFs were differentially expressed in different developmental stages or diverse conditions, which were involved in various biological processes (Gupta et al., 2012). In this study, 328 TDFs were identified by differential expression in the early seedling stage under salt stress, and the GO enrichment and KEGG analysis showed that most of the TDFs were predominantly involved in transport, transcription, translation, carbohydrate transport, metabolism, protein synthesis, and protein fate (Figure 12), which are consistent with those found in Oncidium milliongolds and Lavandula angustifolia (Qian et al., 2014; Banikamali et al., 2020). It has been reported that genes involved in metabolism and energy are activated early in seedlings and persist in functioning during growth (Guo et al., 2009). Our results suggest that these TDF-related genes should be necessary for the regulation of various biological processes and synthesis of numerous substances in sustaining cotyledon growth and response to stress (Street et al., 2008; Liu et al., 2017; Wang et al., 2019b). Corresponding to the unique SC C4 system and high photosynthetic efficiency in S. aralocaspica, in this study, ~5% of the 328 TDFs were photosynthesis-related in the cotyledon, and half of which were associated with the C4 photosynthetic pathway. Our data indicate that active molecular biological events occur in cotyledons of early seedlings as requirements for development and photosynthesis.
To validate the expression profiles, we selected some TDFs from different functional groups as representatives for qRT-PCR analysis (Yang et al., 2012). Three sets of genes involved in development, photosynthesis, and salt stress response were analyzed, and the results revealed that most of the development-related (DR) genes were upregulated in seedling growth compared with those related to response to stress, especially the SYNC1 (asparaginyl-tRNA synthetase) and LST8 (lethal with SEC13 protein 8) genes, which play important roles in regulating amino acid accumulation and nutrition synthesis (Moreau et al., 2012; Arifin et al., 2019), whereas another DR gene, XDH1 (xanthine dehydrogenase) was downregulated in cotyledon development, which is generally upregulated in leaf senescence (Nakagawa et al., 2007), suggesting that cotyledons of S. aralocaspica seedlings are not aging. Genes associated with salt tolerance may limit the uptake and transport of ions and adjust osmotic and ionic balance (Qiu et al., 2017; Riccio-Rengifo et al., 2021). In this study, most of the salt response genes, e.g., BADH4 (betaine aldehyde dehydrogenase 4) and SKD1 (suppressor of K+ transport growth defect 1, encoding for salt-induced AAA-Type ATPase), were upregulated with increasing NaCl concentration. BADH is multifunctional and can increase salt tolerance by protecting the photosynthetic apparatus (Niazian et al., 2021). Salinity stress activates the SKD1 expression of halophytes in the regulation of ion transporters to avoid ion toxicity; whereas reduced SKD1 expression renders Arabidopsis more salt-sensitive (Jou et al., 2004, 2013). Photosynthesis efficiency is an important index for seedling development and measuring physiological sensitivity to abiotic stress (Wang et al., 2019b). In this study, most of the photosynthetic genes were upregulated in seedling growth, e.g., PEPC (phosphoenolpyruvate carboxylase), PPDK (pyruvate orthophosphate dikinase), and MDH (malate dehydrogenase), which are all essentially rate-limiting enzymes in the C4 photosynthetic pathway (Doubnerová et al., 2014; Cheng et al., 2016). Among them, PEPC was remarkably active (~9-fold more on day 33 than on day 3) in the cotyledons, suggesting that the C4 photosynthetic pathway is working with high efficiency in early seedling development even without true leaf growth. Taken together, our data revealed that the majority of tested genes could positively respond to development or salt stress, and that the seedlings of S. aralocaspica could normally grow and develop with only two cotyledons in their natural habitat.
Conclusions
The seedling stage is the most vulnerable time in the life cycle of a plant, especially under abiotic stress. In this study, we comprehensively investigated the strategy of early seedling (actually with two cotyledons only) survival under harsh conditions in S. aralocaspica from the perspective of morphological and physiological performances and gene expression. It was found that the seedlings of S. aralocaspica could survive for about 2 months with only two cotyledons, which grew slowly but continuously with a bigger size rather than becoming senesced in the early developmental stage. Meanwhile, the seedlings could positively respond to stress by accumulation of primary metabolites and activation of related gene expression. Such a delayed but persistently developmental morphology was in accordance with the anatomic structure (i.e., flat AGP for a longer time), physiological performance (e.g., homeostasis of phytohormone change, ion absorption and osmolyte accumulation, and active function of PEPC), and the upregulated genes in response to development, photosynthesis, and stress tolerance. This unique developmental pattern may be a smart strategy for S. aralocaspica seedlings to maintain the balance between development and energy consumption, allowing the early seedling to survive the harsh environment with higher probability. Our findings should contribute to further understanding of the mechanism for early seedlings to adapt to the heterogeneous habitat in desert plant species.
Data Availability Statement
All datasets generated for this study are included in the article/Supplementary Material, further inquiries can be directed to the corresponding author.
Author Contributions
HL, JC, and XL designed the experiments and methodology. JC and HL wrote the manuscript. JC, XL, LC, and MH conducted the experiments and collected the data. JC analyzed the data. All authors contributed critically to the manuscript and gave final approval for publication.
Funding
This study was supported by National Natural Science Foundation of China (32160057 and 31260037).
Conflict of Interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher's Note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
The authors thank the reviewers and all the editors for their helpful comments and suggestions on this manuscript.
Supplementary Material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fpls.2022.844430/full#supplementary-material
A representative result of the silver-stained cDNA-AFLP gel showing differentially expressed genes at different developmental stages in seedlings of S. aralocaspica under salt stress. (A) Seedlings from brown seeds, (B) seedlings from black seeds. 3 d, 18 d, and 33 d represent 3, 18, and 33 days after seedling emergence, respectively. 0, 1, 3, 5 represent 0, 100, 300, 500 mM NaCl treatments. The red arrowheads indicate TDFs with a differentially expressed pattern. Each group was amplified with a different pair of selective AFLP primers.
Relevant parameters of four phytohormones.
Sequences of the adaptors and primers used in cDNA-AFLP analysis.
qRT-PCR primer sequences used in the present study.
The 328 differential expression TDFs identified by cDNA-AFLP against known gene sequences using BLAST.
ReferencesAbidM.AliS.QiL. K.ZahoorR.TianZ.JiangD.. (2018). Physiological and biochemical changes during drought and recovery periods at tillering and jointing stages in wheat (Triticum aestivum L.). Sci. Rep. 8, 4615. 10.1038/s41598-018-21441-729545536Al-AshkarI.RomdhaneW. B.El-SaidR. A.GhazyA.AttiaK.Al-DossA. (2021). Agro-physiologic responses and stress-related gene expression of four doubled haploid wheat lines under salinity stress conditions. Biology10, 56. 10.3390/biology1001005633466713AldonD.MbengueM.MazarsC.GalaudJ. P. (2018). Calcium signalling in plant biotic interactions. Int. J. Mol. Sci.19, 665. 10.3390/ijms1903066529495448AminiS.Maali-AmiriR.MohammadiR.Kazemi-ShahandshtiS. S. (2017). cDNA-AFLP analysis of transcripts induced in chickpea plants by TiO2 nanoparticles during cold stress. Plant Physiol. Biochem.111, 39–49. 10.1016/j.plaphy.2016.11.01127907856AmpofoS. T.MooreK. G.LovellP. H. (1976a). Cotyledon photosynthesis during seedling development in acer. New Phytol. 76, 41–52. 10.1111/j.1469-8137.1976.tb01436.xAmpofoS. T.MooreK. G.LovellP. H. (1976b). The influence of leaves on cotyledon photosynthesis and export during seedling development in Acer. New Phytol. 76, 247–255. 10.1111/j.1469-8137.1976.tb01458.xArifinA. G.GondoT.AkashiR. (2019). Asparaginyl-tRNA synthetase gene (SYNC1) characterized by Lotus corniculatus FOX-superroot lines has effects on plant morphology and amino acid contents of seed in soybean. Plant Biotechnol. 36, 233–240. 10.5511/plantbiotechnology.19.1016a31983877ArnonD. I.HoaglandD. R. (1940). Crop production in artificial culture solutions and in soils with special reference to factors influencing yields and absorption of inorganic nutrients. Soil Sci.50, 463–485.BachemC. W. B.VanderhoevenR. S.DebruijnS. M.VreugdenhilD.ZabeauM. (1996). Visualization of differential gene expression using a novel method of RNA fingerprinting based on AFLP: analysis of gene expression during potato tuber development. Plant J.9, 745–753. 10.1046/j.1365-313X.1996.9050745.x8653120BaisakhN.SubudhiP. K.ParamiN. P. (2006). cDNA-AFLP analysis reveals differential gene expression in response to salt stress in a halophyte Spartina alterniflora Loisel. Plant Sci. 170, 1141–1149. 10.1016/j.plantsci.2006.02.001BanikamaliM.SoltanlooH.RamezanpourS. S.YamchiA.SorahinobarM. (2020). Identification of salinity responsive genes in lavender through cDNA-AFLP. Biotechnol. Rep.28, e00520. 10.1016/j.btre.2020.e0052032963973BradfordM. M. (1976). A rapid and sensitive method for the quantification ofmicrogram quantities of protein utilizing the principle of protein-dye binding. Anal. Biochem.72, 248–254. 10.1016/0003-2697(76)90527-3BreyneP.DreesenR.CannootB.RombautD.VandepoeleK.RombautsS.. (2003). Quantitative cDNA-AFLP analysis for genome-wide expression studies. Mol. Gen. Genomics269, 173–179. 10.1007/s00438-003-0830-612756529BrownA.V.HudsonK.A. (2015). Developmental profiling of gene expression in soybean trifoliate leaves and cotyledons. BMC Plant Biol.15, 169. 10.1186/s12870-015-0553-y26149852CaoJ.ChenL.WangJ.XingJ. J.LvX. Y.MaimaitijiangT.. (2020). Effects of genetic and environmental factors on variations of seed heteromorphism in Suaeda aralocaspica. AoB Plants12, plaa044. 10.1093/aobpla/plaa04433072248CaoJ.LvX. Y.ChenL.XingJ. J.LanH.Y. (2015). Effects of salinity on the growth, physiology and relevant gene expression of an annual halophyte grown from heteromorphic seeds. AoB Plants7, plv112. 10.1093/aobpla/plv11226386128CaoJ.WangL.LanH. Y. (2016). Validation of reference genesfor quantitative RT-RCR normalization in Suaeda aralocaspica, an annualhalophyte with heteromorphism and C4 pathway without Kranz anatomy. Peer J. 4, e1697. 10.7717/peerj.169726893974CastroR. I.Morales-QuintanaL.AlvaradoN.GuzmánL.Forero-DoriaO.Valenzuela-RiffoF.. (2021). Design and optimization of a self-assembling complex based on microencapsulated calcium alginate and glutathione (CAG) using response surface methodology. Polymers13, 2080. 10.3390/polym1313208034202791ChaudharyS.DeviP.BhardwajA.JhaU. C.SharmaK. D.PrasadP. V.V.. (2020). Identification and characterization of contrasting genotypes/cultivars for developing heat tolerance in agricultural crops: current status and prospects. Front. Plant Sci.11, 587264. 10.3389/fpls.2020.58726433193540ChengG.WangL.LanH.Y. (2016). Cloning of PEPC-1 from a C4 halophyte Suaeda aralocaspica without Kranz anatomy and its recombinant enzymatic activity in responses to abiotic stresses. Enzyme Microb. Tech.83, 57–67. 10.1016/j.enzmictec.2015.11.00626777251ChrysargyrisA.MichailidiE.TzortzakisN. (2018). Physiological and biochemical responses of Lavandula angustifolia to salinity under mineral foliar application. Front. Plant Sci. 9, 489. 10.3389/fpls.2018.0048929731759CnuddeF.MorettiC.PorcedduA.PezzottiM.GeratsT. (2003). Transcript profiling on developing Petunia hybrida floral organs. Sex. Plant Reprod.16, 77–85. 10.1007/s00497-003-0178-4Commissione Redactorum Florae Xinjiangensis (1994). Flora Xinjiangensis. Urumchi: Xinjiang Science & Technology & Hygiene Publishing House.CresteS.NetoA.T.FigueiraA. (2001). Detection of single sequence repeat polymorphisms in denaturing polyacrylamide sequencing gels by silver staining. Plant Mol. Biol. Rep.19, 299–306. 10.1007/BF0277282819442642CuiG.LiB.HeW.YinX.LiuS.LianL.. (2018). Physiological analysis of the effect of altitudinal gradients onLeymus secalinus on the Qinghai-Tibetan Plateau. PLoS ONE13:e0202881. 10.1371/journal.pone.020288130183737DoubnerováH. V.MiedzińskaL.DobráJ.VankovaR.RyšlaváH. (2014). Phosphoenolpyruvate carboxylase, NADP-malic enzyme, and pyruvate, phosphate dikinase are involved in the acclimation of Nicotiana tabacum L. to drought stress. J Plant Physiol.171, 19–25. 10.1016/j.jplph.2013.10.01724484954DvorákP.KrasylenkoY.ZeinerA.ŠamajJ.TakáčT. (2021). Signaling toward reactive oxygen species-scavenging enzymes in plants. Front. Plant Sci.11, 618835. 10.3389/fpls.2020.61883533597960EdwardsG. E.FranceschiV. R.VoznesenskayaE. V. (2004). Single-cell C4 photosynthesis versus the dual-cell (Kranz) paradigm. Annu. Rev. Plant Biol.55, 173–196. 10.1146/annurev.arplant.55.031903.14172515377218ElSayedA. I.BoulilaM.RafudeenM. S.MohamedA. H.SenguptaS.RadyM.. (2020). Melatonin regulatory mechanisms and phylogenetic analyses of melatonin biosynthesis related genes extracted from peanut under salinity stress. Plants9, 854. 10.3390/plants907085432640740FatemiF.Hashemi-PetroudiS. H.NematzadehG.AskariH.AbdollahiM. R. (2019). Exploiting differential gene expression to discover ionic and osmotic-associated transcripts in the halophyte grass Aeluropuslittoralis. Biol. Proced. Online21, 14. 10.1186/s12575-019-0103-331337987FennerM. (1987). Seedlings. New Phytol.106, 35–47.FosterK. J.MiklavcicS. J. (2019). A comprehensive biophysical model of ion and water transport in plant roots. II. Clarifying the roles of SOS1 in the salt-stress response in Arabidopsis. Front. Plant Sci.10, 1121. 10.3389/fpls.2019.0112131620152GaoM.SunH.ShiM.WuQ.JiD.WangB.. (2021). 2-Keto-L-gulonic acid improved the salt stress resistance of non-heading chinese cabbage by increasing L-ascorbic acid accumulation. Front. Plant Sci. 12, 697184. 10.3389/fpls.2021.69718434804078GaoS.LiuY. Y.JiangJ. Y.FuY.ZhaoL. X.LiC. Y.. (2019). Protective responses induced by chiral 3-dichloroacetyl oxazolidine safeners in maize (Zea mays L.) and the detoxification mechanism. Molecules24, 3060. 10.3390/molecules2417306031443550GrimeJ.P. (2001). Plant Strategies, Vegetation Processes, and Ecosystem Properties. Chichester: Wiley.GuoH.ZhangL.CuiY. N.WangS. M.BaoA. K. (2019). Identification of candidate genes related to salt tolerance of the secretohalophyteAtriplex canescens by transcriptomic analysis. BMC Plant Biol.19, 213. 10.1186/s12870-019-1827-631117942GuoY.Q.TianZ.Y.QinG.Y.YanD.L.ZhangJ.ZhouW.Z.. (2009). Gene expression of halophyte Kosteletzkya virginica seedlings under salt stress at early stage. Genetica137, 189–199. 10.1007/s10709-009-9384-919588254GuptaN.NaikP. K.ChauhanR. S. (2012). Differential transcript profiling through cDNA-AFLP showed complexity of rutin biosynthesis and accumulation in seeds of a nutraceutical food crop (Fagopyrum spp.). BMC Genom. 13, 231. 10.1186/1471-2164-13-23122686486HanleyM. E. (1998). Seedling herbivory, community composition and plant life history traits. Perspect. Plant Ecol.1, 191–205. 10.1078/1433-8319-0005829491176HanleyM. E.FennerM.WhibleyH.DarvillB. (2004). Early plant growth: identifying the end point of the seedling phase. New Phytol. 163, 61–66. 10.1111/j.1469-8137.2004.01094.x33873793HanleyM. E.MayO. C. (2006). Cotyledon damage at the seedling stage affects growth and flowering potential in mature plants. New Phytol.169, 243–250. 10.1111/j.1469-8137.2005.01578.x16411928ImranM.AaqilK. M.ShahzadR.BilalS.KhanM.YunB. W.. (2021). Melatonin ameliorates thermotolerance in soybean seedling through balancing redox homeostasis and modulating antioxidant defense, phytohormones and polyamines biosynthesis. Molecules26, 5116. 10.3390/molecules2617511634500550JayaramanA.PuranikS.RaiN. K.VidapuS.SahuP. P.LataC.. (2008). cDNA-AFLP analysis reveals differential gene expression in response to salt stress in foxtail millet (Setariaitalica L.). Mol. Biotechnol.40, 241–251. 10.1007/s12033-008-9081-418592419JinX.CaoD.WangZ.MaL.TianK.LiuY.. (2019). Genome-wide identification and expression analyses of the LEA protein gene family in tea plant reveal their involvement in seed development and abiotic stress responses. Sci. Rep.9, 14123. 10.1038/s41598-019-50645-831575979JouY.ChiangC. P.YenH. E. (2013). Changes in cellular distribution regulate SKD1 ATPase activity in response to a sudden increase in environmental salinity in halophyte ice plant. Plant Signal. Behav.8, e27433. 10.4161/psb.2743324390077JouY.ChouP. H.HeM.HungY.YenH. E. (2004). Tissue-specific expression and functional complementation of a yeast potassium-uptake mutant by a salt-induced ice plant gene mcSKD1. Plant Mol. Biol.54, 881–893. 10.1007/s11103-004-0335-715604658JucknischkeA.KutscheraU. (1998). The role of the cotyledons and primary leaves during seedling etablishment in sunflower. J. Plant Physiol.153, 700–705.KeighleyJ. L. S. (1984). Studies on the Photosynthetic Role of the Cotyledons in the Intact Soybean (Glycine max) Seedling. Dissertations & Theses-Gradworks.KhanA.AnwarY.HasanM. M.IqbalA.AliM.AlharbyH. F.. (2017). Attenuation of drought stress in Brassica seedlings with exogenous application of Ca2+ and H2O2. Plants6, 20. 10.3390/plants602002028505096KitajimaK.FennerM. (2000). Ecology of seedling regeneration in Seeds the Ecology of Regeneration in Plant Communities, eds FennerM.(London: CABI), 331–359.KongL.GXieY.HuL.SiJ.WangZ. (2017). Excessive nitrogen application dampens antioxidant capacity and grain filling in wheat as revealed by metabolic and physiological analyses. Sci. Rep.7, 43363. 10.1038/srep4336328233811KoteyevaN. K.VoznesenskayaE. V.BerryJ. O.CousinsA. B.an EdwardsG. E. (2016). The unique structural and biochemical development of single cell C4 photosynthesis along longitudinal leaf gradients in Bienertiasinuspersici and Suaeda aralocaspica (Chenopodiaceae). J. Exp. Bot.67, 2587–2601. 10.1093/jxb/erw08226957565KoteyevaN. K.VoznesenskayaE. V.CousinsA. B.EdwardsG. E. (2014). Differentiation of C4 photosynthesis along a leaf developmental gradient in two Cleome species having different forms of Kranz anatomy. J. Exp. Bot.65, 3525–3541. 10.1093/jxb/eru04224550438LamontB. B.GroomP. K. (2002). Green cotyledons of two Hakea species control seedling mass and morphology by supplying mineral nutrients rather than organic compounds. New Phytol.153, 101–110. 10.1046/j.0028-646X.2001.00300.xLangS. R.LiuX. X.MaG.LanQ. Y.WangX. F. (2014). Identification of desiccation tolerance transcripts potentially involved in rape (Brassica napus L.) seeds development and germination. Plant Physiol. Bioch.83, 316–326. 10.1016/j.plaphy.2014.08.00125221920LiC. X.WeiX. W.LanH. Y. (2021). The environmental adaptation strategy of seed germination, and roles of the seed pappus on dispersal and hypocotyl hairs on seedling anchorage in Tamarixramosissima. AoB Plants13, plab065. 10.1093/aobpla/plab06534804468LiM.HuW.RenL.JiaC.LiuJ.MiaoH.. (2020). Identification, expression, and interaction network analyses of the CDPK gene family reveal their involvement in the development, ripening, and abiotic stress response in banana. Biochem. Genet.58, 40–62. 10.1007/s10528-019-09916-231144068LiX.SuQ.ZhengR.LiuG.LuY.BianL.. (2015). ClRTL1 encodes a Chinese Fir RNase III-like protein involved in regulating shoot branching. Int. J. Mol. Sci.16, 25691–25710. 10.3390/ijms16102569126516842LiuF.WangY.DingZ.T.ZhaoL.XiaoJ.WangL.. (2017). Transcriptomic analysis of flower development in tea (Camellia sinensis (L.)). Gene631, 39–51. 10.1016/j.gene.2017.08.01328844668LiuH.SunZ.HuL.LiC.WangX.YueZ.. (2021). Comparative transcriptome analysis of male sterile anthers induced by high temperature in wheat (Triticum aestivum L.). Front. Plant Sci. 12, 727966. 10.3389/fpls.2021.72796634759937LiuJ.GuoW.Q.ShiD.C. (2010). Seed germination, seedling survival, and physiological response of sunflowers under saline and alkaline conditions. Photosynthetica48, 278–286. 10.1007/s11099-010-0034-3LiuY. X.MaimaitijiangT.ZhangJ. H.MaY. L.LanH. Y. (2020). The developmental enhancement of a C4 system with non-typical C4 physiologicalcharacteristics in Salsola ferganica (kranz anatomy), an annual desert halophyte. Front. Plant Sci.11, 152. 10.3389/fpls.2020.0015232210984LohmanK.N.GanS.JohnM.C.AmasinoR.M. (1994). Molecular analysis of natural leaf senescence in Arabidopsis thaliana. Physiol. Plant.92, 322–328. 10.1111/j.1399-3054.1994.tb05343.xLovellP.MooreK. (1971). A comparative study of the role of the cotyledon in seedling development. J. Exp. Bot.22, 153–162. 10.1093/jxb/22.1.153LyndonR. F. (1998). The shoot apical meristem: its growth and development. Q. Rev. Biol.144, 51–54.23110895MaY. L.WangJ.ZhangJ. HZhangS. Y.LiuY. X.LanH.Y. (2018). Seed heteromorphism and effects of light and abiotic stress on germination of a typical annual halophyte Salsola ferganica in cold desert. Front. Plant Sci. 8, 2257. 10.3389/fpls.2017.0225729387073MacoveiA.TutejaN. (2012). microRNAs targeting DEAD-box helicases are involved in salinity stress response in rice (Oryza sativa L.). BMC Plant Biol. 12, 183. 10.1186/1471-2229-12-18323043463MarekL. F.StewartC. R. (1992). Photosynthesis and photorespiration in presenescent, senescent, and rejuvenated soybean cotyledons. Plant Physiol. 98, 694–699. 10.1104/pp.98.2.69416668697McDonaldJ. H. (2014). Handbook of Biological Statistics, 3rd Edn. Baltimore, MA: Sparky House Publishing.MeyerG.A. (1998). Pattern of defoliation and its effect on photosynthesis and growth of goldenrod. Funct. Ecol.12, 270–279. 10.1046/j.1365-2435.1998.00193.xMilbergP.LamontB.B. (1997). Seed/cotyledon size and nutrient content play a major role in early performance of species in nutrient poor soils. New Phytol. 137, 665–672. 10.1046/j.1469-8137.1997.00870.xMoreauM.AzzopardiM.ClémentG.DobrenelT.MarchiveC.RenneC.. (2012). Mutations in the Arabidopsis homolog of LST8/GβL, a partner of the target of Rapamycin kinase, impair plant growth, flowering, and metabolic adaptation to long days. Plant Cell. 24, 463–481. 10.1105/tpc.111.09130622307851MotmainnaM.JuraimiA. S.UddinM. K.AsibN. B.IslamAKMMAhmad-HamdaniM. S.BerahimZ.. (2021). Physiological and Biochemical Responses of Ageratum conyzoides, Oryza sativa f. spontanea (Weedy Rice) and Cyperus iria to Parthenium hysterophorus Methanol Extract. Plants10, 1205. 10.3390/plants1006120534198474NakagawaA.SakamotoS.TakahashiM.MorikawaH.SakamotoA. (2007). The RNAi-mediated silencing of xanthine dehydrogenase impairs growth and fertility and accelerates leaf senescence in transgenic Arabidopsis plants. Plant Cell Physiol.48, 1484–1495. 10.1093/pcp/pcm11917872919NiazianM.Sadat-NooriS. A.TohidfarM.MortazavianS. M.M.SabbatiniP. (2021). Betaine aldehyde dehydrogenase (BADH) vs. flavodoxin (Fld): two important genes for enhancing plants stress tolerance and productivity. Front. Plant Sci.12, 650215. 10.3389/fpls.2021.65021533868350OjedaV.NájeraV.A.GonzálezM.Pérez-RuizJ.M.CejudoF.J. (2017). Photosynthetic activity of cotyledons is critical during post-germinative growth and seedling establishment. Plant Signal. Behav.12, e1347244. 10.1080/15592324.2017.134724428692378PanX.WeltiR.WangX. (2010). Quantitative analysis of major plant hormones in crude plant extracts by high-performance liquid chromatography-mass spectrometry. Nat. Protoc.5, 986–992. 10.1038/nprot.2010.3720448544ParmarK. M.GaikwadS. L.DhakephalkarP. K.KothariR.SinghR. P. (2017). Intriguing interaction of bacteriophage-host association: an understanding in the era of omics. Front. Microbiol.8, 559. 10.3389/fmicb.2017.0055928439260QianX.GongM.J.WangC.X.TianM. (2014). cDNA-AFLP transcriptional profiling reveals genes expressed during flower development in OncidiumMilliongolds. Genet. Mol. Res.13, 6303–6315. 10.4238/2014QiuN.LiuQ.LiJ.ZhangY.WangF.GaoJ. (2017). Physiological and transcriptomic responses of Chinese cabbage (Brassica rapa L. ssp. pekinensis) to salt stress. Int J Mol Sci.18, 1953. 10.3390/ijms1809195328895882RabeloM. C.BangW. Y.NairV.AlvesR. E.Jacobo-VelázquezD. A.SreedharanS.. (2020). UVC light modulates vitamin C and phenolic biosynthesis in acerola fruit: role of increased mitochondria activity and ROS production. Sci. Rep.10, 21972. 10.1038/s41598-020-78948-133319819ReijansM.LascarisR.GroenegerA. O.WittenbergA.WesselinkE.OeverenJ. V.. (2003). Quantitative comparison of cDNA-AFLP, microarrays, and genechip expression data in Saccharomyces cerevisiae. Genomics82, 606–618. 10.1016/S0888-7543(03)00179-414611802Riccio-RengifoC.FinkeJ.RochaC. (2021). Identifying stress responsive genes using overlapping communities in co-expression networks. BMC Bioinform.22, 541. 10.1186/s12859-021-04462-434743699RuffinoA. M. C.RosaM.HilalM.GonzálezJ. A.PradoF. E. (2010). The role of cotyledon metabolism in the establishment of quinoa (Chenopodium quinoa) seedlings growing under salinity. Plant Soil326, 213–224. 10.1007/s11104-009-9999-8SadakM. S.BakryB. A. (2020). Alleviation of drought stress by melatonin foliar treatment on two flax varieties under sandy soil. Physiol. Mol. Biol. Plants26, 907–919. 10.1007/s12298-020-00789-z32377041ŠamecD.LjubejV.Redovnikovi,ćI. R.Fistani,ćS.Salopek-SondiB. (2022). Low temperatures affect the physiological status and phytochemical content of flat leaf kale (Brassica oleracea var. acephala) sprouts. Foods11, 264. 10.3390/foods1103026435159416SaroshB. R.MeijerJ. (2007). Transcriptional profiling by cDNA-AFLP reveals novel insights during methyl jasmonate, wounding and insect attack in Brassica napus. Plant Mol. Biol.64, 425–438. 10.1007/s11103-007-9164-917401749ShazadeeH.KhanN.WangJ.WangC. C.ZengJ. G.HuangZ.Y.. (2019). Identification and expression profiling of protein phosphatases (PP2C) gene family in Gossypium hirsutum L. Int. J. Mol. Sci.20, 1395. 10.3390/ijms2006139530897702ShiR.ChiangV. L. (2005). Facile means for quantifying microRNA expression by real-time PCR. BioTechniques.39, 519–525. 10.2144/00011201016235564ShiS.AnL.MaoJ.AlukoO. O.UllahZ.XuF.. (2021). The CBL-interacting protein kinase NtCIPK23 positively regulates seed germination and early seedling development in tobacco (Nicotiana tabacum L.). Plants10, 323. 10.3390/plants1002032333567573SmithM. E.KoteyevaN. K.VoznesenskayaE. V.OkitaT. W.EdwardsG. E. (2009). Photosynthetic features of non-Kranz type C4 versus Kranz type C4 and C3 species in subfamily Suaedoideae (Chenopodiaceae). Funct. Plant Biol.36, 770–782. 10.1071/FP0912032688687StreetN. R.SjödinA.Bylesj,öM.GustafssonP.TryggJ.JanssonS. (2008). A cross-species transcriptomics approach to identify genes involved in leaf development. BMC Genom.9, 589. 10.1186/1471-2164-9-58919061504SzékelyG.AbrahámE.CséploA.Rig,óG.ZsigmondL.CsiszárJ.. (2008). Duplicated P5CS genes of Arabidopsis play distinct roles in stress regulation and developmental control of proline biosynthesis. Plant J.53, 11–28. 10.1111/j.1365-313X.2007.03318.x17971042TangH.NiuL.WeiJ.ChenX.ChenY. (2019). Phosphorus limitation improved salt tolerance in maize through tissue mass density increase, osmolytes accumulation, and Na+ uptake inhibition. Front. Plant Sci.10, 856. 10.3389/fpls.2019.0085631333699UrbanavičiuteI.BonfiglioliL.PagnottaM. A. (2021). One hundred candidate genes and their roles in drought and salt tolerance in wheat. Int. J. Mol. Sci.22, 6378. 10.3390/ijms2212637834203629UzildayB.OzgurR.YalcinkayaT.TurkanI.SekmenA. H. (2018). Changes in redox regulation during transition from C3 to single cell C4 photosynthesis in Bienertia sinuspersici. J. Plant Physiol. 220, 1–10. 10.1016/j.jplph.2017.10.00629128610UzildayB.TurkanI.OzgurR.SekmenA. H. (2014). Strategies of ROS regulation and antioxidant defense during transition from C3 to C4 photosynthesis in the genus Flaveria under PEG-induced osmotic stress. J. Plant Physiol. 171, 65–75. 10.1016/j.jplph.2013.06.01623920414UzildayB.TurkanI.SekmenA. H.OzgurR.KarakayaH.C. (2012). Comparison of ROS formation and antioxidant enzymes in Cleome gynandra (C4) and Cleome spinosa (C3) under drought stress. Plant Sci. 182, 59–70. 10.1016/j.plantsci.2011.03.01522118616VoznesenskayaE. V.EdwardsG. E.KiiratsO.ArtyushevaE. G.FranceschiV. R. (2003). Development of biochemical specialization and organelle partitioning in the single celled C4 system in leaves of Borszczowiaaralocaspica (Chenopodiaceae). Am. J. Bot.90, 1669–1680. 10.3732/ajb.90.12.166921653343VoznesenskayaE. V.FranceschiV. R.EdwardsG. E. (2004). Light-dependent development of single cell C4 photosynthesis in cotyledons of Borszczowiaaralocaspica (Chenopodiaceae) during transformation from a storage to a photosynthetic organ. Ann. Bot.93, 177–187. 10.1093/aob/mch02614707001VoznesenskayaE. V.FranceschiV. R.KiiratsO.FreitagH.EdwardsG. E. (2001). Kranz anatomy is not essential for terrestrial C4 plant photosynthesis. Nature414, 543–546. 10.1038/3510707311734854WangH.WangS.FanM. M.ZhangS. H.SunL. L.ZhaoZ. Y. (2021a). Metabolomic insights into the browning of the peel of bagging 'Rui Xue' apple fruit. BMC Plant Biol.21, 209. 10.1186/s12870-021-02974-y33964877WangJ.ChengG.WangC.HeZ. Z.LanX. X.ZhangS. Y.. (2017). The bHLH transcription factor CgbHLH001 is a potential interaction partner of CDPK in halophyte Chenopodium glaucum. Sci. Rep.7, 8441. 10.1038/s41598-017-06706-x30573848WangL.HuangZ. Y.BaskinC. C.BaskinJ. M.DongM. (2008). Germination of dimorphic seeds of the desert annual halophyte Suaeda aralocaspica (Chenopodiaceae), a C4 plant without Kranz anatomy. Ann. Bot.102, 757–769. 10.1093/aob/mcn15818772148WangL.MaG. L.WangH. L.ChengC.MuS. Y.QuanW. L.. (2019a). A draft genome assembly of halophyte Suaeda aralocaspica, a plant that performs C4 photosynthesis within individual cells. Gigascience8, giz116. 10.1093/gigascience/giz11631513708WangP.LuS.ZhangX.HydenB.QinL.LiuL.. (2021b). Double NCED isozymes control ABA biosynthesis for ripening and senescent regulation in peach fruits. Plant Sci.304, 110739. 10.1016/j.plantsci.2020.11073933568291WangY.PengX.SalvatoF.WangY.YanX.ZhouZ.. (2019b). Salt-adaptive strategies in oil seed crop Ricinus communis early seedlings (cotyledon vs. true leaf) revealed from proteomics analysis. Ecotoxicol. Environ. Saf.171, 12–25. 10.1016/j.ecoenv.2018.12.04630593996WuJ.ZhaoG.YangY. N.LeW. Q.KhanM. A.ZhangS. L.. (2013). Identification of differentially expressed genes related to coloration in red/green mutant pear (Pyrus communis L.). Tree Genet. Genom.9, 75–83. 10.1007/s11295-012-0534-3XuX.ZhangL.CaoX.LiuL.JiangB.ZhangC.. (2021). Cotyledons facilitate the adaptation of early-maturing soybean varieties to high-latitude long-day environments. Plant Cell Environ. 44, 2551–2564. 10.1111/pce.1412034050544XueR. F.WuJ.ChenM. L.ZhuZ. D.WangL. F.WangX. M.. (2014). Cloning and characterization of a novel secretory root-expressed peroxidase gene from common bean (Phaseolus vulgaris L.) infected with Fusarium oxysporum f. sp. Phaseoli. Mol. Breeding34, 855–870. 10.1007/s11032-014-0080-9XueR. F.WuJ.ZhuZ. D.WangL. F.WangX. M.WangS. M.. (2015). Differentially expressed genes in resistant and susceptible common bean (Phaseolus vulgaris L.) genotypes in response to Fusarium oxysporum f. sp. phaseoli. PLoS ONE10, e0127698. 10.1371/journal.pone.012769826030070YangD. Q.LiuY. J.SunM. L.ZhaoL.WangY. S.ChenX. T.. (2012). Differential gene expression in tea (Camellia sinensis L.) calli with different morphologies and catechin contents. J. Plant Physiol.169, 163–175. 10.1016/j.jplph.2011.08.01522078995YiX.WangZ.LiuC.LiuG.ZhangM. (2015). Acorn cotyledons are larger than their seedlings' need: evidence from artificial cutting experiments. Sci. Rep.5, 8112. 10.1038/srep0811225630843YouL.ZhangJ.LiL.XiaoC.FengX.ChenS.. (2020). Involvement of abscisic acid, ABI5, and PPC2 in plant acclimation to low CO2. J. Exp. Bot.71, 4093–4108. 10.1093/jxb/eraa14832206789ZhangS. T.ZhaoC.LambE. G. (2011). Cotyledon damage affects seed number through final plant size in the annual grassland species Medicago lupulina. Ann. Bot.107, 437–442. 10.1093/aob/mcq25921196450ZhaoQ.SuoJ.ChenS.JinY.MaX.YinZ.. (2016). Na2CO3-responsive mechanisms in halophyte Puccinelliatenuiflora roots revealed by physiological and proteomic analyses. Sci. Rep. 6, 32717. 10.1038/srep3271727596441ZhengW.YangJ. Y.LuanZ. H.WangP.ZhangH. X.ZhouD. W. (2011). Compensatory growth and photosynthetic responses of Pharbitis purpurea seedlings to clipped cotyledon and second leaf. Photosynthetica49, 21–28. 10.1007/s11099-011-0004-4ZhouC. P.QiY. P.YouX.YangL. T.GuoP.YeX.. (2013). Leaf cDNA-AFLP analysis of two citrus species differing in manganese tolerance in response to long-term manganese-toxicity. BMC Genet.14, 621. 10.1186/1471-2164-14-62124034812ZhouJ. C.ZhaoW. W.YinC. H.SongJ.WangB. S.FanJ. L.. (2014). The role of cotyledons in the establishment of Suaeda physophoraseedlings. Plant Biosyst.148, 584–590. 10.1080/11263504.2013.78857433534902