Differences in gene expression within a striking phenotypic mosaic Eucalyptus tree that varies in susceptibility to herbivory

BACKGROUND Long-lived trees can accumulate mutations throughout their lifetimes that may influence biotic and abiotic interactions. For example, some Eucalyptus trees display marked variation in herbivore defence within a single canopy. These "mosaic" trees support foliage with distinct chemotypes which are differentially favoured by insect and vertebrate herbivores, resulting in susceptible and resistant branches within a single canopy. These mosaic trees provide a unique opportunity to explore the biosynthesis and genetic regulation of chemical defences in the foliage. The biosynthesis of the principal defence compounds, terpenoid-dominated essential oils, is well understood. However, the regulation of the genes involved and thus the control of phenotypic variation within a single tree canopy remains a mystery. RESULTS We sequenced the transcriptomes of the leaves of the two different chemotypes of a chemically mosaic Eucalyptus melliodora tree using 454 pyrosequencing technology. We used gene set enrichment analysis to identify differentially expressed transcripts and found the proportion of differentially expressed genes in the resistant and susceptible foliage similar to the transcript difference between functionally distinct tissues of the same organism, for example roots and leaves. We also investigated sequence differences in the form of single nucleotide polymorphisms and found 10 nucleotides that were different between the two branches. These are likely true SNPs and several occur in regulatory genes. CONCLUSION We found three lines of evidence that suggest changes to a 'master switch' can result in large scale phenotypic changes: 1. We found differential expression of terpene biosynthetic genes between the two chemotypes that could contribute to chemical variation within this plant. 2. We identified many genes that are differentially expressed between the two chemotypes, including some unique genes in each branch. These genes are involved in a variety of processes within the plant and many could contribute to the regulation of secondary metabolism, thus contributing to the chemical variation. 3. We identified 10 SNPs, some of which occur in regulatory genes that could influence secondary metabolism and thus contribute to chemical variation. Whilst this research is inherently limited by sample size, the patterns we describe could be indicative of other plant genetic mosaics.This work was supported by a Discovery grant from the Australian Research Council to WJF (DP0877063)

Accuracy of Genomic Prediction for Foliar Terpene Traits in Eucalyptus polybractea

Unlike agricultural crops, most forest species have not had millennia of improvement through phenotypic selection, but can contribute energy and material resources and possibly help alleviate climate change. Yield gains similar to those achieved in agricultural crops over millennia could be made in forestry species with the use of genomic methods in a much shorter time frame. Here we compare various methods of genomic prediction for eight traits related to foliar terpene yield in Eucalyptus polybractea, a tree grown predominantly for the production of Eucalyptus oil. The genomic markers used in this study are derived from shallow whole genome sequencing of a population of 480 trees. We compare the traditional pedigree-based additive best linear unbiased predictors (ABLUP), genomic BLUP (GBLUP), BayesB genomic prediction model, and a form of GBLUP based on weighting markers according to their influence on traits (BLUP|GA). Predictive ability is assessed under varying marker densities of 10,000, 100,000 and 500,000 SNPs. Our results show that BayesB and BLUP|GA perform best across the eight traits. Predictive ability was higher for individual terpene traits, such as foliar α-pinene and 1,8-cineole concentration (0.59 and 0.73, respectively), than aggregate traits such as total foliar oil concentration (0.38). This is likely a function of the trait architecture and markers used. BLUP|GA was the best model for the two biomass related traits, height and 1 year change in height (0.25 and 0.19, respectively). Predictive ability increased with marker density for most traits, but with diminishing returns. The results of this study are a solid foundation for yield improvement of essential oil producing eucalypts. New markets such as biopolymers and terpene-derived biofuels could benefit from rapid yield increases in undomesticated oil-producing species.Funding for this project was provided by the Australian Research Council Linkage Program (LP110100184) toWJF, the Rural Industries Research and Development Corporation (RIRDC), Australia. Support was also provided by the Center for BioEnergy Innovation (CBI), a U.S DOE Bioenergy Research Center supported by the DOE office of science

The Yield of Essential Oils in Melaleuca alternifolia (Myrtaceae) Is Regulated through Transcript Abundance of Genes in the MEP Pathway

