Cytological characteristics of blueberry fruit development

Using the blueberry cultivar "Powderblue" after pollination, fruits at different developmental stages were collected for study. The transverse and longitudinal diameters, individual fruit weight, and fruit water content were measured during their development. Employing tissue sectioning and microscopy techniques, we systematically studied the morphological features and anatomical structures of the fruits and seeds at various developmental stages, aiming to elucidate the cytological patterns during blueberry fruit development. The results of our study revealed that the "Powderblue" blueberry fruit growth and development followed a double "S" curve. Mature "Powderblue" blueberries were blue-black in color, elliptical in shape, with five locules, an inferior ovary, and an average fruit weight of 1.73 ± 0.17 g, and a moisture content of 78.865 ± 0.9%. Blueberry fruit flesh cells were densely arranged with no apparent intercellular spaces, and mesocarp cells accounted for 52.06 ± 7.4% of fruit cells. In the early fruit development stages, the fruit flesh cells were rapidly dividing, significantly increasing in number but without greatly affecting the fruit's morphological characteristics. During the later stages of fruit development, the expansion of the fruit flesh cells became prominent, resulting in a noticeable increase in the fruit's dimensions. Except for the epidermal cells, cells in all fruit tissues showed varying degrees of rupture as fruit development progressed, with the extent of cell rupture increasing, becoming increasingly apparent as the fruit gradually softened. Additionally, numerous brachysclereids (stone cells) appeared in the fruit flesh cells. Stone cells are mostly present individually in the fruit flesh tissue, while in the placental tissue, they often group together. The "Powderblue" blueberry seeds were light brown, 4.13 ± 0.42 mm long, 2.2 ± 0.14 mm wide, with each fruit containing 50-60 seeds. The "Powderblue" seeds mainly consisted of the seed coat, endosperm, and embryo. The embryo was located at the chalazal end in the center of the endosperm and was spatially separated. The endosperm, occupying the vast majority of the seed volume, comprised both the chalazal and outer endosperm, and the endosperm developed and matured before the embryo. As the seed developed, the seed coat was gradually lignified and consisted of palisade-like stone cells externally and epidermal layer cells internally.

Screening and identification of photoresponse factors in kiwifruit (Actinidia arguta) development

BACKGROUND

Light is essential for kiwifruit development, in which photoresponse factors contributes greatly to the quality formation. 'Light sensitive hypocotyls, also known as light-dependent short hypocotyls' (LSH) gene family can participate in fruit development as photoresponse factor. However, the key LSH gene that determine kiwifruit development remains unclear. This study aim to screen and identify the key gene AaLSH9 in A. arguta.

MATERIALS AND METHODS

Genome-wide identification of the LSH gene family was used to analyse LSH genes in kiwifruit. Homologous cloning was used to confirm the sequence of candidate LSH genes. qRT-PCR and cluster analysis of expression pattern were used to screen the key AaLSH9 gene. Subcellular localization of AaLSH9 in tobacco leaves and overexpression of AaLSH9 in Arabidopsis thaliana hy5 mutant plants were used to define the acting place in cell and identify molecular function, respectively.

RESULTS

We identified 15 LSH genes, which were divided into two sub-families namely A and B. Domain analysis of A and B showed that they contained different domain organizations, which possibly played key roles in the evolution process. Three LSH genes, AaLSH2, AaLSH9, and AaLSH11, were successfully isolated from Actinidia arguta. The expression pattern and cluster analysis of these three AaLSH genes suggested AaLSH9 might be a key photoresponse gene participating in fruit development in A. arguta. Subcellular localization showed AaLSH9 protein was located in the nucleus. The overexpression of AaLSH9 gene in Arabidopsis thaliana hy5 mutant plants partially complemented the long hypocotyls of hy5 mutant, implying AaLSH9 played a key role as photoresponse factor in cells. In addition, the seed coat color of A. thaliana over-expressing AaLSH9 became lighter than the wide type A.thaliana. Finally, AaCOP1 was confirmed as photoresponse factor to participate in developmental process by stable transgenic A. thaliana.

CONCLUSIONS

AaLSH9 can be involved in kiwifruit (A. arguta) development as key photoresponse factor. Our results not only identified the photoresponse factors AaLSH9 and AaCOP1 but also provided insights into their key role in fruit quality improvement in the process of light response.

Genome-wide identification and expression analysis of the JMJ-C gene family in melon (Cucumis melo L.) reveals their potential role in fruit development

BACKGROUND

Proteins with the jumonji (JMJ)-C domain belong to the histone demethylase family and contribute to reverse histone methylation. Although JMJ-C family genes have an essential role in regulating plant growth and development, the characterization of the JMJ-C family genes in melon has not been uncovered.

RESULTS

In this study, a total of 17 JMJ-C proteins were identified in melon (Cucumis melo L.). CmJMJs were categorized into five subfamilies based on the specific conserved domain: KDM4/JHDM3, KDM5/JARID1, JMJD6, KDM3/JHDM2, and JMJ-C domain-only. The chromosome localization analyses showed that 17 CmJMJs were distributed on nine chromosomes. Cis-acting element analyses of the 17 CmJMJ genes showed numerous hormone, light, and stress response elements distributed in the promoter region. Covariance analysis revealed one pair of replicated fragments (CmJMJ3a and CmJMJ3b) in 17 CmJMJ genes. We investigated the expression profile of 17 CmJMJ genes in different lateral organs and four developmental stages of fruit by RNA-seq transcriptome analysis and RT-qPCR. The results revealed that most CmJMJ genes were prominently expressed in female flowers, ovaries, and developing fruits, suggesting their active role in melon fruit development. Subcellular localization showed that the fruit-related CmJMJ5a protein is specifically localized in the cell nucleus.

CONCLUSIONS

This study provides a comprehensive understanding of the gene structure, classification, and evolution of JMJ-C in melon and supports the clarification of the JMJ-C functions in further research.

Loss of Diel Circadian Clock Gene Cycling Is a Part of Grape Berry Ripening

Diel cycles of gene expression are thought to adapt plants to 24-h changes in environmental conditions. The circadian clock contributes to this process, but less is known about circadian programs in developing reproductive organs. While model plants and controlled conditions have contributed greatly to our knowledge of circadian clock function, there is a need to better understand its role in crop plants under field conditions with fluctuating light and temperature. In this study, we investigated changes in the circadian clock during the development of grape berries of Vitis vinifera L. We found that the transcripts of circadian clock homologs had high-amplitude oscillations prior to, but not during, ripening. As ripening progressed, the amplitude and rhythmicity of the diel oscillations decreased until most transcripts tested had no significant fluctuation over the 24-h cycle. Despite this loss of rhythmicity, the majority of circadian clock genes investigated were expressed at or near their abundance at the nadir of their pre-ripening oscillation although the berries remained transcriptionally active. From this, it can be concluded that cycling of the canonical circadian clock appears unnecessary for berry ripening. Our data suggest that changes in circadian clock dynamics during reproductive organ development may have important functional consequences.

Identification of watermelon H3K4 and H3K27 genes and their expression profiles during watermelon fruit development

BACKGROUND

The ClaH3K4s and ClaH3K27s gene families are subfamilies of the SET family, each with a highly conserved SET structure domain and a PHD structural domain. Both participate in histone protein methylation, which affects the chromosome structure and gene expression, and is essential for fruit growth and development.

METHODS AND RESULTS

In order to demonstrate the structure and expression characteristics of ClaH3K4s and ClaH3K27s in watermelon, members of the watermelon H3K4 and H3K27 gene families were identified, and their chromosomal localization, gene structure, and protein structural domains were analyzed. The phylogeny and covariance of the gene families with other species were subsequently determined, and the expression profiles were obtained by performing RNA-Seq and qRT-PCR. The watermelon genome had five H3K4 genes with 3207-8043 bp nucleotide sequence lengths and four H3K27 genes with a 1107-5499 bp nucleotide sequence. Synteny analysis revealed the close relationship between watermelon and cucumber, with the majority of members displaying a one-to-one covariance. Approximately half of the 'Hua-Jing 13 watermelon' ClaH3K4s and ClaH3K27s genes were expressed more in the late fruit development stages, while the changes were minimal for the remaining half. H3K4-2 expression was observed to be slightly greater on day 21 compared to other periods. Moreover, ClaH3K27-1 and ClaH3K27-2 were hardly expressed throughout the developing period, and ClaH3K27-4 exhibited the highest expression.

CONCLUSION

These results serve as a basis for further functional characterization of the H3K4 and H3K27 genes in the fruit development of watermelon.

Brassinosteroid catabolic enzyme CYP734A129 regulates the morphologies of leaves and floral organs in woodland strawberry

Brassinosteroids (BRs) play critical roles in plant growth and development and regulate many important agronomic traits. However, the functions of BRs in strawberry are unclear. This study identified two mutants, named P6 and R87, in woodland strawberry (Fragaria vesca) from EMS mutagenesis populations that exhibit narrow leaves, petals and sepals. Mapping by sequencing and genetic studies revealed that the F. vesca CYP734A129, encoding a putative BR catabolic enzyme, is the causative gene for both P6 and R87. Overexpression of CYP734A129 in both F. vesca and Arabidopsis causes a severe dwarf phenotype, and the BRI1-EMS-SUPPRESSOR 1 (BES1) protein is less abundant in the CYP734A129-overexpressing Arabidopsis seedlings. This suggests that CYP734A129 is functionally conserved with CYP734A1, as a BR-inactivating enzyme. Transcriptome analysis of young leaves revealed that four BR biosynthetic genes were significantly downregulated in P6 (cyp734a129), and photosynthesis-related genes were highly enriched among the up-regulated genes in P6 compared to the wild type. This further supports that CYP734A129 inactivates BRs in F. vesca. Furthermore, we showed that mutations in CYP734A129 do not affect fruit shape and color during ripening in strawberry. Overall, our results suggest that F. vesca CYP734A129 is a BR catabolic enzyme, and provide insights into the roles of CYP734A129 in strawberry.

