Comparative population genomics dissects the genetic basis of seven domestication traits in jujube

Haplotype-resolved T2T reference genomes for wild and domesticated accessions shed new insights into the domestication of jujube

Chinese jujube (Ziziphus jujuba Mill.) is one of the most important deciduous tree fruits in China, with substantial economic and nutritional value. Jujube was domesticated from its wild progenitor, wild jujube (Z. jujuba var. spinosa), and both have high medicinal value. Here we report the 767.81- and 759.24-Mb haplotype-resolved assemblies of a dry-eating 'Junzao' jujube (JZ) and a wild jujube accession (SZ), using a combination of multiple sequencing strategies. Each assembly yielded two complete haplotype-resolved genomes at the telomere-to-telomere (T2T) level, and ~81.60 and 69.07 Mb of structural variations were found between the two haplotypes within JZ and SZ, respectively. Comparative genomic analysis revealed a large inversion on each of chromosomes 3 and 4 between JZ and SZ, and numerous genes were affected by structural variations, some of which were associated with starch and sucrose metabolism. A large-scale population analysis of 672 accessions revealed that wild jujube originated from the lower reaches of the Yellow River and was initially domesticated at local sites. It spread widely and was then independently domesticated at the Shanxi-Shaanxi Gorge of the middle Yellow River. In addition, we identified some new selection signals regions on genomes, which are involved in the tissue development, pollination, and other aspects of jujube tree morphology and fertilization domestication. In conclusion, our study provides high-quality reference genomes of jujube and wild jujube and new insights into the domestication history of jujube.

A natural mutation in the promoter of the aconitase gene ZjACO3 influences fruit citric acid content in jujube

Jujube (Ziziphus jujuba Mill.) is the most economically important fruit tree of the Rhamnaceae and was domesticated from wild or sour jujube (Z. jujuba Mill. var. spinosa Hu). During the process of domestication, there was a substantial reduction in the content of organic acids, particularly malate and citrate, which greatly influence the taste and nutritional value of the fruit. We previously demonstrated that ZjALMT4 is crucial for malate accumulation. However, the mechanism of citrate degradation in jujube remains poorly understood. In the present study, aconitase ZjACO3 was shown to participate in citric acid degradation in the cytoplasm through the GABA pathway. Interestingly, we discovered an E-box mutation in the ZjACO3 promoter (-484A > G; CAAGTG in sour jujube mutated to CAGGTG in cultivated jujube) that was strongly correlated with fruit citrate content; 'A' represented a high-citrate genotype and 'G' represented a low-citrate genotype. We developed and validated an ACO-based Kompetitive allele-specific PCR (KASP) marker for determining citric acid content. Yeast one-hybrid screening, transient dual-luciferase assays, and overexpression analyses showed that the transcription factor ZjbHLH113 protein directly binds to CAGGTG in the promoter of ZjACO3 in cultivated jujube plants, transcriptionally activating ZjACO3 expression, and enhancing citric acid degradation. Conversely, binding ability of the ZjbHLH113 protein to CAAGTG was weakened in sour jujube, thereby promoting citrate accumulation in the fruit. These findings will assist in elucidating the mechanism by which ZjACO3 modulates citrate accumulation in sour jujube and its cultivars.

Comparative analysis of jujube and sour jujube gave insight into their difference in genetic diversity and suitable habitat

Jujube (Ziziphus jujuba var. jujuba Mill.) and sour jujube (Z. jujuba var. spinosa (Bunge) Hu ex H.F.Chow.) are economically, nutritionally, and ecologically significant members of the Rhamnaceae family. Despite their importance, insufficient research on their genetics and habitats has impeded effective conservation and utilization. To address this knowledge gap, we conducted plastome sequencing, integrated distribution data from China, and assessed genetic diversity and suitable habitat. The plastomes of both species exhibited high conservation and low genetic diversity. A new-found 23 bp species-specific Indel in the petL-petG enabled us to develop a rapid Indel-based identification marker for species discrimination. Phylogenetic analysis and dating illuminated their genetic relationship, showing speciation occurred 6.9 million years ago, in a period of dramatic global temperature fluctuations. Substantial variations in suitable climatic conditions were observed, with the mean temperature of the coldest quarter as the primary factor influencing distributions (-3.16°C-12.73°C for jujube and -5.79°C to 4.11°C for sour jujube, suitability exceeding 0.6). Consequently, distinct conservation strategies are warranted due to differences in suitable habitats, with jujube having a broader distribution and sour jujube concentrated in Northern China. In conclusion, disparate habitats and climatic factors necessitate tailored conservation approaches. Comparing genetic diversity and developing rapid species-specific primers will further enhance the sustainable utilization of these valuable species.

