Mitochondrial genome features and systematic evolution of diospyros kaki thunb 'Taishuu'

BACKGROUND

'Taishuu' has a crisp texture, abundant juice, and sweet flavor with hints of cantaloupe. The availability of mitochondrial genome data of Diospyros species is far from the known number of species.

RESULTS

The sequencing data were assembled into a closed circular mitochondrial chromosome with a 421,308 bp length and a 45.79% GC content. The mitochondrial genome comprised 40 protein-coding, 24 tRNA, and three rRNA genes. The most common codons for arginine (Arg), proline (Pro), glycine (Gly), tryptophan (Trp), valine (Val), alanine (Ala), and leucine (Leu) were AGA, CCA, GGA, UGG, GUA, GCA, and CUA, respectively. The start codon for cox1 and nad4L protein-coding genes was ACG (ATG), whereas the remaining protein-coding genes started with ATG. There are four types of stop codons: CGA, TAA, TAG, and TGA, with TAA being the most frequently used stop codon (45.24%). In the D. kaki Thunb. 'Taishuu' mitochondrial genome, a total of 645 repeat sequences were identified, including 125 SSRs, 7 tandem repeats, and 513 dispersed repeats. Collinearity analysis revealed a close relationship between D. kaki Thunb. 'Taishuu' and Diospyros oleifera, with conserved homologous gene fragments shared among these species in large regions of the mitochondrial genome. The protein-coding genes ccmB and nad4L were observed to undergo positive selection. Analysis of homologous sequences between chloroplasts and mitochondria identified 28 homologous segments, with a total length of 24,075 bp, accounting for 5.71% of the mitochondrial genome. These homologous segments contain 8 annotated genes, including 6 tRNA genes and 2 protein-coding genes (rrn18 and ccmC). There are 23 homologous genes between chloroplasts and nuclei. Mitochondria, chloroplasts, and nuclei share two homologous genes, which are trnV-GAC and trnW-CCA.

CONCLUSION

In conclusion, a high-quality chromosome-level draft genome for D. kaki was generated in this study, which will contribute to further studies of major economic traits in the genus Diospyros.

Genetic insights into the dissolution of dioecy in diploid persimmon Diospyros oleifera Cheng

BACKGROUND

Dioecy, a sexual system of single-sexual (gynoecious/androecious) individuals, is rare in flowering plants. This rarity may be a result of the frequent transition from dioecy into systems with co-sexual individuals.

RESULTS

In this study, co-sexual expression (monoecy and hermaphroditic development), previously thought to be polyploid-specific in Diospyros species, was identified in the diploid D. oleifeara historically. We characterized potential genetic mechanisms that underlie the dissolution of dioecy to monoecy and andro(gyno)monoecy, based on multiscale genome-wide investigations of 150 accessions of Diospyros oleifera. We found all co-sexual plants, including monoecious and andro(gyno)monoecious individuals, possessed the male determinant gene OGI, implying the presence of genetic factors controlling gynoecia development in genetically male D. oleifera. Importantly, discrepancies in the OGI/MeGI module were found in diploid monoecious D. oleifera compared with polyploid monoecious D. kaki, including no Kali insertion on the promoter of OGI, no different abundance of smRNAs targeting MeGI (a counterpart of OGI), and no different expression of MeGI between female and male floral buds. On the contrary, in both single- and co-sexual plants, female function was expressed in the presence of a genome-wide decrease in methylation levels, along with sexually distinct regulatory networks of smRNAs and their targets. Furthermore, a genome-wide association study (GWAS) identified a genomic region and a DUF247 gene cluster strongly associated with the monoecious phenotype and several regions that may contribute to andromonoecy.

CONCLUSIONS

Collectively, our findings demonstrate stable breakdown of the dioecious system in D. oleifera, presumably also a result of genomic features of the Y-linked region.

Comparative Metabolomic and Transcriptomic Analyses Reveal Distinct Ascorbic Acid (AsA) Accumulation Patterns between PCA and PCNA Persimmon Developing Fruit

Persimmon fruit has a high nutritional value and significantly varies between pollination-constant astringent (PCA) and pollination-constant non-astringent (PCNA) persimmons. The astringency type affects sugar, flavonoids, and tannin accumulation and is well known in persimmon fruit. However, the impact of the fruit astringency type on ascorbic acid (AsA) accumulation is limited. In this study, typical PCA varieties ('Huojing' and 'Zhongshi5') and PCNA varieties ('Yohou' and 'Jiro') of persimmon fruit were sampled at four developing stages (S1-S4) to provide valuable information on AsA content variation in PCA and PCNA persimmon. Persimmon fruit is rich in ascorbic acid; the AsA content of the four varieties 'Zhongshi5', 'Huojing', 'Jiro', and 'Youhou' mature fruit reached 104.49, 48.69, 69.69, and 47.48 mg/100 g. Fruit of the same astringency type persimmon showed a similar AsA accumulation pattern. AsA content was significantly higher in PCA than PCNA fruit at S1-S3. The initial KEGG analysis of metabolites showed that galactose metabolism is the major biosynthetic pathway of AsA in persimmon fruit. There were significant differences in galactose pathway-related metabolite content in developing PCA and PCNA fruit, such as Lactose, D-Tagatose, and D-Sorbitol content in PCA being higher than that of PCNA. Combined gene expression and WGCNA analyses showed that the expression of the GME (evm.TU.contig4144.37) gene was higher in PCA-type than in PCNA-type fruit in S1-S3 and exhibited the highest correlation with AsA content (r = 690 **, p < 0.01). Four hub genes, including the DNA methylation gene, methyltransferase gene, F-box, and Actin-like Protein, were identified as potential regulators of the GME gene. These results provide basic information on how astringency types affect AsA accumulation and will provide valuable information for further investigation on AsA content variation in persimmon fruit.

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.

Genetic basis of lineage-specific evolution of fruit traits in hexaploid persimmon

Frequent polyploidization events in plants have led to the establishment of many lineage-specific traits representing each species. Little is known about the genetic bases for these specific traits in polyploids, presumably due to plant genomic complexity and their difficulties in applying genetic approaches. Hexaploid Oriental persimmon (Diospyros kaki) has evolved specific fruit characteristics, including wide variations in fruit shapes and astringency. In this study, using whole-genome diploidized/quantitative genotypes from ddRAD-Seq data of 173 persimmon cultivars, we examined their population structures and potential correlations between their structural transitions and variations in nine fruit traits. The population structures of persimmon cultivars were highly randomized and not substantially correlated with the representative fruit traits focused on in this study, except for fruit astringency. With genome-wide association analytic tools considering polyploid alleles, we identified the loci associated with the nine fruit traits; we mainly focused on fruit-shape variations, which have been numerically characterized by principal component analysis of elliptic Fourier descriptors. The genomic regions that putatively underwent selective sweep exhibited no overlap with the loci associated with these persimmon-specific fruit traits. These insights will contribute to understanding the genetic mechanisms by which fruit traits are independently established, possibly due to polyploidization events.

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.

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.

DkmiR397 Regulates Proanthocyanidin Biosynthesis via Negative Modulating DkLAC2 in Chinese PCNA Persimmon

Persimmon fruits accumulate a large amount of proanthocyanidins (PAs), which makes an astringent sensation. Proanthocyanidins (PAs) are the polymers of flavan-3-ols stored in plant vacuoles under laccase activation. A laccase gene, DkLAC2, is putatively involved in PAs biosynthesis and regulated by microRNA (DkmiR397) in persimmon. However, the polymerization of PAs in association with miRNA397 still needs to be explored in persimmon. Here, we identified pre-DkmiR397 and its target gene DkLAC2 in ‘Eshi 1’ persimmon. Histochemical staining with GUS and dual luciferase assay both confirmed DkmiR397-DkLAC2 binding after co-transformation in tobacco leaves. Diverse expression patterns of DkLAC2 and DkmiR397 were exhibited during persimmon fruit development stages. Moreover, a contrasting expression pattern was also observed after the combined DkLAC2-miR397 transformation in persimmon leaves, suggesting that DkmiR397 might be a negative regulator of DkLAC2. Similarly, the transient transformation of DkmiR397 in persimmon fruit discs in vitro also reduced PA accumulation by repressing DkLAC2, whereas the up-regulation of DkLAC2 increased the accumulation of PAs by short tandem target mimic STTM-miR397. A similar expression pattern was observed when overexpressing of DkLAC2 in Arabidopsis wild type (WT) and overexpression of DkLAC2, DkmiR397 in persimmon leaf callus. Our results revealed that the role of DkmiR397 repressed the expression of DkLAC2 concerning PA biosynthesis, providing a potential target for the manipulation of PAs metabolism in persimmon.

