The persimmon (Diospyros oleifera Cheng) genome provides new insights into the inheritance of astringency and ancestral evolution

The complex hexaploid oil-Camellia genome traces back its phylogenomic history and multi-omics analysis of Camellia oil biosynthesis

Oil-Camellia (Camellia oleifera), belonging to the Theaceae family Camellia, is an important woody edible oil tree species. The Camellia oil in its mature seed kernels, mainly consists of more than 90% unsaturated fatty acids, tea polyphenols, flavonoids, squalene and other active substances, which is one of the best quality edible vegetable oils in the world. However, genetic research and molecular breeding on oil-Camellia are challenging due to its complex genetic background. Here, we successfully report a chromosome-scale genome assembly for a hexaploid oil-Camellia cultivar Changlin40. This assembly contains 8.80 Gb genomic sequences with scaffold N50 of 180.0 Mb and 45 pseudochromosomes comprising 15 homologous groups with three members each, which contain 135 868 genes with an average length of 3936 bp. Referring to the diploid genome, intragenomic and intergenomic comparisons of synteny indicate homologous chromosomal similarity and changes. Moreover, comparative and evolutionary analyses reveal three rounds of whole-genome duplication (WGD) events, as well as the possible diversification of hexaploid Changlin40 with diploid occurred approximately 9.06 million years ago (MYA). Furthermore, through the combination of genomics, transcriptomics and metabolomics approaches, a complex regulatory network was constructed and allows to identify potential key structural genes (SAD, FAD2 and FAD3) and transcription factors (AP2 and C2H2) that regulate the metabolism of Camellia oil, especially for unsaturated fatty acids biosynthesis. Overall, the genomic resource generated from this study has great potential to accelerate the research for the molecular biology and genetic improvement of hexaploid oil-Camellia, as well as to understand polyploid genome evolution.

Main flavor compounds and molecular regulation mechanisms in fruits and vegetables

Fruits and vegetables (F&V) are an indispensable part of a healthy diet. The volatile and nonvolatile compounds present in F&V constitute unique flavor substances. This paper reviews the main flavor substances present in F&V, as well as the biosynthetic pathways and molecular regulation mechanisms of these compounds. A series of compounds introduced include aromatic substances, soluble sugars and organic acids, which constitute the key flavor substances of F&V. Esters, phenols, alcohols, amino acids and terpenes are the main volatile aromatic substances, and nonvolatile substances are represented by amino acids, fatty acids and carbohydrates; The combination of these ingredients is the cause of the sour, sweet, bitter, astringent and spicy taste of these foods. This provides a theoretical basis for the study of the interaction between volatile and nonvolatile substances in F&V, and also provides a research direction for the healthy development of food in the future.

Preferential transport activity of DkDTX5/MATE5 affects the formation of different astringency in persimmon

Proanthocyanidins (PAs) are specialized metabolites that influence persimmon fruit quality. Normal astringent (A)-type and non-astringent (NA)-type mutants show significant variation in PA accumulation, but the influencing mechanism remains unclear. In this study, among the six identified DTXs/MATEs proteins associated with PA accumulation, we observed that allelic variation and preferential transport by DkDTX5/MATE5 induced variation in PA accumulation for A-type and NA-type fruit. The expression pattern of DkDTX5/MATE5 was correlated with PA accumulation in NA-type fruit. Upregulation and downregulation of DkDTX5/MATE5 promoted and inhibited PA accumulation, respectively, in the NA-type fruit. Interestingly, transporter assays of Xenopus laevis oocytes indicated that DkDTX5/MATE5 preferentially transported the PA precursors catechin, epicatechin, and epicatechin gallate, resulting in their increased ratios relative to the total PAs, which was the main source of variation in PA accumulation between the A-type and NA-type. The allele lacking Ser-84 in DkDTX5/MATE5 was identified as a dominantly expressed gene in the A-type and lost its transport function. Site-directed mutagenesis revealed that DkDTX5/MATE5 binds to PA precursors via Ser-84. These findings clarify the association between the transporter function of DkDTX5/MATE5 and PA variation, and can contribute to the breeding of new cultivars with improved fruit quality.

Exogenous gibberellin delays maturation in persimmon fruit through transcriptional activators and repressors

