Numerous studies have reported the identification of functional genes that are relevant to plant stress tolerance and that can be used for crop improvement[18-21]. To aid efforts to breed salinity tolerance in soybean, we analyzed the orthologs of plasma membrane intrinsic proteins (PIPs) in soybean and tested their interactions using a yeast two-hybrid system. The interaction between PIP1 and PIP2 proteins functions as the main signal for cell membrane water and salt exchange and PIPs have been considered as functional units that perform their physiological roles under different environmental stresses, such as salt and drought stress. Aquaporins play an important role in growth regulation of plants by influencing root water uptake and leaf gas exchange. The soybean PIP1, GmPIP1;6, has previously been well characterized. Its function in growth regulation and salt tolerance was analysed by constitutive overexpression [22]. A separate study highlighted the involvement of several aquaporin homologs in response to a variety of environmental stressors that interrupt plant cell osmotic balance [23].
PIP1 and PIP2 have highly conserved peptide sequences, and the main differences between them are the lengths of their N and C terminal ends [5]. In soybean, the N terminal ends of GmPIP1 are longer than that of GmPIP2 (approximately 15 amino acids), however, the C terminal ends of GmPIP1 type are shorter than that of GmPIP2 type (approximately 8 amino acids). Interestingly, the transmembrane region of both GmPIP1 and GmPIP2 nearly share the same section (Supplementary Figure 1).
Ispolatov et al. [24] proposed that duplicated proteins were more likely to interact among themselves than with other proteins, and that paralogous interactions were inherited from ancient homo-dimeric proteins, rather than established de novo after gene duplication. In evolution progress, gene duplication events increase gene number by tandem- and segmental-duplication [25]. Soybean (Glycine max (L.) Merr.) is an important crop and well-studied. Previous research reported that soybean is a paleopolyploid, and at least two rounds of large-scale duplication occurred in its ancestral genome at approximately 14- to 42- million ago [26]. In this study, we searched the database PGDD (http://chibba.agtec.uga.edu/duplication/), and found many duplication events in the GmPIP family. Expansion of aquaporin gene families via genome duplication events have been reported in other plants [27].
The transcriptional profiles of PIPs may provide evidence for PIP1–PIP2 interactions. For example, a joint increase (or decrease) in the expression of specific PIP1–PIP2 pairs in plants under stress may indicate shared functionality. Thus, the formation of hetero-tetramers composed of specific PIP1s and PIP2s could be affected by their mRNA abundance [5]. Transcriptional profiles in rice, maize, and Arabidopsis indicate interactions between PIP1–PIP2 pairs in these species [5]. Additionally, Zargar et al. [28] developed a gene co-expression network of rice aquaporin genes (OsPIPs) and tonoplast intrinsic proteins (OsTIPs) using Rice Friend server (http://ricefrend.dna.affrc.go.jp). They found co-expression of PIP1–PIP2 pairs, indicating likely physical interaction between these proteins. In this study, we found physical interactions among GmPIPs using yeast two-hybrid assays. We detected both homotetramers and heterotetramers among these proteins. Salt, but not mannitol, enhanced these interactions. Only GmPIP2;9 could form homotetramers, and this interaction was enhanced by salt stress, but weakened by osmotic stress. Our results corroborate those of Bienert et al. [29], who reported PIP heterotetramerization under salt stress in Selaginella moellendorffii.
The TMDs of GmPIPs were predicted by SMART software. All GmPIPs contained six TMDs. GmPIP1;5 and GmPIP1;6 differed by two amino acids in TMD2 (D/Y) and TMD6 (H/Q). However, there were many differences in TMD2, TMD4, and TMD6 among GmPIP2;4, GmPIP2;6, GmPIP2;8, GmPIP2;9, GmPIP2;10, and GmPIP2;11. This indicates that the highly conserved sequences in TMD1, TMD3, and TMD5 may play a crucial role in the formation of PIP1–PIP2 pairs between GmPIP1;5 and GmPIP1;6, and among GmPIP2;4–GmPIP2;11. Using extensive amino acid substitution mutagenesis, Yoo et al. (2016) studied tetramer formation in ArabidopsisAtPIP2;1. They demonstrated that TMD1, TMD2, and TMD5 contained essential amino acid residues key to tetramer formation.
In addition, the expression profiles of GmPIPs under salt stress also showed similar expression patterns in GmPIP1;5 and GmPIP1;6, and in GmPIP2;3, GmPIP2;4, GmPIP2;5, GmPIP2;6; GmPIP2;9, GmPIP2;10, GmPIP2;11; GmPIP2;13 and GmPIP2;14 respectively, which is partly consistent with the sequence similarities and interaction patterns among GmPIPs. The expression of GmPIP1;3, GmPIP1;4, GmPIP2;1, GmPIP2;8, GmPIP2;9, GmPIP2;10, GmPIP2;11 was significantly up regulated by salt stress. However, the other GmPIPs were significantly downregulated by salt stress, except GmPIP1;9, 1;10, 2;7 and 2;12.
Aquaporins are implicated in a variety of stress responses that disturb plant cell osmotic balance and nutrient homeostasis (Yaneff et al. 2015)[5]. They are involved in the Arabidopsis response to drought stress (Afzal et al., 2016)[23], in leaves and roots of sugar beet under salt stress (Lv et al., 2018)[30], and in rice tolerance to salt stress and cold stress (Qiang et al. 2015)[15]. We have discovered that some GmPIPs are significantly up- or down-regulated in response to salt stress. GmPIP1;6 is the closest ortholog to ArabidopsisAtPIP1;2 (with an amino acid sequence identity of 84.8%), which localizes to the Golgi apparatus and the membrane system. AtPIP1;2 was considered as a functional water channel when it was expressed alone in Xenopus oocytes [31]. A previous study also suggested that many GmPIPs have high sequence similarity but diverse functions (Zhang et al. 2013)[2]. GmPIP1;6 could interact with other GmPIP2 type aquaporins, which implied they may play crucial roles in aquaporin trafficking from the Golgi apparatus to the membrane system in plants. The functions and molecular mechanisms of the diverse families of plant PIPs still need further study.