Reinvention of hermaphroditism via activation of a RADIALIS-like gene in hexaploid persimmon

In flowering plants, different lineages have independently transitioned from the ancestral hermaphroditic state into and out of various sexual systems¹. Polyploidizations are often associated with this plasticity in sexual systems^2,3. Persimmons (the genus Diospyros) have evolved dioecy via lineage-specific palaeoploidizations. More recently, hexaploid D. kaki has established monoecy and also exhibits reversions from male to hermaphrodite flowers in response to natural environmental signals (natural hermaphroditism, NH), or to artificial cytokinin treatment (artificial hermaphroditism, AH). We sought to identify the molecular pathways underlying these polyploid-specific reversions to hermaphroditism. Co-expression network analyses identified regulatory pathways specific to NH or AH transitions. Surprisingly, the two pathways appeared to be antagonistic, with abscisic acid and cytokinin signalling for NH and AH, respectively. Among the genes common to both pathways leading to hermaphroditic flowers, we identified a small-Myb RADIALIS-like gene, named DkRAD, which is specifically activated in hexaploid D. kaki. Consistently, ectopic overexpression of DkRAD in two model plants resulted in hypergrowth of the gynoecium. These results suggest that production of hermaphrodite flowers via polyploidization depends on DkRAD activation, which is not associated with a loss-of-function within the existing sex determination pathway, but rather represents a new path to (or reinvention of) hermaphroditism.

Genome Sequence Resource of the Causal Agent of Persimmon Anthracnose, Colletotrichum horii Strain SD010 from China

Colletotrichum horii is a main causal agent of persimmon (Diospyros kaki) anthracnose and is distributed widely in persimmon-producing areas of the world. Here, we report the first high-quality draft genome sequence of C. horii strain SD010. This will provide a reference for understanding adaptive evolution of genome structure, genes, and population diversity among members of the C. gloeosporioides species complex, and also help in understanding the mechanisms of host-pathogen interactions and improve management strategies of anthracnose.

The regulation of DKGA2ox1 and miR171f_3 in scion dwarfing with Diospyros kaki Thunb. cv. 'Nan-tong-xiao-fang-shi' as interstocks

High-throughput sequencing and yeast one and two-hybrid library screening reveal that DKGA2ox1 and miR171f_3 are involved in the regulation of scion dwarfing with 'Nan-tong-xiao-fang-shi' as interstocks. Diospyros kaki Thunb. cv. Nan-tong-xiao-fang-shi ('Nan-tong-xiao-fang-shi') interstocks play a critical role in the scion dwarfing. However, the understanding of the molecular signaling pathways that regulate the scion dwarfing with 'Nan-tong-xiao-fang-shi' as interstocks remain unclear. In this work, the regulatory network in the scion dwarfing with 'Nan-tong-xiao-fang-shi' as interstocks was identified. Using a yeast one-hybrid library screening, luciferase activity analysis, luciferase complementation imaging assays and GFP signal detection, a SPL transcription factor named Diospyros kaki SPL (DKSPL), potentially functioning as a transcriptional activator of the Diospyros kaki GA2ox1 (DKGA2ox1) gene, was identified as a key stimulating factor in the persimmon growth and development. The DKSPL was found in the nucleus, and might play a role in the transcriptional regulation system. A microRNA named miR171f_3 was identified, which might act as a negative regulator of Diospyros kaki SCR (DKSCR) in persimmon. The interactions between DKSCR and seven proteins were experimentally validated with a yeast two-hybrid library screening. Compared to the non-grafted wildtype persimmon, the tissue section of graft union healed well due to the increased expression of cinnamyl-alcohol dehydrogenase. These results indicate that DKGA2ox1 and miR171f_3 may co-promote the scion dwarfing by plant hormone signal transduction pathways.

The whole-genome sequence of the novel yeast species Metschnikowia persimmonesis isolated from medicinal plant Diospyros kaki Thunb

The new yeast Metschnikowia persimmonesis KCTC 12991BP (KIOM G15050 strain) exhibits strong antimicrobial activity against some pathogens. This activity may be related to the medicinal profile of secondary metabolites that could be found in the genome of this species. Therefore, to explore its future possibility of producing some beneficial activities, including medicinal ability, we report high-quality whole-genome assembly of M. persimmonesis produced by PacBio RS II sequencer. The final draft assembly consisted of 16 scaffolds with GC content of 45.90% and comprised a fairly complete set (82.8%) of BUSCO result using Saccharomycetales lineage data set. The total length of the genome was 16.473 Mb, with a scaffold N50 of 1.982 Mb. Annotation of the M. persimmonesis genome revealed presence of 7029 genes and 6939 functionally annotated proteins. Based on the analysis of phylogenetic relationship and the average nucleotide identities, M. persimmonesis was proved to a novel species within the Metschnikowia genus. This finding is expected to significantly contribute to the discovery of high-value natural products from M. persimmonesis as well as for genome biology and evolution comparative analysis within Metschnikowia species.

Transcription factors DkBZR1/2 regulate cell wall degradation genes and ethylene biosynthesis genes during persimmon fruit ripening

BRASSINAZOLE RESISTANT (BZR) transcription factors are critical components of the brassinosteroid signalling pathway, but their possible roles in fruit ripening have rarely been reported. In this study, four BZR sequences were isolated from persimmon fruit. Among the four BZR genes, DkBZR1/2 were expressed in persimmon fruit; DkBZR1 protein amount decreased and dephosphorylated DkBZR2 gradually accumulated during the storage period. DkBZR1/2 proteins were localized in both the nucleus and cytoplasm and accumulated in the nucleus after 24-epibrassinolide treatment. DkBZR1 suppressed the transcription of Diospyros kaki endo-1,4-betaglucanase 1 (DkEGase1) and 1-aminocyclopropane-1-carboxylate synthase 1 (DkACS1) by binding to the BR response element (BRRE) in their promoters, and DkBZR2 activated the transcription of pectate lyase 1 (DkPL1) and 1-aminocyclopropane-1-carboxylate oxidase 2 (DkACO2) by binding to the E-box motif in their promoters. Transient overexpression of DkBZR2 promoted the conversion of acid-soluble pectin to water-soluble pectin and increased ethylene production in persimmon fruit. Our findings indicate that DkBZR1 and DkBZR2 serve as repressors and activators of persimmon fruit ripening, respectively.

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.

Comparative transcriptome analysis reveals regulatory network and regulators associated with proanthocyanidin accumulation in persimmon

BACKGROUND

Proanthocyanidins (PAs) are important plant secondary metabolites that confer flavor, nutritional value, and resistance to pathogens. Persimmon is one of the PA richest crops. Mature fruits can be inedible because of the astringency caused by high PA levels and need to go through a de-astringency treatment before consumption. The molecular basis for PA accumulation is poorly known, particularly transcriptional regulators. We characterised three genotypes ('Luotiantianshi' (LT), 'Mopanshi' (MP), and 'Youhou' (YH)) with different PA accumulation patterns using an approach that combined PacBio full-length sequencing and Illumina-based RNA sequencing to build high-quality full-length transcriptomes. Additionally, we analysed transcriptome dynamics of the three genotypes (LT, MP, and YH) at four key fruit developmental stages.

RESULTS

A total of 96,463 transcripts were obtained. We identified 80,075 protein-coding sequences (CDSs), 71,137 simple sequence repeats (SSRs), and 27,845 long noncoding RNAs (lncRNAs). Pearson correlation coefficient (PCC), principal component analysis (PCA), and differentially expressed transcripts (DETs) analyses indicated that the four different developmental stages within a genotype exhibited similar transcriptome activities. A total of 2,164 transcripts specific to each fruit developmental stage were detected. The transcripts specific to early stages were attributed to phenylpropanoid and flavonoid biosynthesis. Co-expression network analyses revealed MEbrown and MEblue modules were strongly associated to PA accumulation. From these two modules, 20 hub TFs are potential regulators for PA accumulation. Among them, Cluster_78388 (SBP protein), Cluster_63454 (bZIP protein), and Cluster_66595 (MYB protein) appear to involve in the PA biosynthesis in Chinese genotypes.

CONCLUSIONS

This is the first high-quality reference transcriptome for commercial persimmon. Our work provides insights into the molecular pathways underlying PA accumulation and enhances our global understanding of transcriptome dynamics throughout fruit development.

DkMYB14 is a bifunctional transcription factor that regulates the accumulation of proanthocyanidin in persimmon fruit

