Patterns of Metabolite Changes Identified from Large-Scale Gene Perturbations in Arabidopsis Using a Genome-Scale Metabolic Network

Metabolomics enables quantitative evaluation of metabolic changes caused by genetic or environmental perturbations. However, little is known about how perturbing a single gene changes the metabolic system as a whole and which network and functional properties are involved in this response. To answer this question, we investigated the metabolite profiles from 136 mutants with single gene perturbations of functionally diverse Arabidopsis (Arabidopsis thaliana) genes. Fewer than 10 metabolites were changed significantly relative to the wild type in most of the mutants, indicating that the metabolic network was robust to perturbations of single metabolic genes. These changed metabolites were closer to each other in a genome-scale metabolic network than expected by chance, supporting the notion that the genetic perturbations changed the network more locally than globally. Surprisingly, the changed metabolites were close to the perturbed reactions in only 30% of the mutants of the well-characterized genes. To determine the factors that contributed to the distance between the observed metabolic changes and the perturbation site in the network, we examined nine network and functional properties of the perturbed genes. Only the isozyme number affected the distance between the perturbed reactions and changed metabolites. This study revealed patterns of metabolic changes from large-scale gene perturbations and relationships between characteristics of the perturbed genes and metabolic changes.Rational and quantitative assessment of metabolic changes in response to genetic modification (GM) is an open question and in need of innovative solutions. Nontargeted metabolite profiling can detect thousands of compounds, but it is not easy to understand the significance of the changed metabolites in the biochemical and biological context of the organism. To better assess the changes in metabolites from nontargeted metabolomics studies, it is important to examine the changed metabolites in the context of the genome-scale metabolic network of the organism.Metabolomics is a technique that aims to quantify all the metabolites in a biological system (Nikolau and Wurtele, 2007; Nicholson and Lindon, 2008; Roessner and Bowne, 2009). It has been used widely in studies ranging from disease diagnosis (Holmes et al., 2008; DeBerardinis and Thompson, 2012) and drug discovery (Cascante et al., 2002; Kell, 2006) to metabolic reconstruction (Feist et al., 2009; Kim et al., 2012) and metabolic engineering (Keasling, 2010; Lee et al., 2011). Metabolomic studies have demonstrated the possibility of identifying gene functions from changes in the relative concentrations of metabolites (metabotypes or metabolic signatures; Ebbels et al., 2004) in various species including yeast (Saccharomyces cerevisiae; Raamsdonk et al., 2001; Allen et al., 2003), Arabidopsis (Arabidopsis thaliana; Brotman et al., 2011), tomato (Solanum lycopersicum; Schauer et al., 2006), and maize (Zea mays; Riedelsheimer et al., 2012). Metabolomics has also been used to better understand how plants interact with their environments (Field and Lake, 2011), including their responses to biotic and abiotic stresses (Dixon et al., 2006; Arbona et al., 2013), and to predict important agronomic traits (Riedelsheimer et al., 2012). Metabolite profiling has been performed on many plant species, including angiosperms such as Arabidopsis, poplar (Populus trichocarpa), and Catharanthus roseus (Sumner et al., 2003; Rischer et al., 2006), basal land plants such as Selaginella moellendorffii and Physcomitrella patens (Erxleben et al., 2012; Yobi et al., 2012), and Chlamydomonas reinhardtii (Fernie et al., 2012; Davis et al., 2013). With the availability of whole genome sequences of various species, metabolomics has the potential to become a useful tool for elucidating the functions of genes using large-scale systematic analyses (Fiehn et al., 2000; Saito and Matsuda, 2010; Hur et al., 2013).Although metabolomics data have the potential for identifying the roles of genes that are associated with metabolic phenotypes, the biochemical mechanisms that link functions of genes with metabolic phenotypes are still poorly characterized. For example, we do not yet know the principles behind how perturbing the expression of a single gene changes the metabolic system as a whole. Large-scale metabolomics data have provided useful resources for linking phenotypes to genotypes (Fiehn et al., 2000; Roessner et al., 2001; Tikunov et al., 2005; Schauer et al., 2006; Lu et al., 2011; Fukushima et al., 2014). For example, Lu et al. (2011) compared morphological and metabolic phenotypes from more than 5,000 Arabidopsis chloroplast mutants using gas chromatography (GC)- and liquid chromatography (LC)-mass spectrometry (MS). Fukushima et al. (2014) generated metabolite profiles from various characterized and uncharacterized mutant plants and clustered the mutants with similar metabolic phenotypes by conducting multidimensional scaling with quantified metabolic phenotypes. Nonetheless, representation and analysis of such a large amount of data remains a challenge for scientific discovery (Lu et al., 2011). In addition, these studies do not examine the topological and functional characteristics of metabolic changes in the context of a genome-scale metabolic network. To understand the relationship between genotype and metabolic phenotype, we need to investigate the metabolic changes caused by perturbing the expression of a gene in a genome-scale metabolic network perspective, because metabolic pathways are not independent biochemical factories but are components of a complex network (Berg et al., 2002; Merico et al., 2009).Much progress has been made in the last 2 decades to represent metabolism at a genome scale (Terzer et al., 2009). The advances in genome sequencing and emerging fields such as biocuration and bioinformatics enabled the representation of genome-scale metabolic network reconstructions for model organisms (Bassel et al., 2012). Genome-scale metabolic models have been built and applied broadly from microbes to plants. The first step toward modeling a genome-scale metabolism in a plant species started with developing a genome-scale metabolic pathway database for Arabidopsis (AraCyc; Mueller et al., 2003) from reference pathway databases (Kanehisa and Goto, 2000; Karp et al., 2002; Zhang et al., 2010). Genome-scale metabolic pathway databases have been built for several plant species (Mueller et al., 2005; Zhang et al., 2005, 2010; Urbanczyk-Wochniak and Sumner, 2007; May et al., 2009; Dharmawardhana et al., 2013; Monaco et al., 2013, 2014; Van Moerkercke et al., 2013; Chae et al., 2014; Jung et al., 2014). Efforts have been made to develop predictive genome-scale metabolic models using enzyme kinetics and stoichiometric flux-balance approaches (Sweetlove et al., 2008). de Oliveira Dal’Molin et al. (2010) developed a genome-scale metabolic model for Arabidopsis and successfully validated the model by predicting the classical photorespiratory cycle as well as known key differences between redox metabolism in photosynthetic and nonphotosynthetic plant cells. Other genome-scale models have been developed for Arabidopsis (Poolman et al., 2009; Radrich et al., 2010; Mintz-Oron et al., 2012), C. reinhardtii (Chang et al., 2011; Dal’Molin et al., 2011), maize (Dal’Molin et al., 2010; Saha et al., 2011), sorghum (Sorghum bicolor; Dal’Molin et al., 2010), and sugarcane (Saccharum officinarum; Dal’Molin et al., 2010). These predictive models have the potential to be applied broadly in fields such as metabolic engineering, drug target discovery, identification of gene function, study of evolutionary processes, risk assessment of genetically modified crops, and interpretations of mutant phenotypes (Feist and Palsson, 2008; Ricroch et al., 2011).Here, we interrogate the metabotypes caused by 136 single gene perturbations of Arabidopsis by analyzing the relative concentration changes of 1,348 chemically identified metabolites using a reconstructed genome-scale metabolic network. We examine the characteristics of the changed metabolites (the metabolites whose relative concentrations were significantly different in mutants relative to the wild type) in the metabolic network to uncover biological and topological consequences of the perturbed genes.     

Nontransgenic Genome Modification in Plant Cells

