In planta and in silico characterization of five sesquiterpene synthases from Vitis vinifera (cv. Shiraz) berries.

MAIN CONCLUSION: Five Vitis vinifera sesquiterpene synthases were characterized, two was previously uncharacterized, one being a caryophyllene/cubebene synthase and the other a cadinene synthase. Residue differences with other Vitis sesquiterpene synthases are described. The biochemical composition of grape berries at harvest can have a profound effect on the varietal character of the wine produced. Sesquiterpenes are an important class of volatile compounds produced in grapes that contribute to the flavor and aroma of wine, making the elucidation of their biosynthetic origin an important field of research. Five cDNAs corresponding to sesquiterpene synthase genes (TPSs) were isolated from Shiraz berries and expressed in planta in Nicotiana benthamiana followed by chemical characterization by GC-MS. Three of the TPS cDNAs were isolated from immature berries and two were isolated from ripe Shiraz berries. Two of the investigated enzymes, TPS26 and TPS27, have been previously investigated by expression in E. coli, and the in planta products generally correspond to these previous studies. The enzyme TPS07 differed by eight amino acids (none of which are in the active site) from germacrene B and D synthase isolated from Gewürztraminer grapes and characterized in vitro. Here in planta characterization of VvShirazTPS07 yielded ylangene, germacrene D and several minor products. Two of the enzymes isolated from immature berries were previously uncharacterized enzymes. VvShirazTPS-Y1 produced cadinene as a major product and at least 17 minor sesquiterpenoid skeletons. The second, VvShirazTPS-Y2, was characterized as a caryophyllene/cubebene synthase, a combination of products not previously reported from a single enzyme. Using in silico methods, we identified residues that could play key roles regarding differences in product formation of these enzymes. The first ring closure that is either a 1,10- or 1,11-ring closure is likely controlled by three neighboring amino acids in helices G1, H2, and J. As for many other investigated TPS enzymes, we also observe that only a few residues can account for radical changes in product formation.

Two key polymorphisms in a newly discovered allele of the Vitis vinifera TPS24 gene are responsible for the production of the rotundone precursor α-guaiene.

Rotundone was initially identified as a grape-derived compound responsible for the peppery aroma of Shiraz wine varieties. It has subsequently been found in black and white pepper and several other spices. Because of its potent aroma, the molecular basis for rotundone formation is of particular relevance to grape and wine scientists and industry. We have identified and functionally characterized in planta a sesquiterpene synthase, VvGuaS, from developing grape berries, and have demonstrated that it produces the precursor of rotundone, α-guaiene, as its main product. The VvGuaS enzyme is a novel allele of the sesquiterpene synthase gene, VvTPS24, which has previously been reported to encode VvPNSeInt, an enzyme that produces a variety of selinene-type sesquiterpenes. This newly discovered VvTPS24 allele encodes an enzyme 99.5% identical to VvPNSeInt, with the differences comprising just 6 out of the 561 amino acid residues. Molecular modelling of the enzymes revealed that two of these residues, T414 and V530, are located in the active site of VvGuaS within 4 Å of the binding-site of the substrate, farnesyl pyrophosphate. Mutation of these two residues of VvGuaS into the corresponding polymorphisms in VvPNSeInt results in a complete functional conversion of one enzyme into the other, while mutation of each residue individually produces an intermediate change in the product profile. We have therefore demonstrated that VvGuaS, an enzyme responsible for production of the rotundone precursor, α-guaiene, is encoded by a novel allele of the previously characterized grapevine gene VvTPS24 and that two specific polymorphisms are responsible for functional differences between VvTPS24 alleles.

Evolution of substrate recognition sites (SRSs) in cytochromes P450 from Apiaceae exemplified by the CYP71AJ subfamily.

