Genome-wide identification and characterization of aquaporin gene family in common bean (Phaseolus vulgaris L.)

Identification of VrNIP2-1 aquaporin with novel selective filter regulating the transport of beneficial as well as hazardous metalloids in mungbean (Vigna radiata L.)

Nodulin 26-like intrinsic protein (NIP) subfamily of aquaporins (AQPs) in plants, is known to be involved in the uptake of metalloids including boron, germanium (Ge), arsenic (As), and silicon (Si). In the present study, a thorough evaluation of 55 AQPs found in the mungbean genome, including phylogenetic distribution, sequence homology, expression profiling, and structural characterization, contributed to the identification of VrNIP2-1 as a metalloid transporter. The pore-morphology of VrNIP2-1 was studied using molecular dynamics simulation. Interestingly, VrNIP2-1 was found to harbor an aromatic/arginine (ar/R) selectivity filter formed with ASGR amino acids instead of GSGR systematically reported in metalloid transporters (NIP2s) in higher plants. Evaluation of diverse cultivars showed a high level of Si accumulation in leaves indicating functional Si transport in mungbean. In addition, heterologous expression of VrNIP2-1 in yeast revealed As(III) and GeO2 transport activity. Similarly, VrNIP2-1 expression in Xenopus oocytes confirmed its Si transport ability. The metalloid transport activity with unique structural features will be helpful to better understand the solute specificity of NIP2s in mungbean and related pulses. The information provided here will also serve as a basis to improve Si uptake while restricting hazardous metalloids like As in plants.

Selection signatures and population dynamics of transposable elements in lima bean

The domestication process in lima bean (Phaseolus lunatus L.) involves two independent events, within the Mesoamerican and Andean gene pools. This makes lima bean an excellent model to understand convergent evolution. The mechanisms of adaptation followed by Mesoamerican and Andean landraces are largely unknown. Genes related to these adaptations can be selected by identification of selective sweeps within gene pools. Previous genetic analyses in lima bean have relied on Single Nucleotide Polymorphism (SNP) loci, and have ignored transposable elements (TEs). Here we show the analysis of whole-genome sequencing data from 61 lima bean accessions to characterize a genomic variation database including TEs and SNPs, to associate selective sweeps with variable TEs and to predict candidate domestication genes. A small percentage of genes under selection are shared among gene pools, suggesting that domestication followed different genetic avenues in both gene pools. About 75% of TEs are located close to genes, which shows their potential to affect gene functions. The genetic structure inferred from variable TEs is consistent with that obtained from SNP markers, suggesting that TE dynamics can be related to the demographic history of wild and domesticated lima bean and its adaptive processes, in particular selection processes during domestication.

Genome-Wide Identification and Gene Expression Analysis of Sweet Cherry Aquaporins (Prunus avium L.) under Abiotic Stresses

Aquaporins (AQPs) are integral transmembrane proteins well known as channels involved in the mobilization of water, small uncharged molecules and gases. In this work, the main objective was to carry out a comprehensive study of AQP encoding genes in Prunus avium (cv. Mazzard F12/1) on a genome-wide scale and describe their transcriptional behaviors in organs and in response to different abiotic stresses. A total of 28 non-redundant AQP genes were identified in Prunus spp. Genomes, which were phylogenetically grouped into five subfamilies (seven PIPs, eight NIPs, eight TIPs, three SIPs and two XIPs). Bioinformatic analyses revealed a high synteny and remarkable conservation of structural features among orthologs of different Prunus genomes. Several cis-acting regulatory elements (CREs) related to stress regulation were detected (ARE, WRE3, WUN, STRE, LTR, MBS, DRE, AT-rich and TC-rich). The above could be accounting for the expression variations associated with plant organs and, especially, each abiotic stress analyzed. Gene expressions of different PruavAQPs were shown to be preferentially associated with different stresses. PruavXIP2;1 and PruavXIP1;1 were up-regulated in roots at 6 h and 72 h of hypoxia, and in PruavXIP2;1 a slight induction of expression was also detected in leaves. Drought treatment strongly down-regulated PruavTIP4;1 but only in roots. Salt stress exhibited little or no variation in roots, except for PruavNIP4;1 and PruavNIP7;1, which showed remarkable gene repression and induction, respectively. Interestingly, PruavNIP4;1, the AQP most expressed in cherry roots subjected to cold temperatures, also showed this pattern in roots under high salinity. Similarly, PruavNIP4;2 consistently was up-regulated at 72 h of heat and drought treatments. From our evidence is possible to propose candidate genes for the development of molecular markers for selection processes in breeding programs for rootstocks and/or varieties of cherry.

A transcriptional complex of FtMYB102 and FtbHLH4 coordinately regulates the accumulation of rutin in Fagopyrum tataricum

Tartary buckwheat is rich in flavonoids, which not only play an important role in the plant-environment interaction, but are also beneficial to human health. Rutin is a therapeutic flavonol which is massively accumulated in Tartary buckwheat. It has been demonstrated that transcription factors control rutin biosynthesis. However, the transcriptional regulatory network of rutin is not fully clear. In this study, through transcriptome and target metabolomics, we validated the role of FtMYB102 and FtbHLH4 TFs at the different developmental stages of Tartary buckwheat. The elevated accumulation of rutin in the sprout appears to be closely associated with the expression of FtMYB102 and FtbHLH4. Yeast two-hybrid, transient luciferase activity and co-immunoprecipitation demonstrated that FtMYB102 and FtbHLH4 can interact and form a transcriptional complex. Moreover, yeast one-hybrid showed that both FtMYB102 and FtbHLH4 directly bind to the promoter of chalcone isomerase (CHI), and they can coordinately induce CHI expression as shown by transient luciferase activity assay. Finally, we transferred FtMYB102 and FtbHLH4 into the hairy roots of Tartary buckwheat and found that they both can promote the accumulation of rutin. Our results indicate that FtMYB102 and FtbHLH4 can form a transcriptional complex by inducing CHI expression to coordinately promote the accumulation of rutin.

Genome-Wide Identification of Aquaporin Genes in Adzuki Bean (Vigna angularis) and Expression Analysis under Drought Stress

The adzuki bean Vigna angularis (Wild.) is an important leguminous crop cultivated mainly for food purposes in Asian countries; it represents a source of carbohydrates, digestible proteins, minerals, and vitamins. Aquaporins (AQPs) are crucial membrane proteins involved in the transmembrane diffusion of water and small solutes in all living organisms, including plants. In this study, we used the whole genome sequence of the adzuki bean for in silico analysis to comprehensively identify 40 Vigna angularis aquaporin (VaAQP) genes and reveal how these plants react to drought stress. VaAQPs were compared with AQPs from other closely-related leguminous plants, and the results showed that mustard (Brassica rapa) (59), barrel medic (Medicago truncatula) (46), soybean (Glycine max) (66), and common bean (Phaseolus vulgaris L.) (41) had more AQP genes. Phylogenetic analysis revealed that forty VaAQPs belong to five subfamilies, with the VaPIPs (fifteen) subfamily the largest, followed by the VaNIPs (ten), VaTIPs (ten), VaSIPs (three), and VaXIPs (two) subfamilies. Furthermore, all AQP subcellular locations were found at the plasma membrane, and intron-exon analysis revealed a relationship between the intron number and gene expression, duplication, evolution, and diversity. Among the six motifs identified, motifs one, two, five, and six were prevalent in VaTIP, VaNIP, VaPIP, and VaXIP, while motifs one, three, and four were not observed in VaPIP1-3 and VaPIP1-4. Under drought stress, two of the VaAQPs (VaPIP2-1 and VaPIP2-5) showed significantly higher expression in the root tissue while the other two genes (VaPIP1-1 and VaPIP1-7) displayed variable expression in leaf tissue. This finding revealed that the selected VaAQPs might have unique molecular functions linked with the uptake of water under drought stress or in the exertion of osmoregulation to transport particular substrates rather than water to protect plants from drought. This study presents the first thorough investigation of VaAQPs in adzuki beans, and it reveals the transport mechanisms and related physiological processes that may be utilized for the development of drought-tolerant adzuki bean cultivars.

Genome-Wide Identification of Common Bean PvLTP Family Genes and Expression Profiling Analysis in Response to Drought Stress

Common bean is one of the most important legume crops for human consumption. Its yield is adversely affected by environmental stress. Plant non-specific lipid transfer proteins (nsLTPs) are essential for plant growth, development, and resistance to abiotic stress, such as salt, drought, and alkali. However, changes in nsLTP family genes responding to drought stress are less known. The PvLTP gene family in the common bean was identified by a comprehensive genome-wide analysis. Molecular weights, theoretical isoelectric points, phylogenetic tree, conserved motifs, gene structures, gene duplications, chromosome localization, and expression profiles were analyzed by SignalP 5.0, ExPASy, ClustalX 2.1, MEGA 7.0, NCBI-CDD, MEME, Weblogo, and TBtools 1.09876, respectively. Heatmap and qRT-PCR analyses were performed to validate the expression profiles of PvLTP genes in different organs. In addition, the expression patterns of nine PvLTP genes in common beans treated with drought stress were investigated by qRT-PCR. We obtained 58 putative PvLTP genes in the common bean genome via genome-wide analyses. Based on the diversity of the eight-cysteine motif (ECM), these genes were categorized into five types (I, II, IV, V, and VIII). The signal peptides of the PvLTP precursors were predicted to be from 16 to 42 amino acid residues. PvLTPs had a predicated theoretical isoelectric point of 3.94-10.34 and a molecular weight of 7.15-12.17 kDa. The phylogenetic analysis showed that PvLTPs were closer to AtLTPs than OsLTPs. Conserved motif and gene structure analyses indicated that PvLTPs were randomly distributed on all chromosomes except chromosome 9. In addition, 23 tandem duplicates of PvLTP genes were arranged in 10 gene clusters on chromosomes 1 and 2. The heatmap and qRT-PCR showed that PvLTP expression significantly varied in different tissues. Moreover, 9 PvLTP genes were up-regulated under drought treatment. Our results reveal that PvLTPs play potentially vital roles in plants and provide a comprehensive reference for studies on PvLTP genes and a theoretical basis for further analysis of regulatory mechanisms influencing drought tolerance in the common bean.

