Comparative Genomics Analyses Reveal Extensive Chromosome Colinearity and Novel Quantitative Trait Loci in Eucalyptus

Dense genetic maps, along with quantitative trait loci (QTLs) detected on such maps, are powerful tools for genomics and molecular breeding studies. In the important woody genus Eucalyptus, the recent release of E. grandis genome sequence allows for sequence-based genomic comparison and searching for positional candidate genes within QTL regions. Here, dense genetic maps were constructed for E. urophylla and E. tereticornis using genomic simple sequence repeats (SSR), expressed sequence tag (EST) derived SSR, EST-derived cleaved amplified polymorphic sequence (EST-CAPS), and diversity arrays technology (DArT) markers. The E. urophylla and E. tereticornis maps comprised 700 and 585 markers across 11 linkage groups, totaling at 1,208.2 and 1,241.4 cM in length, respectively. Extensive synteny and colinearity were observed as compared to three earlier DArT-based eucalypt maps (two maps with E. grandis × E. urophylla and one map of E. globulus) and with the E. grandis genome sequence. Fifty-three QTLs for growth (10-56 months of age) and wood density (56 months) were identified in 22 discrete regions on both maps, in which only one colocalizaiton was found between growth and wood density. Novel QTLs were revealed as compared with those previously detected on DArT-based maps for similar ages in Eucalyptus. Eleven to 585 positional candidate genes were obained for a 56-month-old QTL through aligning QTL confidence interval with the E. grandis genome. These results will assist in comparative genomics studies, targeted gene characterization, and marker-assisted selection in Eucalyptus and the related taxa.

Concerted transcription of auxin and carbohydrate homeostasis-related genes underlies improved adventitious rooting of microcuttings derived from far-red treated Eucalyptus globulus Labill mother plants

Economically important plant species, such as Eucalyptus globulus, are often rooting recalcitrant. We have previously shown that far-red light enrichment applied to E. globulus donor-plants improved microcutting rooting competence and increased rooting zone/shoot carbohydrate ratio. To better understand this developmental response, the relative expression profiles of genes involved in auxin signaling (ARF6, ARF8, AGO1), biosynthesis (YUC3) and transport (AUX1, PIN1, PIN2); sucrose cleavage (SUS1, CWINV1), transport (SUC5), hexose phosphorylation (HXK1, FLN1) and starch biosynthesis (SS3) were quantified during adventitious rooting of E. globulus microcuttings derived from donor plants exposed to far-red or white light. Expression of auxin transport-related genes increased in the first days of root induction. Far-red enrichment of donor plants induced ARF6, ARF8 and AGO1 in microcuttings. The first two gene products could activate GH3 and other rooting related genes, whereas AGO1 deregulation of the repressor ARF17 may relief rooting inhibition. Increased sink strength at the basal stem with sucrose unloading in root tissue mediated by SUC and subsequent hydrolysis by SUS1 were also supported by gene expression profile. Fructose phosphorylation and starch biosynthesis could also contribute to proper carbon allocation at the site of rooting, as evidenced by increased expression of related genes. These data are in good agreement with increased contents of hexoses and starch at the cutting base severed from far-red exposed donor plants. To sum up, pathways integrating auxin and carbohydrate metabolism were activated in microcuttings derived from donor plants exposed to far red light enrichment, thereby improving rooting response in E. globulus.

A climate change context for the decline of a foundation tree species in south-western Australia: insights from phylogeography and species distribution modelling

BACKGROUND AND AIMS

A worldwide increase in tree decline and mortality has been linked to climate change and, where these represent foundation species, this can have important implications for ecosystem functions. This study tests a combined approach of phylogeographic analysis and species distribution modelling to provide a climate change context for an observed decline in crown health and an increase in mortality in Eucalyptus wandoo, an endemic tree of south-western Australia.

METHODS

Phylogeographic analyses were undertaken using restriction fragment length polymorphism analysis of chloroplast DNA in 26 populations across the species distribution. Parsimony analysis of haplotype relationships was conducted, a haplotype network was prepared, and haplotype and nucleotide diversity were calculated. Species distribution modelling was undertaken using Maxent models based on extant species occurrences and projected to climate models of the last glacial maximum (LGM).

KEY RESULTS

A structured pattern of diversity was identified, with the presence of two groups that followed a climatic gradient from mesic to semi-arid regions. Most populations were represented by a single haplotype, but many haplotypes were shared among populations, with some having widespread distributions. A putative refugial area with high haplotype diversity was identified at the centre of the species distribution. Species distribution modelling showed high climatic suitability at the LGM and high climatic stability in the central region where higher genetic diversity was found, and low suitability elsewhere, consistent with a pattern of range contraction.

CONCLUSIONS

Combination of phylogeography and paleo-distribution modelling can provide an evolutionary context for climate-driven tree decline, as both can be used to cross-validate evidence for refugia and contraction under harsh climatic conditions. This approach identified a central refugial area in the test species E. wandoo, with more recent expansion into peripheral areas from where it had contracted at the LGM. This signature of contraction from lower rainfall areas is consistent with current observations of decline on the semi-arid margin of the range, and indicates low capacity to tolerate forecast climatic change. Identification of a paleo-historical context for current tree decline enables conservation interventions to focus on maintaining genetic diversity, which provides the evolutionary potential for adaptation to climate change.

Genome-wide patterns of recombination, linkage disequilibrium and nucleotide diversity from pooled resequencing and single nucleotide polymorphism genotyping unlock the evolutionary history of Eucalyptus grandis

We used high-density single nucleotide polymorphism (SNP) data and whole-genome pooled resequencing to examine the landscape of population recombination (ρ) and nucleotide diversity (ϴw ), assess the extent of linkage disequilibrium (r(2) ) and build the highest density linkage maps for Eucalyptus. At the genome-wide level, linkage disequilibrium (LD) decayed within c. 4-6 kb, slower than previously reported from candidate gene studies, but showing considerable variation from absence to complete LD up to 50 kb. A sharp decrease in the estimate of ρ was seen when going from short to genome-wide inter-SNP distances, highlighting the dependence of this parameter on the scale of observation adopted. Recombination was correlated with nucleotide diversity, gene density and distance from the centromere, with hotspots of recombination enriched for genes involved in chemical reactions and pathways of the normal metabolic processes. The high nucleotide diversity (ϴw = 0.022) of E. grandis revealed that mutation is more important than recombination in shaping its genomic diversity (ρ/ϴw = 0.645). Chromosome-wide ancestral recombination graphs allowed us to date the split of E. grandis (1.7-4.8 million yr ago) and identify a scenario for the recent demographic history of the species. Our results have considerable practical importance to Genome Wide Association Studies (GWAS), while indicating bright prospects for genomic prediction of complex phenotypes in eucalypt breeding.

Comparative transcriptional analysis provides new insights into the molecular basis of adventitious rooting recalcitrance in Eucalyptus

Adventitious rooting (AR) is essential in clonal propagation. Eucalyptus globulus is relevant for the cellulose industry due to its low lignin content. However, several useful clones are recalcitrant to AR, often requiring exogenous auxin, adding cost to clonal garden operations. In contrast, E. grandis is an easy-to-root species widely used in clonal forestry. Aiming at contributing to the elucidation of recalcitrance causes in E. globulus, we conducted a comparative analysis with these two species differing in rooting competence, combining gene expression and anatomical techniques. Recalcitrance in E. globulus is reversed by exposure to exogenous indole-3-acetic acid (IAA), which promotes important gene expression modifications in both species. The endogenous content of IAA was significantly higher in E. grandis than in E. globulus. The cambium zone was identified as an active area during AR, concentrating the first cell divisions. Immunolocalization assay showed auxin accumulation in cambium cells, further indicating the importance of this region for rooting. We then performed a cambium zone-specific gene expression analysis during AR using laser microdissection. The results indicated that the auxin-related genes TOPLESS and IAA12/BODENLOS and the cytokinin-related gene ARR1may act as negative regulators of AR, possibly contributing to the hard-to-root phenotype of E. globulus.

Assessment of physiological and biochemical responses, metal tolerance and accumulation in two eucalypt hybrid clones for phytoremediation of cadmium-contaminated waters

Eucalyptus is a promising species for ecological restoration but plant performances under environmental constraints need to be better investigated. In particular, the toxic effects of metals on this plant species are poorly described in the literature. In this work, morpho-physiological and biochemical responses to cadmium were analysed in two eucalypt genotypes (hybrid clones of Eucalyptus camaldulensis × Eucalyptus globulus ssp. bicostata J.B. Kirkp named Velino ex 7 and Viglio ex 358) exposed for 3 weeks to 50 μM CdSO4 under hydroponics. The two eucalypt clones showed a different sensitivity to the metal. The growth reduction caused by cadmium was less than 30% in clone Velino and about 50% in clone Viglio. Cadmium mostly accumulated in plant roots and, to a lesser extent, in stem, as highlighted by the low translocation factor (Tf) measured in both clones. Net photosynthesis measurement, chlorophyll fluorescence images, transpiration values and chlorophyll content revealed a cadmium-induced impairment of physiological processes at the leaf level, which was more evident in clone Viglio. Metal binding and antioxidative compound content was differentially affected by cadmium exposure in the two eucalypt clones. Particularly, the content of thiols like cysteine and glutathione, organic acids like oxalate and citrate, and polyamines were markedly modulated in plant organs by metal treatment and highlighted different defence responses between the clones. Cadmium tolerance and accumulation ability of the eucalypt clones were evaluated and the potential of E. camaldulensis for the reclamation of metal polluted-waters is discussed.

Transcriptionally active LTR retrotransposons in Eucalyptus genus are differentially expressed and insertionally polymorphic

BACKGROUND

In Eucalyptus genus, studies on genome composition and transposable elements (TEs) are particularly scarce. Nearly half of the recently released Eucalyptus grandis genome is composed by retrotransposons and this data provides an important opportunity to understand TE dynamics in Eucalyptus genome and transcriptome.

RESULTS

We characterized nine families of transcriptionally active LTR retrotransposons from Copia and Gypsy superfamilies in Eucalyptus grandis genome and we depicted genomic distribution and copy number in two Eucalyptus species. We also evaluated genomic polymorphism and transcriptional profile in three organs of five Eucalyptus species. We observed contrasting genomic and transcriptional behavior in the same family among different species. RLC_egMax_1 was the most prevalent family and RLC_egAngela_1 was the family with the lowest copy number. Most families of both superfamilies have their insertions occurring <3 million years, except one Copia family, RLC_egBianca_1. Protein theoretical models suggest different properties between Copia and Gypsy domains. IRAP and REMAP markers suggested genomic polymorphisms among Eucalyptus species. Using EST analysis and qRT-PCRs, we observed transcriptional activity in several tissues and in all evaluated species. In some families, osmotic stress increases transcript values.

CONCLUSION

Our strategy was successful in isolating transcriptionally active retrotransposons in Eucalyptus, and each family has a particular genomic and transcriptional pattern. Overall, our results show that retrotransposon activity have differentially affected genome and transcriptome among Eucalyptus species.

Mating system and early viability resistance to habitat fragmentation in a bird-pollinated eucalypt

Habitat fragmentation has been shown to disrupt ecosystem processes such as plant-pollinator mutualisms. Consequently, mating patterns in remnant tree populations are expected to shift towards increased inbreeding and reduced pollen diversity, with fitness consequences for future generations. However, mating patterns and phenotypic assessments of open-pollinated progeny have rarely been combined in a single study. Here, we collected seeds from 37 Eucalyptus incrassata trees from contrasting stand densities following recent clearance in a single South Australian population (intact woodland=12.6 trees ha(-1); isolated pasture=1.7 trees ha(-1); population area=10 km(2)). 649 progeny from these trees were genotyped at eight microsatellite loci. We estimated genetic diversity, spatial genetic structure, indirect contemporary pollen flow and mating patterns for adults older than the clearance events and open-pollinated progeny sired post-clearance. A proxy of early stage progeny viability was assessed in a common garden experiment. Density had no impact on mating patterns, adult and progeny genetic diversity or progeny growth, but was associated with increased mean pollen dispersal. Weak spatial genetic structure among adults suggests high historical gene flow. We observed preliminary evidence for inbreeding depression related to stress caused by fungal infection, but which was not associated with density. Higher observed heterozygosities in adults compared with progeny may relate to weak selection on progeny and lifetime-accumulated mortality of inbred adults. E. incrassata appears to be resistant to the negative mating pattern and fitness changes expected within fragmented landscapes. This pattern is likely explained by strong outcrossing and regular long-distance pollen flow.

Mating patterns and pollinator mobility are critical traits in forest fragmentation genetics

Most woody plants are animal-pollinated, but the global problem of habitat fragmentation is changing the pollination dynamics. Consequently, the genetic diversity and fitness of the progeny of animal-pollinated woody plants sired in fragmented landscapes tend to decline due to shifts in plant-mating patterns (for example, reduced outcrossing rate, pollen diversity). However, the magnitude of this mating-pattern shift should theoretically be a function of pollinator mobility. We first test this hypothesis by exploring the mating patterns of three ecologically divergent eucalypts sampled across a habitat fragmentation gradient in southern Australia. We demonstrate increased selfing and decreased pollen diversity with increased fragmentation for two small-insect-pollinated eucalypts, but no such relationship for the mobile-bird-pollinated eucalypt. In a meta-analysis, we then show that fragmentation generally does increase selfing rates and decrease pollen diversity, and that more mobile pollinators tended to dampen these mating-pattern shifts. Together, our findings support the premise that variation in pollinator form contributes to the diversity of mating-pattern responses to habitat fragmentation.

Parentage Reconstruction in Eucalyptus nitens Using SNPs and Microsatellite Markers: A Comparative Analysis of Marker Data Power and Robustness

Pedigree reconstruction using molecular markers enables efficient management of inbreeding in open-pollinated breeding strategies, replacing expensive and time-consuming controlled pollination. This is particularly useful in preferentially outcrossed, insect pollinated Eucalypts known to suffer considerable inbreeding depression from related matings. A single nucleotide polymorphism (SNP) marker panel consisting of 106 markers was selected for pedigree reconstruction from the recently developed high-density Eucalyptus Infinium SNP chip (EuCHIP60K). The performance of this SNP panel for pedigree reconstruction in open-pollinated progenies of two Eucalyptus nitens seed orchards was compared with that of two microsatellite panels with 13 and 16 markers respectively. The SNP marker panel out-performed one of the microsatellite panels in the resolution power to reconstruct pedigrees and out-performed both panels with respect to data quality. Parentage of all but one offspring in each clonal seed orchard was correctly matched to the expected seed parent using the SNP marker panel, whereas parentage assignment to less than a third of the expected seed parents were supported using the 13-microsatellite panel. The 16-microsatellite panel supported all but one of the recorded seed parents, one better than the SNP panel, although there was still a considerable level of missing and inconsistent data. SNP marker data was considerably superior to microsatellite data in accuracy, reproducibility and robustness. Although microsatellites and SNPs data provide equivalent resolution for pedigree reconstruction, microsatellite analysis requires more time and experience to deal with the uncertainties of allele calling and faces challenges for data transferability across labs and over time. While microsatellite analysis will continue to be useful for some breeding tasks due to the high information content, existing infrastructure and low operating costs, the multi-species SNP resource available with the EuCHIP60k, opens a whole new array of opportunities for high-throughput, genome-wide or targeted genotyping in species of Eucalyptus.

The Transcriptome and Terpene Profile of Eucalyptus grandis Reveals Mechanisms of Defense Against the Insect Pest, Leptocybe invasa

Plants have evolved complex defenses that allow them to protect themselves against pests and pathogens. However, there is relatively little information regarding the Eucalyptus defensome. Leptocybe invasa is one of the most damaging pests in global Eucalyptus forestry, and essentially nothing is known regarding the molecular mechanisms governing the interaction between the pest and host. The aim of the study was to investigate changes in the transcriptional landscape and terpene profile of a resistant and susceptible Eucalyptus genotype in an effort to improve our understanding of this interaction. We used RNA-seqencing to investigate transcriptional changes following L. invasa oviposition. Expression levels were validated using real-time quantitative PCR. Terpene profiles were investigated using gas chromatography coupled to mass spectometry on uninfested and oviposited leaves. We found 698 and 1,115 significantly differentially expressed genes from the resistant and susceptible interactions, respectively. Gene Ontology enrichment and Mapman analyses identified putative defense mechanisms including cell wall reinforcement, protease inhibitors, cell cycle suppression and regulatory hormone signaling pathways. There were significant differences in the mono- and sesquiterpene profiles between genotypes and between control and infested material. A model of the interaction between Eucalyptus and L. invasa was proposed from the transcriptomic and chemical data.

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.

The fasciclin-like arabinogalactan protein family of Eucalyptus grandis contains members that impact wood biology and biomechanics

Fasciclin-like arabinogalactan protein (FLA) families have been identified and characterised in key plant species, with some members exhibiting functional specialization. Here we identify the FLA family of Eucalyptus grandis, and investigate the roles of three single-FAS domain FLAs, with particular focus on secondary cell-wall formation and wood properties. We use various in-silico approaches to identify and characterise E. grandis genome FLAs, and perform phylogenetic comparisons with other species. For three key FLAs, we perform functional testing including promoter-reporter and overexpression transgenic approaches using eucalypts, poplar and tobacco. Of the 18 eucalypt FLAs identified, several were specifically and highly expressed in stems. The specificity to stem xylem vessel and fibre development was demonstrated with EniFLA1promoter:GUS studies in several species. Testing of select eucalypt FLAs resulted in altered wood development and properties, for example 35S:EgrFLA2 led to a 3 degree reduction in cellulose microfibril angle in eucalypt xylem fibres, and 35S:EgrFLA3 to a reduction in tobacco stem flexural strength. These results indicate that the eucalypt FLA family contains diverse members, and particular members with single FAS domains that are functionally specialized for secondary cell wall growth and properties.

Genome-wide analysis of the lignin toolbox of Eucalyptus grandis

Lignin, a major component of secondary cell walls, hinders the optimal processing of wood for industrial uses. The recent availability of the Eucalyptus grandis genome sequence allows comprehensive analysis of the genes encoding the 11 protein families specific to the lignin branch of the phenylpropanoid pathway and identification of those mainly involved in xylem developmental lignification. We performed genome-wide identification of putative members of the lignin gene families, followed by comparative phylogenetic studies focusing on bona fide clades inferred from genes functionally characterized in other species. RNA-seq and microfluid real-time quantitative PCR (RT-qPCR) expression data were used to investigate the developmental and environmental responsive expression patterns of the genes. The phylogenetic analysis revealed that 38 E. grandis genes are located in bona fide lignification clades. Four multigene families (shikimate O-hydroxycinnamoyltransferase (HCT), p-coumarate 3-hydroxylase (C3H), caffeate/5-hydroxyferulate O-methyltransferase (COMT) and phenylalanine ammonia-lyase (PAL)) are expanded by tandem gene duplication compared with other plant species. Seventeen of the 38 genes exhibited strong, preferential expression in highly lignified tissues, probably representing the E. grandis core lignification toolbox. The identification of major genes involved in lignin biosynthesis in E. grandis, the most widely planted hardwood crop world-wide, provides the foundation for the development of biotechnology approaches to develop tree varieties with enhanced processing qualities.

Protein domain evolution is associated with reproductive diversification and adaptive radiation in the genus Eucalyptus

Eucalyptus is a pivotal genus within the rosid order Myrtales with distinct geographic history and adaptations. Comparative analysis of protein domain evolution in the newly sequenced Eucalyptus grandis genome and other rosid lineages sheds light on the adaptive mechanisms integral to the success of this genus of woody perennials. We reconstructed the ancestral domain content to elucidate the gain, loss and expansion of protein domains and domain arrangements in Eucalyptus in the context of rosid phylogeny. We used functional gene ontology (GO) annotation of genes to investigate the possible biological and evolutionary consequences of protein domain expansion. We found that protein modulation within the angiosperms occurred primarily on the level of expansion of certain domains and arrangements. Using RNA-Seq data from E. grandis, we showed that domain expansions have contributed to tissue-specific expression of tandemly duplicated genes. Our results indicate that tandem duplication of genes, a key feature of the Eucalyptus genome, has played an important role in the expansion of domains, particularly in proteins related to the specialization of reproduction and biotic and abiotic interactions affecting root and floral biology, and that tissue-specific expression of proteins with expanded domains has facilitated subfunctionalization in domain families.

Genetic architecture of carbon isotope composition and growth in Eucalyptus across multiple environments

In the context of climate change, the water-use efficiency (WUE) of highly productive tree varieties, such as eucalypts, has become a major issue for breeding programmes. This study set out to dissect the genetic architecture of carbon isotope composition (δ(13) C), a proxy of WUE, across several environments. A family of Eucalyptus urophylla × E. grandis was planted in three trials and phenotyped for δ(13) C and growth traits. High-resolution genetic maps enabled us to target genomic regions underlying δ(13) C quantitative trait loci (QTLs) on the E. grandis genome. Of the 15 QTLs identified for δ(13) C, nine were stable across the environments and three displayed significant QTL-by-environment interaction, suggesting medium to high genetic determinism for this trait. Only one colocalization was found between growth and δ(13) C. Gene ontology (GO) term enrichment analysis suggested candidate genes related to foliar δ(13) C, including two involved in the regulation of stomatal movements. This study provides the first report of the genetic architecture of δ(13) C and its relation to growth in Eucalyptus. The low correlations found between the two traits at phenotypic and genetic levels suggest the possibility of improving the WUE of Eucalyptus varieties without having an impact on breeding for growth.

Genomic patterns of species diversity and divergence in Eucalyptus

We examined genome-wide patterns of DNA sequence diversity and divergence among six species of the important tree genus Eucalyptus and investigated their relationship with genomic architecture. Using c. 90 range-wide individuals of each Eucalyptus species (E. grandis, E. urophylla, E. globulus, E. nitens, E. dunnii and E. camaldulensis), genetic diversity and divergence were estimated from 2840 polymorphic diversity arrays technology markers covering the 11 chromosomes. Species differentiating markers (SDMs) identified in each of 15 pairwise species comparisons, along with species diversity (HHW ) and divergence (FST ), were projected onto the E. grandis reference genome. Across all species comparisons, SDMs totalled 1.1-5.3% of markers and were widely distributed throughout the genome. Marker divergence (FST and SDMs) and diversity differed among and within chromosomes. Patterns of diversity and divergence were broadly conserved across species and significantly associated with genomic features, including the proximity of markers to genes, the relative number of clusters of tandem duplications, and gene density within or among chromosomes. These results suggest that genomic architecture influences patterns of species diversity and divergence in the genus. This influence is evident across the six species, encompassing diverse phylogenetic lineages, geography and ecology.

High-resolution genetic maps of Eucalyptus improve Eucalyptus grandis genome assembly

Genetic maps are key tools in genetic research as they constitute the framework for many applications, such as quantitative trait locus analysis, and support the assembly of genome sequences. The resequencing of the two parents of a cross between Eucalyptus urophylla and Eucalyptus grandis was used to design a single nucleotide polymorphism (SNP) array of 6000 markers evenly distributed along the E. grandis genome. The genotyping of 1025 offspring enabled the construction of two high-resolution genetic maps containing 1832 and 1773 markers with an average marker interval of 0.45 and 0.5 cM for E. grandis and E. urophylla, respectively. The comparison between genetic maps and the reference genome highlighted 85% of collinear regions. A total of 43 noncollinear regions and 13 nonsynthetic regions were detected and corrected in the new genome assembly. This improved version contains 4943 scaffolds totalling 691.3 Mb of which 88.6% were captured by the 11 chromosomes. The mapping data were also used to investigate the effect of population size and number of markers on linkage mapping accuracy. This study provides the most reliable linkage maps for Eucalyptus and version 2.0 of the E. grandis genome.

Comparative interrogation of the developing xylem transcriptomes of two wood-forming species: Populus trichocarpa and Eucalyptus grandis

Wood formation is a complex developmental process governed by genetic and environmental stimuli. Populus and Eucalyptus are fast-growing, high-yielding tree genera that represent ecologically and economically important species suitable for generating significant lignocellulosic biomass. Comparative analysis of the developing xylem and leaf transcriptomes of Populus trichocarpa and Eucalyptus grandis together with phylogenetic analyses identified clusters of homologous genes preferentially expressed during xylem formation in both species. A conserved set of 336 single gene pairs showed highly similar xylem preferential expression patterns, as well as evidence of high functional constraint. Individual members of multi-gene orthologous clusters known to be involved in secondary cell wall biosynthesis also showed conserved xylem expression profiles. However, species-specific expression as well as opposite (xylem versus leaf) expression patterns observed for a subset of genes suggest subtle differences in the transcriptional regulation important for xylem development in each species. Using sequence similarity and gene expression status, we identified functional homologs likely to be involved in xylem developmental and biosynthetic processes in Populus and Eucalyptus. Our study suggests that, while genes involved in secondary cell wall biosynthesis show high levels of gene expression conservation, differential regulation of some xylem development genes may give rise to unique xylem properties.

The effect of elevated carbon dioxide on the interaction between Eucalyptus grandis and diverse isolates of Pisolithus sp. is associated with a complex shift in the root transcriptome

Using the newly available genome for Eucalyptus grandis, we sought to determine the genome-wide traits that enable this host to form mutualistic interactions with ectomycorrhizal (ECM) Pisolithus sp. and to determine how future predicted concentrations of atmospheric carbon dioxide (CO2 ) will affect this relationship. We analyzed the physiological and transcriptomic responses of E. grandis during colonization by different Pisolithus sp. isolates under conditions of ambient (400 ppm) and elevated (650 ppm) CO2 to tease out the gene expression profiles associated with colonization status. We demonstrate that E. grandis varies in its susceptibility to colonization by different Pisolithus isolates in a manner that is not predictable by geographic origin or the internal transcribed spacer (ITS)-based phylogeny of the fungal partner. Elevated concentrations of CO2 alter the receptivity of E. grandis to Pisolithus, a change that is correlated to a dramatic shift in the transcriptomic profile of the root. These data provide a starting point for understanding how future environmental change may alter the signaling between plants and their ECM partners and is a step towards determining the mechanism behind previously observed shifts in Eucalypt-associated fungal communities exposed to elevated concentrations of atmospheric CO2 .

The floral transcriptome of Eucalyptus grandis

As a step toward functional annotation of genes required for floral initiation and development within the Eucalyptus genome, we used short read sequencing to analyze transcriptomes of floral buds from early and late developmental stages, and compared these with transcriptomes of diverse vegetative tissues, including leaves, roots, and stems. A subset of 4807 genes (13% of protein-coding genes) were differentially expressed between floral buds of either stage and vegetative tissues. A similar proportion of genes were differentially expressed among all tissues. A total of 479 genes were differentially expressed between early and late stages of floral development. Gene function enrichment identified 158 gene ontology classes that were overrepresented in floral tissues, including 'pollen development' and 'aromatic compound biosynthetic process'. At least 40 floral-dominant genes lacked functional annotations and thus may be novel floral transcripts. We analyzed several genes and gene families in depth, including 49 putative biomarkers of floral development, the MADS-box transcription factors, 'S-domain'-receptor-like kinases, and selected gene family members with phosphatidylethanolamine-binding protein domains. Expanded MADS-box gene subfamilies in Eucalyptus grandis included SUPPRESSOR OF OVEREXPRESSION OF CO 1 (SOC1), SEPALLATA (SEP) and SHORT VEGETATIVE PHASE (SVP) Arabidopsis thaliana homologs. These data provide a rich resource for functional and evolutionary analysis of genes controlling eucalypt floral development, and new tools for breeding and biotechnology.

Structural, evolutionary and functional analysis of the NAC domain protein family in Eucalyptus

NAC domain transcription factors regulate many developmental processes and stress responses in plants and vary widely in number and family structure. We analysed the characteristics and evolution of the NAC gene family of Eucalyptus grandis, a fast-growing forest tree in the rosid order Myrtales. NAC domain genes identified in the E. grandis genome were subjected to amino acid sequence, phylogenetic and motif analyses. Transcript abundance in developing tissues and abiotic stress conditions in E. grandis and E. globulus was quantified using RNA-seq and reverse transcription quantitative PCR (RT-qPCR). One hundred and eighty-nine E. grandis NAC (EgrNAC) proteins, arranged into 22 subfamilies, are extensively duplicated in subfamilies associated with stress response. Most EgrNAC genes form tandem duplicate arrays that frequently carry signatures of purifying selection. Sixteen amino acid motifs were identified in EgrNAC proteins, eight of which are enriched in, or unique to, Eucalyptus. New candidates for the regulation of normal and tension wood development and cold responses were identified. This first description of a Myrtales NAC domain family reveals an unique history of tandem duplication in stress-related subfamilies that has likely contributed to the adaptation of eucalypts to the challenging Australian environment. Several new candidates for the regulation of stress, wood formation and tree-specific development are reported.

A flexible multi-species genome-wide 60K SNP chip developed from pooled resequencing of 240 Eucalyptus tree genomes across 12 species

We used whole genome resequencing of pooled individuals to develop a high-density single-nucleotide polymorphism (SNP) chip for Eucalyptus. Genomes of 240 trees of 12 species were sequenced at 3.5× each, and 46 997 586 raw SNP variants were subject to multivariable filtering metrics toward a multispecies, genome-wide distributed chip content. Of the 60 904 SNPs on the chip, 59 222 were genotyped and 51 204 were polymorphic across 14 Eucalyptus species, providing a 96% genome-wide coverage with 1 SNP/12-20 kb, and 47 069 SNPs at ≤ 10 kb from 30 444 of the 33 917 genes in the Eucalyptus genome. Given the EUChip60K multi-species genotyping flexibility, we show that both the sample size and taxonomic composition of cluster files impact heterozygous call specificity and sensitivity by benchmarking against 'gold standard' genotypes derived from deeply sequenced individual tree genomes. Thousands of SNPs were shared across species, likely representing ancient variants arisen before the split of these taxa, hinting to a recent eucalypt radiation. We show that the variable SNP filtering constraints allowed coverage of the entire site frequency spectrum, mitigating SNP ascertainment bias. The EUChip60K represents an outstanding tool with which to address population genomics questions in Eucalyptus and to empower genomic selection, GWAS and the broader study of complex trait variation in eucalypts.

The Eucalyptus grandis R2R3-MYB transcription factor family: evidence for woody growth-related evolution and function

The R2R3-MYB family, one of the largest transcription factor families in higher plants, controls a wide variety of plant-specific processes including, notably, phenylpropanoid metabolism and secondary cell wall formation. We performed a genome-wide analysis of this superfamily in Eucalyptus, one of the most planted hardwood trees world-wide. A total of 141 predicted R2R3-MYB sequences identified in the Eucalyptus grandis genome sequence were subjected to comparative phylogenetic analyses with Arabidopsis thaliana, Oryza sativa, Populus trichocarpa and Vitis vinifera. We analysed features such as gene structure, conserved motifs and genome location. Transcript abundance patterns were assessed by RNAseq and validated by high-throughput quantitative PCR. We found some R2R3-MYB subgroups with expanded membership in E. grandis, V. vinifera and P. trichocarpa, and others preferentially found in woody species, suggesting diversification of specific functions in woody plants. By contrast, subgroups containing key genes regulating lignin biosynthesis and secondary cell wall formation are more conserved across all of the species analysed. In Eucalyptus, R2R3-MYB tandem gene duplications seem to disproportionately affect woody-preferential and woody-expanded subgroups. Interestingly, some of the genes belonging to woody-preferential subgroups show higher expression in the cambial region, suggesting a putative role in the regulation of secondary growth.