PKS2: A link between phototropin signalling and auxin transport - a study on how plants sense and respond to light.

The blue light photoreceptors phototropins (phot1 and phot2 in Arabidopsis thaliana (L.)) carry out various light responses of great adaptive value that optimize plant growth. These processes include phototropism (the bending of an organ induced by unequal light distribution), chloroplast movements, stomatal opening, leaf flattening and solar tracking. The biochemical pathways controlling these important blue light responses are just starting to be elucidated. The PHYTOCHROME KINASE SUBSTRATE (PKS1-4) proteins - the subject of this research - have recently been identified as novel phototropism signalling components. PKS1 (the founding member of this family) interacts in a same complex in vivo with phot1 and the important phot1 signalling element NON-PHOTOTROPIC HYPOCOTYL 3 (NPH3). This suggested that the PKS may act as early components of phot signalling. This work further investigates the role of this protein family during phototropin signalling Genetic experiments clearly showed that the PKS do not control chloroplast movements or stomatal opening. However, PKS2 plays a critical role with NPH3 during leaf flattening and solar tracking. Epistasis data indicated that both proteins act in phot1 and phot2 pathways, which is consistent with their in vivo interaction with both phototropins. Because phototropism, leaf flattening and solar tracking are developmental processes regulated by the hormone auxin, the role of PKS2 and NPH3 during auxin homeostasis was also investigated. Interestingly, PKS2 loss-of-function restores leaf flattening in the auxin transporter mutant aux1. Moreover, PKS2 and NPH3 are found in a same complex with AUX1 in vivo. Taken together, these results suggest that PKS2 may act with NPH3 as a connecting point between phot signalling and auxin transport. Further experiments were performed to explore the molecular mode of action of PKS2 and NPH3 in this process. The significance of these results is discussed.

Peroxisomal polyhydroxyalkanoate biosynthesis is a promising strategy for bioplastic production in high biomass crops.

Polyhydroxyalkanoates (PHAs) are bacterial carbon storage polymers with diverse plastic-like properties. PHA biosynthesis in transgenic plants is being developed as a way to reduce the cost and increase the sustainability of industrial PHA production. The homopolymer polyhydroxybutyrate (PHB) is the simplest form of these biodegradable polyesters. Plant peroxisomes contain the substrate molecules and necessary reducing power for PHB biosynthesis, but peroxisomal PHB production has not been explored in whole soil-grown transgenic plants to date. We generated transgenic sugarcane (Saccharum sp.) with the three-enzyme Ralstonia eutropha PHA biosynthetic pathway targeted to peroxisomes. We also introduced the pathway into Arabidopsis thaliana, as a model system for studying and manipulating peroxisomal PHB production. PHB, at levels up to 1.6%-1.8% dry weight, accumulated in sugarcane leaves and A. thaliana seedlings, respectively. In sugarcane, PHB accumulated throughout most leaf cell types in both peroxisomes and vacuoles. A small percentage of total polymer was also identified as the copolymer poly (3-hydroxybutyrate-co-3-hydroxyvalerate) in both plant species. No obvious deleterious effect was observed on plant growth because of peroxisomal PHA biosynthesis at these levels. This study highlights how using peroxisomal metabolism for PHA biosynthesis could significantly contribute to reaching commercial production levels of PHAs in crop plants.

Phytochrome interacting factors 4 and 5 control seedling growth in changing light conditions by directly controlling auxin signaling.

Plant growth is strongly influenced by the presence of neighbors that compete for light resources. In response to vegetational shading shade-intolerant plants such as Arabidopsis display a suite of developmental responses known as the shade-avoidance syndrome (SAS). The phytochrome B (phyB) photoreceptor is the major light sensor to mediate this adaptive response. Control of the SAS occurs in part with phyB, which controls protein abundance of phytochrome-interacting factors 4 and 5 (PIF4 and PIF5) directly. The shade-avoidance response also requires rapid biosynthesis of auxin and its transport to promote elongation growth. The identification of genome-wide PIF5-binding sites during shade avoidance revealed that this bHLH transcription factor regulates the expression of a subset of previously identified SAS genes. Moreover our study suggests that PIF4 and PIF5 regulate elongation growth by controlling directly the expression of genes that code for auxin biosynthesis and auxin signaling components.

Reduced peroxisomal citrate synthase activity increases substrate availability for polyhydroxyalkanoate biosynthesis in plant peroxisomes.