Identification and Characterization of Salt- and Drought-Responsive AQP Family Genes in Medicagosativa L

Aquaporins (AQP) are distributed ubiquitously in plants, and they play important roles in multiple aspects of plant growth and development, as well as in plant resistance to various environmental stresses. In this study, 43 MsAQP genes were identified in the forage crop Medicago sativa. All the MsAQP proteins were clustered into four subfamilies based on sequence similarity and phylogenetic relationship, including 17 TIPs, 14 NIPs, 9 PIPs and 3 SIPs. Analyses of gene structure and conserved domains indicated that the majority of the deduced MsAQP proteins contained the signature transmembrane domains and the NPA motifs. Analyses on cis-acting elements in the promoter region of MsAQP genes revealed the presence of multiple and diverse stress-responsive and hormone-responsive cis-acting elements. In addition, by analyzing the available and comprehensive gene expression data of M. truncatula, we screened ten representative MtAQP genes that were responsive to NaCl or drought stress. By analyzing the sequence similarity and phylogenetic relationship, we finally identified the corresponding ten salt- or drought-responsive AQP genes in M. sativa, including three MsTIPs, three MsPIPs and four MsNIPs. The qPCRs showed that the relative expression levels of these ten selected MsAQP genes responded differently to NaCl or drought treatment in M. sativa. Gene expression patterns showed that most MsAQP genes were preferentially expressed in roots or in leaves, which may reflect their tissue-specific functions associated with development. Our results lay an important foundation for the future characterization of the functions of MsAQP genes, and provide candidate genes for stress resistance improvement through genetic breeding in M. sativa.

The aquaporin gene PvXIP1;2 conferring drought resistance identified by GWAS at seedling stage in common bean

A whole-genome resequencing-derived SNP dataset used for genome-wide association analysis revealed 12 loci significantly associated with drought stress based on survival rate after drought stress at seedling stage. We further confirmed the drought-related function of an aquaporin gene (PvXIP1;2) located at Locus_10. A variety of adverse conditions, including drought stress, severely affect common bean production. Molecular breeding for drought resistance has been proposed as an effective and practical way to improve the drought resistance of common bean. A genome-wide association analysis was conducted to identify drought-related loci based on survival rates at the seedling stage using a natural population consisting of 400 common bean accessions and 3,832,340 SNPs. The coefficient of variation ranged from 40.90 to 56.22% for survival rates in three independent experiments. A total of 12 associated loci containing 89 significant SNPs were identified for survival rates at the seedling stage. Four loci overlapped in the region of the QTLs reported to be associated with drought resistance. According to the expression profiles, gene annotations and references of the functions of homologous genes in Arabidopsis, 39 genes were considered potential candidate genes selected from 199 genes annotated within all associated loci. A stable locus (Locus_10) was identified on chromosome 11, which contained LEA, aquaporin, and proline-rich protein genes. We further confirmed the drought-related function of an aquaporin (PvXIP1;2) located at Locus_10 by expression pattern analysis, phenotypic analysis of PvXIP1;2-overexpressing Arabidopsis and Agrobacterium rhizogenes-mediated hairy root transformation systems, indicating that the association results can facilitate the efficient identification of genes related to drought resistance. These loci and their candidate genes provide a foundation for crop improvement via breeding for drought resistance in common bean.

Genome-Wide Identification and Function of Aquaporin Genes During Dormancy and Sprouting Periods of Kernel-Using Apricot (Prunus armeniaca L.)

Aquaporins (AQPs) are essential channel proteins that play a major role in plant growth and development, regulate plant water homeostasis, and transport uncharged solutes across biological membranes. In this study, 33 AQP genes were systematically identified from the kernel-using apricot (Prunus armeniaca L.) genome and divided into five subfamilies based on phylogenetic analyses. A total of 14 collinear blocks containing AQP genes between P. armeniaca and Arabidopsis thaliana were identified by synteny analysis, and 30 collinear blocks were identified between P. armeniaca and P. persica. Gene structure and conserved functional motif analyses indicated that the PaAQPs exhibit a conserved exon-intron pattern and that conserved motifs are present within members of each subfamily. Physiological mechanism prediction based on the aromatic/arginine selectivity filter, Froger's positions, and three-dimensional (3D) protein model construction revealed marked differences in substrate specificity between the members of the five subfamilies of PaAQPs. Promoter analysis of the PaAQP genes for conserved regulatory elements suggested a greater abundance of cis-elements involved in light, hormone, and stress responses, which may reflect the differences in expression patterns of PaAQPs and their various functions associated with plant development and abiotic stress responses. Gene expression patterns of PaAQPs showed that PaPIP1-3, PaPIP2-1, and PaTIP1-1 were highly expressed in flower buds during the dormancy and sprouting stages of P. armeniaca. A PaAQP coexpression network showed that PaAQPs were coexpressed with 14 cold resistance genes and with 16 cold stress-associated genes. The expression pattern of 70% of the PaAQPs coexpressed with cold stress resistance genes was consistent with the four periods [Physiological dormancy (PD), ecological dormancy (ED), sprouting period (SP), and germination stage (GS)] of flower buds of P. armeniaca. Detection of the transient expression of GFP-tagged PaPIP1-1, PaPIP2-3, PaSIP1-3, PaXIP1-2, PaNIP6-1, and PaTIP1-1 revealed that the fusion proteins localized to the plasma membrane. Predictions of an A. thaliana ortholog-based protein-protein interaction network indicated that PaAQP proteins had complex relationships with the cold tolerance pathway, PaNIP6-1 could interact with WRKY6, PaTIP1-1 could interact with TSPO, and PaPIP2-1 could interact with ATHATPLC1G. Interestingly, overexpression of PaPIP1-3 and PaTIP1-1 increased the cold tolerance of and protein accumulation in yeast. Compared with wild-type plants, PaPIP1-3- and PaTIP1-1-overexpressing (OE) Arabidopsis plants exhibited greater tolerance to cold stress, as evidenced by better growth and greater antioxidative enzyme activities. Overall, our study provides insights into the interaction networks, expression patterns, and functional analysis of PaAQP genes in P. armeniaca L. and contributes to the further functional characterization of PaAQPs.

Plant aquaporins: A frontward to make crop plants drought resistant

Drought stress alters gene expression and causes cellular damage in crop plants. Drought inhibits photosynthesis by reducing the content and the activity of the photosynthetic carbon reduction cycle, ultimately decreasing the crop yield. The role of aquaporins (AQP) in improving the growth and adaptation of crop plants under drought stress is of importance. AQP form channels and control water transport in and out of the cells and are associated with drought tolerance mechanisms. The current review addresses: (1) the evolution of AQPs in plants, (2) the classification of plant AQPs, (3) the role of AQPs in drought alleviation in crop plants, and (4) the phytohormone crosstalk with AQPs in crops exposed to drought stress.

Plant aquaporins alleviate drought tolerance in plants by modulating cellular biochemistry, root-architecture, and photosynthesis

Water is a vital resource for plants to grow, thrive, and complete their life cycle. In recent years, drastic changes in the climate, especially drought frequency and severity, have increased, which reduces agricultural productivity worldwide. Aquaporins are membrane channels belonging to the major intrinsic protein superfamily, which play an essential role in cellular water and osmotic homeostasis of plants under both control and water deficit conditions. A genome-wide search reveals the vast availability of aquaporin isoforms, phylogenetic relationships, different families, conserved residues, chromosomal locations, and gene structure of aquaporins. Furthermore, aquaporins gating and subcellular trafficking are commonly controlled by phosphorylation, cytosolic pH, divalent cations, reactive oxygen species, and stoichiometry. Researchers have identified their involvement in regulating hydraulic conductance, root system architecture, modulation of abiotic stress-related genes, seed viability and germination, phloem loading, xylem water exit, photosynthetic parameters, and post-drought recovery. Remarkable effects following the change in aquaporin activity and/or gene expression have been observed on root water transport properties, nutrient acquisition, physiology, transpiration, stomatal aperture, gas exchange, and water use efficiency. The present review highlights the role of different aquaporin homologs under water-deficit stress condition in model and crop plants. Moreover, the opportunity and challenges encountered to explore aquaporins for engineering drought-tolerant crop plants are also discussed here.

Root-omics for drought tolerance in cool-season grain legumes

Root traits can be exploited to increase the physiological efficiency of crop water use under drought. Root length, root hairs, root branching, root diameter, and root proliferation rate are genetically defined traits that can help to improve the water productivity potential of crops. Recently, high-throughput phenotyping techniques/platforms have been used to screen the germplasm of major cool-season grain legumes for root traits and their impact on different physiological processes, including nutrient uptake and yield potential. Advances in omics approaches have led to the dissection of genomic, proteomic, and metabolomic structures of these traits. This knowledge facilitates breeders to improve the water productivity and nutrient uptake of cultivars under limited soil moisture conditions in major cool-season grain legumes that usually face terminal drought. This review discusses the advances in root traits and their potential for developing drought-tolerant cultivars in cool-season grain legumes.