NAC transcription factor NOR-like1 regulates tomato fruit size

MAIN CONCLUSION

NOR-like1 regulates tomato fruit size by targeting SlARF9, SlGRAS2, SlFW3.2, and SlFW11.3 genes involved in cell division and cell expansion. Fruit size is an important agricultural character that determines the yield of crops. Here, we found that NAC transcription factor NOR-like1 regulated fruit size by regulating cell layer number and cell area in tomato. Over-expressing NOR-like1 gene in tomato reduced fruit weight and size, whereas the knock-out of NOR-like1 increased fruit weight and size. At the molecular level, NOR-like1 binds to the promoter of SlGRAS2, SlFW3.2, and SlFW11.3 to repress their transcription, while it also binds to the promoter of ARF9 to activate its transcription. Overall, these results expand the biological function of NOR-like1 and deepen our understanding of the transcriptional network that regulates tomato fruit size.

Genome-wide characterization of the cellulose synthase gene family in Ziziphus jujuba reveals its function in cellulose biosynthesis during fruit development

The cellulose synthase (Ces/Csl) is a key enzyme in plant cellulose synthesis. Jujube fruits are rich in cellulose. 29 ZjCesA/Csl genes were identified in jujube genome and showed tissue-specific expression. 13 genes highly expressed in jujube fruit exhibited obviously sequential expressions during the fruit development, indicating that they might play distinct roles during the process. Meanwhile, the correlation analysis showed the expressions of ZjCesA1 and ZjCslA1 were significant positive related to the cellulose synthase activities. Furthermore, transient overexpressions of ZjCesA1 or ZjCslA1 in jujube fruits significantly increased cellulose synthase activities and contents, whereas silencing of ZjCesA1 or ZjCslA1 in jujube seedlings obviously reduced cellulose levels. Moreover, the Y2H assays verified that ZjCesA1 and ZjCslA1 may participate in cellulose synthesis by forming protein complexes. The study not only reveals the bioinformatics characteristics and functions of cellulose synthase genes in jujube, but also provides clues for studying cellulose synthesis in other fruits.

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.

Analysis of flavonoid metabolism during fruit development of Lycium chinense

Lycium chinense is an important medicinal plant in the northwest of China. Flavonoids are the major pharmacological components of L. chinense fruits. However, flavonoid metabolism during fruit development of L. chinense remains to be studied. Here, we analyzed the change of flavonoid contents, enzyme activity, and gene expression during fruit development of L. chinense. We found that flavonoids, anthocyanins, and catechins are the most important components of L. chinense fruits. Flavonoid content was increased with fruit development and was high at the late developmental stage. PAL, CHS, and F3H enzymes played a significant role in flavonoid accumulation in fruits. Transcriptomic analysis showed that anthocyanin pathway, flavonol pathway, flavonoid biosynthesis, and phenylpropanoid synthesis pathway were the major pathways involved in flavonoid metabolism in L. chinense. Gene expression analysis indicated that PAL1 and CHS2 genes were critical for flavonoid metabolism in L. chinense fruits. These discoveries help us understand the dynamic changes in flavonoids during fruit development and enhance the use of L. chinense fruits.

Active compound analysis of Ziziphus jujuba cv. Jinsixiaozao in different developmental stages using metabolomic and transcriptomic approaches

Jujube (Ziziphus jujuba Mill.) is a popular fruit with health benefits ascribed to its various metabolites. These metabolites determine the flavors and bioactivities of the fruit, as well as their desirability. However, the dynamics of the metabolite composition and the underlying gene expression that modulate the overall flavor and accumulation of active ingredients during fruit development remain largely unknown. Therefore, we conducted an integrated metabolomic and transcriptomic investigation covering various developmental stages in the jujube cultivar Z. jujuba cv. Jinsixiaozao, which is famous for its nutritional and bioactive properties. A total of 407 metabolites were detected by non-targeted metabolomics. Metabolite accumulation during different jujube developmental stages was examined. Most nucleotides and amino acids and their derivatives accumulated during development, with cAMP increasing notably during ripening. Triterpenes gradually accumulated and were maintained at high concentrations during ripening. Many flavonoids were maintained at relatively high levels in early development, but then rapidly decreased later. Transcriptomic and metabolomic analyses revealed that chalcone synthase (CHS), chalcone isomerase (CHI), flavonol synthase (FLS), and dihydroflavonol 4-reductase (DFR) were mainly responsible for regulating the accumulation of flavonoids. Therefore, the extensive downregulation of these genes was probably responsible for the decreases in flavonoid content during fruit ripening. This study provide an overview of changes of active components in 'Jinsixiaozao' during development and ripening. These findings enhance our understanding of flavor formation and will facilitate jujube breeding for improving both nutrition and function.

Identification and expression analysis of BURP domain-containing genes in jujube and their involvement in low temperature and drought response

Background

Plant-specific BURP domain-containing genes are involved in plant development and stress responses. However, the role of BURP family in jujube (Ziziphus jujuba Mill.) has not been investigated.

Results

In this study, 17 BURP genes belonging to four subfamilies were identified in jujube based on homology analysis, gene structures, and conserved motif confirmation. Gene duplication analysis indicated both tandem duplication and segmental duplication had contributed to ZjBURP expansion. The ZjBURPs were extensively expressed in flowers, young fruits, and jujube leaves. Transcriptomic data and qRT-PCR analysis further revealed that ZjBURPs also significantly influence fruit development, and most genes could be induced by low temperature, salinity, and drought stresses. Notably, several BURP genes significantly altered expression in response to low temperature (ZjPG1) and drought stresses (ZjBNM7, ZjBNM8, and ZjBNM9).

Conclusions

These results provided insights into the possible roles of ZjBURPs in jujube development and stress response. These findings would help selecting candidate ZjBURP genes for cold- and drought-tolerant jujube breeding.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-022-08907-9.

Identification and characterization of auxin response factor (ARF) family members involved in fig (Ficus carica L.) fruit development

The auxin response factor (ARF) combines with AuxREs cis-acting elements in response to auxin to regulate plant development. To date, no comprehensive analysis of ARF genes expressed during fruit development has been conducted for common fig (Ficus carica L.). In this study, members of the FcARF gene family were screened, identified in the fig genome database and their features characterized using bioinformatics. Twenty FcARF genes were clustered into three classes, with almost similar highly conserved DBD (B3-like DNA binding domain), AUX/IAA (auxin/indole-3-acetic acid gene family) and MR domain structure among class members. Analysis of amino acid species in MR domain revealed 10 potential transcription activators and 10 transcription inhibitors, and 17 FcARF members were predicted to be located in the nucleus. DNA sequence analysis showed that the ARF gene family consisted of 4-25 exons, and the promoter region contained 16 cis-acting elements involved in stress response, hormone response and flavonoid biosynthesis. ARF genes were expressed in most tissues of fig, especially flower and peel. Transcriptomics analysis results showed that FcARF2, FcARF11 and FcARF12, belonging to class-Ia, were stably and highly expressed in the early development stage of flower and peel of 'Purple peel' fig. However, their expression levels decreased after maturity. Expression of class-Ic member FcARF3 conformed to the regularity of fig fruit development. These four potential transcription inhibitors may regulate fruit growth and development of 'Purple Peel' fig. This study provides comprehensive information on the fig ARF gene family, including gene structure, chromosome position, phylogenetic relationship and expression pattern. Our work provides a foundation for further research on auxin-mediated fig fruit development.

Comparative transcriptomics identifies candidate genes involved in the evolutionary transition from dehiscent to indehiscent fruits in Lepidium (Brassicaceae)

BACKGROUND

Fruits are the seed-bearing structures of flowering plants and are highly diverse in terms of morphology, texture and maturation. Dehiscent fruits split open upon maturation to discharge their seeds while indehiscent fruits are dispersed as a whole. Indehiscent fruits evolved from dehiscent fruits several times independently in the crucifer family (Brassicaceae). The fruits of Lepidium appelianum, for example, are indehiscent while the fruits of the closely related L. campestre are dehiscent. Here, we investigate the molecular and genetic mechanisms underlying the evolutionary transition from dehiscent to indehiscent fruits using these two Lepidium species as model system.

RESULTS

We have sequenced the transcriptomes and small RNAs of floral buds, flowers and fruits of L. appelianum and L. campestre and analyzed differentially expressed genes (DEGs) and differently differentially expressed genes (DDEGs). DEGs are genes that show significantly different transcript levels in the same structures (buds, flowers and fruits) in different species, or in different structures in the same species. DDEGs are genes for which the change in expression level between two structures is significantly different in one species than in the other. Comparing the two species, the highest number of DEGs was found in flowers, followed by fruits and floral buds while the highest number of DDEGs was found in fruits versus flowers followed by flowers versus floral buds. Several gene ontology terms related to cell wall synthesis and degradation were overrepresented in different sets of DEGs highlighting the importance of these processes for fruit opening. Furthermore, the fruit valve identity genes FRUITFULL and YABBY3 were among the DEGs identified. Finally, the microRNA miR166 as well as the TCP transcription factors BRANCHED1 (BRC1) and TCP FAMILY TRANSCRIPTION FACTOR 4 (TCP4) were found to be DDEGs.

