Skip to main content
  • Research article
  • Open access
  • Published:

New subfamilies of major intrinsic proteins in fungi suggest novel transport properties in fungal channels: implications for the host-fungal interactions



Aquaporins (AQPs) and aquaglyceroporins (AQGPs) belong to the superfamily of Major Intrinsic Proteins (MIPs) and are involved in the transport of water and neutral solutes across the membranes. MIP channels play significant role in plant-fungi symbiotic relationship and are believed to be important in host-pathogen interactions in human fungal diseases. In plants, at least five major MIP subfamilies have been identified. Fungal MIP subfamilies include orthodox aquaporins and five subgroups within aquaglyceroporins. XIP subfamily is common to both plants and fungi. In this study, we have investigated the extent of diversity in fungal MIPs and explored further evolutionary relationships with the plant MIP counterparts.


We have extensively analyzed the available fungal genomes and examined nearly 400 fungal MIPs. Phylogenetic analysis and homology modeling exhibit the existence of a new MIP cluster distinct from any of the known fungal MIP subfamilies. All members of this cluster are found in microsporidia which are unicellular fungal parasites. Members of this family are small in size, charged and have hydrophobic residues in the aromatic/arginine selectivity filter and these features are shared by small and basic intrinsic proteins (SIPs), one of the plant MIP subfamilies. We have also found two new subfamilies (δ and γ2) within the AQGP group. Fungal AQGPs are the most diverse and possess the largest number of subgroups. We have also identified distinguishing features in loops E and D in the newly identified subfamilies indicating their possible role in channel transport and gating.


Fungal SIP-like MIP family is distinct from any of the known fungal MIP families including orthodox aquaporins and aquaglyceroporins. After XIPs, this is the second MIP subfamily from fungi that may have possible evolutionary link with a plant MIP subfamily. AQGPs in fungi are more diverse and possess the largest number of subgroups. The aromatic/arginine selectivity filter of SIP-like fungal MIPs and the δ AQGPs are unique, hydrophobic in nature and are likely to transport novel hydrophobic solutes. They can be attractive targets for developing anti-fungal drugs. The evolutionary pattern shared with their plant counterparts indicates possible involvement of new fungal MIPs in plant-fungi symbiosis and host-pathogen interactions.


The superfamily of Major Intrinsic Proteins (MIPs) contains channel proteins that transport water and other neutral solutes [1]–[4]. These integral membrane proteins are found from bacteria to humans and are abundantly present in plants [5],[6]. Aquaporin (AQP) and aquaglyceroporin (AQGP) are the prototype members of the MIP superfamily [7],[8]. While in mammals three major subfamilies are detected [9], plant MIPs are found to have at least five subfamilies [10]–[12]. In humans, MIP members play significant role in kidney nephron, epithelial fluid secretion, maintaining brain water balance, cell migration, skin hydration, adipocyte metabolism and neuroexcitation [9],[13]–[15]. They are implicated in various human diseases such as glaucoma, epilepsy, cancer, obesity, nephrogenic diabetes insipidus and neuromyelitis optica [15]–[22]. In plants, multiple roles have been recognized in terms of plant development, growth and physiology [2],[23],[24]. They have been shown to be important in stress tolerance in plants and they play vital role in their interactions with soil microorganisms, fungi and pathogens [25]–[29]. In addition to water and glycerol, expression studies in Xenopus oocytes showed that MIPs transport diverse neutral solutes and gases such as ammonia and CO2[1],[3],[25]. MIP members also serve as an important component in host-parasite interactions and several protozoan parasite aquaporins have been identified including those from Plasmodium, Trypanosoma and Leishmania species [30]–[34]. Functional studies suggest that protozoan aquaporins are likely to transport both water and glycerol efficiently [31]. These parasite channel proteins could either serve as potential drug targets or vehicles for transporting cytotoxic compounds [30],[31],[35]. We have recently identified more than 1000 MIP genes from more than 340 organisms and details about the genes, protein products and structural models are available in MIPModDB database [6].

Symbiotic relationship between plants and fungi occurs through mycorrhiza and these interactions help in translocating soil nutrients to the host plants [36],[37]. This mutualistic association also aids in the transfer of organic carbon from plants to fungal partners [38]. In the mycorrhized plants, movement of water and other nutrients are greatly influenced by the expression of both plant and fungal aquaporins and this plays a major role in the drought resistance of plants [39],[40]. The expression of plant aquaporins in mycorrhized and non-mycorrhized forms have been investigated in several plants under different stress conditions [27],[28],[41]–[43]. On the contrary, very few studies have been carried out on the role of fungal aquaporins in mycorrhizal symbiosis [44],[45]. In addition to their role in root water transport in plants through symplastic pathway, fungal aquaporins present in pathogenic fungal species may act as attractive targets for antifungal drugs [16],[31],[46]. With increasing drug resistance reported in human fungal pathogens [47],[48], it has become even more important to identify new drug targets in fungal organisms.

Fungal aquaporins have been identified in individual species such as Laccaria bicolor (an ectomycorrhizal fungus) [44]Glomus intraradices (an arbuscular mycorrhizal fungus) [39] and Encephalitozoon cuniculi (a microsporidia pathogenic to humans) [49]. Their transport properties, gene expression profiles and role in mycorrhizia or pathogenesis have been investigated. Available literature shows that only a few fungal MIPs have been studied in detail. Fungal species are diverse and form a large group of eukaryotic organisms and hence it is anticipated that the MIP channels present in fungi also could be diverse with different transport specificities and unique regulatory mechanisms and/or physiological functions. In an earlier study, Pettersson et al. [50] analyzed genomes of yeast and filamentous fungi. Their phylogenetic analysis with a small number of 55 MIPs revealed three different subgroups within fungal AQGPs in addition to members belonging to orthodox AQP family. A new subfamily of fungal aquaporins, called X-intrinsic proteins (XIPs), has been discovered in our laboratory which is found to be common between fungi and plants [12]. More recently, Dietz et al. [44] and Zwiazek and coworkers [51] have identified even a larger number of 135 and 229 fungal MIPs respectively from database searches. Both their phylogenetic studies indicated at least four clusters including the orthodox AQPs, XIPs and two to four distinct subgroups of AQGPs. With the improved technologies and falling costs of sequencing, fungal genomes are rapidly sequenced and they are available in public databases. In this study, we have systematically searched completely and partially sequenced fungal genomes and identified a large number of fungal MIPs. With a comprehensive analysis of these MIP sequences using different bioinformatics tools, we have identified a new group of fungal MIPs which forms a distinct cluster separately from the orthodox AQPs and AQGPs. Members of this new family share several features with plant SIP (small and basic intrinsic proteins) subfamily. Our analysis also revealed two new subfamilies of aquaglyceroporins which are separate from the other subgroups of AQGPs. Distinct features of these three new families of fungal MIPs, their possible evolutionary relationship with plants and lower organisms and the nature of solutes that are likely to be transported across these channels are presented in detail.


Identification of fungal MIPs and their taxonomic distributions