Zinc finger nucleases (ZFNs) are a powerful tool for genome editing in eukaryotic cells. ZFNs have been used for targeted mutagenesis in model and crop species. In animal and human cells, transient ZFN expression is often achieved by direct gene transfer into the target cells. Stable transformation, however, is the preferred method for gene expression in plant species, and ZFN-expressing transgenic plants have been used for recovery of mutants that are likely to be classified as transgenic due to the use of direct gene-transfer methods into the target cells. Here we present an alternative, nontransgenic approach for ZFN delivery and production of mutant plants using a novel Tobacco rattle virus (TRV)-based expression system for indirect transient delivery of ZFNs into a variety of tissues and cells of intact plants. TRV systemically infected its hosts and virus ZFN-mediated targeted mutagenesis could be clearly observed in newly developed infected tissues as measured by activation of a mutated reporter transgene in tobacco (Nicotiana tabacum) and petunia (Petunia hybrida) plants. The ability of TRV to move to developing buds and regenerating tissues enabled recovery of mutated tobacco and petunia plants. Sequence analysis and transmission of the mutations to the next generation confirmed the stability of the ZFN-induced genetic changes. Because TRV is an RNA virus that can infect a wide range of plant species, it provides a viable alternative to the production of ZFN-mediated mutants while avoiding the use of direct plant-transformation methods.Methods for genome editing in plant cells have fallen behind the remarkable progress made in whole-genome sequencing projects. The availability of reliable and efficient methods for genome editing would foster gene discovery and functional gene analyses in model plants and the introduction of novel traits in agriculturally important species (Puchta, 2002; Hanin and Paszkowski, 2003; Reiss, 2003; Porteus, 2009). Genome editing in various species is typically achieved by integrating foreign DNA molecules into the target genome by homologous recombination (HR). Genome editing by HR is routine in yeast (Saccharomyces cerevisiae) cells (Scherer and Davis, 1979) and has been adapted for other species, including Drosophila, human cell lines, various fungal species, and mouse embryonic stem cells (Baribault and Kemler, 1989; Venken and Bellen, 2005; Porteus, 2007; Hall et al., 2009; Laible and Alonso-González, 2009; Tenzen et al., 2009). In plants, however, foreign DNA molecules, which are typically delivered by direct gene-transfer methods (e.g. Agrobacterium and microbombardment of plasmid DNA), often integrate into the target cell genome via nonhomologous end joining (NHEJ) and not HR (Ray and Langer, 2002; Britt and May, 2003).Various methods have been developed to indentify and select for rare site-specific foreign DNA integration events or to enhance the rate of HR-mediated DNA integration in plant cells. Novel T-DNA molecules designed to support strong positive- and negative-selection schemes (e.g. Thykjaer et al., 1997; Terada et al., 2002), altering the plant DNA-repair machinery by expressing yeast chromatin remodeling protein (Shaked et al., 2005), and PCR screening of large numbers of transgenic plants (Kempin et al., 1997; Hanin et al., 2001) are just a few of the experimental approaches used to achieve HR-mediated gene targeting in plant species. While successful, these approaches, and others, have resulted in only a limited number of reports describing the successful implementation of HR-mediated gene targeting of native and transgenic sequences in plant cells (for review, see Puchta, 2002; Hanin and Paszkowski, 2003; Reiss, 2003; Porteus, 2009; Weinthal et al., 2010).HR-mediated gene targeting can potentially be enhanced by the induction of genomic double-strand breaks (DSBs). In their pioneering studies, Puchta et al. (1993, 1996) showed that DSB induction by the naturally occurring rare-cutting restriction enzyme I-SceI leads to enhanced HR-mediated DNA repair in plants. Expression of I-SceI and another rare-cutting restriction enzyme (I-CeuI) also led to efficient NHEJ-mediated site-specific mutagenesis and integration of foreign DNA molecules in plants (Salomon and Puchta, 1998; Chilton and Que, 2003; Tzfira et al., 2003). Naturally occurring rare-cutting restriction enzymes thus hold great promise as a tool for genome editing in plant cells (Carroll, 2004; Pâques and Duchateau, 2007). However, their wide application is hindered by the tedious and next to impossible reengineering of such enzymes for novel DNA-target specificities (Pâques and Duchateau, 2007).A viable alternative to the use of rare-cutting restriction enzymes is the zinc finger nucleases (ZFNs), which have been used for genome editing in a wide range of eukaryotic species, including plants (e.g. Bibikova et al., 2001; Porteus and Baltimore, 2003; Lloyd et al., 2005; Urnov et al., 2005; Wright et al., 2005; Beumer et al., 2006; Moehle et al., 2007; Santiago et al., 2008; Shukla et al., 2009; Tovkach et al., 2009; Townsend et al., 2009; Osakabe et al., 2010; Petolino et al., 2010; Zhang et al., 2010). Here too, ZFNs have been used to enhance DNA integration via HR (e.g. Shukla et al., 2009; Townsend et al., 2009) and as an efficient tool for the induction of site-specific mutagenesis (e.g. Lloyd et al., 2005; Zhang et al., 2010) in plant species. The latter is more efficient and simpler to implement in plants as it does not require codelivery of both ZFN-expressing and donor DNA molecules and it relies on NHEJ—the dominant DNA-repair machinery in most plant species (Ray and Langer, 2002; Britt and May, 2003).ZFNs are artificial restriction enzymes composed of a fusion between an artificial Cys₂His₂ zinc-finger protein DNA-binding domain and the cleavage domain of the FokI endonuclease. The DNA-binding domain of ZFNs can be engineered to recognize a variety of DNA sequences (for review, see Durai et al., 2005; Porteus and Carroll, 2005; Carroll et al., 2006). The FokI endonuclease domain functions as a dimer, and digestion of the target DNA requires proper alignment of two ZFN monomers at the target site (Durai et al., 2005; Porteus and Carroll, 2005; Carroll et al., 2006). Efficient and coordinated expression of both monomers is thus required for the production of DSBs in living cells. Transient ZFN expression, by direct gene delivery, is the method of choice for targeted mutagenesis in human and animal cells (e.g. Urnov et al., 2005; Beumer et al., 2006; Meng et al., 2008). Among the different methods used for high and efficient transient ZFN delivery in animal and human cell lines are plasmid injection (Morton et al., 2006; Foley et al., 2009), direct plasmid transfer (Urnov et al., 2005), the use of integrase-defective lentiviral vectors (Lombardo et al., 2007), and mRNA injection (Takasu et al., 2010).In plant species, however, efficient and strong gene expression is often achieved by stable gene transformation. Both transient and stable ZFN expression have been used in gene-targeting experiments in plants (Lloyd et al., 2005; Wright et al., 2005; Maeder et al., 2008; Cai et al., 2009; de Pater et al., 2009; Shukla et al., 2009; Tovkach et al., 2009; Townsend et al., 2009; Osakabe et al., 2010; Petolino et al., 2010; Zhang et al., 2010). In all cases, direct gene-transformation methods, using polyethylene glycol, silicon carbide whiskers, or Agrobacterium, were deployed. Thus, while mutant plants and tissues could be recovered, potentially without any detectable traces of foreign DNA, such plants were generated using a transgenic approach and are therefore still likely to be classified as transgenic. Furthermore, the recovery of mutants in many cases is also dependent on the ability to regenerate plants from protoplasts, a procedure that has only been successfully applied in a limited number of plant species. Therefore, while ZFN technology is a powerful tool for site-specific mutagenesis, its wider implementation for plant improvement may be somewhat limited, both by its restriction to certain plant species and by legislative restrictions imposed on transgenic plants.Here we describe an alternative to direct gene transfer for ZFN delivery and for the production of mutated plants. Our approach is based on the use of a novel Tobacco rattle virus (TRV)-based expression system, which is capable of systemically infecting its host and spreading into a variety of tissues and cells of intact plants, including developing buds and regenerating tissues. We traced the indirect ZFN delivery in infected plants by activation of a mutated reporter gene and we demonstrate that this approach can be used to recover mutated plants.     