CONCLUSIONS

Our study reveals differences in gene expression between dehiscent and indehiscent fruits and uncovers miR166, BRC1 and TCP4 as candidate genes for the evolutionary transition from dehiscent to indehiscent fruits in Lepidium.

Proteome Analysis of Vacuoles Isolated from Fig (Ficus carica L.) Flesh during Fruit Development

Fruit flesh cell vacuoles play a pivotal role in fruit growth and quality formation. In the present study, intact vacuoles were carefully released and collected from protoplasts isolated from flesh cells at five sampling times along fig fruit development. Label-free quantification and vacuole proteomic analysis identified 1,251 proteins, 1,137 of which were recruited as differentially abundant proteins (DAPs) by fold change ≥ 1.5, P < 0.05. DAPs were assigned to 10 functional categories; among them, 238, 186, 109, 93 and 90 were annotated as metabolism, transport proteins, membrane fusion or vesicle trafficking, protein fate and stress response proteins, respectively. Decreased numbers of DAPs were uncovered along fruit development. The overall changing pattern of DAPs revealed two major proteome landscape conversions in fig flesh cell vacuoles: the first occurred when fruit developed from late-stage I to mid-stage II, and the second occurred when the fruit started ripening. Metabolic proteins related to glycosidase, lipid and extracellular proteins contributing to carbohydrate storage and vacuole expansion, and protein-degrading proteins determining vacuolar lytic function were revealed. Key tonoplast proteins contributing to vacuole expansion, cell growth and fruit quality formation were also identified. The revealed comprehensive changes in the vacuole proteome during flesh development were compared with our previously published vacuole proteome of grape berry. The information expands our knowledge of the vacuolar proteome and the protein basis of vacuole functional evolution during fruit development and quality formation.

Genome-wide characterization and expression analysis of SAUR gene family in Melon (Cucumis melo L.)

We identified 66 melon SAUR genes by bioinformatic analyses. CmSAUR19, 38, 58, 62 genes are specifically expressed in different stages of fruit growth, suggesting their participation in regulating fruit development. Auxin plays a crucial role in plant growth by regulating the multiple auxin response genes. However, in melon (Cucumis melo L.), the functions of the auxin early response gene family SAUR (Small auxin up RNA) genes in fruit development are still poorly understood. Through genome-wide characterization of CmSAUR family in melon, we identified a total of 66 CmSAUR genes. The open reading frames of the CmSAUR genes ranged from 234 to 525 bp, containing only one exon and lacking introns. Chromosomal position and phylogenetic tree analyses found that the two gene clusters in the melon chromosome are highly homologous in the Cucurbitaceae plants. Among the four conserved motifs in CmSAUR proteins, motif 1, motif 2, and motif 3 located in most of the family protein sequences, and motif 4 showed a close correlation with the two gene clusters. The CmSAUR28 and CmSAUR58 genes have auxin response elements located in the promoters, suggesting they may be involved in the auxin signaling pathway to regulate fruit development. Through transcriptomic profiling in the four developmental stages of fruit and different lateral organs, we selected 16 differentially-expressed SAUR genes for performing further expression analyses. qRT-PCR results showed that five SAUR genes are specifically expressed in flower organs and ovaries. CmSAUR19 and CmSAUR58 were significantly accumulated in the early developmental stage of the fruit. CmSAUR38 and CmAUR62 showed high expression in the climacteric and post-climacteric stages, suggesting their specific role in controlling fruit ripening. This work provides a foundation for further exploring the function of the SAUR gene in fruit development.

Identification of ARF family in blueberry and its potential involvement of fruit development and pH stress response

Background

Auxin responsive factor (ARF) family is one of core components in auxin signalling pathway, which governs diverse developmental processes and stress responses. Blueberry is an economically important berry-bearing crop and prefers to acidic soil. However, the understandings of ARF family has not yet been reported in blueberry.

Results

In the present study, 60 ARF genes (VcARF) were identified in blueberry, and they showed diverse gene structures and motif compositions among the groups and similar within each group in the phylogenetic tree. Noticeably, 9 digenic, 5 trigenic and 6 tetragenic VcARF pairs exhibited more than 95% identity to each other. Computational analysis indicated that 23 VcARFs harbored the miRNA responsive element (MRE) of miR160 or miR167 like other plant ARF genes. Interestingly, the MRE of miR156d/h-3p was observed in the 5’UTR of 3 VcARFs, suggesting a potentially novel post-transcriptional control. Furthermore, the transcript accumulations of VcARFs were investigated during fruit development, and three categories of transcript profiles were observed, implying different functional roles. Meanwhile, the expressions of VcARFs to different pH conditions (pH4.5 and pH6.5) were surveyed in pH-sensitive and tolerant blueberry species, and a number of VcARFs showed different transcript accumulations. More importantly, distinct transcriptional response to pH stress (pH6.5) were observed for several VcARFs (such as VcARF6s and VcARF19-3/19–4) between pH-sensitive and tolerant species, suggesting their potential roles in adaption to pH stress.

Conclusions

Sixty VcARF genes were identified and characterized, and their transcript profiles were surveyed during fruit development and in response to pH stress. These findings will contribute to future research for eliciting the functional roles of VcARFs and regulatory mechanisms, especially fruit development and adaption to pH stress.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-022-08556-y.

"Phoenix in Flight": an unique fruit morphology ensures wind dispersal of seeds of the phoenix tree (Firmiana simplex (L.) W. Wight)

BACKGROUND

Many seed plants produce winged diaspores that use wind to disperse their seeds. The morphology of these diaspores is directly related to the seed dispersal potential. The majority of winged diaspores have flat wings and only seeds; however, some angiosperms, such as Firmiana produce winged fruit with a different morphology, whose seed dispersal mechanisms are not yet fully understood. In this study, we observed the fruit development of F. simplex and determined the morphological characteristics of mature fruit and their effects on the flight performance of the fruit.

RESULTS

We found that the pericarp of F. simplex dehisced early and continued to unfold and expand during fruit development until ripening, finally formed a spoon-shaped wing with multiple alternate seeds on each edge. The wing caused mature fruit to spin stably during descent to provide a low terminal velocity, which was correlated with the wing loading and the distribution of seeds on the pericarp. When the curvature distribution of the pericarp surface substantially changed, the aerodynamic characteristics of fruit during descent altered, resulting in the inability of the fruit to spin.

CONCLUSIONS

Our results suggest that the curved shape and alternate seed distribution are necessary for the winged diaspore of F. simplex to stabilize spinning during wind dispersal. These unique morphological characteristics are related to the early cracking of fruits during development, which may be an adaptation for the wind dispersal of seeds.

Genome-wide identification and expression profile of YABBY genes in Averrhoa carambola

Background

Members of the plant-specific YABBY gene family are thought to play an important role in the development of leaf, flower, and fruit. The YABBY genes have been characterized and regarded as vital contributors to fruit development in Arabidopsis thaliana and tomato, in contrast to that in the important tropical economic fruit star fruit (Averrhoa carambola), even though its genome is available.

Methods

In the present study, a total of eight YABBY family genes (named from AcYABBY1 to AcYABBY8) were identified from the genome of star fruit, and their phylogenetic relationships, functional domains and motif compositions, physicochemical properties, chromosome locations, gene structures, protomer elements, collinear analysis, selective pressure, and expression profiles were further analyzed.

Results

Eight AcYABBY genes (AcYABBYs) were clustered into five clades and were distributed on five chromosomes, and all of them had undergone negative selection. Tandem and fragment duplications rather than WGD contributed to YABBY gene number in the star fruit. Expression profiles of AcYABBYs from different organs and developmental stages of fleshy fruit indicated that AcYABBY4 may play a specific role in regulating fruit size. These results emphasize the need for further studies on the functions of AcYABBYs in fruit development.

Long-read transcriptome sequencing provides insight into lignan biosynthesis during fruit development in Schisandra chinensis

BACKGROUND

Schisandra chinensis, an ancient member of the most basal angiosperm lineage which is known as the ANITA, is a fruit-bearing vine with the pharmacological effects of a multidrug system, such as antioxidant, anti-inflammatory, cardioprotective, neuroprotective, anti-osteoporosis effects. Its major bioactive compound is represented by lignans such as schisandrin. Molecular characterization of lignan biosynthesis in S. chinensis is of great importance for improving the production of this class of active compound. However, the biosynthetic mechanism of schisandrin remains largely unknown.

RESULTS

To understand the potential key catalytic steps and their regulation of schisandrin biosynthesis, we generated genome-wide transcriptome data from three different tissues of S. chinensis cultivar Cheongsoon, including leaf, root, and fruit, via long- and short-read sequencing technologies. A total of 132,856 assembled transcripts were generated with an average length of 1.9 kb and high assembly completeness. Overall, our data presented effective, accurate gene annotation in the prediction of functional pathways. In particular, the annotation revealed the abundance of transcripts related to phenylpropanoid biosynthesis. Remarkably, transcriptome profiling during fruit development of S. chinensis cultivar Cheongsoon revealed that the phenylpropanoid biosynthetic pathway, specific to coniferyl alcohol biosynthesis, showed a tendency to be upregulated at the postfruit development stage. Further the analysis also revealed that the pathway forms a transcriptional network with fruit ripening-related genes, especially the ABA signaling-related pathway. Finally, candidate unigenes homologous to isoeugenol synthase 1 (IGS1) and dirigent-like protein (DIR), which are subsequently activated by phenylpropanoid biosynthesis and thus catalyze key upstream steps in schisandrin biosynthesis, were identified. Their expression was increased at the postfruit development stage, suggesting that they may be involved in the regulation of schisandrin biosynthesis in S. chinensis.

