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

3-ketoacyl-CoA synthase 7 from Xanthoceras sorbifolium seeds is a crucial regulatory enzyme for nervonic acid biosynthesis

Nervonic acid (C24:1) is a very-long-chain fatty acid that plays an imperative role in human brain development and other health benefits. In plants, 3-ketoacyl-CoA synthase (KCS) is the key rate-limiting enzyme for C24:1 biosynthesis. Xanthoceras sorbifolium is a valuable oil-producing economic woody species with abundant C24:1 in seed oils, but the key KCS gene responsible for C24:1 accumulation remains unknown. In this work, a correlation analysis between the transcript profiles of KCS and dynamic change of C24:1 content in developing seeds of X. sorbifolium were conducted to screen out three members of KCS, namely XsKCS4, XsKCS7 and XsKCS8, potentially involved in C24:1 biosynthesis. Of which, the XsKCS7 was highly expressed in developing seeds, while XsKCS4 and XsKCS8 displayed the highest expression in fruits and flowers, respectively. Overexpression of XsKCS4, XsKCS7 and XsKCS8 in yeast Saccharomyces cerevisiae and plant Arabidopsis thaliana indicated that only XsKCS7 possessed the ability to facilitate the biosynthesis of C24:1. These findings collectively suggested that XsKCS7 played a crucial role in specific regulation of C24:1 biosynthesis in X. sorbifolium seeds.

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.

Genome-wide association analysis identifies a candidate gene controlling seed size and yield in Xanthoceras sorbifolium Bunge

Yellow horn (Xanthoceras sorbifolium Bunge) is a woody oilseed tree species whose seed oil is rich in unsaturated fatty acids and rare neuronic acids, and can be used as a high-grade edible oil or as a feedstock for biodiesel production. However, the genetic mechanisms related to seed yield in yellow horn are not well elucidated. This study identified 2 164 863 SNP loci based on 222 genome-wide resequencing data of yellow horn germplasm. We conducted genome-wide association study (GWAS) analysis on three core traits (hundred-grain weight, single-fruit seed mass, and single-fruit seed number) that influence seed yield for the years 2022 and 2020, and identified 399 significant SNP loci. Among these loci, the Chr10_24013014 and Chr10_24012613 loci caught our attention due to their consistent associations across multiple analyses. Through Sanger sequencing, we validated the genotypes of these two loci across 16 germplasms, confirming their consistency with the GWAS analysis results. Downstream of these two significant loci, we identified a candidate gene encoding an AP2 transcription factor protein, which we named XsAP2. RT-qPCR analysis revealed high expression of the XsAP2 gene in seeds, and a significant negative correlation between its expression levels and seed hundred-grain weight, as well as single-fruit seed mass, suggesting its potential role in the normal seed development process. Transgenic Arabidopsis lines with the overexpressed XsAP2 gene exhibited varying degrees of reduction in seed size, number of seeds per silique, and number of siliques per plant compared with wild-type Arabidopsis. Combining these results, we hypothesize that the XsAP2 gene may have a negative regulatory effect on seed yield of yellow horn. These results provide a reference for the molecular breeding of high-yielding yellow horn.

Effects of drought stress on photosynthetic physiological characteristics, leaf microstructure, and related gene expression of yellow horn

Yellow horn grows in northern China and has a high tolerance to drought and poor soil. Improving photosynthetic efficiency and increasing plant growth and yield under drought conditions have become important research content for researchers worldwide. Our study goal is to provide comprehensive information on photosynthesis and some candidate genes breeding of yellow horn under drought stress. In this study, seedlings' stomatal conductance, chlorophyll content, and fluorescence parameters decreased under drought stress, but non-photochemical quenching increased. The leaf microstructure showed that stomata underwent a process from opening to closing, guard cells from complete to dry, and surrounding leaf cells from smooth to severe shrinkage. The chloroplast ultrastructure showed that the changes of starch granules were different under different drought stress, while plastoglobules increased and expanded continuously. In addition, we found some differentially expressed genes related to photosystem, electron transport component, oxidative phosphate ATPase, stomatal closure, and chloroplast ultrastructure. These results laid a foundation for further genetic improvement and deficit resistance breeding of yellow horn under drought stress.

Identification of Functional Brassinosteroid Receptor Genes in Oaks and Functional Analysis of QmBRI1