Multiscale Metabolic Modeling: Dynamic Flux Balance Analysis on a Whole-Plant Scale

Plant metabolism is characterized by a unique complexity on the cellular, tissue, and organ levels. On a whole-plant scale, changing source and sink relations accompanying plant development add another level of complexity to metabolism. With the aim of achieving a spatiotemporal resolution of source-sink interactions in crop plant metabolism, a multiscale metabolic modeling (MMM) approach was applied that integrates static organ-specific models with a whole-plant dynamic model. Allowing for a dynamic flux balance analysis on a whole-plant scale, the MMM approach was used to decipher the metabolic behavior of source and sink organs during the generative phase of the barley (Hordeum vulgare) plant. It reveals a sink-to-source shift of the barley stem caused by the senescence-related decrease in leaf source capacity, which is not sufficient to meet the nutrient requirements of sink organs such as the growing seed. The MMM platform represents a novel approach for the in silico analysis of metabolism on a whole-plant level, allowing for a systemic, spatiotemporally resolved understanding of metabolic processes involved in carbon partitioning, thus providing a novel tool for studying yield stability and crop improvement.Plants are of vital significance as a source of food (Grusak and DellaPenna, 1999; Rogalski and Carrer, 2011), feed (Lu et al., 2011), energy (Tilman et al., 2006; Parmar et al., 2011), and feedstocks for the chemical industry (Metzger and Bornscheuer, 2006; Kinghorn et al., 2011). Given the close connection between plant metabolism and the usability of plant products, there is a growing interest in understanding and predicting the behavior and regulation of plant metabolic processes. In order to increase crop quality and yield, there is a need for methods guiding the rational redesign of the plant metabolic network (Schwender, 2009).Mathematical modeling of plant metabolism offers new approaches to understand, predict, and modify complex plant metabolic processes. In plant research, the issue of metabolic modeling is constantly gaining attention, and different modeling approaches applied to plant metabolism exist, ranging from highly detailed quantitative to less complex qualitative approaches (for review, see Giersch, 2000; Morgan and Rhodes, 2002; Poolman et al., 2004; Rios-Estepa and Lange, 2007).A widely used modeling approach is flux balance analysis (FBA), which allows the prediction of metabolic capabilities and steady-state fluxes under different environmental and genetic backgrounds using (non)linear optimization (Orth et al., 2010). Assuming steady-state conditions, FBA has the advantage of not requiring the knowledge of kinetic parameters and, therefore, can be applied to model detailed, large-scale systems. In recent years, the FBA approach has been applied to several different plant species, such as maize (Zea mays; Dal’Molin et al., 2010; Saha et al., 2011), barley (Hordeum vulgare; Grafahrend-Belau et al., 2009b; Melkus et al., 2011; Rolletschek et al., 2011), rice (Oryza sativa; Lakshmanan et al., 2013), Arabidopsis (Arabidopsis thaliana; Poolman et al., 2009; de Oliveira Dal’Molin et al., 2010; Radrich et al., 2010; Williams et al., 2010; Mintz-Oron et al., 2012; Cheung et al., 2013), and rapeseed (Brassica napus; Hay and Schwender, 2011a, 2011b; Pilalis et al., 2011), as well as algae (Boyle and Morgan, 2009; Cogne et al., 2011; Dal’Molin et al., 2011) and photoautotrophic bacteria (Knoop et al., 2010; Montagud et al., 2010; Boyle and Morgan, 2011). These models have been used to study different aspects of metabolism, including the prediction of optimal metabolic yields and energy efficiencies (Dal’Molin et al., 2010; Boyle and Morgan, 2011), changes in flux under different environmental and genetic backgrounds (Grafahrend-Belau et al., 2009b; Dal’Molin et al., 2010; Melkus et al., 2011), and nonintuitive metabolic pathways that merit subsequent experimental investigations (Poolman et al., 2009; Knoop et al., 2010; Rolletschek et al., 2011). Although FBA of plant metabolic models was shown to be capable of reproducing experimentally determined flux distributions (Williams et al., 2010; Hay and Schwender, 2011b) and generating new insights into metabolic behavior, capacities, and efficiencies (Sweetlove and Ratcliffe, 2011), challenges remain to advance the utility and predictive power of the models.Given that many plant metabolic functions are based on interactions between different subcellular compartments, cell types, tissues, and organs, the reconstruction of organ-specific models and the integration of these models into interacting multiorgan and/or whole-plant models is a prerequisite to get insight into complex plant metabolic processes organized on a whole-plant scale (e.g. source-sink interactions). Almost all FBA models of plant metabolism are restricted to one cell type (Boyle and Morgan, 2009; Knoop et al., 2010; Montagud et al., 2010; Cogne et al., 2011; Dal’Molin et al., 2011), one tissue or one organ (Grafahrend-Belau et al., 2009b; Hay and Schwender, 2011a, 2011b; Pilalis et al., 2011; Mintz-Oron et al., 2012), and only one model exists taking into account the interaction between two cell types by specifying the interaction between mesophyll and bundle sheath cells in C4 photosynthesis (Dal’Molin et al., 2010). So far, no model representing metabolism at the whole-plant scale exists.Considering whole-plant metabolism raises the problem of taking into account temporal and environmental changes in metabolism during plant development and growth. Although classical static FBA is unable to predict the dynamics of metabolic processes, as the network analysis is based on steady-state solutions, time-dependent processes can be taken into account by extending the classical static FBA to a dynamic flux balance analysis (dFBA), as proposed by Mahadevan et al. (2002). The static (SOA) and dynamic optimization approaches introduced in this work provide a framework for analyzing the transience of metabolism by integrating kinetic expressions to dynamically constrain exchange fluxes. Due to the requirement of knowing or estimating a large number of kinetic parameters, so far dFBA has only been applied to a plant metabolic model once, to study the photosynthetic metabolism in the chloroplasts of C3 plants by a simplified model of five biochemical reactions (Luo et al., 2009). Integrating a dynamic model into a static FBA model is an alternative approach to perform dFBA.In this study, a multiscale metabolic modeling (MMM) approach was applied with the aim of achieving a spatiotemporal resolution of cereal crop plant metabolism. To provide a framework for the in silico analysis of the metabolic dynamics of barley on a whole-plant scale, the MMM approach integrates a static multiorgan FBA model and a dynamic whole-plant multiscale functional plant model (FPM) to perform dFBA. The performance of the novel whole-plant MMM approach was tested by studying source-sink interactions during the seed developmental phase of barley plants.     

Tyloses and Phenolic Deposits in Xylem Vessels Impede Water Transport in Low-Lignin Transgenic Poplars: A Study by Cryo-Fluorescence Microscopy