Medicinal tea tree (Melaleuca alternifolia) leaves contain large amounts of an essential oil, dominated by monoterpenes. Several enzymes of the chloroplastic methylerythritol phosphate (MEP) pathway are hypothesised to act as bottlenecks to the production of monoterpenes. We investigated, whether transcript abundance of genes encoding for enzymes of the MEP pathway were correlated with foliar terpenes in M. alternifolia using a population of 48 individuals that ranged in their oil concentration from 39 -122 mg x g DM(-1). Our study shows that most genes in the MEP pathway are co-regulated and that the expression of multiple genes within the MEP pathway is correlated with oil yield. Using multiple regression analysis, variation in expression of MEP pathway genes explained 87% of variation in foliar monoterpene concentrations. The data also suggest that sesquiterpenes in M. alternifolia are synthesised, at least in part, from isopentenyl pyrophosphate originating from the plastid via the MEP pathway.The work was funded by an Australian Research Council Linkage grant to W.J.F. (LP110100184) with the active partnership of the Australian Tea Tree Industry Association (ATTIA) and a supplementary grant from the Rural Industries Research and Development Corporation. The funders had no role is study design, data collection and analysis, decision to publish, or preparation of the manuscript

Comparative SNP diversity among four Eucalyptus species for genes from secondary metabolite biosynthetic pathways

<p>Abstract</p> <p>Background</p> <p>There is little information about the DNA sequence variation within and between closely related plant species. The combination of re-sequencing technologies, large-scale DNA pools and availability of reference gene sequences allowed the extensive characterisation of single nucleotide polymorphisms (SNPs) in genes of four biosynthetic pathways leading to the formation of ecologically relevant secondary metabolites in <it>Eucalyptus</it>. With this approach the occurrence and patterns of SNP variation for a set of genes can be compared across different species from the same genus.</p> <p>Results</p> <p>In a single GS-FLX run, we sequenced over 103 Mbp and assembled them to approximately 50 kbp of reference sequences. An average sequencing depth of 315 reads per nucleotide site was achieved for all four eucalypt species, <it>Eucalyptus globulus</it>, <it>E. nitens</it>, <it>E. camaldulensis </it>and <it>E. loxophleba</it>. We sequenced 23 genes from 1,764 individuals and discovered 8,631 SNPs across the species, with about 1.5 times as many SNPs per kbp in the introns compared to exons. The exons of the two closely related species (<it>E. globulus </it>and <it>E. nitens</it>) had similar numbers of SNPs at synonymous and non-synonymous sites. These species also had similar levels of SNP diversity, whereas <it>E. camaldulensis </it>and <it>E. loxophleba </it>had much higher SNP diversity. Neither the pathway nor the position in the pathway influenced gene diversity. The four species share between 20 and 43% of the SNPs in these genes.</p> <p>Conclusion</p> <p>By using conservative statistical detection methods, we were confident about the validity of each SNP. With numerous individuals sampled over the geographical range of each species, we discovered one SNP in every 33 bp for <it>E. nitens </it>and one in every 31 bp in <it>E. globulus</it>. In contrast, the more distantly related species contained more SNPs: one in every 16 bp for <it>E. camaldulensis </it>and one in 17 bp for <it>E. loxophleba</it>, which is, to the best of our knowledge, the highest frequency of SNPs described in woody plant species.</p

High marker density GWAS provides novel insights into the genomic architecture of terpene oil yield in Eucalyptus

Terpenoid based essential oils are economically important commodities, yet beyond their biosynthetic pathways, little is known about the genetic architecture of terpene oil yield from plants. Transport, storage, evaporative loss, transcriptional regulation and precursor competition may be important contributors to this complex trait. Here, we associate 2.39 M single nucleotide polymorphisms derived from shallow whole genome sequencing of 468 Eucalyptus polybractea individuals with 12 traits related to the overall terpene yield, eight direct measures of terpene concentration and four biomass‐related traits. Our results show that in addition to terpene biosynthesis, development of secretory cavities where terpenes are both synthesised and stored, and transport of terpenes were important components of terpene yield. For sesquiterpene concentrations, the availability of precursors in the cytosol was important. Candidate terpene synthase genes for the production of 1,8‐cineole and α‐pinene, and β‐pinene, (which made up more than 80% of the total terpenes) were functionally characterised as a 1,8‐cineole synthase and a β / α‐pinene synthase. Our results provide novel insights of the genomic architecture of terpene yield and we provide candidate genes for breeding or engineering of crops for biofuels or the production of industrially valuable terpenes

The yield of essential oils in Melaleuca alternifolia (Myrtaceae) is regulated through transcript abundance of genes in the MEP pathway

