Figures (5)  Tables (2)
    • Figure 1. 

      cDNA-AFLP profiles using one EcoR I selective primer and eight Mes I selective primers. One EcoR I selective primer: EcoR I-AC; Eight Mes I selective primers: Mes I-AA, AG, AC, AT, CC, CG, CT, and CA; Lane 1: Normal juice sacs of C. grandis; Lane 2: Granulated juice sacs of C. grandis; Arrows indicate differentially expressed transcript-derived fragments.

    • Figure 2. 

      Functional classification of differentially expressed transcript-derived fragments (TDFs) in normal and granulated juice sacs of C. grandis. Functional classification was performed based on information reported for each sequence in the NCBI database (https://blast.ncbi.nlm.nih.gov/Blast.cgi).

    • Figure 3. 

      Relative expression levels of transcript-derived fragments (TDFs) in C. grandis normal and granulated juice sacs. (a) Relative expression levels of genes encoding β-amylase 4 (TDF #64-2); cytokinin-O-glucosyltransferase 1 (TDF #94-1); gibberellin 20 oxidase (TDF #13-2); galactose-1-phosphate guanylyltransferases (TDF #200-1); α-galactosidase precursor (TDF #160-1); cytochrome b5 (TDF #123-2); ADP-ribosylation factor 3 (TDF #246-4); 2-oxoglutarate-dependent dioxygenase (TDF #13-1); 1,4-alpha-glucan-maltohydrolase (TDF #20-3); ethylene insensitive 3-like protein (TDF #147-4). (b) Relative expression levels of genes encoding translation initiation factor 4A2 (TDF #14-3); cytochrome P450 (TDF #48-3); auxin down-regulated-like protein (TDF #5-3); cellulose synthase (TDF #249-3); transport protein SEC31 (TDF #34-2); Ca2+-transporting ATPase (TDF #252-3); heat shock protein (TDF #252-1); senescence-associated protein (TDF #5-1); dicer-like protein 4 (TDF #244-6); cell wall-associated hydrolase (TDF #216-3). Results are shown as the mean ± SD of three independent experiments. Different letters above bars indicate significant differences at P < 0.05.

    • Figure 4. 

      Granulation led to alterations in hormone content in C. grandis juice sacs. The content of five hormones were measured by ultra-performance liquid chromatography mass-spectrometry (UPLC-MS). Normal juice sacs were considered as negative control. GA3: gibberellin 3; GA20: gibberellin 20, CK: cytokinins; ZT: zeatin; and IAA: indole-3-acetic acid. Error bars represent standard deviations calculated from three biological replicates. Different letters above bars indicate significant differences at P < 0.05.

    • Figure 5. 

      Proposed regulatory network for the granulation phenomenon in C. grandis juice sacs. Red arrows indicate upregulated genes. HSP: heat-shock protein.

    • Found only in normalFound only in granulatedFound in both juice sacsTotal
      Total TDFs detected5368743,0144,424
      Total differentially expressed TDFs detected688926183
      TDFs produced useable sequence data416114116
      TDFs encoding known or putative proteins385711106
      TDFs encoding predicted, uncharacterized89219
      TDFs without matches in the database1126441

      Table 1. 

      Summary of transcript-derived fragments (TDFs) in normal and granulated juice sacs of C. grandis.

    • TDF numberSize
      (bp)
      Homologous proteinOrganism of originE-valueDegree of similarity (%)Genbank IDFold change
      Hormone and biological regulation
      48-2292CAAX amino terminal protease familyCucumis melo7E-1890ADN33781.1+
      5-3259Auxin down-regulated-like protein, partialPicea sitchensis5E-0555ADM77850.1+
      6-2158Protein embryo defective 2752Arabidopsis thaliana2E-1261NP_567830.10
      118-3290Expressed proteinOryza sativa5E-2763ABF95726.10
      94-1358Cytokinin-O-glucosyltransferase 1Aegilops tauschii9E-2548EMT28784.1+
      103-1211Nitrate reductaseCucumis sativus1E-3089ADK77877+
      219-5244ARF-GAP domain 2Arabidopsis lyrata subsp. lyrata5E-0943NP_176283.10
      13-2269Gibberellin 20 oxidaseMedicago truncatula4E-4281AES62614.2+
      147-4271EIN3-like (Ethylene insensitive 3) proteinC. melo8E-41100BAB64345+
      195-6197MATE efflux family proteinTheobroma cacao2E-1979EOX90702.10
      115-1355Transmembrane emp24 domain-containing protein p24delta9-likeCrocus sativus3E-6098XP_004143772.1+
      195-1455Multicatalytic endopeptidase complexA. thaliana8E-6984CAA74030.1+
      Carbohydrate and energy metabolism
      134-1192ATP-binding protein-likeA. thaliana2E-1075BAB09414.1+
      249-1311UDP-glucosyltransferase family 1 proteinCamellia sinensis2E-2593ACS87991.1+
      6-1214Fructose-bisphosphate aldolaseLemna minor1E-3394ACD10928.10
      5-4148ATP synthase subunit betaMedicago truncatula2E-1395XP_003627732.1+
      118-2304Glyceraldehyde-3-phosphate dehydrogenase, partialVernicia fordii4E-1194AFJ04516.10
      118-1441Glycosyltransferase, CAZy family GT8Populus trichocarpa5E-7590XP_002312381.1+
      64-2220β-amylase 4Citrus trifoliata2E-2269AFQ33616+
      249-3217Cellulose synthasePopulus tremula × Populus tremuloides5E-1291AAT09895.1+
      141-3499Mitochondrial benzaldehyde dehydrogenaseAntirrhinum majus8E-8980ACM89738.1+
      160-1174α-galactosidase precursorCoffea arabica2E-0462CAJ40777.1+
      13-12692-oxoglutarate-dependent dioxygenasePopulus trichocarpa8E-4889XP_002330269.17.07 ± 0.52
      200-1389Galactose-1-phosphate guanylyltransferasesT. cacao1E-6583EOY12255.1+
      20-32201,4-alpha-glucan-maltohydrolaseSolanum lycopersicum1E-1460NP_001234052.1+
      246-3192ABC transporter family proteinP. trichocarpa1E-1280XP_002310031.20
      221-3123Diacylglycerol kinase-like proteinA. thaliana6E-1169BAB09587.1+
      246-4172ADP-ribosylation factor 3A. lyrata subsp. lyrata4E-9896XP_002869315.10
      65-1363Methylenetetrahydrofolate reductase family protein isoform 3T. cacao7E-4972EOY04345.10
      Protein and nucleic acid metabolism
      228-1243Ribonucleoside-diphosphate reductase
      subunit M1
      A. thaliana2E-3389AEC07222.10
      251-1188BET1P/SFT1P-like protein 14AA. thaliana6E-0778NP_191376.10
      253-2306Class II aaRS and biotin synthetases
      superfamily protein
      A. thaliana5E-4183NP_186925.4+
      130-3193Ribosomal protein L5Citrullus lanatus1E-3096YP_003587255.10
      119-3254Nuclear transport factor 2 family proteinT. cacao1E-2262EOY06196.1+
      27-3224Glutathione S-transferase family proteinT. cacao1E-0762EOY27562.1+
      195-5268BRCA1-associated proteinM. truncatula8E-3973XP_003609376.20
      54-2169Pre-mRNA splicing factor PRP38 family proteinP.trichocarpa2E-2292ERP53525.10
      151-4347Ribonuclease II family proteinA. thaliana2E-5180NP_565418.1+
      127-1180Mitochondrial substrate carrier family protein isoform 2T. cacao2E-2861EOY07093.1+
      15-2349Nuclear transcription factor Y subunit B18M. truncatula1E-0968AFK49668.1+
      15-1392RRNA intron-encoded homing endonucleaseM. truncatula6E-1388XP_003614385.1+
      20-1280Solute carrier family 25 memberM. truncatula6E-2363XP_003615848.20
      221-1256TPA: heterogeneous nuclear ribonucleoprotein A3-like protein 2 isoform 1Zea mays1E-2982DAA58966.10.56 ± 0.06
      217-5194Adenine nucleotide alpha hydrolases-like superfamily proteinT. cacao2E-0696EOY06709.11.86 ± 0.34
      215-1363Transcription regulatorA. thaliana1E-0688NP_171710.4+
      195-4333PLP-dependent transferases superfamily proteinA. thaliana2E-3760NP_191772.1+
      21-221060S ribosomal protein L24-1T. cacao1E-3797EOY23121.1+
      155-3187Valyl-tRNA synthetase/valine-tRNA ligaseT. cacao3E-19
      75
      EOY31957.1
      2.23 ± 0.11
      252-2299ACT domain-containing protein ACR8Ricinus communis2E-3086XP_002509632.10.79 ± 0.03
      151-5
      347Exosome complex exonuclease RRP44
      homolog A
      R. communis5E-5382XP_002521738.1+
      217-2369Transcription initiation factor TFIID subunit AA. thaliana1E-2277NP_564023.1+
      218-1239Chaperonin 60 alpha subunitArachis diogoi1E-2086ACA23472.10
      31-1293Aspartyl protease family proteinA. thaliana3E-2357XP_002891474.1+
      236-3251Protein kinase domain-containing proteinA. thaliana3E-3080AEE27605.10
      247-1324Importin beta-2 subunit family proteinA. thaliana2E-1676XP_002867489.1+
      253-1263Serine/threonine protein kinase TNNI3KM. truncatula2E-0976XP_003601186.10
      5-2151Amino acid adenylation proteinCalothrix sp. PCC 63039.940YP_007137552.1+
      119-1290Ricin B-like lectin EULS3A. lyrata subsp. lyrata4E-3673XP_002862306.2+
      119-2190Spl1-Related 2 proteinA. thaliana2E-1876CAB56773.1+
      128-1286Chloroplast elongation factor TuB (EF-TuB)Nicotiana sylvestris7E-1290BAA01975.1+
      94-2242Clone 6F8 eukaryotic initiation factor 4A-14 geneNicotiana benthamiana2E-4789JN688263.1+
      34-2343Transport protein SEC31T. cacao2E-0681EOY23302.1+
      91-1199Phosphatase 2C family protein isoform 2T. cacao2E-8387EOY06499.1+
      160-6271Kinase superfamily protein isoform 1T. cacao3E-3478EOY06443.10
      209-3337Ser/Thr phosphatase-containing Kelch repeat domain protein, partialN. benthamiana1E-5090AFN44702.10
      220-3208Pentatricopeptide repeat (PPR) superfamily protein isoform 2T. cacao5E -0851EOY04957.12.25 ± 0.22
      197-1398Ubiquitin-specific protease family C19-related proteinA. thaliana9E-6879NP_564009.1+
      87-2220Tetratricopeptide repeat (TPR)-like superfamily protein isoform 1T. cacao3E-2786EOY33236.12.55 ± 0.06
      209-5319Serine/threonine-protein phosphatase, partialGenlisea aurea3E-5288EPS64063.10
      194-1474Proteasome subunit beta type-4Solanum nigrum6E-7391ADW66147.10
      150-9199Serine/threonine-protein kinase AtPK2/AtPK19R. communis2.851XP_002528702.1+
      123-4246IFA binding proteinLilium longiflorum2E-3774ABM68547.10
      123-3302Dual specificity kinase 1Arabidopsis thaliana3E-2793AEE276350
      26-2324Dephospho-CoA kinaseA. thaliana1E-4369NP_180318.10
      14-3222Translation initiation factor 4A2Z. mays7E-3797AAD20980.10.65 ± 0.01
      244-1302Translation elongation factor, partialAmmopiptanthus mongolicus1E-0588AFC01200.1+
      246-2165Eukaryotic translation initiation factor 5
      isoform 2
      T. cacao3E-0552EOX90767.10
      Lipid metabolism
      123-1278Patellin-5A. lyrata subsp. lyrata3E-1367XP_002872438.10
      141-4209Patellin-5Triticum urartu8E-2771EMS62735
      XP_003623
      0.27 ± 0.02
      195-3347Non-specific lipid-transfer proteinM. truncatula1E-0592596.3+
      197-3327Pleckstrin (PH) and lipid-binding START domains-containing protein isoform 2T. cacao1E-4778EOY34386.1+
      125-2278Glyoxylate/hydroxypyruvate reductase
      A HPR2
      A. lyrata subsp. lyrata1E-2970XP_002889322.1+
      Stress response and defense
      48-1263Trehalose 6-phosphate synthaseNicotiana tabacum3E-1891BAI99252.1+
      220-2208Transcription factor bHLH130M. truncatula4E-2275XP_003590427.11.41 ± 0.06
      123-2189Cytochrome b5N. tabacum9E-2680CAA50575 +
      48-3268Cytochrome P450Citrus sinensis2E-2895AAL24049.1+
      64-1297Cytochrome P450A. thaliana1E-2381NP_176086.1+
      221-2451Cytochrome oxidase subunit 1Curcuma longa1E-1056ABY83898.1+
      218-2341DNA damage-binding protein, partialM. truncatula3E-6467XP003638090.10
      252-1353Heat shock proteinM. truncatula6E-0745XP_003621962.10.36 ± 0.04
      115-2290Stress responsive proteinZ. mays5E-3264NP_001149550.1+
      5-1230Senescence-associated protein
      Picea abies2E-4596ACA04850.1+
      154-5184Dehydration-induced 19-like proteinGossypium hirsutum8E-0556ADP30960.1+
      27-1271B-box zinc finger proteinBambusa oldhamii7E-0957ACF35275.10
      244-6248Dicer-like protein 4A. thaliana9E-0963NP_197532.3+
      Cell metabolism
      252-3355Ca2+-transporting ATPaseA. thaliana2E-0472NP_195479+
      130-1240Plasma membrane isoform 4T. cacao3E-2983EOY10146.10
      244-2287Cinnamyl-alcohol dehydrogenaseA. thaliana3E-0591AAA99511.1+
      40-1163Clathrin adaptor complexes medium subunit family proteinA. lyrata subsp. lyrata3E-2293XP_002886592.10
      154-3429RAB GTPase homolog A5AA. thaliana1E-1790NP_199563.10.67 ± 0.01
      217-6323Receptor-like kinase binding proteinP. trichocarpa4E-2858XP_002325092.10.87 ± 0.05
      147-3168Calreticulin-like protein
      Solanum melongena6E-2288BAA85118.13.35 ± 0.12
      16-3238Ycf2 (chloroplast)Gossypium raimondii3E-3796YP_005087735.1+
      143-4350Tetraspanin8T. cacao9E-2469EOY31574.1+
      216-3443Cell wall-associated hydrolaseVibrio cholerae5E-3976ACX81677.10
      200-2283Nucleic acid binding proteinZ. mays5E-4389NP_001152488.1+
      Unknown biological processes
      217-3323Hypothetical protein AT5G07270A. thaliana4E-0959NP_1963440
      236-1256Choline/ethanolamine kinase, putativeRicinus communis4E-3583XP_002525542.1+
      236-2322Uncharacterized protein LOC8268581R. communis3E-3182XP_002530954.10
      249-2196Amino acid transporter, putativeR. communis3E-0339XP_002531860.1+
      8-2168Predicted: monoacylglycerol lipase abhd6-B-likeFragaria vesca subsp. vesca1E-1880XP_004303453.10
      119-4310Domain of uncharacterized protein function 724 6, putative isoform 1T. cacao2E-2268EOX95351.10
      57-1205Predicted: Vitis vinifera peroxidase 3-like mRNAVitis vinifera2E-2191XM_002280238.4+
      100-1318Hypothetical protein CICLE_v10006049mgCitrus clementina6E-14100ESR32793.10
      17-3206Hypothetical protein MTR_2g077840M. truncatula2E-376XP_003596462.10
      197-2335Mitochondrial protein, putativeM. truncatula3E-0488XP_003588355.10
      217-4261Putative ATP synthetase alpha chainOryza sativa subsp. japonica3E-1363AAO72570.10
      154-1202Hypothetical protein CICLE_v10022616mgCitrus clementina7E-2295ESR54213.10
      89-3215Hypothetical protein CICLE_v10033239mgC. clementina3E-19100ESR51519.1+
      TDFs: Transcript-derived fragments. Results are shown as the mean ± SD of at least three independent experiments. Fold change: 0 indicates TDFs only detected in normal juice sacs; + indicates TDFs only detected in granulated juice sacs. Relative expression ratio was obtained by analyzing gel images using PDQuest version 8.0.1 (Bio-Rad, Hercules, CA, USA).

      Table 2. 

      Homology of differentially expressed cDNA-AFLP fragments with known gene sequences in the database using BLASTX algorithm along with their expression patterns in granulation juice sacs of C. grandis.