Brassinosteroids (BRs) play important regulatory roles in plant growth and development, with functional BR receptors being crucial for BR recognition or signaling. Although functional BR receptors have been extensively studied in herbaceous plants, they remain largely under-studied in forest tree species. In this study, nine BR receptors were identified in three representative oak species, of which BRI1s and BRL1s were functional BR receptors. Dispersed duplications were a driving force for oak BR receptor expansion, among which the Brassinosteroid-Insensitive-1 (BRI1)-type genes diverged evolutionarily from most rosids. In oak BRI1s, we identified that methionine in the conserved Asn-Gly-Ser-Met (NGSM) motif was replaced by isoleucine and that the amino acid mutation occurred after the divergence of Quercus and Fagus. Compared with QmBRL1, QmBRI1 was relatively highly expressed during BR-induced xylem differentiation and in young leaves, shoots, and the phloem and xylem of young stems of Quercus mongolica. Based on Arabidopsis complementation experiments, we proved the important role of QmBRI1 in oak growth and development, especially in vascular patterning and xylem differentiation. These findings serve as an important supplement to the findings of the structural, functional and evolutionary studies on functional BR receptors in woody plants and provide the first example of natural mutation occurring in the conserved BR-binding region (NGSM motif) of angiosperm BRI1s.

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.

Characterization of the horse chestnut genome reveals the evolution of aescin and aesculin biosynthesis

Horse chestnut (Aesculus chinensis) is an important medicinal tree that contains various bioactive compounds, such as aescin, barrigenol-type triterpenoid saponins (BAT), and aesculin, a glycosylated coumarin. Herein, we report a 470.02 Mb genome assembly and characterize an Aesculus-specific whole-genome duplication event, which leads to the formation and duplication of two triterpenoid biosynthesis-related gene clusters (BGCs). We also show that AcOCS6, AcCYP716A278, AcCYP716A275, and AcCSL1 genes within these two BGCs along with a seed-specific expressed AcBAHD6 are responsible for the formation of aescin. Furthermore, we identify seven Aesculus-originated coumarin glycoside biosynthetic genes and achieve the de novo synthesis of aesculin in E. coli. Collinearity analysis shows that the collinear BGC segments can be traced back to early-diverging angiosperms, and the essential gene-encoding enzymes necessary for BAT biosynthesis are recruited before the splitting of Aesculus, Acer, and Xanthoceras. These findings provide insight on the evolution of gene clusters associated with medicinal tree metabolites.

Genome-wide analysis of the heat shock transcription factor family reveals saline-alkali stress responses in Xanthoceras sorbifolium

The heat shock transcription factor (HSF) family is involved in regulating growth, development, and abiotic stress. The characteristics and biological functions of HSF family member in X. sorbifolium, an important oil and ornamental plant, have never been reported. In this study, 21 XsHSF genes were identified from the genome of X. sorbifolium and named XsHSF1-XsHSF21 based on their chromosomal positions. Those genes were divided into three groups, A, B, and C, containing 12, one, and eight genes, respectively. Among them, 20 XsHSF genes are located on 11 chromosomes. Protein structure analysis suggested that XsHSF proteins were conserved, displaying typical DNA binding domains (DBD) and oligomerization domains (OD). Moreover, HSF proteins within the same group contain specific motifs, such as motif 5 in the HSFC group. All XsHSF genes have one intron in the CDS region, except XsHSF1 which has two introns. Promoter analysis revealed that in addition to defense and stress responsiveness elements, some promoters also contained a MYB binding site and elements involved in multiple hormones responsiveness and anaerobic induction. Duplication analysis revealed that XsHSF1 and XsHSF4 genes were segmentally duplicated while XsHSF2, XsHSF9, and XsHSF13 genes might have arisen from transposition. Expression pattern analysis of leaves and roots following salt-alkali treatment using qRT-PCR indicated that five XsHSF genes were upregulated and one XsHSF gene was downregulated in leaves upon NaCl treatment suggesting these genes may play important roles in salt response. Additionally, the expression levels of most XsHSFs were decreased in leaves and roots following alkali-induced stress, indicating that those XsHSFs may function as negative regulators in alkali tolerance. MicroRNA target site prediction indicated that 16 of the XsHSF genes may be regulated by multiple microRNAs, for example XsHSF2 might be regulated by miR156, miR394, miR395, miR408, miR7129, and miR854. And miR164 may effect the mRNA levels of XsHSF3 and XsHSF17, XsHSF9 gene may be regulated by miR172. The expression trends of miR172 and miR164 in leaves and roots on salt treatments were opposite to the expression trend of XsHSF9 and XsHSF3 genes, respectively. Promoter analysis showed that XsHSFs might be involved in light and hormone responses, plant development, as well as abiotic stress responses. Our results thus provide an overview of the HSF family in X. sorbifolium and lay a foundation for future functional studies to reveal its roles in saline-alkali response.

