A phylogenomic approach reveals a low somatic mutation rate in a long-lived plant.

Somatic mutations can have important effects on the life history, ecology, and evolution of plants, but the rate at which they accumulate is poorly understood and difficult to measure directly. Here, we develop a method to measure somatic mutations in individual plants and use it to estimate the somatic mutation rate in a large, long-lived, phenotypically mosaic Eucalyptus melliodora tree. Despite being 100 times larger than Arabidopsis, this tree has a per-generation mutation rate only ten times greater, which suggests that this species may have evolved mechanisms to reduce the mutation rate per unit of growth. This adds to a growing body of evidence that illuminates the correlated evolutionary shifts in mutation rate and life history in plants.

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 million 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 synthesized 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 comprised > 80% of the total terpenes) were functionally characterized as a 1,8-cineole synthase and a ß/α-pinene synthase. Our results provide novel insights into 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.

Four terpene synthases contribute to the generation of chemotypes in tea tree (Melaleuca alternifolia).

BACKGROUND: Terpene rich leaves are a characteristic of Myrtaceae. There is significant qualitative variation in the terpene profile of plants within a single species, which is observable as "chemotypes". Understanding the molecular basis of chemotypic variation will help explain how such variation is maintained in natural populations as well as allowing focussed breeding for those terpenes sought by industry. The leaves of the medicinal tea tree, Melaleuca alternifolia, are used to produce terpinen-4-ol rich tea tree oil, but there are six naturally occurring chemotypes; three cardinal chemotypes (dominated by terpinen-4-ol, terpinolene and 1,8-cineole, respectively) and three intermediates. It has been predicted that three distinct terpene synthases could be responsible for the maintenance of chemotypic variation in this species. RESULTS: We isolated and characterised the most abundant terpene synthases (TPSs) from the three cardinal chemotypes of M. alternifolia. Functional characterisation of these enzymes shows that they produce the dominant compounds in the foliar terpene profile of all six chemotypes. Using RNA-Seq, we investigated the expression of these and 24 additional putative terpene synthases in young leaves of all six chemotypes of M. alternifolia. CONCLUSIONS: Despite contributing to the variation patterns observed, variation in gene expression of the three TPS genes is not enough to explain all variation for the maintenance of chemotypes. Other candidate terpene synthases as well as other levels of regulation must also be involved. The results of this study provide novel insights into the complexity of terpene biosynthesis in natural populations of a non-model organism.

Transcriptome analysis of terpene chemotypes of Melaleuca alternifolia across different tissues.

Plant chemotypes or chemical polymorphisms are defined by discrete variation in secondary metabolites within a species. This variation can have consequences for ecological interactions or the human use of plants. Understanding the molecular basis of chemotypic variation can help to explain how variation of plant secondary metabolites is controlled. We explored the transcriptomes of the 3 cardinal terpene chemotypes of Melaleuca alternifolia in young leaves, mature leaves, and stem and compared transcript abundance to variation in the constitutive profile of terpenes. Leaves from chemotype 1 plants (dominated by terpinen-4-ol) show a similar pattern of gene expression when compared to chemotype 5 plants (dominated by 1,8-cineole). Only terpene synthases in young leaves were differentially expressed between these chemotypes, supporting the idea that terpenes are mainly synthetized in young tissue. Chemotype 2 plants (dominated by terpinolene) show a greater degree of differential gene expression compared to the other chemotypes, which might be related to the isolation of plant populations that exhibit this chemotype and the possibility that the terpinolene synthase gene in M. alternifolia was derived by introgression from a closely related species, Melaleuca trichostachya. By using multivariate analyses, we were able to associate terpenes with candidate terpene synthases.

The Eucalyptus terpene synthase gene family.

BACKGROUND: Terpenoids are abundant in the foliage of Eucalyptus, providing the characteristic smell as well as being valuable economically and influencing ecological interactions. Quantitative and qualitative inter- and intra- specific variation of terpenes is common in eucalypts. RESULTS: The genome sequences of Eucalyptus grandis and E. globulus were mined for terpene synthase genes (TPS) and compared to other plant species. We investigated the relative expression of TPS in seven plant tissues and functionally characterized five TPS genes from E. grandis. Compared to other sequenced plant genomes, Eucalyptus grandis has the largest number of putative functional TPS genes of any sequenced plant. We discovered 113 and 106 putative functional TPS genes in E. grandis and E. globulus, respectively. All but one TPS from E. grandis were expressed in at least one of seven plant tissues examined. Genomic clusters of up to 20 genes were identified. Many TPS are expressed in tissues other than leaves which invites a re-evaluation of the function of terpenes in Eucalyptus. CONCLUSIONS: Our data indicate that terpenes in Eucalyptus may play a wider role in biotic and abiotic interactions than previously thought. Tissue specific expression is common and the possibility of stress induction needs further investigation. Phylogenetic comparison of the two investigated Eucalyptus species gives insight about recent evolution of different clades within the TPS gene family. While the majority of TPS genes occur in orthologous pairs some clades show evidence of recent gene duplication, as well as loss of function.

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.