Proanthocyanidins (PAs) are phenolic secondary metabolites that contribute to the protection of plant and human health. Persimmon (Diospyros kaki Thunb.) can accumulate abundant PAs in fruit, which cause a strong sensation of astringency. Proanthocyanidins can be classified into soluble and insoluble PAs; the former cause astringency but the latter do not. Soluble PAs can be converted into insoluble PAs upon interacting with acetaldehydes. We demonstrate here that DkMYB14, which regulates the accumulation of PA in persimmon fruit flesh, is a bifunctional transcription factor that acts as a repressor in PA biosynthesis but becomes an activator when involved in acetaldehyde biosynthesis. Interestingly, both functions contribute to the elimination of astringency by decreasing PA biosynthesis and promoting its insolubilization. We show that the amino acid Gly39 in the R2 domain and the ethylene response factor-associated amphiphilic repression-like motif in the C-terminal of DkMYB14 are essential for the regulation of both PA and acetaldehyde synthesis. The repressive function of DkMYB14 was lost after the mutation of either motif, and all activities of DkMYB14 were eliminated following the mutation of both motifs. Our results demonstrate that DkMYB14 functions as both a transcriptional activator and a repressor, directly repressing biosynthesis of PA and promoting its insolubilization, resulting in non-astringency in persimmon.

Phytohormone and integrated mRNA and miRNA transcriptome analyses and differentiation of male between hermaphroditic floral buds of andromonoecious Diospyros kaki Thunb

BACKGROUND

Persimmon (Diospyros kaki Thunb.) has various labile sex types, and studying its sex differentiation can improve breeding efficiency. However, studies on sexual regulation patterns in persimmon have focused mainly on monoecy and dioecy, whereas little research has been published on andromonoecy. In order to reveal the sex differentiation regulation mechanism of andromonoecious persimmon, we performed histological and cytological observations, evaluated OGI and MeGI expression and conducted phytohormones assays and mRNA and small RNA transcriptome analyses of the male and hermaphroditic floral buds of the andromonoecious persimmon 'Longyanyeshi 1'.

RESULTS

Stages 2 and 4 were identified as the critical morphological periods for sex differentiation of 'Longyanyeshi 1' by histological and cytological observation. At both stages, OGI was differentially expressed in male and hermaphroditic buds, but MeGI was not. This was different from their expressions in dioecious and monoecious persimmons. Meantime, the results of phytohormones assays showed that high IAA, ABA, GA3, and JA levels at stage 2 may have promoted male floral bud differentiation. However, high JA levels at stage 4 and high ZT levels at stages 2 and 4 may have promoted hermaphroditic floral bud differentiation. In these phytohormone biosynthesis and signaling pathways, 52 and 54 differential expression genes (including Aux/IAA, ARFs, DELLA, AHP, A-ARR, B-ARR, CYP735A, CRE1, PP2C, JAZ, MYC2, COI1, CTR1, SIMKK, ACO, and MPK6) were identified, respectively. During the development of male floral buds, five metacaspases genes may have been involved in pistil abortion. In addition, MYB, FAR1, bHLH, WRKY, and MADS transcription factors might play important roles in persimmon floral bud sex differentiation. Noteworthy, miR169v_1, miR169e_3, miR319_1, and miR319 were predicted to contribute to phytohormone biosynthesis and signaling pathways and floral organogenesis and may also regulate floral bud sex differentiation.

CONCLUSION

The present study revealed the differences in morphology and phytohormones content between male and hermaphroditic floral buds of 'Longyanyeshi 1' during the process of sex differentiation, and identified a subset of candidate genes and miRNAs putatively associated with its sex differentiation. These findings can provide a foundation for molecular regulatory mechanism researching on andromonoecious persimmon.

Generating long-read sequences using Oxford Nanopore Technology from Diospyros celebica genomic DNA

OBJECTIVES

Development of sequencing technology has opened up vast opportunities for tree genomic research in the tropics. One of the aforesaid technologies named ONT (Oxford Nanopore Technology) has attracted researchers in undertaking testings and experiments due to its affordability and accessibility. To the best of our knowledge, there has been no published reports on the use of ONT for genomic analysis of Indonesian tree species. This progress is promising for further improvement in order to acquire more genomic data for research purposes. Therefore, the present study was carried out to determine the effectiveness of ONT in generating long-read DNA sequences using DNA isolated from leaves and wood cores of Macassar ebony (Diospyros celebica Bakh.).

DATA DESCRIPTION

Long-read sequences data of leaves and wood cores of Macassar ebony were generated by using the MinION device and MinKnow v3.6.5 (ONT). The obtained data, as the first long-read sequence dataset for Macassar ebony, is of great importance to conserve the genetic diversity, understanding the molecular mechanism, and sustainable use of plant genetic resources for downstream applications.

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.

Overexpression of the persimmon abscisic acid β-glucosidase gene (DkBG1) alters fruit ripening in transgenic tomato

β-Glucosidases (BG) are present in many plant tissues. Among these, abscisic acid (ABA) β-glucosidases are thought to take part in the adjustment of cellular ABA levels, however the role of ABA-BG in fruits is still unclear. In this study, through RNA-seq analysis of persimmon fruit, 10 full-length DkBG genes were isolated and were all found to be expressed. In particular, DkBG1 was highly expressed in persimmon fruits with a maximum expression 95 days after full bloom (DAFD). We verified that, in vitro, DkBG1 protein can hydrolyze ABA-glucose ester (ABA-GE) to release free ABA. Compared with wild-type, tomato plants that overexpressed DkBG1 significantly upregulated the expression of ABA receptor PYL3/7 genes and showed typical symptoms of ABA hypersensitivity in fruits. DkBG1 overexpression (DkBG1-OE) accelerated fruit ripening onset by 3-4 days by increasing ABA levels at the pre-breaker stage and induced early ethylene release compared with wild-type fruits. DkBG1-OE altered the expression of ripening regulator NON-RIPENING (NOR) and its target genes; this in turn altered fruit quality traits such as coloration. Our results demonstrated that DkBG1 plays an important role in fruit ripening and quality by adjusting ABA levels via hydrolysis of ABA-GE.

High CO2/hypoxia-induced softening of persimmon fruit is modulated by DkERF8/16 and DkNAC9 complexes

Most persimmon (Diospyros kaki) cultivars are astringent and require post-harvest deastringency treatments such as 95% CO2 (high-CO2 treatment) to make them acceptable to consumers. High-CO2 treatment can, however, also induce excessive softening, which can be reduced by adding 1-methylcyclopropene (1-MCP). Previous studies have shown that genes encoding the ETHYLENE RESPONSE FACTORS (ERFs) DkERF8/16/19 can trans-activate xyloglucan endotransglycosylase/hydrolase (DkXTH9), which encodes the cell wall-degrading enzyme associated with persimmon fruit softening. In this study, RNA-seq data between three treatments were compared, namely high-CO2, high-CO2+1-MCP, and controls. A total of 227 differentially expressed genes, including 17 transcription factors, were predicted to be related to persimmon post-deastringency softening. Dual-luciferase assays indicated that DkNAC9 activated the DkEGase1 promoter 2.64-fold. Synergistic effects on transcription of DkEGase1 that involved DkNAC9 and the previously reported DkERF8/16 were identified. Electrophoretic mobility shift assay indicated that DkNAC9 could physically bind to the DkEGase1 promoter. Bimolecular fluorescence complementation and firefly luciferase complementation imaging assays indicated protein-protein interactions between DkNAC9 and DkERF8/16. Based on these findings, we conclude that DkNAC9 is a direct transcriptional activator of DkEGase1 that can co-operate with DkERF8/16 to enhance fruit post-deastringency softening.

The persimmon genome reveals clues to the evolution of a lineage-specific sex determination system in plants

Most angiosperms bear hermaphroditic flowers, but a few species have evolved outcrossing strategies, such as dioecy, the presence of separate male and female individuals. We previously investigated the mechanisms underlying dioecy in diploid persimmon (D. lotus) and found that male flowers are specified by repression of the autosomal gene MeGI by its paralog, the Y-encoded pseudo-gene OGI. This mechanism is thought to be lineage-specific, but its evolutionary path remains unknown. Here, we developed a full draft of the diploid persimmon genome (D. lotus), which revealed a lineage-specific whole-genome duplication event and provided information on the architecture of the Y chromosome. We also identified three paralogs, MeGI, OGI and newly identified Sister of MeGI (SiMeGI). Evolutionary analysis suggested that MeGI underwent adaptive evolution after the whole-genome duplication event. Transformation of tobacco plants with MeGI and SiMeGI revealed that MeGI specifically acquired a new function as a repressor of male organ development, while SiMeGI presumably maintained the original function. Later, a segmental duplication event spawned MeGI’s regulator OGI on the Y-chromosome, completing the path leading to dioecy, and probably initiating the formation of the Y-chromosome. These findings exemplify how duplication events can provide flexible genetic material available to help respond to varying environments and provide interesting parallels for our understanding of the mechanisms underlying the transition into dieocy in plants.