High-quality genome assembly and comparative genomic profiling of yellowhorn (Xanthoceras sorbifolia) revealed environmental adaptation footprints and seed oil contents variations

Yellowhorn (Xanthoceras sorbifolia) is a species of deciduous tree that is native to Northern and Central China, including Loess Plateau. The yellowhorn tree is a hardy plant, tolerating a wide range of growing conditions, and is often grown for ornamental purposes in parks, gardens, and other landscaped areas. The seeds of yellowhorn are edible and contain rich oil and fatty acid contents, making it an ideal plant for oil production. However, the mechanism of its ability to adapt to extreme environments and the genetic basis of oil synthesis remains to be elucidated. In this study, we reported a high-quality and near gap-less yellowhorn genome assembly, containing the highest genome continuity with a contig N50 of 32.5 Mb. Comparative genomics analysis showed that 1,237 and 231 gene families under expansion and the yellowhorn-specific gene family NB-ARC were enriched in photosynthesis and root cap development, which may contribute to the environmental adaption and abiotic stress resistance of yellowhorn. A 3-ketoacyl-CoA thiolase (KAT) gene (Xso_LG02_00600) was identified under positive selection, which may be associated with variations of seed oil content among different yellowhorn cultivars. This study provided insights into environmental adaptation and seed oil content variations of yellowhorn to accelerate its genetic improvement.

Yellowhorn Xso-miR5149-XsGTL1 enhances water-use efficiency and drought tolerance by regulating leaf morphology and stomatal density

Yellowhorn (Xanthoceras sorbifolium) is a unique edible woody oil tree species in China. Drought stress is the major yield-limiting factor of yellowhorn. MicroRNAs play an important role in regulating the response of woody plants to drought stress. However, the regulatory function of miRNAs in yellowhorn remains unclear. Here, we first constructed coregulatory networks integrated with miRNAs and their target genes. According to GO function and expression pattern analysis, we selected the Xso-miR5149-XsGTL1 module for further study. Xso-miR5149 is a key regulator of leaf morphology and stomatal density by directly mediating the expression of the transcription factor XsGTL1. Downregulation of XsGTL1 in yellowhorn led to increased leaf area and reduced stomatal density. RNA-seq analysis indicated that downregulation of XsGTL1 increased the expression of genes involved in the negative control of stomatal density, leaf morphology, and drought tolerance. After drought stress treatments, the XsGTL1-RNAi yellowhorn plants were less damaged and had higher water-use efficiency than the WT plants, while destruction of Xso-miR5149 or overexpression of XsGTL1 had the opposite effect. Our findings indicated that the Xso-miR5149-XsGTL1 regulatory module plays a critical role in controlling leaf morphology and stomatal density; hence, it's a potential candidate module for engineering enhanced drought tolerance in yellowhorn.

The double flower variant of yellowhorn is due to a LINE1 transposon-mediated insertion

As essential organs of reproduction in angiosperms, flowers, and the genetic mechanisms of their development have been well characterized in many plant species but not in the woody tree yellowhorn (Xanthoceras sorbifolium). Here, we focused on the double flower phenotype in yellowhorn, which has high ornamental value. We found a candidate C-class gene, AGAMOUS1 (XsAG1), through bovine serum albumin sequencing and genetics analysis with a Long Interpersed Nuclear Elements 1 (LINE1) transposable element fragment (Xsag1-LINE1-1) inserted into its second intron that caused a loss-of-C-function and therefore the double flower phenotype. In situ hybridization of XsAG1 and analysis of the expression levels of other ABC genes were used to identify differences between single- and double-flower development processes. These findings enrich our understanding of double flower formation in yellowhorn and provide evidence that transposon insertions into genes can reshape plant traits in forest trees.

The core triacylglycerol toolbox in woody oil plants reveals targets for oil production bioengineering

