The genome of Eucalyptus grandis

Genomic determinants, architecture, and constraints in drought-related traits in Corymbia calophylla

BACKGROUND

Drought adaptation is critical to many tree species persisting under climate change, however our knowledge of the genetic basis for trees to adapt to drought is limited. This knowledge gap impedes our fundamental understanding of drought response and application to forest production and conservation. To improve our understanding of the genomic determinants, architecture, and trait constraints, we assembled a reference genome and detected ~ 6.5 M variants in 432 phenotyped individuals for the foundational tree Corymbia calophylla.

RESULTS

We found 273 genomic variants determining traits with moderate heritability (h2SNP = 0.26-0.64). Significant variants were predominantly in gene regulatory elements distributed among several haplotype blocks across all chromosomes. Furthermore, traits were constrained by frequent epistatic and pleiotropic interactions.

CONCLUSIONS

Our results on the genetic basis for drought traits in Corymbia calophylla have several implications for the ability to adapt to climate change: (1) drought related traits are controlled by complex genomic architectures with large haplotypes, epistatic, and pleiotropic interactions; (2) the most significant variants determining drought related traits occurred in regulatory regions; and (3) models incorporating epistatic interactions increase trait predictions. Our findings indicate that despite moderate heritability drought traits are likely constrained by complex genomic architecture potentially limiting trees response to climate change.

Rare but diverse off-target and somatic mutations found in field and greenhouse grown trees expressing CRISPR/Cas9

Introduction: CRISPR gene editing, while highly efficient in creating desired mutations, also has the potential to cause off-target mutations. This risk is especially high in clonally propagated plants, where editing reagents may remain in the genome for long periods of time or in perpetuity. We studied a diverse population of Populus and Eucalyptus trees that had CRISPR/Cas9-containing transgenes that targeted one or two types of floral development genes, homologs of LEAFY and AGAMOUS. Methods: Using a targeted sequence approach, we studied approximately 20,000 genomic sites with degenerate sequence homology of up to five base pairs relative to guide RNA (gRNA) target sites. We analyzed those sites in 96 individual tree samples that represented 37 independent insertion events containing one or multiples of six unique gRNAs. Results: We found low rates of off-target mutations, with rates of 1.2 × 10-9 in poplar and 3.1 × 10-10 in eucalypts, respectively, comparable to that expected due to sexual reproduction. The rates of mutation were highly idiosyncratic among sites and not predicted by sequence similarity to the target sites; a subset of two gRNAs showed off-target editing of four unique genomic sites with up to five mismatches relative to the true target sites, reaching fixation in some gene insertion events and clonal ramets. The location of off-target mutations relative to the PAM site were essentially identical to that seen with on-target CRISPR mutations. Discussion: The low rates observed support many other studies in plants that suggest that the rates of off-target mutagenesis from CRISPR/Cas9 transgenes are negligible; our study extends this conclusion to trees and other long-lived plants where CRISPR/Cas9 transgenes were present in the genome for approximately four years.

Genome-wide identification and characterization of SNW/SKIP domain-containing proteins in plants

Sessile organisms, such as plants, developed various ways to sense and respond to external and internal stimuli to maximize their fitness through evolutionary time. Transcripts and protein regulation are, among many, the main mechanisms that plants use to respond to environmental changes. SKIP protein is one such, presenting an SNKW interacting domain, which is highly conserved among eukaryotes, where SKI interacting protein acts in regulating key processes. In the present work, many bioinformatics tools, such as phylogenetic relationships, gene structure, physical-chemical properties, conserved motifs, prediction of regulatory cis-elements, chromosomal localization, and protein-protein interaction network, were used to better understand the genome-wide SNW/SKIP domain-containing proteins. In total, 28 proteins containing the SNW/SKIP domain were identified in different plant species, including plants of agronomic interest. Two main protein clusters were formed in phylogenetic analysis, and gene structure analysis revealed that, in general, the coding region had no introns. Also, expression of these genes is possibly induced by abiotic stress stimuli. Primary structure analysis of the proteins revealed the existence of an evolutionarily conserved functional unit. But physicochemical properties show that proteins containing the SNW/SKIP domain are commonly unstable under in vivo conditions. In addition, the protein network, demonstrated that SKIP homologues could act by modulating plant fitness through gene expression regulation at the transcriptional and post-transcriptional levels. This could be corroborated by the expression number of gene copies of SKIP proteins in many species, highlighting it's crucial role in plant development and tolerance through the course of evolution.

Genome-wide identification and expression analysis of GRAS gene family in Eucalyptus grandis

BACKGROUND

The GRAS gene family is a class of plant-specific transcription factors with important roles in many biological processes, such as signal transduction, disease resistance and stress tolerance, plant growth and development. So far, no information available describes the functions of the GRAS genes in Eucalyptus grandis.

RESULTS

A total of 82 GRAS genes were identified with amino acid lengths ranging from 267 to 817 aa, and most EgrGRAS genes had one exon. Members of the GRAS gene family of Eucalyptus grandis are divided into 9 subfamilies with different protein structures, while members of the same subfamily have similar gene structures and conserved motifs. Moreover, these EgrGRAS genes expanded primarily due to segmental duplication. In addition, cis-acting element analysis showed that this family of genes was involved involved in the signal transduction of various plant hormones, growth and development, and stress response. The qRT-PCR data indicated that 18 EgrGRAS genes significantly responded to hormonal and abiotic stresses. Among them, the expression of EgrGRAS13, EgrGRAS68 and EgrGRAS55 genes was significantly up-regulated during the treatment period, and it was hypothesised that members of the EgrGRAS family play an important role in stress tolerance.

CONCLUSIONS

In this study, the phylogenetic relationship, conserved domains, cis-elements and expression patterns of GRAS gene family of Eucalyptus grandis were analyzed, which filled the gap in the identification of GRAS gene family of Eucalyptus grandis and laid the foundation for analyzing the function of EgrGRAS gene in hormone and stress response.

An improved genome assembly of Chrysanthemum nankingense reveals expansion and functional diversification of terpene synthase gene family

BACKGROUND

Terpenes are important components of plant aromas, and terpene synthases (TPSs) are the key enzymes driving terpene diversification. In this study, we characterized the volatile terpenes in five different Chrysanthemum nankingense tissues. In addition, genome-wide identification and expression analysis of TPS genes was conducted utilizing an improved chromosome-scale genome assembly and tissue-specific transcriptomes. The biochemical functions of three representative TPSs were also investigated.

RESULTS

We identified tissue-specific volatile organic compound (VOC) and volatile terpene profiles. The improved Chrysanthemum nankingense genome assembly was high-quality, including a larger assembled size (3.26 Gb) and a better contig N50 length (3.18 Mb) compared to the old version. A total of 140 CnTPS genes were identified, with the majority representing the TPS-a and TPS-b subfamilies. The chromosomal distribution of these TPS genes was uneven, and 26 genes were included in biosynthetic gene clusters. Closely-related Chrysanthemum taxa were also found to contain diverse TPS genes, and the expression profiles of most CnTPSs were tissue-specific. The three investigated CnTPS enzymes exhibited versatile activities, suggesting multifunctionality.

CONCLUSIONS

We systematically characterized the structure and diversity of TPS genes across the Chrysanthemum nankingense genome, as well as the potential biochemical functions of representative genes. Our results provide a basis for future studies of terpene biosynthesis in chrysanthemums, as well as for the breeding of improved chrysanthemum varieties.

Comparative genomic analysis reveals expansion of the DnaJ gene family in Lagerstroemia indica and its members response to salt stress

DnaJs/Hsp40s/JPDs are obligate co-chaperones of heat shock proteins (Hsp70), performing crucial biological functions within organisms. A comparative genome analysis of four genomes (Vitis vinifera, Eucalyptus grandis, Lagerstroemia indica, and Punica granatum) revealed that the DnaJ gene family in L. indica has undergone expansion, although not to the extent observed in P. granatum. Inter-genome collinearity analysis of four plants indicates that members belonging to Class A and B are more conserved during evolution. In L. indica, the expanded members primarily belong to Class-C. Tissue expression patterns and the biochemical characterization of LiDnaJs further suggested that DnaJs may be involved in numerous biological processes in L. indica. Transcriptome and qPCR analyses of salt stressed leaves identified at least ten LiDnaJs that responded to salt stress. In summary, we have elucidated the expansion mechanism of the LiDnaJs, which is attributed to a recent whole-genome triplication. This research laid the foundation for functional analysis of LiDnaJs and provides gene resources for breeding salt-tolerant varieties of L. indica.

De novo genome assembly and population genomics of a shrub tree Barthea barthei (Hance) krass provide insights into the adaptive color variations

