The genome of Eucalyptus grandis

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.

Gap-free genome assembly and comparative analysis reveal the evolution and anthocyanin accumulation mechanism of Rhodomyrtus tomentosa

Rhodomyrtus tomentosa is an important fleshy-fruited tree and a well-known medicinal plant of the Myrtaceae family that is widely cultivated in tropical and subtropical areas of the world. However, studies on the evolution and genomic breeding of R. tomentosa were hindered by the lack of a reference genome. Here, we presented a chromosome-level gap-free T2T genome assembly of R. tomentosa using PacBio and ONT long read sequencing. We assembled the genome with size of 470.35 Mb and contig N50 of ~43.80 Mb with 11 pseudochromosomes. A total of 33 382 genes and 239.31 Mb of repetitive sequences were annotated in this genome. Phylogenetic analysis elucidated the independent evolution of R. tomentosa starting from 14.37MYA and shared a recent WGD event with other Myrtaceae species. We identified four major compounds of anthocyanins and their synthetic pathways in R. tomentosa. Comparative genomic and gene expression analysis suggested the coloring and high anthocyanin accumulation in R. tomentosa tends to be determined by the activation of anthocyanin synthesis pathway. The positive selection and up-regulation of MYB transcription factors were the implicit factors in this process. The copy number increase of downstream anthocyanin transport-related OMT and GST gene were also detected in R. tomentosa. Expression analysis and pathway identification enriched the importance of starch degradation, response to stimuli, effect of hormones, and cell wall metabolism during the fleshy fruit development in Myrtaceae. Our genome assembly provided a foundation for investigating the origins and differentiation of Myrtaceae species and accelerated the genetic improvement of R. tomentosa.

Interspecies genome divergence is predominantly due to frequent small scale rearrangements in Eucalyptus

Synteny, the ordering of sequences within homologous chromosomes, must be maintained within the genomes of sexually reproducing species for the sharing of alleles and production of viable, reproducing offspring. However, when the genomes of closely related species are compared, a loss of synteny is often observed. Unequal homologous recombination is the primary mechanism behind synteny loss, occurring more often in transposon rich regions, and resulting in the formation of chromosomal rearrangements. To examine patterns of synteny among three closely related, interbreeding, and wild Eucalyptus species, we assembled their genomes using long-read DNA sequencing and de novo assembly. We identify syntenic and rearranged regions between these genomes and estimate that ~48% of our genomes remain syntenic while ~36% is rearranged. We observed that rearrangements highly fragment microsynteny. Our results suggest that synteny between these species is primarily lost through small-scale rearrangements, not through sequence loss, gain, or sequence divergence. Further examination of identified rearrangements suggests that rearrangements may be altering the phenotypes of Eucalyptus species. Our study also underscores that the use of single reference genomes in genomic variation studies could lead to reference bias, especially given the scale at which we show potentially adaptive loci have highly diverged, deleted, duplicated and/or rearranged. This study provides an unbiased framework to look at potential speciation and adaptive loci among a rapidly radiating foundation species of woodland trees that are free from selective breeding seen in most crop species.

Genome-wide analysis of R2R3-MYB transcription factors reveals their differential responses to drought stress and ABA treatment in desert poplar (Populus euphratica)

The R2R3-MYB transcription factors are widely involved in the regulation of plant growth, biotic and abiotic stress responses. Meanwhile, seed germination, which is stimulated by internal and external environments, is a critical stage in the plant life cycle. However, the identification, characterization, and expression profiling of the Populus euphratica R2R3-MYB family in drought response during seed germination have been unknown. Our study attempted to identify and characterize the R2R3-MYB genes in P. euphratica (PeR2R3-MYBs) and explore how R2R3-MYBs trigger the drought and abscisic acid (ABA) response mechanism in its seedlings. Based on the analysis of comparative genomics, 174 PeR2R3-MYBs were identified and expanded driven by whole genome duplication or segment duplication events. The analysis of Ka/Ks ratios showed that, in contrast to most PeR2R3-MYBs, the other PeR2R3-MYBs were subjected to positive selection in P. euphratica. Further, the expression data of PeR2R3-MYBs under drought stress and ABA treatment, together with available functional data for Arabidopsis thaliana MYB genes, supported the hypothesis that PeR2R3-MYBs involved in response to drought are dependent or independent on ABA signaling pathway during seed germination, especially PeR2R3-MYBs with MYB binding sites (MBS) cis-element and/or tandem duplication. This study is the first report on the genome-wide analysis of PeR2R3-MYBs, as well as the other two Salicaceae species. The duplication events and differential expressions of PeR2R3-MYBs play important roles in enhancing the adaptation to drought desert environment. Our results provide a reference for prospective functional studies of R2R3-MYBs of poplars and lay the foundation for new breeding strategies to improve the drought tolerance of P. euphratica.

Comparative phylogenomic and structural analysis of canonical secretory PLA2 and novel PLA2-like family in plants

Plant secretory phospholipase A₂ (sPLA₂) is a family of lipolytic enzymes involved in the sn-2 hydrolysis of phospholipid carboxyester bonds, characterized by the presence of a conserved PA2c domain. PLA₂ produces free fatty acids and lysophospholipids, which regulate several physiological functions, including lipid metabolism, plant growth and development, signal transduction, and response to various environmental stresses. In the present work, we have performed a comparative analysis of PA2c domain-containing genes across plants, focusing on gene distribution, phylogenetic analysis, tissue-specific expression, and homology modeling. Our data revealed the widespread occurrence of multiple sPLA₂ in most land plants and documented single sPLA₂ in multiple algal groups, indicating an ancestral origin of sPLA₂. We described a novel PA2c-containing gene family present in all plant lineages and lacking secretory peptide, which we termed PLA₂-like. Phylogenetic analysis revealed two independent clades in canonical sPLA₂ genes referred to as α and β clades, whereas PLA₂-like genes clustered independently as a third clade. Further, we have explored clade-specific gene expressions showing that while all three clades were expressed in vegetative and reproductive tissues, only sPLA₂-β and PLA₂-like members were expressed in the pollen and pollen tube. To get insight into the conservation of the gene regulatory network of sPLA₂ and PLA₂-like genes, we have analyzed the occurrence of various cis-acting promoter elements across the plant kingdom. The comparative 3D structure analysis revealed conserved and unique features within the PA2c domain for the three clades. Overall, this study will help to understand the evolutionary significance of the PA2c family and lay the foundation for future sPLA₂ and PLA₂-like characterization in plants.

Recruitment of distinct UDP-glycosyltransferase families demonstrates dynamic evolution of chemical defense within Eucalyptus L'Hér

The economic and ecologically important genus Eucalyptus is rich in structurally diverse specialized metabolites. While some specialized metabolite classes are highly prevalent across the genus, the cyanogenic glucoside prunasin is only produced by c. 3% of species. To investigate the evolutionary mechanisms behind prunasin biosynthesis in Eucalyptus, we compared de novo assembled transcriptomes, together with online resources between cyanogenic and acyanogenic species. Identified genes were characterized in vivo and in vitro. Pathway characterization of cyanogenic Eucalyptus camphora and Eucalyptus yarraensis showed for the first time that the final glucosylation step from mandelonitrile to prunasin is catalyzed by a novel UDP-glucosyltransferase UGT87. This step is typically catalyzed by a member of the UGT85 family, including in Eucalyptus cladocalyx. The upstream conversion of phenylalanine to mandelonitrile is catalyzed by three cytochrome P450 (CYP) enzymes from the CYP79, CYP706, and CYP71 families, as previously shown. Analysis of acyanogenic Eucalyptus species revealed the loss of different ortholog prunasin biosynthetic genes. The recruitment of UGTs from different families for prunasin biosynthesis in Eucalyptus demonstrates important pathway heterogeneities and unprecedented dynamic pathway evolution of chemical defense within a single genus. Overall, this study provides relevant insights into the tremendous adaptability of these long-lived trees.