Woody oil plants are the most productive oil-bearing species that produce seeds with high levels of valuable triacylglycerols (TAGs). TAGs and their derivatives are the raw materials for many macromolecular bio-based products, such as nylon precursors, and biomass-based diesel. Here, we identified 280 genes encoding seven distinct classes of enzymes (i.e., G3PAT, LPAAT, PAP, DGAT, PDCT, PDAT, and CPT) involved in TAGs-biosynthesis. Several multigene families are expanded by large-scale duplication events, such as G3PATs, and PAPs. RNA-seq was used to survey the expression profiles of these TAG pathway-related genes in different tissues or development, indicating functional redundancy for some duplicated genes originated from the large-scale duplication events, and neo-functionalization or sub-functionalization for some of them. Sixty-two genes showed strong, preferential expression during the period of rapid seed lipid synthesis, suggesting that their might represented the core TAG-toolbox. We also revealed for the first time that there is no PDCT pathway in Vernicia fordii and Xanthoceras sorbifolium. The identification of key genes involved in lipid biosynthesis will be the foundation to plan strategies to develop woody oil plant varieties with enhanced processing properties and high oil content.

Identification and Characterization of AP2/ERF Transcription Factors in Yellow Horn

The AP2/ERF gene family involves numerous plant processes, including growth, development, metabolism, and various plant stress responses. However, several studies have been conducted on the AP2/ERF gene family in yellow horn, a new type of oil woody crop and an essential oil crop in China. According to sequence alignment and phylogenetic analyses, one hundred and forty-five AP2/ERF genes were detected from the yellow horn genome. They were divided into four relatively conserved subfamilies, including 21 AP2 genes, 119 ERBP genes, 4 RAV genes, and 1 Soloist gene. Gene analysis of XsAP2/ERF TFs showed 87 XsAP2/ERF TFs lacked introns. There were 75 pairs of collinearity relationships between X. sorbifolium and Arabidopsis, indicating a close similarity. In addition, the expression patterns of XsAP2/ERF TFs under cold treatments confirmed that the XsAP2/ERF TFs play essential roles in abiotic stress response. The expression of eight XsAP2/ERF transcription factors was verified in different tissues and under various stress treatments using RT-qPCR. This study establishes a starting point for further research to explore the potential mechanisms of identifying candidate AP2/ERF TFs that could respond to the abiotic stress of yellow horn.

A naturally-occurring phenomenon of flower color change during flower development in Xanthoceras sorbifolium

Xanthoceras sorbifolium (yellowhorn) is originated in China and is a unique tree in northern China. Yellowhorn is very popular because of the gradual color change of its flower at different flower developmental stages. During flowering development, the color at the base of yellowhorn flower petals gradually changes from yellow to purple. The mechanism of this miraculous phenomenon is still unclear. Here we show that anthocyanin accumulation during flowering development is the main reason for this color change. RT-PCR results show that the expression level of a variety of anthocyanin biosynthesis genes changes in different flower developmental stages. Realtime results show that the expression changes of these anthocyanin biosynthesis genes are positively regulated by a cluster of R2R3-MYB transcription factor genes, XsMYB113s. Furthermore, the DNA methylation analysis showed that CHH methylation status on the transposon element near the XsMYB113-1 influence its transcript level during flowering development. Our results suggest that dynamic epigenetic regulation of the XsMYB113-1 leads to the accumulation of anthocyanins during yellowhorn flower color change. These findings reemphasize the role of epigenetic regulation in flower development and provide a foundation for further studies of epigenetic regulation in long-lived woody perennials.

Genome sequencing and comparative analysis of Ficus benghalensis and Ficus religiosa species reveal evolutionary mechanisms of longevity

Ficus benghalensis and Ficus religiosa are large woody trees well known for their long lifespan, ecological and traditional significance, and medicinal properties. To understand the genomic and evolutionary aspects of these characteristics, the whole genomes of these Ficus species were sequenced using 10x Genomics linked reads and Oxford Nanopore long reads. The draft genomes of F. benghalensis and F. religiosa comprised of 392.89 Mbp and 332.97 Mbp, respectively. We established the genome-wide phylogenetic positions of the two Ficus species with respect to 50 other Angiosperm species. Comparative evolutionary analyses with other phylogenetically closer Eudicot species revealed adaptive evolution in genes involved in key cellular mechanisms associated with prolonged survival including phytohormones signaling, senescence, disease resistance, and abiotic stress tolerance, which provide genomic insights into the mechanisms conferring longevity and suggest that longevity is a multifaceted phenomenon. This study also provides clues on the existence of CAM pathway in these Ficus species.