Plant sexuality has fascinated scientists for decades. Most plants can self-reproduce but not all. For example, a small subset of species have evolved a system called dioecy, with separate male and female individuals. Dioecy has evolved multiple times independently and, while we do not understand the molecular mechanisms underlying dioecy in many of these species yet, a picture is starting to emerge with recent progress in several dioecious species. Here, we focused on the evolutionary events leading to dioecy in persimmon. Our previous work had identified a pair of genes regulating sex in this species, called OGI and MeGI. We drafted the whole genome sequence of diploid persimmon to investigate their evolutionary history. We discovered a lineage-specific whole-genome duplication event, and observed that MeGI underwent adaptive evolution after this event. Transgenic analyses validated that MeGI newly acquired a male-suppressor function, while the other copy of this gene, SiMeGI, did not. The regulator of MeGI, OGI, resulted from a second smaller-scale segmental duplication event, finalizing the system. This study sheds light on the role of duplication as a mechanism that promote flexible genes functions, and how it can affect important biological functions, such as the establishment of a new sexual system.

Epigenetic Flexibility Underlies Somaclonal Sex Conversions in Hexaploid Persimmon

Epigenetic regulation adds a flexible layer to genetic variations, potentially enabling long-term, but reversible, changes to a trait, while maintaining genetic information. In the hexaploid Oriental persimmon (Diospyros kaki), genetically monoecious cultivars bearing male flowers require the Y-encoded small RNA (smRNA) gene, OGI. This gene represses the expression of its autosomal counterpart gene, MeGI, as part of the canonical male production system. However, a D. kaki cultivar, Saijo, which lacks the OGI gene and originally bears only female flowers, occasionally produces somaclonal mutant male and revertant female (RF) branches. In this study, we investigated the mechanisms underlying these somaclonal sex conversions in persimmon. Specifically, we aimed to unravel how a genetically female tree without the OGI gene can produce male flowers and RF flowers. Applying multi-omics approaches, we revealed that this noncanonical male production system is basically consistent with the canonical system, in which the accumulation of smRNA targeting MeGI and the considerable DNA methylation of MeGI are involved. The epigenetic status of MeGI on CGN and CHG was synchronized to the genome-wide methylation patterns, both in transition to and from the male production system. These results suggest that the somaclonal sex conversions in persimmon are driven by the genome-wide epigenetic regulatory activities. Moreover, flexibility in the epigenetic layers of long-lived plant species (e.g. trees) is important for overcoming genetic robustness.

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.

DkWRKY interacts with pyruvate kinase gene DkPK1 and promotes natural deastringency in C-PCNA persimmon

PAs, also known as condensed tannins, cause the astringency sensation in the persimmon fruit. The astringency of Chinese pollination-constant non-astringent (C-PCNA) persimmon (Diospyros kaki Thunb.) can be naturally removed on the tree, but the regulatory mechanisms of deastringency remain to be elucidated. In our previous research, DkPK1 was shown to be involved in the natural loss of astringency of C-PCNA persimmon fruit. In the present study, yeast one-hybrid (Y1H) library screening using the DkPK1 promoter as baits identified two DkWRKY transcription factor genes (DkWRKY3 and -15). The transcript levels of both DkWRKY3 and -15 exhibited a positive correlation with the decrease in soluble proanthocyanidin (PA) content during the last developmental stage in C-PCNA persimmon. Multiple sequence analysis and subcellular localization confirmed that DkWRKY3 and -15 belonging to the group II and I families, respectively, were both located in the nucleus. Dual-luciferase and Y1H assays demonstrated that DkWRKY3 and -15 can transactivate the DkPK1 promoters. The combination of DkWRKY3 and -15 most likely produced an additive activation effect compared to a single activator on DkPK1, although the two transcriptional activators were not capable of interacting. Notably, DkWRKY3 and -15 showed ubiquitous expression in various organs and abundant upregulation in seeds. Furthermore, transient overexpression of both DkWRKY3 and -15 in persimmon leaves led to a significant decrease in the content of soluble PAs but a significant increase in the expression levels of the acetaldehyde metabolism-related DkPK, DkPDC and DkADH genes. Thus, we suggest that DkWRKY3 and -15 are the upstream regulators of DkPK1 and positively regulate the natural deastringency in C-PCNA persimmon.

C2H2-Type Zinc Finger Proteins (DkZF1/2) Synergistically Control Persimmon Fruit Deastringency

Hypoxic environments are generally undesirable for most plants, but for astringent persimmon, high CO₂ treatment (CO₂ > 90%), also termed artificial high-CO₂ atmosphere (AHCA), causes acetaldehyde accumulation and precipitation of soluble tannins and could remove astringency. The multiple transcriptional regulatory linkages involved in persimmon fruit deastringency have been advanced significantly by characterizing the ethylene response factors (ERFs), WRKY and MYB; however, the involvement of zinc finger proteins for deastringency has not been investigated. In this study, five genes encoding C2H2-type zinc finger proteins were isolated and designed as DkZF1-5. Phylogenetic and sequence analyses suggested the five DkZFs could be clustered into two different subgroups. qPCR analysis indicated that transcript abundances of DkZF1/4 were significantly upregulated during AHCA treatment (1% O₂ and 95% CO₂) at day 1, DkZF2/5 at both day 1 and 2, while DkZF3 at day 2. Dual-luciferase assay indicated DkZF1 and DkZF2 as the activators of deastringency-related structural genes (DkPDC2 and DkADH1) and transcription factors (DkERF9/10). Moreover, combinative effects between various transcription factors were investigated, indicating that DkZF1 and DkZF2 synergistically showed significantly stronger activations on the DkPDC2 promoter. Further, both bimolecular fluorescence complementation (BiFC) and yeast two hybrid (Y2H) assays confirmed that DkZF2 had protein–protein interactions with DkZF1. Thus, these findings illustrate the regulatory mechanisms of zinc finger proteins for persimmon fruit deastringency under AHCA.

Evolution of Lineage-Specific Gene Networks Underlying the Considerable Fruit Shape Diversity in Persimmon

The shapes of plant organs reflect the evolution of each lineage and have been diversified according to lineage-specific adaptations to environment. Research on the molecular pathways responsible for organ shapes has traditionally been focused mainly on leaves or flowers. Thus, little is known about the pathways controlling fruit shapes, despite their diversity in some plant species. In this study, we analyzed oriental persimmon (Diospyros kaki), which exhibits considerable diversity in fruit shapes among cultivars, to elucidate the underlying molecular mechanism using transcriptomic data and quantitative evaluation. First, to filter the candidate genes associated with persimmon fruit shapes, the whole gene expression patterns obtained using mRNA-Seq analysis from 100 individuals, including a segregated population and various cultivars, were assessed to detect correlations with principal component scores for fruit shapes characterized with elliptic Fourier descriptors. Next, a gene co-expression network analysis with weighted gene co-expression network analysis (WGCNA) package revealed that class 1 KNOX family genes and SEEDSTICK function as integrators along with some phytohormone-related genes, to regulate the fruit shape diversity. On the other hand, the OVATE family genes also contribute to fruit shape diversity, of which pathway would be potentially shared with other plant species. Evolutionary aspects suggest that acquisition of a high lineage-specific and variable expression of class 1 KNOX gene, knotted-like homeobox of Arabidopsis thaliana 1 (KNAT1), in young fruit is important for establishing the persimmon-specific mechanism that determines fruit shape diversity.

Comparative transcriptomic analysis reveals genetic divergence and domestication genes in Diospyros

BACKGROUND

Persimmon (Diospyros kaki) is the most economically cultivated species belonging to the genus Diospyros. However, little is known about the interspecific diversity and mechanism of domestication, partly due to the lack of genomic information that is available for closely related species of D. kaki (DK). Here, we performed transcriptome sequencing on nine samples, including DK, a variety of DK and seven closely related species, to evaluate the interspecific genetic divergence and to identify candidate genes involved in persimmon domestication.

RESULTS

