Search

A multi-tiered transcriptional network regulates xylem differentiation and secondary cell wall (SCW) formation in plants, with evidence of both conserved and lineage-specific SCW network architecture. We aimed to elucidate the roles of selected R2R3-MYB transcription factors (TFs) linked to Eucalyptus wood formation by identifying genome-wide TF binding sites and direct target genes through an improved DAP-seq protocol combined with machine learning for target gene assignment (DAP-seq-ML). We applied this to five TFs including a well-studied SCW master regulator (EgrMYB2; homolog of AtMYB83), a repressor of lignification (EgrMYB1; homolog of AtMYB4), a TF affecting SCW thickness and vessel density (EgrMYB137; homolog of PtrMYB074) and two TFs with unclear roles in SCW regulation (EgrMYB135 and EgrMYB122). Each DAP-seq TF peak set (average 12,613 peaks) was enriched for canonical R2R3-MYB binding motifs. To improve the reliability of target gene assignment to peaks, a random forest classifier was developed from Arabidopsis DAP-seq, RNA-seq, chromatin, and conserved noncoding sequence data which demonstrated significantly higher precision and recall to the baseline method of assigning genes to proximal peaks. EgrMYB1, EgrMYB2 and EgrMYB137 predicted targets showed clear enrichment for SCW-related biological processes. As validation, EgrMYB137 overexpression in transgenic Eucalyptus hairy roots increased xylem lignification, while its dominant repression in transgenic Arabidopsis and Populus reduced xylem lignification, stunted growth, and caused downregulation of SCW genes. EgrMYB137 targets overlapped significantly with those of EgrMYB2, suggesting partial functional redundancy. Our results show that DAP-seq-ML identified biologically relevant R2R3-MYB targets supported by the finding that EgrMYB137 promotes SCW lignification in planta., (© 2023. The Author(s), under exclusive licence to Springer Nature B.V.)

To improve our understanding of genetic mechanisms underlying complex traits in plants, a comprehensive analysis of gene variants is required. Eucalyptus is an important forest plantation genus that is highly outbred. Trait dissection and molecular breeding in eucalypts currently relies on biallelic single-nucleotide polymorphism (SNP) markers. These markers fail to capture the large amount of haplotype diversity in these species, and thus multi-allelic markers are required. We aimed to develop a gene-based haplotype mining panel for Eucalyptus species. We generated 17 999 oligonucleotide probe sets for targeted sequencing of selected regions of 6293 genes implicated in growth and wood properties, pest and disease resistance, and abiotic stress responses. We identified and phased 195 834 SNPs using a read-based phasing approach to reveal SNP-based haplotypes. A total of 8915 target regions (at 4637 gene loci) passed tests for Mendelian inheritance. We evaluated the haplotype panel in four Eucalyptus species (E. grandis, E. urophylla, E. dunnii and E. nitens) to determine its ability to capture diversity across eucalypt species. This revealed an average of 3.13-4.52 haplotypes per target region in each species, and 33.36% of the identified haplotypes were shared by at least two species. This haplotype mining panel will enable the analysis of haplotype diversity within and between species, and provide multi-allelic markers that can be used for genome-wide association studies and gene-based breeding approaches., (© 2022 The Authors. The Plant Journal published by Society for Experimental Biology and John Wiley & Sons Ltd.)

Background: De novo phased (haplo)genome assembly using long-read DNA sequencing data has improved the detection and characterization of structural variants (SVs) in plant and animal genomes. Able to span across haplotypes, long reads allow phased, haplogenome assembly in highly outbred organisms such as forest trees. Eucalyptus tree species and interspecific hybrids are the most widely planted hardwood trees with F1 hybrids of Eucalyptus grandis and E. urophylla forming the bulk of fast-growing pulpwood plantations in subtropical regions. The extent of structural variation and its effect on interspecific hybridization is unknown in these trees. As a first step towards elucidating the extent of structural variation between the genomes of E. grandis and E. urophylla, we sequenced and assembled the haplogenomes contained in an F1 hybrid of the two species., Findings: Using Nanopore sequencing and a trio-binning approach, we assembled the separate haplogenomes (566.7 Mb and 544.5 Mb) to 98.0% BUSCO completion. High-density SNP genetic linkage maps of both parents allowed scaffolding of 88.0% of the haplogenome contigs into 11 pseudo-chromosomes (scaffold N50 of 43.8 Mb and 42.5 Mb for the E. grandis and E. urophylla haplogenomes, respectively). We identify 48,729 SVs between the two haplogenomes providing the first detailed insight into genome structural rearrangement in these species. The two haplogenomes have similar gene content, 35,572 and 33,915 functionally annotated genes, of which 34.7% are contained in genome rearrangements., Conclusions: Knowledge of SV and haplotype diversity in the two species will form the basis for understanding the genetic basis of hybrid superiority in these trees., (© The Author(s) 2023. Published by Oxford University Press GigaScience.)

From its origins in Australia, Eucalyptus grandis has spread to every continent, except Antarctica, as a wood crop. It has been cultivated and bred for over 100 yr in places such as South Africa. Unlike most annual crops and fruit trees, domestication of E. grandis is still in its infancy, representing a unique opportunity to interrogate the genomic consequences of artificial selection early in the domestication process. To determine how a century of artificial selection has changed the genome of E. grandis, we generated single nucleotide polymorphism genotypes for 1080 individuals from three advanced South African breeding programmes using the EUChip60K chip, and investigated population structure and genome-wide differentiation patterns relative to wild progenitors. Breeding and wild populations appeared genetically distinct. We found genomic evidence of evolutionary processes known to have occurred in other plant domesticates, including interspecific introgression and intraspecific infusion from wild material. Furthermore, we found genomic regions with increased linkage disequilibrium and genetic differentiation, putatively representing early soft sweeps of selection. This is, to our knowledge, the first study of genomic signatures of domestication in a timber species looking beyond the first few generations of cultivation. Our findings highlight the importance of intra- and interspecific hybridization during early domestication., (© 2022 The Authors. New Phytologist © 2022 New Phytologist Foundation.)

Terpenes are an important group of plant specialized metabolites influencing, amongst other functions, defence mechanisms against pests. We used a genome-wide association study to identify single nucleotide polymorphism (SNP) markers and putative candidate genes for terpene traits. We tested 15,387 informative SNP markers derived from genotyping 416 Eucalyptus grandis individuals for association with 3 terpene traits, 1,8-cineole, γ-terpinene, and p-cymene. A multilocus mixed model analysis identified 21 SNP markers for 1,8-cineole on chromosomes 2, 4, 6, 7, 8, 9, 10, and 11, that individually explained 3.0%-8.4% and jointly 42.7% of the phenotypic variation. Association analysis of γ-terpinene found 32 significant SNP markers on chromosomes 1, 2, 4, 5, 6, 9, and 11, explaining 3.4-15.5% and jointly 54.5% of phenotypic variation. For p-cymene, 28 significant SNP markers were identified on chromosomes 1, 2, 3, 5, 6, 7, 10, and 11, explaining 3.4-16.1% of the phenotypic variation and jointly 46.9%. Our results show that variation underlying the 3 terpene traits is influenced by a few minor loci in combination with a few major effect loci, suggesting an oligogenic nature of the traits., (© The Author(s) 2022. Published by Oxford University Press on behalf of Genetics Society of America.)

We performed gene and genome targeted SNP discovery towards the development of a genome-wide, multispecies genotyping array for tropical pines. Pooled RNA-seq data from shoots of seedlings from five tropical pine species was used to identify transcript-based SNPs resulting in 1.3 million candidate Affymetrix SNP probe sets. In addition, we used a custom 40 K probe set to perform capture-seq in pooled DNA from 81 provenances representing the natural ranges of six tropical pine species in Mexico and Central America resulting in 563 K candidate SNP probe sets. Altogether, 300 K RNA-seq (72%) and 120 K capture-seq (28%) derived SNP probe sets were tiled on a 420 K screening array that was used to genotype 576 trees representing the 81 provenances and commercial breeding material. Based on the screening array results, 50 K SNPs were selected for commercial SNP array production including 20 K polymorphic SNPs for P. patula, P. tecunumanii, P. oocarpa and P. caribaea, 15 K for P. greggii and P. maximinoi, 13 K for P. elliottii and 8K for P. pseudostrobus. We included 9.7 K ancestry informative SNPs that will be valuable for species and hybrid discrimination. Of the 50 K SNP markers, 25% are polymorphic in only one species, while 75% are shared by two or more species. The Pitro50K SNP chip will be useful for population genomics and molecular breeding in this group of pine species that, together with their hybrids, represent the majority of fast-growing tropical and subtropical pine plantations globally., (© 2021 John Wiley & Sons Ltd.)

Background: Affordable high-throughput DNA and RNA sequencing technologies are allowing genomic analysis of plant and animal populations and as a result empowering new systems genetics approaches to study complex traits. The availability of intuitive tools to browse and analyze the resulting large-scale genetic and genomic datasets remain a significant challenge. Furthermore, these integrative genomics approaches require innovative methods to dissect the flow and interconnectedness of biological information underlying complex trait variation. The Plant Genome Integrative Explorer (PlantGenIE.org) is a multi-species database and domain that houses online tools for model and woody plant species including Eucalyptus. Since the Eucalyptus Genome Integrative Explorer (EucGenIE) is integrated within PlantGenIE, it shares genome and expression analysis tools previously implemented within the various subdomains (ConGenIE, PopGenIE and AtGenIE). Despite the success in setting up integrative genomics databases, online tools for systems genetics modelling and high-resolution dissection of complex trait variation in plant populations have been lacking., Results: We have developed qtlXplorer ( https://eucgenie.org/QTLXplorer ) for visualizing and exploring systems genetics data from genome-wide association studies including quantitative trait loci (QTLs) and expression-based QTL (eQTL) associations. This module allows users to, for example, find co-located QTLs and eQTLs using an interactive version of Circos, or explore underlying genes using JBrowse. It provides users with a means to build systems genetics models and generate hypotheses from large-scale population genomics data. We also substantially upgraded the EucGenIE resource and show how it enables users to combine genomics and systems genetics approaches to discover candidate genes involved in biotic stress responses and wood formation by focusing on two multigene families, laccases and peroxidases., Conclusions: qtlXplorer adds a new dimension, population genomics, to the EucGenIE and PlantGenIE environment. The resource will be of interest to researchers and molecular breeders working in Eucalyptus and other woody plant species. It provides an example of how systems genetics data can be integrated with functional genetics data to provide biological insight and formulate hypotheses. Importantly, integration within PlantGenIE enables novel comparative genomics analyses to be performed from population-scale data., (© 2021. The Author(s).)

Eucalyptus is among the most widely planted taxa of forest trees worldwide. However, its spread as an exotic or genetically engineered form can create ecological and social problems. To mitigate gene flow via pollen and seeds, we mutated the Eucalyptus orthologue of LEAFY (LFY) by transforming a Eucalyptus grandis × urophylla wild-type hybrid and two Flowering Locus T (FT) overexpressing (and flowering) lines with CRISPR Cas9 targeting its LFY orthologue, ELFY. We achieved high rates of elfy biallelic knockouts, often approaching 100% of transgene insertion events. Frameshift mutations and deletions removing conserved amino acids caused strong floral alterations, including indeterminacy in floral development and an absence of male and female gametes. These mutants were otherwise visibly normal and did not differ statistically from transgenic controls in juvenile vegetative growth rate or leaf morphology in greenhouse trials. Genes upstream or near to ELFY in the floral development pathway were overexpressed, whereas floral organ identity genes downstream of ELFY were severely depressed. We conclude that disruption of ELFY function appears to be a useful tool for sexual containment, without causing statistically significant or large adverse effects on juvenile vegetative growth or leaf morphology., (© 2021 The Authors. Plant Biotechnology Journal published by Society for Experimental Biology and The Association of Applied Biologists and John Wiley & Sons Ltd.)

Eucalyptus grandis is one of the most important species for hardwood plantation forestry around the world. At present, its commercial deployment is in decline because of pests and pathogens such as Leptocybe invasa gall wasp ( Lepto ), and often co-occurring fungal stem diseases such as Botryosphaeria dothidea and Teratosphaeria zuluensis ( BotryoTera ). This study analyzed Lepto , BotryoTera , and stem diameter growth in an E. grandis multi-environmental, genetic trial. The study was established in three subtropical environments. Diameter growth and BotryoTera incidence scores were assessed on 3,334 trees, and Lepto incidence was assessed on 4,463 trees from 95 half-sib families. Using the Eucalyptus EUChip60K SNP chip, a subset of 964 trees from 93 half-sib families were genotyped with 14,347 informative SNP markers. We employed single-step genomic BLUP (ssGBLUP) to estimate genetic parameters in the genetic trial. Diameter and Lepto tolerance showed a positive genetic correlation (0.78), while BotryoTera tolerance had a negative genetic correlation with diameter growth (-0.38). The expected genetic gains for diameter growth and Lepto and BotryoTera tolerance were 12.4, 10, and -3.4%, respectively. We propose a genomic selection breeding strategy for E. grandis that addresses some of the present population structure problems., Competing Interests: MM was employed by the company Mondi South Africa (Pty) Ltd. The authors declare that this study received funding from Mondi South Africa (Pty) Ltd. The funder was not involved in the study design, collection, analysis, and interpretation of data, the writing of this article, or the decision to submit it for publication., (Copyright © 2021 Mphahlele, Isik, Hodge and Myburg.)

The genetic consequences of adaptation to changing environments can be deciphered using population genomics, which may help predict species' responses to global climate change. Towards this, we used genome-wide SNP marker analysis to determine population structure and patterns of genetic differentiation in terms of neutral and adaptive genetic variation in the natural range of Eucalyptus grandis, a widely cultivated subtropical and temperate species, serving as genomic reference for the genus. We analysed introgression patterns at subchromosomal resolution using a modified ancestry mapping approach and identified provenances with extensive interspecific introgression in response to increased aridity. Furthermore, we describe potentially adaptive genetic variation as explained by environment-associated SNP markers, which also led to the discovery of what is likely a large structural variant. Finally, we show that genes linked to these markers are enriched for biotic and abiotic stress responses., (© 2020 John Wiley & Sons Ltd.)

Gall-inducing insects and their hosts present some of the most intricate plant-herbivore interactions. Oviposition on the host is often the first cue of future herbivory and events at this early time point can affect later life stages. Many gallers are devastating plant pests, yet little information regarding the plant-insect molecular interplay exists, particularly following egg deposition. We studied the physiological and transcriptional responses of Eucalyptus following oviposition by the gall-inducing wasp, Leptocybe invasa, to explore potential mechanisms governing defence responses and gall development. RNA sequencing and microscopy were used to explore a susceptible Eucalyptus-L. invasa interaction. Infested and control material was compared over time (1-3, 7 and 90 days post oviposition) to examine the transcriptional and morphological changes. Oviposition induces accumulation of reactive oxygen species and phenolics which is reflected in the transcriptome analysis. Gene expression supports phytohormones and 10 transcription factor subfamilies as key regulators. The egg and oviposition fluid stimulate cell division resulting in gall development. Eucalyptus responses to oviposition are apparent within 24 hr. Putative defences include the oxidative burst and barrier reinforcement. However, egg and oviposition fluid stimuli may redirect these responses towards gall development., (© 2020 John Wiley & Sons Ltd.)

Human life intimately depends on plants for food, biomaterials, health, energy, and a sustainable environment. Various plants have been genetically improved mostly through breeding, along with limited modification via genetic engineering, yet they are still not able to meet the ever-increasing needs, in terms of both quantity and quality, resulting from the rapid increase in world population and expected standards of living. A step change that may address these challenges would be to expand the potential of plants using biosystems design approaches. This represents a shift in plant science research from relatively simple trial-and-error approaches to innovative strategies based on predictive models of biological systems. Plant biosystems design seeks to accelerate plant genetic improvement using genome editing and genetic circuit engineering or create novel plant systems through de novo synthesis of plant genomes. From this perspective, we present a comprehensive roadmap of plant biosystems design covering theories, principles, and technical methods, along with potential applications in basic and applied plant biology research. We highlight current challenges, future opportunities, and research priorities, along with a framework for international collaboration, towards rapid advancement of this emerging interdisciplinary area of research. Finally, we discuss the importance of social responsibility in utilizing plant biosystems design and suggest strategies for improving public perception, trust, and acceptance., Competing Interests: The authors declare that they have no conflicts of interest regarding the publication of this article., (Copyright © 2020 Xiaohan Yang et al.)

The galling insect, Leptocybe invasa, causes significant losses in plantations of various Eucalyptus species and hybrids, threatening its economic viability. We applied a genome-wide association study (GWAS) to identify single-nucleotide polymorphism (SNP) markers associated with resistance to L. invasa. A total of 563 insect-challenged Eucalyptus grandis trees, from 61 half-sib families, were genotyped using the EUChip60K SNP chip, and we identified 15,445 informative SNP markers in the test population. Multi-locus mixed-model (MLMM) analysis identified 35 SNP markers putatively associated with resistance to L. invasa based on four discreet classes of insect damage scores: (0) not infested, (1) infested showing evidence of oviposition but no gall development, (2) infested with galls on leaves, midribs or petioles and (3) stunting and lethal gall formation. MLMM analysis identified three associated genomic regions on chromosomes 3, 7 and 8 jointly explaining 17.6% of the total phenotypic variation. SNP analysis of a validation population of 494 E. grandis trees confirmed seven SNP markers that were also detected in the initial association analysis. Based on transcriptome profiles of resistant and susceptible genotypes from an independent experiment, we identified several putative candidate genes in associated genomic loci including Nucleotide-binding ARC- domain (NB-ARC) and toll-interleukin-1-receptor-Nucleotide binding signal- Leucine rich repeat (TIR-NBS-LRR) genes. Our results suggest that Leptocybe resistance in E. grandis may be influenced by a few large-effect loci in combination with minor effect loci segregating in our test and validation populations., (© The Author(s) 2020. Published by Oxford University Press on behalf of Japanese Society of Plant Physiologists. All rights reserved. For permissions, please email: journals.permissions@oup.com.)

SECONDARY WALL-ASSOCIATED NAC DOMAIN1 ( SND1 ) is a master regulator of fibre secondary wall deposition in Arabidopsis thaliana (Arabidopsis), with homologs in other angiosperms and gymnosperms. However, it is poorly understood to what extent the fibre-specific regulation of the SND1 promoter, and that of its orthologs, is conserved between diverged herbaceous and woody lineages. We performed a reciprocal reporter gene analysis of orthologous SND1 promoters from Arabidopsis ( AthSND1 ), Eucalyptus grandis ( EgrNAC61 ) and Populus alba × P. grandidentata ( PagWND1A ) relative to secondary cell wall-specific Cellulose Synthase4 (CesA4 ) and CesA7 promoters, in both a non-woody (Arabidopsis) and a woody (poplar) system. β-glucuronidase (GUS) reporter analysis in Arabidopsis showed that the SND1 promoter was active in vascular tissues as previously reported and showed interfascicular and xylary fibre-specific expression in inflorescence stems, while reporter constructs of the woody plant-derived promoters were partial to the (pro)cambium-phloem and protoxylem. In transgenic P. tremula × P. alba plants, all three orthologous SND1 promoters expressed the GUS reporter similarly and preferentially in developing secondary xylem, ray parenchyma and cork cambium. Ours is the first study to reciprocally test orthologous SND1 promoter specificity in herbaceous and woody species, revealing diverged regulatory functions in the herbaceous system.

Fusarium circinatum poses a serious threat to many pine species in both commercial and natural pine forests. Knowledge regarding the molecular basis of pine- F. circinatum host-pathogen interactions could assist efforts to produce more resistant planting stock. This study aimed to identify molecular responses underlying resistance against F. circinatum . A dual RNA-seq approach was used to investigate host and pathogen expression in F. circinatum challenged Pinus tecunumanii (resistant) and Pinus patula (susceptible), at three- and seven-days post inoculation. RNA-seq reads were mapped to combined host-pathogen references for both pine species to identify differentially expressed genes (DEGs). F. circinatum genes expressed during infection showed decreased ergosterol biosynthesis in P. tecunumanii relative to P. patula . For P. tecunumanii , enriched gene ontologies and DEGs indicated roles for auxin-, ethylene-, jasmonate- and salicylate-mediated phytohormone signalling. Correspondingly, key phytohormone signaling components were down-regulated in P. patula . Key F. circinatum ergosterol biosynthesis genes were expressed at lower levels during infection of the resistant relative to the susceptible host. This study further suggests that coordination of phytohormone signaling is required for F. circinatum resistance in P. tecunumanii , while a comparatively delayed response and impaired phytohormone signaling contributes to susceptibility in P. patula .

Accessible chromatin changes dynamically during development and harbours functional regulatory regions which are poorly understood in the context of wood development. We explored the importance of accessible chromatin in Eucalyptus grandis in immature xylem generally, and MYB transcription factor-mediated transcriptional programmes specifically. We identified biologically reproducible DNase I Hypersensitive Sites (DHSs) and assessed their functional significance in immature xylem through their associations with gene expression, epigenomic data and DNA sequence conservation. We identified in vitro DNA binding sites for six secondary cell wall-associated Eucalyptus MYB (EgrMYB) transcription factors using DAP-seq, reconstructed protein-DNA networks of predicted targets based on binding sites within or outside DHSs and assessed biological enrichment of these networks with published datasets. 25 319 identified immature xylem DHSs were associated with increased transcription and significantly enriched for various epigenetic signatures (H3K4me3, H3K27me3, RNA pol II), conserved noncoding sequences and depleted single nucleotide variants. Predicted networks built from EgrMYB binding sites located in accessible chromatin were significantly enriched for systems biology datasets relevant to wood formation, whereas those occurring in inaccessible chromatin were not. Our study demonstrates that DHSs in E. grandis immature xylem, most of which are intergenic, are of functional significance to gene regulation in this tissue., (© 2019 University of Pretoria New Phytologist © 2019 New Phytologist Trust.)

Acetyl- and methylglucuronic acid decorations of xylan, the dominant hemicellulose in secondary cell walls (SCWs) of woody dicots, affect its interaction with cellulose and lignin to determine SCW structure and extractability. Genes and pathways involved in these modifications may be targets for genetic engineering; however, little is known about the regulation of xylan modifications in woody plants. To address this, we assessed genetic and gene expression variation associated with xylan modification in developing xylem of Eucalyptus grandis × Eucalyptus urophylla interspecific hybrids. Expression quantitative trait locus (eQTL) mapping identified potential regulatory polymorphisms affecting gene expression modules associated with xylan modification. We identified 14 putative xylan modification genes that are members of five expression modules sharing seven trans-eQTL hotspots. The xylan modification genes are prevalent in two expression modules. The first comprises nucleotide sugar interconversion pathways supplying the essential precursors for cellulose and xylan biosynthesis. The second contains genes responsible for phenylalanine biosynthesis and S-adenosylmethionine biosynthesis required for glucuronic acid and monolignol methylation. Co-expression and co-regulation analyses also identified four metabolic sources of acetyl coenxyme A that appear to be transcriptionally coordinated with xylan modification. Our systems genetics analysis may provide new avenues for metabolic engineering to alter wood SCW biology for enhanced biomass processability., (© 2019 The Authors. New Phytologist © 2019 New Phytologist Trust.)

Fast-growing forest plantations are sustainable feedstocks of plant biomass that can serve as alternatives to fossil carbon resources for materials, chemicals, and energy. Their ability to efficiently harvest light energy and carbon from the atmosphere and sequester this into metabolic precursors for lignocellulosic biopolymers and a wide range of plant specialized metabolites make them excellent biochemical production platforms and living biorefineries. Their large sizes have facilitated multi-omics analyses and systems modeling of key biological processes such as lignin biosynthesis in trees. High-throughput 'omics' approaches have also been applied in segregating tree populations where genetic variation creates abundant genetic perturbations of system components allowing construction of systems genetics models linking genes and pathways to complex trait variation. With this information in hand, it is now possible to start using synthetic biology and genome editing techniques in a bioengineering approach based on a deeper understanding and rational design of biological parts, devices, and integrated systems. However, the complexity of the biology and interacting components will require investment in big data informatics, machine learning, and intuitive visualization to fully explore multi-dimensional patterns and identify emergent properties of biological systems. Predictive systems models could be tested rapidly through high-throughput synthetic biology approaches and multigene editing. Such a bioengineering paradigm, together with accelerated genomic breeding, will be crucial for the development of a new generation of woody biorefinery crops.

Subcellular compartmentation of plant biosynthetic pathways in the mitochondria and plastids requires coordinated regulation of nuclear encoded genes, and the role of these genes has been largely ignored by wood researchers. In this study, we constructed a targeted systems genetics coexpression network of xylogenesis in Eucalyptus using plastid and mitochondrial carbon metabolic genes and compared the resulting clusters to the aspen xylem developmental series. The constructed network clusters reveal the organization of transcriptional modules regulating subcellular metabolic functions in plastids and mitochondria. Overlapping genes between the plastid and mitochondrial networks implicate the common transcriptional regulation of carbon metabolism during xylem secondary growth. We show that the central processes of organellar carbon metabolism are distinctly coordinated across the developmental stages of wood formation and are specifically associated with primary growth and secondary cell wall deposition. We also demonstrate that, during xylogenesis, plastid-targeted carbon metabolism is partially regulated by the central clock for carbon allocation towards primary and secondary xylem growth, and we discuss these networks in the context of previously established associations with wood-related complex traits. This study provides a new resolution into the integration and transcriptional regulation of plastid- and mitochondrial-localized carbon metabolism during xylogenesis., (© 2019 University of Pretoria New Phytologist © 2019 New Phytologist Trust.)

Lignocellulosic biomass, encompassing cellulose, lignin and hemicellulose in plant secondary cell walls (SCWs), is the most abundant source of renewable materials on earth. Currently, fast-growing woody dicots such as Eucalyptus and Populus trees are major lignocellulosic (wood fiber) feedstocks for bioproducts such as pulp, paper, cellulose, textiles, bioplastics and other biomaterials. Processing wood for these products entails separating the biomass into its three main components as efficiently as possible without compromising yield. Glucuronoxylan (xylan), the main hemicellulose present in the SCWs of hardwood trees carries chemical modifications that are associated with SCW composition and ultrastructure, and affect the recalcitrance of woody biomass to industrial processing. In this review we highlight the importance of xylan properties for industrial wood fiber processing and how gaining a greater understanding of xylan biosynthesis, specifically xylan modification, could yield novel biotechnology approaches to reduce recalcitrance or introduce novel processing traits. Altering xylan modification patterns has recently become a focus of plant SCW studies due to early findings that altered modification patterns can yield beneficial biomass processing traits. Additionally, it has been noted that plants with altered xylan composition display metabolic differences linked to changes in precursor usage. We explore the possibility of using systems biology and systems genetics approaches to gain insight into the coordination of SCW formation with other interdependent biological processes. Acetyl-CoA, s-adenosylmethionine and nucleotide sugars are precursors needed for xylan modification, however, the pathways which produce metabolic pools during different stages of fiber cell wall formation still have to be identified and their co-regulation during SCW formation elucidated. The crucial dependence on precursor metabolism provides an opportunity to alter xylan modification patterns through metabolic engineering of one or more of these interdependent pathways. The complexity of xylan biosynthesis and modification is currently a stumbling point, but it may provide new avenues for woody biomass engineering that are not possible for other biopolymers.

Re-engineering of transcriptional networks regulating secondary cell wall formation may allow the improvement of plant biomass in widely grown plantation crops such as Eucalyptus. However, there is currently a scarcity of freely available standardized biological parts (e.g., Phytobricks) compatible with Type IIS assembly approaches from forest trees, and there is a need to accelerate transcriptional network inference in nonmodel biomass crops. Here we describe the design and synthesis of a versatile three-panel biological parts collection of 221 secondary cell wall-related Eucalyptus grandis transcription factor coding sequences and 65 promoters that are compatible with GATEWAY, Golden Gate, MoClo, and GoldenBraid DNA assembly methods and generally conform to accepted Phytobrick syntaxes. This freely available resource is intended to accelerate synthetic biology applications in multiple plant biomass crops and enable reconstruction of secondary cell wall transcriptional networks using high-throughput assays such as DNA affinity purification sequencing (DAP-seq) and enhanced yeast one-hybrid (eY1H) screening.

Background: Land plant organellar genomes have significant impact on metabolism and adaptation, and as such, accurate assembly and annotation of plant organellar genomes is an important tool in understanding the evolutionary history and interactions between these genomes. Intracellular DNA transfer is ongoing between the nuclear and organellar genomes, and can lead to significant genomic variation between, and within, species that impacts downstream analysis of genomes and transcriptomes., Results: In order to facilitate further studies of cytonuclear interactions in Eucalyptus, we report an updated annotation of the E. grandis plastid genome, and the second sequenced and annotated mitochondrial genome of the Myrtales, that of E. grandis. The 478,813 bp mitochondrial genome shows the conserved protein coding regions and gene order rearrangements typical of land plants. There have been widespread insertions of organellar DNA into the E. grandis nuclear genome, which span 141 annotated nuclear genes. Further, we identify predicted editing sites to allow for the discrimination of RNA-sequencing reads between nuclear and organellar gene copies, finding that nuclear copies of organellar genes are not expressed in E. grandis., Conclusions: The implications of organellar DNA transfer to the nucleus are often ignored, despite the insight they can give into the ongoing evolution of plant genomes, and the problems they can cause in many applications of genomics. Future comparisons of the transcription and regulation of organellar genes between Eucalyptus genotypes may provide insight to the cytonuclear interactions that impact economically important traits in this widely grown lignocellulosic crop species.

Background: Fusarium circinatum is a pressing threat to the cultivation of many economically important pine tree species. Efforts to develop effective disease management strategies can be aided by investigating the molecular mechanisms involved in the host-pathogen interaction between F. circinatum and pine species. Pinus tecunumanii and Pinus patula are two closely related tropical pine species that differ widely in their resistance to F. circinatum challenge, being resistant and susceptible respectively, providing the potential for a useful pathosystem to investigate the molecular responses underlying resistance to F. circinatum. However, no genomic resources are available for P. tecunumanii. Pathogenesis-related proteins are classes of proteins that play important roles in plant-microbe interactions, e.g. chitinases; proteins that break down the major structural component of fungal cell walls. Generating a reference sequence for P. tecunumanii and characterizing pathogenesis related gene families in these two pine species is an important step towards unravelling the pine-F. circinatum interaction., Results: Eight reference based and 12 de novo assembled transcriptomes were produced, for juvenile shoot tissue from both species. EvidentialGene pipeline redundancy reduction, expression filtering, protein clustering and taxonomic filtering produced a 50 Mb shoot transcriptome consisting of 28,621 contigs for P. tecunumanii and a 72 Mb shoot transcriptome consisting of 52,735 contigs for P. patula. Predicted protein sequences encoded by the assembled transcriptomes were clustered with reference proteomes from 92 other species to identify pathogenesis related gene families in P. patula, P. tecunumanii and other pine species., Conclusions: The P. tecunumanii transcriptome is the first gene catalogue for the species, representing an important resource for studying resistance to the pitch canker pathogen, F. circinatum. This study also constitutes, to our knowledge, the largest index of gymnosperm PR-genes to date.

Leptocybe invasa is an insect pest causing gall formation on oviposited shoot tips and leaves of Eucalyptus trees leading to leaf deformation, stunting, and death in severe cases. We previously observed different constitutive and induced terpenes, plant specialized metabolites that may act as attractants or repellents to insects, in a resistant and susceptible clone of Eucalyptus challenged with L. invasa. We tested the hypothesis that specific terpenes are associated with pest resistance in a Eucalyptus grandis half-sib population. Insect damage was scored over 2 infestation cycles, and leaves were harvested for near-infrared reflectance (NIR) and terpene measurements. We used Bayesian model averaging for terpene selection and obtained partial least squares NIR models to predict terpene content and L. invasa infestation damage. In our optimal model, 29% of the phenotypic variation could be explained by 7 terpenes, and the monoterpene combination, limonene, α-terpineol, and 1,8-cineole, could be predicted with an NIR prediction ability of  .67. Bayesian model averaging supported α-pinene, γ-terpinene, and iso-pinocarveol as important for predicting L. invasa infestation. Susceptibility was associated with increased γ-terpinene and α-pinene, which may act as a pest attractant, whereas reduced susceptibility was associated with iso-pinocarveol, which may act to recruit parasitoids or have direct toxic effects., (© 2018 John Wiley & Sons Ltd.)

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.

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, EgLBD29 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., (© 2017 The Authors. Plant Biotechnology Journal published by Society for Experimental Biology and The Association of Applied Biologists and John Wiley & Sons Ltd.)

Lignocellulosic biomass is a promising feedstock for the manufacture of biodegradable and renewable bioproducts. However, the complex lignocellulosic polymeric structure of woody tissue is difficult to access without extensive industrial pre-treatment. Enzyme processing of partly depolymerised biomass is an established technology, and there is evidence that high temperature (extremely thermophilic) lignocellulose degrading enzymes [carbohydrate active enzymes (CAZymes)] may enhance processing efficiency. However, wild-type thermophilic CAZymes will not necessarily be functionally optimal under industrial pre-treatment conditions. With recent advances in synthetic biology, it is now potentially possible to build CAZyme constructs from individual protein domains, tailored to the conditions of specific industrial processes. In this review, we identify a 'toolbox' of thermostable CAZyme domains from extremely thermophilic organisms and highlight recent advances in CAZyme engineering which will allow for the rational design of CAZymes tailored to specific aspects of lignocellulose digestion.

Key Message: QTL consistent across seasons were detected for resistance to cassava brown streak disease induced root necrosis and foliar symptoms. The CMD2 locus was detected in an East African landrace, and comprised two QTL. Cassava production in Africa is compromised by cassava brown streak disease (CBSD) and cassava mosaic disease (CMD). To reduce costs and increase the precision of resistance breeding, a QTL study was conducted to identify molecular markers linked to resistance against these diseases. A bi-parental F 1 mapping population was developed from a cross between the Tanzanian farmer varieties, Namikonga and Albert. A one-step genetic linkage map comprising 943 SNP markers and 18 linkage groups spanning 1776.2 cM was generated. Phenotypic data from 240 F 1 progeny were obtained from two disease hotspots in Tanzania, over two successive seasons, 2013 and 2014. Two consistent QTLs linked to resistance to CBSD-induced root necrosis were identified in Namikonga on chromosomes II (qCBSDRNFc2Nm) and XI (qCBSDRNc11Nm) and a putative QTL on chromosome XVIII (qCBSDRNc18Nm). qCBSDRNFc2Nm was identified at Naliendele in both seasons. The same QTL was also associated with CBSD foliar resistance. qCBSDRNc11Nm was identified at Chambezi in both seasons, and was characterized by three peaks, spanning a distance of 253 kb. Twenty-seven genes were identified within this region including two LRR proteins and a signal recognition particle. In addition, two highly significant CMD resistance QTL (qCMDc12.1A and qCMDc12.2A) were detected in Albert, on chromosome 12. Both qCMDc12.1A and qCMDc12.2A lay within the range of markers reported earlier, defining the CMD2 locus. This is the first time that two loci have been identified within the CMD2 QTL, and in germplasm of apparent East African origin. Additional QTLs with minor effects on CBSD and CMD resistance were also identified.

Conversion of lignocellulosic biomass to biofuels and biomaterials suffers from high production costs associated with biomass pretreatment and enzymatic hydrolysis. In-planta expression of lignocellulose-digesting enzymes is a promising approach to reduce these cost elements. However, this approach faces a number of challenges, including auto-hydrolysis of developing cell walls, plant growth and yield penalties, low expression levels and the limited stability of expressed enzymes at the high temperatures generally used for biomass processing to release fermentable sugars. To overcome these challenges we expressed codon-optimized recombinant hyperthermophilic endoglucanase (EG) and xylanase (Xyn) genes in A. thaliana. Transgenic Arabidopsis lines expressing EG and Xyn enzymes at high levels without any obvious plant growth or yield penalties were selected for further analysis. The highest enzyme activities were observed in the dry stems of transgenic lines, indicating that the enzymes were not degraded during stem senescence and storage. Biomass from transgenic lines exhibited improved saccharification efficiency relative to WT control plants. We conclude that the expression of hyperthermophilic enzymes in plants is a promising approach for combining pretreatment and enzymatic hydrolysis processes in lignocellulosic digestion. This study provides a valid foundation for further studies involving in planta co-expression of core and accessory lignocellulose-digesting enzymes.

A time-course transcriptome analysis of two cassava varieties that are either resistant or susceptible to cassava brown streak disease (CBSD) was conducted using RNASeq, after graft inoculation with Ugandan cassava brown streak virus (UCBSV). From approximately 1.92 billion short reads, the largest number of differentially expressed genes (DEGs) was obtained in the resistant (Namikonga) variety at 2 days after grafting (dag) (3887 DEGs) and 5 dag (4911 DEGs). At the same time points, several defense response genes (encoding LRR-containing, NBARC-containing, pathogenesis-related, late embryogenesis abundant, selected transcription factors, chaperones, and heat shock proteins) were highly expressed in Namikonga. Also, defense-related GO terms of 'translational elongation', 'translation factor activity', 'ribosomal subunit' and 'phosphorelay signal transduction', were overrepresented in Namikonga at these time points. More reads corresponding to UCBSV sequences were recovered from the susceptible variety (Albert) (733 and 1660 read counts per million (cpm)) at 45 dag and 54 dag compared to Namikonga (10 and 117 cpm respectively). These findings suggest that Namikonga's resistance involves restriction of multiplication of UCBSV within the host. These findings can be used with other sources of evidence to identify candidate genes and biomarkers that would contribute substantially to knowledge-based resistance breeding.

Genetic mapping of quantitative trait loci (QTL) for resistance to cassava brown streak disease (CBSD), cassava mosaic disease (CMD), and cassava green mite (CGM) was performed using an F 1 cross developed between the Tanzanian landrace, Kiroba, and a breeding line, AR37-80. The population was evaluated for two consecutive years in two sites in Tanzania. A genetic linkage map was derived from 106 F 1 progeny and 1,974 SNP markers and spanned 18 chromosomes covering a distance of 1,698 cM. Fifteen significant QTL were identified; two are associated with CBSD root necrosis only, and were detected on chromosomes V and XII, while seven were associated with CBSD foliar symptoms only and were detected on chromosomes IV, VI, XVII, and XVIII. QTL on chromosomes 11 and 15 were associated with both CBSD foliar and root necrosis symptoms. Two QTL were found to be associated with CMD and were detected on chromosomes XII and XIV, while two were associated with CGM and were identified on chromosomes V and X. There are large Manihot glaziovii introgression regions in Kiroba on chromosomes I, XVII, and XVIII. The introgression segments on chromosomes XVII and XVIII overlap with QTL associated with CBSD foliar symptoms. The introgression region on chromosome I is of a different haplotype to the characteristic "Amani haplotype" found in the landrace Namikonga and others, and unlike some other genotypes, Kiroba does not have a large introgression block on chromosome IV. Kiroba is closely related to a sampled Tanzanian "tree cassava." This supports the observation that some of the QTL associated with CBSD resistance in Kiroba are different to those observed in another variety, Namikonga.

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.

Eucalyptus species are cultivated for forestry and are of economic importance. The fungal stem canker pathogen Chrysoporthe austroafricana causes disease of varying severity on E. grandis. The Eucalyptus grandis-Chrysoporthe austroafricana interaction has been established as a model system for studying Eucalyptus antifungal defence. Previous studies revealed that the phytohormone salicylic acid (SA) affects the levels of resistance in highly susceptible (ZG14) and moderately resistant (TAG5) clones. The aims of this study were to examine histochemical changes in response to wounding and inoculation as well as host responses at the protein level. The anatomy and histochemical changes induced by wounding and inoculation were similar between the clones, suggesting that anatomical differences do not underlie their different levels of resistance. Tyloses and gum-like substances were present after inoculation and wounding, but cell death occurred only after inoculation. Hyphae of C. austroafricana were observed inside dead and living cells, suggesting that the possibility of a hemibiotrophic interaction requires further investigation. Proteomics analysis revealed the possible involvement of proteins associated with cell death, SA signalling and systemic resistance. In combination with previous information, this study forms a basis for future functional characterisation of candidate genes involved in resistance of E. grandis to C. austroafricana.

We used a systems genetics approach to elucidate the molecular mechanisms of the responses of maize to grey leaf spot (GLS) disease caused by Cercospora zeina, a threat to maize production globally. Expression analysis of earleaf samples in a subtropical maize recombinant inbred line population (CML444 × SC Malawi) subjected in the field to C. zeina infection allowed detection of 20 206 expression quantitative trait loci (eQTLs). Four trans-eQTL hotspots coincided with GLS disease QTLs mapped in the same field experiment. Co-expression network analysis identified three expression modules correlated with GLS disease scores. The module (GY-s) most highly correlated with susceptibility (r = 0.71; 179 genes) was enriched for the glyoxylate pathway, lipid metabolism, diterpenoid biosynthesis and responses to pathogen molecules such as chitin. The GY-s module was enriched for genes with trans-eQTLs in hotspots on chromosomes 9 and 10, which also coincided with phenotypic QTLs for susceptibility to GLS. This transcriptional network has significant overlap with the GLS susceptibility response of maize line B73, and may reflect pathogen manipulation for nutrient acquisition and/or unsuccessful defence responses, such as kauralexin production by the diterpenoid biosynthesis pathway. The co-expression module that correlated best with resistance (TQ-r; 1498 genes) was enriched for genes with trans-eQTLs in hotspots coinciding with GLS resistance QTLs on chromosome 9. Jasmonate responses were implicated in resistance to GLS through co-expression of COI1 and enrichment of genes with the Gene Ontology term 'cullin-RING ubiquitin ligase complex' in the TQ-r module. Consistent with this, JAZ repressor expression was highly correlated with the severity of GLS disease in the GY-s susceptibility network., (© 2016 The Authors The Plant Journal © 2016 John Wiley & Sons Ltd.)

As a consequence of their remarkable adaptability, fast growth, and superior wood properties, eucalypt tree plantations have emerged as key renewable feedstocks (over 20 million ha globally) for the production of pulp, paper, bioenergy, and other lignocellulosic products. However, most biomass properties such as growth, wood density, and wood chemistry are complex traits that are hard to improve in long-lived perennials. Systems genetics, a process of harnessing multiple levels of component trait information (e.g., transcript, protein, and metabolite variation) in populations that vary in complex traits, has proven effective for dissecting the genetics and biology of such traits. We have applied a network-based data integration (NBDI) method for a systems-level analysis of genes, processes and pathways underlying biomass and bioenergy-related traits using a segregating Eucalyptus hybrid population. We show that the integrative approach can link biologically meaningful sets of genes to complex traits and at the same time reveal the molecular basis of trait variation. Gene sets identified for related woody biomass traits were found to share regulatory loci, cluster in network neighborhoods, and exhibit enrichment for molecular functions such as xylan metabolism and cell wall development. These findings offer a framework for identifying the molecular underpinnings of complex biomass and bioprocessing-related traits. A more thorough understanding of the molecular basis of plant biomass traits should provide additional opportunities for the establishment of a sustainable bio-based economy., Competing Interests: The authors declare no conflict of interest.

Chrysoporthe austroafricana is a fungal pathogen that causes the development of stem cankers on susceptible Eucalyptus grandis trees. Clones of E. grandis that are partially resistant and highly susceptible have been identified based on the extent of lesion formation on the stem upon inoculation with C. austroafricana. These interactions have been used as a model pathosystem to enhance our understanding of interactions between pathogenic fungi and woody hosts, which may be different to herbaceous hosts. In previous research, transcriptomics of host responses in these two clones to C. austroafricana suggested roles for salicylic acid and gibberellic acid phytohormone signaling in defense. However, it is unclear how the pathogen infiltrates host tissue and which pathogenicity factors facilitate its spread in the two host genotypes. The aim of this study was to investigate these two aspects of the E. grandis-C. austroafricana interaction and to test the hypothesis that the pathogen possesses mechanisms to modulate the tree phytohormone-mediated defenses. Light microscopy showed that the pathogen occurred in most cell types and structures within infected E. grandis stem tissue. Notably, the fungus appeared to spread through the stem by penetrating cell wall pits. In order to understand the molecular interaction between these organisms and predict putative pathogenicity mechanisms of C. austroafricana , fungal gene expression was studied in vitro and in planta . Fungal genes associated with cell wall degradation, carbohydrate metabolism and phytohormone manipulation were expressed in planta by C. austroafricana . These genes could be involved in fungal spread by facilitating cell wall pit degradation and manipulating phytohormone mediated defense in each host environment, respectively. Specifically, the in planta expression of an ent-kaurene oxidase and salicylate hydroxylase in C. austroafricana suggests putative mechanisms by which the pathogen can modulate the phytohormone-mediated defenses of the host. These mechanisms have been reported in herbaceous plant-pathogen interactions, supporting the notion that these aspects of the interaction are similar in a woody species. This study highlights ent-kaurene oxidase and salicylate hydroxylase as candidates for further functional characterization.

Gall-inducing insects are capable of exerting a high level of control over their hosts' cellular machinery to the extent that the plant's development, metabolism, chemistry, and physiology are all altered in favour of the insect. Many gallers are devastating pests in global agriculture and the limited understanding of their relationship with their hosts prevents the development of robust management strategies. Omics technologies are proving to be important tools in elucidating the mechanisms involved in the interaction as they facilitate analysis of plant hosts and insect effectors for which little or no prior knowledge exists. In this review, we examine the mechanisms behind insect gall development using evidence from omics-level approaches. The secretion of effector proteins and induced phytohormonal imbalances are highlighted as likely mechanisms involved in gall development. However, understanding how these components function within the system is far from complete and a number of questions need to be answered before this information can be used in the development of strategies to engineer or breed plants with enhanced resistance., Competing Interests: The authors declare no conflict of interest.

Many plant genes are known to be involved in the development of cambium and wood, but how the expression and functional interaction of these genes determine the unique biology of wood remains largely unknown. We used the soc1ful loss of function mutant - the woodiest genotype known in the otherwise herbaceous model plant Arabidopsis - to investigate the expression and interactions of genes involved in secondary growth (wood formation). Detailed anatomical observations of the stem in combination with mRNA sequencing were used to assess transcriptome remodeling during xylogenesis in wild-type and woody soc1ful plants. To interpret the transcriptome changes, we constructed functional gene association networks of differentially expressed genes using the STRING database. This analysis revealed functionally enriched gene association hubs that are differentially expressed in herbaceous and woody tissues. In particular, we observed the differential expression of genes related to mechanical stress and jasmonate biosynthesis/signaling during wood formation in soc1ful plants that may be an effect of greater tension within woody tissues. Our results suggest that habit shifts from herbaceous to woody life forms observed in many angiosperm lineages could have evolved convergently by genetic changes that modulate the gene expression and interaction network, and thereby redeploy the conserved wood developmental program., (© 2016 The Authors. The Plant Journal published by Society for Experimental Biology and John Wiley & Sons Ltd.)

In woody plants, xylogenesis is an exceptionally strong carbon sink requiring robust transcriptional control and dynamic coordination of cellular and metabolic processes directing carbon allocation and partitioning into secondary cell wall biosynthesis. As a biological process, wood formation is an excellent candidate for systems modeling due to the strong correlation patterns and interconnectedness observed for transcriptional and metabolic component traits contributing to complex phenotypes such as cell wall chemistry and ultrastructure. Genetic variation in undomesticated tree populations provides abundant perturbation of systems components, adding another dimension to plant systems biology (besides spatial and temporal variation). High-throughput analysis of molecular component traits in adult trees has provided the first insights into the systems genetics of wood, an important renewable feedstock for biomaterials and bioenergy., (Copyright © 2016 Elsevier Ltd. All rights reserved.)

Damage caused by Phytophthora cinnamomi Rands remains an important concern on forest tree species. The pathogen causes root and collar rot, stem cankers, and dieback of various economically important Eucalyptus spp. In South Africa, susceptible cold tolerant Eucalyptus plantations have been affected by various Phytophthora spp. with P. cinnamomi considered one of the most virulent. The molecular basis of this compatible interaction is poorly understood. In this study, susceptible Eucalyptus nitens plants were stem inoculated with P. cinnamomi and tissue was harvested five days post inoculation. Dual RNA-sequencing, a technique which allows the concurrent detection of both pathogen and host transcripts during infection, was performed. Approximately 1% of the reads mapped to the draft genome of P. cinnamomi while 78% of the reads mapped to the Eucalyptus grandis genome. The highest expressed P. cinnamomi gene in planta was a putative crinkler effector (CRN1). Phylogenetic analysis indicated the high similarity of this P. cinnamomi CRN1 to that of Phytophthora infestans. Some CRN effectors are known to target host nuclei to suppress defense. In the host, over 1400 genes were significantly differentially expressed in comparison to mock inoculated trees, including suites of pathogenesis related (PR) genes. In particular, a PR-9 peroxidase gene with a high similarity to a Carica papaya PR-9 ortholog previously shown to be suppressed upon infection by Phytophthora palmivora was down-regulated two-fold. This PR-9 gene may represent a cross-species effector target during P. cinnamomi infection. This study identified pathogenicity factors, potential manipulation targets, and attempted host defense mechanisms activated by E. nitens that contributed to the susceptible outcome of the interaction.

DUF1218 is a land plant-specific innovation and has previously been shown to be associated with cell wall biology, vasculature patterning and abiotic/biotic stress response. The Arabidopsis genome encodes 15 members, two of which (At1g31720 and At4g27435) are preferentially expressed in the secondary cell wall depositing inflorescence stems. To further our understanding of the roles of DUF1218-containing proteins in secondary cell wall biology, we functionally characterized At1g31720 (herein referred to as MODIFYING WALL LIGNIN-1 or MWL-1). Since related gene family members may contribute to functional redundancy, we also characterized At4g19370 (MWL-2), the most closely related gene to MWL-1 in the protein family. Subcellular localization revealed that both Arabidopsis proteins are targeted to the cell periphery. The single T-DNA knockout lines, mwl-1 and mwl-2, and independent overexpression lines showed no significant differences in plant growth or changes in total lignin content relative to wild-type (WT) control plants. However, the double homozygous mutant, mwl-1/mwl-2, had smaller rosettes with a significant decrease in rosette fresh weight and stem height relative to the WT control at four weeks and six weeks, respectively. Moreover, mwl-1/mwl-2 showed a significant reduction in total lignin content (by ca. 11% relative to WT) and an increase in syringyl/guaiacyl (S/G) monomer ratio relative to the control plants. Our study has identified two additional members of the DUF1218 family in Arabidopsis as novel contributors to secondary cell wall biology, specifically lignin biosynthesis, and these proteins appear to function redundantly.

Background: Pines are the most important tree species to the international forestry industry, covering 42 % of the global industrial forest plantation area. One of the most pressing threats to cultivation of some pine species is the pitch canker fungus, Fusarium circinatum, which can have devastating effects in both the field and nursery. Investigation of the Pinus-F. circinatum host-pathogen interaction is crucial for development of effective disease management strategies. As with many non-model organisms, investigation of host-pathogen interactions in pine species is hampered by limited genomic resources. This was partially alleviated through release of the 22 Gbp Pinus taeda v1.01 genome sequence ( http://pinegenome.org/pinerefseq/ ) in 2014. Despite the fact that the fragmented state of the genome may hamper comprehensive transcriptome analysis, it is possible to leverage the inherent redundancy resulting from deep RNA sequencing with Illumina short reads to assemble transcripts in the absence of a completed reference sequence. These data can then be integrated with available genomic data to produce a comprehensive transcriptome resource. The aim of this study was to provide a foundation for gene expression analysis of disease response mechanisms in Pinus patula through transcriptome assembly., Results: Eighteen de novo and two reference based assemblies were produced for P. patula shoot tissue. For this purpose three transcriptome assemblers, Trinity, Velvet/OASES and SOAPdenovo-Trans, were used to maximise diversity and completeness of assembled transcripts. Redundancy in the assembly was reduced using the EvidentialGene pipeline. The resulting 52 Mb P. patula v1.0 shoot transcriptome consists of 52 112 unigenes, 60 % of which could be functionally annotated., Conclusions: The assembled transcriptome will serve as a major genomic resource for future investigation of P. patula and represents the largest gene catalogue produced to date for this species. Furthermore, this assembly can help detect gene-based genetic markers for P. patula and the comparative assembly workflow could be applied to generate similar resources for other non-model species.

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 Author 2015. Published by Oxford University Press on behalf of Japanese Society of Plant Physiologists. All rights reserved. For permissions, please email: journals.permissions@oup.com.)

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.

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., (© 2014 The Authors. New Phytologist © 2014 New Phytologist Trust.)

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., (© 2014 The Authors. New Phytologist © 2014 New Phytologist Trust.)

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., (© 2015 University of Tasmania New Phytologist © 2015 New Phytologist Trust.)

Horsetails represent an enigmatic clade within the land plants. Despite consisting only of one genus (Equisetum) that contains 15 species, they are thought to represent the oldest extant genus within the vascular plants dating back possibly as far as the Triassic. Horsetails have retained several ancient features and are also characterized by a particularly high chromosome count (n = 108). Whole-genome duplications (WGDs) have been uncovered in many angiosperm clades and have been associated with the success of angiosperms, both in terms of species richness and biomass dominance, but remain understudied in nonangiosperm clades. Here, we report unambiguous evidence of an ancient WGD in the fern lineage, based on sequencing and de novo assembly of an expressed gene catalog (transcriptome) from the giant horsetail (Equisetum giganteum). We demonstrate that horsetails underwent an independent paleopolyploidy during the Late Cretaceous prior to the diversification of the genus but did not experience any recent polyploidizations that could account for their high chromosome number. We also discuss the specific retention of genes following the WGD and how this may be linked to their long-term survival., (© 2015 American Society of Plant Biologists. All rights reserved.)

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., (© 2014 University of Pretoria New Phytologist © 2014 New Phytologist Trust.)