CONCLUSIONS

Our results provide new insights into the production and accumulation of schisandrin in S. chinensis berries and will be utilized as a valuable transcriptomic resource for improving the schisandrin content.

Tomato SlPP2C5 Is Involved in the Regulation of Fruit Development and Ripening

Abscisic acid (ABA) regulates plant development mainly through its signaling, in which ABA binds to receptors to inhibit type 2C protein phosphatases (PP2Cs). The exact roles of PP2Cs in fruit development are still unclear. In this work, we verify that tomato SlPP2C5 works as a negative regulator in ABA signaling during fruit development. SlPP2C5 was inhibited by both monomeric and dimeric ABA receptors SlPYLs through ABA dose-dependent way, and it interacted physically with SlPYLs and SlSnRK2s. SlPP2C5 was highly expressed in fruits induced by exogenous ABA. Plants with overexpressed SlPP2C5 had lower sensitivity to ABA, which showed faster seed germination and primary root growth compared to Wild type (WT), while SlPP2C5-suppressed plants were more sensitive to ABA. SlPP2C5-over-expression (OE) delayed fruit ripening onset, while SlPP2C5-RNAi advanced fruit ripening. Alteration of SlPP2C5 expression impacts fruit quality parameters as well, including pericarp thickness, fruit shape index, seed number and weight and the soluble solid content. RNA-seq analysis revealed that there were significant expression differences of genes related to ethylene release and lycopene synthesis between WT and both SlPP2C5-OE and SlPP2C5-RNAi lines with an inversed variation. Taken together, our findings demonstrate that SlPP2C5 plays an important role in the regulation of fruit development, ripening and quality.

Plasmodesmata and their role in the regulation of phloem unloading during fruit development

Fruit consumption is fundamental to a balanced diet. The contemporary challenge of maintaining a steady food supply to meet the demands of a growing population is driving the development of strategies to improve the production and nutritional quality of fruit. Plasmodesmata, the structures that mediate symplasmic transport between plant cells, play an important role in phloem unloading and distribution of sugars and signalling molecules into developing organs. Targeted modifications to the structures and functioning of plasmodesmata have the potential to improve fruit development; however, knowledge on the mechanisms underpinning plasmodesmata regulation in this context is scarce. In this review, we have compiled current knowledge on plasmodesmata and their structural characterisation during the development of fruit organs. We discuss key questions on phloem unloading, including the pathway shift from symplasmic to apoplastic that takes place during the onset of ripening as potential targets for improving fruit quality.

Genome-wide investigation of the GRAS transcription factor family in foxtail millet (Setaria italica L.)

BACKGROUND

GRAS transcription factors perform indispensable functions in various biological processes, such as plant growth, fruit development, and biotic and abiotic stress responses. The development of whole-genome sequencing has allowed the GRAS gene family to be identified and characterized in many species. However, thorough in-depth identification or systematic analysis of GRAS family genes in foxtail millet has not been conducted.

RESULTS

In this study, 57 GRAS genes of foxtail millet (SiGRASs) were identified and renamed according to the chromosomal distribution of the SiGRAS genes. Based on the number of conserved domains and gene structure, the SiGRAS genes were divided into 13 subfamilies via phylogenetic tree analysis. The GRAS genes were unevenly distributed on nine chromosomes, and members of the same subfamily had similar gene structures and motif compositions. Genetic structure analysis showed that most SiGRAS genes lacked introns. Some SiGRAS genes were derived from gene duplication events, and segmental duplications may have contributed more to GRAS gene family expansion than tandem duplications. Quantitative polymerase chain reaction showed significant differences in the expression of SiGRAS genes in different tissues and stages of fruits development, which indicated the complexity of the physiological functions of SiGRAS. In addition, exogenous paclobutrazol treatment significantly altered the transcription levels of DELLA subfamily members, downregulated the gibberellin content, and decreased the plant height of foxtail millet, while it increased the fruit weight. In addition, SiGRAS13 and SiGRAS25 may have the potential for genetic improvement and functional gene research in foxtail millet.

CONCLUSIONS

Collectively, this study will be helpful for further analysing the biological function of SiGRAS. Our results may contribute to improving the genetic breeding of foxtail millet.

Genome-wide analysis of the NF-Y gene family and their roles in relation to fruit development in Tartary buckwheat (Fagopyrum tataricum)

Nuclear factor Y (NF-Y) is a heterotrimeric transcription factor playing crucial roles in various biological process in plant. However, thorough research on NF-Y gene family of Tartary buckwheat (Fagopyrum tataricum) is little. In this study, 38 FtNF-Y genes (12 FtNF-YAs, 17 FtNF-YBs, and 9 FtNF-YCs) were identified and renamed on the basis of their subfamily and chromosomal location. Their gene structure, genomic mapping, motif composition, conserved domain, phylogenetic relationships, cis-acting elements and gene expression were investigated. Illustration of gene structures and conserved domains of FtNF-Ys revealed their functional conservation and specificity. Construction of phylogenetic trees of NF-Ys in Tartary buckwheat, Arabidopsis, tomato, rice and banana, allowed us to predict functional similarities among NF-Ys from different species. Gene expression analysis displayed that twenty-four FtNF-Ys were expressed in all the tissues and the transcript levels of them were different, suggesting their function varieties. Moreover, expression profiles of twenty FtNF-Ys along five different fruit development stages acquired by real-time quantitative PCR (RT-qPCR) demonstrated distinct abundance diversity at different stages, providing some clues of potential fruit development regulators. Our study could provide helpful reference information for further function characterization of FtNF-Ys and for the fruit quality enhancement of Tartary buckwheat.

Alternative splicing regulation appears to play a crucial role in grape berry development and is also potentially involved in adaptation responses to the environment

BACKGROUND

Alternative splicing (AS) produces transcript variants playing potential roles in proteome diversification and gene expression regulation. AS modulation is thus essential to respond to developmental and environmental stimuli. In grapevine, a better understanding of berry development is crucial for implementing breeding and viticultural strategies allowing adaptation to climate changes. Although profound changes in gene transcription have been shown to occur in the course of berry ripening, no detailed study on splicing modifications during this period has been published so far. We report here on the regulation of gene AS in developing berries of two grapevine (Vitis vinifera L.) varieties, Gewurztraminer (Gw) and Riesling (Ri), showing distinctive phenotypic characteristics. Using the software rMATS, the transcriptomes of berries at four developmental steps, from the green stage to mid-ripening, were analysed in pairwise comparisons between stages and varieties.

RESULTS

A total of 305 differential AS (DAS) events, affecting 258 genes, were identified. Interestingly, 22% of these AS events had not been reported before. Among the 80 genes that underwent the most significant variations during ripening, 22 showed a similar splicing profile in Gw and Ri, which suggests their involvement in berry development. Conversely, 23 genes were subjected to splicing regulation in only one variety. In addition, the ratios of alternative isoforms were different in Gw and Ri for 35 other genes, without any change during ripening. This last result indicates substantial AS differences between the two varieties. Remarkably, 8 AS events were specific to one variety, due to the lack of a splice site in the other variety. Furthermore, the transcription rates of the genes affected by stage-dependent splicing regulation were mostly unchanged, identifying AS modulation as an independent way of shaping the transcriptome.

CONCLUSIONS

The analysis of AS profiles in grapevine varieties with contrasting phenotypes revealed some similarity in the regulation of several genes with developmental functions, suggesting their involvement in berry ripening. Additionally, many splicing differences were discovered between the two varieties, that could be linked to phenotypic specificities and distinct adaptive capacities. Together, these findings open perspectives for a better understanding of berry development and for the selection of grapevine genotypes adapted to climate change.

Multiple and integrated functions of floral C-class MADS-box genes in flower and fruit development of Physalis floridana

This work reveals potentially multiple and integrated roles in flower and fruit development of floral C-class MADS-box genes in Physalis. The Physalis fruit features a morphological novelty, the Chinese lantern. Floral C-class MADS-domain AGAMOUS-like (AG-like) proteins can interact with the identified regulators of this novel structure. However, the developmental role of the floral C-class genes is unknown in Physalis. Here, we characterized two AG-like genes from Physalis floridana, designated PFAG1 and PFAG2. The two paralogous genes shared around 61.0% of sequence identity and had similar expression domains, with different expression levels in the floral and berry development. However, the genes had distinct expression patterns in leaf and calyx development. Protein-protein interaction analyses revealed that PFAG1 and PFAG2 could commonly or specifically dimerize with certain floral MADS-domain proteins as well as non-MADS-domain proteins involved in various floral developmental processes. Gene downregulation analyses demonstrated that PFAG1 may repress PFAG2, but PFAG2 did not affect PFAG1. Downregulating PFAG1 led to incomplete floral homeotic variation in the stamens and carpels, and alteration of petal coloration pattern, while downregulating PFAG2 did not result in any floral homeotic variation. PFAG1 affected pollen maturation, while PFAG2 affected female fertility. However, simultaneously downregulating PFAG1 and PFAG2 caused loss of the complete C-function, indicating that the two PFAG genes interact to determine the identity and functionality of androecia and gynoecia organs. Their potential roles in regulating fruit size and the Chinese lantern are also discussed. Our results reveal functional divergence of floral C-class MADS-box genes in Physalis, demonstrating that they may play multiple and integrated roles in flower and fruit development.