BACKGROUND: Large proliferations of cytochrome P450 encoding genes resulting from gene duplications can be termed as 'blooms', providing genetic material for the genesis and evolution of biosynthetic pathways. Furanocoumarins are allelochemicals produced by many of the species in Apiaceaous plants belonging to the Apioideae subfamily of Apiaceae and have been described as being involved in the defence reaction against phytophageous insects. RESULTS: A bloom in the cytochromes P450 CYP71AJ subfamily has been identified, showing at least 2 clades and 6 subclades within the CYP71AJ subfamily. Two of the subclades were functionally assigned to the biosynthesis of furanocoumarins. Six substrate recognition sites (SRS1-6) important for the enzymatic conversion were investigated in the described cytochromes P450 and display significant variability within the CYP71AJ subfamily. Homology models underline a significant modification of the accession to the iron atom, which might explain the difference of the substrate specificity between the cytochromes P450 restricted to furanocoumarins as substrates and the orphan CYP71AJ. CONCLUSION: Two subclades functionally assigned to the biosynthesis of furanocoumarins and four other subclades were identified and shown to be part of two distinct clades within the CYP71AJ subfamily. The subclades show significant variability within their substrate recognition sites between the clades, suggesting different biochemical functions and providing insights into the evolution of cytochrome P450 'blooms' in response to environmental pressures.

Engineering mammalian mucin-type O-glycosylation in plants.

Mucin-type O-glycosylation is an important post-translational modification that confers a variety of biological properties and functions to proteins. This post-translational modification has a particularly complex and differentially regulated biosynthesis rendering prediction and control of where O-glycans are attached to proteins, and which structures are formed, difficult. Because plants are devoid of GalNAc-type O-glycosylation, we have assessed requirements for establishing human GalNAc O-glycosylation de novo in plants with the aim of developing cell systems with custom-designed O-glycosylation capacity. Transient expression of a Pseudomonas aeruginosa Glc(NAc) C4-epimerase and a human polypeptide GalNAc-transferase in leaves of Nicotiana benthamiana resulted in GalNAc O-glycosylation of co-expressed human O-glycoprotein substrates. A chimeric YFP construct containing a 3.5 tandem repeat sequence of MUC1 was glycosylated with up to three and five GalNAc residues when co-expressed with GalNAc-T2 and a combination of GalNAc-T2 and GalNAc-T4, respectively, as determined by mass spectrometry. O-Glycosylation was furthermore demonstrated on a tandem repeat of MUC16 and interferon α2b. In plants, prolines in certain classes of proteins are hydroxylated and further substituted with plant-specific O-glycosylation; unsubstituted hydroxyprolines were identified in our MUC1 construct. In summary, this study demonstrates that mammalian type O-glycosylation can be established in plants and that plants may serve as a host cell for production of recombinant O-glycoproteins with custom-designed O-glycosylation. The observed hydroxyproline modifications, however, call for additional future engineering efforts.

VTCdb: a gene co-expression database for the crop species Vitis vinifera (grapevine).

