Genome-wide identification of rubber tree (Hevea brasiliensis Muell. Arg.) aquaporin genes and their response to ethephon stimulation in the laticifer, a rubber-producing tissue

TIP aquaporins in Cyperus esculentus: genome-wide identification, expression profiles, subcellular localizations, and interaction patterns

BACKGROUND

Tonoplast intrinsic proteins (TIPs), which typically mediate water transport across vacuolar membranes, play an essential role in plant growth, development, and stress responses. However, their characterization in tigernut (Cyperus esculentus L.), an oil-bearing tuber plant of the Cyperaceae family, is still in the infancy.

RESULTS

In this study, a first genome-wide characterization of the TIP subfamily was conducted in tigernut, resulting in ten members representing five previously defined phylogenetic groups, i.e., TIP1-5. Although the gene amounts are equal to that present in two model plants Arabidopsis and rice, the group composition and/or evolution pattern were shown to be different. Except for CeTIP1;3 that has no counterpart in both Arabidopsis and rice, complex orthologous relationships of 1:1, 1:2, 1:3, 2:1, and 2:2 were observed. Expansion of the CeTIP subfamily was contributed by whole-genome duplication (WGD), transposed, and dispersed duplications. In contrast to the recent WGD-derivation of CeTIP3;1/-3;2, synteny analyses indicated that TIP4 and - 5 are old WGD repeats of TIP2, appearing sometime before monocot-eudicot divergence. Expression analysis revealed that CeTIP genes exhibit diverse expression profiles and are subjected to developmental and diurnal fluctuation regulation. Moreover, when transiently overexpressed in tobacco leaves, CeTIP1;1 was shown to locate in the vacuolar membrane and function in homo/heteromultimer, whereas CeTIP2;1 is located in the cell membrane and only function in heteromultimer. Interestingly, CeTIP1;1 could mediate the tonoplast-localization of CeTIP2;1 via protein interaction, implying complex regulatory patterns.

CONCLUSIONS

Our findings provide a global view of CeTIP genes, which provide valuable information for further functional analysis and genetic improvement through manipulating key members in tigernut.

Integrative Analysis of Oleosin Genes Provides Insights into Lineage-Specific Family Evolution in Brassicales

Oleosins (OLEs) are a class of small but abundant structural proteins that play essential roles in the formation and stabilization of lipid droplets (LDs) in seeds of oil crops. Despite the proposal of five oleosin clades (i.e., U, SL, SH, T, and M) in angiosperms, their evolution in eudicots has not been well-established. In this study, we employed Brassicales, an economically important order of flowering plants possessing the lineage-specific T clade, as an example to address this issue. Three to 10 members were identified from 10 species representing eight plant families, which include Caricaceae, Moringaceae, Akaniaceae, Capparaceae, and Cleomaceae. Evolutionary and reciprocal best hit-based homologous analyses assigned 98 oleosin genes into six clades (i.e., U, SL, SH, M, N, and T) and nine orthogroups (i.e., U1, U2, SL, SH1, SH2, SH3, M, N, and T). The newly identified N clade represents an ancient group that has already appeared in the basal angiosperm Amborella trichopoda, which are constitutively expressed in the tree fruit crop Carica papaya, including pulp and seeds of the fruit. Moreover, similar to Clade N, the previously defined M clade is actually not Lauraceae-specific but an ancient and widely distributed group that diverged before the radiation of angiosperm. Compared with A. trichopoda, lineage-specific expansion of the family in Brassicales was largely contributed by recent whole-genome duplications (WGDs) as well as the ancient γ event shared by all core eudicots. In contrast to the flower-preferential expression of Clade T, transcript profiling revealed an apparent seed/embryo/endosperm-predominant expression pattern of most oleosin genes in Arabidopsis thaliana and C. papaya. Moreover, the structure and expression divergence of paralogous pairs was frequently observed, and a good example is the lineage-specific gain of an intron. These findings provide insights into lineage-specific family evolution in Brassicales, which facilitates further functional studies in nonmodel plants such as C. papaya.

Analysis of Carica papaya Informs Lineage-Specific Evolution of the Aquaporin (AQP) Family in Brassicales