Of 14 transgenic poplar genotypes (Populus tremula × Populus alba) with antisense 4-coumarate:coenzyme A ligase that were grown in the field for 2 years, five that had substantial lignin reductions also had greatly reduced xylem-specific conductivity compared with that of control trees and those transgenic events with small reductions in lignin. For the two events with the lowest xylem lignin contents (greater than 40% reduction), we used light microscopy methods and acid fuchsin dye ascent studies to clarify what caused their reduced transport efficiency. A novel protocol involving dye stabilization and cryo-fluorescence microscopy enabled us to visualize the dye at the cellular level and to identify water-conducting pathways in the xylem. Cryo-fixed branch segments were planed in the frozen state on a sliding cryo-microtome and observed with an epifluorescence microscope equipped with a cryo-stage. We could then distinguish clearly between phenolic-occluded vessels, conductive (stain-filled) vessels, and nonconductive (water- or gas-filled) vessels. Low-lignin trees contained areas of nonconductive, brown xylem with patches of collapsed cells and patches of noncollapsed cells filled with phenolics. In contrast, phenolics and nonconductive vessels were rarely observed in normal colored wood of the low-lignin events. The results of cryo-fluorescence light microscopy were supported by observations with a confocal microscope after freeze drying of cryo-planed samples. Moreover, after extraction of the phenolics, confocal microscopy revealed that many of the vessels in the nonconductive xylem were blocked with tyloses. We conclude that reduced transport efficiency of the transgenic low-lignin xylem was largely caused by blockages from tyloses and phenolic deposits within vessels rather than by xylem collapse.Secondary xylem in woody plants is a complex vascular tissue that functions in mechanical support, conduction, storage, and protection (Carlquist, 2001; Tyree and Zimmermann, 2002). The xylem must provide a sufficient and safe water supply throughout the entire pathway from roots to leaves for transpiration and photosynthesis. It is well established that enhanced water conductivity of xylem can increase total plant carbon gain (Domec and Gartner, 2003; Santiago et al., 2004; Brodribb and Holbrook, 2005a). According to the Hagen-Poiseuille equation, xylem conductivity should scale with vessel lumen diameter to the fourth power (Tyree and Zimmermann, 2002). Indeed, xylem conductivity largely depends on anatomical features, including conduit diameters and frequencies (Salleo et al., 1985; McCulloh and Sperry, 2005). However, there are hydraulic limits to maximum vessel diameters, because xylem conduits have to withstand the strong negative pressures of the transpiration stream that could cause cell collapse or embolisms within vessels that are structurally inadequate to withstand these forces (Tyree and Sperry, 1989; Lo Gullo et al., 1995; Hacke et al., 2000). To some extent, stomatal regulation of transpiration limits the negative pressures that the xylem experiences (Tardieu and Davies, 1993; Cochard et al., 2002; Meinzer, 2002; Brodribb and Holbrook, 2004; Buckley, 2005; Franks et al., 2007; Woodruff et al., 2007). Nevertheless, plants rely on an array of structural reinforcements of xylem to ensure the safety of water transport. The size of xylem elements, vessel redundancy, intervessel pit and membrane geometries, and the thickness, microstructure, and chemical composition of cell walls are among the features that regulate tradeoffs between efficiency and safety of xylem water transport (Baas and Schweingruber, 1987; Hacke et al., 2001; Domec et al., 2006; Ewers et al., 2007; Choat et al., 2008).The xylem cell wall is made up of cellulose bundles that are hydrogen bonded with hemicelluloses, which are in turn embedded within a lignin matrix (Mansfield, 2009; Salmén and Burgert, 2009). Besides providing this matrix for the cell wall itself, lignin is thought to contribute to many of the mechanical and physical characteristics of wood as well as conferring passive resistance to the spread of pathogens within a plant (Niklas, 1992; Boyce et al., 2004; Davin et al., 2008). Lignin typically represents 20% to 30% of the dry mass of wood and therefore is among the most abundant stores of carbon in the biosphere (Zobel and van Buijtenen, 1989). The complex molecular structure and biosynthetic pathway of various types of lignins have been studied extensively (Boerjan et al., 2003; Ralph et al., 2004, 2007; Higuchi, 2006; Boudet, 2007; Davin et al., 2008). The monomeric composition of lignin varies between different cell types of the same species depending on the functional specialization of the cell (Yoshinaga et al., 1992; Watanabe et al., 2004; Xu et al., 2006). The composition and amount of lignin in wild plants varies in response to climatic conditions (Donaldson, 2002) or gravitational and mechanical demands (Pruyn et al., 2000; Kern et al., 2005; Rüggeberg et al., 2008). It is clear that plants are capable of regulating the lignification pattern in differentiating cells, which provides them with flexibility for responding to environmental stresses (Donaldson, 2002; Koehler and Telewski, 2006; Ralph et al., 2007; for review, see Vanholme et al., 2008).Whereas some level of lignin is a requisite for all vascular plants, it is often an unwanted product in the pulp and paper industry because it increases the costs of paper production and associated water treatments necessary for environmental protection (Chen et al., 2001; Baucher et al., 2003; Peter et al., 2007). Reducing the lignin content of the raw biomass material may allow more efficient hydrolysis of polysaccharides in biomass and thus facilitate the production of biofuel (Chen and Dixon, 2007). With the ultimate goal of development of wood for more efficient processing, much research has been aimed at the production of genetically modified trees with altered lignin biosynthesis (Boerjan et al., 2003; Boudet et al., 2003; Li et al., 2003; Halpin, 2004; Ralph et al., 2004, 2008; Chiang, 2006; Coleman et al., 2008a, 2008b; Vanholme et al., 2008; Wagner et al., 2009). It is now technically possible to achieve more than 50% reductions of lignin content in xylem of poplar (Populus spp.; Leplé et al., 2007; Coleman et al., 2008a, 2008b), but the consequences of such reduction on plant function have received relatively little attention (Koehler and Telewski, 2006). In-depth studies on the xylem structure and functional performance of transgenic plants with low lignin are limited, despite the need to assess their long-term sustainability for large-scale production (Anterola and Lewis, 2002; Hancock et al., 2007; Coleman et al., 2008b, Voelker, 2009; Horvath et al., 2010).Genetically modified plants are suitable models for studying fundamental questions of the physiological role of lignin because of the possibility of controlling lignification without the confounding effects encountered when comparing across plant tissues or stages of development (Koehler and Telewski, 2006; Leplé et al., 2007; Coleman et al., 2008b). Research on Arabidopsis (Arabidopsis thaliana) and tobacco (Nicotiana tabacum) has shown that down-regulation of lignin biosynthesis can have diverse effects on plant metabolism and structure, including changes in the lignin amount and composition (p-hydroxyphenyl/guaiacyl/syringyl units ratio) as well as the collapse of xylem vessel elements (Lee et al., 1997; Sewalt et al., 1997; Piquemal et al., 1998; Chabannes et al., 2001; Jones et al., 2001; Franke et al., 2002; Dauwe et al., 2007). Among temperate hardwoods, poplar has been established as a model tree for genetic manipulations because of its ecological and economic importance, fast growth, ease of vegetative propagation, and its widespread use in traditional breeding programs (Bradshaw et al., 2001; Brunner et al., 2004). The question of how manipulation of lignin can affect the anatomy and physiological function of xylem in poplar has been addressed in part by several research groups (Anterola and Lewis, 2002; Boerjan et al., 2003; Leplé et al., 2007; Coleman et al., 2008b). Some studies that involved large lignin reductions reported no significant alterations in the xylem anatomy (Hu et al., 1999; Li et al., 2003). However, in many other experiments, reduced total lignin content was associated with significant growth retardation, alterations in the lignin monomer composition, irregularities in the xylem structure (Anterola and Lewis, 2002; Leplé et al., 2007; Coleman et al., 2008b), and the patchy occurrence of collapsed xylem cells (Coleman et al., 2008b; Voelker, 2009). Furthermore, severely down-regulated lignin biosynthesis has resulted in greatly reduced xylem water-transport efficiency (Coleman et al., 2008b; Lachenbruch et al., 2009; Voelker, 2009). It is generally assumed that the reduced water transport ability of xylem with very low lignin contents is caused by collapsed conduits and/or increased embolism due to the entry of air bubbles into the water-conducting cells (Coleman et al., 2008b; Wagner et al., 2009), but detailed anatomical investigations of the causes of impaired xylem conductivity of low-lignin trees are lacking. Analysis of the anatomical basis for the properties of xylem conduits in plants with genetically manipulated amounts and composition of lignin can provide a deeper understanding of the physiological role of lignin as well as the lower limit of down-regulation of lignin biosynthesis at which trees can still survive within natural environments.One of the approaches for the suppression of lignin biosynthesis is down-regulation of 4-coumarate:coenzyme A ligase (4CL), an enzyme that functions in phenylpropanoid metabolism by producing the monolignol precursor p-coumaroyl-CoA (Kajita et al.,1997; Allina et al., 1998; Hu et al., 1998; Harding et al., 2002; Jia et al., 2004; Costa et al., 2005; Friedmann et al., 2007; Wagner et al., 2009). In a 2-year-long field trial on the physiological performance of poplar (Populus tremula × Populus alba) transgenic clones, out of 14 genotypes with altered lignin biosynthesis (down-regulated 4CL), five showed dramatically reduced wood-specific conductivity (k_s) compared with that of control trees (Voelker, 2009). Those mutants with the severely reduced k_s were also characterized by having the lowest wood lignin contents (up to an approximately 40% reduction) in the study. Trees with transgenic events characterized by the formation of abnormally brown wood exhibited regular branch dieback at the end of the growing season, despite having been regularly watered (Voelker, 2009). Our objective was to identify the structural features responsible for reduced transport efficiency in the xylem of transgenic poplars with extremely low lignin contents. We employed fluorescence and laser scanning confocal microscopy for anatomical analyses of xylem structure as well as dye-flow experiments followed by cryo-fluorescence microscopy to visualize the functioning water-conductive pathways in xylem at the cellular level. We report the frequent occurrence of tyloses and phenolic depositions in xylem vessels of strongly down-regulated trees that may be the cause of their reduced xylem conductivity.     