Medicinal tea tree (Melaleuca alternifolia) leaves contain large amounts of an essential oil, dominated by monoterpenes. Several enzymes of the chloroplastic methylerythritol phosphate (MEP) pathway are hypothesised to act as bottlenecks to the production of monoterpenes. We investigated, whether transcript abundance of genes encoding for enzymes of the MEP pathway were correlated with foliar terpenes in M. alternifolia using a population of 48 individuals that ranged in their oil concentration from 39 -122 mg x g DM(-1). Our study shows that most genes in the MEP pathway are co-regulated and that the expression of multiple genes within the MEP pathway is correlated with oil yield. Using multiple regression analysis, variation in expression of MEP pathway genes explained 87% of variation in foliar monoterpene concentrations. The data also suggest that sesquiterpenes in M. alternifolia are synthesised, at least in part, from isopentenyl pyrophosphate originating from the plastid via the MEP pathway

The significance of feedback de-excitation

During photosynthesis sunlight is absorbed by photosynthetic pigments and converted into organic compounds, such as carbohydrates. Photosynthesis needs to be highly regulated, since both too much and too little light are harmful to plant. If too little light is absorbed, a plant cannot store enough energy, which will have effects on growth and fitness of the plant. With too much light absorbed, a dangerous side reaction of photosynthesis, the production of reactive oxygen species can happen. These reactive oxygen species can damage the proteins in the chloroplast and the lipids of the chloroplast. To avoid the production of reactive oxygen species, plants have evolved many mechanisms, which act on different time-scales and different levels of organization. As a first measure, when the absorbed light is exceeding the capacity for its utilization, is to switch the light-harvesting antenna from efficient light harvesting to energy dissipation. This process is called feedback de-excitation (FDE). The protein PsbS is essential for this process as well as a functioning xanthophylls cycle with the enzyme violaxanthin de-epoxidase (VDE). I have investigated the effects of plants with changes in their ability to dissipate excess excitation energy in the model plants species Arabidopsis thaliana. Three genotypes with either increased or decreased capacity for FDE were used during my experiments. The first genotype over-expresses the PsbS gene, having approximately two-fold increased amounts of PsbS and FDE. The second is a PsbS deletion mutant with no PsbS protein and no FDE. The third genotype cannot perform the conversion of violaxanthin to zeaxanthin, because the enzyme VDE is missing. This mutant has some FDE left. Arabidopsis thaliana is an annual plant, which flowers only once in its lifetime. Therefore, when counting the seeds produced an estimation of fitness can be made from the amount of seeds produced. This was done during my experiments and shown that FDE is a trait and that plants with increased FDE have a higher fitness and vice versa. This was also the case for a collection of plants lacking a single protein from the light harvesting antenna. All of these genotypes had a fitness reduction, proving that their function is not redundant. In an attempt to explain why the fitness is reduced in plants with altered FDE, photosynthetic measurements, as well as a determination of the transcriptome and the metabolome was performed. Plants lacking FDE had higher levels of photoinhibition, leading both to lower rates of photosynthesis and to higher repair cost. This could in part explain the reduction in fitness. These plants also had major changes in their transcriptome and their metabolome. Primary metabolism was most effected, for example carbohydrate and amino acid metabolism. But there were also changes in secondary metabolism such as an up regulation of the biosynthesis of anthocyanins

The genetic basis of foliar terpene yield: implications for breeding and profitability of Australian essential oil crops

The family Myrtaceae is known for its high foliar terpene concentrations as well as significant qualitative and quantitative variation in foliar terpenes between taxa, populations and individuals. To date, few studies have investigated the genetic and biochemical processes, which underlie this variation, much of which is known to be under genetic control. Differences in yield are both ecologically and commercially important and a better understanding of its basis will allow a greater understanding of Australian ecosystems as well as improve commercial viability of essential oil industries. Over the past decade a good understanding of the genes involved in terpene biosynthesis has developed in other species and several important regulatory steps have been identified. Much of this work has been done in transgenic plants, so our understanding at a molecular level is strong. Nonetheless, it remains unclear if these processes are transferrable to wild populations, or indeed how ecologically important quantitative variation in terpenoids arise and are maintained in natural ecosystems. In this review we will summarize what is known about terpene biosynthesis and the control of flux through the terpene biosynthetic pathways. We will then argue that this platform of work provides a great resource for Myrtaceae, as well as other plants, to identify candidate genes that control flux through the biosynthetic pathways and how this will inform further studies into the ecological implications of quantitative variation of terpenes. Work into terpene biosynthesis would also provide a framework to improve the profitability of essential oil crops