The molecular basis of quantitative variation in foliar secondary metabolites in Eucalyptus globulus.

Eucalyptus is characterized by high foliar concentrations of plant secondary metabolites with marked qualitative and quantitative variation within a single species. Secondary metabolites in eucalypts are important mediators of a diverse community of herbivores. We used a candidate gene approach to investigate genetic associations between 195 single nucleotide polymorphisms (SNPs) from 24 candidate genes and 33 traits related to secondary metabolites in the Tasmanian Blue Gum (Eucalyptus globulus). We discovered 37 significant associations (false discovery rate (FDR) Q < 0.05) across 11 candidate genes and 19 traits. The effects of SNPs on phenotypic variation were within the expected range (0.018 < r(2) < 0.061) for forest trees. Whereas most marker effects were nonadditive, two alleles from two consecutive genes in the methylerythritol phosphate pathway (MEP) showed additive effects. This study successfully links allelic variants to ecologically important phenotypes which can have a large impact on the entire community. It is one of very few studies to identify the genetic variants of a foundation tree that influences ecosystem function.

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

BACKGROUND: 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 Eucalyptus. 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. RESULTS: 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, Eucalyptus globulus, E. nitens, E. camaldulensis and E. loxophleba. 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 (E. globulus and E. nitens) had similar numbers of SNPs at synonymous and non-synonymous sites. These species also had similar levels of SNP diversity, whereas E. camaldulensis and E. loxophleba 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. CONCLUSION: 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 E. nitens and one in every 31 bp in E. globulus. In contrast, the more distantly related species contained more SNPs: one in every 16 bp for E. camaldulensis and one in 17 bp for E. loxophleba, which is, to the best of our knowledge, the highest frequency of SNPs described in woody plant species.

Gene expression in Eucalyptus branch wood with marked variation in cellulose microfibril orientation and lacking G-layers.

In response to gravitational stresses, angiosperm trees form tension wood in the upper sides of branches and leaning stems in which cellulose content is higher, microfibrils are typically aligned closely with the fibre axis and the fibres often have a thick inner gelatinous cell wall layer (G-layer). Gene expression was studied in Eucalyptus nitens branches oriented at 45 degrees using microarrays containing 4900 xylem cDNAs, and wood fibre characteristics revealed by X-ray diffraction, chemical and histochemical methods. Xylem fibres in tension wood (upper branch) had a low microfibril angle, contained few fibres with G-layers and had higher cellulose and decreased Klason lignin compared with lower branch wood. Expression of two closely related fasciclin-like arabinogalactan proteins and a beta-tubulin was inversely correlated with microfibril angle in upper and lower xylem from branches. Structural and chemical modifications throughout the secondary cell walls of fibres sufficient to resist tension forces in branches can occur in the absence of G-layer enriched fibres and some important genes involved in responses to gravitational stress in eucalypt xylem are identified.

Identification of quantitative trait loci influencing foliar concentrations of terpenes and formylated phloroglucinol compounds in Eucalyptus nitens.

Leaves of eucalypt species contain a variety of plant secondary metabolites, including terpenoids and formylated phloroglucinol compounds (FPCs). Both terpene and FPC concentrations are quantitative traits that can show large variation within a population and have been shown to be heritable. The molecular genetic basis of this variation is currently unknown. Progeny from a field trial of a three-generation mapping pedigree of Eucalyptus nitens were assayed for terpenes and FPCs. Quantitative trait loci (QTL) analyses were conducted using a map constructed from 296 markers to locate regions of the genome influencing foliar concentrations of these plant secondary compounds. A large number of significant QTL for 14 traits were located across nine linkage groups, with significant clustering of QTL on linkage groups 7, 8 and 9. As expected, QTL for biosynthetically related compounds commonly colocated, but QTL for unrelated monterpenes and FPCs also mapped closely together. Colocation of these QTL with mapped candidate genes from the various biosynthetic pathways, and subsequent use of these genes in association mapping, will assist in determining the causes of variation in plant secondary metabolites in eucalypts.

beta-tubulin affects cellulose microfibril orientation in plant secondary fibre cell walls.

