Next Article in Journal
Intrinsically Connected: Therapeutically Targeting the Cathepsin Proteases and the Bcl-2 Family of Protein Substrates as Co-regulators of Apoptosis
Previous Article in Journal
The Novel Role of PGC1α in Bone Metabolism
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
IJMS
Volume 22
Issue 9
10.3390/ijms22094671
ijms-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Identification and Expression Analysis of SLAC/SLAH Gene Family in Brassica napus L.

by
Yunyou Nan
1,
Yuyu Xie
1,
Ayub Atif
1,
Xiaojun Wang
1,
Yanfeng Zhang
2,
Hui Tian
1,* and
Yajun Gao
1,3,*
1
College of Natural Resource and Environment, Northwest A&F University, Yangling 712100, China
2
Hybrid Rapeseed Research Center of Shaanxi Province, Yangling 712100, China
3
Key Laboratory of Plant Nutrition and the Agri-Environment in Northwest China, Ministry of Agriculture, Yangling 712100, China
*
Authors to whom correspondence should be addressed.
Int. J. Mol. Sci. 2021, 22(9), 4671; https://doi.org/10.3390/ijms22094671
Submission received: 4 March 2021 / Revised: 16 April 2021 / Accepted: 21 April 2021 / Published: 28 April 2021
(This article belongs to the Section Molecular Plant Sciences)
Browse Figures
Review Reports Versions Notes

Abstract

:
Slow type anion channels (SLAC/SLAHs) play important roles during anion transport, growth and development, abiotic stress responses and hormone responses in plants. However, there is few report on SLAC/SLAHs in rapeseed (Brassica napus). Genome-wide identification and expression analysis of SLAC/SLAH gene family members were performed in B. napus. A total of 23 SLAC/SLAH genes were identified in B. napus. Based on the structural characteristics and phylogenetic analysis of these members, the SLAC/SLAHs could be classified into three main groups. Transcriptome data demonstrated that BnSLAH3 genes were detected in various tissues of the rapeseed and could be up-regulated by low nitrate treatment in roots. BnSLAC/SLAHs were exclusively localized on the plasma membrane in transient expression of tobacco leaves. These results will increase our understanding of the evolution and expression of the SLAC/SLAHs and provide evidence for further research of biological functions of candidates in B. napus.

1. Introduction

Slow and rapid types of anion channels have been found in plants, which are expressed in guard cells and mediate slow and fast anion currents [1]. Slow type anion channel proteins (SLAC/SLAH) are one of the subfamilies of nitrate transporters. There are five SLAC/SLAH members in Arabidopsis (Arabidopsis thaliana) and they play an important role in stress signaling, growth and development, and hormonal response [2]. The interaction of AtSLAC1 and AtSLAH3 with different kinase phosphatases is associated with water stress signals [3,4,5]. In addition, SLAC1 is predominantly distributed in guard cells and is phosphorylated by Open stomata 1 (OST1) kinase, resulting in anion efflux from guard cells mediating stomatal closure [6,7]. In terms of hormonal signaling, abscisic acid (ABA) activates SLAC1 from direct interactions with OST1 and Calcium-dependent protein kinases (CPKs), and contributes to increase drought tolerance by regulating stomatal closure [5]. SLAC/SLAHs are also involved in the process of plant growth and development. For instance, SLAH3 phosphorylated by calcium-dependent protein kinase such as CPK2 and CPK20 in Arabidopsis regulated pollen tube growth [8]. In Arabidopsis, SLAC1 is mainly expressed in guard cells but weakly expressed in root/stamen/young silique, whereas SLAH3 is weakly expressed in guard cells and stamens, but strongly expressed in root [9]. SLAH2, the closest homologous protein to SLAH3, exclusively absorbs nitrate which is different from other SLAC/SLAH members with both absorption to nitrate and chloride [10]. SLAH1 is expressed in root, hypocotyl and stamen, and controls Cl− transport from root to shoot without affecting the absorption of NO3−, which makes SLAH1 a promising candidate gene for salt-tolerant plant biotechnology breeding [11,12]. SLAH4 is strongly expressed near the root tip, [13]. SLAC1 and SLAH3 have permeability ratios of 10 and 20 to NO3−/Cl−, respectively, when they are expressed in Xenopus oocytes [13]. SLAH3 shows a higher preference for nitrate than SLAC1, which is supposed to be a nitrate channel protein and functions in nitrate-dependent alleviation of ammonium toxicity [9]. The latest studies showed that sucrose non-fermenting 1-related protein kinase 1.1 (SnRK1.1) phosphorylated the C-terminal of SLAH3 and negatively regulated the process of SLAH3 involved in nitrate-dependent ammonium toxicity alleviation [14].
In rice (Oryza sativa), a total of 9 SLAC/SLAH members have been identified. Protein kinase OsSAPK8 is an important activator of OsSLAC1, and OsSLAC1 has nitrate selectivity [15]. In barley (Hordeum vulgare), HvSLAC1 requires nitrate to achieve ABA-induced stomatal closure. HvSLAC1 and other slow type anion channels in monocot have nitrate-dependent gating characteristics after the separation of monocot and dicot by the evolution of TMD3 series motifs in the trans-membrane domain [16]. ZmSLAC1 is a nitrate selective anion channel without obvious permeability to chloride, sulfate and malatein maize [17]. In pear (Pyrus bretschneideri), PbrSLAH3 are localized in the plasma membrane without expression in flower. Moreover, PbrSLAH3 has a strong selective absorption to nitrate and no permeability to chlorine, which is similar to ZmSLAC1. PbrSLAH3 regulated by PbrCPK32 is involved in the absorption and transport of nitrate in pear roots [18]. PttSLAH3 of poplar is not activated by protein kinase phosphorylation to absorb nitrate and chloride ions [19]. These studies indicate that SLAC/SLAHs has the nitrate selectivity and permeability and can regulate stomatal closure by regulating anion efflux from guard cells in response to ABA signals.
At present, more and more studies have focused on the slow anion channel proteins, especially those involved in the absorption and transport of nitrate. SLAC/SLAH family members have been identified and studied in several species such as Arabidopsis [9], rice [15], maize [17], barley [16], tobacco [20], poplar [19], and pear [21]. However, Brassica napus (B. napus), an important oil crop, has not been reported on SLAC/SLAHs. B. napus has a higher nitrogen requirement than other crops [22]. It was reported earlier that SLAC/SLAHs were highly expressed in roots and had the nitrate-permeable channel activity [9], we thus hypothesized that SLAC/SLAH family members in B. napus might be involved in nitrate uptake and transport. Therefore, the main objective of this study is to identify the SLAC/SLAH gene family of B. napus and analyze the structure, physicochemical properties, evolutionary relationships and expression profiles of corresponding genes through bioinformatics analysis. It is hoped to provide reference for further research on the biological functions of key candidate genes in nitrogen uptake and utilization of rapeseed.

2. Results

2.1. Identification and Classification of SLAC/SLAH Genes in B. napus

To identify the members of SLAC/SLAH family in B. napus, local protein BLAST using SLAC/SLAH protein sequences from Arabidopsis and Hidden Markov Model (HMM) search with conserved model (SLAC1.hmm) as query were performed. A total of 23 genes was identified in B. napus after several sequences redundant or without conserved domain (TDT_SLAC1_like) were removed (Table 1). Meanwhile, eleven and 12 SLAC/SLAH family members were identified in Brassica rapa (B. rapa) and Brassica oleracea (B. oleracea), respectively (Tables S1 and S2). The SLAC/SLAHs in B. napus were named according to the order of closest orthologues with Arabidopsis. B. napus SLAC/SLAH genes (hereafter termed as BnSLAC/BnSLAH) mainly included 2 to 6 exons encoding proteins of 329 to 640 residues and about 37.57 to 72.59 kDa molecular weight and 6 to 9.7 Isoelectric point (pI) value. These proteins were predicted to localize in plasma-membrane with 6 to 10 putative transmembrane domains (TMDs).

2.2. Phylogenetic Analysis of the SLAC/SLAHs between B. napus and Other Species and Identification of Evolutionary Selection Pressure on BnSLAC/SLAHs

