Next Article in Journal
Separation and Paired Proteome Profiling of Plant Chloroplast and Cytoplasmic Ribosomes
Previous Article in Journal
Identification and Expression Analysis of Cold Shock Protein 3 (BcCSP3) in Non-Heading Chinese Cabbage (Brassica rapa ssp. chinensis)
Previous Article in Special Issue
Spinach Plants Favor the Absorption of K+ over Na+ Regardless of Salinity, and May Benefit from Na+ When K+ is Deficient in the Soil
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

RNAseq Analysis Reveals Altered Expression of Key Ion Transporters Causing Differential Uptake of Selective Ions in Canola (Brassica napus L.) Grown under NaCl Stress

1
Department of Botany, Government College University, Lahore 54000, Pakistan
2
Department of Botany, Lahore College for Women University, Lahore 54000, Pakistan
3
Institute of Pure and Applied Biology, Bhauddin Zakria University, Multan 66000, Pakistan
4
Department of Plant Physiology, Institute of Biology, Warsaw University of Life Sciences SGGW, Nowoursynowska 159, 02-776 Warsaw, Poland
*
Authors to whom correspondence should be addressed.
Submission received: 15 May 2020 / Revised: 5 July 2020 / Accepted: 9 July 2020 / Published: 14 July 2020
(This article belongs to the Special Issue Mechanism of Salinity Tolerance in Plants)

Abstract

:
Salinity is one of the major abiotic stresses prevailing throughout the world that severely limits crop establishment and production. Every crop has an intra-specific genetic variation that enables it to cope with variable environmental conditions. Hence, this genetic variability is a good tool to exploit germplasms in salt-affected areas. Further, the selected cultivars can be effectively used by plant breeders and molecular biologists for the improvement of salinity tolerance. In the present study, it was planned to identify differential expression of genes associated with selective uptake of different ions under salt stress in selected salt-tolerant canola (Brassica napus L.) cultivar. For the purpose, an experiment was carried out to evaluate the growth response of different salt-sensitive and salt-tolerant canola cultivars. Plants were subjected to 200 mM NaCl stress. Canola cultivars—Faisal Canola, DGL, Dunkled, and CON-II—had higher growth than in cvs Cyclone, Ac-EXcel, Legend, and Oscar. Salt-tolerant cultivars were better able to maintain plant water status probably through osmotic adjustment as compared to salt-sensitive cultivars. Although salt stress increased shoot Na+ and shoot Cl contents in all canola cultivars, salt-tolerant cultivars had a lower accumulation of these toxic nutrients. Similarly, salt stress reduced shoot K+ and Ca2+ contents in all canola cultivars, while salt-tolerant cultivars had a higher accumulation of K+ and Ca2+ in leaves, thereby having greater shoot K+/Na+ and Ca2+/Na+ ratios. Nutrient utilization efficiency decreased significantly in all canola cultivars due to the imposition of salt stress; however, it was greater in salt-tolerant cultivars—Faisal Canola, DGL, and Dunkled. Among four salt-tolerant canola cultivars, cv Dunkled was maximal in physiological attributes, and thus differentially expressed genes (DEGs) were assessed in it by RNA-seq analysis using next-generation sequencing (NGS) techniques. The differentially expressed genes (DEG) in cv Dunkled under salt stress were found to be involved in the regulation of ionic concentration, photosynthesis, antioxidants, and hormonal metabolism. However, the most prominent upregulated DEGs included Na/K transporter, HKT1, potassium transporter, potassium channel, chloride channel, cation exchanger, Ca channel. The RNA-seq data were validated through qRT-PCR. It was thus concluded that genes related to the regulation of ionic concentrate are significantly upregulated and expressed under salt stress, in the cultivar Dunkled.

Graphical Abstract

1. Introduction

Salinity stress negatively affects the plant growth and yield of all glycophytic plants and cause multi-million dollar losses worldwide [1]. Growth reduction in plants occurs in salt-affected soil due to osmotic and ion toxic effects and limited uptake of nutrients [2]. However, due to the unique genetic makeup of every plant species, the extent of inhibitory effects varies greatly in all plant species, and some of them can cope with this problem better as compared to others. Canola (Brassica napus L.) has been reported to have the potential to cope with the toxicity of excess salt in the soil, and hence it can be successfully grown on salt-affected areas [3]. In an earlier study, screening was carried out for salt stress tolerance of 34 accessions of B. napus, and the selection of physiological parameters was done. Cultivar Dunkled was ranked salt stress-tolerant [4]. Brassica napus L. (Canola) ranks second after soybean, among oilseed crops, and provides 13% of the total oil supply in the world. However, further crop improvement up to 60% using current approaches by 2050 in this potential oilseed crop seems to be difficult. Thus, some novel sources of salt tolerance in this potential oilseed crop need to be explored. However, the problem is a lack of understanding of the detailed mechanism of salt tolerance or exact cellular and metabolic sites of salt-induced damages [5,6,7]. Several scientists have suggested that understanding how plants sense salt stress and respond to it through the various physiological processes will help us in devising a strategy to develop salt-tolerant plants [8]. For example, when plants are exposed to high salinity stress, membrane proteins at root hairs or root epidermal cells sense the extent of salt stress and communicate to whole plant body and alter different biochemical pathways, resulting in adjustments in ion homeostasis, detoxification of reactive oxygen species, changes in CO2 assimilation hormones, and growth regulation [9,10]. All these changes in the physiological process in response to salt stress are governed by the expression of salt stress-responsive genes or trigger the expression of some other related genes [11]. Salt overly sensitive (SOS) pathway is known to control ionic homeostasis in association with calcium-binding protein sensors that are important in signal transduction under salt stress [12,13,14]. The K+/Na+ balance under salt stress is also decisive in studying growth responses of plants under stress tolerance, and several transporters and their genes have been identified [14]. Serine/threonine type protein kinase CIPK24 is observed to be interacting with the Ca2+ sensor CBL4 and regulates the Na+/H+ exchangers and deliberates salinity tolerance [12]. These reports have suggested that genes can be grouped in sensing and signaling, ion transporters, and salt stress-related gene and transcriptional regulators. Identification of these genes playing their role in salt tolerance is necessary to develop salt-tolerant crops either through breeding or genetic engineering techniques.
In recent years, physiological processes complemented with proteome profiling and/or transcriptome mapping in plants have gained popularity to explore mechanisms of salt tolerance [15,16,17,18]. Next-generation sequencing (NGS) has become a potential technique to evaluate molecular profiles for crop plants [19]. It is being applied successfully throughout the world to identify transcriptome variation in plants [17,20]. Few studies on B. napus are available for the identification of salt-tolerant genes using the RNA-seq approach [16]. Keeping in view the above-mentioned facts, the first part of the present study was planned to evaluate salt tolerance in canola (Brassica napus L.) cultivars differing in salinity tolerance using the physiological and biochemical approach and to reconfirm behavior of cultivar Dunkled. In the second part, the transcriptome analysis (RNA-seq approach) was carried out in selected salt-tolerant cultivar Dunkled to identify key pathways and genes responsible for salt tolerance. Transcriptome analysis was further validated through gene expression of selected genes/salinity responsive transcripts. Our findings would help plant breeders and molecular biologists to improve cultivated canola varieties for better use of salt-affected areas all over the world.

