Subcellular Localization Properties Research Articles

Human CKAP2L shows a cell cycle-dependent expression pattern and exhibits microtubule-stabilizing properties.

Cytoskeleton-associated protein 2-like (CKAP2L) is a paralogue of cytoskeleton-associated protein 2 (CKAP2). We characterized the expression pattern, subcellular localization, and microtubule-stabilizing properties of human CKAP2L. The levels of both CKAP2L transcript and protein were cell cycle phase-dependent, peaking during the G2/M phase and relatively high in certain human tissues, including testis, intestine, and spleen. CKAP2L protein was detectable in all human cancer cell lines we tested. CKAP2L localized to the mitotic spindle apparatus during mitosis, as reported previously. During interphase, however, CKAP2L localized mainly to the nucleus. Ectopic overexpression of CKAP2L resulted in 'microtubule bundling', and, consequently, an elevated CKAP2L level led to prolonged mitosis. These findings support the mitotic role of CKAP2L during the human cell cycle.

Stochastic RNA editing of the Complexin C-terminus within single neurons regulates neurotransmitter release

Neurotransmitter release requires assembly of the SNARE complex fusion machinery, with multiple SNARE-binding proteins regulating when and where synaptic vesicle fusion occurs. The presynaptic protein Complexin (Cpx) controls spontaneous and evoked neurotransmitter release by modulating SNARE complex zippering. Although the central SNARE-binding helix is essential, post-translational modifications to Cpx's C-terminal membrane-binding amphipathic helix regulate its ability to control synaptic vesicle fusion. Here, we demonstrate that RNA editing of the Cpx C-terminus modifies its ability to clamp SNARE-mediated fusion and alters presynaptic output. RNA editing of Cpx across single neurons is stochastic, generating up to eight edit variants that fine tune neurotransmitter release by altering the subcellular localization and clamping properties of the protein. Similar stochastic editing rules for other synaptic genes were observed, indicating editing variability at single adenosines and across multiple mRNAs generates unique synaptic proteomes within the same population of neurons to fine tune presynaptic output.

Purification and biochemical analysis of native AMPA receptors from three different mammalian species.

The majority of fast, excitatory synaptic transmission in the central nervous system (CNS) is mediated by α-amino-3-hydroxy-5-methyl-4-isoxazolepropionic acid receptors (AMPARs), which are glutamate-activated ion channels integral to synaptic plasticity, motor coordination, learning, and memory. Native AMPARs are multiprotein assemblies comprised of a tetrameric receptor core that co-assembles with a broad range of peripheral auxiliary proteins which shape subcellular localization and signaling properties of the resulting complexes. Structure determination of AMPARs has traditionally relied on recombinant expression systems; however, these methods are not well suited to elucidate the diverse array of AMPAR assemblies that are differentially expressed in mammalian brains. While recent studies of native receptor complexes have advanced our understanding of endogenous assemblies, receptors thus far have only been isolated from rodent brain tissue. Here, we employed an immunoaffinity purification strategy to isolate native AMPARs from the brains of three different mammals-pigs, sheep, and cows. Compared to rodents, pigs, sheep, and cows are ungulate mammals, animals with closer genomic identity with humans. Here we determined the molecular size, overall yield, and purity of native AMPARs isolated from these three mammals, thereby demonstrating that structural determination and biochemical analysis is possible from a clade of mammals evolutionarily distinct from rodents.

Regulation of Neurotransmitter Release byK+ Channels.

K+ channels play potent roles in the process of neurotransmitter release by influencing the action potential waveform and modulating neuronal excitability and release probability. These diverse effects of K+ channel activation are ensured by the wide variety of K+ channel genes and their differential expression in different cell types. Accordingly, a variety of K+ channels have been implicated in regulating neurotransmitter release, including the Ca2+- and voltage-gated K+ channel Slo1 (also known as BK channel), voltage-gated K+ channels of the Kv3 (Shaw-type), Kv1 (Shaker-type), and Kv7 (KCNQ) families, G-protein-gated inwardly rectifying K+ (GIRK) channels, and SLO-2 (a Ca2+-. Cl-, and voltage-gated K+ channel in C. elegans). These channels vary in their expression patterns, subcellular localization, and biophysical properties. Their roles in neurotransmitter release may also vary depending on the synapse and physiological or experimental conditions. This chapter summarizes key findings about the roles of K+ channels in regulating neurotransmitter release.