Flower color is a classic example of an ecologically important trait under selection in plants. Understanding the genetic mechanisms underlying shifts in flower color can provide key insights into ecological speciation. In this study, we investigated the genetic basis of flower color divergence in Barthea barthei, a shrub tree species exhibiting natural variation in flower color. We assembled a high-quality genome assembly for B. barthei with a contig N50 of 2.39 Mb and a scaffold N50 of 16.21 Mb. The assembly was annotated with 46,430 protein-coding genes and 1,560 non-coding RNAs. Genome synteny analysis revealed two recent tetraploidization events in B. barthei, estimated to have occurred at approximately 17 and 63 million years ago. These tetraploidization events resulted in massive duplicated gene content, with over 70% of genes retained in collinear blocks. Gene family members of the core regulators of the MBW complex were significantly expanded in B. barthei compared to Arabidopsis, suggesting that these duplications may have provided raw genetic material for the evolution of novel regulatory interactions and the diversification of anthocyanin pigmentation. Transcriptome profiling of B. barthei flowers revealed differential expression of 9 transcription factors related to anthocyanin biosynthesis between the two ecotypes. Six of these differentially expressed transcription factors were identified as high-confidence candidates for adaptive evolution based on positive selection signals. This study provides insights into the genetic basis of flower color divergence and the evolutionary mechanisms underlying ecological adaptation in plants.

Fungal metabolism and free amino acid content may predict nitrogen transfer to the host plant in the ectomycorrhizal relationship between Pisolithus spp. and Eucalyptus grandis

Ectomycorrhizal (ECM) fungi are crucial for tree nitrogen (N) nutrition; however, mechanisms governing N transfer from fungal tissues to the host plant are not well understood. ECM fungal isolates, even from the same species, vary considerably in their ability to support tree N nutrition, resulting in a range of often unpredictable symbiotic outcomes. In this study, we used isotopic labelling to quantify the transfer of N to the plant host by isolates from the ECM genus Pisolithus, known to have significant variability in colonisation and transfer of nutrients to a host. We considered the metabolic fate of N acquired by the fungi and found that the percentage of plant N acquired through symbiosis significantly correlated to the concentration of free amino acids in ECM extra-radical mycelium. Transcriptomic analyses complemented these findings with isolates having high amino acid content and N transfer showing increased expression of genes related to amino acid transport and catabolic pathways. These results suggest that fungal N metabolism impacts N transfer to the host plant in this interaction and that relative N transfer may be possible to predict through basic biochemical analyses.

High-density genetic linkage mapping in Sitka spruce advances the integration of genomic resources in conifers

In species with large and complex genomes such as conifers, dense linkage maps are a useful resource for supporting genome assembly and laying the genomic groundwork at the structural, populational, and functional levels. However, most of the 600+ extant conifer species still lack extensive genotyping resources, which hampers the development of high-density linkage maps. In this study, we developed a linkage map relying on 21,570 single nucleotide polymorphism (SNP) markers in Sitka spruce (Picea sitchensis [Bong.] Carr.), a long-lived conifer from western North America that is widely planted for productive forestry in the British Isles. We used a single-step mapping approach to efficiently combine RAD-seq and genotyping array SNP data for 528 individuals from 2 full-sib families. As expected for spruce taxa, the saturated map contained 12 linkages groups with a total length of 2,142 cM. The positioning of 5,414 unique gene coding sequences allowed us to compare our map with that of other Pinaceae species, which provided evidence for high levels of synteny and gene order conservation in this family. We then developed an integrated map for P. sitchensis and Picea glauca based on 27,052 markers and 11,609 gene sequences. Altogether, these 2 linkage maps, the accompanying catalog of 286,159 SNPs and the genotyping chip developed, herein, open new perspectives for a variety of fundamental and more applied research objectives, such as for the improvement of spruce genome assemblies, or for marker-assisted sustainable management of genetic resources in Sitka spruce and related species.

Sesquiterpenes of the ectomycorrhizal fungus Pisolithus microcarpus alter root growth and promote host colonization

Trees form symbioses with ectomycorrhizal (ECM) fungi, maintained in part through mutual benefit to both organisms. Our understanding of the signaling events leading to the successful interaction between the two partners requires further study. This is especially true for understanding the role of volatile signals produced by ECM fungi. Terpenoids are a predominant class of volatiles produced by ECM fungi. While several ECM genomes are enriched in the enzymes responsible for the production of these volatiles (i.e., terpene synthases (TPSs)) when compared to other fungi, we have limited understanding of the biochemical products associated with each enzyme and the physiological impact of specific terpenes on plant growth. Using a combination of phylogenetic analyses, RNA sequencing, and functional characterization of five TPSs from two distantly related ECM fungi (Laccaria bicolor and Pisolithus microcarpus), we investigated the role of these secondary metabolites during the establishment of symbiosis. We found that despite phylogenetic divergence, these TPSs produced very similar terpene profiles. We focused on the role of P. microcarpus terpenes and found that the fungus expressed a diverse array of mono-, di-, and sesquiterpenes prior to contact with the host. However, these metabolites were repressed following physical contact with the host Eucalyptus grandis. Exposure of E. grandis to heterologously produced terpenes (enriched primarily in γ -cadinene) led to a reduction in the root growth rate and an increase in P. microcarpus-colonized root tips. These results support a very early putative role of fungal-produced terpenes in the establishment of symbiosis between mycorrhizal fungi and their hosts.

Profiling of the gene expression and alternative splicing landscapes of Eucalyptus grandis