Aquaporins (AQPs), a type of intrinsic membrane proteins that transport water and small solutes across biological membranes, play crucial roles in plant growth and development. This study presents a first genome-wide identification and comparative analysis of the AQP gene family in papaya (Carica papaya L.), an economically and nutritionally important fruit tree of tropical and subtropical regions. A total of 29 CpAQP genes were identified, which represent five subfamilies, i.e., nine plasma intrinsic membrane proteins (PIPs), eight tonoplast intrinsic proteins (TIPs), seven NOD26-like intrinsic proteins (NIPs), two X intrinsic proteins (XIPs), and three small basic intrinsic proteins (SIPs). Although the family is smaller than the 35 members reported in Arabidopsis, it is highly diverse, and the presence of CpXIP genes as well as orthologs in Moringa oleifera and Bretschneidera sinensis implies that the complete loss of the XIP subfamily in Arabidopsis is lineage-specific, sometime after its split with papaya but before Brassicaceae-Cleomaceae divergence. Reciprocal best hit-based sequence comparison of 530 AQPs and synteny analyses revealed that CpAQP genes belong to 29 out of 61 identified orthogroups, and lineage-specific evolution was frequently observed in Brassicales. Significantly, the well-characterized NIP3 group was completely lost; lineage-specific loss of the NIP8 group in Brassicaceae occurred sometime before the divergence with Cleomaceae, and lineage-specific loss of NIP2 and SIP3 groups in Brassicaceae occurred sometime after the split with Cleomaceae. In contrast to a predominant role of recent whole-genome duplications (WGDs) on the family expansion in B. sinensis, Tarenaya hassleriana, and Brassicaceae plants, no recent AQP repeats were identified in papaya, and ancient WGD repeats are mainly confined to the PIP subfamily. Subfamily even group-specific evolution was uncovered via comparing exon-intron structures, conserved motifs, the aromatic/arginine selectivity filter, and gene expression profiles. Moreover, down-regulation during fruit ripening and expression divergence of duplicated CpAQP genes were frequently observed in papaya. These findings will not only improve our knowledge on lineage-specific family evolution in Brassicales, but also provide valuable information for further studies of AQP genes in papaya and species beyond.

Molecular characterization of oleosin genes in Cyperus esculentus, a Cyperaceae plant producing oil in underground tubers

KEY MESSAGE

CeOLE genes exhibit a tuber-predominant expression pattern and their mRNA/protein abundances are positively correlated with oil accumulation during tuber development. Overexpression could significantly increase the oil content of tobacco leaves. Oleosins (OLEs) are abundant structural proteins of lipid droplets (LDs) that function in LD formation and stabilization in seeds of oil crops. However, little information is available on their roles in vegetative tissues. In this study, we present the first genome-wide characterization of the oleosin family in tigernut (Cyperus esculentus L., Cyperaceae), a rare example accumulating high amounts of oil in underground tubers. Six members identified represent three previously defined clades (i.e. U, SL and SH) or six out of seven orthogroups (i.e. U, SL1, SL2, and SH1-3) proposed in this study. Comparative genomics analysis reveals that lineage-specific expansion of Clades SL and SH was contributed by whole-genome duplication and dispersed duplication, respectively. Moreover, presence of SL2 and SH3 in Juncus effuses implies their appearance sometime before Cyperaceae-Juncaceae divergence, whereas SH2 appears to be Cyperaceae specific. Expression analysis showed that CeOLE genes exhibit a tuber-predominant expression pattern and transcript levels are considerably more abundant than homologs in the close relative Cyperus rotundus. Moreover, CeOLE mRNA and protein abundances were shown to positively correlate with oil accumulation during tuber development. Additionally, two dominant isoforms (i.e. CeOLE2 and -5) were shown to locate in LDs as well as the endoplasmic reticulum of tobacco (Nicotiana benthamiana) leaves, and are more likely to function in homo and heteromultimers. Furthermore, overexpression of CeOLE2 and -5 in tobacco leaves could significantly increase the oil content, supporting their roles in oil accumulation. These findings provide insights into lineage-specific family evolution and putative roles of CeOLE genes in oil accumulation of vegetative tissues, which facilitate further genetic improvement for tigernut.

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.

Analysis of Lhc family genes reveals development regulation and diurnal fluctuation expression patterns in Cyperus esculentus, a Cyperaceae plant

Sixteen Lhc genes representing 13 phylogenetic groups were identified from the full-length transcriptome of tigernut, exhibiting development regulation and diurnal fluctuation expression patterns in leaves. Nuclear encoded light-harvesting chlorophyll a/b-binding (Lhc) proteins play indispensable roles in oxygenic photosynthesis. In this study, we present the first transcriptome-based characterization of Lhc family genes in tigernut (Cyperus esculentus L.), a Cyperaceae C₄ plant producing oil in underground tubers. A number of 16 Lhc genes representing 13 phylogenetic groups identified from the full-length tigernut transcriptome are equal to that found in both Carex littledalei (another Cyperaceae plant) and papaya, slightly more than 15 members present in both rice and jatropha, but relatively less than 18, 20, and 21 members present in sorghum, cassava, and Arabidopsis, respectively. Nevertheless, nearly one-vs-one orthologous relationship was observed in most groups, though some of them are no longer located in syntenic blocks and species-specific expansion was frequently found in Lhcb1. Comparative genomics analysis revealed that the loss of two groups (i.e., Lhca2 and Lhca5) in C. littledalei is species-specific, sometime after the split with tigernut, and the expansion of Lhcb1 was mainly contributed by tandem duplication as observed in most species. Interestingly, a transposed duplication, which appears to be shared by monocots, was also identified in Lhcb1. Further transcriptome profiling revealed a predominant expression pattern of most CeLhc family genes in photosynthetic tissues and enhanced transcription during leaf maturation, reflecting their key roles in light absorption. Moreover, qRT-PCR analysis revealed an apparent diurnal fluctuation expression pattern of 11 dominant CeLhc genes. These findings not only highlight species-specific evolution of Lhc genes in the Cyperaceae family as well as the monocot lineage, but also provide valuable information for further functional analysis and genetic improvement in tigernut.