Genome-Wide Identification and Analysis of the R2R3-MYB Gene Family in Theobroma cacao.

The MYB gene family is involved in the regulation of plant growth, development and stress responses. In this paper, to identify Theobroma cacao R2R3-MYB (TcMYB) genes involved in environmental stress and phytohormones, we conducted a genome-wide analysis of the R2R3-MYB gene family in Theobroma cacao (cacao). A total of 116 TcMYB genes were identified, and they were divided into 23 subgroups according to the phylogenetic analysis. Meanwhile, the conserved motifs, gene structures and cis-acting elements of promoters were analyzed. Moreover, these TcMYB genes were distributed on 10 chromosomes. We conducted a synteny analysis to understand the evolution of the cacao R2R3-MYB gene family. A total of 37 gene pairs of TcMYB genes were identified through tandem or segmental duplication events. Additionally, we also predicted the subcellular localization and physicochemical properties. All the studies showed that TcMYB genes have multiple functions, including responding to environmental stresses. The results provide an understanding of R2R3-MYB in Theobroma cacao and lay the foundation for a further functional analysis of TcMYB genes in the growth of cacao.

Genome-wide identification and expression analyses of late embryogenesis abundant (LEA) gene family in tobacco (Nicotiana tabacum L.) reveal their function in abiotic stress responses

Late embryogenesis abundant (LEA) proteins play an important role in plant growth and response to abiotic stresses. However the late embryogenesis abundant (LEA) gene family in Nicotiana tabacum has not been systematically studied. In this study, 123 NtLEA genes were identified in Nicotiana tabacum, and divided into 8 groups, including LEA_1, LEA_2, LEA_3, LEA_4, LEA_5, LEA_6, DHN (dehydratin) and SMP (Seed Maturation Protein). The LEA_2 group is the most abundant of the NtLEA family. The gene structure, conserved motifs, subcellular localization and physicochemical properties of the NtLEA genes were analyzed. RNA-seq and qPCR analyses showed that the NtLEA genes were significantly induced under two different abiotic stresses and showed different expression patterns. The expression patterns of 35 NtLEA genes responding to ABA and 3 NtLEA genes responding to NaCl abiotic stress, respectively, were characterized. The protein–protein interaction network revealed that most NtLEA proteins (>78%) had the potential function to enhance tobacco resistance to abiotic stress. The transcriptional regulatory network showed that 21 transcription factor families were involved in regulating the expression of the NtLEA genes. These results are beneficial for future studies of the function of the NtLEA genes.

Sialidase NEU3 and its pathological significance.

Sialidases (EC 3.2.1.18, also called neuraminidases) catalyze the removal of α-glycosidically linked sialic acid residues from glycoproteins and glycolipids; this is the initial step in the degradation of these glycoconjugates. Sialidases of mammalian origin have been implicated in not only lysosomal catabolism but also the modulation of functional molecules involved in many biological processes. To date, four types of mammalian sialidases have been cloned and designated as Neu1, Neu2, Neu3 and Neu4. These sialidases differ in their subcellular localization and enzymatic properties, as well as their chromosomal localization, and they are expressed in a tissue-specific manner. Among the sialidases, the plasma membrane-associated sialidase Neu3 appears to play particular roles in controlling transmembrane signaling through the modulation of gangliosides, and its aberrant expression is closely related to various pathogeneses, including that of cancer. Interestingly, the human orthologue NEU3 acts in two ways, catalytic hydrolysis of gangliosides and protein interactions with other signaling molecules. Aberrant NEU3 expression can induce various pathological conditions. This review briefly summarizes recent studies, focusing on the involvement of NEU3 in various pathological phenomena.

In Silico Characterization of a Hypothetical Protein from Shigella dysenteriae ATCC 12039 Reveals a Pathogenesis-Related Protein of the Type-VI Secretion System.