ZjWRKY23 and ZjWRKY40 Promote Fruit Size Enlargement by Targeting and Downregulating Cytokinin Oxidase/Dehydrogenase 5 Expression in Chinese Jujube

Fruit size is crucial for fruit trees, as it contributes to both quality and yield. However, the underlying mechanism of fruit size regulation remains largely unknown. Taking advantage of using a fruit double-sized bud mutant of Chinese jujube, "Jinkuiwang" and its wild type, "Jinsixiaozao", we carried out a comprehensive study on the mechanism of fruit size development in jujube. Using weighted gene coexpression network analyses, a number of candidate regulators for fruit size including those involved in hormonal signaling pathways, transcription factors, and heat shock proteins were identified. A hub gene named cytokinin oxidase/dehydrogenase 5 (ZjCKX5), responsible for cytokinin degradation, was found to play a negative role in regulating fruit size development, and overexpressing ZjCKX5 in tomato and Arabidopsis resulted in much smaller fruits and dwarf plants. Furthermore, another two hub genes, ZjWRKY23 and ZjWRKY40 transcription factors, were found to participate in fruit size regulation by targeting and downregulating the ZjCKX5 expression. Overexpressing ZjWRKY23 or ZjWRKY40 in tomato led to much larger fruits and promoted plant architecture. Based on these results, a molecular framework for jujube fruit size regulation, namely, ZjWRKY-ZjCKX5 module, was proposed. This study provides a new insight into the molecular networks underlying fruit size regulation.

Large-scale population structure and genetic architecture of agronomic traits of garlic

Garlic, an asexually propagated crop, is the second important bulb crop after the onion and is used as a vegetable and medicinal plant. Abundant and diverse garlic resources have been formed over thousands of years of cultivation. However, genome variation, population structure and genetic architecture of garlic agronomic traits were still not well elucidated. Here, 1 100 258 single nucleotide polymorphisms (SNPs) were identified using genotyping-by-sequencing in 606 garlic accessions collected from 43 countries. Population structure, principal component and phylogenetic analysis showed that these accessions were divided into five subpopulations. Twenty agronomic traits, including above-ground growth traits, bulb-related and bolt-related traits in two consecutive years were implemented in a genome-wide association study. In total, 542 SNPs were associated with these agronomic traits, among which 188 SNPs were repeatedly associated with more than two traits. One SNP (chr6: 1896135972) was repeatedly associated with ten traits. These associated SNPs were located within or near 858 genes, 56 of which were transcription factors. Interestingly, one non-synonymous SNP (Chr4: 166524085) in ribosomal protein S5 was repeatedly associated with above-ground growth and bulb-related traits. Additionally, gene ontology enrichment analysis of candidate genes for genomic selection regions between complete-bolting and non-bolting accessions showed that these genes were significantly enriched in 'vegetative to reproductive phase transition of meristem', 'shoot system development', 'reproductive process', etc. These results provide valuable information for the reliable and efficient selection of candidate genes to achieve garlic genetic improvement and superior varieties.

A comparative plastomic analysis of Ziziphus jujuba var. spinosa (Bunge) Hu ex H. F. Chow and implication of the origin of Chinese jujube