Polyhydroxyalkanoates (PHAs) are bacterial carbon storage polymers used as renewable, biodegradable plastics. PHA production in plants may be a way to reduce industrial PHA production costs. We recently demonstrated a promising level of peroxisomal PHA production in the high biomass crop species sugarcane. However, further production strategies are needed to boost PHA accumulation closer to commercial targets. Through exogenous fatty acid feeding of Arabidopsis thaliana plants that contain peroxisome-targeted PhaA, PhaB and PhaC enzymes from Cupriavidus necator, we show here that the availability of substrates derived from the β-oxidation cycle limits peroxisomal polyhydroxybutyrate (PHB) biosynthesis. Knockdown of peroxisomal citrate synthase activity using artificial microRNA increased PHB production levels approximately threefold. This work demonstrates that reduction of peroxisomal citrate synthase activity may be a valid metabolic engineering strategy for increasing PHA production in other plant species.

A hormonal regulatory module that provides flexibility to tropic responses.

Plants orient their growth depending on directional stimuli such as light and gravity, in a process known as tropic response. Tropisms result from asymmetrical accumulation of auxin across the responding organ relative to the direction of the stimulus, which causes differential growth rates on both sides of the organ. Here, we show that gibberellins (GAs) attenuate the gravitropic reorientation of stimulated hypocotyls of dark-grown Arabidopsis (Arabidopsis thaliana) seedlings. We show that the modulation occurs through induction of the expression of the negative regulator of auxin signaling INDOLE-3-ACETIC ACID INDUCIBLE19/MASSUGU2. The biological significance of this regulatory mechanism involving GAs and auxin seems to be the maintenance of a high degree of flexibility in tropic responses. This notion is further supported by observations that GA-deficient seedlings showed a much lower variance in the response to gravity compared to wild-type seedlings and that the attenuation of gravitropism by GAs resulted in an increased phototropic response. This suggests that the interplay between auxin and GAs may be particularly important for plant orientation under competing tropic stimuli.

A molecular framework for light and gibberellin control of cell elongation.

Cell elongation during seedling development is antagonistically regulated by light and gibberellins (GAs). Light induces photomorphogenesis, leading to inhibition of hypocotyl growth, whereas GAs promote etiolated growth, characterized by increased hypocotyl elongation. The mechanism underlying this antagonistic interaction remains unclear. Here we report on the central role of the Arabidopsis thaliana nuclear transcription factor PIF4 (encoded by PHYTOCHROME INTERACTING FACTOR 4) in the positive control of genes mediating cell elongation and show that this factor is negatively regulated by the light photoreceptor phyB (ref. 4) and by DELLA proteins that have a key repressor function in GA signalling. Our results demonstrate that PIF4 is destabilized by phyB in the light and that DELLAs block PIF4 transcriptional activity by binding the DNA-recognition domain of this factor. We show that GAs abrogate such repression by promoting DELLA destabilization, and therefore cause a concomitant accumulation of free PIF4 in the nucleus. Consistent with this model, intermediate hypocotyl lengths were observed in transgenic plants over-accumulating both DELLAs and PIF4. Destabilization of this factor by phyB, together with its inactivation by DELLAs, constitutes a protein interaction framework that explains how plants integrate both light and GA signals to optimize growth and development in response to changing environments.

Cathepsin B-like and cell death in the unicellular human pathogen Leishmania.

In several studies reporting cell death (CD) in lower eukaryotes and in the human protozoan parasite Leishmania, proteolytic activity was revealed using pan-caspase substrates or inhibitors such as carbobenzoxy-valyl-alanyl-aspartyl-[O-methyl]-fluoromethylketone (Z-VAD-FMK). However, most of the lower eukaryotes do not encode caspase(s) but MCA, which differs from caspase(s) in its substrate specificity and cannot be accountable for the recognition of Z-VAD-FMK. In the present study, we were interested in identifying which enzyme was capturing the Z-VAD substrate. We show that heat shock (HS) induces Leishmania CD and leads to the intracellular binding of Z-VAD-FMK. We excluded binding and inhibition of Z-VAD-FMK to Leishmania major metacaspase (LmjMCA), and identified cysteine proteinase C (LmjCPC), a cathepsin B-like (CPC) enzyme, as the Z-VAD-FMK binding enzyme. We confirmed the specific interaction of Z-VAD-FMK with CPC by showing that Z-VAD binding is absent in a Leishmania mexicana strain in which the cpc gene was deleted. We also show that parasites exposed to various stress conditions release CPC into a soluble fraction. Finally, we confirmed the role of CPC in Leishmania CD by showing that, when exposed to the oxidizing agent hydrogen peroxide (H(2)O(2)), cpc knockout parasites survived better than wild-type parasites (WT). In conclusion, this study identified CPC as the substrate of Z-VAD-FMK in Leishmania and as a potential additional executioner protease in the CD cascade of Leishmania and possibly in other lower eukaryotes.