Shigellosis caused by Shigella dysenteriae is a major public health concern worldwide, particularly in developing countries. The bacterial genome is known, but there are many hypothetical proteins whose functions are yet to be discovered. A hypothetical protein (accession no. WP_128879999.1, 161 residues) of S. dysenteriae ATCC 12039 strain was selected in this study for comprehensive structural and functional analysis. Subcellular localization and different physicochemical properties of this hypothetical protein were estimated indicating it as a stable, soluble, and extracellular protein. Functional annotation tools, such as NCBI-CD Search, Pfam, and InterProScan, predicted our target protein to be an amidase effector protein 4 (Tae4) of type-VI secretion system (T6SS). Multiple sequence alignment of the homologous sequences coincided with previous findings. Random coil was found to be predominant in secondary structure. Three-dimensional (3D) structure of the protein was obtained using homology modeling method by SWISS-MODEL server using a template protein (PDB ID: 4J30) of 80.12% sequence identity. The 3D structure became more stable after YASARA energy minimization and was validated by several quality assessment tools like PROCHECK, QMEAN, Verify3D, and ERRAT. Superimposition of the target with the template protein by UCSF Chimera generated RMSD value of 0.115 Å, suggesting a reliable 3D structure. The active site of the modeled structure was predicted and visualized by CASTp server and PyMOL. Interestingly, similar binding affinity and key interacting residues were found for the target protein and a Salmonella enterica Tae4 protein with the ligand L-Ala D-Glu-mDAP by molecular docking analysis. Protein-protein docking was also performed between the target protein and hemolysin coregulated protein 1 of T6SS. Finally, the protein was found to be a unique protein of S. dysenteriae nonhomologous to human by comparative genomics approach indicating a potential therapeutic target. Most pathogens harboring T6SS in their system pose a significant threat to the human health. Many T6SSs and their effectors are associated with interbacterial competition, pathogenesis, and virulency; however, relationships between these effectors and pathogenicity of S. dysenteriae are yet to be determined. The study findings provide a lucrative platform for future antibacterial treatment.

Proteomic Approaches to Investigate Regulated Trafficking and Signaling of G Protein-Coupled Receptors.

Advances in proteomic methodologies based on quantitative mass spectrometry are now transforming pharmacology and experimental biology more broadly. The present review will discuss several examples based on work in the author's laboratory, which focuses on delineating relationships between G protein-coupled receptor signaling and trafficking in the endocytic network. The examples highlighted correspond to those discussed in a talk presented at the 2019 EB/ASPET meeting, which was organized by Professor Joe Beavo to commemorate his receipt of the Julius Axelrod Award. SIGNIFICANCE STATEMENT: GPCRs are allosteric machines that signal by interacting with other cellular proteins, and this, in turn, is determined by a complex interplay between the biochemical, subcellular localization, and membrane trafficking properties of receptors relative to transducer and regulatory proteins. The present minireview highlights recent advances and challenges in elucidating this dynamic cell biology and toward delineating the cellular basis of drug action at the level of defined GPCR interaction networks using proteomic approaches enabled by quantitative mass spectrometry.

Nucleocytoplasmic shuttling of the glucocorticoid receptor is influenced by tetratricopeptide repeat-containing proteins.

It has been demonstrated that tetratricopeptide-repeat (TPR) domain proteins regulate the subcellular localization of glucocorticoid receptor (GR). This study analyses the influence of the TPR domain of high molecular weight immunophilins in the retrograde transport and nuclear retention of GR. Overexpression of the TPR peptide prevented efficient nuclear accumulation of the GR by disrupting the formation of complexes with the dynein-associated immunophilin FKBP52 (also known as FKBP4), the adaptor transporter importin-β1 (KPNB1), the nuclear pore-associated glycoprotein Nup62 and nuclear matrix-associated structures. We also show that nuclear import of GR was impaired, whereas GR nuclear export was enhanced. Interestingly, the CRM1 (exportin-1) inhibitor leptomycin-B abolished the effects of TPR peptide overexpression, although the drug did not inhibit GR nuclear export itself. This indicates the existence of a TPR-domain-dependent mechanism for the export of nuclear proteins. The expression balance of those TPR domain proteins bound to the GR-Hsp90 complex may determine the subcellular localization and nucleocytoplasmic properties of the receptor, and thereby its pleiotropic biological properties in different tissues and cell types.

In silico cloning and comparative analysis of NAD1 gene in three common human parasites