Genomic Analysis of LEA Genes in Carica papaya and Insight into Lineage-Specific Family Evolution in Brassicales

Late embryogenesis abundant (LEA) proteins comprise a diverse superfamily involved in plant development and stress responses. This study presents a first genome-wide analysis of LEA genes in papaya (Carica papaya L., Caricaceae), an economically important tree fruit crop widely cultivated in the tropics and subtropics. A total of 28 members were identified from the papaya genome, which belong to eight families with defined Pfam domains, i.e., LEA_1 (3), LEA_2 (4), LEA_3 (5), LEA_4 (5), LEA_5 (2), LEA_6 (2), DHN (4), and SMP (3). The family numbers are comparable to those present in Ricinus communis (Euphorbiaceae, 28) and Moringa oleifera (Moringaceae, 29), but relatively less than that found in Moringa oleifera (Cleomaceae, 39) and Arabidopsis thaliana (Brassicaceae, 51), implying lineage-specific evolution in Brassicales. Indeed, best-reciprocal-hit-based sequence comparison and synteny analysis revealed the presence of 29 orthogroups, and significant gene expansion in Tarenaya and Arabidopsis was mainly contributed by whole-genome duplications that occurred sometime after their split with the papaya. Though a role of transposed duplication was also observed, tandem duplication was shown to be a key contributor in gene expansion of most species examined. Further comparative analyses of exon-intron structures and protein motifs supported fast evolution of this special superfamily, especially in Arabidopsis. Transcriptional profiling revealed diverse expression patterns of CpLEA genes over various tissues and different stages of developmental fruit. Moreover, the transcript level of most genes appeared to be significantly regulated by drought, cold, and salt stresses, corresponding to the presence of cis-acting elements associated with stress response in their promoter regions. These findings not only improve our knowledge on lineage-specific family evolution in Brassicales, but also provide valuable information for further functional analysis of LEA genes in papaya.

Genomic insights into lineage-specific evolution of the oleosin family in Euphorbiaceae

BACKGROUND

Lipid droplets (LDs) present in land plants serve as an essential energy and carbon reserve for seed germination and seedling development. Oleosins, the most abundant structural proteins of LDs, comprise a small family involved in LD formation, stabilization and degradation. Despite their importance, our knowledge on oleosins is still poor in Euphorbiaceae, a large plant family that contains several important oil-bearing species.

RESULTS

To uncover lineage-specific evolution of oleosin genes in Euphorbiaceae, in this study, we performed a genome-wide identification and comprehensive comparison of the oleosin family in Euphorbiaceae species with available genome sequences, i.e. castor bean (Ricinus communis), physic nut (Jatropha curcas), tung tree (Vernicia fordii), Mercurialis annua, cassava (Manihot esculenta) and rubber tree (Hevea brasiliensis), and a number of five, five, five, five, eight and eight members were found, respectively. Synteny analysis revealed one-to-one collinear relationship of oleosin genes between the former four (i.e. castor bean, physic nut, tung tree and M. annua) as well as latter two species (i.e. cassava and rubber tree), whereas one-to-one and one-to-two collinear relationships were observed between physic nut and cassava, reflecting the occurrence of one recent whole-genome duplication (WGD) in the last common ancestor of cassava and rubber tree. The presence of five ortholog groups representing three previously defined clades (i.e. U, SL and SH) dates back at least to the Malpighiales ancestor, because they are also conserved in poplar (Populus trichocarpa), a tree having experienced one Salicaceae-specific recent WGD. As observed in poplar, WGD was shown to be the main driver for the family expansion in both cassava and rubber tree. Nevertheless, same retention patterns of WGD-derived duplicates observed in cassava and rubber tree are somewhat different from that of poplar, though certain homologous fragments are still present in rubber tree. Further transcriptional profiling revealed an apparent seed-predominant expression pattern of oleosin genes in physic nut, castor bean and rubber tree. Moreover, structure and expression divergence of paralogous pairs were also observed in both cassava and rubber tree.

CONCLUSION

Comparative genomics analysis of oleosin genes reported in this study improved our knowledge on lineage-specific family evolution in Euphorbiaceae, which also provides valuable information for further functional analysis and utilization of key members and their promoters.

Transcriptome analysis of Pará rubber tree (H. brasiliensis) seedlings under ethylene stimulation

BACKGROUND