Draft genome of the medicinal tea tree Melaleuca alternifolia

BACKGROUND

Melaleuca alternifolia is a commercially important medicinal tea tree native to Australia. Tea tree oil, the essential oil distilled from its branches and leaves, has broad-spectrum germicidal activity and is highly valued in the pharmaceutical and cosmetic industries. Thus, the study of genome, which can provide reference for the investigation of genes involved in terpinen-4-ol biosynthesis, is quite crucial for improving the productivity of Tea tree oil.

METHODS AND RESULTS

In our study, the next-generation sequencing was used to investigate the whole genome of Melaleuca alternifolia. About 114 Gb high quality sequence data were obtained and assembled into 1,838,159 scafolds with an N50 length of 1021 bp. The assembled genome size is about 595 Mb, twice of that predicted by flow cytometer (300 Mb) and k-mer analysis (345 Mb). Benchmarking Universal Single-Copy Orthologs analyses indicated that only 11.3% of the conserved single-copy genes were miss. Repetitive regions cover over 40.43% of the genome. A total of 44,369 protein-coding genes were predicted and annotated against Nr, Swissprot, Refseq, COG, KOG, and KEGG database. Among these genes, 32,909 and 16,241 genes were functionally annotated in Nr and KEGG, respectively. Moreover, 29,411 and 14,435 genes were functionally annotated in COG and KOG. Additionally, 457,661 simple sequence repeats and 1109 transcription factors (TFs) form 67 TF families were identified in the assembled genome.

CONCLUSION

Our findings provide a draft genome sequencing of M. alternifolia which can act as a reference for the deep sequencing strategies, and are useful for future functional and comparative genomics analyses.

Chromosomal-level genome assembly of Melastoma candidum provides insights into trichome evolution

Melastoma, consisting of ~100 species diversified in tropical Asia and Oceania in the past 1-2 million years, represents an excellent example of rapid speciation in flowering plants. Trichomes on hypanthia, twigs and leaves vary markedly among species of this genus and are the most important diagnostic traits for species identification. These traits also play critical roles in contributing to differential adaptation of these species to their own habitats. Here we sequenced the genome of M. candidum, a common, erect-growing species from southern China, with the aim to provide genomic insights into trichome evolution in this genus. We generated a high-quality, chromosome-level genome assembly of M. candidum, with the genome size of 256.2 Mb and protein-coding gene number of 40,938. The gene families specific to, and significantly expanded in Melastoma are enriched for GO terms related to trichome initiation and differentiation. We provide evidence that Melastoma and its sister genus Osbeckia have undergone two whole genome duplications (WGDs) after the triplication event (γ) shared by all core eudicots. Preferential retention of trichome development-related transcription factor genes such as C2H2, bHLH, HD-ZIP, WRKY, and MYB after both WGDs might provide raw materials for trichome evolution and thus contribute to rapid species diversification in Melastoma. Our study provides candidate transcription factor genes related to trichome evolution in Melastoma, which can be used to evolutionary and functional studies of trichome diversification among species of this genus.

Genomic and epigenomic variation in Psidium species and their outcome under the yield and composition of essential oils

Diploid and polyploid species derived from the euploid series x = 11 occur in the genus Psidium, as well as intraspecific cytotypes. Euploidy in the genus can alter the gene copy number, resulting in several "omics" variations. We revisited the euploidy, reported genomic (nuclear 2C value, GC%, and copy number of secondary metabolism genes) and epigenomic (5-mC%) differences in Psidium, and related them to essential oil yield and composition. Mean 2C values ranged from 0.90 pg (P. guajava) to 7.40 pg (P. gaudichaudianum). 2C value is intraspecifically varied in P. cattleyanum and P. gaudichaudianum, evidencing cytotypes that can be formed from euploid (non-reduced) and/or aneuploid reproductive cells. GC% ranged from 34.33% (P. guineense) to 48.95% (P. myrtoides), and intraspecific variations occurred even for species without 2C value intraspecific variation. Essential oil yield increased in relation to 2C value and to GC%. We showed that P. guajava (diploid) possesses two and P. guineense (tetraploid) four copies of the one specific TPS gene, as well as eight and sixteen copies respectively of the conserved regions that occur in eight TPS genes. We provide a wide "omics'' characterization of Psidium and show the outcome of the genome and epigenome variation in secondary metabolism.

Flavan-3-ol Galloylation-Related Functional Gene Cluster and the Functional Diversification of SCPL Paralogs in Camellia sp

The high accumulation of galloylated flavan-3-ols in Camellia sp. is a noteworthy phenomenon. We identified a flavan-3-ol galloylation-related functional gene cluster in tannin-rich plant Camellia sp., which included UGT84A22 and SCPL-AT gene clusters. We investigated the possible correlation between the accumulation of metabolites and the expression of SCPL-ATs and UGT84A22. The results revealed that C. sinensis, C. ptilophylla, and C. oleifera accumulated galloylated cis-flavan-3-ols (EGCG), galloylated trans-flavan-3-ols (GCG), and hydrolyzed tannins, respectively; however, C. nitidissima did not accumulate any galloylated compounds. C. nitidissima exhibited no expression of SCPL-AT or UGT84A22, whereas the other three species of Camellia exhibited various expression patterns. This indicated that the functions of the paralogs of SCPL-AT vary. Enzymatic analysis revealed that SCPL5 was neofunctionalized as a noncatalytic chaperone paralog, a type of chaerone-like protein, associating with flavan-3-ol galloylation; moreover, CsSCPL4 was subfunctionalized in association with the galloylation of cis- and trans-flavan-3-ols. In C. nitidissima, an SCPL4 homolog was noted with mutations in two cysteine residues forming a disulfide bond, which suggested that this homolog was defunctionalized. The findings of this study improve our understanding of the functional diversification of SCPL paralogs in Camellia sp.

Transcriptome Profiling of Gossypium anomalum Seedlings Reveals Key Regulators and Metabolic Pathways in Response to Drought Stress