Peroxisomal Delta(3),Delta(2)-enoyl CoA isomerases and evolution of cytosolic paralogues in embryophytes

Delta(3),Delta(2)-enoyl CoA isomerase (ECI) is an enzyme that participates in the degradation of unsaturated fatty acids through the beta-oxidation cycle. Three genes encoding Delta(3),Delta(2)-enoyl CoA isomerases and named AtECI1, AtECI2 and AtECI3 have been identified in Arabidopsis thaliana. When expressed heterologously in Saccharomyces cerevisiae, all three ECI proteins were targeted to the peroxisomes and enabled the yeast Deltaeci1 mutant to degrade 10Z-heptadecenoic acid, demonstrating Delta(3),Delta(2)-enoyl CoA isomerase activity in vivo. Fusion proteins between yellow fluorescent protein and AtECI1 or AtECI2 were targeted to the peroxisomes in onion epidermal cells and Arabidopsis root cells, but a similar fusion protein with AtECI3 remained in the cytosol for both tissues. AtECI3 targeting to peroxisomes in S. cerevisiae was dependent on yeast PEX5, while expression of Arabidopsis PEX5 in yeast failed to target AtECI3 to peroxisomes. AtECI2 and AtECI3 are tandem duplicated genes and show a high level of amino acid conservation, except at the C-terminus; AtECI2 ends with the well conserved peroxisome targeting signal 1 (PTS1) terminal tripeptide PKL, while AtECI3 possesses a divergent HNL terminal tripeptide. Evolutionary analysis of ECI genes in plants revealed several independent duplication events, with duplications occurring in rice and Medicago truncatula, generating homologues with divergent C-termini and no recognizable PTS1. All plant ECI genes analyzed, including AtECI3, are under negative purifying selection, implying functionality of the cytosolic AtECI3. Analysis of the mammalian and fungal genomes failed to identify cytosolic variants of the Delta(3),Delta(2)-enoyl CoA isomerase, indicating that evolution of cytosolic Delta(3),Delta(2)-enoyl CoA isomerases is restricted to the plant kingdom

Quantifying degradation rates of transmembrane receptor kinases.

Transmembrane receptor-kinases are widespread throughout eukaryotes and their activities are known to regulate all kinds of cellular responses in diverse organs and cell types. In order to guarantee the correct amplitude and duration of signals, receptor levels at the cellular surface need to be tightly controlled. The regulation of receptor degradation is the most direct way to achieve this and elaborate mechanisms are in place to control this process. Therefore, the rate of receptor degradation is a parameter of central importance for understanding the dynamics of a signal transduction cascade. Unfortunately, degradation of transmembrane receptors is a complicated multistep process that involves internalization from the plasma membrane, invagination into the lumen of endosomal compartments, and finally fusion with the vacuole for degradation by vacuolar proteases. Therefore, degradation should be measured in an as noninvasive way as possible, such as not to interfere with the complicated transport processes. Here, a method for minimally invasive, in vivo turn-over measurements in intact organs is provided. This technique was used for quantifying the turn-over rates of the Brassinosteroid receptor kinase BRI1 (BRASSINOSTEROID INSENSITIVE 1) in Arabidopsis thaliana root meristems. Pulse-chase expression of a fluorescently labeled BRI1 variant was used and its turn-over rate was determined by quantitative confocal microscopy. This method is well suited to measure turn-over of transmembrane kinases, but can evidently be extended to measure turn-over of any types of transmembrane proteins.

The fou2 mutation in the major vacuolar cation channel TPC1 confers tolerance to inhibitory luminal calcium.