As the harvest season of most fruit is concentrated, fruit maturation manipulation is essential for the fresh fruit industry to prolong sales time. Gibberellin (GA), an important phytohormone necessary for plant growth and development, has also shown a substantial regulatory effect on fruit maturation; however, its regulatory mechanisms remain inconclusive. In this research, preharvest GA3 treatment effectively delayed fruit maturation in several persimmon (Diospyros kaki) cultivars. Among the proteins encoded by differentially expressed genes, 2 transcriptional activators (NAC TRANSCRIPTION FACTOR DkNAC24 and ETHYLENE RESPONSIVE FACTOR DkERF38) and a repressor (MYB-LIKE TRANSCRIPTION FACTOR DkMYB22) were direct regulators of GERANYLGERANYL DIPHOSPHATE SYNTHASE DkGGPS1, LYSINE HISTIDINE TRANSPORTER DkLHT1, and FRUCTOSE-BISPHOSPHATE ALDOLASE DkFBA1, respectively, resulting in the inhibition of carotenoid synthesis, outward transport of an ethylene precursor, and consumption of fructose and glucose. Thus, the present study not only provides a practical method to prolong the persimmon fruit maturation period in various cultivars but also provides insights into the regulatory mechanisms of GA on multiple aspects of fruit quality formation at the transcriptional regulation level.

Activity and potential mechanisms of action of persimmon tannins according to their structures: A review

One distinguishing feature of the persimmon, that differentiates it from other fruits, is its high proanthocyanidins content, known as persimmon tannin (PT). Despite the poor absorption of PT in the small intestine, results from animal studies demonstrate that PT has many health benefits. Our goal in this review is to summarize the literature that elucidates the relationship between PT structure and activity. In addition, we also summarize the potential mechanisms underlying the health benefits that result from PT consumption; this includes the hypolipidemic, hypoglycemic, antioxidant, anti-inflammatory, antiradiation, antibacterial and antiviral, detoxification effects on snake venom, and the absorption of heavy metals and dyes. Studies show that PT is a structurally distinct proanthocyanidins that exhibits a high degree of polymerization. It is galloylation-rich and possesses unique A-type interflavan linkages in addition to the more common B-type interflavan bonds. Thus, PT is converted into oligomeric proanthocyanidins by depolymerization strategies, including the nucleophilic substitution reaction, acid hydrolysis, and hydrogenolysis. In addition, multiple health benefits exerted by PT mainly involve the inactivation of lipogenic and intracellular inflammatory signaling pathways, activation of the fatty acid oxidation signaling pathway, regulation of gut microbiota, and highly absorptive properties.

Population Genetic Analysis in Persimmons (Diospyros kaki Thunb.) Based on Genome-Wide Single-Nucleotide Polymorphisms

This study investigated the genetic diversity and population structure of a persimmon (Diospyros kaki Thunb., 2n = 6x = 90) collection in South Korea by evaluating 9751 genome-wide single-nucleotide polymorphisms (SNPs) detected using genotyping-by-sequencing in 93 cultivars. The results of neighbor-joining clustering, principal component analysis, and STRUCTURE analysis based on SNPs indicated clear separation between cultivar groups (pollination-constant nonastringent (PCNA, 40 cultivars), pollination-constant astringent (PCA, 19), pollination-variant nonastringent (PVNA, 23), and the pollination-variant astringent type (PVA, 9)) based on the astringency types, while separation between PVA and PVNA-type cultivars was unclear. Population genetic diversity based on SNPs showed that the proportions of polymorphic SNPs within each group ranged from 99.01% (PVNA) to 94.08% (PVA), and the PVNA group exhibited the highest genetic diversity (He = 3.86 and uHe = 0.397). F (fixation index) values were low ranging from -0.024 (PVA) to 0.176 (PCA) with an average of 0.089, indicating a deficiency of heterozygosity. Analysis of molecular variance (AMOVA) and Fst among cultivar groups indicated that variation within individuals was higher than that among the groups. Pairwise Fst values among the groups ranged from 0.01566 (between PVA and PVNA) to 0.09416 (between PCA and PCNA), indicating a low level of cultivar type differentiation. These findings highlight the potential application of biallelic SNPs in population genetics studies of allopolyploids species and provide valuable insights that may have significant implications for breeding and cultivar identification in persimmon.

Allele-aware chromosome-level genome assembly of the autohexaploid Diospyros kaki Thunb

