Search

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

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

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

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

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

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

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

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

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.

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.