Understanding aquaporin transport system, silicon and other metalloids uptake and deposition in bottle gourd (Lagenaria siceraria)

Aquaporins (AQPs) facilitates the transport of small solutes like water, urea, carbon dioxide, boron, and silicon (Si) and plays a critical role in important physiological processes. In this study, genome-wide characterization of AQPs was performed in bottle gourd. A total of 36 AQPs were identified in the bottle gourd, which were subsequently analyzed to understand the pore-morphology, exon-intron structure, subcellular-localization. In addition, available transcriptome data was used to study the tissue-specific expression. Several AQPs showed tissue-specific expression, more notably the LsiTIP3-1 having a high level of expression in flowers and fruits. Based on the in-silico prediction of solute specificity, LsiNIP2-1 was predicted to be a Si transporter. Silicon was quantified in different tissues, including root, young leaves, mature leaves, tendrils, and fruits of bottle gourd plants. More than 1.3% Si (d.w.) was observed in bottle gourd leaves, testified the in-silico predictions. Silicon deposition evaluated with an energy-dispersive X-ray coupled with a scanning electron microscope showed a high Si accumulation in the shaft of leaf trichomes. Similarly, co-localization of Si with arsenic and antimony was observed. Expression profiling performed with real-time quantitative PCR showed differential expression of AQPs in response to Si supplementation. The information provided in the present study will be helpful to better understand the AQP transport mechanism, particularly Si and other metalloids transport and localization in plants.

Targeting aquaporins to alleviate hazardous metal(loid)s imposed stress in plants

Uptake of hazardous metal(loid)s adversely affects plants and imposes a threat to the entire food chain. Here, the role of aquaporins (AQPs) providing tolerance against hazardous metal(loid)s in plants is discussed to provide a perspective on the present understanding, knowledge gaps, and opportunities. Plants adopt complex molecular and physiological mechanisms for better tolerance, adaptability, and survival under metal(loid)s stress. Water conservation in plants is one such primary strategies regulated by AQPs, a family of channel-forming proteins facilitating the transport of water and many other solutes. The strategy is more evident with reports suggesting differential expression of AQPs adopted by plants to cope with the heavy metal stress. In this regard, numerous studies showing enhanced tolerance against hazardous elements in plants due to AQPs activity are discussed. Consequently, present understanding of various aspects of AQPs, such as tertiary-structure, transport activity, solute-specificity, differential expression, gating mechanism, and subcellular localization, are reviewed. Similarly, various tools and techniques are discussed in detail aiming at efficient utilization of resources and knowledge to combat metal(loid)s stress. The scope of AQP transgenesis focusing on heavy metal stresses is also highlighted. The information provided here will be helpful to design efficient strategies for the development of metal(loid)s stress-tolerant crops.

Transcriptome sequencing revealed the mechanism of promoting floret opening by exogenous methyl jasmonate in sorghum

Flowering time is a critical trait reflecting the adaptation of plants to their environments. Our initial research has shown that exogenous methyl jasmonate (MeJA) significantly promoted the floret opening of sorghum. To better understand the mechanism of this phenomenon in sorghum, the comparative transcriptome analysis was performed. Transcriptomic analysis showed that the most number of differentially expressed genes was presented between control plants and plants treated with 2.0 mM exogenous MeJA in 2.5 h. A large number of differentially expressed genes were assigned to the subcategory of carbohydrate metabolism and lipid metabolism. The transcriptomic analysis of differentially expressed genes involved in glycolysis/gluconeogenesis and tricarboxylic acid cycle indicated a close relationship between carbohydrates metabolism and flowering. In addition, potassium uptake proteins and aquaporins also played important role in response to the exogenous MeJA in the flowering process. These results provide insights into the effect of MeJA on flowering time and explore the possible molecular mechanism of advancing the flowering period by spraying MeJA.

SUPPLEMENTARY INFORMATION

The online version contains supplementary material available at 10.1007/s13205-021-02743-6.

Drought stress-induced physiological mechanisms, signaling pathways and molecular response of chloroplasts in common vegetable crops

Drought stress is one of the most adverse abiotic stresses that hinder plants' growth and productivity, threatening sustainable crop production. It impairs normal growth, disturbs water relations and reduces water-use efficiency in plants. However, plants have evolved many physiological and biochemical responses at the cellular and organism levels, in order to cope with drought stress. Photosynthesis, which is considered one of the most crucial biological processes for survival of plants, is greatly affected by drought stress. A gradual decrease in CO₂ assimilation rates, reduced leaf size, stem extension and root proliferation under drought stress, disturbs plant water relations, reducing water-use efficiency, disrupts photosynthetic pigments and reduces the gas exchange affecting the plants adversely. In such conditions, the chloroplast, organelle responsible for photosynthesis, is found to counteract the ill effects of drought stress by its critical involvement as a sensor of changes occurring in the environment, as the first process that drought stress affects is photosynthesis. Beside photosynthesis, chloroplasts carry out primary metabolic functions such as the biosynthesis of starch, amino acids, lipids, and tetrapyroles, and play a central role in the assimilation of nitrogen and sulfur. Because the chloroplasts are central organelles where the photosynthetic reactions take place, modifications in their physiology and protein pools are expected in response to the drought stress-induced variations in leaf gas exchanges and the accumulation of ROS. Higher expression levels of various transcription factors and other proteins including heat shock-related protein, LEA proteins seem to be regulating the heat tolerance mechanisms. However, several aspects of plastid alterations, following a water deficit environment are still poorly characterized. Since plants adapt to various stress tolerance mechanisms to respond to drought stress, understanding mechanisms of drought stress tolerance in plants will lead toward the development of drought tolerance in crop plants. This review throws light on major droughts stress-induced molecular/physiological mechanisms in response to severe and prolonged drought stress and addresses the molecular response of chloroplasts in common vegetable crops. It further highlights research gaps, identifying unexplored domains and suggesting recommendations for future investigations.

Two Aquaporin Genes, GhPIP2;7 and GhTIP2;1, Positively Regulate the Tolerance of Upland Cotton to Salt and Osmotic Stresses

Aquaporins (AQPs) facilitate the transport of water and small molecules across intrinsic membranes and play a critical role in abiotic stresses. In this study, 111, 54, and 56 candidate AQP genes were identified in Gossypium hirsutum (AD₁), Gossypium arboreum (A₂), and Gossypium raimondii (D₅), respectively, and were further classified into five subfamilies, namely, plasma intrinsic protein (PIP), tonoplast intrinsic protein (TIP), nodulin 26-like intrinsic protein (NIP), small basic intrinsic protein (SIP), and uncategorized X intrinsic protein (XIP). Transcriptome analysis and quantitative real-time PCR (qRT-PCR) revealed some high-expression GhPIPs and GhTIPs (PIP and TIP genes in G. hirsutum, respectively) in drought and salt stresses. GhPIP2;7-silenced plants decreased in the chlorophyll content, superoxide dismutase (SOD) activity, and peroxidase (POD) activity comparing the mock control (empty-vector) under 400 mM NaCl treatment, which indicated a positive regulatory role of GhPIP2;7 in salt tolerance of cotton. The GhTIP2;1-silenced cotton plants were more sensitive to osmotic stress. GhTIP2;1-overexpressed plants exhibited less accumulation of H₂O₂ and malondialdehyde but higher proline content under osmotic stress. In summary, our study elucidates the positive regulatory roles of two GhAQPs (GhPIP2;7 and GhTIP2;1) in salt and osmotic stress responses, respectively, and provides a new gene resource for future research.

Comparative transcriptome and co-expression analysis reveal key genes involved in leaf margin serration in Perilla frutescens

Objective

In this study, we aimed to identify the genes involved in leaf margin serration in Perilla frutescens. P. frutescens (Family: Lamiaceae) is widely grown in Asian countries. Perilla leaf is the medicinal part stipulated in the Chinese Pharmacopoeia. There are mainly two types of perilla leaves: one with serrated leaf margin which is the phenotype described in the pharmacopoeia and the other with smooth leaf margin.

Methods

Transcriptome sequencing, co-expression analysis, and qRT-PCR analysis of six perilla tissues sampled from two different phenotypes (serrated and smooth leaves) were performed.

Results

Forty-three differentially expressed genes (DEGs), which may potentially regulate leaf shape, were identified through de novo transcriptome sequencing between the two groups. Genes involved in leaf shape regulation were identified. Simultaneously, we validated five DEGs by qRT-PCR, and the results were consistent with the transcriptome data. In addition, 1186 transcription factors (TFs) belonging to 45 TF families were identified. Moreover, the co-expression network of DEGs was constructed.

Conclusion

The study identified the key genes that control leaf shape by comparing the transcriptomes. Our findings also provide basic data for further exploring P. frutescens, which can help study the mechanism of leaf shape development and molecular breeding.

Exploration of the Yield Potential of Mesoamerican Wild Common Beans From Contrasting Eco-Geographic Regions by Nested Recombinant Inbred Populations