Artificially improving persimmon (Diospyros kaki Thunb.), one of the most important fruit trees, remains challenging owing to the lack of reference genomes. In this study, we generated an allele-aware chromosome-level genome assembly for the autohexaploid persimmon 'Xiaoguotianshi' (Chinese-PCNA type) using PacBio CCS and Hi-C technology. The final assembly contained 4.52 Gb, with a contig N50 value of 5.28 Mb and scaffold N50 value of 44.01 Mb, of which 4.06 Gb (89.87%) of the assembly were anchored onto 90 chromosome-level pseudomolecules comprising 15 homologous groups with 6 allelic chromosomes in each. A total of 153,288 protein-coding genes were predicted, of which 98.60% were functionally annotated. Repetitive sequences accounted for 64.02% of the genome; and 110,480 rRNAs, 12,297 tRNAs, 1,483 miRNAs, and 3,510 snRNA genes were also identified. This genome assembly fills the knowledge gap in the autohexaploid persimmon genome, which is conducive in the study on the regulatory mechanisms underlying the major economically advantageous traits of persimmons and promoting breeding programs.

Full-length transcriptome profiling for fruit development in Diospyros oleifera using nanopore sequencing

OBJECTIVES

Diospyros oleifera, one of the most economically important Diospyros species, is an ideal model for studying the fruit development of persimmon. While, the lack of whole-transcriptome has hindered the complex transcriptional regulation mechanisms of sugar and tannin during fruit development.

DATA DESCRIPTION

We applied Oxford Nanopore Technologies to six developmental stage of fruit from D. oleifera for use in transcriptome sequencing. As a result of full-length transcriptome sequencing, 55.87 Gb of clean data were generated. After mapping onto the reference genome of D. oleifera, 51,588 full-length collapsing transcripts, including 2,727 new gene loci and 43,223 transcripts, were obtained. Comprehensively annotated, 38,086 of new transcripts were functional annotation, and 972 lncRNAs, 7,159 AS events were predicted. Here, we released the transcriptome database of D. oleifera at different stage of fruit development,which will provide a fundamention of to investigatethe transcript structure, variants and evolution of persimmon.

Investigating phenotypic relationships in persimmon accessions through integrated proteomic and metabolomic analysis of corresponding fruits

Together with phenological and genomic approaches, gel-based and label-free proteomic as well metabolomic procedures were separately applied to plants to highlight differences between ecotypes, to estimate genetic variability within/between organism populations, or to characterize specific mutants/genetically modified lines at metabolic level. To investigate the possible use of tandem mass tag (TMT)-based quantitative proteomics in the above-mentioned contexts and based on the absence of combined proteo-metabolomic studies on Diospyros kaki cultivars, we here applied integrated proteomic and metabolomic approaches to fruits from Italian persimmon ecotypes with the aim to characterize plant phenotypic diversity at molecular level. We identified 2255 proteins in fruits, assigning 102 differentially represented components between cultivars, including some related to pomological, nutritional and allergenic characteristics. Thirty-three polyphenols were also identified and quantified, which belong to hydroxybenzoic acid, flavanol, hydroxycinnamic acid, flavonol, flavanone and dihydrochalcone sub-classes. Heat-map representation of quantitative proteomic and metabolomic results highlighted compound representation differences in various accessions, whose elaboration through Euclidean distance functions and other linkage methods defined dendrograms establishing phenotypic relationships between cultivars. Principal component analysis of proteomic and metabolomic data provided clear information on phenotypic differences/similarities between persimmon accessions. Coherent cultivar association results were observed between proteomic and metabolomic data, emphasizing the utility of integrating combined omic approaches to identify and validate phenotypic relationships between ecotypes, and to estimate corresponding variability and distance. Accordingly, this study describes an original, combined approach to outline phenotypic signatures in persimmon cultivars, which may be used for a further characterization of other ecotypes of the same species and an improved description of nutritional characteristics of corresponding fruits.

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

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

The chromosome-level genome for Toxicodendron vernicifluum provides crucial insights into Anacardiaceae evolution and urushiol biosynthesis

The lacquer tree (Toxicodendron vernicifluum (Stokes) F.A. Barkley) is an important tree with economic, industrial, and medicinal values. Here, we generated the reference genome of T. vernicifluum at the chromosome level with 491.93 Mb in size, in which 98.26% of the assembled contigs were anchored onto 15 pseudochromosomes with the scaffold N50 of 32.97 Mb. Comparative genomic analysis revealed the gene families related to urushiol biosynthesis were expanded, contributing to the ecological fitness and biological adaptability of the lacquer tree. We combined multi-omics data to identify genes that encode key enzymes in the T. vernicifluum urushiol and lignin biosynthetic pathways. Furthermore, the unique active metabolites, such as butin and fisetin, in cultivar lacquers were identified by metabolism profiling. Our work would provide crucial insights into metabolite synthesis such as urushiol and lignin, meanwhile offer a basis for further exploration of the cultivation and breeding of T. vernicifluum and other Anacardiaceae members.