•

First whole genome assemblies of Ficus benghalensis and Ficus religiosa trees

•

Genome-wide phylogeny with 50 other Angiosperm species

•

Evolution of genes in phytohormone signaling, senescence, and stress tolerance

•

Genomic insights into longevity-regulating mechanisms

Genomic and transcriptomic analyses provide insights into valuable fatty acid biosynthesis and environmental adaptation of yellowhorn

Yellowhorn (Xanthoceras sorbifolium) is an oil-bearing tree species growing naturally in poor soil. The kernel of yellowhorn contains valuable fatty acids like nervonic acid. However, the genetic basis underlying the biosynthesis of valued fatty acids and adaptation to harsh environments is mainly unexplored in yellowhorn. Here, we presented a haplotype-resolved chromosome-scale genome assembly of yellowhorn with the size of 490.44 Mb containing scaffold N50 of 34.27 Mb. Comparative genomics, in combination with transcriptome profiling analyses, showed that expansion of gene families like long-chain acyl-CoA synthetase and ankyrins contribute to yellowhorn fatty acid biosynthesis and defense against abiotic stresses, respectively. By integrating genomic and transcriptomic data of yellowhorn, we found that the transcription of 3-ketoacyl-CoA synthase gene XS04G00959 was consistent with the accumulation of nervonic and erucic acid biosynthesis, suggesting its critical regulatory roles in their biosynthesis. Collectively, these results enhance our understanding of the genetic basis underlying the biosynthesis of valuable fatty acids and adaptation to harsh environments in yellowhorn and provide foundations for its genetic improvement.

The yellowhorn AGL transcription factor gene XsAGL22 contributes to ABA biosynthesis and drought tolerance in poplar

Regulation of abscisic acid (ABA) biosynthesis helps plants adapt to drought stress, but the underlying molecular mechanisms are largely unclear. Here, a drought-induced transcription factor XsAGL22 was isolated from yellowhorn (Xanthoceras sorbifolium Bunge). Yeast one-hybrid and electrophoretic mobility shift assays indicated that XsAGL22 can physically bind to the promoters of the ABA biosynthesis-related genes XsNCED6 and XsBG1, and a dual-luciferase assay showed that XsAGL22 activates the promoters of the later two genes. Transient overexpression of XsAGL22 in yellowhorn leaves also increased the expression of XsNCED6 and XsBG1 and increased cellular ABA levels. Finally, heterologous overexpression of XsAGL22 in poplar increased ABA content, reduced stomatal aperture and increased drought resistance. Our results suggest that XsAGL22 is a powerful regulator of ABA biosynthesis and plays a critical role in drought resistance in plants.

Editing of the starch branching enzyme gene SBE2 generates high-amylose storage roots in cassava

The production of high-amylose cassava through CRISPR/Cas9-mediated mutagenesis of the starch branching enzyme gene SBE2 was firstly achieved. High-amylose cassava (Manihot esculenta Crantz) is desirable for starch industrial applications and production of healthier processed food for human consumption. In this study, we report the production of high-amylose cassava through CRISPR/Cas9-mediated mutagenesis of the starch branching enzyme 2 (SBE2). Mutations in two targeted exons of SBE2 were identified in all regenerated plants; these mutations, which included nucleotide insertions, and short or long deletions in the SBE2 gene, were classified into eight mutant lines. Three mutants, M6, M7 and M8, with long fragment deletions in the second exon of SBE2 showed no accumulation of SBE2 protein. After harvest from the field, significantly higher amylose (up to 56% in apparent amylose content) and resistant starch (up to 35%) was observed in these mutants compared with the wild type, leading to darker blue coloration of starch granules after quick iodine staining and altered starch viscosity with a higher pasting temperature and peak time. Further ¹H-NMR analysis revealed a significant reduction in the degree of starch branching, together with fewer short chains (degree of polymerization [DP] 15-25) and more long chains (DP>25 and especially DP>40) of amylopectin, which indicates that cassava SBE2 catalyzes short chain formation during amylopectin biosynthesis. Transition from A- to B-type crystallinity was also detected in the starches. Our study showed that CRISPR/Cas9-mediated mutagenesis of starch biosynthetic genes in cassava is an effective approach for generating novel varieties with valuable starch properties for food and industrial applications.