2. Results

2.1. Evaluation of Different Cultivars of Canola (B. napus L.) for Salt Stress Tolerance

2.1.1. Growth Attributes

Analysis of variance of the data for shoot fresh and dry weights of eight canola cultivars differing in salinity tolerance revealed that the imposition of salt stress caused a significant reduction in shoot fresh and dry weight (Table 1). Canola cultivars were significantly different in their fresh and dry weights under normal or saline (200 mM NaCl) conditions (Figure 1). As expected, salt-tolerant canola cultivars (DGL, Dunkled, Faisal Canola, and CON-II) had greater fresh and dry weights than those of salt-sensitive canola cultivars (Legend, Oscar, Cyclone, and Ac-Excel). However, Dunkled had maximum growth potential under saline conditions. However, cultivars—Legend and Oscar—were highly salt-sensitive being lowest in their fresh and dry weights of shoots.

2.1.2. Relative Water Contents (RWC %)

Data for shoot dry weights of all canola cultivars under the study showed that the imposition of salt stress caused a significant reduction in shoot dry weights, as was evident from the analysis of variance (Table 1). Maximum dry weights under control conditions were observed in Dunkled (5.50 g), while minimum in Legend and Oscar. Under the influence of the saline environment, maximum dry weight was produced in Faisal Canola, DGL, Dunkled, and CON-II (4.19, 4.35, 4.61, and 4.21 g/plant, respectively) and statistically ranked same as the apparent differences had no statistical significance. Legend and Oscar were at the bottom end (Figure 1).

2.1.3. Osmotic Potential (−MPa)

Leaf osmotic potential decreased in the plants of all canola cultivars when grown under saline condition (Table 1; Figure 2). A maximum decrease in leaf osmotic potential was found in salt-tolerant cv Dunkled, whereas the lowest decrease in leaf osmotic potential (more negative) was observed in salt-sensitive canola cultivars—Legend, Cyclone, AC-Excel, and Oscar.

2.1.4. Leaf Mineral Nutrients (Leaf Na+, K+, Ca2+, Cl, K+/Na+ Ratio, Ca2+/Na+ Ratio)

Leaf Na+ and leaf Cl markedly increased in salt-stressed plants of all canola cultivars, whereas the accumulation of K+ and Ca2+ in the leaves decreased significantly due to the imposition of salt stress (Table 1; Figure 3). Moreover, the accumulation of Na+ in the leaves was lower in all four salt-tolerant cultivars than in all four salt-sensitive canola cultivars. Similarly, salt-tolerant canola cultivars had greater leaf K+ under saline conditions. A maximum increase in leaf K+ was observed in DGL and Dunkled, and minimum K+ uptake was observed in the cultivar Oscar (Figure 3). Salt-tolerant cv. Dunkled had significantly greater shoot Ca2+ contents under saline conditions as compared to those of all other cultivars (Figure 3). Whereas, salt-sensitive cultivars—AC-Excel and Oscar—showed minimum shoot Ca2+ contents. Maximum shoot Cl contents were observed in salt-sensitive cultivars—Legend, Cyclone, AC-Excel, and Oscar (Figure 3). In contrast, the minimum increase in shoot Cl contents was recorded in salt-sensitive cultivar Faisal Canola under salt stress conditions.

2.1.5. Shoot K+/Na+ and Ca2+/Na+ Ratio

Salt stress decreased the shoot K+/Na+ ratio in all canola cultivars (Table 1). The minimum K+/Na+ ratio was observed in salt-sensitive cultivar Oscar, while the highest K+/Na+ ratio was observed in salt-tolerant cultivar Dunkled. In contrast, shoot Ca2+/Na+ ratio remained the same in salt-tolerant cultivars, whereas it increased in salt-stressed plants of salt-sensitive cultivars of canola (Table 1; Figure 4).

2.1.6. Nutrient Utilization Efficiency (mg2/µg)

It was evident from the statistical analysis of data regarding nutrient utilization efficiency that the imposition of salt stress in the rooting medium had a significant effect on nutrient utilization efficiency of canola cultivars (Table 1). Nutrient utilization efficiency decreased significantly in all the cultivars on the imposition of salt stress in the rooting medium (Figure 4). Under saline conditions, maximum nutrient utilization efficiency was observed in cultivars DGL and Dunkled. However, under saline conditions, minimum nutrient utilization efficiency was recorded in cultivar Oscar.

2.2. RNASeq Analysis and Differential Expression of Genes in Salt-Tolerant Cultivar Dunkled