We obtained a total of 483,421 unigenes with N50 at 1490 bp in the nine Diospyros samples and identified 2603 orthogroups that were shared among all the samples using OrthoMCL analysis. A phylogenetic tree was established based on the tandem 2603 one-to-one single copy gene alignments, showing that DK was closely related to D. kaki var. silvestris (DKV) and that it clustered with the clade of D. deyangnsis (DD) and was farthest from the D. cathayensis (DC) species. The nonsynonymous substitutions (Ka), via synonymous substitution (Ks) ratios, was directly proportional to the genetic relationship of the different species. The higher the Ka/Ks ratios, the longer the distance was. Moreover, 31 positively selected genes (PSGs) involved in carbohydrate metabolism and phenolic metabolism were identified and isolated, and nearly all PSGs except the MATE gene had a high expression in the DK or DKV species. It was hypothesized that these genes might contribute to the domestication of the DK species. Finally, we developed the expressed sequence tag-simple sequence repeat (EST-SSR) and identified 2 unique amplicons DKSSR10 and DKSSR39: the former was absent in the DC species but was present in the other species, the latter had a long amplification product in the DJ species.

CONCLUSION

This study presents the first transcriptome resources for the closely related species of persimmon and reveals interspecific genetic divergence. It is speculated that DK is derived from the hybridization of DD and DO species. Furthermore, our analysis suggests candidate PSGs that may be crucial for the adaptation, domestication, and speciation of persimmon relatives and suggests that DKSSR10 and DKSSSR39 could potentially serve as species-specific molecular markers.

High-CO2/Hypoxia-Responsive Transcription Factors DkERF24 and DkWRKY1 Interact and Activate DkPDC2 Promoter

Identification and functional characterization of hypoxia-responsive transcription factors is important for understanding plant responses to natural anaerobic environments and during storage and transport of fresh horticultural products. In this study, yeast one-hybrid library screening using the persimmon (Diospyros kaki) pyruvate decarboxylase (DkPDC2) promoter identified three ethylene response factor (ERF) genes (DkERF23/DkERF24/DkERF25) and four WRKY transcription factor genes (DkWRKY/DdkWRKY5/DkWRKY6/DkWRKY7) that were differentially expressed in response to high CO2 (95%, with 4% N2 and 1% oxygen) and high N2 (99% N2 and 1% oxygen). Yeast one-hybrid assays and electrophoretic mobility shift assays indicated that DkERF23, DkERF24, DkERF25, DkWRKY6, and DkWRKY7 could directly bind to the DkPDC2 promoter. Dual-luciferase assays confirmed that these transcription factors were capable of transactivating the DkPDC2 promoter. DkERF24 and DkWRKY1 in combination synergistically transactivated the DkPDC2 promoter, and yeast two-hybrid and bimolecular fluorescence complementation assays confirmed protein-protein interaction between DkERF24 and DkWRKY1. Transient overexpression of DkERF24 and DkWRKY1 separately and in combination in persimmon fruit discs was effective in maintaining insolubilization of tannins, concomitantly with the accumulation of DkPDC2 transcripts. Studies with Arabidopsis (Arabidopsis thaliana) homologs AtERF1 and AtWRKY53 indicated that similar protein-protein interactions and synergistic regulatory effects also occur with the DkPDC2 promoter. We propose that an ERF and WRKY transcription factor complex contributes to responses to hypoxia in both persimmon fruit and Arabidopsis, and the possibility that this is a general plant response requires further investigation.

Gene networks orchestrated by MeGI: a single-factor mechanism underlying sex determination in persimmon

Separating male and female sex organs is one of the main strategies used to maintain genetic diversity within a species. However, the genetic determinants and their regulatory mechanisms have been identified in only a few species. In dioecious persimmons, the homeodomain transcription factor, MeGI, which is the target of a Y chromosome-encoded small-RNA, OGI, can determine floral sexuality. The basic features of this system are conserved in the monoecious hexaploid Oriental persimmon, in which an additional epigenetic regulation of MeGI determines floral sexuality. The downstream regulatory pathways of MeGI remain uncharacterized. In this study, we examined transcriptomic data for male and female flowers from monoecious persimmon cultivars to unveil the gene networks orchestrated by MeGI. A network visualization and cistrome assessment suggested that class-1 KNOTTED-like homeobox (KNOX)/ovate family protein (OFP)/growth regulating factors (GRFs) and short vegetative phase (SVP) genes mediate the differences in gynoecium and androecium development between male and female flowers, respectively. The expression of these genes is directly controlled by MeGI. The gene networks also suggested that some cytokinin, auxin, and gibberellin signaling genes function cooperatively in the KNOX/OFP/GRF pathway during gynoecium differentiation. Meanwhile, SVP may repress PI expression in developing androecia. Overall, our results suggest that MeGI evolved the ability to promote gynoecium development and suppress androecium development by regulating KNOX/OFP/GRF and SVP expression levels, respectively. These insights may help to clarify the molecular mechanism underlying the production of unisexual flowers, while also elucidating the physiological background enabling a single-factor system to establish dioecy in plants.

Interspecific chloroplast genome sequence diversity and genomic resources in Diospyros

BACKGROUND

Fruits of persimmon plants are traditional healthy food in China, Korea and Japan. However, due to the shortage of morphological and DNA markers, the development of persimmon industry has been heavily inhibited.

RESULTS

Chloroplast genomes of Diospyros cathayensis, D. virginiana, D. rhombifolia and D. deyangensis were newly sequenced. Comparative analyses of ten chloroplast genomes including six previously published chloroplast genomes of Diospyros provided new insights into the genome sequence diversity and genomic resources of the genus. Eight hyper-variable regions, trnH-psbA, rps16-trnQ, rpoB-trnC, rps4-trnT-trnL, ndhF, ndhF-rpl32-trnL, ycf1a, and ycf1b, were discovered and can be used as chloroplast DNA markers at/above species levels. The complete chloroplast genome sequences provided the best resolution at inter-specific level in comparison with different chloroplast DNA sequence datasets.

CONCLUSION

Diospyros oleifera, D. deyangensis, D. virginiana, D. glaucifolia, D. lotus and D. jinzaoshi are important wild species closely related to the cultivated persimmon D. kaki. The hyper-variable regions can be used as DNA markers for global genetic diversity detection of Diospyros. Deeper study on these taxa would be helpful for elucidating the origin of D. kaki.

Functional characterization of persimmon β-galactosidase gene DkGAL1 in tomato reveals cell wall modification related to fruit ripening and radicle elongation

Cell wall metabolism during fruit ripening is a highly organized process that involves complex interplay among various cell wall hydrolases. Among these cell wall hydrolases, β-galactosidase has been identified to participate in cell wall metabolism via its ability to catalyze galactosyl metabolism from the large and complex side chains of cell walls. In this study, the galactose content in the pericarp increased during persimmon fruit ripening, but cell wall galactosyl residues decreased, indicating a relationship between galactose metabolism and persimmon fruit ripening. Expression of a previously isolated β-galactosidase gene, DkGAL1, increased 25.01-fold during fruit ripening. Heterologous expression of DkGAL1 under the CaMV 35S promoter in tomato accelerated on-plant and postharvest fruits ripening. The fruit firmness of one of transgenic line, OE-18, was 23.83% lower than that of WT at the breaker stage. The transgenic fruits produced more ethylene by promoting the expression of ethylene synthesis-related genes and cell wall degradation-related genes. Overexpression of DkGAL1 in tomato also reduced cell-to-cell adhesion and promoted both wider intercellular spaces and less cell compaction in transgenic fruit structures. Moreover, DkGAL1 was involved in seed germination and radicle elongation in transgenic tomato seeds. These results confirm the role of DkGAL1 in fruit ripening and suggest that this gene alters galactose metabolism in the fruit, which can promote ripening and reduce cellular adhesion. In addition, the role of DkGAL1 is not limited to fruit softening; DkGAL1 was also involved in seed germination and radicle elongation in transgenic tomato seeds.

A transcription factor network responsive to high CO2/hypoxia is involved in deastringency in persimmon fruit

Plant responses to anaerobic environments are regulated by ethylene-response factors (ERFs) in both vegetative and productive organs, but the roles of other transcription factors (TFs) in hypoxia responses are poorly understood. In this study, eight TFs (DkbHLH1, DkMYB9/10/11, DkRH2-1, DkGT3-1, DkAN1-1, DkHSF1) were shown to be strongly up-regulated by an artificial high-CO2 atmosphere (1% O2 and 95% CO2). Dual-luciferase assays indicated that some TFs were activators of previously characterized DkERFs, including DkMYB10 for the DkERF9 promoter, DkERF18/19 and DkMYB6 for the DkERF19 promoter, and DkERF21/22 for the DkERF10 promoter. Yeast one-hybrid and cis-element mutagenesis confirmed these physical interactions with one exception. The potential roles of these TFs in persimmon fruit deastringency were analysed by investigating their transient over-expression (TOX) in persimmon fruit discs, which indicated that DkMYB6TOX, DkMYB10TOX, DkERF18TOX, and DkERF19TOX were all effective in causing insolubilization of tannins, concomitantly with the up-regulation of the corresponding genes. These results indicated that multiple TFs of different classes are responsive to high-CO2/hypoxia in fruit tissues, and that a TF-TF regulatory cascade is involved in the hypoxia responses involving the Group VII DkERF10, and DkERFs and DkMYBs.