The SV channel encoded by the TPC1 gene represents a Ca(2+)- and voltage-dependent vacuolar cation channel. Point mutation D454N within TPC1, named fou2 for fatty acid oxygenation upregulated 2, results in increased synthesis of the stress hormone jasmonate. As wounding causes Ca2+ signals and cytosolic Ca2+ is required for SV channel function, we here studied the Ca(2+)-dependent properties of this major vacuolar cation channel with Arabidopsis thaliana mesophyll vacuoles. In patch clamp measurements, wild-type and fou2 SV channels did not exhibit differences in cytosolic Ca2+ sensitivity and Ca2+ impermeability. K+ fluxes through wild-type TPC1 were reduced or even completely faded away when vacuolar Ca2+ reached the 0.1-mm level. The fou2 protein under these conditions, however, remained active. Thus, D454N seems to be part of a luminal Ca2+ recognition site. Thereby the SV channel mutant gains tolerance towards elevated luminal Ca2+. A three-fold higher vacuolar Ca/K ratio in the fou2 mutant relative to wild-type plants seems to indicate that fou2 can accumulate higher levels of vacuolar Ca(2+) before SV channel activity vanishes and K(+) homeostasis is impaired. In response to wounding fou2 plants might thus elicit strong vacuole-derived cytosolic Ca2+ signals resulting in overproduction of jasmonate.

Systems biology of plants : from intraspecific genome variation to computational morphodynamics

Recently, the introduction of second generation sequencing and further advance-ments in confocal microscopy have enabled system-level studies for the functional characterization of genes. The degree of complexity intrinsic to these approaches needs the development of bioinformatics methodologies and computational models for extracting meaningful biological knowledge from the enormous amount of experi¬mental data which is continuously generated. This PhD thesis presents several novel bioinformatics methods and computational models to address specific biological questions in Plant Biology by using the plant Arabidopsis thaliana as a model system. First, a spatio-temporal qualitative analysis of quantitative transcript and protein profiles is applied to show the role of the BREVIS RADIX (BRX) protein in the auxin- cytokinin crosstalk for root meristem growth. Core of this PhD work is the functional characterization of the interplay between the BRX protein and the plant hormone auxin in the root meristem by using a computational model based on experimental evidence. Hyphotesis generated by the modelled to the discovery of a differential endocytosis pattern in the root meristem that splits the auxin transcriptional response via the plasma membrane to nucleus partitioning of BRX. This positional information system creates an auxin transcriptional pattern that deviates from the canonical auxin response and is necessary to sustain the expression of a subset of BRX-dependent auxin-responsive genes to drive root meristem growth. In the second part of this PhD thesis, we characterized the genome-wide impact of large scale deletions on four divergent Arabidopsis natural strains, through the integration of Ultra-High Throughput Sequencing data with data from genomic hybridizations on tiling arrays. Analysis of the identified deletions revealed a considerable portion of protein coding genes affected and supported a history of genomic rearrangements shaped by evolution. In the last part of the thesis, we showed that VIP3 gene in Arabidopsis has an evo-lutionary conserved role in the 3' to 5' mRNA degradation machinery, by applying a novel approach for the analysis of mRNA-Seq data from random-primed mRNA. Altogether, this PhD research contains major advancements in the study of natural genomic variation in plants and in the application of computational morphodynamics models for the functional characterization of biological pathways essential for the plant. - Récemment, l'introduction du séquençage de seconde génération et les avancées dans la microscopie confocale ont permis des études à l'échelle des différents systèmes cellulaires pour la caractérisation fonctionnelle de gènes. Le degrés de complexité intrinsèque à ces approches ont requis le développement de méthodologies bioinformatiques et de modèles mathématiques afin d'extraire de la masse de données expérimentale générée, des information biologiques significatives. Ce doctorat présente à la fois des méthodes bioinformatiques originales et des modèles mathématiques pour répondre à certaines questions spécifiques de Biologie Végétale en utilisant la plante Arabidopsis thaliana comme modèle. Premièrement, une analyse qualitative spatio-temporelle de profiles quantitatifs de transcripts et de protéines est utilisée pour montrer le rôle de la protéine BREVIS RADIX (BRX) dans le dialogue entre l'auxine et les cytokinines, des phytohormones, dans la croissance du méristème racinaire. Le noyau de ce travail de thèse est la caractérisation fonctionnelle de l'interaction entre la protéine BRX et la phytohormone auxine dans le méristème de la racine en utilisant des modèles informatiques basés sur des preuves expérimentales. Les hypothèses produites par le modèle ont mené à la découverte d'un schéma différentiel d'endocytose dans le méristème racinaire qui divise la réponse transcriptionnelle à l'auxine par le partitionnement de BRX de la membrane plasmique au noyau de la cellule. Cette information positionnelle crée une réponse transcriptionnelle à l'auxine qui dévie de la réponse canonique à l'auxine et est nécessaire pour soutenir l'expression d'un sous ensemble de gènes répondant à l'auxine et dépendant de BRX pour conduire la croissance du méristème. Dans la seconde partie de cette thèse de doctorat, nous avons caractérisé l'impact sur l'ensemble du génome des délétions à grande échelle sur quatre souches divergentes naturelles d'Arabidopsis, à travers l'intégration du séquençage à ultra-haut-débit avec l'hybridation génomique sur puces ADN. L'analyse des délétions identifiées a révélé qu'une proportion considérable de gènes codant était affectée, supportant l'idée d'un historique de réarrangement génomique modelé durant l'évolution. Dans la dernière partie de cette thèse, nous avons montré que le gène VÏP3 dans Arabidopsis a conservé un rôle évolutif dans la machinerie de dégradation des ARNm dans le sens 3' à 5', en appliquant une nouvelle approche pour l'analyse des données de séquençage d'ARNm issue de transcripts amplifiés aléatoirement. Dans son ensemble, cette recherche de doctorat contient des avancées majeures dans l'étude des variations génomiques naturelles des plantes et dans l'application de modèles morphodynamiques informatiques pour la caractérisation de réseaux biologiques essentiels à la plante. - Le développement des plantes est écrit dans leurs codes génétiques. Pour comprendre comment les plantes sont capables de s'adapter aux changements environnementaux, il est essentiel d'étudier comment leurs gènes gouvernent leur formation. Plus nous essayons de comprendre le fonctionnement d'une plante, plus nous réalisons la complexité des mécanismes biologiques, à tel point que l'utilisation d'outils et de modèles mathématiques devient indispensable. Dans ce travail, avec l'utilisation de la plante modèle Arabidopsis thalicinci nous avons résolu des problèmes biologiques spécifiques à travers le développement et l'application de méthodes informatiques concrètes. Dans un premier temps, nous avons investigué comment le gène BREVIS RADIX (BRX) régule le développement de la racine en contrôlant la réponse à deux hormones : l'auxine et la cytokinine. Nous avons employé une analyse statistique sur des mesures quantitatives de transcripts et de produits de gènes afin de démontrer que BRX joue un rôle antagonisant dans le dialogue entre ces deux hormones. Lorsque ce-dialogue moléculaire est perturbé, la racine primaire voit sa longueur dramatiquement réduite. Pour comprendre comment BRX répond à l'auxine, nous avons développé un modèle informatique basé sur des résultats expérimentaux. Les simulations successives ont mené à la découverte d'un signal positionnel qui contrôle la réponse de la racine à l'auxine par la régulation du mouvement intracellulaire de BRX. Dans la seconde partie de cette thèse, nous avons analysé le génome entier de quatre souches naturelles d'Arabidopsis et nous avons trouvé qu'une grande partie de leurs gènes étaient manquant par rapport à la souche de référence. Ce résultat indique que l'historique des modifications génomiques conduites par l'évolution détermine une disponibilité différentielle des gènes fonctionnels dans ces plantes. Dans la dernière partie de ce travail, nous avons analysé les données du transcriptome de la plante où le gène VIP3 était non fonctionnel. Ceci nous a permis de découvrir le rôle double de VIP3 dans la régulation de l'initiation de la transcription et dans la dégradation des transcripts. Ce rôle double n'avait jusqu'alors été démontrée que chez l'homme. Ce travail de doctorat supporte le développement et l'application de méthodologies informatiques comme outils inestimables pour résoudre la complexité des problèmes biologiques dans la recherche végétale. L'intégration de la biologie végétale et l'informatique est devenue de plus en plus importante pour l'avancée de nos connaissances sur le fonctionnement et le développement des plantes.

