Introduction
============

Bread wheat, an allohexaploid plant also known as *Triticum aestivum* Linnaeus (L.), is cultivated in the Southern and Northern hemispheres and is one of the most important cereal crops next to rice in the world in terms of planting area. Bread wheat contributes approximately 95% towards total production and it occupies a central position in agricultural policies. Extreme environmental conditions such as rise in temperature and drought are factors that affect its production. [@b1-ebo-9-2013-437] Plants have an adaptive mechanism to regulate water balance in response to a variety of challenging environmental conditions. To withstand the drought stress, it induces a number of proteins, such as *late embryogen esis abundant (LEA), aquaporin (AQP),* and *proline synthetase*.[@b2-ebo-9-2013-437] Aquaporins are \~26--30 kDa water channel proteins, belongs to a highly conserved group of membrane proteins called the major intrinsic proteins (MIPs) that form a large family. MIPs have been reported in animals, insects, fungi, protozoa, bacteria, archaea, and plants and play a key role in transport of water molecules across cell membrane.[@b3-ebo-9-2013-437] The molecular and functional characterization of aquaporins has highlighted the importance of their regulation in response to environmental stimuli.[@b4-ebo-9-2013-437] In plants, aquaporins are most abundantly present in the plasma membrane and in the vacuolar membrane, and possibly in other internal membranes.[@b5-ebo-9-2013-437] Considering the importance of water transport across the membrane, studies on plant AQPs and plant water relations have been well documented. Earlier research has shown that AQPs assist with transportation through the membrane for various physiologically important molecules such as water, glycerol,[@b6-ebo-9-2013-437] CO~2~, H~2~O~2~, boron and silicon,[@b7-ebo-9-2013-437] as well as assist with various physiological processes such as drought avoidance,[@b8-ebo-9-2013-437] salt tolerance, [@b9-ebo-9-2013-437] and the chilling response.[@b10-ebo-9-2013-437] This implies that AQPs have vital roles in transporting a large volume of water with minimal energy expenditure and appear to regulate the trans-cellular route of water.[@b10-ebo-9-2013-437]

On the basis of sequence similarities, plant aquaporins have been classified into 4 subfamilies: the plasma membrane intrinsic proteins (PIPs); tonoplast intrinsic proteins (TIPs); nodulin26-like membrane intrinsic protein (NIPs); small membrane intrinsic proteins (SIPs) and basic membrane intrinsic proteins (BIPs).[@b11-ebo-9-2013-437],[@b12-ebo-9-2013-437] There are a number of studies focused on the *in silico* prediction of aquaporin gene in plants. For instance, 35 aquaporin genes were identified in *Arabidopsis*,[@b13-ebo-9-2013-437],[@b14-ebo-9-2013-437] 33 from *Oryza sativa* L*.,*[@b15-ebo-9-2013-437] 28 from *Vitis vinifera* L.,[@b16-ebo-9-2013-437] and 23 from a moss called *Physcomitrella patens*[@b17-ebo-9-2013-437] using whole genome sequences. Recently, 55 expressed aquaporin genes were detected in *Populus trichocarpa*[@b18-ebo-9-2013-437] genome sequence data while 33 aquaporin genes from *Zea mays*[@b19-ebo-9-2013-437] were determined using expressed sequence tag (EST) data analysis.

Aquaporins exhibit a typically conserved structure with 6 transmembrane helices (TMHs; H1--H6), connected by 5 loops (loops A--E), and two Asn-Pro- Ala amino acid motif. This conserved NPA motif is found to confer substrate selectivity for molecularw transport.[@b13-ebo-9-2013-437] The expression profile of each aquaporin gene family is regulated differentially. Javot and Maurel's[@b20-ebo-9-2013-437] studies revealed that AQPs expression level was highly abundant in roots during soil water uptake.

Promising observations from the evaluation ofw aquaporins among stress resistance or sensitive plants, such as drought susceptible and drought tolerant wheat cultivars,[@b22-ebo-9-2013-437] other crop cultivars,[@b23-ebo-9-2013-437] or stressed EST libraries,[@b24-ebo-9-2013-437] clearly indicate that aquaporins would be important for water uptake, transport and identification or development of any stress tolerant genotypes of crop species. The present work tackles challenges for development of functional marker (FMs) from sequence polymorphisms present in allelic variants of an aquaporin gene. FMs precisely distinguish alleles of a targeted gene, and are modern molecular markers for marker-assisted selection in wheat breeding.[@b25-ebo-9-2013-437] Improvement of productivity of wheat cultivars under drought conditions has become one of the important breeding program objectives in wheat. The performance of genotypes under drought conditions is largely attributed to genetic variations, mostly at the single nucleotide level. Therefore, an *in silico* analysis was conducted with the aim to identify all possible SNPs within aquaporin ESTs and to demonstrate the utility of SNP in genetic mapping and genetic diversity applications.