Role of Brassinosteroids in Persimmon ( Diospyros kaki L.) Fruit Ripening

Brassinosteroids (BRs) are phytohormones that regulate numerous processes including fruit ripening. In this study, persimmon ( Diospyros kaki L.) fruits were treated with 24-epibrassinolide (EBR) or brassinazole (Brz, a BR biosynthesis inhibitor) and then stored at ambient temperature. The results show that endogenous BR contents gradually increased during persimmon fruit ripening. EBR treatment significantly increased both the content of water-soluble pectin and the activities of polygalacturonase, pectate lyase, and endo-1,4-beta-glucanase but significantly reduced the content of acid-soluble pectin and cellulose, resulting in rapid fruit softening. The EBR treatment also promoted ethylene production and respiration rate. In contrast, Brz treatment delayed persimmon fruit ripening. qRT-PCR analysis showed that DkPG1, DkPL1, DkPE2, DkEGase1, DkACO2, DkACS1, and DkACS2 were up-regulated (especially a 38-fold increase in DkEGase1) in the fruit of the EBR-treated group. These results suggest that BRs are involved in persimmon fruit ripening by influencing cell-wall-degrading enzymes and ethylene biosynthesis.

DkNAC7, a novel high-CO2/hypoxia-induced NAC transcription factor, regulates persimmon fruit de-astringency

Artificial high-CO₂ atmosphere (AHCA, 95% CO₂ and 1% O₂) has been widely applied as a postharvest de-astringency treatment for persimmon fruit. AHCA increases expression of transcription factors, including ethylene response factors (DkERF), that target de-astringency genes. Here, the promoter of DkERF9, a previously characterized AHCA-inducible and de-astringency regulator, was utilized to screen a cDNA library by yeast one hybrid assay. A novel NAC transcription factor, named DkNAC7, was identified. Dual-luciferase assay indicated that DkNAC7 could not only trans-activate the promoter of DkERF9, but also activated the previously identified deastringency-related gene DkPDC2. Real-time PCR analysis showed that DkNAC7 was up-regulated by AHCA treatment, in concert with the removal of astringency from persimmon fruit and subcellular localization showed DkNAC7 was located in the nucleus. Thus, these results indicate that DkNAC7 is a putative transcriptional activator involved in regulating persimmon fruit deastringency by trans-activition on both DkERF9 and DkPDC2, which encodes pyruvate decarboxylase.

Characterization of a gene regulatory network underlying astringency loss in persimmon fruit

Transcriptome analysis of a persimmon population segregating for an astringency trait in fruit suggested central roles for a limited number of transcriptional regulators in the loss of proanthocyanidin accumulation. Persimmon (Diospyros kaki; 2n = 6x = 90) accumulates a large amount of proanthocyanidins (PAs) in its fruit, resulting in an astringent taste. Persimmon cultivars are classified into four types based on the nature of astringency loss and the amount of PAs at maturity. Pollination constant and non-astringent (PCNA)-type cultivars stop accumulating PAs in the early stages of fruit development and their fruit can be consumed when still firm without the need for artificial deastringency treatments. While the PCNA trait has been shown to be conferred by a recessive allele at a single locus (ASTRINGENCY; AST), the exact genetic determinant remains unidentified. Here, we conducted transcriptome analyses to elucidate the regulatory mechanism underlying this trait using developing fruits of an F₁ population segregating for the PCNA trait. Comparisons of the transcriptomes of PCNA and non-PCNA individuals and hierarchical clustering revealed that genes related to the flavonoid pathway and to abiotic stress responses involving light stimulation were expressed coordinately with PA accumulation. Furthermore, coexpression network analyses suggested that three putative transcription factors were central to the PA regulatory network and that at least DkMYB4 and/or DkMYC1, which have been reported to form a protein complex with each other for PA regulation, may have a central role in the differential expression of PA biosynthetic pathway genes between PCNA and non-PCNA.

Hypoxia-responsive ERFs involved in postdeastringency softening of persimmon fruit

Removal of astringency by endogenously formed acetaldehyde, achieved by postharvest anaerobic treatment, is of critical importance for many types of persimmon fruit. Although an anaerobic environment accelerates de-astringency, it also has the deleterious effect of promoting excessive softening, reducing shelf life and marketability. Some hypoxia-responsive ethylene response factors (ERFs) participate in anaerobic de-astringency, but their role in accelerated softening was unclear. Undesirable rapid softening induced by high CO₂ (95%) was ameliorated by adding the ethylene inhibitor 1-MCP (1 μL/L), resulting in reduced astringency while maintaining firmness, suggesting that CO₂ -induced softening involves ethylene signalling. Among the hypoxia-responsive genes, expression of eight involved in fruit cell wall metabolism (Dkβ-gal1/4, DkEGase1, DkPE1/2, DkPG1, DkXTH9/10) and three ethylene response factor genes (DkERF8/16/19) showed significant correlations with postdeastringency fruit softening. Dual-luciferase assay indicated that DkERF8/16/19 could trans-activate the DkXTH9 promoter and this interaction was abolished by a mutation introduced into the C-repeat/dehydration-responsive element of the DkXTH9 promoter, supporting the conclusion that these DkERFs bind directly to the DkXTH9 promoter and regulate this gene, which encodes an important cell wall metabolism enzyme. Some hypoxia-responsive ERF genes are involved in deastringency and softening, and this linkage was uncoupled by 1-MCP. Fruit of the Japanese cultivar 'Tonewase' provide a model for altered anaerobic response, as they lost astringency yet maintained firmness after CO₂ treatment without 1-MCP and changes in cell wall enzymes and ERFs did not occur.

Overexpression of persimmon DkXTH1 enhanced tolerance to abiotic stress and delayed fruit softening in transgenic plants

DkXTH1 promoted cell elongation and more strength to maintain structural integrity by involving in cell wall assembly, thus enhanced tolerance to abiotic stress with broader phenotype in transgenic plants. Xyloglucan endotransglucosylase/hydrolase (XTH) is thought to play a key role in cell wall modifications by cleaving and re-joining xyloglucan, and participates in the diverse physiological processes. DkXTH1 was found to peak in immature expanding persimmon fruit, and its higher expression level exhibited along with firmer fruit during storage. In the present study, transgenic Arabidopsis and tomato plants were generated with DkXTH1 constitutively expressed. Overexpression of DkXTH1 enhanced tolerance to salt, ABA and drought stresses in transgenic Arabidopsis plants with respect to root and leaf growth, and survival. Transgenic tomatoes collected at the mature green stage, presented delayed fruit softening coupled with postponed color change, a later and lower ethylene peak, and higher firmness in comparison with the wild-type tomatoes during storage. Furthermore, broader leaves and tomato fruit with larger diameter were gained in transgenic Arabidopsis and tomato, respectively. Most importantly, transgenic plants exhibited more large and irregular cells with higher density of cell wall and intercellular spaces, resulting from the overactivity of XET enzymes involving in cell wall assembly. We suggest that DkXTH1 expression resulted in cells with more strength and thickness to maintain structural integrity, and thus enhanced tolerance to abiotic stress and delayed fruit softening in transgenic plants.

DkXTH8, a novel xyloglucan endotransglucosylase/hydrolase in persimmon, alters cell wall structure and promotes leaf senescence and fruit postharvest softening

Fruit softening is mainly associated with cell wall structural modifications, and members of the xyloglucan endotransglucosylase/hydrolase (XTH) family are key enzymes involved in cleaving and re-joining xyloglucan in the cell wall. In this work, we isolated a new XTH gene, DkXTH8, from persimmon fruit. Transcriptional profiling revealed that DkXTH8 peaked during dramatic fruit softening, and expression of DkXTH8 was stimulated by propylene and abscisic acid but suppressed by gibberellic acid and 1-MCP. Transient expression assays in onion epidermal cells indicated direct localization of DkXTH8 to the cell wall via its signal peptide. When expressed in vitro, the recombinant DkXTH8 protein exhibited strict xyloglucan endotransglycosylase activity, whereas no xyloglucan endohydrolase activity was observed. Furthermore, overexpression of DkXTH8 resulted in increased leaf senescence coupled with higher electrolyte leakage in Arabidopsis and faster fruit ripening and softening rates in tomato. Most importantly, transgenic plants overexpressing DkXTH8 displayed more irregular and twisted cells due to cell wall restructuring, resulting in wider interstitial spaces with less compact cells. We suggest that DkXTH8 expression causes cells to be easily destroyed, increases membrane permeability and cell peroxidation, and accelerates leaf senescence and fruit softening in transgenic plants.