Insect eggs suppress plant defence against chewing herbivores.

Plants activate direct and indirect defences in response to insect egg deposition. However, whether eggs can manipulate plant defence is unknown. In Arabidopsis thaliana, oviposition by the butterfly Pieris brassicae triggers cellular and molecular changes that are similar to the changes caused by biotrophic pathogens. In the present study, we found that the plant defence signal salicylic acid (SA) accumulates at the site of oviposition. This is unexpected, as the SA pathway controls defence against fungal and bacterial pathogens and negatively interacts with the jasmonic acid (JA) pathway, which is crucial for the defence against herbivores. Application of P. brassicae or Spodoptera littoralis egg extract onto leaves reduced the induction of insect-responsive genes after challenge with caterpillars, suggesting that egg-derived elicitors suppress plant defence. Consequently, larval growth of the generalist herbivore S. littoralis, but not of the specialist P. brassicae, was significantly higher on plants treated with egg extract than on control plants. In contrast, suppression of gene induction and enhanced S. littoralis performance were not seen in the SA-deficient mutant sid2-1, indicating that it is SA that mediates this phenomenon. These data reveal an intriguing facet of the cross-talk between SA and JA signalling pathways, and suggest that insects have evolved a way to suppress the induction of defence genes by laying eggs that release elicitors. We show here that egg-induced SA accumulation negatively interferes with the JA pathway, and provides an advantage for generalist herbivores.