Natural rubber (cis-1,4-polyioprene, NR) is an indispensable industrial raw material obtained from the Pará rubber tree (H. brasiliensis). Natural rubber cannot be replaced by synthetic rubber compounds because of the superior resilience, elasticity, abrasion resistance, efficient heat dispersion, and impact resistance of NR. In NR production, latex is harvested by periodical tapping of the trunk bark. Ethylene enhances and prolongs latex flow and latex regeneration. Ethephon, which is an ethylene-releasing compound, applied to the trunk before tapping usually results in a 1.5- to 2-fold increase in latex yield. However, intense mechanical damage to bark tissues by excessive tapping and/or over-stimulation with ethephon induces severe oxidative stress in laticifer cells, which often causes tapping panel dryness (TPD) syndrome. To enhance NR production without causing TPD, an improved understanding of the molecular mechanism of the ethylene response in the Pará rubber tree is required. Therefore, we investigated gene expression in response to ethephon treatment using Pará rubber tree seedlings as a model system.

RESULTS

After ethephon treatment, 3270 genes showed significant differences in expression compared with the mock treatment. Genes associated with carotenoids, flavonoids, and abscisic acid biosynthesis were significantly upregulated by ethephon treatment, which might contribute to an increase in latex flow. Genes associated with secondary cell wall formation were downregulated, which might be because of the reduced sugar supply. Given that sucrose is an important molecule for NR production, a trade-off may arise between NR production and cell wall formation for plant growth and for wound healing at the tapping panel.

CONCLUSIONS

Dynamic changes in gene expression occur specifically in response to ethephon treatment. Certain genes identified may potentially contribute to latex production or TPD suppression. These data provide valuable information to understand the mechanism of ethylene stimulation, and will contribute to improved management practices and/or molecular breeding to attain higher yields of latex from Pará rubber trees.

Genome-Wide Identification, Structure Characterization, and Expression Pattern Profiling of the Aquaporin Gene Family in Betula pendula

Aquaporin water channels (AQPs) constitute a large family of transmembrane proteins present throughout all kingdoms of life. They play key roles in the flux of water and many solutes across the membranes. The AQP diversity, protein features, and biological functions of silver birch are still unknown. A genome analysis of Betula pendula identified 33 putative genes encoding full-length AQP sequences (BpeAQPs). They are grouped into five subfamilies, representing ten plasma membrane intrinsic proteins (PIPs), eight tonoplast intrinsic proteins (TIPs), eight NOD26-like intrinsic proteins (NIPs), four X intrinsic proteins (XIPs), and three small basic intrinsic proteins (SIPs). The BpeAQP gene structure is conserved within each subfamily, with exon numbers ranging from one to five. The predictions of the aromatic/arginine selectivity filter (ar/R), Froger's positions, specificity-determining positions, and 2D and 3D biochemical properties indicate noticeable transport specificities to various non-aqueous substrates between members and/or subfamilies. Nevertheless, overall, the BpePIPs display mostly hydrophilic ar/R selective filter and lining-pore residues, whereas the BpeTIP, BpeNIP, BpeSIP, and BpeXIP subfamilies mostly contain hydrophobic permeation signatures. Transcriptional expression analyses indicate that 23 BpeAQP genes are transcribed, including five organ-related expressions. Surprisingly, no significant transcriptional expression is monitored in leaves in response to cold stress (6 °C), although interesting trends can be distinguished and will be discussed, notably in relation to the plasticity of this pioneer species, B. pendula. The current study presents the first detailed genome-wide analysis of the AQP gene family in a Betulaceae species, and our results lay a foundation for a better understanding of the specific functions of the BpeAQP genes in the responses of the silver birch trees to cold stress.

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.

A predominant isoform of fructokinase, HbFRK2, is involved in Hevea brasiliensis (para rubber tree) latex yield and regeneration

Fructokinase (FRK) mediates fructose phosphorylation to regulate the carbon flow and its assignment to sink tissues. Out of five HbFRKs in the genome of the rubber tree, three (HbFRK1-3) that were highly expressed in latex (cytoplasm of laticifers) were isolated and examined. According to phylogenetic analysis and intracellular location experiment, both HbFRK2 and HbFRK3 were highly possible to be expressed in cytosol, while HbFRK1 was in plastid. As the predominant isoform in laticifers, HbFRK2 had the highest transcripts, followed by HbFRK3 and HbFRK1. In enzymatic function, HbFRK2 also showed the highest affinity for fructose. To examine the roles of FRKs in latex yield and regeneration, changes in HbFRKs were examined when latex outflow from the trees were increased through two experimental interventions. In the first approach, tapping was initiated on previously untapped trees, resulting in latex yield increasing with consecutive tapping at the initial stage before it stabilized. In the second approach, latex yield from trees that were already in regular tapping was stimulated by treatment with the ethylene-based yield stimulant, ethephon. Using either method to induce an increase in latex yield, the abundance of HbFRK2 and HbFRK3 in transcripts, was increased. This development, which was especially marked in HbFRK2, may reflect a strengthening of glycolysis to meet the carbon flux and energy demands for increased rubber biosynthesis to replace rubber lost in the increased latex yield. Our results, therefore, suggest that HbFRK2 plays a critical role in fructose catabolism to facilitate rubber regeneration in the commercially exploited rubber tree.

