Qian Li1,2,3†
Jiaxing Song
Yongzhen Pang
Bo Zhang- 1Key Laboratory of Grassland Resources and Ecology of Western Arid Region, Ministry of Education, College of Grassland Science, Xinjiang Agricultural University, Urumqi, China
- 2Key Laboratory of Grassland Resources and Ecology of Xinjiang, College of Grassland Science, Xinjiang Agricultural University, Urumqi, China
- 3Institute of Animal Science, Chinese Academy of Agricultural Sciences, Beijing, China
- 4College of Grassland Agriculture, Northwest A&F University, Yangling, China
- 5School of Agriculture Food and Wine, The University of Adelaide, Urrbrae, SA, Australia
Bromus inermis L. (commonly known as smooth bromegrass) is a grass species with high nutritional value, great palatability, cold tolerance, and grazing resistance, which has been widely cultivated for pasture and sand fixation in northern and northwestern China. Salt stress is a main environmental factor limiting growth and production of smooth bromegrass. In this study, we performed PacBio Iso-Seq to construct the first full-length transcriptome database for smooth bromegrass under 300 mM NaCl treatment at different time points. Third-generation full-length transcriptome sequencing yielded 19.67 G polymerase read bases, which were assembled into 355,836 full-length transcripts with an average length of 2,542 bp. A total of 116,578 differentially expressed genes were obtained by comparing the results of third-generation sequencing and second-generation sequencing. GO and KEGG enrichment analyses revealed that multiple pathways were differently activated in leaves and roots. In particular, a number of genes participating in the molecular network of plant signal perception, signal transduction, transcription regulation, antioxidant defense, and ion regulation were affected by NaCl treatment. In particular, the CBL-CIPK, MAPK, ABA signaling network, and SOS core regulatory pathways of Ca2+ signal transduction were activated to regulate salt stress response. In addition, the expression patterns of 10 salt-responsive genes were validated by quantitative real-time PCR, which were consistent with those detected by RNA-Seq. Our results reveal the molecular regulation of smooth bromegrass in response to salt stress, which are important for further investigation of critical salt responsive genes and molecular breeding of salt-tolerant smooth bromegrass.
Introduction
The salt content of the soil surface is high in arid and semi-arid areas due to low precipitation, high evaporation, and inadequate water resource management (Rengasamy, 2010). Soil salinity exists in almost every irrigated area in the world, as well as in non-irrigated farmland and pasture (Ivushkin et al., 2019; Kearl et al., 2019). High content of salt in the soil causes dehydration of plant roots, imbalance of intracellular osmotic regulation, and osmotic stress. In addition, salt stress affects biological functions in plant cells, for example, inhibiting its ability to uptake water and nutrients, accelerating chlorophyll degradation, inhibiting photosynthesis, and blocking various metabolic reactions (Karan and Subudhi, 2012; Deinlein et al., 2014; Al-Farsi et al., 2020). The most direct damage is high osmotic pressure resulting from low water use efficiency and ionic toxicity (mainly Na+) (Chinnusamy et al., 2006). Plants generally respond to salt stress through various mechanisms, including osmotic adjustment, ion transmembrane transport, ion compartmentalization, and active oxygen scavenging (ROS) (Liu et al., 2021).
In response to high salt stress, plants have evolved a number of signal transduction pathways involving regulatory genes and proteins for stress perception, stress signaling, and downstream metabolites (Bhattarai et al., 2020). Under salt stress, various stress signal receptors located on the cell membrane can quickly perceive changes in the external environment, and second messengers such as Ca2+, ROS, and phytohormone are rapidly generated in the cytoplasm to gradually decode and amplify salt stress signals (Fahad et al., 2014; Steinhorst and Kudla, 2019). These signals regulate downstream TFs through a cascade reaction to further alter transcript levels of many TF genes, such as AP2/ERF, bZIP, bHLH, WRKY, MYB, NAC, and DREB (Zhao et al., 2020). Thereafter, the expression of various osmotic stress-responsive genes (e.g., P5CSs and P5CRs) and ionic stress-responsive genes (e.g., HKTs, AKTs, and NSCCs) is eventually affected and ultimately contributes to plant salt tolerance (Zhu, 2016). Several signaling pathways involved in salt stress response have been revealed, including calcineurin B-like protein (CBL)–interacting protein kinase (CBL-CIPK), mitogen-activated protein kinase (MAPK), SOS, plant hormone (e.g., ABA, Eth, BR, and GAs), and calcium-dependent protein kinase (CDPK) pathways (Zhang et al., 2020; Zhao et al., 2020).
Smooth bromegrass is a perennial fine pasture of the grass family, which is widely used as hay, silage, and pasture for ruminants and dairy production (Ferdinandez and Coulman, 2001; Smart et al., 2006). It is mainly cultivated in northern China. However, due to climatic factors and poor agricultural management, soil salinization in these regions is becoming worse, which limits the production and utilization of smooth bromegrass (Wang and Yan, 2013 Wang et al., 2021). Therefore, it is important to explore molecular mechanisms of smooth bromegrass under salt stress and identify valuable resources for further breeding programs.
Currently, next-generation sequencing (NGS) provides an accurate and comprehensive analysis of differentially expressed genes and has made significant progress in understanding plant responses to drought and salt stress. NGS based on Illumina sequencing technology has been applied in Spartina pectinata (Robben and Gonzalez, 2016), Hordeum vulgare (Zhang et al., 2021), Spartina alterniflora (Hana, 2015), and Sorghum bicolor (Cui et al., 2018) to explore the molecular response mechanism under salt stress. Meanwhile, the single-molecule real-time (SMRT) sequencing provides third-generation sequencing (TGS) technology, which outperforms NGS technology in terms of read length (Rhoads and Au, 2015). In comparison to standard NGS technology, PacBio RSII TGS employs SMRT isoform sequencing (Iso-Seq), which provides longer read length, uniform coverage, and high accuracy (Dong et al., 2015). Therefore, SMRT sequencing technology has been used in Zea mays (Wang et al., 2016a), Sorghum bicolor (Abdel-Ghany et al., 2016), Medicago sativa (Luo et al., 2019), and Arabidopsis pumila (Yang et al., 2018) to explore how plants respond to stress.
In this study, we combined NGS and SMRT sequencing to generate a full-length transcriptome in smooth bromegrass. We analyzed the expression profile of differentially expressed genes in leaves and roots of smooth bromegrass under salt stress. In addition, we also identified a set of key genes that are likely involved in salinity adaptation in smooth bromegrass. These sequencing data will facilitate future research on gene function characterization and salinity adaptation mechanism in smooth bromegrass.
Materials and Methods
Plant Material and Growth
Bromus inermis cv. Wusu No. 1 is a variety that is widely cultivated at Qitai Grassland Station in Xinjiang. It possesses characteristics of high yield, superior drought tolerance, cold resistance, rapid regeneration, early greening, and disease resistance. The full-grained seeds of Wusu No. 1 smooth bromegrass were sterilized and placed into a 26°C/16°C (light/dark, 8h/16h) light incubator for germination. At the two-leaf stage, the seedlings were transplanted into soil containing vermiculite and perlite at a ratio of 3:1 and grew in a greenhouse. The plants were irrigated with a fresh 1/2 Hoagland solution every 3 days. When the seedlings grew to the 4-leaf stage, they were selected and carefully transplanted into a triangular flask containing a 1/2 Hoagland solution to recover for 3 days. These seedlings were then treated with a 1/2 Hoagland solution containing 300 mM NaCl, and leaves and roots from 10 seedlings (with a triplicate) were sampled at 0 h, 12 h, 24 h, 36 h, and 48 h, frozen in liquid nitrogen, and stored at −80 °C.
Analysis on Physiological Index
Chlorophyll was measured with a SPAD-502 chlorophyll meter, and the SPAD value was recorded with 10 repeats. A handy PEA portable plant fluorometer was used to determine the maximum photochemical efficiency. A DDS-11A conductivity meter was utilized to measure relative conductivity, and the relative water content method was used to measure fresh weight (FW) immediately after harvest. The plants were placed in deionized water at 4°C for 8 h. The weight (TW) saturated with water absorption was measured, dried at 105°C for 10 min, and moved to 80°C for 24 h. The dry weight (DW) of the plants was measured, and the RWC was calculated using the following formula: RWC (%) = (FW-DW)/(TW-DW) × 100. The parameters of MDA, proline, SOD, and GSH content were determined by using kits developed by Nanjing Jiancheng Bioengineering Research Institute Co., Ltd. Measurement methods were referred to the instructions of the plant malondialdehyde (MDA) assay kit (A003-3-1), proline assay kit (A107-1-1), superoxide dismutase (SOD) assay kit (A001-3-2), and reduced glutathione (GSH) assay kit (A006-2-1), respectively, as previously reported (Wang et al., 2016b, 2019; Meng et al., 2020).
Sequencing, Assembly, and Annotation of SMRT Library
RNAs were extracted by using an EasyPure Plant RNA Kit (No. ER301), and a Nanodrop 2000 spectrophotometer was used to detect the purity of RNA. A SMARTer PCR cDNA Synthesis Kit was employed for reverse transcription. The PacBio Sequel platform (Pacific Biosciences: http://www.pacb.com) was used for sequencing, and the PacBio official software package SMRTlink was used to process the original offline data. The subread sequence and the circular consensus sequence (CCS) were produced by correcting the subreads, and the sequences were then classified as full-length or non-full-length based on the presence of 5'-end primer, 3'-end primer, or poly-A tail. The full-length sequences were then clustered using ICE to get the cluster consensus sequence. In order to measure the correlation and similarity among samples, the PCA and Pearson correlation coefficient (R2) were analyzed by using the R package.
We used LoRDEC software (Leena and Eric, 2014) to analyze and correct the second-generation data and third-generation PacBio data. The corrected transcripts were sequenced and clustered by CD-HIT software (Fu et al., 2012), and redundant and similar sequences were removed. CDS predictive analysis was performed using ANGEL software (Shimizu et al., 2006), and the fault-tolerant mode was adopted by default. We compared all predicted protein coding sequences with protein and nucleotide databases using BLASTX with the following databases: NR, Nt, Pfam, KOG/COG, SWISS-PROT, KEGG, and GO. The cutoff E-value for NCBI non-redundant protein (NR) and NCBI non-redundant nucleotide (NT) was ≤ 10−5, and that for the protein family was E-value ≤ 10−5. The cutoff E-values were all ≤ 10−5 for KOG (http://www.ncbi.nlm.nih.gov/COG/), SWISS-PROT (Amos and Rolf, 2000), and KEGG analyses (http://www.genome.jp/kegg/) (Minoru et al., 2004). The setting parameters for Gene Ontology were referred to Ashburner et al. (http://www.geneontology.org/) (Ashburner et al., 2000), and GO enrichment analysis was performed by using agriGO v2.0 (p-value ≤ 0.05) (Tian et al., 2017). All differentially expressed genes were analyzed by using iTAK for online plant transcription factor prediction (http://itak.feilab.net/cgi-bin/itak/index.cgi) (Yi et al., 2016).
Illumina Sequencing and Sequence Assembly
Magnetic beads with oligo (dT) were used to enrich mRNA, and six-base random primers (random hexamers) were employed to synthesize cDNA. The double-stranded cDNAs were purified with AMPure XP beads, and PCR amplification was performed to construct the cDNA library. An Agilent 2,100 Bioanalyzer was used to detect the insert size of the library. After qualification, the Illumina second-generation high-throughput sequencing platform was used, and raw reads were filtered to obtain clean reads. Trinity (Grabherr et al., 2011) was used to assemble the transcripts of clean reads, and software RSEM (Dewey and Li, 2011) was used to calculate the gene expression using the FPKM method (Trapnell et al., 2010). Based on the average FPKM values under each treatment in roots and leaves, both the absolute values of the |log2 (fold change)| ≥1 and the adjusted P-value < 0.05 were used as thresholds to identify DEGs.
RT-qPCR Analysis
In order to verify the accuracy of the sequencing results, 10 differential genes were selected for RT-qPCR analysis. The TransScript One-Step gDNA Removal and cDNA Synthesis Super-Mix Reverse Transcription Kit were used for reverse transcription. Total RNAs extracted from leaves and roots from the control and salt treatment groups were used to synthesize cDNA. Primer Premier 5 software was used to design primers (Supplementary Table S1) (Regina et al., 2010; Bahrini et al., 2011). Real-time PCR was performed using the Applied Biosystems 7500/7500 fast real-time PCR system (ABI, Foster City, California, USA) and the SYBR Green PCR Master Mix system (Takara). We referred to the Trans-Start Top Green qPCR Super-Mix manual for RT-qPCR with triplicates. The 2−ΔΔCt method (Livak and Schmittgen, 2001) was used to quantitatively calculate the relative expression level of candidate genes.
Statistical Analysis
Statistical data analysis was performed using Excel 2016 software and the SPSS software package (ver. 22.0; SPSS Inc., Chicago, IL, USA). MEV 4.9 software (https://sourceforge.net/pro-jects/mev-tm4/files/mev-tm4/) was used for cluster analysis and expression pattern analysis.
Results
Effects of Salt Treatment on Physiological Changes of Smooth Bromegrass
Under the treatment of 300 mM NaCl, seedlings of smooth bromegrass gradually withered from 12 h to 48 h (Supplementary Figure S1). Compared with the control group (CK), the relative water content gradually decreased under NaCl treatment, while the relative conductivity gradually increased and reached the peak at 48 h (Figures 1A,B). No significant difference in the chlorophyll SPAD value was found for the treatment group when compared with the control group, but the photochemical efficiency was reduced by 10.98% after 48-h treatment compared with the CK (Figures 1C,D). Proline and Na+ contents increased during treatment time (Figures 1E,F), while the K+ content increased slightly at 12 h and decreased gradually from 36 h to 48 h (Figure 1G).
Figure 1. Analyses of dynamic physiological effects under 300 mM NaCl stress in the leaves. CK: untreated control; L1: salt stress treatment for 12h; L2: 24 h; L3: 36 h; L4: 48 h. (A) rRelative water content; (B) relative electrolyte leakage; (C) chlorophyll content (SPAD); (D) photochemical efficiency; (E) proline content; (F) Na+ content; (G) K+ content; (H) MDA content; (I) superoxide dismutase; (J) GSH content.
The malondialdehyde content can be used as an indicator of membrane lipid peroxidation. The malondialdehyde level in leaves of smooth bromegrass was considerably elevated by 37.19% at 48 h in the NaCl treatment group compared to the control (Figure 1H). In order to investigate the antioxidant defense system of smooth bromegrass cells under NaCl treatment, activities of key enzymes such as superoxide dismutase (SOD) and reduced glutathione (GSH) were tested in both control and treatment groups. It showed that SOD activity reached its peak after a 12-h treatment and started to decrease from 24 h to 48 h (Figure 1I). Unlike SOD, the GSH content increased considerably to a maximum at 24 h and then decreased from 36 h to 48 h (Figure 1J).
Statistical Analysis of Transcriptome Sequencing Data
The leaves and roots of smooth bromegrass treated with 300 mM NaCl at five time points were separately sampled and subjected to NGS sequencing. In total, about 226.11 G clean reads were generated by Illumina sequencing (Supplementary Table S2). The PacBio RS II platform was utilized to create the Iso-Seq library and single-molecule sequencing. In total, 19.67 G polymerase read bases and 516,677 polymerase reads were obtained. It showed that the final insert subread is 19.06 G, and the number of final insert subreads is 8,354,897 (Table 1).
Table 1. Full-length transcriptome sequencing results.
In total, 369,646 full-length (FL) reads, 41,968 non-full-length (NFL) reads, and 355,836 full-length non-chimeric (FLNC) reads (Supplementary Table S3) were obtained. To eliminate redundancy, we used the ICE algorithm to cluster FLNC sequences with the same transcript, and a total of 202,837 consensus numbers were obtained. The length of the polished consensus sequence for each sample and the length distribution map are shown in Supplementary Figures S2A,B. In order to test the repeatability between samples, both principal component analysis and Pearson correlation were performed for a total of 30 samples. The principal component analysis showed the samples for L0 and R0 from leaves and roots were grouped separately, which was far away from the samples treated with salt (Supplementary Figure S3). Correlation analysis showed the Pearson correlation coefficient between samples of the same treatment was greater than 0.8 (Supplementary Figure S4), indicating that sampling was relatively accurate with small error, and the data analyzed later are also relatively reliable. Furthermore, the consensus sequences obtained from the subreads yielded 413,135 circular consensus sequences (CCS), and the average number of CCS was 17. FLNC reads accounts for 86.13% of the CCS (Supplementary Table S5).
We used LoRDEC software (Leena and Eric, 2014) to correct and analyze the Illumina data and SMRT PacBio data, and CD-HIT software (Fu et al., 2012) to remove redundant and similar sequences, and finally, we obtained 202,837 transcripts (Table 2). Among all these transcripts, only 1.66% of them were >500 bp, 37.77% of them were 2-3kb, and 31.57% of them were longer than 3kb (Table 2). Our results confirmed that SMRT sequencing provides a large number of full-length and high-quality transcripts and that using Illumina data to correct low-quality SMRT reads improved the accuracy of long reads of PacBio (Figure 2A). The average length of genes detected by SMRT is more than that detected by Illumina. Finally, 116,578 unigenes were obtained after the removal of redundant sequences, and 1.11% of the retrieved unigenes were >500 bp, whereas 39.10% of them were longer than 3 kb (Table 2, Figure 2B).
Table 2. Statistical table of length frequency distribution before and after transcript de-redundancy.
Figure 2. Unigene length distribution map and gene function annotation Venn diagram. (A) Illumina and SMRT measured the gene count distribution map. The abscissa is the length (nt) of the predicted gene; the ordinate is the number of transcripts reads of the gene; (B) Illumina and SMRT measured gene length distribution box diagram; (C) seven database annotation results statistical graph; (D) gene function annotation Venn diagram.
Analysis of Gene Function Annotation
For more accurate gene analysis, we used ANGEL software (Shimizu et al., 2006) to predict CDS regions and obtained the distribution map of CDS (Supplementary Figure S3A). Unigenes were used for gene function annotation, with the following databases: Nr, Nt, Pfam, KOG/COG, SWISS-PROT, KEGG, and GO. In total, 66,847 (57.34%) unigenes had their annotation based on the GO database, and 104,784 (89.88%) and 107,022 (91.80%) unigenes based on KEGG and Nr databases, respectively (Figure 2C; Supplementary Table S4). Based on the plant species in the Nr database, a large number of unigenes were annotated to Hordeum vulgare, indicating that smooth bromegrass is closely related to Hordeum vulgare (Supplementary Figure S3B). Five major databases (viz., Nr, GO, KEGG, KOG, and Nt) jointly annotated 45,200 shared unigenes, and 523 were solely found in the NR database, 486 in the GO database, 32 in KEGG, 41 in KOG, and 3,679 in NT (Figure 2D).
Analysis of Differentially Expressed Genes
Compared with the control, 1,727 unigenes were specifically expressed in leaves under NaCl treatment at 12 h (L1), 7,523 at 24 h (L2), 1,458 at 36 h (L3), and 1,068 at 48 h (L4), and 4,753 of them were shared by these four comparison pairs (Figure 3A). Among differentially expressed genes, 7,954, 12,247, 6,283, and 5,552 unigenes were upregulated, and 4,845, 8,418, 4,913, and 4,224 unigenes were downregulated for the four comparison pairs, respectively (Figure 3B).
Figure 3. Distribution of DEGs in different periods. (A,B) Numbers of DEGs in leaves and roots of smooth bromegrass under different stress times compared with the control; (C,D) Summary of the numbers of up- and downregulated DEGs in leaves and roots of smooth bromegrass.
In roots, 1,162 unigenes were specifically expressed at the R1 stage under NaCl treatment in comparison with the control group, 1,018 at the R2 stage, 772 at the R3 stage, and 768 at the R4 stage, and a total of 630 unigenes were shared by the four comparison pairs (Figure 3C). Meanwhile, 3,455, 2,937, 2,764, and 2,399 unigenes were found to be differentially expressed in roots of smooth bromegrass at four different periods in comparison with the control group. Among them, 2,662, 2,051, 1,674, and 1,353 unigenes were upregulated, and 793, 886, 1,090, and 1,046 unigenes were downregulated at the four stages, respectively (Figure 3D).
Analysis of GO Annotation
The GO database was used for the annotation of DEGs between the control and treatment groups. The heat map for GO of DEGs was generated based on the calculated -log2(p-value), and a total of 17 biological processes were identified (Supplementary Table S5). Among them, the carbohydrate metabolism process, cellular amine metabolism process, and response to oxidative stress process showed higher expression levels in each stage under salt stress (Figure 4A). These data were also analyzed by using a direct acyclic graph (DAG) tree, resulting in the identification of response to oxidative stress, carbohydrate metabolic process, cellular amine metabolic process, and other important processes (Figure 4B).
Figure 4. GO enrichment results of DEGs after NaCl treatment. (A) GO annotation results of DEGs; (B–D) GO term for DEGs that respond to salt stress. Each box represents a GO term; the depth of the color represents the degree of enrichment, and the darker the color, the higher the degree of enrichment.
Furthermore, we found seven processes related to cell components (Supplementary Table S5), which were expressed differently in both leaves and roots. The DAG tree analysis also revealed important processes of the polyphosphate kinase complex and membrane (Figure 4C). By analyzing these data, we also found 18 processes related to molecular functions (Supplementary Table S5), and these processes varied in different stages under NaCl treatment. In particular, the catalytic activity process was significantly expressed in leaves of smooth bromegrass at all treatment time points (Figure 4A). The oxidoreductase activity process exhibited relatively high expression in various processing stages in both leaves and roots. In the process of peroxidase activity, catalase activity and antioxidant activity showed higher expression in roots. Moreover, several important molecular functional processes were also detected by using DAG analysis, such as catalytic activity, antioxidant activity, oxidoreductase activity, peroxidase activity, catalase activity, and iron ion binding (Figure 4D).
Analysis of KEGG Pathway
Analysis of the KEGG metabolic pathway was performed with DEGs between the control group and the four treatment groups in both leaves and roots, and 14 important pathways were significantly enriched (Supplementary Table S6), including ascorbate and aldarate metabolism, oxidative phosphorylation, glutathione metabolism, plant hormone signal transduction, phenylpropanoid biosynthesis, and starch and sucrose metabolism (Figure 5). These pathways were affected in different time periods under salt stress. The phenylpropanoid biosynthesis pathway, oxidative phosphorylation, phosphatidylinositol signaling system, and sugar metabolism were highly expressed, with more genes in leaves than in roots. Interestingly, the pathways of riboflavin metabolism only occurred in leaves but not in roots, while the carbon metabolism pathway, carbon fixation in photosynthetic organisms, and ascorbate and aldarate metabolism only occurred in roots but not in leaves, indicating they are tissue-specific pathways in smooth bromegrass (Figure 5).
Figure 5. Scatter plot of the enriched KEGG pathway of DEGs under NaCl stress.
Analyses of Gene Expression Levels Involved in Specific Pathways
In this study, we also analyzed ABA signaling pathway genes and identified 56 genes encoding signal receptor sensors related to NaCl stress through DEGs annotation analysis (Supplementary Figure S6A and Supplementary Table S7). In addition, 3 RLKs, 8 CDPKs, 20 CIPKs, 15 CMLs, and 10 MAPKs in the signal sensor pathway were also identified (Supplementary Figure S6B and Supplementary Table S8).
Meanwhile, we also found that most of the DEGs were associated with antioxidant enzymes (SODs, PODs, GPXs, CATs, and APXs), non-enzymatic antioxidants (GSTs, GSHs, and GSRs), and proline synthase (P5CSs, and P5CRs), and their corresponding coding genes were all regulated by NaCl stress at various levels (Figure 6A; Supplementary Table S9). We found that the expression levels of 9 AKTs, 11 AVPs, 7 CHXs, 6 CNGCs, 3 HKTs, 3 KATs, and 6 PMAs genes were also affected under NaCl treatment in both leaves and roots (Figure 6B; Supplementary Table S10). Furthermore, we identified 12 NHX-, 10 SOS1-, 5 SOS2-, and 4 SOS3-related genes (Figure 6C, Supplementary Table S11), and these genes were upregulated after salt stress, indicating that Na+ regulatory mechanism in smooth bromegrass was activated.
Figure 6. Heat map plot of the expression levels of the DEGs involved in antioxidative, SOS, and ion transporter pathway under NaCl stress. (A) Heat map plot of the antioxidative defense-related DEGs by NaCl stress; (B,C) heat map plot of the DEGs of SOS pathway and ion transporters by NaCl stress.
Identification of Salt-Responsive Transcription Factor Genes
In this study, a total of 1,223 transcription factor genes were identified and classified into 53 families. Among them, the bZIP family is the largest family with 127 genes, followed by 96 NAC-, 76 bHLH-, 73 AP2/ERF-, and 70 MYB-related genes (Figure 7; Supplementary Table S12). In order to investigate the expression of these TF genes, they were all analyzed at all the treatment time points. It was found that the expression levels of many transcription factor genes were highly upregulated at different time points after salt stress, including several main TF family genes AP2/ERF, bHLH, bZIP, GRAS, MYB, and WRKY (Supplementary Figure S7). In particular, many of them showed distinct expression profiles between leaves and roots (Supplementary Figure S7). The dynamic changes in the expression level of these TF genes implied their important roles in regulating plant salt tolerance.
Figure 7. Statistics of transcription factors under NaCl stress.
Validation of NaCl Stress-Related Genes by RT-qPCR
To validate the transcriptome data, 10 DEGs were randomly selected, and their expression levels were analyzed in the control and treatment groups by RT-qPCR (Supplementary Figure S8A), with actin gene as an internal reference. These genes included the putative peroxide gene, MYB39, and CYP73A2 (Supplementary Figure S8A). By comparing and analyzing these results with the RNA-Seq data (Supplementary Figure S8B), it was found that the expression of these 10 selected genes showed similar expression profiles in the leaves and roots at four time points (Supplementary Figure S8A). Further analyses showed that RT-qPCR and RNA-Seq data were correlated linearly either in leaves or in roots at 12 h, 24 h, 36 h, and 48 h (Supplementary Figure S8B–E), indicating RNA-Seq data were reliable for further gene screening and expression level comparison.
Discussion
At present, the third-generation sequencing technology (e.g., PacBio SMRT) has greatly facilitated the de novo assembly of transcriptomes of many higher organisms and has helped to overcome the problems of short splicing and incomplete information for the species without reference genomes (Cheng et al., 2017; Chao et al., 2019; Kan et al., 2020; Shen et al., 2021). Meanwhile, second-generation sequencing (e.g., Illumina platform) has been widely used to obtain more comprehensive annotation information of many plant species (Yang et al., 2018; Luo et al., 2019). In this study, for the first time, the Illumina platform was used to study the transcriptome changes of global genes in response to salt stress in the leaves and roots of smooth bromegrass, and the reference sequence was obtained from its full-length transcript database using PacBio SMRT sequencing technology. After screening, we acquired 19.67 G polymerase read bases (Table 1) and 116,578 distinct genes in 30 sample libraries (Table 2). Of these genes, 95.82% of DEGs were annotated in at least one database (Figure 3C), and five major databases (viz., Nr, GO, KEGG, KOG, and Nt) jointly annotated 45,200 unigenes (Figure 3D), which is higher than Elymus sibiricus (79.81%) (Zhou et al., 2016); however, it is lower than alfalfa (99%) (Luo et al., 2019). Compared with the control, 12,799, 20,665,11,196, and 9,776 DEGs were specifically expressed in leaves under NaCl treatment at 12 h, 24 h, 36 h, and 48 h, respectively (Figure 3B), and in roots, 3,455, 2,937, 2764, and 2,399 DEGs were identified in four different periods (Figure 3D). The candidate genes were analyzed in the control and treatment groups by RT-qPCR (Supplementary Figure S8A), which were consistent with the RNA-Seq data, implying our sequencing data were reliable (Supplementary Figure S8). Our study provided valuable and nearly complete sequence information on smooth bromegrass. Furthermore, we depicted a global molecular mechanism model of salt response in smooth bromegrass (Figure 8), with all these transcriptome information.
Figure 8. Model describing the signal transduction pathways involved in the acquisition of salt tolerance in smooth bromegrass. Red dashed lines indicate proposed regulation, solid arrowheads indicate activation, and vertical solid lines indicate suppression.
Under salt stress conditions, plants produced a large number of genes related to osmotic regulation, which regulate the accumulation of primary metabolites in plants such as proline, betaine, and sugar alcohol to maintain cell osmotic potential and improve their own salt tolerance (Xiong et al., 2017; Singh et al., 2018; Gao et al., 2019). Our studies showed that proline and soluble sugar contents increased in smooth bromegrass leaves under salt treatment (Figure 1), which was consistent with other studies in Arabidopsis thaliana (Yang et al., 2018), Medicago sativa (Luo et al., 2019), Sorghum bicolor (Cui et al., 2018), Arundo donax (Sicilia et al., 2020), and Zea mays (Kakumanu et al., 2012), indicating that plants share similar salt tolerance mechanisms. At the molecular level, the GO categories “response to oxidative stress,” “ion transport,” “catalytic activity,” “antioxidant activity,” and “hormone activity”, and the KEGG pathway categories “plant hormone signal transduction,” “phenylpropanoid biosynthesis,” and “flavonoid biosynthesis” were significantly enriched under salt stress (Figures 4, 5), which were consistent with rice (Oryza sativa) (Rabbani et al., 2003), maize (Zea mays) (Sun et al., 2015), and Sophora alopecuroides (Zhu et al., 2021). These results implied secondary metabolic processes and antioxidant systems are activated through secondary metabolites (phenylpropane and flavonoids), hormones, ion transport (Na+ and K+), antioxidant enzymes (SOD and GSH), and second messengers (such as ROS and Ca2+) in response to salt stress. Consistently, our physiological work also revealed that SOD and GSH contents increased at the physiological level by NaCl stress (Figures 1I,J), and increased levels of Na+ and K+ contents were consistent with the ion transport pathway (Figures 1F,G). The physiological regulation mechanism of salt stress usually includes osmotic regulation, ion transmembrane transport, ion compartmentalization, and active oxygen scavenging. Plant hormones and transcription factors also respond to salt stress (Zhu, 2016; Geng et al., 2020; Luo et al., 2020), which enabled smooth bromegrass to respond to salt stress and further adapt to changes in the external environment by ROS accumulation.
Under NaCl stress, the osmotic stress caused by the increase in the Na+ concentration in the rhizosphere of plants reduces the water potential. On the other hand, this ion stress induces a nutritional imbalance in plants (Huang et al., 2012; Al-Farsi et al., 2020). The active oxygen signal pathway is activated in response to osmotic stress, and excessive ROS causes serious damage to plant cell membranes, proteins, nucleic acids, etc. (Yun et al., 2011; Zhang et al., 2016). Plants produce redox-sensitive TFs and other molecular sensors to perceive stress signals, thus increasing the ROS concentration and activating ROS defense/metabolic pathways for ROS clearance (Chawla et al., 2013; Fidalgo et al., 2015; Mahmoud et al., 2020). Physiological studies on smooth bromegrass leaves revealed that NaCl produced oxidative damage, resulting in a substantial reduction in the chlorophyll content, an elevation in relative electrical conductivity, and a reduction in the relative water content (Figure 1). Smooth bromegrass activates the antioxidant defense system by adjusting the activity of antioxidants (POD, CAT, and SOD) and the content of osmotic regulators (PRO) (Supplementary Table S9). The same transcriptome data also revealed that NaCl stress activates the antioxidant defense system at the molecular level. When exposed to NaCl treatment GO-type oxidoreductase activity and KEGG pathways peroxisome and glutathione metabolism, as well as hormone signal transduction, were all enriched, as illustrated in Figures 4, 5. In addition, most of the DEGs associated with antioxidant enzymes (SODs, PODs, GPXs, CATs, and APXs), non-enzymatic antioxidants (GSTs, GSHs, and GSRs), and proline synthase (P5CSs and P5CRs) were all regulated by NaCl stress (Supplementary Table S9). These physiological indicators are consistent with the results of the transcriptome data, which strongly indicate the important role of these DEGs in responding to osmotic stress and therefore protecting smooth bromegrass from ROS damage.
Plants respond to external stress by activating signal pathways, including the ROS signal pathway, the MAPK cascade signal system, the ABA signal response pathway, and Ca2+-dependent proteins (Wang et al., 2016b; Danquah et al., 2014; Steinhorst and Kudla, 2019). Osmotic stress signals are sensed by protein receptors on the cell membrane. These stress signal receptors can quickly sense changes in the external environment. Subsequently, many second messengers such as Ca2+, inositol phosphate, ROS, and ABA signals will be rapidly produced in the cytoplasm to further decode and amplify the salt stress signal (Mehlmer et al., 2010; Zhu, 2016; Cao et al., 2020). The second messenger stimulates downstream signals through protein phosphorylation cascade changes, such as CDPKs, CBLs, CIPKs, CMLs, and MAPKs, linking the perception of external stimuli with cellular responses (Danquah et al., 2014; Bakshi et al., 2019). In this study, we identified a number of genes encoding signal receptor sensors related to NaCl stress through DEG annotation analysis (Supplementary Figure S4B), including CDPKs, CIPKs, CMLs, and MAPKs (Supplementary Table S8). Since roots are directly exposed to salt stress, the gene expression of smooth bromegrass exhibits distinct expression profiles in leaves and roots.
The characteristic of the CBL-CIPK regulatory network provides a great contribution to the efficient transduction of Ca2+ signals and also plays an important role in regulating the response of plants to salt stress (Thoday-Kennedy et al., 2015; Liu et al., 2021). The CBL-CIPK module can interact with many ion transporters and regulate their activity, especially Na+/K+ transporters, which are essential for ion homeostasis (Pandey et al., 2015; Fan et al., 2016). Studies have reported that cyclic nucleotide-gated channels (CNGCs) can mediate the flow of Ca2+ into the cytoplasm under stress and maintain ion balance in the cell (Kugler et al., 2009; Jin et al., 2015). In many cases, these ion transporters interact with the SOS signaling network in a synergistic or antagonistic manner to maintain ion “balance” in plants under salt stress conditions. In our study, the GO enrichment of DEGs revealed genes involved in ion homeostasis, iron ion transport, ion channel activity, ion transmembrane transporter activity, and iron ion binding, as well as the ion transporter (Figure 4). Meanwhile, genes involved in this signaling pathway were identified, including AKTs, AVPs, HKTs, and PMAs (Figure 6B; Supplementary Table S10), indicating that these DEGs may directly or indirectly promote the tolerance mechanism of smooth bromegrass to salt stress. The SOS signal pathway is an important signal transduction pathway for plant salt tolerance, and it plays an important role in regulating Na+ and K+ homeostasis and salt tolerance (Deinlein et al., 2014; Steinhorst and Kudla, 2019). In total, 31 genes related to NHXs and SOS pathways were identified (Supplementary Table S11), and they were largely upregulated by salt stress (Figure 6B), indicating that Na+ regulatory mechanism was activated in smooth bromegrass, which was consistent with an increase in the Na+ content (Figure 1F).
ABA is an important plant hormone that plays a key role in regulating plant growth and resisting adversity stress. The ABA signaling pathway is also activated in the early stage of osmotic stress (Raghavendra et al., 2010; Komatsu et al., 2020). Osmotic stress causes plants to accumulate a significant quantity of active oxygen, and NADPH oxidase in the ABA pathway has a certain mitigating effect on alleviating osmotic stress (Zhao et al., 2020). In this study, 36 ABA signaling pathway-related genes were analyzed, including 5 PYR/PYLs, 19 PP2Cs, and 12 SnRK2s (Supplementary Figure S6A; Supplementary Table S7), and the expression level of these genes increased significantly after salt stress treatment in smooth bromegrass.
TFs are key components in the osmotic stress signal pathway, and they are involved in the signal perception of salt stress and the expression of downstream key genes in response to salt stress. They participate in a variety of biological processes and play an important role in regulating plant growth and stress response (Hartmann et al., 2015; Yang et al., 2020). The transcriptome of smooth bromegrass revealed that a number of different transcription factor families were affected by salt stress, including bZIP, NAC, bHLH, WRKY, AP2/ERF, MYB-related, C3H, C2H2, GRAS, and MYB (Figure 7), which was similar to those major TF families as in wheat under salt stress (Malik et al., 2020). These transcription factor genes have also been found in response to salt stress of Arabidopsis thaliana (Yang et al., 2018), Caragana korshinskii (Li et al., 2016), and Cynodon dactylon (Hu et al., 2015). Transcription factors such as AP2/ERF, bHLH, bZIP, C2H2, C3H, GRAS, NAC, MYB-related, and WRKY showed different expression patterns in leaves and roots after NaCl stress (Supplementary Figure S7), indicating that NaCl had different effects on the complex transcriptional regulation of smooth bromegrass adaption to salt stress.
Data Availability Statement
The datasets presented in this study can be found in online repositories. All sequencing data of PacBio Iso-Seq and Illumina SGS RNA-Seq in this article can be obtained in NCBI SRA (http://www.ncbi.nlm.nih.gov/sra): SRR15634023 and SRR15633796-SRR15633815.
Author Contributions
JS and QL designed this experiment, performed the experiments, and drafted the manuscript. JS, YC, and LZ analyzed experimental data and visualized it. BZ, YP, and YC revised the manuscript and directed the study. All authors have read and agreed to the published version of the manuscript.
Funding
This research was funded by the China Agriculture Research System of MOF and MARA (CARS-34).
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.
Acknowledgments
We sincerely thank Drs. Rongchen Yang and Jing Yang at Northwest A&F University for their help and suggestions for this experiment and manuscript. We also thank Drs. Shuzhen Zhang, Yanhui Zhang, Yuxiang Wang, Lili Gu, and Chenjian Li for their help and advice during the experiment.
Supplementary Material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fpls.2022.917338/full#supplementary-material
References
Abdel-Ghany, S. E., Hamilton, M., Jacobi, J. L., Ngam, P., Devitt, N., Schilkey, F., et al. (2016). A survey of the sorghum transcriptome using single-molecule long reads. Nat. Commun. 7, 11706. doi: 10.1038/ncomms11706
Al-Farsi, S. M., Nawaz, A., Anees Ur, R., Nadaf, S. K., Al-Sadi, A. M., Siddique, K. H. M., et al. (2020). Effects, tolerance mechanisms and management of salt stress in lucerne (Medicago sativa). Crop. Pasture. Sci. 71, 411–428. doi: 10.1071/CP20033
Amos, B., and Rolf, A. (2000). The SWISS-PROT protein sequence database and its supplement TrEMBL in 2000. Nucl. Acids Res. 28, 45–48. doi: 10.1093/nar/28.1.45
Ashburner, M., Ball, C. A., Blake, J. A., Botstein, D., Butler, H., Cherry, J. M., et al. (2000). Gene Ontology: tool for the unification of biology. Nat. Genet. 25, 25–29. doi: 10.1038/75556
Bahrini, I., Ogawa, T., Kobayashi, F., Kawahigashi, H., and Handa, H. (2011). Overexpression of the pathogen-inducible wheat TaWRKY45 gene confers disease resistance to multiple fungi in transgenic wheat plants. Breed 61, 319–326. doi: 10.1270/jsbbs.61.319
Bakshi, P., Handa, N., Gautam, V., Kaur, P., Sareen, S., Mir, B., et al. (2019). Role and Regulation of Plant Hormones as a Signal Molecule in Response to Abiotic Stresses. Cambridge: Woodhead Publ Ltd. p. 303–317. doi: 10.1016/B978-0-12-816451-8.00018-6
Bhattarai, S., Biswas, D., Fu, Y. B., and Biligetu, B. (2020). Morphological, physiological, and genetic responses to salt stress in alfalfa: a review. Agronomy-Basel. 10, 577. doi: 10.3390/agronomy10040577
Cao, Y., Zhang, M., Liang, X., Li, F., and Jiang, C. (2020). Natural variation of an EF-hand Ca2+-binding-protein coding gene confers saline-alkaline tolerance in maize. Nat. Commun. 11, 186. doi: 10.1038/s41467-019-14027-y
Chao, Y. H., Yuan, J. B., Guo, T., Xu, L. X., Mu, Z. Y., and Han, L. B. (2019). Analysis of transcripts and splice isoforms in Medicago sativa L. by single-molecule long-read sequencing. Plant Mol. Biol. 99, 219–235. doi: 10.1007/s11103-018-0813-y
Chawla, S., Jain, S., and Jain, V. (2013). Salinity induced oxidative stress and antioxidant system in salt-tolerant and salt-sensitive cultivars of rice (Oryza Sativa L.). J. Plant Biochem. Biot. 22, 27–34. doi: 10.1007/s13562-012-0107-4
Cheng, B., Furtado, A., and Henry, R. J. (2017). Long-read sequencing of the coffee bean transcriptome reveals the diversity of full-length transcripts. Gigascience. 6, 13. doi: 10.1093/gigascience/gix086
Chinnusamy, V., Zhu, J., and Zhu, J. K. (2006). Salt stress signaling and mechanisms of plant salt tolerance. Genet. Eng. 27, 141–177. doi: 10.1007/0-387-25856-6_9
Cui, J., Ren, G., Qiao, H., Xiang, X., Huang, L., and Chang, J. (2018). Comparative transcriptome analysis of seedling stage of two Sorghum cultivars under salt stress. J. Plant Growth Regul. 37, 986–998. doi: 10.1007/s00344-018-9796-9
Danquah, A., Zelicourt, A. D., Colcombet, J., and Hirt, H. (2014). The role of ABA and MAPK signaling pathways in plant abiotic stress responses. Biotechnol. Adv. 32, 40–52. doi: 10.1016/j.biotechadv.2013.09.006
Deinlein, U., Stephan, A. B., Horie, T., Luo, W., and Schroeder, J. I. (2014). Plant salt-tolerance mechanisms. Trends Plant Sci. 19, 371–379. doi: 10.1016/j.tplants.2014.02.001
Dewey, C. N., and Li, B. (2011). RSEM: accurate transcript quantification from RNA-Seq data with or without a reference genome. BMC Bioinform. 12, 323–323. doi: 10.1186/1471-2105-12-323
Dong, L., Liu, H., Zhang, J., Yang, S., Kong, G., Chu, J. S. C., et al. (2015). Single-molecule real-time transcript sequencing facilitates common wheat genome annotation and grain transcriptome research. BMC Genom. 16, 1039. doi: 10.1186/s12864-015-2257-y
Fahad, S., Hussain, S., Matloob, A., Khan, F. A., and Huang, J. (2014). Phytohormones and plant responses to salinity stress: a review. Plant Growth Regul. 75, 391–404. doi: 10.1007/s10725-014-0013-y
Fan, Z., Li, L., Jiao, Z., Chen, Y., and Zheng, J. (2016). Characterization of the Calcineurin B-Like (CBL) gene family in maize and functional analysis of ZmCBL9 under abscisic acid and abiotic stress treatments. Plant Sci. 253, 118–129. doi: 10.1016/j.plantsci.2016.09.011
Ferdinandez, Y., and Coulman, B. E. (2001). Nutritive values of Smooth Bromegrass, Meadow Bromegrass, and Meadow × Smooth Bromegrass hybrids for different plant parts and growth stages. Crop Sci. 41, 473–478. doi: 10.2135/cropsci2001.412473x
Fidalgo, F., Santos, A., Santos, I., and Salema, R. (2015). Effects of long-term salt stress on antioxidant defence systems, leaf water relations and chloroplast ultrastructure of potato plants. Ann. Appl. Biol. 145, 185–192. doi: 10.1111/j.1744-7348.2004.tb00374.x
Fu, L., Niu, B., Zhu, Z., Wu, S., and Li, W. (2012). CD-HIT: accelerated for clustering the next-generation sequencing data. Bioinformatics. 28, 3150–3152. doi: 10.1093/bioinformatics/bts565
Gao, Y. L., Cui, Y. J., Long, R. C., Sun, Y., Zhang, T. J., Yang, Q. C., et al. (2019). Salt-stress induced proteomic changes of two contrasting alfalfa cultivars during germination stage. J. Sci. Food Agric. 99, 1384–1396. doi: 10.1002/jsfa.9331
Geng, W., Li, Z., Hassan, M. J., and Peng, Y. (2020). Chitosan regulates metabolic balance, polyamine accumulation, and Na+ transport contributing to salt tolerance in creeping bentgrass. BMC Plant Biol. 20, 506. doi: 10.1186/s12870-020-02720-w
Grabherr, M. G., Haas, B. J., Yassour, M., Levin, J. Z., Thompson, D. A., Amit, I., et al. (2011). Full-length transcriptome assembly from RNA-Seq data without a reference genome. Nat. Biotech. 29, 644–652. doi: 10.1038/nbt.1883
Hana, Z., karimi Renesh, B.edre, Julio, S.olis, Venkata, M., angu, Niranjan, and Baisakh (2015). Sequencing and expression analysis of salt-responsive miRNAs and target genes in the halophyte smooth cordgrass (Spartina alternifolia L.). Mol. Biol. Rep. 42, 1341–1350. doi: 10.1007/s11033-015-3880-z
Hartmann, L., Pedrotti, L., Weiste, C., Fekete, A., Schierstaedt, J., GöTtler, J., et al. (2015). Crosstalk between two bZIP signaling pathways orchestrates salt-induced metabolic reprogramming in Arabidopsis roots. Plant Cell. 27, 2244–2260. doi: 10.1105/tpc.15.00163
Hu, L. X., Li, H. Y., Chen, L., Lou, Y. H., and Amombo, E. (2015). RNA-seq for gene identification and transcript profiling in relation to root growth of bermudagrass (Cynodon dactylon) under salinity stress. BMC Genomics. 16, 1–12. doi: 10.1186/s12864-015-1799-3
Huang, G. T., Ma, S. L., Bai, L. P., Li, Z., Hui, M., Jia, P., et al. (2012). Signal transduction during cold, salt, and drought stresses in plants. Mol. Biol. Rep. 39, 969–987. doi: 10.1007/s11033-011-0823-1
Ivushkin, K., Bartholomeus, H., Bregt, A. K., Pulatov, A., and Sousa, L. D. (2019). Global mapping of soil salinity change. Remote Sens. Environ. 231, 111260. doi: 10.1016/j.rse.2019.111260
Jin, Y., Jing, W., Zhang, Q., and Zhang, W. (2015). Cyclic nucleotide gated channel 10 negatively regulates salt tolerance by mediating Na+ transport in Arabidopsis. J. Plant Res. 128, 211–220. doi: 10.1007/s10265-014-0679-2
Kakumanu, A., Ambavaram, M. M. R., Klumas, C., Krishnan, A., Batlang, U., Myers, E., et al. (2012). Effects of drought on gene expression in maize reproductive and leaf meristem tissue revealed by RNA-Seq. Plant Physiol. 160, 846–867. doi: 10.1104/pp.112.200444
Kan, L., Liao, Q. C., Su, Z. Y., Tan, Y. S., Wang, S. Y., and Zhang, L. (2020). Single-Molecule Real-Time sequencing of the Madhuca pasquieri (Dubard) Lam. transcriptome reveals the diversity of full-length transcripts. Forests 11, 866. doi: 10.3390/f11080866
Karan, R., and Subudhi, P. K. (2012). Approaches to Increasing Salt Tolerance in Crop Plants. New York, NY: Springer. p. 63–88. doi: 10.1007/978-1-4614-0634-1_4
Kearl, J., Mcnary, C., Lowman, J. S., Mei, C., and Nielsen, B. L. (2019). Salt-tolerant halophyte rhizosphere bacteria stimulate growth of alfalfa in salty soil. Front. Microbiol. 10, 1849. doi: 10.3389/fmicb.2019.01849
Komatsu, K., Takezawa, D., and Sakata, Y. (2020). Decoding ABA and osmostress signalling in plants from an evolutionary point of view. Plant. Cell Environ. 43, 2894–2911. doi: 10.1111/pce.13869
Kugler, A., K?Hler, B., Palme, K., Wolff, P., and Dietrich, P. (2009). Salt-dependent regulation of a CNG channel subfamily in Arabidopsis. BMC Plant Biol. 9, 140. doi: 10.1186/1471-2229-9-140
Leena, S., and Eric, R. (2014). LoRDEC: accurate and efficient long read error correction. Bioinformatics 30, 3506–3514. doi: 10.1093/bioinformatics/btu538
Li, S., Fan, C., Li, Y., Zhang, J., Sun, J., Chen, Y., et al. (2016). Effects of drought and salt-stresses on gene expression in Caragana korshinskii seedlings revealed by RNA-Seq. BMC Genom. 17, 1–19. doi: 10.1186/s12864-016-2562-0
Liu, C., Mao, B., Yuan, D., Chu, C., and Duan, M. (2021). Salt tolerance in rice: physiological responses and molecular mechanisms. The Crop J. 10, 13–25 doi: 10.1016/j.cj.2021.02.010
Livak, K. J., and Schmittgen, T. D. (2001). Analysis of relative gene expression data using real-time quantitative PCR and the 2−ΔΔCT Method. Methods. 25, 402–408. doi: 10.1006/meth.2001.1262
Luo, D., Zhou, Q., Wu, Y., Chai, X., Liu, W., Wang, Y., et al. (2019). Full-length transcript sequencing and comparative transcriptomic analysis to evaluate the contribution of osmotic and ionic stress components towards salinity tolerance in the roots of cultivated alfalfa (Medicago sativa L.). BMC Plant Biol. 19, 32. doi: 10.1186/s12870-019-1630-4
Luo, H., Zhou, Z., Song, G., Yao, H., and Han, L. (2020). Antioxidant enzyme activity and microRNA are associated with growth of Poa pratensis callus under salt stress. Plant Biotechnol. Rep. 14, 429–438. doi: 10.1007/s11816-020-00620-x
Mahmoud, Z. S., El-Sayed, M. D., Ali, O., and Mostafa, M. R. (2020). Pumpkin seed protein hydrolysate treatment alleviates salt stress effects on Phaseolus vulgaris by elevating antioxidant capacity and recovering ion homeostasis. Sci. Hortic. 271, 109495, doi: 10.1016/j.scienta.2020.109495
Malik, A., Gul, A., Hanif, U., Kubra, G., Bibi, S., Ali, M., et al. (2020). “Salt responsive transcription factors in wheat”, in Climate Change and Food Security with Emphasis on Wheat. Elsevier, Pakistan. 107–127. doi: 10.1016/B978-0-12-819527-7.00007-8
Mehlmer, N., Wurzinger, B., Stael, S., Hofmann-Rodrigues, D., Csaszar, E., Pfister, B., et al. (2010). The Ca2+-dependent protein kinase CPK3 is required for MAPK-independent salt-stress acclimation in Arabidopsis. Plant J. 63, 484–498. doi: 10.1111/j.1365-313X.2010.04257.x
Meng, X., Liu, S., Dong, T., Xu, T., and Zhu, M. (2020). Comparative transcriptome and proteome analysis of salt-tolerant and salt-sensitive sweet potato and overexpression of IbNAC7 confers salt tolerance in Arabidopsis. Front. Plant Sci. 11, 572540. doi: 10.3389/fpls.2020.572540
Minoru, K., Susumu, G., Shuichi, K., Yasushi, O., and Masahiro, H. (2004). The KEGG resource for deciphering the genome. Nucl. Acids Res. 32, 277–280. doi: 10.1093/nar/gkh063
Pandey, G. K., Kanwar, P., Singh, A., Steinhorst, L., Pandey, A., Yadav, A. K., et al. (2015). Calcineurin B-like protein-interacting protein kinase CIPK21 regulates osmotic and salt stress responses in Arabidopsis. Plant Physiol. 169, 780–792. doi: 10.1104/pp.15.00623
Rabbani, M. A., Maruyama, K., Abe, H., Khan, M. A., Katsura, K., Ito, Y., et al. (2003). Monitoring expression profiles of rice genes under cold, drought, and high salinity stresses and abscisic acid application using cDNA microarray and RNA gel blot analyses. Plant Physiol. 133, 1755–67. doi: 10.1104/pp.103.025742
Raghavendra, A. S., Gonugunta, V. K., Christmann, A., and Grill, E. (2010). ABA perception and signalling. Trends Plant Sci. 15, 395–401. doi: 10.1016/j.tplants.2010.04.006
Regina, A., Kosar-Hashemi, B., Ling, S., Li, Z., Rahman, S., and Morell, M. (2010). Control of starch branching in barley defined through differential RNAi suppression of starch branching enzyme IIa and IIb. J. Exper. Bota. 61, 1469–1482. doi: 10.1093/jxb/erq011
Rengasamy, P. (2010). Soil processes affecting crop production in salt-affected soils. Plant Biol. 37, 613–620. doi: 10.1071/FP09249
Rhoads, A., and Au, K. F. (2015). Sequencing and its applications. Genomics, Proteom. Bioinform. 13, 278–289. doi: 10.1016/j.gpb.2015.08.002
Robben, M., and Gonzalez, J. (2016). RNA-Seq analysis of the salt tolerance response in prairie cordgrass. The FASEB J. 30, lb165–lb165. doi: 10.1096/fasebj.30.1_supplement.lb165
Shen, G. M., Ou, S. Y., He, C., Liu, J., and He, L. (2021). Full length sequencing reveals novel transcripts of detoxification genes along with related alternative splicing events and lncRNAs in Phyllotreta striolata. PLoS ONE 16, e0248749. doi: 10.1371/journal.pone.0248749
Shimizu, K., Adachi, J., and Muraoka, Y. (2006). ANGLE: a sequencing errors resistant program for predicting protein coding regions in unfinished cDNA. J. Bioinform. Comput. Biol. 4, 649–664. doi: 10.1142/S0219720006002260
Sicilia, A., Santoro, D. F., Testa, G., Cosentino, S. L., and Piero, A. (2020). Transcriptional response of giant reed (Arundo donax L.) low ecotype to long-term salt stress by unigene-based RNA-seq. Phytochemist. 177, 112436. doi: 10.1016/j.phytochem.2020.112436
Singh, V., Singh, A. P., Bhadoria, J., Giri, J., Singh, J.V., et al. (2018). Differential expression of salt-responsive genes to salinity stress in salt-tolerant and salt-sensitive rice (Oryza sativa L.) at seedling stage. Protoplasma. 255, 1667–1681. doi: 10.1007/s00709-018-1257-6
Smart, A. J., Schacht, W. H., Volesky, J. D., and Moser, L. E. (2006). Seasonal changes in dry matter partitioning, yield, and crude protein of intermediate wheatgrass and Smooth Bromegrass. Agron. J. 98, 986–991. doi: 10.2134/agronj2005.0233
Steinhorst, L., and Kudla, J. (2019). How plants perceive salt. Nature 572, 318–320. doi: 10.1038/d41586-019-02289-x
Sun, C., Gao, X., Fu, J., Zhou, J., and Wu, X. (2015). Metabolic response of maize (Zea mays L.) plants to combined drought and salt stress. Plant Soil. 388:99–117. doi: 10.1007/s11104-014-2309-0
Thoday-Kennedy, E. L., Jacobs, A. K., and Roy, S. J. (2015). The role of the CBL-CIPK calcium signalling network in regulating ion transport in response to abiotic stress. Plant Growth Regul. 76, 3–12. doi: 10.1007/s10725-015-0034-1
Tian, T., Liu, Y., Yan, H., You, Q., Yi, X., Du, Z., et al. (2017). agriGO v2.0: a GO analysis toolkit for the agricultural community, 2017 update. Nucl. Acids Res. 45, W122–W129. doi: 10.1093/nar/gkx382
Trapnell, C., Williams, B. A., Pertea, G., Mortazavi, A., Kwan, G., Baren, M., et al. (2010). Transcript assembly and quantification by RNA-Seq reveals unannotated transcripts and isoform switching during cell differentiation. Nat. Biotechnol. 28, 511–515. doi: 10.1038/nbt.1621
Wang, B., Tseng, E., Regulski, M., Clark, T. A., Hon, T., Jiao, Y., et al. (2016a). Unveiling the complexity of the maize transcriptome by single-molecule long-read sequencing. Nat. Commun. 7, 11708. doi: 10.1038/ncomms11708
Wang, X., Gao, F., Bing, J., Sun, W., Feng, X., Ma, X. F., et al. (2019). Overexpression of the Jojoba Aquaporin gene, ScPIP1, enhances drought and salt tolerance in transgenic Arabidopsis. Int. J. Mol. Sci. 20, 153. doi: 10.3390/ijms20010153
Wang, X., Li, X., Cai, D., Lou, J., and Liu, F. (2021). Salinification and salt transports under aeolian processes in potential desertification regions of China. Sci. Total Environ. 782, 146832. doi: 10.1016/j.scitotenv.2021.146832
Wang, X. L., Li, X., Zhang, S., Korpelainen, H., and Li, C. Y. (2016b). Physiological and transcriptional responses of two contrasting Populus clones to nitrogen stress. Tree Physiol. 36, 628–644. doi: 10.1093/treephys/tpw019
Wang, Y., and Yan, L. (2013). Land exploitation resulting in soil salinization in a desertoasis ecotone. Catena 100, 50–56. doi: 10.1016/j.catena.2012.08.005
Xiong, J. B., Sun, Y., Yang, Q. C., Tian, H., Zhang, H. S., Liu, Y., et al. (2017). Proteomic analysis of early salt stress responsive proteins in alfalfa roots and shoots. Proteome Sci.15, 1–19. doi: 10.1186/s12953-017-0127-z
Yang, L., Jin, Y., Huang, W., Sun, Q., Liu, F., and Huang, X. (2018). Full-length transcriptome sequences of ephemeral plant Arabidopsis pumila provides insight into gene expression dynamics during continuous salt stress. BMC Genomics 19, 717–731. doi: 10.1186/s12864-018-5106-y
Yang, Y., Yu, T. F., Ma, J., Chen, J., and Xu, Z. S. (2020). The Soybean bZIP transcription factor gene GmbZIP2 confers drought and salt resistances in transgenic plants. Int. J. Mol. Sci. 21, 670. doi: 10.3390/ijms21020670
Yi, Z., Chen, J., Sun, H., Rosli, H. G., Pombo, M. A., Zhang, P., et al. (2016). iTAK: a program for genome-wide prediction and classification of plant transcription factors, transcriptional regulators, and protein kinases. Mol. Plant 9, 1667–1670. doi: 10.1016/j.molp.2016.09.014
Yun, S. B., Oh, H., Rhee, S. G., and Yoo, Y. D. (2011). Regulation of reactive oxygen species generation in cell signaling. Mol. Cells 32, 491–509. doi: 10.1007/s10059-011-0276-3
Zhang, H., Zhao, Y., and Zhu, J. K. (2020). Thriving under stress: how plants balance growth and the stress response. Dev. Cell 55, 529–543. doi: 10.1016/j.devcel.2020.10.012
Zhang, M., Smith, J. A. C., Harberd, N. P., and Jiang, C. F. (2016). The regulatory roles of ethylene and reactive oxygen species (ROS) in plant salt stress responses. Plant Mol. Biol. 91, 651–659. doi: 10.1007/s11103-016-0488-1
Zhang, X., Zhang, L., Chen, Y., Wang, S., Fang, Y., Zhang, X., et al. (2021). Genome-wide identification of the SOD gene family and expression analysis under drought and salt stress in barley. Plant Growth Regul. 94, 49–60. doi: 10.1007/s10725-021-00695-8
Zhao, C., Zhang, H., Song, C., Zhu, J. K., and Shabala, S. (2020). Mechanisms of plant responses and adaptation to soil Salinity. Innovat. 1, 100017. doi: 10.1016/j.xinn.2020.100017
Zhou, Q., Luo, D., Ma, L., Xie, W., Wang, Y., Wang, Y., et al. (2016). Development and cross-species transferability of EST-SSR markers in Siberian wildrye (Elymus sibiricus L.) using Illumina sequencing. Sci. Rep. 6, 20549. doi: 10.1038/srep20549
Zhu, J. K. (2016). Abiotic stress signaling and responses in plants. CELL 167, 313–324. doi: 10.1016/j.cell.2016.08.029
Keywords: smooth bromegrass, full-length transcriptomics, salt tolerance, ion signal transduction, SOS regulatory pathways
Citation: Li Q, Song J, Zhou Y, Chen Y, Zhang L, Pang Y and Zhang B (2022) Full-Length Transcriptomics Reveals Complex Molecular Mechanism of Salt Tolerance in Bromus inermis L. Front. Plant Sci. 13:917338. doi: 10.3389/fpls.2022.917338
Received: 11 April 2022; Accepted: 09 May 2022;
Published: 09 June 2022.
Edited by:
Jiyu Zhang, Lanzhou University, ChinaReviewed by:
Jianxiu Liu, Institute of Botany, Jiangsu Province (CAS), ChinaMarcelo Nogueira Do Amaral, Federal University of Pelotas, Brazil
Copyright © 2022 Li, Song, Zhou, Chen, Zhang, Pang and Zhang. 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.
*Correspondence: Bo Zhang, xjauzb@126.com; Yongzhen Pang, pangyongzhen@caas.cn
†These authors have contributed equally to this work