The cuticle : not only a barrier for plant defence: A novel defence syndrome in plants with cuticular defects.

The cuticle is a physical barrier that prevents water loss and protects against irradiation, xenobiotics and pathogens. This classic textbook statement has recently been revisited and several observations were made showing that this dogma falls short of being universally true. Both transgenic Arabidopsis thaliana lines expressing cell wall-targeted fungal cutinase (so-called CUTE plants) or lipase as well as several A. thaliana mutants with altered cuticular structure remained free of symptoms after an inoculation with Botrytis cinerea. The alterations in cuticular structure lead to the release of fungitoxic substances and changes in gene expression that form a multifactorial defence response. Several models to explain this syndrome are discussed.

Plasma membrane cyclic nucleotide gated calcium channels control land plant thermal sensing and acquired thermotolerance.

Typically at dawn on a hot summer day, land plants need precise molecular thermometers to sense harmless increments in the ambient temperature to induce a timely heat shock response (HSR) and accumulate protective heat shock proteins in anticipation of harmful temperatures at mid-day. Here, we found that the cyclic nucleotide gated calcium channel (CNGC) CNGCb gene from Physcomitrella patens and its Arabidopsis thaliana ortholog CNGC2, encode a component of cyclic nucleotide gated Ca(2+) channels that act as the primary thermosensors of land plant cells. Disruption of CNGCb or CNGC2 produced a hyper-thermosensitive phenotype, giving rise to an HSR and acquired thermotolerance at significantly milder heat-priming treatments than in wild-type plants. In an aequorin-expressing moss, CNGCb loss-of-function caused a hyper-thermoresponsive Ca(2+) influx and altered Ca(2+) signaling. Patch clamp recordings on moss protoplasts showed the presence of three distinct thermoresponsive Ca(2+) channels in wild-type cells. Deletion of CNGCb led to a total absence of one and increased the open probability of the remaining two thermoresponsive Ca(2+) channels. Thus, CNGC2 and CNGCb are expected to form heteromeric Ca(2+) channels with other related CNGCs. These channels in the plasma membrane respond to increments in the ambient temperature by triggering an optimal HSR, leading to the onset of plant acquired thermotolerance.

Characterization of the rice PHO1 gene family reveals a key role for OsPHO1;2 in phosphate homeostasis and the evolution of a distinct clade in dicotyledons.

Phosphate homeostasis was studied in a monocotyledonous model plant through the characterization of the PHO1 gene family in rice (Oryza sativa). Bioinformatics and phylogenetic analysis showed that the rice genome has three PHO1 homologs, which cluster with the Arabidopsis (Arabidopsis thaliana) AtPHO1 and AtPHO1;H1, the only two genes known to be involved in root-to-shoot transfer of phosphate. In contrast to the Arabidopsis PHO1 gene family, all three rice PHO1 genes have a cis-natural antisense transcript located at the 5 ' end of the genes. Strand-specific quantitative reverse transcription-PCR analyses revealed distinct patterns of expression for sense and antisense transcripts for all three genes, both at the level of tissue expression and in response to nutrient stress. The most abundantly expressed gene was OsPHO1;2 in the roots, for both sense and antisense transcripts. However, while the OsPHO1;2 sense transcript was relatively stable under various nutrient deficiencies, the antisense transcript was highly induced by inorganic phosphate (Pi) deficiency. Characterization of Ospho1;1 and Ospho1;2 insertion mutants revealed that only Ospho1;2 mutants had defects in Pi homeostasis, namely strong reduction in Pi transfer from root to shoot, which was accompanied by low-shoot and high-root Pi. Our data identify OsPHO1;2 as playing a key role in the transfer of Pi from roots to shoots in rice, and indicate that this gene could be regulated by its cis-natural antisense transcripts. Furthermore, phylogenetic analysis of PHO1 homologs in monocotyledons and dicotyledons revealed the emergence of a distinct clade of PHO1 genes in dicotyledons, which include members having roles other than long-distance Pi transport.