Genes encoding light-harvesting chlorophyll a/b-binding proteins in papaya (Carica papaya L.) and insight into lineage-specific evolution in Brassicaceae

Light-harvesting chlorophyll a/b-binding (Lhc) proteins comprise a plant-specific superfamily involved in photosynthesis and stress responses. Despite their importance, little is known in papaya (Carica papaya), an economically important tree fruit crop as well as a species close to the model plant arabidopsis (Arabidopsis thaliana). This study reports a first genome-wide analysis of Lhc superfamily genes in papaya, and a total of 28 members that represent four defined families or 26 orthologous groups were identified from the papaya genome. The superfamily number is comparable to 28 or 27 reported in castor (Ricinus communis) and jatropha (Jatropha curcas), respectively, two Euphorbiaceous plants also without any recent whole-genome duplication (WGD), but relatively less than 35, 34, 32, 32, 37, 30 or 32 present in cassava (Manihot esculenta), arabidopsis, A. lyrata, A. halleri, Capsella rubella, C. grandiflora, and Eutrema salsugineum, respectively, representative species having experienced one or two recent WGDs. Local duplication was shown to play a predominant role in gene expansion in papaya, castor, and jatropha, which is only confined to the Lhcb1 group. By contrast, WGD plays a relatively more important role in cassava, arabidopsis, and other Brassicaceous plants. Further comparison of Brassicaceous plants revealed that loss of the SEP6 group in arabidopsis is lineage-specific, occurring sometime after papaya-arabidopsis divergence but before the radiation of Brassicaceous plants. Transcriptional profiling revealed a leaf-preferential expression pattern of most CpLhc superfamily genes and their transcript levels were markedly regulated by three abiotic stresses, i.e., mimicking drought, cold, and high salt. These findings not only facilitate further functional studies in papaya, but also improve our knowledge on lineage-specific evolution of this special gene superfamily in Brassicaceae.

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.

Gene-Wide Analysis of Aquaporin Gene Family in Malus domestica and Heterologous Expression of the Gene MpPIP2;1 Confers Drought and Salinity Tolerance in Arabidposis thaliana

The aquaporins (AQPs) are a family of integral membrane proteins involved in the transcellular membrane transport of water and other small molecules. A scan of the apple (Malus domestica) genome revealed the presence of 42 genes encoding putative AQPs. Based on a phylogenetic analysis of the deduced peptide sequences of the AQPs generated by Arabidopsis thaliana, poplar (Populus trichocarpa), and rubber (Hevea brasiliensis), the apple AQPs were each assigned membership of the five established AQP subfamilies, namely the PIPs (eleven members), the TIPs (thirteen members), the NIPs (eleven members), the SIPs (five members), and the XIPs (two members). The apple AQPs included asparagine-proline-alanine (NPA) motifs, an aromatic/arginine (ar/R) selectivity filter, and the Froger’s positions. The heterologous expression of MpPIP2;1 in A. thaliana was shown to enhance the level of tolerance exhibited against both drought and salinity.

Genomics analysis of genes encoding respiratory burst oxidase homologs (RBOHs) in jatropha and the comparison with castor bean

Respiratory burst oxidase homologs (RBOHs), which catalyze the production of superoxide from oxygen and NADPH, play key roles in plant growth and development, hormone signaling, and stress responses. Compared with extensive studies in model plants arabidopsis and rice, little is known about RBOHs in other species. This study presents a genome-wide analysis of Rboh family genes in jatropha (Jatropha curcas) as well as the comparison with castor bean (Ricinus communis), another economically important non-food oilseed crop of the Euphorbiaceae family. The family number of seven members identified from the jatropha genome is equal to that present in castor bean, and further phylogenetic analysis assigned these genes into seven groups named RBOHD, -C, -B, -E, -F, -N, and -H. In contrast to a high number of paralogs present in arabidopsis and rice that experienced several rounds of recent whole-genome duplications, no duplicate was identified in both jatropha and castor bean. Conserved synteny and one-to-one orthologous relationship were observed between jatropha and castor bean Rboh genes. Although exon-intron structures are usually highly conserved between orthologs, loss of certain introns was observed for JcRbohB, JcRbohD, and RcRbohN, supporting their divergence. Global gene expression profiling revealed diverse patterns of JcRbohs over various tissues. Moreover, expression patterns of JcRbohs during flower development as well as various stresses were also investigated. These findings will not only improve our knowledge on species-specific evolution of the Rboh gene family, but also provide valuable information for further functional analysis of Rboh genes in jatropha.