Wettability,Polarity, and Water Absorption of Holm Oak Leaves: Effect of Leaf Side and Age

Plant trichomes play important protective functions and may have a major influence on leaf surface wettability. With the aim of gaining insight into trichome structure, composition, and function in relation to water-plant surface interactions, we analyzed the adaxial and abaxial leaf surface of holm oak (Quercus ilex) as a model. By measuring the leaf water potential 24 h after the deposition of water drops onto abaxial and adaxial surfaces, evidence for water penetration through the upper leaf side was gained in young and mature leaves. The structure and chemical composition of the abaxial (always present) and adaxial (occurring only in young leaves) trichomes were analyzed by various microscopic and analytical procedures. The adaxial surfaces were wettable and had a high degree of water drop adhesion in contrast to the highly unwettable and water-repellent abaxial holm oak leaf sides. The surface free energy and solubility parameter decreased with leaf age, with higher values determined for the adaxial sides. All holm oak leaf trichomes were covered with a cuticle. The abaxial trichomes were composed of 8% soluble waxes, 49% cutin, and 43% polysaccharides. For the adaxial side, it is concluded that trichomes and the scars after trichome shedding contribute to water uptake, while the abaxial leaf side is highly hydrophobic due to its high degree of pubescence and different trichome structure, composition, and density. Results are interpreted in terms of water-plant surface interactions, plant surface physical chemistry, and plant ecophysiology.Plant surfaces have an important protecting function against multiple biotic and abiotic stress factors (Riederer, 2006). They may, for example, limit the attack of insects (Eigenbrode and Jetter, 2002) or pathogenic fungi (Gniwotta et al., 2005; Łaźniewska et al., 2012), avoid damage caused by high intensities of UV and visible radiation (Reicosky and Hanover, 1978; Karabourniotis and Bormann, 1999), help to regulate leaf temperature (Ehleringer and Björkman, 1978; Ripley et al., 1999), and chiefly prevent plant organs from dehydration (Riederer and Schreiber, 2001).The epidermis of plants has been found to have a major degree of physical and chemical variability and may often contain specialized cells such as trichomes or stomata (Roth-Nebelsick et al., 2009; Javelle et al., 2011). Most aerial organs are covered with an extracellular and generally lipid-rich layer named the cuticle, which is typically composed of waxes embedded in (intracuticular waxes) or deposited on (epicuticular waxes) a biopolymer matrix of cutin, forming a network of cross-esterified hydroxy C₁₆ and/or C₁₈ fatty acids, and/or cutan, with variable amounts of polysaccharides and phenolics (Domínguez et al., 2011; Yeats and Rose, 2013). Different nano- and/or microscale levels of plant surface sculpturing have been observed by scanning electron microscopy (SEM), generally in relation to the topography of epicuticular waxes, cuticular folds, and epidermal cells (Koch and Barthlott, 2009). Such surface features together with their chemical composition (Khayet and Fernández, 2012) may lead to a high degree of roughness and hydrophobicity (Koch and Barthlott, 2009; Konrad et al., 2012). The interactions of plant surfaces with water have been addressed in some investigations (Brewer et al., 1991; Brewer and Smith, 1997; Pandey and Nagar, 2003; Hanba et al., 2004; Dietz et al., 2007; Holder, 2007a, 2007b; Fernández et al., 2011, 2014; Roth-Nebelsick et al., 2012; Wen et al., 2012; Urrego-Pereira et al., 2013) and are a topic of growing interest for plant ecophysiology (Helliker and Griffiths, 2007; Aryal and Neuner, 2010; Limm and Dawson, 2010; Kim and Lee, 2011; Berry and Smith, 2012; Berry et al., 2013; Rosado and Holder, 2013; Helliker, 2014). On the other hand, the mechanisms of foliar uptake of water and solutes by plant surfaces are still not fully understood (Fernández and Eichert, 2009; Burkhardt and Hunsche, 2013), but they may play an important ecophysiological role (Limm et al., 2009; Johnstone and Dawson, 2010; Adamec, 2013; Berry et al., 2014).The importance of trichomes and pubescent layers on water drop-plant surface interactions and on the subsequent potential water uptake into the organs has been analyzed in some investigations (Fahn, 1986; Brewer et al., 1991; Grammatikopoulos and Manetas, 1994; Brewer and Smith, 1997; Pierce et al., 2001; Kenzo et al., 2008; Fernández et al., 2011, 2014; Burrows et al., 2013). Trichomes are unicellular or multicellular and glandular or nonglandular appendages, which originate from epidermal cells only and develop outwards on the surface of plant organs (Werker, 2000). Nonglandular trichomes are categorized according to their morphology and exhibit a major variability in size, morphology, and function. On the other hand, glandular trichomes are classified by the secretory materials they excrete, accumulate, or absorb (Johnson, 1975; Werker, 2000; Wagner et al., 2004). Trichomes can be often found in xeromorphic leaves and in young organs (Fahn, 1986; Karabourniotis et al., 1995). The occurrence of protecting leaf trichomes has been also reported for Mediterranean species such as holm oak (Quercus ilex; Karabourniotis et al., 1995, 1998; Morales et al., 2002; Karioti et al., 2011; Camarero et al., 2012). There is limited information about the nature of the surface of trichomes, but they are also covered with a cuticle similarly to other epidermal cell types (Fernández et al., 2011, 2014).In this study and using holm oak as a model, we assessed, for the first time, the leaf surface-water relations of the abaxial (always pubescent) versus the adaxial (only pubescent in developing leaves and for a few months) surface, including their capacity to absorb surface-deposited water drops. Based on membrane science methodologies (Fernández et al., 2011; Khayet and Fernández, 2012) and following a new integrative approach, the chemical, physical, and anatomical properties of holm oak leaf surfaces and trichomes were analyzed, with the aim of addressing the following questions. Are young and mature adaxial and abaxial leaf surfaces capable of absorbing water deposited as drops on to the surfaces? Are young and mature abaxial and adaxial leaf surfaces similar in relation to their wettability, hydrophobicity, polarity, work of adhesion (W_a) for water, solubility parameter (δ), and surface free energy (γ)? What is the physical and chemical nature of the adaxial versus the abaxial trichomes, chiefly in relation to young leaves?     