Drought stress is a key limiting factor for cotton (Gossypium spp.) growth, production, development, and production worldwide. Some wild diploid cotton species are remarkably tolerant of water deficit and constitute an important reservoir for understanding the molecular mechanisms of Gossypium spp. drought tolerance and improving cultivated upland cotton. Here, we utilized RNA-Seq technology to characterize the leaf transcriptomes of a wild African diploid cotton species, Gossypium anomalum, under drought stress. A total of 12,322 differentially expressed genes (DEGs) were identified after mapping valid clean reads to the reference genome of G. anomalum, of which 1243 were commonly differentially expressed at all stages of drought stress. These genes were significantly enriched for molecular functions Gene Ontology terms related to cytoskeleton, hydrolase activity, cellular redox, and binding. Additionally, a substantial proportion of enriched biological process terms concerned cell or subcellular processes, while most in the cellular components category concerned membrane function and photosynthesis. An enrichment analysis against the Kyoto Encyclopedia of Genes and Genomes showed the top significantly enriched pathways to be photosynthesis-antenna proteins, amino sugar and nucleotide sugar metabolism, starch and sucrose metabolism, MAPK signaling pathway, glutathione metabolism, and plant hormone signal transduction. The DEGs also exhibited interestingly significant enrichments for drought stress-induced tandemly repeated genes involved in iron ion binding, oxidoreductase activity, heme binding, and other biological processes. A large number of genes encoding transcription factors, such as MYB, bHLH, ERF, NAC, WRKY, and bZIP, were identified as playing key roles in acclimatizing to drought stress. These results will provide deeper insights into the molecular mechanisms of drought stress adaptation in Gossypium spp.

Tomato POLLEN DEFICIENT 2 encodes a G-type lectin receptor kinase required for viable pollen grain formation

Pollen development is a crucial biological process indispensable for seed set in flowering plants and for successful crop breeding. However, little is known about the molecular mechanisms regulating pollen development in crop species. This study reports a novel male-sterile tomato mutant, pollen deficient 2 (pod2), characterized by the production of non-viable pollen grains and resulting in the development of small parthenocarpic fruits. A combined strategy of mapping-by-sequencing and RNA interference-mediated gene silencing was used to prove that the pod2 phenotype is caused by the loss of Solanum lycopersicum G-type lectin receptor kinase II.9 (SlG-LecRK-II.9) activity. In situ hybridization of floral buds showed that POD2/SlG-LecRK-II.9 is specifically expressed in tapetal cells and microspores at the late tetrad stage. Accordingly, abnormalities in meiosis and tapetum programmed cell death in pod2 occurred during microsporogenesis, resulting in the formation of four dysfunctional microspores leading to an aberrant microgametogenesis process. RNA-seq analyses supported the existence of alterations at the final stage of microsporogenesis, since we found tomato deregulated genes whose counterparts in Arabidopsis are essential for the normal progression of male meiosis and cytokinesis. Collectively, our results revealed the essential role of POD2/SlG-LecRK-II.9 in regulating tomato pollen development.

Proteomics research in forest trees: A 2012-2022 update

This review is a compilation of proteomic studies on forest tree species published in the last decade (2012-2022), mostly focused on the most investigated species, including Eucalyptus, Pinus, and Quercus. Improvements in equipment, platforms, and methods in addition to the increasing availability of genomic data have favored the biological knowledge of these species at the molecular, organismal, and community levels. Integration of proteomics with physiological, biochemical and other large-scale omics in the direction of the Systems Biology, will provide a comprehensive understanding of different biological processes, from growth and development to responses to biotic and abiotic stresses. As main issue we envisage that proteomics in long-living plants will thrive light on the plant responses and resilience to global climate change, contributing to climate mitigation strategies and molecular breeding programs. Proteomics not only will provide a molecular knowledge of the mechanisms of resilience to either biotic or abiotic stresses, but also will allow the identification on key gene products and its interaction. Proteomics research has also a translational character being applied to the characterization of the variability and biodiversity, as well as to wood and non-wood derived products, traceability, allergen and bioactive peptides identification, among others. Even thought, the full potential of proteomics is far from being fully exploited in forest tree research, with PTMs and interactomics being reserved to plant model systems. The most outstanding achievements in forest tree proteomics in the last decade as well as prospects are discussed.

Haplotype mining panel for genetic dissection and breeding in Eucalyptus

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

Compilation of reaction kinetics parameters determined in the Key Development Project for Air Pollution Formation Mechanism and Control Technologies in China

A compilation of new advances made in the research field of laboratory reaction kinetics in China's Key Development Project for Air Pollution Formation Mechanism and Control Technologies was presented. These advances are grouped into six broad, interrelated categories, including volatile organic compound (VOC) oxidation, secondary organic aerosol (SOA) formation, new particle formation (NPF) and gas-particle partitioning, ozone chemistry, model parameters, and secondary inorganic aerosol (SIA) formation, highlighting the laboratory work done by Chinese researchers. For smog chamber applications, the current knowledge gained from laboratory studies is reviewed, with emphasis on summarizing the oxidation mechanisms of long-chain alkanes, aromatics, alkenes, aldehydes/ketones in the atmosphere, SOA formation from anthropogenic emission sources, and oxidation of aromatics, isoprene, and limonene, as well as SIA formation. For flow tube applications, atmospheric oxidation mechanisms of toluene and methacrolein, SOA formation from limonene oxidation by ozone, gas-particle partitioning of peroxides, and sulfuric acid-water (H₂SO₄-H₂O) binary nucleation, methanesulfonic acid-water (MSA-H₂O) binary nucleation, and sulfuric acid-ammonia-water (H₂SO₄-NH₃-H₂O) ternary nucleation are discussed.

The genome of the king protea, Protea cynaroides

The king protea (Protea cynaroides), an early-diverging eudicot, is the most iconic species from the Megadiverse Cape Floristic Region, and the national flower of South Africa. Perhaps best known for its iconic flower head, Protea is a key genus for the South African horticulture industry and cut-flower market. Ecologically, the genus and the family Proteaceae are important models for radiation and adaptation, particularly to soils with limited phosphorus bio-availability. Here, we present a high-quality chromosome-scale assembly of the P. cynaroides genome as the first representative of the fynbos biome. We reveal an ancestral whole-genome duplication event that occurred in the Proteaceae around the late Cretaceous that preceded the divergence of all crown groups within the family and its extant diversity in all Southern continents. The relatively stable genome structure of P. cynaroides is invaluable for comparative studies and for unveiling paleopolyploidy in other groups, such as the distantly related sister group Ranunculales. Comparative genomics in sequenced genomes of the Proteales shows loss of key arbuscular mycorrhizal symbiosis genes likely ancestral to the family, and possibly the order. The P. cynaroides genome empowers new research in plant diversification, horticulture and adaptation, particularly to nutrient-poor soils.

Genome-wide identification of GRF gene family and their contribution to abiotic stress response in pitaya (Hylocereus polyrhizus)

Growth-regulating factors (GRFs) are plant-specific transcription factors identified in many land plants. Recently, their indispensable roles in stress response are highlighted. In present work, 11 HpGRFs were cloned in pitaya. Segmental duplication is considered essential for the expansion of HpGRFs. A phylogenetic tree suggested that GRFs could be divided into eight categories, among which G-I was a Caryophyllales-specific one. The categorization was further evidenced by differences in the gene structure, collinearity, protein domain of HpGRFs. Five miR396 hairpins giving rise to two types of matured miR396s were identified in pitaya via sRNA-Seq in combination with bioinformatic analysis. Parallel analysis of RNA ends proved that HpGRFs except HpGRF5 were degraded by miR396-directed cleavages at the regions which code the conserved WRC motifs of HpGRFs. Multiple cis-regulatory elements were discovered in the promoters of HpGRFs. Among the elements, most are involved in stress and phytohormone response as well as plant growth, indicating a crosstalk between them. Expression analysis showed the responsive patterns of the miR396-GRF module under abiotic stresses. To conclude, our work systematically identified the miR396-targeted HpGRFs in pitaya and confirmed their involvement in stress response, providing novel insights into the comprehensive understanding of the stress resistance of pitaya.