Genetic analyses and utilization of wild genetic variation for crop improvement in common bean (Phaseolus vulgaris L.) have been hampered by yield evaluation difficulties, identification of advantageous variation, and linkage drag. The lack of adaptation to cultivation conditions and the existence of highly structured populations make association mapping of diversity panels not optimal. Joint linkage mapping of nested populations avoids the later constraint, while populations crossed with a common domesticated parent allow the evaluation of wild variation within a more adapted background. Three domesticated by wild backcrossed-inbred-line populations (BC₁S₄) were developed using three wild accessions representing the full range of rainfall of the Mesoamerican wild bean distribution crossed to the elite drought tolerant domesticated parent SEA 5. These populations were evaluated under field conditions in three environments, two fully irrigated trials in two seasons and a simulated terminal drought in the second season. The goal was to test if these populations responded differently to drought stress and contained progenies with higher yield than SEA 5, not only under drought but also under water-watered conditions. Results revealed that the two populations derived from wild parents of the lower rainfall regions produced lines with higher yield compared to the domesticated parent in the three environments, i.e., both in the drought-stressed environment and in the well-watered treatments. Several progeny lines produced yields, which on average over the three environments were 20% higher than the SEA 5 yield. Twenty QTLs for yield were identified in 13 unique regions on eight of the 11 chromosomes of common bean. Five of these regions showed at least one wild allele that increased yield over the domesticated parent. The variation explained by these QTLs ranged from 0.6 to 5.4% of the total variation and the additive effects ranged from -164 to 277 kg ha^-1, with evidence suggesting allelic series for some QTLs. Our results underscore the potential of wild variation, especially from drought-stressed regions, for bean crop improvement as well the identification of regions for efficient marker-assisted introgression.

Phosphorylation of TIP3 Aquaporins during Phaseolus vulgaris Embryo Development

The membrane phosphoproteome in plant seed changes dynamically during embryo development. We examined the patterns of Phaseolus vulgaris (common bean) seed membrane protein phosphorylation from the mid-maturation stage until two days after germination. Serine and threonine phosphorylation declined during seed maturation while tyrosine phosphorylation remained relatively constant. We discovered that the aquaporin PvTIP3;1 is the primary seed membrane phosphoprotein, and PvTIP3;2 shows a very low level of expression. The level of phosphorylated Ser7 in PvTIP3;1 increased four-fold after seed maturation. Since phosphorylation increases water channel activity, we infer that water transport by PvTIP3;1 is highest in dry and germinating seeds, which would be optimal for seed imbibition. By the use of isoform-specific, polyclonal peptide antibodies, we found that PvTIP3;2 is expressed in a developmental pattern similar to PvTIP3;1. Unexpectedly, PvTIP3;2 is tyrosine phosphorylated following seed maturation, which may suggest a mechanism for the regulation of PvTIP3;2 following seed germination. Analysis of protein secondary structure by circular dichroism spectroscopy indicated that the amino-terminal domain of PvTIP3;1 is generally unstructured, and phosphorylation increases polyproline II (PPII) helical structure. The carboxy-terminal domain also gains PPII character, but in a pH-dependent manner. These structural changes are a first step to understand TIP3 aquaporin regulation.

Plant Aquaporins: Diversity, Evolution and Biotechnological Applications

The plasma membrane forms a permeable barrier that separates the cytoplasm from the external environment, defining the physical and chemical limits in each cell in all organisms. The movement of molecules and ions into and out of cells is controlled by the plasma membrane as a critical process for cell stability and survival, maintaining essential differences between the composition of the extracellular fluid and the cytosol. In this process aquaporins (AQPs) figure as important actors, comprising highly conserved membrane proteins that carry water, glycerol and other hydrophilic molecules through biomembranes, including the cell wall and membranes of cytoplasmic organelles. While mammals have 15 types of AQPs described so far (displaying 18 paralogs), a single plant species can present more than 120 isoforms, providing transport of different types of solutes. Such aquaporins may be present in the whole plant or can be associated with different tissues or situations, including biotic and especially abiotic stresses, such as drought, salinity or tolerance to soils rich in heavy metals, for instance. The present review addresses several aspects of plant aquaporins, from their structure, classification, and function, to in silico methodologies for their analysis and identification in transcriptomes and genomes. Aspects of evolution and diversification of AQPs (with a focus on plants) are approached for the first time with the aid of the LCA (Last Common Ancestor) analysis. Finally, the main practical applications involving the use of AQPs are discussed, including patents and future perspectives involving this important protein family.

Genome-wide characterization and expression analysis of aquaporins in salt cress (Eutrema salsugineum)

Aquaporins (AQPs) serve as water channel proteins and belong to major intrinsic proteins (MIPs) family, functioning in rapidly and selectively transporting water and other small solutes across biological membranes. Importantly, AQPs have been shown to play a critical role in abiotic stress response pathways of plants. As a species closely related to Arabidopsis thaliana, Eutrema salsugineum has been proposed as a model for studying salt resistance in plants. Here we surveyed 35 full-length AQP genes in E. salsugineum, which could be grouped into four subfamilies including 12 plasma membrane intrinsic proteins (PIPs), 11 tonoplast intrinsic proteins (TIPs), nine NOD-like intrinsic proteins (NIPs), and three small basic intrinsic proteins (SIPs) by phylogenetic analysis. EsAQPs were comprised of 237-323 amino acids, with a theoretical molecular weight (MW) of 24.31-31.80 kDa and an isoelectric point (pI) value of 4.73-10.49. Functional prediction based on the NPA motif, aromatic/arginine (ar/R) selectivity filter, Froger's position and specificity-determining position suggested quite differences in substrate specificities of EsAQPs. EsAQPs exhibited global expressions in all organs as shown by gene expression profiles and should be play important roles in response to salt, cold and drought stresses. This study provides comprehensive bioinformation on AQPs in E. salsugineum, which would be helpful for gene function analysis for further studies.

Genome-wide identification and expression analysis of aquaporin gene family related to abiotic stress in watermelon

The plant aquaporins (AQPs) are highly conserved integral membrane proteins that participate in multiple developmental processes and responses to various stresses. In this study, a total of 35 AQP genes were identified in the watermelon genome. The phylogenetic analysis showed that these AQPs can be divided into five types, including 16 plasma membrane intrinsic proteins (PIPs), eight tonoplast intrinsic proteins (TIPs), eight nodulin 26-like intrinsic proteins (NIPs), two small basic intrinsic proteins (SIPs), and one uncategorized X intrinsic protein (XIP). A number of cis-elements related to plant responses to hormones and stresses were detected in the promoter sequences of ClAQP genes. Chromosome distribution analysis revealed that the genes are unevenly distributed on eight chromosomes, with chromosomes 1 and 4 possessing the most genes. Expression analysis at different developmental stages in flesh and rind indicated that most of ClAQPs have tissue-specific expression. Meanwhile, some other AQP genes showed differential expression in response to cold, salt, and ABA treatments, which is consistent with the organization of the stress-responsive cis-elements detected in the promoter regions. Our results lay a foundation for understanding the specific functions of ClAQP genes to help the genetic improvement of watermelon.

Genome-wide identification, structure characterization, and expression pattern profiling of aquaporin gene family in cucumber

BACKGROUND

Aquaporin (AQP) proteins comprise a group of membrane intrinsic proteins (MIPs) that are responsible for transporting water and other small molecules, which is crucial for plant survival under stress conditions including salt stress. Despite the vital role of AQPs, little is known about them in cucumber (Cucumis sativus L.).

RESULTS

In this study, we identified 39 aquaporin-encoding genes in cucumber that were separated by phylogenetic analysis into five sub-families (PIP, TIP, NIP, SIP, and XIP). Their substrate specificity was then assessed based on key amino acid residues such as the aromatic/Arginine (ar/R) selectivity filter, Froger's positions, and specificity-determining positions. The putative cis-regulatory motifs available in the promoter region of each AQP gene were analyzed and results revealed that their promoter regions contain many abiotic related cis-regulatory elements. Furthermore, analysis of previously released RNA-seq data revealed tissue- and treatment-specific expression patterns of cucumber AQP genes (CsAQPs). Three aquaporins (CsTIP1;1, CsPIP2;4, and CsPIP1;2) were the most transcript abundance genes, with CsTIP1;1 showing the highest expression levels among all aquaporins. Subcellular localization analysis in Nicotiana benthamiana epidermal cells revealed the diverse and broad array of sub-cellular localizations of CsAQPs. We then performed RNA-seq to identify the expression pattern of CsAQPs under salt stress and found a general decreased expression level of root CsAQPs. Moreover, qRT-PCR revealed rapid changes in the expression levels of CsAQPs in response to diverse abiotic stresses including salt, polyethylene glycol (PEG)-6000, heat, and chilling stresses. Additionally, transient expression of AQPs in N. benthamiana increased leaf water loss rate, suggesting their potential roles in the regulation of plant water status under stress conditions.

CONCLUSIONS

Our results indicated that CsAQPs play important roles in response to salt stress. The genome-wide identification and primary function characterization of cucumber aquaporins provides insight to elucidate the complexity of the AQP gene family and their biological functions in cucumber.

Combining genome-wide and transcriptome-wide analyses reveal the evolutionary conservation and functional diversity of aquaporins in cotton

BACKGROUND

Aquaporins (AQPs) are integral membrane proteins from a larger family of major intrinsic proteins (MIPs) and function in a huge variety of processes such as water transport, plant growth and stress response. The availability of the whole-genome data of different cotton species allows us to study systematic evolution and function of cotton AQPs on a genome-wide level.

RESULTS

Here, a total of 53, 58, 113 and 111 AQP genes were identified in G. arboreum, G. raimondii, G. hirsutum and G. barbadense, respectively. A comprehensive analysis of cotton AQPs, involved in exon/intron structure, functional domains, phylogenetic relationships and gene duplications, divided these AQPs into five subfamilies (PIP, NIP, SIP, TIP and XIP). Comparative genome analysis among 30 species from algae to angiosperm as well as common tandem duplication events in 24 well-studied plants further revealed the evolutionary conservation of AQP family in the organism kingdom. Combining transcriptome analysis and Quantitative Real-time PCR (qRT-PCR) verification, most AQPs exhibited tissue-specific expression patterns both in G. raimondii and G. hirsutum. Meanwhile, a bias of time to peak expression of several AQPs was also detected after treating G. davidsonii and G. hirsutum with 200 mM NaCl. It is interesting that both PIP1;4 h/i/j and PIP2;2a/e showed the highly conserved tandem structure, but differentially contributed to tissue development and stress response in different cotton species.