To classify the SLAC/SLAHs and investigate their evolutionary relationships identified in B. napus, B. rapa and B. oleracea, phylogenetic analysis were performed to infer a functional relationship among them. A phylogenetic tree based on SLAC/SLAH protein sequences of A. thaliana, O. sativa, B. napus, B. rape and B. oleracea was constructed by the Neighbor-joining (NJ) method, which showed that the SLAC/SLAHs from the 5 species could be divided into three main groups, including SLAC, SLAH2/3 and SLAH1/4 (Figure 1). In phylogeny, group SLAH2/3 formed the largest cluster with 11 members in B. napus, 6 members in B. rapa, 5 members in B. oleracea, 2 members in A. thaliana and 6 members in O. sativa.
To characterize the selection pressure on the BnSLAC/SLAHs during the evolutionary process, we used the orthologous SLAC/SLAH gene pairs between B. napus and A. thaliana to determine the values of synonymous (Ks) and non-synonymous (Ka) nucleotide substitution rates, and Ka/Ks (Table 2). The Ka values of BnSLAC/SLAHs ranged from 0.0567 (BnSLAH2-1) to 0.2183 (BnSLAH4-5 and BnSLAH4-6) with an average of 0.1063, and the Ks values of BnSLAC/SLAHs ranged from 0.0677 (BnSLAH4-1) to 0.4910 (BnSLAH3-6) with an average of 0.3518. Further, we found that the Ka/Ks values of most BnSLAC/SLAHs were less than 1.0 except that of BnSLAH4 subfamily members were more than 1.0. Therefore, we presumed that the BnSLAC1s, BnSLAH1s, BnSLAH2s, BnSLAH3s might have experienced a very strong purify selection to preserve their function while BnSLAH4s were positively selected.

2.3. Chromosomal Location and Duplication Patterns

B. napus (AACC, 2n = 38) is an allotetraploid hybrid species of B. rapa (AA, 2n = 20) and B. oleracea (CC, 2n = 18) [23]. All the BnSLAC/SLAH genes were unevenly distributed on the rest 13 of 19 chromosomes according to the position of the start and stop codons (Figure 2) except for chromosomes A03-A05, A8, A10 and C06. Among them, 9 BnSLAC/SLAH genes were mapped to the chromosomes of A genomes while 14 members were found in the C genomes. Chromosome A09 and C01 included a maximum of 3 BnSLAC/SLAH genes, respectively. However, Chromosome A02, A07, C02, C03 and C05 only had a single BnSLAC/SLAH member. Genomic duplications in plants could be either whole-genome duplications (WGD) as large scale or tandem and segmental duplications as small scale. Tandem duplications can form on the same chromosome while segmental duplications events often occur on the different chromosomes. The number of BnSLAC/SLAHs in B. napus is different from that in B. rapa and in B. oleracea genomes, indicating that the slow type anion channels in B. napus underwent expansion and evolution in comparison with its ancestor species, B. rapa and B. oleracea. This study identified several tandem and segmental duplications between some BnSLAC/SLAH gene pairs. BnSLAH4-5 and BnSLAH4-6 gene pair located on the chromosome C01 generated from tandem duplications was found whereas 2 BnSLAC1s (BnSLAC1-1 and BnSLAC1-2), 2 BnSLAH1s (BnSLAH1-2 and BnSLAH1-3), 2 BnSLAH2s (BnSLAH2-1 and BnSLAH2-2), 5 BnSLAH3s (BnSLAH3-2 to BnSLAH3-6) and 4 BnSLAH4s (BnSLAH4-1 to BnSLAH4-4) were identified to be segmental duplicated genes (Figure 2). These results indicated that segmental duplication events functioned as a major driving force for the expansion of SLAC/SLAHs in B. napus.

2.4. Gene Structure and Conserved Motif Analysis of the SLAC/SLAH Family Members in B. napus

The phylogenetic tree of BnSLAC/SLAH protein sequences illustrated the evolutionary relationships of these members in B. napus (Figure 3). All the coding sequences of the BnSLAC/SLAH genes were spaced by introns with numbers varying from 1 to 5. The number of introns and the lengths of individual exons had the high similarity within the same subfamily, except that BnSLAH4 genes had different numbers of exons ranging from 2 to 5 (Figure 3 middle). Motif analysis was performed for the most conserved 10 motifs sequences in 23 BnSLAC/SLAH proteins using MEME tool. Most of the closely related members in the phylogenetic tree shared common motif compositions (Figure 3 right). BnSLAH3 subfamily contained motif 1 to motif 10 while BnSLAH2 mainly had motif 1-motif 5 and motif 7. Motif 1 to motif 7 were widely distributed in the BnSLAH1, BnSLAH4 and BnSLAC1 subfamilies. These conserved motifs are considered to possibly have functional and/or structural roles in active proteins.

2.5. Promoter cis-Acting Element Analysis and Interaction Protein Prediction

Attempting to understand the transcriptional regulation mechanisms of the identified BnSLAC/SLAH genes, it was essential to have insights about the upstream promoter region (Ravel et al. 2015). Two-kb upstream region from the initiation codon of BnSLAC/SLAH genes was obtained and supplied to the Plant CARE database to investigate of cis-regulatory elements. A total of 42 different cis-elements associated with light responsive, stress responsive, the phyto-hormone responsive and growth regulation have been identified in upstream regions of 23 BnSLAC/SLAH genes. A heatmap was constructed based on presence or absence of regulatory elements in each corresponding gene to better characterize the large number of motifs (Figure 4). These results indicated that complex regulatory networks may be implicated in the transcriptional regulation of BnSLAC/SLAHs. Cis-regulatory elements, CAAT-box and TATA-box, were commonly shared by all BnSLAC/SLAH genes. Obviously, CAAT- and TATA boxes are two common cis-regulatory elements in upstream regions of eukaryotic genes. Box4, G-box and TCT-motif elements responding to light existed in the 2-kb upstream region of more than approximately 70% BnSLAC/SLAH genes. Most BnSLAC/SLAH genes contained ABRE element (ABA responsive) except for BnSLAH2-1, BnSLAH2-3 and BnSLAH3-2. BnSLAH2 and BnSLAH3 subfamily genes except for BnSLAH2-1 contains TGACG-motif element involved in MeJA response. Moreover, eleven BnSLAC/SLAHs harbored drought responsive cis-elements (MBS) while 8 BnSLAC/SLAHs contained low-temperature responsive cis-element (LTR).
To further identify the proteins potentially interacting with the SLAC/SLAH family members, we constructed a protein interaction network involving direct (physical) and indirect (function) association by using the STRING database based on either known experimental or predicted interactions. All the SLAC/SLAH proteins consistently interacted with OST1. Besides, CPK6, CPK 21 and CPK23 were also observed to interact with the SLAC/SLAH proteins (Figure S1). Interestingly, BnSLAH1 and BnSLAH3 were predicted to interact with nitrate excretion transporter 1 (NAXT1) which was one of the nitrate transporters.

2.6. Synteny Analysis

In the present study, we also identified 11 and 12 SLAC/SLAH genes in the B. rapa and B. oleracea genomes, respectively. To further infer the phylogenetic mechanisms of the BnSLAC/SLAHs, we constructed a comparative syntenic map of B. napus and its ancestors (Arabidopsis, B. rapa and B. oleracea). Collinearity analysis revealed that there were strong orthologs of SLAC/SLAH genes between B. napus and the other three ancestral species (Figure 5). Nine and 4 of the SLAC/SLAH genes in the A subgenome of B. napus showed syntenic relationships with 9 and 4 SLAC/SLAH genes in the B. rapa and Arabidopsis genomes, respectively. In contrast, eight and 2 of the genes in the B. napus C subgenome were syntenic with 9 and 2 of the B. oleracea and Arabidopsis genomes, respectively. Additionally, 10 pairs of BnSLAC/SLAH genes are paralogs to each other in B. napus. The fact that nearly all of the homologous BrSLAC/SLAHs and BoSLAC/SLAHs maintained a syntenic relationship with BnSLAC/SLAHs suggested that whole-genome duplication (polyploidy) also played a major driving force for BnSLAC/SLAHs evolution besides segmental duplication.

2.7. Subcellular Localization of the SLAC/SLAH Genes of B. napus