Characterization and phylogenetic analysis of the complete mitochondrial genome sequence of Diospyros oleifera, the first representative from the family Ebenaceae

Plant mitochondrial genomes are a valuable source of genetic information for a better understanding of phylogenetic relationships. However, no mitochondrial genome of any species in Ebenaceae has been reported. In this study, we reported the first mitochondrial genome of an Ebenaceae model plant Diospyros oleifera. The mitogenome was 493,958 bp in length, contained 39 protein-coding genes, 27 transfer RNA genes, and 3 ribosomal RNA genes. The rps2 and rps11 genes were missing in the D. oleifera mt genome, while the rps10 gene was identified. The length of the repetitive sequence in the D. oleifera mt genome was 31 kb, accounting for 6.33%. A clear bias in RNA-editing sites were found in the D. oleifera mt genome. We also detected 28 chloroplast-derived fragments significantly associated with D. oleifera mt genes, indicating intracellular tRNA genes transferred frequently from chloroplasts to mitochondria in D. oleifera. Phylogenetic analysis based on the mt genomes of D. oleifera and 27 other taxa reflected the exact evolutionary and taxonomic status of D. oleifera. Ka/Ks analysis revealed that 95.16% of the protein-coding genes in the D. oleifera mt genome had undergone negative selections. But, the rearrangement of mitochondrial genes has been widely occur among D. oleifera and these observed species. These results will lay the foundation for identifying further evolutionary relationships within Ebenaceae.

The Chromosome-Level Genome of Miracle Fruit (Synsepalum dulcificum) Provides New Insights Into the Evolution and Function of Miraculin

Miracle fruit (Synsepalum dulcificum) is a rare valuable tropical plant famous for a miraculous sweetening glycoprotein, miraculin, which can modify sour flavors to sweet flavors tasted by humans. Here, we present a chromosome-level high-quality genome of S. dulcificum with an assembly genome size of ∼550 Mb, contig N50 of ∼14.14 Mb, and 37,911 annotated protein-coding genes. Phylogenetic analysis revealed that S. dulcificum was most closely related to Camellia sinensis and Diospyros oleifera, and that S. dulcificum diverged from the Diospyros genus ∼75.8 million years ago (MYA), and that C. sinensis diverged from Synsepalum ∼63.5 MYA. Ks assessment and collinearity analysis with S. dulcificum and other species suggested that a whole-genome duplication (WGD) event occurred in S. dulcificum and that there was good collinearity between S. dulcificum and Vitis vinifera. On the other hand, transcriptome and metabolism analysis with six tissues containing three developmental stages of fleshes and seeds of miracle fruit revealed that Gene Ontology (GO) terms and metabolic pathways of "cellular response to chitin," "plant-pathogen interaction," and "plant hormone signal transduction" were significantly enriched during fruit development. Interestingly, the expression of miraculin (Chr10G0299340) progressively increased from vegetative organs to reproductive organs and reached an incredible level in mature fruit flesh, with an fragments per kilobase of transcript per million (FPKM) value of ∼113,515, which was the most highly expressed gene among all detected genes. Combining the unique signal peptide and the presence of the histidine-30 residue together composed the main potential factors impacting miraculin's unique properties in S. dulcificum. Furthermore, integrated analysis of weighted gene coexpression network analysis (WGCNA), enrichment and metabolite correlation suggested that miraculin plays potential roles in regulating plant growth, seed germination and maturation, resisting pathogen infection, and environmental pressure. In summary, valuable genomic, transcriptomic, and metabolic resources provided in this study will promote the utilization of S. dulcificum and in-depth research on species in the Sapotaceae family.

UV-B induces the expression of flavonoid biosynthetic pathways in blueberry (Vaccinium corymbosum) calli