CELLULOSE SYNTHASE-LIKE A2, a Glucomannan Synthase,Is Involved in Maintaining Adherent Mucilage Structure in Arabidopsis Seed

Mannans are hemicellulosic polysaccharides that are considered to have both structural and storage functions in the plant cell wall. However, it is not yet known how mannans function in Arabidopsis (Arabidopsis thaliana) seed mucilage. In this study, CELLULOSE SYNTHASE-LIKE A2 (CSLA2; At5g22740) expression was observed in several seed tissues, including the epidermal cells of developing seed coats. Disruption of CSLA2 resulted in thinner adherent mucilage halos, although the total amount of the adherent mucilage did not change compared with the wild type. This suggested that the adherent mucilage in the mutant was more compact compared with that of the wild type. In accordance with the role of CSLA2 in glucomannan synthesis, csla2-1 mucilage contained 30% less mannosyl and glucosyl content than did the wild type. No appreciable changes in the composition, structure, or macromolecular properties were observed for nonmannan polysaccharides in mutant mucilage. Biochemical analysis revealed that cellulose crystallinity was substantially reduced in csla2-1 mucilage; this was supported by the removal of most mucilage cellulose through treatment of csla2-1 seeds with endo-β-glucanase. Mutation in CSLA2 also resulted in altered spatial distribution of cellulose and an absence of birefringent cellulose microfibrils within the adherent mucilage. As with the observed changes in crystalline cellulose, the spatial distribution of pectin was also modified in csla2-1 mucilage. Taken together, our results demonstrate that glucomannans synthesized by CSLA2 are involved in modulating the structure of adherent mucilage, potentially through altering cellulose organization and crystallization.Mannan polysaccharides are a complex set of hemicellulosic cell wall polymers that are considered to have both structural and storage functions. Based on the particular chemical composition of the backbone and the side chains, mannan polysaccharides are classified into four types: pure mannan, glucomannan, galactomannan, and galactoglucomannan (Moreira and Filho, 2008; Wang et al., 2012; Pauly et al., 2013). Each of these polysaccharides is composed of a β-1,4-linked backbone containing Man or a combination of Glc and Man residues. In addition, the mannan backbone can be substituted with side chains of α-1,6-linked Gal residues. Mannan polysaccharides have been proposed to cross link with cellulose and other hemicelluloses via hydrogen bonds (Fry, 1986; Iiyama et al., 1994; Obel et al., 2007; Scheller and Ulvskov, 2010). Furthermore, it has been reported that heteromannans with different levels of substitution can interact with cellulose in diverse ways (Whitney et al., 1998). Together, these observations indicate the complexity of mannan polysaccharides in the context of cell wall architecture.CELLULOSE SYNTHASE-LIKE A (CSLA) enzymes have been shown to have mannan synthase activity in vitro. These enzymes polymerize the β-1,4-linked backbone of mannans or glucomannans, depending on the substrates (GDP-Man and/or GDP-Glc) provided (Richmond and Somerville, 2000; Liepman et al., 2005, 2007; Pauly et al., 2013). In Arabidopsis (Arabidopsis thaliana), nine CSLA genes have been identified; different CSLAs are responsible for the synthesis of different mannan types (Liepman et al., 2005, 2007). CSLA7 has mannan synthase activity in vitro (Liepman et al., 2005) and has been shown to synthesize stem glucomannan in vivo (Goubet et al., 2009). Disrupting the CSLA7 gene results in defective pollen growth and embryo lethality phenotypes in Arabidopsis, indicating structural or signaling functions of mannan polysaccharides during plant embryo development (Goubet et al., 2003). A mutation in CSLA9 results in the inhibition of Agrobacterium tumefaciens-mediated root transformation in the rat4 mutant (Zhu et al., 2003). CSLA2, CSLA3, and CSLA9 are proposed to play nonredundant roles in the biosynthesis of stem glucomannans, although mutations in CSLA2, CSLA3, or CSLA9 have no effect on stem development or strength (Goubet et al., 2009). All of the Arabidopsis CSLA proteins have been shown to be involved in the biosynthesis of mannan polysaccharides in the plant cell wall (Liepman et al., 2005, 2007), although the precise physiological functions of only CSLA7 and CSLA9 have been conclusively demonstrated.In Arabidopsis, when mature dry seeds are hydrated, gel-like mucilage is extruded to envelop the entire seed. Ruthenium red staining of Arabidopsis seeds reveals two different mucilage layers, termed the nonadherent and the adherent mucilage layers (Western et al., 2000; Macquet et al., 2007a). The outer, nonadherent mucilage is loosely attached and can be easily extracted by shaking seeds in water. Compositional and linkage analyses suggest that this layer is almost exclusively composed of unbranched rhamnogalacturonan I (RG-I) (>80% to 90%), with small amounts of branched RG-I, arabinoxylan, and high methylesterified homogalacturonan (HG). By contrast, the inner, adherent mucilage layer is tightly attached to the seed and can only be removed by strong acid or base treatment, or by enzymatic digestion (Macquet et al., 2007a; Huang et al., 2011; Walker et al., 2011). As with the nonadherent layer, adherent mucilage is also mainly composed of unbranched RG-I, but with small numbers of arabinan and galactan ramifications (Penfield et al., 2001; Willats et al., 2001; Dean et al., 2007; Macquet et al., 2007a, 2007b; Arsovski et al., 2009; Haughn and Western, 2012). There are also minor amounts of pectic HG in the adherent mucilage, with high methylesterified HG in the external domain compared with the internal domain of the adherent layer (Willats et al., 2001; Macquet et al., 2007a; Rautengarten et al., 2008; Sullivan et al., 2011; Saez-Aguayo et al., 2013). In addition, the adherent mucilage contains cellulose (Blake et al., 2006; Macquet et al., 2007a), which is entangled with RG-I and is thought to anchor the pectin-rich mucilage onto seeds (Macquet et al., 2007a; Harpaz-Saad et al., 2011, 2012; Mendu et al., 2011; Sullivan et al., 2011). As such, Arabidopsis seed mucilage is considered to be a useful model for investigating the biosynthesis of cell wall polysaccharides and how this process is regulated in vivo (Haughn and Western, 2012).Screening for altered seed coat mucilage has led to the identification of several genes encoding enzymes that are involved in the biosynthesis or modification of mucilage components. RHAMNOSE SYNTHASE2/MUCILAGE-MODIFIED4 (MUM4) is responsible for the synthesis of UDP-l-Rha (Usadel et al., 2004; Western et al., 2004; Oka et al., 2007). The putative GALACTURONSYLTRANSFERASE11 can potentially synthesize mucilage RG-I or HG pectin from UDP-d-GalUA (Caffall et al., 2009). GALACTURONSYLTRANSFERASE-LIKE5 appears to function in the regulation of the final size of the mucilage RG-I (Kong et al., 2011, 2013). Mutant seeds defective in these genes display reduced thickness of the extruded mucilage layer compared with wild-type Arabidopsis seeds.RG-I deposited in the apoplast of seed coat epidermal cells appears to be synthesized in a branched form that is subsequently modified by enzymes in the apoplast. MUM2 encodes a β-galactosidase that removes Gal residues from RG-I side chains (Dean et al., 2007; Macquet et al., 2007b). β-XYLOSIDASE1 encodes an α-l-arabinfuranosidase that removes Ara residues from RG-I side chains (Arsovski et al., 2009). Disruptions of these genes lead to defective hydration properties and affect the extrusion of mucilage. Furthermore, correct methylesterification of mucilage HG is also required for mucilage extrusion. HG is secreted into the wall in a high methylesterified form that can then be enzymatically demethylesterified by pectin methylesterases (PMEs; Bosch and Hepler, 2005). PECTIN METHYLESTERASE INHIBITOR6 (PMEI6) inhibits PME activities (Saez-Aguayo et al., 2013). The subtilisin-like Ser protease (SBT1.7) can activate other PME inhibitors, but not PMEI6 (Rautengarten et al., 2008; Saez-Aguayo et al., 2013). Disruption of either PMEI6 or SBT1.7 results in the delay of mucilage release.Although cellulose is present at low levels in adherent mucilage, it plays an important adhesive role for the attachment of mucilage pectin to the seed coat epidermal cells. The orientation and amount of pectin associated with the cellulose network is largely determined by cellulose conformation properties (Macquet et al., 2007a; Haughn and Western, 2012). Previous studies have demonstrated that CELLULOSE SYNTHASE A5 (CESA5) is required for the production of seed mucilage cellulose and the adherent mucilage in the cesa5 mutant can be easily extracted with water (Harpaz-Saad et al., 2011, 2012; Mendu et al., 2011; Sullivan et al., 2011).Despite all of these discoveries, large gaps remain in the current knowledge of the biosynthesis and functions of mucilage polysaccharides in seed coats. In this study, we show that CSLA2 is involved in the biosynthesis of mucilage glucomannan. Furthermore, we show that CSLA2 functions in the maintenance of the normal structure of the adherent mucilage layer through modifying the mucilage cellulose ultrastructure.     