Genomic technologies for Hevea breeding

The commercial production of high quality natural rubber (NR) solely depends on Hevea brasiliensis Muell. Arg, (Para rubber tree) and accounts for >98% of total production worldwide. NR with its unique properties is an essential commodity for the automobile industry and its synthetic counterparts are in no way substitute to it. The rubber tree genome is very complex and plays an important role in delivering the unique properties of Hevea. But a lack of knowledge on the molecular mechanisms of rubber biosynthesis, disease resistance, etc., in elite clones of rubber still persists. Marker-assisted selection and transgenic techniques were proved to be advantageous in improving the breeding efficiency for latex yield, disease resistance, etc. The suppression subtractive hybridization (SSH), in the form of subtracted cDNA libraries and microarrays, can assist in searching the functions of expressed genes (candidate gene approach). Expressed sequence tags (ESTs) related to various metabolic aspects are well utilized to create EST banks that broadly represent the genes expressed in one tissue, such as latex cells, that assists in the study of gene function and regulation. Transcriptome analysis and gene mapping have been accomplished in Hevea at various stages. However, a selection criterion to delineate high yielding genotypes at the juvenile stage has not been accomplished so far. This is the main pit fall for rubber breeding apart from stock-scion interactions leading to yield differences among a clonally multiplied population. At least four draft genome sequences have been published on Hevea rubber, and all give different genome size and contig lengths-a comprehensive and acceptable genomic map remains unfulfilled. The progress made in molecular markers, latex biosynthesis genes, transcriptome analysis, chloroplast and mitochondrial DNA diversity, paternity identification through Breeding without Breeding (BwB), stimulated latex production and its molecular intricacies, molecular biology of tapping panel dryness, genomics for changed climates and genome mapping are discussed in this review. These information can be utilized to improvise the molecular breeding programs of Hevea in future.

Genome-wide comparison reveals divergence of cassava and rubber aquaporin family genes after the recent whole-genome duplication

Background

Aquaporins (AQPs) are a class of integral membrane proteins that facilitate the passive transport of water and other small solutes across biological membranes. Despite their importance, little information is available in cassava (Manihot esculenta), a perennial shrub of the Euphorbiaceae family that serves the sixth major staple crop in the world.

Results

This study presents a genome-wide analysis of the AQP gene family in cassava. The family of 42 members in this species could be divided into five subfamilies based on phylogenetic analysis, i.e., 14 plasma membrane intrinsic proteins (PIPs), 13 tonoplast intrinsic proteins (TIPs), nine NOD26-like intrinsic proteins (NIPs), four X intrinsic proteins (XIPs), and two small basic intrinsic proteins (SIPs). Best-reciprocal-hit-based sequence comparison and synteny analysis revealed 34 orthologous groups (OGs) present in the Euphorbiaceae ancestor, and nearly one-to-one or two-to-one orthologous relationships were observed between cassava with rubber/physic nut, reflecting the occurrence of one so-called ρ recent whole-genome duplication (WGD) in the last common ancestor of cassava and rubber. In contrast to a predominant role of the ρ WGD on family expansion in rubber, cassava AQP duplicates were derived from the WGD as well as local duplication. Species-specific gene loss was also observed in cassava, which includes the entire NIP4 group and/or six OGs. Comparison of conserved motifs and gene expression profiles revealed divergence of paralogs in cassava as observed in rubber.

Conclusions

Our findings will not only improve our knowledge on family evolution in Euphorbiaceae, but also provide valuable information for further functional analysis of AQP genes in cassava and rubber.

Electronic supplementary material

The online version of this article (10.1186/s12864-019-5780-4) contains supplementary material, which is available to authorized users.

Genome-wide identification and comparative evolutionary analysis of the Dof transcription factor family in physic nut and castor bean

DNA-binding with one finger (Dof) proteins comprise a plant-specific transcription factor family involved in plant growth, development and stress responses. This study presents a genome-wide comparison of Dof family genes in physic nut (Jatropha curcas) and castor bean (Ricinus communis), two Euphorbiaceae plants that have not experienced any recent whole-genome duplication. A total of 25 or 24 Dof genes were identified from physic nut and castor genomes, respectively, where JcDof genes are distributed across nine out of 11 chromosomes. Phylogenetic analysis assigned these genes into nine groups representing four subfamilies, and 24 orthologous groups were also proposed based on comparison of physic nut, castor, Arabidopsis and rice Dofs. Conserved microsynteny was observed between physic nut and castor Dof-coding scaffolds, which allowed anchoring of 23 RcDof genes to nine physic nut chromosomes. In contrast to how no recent duplicate was present in castor, two tandem duplications and one gene loss were found in the Dof gene family of physic nut. Global transcriptome profiling revealed diverse patterns of Jc/RcDof genes over various tissues, and key Dof genes involved in flower development and stress response were also identified in physic nut. These findings provide valuable information for further studies of Dof genes in physic nut and castor.