To in silico clone the NAD1 gene of three common parasites and analyze their bioinformatics, so as to lay the foundation for further research on the NAD gene. By using the in silico cloning method, the full length cDNA (s) of NAD 1 genes of Clonorchis sinensis, Ascaris lumbricoides and Schistosoma japonicum were got, then their physical and chemical properties, compositions of amino acids, subcellular localizations, binary and ternary structures were contrastively analyzed. The three kinds of NAD1 proteins were similar in the relative molecular weight, subcellular localization, and physical and chemical properties. The NAD1 proteins were highly similar in binary and ternary structures of A. lumbricoides and S. japonicum. The phylogenetic analysis showed that C. sinensis, A. lumbricoides and S. japonicum belonged to the different evolutionary branches with a certain of genetic distance. The three NAD1 genes got from C. sinensis, A. lumbricoides and S. japonicum by in silico cloning belong to the same gene of different species, which can be widely used in the researches of heritable variation of parasites.

Improving prediction of heterodimeric protein complexes using combination with pairwise kernel

BackgroundSince many proteins become functional only after they interact with their partner proteins and form protein complexes, it is essential to identify the sets of proteins that form complexes. Therefore, several computational methods have been proposed to predict complexes from the topology and structure of experimental protein-protein interaction (PPI) network. These methods work well to predict complexes involving at least three proteins, but generally fail at identifying complexes involving only two different proteins, called heterodimeric complexes or heterodimers. There is however an urgent need for efficient methods to predict heterodimers, since the majority of known protein complexes are precisely heterodimers.ResultsIn this paper, we use three promising kernel functions, Min kernel and two pairwise kernels, which are Metric Learning Pairwise Kernel (MLPK) and Tensor Product Pairwise Kernel (TPPK). We also consider the normalization forms of Min kernel. Then, we combine Min kernel or its normalization form and one of the pairwise kernels by plugging. We applied kernels based on PPI, domain, phylogenetic profile, and subcellular localization properties to predicting heterodimers. Then, we evaluate our method by employing C-Support Vector Classification (C-SVC), carrying out 10-fold cross-validation, and calculating the average F-measures. The results suggest that the combination of normalized-Min-kernel and MLPK leads to the best F-measure and improved the performance of our previous work, which had been the best existing method so far.ConclusionsWe propose new methods to predict heterodimers, using a machine learning-based approach. We train a support vector machine (SVM) to discriminate interacting vs non-interacting protein pairs, based on informations extracted from PPI, domain, phylogenetic profiles and subcellular localization. We evaluate in detail new kernel functions to encode these data, and report prediction performance that outperforms the state-of-the-art.

Biological and Pathological Roles of Ganglioside Sialidases.

Sialidases are glycosidases responsible for the removal of α-glycosidically linked sialic acid residues from carbohydrate portions of glycoproteins and glycolipids, this process being the initial step in the degradation of such glycoconjugates. Sialic acids are considered to play important roles in various biological processes largely in two ways, one related to their hydrophilic and acidic properties exerting physicochemical effects on the glycoconjugates to which they are attached, and the other as recognition sites or in an opposing fashion as masking sites. The removal of sialic acids catalyzed by a sialidase, therefore greatly influences many biological processes through changing the conformation of glycoproteins and through recognition and masking of biological sites of functional molecules. Sialidases are found widely distributed in metazoan animals, from echinoderms to mammals, and are also present in viruses and other microorganisms, including fungi, protozoa, and bacteria even mostly lacking sialic acids. In mammals, there are four forms of sialidase (Neu1, Neu2, Neu3, and Neu4), differing in their major subcellular localization and enzymatic properties. They have been implicated in regulation of various cellular activities, such as cell differentiation, cell growth, and cell adhesion and motility, depending on their particular properties. In contrast, in microorganisms the enzymes appear to play roles limited to nutrition and pathogenesis. In this chapter, the focus is on mammalian sialidases preferentially hydrolyzing gangliosides, mostly Neu3 and Neu4, with an attempt to provide a brief overview of their physiological and pathological roles.

Silencing of D-Lactate Dehydrogenase Impedes Glyoxalase System and Leads to Methylglyoxal Accumulation and Growth Inhibition in Rice.