CONCLUSIONS

These results demonstrated that cotton AQPs were structural conservation while experienced the functional differentiation during the process of evolution and domestication. This study will further broaden our insights into the evolution and functional elucidation of AQP gene family in cotton.

Aquaporin 1 is located on the intestinal basolateral membrane in Toxocara canis and might play a role in drug uptake

Background

Aquaporins (AQPs) are a family of integral membrane channel proteins that facilitate the transport of water and other small solutes across cell membranes. AQPs appear to play crucial roles in parasite survival and represent possible drug targets for novel intervention strategy. In this work, we investigated the tissue distribution and biological roles of an aquaporin TcAQP1 in the neglected parasitic nematode Toxocara canis.

Methods

Recombinant C-terminal hydrophilic domain of AQP1 of T. canis (rTcAQP1c) and polyclonal antibody against rTcAQP1c were produced to analyse the tissue expression of native TcAQP1 in adult (female and male) worms using an immunohistochemical approach. RNA interference (RNAi), quantitative real-time PCR (qRT-PCR) and nematocidal assays were performed to investigate the functional roles of TcAQP1 in the adult stage of T. canis.

Results

Immunofluorescence analysis showed that TcAQP1 was localised predominantly in the epithelial linings of the reproductive tract and basolateral membrane of the intestine in the adult stage (female and male) of T. canis, indicating important roles in reproduction, nutrient absorption and/or osmoregulation. Treatment with silencing RNA for 24 h resulted in a significant reduction of Tc-aqp-1 mRNA level in adult T. canis, though no phenotypical change was observed. The efficient gene knockdown compromised the nematocidal activity of albendazole in vitro, suggesting the role of TcAQP1 in drug uptake.

Conclusions

The findings of this study provide important information about tissue expression and functional roles of TcAQP1 protein in adult T. canis. Understanding the biological functions of this protein in other developmental stages of T. canis and related parasitic nematodes would contribute to the discovery of novel diagnostic or anthelmintic targets.

The light-induced transcription factor FtMYB116 promotes accumulation of rutin in Fagopyrum tataricum

Tartary buckwheat (Fagopyrum tataricum) not only provides a supplement to primary grain crops in China but also has high medicinal value, by virtue of its rich content of flavonoids possessing antioxidant, anti-inflammatory, and anticancer properties. Light is an important environmental factor that can regulate the synthesis of plant secondary metabolites. In this study, we treated tartary buckwheat seedlings with different wavelengths of light and found that red and blue light could increase the content of flavonoids and the expression of genes involved in flavonoid synthesis pathways. Through coexpression analysis, we identified a new MYB transcription factor (FtMYB116) that can be induced by red and blue light. Yeast one-hybrid assays and an electrophoretic mobility shift assay showed that FtMYB116 binds directly to the promoter region of flavonoid-3'-hydroxylase (F3'H), and a transient luciferase activity assay indicated that FtMYB116 can induce F3'H expression. After transforming FtMYB116 into the hairy roots of tartary buckwheat, we observed significant increases in the content of rutin and quercetin. Collectively, our results indicate that red and blue light promote an increase in flavonoid content in tartary buckwheat seedlings; we also identified a new MYB transcription factor, FtMYB116, that promotes the accumulation of rutin via direct activation of F3'H expression.

Identification and characterization of aquaporin genes in Arachis duranensis and Arachis ipaensis genomes, the diploid progenitors of peanut

BACKGROUND

Aquaporins (AQPs) facilitate transport of water and small solutes across cell membranes and play an important role in different physiological processes in plants. Despite their importance, limited data is available about AQP distribution and function in the economically important oilseed crop peanut, Arachis hypogea (AABB). The present study reports the identification and structural and expression analysis of the AQPs found in the diploid progenitor genomes of A. hypogea i.e. Arachis duranensis (AA) and Arachis ipaensis (BB).

RESULTS

Genome-wide analysis revealed the presence of 32 and 36 AQPs in A. duranensis and A. ipaensis, respectively. Phylogenetic analysis showed similar numbers of AQPs clustered in five distinct subfamilies including the plasma membrane intrinsic proteins (PIPs), the tonoplast intrinsic proteins (TIPs), the nodulin 26-like intrinsic proteins (NIPs), the small basic intrinsic proteins (SIPs), and the uncharacterized intrinsic proteins (XIPs). A notable exception was the XIP subfamily where XIP1 group was observed only in A. ipaensis genome. Protein structure evaluation showed a hydrophilic aromatic/arginine (ar/R) selectivity filter (SF) in PIPs whereas other subfamilies mostly contained a hydrophobic ar/R SF. Both genomes contained one NIP2 with a GSGR SF indicating a conserved ability within the genus to uptake silicon. Analysis of RNA-seq data from A. hypogea revealed a similar expression pattern for the different AQP paralogs of AA and BB genomes. The TIP3s showed seed-specific expression while the NIP1s' expression was confined to roots and root nodules.

CONCLUSIONS

The identification and the phylogenetic analysis of AQPs in both Arachis species revealed the presence of all five sub-families of AQPs. Within the NIP subfamily, the presence of a NIP2 in both genomes supports a conserved ability to absorb Si within plants of the genus. The global expression profile of AQPs in A. hypogea revealed a similar pattern of AQP expression regardless of the subfamilies or the genomes. The tissue-specific expression of AQPs suggests an important role in the development and function of the respective organs. The AQPs identified in the present study will serve as a resource for further characterization and possible exploitation of AQPs to understand their physiological role in A. hypogea.

Investigation of the AQP Family in Soybean and the Promoter Activity of TIP2;6 in Heat Stress and Hormone Responses

Aquaporins (AQPs) are one diverse family of membrane channel proteins that play crucial regulatory roles in plant stress physiology. However, the heat stress responsiveness of AQP genes in soybean remains poorly understood. In this study, 75 non-redundant AQP encoding genes were identified in soybean. Multiple sequence alignments showed that all GmAQP proteins possessed the conserved regions, which contained 6 trans-membrane domains (TM1 to TM6). Different GmAQP members consisted of distinct Asn-Pro-Ala (NPA) motifs, aromatic/arginine (ar/R) selectivity filters and Froger's positions (FPs). Phylogenetic analyses distinguished five sub-families within these GmAQPs: 24 GmPIPs, 24 GmTIPs, 17 GmNIPs, 8 GmSIPs, and 2 GmXIPs. Promoter cis-acting elements analyses revealed that distinct number and composition of heat stress and hormone responsive elements existed in different promoter regions of GmAQPs. QRT-PCR assays demonstrated that 12 candidate GmAQPs with relatively extensive expression in various tissues or high expression levels in root or leaf exhibited different expression changes under heat stress and hormone cues (abscisic acid (ABA), l-aminocyclopropane-l-carboxylic acid (ACC), salicylic acid (SA) and methyl jasmonate (MeJA)). Furthermore, the promoter activity of one previously functionally unknown AQP gene-GmTIP2;6 was investigated in transgenic Arabidopsis plants. The beta-glucuronidase (GUS) activity driven by the promoter of GmTIP2;6 was strongly induced in the heat- and ACC-treated transgenic plants and tended to be accumulated in the hypocotyls, vascular bundles, and leaf trichomes. These results will contribute to uncovering the potential functions and molecular mechanisms of soybean GmAQPs in mediating heat stress and hormone signal responses.

Biomimetic preparation of silicon quantum dots and their phytophysiology effect on cucumber seedlings

In this study, a biomimetic synthetic strategy was proposed for a facile preparation of red fluorescent silicon quantum dots (SiQDs) using unicellular algae of diatoms as reaction precursor.

Whole-Genome and Expression Analyses of Bamboo Aquaporin Genes Reveal Their Functions Involved in Maintaining Diurnal Water Balance in Bamboo Shoots

Water supply is essential for maintaining normal physiological function during the rapid growth of bamboo. Aquaporins (AQPs) play crucial roles in water transport for plant growth and development. Although 26 PeAQPs in bamboo have been reported, the aquaporin-led mechanism of maintaining diurnal water balance in bamboo shoots remains unclear. In this study, a total of 63 PeAQPs were identified, based on the updated genome of moso bamboo (Phyllostachys edulis), including 22 PePIPs, 20 PeTIPs, 17 PeNIPs, and 4 PeSIPs. All of the PeAQPs were differently expressed in 26 different tissues of moso bamboo, based on RNA sequencing (RNA-seq) data. The root pressure in shoots showed circadian rhythm changes, with positive values at night and negative values in the daytime. The quantitative real-time PCR (qRT-PCR) result showed that 25 PeAQPs were detected in the base part of the shoots, and most of them demonstrated diurnal rhythm changes. The expression levels of some PeAQPs were significantly correlated with the root pressure. Of the 86 sugar transport genes, 33 had positive co-expression relationships with 27 PeAQPs. Two root pressure-correlated PeAQPs, PeTIP4;1 and PeTIP4;2, were confirmed to be highly expressed in the parenchyma and epidermal cells of bamboo culm, and in the epidermis, pith, and primary xylem of bamboo roots by in situ hybridization. The authors’ findings provide new insights and a possible aquaporin-led mechanism for bamboo fast growth.