A high-quality pseudo-phased genome for Melaleuca quinquenervia shows allelic diversity of NLR-type resistance genes

BACKGROUND

Melaleuca quinquenervia (broad-leaved paperbark) is a coastal wetland tree species that serves as a foundation species in eastern Australia, Indonesia, Papua New Guinea, and New Caledonia. While extensively cultivated for its ornamental value, it has also become invasive in regions like Florida, USA. Long-lived trees face diverse pest and pathogen pressures, and plant stress responses rely on immune receptors encoded by the nucleotide-binding leucine-rich repeat (NLR) gene family. However, the comprehensive annotation of NLR encoding genes has been challenging due to their clustering arrangement on chromosomes and highly repetitive domain structure; expansion of the NLR gene family is driven largely by tandem duplication. Additionally, the allelic diversity of the NLR gene family remains largely unexplored in outcrossing tree species, as many genomes are presented in their haploid, collapsed state.

RESULTS

We assembled a chromosome-level pseudo-phased genome for M. quinquenervia and described the allelic diversity of plant NLRs using the novel FindPlantNLRs pipeline. Analysis reveals variation in the number of NLR genes on each haplotype, distinct clustering patterns, and differences in the types and numbers of novel integrated domains.

CONCLUSIONS

The high-quality M. quinquenervia genome assembly establishes a new framework for functional and evolutionary studies of this significant tree species. Our findings suggest that maintaining allelic diversity within the NLR gene family is crucial for enabling responses to environmental stress, particularly in long-lived plants.

Haplogenome assembly reveals structural variation in Eucalyptus interspecific hybrids

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.

Genome-wide identification of sucrose non-fermenting-1-related protein kinase genes in maize and their responses to abiotic stresses

Introduction

Protein kinases play an important role in plants in response to environmental changes through signal transduction. As a large family of protein kinases, sucrose non-fermenting-1 (SNF1)-related kinases (SnRKs) were found and functionally verified in many plants. Nevertheless, little is known about the SnRK family of Zea mays.

Methods

Evolutionary relationships, chromosome locations, gene structures, conserved motifs, and cis-elements in promoter regions were systematically analyzed. Besides, tissue-specific and stress-induced expression patterns of ZmSnRKs were determined. Finally, functional regulatory networks between ZmSnRKs and other proteins or miRNAs were constructed.

Results and Discussion

In total, 60 SnRK genes located on 10 chromosomes were discovered in maize. ZmSnRKs were classified into three subfamilies (ZmSnRK1, ZmSnRK2, and ZmSnRK3), consisting of 4, 14, and 42 genes, respectively. Gene structure analysis showed that 33 of the 42 ZmSnRK3 genes contained only one exon. Most ZmSnRK genes contained at least one ABRE, MBS, and LTR cis-element and a few ZmSnRK genes had AuxRR-core, P-box, MBSI, and SARE ciselements in their promoter regions. The Ka:Ks ratio of 22 paralogous ZmSnRK gene pairs revealed that the ZmSnRK gene family had experienced a purifying selection. Meanwhile, we analyzed the expression profiles of ZmSnRKs, and they exhibited significant differences in various tissues and abiotic stresses. In addition, A total of eight ZmPP2Cs, which can interact with ZmSnRK proteins, and 46 miRNAs, which can target 24 ZmSnRKs, were identified. Generally, these results provide valuable information for further function verification of ZmSnRKs, and improve our understanding of the role of ZmSnRKs in the climate resilience of maize.

A Eucalyptus Pht1 Family Gene EgPT8 Is Essential for Arbuscule Elongation of Rhizophagus irregularis

The majority of vascular flowering plants can establish arbuscular mycorrhizal (AM) symbiosis with AM fungi. These associations contribute to plant health and plant growth against various environmental stresses. In the mutualistic endosymbiosis, the AM fungi deliver phosphate (Pi) to the host root through highly branched hyphae called arbuscules. The molecular mechanisms of Pi transfer from AM fungi to the plant have been determined, which are dominated by AM-specific Pi transporters belonging to the PHOSPHATE TRANSPORTER 1 (Pht1) family within the subfamily I. However, it is unknown whether Pht1 family proteins are involved in other regulations in AM symbiosis. Here, we report that the expression of EgPT8 is specifically activated by AM fungus Rhizophagus irregularis and is localized in root cortical cells containing arbuscules. Interestingly, knockdown of EgPT8 function does not affect the Eucalyptus grandis growth, total phosphorous (P) concentration, and arbuscule formation; however, the size of mature arbuscules was significantly suppressed in the RNAi-EgPT8 lines. Heterogeneous expression of EgPT4, EgPT5, and EgPT8 in the Medicago truncatula mutant mtpt4-2 indicates that EgPT4 and EgPT5 can fully complement the defects of mutant mtpt4-2 in mycorrhizal Pi uptake and arbuscule formation, while EgPT8 cannot complement the defective AM phenotype of the mtpt4-2 mutant. Based on our results, we propose that the AM fungi-specific subfamily I transporter EgPT8 has novel functions and is essential to arbuscule elongation. IMPORTANCE Arbuscular mycorrhizal (AM) formation in host root cortical cells is initiated by exchanges of diffusible molecules, among which Pi uptake is known as the important feature of AM fungi on symbiosis functioning. Over the last two decades, it has been repeatedly proven that most vascular plants harbor two or more AM-specific Pht1 proteins; however, there is no direct evidence regarding the potential link among these Pi transporters at the symbiotic interface. This work revealed a novel function of a structurally conserved protein involved in lateral arbuscule development. In total, we confirmed that three AM-specific Pht1 family proteins are nonredundant in Eucalyptus grandis and that EgPT8 is responsible for fungal arbuscule elongation of Rhizophagus irregularis.

A box on the river: The phylogenetics and phylogeography of Eucalyptus baueriana (Eucalyptus sect. Adnataria ser. Heterophloiae)

We present a phylogeographic study of the tree species Eucalyptus baueriana Schauer, which occurs in disjunct areas on the near coastal plains and ranges of the south-east Australian mainland. DArTseq data are used to build a phylogeny including E. baueriana and closely related taxa to test its monophyly, test the genetic distinctness of the three subspecies of E. baueriana, and investigate relationships between its disjunct populations. Additionally, we use population structure analysis to investigate the genetic distinctness of populations, and MaxEnt to investigate the environmental factors potentially influencing the species' distribution. We show E. baueriana is monophyletic and most closely related to three other Blue Box eucalypt species: E. conica H.Deane & Maiden, E. dalveenica T.L.Collins, R.L.Andrew & J.J.Bruhl and E. magnificata L.A.S.Johnson & K.D.Hill, with some evidence for genetic introgression between these taxa. Within E. baueriana, the deepest genetic breaks do not correspond with the subspecies classification as the two geographically restricted subspecies, together with samples of the more widespread E. baueriana subsp. baueriana from west of the Gippsland lowlands, form a south-western clade with that is sister to other populations of subsp. baueriana. The oldest genetic break in the species occurs in far eastern Gippsland (Victoria), corresponding to one of the shortest geographic disjunctions in the species' distribution. Genetic breaks in other species have been observed in this region which is broadly referred to as the southern transition zone. Both total annual rainfall and the seasonality of this rainfall are hypothesised to affect the species' distribution; gaps in its distribution are in areas of higher rainfall that support closed forest and in regions with more winter dominated rainfall.