Ultraviolet-B (UV-B) radiation is an environmental signal that affects the accumulation of secondary metabolites in plants. In particular, UV-B promotes flavonoid biosynthesis, leading to improved fruit quality. To explore the underlying molecular mechanism, we exposed blueberry (Vaccinium corymbosum) calli to UV-B radiation and performed a transcriptome deep sequencing (RNA-seq) analysis to identify differentially expressed genes (DEGs). We detected 16,899 DEGs among different treatments, with the largest number seen after 24 h of UV-B exposure relative to controls. Functional annotation and enrichment analysis showed a significant enrichment for DEGs in pathways related to plant hormone signal transduction and phenylpropanoid and flavonoid biosynthesis. In agreement with the transcriptome data, flavonol, anthocyanin and proanthocyanidin accumulated upon UV-B radiation, and most DEGs mapping to the phenylpropanoid and flavonoid biosynthetic pathways using the KEGG mapper tool were upregulated under UV-B radiation. We also performed a weighted gene co-expression network analysis (WGCNA) to explore the relationship among genes involved in plant hormone signal transduction, encoding transcription factors or participating in flavonoid biosynthesis. The transcription factors VcMYBPA1, MYBPA2.1, MYB114, MYBA2, MYBF, and MYB102 are likely activators, whereas MYB20, VcMYB14, MYB44, and VcMYB4a are inhibitors of the flavonoid biosynthetic pathway, as evidenced by the direction of correlation between the expression of these MYBs and flavonoid biosynthesis-related genes. The transcription factors bHLH74 and bHLH25 might interact with MYB repressors or directly inhibited the expression of flavonoid biosynthetic genes to control flavonoid accumulation. We also observed the downregulation of several genes belonging to the auxin, gibberellin and brassinosteroid biosynthetic pathways, suggesting that MYB inhibitors or activators are directly or indirectly regulated to promote flavonoid biosynthesis under UV-B radiation.

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.

The regulatory role of gibberellin related genes DKGA2ox1 and MIR171f_3 in persimmon dwarfism

'Nantongxiaofangshi' (Diospyros kaki Thunb., D. kaki Thunb.) is a local cultivar of persimmon with dwarf-like traits in Jiangsu, China. Closely spaced planting afforded by dwarfism is usually one of the most important ways to promote fruit cultivation and production. However, the understanding of dwarfism in D. kaki Thunb. is very limited at the molecular level, which hinders the further increase of the fruit production. In this work, a persimmon transgenic system was successfully established, and the field experiments of grafting phenotype were carried out. The results showed that D. kaki Thunb. could be used as an interstock to induce dwarfing in grafted scions, and the dwarf character was better when interstock lengths were between 20 and 25 cm. Furthermore, the key genes related to dwarfism in D. kaki Thunb. were screened and verified, and subsequently, the regulatory role of related genes in persimmon dwarfism was figured out. It was found that the gene encoding gibberellin 2-oxidase-1 (DkGA2ox1) involved in GA biosynthesis was associated with the dwarfing in D. kaki Thunb. Overexpression of DkGA2ox1 in Diospyros lotus resulted in a typical dwarf phenotype. Meanwhile, the microRNA data showed that the miR171f_3 demonstrated the active involvement in GA pathway response in persimmon dwarfism. DkGA2ox1 and MIR171f_3, as two highly expressed genes in D. kaki Thunb. interstock, could be used as stimulus signals to affect the content of GA in scion, however, the specific transmission mechanism still needs to be further explored. Ultimately, the bioactive GA level was decreased, resulting in the scion dwarfism.

The chromosome-scale reference genome of Rubus chingii Hu provides insight into the biosynthetic pathway of hydrolyzable tannins

Rubus chingii Hu (Fu-Pen-Zi), a perennial woody plant in the Rosaceae family, is a characteristic traditional Chinese medicinal plant because of its unique pharmacological effects. There are abundant hydrolyzable tannin (HT) components in R. chingii that provide health benefits. Here, an R. chingii chromosome-scale genome and related functional analysis provide insights into the biosynthetic pathway of HTs. In total, sequence data of 231.21 Mb (155 scaffolds with an N50 of 8.2 Mb) were assembled into seven chromosomes with an average length of 31.4 Mb, and 33 130 protein-coding genes were predicted, 89.28% of which were functionally annotated. Evolutionary analysis showed that R. chingii was most closely related to Rubus occidentalis, from which it was predicted to have diverged 22.46 million years ago (Table S8). Comparative genomic analysis showed that there was a tandem gene cluster of UGT, carboxylesterase (CXE) and SCPL genes on chromosome 02 of R. chingii, including 11 CXE, eight UGT, and six SCPL genes, which may be critical for the synthesis of HTs. In vitro enzyme assays indicated that the proteins encoded by the CXE (LG02.4273) and UGT (LG02.4102) genes have tannin hydrolase and gallic acid glycosyltransferase functions, respectively. The genomic sequence of R. chingii will be a valuable resource for comparative genomic analysis within the Rosaceae family and will be useful for understanding the biosynthesis of HTs.

