Structural module of the two rice SIP-like proteins
Previous studies have demonstrated that the OsSIPs in rice(Oryza sativa) consists of two members OsSIP1 and OsSIP2 [17]. By constructing evolutionary trees for multiple species and analyzing amino acid sequence homology and identity, it was revealed that they belong to two different subfamilies, SIP1s and SIP2s, respectively (Fig. 1). Topological analysis of OsSIP1 and OsSIP2 revealed that both have a short N-terminus, 6 transmembrane structural domains and five loops that connect the TMs. The C-terminus, N-terminus, LB, and LD loops are intracellularly held, while the LA, LC, and LE loops are anchored extracellularly, Additionally, two NPA motifs, and an ar/R-restricted region are present (Fig. 2A, B).
Unlike typical AQPs, the NPA regions conserved between TM2 and TM3 of OsSIP1 and OsSIP2 changed to NPT and NPL, correspondingly. The ar/R-restricted regions, which govern their permeability to different substrates, are shown in (Fig. 2C, D). The water permeability of Arabidopsis AtSIPs and grape VvSIP1;1 has been confirmed by previous studies [24, 28]. Further research has found that the Homo sapiens hAQP11, which exhibits a closer homology to the plant SIPs family (Fig. 1), can transport water, H2O2 and glycerol [40, 41]. A comparison of the two conserved constrictions of OsSIPs with those of other SIP orthologs and other typical AQPs members that have been studied so far suggested that the OsSIP2 is likely more restrict in the pore layouts owing to the projected side strains of the Leu59 and His179 in the ar/R region (Fig. 2D). Based on the predicted structural modules (Fig. 2C, D), we suppose that OsSIPs may behavior more like a water channel rather than a glycerol channel. Finaly, putative phosphorylation sites are presented in the cytoplasmic LB loop and N/C-terminus OsSIPs (Fig. 2A, B). As phosphorylation at several sites in the LB loop and C-terminus of AQPs has been found to control their water transport activity in plants [42–45], these predicted sites are potentially regulatory in nature.
Phenotypical analysis of yeast cells expressing OsSIP1 and OsSIP2
A S. cerevisiae fps1 mutant, BY4742 (Δfps1), which lacks the endogenous aquaglyceroporin Fps1, was successfully used for functional complementation analysis to investigate the permeability of different aquaporins to several substrates [36, 46–50]. To investigate the impact of OsSIP1 and OsSIP2 on the transport of substrate, we assessed the growth of yeast BY4742 (Δfps1) cells transformed with pRS426met25-OsSIPs-GFP (C-terminal GFP fusion) or pRS426met25-GFP (negative control) on plates with or without a specified substrate. Furthermore, we employed pRS426met25-GFP transformed wild-type yeast strain BY4742 as a positive control. Fluorescence of the OsSIPs-GFP fusion proteins or GFP indicated successful expression of the proteins (Fig. 3A), which was further verified by western blot assays (representative as in Fig. 3B). The expression of OsSIP1-GFP and OsSIP2-GFP exhibited no toxicity to the yeast cells, as evidenced by comparable growth on plates without substrate addition that was identical to that of the cells transformed with the control vectors (Fig. 3).
When challenged by hyperosmotic stress imposed by high concentration of either KCl or sorbitol, the yeast cell accumulates glycerol intracellularly and closes the aquaglycerol channel Fps1 to survival, and the fps1 mutant is more resistant due to a lack of expressed water and/or glycerol channel [49]. Expression of both OsSIPs made the yeast cells more sensitive than the wild-type yeast and the fps1 mutant (Fig. 3C, D), indicating that a facilitated efflux of water could attribute to the defect growth. Thus, according to their phenotype, OsSIP1 displays higher water or glycerol permeability than OsSIP2.
To test whether the OsSIP1 and OsSIP2 could transport H2O2, growth assays were carried out in agar medium (SD-Ura, pH 5.5) with increasing concentrations of H2O2. The wild-type cell grew slightly better than the fps1 mutant under 2 mM H2O2 but could not withstand at 2.5 mM H2O2 (Fig. 3E), consistent with previously reported notion that oxidative stress rapidly induces Fps1 degradation at low H2O2 concentration [50]. Expression of OsSIP1 and OsSIP2 in the mutant resulted in a significant reduction of growth under increasing H2O2 concentration (Fig. 3E). Like the situation under hyperosmotic stress, OsSIP1 seemed to cause the yeast cell to be slightly more sensitive to H2O2 than did the OsSIP2 (Fig. 3E). The sensitivity to H2O2 is presumably contributed by the SIP proteins sit at the plasma membrane although their GFP fluorescent signal is predominant at the ER (Fig. 3E).
To determine ammonia permeation, the agar medium (SD-Ura, pH 5.5) was supplemented with the toxic ammonia analog MeA, which was taken up by the MEP proteins and accumulated in the cytosol causing detrimental effects to the cells. Detoxification could be achieved by the presence of a functional channel that permits the release of MeA outside the cells [36]. As expected, the fps1 mutant was unable to grow in MeA containing medium, while cells with FPS1 permitted sustainable growth in medium with up to 100 mM MeA (Fig. 3F). The expression of OsSIP1 in the mutant yeast rescues only weak growth whereas that of OsSIP2 had no effect at all (Fig. 3F). This observation demonstrates that OsSIP1 may have limited methylamine permeability and OsSIP2 is likely a strict water channel. Another possibility was not excluded that both SIP proteins are not sufficiently localized to the plasma membrane.
Expression pattern of OsSIP1 and OsSIP2
Tissue expressing the highest level of OsSIP1 was found in leaves of seedling and at reproductive stage, with moderately reduced level at tillering stage. OsSIP1 transcripts were detected at roughly similar levels in most tissues, but were reduced in root (Fig. 5A). Much lower expression of OsSIP2 was found in almost all tissues except anther, where it reached maximum transcript levels similar to those of OsSIP1 in young leaves (Fig. 5A).
To further investigate the tissue-specific expression of OsSIP1 and OsSIP2 in rice, we generated transgenic expressing the β-glucuronidase (GUS) reporter gene under the control of the native promoters (pOsSIP1::GUS and pOsSIP2::GUS) (Fig. S1). Histochemical GUS staining was carried out on germinating seeds, 2-week-old seedlings, and reproductive glumes of pOsSIP1::GUS and pOsSIP2::GUS transgenic plants (homozygous T2 progeny), and the wild-type plant served as a negative control of staining (Fig. S2). It was observed that OsSIP1 displays stronger GUS signals than OsSIP2 in all rice tissues tested, excepting the anthers. Stronger GUS staining signals were found in both OsSIP1 and OsSIP2 at the shoot tips of germinating seeds (Fig. 5B a, m) and at the stem tips of seedlings (Fig. 5B b, n). Additionally, expression of OsSIP1 was detected in adventitious roots of seedlings (Fig. 5B c), particularly at root hairs and root tips (Fig. 5B d, e). On contrary, no OsSIP2 expression signals were detected in the roots (Fig. 5B n, o, p), suggesting that OsSIP2 was hardly expressed in rice roots. In stem, leaf sheath and glumes, both OsSIP1 and OsSIP2 produced visible staining signals, with OsSIP1 displaying much higher intensity than OsSIP2 (Fig. 5B f, g, h, l, and r, s, t, x). In addition, less staining signals were detected in leaf blades and veins for OsSIP1, but not for OsSIP2 (Fig. 5B i, j, u, v). Interestingly, OsSIP1 expression signal was detected at the base of filaments and pistils (Fig. 5B k), while OsSIP2 displayed a robust GUS staining signal in anthers (Fig. 5B w). Overall, OsSIP1 showed variable expression in all tissues, whereas OsSIP2 showed high expression specific to anthers. The result of the promoter-GUS experiment was in good agreement with that obtained by qRT-PCR detection.
Expression profiles of OsSIPs during biotic and abiotic stresses, as well as exogenous phytohormone treatment
AQPs play a critical role in plant responses to various stressors, including drought, salt, temperature, redox, pathogen attack, as well as hormonal stimuli [4, 15, 20, 21]. To determine the putative cis-acting elements in promoter of the OsSIP genes, analysis of the 2000 bp region upstream of the initiation codon (ATG) of the two OsSIP genes were performed using the PlantCARE database (http://bioinformatics.psb.ugent.be/webtools/plantcare/html/). The identified elements are summarized in Table. S1. In addition to transcription factor recognition sites, such as those for MYC and MYB, an extensive repertoire of light-responsive elements, meristematic expression elements (CAT-box, CCGTCC-box), gibberellin-responsive elements (P-box, TATC-box), anaerobic induction elements (ARE), drought stress-responsive elements (MBS), and several stress-responsive elements (STRE, WRE3, W-box) are predicted in the promoters of both OsSIPs genes. The promoter of the OsSIP1 gene also contains response elements for a salicylic acid (SA), methyl jasmonate (MeJA), ethylene, and low temperature (LTR). Similarly, the promoter of the OsSIP2 gene shows the presence of response elements for abscisic acid (ABA) and auxin (Table. S1). Overall, OsSIP1 contains a greater number of response elements in its promoter than OsSIP2.
We then conducted experiments using 2-week-old rice seedling challenged by salt and osmotic stresses and showed that both OsSIP1 and OsSIP2 were up-regulated (Fig. 6). However, the expression of OsSIP1 and OsSIP2 was significantly reduced in the condition in which water was withdrawn to cause dehydration. OsP5CS1 and OsMYB48-1, which are known to respond positively to dehydration, salt stress, and osmotic stress [52, 53], were used as treatment controls (Fig. 6). On the other hand, cold and high temperature stresses caused a shift in expression of OsSIP1 and OsSIP2 during the first three hours after treatment, which was similar to that of the known stress-responsive marker genes, namely OsCBF3 and OsHSP14.7 [54, 55] (Fig. 7). In addition, oxidative stress imposed by hydrogen peroxide treatment up-regulated the expression of OsSIP1 and OsSIP2 as well as an oxidative stress marker gene OsNAC066 [56] (Fig. 8). Similarly, under DTT-induced ER stress, the expression of both OsSIP1 and OsSIP2 was slightly increased, in contrast to the strong response of the marker gene OsbZIP50 [57] (Fig. 8). It was concluded that, in terms of quantitative comparison, OsSIP1 was most responsive to H2O2 and salt treatments while OsSIP2 to heat and DTT stresses.
When seedlings were treated with gibberellin (GA), abscisic acid (ABA), methyl jasmonate (MeJA) or salicylic acid (SA), both OsSIP1 and OsSIP2 showed up-regulated expression to some extent (Fig. 9). We also checked the expression of the rice GA synthesis gene OsGA20ox2 [58], and marker genes for hormone signaling response, namely OsABF1, OsLOX2 and X58877 [59–61] for ABA, MeJA and SA, under the same conditions (Fig. 9). These results suggest that OsSIP1 and OsSIP2 are little regulated by the tested phytohormone signaling pathways.
Finally, to explore the potential involvement of OsSIPs genes in the response to pathogen infestation, we spray-inoculated Magnaporthe oryzae strain Guy11 on rice seedlings at the three-leaf stage and collected samples of diseased leaves every 24 hours over a period of 5 days (Fig. 10). Both OsSIP1 and OsSIP2 showed moderate upregulation in response to M. oryzae infection, with OsSIP1 being stronger. The OsNAC4 gene, used as a control marker, showed increased transcript levels in the infected leaves (Fig. 10).