Scholarly article on topic 'Divergence in function and expression of the NOD26-like intrinsic proteins in plants'

Divergence in function and expression of the NOD26-like intrinsic proteins in plants Academic research paper on "Biological sciences"

Share paper
Academic journal
BMC Genomics
OECD Field of science

Academic research paper on topic "Divergence in function and expression of the NOD26-like intrinsic proteins in plants"

BMC Genomics

BioMed Central

Open Access

Research article

Divergence in function and expression of the NOD26-like intrinsic proteins in plants

Qingpo Liu*1, Huasen Wang1, Zhonghua Zhang2, Jiasheng Wu3, Ying Feng4 and Zhujun Zhu*1

Address: School of Agriculture and Food Science, Zhejiang Forestry University, Lin'an, Hangzhou 311300, PR China, 2Institute of Vegetables and Flowers, Chinese Academy of Agricultural Sciences, Beijing 100081, PR China, 3School of Forestry and Biotechnology, Zhejiang Forestry University, Lin'an, Hangzhou 311300, PR China and 4College of Environmental and Resources Science, Zhejiang University, Hangzhou 310029, PR China

Email: Qingpo Liu* -; Huasen Wang -; Zhonghua Zhang -; Jiasheng Wu -; Ying Feng -; Zhujun Zhu* - * Corresponding authors

Published: 15 July 2009 Received: 19 March 2009

BMC Genomics 2009, 10:313 doi:l0.ll86/l47l-2l64-l0-3l3 Acœpted I5 Ju|y 2009

This article is available from: © 2009 Liu et al; licensee BioMed Central Ltd.