Population Genetics and Development of a Core Collection from Elite Germplasms of Xanthoceras sorbifolium Based on Genome-Wide SNPs

Xanthoceras sorbifolium is one of the most important species of woody oil. In this study, whole genome re-sequencing of 119 X. sorbifolium germplasms was conducted and, after filtering, 105,685,557 high-quality SNPs were identified, which were used to perform population genetics and core collection development analyses. The results from the phylogenetic, population structure, and principal component analyses showed a high level of agreement, with 119 germplasms being classified into three main groups. The germplasms were not completely classified based on their geographical origins and flower colors; furthermore, the genetic backgrounds of these germplasms were complex and diverse. The average polymorphsim information content (PIC) values for the three inferred groups clustered by structure analysis and the six classified color groups were 0.2445 and 0.2628, respectively, indicating a low to medium informative degree of genetic diversity. Moreover, a core collection containing 29.4% (35) out of the 119 X. sorbifolium germplasms was established. Our results revealed the genetic diversity and structure of X. sorbifolium germplasms, and the development of a core collection will be useful for the efficient improvement of breeding programs and genome-wide association studies.

Two divergent haplotypes from a highly heterozygous lychee genome suggest independent domestication events for early and late-maturing cultivars

Lychee is an exotic tropical fruit with a distinct flavor. The genome of cultivar ‘Feizixiao’ was assembled into 15 pseudochromosomes, totaling ~470 Mb. High heterozygosity (2.27%) resulted in two complete haplotypic assemblies. A total of 13,517 allelic genes (42.4%) were differentially expressed in diverse tissues. Analyses of 72 resequenced lychee accessions revealed two independent domestication events. The extremely early maturing cultivars preferentially aligned to one haplotype were domesticated from a wild population in Yunnan, whereas the late-maturing cultivars that mapped mostly to the second haplotype were domesticated independently from a wild population in Hainan. Early maturing cultivars were probably developed in Guangdong via hybridization between extremely early maturing cultivar and late-maturing cultivar individuals. Variable deletions of a 3.7 kb region encompassed by a pair of CONSTANS-like genes probably regulate fruit maturation differences among lychee cultivars. These genomic resources provide insights into the natural history of lychee domestication and will accelerate the improvement of lychee and related crops.

Two divergent haplotypes from a highly heterozygous lychee genome of the cultivar ‘Feizixiao’ and resequencing of 72 lychee accessions provide insights into the genome evolution and domestication history of lychee.

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.

Advances in genomics and genome editing for breeding next generation of fruit and nut crops

Fruit tree crops are an essential part of the food production systems and are key to achieve food and nutrition security. Genetic improvement of fruit trees by conventional breeding has been slow due to the long juvenile phase. Advancements in genomics and molecular biology have paved the way for devising novel genetic improvement tools like genome editing, which can accelerate the breeding of these perennial crops to a great extent. In this article, advancements in genomics of fruit trees covering genome sequencing, transcriptome sequencing, genome editing technologies (GET), CRISPR-Cas system based genome editing, potential applications of CRISPR-Cas9 in fruit tree crops improvement, the factors influencing the CRISPR-Cas editing efficiency and the challenges for CRISPR-Cas9 applications in fruit tree crops improvement are reviewed. Besides, base editing, a recently emerging more precise editing system, and the future perspectives of genome editing in the improvement of fruit and nut crops are covered.

Transcriptome and physiological analyses provide insights into the leaf epicuticular wax accumulation mechanism in yellowhorn