Genome-wide identification and characterization of the Dof gene family in cassava (Manihot esculenta)

Dof proteins comprise a plant-specific transcription factor family involved in regulation of many plant processes. Despite their diverse roles, little information is available in cassava (Manihot esculenta, Euphorbiaceae), the sixth major staple crop in the world. This study presents the first genome-wide identification of Dof family genes in cassava, resulting in 45 family members that are unevenly distributed across 17 out of the 18 chromosomes. Phylogenetic analysis assigned these genes into eight groups representing four subfamilies. Synteny analysis and BRH (Best Reciprocal Hit) method were also adopted to study evolutionary relationships. Among 21 MeDof duplicates identified, the majority of them were derived from the recent whole-genome duplication (WGD), whereas only one from local duplication. Transcriptional profiling revealed diverse expression patterns among various tissues, and further comparison of conserved motifs and expression profiles revealed divergence of MeDof duplicates. These findings provide valuable information for functional analysis of Dof genes in cassava.

Genome-Wide Identification and Characterization of Aquaporins and Their Role in the Flower Opening Processes in Carnation (Dianthus caryophyllus)

Aquaporins (AQPs) are associated with the transport of water and other small solutes across biological membranes. Genome-wide identification and characterization will pave the way for further insights into the AQPs’ roles in the commercial carnation (Dianthus caryophyllus). This study focuses on the analysis of AQPs in carnation (DcaAQPs) involved in flower opening processes. Thirty DcaAQPs were identified and grouped to five subfamilies: nine PIPs, 11 TIPs, six NIPs, three SIPs, and one XIP. Subsequently, gene structure, protein motifs, and co-expression network of DcaAQPs were analyzed and substrate specificity of DcaAQPs was predicted. qRT-PCR, RNA-seq, and semi-qRTRCR were used for DcaAQP genes expression analysis. The analysis results indicated that DcaAQPs were relatively conserved in gene structure and protein motifs, that DcaAQPs had significant differences in substrate specificity among different subfamilies, and that DcaAQP genes’ expressions were significantly different in roots, stems, leaves and flowers. Five DcaAQP genes (DcaPIP1;3, DcaPIP2;2, DcaPIP2;5, DcaTIP1;4, and DcaTIP2;2) might play important roles in flower opening process. However, the roles they play are different in flower organs, namely, sepals, petals, stamens, and pistils. Overall, this study provides a theoretical basis for further functional analysis of DcaAQPs.

Genome-wide comparative analysis of papain-like cysteine protease family genes in castor bean and physic nut

Papain-like cysteine proteases (PLCPs) are a class of proteolytic enzymes involved in many plant processes. Compared with the extensive research in Arabidopsis thaliana, little is known in castor bean (Ricinus communis) and physic nut (Jatropha curcas), two Euphorbiaceous plants without any recent whole-genome duplication. In this study, a total of 26 or 23 PLCP genes were identified from the genomes of castor bean and physic nut respectively, which can be divided into nine subfamilies based on the phylogenetic analysis: RD21, CEP, XCP, XBCP3, THI, SAG12, RD19, ALP and CTB. Although most of them harbor orthologs in Arabidopsis, several members in subfamilies RD21, CEP, XBCP3 and SAG12 form new groups or subgroups as observed in other species, suggesting specific gene loss occurred in Arabidopsis. Recent gene duplicates were also identified in these two species, but they are limited to the SAG12 subfamily and were all derived from local duplication. Expression profiling revealed diverse patterns of different family members over various tissues. Furthermore, the evolution characteristics of PLCP genes were also compared and discussed. Our findings provide a useful reference to characterize PLCP genes and investigate the family evolution in Euphorbiaceae and species beyond.

Transcript Profiling of Hevea brasiliensis during Latex Flow

Latex exploitation enhances latex regeneration in rubber trees. The latex exploitation-caused latex flow lasts from 10 min to a few hours, which is convenient for exploring the transcript profiling of latex metabolism-related genes at the different stages of latex flow. In the present study, the expression pattern of 62 latex metabolism-related genes involved in water transportation, carbohydrate metabolism, natural rubber biosynthesis, hormone signaling, ROS generation and scavenging, and latex coagulum across three stages of latex flow between rubber tree clones CATAS7-33-97 and CATAS8-79 were comparatively analyzed by quantitative real-time PCR. The two clones show differences in latex regeneration and have a different duration of latex flow. The results showed that the expression levels of 38 genes were significantly higher in CATAS8-79 latex than in CATAS7-33-97 during latex regeneration, while 45 genes had a notably higher expression level in CATAS8-79 latex during latex flow. Together with the activation of the MEP pathway and jasmonate pathway in CATAS8-79 latex, HbPIP1;3, HbPIP1;4, HbSUT3, HbSus3, HbHMGS1-2, HbMK should contribute to the high latex regeneration ability. The up-regulation of ethylene signaling and Hb44KD and the down-regulation of latex coagulation-related genes in CATAS8-79 latex might contribute to its longer latex flow duration. This study provides some cues for revealing the regulation of latex metabolism in rubber trees.