Cellulose microfibrils are the major structural component of plant secondary cell walls. Their arrangement in plant primary cell walls, and its consequent influence on cell expansion and cellular morphology, is directed by cortical microtubules; cylindrical protein filaments composed of heterodimers of alpha- and beta-tubulin. In secondary cell walls of woody plant stems the orientation of cellulose microfibrils influences the strength and flexibility of wood, providing the physical support that has been instrumental in vascular plant colonization of the troposphere. Here we show that a Eucalyptus grandisbeta-tubulin gene (EgrTUB1) is involved in determining the orientation of cellulose microfibrils in plant secondary fibre cell walls. This finding is based on RNA expression studies in mature trees, where we identified and isolated EgrTUB1 as a candidate for association with wood-fibre formation, and on the analysis of somatically derived transgenic wood sectors in Eucalyptus. We show that cellulose microfibril angle (MFA) is correlated with EgrTUB1 expression, and that MFA was significantly altered as a consequence of stable transformation with EgrTUB1. Our findings present an important step towards the production of fibres with altered tensile strength, stiffness and elastic properties, and shed light on one of the molecular mechanisms that has enabled trees to dominate terrestrial ecosystems.

Polymorphisms in cinnamoyl CoA reductase (CCR) are associated with variation in microfibril angle in Eucalyptus spp.

Linkage disequilibrium (LD) mapping using natural populations results in higher resolution of marker-trait associations compared to family-based quantitative trait locus (QTL) studies. Depending on the extent of LD, it is possible to identify alleles within candidate genes associated with a trait. Analysis of a natural mutant in Arabidopsis has shown that mutations in cinnamoyl CoA reductase (CCR), a key lignin gene, affect physical properties of the secondary cell wall such as stiffness and strength. Using this gene, we tested whether LD mapping could identify alleles associated with microfibril angle (MFA), a wood quality trait affecting stiffness and strength of wood. We identified 25 common single-nucleotide polymorphism (SNP) markers in the CCR gene in Eucalyptus nitens. Using single-marker and haplotype analyses in 290 trees from a E. nitens natural population, two haplotypes significantly associated with MFA were found. These results were confirmed in two full-sib families of E. nitens and Eucalyptus globulus. In an effort to understand the functional significance of the SNP markers, we sequenced the cDNA clones and identified an alternatively spliced variant from the significant haplotype region. This study demonstrates that LD mapping can be used to identify alleles associated with wood quality traits in natural populations of trees.

Effects of native forest regeneration practices on genetic diversity in Eucalyptus consideniana.

Impacts of forest harvesting and regeneration practices on genetic diversity in the Australian native forest species Eucalyptus consideniana Maiden (yertchuk) were examined using 29 Mendelian DNA markers (18 RFLPs and 11 microsatellites). Two replicate logging coupes were studied from each of the two most commonly employed silvicultural treatments: clear felling with aerial re-sowing and the seed tree system. For each coupe, genetic diversity measures were compared between a sample of the sapling regeneration and a corresponding control sample from bordering unharvested trees. When calculations were performed over all 29 loci, significant reductions of allelic richness (AR), effective number of alleles (AE) and/or expected heterozygosity (HE) were detected on one or both of seed tree coupes, but on neither of the clear falls. When calculations were performed over the 11 microsatellites alone, all three measures, AR, AE and HE, were significantly reduced on both of the seed replicates but on neither of the two clear falls. In contrast, when the RFLPs were examined separately, there were no significant reductions of diversity on either of the two seed tree coupes or on the two clear falls. These results suggest that genetic erosion is more likely under the seed tree system than under clear-felling with aerial re-sowing and that there is greater statistical power to detect it with microsatellites than with RFLPs. A Monte Carlo simulation to test the statistical significance of the number of apparently lost or gained alleles showed that significant losses of alleles above specified threshold frequencies occurred only in the two seed tree replicates. Three of the four control and regeneration population pairs were significantly differentiated, as indicated by exact tests or by pairwise FST estimates. Comparisons of CONTML dendrograms, constructed for the regeneration populations only versus the control populations only, indicated that genetic drift was significantly promoted under forest management. No significant decreases in observed heterozygosity, or increases in the panmictic index (f), were observed in any of the comparisons suggesting that inbreeding was not promoted by a single rotation of forest management.