Arbuscular Mycorrhizal Symbiosis Leads to Differential Regulation of Drought-Responsive Genes in Tissue-Specific Root Cells of Common Bean

Arbuscular mycorrhizal fungi (AMF) colonization in plants promotes both local and systemic changes in the gene expression profiles of the host that might be relevant for drought-stress perception and response. Drought-tolerant common bean plants (cv. BAT 477), colonized by a mixture of AMF (Glomus clarum, Acaulospora scrobiculata, and Gigaspora rosea), were exposed to a water deprivation regime of 96 h during pre-flowering. Root transcriptomes were accessed through RNA-Seq revealing a set of 9,965 transcripts with significant differential regulation in inoculated plants during a water deficit event, and 10,569 in non-inoculated. These data include 1,589 transcripts that are exclusively regulated by AMF-inoculation, and 2,313 under non-inoculation conditions. Relative gene expression analyses of nine aquaporin-related transcripts were performed in roots and leaves of plants harvested at initial stages of treatment. Significant shifts in gene expression were detected in AM water deficit-treated roots, in relation to non-inoculated, between 48 and 72 h. Leaves also showed significant mycorrhizal influence in gene expression, especially after 96 h. Root cortical cells, harboring or not arbuscules, were collected from both inoculation treatments through a laser microdissection-based technique. This allowed the identification of transcripts, such as the aquaporin PvPIP2;3 and Glucan 1,3 β-Glucosidase, that are unique to arbuscule-containing cells. During the water deficit treatment, AMF colonization exerted a fine-tune regulation in the expression of genes in the host. That seemed to initiate in arbuscule-containing cells and, as the stressful condition persisted, propagated to the whole-plant through secondary signaling events. Collectively, these results demonstrate that arbuscular mycorrhization leads to shifts in common bean's transcriptome that could potentially impact its adaptation capacity during water deficit events.

Red and Blue Light Promote the Accumulation of Artemisinin in Artemisia Annua L

Artemisinin, which has been isolated from Artemisia annua L., is the most effective antimalarial drug and has saved millions of lives. In addition, artemisinin and its derivatives have anti-tumor, anti-parasitic, anti-fibrosis, and anti-arrhythmic properties, which enhances the demand for these compounds. Improving the content of artemisinin in A. annua is therefore becoming an increasing research interest, as the chemical synthesis of this metabolite is not viable. Ultraviolet B and C irradiation have been reported to improve the artemisinin content in A. annua, but they are harmful to plant growth and development. Therefore, we screened other light sources to examine if they could promote artemisinin content without affecting plant growth and development. We found that red and blue light could enhance artemisinin accumulation by promoting the expression of the genes that were involved in artemisinin biosynthesis, such as amorpha-4,11-diene synthase (ADS) and cytochrome P450 monooxygenase (CYP71AV1) genes. Thus, in addition to being the main light sources for photosynthesis, red and blue light play a key role in plant secondary metabolism, and optimizing the combination of these light might allow for the productionof artemisinin-rich A. annua.

Identification and substrate prediction of new Fragaria x ananassa aquaporins and expression in different tissues and during strawberry fruit development

The newly identified aquaporin coding sequences presented here pave the way for further insights into the plant-water relations in the commercial strawberry (Fragaria x ananassa). Aquaporins are water channel proteins that allow water to cross (intra)cellular membranes. In Fragaria x ananassa, few of them have been identified hitherto, hampering the exploration of the water transport regulation at cellular level. Here, we present new aquaporin coding sequences belonging to different subclasses: plasma membrane intrinsic proteins subtype 1 and subtype 2 (PIP1 and PIP2) and tonoplast intrinsic proteins (TIP). The classification is based on phylogenetic analysis and is confirmed by the presence of conserved residues. Substrate-specific signature sequences (SSSSs) and specificity-determining positions (SDPs) predict the substrate specificity of each new aquaporin. Expression profiling in leaves, petioles and developing fruits reveals distinct patterns, even within the same (sub)class. Expression profiles range from leaf-specific expression over constitutive expression to fruit-specific expression. Both upregulation and downregulation during fruit ripening occur. Substrate specificity and expression profiles suggest that functional specialization exists among aquaporins belonging to a different but also to the same (sub)class.

RNA Sequencing and Coexpression Analysis Reveal Key Genes Involved in α-Linolenic Acid Biosynthesis in Perilla frutescens Seed

Perilla frutescen is used as traditional food and medicine in East Asia. Its seeds contain high levels of α-linolenic acid (ALA), which is important for health, but is scarce in our daily meals. Previous reports on RNA-seq of perilla seed had identified fatty acid (FA) and triacylglycerol (TAG) synthesis genes, but the underlying mechanism of ALA biosynthesis and its regulation still need to be further explored. So we conducted Illumina RNA-sequencing in seven temporal developmental stages of perilla seeds. Sequencing generated a total of 127 million clean reads, containing 15.88 Gb of valid data. The de novo assembly of sequence reads yielded 64,156 unigenes with an average length of 777 bp. A total of 39,760 unigenes were annotated and 11,693 unigenes were found to be differentially expressed in all samples. According to Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis, 486 unigenes were annotated in the “lipid metabolism” pathway. Of these, 150 unigenes were found to be involved in fatty acid (FA) biosynthesis and triacylglycerol (TAG) assembly in perilla seeds. A coexpression analysis showed that a total of 104 genes were highly coexpressed (r > 0.95). The coexpression network could be divided into two main subnetworks showing over expression in the medium or earlier and late phases, respectively. In order to identify the putative regulatory genes, a transcription factor (TF) analysis was performed. This led to the identification of 45 gene families, mainly including the AP2-EREBP, bHLH, MYB, and NAC families, etc. After coexpression analysis of TFs with highly expression of FAD2 and FAD3 genes, 162 TFs were found to be significantly associated with two FAD genes (r > 0.95). Those TFs were predicted to be the key regulatory factors in ALA biosynthesis in perilla seed. The qRT-PCR analysis also verified the relevance of expression pattern between two FAD genes and partial candidate TFs. Although it has been reported that some TFs are involved in seed development, more direct evidence is still needed to verify their function. However, these findings can provide clues to reveal the possible molecular mechanisms of ALA biosynthesis and its regulation in perilla seed.

Prediction of arsenic and antimony transporter major intrinsic proteins from the genomes of crop plants