On average, the control sample produced 58.97 million (58,969,962) reads, and total read bases were 8.9 G bp. The GC content was 49.15%, and Q30 was 94.06%. For the stress sample, 45.04 million (45,035,200) reads were produced with a total read base of 6.8 G bp. The GC content was 47.87%, and Q30 was 93.99%. Box plot showing raw expression (log2) is presented in (Figure S1). For the construction of contigs, raw reads were processed, as described in the Materials and Methods section. A total of 76,181 gene IDs were mapped on the reference genome (Table S1a). Differentially expressed genes (DEGs) were analyzed using R software. The upregulated genes were observed to be 29,187, and salt-stressed-induced downregulated genes were observed to be 29,291 (Table S1b). The analysis using DAVID revealed 17,912 characterized upregulated DEGs in control, and 20,931 characterized downregulated DEGs in the saline sample. Significant DEGs included sodium hydrogen exchanger (NHX), sodium transporter (HKT), potassium transporter (POT), Na-K-Cl co-transporter (NCKK1), cyclic nucleotide-gated ion channel 1 (CNGC1), mechanosensitive ion channel (MSL), potassium channel (KOR), chloride channel (ClCa), calmodulin (Calm), calmodulin binding transcriptional activator (CBTA), calcium transporting ATPase (Ca-ATPase), vacuolar ATPase (V-ATPase), heat shock proteins (HSP), late embryogenesis abundant proteins (LEA), Fe-superoxide dismutase (Fe-SOD), Cu-superoxide dismutase (Cu-SOD), Mn-superoxide dismutase (Mn-SOD), catalase (CAT), peroxidase (POD), tonoplast intrinsic protein (TIP), plasma-membrane intrinsic protein (PIP), nucleoplasm intrinsic protein (NIP), expansins (EXP), cell wall integrity and stress response component (WSC), NAC domain containing transcription factor (NAC), ethylene responsive transcription factor (ERF), MYB domain containing transcription factor (MYB), bZIP transcription factor (bZIP), 1-aminocyclopropane-1-carboxylic acid transcription factor (ACC), heat shock transcription factor (HSP) (Table S2a,b).
After selecting DEGs with fold change >5 and −5 and excluding DEGs other than protein-coding, DEGs were selected for gene function enrichment analysis by g:Profiler (https://biit.cs.ut.ee/gprofiler/gost). A bubble plot of statistically significant enriched terms is presented in Figure 5.
For upregulated DEGs, significant GO terms of the “Molecular Function” category included antioxidant activity (GO:0016209), DNA binding (GO:0003677), and peroxidase activity (GO:0004601). Significant GO terms of the “Biological Process” category included regulation of cellular process (GO:0050794), regulation of biological process (GO:0050789), homeostatic process (GO:0042592), hydrogen peroxide metabolic process (GO:0042743), cellular detoxification (GO:1990748), regulation of gene expression (GO:0010468), response to oxidative stress (GO:0006979), transmembrane transport (GO:0055085), and potassium ion transmembrane transport (GO:0071805). Significant GO terms of the “Cellular Component” category included intracellular membrane-bounded organelle (GO:0043231), cellular anatomical entity (GO:0110165), and cellular component (GO:0005575) (Table 2). For downregulated DEGs, significant GO terms of the molecular function category included MAP kinase activity (GO:0004707) and proline dehydrogenase activity (GO:0004657). Significant GO terms of the “Biological Process” category included MAPK cascade (GO:0000165), intracellular signal transduction (GO:0035556), and proline metabolic process (GO:0006560). Significant GO terms of the cellular component category included photosynthetic membrane (GO:0034357) and membrane protein complex (GO:0098796) (Figure 5).
Annotation by BlastKOALA resulted in 30 distinct terms with unique K numbers that were related to various transporter proteins, channels, and pumps (Table 2). It was inferred from these findings that salt stress caused over-expressed membrane-bound specialized proteins causing uptake of Na+, K+, Ca2+, and Cl ions, signal transduction, and triggering of differential gene expression by activation of various transcription factors (Table S4).

2.3. qRT-PCR

Transcriptome data was validated through qRT-PCR analysis. Among the upregulated genes, homologous genes of AtNKX and AtHKT1 were observed to be over-expressed in salt-stressed Dunkled plants when exposed to 200 mM NaCl (Figure 6).

3. Discussion

Plant growth responses to salt stress differ widely among glycophytes, reflecting differences in their genetic make-ups and respective potentials to withstand various stress levels [10]. Canola (Brassica napus L.) is moderately salt-tolerant with wide interspecific variation [1,21]. In the first part of the present study, four salt-tolerant and four salt-sensitive canola cultivars were used. Vegetative growth and development are some of the potential determinants of salt tolerance [22]. Plant growth results verified the degree of salt tolerance in canola cultivars examined in the present study. The existence of genetic diversity for abiotic stress tolerance in any crop species is prime for developing stress-tolerant cultivars through breeding programs [3]. Thus, studying physiological and biochemical responses in cultivars with diverse genetic background help in understanding the detailed mechanism of stress tolerance as well [10]. Although various physiological and biochemical processes are involved in salt stress tolerance, their contribution varies with species and type of cultivar. The well-known mechanisms can be categorized as osmotic tolerance, ion exclusion, and tissue tolerance [3,10,23,24]. However, several researchers have suggested that physiological traits contributing to these mechanisms should be evaluated individually as the ability of plants to maintain them under saline conditions [6,25,26,27]. Osmoregulation is an important phenomenon that is essential for normal cellular metabolism [28]. Under saline conditions, plants can accumulate organic and/or inorganic solutes, e.g., Glycine betain, proline, soluble sugars, Na+, K+, etc., to decrease its leaf osmotic potential, thereby increasing leaf turgor potential. Thus, if we establish a relationship between leaf relative water contents and osmotic potential, it becomes clear that leaf osmotic potential decreases with the decrease in relative water contents in plants grown under salt stress. Hence, an increase in turgor pressure is associated with a decrease in osmotic potential [29]. However, turgor by itself has no direct control on plant growth, rather it is essential for plant growth, as it is employed as an extending force required for the expansion of the cell wall [26]. In the present study, the leaf osmotic potential of all the cultivars decreased under the influence of saline conditions in the rooting medium. In plants grown in saline solution, lowest osmotic potential (more negative) was observed in Dunkled, while the highest was observed in Legend, Cyclone, AC-Excel, and Oscar. It is suggested that salt-tolerant canola cultivars maintained their water status through osmotic adjustment (either through the accumulation of organic solutes or inorganic osmotica) [26,30,31]. Leaf relative water contents of all the cultivars decreased under salinity stress. There was a marked increase in shoot Na+ and a decrease in shoot K+ contents in all canola cultivars under study, which resulted in a decrease in shoot K+/Na+ ratio. These results suggested that the accumulation of Na+ or K+ was not the main contributor to osmotic adjustment [24]. However, the K+/Na+ ratio was observed to show a positive relationship with salt tolerance, suggesting maintenance of ion homeostasis is more important in this regard [7].
Selection for salt-tolerant cultivars or identification of the mechanism of salt tolerance has been carried out through a variety of approaches, such as metabolite analysis, the gene for these metabolites through transcriptome analysis, miRNAs regulating genes responsible for metabolites and metabolic pathways, proteins responsible for specific pathways, or metabolites through proteome analysis, etc. [32,33,34,35]. Transcriptome studies to explore the mechanism of salt tolerance is gaining significant ground [16,17,34]. Thus, leaf transcriptome of salt-tolerant canola cultivar Dunkled was analyzed. Though roots are directly damaged by salt stress, growth inhibition of leaves in salt stress plants is more than that of roots. In the present study, differentially expressed genes (DEGs) were mainly involved in the regulation of ionic concentration in salt-tolerant cultivar Dunkled. Potassium transporters and channels upregulated under saline conditions, including two HKT1, two potassium transporter 9-like, one potassium transporter 4-like, one potassium channel KOR1 like, and one sodium/potassium/calcium exchanger. Greater salt tolerance in cv Dunkled was rendered to upregulation of transport proteins associated with sodium exclusion and preferential uptake of potassium. It was mediated either through adjusting K+ homeostasis by inducing the expression of potassium transporters or through vacuolar sodium sequestration by NHX, which is supported by the fact that salt stress upregulates vacuolar ATPases [34]. High Na+ concentration is reported to block high-affinity K+ transporters that result in reduced K+ influx and elevated Na+ influx [7,14,34].
Increased uptake of Ca2+ in the shoot was also reported to occur under salt stress, as was earlier observed in barley [36]. Nine calcium-sensing proteins were also upregulated, including calcium-sensing receptor chloroplastic like, calmodulin-like, and calmodulin-binding transcription activator proteins. Various reports have suggested that calcium sensory proteins activate sodium exclusion by activating plasma membrane or vacuolar Na+/H+ antiporters [14,24]. Downstream to calcium signaling MAP kinases may transduce the salt-induced osmotic stress signal and activate transcription factors by binding calmodulin and calmodulin-binding transcriptional activators [37]. Thus, the DEGs, like NHX, HKT, POT, NCKK1, CNGC1, MSL, KOR, ClCa, Calm, CBTA, Ca-ATPase, V-ATPase, HSPs, LEA, Fe-SOD, Cu-SOD, Mn-SOD, CAT, POD, TIP, PIP, NIP, EXP, WSC, NAC, ERF, MYB, bZIP, ACC, HSP, were inferred. Over-expression of AtNKX and AtHKT1 using qRT-PCR in Dunkled under salt stress validated that sodium exclusion and maintenance of potassium in leaf were mainly regulated through selective channels and transporters.

4. Materials and Methods

4.1. Evaluation of Different Cultivars of Canola (B. napus L.) for Salt Stress Tolerance

Canola (Brassica napus L.) plants from four salt-tolerant cultivars and four salt-sensitive cultivars were grown (25 plants/pot) in pots filled with sand. The sand was air-dried for 72 h, cleaned from straws, stones, and passed through 5 mm sieve in order to have uniformity. Pots were filled with 4.75 kg of sand. Thinning was done 5 days after completion of germination, leaving 6 plants/pot. Two-week-old plants were divided into two treatments, i.e., control (0 mM NaCl) and saline (200 mM NaCl). The salt solution was applied to target plant pots in such a quantity that the solution being eluted from the pot had the same EC (Electrical conductivity) as that of the irrigated salt solution, in order to attain a reliable stress level. Plants’ responses were observed and recorded for growth and physiological parameters after 5 weeks from the application of stress.

4.1.1. Growth Attributes

Plants were harvested after five weeks of salt stress and separated into shoots and roots. Fresh weights of shoots and roots were recorded. In addition, shoot and root lengths were also recorded.

4.1.2. Relative Water Contents (%)

Fully expanded youngest leaf (mostly 3rd leaf from the top) was cut from each replicate (one leaf per plant). Leaves were numbered, weighed, and the respective fresh weights were recorded. Leaves were transferred to a tray filled with water and kept for 60 min. Afterward, leaves were taken out, surface blotted, and weighed for turgid weights. Leaves were then kept in the oven for seven days at 65 °C, and dry weights were noted. Relative water contents were measured as follows:
R W C = ( F r e s h   w e i g h t D r y   w e i g h t ) × 100 ( T u r g i d   w e i g h t D r y   w e i g h t )

4.1.3. Leaf Osmotic Potential (−MPa)

The third leaf from the top was excised from each replicate. Leaves were stored at −20 °C for a week. Frozen leaf material was allowed to thaw and thoroughly pressed with a glass rod to extract sap. The osmotic potential was determined using a vapor pressure osmometer.

4.1.4. Mineral Contents (K+, Na+, Cl, Ca2+)

The 0.1 g of powdered dry shoot material was digested using H2SO4 and H2O2, according to the method of [38]. The amount of Na+, K+, and Ca2+ ions in the digested shoots and roots were measured using a flame photometer (PFP 7, Jenway, UK). However, Cl ions were extracted in distilled water (10 mL). The mixture was continuously heated (80 °C) until the volume of the mixture was reduced to half. The volume of the mixture was restored to 10 mL with distilled water. Chloride meter was used to determine Cl contents in the samples.

4.1.5. Nutrient Utilization Efficiency (mg2/µg)

Nutrient utilization efficiency was determined on the basis of shoot dry matter of the plant, following [39]:
N u t r i e n t   u t i l i z a t i o n   e f f i c i e n c y   ( m g 2 µ g ) = 1 × s h o o t   d r y   w e i g h t   ( m g ) S h o o t   n u t r i e n t   c o n c e n t r a t i o n   ( µ g / m g )

4.2. RNASeq Analysis and Differential Expression of Genes in Salt-Tolerant Cultivar Dunkled

4.2.1. Isolation of Total RNA

Samples of leaves were harvested after 24 h of salt stress, immediately immersed in liquid nitrogen, and stored in a plastic zipper bag at −80 °C. Later on, leaves were ground in liquid nitrogen, and total RNA was isolated by an optimized protocol based on the hot borate method by Wan and Wilkins [40].

4.2.2. Next-Generation Sequencing (NGS)

RNA was subjected to next-generation sequencing (Mcrogen, Seoul, Korea). Quality control (QC) was performed, and qualified samples proceeded to library construction. TruSeq Stranded mRNA LT Sample Prep Kit (Illumina, San Diego, CA, USA) was used using TruSeq stranded mRNA sample preparation guide, part # 15031047 Rev. E. The sequencing library was prepared by random fragmentation of the cDNA samples, followed by 5′ and 3′ adapter ligation. Alternatively, “tagmentation” combined the fragmentation and ligation reactions into a single step to increase the efficiency of the library preparation process. Adapter-ligated fragments were then PCR amplified and purified.

4.2.3. Sequencing for Cluster Generation

The library was loaded into a flow cell where fragments were captured on a lawn of surface-bound oligos complementary to the library adapters. Each fragment was then amplified into distinct, clonal clusters through bridge amplification. When cluster generation was complete, the templates were ready for sequencing. Illumina SBS (sequencing by synthesis) technology utilized a proprietary reversible terminator-based method that detected single bases as they were incorporated into DNA template strands. As all 4 reversible, terminator-bound dNTPs were present during each sequencing cycle, natural competition minimized incorporation bias and greatly reduced raw error rates compared to other technologies. The result was highly accurate base-by-base sequencing that virtually eliminated sequence-context-specific errors, even within repetitive sequence regions and homopolymers.

4.2.4. Generation of Raw Data

The Illumina sequencer generated raw images, utilizing sequencing control software for system control and base calling through integrated primary analysis software called RTA (real-time analysis). The BCL (base calls) binary was converted into FASTQ, utilizing Illumina package bcl2fastq.

4.2.5. Bioinformatics Analysis of Raw Data

Data quality was assessed using the FastQC tool. Data contained “Nextera Transpose Sequence” adaptor content. Adaptor sequences were removed using fastp [41] software. High-quality reads (contigs) were assembled using TRINITY software. For functional annotation, assembled contigs were annotated using reference Brassica napus transcriptome data (https://0-www-ncbi-nlm-nih-gov.brum.beds.ac.uk/assembly/GCF_000686985.2/) through BWA [42] software (parameters: mem -t 4 -k 32 -M). Transcript counts were calculated using feature-count [43] software using -gene parameter.

4.2.6. Bioinformatics Analysis of Feature Count Matrix

Differential expression of genes was determined using Bioconductor’s edgeR [44] package in R software (version: 3.5.3), yielding upregulated and downregulated DEGs. DEGs were converted to DAVID IDs online (https://david.ncifcrf.gov/conversion.jsp) using the gene ID conversion tool, selecting the official gene symbol as an identifier and Brassica napus as selected species. Only protein-coding DEGs were selected for further analysis. Selected DEGs with fold change (log2FC) greater than 5 (for upregulated DEGs) and −5 (for downregulated DEGs) were subjected to functional enrichment analysis of significant GO terms using g:Profiler (version e99_eg46_p14_f929183). Benjamini–Hochberg false discovery rate (FDR) threshold was set to 0.05. A bubble plot was made with the z-score (x-axis) against the negative logarithm of the adjusted p-value (−log Padj) (y-axis). Further annotation was carried out by KOALA (KEGG Orthology And Links Annotation) of KEGG GENES using BlastKOALA (https://www.kegg.jp/blastkoala/), and K numbers were assigned to the respective amino acid sequence data.

4.2.7. Validation of NGS Data by qRT-PCR

The expression of two genes, i.e., AtNKX and AtHKT1, was assessed using RT-PCR. These genes were selected on the basis of having a documented role in salt tolerance mechanism [7,14,34]. RNA was converted to cDNA using Vivantis cDSK01-050 cDNA synthesis kit (Vivantis Technologies, 40170 Shah Alam, Selangor Darul Ehsan, Malaysia). Real-time polymerase chain reactions were performed on Mic PCR (Bio Molecular Systems, Brisbane Queensland, Australia). Beta-ACTIN was used as an internal control. Primer sequences are given in Table S5.

5. Conclusions

Salt tolerance in salt-tolerant canola cultivars—Faisal Canola, DGL, Dunkled, and CON-II—was found to be associated with maintenance of plant water status and salt exclusion. Salt-tolerant canola cultivars were better able to osmotically adjust and exclude toxic salts. Moreover, RNA-seq data verified that DEGs in salt-tolerant canola cultivar Dunkled were mainly related to membrane proteins related to water and ion transport (Na+, K+, Ca2+, Cl).
It is thus concluded that salt stress upregulated various ion transporters while downregulated regulatory genes of signaling pathways. Among upregulated ion transporters, vacuolar and plasma membrane Na+/H+ antiporters, potassium transporters and channels (HKT1, KOR), calcium transporters, ABC transporters, cyclic nucleotide-gated ion channels (CNGCs), mechanosensitive ion channels, and chloride channels were included. Among the downregulated genes, transcripts related to MAPKs and senescence-associated genes were included. Therefore, enhanced salt stress tolerance observed in canola cultivar could be associated with water uptake, ion exclusion, and osmotic adjustment. The overexpression results of qRT-PCR with documented homologous Arabidopsis thaliana genes (AtNHX and AtHKT1) in salt-tolerant canola cultivar Dunkled further confirmed our findings.

Supplementary Materials

The following are available online at https://0-www-mdpi-com.brum.beds.ac.uk/2223-7747/9/7/891/s1, Figure S1: Box plot showing raw expression (Log2 scale) for control (0 mM NaCl) and saline (200 mM NaCl) plant samples from Brassica napus L. cultivar Dunkled, Table S1a: Gene count data for RNASeq from Canola (Brassica napus L.) plants grown under Control (0 mmol NaCl) or Saline (200 mmol NaCl) conditions, Table S1b: Fold Change (log2FC) for DEGs in Canola (Brassica napus L.) plants grown under Control (0 mmol NaCl) or Saline (200 mmol NaCl) conditions, Table S2a: List of selective upregulated DEGs (Fold change >5), Table S2b: List of selective downregulated DEGs (Fold change > −5), Table S3a: Functional enrichment analysis of selective upregulated DEGs (Fold change >5), Table S3b: Functional enrichment analysis of selective downregulated DEGs (Fold change > −5), Table S4: Annotation data for ion transporters after BlastKOALA, Table S5: List of primers used for qRT-PCR.

Author Contributions

Conceptualization, H.-u.-R.A.; Formal analysis, M.U.; Investigation, M.U.; Methodology, M.U.; Supervision, H.-u.-R.A. and Z.-d.K.; Writing—original draft, M.U.; Writing—review and editing, H.M.K. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Athar, H.R.; Ashraf, M.; Wahid, A.; Jamil, A. Inducing salt tolerance in canola (Brassica napus L.) by exogenous application of glycinebetaine and proline: Response at the initial growth stages. Pak. J. Bot. 2009, 41, 1311–1319. [Google Scholar]
  2. Munns, R. Comparative physiology of salt and water stress. Plant Cell Environ. 2002, 25, 239–250. [Google Scholar] [CrossRef]
  3. Ashraf, M.; McNeilly, T. Salinity tolerance in brassica oilseeds. Crit. Rev. Plant Sci. 2004, 23, 157–174. [Google Scholar] [CrossRef]
  4. Ulfat, M.; Athar, H.; Ashraf, M.; Akram, N.A.; Jamil, A. Appraisal of physiological and biochemical selection criteria for evaluation of salt tolerance in canola (Brassica napus L.). Pak. J. Bot. 2007, 39, 1593–1608. [Google Scholar]
  5. Ashraf, M.; Athar, H.; Harris, P.; Kwon, T. Some prospective strategies for improving crop salt tolerance. Adv. Agron. 2008, 97, 45–110. [Google Scholar] [CrossRef]
  6. Morton, M.J.L.; Awlia, M.; Al-Tamimi, N.; Saade, S.; Pailles, Y.; Negrão, S.; Tester, M. Salt stress under the scalpel—Dissecting the genetics of salt tolerance. Plant J. 2019, 97, 148–163. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  7. Van Zelm, E.; Zhang, Y.; Testerink, C. Salt tolerance mechanisms of plants. Annu. Rev. Plant Biol. 2020, 71, 403–433. [Google Scholar] [CrossRef] [Green Version]
  8. Munns, R.; James, R.; Gilliham, M.; Flowers, T.J.; Colmer, T.D. Tissue tolerance: An essential but elusive trait for salt-tolerant crops. Funct. Plant Biol. 2016, 43, 1103. [Google Scholar] [CrossRef] [Green Version]
  9. Hasegawa, P.M.; Bressan, R.A.; Zhu, J.K.; Bohnert, H.J. Plant cellular and molecular responses to high salinity. Annu. Rev. Plant Physiol. Plant Mol. Biol. 2000, 51, 463–499. [Google Scholar] [CrossRef] [Green Version]
  10. Munns, R.; Tester, M. Mechanisms of salinity tolerance. Annu. Rev. Plant Biol. 2008, 59, 651–681. [Google Scholar] [CrossRef] [Green Version]
  11. Munns, R. Genes and salt tolerance: Bringing them together. New Phytol. 2005, 167, 645–663. [Google Scholar] [CrossRef]
  12. Ishitani, M.; Liu, J.P.; Halfter, U.; Kim, C.S.; Shi, W.M.; Zhu, J.K. SOS3 function in plant salt tolerance requires N-myristoylation and calcium binding. Plant Cell 2000, 12, 1667–1677. [Google Scholar] [CrossRef] [Green Version]
  13. Zhu, J.-K. Genetic analysis of plant salt tolerance using Arabidopsis. Plant Physiol. 2000, 124, 941–948. [Google Scholar] [CrossRef] [Green Version]
  14. Ismail, A.M.; Horie, T. Genomics, physiology, and molecular breeding approaches for improving salt tolerance. Annu. Rev. Plant Biol. 2017, 68, 405–434. [Google Scholar] [CrossRef] [Green Version]
  15. Guo, J.; Shi, G.; Guo, X.; Zhang, L.; Xu, W.; Wang, Y.; Su, Z.; Hua, J. Transcriptome analysis reveals that distinct metabolic pathways operate in salt-tolerant and salt-sensitive upland cotton varieties subjected to salinity stress. Plant Sci. 2015, 238, 33–45. [Google Scholar] [CrossRef]
  16. Long, W.; Zou, X.; Zhang, X. Transcriptome analysis of canola (Brassica napus) under salt stress at the germination stage. PLoS ONE 2015, 10, e0116217. [Google Scholar] [CrossRef]
  17. Lin, J.; Li, J.; Yuan, F.; Yang, Z.; Wang, B.; Chen, M. Transcriptome profiling of genes involved in photosynthesis in Elaeagnus angustifolia L. under salt stress. Photosynthesis 2018, 56, 998–1009. [Google Scholar] [CrossRef]
  18. Iqbal, M.; Athar, H.-U.-R.; Ibrahim, M.; Javed, M.; Zafar, Z.U.; Ashraf, A. Leaf proteome analysis signified that photosynthesis and antioxidants are key indicators of salinity tolerance in canola (Brassica napus L.). Pak. J. Bot. 2019, 51, 52. [Google Scholar] [CrossRef] [Green Version]
  19. Sicilia, A.; Testa, G.; Santoro, D.F.; Cosentino, S.L.; Piero, A.R.L. RNASeq analysis of giant cane reveals the leaf transcriptome dynamics under long-term salt stress. BMC Plant Biol. 2019, 19, 355. [Google Scholar] [CrossRef] [Green Version]
  20. Du, X.; Wang, G.; Ji, J.; Shi, L.; Guan, C.; Jin, C. Comparative transcriptome analysis of transcription factors in different maize varieties under salt stress conditions. Plant Growth Regul. 2016, 81, 183–195. [Google Scholar] [CrossRef]
  21. Ashraf, M.; Oztürk, M.A.; Athar, H.R. Salinity and Water Stress: Improving Crop Efficiency; Springer: Berlin, Germany, 2009; ISBN 978-1-4020-9065-3. [Google Scholar]
  22. Lauchli, A.; Grattan, S. Plant growth and development under salinity stress. In Advances in Molecular Breeding Toward Drought and Salt Tolerant Crops; Jenks, M.A., Hasegawa, P.M., Jain, S.M., Eds.; Springer: Dordrecht, The Netherlands, 2007; pp. 1–32. ISBN 978-1-4020-5578-2. [Google Scholar]
  23. Roy, S.; Chakraborty, U. Salt tolerance mechanisms in Salt Tolerant Grasses (STGs) and their prospects in cereal crop improvement. Bot. Stud. 2014, 55, 31. [Google Scholar] [CrossRef] [Green Version]
  24. Roy, S.; Negrão, S.; Tester, M. Salt resistant crop plants. Curr. Opin. Biotechnol. 2014, 26, 115–124. [Google Scholar] [CrossRef]
  25. Cuartero, J.; Romero-Aranda, R.; Yeo, A.; Flowers, T. Variability for some physiological characters affecting salt tolerance in tomato. Acta Hortic. 2002, 573, 435–441. [Google Scholar] [CrossRef]
  26. Ashraf, M. Some important physiological selection criteria for salt tolerance in plants. Flora Morphol. Distrib. Funct. Ecol. Plants 2004, 199, 361–376. [Google Scholar] [CrossRef]
  27. Negrão, S.; Schmöckel, S.M.; Tester, M. Evaluating physiological responses of plants to salinity stress. Ann. Bot. 2016, 119, 1–11. [Google Scholar] [CrossRef] [Green Version]
  28. Taiz, L.; Zeiger, E.; Moller, I.S.; Murphy, A. Plant Physiology and Development; Sinauer Associates Inc.: Boston, MA, USA, 2015. [Google Scholar]
  29. Ali, Q.; Athar, H.R.; Ashraf, M. Ion transport in four canola cultivars as influenced by salt stress. Pak. J. Bot. 2006, 8, 1703–1708. [Google Scholar]
  30. Akram, M.S.; Ali, Q.; Athar, H.R.; Bhatti, A.S. Ion uptake and distribution in Panicum antidotale Retz. under salt stress. Pak. J. Bot. 2006, 38, 1661–1669. [Google Scholar]
  31. Munns, R.; Passioura, J.; Colmer, T.D.; Byrt, C. Osmotic adjustment and energy limitations to plant growth in saline soil. New Phytol. 2019, 225, 1091–1096. [Google Scholar] [CrossRef] [Green Version]
  32. Winicov, I.; Bastola, D.R. Transgenic overexpression of the transcription factor Alfin1 enhances expression of the endogenous MsPRP2 gene in Alfalfa and improves salinity tolerance of the plants. Plant Physiol. 1999, 120, 473–480. [Google Scholar] [CrossRef] [Green Version]
  33. Postnikova, O.A.; Shao, J.; Nemchinov, L.G. Analysis of the Alfalfa root transcriptome in response to salinity stress. Plant Cell Physiol. 2013, 54, 1041–1055. [Google Scholar] [CrossRef]
  34. Lei, Y.; Xu, Y.; Hettenhausen, C.; Lu, C.; Shen, G.; Zhang, C.; Li, J.; Song, J.; Lin, H.-H.; Wu, J. Comparative analysis of alfalfa (Medicago sativa L.) leaf transcriptomes reveals genotype-specific salt tolerance mechanisms. BMC Plant Biol. 2018, 18, 1–14. [Google Scholar] [CrossRef] [Green Version]
  35. Long, R.; Li, M.; Kang, J.; Zhang, T.; Sun, Y.; Yang, Q. Small RNA deep sequencing identifies novel and salt-stress-regulated microRNAs from roots of Medicago sativa and Medicago truncatula. In The Model Legume Medicago truncatula; Bruijn, F.D., Ed.; Wiley online Library: Hoboken, NJ, USA, 2020; pp. 963–976. [Google Scholar]
  36. Ageeva, I.F.; Prusakova, L.D.; Chizhova, S.I. Influence of brassinosteroids on stem formation and potassium and calcium ion contents in spring barley plants. Agrokhimiya 2001, 6, 49–55. [Google Scholar]
  37. Deinlein, U.; Stephan, A.B.; Horie, T.; Luo, W.; Xu, G.; I Schroeder, J. Plant salt-tolerance mechanisms. Trends Plant Sci. 2014, 19, 371–379. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  38. Wolf, B. A comprehensive system of leaf analyses and its use for diagnosing crop nutrient status. Commun. Soil Sci. Plant Anal. 1982, 13, 1035–1059. [Google Scholar] [CrossRef]
  39. Glass, A.D.M. Regulation of Ion Transport. Annu. Rev. Plant Physiol. 1983, 34, 311–326. [Google Scholar] [CrossRef]
  40. Wan, C.; Wilkins, T. A modified hot borate method significantly enhances the yield of high-quality RNA from cotton (Gossypium hirsutum L.). Anal. Biochem. 1994, 223, 7–12. [Google Scholar] [CrossRef]
  41. Chen, S.; Zhou, Y.; Chen, Y.; Gu, J. Fastp: An ultra-fast all-in-one FASTQ preprocessor. Bioinformatics 2018, 34, i884–i890. [Google Scholar] [CrossRef]
  42. Li, H.; Durbin, R. Fast and accurate short read alignment with Burrows-Wheeler transform. Bioinformatics 2009, 25, 1754–1760. [Google Scholar] [CrossRef] [Green Version]
  43. Liao, Y.; Smyth, G.K.; Shi, W. FeatureCounts: An efficient general purpose program for assigning sequence reads to genomic features. Bioinformatics 2013, 30, 923–930. [Google Scholar] [CrossRef] [Green Version]
  44. Robinson, M.D.; McCarthy, D.J.; Smyth, G.K. EdgeR: A Bioconductor package for differential expression analysis of digital gene expression data. Bioinformatics 2009, 26, 139–140. [Google Scholar] [CrossRef] [Green Version]
Figure 1. Shoot fresh (A) and dry weights (B) (g/plant) of eight cultivars of canola (Brassica napus L.) when two weeks old plants were grown under control (0 mM NaCl) or saline (200 mM NaCl) conditions for further five weeks (n = 3); Error bars are representing standard error.
Figure 1. Shoot fresh (A) and dry weights (B) (g/plant) of eight cultivars of canola (Brassica napus L.) when two weeks old plants were grown under control (0 mM NaCl) or saline (200 mM NaCl) conditions for further five weeks (n = 3); Error bars are representing standard error.
Plants 09 00891 g001
Figure 2. Relative water contents (A) (%) and osmotic potential (B) (−MPa) of eight cultivars of canola (Brassica napus L.) when two weeks old plants were grown under control (0 mM NaCl) or saline (200 mM NaCl) conditions for further five weeks (n = 3): Error bars are representing standard error.
Figure 2. Relative water contents (A) (%) and osmotic potential (B) (−MPa) of eight cultivars of canola (Brassica napus L.) when two weeks old plants were grown under control (0 mM NaCl) or saline (200 mM NaCl) conditions for further five weeks (n = 3): Error bars are representing standard error.
Plants 09 00891 g002
Figure 3. Shoot Na+ (A), K+ (B), Ca2+ (C), and Cl (D) contents (mg/g dry wt.) of eight cultivars of canola (Brassica napus L.) when two weeks old plants were grown under control (0 mM NaCl) or saline (200 mM NaCl) conditions for further five weeks (n = 3); Error bars are representing standard error.
Figure 3. Shoot Na+ (A), K+ (B), Ca2+ (C), and Cl (D) contents (mg/g dry wt.) of eight cultivars of canola (Brassica napus L.) when two weeks old plants were grown under control (0 mM NaCl) or saline (200 mM NaCl) conditions for further five weeks (n = 3); Error bars are representing standard error.
Plants 09 00891 g003
Figure 4. Shoot K+/Na+ ratio (A), Ca2+/Na+ ratio (B), and nutrient utilization efficiency (C) (mg2/µg) of eight cultivars of canola (Brassica napus L.) when two weeks old plants were grown under control (0 mM NaCl) or saline (200 mM NaCl) conditions for further five weeks (n = 3); Error bars are representing standard error.
Figure 4. Shoot K+/Na+ ratio (A), Ca2+/Na+ ratio (B), and nutrient utilization efficiency (C) (mg2/µg) of eight cultivars of canola (Brassica napus L.) when two weeks old plants were grown under control (0 mM NaCl) or saline (200 mM NaCl) conditions for further five weeks (n = 3); Error bars are representing standard error.
Plants 09 00891 g004
Figure 5. Bubble plot to represent functional profiling of significant gene ontology (GO) terms of selective upregulated (A) and downregulated (B) differentially expressed genes (DEGs) in canola (Brassica napus L.) grown under salt stress (Benjamini-Hochberg FDR threshold 0.05) using g:Profiler (version e99_eg46_p14_f929183). The x-axis represents z-score, whereas the y-axis represents the negative logarithm of the adjusted p-value. The area of the circles represents the number of genes assigned to the particular term. Detailed results are presented in Table S3a,b.
Figure 5. Bubble plot to represent functional profiling of significant gene ontology (GO) terms of selective upregulated (A) and downregulated (B) differentially expressed genes (DEGs) in canola (Brassica napus L.) grown under salt stress (Benjamini-Hochberg FDR threshold 0.05) using g:Profiler (version e99_eg46_p14_f929183). The x-axis represents z-score, whereas the y-axis represents the negative logarithm of the adjusted p-value. The area of the circles represents the number of genes assigned to the particular term. Detailed results are presented in Table S3a,b.
Plants 09 00891 g005
Figure 6. Expression (fold change) as depicted by qRT-PCR for two DEGs in leaf tissue of canola (Brassica napus L.) salt-tolerant cultivar Dunkled grown under saline (200 mM NaCl) conditions. Primer sequences are given in Table S5.
Figure 6. Expression (fold change) as depicted by qRT-PCR for two DEGs in leaf tissue of canola (Brassica napus L.) salt-tolerant cultivar Dunkled grown under saline (200 mM NaCl) conditions. Primer sequences are given in Table S5.
Plants 09 00891 g006
Table 1. Mean squares from the analysis of variance (ANOVA) of the data regarding eight cultivars of canola (Brassica napus L.) grown under control (0 mM NaCl) or saline (200 mM NaCl) conditions (n = 3).
Table 1. Mean squares from the analysis of variance (ANOVA) of the data regarding eight cultivars of canola (Brassica napus L.) grown under control (0 mM NaCl) or saline (200 mM NaCl) conditions (n = 3).
Source of VariationCultivar (df = 7)Salinity (df = 1)Cultivar × Salinity (df = 7)Error (df = 32)
Shoot fresh weight257.58 ***2254.50 ***26.85 ns11.75
Shoot dry weight0.79 **10.71 ***0.11 ns0.18
Osmotic potential0.16 ***1.35 ***0.07 ***0.01
Relative water contents15.44 ns685.35 ***32.28 *10.42
Shoot Na+ contents1.58 ns3579.38 ***36.46 ***6.12
Shoot K+ contents209.51 ***3291.80 ***45.73 ***7.57
Shoot Ca2+ contents21.42 ***643.15 ***2.15 ns3.85
Shoot Cl contents9.13 ns23919.01 ***28.71 **6.97
Shoot Na+/K+ ratio0.18 ***23.13 ***0.10 *0.03
Shoot Ca2+/Na+ ratio0.01 ns8.46 ***0.03 ns0.02
Nutrient utilization efficiency3.12e−5 *0.004 ***7.53e−5 ***1.28e−5
*** = significant at 0.001, ** = significant at 0.025, * = significant at 0.05, ns = non-significant. df = degree of freedom.
Table 2. KEGG (Kyoto Encyclopedia of Genes and Genomes) Orthology terms concluded from BlastKOALA for nucleotide sequences of the selective differentially expressed genes (DEGs) from plants of B. napus cultivar Dunkled grown under salt stress. (Detailed results are given in Table S4).
Table 2. KEGG (Kyoto Encyclopedia of Genes and Genomes) Orthology terms concluded from BlastKOALA for nucleotide sequences of the selective differentially expressed genes (DEGs) from plants of B. napus cultivar Dunkled grown under salt stress. (Detailed results are given in Table S4).
S.NoKEGG OrthologyK No.Protein IDTerm Size
1.calcium_sensing_receptor,_chloroplastic-likeK01013XP_013677299114
2.calcium-transporting_ATPase_10,_plasma_membrane-type-likeK01537XP_013668846 345
3.V-type_proton_ATPase_subunit_H-likeK02144XP_022557635 163
4.V-type_proton_ATPase_subunit_G2K02152XP_013674837 79
5.calmodulin-like_protein_12K02183XP_013736721 447
6.potassium_transporter_4-likeK03549XP_013689829 783
7.mitogen-activated_protein_kinase_19K04371XP_013644611603
8.probable_cyclic_nucleotide-gated_ion_channel_14K05391XP_013660045733
9.ABC_transporter_B_family_member_13-likeK05658XP_013650695328
10.ABC_transporter_C_family_member_7K05666XP_0136682551477
11.AP2-like_ethylene-responsive_transcription_factor_AIL1K09285XP_013746829457
12.ethylene-responsive_transcription_factor_ERF056-likeK09286XP_022545196151
13.dehydration-responsive_element-binding_protein_2B-likeK09287XP_022555207394
14.transcription_factor_MYB35-likeK09422XP_013655211310
15.aquaporin_TIP3-1K09873XP_013649603265
16.aquaporin_NIP6-1_XPK09874013725889305
17.V-type_proton_ATPase_subunit_G3K10604XP_013669950108
18.mitogen-activated_protein_kinase_kinase_4-likeK13413XP_013655549353
19.calmodulin-like_protein_8K13448XP_013710349153
20.sodium/potassium/calcium_exchanger_1-likeK13749XP_013717164288
21.abscisic_acid_receptor_PYL10-likeK14496XP_013741522222
22.mitogen-activated_protein_kinase_6-likeK14512XP_013656013392
23.mitogen-activated_protein_kinase_1K20535XP_013640933369
24.mitogen-activated_protein_kinase_3K20536NP_001303218370
25.mitogen-activated_protein_kinase_7K20537NP_001303162368
26.mitogen-activated_protein_kinase_kinase_9K20604XP_013648942308
27.mitogen-activated_protein_kinase_kinase_kinase_ANP1-like_isoform_X1K20606XP_013641030668
28.mitogen-activated_protein_kinase_kinase_kinase_18-likeK20716XP_013651099456
29.potassium_channel_KOR1-likeK21867XP_022549779632
30.mechanosensitive_ion_channel_protein_9-likeK22048XP_013659826738

Share and Cite

MDPI and ACS Style

Ulfat, M.; Athar, H.-u.-R.; Khan, Z.-d.; Kalaji, H.M. RNAseq Analysis Reveals Altered Expression of Key Ion Transporters Causing Differential Uptake of Selective Ions in Canola (Brassica napus L.) Grown under NaCl Stress. Plants 2020, 9, 891. https://0-doi-org.brum.beds.ac.uk/10.3390/plants9070891

AMA Style

Ulfat M, Athar H-u-R, Khan Z-d, Kalaji HM. RNAseq Analysis Reveals Altered Expression of Key Ion Transporters Causing Differential Uptake of Selective Ions in Canola (Brassica napus L.) Grown under NaCl Stress. Plants. 2020; 9(7):891. https://0-doi-org.brum.beds.ac.uk/10.3390/plants9070891

Chicago/Turabian Style

Ulfat, Mobina, Habib-ur-Rehman Athar, Zaheerud-din Khan, and Hazem M. Kalaji. 2020. "RNAseq Analysis Reveals Altered Expression of Key Ion Transporters Causing Differential Uptake of Selective Ions in Canola (Brassica napus L.) Grown under NaCl Stress" Plants 9, no. 7: 891. https://0-doi-org.brum.beds.ac.uk/10.3390/plants9070891

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop