Search

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.)

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

Tension wood has distinct physical and chemical properties, including altered fibre properties, cell wall composition and ultrastructure. It serves as a good system for investigating the genetic regulation of secondary cell wall biosynthesis and wood formation. The reference genome sequence for Eucalyptus grandis allows investigation of the global transcriptional reprogramming that accompanies tension wood formation in this global wood fibre crop. We report the first comprehensive analysis of physicochemical wood property changes in tension wood of Eucalyptus measured in a hybrid (E. grandis × Eucalyptus urophylla) clone, as well as genome-wide gene expression changes in xylem tissues 3 wk post-induction using RNA sequencing. We found that Eucalyptus tension wood in field-grown trees is characterized by an increase in cellulose, a reduction in lignin, xylose and mannose, and a marked increase in galactose. Gene expression profiling in tension wood-forming tissue showed corresponding down-regulation of monolignol biosynthetic genes, and differential expression of several carbohydrate active enzymes. We conclude that alterations of cell wall traits induced by tension wood formation in Eucalyptus are a consequence of a combination of down-regulation of lignin biosynthesis and hemicellulose remodelling, rather than the often proposed up-regulation of the cellulose biosynthetic pathway., (© 2014 University of Pretoria New Phytologist © 2014 New Phytologist Trust.)

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

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

Background: Carbohydrate metabolism is a key feature of vascular plant architecture, and is of particular importance in large woody species, where lignocellulosic biomass is responsible for bearing the bulk of the stem and crown. Since Carbohydrate Active enZymes (CAZymes) in plants are responsible for the synthesis, modification and degradation of carbohydrate biopolymers, the differences in gene copy number and regulation between woody and herbaceous species have been highlighted previously. There are still many unanswered questions about the role of CAZymes in land plant evolution and the formation of wood, a strong carbohydrate sink., Results: Here, twenty-two publically available plant genomes were used to characterize the frequency, diversity and complexity of CAZymes in plants. We find that a conserved suite of CAZymes is a feature of land plant evolution, with similar diversity and complexity regardless of growth habit and form. In addition, we compared the diversity and levels of CAZyme gene expression during wood formation in trees using mRNA-seq data from two distantly related angiosperm tree species Eucalyptus grandis and Populus trichocarpa, highlighting the major CAZyme classes involved in xylogenesis and lignocellulosic biomass production., Conclusions: CAZyme domain ratio across embryophytes is maintained, and the diversity of CAZyme domains is similar in all land plants, regardless of woody habit. The stoichiometric conservation of gene expression in woody and non-woody tissues of Eucalyptus and Populus are indicative of gene balance preservation.

Background: Histone 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., Results: Plant 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., Conclusions: In 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.

Background: Eucalyptus species and interspecific hybrids exhibit valuable growth and wood properties that make them a highly desirable commodity. However, these trees are challenged by a wide array of biotic stresses during their lifetimes. The Eucalyptus grandis reference genome sequence provides a resource to study pest and pathogen defence mechanisms in long-lived woody plants. E. grandis trees are generally susceptible to Chrysoporthe austroafricana, a causal agent of stem cankers on eucalypts. The aim of this study was to characterize the defence response of E. grandis against C. austroafricana., Results: Hormone profiling of susceptible and moderately resistant clonal E. grandis genotypes indicated a reduction in salicylic acid and gibberellic acid levels at 3 days post inoculation. We hypothesized that these signaling pathways may facilitate resistance. To further investigate other defence mechanisms at this time point, transcriptome profiling was performed. This revealed that cell wall modifications and response to oxidative stress form part of the defence responses common to both genotypes, whilst changes in the hormone signaling pathways may contribute to resistance. Additionally the expression of selected candidate defence response genes was induced earlier in moderately resistant trees than in susceptible trees, supporting the hypothesis that a delayed defence response may occur in the susceptible interaction., Conclusion: The ability of a host to fine-tune its defence responses is crucial and the responses identified in this study extends our understanding of plant defence, gained from model systems, to woody perennials.

Auxin plays a pivotal role in various plant growth and development processes, including vascular differentiation. The modulation of auxin responsiveness through the auxin perception and signaling machinery is believed to be a major regulatory mechanism controlling cambium activity and wood formation. To gain more insights into the roles of key Aux/IAA gene regulators of the auxin response in these processes, we identified and characterized members of the Aux/IAA family in the genome of Eucalyptus grandis, a tree of worldwide economic importance. We found that the gene family in Eucalyptus is slightly smaller than that in Populus and Arabidopsis, but all phylogenetic groups are represented. High-throughput expression profiling of different organs and tissues highlighted several Aux/IAA genes expressed in vascular cambium and/or developing xylem, some showing differential expression in response to developmental (juvenile vs. mature) and/or to environmental (tension stress) cues. Based on the expression profiles, we selected a promising candidate gene, EgrIAA4, for functional characterization. We showed that EgrIAA4 protein is localized in the nucleus and functions as an auxin-responsive repressor. Overexpressing a stabilized version of EgrIAA4 in Arabidopsis dramatically impeded plant growth and fertility and induced auxin-insensitive phenotypes such as inhibition of primary root elongation, lateral root emergence and agravitropism. Interestingly, the lignified secondary walls of the interfascicular fibers appeared very late, whereas those of the xylary fibers were virtually undetectable, suggesting that EgrIAA4 may play crucial roles in fiber development and secondary cell wall deposition., (© 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.)

Plant organisms contain a large number of genes belonging to numerous multigenic families whose evolution size reflects some functional constraints. Sequences from eight multigenic families, involved in biotic and abiotic responses, have been analyzed in Eucalyptus grandis and compared with Arabidopsis thaliana. Two transcription factor families APETALA 2 (AP2)/ethylene responsive factor and GRAS, two auxin transporter families PIN-FORMED and AUX/LAX, two oxidoreductase families (ascorbate peroxidases [APx] and Class III peroxidases [CIII Prx]), and two families of protective molecules late embryogenesis abundant (LEA) and DNAj were annotated in expert and exhaustive manner. Many recent tandem duplications leading to the emergence of species-specific gene clusters and the explosion of the gene numbers have been observed for the AP2, GRAS, LEA, PIN, and CIII Prx in E. grandis, while the APx, the AUX/LAX and DNAj are conserved between species. Although no direct evidence has yet demonstrated the roles of these recent duplicated genes observed in E. grandis, this could indicate their putative implications in the morphological and physiological characteristics of E. grandis, and be the key factor for the survival of this nondormant species. Global analysis of key families would be a good criterion to evaluate the capabilities of some organisms to adapt to environmental variations., (© The Author(s) 2015. Published by Oxford University Press on behalf of the Society for Molecular Biology and Evolution.)

Auxin is a central hormone involved in a wide range of developmental processes including the specification of vascular stem cells. Auxin Response Factors (ARF) are important actors of the auxin signalling pathway, regulating the transcription of auxin-responsive genes through direct binding to their promoters. The recent availability of the Eucalyptus grandis genome sequence allowed us to examine the characteristics and evolutionary history of this gene family in a woody plant of high economic importance. With 17 members, the E. grandis ARF gene family is slightly contracted, as compared to those of most angiosperms studied hitherto, lacking traces of duplication events. In silico analysis of alternative transcripts and gene truncation suggested that these two mechanisms were preeminent in shaping the functional diversity of the ARF family in Eucalyptus. Comparative phylogenetic analyses with genomes of other taxonomic lineages revealed the presence of a new ARF clade found preferentially in woody and/or perennial plants. High-throughput expression profiling among different organs and tissues and in response to environmental cues highlighted genes expressed in vascular cambium and/or developing xylem, responding dynamically to various environmental stimuli. Finally, this study allowed identification of three ARF candidates potentially involved in the auxin-regulated transcriptional program underlying wood formation.

Long-lived tree species are subject to attack by various pests and pathogens during their lifetime. This problem is exacerbated by climate change, which may increase the host range for pathogens and extend the period of infestation by pests. Plant defences may involve preformed barriers or induced resistance mechanisms based on recognition of the invader, complex signalling cascades, hormone signalling, activation of transcription factors and production of pathogenesis-related (PR) proteins with direct antimicrobial or anti-insect activity. Trees have evolved some unique defence mechanisms compared with well-studied model plants, which are mostly herbaceous annuals. The genome sequence of Eucalyptus grandis W. Hill ex Maiden has recently become available and provides a resource to extend our understanding of defence in large woody perennials. This review synthesizes existing knowledge of defence mechanisms in model plants and tree species and features mechanisms that may be important for defence in Eucalyptus, such as anatomical variants and the role of chemicals and proteins. Based on the E. grandis genome sequence, we have identified putative PR proteins based on sequence identity to the previously described plant PR proteins. Putative orthologues for PR-1, PR-2, PR-4, PR-5, PR-6, PR-7, PR-8, PR-9, PR-10, PR-12, PR-14, PR-15 and PR-17 have been identified and compared with their orthologues in Populus trichocarpa Torr. & A. Gray ex Hook and Arabidopsis thaliana (L.) Heynh. The survey of PR genes in Eucalyptus provides a first step in identifying defence gene targets that may be employed for protection of the species in future. Genomic resources available for Eucalyptus are discussed and approaches for improving resistance in these hardwood trees, earmarked as a bioenergy source in future, are considered., (© The Author 2014. Published by Oxford University Press. All rights reserved. For Permissions, please email: journals.permissions@oup.com.)

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.

Lignocellulosic biomass is an important feedstock for the pulp and paper industry as well as emerging biofuel and biomaterial industries. However, the recalcitrance of the secondary cell wall to chemical or enzymatic degradation remains a major hurdle for efficient extraction of economically important biopolymers such as cellulose. It has been estimated that approximately 10-15% of about 27,000 protein-coding genes in the Arabidopsis genome are dedicated to cell wall development; however, only about 130 Arabidopsis genes thus far have experimental evidence validating cell wall function. While many genes have been implicated through co-expression analysis with known genes, a large number are broadly classified as proteins of unknown function (PUFs). Recently the functionality of some of these unknown proteins in cell wall development has been revealed using reverse genetic approaches. Given the large number of cell wall-related PUFs, how do we approach and subsequently prioritize the investigation of such unknown genes that may be essential to or influence plant cell wall development and structure? Here, we address the aforementioned question in two parts; we first identify the different kinds of PUFs based on known and predicted features such as protein domains. Knowledge of inherent features of PUFs may allow for functional inference and a concomitant link to biological context. Secondly, we discuss omics-based technologies and approaches that are helping identify and prioritize cell wall-related PUFs by functional association. In this way, hypothesis-driven experiments can be designed for functional elucidation of many proteins that remain missing links in our understanding of plant cell wall biosynthesis., (© The Author 2014. 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: Gray leaf spot (GLS) is a globally important foliar disease of maize. Cercospora zeina, one of the two fungal species that cause the disease, is prevalent in southern Africa, China, Brazil and the eastern corn belt of the USA. Identification of QTL for GLS resistance in subtropical germplasm is important to support breeding programmes in developing countries where C. zeina limits production of this staple food crop., Results: A maize RIL population (F7:S6) from a cross between CML444 and SC Malawi was field-tested under GLS disease pressure at five field sites over three seasons in KwaZulu-Natal, South Africa. Thirty QTL identified from eleven field trials (environments) were consolidated to seven QTL for GLS resistance based on their expression in at least two environments and location in the same core maize bins. Four GLS resistance alleles were derived from the more resistant parent CML444 (bin 1.10, 4.08, 9.04/9.05, 10.06/10.07), whereas the remainder were from SC Malawi (bin 6.06/6.07, 7.02/7.03, 9.06). QTLs in bin 4.08 and bin 6.06/6.07 were also detected as joint QTLs, each explained more than 11% of the phenotypic variation, and were identified in four and seven environments, respectively. Common markers were used to allocate GLS QTL from eleven previous studies to bins on the IBM2005 map, and GLS QTL "hotspots" were noted. Bin 4.08 and 7.02/7.03 GLS QTL from this study overlapped with hotspots, whereas the bin 6.06/6.07 and bin 9.06 QTLs appeared to be unique. QTL for flowering time (bin 1.07, 4.09) in this population did not correspond to QTL for GLS resistance., Conclusions: QTL mapping of a RIL population from the subtropical maize parents CML444 and SC Malawi identified seven QTL for resistance to gray leaf spot disease caused by C. zeina. These QTL together with QTL from eleven studies were allocated to bins on the IBM2005 map to provide a basis for comparison. Hotspots of GLS QTL were identified on chromosomes one, two, four, five and seven, with QTL in the current study overlapping with two of these. Two QTL from this study did not overlap with previously reported QTL.

Plant biomass, as an abundant renewable carbon source, is a promising alternative to fossil fuels. However, the enzymes most commonly used for depolymerization of lignocellulosic biomass are expensive, and the development of cost-effective alternative conversion technologies would be desirable. One possible option is the heterologous expression of genes encoding lignocellulose-digesting enzymes in plant tissues. To overcome simultaneously issues of toxicity and incompatibility with high-temperature steam explosion processes, the use of heterologous genes encoding hyperthermophilic enzymes may be an attractive alternative. This approach could reduce the need for exogenous enzyme additions prior to fermentation, reducing the cost of the complete processing operation. This review highlights recent advances and future prospects for using hyperthermophilic enzymes in the biofuels industry., (Copyright © 2014 Elsevier Ltd. All rights reserved.)

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.

Two important role players in plant defence response are the phytohormones salicylic acid (SA) and jasmonic acid (JA); both of which have been well described in model species such as Arabidopsis thaliana. Several pathogenesis related (PR) genes have previously been used as indicators of the onset of SA and JA signaling in Arabidopsis. This information is lacking in tree genera such as Eucalyptus. The aim of this study was to characterize the transcriptional response of PR genes (EgrPR2, EgrPR3, EgrPR4, EgrPR5, and EgrLOX) identified in Eucalyptus grandis to SA and methyl jasmonate (MeJA) treatment as well as to qualify them as diagnostic for the two signaling pathways. Using the genome sequence of E. grandis, we identified candidate Eucalyptus orthologs EgrPR2, EgrPR3, EgrPR4, EgrPR5, and EgrLOX based on a co-phylogenetic approach. The expression of these genes was investigated after various doses of SA and MeJA (a derivative of JA) treatment as well as at various time points. The transcript levels of EgrPR2 were decreased in response to high concentrations of MeJA whereas the expression of EgrPR3 and EgrLOX declined as the concentrations of SA treatment increased, suggesting an antagonistic relationship between SA and MeJA. Our results support EgrPR2 as potentially diagnostic for SA and EgrPR3, EgrPR4, and EgrLOX as indicators of MeJA signaling. To further validate the diagnostic potential of the PR genes we challenged E. grandis clones with the fungal necrotrophic pathogen Chrysoporthe austroafricana. The tolerant clone showed high induction of EgrPR2 and decreased transcript abundance of EgrPR4. Pre-treatment of the susceptible genotype with 5 mM SA resulted in lesion lengths comparable to the tolerant genotype after artificial inoculation with C. austroafricana. Thus expression profiling of EgrPR2 and EgrPR4 genes could serve as a useful diagnostic approach to determine which of the two signaling pathways are activated against various pathogens in Eucalyptus.

The increasing focus on plantation forestry as a renewable source of cellulosic biomass has emphasized the need for tools to study the unique biology of woody genera such as Eucalyptus, Populus and Pinus. The domestication of these woody crops is hampered by long generation times, and breeders are now looking to molecular approaches such as marker-assisted breeding and genetic modification to accelerate tree improvement. Much of what is known about genes involved in the growth and development of plants has come from studies of herbaceous models such as Arabidopsis and rice. However, transferring this information to woody plants often proves difficult, especially for genes expressed in woody stems. Here we report the use of induced somatic sector analysis (ISSA) for characterization of promoter expression patterns directly in the stems of Populus and Eucalyptus trees. As a case study, we used previously characterized primary and secondary cell wall-related cellulose synthase (CesA) promoters cloned from Eucalyptus grandis. We show that ISSA can be used to elucidate the phloem and xylem expression patterns of the CesA genes in Eucalyptus and Populus stems and also show that the staining patterns differ in Eucalyptus and Populus stems. These findings show that ISSA is an efficient approach to investigate promoter function in the developmental context of woody plant tissues and raise questions about the suitability of heterologous promoters for genetic manipulation in plant species.

Background: F1 hybrid clones of Eucalyptus grandis and E. urophylla are widely grown for pulp and paper production in tropical and subtropical regions. Volume growth and wood quality are priority objectives in Eucalyptus tree improvement. The molecular basis of quantitative variation and trait expression in eucalypt hybrids, however, remains largely unknown. The recent availability of a draft genome sequence (http://www.phytozome.net) and genome-wide genotyping platforms, combined with high levels of genetic variation and high linkage disequilibrium in hybrid crosses, greatly facilitate the detection of quantitative trait loci (QTLs) as well as underlying candidate genes for growth and wood property traits. In this study, we used Diversity Arrays Technology markers to assess the genetic architecture of volume growth (diameter at breast height, DBH) and wood basic density in four-year-old progeny of an interspecific backcross pedigree of E. grandis and E. urophylla. In addition, we used Illumina RNA-Seq expression profiling in the E. urophylla backcross family to identify cis- and trans-acting polymorphisms (eQTLs) affecting transcript abundance of genes underlying QTLs for wood basic density., Results: A total of five QTLs for DBH and 12 for wood basic density were identified in the two backcross families. Individual QTLs for DBH and wood basic density explained 3.1 to 12.2% of phenotypic variation. Candidate genes underlying QTLs for wood basic density on linkage groups 8 and 9 were found to share trans-acting eQTLs located on linkage groups 4 and 10, which in turn coincided with QTLs for wood basic density suggesting that these QTLs represent segregating components of an underlying transcriptional network., Conclusion: This is the first demonstration of the use of next-generation expression profiling to quantify transcript abundance in a segregating tree population and identify candidate genes potentially affecting wood property variation. The QTLs identified in this study provide a resource for identifying candidate genes and developing molecular markers for marker-assisted breeding of volume growth and wood basic density. Our results suggest that integrated analysis of transcript and trait variation in eucalypt hybrids can be used to dissect the molecular basis of quantitative variation in wood property traits.

Background: Genetic linkage maps are invaluable resources in plant research. They provide a key tool for many genetic applications including: mapping quantitative trait loci (QTL); comparative mapping; identifying unlinked (i.e. independent) DNA markers for fingerprinting, population genetics and phylogenetics; assisting genome sequence assembly; relating physical and recombination distances along the genome and map-based cloning of genes. Eucalypts are the dominant tree species in most Australian ecosystems and of economic importance globally as plantation trees. The genome sequence of E. grandis has recently been released providing unprecedented opportunities for genetic and genomic research in the genus. A robust reference linkage map containing sequence-based molecular markers is needed to capitalise on this resource. Several high density linkage maps have recently been constructed for the main commercial forestry species in the genus (E. grandis, E. urophylla and E. globulus) using sequenced Diversity Arrays Technology (DArT) and microsatellite markers. To provide a single reference linkage map for eucalypts a composite map was produced through the integration of data from seven independent mapping experiments (1950 individuals) using a marker-merging method., Results: The composite map totalled 1107 cM and contained 4101 markers; comprising 3880 DArT, 213 microsatellite and eight candidate genes. Eighty-one DArT markers were mapped to two or more linkage groups, resulting in the 4101 markers being mapped to 4191 map positions. Approximately 13% of DArT markers mapped to identical map positions, thus the composite map contained 3634 unique loci at an average interval of 0.31 cM., Conclusion: The composite map represents the most saturated linkage map yet produced in Eucalyptus. As the majority of DArT markers contained on the map have been sequenced, the map provides a direct link to the E. grandis genome sequence and will serve as an important reference for progressing eucalypt research.

Fast-growing, short-rotation forest trees, such as Populus and Eucalyptus, produce large amounts of cellulose-rich biomass that could be utilized for bioenergy and biopolymer production. Major obstacles need to be overcome before the deployment of these genera as energy crops, including the effective removal of lignin and the subsequent liberation of carbohydrate constituents from wood cell walls. However, significant opportunities exist to both select for and engineer the structure and interaction of cell wall biopolymers, which could afford a means to improve processing and product development. The molecular underpinnings and regulation of cell wall carbohydrate biosynthesis are rapidly being elucidated, and are providing tools to strategically develop and guide the targeted modification required to adapt forest trees for the emerging bioeconomy. Much insight has already been gained from the perturbation of individual genes and pathways, but it is not known to what extent the natural variation in the sequence and expression of these same genes underlies the inherent variation in wood properties of field-grown trees. The integration of data from next-generation genomic technologies applied in natural and experimental populations will enable a systems genetics approach to study cell wall carbohydrate production in trees, and should advance the development of future woody bioenergy and biopolymer crops.

Diversity Arrays Technology (DArT) provides a robust, high throughput, cost-effective method to query thousands of sequence polymorphisms in a single assay. Despite the extensive use of this genotyping platform for numerous plant species, little is known regarding the sequence attributes and genome-wide distribution of DArT markers. We investigated the genomic properties of the 7,680 DArT marker probes of a Eucalyptus array, by sequencing them, constructing a high density linkage map and carrying out detailed physical mapping analyses to the Eucalyptus grandis reference genome. A consensus linkage map with 2,274 DArT markers anchored to 210 microsatellites and a framework map, with improved support for ordering, displayed extensive collinearity with the genome sequence. Only 1.4 Mbp of the 75 Mbp of still unplaced scaffold sequence was captured by 45 linkage mapped but physically unaligned markers to the 11 main Eucalyptus pseudochromosomes, providing compelling evidence for the quality and completeness of the current Eucalyptus genome assembly. A highly significant correspondence was found between the locations of DArT markers and predicted gene models, while most of the 89 DArT probes unaligned to the genome correspond to sequences likely absent in E. grandis, consistent with the pan-genomic feature of this multi-Eucalyptus species DArT array. These comprehensive linkage-to-physical mapping analyses provide novel data regarding the genomic attributes of DArT markers in plant genomes in general and for Eucalyptus in particular. DArT markers preferentially target the gene space and display a largely homogeneous distribution across the genome, thereby providing superb coverage for mapping and genome-wide applications in breeding and diversity studies. Data reported on these ubiquitous properties of DArT markers will be particularly valuable to researchers working on less-studied crop species who already count on DArT genotyping arrays but for which no reference genome is yet available to allow such detailed characterization.

Background: NAC 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., Results: We 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., Conclusions: This 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.

Background: Microarray technology has matured over the past fifteen years into a cost-effective solution with established data analysis protocols for global gene expression profiling. The Agilent-016047 maize 44 K microarray was custom-designed from EST sequences, but only reporter sequences with EST accession numbers are publicly available. The following information is lacking: (a) reporter - gene model match, (b) number of reporters per gene model, (c) potential for cross hybridization, (d) sense/antisense orientation of reporters, (e) position of reporter on B73 genome sequence (for eQTL studies), and (f) functional annotations of genes represented by reporters. To address this, we developed a strategy to annotate the Agilent-016047 maize microarray, and built a publicly accessible annotation database., Description: Genomic annotation of the 42,034 reporters on the Agilent-016047 maize microarray was based on BLASTN results of the 60-mer reporter sequences and their corresponding ESTs against the maize B73 RefGen v2 "Working Gene Set" (WGS) predicted transcripts and the genome sequence. The agreement between the EST, WGS transcript and gDNA BLASTN results were used to assign the reporters into six genomic annotation groups. These annotation groups were: (i) "annotation by sense gene model" (23,668 reporters), (ii) "annotation by antisense gene model" (4,330); (iii) "annotation by gDNA" without a WGS transcript hit (1,549); (iv) "annotation by EST", in which case the EST from which the reporter was designed, but not the reporter itself, has a WGS transcript hit (3,390); (v) "ambiguous annotation" (2,608); and (vi) "inconclusive annotation" (6,489). Functional annotations of reporters were obtained by BLASTX and Blast2GO analysis of corresponding WGS transcripts against GenBank.The annotations are available in the Maize Microarray Annotation Database http://MaizeArrayAnnot.bi.up.ac.za/, as well as through a GBrowse annotation file that can be uploaded to the MaizeGDB genome browser as a custom track.The database was used to re-annotate lists of differentially expressed genes reported in case studies of published work using the Agilent-016047 maize microarray. Up to 85% of reporters in each list could be annotated with confidence by a single gene model, however up to 10% of reporters had ambiguous annotations. Overall, more than 57% of reporters gave a measurable signal in tissues as diverse as anthers and leaves., Conclusions: The Maize Microarray Annotation Database will assist users of the Agilent-016047 maize microarray in (i) refining gene lists for global expression analysis, and (ii) confirming the annotation of candidate genes before functional studies.

Fusarium circinatum and Fusarium subglutinans are two distinct species in the Gibberella fujikuroi species complex. A genetic linkage map produced from an interspecific cross between these species was used to identify quantitative trait loci (QTLs) associated with variation in mycelial growth and morphology of colony margins (CMs) in the 94 F(1) progeny. Mycelial growth was assessed by measuring culture size at 25°C and 30°C, while CM morphology was characterized in the parents and assessed in their F(1) progeny. In order to test the pathogenicity of the progeny, Pinus patula seedlings were inoculated and lesion lengths were measured after 3weeks. Seven putative QTLs were associated with mycelial growth, three for growth at 25°C and four at 30°C. One highly significant QTL (P<0.001) was present at both growth temperatures. For CM morphology, a QTL was identified at the same position (P<0.001) as the QTL responsible for growth at the two temperatures. The putative QTLs accounted for 45 and 41% of the total mycelial growth variation at 25°C and 30°C, respectively, and for 21% of the variation in CM morphology. Only one of the 94 F(1) progeny was pathogenic on P. patula seedlings. This observation could be explained by the genetic constitution of this F(1) isolate, namely that ∼96% of its genome originated from the F. circinatum parent. This F(1) individual also grew significantly faster at 25°C than the F. circinatum parent (P<0.05), as well as more rapidly than the average growth for the remaining 93 F(1) progeny (P<0.05). However, no association was found between mycelial growth and pathogenicity at 25°C. The highly significant QTL associated with growth at two temperatures, suggests that this is a principal genomic region involved in mycelial growth at both temperatures, and that the same region is also responsible for CM morphology., (Copyright © 2011 British Mycological Society. Published by Elsevier Ltd. All rights reserved.)

A set of over 8000 Diversity Arrays Technology (DArT) markers was tested for its utility in high-resolution population and phylogenetic studies across a range of Eucalyptus taxa. Small-scale population studies of Eucalyptus camaldulensis, Eucalyptus cladocalyx, Eucalyptus globulus, Eucalyptus grandis, Eucalyptus nitens, Eucalyptus pilularis and Eucalyptus urophylla demonstrated the potential of genome-wide genotyping with DArT markers to differentiate species, to identify interspecific hybrids and to resolve biogeographic disjunctions within species. The population genetic studies resolved geographically partitioned clusters in E. camaldulensis, E. cladocalyx, E. globulus and E. urophylla that were congruent with previous molecular studies. A phylogenetic study of 94 eucalypt species provided results that were largely congruent with traditional taxonomy and ITS-based phylogenies, but provided more resolution within major clades than had been obtained previously. Ascertainment bias (the bias introduced in a phylogeny from using markers developed in a small sample of the taxa that are being studied) was not detected. DArT offers an unprecedented level of resolution for population genetic, phylogenetic and evolutionary studies across the full range of Eucalyptus species., (Copyright © 2011. Published by Elsevier Inc.)

Background: De novo assembly of transcript sequences produced by short-read DNA sequencing technologies offers a rapid approach to obtain expressed gene catalogs for non-model organisms. A draft genome sequence will be produced in 2010 for a Eucalyptus tree species (E. grandis) representing the most important hardwood fibre crop in the world. Genome annotation of this valuable woody plant and genetic dissection of its superior growth and productivity will be greatly facilitated by the availability of a comprehensive collection of expressed gene sequences from multiple tissues and organs., Results: We present an extensive expressed gene catalog for a commercially grown E. grandis × E. urophylla hybrid clone constructed using only Illumina mRNA-Seq technology and de novo assembly. A total of 18,894 transcript-derived contigs, a large proportion of which represent full-length protein coding genes were assembled and annotated. Analysis of assembly quality, length and diversity show that this dataset represent the most comprehensive expressed gene catalog for any Eucalyptus tree. mRNA-Seq analysis furthermore allowed digital expression profiling of all of the assembled transcripts across diverse xylogenic and non-xylogenic tissues, which is invaluable for ascribing putative gene functions., Conclusions: De novo assembly of Illumina mRNA-Seq reads is an efficient approach for transcriptome sequencing and profiling in Eucalyptus and other non-model organisms. The transcriptome resource (Eucspresso, http://eucspresso.bi.up.ac.za/) generated by this study will be of value for genomic analysis of woody biomass production in Eucalyptus and for comparative genomic analysis of growth and development in woody and herbaceous plants.

Background: A number of molecular marker technologies have allowed important advances in the understanding of the genetics and evolution of Eucalyptus, a genus that includes over 700 species, some of which are used worldwide in plantation forestry. Nevertheless, the average marker density achieved with current technologies remains at the level of a few hundred markers per population. Furthermore, the transferability of markers produced with most existing technology across species and pedigrees is usually very limited. High throughput, combined with wide genome coverage and high transferability are necessary to increase the resolution, speed and utility of molecular marker technology in eucalypts. We report the development of a high-density DArT genome profiling resource and demonstrate its potential for genome-wide diversity analysis and linkage mapping in several species of Eucalyptus., Findings: After testing several genome complexity reduction methods we identified the PstI/TaqI method as the most effective for Eucalyptus and developed 18 genomic libraries from PstI/TaqI representations of 64 different Eucalyptus species. A total of 23,808 cloned DNA fragments were screened and 13,300 (56%) were found to be polymorphic among 284 individuals. After a redundancy analysis, 6,528 markers were selected for the operational array and these were supplemented with 1,152 additional clones taken from a library made from the E. grandis tree whose genome has been sequenced. Performance validation for diversity studies revealed 4,752 polymorphic markers among 174 individuals. Additionally, 5,013 markers showed segregation when screened using six inter-specific mapping pedigrees, with an average of 2,211 polymorphic markers per pedigree and a minimum of 859 polymorphic markers that were shared between any two pedigrees., Conclusions: This operational DArT array will deliver 1,000-2,000 polymorphic markers for linkage mapping in most eucalypt pedigrees and thus provide high genome coverage. This array will also provide a high-throughput platform for population genetics and phylogenetics in Eucalyptus. The transferability of DArT across species and pedigrees is particularly valuable for a large genus such as Eucalyptus and will facilitate the transfer of information between different studies. Furthermore, the DArT marker array will provide a high-resolution link between phenotypes in populations and the Eucalyptus reference genome, which will soon be completed.

* The cellulose synthase (CesA) gene family encodes the catalytic subunits of a large protein complex responsible for the deposition of cellulose into plant cell walls. Early in vascular plant evolution, the gene family diverged into distinct members with conserved structures and functions (e.g. primary or secondary cell wall biosynthesis). Although the functions and expression domains of CesA genes have been extensively studied in plants, little is known about transcriptional regulation and promoter evolution in this gene family. * Here, comparative sequence analysis of orthologous CesA promoters from three angiosperm genera, Arabidopsis, Populus and Eucalyptus, was performed to identify putative cis-regulatory sequences. The promoter sequences of groups of Arabidopsis genes that are co-expressed with the primary or secondary cell wall-related CesA genes were also analyzed. * Reporter gene analysis of newly isolated promoter regions of six E. grandis CesA genes in Arabidopsis revealed the conserved functionality of the promoter sequences. Comparative sequence analysis identified 71 conserved sequence motifs, of which 66 were significantly over-represented in either primary or secondary wall-associated promoters. * The presence of conserved cis-regulatory elements in the evolutionary distant CesA promoters of Arabidopsis, Populus and Eucalyptus suggests an ancient transcriptional network regulating cellulose biosynthesis in vascular plants.

The Gibberella fujikuroi complex includes many plant pathogens of agricultural crops and trees, all of which have anamorphs assigned to the genus Fusarium. In this study, an interspecific hybrid cross between Gibberella circinata and Gibberella subglutinans was used to compile a genetic linkage map. A framework map was constructed using a total of 578 AFLP markers together with the mating type (MAT-1 and MAT-2) genes and the histone (H3) gene. Twelve major linkage groups were identified (n=12). Fifty percent of the markers showed significant deviation from the expected 1:1 transmission ratio in a haploid F(1) cross (P <0.05). The transmission of the markers on the linkage map was biased towards alleles of the G. subglutinans parent, with an estimated 60% of the genome of F(1) individuals contributed by this parent. This map will serve as a powerful tool to study the genetic architecture of interspecific differentiation and pathogenicity in the two parental genomes.

Knowledge of genetic diversity (GD) and relationships among maize inbred lines is indispensable in a breeding program. Our objectives were to (1) investigate the level of genetic diversity among maize inbred lines and (2) assess their genetic structures by applying simple sequence repeat (SSR) markers. Fifty-six highland and mid-altitude maize inbred lines obtained from CIMMYT programs in Ethiopia and Zimbabwe were genotyped using 27 SSR loci. All of the genotypes studied could unequivocally be distinguished with the combination of the SSRs used. In total, 104 SSR alleles were identified, with a mean of 3.85 alleles per locus. The average polymorphism information content (PIC) was 0.58. GD expressed as Euclidean distance, varied from 0.28 to 0.73 with an average of 0.59. Cluster analysis using unweighted pair group method with arithmetic average (UPGMA) suggested five groups among the inbred lines. Most of the inbred lines adapted to the highlands and the mid-altitudes were positioned in different clusters with a few discrepancies. The pattern of groupings of the inbred lines was mostly consistent with available pedigree information. The variability detected using SSR markers could potentially contribute towards effective utilization of the inbred lines for the exploitation of heterosis and formation of genetically diverse source populations in Ethiopian maize improvement programs.

Higher plants contain a family of cellulose synthase catalytic subunit (CesA) genes that encode components of an enzyme complex embedded in the cell membrane. Recent studies in several higher plant species have demonstrated that two groups of CesA genes exist, associated with either primary or secondary cell wall deposition. We cloned six full-length CesA cDNAs from Eucalyptus grandis W. Hill ex Maiden (EgCesA1 through 6) and determined their expression patterns in a variety of organs from an adult tree. The six EgCesA genes encode predicted proteins of 978 to 1097 amino acid residues, each of which contains all of the key regions and motifs characteristic of functional CESA proteins. The predicted proteins share limited amino acid identity with each other, ranging from 61 to 70%. In contrast, similar CESA proteins from higher plant species exhibit 81 to 90% identity with the six EgCESAs. Gene expression analysis using quantitative reverse-transcription polymerase chain reaction indicated that transcripts of EgCesA1 to 3 were abundant in tissues enriched with cells laying down secondary cell walls (e.g., xylem), but were weakly expressed in tissues undergoing primary growth (e.g., unfolding leaves). Expression of EgCesA4 and EgCesA5 was upregulated in tissues rich in rapidly dividing cells undergoing primary wall synthesis, whereas EgCesA6 was weakly expressed in all tissues analyzed. These results suggest that Eucalyptus, like other higher plants, expresses two contrasting groups of apparently co-regulated CesAs involved in either primary or secondary cell wall biosynthesis.

Despite the availability of high-throughput transcript profiling technology, little is known about tissue-specific gene expression patterns in the wood-forming tissues of Eucalyptus plantation tree species. We used cDNA-amplified fragment length polymorphism (AFLP) analysis in combination with infrared fragment detection and semi-automated band quantification to profile gene expression in a 6-year-old, fast- growing Eucalyptus tree. The expression profiles of 6385 transcript-derived fragments (TDFs) were analyzed across four major woody tissues (mature xylem, immature xylem, phloem and cork) collected from two stem positions, to provide a global view of transcript abundance and variability in the Eucalyptus stem. About 21% of the TDFs were differentially expressed and could be grouped into clusters representing co- expressed genes. A total of 71 TDFs representing different gene clusters were isolated and characterized. These included genes implicated in cell fate, signal transduction and cell wall biosynthesis, processes closely associated with xylogenesis. Analysis of the expression levels of selected TDFs by quantitative RT-PCR corroborated the TDF quantification and confirmed that cDNA-AFLP analysis is a highly efficient and accurate tool for transcript profiling and gene discovery in wood-forming tissues of tree species.

Species diversity may have evolved by differential regulation of a similar set of genes. To analyze and compare the genetic architecture of transcript regulation in different genetic backgrounds of Eucalyptus, microarrays were used to examine variation in mRNA abundance in the differentiating xylem of a E. grandis pseudobackcross population [E. grandis x F(1) hybrid (E. grandis x E. globulus)]. Least-squares mean estimates of transcript levels were generated for 2608 genes in 91 interspecific backcross progeny. The quantitative measurements of variation in transcript abundance for specific genes were mapped as expression QTL (eQTL) in two single-tree genetic linkage maps (F(1) hybrid paternal and E. grandis maternal). EQTL were identified for 1067 genes in the two maps, of which 811 were located in the F(1) hybrid paternal map, and 451 in the E. grandis maternal map. EQTL for 195 genes mapped to both parental maps, the majority of which localized to nonhomologous linkage groups, suggesting trans-regulation by different loci in the two genetic backgrounds. For 821 genes, a single eQTL that explained up to 70% of the transcript-level variation was identified. Hotspots with colocalized eQTL were identified in both maps and typically contained genes associated with specific metabolic and regulatory pathways, suggesting coordinated genetic regulation.

Development of improved Eucalyptus genotypes involves the routine identification of breeding stock and superior clones. Currently, microsatellites and random amplified polymorphic DNA markers are the most widely used DNA-based techniques for fingerprinting of these trees. While these techniques have provided rapid and powerful fingerprinting assays, they are constrained by their reliance on gel or capillary electrophoresis, and therefore, relatively low throughput of fragment analysis. In contrast, recently developed microarray technology holds the promise of parallel analysis of thousands of markers in plant genomes. The aim of this study was to develop a DNA fingerprinting chip for Eucalyptus grandis and to investigate its usefulness for fingerprinting of eucalypt trees. A prototype chip was prepared using a partial genomic library from total genomic DNA of 23 E. grandis trees, of which 22 were full siblings. A total of 384 cloned genomic fragments were individually amplified and arrayed onto glass slides. DNA fingerprints were obtained for 17 individuals by hybridizing labeled genome representations of the individual trees to the 384-element chip. Polymorphic DNA fragments were identified by evaluating the binary distribution of their background-corrected signal intensities across full-sib individuals. Among 384 DNA fragments on the chip, 104 (27%) were found to be polymorphic. Hybridization of these polymorphic fragments was highly repeatable (R2>0.91) within the E. grandis individuals, and they allowed us to identify all 17 full-sib individuals. Our results suggest that DNA microarrays can be used to effectively fingerprint large numbers of closely related Eucalyptus trees.

Phenotypic, genotypic, and transcript level (microarray) data from an interspecific backcross population of Eucalyptus grandis and Eucalyptus globulus were integrated to dissect the genetic and metabolic network underlying growth variation. Transcript abundance, measured for 2,608 genes in the differentiating xylem of a 91 (E. grandis x E. globulus) x E. grandis backcross progeny was correlated with diameter variation, revealing coordinated down-regulation of genes encoding enzymes of the lignin biosynthesis and associated methylation pathways in fast growing individuals. Lignin analysis of wood samples confirmed the content and quality predicted by the transcript levels measured on the microarrays. Quantitative trait locus (QTL) analysis of transcript levels of lignin-related genes showed that their mRNA abundance is regulated by two genetic loci, demonstrating coordinated genetic control over lignin biosynthesis. These two loci colocalize with QTLs for growth, suggesting that the same genomic regions are regulating growth, and lignin content and composition in the progeny. Genetic mapping of the lignin genes revealed that most of the key biosynthetic genes do not colocalize with growth and transcript level QTLs, with the exception of the locus encoding the enzyme S-adenosylmethionine synthase. This study illustrates the power of integrating quantitative analysis of gene expression data and genetic map information to discover genetic and metabolic networks regulating complex biological traits. (Sequence data for this article have been deposited with the EMBL/GenBank data libraries under accession numbers CB 967505 to CB 968059; CD 667988 to CD 670002; CD 670004; CD 670097; CD 670101 to CD 670112; and CD 670114 to CD 670137.)

The genetic architecture of hybrid fitness characters can provide valuable insights into the nature and evolution of postzygotic reproductive barriers in diverged species. We determined the genome-wide distribution of barriers to introgression in an F(1) hybrid of two Eucalyptus tree species, Eucalyptus grandis (W. Hill ex Maiden.) and E. globulus (Labill.). Two interspecific backcross families (N = 186) were used to construct comparative, single-tree, genetic linkage maps of an F(1) hybrid individual and two backcross parents. A total of 1354 testcross AFLP marker loci were evaluated in the three parental maps and a substantial proportion (27.7% average) exhibited transmission ratio distortion (alpha = 0.05). The distorted markers were located in distinct regions of the parental maps and marker alleles within each region were all biased toward either of the two parental species. We used a Bayesian approach to estimate the position and effect of transmission ratio distorting loci (TRDLs) in the distorted regions of each parental linkage map. The relative viability of TRDL alleles ranged from 0.20 to 0.72. Contrary to expectation, heterospecific (donor) alleles of TRDLs were favored as often as recurrent alleles in both backcrosses, suggesting that positive and negative heterospecific interactions affect introgression rates in this wide interspecific pedigree.

Comparative genetic mapping in interspecific pedigrees presents a powerful approach to study genetic differentiation, genome evolution and reproductive isolation in diverging species. We used this approach for genetic analysis of an F(1) hybrid of two Eucalyptus tree species, Eucalyptus grandis (W. Hill ex Maiden.) and Eucalyptus globulus (Labill.). This wide interspecific cross is characterized by hybrid inviability and hybrid abnormality. Approximately 20% of loci in the genome of the F(1) hybrid are expected to be hemizygous due to a difference in genome size between E. grandis (640 Mbp) and E. globulus (530 Mbp). We investigated the extent of colinearity between the two genomes and the distribution of hemizygous loci in the F(1) hybrid using high-throughput, semi-automated AFLP marker analysis. Two pseudo-backcross families (backcrosses of an F(1) individual to non-parental individuals of the parental species) were each genotyped with more than 800 AFLP markers. This allowed construction of de novo comparative genetic linkage maps of the F(1) hybrid and the two backcross parents. All shared AFLP marker loci in the three single-tree parental maps were found to be colinear and little evidence was found for gross chromosomal rearrangements. Our results suggest that hemizygous AFLP loci are dispersed throughout the E. grandis chromosomes of the F(1) hybrid.

Amplified fragment length polymorphism (AFLP) analysis is currently the most powerful and efficient technique for the generation of large numbers of anonymous DNA markers in plant and animal genomes. We have developed a protocol for high-throughput AFLP analysis that allows up to 70,000 polymorphic marker genotype determinations per week on a single automated DNA sequencer. This throughput is based on multiplexed PCR amplification of AFLP fragments using two different infrared dyelabeled primer combinations. The multiplexed AFLPs are resolved on a two-dye, model 4200 LI-COR automated DNA sequencer, and the digital images are scored using semi-automated scoring software specifically designed for complex AFLP banding patterns (AFLP-Quantar). Throughput is enhanced by using high-quality genomic DNA templates obtained by a 96-well DNA isolation procedure.