Is speciation an unrelenting march to reproductive isolation?

Speciation is often portrayed as an "incomplete" or "incipient" process if two groups of organisms, technically distinguishable either by morphology or genetics, can exchange genes. The ultimate outcome of diversification, given this perspective, is complete reproductive isolation. But an increasing amount of evidence suggests that speciation is rarely complete and inter-fertility between different taxonomically accepted species is consistently maintained. In this issue of Molecular Ecology, Linan et al. (2021) provide results that bridge evolutionary processes from populations to phylogenies that indicate suites of closely related tree species in the Mascarene Islands actively exchange genes, evolving as a nested set of syngameons with a hierarchical pattern of interfertility. The deep insight into diversification provided by this study is particularly powerful because of the genomic scale of the data and the complete taxonomic sampling of an island clade evolving in situ. The prevalence of syngameon dynamics in a broad range of organisms indicates that we should adopt a fluid and comprehensive approach to defining evolutionary units for conservation and research. We should move beyond focusing on single endangered species in evolutionary and ecological isolation from other species but consider the entire network of potentially interfertile species and the potential for future adaptation and innovation, particularly in a human dominated world.

Transcriptome Analysis Revealed the Roles of Carbohydrate Metabolism on Differential Acetaldehyde Production Capacity in Persimmon Fruit in Response to High-CO2 Treatment

Persimmon (Diospyros kaki Thunb.) fruit is unique due to the continuous accumulation of soluble tannins during fruit development in most cultivars, which causes undesired astringency. High-CO₂ treatment was the most effective widely used method for astringency removal. However, differential effects of high-CO₂ treatment between cultivars were observed and the molecular basis remained inclusive. Previously, one cultivar ("Luoyangfangtianshengshi," LYFTSS) showed rapid deastringency, while two cultivars ("Shijiazhuanglianhuashi," SJZLHS; "Laopige," LPG) showed slow deastringency in response to high-CO₂ (95% CO₂) treatment. In this study, the metabolites (acetaldehyde and ethanol) related to deastringency were further analyzed and both acetaldehyde and ethanol were higher in SJZLHS and LYFTSS than that in LPG, where acetaldehyde was undetectable. Based on the RNA-seq data, the weighted gene coexpression network analysis (WGCNA) revealed that one module, comprised of 1773 unigenes, significantly correlated with the contents of acetaldehyde and ethanol (P < 0.001). Further analysis based on the acetaldehyde metabolism pathway indicated that the differentially expressed structural genes, including previously characterized DkADH and DkPDC and also their upstream members (e.g., PFK, phosphofructokinase), showed positive correlations with acetaldehyde production. Quantitative analysis of the precursor substances indicated that sucrose, glucose, and fructose exhibited limited differences between cultivar except for malic acid. However, the content of malic acid is much less than the total soluble sugar content. To verify the correlations between these genes and acetaldehyde production, the fruit from 14 more cultivars were collected and treated with high CO₂. After the treatment, acetaldehyde contents in different cultivars ranked in 30.4-255.5 μg/g FW. Real-time polymerase chain reaction (PCR) and correlation analysis indicated that the EVM0002315 (PFK) gene, belonging to carbohydrate metabolism, was significantly correlated with acetaldehyde content in fruit. Thus, it could be proposed that the differentially expressed carbohydrate metabolism related genes (especially PFK) are the basis for the variance of acetaldehyde production among different persimmon cultivars.

lepidium-like, a Naturally Occurring Mutant of Capsella bursa-pastoris, and Its Implications on the Evolution of Petal Loss in Cruciferae