Comparative study on fruit development and oil synthesis in two cultivars of Camellia oleifera

BACKGROUND

The oil-tea tree (Camellia oleifera Abel.) is a woody tree species that produces edible oil in the seed. C. oleifera oil has high nutritional value and is also an important raw material for medicine and cosmetics. In China, due to the uncertainty on maturity period and oil synthesis mechanism of many C. oleifera cultivars, growers may harvest fruits prematurely, which could not maximize fruit and oil yields. In this study, our objective was to explore the mechanism and differences of oil synthesis between two Camellia oleifera cultivars for a precise definition of the fruit ripening period and the selection of appropriate cultivars.

RESULTS

The results showed that 'Huashuo' had smaller fruits and seeds, lower dry seed weight and lower expression levels of fatty acid biosynthesis genes in July. We could not detect the presence of oil and oil bodies in 'Huashuo' seeds until August, and oil and oil bodies were detected in 'Huajin' seeds in July. Moreover, 'Huashuo' seeds were not completely blackened in October with up to 60.38% of water and approximately 37.98% of oil in seed kernels whose oil content was much lower than normal mature seed kernels. The oil bodies in seed endosperm cells of 'Huajin' were always higher than those of 'Huashuo' from July to October.

CONCLUSION

Our results confirmed that C. oleifera 'Huashuo' fruits matured at a lower rate compared to 'Huajin' fruits and that 'Huajin' seeds entered the oil synthesis period earlier than 'Huashuo' seeds. Moreover, 'Huashuo' fruits did not mature during the Frost's Descent period (October 23-24 each year).

Transcriptomic analysis of temporal shifts in berry development between two grapevine cultivars of the Pinot family reveals potential genes controlling ripening time

BACKGROUND

Grapevine cultivars of the Pinot family represent clonally propagated mutants with major phenotypic and physiological differences, such as different colour or shifted ripening time, as well as changes in important viticultural traits. Specifically, the cultivars 'Pinot Noir' (PN) and 'Pinot Noir Precoce' (PNP, early ripening) flower at the same time, but vary in the beginning of berry ripening (veraison) and, consequently, harvest time. In addition to genotype, seasonal climatic conditions (i.e. high temperatures) also affect ripening times. To reveal possible regulatory genes that affect the timing of veraison onset, we investigated differences in gene expression profiles between PN and PNP throughout berry development with a closely meshed time series and over two separate years.

RESULTS

The difference in the duration of berry formation between PN and PNP was quantified to be approximately two weeks under the growth conditions applied, using plant material with a proven PN and PNP clonal relationship. Clusters of co-expressed genes and differentially expressed genes (DEGs) were detected which reflect the shift in the timing of veraison onset. Functional annotation of these DEGs fit to observed phenotypic and physiological changes during berry development. In total, we observed 3,342 DEGs in 2014 and 2,745 DEGs in 2017 between PN and PNP, with 1,923 DEGs across both years. Among these, 388 DEGs were identified as veraison-specific and 12 were considered as berry ripening time regulatory candidates. The expression profiles revealed two candidate genes for ripening time control which we designated VviRTIC1 and VviRTIC2 (VIT_210s0071g01145 and VIT_200s0366g00020, respectively). These genes likely contribute the phenotypic differences observed between PN and PNP.

CONCLUSIONS

Many of the 1,923 DEGs show highly similar expression profiles in both cultivars if the patterns are aligned according to developmental stage. In our work, putative genes differentially expressed between PNP and PN which could control ripening time as well as veraison-specific genes were identified. We point out connections of these genes to molecular events during berry development and discuss potential candidate genes which may control ripening time. Two of these candidates were observed to be differentially expressed in the early berry development phase. Several down-regulated genes during berry ripening are annotated as auxin response factors / ARFs. Conceivably, general changes in auxin signaling may cause the earlier ripening phenotype of PNP.

Comparative anatomical and transcriptomic insights into Vaccinium corymbosum flower bud and fruit throughout development

BACKGROUND

Blueberry (Vaccinium spp.) is characterized by the production of berries that are smaller than most common fruits, and the underlying mechanisms of fruit size in blueberry remain elusive. V. corymbosum 'O'Neal' and 'Bluerain' are commercial southern highbush blueberry cultivars with large- and small-size fruits, respectively, which mature 'O'Neal' fruits are 1 ~ 2-fold heavier than those of 'Bluerain'. In this study, the ontogenetical patterns of 'O'Neal' and 'Bluerain' hypanthia and fruits were compared, and comparative transcriptomic analysis was performed during early fruit development.

RESULTS

V. corymbosum 'O'Neal' and 'Bluerain' hypanthia and fruits exhibited intricate temporal and spatial cell proliferation and expansion patterns. Cell division before anthesis and cell expansion after fertilization were the major restricting factors, and outer mesocarp was the key tissue affecting fruit size variation among blueberry genotypes. Comparative transcriptomic and annotation analysis of differentially expressed genes revealed that the plant hormone signal transduction pathway was enriched, and that jasmonate-related TIFYs genes might be the key components orchestrating other phytohormones and influencing fruit size during early blueberry fruit development.

CONCLUSIONS

These results provided detailed ontogenetic evidence for determining blueberry fruit size, and revealed the important roles of phytohormone signal transductions involving in early fruit development. The TIFY genes could be useful as markers for large-size fruit selection in the current breeding programs of blueberry.

Downstream of GA4, PbCYP78A6 participates in regulating cell cycle-related genes and parthenogenesis in pear (Pyrus bretshneideri Rehd.)

BACKGROUND

Parthenocarpy results in traits attractive to both consumers and breeders, and it overcomes the obstacle of self-incompatibility in the fruit set of horticultural crops, including pear (Pyrus bretshneider). However, there is limited knowledge regarding the genetic and molecular mechanisms that regulate parthenogenesis.

RESULTS

Here, in a transcriptional comparison between pollination-dependent fruit and GA₄-induced parthenocarpy, PbCYP78A6 was identified and proposed as a candidate gene involved in parthenocarpy. PbCYP78A6 is similar to Arabidopsis thaliana CYP78A6 and highly expressed in pear hypanthia. The increased PbCYP78A6 expression, as assessed by RT-qPCR, was induced by pollination and GA₄ exposure. The ectopic overexpression of PbCYP78A6 contributed to parthenocarpic fruit production in tomato. The PbCYP78A6 expression coincided with fertilized and parthenocarpic fruitlets development and the expression of fruit development-related genes as assessed by cytological observations and RT-qPCR, respectively. PbCYP78A6 RNA interference and overexpression in pear calli revealed that the gene is an upstream regulator of specific fruit development-related genes in pear.

CONCLUSIONS

Our findings indicate that PbCYP78A6 plays a critical role in fruit formation and provide insights into controlling parthenocarpy.

Identification and Expression Analyses of the Special 14-3-3 Gene Family in Papaya and its Involvement in Fruit Development, Ripening, and Abiotic Stress Responses

Plant 14-3-3 proteins play key roles in regulating growth, development, and stress responses. However, little is known about this gene family in papaya (Carica papaya L.). We characterized eight 14-3-3 genes from the papaya genome and designed them as CpGRF1-8. Based on phylogenetic, conserved motif, and gene structure analyses, papaya CpGRFs were divided into ε and non-ε groups. Expression analysis showed differential and class-specific transcription patterns in different organs. Quantitative real-time polymerase chain reaction analysis showed that most CpGRFs had large changes in expression during fruit development and ripening. This indicated that the CpGRFs were involved in regulating fruit development and ripening. Significant expression changes occurred after cold, salt, and drought treatments in papaya seedlings, indicating that CpGRFs were also involved in signaling responses to abiotic stress. These results provide a transcription profile of 14-3-3 genes in organs, during fruit development and ripening and in response to stress. Some highly expressed, fruit-specific, and stress-responsive candidate CpGRFs will be identified for further genetic improvement of papayas.

CcBLH6, a bell-like homeodomain-containing transcription factor, regulates the fruit lignification pattern

CcBLH6 is a bell-like homeodomain-containing transcription factor that plays an important role of lignin biosynthesis in the control of fruit lignification pattern in Camellia chekiangoleosa. The fruit of Camellia chekiangoleosa has a unique lignification pattern that features with a thick pericarp containing a low level of lignification. Yet the fruit lignification pattern and the regulatory network of responsible gene transcription are poorly understood. Here, we characterized a bell-like homeodomain-containing (BLH) transcription factor from C. chekiangoleosa, CcBLH6, in the control of fruit lignification. CcBLH6 expression was highly correlated with the unique lignification pattern during fruit development. The ectopic expression of CcBLH6 promoted the lignification process of stem and root in Arabidopsis. We found that expression of genes related to lignin biosynthesis and its transcriptional regulation was altered in transgenic lines. In a Camellia callus-transformation system, overexpression of CcBLH6 greatly enhanced the expression of genes related to lignin biosynthesis and its transcriptional regulation was altered in transgenic lines. In the callus-transformation system, overexpression of CcBLH6 greatly enhanced the lignification of parenchymal cells, and the regulation of several genes involved in lignin accumulation was largely consistent between Arabidopsis and Camellia. Our study reveals a positive role of CcBLH6 in the regulation of lignin biosynthesis during fruit lignification in Camellia.

Developmental and water deficit-induced changes in hydraulic properties and xylem anatomy of tomato fruit and pedicels