Ten MIP sequences representing different subfamilies and organism groups were used as query sequences in the PSI-BLAST and BLAST searches to search the non-redundant database in NCBI ( and UniProt ( database [52] respectively. Although the MIP sequences are diverse, they all adopt a characteristic hour-glass helical fold consisting of six transmembrane helical segments (TM1 to TM6) and two functionally important loops with half-helices (LB and LE) which meet in the middle of the membrane (Figure 1). In addition to the above sequences, we also used sequences belonging to fungal XIP subfamily [12] as query sequences. Among the 8 phyla present in the kingdom of fungi [53],[54], at the time of database search, genome sequences of 3 phyla were available either in completed form (35 fungal organisms) or in the form of whole genome shotgun sequences (241 fungi). Our search, as described in the Methods section, yielded 326 fungal MIP sequences. Additional 69 sequences were taken from the studies of Xu et al. [51] and many are found in the database of Joint Genome Institute (JGI; The phylogenetic analysis of all 395 fungal MIP sequences is presented in Figure 2a. It is clearly evident that there are two major clusters representing aquaporins and aquaglyceroporins. The dominant AQP and AQGP clusters contain 163 and 199 members respectively. The third cluster containing the XIP sequences was identified in our earlier studies [12]. In the present study, we report 17 XIP sequences and this group is common to both fungi and plants. This analysis shows another cluster with 16 sequences distinct from the above mentioned three groups. We name this cluster as SIP-like MIP sequences. SIPs form one of the subgroups within plant MIPs [11] and they are small and basic intrinsic proteins. The reasons for naming this fungal MIP cluster as SIP-like sequences are discussed in the following sections. In the case of AQGPs, the subgroups identified in earlier studies [44],[50],[51] are present in our analysis also. We have also identified two additional AQGP subgroups in our analysis (Figure 2b). Sequences, structural models and other associated details of all 395 sequences from 172 different fungal organisms are available in our MIPModDB database ( [6]. The accession codes of these sequences from GenBank or JGI are provided in the Additional file 1: Table S1.

Figure 1
figure 1

Topology, structure and selectivity filter of MIP channels. (a) Topology diagram of a MIP channel exhibiting six transmembrane segments (TM1 to TM6) and the two half-helices formed by the functionally important loops LB and LE. (b) Superposition of three MIP channel structures. Water-transporting AQP1 (brown; PDB ID: 1J4N) and glycerol-specific GlpF (green; PDB ID: 1FX8) are superposed on a modeled fungal MIP structure (purple). Only the helical backbone of transmembrane segments and the loop LB and LE are shown for clarity. (c) Superposition of residues forming the aromatic/arginine selectivity filter from the water-transporting AQP1 (brown) and a model from a representative example of SIP-like fungal MIP channel (purple). Amino acids in one letter codes are labeled in the respective color of each structure. The contributing transmembrane segments (TM2 or TM5) and the loop positions (LE1 or LE2) are also indicated.

Figure 2
figure 2

Phylogenetic analysis of fungal MIP channels. (a) Phylogenetic analysis of 395 fungal MIP channels. The newly identified “SIP-like” fungal MIP channels are distinct from orthodox AQPs, AQGPs and XIPs. (b) Phylogenetic analysis of all fungal AQGPs. The newly identified δ-cluster forms one of the major subfamilies along with α, β, Fps1-like and Yfl054-like AQGPs. γ1 and γ2 form two small clusters. The color of the branches indicates the taxonomy as mentioned in the legends.

Taxonomical distribution of fungal MIPs

Taxonomical distribution of different clusters and the subgroups are also indicated in Figure 2. The orthodox fungal AQPs are present mainly in the species of Ascomycota (yeast and filamentous) and the rest are found in filamentous Basidiomycota. Only two examples are found in Glomeromycota and a lone representative is from Mucoromycotina. All XIPs are found in Ascomycota filamentous fungi. Among the fungal AQGPs, we have identified two additional subgroups in this study apart from the known subgroups. They have different taxonomic distributions. The phylogenetic tree of all fungal AQGPs is presented in Figure 2b. The Fps1-like AQGPs are found invariably in yeast Ascomycota while the Yfl054-like AQGPs, also known as facultative aquaporins [51], are predominantly identified in yeast and filamentous Ascomycota. The other major subgroups classified as α and β are mainly found in Ascomycota (filamentous) and Basidomycota (filamentous) respectively. There are two small subgroups which we have designated as γ1 and γ2. The γ1 cluster has been previously recognized [51] and is mainly found in the species of Mucoromycotina. The second cluster γ2 is newly recognized in this study and is found in filamentous Ascomycota. The major fungal AQGP subgroup newly identified from this analysis consists of 49 members and the members of this group are found in filamentous Ascomycota and Basidiomycota (yeast & filamentous). We have named this cluster as δ. The new fungal MIP group, which is distinct from orthodox AQPs and AQGPs, are found exclusively in Microsporidia and these parasitic species are known to infect their animal hosts. In this context, it is even more important to understand and characterize this MIP group. We have designated this group as “SIP-like” fungal MIPs since they possess several characteristics of plant SIP subgroup [11] (see below). From this analysis, it appears that some fungal MIP subfamilies such as XIPs and the newly identified SIP-like members are exclusive to a single fungal phylum (Figure 2). However, one has to keep in mind that the present analysis was carried out with limited number of available fungal genome sequences. When more fungal genome sequences are available, it is possible that MIP members from these subfamilies may be found in other fungal phyla also.

Characterization of new fungal MIP cluster that is distinct from AQPs, AQGPs and XIPs

Sixteen fungal MIPs form a cluster separate from AQPs, AQGPs and XIPs. Average pairwise sequence identity of all these 16 sequences is about 54% and similarity is close to 70%. We also calculated the average pairwise sequence identity between the sequences from the newly identified cluster and all the sequences belonging to the clusters of orthodox AQPs and XIPs. A similar exercise was carried out for all the subgroups of AQGPs also. The results of pairwise sequence analysis are given in Additional file 2: Table S2. Inter-group average sequence identity varies from 16 to 27% and similarity varies from 32 to 44% indicating that this group of sequences exhibits a large degree of divergence from all other fungal MIP sequences.

MIP subfamilies within plants have been shown to have different biochemical properties [11]. To find out whether the newly identified fungal MIP group exhibits differential properties, we plotted isoelectric point versus molecular weight for each fungal MIP groups and subgroups. We also carried out comparative studies of fungal MIP groups with plant MIP subfamilies (Figure 3). Our analysis shows that members of the newly identified fungal MIP cluster are smaller in size and have molecular weight similar to that found for the plant TIP and SIP subfamilies. While SIP members are more basic, the new fungal MIP sequences have almost equal distribution of acidic and basic residues.

Figure 3
figure 3

Relationship between isoelectric point and molecular weight for the newly identified fungal MIP cluster and selected plant MIP subfamilies. Plot showing isoelectric point versus molecular weight for SIP-like MIPs (red squares), plant TIPs (blue inverted triangles) and plant SIPs (green triangles). One hundred eighty one plant TIPs from 45 species and 25 SIPs from 10 different plant species were considered for this analysis. See also Table 2 for details of the species.

Analysis of Ar/R selectivity filter residues

To further characterize the new subfamily, we constructed homology models of all fungal MIP sequences as described in the Methods section. We particularly focused on the aromatic/arginine selectivity filter residues that have been shown to play significant role in determining the selectivity of the solutes to be transported [55]–[57]. The four residues are contributed by TM2, TM5 and the two loop residues from loop E, LE1 and LE2 (Figure 1). As the name suggests, usually the TM2 position is predominantly occupied by an aromatic residue and the LE2 position is almost always an arginine residue. This is true for majority of other groups of fungal MIP sequences except the newly identified family. It has neither aromatic nor arginine residue in its selectivity filter. Both the TM2 and LE2 positions are occupied by bulky hydrophobic residues and small residues (A/S/G) are found in the other two positions TM5 and LE1 (Table 1). Sequence Logo [58] produced for loop E region also clearly highlights this fact (Figure 4). While all the AQGP subfamilies and AQP and XIP clusters possess the highly conserved arginine which forms part of the selectivity filter, only the newly identified fungal MIP cluster is devoid of this arginine residue in its selectivity filter. Since no fungal MIP subgroup has this feature, we looked at the selectivity filter residues of plant SIP and TIP subfamilies. In both cases, we found members in which aromatic or arginine residue is absent in the selectivity filter (Table 2). A significant number of TIPs do not have arginine as part of the selectivity filter. However, there is a histidine in TM2 position which could provide necessary basic character in the place of arginine. Analysis of 25 plant SIPs from different species reveals that none of them have the arginine residue. Although in some cases an aromatic residue is found in TM2/TM5/LE2 position, the selectivity filter of plant SIP members is in general more hydrophobic as in the newly identified fungal MIP group (Table 2).

Table 1 Aromatic/arginine selectivity filter residues and substitutions in the conserved NPA motifs in different fungal MIP subgroups
Figure 4
figure 4

Sequence logos produced for the loop E region of different fungal MIP groups. The positions corresponding to the conserved NPA motif are used as the reference point and the position of Asn in the NPA motif is designated as ‘0”. All other amino acids are marked relative to this position. The arrows shown in cyan point to the residues from LE1 and LE2 positions that participate in the formation of the narrow selectivity filter. The acidic and basic residues that can form intra-helical salt-bridge interaction are indicated in the figure. The sequences logos were produced from the web server

Table 2 Aromatic/arginine selectivity filter residues and substitutions in the conserved NPA motifs in plant subfamilies TIPs and SIPs

Substitutions in the highly conserved NPA boxes

The third factor we examined is the substitution in the highly conserved NPA motifs. In majority of the members belonging to the newly identified fungal MIP cluster, alanine residue of NPA motif from loop E has been substituted (Table 1). Substitution of conserved NPA motif is also common in plant SIP subfamily (Table 2). While in almost all plant TIPs both NPA motifs are strictly conserved, substitution in NPA motif of loop B is observed in 21 out of 25 SIP members. Thus the newly identified fungal MIP group resembles more like the plant SIP subfamily in terms of the molecular weight, nature of selectivity filter residues and substitution in the highly conserved NPA motifs. Hence, we have classified this new fungal MIP cluster as “SIP-like” subfamily taking all these factors into account. In our previous studies we have demonstrated that MIP sequences show high group conservation of residues at 17 positions which occur at the helix-helix interface of the hour-glass helical fold [12],[59]. Analysis of the SIP-like members in fungi reveals that 14 out of 17 positions exhibit very high group conservation confirming that the new members have characteristic features of MIP superfamily (Additional file 3: Table S3). Further analysis shows that the SIP-like subfamily members have small N- and C-termini regions and the loop D connecting TM4 and TM5 is the longest in this family compared to all other fungal MIP clusters (Figure 5). This loop in SIP-like MIPs also possesses highest number of charged residues among all fungal MIP subgroups.

Figure 5
figure 5

Analysis of loop D in various fungal MIP groups. (a) Average length and (b) Average number of charged residues of loop D connecting the transmembrane segments TM4 and TM5 are shown for different fungal MIP subfamilies.

Characterization of δ subfamily of fungal AQGPs

The newly identified SIP-like fungal MIP family is distinct from AQPs, AQGPs and XIPs within fungal MIP family. Within the fungal AQGP group, five different subgroups have already been identified in a recently published study [51]. In the current study, we have also identified two additional subgroups distinct from the previously identified fungal AQGP clusters. The biggest cluster is δ with 49 members (Figure 2b). The average pairwise sequence identity and similarity within this group are about 51% and 67% respectively. This indicates that the sequences within the group are closely related. However, the same cannot be said for AQGP members from other subgroups. The inter-group average sequence identity between δ cluster and other AQGP groups varies from 24% to 30% indicating that δ subgroup members have diverged significantly and are distantly related to other fungal AQGP groups (Additional file 2: Table S2). Hence, our study mainly focused on features that are unique to δ cluster AQGPs and hence, we present our analysis by comparing δ subgroup with other fungal AQGP groups.

We have analyzed the molecular weight of these sequences as a function of isoelectric point and compared them with other fungal AQGP subgroups. We did not find any feature in this plot that distinguishes δ from other fungal AQGP clusters (Additional file 4: Figure S1). We then analyzed the selectivity filter residues from the homology models we have generated. All the members from this subgroup contain arginine in the LE2 position and the TM2 position is predominantly occupied by an aromatic residue (Table 1). The other two positions from TM5 and LE1 are preferred by bulky aliphatic residues Ile, Leu and Val (Figure 6). This results in a relatively more hydrophobic environment in the selectivity filter region. This also distinguishes δ from α, β, and γ2 members in which two aromatic residues are found in TM2 and LE1 positions and usually a small residue is found in TM5 position (Table 1). In the subgroups Fps1-like and facultative aquaporins, at least one position is occupied by a small residue. Thus it is clear that the δ subgroup MIP members have selectivity filter distinct from other fungal AQGPs. We anticipate that members from this group are likely to transport relatively more hydrophobic solutes.

Figure 6
figure 6

Aromatic/arginine selectivity filter is novel in δ AQGPs. Superposition of selectivity filter residues from glycerol-transporting GlpF (brown; PDB ID: 1FX8) and one of the modeled fungal AQGP member from the δ subgroup (purple).

We have examined the conservation of the signature NPA motifs in loops LB and LE and compared the pattern with other AQGPs. While the NPA motif of loop LB is substituted in some members of δ cluster, the proline residue is invariably replaced in almost every member in the second NPA motif that occurs in LE (Table 1 and Figure 4). This is also observed in α cluster and Fps1-like members. Analysis of sequences in the loop B region reveals that the highly conserved Asn in the NPA motif is substituted by Ser or His in some of δ AQGPs indicating that side-chains of these residues are likely to play a similar role as that of Asn (Additional file 5: Figure S2). Sequence logo of loop LE region reveals strict conservation of an aspartate and an arginine residue near the NPA motif of loop LE in all AQGP members (Figure 4). These two residues are separated by four positions enabling them to form a salt-bridge interaction in the half-helix formed by the LE loop (Figure 7). Such an interaction gives additional stability to the helical region within LE and the residues forming these interactions do not face the channel interior. Such a stabilizing interaction seems to be unique to AQGPs from fungi and other organisms (R. K. Verma, N. D. Prabh and R. Sankararamakrishnan, Unpublished results). The acidic and basic residues are not conserved in orthodox AQPs from fungi, XIPs and the newly identified SIP-like MIPs (Figure 4). Molecular dynamics simulations of mammalian AQP, E. coli AQGP and Plasmodium AQGP reveal that while the half-helix LE is stable in E.coli and Plasmodium AQGPs, unwinding of LE half-helix is observed in mammalian AQP [60]. In AQGPs from E. coli and Plasmodium, the intra-helical salt-bridge interaction gives additional stability to the half-helix in loop LE. Since loop LE contributes two out of four residues (LE1 and LE2) for the Ar/R selectivity filter, we speculate that the unwinding character may help in regulating the channel transport in orthodox AQPs that lack the intra-helical salt-bridge interaction. However in the case of fungal AQGPs, presence of this interaction seems to maintain the helix stable and there could be other mechanisms to regulate the channel function. In this respect, it is important to note the two highly conserved residues, glycine and proline, are present between the acidic and basic residue positions in all fungal AQGP groups except the δ cluster of AQGPs (Figure 4). Both these residues are known to be helix breakers. It appears that the helix stability of loop E half-helix is balanced by two opposing forces, one with stabilizing interaction and the other with the tendency to break helices. In the case of δ cluster, the proline residue is not conserved (Figure 4). Hence, we can anticipate that the half-helix in loop E in δ cluster will be relatively more stable than other fungal AQGPs.

Figure 7
figure 7

Intra-helical salt-bridge in the loop E half-helix of fungal AQGPs. MIP channels from representative members from (a) SIP-like, (b) δ AQGP, (c) α AQGP and (d) orthodox AQP clusters. The view is down the channel axis from the extracellular side. The residues forming the selectivity filter are shown in green in stick representation. Loop E is shown in dark brown. Acidic and basic residues forming the intra-helical salt-bridge are shown for the δ and α members. The equivalent positions are shown in SIP-like and orthodox AQP members. With the exception of the δ-subgroup, a proline residue occurs at the intervening positions of the acidic and basic residues and it is highly conserved in all other AQGP subgroups and also in most of the orthodox AQPs. This residue is shown in blue and it is substituted by Ala in the representative δ MIP model shown here.

γ2 cluster of AQGPs

This small group contains only 5 members forming a separate cluster within the fungal AQGPs (Figure 2b) and all of them belong to filamentous Ascomycota. Their selectivity filter contains a Tyr residue at LE1 position which seems to be a feature shared by many members of α and β clusters (Table 1). The TM2 position is preferred by a bulky residue (Trp/Met/Ile) and LE2 position is occupied by Arg. There is no strikingly different feature observed in this small group. Perhaps, a larger dataset would help to identify unique characteristics of small clusters like γ1 and γ2.


Fungal MIPs are as diverse as plant MIPs

The diversity of plant MIPs has been established in many earlier studies [10]–[12]. Only few studies identified and addressed the extent of diversity in fungal MIPs [44],[50],[51]. In the present study, by identifying two major clusters and one minor cluster within the fungal MIP superfamily, we have determined that fungal MIPs are as diverse as their plant counterparts. Till date, the only MIP subgroup common to both plant and fungal species has been found to be XIPs [12]. However, by systematically searching the genome sequences of fungal species, we have recognized another group in fungal MIPs which has common features with plant SIP subfamily. Although phylogenetic analysis cluster the newly identified “SIP-like” fungal MIPs separately from the plant SIPs, several features are shared by the two MIP families from the two different species group. The small size of the channel protein, the nature of selectivity filter residues and substitutions in the conserved NPA boxes indicate that the new fungal MIP cluster shares some characteristic features with plant SIPs.

Unique selectivity filter suggests that SIP-like MIPs are likely to transport larger hydrophobic solutes

Very few functional studies have been carried on fungal MIPs [39],[41],[44],[45]. Experimental, structural and simulation studies have emphasized the role of the aromatic/arginine selectivity filter in MIP channel’s transport and selectivity [55]–[57],[61]–[64]. In this context, it would be interesting to at least make a speculation regarding the transport properties of SIP-like fungal MIPs based on the selectivity filter residues. Even for plant SIPs, not many studies have been reported regarding their function. When Arabidopsis SIPs (AtSIPs) were expressed in yeast, at least one of the two members displayed water channel activity [65]. This water transport has been reported in spite of the fact that the aromatic/arginine selectivity filter of AtSIPs lacks the arginine residue. In a recent study [34], an aquaporin AQP2 from the protozoan parasite Trypanosoma brucei gambiense or Trypanosoma brucei rhodesiense has been shown to have an unusual selectivity filter comprising of large hydrophobic residues IVLL. T. brucei causes sleeping sickness in humans and with all hydrophobic residues in the selectivity filter without any aromatic and arginine residue, this aquaporin has been shown to play a vital role in increased cross resistance to the drugs melarsoprol and pentamidine [66]. Based on their studies, Baker et al. [34] have suggested that T. brucei aquaporin AQP2 may function as a transporter of pentamidine and melarsoprol drugs. Both drugs with molecular mass above 330 Da are larger molecules compared to glycerol (92 Da). It is suggested that the unique selectivity filter of AQP2 from T. Brucei with bulky hydrophobic residues may facilitate the transport of large molecules.

We have also examined all non-plant MIPs whose selectivity filter lacks arginine residue. In MIPModDB [6] which has more than 1000 MIPs, we could find only 52 examples in which arginine is absent in the selectivity filter. Among them, only 17 of them including mammalian AQP12 have features similar to that found in SIP-like fungal MIPs (two bulky hydrophobic residues and two small residues in the selectivity filter). We found that in mammalian AQP12, both TM2 and LE2 positions are occupied by Leu residues and Ala is found in TM5 and LE1 positions. Although no direct functional studies are reported for AQP12, it has been speculated that AQP12 may regulate proper secretion of pancreatic fluid or digestive enzymes [67],[68]. Thus, plant SIPs, AQP12 and AQP2 of T. Brucei also have selectivity filters similar to fungal SIP-like MIPs with all four positions occupied by hydrophobic residues and arginine is absent. Based on the available literature on other MIPs with similar features, we conclude that the SIP-like fungal MIPs are likely to transport solutes which are novel, hydrophobic and larger in size. The same argument may hold good for the plant SIPs also.

Loop D in SIP-like members may play a role in regulating the channel function

The SIP-like fungal MIPs exhibit some characteristic features in loop D which connects the TM4 and TM5 helices. The average length of loop D is the longest in SIP-like cluster compared to all other fungal MIP subgroups (Figure 5). Moreover, large number of charged residues is found in loop D in SIP-like members. This observation indicates that loop D in SIP-like MIPs can play a special role in the function of these channels. Loop D in some mammalian and plant aquaporin channels have been shown to be important in channel gating and regulation of channel function. Experimental and computational studies in AQP1, spinach aquaporin SoPIP2;1 and AQP4 have suggested a role for loop D in channel gating [69],[70] or binding site for channel agonists [71] or antagonists [72]. Specific residues in loop D which respond to pH have been shown to be responsible for the closed state in AQP1 channel [73]. The arginine-rich loop D in AQP1 has also been shown to bind cGMP which seems to initiate channel opening [74]. Loop D in AQP4 serves as a possible region for metal-binding activity and has also been implicated in regulating the channel activity [75],[76]. Structural studies on spinach aquaporin characterized the open and closed states and exhibited two different conformations of loop D revealing that the gating mechanism at molecular level involves a large movement of this loop [77]. The available experimental and simulation studies provide a clue for the role of loop D in fungal SIP-like MIPs. The charged residues in this loop could be involved in binding ions or molecules that could trigger the channel opening and closing. The longer loop in this subfamily can easily adopt different conformations that can regulate the channel transport. Both possibilities can also simultaneously exist implying that the loop D in SIP-like channels should be given specific attention while investigating the structure-function relationship of these subgroup members.

SIP-like channels could be attractive targets for some anti-fungal diseases

It is interesting to note that all the SIP-like members are from the lower fungi Microsporidia. We did not find even a single example from higher fungi. Although Micorsporidia were earlier classified as part of the kingdom protozoa, these spore-forming unicellular parasites are now known to be fungi. They infect both humans and animals and cause microsporidiosis in humans [78]. As human pathogens they cause diarrhea and infections in immune-compromised individuals (Additional file 6: Table S4). Hence, understanding the host-pathogen interactions and the infection mechanism is extremely important that will aid in developing anti-fungal drugs. SIP-like members in Microsporidia are unique in several respects. They are small in size and the channels formed by them have hydrophobic selectivity filters. The loop D in these channels is longer and is highly charged. The nature of molecules that are transported through these channels and the role of loop D in regulating the channel function have to be first experimentally investigated. Although human AQP12 has similar selectivity filter, phylogenetic analysis clearly clusters them into separate clades. Hence, SIP-like MIP channels in disease causing Microsporidial fungi can be considered as attractive drug targets.

Fungal AQGPs have the largest number of subfamilies

With the identification of δ and γ2 AQGP subfamilies in this study, fungal AQGPs can be classified into seven different subgroups (Fps1-like, facultative, α, β, γ1, γ2 and δ). To our knowledge, the kingdom of fungi possesses the largest subfamilies within the major group of AQGPs. This makes the fungal AQGPs as the most diverse among all known AQGPs compared to any major organism groups. In majority of fungal AQGPs, Trp is found to be preferred in the TM2 position of the selectivity filter and Arg is absolutely conserved in LE2 position of all fungal AQGPs (Table 1). Most of the variations are seen in TM5 and LE1 positions. Members of γ1 AQGPs and the newly identified δ subgroup of AQGPs defy the trend observed in TM2 position. While Trp is absent in members of both the clusters, a small polar residue (Thr, Ser or Asn) is found in γ1 at TM2 position. However in δ subgroup, bulky residues of hydrophobic nature (Phe, Tyr, Met, Val, Ile and Leu) are found in TM2, TM5 and LE1 positions. As in SIP-like MIPs, the selectivity filter of δ AQGPs appears to be unique. We searched MIPModDB database with more than 1000 MIP members looking for examples with similar features in their selectivity filters. We could find only four MIPs in which the channel selectivity filter is formed by three bulky hydrophobic residues and the fourth position is occupied by an arginine residue.

As mentioned earlier, δ subgroup of AQGPs has been found in filamentous Ascomycota and yeast and filamentous Basidiomycota. Many examples from this species group are known to be plant pathogens [79] (Additional file 7: Table S5). Members belonging to this group have been identified to cause serious diseases in wide range of plants including agriculturally important cereals such as rice, wheat, rye and barley. They also affect economically important crops like cotton. These pathogens display resistance to fungicides and δ subgroup of AQGPs can be considered as an important target to contain these plant pathogens.

Loop E as a marker to study the evolution of MIPs

Loop E provides two residues for the selectivity filter. It also possesses one of the two conserved NPA motifs. While the other conserved NPA motif resides in loop B, majority of the positions in this loop show variations (Additional file 5: Figure S2). In the case of loop E, this region seems to show some distinctive features depending upon the subgroups. Six out of seven AQGPs exhibit a conserved motif RDxGPR next to the NPA motif in loop E (Figure 4). A very high conservation of acidic and basic residues separated by four positions will enable the formation of salt-bridge interaction in the half-helix. However, as mentioned earlier, the presence of two helix destabilizing residues Gly and Pro will weaken the stability of the same helix. In δ AQGPs, the motif present is RDxGxR and the Pro residue next to Gly is mostly substituted by other residues. Hence, we can think of the half-helix in δ AQGPs as relatively more stable compared to all other fungal AQGP subgroups. In XIPs, the acidic and basic residues are replaced by a strictly conserved cysteine and mostly bulky residues (Tyr, Phe or Met) respectively. The motif in XIPs is found to be RCx[G/A]xx. SIP-like members exhibit no such conserved motif in this region. In orthodox AQPs, the acidic and basic residues are replaced by small residues and the conservation pattern in the same region is R[S/A]xG[P/A][A/C/D/S]. It should be noted that the absence of intra-helical salt-bridge and the conservation of Gly and to some extent Pro in AQPs will render this half-helix relatively less stable compared to the same region in AQGPs in general and δ-AQGPs in particular. An extensive bioinformatics analysis of MIPs from microbial organisms to mammals in loop E is necessary and this may spring some surprising insights into the evolution of MIPs.


The role of different fungal proteins in the symbiotic relationship with plants needs to be clearly established. Similarly, host-pathogen interactions at the time of fungi-induced infection require understanding at the molecular level. In this context, understanding the evolution, function and diversity of fungal MIP channels has become very significant. Experimental studies clearly reveal that MIP channels are important players in plant-fungi interactions. Identification of new targets for anti-fungal drugs is a major goal for several human diseases caused by fungi and MIP channels could be considered as attractive anti-fungal targets. In the present study, analysis of fungal genome sequences has identified additional MIP channels in different species groups within the kingdom of fungi. Phylogenetic analysis of nearly 400 fungal MIP channels has revealed the existence of a new MIP cluster completely distinct from the orthodox AQP and AQGP channels. Further sequence analysis and homology modeling studies indicate features that are shared between the new fungal MIP cluster and the plant SIP subfamily. The size of the protein, chemical nature of the residues that form the narrow aromatic/arginine selectivity filter and the substitutions found in the conserved NPA motifs are the common characteristics between the new family and the plant SIP channels. Hence in addition to the XIPs, this “SIP-like” fungal MIP channels can possibly be another evolutionary link between the plants and fungi. Since SIP-like channels are observed only in Microsporidia which are the unicellular parasites, they can also be considered as drug targets for developing anti-fungal drugs in human infections caused by fungi.

We have also identified one major subgroup and another minor cluster within the fungal AQGP family. The δ-subgroup with 49 members have unique selectivity filter very rarely found in other MIP channels. With the discovery of these two new AQGP subfamilies, it appears that fungal AQGPs are the most diverse among all known AQGPs. Many fungal species possessing δ-AQGPs also are known to act as plant pathogens. The newly identified δ-group MIP channels can be exploited to contain some of the serious infections affecting agriculturally and economically important crops.

With very few functional studies available in fungal MIP channels, the present study has given a picture of the diverse fungal organisms evolved with two new MIP subfamilies with unique selectivity filter residues. The possible solutes that are transported through these narrow regions within the channel interior, the role of loop D in SIP-like channels and the importance of these channels in the fungal life cycles are some of the immediate questions that have to be addressed by researchers in this field.


Identification of fungal MIP sequences

The non-redundant database of NCBI ( and the UniProtKB ( database [52] were searched using known AQP and AQGP sequences as query. The UniProt accession IDs of these sequences are P0AER0, P47865, P60844, Q6J8I9, Q9C4Z5, Q41372, Q8WPZ6, P55064, P55087 and F2QVG4. These sequences belong to different species and organism groups and they come from bacteria (Escherichia coli), mammals (Bos taurus, Ovis aries and Homo sapiens), archaea (Methanothermobacter marburgensis), plant (Spinacia oleracea), protozoan parasite (Plasmodium falciparum) and yeast (Komagataella pastoris). The three-dimensional structures of all these channels have been determined. In addition to these sequences, we have also considered the sequences belonging to the XIP subfamily from fungal species as query sequences previously identified in our laboratory [12]. We used PSI-BLAST [80] on NCBI non-redundant database to identify fungal MIP sequences. For each sequence, three iterations were performed with threshold e-value 0.001. When going from previous iteration to the next iteration, we included only those hits whose expect value was 0.001 or less to generate the position-specific scoring matrix that can be used in the next iteration. This is to make sure that the PSSM is not corrupted and at the same time diverse MIP sequences can be obtained as hits. BLASTp was used with UniProt database with the same query sequences. Additionally, we also used tBLASTn against NCBI GenBank database to retrieve additional fungal MIP sequences. The newly identified sequences from PSI-BLAST and tBLASTn searches were again used as query sequences and second round of search was carried out using the same tools. This exercise was repeated until no new MIP sequences were found.

We then used the program CD-HIT [81],[82] on the set all fungal MIP sequences thus obtained to remove any redundancy. It is possible that some of the MIP sequences in this non-redundant set may be partial or may not have all the features associated with a MIP channel protein. In order to find out the partial sequences or sequences that don’t have all MIP features, we adopted the following strategy. We used a multiple sequence alignment program PRALINE ( to align all the fungal MIP sequences [83]. There are advantages using PRALINE. It uses profile-based approach to align sequences that have low sequence identity. It also guides the sequence alignment using secondary structure information by employing different scoring matrices for helix, strand or coil. For each MIP sequence, we used TMHMM [84] and PSIPRED [85] which are available as part of the PRALINE toolkit to predict transmembrane regions and secondary structures respectively. The multiple sequence alignment thus created was examined to find out the following features specific to MIP family sequences. These are (a) presence of two NPA or NPA-like motifs, (b) presence of six transmembrane segments and two functionally important loops possessing the features characteristically present in MIP channels and (c) group-based conservation of small and weakly polar residues in most of the 17 positions that typically occur in the helix-helix interface of MIP hour-glass helical fold [59]. Only those sequences which satisfy all the three conditions were considered further for analysis.

We identified 326 fungal MIP sequences from the database search described above. We also included additional 69 fungal MIP sequences from the recently published work of Xu et al. [51] which fulfilled all our criteria. In total, 395 fungal MIP sequences were considered for phylogenetic analysis.

Phylogenetic analysis

To understand the fungal MIP diversity, the evolution of MIP family and to enable comparison with other groups of species, we performed phylogenetic analysis on all fungal MIP sequences using MEGA 5.0 suite of software [86]. Two different clustering algorithms, namely neighbor-joining and maximum parsimony methods, were used to derive a phylogenetic tree. Reliability of individual branches of the tree was estimated by performing bootstrapping with 1000 replicates. We applied 50% majority rule so that branches with less than 50% confidence level were collapsed to get the final tree topology. The fungal MIP subgroups obtained from this tree was considered for further analysis. The groupings obtained from this tree were validated by comparing the trees obtained from the two different clustering methods. We also used different multiple sequence alignments as input for MEGA 5.0 obtained from various MSA tools like Clustal-W [87], Clustal-Ω [88] and MUSCLE [89] and the resulting trees generated from the newly obtained multiple sequence alignments were compared.

Homology modeling of fungal MIP sequences

Three-dimensional structures of all fungal MIP sequences were modeled using the protocol developed in our laboratory earlier [12],[59]. Briefly, the software MODELLER version 9.10 [90],[91] was used to build the structure and the experimentally determined high-resolution MIP structures from mammalian AQP1 (PDB ID:1J4N) [92], bacterial GlpF (PDB ID: 1FX8) [93] and archael AQPM (PDB ID: 2F2B) [94] were used as templates. The target-template sequence alignment was manually examined to find out whether there are gaps in the transmembrane helical regions or in the functionally important loop regions which possess the highly conserved NPA motifs. We also looked at the high conservation of at least one residue in each transmembrane helical segment (see above). Among the ten models built for each fungal MIP sequence, the one with the optimum MODELLER objective function was selected. Loops and side-chain conformations of non-conserved residues were further refined using the MODELLER’s loop optimization protocol and the program SCWRL3 [95] respectively. The resultant model was minimized using GROMACS ver. 4.5 [96] and the quality of the model was examined using PROCHECK [97].

No human or animal experiments were carried out in this study.

Additional files



Major intrinsic proteins






X-intrinsic proteins


Small and basic intrinsic proteins


Tonoplast intrinsic proteins


  1. Gomes D, Agasse A, Thiebaud P, Delrot S, Geros H, Chaumont F: Aquaporins are multifunctional water and solute transporters highly divergent in living organisms. Biochim Biophys Acta. 2009, 1788: 1213-1228. 10.1016/j.bbamem.2009.03.009.

    Article  PubMed  CAS  Google Scholar 

  2. Maurel C: Plant aquaporins: Novel functions and regulation properties. FEBS Lett. 2007, 581: 2227-2236. 10.1016/j.febslet.2007.03.021.

    Article  PubMed  CAS  Google Scholar 

  3. Wu B, Beitz E: Aquaporins with selectivity for unconventional permeants. Cell Mol Life Sci. 2007, 64: 2413-2421. 10.1007/s00018-007-7163-2.

    Article  PubMed  CAS  Google Scholar 

  4. Bienert GP, Schussler MD, Jahn TP: Metalloids: essential, beneficial or toxic? Major intrinsic proteins sort it out. Trends Biochem Sci. 2008, 33: 20-26. 10.1016/j.tibs.2007.10.004.

    Article  PubMed  CAS  Google Scholar 

  5. Abascal F, Irisarri I, Zardoya R: Diversity and evolution of membrane intrinsic proteins. Biochim Biophys Acta. 2014, 1840: 1468-1481. 10.1016/j.bbagen.2013.12.001.

    Article  PubMed  CAS  Google Scholar 

  6. Gupta AB, Verma RK, Agarwal V, Vajpai M, Bansal V, Sankararamakrishnan R: MIPModDB: A central resource for the superfamily of major intrinsic proteins. Nucleic Acids Res. 2012, 40: D362-D369. 10.1093/nar/gkr914.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  7. Borgnia MJ, Agre P: Reconstitution and functional comparison of purified GlpF and AqpZ, the glycerol and water channels from Escherichia coli. Proc Natl Acad Sci U S A. 2001, 98: 2888-2893. 10.1073/pnas.051628098.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  8. Rojek A, Praetorius J, Frokiaer J, Nielsen S, Fenton RA: A current view of the mammalian aquaglyceroporins. Annu Rev Physiol. 2008, 70: 301-327. 10.1146/annurev.physiol.70.113006.100452.

    Article  PubMed  CAS  Google Scholar 

  9. Ishibashi K, Hara S, Kondo S: Aquaporin water channels in mammals. Clin Exp Nephrol. 2009, 13: 107-117. 10.1007/s10157-008-0118-6.

    Article  PubMed  CAS  Google Scholar 

  10. 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. 10.1104/pp.126.4.1358.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  11. Johanson U, Gustavsson S: A new subfamily of major intrinsic proteins in plants. Mol Biol Evol. 2002, 19: 456-461. 10.1093/oxfordjournals.molbev.a004101.

    Article  PubMed  CAS  Google Scholar 

  12. Gupta AB, Sankararamakrishnan R: Genome-wide analysis of major intrinsic proteins in the tree plant Populus trichocarpa: Characterization of XIP subfamily of aquaporins from evolutionary perspective.BMC Plant Biol 2009, 9:Art. no. 134.,

  13. Nielsen S, Frokiaer J, Marples D, Kwon TH, Agre P, Knepper MA: Aquaporins in kidney: From molecules to medicine. Physiol Rev. 2002, 82: 205-244.

    Article  PubMed  CAS  Google Scholar 

  14. Papadopoulos MC, Verkman AS: Aquaporin water channels in the nervous system. Nature Rev Neurosci. 2013, 14: 265-277. 10.1038/nrn3468.

    Article  CAS  Google Scholar 

  15. Verkman AS: Aquaporins in clinical medicine. Annu Rev Med. 2012, 63: 303-316. 10.1146/annurev-med-043010-193843.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  16. Agre P, King LS, Yasui M, Guggino WB, Otterson OP, Fujiyoshi Y, Engel A, Nielsen S: Aquaporin water channels - from atomic structure to clinical medicine. J Physiol London. 2002, 542: 3-16. 10.1113/jphysiol.2002.020818.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  17. Lennon VA, Kryzer TJ, Pittock SJ, Verkman AS, Hinson SR: IgG marker of optic-spinal multiple sclerosis binds to the aquaporin-4 water channel. J Exp Med. 2005, 202: 473-477. 10.1084/jem.20050304.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  18. Verkman AS, Hara-Chikuma M, Papadopoulos MC: Aquaporins - new players in cancer biology. J Mol Med. 2008, 86: 523-529. 10.1007/s00109-008-0303-9.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  19. Verkman AS: Role of aquaporin water channels in eye function. Exp Eye Res. 2003, 76: 137-143. 10.1016/S0014-4835(02)00303-2.

    Article  PubMed  CAS  Google Scholar 

  20. Binder DK, Nagelhus EA, Ottersen OP: Aquaporin-4 and epilepsy. Glia. 2012, 60: 1203-1214. 10.1002/glia.22317.

    Article  PubMed  Google Scholar 

  21. Hibuse T, Maeda N, Funahashi T, Yamamoto K, Nagasawa A, Mizunoya W, Kishida K, Inoue K, Kuriyama H, Nakamura T, Fushiki T, Kihara S, Shimomura L: Aquaporin 7 deficiency is associated with development of obesity through activation of adipose glycerol kinase. Proc Natl Acad Sci U S A. 2005, 102: 10993-10998. 10.1073/pnas.0503291102.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  22. Loonen AJM, Knoers NVAM, van Os CH, Deen PMT: Aquaporin 2 mutations in nephrogenic diabetes insipidus. Seminars Nephrol. 2008, 26: 252-265. 10.1016/j.semnephrol.2008.03.006.

    Article  Google Scholar 

  23. Maurel C, Verdoucq L, Luu DT, Santoni V: Plant aquaporins: Membrane channels with multiple integrated functions. Annu Rev Plant Biol. 2008, 59: 595-624. 10.1146/annurev.arplant.59.032607.092734.

    Article  PubMed  CAS  Google Scholar 

  24. 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, 1465: 324-342. 10.1016/S0005-2736(00)00147-4.

    Article  PubMed  CAS  Google Scholar 

  25. Tyerman SD, Niemietz CM, Bramley H: Plant aquaporins: multifunctional water and solute channels with expanding roles. Plant Cell Environ. 2002, 25: 173-194. 10.1046/j.0016-8025.2001.00791.x.

    Article  PubMed  CAS  Google Scholar 

  26. Sade N, Vinocur BJ, Diber A, Shatil A, Ronen G, Nissan H, Wallach R, Karchi H, Moshelion M: Improving plant stress tolerance and yield production: Is the tonoplat aquaporin SITIP2;2 a key to isohydric to anisohydric conversion?. New Phytol. 2009, 181: 651-661. 10.1111/j.1469-8137.2008.02689.x.

    Article  PubMed  CAS  Google Scholar 

  27. Marjanovic Z, Uehlein N, Kaldenhoff R, Zwiazek JJ, Weiss M, Hampp R, Nehls U: Aquaporins in poplar: What a difference a symbiont makes!. Planta. 2005, 222: 258-268. 10.1007/s00425-005-1539-z.

    Article  PubMed  CAS  Google Scholar 

  28. Uehlein N, Fileschi K, Eckert M, Bienert GP, Bertl A, Kaldenhoff R: Arbuscular mycorrhizal symbiosis and plant aquaporin expression. Phytochemistry. 2007, 68: 122-129. 10.1016/j.phytochem.2006.09.033.

    Article  PubMed  CAS  Google Scholar 

  29. Dowd C, Wilson LW, McFadden H: Gene expression profile changes in cotton root and hypocotyl tissues in response to infection with Fusarium oxysporum f.sp vasinfectum. Mol Plant Microbe Interact. 2004, 17: 654-667.

    Article  PubMed  CAS  Google Scholar 

  30. Fadiel A, Isokpehi RD, Stambouli N, Hamza A, Benammar-Elgaaied A, Scalise TJ: Protozoan parasite aquaporins. Expert Rev Proteomics. 2009, 6: 199-211. 10.1586/epr.09.10.

    Article  PubMed  CAS  Google Scholar 

  31. Beitz E: Aquaporins from pathogenic protozoan parasites: structure, function and potential for chemotherapy. Biol Cell. 2005, 97: 373-383. 10.1042/BC20040095.

    Article  PubMed  CAS  Google Scholar 

  32. Hansen M, Kun JFJ, Schultz JE, Beitz E: A single, bi-functional aquaglyceroporin in blood-stage Plasmodium falciparum malaria parasites. J Biol Chem. 2002, 277: 4874-4882. 10.1074/jbc.M110683200.

    Article  PubMed  CAS  Google Scholar 

  33. Montalvetti A, Rohloff P, Docampo R: A functional aquaporin co-localizes with the vacuolar proton pyrophosphatase to acidocalcisomes and the contractile vacuole complex of Trypanosoma cruzi. J Biol Chem. 2004, 279: 38673-38682. 10.1074/jbc.M406304200.

    Article  PubMed  CAS  Google Scholar 

  34. Baker N, Glover L, Munday JC, Andres DA, Barrett MP, de Koning HP, Horn D: Aquaglyceroporin 2 controls susceptibility to melarsoprol and pentamidine in African tyrpanosomes. Proc Natl Acad Sci U S A. 2012, 109: 10996-11001. 10.1073/pnas.1202885109.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  35. Kun JF, de Carvalho EG: Novel therapeutic targets in Plasmodium falciparum: aquaglyceroporins. Expert Opin Ther Targets. 2009, 13: 385-394. 10.1517/14728220902817839.

    Article  PubMed  CAS  Google Scholar 

  36. Bonfante P, Genre A: Mechanisms underlying beneficial plant-fungus interactions in mycorrhizal symbiosis.Nature Commn 2010, 1:Art. No. 48.,

  37. Casieri L, Lahmidi NA, Doidy J, Veneault-Fourrey C, Migeon A, Bonneau L, Courty PE, Garcia K, Charbonnier M, Delteil A, Brun A, Zimmermann S, Plassard C, Wipf D: Biotrophic transportome in mutualistic plant-fungal interactions. Mycorrhiza. 2013, 23: 597-625. 10.1007/s00572-013-0496-9.

    Article  PubMed  CAS  Google Scholar 

  38. Parniske M: Arbuscular mycorrhiza: the mother of plant root endosymbioses. Nature Rev Microbiol. 2008, 6: 763-775. 10.1038/nrmicro1987.

    Article  CAS  Google Scholar 

  39. Li T, Hu YJ, Hao ZP, Li H, Wang YS, Chen BD: First cloning and characterization of two functional aquaporin genes from an arbuscular mycorrhizal fungus Glomus intraradices. New Phytol. 2013, 197: 617-630. 10.1111/nph.12011.

    Article  PubMed  CAS  Google Scholar 

  40. Giovannetti M, Balestrini R, Volpe V, Guether M, Straub D, Costa A, Ludewig U, Bonfante P: Two putative aquaporin genes are differently expressed during arbuscular mycorrhizal symbiosis inLotus japonicus.BMC Plant Biol 2012, 12:Art. No: 186.,

  41. Navarro-Rodenas A, Ruiz-Lozano JM, Kaldenhoff R, Morte A: The aquaporin TcAQP1 of the desert truffle Terfezia claveryi is a membrane pore for water and CO2 transport. Mol Plant-Microbe Interact. 2012, 25: 259-266. 10.1094/MPMI-07-11-0190.

    Article  PubMed  CAS  Google Scholar 

  42. Lee SH, Calvo-Polanco M, Chung GC, Zwiazek JJ: Role of aquaporins in root water transport of ectomycorrhizal jack pine (Pinus banksiana) seedlings exposed to NaCl and fluoride. Plant Cell Environ. 2010, 33: 769-780.

    PubMed  CAS  Google Scholar 

  43. Aroca R, Porcel R, Ruiz-Lozano JM: How does arbuscular mycorrhizal symbiosis regulate root hydraulic properties and plasma membrane aquaporins in Phaseolus vulgaris under drought, cold or salinity stresses?. New Phytol. 2007, 173: 808-816. 10.1111/j.1469-8137.2006.01961.x.

    Article  PubMed  CAS  Google Scholar 

  44. Dietz S, von Bulow J, Baker N, Nehls U: The aquaporin gene family of the ectomycorrhizal fungus Laccaria bicolor: lessons for symbiotic functions. New Phytol. 2011, 190: 927-940. 10.1111/j.1469-8137.2011.03651.x.

    Article  PubMed  CAS  Google Scholar 

  45. Aroca R, Bago A, Sutka M, Paz JA, Cano C, Amodeo G, Ruiz-Lozano JM: Expression analysis of the first arbuscular mycorrhizal fungi aquaporin described reveals concerted gene expression between salt-stressed and nonstressed mycelium. Mol Plant Microbe Interact. 2009, 22: 1169-1178.

    Article  PubMed  CAS  Google Scholar 

  46. Castle NA: Aquaporins as targets for drug discovery. Drug Disc Today. 2005, 10: 485-493. 10.1016/S1359-6446(05)03390-8.

    Article  CAS  Google Scholar 

  47. Kontoyiannis DP, Lewis RE: Antifungal drug resistance of pathogenic fungi. Lancet. 2002, 359: 1135-1144. 10.1016/S0140-6736(02)08162-X.

    Article  PubMed  CAS  Google Scholar 

  48. Sanglard D: Resistance of human fungal pathogens to antifungal drugs. Curr Opin Microbiol. 2002, 5: 379-385. 10.1016/S1369-5274(02)00344-2.

    Article  PubMed  CAS  Google Scholar 

  49. Ghosh K, Cappiello CD, McBride SM, Occi JL, Cali A, Takvorian PM, McDonald TV, Weiss LM: Functional characterization of a putative aquaporin from Encephalitozoon cuniculi, a microsporidia pathogenic to humans. Int J Parasitol. 2006, 36: 57-62. 10.1016/j.ijpara.2005.08.013.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  50. Pettersson N, Filipsson C, Becit E, Brive L, Hohmann S: Aquaporins in yeasts and filamentous fungi. Biol Cell. 2005, 97: 487-500. 10.1042/BC20040144.

    Article  PubMed  CAS  Google Scholar 

  51. Xu H, Cooke JEK, Zwiazek JJ: Phylogenetic analysis of fungal aquaporins provide insight into their possible role in water transport of mycorrhizal associations. Botany-Botanique. 2013, 91: 495-504. 10.1139/cjb-2013-0041.

    Article  CAS  Google Scholar 

  52. Consortium TU: Activities at the Universal Protein Resource (UniProt). Nucleic Acids Res. 2014, 42: D191-D198. 10.1093/nar/gkt1140.

    Article  Google Scholar 

  53. Kirk PM, Cannon PF, Minter DW, Stalpers JA: Dictionary of the fungi. 2008, CAB International, Wallingford, U. K.

    Google Scholar 

  54. Jones MDM, Forn I, Gadelha C, Egan MJ, Bass D, Massana R, Richards TA: Discovery of novel intermediate forms redefines the fungal tree of life. Nature. 2011, 474: 200-203. 10.1038/nature09984.

    Article  PubMed  CAS  Google Scholar 

  55. Azad AK, Yoshikawa N, Ishikawa T, Sawa Y, Shibata H: Substitution of a single amino acid residue in the aromatic/arginine selectivity filter alters the transport profiles of tonoplast aquaporin homologs. Biochim Biophys Acta. 1818, 2012: 1-11.

    Google Scholar 

  56. Mitani-Ueno N, Yamaji N, Zhao FJ, Ma JF: The aromatic/arginine selectivity filter of NIP aquaporins plays a critical role in substrate selectivity for silicon, boron, and arsenic. J Exp Botany. 2011, 62: 4391-4398. 10.1093/jxb/err158.

    Article  CAS  Google Scholar 

  57. Beitz E, Wu B, Holm LM, Schultz JE, Zeuthen T: Point mutations in the aromatic/arginine region in aquaporin 1 allow passage of urea, glycerol, ammonia, and protons. Proc Natl Acad Sci U S A. 2006, 103: 269-274. 10.1073/pnas.0507225103.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  58. Crooks GE, Hon G, Chandonia JM, Brenner SE: WebLogo: A sequence logo generator. Genome Res. 2004, 14: 1186-1190. 10.1101/gr.849004.

    Article  Google Scholar 

  59. 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:Art. No. 27.,

  60. Verma RK, Jain A, Sankararamakrishnan R: Distinguishing features of aquaglyceroporin in Plasmodium falciparum: Comparative molecular dynamics simulations of three aquaporins.Biophys J 2012, 102:452A.,

  61. Li H, Chen HN, Steinbronn C, Wu BH, Beitz E, Zeuthen T, Voth GA: Enhancement of proton conductance by mutations of the selectivity filter of aquaporin-1. J Mol Biol. 2011, 407: 607-620. 10.1016/j.jmb.2011.01.036.

    Article  PubMed  CAS  Google Scholar 

  62. Savage DF, O’Connell JD, Miercke LJW, Finer-Moore J, Stroud RM: Structural context shapes the aquaporin selectivity filter. Proc Natl Acad Sci U S A. 2010, 107: 17164-17169. 10.1073/pnas.1009864107.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  63. Bienert GP, Bienert MD, Jahn TP, Boutry M, Chaumont F: Solanaceae XIPs are plasma membrane aquaporins that facilitate the transport of many uncharged substrates. Plant J. 2011, 66: 306-317. 10.1111/j.1365-313X.2011.04496.x.

    Article  PubMed  CAS  Google Scholar 

  64. Hub JS, de Groot BL: Does CO2 permeate through aquaporin-1?. Biophys J. 2006, 91: 842-848. 10.1529/biophysj.106.081406.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  65. Ishikawa F, Suga S, Uemura T, Sato MH, Maeshima M: Novel type aquaporin SIPs are mainly localized to the ER membrane and show cell-specific expression in Arabidopsis thaliana. FEBS Lett. 2005, 579: 5814-5820. 10.1016/j.febslet.2005.09.076.

    Article  PubMed  CAS  Google Scholar 

  66. Graf FE, Ludin P, Wenzler T, Kaiser M, Brun R, Pyana PP, Buscher P, de Koning HP, Horn D, Maser P: Aquaporin 2 mutations in Trypanosoma brucei gambiense field isolates correlate with decreased susceptibility to pentamidine and melarsoprol.PLoS Neg Trop Dis 2013, 7:Art. No. e2475.,

  67. Ohta E, Itoh T, Nemoto T, Kumagai J, Ko SBH, Ishibashi K, Ohno M, Uchida K, Ohta A, Sohara E, Uchida S, Sasaki S, Rai T: Pancreas-specific aquaporin 12 null mice showed increased susceptibility to caerulein-induced acute pancreatitis. Am J Physiol Cell Physiol. 2009, 297: C1368-C1378.

    Article  PubMed  CAS  Google Scholar 

  68. Itoh T, Rai T, Kuwahara M, Ko SBH, Uchida S, Sasaki S, Ishibashi K: Identification of a novel aquaporin, AQP12, expressed in pancreatic acinar cells. Biochem Biophys Res Commn. 2005, 330: 832-838. 10.1016/j.bbrc.2005.03.046.

    Article  CAS  Google Scholar 

  69. Nyblom M, Frick A, Wang Y, Ekvall M, Hallgren K, Hedfalk K, Neutze R, Tajkhorshid E, Tornroth-Horsefield S: Structural and functional analysis of SoPIP2;1 mutants adds insight into plant aquaporin gating. J Mol Biol. 2009, 387: 653-668. 10.1016/j.jmb.2009.01.065.

    Article  PubMed  CAS  Google Scholar 

  70. Yukutake Y, Tsuji S, Hirano Y, Adachi T, Takahashi T, Fujihara K, Agre P, Yasui M, Suematsu M: Mercury chloride decreases the water permeability of aquaporin-4-reconstituted proteoliposomes. Biol Cell. 2008, 100: 355-363. 10.1042/BC20070132.

    Article  PubMed  CAS  Google Scholar 

  71. Yool AJ, Morelle J, Cnops Y, Verbavatz JM, Campbell EM, Beckett EAH, Booker GW, Flynn G, Devuyst O: AqF206 is a pharmacologic agonist of the water channel aquaporin-1. J Am Soc Nephrol. 2013, 24: 1045-1052. 10.1681/ASN.2012080869.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  72. Yool AJ, Brown EA, Flynn GA: Roles for novel pharmacological blockers of aquaporins in the treatment of brain oedema and cancer. Clin Exp Pharmacol Physiol. 2010, 37: 403-409. 10.1111/j.1440-1681.2009.05244.x.

    Article  PubMed  CAS  Google Scholar 

  73. Frick A, Jarva M, Tornroth-Horsefield S: Structural basis for pH gating of plant aquaporins. FEBS Lett. 2013, 587: 989-993. 10.1016/j.febslet.2013.02.038.

    Article  PubMed  CAS  Google Scholar 

  74. Yu J, Yool AJ, Schulten K, Tajkhorshid E: Mechanism of gating and ion conductivity of a possible tetrameric pore in aquaporin-1. Structure. 2006, 14: 1411-1423. 10.1016/j.str.2006.07.006.

    Article  PubMed  CAS  Google Scholar 

  75. Yukutake Y, Yasui M: Regulation of water permeability through aquaporin-4. Neuroscience. 2010, 168: 885-891. 10.1016/j.neuroscience.2009.10.029.

    Article  PubMed  CAS  Google Scholar 

  76. Yukutake Y, Hirano Y, Suematsu M, Yasui M: Rapid and reversible inhibition of aquaporin-4 by zinc. Biochemistry. 2009, 48: 12059-12061. 10.1021/bi901762y.

    Article  PubMed  CAS  Google Scholar 

  77. Tornroth-Horsefield S, Wang Y, Hedfalk K, Johanson U, Karlsson M, Tajkhorshid E, Neutze R, Kjellbom P: Structural mechanism of plant aquaporin gating. Nature. 2006, 439: 688-694. 10.1038/nature04316.

    Article  PubMed  Google Scholar 

  78. Didier ES: Microsporidiosis: An emerging and opportunistic infection in humans and animals. Acta Trop. 2005, 94: 61-76. 10.1016/j.actatropica.2005.01.010.

    Article  PubMed  CAS  Google Scholar 

  79. Oide S, Moeder W, Krasnoff S, Gibson D, Haas H, Yoshioka K, Turgeon BG: NPS6, encoding a nonribosomal pepide synthetase involved in siderophore-mediated iron metabolism, is a conserved virulence determinant of plant pathogenic ascomycetes. Plant Cell. 2006, 18: 2836-2853. 10.1105/tpc.106.045633.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  80. Altschul SF, Madden TL, Schaffer AA, Zhang JH, Zhang Z, Miller W, Lipman DJ: Gapped BLAST and PSI-BLAST: a new generation of protein database search programs. Nucleic Acids Res. 1997, 25: 3389-3402. 10.1093/nar/25.17.3389.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  81. Li WZ, Godzik A: Cd-hit: a fast program for clustering and comparing large sets of protein or nucleotide sequences. Bioinformatics. 2006, 22: 1658-1659. 10.1093/bioinformatics/btl158.

    Article  PubMed  CAS  Google Scholar 

  82. Huang Y, Niu BF, Gao Y, Fu LM, Li WZ: CD-HIT suite: a web server for clustering and comparing biological sequences. Bioinformatics. 2010, 26: 680-682. 10.1093/bioinformatics/btq003.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  83. Simossis VA, Heringa J: PRALINE: a multiple sequence alignment toolbox that integrates homology-extended and secondary structure information. Nucleic Acids Res. 2005, 33: W289-W294. 10.1093/nar/gki390.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  84. Krogh A, Larsson B, von Heijne G, Sonnhammer ELL: Predicting transmembrane protein topology with a hidden Markov model: Application to complete genomes. J Mol Biol. 2001, 305: 567-580. 10.1006/jmbi.2000.4315.

    Article  PubMed  CAS  Google Scholar 

  85. McGuffin LJ, Bryson K, Jones DT: The PSIPRED protein structure prediction server. Bioinformatics. 2000, 16: 404-405. 10.1093/bioinformatics/16.4.404.

    Article  PubMed  CAS  Google Scholar 

  86. Tamura K, Peterson D, Peterson N, Stecher G, Nel M, Kumar S: MEGA5: Molecular evolutionary genetics analysis using maximum likelihood, evolutionary distance, and maximum parsimony methods. Mol Biol Evol. 2011, 28: 2731-2739. 10.1093/molbev/msr121.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  87. Chenna R, Sugawara H, Koike T, Lopez R, Gibson TJ, Higgins DG, Thompson JD: Multiple sequence alignment with the Clustal series of programs. Nucleic Acids Res. 2003, 31: 3497-3500. 10.1093/nar/gkg500.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  88. Sievers F, Wilm A, Dineen D, Gibson TJ, Karplus K, Li W, Lopez R, McWilliam H, Remmert M, Soding J, Thompson JD, Higgins DG: Fast, scalable generation of high-quality protein multiple sequence alignments using Clustal Omega.Mol Systems Biol 2011, 7:Art. No. 539.,

  89. Edgar RC: MUSCLE: a multiple sequence alignment method with reduced time and space complexity.BMC Bioinformatics 2004, 5:Art. No. 113.,

  90. Sali A, Blundell TL: Comparative protein modelling by satisfaction of spatial restraints. J Mol Biol. 1993, 234: 779-815. 10.1006/jmbi.1993.1626.

    Article  PubMed  CAS  Google Scholar 

  91. Eswar N, Eramian D, Webb B, Shen M-Y, Sali A: Protein structure modeling with MODELLER. Methods Mol Biol. 2008, 426: 145-159. 10.1007/978-1-60327-058-8_8.

    Article  PubMed  CAS  Google Scholar 

  92. Sui H, Han BG, Lee JK, Walian P, Jap BK: Structural basis of water-specific transport through the AQP1 water channel. Nature. 2001, 414: 872-878. 10.1038/414872a.

    Article  PubMed  CAS  Google Scholar 

  93. Fu D, Libson A, Miercke LJ, Weitzman C, Nollert P, Krucinski J, Stroud RM: Structure of a glycerol-conducting channel and the basis for its selectivity. Science. 2000, 290: 481-486. 10.1126/science.290.5491.481.

    Article  PubMed  CAS  Google Scholar 

  94. Lee JK, Kozono D, Remis J, Kitagawa Y, Agre P, Stroud RM: Structural basis for conductance by the archaeal aquaporin AqpM at 1.68 A. Proc Natl Acad Sci U S A. 2005, 102: 18932-18937. 10.1073/pnas.0509469102.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  95. Canutescu AA, Shelenkov AA, Dunbrack RL: A graph-theory algorithm for rapid protein side-chain prediction. Proc Natl Acad Sci U S A. 2003, 12: 2001-2014.

    CAS  Google Scholar 

  96. Hess B, Kutzner C, van der Spoel D, Lindahl E: GROMACS 4: Algorithms for highly efficient, load-balanced, and scalable molecular simulation. J Chem Theoy Comput. 2008, 4: 435-447.

    Article  CAS  Google Scholar 

  97. Laskowski RA, MacArthur MW, Moss DS, Thornton JM: PROCHECK: A program to check the stereochemical quality of protein structures. J Appl Cryst. 1993, 26: 283-291. 10.1107/S0021889892009944.

    Article  CAS  Google Scholar 

Download references


We gratefully acknowledge the support received from the Centre of Excellence for Chemical Biology at IIT-Kanpur funded by MHRD, Government of India. RS is holding USV Chair Professorship at IIT-Kanpur. RKV thanks CSIR, Government of India for a Senior Research Fellowship. NDP thanks IIT-Kanpur for a fellowship. We acknowledge Dr. Anjali Bansal Gupta and Mr. Manu Vajpai for their help. We thank our lab members for useful discussions.

Author information

Authors and Affiliations


Corresponding author

Correspondence to Ramasubbu Sankararamakrishnan.

Additional information

Competing interests

The authors declare that they have no competing interests.

Authors’ contributions

RS conceived the study. RKV carried out the sequence search, phylogeneic analysis and homology modeling on fungal channels. NDP and RKV performed the studies on plant channels. RKV, NDP and RS analyzed the results. RS wrote the manuscript. All authors read and approved the final manuscript.

Electronic supplementary material

Additional file 1: Table S1.: Contains accession IDs of all fungal MIPs belonging to different subgroups. (DOC 386 KB)


Additional file 2: Table S2.: Contains inter- and intra-group average pairwise sequence identities and similarities of all fungal MIP families. (DOC 42 KB)


Additional file 3: Table S3.: Shows group conservation of small and weakly polar residues at the helix-helix interface calculated for all fungal MIP subgroups. (DOC 54 KB)


Additional file 4: Figure S1.: Relationship between isoelectric point and molecular weight for all the fungal AQGP groups. (TIFF 837 KB)


Additional file 5: Figure S2.: Contains sequence logos produced for Loop B region for different fungal MIP groups. (TIFF 2 MB)


Additional file 6: Table S4.: Contains list of pathogenic fungi that have at least one SIP-like MIP channel. (DOC 114 KB)


Additional file 7: Table S5.: Contains list of plant pathogenic fungi that have at least one member from the δ subgroup of AQGPs. (DOC 97 KB)

Authors’ original submitted files for images

Rights and permissions

Open Access  This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made.

The images or other third party material in this article are included in the article’s Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder.

To view a copy of this licence, visit

The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Verma, R.K., Prabh, N.D. & Sankararamakrishnan, R. New subfamilies of major intrinsic proteins in fungi suggest novel transport properties in fungal channels: implications for the host-fungal interactions. BMC Evol Biol 14, 173 (2014).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: