Steven G. Hussey

The genome of Eucalyptus grandis

Eucalypts are the world’s most widely planted hardwood trees. Their outstanding diversity, adaptability and growth have made them a global renewable resource of fibre and energy. We sequenced and assembled >94% of the 640-megabase genome of Eucalyptus grandis. Of 36,376 predicted protein-coding genes, 34% occur in tandem duplications, the largest proportion thus far in plant genomes. Eucalyptus also shows the highest diversity of genes for specialized metabolites such as terpenes that act as chemical defence and provide unique pharmaceutical oils. Genome sequencing of the E. grandis sister species E. globulus and a set of inbred E. grandis tree genomes reveals dynamic genome evolution and hotspots of inbreeding depression. The E. grandis genome is the first reference for the eudicot order Myrtales and is placed here sister to the eurosids. This resource expands our understanding of the unique biology of large woody perennials and provides a powerful tool to accelerate comparative biology, breeding and biotechnology.

SND2, a NAC transcription factor gene, regulates genes involved in secondary cell wall development in Arabidopsis fibres and increases fibre cell area in Eucalyptus

BackgroundNAC domain transcription factors initiate secondary cell wall biosynthesis in Arabidopsis fibres and vessels by activating numerous transcriptional regulators and biosynthetic genes. NAC family member SND2 is an indirect target of a principal regulator of fibre secondary cell wall formation, SND1. A previous study showed that overexpression of SND2 produced a fibre cell-specific increase in secondary cell wall thickness in Arabidopsis stems, and that the protein was able to transactivate the cellulose synthase8 (CesA8) promoter. However, the full repertoire of genes regulated by SND2 is unknown, and the effect of its overexpression on cell wall chemistry remains unexplored.ResultsWe overexpressed SND2 in Arabidopsis and analyzed homozygous lines with regards to stem chemistry, biomass and fibre secondary cell wall thickness. A line showing upregulation of CesA8 was selected for transcriptome-wide gene expression profiling. We found evidence for upregulation of biosynthetic genes associated with cellulose, xylan, mannan and lignin polymerization in this line, in agreement with significant co-expression of these genes with native SND2 transcripts according to public microarray repositories. Only minor alterations in cell wall chemistry were detected. Transcription factor MYB103, in addition to SND1, was upregulated in SND2-overexpressing plants, and we detected upregulation of genes encoding components of a signal transduction machinery recently proposed to initiate secondary cell wall formation. Several homozygous T4 and hemizygous T1 transgenic lines with pronounced SND2 overexpression levels revealed a negative impact on fibre wall deposition, which may be indirectly attributable to excessive overexpression rather than co-suppression. Conversely, overexpression of SND2 in Eucalyptus stems led to increased fibre cross-sectional cell area.ConclusionsThis study supports a function for SND2 in the regulation of cellulose and hemicellulose biosynthetic genes in addition of those involved in lignin polymerization and signalling. SND2 seems to occupy a subordinate but central tier in the secondary cell wall transcriptional network. Our results reveal phenotypic differences in the effect of SND2 overexpression between woody and herbaceous stems and emphasize the importance of expression thresholds in transcription factor studies.

Navigating the transcriptional roadmap regulating plant secondary cell wall deposition

The current status of lignocellulosic biomass as an invaluable resource in industry, agriculture, and health has spurred increased interest in understanding the transcriptional regulation of secondary cell wall (SCW) biosynthesis. The last decade of research has revealed an extensive network of NAC, MYB and other families of transcription factors regulating Arabidopsis SCW biosynthesis, and numerous studies have explored SCW-related transcription factors in other dicots and monocots. Whilst the general structure of the Arabidopsis network has been a topic of several reviews, they have not comprehensively represented the detailed protein–DNA and protein–protein interactions described in the literature, and an understanding of network dynamics and functionality has not yet been achieved for SCW formation. Furthermore the methodologies employed in studies of SCW transcriptional regulation have not received much attention, especially in the case of non-model organisms. In this review, we have reconstructed the most exhaustive literature-based network representations to date of SCW transcriptional regulation in Arabidopsis. We include a manipulable Cytoscape representation of the Arabidopsis SCW transcriptional network to aid in future studies, along with a list of supporting literature for each documented interaction. Amongst other topics, we discuss the various components of the network, its evolutionary conservation in plants, putative modules and dynamic mechanisms that may influence network function, and the approaches that have been employed in network inference. Future research should aim to better understand network function and its response to dynamic perturbations, whilst the development and application of genome-wide approaches such as ChIP-seq and systems genetics are in progress for the study of SCW transcriptional regulation in non-model organisms.

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.

Genome-wide mapping of histone H3 lysine 4 trimethylation in Eucalyptus grandis developing xylem

BackgroundHistone modifications play an integral role in plant development, but have been poorly studied in woody plants. Investigating chromatin organization in wood-forming tissue and its role in regulating gene expression allows us to understand the mechanisms underlying cellular differentiation during xylogenesis (wood formation) and identify novel functional regions in plant genomes. However, woody tissue poses unique challenges for using high-throughput chromatin immunoprecipitation (ChIP) techniques for studying genome-wide histone modifications in vivo. We investigated the role of the modified histone H3K4me3 (trimethylated lysine 4 of histone H3) in gene expression during the early stages of wood formation using ChIP-seq in Eucalyptus grandis, a woody biomass model.ResultsPlant chromatin fixation and isolation protocols were optimized for developing xylem tissue collected from field-grown E. grandis trees. A “nano-ChIP-seq” procedure was employed for ChIP DNA amplification. Over 9 million H3K4me3 ChIP-seq and 18 million control paired-end reads were mapped to the E. grandis reference genome for peak-calling using Model-based Analysis of ChIP-Seq. The 12,177 significant H3K4me3 peaks identified covered ~1.5% of the genome and overlapped some 9,623 protein-coding genes and 38 noncoding RNAs. H3K4me3 library coverage, peaking ~600 - 700 bp downstream of the transcription start site, was highly correlated with gene expression levels measured with RNA-seq. Overall, H3K4me3-enriched genes tended to be less tissue-specific than unenriched genes and were overrepresented for general cellular metabolism and development gene ontology terms. Relative expression of H3K4me3-enriched genes in developing secondary xylem was higher than unenriched genes, however, and highly expressed secondary cell wall-related genes were enriched for H3K4me3 as validated using ChIP-qPCR.ConclusionsIn this first genome-wide analysis of a modified histone in a woody tissue, we optimized a ChIP-seq procedure suitable for field-collected samples. In developing E. grandis xylem, H3K4me3 enrichment is an indicator of active transcription, consistent with its known role in sustaining pre-initiation complex formation in yeast. The H3K4me3 ChIP-seq data from this study paves the way to understanding the chromatin landscape and epigenomic architecture of xylogenesis in plants, and complements RNA-seq evidence of gene expression for the future improvement of the E. grandis genome annotation.

Characterising the role of the Eucalyptus grandis SND2 promoter in secondary cell wall biosynthesis

Background NAC and MYB transcription factors (TFs) have been shown to play prominent roles in the regulation of plant developmental processes. Two Arabidopsis thaliana NAC domain TFs (AtSND2, AtSND3) and one MYB domain TF (AtMYB103) were shown to be downstream targets of two master regulators of xylem fibre cell development, NST1 and SND1 [1,2]. These TFs were able to induce the expression of a GUS reporter gene under the control of the AtCesA8 promoter [3], indicating that they may be involved in the regulation of cellulose biosynthesis in the secondary cell walls of A. thaliana xylem fibres. It is hypothesized that putative orthologs of these TFs will also play important roles in regulating fibre secondary cell wall biosynthesis in woody plants such as Eucalyptus grandis. The transcriptional network regulating wood fibre development is uncharacterized in E. grandis, therefore it would be beneficial to identify upstream components of this transcriptional network in E. grandis. In this ongoing study we aim to identify TFs which bind to the promoter of the putative ortholog of AtSND2 in E. grandis (EgSND2) and possibly also to the promoters of the orthologs of AtSND3 and AtMYB103. In parallel, we are characterizing the detailed heterologous expression patterns of the EgSND2, EgSND3 and MYB103 promoter regions in A. thaliana plants. This work forms part of an effort to elucidate the transcriptional network regulating wood fibre development in E. grandis. Methods A reverse BLAST approach was used to identify candidate orthologs of AtSND2, AtSND3 and AtMYB103 in Eucalyptus, Populus and Vitis. In silico cis-element analysis was performed on the promoter regions of the putative orthologs to identify previously characterised and novel cis-elements. The 1.5 kb regions upstream of the translational start site (TSS) of EgSND2, EgSND3 and EgMYB103 were isolated from E. grandis genomic DNA. The amplified fragments were cloned into pMDC162, a GUS reporter vector, and introduced into A. thaliana Col-O plants for heterologous GUS expression analysis. The same reporter constructs were also transformed directly into the vascular cambium of potted E. grandis plants to determine endogenous promoter activity by Induced Somatic Sector Analysis (ISSA, [4]). Qualitative GUS reporter analyses were performed on the EgSND2promoter::GUS, EgSND3promoter::GUS and EgMYB103promoter::GUS constructs in Arabidopsis plants at 1, 3 and 6 weeks. A 500 bp truncation upstream of the TSS of EgSND2 was generated and cloned into the pHIS2.1 vector along with the 1.5 kb EgSND2 promoter sequence for use in yeast onehybrid (Y1-H) screening. Candidate proteins which bind to the cloned promoter sequences will be identified using Y1-H analysis and further characterised.

Integrated analysis and transcript abundance modelling of H3K4me3 and H3K27me3 in developing secondary xylem

Despite the considerable contribution of xylem development (xylogenesis) to plant biomass accumulation, its epigenetic regulation is poorly understood. Furthermore, the relative contributions of histone modifications to transcriptional regulation is not well studied in plants. We investigated the biological relevance of H3K4me3 and H3K27me3 in secondary xylem development using ChIP-seq and their association with transcript levels among other histone modifications in woody and herbaceous models. In developing secondary xylem of the woody model Eucalyptus grandis, H3K4me3 and H3K27me3 genomic spans were distinctly associated with xylogenesis-related processes, with (late) lignification pathways enriched for putative bivalent domains, but not early secondary cell wall polysaccharide deposition. H3K27me3-occupied genes, of which 753 (~31%) are novel targets, were enriched for transcriptional regulation and flower development and had significant preferential expression in roots. Linear regression models of the ChIP-seq profiles predicted ~50% of transcript abundance measured with strand-specific RNA-seq, confirmed in a parallel analysis in Arabidopsis where integration of seven additional histone modifications each contributed smaller proportions of unique information to the predictive models. This study uncovers the biological importance of histone modification antagonism and genomic span in xylogenesis and quantifies for the first time the relative correlations of histone modifications with transcript abundance in plants.

The role of SND2 in the regulation of Arabidopsis fibre secondary cell wall formation

Background Transcription factors (TFs) play important roles in the regulation of secondary cell wall (SCW) biosynthesis in herbaceous and woody plants. In Arabidopsis, the onset of SCW deposition is initiated by a nexus of NAC, MYB, homeodomain and several other families of TFs, which function in a transcriptional network regulating SCW biosynthetic genes. NAC family members SND1/ NST1 and VND6/VND7 have been identified as functionally redundant master regulators of SCW formation in fibres and vessels, respectively [1,2]. Arabidopsis plants overexpressing SND2, an indirect target of fibre master regulator SND1, exhibited increased SCW thickness in inflorescence stem fibres, whilst dominant repression lines exhibited a decrease in fibre SCW thickness associated with a reduction in glucose and xylose cell wall sugar content [3]. The ability of SND2 to transactivate the CesA8 promoter [3] suggested that SND2 may regulate cellulose biosynthetic genes during fibre SCW formation. The evaluation of this hypothesis necessitates the identification of all downstream genes potentially regulated by SND2, and the analysis of SCW chemistry and morphology in overexpression lines. The aim of this ongoing study is to further elucidate the role of SND2 in fibre SCW formation through microarray analysis of overexpression lines and by independent confirmation of the effect of SND2 overexpression on fibre SCW chemistry and morphology in Arabidopsis plants.

Temporal analysis of Arabidopsis genes activated by Eucalyptus grandis NAC transcription factors associated with xylem fibre and vessel development

Secondary cell wall (SCW) deposition in Arabidopsis is regulated among others by NAC transcription factors, where SND1 chiefly initiates xylem fibre differentiation while VND6 controls metaxylem vessel SCW development, especially programmed cell death and wall patterning. The translational relevance of Arabidopsis SCW regulation theory and the utility of characterized transcription factors as modular synthetic biology tools for improving commercial fibre crops is unclear. We investigated inter-lineage gene activation dynamics for potential fibre and vessel differentiation regulators from the widely grown hardwood Eucalyptus grandis (Myrtales). EgrNAC26, a VND6 homolog, and EgrNAC61, an SND1 homolog, were transiently expressed in Arabidopsis mesophyll protoplasts in parallel to determine early and late (i.e. 7 and 14 hours post-transfection) gene targets. Surprisingly, across the time series EgrNAC26 activated only a subset of SCW-related transcription factors and biosynthetic genes activated by EgrNAC61, specializing instead in targeting vessel-specific wall pit and programmed cell death markers. Promoters of EgrNAC26 and EgrNAC61 both induced reporter gene expression in vessels of young Arabidopsis plants, with EgrNAC61 also conferring xylem- and cork cambium-preferential expression in Populus. Our results demonstrate partial conservation, with notable exceptions, of SND1 and VND6 homologs in Eucalyptus and a first report of cork cambium expression for EgrNAC61.

Genome-wide analysis of the lateral organ boundaries domain (LBD) gene family in Eucalyptus grandis reveals members that differentially impact secondary growth.

Summary Lateral Organ Boundaries Domain (LBD) proteins are plant‐specific transcription factors playing crucial roles in growth and development. However, the function of LBD proteins in Eucalyptus grandis remains largely unexplored. In this study, LBD genes in E. grandis were identified and characterized using bioinformatics approaches. Gene expression patterns in various tissues and the transcriptional responses of EgLBDs to exogenous hormones were determined by qRT‐PCR. Functions of the selected EgLBDs were studied by ectopically overexpressing in a hybrid poplar (Populus alba × Populus glandulosa). Expression levels of genes in the transgenic plants were investigated by RNA‐seq. Our results showed that there were forty‐six EgLBD members in the E. grandis genome and three EgLBDs displayed xylem‐ (EgLBD29) or phloem‐preferential expression (EgLBD22 and EgLBD37). Confocal microscopy indicated that EgLBD22, EgLBD 29 and EgLBD37 were localized to the nucleus. Furthermore, we found that EgLBD22, EgLBD29 and EgLBD37 were responsive to the treatments of indol‐3‐acetic acid and gibberellic acid. More importantly, we demonstrated EgLBDs exerted different influences on secondary growth. Namely, 35S::EgLBD37 led to significantly increased secondary xylem, 35S::EgLBD29 led to greatly increased phloem fibre production, and 35S::EgLBD22 showed no obvious effects. We revealed that key genes related to gibberellin, ethylene and auxin signalling pathway as well as cell expansion were significantly up‐ or down‐regulated in transgenic plants. Our new findings suggest that LBD genes in E. grandis play important roles in secondary growth. This provides new mechanisms to increase wood or fibre production.