Comparative plastomics can be used to explicitly dissect various types of plastome variation. In the present study, the plastome variation pattern of Ziziphus jujuba var. spinosa (also called sour jujube) and its phylogenomic relationship with Chinese jujube were investigated. Plastomes of 21 sour jujube individuals were sequenced and assembled. The length of the sour jujube plastomes ranged between 159399 and 161279 bp. The plastomes exhibited collinearity of structure, gene order and content. The most divergent regions were located in the intergenic spacers, such as trnR-UCU-atpA and psbZ-trnG-UCC. Sliding window analysis demonstrated that the sequence variation among the sour jujube plastomes was relatively low. Sixty-two to 76 SSRs with 4 motif types were identified in the sour jujube plastomes with a predominant motif type of A/T. Three protein-coding genes exhibited higher nonsynonymous/synonymous substitution ratios, indicating that these genes may undergo positive selection. A total of 80 SNPs were detected and 1266 potential RNA editing sites of 23 protein-coding genes were predicted. In the phylogenomic tree constructed, sour jujube has a sister relationship to Chinese jujube, which indicates that Chinese jujube may have originated or been domesticated from sour jujube. The present study explicitly investigated the individual-level plastome variation of sour jujube and provides potential valuable molecular markers for future genetic-related study of this lineage.

Genome-wide association study of 23 flowering phenology traits and 4 floral agronomic traits in tree peony (Paeonia section Moutan DC.) reveals five genes known to regulate flowering time

Tree peony is a unique traditional flower in China, with large, fragrant, and colorful flowers. However, a relatively short and concentrated flowering period limits the applications and production of tree peony. A genome-wide association study (GWAS) was conducted to accelerate molecular breeding for the improvement of flowering phenology traits and ornamental phenotypes in tree peony. A diverse panel of 451 tree peony accessions was phenotyped for 23 flowering phenology traits and 4 floral agronomic traits over 3 years. Genotyping by sequencing (GBS) was used to obtain a large number of genome-wide single-nucleotide polymorphisms (SNPs) (107 050) for the panel genotypes, and 1047 candidate genes were identified by association mapping. Eighty-two related genes were observed during at least 2 years for flowering, and seven SNPs repeatedly identified for multiple flowering phenology traits over multiple years were highly significantly associated with five genes known to regulate flowering time. We validated the temporal expression profiles of these candidate genes and highlighted their possible roles in the regulation of flower bud differentiation and flowering time in tree peony. This study shows that GWAS based on GBS can be used to identify the genetic determinants of complex traits in tree peony. The results expand our understanding of flowering time control in perennial woody plants. Identification of markers closely related to these flowering phenology traits can be used in tree peony breeding programs for important agronomic traits.

Integrated microRNA and transcriptome profiling reveals the regulatory network of embryo abortion in jujube

Hybridization is an important approach to the production of new varieties with exceptional traits. Although the kernel rate of wild jujube (Ziziphus jujuba Mill. var. spinosa Hu.) is generally high, that of cultivated jujube (Z. jujuba Mill.) is low, greatly hampering the jujube breeding process. However, the mechanism by which this trait changed during jujube domestication remains unclear. Here, we explored the potential regulatory network that governs jujube embryo abortion using correlation analysis of population traits, artificial pollination, sugar content measurements and multi-omics analysis. The results showed that embryo abortion was an important reason for the low kernel rate of cultivated jujube, and kernel rate was negatively correlated with edible rate. Twenty-one days after pollination was a critical period for embryo abortion. At this time, the sugar content of cultivated 'Junzao' kernels decreased significantly compared with that of the pulp, but sugar content remained relatively stable in kernels of wild 'Suanzao'. A total of 1142 differentially expressed genes targeted by 93 microRNAs (miRNAs) were identified by transcriptome, miRNA and degradome sequencing, and may be involved in the regulation of embryo abortion during kernel development. Among them, DELLA protein, TCP14 and bHLH93 transcription factors have been shown to participate in the regulation of embryonic development. Our findings suggest that carbohydrate flow between different tissues of cultivated jujube exhibits a bias toward the pulp at 21 days after pollination, thereby restricting the process of kernel development. This information enhances our understanding of the embryo abortion process and reveals miRNA-target gene pairs that may be useful for molecular-assisted breeding.

Physico-chemical properties, antioxidant activity, and ACE inhibitory activity of protein hydrolysates from wild jujube seed