D-Lactate is oxidized by two classes of D-lactate dehydrogenase (D-LDH), namely, NAD-dependent and NAD-independent D-LDHs. Little is known about the characteristics and biological functions of D-LDHs in rice. In this study, a functional NAD-independent D-LDH (LOC_Os07g06890) was identified in rice, as a result of alternative splicing events. Characterization of the expression profile, subcellular localization, and enzymatic properties of the functional OsD-LDH revealed that it is a mitochondrial cytochrome-c-dependent D-LDH with high affinity and catalytic efficiency. Functional analysis of OsD-LDH RNAi transgenic rice demonstrated that OsD-LDH participates in methylglyoxal metabolism by affecting the activity of the glyoxalase system and aldo-keto reductases. Under methylglyoxal treatment, silencing of OsD-LDH in rice resulted in the accumulation of methylglyoxal and D-lactate, the decrease of reduced glutathione in leaves, and ultimately severe growth inhibition. Moreover, the detached leaves of OsD-LDH RNAi plants were more sensitive to salt stress. However, the silencing of OsD-LDH did not affect the growth under photorespiration conditions. Our results provide new insights into the role of NAD-independent D-LDHs in rice.

Organelle-specific single-molecule imaging of α4β2 nicotinic receptors reveals the effect of nicotine on receptor assembly and cell-surface trafficking

Nicotinic acetylcholine receptors (nAChRs) assemble in the endoplasmic reticulum (ER) and traffic to the cell surface as pentamers composed of α and β subunits. Many nAChR subtypes can assemble with varying subunit ratios, giving rise to multiple stoichiometries exhibiting different subcellular localization and functional properties. In addition to the endogenous neurotransmitter acetylcholine, nicotine also binds and activates nAChRs and influences their trafficking and expression on the cell surface. Currently, no available technique can specifically elucidate the stoichiometry of nAChRs in the ER versus those in the plasma membrane. Here, we report a method involving single-molecule fluorescence measurements to determine the structural properties of these membrane proteins after isolation in nanoscale vesicles derived from specific organelles. These cell-derived nanovesicles allowed us to separate single membrane receptors while maintaining them in their physiological environment. Sorting the vesicles according to the organelle of origin enabled us to determine localized differences in receptor structural properties, structural influence on transport between organelles, and changes in receptor assembly within intracellular organelles. These organelle-specific nanovesicles revealed that one structural isoform of the α4β2 nAChR was preferentially trafficked to the cell surface. Moreover, nicotine altered nAChR assembly in the ER, resulting in increased production of the receptor isoform that traffics more efficiently to the cell surface. We conclude that the combined effects of the increased assembly of one nAChR stoichiometry and its preferential trafficking likely drive the up-regulation of nAChRs on the cell surface upon nicotine exposure.

Localization and functional characterization of the p.Asn965Ser (N965S) ABCA4 variant in mice reveal pathogenic mechanisms underlying Stargardt macular degeneration.

ABCA4 is a member of the superfamily of ATP-binding cassette (ABC) proteins that transports N-retinylidene-phosphatidylethanolamine (N-Ret-PE) across outer segment disc membranes thereby facilitating the removal of potentially toxic retinoid compounds from photoreceptor cells. Mutations in the gene encoding ABCA4 are responsible for Stargardt disease (STGD1), an autosomal recessive retinal degenerative disease that causes severe vision loss. To define the molecular basis for STGD1 associated with the p.Asn965Ser (N965S) mutation in the Walker A motif of nucleotide binding domain 1 (NBD1), we generated a p.Asn965Ser knockin mouse and compared the subcellular localization and molecular properties of the disease variant with wild-type (WT) ABCA4. Here, we show that the p.Asn965Ser ABCA4 variant expresses at half the level of WT ABCA4, partially mislocalizes to the endoplasmic reticulum (ER) of photoreceptors, is devoid of N-Ret-PE activated ATPase activity, and causes an increase in autofluorescence and the bisretinoid A2E associated with lipofuscin deposits in retinal pigment epithelial cells as found in Stargardt patients and Abca4 knockout mice. We also show for the first time that a significant fraction of WT ABCA4 is retained in the inner segment of photoreceptors. On the basis of these studies we conclude that loss in substrate-dependent ATPase activity and protein misfolding are mechanisms underlying STGD1 associated with the p.Asn965Ser mutation in ABCA4. Functional and molecular modeling studies further suggest that similar pathogenic mechanisms are responsible for Tangiers disease associated with the p.Asn935Ser (N935S) mutation in the NBD1 Walker A motif of ABCA1.

The glycine hinge of transmembrane segment 2 modulates the subcellular localization and gating properties in TREK channels