Arbuscular mycorrhiza mitigates zinc stress on Eucalyptus grandis through regulating metal tolerance protein gene expression and ionome uptake

Arbuscular mycorrhizal (AM) fungi are symbionts of most terrestrial plants and enhance their adaptability in metal-contaminated soils. In this study, mycorrhized and non-mycorrhized Eucalyptus grandis were grown under different Zn treatments. After 6 weeks of treatment, the growing status and ionome content of plants as well as the expression patterns of metal tolerance proteins and auxin biosynthesis-related genes were measured. In this study, mycorrhized E. grandis showed higher biomass and height at a high level of Zn compared with non-mycorrhized plants. In addition, AM plants accumulated P, Mg, and Mn in roots and P, Fe, and Cu in shoots, which indicate that AM fungi facilitate the uptake of ionome nutrients to promote plant growth. In addition, mycorrhiza upregulated the expression of EgMTP1 and EgMTP7, whose encoding proteins were predicted to be located at the vacuolar membrane. Meanwhile, Golgi membrane transporter EgMTP5 was also induced in AM shoot. Our results suggest that AM likely mitigates Zn toxicity through sequestrating excess Zn into vacuolar and Golgi. Furthermore, the expression of auxin biosynthesis-related genes was facilitated by AM, and this is probably another approach for Zn tolerance.

Genetic associations with resistance to Meloidogyne enterolobii in guava (Psidium sp.) using cross-genera SNPs and comparative genomics to Eucalyptus highlight evolutionary conservation across the Myrtaceae

Tropical fruit tree species constitute a yet untapped supply of outstanding diversity of taste and nutritional value, barely developed from the genetics standpoint, with scarce or no genomic resources to tackle the challenges arising in modern breeding practice. We generated a de novo genome assembly of the Psidium guajava, the super fruit “apple of the tropics”, and successfully transferred 14,268 SNP probesets from Eucalyptus to Psidium at the nucleotide level, to detect genomic loci linked to resistance to the root knot nematode (RKN) Meloidogyne enterolobii derived from the wild relative P. guineense. Significantly associated loci with resistance across alternative analytical frameworks, were detected at two SNPs on chromosome 3 in a pseudo-assembly of Psidium guajava genome built using a syntenic path approach with the Eucalyptus grandis genome to determine the order and orientation of the contigs. The P. guineense-derived resistance response to RKN and disease onset is conceivably triggered by mineral nutrients and phytohormone homeostasis or signaling with the involvement of the miRNA pathway. Hotspots of mapped resistance quantitative trait loci and functional annotation in the same genomic region of Eucalyptus provide further indirect support to our results, highlighting the evolutionary conservation of genomes across genera of Myrtaceae in the adaptation to pathogens. Marker assisted introgression of the resistance loci mapped should accelerate the development of improved guava cultivars and hybrid rootstocks.

Chromosome-level genome assembly and functional characterization of terpene synthases provide insights into the volatile terpenoid biosynthesis of Wurfbainia villosa

Wurfbainia villosa is a well-known medicinal and edible plant that is widely cultivated in the Lingnan region of China. Its dried fruits (called Fructus Amomi) are broadly used in traditional Chinese medicine for curing gastrointestinal diseases and are rich in volatile terpenoids. Here, we report a high-quality chromosome-level genome assembly of W. villosa with a total size of approximately 2.80 Gb, 42 588 protein-coding genes, and a very high percentage of repetitive sequences (87.23%). Genome analysis showed that W. villosa likely experienced a recent whole-genome duplication event prior to the W. villosa-Zingiber officinale divergence (approximately 11 million years ago), and a recent burst of long terminal repeat insertions afterward. The W. villosa genome enabled the identification of 17 genes involved in the terpenoid skeleton biosynthesis pathway and 66 terpene synthase (TPS) genes. We found that tandem duplication events have an important contribution to the expansion of WvTPSs, which likely drove the production of volatile terpenoids. In addition, functional characterization of 18 WvTPSs, focusing on the TPS-a and TPS-b subfamilies, showed that most of these WvTPSs are multi-product TPS and are predominantly expressed in seeds. The present study provides insights into the genome evolution and the molecular basis of the volatile terpenoids diversity in W. villosa. The genome sequence also represents valuable resources for the functional gene research and molecular breeding of W. villosa.

Multiple PHT1 family phosphate transporters are recruited for mycorrhizal symbiosis in Eucalyptus grandis and conserved PHT1;4 is a requirement for the arbuscular mycorrhizal symbiosis

Eucalypts engage in a mutualistic endosymbiosis with arbuscular mycorrhizal (AM) fungi to acquire mineral nutrients from soils, particularly inorganic phosphate (Pi). In return, the host plant provides organic carbons to its fungal partners. However, the mechanism by which the Eucalyptus plants acquire Pi released from the AM fungi has remained elusive. In this study, we investigated the characterization of potential PHOSPHATE TRANSPORTER1 (PHT1) family Pi transporters in AM symbiosis in Eucalyptus grandis W. Hill ex Maiden. We show that multiple PHT1 family Pi transporters were recruited for AM symbiosis in E. grandis. We further report that EgPT4, an E. grandis member of the PHT1 family, is conserved across angiosperms and is exclusively expressed in AM roots with arbuscule-containing cells and localizes to the periarbuscular membrane (PAM). EgPT4 was able to complement a yeast mutant strain defective in all inorganic Pi transporters and mediate Pi uptake. Importantly, EgPT4 is essential for improved E. grandis growth, total phosphorus concentration and arbuscule development during symbiosis. Moreover, silencing of EgPT4 led to the induction of polyphosphate accumulation relevant genes of Rhizophagus irregularis DAOM 197198. Collectively, our results unravel a pivotal role for EgPT4 in symbiotic Pi transport across the PAM required for arbuscule development in E. grandis.

A high-quality Buxus austro-yunnanensis (Buxales) genome provides new insights into karyotype evolution in early eudicots

BACKGROUND

Eudicots are the most diverse group of flowering plants that compromise five well-defined lineages: core eudicots, Ranunculales, Proteales, Trochodendrales, and Buxales. However, the phylogenetic relationships between these five lineages and their chromosomal evolutions remain unclear, and a lack of high-quality genome analyses for Buxales has hindered many efforts to address this knowledge gap.

RESULTS

Here, we present a high-quality chromosome-level genome of Buxus austro-yunnanensis (Buxales). Our phylogenomic analyses revealed that Buxales and Trochodendrales are genetically similar and classified as sisters. Additionally, both are sisters to the core eudicots, while Ranunculales was found to be the first lineage to diverge from these groups. Incomplete lineage sorting and hybridization were identified as the main contributors to phylogenetic discordance (34.33%) between the lineages. In fact, B. austro-yunnanensis underwent only one whole-genome duplication event, and collinear gene phylogeny analyses suggested that separate independent polyploidizations occurred in the five eudicot lineages. Using representative genomes from these five lineages, we reconstructed the ancestral eudicot karyotype (AEK) and generated a nearly gapless karyotype projection for each eudicot species. Within core eudicots, we recovered one common chromosome fusion event in asterids and malvids, respectively. Further, we also found that the previously reported fused AEKs in Aquilegia (Ranunculales) and Vitis (core eudicots) have different fusion positions, which indicates that these two species have different karyotype evolution histories.