Mosaic eucalypt trees suggest genetic control at a point that influences several metabolic pathways.

Mosaic trees contain more than one phenotype. The two Eucalyptus mosaic trees studied here (E. melliodora and E. sideroxylon) are predominantly susceptible to insect herbivory, with the leaves on a single large branch on each tree resisting herbivory. We used gas chromatography-mass spectrometry and high-pressure liquid chromatography to analyze the chemical profile of leaves of the mosaic trees, as well as leaves of adjacent non-mosaic con-specifics. We show that the leaves of the two phenotypes are distinctly different. Compared to the susceptible (S) leaves on the same tree, the resistant (R) leaves on the mosaic trees had low concentrations of sesquiterpenes (E. melliodora: 2 vs. 24 mg·g(-1) dry matter; E. sideroxylon: 8 vs. 22 mg·g(-1) dry matter), high concentrations of FPCs (E. melliodora: 5.4 vs. 0.3 mg·g(-1) dry matter; E. sideroxylon: 9.8 vs. 0.2 mg·g(-1) dry matter) but similar concentrations of nitrogen (E. melliodora: 15.4 vs. 16.8 mg·g(-1) dry matter; E. sideroxylon: 13.1 vs. 14.0 mg·g(-1) dry matter). The only difference between the two mosaic trees was in the levels of monoterpenes. The R leaves from the mosaic E. melliodora contained higher concentrations of monoterpenes compared to the S leaves (12 vs. 6 mg·g(-1) dry matter). In contrast, the leaves from the E. sideroxylon mosaic contained much higher concentrations of monoterpenes with a reversed pattern (R: 26 vs. S: 45 mg·g(-1) dry matter). There were qualitative differences too on the mosaic trees. The R leaves of both species contained much higher concentrations of the monoterpene, 1,8-cineole, whereas the S chemotype of E. sideroxylon only had high concentrations of phellandrenes. Furthermore, the chemical differences between leaves on the R and S branches of the mosaic trees resembled those between the leaves of R and S con-specific trees in the same population. We use these data and knowledge of secondary metabolite biosynthesis to propose that high-level transcriptional differences may control the profile of specialized metabolites in eucalypts.

Global population genomic signature of Spodoptera frugiperda (fall armyworm) supports complex introduction events across the Old World.

Native to the Americas, the invasive Spodoptera frugiperda (fall armyworm; FAW) was reported in West Africa in 2016, followed by its chronological detection across the Old World and the hypothesis of an eastward Asia expansion. We explored population genomic signatures of American and Old World FAW and identified 12 maternal mitochondrial DNA genome lineages across the invasive range. 870 high-quality nuclear single nucleotide polymorphic DNA markers identified five distinct New World population clusters, broadly reflecting FAW native geographical ranges and the absence of host-plant preferences. We identified unique admixed Old World populations, and admixed and non-admixed Asian FAW individuals, all of which suggested multiple introductions underpinning the pest's global spread. Directional gene flow from the East into eastern Africa was also detected, in contrast to the west-to-east spread hypothesis. Our study demonstrated the potential of population genomic approaches via international partnership to address global emerging pest threats and biosecurity challenges.

Whole-genome sequencing to detect mutations associated with resistance to insecticides and Bt proteins in Spodoptera frugiperda.