An enormous deal is known about aquaporin proteins in several plant species; however, little information is available about the aquaporin gene family in wheat. This may be due to the unavailability of its complete genome sequence, as well as the allohexaploid nature of the genome, which is proving to be for gene identification and analysis. Previously, 35 PIP and TIP aquaporin genes have been identified in wheat as reported by Forrest and Bhave,[@b11-ebo-9-2013-437] as well as Yousif and Bhave.[@b12-ebo-9-2013-437] Determining and understanding of the molecular mechanisms underlying the response to abiotic stress responses is required for genetic improvement of wheat for stress tolerance.[@b23-ebo-9-2013-437] Although *AQP* genes respond to multiple stresses, their precise role in abiotic stress tolerance remains unclear. As reported by Forrest and Bhave,[@b11-ebo-9-2013-437] more MIPs may exist in wheat due to presence of 3 homologous genomes in wheat. The aim of this work was to identify and characterize more wheat aquaporin gene encoding AQPs that are potentially important for the regulation of root water flow. In this regard, we developed a bioinformatics pipeline by combining an in-house PERL script with openly available analysis tools to identify new wheat aquaporin (*TaAQP*) genes from a publicly available wheat EST (wEST) database. Here, identification of 13 MIP genes in wheat provides extensive information for functional studies and the development of markers for stress tolerance. Phylogenetic analysis revealed how the evolutionarily-conserved multigene aquaporin family behaves in a polyploid species such as wheat. Conserved consensus motifs exist across the aquaporin subfamilies to support their association. Using the available EST information as a source of expression profile, aquaporin genes from wheat were detected in 9 different tissues. A 3- dimensional (3D) model of the *T. aestivum* aquaporin (PIP1:2) has been developed. This is the first report of the protein structure analysis of *T. aestivum PIP* gene and development of *PIP* gene specific marker. Identification of the functional genes in wheat is becoming an important line of research, which could help us to understand the molecular genetic basis for the wheat genetic improvement and also provide the functional genetic resources for transgenic research.

Materials and Methods
=====================

Database sources
----------------