All of the BnSLAC/SLAHs were predicted to localize on the plasma membrane. To confirm the above prediction, three BnSLAC/SLAH genes were cloned and fused to the N-terminus with the GFP protein driven by the CaMV 35S promoter. The fusion proteins (BnSLAC/SLAH-GFP) and control (empty vector) were transiently transformed into tobacco leaves. Microscopic visualization showed that green fluorescence was distributed throughout the whole cell when the empty vector was used. The green fluorescence was exclusively detected on the plasma membrane by confocal microscopic when the vectors contained BnSLAH1-1, BnSLAH3-2 and BnSLAH3–3 while that was detected both on the plasma membrane and nucleus in the control (Figure 6). Thus, the BnSLAC/SLAHs might have a similar subcellular localization pattern that they were all membrane protein.

2.8. Expression Level of the SLAC/SLAH Genes in B. napus

The expression level of genes usually affects their biological function. In order to better investigate the functions of the BnSLAC/SLAH genes, we analyzed the expression patterns of the BnSLAC/SLAH genes in 12 different tissues based on the RNA-seq data from the B. napus cv. ‘Zhongshuang 11′ (BioProject ID PRJNA394926) (Figure 7, Table S4). According to the results, several genes belong to the BnSLAH3, BnSLAH1 and BnSLAC1 subfamilies had a higher expression in most tissues than others, which shared similar expression patterns in specific tissue. Especially, BnSLAH3-2, BnSLAH3-3 and BnSLAH3-4 were highly expressed in the pistil, root, stamen, silique and leaf. The transcripts of the BnSLAH4 and BnSLAH2 subfamilies genes were weakly expressed in the 12 different tissues. Moreover, BnSLAH1-1 was specifically expressed in the stamen, which was similar to BnSLAH1-2.
To investigate the mechanisms of nitrate transport and uptake of slow anion channel in B. napus, we selected several BnSLAC/SLAH genes for further study. The expression of BnSLAC/SLAH genes in leaves and roots of rapeseed seedlings under high nitrate (HN 7.5 mM) and low nitrate (LN 0.19 mM) concentrations was determined as was shown in Figure 8. In leaves, the expression level of the BnSLAC1-1 and BnSLAH3-2 increased about 10 and 5 folds after 7 d LN treatment. In addition, the expression level of BnSLAH1-1 and BnSLAH3-3 evaluated about 2 to 4 folds after 3 h and 12 h LN treatment while that of BnSLAH3-4 had no obvious change. In roots, the expression of BnSLAC1-1 was increased 3.5 times at 24 h after LN treatment but not significantly changed at other treatment time points. The expression levels of BnSLAH1-1 were significantly up-regulated especially at 3 h after low nitrogen treatment. In addition, BnSLAH3-2 was significantly up-regulated 2 to 3 times under 12 h and 24 h low N treatment. Besides, 3 h and 12 h LN treatment induced a 2 to 4-fold up-regulation of BnSLAH3-3 and BnSLAH3-4 expression. Therefore, the slow anion channel proteins are able to respond to low nitrogen signal in a short time, which may be involved in nitrate uptake and transport in B. napus.
ABA is a phyto-hormone involved in a variety of abiotic stress response. By analyzing the Cis-acting elements of the BnSLAC/SLAH promoter, we found that most members contained ABA response elements (ABRE). Therefore, we detected the change of gene expression in rapeseed seedlings treated with 100 uM ABA (Figure 9). The results showed that the expressions of BnSLAC1-1 and BnSLAH3-4 genes were significantly decreased after 12 h ABA treatment in leaves. BnSLAH3-2 was up-regulated 3–4 folds at 3 h and 24 h after ABA treatment, while the expression of BnSLAH3-3 displayed 58-fold up-regulation at 12 h after ABA treatment in leaves. In roots, BnSLAC1-1 was up-regulated about 2.5 folds after 24 h ABA treatment, while BnSLAH1-1, BnSLAH3-2, BnSLAH3-3 and BnSLAH3-4 showed down-regulation trend after 3 h, 12 h and 24 h ABA treatment.

3. Discussion

In this study, 23, 11 and 12 SLAC/SLAH family members were identified in B. napus, B. rapa and B. oleracea respectively using 5 AtSLAC/SLAH sequences as queries. B. napus is a heterotetraploid formed by natural multiplication and interspecific hybridization between B. rapa and B. oleracea [23]. However, 9 BnSLAC/SLAH members are located in the A sub-genome of B. napus, which is less than the 11 BrSLAC/SLAH members identified in B. rapa. In addition, 14 BnSLAC/SLAH members are located in the C sub-genome of B. napus, which is more than 12 BoSLAC/SLAH members of B. oleracea. This may be due to the occurrence of genome-wide duplication events, tandem and segmental duplication events, resulting in partial gene deletions and increases during evolution. To explore the evolutionary process of BnSLAC/SLAHs, we investigated phylogeny, gene duplication events and syntenic relationships between B. napus and other species. These SLAC/SLAHs of Arabidopsis, rice, B. napus, B. rapa, and B. oleracea were divided into three groups, including SLAC1, SLAH2/3 and SLAH1/4, which was consistent with the classification of SLAC/SLAHs in Arabidopsis [9]. We found that expansion of SLAC/SLAH genes was attributed to WGD and segmental duplication event, and Ka/Ks analysis indicated that all the SLAC/SLAHs except BnSLAH4s were purified to evolve. However, previous study indicated that SLAC/SLAH developed along with the stoma evolution [24]. These results implies that SLAC/SLAH plays a crucial role in regulating stomatal movement. We found that BnSLAH3 members contain 5 exons and 10 conserved motifs, while BnSLAH1 genes only have 2 exons, which indicate BnSLAH3 may have functional diversity during growth and development in plants. SLAC/SLAH genes have a wide range of functions in plant growth and development. In addition, the result of promoter analysis showed that 2-kb upstream region of BnSLAC/SLAH genes contained different kinds of cis-elements involved in the light response, phyto-hormone response (including ABA, IAA, GA, MeJA and SA), drought response, low temperature response and growth regulation. BnSLAC/SLAH may function in the abiotic stress and growth regulation. AtSLAC1 encodes a plasma membrane-localized protein that is highly permeable to malate and chloride in Arabidopsis [25]. In this study, several BnSLAC/SLAH genes were selected for subcellular experiment. The observation of confocal microscopy suggested that BnSLAH1-1, BnSLAH3-2 and BnSLAH3-3 were localized on the plasma membrane, which is consistent with the results observed in Arabidopsis and pear as well as the prediction of subcellular localization in B. napus.
Five slow anion channel proteins, SLAC1 and its homologues SLAH1-SLAH4, were found in the model plant Arabidopsis, which were differentially expressed in leaves, stamens and roots [9]. RNA-seq results in B. napus showed that BnSLAH3-2, BnSLAH3-3 and BnSLAH3-4 genes were highly expressed in stamen, pistil, root, leaf, silique and sepal, while BnSLAC1-1 and BnSLAC1-2 genes showed a high expression level in the sepal, which is different from the expression patterns of AtSLAC1. This suggests that several BnSLAH3 genes may play an important role in the growth and development of rapeseed in comparison with other BnSLAC/SLAH members. Previous studies have shown that the plant slow anion channel proteins are involved in regulating the electrical conductivity of various anions such as chloride and nitrate [25], and function in stress signal, growth and development, and hormonal response [2]. Among them, more studies have focused on SLAC1 and SLAH3. AtSLAH3 was reported to be highly expressed in root besides leaves and stamens [9]. The absorption ratio of NO3−/Cl− of AtSLAH3 was 20 in Xenopus oocytes and showed a higher preference for NO3−, which was consequently considered as a nitrate channel protein [5,6,26]. According to the previous studies and RNA-seq analysis of different tissues in B. napus, several BnSLAH/SLAH genes were selected for qRT-PCR analysis. The results showed that BnSLAH3-2 and BnSLAH3-3 and BnSLAH3-4 in root could be up-regulated at 12 h after low nitrate treatment (0.19 mM), indicating that the above BnSLAH3 genes could respond to low nitrate stress in a short time and they might promote nitrate uptake and transport in rapeseed roots. However, There was a report suggesting that high concentration (64 mM) of nitrate induced the expression of SLAC/SLAH genes in pear [21]. It may be attributed to different species, nitrate concentration and treatment time. Accumulation of ABA in wheat induces the expression of NRT2 gene to promote nitrate uptake [27]. Exogenous ABA treatment can improve the yield and nitrogen use efficiency of B. napus [28]. BnSLAH3-3 were significantly upregulated about sixty times at 12 h in leaves while BnSLAH3-2, BnSLAH3-3 and BnSLAH3-3 were down-regulated in roots after ABA treatment. It is suggested that BnSLAH3s in leaves stimulated by exogenous ABA signals may be mainly involved in stomatal movement rather than nitrate uptake in roots.
AtSLAC1 and AtSLAH3 were activated by the calcium-independent protein kinase OST1 as well as by CPKs [29,30,31]. PbrSLAH3 showed a high selectivity for nitrate over chloride and was confirmed to interact with PbrCPK32 through yeast two-hybrid and bimolecular fluorescence complementation assays transport [18]. In addition, the prediction of protein interaction displayed that NAXT1, one of the NRT1 family members, was capable of interacting with SLAH3. Nitrogen is an essential nutrient element for plant. Adequate nitrogen uptake from the external environment is a crucial process for plant growth and development [26]. Rapeseed is the third largest oil crop worldwide [32]. B.napus, main cultivar of rapeseed, requires more nitrogen fertilizer with low nitrogen utilization efficiency [33,34]. Based on previous reports on the high expression level of SLAH3 gene in the root system, the strong permeable to nitrate and the interaction between SLAH3 and NAXT1, it was presumed that BnSLAH3s may be involved in the absorption of nitrate in B. napus. In this study, 23 BnSLAC/SLAHs were identified, and the results showed that BnSLAH3 genes were expressed in multiple tissues and could respond to low nitrate in a short time. Therefore, further study will be focus on the biological function of BnSLAH3 genes and the regulatory mechanism of nitrate absorption and utilization in B. napus.

