- Open Access
Integrative transcriptome and proteome analyses provide new insights into different stages of Akebia trifoliata fruit cracking during ripening
Biotechnology for Biofuels volume 13, Article number: 149 (2020)
Akebia trifoliata (Thunb.) Koidz may have applications as a new potential source of biofuels owing to its high seed count, seed oil content, and in-field yields. However, the pericarp of A. trifoliata cracks longitudinally during fruit ripening, which increases the incidence of pests and diseases and can lead to fruit decay and deterioration, resulting in significant losses in yield. Few studies have evaluated the mechanisms underlying A. trifoliata fruit cracking.
In this study, by observing the cell wall structure of the pericarp, we found that the cell wall became thinner and looser and showed substantial breakdown in the pericarp of cracking fruit compared with that in non-cracking fruit. Moreover, integrative analyses of transcriptome and proteome profiles at different stages of fruit ripening demonstrated changes in the expression of various genes and proteins after cracking. Furthermore, the mRNA levels of 20 differentially expressed genes were analyzed, and parallel reaction monitoring analysis of 20 differentially expressed proteins involved in cell wall metabolism was conducted. Among the molecular targets, pectate lyases and pectinesterase, which are involved in pentose and glucuronate interconversion, and β-galactosidase 2, which is involved in galactose metabolism, were significantly upregulated in cracking fruits than in non-cracking fruits. This suggested that they might play crucial roles in A. trifoliata fruit cracking.
Our findings provided new insights into potential genes influencing the fruit cracking trait in A. trifoliata and established a basis for further research on the breeding of cracking-resistant varieties to increase seed yields for biorefineries.
Akebia trifoliata (Thunb.) Koidz is a perennial, wild, woody vine, belonging to the family Lardizabalaceae and the subgenus Akebia Decne . This species is characterized by numerous seeds (300–800), high seed yield (400 kg/667 m2), and high seed oil contents (47.62%) [2,3,4]. The physicochemical properties of seeds, such as water content, iodine value, acid value, peroxide value, cetane number, flash point, and cold filter plugging point of the oil, meet the biodiesel standard, implying applications in biodiesel production for the development and utilization of liana oil crops [3, 5]. Similar in structure to petrol-diesel, biodiesel prepared with C15–C18 alkane fractions with a high cetane number of 91 can be obtained from A. trifoliata seed oil. This cetane number is similar to that of biodiesel prepared using Sapium sebiferum oil (40.2) and Vernicia fordii (53) [6,7,8]. Moreover, the fruits of A. trifoliata have high sugar content, which is important in industrial bioethanol production because it can boost the amount of ethanol produced at the end of fermentation [9, 10]. Enzymatic saccharification and ethanol yields of 98.0% and 100%, respectively, were achieved from the solid residues of Akebia after CHCl-formic acid pretreatment, suggesting potential applications as a feedstock for cellulosic ethanol production . The advantages of A. trifoliata, including high seed count, high seed oil contents, abundant yields from the field, high total sugar contents of fruit, wide adaptability, tolerance to both drought and heat, and ease of management, could provide new resources for biodiesels and feedstocks for cellulosic ethanol production [4, 5, 11]. A. trifoliata can also be exploited as a new high-value fruit crop owing to its high nutritional, ornamental, economic, medicinal, and potential development value [12, 13].
However, the pericarp of A. trifoliata cracks longitudinally along the ventral suture when matured. Fruit cracking is a serious problem that increases the incidence of pests and diseases, leading to fruit decay and deterioration, affecting the utilization rate of fruit and seeds and causing significant losses in yields and commercial value [14, 15]. Studies have indicated that the dry mass and oil contents of insect-infested fruits (33.1 mg/fruit and 3.9%, respectively) are lower than those of healthy fruits (67.4 mg/fruit and 39.9%, respectively); moreover, the observed fruit and oil yields (2.9 and 0.6 kg/tree, respectively) are lower than the excepted yields (4.7 and 1.9 kg/tree, respectively) . Fruit ripening and cracking is a complex, genetically programmed process that is accompanied by the synthesis of large amounts of proteins and the transcription of many genes [17, 18]. Studies have indicated that different genetic accessions exhibit major differences in cracking resistance, suggesting that genetic factors play important roles in sweet cherry fruit cracking resistance . Studies on fruit cracking have been performed in many species, including tomatoes, litchis, durians, and apples; these studies have suggested that cell wall-modifying proteins, such as polygalacturonases (PGs), pectinesterase (PE), β-galactosidases (β-GALs), expansins (EXPs), and xyloglucan endotransglycosylase proteins [14, 20,21,22], are associated with fruit cracking. The basic helix-loop-helix (bHLH) gene INDEHISCENT regulates Lepidium campestre fruit dehiscence . Additionally, Dong et al.  found that pod shattering resistance in soybeans was mediated by the NAC gene. Sorefan et al.  indicated that a regulated auxin minimum was required for seed dispersal in Arabidopsis. Consequently, a better understanding of the genetic analysis of fruit cracking is necessary to prevent cracking phenotypes.
Although many studies have evaluated fruit cracking, little progress has been made with regard to our understanding of the molecular mechanisms underlying A. trifoliata fruit cracking, and only two studies have performed transcriptome sequencing in this species [1, 26]. Our limited knowledge of the molecular characteristics of A. trifoliata has made it difficult to recommend preventive measures for fruit cracking. Next-generation sequencing methods, such as transcriptome and proteome technologies for measuring gene expression and protein abundance, have become powerful tools for the discovery of novel genes and their functions in regulating fruit ripening and cracking [27, 28].
Therefore, in this study, an integrative analysis of the transcriptome and proteome was performed to enhance our understanding of A. trifoliata fruit cracking during ripening at the molecular level. Our data provide important insights into omics resources and candidate genes responsible for fruit cracking traits in A. trifoliata.
Changes in pericarp structure
First, we evaluated the dynamic structures of fruit pericarps during different stages of ripening (Fig. 1). In the non-cracking stage (PS), the pericarp cells and cuticles were densely arranged, with small intercellular spaces and continuous distribution (Fig. 1a, d, g). However, in the initial cracking stage, the cell wall thinned, the cell volume increased, the number of cell layers decreased, and the cells were loosely arranged with poor integrity; moreover, the spacing between the cells increased, and the cell of the exocarp and mesocarp began to degrade in the initial cracking stage (PM; Fig. 1b, e, h). Irregularly arranged layers with continuously reduced numbers and larger spaces were observed as an evidence of cell degradation throughout the cracking stages (PL; Fig. 1c, f, i).
Transcriptomic analysis overview
Nine cDNA libraries were constructed, and 47.05, 46.92, and 54.00 million raw sequence reads were generated from the PS, PM, and PL libraries, respectively. After removing low-quality reads and adaptor sequences, 46.45, 46.35, and 53.46 million clean reads with 91.58–93.75% Q30 bases and 45.94–48.48% GC content, respectively, were obtained (Additional file 1: Table S1). The resulting A. trifoliata transcriptome contained 241,376 transcripts, ranging from 201 to 2000 bp, and 186,054 unigenes (> 200 bp; Table 1 and Additional file 1: Fig. S1).
All unigenes were annotated using Basic Local Alignment Search Tool (BLAST) searches against the following five databases: National Center for Biotechnology Information nonredundant protein sequences database (NR; 100,329; 53.9% and 41.57% of all identified unigenes and transcripts, respectively), SwissProt (56,346; 30.3% and 23.3%), Protein Families database (Pfam; 34,428; 18.5% and 14.3%), Gene Ontology database (GO; 44,558; 23.9% and 18.5%), and Kyoto Encyclopedia of Genes and Genomes pathway database (KEGG; 30,298; 16.3% and 12.6%). This indicated that the NR database provided the largest number of annotations, suggesting that 100,329 unigenes corresponded with sequences from at least one of the public databases, and 7283 unigenes were annotated to all databases. Moreover, 17,601, 19,281, and 13,525 unique unigenes and 45,866, 48,104, and 43,303 absent unigenes were identified in PS, PM, and PL, respectively. Most of these were uncharacterized, unknown, and hypothetical proteins in the annotation results (Additional file 1: Tables S2–S3).
Among these unigenes, 9301 were identified as differently expressed genes (DEGs), including 2703 (779 up- and 1924 downregulated), 4694 (3815 up- and 879 downregulated) and 1904 co-expressed in PM versus PS and PL versus PM groups, respectively (Table 2), which are presented in a volcano plot (Additional file 1: Fig. S2a–b). The summary information of these DEGs is shown in Additional file 1: Table S4.
Functional classification of identified DEGs
Bioinformatics analysis indicated that most DEGs in the GO terms of biological process (BP) were involved in cellular amide metabolic processes and amide metabolic processes. Genes related to structural molecular activity and oxidoreductase activity accounted for the highest proportions of DEGs in the molecular function (MF) category in both the PM_PS and PL_PM groups. Cytoplasmic parts and intracellular ribonucleoprotein complexes comprised the highest proportions of DEGs in the cell components (CCs) in the PM_PS and PL_PM groups, respectively (Fig. 2a, b).
KEGG pathway analysis indicated that many DEGs were enriched in metabolic pathways and ribosomes in the PM_PS comparison and in the biosynthesis of secondary metabolites and metabolic pathways in the PL_PM comparison (Fig. 2e, f). Cluster analysis showed that cell wall-related DEGs (285), including pentose and glucuronate interconversion, the phenylpropanoid pathway, galactose metabolism, starch and sucrose metabolism, amino sugar and nucleotide sugar metabolism, and transcription factors, were clustered closely in the PM_PS and PL_PM groups. Most DEGs were downregulated in the PM_PS group, but these were upregulated in the PL_PM group (Fig. 3).
Quantitative proteome analysis
In total, 812,625 spectra, 68,151 identified spectra, 12,456 peptides, 10,572 unique peptides, and 2839 proteins were determined via proteomic analysis (Table 1 and Additional file 1: Table S3). In terms of protein mass distribution, proteins with molecular weights greater than 9 kDa had a wide range and good coverage, with a maximum distribution area of 10–40 kDa. Peptide quantitative analysis of the proteins showed that protein quantity decreased with an increase in matching peptides (Additional file 1: Fig. S3).
Among these 2839 proteins, 223 were identified as differentially abundant proteins (DAPs), including 173 (75 up- and 98 downregulated), 33 (20 up- and 13 downregulated), and 17 co-expressed proteins in the PM_PS and PL_PM groups, respectively (Table 2), which are presented in a volcano plot (Additional file 1: Fig. S2c, d). The summary information of these DAPs is shown in Additional file 1: Table S4.
Functional classification of the identified DAPs
GO analysis showed that most DAPs in the BP category were involved in cellular responses to chemical stimulus and cellular oxidant detoxification processes in the PM_PS group and in metabolic processes and macromolecule metabolic processes in the PL_PM group. The highest proportions of DAPs in the MF category were involved in oxidoreductase activity and antioxidant activity in the PM_PS group and structural constituent of ribosome and structural molecule activity in the PL_PM group. The extracellular region in the PM_PS group and cytoplasmic parts in the PL _PM group showed the highest portions of DAPs in the CC category (Fig. 2c, d). KEGG pathway analyses indicated that many proteins were enriched in the two-component system and ribosome pathways in the PM_PS and PL_PM groups, respectively (Fig. 2g, h).
Hierarchical cluster analysis showed that 40 cell wall-related DAPs, including pentose and glucuronate interconversions, the phenylpropanoid pathway, galactose metabolism, starch and sucrose metabolism, amino sugar and nucleotide sugar metabolism, and other cell wall metabolism-related proteins, were clustered closely in the PM_PS and PL_PM groups. Notably, most DAPs were upregulated in both the PM_PS and PL_PM groups. Alternatively, DAPs involved in phenylpropanoid pathways and galactose metabolism were downregulated in the PM_PS group but upregulated in the PL_PM group (Fig. 3g–l). Moreover, the protein–protein interaction (PPI) network analysis indicated that 85 DAPs, including 29 upregulated (28 and 1) and 56 downregulated DAPs (46 and 10), were involved in the interaction networks in the PM_PS and PL_PM groups, respectively (Fig. 4).
Comparative analysis between protein abundance and gene expression
There were more DEGs (11,205) than DAPs (240) and shared DEGs (1904) than shared DAPs (17) in cracked fruit than in non-cracking fruit. Most were downregulated in the PM_PS group but upregulated in the PL_PM group (Table 2). Moreover, integration of the proteome and transcriptome data showed that 14 and 4 DAPs were matched with their DEGs. In addition, 12 (4 up- and 8 downregulated) and 4 DAPs (2 up- and 2 downregulated) showed the same tendency as DEGs in the PM_PS and PL_PM groups, respectively, with two showing the opposite tendency as DEGs in the PM_PS group (Fig. 5a, b). Furthermore, the fold-changes in DAPs indicated differentially positive correlations with their corresponding DEGs based on Pearson’s correlation tests. A limited correlation (r = 0.03 and 0.11) was detected between the proteome and transcriptome, and a relatively higher positive correlation was identified with the same trend (r = 0.9161 and 0.8) for DEGs and DAPs in the PM_PS and PL_PM groups, respectively (Fig. 5c–f).
Identification of DAPs and DEGs associated with candidate pathways
The correlation of proteome and transcriptome GO enrichment indicated that the largest groups of those DEGs and DAPs within the BP, MF, and CC categories were linked to metabolic and cellular processes, catalytic activity and binding, and cells and cell parts, respectively, in both the PM_PS and PL_PM groups (Additional file 1: Fig. S4a–b). KEGG enrichment showed that most were significantly enriched in phenylpropanoid biosynthesis, pentose and glucuronate interconversions, amino sugar and nucleotide sugar metabolism, and galactose metabolism in the PM_PS group. Most were significantly enriched in pentose and glucuronate interconversions and the galactose metabolism pathway in the PL_PM group for both the proteome and transcriptome (Additional file 1: Fig S4c–d). Comparative analysis showed that two pathways, i.e., pentose and glucuronate interconversion and galactose metabolism pathways, were shared in both the PM_PS and PL_PM groups. In contrast, the phenylpropanoid biosynthesis pathway was only shared in the PM_PS group.
Validation of data reliability through reverse transcription real-time quantitative PCR (qPCR) and parallel reaction monitoring (PRM)
qPCR experiments were performed for 20 selected DEGs, and the results are shown in Fig. 6. In the PM_PS group, the phenylpropanoid pathway-related genes 4-coumarate-COA-ligase (4CL), peroxidase (PRX), and PRX2 were downregulated, whereas cinnamyl-alcohol dehydrogenase (CAD) and shikimate O-hydroxycinnamoyltransferase (HCT) were upregulated. The galactose metabolism-related genes β-GAL1 and β-GAL2; amino sugar and nucleotide sugar metabolism-related gene beta-d-xylosidase (BXL); and starch and sucrose metabolism-related genes cellulase (CEL), cellulose synthase-like protein (CSLG), and glucan endo-1,3-beta-d-glucosidase (ENDOB) were downregulated. The cell wall metabolism genes NAC, NAC-like, and EXP1 were downregulated, whereas the BHLH transcription factor and dirigent protein (DIR2) were upregulated. The pentose and glucuronate interconversion-related genes PL, PG, and PE were upregulated. Most of these genes, except 4CL, CAD, β-GAL, and EXP1, were significantly upregulated in the PL_PM group. Moreover, the expression of 13 candidate genes, including DIR2, NAC-like, EXP1, CAD, β-GAL1, β-GAL2, 4CL, ENDOB, PE, BHLH, PG3, CEL, and PRX2, showed strong correlations and that of seven genes showed poor correlations with the corresponding RNA-seq data (Fig. 6; Table 3).
Moreover, 20 DAPs were selected for PRM analysis, of which 18 exhibited significantly different levels. Of these 18 DAPs, 14 (77.8%) showed the same trends in abundance between PRM and TMT quantification, including PE, PL, PG2, furostanol glycoside 26-O-beta-glucosidase (F26G), β-GAL2, auxin efflux carrier, alpha/beta hydrolase (α-HY), PRX2, PG4, PRX5, PRX3, endoglucanase 19, endoglucanase 8, and DIR1. Additionally, four genes (PRX, beta-fructofuranosidase, BXL, and BGLU33) showed inconsistent abundance compared with protein levels quantified by TMT (Table 4). In general, the trends in the expression changes measured by PRM and TMT were consistent. Notably, in the PM_PS and PL_PM groups, several genes involved in cell wall metabolism pathways showed consistent upregulation/downregulation in the transcriptome and proteome (Fig. 7).
Structural changes in the pericarp cell wall might affect A. trifoliata fruit cracking
Fruit cracking occurs when the stress exerted on the pericarp from the enlarged aril is greater than the strength of the fruit skin, and the mechanical strength of the pericarp depends largely on its cell wall . Studies have shown that jujube fruit cracking may be related to structural changes and rearrangement of the cell wall during the later stages of fruit ripening . Arrangement of the subcutaneous layers of cells was found to be relatively regular, and cell layers had a closer arrangement in the cracking-resistant tomato genotype . In this study, the cell wall of the pericarp had poor integrity, loose structures, deformed and reduced cell layers, and larger spaces and began to degrade in cracked fruits during ripening, consistent with previous results in grapes and oranges [32, 33]. This suggests that structural changes in the cell wall might play key roles in the occurrence of A. trifoliata fruit cracking.
General features of the transcriptomes and proteomes of different A. trifoliata pericarps
Fruit cracking is a key factor affecting seed yield. Elucidating the molecular mechanisms underlying fruit cracking could aid the utilization of A. trifoliata seeds for biofuels. However, the underlying mechanisms of A. trifoliata fruit cracking remain largely unknown. In this study, differences in the transcriptome and proteome were investigated based on RNA-seq and TMT data during different stages. As the transcriptome database was used for protein identification, the quality of the sequencing and assembly of the transcriptome data was crucial for subsequent analyses. In total, 186,054 unigenes (> 200 bp) were assembled in the transcriptome of A. trifoliata pericarp, and this number is much higher than that previously obtained for this species and other Ranunculales spp., e.g., A. trifoliata (11,749 by Yang et al. ; 65,757 by Niu et al. ), Dysosma aurantiocaulis (53,929) , and D. versipellis (44,855) . The obtained percentages of the Q30 bases (91.58–93.75%) and GC contents (45.94–48.48%) were similar to those reported in other studies of the A. trifoliata transcriptome, e.g., 89.06–93.33% and 43.20–43.93%  and 96.31% and 45.10% , respectively. In total, 100,329 (53.9%) and 56,346 (30.3%) of all unigenes identified in the current study matched those in the NR and SwissProt databases, respectively. These results are similar to values reported for D. aurantiocaulis (29,497 and 18,029; 54.70% and 33.43%)  and were higher than those presented in other similar A. trifoliata transcriptome studies, including 34,245 and 23,352 (32.14% and 21.91%) in the NR and SwissProt databases, respectively , and 19,096 (29.95%) in SwissProt . Thus, our results provided extensive sequence and unigene resources for A. trifoliata.
Moreover, 812,625 spectra, 10,572 unique peptides, and 2839 proteins were identified based on transcriptome data from A. trifoliata pericarp. The RNA-seq and protein sequencing methods identified and annotated many genes and proteins, providing a basis for a more precise and detailed descriptions of molecular processes and for the elucidation of complex physiological processes and their genetic regulation . The data presented here were sufficiently high and provided useful tools for future genetic research on fruit cracking in A. trifoliata and other Lardizabalaceae species. Additionally, there were many “uncharacterized,” “predicted,” or “putative” transcripts and proteins in the annotation results, indicating the limited nature of the analysis owing to a lack of genomic information . Thus, the roles of these unknown or uncharacterized genes and proteins in A. trifoliata fruit cracking should be investigated in future studies.
Translational and post-translational regulation of A. trifoliata fruit cracking
Some inevitable passive processes will occur in the cracked pericarp following fruit cracking, such as oxidative stress and microbial invasion . Therefore, differences in the expression of genes and proteins among PS, PM, and PL alone cannot accurately reflect the cause of fruit cracking. In this study, our comparative analysis of protein and gene expression levels indicated that more DEGs and shared DEGs than DAPs were identified in cracked fruits than in non-cracking fruits. This could be explained by the technical limitations of MS-based proteomics, such as low amount of readily available samples and low MS scanning rate for quantitative data acquisition, and the need for extensive fractionation, limiting the capacity to identify proteins . The PPI analysis indicated that the interactions of proteins were infrequent and weak, which was associated with the results that the interactions are often weak for many cellular processes, which are regulated by post-translational modifications that are recognized by specific domains in protein binding partners . Moreover, correlation analyses showed negative correlations between the proteome and transcriptome, indicating a discordance between the transcript levels and protein abundance. These results were similar to those observed in previous reports, suggesting that post-transcriptional and post-translational regulation, reversible phosphorylation, splicing events in cells, and translation efficiency play key roles in the regulation of fruit ripening [40, 41]. Thus, gene translation and post-translation processes could be important regulatory methods in fruit ripening and cracking.
Potential regulators and metabolic pathways during fruit cracking
GO and KEGG functional enrichment provide prediction information of inner-cell metabolic pathways and of the genetic and biologic behaviors of genes. GO terms associated with oxidoreductase activity and structural molecule activity were mainly enriched in the PM_PS and PL_PM groups, respectively, similar to previous studies showing that the target genes could be classified into different categories based on their functions, such as peroxidase and cell wall polysaccharide [20, 42]. Cell polysaccharides are degraded by cell wall hydrolases, and whereas phenolic crosslinking of cell wall structural components is catalyzed by cell wall peroxidase. These modifications reduce the strength of the fruit pericarp, resulting in changes in pericarp development and fruit cracking . KEGG analysis indicated that cell wall-related pathways, including pentose and glucuronate interconversions, the galactose metabolism pathway, and the phenylpropanoid biosynthesis pathway, were common pathways shared by DEGs and DAPs, suggesting that cell well metabolism may have important roles in A. trifoliata fruit ripening and cracking. Moreover, these results revealed that the proteomic data and transcriptome data were complementary, and that the proteome could confirm the transcriptome data; in addition, genes perform the same function at the transcriptome and proteome levels . Functional classification of the transcriptome and proteome could improve our understanding of the molecular physiology of fruit ripening and cracking.
Candidate DEGs and DAPs may play key roles in fruit pericarp cracking
Fruit cracking is a complex physiological process that is controlled by numerous genes working together, rather than a single gene directly controlling the process [18, 45]. Researchers have found that several cell wall modification genes, including β-GAL, β-GLU, PE and PG are differentially expressed in cracked fruits compared with levels in non-cracking litchi fruits . Downregulation of PpBGAL can delay peach fruit softening by reducing PG and pectin methylesterase activity . Antisense inhibition of PE and PG activity in tomatoes was found to reduce fruit cracking . Additionally, silencing the SIPL gene can enhance fruit firmness and reduce the content of pectin, suggesting that this gene participates in pericarp cell wall rearrangement during fruit softening . Fruit softening is mainly caused by hemicellulose and pectin degradation proteins, such as XYL, BGAL, PE, and PG [49, 50]. As a major cellulose degradation enzyme, downregulation of BGLU protein in strawberries can delay fruit maturation . PRXs are also involved in rearrangement of cell wall polysaccharides during plant development .
In this study, 13 DEGs and 14 DAPs involved in cell wall metabolism showed strong correlations with RNA-seq data and protein expression levels, consistent with results showing that cell polysaccharide metabolism might play key roles in fruit ripening and cracking [42, 43]. Notably, the expression levels of three proteins (β-GAL2, PE, and PL) associated with galactose metabolism and the pentose and glucuronate interconversion pathways showed higher protein and gene expression in cracking fruit than in non-cracking fruit, suggesting that PE, PL, and β-GAL2 may play important roles in the regulation of A. trifoliata fruit cracking.
In addition, F26G, PG, PG3, XYL, PRX3, and PRX5 also showed higher protein and gene expression in cracking fruit than in non-cracking fruit. The significantly increased expression of these cell wall metabolism proteins indicated a network of cell wall cellulose, hemicellulose, and pectin in fruit cracking during ripening [50, 53].
Our findings reveal the structural changes in the cell wall during different stages of A. trifoliata pericarp development, suggesting that structural changes between unripe and ripe fruits are important factors affecting the cracking tendency of fruits. To strengthen the structure of the pericarp, bagging and 1-MCP treatment, which reduce fruit cracking, could be performed [54, 55]. Comparative transcriptomic analysis and proteomic profiling were applied to identify the potential genes/proteins that participated in fruit cracking. Various genes and proteins were found to be differentially expressed after cracking. Three co-expressed genes (PE, PL, and β-GAL2) were implicated in various signaling pathways in fruit development, cell wall degradation, and fruit softening. These omics data provide a new perspective to understand the process of fruit cracking in A. trifoliate. In future works, transgenic plants overexpressing candidate genes are needed to validate the functions of these genes. Molecular markers for fruit cracking will be developed to utilize this trait in breeding and to increase seed yields for biorefineries.
The wild germplasm Nong No. 8, which had been transplanted 9 years prior in the nursery at Hunan Academy of Agricultural Sciences, Changsha, P. R. China, was used as the research material in this study. According to our observations, in Changsha, Hunan Province, the blooming stage for the germplasm of Nong No. 8 was in early April, when 50% of the A. trifoliata flowers were in bloom with a flowering period of approximately 30–45 days. Fruit development lasts for 5 months, from ovary inflation through fruit setting, and longitudinal and transverse elongation, flesh softening, and peel cracking complete the development process . The Nong No. 8 variety usually ripens in early October. Because the fruit of A. trifoliata do not develop uniformly, fruits were harvested from different stages of the same Nong No. 8 tree at the same time and were then sorted according to their developmental stage, considered to be a sample. Different fruits, including non-cracking fruits (PS), initially cracking fruits (PM), and total cracking fruits (PL) were randomly collected every 10 days during the ripening stage (September 18, September 28, and October 8, 2018) and three fruits were mixed into one biological replicate for further analysis. In total, three biological replicates were collected for each stage. Sampled pericarps were rapidly collected and immediately frozen in liquid nitrogen and stored at − 80 °C until use for transcriptome, proteome, qPCR, and PRM analyses.
Anatomical structure of the pericarp
The anatomy of pericarp samples collected from the Nong No. 8 fruits at different development stages, including PS, PM, and PL, were prepared for paraffin sectioning, and scanning electron microscopy (SEM) was carried out according to previous studies [57, 58]. Pericarp samples were fixed directly in the field using FAA (70% ethyl alcohol + 38% methyl aldehyde + 25% acetic acid [16:1:1]). The tissues were subsequently dehydrated through an ethanol series with increasing ethanol concentrations and embedded in paraffin. Subsequently, paraffin sections were stained with Safranin O and Fast Green staining and observed with an Axio Imager (Zeiss, Oberkochen, Germany); upright microscopy images were displayed using Image-pro Plus 6.0 software.
Pericarp samples were fixed in 2.5% glutaraldehyde (pH 7.4) for 4 h under syringe suction and washed subsequently in phosphate buffered-saline (PBS, 0.1 M, pH 7.0). After post-fixing with 1% osmium for 1–2 h and washing three times in PBS, samples were dehydrated with 30%, 50%, and 70% ethanol for 20 min and then with ethanol and iso-amyl acetate (V:V = 1:1) for 30 min. Samples were finally dried in liquid carbon dioxide. Dried samples were coated with gold–palladium in a Hitachi Model E-1010 ion sputter for 4–5 min and observed via SEM with a Hitachi Model SU-8010 (Hitachi, SU-8010, Japan).
RNA isolation, library construction, and sequencing
Total RNA used for the RNA-seq assays was isolated from three independent replicates of pericarp in the PS, PM, and PL stages, as described by Tao et al. . RNA samples were detected based on the A260/A280 absorbance ratio with a Nanodrop ND-1000 system (Thermo Scientific, USA). Paired-end libraries were prepared using a NEBNext UltraTM RNA Library Prep Kit for Illumina (NEB, USA) following the manufacturer’s instructions. The mRNA was purified from 3 µg total RNA using oligo (dT) magnetic beads followed by fragmentation carried out using divalent cations at elevated temperatures in NEBNext First Strand Synthesis Reaction Buffer. Subsequently, first-strand cDNAs were synthesized with random hexamer primers and Reverse Transcriptase (RNase H-) using mRNA fragments as templates, followed by second-strand cDNA synthesis using DNA polymerase I, RNAseH, buffer, and dNTPs. The synthesized double-stranded cDNA fragments were then purified with an AMPure XP system (Beckman Coulter, Beverly, MA, USA). The purified double-stranded cDNA was polyadenylated and adapter-ligated for preparation of the paired-end library. Adaptor-ligated cDNA and adaptor primers were used for PCR amplification. PCR products were purified (AMPure XP system) and library quality was assessed on an Agilent Bioanalyzer 2100 system. Finally, sequencing was performed with an Illumina HiSeq2500 instrument by Shanghai Applied Protein Technology (Shanghai, China).
Quality control and transcriptome assembly
The raw paired-end reads in fastq format produced from the sequencing were first processed with the in-house Perl scripts. Reads containing adapters, excess “N” nucleotides with more than 10% of the bases and reads of low-quality (reads with quality values ≤ 10) were removed by filter_fq software. The Q20, Q30, GC content, and sequence duplication levels of the obtained clean reads were calculated. The assembly of clean reads to unigene collections was performed using the Trinity software package (https://github.com/trinityrnaseq/trinityrnaseq/releases) . Trinity software consists of three independent software modules, including Inchworm, Chrysalis, and Butterfly, and transcripts less than 200 bp in length were discarded. Sequences containing the longest cluster transcripts without redundancy extracted from transcripts were considered unigenes.
The de novo assembled unigenes were annotated in five databases, including NR, Pfam, SwissProt, GO, and KEGG pathway databases, based on a BLAST search with an E-value threshold of 1E−5. Moreover, to further analyze the annotation results, GO and KEGG results with E-values of 1E−5 were used for functional gene annotation. GO terms could be classified into three categories, including BP, MF, and CC. In addition to the GO terms, the pathway maps were determined based on the KEGG database.
The normalized transcript abundances of the genes were estimated using the fragments per kilobase per million reads mapped based on the length of the gene and read counts mapped to this gene. DESeq2 R package (1.16.1) software was used to identify DEGs, and the false discovery rate (FDR) was controlled using the Benjamini and Hochberg’s approach to adjust the P value. Genes with an adjusted P value of less than 0.05 and absolute fold-change of 2 were deemed to be differentially expressed between the two samples. In addition, GO and KEGG pathway enrichment analysis of DEGs was implemented with the clusterProfiler R package. Transcription factor analysis of DEGs was performed against the PlantTFDB database (https://planttfdb.cbi.pku.edu.cn/). The heat map was visualized using heatmap 2.0 in the gplot R package.
Protein extraction from A. trifoliata pericarp was performed from each sample as described previously . The samples were frozen in liquid nitrogen and ground into powder. A five-fold volume of TCA/acetone (1:9) was added, and the sample was vortexed, mixed, and incubated at − 20℃ for 4 h. Samples were then then centrifuged at 6000×g at 4℃ for 40 min. The supernatant was discarded, and the precipitate was washed three times with pre-cooled acetone. After precipitation and air drying, the precipitate was redissolved in buffer (4% sodium dodecyl sulfate [SDS], 100 mM Tris–HCl and 1 mM DTT; pH 7.6). After sonication and boiling for 15 min, the lysate was centrifuged for 40 min, and the supernatant was filtered and quantified using the BCA Protein Assay Kit (Bio-Rad, USA).
Trypsin digestion and TMT labeling
For digestion, the samples were added to the buffer (4% SDS, 100 mM DTT, 150 mM tris–HCl, pH 8.0) and UA buffer (8 M urea, 150 mM Tris–HCl pH 8.0) by repeated ultrafiltration (Microcon units, 10 kDa). Then, iodoacetamide (100 mM IAA in UA buffer) was added to the samples to block reduced cysteine residues, and samples were incubated for 30 min in the dark at room temperature. After the filters were washed with UA buffer and triethylamine borane (TEAB) buffer in turn, the suspensions were digested with trypsin (Promega, Madison, WI, USA) in TEAB buffer overnight. After trypsin digestion, the samples (100 μg protein) were categorized for labeling with 129-tag (PS), 130-tag (PM), and 131-tag (PL; Thermo Fisher Scientific, Waltham, MA, USA). Finally, TMT-labeled peptide aliquots were pooled for subsequent fractionation using a Pierce high pH reversed-phase fractionation kit (Thermo scientific, USA).
High-performance liquid chromatography fractionation and liquid chromatography-tandem mass spectrometry (LC–MS/MS) analysis
For the fractionation of labeled peptides, samples were loaded onto a reverse-phase trap column (Thermo Scientific Acclaim PepMap100, 100 μm × 2 cm, nanoViper C18) connected to the C18 reverse-phase analytical column (Thermo Scientific Easy Column, 10 cm long, 75 μm inner diameter, 3 μm resin) in buffer A (0.1% formic acid) and separated with a linear gradient of buffer B (84% acetonitrile and 0.1% formic acid) at a flow rate of 300 nL/min controlled by IntelliFlow technology. The resulting peptides were further processed using a Q Exactive mass spectrometer (Thermo Scientific, USA) coupled with an Easy nLC (Thermo Fisher Scientific). MS analysis was performed in positive ion mode, and MS data were acquired using a data-dependent top10 method dynamically choosing the most abundant precursor ions from the survey scan (300–1800 m/z) using the higher-energy collision dissociation (HCD) fragmentation method. Automatic gain control (AGC) was set as 3E6 and the maximum inject time was 10 ms, with the following parameters: 40.0-s dynamic exclusion duration; 70,000 resolutions with survey scans at m/z 200 and resolution for HCD spectra at 17,500 at m/z 200; 2 m/z of isolation width; 30 eV of normalized collision energy and the underfill ratio was defined as 0.1%.
Sequence database search and data analysis
The raw data were processed by the MASCOT engine (Matrix Science, London, UK; version 2.2), and Proteome Discoverer 1.4 software was used to process MS/MS spectra. The search was performed using the following settings based on the A. trifoliata database: trypsin for the enzyme and 2 as the maximum missed cleavage allowed; fixed modifications of carbamidomethyl (C), TMT-6plex (N-term), and TMT-6plex (K), variable modification of oxidation (M); mass tolerance for fragment ions of 0.1 Da, and 20 ppm for peptide ions, as well as both peptide and protein levels of FDR less than 0.01, and only unique peptides of the protein were employed for the protein identification and quantification. Proteins with a P value less than 0.05 and fold-change greater than or equal to 1.2 or less than or equal to 0.83 within a comparison were recognized as DAPs.
Functional categorization was performed using GO and KEGG pathway databases with P values less than or equal to 0.05. The protein functional network was performed with STRING 9.0 software (https://string-db.org). Clustering analysis of the DEGs and DAPs was performed using Cluster 3.0 (https://bioservices.capitalbio.com/xzzq/rj/3885.shtml) and the Java Treeview software (https://jtreeview.sourceforge.net). Correlations were analyzed based on the DEGs and DAPs, and Person correlation tests were conducted for each comparison group, including PM versus PS and PL versus PM.
Reverse transcription qPCR
The method for total RNA extraction and synthesis of cDNA was described previously. A Bio-Rad CFX96 Touch detection system (Bio-Rad Laboratories, Richmond, CA, USA) with SYBR Green PCR master mix (Aidlab Biotechnologies, Co., Ltd.) was used for qPCR for each sample. In this study, the EF-1α gene was used as the internal control gene, which was detected by the de novo transcriptome sequencing of A. trifoliata . Primers for qPCR experiments were designed using Primer 5.0 software (Additional file 1: Table S5), and the gene sequences were blasted against the NCBI database. The amplification reactions contained 12.5 μL SYBR Green PCR master mix, 1 μL cDNA, and 0.5 μL of each primer in a final reaction volume of 25 μL. The thermal cycling program began with 3 min at 95 °C, followed by 40 cycles of 95 °C for 10 s and 55 °C for 30 s, with a melt curve from 65 to 95 °C based on increments of 0.5 °C for 5 s. After PCR amplification, quantitative variations were analyzed using the Delta Ct method, and the analysis of statistically significant differences in gene expression was performed by independent samples t-test analysis at P < 0.05 using GraphPad Prism 8 software. Correlation analysis between cell wall-related gene and protein expression was performed by Pearson’s correlation coefficient analysis .
To verify the protein expression levels obtained by TMT analysis, the expression levels of selected proteins were further quantified by PRM analysis. The AQUA stable isotope peptide was spiked in each sample as an internal standard reference. The tryptic peptides were directly loaded on C18 STAGE-tips for desalting prior to reverse-phase chromatography on an Easy nLC-1200 system (Thermo Scientific, USA). Acetonitrile was used for 45 min, increasing from 5 to 35% for a 1-h liquid chromatography gradient, and PRM analysis was performed on a Q Exactive Plus mass spectrometer (Thermo Scientific). The mass spectrometer was operated in positive ion mode, and the full MS1 scan was acquired with a resolution of 60,000 (at 200 m/z). AGC target values and maximum ion injection times were set at 3e6 and 200 ms, respectively; full MS scans were followed by 20 PRMs (MS2 scans) at 30,000 resolution (at m/z 200) with AGC of 3e6 and the maximum injection time set as 120 ms. The targeted peptides were isolated with a 2-Th window. Ion activation/dissociation was performed at a normalized collision energy of 27 in the HCD collision cell. Skyline version 3.7.0 was used to analyze the MS data where signal intensities for individual peptide sequences for each of the significantly altered proteins were quantified relative to each sample and normalized to a standard reference .
Availability of data and materials
The raw sequence data have been deposited in the National Center for Biotechnology Information (NCBI) database (https://www.ncbi.nlm.nih.gov/sra/SRP246017). The mass spectrometry proteomics data have been deposited to the PRIDE partner repository with the dataset identifier PXD017282.
- Akebia trifoliata :
Tandem mass tag
The non-cracking stage
The initial cracking stage
The total cracking stage
Basic Local Alignment Search Tool
Nonredundant protein sequence database
Protein families Fdatabase
Gene Ontology database
Kyoto Encyclopedia of Genes and Genomes pathway database
Differently expressed gene
Real-time quantitative polymerase chain reaction
Shikimate O-hydroxycinnamoyl transferase
Furostanol glycoside 26-O-beta-glucosidase
Differentially abundant protein
Niu J, Wang YJ, Shi YL, Wang XF, Sun ZM, Huang KY, Gong C, Luan MB, Chen JH. Development of SSR markers via de novo transcriptome assembly in Akebia trifoliata (Thunb.) Koidz. Genome. 2019;999:1–15.
Zhong CH, Bu FW, Wang ZY, Peng DF. Fruit development and biological characteristics of seedling progenies of Akebia trifoliata. Hunan Agr Sci. 2006;1:27–9.
Luo LP, Lan YH, Hu BF, Zhang X, Guo XL. Physicochemical properties and biodiesel preparation of Akebia trifoliate seed oil. Nat Prod Res Dev. 2013;25(8):1095–100.
Wang YJ, Ao WC, He XS, Duan WH, Xin WN, Cong C. Comparative study of physical and chemical components of Akebia trifoliata from 9 places of production in China. J West China Forestry Sci. 2016;45(6):45–8.
Yang BY, Zhang SH, Liu MY, Yang ZP, Li W. Exploitation and utilization of new-type liana oil plant Akebia trifoliata ssp. australis. Acta Agriculturae Jiangxi. 2014;26(9):20–2.
Liu C, Jiang SH, Xu LL, Li TJ, Han XJ, Wang JX, Liao L, Rong L. Biodiesel prepared by hydrocracking Akebia trifoliata oil. CCOA. 2015;30(1):76–80.
Wang R, Hanna MA, Zhou WW, Bhadury PS, Chen Q, Song BA. Production and selected fuel properties of biodiesel from promising non-edible oils: Euphorbia lathyris L. Sapium sebiferum L. and Jatropha curcas L. Bioresource Technol. 2011;102(2):1194–9.
Chung KH. Transesterification of Camellia japonica and Vernicia fordii seed oils on alkali catalysts for biodiesel production. J Ind Eng Chem. 2010;16(4):506–9.
Wang DZ, Li F, Yan J, Zhong HM. Study and application of nutritional components of wild plant Var australis (Diels) Rehd. Amino Acids Biotic Resour. 2004;26:16–7.
Ishmayana S, Learmonth RP, Kennedy UJ. Fermentation performance of the yeast Saccharomyces cerevisiae in media with high sugar concentration. In: The 2nd international seminar on chemistry; 2011, p. 379–85.
Yu Q, Zhang A, Wang W, Chen L, Ruxue B, Zhuang XS, Wang Q, Wang ZM, Yuan ZH. Deep eutectic solvents from hemicellulose-derived acids for the cellulosic ethanol refining of Akebia’ herbal residues. Bioresource Technol. 2017;247:705–10.
Li L, Yao X, Zhong C, Chen X, Huang H. Akebia: a potential new fruit crop in China. Hortic Sci. 2010;45:4–10.
Zou S, Yao X, Zhong C, Zhao T, Huang H. Effectiveness of recurrent selection in Akebia trifoliata (Lardizabalaceae) breeding. Sci Hortic. 2019;246:79–85.
Joshi M, Baghel RS, Fogelman E, Stern RA, Ginzberg I. Identification of candidate genes mediating apple fruit-cracking resistance following the application of gibberellic acids 4+7 and the cytokinin 6-benzyladenine. Plant Physiol Bioch. 2018;127:436–45.
Huang X, Wang HC, Li J, Yin J, Yuan W, Lu J, Huang HB. An overview of calcium’s role in lychee fruit cracking. Acta Hortic. 2003;665:1026–35.
Liu QX, Wen LZ, Zhou HJ, Wu Q, Marvin KH, Xiao ZS. Effects of seed wasps and seedless fruits on fruit and oil yields of Pistacia chinensis as a biofuel tree. Acta Ent Sinica. 2011;54(10):1133–9.
Wu J, Xu Z, Zhang Y, Chai L, Yi H, Deng X. An integrative analysis of the transcriptome and proteome of the pulp of a spontaneous late-ripening sweet orange mutant and its wild type improves our understanding of fruit ripening in citrus. J Exp Bot. 2014;65(6):1651–71.
Mustafa M, Syukur M, Sutjahjo SH. Inheritance of radial fruit cracking resistance in tomatoes (Solanum lycopersicum L.). IOP Conf Ser Earth Environ Sci. 2019;270(1):012032.
Correia S, Schouten R, Silva AP, Goncalves B. Sweet cherry fruit cracking mechanisms and prevention strategies: a review. Sci Hortic. 2018;240:369–77.
Jiang F, Lopez A, Jeon S, de Freitas ST, Yu Q, Wu Z, Mitcham E. Disassembly of the fruit cell wall by the ripening-associated polygalacturonase and expansin influences tomato cracking. Hortic Res. 2019;6(1):1–15.
Li WC, Wu JY, Zhang HN, Shi SY, Liu LQ, Shu B, Wei YZ. De novo assembly and characterization of pericarp transcriptome and identification of candidate genes mediating fruit cracking in Litchi chinensis Sonn. Int J Mol Sci. 2014;15(10):17667–85.
Khurnpoon L, Siriphanich J, Labavitch JM. Cell wall metabolism during durian fruit dehiscence. Postharvest Biol Tec. 2008;48(3):391–401.
Lenser T, Theißen G. Conservation of fruit dehiscence pathways between Lepidium campestre and Arabidopsis thaliana sheds light on the regulation of indehiscent. Plant J. 2013;76(4):545–56.
Dong Y, Yang X, Liu J, Wang BH, Liu BL, Wang YZ. Pod shattering resistance associated with domestication is mediated by a NAC gene in soybean. Nat Commun. 2014;5:3352.
Sorefan K, Girin T, Liljegren SJ, Ljung K, Robles P, Galván-Ampudia CS, Østergaard L. A regulated auxin minimum is required for seed dispersal in Arabidopsis. Nat. 2009;459(7246):583.
Yang H, Liu H, Shi X, Ge F, Zhao Z, Luo H. Transcriptome analysis of Akebia trifoliate in ripening process. J Mountain Agrie Biol. 2016;35(2):46–51.
Wang J, Gao X, Ma Z, Chen J, Liu Y. Analysis of the molecular basis of fruit cracking susceptibility in Litchi chinensis cv. Baitangying by transcriptome and quantitative proteome profiling. J Plant Physiol. 2019;234:106–16.
Voelckel C, Gruenheit N, Lockhart P. Evolutionary transcriptomics and proteomics: insight into plant adaptation. Trends Plant Sci. 2017;22(6):462–71.
Huang XM, Yuang WQ, Wang HC, Li JG, Huang HB, Luo S, Yin JH. Linking cracking resistance and fruit desiccation rate to pericarp structure in litchi (Litchi Chinensis Sonn.). J Hortic Sci Biotech. 2004;79:897–905.
Wang BM, Ding GX, Wang XY, Fu CB, Qin GJ, Yang JQ, Wen PF. Changes of histological structure and water potential of huping jujube fruit cracking. Scientia Agr Sinica. 2013;46(21):4558–68.
Yang Z, Wu Z, Zhang C, Hu E, Zhou R, Jiang F. The composition of pericarp, cell aging, and changes in water absorption in two tomato genotypes: mechanism, factors, and potential role in fruit cracking. Acta Physiol Plant. 2016;38(9):215.
Li J, Chen JZ, Wan JF, Huang YJ, Yao Q. The correlation of cracking fruits types and peel structure between ‘Washington’ and ‘Bonanza’ navel oranges. Chin J Trop Crop. 2011;32:921–5.
Zhang C, Guan L, Fan X, Zheng T, Fang J. Anatomical characteristics associated with different degrees of berry cracking in grapes. Sci Hortic. 2020;261:108992.
Mao Y, Zhang Y, Xu C, Qiu Y. Comparative transcriptome resources of two Dysosma species (Berberidaceae) and molecular evolution of the CYP 719A gene in Podophylloideae. Mol Ecol Resour. 2016;16(1):228–41.
Guo R, Mao YR, Cai JR, Wang JY, Wu J, Qiu YX. Characterization and cross-species transferability of EST–SSR markers developed from the transcriptome of Dysosma versipellis (Berberidaceae) and their application to population genetic studies. Mol Breeding. 2014;34(4):1733–46.
Blencowe BJ, Ahmad S, Lee LJ. Current-generation high-throughput sequencing: deepening insights into mammalian transcriptomes. Genes Dev. 2009;23:1379–86.
Cramaro WJ, Revets D, Hunewald OE, Sinner R, Reye A, Muller CP. Integration of Ixodes ricinus genome sequencing with transcriptome and proteome annotation of the nave midgut. BMC Genomics. 2015;16(1):871.
Zhou Y, Meng Z, Edman-Woolcott M, Hamm-Alvarez SF, Zandi E. Multidimensional separation using HILIC and SCX pre-fractionation for RP LC-MS/MS platform with automated exclusion list-based MS data acquisition with increased protein quantification. J Proteomics Bioinform. 2015;8(11):260–5.
Tripsianes K, Chu NK, Friberg A, Sattler M, Becker CF. Studying weak and dynamic interactions of posttranslationally modified proteins using expressed protein ligation. ACS Chem Biol. 2014;9(2):347–52.
Luo X, Cao D, Li H, Zhao D, Xue H, Niu J, Cao S. Complementary iTRAQ-based proteomic and RNA sequencing-based transcriptomic analyses reveal a complex network regulating pomegranate (Punica granatum L.) fruit peel colour. Sci Rep. 2018;8(1):12362.
Chen X, Tao Y, Ali A, Zhuang Z, Guo D, Guo Q, Wang J. Transcriptome and proteome profiling of different colored rice reveals physiological dynamics involved in the flavonoid pathway. Int J Mol Sci. 2019;20(10):2463.
Xue LZ, Sun MT, Wu Z, Yu L, Yu QH, Tang YP, Jiang FL. LncRNA regulates tomato fruit cracking by coordinating gene expression via a hormone-redox-cell wall network. BMC Plant Biol. 2020;20:162.
Li J, Chen JZ. Citrus fruit-cracking: causes and occurrence. HPJ. 2017;3(6):255–60.
Muers M. Gene expression: transcriptome to proteome and back to genome. Nat Rev Genet. 2011;12(8):518.
Chen J, Duan Y, Hu Y, Li W, Xie J. Transcriptome analysis of atemoya pericarp elucidates the role of polysaccharide metabolism in fruit ripening and cracking after harvest. BMC Plant Biol. 2019;19(1):2–19.
Liu H, Qian M, Song CH, Li JJ, Zhao CP, Li GF, Wang AZ, Han MY. Down-regulation of PpGAL10 and PpGAL16 delays fruit softening in peach by reducing polygalacturonase and pectin methylesterase activity. Front Plant Sci. 2018;9:1015.
Brummell DA, Harpster MH. Cell wall metabolism in fruit softening and quality and its manipulation in transgenic plants. Plant Mol Biol. 2001;47:311–8.
Yang L, Huang W, Xiong F, Xian Z, Su D, Ren M, et al. Silencing of SlPL, which encodes a pectate lyase in tomato, confers enhanced fruit firmness, prolonged shelf-life, and reduced susceptibility to gray mold. Plant Biotechnol J. 2017;15(9):1544–55.
Bustamante CA, Civello PM, Martínez GA. Cloning of the promoter region of β-xylosidase (FaXyl1) gene and effect of plant growth regulators on the expression of FaXyl1 in strawberry fruit. Plant Sci. 2009;177:49–56.
Song XX, Hu Z, Bi Y, Jiang AL. Research progress in softening mechanism and regulation of fresh-cut fruits and vegetables. Sci Technol Food Ind. 2013;34(4):397–9.
Sara F, Simona M, Chiara M. Fruit ripening: the role of hormones, cell wall modifications, and their relationship with pathogens. J Exp Bot. 2019;70(11):2993–3006.
Francoz E, Ranocha P, Nguyen-Kim H, Jamet E, Burlat V, Dunand C. Roles of cell wall peroxidases in plant development. Phytochemistry. 2015;112:15–211.
Huang XM, Wang HC, Lu XJ, Yuan WQ, Lu JM, Li JG, Huang HB. Cell wall modifications in the pericarp of litchi (Litchi chinensis sonn.) cultivars that differ in their resistance to cracking. J Hortic Sci Biotechnol. 2006;81(2):231–7.
Son IC, Kim DI. Effects of bagging periods on pericarp characteristics and berry cracking in 'Kyoho' Grape (Vitis sp.). Korean J Hortic Sci Technol. 2010;28(3):381–6.
Jeong J, Huber DJ. Suppression of avocado (Persea americana Mill.) fruit softening and changes in cell wall matrix polysaccharides and enzyme activities: differential responses to 1-MCP and delayed ethylene application. J Am Soc Hortic Sci. 2004;129:752–9.
Zhang XF. Study on the physiological and ecological characteristics and quality of wild Akebia trifoliate (Thunb.) Koidz. Guiyang: Guizhou University; 2017.
Zhang EZ, Liao F, Yang RG, He QG, Huang MH, Xin M, Huang ZY. Comparison of cell structure of the fruit peel in seventeen varieties of Mangifera indica L during different postharvest stages. Chin J Tropical Crops. 2016;37(11):2238–43.
Chen L, Zhang J, Li H, Niu J, Xue H, Liu B, Cao S. Transcriptomic analysis reveals candidate genes for female sterility in pomegranate flowers. Front Plant Sci. 2017;8:1430.
Tao X, Fang Y, Xiao Y, Jin YL, Ma XR, Zhao Y, Wang HY. Comparative transcriptome analysis to investigate the high starch accumulation of duckweed (Landoltia punctata) under nutrient starvation. Biotechnol Biofuels. 2013;6(1):72.
Grabherr MG, Haas BJ, Yassour M, Levin JZ, Thompson DA, Amit I, Chen Z. Full-length transcriptome assembly from RNA-Seq data without a reference genome. Nat Biotechnol. 2011;29(7):644.
Liu Y, Cao D, Ma L, Jin X, Yang P, Ye F, Wei C. TMT-based quantitative proteomics analysis reveals the response of tea plant (Camellia sinensis) to fluoride. J Proteomics. 2018;176:71–81.
Yang H. Analysis of gene expression on ethylene biosynthesis of Akebia trifoliate (Thunb.) Koidz in ripping process. Guiyang: Guizhou University; 2016.
Sato A, Zhang T, Yonekura L, Tamura H. Antiallergic activities of eleven onions (Allium cepa) were attributed to quercetin 4'-glucoside using QuEChERS method and Pearson's correlation coefficient. J Funct Foods. 2015;14:581–9.
MacLean B, Tomazela DM, Shulman N, Chambers M, Finney GL, Frewen B, Kern R, Tabb DL, Liebler DC, MacCoss MJ. Skyline: an open source document editor for creating and analyzing targeted proteomics experiments. Bioinformatics. 2010;26(7):966–8.
This research was financially supported by the Central Public- interest Scientific Institution Basal Research Fund (1610242020010) and Agricultural Science, Technology Innovation Program (ASTIP) of CAAS (Grant No.2017IBFC).
Ethics approval and consent to participate
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Additional file 1: Fig. S1.
The size distribution of the assembled transcript and unigene sequences of A. trifoliate. Fig. S2. Volcano plot depicting the transcriptome and proteome data of A. trifoliate. a, c Volcano plot depicting the transcriptome and proteome data in PM_PS. b, d Volcano plot depicting the transcriptome and proteome data in PL_PM. Absolute log10 and log2 fold changes are plotted on the y-axis and x-axis, respectively. Horizontal dotted line presents p values of 0.05 cut-off position while the vertical dotted lines discriminate between genes and proteins having absolute log2 fold change of 1. Red dots represent significantly different upregulated genes. Green dots represent significantly different downregulated genes. Pink dots represent a log2 fold change >1 with p < 0.05 in protein expression. Black dots indicate no difference in gene and protein expression. Fig. S3. Molecular weight and peptide count distribution were identified from TMT proteomics by searching against the database. a Distribution of the proteins that were identified among different molecular weights. b Distribution of peptide count of the proteins were identified from TMT data. Fig. S4. GO and KEGG pathway functional enrichment analysis of co-regulated genes and proteins in A. trifoliate. a, b GO enrichment analysis of co-regulated genes and proteins in PM_PS and PL_PM, respectively. c, d KEGG pathways enrichment analysis of co-regulated genes and proteins in PM_PS and PL_PM, respectively. Table S1. Sequencing statistics for A. trifoliate. Table S2. Statistical data of unigenes annotation. Table S3. The summary of the total number of transcripts and proteins identified from different stages and replicates. Table S4. The summary of the DEGs and DAPs identified from different stages and replicates. Table S5. Sequences of specific primers used for qPCR experiment.
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.
About this article
Cite this article
Niu, J., Shi, Y., Huang, K. et al. Integrative transcriptome and proteome analyses provide new insights into different stages of Akebia trifoliata fruit cracking during ripening. Biotechnol Biofuels 13, 149 (2020). https://doi.org/10.1186/s13068-020-01789-7
- Akebia trifoliata
- Pericarp structure
- Fruit cracking
- Pentose and glucuronate interconversions
- Galactose metabolism