CONCLUSIONS

Based on our phylogenomic and karyotype evolution analyses, we revealed the likely relationships and evolutionary histories of early eudicots. Ultimately, our study expands genomic resources for early-diverging eudicots.

Terpenoids and their gene regulatory networks in Opisthopappus taihangensis 'Taihang Mingzhu' as detected by transcriptome and metabolome analyses

'Taihang Mingzhu' is the hybrid offspring of Opisthopappus taihangensis, and it has an excellent characteristic of whole-plant fragrance. At present, the genes and metabolites involved in the synthesis of its aromatic compounds are unknown because of the paucity of molecular biology studies on flowering in the Opisthopappus Shih genus. To elucidate the biosynthetic pathways of terpenoids, the main aromatic compounds in 'Taihang Mingzhu', we conducted transcriptome and metabolite analyses on its leaves and bud, inflorescences at the color-development, flowering, and full-bloom stages. A total of 82,685 unigenes were obtained, of which 43,901 were annotated on the basis of information at the protein databases Nr, SwissProt, KEGG, and COG/KOG (e-value<0.00001). Using gas headspace solid-phase microextraction chromatography - mass spectrometry (HS-SPME-GC/MS), 1350 metabolites were identified, the most abundant of which were terpenoids (302 metabolites). Analyses of the gene regulatory network of terpenoids in 'Taihang Mingzhu' identified 52 genes potentially involved in the regulation of terpenoid synthesis. The correlations between genes related to terpenoid metabolism/regulation and metabolite abundance were analyzed. We also extracted the essential oil from the leaves of 'Taihang Mingzhu' by hydrodistillation, and obtained 270 aromatic compounds. Again, the most abundant class was terpenoids. These results provide guidance for the extraction of essential oil from 'Taihang Mingzhu' leaves and flowers. In addition, our analyses provide valuable genetic resources to identify genetic targets to manipulate the aromatic profiles of this plant and other members the Opisthopappus Shih genus by molecular breeding.

Genome-wide identification of calcineurin B-like protein-interacting protein kinase gene family reveals members participating in abiotic stress in the ornamental woody plant Lagerstroemia indica

Calcineurin B-like protein-interacting protein kinases (CIPKs) play important roles in plant responses to stress. However, their function in the ornamental woody plant Lagerstroemia indica is remains unclear. In this study, the LiCIPK gene family was analyzed at the whole genome level. A total of 37 LiCIPKs, distributed across 17 chromosomes, were identified. Conserved motif analysis indicated that all LiCIPKs possess a protein kinase motif (S_TKc) and C-terminal regulatory motif (NAF), while seven LiCIPKs lack a protein phosphatase interaction (PPI) motif. 3D structure analysis further revealed that the N-terminal and C-terminal 3D-structure of 27 members are situated near to each other, while 4 members have a looser structure, and 6 members lack intact structures. The intra- and interspecies collinearity analysis, synonymous substitution rate (K _s ) peaks of duplicated LiCIPKs, revealed that ∼80% of LiCIPKs were retained by the two whole genome duplication (WGD) events that occurred approximately 56.12-61.16 million year ago (MYA) and 16.24-26.34 MYA ago. The promoter of each LiCIPK contains a number of auxin, abscisic acid, gibberellic acid, salicylic acid, and drought, anaerobic, defense, stress, and wound responsive cis-elements. Of the 21 members that were successfully amplified by qPCR, 18 LiCIPKs exhibited different expression patterns under NaCl, mannitol, PEG8000, and ABA treatments. Given that LiCIPK30, the AtSOS2 ortholog, responded to all four types of stress it was selected for functional verification. LiCIPK30 complements the atsos2 phenotype in vivo. 35S:LiCIPK-overexpressing lines exhibit increased leaf area increment, chlorophyll a and b content, reactive oxygen species scavenging enzyme activity, and expression of ABF3 and RD22, while the degree of membrane lipid oxidation decreases under NaCl treatment compared to WT. The evolutionary history, and potential mechanism by which LiCIPK30 may regulate plant tolerance to salt stress were also discussed. In summary, we identified LiCIPK members involved in abiotic stress and found that LiCIPK30 transgenic Arabidopsis exhibits more salt and osmotic stress tolerance than WT. This research provides a theoretical foundation for further investigation into the function of LiCIPKs, and for mining gene resources to facilitate the cultivation and breeding of new L. indica varieties in coastal saline-alkali soil.

Reference quality genome sequence of Indian pomegranate cv. 'Bhagawa' (Punica granatum L.)

Pomegranate is an important fruit crop for ensuring livelihood and nutrition security in fragile semi-arid regions of the globe having limited irrigation resources. This is a high-value, nutritionally rich, and export-oriented agri-commodity that ensures high returns on investment to growers across the world. Although it is a valuable fruit crop, it has received only a limited genomics research outcome. To fast-track the pomegranate improvement program, de novo whole-genome sequencing of the main Indian cultivar 'Bhagawa' was initiated by the Indian Council of Agricultural Research-National Research Center on Pomegranate (ICAR-NRCP). We have demonstrated that a combination of commercially available technologies from Illumina, PacBio, 10X Genomics, and BioNano Genomics could be used efficiently for sequencing and reference-grade de novo assembly of the pomegranate genome. The research led to a final reference-quality genome assembly for 'Bhagawa' of 346.08 Mb in 342 scaffolds and an average N50 of 16.12 Mb and N90 of 1088.62 Kb. This assembly covered more than 98% of the estimated pomegranate genome size, 352.54 Mb. The LTR assembly index (LAI) value of 10 and 93.68% Benchmarking Universal Single-Copy Orthologs (BUSCO) completeness score over the 1,440 ortholog genes of the completed pomegranate genome indicates the quality of the assembled pomegranate genome. Furthermore, 29,435 gene models were discovered with a mean transcript length of 2,954 bp and a mean coding sequence length 1,090 bp. Four transcript data samples of pomegranate tissues were mapped over the assembled 'Bhagawa' genome up to 95% significant matches, indicating the high quality of the assembled genome. We have compared the 'Bhagawa' genome with the genomes of the pomegranate cultivars 'Dabenzi' and 'Taishanhong.' We have also performed whole-genome phylogenetic analysis using Computational Analysis of Gene Family Evolution (CAFE) and found that Eucalyptus grandis and pomegranate diverged 64 (60-70) million years ago. About 1,573 protein-coding resistance genes identified in the 'Bhagawa' genome were classified into 32 domains. In all, 314 copies of miRNA belonging to 26 different families were identified in the 'Bhagawa' genome. The reference-quality genome assembly of 'Bhagawa' is certainly a significant genomic resource for accelerated pomegranate improvement.

Low phosphorus induces differential metabolic responses in eucalyptus species improving nutrient use efficiency