Plantations and production of yellowhorn, one of the most important woody oil and urban greening trees widely cultivated in northern China, have gradually become limited by drought stress. The epicuticular wax layer plays a key role in the protection of yellowhorn trees from drought and other stresses. However, there is no research on the mechanism of wax loading in yellowhorn trees. In this study, we investigated the anatomical and physiological characteristics of leaves from different germplasm resources and different parts of the same tree and compared their cuticle properties. In addition, the different expression patterns of genes involved in wax accumulation were analyzed, and a coexpression network was built based on transcriptome sequencing data. Morphological and physiological comparisons found that the sun leaves from the outer part of the crown had thicker epicuticular wax, which altered the permeability and improved the drought resistance of leaves, than did shade leaves. Based on transcriptome data, a total of 3008 and 1324 differentially expressed genes (DEGs) were identified between the sun leaves and shade leaves in glossy- and non-glossy-type germplasm resources, respectively. We identified 138 DEGs involved in wax biosynthesis and transport, including structural genes (such as LACS8, ECH1, and ns-LTP) and transcription factors (such as MYB, WRKY, and bHLH transcription factor family proteins). The coexpression network showed a strong correlation between these DEGs. The differences in gene expression patterns between G- and NG-type germplasm resources under different light conditions were very clear. These results not only provide a theoretical basis for screening and developing drought-resistant yellowhorn germplasm resources but also provide a data platform to reveal the wax accumulation process of yellowhorn leaves.

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.

Quantitative proteomic analysis of Xanthoceras sorbifolium Bunge seedlings in response to drought and heat stress

Yellowhorn (Xanthoceras sorbifolium Bunge) is a woody oil species that is widely distributed in northwestern China. To investigate the molecular mechanisms underlying the drought and heat tolerance response of yellowhorn seedlings, changes in protein abundance were analyzed via comparative proteomics. Drought and heat treatment of seedlings was applied in growth chamber, and the leaves were harvested after 7 days of treatment. The total protein was extracted, and comparative proteomic analysis was performed via isobaric tag for relative and absolute quantitation (iTRAQ). The abundance of most of the proteins associated with oxidative phosphorylation, NADH dehydrogenase and superoxide dismutase (SOD) was reduced. The differential proteins associated with photosynthesis enzymes indicated that stress had different effects on photosystem I (PSI) and photosystem II (PSII). After comprehensively analyzing the results, we speculated that drought and heat stress could hinder the synthesis of riboflavin, reducing NADH dehydrogenase content, which might further have an impact on energy utilization. Yellowhorn seedlings relied on Fe-Mn SOD enzymes rather than Cu/Zn SOD enzymes to remove reactive oxygen species (ROS). In addition, heat-shock proteins (HSPs) had significant increase and played a key role in stress response, which could be divided into two categories according to their transcription and translation efficiency. Over all, the results can provide a basis for understanding the molecular mechanism underlying resistance to drought and heat stress in yellowhorn and for subsequent research of posttranslational modification-related omics of key proteins.

The changing face of genome assemblies: Guidance on achieving high-quality reference genomes

The quality of genome assemblies has improved rapidly in recent years due to continual advances in sequencing technology, assembly approaches, and quality control. In the field of molecular ecology, this has led to the development of exceptional quality genome assemblies that will be important long-term resources for broader studies into ecological, conservation, evolutionary, and population genomics of naturally occurring species. Moreover, the extent to which a single reference genome represents the diversity within a species varies: pan-genomes will become increasingly important ecological genomics resources, particularly in systems found to have considerable presence-absence variation in their functional content. Here, we highlight advances in technology that have raised the bar for genome assembly and provide guidance on standards to achieve exceptional quality reference genomes. Key recommendations include the following: (a) Genome assemblies should include long-read sequencing except in rare cases where it is effectively impossible to acquire adequately preserved samples needed for high molecular weight DNA standards. (b) At least one scaffolding approach should be included with genome assembly such as Hi-C or optical mapping. (c) Genome assemblies should be carefully evaluated, this may involve utilising short read data for genome polishing, error correction, k-mer analyses, and estimating the percent of reads that map back to an assembly. Finally, a genome assembly is most valuable if all data and methods are made publicly available and the utility of a genome for further studies is verified through examples. While these recommendations are based on current technology, we anticipate that future advances will push the field further and the molecular ecology community should continue to adopt new approaches that attain the highest quality genome assemblies.

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.

Distinct Evolutionary Patterns of NBS-Encoding Genes in Three Soapberry Family (Sapindaceae) Species