TWIK-Related K+ channels (TREK), including TREK-1 and TREK-2, belong to the TREK/TRAAK subclass of two-pore domain K+ (K2P) family. The important functions of transmembrane segment 4 (M4)-glycine hinge in TREK channel gating have been characterized, but the roles of M2-hinge (the equivalent residue of M4-hinge) remain unclear. Here, by characterizing the macroscopic currents, subcellular localization and gating properties of their M2-hinge mutants (G166A for TREK-1 and G196A for TREK-2), we investigated the functions of M2-hinge. G166A displayed decreased whole-cell currents, whereas no current was produced by G196A. Subcellular analysis indicated that both mutants were aggregated near the perinuclear region, and most of them were retented within the endoplasmic reticulum (ER). Next, to explore the roles of M2-hinge in the gating mechanism, we tested the responses of the related M2-hinge mutants to 2-Aminoethoxydiphenyl borate (2-APB) and extracellular pH alteration (ΔpHo). TREK-1mut7-G166A displayed reduced sensitivity to 2-APB activation, but similar sensitivity to ΔpHo, when compared with TREK-1mut7. WT-ΔpCt, a TREK-2 tandom dimer, was used to assess the function of M2-hinge in the cis-type gating of TREK-2. The sensitivities of G196A-ΔpCt to both 2-APB and ΔpHo decreased compared with WT-ΔpCt. Taken together, our results reveal that the M2-hinge of TREK channels control their macroscopic current, subcellular localization and gating process.

Optimized FISH methods for visualizing RNA localization properties in Drosophila and human tissues and cultured cells.

Eukaryotic gene expression is orchestrated by a large number of regulatory steps to modulate the synthesis, maturation and fate of various families of protein-coding and non-coding RNA molecules. Defining the subcellular localization properties of an RNA molecule is thus of considerable importance for gleaning its function(s) and for elucidating post-transcriptional gene regulation pathways. For decades, fluorescent In Situ hybridization (FISH) has constituted the gold-standard technique for assessing RNA expression and distribution properties in cultured cells, tissue specimens, and whole mount organisms. Recently, several attempts aimed at advancing multiplex RNA-FISH experiments have been published. However, these procedures are both financially demanding and technically challenging, while their full potential remains unexploited. Here we describe an optimized RNA-FISH method employing the Tyramide Signal Amplification system that robustly enhances resolution and sensitivity needed for exploring RNA localization in Drosophila embryos, tissues and commonly cultured human and insect cell lines. Methodological details and key parameters are outlined for high-throughput analyses conducted in 96-well plate format.

Structure-function analysis of human stomatin: A mutation study.

Stomatin is an ancient, widely expressed, oligomeric, monotopic membrane protein that is associated with cholesterol-rich membranes/lipid rafts. It is part of the SPFH superfamily including stomatin-like proteins, prohibitins, flotillin/reggie proteins, bacterial HflK/C proteins and erlins. Biochemical features such as palmitoylation, oligomerization, and hydrophobic “hairpin” structure show similarity to caveolins and other integral scaffolding proteins. Recent structure analyses of the conserved PHB/SPFH domain revealed amino acid residues and subdomains that appear essential for the structure and function of stomatin. To test the significance of these residues and domains, we exchanged or deleted them, expressed respective GFP-tagged mutants, and studied their subcellular localization, molecular dynamics and biochemical properties. We show that stomatin is a cholesterol binding protein and that at least two domains are important for the association with cholesterol-rich membranes. The conserved, prominent coiled-coil domain is necessary for oligomerization, while association with cholesterol-rich membranes is also involved in oligomer formation. FRAP analyses indicate that the C-terminus is the dominant entity for lateral mobility and binding site for the cortical actin cytoskeleton.

The Interactions of Aquaporins and Mineral Nutrients in Higher Plants.

Aquaporins, major intrinsic proteins (MIPs) present in the plasma and intracellular membranes, facilitate the transport of small neutral molecules across cell membranes in higher plants. Recently, progress has been made in understanding the mechanisms of aquaporin subcellular localization, transport selectivity, and gating properties. Although the role of aquaporins in maintaining the plant water status has been addressed, the interactions between plant aquaporins and mineral nutrients remain largely unknown. This review highlights the roles of various aquaporin orthologues in mineral nutrient uptake and transport, as well as the regulatory effects of mineral nutrients on aquaporin expression and activity, and an integrated link between aquaporins and mineral nutrient metabolism was identified.