Transcriptomic profiling of genes in matured dimorphic seeds of euhalophyte Suaeda salsa
BMC Genomics volume 18, Article number: 727 (2017)
Suaeda salsa (S. salsa) is a euhalophyte with high economic value. S. salsa can produce dimorphic seeds. Brown seeds are more salt tolerant, can germinate quickly and maintain the fitness of the species under high saline conditions. Black seeds are less salt tolerant, may become part of the seed bank and germinate when soil salinity is reduced. Previous reports have mainly focused on the ecophysiological traits of seed germination and production under saline conditions in this species. However, there is no information available on the molecular characteristics of S. salsa dimorphic seeds.
In the present study, a total of 5825 differentially expressed genes were obtained; and 4648 differentially expressed genes were annotated based on a sequence similarity search, utilizing five public databases by transcriptome analysis. The different expression of these genes may be associated with embryo development, fatty acid, osmotic regulation substances and plant hormones in brown and black seeds. Compared to black seeds, most genes may relate to embryo development, and various genes that encode fatty acid desaturase and are involved in osmotic regulation substance synthesis or transport are upregulated in brown seeds. A large number of differentially expressed genes related to plant hormones were found in brown and black seeds, and their possible roles in regulating seed dormancy/germination were discussed.
Upregulated genes involved in seed development and osmotic regulation substance accumulation may relate to bigger seed size and rapid seed germination in brown seeds, compared to black seeds. Differentially expressed genes of hormones may relate to seed dormancy/germination and the development of brown and black seeds. The transcriptome dataset will serve as a valuable resource to further understand gene expression and functional genomics in S. salsa dimorphic seeds.
At present, irreversible soil salinization has been a global issue, which not only affects the sustainable development of agriculture, but also becomes a potential threat of desertification . Salt stress is harmful to plant growth and yield in crops . Therefore, it is essential to develop and utilize saline-alkali land by developing salt-tolerant crops (salt-tolerant germplasm), thereby improving saline soil .
Halophytes, which are defined as plants that naturally inhabit saline environments and benefit from salt in their growth cycle, could be the best choice to develop saline agriculture [4, 5]. Suaeda salsa L. (S. salsa) is an annual euhalophytic herb with succulent leaves, which can be utilized as food, medicine, forage and bioenergy, and plays an important role in the restoration of contaminated land . The species could be a promising model to understand salt tolerance and develop saline agriculture in the future . Seed dimorphism is a powerful germination strategy in unpredictable environments such as deserts and high salt soil areas . S. salsa produces dimorphic seeds in one plant; that is, brown seeds with a soft, semi-transparent outer seed coat and black seeds with a hard, non-transparent outer seed coat [8, 9]. Brown seeds are bigger, more salt tolerant and can absorb water more quickly than black seeds . It has been suggested that brown seeds could germinate in spring under high saline conditions, and black seeds may germinate in late spring and summer when soil-salt content is low due to rainfall [8, 9]. In addition, black seeds have a relatively long period of dormancy, and light is necessary for black seeds under high saline conditions during germination . By producing dimorphic seeds, S. salsa selects the “hedge-betting” strategy to adapt to saline environments during germination .
Seed germination, which determines whether the plants can successfully establish their progeny, is critical in the life cycle of halophytes in saline environments . Seed germination will be delayed or prohibited under salt stress for both halophytes and non-halophytes in general conditions [9, 13, 14]. Seed germination is a complicated process influenced by a large number of genes and environmental factors . The characteristics of the seed itself perform a pivotal role in seed germination such as seed structures, developmental status of the embryo, compounds that are imported from the mother plant, and some factors that are produced by the embryo itself, including plant hormones . In previous studies, a number of experiments on the germination of S. salsa have been performed to investigate the difference between black seeds and brown seeds. However, studies of molecular biology on the difference of dimorphic seeds themselves in the regulation of dormancy and germination have seldom been reported. Furthermore, there is no report on the functional genomics of S. salsa seeds, especially for differentially expressed genes (DEGs) in dimorphic seeds by transcriptome analysis.
In the present study, the transcriptome of brown and black seeds were analyzed by high-throughput Illumina RNA sequencing (RNA-seq). By comparing the transcriptome of black and brown seeds, large numbers of DEGs were identified. The results of the present study may provide further insight into this complex regulatory network, which may be related to seed development, dormancy and seed germination under saline conditions in S. salsa.
Transcriptome sequencing and de novo assembly
A total of 194.87 million clean reads were obtained for assembly, and further analysis after trimming adapters were performed to filter out low quality reads. The reads from S. salsa were assembled into 11,824,856 contigs, with a mean length of 47.91 bp and a N50 length of 48 bp. By performing pair-end joining and clustering, contig libraries were assembled into 208,588 transcripts with a mean length of 1055.93 bp and N50 lengths of 1751 bp. Then, these transcripts were further clustered into unigenes. A total of 91,753 unigenes were obtained from S. salsa with a mean length of 701.53 bp and a N50 length of 1099 bp. The size distribution of these unigenes is shown in Additional file 1: Figure S1 and Table 1. All raw transcriptome data have been deposited in the National Center for Biotechnology Information (NCBI) SRA database, and the accession number can be found in the “Availability of supporting data”.
Selection and identification of DEGs
A total of 36,318 unigenes (39.58%) were identified with a significance threshold. A total of 5825 genes were found to be differentially expressed in black and brown seeds. Among them, 2919 genes were upregulated and 2906 genes were downregulated in brown seeds, compared to black seeds (Fig. 1).
Functional annotation of DEGs
A total of 4648 DEGs were annotated against five public databases, and allowed the assignment of several functional annotations. DEG sequences were searched against the NCBI non-redundant (Nr) database  and Kyoto Encyclopedia of Genes and Genomes (KEGG) database , and 4641 and 1102 DEGs were annotated, respectively. In order to obtain more detailed protein annotations, the sequences were also blasted against the Swiss-Prot databases , and subsequently, 3410 DEGs were annotated. Furthermore, 1900 DEGs were annotated in the Clusters of Orthologous Groups (COG) databases , and 3659 DEGs were annotated into the Gene Ontology (GO) database .
Functional classification by GO
GO analysis provides a structured, precisely defined, dynamic, controlled vocabulary for describing the roles of genes and gene products in any organism with a biological process, molecular function and cellular component . A total of 3659 DEGs were assigned to 53 level-2 GO terms  according to three gene ontology classes, including 25 for biological process, 14 for cellular components, and 14 for molecular function (Fig. 2). For the cellular component group, the most represented category was cell part, cell and organelle. For biological process, metabolic process and cellular process were highly represented. Regarding molecular function, the category of catalytic activity was the most represented GO term, followed second by the category of binding.
Functional classification by COG
COG classification, which is a tool for the genome-scale analysis of protein function and evolution, was carried out to assign functional information to DEGs . Among the 25 COG categories, the cluster for “general function prediction only” had the largest number, followed by “replication, recombination and repair”, “carbohydrate transport and metabolism”, “transcription”, and “signal transduction mechanisms” (Fig. 3). There were no genes assigned in the “nuclear structure” and “intracellular trafficking, secretion, and vesicular transport” (Fig. 3).
Functional classification by KEGG
A total of 1102 DEGs were assigned to 111 KEGG pathways. Most of these DEGs were mapped to “oxidative phosphorylation”, “ribosome”, “glycolysis/gluconeogenesis”, “starch and sucrose metabolism”, “amino sugar and nucleotide sugar metabolism”, “protein processing in endoplasmic reticulum”, and “plant hormone signal transduction” (Additional file 2: Table S1), which indicate that hormone content and carbohydrate metabolism may lead to different germination characteristics between brown and black seeds.
Genes involved in fatty acid and embryo development
Through transcriptome sequencing, 75 genes may be related to embryo development (Additional file 3: Table S2) and 90 genes related to the synthesis and metabolism of fatty acid were also obtained (Additional file 4: Table S3). Most of these genes, which were involved in fatty acid and embryo development, were also upregulated in brown seeds, compared to black seeds (Fig. 4). A total of 10 DEGs were annotated to “biosynthesis of unsaturated fatty acids” by KEGG analysis (Additional file 5: Table S4, Fig. 5). For the biosynthesis of unsaturated fatty acids, genes encoding 3-oxoacyl-[acyl-carrier protein] reductase, acyl-[acyl-carrier-protein] desaturase and stearoyl-CoA desaturase (Delta-9 desaturase) were downregulated. Genes related to very-long-chain enoyl-CoA reductase, acyl-lipid omega-3 desaturase (Delta-15 desaturase), acyl-lipid omega-6 desaturase (Delta-12 desaturase), and ACAA1 (acetyl-CoA acyltransferase 1) were upregulated (Additional file 5: Table S4, Fig. 5).
Genes involved in betaine, proline, sodium, chloride, potassium and calcium accumulation
Numerous DEGs were annotated to the GO classification, and may be related to betaine, proline, sodium, chloride, potassium, and calcium accumulation, which were 2, 12, 11, 3, 18, and 93, respectively (Fig. 6). Most DEGs such as genes related to sodium ion transport (GO: 0006814), chloride transport (GO: 0006821), voltage-gated chloride channel activity (GO: 0005247), potassium ion transmembrane transport (GO: 0071805), cellular potassium ion homeostasis (GO: 0030007), potassium ion transmembrane transporter activity (GO: 0015079), and inward rectifier potassium channel activity (GO: 0005242) were upregulated in brown seeds, compared to black seeds (Additional file 6: Table S5).
Genes involved in auxin (IAA), abscisic acid (ABA), cytokinins (CTK), salicylic acid (SA), jasmonic acid (JA), brassinosteroid (BR), ethylene (ETH) and gibberellins (GAs)
According to the GO annotations, 367 genes were found, which may be related to endogenous hormones. These genes may participate in the hormone biosynthetic process, metabolic or process signaling pathway, and polar transport from all DEGs (Additional file 7: Table S6). Among these, 283 genes were upregulated and 84 genes were downregulated in brown seeds, compared to black seeds (Fig. 7). Most of these genes were related to abscisic acid and auxin (Fig. 7). For example, 38 genes were annotated in “plant hormone signal transduction” by KEGG analysis (Fig. 8, Additional file 8: Table S7), and 11 genes were annotated to the auxin signal transduction process, in which most of genes were upregulated (the proportion was 81.82%). These genes were auxin influx carrier (AUX1), transport inhibitor response-1 (TIR1), auxin responsive GH3 gene family (GH3), and SAUR family protein (SAUR). Furthermore, AUX1, TIR1 and GH3 were upregulated. There were five genes (three genes were upregulated and two genes were downregulated) encoding SAUR. In the cytokinin signal transduction process, A-ARR (two-component response regulator ARR-A family: ARR4 and ARR9) was upregulated, and B-ARR (two-component response regulator ARR-B family: ARR11) was downregulated. In the gibberellin signal transduction process, GID1 (gibberellin-insensitive dwarf 1: GID1C) was downregulated. In the abscisic acid signal transduction process, SNRK2 (serine/threonine-protein kinase SRK2) was upregulated and ABF (ABA responsive element binding factor) was downregulated. There change to were three genes (two genes were upregulated and one gene was downregulated) encoding PYR/PYL (abscisic acid receptor PYR/PYL family), and two genes (one gene was upregulated and one gene was downregulated) encoding PP2C (protein phosphatase 2C). In the ethylene signal transduction process, ethylene receptor (ETR) was upregulated, while mitogen-activated protein kinase 6 (MPK6) and ethylene-insensitive protein 3 (EIN3) were downregulated. In the brassinosteroid signal transduction process, BRI1 kinase inhibitor 1 (BKI1), BR-signaling kinase (BSK), protein brassinosteroid insensitive 2 (BIN2), BZR1/2 (brassinosteroid resistant 1/2: BZR1 and BZR2/BES1), and cyclin D3 (CYCD3) were upregulated. There are four genes (two genes were upregulated and two genes were downregulated) encoding TGA (transcription factor TGA) in the salicylic acid signal transduction process. There were no DEGs in the jasmonic acid signal transduction process (Fig. 8).
Contents of endogenous hormones
Concentration of Ca2+
The concentration of Ca2+ in brown seeds was much higher than that in black seeds; that is, the Ca2+ concentration in brown seeds was 5.3 times of that in black seeds (Fig. 10).
qRT-PCR validation of RNA-Seq results
Based on the results mentioned above, 20 pairs of primers were designed (Additional file 9: Table S8), and qRT-PCR was applied to validate the correlation against RNA-Seq. Results were basically identical with the sequencing results (R 2 = 0.97, Fig. 11), which indicates that the transcriptome sequencing results are reliable.
In the present study, a large number of differentially expressed genes were identified between brown and black seeds. Brown seeds of S. salsa were much heavier (bigger) and better developed than black seeds . DEGs, according to the GO annotation of the transcriptome sequencing, may participate in the biological process of embryo development (the proportion of upregulated genes in brown seeds compared to black seeds was 61.33%); that is, they may play a positive role in embryo development. Seeds of S. salsa can produce oil, which contains several kinds of fatty acid with sufficient unsaturated fatty acids, especially linoleic acid . In the biosynthesis of the unsaturated fatty acid signal transduction process, fatty acid desaturase, which almost exists in all organisms, performed a vital role . Delta-9 desaturase can convert the saturated fatty acids into oleic acid by introducing the first double bond, and the delta-12 desaturase catalyze oleic acid into linoleic acid, which may further be converted by delta-15 desaturase into a-linolenic acid . Genes that encode delta-12 desaturase had a relatively higher expression in brown seeds than in black seeds. It is possible that better developed embryos and more unsaturated fatty acids can provide more nutrition during rapid seed germination. This deduction needs to be further verified.
In a saline environment, seed germination was inhibited to a certain extent, and was affected by specific-ion toxicity or/and osmotic stress [26,27,28]. High concentration of inorganic ions such as Na+, K+ and Cl− can lower water potential and help seeds absorb water from saline soils during germination [29, 30]. Liu et al. found that the concentration of Na+, K+, Cl− and total ions of brown seeds was much higher than those of black seeds . In addition, the concentration of Ca2+ revealed a similar tendency (Fig. 10). These present results of transcriptome sequencing indicate that most genes related to Na+, K+, Ca2+ and Cl− accumulation were upregulated in brown seeds compared to black seeds (Fig. 6). Vacuolar cation/proton exchanger 3 (CAX3) is localized in the tonoplast, and is involved in ion vacuolar compartmentalization [32,33,34]. Punshon et al. found that the CAX3 transcript was abundant in dry seeds of Arabidopsis . CAX3 is strongly expressed during seed germination [36, 37]. In the present study, CAX3 was upregulated in brown seeds compared to black seeds of S. salsa (Additional file 6: Table S5). Therefore, CAX3 may play an important role in nutrient acquisition in brown seeds. In addition, the high expression of CAX3 may be contributed to the higher salt resistance of brown seeds compared to black seeds in S. salsa. Proteins of the CLC (chloride channel) family, which participate in Cl− balance and prevent Cl− toxicity under salt stress, has a potential physiology function in the salt tolerance of plants [38, 39]. In the present study, all genes related to chloride transport and voltage-gated chloride channel activity (CLC-c, CLC-f, and CLC-g) were upregulated in brown seeds, compared to black seeds of S. salsa (Additional file 6: Table S5). Genes that encode putative potassium transporter 12, two-pore potassium channel 1, and 1-aminocyclopropane-1-carboxylate are significantly upregulated, while the expression of genes that encode probable potassium transporter 13, potassium channel AKT1, and two-pore potassium channel 3 exhibited an opposite trend (Additional file 6: Table S5). The expression of genes may relate to Na+, K+, Ca2+ and Cl− transport carriers or channels that can help regulate salt ions in and out, and reduce salt ion poison under salt stress . It is possible that the interaction of these genes can maintain ion homeostasis in S. salsa dimorphic seeds during germination under salt stress.
In addition to the large number of inorganic ions, some organic osmotic regulation substances were also synthesized to lower the water potential of plants under salt stress [41,42,43]. Proline and betaine, which can adjust osmotic potential and stabilize biological macromolecules, are two important osmotic regulation substances [41, 44]. Compared to black seeds, genes that could encode betaine aldehyde dehydrogenase (BADH) were upregulated (Additional file 6: Table S5). Li et al. found that overexpressing the BADH gene, which was cloned from Suaeda liaotungensis in tobacco, could improve the salt tolerance of transgenic tobacco by producing betaine . Therefore, the upregulation of BADH (Additional file 6: Table S5) may help brown seeds accumulate more betaine than that of black seeds. Yokoishi and Tanimoto found that an increase in betaine appeared to secure seed germination under salt stress in halophyte Suaeda japonica . Genes (fis1) that may encode aldehyde dehydrogenase, which are most likely the P5CDH (delta-1-pyrroline-5-carboxylate dehydrogenase) enzyme related to the proline catabolic process to glutamate, were upregulated in brown seeds, compared to black seeds . Poljakoff-Mayber et al. found that proline concentration was low in dry Kosteletzkya virginica seeds, and proline content increased when seeds were germinated under salt stress . It is possible that proline keeps the developing embryonic axis in a resting state, and the synthesis of proline may be an adaptive strategy of the plant to prevent germination under stressful environments . Therefore, the upregulation of P5CDH may not be related to the high ability of osmotic adjustment in brown seeds, but plays a protective role in brown seed germination under saline conditions.
Plant hormones have important effects on plant development with vanishingly low concentrations, as well as in seed dormancy and germination [50, 51]. Germination is crucial in the whole life cycle of plants, especially for halophytes in variable saline environments [12, 52]. Seed dormancy is a mechanism that can prohibit seed germination in unfavorable conditions . There were large numbers of DEGs probably related to hormone (GAs, ABA, CTK, IAA, ETH, SA, JA, and BR) synthesis, transport and metabolism in the transcriptome sequencing of S. salsa dimorphic seeds (Additional file 7: Table S6). Except for genes for GAs, most genes may be related to other hormones in brown seeds were upregulated, compared to black seeds (Fig. 7). Correspondingly, the content of these hormones in brown seeds were higher than in black seeds, except for ZR and BR (Fig. 9). It is possible that the interactions between these hormones (hormone cross-talk) play profound effects on seed germination. For example, ABA can inhibit germination and induce dormancy, but GAs can stimulate germination . Ethylene promotes seed germination possibly through the antagonism of ABA signaling . JA and SA play a negative regulation of seed germination [54, 55]. It is possible that JA and ABA play a synergistic role during seed germination, but both synergism and antagonism can occur between the JA and ethylene pathways . The auxin influence on seed germination is based on GAs, ABA and ETH regulations . The signals of BR decrease the sensitivity of ABA, and promote seed germination . In a controlled experiment, ZR and IAA content was much higher in brown seeds than in black seeds of S. salsa, which indicate that ZR and IAA may take important functions as growth regulators, and may promote the development of brown seeds, compared to black seeds . CTK and IAA can also increase the salt tolerance of seeds in certain plant species . Therefore, those upregulated genes related to CTK and IAA may promote the development of brown seeds. As a result, brown seeds are much bigger and have a higher ability of salt tolerance than black seeds . Seed germination/dormancy appears to be controlled by the GA/ABA ratio rather than by absolute hormone levels [22, 61]. Most of the genes involved in the gibberellins biosynthetic process (GO: 0009686) were upregulated (Additional file 7: Table S6). Genes that encode GA 20-oxidase (c241829.graph_c0, c223973.graph_c1) (Additional file 7: Table S6), which is a key enzyme that catalyses the last three steps related to GA biosynthesis, were also highly upregulated . It is possible that the upregulation of these genes positively accumulates GAs and further affects the germination/dormancy of brown seeds. ABFs, which belong to a distinct subfamily of bZIP proteins, are induced by ABA and various stress treatments . SnRK2s can phosphorylate ABFs to activate ABA-responsive gene expression . In these present results, ABF was downregulated in brown seeds compared to black seeds (Fig. 8). Therefore, ABFs may increase the degree of dormancy in black seeds of S. salsa. The vital role of ethylene is performed in the dormancy release and germination to non-dormant seeds . ETR is a ethylene receptor that negatively regulates ethylene signaling . The germination of ethylene-insensitive mutant (etr mutant) was inhibited at the same conditions compared with wild-type seeds of Arabidopsis [67, 68]. These present results indicate that the upregulation of ETR in brown seeds (Fig. 8) is probably one of reasons for the quick and high germination of brown seeds. BR is considered to have performed a vital role in antagonizing seed dormancy and stimulating germination . In the present study, many genes related to BR were upregulated in brown seeds, compared with black seeds (Fig. 7, Additional file 7: Table S6), which indicate that BR may also play a role in the quick germination of brown seeds in S. salsa.
In conclusion, the present study reported a transcriptome view of the complexity of gene expression in matured dimorphic seeds of S. salsa. This can help to further investigate the adaptation of seed dimorphic plants to variable saline environments. Further genetic and biochemical analysis should be performed to detect the functions of key genes related to embryo development, inorganic and organic substance accumulation, dormancy and germination in S. salsa dimorphic seeds.
Plant material and RNA extraction
Approximately 50 plants of Suaeda salsa L. from saline inland (N37°20′, E118°36′) were randomly collected in the Yellow River Delta in Shandong province, China, when the seeds became completely mature. After seeds were air-dried, the brown and black seeds were separated; and the seeds were immediately frozen in liquid nitrogen and stored at −80 °C. Total RNA was extracted from pooled samples of brown or black seeds using a total RNA kit (E.Z.N.A.TM Plant RNA Kit R6827-01; OMEGA, Beijing, China), according to manufacturer’s instructions. Then, these were examined using a NanoDrop ND1000 spectrophotometer (NanoDrop Technologies, Wilmington, DE, USA) and characterized on 1% agarose gel electrophoresis. The standards applied were OD260/OD230 ≥ 1.0, and 1.8 ≤ OD260/OD280 ≤ 2.2. The RNA integrity number (RIN) values (>8.0) of these samples were assessed using an Agilent 2100 Bioanalyzer (Santa Clara, CA, USA). At least 20 μg of RNA was pooled in an equimolar fashion from each of the samples. Two replicates of black seeds (T1 and T2) or brown seeds (T3 and T4) were performed in this study, and each sample had 0.5 g of seeds.
Illumina sequencing and de novo assembly
In this experiment, mRNA was enriched and purified with oligo (dT)-rich magnetic beads, and the cDNA libraries were constructed. We tested the library insert fragment size in 1.8% agarose gel electrophoresis. Then, QPCR was performed to enrich the purified cDNA template. Finally, the libraries were sequenced using Illumina HiSeq™2500 with a paired-end strategy at the Illumina Sequencing Services of Biomarker Technologies Co, LTD. (Beijing, China). The original reads obtained from the sequencing data was evaluated and filtered. The Trinity software (http://trinityrnaseq.sourceforge.net/) was adopted to assemble the high quality data after these were filtered, which did not include the Gap in the final length of the constructed transcription. Then, non-redundant unigenes were generated by de novo assembly.
Gene expression analysis
RPKM (Reads per Kilobase per Million mapped reads) values  were used to estimate all unigenes of the gene expression in the brown and black seeds of S. salsa. DEGs in the brown and black seeds of S. salsa were filtered using the DESeq software . Differentially expressed genes were selected on the basis of a threshold of false discovery rate (FDR) <0.01 and an absolute value of log2 Fold Change (log2FC) >1.
Determination of hormone contents
The contents of endogenous IAA, ABA, ZR, JA, BR and GAs in black and brown seeds were analyzed by enzyme linked immunosorbent assay (ELISA) at the College of Crop Sciences, China Agricultural University [59, 72, 73]. Three replicates were used for brown and black seeds.
Measurement of Ca2+ concentration
The sample (0.3 g of dry black or brown seeds for each sample) was prepared by microwave digestion in HNO3-H2O2, and ultrapure water was added to obtain a final volume of 25 mL . Then, these were analyzed using an inductively coupled plasma (ICP-OES) spectrometer (Optima 7300 DV, Perkin-Elmer, USA). Three replicates were used for brown and black seeds.
In order to test the expression of the target gene transcript, a total of 19 DEGs were randomly selected for real-time PCR analysis. All primers were designed using the Beacon Designer software (version 7.9). The housekeeping gene actin (GenBank ID: EU429457) was used as an internal standard . Real-time PCR was carried out in a fluorometric thermal cycler (Bio-Rad CFX96TM Real-time PCR System, California, USA) using SYBR Premix Ex TaqTM (TaKaRa) in 20-μL reaction mixtures containing 10 μL of 2 × SYBR Premix Ex TaqTM (TaKaRa), 50 ng of cDNA and 0.8 μM of primer pairs (actin and target gene). The relative expression of each gene was calculated by 2− ΔΔCT method . Three replicates were performed.
Data for concentrations of endogenous hormones and Ca2+ were analyzed by one-way ANOVA. SAS™ software  was used for statistical analysis.
Two-component response regulator ARR-A family
ABA responsive element binding factor
Acetyl-CoA acyltransferase 1
Auxin influx carrier
Betaine aldehyde dehydrogenase
Two-component response regulator ARR-B family
Protein brassinosteroid insensitive 2
BRI1 kinase inhibitor 1
Brassinosteroid resistant 1/2: BZR1 and BZR2/BES1
Vacuolar cation/proton exchanger 3
Clusters of orthologous groups
Differentially expressed genes
Ethylene-insensitive protein 3
False discovery rate
Auxin responsive GH3 gene family
Gibberellin-insensitive dwarf 1
Kyoto encyclopedia of genes and genomes
Mitogen-activated protein kinase 6
Protein phosphatase 2C
Abscisic acid receptor PYR/PYL family
RNA integrity number
Reads per KB per million
SAUR family protein
Serine/threonine-protein kinase SRK2
Transcription factor TGA
Transport inhibitor response 1
Munns R, Tester M. Mechanisms of salinity tolerance. Annu Rev Plant Biol. 2008;59:651–81.
Paul D, Lade H. Plant-growth-promoting rhizobacteria to improve crop growth in saline soils: a review. Agron Sustain Dev. 2014;34(4):737–52.
Rozema J, Flowers T. Crops for a salinized world. Science. 2008;322(5907):1478–80.
Shabala S. Learning from halophytes: physiological basis and strategies to improve abiotic stress tolerance in crops. Ann Bot. 2013;112(7):1209–21.
Flowers TJ, Colmer TD. Plant salt tolerance: adaptations in halophytes. Ann Bot. 2015;115(3):327–31.
Song J, Wang BS. Using euhalophytes to understand salt tolerance and to develop saline agriculture: Suaeda salsa as a promising model. Ann Bot. 2015;115(3):541–53.
Wang L, Dong M, Huang ZY. Review of research on seed heteromorphism and its ecological significance. Chin J Plant Ecol. 2010;34(5):578–90.
Li W, Liu X, Khan MA, Yamaguchi S. The effect of plant growth regulators, nitric oxide, nitrate, nitrite and light on the germination of dimorphic seeds of Suaeda salsa under saline conditions. J Plant Res. 2005;118(3):207–14.
Song J, Fan H, Zhao YY, Jia YH, Du XH, Wang BS. Effect of salinity on germination, seedling emergence, seedling growth and ion accumulation of a euhalophyte Suaeda salsa in an intertidal zone and on saline inland. Aquat Bot. 2008;88(4):331–7.
Li WQ, An P, Liu XJ. Effect of storage, stratification, temperature and gibberellins on germination of dimorphic seeds of Suaeda salsa under saline conditions. Seed Sci Technol. 2008;36(1):122–32.
Kranner I, Seal CE. Salt stress, signalling and redox control in seeds. Funct Plant Biol. 2013;40(9):848–59.
Song J, Shi WW, Liu RR, Xu YG, Sui N, Zhou JC, et al. The role of the seed coat in adaptation of dimorphic seeds of the euhalophyte Suaeda salsa to salinity. Plant Spec Biol. 2017;32:107–14.
Tessier M, Gloaguen JC, Lefeuvre JC. Factors affecting the population dynamics of Suaeda maritime at initial stages of development. Plant Ecol. 2000;147(2):193–203.
Khajeh-Hosseini M, Powell AA, Bingham IJ. Experimental approach influences soybean (Glycine max L.) seed and seedling response to salinity. Seed Sci Technol. 2005;33(3):629–38.
Bentsink L, Koornneef M. Seed dormancy and germination. Arabidopsis Book. 2008;6(6):e0119.
Koornneef M, Bentsink L, Hilhorst H. Seed dormancy and germination. Curr Opin Plant Biol. 2002;5(1):33–6.
Pruitt KD, Tatusova T, Maglott DR. NCBI Reference Sequence (RefSeq): a curated non-redundant sequence database of genomes, transcripts and proteins. Nucleic Acids Res. 2005;33(suppl 1):D501–4.
Kanehisa M, Goto S, Kawashima S, Okuno Y, Hattori M. The KEGG resource for deciphering the genome. Nucleic Acids Res. 2004;32(suppl 1):D277–80.
Apweiler R, Bairoch A, Wu CH, Barker WC, Boeckmann B, Ferro S, et al. UniProt: the universal protein knowledgebase. Nucleic Acids Res. 2004;32(suppl 1):D115–9.
Tatusov RL, Galperin MY, Natale DA, Koonin EV. The COG database: a tool for genome-scale analysis of protein functions and evolution. Nucleic Acids Res. 2000;28(1):33–6.
Ashburner M, Ball CA, Blake JA, Botstein D, Butler H, Cherry JM, et al. Gene ontology: tool for the unification of biology. Nat Genet. 2000;25(1):25–9.
Conesa A, Götz S, García-Gómez JM, Terol J, Talón M, Robles M. Blast2GO: A universal tool for annotation, visualization and analysis in functional genomics research. Bioinformatics (Oxford, England). 2005;21(18):3674–6.
Li HS, Fan YX. Extraction and characteristics analysis of Suaeda salsa seed oil. China Oils Fats. 2010;35(1):74–6.
Zeng SS, Jiang LM, Yuan DJ. Research and dvelopment of fatty acid desaturase. Life Sci. 2008;20(5):816–20.
Certik M, Shimizu S. Biosynthesis and regulation of microbial polyunsaturated fatty acid production. J Biosci Bioeng. 1999;87(1):1–14.
Gulzar S, Khan MA. Seed germination of a halophytic grass Aeluropus lagopoides. Ann Bot. 2001;87(3):319–24.
Gul B, Ansari R, Flowers TJ, Khan MA. Germination strategies of halophyte seeds under salinity. Environ Exp Bot. 2013;92:4–18.
Hameed A, Rasheed A, Gul B, Khan MA. Salinity inhibits seed germination of perennial halophytes Limonium stocksii and Suaeda fruticosa by reducing water uptake and ascorbate dependent antioxidant system. Environ Exp Bot. 2014;107:32–8.
Zhao KF, Fan H, Zhou S, Song J. Study on the salt and drought tolerance of Suaeda salsa and Kalanchoe claigremontiana under iso-osmotic salt and water stress. Plant Sci. 2003;165(4):837–44.
Song J, Feng G, Tian CY, Zhang FS. Strategies for adaptation of Suaeda physophora, Haloxylon ammodendron and Haloxylon persicum to a saline environment during seed-germination stage. Ann Bot. 2005;96(3):399–405.
Liu Y, Zhou JC, Zhang XD, Li X, Fan H, Wang BS, et al. Adaptations of dimorphic seeds and seedlings of Suaeda salsa to saline environments. Acta Ecol Sin. 2013;33(17):5162–8.
Shigaki T, Cheng NH, Pittman JK, Hirschi K. Structural determinants of Ca2+ transport in the Arabidopsis H+/Ca2+ antiporter CAX1. J Biol Chem. 2001;276(46):43152–9.
Cheng NH, Pittman JK, Shigaki T, Lachmansingh J, LeClere S, Lahner B, et al. Functional association of Arabidopsis CAX1 and CAX3 is required for normal growth and ion homeostasis. Plant Physiol. 2005;138(4):2048–60.
Martinoia E, Maeshima M, Neuhaus HE. Vacuolar transporters and their essential role in plant metabolism. J Exp Bot. 2007;58(1):83–102.
Punshon T, Hirschi K, Yang J, Lanzirotti A, Lai B, Guerinot ML. The role of CAX1 and CAX3 in elemental distribution and abundance in Arabidopsis seed. Plant Physiol. 2012;158(1):352–62.
Zhao J, Barkla BJ, Marshall J, Pittman JK, Hirschi KD. The Arabidopsis cax3 mutants display altered salt tolerance, pH sensitivity and reduced plasma membrane H+-ATPase activity. Planta. 2007;227(3):659–69.
Zhao J, Shigaki T, Mei H, Guo YQ, Cheng NH, Hirschi KD. Interaction between Arabidopsis Ca2+/H+ exchangers CAX1 and CAX3. J Biol Chem. 2009;284(7):4605–15.
De Angeli A, Monachello D, Ephritikhine G, Frachisse JM, Thomine S, Gambale F, et al. CLC-mediated anion transport in plant cells. Philos T R Soc B. 2009;364(1514):195–201.
Zifarelli G, Pusch M. CLC transport proteins in plants. FEBS Lett. 2010;584(10):2122–7.
Zhu JK. Salt and drought stress signal transduction in plants. Annu Rev Plant Biol. 2002;53(1):247–73.
Flowers TJ, Colmer TD. Salinity tolerance in halophytes. New Phytol. 2008;179(4):945–63.
Guo JR, Suo SS, Wang BS. Sodium chloride improves seed vigour of the euhalophyte Suaeda salsa. Seed Sci Res. 2015;25(3):335–44.
Zhou JC, Fu TT, Sui N, Guo JR, Feng G, Fan JL, et al. The role of salinity in seed maturation of the euhalophyte Suaeda salsa. Plant Biosyst. 2014;150(1):83–90.
Yeo A. Molecular biology of salt tolerance in the context of whole-plant physiology. J Exp Bot. 1998;49(323):915–29.
Li QL, Gao XR, Yu XH, Wang XZ, An LJ. Molecular cloning and characterization of betaine aldehyde dehydrogenase gene from Suaeda liaotungensis and its use in improved tolerance to salinity in transgenic tobacco. Biotechnol Lett. 2003;25(17):1431–6.
Yokoishi T, Tanimoto S. Seed germination of the halophyte Suaeda japonica under salt stress. J Plant Res. 1994;107(4):385–8.
Mitchell HJ, Ayliffe MA, Rashid KY, Pryor AJ. A rust-inducible gene from flax (fis1) is involved in proline catabolism. Planta. 2006;223(2):213–22.
Poljakoff-Mayber A, Somers GF, Werker E, Gallagher JL. Seeds of Kosteletzkya virginica (Malvaceae): their structure, germination, and salt tolerance. II. Germination and salt tolerance. Am J Bot. 1994;81(1):54–9.
Thakur M, Sharma AD. Salt-stress-induced proline accumulation in germinating embryos: evidence suggesting a role of proline in seed germination. J Arid Environ. 2005;62(3):517–23.
Kucera B, Cohn MA, Leubner-Metzger G. Plant hormone interactions during seed dormancy release and germination. Seed Sci Res. 2005;15(4):281–307.
Miransari M, Smith DL. Plant hormones and seed germination. Environ Exp Bot. 2014;99:110–21.
Wei Y, Dong M, Huang ZY, Tan DY. Factors influencing seed germination of Salsola affinis (Chenopodiaceae), a dominant annual halophyte inhabiting the deserts of Xinjiang, China. Flora. 2008;203(2):134–40.
Finkelstein R, Reeves W, Ariizumi T, Steber C. Molecular aspects of seed dormancy. Annu Rew Plant Biol. 2008;59(1):387–415.
Rajjou L, Gallardo K, Job C, Job D. Proteome analysis for the study of developmental processes in plants. Plant Proteomics Ann Plant Rev. 2008;28:151–84.
Dave A, Hernandez ML, He Z, Andriotis VM, Vaistij FE, Larson TR, et al. 12-oxo-phytodienoic acid accumulation during seed development represses seed germination in Arabidopsis. Plant Cell. 2011;23(2):583–99.
Linkies A, Leubner-Metzger G. Beyond gibberellins and abscisic acid: how ethylene and jasmonates control seed germination. Plant Cell Rep. 2012;31(2):253–70.
Carrera E, Holman T, Medhurst A, Dietrich D, Footitt S, Theodoulou FL, et al. Seed after-ripening is a discrete developmental pathway associated with specific gene networks in Arabidopsis. Plant J. 2008;53(2):214–24.
Steber CM, McCourt P. A role for brassinosteroids in germination in Arabidopsis. Plant Physiol. 2001;125(2):763–9.
Wang FX, Xu YG, Wang BS, Shi WW, Liu RR, Feng G, et al. Salinity affects production and salt tolerance of dimorphic seeds of Suaeda salsa. Plant Physiol Bioch. 2015;95:41–8.
Park J, Kim YS, Kim SG, Jung JH, Woo JC, Park CM. Integration of auxin and salt signals by the NAC transcription factor NTM2 during seed germination in Arabidopsis. Plant Physiol. 2011;156(2):537–49.
Shu K, Liu XD, Xie Q, He ZH. Two faces of one seed: hormonal regulation of dormancy and germination. Mol Plant. 2016;9(1):34–45.
Calvo AP, Nicolas C, Nicolas G, Rodríguez D. Evidence of a cross-talk regulation of a GA 20-oxidase (FsGA20ox1) by gibberellins and ethylene during the breaking of dormancy in Fagus sylvatica seeds. Physiol Plant. 2004;120(4):623–30.
Choi HI, Hong JH, Ha JO, Kang JY, Kim SY. ABFs, a family of ABA-responsive element binding factors. J Biol Chem. 2000;275(3):1723–30.
Nakashima K, Yamaguchi-Shinozaki K. ABA signaling in stress-response and seed development. Plant Cell Rep. 2013;32(7):959–70.
Kepczynski J, Kepczynska E. Ethylene in seed dormancy and germination. Physiol Plant. 1997;101(4):720–6.
Hua J, Meyerowitz EM. Ethylene responses are negatively regulated by a receptor gene family in Arabidopsis thaliana. Cell. 1998;94(2):261–71.
Bleecker AB, Estelle MA, Somerville C, Kende H. Insensitivity to ethylene conferred by a dominant mutation in Arabidopsis thaliana. Science. 1988;241(4869):1086–9.
Lashbrook CC, Tieman DM, Klee HJ. Differential regulation of the tomato ETR gene family throughout plant development. Plant J. 1998;15(2):243–52.
Mortazavi A, Williams BA, McCue K, Schaeffer L, Wold B. Mapping and quantifying mammalian transcriptomes by RNA-Seq. Nat Methods. 2008;5(7):621–8.
Anders S, Huber W. Differential expression analysis for sequence count data. Genome Boil. 2010;11(10):R106.
Altschul SF, Madden TL, Schaffer AA, Zhang JH, Zhang Z, Miller W, et al. Gapped BLAST and PSI-BLAST: a new generation of protein database search programs. Nucleic Acids Res. 1997;25(17):3389–402.
Yang YM, Xu CN, Wang BM, Jia JZ. Effects of plant growth regulators on secondary wall thickening of cotton fibres. Plant Growth Regul. 2001;35(3):233–7.
Xu YG, Liu RR, Sui N, Shi WW, Wang L, Tian CY, et al. Changes in endogenous hormones and seed coat phenolics during seed storage of two Suaeda salsa populations. Aust J Bot. 2016;64(4):325–32.
Zhang SB, Ji XW, Liu CP. Study on multi-elements in traditional chinese medicines Ophiopogon japonicus and Lotus seeds by the ICP-AES with microwave digestion. Adv Mater Res Trans Tech Publ. 2012;535:1126–32.
Ma Q, Zhou XR, Wu GQ, Wang SM. Cloning and sequence analysis of actin gene fragment from the halophyte Suaeda glauca. Biotechnology. 2009;19(1):1–3.
Livak KJ, Schmittgen TD. Analysis of relative gene expression data using real-time quantitative PCR and the 2− ΔΔCT method. Methods. 2001;25(4):402–8.
SAS Institute Inc. SAS/STAT user's guide. Cary: SAS Institute; 1989.
The financial support from National Natural Science Research Foundation of China (31370420, 31570392) is greatly appreciated.
The financial of this research was supported by National Natural Science Research Foundation of China (31370420, 31570392).
Availability of data and materials
The datasets supporting the results of this article are included within the article and its additional files. All raw transcriptome data have been deposited in the National Center for Biotechnology Information (NCBI) SRA database with the accession number of SRR5877206, and SRR5877205 for black seeds; SRR5877204, and SRR5877203 for brown seeds. Access to the data is available upon publication at http://www.ncbi.nlm.nih.gov/bioproject/395283.
Ethics approval and consent to participate
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Length distribution of unigenes identified in the transcriptome of S. salsa dimorphic seeds. (DOC 141 kb)
Functional classification by KEGG pathways in S. salsa dimorphic seeds. (XLS 20 kb)
DEGs related to embryo development according to GO function annotation. (XLS 27 kb)
DEGs related to fatty acid according to GO function annotation. (XLS 31 kb)
DEGs in the biosynthesis of the unsaturated fatty acid pathway. (XLS 10 kb)
DEGs related to osmotic regulation substances according to GO function annotation. (XLS 57 kb)
DEGs related to plant hormones according to GO function annotation. (XLS 118 kb)
DEGs in the plant hormone signal transduction pathway. (XLS 23 kb)
Primer pairs for real-time quantitative PCR. (XLS 17 kb)
About this article
Cite this article
Xu, Y., Zhao, Y., Duan, H. et al. Transcriptomic profiling of genes in matured dimorphic seeds of euhalophyte Suaeda salsa . BMC Genomics 18, 727 (2017). https://doi.org/10.1186/s12864-017-4104-9