Research article | Open | Published:
A new subfamily LIP of the major intrinsic proteins
BMC Genomicsvolume 15, Article number: 173 (2014)
Proteins of the major intrinsic protein (MIP) family, or aquaporins, have been detected in almost all organisms. These proteins are important in cells and organisms because they allow for passive transmembrane transport of water and other small, uncharged polar molecules.
We compared the predicted amino acid sequences of 20 MIPs from several algae species of the phylum Heterokontophyta (Kingdom Chromista) with the sequences of MIPs from other organisms. Multiple sequence alignments revealed motifs that were homologous to functionally important NPA motifs and the so-called ar/R-selective filter of glyceroporins and aquaporins. The MIP sequences of the studied chromists fell into several clusters that belonged to different groups of MIPs from a wide variety of organisms from different Kingdoms. Two of these proteins belong to Plasma membrane intrinsic proteins (PIPs), four of them belong to GlpF-like intrinsic proteins (GIPs), and one of them belongs to a specific MIPE subfamily from green algae. Three proteins belong to the unclassified MIPs, two of which are of bacterial origin. Eight of the studied MIPs contain an NPM-motif in place of the second conserved NPA-motif typical of the majority of MIPs. The MIPs of heterokonts within all detected clusters can differ from other MIPs in the same cluster regarding the structure of the ar/R-selective filter and other generally conserved motifs.
We proposed placing nine MIPs from heterokonts into a new group, which we have named the LIPs (large intrinsic proteins). The possible substrate specificities of the studied MIPs are discussed.
The major intrinsic proteins (MIPs) , or aquaporins , allow for the passive transmembrane transport of water and other small, uncharged polar molecules . Glyceroporin (GlpF) from Escherichia coli and aquaporin 1 (AQP1) from bovine  were the first MIPs for which the 3D structures were established through X-ray crystallographic analysis. The similarities between proteins of the MIP family suggest that they have a common origin . Plant aquaporins comprise a large protein family [7–9]. The topology of MIPs resembles a sandwich, consisting of six transmembrane α-helical strands (denoted H1 through H6). These strands are connected to each other by five loops (denoted LA through LE). The LB and LE loops each consist of a short α-helix connected by highly conserved NPA-motifs, and these loops are partly located within the membrane . Certain amino acid residues (a.a.) of the H2 and H5 strands, together with two a.a. of the LE site in the same plane, form the so-called ar/R-filter (aromatic/arginine), which determines the substrate specificity of the protein [4, 5, 11]. Certain sub-families of MIPs contain conserved a.a. within the ar/R-filter [12, 13]: e.g. F56, H180, C189, and R195 in the aquaporin HsAQP1 ; and W48, G191, F200, and R205 in the glyceroporin EcGlpF . The pore diameter of the latter is larger than the former . It was recently shown that substitution of a.a. within the ar/R-filter results in a change in substrate specificity or in a loss of function . However, the design of the ar/R filter is not the only determinant of specificity. Determinants other than the ar/R filter before experimental studies cannot be identified by theoretical analysis for a majority of MIPs which are considered in the present study.
Diatoms are unicellular, phototrophic, eukaryotic organisms that are present in all marine and freshwater habitats. They originated as a result of double endosymbiosis followed by long-term (240 million years) evolution [16–18], which resulted in the migration of many genes, such as bacterial genes, into the diatom nuclear genome [19, 20]. A general feature of diatoms is the presence of an intricately ornamented cell wall, known as a frustule, which consists of silica. The synthesis of the solid and nearly anhydrous elements of the frustule takes place within specialised sub-cellular vesicles (silica deposition vesicules, SDVs) [21, 22]. Maturation of the frustule requires the removal of water from the SDV. It has been proposed that this process is mediated by aquaporins .
In this study we investigate MIPs from Chromista, phylum (Infrakindom) Heterokontophyta (Additional file 1). It was proposed that all pigmented heterokonts appeared due to double symbiosis and the simultaneous appearance of the ability to build cell walls with silica . We used taxonomy from algae base .
Ten MIP genes were found in the complete genome sequences for the diatoms Thalassiosira pseudonana TpMIP1, TpMIP2  and T. oceanica ToMIP1, ToMIP2  (class Coscinodiscophyceae), Phaeodactylum tricornutum PtMIP1, PtMIP2, PtMIP3, PtMIP4, PtMIP5  (class Bacillariophyceae), and Nannochloropsis gaditana NgMIP  (class Eustigmatophyceae). A MIP gene was also recently discovered in the genome of the freshwater araphid pennate diatom Synedra acus subsp. radians (class Fragillaryophyceae). The length of the predicted SarMIP is 286 a.a. . MIP genes were found in the genomes of the diatoms Pseudo-nitzschia multiseries and Fragilariopsis cylindrus (class Bacillariophyceae).
The database dedicated to Major Intrinsic Proteins MIPdb  contains seven MIP sequences from Ectocarpus siliculosus, EsAQP, EsPIP, and EsMIP (class Phaeophyceae), as well as from Aureococcus anophagefferens, AaMIP1, AaMIP2, AaMIP3, and AaMIP4 (class Pelagophyceae).
The purpose of the present study was to compare the predicted a.a. sequences of 20 MIPs from these algae with MIP sequences from a wide variety of organisms. We show that MIPs from heterokonts belong to different subfamilies, and nine of them merge into a new Large intrinsic protein (LIP) subfamily, which is closely related with the SIP subfamily  and the MIPC subfamily .
Search for homologues in the MIPdb
We used the MIPdb to find homologues of the 20 MIPs from heterokonts. This database contains 8429 MIPs belonging to 11 groups (subfamilies): 577 AQPes, 1150 AQPps, 363 GLAes, 1827 GLAps, 1052 GLPps, 192 NIPs, 661 PIPs, 42 SIPs, and 375 TIPs, as well as 1053 predicted MIPs. There are 1137 sequences that are unclassified MIPs. MIPdb contains 16 MIPs from heterokonts, 15 of them belonging to the unclassified group and one of them (EsPIP) belonging to AQPe (Additional file 1).
To find closely related sequences, we used the phmmer procedure in the HMMER3-package. We selected the proteins with the smallest e-values from each of the MIPdb groups. These proteins are hypothesised to be related to the 20 MIPs studied.
To classify the MIPs from heterokonts, their a.a. sequences were subjected to a phylogenetic analysis together with proteins selected based on the results of a search for homologues and reference sequences, which included a typical aquaporin (HsAQP1), a glyceroporin (EcGlpF), and NIPs from rice and maize (OsNIP21, OsNIP22, ZmNIP22, and ZmNIP23). It has been showed that these four NIPs are transporters of silicic acid [34, 35]. Among the 22 analysed MIPs from green algae, we identified seven different groups, including PIPs, GlpF-like intrinsic proteins, and MIPs unique to green algae (MIPA to MIPE) . These proteins have also been analysed in order to compare the results. The total number of sequences subjected to phylogenetic analysis was 212. Alignment was performed using the profile of the MIP PF00230 in the Pfam database.
The results of the phylogenetic analysis (Figure 1) demonstrated that the topology of the tree constructed from the selected sequences corresponds to the modern understanding of the phylogeny of MIPs . MIPs from heterokonts fall into different clades.
Three sequences (PtMIP1, PtMIP2, and PtMIP3) from P. tricornutum, and AaMIP2 from the A. anophagefferens clustered within a large clade that includes GLAp, GLAe, and GLPp with a bootstrap support of 58%.
The sequences of TpMIP1 from T. pseudonana, ToMIP2 from T. oceanic, and AaMIP3 from A. anophagefferens clustered with unclassified MIPs from bacteria with bootstrap supports of 52% and 100%, respectively.
Two sequences (EsAQP and EsPIP) from the brown alga E. siliculosus clustered with plant PIPs with a bootstrap support of 79%. Nine sequences (EsMIP, NgMIP, AaMIP1, PtMIP5, PnmMIP, FcMIP, SarMIP, TpMIP2, and ToMIP1) constituted a separate clade with a bootstrap support of 98%. Their closest relatives are the plant SIPs and MIPC from green algae. AaMIP4 from A. anophagefferens clustered with MIPE from green algae with a bootstrap support of 71%. Only one sequence of PtMIP4 did not reliably belong to any clade. EsAQP and EsPIP clustered with CcPIP4;1 and CcPIP4;2 from green algae. AaMIP2, PtMIP1, PtMIP2, and PtMIP3 clustered with GIPs from green algae as a sister clade with a bootstrap support of 58%.
Peculiarities of the structures of MIPs from heterokonts
Of the 20 heterokonts’ MIPs studied (Figure 2), only eight contain a pair of conserved NPA motifs. It worth noting that the EsMIP, NgMIP, PtMIP5, PnmMIP, FcMIP, SarMIP, TpMIP2, and ToMIP1 proteins, which are the closest to SIPs, all have an NPM motif in place of the second NPA. However, the second motif of SIPs is an NPA and the first motif is variable in the third a.a. position. AaMIP4 has NGA instead the first NPA motif. The MIPC sequences contain the first modified motif NP[T/V].
The amino acids belonging to the ar/R filters are shown in Additional file 1 and Table 1 and Figure 2. MIPs 1, 2, and 3 from P. tricornutum have a.a. compositions in the ar/R filter that are similar to the GLAp MIP, whereas the ar/R filter of AaMIP2, which was in the same clade as GLAp, GLAe, and GLPp, has the same a.a. composition as human aquaporin HsAQP1 at positions H2, LE1, and LE2.
The amino acid composition of the ar/R filters of PtMIP4, AaMIP4, and EsPIP was identical to that of the ar/R filter of human aquaporin HsAQP1. The ar/R filter of EsAQP differs from the ar/R filter of PtMIP4, AaMIP4, and EsPIP by the presence of an alanine (A) in place of the cysteine (C) at position LE1. The LE1 position in the filters of EsMIP, NgMIP, PtMIP5, PnmMIP, FcMIP, AaMIP1, SarMIP, TpMIP2, and ToMIP1, as well as a majority of the filters of the sister clade, was occupied by a proline (P). Unlike the other MIPs, their LE2 positions are occupied by leucine (L) or isoleucine (I). Hence, EsMIP, NgMIP, PtMIP5, PnmMIP, FcMIP, AaMIP1, SarMIP, TpMIP2, and ToMIP1 differed from the other MIPs in this feature. The H2 strand of the filter includes tryptophan (W)/tyrosine (Y) in these nine proteins, which is a more typical feature of GLPp, GLAp, and NIPs that are not able to transport silicic acid (Additional file 1).
Phylogenetic analysis of MIP a.a. sequences predicted from the nucleotide sequences of the respective genes, as well as a comparison of the ar/R filters of 20 MIPs from heterokonts, revealed that some of these proteins have very close homologues among the 8429 proteins documented in the MIPdb. On the basis of the ar/R filter, substrate specificity could be suggested.
Three proteins (PtMIP1, PtMIP2, and PtMIP3) from the diatom P. tricornutum were classified as GlpF-like intrinsic proteins. Their closest homologues are DaMIP, which is from the sulphate reducing anaerobic proteobacteria Desulfuromonas acetoxidans, and GhMIP, which is from the Gram-positive coccus Gemella haemolysans. The ar/R filters of PtMIP1, PtMIP2, and PtMIP3 are identical to those of NIPs of subgroup I, which have glycerol permease activity (Table 1) [12, 37].
AaMIP2 from A. anophagefferens also belongs to the GlpF-like intrinsic proteins, as revealed by the phylogenetic analysis (Figure 1). However, three a.a. in the ar/R filter in the H2, LE1, and LE2 of AaMIP2 are identical to those found in the human aquaporin HsAQP1.
Two proteins, EsAQP and EsPIP, of the E. siliculosus belong to the PIP subfamily, as revealed by the phylogenetic analysis (Figure 1). Unlike EsAQP and EsPIP, all other PIPs are highly conserved [38, 33]. The composition of the ar/R filter of EsAQP and EsPIP is different from that of the ar/R filter of PIPs at some positions (Additional file 1). However, the composition of the ar/R filters of EsPIP and HsAQP1 is identical. We propose that, based on their sequence composition, these two proteins from brown alga E. siliculosus are intermediate forms between the human HsAQP1 and plant PIPs.
AaMIP4 from A. anophagefferens clusters with a specific subfamily MIPE from green algae on the phylogenetic tree. The ar/R filters of AaMIP4 and MIPEs are identical to those of HsAQP1. This similarity suggests these proteins have specificities for water.
A phylogenetic analysis revealed that PtMIP4 from the diatom P. tricornutum does not have close relatives. The sequence of this protein differs from classical aquaporins in that the first NPA motif in PtMIP4 is transformed into NPG, but all residues of the ar/R filter are identical to those of the HsAQP1 and MIPEs. TpMIP1 and ToMIP2 from diatoms and AaMIP3 from A. anophagefferens cluster with bacterial MIPs on the phylogenetic tree. Multiple alignments have shown that the ar/R filters of bacterial MIPs, TpMIP1, ToMIP2, and AaMIP3 differ from those found in the proteins from other subfamilies (Additional file 1). The functions of these bacterial MIPs are not yet known. Therefore, no function can be proposed at this time for TpMIP1, ToMIP2, or AaMIP3. Interestingly, the aquaporin TcAQPe of the parasitic trypanosome Trypanosoma cruzi falls into the same clade, although its ar/R filter  is different from that of typical aquaporins at all four positions.
Nine proteins (EsMIP, NgMIP, PtMIP5, PnmMIP, FcMIP, AaMIP1, SarMIP, TpMIP2, and ToMIP1) form a separate clade adjacent to the SIP clade of plants and the MIPC clade of green algae. Ishibashi et al.  concluded that during their evolution, SIPs and XIPs lost conservation of the NPA motifs. The first motif that replaced NPA in SIPs was NP[T/L/S/I], while in MIPC it was NP[T/V]. The motif that is found in place of the second NPA in these nine proteins is NPM. We found four other proteins with NPM in place of the second NPA in the MIPdb. One of these proteins belongs to the NIP subfamily, and the other three belong to an uncharacterised group. However, arginine (R) is C–terminal to NPM in the LE2 position in these four proteins, whereas leucine (L) or isoleucine (I) are C-terminal to NPM in MIPs from heterokonts. The LE1 position in the ar/R filters of EsMIP, NgMIP, PtMIP5, PnmMIP, FcMIP, AaMIP1, SarMIP, TpMIP2, and ToMIP1 are occupied by the same a.a. as those in these positions in SIPs. However, the sites at positions H2 and LE2 are occupied by different a.a. from those found in the same positions in the SIPs. Remarkably, a tryptophan (W) at position H2 occurs in all glyceroporins, as well as in NIPs that are not able to transport silicic acid (Additional file 1). Residues of the ar/R filter of MIPC do not match any one a.a. position of the ar/R filter of these nine proteins. Similarities were revealed in the terminal a.a. of lysine (K) in SIPs and EsMIP, NgMIP, PtMIP5, PnmMIP, FcMIP, SarMIP, TpMIP2, and ToMIP1 (Figure 2). Of all the MIPCs, only one sequence (MrMIPC1;1) contains a terminal lysine (K).
On the basis of the above evidence, we suggest a new phylogenetic clade, LIPs, which includes nine proteins: EsMIP, NgMIP, PtMIP5, PnmMIP, FcMIP, AaMIP1, SarMIP, TpMIP2, and ToMIP1 (Figure 1). This new clade has high bootstrap support. Indeed, seven of the nine proteins of the subfamily are large (280 to 317 a.a.), with EsMIP consisting of 225 a.a. and NgMIP consisting of 230 a.a. (Additional file 1). However, EsMIP and NgMIP are similar to other LIPs in terms of phylogeny (Figure 1) and the structure of the ar/R filters (Figure 2).
Heterokonts are thought to be derived from a secondary endosymbiotic process between a red alga and a heterotrophic eukaryote . Recent studies of the genomes of the diatoms have revealed the participation of green algae in the origin of some membrane transporters . We showed that of the MIPs from heterokonts, one protein (AaMIP4) has a relationship with a specific subfamily MIPE from green algae.
According to our analysis, none of the 20 analysed MIPs from heterokonts are relatives to NIPs, which transports silicic acid, on the phylogeny (Figure 1), and have dissimilar ar/R filter (Additional file 1 and Table 1).
Heterokonts, like other organisms, contain a variety of MIPs, which could allow for the transport of substances, such as water, glycerol, urea, carbon dioxide, etc. We found that heterokonts contain MIPs that belong to different subfamilies, such as PIP, GIP, and MIPE. The most surprising finding is that during their evolution, heterokonts acquired unique genes, such as those that encode the MIPs of the LIP subfamily. These unusual proteins encoded by these genes are only distantly related to typical aqua- or glyceroporins, and are characterised by a specific motif and the composition of the ar/R filter. Notably, none of MIPs from heterokonts have any similarities with NIPs that are responsible for transporting silicic acid.
Search for homology
A search for closely related MIPs was carried out with MIPdb, which is a motif-oriented database that allows for analyses on the biological, structural, and functional levels and is used to identify highly specific domains of unknown proteins. To analyse the similarities between MIPs from heterokonts and MIP sequences from a wide variety of organisms, we used HMMER3  with the procedure phhmer to carry out a BLAST-like search for a specified sequence in the database.
Alignments and phylogenetic analysis
Multiple sequence alignments of aquaporin amino acids on the Pfam profile of the MIP family PF00230 was carried out using HMMER3 with the procedure hmmalign. The resulting alignment was edited in the JalView program  to remove non-informative C- and N-termini. The phylogenetic trees were constructed using MEGA5 v5.1  using the bootstrap neighbour-joining (NJ) method with 1000 replicates and the Jones-Taylor-Thornton (JTT) model. For tree visualisation, we used iTOL . Multiple sequence alignments of the MIPs of heterokonts were carried out using the program Muscle v3.8.31 .
Availability of supporting data
The data sets supporting the results of this article are available in the DRYAD repository, http://dx.doi.org/10.5061/dryad.8f61v.
Aquaporin of eukaryotes
Aquaporin of prokaryotes
GlpF-like intrinsic protein
Aquaglyceroporin of eukaryotes
Aquaglyceroporin of prokaryotes
Glycerol uptake facilitator of prokaryotes
Glycerol uptake facilitator protein
Major intrinsic protein
Large intrinsic protein
MIP data base
Nodulin 26-like intrinsic protein
Plasma membrane intrinsic protein
Small basic intrinsic protein
Silica deposition vesicle
Tonoplast intrinsic protein
X intrinsic protein.
Gorin MB, Yancey SB, Cline J, Revel J-P, Horwitz J: The major intrinsic protein (MIP) of the bovine lens fiber membrane: characterization and structure based on cDNA clining. Cell. 1984, 39: 49-59. 10.1016/0092-8674(84)90190-9.
Agre P, Sasaki S, Chrispeels MJ: Aquaporins: a family of membrane water channels. Am J Physiol. 1993, 265: F461-
Ishibashi K, Sasaki S, Fushimi K, Uchida S, Kuwahara M, Saito H, Furukawa T, Nakajima K, Yamaguchi Y, Gojobori T: Molecular cloning and expression of a member of the aquaporin family with permeability to glycerol and urea in addition to water expressed at the basolateral membrane of kidney collecting duct cells. Proc Natl Acad Sci USA. 1994, 91: 6269-6273. 10.1073/pnas.91.14.6269.
Fu D, Libson A, Miercke LJW, Weitzman C, Nollert P, Krucinski K, Stroud R: Structure of a glycerol conducting channel and the basis for its selectivity. Science. 2000, 290: 481-486. 10.1126/science.290.5491.481.
Sui H, Han B-G, 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.
Reizer J, Reizer A, Saier MH: The MIP family of integral membrane channel proteins: sequence comparisons, evolutionary relationships, reconstructed pathway of evolution, and proposed functional differentiation of the two repeated halves of the proteins. Crit Rev Biochem Mol Biol. 1993, 28: 235-257. 10.3109/10409239309086796.
Johanson U, Karlsson M, Johansson I, Gustavsson S, Sjövall 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.
Quigley F, Rosenberg JM, Shachar-Hill Y, Bohnert HJ: From genome to function: the Arabidopsis aquaporins. Genome Biol. 2001, 3 (1): research0001.1–0001.17-
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. 10.1093/pcp/pci172.
Tajkhorshid E, Nollert P, Jensen M, Miercke L, O’Connell J, Stroud RM, Schulten K: Control of the selectivity of the aquaporin water channel family by global orientation tuning. Science. 2002, 296: 525-530. 10.1126/science.1067778.
Thomas D, Bron P, Ranchy G, Duchesne L, Cavalier A, Rolland JP, Raguénès-Nicol C, Hubert JF, Haase W, Delamarche C: Aquaglyceroporins, one channel for two molecules. Biochim Biophys Acta. 2002, 1555: 181-186. 10.1016/S0005-2728(02)00275-X.
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-1068. 10.1104/pp.103.033415.
Hub JS, de Groot BL: Mechanism of selectivity in aquaporins and aquaglyceroporins. Proc Natl Acad Sci USA. 2007, 105 (4): 1198-1203.
Wspalz T, Fujiyoshi Y, Engel A: The AQP Structure and Functional Implications. Handbool of Experimental Pharmacology. Aquaporins. Edited by: Beitz E. 2009, Berlin, Heidelberg: Springer-Verlag, 31-56.
Mitani-Ueno N, Yamaji N, Zhao FJ, Ma J: The aromatic/arginine selectivity filter of NIP aquaporins plays a critical role in substrate selectivity for silicon, boron, and arsenic. J Exp Bot. 2011, 62 (12): 4391-4398. 10.1093/jxb/err158.
Kooistra WHCF, Medlin LK: Evolution of the diatoms (Bacillariophyta). IV. A reconstruction of their age from small subunit rRNA coding regions and the fossil record. Mol Phyl Evol. 1996, 6: 391-407. 10.1006/mpev.1996.0088.
Medlin LK, Kooistra WHCF, Potter D, Saunders GW, Anderson RA: Phylogenetic relationships of the “golden algae” (haptophytes, heterokonts, chrysophytes) and their plastids. Plant Syst Evol. 1997, 11 (Suppl): 187-210.
Medlin LK, Kaczmarska I: Evolution of the diatoms: V. Morphological and cytological support for the major clades and a taxonomic revision. Phycologia. 2004, 43 (3): 245-270. 10.2216/i0031-8884-43-3-245.1.
Armbrust EV, Berges JA, Bowler C, Green BR, Martinez D, Putnam NH, Zhou S, Allen AE, Apt KE, Bechner M, Brzezinski MA, Chaal BK, Chiovitti A, Davis AK, Demarest MS, Detter JC, Glavina T, Goodstein D, Hadi MZ, Hellsten U, Hildebrand M, Jenkins BD, Jurka J, Kapitonov VV, Kröger N, Lau WW, Lane TW, Larimer FW, Lippmeier JC, Lucas S, et al: The Genome of the diatom Thalassiosira pseudonana: ecology, evolution, and metabolism. Science. 2004, 306: 79-86. 10.1126/science.1101156.
Bowler C, Allen AE, Badger JH, Grimwood J, Jabbari K, Kuo A, Maheswari U, Martens C, Maumus F, Otillar RP, Rayko E, Salamov A, Vandepoele K, Beszteri B, Gruber A, Heijde M, Katinka M, Mock T, Valentin K, Verret F, Berges JA, Brownlee C, Cadoret JP, Chiovitti A, Choi CJ, Coesel S, De Martino A, Detter JC, Durkin C, Falciatore A, et al: The Phaeodactylum genome reveals the evolutionary history of diatom genomes. Nature. 2008, 456: 239-244. 10.1038/nature07410.
Drum RW, Pankratz HS: Post mitotic fine structure of Gomphinema parvulum. J Ultrastruct Res. 1964, 10: 217-223. 10.1016/S0022-5320(64)80006-X.
Reimann BEF: Deposition of silica inside a diatom cell. Exp Cell Res. 1964, 34: 605-608. 10.1016/0014-4827(64)90248-4.
Grachev MA, Annenkov VV, Likhoshway YV: Silicon nanotechnologies of pigmented heterokonts. BioEssays. 2008, 30: 328-337. 10.1002/bies.20731.
Medlin LK, Koositra WHCF, Schmid A-MM: A review of the evolution of the diatoms: a total approach using molecules, morphology an geology. The Origin and Early Evolution of the Diatoms: Fossil, Molecular and Biogeographical Approaches. Edited by: Witkowski A, Sieminska J. 1999, Cracow: W. Szafer Institute of Botany, Polish Academy of Sciences, 13-34.
Global species database of information on all groups of algae. http://www.algaebase.org,
Lommer M, Specht M, Roy A-S, Kraemer L, Andreson R, Gutowska MA, Wolf J, Bergner SV, Schilhabel MB, Klostermeier UC, Beiko RG, Rosenstiel P, Hippler M, LaRoche J: Genome and low-iron response of an oceanic diatom adapted to chronic iron limitation. Genome Biol. 2012, 12 (7): R66-
Radakovits R, Jinkerson R, Fuerstenberg SI, Tae H, Settlage RE, Boore JL, Posewitz MC: Draft genome sequence and genetic transformation of the oleaginous alga Nannochloropis gaditana. Nat. Commun. 2012, doi:10.1038/ncomms1688
Petrova DP, Khabudaev KV, Marchenkov AM, Galachyants YP, Kaluyzhnaya OV, Zakharova YR, Likhoshway YV, Grachev MA: The aquaporin-like protein of the diatom Synedra acus. Dokl Biochem Biophys. 2013, 448 (2): 1-4.
JGI P. multiseries v1. http://genome.jgi.doe.gov/Psemu1/Psemu1.home.html,
JGI F. cylindrus v1. http://genome.jgi-psf.org/Fracy1/Fracy1.home.html,
A database dedicated to major intrinsic proteins. http://mipdb.genouest.org,
Johanson U, Gustavsson S: A new subfamily of major intrinsic proteins in plants. Mol Evol Biol. 2002, 19 (4): 456-461. 10.1093/oxfordjournals.molbev.a004101.
Anderberg HI, Danielson JA, Johansos U: Algal MIPs, high diversity and conserved motifs. BMC Evol Biol. 2011, 11: 110-10.1186/1471-2148-11-110.
Ma JF, Tanai K, Yamaji N, Mutani N, Konishi S, Katsuhara M, Ishiguro M, Murata Y, Yano M: A silicon transporter in rice. Nature. 2006, 440 (30): 688-691.
Mitani N, Yamaji N, Ma JF: Identification of maize silicon influx transporters. Plant Cell Physiol. 2008, 50 (1): 5-12.
Danielson JH, Johanson U: MIPs and their role in the exchange of metalloids. Advances in Experimental Medicine and Biology. Edited by: Thomas PJ, Gerd PB. 2010, New-York: Springer Science & Business Media, LLC dual imprint, 19-27.
Weig AR, Jakob C: Functional identification of the glycerol permease activity of Arabidopsis thaliana NLM1 and NLM2 proteins by heterologous expression in Saccharomyces cerevisiae. FEBS Lett. 2000, 481: 293-298. 10.1016/S0014-5793(00)02027-5.
Zardoya R: Phylogeny and evolution of the major intrinsic protein family. Biol Cell. 2005, 97: 397-414. 10.1042/BC20040134.
Ishibashi K, Kondo S, Hara S, Morishita Y: The evolutionary aspects of aquaporin family. Am J Physiol Regul Integr Comp Physiol. 2011, 300: 566-576. 10.1152/ajpregu.90464.2008.
Chan CX, Reyes-Prieto A, Bhattacharya D: Red and green algal origin of diatom membrane transporters: insights into environmental adaptation and cell evolution. PLoS One. 2011, 6 (12): e29138-10.1371/journal.pone.0029138. doi: 10.1371/journal.pone.0029138
Wistrand M, Sonnhammer EL: Improved profile HMM performance by assessment of critical algorithmic features in SAM and HMMER. BMC Bioinforma. 2005, 6: 99-10.1186/1471-2105-6-99. doi: 10.1186/1471-2105-6-99
Waterhouse A, Procter J, Martin D, Clamp M, Barton G: Jalview Version 2 - a multiple sequence alignment editor and analysis workbench. Bioinformatics. 2009, 25 (9): 1189-1191. 10.1093/bioinformatics/btp033.
Tamura K, Peterson D, Peterson N, Stecher G, Nei M, Kumar S: MEGA5: molecular evolutionary genetics analysis using maximum likelihood, evolutionary distance, and maximum parsimony methods. Mol Biol Evol. 2001, 28: 2731-2739.
Letunic I, Bork P: Interactive Tree of Life (iTOL): an online tool for phylogenetic tree display and annotation. Bioinformatics. 2007, 23: 127-128. 10.1093/bioinformatics/btl529.
Edgar RC: MUSCLE: multiple sequence alignment with high accuracy and high throughput. Nucleic Acids Res. 2004, 32 (5): 1792-1797. 10.1093/nar/gkh340.
The authors are thankful to Yuri Galachyants for his valuable comments on the manuscript and Evgeniy Pozdnyak for pre-analysis of the data.
This work was supported by the Program of the Russian Academy Presidium “Molecular and Cell Biology” [grant 6.9].
The authors declare that they have no competing interests.
KVK and DPP participated in the sequence alignment and studied the homology of MIPs. KVK selected MIPs most suitable for phylogenetic analysis, and drafted the manuscript. MAG participated in the design of the study and edited the manuscript. YVL conceived of the study, participated in its design, supervised phylogenetic analysis and drafted the manuscript. All authors read and approved the final manuscript.