The fall armyworm (FAW), Spodoptera frugiperda, is a major pest native to the Americas that has recently invaded the Old World. Point mutations in the target-site proteins acetylcholinesterase-1 (ace-1), voltage-gated sodium channel (VGSC) and ryanodine receptor (RyR) have been identified in S. frugiperda as major resistance mechanisms to organophosphate, pyrethroid and diamide insecticides respectively. Mutations in the adenosine triphosphate-binding cassette transporter C2 gene (ABCC2) have also been identified to confer resistance to Cry1F protein. In this study, we applied a whole-genome sequencing (WGS) approach to identify point mutations in the target-site genes in 150 FAW individuals collected from China, Malawi, Uganda and Brazil. This approach revealed three amino acid substitutions (A201S, G227A and F290V) of S. frugiperda ace-1, which are known to be associated with organophosphate resistance. The Brazilian population had all three ace-1 point mutations and the 227A allele (mean frequency = 0.54) was the most common. Populations from China, Malawi and Uganda harbored two of the three ace-1 point mutations (A201S and F290V) with the 290V allele (0.47-0.58) as the dominant allele. Point mutations in VGSC (T929I, L932F and L1014F) and RyR (I4790M and G4946E) were not detected in any of the 150 individuals. A novel 12-bp insertion mutation in exon 15 of the ABCC2 gene was identified in some of the Brazilian individuals but absent in the invasive populations. Our results not only demonstrate robustness of the WGS-based genomic approach for detection of resistance mutations, but also provide insights for improvement of resistance management tactics in S. frugiperda.

Genetic evidence supports three previously described species of greater glider, Petauroides volans, P. minor, and P. armillatus.

The identification and classification of species are essential for effective conservation management. This year, Australia experienced a bushfire season of unprecedented severity, resulting in widespread habitat loss and mortality. As a result, there has been an increased focus on understanding genetic diversity and structure across the range of individual species to protect resilience in the face of climate change. The greater glider (Petauroides volans) is a large, gliding eucalypt folivore. This nocturnal arboreal marsupial has a wide distribution across eastern Australia and is considered the sole extant member of the genus Petauroides. Differences in morphology have led to suggestions that the one accepted species is actually three. This would have substantial impacts on conservation management, particularly given a recent history of declining populations, coupled with extensive wildfires. Until now, genetic evidence to support multiple species has been lacking. For the first time, we used DArT sequencing on greater glider tissue samples from multiple regions and found evidence of three operational taxonomic units (OTUs) representing northern, central and southern groups. The three OTUs were also supported by our morphological data. These findings have important implications for greater glider management and highlight the role of genetics in helping to assess conservation status.

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.

Transcriptome sequencing of two phenotypic mosaic Eucalyptus trees reveals large scale transcriptome re-modelling.

Phenotypic mosaic trees offer an ideal system for studying differential gene expression. We have investigated two mosaic eucalypt trees from two closely related species (Eucalyptus melliodora and E. sideroxylon), which each support two types of leaves: one part of the canopy is resistant to insect herbivory and the remaining leaves are susceptible. Driving this ecological distinction are differences in plant secondary metabolites. We used these phenotypic mosaics to investigate genome wide patterns of foliar gene expression with the aim of identifying patterns of differential gene expression and the somatic mutation(s) that lead to this phenotypic mosaicism. We sequenced the mRNA pool from leaves of the resistant and susceptible ecotypes from both mosaic eucalypts using the Illumina HiSeq 2000 platform. We found large differences in pathway regulation and gene expression between the ecotypes of each mosaic. The expression of the genes in the MVA and MEP pathways is reflected by variation in leaf chemistry, however this is not the case for the terpene synthases. Apart from the terpene biosynthetic pathway, there are several other metabolic pathways that are differentially regulated between the two ecotypes, suggesting there is much more phenotypic diversity than has been described. Despite the close relationship between the two species, they show large differences in the global patterns of gene and pathway regulation.

The molecular basis of host plant selection in Melaleuca quinquenervia by a successful biological control agent.

Melaleuca quinquenervia (Cav) S.T. Blake (broadleaf paperbark) is an Australian tree that has become a serious weed in many places around the world. Two insects Oxyops vitiosa (the melaleuca weevil), and Boreioglycaspis melaleucae (the melaleuca psyllid), which were introduced to Florida as part of a biological control programme, have been very effective in reducing survival and reproduction of this weed. There are two terpene chemotypes of M. quinquenervia; one rich in the sesquiterpene E-nerolidol whereas the other is rich in viridiflorol. Viridiflorol is a strong feeding deterrent for the melaleuca weevil and retards larval development. The larvae therefore avoid the viridiflorol-rich chemotype, in contrast, female melaleuca psyllids prefer to oviposit on these leaves. To identify the molecular basis of these preferences, we isolated and characterised two terpene synthases from the viridiflorol-rich chemotype, both of which utilise farnesyl pyrophosphate and have the same product profile. Chemotypic variation in terpenes in M. quinquenervia is under strong genetic control and the reproductive potential of each chemotype is limited by a different insect. These insects could, therefore, be selective agents for the maintenance of chemotypic variation in M. quinquenervia.