Distribution and in situ conservation of a relic Chinese oil woody species Xanthoceras sorbifolium (yellowhorn)

Agronomic, physiological and transcriptional characteristics provide insights into fatty acid biosynthesis in yellowhorn (Xanthoceras sorbifolium Bunge) during fruit ripening

Yellowhorn (Xanthoceras sorbifolium Bunge) is an oil-bearing tree species in northern China. In this study, we used yellowhorn from Heilongjiang to analyze the morphological and physiological changes of fruit development and conducted transcriptome sequencing. The results showed that the fruit experienced relatively slow growth from fertilization to DAF20 (20 days after flowering). From DAF40 to DAF60, the fruit entered an accelerated development stage, with a rapid increase in both transverse and longitudinal diameters, and the kernel contour developed completely at DAF40. From DAF60 to DAF80, the transverse and vertical diameters of the fruit developed slowly, and the overall measures remained stable until maturity. The soluble sugar, starch, and anthocyanin content gradually accumulated until reaching a peak at DAF80 and then rapidly decreased. RNA-seq analysis revealed differentially expressed genes (DEGs) in the seed coat and kernel, implying that seed components have different metabolite accumulation mechanisms. During the stages of seed kernel development, k-means clustering separated the DEGs into eight sub-classes, indicating gene expression shifts during the fruit ripening process. In subclass 8, the fatty acid biosynthesis pathway was enriched, suggesting that this class was responsible for lipid accumulation in the kernel. WGCNA revealed ten tissue-specific modules for the 12 samples among 20 modules. We identified 54 fatty acid biosynthesis pathway genes across the genome, of which 14 was quantified and confirmed by RT-qPCR. Most genes in the plastid synthesis stage showed high expression during the DAF40-DAF60 period, while genes in the endoplasmic reticulum synthesis stage showed diverse expression patterns. EVM0012847 (KCS) and EVM0002968 (HCD) showed similar high expression in the early stages and low expression in the late stages. EVM0022385 (HCD) exhibited decreased expression from DAF40 to DAF60 and then increased from DAF60 to DAF100. EVM0000575 (KCS) was increasingly expressed from DAF40 to DAF60 and then decreased from DAF60 to DAF100. Finally, we identified transcription factors (TFs) (HB-other, bHLH and ARF) that were predicted to bind to fatty acid biosynthesis pathway genes with significant correlations. These results are conducive to promoting the transcriptional regulation of lipid metabolism and the genetic improvement in terms of high lipid content of yellowhorn.

Differentiated responses of the phyllosphere bacterial community of the yellowhorn tree to precipitation and temperature regimes across Northern China

Introduction

As an ephemeral and oligotrophic environment, the phyllosphere harbors many highly diverse microorganisms. Importantly, it is known that their colonization of plant leaf surfaces is considerably influenced by a few abiotic factors related to climatic conditions. Yet how the dynamics of phyllosphere bacterial community assembly are shaped by detailed climatological elements, such as various bioclimatic variables, remains poorly understood.

Methods

Using high-throughput 16S rRNA gene amplicon sequencing technology, we analyzed the bacterial communities inhabiting the leaf surfaces of an oilseed tree, yellowhorn (Xanthoceras sorbifolium), grown at four sites (Yinchuan, Otogqianqi, Tongliao, and Zhangwu) whose climatic status differs in northern China.

Results and Discussion

We found that the yellowhorn phyllosphere's bacterial community was generally dominated by four phyla: Proteobacteria, Firmicutes, Actinobacteria, and Bacteroidetes. Nevertheless, bacterial community composition differed significantly among the four sampled site regions, indicating the possible impact of climatological factors upon the phyllosphere microbiome. Interestingly, we also noted that the α-diversities of phyllosphere microbiota showed strong positive or negative correlation with 13 bioclimatic factors (including 7 precipitation factors and 6 temperature factors). Furthermore, the relative abundances of 55 amplicon sequence variants (ASVs), including three ASVs representing two keystone taxa (the genera Curtobacterium and Streptomyces), exhibited significant yet contrary responses to the precipitation and temperature climatic variables. That pattern was consistent with all ASVs' trends of possessing opposite correlations to those two parameter classes. In addition, the total number of links and nodes, which conveys community network complexity, increased with rising values of most temperature variables. Besides that, remarkably positive relevance was found between average clustering coefficient and most precipitation variables. Altogether, these results suggest the yellowhorn phyllosphere bacterial community is capable of responding to variation in rainfall and temperature regimes in distinctive ways.

Centromere-Specific Retrotransposons and Very-Long-Chain Fatty Acid Biosynthesis in the Genome of Yellowhorn (Xanthoceras sorbifolium, Sapindaceae), an Oil-Producing Tree With Significant Drought Resistance

In-depth genome characterization is still lacking for most of biofuel crops, especially for centromeres, which play a fundamental role during nuclear division and in the maintenance of genome stability. This study applied long-read sequencing technologies to assemble a highly contiguous genome for yellowhorn (Xanthoceras sorbifolium), an oil-producing tree, and conducted extensive comparative analyses to understand centromere structure and evolution, and fatty acid biosynthesis. We produced a reference-level genome of yellowhorn, ∼470 Mb in length with ∼95% of contigs anchored onto 15 chromosomes. Genome annotation identified 22,049 protein-coding genes and 65.7% of the genome sequence as repetitive elements. Long terminal repeat retrotransposons (LTR-RTs) account for ∼30% of the yellowhorn genome, which is maintained by a moderate birth rate and a low removal rate. We identified the centromeric regions on each chromosome and found enrichment of centromere-specific retrotransposons of LINE1 and Gypsy in these regions, which have evolved recently (∼0.7 MYA). We compared the genomes of three cultivars and found frequent inversions. We analyzed the transcriptomes from different tissues and identified the candidate genes involved in very-long-chain fatty acid biosynthesis and their expression profiles. Collinear block analysis showed that yellowhorn shared the gamma (γ) hexaploidy event with Vitis vinifera but did not undergo any further whole-genome duplication. This study provides excellent genomic resources for understanding centromere structure and evolution and for functional studies in this important oil-producing plant.

Small RNA profiling for identification of microRNAs involved in regulation of seed development and lipid biosynthesis in yellowhorn

BACKGROUND

Yellowhorn (Xanthoceras sorbifolium), an endemic woody oil-bearing tree, has become economically important and is widely cultivated in northern China for bioactive oil production. However, the regulatory mechanisms of seed development and lipid biosynthesis affecting oil production in yellowhorn are still elusive. MicroRNAs (miRNAs) play crucial roles in diverse aspects of biological and metabolic processes in seeds, especially in seed development and lipid metabolism. It is still unknown how the miRNAs regulate the seed development and lipid biosynthesis in yellowhorn.

RESULTS

Here, based on investigations of differences in the seed growth tendency and embryo oil content between high-oil-content and low-oil-content lines, we constructed small RNA libraries from yellowhorn embryos at four seed development stages of the two lines and then profiled small RNA expression using high-throughput sequencing. A total of 249 known miRNAs from 46 families and 88 novel miRNAs were identified. Furthermore, by pairwise comparisons among the four seed development stages in each line, we found that 64 miRNAs (53 known and 11 novel miRNAs) were differentially expressed in the two lines. Across the two lines, 15, 11, 10, and 7 differentially expressed miRNAs were detected at 40, 54, 68, and 81 days after anthesis, respectively. Bioinformatic analysis was used to predict a total of 2654 target genes for 141 differentially expressed miRNAs (120 known and 21 novel miRNAs). Most of these genes were involved in the fatty acid biosynthetic process, regulation of transcription, nucleus, and response to auxin. Using quantitative real-time PCR and an integrated analysis of miRNA and mRNA expression, miRNA-target regulatory modules that may be involved in yellowhorn seed size, weight, and lipid biosynthesis were identified, such as miR172b-ARF2 (auxin response factor 2), miR7760-p3_1-AGL61 (AGAMOUS-LIKE 61), miR319p_1-FAD2-2 (omega-6 fatty acid desaturase 2-2), miR5647-p3_1-DGAT1 (diacylglycerol acyltransferase 1), and miR7760-p5_1-MED15A (Mediator subunit 15a).

CONCLUSIONS

This study provides new insights into the important regulatory roles of miRNAs in the seed development and lipid biosynthesis in yellowhorn. Our results will be valuable for dissecting the post-transcriptional and transcriptional regulation of seed development and lipid biosynthesis, as well as improving yellowhorn in northern China.

Role of Xanthoceras sorbifolium MYB44 in tolerance to combined drought and heat stress via modulation of stomatal closure and ROS homeostasis

Yellowhorn (Xanthoceras sorbifolium) is an important edible woody oil tree species that is endemic to China. Drought and heat stresses are factors severely limiting the high-quality development of the yellowhorn industry. Transcription factors (TFs) play critical roles in regulating the response of woody plant species to water deficit or high temperature. However, the MYB TFs that respond to combined drought and heat stress in yellowhorn remain unclear. Here, we first investigated the physiological changes in 5 yellowhorn varieties in response to combined stress treatments. We observed significant changes in antioxidant enzyme activities and photosynthesis. The Maigaiti variety yielded the best results and was selected for subsequent experiments. An R2R3-type MYB TF, designated XsMYB44, was isolated from the leaves of yellowhorn. XsMYB44 expression was strongly induced by combined stress. Suppression of XsMYB44 expression via virus-induced gene silencing weakened yellowhorn tolerance to both individual and combined drought and heat stress, and the increased susceptibility was coupled with decreased plant height, fresh weight and relative water content and inhibited stomatal closure. Moreover, compared with the individual stresses, the combined stress caused increased reactive oxygen species levels and decreased antioxidant enzyme activities and proline content in XsMYB44-silenced plants. Furthermore, the expression levels of several defense-related genes were reduced in the XsMYB44-silenced plants. Overall, we studied the physiological characteristics of 5 yellowhorn varieties, and the results demonstrated that XsMYB44 acts as a positive regulator in the yellowhorn response to combined stress by triggering stomatal closure to maintain water levels and by modulating ROS homeostasis.

Yellowhorn drought-induced transcription factor XsWRKY20 acts as a positive regulator in drought stress through ROS homeostasis and ABA signaling pathway

Yellowhorn (Xanthoceras sorbifolium) is a peculiar woody edible oil-bearing tree in China. WRKY transcription factors have specific roles in plant multiple abiotic stress responses. However, it is still not clear that the molecular mechanisms of WRKYs involve in drought tolerance in yellowhorn. In this study, we isolated a drought-induced group I WRKY gene from yellowhorn, designated as XsWRKY20. Expression of XsWRKY20 was strongly induced by PEG6000, NaCl, ABA and SA. Virus-induced gene silencing (VIGS) of XsWRKY20 reduced tolerance to drought stress in yellowhorn, as determined through physiological analyses of POD activity, SOD activity and proline content. This susceptibility was coupled with decreased expression of stress-related genes. In contrast, overexpression of XsWRKY20 in tobacco notably improved drought tolerance. Compared with the WT plants, the XsWRKY20-transgenic lines exhibited lower ROS and MDA content and higher antioxidant enzyme activity and proline content after drought treatment. Moreover, overexpression of XsWRKY20 enhanced the expression of several genes associated with encoding these antioxidant enzymes, proline biosynthesis and ABA signaling pathway. Taken together, XsWRKY20 functions as a positive regulator contributing to drought stress tolerance through either ROS homeostasis by antioxidant systems or ABA-dependent/independent gene expression pathway.

Transcriptome and Co-Expression Network Analyses Identify the Molecular Signatures Underlying Drought Resistance in Yellowhorn

Drought is a key factor that limits plant growth and yield. Yellowhorn is an important and vigorously promoted oil tree in China. It can survive under certain drought conditions, but a lack of water severely restricts its growth and results in yield losses in arid and semi-arid areas. Therefore, it is important to identify the key pathways and genes to understand the mechanisms of its drought resistance. In this study, we evaluated drought resistance in four types of yellowhorn, and obtained 2669 and 2451 differentially expressed genes (DEGs) via the transcriptome analysis of the comparison of water-saving/water-consuming and fast-growing/slow-growing yellowhorn, respectively, under long-term drought conditions. The gene ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) analyses of DEGs showed the key biological processes and metabolic pathways involved in drought resistance, which demonstrated that there are both the same and different biological processes involved in regulating water use efficiency (WUE) and growth in response to drought stress. Furthermore, the network analysis indicated hub genes (especially seven protein kinases) and potential co-expressed gene clusters in a greenyellow module associated with WUE and a blue module associated with growth. These identified hub genes and key biological processes can significantly enhance our knowledge about the molecular mechanisms of drought resistance in yellowhorn.

Pseudomolecule-level assembly of the Chinese oil tree yellowhorn (Xanthoceras sorbifolium) genome

BACKGROUND