Papain-like cysteine protease encoding genes in rubber (Hevea brasiliensis): comparative genomics, phylogenetic, and transcriptional profiling analysis

43 HbPLCPs representing nine subfamilies or 20 orthologous groups were found in rubber, where paralogs were resulted from the recent WGD and local duplication. Several senescence-associated genes were also identified. Papain-like cysteine proteases (PLCPs) comprise a large family of proteolytic enzymes involved in plant growth and development, seed germination, organ senescence, immunity, and stress response. Despite their importance and the extensive research in the model plant Arabidopsis thaliana, little information is available on rubber tree (Hevea brasiliensis), a rubber-producing plant of the Euphorbiaceae family. This study performed a genome-wide identification of PLCP family genes in rubber, resulting in a relatively high number of 43 members. The phylogenetic analysis assigned these genes into nine subfamilies, i.e., RD21 (6), CEP (4), XCP (4), XBCP3 (2), THI (1), SAG12 (18), RD19 (4), ALP (2), and CTB (2). Most of them were shown to have orthologs in Arabidopsis; however, several members in SAG12, CEP and XBCP3 subfamilies form new groups as observed in other core eudicots such as Manihot esculenta, Ricinus communis, Populus trichocarpa, and Vitis vinifera. Based on an expert sequence comparison, 20 orthologous groups (OGs) were proposed for core eudicots, and rubber paralogs were shown to be resulted from the recent whole-genome duplication (WGD) as well as local duplication. Transcriptional profiling showed distinct expression pattern of different members across various tissues, e.g., root, leaf, bark, laticifer, flower, and seed. By using the senescence-specific HbSAG12H1 as the indicator, the transcriptome of senescent rubber leaves was deeply sequenced and several senescence-associated PLCP genes were identified. Results obtained from this study provide valuable information for future functional analysis and utilization of PLCP genes in Hevea and other species.

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.

Survey of the rubber tree genome reveals a high number of cysteine protease-encoding genes homologous to Arabidopsis SAG12

Arabidopsis thaliana SAG12, a senescence-specific gene encoding a cysteine protease, is widely used as a molecular marker for the study of leaf senescence. To date, its potential orthologues have been isolated from several plant species such as Brassica napus and Nicotiana tabacum. However, little information is available in rubber tree (Hevea brasiliensis), a rubber-producing plant of the Euphorbiaceae family. This study presents the identification of SAG12-like genes from the rubber tree genome. Results showed that an unexpected high number of 17 rubber orthologues with a single intron were found, contrasting the single copy with two introns in Arabidopsis. The gene expansion was also observed in another two Euphorbiaceae plants, castor bean (Ricinus communis) and physic nut (Jatropha curcas), both of which contain 8 orthologues. In accordance with no occurrence of recent whole-genome duplication (WGD) events, most duplicates in castor and physic nut were resulted from tandem duplications. In contrast, the duplicated HbSAG12H genes were derived from tandem duplications as well as the recent WGD. Expression analysis showed that most HbSAG12H genes were lowly expressed in examined tissues except for root and male flower. Furthermore, HbSAG12H1 exhibits a strictly senescence-associated expression pattern in rubber tree leaves, and thus can be used as a marker gene for the study of senescence mechanism in Hevea.

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.

The Eucalyptus Tonoplast Intrinsic Protein (TIP) Gene Subfamily: Genomic Organization, Structural Features, and Expression Profiles

Plant aquaporins are water channels implicated in various physiological processes, including growth, development and adaptation to stress. In this study, the Tonoplast Intrinsic Protein (TIP) gene subfamily of Eucalyptus, an economically important woody species, was investigated and characterized. A genome-wide survey of the Eucalyptus grandis genome revealed the presence of eleven putative TIP genes (referred as EgTIP), which were individually assigned by phylogeny to each of the classical TIP1-5 groups. Homology modeling confirmed the presence of the two highly conserved NPA (Asn-Pro-Ala) motifs in the identified EgTIPs. Residue variations in the corresponding selectivity filters, that might reflect differences in EgTIP substrate specificity, were observed. All EgTIP genes, except EgTIP5.1, were transcribed and the majority of them showed organ/tissue-enriched expression. Inspection of the EgTIP promoters revealed the presence of common cis-regulatory elements implicated in abiotic stress and hormone responses pointing to an involvement of the identified genes in abiotic stress responses. In line with these observations, additional gene expression profiling demonstrated increased expression under polyethylene glycol-imposed osmotic stress. Overall, the results obtained suggest that these novel EgTIPs might be functionally implicated in eucalyptus adaptation to stress.

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.

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.