Epigenetic Regulation of the Sex Determination Gene MeGI in Polyploid Persimmon

Epigenetic regulation can add a flexible layer to genetic variation, potentially enabling long-term but reversible cis-regulatory changes to an allele while maintaining its DNA sequence. Here, we present a case in which alternative epigenetic states lead to reversible sex determination in the hexaploid persimmon Diospyros kaki Previously, we elucidated the molecular mechanism of sex determination in diploid persimmon and demonstrated the action of a Y-encoded sex determinant pseudogene called OGI, which produces small RNAs targeting the autosomal gene MeGI, resulting in separate male and female individuals (dioecy). We contrast these findings with the discovery, in hexaploid persimmon, of an additional layer of regulation in the form of DNA methylation of the MeGI promoter associated with the production of both male and female flowers in genetically male trees. Consistent with this model, developing male buds exhibited higher methylation levels across the MeGI promoter than developing female flowers from either monoecious or female trees. Additionally, a DNA methylation inhibitor induced developing male buds to form feminized flowers. Concurrently, in Y-chromosome-carrying trees, the expression of OGI is silenced by the presence of a SINE (short interspersed nuclear element)-like insertion in the OGI promoter. Our findings provide an example of an adaptive scenario involving epigenetic plasticity.

Molecular cloning and functional characterization of DkMATE1 involved in proanthocyanidin precursor transport in persimmon (Diospyros kaki Thunb.) fruit

Persimmon fruits accumulate a large amount of proanthocyanidins (PAs) in "tannin cells" during development that cause the sensation of astringency due to coagulation of oral proteins. Pollination-constant non-astringent (PCNA) is a spontaneous mutant persimmon phenotype that loses its astringency naturally on the tree at maturity; while the more common non-PCNA fruits remain rich in PAs until they are fully ripened. Here, we isolated a DkMATE1 gene encoding a Multidrug And Toxic Compound Extrusion (MATE) family protein from the Chinese PCNA (C-PCNA) 'Eshi 1'. Expression patterns of DkMATE1 were positively correlated with the accumulation of PAs in different types of persimmons fruits during fruit development. An analysis of the inferred amino acid sequences and phylogenetic relationships indicated that DkMATE1 is a putative PA precursor transporter, and subcellular localization assays revealed that DkMATE1 is localized in the vacuolar membrane. Ectopic expression of the DkMATE1 in Arabidopsis tt12 mutant supported that DkMATE1 could complement its biological function in transporting epicatechin 3'-O-glucoside as a PAs precursor from the cytoplasm to vacuole. Furthermore, the transient over-expression and silencing of DkMATE1 in 'Mopanshi' persimmon leaves resulted in a significant increase and a decrease in PA content, respectively. The analysis of cis-elements in DkMATE1 promoter regions indicated that DkMATE1 might be regulated by DkMYB4, another well-known structural gene in persimmon. Overall, our results show that DkMATE1 may be an essential PA precursor membrane transporter that plays an important role in PA biosynthesis in persimmon.

Selection of Suitable Reference Genes for RT-qPCR Normalization under Abiotic Stresses and Hormone Stimulation in Persimmon (Diospyros kaki Thunb)

The success of quantitative real-time reverse transcription polymerase chain reaction (RT-qPCR) to quantify gene expression depends on the stability of the reference genes used for data normalization. To date, systematic screening for reference genes in persimmon (Diospyros kaki Thunb) has never been reported. In this study, 13 candidate reference genes were cloned from 'Nantongxiaofangshi' using information available in the transcriptome database. Their expression stability was assessed by geNorm and NormFinder algorithms under abiotic stress and hormone stimulation. Our results showed that the most suitable reference genes across all samples were UBC and GAPDH, and not the commonly used persimmon reference gene ACT. In addition, UBC combined with RPII or TUA were found to be appropriate for the "abiotic stress" group and α-TUB combined with PP2A were found to be appropriate for the "hormone stimuli" group. For further validation, the transcript level of the DkDREB2C homologue under heat stress was studied with the selected genes (CYP, GAPDH, TUA, UBC, α-TUB, and EF1-α). The results suggested that it is necessary to choose appropriate reference genes according to the test materials or experimental conditions. Our study will be useful for future studies on gene expression in persimmon.

Five Complete Chloroplast Genome Sequences from Diospyros: Genome Organization and Comparative Analysis

Diospyros is the largest genus in Ebenaceae, comprising more than 500 species with remarkable economic value, especially Diospyros kaki Thunb., which has traditionally been an important food resource in China, Korea, and Japan. Complete chloroplast (cp) genomes from D. kaki, D. lotus L., D. oleifera Cheng., D. glaucifolia Metc., and Diospyros 'Jinzaoshi' were sequenced using Illumina sequencing technology. This is the first cp genome reported in Ebenaceae. The cp genome sequences of Diospyros ranged from 157,300 to 157,784 bp in length, presenting a typical quadripartite structure with two inverted repeats each separated by one large and one small single-copy region. For each cp genome, 134 genes were annotated, including 80 protein-coding, 31 tRNA, and 4 rRNA unique genes. In all, 179 repeats and 283 single sequence repeats were identified. Four hypervariable regions, namely, intergenic region of trnQ_rps16, trnV_ndhC, and psbD_trnT, and intron of ndhA, were identified in the Diospyros genomes. Phylogenetic analyses based on the whole cp genome, protein-coding, and intergenic and intron sequences indicated that D. oleifera is closely related to D. kaki and could be used as a model plant for future research on D. kaki; to our knowledge, this is proposed for the first time. Further, these analyses together with two large deletions (301 and 140 bp) in the cp genome of D. 'Jinzaoshi', support its placement as a new species in Diospyros. Both maximum parsimony and likelihood analyses for 19 taxa indicated the basal position of Ericales in asterids and suggested that Ebenaceae is monophyletic in Ericales.

Biochemical Constituents and in Vitro Antioxidant and Anticholinesterase Potential of Seeds from Native Korean Persimmon Genotypes

In the current study, the functional and biochemical potential of the seeds of four persimmon cultivars (PC1, PC2, PC3 and PC4) and their role against oxidative stress and acetylcholinesterase (AChE) inhibition were evaluated. In terms of biochemical compositions, free amino acids, fatty acids and organic acids analysis was performed. The free amino acids ranged from 2617.31 (PC2) to 3773.01 μg∙g(-1) dry weight (PC4). Oleic acid and linoleic acid were the principal fatty acids, which were significantly higher in PC4 and PC1, respectively. PC4 presented the highest amount of organic acid content (4212 mg∙kg(-1)), whereas PC2 presented the lowest (2498 mg∙kg(-1)). PC2 contained higher total phenolic content and flavonoid content, whereas PC3 had the lowest amount as compared to other cultivars. The in vitro DPPH, ABTS and superoxide anion radicals scavenging activity increased in a dose-dependent manner, whereas PC2 showed significantly higher scavenging activities as compared to PC1, PC2 and PC4 types. In the case of AChE inhibition, PC4 showed a moderate activity (67.34% ± 1.8%). In conclusion, the current findings reveal that the studied persimmon seeds cultivars are a source of bioactive natural antioxidants and AChE inhibitors. Such natural products could be employed in pharmaceutical and food industries, whilst can also be considered for the treatment of neurodegenerative diseases such as Alzheimer's.

Sequencing of whole plastid genomes and nuclear ribosomal DNA of Diospyros species (Ebenaceae) endemic to New Caledonia: many species, little divergence

BACKGROUND AND AIMS

Some plant groups, especially on islands, have been shaped by strong ancestral bottlenecks and rapid, recent radiation of phenotypic characters. Single molecular markers are often not informative enough for phylogenetic reconstruction in such plant groups. Whole plastid genomes and nuclear ribosomal DNA (nrDNA) are viewed by many researchers as sources of information for phylogenetic reconstruction of groups in which expected levels of divergence in standard markers are low. Here we evaluate the usefulness of these data types to resolve phylogenetic relationships among closely related Diospyros species.

METHODS

Twenty-two closely related Diospyros species from New Caledonia were investigated using whole plastid genomes and nrDNA data from low-coverage next-generation sequencing (NGS). Phylogenetic trees were inferred using maximum parsimony, maximum likelihood and Bayesian inference on separate plastid and nrDNA and combined matrices.

KEY RESULTS

The plastid and nrDNA sequences were, singly and together, unable to provide well supported phylogenetic relationships among the closely related New Caledonian Diospyros species. In the nrDNA, a 6-fold greater percentage of parsimony-informative characters compared with plastid DNA was found, but the total number of informative sites was greater for the much larger plastid DNA genomes. Combining the plastid and nuclear data improved resolution. Plastid results showed a trend towards geographical clustering of accessions rather than following taxonomic species.