Phosphorus (P) is a vital nutrient for plant growth. P availability is generally low in soils, and plant responses to low P availability need to be better understood. In a previous study, we studied the growth and physiological responses of 24 species to low P availability in the soil and verified of eucalypts, five (Eucalyptus acmenoides, E. grandis, E. globulus, E. tereticornis, and Corymbia maculata) contrasted regarding their efficiency and responsiveness to soil P availability. Here, we obtained the metabolomic and lipidomic profile of leaves, stems, and roots from these species growing under low (4.5 mg dm^-3) and sufficient (10.8 mg dm^-3) P in the soil. Disregarding the level of P in the soils, P allocation was always higher in the stems. However, when grown in the P-sufficient soil, the stems steadily were the largest compartment of the total plant P. Under low P, the relative contents of primary metabolites, such as amino acids, TCA cycle intermediates, organic acids and carbohydrates, changed differently depending on the species. Additionally, phosphorylated metabolites showed enhanced turnover or reductions. While photosynthetic efficiencies were not related to higher biomass production, A/Ci curves showed that reduced P availability increased the eucalypt species' Vcmax, Jmax and photosynthetic P-use efficiency. Plants of E. acmenoides increased galactolipids and sulfolipids in leaves more than other eucalypt species, suggesting that lipid remodelling can be a strategy to cope with the P shortage in this species. Our findings offer insights to understand genotypic efficiency among eucalypt species to accommodate primary metabolism under low soil P availability and eventually be used as biochemical markers for breeding programs.

Genomic insights into rapid speciation within the world's largest tree genus Syzygium

Species radiations, despite immense phenotypic variation, can be difficult to resolve phylogenetically when genetic change poorly matches the rapidity of diversification. Genomic potential furnished by palaeopolyploidy, and relative roles for adaptation, random drift and hybridisation in the apportionment of genetic variation, remain poorly understood factors. Here, we study these aspects in a model radiation, Syzygium, the most species-rich tree genus worldwide. Genomes of 182 distinct species and 58 unidentified taxa are compared against a chromosome-level reference genome of the sea apple, Syzygium grande. We show that while Syzygium shares an ancient genome doubling event with other Myrtales, little evidence exists for recent polyploidy events. Phylogenomics confirms that Syzygium originated in Australia-New Guinea and diversified in multiple migrations, eastward to the Pacific and westward to India and Africa, in bursts of speciation visible as poorly resolved branches on phylogenies. Furthermore, some sublineages demonstrate genomic clines that recapitulate cladogenetic events, suggesting that stepwise geographic speciation, a neutral process, has been important in Syzygium diversification.

Exploring leaf hydraulic traits to predict drought tolerance of Eucalyptus clones

Ongoing changes in climate, and the consequent mortality of natural and cultivated forests across the globe, highlight the urgent need to understand the plant traits associated with greater tolerance to drought. Here, we aimed at assessing key foliar traits, with a focus on the hydraulic component, that could confer a differential ability to tolerate drought in three commercial hybrids of the most important Eucalyptus species utilized in tropical silviculture: E. urophyla, E. grandis and E. camaldulensis. All genotypes exhibited similar water potential when the 90% stomatal closure (Ψgs90) occurs with Ψgs90 always preceding the start of embolism events. The drought-tolerant hybrid showed a higher leaf resistance to embolism, but the leaf hydraulic efficiency was similar among all genotypes. Other traits presented by the drought-tolerant hybrid were a higher cell wall reinforcement, lower value of osmotic potential at full turgor and greater bulk modulus of elasticity. We also identified that the leaf capacitance after the turgor loss, the ratio between cell wall thickness (t) and lumen breadth (b) ratio (t/b)3, and the minimal conductance might be good proxies for screening drought-tolerant Eucalyptus genotypes. Our findings suggest that xylem resistance to embolism can be an important component of drought tolerance in Eucalyptus in addition to other traits aimed at delaying the development of high tensions in the xylem. Highlight Drought tolerance in tropical Eucalyptus hybrids encompasses a high leaf resistance to embolism and a suite of traits aimed at delaying the development of high tensions in the xylem.

Genomic consequences of artificial selection during early domestication of a wood fibre crop

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

The Eucalyptus grandis chloroplast proteome: Seasonal variations in leaf development

Chloroplast metabolism is very sensitive to environmental fluctuations and is intimately related to plant leaf development. Characterization of the chloroplast proteome dynamics can contribute to a better understanding on plant adaptation to different climate scenarios and leaf development processes. Herein, we carried out a discovery-driven analysis of the Eucalyptus grandis chloroplast proteome during leaf maturation and throughout different seasons of the year. The chloroplast proteome from young leaves differed the most from all assessed samples. Most upregulated proteins identified in mature and young leaves were those related to catabolic-redox signaling and biogenesis processes, respectively. Seasonal dynamics revealed unique proteome features in the fall and spring periods. The most abundant chloroplast protein in humid (wet) seasons (spring and summer) was a small subunit of RuBisCO, while in the dry periods (fall and winter) the proteins that showed the most pronounced accumulation were associated with photo-oxidative damage, Calvin cycle, shikimate pathway, and detoxification. Our investigation of the chloroplast proteome dynamics during leaf development revealed significant alterations in relation to the maturation event. Our findings also suggest that transition seasons induced the most pronounced chloroplast proteome changes over the year. This study contributes to a more comprehensive understanding on the subcellular mechanisms that lead to plant leaf adaptation and ultimately gives more insights into Eucalyptus grandis phenology.

Genome-Wide Identification of Eucalyptus Heat Shock Transcription Factor Family and Their Transcriptional Analysis under Salt and Temperature Stresses

Heat shock transcription factors (HSFs) activate heat shock protein gene expression by binding their promoters in response to heat stress and are considered to be pivotal transcription factors in plants. Eucalyptus is a superior source of fuel and commercial wood. During its growth, high temperature or other abiotic stresses could impact its defense capability and growth. Hsf genes have been cloned and sequenced in many plants, but rarely in Eucalyptus. In this study, we used bioinformatics methods to analyze and identify Eucalyptus Hsf genes, their chromosomal localization and structure. The phylogenetic relationship and conserved domains of their encoded proteins were further analyzed. A total of 36 Hsf genes were identified and authenticated from Eucalyptus, which were scattered across 11 chromosomes. They could be classified into three classes (A, B and C). Additionally, a large number of stress-related cis-regulatory elements were identified in the upstream promoter sequence of HSF, and cis-acting element analysis indicated that the expression of EgHsf may be regulated by plant growth and development, metabolism, hormones and stress responses. The expression profiles of five representative Hsf genes, EgHsf4, EgHsf9, EgHsf13, EgHsf24 and EgHsf32, under salt and temperature stresses were examined by qRT-PCR. The results show that the expression pattern of class B genes (EgHsf4, EgHsf24 and EgHsf32) was more tolerant to abiotic stresses than that of class A genes (EgHsf9 and EgHsf13). However, the expressions of all tested Hsf genes in six tissues were at different levels. Finally, we investigated the network of interplay between genes, and the results suggest that there may be synergistic effects between different Hsf genes in response to abiotic stresses. We conclude that the Hsf gene family played an important role in the growth and developmental processes of Eucalyptus and could be vital for maintaining cell homeostasis against external stresses. This study provides basic information on the members of the Hsf gene family in Eucalyptus and lays the foundation for the functional identification of related genes and the further investigation of their biological functions in plant stress regulation.

Genetic control of the operculum and capsule morphology of Eucalyptus globulus

BACKGROUND AND AIMS

The petaline operculum that covers the inner whorls until anthesis and the woody capsule that develops after fertilization are reproductive structures of eucalypts that protect the flower and seeds. Although they are distinct organs, they both develop from flower buds and this common ontogeny suggests shared genetic control. In Eucalyptus globulus their morphology is variable and we aimed to identify the quantitative trait loci (QTL) underlying this variation and determine whether there is common genetic control of these ecologically and taxonomically important reproductive structures.