Nucleotide-binding site (NBS)-type disease resistance genes (R genes) play key roles in plant immune responses and have co-evolved with pathogens over the course of plant lifecycles. Comparative genomic studies tracing the dynamic evolution of NBS-encoding genes have been conducted using many important plant lineages. However, studies on Sapindaceae species have not been performed. In this study, a discrepant number of NBS-encoding genes were identified in the genomes of Xanthoceras sorbifolium (180), Dinnocarpus longan (568), and Acer yangbiense (252). These genes were unevenly distributed and usually clustered as tandem arrays on chromosomes, with few existed as singletons. The phylogenetic analysis revealed that NBS-encoding genes formed three monophyletic clades, RPW8-NBS-LRR (RNL), TIR-NBS-LRR (TNL), and CC-NBS-LRR (CNL), which were distinguished by amino acid motifs. The NBS-encoding genes of the X. sorbifolium, D. longan, and A. yangbiense genomes were derived from 181 ancestral genes (three RNL, 23 TNL, and 155 CNL), which exhibited dynamic and distinct evolutionary patterns due to independent gene duplication/loss events. Specifically, X. sorbifolium exhibited a “first expansion and then contraction” evolutionary pattern, while A. yangbiense and D. longan exhibited a “first expansion followed by contraction and further expansion” evolutionary pattern. However, further expansion in D. longan was stronger than in A. yangbiense after divergence, suggesting that D. longan gained more genes in response to various pathogens. Additionally, the ancient and recent expansion of CNL genes generated the dominance of this subclass in terms of gene numbers, while the low copy number status of RNL genes was attributed to their conserved functions.

Integrated transcriptomic and metabolomic analyses of yellow horn (Xanthoceras sorbifolia) in response to cold stress

Xanthoceras sorbifolia, a medicinal and oil-rich woody plant, has great potential for biodiesel production. However, little study explores the link between gene expression level and metabolite accumulation of X. sorbifolia in response to cold stress. Herein, we performed both transcriptomic and metabolomic analyses of X. sorbifolia seedlings to investigate the regulatory mechanism of resistance to low temperature (4 °C) based on physiological profile analyses. Cold stress resulted in a significant increase in the malondialdehyde content, electrolyte leakage and activity of antioxidant enzymes. A total of 1,527 common differentially expressed genes (DEGs) were identified, of which 895 were upregulated and 632 were downregulated. Annotation of DEGs revealed that amino acid metabolism, glycolysis/gluconeogenesis, starch and sucrose metabolism, galactose metabolism, fructose and mannose metabolism, and the citrate cycle (TCA) were strongly affected by cold stress. In addition, DEGs within the plant mitogen-activated protein kinase (MAPK) signaling pathway and TF families of ERF, WRKY, NAC, MYB, and bHLH were transcriptionally activated. Through metabolomic analysis, we found 51 significantly changed metabolites, particularly with the analysis of primary metabolites, such as sugars, amino acids, and organic acids. Moreover, there is an overlap between transcript and metabolite profiles. Association analysis between key genes and altered metabolites indicated that amino acid metabolism and sugar metabolism were enhanced. A large number of specific cold-responsive genes and metabolites highlight a comprehensive regulatory mechanism, which will contribute to a deeper understanding of the highly complex regulatory program under cold stress in X. sorbifolia.

Retracted Article: Investigation of yellow horn (Xanthoceras sorbifolia Bunge) transcriptome in response to different abiotic stresses: a comparative RNA-Seq study

Yellow horn (Xanthoceras sorbifolia Bunge) is a well-known oil-rich seed shrub which can grow well in barren and arid environments in the northern part of China. Yellow horn has received worldwide attention because of its excellent economic and environmental value. However, because of its limited genetic data, little information can be found regarding the molecular defense mechanisms of yellow horn exposed to various abiotic stresses. In view of this, the current study aims to investigate the impact of different abiotic stresses (i.e. NaCl, ABA and low temperature) on the transcriptome of yellow horn using RNA-Seq. Based on the transcriptome sequencing data, approximately 27% to 45% of stress-responsive genes were found highly expressed after stress treatment for 24 h. In addition, these genes were found to be still expressed after stress treatment for 48 h. However, many additional genes were stress-regulated after 48 h treatment compared with the 24 h treatment. GO enrichment analysis revealed that the expression patterns of the stress-responsive, type-specific terms were generally down-regulated. Most shared GO terms were primarily involved in protein folding, unfolding protein binding, protein transport and protein modification. Further, transcription factors (TFs), such as ERFs, bHLH, GRAS and NAC, were found to be enriched only in the low temperature treatment group, particularly the ERF TFs families. These combined results suggested that yellow horn may have developed specific molecular defense systems against diverse abiotic stresses.

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.