Xylem water transport from the parent plant plays a crucial role in fruit growth, development, and the determination of quality. Attempts have been made to partition the hydraulic resistance of the pathway over the course of development, but no consensus has been reached. Furthermore, the issue has not been addressed in the context of changing plant and fruit water status under water deficit conditions. In this study, we have conducted a rigorous investigation into the developmental changes that occur in the hydraulic properties of tomato fruits and their pedicels under well-irrigated and water deficit conditions, based on hydraulic measurements, fruit rehydration, dye-tracing, light and electron microscopy, and flow modeling. We found that a decline in water transport capacity during development did not occur in the xylem pathway leading up to the fruit, but within the fruit itself, where the effect might reside either inside or outside of the xylem pathway. The developmental pattern of the hydraulic resistance of the xylem pathway was not significantly influenced by water deficit. The changes in xylem flow between the fruit and the parent plant resulting from the reduced driving force under water deficit could explain the reduced accumulation of water in the fruit. This study provides new insights that aid our understanding of xylem water transport in fleshy fruits and its sensitivity to water deficit from a hydraulic perspective.

Genome-wide identification and comprehensive analysis of NAC family genes involved in fruit development in kiwifruit (Actinidia)

BACKGROUND

NAC transcription factors (TFs) are plant-specific proteins encoded by a large gene family. They play important roles in diverse biological processes, such as plant growth and development, leaf senescence, and responses to biotic or abiotic stresses. Functions of a number of NAC TFs have been identified mainly in model plants. However, very few studies on NAC TFs have been conducted in the fruit tree of kiwifruit.

RESULTS

Genome-wide NAC genes were identified and their phylogeny, genomic structure, chromosomal location, synteny relationships, protein properties and conserved motifs were analyzed. In addition, the fruit developmental process was evaluated in a new kiwifruit cultivar of Actinidia eriantha 'Ganlu 1'. And expressions for all those NAC genes were analyzed by quantitative real-time PCR method in fruits of 'Ganlu 1' during its developmental process. Our research identified 142 NAC TFs which could be phylogenetically divided into 23 protein subfamilies. The genomic structures of those NAC genes indicated that their exons were between one and ten. Analysis of chromosomal locations suggested that 116 out of 142 NACs distributed on all the 29 kiwifruit chromosomes. In addition, genome-wide gene expression analysis showed that expressions of 125 out of 142 NAC genes could be detected in fruit samples.

CONCLUSION

Our comprehensive study provides novel information on NAC genes and expression patterns in kiwifruit fruit. This research would be helpful for future functional identification of NAC genes involved in kiwifruit fruit development.

Transcriptome profiling provides insights into the fruit color development of wild Lycium ruthenicum Murr. from Qinghai-Tibet Plateau

Lycium ruthenicum Murr. is an important ecological and economic species in the Qaidam Basin of Qinghai-Tibet Plateau. Its black fruits (BF) are rich in anthocyanins, which have health-promoting properties for humans and thus provide nutritional benefits for this plant. Although the fruit quality of natural white fruit (WF) is affected by the disappearance of pigmentation in phenotypes, this phenomenon provides an opportunity to unravel the complex color metabolic networks. In this study, anthocyanin profiling confirmed that WF was formed due to anthocyanin loss. Transcriptome analysis of BF and WF revealed 101,466 unigenes, 261 of which were identified as the putative homologs of color-related genes in other species. Genes encoding the enzymes involved in flavonoid biosynthesis were also identified systematically. The structural gene expression levels of chalcone synthase (CHS), chalcone isomerase (CHI), flavonoid 3'5'-hydroxylase (F3'5'H), dihydroflavonol 4-reductase (DFR), anthocyanidin synthase (ANS), and anthocyanidin 3-O-glucosyltransferase (UFGT) were highly similar and significantly positively correlated with anthocyanin accumulation rate in BF. In particular, F3'5'H, UFGT, ANS, and DFR expression levels in BF were 2391, 119, 96, and 85 times higher than those in WF at S3 (35 days after anthesis), respectively. This result strongly suggests that the low expression of these genes in WF is responsible for the anthocyanin loss. Meanwhile, the expression patterns of the anthocyanin regulatory genes were also investigated by qRT-PCR. Mass sequencing data were obtained and annotated by deep sequencing and provided a platform for future function and molecular biological research on L. ruthenicum Murr.

makuwa)

Rootstocks are among the primary factors that influence fruit yield and quality as well as melon development. To understand the differences in the molecular mechanisms and gene expression networks of fruit development between grafted and nongrafted plants in oriental melon, we performed a comprehensive analysis of the transcriptome and proteome dynamic gene/protein expression profiles during fruit development in oriental melon (Cucumis melo L. var. makuwa). Using pairwise comparisons between grafted and nongrafted samples by transcriptome analysis, we identified a large number of candidate genes involved in hormonal signaling pathways, transcription factors, resistance-related biosynthetic pathways and photosynthesis-related metabolic pathways. Many transcription factor-encoded genes were significantly more strongly expressed in the grafted samples, for example, AP2/ERF, C2H2, MYB, bHLH, and AUX/IAA, which are well-known participants in the regulation of developmental processes and hormonal signaling metabolism. Some differentially expressed genes (DEGs) were enriched in flavonoid biosynthesis and phenylpropanoid biosynthesis and determined plant resistance. In addition, some differentially expressed proteins (DEPs) were enriched in photosynthesis-related pathways, which could improve fruit quality and yield. Moreover, through weighted gene coexpression network analyses, we identified modules of coexpressed genes and hub genes specifically related to grafting for different fruit developmental stages. The results suggested that graft-related modules and hub genes were primarily associated with photosynthate metabolism and hormonal signaling pathways. The results obtained in this study provide a valuable resource for dissecting the role of candidate genes governing graft-related metabolism in oriental melon fruit, suggesting an interesting correlation with the effects of rootstock on fruit development.

Anthocyanin accumulation is initiated by abscisic acid to enhance fruit color during fig (Ficus carica L.) ripening

Fig fruit is well-known for its attractive flavor, color, and nutritional and medicinal value. Anthocyanin contributes to the fruit's color and constitutes a high percentage of the total antioxidant content of the fig fruit. We quantified the major anthocyanins and characterized the expression levels of anthocyanin-biosynthesis and transcription factor genes in fruit treated on-tree with exogenous abscisic acid (ABA) or ethephon, or the ABA inhibitors nordihydroguaiaretic acid (NDGA) or fluridone. The major anthocyanins cyanidin 3-O-glucoside and cyanidin 3-O-rutinoside were found in significantly higher quantities in exogenous ABA- and ethephon-treated fruit, with early dark purple color compared to the controls. On the other hand, NDGA- and fluridone-treated fruit had significantly lower amounts of anthocyanins, with less purple color coverage than controls. Expression levels of the anthocyanin-biosynthesis genes FcPAL, FcCHS2, FcCHI, FcF3H, FcDFR, FcANS, FcUFGT and Fc3RT were upregulated by exogenous ABA and ethephon treatment, and downregulated by NDGA and fluridone treatment. The MYB-bHLH-WD40 complex-related genes of ripe fig fruit were identified. In particular, FcMYB113 was strongly upregulated by exogenous ABA and ethephon, and strongly downregulated by NDGA and fluridone. In addition, moderate upregulation of FcGL3 and FcWD40 was observed with exogenous ABA and ethephon treatment, and moderate downregulation in NDGA- and fluridone-treated fruit. These results indicate that ABA can initiate anthocyanin biosynthesis, which ultimately improves the color and nutritional value of fig fruit, enhancing their attractiveness to consumers.

Genome-wide analysis of the bZIP gene family in Chinese jujube (Ziziphus jujuba Mill.)

BACKGROUND

Among several TF families unique to eukaryotes, the basic leucine zipper (bZIP) family is one of the most important. Chinese jujube (Ziziphus jujuba Mill.) is a popular fruit tree species in Asia, and its fruits are rich in sugar, vitamin C and so on. Analysis of the bZIP gene family of jujube has not yet been reported. In this study, ZjbZIPs were identified firstly, their expression patterns were further studied in different tissues and in response to various abiotic and phytoplasma stresses, and their protein-protein interactions were also analyzed.

RESULTS

At the whole genome level, 45 ZjbZIPs were identified and classified into 14 classes. The members of each class of bZIP subfamily contain a specific conserved domain in addition to the core bZIP conserved domain, which may be related to its biological function. Relative Synonymous Codon Usage (RSCU) analysis displayed low values of NTA and NCG codons in ZjbZIPs, which would be beneficial to increase the protein production and also indicated that ZjbZIPs were at a relative high methylation level. The paralogous and orthologous events occurred during the evolutionary process of ZjbZIPs. Thirty-four ZjbZIPs were mapped to but not evenly distributed among 10 pseudo- chromosomes. 30 of ZjbZIP genes showed diverse tissue-specific expression in jujube and wild jujube trees, indicating that these genes may have multiple functions. Some ZjbZIP genes were specifically analyzed and found to play important roles in the early stage of fruit development. Moreover, some ZjbZIPs that respond to phytoplasma invasion and abiotic stress environmental conditions, such as salt and low temperature, were found. Based on homology comparisons, prediction analysis and yeast two-hybrid, a protein interaction network including 42 ZjbZIPs was constructed.

CONCLUSIONS

The bioinformatics analyses of 45 ZjbZIPs were implemented systematically, and their expression profiles in jujube and wild jujube showed that many genes might play crucial roles during fruit ripening and in the response to phytoplasma and abiotic stresses. The protein interaction networks among ZjbZIPs could provide useful information for further functional studies.