Major intrinsic proteins (MIPs), commonly known as aquaporins, transport water and non-polar small solutes. Comparing the 3D models and the primary selectivity-related motifs (two Asn-Pro-Ala (NPA) regions, the aromatic/arginine (ar/R) selectivity filter, and Froger's positions (FPs)) of all plant MIPs that have been experimentally proven to transport arsenic (As) and antimony (Sb), some substrate-specific signature sequences (SSSS) or specificity determining sites (SDPs) have been predicted. These SSSS or SDPs were determined in 543 MIPs found in the genomes of 12 crop plants; the As and Sb transporters were predicted to be distributed in noduline-26 like intrinsic proteins (NIPs), and every plant had one or several As and Sb transporter NIPs. Phylogenetic grouping of the NIP subfamily based on the ar/R selectivity filter and FPs were linked to As and Sb transport. We further determined the group-wise substrate selectivity profiles of the NIPs in the 12 crop plants. In addition to two NPA regions, the ar/R filter, and FPs, certain amino acids especially in the pore line, loop D, and termini contribute to the functional distinctiveness of the NIP groups. Expression analysis of transcripts in different organs indicated that most of the As and Sb transporter NIPs were expressed in roots.

Genome-wide identification and characterization of aquaporin gene family in Beta vulgaris

Aquaporins (AQPs) are essential channel proteins that execute multi-functions throughout plant growth and development, including water transport, uncharged solutes uptake, stress response, and so on. Here, we report the first genome-wide identification and characterization AQP (BvAQP) genes in sugar beet (Beta vulgaris), an important crop widely cultivated for feed, for sugar production and for bioethanol production. Twenty-eight sugar beet AQPs (BvAQPs) were identified and assigned into five subfamilies based on phylogenetic analyses: seven of plasma membrane (PIPs), eight of tonoplast (TIPs), nine of NOD26-like (NIPs), three of small basic (SIPs), and one of x-intrinsic proteins (XIPs). BvAQP genes unevenly mapped on all chromosomes, except on chromosome 4. Gene structure and motifs analyses revealed that BvAQP have conserved exon-intron organization and that they exhibit conserved motifs within each subfamily. Prediction of BvAQPs functions, based on key protein domains conservation, showed a remarkable difference in substrate specificity among the five subfamilies. Analyses of BvAQPs expression, by mean of RNA-seq, in different plant organs and in response to various abiotic stresses revealed that they were ubiquitously expressed and that their expression was induced by heat and salt stresses. These results provide a reference base to address further the function of sugar beet aquaporins and to explore future applications for plants growth and development improvements as well as in response to environmental stresses.

Carbon dioxide and water transport through plant aquaporins

Aquaporins are channel proteins that function to increase the permeability of biological membranes. In plants, aquaporins are encoded by multigene families that have undergone substantial diversification in land plants. The plasma membrane intrinsic proteins (PIPs) subfamily of aquaporins is of particular interest given their potential to improve plant water relations and photosynthesis. Flowering plants have between 7 and 28 PIP genes. Their expression varies with tissue and cell type, through development and in response to a variety of factors, contributing to the dynamic and tissue specific control of permeability. There are a growing number of PIPs shown to act as water channels, but those altering membrane permeability to CO₂ are more limited. The structural basis for selective substrate specificities has not yet been resolved, although a few key amino acid positions have been identified. Several regions important for dimerization, gating and trafficking are also known. PIP aquaporins assemble as tetramers and their properties depend on the monomeric composition. PIPs control water flux into and out of veins and stomatal guard cells and also increase membrane permeability to CO₂ in mesophyll and stomatal guard cells. The latter increases the effectiveness of Rubisco and can potentially influence transpiration efficiency.

Characterization and Regulation of Aquaporin Genes of Sorghum [Sorghum bicolor (L.) Moench] in Response to Waterlogging Stress

Waterlogging is a significant environmental constraint to crop production, and a better understanding of plant responses is critical for the improvement of crop tolerance to waterlogged soils. Aquaporins (AQPs) are a class of channel-forming proteins that play an important role in water transport in plants. This study aimed to examine the regulation of AQP genes under waterlogging stress and to characterize the genetic variability of AQP genes in sorghum (Sorghum bicolor). Transcriptional profiling of AQP genes in response to waterlogging stress in nodal root tips and nodal root basal regions of two tolerant and two sensitive sorghum genotypes at 18 and 96 h after waterlogging stress imposition revealed significant gene-specific pattern with regard to genotype, root tissue sample, and time point. For some tissue sample and time point combinations, PIP2-6, PIP2-7, TIP2-2, TIP4-4, and TIP5-1 expression was differentially regulated in tolerant compared to sensitive genotypes. The differential response of these AQP genes suggests that they may play a tissue specific role in mitigating waterlogging stress. Genetic analysis of sorghum revealed that AQP genes were clustered into the same four subfamilies as in maize (Zea mays) and rice (Oryza sativa) and that residues determining the AQP channel specificity were largely conserved across species. Single nucleotide polymorphism (SNP) data from 50 sorghum accessions were used to build an AQP gene-based phylogeny of the haplotypes. Phylogenetic analysis based on single nucleotide polymorphisms of sorghum AQP genes placed the tolerant and sensitive genotypes used for the expression study in distinct groups. Expression analyses suggested that selected AQPs may play a pivotal role in sorghum tolerance to water logging stress. Further experimentation is needed to verify their role and to leverage phylogenetic analyses and AQP expression data to improve waterlogging tolerance in sorghum.

Genome-wide identification, characterization, and expression profile of aquaporin gene family in flax (Linum usitatissimum)

Membrane intrinsic proteins (MIPs) form transmembrane channels and facilitate transport of myriad substrates across the cell membrane in many organisms. Majority of plant MIPs have water transporting ability and are commonly referred as aquaporins (AQPs). In the present study, we identified aquaporin coding genes in flax by genome-wide analysis, their structure, function and expression pattern by pan-genome exploration. Cross-genera phylogenetic analysis with known aquaporins from rice, arabidopsis, and poplar showed five subgroups of flax aquaporins representing 16 plasma membrane intrinsic proteins (PIPs), 17 tonoplast intrinsic proteins (TIPs), 13 NOD26-like intrinsic proteins (NIPs), 2 small basic intrinsic proteins (SIPs), and 3 uncharacterized intrinsic proteins (XIPs). Amongst aquaporins, PIPs contained hydrophilic aromatic arginine (ar/R) selective filter but TIP, NIP, SIP and XIP subfamilies mostly contained hydrophobic ar/R selective filter. Analysis of RNA-seq and microarray data revealed high expression of PIPs in multiple tissues, low expression of NIPs, and seed specific expression of TIP3 in flax. Exploration of aquaporin homologs in three closely related Linum species bienne, grandiflorum and leonii revealed presence of 49, 39 and 19 AQPs, respectively. The genome-wide identification of aquaporins, first in flax, provides insight to elucidate their physiological and developmental roles in flax.

Aquaporins as potential drought tolerance inducing proteins: Towards instigating stress tolerance

Aquaporins (AQPs) are primarily involved in maintaining cellular water homeostasis. Their role in diverse physiological processes has fascinated plant scientists for more than a decade, particularly concerning abiotic stresses. Increasing examples of evidence in various crop plants indicate that the AQPs are responsible for precise regulation of water movement and consequently play a crucial role in the drought stress tolerance. Since drought is one of the major abiotic stresses affecting agricultural production worldwide, it has become a critical agenda to focus research on the development of drought tolerant crop plants. AQPs can act as key candidate molecules to confront this issue. Hence, there is an important need to explore the potential of AQPs by understanding the molecular mechanisms and pathways through which they induce drought tolerance. Moreover, the signalling network/s involved in such pathways needs to be mined and understood correctly, and that may lead to the development of drought tolerance in crop plants. In the present review, opportunity and challenges regarding the efficient utilization of AQP-related information is presented and discussed. The complied information and the discussion will be helpful for designing future experiments and to set the specific goals for the enhancement of drought tolerance in crop plants. Biological Significance Knowledge on the role of AQPs in maintaining cellular water homeostasis has given new hope for developing drought tolerance in crop plants. Since drought is one of the major abiotic stresses affecting agricultural production worldwide, it has become a critical agenda to focus research on the development of drought-tolerant crop plants. AQPs can act as key candidate molecules to solve this problem through genetic engineering. For this, it is important to understand the molecular mechanisms and inter-related pathways through which AQPs induce drought tolerance and to explore the signaling network/s involved in such pathways.

Single amino acid substitutions in the selectivity filter render NbXIP1;1α aquaporin water permeable

BACKGROUND

Aquaporins (AQPs) are integral membrane proteins that facilitate transport of water and/or other small neutral solutes across membranes in all forms of life. The X Intrinsic Proteins (XIPs) are the most recently recognized and the least characterized aquaporin subfamily in higher plants. XIP1s have been shown to be impermeable to water but permeable to boric acid, glycerol, hydrogen peroxide and urea. However, uncertainty regarding the determinants for selectivity and lack of an activity that is easy to quantify have hindered functional investigations. In an effort to resolve these issues, we set out to introduce water permeability in Nicotiana benthamiana XIP1;1α (NbXIP1;1α), by exchanging amino acid residues of predicted alternative aromatic/arginine (ar/R) selectivity filters of NbXIP1;1α for residues constituting the water permeable ar/R selectivity filter of AtTIP2;1.

RESULTS

Here, we present functional results regarding the amino acid substitutions in the putative filters as well as deletions in loops C and D of NbXIP1;1α. In addition, homology models were created based on the high resolution X-ray structure of AtTIP2;1 to rationalize the functional properties of wild-type and mutant NbXIP1;1α. Our results favour Thr 246 rather than Val 242 as the residue at the helix 5 position in the ar/R filter of NbXIP1;1α and indicate that the pore is not occluded by the loops when heterologously expressed in Pichia pastoris. Moreover, our results show that a single amino acid substitution in helix 1 (L79G) or in helix 2 (I102H) is sufficient to render NbXIP1;1α water permeable. Most of the functional results can be rationalized from the models based on a combination of aperture and hydrophobicity of the ar/R filter.

CONCLUSION

The water permeable NbXIP1;1α mutants imply that the heterologously expressed proteins are correctly folded and offer means to explore the structural and functional properties of NbXIP1;1α. Our results support that Thr 246 is part of the ar/R filter. Furthermore, we suggest that a salt bridge to an acidic residue in helix 1, conserved among the XIPs in clade B, directs the orientation of the arginine in the ar/R selectivity filter and provides a novel approach to tune the selectivity of AQPs.

Genome-Wide Identification and Characterization of the Aquaporin Gene Family and Transcriptional Responses to Boron Deficiency in Brassica napus

Aquaporins (AQPs) are an abundant protein family and play important roles to facilitate small neutral molecule transport across membranes. Oilseed rape (Brassica napus L.) is an important oil crop in China and elsewhere in the world, and is very sensitive to low boron (B) stress. Several AQP family genes have been reported to be involved in B transport across plasma membranes in plants. In this study, a total of 121 full-length AQPs were identified and characterized in B. napus (AC genome), and could be classified into four sub-families, including 43 PIPs (plasma membrane intrinsic proteins), 35 TIPs (tonoplast intrinsic proteins), 32 NIPs (NOD26-like intrinsic proteins), and 11 SIPs (small basic intrinsic proteins). The gene characteristics of BnaAQPs were similar to those of BraAQPs (A genome) and BolAQPs (C genome) including the composition of each sub-family, gene structure, and substrate selectivity filters. The BnaNIP was the most complex AQP sub-family, reflecting the composition of substrate selectivity filter structures which affect the permeation of solution molecules. In this study, the seedlings of both B-efficient (QY10) and B-inefficient (W10) cultivars were treated with two boron (B) levels: deficient (0.25 μM B) and sufficient (25 μM B). The transcription of AQP genes in root (R), juvenile leaf (JL), and old leaf (OL) tissues of both cultivars was investigated under B deficient and sufficient conditions. Transcription of most BnaPIPs and BnaTIPs was significantly increased compared with other BnaAQPs in all the three tissues, especially in the roots, of both B-efficient and B-inefficient cultivars under both B conditions. With B deprivation, the expression of the majority of the BnaPIPs and BnaTIPs was down-regulated in the roots. However, the BnaNIPs were up-regulated. In addition, the BnaCnn_random.PIP1;4b, BnaPIP2;4s, BnaC04.TIP4;1a, BnaAnn_random.TIP1;1b, and BnaNIP5;1s (except for BnaA07.NIP5;1c and BnaC06.NIP5;1c) exhibited obvious differences at low B between B-efficient and B-inefficient cultivars. These results will help us to understand boron homeostasis in B. napus.

Characterization and differential expression analysis of Toxocara canis aquaporin-1 gene

Toxocara canis is an intestinal nematode of canids with a worldwide distribution, causing an important but neglected parasitic zoonosis in humans. Aquaporins (AQP) are a family of water channel proteins, which function as membrane channels to regulate water homeostasis. In this study, the coding sequence of aquaporin-1 gene of T. canis (Tc-aqp-1) was cloned and characterized. The obtained Tc-aqp-1 coding sequence was 933 bp in length, which predicted to encode 311 amino acids. Two conserved asparagine-proline-alanine (NPA) motifs were identified in the multiple sequence alignments. Phylogenetic analysis revealed the closest relationship between T. canis and Opisthorchis viverrini based on aquaporin-1 amino acid sequence. A structure was predicted with ligand binding sites predicted at H93, N95, N226, L94, I79, and I210 and with active sites predicted at I256 and G207. Gene Ontology (GO) annotations predicted its cellular component term of integral component of plasma membrane (GO: 0005887), molecular function term of channel activity (GO: 0015250), and biological process term of water transport (GO: 0006833). Tissue expression analysis revealed that the Tc-aqp-1 was highly expressed in the intestine of adult male. The findings of the present study provide the basis for further functional studies of T. canis aquaporin-1.

Genome-Wide Identification of Jatropha curcas Aquaporin Genes and the Comparative Analysis Provides Insights into the Gene Family Expansion and Evolution in Hevea brasiliensis

Aquaporins (AQPs) are channel-forming integral membrane proteins that transport water and other small solutes across biological membranes. Despite the vital role of AQPs, to date, little is known in physic nut (Jatropha curcas L., Euphorbiaceae), an important non-edible oilseed crop with great potential for the production of biodiesel. In this study, 32 AQP genes were identified from the physic nut genome and the family number is relatively small in comparison to 51 in another Euphorbiaceae plant, rubber tree (Hevea brasiliensis Muell. Arg.). Based on the phylogenetic analysis, the JcAQPs were assigned to five subfamilies, i.e., nine plasma membrane intrinsic proteins (PIPs), nine tonoplast intrinsic proteins (TIPs), eight NOD26-like intrinsic proteins (NIPs), two X intrinsic proteins (XIPs), and four small basic intrinsic proteins (SIPs). Like rubber tree and other plant species, functional prediction based on the aromatic/arginine selectivity filter, Froger's positions, and specificity-determining positions showed a remarkable difference in substrate specificity among subfamilies of JcAQPs. Genome-wide comparative analysis revealed the specific expansion of PIP and TIP subfamilies in rubber tree and the specific gene loss of the XIP subfamily in physic nut. Furthermore, by analyzing deep transcriptome sequencing data, the expression evolution especially the expression divergence of duplicated HbAQP genes was also investigated and discussed. Results obtained from this study not only provide valuable information for future functional analysis and utilization of Jc/HbAQP genes, but also provide a useful reference to survey the gene family expansion and evolution in Euphorbiaceae plants and other plant species.

Genome-wide identification and characterization of aquaporin gene family in moso bamboo (Phyllostachys edulis)

Aquaporins (AQPs) are known to play a major role in maintaining water and hydraulic conductivity balance in the plant system. Numerous studies have showed AQPs execute multi-function throughout plant growth and development, including water transport, nitrogen, carbon, and micronutrient acquisition etc. However, little information on AQPs is known in bamboo. In this study, we present the first genome-wide identification and characterization of AQP genes in moso bamboo (Phyllostachys edulis) using bioinformatics. In total, 26 AQP genes were identified by homologous analysis, which were divided into four groups (PIPs, TIPs, NIPs, and SIPs) based on the phylogenetic analysis. All the genes were located on 26 different scaffolds respectively on basis of the gene mapped to bamboo genome. Evolutionary analysis indicated that Ph. edulis was more close to Oryza sativa than Zea mays in the genetic relationship. Besides, qRT-PCR was used to analyze gene expression profiles, which revealed that AQP genes were expressed constitutively in all the detected tissues, and were all responsive to the environmental cues such as drought, water, and NaCl stresses. This data suggested that AQPs may play fundamental roles in maintaining normal growth and development of bamboo, which would contribute to better understanding for the complex regulation mechanism involved in the fast-growing process of bamboo. Furthermore, the result could provide valuable information for further research on bamboo functional genomics.

Plant Aquaporins: Genome-Wide Identification, Transcriptomics, Proteomics, and Advanced Analytical Tools

Aquaporins (AQPs) are channel-forming integral membrane proteins that facilitate the movement of water and many other small molecules. Compared to animals, plants contain a much higher number of AQPs in their genome. Homology-based identification of AQPs in sequenced species is feasible because of the high level of conservation of protein sequences across plant species. Genome-wide characterization of AQPs has highlighted several important aspects such as distribution, genetic organization, evolution and conserved features governing solute specificity. From a functional point of view, the understanding of AQP transport system has expanded rapidly with the help of transcriptomics and proteomics data. The efficient analysis of enormous amounts of data generated through omic scale studies has been facilitated through computational advancements. Prediction of protein tertiary structures, pore architecture, cavities, phosphorylation sites, heterodimerization, and co-expression networks has become more sophisticated and accurate with increasing computational tools and pipelines. However, the effectiveness of computational approaches is based on the understanding of physiological and biochemical properties, transport kinetics, solute specificity, molecular interactions, sequence variations, phylogeny and evolution of aquaporins. For this purpose, tools like Xenopus oocyte assays, yeast expression systems, artificial proteoliposomes, and lipid membranes have been efficiently exploited to study the many facets that influence solute transport by AQPs. In the present review, we discuss genome-wide identification of AQPs in plants in relation with recent advancements in analytical tools, and their availability and technological challenges as they apply to AQPs. An exhaustive review of omics resources available for AQP research is also provided in order to optimize their efficient utilization. Finally, a detailed catalog of computational tools and analytical pipelines is offered as a resource for AQP research.

The Aquaporin Splice Variant NbXIP1;1α Is Permeable to Boric Acid and Is Phosphorylated in the N-terminal Domain

Aquaporins (AQPs) are membrane channel proteins that transport water and uncharged solutes across different membranes in organisms in all kingdoms of life. In plants, the AQPs can be divided into seven different subfamilies and five of these are present in higher plants. The most recently characterized of these subfamilies is the XIP subfamily, which is found in most dicots but not in monocots. In this article, we present data on two different splice variants (α and β) of NbXIP1;1 from Nicotiana benthamiana. We describe the heterologous expression of NbXIP1;1α and β in the yeast Pichia pastoris, the subcellular localization of the protein in this system and the purification of the NbXIP1;1α protein. Furthermore, we investigated the functionality and the substrate specificity of the protein by stopped-flow spectrometry in P. pastoris spheroplasts and with the protein reconstituted in proteoliposomes. The phosphorylation status of the protein and localization of the phosphorylated amino acids were verified by mass spectrometry. Our results show that NbXIP1;1α is located in the plasma membrane when expressed in P. pastoris, that it is not permeable to water but to boric acid and that the protein is phosphorylated at several amino acids in the N-terminal cytoplasmic domain of the protein. A growth assay showed that the yeast cells expressing the N-terminally His-tagged NbXIP1;1α were more sensitive to boric acid as compared to the cells expressing the C-terminally His-tagged isoform. This might suggest that the N-terminal His-tag functionally mimics the phosphorylation of the N-terminal domain and that the N-terminal domain is involved in gating of the channel.

Genome-Wide Analysis of the Aquaporin Gene Family in Chickpea (Cicer arietinum L.)

Aquaporins (AQPs) are essential membrane proteins that play critical role in the transport of water and many other solutes across cell membranes. In this study, a comprehensive genome-wide analysis identified 40 AQP genes in chickpea (Cicer arietinum L.). A complete overview of the chickpea AQP (CaAQP) gene family is presented, including their chromosomal locations, gene structure, phylogeny, gene duplication, conserved functional motifs, gene expression, and conserved promoter motifs. To understand AQP's evolution, a comparative analysis of chickpea AQPs with AQP orthologs from soybean, Medicago, common bean, and Arabidopsis was performed. The chickpea AQP genes were found on all of the chickpea chromosomes, except chromosome 7, with a maximum of six genes on chromosome 6, and a minimum of one gene on chromosome 5. Gene duplication analysis indicated that the expansion of chickpea AQP gene family might have been due to segmental and tandem duplications. CaAQPs were grouped into four subfamilies including 15 NOD26-like intrinsic proteins (NIPs), 13 tonoplast intrinsic proteins (TIPs), eight plasma membrane intrinsic proteins (PIPs), and four small basic intrinsic proteins (SIPs) based on sequence similarities and phylogenetic position. Gene structure analysis revealed a highly conserved exon-intron pattern within CaAQP subfamilies supporting the CaAQP family classification. Functional prediction based on conserved Ar/R selectivity filters, Froger's residues, and specificity-determining positions suggested wide differences in substrate specificity among the subfamilies of CaAQPs. Expression analysis of the AQP genes indicated that some of the genes are tissue-specific, whereas few other AQP genes showed differential expression in response to biotic and abiotic stresses. Promoter profiling of CaAQP genes for conserved cis-acting regulatory elements revealed enrichment of cis-elements involved in circadian control, light response, defense and stress responsiveness reflecting their varying pattern of gene expression and potential involvement in biotic and abiotic stress responses. The current study presents the first detailed genome-wide analysis of the AQP gene family in chickpea and provides valuable information for further functional analysis to infer the role of AQP in the adaptation of chickpea in diverse environmental conditions.