Yellowhorn (Xanthoceras sorbifolium) is a species of the Sapindaceae family native to China and is an oil tree that can withstand cold and drought conditions. A pseudomolecule-level genome assembly for this species will not only contribute to understanding the evolution of its genes and chromosomes but also bring yellowhorn breeding into the genomic era.

FINDINGS

Here, we generated 15 pseudomolecules of yellowhorn chromosomes, on which 97.04% of scaffolds were anchored, using the combined Illumina HiSeq, Pacific Biosciences Sequel, and Hi-C technologies. The length of the final yellowhorn genome assembly was 504.2 Mb with a contig N50 size of 1.04 Mb and a scaffold N50 size of 32.17 Mb. Genome annotation revealed that 68.67% of the yellowhorn genome was composed of repetitive elements. Gene modelling predicted 24,672 protein-coding genes. By comparing orthologous genes, the divergence time of yellowhorn and its close sister species longan (Dimocarpus longan) was estimated at ∼33.07 million years ago. Gene cluster and chromosome synteny analysis demonstrated that the yellowhorn genome shared a conserved genome structure with its ancestor in some chromosomes.

CONCLUSIONS

This genome assembly represents a high-quality reference genome for yellowhorn. Integrated genome annotations provide a valuable dataset for genetic and molecular research in this species. We did not detect whole-genome duplication in the genome. The yellowhorn genome carries syntenic blocks from ancient chromosomes. These data sources will enable this genome to serve as an initial platform for breeding better yellowhorn cultivars.

The genome assembly and annotation of yellowhorn (Xanthoceras sorbifolium Bunge)

BACKGROUND

Yellowhorn (Xanthoceras sorbifolium Bunge), a deciduous shrub or small tree native to north China, is of great economic value. Seeds of yellowhorn are rich in oil containing unsaturated long-chain fatty acids that have been used for producing edible oil and nervonic acid capsules. However, the lack of a high-quality genome sequence hampers the understanding of its evolution and gene functions.

FINDINGS

In this study, a whole genome of yellowhorn was sequenced and assembled by integration of Illumina sequencing, Pacific Biosciences single-molecule real-time sequencing, 10X Genomics linked reads, Bionano optical maps, and Hi-C. The yellowhorn genome assembly was 439.97 Mb, which comprised 15 pseudo-chromosomes covering 95.42% (419.84 Mb) of the assembled genome. The repetitive fractions accounted for 56.39% of the yellowhorn genome. The genome contained 21,059 protein-coding genes. Of them, 18,503 (87.86%) genes were found to be functionally annotated with ≥1 "annotation" term by searching against other databases. Transcriptomic analysis showed that 341, 135, 125, 113, and 100 genes were specifically expressed in hermaphrodite flower, staminate flower, young fruit, leaf, and shoot, respectively. Phylogenetic analysis suggested that yellowhorn and Dimocarpus longan diverged from their most recent common ancestor ∼46 million years ago.

CONCLUSIONS

The availability and subsequent annotation of the yellowhorn genome, as well as the identification of tissue-specific functional genes, provides a valuable reference for plant comparative genomics, evolutionary studies, and molecular design breeding.

Conservation of genetic diversity hotspots of the high-valued relic yellowhorn (Xanthoceras sorbifolium) considering climate change predictions

Genetic structure and major climate factors may contribute to the distribution of genetic diversity of a highly valued oil tree species Xanthoceras sorbifolium (yellowhorn). Long-term over utilization along with climate change is affecting the viability of yellowhorn wild populations. To preserve the species known and unknown valuable gene pools, the identification of genetic diversity "hotspots" is a prerequisite for their consideration as in situ conservation high priority. Chloroplast DNA (cpDNA) diversity was high among 38 natural populations (H d = 0.717, K = 4.616, Tajmas' D = -0.22) and characterized by high genetic divergence (F ST = 0.765) and relatively low gene flow (N m = 0.03), indicating populations isolation reflecting the species' habitat fragmentation and inbreeding depression. Six out of the studied 38 populations are defined as genetic diversity "hotspots." The number and geographic direction of cpDNA mutation steps supported the species southwest to northeast migration history. Climatic factors such as extreme minimum temperature over 30 years indicated that the identified genetic "hotspots" are expected to experience 5°C temperature increase in next following 50 years. The results identified vulnerable genetic diversity "hotspots" and provided fundamental information for the species' future conservation and breeding activities under the anticipated climate change. More specifically, the role of breeding as a component of a gene resource management strategy aimed at fulfilling both utilization and conservation goals.