Dynamics of ascorbic acid content in apple (Malus x domestica) during fruit development and storage

Vitamin C is a crucial antioxidant and cofactor for both plants and humans. Apple fruits generally contain low levels of vitamin C, making vitamin C content an interesting trait for apple crop improvement. With the aim of breeding high vitamin C apple cultivars it is important to get an insight in the natural biodiversity of vitamin C content in apple fruits. In this study, quantification of ascorbic acid (AsA), dehydroascorbic acid (DHA), and total AsA (AsA + DHA) in apple pulp of 79 apple accessions at harvest revealed significant variation, indicating a large genetic biodiversity. High density genotyping using an 8 K SNP array identified 21 elite and 58 local cultivars in this germplasm, with local accessions showing similar levels of total AsA but higher amounts of DHA compared to elite varieties. Out of the 79 apple cultivars screened, ten genotypes with either the highest or the lowest concentration of total AsA at harvest were used for monitoring vitamin C dynamics during fruit development and storage. For all these cultivars, the AsA/DHA ratio in both apple pulp and peel increased throughout fruit development, whereas the AsA/DHA balance always shifted towards the oxidized form during storage and shelf life, putatively reflecting an abiotic stress response. Importantly, at any point during apple fruit development and storage, the apple peel contained a higher level of vitamin C compared to the pulp, most likely because of its direct exposure to abiotic and biotic stresses.

Genome-wide characterization and expression analysis of the auxin response factor (ARF) gene family during melon (Cucumis melo L.) fruit development

ARFs in plants mediate auxin signaling transduction and regulate growth process. To determine genome-wide characterization of ARFs family in melon (Cucumis melo L.), ARFs were identified via analysis of information within the melon genomic database, and bioinformatic analyses were performed using various types of software. Based on different treatment methods involving dipping with the growth regulator Fengchanji No. 2 and artificial pollination, Jingmi No. 11 melon was used as the test material, and melon plants with unpollinated ovaries served as controls. The expression of ARFs during the early development of melon was analyzed via qRT-PCR. Seventeen genes that encode ARF proteins were identified in the melon genome for the first time. The expression of these ARFs differed in different tissues. The expression levels of CmARF2, CmARF16-like, CmARF18-like2, and CmARF19-like were especially high in melon fruits. The expression of ARFs during the early development of melon fruits differed in response to the different treatments, which suggested that CmARF9, CmARF16-like, CmARF19-like, CmARF19, CmARF1, CmARF2, CmARF3, and CmARF5 may be associated with melon fruit growth during early development. Interestingly, the increase in the transverse diameter of fruits treated with growth regulators was significantly greater than that of fruits resulting from artificial pollination, while the increase in the longitudinal diameter of the fruits resulting from artificial pollination was significantly greater.

'Fujiminori')

The fruit is the most important economical organ in the grape; accordingly, to investigate the grapevine genomic methylation landscape and examine its functional significance during fruit development, we generated whole genome DNA methylation maps for various developmental stages in the fruit of grapevine. In this study, thirteen DNA methylation-related genes and their expression profiles were identified and analyzed. The methylation levels for mC, mCG, mCHG, and mCHH contexts in 65 days after flowering (65DAF) fruit (véraison stage) were higher than those in 40DAF (green stage) and 90DAF (mature stage) fruits. Relative to methylation in the mC context, methylation levels in the mCHH context were higher than those of mCG and mCHG. The DNA methylation level in the ncRNA regions was significantly higher than that in exon, gene, intron, and mRNA regions. The differentially methylated regions (DMRs) and differentially methylated promoters (DMPs) in 65DAF_vs_40DAF were both higher than those in 90DAF_vs_65DAF and 90DAF_vs_40DAF. Most DMRs (or DMPs) were involved in metabolic processes and cell processes, binding, and catalytic activity. These results indicated that DNA methylation represses gene expression during grape fruit development, and it broadens our understanding of the landscape and function of DNA methylation in grapevine genomes.

Comprehensive transcriptome analysis of reference genes for fruit development of Euscaphis konishii

Background

Quantitativereal-time reverse transcriptase polymerase chain reaction is the common method to quantify relative gene expression. Normalizating using reliable genes is critical in correctly interpreting expression data from qRT-PCR. Euscaphis konishii is a medicinal plant with a long history in China, which has various chemical compounds in fruit. However, there is no report describing the selection of reference genes in fruit development of Euscaphis konishii.

Methods

We selected eight candidate reference genes based on RNA-seq database analysis, and ranked expression stability using statistical algorithms GeNorm, NormFinder, BestKeeper and ReFinder. Finally, The nine genes related to the anthocyanin synthesis pathway of Euscaphis konishii were used to verify the suitability of reference gene.

Results

The results showed that the stability of EkUBC23, EkCYP38 and EkGAPDH2 was better, and the low expression reference genes (EkUBC23 and EkCYP38) were favourable for quantifying low expression target genes, while the high expression reference gene (EkGAPDH2) was beneficial for quantifying high expression genes. In this study, we present the suitable reference genes for fruit development of Euscaphis konishii based on transcriptome data, our study will contribute to further studies in molecular biology and gene function on Euscaphis konishii and other closely related species.

Hormonal and transcriptional analyses of fruit development and ripening in different varieties of black pepper (Piper nigrum)

Black pepper (Piper nigrum L.) is one of the most popular and oldest spices in the world with culinary uses and various pharmacological properties. In order to satisfy the growing worldwide demand for black pepper, improved productivity of pepper is highly desirable. A primary constraint in black pepper production is the non-synchronous nature of flower development and non-uniform fruit ripening within a spike. The uneven ripening of pepper berries results in a high labour requirement for selective harvesting contributes to low productivity and affects the quality of the pepper products. In Malaysia, there are a few recommended varieties for black pepper planting, each having some limitations in addition to the useful characteristics. Therefore, a comparative study of different black pepper varieties will provide a better understanding of the mechanisms regulates fruit development and ripening. Plant hormones are known to influence the fruit development process and their roles in black pepper flower and fruit development were inferred based on the probe-based gene expression analysis and the quantification of the multiple plant hormones using high-performance liquid chromatography coupled to tandem mass spectrometry (HPLC-MS/MS). In this study, jasmonic acid and salicylic acid were found to play roles in flowering and fruit setting, whereas auxin, gibberellin and cytokinins are important for fruit growth. Abscisic acid has positive role in fruit maturation and ripening in the development process. Distinct pattern of plant hormones related gene expression profiles with the hormones accumulation profiles suggested a complex network of regulation is involved in the signaling process and crosstalk between plant hormones was another layer of regulation in the black pepper fruit development mechanisms. The current study provides clues to help in elucidating the timing of the action of each specific plant hormone during fruit development and ripening which could be applied to enhance our ability to control the ripening process, leading to improving procedures for the production and post-harvest handling of pepper fruits.

The power of model-to-crop translation illustrated by reducing seed loss from pod shatter in oilseed rape

Elucidation of key regulators in Arabidopsis fruit patterning has facilitated knowledge-translation into crop species to address yield loss caused by premature seed dispersal (pod shatter). In the 1980s, plant scientists descended on a small weed Arabidopsis thaliana (thale cress) and developed it into a powerful model system to study plant biology. The massive advances in genetics and genomics since then have allowed us to obtain incredibly detailed knowledge on specific biological processes of Arabidopsis growth and development, its genome sequence and the function of many of the individual genes. This wealth of information provides immense potential for translation into crops to improve their performance and address issues of global importance such as food security. Here, we describe how fundamental insight into the genetic mechanism by which seed dispersal occurs in members of the Brassicaceae family can be exploited to reduce seed loss in oilseed rape (Brassica napus). We demonstrate that by exploiting data on gene function in model species, it is possible to adjust the pod-opening process in oilseed rape, thereby significantly increasing yield. Specifically, we identified mutations in multiple paralogues of the INDEHISCENT and GA4 genes in B. napus and have overcome genetic redundancy by combining mutant alleles. Finally, we present novel software for the analysis of pod shatter data that is applicable to any crop for which seed dispersal is a serious problem. These findings highlight the tremendous potential of fundamental research in guiding strategies for crop improvement.

The nutrient distribution in the continuum of the pericarp, seed coat, and kernel during Styrax tonkinensis fruit development

Background

Styrax tonkinensis is a great potential biofuel as the species contains seeds with a particularly high oil content. Understanding the nutrient distribution in different parts of the fruit is imperative for the development and enhancement of S. tonkinensis as a biodiesel feedstock.

Methods

From 30 to 140 days after flowering (DAF), the development of S. tonkinensis fruit was tracked. The morphology change, nutrient content, and activity of associated enzymes in the continuum of the pericarp, seed coat, and kernel were analyzed.

Results