Wild jujube seed protein (WJSP) as one kind of functional food material has attracted much attention due to its highly nutritive and medicinal value in anti-inflammatory and improving immunomodulatory ability. However, owing to its large molecular weight and complex structure, biological activities of WJSP were greatly limited and cannot be fully utilized by the human body. Therefore, how to improve the bioavailability of WJSP and develop promising WJSP nutritious materials is a great challenge. In this work, wild jujube seed protein hydrolysates (WJSPHs) were prepared from WJSP via enzymatic hydrolysis method, and their physico-chemical properties, antioxidant activity, and angiotensin converting enzyme (ACE) inhibitory activity in vitro have been investigated for the first time. SDS-PAGE electrophoresis and size-exclusion chromatographic results indicate that WJSPHs have lower molecular weight distribution (< 5,000 Da) than WJSP. Circular dichroism (CD) spectroscopy and Fourier transform infrared spectroscopy (FTIR) results illustrated that random coil is the main secondary structure of WJSPHs. Antioxidant experiments indicate that WJSPHs exhibit high radicals-scavenging ability of 1,1-diphenyl-2-picrylhydrazyl (DPPH) radicals (94.60%), 2,2'-azinobis-(3-ethylbenzthiazoline-6-sulfonate) (ABTS⁺ ) radicals (90.84%), superoxide radicals (44.77%), and hydroxyl radicals (47.77%). In vitro, WJSPHs can significantly decrease the accumulation of reactive oxygen species (ROS) and malondialdehyde (MDA), and increase the activity of superoxide dismutase (SOD) and glutathione peroxidase (GSH-Px) in HepG2 cells. Moreover, ACE activity was found that can be significantly inhibited by WJSPHs (73.02%). Therefore, all previously mentioned results suggest that WJSPHs may be a promising antioxidant food to prevent oxidative-related diseases in future. PRACTICAL APPLICATION: This study shows that WJSPHs exhibit high antioxidant activity and ACE inhibitory activity in vitro, which provide potential application value as antioxidant peptides to prevent oxidative-related diseases.

Distribution of polyphenolic compounds, antioxidant potential, and free amino acids in Ziziphus fruits extract; a study for determining the influence of wider geography

Ziziphus fruits have attracted much attention within the field of medicine due to their high potential against central nervous system disorders. Abundance of secondary metabolites and their composition is key to the pharmaceutical potential and commercial qualities of plants. The in vitro antioxidant activities of Ziziphus nummularia (Burm. f.) and Ziziphus oxyphylla Edgew fruit extract were analyzed using 2,2‐diphenil‐1‐pycrilhydrazyl (DPPH) and 2,2′‐azino‐bis (3‐ethylbenzothiazoline)‐6‐sulfonic acid (ABTS) free radical scavenging assay methods. Phenolic profiles were explored using high‐performance liquid chromatography‐diode array detector (HPLC‐DAD). The result revealed high concentration of polyphenols and their antioxidant potential. In Z. nummularia, the total phenolic content (TPC) (80.270 ± 0.422 μg/ml), DPPH (62.03 ± 0.98 μg/ml), ABTS (66.32 ± 0.73 μg/ml), and TFC (90.683 ± 0.274 μg/ml) were recorded. However, in Z. oxyphylla, DPPH and ABTS values were 60.66 ± 0.56 μg/ml and 61.55 ± 0.77 μg/ml, respectively, indicative of the impacts of climate and soil nutrients. The overall screening of phytochemicals revealed that both the Ziziphus species contain diverse bioactive compounds, including spinacetine‐3‐O‐(2 feruloyl glucopyranosyl)‐glucopyranoside, kaempferol‐3‐O‐glucoside‐7‐O‐glucoside, and caffeic acid; p‐hydroxybenzoyl hexose, p‐coumaric acid, salicylic acid, and ellagic acid pentoxide. Additionally, the highest concentrated amino acid noted was of Lue 0.19 g/100 g with 596.00 retention time (RT), followed by Thr>Ale>Isl>Phya>Val in Z. nummularia. Similarly, the highest concentration of Lue amino acid was recorded as 0.18/100 g with 564.52 RT followed by Pr>Thr>Ale>Lue>Isl>Phya>Val in all genotypes of Z. oxyphylla. Reporting of polyphenols rich and stable species along with identification of favorable regions of cultivation for amino acid, polyphenols, and higher antioxidant potential may lead the way for the identification of elite clones of the species as well as may result in new drug discovery.