4. Materials and Methods

4.1. Identification of SLAC/SLAH Gene Family Members in B. napus

To identify the SLAC/SLAH genes in B. napus, various database searches were conducted. The SLAC/SLAH protein sequences of Arabidopsis, including AtSLAC1 (AT1G12480), AtSLAH1 (AT1G62280), AtSLAH2 (AT4G27970), AtSLAH3 (AT5G24030) and AtSLAH4 (AT1G62262), were downloaded from the TAIR (http://www.arabidopsis.org/, accessed on 26 January 2021) [35]. In addition, the SLAC/SLAH protein sequences of rice, such as Os04g48530.1, Os01g43460.1, Os05g13320.1, Os01g28840.1, Os01g12680.1, Os07g08350.1, Os05g50770.2, Os01g14520.1 and Os05g18670.1, were downloaded from Phytozome (http://phytozome.jgi.doe.gov/pz/portal.html#, accessed on 26 January 2021). The protein sequences of B. napus, B. rapa and B. oleracea were download from NCBI (https://www.ncbi.nlm.nih.gov/datasets/, accessed on 26 January 2021). The candidate SLAC/SLAH genes in B. napus were identified by using protein BLAST (https://www.ncbi.nlm.nih.gov/, accessed on 26 January 2021) to search the local Brassica database and using Arabidopsis SLAC/SLAH protein sequences as queries and an E-value ≤ 1 × 10−20 as a threshold. The HMM profiles of the SLAC1.hmm downloaded from the HMM database (https://www.ebi.ac.uk/Tools/hmmer/, accessed on 26 January 2021) were applied to search for the SLAC1-type sequences in the local protein databases of B. napus by performing local HMMER search program with an E-value ≤ 1 × 10−20 as a threshold. Furthermore, all the obtained SLAC/SLAH protein sequences were analyzed against the Pfam database (http://pfam.xfam.org/, accessed on 26 January 2021) to verify the presence of SLAC1 domains (PF03595). Protein sequences lacking the SLAC1 domain or having E-values more than 1 × 10−20 were removed. The same method was used to identify the SLAC/SLAHs in B. rapa and B. oleracea.

4.2. Phylogenetic and Synteny Analysis of SLAC/SLAH

The phylogenetic trees were constructed based on the SLAC/SLAH protein sequences from Arabidopsis, O. sativa, B. rapa, B. oleracea and B. napus using the Neighbor-Joining (NJ) method in MEGA6.0 (http://www.megasoftware.net/, accessed on 26 January 2021) [36], and the bootstrap test was carried out with 1000 replicates. DnaSP (DNA Sequence Polymorphism) v5 [37] was used to calculate the ratio of the non-synonymous substitution rate (Ka) to the synonymous substitution rate (Ks) and the Ka/Ks value between paralogous gene pairs.
MCScanX was used to identify syntenic chains according to the homologous pairs as input [38]. The identification of whole-genome (WGD)/segmental, tandem, proximal and dispersed duplications in the SLAC/SLAH gene family were performed by downstream analysis tools in the MCScanX package. The results were displayed using TBtools [39].

4.3. Sequence Analysis of SLAC/SLAH Genes/Proteins in B. napus

ProtParam tool was used to calculate the Physicochemical features of SLAC/SLAHs (http://web.expasy.org/protparam/, accessed on 26 January 2021). The structures of the SLAC/SLAH genes was analyzed using GSDS server (gsds.cbi.pku.edu.cn/, accessed on 26 January 2021) [40]. Trans-membrane domains of SLAC/SLAHs were predicted by SCAMPI server (http://scampi.cbr.su.se/, accessed on 26 January 2021; Bernsel et al. 2008). The subcellular localization of SLAC/SLAH proteins in B. napus were predicted with the WoLF PSORT online tool (https://www.genscript.com/wolf-psort.html, accessed on 26 January 2021) [41]. Conserved motif of the SLAC/SLAH protein sequences were identified by MEME suite with parameters, where maximum number of motifs was 10, minimum motif width was 6 and maximum motif width was 50 (http://meme-suite.org/, accessed on 26 January 2021) [42]. The phylogenetic tree of SLAC/SLAH in B. napus was drawn with the neighbor-joining way by MEGA 6.0 [36].

4.4. Chromosomal Location and Duplication of SLAC/SLAH Genes in B. napus

Characterization of the chromosomal location of each SLAC/SLAH gene was acquired from NCBI (https://www.ncbi.nlm.nih.gov/datasets/, accessed on 26 January 2021). Based on the position information of SLAC/SLAH genes on the chromosomes, all of them were located on the A and C genome chromosomes. SLAC/SLAH genes were mapped on the B. napus chromosome using Mapchart. Genes located surrounding and parted by less than five genes can be as tandem duplications. The segmental duplications were analyzed by the Maher’s study [43].

4.5. Promoter cis-Acting Element Analysis and Interaction Protein Prediction

Two thousand bp upstream regions from the start codon site of SLAC/SLAH genes were retrieved from NCBI database and supplied to Plant CARE database for promoter analysis (http://bioinformatics.psb.ugent.be/webtools/plantcare/html/, accessed on 26 January 2021) [44]. STRING (Search Tool for Recurring Instances of Neighboring Genes) v 11.0 (https://string-db.org, accessed on 26 January 2021) web-server was selected to predict and display protein interaction networks of the SLAC/SLAH proteins in B. napus.

4.6. Subcellular Localization of the SLAC/SLAHs in B. napus

The full-length coding sequences without termination codon of the selected SLAC/SLAH genes of B. napus were isolated and linked into the pEGFP vector containing GFP reporter (saved in our laboratory). The competent cells of Escherichia coli (DH5α) and Agrobacterium (LBA4404) were used for the transformation of recombinants. Primers used for gene cloning and vector construction are shown in Supplementary Table S3. Agrobacterium-mediated transient expression in tobacco (Nicotiana benthamiana) leaves was performed as previously described [18]. Images were processed using the laser scanning confocal microscopy (Olympus FV3000, Tokyo, Japan).

4.7. Plant Material and Stress Treatment

According to the previous study [45], B. napus cv. “Zhongshuang 11” (ZS11) was planted and treated for qRT-PCR analysis. The seeds were sterilized and placed a wet filter paper for germination. Ten days later, seedlings were cultured in hydroponics with modified Hoagland’s nutrient solution (NO3−, 7.5 mM). After two weeks, they were used for low nitrogen treatment (NO3−, 0.19 mM) and ABA treatment (100 uM) under the condition of 16 h light, 8h dark and 5000 Lux illumination intensity at 25 °C. Roots and leaves were sampled after 3 h, 12 h, 24 h, 3 d and 7 d nitrogen treatment from high nitrogen level (NO3−, 7.5 mM) and low nitrogen level (NO3−, 0.19 mM). Similarly, Roots and leaves were sampled after 3 h, 12 h and 24 h treatment with or without 100 uM ABA. All samples were collected and stored in the −80 °C.

4.8. Gene Expression Analysis of SLAC/SLAHs in B. napus

Twelve different tissues and organs based on the RNA-seq data from the B. napus ‘Zhongshuang 11′ (BioProject ID PRJNA394926) were downloaded from NCBI-SRA database (https://www.ncbi.nlm.nih.gov/sra/, accessed on 26 January 2021) for further analysis. Total RNA was isolated from the collected tissues using a total RNA extracting kit (Bioteke, Beijing, China) according to the manufacturer’s instructions. Gel electrophoresis was used to determine the quality of RNA samples (Figure S2). Then 1 μg RNA of each sample was used to synthesize cDNA with the PrimeScript RT reagent Kit (Trans Gen, Beijing, China). The specific primers of the BnSLAC/SLAH genes and the housekeeping gene (Bnactin7) were designed using the Primer Premier 5.0 software (Supplementary Table S3). The NCBI-Primer BLAST against the rapeseed cDNA database and the melt curves (Figure S3) corresponding primers were examined to verify the specificity. The qRT-PCR assays were performed with three replicates. All of the real-time PCR reactions included 2 μL diluted cDNA, 200 nM of each primer, 2× SYBER GREEN Master Mix (Bioteke, Beijing, China) and sterile water, for a final volume of 20 μL. The following thermal cycle conditions were used with the Applied Biosystems (ABI Quantstudio 5, Foster city, CA, USA): pre-incubation at 94 °C for 2 min, then 45 cycles of 94 °C for 15 s, 60 °C for 15 s and 72 °C for 30 s, and with a default procedure for melt curve. Additionally, the relative expression levels were calculated with the 2−ΔΔCt method. Data was analyzed using the Office 2010 software, and statistical analyses were conducted with SPSS17.0 software using Student’s t test at the p < 0.05 level of significance.

5. Conclusions

In conclusion, 23 full-length SLAC/SLAH genes were identified in B. napus genomes. The subcellular localizations of BnSLAC/SLAHs provided evidence regarding their functions across the plasma membrane. Moreover, RNA-seq and qRT-PCR analysis showed that BnSLAH3-2, BnSLAH3-3 and BnSLAH3-4 genes were highly expressed in stamens, pistils, roots, leaves, silique and sepals, and they may be involved nitrate uptake and transport in B. napus. The characteristics and analysis of gene structure, physic-chemical properties, homologs phylogeny and experimental data provide a framework for further analysis of the BnSLAC/SLAH genes to define their biological functions during stress responses, as well as growth and development in B. napus.

Supplementary Materials

The following are available online at https://www.mdpi.com/article/10.3390/ijms22094671/s1, Table S1: SLAC/SLAH genes in the Brassica rapa genome, Table S2: SLAC/SLAH genes in the Brassica oleracea genome, Table S3: Primers used for gene cloning and qRT-PCR, Table S4: The FPKM values of expression levels of BnSLAC/SLAH genes in different tissues, Figure S1: Prediction of SLAC/SLAH protein interaction in Arabidopsis, Figure S2: Electropherograms of RNA samples from high nitrogen level (HN, 7.5 mM), low nitrogen level (LN, 0.19 mM) and ABA treatment (100 μM), Figure S3: Melt curves of primers used for qRT-PCR.

Author Contributions

Conceptualization, Y.N. and Y.G.; methodology, Y.X.; software, Y.N.; validation, Y.N., Y.X. and X.W.; formal analysis, A.A.; writing—original draft preparation, Y.N.; writing—review and editing, Y.Z., H.T. and Y.G. All authors have read and agreed to the published version of the manuscript.

Funding

This research was supported by grants from the National Key R&D Program of China (2018YFD0200907), grants from the Special Fund for Agro-scientific Research in the Public Interest (201503124), and the Innovative Research Team Plan of the Agriculture Ministry.

Institutional Review Board Statement

This study did not involve humans or animals.

Informed Consent Statement

This study did not involve humans.

Data Availability Statement

The RNA-Seq data used in this study are available in the Sequence Read Archive (SRA) at NCBI (SRA accession: PRJNA394926) repository.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Roelfsema, M.R.G.; Hedrich, R.; Geiger, D. Anion channels: Master switches of stress responses. Trends Plant Sci. 2012, 17, 221–229. [Google Scholar] [CrossRef] [PubMed]
  2. Vahisalu, T.; Kollist, H.; Wang, Y.F.; Nishimura, N.; Chan, W.Y.; Valerio, G.; Lamminmaki, A.; Brosche, M.; Moldau, H.; Desikan, R.; et al. SLAC1 is required for plant guard cell S-type anion channel function in stomatal signalling. Nature 2008, 452, 487-U415. [Google Scholar] [CrossRef]
  3. Brandt, B.; Brodsky, D.E.; Xue, S.W.; Negi, J.; Iba, K.; Kangasjarvi, J.; Ghassemian, M.; Stephan, A.B.; Hu, H.H.; Schroeder, J.I. Reconstitution of abscisic acid activation of SLAC1 anion channel by CPK6 and OST1 kinases and branched ABI1 PP2C phosphatase action. Proc. Natl. Acad. Sci. USA 2012, 109, 10593–10598. [Google Scholar] [CrossRef] [Green Version]
  4. Geiger, D.; Scherzer, S.; Mumm, P.; Marten, I.; Ache, P.; Matschi, S.; Liese, A.; Wellmann, C.; Al-Rasheid, K.A.S.; Grill, E.; et al. Guard cell anion channel SLAC1 is regulated by CDPK protein kinases with distinct Ca2+ affinities. Proc. Natl. Acad. Sci. USA 2010, 107, 8023–8028. [Google Scholar] [CrossRef] [Green Version]
  5. Geiger, D.; Maierhofer, T.; AL-Rasheid, K.A.S.; Scherzer, S.; Mumm, P.; Liese, A.; Ache, P.; Wellmann, C.; Marten, I.; Grill, E.; et al. Stomatal Closure by Fast Abscisic Acid Signaling Is Mediated by the Guard Cell Anion Channel SLAH3 and the Receptor RCAR1. Sci. Signal 2011, 4, ra32. [Google Scholar] [CrossRef] [PubMed]
  6. Vahisalu, T.; Puzorjova, I.; Brosche, M.; Valk, E.; Lepiku, M.; Moldau, H.; Pechter, P.; Wang, Y.S.; Lindgren, O.; Salojarvi, J.; et al. Ozone-triggered rapid stomatal response involves the production of reactive oxygen species, and is controlled by SLAC1 and OST1. Plant J. 2010, 62, 442–453. [Google Scholar] [CrossRef] [PubMed]
  7. Schafer, N.; Maierhofer, T.; Herrmann, J.; Jorgensen, M.E.; Lind, C.; von Meyer, K.; Lautner, S.; Fromm, J.; Felder, M.; Hetherington, A.M.; et al. A Tandem Amino Acid Residue Motif in Guard Cell SLAC1 Anion Channel of Grasses Allows for the Control of Stomatal Aperture by Nitrate. Curr. Biol. 2018, 28, 1370-U1145. [Google Scholar] [CrossRef] [Green Version]
  8. Gutermuth, T.; Lassig, R.; Portes, M.T.; Maierhofer, T.; Romeis, T.; Borst, J.W.; Hedrich, R.; Feijo, J.A.; Konrad, K.R. Pollen Tube Growth Regulation by Free Anions Depends on the Interaction between the Anion Channel SLAH3 and Calcium-Dependent Protein Kinases CPK2 and CPK20. Plant Cell 2013, 25, 4525–4543. [Google Scholar] [CrossRef] [Green Version]
  9. Zheng, X.J.; He, K.; Kleist, T.; Chen, F.; Luan, S. Anion channel SLAH3 functions in nitrate-dependent alleviation of ammonium toxicity in Arabidopsis. Plant Cell Environ. 2015, 38, 474–486. [Google Scholar] [CrossRef] [PubMed]
  10. Maierhofer, T.; Lind, C.; Huttl, S.; Scherzer, S.; Papenfuss, M.; Simon, J.; Al-Rasheid, K.A.S.; Ache, P.; Rennenberg, H.; Hedrich, R.; et al. A Single-Pore Residue Renders the Arabidopsis Root Anion Channel SLAH2 Highly Nitrate Selective. Plant Cell 2014, 26, 2554–2567. [Google Scholar] [CrossRef] [Green Version]
  11. Qiu, J.E.; Henderson, S.W.; Tester, M.; Roy, S.J.; Gilliham, M. SLAH1, a homologue of the slow type anion channel SLAC1, modulates shoot Cl- accumulation and salt tolerance in Arabidopsis thaliana. J. Exp. Bot. 2016, 67, 4495–4505. [Google Scholar] [CrossRef] [Green Version]
  12. Cubero-Font, P.; Maierhofer, T.; Jaslan, J.; Rosales, M.A.; Espartero, J.; Diaz-Rueda, P.; Muller, H.M.; Hurter, A.L.; AL-Rasheid, K.A.S.; Marten, I.; et al. Silent S-Type Anion Channel Subunit SLAH1 Gates SLAH3 Open for Chloride Root-to-Shoot Translocation. Curr. Biol. 2016, 26, 2213–2220. [Google Scholar] [CrossRef] [Green Version]
  13. Geiger, D.; Scherzer, S.; Mumm, P.; Stange, A.; Marten, I.; Bauer, H.; Ache, P.; Matschi, S.; Liese, A.; Al-Rasheid, K.A.S.; et al. Activity of guard cell anion channel SLAC1 is controlled by drought-stress signaling kinase-phosphatase pair. Proc. Natl. Acad. Sci. USA 2009, 106, 21425–21430. [Google Scholar] [CrossRef] [Green Version]
  14. Sun, D.; Fang, X.; Xiao, C.; Ma, Z.; Huang, X.; Su, J.; Li, J.; Wang, J.; Wang, S.; Luan, S.; et al. Kinase SnRK1.1 Regulates nitrate channel SLAH3 Engaged in Nitrate-Dependent Alleviation of Ammonium Toxicity. Plant Physiol. 2021. [Google Scholar] [CrossRef]
  15. Sun, S.J.; Qi, G.N.; Gao, Q.F.; Wang, H.Q.; Yao, F.Y.; Hussain, J.; Wang, Y.F. Protein kinase OsSAPK8 functions as an essential activator of S-type anion channel OsSLAC1, which is nitrate-selective in rice. Planta 2016, 243, 489–500. [Google Scholar] [CrossRef]
  16. Liu, X.H.; Mak, M.; Babla, M.; Wang, F.F.; Chen, G.; Veljanoski, F.; Wang, G.; Shabala, S.; Zhou, M.X.; Chen, Z.H. Linking stomatal traits and expression of slow anion channel genes HvSLAH1 and HvSLAC1 with grain yield for increasing salinity tolerance in barley. Front. Plant Sci. 2014, 5, 634. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  17. Qi, G.N.; Yao, F.Y.; Ren, H.M.; Sun, S.J.; Tan, Y.Q.; Zhang, Z.C.; Qiu, B.S.; Wang, Y.F. The S-Type Anion Channel ZmSLAC1 Plays Essential Roles in Stomatal Closure by Mediating Nitrate Efflux in Maize. Plant Cell Physiol. 2018, 59, 614–623. [Google Scholar] [CrossRef] [PubMed]
  18. Chen, G.D.; Wang, L.; Chen, Q.; Qi, K.J.; Yin, H.; Cao, P.; Tang, C.; Wu, X.; Zhang, S.L.; Wang, P.; et al. PbrSLAH3 is a nitrate-selective anion channel which is modulated by calcium-dependent protein kinase 32 in pear. BMC Plant Biol. 2019, 19, 1–12. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  19. Jaborsky, M.; Maierhofer, T.; Olbrich, A.; Escalante-Perez, M.; Muller, H.M.; Simon, J.; Krol, E.; Cuin, T.A.; Fromm, J.; Ache, P.; et al. SLAH3-type anion channel expressed in poplar secretory epithelia operates in calcium kinase CPK-autonomous manner. New Phytol. 2016, 210, 922–933. [Google Scholar] [CrossRef] [Green Version]
  20. Kurusu, T.; Saito, K.; Horikoshi, S.; Hanamata, S.; Negi, J.; Yagi, C.; Kitahata, N.; Iba, K.; Kuchitsu, K. An S-Type Anion Channel SLAC1 Is Involved in Cryptogein-Induced Ion Fluxes and Modulates Hypersensitive Responses in Tobacco BY-2 Cells. PLoS ONE 2013, 8, e70623. [Google Scholar] [CrossRef]
  21. Chen, G.D.; Li, X.L.; Qiao, X.; Li, J.M.; Wang, L.; Kou, X.B.; Wu, X.; Wang, G.M.; Yin, H.; Wang, P.; et al. Genome-wide survey and expression analysis of the SLAC/SLAH gene family in pear (Pyrus bretschneideri) and other members of the Rosaceae. Genomics 2019, 111, 1097–1107. [Google Scholar] [CrossRef]
  22. Bouchet, A.S.; Laperche, A.; Bissuel-Belaygue, C.; Baron, C.; Morice, J.; Rousseau-Gueutin, M.; Dheu, J.E.; George, P.; Pinochet, X.; Foubert, T.; et al. Genetic basis of nitrogen use efficiency and yield stability across environments in winter rapeseed. BMC Genet. 2016, 17. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  23. Chalhoub, B.; Denoeud, F.; Liu, S.Y.; Parkin, I.A.P.; Tang, H.B.; Wang, X.Y.; Chiquet, J.; Belcram, H.; Tong, C.B.; Samans, B.; et al. Early allopolyploid evolution in the post-Neolithic Brassica napus oilseed genome. Science 2014, 345, 950–953. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  24. Hedrich, R.; Geiger, D. Biology of SLAC1-type anion channels—From nutrient uptake to stomatal closure. New Phytol. 2017, 216, 46–61. [Google Scholar] [CrossRef] [Green Version]
  25. Negi, J.; Matsuda, O.; Nagasawa, T.; Oba, Y.; Takahashi, H.; Kawai-Yamada, M.; Uchimiya, H.; Hashimoto, M.; Iba, K. CO2 regulator SLAC1 and its homologues are essential for anion homeostasis in plant cells. Nature 2008, 452, 483-U413. [Google Scholar] [CrossRef] [PubMed]
  26. Wang, Y.Y.; Hsu, P.K.; Tsay, Y.F. Uptake, allocation and signaling of nitrate. Trends Plant Sci. 2012, 17, 458–467. [Google Scholar] [CrossRef]
  27. Wang, M.; Zhang, P.L.; Liu, Q.; Li, G.J.; Di, D.W.; Xia, G.M.; Kronzucker, H.J.; Fang, S.; Chu, J.F.; Shi, W.M. TaANR1-TaBG1 and TaWabi5-TaNRT2s/NARs Link ABA Metabolism and Nitrate Acquisition in Wheat Roots(1) ([OPEN]). Plant Physiol. 2020, 182, 1440–1453. [Google Scholar] [CrossRef] [Green Version]
  28. Han, Y.L.; Liao, J.Y.; Yu, Y.; Song, H.X.; Rong, N.; Guan, C.Y.; Lepo, J.E.; Ismail, A.M.; Zhang, Z.H. Exogenous abscisic acid promotes the nitrogen use efficiency of Brassica napus by increasing nitrogen remobilization in the leaves. J. Plant Nutr. 2017, 40, 2540–2549. [Google Scholar] [CrossRef]
  29. Lee, T.H.; Tang, H.B.; Wang, X.Y.; Paterson, A.H. PGDD: A database of gene and genome duplication in plants. Nucleic Acids Res. 2013, 41, D1152–D1158. [Google Scholar] [CrossRef] [PubMed]
  30. Krzywinski, M.; Schein, J.; Birol, I.; Connors, J.; Gascoyne, R.; Horsman, D.; Jones, S.J.; Marra, M.A. Circos: An information aesthetic for comparative genomics. Genome Res. 2009, 19, 1639–1645. [Google Scholar] [CrossRef] [Green Version]
  31. Drechsel, G.; Bergler, J.; Wippel, K.; Sauer, N.; Vogelmann, K.; Hoth, S. C-terminal armadillo repeats are essential and sufficient for association of the plant U-box armadillo E3 ubiquitin ligase SAUL1 with the plasma membrane. J. Exp. Bot. 2011, 62, 775–785. [Google Scholar] [CrossRef] [Green Version]
  32. Allender, C.J.; King, G.J. Origins of the amphiploid species Brassica napus L. investigated by chloroplast and nuclear molecular markers. BMC Plant Biol. 2010, 10, 1–9. [Google Scholar] [CrossRef] [Green Version]
  33. Aufhammer, W.; Kubler, E.; Bury, M. Nitrogen Uptake and Nitrogen Residuals of Winter Oilseed Rape and Fallout Rape. J. Agron. Crop Sci. 1994, 172, 255–264. [Google Scholar] [CrossRef]
  34. Rathke, G.W.; Christen, O.; Diepenbrock, W. Effects of nitrogen source and rate on productivity and quality of winter oilseed rape (Brassica napus L.) grown in different crop rotations. Field Crop. Res. 2005, 94, 103–113. [Google Scholar] [CrossRef]
  35. Swarbreck, D.; Wilks, C.; Lamesch, P.; Berardini, T.Z.; Garcia-Hernandez, M.; Foerster, H.; Li, D.; Meyer, T.; Muller, R.; Ploetz, L.; et al. The Arabidopsis Information Resource (TAIR): Gene structure and function annotation. Nucleic Acids Res. 2008, 36, D1009–D1014. [Google Scholar] [CrossRef] [PubMed]
  36. Tamura, K.; Stecher, G.; Peterson, D.; Filipski, A.; Kumar, S. MEGA6: Molecular Evolutionary Genetics Analysis Version 6.0. Mol. Biol. Evol. 2013, 30, 2725–2729. [Google Scholar] [CrossRef] [Green Version]
  37. Librado, P.; Rozas, J. DnaSP v5: A software for comprehensive analysis of DNA polymorphism data. Bioinformatics 2009, 25, 1451–1452. [Google Scholar] [CrossRef] [Green Version]
  38. Tang, H.B.; Wang, X.Y.; Bowers, J.E.; Ming, R.; Alam, M.; Paterson, A.H. Unraveling ancient hexaploidy through multiply-aligned angiosperm gene maps. Genome Res. 2008, 18, 1944–1954. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  39. Leblanc, A.; Renault, H.; Lecourt, J.; Etienne, P.; Deleu, C.; Le Deunff, E. Elongation changes of exploratory and root hair systems induced by aminocyclopropane carboxylic acid and aminoethoxyvinylglycine affect nitrate uptake and BnNrt2.1 and BnNrt1.1 transporter gene expression in oilseed rape. Plant Physiol. 2008, 146, 1928–1940. [Google Scholar] [CrossRef] [Green Version]
  40. Hu, B.; Jin, J.P.; Guo, A.Y.; Zhang, H.; Luo, J.C.; Gao, G. GSDS 2.0: An upgraded gene feature visualization server. Bioinformatics 2015, 31, 1296–1297. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  41. Zhao, H.M.; Ma, H.L.; Yu, L.; Wang, X.; Zhao, J. Genome-Wide Survey and Expression Analysis of Amino Acid Transporter Gene Family in Rice (Oryza sativa L.). PLoS ONE 2012, 7, e49210. [Google Scholar] [CrossRef] [PubMed]
  42. Bailey, T.L.; Boden, M.; Buske, F.A.; Frith, M.; Grant, C.E.; Clementi, L.; Ren, J.Y.; Li, W.W.; Noble, W.S. MEME SUITE: Tools for motif discovery and searching. Nucleic Acids Res. 2009, 37, W202–W208. [Google Scholar] [CrossRef]
  43. Maher, C.; Stein, L.; Ware, D. Evolution of Arabidopsis microRNA families through duplication events. Genome Res. 2006, 16, 510–519. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  44. Lescot, M.; Dehais, P.; Thijs, G.; Marchal, K.; Moreau, Y.; Van de Peer, Y.; Rouze, P.; Rombauts, S. PlantCARE, a database of plant cis-acting regulatory elements and a portal to tools for in silico analysis of promoter sequences. Nucleic Acids Res. 2002, 30, 325–327. [Google Scholar] [CrossRef]
  45. Qin, L.; Walk, T.C.; Han, P.P.; Chen, L.Y.; Zhang, S.; Li, Y.S.; Hu, X.J.; Xie, L.H.; Yang, Y.; Liu, J.P.; et al. Adaption of Roots to Nitrogen Deficiency Revealed by 3D Quantification and Proteomic Analysis. Plant Physiol. 2019, 179, 329–347. [Google Scholar] [CrossRef] [PubMed] [Green Version]
Ijms 22 04671 g001 550
Figure 1. Phylogenetic tree of SLAC/SLAHs among B. napus, B. rapa, B. oleracea, A. thaliana and O. sativa. The tree was generated using MEGA 6.0 with the Neighbor-Joining method. The proteins clustered into three subgroups. Yellow, blue and green sections indicate the three subfamilies of the SLAC/SLAH proteins. B. napus SLAC/SLAHs were marked with red dot.
Figure 1. Phylogenetic tree of SLAC/SLAHs among B. napus, B. rapa, B. oleracea, A. thaliana and O. sativa. The tree was generated using MEGA 6.0 with the Neighbor-Joining method. The proteins clustered into three subgroups. Yellow, blue and green sections indicate the three subfamilies of the SLAC/SLAH proteins. B. napus SLAC/SLAHs were marked with red dot.
Ijms 22 04671 g001
Ijms 22 04671 g002 550
Figure 2. Chromosomal localizations and gene duplication of BnSLAC/SLAH genes. Tandem duplication genes were circled with red line, and segmental duplication were linked with green lines.
Figure 2. Chromosomal localizations and gene duplication of BnSLAC/SLAH genes. Tandem duplication genes were circled with red line, and segmental duplication were linked with green lines.
Ijms 22 04671 g002
Ijms 22 04671 g003 550
Figure 3. The phylogenetic relationships, gene structures and conserved motifs of SLAC/SLAHs in B. napus.
Figure 3. The phylogenetic relationships, gene structures and conserved motifs of SLAC/SLAHs in B. napus.
Ijms 22 04671 g003
Ijms 22 04671 g004 550
Figure 4. Cis-elements in the 2.0-kb upstream regions of the BnSLAC/SLAH family genes in B. napus. On the left, five different color (brown, blue, yellow, green and purple) represent various types of cis-elements, including core promoter elements, light responsive, phytohormone responsive, abiotic stress responsive and others related to growth. On the right, Red square shows the presence of cis-element while the green indicates the absence of cis-element
Figure 4. Cis-elements in the 2.0-kb upstream regions of the BnSLAC/SLAH family genes in B. napus. On the left, five different color (brown, blue, yellow, green and purple) represent various types of cis-elements, including core promoter elements, light responsive, phytohormone responsive, abiotic stress responsive and others related to growth. On the right, Red square shows the presence of cis-element while the green indicates the absence of cis-element
Ijms 22 04671 g004
Ijms 22 04671 g005 550
Figure 5. Syntenic relationships of SLAC/SLAH genes among B. napus, B. rapa, B. oleracea and Arabidopsis. The circular forms of Brassica species and Arabidopsis chromosomes are shown in different colors. The approximate positions of the SLAC/SLAH genes are marked with short black lines on the circles. Gene pairs with syntenic relationships are joined by the colored lines.
Figure 5. Syntenic relationships of SLAC/SLAH genes among B. napus, B. rapa, B. oleracea and Arabidopsis. The circular forms of Brassica species and Arabidopsis chromosomes are shown in different colors. The approximate positions of the SLAC/SLAH genes are marked with short black lines on the circles. Gene pairs with syntenic relationships are joined by the colored lines.
Ijms 22 04671 g005
Ijms 22 04671 g006 550
Figure 6. Subcellular localization of BnSLAC/SLAHs in B. napus. The selected BnSLAC/SLAH-GFP fusion proteins (BnSLAH1-1-GFP, BnSLAH3-2-GFP and PbrSLAH3-3-GFP) as well as 35S-GFP as the control were independently transiently expressed in tobacco leaves and imaged under a confocal microscope. Bars = 50 μm.
Figure 6. Subcellular localization of BnSLAC/SLAHs in B. napus. The selected BnSLAC/SLAH-GFP fusion proteins (BnSLAH1-1-GFP, BnSLAH3-2-GFP and PbrSLAH3-3-GFP) as well as 35S-GFP as the control were independently transiently expressed in tobacco leaves and imaged under a confocal microscope. Bars = 50 μm.
Ijms 22 04671 g006
Ijms 22 04671 g007 550
Figure 7. Expression patterns of BnSLAC/SLAH family genes in twelve tissues. The expression data were gained from the RNA-seq data and calculated by fragments per kilobase of exon model per million (FPKM) values. The label above the heatmap represents the different tissues of B. napus ZS11, the right side of the heatmap represents different BnSLAC/SLAH genes. The colour bar represents log2 FPKM values.
Figure 7. Expression patterns of BnSLAC/SLAH family genes in twelve tissues. The expression data were gained from the RNA-seq data and calculated by fragments per kilobase of exon model per million (FPKM) values. The label above the heatmap represents the different tissues of B. napus ZS11, the right side of the heatmap represents different BnSLAC/SLAH genes. The colour bar represents log2 FPKM values.
Ijms 22 04671 g007
Ijms 22 04671 g008 550
Figure 8. The expression patterns of selected BnSLAC/SLAH genes in leaves (A) and roots (B) after 3 h, 12 h, 24 h, 3 d and 7 d HN (7.5 mM) and LN (0.19 mM) treatment. The asterisk represents statistical significance (Student’s t-test, “*” p < 0.05, “**” p < 0.01) in comparison with control. The data are shown as mean values ± SE (n = 3).
Figure 8. The expression patterns of selected BnSLAC/SLAH genes in leaves (A) and roots (B) after 3 h, 12 h, 24 h, 3 d and 7 d HN (7.5 mM) and LN (0.19 mM) treatment. The asterisk represents statistical significance (Student’s t-test, “*” p < 0.05, “**” p < 0.01) in comparison with control. The data are shown as mean values ± SE (n = 3).
Ijms 22 04671 g008
Ijms 22 04671 g009 550
Figure 9. The expression patterns of selected BnSLAC/SLAH genes in leaves (A) and roots (B) after 3 h, 12 h and 24 h ABA (100 uM) treatment. The asterisk represents statistical significance (Student’s t-test, “*” p < 0.05, “**” p < 0.01) in comparison with control. The data are shown as mean values ± SE (n = 3).
Figure 9. The expression patterns of selected BnSLAC/SLAH genes in leaves (A) and roots (B) after 3 h, 12 h and 24 h ABA (100 uM) treatment. The asterisk represents statistical significance (Student’s t-test, “*” p < 0.05, “**” p < 0.01) in comparison with control. The data are shown as mean values ± SE (n = 3).
Ijms 22 04671 g009
Table
Table 1. Characterization of 23 SLAC/SLAH members identified in B. napus.
Table 1. Characterization of 23 SLAC/SLAH members identified in B. napus.
Gene FeaturesProtein Features
Putative Gene Name Gene IDExon No.Chr. LocationProtein Length (aa)MW (KDa)pIGRAVYTMD
BnSLAC1-1LOC1064455773A06 (5554449..5557045)58866.649.22−0.0468
BnSLAC1-2LOC1112062993C05(7993671..7996097)58666.189.30.038
BnSLAH1-1LOC1063652682A09(7429214..7430811)38743.409.280.58
BnSLAH1-2LOC1063736082C04(33942406..33943955)38743.819.290.4558
BnSLAH1-3LOC1063730082A09 (9740564..9742068)38142.999.270.5668
BnSLAH1-4LOC1064000982C09 (12333910..12335636)38143.049.110.5518
BnSLAH2-1LOC1064098454C07 (54435207..54437674)49856.269.30.23110
BnSLAH2-2LOC1063598354A07(11635779..11638473)51557.789.650.23410
BnSLAH2-3LOC1064138964C03 (16311421..16314068)51557.809.541.23510
BnSLAH2-4LOC1064527143A01 (12909477..12912194)35940.319.440.6269
BnSLAH2-5LOC1063761453C01 (21092721..21095213)35940.369.440.6139
BnSLAH3-1LOC1063482905A06 (25034448..25039244)63672.428.34−0.0929
BnSLAH3-2LOC1064364655C07 (46161649..46166453)63672.418.39−0.11110
BnSLAH3-3LOC1064284595C09 (3799924..3805933)62371.086.88−0.1059
BnSLAH3-4LOC1063660195A01(11440762..11445373)62371.077.07−0.119
BnSLAH3-5LOC1064152535A02 (32761829..32765883)64072.598.95−0.079
BnSLAH3-6LOC1063789925C02 (66036653..66039643)63672.248.88−0.0749
BnSLAH4-1LOC1064179793A09 (9746804..9754357)47654.279.260.30810
BnSLAH4-2LOC1063977063C09 (12312761..12318196)47053.589.220.3210
BnSLAH4-3LOC1064122013C08(29119049..29124524)47053.609.210.32410
BnSLAH4-4LOC1063953192C04 (33847294..33848994)32937.578.970.6496
BnSLAH4-5LOC1063788466C01(36168496..36175425)53360.566.040.198
BnSLAH4-6LOC1063788474C01(36145965..36151210)44751.016.430.2778
Table
Table 2. The synonymous substitution rates (Ks) and non-synonymous substitution rates (Ka) of the SLAC/SLAHs in A. thaliana and B. napus.
Table 2. The synonymous substitution rates (Ks) and non-synonymous substitution rates (Ka) of the SLAC/SLAHs in A. thaliana and B. napus.
Gene Name in A. thalianaGene Name in B. napusKaKsKa/Ks
AtSLAC1BnSLAC1-10.06160.45410.1356529
AtSLAC1BnSLAC1-20.06340.42920.1477167
AtSLAH1BnSLAH1-10.08670.46710.1856134
AtSLAH1BnSLAH1-20.09580.47020.2037431
AtSLAH1BnSLAH1-30.09290.43150.2152955
AtSLAH1BnSLAH1-40.09550.41820.2283596
AtSLAH2BnSLAH2-10.05670.37620.1507177
AtSLAH2BnSLAH2-20.06860.40210.1706043
AtSLAH2BnSLAH2-30.07060.41370.1706551
AtSLAH2BnSLAH2-40.07620.44820.1700134
AtSLAH2BnSLAH2-50.07410.460.161087
AtSLAH3BnSLAH3-10.06940.43530.1594303
AtSLAH3BnSLAH3-20.06860.38960.176078
AtSLAH3BnSLAH3-30.06380.44190.1443766
AtSLAH3BnSLAH3-40.06350.46990.1351351
AtSLAH3BnSLAH3-50.10840.47670.2273967
AtSLAH3BnSLAH3-60.10920.4910.2224033
AtSLAH4BnSLAH4-10.1720.06772.5406204
AtSLAH4BnSLAH4-20.17570.06792.5876289
AtSLAH4BnSLAH4-30.17560.0682.5823529
AtSLAH4BnSLAH4-40.15950.09891.6127401
AtSLAH4BnSLAH4-50.21830.15751.3860317
AtSLAH4BnSLAH4-60.21830.15751.3860317
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Nan, Y.; Xie, Y.; Atif, A.; Wang, X.; Zhang, Y.; Tian, H.; Gao, Y. Identification and Expression Analysis of SLAC/SLAH Gene Family in Brassica napus L. Int. J. Mol. Sci. 2021, 22, 4671. https://doi.org/10.3390/ijms22094671

AMA Style

Nan Y, Xie Y, Atif A, Wang X, Zhang Y, Tian H, Gao Y. Identification and Expression Analysis of SLAC/SLAH Gene Family in Brassica napus L. International Journal of Molecular Sciences. 2021; 22(9):4671. https://doi.org/10.3390/ijms22094671

Chicago/Turabian Style

Nan, Yunyou, Yuyu Xie, Ayub Atif, Xiaojun Wang, Yanfeng Zhang, Hui Tian, and Yajun Gao. 2021. "Identification and Expression Analysis of SLAC/SLAH Gene Family in Brassica napus L." International Journal of Molecular Sciences 22, no. 9: 4671. https://doi.org/10.3390/ijms22094671

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