Eucalyptus is a widely planted hardwood tree species due to its fast growth, superior wood properties and adaptability. However, the post-transcriptional regulatory mechanisms controlling tissue development and stress responses in Eucalyptus remain poorly understood. In this study, we performed a comprehensive analysis of the gene expression profile and the alternative splicing (AS) landscape of E. grandis using strand-specific RNA-Seq, which encompassed 201 libraries including different organs, developmental stages, and environmental stresses. We identified 10 416 genes (33.49%) that underwent AS, and numerous differentially expressed and/or differential AS genes involved in critical biological processes, such as primary-to-secondary growth transition of stems, adventitious root formation, aging and responses to phosphorus- or boron-deficiency. Co-expression analysis of AS events and gene expression patterns highlighted the potential upstream regulatory role of AS events in multiple processes. Additionally, we highlighted the lignin biosynthetic pathway to showcase the potential regulatory functions of AS events in the KNAT3 and IRL3 genes within this pathway. Our high-quality expression atlas and AS landscape serve as valuable resources for unravelling the genetic control of woody plant development, long-term adaptation, and understanding transcriptional diversity in Eucalyptus. Researchers can conveniently access these resources through the interactive ePlant browser (https://bar.utoronto.ca/eplant_eucalyptus).

Plasticity of repetitive sequences demonstrated by the complete mitochondrial genome of Eucalyptus camaldulensis

The tree Eucalyptus camaldulensis is a ubiquitous member of the Eucalyptus genus, which includes several hundred species. Despite the extensive sequencing and assembly of nuclear genomes from various eucalypts, the genus has only one fully annotated and complete mitochondrial genome (mitogenome). Plant mitochondria are characterized by dynamic genomic rearrangements, facilitated by repeat content, a feature that has hindered the assembly of plant mitogenomes. This complexity is evident in the paucity of available mitogenomes. This study, to the best of our knowledge, presents the first E. camaldulensis mitogenome. Our findings suggest the presence of multiple isomeric forms of the E. camaldulensis mitogenome and provide novel insights into minor rearrangements triggered by nested repeat sequences. A comparative sequence analysis of the E. camaldulensis and E. grandis mitogenomes unveils evolutionary changes between the two genomes. A significant divergence is the evolution of a large repeat sequence, which may have contributed to the differences observed between the two genomes. The largest repeat sequences in the E. camaldulensis mitogenome align well with significant yet unexplained structural variations in the E. grandis mitogenome, highlighting the adaptability of repeat sequences in plant mitogenomes.

Comparison of ddRADseq and EUChip60K SNP genotyping systems for population genetics and genomic selection in Eucalyptus dunnii (Maiden)

Eucalyptus dunnii is one of the most important Eucalyptus species for short-fiber pulp production in regions where other species of the genus are affected by poor soil and climatic conditions. In this context, E. dunnii holds promise as a resource to address and adapt to the challenges of climate change. Despite its rapid growth and favorable wood properties for solid wood products, the advancement of its improvement remains in its early stages. In this work, we evaluated the performance of two single nucleotide polymorphism, (SNP), genotyping methods for population genetics analysis and Genomic Selection in E. dunnii. Double digest restriction-site associated DNA sequencing (ddRADseq) was compared with the EUChip60K array in 308 individuals from a provenance-progeny trial. The compared SNP set included 8,011 and 19,008 informative SNPs distributed along the 11 chromosomes, respectively. Although the two datasets differed in the percentage of missing data, genome coverage, minor allele frequency and estimated genetic diversity parameters, they revealed a similar genetic structure, showing two subpopulations with little differentiation between them, and low linkage disequilibrium. GS analyses were performed for eleven traits using Genomic Best Linear Unbiased Prediction (GBLUP) and a conventional pedigree-based model (ABLUP). Regardless of the SNP dataset, the predictive ability (PA) of GBLUP was better than that of ABLUP for six traits (Cellulose content, Total and Ethanolic extractives, Total and Klason lignin content and Syringyl and Guaiacyl lignin monomer ratio). When contrasting the SNP datasets used to estimate PAs, the GBLUP-EUChip60K model gave higher and significant PA values for six traits, meanwhile, the values estimated using ddRADseq gave higher values for three other traits. The PAs correlated positively with narrow sense heritabilities, with the highest correlations shown by the ABLUP and GBLUP-EUChip60K. The two genotyping methods, ddRADseq and EUChip60K, are generally comparable for population genetics and genomic prediction, demonstrating the utility of the former when subjected to rigorous SNP filtering. The results of this study provide a basis for future whole-genome studies using ddRADseq in non-model forest species for which SNP arrays have not yet been developed.

Genetically engineered eucalyptus expressing pesticidal proteins from Bacillus thuringiensis for insect resistance: a risk assessment evaluation perspective

Eucalyptus covers approximately 7.5 million hectares in Brazil and serves as the primary woody species cultivated for commercial purposes. However, native insects and invasive pests pose a significant threat to eucalyptus trees, resulting in substantial economic losses and reduced forest productivity. One of the primary lepidopteran pests affecting eucalyptus is Thyrinteina arnobia (Stoll, 1782) (Lepidoptera: Geometridae), commonly referred to as the brown looper caterpillar. To address this issue, FuturaGene, the biotech division of Suzano S.A., has developed an insect-resistant (IR) eucalyptus variety, which expresses Cry pesticidal proteins (Cry1Ab, Cry1Bb, and Cry2Aa), derived from Bacillus thuringiensis (Bt). Following extensive safety assessments, including field trials across various biomes in Brazil, the Brazilian National Technical Commission of Biosafety (CTNBio) recently approved the commercialization of IR eucalyptus. The biosafety assessments involved the analysis of molecular genomics, digestibility, thermostability, non-target organism exposure, degradability in the field, and effects on soil microbial communities and arthropod communities. In addition, in silico studies were conducted to evaluate allergenicity and toxicity. Results from both laboratory and field studies indicated that Bt eucalyptus is as safe as the conventional eucalyptus clone for humans, animals, and the environment, ensuring the secure use of this insect-resistant trait in wood production.

Genomic and transcriptomic studies on flavonoid biosynthesis in Lagerstroemia indica

BACKGROUND

Lagerstroemia indica is a widely cultivated ornamental woody shrub/tree of the family Lythraceae that is used as a traditional medicinal plant in East Asia and Egypt. However, unlike other ornamental woody plants, its genome is not well-investigated, which hindered the discovery of the key genes that regulate important traits and the synthesis of bioactive compounds.

RESULTS

In this study, the genomic sequences of L. indica were determined using several next-generation sequencing technologies. Altogether, 324.01 Mb sequences were assembled and 98.21% (318.21 Mb) of them were placed in 24 pseudo-chromosomes. The heterozygosity, repeated sequences, and GC residues occupied 1.65%, 29.17%, and 38.64% of the genome, respectively. In addition, 28,811 protein-coding gene models, 327 miRNAs, 552 tRNAs, 214 rRNAs, and 607 snRNAs were identified. The intra- and interspecies synteny and Ks analysis revealed that L. indica exhibits a hexaploidy. The co-expression profiles of the genes involved in the phenylpropanoid (PA) and flavonoid/anthocyanin (ABGs) pathways with the R2R3 MYB genes (137 members) showed that ten R2R3 MYB genes positively regulate flavonoid/anthocyanin biosynthesis. The colors of flowers with white, purple (PB), and deep purplish pink (DPB) petals were found to be determined by the levels of delphinidin-based (Dp) derivatives. However, the substrate specificities of LiDFR and LiOMT probably resulted in the different compositions of flavonoid/anthocyanin. In L. indica, two LiTTG1s (LiTTG1-1 and LiTTG1-2) were found to be the homologs of AtTTG1 (WD40). LiTTG1-1 was found to repress anthocyanin biosynthesis using the tobacco transient transfection assay.

CONCLUSIONS

This study showed that the ancestor L. indica experienced genome triplication approximately 38.5 million years ago and that LiTTG1-1 represses anthocyanin biosynthesis. Furthermore, several genes such as LiDFR, LiOMTs, and R2R3 LiMYBs are related to anthocyanin biosynthesis. Further studies are required to clarify the mechanisms and alleles responsible for flower color development.

Wood of trees: Cellular structure, molecular formation, and genetic engineering

Wood is an invaluable asset to human society due to its renewable nature, making it suitable for both sustainable energy production and material manufacturing. Additionally, wood derived from forest trees plays a crucial role in sequestering a significant portion of the carbon dioxide fixed during photosynthesis by terrestrial plants. Nevertheless, with the expansion of the global population and ongoing industrialization, forest coverage has been substantially decreased, resulting in significant challenges for wood production and supply. Wood production practices have changed away from natural forests toward plantation forests. Thus, understanding the underlying genetic mechanisms of wood formation is the foundation for developing high-quality, fast-growing plantation trees. Breeding ideal forest trees for wood production using genetic technologies has attracted the interest of many. Tremendous studies have been carried out in recent years on the molecular, genetic, and cell-biological mechanisms of wood formation, and considerable progress and findings have been achieved. These studies and findings indicate enormous possibilities and prospects for tree improvement. This review will outline and assess the cellular and molecular mechanisms of wood formation, as well as studies on genetically improving forest trees, and address future development prospects.

Genomic evidence for rediploidization and adaptive evolution following the whole-genome triplication

Whole-genome duplication (WGD), or polyploidy, events are widespread and significant in the evolutionary history of angiosperms. However, empirical evidence for rediploidization, the major process where polyploids give rise to diploid descendants, is still lacking at the genomic level. Here we present chromosome-scale genomes of the mangrove tree Sonneratia alba and the related inland plant Lagerstroemia speciosa. Their common ancestor has experienced a whole-genome triplication (WGT) approximately 64 million years ago coinciding with a period of dramatic global climate change. Sonneratia, adapting mangrove habitats, experienced extensive chromosome rearrangements post-WGT. We observe the WGT retentions display sequence and expression divergence, suggesting potential neo- and sub-functionalization. Strong selection acting on three-copy retentions indicates adaptive value in response to new environments. To elucidate the role of ploidy changes in genome evolution, we improve a model of the polyploidization-rediploidization process based on genomic evidence, contributing to the understanding of adaptive evolution during climate change.

Transcriptome analysis of the transition from primary to secondary growth of vertical stem in Eucalyptus grandis

Eucalyptus was one of the most cultivated hardwood species worldwide, with rapid growth, good wood properties and a wide range of adaptability. Eucalyptus stem undergoes primary growth (longitudinal growth) followed by secondary growth (radial growth), which produces biomass that is an important source of energy worldwide. In order to better understand the genetic regulation of secondary growth in Eucalyptus grandis, Transcriptome analyses in stem segments along a developmental gradient from the third internode to the eleventh internode of E. grandis that spanned primary to secondary growth were carried out. 5,149 genes that were differentially expressed during stem development were identified. Combining the trend analysis by the Mfuzz method and the module-trait correlation analysis by the Weighted Gene Co-expression Network Analysis method, a total of 70 differentially expressed genes (DEGs) selected from 868 DEGs with high connectivity were found to be closely correlated with secondary growth. Results revealed that the differential expression of these DEGs suggests that they may involve in the primary growth or secondary growth. AP1, YAB2 TFs and EXP genes are highly expressed in the IN3, whereas NAC, MYB TFs are likely to be important for secondary growth. These results will expand our understanding of the complex molecular and cellular events of secondary growth and provide a foundation for future studies on wood formation in Eucalyptus.

Cytokinin promotes anthocyanin biosynthesis via regulating sugar accumulation and MYB113 expression in Eucalyptus

Anthocyanins are flavonoid-like substances that play important roles in plants' adaptation to various environmental stresses. In this research, we discovered that cytokinin (CK) alone could effectively induce the anthocyanin biosynthesis in Eucalyptus and many other perennial woody plant species, but not in tobacco and Arabidopsis, suggesting a diverse role of CK in regulating anthocyanin biosynthesis in different species. Transcriptomic and metabolomic strategies were used to further clarify the specific role of CK in regulating anthocyanin biosynthesis in Eucalyptus. The results showed that 801 and 2241 genes were differentially regulated at 6 and 24 h, respectively, after CK treatment. Pathway analysis showed that most of the differentially expressed genes were categorized into pathways related to cellular metabolism or transport of metabolites, including amino acids and sugars. The metabolomic results well supported the transcriptome data, which showed that most of the differentially regulated metabolites were related to the metabolism of sugar, amino acids and flavonoids. Moreover, CK treatment significantly induced the accumulation of sucrose in the CK-treated leaves, while sugar starvation mimicked by either defoliation or shading treatment of the basal leaves significantly reduced the sugar increase of the CK-treated leaves and thus inhibited CK-induced anthocyanin biosynthesis. The results of in vitro experiment also suggested that CK-induced anthocyanin in Eucalyptus was sugar-dependent. Furthermore, we identified an early CK-responsive transcription factor MYB113 in Eucalyptus, the expression of which was significantly upregulated by CK treatment in Eucalyptus, but was inhibited in Arabidopsis. Importantly, the overexpression of EgrMYB113 in the Eucalyptus hairy roots was associated with significant anthocyanin accumulation and upregulation of most of the anthocyanin biosynthetic genes. In conclusion, our study demonstrates a key role of CK in the regulation of anthocyanin biosynthesis in Eucalyptus, providing a molecular basis for further understanding the regulatory mechanism and diversity of hormone-regulated anthocyanin biosynthesis in different plant species.

Reference genes for Eucalyptus spp. under Beauveria bassiana inoculation and subsequently infestation by the galling wasp Leptocybe invasa

Relative gene expression analysis through RT-qPCR is an important molecular technique that helps understanding different molecular mechanisms, such as the plant defense response to insect pests. However, the use of RT-qPCR for gene expression analysis can be affected by factors that directly affect the reliability of the results. Among these factors, the appropriate choice of reference genes is crucial and can strongly impact RT-qPCR relative gene expression analyses, highlighting the importance in correctly choosing the most suitable genes for the success of the analysis. Thus, this study aimed to select and validate reference genes for relative gene expression studies through RT-qPCR in hybrids of Eucalyptus tereticornis × Eucalyptus camaldulensis (drought tolerant and susceptible to Leptocybe invasa) under conditions of inoculation by the Beauveria bassiana fungus and subsequent infestation by L. invasa. The expression level and stability of eleven candidate genes were evaluated. Stability was analyzed using the RefFinder tool, which integrates the geNorm, NormFinder, BestKeeper, and Delta-Ct algorithms. The selected reference genes were validated through the expression analysis of the transcriptional factor EcDREB2 (dehydration-responsive element-binding protein 2). For all treatments evaluated, EcPTB, EcPP2A-1, and EcEUC12 were the best reference genes. The triplets EcPTB/EcEUC12/EcUBP6, EcPP2A-1/EcEUC12/EcPTB, EcIDH/EcSAND/Ecα-TUB, EcPP2A-1/Ecα-TUB/EcPTB, and EcPP2A-1/EcUPL7/EcSAND were the best reference genes for the control plants, mother plants, plants inoculated with B. bassiana, plants infested with L. invasa, and plants inoculated with B. bassiana and subsequently infested with L. invasa, respectively. The best determined reference genes were used to normalize the RT-qPCR expression data for each experimental condition evaluated. The results emphasize the importance of this type of study to ensure the reliability of relative gene expression analyses. Furthermore, the findings of this study can be used as a basis for future research, comprising gene expression analysis of different eucalyptus metabolic pathways.

Phylogenetic and functional analysis of tiller angle control homeologs in allotetraploid cotton

Introduction

Plants can adapt their growth to optimize light capture in competitive environments, with branch angle being a crucial factor influencing plant phenotype and physiology. Decreased branch angles in cereal crops have been shown to enhance productivity in high-density plantings. The Tiller Angle Control (TAC1) gene, known for regulating tiller inclination in rice and corn, has been found to control branch angle in eudicots. Manipulating TAC1 in field crops like cotton offers the potential for improving crop productivity.

Methods

Using a homolog-based methodology, we examined the distribution of TAC1-related genes in cotton compared to other angiosperms. Furthermore, tissue-specific qPCR analysis unveiled distinct expression patterns of TAC1 genes in various cotton tissues. To silence highly expressed specific TAC1 homeologs in the stem, we applied CRISPR-Cas9 gene editing and Agrobacterium-mediated transformation, followed by genotyping and subsequent phenotypic validation of the mutants.

Results

Gene duplication events of TAC1 specific to the Gossypium lineage were identified, with 3 copies in diploid progenitors and 6 copies in allotetraploid cottons. Sequence analysis of the TAC1 homeologs in Gossypium hirsutum revealed divergence from other angiosperms with 1-2 copies, suggesting possible neo- or sub-functionalization for the duplicated copies. These TAC1 homeologs exhibited distinct gene expression patterns in various tissues over developmental time, with elevated expression of A11G109300 and D11G112200, specifically in flowers and stems, respectively. CRISPR-mediated loss of these TAC1 homeologous genes resulted in a reduction in branch angle and altered petiole angles, and a 5 to 10-fold reduction in TAC1 expression in the mutants, confirming their role in controlling branch and petiole angles. This research provides a promising strategy for genetically engineering branch and petiole angles in commercial cotton varieties, potentially leading to increased productivity.

E. globulus leaf EO exhibits anti-inflammatory effects by regulating GSDMD-mediated pyroptosis, thereby alleviating neurological impairment and neuroinflammation in experimental stroke mice

ETHNOPHARMACOLOGICAL RELEVANCE

Aromatic and medicinal plants continue to be a major component of alternative and traditional medicine in the developing countries. Eucalyptus globulus (Labill.) is being employed to cultivation and production in China. However, few studies have reported the chemical composition and anti-inflammatory activity of Eucalyptus globulus (Labill.) leaf essential oil (E. globulus leaf EO) extracted from Eucalyptus globulus.

AIM OF THE STUDY

This study aimed to assess the composition of E. globulus leaf EO and identify its bacteriostatic action as well as anti-inflammatory activity. Importantly, we evaluated the effect of E. globulus leaf EO on neurological impairment and neuroinflammation in experimental stroke mice.

MATERIALS AND METHODS

Gas Chromatography-Mass Spectrometer (GC-MS) analyses was employed to evaluate the chemical components of E. globulus leaf EO, and the relative content of each component was determined by area normalization method. The antimicrobial activity of E. globulus leaf EO was determined by Oxford cup method and microbroth dilution assay. Cytotoxic activity of E. globulus leaf EO on THP-1 cells or BV2 cells in vitro was determined by CCK8 assay. In addition, the lipopolysaccharide (LPS)/ATP-induced inflammation model in THP-1 cells or BV2 cells were established, and the relative expression of TNF-α, IL-1β, MCP-1and IL-6 were confirmed by RT-PCR. Furthermore, the expression of protein GSDMD, IL-lβ, NLRP3 and NFκB signaling pathway were assessed by immunoblotting. In vivo,the experimental stroke model constructed by middle cerebral artery occlusion/reperfusion (MCAO/R) in mice was employed and subsequently treated with E. globulus leaf EO (50,100 mg/kg, subcutaneous injection) for 3 days to assess neurological impairment and neuroinflammation. Behavioral and neuronal damage were assessed using grip strength test, rod trarod test, and Nissl staining. Pro-inflammatory factors in serum or ischemic brain tissue was detected by ELISA kits.

RESULTS

GC-MS analyses revealed that the major compound in E. globulus leaf EO was eudesmol (71.967%). E. globulus leaf EO has antimicrobial activity against Staphylococcus aureus (gram positive bacteria, MIC = 0.0625 mg/mL), Escherichia coli (gram negative bacteria, MIC = 1 mg/mL), and Candida albicans (MIC = 4 mg/mL). E. globulus leaf EO (0.5312, 1.0625, and 2.15 mg/mL) significantly decreased the expression of inflammation-related genes, including IL-1β, TNF-α, MCP-1, and IL-6. Furthermore, reduced levels of TLR4, Myd88, phosphorylated NF-κB P65, and IκBα were found in the E. globulus leaf EO group for BV2 cells (1.025, and 2.125 mg/mL). In addition, the expression levels of GSDMD, NLRP3, IL-1β and AIM2 were significantly decreased in the E. globulus leaf EO group when compared with the LPS -stimulated group, regulating GSDMD-mediated pyroptosis. In vivo, E. globulus leaf EO improved neurological functional deficits, inhibited excessive activation of microglia, and reduced the secretion of pro-inflammatory factors IL-1β, TNF-α in the ischemic tissue and serum after MCAO/R.

CONCLUSION

E. globulus leaf EO has strong antibacterial and anti-inflammatory activity, which has been implicated in blocking GSDMD-mediated pyroptosis. Moreover, E. globulus leaf EO could exert neuroprotective effect on cerebral ischemia-reperfusion injury.

Genomic evolution and insights into agronomic trait innovations of Sesamum species

Sesame is an ancient oilseed crop with high oil content and quality. However, the evolutionary history and genetic mechanisms of its valuable agronomic traits remain unclear. Here, we report chromosome-scale genomes of cultivated sesame (Sesamum indicum L.) and six wild Sesamum species, representing all three karyotypes within this genus. Karyotyping and genome-based phylogenic analysis revealed the evolutionary route of Sesamum species from n = 13 to n = 16 and revealed that allotetraploidization occurred in the wild species Sesamum radiatum. Early divergence of the Sesamum genus (48.5-19.7 million years ago) during the Tertiary period and its ancient phylogenic position within eudicots were observed. Pan-genome analysis revealed 9164 core gene families in the 7 Sesamum species. These families are significantly enriched in various metabolic pathways, including fatty acid (FA) metabolism and FA biosynthesis. Structural variations in SiPT1 and SiDT1 within the phosphatidyl ethanolamine-binding protein gene family lead to the genomic evolution of plant-architecture and inflorescence-development phenotypes in Sesamum. A genome-wide association study (GWAS) of an interspecific population and genome comparisons revealed a long terminal repeat insertion and a sequence deletion in DIR genes of wild Sesamum angustifolium and cultivated sesame, respectively; both variations independently cause high susceptibility to Fusarium wilt disease. A GWAS of 560 sesame accessions combined with an overexpression study confirmed that the NAC1 and PPO genes play an important role in upregulating oil content of sesame. Our study provides high-quality genomic resources for cultivated and wild Sesamum species and insights that can improve molecular breeding strategies for sesame and other oilseed crops.

Exploring the versatility of sesquiterpene biosynthesis in guava plants: a comparative genome-wide analysis of two cultivars

Psidium guajava L., a fruit crop belonging to the Myrtaceae family, is highly valued for its nutritional and medicinal properties. The family exhibits a diverse chemical profile of essential oils and serves as a valuable resource due to its ecological interactions, adaptability, and dispersal capacity. The Myrtaceae family has been extensively studied for its terpenoids. Genetic studies have focused on foliar terpene yield in species from the Eucalypteae and Melaleucaceae tribes. To understand the evolutionary trends in guava breeding, this study predicted terpene synthase genes (TPS) from different cultivars. Through this analysis, 43 full-length TPS genes were identified, and approximately 77% of them exhibited relative expression in at least one of the five investigated plant tissues (root, leaf, bud, flower, and fruit) of two guava cultivars. We identified intra-species variation in the terpene profile and single nucleotide polymorphisms (SNPs) in twelve TPS genes, resulting in the clustering of 62 genotypes according to their essential oil chemotypes. The high concentration of sesquiterpenes is supported by the higher number of TPS-a genes and their expression. The expansion for TPS sub-families in P. guajava occurred after the expansion of other rosids species. Providing insight into the origin of structural diversification and expansion in each clade of the TPS gene family within Myrtaceae. This study can provide insights into the diversity of genes for specialized metabolites such as terpenes, and their regulation, which can lead to a diverse chemotype of essential oil in different tissues and genotypes. This suggests a mode of enzymatic evolution that could lead to high sesquiterpene production, act as a chemical defense and contribute to the adaptive capacity of this species to different habitats.

A two-sequence motif-based method for the inventory of gene families in fragmented and poorly annotated genome sequences

Databases of genome sequences are growing exponentially, but, in some cases, assembly is incomplete and genes are poorly annotated. For evolutionary studies, it is important to identify all members of a given gene family in a genome. We developed a method for identifying most, if not all, members of a gene family from raw genomes in which assembly is of low quality, using the P-type ATPase superfamily as an example. The method is based on the translation of an entire genome in all six reading frames and the co-occurrence of two family-specific sequence motifs that are in close proximity to each other. To test the method's usability, we first used it to identify P-type ATPase members in the high-quality annotated genome of barley (Hordeum vulgare). Subsequently, after successfully identifying plasma membrane H+-ATPase family members (P3A ATPases) in various plant genomes of varying quality, we tested the hypothesis that the number of P3A ATPases correlates with the ability of the plant to tolerate saline conditions. In 19 genomes of glycophytes and halophytes, the total number of P3A ATPase genes was found to vary from 7 to 22, but no significant difference was found between the two groups. The method successfully identified P-type ATPase family members in raw genomes that are poorly assembled.

Exogenous applications of brassinosteroids promote secondary xylem differentiation in Eucalyptus grandis

Brassinosteroids (BRs) play many pivotal roles in plant growth and development, especially in cell elongation and vascular development. Although its biosynthetic and signal transduction pathway have been well characterized in model plants, their biological roles in Eucalyptus grandis, a major hardwood tree providing fiber and energy worldwide, remain unclear. Here, we treated E. grandis plantlets with 24-epibrassinolide (EBL), the most active BR and/or BR biosynthesis inhibitor brassinazole. We recorded the plant growth and analyzed the cell structure of the root and stem with histochemical methods; then, we performed a secondary growth, BR synthesis, and signaling-related gene expression analysis. The results showed that the BRs dramatically increased the shoot length and diameter, and the exogenous BR increased the xylem area of the stem and root. In this process, EgrBRI1, EgrBZR1, and EgrBZR2 expression were induced by the BR treatment, and the expressions of HD-ZIPIII and cellulose synthase genes were also altered. To further verify the effect of BRs in secondary xylem development in Eucalyptus, we used six-month-old plants as the material and directly applied EBL to the xylem and cambium of the vertical stems. The xylem area, fiber cell length, and cell numbers showed considerable increases. Several key BR-signaling genes, secondary xylem development-related transcription factor genes, and cellulose and lignin biosynthetic genes were also considerably altered. Thus, BR had regulatory roles in secondary xylem development and differentiation via the BR-signaling pathway in this woody plant.

Automatic annotation of the bHLH gene family in plants

BACKGROUND

The bHLH transcription factor family is named after the basic helix-loop-helix (bHLH) domain that is a characteristic element of their members. Understanding the function and characteristics of this family is important for the examination of a wide range of functions. As the availability of genome sequences and transcriptome assemblies has increased significantly, the need for automated solutions that provide reliable functional annotations is emphasised.

RESULTS

A phylogenetic approach was adapted for the automatic identification and functional annotation of the bHLH transcription factor family. The bHLH_annotator, designed for the automated functional annotation of bHLHs, was implemented in Python3. Sequences of bHLHs described in literature were collected to represent the full diversity of bHLH sequences. Previously described orthologs form the basis for the functional annotation assignment to candidates which are also screened for bHLH-specific motifs. The pipeline was successfully deployed on the two Arabidopsis thaliana accessions Col-0 and Nd-1, the monocot species Dioscorea dumetorum, and a transcriptome assembly of Croton tiglium. Depending on the applied search parameters for the initial candidates in the pipeline, species-specific candidates or members of the bHLH family which experienced domain loss can be identified.

CONCLUSIONS

The bHLH_annotator allows a detailed and systematic investigation of the bHLH family in land plant species and classifies candidates based on bHLH-specific characteristics, which distinguishes the pipeline from other established functional annotation tools. This provides the basis for the functional annotation of the bHLH family in land plants and the systematic examination of a wide range of functions regulated by this transcription factor family.

Reference genome of the nutrition-rich orphan crop chia (Salvia hispanica) and its implications for future breeding

Chia (Salvia hispanica L.) is one of the most popular nutrition-rich foods and pseudocereal crops of the family Lamiaceae. Chia seeds are a rich source of proteins, polyunsaturated fatty acids (PUFAs), dietary fibers, and antioxidants. In this study, we present the assembly of the chia reference genome, which spans 303.6 Mb and encodes 48,090 annotated protein-coding genes. Our analysis revealed that ~42% of the chia genome harbors repetitive content, and identified ~3 million single nucleotide polymorphisms (SNPs) and 15,380 simple sequence repeat (SSR) marker sites. By investigating the chia transcriptome, we discovered that ~44% of the genes undergo alternative splicing with a higher frequency of intron retention events. Additionally, we identified chia genes associated with important nutrient content and quality traits, such as the biosynthesis of PUFAs and seed mucilage fiber (dietary fiber) polysaccharides. Notably, this is the first report of in-silico annotation of a plant genome for protein-derived small bioactive peptides (biopeptides) associated with improving human health. To facilitate further research and translational applications of this valuable orphan crop, we have developed the Salvia genomics database (SalviaGDB), accessible at https://salviagdb.org.

Chromosome-scale genomes of commercially important mahoganies, Swietenia macrophylla and Khaya senegalensis

Mahogany species (family Meliaceae) are highly valued for their aesthetic and durable wood. Despite their economic and ecological importance, genomic resources for mahogany species are limited, hindering genetic improvement and conservation efforts. Here we perform chromosome-scale genome assemblies of two commercially important mahogany species: Swietenia macrophylla and Khaya senegalensis. By combining 10X sequencing and Hi-C data, we assemble high-quality genomes of 274.49 Mb (S. macrophylla) and 406.50 Mb (K. senegalensis), with scaffold N50 lengths of 8.51 Mb and 7.85 Mb, respectively. A total of 99.38% and 98.05% of the assembled sequences are anchored to 28 pseudo-chromosomes in S. macrophylla and K. senegalensis, respectively. We predict 34,129 and 31,908 protein-coding genes in S. macrophylla and K. senegalensis, respectively, of which 97.44% and 98.49% are functionally annotated. The chromosome-scale genome assemblies of these mahogany species could serve as a vital genetic resource, especially in understanding the properties of non-model woody plants. These high-quality genomes could support the development of molecular markers for breeding programs, conservation efforts, and the sustainable management of these valuable forest resources.

Pangenome of water caltrop reveals structural variations and asymmetric subgenome divergence after allopolyploidization

Water caltrop (Trapa spp., Lythraceae) is a traditional but currently underutilized non-cereal crop. Here, we generated chromosome-level genome assemblies for the two diploid progenitors of allotetraploid Trapa. natans (4x, AABB), i.e., diploid T. natans (2x, AA) and Trapa incisa (2x, BB). In conjunction with four published (sub)genomes of Trapa, we used gene-based and graph-based pangenomic approaches and a pangenomic transposable element (TE) library to develop Trapa genomic resources. The pangenome displayed substantial gene-content variation with dispensable and private gene clusters occupying a large proportion (51.95%) of the total cluster sets in the six (sub)genomes. Genotyping of presence-absence variation (PAVs) identified 40 453 PAVs associated with 2570 genes specific to A- or B-lineages, of which 1428 were differentially expressed, and were enriched in organ development process, organic substance metabolic process and response to stimulus. Comparative genome analyses showed that the allotetraploid T. natans underwent asymmetric subgenome divergence, with the B-subgenome being more dominant than the A-subgenome. Multiple factors, including PAVs, asymmetrical amplification of TEs, homeologous exchanges (HEs), and homeolog expression divergence, together affected genome evolution after polyploidization. Overall, this study sheds lights on the genome architecture and evolution of Trapa, and facilitates its functional genomic studies and breeding program.

Chromosome-scale assemblies of S. malaccense, S. aqueum, S. jambos, and S. syzygioides provide insights into the evolution of Syzygium genomes

Syzygium is a large and diverse tree genus in the Myrtaceae family. Genome assemblies for clove (Syzygium aromaticum, 370 Mb) and sea apple (Syzygium grande, 405 Mb) provided the first insights into the genomic features and evolution of the Syzygium genus. Here, we present additional de novo chromosome-scale genome assemblies for Syzygium malaccense, Syzygium aqueum, Syzygium jambos, and Syzygium syzygioides. Genome profiling analyses show that S. malaccense, like S. aromaticum and S. grande, is diploid (2n = 2x = 22), while the S. aqueum, S. jambos, and S. syzygioides specimens are autotetraploid (2n = 4x = 44). The genome assemblies of S. malaccense (430 Mb), S. aqueum (392 Mb), S. jambos (426 Mb), and S. syzygioides (431 Mb) are highly complete (BUSCO scores of 98%). Comparative genomics analyses showed conserved organization of the 11 chromosomes with S. aromaticum and S. grande, and revealed species-specific evolutionary dynamics of the long terminal repeat retrotransposon elements belonging to the Gypsy and Copia lineages. This set of Syzygium genomes is a valuable resource for future structural and functional comparative genomic studies on Myrtaceae species.

Functional investigation of five R2R3-MYB transcription factors associated with wood development in Eucalyptus using DAP-seq-ML

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

Cell-type-specific control of secondary cell wall formation by Musashi-type translational regulators in Arabidopsis

Deciphering the mechanism of secondary cell wall/SCW formation in plants is key to understanding their development and the molecular basis of biomass recalcitrance. Although transcriptional regulation is essential for SCW formation, little is known about the implication of post-transcriptional mechanisms in this process. Here we report that two bonafide RNA-binding proteins homologous to the animal translational regulator Musashi, MSIL2 and MSIL4, function redundantly to control SCW formation in Arabidopsis. MSIL2/4 interactomes are similar and enriched in proteins involved in mRNA binding and translational regulation. MSIL2/4 mutations alter SCW formation in the fibers, leading to a reduction in lignin deposition, and an increase of 4-O-glucuronoxylan methylation. In accordance, quantitative proteomics of stems reveal an overaccumulation of glucuronoxylan biosynthetic machinery, including GXM3, in the msil2/4 mutant stem. We showed that MSIL4 immunoprecipitates GXM mRNAs, suggesting a novel aspect of SCW regulation, linking post-transcriptional control to the regulation of SCW biosynthesis genes.

Exploring Pan-Genomes: An Overview of Resources and Tools for Unraveling Structure, Function, and Evolution of Crop Genes and Genomes

The availability of multiple sequenced genomes from a single species made it possible to explore intra- and inter-specific genomic comparisons at higher resolution and build clade-specific pan-genomes of several crops. The pan-genomes of crops constructed from various cultivars, accessions, landraces, and wild ancestral species represent a compendium of genes and structural variations and allow researchers to search for the novel genes and alleles that were inadvertently lost in domesticated crops during the historical process of crop domestication or in the process of extensive plant breeding. Fortunately, many valuable genes and alleles associated with desirable traits like disease resistance, abiotic stress tolerance, plant architecture, and nutrition qualities exist in landraces, ancestral species, and crop wild relatives. The novel genes from the wild ancestors and landraces can be introduced back to high-yielding varieties of modern crops by implementing classical plant breeding, genomic selection, and transgenic/gene editing approaches. Thus, pan-genomic represents a great leap in plant research and offers new avenues for targeted breeding to mitigate the impact of global climate change. Here, we summarize the tools used for pan-genome assembly and annotations, web-portals hosting plant pan-genomes, etc. Furthermore, we highlight a few discoveries made in crops using the pan-genomic approach and future potential of this emerging field of study.

High-quality genome of a pioneer mangrove Laguncularia racemosa explains its advantages for intertidal zone reforestation

Ecological restoration of mangrove ecosystems that became susceptible to recent habitat perturbations is crucial for tropical coast conservation. The white mangrove Laguncularia racemosa, a pioneer species inhabiting intertidal environments of the Atlantic East Pacific (AEP) region, has been used for reforestation in China for decades. However, the molecular mechanisms underlying its fast growth and high adaptive potential remain unknown. Using PacBio single-molecule real-time sequencing, we completed a high-quality L. racemosa genome assembly covering 1105 Mb with scaffold N50 of 3.46 Mb. Genomic phylogeny shows that L. racemosa invaded intertidal zones during a period of global warming. Multi-level genomic convergence analyses between L. racemosa and three native dominant mangrove clades show that they experienced convergent changes in genes involved in nutrient absorption and high salinity tolerance. This may explain successful L. racemosa adaptation to stressful intertidal environments after introduction. Without recent whole-genome duplications or activated transposable elements, L. racemosa has retained many tandem gene duplications. Some of them are involved in auxin biosynthesis, intense light stress and cold stress response pathways, associated with L. racemosa's ability to grow fast under high light or cold conditions when used for reforestation. In summary, our study identifies shared mechanisms of intertidal environmental adaptation and unique genetic changes underlying fast growth in mangrove-unfavourable conditions and sheds light on the molecular mechanisms of the white mangrove utility in ecological restoration.

Genome assembly and resequencing analyses provide new insights into the evolution, domestication and ornamental traits of crape myrtle

Crape myrtle (Lagerstroemia indica) is a globally used ornamental woody plant and is the representative species of Lagerstroemia. However, studies on the evolution and genomic breeding of L. indica have been hindered by the lack of a reference genome. Here we assembled the first high-quality genome of L. indica using PacBio combined with Hi-C scaffolding to anchor the 329.14-Mb genome assembly into 24 pseudochromosomes. We detected a previously undescribed independent whole-genome triplication event occurring 35.5 million years ago in L. indica following its divergence from Punica granatum. After resequencing 73 accessions of Lagerstroemia, the main parents of modern crape myrtle cultivars were found to be L. indica and L. fauriei. During the process of domestication, genetic diversity tended to decrease in many plants, but this was not observed in L. indica. We constructed a high-density genetic linkage map with an average map distance of 0.33 cM. Furthermore, we integrated the results of quantitative trait locus (QTL) using genetic mapping and bulk segregant analysis (BSA), revealing that the major-effect interval controlling internode length (IL) is located on chr1, which contains CDL15, CRG98, and GID1b1 associated with the phytohormone pathways. Analysis of gene expression of the red, purple, and white flower-colour flavonoid pathways revealed that differential expression of multiple genes determined the flower colour of L. indica, with white flowers having the lowest gene expression. In addition, BSA of purple- and green-leaved individuals of populations of L. indica was performed, and the leaf colour loci were mapped to chr12 and chr17. Within these intervals, we identified MYB35, NCED, and KAS1. Our genome assembly provided a foundation for investigating the evolution, population structure, and differentiation of Myrtaceae species and accelerating the molecular breeding of L. indica.

Chromosome-scale genomes of commercial timber trees (Ochroma pyramidale, Mesua ferrea, and Tectona grandis)

Wood is the most important natural and endlessly renewable source of energy. Despite the ecological and economic importance of wood, many aspects of its formation have not yet been investigated. We performed chromosome-scale genome assemblies of three timber trees (Ochroma pyramidale, Mesua ferrea, and Tectona grandis) which exhibit different wood properties such as wood density, hardness, growth rate, and fiber cell wall thickness. The combination of 10X, stLFR, Hi-Fi sequencing and HiC data led us to assemble high-quality genomes evident by scaffold N50 length of 55.97 Mb (O. pyramidale), 22.37 Mb (M. ferrea) and 14.55 Mb (T. grandis) with >97% BUSCO completeness of the assemblies. A total of 35774, 24027, and 44813 protein-coding genes were identified in M. ferrea, T. grandis and O. pyramidale, respectively. The data generated in this study is anticipated to serve as a valuable genetic resource and will promote comparative genomic analyses, and it is of practical importance in gaining a further understanding of the wood properties in non-model woody species.

Multi-omics reveal key enzymes involved in the formation of phenylpropanoid glucosides in Artemisia annua

Although mainly known for producing artemisinin, Artemisia annua is enriched in phenylpropanoid glucosides (PGs) with significant bioactivities. However, the biosynthesis of A. annua PGs is insufficiently investigated. Different A. annua ecotypes from distinct growing environments accumulate varying amounts of metabolites, including artemisinin and PGs such as scopolin. UDP-glucose:phenylpropanoid glucosyltransferases (UGTs) transfers glucose from UDP-glucose in PG biosynthesis. Here, we found that the low-artemisinin ecotype GS produces a higher amount of scopolin, compared to the high-artemisinin ecotype HN. By combining transcriptome and proteome analyses, we selected 28 candidate AaUGTs from 177 annotated AaUGTs. Using AlphaFold structural prediction and molecular docking, we determined the binding affinities of 16 AaUGTs. Seven of the AaUGTs enzymatically glycosylated phenylpropanoids. AaUGT25 converted scopoletin to scopolin and esculetin to esculin. The lack of accumulation of esculin in the leaf and the high catalytic efficiency of AaUGT25 on esculetin suggest that esculetin is methylated to scopoletin, the precursor of scopolin. We also discovered that AaOMT1, a previously uncharacterized O-methyltransferase, converts esculetin to scopoletin, suggesting an alternative route for producing scopoletin, which contributes to the high-level accumulation of scopolin in A. annua leaves. AaUGT1 and AaUGT25 responded to induction of stress-related phytohormones, implying the involvement of PGs in stress responses.

Integrated transcriptomic and metabolomic analysis reveals the effects of polyploidization on the lignin content and metabolic pathway in Eucalyptus

BACKGROUND

Lignin is a major restriction factor for the industrial production of biomass resources, such as pulp and bioenergy. Eucalyptus is one of the most important sources of pulp and bioenergy. After polyploidization, the lignin content of forest trees is generally reduced, which is considered a beneficial genetic improvement. However, the differences in the lignin content between triploid and diploid Eucalyptus and the underlying regulatory mechanism are still unclear.

RESULTS

We conducted a comprehensive analysis at the phenotypic, transcriptional and metabolite levels between Eucalyptus urophylla triploids and diploids to reveal the effects of polyploidization on the lignin content and lignin metabolic pathway. The results showed that the lignin content of Eucalyptus urophylla triploid stems was significantly lower than that of diploids. Lignin-related metabolites were differentially accumulated between triploids and diploids, among which coniferaldehyde, p-coumaryl alcohol, sinapaldehyde and coniferyl alcohol had significant positive correlations with lignin content, indicating that they might be primarily contributing metabolites. Most lignin biosynthetic genes were significantly downregulated, among which 11 genes were significantly positively correlated with the lignin content and above metabolites. Furthermore, we constructed a co-expression network between lignin biosynthetic genes and transcription factors based on weighted gene co-expression network analysis. The network identified some putative orthologues of secondary cell wall (SCW)-related transcription factors, among which MYB52, MYB42, NAC076, and LBD15 were significantly downregulated in Eucalyptus urophylla triploids. In addition, potential important transcription factors, including HSL1, BEE3, HHO3, and NAC046, also had high degrees of connectivity and high edge weights with lignin biosynthetic genes, indicating that they might also be involved in the variation of lignin accumulation between triploid and diploid Eucalyptus urophylla.

CONCLUSIONS

The results demonstrated that some lignin-related metabolites, lignin biosynthetic genes and transcription factors in Eucalyptus urophylla triploids may be relatively sensitive in response to the polyploidization effect, significantly changing their expression levels, which ultimately correlated with the varied lignin content. The analysis of the underlying formation mechanism could provide beneficial information for the development and utilization of polyploid biomass resources, which will be also valuable for genetic improvement in other bioenergy plants.

Knockout of floral and meiosis genes using CRISPR/Cas9 produces male-sterility in Eucalyptus without impacts on vegetative growth

Eucalyptus spp. are widely cultivated for the production of pulp, energy, essential oils, and as ornamentals. However, their dispersal from plantings, especially when grown as an exotic, can cause ecological disruptions. To provide new tools for prevention of sexual dispersal by pollen as well as to induce male-sterility for hybrid breeding, we studied the clustered regularly interspaced short palindromic repeats (CRISPR)/Cas9-mediated knockout of three floral genes in both FT-expressing (early-flowering) and non-FT genotypes. We report male-sterile phenotypes resulting from knockout of the homologs of all three genes, including one involved in meiosis and two regulating early stages of pollen development. The targeted genes were Eucalyptus homologs of REC8 (EREC8), TAPETAL DEVELOPMENT AND FUNCTION 1 (ETDF1), and HECATE3 (EHEC3-like). The erec8 knockouts yielded abnormal pollen grains and a predominance of inviable pollen, whereas the etdf1 and ehec3-like knockouts produced virtually no pollen. In addition to male-sterility, both erec8 and ehec3-like knockouts may provide complete sterility because the failure of erec8 to undergo meiosis is expected to be independent of sex, and ehec3-like knockouts produce flowers with shortened styles and no visible stigmas. When comparing knockouts to controls in wild-type (non-early-flowering) backgrounds, we did not find visible morphological or statistical differences in vegetative traits, including average single-leaf mass, stem volume, density of oil glands, or chlorophyll in leaves. Loss-of-function mutations in any of these three genes show promise as a means of inducing male- or complete sterility without impacting vegetative development.

Duplicate Genes Contribute to Variability in Abiotic Stress Resistance in Allopolyploid Wheat

Gene duplication is a universal biological phenomenon that drives genomic variation and diversity, plays a crucial role in plant evolution, and contributes to innovations in genetic engineering and crop development. Duplicated genes participate in the emergence of novel functionality, such as adaptability to new or more severe abiotic stress resistance. Future crop research will benefit from advanced, mechanistic understanding of the effects of gene duplication, especially in the development and deployment of high-performance, stress-resistant, elite wheat lines. In this review, we summarize the current knowledge of gene duplication in wheat, including the principle of gene duplication and its effects on gene function, the diversity of duplicated genes, and how they have functionally diverged. Then, we discuss how duplicated genes contribute to abiotic stress response and the mechanisms of duplication. Finally, we have a future prospects section that discusses the direction of future efforts in the short term regarding the elucidation of replication and retention mechanisms of repetitive genes related to abiotic stress response in wheat, excellent gene function research, and practical applications.

Chromosome-level assemblies of cultivated water chestnut Trapa bicornis and its wild relative Trapa incisa

Water chestnut (Trapa L.) is a floating-leaved aquatic plant with high edible and medicinal value. In this study, we presented chromosome-level genome assemblies of cultivated large-seed species Trapa bicornis and its wild small-seed relative Trapa incisa by using PacBio HiFi long reads and Hi-C technology. The T. bicornis and T. incisa assemblies consisted of 479.90 Mb and 463.97 Mb contigs with N50 values of 13.52 Mb and 13.77 Mb, respectively, and repeat contents of 62.88% and 62.49%, respectively. A total of 33,306 and 33,315 protein-coding genes were predicted in T. bicornis and T. incisa assemblies, respectively. There were 159,232 structural variants affecting more than 11 thousand genes detected between the two genomes. The phylogenetic analysis indicated that the lineage leading to Trapa was diverged from the lineage to Sonneratia approximately 23 million years ago. These two assemblies provide valuable resources for future evolutionary and functional genomic research and molecular breeding of water chestnut.

Using knowledge graphs to infer gene expression in plants

Introduction

Climate change is already affecting ecosystems around the world and forcing us to adapt to meet societal needs. The speed with which climate change is progressing necessitates a massive scaling up of the number of species with understood genotype-environment-phenotype (G×E×P) dynamics in order to increase ecosystem and agriculture resilience. An important part of predicting phenotype is understanding the complex gene regulatory networks present in organisms. Previous work has demonstrated that knowledge about one species can be applied to another using ontologically-supported knowledge bases that exploit homologous structures and homologous genes. These types of structures that can apply knowledge about one species to another have the potential to enable the massive scaling up that is needed through in silico experimentation.

Methods

We developed one such structure, a knowledge graph (KG) using information from Planteome and the EMBL-EBI Expression Atlas that connects gene expression, molecular interactions, functions, and pathways to homology-based gene annotations. Our preliminary analysis uses data from gene expression studies in Arabidopsis thaliana and Populus trichocarpa plants exposed to drought conditions.

Results

A graph query identified 16 pairs of homologous genes in these two taxa, some of which show opposite patterns of gene expression in response to drought. As expected, analysis of the upstream cis-regulatory region of these genes revealed that homologs with similar expression behavior had conserved cis-regulatory regions and potential interaction with similar trans-elements, unlike homologs that changed their expression in opposite ways.

Discussion

This suggests that even though the homologous pairs share common ancestry and functional roles, predicting expression and phenotype through homology inference needs careful consideration of integrating cis and trans-regulatory components in the curated and inferred knowledge graph.

Transcriptional reprogramming during recovery from drought stress in Eucalyptus grandis

The importance of drought as a constraint to agriculture and forestry is increasing with climate change. Genetic improvement of plants' resilience is one of the mitigation strategies to curb this threat. Although recovery from drought stress is important to long-term drought adaptation and has been considered as an indicator of dehydration tolerance in annual crops, this has not been well explored in forest trees. Thus, we aimed to investigate the physiological and transcriptional changes during drought stress and rewatering in Eucalyptus grandis W. Hill ex Maiden. We set up a greenhouse experiment where we imposed drought stress on 2-year-old seedlings and rewatered the recovery group after 17 days of drought. Our measurement of leaf stomatal conductance (gs) showed that, while gs was reduced by drought stress, it fully recovered after 5 days of rewatering. The RNA-seq analysis from stem samples revealed that genes related to known stress responses such as phytohormone and reactive oxygen species signaling were upregulated, while genes involved in metabolism and growth were downregulated due to drought stress. We observed reprogramming of signal transduction pathways and metabolic processes at 1 day of rewatering, indicating a quick response to rewatering. Our results suggest that recovery from drought stress may entail alterations in the jasmonic acid, salicylic acid, ethylene and brassinosteroid signaling pathways. Using co-expression network analysis, we identified hub genes, including the putative orthologs of ABI1, ABF2, ABF3, HAI2, BAM1, GolS2 and SIP1 during drought and CAT2, G6PD1, ADG1 and FD-1 during recovery. Taken together, by highlighting the molecular processes and identifying key genes, this study gives an overview of the mechanisms underlying the response of E. grandis to drought stress and recovery that trees may face repeatedly throughout their long life cycle. This provides a useful reference to the identification and further investigation of signaling pathways and target genes for future tree improvement.

The In Planta Gene Expression of Austropuccinia psidii in Resistant and Susceptible Eucalyptus grandis

Austropuccinia psidii, commonly known as myrtle rust, is an obligate, biotrophic rust pathogen that causes rust disease in a broad host range of Myrtaceae species. Eucalyptus grandis, a widely cultivated hardwood Myrtaceae species, is susceptible to A. psidii infection, with this pathogen threatening both their natural range and various forest plantations across the world. This study aimed to investigate the A. psidii transcriptomic responses in resistant and susceptible E. grandis at four time points. RNA-seq reads were mapped to the A. psidii reference genome to quantify expressed genes at 12 h postinoculation and 1, 2, and 5 days postinoculation (dpi). A total of eight hundred and ninety expressed genes were found, of which 43 were candidate effector protein genes. These included rust transferred protein 1 (RTP1), expressed in susceptible hosts at 5 dpi, and a hydrolase protein gene expressed in both resistant and susceptible hosts over time. Functional categorization of expressed genes revealed processes enriched in susceptible hosts, including malate metabolic and malate dehydrogenase activity, implicating oxalic acid in disease susceptibility. These results highlight putative virulence or pathogenicity mechanisms employed by A. psidii to cause disease, and they provide the first insight into the molecular responses of A. psidii in E. grandis over time.

E. urophylla × E. grandis high-quality genome and comparative genomics provide insights on evolution and diversification of eucalyptus

BACKGROUND

Eucalyptus urophylla × Eucalyptus grandis, an economically important forest tree, provides important raw material for energy and reduces damage to native forests. However, the absence of a high-quality E. urophylla × E. grandis reference genome has significantly hindered its evolution and genetic analysis.

RESULTS

We successfully presented a high-quality reference genome of E. urophylla × E. grandis (545.75 Mb; scaffold N50, 51.62 Mb) using a combination of the Illumina, PacBio HiFi, and Hi-C sequencing platforms. A total of 34,502 genes and 58.56% of the repetitive sequences in this genome were annotated. Using genome evolution analyses, we identified a recent whole-genome duplication (WGD) event in E. urophylla × E. grandis. We further found that gene families associated with starch and sucrose metabolism, flavonoid biosynthesis, and plant-pathogen interaction were significantly expanded in E. urophylla × E. grandis. Moreover, comparative genomic and evolutionary analyses showed large structural variations among the different chromosomes of the 34 Eucalyptus accessions, which were divided into six clades.

CONCLUSIONS

Overall, our findings provide a valuable resource for expanding our understanding of the E. urophylla × E. grandis genome evolution, genetic improvement, and its comparative biology.

Transcriptome-Based Construction of the Gibberellin Metabolism and Signaling Pathways in Eucalyptus grandis × E. urophylla, and Functional Characterization of GA20ox and GA2ox in Regulating Plant Development and Abiotic Stress Adaptations

Gibberellins (GAs) are the key regulators controlling plant growth, wood production and the stress responses in perennial woody plants. The role of GA in regulating the above-mentioned processes in Eucalyptus remain largely unclear. There is still a lack of systematic identification and functional characterization of GA-related genes in Eucalyptus. In this study, a total of 59,948 expressed genes were identified from the major vegetative tissues of the E. grandis × E. urophylla using transcriptome sequencing. Then, the key gene families in each step of GA biosynthesis, degradation and signaling were investigated and compared with those of Arabidopsis, rice, and Populus. The expression profile generated using Real-time quantitative PCR showed that most of these genes exhibited diverse expression patterns in different vegetative organs and in response to abiotic stresses. Furthermore, we selectively overexpressed EguGA20ox1, EguGA20ox2 and EguGA2ox1 in both Arabidopsis and Eucalyptus via Agrobacterium tumefaciens or A. rhizogenes-mediated transformation. Though both Arabidopsis EguGA20ox1- and EguGA20ox2-overexpressing (OE) lines exhibited better vegetative growth performance, they were more sensitive to abiotic stress, unlike EguGA2ox1-OE plants, which exhibited enhanced stress resistance. Moreover, overexpression of EguGA20ox in Eucalyptus roots caused significantly accelerated hairy root initiation and elongation and improved root xylem differentiation. Our study provided a comprehensive and systematic study of the genes of the GA metabolism and signaling and identified the role of GA20ox and GA2ox in regulating plant growth, stress tolerance, and xylem development in Eucalyptus; this could benefit molecular breeding for obtaining high-yield and stress-resistant Eucalyptus cultivars.

Insights on the Impact of Arbuscular Mycorrhizal Symbiosis on Eucalyptus grandis Tolerance to Drought Stress

Drought stress has a negative impact on plant growth and production. Arbuscular mycorrhizal (AM) fungi, which establish symbioses with most terrestrial vascular plant species, play important roles in improving host plant mineral nutrient acquisition and resistance to drought. However, the physiological and molecular regulation mechanisms occurring in mycorrhizal Eucalyptus grandis coping with drought stress remain unclear. Here, we studied the physiological changes and mitogen-activated protein kinase (MAPK) cascade gene expression profiles of E. grandis associated with AM fungi under drought stress. The results showed that colonization by AM fungi significantly enhanced plant growth, with higher plant biomass, shoot height, root length, and relative water content (RWC) under drought conditions. Mycorrhizal plants had lower levels of accumulation of proline, malondialdehyde (MDA), H2O2, and O2·- than seedlings not colonized with AM fungi. In addition, mycorrhizal E. grandis also had higher peroxidase (POD), superoxide dismutase (SOD), and catalase (CAT) activities under drought conditions, improving the antioxidant system response. Eighteen MAPK cascade genes were isolated from E. grandis, and the expression levels of the MAPK cascade genes were positively induced by symbiosis with AM fungi, which was correlated with changes in the proline, MDA, H2O2, and O2·- contents and POD, SOD, and CAT activities. In summary, our results showed that AM symbiosis enhances E. grandis drought tolerance by regulating plant antioxidation abilities and MAPK cascade gene expression. IMPORTANCE Arbuscular mycorrhizal (AM) fungi play an important role in improving plant growth and development under drought stress. The MAPK cascade may regulate many physiological and biochemical processes in plants in response to drought stress. Previous studies have shown that there is a complex regulatory network between the plant MAPK cascade and drought stress. However, the relationship between the E. grandis MAPK cascade and AM symbiosis in coping with drought remains to be investigated. Our results suggest that AM fungi could improve plant drought tolerance mainly by improving the antioxidant ability to protect plants from reactive oxygen species (ROS) and alleviate oxidative stress damage. The expression of the MAPK cascade genes was induced in mycorrhizal E. grandis seedlings under drought stress. This study revealed that MAPK cascade regulation is of special significance for improving the drought tolerance of E. grandis. This study provides a reference for improving mycorrhizal seedling cultivation under stress.

Genome-wide identification, expression profile and evolutionary relationships of TPS genes in the neotropical fruit tree species Psidium cattleyanum

Terpenoids are essential for plant growth, development, defense, and adaptation mechanisms. Psidium cattleyanum (Myrtaceae) is a fleshy fruit tree species endemics from Atlantic Forest, known for its pleasant fragrance and sweet taste, attributed to terpenoids in its leaves and fruits. In this study, we conducted genome-wide identification, evolutionary and expression analyses of the terpene synthase gene (TPS) family in P. cattleyanum red guava (var. cattleyanum), and yellow guava (var. lucidum Hort.) morphotypes. We identified 32 full-length TPS in red guava (RedTPS) and 30 in yellow guava (YlwTPS). We showed different expression patterns of TPS paralogous in the two morphotypes, suggesting the existence of distinct gene regulation mechanisms and their influence on the final essential oil content in both morphotypes. Moreover, the oil profile of red guava was dominated by 1,8-cineole and linalool and yellow guava was enriched in α-pinene, coincident in proportion to TPS-b1 genes, which encode enzymes that produce cyclic monoterpenes, suggesting a lineage-specific subfamily expansion of this family. Finally, we identified amino acid residues near the catalytic center and functional areas under positive selection. Our findings provide valuable insights into the terpene biosynthesis in a Neotropical Myrtaceae species and their potential involvement in adaptation mechanisms.