Jujube metabolome selection determined the edible properties acquired during domestication

Plants supply both food and medicinal compounds, which are ascribed to diverse metabolites produced by plants. However, studies on domestication-driven changes in the metabolome and genetic basis of bioactive molecules in perennial fruit trees are generally lacking. Here, we conducted multidimensional analyses revealing a singular domestication event involving the genomic and metabolomic selection of jujube trees (Ziziphus jujuba Mill.). The genomic selection for domesticated genes was highly enriched in metabolic pathways, including carbohydrates and specialized metabolism. Domesticated metabolome profiling indicated that 187 metabolites exhibited significant divergence as a result of directional selection. Malic acid was directly selected during domestication, and the simultaneous selection of specialized metabolites, including triterpenes, consequently lead to edible properties. Cyclopeptide alkaloids (CPAs) were specifically targeted for the divergence between dry and fresh cultivars. We identified 1080 significantly associated loci for 986 metabolites. Among them, 15 triterpenes were directly selected at six major loci, allowing the identification of a homologous cluster containing seven 2,3-oxidosqualene cyclases (OSCs). An OSC gene was found to contribute to the reduction in the content of triterpenes during domestication. The complete pathway for synthesizing ursolic acid was dissected by integration of the metabolome and transcriptome. Additionally, an N-methyltransferase involved in the biosynthesis of CPA and responsible for inter-cultivar content variation was identified. The present study promotes our understanding of the selection process of the global metabolome subsequent to fruit tree domestication and facilitates the genetic manipulation of specialized metabolites to enhance their edible traits.

Chromosome-Scale Genome Assembly for Chinese Sour Jujube and Insights Into Its Genome Evolution and Domestication Signature

Sour or wild jujube fruits and dried seeds are popular food all over the world. In this study, we reported a high-quality genome assembly of sour jujube (Ziziphus jujuba Mill. var. spinosa), with a size of 406 Mbp and scaffold N50 of 30.3 Mbp, which experienced only γ hexaploidization event, without recent genome duplication. Population structure analysis identified four jujube subgroups (two domesticated ones, i.e., D1 in West China and D2 in East/SouthEast China, semi-wild, and wild), which underwent an evolutionary history of a significant decline of effective population size during the Last Glacial Period. The respective selection signatures of three subgroups were discovered, such as strong peaks on chromosomes #3 in D1, #1 in D2, and #4 in wild. Genes under the most significant selection on chromosomes #4 in wild were confirmed to be involved in fruit variations among jujube accessions, in transcriptomic analysis. Our study offered novel insights into the jujube population structure and domestication and provided valuable genomic resources for jujube improvement in stress response and fruit flavor in the future.

Jujube witches' broom phytoplasma effectors SJP1 and SJP2 induce lateral bud outgrowth by repressing the ZjBRC1-controlled auxin efflux channel

Comprehensively controlling phytoplasma-associated jujube witches' broom (JWB) disease is extremely challenging for the jujube industry. Although the pathogenesis of phytoplasma disease has been highlighted in many plant species, the release of lateral buds from dormancy under JWB phytoplasma infection has not been characterized in woody perennial jujube. Here, two 16SrV-B group phytoplasma effectors, SJP1 and SJP2, were experimentally determined to induce witches' broom with increased lateral branches. In vivo interaction and subcellular localization analyses showed that both SJP1 and SJP2 were translocated from the cytoplasm to the nucleus to target the CYC/TB1-TCP transcription factor ZjBRC1. The N- and C-terminal coiled-coil domains of SJP1 and SJP2 were required for the TCP-binding ability. ZjBRC1 bound directly to the auxin efflux carrier ZjPIN1c/3 promoters and down-regulated their expression to promote the accumulation of endogenous auxin indole-3-acetic acid in jujube calli. Furthermore, JWB phytoplasma infection suppressed ZjBRC1 accumulation and induced ZjPIN1c/3 expression to stimulate lateral bud outgrowth. Therefore, SJP1 and SJP2 stimulate lateral bud outgrowth, at least partly, by repressing the ZjBRC1-controlled auxin efflux channel in jujube, representing a potential strategy for comprehensive phytoplasma-associated disease control and a resource for gene editing breeding to create new cultivars with varying degrees of shoot branching.