Engineering Triterpene and Methylated Triterpene Production in Plants Provides Biochemical and Physiological Insights into Terpene Metabolism

Linear, branch-chained triterpenes, including squalene (C30), botryococcene (C30), and their methylated derivatives (C31–C37), generated by the green alga Botryococcus braunii race B have received significant attention because of their utility as chemical and biofuel feedstocks. However, the slow growth habit of B. braunii makes it impractical as a production system. In this study, we evaluated the potential of generating high levels of botryococcene in tobacco (Nicotiana tabacum) plants by diverting carbon flux from the cytosolic mevalonate pathway or the plastidic methylerythritol phosphate pathway by the targeted overexpression of an avian farnesyl diphosphate synthase along with two versions of botryococcene synthases. Up to 544 µg g⁻¹ fresh weight of botryococcene was achieved when this metabolism was directed to the chloroplasts, which is approximately 90 times greater than that accumulating in plants engineered for cytosolic production. To test if methylated triterpenes could be produced in tobacco, we also engineered triterpene methyltransferases (TMTs) from B. braunii into wild-type plants and transgenic lines selected for high-level triterpene accumulation. Up to 91% of the total triterpene contents could be converted to methylated forms (C31 and C32) by cotargeting the TMTs and triterpene biosynthesis to the chloroplasts, whereas only 4% to 14% of total triterpenes were methylated when this metabolism was directed to the cytoplasm. When the TMTs were overexpressed in the cytoplasm of wild-type plants, up to 72% of the total squalene was methylated, and total triterpene (C30+C31+C32) content was elevated 7-fold. Altogether, these results point to innate mechanisms controlling metabolite fluxes, including a homeostatic role for squalene.Terpenes and terpenoids represent a distinct class of natural products (Buckingham, 2003) that are derived from two universal five-carbon precursors: isopentenyl diphosphate (IPP) and dimethylallyl diphosphate (DMAPP). In eukaryotic fungi and animals, IPP and DMAPP are synthesized via the mevalonate (MVA) pathway, whereas in prokaryotes, they are synthesized via the methylerythritol phosphate (MEP) pathway. In higher plants, the pathways are present in separate compartments and are believed to operate independently. The MVA pathway in the cytoplasm is predominantly responsible for sesquiterpene (C15), triterpene (C30), and polyprenol (greater than C45) biosynthesis and associated with the endoplasmic reticulum (ER) system. The MEP pathway resides in plastids and is dedicated to monoterpenes (C10), diterpenes (C20), carotenoids (C40), and long-chain phytol biosynthesis. All these compounds are usually produced by plants for a variety of physiological (i.e. hormones, aliphatic membrane anchors, and maintaining membrane structure) and ecological (i.e. defense compounds and insect/animal attractants) roles (Kempinski et al., 2015). Terpenes are also important for various industrial applications, ranging from flavors and fragrances (Schwab et al., 2008) to medicines (Dewick, 2009; Niehaus et al., 2011; Shelar, 2011).The utility of terpenes as chemical and biofuel feedstocks has also received considerable attention recently. Isoprenoid-derived biofuels include farnesane (Renninger and McPhee, 2008; Rude and Schirmer, 2009), bisabolene (Peralta-Yahya et al., 2011), pinene dimers (Harvey et al., 2010), isopentenal (Withers et al., 2007), and botryococcene (Moldowan and Seifert, 1980; Hillen et al., 1982; Glikson et al., 1989; Mastalerz and Hower, 1996). The richness of branches within these hydrocarbon scaffolds correlate with their high-energy content, which enables them to serve as suitable alternatives to crude petroleum (Peralta-Yahya and Keasling, 2010). Indeed, some of them are already major contributors to current-day petroleum-based fuels. One of the best examples of this is the triterpene oil accumulating in the green alga Botryococcus braunii race B, which is considered a major progenitor to oil and coal shale deposits (Moldowan and Seifert, 1980). This alga has been well studied, and the major constituents of its prodigious hydrocarbon oil are a group of triterpenes including squalene (C30), organism-specific botryococcene (C30), methylated squalene (C31–C34), and methylated botryococcene (C31–C37; Metzger et al., 1988; Huang and Poulter, 1989; Okada et al., 1995), which can be readily converted into all classes of combustible fuels under hydrocracking conditions (Hillen et al., 1982).The unique biosynthetic mechanism for the triterpenes in B. braunii was recently described by Niehaus et al. (2011), and a series of novel squalene synthase-like genes were identified (Fig. 1). In short, squalene synthase-like enzyme, SSL-1, performs a head-to-head condensation of two farnesyl diphosphate (FPP) molecules into presqualene diphosphate, followed by a reductive rearrangement to yield squalene (C30) by the enzyme SSL-2, or is converted by SSL-3 to form botryococcene through a different reductive rearrangement (Niehaus et al., 2011). Methylated derivatives are the dominant triterpene species generated by B. braunii race B (Metzger, 1985; Metzger et al., 1988), and these derivatives are known to yield higher quality fuels due to their high energy content and the hydrocracking products derived by virtue of having more hydrocarbon branches. Triterpene methyltransferases (TMTs) that can methylate squalene and botryococcene have been successfully characterized by Niehaus et al. (2012). TRITERPENE METHYLTRANSFERASE1 (TMT-1) and TMT-2 prefer squalene C30 as their substrate for the production of monomethylated (C31) or dimethylated (C32) squalene, while TMT-3 prefers botryococcene as its substrate for the biosynthesis of monomethylated (C31) or dimethylated (C32) botryococcene (Fig. 1). These TMTs are believed to be insoluble enzymes; they exhibit large hydrophobic areas, and their activities were only observed in vitro using yeast microsomal preparations (no activity was observed when expressed in bacteria; Niehaus et al., 2012).Open in a separate windowFigure 1.Depiction of the catalytic roles of novel SSL and TMT enzymes in B. braunii race B and their putative contributions to the triterpene constituents (Niehaus et al., 2011; Niehaus et al., 2012). SSL-1 catalyzes the condensation of two farnesyl diphosphate (FPP) molecules to presqualene diphosphate (PSPP), which is converted to either squalene or botryococcene by SSL-2 or SSL-3, respectively. Squalene can also be synthesized directly from the condensation of two FPP molecules catalyzed by squalene synthase (SQS). TMT-1 and TMT-2 transfer the methyl donor group from S-adenosylmethionine (SAM) to squalene to form monomethylated and dimethylated squalene, whereas TMT-3 acts on botryococcene to form monomethylated and dimethylated botryococcene (Niehaus et al., 2012).Like the majority of identified methyltransferases, these TMTs utilize the methyl donor S-adenosyl methionine (SAM), which is ubiquitous in prokaryotes and eukaryotes (Scheer et al., 2011; Liscombe et al., 2012). In plants, SAM is one of the most abundant cofactors (Fontecave et al., 2004; Sauter et al., 2013) and is synthesized exclusively in the cytosol (Wallsgrove et al., 1983; Ravanel et al., 1998, 2004; Bouvier et al., 2006). While it is used predominantly as a methyl donor in the methylation reaction (Ravanel et al., 2004), it also serves as the primary precursor for the biosynthesis of ethylene (Wang et al., 2002b), polyamines (Kusano et al., 2008), and nicotianamine (Takahashi et al., 2003), which play a variety of important roles for plant growth and development (Huang et al., 2012; Sauter et al., 2013). The SAM present in organelles, like the chloroplast, appears to be imported from the cytosol by specific SAM/S-adenosylhomocysteine exchange transporters that reside on the envelope membranes of plastids (Ravanel et al., 2004; Bouvier et al., 2006). The imported SAM is involved in the biogenesis of Asp-derived amino acids (Curien et al., 1998; Jander and Joshi, 2009; Sauter et al., 2013) and serves as the methyl donor for the methylation of macromolecules, such as plastid DNA (Nishiyama et al., 2002; Ahlert et al., 2009) and proteins (Houtz et al., 1989; Niemi et al., 1990; Ying et al., 1999; Trievel et al., 2003; Alban et al., 2014), and small molecule metabolites, such as prenylipids (e.g. plastoquinone, tocopherol, chlorophylls, and phylloquinone; Bouvier et al., 2005, 2006; DellaPenna, 2005).Although plants and microbes are the natural sources for useful terpenes, most of them are produced in very small amounts and often as complex mixtures. In contrast, B. braunii produces large quantities of triterpenes, but its slow growth makes it undesirable as a viable production platform (Niehaus et al., 2011). Nevertheless, metabolic engineering and synthetic biology offer many strategies to manipulate terpene metabolism in various biological systems to achieve high-value terpene production with high yield and high fidelity for particular practical applications (Nielsen and Keasling, 2011). Many successes have been achieved in engineering valuable terpenes in heterotrophic microbes, such as Escherichia coli (Nishiyama et al., 2002; Martin et al., 2003; Ajikumar et al., 2010) and Saccharomyces cerevisiae (Ro et al., 2006; Takahashi et al., 2007; Westfall et al., 2012; Zhuang and Chappell, 2015). The strategies developed in these efforts usually take advantage of specific microbe strains whose innate biosynthetic machinery is genetically modified to accumulate certain prenyldiphosphate precursors (e.g. IPP or FPP), which can be utilized by other introduced terpene synthase(s) for the production of the desired terpene(s). For example, greater than 900 mg L⁻¹ bisabolene was produced when bisabolene synthase genes from plants were introduced into FPP-overproducing E. coli or S. cerevisiae strains (Peralta-Yahya et al., 2011). High levels of farnesane production for diesel fuels were also achieved by reductive hydrogenation of its precursor farnesene, which was generated from a genetically engineered yeast (e.g. Saccharomyces cerevisiae) strain using plant farnesene synthases (Renninger and McPhee, 2008; Ubersax and Platt, 2010). However, terpene production using microbial platforms is still dependent on exogenous feedstocks (i.e. sugars) and elaborate production facilities, both of which add significantly to their production costs.Compared with microbial systems, engineering terpene production in plant systems seems like an attractive target as well. This is because plants can take advantage of photosynthesis by using atmospheric CO₂ as their carbon resource instead of relying on exogenous carbon feedstocks. Moreover, crop plants such as tobacco (Nicotiana tabacum) can generate a large amount of green tissues efficiently when grown for biomass production (Schillberg et al., 2003; Andrianov et al., 2010), making them a robust, sustainable, and scalable platform for large-scale terpene production. Nonetheless, compared with microbial platforms, there are only a few examples of elevating terpene production in bioengineered plants. This is due partly to higher plants being complex multicellular organisms, in which terpene metabolism generally utilizes more complex innate machinery that can be compartmentalized intracellularly and to cell/tissue specificities (Lange and Ahkami, 2013; Kempinski et al., 2015). Significant efforts have been made to overcome these obstacles to improve the production of valuable terpenes in plants, including monoterpenes (Lücker et al., 2004; Ohara et al., 2010; Lange et al., 2011), sesquiterpenes (Aharoni et al., 2003; Kappers et al., 2005; Wu et al., 2006; Davidovich-Rikanati et al., 2008), diterpenes (Besumbes et al., 2004; Anterola et al., 2009), and triterpenes (Inagaki et al., 2011; Wu et al., 2012). Among these, engineering terpene metabolism into a subcellular organelle, where the engineered enzymes/pathways can utilize unlimited/unregulated precursors as substrates, appears most successful. For example, Wu et al. (2006, 2012) expressed an avian farnesyl diphosphate synthase (FPS) with foreign sesquiterpene/triterpene synthases targeted to the plastid to divert the IPP/DMAPP pool from the plastidic MEP pathway to synthesize high levels of the novel sesquiterpenes patchoulol and amorpha-4,11-diene up to 30 µg g⁻¹ fresh weight and the triterpene squalene up to 1,000 µg g⁻¹ fresh weight. This strategy appears to be particularly robust because it avoids possible endogenous regulation of sesquiterpene and triterpene biosynthesis, which occurs normally in the cytoplasm, and relies upon more plastic precursor pools of IPP/DMAPP inherent in the plastid, which are primarily derived from the local CO₂ fixation (Wright et al., 2014).The goal of this study was to evaluate the prospects for engineering advanced features of triterpene metabolism from B. braunii into tobacco and, thus, to probe the innate intricacies of isoprenoid metabolism in plants. In order to achieve this, we first introduced the key steps of botryococcene biosynthesis into specific subcellular compartments of tobacco cells under the direction of constitutive or trichome-specific promoters. The transgenic lines expressing the enzymes in the chloroplast were found to accumulate the highest levels of botryococcene. Triterpene methyltransferases were next introduced into the same intracellular compartments of selected high-triterpene-accumulating lines. A high yield of methylated triterpenes was also achieved in transgenic lines when the TMTs were targeted to the chloroplast. Through careful comparison of the levels of triterpenes and the methylated triterpene products in the various transgenic lines, we have also gained a deeper insight into the subcellular distribution of the triterpene products in these transgenic lines as well as a better understanding of methylation metabolism for specialized metabolites in particular compartments. These findings all contribute to our understanding of the regulatory elements that control carbon flux through the innate terpene biosynthetic pathways operating in plants.