Naturally occurring mutants whose phenotype recapitulates the changes that distinguish closely related species are of special interest from the evolutionary point of view. They can give a key about the genetic control of the changes that led to speciation. In this study, we described lepidium-like (lel), a naturally occurring variety of an allotetraploid species Capsella bursa-pastoris that is characterized by the typical loss of all four petals. In some cases, one or two basal flowers in the raceme had one or two small petals. The number and structure of other floral organs are not affected. Our study of flower development in the mutant showed that once initiated, petals either cease further development and cannot be traced in anthetic flowers or sometimes develop to various degrees. lel plants showed an earlier beginning of floral organ initiation and delayed petal initiation compared to the wild-type plants. lel phenotype has a wide geographical distribution, being found at the northern extremity of the species range as well as in the central part. The genetic analysis of inheritance demonstrated that lel phenotype is controlled by two independent loci. While the flower in the family Cruciferae generally has a very stable structure (i.e., four sepals, four petals, six stamens, and two carpels), several deviations from this ground plan are known, in particular in the genus Lepidium, C. bursa-pastoris is an emerging model for the study of polyploidy (which is also very widespread in Cruciferae); the identification and characterization of the apetalous mutant lays a foundation for further research of morphological evolution in polyploids.

Genome-Wide Identification of Tannase Genes and Their Function of Wound Response and Astringent Substances Accumulation in Juglandaceae

Tannins are important polyphenol compounds with different component proportions in different plant species. The plants in the Juglandaceae are rich in tannins, including condensed tannins and hydrolyzable tannins. In this study, we identified seven tannase genes (TAs) responsible for the tannin metabolism from walnut, pecan, and Chinese hickory, and three nut tree species in the Juglandaceae, which were divided into two groups. The phylogenetic and sequence analysis showed that TA genes and neighboring clade genes (TA-like genes) had similar sequences compared with other carboxylesterase genes, which may be the origin of TA genes produced by tandem repeat. TA genes also indicated higher expressions in leaf than other tissues and were quickly up-regulated at 3 h after leaf injury. During the development of the seed coat, the expression of the synthesis-related gene GGTs and the hydrolase gene TAs was continuously decreased, resulting in the decrease of tannin content in the dry sample of the seed coat of Chinese hickory. However, due to the reduction in water content during the ripening process, the tannin content in fresh sample increased, so the astringent taste was obvious at the mature stage. In addition, the CcGGTs' expression was higher than CiGGTs in the initiation of development, but CcTAs continued to be down-regulated while CiTA2a and CiTA2b were up-regulated, which may bring about the significant differences in tannin content and astringent taste between Chinese hickory and pecan. These results suggested the crucial role of TAs in wound stress of leaves and astringent ingredient accumulation in seed coats of two nut tree species in the Juglandaceae.

A high-quality chromosomal genome assembly of Diospyros oleiferaCheng

BACKGROUND

Diospyros oleifera Cheng, of the family Ebenaceae, is an economically important tree. Phylogenetic analyses indicate that D. oleifera is closely related to Diospyros kaki Thunb. and could be used as a model plant for studies of D. kaki. Therefore, development of genomic resources of D. oleifera will facilitate auxiliary assembly of the hexaploid persimmon genome and elucidate the molecular mechanisms of important traits.

FINDINGS

The D. oleifera genome was assembled with 443.6 Gb of raw reads using the Pacific Bioscience Sequel and Illumina HiSeq X Ten platforms. The final draft genome was ∼812.3 Mb and had a high level of continuity with N50 of 3.36 Mb. Fifteen scaffolds corresponding to the 15 chromosomes were assembled to a final size of 721.5 Mb using 332 scaffolds, accounting for 88.81% of the genome. Repeat sequences accounted for 54.8% of the genome. By de novo sequencing and analysis of homology with other plant species, 30,530 protein-coding genes with an average transcript size of 7,105.40 bp were annotated; of these, 28,580 protein-coding genes (93.61%) had conserved functional motifs or terms. In addition, 171 candidate genes involved in tannin synthesis and deastringency in persimmon were identified; of these chalcone synthase (CHS) genes were expanded in the D. oleifera genome compared with Diospyros lotus, Camellia sinensis, and Vitis vinifera. Moreover, 186 positively selected genes were identified, including chalcone isomerase (CHI) gene, a key enzyme in the flavonoid-anthocyanin pathway. Phylogenetic tree analysis indicated that the split of D. oleifera and D. lotus likely occurred 9.0 million years ago. In addition to the ancient γ event, a second whole-genome duplication event occurred in D. oleifera and D. lotus.

CONCLUSIONS

We generated a high-quality chromosome-level draft genome for D. oleifera, which will facilitate assembly of the hexaploid persimmon genome and further studies of major economic traits in the genus Diospyros.