Resequencing and transcriptomic analysis reveal differences in nitrite reductase in jujube fruit (Ziziphus jujuba Mill.)

BACKGROUND

Jujube is a typical fruit tree species from China. 'Muzao', a cracking-susceptible cultivar, and 'Linhuang No. 1', a cracking-resistant cultivar, were selected in a previous study as contrasting research materials. Whole-genome resequencing and transcriptomic analysis of 'Linhuang No. 1' and 'Muzao' allowed the identification of differentially expressed genes with different gene structures between the two cultivars and could be helpful in explaining the differences and similarities between the two cultivars.

RESULTS

Resequencing identified 664,129 polymorphic variable sites between 'Linhuang No. 1' and 'Muzao'. To determine the genetic relationship among 'Linhuang No. 1', 'Muzao' and the jujube genome reference cultivar 'Dongzao', the characteristic polymorphic variable sites were analysed by principal component analysis. The genetic relationship between 'Linhuang No. 1' and 'Muzao' was closer than that of either variety and 'Dongzao'. Nineteen differentially expressed genes were identified by combining transcriptomic analysis with resequencing analysis. LOC107427052 (encoding a nitrite reductase) was identified by Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analysis for further study. The identified insertion was not in the domain region of the LOC107427052 gene coding sequence (CDS) region and was verified by the finding that the insertion did not affect translation of the protein. The LOC107427052 gene expression levels, nitrite reductase activities and nitrite contents of 'Muzao' were significantly higher than the corresponding values of 'Linhuang No. 1' at the young fruit stage. There was no significant difference in the quantity of the product of nitrite reductase, namely, ammonia, between the two cultivars.

CONCLUSIONS

The present study was the first to explore the differences between different jujube cultivars ('Linhuang No. 1' and 'Muzao') by combining genome resequencing and transcriptomics. LOC107427052 (encoding a nitrite reductase) was characterized by KEGG enrichment analysis. The insertion in the CDS region of the LOC107427052 gene provides a new direction for the study of nitrogen metabolism in jujube. Our study has laid a foundation for the comparative analysis of nitrite metabolism between the jujube cultivars 'Linhuang No. 1' and 'Muzao'.

Genomic analyses of diverse wild and cultivated accessions provide insights into the evolutionary history of jujube

The Chinese jujube (Ziziphus jujuba Mill.), a member of the Rhamnaceae family, is an important perennial fruit tree crop of substantial economic, ecological and nutritional value, and is also used as a traditional herbal medicine. Here, we report the resequencing of 493 jujube accessions, including 202 wild and 291 cultivated accessions at >16× depth. Our population genomic analyses revealed that the Shanxi-Shaanxi area of China was jujube's primary domestication centre and that jujube was then disseminated into East China before finally extending into South China. Divergence events analysis indicated that Ziziphus acidojujuba and Ziziphus jujuba diverged around 2.7 Mya, suggesting the interesting possibility that a long pre-domestication period may have occurred prior to human intervention. Using the large genetic polymorphism data set, we identified a 15-bp tandem insertion in the promoter of the jujube ortholog of the POLLEN DEFECTIVE IN GUIDANCE 1 (POD1) gene, which was strongly associated with seed-setting rate. Integrating genome-wide association study (GWAS), transcriptome data, expression analysis and transgenic validation in tomato, we identified a DA3/UBIQUITIN-SPECIFIC PROTEASE 14 (UBP14) ortholog, which negatively regulate fruit weight in jujube. We also identified candidate genes, which have likely influenced the selection of fruit sweetness and crispness texture traits among fresh and dry jujubes. Our study not only illuminates the genetic basis of jujube evolution and domestication and provides a deep and rich genomic resource to facilitate both crop improvement and hypothesis-driven basic research, but also identifies multiple agriculturally important genes for this unique perennial tree fruit species.