The *T. aestivum* EST sequences were downloaded from the National Center for Biotechnology Information (NCBI) in the form of UniGene (Ta.seq.all.gz). Rice aquaporins family sequences were retrieved from Uniprot database (<http://www.uniprot.org/>) and used as query. To avoid false positives, 50% identity and e-value of 1e-3 was taken as the threshold for the sequences obtained from BLAST analysis. The transcripts obtained from the search were analyzed using the CAP3 Sequence Assembly Program.[@b26-ebo-9-2013-437]

Bioinformatics analysis
-----------------------

Open reading frames (ORFs) were generated by ORF finder (<http://www.ncbi.nlm.nih.gov/gorf>). Sub-cellular localization prediction of each predicted TaAQP was carried out using WoLF PSORT (<http://wolfpsort.seqcbrc.jp>) and TargetP 1.1 (<http://www.cbs.dtu.dk/services/TargetP>). 2D and 3D structure alignment was performed using MATRAS (<http://strcomp.protein.osaka-u.ac.jp/matras/>). Conserved domains and signature pattern analysis was performed using SMART (<http://smart.embl-heidelberg.de/>) and Scanprosite (<http://prosite.expasy.org/scanprosite/>) programs.

Phylogenetic tree construction
------------------------------

Amino acid sequences of 67 gramineous PIPs, TIPs, SIPs and TIPs, including 20 from rice (*Oryza sativa*), 29 from maize (*Zea mays*), 11 from wheat, and 7 from barley (*Hordeum vulgare*), and the remainder from from dicot species like *Arabidopsis thaliana* and *Brachypodium distachyon* were aligned with 13 *TaAQP* identified in this study using ClustalW program (<http://www.ebi.ac.uk/clustalw/>). To investigate the evolutionary relationship among aquaporins proteins, a phylogenetic tree was constructed by employing the minimal evolution (ME) method and the neighbor-joining (NJ) method using the Bootstrap value of 1,000 replicates wrapped in the MEGA4 (<http://www.megasoftware.net/mega.html>)[@b27-ebo-9-2013-437] software suite.

Analysis of expression profiles
-------------------------------

The expression profile was determined by analyzing the EST counts based on UniGene of *T. aestivum* for the various tissues (<http://www.ncbi.nlm.nih.gov/UniGene>). The EST expression profile was calculated by

Transcript per million

=

Gene EST

/

Σ

Total EST in the pool

*In silico* mining of SNP
-------------------------

Aquaporin ESTs were downloaded from Unigene and were cleaned to remove contaminating sequences. Vector sequences and other contaminations were identified by using the VecScreen web server (<http://www.ncbi.nlm.nih.gov/VecScreen/VecScreen.html>). Poly-A/T tails were completely trimmed by the EST trimmer Perl program. After pre-cleaning, EST sequences shorter than 50 bases were discarded. Furthermore, low complexity regions were masked by using Repeat Masker (<http://www.repeatmasker.org/>). ESTs clustering and SNP identification as performed on the Seqman module of DNAstar software (DNASTAR, USA). A Perl script was written to parse Indel data from output results.

Validation of SNPs
------------------

The allele-specific primers were designed from the position of base transition. The rationale in designing the primers was based on the premise that the 3′-terminal positions ought to be unique among the known wheat genomic sequences. The primers were designed from consensus sequence of contigs which have the candidates SNPs by Primer3 (<http://frodo.wi.mit.edu/cgibin/primer3/primer3_www.cgi>) with Tm 55 °C--65 °C, and PCR product size 175 bp. Polymerase chain reaction (PCR) was performed in a 25 μL volume containing 100 ng of genomic DNA, 2.5 μL of 10 xPCR buffer, 200 μM of each dNTP, 0.2 μM of each primer and 1.0 unit of *Taq* DNA polymerase (Bangalore Genei, India). The thermocycling program consisted of an initial denaturation at 94 °C for 4 min, followed by 30 cycles of 45 sec at 94 °C, 45 sec at 64 °C, 60 sec at 72 °C and a final cycle of 5 min at 72 °C in S-1000 Thermal Cycler (Bio-Rad). The amplified product was separated by electrophoresis on 3% (w/v) agarose gel and stained with ethidium bromide.

3-D model generation and validation
-----------------------------------

The model of *PIP* gene (Accession no AF366564) was built by Modeller 9v8 (<http://www.salilab.org/modeller/>), using a template structure (2B5F: spinach aquaporin) with 74% identity. The model with the lowest discrete optimized protein energy (DOPE) score was selected. Finally, the model was minimized using GROMACS[@b28-ebo-9-2013-437] and its stereochemical quality was evaluated using PROCHECK (<http://www.ebi.ac.uk/thornton-srv/software/PROCHECK/>) RSMD tertiary structures of both template and modeled protein was calculated by SuperPose (<http://wishart.biology.ualberta.ca/SuperPose/>). The final 3D structure of *T. aestivum* aquaporin was submitted to the Protein Model Database (PMDB; <http://mi.caspur.it/PMDB/>).

Results and Discussion
======================

Identification and characterization of aquaporin genes in wheat
---------------------------------------------------------------

ESTs have proven to be an excellent source for gene discovery, molecular marker discovery and gene expression analysis. To ascertain aquaporin genes in the wheat genome, wheat's UniGene database, Ta.seq. all, which contains 56,943 entries, was obtained. Aquaporin families from *O. sativa* were employed against the UniGene database of wheat by BLAST application. As a result 349 selected ESTs were assembled into continuous sequences of length ranged from 292--128 bp, encoding 13 aquaporin genes. The BLAST identity ranged from 51%--90% by E-value of 0 to 8e-10^9^ ([Table 1](#t1-ebo-9-2013-437){ref-type="table"}). The total number of ESTs used for the assembly was highest in TaPIP1-1 ([Table 1](#t1-ebo-9-2013-437){ref-type="table"}). All AQP existence was supported by EST hits from SwissProt or GenBank database ([Table 1](#t1-ebo-9-2013-437){ref-type="table"}). Since there was standard nomenclature established for AQPs in *Arabidopsis* and rice, we named the predicted MIPs accordingly. In this study we noticed that all the predicted AQPs had conserved transmemebrane helix domains, which is the backbone for the AQP family. Based on homologies between DNA sequences, amino acid sequences, structure domain, and the phylogenetic analysis in comparisons with rice, 6 were confirmed as TaPIP subfamily as they exhibited closest sequence identity (61%--90%) to rice PIP cDNAs. Likewise, the others were confirmed as TaNIP based upon grouping and closest identity to rice NIP cDNAs (73%--86%), followed by TaSIP (51%--88% identity) whereas TaTIP showed 86% identity with rice TIP cDNAs (for detail see [Table 1](#t1-ebo-9-2013-437){ref-type="table"}). The identified TaPIP subfamily can be further divided into 2 subfamilies named PIP1 and PIP2 (PIP2-1, PIP2-3, PIP1-3, PIP1-1, PIP2-2, PIP2-4). Similarly TaNIP can be further divided into 3 subgroups TaNIP1, TaNIP2, and TaNIP 4, (NIP4-1, NIP1-3, NIP1-1, NIP2-1) and SIP into 2 subgroups TaSIP1 and TaSIP2 (TaSIP1-2, TaSIP1-1), respectively. In the light of the large numbers of AQPs genes in rice, maize and *Arabidopsis*, it was possible that additional genes may exist in wheat. Our results are addition on the AQP identified by Forrest and Bhave[@b11-ebo-9-2013-437],[@b12-ebo-9-2013-437] as they do not identify SIP and NIP in their investigation under a genome-wide analysis of the aquaporin gene family in wheat. Thus, PIP and TIP genes with a new set of wEST followed by new members of NIP and SIP were identified in the wheat genome, which will add further diversity to the wheat MIP super family and provide valuable information for understanding the whole AQP family of common wheat.

Comparisons of aquaporin genes in selected plant species
--------------------------------------------------------

For comparative analysis, we used aquaporin genes of rice, barley, wheat, *Arabidopsis*, maize and grapes available in the public domain ([Table 2](#t2-ebo-9-2013-437){ref-type="table"}). The genome of *Arabidopsis* encodes the largest number of homologous aquaporin,[@b13-ebo-9-2013-437] followed by *O. sativa* L.[@b15-ebo-9-2013-437] and *Z. mays* L.[@b19-ebo-9-2013-437] It is interesting to note here that wheat, rice and barley belong to the same grass family so their overall aquaporin gene number (identified aquaporin) proportions were very close.

Sequence alignment and conserved motifs in wheat and rice
---------------------------------------------------------

The sequence and structural alignment of deduced amino acid sequences of wheat and some well-described AQPs from selected plant species showed a high degree of homology ([Fig. 1](#f1-ebo-9-2013-437){ref-type="fig"}). Considering the high homology between aquaporin subgroups irrespective of species, these groups may have evolved before the divergence of higher plants. It was observed that all the aquaoprin subfamilies have a tendency to have alpha (α) helical structure followed by a beta (β) sheet. PIP2 genes are highly homologous as compared to PIP1 and other aquaporin subgroups ([Fig. 1B and C](#f1-ebo-9-2013-437){ref-type="fig"}). Functional studies of AQPs showed that PIP2s reveal high osmotic water permeability in contrast with PIP1 members that show lower or no water permeability when expressed in *Xenopus oocytes,* maize,[@b19-ebo-9-2013-437] poplar,[@b29-ebo-9-2013-437] gravepine[@b16-ebo-9-2013-437] and wheat.[@b30-ebo-9-2013-437] A high degree of divergence was observed in NIPs ([Fig. 1B](#f1-ebo-9-2013-437){ref-type="fig"}). Variations were observed between *A. thaliana* and *G. max* in SIP1-1, whereas in SIP1-2 subfamily, no conserved residues were detected ([Fig. 1D](#f1-ebo-9-2013-437){ref-type="fig"}). Because of the high amino acid sequence homology and functional diversity of AQPs in plants, a systematic analysis of signature sequences and residues is required to precisely identify and classify newly predicted AQPs. The translated amino acid sequences of TaAQP were scanned for motifs, signature pattern and conserved secondary structure to obtain a detailed functional analysis.

Computational prediction revealed up to 6 transmembrane helices with 21 amino acid residues in each helix among all 13 deduced amino acid sequences of TaAQPs as shown in [Supplementary Figure 1](#s1-ebo-9-2013-437){ref-type="supplementary-material"}. Predicted amino acid sequences of the aquaporin gene family revealed an NPA motif (Asn-Pro-Ala) that was extremely conserved.[@b11-ebo-9-2013-437],[@b12-ebo-9-2013-437] Two NPA motifs were predicted in PIPs, NIPs, TIPs, and one in SIPs ([Supplementary Fig. 2](#s1-ebo-9-2013-437){ref-type="supplementary-material"}). The NPA boxes and adjacent residues are considered essential for water transport activity.[@b31-ebo-9-2013-437] The first NPA box is located in the first cytosolic loop and the second NPA box is located in the third extracellular (or vacuolar for TIPs) loop. Substrate selectivity for molecular transport is recognized by the presence of an NPA box.[@b32-ebo-9-2013-437] The helical region is of particular importance for AQP structure because they contain the conserved NPA motifs that are functionally important for water channels. The C-terminal region of TaPIPs and NIPs encompasses conserved phosphorylation sites KXSXXR/K ([Supplementary Fig. 2](#s1-ebo-9-2013-437){ref-type="supplementary-material"}). Phosphorylation of the C-terminal serine may regulate aquaporin activity in response to an osmotic signal.[@b18-ebo-9-2013-437] An increase in aquaporin activity following phosphorylation was demonstrated in α-TIP and PM28A in oocytes.[@b33-ebo-9-2013-437] AQPs belonging to PIP1 subfamily encode polypeptides of 244--292 amino acids in length with 92.47% (TaPIP1-1 and TaPIP1-3) shared sequence identity in wheat. The length of the PIP2 (PIP2-3) ORF ranged from 170--320 amino acids, which had a relatively high sequence similarity of 88% with PIP2-2. In addition, PIP2-2 and PIP2-4 shared 80% amino acid sequence identity with each other. In general, lengths of estimated ORFs were quite different between PIP1 and PIP2. The predicted single polypeptide of TIP1 was 252 amino acids in length. However, predicted SIP1-1 polypeptide ORFs ranged from 128 to 244 amino acids that shared 55% sequence identity (see [Supplementary Table 1](#s1-ebo-9-2013-437){ref-type="supplementary-material"}). NIPs were divided phylogenetically into 4 different subgroups (NIP1-2, NIP1-3, NIP2-1 and NIP4-1) with ORF lengths ranging from 245 to 250 amino acids. However, TaNIP2-1 and TaNIP4-1 (26%) showed minimum identity. After detailed fingerprinting analysis of different classes of putative aquaporin, we found that these AQP proteins have a signature motif identical to phosphokinase C, tyrosine kinase and casein kinase II proteins, as well as N-myristoylation site and MIP family signature sequences SGxHxNPA shown in [Table 3](#t3-ebo-9-2013-437){ref-type="table"}. Specific signature pattern scans were identified to be PKC_PHOSHO_SITE, N-Myristoylation site, MIP signature and CK2_PHOSPHO_SITE families of 3, 6, 4 and 9 amino acids in length and present in all 4 classes of aquaporin. These families have a more or less significant relationship with regulation and metabolism. Plant protein kinase C (PKC) cascades are likely to be involved in mitogen- activated signaling pathways, cellular regulation and metabolism, gene expression, metabolism, motility, membrane transport, and apoptosis.[@b34-ebo-9-2013-437] Previous studies have reported that Casein kinase II, a selective protein kinase, is possibly involved in circadian clock regulation, photoperiod sensitivity[@b35-ebo-9-2013-437] and various regulatory processes of plants. In eukaryotes, N-myristoylation (N-MYR) plays an important role in cell physiology, alters the lipophilicity[@b36-ebo-9-2013-437] of the target protein and facilitates its interaction with membranes, thereby affecting its subcellular localization. N-MYR in plant cells has a critical role in controlling membrane-signaling pathways that lead to specific plant immunity.[@b37-ebo-9-2013-437] For instance, in *Arabidopsis,* N-MYR-related protein mediates plant rescue metabolism to survive damage created by external uneven conditions. This provides evidence of structural and functional correlation of the aquaporin protein class in relation to plant metabolism. We found that functional amino acid sequence motifs are well conserved among members of a subgroup and they are likely to have similar functions. Typical features were revealed by amino acids of the predicted aquaporin proteins such as 6 transmembrane α-helices, NPA box, and homologies with other plant proteins, which provided strong indications that the identified, putative wheat AQPs exhibit water channel activity. All these highly-conserved residues have been shown functional importance for substrate filtering and gating of aquaporin channel proteins.[@b38-ebo-9-2013-437]

Sub-cellular location of the predicted AQPs
-------------------------------------------

AQPs were characterized in this study and a link with AQPs subcellular localization was established. Identified TaPIPs, TaNIPs, TaTIPs and TaSIPs were figured out, to contain a signal peptide localized in various places; ie, mitochondria (TaPIP1-1, TaPIP1-3 and TaPIP2-1), tonoplast (TaTIP4-1), plasma membrane (PM) and mitochondria (TaNIP1-1, TaNIP1- 3, TaNIP2-1, TaNIP4-1, TaPIP2-2, TaPIP2-3 and TaPIP2-4), and endoplasmic reticulum (TaSIP1-1, and TaSIP1-1; [Table 1](#t1-ebo-9-2013-437){ref-type="table"}). Proper sub-cellular localization represents a mechanistic regulation of aquaporin activity that possibly relies on its ability to form multimers between members of different subgroups. [@b39-ebo-9-2013-437] In consistent with previous reports, TaPIPs proteins were found in PM and mitochondrial subcellular compartments, which are identical to those predicted in ice plants (*Mesembryanthemum crystallinum*). [@b9-ebo-9-2013-437] However, most of the cotton PIPs were also found to localize in PM, and these would be important for movement of water and other nonpolar small molecules.[@b21-ebo-9-2013-437],[@b40-ebo-9-2013-437] In rice, sub-cellular localization of all PIP members were predicted to exist in the PM in addition to chloroplast and mitochondria or tonoplast. They could also exist in both PM and tonoplast, while TIPs are localized either to PM or tonoplast or both.[@b11-ebo-9-2013-437]

Phylogenetic relationships of aquaporin family genes in different plant species
-------------------------------------------------------------------------------

To understand the evolution and conservation of crop species, elucidation of the evolutionary relationships is a crucial step.[@b41-ebo-9-2013-437] Earlier, Zardoya et al[@b42-ebo-9-2013-437] established a phylogenetic framework for the aquaporin family in eukaryotes. In order to systematically classify the putative TaAQP, a phylogenetic tree was constructed using bootstrap analysis based on multiple sequence alignments of the proteins sequences of *Vitis vinifera*, *Glycine max, B. distortion, Z. mays, O. sativa, A. thaliana* and *H. vulgar.* Phylogenetic analysis confirmed that the identified wheat AQPs could be classified into 4 large and highly similar orthologous groups of aquaporin subfamilies; ie, PIP, TIP, SIP and NIP ([Fig. 2](#f2-ebo-9-2013-437){ref-type="fig"}). The NIP subfamily exhibits 3 distinct clusters corresponding to the NIP1, NIP2, NIP4 subgroups. The PIP subfamily is divided into PIP1 and PIP2 subgroups but has slight pair wise-sequence divergence when compared with NIPs, TIPs and SIPs. This may be due to the slower rate of evolution of PIPs compared to other subgroups. There were 6 AQPs from rice, with *Arabidopsis* and wheat grouped under NIP subfamily, 4 AQPs from wheat, 8 from *Arabidopsis* and 5 from rice grouped into the TIP subfamily. There were only 2 AQPs from wheat, 2 from rice and 3 from *Arabidopsis* were grouped under SIP family, which was confirmed to be a smaller subfamily of AQPs. PIP proved to be largest subfamily among all AQPs, with 10 PIPs from wheat, 10 from rice and 12 from *Arabidopsis,* respectively. PIPs, SIPs, are closely related. As the results indicated, wheat, rice and *Arabidopsis* show similarity in the number of AQPs in each subgroup. All branches were supported by bootstrap value. Branches corresponding to less than 50% bootstrap replicates are collapsed. Phylogenetic analysis performed based upon AQPs homologous revealed that as many as 13 identified AQP genes were clearly divided into 4 different subfamilies (PIP, SIP, NIP and TIP). Thus, the phylogenetic clustering pattern analysis obtained in this study is similar to the findings reported in other plants previously.[@b8-ebo-9-2013-437],[@b15-ebo-9-2013-437]

Expression profile analysis of aquaporin genes
----------------------------------------------

It has been suggested that the aquaporin gene family expression pattern in response to abiotic stress in different plant tissues at various stages has been detected in various plants. To demonstrate the utility of digital expression analysis of numerous genes across various tissues, we performed gene expression analysis of the wEST collection.[@b43-ebo-9-2013-437] The expression profiles varied with aquaporins subfamilies. 5 genes from PIP (PIP1-1, PIP1, PIP2, PIP2-3, and PIP2-7), 6 genes from TIP (TIP2-2, TIP2-1, TIP2-3, TIP1-2, TIP3-1 and TIP4-3), 3 each from NIP (NIP1-1, NIP1-3, and NIP3-2) and AQPs (AQP5, AQP3, AQP4, AQP2, AQP7) were analyzed to compare the abundance of mRNA transcripts in various tissues. The aquaporin family genes from *T. aestivum* were found in the callus, crown, flower, inflorescence, leaf, root, seed, sheath and stem. [Figure 3](#f3-ebo-9-2013-437){ref-type="fig"} shows aquaporin gene transcripts expressed at relatively high levels in crown and root tissue (36%), followed by stem tissue (13%), crown (12%), inflorescences and leaf (8%) flower (10%), and, least abundantly, in cell culture and callus (3%). We performed gene expression analysis in as many as 9 wheat tissues; ie, callus, cell culture, crown, stem, inflorescence, leaf, flower, seed and root ([Fig. 3](#f3-ebo-9-2013-437){ref-type="fig"}). Our results showed that PIP1-2, PIP1, TIP2-1 and APQ3 genes were detected in most of the tissues. It is noteworthy that a relatively very high expression level of TIP2-1 and APQ3 was found in the root and PIP1-2, and PIP1 in crown tissues ([Fig. 4](#f4-ebo-9-2013-437){ref-type="fig"}). However, the level of expression pattern was observed to be low in the majority of tissues examined for PIP2-3, PIP2-7, NIP1-3, NIP3-2, TIP1-2 and TIP4-2. In order to determine the functionality of this new MIP subfamily, advance studies are required to investigate expression pattern, localization and substrate specificity. Understanding expression profiling of aquaporin genes will help to develop *T. aestivum* cultivars that may be adaptable to drought stress conditions and to support their potential role in the development of the respective tissue or process in these tissues.[@b44-ebo-9-2013-437]

*In silico* identification and validation of SNP
------------------------------------------------

The availability of huge plant sequence data in the public domain is a rich resource for discovery of high-quality SNPs using bioinformatics pipeline. EST databases are an abundant source of SNP markers.[@b45-ebo-9-2013-437] Recently, Mondini et al[@b46-ebo-9-2013-437] identified SNPs involved in drought and salt stress tolerance in durum wheat. In this study, 1381 aquaporin ESTs were assembled into 32 contigs. Assembly of these cleaned ESTs was done with stringent parameters (match size = 40, sequence length = 100, maximum expected coverage = 40 and match percentage = 95). Only contigs containing a minimum of 4 and maximum of 80 members with 4 different cultivars were analyzed further. Not enough information can be gained from contigs \< 4 EST and it becomes difficult to view and edit contigs \> 80 ESTs. SNPs were declared only when there was no mismatch, no gaps before and after putative SNP site; in addition, the alternative base in the consensus sequence was present at least more than twice in an alignment. 6 SNPs from 9 contigs were taken into consideration with SNP score ≥ 40 and SNPs in the start and end of the alignment were ignored ([Table 4](#t4-ebo-9-2013-437){ref-type="table"}). In current investigation, we found only a few SNPs because less EST sequences were available for the aquaporin gene. The allele-specific marker ([Supplementary Table 2](#s1-ebo-9-2013-437){ref-type="supplementary-material"}) was used for amplification of DNA fragment in 38 wheat genotypes ([Supplementary Table 3](#s1-ebo-9-2013-437){ref-type="supplementary-material"}), representing a core set.[@b47-ebo-9-2013-437] 20 genotypes resulted in the presence of a band of \~170 bp fragments ([Fig. 5](#f5-ebo-9-2013-437){ref-type="fig"}), and the remaining 18 genotypes did not yield any amplification. 1 SNP was found between these genotypes. The primer is amplified mostly in drought tolerant genotypes, based on the reported phenotypic data analysis (BS Tyagi, personal communication). This is the first report of a PIP-derived SNP marker in wheat that can be used by the breeders for improving tolerance to drought and high temperatures in wheat breeding programs.

3D model prediction and validation
----------------------------------

After establishing functional annotation, we have now extended research on the structural background of aquaporin protein in order to gain more insights. Lack of theoretical 3D structure of TaAQP protein motivated us to carry out this study. Therefore, we modeled wheat PIP (Accession no: AAM00368) using BLASTP identified the crystal structure with PDB ID: 2B5F as templates. Query sequences showed 74% sequence identity with the template with an E-value of 1e-59 ([Fig. 6A](#f6-ebo-9-2013-437){ref-type="fig"}). Only 1 (0.4%) out of 292 residues was present in the disallowed region whereas another 4 (1.7%) residues were present in the generously allowed regions of the Ramachandran plot, respectively ([Supplementary Fig. 3](#s1-ebo-9-2013-437){ref-type="supplementary-material"}). Similarly, despite of having \~74% sequence homology amongst template and modeled proteins, the tertiary structure was also found to be comparable as indicated by a root mean square deviation of 0.60 Å. Holistically, the modeled tertiary structure of *T. aestivum* PIP is structurally very similar to the spinach aquaporin transport protein. The modeled TaPIP protein comprised of 6 helices, 18 helix-helix interactions, 14 beta turns, and 3 gamma turns were predicted in the 3D structure of aquaporin, implying a helix-rich structure ([Fig. 6B](#f6-ebo-9-2013-437){ref-type="fig"}). The final protein structure was deposited in PMDB and is available under the Accession ID: PM0078491.

We believe that the present findings illustrate more insights into the structure-function role of AQPs protein in molecular terms. In order to understand the function of individual MIPs in maintaining water homeostasis, it is necessary to carry out knock-out experiments and promoter analyses, as well as substrate specificities under various physiological conditions in relation to water balance and nutrient uptake in wheat and other plant systems. Current biotechnology and bioinformatics tools may identify and characterize genes in their respective subclasses. As a model plant, and having a great synteny with the grass family with respect to gene structure, the information generated about the aquaoporin gene family in wheat will also provide a platform for predicting the function of genes of crops whose genome sequences are in their infancy.

Conclusions
===========

The significance of the multigene family of aquaporin transmembrane proteins is emerging from studies aimed at optimizing water and nutrient use efficiency. Complete set of AQPs have been identified and classified in some agriculturally important crops with well-determined genomic sequences. Though the global importance of wheat as a most important cereal grain in world trade is well established, its AQPs remain less studied to date because of the lack of wheat genomic sequence information. The goal of this study was to identify new members of the aquaporin family in the wheat genome. Our combined approach enabled us to identify a total of 13 new aquaporin genes in wheat that showed significant sequence identity (\>50% identity) with those from rice. A number of motifs and signature patterns were identified that are related to subcellular localization and functional annotation. Characterization of SNPs in candidate genes for drought tolerance such as aquaporin is a promising approach for identifying alleles that are associated with drought phenotypes. The aquaporin gene obtained in this work thus provides further tools for the physical and genetic mapping of these important genes, for identifying their chromosomal locations or genetic linkage to water homeostasis-related traits, respectively.

Supplementary Materials
=======================

###### 

Supplementary Figure 1

Transmemebrane α helix from the amino acid sequence of wheat AQPs based on SMART analysis.

Supplementary Figure 2

Multiple sequence alignment of wheat aquaporins. Deduced amino acid sequences were aligned using the ClustalW program. NPA motifs, black box. The MIP family signature sequences are in red box and the consensus phosphorylation site located in the C-terminal region is shaded yellow.

Supplementary Figure 3

Model structure of PIP evaluated by Ramachandran's plot.

Supplementary Table 1

The ClustalW multiple sequence alignment of identified wheat aquaporin protein sequences.

Supplementary Table 2

A list of the EST based SNP primers designed for aquaporin.

Supplementary Table 3

List of released Indian wheat genotypes used in this study.

This study is part of the requirements for BP's doctoral degree in Biotechnology. The financial support for Agri-Bioinformatics Promotion Program provided by Bioinformatics Initiative Division, Department of Information Technology, Ministry of Communications and Information Technology, Government of India, New Delhi is gratefully acknowledged.

**Author Contributions**

Conceived and designed the experiments: PS, BP. Analyzed the data: BP. Wrote the first draft of the manuscript: BP, PS. Contributed to the writing of the manuscript: BP, PS. Agree with manuscript results and conclusions: BP, DMP, PS, RC, IS. Jointly developed the structure and arguments for the paper: BP, PS. Made critical revisions and approved final version: BP, PS, DMP, RC, IS. All authors reviewed and approved of the final manuscript.

**Competing Interests**

Author(s) disclose no potential conflicts of interest.

**Disclosures and Ethics**

As a requirement of publication the authors have provided signed confirmation of their compliance with ethical and legal obligations including but not limited to compliance with ICMJE authorship and competing interests guidelines, that the article is neither under consideration for publication nor published elsewhere, of their compliance with legal and ethical guidelines concerning human and animal research participants (if applicable), and that permission has been obtained for reproduction of any copyrighted material. This article was subject to blind, independent, expert peer review. The reviewers reported no competing interests.

**Funding**

This work was partially funded by Indian Council of Agricultural Research (ICAR) Grant-in-Aid DWR/RP/2010-5.3 to PS and support from Agri-Bioinformatics Promotion Program of Ministry of Communications and Information Technology, New Delhi to RC.

###### 

Multiple sequence alignment using amino acid sequence of (**A**) TaNIP1-1, TaNIP1-3, TaNIP2-1, TaNIP4-1 (**B**) TaPIP1-1 and TaPIP1-3 (**C**) TaPIP2-2, TaPIP2-3 and TaPIP2-4 (**D**) TaTIP1-4, TaSIP1-1 and TaSIP1-2 from different cereal crops. Conserved amino acids are shown in black boxes. Protein structural features are indicated above the alignment. Alpha helix and beta strands are elements of secondary structure represented as rods and arrows.

![](ebo-9-2013-437f1a)

![](ebo-9-2013-437f1b)

![](ebo-9-2013-437f1c)

![](ebo-9-2013-437f1d)

![Phylogenetic tree was constructed using the neighbor-joining method and diagrams drawn with MEGA4. NJ method used the multiple sequence alignment generated by ClustalW to generate the tree. Bootstrap values are indicated against each branch. Phylogenetic tree showing the 4 clusters of PIPs, TIPs, NIPs and SIPs. The 13 wheat AQPs are compared with all the PIPs as well as the TIPs, NIPs and SIPs from *A. thaliana*, *G. max,* and *B. distachyon, H. vulgare, V. vinifera* and *Z. mays.* For visibility reasons identified wheat aquaporin are indicated by a red triangle. Branch lengths are proportional to evolutionary distance.](ebo-9-2013-437f2){#f2-ebo-9-2013-437}

![Distribution of *T. aestivum* aquaporin family genes in various tissues. The expression profile was determined by analyzing the EST counts of 9 different kinds of tissues.](ebo-9-2013-437f3){#f3-ebo-9-2013-437}

![Tissue specific expression of the aquaporin genes profile was determined by analyzing the EST counts based on UniGene.](ebo-9-2013-437f4){#f4-ebo-9-2013-437}

![PCR analysis using the SNP marker. List of genotype was mentioned in [Supplementary Table 2](#s1-ebo-9-2013-437){ref-type="supplementary-material"}. The size of the amplification product is shown on the left. Absence of a band indicates the specific sequence is absent in wheat genotype. Lane 39 is a negative control, which used H~2~O as the template.](ebo-9-2013-437f5){#f5-ebo-9-2013-437}

![(**A**) Pairwise alignment of aquaporin family gene from *Triticum aestivum* and the template (PDB ID: 1N0J). The dash represents insertion and deletion; conserved residues showed in shaded region, NPA box is labeled and the aquaporin signature is shown in red (**B**) Aquaporin structure model produced using Modeller9.](ebo-9-2013-437f6){#f6-ebo-9-2013-437}

###### 

Blast results of identified aquaporin to classify them accordingly.

  Aquaporin identified            No. of ESTs in TC   TC length   Rice gene   Chromosome no   Cellular localization   Bit score   \% identity   E-value
  ------------------------------- ------------------- ----------- ----------- --------------- ----------------------- ----------- ------------- ----------
  **Blast sequences alignment**                                                                                                                 
  TaNIP1-1                        24                  278         OsNIP1-1    2               PM and M                466         86            4e-165
  TaNIP1-3                        24                  153         OsNIP1-3    --              PM and M                459         76            7e-68
  TaNIP2-1                        41                  295         OsNIP2-1    2               PM and M                510         86            0
  TaNIP4-1                        6                   290         OsNIP4-1    1               PM and M                323         73            8e-109
  TaPIP1-1                        64                  290         OsPIP1-1    3               M                       483         87            2e-171
  TaPIP1-3                        32                  292         OsPIP1-3    --              M                       535         90            0
  TaPIP2-1                        25                  353         OsPIP2-1    7               M                       250         61            3.00E-77
  TaPIP2-2                        39                  290         OsPIP2-2    2               PM and M                489         83            2e-175
  TaPIP2-3                        8                   286         OsPIP2-4    --              PM and M                494         86            3e-175
  TaPIP2-4                        3                   286         OsPIP2-6    4               PM and M                526         90            0
  TaSIP2-1                        38                  244         OsSIP2-1    --              ER                      218         51            2e-70
  TaSIP1-1                        38                  128         OsSIP1-1    --              ER                      207         88            2e-66
  TaTIP4-1                        7                   252         OsTIP4-1    5               T                       370         86            3e-130

**Abbreviations:** TC, tentative consensus sequence; PM, plasma membrane; M, mitochondria; T, tonoplast; ER, endoplasmic reticulum.

###### 

Summary of the aquaporin family in different plant species.

  Plant       *A. thaliana*   *O. sativa*   *Z. mays*   *V. vinifera*   *H. vulgare*   *T. aestivum*
  ----------- --------------- ------------- ----------- --------------- -------------- --------------------
  **Class**                                                                            *Newly identified*
   NIP                                                                                 
    NIP1-1    3               1             --          1               1              1
    NIP1-3    --              1             --          --              --             1
    NIP2-1    3               1             --          1               1              1
    NIP4-1    1               1             --          --              1              1
   PIP                                                                                 
    PIP1-1    5               1             1           1               1              1
    PIP1-3    5               1             4           1               1              1
    PIP2-1    7               1             4           5               1              1
    PIP2-2    3               1             2           1               1              1
    PIP2-3    3               1             1           1               1              1
    PIP2-4    3               1             2           1               1              1
   SIP                                                                                 
    SIP2-1    1               1             1           --              --             1
    SIP1-1    3               1             --          --              1              1
   TIP                                                                                 
    TIP4-1    3               1             --          1               1              1

###### 

Fingerprint analysis result for different aquaporin classes.

  MIP class   MIP subclass   MIP family signature   N-myristoylation site                                                                                               PKC_PHOSPHO_SITE               CK2_PHOSPHO_SITE
  ----------- -------------- ---------------------- ------------------------------------------------------------------------------------------------------------------- ------------------------------ ------------------------------------
  NIP         NIP1-1         HFNPAVTFA              \[G\]\[A/L/T/C/G\]\[V/A/I/T/\]\[L/V/P/A\]\[G/T/A\]\[S/G/A\]                                                         \[T\]\[V/I/N\]\[R/K\]          TGSD
              NIP1-3         HMNPAVTLA              \[G\]\[A/L/T/A/C/G\]\[V/A/I/T/G/D\]\[A/V/P/I/N/D\]\[G/T/A/S/N\]\[A/G/S/I/Q\]                                        \[T\]\[I/N\]\[K/R\]            TGSD
              NIP2-1         HMNPAVTLA              \[G\]\[A/L/T/A/C/G\]\[V/A/I/T/G/D\]\[A/V/P/I/N/D\]\[G/T/A/S/N\]\[A/G/S/I/Q\]                                        \[T/S\]\[E/N/F\]\[R/K\]        \[S\]\[I/A\]\[A/H\]\[D\]
              NIP4-1         HFNPAVTFA              \[G\]\[S/A/L/M/V/G\]\[R/P/A/N/I\]\[G/K/V/I/N/Q\]\[T/A/S/N/G\]\[A/Q/H/T/R/D/S/N\]                                    \[T\]\[V/I/N\]\[R/K\]          TGSD
  PIP         PIP1-1         HINPAVTFG              \[G\]\[G/I/L/A/T\]\[D/A/S/F/M/V/N/H\]\[F/W/G/I/C/K/N/P/Q\]\[A/G/N/T\]\[A/F/H/R/V/T/I/G/L\]                          SAR                            \[T/S\]\[K/A/P\]\[G/T/R/Q\]\[D/E\]
              PIP1-3         HINPAVTFG              \[G\]\[N/S/I/L/A/V/N/G\]\[P/K/A/S/F/I/V/N\]\[P/K/A/S/F/I/V/G/N\]\[S/C/W/G/I/C/K/N/P\]\[G/S/A/N\]\[S/T/F/H/R/A/V\]   \[S\]\[A/E\]\[R\]              \[T/S\]\[K/A\]\[G/T/R\]\[D\]
              PIP2-1         HINPAVTFG              \[G\]\[G/I/L/A/T\]\[D/A/S/F/M/V/N/H\]\[F/W/G/I/C/K/N/P/Q\]\[A/G/N/T\]\[A/F/H/R/V/T/I/G/L\]                          \[S/T\]\[A/M\]\[R\]            \[T/S\]\[A\]\[A/T/R\]\[D\]
              PIP2-2         HINPAVTFG              \[G\]\[G/I/L/A/T\]\[D/A/S/F/M/V/N/H\]\[F/W/G/I/C/K/N/P/Q\]\[A/G/N/T\]\[A/F/H/R/V/T/I/G/L\]                          \[S\]\[A/F\]\[R\]              \[S\]\[A\]\[T/R\]\[D\]
              PIP2-3         HINPAVTFG              \[G\]\[N/S/I/L/A/V/N/G\]\[P/K/A/S/F/I/V/N\]\[P/K/A/S/F/I/V/G/N\]\[S/C/W/G/I/C/K/N/P\]\[G/S/A/N\]\[S/T/F/H/R/A/V\]   \[S\]\[K/A/F\]\[R\]            \[S\]\[G/A\]\[P/T\]\[D\]
              PIP2-4         HINPAVTFG              \[G\]\[A/L/T/C/G\]\[V/A/I/T/G/D\]\[A/V/P/I/N/D\]\[G/T/A/S/N\]\[A/G/S/L/I/Q\]                                        RGRR                           TAAD
  SIP         SIP1-1         HINPAVTFG              \[G\]\[A/I/G\]\[S/H/V/I\]\[T/Y/G/I\]\[A/G/T\]\[A/L/I/G\]                                                            \[S/T\]\[I/F/K\]\[K/R\]        SLHE
              SIP2-1         HINPAVTFG              \[G\]\[A\]\[I/A\]\[Y/L\]\[T/A\]\[G/A\]                                                                              \[S/T\]\[F/L\]\[K/R\]          --
  TIP         TIP4-1         HLNPAVTVA              \[G\]\[V/L/N/G\]\[P/W/F/N/I/A\]\[I/V/F/S/T/A\]\[A/T/G/S/N\]\[Q/G/A/H/G\]                                            \[S/T\]\[S/I\]\[P/I\]\[D/E\]   --

###### 

Details of ESTs and predicted SNPs in aquaporin genes.

  Gene name   No. of EST   Unigene no.   Contigs   Singletons   SNP (\>40%)
  ----------- ------------ ------------- --------- ------------ -------------
  PIP1:2      886          Ta.50448      6         --           4
  PIP1        336          Ta.23833      3         --           2