This is an Open Access article distributed under the terms of the Creative Commons Attribution License (, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.


Background: NOD26-like intrinsic proteins (NIPs) that belong to the aquaporin superfamily are plant-specific and exhibit a similar three-dimensional structure. Experimental evidences however revealed that functional divergence should have extensively occurred among NIP genes. It is therefore intriguing to further investigate the evolutionary mechanisms being responsible for the functional diversification of the NIP genes. To better understand this process, a comprehensive analysis including the phylogenetic, positive selection, functional divergence, and transcriptional analysis was carried out.

Results: The origination of NIPs could be dated back to the primitive land plants, and their diversification would be no younger than the emergence time of the moss P. patens. The rapid proliferation of NIPs in plants may be primarily attributed to the segmental chromosome duplication produced by polyploidy and tandem duplications. The maximum likelihood analysis revealed that NIPs should have experienced strong selective pressure for adaptive evolution after gene duplication and/or speciation, prompting the formation of distinct NIP groups. Functional divergence analysis at the amino acid level has provided strong statistical evidence for shifted evolutionary rate and/or radical change of the physiochemical properties of amino acids after gene duplication, and DIVERGE2 has identified the critical amino acid sites that are thought to be responsible for the divergence for further investigation. The expression of plant NIPs displays a distinct tissue-, cell-type-, and developmental specific pattern, and their responses to various stress treatments are quite different also. The differences in organization of c/s-acting regulatory elements in the promoter regions may partially explain their distinction in expression.

Conclusion: A number of analyses both at the DNA and amino acid sequence levels have provided strong evidences that plant NIPs have suffered a high divergence in function and expression during evolution, which is primarily attributed to the strong positive selection or a rapid change of evolutionary rate and/or physiochemical properties of some critical amino acid sites.


The aquaporins are a family of small pore-forming integral membrane proteins, which contain six membrane-spanning a-helices, with the N- and C-termini facing the cytosol [1]. On the basis of the first sequenced member -the Major Intrinsic Protein of bovine lens cells [2], this family is also referred to as the MIP superfamily. The MIP proteins are found to be present in all kingdoms from archaea to plants and animals [1]. However, MlPs constitute a larger and more diverse family in plants than in animals. There are 35 and 39 MIP genes in the genomes of the model plants Arabidopsis [3] and rice [4] respectively. By contrast, mammals only possess 13 distinct MIPs [5]. According to sequence similarities, plant aquaporins are clearly classified into five major subfamilies: plasmamembrane intrinsic proteins (PIPs), tonoplast intrinsic proteins (TIPs), NOD26-like intrinsic proteins (NIPs), small basic intrinsic proteins (SIPs), and the GlpF-like intrinsic proteins (GIPs). The divergence of plant aquaporins into five subfamilies had been already established as early as the time of emergence of primitive terrestrial plants [6,7]. The last subfamily of MIP genes has so far only been identified in two mosses [8]. Within each of the other subfamilies, MIP genes can be further subdivided into more than two groups [9] that may correspond to their localization and transport selectivity [10].

There are strong evidences that aquaporins are central components in plant water relations [10,11]. The significances of plant aquaporins functioning in mediating water and/or other small solutes across biomembrane and response to external environmental stresses have been widely reported [12-14]. It is likely that the particular abundance of MIP genes in plants may be attributed to the higher degree of compartmentalization of plant cells and their greater necessity for fine-tuned water control [3]. Alternatively, the extensive proliferation of aquaporin iso-forms may offer an adaptive advantage for plants to grow in different environmental conditions, possibly as a result of divergent transport selectivities or regulatory mechanisms [11].

NIPs that were defined as NOD26-like intrinsic proteins on the basis of the archetype nodulin26 protein identified firstly in soybean [15], are unique to plants. These proteins are presumed to be involved in exchange of metabolites between the host and the symbiont [1]. Nevertheless, NIPs are widely distributed in both leguminous and non-leguminous plants, indicating that plant NIP function is not limited to the role that they play in nodule symbiosis [16]. Of plant aquaporins, only proteins belonging to this subfamily have glycerol transport activity [12]. Thus this may suggest that the common ancestor of plant aquaporins had lacked the ability to transport glycerol and later on NIPs had acquired this transport activity during evolution

to compensate for the absence of GLPs in plants [17]. By reconstructing phylogenetic trees, Zardoya et al. [17] suggested that NIPs might originate from bacteria at the origin of plants through a single horizontal gene transfer event. Zardoya [9] confirmed the functional recruitment of NIPs to glycerol transport, however the gene horizontal transfer origination of NIPs could not be effectively recovered because of the general lack of resolution of deeper nodes.

The Arabidopsis thaliana, Oryza sativa, and Physcomitrella patens genomes encode 9, 13, and at least 5 NIP proteins respectively [3,4,7]. The occurrence and increase of NIP genes in plants indicate a wider range of function that may include a greater range in selectivity [12,18]. Functional analysis of NIPs indeed has revealed diversity in their transport substrates. The soybean nodulin26 protein not only transports water but also glycerol, formamide, malat, and NH3 [12,19,20]. In Arabidopsis, the AtNIP2;1 shows minimal water and glycerol transport, but displays transport of lactic acid with a preference under the anaerobic condition [21], whereas AtNIP5;1 is essential for efficient boron uptake and plant development under boron limitation [22]. Two NIPs (OsNIP2;1 and OsNIP2;2) in rice show a transport activity of larger solute silicic acid that enhances resistance of plants to biotic and abiotic stresses

[23], whereas OsNIP1;1 and OsNIP3;1 do not; furthermore, the OsNIP2;1 can be permeable to water, urea, as well as boric acid, but not glycerol [24]. Therefore, the substrate specificity is NIP-dependent. Nonetheless, other unknown factors or structural features may also be involved in the process of efficient substrate recognition

[24]. Very recently, some NIPs, including AtNIP1;1, AtNIP1;2, AtNIP5;1, AtNIP6;1, AtNIP7;1, OsNIP2;1, OsNIP2;2, OsNIP3;2, and LjNIP5;1, LjNIP6;1 were found to be responsible for the permeability to arsenite [25-28]. It can be proposed thus that arsenite transport through NIPs should be a conserved and ancient feature.

The molecular basis of aquaporins selectivity is pivotally due to two filters within the pore; the first one is formed by the conserved dual NPA motif, and the second one formed by a constriction region that is also called the ar/R (aromatic/arginine) filter [11,29,30]. However, it seems that the NPA motif is not crucial for [30,31], while the ar/ R filter plays an important role in determining the substrate selectivity for the NIP subfamily [24]. The ar/R filter is located in the narrowest region on the extra-membrane mouth of the pore, and formed by four residues, one each from helix 2 (H2) and helix 5 (H5), as well as two residues from loop E (LE1 and LE2) [16]. It appears that the properties of the four residues making up the ar/R selectivity filter govern the substrate specificity of the pore [23], and are thought to be useful for predicting the function of the proteins [30]. Based on the ar/R regions of aquaporins,

NIPs can be divided into three distinct groups [32]. NIP I proteins in Arabidopsis have been reported to transport water, glycerol, and lactic acid, whose ar/R region is characterized by Trp (W), Val (V)/Ile (I), Ala (A), and Arg (R). The consensus of the ar/R region for NIP II proteins, which are permeable to larger solutes than NIP I protein, is Thr (T)/Ala (A), Ala (A)/Ile (I)/Val (V), Gly (G)/Ala (A), and R [23,24]. OsNIP2;1 and OsNIP2;2 belong to the third group (NIP III) [24]. The ar/R region of NIP III consists of Gly (G), Ser (S), Gly (G), and Arg (R), forming a larger constriction size (>6A) compared with other NIP groups (<5A and 3.5A), which allows solutes like silicic acid with a larger diameter (4.38A) to permeate [24,3133]. The critical role of the ar/R region is further demonstrated by the observations that substitution of Trp (W) with His (H) at the position H2 in LIMP2 abolished its glycerol transport when expressed in Xenopus oocytes [34]. Moreover, mutation at the Arg (R) residue in LE2 can cause human disease [35]. This residue is strictly conserved in NIPs, and thought to be important for providing hydrogen bonds for transport of water or glycerol molecules and to repel cations from the pore [36].

With a few of exceptions, the subcellular location of most members from the NIP subfamily is still uncertain. The archetype NIP, nodulin26 is located in the peribacteroid membrane (PBM) of soybean nodule cells [37]. The Lotus japonicus NIP gene, LIMP2 is also probably located in the PBM [38], whereas Arabidopsis AtNIP1;1 [28], AtNIP5;1 [22], rice NIP2 genes (OsNIP2;1 and OsNIP2;2) and the barley HvLsi1 are found to be localized on the plasma membrane [39-41]. In addition, it was found that the Arabidopsis AtNIP2;1 is predominantly expressed in young roots and is mainly located to the endoplasmic reticulum membrane [42]. Accordingly, the subcellular localization of each NIP protein may be diverse.

It has been reported that the activity of both plant and animal aquaporins may be regulated by phosphorylation [11]. The soybean NOD26 can be phosphorylated by a calcium-dependent protein kinase (CDPK) [37]; phosphorylation of NOD26 on Ser262 enhanced its water permeability [43]. The phosphorylated Ser262 of NOD26 is conserved in most but not all NIPs from Arabidopsis [44]. Experiments with the C-terminal extension of the Arabi-dopsis AtNIP7;1 showed that this gene could be phosphorylated by activated AtMPK4 in vitro [16]. Mitani et al. [24] however found no evidence for the involvement of phosphorylation in the regulation of OsNIP2;1, because neither okadaic acid nor K252a affected the transport activity of this protein for silicic acid in oocytes.

Although NIP proteins exhibit a similar three-dimensional structure [32], functional divergence has extensively occurred among NIP groups. It is therefore

intriguing to further investigate the evolutionary mechanisms driving the functional diversification of the NIP groups at the nucleotide and amino acid sequence level respectively. Here we showed that strong positive selection had occurred after gene duplication and/or specia-tion. Functional divergence analysis provided convincing evidence for shifted evolutionary rate and/or rapid changes of amino acid properties between NIP groups, and identified some critical amino acid sites that are thought to be significantly functional divergence related. Further, we investigated the expression pattern and the cis-acting regulatory element organization in plant NIPs also.

Results and discussion

Phylogenetic and sequence character analysis of the plant NIP subfamily

On the basis of sequence similarity and the conserved MIP domain, we have identified 6, 8, 9, 11, 6, and 7 NIP genes from the Vitis vinifera, Populus trichocarpa, Sorghum bicolor, Glycine max, Cucumis sativus, and Pinus taeda genomes respectively (Additional file 1). Danielson and Johanson [7] revealed that the moss Physcomitrella patens encodes at least five NIP proteins, of which PpNIP6;1 is partial and excluded from the phylogenetic analysis. Furthermore, although Chlamydomonas reinhardtii might lack NIP proteins [45], one NIP homolog was identified from the green alga Ostreococcus lucimarinus in JGI (protein ID 25291), indicating that the NIP subfamily would be more ancient.

More than eighty NIP sequences were collected, aligned, and used to reconstruct phylogenetic trees. Figure 1 shows that plant NIPs are clearly divided into three groups, which supports the classification of NIPs based on the ar/ R region [24,32]. It is apparent that each NIP group possesses at least one NIP protein each from P. taeda and P. patens (Figure 1), with the exception of the NIP I group where no P. patens NIPs are represented thus far, suggesting that the most recent common ancestor (MRCA) of P. patens and higher plants should have three ancestral NIP genes corresponding to the 3 groups. The diversification of the three ancient progenitor NIPs should have occurred before the emergence of the moss P. patens but after the green alga O. lucimarinus.

The diverse NIPs should have arisen from three ancestral genes, and subsequently species-specific expansion of this subfamily has occurred to a great extent. For example, three P. patens NIP5 (PpNIP5;1 - PpNIP5;3) and four P. taeda NIP1 isoforms (PitNIP1;1 - PitNIP1;4) were tightly clustered together in the NIP III and NIP I group respectively (Figure 1). Within each group, the monocot and dicot NIPs have formed distinct clades, indicating that the extensive proliferation of NIPs should have occurred after the monocot-dicot split. This proliferation may be attrib-

Figure 1

Phylogenetic tree of NIP proteins in plants. The

number beside the branches represents bootstrap values > 60% based on 1000 resamplings. The scale is in amino acid substitutions per site. The two NPA motifs and the four residues making up the ar/R filter are indicated in the figure. To identify the species of origin for each NIP gene, a species acronym is included before the gene name: An, Atriplex num-mularia; At, Arabidopsis thaliana; Ca, Cicer arietinum; Cp, Cucurbita pepo; Cs, Cucumis sativus; Gm, Glycine max; Lj, Lotus japonicus; Mt, Medicago truncatula; Os, Oryza sativa; Pit, Pinus taeda; Pp, Physcomitrella patens; Pt, Populus trichocarpa; Sb, Sorghum bicolor; Vu, Vigna unguiculata; Vv, Vitis vinifera; Zm, Zea mays.

uted to the duplication of individual genes, or even the entire genome. The rice and Arabidopsis genomes have undergone two and three rounds of whole genome duplications (WGD) respectively [46,47], which may be responsible for the expansion of the NIP subfamily in the two model plants. We searched against the segmentally duplicated chromosome segments identified based on the TIGR rice ( segmental_dup/) and Arabidopsis ( tdb/e2k1/ath1 /Arabidopsis_genome_duplication.shtml) genome annotation release version 5.0 respectively. As expected, three gene pairs including OsNIP1;1/OsNIP1;4, OsNIP1;2/OsNIP1;3, and OsNIP2;1 /OsNIP2;2 were found to be located into chromosomal regions that were supposed to have undergone large-scale segmental duplications. Moreover, three genes (OsNIP3;2, OsNIP3;3, and OsNIP3;5) were arranged as tandem on chromosome 8. As observed in rice, we also found evidence for segmental duplication (AtNIP3;1) and tandem duplications (AtNIP1;1/AtNIP1;2, and AtNIP4;1 / AtNIP4;2) in Arabidopsis. The duplicated copy for the AtNIP3;1 gene however should have lost during evolution in that only one copy of it was found in the corresponding segmental duplicated regions. Thus, the explosive expansion of rice and Arabi-dopsis NIPs would attribute primarily to the segmental duplications and tandem duplications.

Furthermore, the exon/intron lengths and gene structure for each NIP gene were examined (Additional file 1), since the exon/intron structures were terribly important for investigating the evolution of any gene family. It was found that most of the NIP genes possess four or five exons, and only a few of them have 2, 3, or 6 exons. Moreover, we observed that with one exception (OsNIP1;2/ OsNIP1;3), the gene pairs predicted to be produced by duplications showed similar gene structures, even though the lengths of introns might be different. For example, both AtNIP4;1 and AtNIP4;2 have 5 exons, and the length is the same for each of the five exons too. However, the lengths of the four introns vary extensively (Additional file 1). The similar cases were also found in the OsNIP2;1/ OsNIP2;2, OsNIP1;1/OsNIP1;4 gene pairs, and in the OsNIP3;2/OsNIP3;3/OsNIP3;5 gene cluster. Although OsNIP1;3 lacks one exon than OsNIP1;2, the length of the second exon of the former gene is nearly equal to the total length of the second and third exons of the latter one. Besides, the lengths of the other three exons are similar to each other too (Additional file 1). It could be thus inferred that OsNIP1;2 should have independently gained an extra intron during evolution, because almost all NIP I genes in monocot plants possess four exons.

The two major constriction filters corresponding to the two NPA motifs and the ar/R filter were shown in Figure

1. It was observed that NIP proteins often have unorthodox NPA motifs, which is different from other MIPs. In the NIP subfamily, the first and second NPA motifs are always replaced by NPS and NPV (NPG), respectively (Figure 1, and Additional file 1). In the NIP I and III groups, the first NPA motif is highly invariant, and the second NPA motif is also conserved besides several exceptions. However, more complicated NPA patterns were present in the NIP II group. There are two types of NPA motifs each in the first (NPA and NPS) and second (NPA and NPV) positions. The monocot NIP3 and dicot NIP7 proteins possess conserved NPA in the second position, while it is replaced by NPV in other NIP II proteins (Figure 1). Wallace and Roberts [31] demonstrated that in AtNIP6;1 substitution of Ala for Val in the NPA2 region did not alter its transport selectivity, indicating that the NPA motifs should have little effect on determining the transport specificity for this protein.

The ar/R region is the second constriction filter. It is obvious that this region is highly subgroup specific for the NIP I and III groups. With one exception (AtNIP3;1), the ar/R region is characterized by W, V, A, R for the NIP I group. The hydrophobic residues W, V, and A, provide a nonpolar surface that hydrophobically interacts with the hydrocarbon skeleton of glycerol, whereas the R residue creates a hydrogen bonding with the hydroxyls of glyc-erol, which facilitates the transit of glycerol and water molecules across the central water channel [32]. In NIP III, monocot and dicot NIPs possess a specific ar/R filter consisting of G, S, G, and R. The tiny character of G and S residues insures a wider aperture for the ar/R filter [32]. Nonetheless, there are some exceptions. Compared with CsNIP2;1, CsNIP2;2 possesses a different ar/R filter where the tiny Gly (G) residue is replace by the bulky Cys (C), and this ar/R filter may offer a narrower aperture than the former gene, indicating that the two genes may perform different transporter functions. The similar case was also found in PitNIP2;1, where the slight larger Ser (S) residue substitutes for the first Gly (G). The moss P. patens PpNIP5s have a unique combination of amino acids at the ar/R filter (Figure 1, see also [7]), where the four residues (F, A, A, and R) constitute a mediate aperture that should be larger than the NIP I protein, but narrower than the NIP III protein. A considerable variation in the ar/R region exists in NIP II proteins. The size and hydrophilic-ity of the four residues are different from each other, suggesting that members belonging to this group should function in different ways. Notwithstanding, the NPA motifs (NPA/NPV) and the ar/R filter (A, I, A, R) in the moss PpNIP3;1 is identical to AtNIP6;1, suggesting that these genes had the signatures before the split of bryo-phytes and vascular plants [7]. However, whether PpNIP3;1 functions as a boric acid transporter like AtNIP6;1 [48] needs further experimental investigation.

In addition, we found that except for the OsNIP3;2/ OsNIP3;3/OsNIP3;5 gene cluster, other putative gene duplicates in rice and Arabidopsis possess the same NPA motifs and the ar/R signature (Additional file 1). The AtNIP4s (4;1 and 4;2) have the similar gene structure and the two constriction filters, suggesting that the two genes should perform same or similar functions in transporting water and glycerol, etc. In rice, both OsNIP2;1 and OsNIP2;2 were demonstrated to be permeable to the larger solute silicic acid [39,40], which supports our speculation on the AtNIP4s. Given that the ar/R regions in OsNIP3;2, OsNIP3;3 and OsNIP3;5 are clearly differentiated, we speculated that functional divergence should have occurred among them. Moreover, we found that the putative orthologs in the NIP I and III groups have the same (or with a slight modification) NPA motifs and the ar/R signature, whereas the orthologs in the NIP II group are different from each other (Figure 1). It was worth to note that the corresponding motifs in the green algae O. lucimarinus, are NPS/NAA, and Y, L, G, R respectively. If this is not the result of sequencing error, it can be inferred that this gene should play a quite different role from other NIPs in moss and higher plants.

In Arabidopsis, as much as seven NIP groups can be defined [3], whereas there are only three groups in rice [4,49]. This is a consequence of the large variation in NIP sequence divergence in dicot species [3]. It is worth noting that there is no Arabidopsis NIPs classified into the NIP III group, although the NIP2 genes of some dicots such as poplar, grape, and cucumber belong to this group (Figure 1). This suggests that Arabidopsis should have lost its counterpart in this group during evolution. Alternatively, it is likely that extensive sequence and structural variation should have occurred in Arabidopsis AtNIP2;1, leading to its functional divergence from other plant NIP2 genes. Accordingly we examined the gene structure and ar/R filter for this gene, and found that AtNIP2;1 possesses four exons and its ar/R filter is W, V, A, R, which is quite different from other NIP2 genes that always have five exons and whose ar/R filter is composed of G, S, G, R (see Additional file 1).

Detection of positive selection signal

The multiple alignments of NIP protein sequences were shown in Additional file 2. The estimation of positive selection was based on the corresponding nucleotide ML tree (Additional file 3), where the sequences in partial and whose length are less than 240 codons were excluded so as to avoid possible analysis biases. To test for variable ra ratios among lineages, we conducted the likelihood ratio test (LRT) to compare the two extreme models: the one-ratio model that assumes a unique rate ratio for all branches, and the free-ratio model that assumes an independent ra ratio for each branch [50,51]. The log likeli-

hood value under the one-ratio model is -23356.6, while the value is -23007.7 for the free-ratio model. Twice the log likelihood difference, 2A€ = 697.8, is strongly statistically significant, revealing a heterogeneous selective pressure among lineages. It is obvious that some branches of the NIP phylogeny include some internal branches having ra > 1, showing strong evidence for adaptive evolution. Moreover, we observed that the branches that were detected to be under positive selection should correspond to gene duplication or speciation events (see Additional file 3), which might act as a major evolutionary force driving the divergence of NIP functions.

The divergence of family members may involve positive selection, as indicated by many typical studies [52]. In a gene family, the fate of new genes produced by duplication would either evolve a new function under positive selection, or be lost during evolution [53]. Plants have evolved more abundant aquaporins with multifunctions, which may reflect the need for plants to better adapt external environmental conditions. Evidence for adaptive evolution thus clearly imply that functional diversification of the NIPs represents an evolutionary advantage for undergoing ecological adaptation to local environment [11]. It thus appears that continued positive selection should have acted on the NIP subfamily during evolution and this selection would be remarkably significant.

Functional divergence analysis (FDA) of plant NIP proteins

Two types of functional divergence (Type-I and Type-II) between gene clusters of the NIP subfamily were estimated by posterior analysis using DIVERGE2 that evaluates shifted evolutionary rate and altered amino acid property after gene duplication [54,55]. The advantage of these methods is that they use amino acid sequences, and thereby is not sensitive to saturation of synonymous sites. The NIP subfamily consists of three major groups (NIP I, II, and III) ([24,32]; see also Figure 1). However, the results of positive selection analysis and primary FDA support the classification of NIP IIA and NIP IIB as two dis-

tinct groups (see Additional file 3). Thus, herein four gene clusters of interests were used as input for the DIVERGE2 analysis. As expected, the null hypothesis (no functional divergence) could be strongly rejected in that the coefficients of Type-I functional divergence (9t) between NIP groups were statistically significant (p < 0.01; Table 1), indicating that significant amino acid site-specific selective constraints operate on different types of NIP members leading to a subgroup-specific functional evolution after their diversification. Further, a functional distance analysis was conducted. We estimated the functional branch length (bF) for each group by employing the least-squares method [54], and found that the level of altered selective constraints of group genes, measured by this index, followed bF (II, 0.48) > bF (IV, 0.32) > bF (I, 0.19) > bF(III, 0.12), suggesting that the NIP II group should be significantly divergent in function from other groups. On the other hand, it was found that the coefficients of Type-II functional divergence (9n) between I/II, I/III, and I/IV were insignificant (p > 0.05). Nevertheless, we found evidences for Type-II functional divergence between three group pairs including II/III, II/IV, and III/IV, indicative of a radical shift of amino acid property [55].

To identify critical amino acid sites that may be responsible for functional divergence between NIP groups, the posterior probability (Qk) of divergence was determined for each site. According to the definition, large Qk indicates a high possibility that the evolutionary rate or amino acid physiochemical property of a site is different between two clusters. DIVERGE2 thus identified some critical amino acid sites (CAASs) that are highly relevant to functional divergence (see Additional file 4, and Table 1). In order to extensively reduce positive false, Qk > 0.8 was empirically used as cutoff to identify the Type-I and Type-II functional divergence-related residues between gene clusters. The results showed that more than 13 CAASs were supposed to be responsible for the functional divergence between NIP II and I, III, IV respectively; whereas there were only two, seven, and one CAAS with Qk > 0.8

Table 1: Functional divergence between groups of the plant NIP subfamily

Group! Group2 Type-I Type-II

9 ± S.E. LRT p Qk > 0.8 Qk > 0.9 9 ± S.E. Qk > 0.8

NIP I II 0.449 ± 0.063 5 .3 <0.0! !7 8 -0.085 ± 0. !63 0

NIP I III 0.190 ± 0.068 7.7 <0.0! 2 0 -0.!05 ± 0. !75 0

NIP I IV 0.429 ± 0.087 24.5 <0.0! 7 3 0.004 ± 0.!78 0

NIP II III 0.486 ± 0.068 50.6 <0.0! 25 !! 0.!06 ± 0.!!2 2

NIP II IV 0.534 ± 0.098 29.4 <0.0! !3 3 0.!3! ± 0.!!5 2

NIP III IV 0.286 ± 0.092 9.7 <0.0! ! ! 0.22! ± 0.!!9 7

Note: 9 and 9N, the coefficients of Type-I and Type-II functional divergence between two gene clusters; LRT, Likelihood Ratio Statistic; Qk, posterior probability. Large Qk value indicates a high possibility that the functional constraint (or the evolutionary rate) or the physicochemical properties of a given amino acid site is different between two clusters.

identified between I/III, I/IV, and III/IV respectively. Interestingly, in contrast to the Type-I functional divergence, two, two, and seven Type-II related CAASs were identified for the II/III, II/IV, and III/IV pairs respectively (Table 1). Compared with only one CAAS for the Type-I functional divergence, there were seven predicted sites for the Type-II functional divergence between III/IV, indicating that the functional divergence between the two group genes was mainly attributed to the rapid changes of amino acid physiochemical property, and secondly to the shifted evolutionary rate. The contrary cases were found for the II/I, II/III, and II/IV pairs. These observations indicated that site-specific shift of evolutionary rate and changes of amino acid property should not uniformly act on the NIP subfamily members during long periods time of evolution. The relative importance of Type-I and Type-II functional divergence may be associated with specific functional classes of this protein family. However, the degree of functional divergence between NIP I/III was not remarkably significant, because there were only two CAASs with posterior probability >0.8. This suggests that genes belonging to these two groups might perform similar functions in some aspects [25,27,28].

These CAASs identified by DIVERGE2 were mapped onto the alignments of protein sequences (Additional file 4). We found that these CAASs were mainly located in four transmembrane regions (TMs 1, 2, 5, and 6), and only a few of them fell into the intra- or extra-cellular loops. We take the II/III pair as an example. There were 25 and 2 CAASs obtained for the two types of functional divergence respectively. Among the 27 predicted CAASs, there were 10 sites located in the loop regions, while 17 sites were in the TMs, with particular abundance in TM6. The similar cases were also found in other group pairs (Additional file 4). Given the observation that no possible CAAS was identified in TM4, it indicated that this region should be more conserved during evolution, and play important roles in maintaining the fundamental function of NIP aquapor-ins. In reality, conserved amino acids located in the helices 4 and 6 are essential to maintain the tetrameric structure of aquaporins [56].

We observed that among the predicted CAASs, two sites for the I/II pair, and one each for the II/III, and II/IV pairs corresponded to the first and/or third residues in the ar/R region (Additional file 4), suggesting that selective forces should have worked on the two residues of the ar/R filter, thereby leading to the diversity of substrate selectivity for different NIP proteins. In NIP I, the ar/R residues include W, V/I, A, and R, while these positions are occupied by A/ G/T/S/N, A/I/V, G/A, and R for the NIP II group. The two identified CAASs in the I/II pair correspond to the residues in the H2 and LE1 positions, which involve the changes of the size and hydrophilicity of the corresponding residues.

For example, the bulky Trp (W) in NIP I was replaced by the tiny A/S/T/N in NIP II; the hydrophilic G/S substituted for the hydrophibic A. These changes play a pivotal role in determining the substrate specificity. It has been reported that NIP I proteins can permeate small solutes such as water and glycerol, whereas NIP II proteins show no measurable water permeability but transport glycerol, for-mamide, as well as some larger uncharged solutes than NIP I protein. Experimental evidences showed that the substitution of a Trp (W) residue for Ala (A) at the H2 position of the ar/R tetrad of AtNIP6;1 results in enhanced water-transport activity, but the permeability of the mutant protein to urea, similar to the NIP I protein soybean nodulin 26, becomes more restrictive [31]. Using site-directed mutagenesis method, Wallace et al. [34] demonstrated that the residue in the H2 position should be a major determinant of glycerol selectivity.

In the II/III pair, the third residue (LE1) of the ar/R filter is predicted to be functional divergence related. This position is invariant Gly (G) in NIP III proteins, whereas the same position in NIP II contains several amino acids (glycine, alanine, or serine) with different properties, such as the nonpolar Ala (A). The size of the four residues (G, S, G, and R) in NIP III is smaller compared to other groups of NIPs, which form a larger constriction size (>6A) and allow for passage of much larger solutes [24]. However, NIP III proteins cannot be permeable to glycerol [57], although the molecular size of glycerol is smaller than that of silicic acid. This suggests that other unknown factors may also be involved in the substrate specificity [24].

To date, the monocot NIP3s (defined as NIP IV in this study) were included in the NIP II group. The present study however provided statistical evidence that monocot NIP3s should have largely diverged from other NIP II proteins in function, and fifteen CAASs including the H2 residue of the ar/R region were predicted to be responsible for the functional divergence (Table 1 and Additional file 4). In NIP IV, the residue in H2 is highly conserved to Ala (A). In contrast, several types of amino acids, such as alanine, asparagines, threonine, and serine, occupy the corresponding position in the NIP II group. However, whether this residue position is determinant for functional differentiation between NIP II and IV groups remains unknown.

Using more than 150 MIP proteins, Froger et al. [58] identified five positions (P1-P5) where the physicochemical properties of the corresponding amino acids are drastically different in aquaporins and glycerol permeases. Of the five positions, P4 and P5 correspond to two consecutive amino acids located in the sixth transmembrane segment (TM6). Lagree et al. [59] demonstrated that mutations of YW to PL at the two positions totally con-

verted the selectivity of the channel from water to glycerol. Thus, it seems that positions P4P5 play a crucial role in transport specificity. The five positions were designated in Additional file 4. We found that one position (P5) each for the I/II and II/III pairs was predicted to be highly functional divergence related, suggesting that the corresponding position should have acted as a determinant factor in subgroup-specific substrate transport.

Moreover, we detected some "hotspot" amino acid sites that are highly contributed to the functional divergence between different NIP subgroups. For example, there are three (I54, Q72, and A245) and four (P44, S85, Q125, and G242) such positions identified in OsNIP1;1 and OsNIP2;1 respectively (Additional file 4). However, because plant NIPs are comparatively less studied, the functional importance of the CAASs identified firstly needs to be further experimentally examined.

Expression analysis of NIP genes

The transcriptional patterns of rice NIP genes in nine tissues and forty cell types were investigated. It is observed that OsNIPs are unevenly expressed in the examined tissues, and exhibit a clearly tissue-specific expression pattern. OsNIP2;1, OsNIP3;1, and OsNIP2;2 are predominantly expressed in root (Figure 2A); their Specificity Measure (SPM) values are 0.994, 0.991, and 0.864 respectively. Similarly, relative to other tissues, OsNIP4;1 and OsNIP3;2 have a much stronger expression level in anther (SPM = 0.828) and suspension cell (SPM = 0.810) respectively, suggesting that they should play specific roles in the corresponding tissues. The similar cases were also found in Arabidopsis (Additional file 5). Consistent with the above postulation, Ma et al. [39] demonstrated that OsNIP2;1 was responsible for the uptake of silicic acid from soil. The poplar PtNIP1;1 was specifically expressed in the suspensor ligament of the embryo, and played important roles in the corresponding process [60]. Further, the microarray data for rice NIP genes in forty cell types were analyzed. It is obvious that OsNIPs represent a cell-type-specific expression pattern (Figure 2B), supporting the idea that NIP transport activities may be prevalent in a more defined set of cells in the plant [16]. This result was further validated by the analysis of Arabidopsis AtNIPs in the root cell-types (Additional file 6).

We further analyzed the expression of OsNIPs during the reproductive development, and found that the OsNIPs could be classified in two groups of expression pattern (Figure 2C). Five OsNIPs (3;3, 3;5, 1;4, 4;1, and 1;1) that are largely expressed in the inflorescence stage, represent the first group, while the other OsNIPs that appear to be root specific, form the second one. Therefore, the expression of NIPs would be developmental-related. Coinci-dently, it was found that PsNIP1;1 was expressed in the

developing seed coats of pea seeds in a developmentally specific manner [61].

The expression of AtNIPs is generally low under the treatment of constant white light, but the responses of AtNIPs are different (Figure 2D). For example, AtNIP3;1 exhibits a higher expression level under 36 hrs constant lighting. Moreover, we observed that most of the OsNIPs expression is down-regulated under a variety of stresses, such as drought, salt, and cold (Figure 2E). However, OsNIP4;1 is highly regulated by drought. These results indicate that NIPs are not exception in responding to a diverse array of stress related signals. Weig et al. [62] found that the expression of AtNIP1;1 was severely down-regulated by the stresses drought, salinity, as well as hormone abscisic acid. Arabidopsis AtNIP2;1 is an anaerobic-induced lactic acid transporter that may play a role in adaptation to lactic fermentation under anaerobic stress [21], which promotes us to further examine the expression profiles of OsNIP genes under the aerobic and anoxic conditions [63]. Interestingly, we observed that in rice coleoptile OsNIPs could be clearly divided into two groups: the first group includes four genes (OsNIPs 1;1, 2;1, 2;2, and 3;1), which are aerobic-dependent; the second group consists of the other seven OsNIPs, whose expression is regulated by anoxic stress (Figure 2F). Overall, the detailed study of NIP expression will refine our understanding of their physiological roles in plants.

Regulatory elements for plant NIP genes

Transcription factors bind to corresponding TFBSs (transcription factor binding sites) upstream from genes of interest, and the profiles of cis-acting elements may thus provide information for understanding the regulatory mechanism of gene expression. A computational tool PlantCARE [64] was adopted to identify putative TFBSs in the 1000 bp DNA sequence upstream of the translation initiation codon of NIP genes in rice, S. bicolor, grape and poplar.

Four types of cis-elements were found to be significantly abundant in the promoter region of plant NIP genes (Additional file 7). The first type of cis-element enriched in the promoter region is the light responsive elements, such as G-Box [65], GAG-motif [66], and Box 4 [67] etc. G-Box is the most abundant cis-element in rice. All but two (OsNIP3;3 and OsNIP4;1) have at least one copy of this element, whereas the Box 4 element appears to be more abundant in grape and poplar. Plant hormone responsive elements, such as ABRE [68], P-box [69], as well as the TCA-element [70], constitute the second class. It seems that ABRE is the most abundant hormone-related cis-element in rice, suggestive of the regulation of the expression of some OsNIPs by abscisic acid (ABA); whereas no such element has been detected in grape and

05 (Q CD CD O

to to £

o ^3 ^3

>s ri 3

g TO n

cr CT-

a- cr to

6 s re

TO > ^

60 > %

OQ C£ 3

3 3 3 &

n to a a ~ 8

O 3 TO 3 OQ

a » » 3 3


3 ^ ho 3 "

ft TO,

„ 3 3. OQ 3 3"

. . TO

3 ci a»

,—, to 3

^ 3 TO &

1—1 C/5 C/3

^—1 • (X> ¡D

to a ^

№ r+

№ cr Hi cr

3 O ^ rr

H-J (t>

£ & M r+ ,—, 3

^ 3 3 ^

3- 1=1 2 §

13 a" 3"

OQ ^ 3 O

Xk cr ^

<" T3 fl>

3 o T3

o Q. 3

« ■■ i5 n rt- !£.

-s n> p "> < 2 ft S-® 3-g 3

(t 3 "

rt (P S,

_ n " 3

S-l" ? *

§1 i/i a. rt o

fl> org — fl)

3" fl> rt W

3- 3 2 f*

- OQ P


7-day-old Seedling 1 7-day-old Seedling 2 7-day-old Seedling 3 Drought stress 1 Drought stress 2 Drought stress 3 Salt stress 1 Salt stress 2 Salt stress 3 Cold stress 1 Cold stress 2 Cold stress 3

_l M _l ftl A (.1 (il

■ N) GJ N) 01 GJ ^ ■

r^TH rrh


to N1 Kl rj -i -i ii. CT1 fj

Root Root Root

Mature leaf Mature leaf Mature leaf

Y Leaf

Y Leaf


Young inflorescence Young inflorescence Young inflorescence Inflorescence (P2 Inflorescence (P2 Inflorescence (P2 Inflorescence (P3 Inflorescence (P3 Inflorescence (P3 Inflorescence (P4 Inflorescence (P4 Inflorescence (P4 Inflorescence (P5 Inflorescence (P5 Inflorescence (P5 Inflorescence (P6 Inflorescence (P6 Inflorescence (P6 Seed (S1 Seed (S1 Seed (S1 Seed (S2 Seed (S2 Seed (S2 Seed (S3 Seed (S3 Seed(S3 Seed(S4 Seed(S4 Seed(S4 Seed(S5 Seed(S5 Seed(S5

(P1 (P1 (P1


cil M -i ùj -i -i Ni N1 "i ùj






Ovary 2







Scutellum (Ohr) Scutellum (12hr) Scutellum (24hr) Coleoptlle (Ohr) Coleoptlle (12hr) Coleoptlle (24hr) Plumule (Ohr) Plumule (12hr) Plumule (24hr) Eplblast (Ohr) Eplblast (12hr) Epl blast (24hr) Radicle (Ohr) Radicle (12hr) Radicle (24hr) Axillary primordium Axillary meristem Apical meristem P1 P2 P3

Seedling blade bulllform Seedling blade stomata Seedling blade long cell Seedling blade mesophyll Seedling blade bundle sheath Seedling blade vein Lateral root cap Root tl p cortex Root tip vascular bundle Root tl p metaxy I em Elongation epidermis Elongation cortex Elongation endodermis Elong. vascular bundle Elongation metaxylem Maturation epidermis Maturation cortex Maturation endodermis Matur vascular bundle Whole root Whole I eat (fresh)

ro o o

S?cî?cî?cï?cï?cï?c7?c7?<ï? 2ZZZZZZZZ


¡=t TD

CD Q. O fD

two distinct expression patterns (Figure 2F). We speculated therefore that the anaerobic regulation of OsNIPs expression should be tissue or developmental-stage dependent. The drought responsive element MBS [72] is also enriched in the promoter. With a few exceptions, OsNIPs contain more than two copies of this element (Additional file 7). Circadian that is involved in circadian control [73] is the fourth type of cis-element found abundant in the promoters of monocot NIP genes. PlantCARE [64] identified one circadian in the rice silicic acid transporter OsNIP2;1 [39], which may be responsible for its distinct diurnal expression pattern [74]. The presence of a diverse of cis-elements in the upstream regions of NIPs indicates that plant NIPs may function in a wider range of ways.

In addition, NIP membership-specific as-elements have been observed. For example, in rice, OsNIP2;1 and OsNIP3;2 each possesses a putative LTR motif that is response to low temperature [75,76]. AC-I, a cis-element conferring enhanced xylem expression [77,78], is specifically present in OsNIP1;1. Furthermore, orthologous NIPs in monocot and dicot plants have different cis-element organization as well (Additional file 7). These results suggest that plant NIPs should have evolved some specific regulatory elements, and thereby leading to the differentiation of expression patterns.


Plant NIP subfamily is more ancient, and their diversification can be placed at the time before the emergence of the moss P. patens. As many typical gene families, NIPs have experienced strong positive selection during evolution. Consistently, the amino acid level analysis suggests that functional divergence has occurred between plant NIP proteins, and identified the critical amino acid sites involved in this divergence for further investigation. Rice and Arabidopsis NIPs exhibit distinct expression pattern. The survey of upstream elements reveals four major classes of cis-elements in the promoter region of NIPs and their distinct organization pattern is interpreted to reflect their varying participation in gene expression regulation. The rapid proliferation and functional diversification of plant NIPs is argued to have partially attributed to the need for plants to better adapt to external different environments. These findings provide new insights into understanding the evolutionary mechanisms of NIP proteins and their functional diversification.

Methods Sequence data

O. sativa, A. thaliana, Z. Mays, and P. patens NIP sequences were downloaded from the GenBank and JGI databases according to the published literatures [3,7,49,79]. The NIPs were used as query to search against the V. vinifera, P. trichocarpa, S. bicolor, G. max, C. sativus, and P. taeda

genomes using the BLASTP and TBLASTN programs respectively. Other sequences including the Cicer arieti-num CaNIP (AN: CAG34223), and Cucurbita pepo CpNIP (AN: CAD67694) were collected from literature [32]. Programs InterProScan [80] and ConPred II [81] were employed to detect the conserved MIP domain and predict the putative transmembrane regions respectively.

Multiple sequence alignment and phylogenetic tree reconstruction

Plant NIP protein sequences were aligned using the program L-INS-i implemented in MAFFT v6.6 [82], with the parameters: Scoring matrix for amino acid sequences, BLOSUM62; Gap opening penalty, 2.0; and Gap extension penalty, 0.2. The resulting protein alignment was subsequently employed to generate the codon-alignment of corresponding coding DNA sequences using a custom PERL script. Maximum likelihood (ML) phylogenies were reconstructed with PHYML v2.4 [83]. The programs PROTTEST [84] and ModelGenerator [85] were utilized to determine the best model for each ML analysis. Here the JTT+I+G model for the protein alignment and the HKY+I+G+F model for the codon alignment were determined respectively. The reliability of interior branches was assessed with 500 bootstrap resamplings. Phyloge-netic trees were displayed using MEGA v4.0 [86].

Test of positive selection

The CODEML program implemented in the PAML v4.0 software package [51] was utilized to test the hypothesis of positive selection in the NIP subfamily during evolution. To test for heterogeneous selective pressure among lineages [50], models of variable ra ratios among lineages were fitted by ML to the NIP sequence alignment. The ratio of nonsynonymous-to-synonymous for each branch under two models (one-ratio and free-ratio for branches) was calculated, and the two models were compared using the LRT test to see whether the ra ratios are different among lineages; that is, positive selection is indicated if the free-ratio model that allows for selection is significantly better than the one-ratio model (no selection) in the LRT analysis.

Estimation of functional divergence

The software DIVERGE2 [55] was used to detect functional divergence between members of the plant NIP protein subfamily, where 4 gene clusters of interests were selected. The coefficients of Type-I and Type-II functional divergence 9I and 9II between any two NIP groups were calculated. If 9I or 9II is significantly greater than 0, it means site-specific altered selective constraints or a radical shift of amino acid physiochemical property after gene duplication [54,55]. Moreover, a site-specific posterior analysis was used to predict amino acid residues that were crucial for functional divergence.

Investigation of transcription patterns

Gene expression microarray datasets (GSE7951, GSE13161, GSE6893, GSE6908, and GSE6901 for rice; GSE680, GSE7641, and GSE8365 for Arabidopsis) were downloaded from the GEO database in NCBI. The micro-array data of rice include the analysis of gene expression profiles in nine tissues [87] and forty cell types; during reproductive development; seven-day-old seedlings under drought, salt, and cold stress treatments [88]; and rice coleoptile under the aerobic and anoxic conditions [63]. In Arabidopsis, the transcript profiles of NIP genes in root cell-types after treatment with salt [89], under constant white light treatment [90], and during the whole plant life cycle, were investigated as well. Program dChip 2008 (5/ 8/08) [91] was used to perform the cluster analysis and display the expression patterns of rice and Arabidopsis NIP genes using microarray data as input. The GEPS software [92] was employed to quantitatively analyze the expression pattern of NIP genes. Specificity Measure (SPM) was used to define the tissue-specific expression pattern of a gene, which may be useful for further understanding its physiological behaviors [92].

Analysis of cis-acting regulatory elements

1000 bp of nucleotide sequences upstream of the translation initiation codon for each NIP gene in four species (rice, S. bicolor, grape, and poplar) was extracted, which were used for the TFBSs analysis. At present, no full-length cDNA sequences for grape were available. In order to facilitate comparison between species, the sequences upstream of translation initiation codon rather than transcription start site were used to screen possible cis-acting regulatory elements. The software PlantCARE [64] was utilized to determine putative plant-specific TFBSs in a given DNA sequence. To avoid biases in analysis, only TFBSs whose matrix score is not less than 6 were considered further.

Authors' contributions

QL and ZZ conceived and designed the experiments. QL, HW, and YF performed the experiments and analyzed the data. ZHZ identified and annotated the exon/intron structures of the NIP genes in cucumber. JW facilitated the usage of the dChip program and performed the statistical analysis. QL wrote the paper. All authors read and approved the final manuscript.

Additional material

Additional file 1

List of NIP genes in plants. For each NIP gene, the information about the accession number, chromosomal localization, gene length, protein length, gene structure, the NPA motif, and the ar/R filter sequence were listed.

Click here for file


Additional file 2

The multiple alignments of plant NIP protein sequences. These sequences were aligned using the program L-INS-i implemented in MAFFT v6.6. The alignments were shaded in the "Quantify Mode", and the residues were displayed in the "Difference Mode" with the "Diff/Con-sensus Line" style. Dots and "-" indicate similar residues and gaps on the alignment respectively. Click here for file


Additional file 3

Phylogenetic tree reconstructed using plant NIP nucleotide sequences. The number beside the branches represents bootstrap values >300 based on 500 resamplings. The scale bar shows total nucleotide distance. The NIP homologue in the green alga Ostreococcus lucimarinus (defined as galgaNIP) is used as outgroup sequence to root the tree. Branches with rates of numbers of nonsynonymous and synonymous substitutions >1, are indicated by red thick lines. Click here for file


Additional file 4

Functional divergence significantly related amino acid site candidates. A site-specific profile based on the posterior probability (Qk) was used to identify critical amino acid sites that were responsible for functional divergence between NIP subfamily members. According to the definition, large Qk indicates a high possibility that the functional constraint (or, the evolutionary rate) or the radical change of amino acid property of a site is different between two clusters. Dots indicate conserved residues with the first protein, such as OsNIP2-1. Amino acids (AAs) with Qk values 0.9>Qk>0.8 and Qk>0.9 are shaded with green and red color respectively. The predicted AAs for Type-II functional divergence are shown in purple. The six predicted transmembrane regions (TMs) are marked above the sequences. Arrows indicate the first and/or the third residue(s) in the ar/R region. The five positions identified by Froger et al. [58] are indicated by stars. (A) I/II; (B) I/III; (C) I/IV; (D) II/III; (E) II/IV; (F) III/IV. Click here for file


Additional file 5

The expression patterns of Arabidopsis AtNIPs during the plant life cycle.

Click here for file


Additional file 6

The diversity of expression profiles of Arabidopsis AtNIP genes in root cell-types after treatment with salt.

Click here for file


Additional file 7

Analysis of cis-acting elements in the 1000 bp sequence upstream of the translation initiation codon in plant NIP genes.

Click here for file



We thank the anonymous reviewers for their valuable and constructive suggestions for this manuscript. This work was supported by the Foundation of Zhejiang Educational Committee, China (No. Y200805229), the intramural fund from Zhejiang Forestry University (to Qingpo Liu), and grants from National Basic Research Program of China (973 program; No. 2007CBI09305) to YF, and the Special Fund for Grade B Innovative Research Team from Zhejiang Forestry University to ZZ.


1. Kruse E, Uehlein N, Kaldenhoff R: The aquaporins. Genome Biol

2006, 7:206.

2. Gorin MB, Yancey SB, Cline J, Revel JP, Horwitz J: The major intrinsic protein (MIP) of the bovine lens fiber membrane: characterization and structure based on cDNA cloning. Cell 1984, 39:49-59.

3. Johanson U, Karlsson M, Johansson I, Gustavsson S, Sjovall S, Fraysse L, Weig AR, Kjellbom P: The complete set of genes encoding major intrinsic proteins in Arabidopsis provides a framework for a new nomenclature for major intrinsic proteins in plants. Plant Physiol 2001, 126:1358-1369.

4. Bansal A, Sankararamakrishnan R: Homology modeling of major intrinsic proteins in rice, maize and Arabidopsis: comparative analysis of transmembrane helix association and aromatic/ arginine selectivity filters. BMC Struct Biol 2007, 7:27.

5. Portincasa P, Palasciano G, Svelto M, Calamita G: Aquaporins in the hepatobiliary tract. Which, where and what they do in health and disease. Eur J Clin Invest 2008, 38:1-10.

6. Borstlap AC: Early diversification of plant aquaporins. Trends Plant Sci 2002, 7:529-530.

7. Danielson JAH, Johanson U: Unexpected complexity of the Aquaporin gene family in the moss Physcomitrella patens . BMC Plant Biol 2008, 8:45.

8. Gustavsson S, Lebrun AS, Norden K, Chaumont F, Johanson U: A novel plant major intrinsic protein in Physcomitrella patens most similar to bacterial glycerol channels. Plant Physiol 2005, 139:287-295.

9. Zardoya R: Phylogeny and evolution of the major intrinsic protein family. Biol Cell 2005, 97:397-414.

10. Tyerman SD, Niemietz CM, Bramley H: Plant aquaporins: multifunctional water and solute channels with expanding roles. Plant Cell Environ 2002, 25:173-194.

11. Chaumont F, Moshelion M, Daniels MJ: Regulation of plant aquaporin activity. Biol Cell 2005, 97:749-764.

12. Kaldenhoff R, Fischer M: Functional aquaporin diversity in plants. Biochim Biophys Acta 2006, 1758: I 134-1141.

13. Jang JY, Lee SH, Rhee JY, Chung GC, Ahn SJ, Kang H: Transgenic Arabidopsis and tobacco plants overexpressing an aquaporin respond differently to various abiotic stresses. Plant Mol Biol

2007, 64:62I -632.

14. Zhang Y, Wang Z, Chai T, Wen Z, Zhang H: Indian mustard aquaporin improves drought and heavy-metal resistance in tobacco. Mol Biotechnol 2008, 40:280-292.

15. Fortin MG, Morrison NA, Verma DP: Nodulin-26, a peribacteroid membrane nodulin is expressed independently of the development of the peribacteroid compartment. Nucleic Acids Res I987, l5:8I3-824.

16. Wallace IS, Choi WG, Roberts DM: The structure, function and regulation of the nodulin26-like intrinsic protein family of plant aquaglyceroporins. Biochim Biophys Acta 2006, 1758: I I65-II75.

17. Zardoya R, Ding X, Kitagawa Y, Chrispeels MJ: Origin of plant glyc-erol transporters by horizontal gene transfer and functional recruitment. Proc Natl Acad Sci USA 2002, 99: I4893-I4896.

18. Zhao CX, Shao HB, Chu LY: Aquaporin structure-function relationships: water flow through plant living cells. Colloids Surf B Biointerfaces 2008, 62: I63-172.

19. Dean RM, Rivers RL, Zeidel ML, Roberts DM: Purification and functional reconstitution of soybean nodulin 26, an aquaporin with water and glycerol transport properties. Biochemistry I999, 38:347-353.

20. Niemietz CM, Tyerman SD: Channel-mediated permeation of ammonia gas through the peribacteroid membrane of soybean nodules. FEBS Lett 2000, 465: II0-II4.

21. Choi WG, Roberts DM: Arabidopsis NIP2;l, a major intrinsic protein transporter of lactic acid induced by anoxic stress. J Biol Chem 2007, 282:24209-24218.

22. Takano J, Wada M, Ludewig U, Schaaf G, von Wiren N, Fujiwara T: The Arabidopsis major intrinsic protein NIP5;l is essential for efficient boron uptake and plant development under boron limitation. Plant Cell 2006, 18: I498-I509.

23. Ma JF, Yamaji N: Functions and transport of silicon in plants. Cell Mol Life Sci 2008, 65:3049-3057.

24. Mitani N, Yamaji N, Ma JF: Characterization of substrate specificity of a rice silicon transporter, Lsil. Pflugers Arch - Eur J Physiol 2008, 456:679-686.

25. Bienert GP, Thorsen M, Schussler MD, Nilsson HR, Wagner A, Tamas MJ, Jahn TP: A subgroup of plant aquaporins facilitate the bidirectional diffusion of As(OH)3 and Sb(OH)3 across membranes. BMC Biol 2008, 6:26.

26. Isayenkov SV, Maathuis FJM: The Arabidopsis thaliana aquaglyc-eroporin AtNIP7;l is a pathway for arsenite uptake. FEBS Lett 2008, 582: I625-1628.

27. Ma JF, Yamaji N, Mitani N, Xu XY, Su YH, McGrath SP, Zhao FJ: Transporters of arsenite in rice and their role in arsenic accumulation in rice grain. Proc Natl Acad Sci USA 2008, l05:993I-9935.

28. Kamiya T, Tanaka M, Mitani N, Ma JF, Maeshima M, Fujiwara T: NIPl;l, an aquaporin homolog, determines the arsenite sensitivity of Arabidopsis thaliana. J Biol Chem 2009, 284:2II4-2120.

29. Wallace IS, Roberts DM: Homology modeling of representative subfamilies of Arabidopsis major intrinsic proteins. Classification based on the aromatic/arginine selectivity filter. Plant Physiol 2004, 135:1059-I068.

30. Forrest KL, Bhave M: Major intrinsic proteins (MIPs) in plants: a complex gene family with major impacts on plant pheno-type. Funct Integr Genomics 2007, 7:263-289.

31. Wallace IS, Roberts DM: Distinct transport selectivity of two structural subclasses of the nodulin-like intrinsic protein family of plant aquaglyceroporin channels. Biochemistry 2005, 44: I6826-I6834.

32. Rouge P, Barre A: A molecular modeling approach defines a new group of nodulin26-like aquaporins in plants. Biochem Biophys Res Commun 2008, 367:60-66.

33. Wu B, Beitz E: Aquaporins with selectivity for unconventional permeants. Cell Mol Life Sci 2007, 64:2413-2421.

34. Wallace IS, Wills DM, Guenther JF, Roberts DM: Functional selectivity for glycerol of the nodulin 26 subfamily of plant membrane intrinsic proteins. FEBS Lett 2002, 523: I09-I I2.

35. Agre P, Kozono D: Aquaporin water channels: molecular mechanisms for human disease. FEBS Lett 2003, 555:72-78.

36. Sui H, Han BG, Lee JK, Walian P, Jap BK: Structural basis of water-specific transport through the AQPl water channel. Nature 200I, 4l4:872-878.

37. Weaver CD, Crombie B, Stacey G, Roberts DM: Calcium-dependent phosphorylation of symbiosome membrane proteins from nitrogen-fixing soybean nodules: evidence for phospho-rylation of nodulin-26. Plant Physiol I99I, 95:222-227.

38. Guenther JF, Roberts DM: Water-selective and multifunctional aquaporins from Lotus japonicus nodules. Planta 2000, 2 l0:74I-748.

39. Ma JF, Tamai K, Yamaji N, Mitani N, Konishi S, Katsuhara M, Ishiguro M, Murata Y, Yano M: A silicon transporter in rice. Nature 2006, 440:688-69I.

40. Yamaji N, Mitatni N, Ma JF: A transporter regulating silicon distribution in rice shoots. Plant Cell 2008, 20: I38I-I389.

41. Chiba Y, Mitani N, Yamaji N, Ma JF: HvLsil is a silicon influx transporter in barley. Plant J 2009, 57:810-818.

42. Mizutani M, Watanabe S, Nakagawa T, Maeshima M: Aquaporin NIP2;l is mainly localized to the ER membrane and shows root-specific accumulation in Arabidopsis thaliana. Plant Cell Physiol 2006, 47: I420-1426.

43. Guenther JF, Chanmanivone N, Galetovic MP, Wallace IS, Cobb JA, Roberts DM: Phosphorylation of soybean nodulin 26 on serine 262 enhances water permeability and is regulated develop-mentally and by osmotic signals. Plant Cell 2003, l5:98I-99I.

44. Johansson I, Karlsson M, Johanson U, Larsson C, Kjellbom P: The role of aquaporins in cellular and whole plant water balance. Biochim Biophys Acta 2000, l465:324-342.

45. Raven JA: Costs and benefits of low osmolarity in cells of freshwater algae. Funct Ecol I995, 9:70I-707.

46. Blanc G, Hokamp K, Wolfe KH: A recent polyploidy superimposed on older large-scale duplications in the Arabidopsis genome. Genome Res 2003, l3: I37-I44.

47. Guyot R, Keller B: Ancestral genome duplication in rice. Genome 2004, 47:610-614.

48. Tanaka M, Wallace IS, Takano J, Roberts DM, Fujiwara T: NIP6;I is a boric acid channel for preferential transport of boron to growing shoot tissues in Arabidopsis. Plant Cell 2008, 20:2860-2875.

49. Sakurai J, Ishikawa F, Yamaguchi T, Uemura M, Maeshima M: Identification of 33 rice aquaporin genes and analysis of their expression and function. Plant Cell Physiol 2005, 46:1568-1577.

50. Yang Z: Likelihood ratio test for detecting positive selection and application to primate lysozyme evolution. Mol Biol Evol 1998, 15:568-573.

51. Yang Z: PAML4: a program package for phylogenetic analysis by maximum likelihood. Mol Biol Evol 2007, 24:1586-1591.

52. Zhu S: Positive selection targeting the cathelin-like domain of the antimicrobial cathelicidin family. Cell Mol Life Sci 2008, 65:1285-1294.

53. Lynch M, Conery JS: The evolutionary fate and consequences of duplicate genes. Science 2000, 290:1151-1155.

54. Gu X: Statistical methods for testing functional divergence after gene duplication. Mol Biol Evol 1999, 16:1664-1674.

55. Gu X: A simple statistical method for estimating type-II (cluster-specific) functional divergence of protein sequences. Mol Biol Evol 2006, 23:1937-1945.

56. Walz T, Hirai T, Murata K, Heymann JB, Mitsuoka K, Fujiyoshi Y, Smith BL, Agre P, Engel A: The three-dimensional structure of aquaporin-I. Nature 1997, 387:624-627.

57. Klebl F, Wolf M, Sauer N: A defect in the yeast plasma membrane urea transporter Dur3p is complemented by CpNIP1, a Nod26-like protein from zucchini (Cucurbita pepo L.), and by Arabidopsis thaliana 5-TIP or y-TIP. FEBS Lett 2003, 547:69-74.

58. Froger A, Tallur B, Thomas D, Delamarche C: Prediction of functional residues in water channels and related proteins. Protein Sci 1998, 7:1458-1468.

59. Lagree V, Froger A, Deschamps S, Hubert JF, Delamarche C, Bonnec G, Thomas D, Gouranton J, Pellerin I: Switch from an aquaporin to a glycerol channel by two amino acids substitution. J Biol Chem 1999, 274:6817-6819.

60. Ciavatta VT, Egertsdotter U, Clapham D, von Arnold S, Cairney J: A promoter from the loblolly pine PtNIPI;I gene directs expression in an early-embryogenesis and suspensor-specific fashion. Planta 2002, 215:694-698.

61. Schuurmans JA, van Dongen JT, Rutjens BP, Boonman A, Pieterse CM, Borstlap AC: Members of the aquaporin family in the developing pea seed coat include representatives of the PIP, TIP, and NIP subfamilies. Plant Mol Biol 2003, 53:633-645.

62. Weig A, Deswarte C, Chrispeels MJ: The major intrinsic protein family of Arabidopsis has 23 members that form three distinct groups with functional aquaporins in each group. Plant Physiol 1997, 114:1347-1357.

63. Lasanthi-Kudahettige R, Magneschi L, Loreti E, Gonzali S, Licausi F, Novi G, Beretta O, Vitulli F, Alpi A, Perata P: Transcript profiling of the anoxic rice coleoptile. Plant Physiol 2007, 144:218-231.

64. Lescot M, Dehais P, Thijs G, Marchal K, Moreau Y, Peer Y Van de, 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.

65. Sommer H, Saedler H: Structure of the chalcone synthase gene of Antirrhinum majus. Mol Gen Genet 1986, 202:429-434.

66. Werneke JM, Chatfield JM, Ogren WL: Alternative mRNA splicing generates the two ribulosebisphosphate carboxylase/ oxygenase activase polypeptides in spinach and Arabidopsis. Plant Cell 1989, 1:815-825.

67. Lois R, Dietrich A, Hahlbrock K, Schulz W: A phenylalanine ammonia-lyase gene from parsley: structure, regulation and identification of elicitor and light responsive cis-acting elements. EMBO J 1989, 8:1641-1648.

68. Simpson SD, Nakashima K, Narusaka Y, Seki M, Shinozaki K, Yamaguchi-Shinozaki K: Two different novel cis-acting elements of erdI, a clpA homologous Arabidopsis gene function in induction by dehydration stress and dark-induced senescence. Plant J 2003, 33:259-270.

69. Kim JK, Cao J, Wu R: Regulation and interaction of multiple protein factors with the proximal promoter regions of a rice high pI alpha-amylase gene. Mol Gen Genet 1992, 232:383-393.

70. Pastuglia M, Roby D, Dumas C, Cock JM: Rapid induction by wounding and bacterial infection of an S gene family receptor-like kinase in Brassica oleracea. Plant Cell. 1997, 9(1):49-60.

71. Manjunath S, Sachs MM: Molecular characterization and promoter analysis of the maize cytosolic glyceraldehyde 3-phos-phate dehydrogenase gene family and its expression during anoxia. Plant Mol Biol 1997, 33:97-1 12.

72. Yamaguchi-Shinozaki K, Shinozaki K: Arabidopsis DNA encoding two desiccation-responsive rd29 genes. Plant Physiol 1993, 101:11 19-1 120.

73. Pichersky E, Bernatzky R, Tanksley SD, Breidenbach RB, Kausch AP, Cashmore AR: Molecular characterization and genetic mapping of two clusters of genes encoding chlorophyll a/b-biding proteins in Lycopersicon esculentum (tomato). Gene 1985, 40:247-258.

74. Yamaji N, Ma JF: Spatial distribution and temporal variation of the rice silicon transporter Lsi1. Plant Physiol 2007, 143:1306-1313.

75. Baker SS, Wilhelm KS, Thomashow MF: The 5'-region of Arabidopsis thaliana core 15a has cis-acting elements that confer cold, drought- and ABA-regulated gene expression. Plant Mol Biol 1994, 24:701-713.

76. White AJ, Dunn MA, Brown K, Hughes MA: Comparative analysis of genomic sequence and expression of a lipid transfer protein gene family in winter barley. J Exp Bot 1994, 45:1885-1892.

77. Hatton D, Sablowski R, Yung MH, Smith C, Schuch W, Bevan M: Two classes of cis sequences contribute to tissue-specific expression of a PAL2 promoter in transgenic tobacco. Plant J 1995, 7:859-876.

78. Patzlaff A, Newman LJ, Dubos C, Whetten RW, Smith C, McInnis S, Bevan MW, Sederoff RR, Campbell MM: Characterization of PtMYBI, an R2R3-MYB from pine xylem. Plant Mol Biol 2003, 53:597-608.

79. Chaumont F, Barrieu F, Wojcik E, Chrispeels MJ, Jung R: Aquaporins constitute a large and highly divergent protein family in maize. Plant Physiol 2001, 125:1206-1215.

80. Quevillon E, Silventoinen V, Pillai S, Harte N, Mulder N, Apweiler R, Lopez R: InterProScan: protein domains identifier. Nucleic Acids Res 2005, 33:WI 16-120.

81. Arai M, Mitsuke H, Ikeda M, Xia JX, Kikuchi T, Satake M, Shimizu T: ConPred II: a consensus prediction method for obtaining transmembrane topology models with high reliability. Nucleic Acids Res 2004, 32:W390-393.

82. Katoh K, Kuma K, Miyata T, Toh H: Improvement in the accuracy of multiple sequence alignment program MAFFT. Genome Inform 2005, 16:22-33.

83. Guindon S, Gascuel O: A simple, fast, and accurate algorithm to estimate large phylogenies by maximum likelihood. Syst Biol 2003, 52:696-704.

84. Abascal F, Zardoya R, Posada D: ProtTest: selection of best-fit models of protein evolution. Bioinformatics 2005, 21:2104-2105.

85. Keane TM, Creevey CJ, Pentony MM, Naughton TJ, Mclnerney JO: Assessment of methods for amino acid matrix selection and their use on empirical data shows that ad hoc assumptions for choice of matrix are not justified. BMC Evol Biol 2006, 6:29.

86. Tamura K, Dudley J, Nei M, Kumar S: MEGA4: molecular evolutionary genetics analysis (MEGA) software version 4.0. Mol Biol Evol 2007, 24:1596-1599.

87. Li M, Xu W, Yang W, Kong Z, Xue Y: Genome-wide gene expression profiling reveals conserved and novel molecular functions of the stigma in rice. Plant Physiol 2007, 144:1797-1812.

88. Jain M, Nijhawan A, Arora R, Agarwal P, Ray S, Sharma P, Kapoor S, Tyagi AK, Khurana JP: Genome-wide analysis, classification, emporal and spatial gene expression during panicle and seed development, and regulation by light and abiotic stress. Plant Physiol 2007, 143:1467-1483.

89. Dinneny JR, Long TA, Wang JY, Jung JW, Mace D, Pointer S, Barron C, Brady SM, Schiefelbein J, Benfey PN: Cell identity mediates the response of Arabidopsis roots to abiotic stress. Science 2008, 320:942-945.

90. Covington MF, Harmer SL: The circadian clock regulates auxin signaling and responses in Arabidopsis. PLOS Biol 2007, 5:e222.

91. Li C: Automaing dChip: toward reproducible sharing of microarray data analysis. BMC Bioinformatics. 2008, 9:231.

92. Wang YP, Liang L, Han BC, Quan Y, Wang X, Tao T, Ji ZL: GEPS: the gene expression pattern scanner. Nucleic Acids Res 2006, 34:W492-497.