-
A total of 24 SWEET homologs were identified based on the genome information of M. alba and named according to their orthologs in A. thaliana, P. trichocarpa or V. vinifera[45]. These MaSWEETs encode proteins ranging from 197 aa to 304 aa with molecular weight from 21.45 to 34.12 kDa and theoretical isoelectric points from 7.16 to 9.61 (Table 1). Subcellular localization prediction of these MaSWEETs showed that most of them (18/24) distributed on membrane structures such as plasma membrane (PM), tonoplast membrane (TM) and chloroplast thylakoid membrane (CTM). Phylogenetic analysis of MaSWEETs and SWEETs from model plants such as A. thaliana, P. trichocarpa and O. sativa showed that four clades were formed by these SWEETs (Fig. 1). According to previous studies, SWEETs in plants were generally classed as four phylogenetic clades which is in agreement with our results[8]. Major MaSWEETs (10/24) together with AtSWEET1, 2 and 3 belong to clade I. Clade II and IV contain five MaSWEETs each and Clade III contains four MaSWEETs (Fig. 1 and Table 1).
Table 1. SWEET gene family in Morus alba.
Clade Gene name Accession no. Gene ID CDS Size
(bp)Protein physicochemical characteristics TMHs Subcellular MtN3/Saliva (PQ-Loop Repeat) Length (aa) MW (kDa) pI Aliphatic index Localization* Domain Position I MaSWEET1a XM_024170697.1-0 M.alba_G0012049 729 242 26.55 9.61 113.51 7 CTM 6-94, 131-209 I MaSWEET1b XM_024170698.1-0 M.alba_G0012049 678 225 21.45 9.51 115.74 6 CTM 1-49, 86-164 I MaSWEET2a XM_024163709.1-0 M.alba_G0019244 708 235 25.94 8.39 129.79 7 CTM 18-104, 137-221 I MaSWEET2b XM_024164777.1-0 M.alba_G0010863 705 234 25.92 8.98 122.86 7 TM 17-101, 137-218 I MaSWEET2c XM_024163707.1-0 M.alba_G0019244 774 257 28.58 8.58 132.33 7 PM 54-126, 159-243 I MaSWEET2d XM_024163712.1-0 M.alba_G0019244 681 226 25.27 8.22 128.89 6 EX 23-95, 128-212 I MaSWEET2e XM_024163708.1-0 M.alba_G0019244 729 242 26.72 8.49 128.06 7 CTM 39-111, 144-228 I MaSWEET2f XM_024163703.1-0 M.alba_G0019244 777 258 28.8 8.49 132.56 8 PM 41-127, 160-244 I MaSWEET2g XM_024163711.1-0 M.alba_G0019244 684 227 25.49 7.62 129.16 7 PM 10-96, 129-213 I MaSWEET3 XM_010099554.2-0 M.alba_G0003063 783 260 29.01 8.89 115.73 7 PM 9-98, 132-216 II MaSWEET4a XM_010091939.1-0 M.alba_G0009276 735 244 27.45 9.28 109.39 7 CTM 10-98, 134-218 II MaSWEET4b XM_010113461.2-0 M.alba_G0001536 738 245 27.4 8.98 122.08 7 EX 11-97, 134-216 II MaSWEET5 XM_024168739.1-0 M.alba_G0018295 711 236 26.64 7.63 120.93 7 CY 10-93, 131-127 II MaSWEET7a XM_010108966.2-0 M.alba_G0005110 774 257 28.32 9.57 128.56 7 CTM 11-95, 134-218 II MaSWEET7b XM_010108964.1-0 M.alba_G0005109 789 262 28.99 9 124.96 7 CTM 10-97, 134-218 III MaSWEET10 XM_010095631.2-0 M.alba_G0018016 888 295 33.12 8.86 120.31 7 CTM 11-96, 132-216 III MaSWEET11a XM_010114440.2-0 M.alba_G0016901 804 267 29.63 9.47 124.08 7 CTM 11-99, 135-220 III MaSWEET11b XM_010095633.2-0 M.alba_G0018015 915 304 34.12 7.57 112.89 7 CTM 12-99, 134-219 III MaSWEET15 XM_010092381.2-0 M.alba_G0006767 885 294 33.42 7.16 109.01 7 CTM 12-99, 133-219 IV MaSWEET16 XM_024167733.1-0 M.alba_G0014617 909 302 33.2 9.08 114.27 7 CTM 20-92, 129-212 IV MaSWEET17a XM_024171451.1-0 M.alba_G0014614 708 235 26.46 8.71 119.87 5 CY 5-78, 116-198 IV MaSWEET17b XM_024167902.1-0 M.alba_G0014613 753 250 28.05 8.71 120.88 6 PM 8-93, 131-213 IV MaSWEET17c XM_024171286.1-0 M.alba_G0008195 720 239 26.72 8.94 111.72 7 CY 6-92, 129-212 IV MaSWEET17d XM_024171287.1-1 M.alba_G0008196 723 240 27 9.43 122.67 7 CTM 6-92, 127-213 * The subcellular localizations were predicted by WoLFPSORT. PM, plasma membrane; EX, extracellular; CY, cytoplasmic; TM, tonoplast membrane; CTM, chloroplast thylakoid membrane. Figure 1.
Phylogenetic relationships of the SWEET family genes in Arabidopsis, Oryza sativa, Populus, Vitis vinifera, and Morus alba. The sequences of the 104 SWEET proteins from the above four plant species were aligned by Clustal Omega, and the phylogenetic tree was constructed by the MEGA 11.0 using the NJ method with 1000 bootstrap replicates. The proteins from Arabidopsis, Oryza sativa, Populus, Vitis vinifera, and Morus alba are indicated with the prefixes of At, Os, Pt, Vv, and Ma, respectively.
Chromosomal location and gene duplication
-
MaSWEETs distributed on 12 chromosomes except chromosome 1 and 3. Chromosome 5 occupied five MaSEETs which formed a gene cluster. Chromosome 5 is the chromosome that had the most SWEETs and the following is chromosome 6 and 8 which had four MaWSEETs each. There was only one MaSWEET locating on chromosome 4 (Fig. 2). In addition, there were three MaSWEETs on chromosome 2 and 7 respectively and two MaSWEETs on chromosome 1 and 5 respectively. Gene duplication including block duplication and tandem duplication is the main cause for gene family expansion. Tandem duplications were found on chromosome 2, 6 and 12 (linked by red lines in Fig. 2). It is also interesting to find several possible gene clusters such as MaSWEET1a-b on chromosome 5, MaSWEET2a-g on chromosome 9, MaSWEET7 a-b on chromosome 13, MaSWEET17a-b on chromosome 12 and MaSWEET17c-d on chromosome 6. Two gene pairs (MaSWEET4b/MaSWEET5 and MaSWEET11b/ MaSWEET15) resulting from block duplications were also marked (linked by black lines in Fig. 2).
Figure 2.
Distribution of MaSWEET genes in Morus alba chromosomes. The tandem gene pairs are linked by red lines. The block duplications gene pairs are marked by black lines. The scale is provided in megabase (Mb).
Sequence analysis of MaSWEETs
-
MaSWEETs always located on membrane structures with transmembrane domains. The prediction results of MaSWEETs using DeepTMHMM showed that most MaSWEETs posed seven types of transmembrane helices (TMH) named TMH1-7 (Table 1, Supplemental Fig. S1). Alignment and conserved motif analysis showed that almost all MaSWEETs kept the conserved TMH and active sites Tyr and Asp indicating by red full triangles (Fig. 3). The active residues Tyr and Asp were reported to be involved in forming hydrogen bonds to ensure sugar transport activity[14]. In addition, all MaSWEETs except MaSWEET4a had a conserved Ser in each triple helix bundle (THB) which can be phosphorylated and is important for SWEET activity (Fig. 3). MaSWEET4a replaced Ser with Thr at the first Ser phosphorylation site between TMH1 and TMH2 which may also retain similar activity as Thr is also a common phosphorylation site. All MaSWEETs retained the conserved second Ser phosphorylation site between TMH5 and TMH6 (Fig. 3).
Figure 3.
Multiple sequence alignment of MaSWEET proteins. The positions of the TMHs are underlined. The positions of the active sites of tyrosine (Y) and aspartic acid (D) are indicated by red triangles. The conserved serine (S) phosphorylation sites are indicated by blue triangles.
Gene organization and promoter analysis of MaSWEETs
-
MaSWEETs gene structures were identified based on the annotation information of the M. alba genome. In summary, there were six MaSWEETs with six introns, 12 MaSWEETs with five introns and five MaSWEETs with four introns (Fig. 4). In addition, genes clustered together based on phylogenetic analysis are likely to show similar gene structures and length. For example, MaSWEET2a, c, d, e, f and g with six introns, MaSWEET7a and MaSWEET7b with four introns, and MaSWEET17a and MaSWEET17b with four introns.
Figure 4.
Gene organization of MaSWEETs and cis-elements in promoter regions of MaSWEETs. (a) Phylogenetic tree using 24 MaSWEETs. (b) Exon/intron structures of Morus alba L. SWEETs. (c) Cis-element distribution in the promoter regions of MaSWEETs.
Promoter region analysis of MaSWEETs indicated the possible function of MaSWEETs in response to hormones and environment stimulus. Among all the 22 types of cis-elements identified, most of them are light response elements accounting for 44% of the total elements (Supplemental Table S2). In addition, hormone response elements were also widely identified in the promoters of MaSWEETs (Fig. 4c). Most MaSWEETs had abscisic acid (ABA), salicylic acid (SA) or methyl jasmonate (MeJA) related response elements in their promoter regions. Especially, MaSWEET1a-b, MaSWEET16, MaSWEET17a-d had cis-elements involved in response to five types of hormones (ABA, SA, MeJA, auxin and gibberellins). Several Myb binding cis-elements were also identified in promoter regions of MaSWEET2a and MaSWEET10 (Fig. 4c, Supplemental Table S2).
Expression profile of MaSWEETs in different tissues of mulberry
-
The tissue or organ expression profiles of MaSWEETs were revealed. The MaSWEETs with high sequence identity (> 91%) are hard to distinguished by RT-qPCR and were determined by common primers to reveal their total transcription levels. MaSWEETs in phylogenetic clade I showed quite similar expression patterns with highest expression levels in leaf and relatively higher expression levels in early stages (S0 and S1) of fruit development except MaSWEET1a/b (Fig. 5a−d). MaSWEET1a/b showed higher expression levels in fruits during whole fruit development with highest expression level in fruit at S0 stage (Fig. 5a). However, MaSWEETs (MaSWEET10, 11a-b, and 15) in phylogenetic clade III showed preferential expression in fruits especially at the late stages (S2 and S3) (Fig. 5j−m). MaSWEETs in phylogenetic clade II had different expression patterns in different tissues or organs. MaSWEET4a and b showed highest expression level in fruit at the S1 stage while MaSWEET5 showed obviously preferential expression in xylem (Fig. 5e−g). MaSWEET7a and b showed similar expression pattern with MaSWEET2 cluster and MaSWEET3 from phylogenetic clade I. MaSWEETs in phylogenetic clade IV showed similar expression pattern with highest expression levels in leaf (Fig. 5b, d, h, i). MaSWEET16 also had higher expression in fruit with similar expression level at four different development stages (Fig. 5n).
Figure 5.
Transcript levels of MaSWEETs in leaves, xylem, phloem, and different development stages of fruit. Three technical replicates were analyzed. Error bars represent SE. Different letters indicate statistically significant differences (Duncan's test, p < 0.05).
Transcription-level responses of MaSWEETs to various stresses
-
Mulberry sclerotiniose is a fungal disease resulting from fungal pathogen infection. Most (20/24) MaSWEETs showed positive or negative responses to the fungal infection. MaSWEET1a/b, MaSWEET2 cluster, MaSWEET4b, and MaSWEET17 a-d showed a significant decrease of expression levels in diseased fruits with sclerotiniose compared with the expression levels in healthy fruits (Fig. 6). In contrast, MaSWEET2b, MaSWEET3, MaSWEET7b, MaSWEET10 and MaSWEET11a-b showed significant increases of expression levels in diseased fruits. MaSWEETs also played roles in response to various abiotic stresses including drought, water logging, cold and high temperature. MaSWEET1a/b showed a positive response to drought with significant increasing expression levels while other clade I MaSWEETs, MaSWEET2b and 3 significantly decreased their expression level under detected abiotic stresses (Fig. 7a−d). In contrast, MaSWEET16 significantly increased its expression level under detected abiotic stresses. MaSWEET4a-b, MaSWEET5 in phylogenetic clade II and MaSWEET11a-b in clade III showed similar response patterns with a significant increase of expression levels in response to low temperature (4 °C), high temperature (40 °C) or drought (Fig. 7e−g). MaSWEET15 showed high sensitivity for drought and significant increase of expression levels under drought stress. MaSWEET17a-d showed a negative response to temperature change with a significant decrease of expression levels under low temperature and high temperature treatments (Fig. 7o−q).
Figure 6.
Expression levels of 24 selected MaSWEET genes in response to fungi stress conditions. Three technical replicates were analyzed. Error bars represent SE. Asterisks indicate significant difference as determined by Student’s t-test (* p < 0.05; ** p < 0.01; *** p < 0.001; **** p < 0.0001).
Figure 7.
Expression levels of 24 selected MaSWEET genes in response to 4 °C, 40 °C , drought, and flood stress conditions. Three technical replicates were analyzed. Error bars represent SE. Asterisks indicate significant difference as determined by Student’s t-test (* p < 0.05 ; ** p < 0.01 ; *** p < 0.001 ; **** p < 0.0001 ).
Functional characterization of MaSWEET1a in response to sclerotiniose infection
-
MaSWEET1a/b is quite different from other MaSWEETs in clade I based on expression profile analysis, which showed preferential expression in fruits and a negative response to sclerotiniose pathogen (Ciboria shiraiana) infection. Our unpublished data indicated MaSWEET1a as key genes involved in the pathogen infection process based on comparative transcriptome analysis. Transient overexpression of MaSWEET1a in tobacco and VIGS knock-down of MaSWEET1a/b in mulberry were performed. RT-qPCR results validated the successful overexpression of MaSWEET1a in tobacco and knock-down of MaSWEET1a/b in mulberry (Fig. 8b, d). The expression level of MaSWEET1a affected the resistance to C. shiraiana infection in both tobacco and mulberry (Fig. 8a, d). Overexpression of MaSWEET1a decreased the resistance to C. shiraiana infection with more severe cell death symptoms observed in OE-line tobacco (Fig. 8a). Knock-down of MaSWEET1a/b in mulberry could increase the resistance to C. shiraiana infection in mulberry (Fig. 8d). These results proved that MaSWEET1a is an important negative regulator of resistance to C. shiraiana infection in mulberry.
Figure 8.
Functional characterization of MaSWEET1a/b in tobacco and mulberry. (a) Damage of tobacco overexpressing MaSWEET1a after infection by C. shiraiana. (b) Expression levels of MaSWEET1a in tobacco. MaSWEET1a-T1, T2, T3 are three independent treatments with transient overexpressed MaSWEET1a. (c) Damage of mulberry with knock-down MaSWEET1a/b after infection by C. shiraiana. (d) The expression level of MaSWEET1a/b in mulberry. MaSWEET1a/b-T1, T3, T4 are three independent treatments with down-regulation of MaSWEET1a/b using VIGS.
-
Functional studies on SWEETs have revealed that SWEET homologs not only act as loading and unloading transporters of sugars but also play critical roles in various biological processes. Typically, angiosperm genomes contain about 15–25 SWEET genes[4]. In the present study, a total of 24 SWEET genes were identified and clustered into four clades corresponding with the knowledge of SWEETs in angiosperm. Several SWEETs including MaSWEET1a/b, MaSWEET4a-b, MaSWEET10, MaSWEET11a-b and MaSWEET15 showed preferential expression in fruits indicating their possible roles in fruit development. Similar expression preference of AtSWEETs was also reported in Arabidopsis[2]. It is noted that fruit-preferential expressed MaSWEETs still showed temporal expression difference during fruit development indicating time-course regulation of MaSWEETs for fruit ripening in mulberry. Early-stage expressed MaSWEET1a/b and MaSWEET4a/b further differed from late-stage expressed MaSWEET10, MaSWEET11a/b and MaSWEET15 in terms of detailed expression patterns during fruit ripening. MaSWEET2b and 2 cluster genes (MaSWEET2a, c-g) showed preferential expression in leaves which is similar with the ortholog ZjSWEET2.2 in Ziziphus jujuba. ZjSWEET2.2 was reported to be involved in mediating sugar loading in leaves[46]. MaSWEET3, 7a-b,16 and 17a-d also showed highest expression levels in leaves indicating their possible roles in sugar source loading or unloading. MaSWEET11b showed higher expression levels in phloem and its ortholog AtSWEET11 in Arabidopsis was reported to be involved in sugar phloem loading[3].
Sugar signal is critical for plants in response to various stresses. Previous studies have shown that SWEETs participated in abiotic and biotic responses in many plant species including arabidopsis and rice[21, 25]. AtSWEET11, 12, 15 and 16 were reported to be involved in affecting cold tolerance in Arabidopsis[19−21]. HfSWEET17 was also reported as a positive regulator of resistance to cold stress in daylily[11]. Cold environment (4 °C) induced expression of MaSWEET4a, 4b, 5, 11a, 11b and 16 in mulberry. Interestingly, these cold-induced MaSWEETs can also be induced by high temperature. MaSWEET15 which is the ortholog of AtSWEET15 can be induced by drought as well as low or high temperature. MaSWEET15, MaSWEET1a/b. 4a, 4b, 5, 7a, 11a, 11b, and 16 also showed positive responses to drought. It is obvious that some SWEET genes can be induced by different stresses. AtSWEET15 was also reported to be induced by osmotic, drought and salinity[21]. MaSWEET4a, 4b, 11a, 11b and 16 can be induced by low or high temperature and drought indicating their important roles in response to various abiotic stresses in mulberry.
SWEETs were generally thought to 'support the enemy' during infection. SWEETs especially those that function as exporters generally facilitate the export of sugars out of host cells, which support pathogen growth in the apoplasm[15,47,48]. Clade III SWEETs including AtSWEET11, 12, OsSWEET11 were characterized as negative regulators of resistance to fungal infection and Clade III SWEETs including OsSWEET11, 13, 14 and GhSWEET10 were characterized as negative regulators of resistance to bacterial pathogen infection[2,24,49,50]. In contrast, clade I AtSWEET2, a glucose importer and clade III IbSWEET10 were reported as positive regulators of resistance to fungal infection[48,51]. Therefore, roles of SWEETs in response to pathogen infection may be quite different. Most MaSWEETs were disturbed in diseased fruits that resulted from sclerotiniose pathogen infection. MaSWEET2b, MaSWEET3 in clade I, MaSWEET7b in clade II, MaSWEET10 and MaSWEET11a-b in clade III showed significant increase in expression levels in diseased fruits while MaSWEET1a/b, MaSWEET2 cluster (MaSWEET2a, c-g) in clade I showed significant decrease of expression levels in diseased fruits. MaSWEET1a was further validated as a negative regulator of resistance to C. shiraiana infection. Given the fact that MaSWEET1a/b was repressed in diseased fruits, it is likely that a possible pathway through repression of MaSWEET1a exists in mulberry to defense pathogen infection.
In conclusion, we have performed a genome-wide investigation of SWEET genes in Morus and a comprehensive analysis including phylogenetic analysis, promoter analysis and expression profile analysis was also carried out. Their possible roles in development and response to abiotic and biotic stresses were addressed. In particular, the functional role of MaSWEET1a in regulation of tolerance to C. shiraiana infection was validated using both VIGS knock-down and transient overexpression in tobacco combined with inoculation of C. shiraiana. The results in this study provides a foundation for studying the function of the SWEET family in mulberry plants and provides a negative regulator of resistance to C. shiraiana infection for further genetic modification.
-
About this article
Cite this article
Kang X, Huang S, Feng Y, Fu R, Tang F, et al. 2023. SWEET transporters and their potential roles in response to abiotic and biotic stresses in mulberry. Beverage Plant Research 3:6 doi: 10.48130/BPR-2023-0006
SWEET transporters and their potential roles in response to abiotic and biotic stresses in mulberry
- Received: 24 November 2022
- Revised: 07 January 2023
- Accepted: 17 January 2023
- Published online: 02 March 2023
Abstract: Mulberry (Morus spp., Moraceae) is a traditional economic crop plant and is also being gradually utilized as a beverage plant. SWEETs (Sugars Will Eventually be Exported Transporter) are important sugar transporters involved in various biological processes and responses to various stresses. However, SWEETs in mulberry are still poorly studied without a comprehensive functional analysis of SWEETs. In the present study, a total of 24 SWEETs were identified using the Morus alba (Ma) genome. Phylogenetic analysis showed that these 24 MaSWEETs were clustered with SWEETs from Arabidopsis, Populus and Oryza and fell into four clades. MaSWEETs in the same clade are likely to pose similar intron/exon patterns. These MaSWEETs distributed on 12 chromosomes and tandem duplication and block duplication were responsible for the expansion of SWEETs in mulberry. Transmembrane domains and conserved active sites of Tyr and Asp were observed in MaSWEETs. Cis-elements in promoter regions of MaSWEETs indicated the possible function of MaSWEETs in response to hormones and environment stimulus. MaSWEETs showed quite different expression preference in tissues and organs indicating the possible function divergence. In addition, most MaSWEETs showed a disturbed expression levels in response to various abiotic stresses and Ciboria shiraiana infection. MaSWEET1a was functionally characterized as a negative regulator of resistance to C. shiraiana infection based on in vivo transient overexpression of MaSWEET1a in tobacco and down-regulation of MaSWEET1a/b in mulberry. Our results provided foundation for further functional dissection of SWEETs in mulberry and a potential regulator for genetic modification.
-
Key words:
- Mulberry /
- SWEET /
- Abiotic stress /
- Fungal pathogen /
- Functional divergence