CONCLUSIONS

In plant groups in which multiple plastid markers are not sufficiently informative, an investigation at the level of the entire plastid genome may also not be sufficient for detailed phylogenetic reconstruction. Sequencing of complete plastid genomes and nrDNA repeats seems to clarify some relationships among the New Caledonian Diospyros species, but the higher percentage of parsimony-informative characters in nrDNA compared with plastid DNA did not help to resolve the phylogenetic tree because the total number of variable sites was much lower than in the entire plastid genome. The geographical clustering of the individuals against a background of overall low sequence divergence could indicate transfer of plastid genomes due to hybridization and introgression following secondary contact.

Involvement of DkTGA1 Transcription Factor in Anaerobic Response Leading to Persimmon Fruit Postharvest De-Astringency

Persimmon fruit are unique in accumulating proanthocyanidins (tannins) during development, which cause astringency in mature fruit. In ‘Mopanshi’ persimmon, astringency can be removed by treatment with 95% CO₂, which increases the concentrations of ethanol and acetaldehyde by glycolysis, and precipitates the soluble tannin. A TGA transcription factor, DkTGA1, belonging to the bZIP super family, was isolated from an RNA-seq database and real-time quantitative PCR indicated that DkTGA1 was up-regulated by CO₂ treatment, in concert with the removal of astringency from persimmon fruit. Dual-luciferase assay revealed that DkTGA1 had a small (less than 2-fold), but significant effect on the promoters of de-astringency-related genes DkADH1, DkPDC2 and DkPDC3, which encode enzymes catalyzing formation of acetaldehyde and ethanol. A combination of DkTGA1 and a second transcription factor, DkERF9, shown previously to be related to de-astringency, showed additive effects on the activation of the DkPDC2 promoter. Yeast one-hybrid assay showed that DkERF9, but not DkTGA1, could bind to the DkPDC2 promoter. Thus, although DkTGA1 expression is positively associated with persimmon fruit de-astringency, trans-activation analyses with DkPDC2 indicates it is likely to act by binding indirectly DkPDC2 promoter, might with helps of DkERF9.

Investigation and Analysis of Genetic Diversity of Diospyros Germplasms Using SCoT Molecular Markers in Guangxi

BACKGROUND

Knowledge about genetic diversity and relationships among germplasms could be an invaluable aid in diospyros improvement strategies.

METHODS

This study was designed to analyze the genetic diversity and relationship of local and natural varieties in Guangxi Zhuang Autonomous Region of China using start codon targeted polymorphism (SCoT) markers. The accessions of 95 diospyros germplasms belonging to four species Diospyros kaki Thunb, D. oleifera Cheng, D. kaki var. silverstris Mak, and D. lotus Linn were collected from different eco-climatic zones in Guangxi and were analyzed using SCoT markers.

RESULTS

Results indicated that the accessions of 95 diospyros germplasms could be distinguished using SCoT markers, and were divided into three groups at similarity coefficient of 0.608; these germplasms that belong to the same species were clustered together; of these, the degree of genetic diversity of the natural D. kaki var. silverstris Mak population was richest among the four species; the geographical distance showed that the 12 natural populations of D. kaki var. silverstris Mak were divided into two groups at similarity coefficient of 0.19. Meanwhile, in order to further verify the stable and useful of SCoT markers in diospyros germplasms, SSR markers were also used in current research to analyze the genetic diversity and relationship in the same diospyros germplasms. Once again, majority of germplasms that belong to the same species were clustered together. Thus SCoT markers were stable and especially useful for analysis of the genetic diversity and relationship in diospyros germplasms.

DISCUSSION

The molecular characterization and diversity assessment of diospyros were very important for conservation of diospyros germplasm resources, meanwhile for diospyros improvement.

Analysis of xyloglucan endotransglycosylase/hydrolase (XTH) genes and diverse roles of isoenzymes during persimmon fruit development and postharvest softening

Xyloglucan endotransglycosylase/hydrolase (XTH) enzymes have played a role in the remodeling of cell wall hemicelluloses. To investigate the function of XTHs in persimmon (Diospyros kaki L.) fruit development and postharvest softening, five cDNAs (DkXTH1 to DkXTH5), whose putative proteins contained the conserved DEIDFEFLG motif of XTH, were cloned. Real time quantitative PCR analysis revealed that DkXTH1, DkXTH4, and DkXTH5 peaked in immature expanding fruit, and their higher expression was observed along with higher fruit firmness in cold-treated fruit or firmer cultivar fruit during storage. The opposite gene expression patterns were observed in DkXTH2 and DkXTH3, which reached maxima concomitance with pronounced fruit softening. Meanwhile, the xyloglucan endotransglycosylase (XET) enzymes play important roles in both the rapid growth and ripening of persimmon fruit. Furthermore, the recombined DkXTH1 and DkXTH2 proteins showed significant XET activity without any detected XEH activity. However, the XET activity of recombined DkXTH2 protein had a higher affinity for small acceptor molecules than that of recombined DkXTH1 protein. The former might prefer to participate in cell wall restructuring, and the latter is more inclined to participate in cell wall assembly. Besides, DKXTH proteins could function by targeting to the cell wall under regulation of a signal peptide. The data suggested that individual DKXTHs could exhibit different patterns of expression, and the encoded products possessed specific enzymatic properties conferring on their respective functions in growth and postharvest softening of persimmon fruit.

Identification and characterization of microRNAs from Chinese pollination constant non-astringent persimmon using high-throughput sequencing

BACKGROUND

microRNAs (miRNAs) have been shown to play key roles in regulating gene expression at post-transcriptional level, but miRNAs associated with natural deastringency of Chinese pollination-constant nonastringent persimmon (CPCNA) have never been identified.

RESULTS

In this study, two small RNA libraries established using 'Eshi No. 1' persimmon (Diospyros kaki Thunb.; CPCNA) fruits collected at 15 and 20 weeks after flowering (WAF) were sequenced through Solexa platform in order to identify miRNAs involved in deastringency of persimmon. A total of 6,258,487 and 7,634,169 reads were generated for the libraries at 15 and 20 WAF, respectively. Based on sequence similarity and hairpin structure prediction, 236 known miRNAs belonging to 65 miRNA families and 33 novel miRNAs were identified using persimmon transcriptome data. Sixty one of the characterized miRNAs exhibited pronounced difference in the expression levels between 15 and 20 WAF, 17 up-regulated and 44 down-regulated. Expression profiles of 12 conserved and 10 novel miRNAs were validated by stem loop qRT-PCR. A total of 198 target genes were predicted for the differentially expressed miRNAs, including several genes that have been reported to be implicated in proanthocyanidins (PAs, or called tannin) accumulation. In addition, two transcription factors, a GRF and a bHLH, were experimentally confirmed as the targets of dka-miR396 and dka-miR395, respectively.

CONCLUSIONS

Taken together, the present data unraveled several important miRNAs in persimmon. Among them, miR395p-3p and miR858b may regulate bHLH and MYB, respectively, which are influenced by SPL under the control of miR156j-5p and in turn regulate the structural genes involved in PA biosynthesis. In addition, dka-miR396g and miR2911a may regulate their target genes associated with glucosylation and insolubilization of tannin precursors. All of these miRNAs might play key roles in the regulation of (de)astringency in persimmon fruits under normal development conditions.

Two novel anoxia-induced ethylene response factors that interact with promoters of deastringency-related genes from persimmon

A hypoxic environment is generally undesirable for most plants and stimulates anaerobic metabolism. It is a beneficial treatment, however, for the removal of astringency from persimmon to improve the fruit quality after harvest. High soluble tannins (SCTs) content is one of most important causes of astringency. High CO2 (95%) treatment effectively reduced SCTs in both "Mopan" and "Gongcheng-shuishi" persimmon fruit by causing increases in acetaldehyde. Using RNA-seq and realtime PCR, twelve ethylene response factor genes (DkERF11-22) were isolated and characterized, to determine those responsive to high CO2 treatment. Only two genes, DkERF19 and DkERF22, showed trans-activation effects on the promoters of deastringency-related genes pyruvate decarboxylase genes (DkPDC2 and DkPDC3) and the transcript levels of these genes was enhanced by hypoxia. Moreover, DkERF19 and the previously isolated DkERF9 had additive effects on activating the DkPDC2 promoter. Taken together, these results provide further evidence that transcriptome changes in the level of DkERF mRNAs regulate deastringency-related genes and their role in the mechanism of persimmon fruit deastringency is discussed.