BACKGROUND: Gene expression datasets in model plants such as Arabidopsis have contributed to our understanding of gene function and how a single underlying biological process can be governed by a diverse network of genes. The accumulation of publicly available microarray data encompassing a wide range of biological and environmental conditions has enabled the development of additional capabilities including gene co-expression analysis (GCA). GCA is based on the understanding that genes encoding proteins involved in similar and/or related biological processes may exhibit comparable expression patterns over a range of experimental conditions, developmental stages and tissues. We present an open access database for the investigation of gene co-expression networks within the cultivated grapevine, Vitis vinifera. DESCRIPTION: The new gene co-expression database, VTCdb (http://vtcdb.adelaide.edu.au/Home.aspx), offers an online platform for transcriptional regulatory inference in the cultivated grapevine. Using condition-independent and condition-dependent approaches, grapevine co-expression networks were constructed using the latest publicly available microarray datasets from diverse experimental series, utilising the Affymetrix Vitis vinifera GeneChip (16 K) and the NimbleGen Grape Whole-genome microarray chip (29 K), thus making it possible to profile approximately 29,000 genes (95% of the predicted grapevine transcriptome). Applications available with the online platform include the use of gene names, probesets, modules or biological processes to query the co-expression networks, with the option to choose between Affymetrix or Nimblegen datasets and between multiple co-expression measures. Alternatively, the user can browse existing network modules using interactive network visualisation and analysis via CytoscapeWeb. To demonstrate the utility of the database, we present examples from three fundamental biological processes (berry development, photosynthesis and flavonoid biosynthesis) whereby the recovered sub-networks reconfirm established plant gene functions and also identify novel associations. CONCLUSIONS: Together, we present valuable insights into grapevine transcriptional regulation by developing network models applicable to researchers in their prioritisation of gene candidates, for on-going study of biological processes related to grapevine development, metabolism and stress responses.

Toward stable genetic engineering of human O-glycosylation in plants.

Glycosylation is the most abundant and complex posttranslational modification to be considered for recombinant production of therapeutic proteins. Mucin-type (N-acetylgalactosamine [GalNAc]-type) O-glycosylation is found in eumetazoan cells but absent in plants and yeast, making these cell types an obvious choice for de novo engineering of this O-glycosylation pathway. We previously showed that transient implementation of O-glycosylation capacity in plants requires introduction of the synthesis of the donor substrate UDP-GalNAc and one or more polypeptide GalNAc-transferases for incorporating GalNAc residues into proteins. Here, we have stably engineered O-glycosylation capacity in two plant cell systems, soil-grown Arabidopsis (Arabidopsis thaliana) and tobacco (Nicotiana tabacum) Bright Yellow-2 suspension culture cells. Efficient GalNAc O-glycosylation of two stably coexpressed substrate O-glycoproteins was obtained, but a high degree of proline hydroxylation and hydroxyproline-linked arabinosides, on a mucin (MUC1)-derived substrate, was also observed. Addition of the prolyl 4-hydroxylase inhibitor 2,2-dipyridyl, however, effectively suppressed proline hydroxylation and arabinosylation of MUC1 in Bright Yellow-2 cells. In summary, stably engineered mammalian type O-glycosylation was established in transgenic plants, demonstrating that plants may serve as host cells for the production of recombinant O-glycoproteins. However, the present stable implementation further strengthens the notion that elimination of endogenous posttranslational modifications may be needed for the production of protein therapeutics.

Identification and characterization of a kunzeaol synthase from Thapsia garganica: implications for the biosynthesis of the pharmaceutical thapsigargin.

Thapsigargin is a major terpenoid constituent of Thapsia garganica root. Owing to its potent antagonistic effect on the sarcoplasmic/endoplasmic reticulum Ca2+-ATPase, thapsigargin has been widely used to study Ca2+ signalling and is also a potential drug for prostate cancer. Despite its importance, thapsigargin biosynthesis in T. garganica remains unknown. In order to decipher thapsigargin biosynthesis, deep transcript sequencing (454 and Illumina) of the T. garganica root was performed, and two terpene synthases (TgTPS1/2) were identified. Functional characterization of their encoded enzymes in a metabolically engineered yeast revealed that TgTPS1 synthesized δ-cadinene, whereas TgTPS2 produced ten distinct terpenoids. However, cultivation of the TgTPS2-expressing yeast in pH-maintained conditions (pH 6-7) yielded one major oxygenated sesquiterpenoid, suggesting that formation of multiple terpenoids was caused by acidity. The major terpene product from TgTPS2 was identified as 6ß-hydroxygermacra-1(10),4-diene (kunzeaol) by mass-fragmentation pattern, retention index, the nature of its acid-induced degradation and NMR. Also, recombinant TgTPS2 efficiently catalysed the synthesis of kunzeaol in vitro from farnesyl diphosphate with a Km of 2.6 µM and a kcat of 0.03 s-1. The present paper is the first report of a kunzeaol synthase, and a mechanism for the transformation of kunzeaol into the thapsigargin backbone is proposed.

Transcriptome analysis of Thapsia laciniata Rouy provides insights into terpenoid biosynthesis and diversity in Apiaceae.

Thapsia laciniata Rouy (Apiaceae) produces irregular and regular sesquiterpenoids with thapsane and guaiene carbon skeletons, as found in other Apiaceae species. A transcriptomic analysis utilizing Illumina next-generation sequencing enabled the identification of novel genes involved in the biosynthesis of terpenoids in Thapsia. From 66.78 million HQ paired-end reads obtained from T. laciniata roots, 64.58 million were assembled into 76,565 contigs (N50: 1261 bp). Seventeen contigs were annotated as terpene synthases and five of these were predicted to be sesquiterpene synthases. Of the 67 contigs annotated as cytochromes P450, 18 of these are part of the CYP71 clade that primarily performs hydroxylations of specialized metabolites. Three contigs annotated as aldehyde dehydrogenases grouped phylogenetically with the characterized ALDH1 from Artemisia annua and three contigs annotated as alcohol dehydrogenases grouped with the recently described ADH1 from A. annua. ALDH1 and ADH1 were characterized as part of the artemisinin biosynthesis. We have produced a comprehensive EST dataset for T. laciniata roots, which contains a large sample of the T. laciniata transcriptome. These transcriptome data provide the foundation for future research into the molecular basis for terpenoid biosynthesis in Thapsia and on the evolution of terpenoids in Apiaceae.

Structural and functional analyses of PpENA1 provide insights into cation binding by type IID P-type ATPases in lower plants and fungi.

PpENA1 is a membrane-spanning transporter from the moss Physcomitrella patens, and is the first type IID P-type ATPase to be reported in the plant kingdom. In Physcomitrella, PpENA1 is essential for normal growth under moderate salt stress, while in yeast, type IID ATPases provide a vital efflux mechanism for cells under high salt conditions by selectively transporting Na+ or K+ across the plasma membrane. To investigate the structural basis for cation-binding within the type IID ATPase subfamily, we used homology modeling to identify a highly conserved cation-binding pocket between membrane helix (MH) 4 and MH 6 of the membrane-spanning pore of PpENA1. Mutation of specific charged and polar residues on MHs 4-6 resulted in a decrease or loss of protein activity as measured by complementation assays in yeast. The E298S mutation on MH 4 of PpENA1 had the most significant effect on activity despite the presence of a serine at this position in fungal type IID ATPases. Activity was partially restored in an inactivated PpENA1 mutant by the insertion of two additional serine residues on MH 4 and one on MH 6 based on the presence of these residues in fungal type IID ATPases. Our results suggest that the residues responsible for cation-binding in PpENA1 are distinct from those in fungal type IID ATPases, and that a fungal-type cation binding site can be successfully engineered into the moss protein.

Transcriptome analysis at four developmental stages of grape berry (Vitis vinifera cv. Shiraz) provides insights into regulated and coordinated gene expression.

BACKGROUND: Vitis vinifera berry development is characterised by an initial phase where the fruit is small, hard and acidic, followed by a lag phase known as veraison. In the final phase, berries become larger, softer and sweeter and accumulate an array of organoleptic compounds. Since the physiological and biochemical makeup of grape berries at harvest has a profound impact on the characteristics of wine, there is great interest in characterising the molecular and biophysical changes that occur from flowering through veraison and ripening, including the coordination and temporal regulation of metabolic gene pathways. Advances in deep-sequencing technologies, combined with the availability of increasingly accurate V. vinifera genomic and transcriptomic data, have enabled us to carry out RNA-transcript expression analysis on a global scale at key points during berry development. RESULTS: A total of 162 million 100-base pair reads were generated from pooled Vitis vinifera (cv. Shiraz) berries sampled at 3-weeks post-anthesis, 10- and 11-weeks post-anthesis (corresponding to early and late veraison) and at 17-weeks post-anthesis (harvest). Mapping reads from each developmental stage (36-45 million) onto the NCBI RefSeq transcriptome of 23,720 V. vinifera mRNAs revealed that at least 75% of these transcripts were detected in each sample. RNA-Seq analysis uncovered 4,185 transcripts that were significantly upregulated at a single developmental stage, including 161 transcription factors. Clustering transcripts according to distinct patterns of transcription revealed coordination in metabolic pathways such as organic acid, stilbene and terpenoid metabolism. From the phenylpropanoid/stilbene biosynthetic pathway at least 46 transcripts were upregulated in ripe berries when compared to veraison and immature berries, and 12 terpene synthases were predominantly detected only in a single sample. Quantitative real-time PCR was used to validate the expression pattern of 12 differentially expressed genes from primary and secondary metabolic pathways. CONCLUSIONS: In this study we report the global transcriptional profile of Shiraz grapes at key stages of development. We have undertaken a comprehensive analysis of gene families contributing to commercially important berry characteristics and present examples of co-regulation and differential gene expression. The data reported here will provide an invaluable resource for the on-going molecular investigation of wine grapes.

A comparison of headspace solid-phase microextraction and classic hydrodistillation for the identification of volatile constituents from Thapsia spp. provides insights into guaianolide biosynthesis in Apiaceae.

INTRODUCTION: Thapsia spp. (Apiaceae) are the major natural source of polyoxygenated guaianolide sesquiterpene lactones known as thapsigargins, which induce apoptosis in mammalian cells via a high affinity inhibition of the sarco/endoplasmic reticulum Ca(2+) ATPase. The mechanism of biosynthesis of thapsigargins has not been elucidated, and probable biochemical precursors such as hydrocarbon or oxygenated sesquiterpenes have not been identified in previous phytochemical analyses of essential oils from this genus. OBJECTIVE: To investigate the utility of solid phase micro-extraction (SPME), when compared with classical essential oil distillates, for identifying potential precursors of guaianolide sesquiterpene lactones from Thapsia garganica L. and Thapsia villosa L. type II. METHODOLOGY: A systematic description of the volatile components of roots, flowers, stems and fruits of T. villosa and of root, flower and fruits of T. garganica was constructed via GC-MS analyses of SPME-adsorbed compounds and of essential oils obtained through hydrodistillation of the same tissues. RESULTS: The sesquiterpenoids δ-cadinene, α- and δ-guaiene, elemol and guaiols were found to be major volatile constituents of the roots of T. garganica and T. villosa trapped using SPME. In contrast, these sesquiterpenoids were not detected or were at negligible levels in essential oils, where sesquiterpenoids are potentially converted to azulenes during hydrodistillation. CONCLUSION: The new data reported in this study demonstrates that SPME is a valuable tool for the identification of volatile sesquiterpenes when compared with analysis of essential oils, and we postulate that guaiene is the likely precursor of guaianolide sesquiterpenes from Thapsia.

Stabilization of dhurrin biosynthetic enzymes from Sorghum bicolor using a natural deep eutectic solvent.

In recent years, ionic liquids and deep eutectic solvents (DESs) have gained increasing attention due to their ability to extract and solubilize metabolites and biopolymers in quantities far beyond their solubility in oil and water. The hypothesis that naturally occurring metabolites are able to form a natural deep eutectic solvent (NADES), thereby constituting a third intracellular phase in addition to the aqueous and lipid phases, has prompted researchers to study the role of NADES in living systems. As an excellent solvent for specialized metabolites, formation of NADES in response to dehydration of plant cells could provide an appropriate environment for the functional storage of enzymes during drought. Using the enzymes catalyzing the biosynthesis of the defense compound dhurrin as an experimental model system, we demonstrate that enzymes involved in this pathway exhibit increased stability in NADES compared with aqueous buffer solutions, and that enzyme activity is restored upon rehydration. Inspired by nature, application of NADES provides a biotechnological approach for long-term storage of entire biosynthetic pathways including membrane-anchored enzymes.

Changes in the Composition of the Gut Microbiota and the Blood Transcriptome in Preterm Infants at Less than 29 Weeks Gestation Diagnosed with Bronchopulmonary Dysplasia.

Bronchopulmonary dysplasia (BPD) is a common chronic lung condition in preterm infants that results in abnormal lung development and leads to considerable morbidity and mortality, making BPD one of the most common complications of preterm birth. We employed RNA sequencing and 16S rRNA gene sequencing to profile gene expression in blood and the composition of the fecal microbiota in infants born at <29 weeks gestational age and diagnosed with BPD in comparison to those of preterm infants that were not diagnosed with BPD. 16S rRNA gene sequencing, performed longitudinally on 255 fecal samples collected from 50 infants in the first months of life, identified significant differences in the relative levels of abundance of Klebsiella, Salmonella, Escherichia/Shigella, and Bifidobacterium in the BPD infants in a manner that was birth mode dependent. Transcriptome sequencing (RNA-Seq) analysis revealed that more than 400 genes were upregulated in infants with BPD. Genes upregulated in BPD infants were significantly enriched for functions related to red blood cell development and oxygen transport, while several immune-related pathways were downregulated. We also identified a gene expression signature consistent with an enrichment of immunosuppressive CD71+ early erythroid cells in infants with BPD. Intriguingly, genes that were correlated in their expression with the relative abundances of specific taxa in the microbiota were significantly enriched for roles in the immune system, suggesting that changes in the microbiota might influence immune gene expression systemically.IMPORTANCE Bronchopulmonary dysplasia (BPD) is a serious inflammatory condition of the lung and is the most common complication associated with preterm birth. A large body of evidence now suggests that the gut microbiota can influence immunity and inflammation systemically; however, the role of the gut microbiota in BPD has not been evaluated to date. Here, we report that there are significant differences in the gut microbiota of infants born at <29 weeks gestation and subsequently diagnosed with BPD, which are particularly pronounced when infants are stratified by birth mode. We also show that erythroid and immune gene expression levels are significantly altered in BPD infants. Interestingly, we identified an association between the composition of the microbiota and immune gene expression in blood in early life. Together, these findings suggest that the composition of the microbiota may influence the risk of developing BPD and, more generally, may shape systemic immune gene expression.

Perspectives on using physcomitrella patens as an alternative production platform for thapsigargin and other terpenoid drug candidates.

To overcome the potential future demand for terpenoids used as drugs, a new production platform is currently being established in our laboratory. The moss Physcomitrella has been chosen as the candidate organism for production of drug candidates based on terpenoids derived from plants, with a primary focus on the sesquiterpene lactone, thapsigargin. This drug candidate and other candidates/drugs with sesquiterpene skeleton are difficult to obtain by chemical synthesis due to their large number of chiral centers. Furthermore, they are not available in sufficient amounts from their original plant. The requirement for a new production system to meet the potential market demand for these compounds is not only obvious, but also essential if sufficient quantities of the drug candidates are to be available for the potential therapeutic use.

Heterologous expression of cDNAs encoding monodehydroascorbate reductases from the moss, Physcomitrella patens and characterization of the expressed enzymes.

Monodehydroascorbate reductase (MDHAR; EC 1.6.5.4) catalyses the reduction of the monodehydroascorbate (MDHA) radical to ascorbate, using NADH or NADPH as an electron donor, and is believed to be involved in maintaining the reactive oxygen scavenging capability of plant cells. This key enzyme in the ascorbate-glutathione cycle has been studied here in the moss Physcomitrella patens, which is tolerant to a range of abiotic stresses and is increasingly used as a model plant. In the present study, three cDNAs encoding different MDHAR isoforms of 47 kDa were identified in P. patens, and found to exhibit enzymic characteristics similar to MDHARs in vascular plants despite low-sequence identity and a distant evolutionary relationship between the species. The three cDNAs for the P. patens MDHAR enzymes were expressed in Escherichia coli and the active enzymes were purified and characterized. Each recombinant protein displayed an absorbance spectrum typical of flavoenzymes and contained a single non-covalently bound FAD coenzyme molecule. The Km and kcat values for the heterologously expressed PpMDHAR enzymes ranged from 8 to 18 microM and 120-130 s(-1), respectively, using NADH as the electron donor. The Km values were at least an order of magnitude higher for NADPH. The Km values for the MDHA radical were approximately 0.5-1.0 microM for each of the purified enzymes, and further kinetic analyses indicated that PpMDHARs follow a 'ping-pong' kinetic mechanism. In contrast to previously published data, site-directed mutagenesis indicated that the conserved cysteine residue is not directly involved in the reduction of MDHA.

Exclusion of Na+ via sodium ATPase (PpENA1) ensures normal growth of Physcomitrella patens under moderate salt stress.

The bryophyte Physcomitrella patens is unlike any other plant identified to date in that it possesses a gene that encodes an ENA-type Na(+)-ATPase. To complement previous work in yeast (Saccharomyces cerevisiae), we determined the importance of having a Na(+)-ATPase in planta by conducting physiological analyses of PpENA1 in Physcomitrella. Expression studies showed that PpENA1 is up-regulated by NaCl and, to a lesser degree, by osmotic stress. Maximal induction is obtained after 8 h at 60 mm NaCl or above. No other abiotic stress tested led to significant increases in PpENA1 expression. In the gametophyte, strong expression was confined to the rhizoids, stem, and the basal part of the leaf. In the protonemata, expression was ubiquitous with a few filaments showing stronger expression. At 100 mm NaCl, wild-type plants were able to maintain a higher K(+)-to-Na(+) ratio than the PpENA1 (ena1) knockout gene, but at higher NaCl concentrations no difference was observed. Although no difference in chlorophyll content was observed between ena1 and wild type at 100 mm NaCl, the impaired Na(+) exclusion in ena1 plants led to an approximately 40% decrease in growth.

Gene structure and expression pattern analysis of three monodehydroascorbate reductase (Mdhar) genes in Physcomitrella patens: implications for the evolution of the MDHAR family in plants.

The ascorbate-glutathione pathway plays a major role in the detoxification of reactive oxygen species (ROS) in vascular plants. One of the key enzymes in this pathway is monodehydroascorbate reductase (MDHAR), a FAD enzyme that catalyses the reduction of the monodehydroascorbate radical. To elucidate the evolution and functional role of MDHAR we identified and characterised MDHARs from the moss Physcomitrella patens. Expressed sequence tag (EST) databases containing approximately 100.000 ESTs from Physcomitrella were searched and three isoforms of monodehydroascorbate reductase (PpMDHAR1, PpMDHAR2 and PpMDHAR3) were identified. In vascular plants MDHAR is found in the cytosol, chloroplast, mitochondria and peroxisome. Surprisingly, all three PpMDHARs resembled the cytosolic isoforms from vascular plants lacking the NH(2)-terminal or COOH-terminal extension found in organelle targeted MDHARs. The number and position of introns was also conserved between PpMDHARs and cytosolic MDHARs from vascular plants. Phylogenetic analysis revealed that cytosolic MDHARs are monophyletic in origin and the ancestral gene evolved before the divergence of bryophytes more than 400 million years ago. Transcript analyses showed that expression of PpMdhar1 and PpMdhar3 was increased up to 5-fold under salt stress, osmotic stress or upon exposure to abscisic acid. In contrast, PpMdhar transcription levels were unchanged upon chilling, UV-B exposure or oxidative stress. The conservation of cytosolic MDHAR in the land-plant lineage and the transcriptional upregulation under water deficiency suggest that the evolution of cytosolic MDHAR played an essential role in stress protection for land plants when they inhabited the dry terrestrial environment.