Between 30 and 70 DAF, the main locus of dry matter deposition shifted from the seed coat to the kernel. The water content within the pericarp remained high throughout development, but at the end (130 DAF later) decreased rapidly. The water content within both the seed coat and the kernel consistently declined over the course of the fruit development (30–110 DAF). Between 70 and 80 DAF, the deposition centers for sugar, starch, protein, potassium, and magnesium was transferred to the kernel from either the pericarp or the seed coat. The calcium deposition center was transferred first from pericarp to the seed coat and then to the kernel before it was returned to the pericarp. The sucrose to hexose ratio in the seed coat increased between 30 and 80 DAF, correlating with the accumulation of total soluble sugar, starch, and protein. In the pericarp, the sucrose to hexose ratio peaked at 40 and 100 DAF, correlating with the reserve deposition in the following 20–30 days. After 30 DAF, the chlorophyll concentration of both the pericarp and the seed coat dropped. The maternal unit (the pericarp and the seed coat) in fruit showed a significant positive linear relationship between chlorophyll b/a and the concentration of total soluble sugar. The potassium content had significant positive correlation with starch (ρ = 0.673, p = 0.0164), oil (ρ = 0.915, p = 0.000203), and protein content (ρ = 0.814, p = 0.00128), respectively. The concentration of magnesium had significant positive correlation with starch (ρ = 0.705, p = 0.0104), oil (ρ = 0.913, p = 0.000228), and protein content (ρ = 0.896, p = 0.0000786), respectively. Calcium content had a significant correlation with soluble sugar content (ρ = 0.585, p = 0.0457).

Conclusions

During the fruit development of S. tonkinensis, the maternal unit, that is, the pericarp and seed coat, may act a nutrient buffer storage area between the mother tree and the kernel. The stage of 70–80 DAF is an important time in the nutrient distribution in the continuum of the pericarp, seed coat, and kernel. Our results described the metabolic dynamics of the continuum of the pericarp, seed coat, and kernel and the contribution that a seed with high oil content offers to biofuel.

Genome-wide identification and analysis of the MADS-box gene family and its potential role in fruit development and ripening in red bayberry (Morella rubra)

The MADS-box gene family encodes transcription factors and plays an important role in plant growth and the development of flower and fruit. A perennial dioecious plant, the red bayberry genome has been published recently, providing the opportunity to analyze the MADS-box gene family and its role in fruit development and ripening. Here, we identified 54 MADS-box genes in the red bayberry genome, and classified them into two types based on phylogenetic analysis. Thirteen Type I MADS-box genes were subdivided into three subfamilies and 41 Type II MADS-box genes into 13 subfamilies. A total of 46 MADS-box genes were distributed across eight red bayberry chromosomes, and the other eight genes were located on the unmapped scaffolds. Transcriptome analysis suggested that the expression of most Type II genes was higher than Type I in five female tissues. Moreover, 26 MADS-box genes were expressed during red bayberry fruit development and ten of them showed high expression. qRT-PCR showed that the expression of MrMADS01 (SEP, MIKC^C), with differences between the pale pink and red varieties, increased significantly at the final ripening stage, suggesting it may participate in ripening as positive regulator and related to anthocyanin biosynthesis. These results provide some clues for future study of MADS-box genes in red bayberry, especially in ripening process.

Genome-wide identification, expression analysis and functional study of the GRAS gene family in Tartary buckwheat (Fagopyrum tataricum)

BACKGROUND

GRAS are plant-specific transcription factors that play important roles in plant growth and development. Although the GRAS gene family has been studied in many plants, there has been little research on the GRAS genes of Tartary buckwheat (Fagopyrum tataricum), which is an important crop rich in rutin. The recently published whole genome sequence of Tartary buckwheat allows us to study the characteristics and expression patterns of the GRAS gene family in Tartary buckwheat at the genome-wide level.

RESULTS

In this study, 47 GRAS genes of Tartary buckwheat were identified and divided into 10 subfamilies: LISCL, HAM, DELLA, SCR, PAT1, SCL4/7, LAS, SHR, SCL3, and DLT. FtGRAS genes were unevenly distributed on 8 chromosomes, and members of the same subfamily contained similar gene structures and motif compositions. Some FtGRAS genes may have been produced by gene duplications; tandem duplication contributed more to the expansion of the GRAS gene family in Tartary buckwheat. Real-time PCR showed that the transcription levels of FtGRAS were significantly different in different tissues and fruit development stages, implying that FtGRAS might have different functions. Furthermore, an increase in fruit weight was induced by exogenous paclobutrazol, and the transcription level of the DELLA subfamily member FtGRAS22 was significantly upregulated during the whole fruit development stage. Therefore, FtGRAS22 may be a potential target for molecular breeding or genetic editing.

CONCLUSIONS

Collectively, this systematic analysis lays a foundation for further study of the functional characteristics of GRAS genes and for the improvement of Tartary buckwheat crops.

Genome-wide identification of the SPL gene family in Tartary Buckwheat (Fagopyrum tataricum) and expression analysis during fruit development stages

BACKGROUND

SPL (SQUAMOSA promoter binding protein-like) is a class of plant-specific transcription factors that play important roles in many growth and developmental processes, including shoot and inflorescence branching, embryonic development, signal transduction, leaf initiation, phase transition, and flower and fruit development. The SPL gene family has been identified and characterized in many species but has not been well studied in tartary buckwheat, which is an important edible and medicinal crop.

RESULTS

In this study, 24 Fagopyrum tataricum SPL (FtSPL) genes were identified and renamed according to the chromosomal distribution of the FtSPL genes. According to the amino acid sequence of the SBP domain and gene structure, the SPL genes were divided into eight groups (group I to group VII) by phylogenetic tree analysis. A total of 10 motifs were detected in the tartary buckwheat SPL genes. The expression patterns of 23 SPL genes in different tissues and fruits at different developmental stages (green fruit stage, discoloration stage and initial maturity stage) were determined by quantitative real-time polymerase chain reaction (qRT-PCR).

CONCLUSIONS

The tartary buckwheat genome contained 24 SPL genes, and most of the genes were expressed in different tissues. qRT-PCR showed that FtSPLs played important roles in the growth and development of tartary buckwheat, and genes that might regulate flower and fruit development were preliminarily identified. This work provides a comprehensive understanding of the SBP-box gene family in tartary buckwheat and lays a significant foundation for further studies on the functional characteristics of FtSPL genes and improvement of tartary buckwheat crops.

Genome-wide investigation of the ZF-HD gene family in Tartary buckwheat (Fagopyrum tataricum)

BACKGROUND

ZF-HD is a family of genes that play an important role in plant growth, development, some studies have found that after overexpression AtZHD1 in Arabidopsis thaliana, florescence advance, the seeds get bigger and the life span of seeds is prolonged, moreover, ZF-HD genes are also participate in responding to adversity stress. The whole genome of the ZF-HD gene family has been studied in several model plants, such as Arabidopsis thaliana and rice. However, there has been little research on the ZF-HD genes in Tartary buckwheat (Fagopyrum tataricum), which is an important edible and medicinal crop. The recently published whole genome sequence of Tartary buckwheat allows us to study the tissue and expression profiles of the ZF-HD gene family in Tartary buckwheat on a genome-wide basis.

RESULTS

In this study, the whole genome and expression profile of the ZF-HD gene family were analyzed for the first time in Tartary buckwheat. We identified 20 FtZF-HD genes and divided them into MIF and ZHD subfamilies according to phylogeny. The ZHD genes were divided into 5 subfamilies. Twenty FtZF-HD genes were distributed on 7 chromosomes, and almost all the genes had no introns. We detected seven pairs of chromosomes with fragment repeats, but no tandem repeats were detected. In different tissues and at different fruit development stages, the FtZF-HD genes obtained by a real-time quantitative PCR analysis showed obvious expression patterns.

CONCLUSIONS

In this study, 20 FtZF-HD genes were identified in Tartary buckwheat, and the structures, evolution and expression patterns of the proteins were studied. Our findings provide a valuable basis for further analysis of the biological function of the ZF-HD gene family. Our study also laid a foundation for the improvement of Tartary buckwheat crops.

Cytokinin treatment modifies litchi fruit pericarp anatomy leading to reduced susceptibility to post-harvest pericarp browning

Litchi (Litchi chinensis Sonn.) is a subtropical fruit known for its attractive red pericarp color, semi-translucent white aril and unique flavor and aroma. Rapid post-harvest pericarp browning strictly limits litchi fruit marketing. In the current research, we hypothesized that modification of litchi fruit pericarp anatomy by hormone application may reduce fruit susceptibility to post-harvest pericarp browning. In this context, we hypothesized that cytokinin treatment, known to induce cell division, may yield fruit with thicker pericarp and reduced susceptibility for fruit surface micro-crack formation, water loss and post-harvest pericarp browning. Exogenous cytokinin treatment was applied at different stages along the course of litchi fruit development and the effect on fruit pericarp anatomy, fruit maturation and postharvest pericarp browning was investigated. Interestingly, cytokinin treatment, applied 4 weeks after full female bloom (WFB), during the phase of pericarp cell division, led to mature fruit with thicker pericarp, reduced rate of post-harvest water loss and reduced susceptibility to post-harvest pericarp browning, as compared to non-treated control fruit. Histological sections ascribe the difference in pericarp anatomy to increased cell proliferation in the parenchymatic tissue and the highly-lignified brachysclereid cell layer. In contrast, exogenous cytokinin treatment applied 7 WFB, following the phase of pericarp cell division, significantly increased epidermal-cell proliferation but had no significant effect on overall fruit pericarp thickness and only minor affect on post-harvest water loss or pericarp browning. Interestingly, the late cytokinin treatment also significantly postponed fruit maturation-associated anthocyanin accumulation and chlorophyll degradation, as previously reported, but had no effect on other parameters of fruit maturation, like total soluble sugars and total titratable acids typically modified during aril maturation. In conclusion, exogenous cytokinin treatment at different stages in fruit development differentially modifies litchi fruit pericarp anatomy by induction of cell-type specific cell proliferation. Early cytokinin treatment during the phase of pericarp cell division may prolong litchi fruit storage by reducing fruit susceptibility to post-harvest water loss and pericarp browning.