METHODS

Samples of opercula and capsules were collected from 206 trees that belong to a large outcrossed F2E. globulus mapping population. The morphological variation in these structures was characterized by measuring six operculum and five capsule traits. QTL analysis was performed using these data and a linkage map consisting of 480 markers.

KEY RESULTS

A total of 27 QTL were detected for operculum traits and 28 for capsule traits, with the logarithm of odds ranging from 2.8 to 11.8. There were many co-located QTL associated with operculum or capsule traits, generally reflecting allometric relationships. A key finding was five genomic regions where co-located QTL affected both operculum and capsule morphology, and the overall trend for these QTL was to affect elongation of both organs. Some of these QTL appear to have a significant effect on the phenotype, with the strongest QTL explaining 26.4 % of the variation in operculum shape and 16.4 % in capsule shape. Flower bud measurements suggest the expression of these QTL starts during bud development. Several candidate genes were found associated with the QTL and their putative function is discussed.

CONCLUSIONS

Variation in both operculum and capsule traits in E. globulus is under strong genetic control. Our results suggest that these reproductive structures share a common genetic pathway during flower bud development.

Xylem Transcriptome Analysis in Contrasting Wood Phenotypes of Eucalyptus urophylla × tereticornis Hybrids

An investigation of the effects of two important post-transcriptional regulatory mechanisms, gene transcription and alternative splicing (AS), on the wood formation of Eucalyptusurophylla × tereticornis, an economic tree species widely planted in southern China, was carried out. We performed RNA-seq on E. urophylla × tereticornis hybrids with highly contrasting wood basic density (BD), cellulose content (CC), hemicellulose content (HC), and lignin content (LC). Signals of strong differentially expressed genes (DEGs) and differentially spliced genes (DSGs) were detected in all four groups of wood properties, suggesting that gene transcription and selective splicing may have important regulatory roles in wood properties. We found that there was little overlap between DEGs and DSGs in groups of the same trait. Furthermore, the key DEGs and DSGs that were detected simultaneously in the four groups tended to be enriched in different Gene Ontology terms, Kyoto Encyclopedia of Genes and Genomes pathways, and transcription factors. These results implied that regulation of gene transcription and AS is controlled by independent regulatory systems in wood formation. Lastly, we detected transcript levels of known wood biosynthetic genes and found that 79 genes encoding mainly enzymes or proteins such as UGT, LAC, CAD, and CESA may be involved in the positive or negative regulation of wood properties. This study reveals potential molecular mechanisms that may regulate wood formation and will contribute to the genetic improvement of Eucalyptus.

The clove (Syzygium aromaticum) genome provides insights into the eugenol biosynthesis pathway

The clove (Syzygium aromaticum) is an important tropical spice crop in global trade. Evolving environmental pressures necessitate modern characterization and selection techniques that are currently inaccessible to clove growers owing to the scarcity of genomic and genetic information. Here, we present a 370-Mb high-quality chromosome-scale genome assembly for clove. Comparative genomic analysis between S. aromaticum and Eucalyptus grandis—both species of the Myrtaceae family—reveals good genome structure conservation and intrachromosomal rearrangements on seven of the eleven chromosomes. We report genes that belong to families involved in the biosynthesis of eugenol, the major bioactive component of clove products. On the basis of our transcriptomic and metabolomic findings, we propose a hypothetical scenario in which eugenol acetate plays a key role in high eugenol accumulation in clove leaves and buds. The clove genome is a new contribution to omics resources for the Myrtaceae family and an important tool for clove research.

A newly assembled clove genome is compared with E. grandis to investigate genome evolution between the two genera of the Myrtaceae family, and putative genes involved in the biosynthesis of eugenol are identified through transcriptomics and metabolomics.

Convergent gene clusters underpin hyperforin biosynthesis in St John's wort

The meroterpenoid hyperforin is responsible for the antidepressant activity of St John's wort extracts, but the genes controlling its biosynthesis are unknown. Using genome mining and biochemical work, we characterize two biosynthetic gene clusters (BGCs) that encode the first three steps in the biosynthesis of hyperforin precursors. The findings of syntenic and phylogenetic analyses reveal the parallel assembly of the two BGCs. The syntenous BGC in Mesua ferrea indicates that the first cluster was assembled before the divergence of the Hypericaceae and Calophyllaceae families. The assembly of the second cluster is the result of a coalescence of genomic fragments after a major duplication event. The differences between the two BGCs - in terms of gene expression, response to methyl jasmonate, substrate specificity and subcellular localization of key enzymes - suggest that the presence of the two clusters could serve to generate separate pools of precursors. The parallel assembly of two BGCs with similar compositions in a single plant species is uncommon, and our work provides insights into how and when these gene clusters form. Our discovery helps to advance our understanding of the evolution of plant specialized metabolism and its genomic organization. Additionally, our results offer a foundation from which hyperforin biosynthesis can be more fully understood, and which can be used in future metabolic engineering applications.

The Larix kaempferi genome reveals new insights into wood properties

Here, through single-molecule real-time sequencing, we present a high-quality genome sequence of the Japanese larch (Larix kaempferi), a conifer species with great value for wood production and ecological afforestation. The assembled genome is 10.97 Gb in size, harboring 45,828 protein-coding genes. Of the genome, 66.8% consists of repeat sequences, of which long terminal repeat retrotransposons are dominant and make up 69.86%. We find that tandem duplications have been responsible for the expansion of genes involved in transcriptional regulation and stress responses, unveiling their crucial roles in adaptive evolution. Population transcriptome analysis reveals that lignin content in L. kaempferi is mainly determined by the process of monolignol polymerization. The expression values of six genes (LkCOMT7, LkCOMT8, LkLAC23, LkLAC102, LkPRX148, and LkPRX166) have significantly positive correlations with lignin content. These results indicated that the increased expression of these six genes might be responsible for the high lignin content of the larches' wood. Overall, this study provides new genome resources for investigating the evolution and biological function of conifer trees, and also offers new insights into wood properties of larches.

Genome structure and evolutionary history of frankincense producing Boswellia sacra

Boswellia sacra Flueck (family Burseraceae) tree is wounded to produce frankincense. We report its de novo assembled genome (667.8 Mb) comprising 18,564 high-confidence protein-encoding genes. Comparing conserved single-copy genes across eudicots suggest >97% gene space assembly of B. sacra genome. Evolutionary history shows B. sacra gene-duplications derived from recent paralogous events and retained from ancient hexaploidy shared with other eudicots. The genome indicated a major expansion of Gypsy retroelements in last 2 million years. The B. sacra genetic diversity showed four clades intermixed with a primary genotype—dominating most resin-productive trees. Further, the stem transcriptome revealed that wounding concurrently activates phytohormones signaling, cell wall fortification, and resin terpenoid biosynthesis pathways leading to the synthesis of boswellic acid—a key chemotaxonomic marker of Boswellia. The sequence datasets reported here will serve as a foundation to investigate the genetic determinants of frankincense and other resin-producing species in Burseraceae.

•

Assembly and architecture of frankincense producing Boswellia sacra Flueck

•

Comparative genomics and evolutionary history of frankincense tree within orders

•

Transcriptome of stem part and gene expression patterns of wounding to the tree

•

Resin biosynthesis pathway and related CYP450 enzymes and gene families