Next Article in Journal
Chromosome Evolution in Marsupials
Previous Article in Journal
Genome-Wide Identification and Expression Analysis of the HD-Zip Gene Family in Wheat (Triticum aestivum L.)
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Transcriptomic Analysis Provides Insights into Grafting Union Development in Pecan (Carya illinoinensis)

1
College of Forestry, Nanjing Forestry University, Nanjing 210037, China
2
Co-Innovation Center for the Sustainable Forestry in Southern China, Nanjing Forestry University, Nanjing 210037, China
*
Author to whom correspondence should be addressed.
Submission received: 30 November 2017 / Revised: 18 January 2018 / Accepted: 26 January 2018 / Published: 5 February 2018
(This article belongs to the Section Plant Genetics and Genomics)

Abstract

:
Pecan (Carya illinoinensis), as a popular nut tree, has been widely planted in China in recent years. Grafting is an important technique for its cultivation. For a successful grafting, graft union development generally involves the formation of callus and vascular bundles at the graft union. To explore the molecular mechanism of graft union development, we applied high throughput RNA sequencing to investigate the transcriptomic profiles of graft union at four timepoints (0 days, 8 days, 15 days, and 30 days) during the pecan grafting process. After de novo assembly, 83,693 unigenes were obtained, and 40,069 of them were annotated. A total of 12,180 differentially expressed genes were identified between by grafting. Genes involved in hormone signaling, cell proliferation, xylem differentiation, cell elongation, secondary cell wall deposition, programmed cell death, and reactive oxygen species (ROS) scavenging showed significant differential expression during the graft union developmental process. In addition, we found that the content of auxin, cytokinin, and gibberellin were accumulated at the graft unions during the grafting process. These results will aid in our understanding of successful grafting in the future.

1. Introduction

Pecan (Carya illinoinensis), a member of the Juglandaceae family, is an economically important nut tree native to North America. It was introduced to China more than 100 years ago; however, for a long time, there was little incentive for its commercial planting, due to an extremely long juvenile stage, with approximately 10 years to maturity. Grafting is an effective approach to shorten the duration of vegetative growth, by which, pecan can start to bear fruits within 5–8 years. For the trees in the Juglandaceae family, grafting is more difficult in comparison to other fruit trees. In recent years, patch budding, one of the most commonly used grafting methods, conducted from July to September in China, has achieved about 90% of grafting success [1], which makes the large-scale cultivation of pecan possible. An in-depth understanding of the mechanism underlying successful grafting will help increase the production efficiency of pecan, as well as other trees in the future.
When grafting is performed, the grafted partners, scion, and rootstock are cut and joined together. Once the scion and rootstock come into intimate contact, an intricate structural and biochemical response occurs at the graft union for a successful graft. For woody trees, following the initial adhesion of grafted partners, the graft union undergoes two essential developmental processes: the formation of callus tissues, and the sufficient connection of functional vascular bundles between the scion and rootstock [2,3]. Therefore, graft union development is a process that involves cell division and differentiation at the graft junction.
Currently, reports regarding the molecular mechanism of graft union development are still limited. A cDNA-AFLP method was applied to investigate the gene expression in the graft process of hickory, and the research obtained 49 differentially expressed genes that were related to signal transduction, auxin transportation, metabolism, cell cycle, wound response, and cell wall synthesis [4]. In the hypocotyl grafts of Arabidopsis, changes in global gene expression were evaluated 24 h after grafting, and graft union development was revealed to involve signal transduction as well as cellular debris elimination [5]. In grapevine autografts, transcriptional changes were examined via whole genome microarray analysis, and the results revealed that graft union development triggered numerous gene expression changes related to wounding, cell wall modification, hormone signaling, and secondary metabolism [6]. A comparison the gene expression between the hetero- and autografts of grapevine indicated that genes involved in stress responses were up-regulated [7]. In recent years, RNA sequencing (RNA-seq) is a rapidly emerging transcriptome technology that can be performed without a reference genome. It has been employed to analyze the expression of mRNA and miRNA in hickory graft process, through which candidate genes involved in the auxin and cytokinin signaling were identified [8]; otherwise, a total of 12 candidate grafting-responsive miRNA were detected [9]. A comparative proteomic analysis of the hickory graft unions revealed that key enzymes involved in flavonoid biosynthesis were up-regulated 7 days after grafting [10].
Previously, we have paid attention to the morphological and proteomic changes in pecan homografts [11]. However, to the best of our knowledge, there have still been no reports describing the genes and gene networks underlying graft union development of pecan. In this study, we applied RNA-seq technology to construct mRNA libraries from the graft unions that were collected at 0, 8, 15, and 30 days after grafting, and analyzed the transcriptomic changes across the graft process.

2. Materials and Methods

2.1. Plant Material and Grafting Procedures

Pecan homografts were made in August using patch budding at the experimental orchard of Nanjing Forestry University (China). Graft unions (approximately 3 cm in length, the budding segment that includes the tissues of scion and the developing xylem of rootstock) were collected at 0, 8, 15, and 30 days after grafting, and immediately frozen in liquid nitrogen. The sampling timepoints were determined according to our histological analysis of the graft union developmental process in pecan homografts. In detail, the samples at 8 days, 15 days, and 30 days were selected for exploring the differentially expressed genes involved in the initial callus proliferation, massive callus proliferation accompanied by cambium establishment, and functional vascular bundles formation, respectively. Samples at 0 days were graft unions collected immediately from the scion and rootstock before grafting, and were used as controls. Three biological replicates were performed for each timepoint.

2.2. RNA Extraction, Library Construction, and Sequencing

Total RNA was isolated from the graft unions using the Universal Plant RNA Kit (BioTeke, Beijing, China) and treated with RNase-free DNase I (Takara) to degrade genomic DNA. RNA quality and quantity were monitored by a Nanodrop 1000 Spectrophotometer (Thermo Fisher Scientific, Wilmington, DE, USA) and Agilent 2100 Bioanalyzer (Agilent Technologies, Santa Clara, CA, USA). For each sample, about 3 μg of the total RNA that passed the quality examinations was used to prepare the cDNA library. Construction of sequencing libraries was performed by NEBNext® UltraTM RNA Library Prep Kit for Illumina® (NEB, Ipswich, MA, USA) according to the protocol. Briefly, the mRNA was enriched by oligo (dT)-attached magnetic beads and fragmented into short pieces, which were taken as templates for the first-strand and second-strand cDNA synthesis. Then, exonuclease/polymerase was used to convert the remaining overhangs into blunt ends. The resulting fragments were end-repaired by inserting an “A” base to the 3’ ends of the cDNA. NEBNext adapters with a hairpin loop structure were then ligated to the fragments. The library fragments were purified by an AMPure XP system (Beckman Coulter, Beverly, MA, USA) to select suitable cDNA fragments. Then, the products were amplified by PCR to create sequencing libraries. The constructed libraries were sequenced by an Illumina HiSeqTM 4000 platform (Biomarker Technology Company, Beijing, China). The sequencing raw data was deposited in the NCBI Sequence Read Archive (SRA) with the accession number SRP118757.

2.3. De Novo Assembly and Functional Annotation

After RNA sequencing, adapter sequences, poly-N reads, and low-quality reads from raw data were removed by in-house perl scripts to obtain clean reads. The resulting clean reads from all the samples were pooled for generating reference genes as far as possible. Trinity software with a k-mer length of 25 and other default parameters were used in the subsequent de novo assembly of transcriptome. Clean reads were assembled into contigs, and then further linked into transcripts through pair-end joining. The produced transcripts were clustered with a TGI clustering tool, and the longest transcripts were recognized as unigenes. For functional annotation, unigenes were compared against the following databases including NCBI non-redundant protein (Nr), Clusters of Orthologous Groups of proteins (COG), euKaryotic Orthologous Groups (KOG), Gene Ontology (GO), Kyoto Encyclopedia of Genes and Genomes (KEGG), Protein family (Pfam) and Swiss-Prot using the BLASTX program with E-value of 10−5.

2.4. Analysis of Differentially Expressed Genes (DEGs)

The clean reads sequenced from each sample were mapped back to the unigene library to calculate the abundance of unigenes. To quantify the gene expression level, FPKM (fragments per kilobase of exon per million mapped reads) was calculated in each sample by RSEM. Differential expression analysis was then performed using the DESeq R package for three comparisons (8 days vs. 0 days, 15 days vs. 0 days, and 30 days vs. 0 days). The false discovery rate (FDR) was applied to identify the p value threshold in multiple test and analysis. Only genes with FDR < 0.01 and more than two-fold change in expression between samples were considered as DEGs. GO enrichment analysis of DEGs was carried out by the topGO R package based on the hypergeometric test. Additionally, we used KOBAS software to test the enriched pathway of DEGs. GO terms and KEGG pathways with corrected p value ≤ 0.01 were recognized as significantly over-represented.

2.5. Validation of RNA-Seq Data by Quantitative Real-Time PCR (qRT-PCR)

RNA preparation with three biological replicates for each sample was conducted as described above. The first-strand cDNA synthesis was performed using a Prime-Script™ II First Strand cDNA synthesis kit (Takara Bio, Dalian, China) according to the manufacturer’s instructions. The primer sets for each unigene were designed by Primer Premier 5.0, and their sequences are listed in Table S1. qRT-PCR was carried out on an ABI 7500 Real-Time PCR System (Thermo Fisher Scientific, Inc., Waltham, MA, USA) with SYBR Premix Ex Taq™ II kit (Takara). Expression was calculated as 2−ΔΔCt and normalized to that of the reference gene Actin.

2.6. Detection of Hormones Content by ELISA

Samples were taken from the graft union at 0, 8, 15, and 30 days after grafting with three biological replicates. The contents of endogenous indole-3-acetic acid (IAA), zeatin riboside (ZR), and gibberellin (GA) were measured with the enzyme linked immunosorbent assay (ELISA). The hormone ELISA kits were developed from China Agricultural University, which have been validated with GC-MS and HPLC method. The determination of hormone content was performed as outlined by [12].

3. Results and Discussion

3.1. De Novo Assembly and Functional Annotation

To gain a comprehensive overview of transcriptome associated with graft union development in pecan, samples at different time points (0 days, 8 days, 15 days, and 30 days after grafting) with three biological replicates were subjected to illuminate sequencing. Raw reads were cleaned to generate a total of 312.08 million high-quality reads, encompassing 93.22 gigabase pairs with an average GC percentage of 46.41%. As a whole, all libraries showed good sequencing quality with Q30 (sequencing error rate less than 0.1%) more than 86.58% (Table S2). After sequence cleaning, reads from all samples were mixed to perform de novo assembly by Trinity software. Short reads were assembled into 140,455 transcripts with N50 length of 1905 bp (50% of total assembly length is contained within unigenes at least 1905 bp) and an average length of 1178 bp. There were 31,127 (22.16%) transcripts in the range between 1000 bp to 2000 bp, and 25,188 (17.93%) with length more than 2000 bp. All transcripts were subsequently clustered to yield 83,693 unigenes with an N50 length of 1350 bp and mean length of 892 bp. Among these, 13,698 (16.37%) unigenes were in the range of 1000–2000 bp, and 8364 (9.99%) exceeded 2000 bp (Table 1).
All the 83,693 unigenes were aligned with available protein databases using the BLASTx algorithm with E-value of 10−5. The results showed that there were 11,762 (14.05%) unigenes matched in the COG database, 23,260 (27.79%) in the GO database, 13,859 (16.56%) in the KEGG database, 21,612 (25.82%) in the KOG database, 25,909 (30.96%) in the Pfam database, 23,243 (27.77%) in the Swiss-Prot database, and 38,793 (46.35%) in the NR database. In total, there were 40,069 unigenes annotated in at least one database, accounting for 47.88% of all unigenes (Table 2). There was a relatively large portion of unigenes that had no significant hits to current known proteins, which might represent novel genes in pecan. For functional classification of the assembled unignenes, COG and GO annotation were carried out to gain the distributions of the functional categories (Figure S1).

3.2. Analysis of DEGs in the Graft Process of Pecan

Clean reads from the 12 libraries were aligned to the obtained unigenes, and quantified to calculate the expression levels by fragments per kilobase of transcript per million fragments mapped reads (FPKM). Based on the FPKM values of all unigenes, the correlations between biological replicates at each time point were analyzed. We found that there were strong correlations between the biological repetitions, with correlation coefficients over 0.90 (Figure S2).
According to the criteria of at least two-fold change and FDR < 0.01, a total of 3470 DEGs were discovered by analyzing 8 days/0 days, with 2154 up-regulated and 1316 down-regulated; 4942 DEGs were identified in the comparison of 15 days/0 days, with 2750 up-regulated and 2192 down-regulated; 9145 DEGs were found by comparing 30 days/0 days, with 3001 up-regulated and 6144 down-regulated. In total, 12,180 DEGs were identified during the grafting process, among which, 1499 genes were detected at all comparisons (Figure 1, Table S3). The number of DEGs in 30 days/0 days was greater than 8 days/0 days and 15 days/0 days, indicating the involvement of complex molecular responses during the developmental stage of vascular tissue formation.

3.3. Gene Ontology and Pathway Enrichment Analyses of DEGs

To elucidate the associated biological processes in which the DEGs were involved, the enrichment of GO terms was analyzed. For the ontology of biological process, there were 44, 28, and 38 GO enriched terms in the comparisons of 8 days/0 days, 15 days/0 days, and 30 days/0 days, respectively (Table S4, Figure 2). We found that “response to hormone”, “response to oxidative stress”, and “regulation of cell cycle” were simultaneously enriched in all the comparisons, suggesting the critical roles of these biological processes for a successful grafting. Interestingly, GO terms related to “plant-type secondary cell wall biogenesis” and “lignin metabolic process” were specially enriched in the 30 days vs. 0 days comparison, which were in good agreement with the developmental stage of vasculature formation at 30 days.
Additionally, KEGG enrichment analysis was performed to reveal the relevant metabolic pathways in which the DEGs participated. We identified 7, 5, and 2 significantly enriched pathways in 8 days/0 days, 15 days/0 days, and 30 days/0 days comparisons, respectively (Figure 3). Among those, the “phenylpropanoid biosynthesis” was the overlapping pathway that was identified in three comparisons, which was consistent with the significant role of this metabolic pathway during the grafting process [10,13].

3.4. Hormones Were Critical Regulators for Graft Union Development

During the grafting process, a block in auxin basipetal transport is produced due to vasculature damage, which leads to auxin accumulation at the graft junction. In our study, the content of auxin was increased distinctly at 8 days, 15 days, and 30 days after grafting (Figure 4A). Correspondingly, all the unigenes except one (c147017.graph_c0) encoding auxin influx carrier, and two unigenes (c129281.graph_c0 and c94763.graph_c0) encoding auxin efflux carrier were significantly up-regulated over the course of graft union development (Figure 5). Similarly, genes responsible for auxin transport were induced during the grafting process of grapevine [6] and hickory [8]. The accumulated auxin is indispensable for the regulation of callus proliferation and cambial activity [14,15]. A high level of auxin would release the transcriptional activity of auxin response factors (ARFs), which would induce the expression of genes that contain auxin responsive elements (AuxREs) in their promoter regions [16,17]. Previous studies have reported that ARF6 and ARF8 mutants showed cell division defects [15], and ARF5 mutants exhibited abnormality in vasculature development [18], indicating auxin signaling via ARFs appeared to be essential for the graft union development. In the present study, six unigenes encoding ARFs were differentially expressed (Figure 5). Intriguingly, three unigenes (c37236.graph_c0, c38752.graph_c0, and c114601.graph_c0) were up-regulated significantly at 8 days or 15 days, indicating a possible role in callus formation; while one unigene (c142339.graph_c0) was greatly up-regulated at 30 days, which suggested it might function in vasculature development.
Additionally, mounting evidence supports the involvement of cytokinin in cell division and vasculature differentiation [19,20,21,22]. Consistent with its role in graft union development, we found that the content of zeatin riboside (ZR), a major form of cytokinin, in graft junction was elevated significantly at 15 days and 30 days (Figure 4B). It has been reported that although auxin is capable of stimulating cell division, cytokinin is required for its full induction [23,24]. Therefore, massive callus proliferation at 15 days in this research might have resulted from the increased cytokinin as well as auxin. Cytokinin signal transduction is mediated via the two-component regulatory pathway to activate type-B ARR transcription factors [25]. Previous studies have showed that triple mutants of TYPE-B ARRs (ARR1, ARR10, ARR12) showed reduced callus formation [26], while overexpression of ARR1 enhanced callus formation [27]. The activated type-B ARRs are likely to be principal regulators of the cytokinin-induced callus proliferation. Three unigenes encoding type-B ARR protein were identified in our DEGs data, and one of them, c150807.graph_c1, was greatly up-regulated at 15 days (Figure 5), which might play an important role in callus formation.
More reports have suggested that gibberellin (GA) could trigger xylogenesis [28,29], which is probably important for the vascular bundle formation cross the grafted partners. This coincided with our biochemical analysis of GA content by ELISA, showing that GA was increased significantly at 30 days, while had no significant difference at other time points (Figure 4C). Accordingly, one GA synthesis gene, GA20ox, was highly up-regulated after grafting, and achieved its peak expression at 30 days, while three GA deactivating genes, GA2ox, were all significantly down-regulated at 30 days (Figure 5). GA signaling could promote the expression of genes involved in cell expansion, as well as secondary wall biosynthesis during xylem differentiation [30,31]. For the genes involved in GA signaling, we found that one unigene encoding GID1, a gibberellin receptor, was strikingly up-regulated at 30 days, and might take part in vasculature formation.

3.5. Genes Responsible for Callus Formation

Callus formation is a basic wound response to grafting, and the lack of callus production at the graft interface can lead to graft failure [32]. Genes involved in cell division are pivotal for callus formation [33]. Genome-wide transcriptomic study of callus initiation in Arabidopsis has revealed the up-regulation of various cell division related genes [34]. Cyclins, together with their catalytically active partners, cyclin dependent kinases (CDKs), are key regulators of cell cycle progression in eukaryotes. Three classes of cyclins (A-, B-, and D-type) exist in plants. Among them, D-type cyclin (CYCD) is a rate limiting factor for the G1/S transition in cell cycle, which plays an important role in driving the entry of cell cycle. It is usually considered as a sensor of external conditions that could be regulated by auxin and cytokinin [20,35,36]. In Arabidopsis, overexpression of a CYCD led to increased callus growth rate and callus induction frequency [37]. In this work, a considerable number of cyclins and CDKs were identified and all of them except one, c123028.graph_c0, were up-regulated across the grafting process (Figure 6), which might facilitate the activation of cellular proliferation.
In addition, we found one E2Fa gene, c146931.graph_c0, and one MYB3R-1 gene, c145971.graph_c0, were up-regulated significantly during the entire period of grafting (Figure 6). E2Fa is a transcription factor that drives the expression of genes required for the S-phase in cell cycle [35]. Transgenic Arabidopsis overexpressing E2FA could induce cell division in tissues already devoid of proliferation [38]. MYB3R-1 is an R1R2R3-type MYB transcription factor that aims at inducing genes required for the M-phase in cell cycle [39]. The up-regulated E2Fa along with MYB3R1 would probably facilitate cell cycle progression in this study. Besides the genes involved in cell cycle, we found that various genes responsible for nucleosome component synthesis (histone), DNA replication (DNA replication licensing factor MCM5 and MCM6), microtubule cytoskeleton organization (microtubule-associated protein RP/EB family and tubulin), and cytokinesis (kinesin-like protein) were generally up-regulated over the course of grafting (Figure 6), which might have participated in callus formation.

3.6. Genes Participated in Vascular Bundle Formation

Production of new vascular tissues permits the long-distance transport of nutrients between the grafting partners, and is recognized as a mark of successful grafting [32]. Based on our morphological observation of pecan graft union development, xylem is the main vascular tissue at 30 days [11]. The sequential developmental processes underlying wood formation include the promotion of vascular cambial activity, xylem differentiation, cell elongation, secondary cell wall thickening, and programmed cell death [40]. Plant hormones, including auxin and cytokinin, have been implicated in stimulating vascular cambium activity [40]. Previous studies have revealed the involvement of class III homeodomain-leucine zipper (HD-ZIP III) proteins in xylem differentiation [41,42]. It has been reported that the levels of HD-ZIP III determine the identity of xylem. Low levels of HD-ZIP III promote protoxylem identity, while high levels of HD-ZIP III promote metaxylem identity [43]. In our study, three down-regulated HD-ZIP III were identified at 30 days (Figure 7), which might function in xylem differentiation. Cells undergo significant enlargement following xylem differentiation. Enzymes such as expansions were required not only for cell growth, but also for the loosening of existing cell wall architecture during cell elongation [44]. As expected, most (six out of seven) unigenes encoding expansion in this study were up-regulated in both the callus proliferative phase and the vasculature formative phase (Figure 7). Tubulin, aside from its role in cell division, also plays a role in cell elongation by guiding nascent cellulose microfibrils deposition [45]. Expression of tubulin genes was also elevated at 30 days in our study (Figure 6). It is presumed that tubulin might also be involved in cell expansion during vascular development.
Following the completion of cell elongation, differentiating vascular cells go through the deposition of cellulose, hemicellulose, and lignin in the secondary cell wall. We identified various genes encoding the key biosynthetic enzymes of secondary cell wall components, and most of those genes were up-regulated, with the highest expression value at 30 days (Figure 7). CELLULOSE SYNTHASE, a gene implicated in cellulose synthesis, was found to be strongly expressed in the developing secondary xylem of Populus [46]. CINNAMOYL-COA REDUCTASE (CCR), CINNAMYL-ALCOHOL DEHYDROGENASE (CAD), and CAFFEOYL COA 3-O-METHYLTRANSFERASE (CCoAOMT) are the genes involved in the phenylpropanoid pathway, all of them critical for monolignol synthesis. The gene product of LACCASE is a polyphenol oxidase enzyme, which plays a critical role in lignin formation through inducing the oxidative polymerization of monolignols [47]. Mutations in LACCASE4 and 17 showed reduced lignin content in Arabidopsis [48]. IRREGULAR XYLEM 9 (IRX9), IRX10, and IRX15 are the genes that participate in hemicellulose synthesis.
Additionally, we identified candidate transcription factors involved in the transcriptional regulation of secondary cell wall deposition. Some NAC transcription factors are master regulators in controlling the entire developmental process of secondary cell wall synthesis [49]. Overexpression of NAC genes in plants induced secondary wall thickening in various tissues, while repression of their function suppressed secondary wall deposition [50,51]. Three particular secondary cell wall related NACs were identified, and all the NACs were strongly up-regulated at 30 days (Figure 7), which might imply an important role during vasculature differentiation. Secondary wall related R2R3-type MYB transcription factors are also important regulators, which have already been identified as transcriptional regulators of phenylpropanoid biosynthesis pathway [52]. In this study, nine candidate R2R3-type MYB genes were found to be differentially expressed (Figure 7). Among them, two unigenes, c128129.graph_c0 and c124957.graph_c0, were annotated as MYB46, which act as direct targets of NAC domain regulator, exhibiting up-regulation throughout the grafting process with peak value at 30 days (Figure 7). The induced MYB46 was not only able to activate synthesis of the lignin, but also the cellulose and hemicellulose [53]. One unigene, c110812.graph_c0, was annotated as MYB4, a gene that negatively regulated secondary cell wall formation [54], showing great down-regulation at 30 days (Figure 7). Ectopic overexpression of poplar PdMYB221, an ortholog of Arabidopsis MYB4, was reported to result in decreased thickness of cell walls [55]. Collectively, these DEGs might indicate the synthesis of secondary cell wall components during vasculature differentiation.
After fulfilling secondary cell wall deposition, developing vascular cells trigger programmed cell death (PCD) to digest the cellular contents. Hydrolytic enzymes, such as aspartic proteinase, cysteine proteinase, and nucleases (endonuclease, exonuclease, and ribonuclease) have been demonstrated to operate during xylogenesis [56,57,58], which were generally detected with increased expression at 30 days in our research (Figure 7). Metacaspases are a class of enzymes with structural similarity to animal caspases that could regulate the process of plant programmed cell death [59]. Analysis of microarray data revealed that the expression of an Arabidopsis METACASPASE 9 (AtMC9) homologue gene in Populus was up-regulated during xylem maturation [57]. In the present study, the great up-regulation of metacaspase at 30 days (Figure 7) might suggest its involvement in vascular bundle differentiation. Otherwise, we found one CYSTEINE PROTEINASE INHIBITOR was down-regulated drastically at 30 days (Figure 7), which probably indicated the strong activity of PCD during vascular bundle formation.

3.7. Genes Involved in Reactive Oxygen Species (ROS) Scavenging

During grafting, mechanical damage inevitably occurs at the graft interface. In higher plants, the production of ROS is a general event following wounding [60]. There are indications that non-successful grafts show signs of excess oxidative stress [61,62,63]. An efficient antioxidant defense system in plants might be an important factor in achieving successful grafting. In the current investigation, we found 19 DEGs that could scavenge the ROS, including genes encoding peroxidase (POD), catalase (CAT), ascorbate peroxidase (APX), cationic peroxidase, and peroxiredoxin (Prx), and most of them (13 out of 15) showed increased expression during the grafting process (Figure 8), which was presumably related to mitigating the ROS toxicity.

3.8. Validation of RNA-Seq Data by Real-Time RT-PCR (qRT-PCR)

We selected twelve genes that were predicted to be associated with hormone signaling, cell division, secondary cell wall formation, programmed cell death, and ROS scavenging to validate the sequencing data. Gene expression at different time points was relative to basal level (0 days). For a specific treatment, the fold changes of some genes in their expression detected by qRT-PCR and sequencing did not match exactly, however, the expression patterns of those selected genes were basically identical between the two data sets (Figure 9).

4. Conclusions

In this work, transcriptomic analysis was applied to explore the differentially expressed genes at the graft union during the pecan homograft process. A total of 12,180 DEGs were identified at the comparisons of 8 days/0 days, 15 days/0 days, and 30 days/0 days. Candidate genes that would participate in successful grafting were further analyzed. Based on our result, a suggested model for depicting the molecular mechanism of graft union development in pecan could be summarized in Figure 10. Upon grafting, signal transduction pathways including hormone (IAA, CK, and GA) signaling and other unknown signaling are activated. The activated signaling might induce the expression of various genes related to ROS scavenging, cell division, vasculature differentiation, cell elongation, secondary cell wall synthesis, and programmed cell death, resulting in a successful graft.

Supplementary Materials

The following are available online at https://0-www-mdpi-com.brum.beds.ac.uk/2073-4425/9/2/71/s1. Figure S1: Functional categories of assembled unigenes. (A) Distribution of number of annotated unigenes in COG database; (B) Distribution of number of annotated unigenes in GO database; Figure S2: Distribution of correlation co-efficiencies between each pair of samples; Table S1: The primer sequences of selected unigenes; Table S2: Summary of sequences analysis in all samples; Table S3: The identified DEGs in the graft process of pecan; Table S4: Significant enriched GO terms of DEGs between different developmental stages.

Acknowledgments

The authors appreciate the financial support from the SanXin project of Jiangsu province (LYSX[2016]44), the state bureau of forestry 948 project (2015-4-16) and the Priority Academic Program Development of Jiangsu Higher Education Institutions (PAPD).

Author Contributions

F.P. conceived and designed the study. Z.M. performed the entire data analysis, carried out qRT-PCR and wrote the manuscript. G.F. and W.S. participated in hormone detection. Z.L. was involved in sample collection. All authors read and approved the final manuscript.

Conflicts of Interest

The authors declare that they have no competing interests.

References

  1. Zhang, R.; Peng, F.; Li, Y. Pecan production in China. Sci. Hortic. 2015, 197, 719–727. [Google Scholar] [CrossRef]
  2. Asante, A.; Barnett, J. Graft union formation in mango (Mangifera indica L.). J. Hortic. Sci. 1997, 72, 781–790. [Google Scholar] [CrossRef]
  3. Yang, Z.J.; Feng, J.L.; Chen, H. Study on the anatomical structures in development of the nurse seed grafted union of Camellia oleifera. Plant Sci. J. 2013, 31, 313–320. [Google Scholar] [CrossRef]
  4. Zheng, B.S.; Chu, H.L.; Jin, S.H.; Huang, Y.J.; Wang, Z.J.; Chen, M.; Huang, J.Q. cDNA-AFLP analysis of gene expression in hickory (Carya carthayensis) during graft process. Tree Physiol. 2009, 30, 297–303. [Google Scholar] [CrossRef] [PubMed]
  5. Yin, H.; Yan, B.; Sun, J.; Jia, P.; Zhang, Z.; Yan, X.; Chai, J.; Ren, Z.; Zheng, G.; Liu, H. Graft-union development: A delicate process that involves cell–cell communication between scion and stock for local auxin accumulation. J. Exp. Bot. 2012, 63, 4219–4232. [Google Scholar] [CrossRef] [PubMed]
  6. Cookson, S.J.; Clemente Moreno, M.J.; Hevin, C.; Nyamba Mendome, L.Z.; Delrot, S.; Trossatmagnin, C.; Ollat, N. Graft union formation in grapevine induces transcriptional changes related to cell wall modification, wounding, hormone signalling, and secondary metabolism. J. Exp. Bot. 2013, 64, 2997–3008. [Google Scholar] [CrossRef] [PubMed]
  7. Cookson, S.J.; Moreno, M.J.C.; Hevin, C.; Mendome, L.Z.N.; Delrot, S.; Magnin, N.; Trossat-Magnin, C.; Ollat, N. Heterografting with nonself rootstocks induces genes involved in stress responses at the graft interface when compared with autografted controls. J. Exp. Bot. 2014, 65, 2473–2481. [Google Scholar] [CrossRef] [PubMed]
  8. Qiu, L.; Jiang, B.; Fang, J.; Shen, Y.; Fang, Z.; RM, S.K.; Yi, K.; Shen, C.; Yan, D.; Zheng, B. Analysis of transcriptome in hickory (Carya carthayensis), and uncover the dynamics in the hormonal signaling pathway during graft process. BMC Genom. 2016, 17, 935. [Google Scholar] [CrossRef] [PubMed]
  9. Sima, X.; Jiang, B.; Fang, J.; He, Y.; Fang, Z.; Saravana Kumar, K.M.; Ren, W.; Qiu, L.; Chen, X.; Zheng, B. Identification by deep sequencing and profiling of conserved and novel hickory micrornas involved in the graft process. Plant Biotechnol. Rep. 2015, 9, 115–124. [Google Scholar] [CrossRef]
  10. Xu, D.; Yuan, H.; Tong, Y.; Zhao, L.; Qiu, L.; Guo, W.; Shen, C.; Liu, H.; Yan, D.; Zheng, B. Comparative proteomic analysis of the graft unions in hickory (Carya carthayensis) provides insights into response mechanisms to grafting process. Front. Plant Sci. 2017, 8. [Google Scholar] [CrossRef] [PubMed]
  11. Mo, Z.; He, H.; Su, W.; Peng, F. Analysis of differentially accumulated proteins associated with graft union formation in pecan (Carya illinoensis). Sci. Hortic. 2017, 224, 126–134. [Google Scholar] [CrossRef]
  12. Yang, Y.M.; Xu, C.N.; Wang, B.M.; Jia, J.Z. Effects of plant growth regulators on secondary wall thickening of cotton fibres. Plant Growth Regul. 2001, 35, 233–237. [Google Scholar]
  13. Pina, A.; Errea, P. Differential induction of phenylalanine ammonia-lyase gene expression in response to in vitro callus unions of Prunus spp. J. Plant Physiol. 2008, 165, 705–714. [Google Scholar] [CrossRef] [PubMed]
  14. Melnyk, C.W.; Schuster, C.; Leyser, O.; Meyerowitz, E.M. A developmental framework for graft formation and vascular reconnection in Arabidopsis thaliana. Curr. Biol. 2015, 25, 1306–1318. [Google Scholar] [CrossRef] [PubMed]
  15. Pitaksaringkarn, W.; Ishiguro, S.; Asahina, M.; Satoh, S. ARF6 and ARF8 contribute to tissue reunion in incised Arabidopsis inflorescence stems. Plant Biotechnol. 2014, 31, 49–53. [Google Scholar] [CrossRef]
  16. Li, S.B.; Xie, Z.Z.; Hu, C.G.; Zhang, J.Z. A review of auxin response factors (ARFs) in plants. Front. Plant Sci. 2016, 7. [Google Scholar] [CrossRef] [PubMed]
  17. Mattsson, J.; Ckurshumova, W.; Berleth, T. Auxin signaling in Arabidopsis leaf vascular development. Plant Physiol. 2003, 131, 1327–1339. [Google Scholar] [CrossRef] [PubMed]
  18. Hardtke, C.S.; Berleth, T. The Arabidopsis gene MONOPTEROS encodes a transcription factor mediating embryo axis formation and vascular development. EMBO J. 1998, 17, 1405–1411. [Google Scholar] [CrossRef] [PubMed]
  19. Aloni, R.; Aloni, E.; Langhans, M.; Ullrich, C. Role of cytokinin and auxin in shaping root architecture: Regulating vascular differentiation, lateral root initiation, root apical dominance and root gravitropism. Ann. Bot. 2006, 97, 883–893. [Google Scholar] [CrossRef] [PubMed]
  20. Ikeuchi, M.; Sugimoto, K.; Iwase, A. Plant callus: Mechanisms of induction and repression. Plant Cell 2013, 25, 3159–3173. [Google Scholar] [CrossRef] [PubMed]
  21. Nieminen, K.; Immanen, J.; Laxell, M.; Kauppinen, L.; Tarkowski, P.; Dolezal, K.; Tähtiharju, S.; Elo, A.; Decourteix, M.; Ljung, K. Cytokinin signaling regulates cambial development in poplar. Proc. Natl. Acad. Sci. USA 2008, 105, 20032–20037. [Google Scholar] [CrossRef] [PubMed]
  22. Hejátko, J.; Ryu, H.; Kim, G.T.; Dobešová, R.; Choi, S.; Choi, S.M.; Souček, P.; Horák, J.; Pekárová, B.; Palme, K. The histidine kinases CYTOKININ-INDEPENDENT1 and ARABIDOPSIS HISTIDINE KINASE2 and 3 regulate vascular tissue development in Arabidopsis shoots. Plant Cell 2009, 21, 2008–2021. [Google Scholar] [CrossRef] [PubMed]
  23. Del Pozo, J.C.; Lopez-Matas, M.; Ramirez-Parra, E.; Gutierrez, C. Hormonal control of the plant cell cycle. Physiol. Plant. 2005, 123, 173–183. [Google Scholar] [CrossRef]
  24. Perrot-Rechenmann, C. Cellular responses to auxin: Division versus expansion. Cold Spring Harb. Perspect. Biol. 2010, 2. [Google Scholar] [CrossRef] [PubMed]
  25. Pils, B.; Heyl, A. Unraveling the evolution of cytokinin signaling. Plant Physiol. 2009, 151, 782–791. [Google Scholar] [CrossRef] [PubMed]
  26. Mason, M.G.; Mathews, D.E.; Argyros, D.A.; Maxwell, B.B.; Kieber, J.J.; Alonso, J.M.; Ecker, J.R.; Schaller, G.E. Multiple type-B response regulators mediate cytokinin signal transduction in Arabidopsis. Plant Cell 2005, 17, 3007–3018. [Google Scholar] [CrossRef] [PubMed]
  27. Sakai, H.; Honma, T.; Aoyama, T.; Sato, S.; Kato, T.; Tabata, S.; Oka, A. ARR1, a transcription factor for genes immediately responsive to cytokinins. Science 2001, 294, 1519–1521. [Google Scholar] [CrossRef] [PubMed]
  28. Mauriat, M.; Moritz, T. Analyses of GA20ox- and GID1-over-expressing aspen suggest that gibberellins play two distinct roles in wood formation. Plant J. 2009, 58, 989–1003. [Google Scholar] [CrossRef] [PubMed]
  29. Ragni, L.; Nieminen, K.; Pacheco-Villalobos, D.; Sibout, R.; Schwechheimer, C.; Hardtke, C.S. Mobile gibberellin directly stimulates Arabidopsis hypocotyl xylem expansion. Plant Cell 2011, 23, 1322. [Google Scholar] [CrossRef] [PubMed]
  30. Guo, H.; Wang, Y.; Liu, H.; Hu, P.; Jia, Y.; Zhang, C.; Wang, Y.; Gu, S.; Yang, C.; Wang, C. Exogenous GA3 application enhances xylem development and induces the expression of secondary wall biosynthesis related genes in Betula platyphylla. Int. J. Mol. Sci. 2015, 16, 22960–22975. [Google Scholar] [CrossRef] [PubMed]
  31. Israelsson, M.; Sundberg, B.; Moritz, T. Tissue-specific localization of gibberellins and expression of gibberellin-biosynthetic and signaling genes in wood-forming tissues in aspen. Plant J. 2005, 44, 494–504. [Google Scholar] [CrossRef] [PubMed]
  32. Pina, A.; Errea, P. A review of new advances in mechanism of graft compatibility–incompatibility. Sci. Hortic. 2005, 106, 1–11. [Google Scholar] [CrossRef]
  33. Chen, C.C.; Fu, S.F.; Lee, Y.I.; Lin, C.Y.; Lin, W.C.; Huang, H.J. Transcriptome analysis of age-related gain of callus-forming capacity in Arabidopsis hypocotyls. Plant Cell Physiol. 2012, 53, 1457–1469. [Google Scholar] [CrossRef] [PubMed]
  34. Xu, K.; Liu, J.; Fan, M.; Xin, W.; Hu, Y.; Xu, C. A genome-wide transcriptome profiling reveals the early molecular events during callus initiation in Arabidopsis multiple organs. Genomics 2012, 100, 116–124. [Google Scholar] [CrossRef] [PubMed]
  35. Fehér, A.; Magyar, Z. Coordination of cell division and differentiation in plants in comparison to animals. Acta Biol. Szeged. 2015, 59, 275–289. [Google Scholar]
  36. Riou-Khamlichi, C.; Huntley, R.; Jacqmard, A.; Murray, J.A. Cytokinin activation of Arabidopsis cell division through a D-type cyclin. Science 1999, 283, 1541–1544. [Google Scholar] [CrossRef] [PubMed]
  37. Cockcroft, C.E.; den Boer, B.G.; Healy, J.S.; Murray, J.A. Cyclin D control of growth rate in plants. Nature 2000, 405, 575–579. [Google Scholar] [PubMed]
  38. De, V.L.; Beeckman, T.; Beemster, G.T.; De, A.E.J.; Ormenese, S.; Maes, S.; Naudts, M.; Van, D.S.E.; Jacqmard, A.; Engler, G. Control of proliferation, endoreduplication and differentiation by the Arabidopsis E2Fa-DPa transcription factor. EMBO J. 2002, 21, 1360–1368. [Google Scholar]
  39. Magyar, Z.; Bögre, L.; Ito, M. Dreams make plant cells to cycle or to become quiescent. Curr. Opin. Plant Biol. 2016, 34, 100–106. [Google Scholar] [CrossRef] [PubMed]
  40. Ye, Z.; Zhong, R. Molecular control of wood formation in trees. J. Exp. Bot. 2015, 66, 4119–4131. [Google Scholar] [CrossRef] [PubMed]
  41. Du, J.; Miura, E.; Robischon, M.; Martinez, C.; Groover, A. The Populus Class III HD ZIP transcription factor Popcorona affects cell differentiation during secondary growth of woody stems. PLoS ONE 2011, 6, e17458. [Google Scholar] [CrossRef] [PubMed]
  42. Ilegems, M.; Douet, V.; Meylanbettex, M.; Uyttewaal, M.; Brand, L.; Bowman, J.L.; Stieger, P.A. Interplay of auxin, KANADI and Class III HD-ZIP transcription factors in vascular tissue formation. Development 2010, 137, 975–984. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  43. De, R.B.; Mähönen, A.P.; Helariutta, Y.; Weijers, D. Plant vascular development: From early specification to differentiation. Nat. Rev. Mol. Cell Biol. 2015, 17, 30–40. [Google Scholar]
  44. Ye, Z.H. Vascular tissue differentiation and pattern formation in plants. Annu. Rev. Plant Biol. 2002, 53, 183–202. [Google Scholar] [CrossRef] [PubMed]
  45. Mendu, N.; Silflow, C.D. Elevated levels of tubulin transcripts accompany the GA3-induced elongation of oat internode segments. Plant Cell Physiol. 1993, 34, 973–983. [Google Scholar]
  46. Suzuki, S.; Li, L.; Sun, Y.H.; Chiang, V.L. The cellulose synthase gene superfamily and biochemical functions of xylem-specific cellulose synthase-like genes in Populus trichocarpa. Plant Physiol. 2006, 142, 1233–1245. [Google Scholar] [CrossRef] [PubMed]
  47. Zhao, Q.; Nakashima, J.; Chen, F.; Yin, Y.; Fu, C.; Yun, J.; Shao, H.; Wang, X.; Wang, Z.-Y.; Dixon, R.A. Laccase is necessary and nonredundant with peroxidase for lignin polymerization during vascular development in Arabidopsis. Plant Cell 2013, 25, 3976–3987. [Google Scholar] [CrossRef] [PubMed]
  48. Berthet, S.; Demont-Caulet, N.; Pollet, B.; Bidzinski, P.; Cézard, L.; Le Bris, P.; Borrega, N.; Hervé, J.; Blondet, E.; Balzergue, S. Disruption of Laccase4 and 17 results in tissue-specific alterations to lignification of Arabidopsis thaliana stems. Plant Cell 2011, 23, 1124–1137. [Google Scholar] [CrossRef] [PubMed]
  49. Zhong, R.; Lee, C.; Zhou, J.; McCarthy, R.L.; Ye, Z.-H. A battery of transcription factors involved in the regulation of secondary cell wall biosynthesis in Arabidopsis. Plant Cell 2008, 20, 2763–2782. [Google Scholar] [CrossRef] [PubMed]
  50. Mitsuda, N.; Iwase, A.; Yamamoto, H.; Yoshida, M.; Seki, M.; Shinozaki, K.; Ohme-Takagi, M. NAC transcription factors, NST1 and NST3, are key regulators of the formation of secondary walls in woody tissues of Arabidopsis. Plant Cell 2007, 19, 270–280. [Google Scholar] [CrossRef] [PubMed]
  51. Zhou, J.; Zhong, R.; Ye, Z. Arabidopsis NAC domain proteins, VND1 to VND5, are transcriptional regulators of secondary wall biosynthesis in vessels. PLoS ONE 2014, 9, e105726. [Google Scholar] [CrossRef] [PubMed]
  52. Stracke, R.; Werber, M.; Weisshaar, B. The R2R3-MYB gene family in Arabidopsis thaliana. Curr. Opin. Plant Biol. 2001, 4, 447–456. [Google Scholar] [CrossRef]
  53. Yang, J.H.; Wang, H. Molecular mechanisms for vascular development and secondary cell wall formation. Front. Plant Sci. 2016, 7. [Google Scholar] [CrossRef] [PubMed]
  54. Zhao, Q.; Dixon, R.A. Transcriptional networks for lignin biosynthesis: More complex than we thought? Trends Plant Sci. 2011, 16, 227–233. [Google Scholar] [CrossRef] [PubMed]
  55. Tang, X.; Zhuang, Y.; Qi, G.; Wang, D.; Liu, H.; Wang, K.; Chai, G.; Zhou, G. Poplar PdMYB221 is involved in the direct and indirect regulation of secondary wall biosynthesis during wood formation. Sci. Rep. 2015, 5. [Google Scholar] [CrossRef] [PubMed]
  56. Iliev, I.; Savidge, R. Proteolytic activity in relation to seasonal cambial growth and xylogenesis in Pinus banksiana. Phytochemistry 1999, 50, 953–960. [Google Scholar] [CrossRef]
  57. Courtois-Moreau, C.L.; Pesquet, E.; Sjödin, A.; Muñiz, L.; Bollhöner, B.; Kaneda, M.; Samuels, L.; Jansson, S.; Tuominen, H. A unique program for cell death in xylem fibers of Populus stem. Plant J. 2009, 58, 260–274. [Google Scholar] [CrossRef] [PubMed]
  58. Iakimova, E.T.; Woltering, E.J. Xylogenesis in zinnia (Zinnia elegans) cell cultures: Unravelling the regulatory steps in a complex developmental programmed cell death event. Planta 2017, 245, 681–705. [Google Scholar] [CrossRef] [PubMed]
  59. Petzold, H.E.; Zhao, M.; Beers, E.P. Expression and functions of proteases in vascular tissues. Physiol. Plant. 2012, 145, 121–129. [Google Scholar] [CrossRef] [PubMed]
  60. Ślesak, I.; Ślesak, H.; Libik, M.; Miszalski, Z. Antioxidant response system in the short-term post-wounding effect in Mesembryanthemum crystallinum leaves. J. Plant Physiol. 2008, 165, 127–137. [Google Scholar] [CrossRef] [PubMed]
  61. Aloni, B.; Karni, L.; Deventurero, G.; Levin, Z.; Cohen, R.; Katzir, N.; Lotan-Pompan, M.; Edelstein, M.; Aktas, H.; Turhan, E. Physiological and biochemical changes at the rootstock-scion interface in graft combinations between Cucurbita rootstocks and a melon scion. J. Hortic. Sci. Biotech. 2008, 83, 777–783. [Google Scholar] [CrossRef]
  62. Irisarri, P.; Binczycki, P.; Errea, P.; Martens, H.J.; Pina, A. Oxidative stress associated with rootstock–scion interactions in pear/quince combinations during early stages of graft development. J. Plant Physiol. 2015, 176, 25–35. [Google Scholar] [CrossRef] [PubMed]
  63. Xu, Q.; Guo, S.R.; Li, H.; Du, N.S.; Shu, S.; Sun, J. Physiological aspects of compatibility and incompatibility in grafted cucumber seedlings. J. Am. Soc. Hort. Sci. 2015, 140, 299–307. [Google Scholar]
Figure 1. The differentially expressed genes (DEGs) in different comparisons (8 days/0 days, 15 days/0 days, and 30 days/0 days) during graft union development in pecan homografts.
Figure 1. The differentially expressed genes (DEGs) in different comparisons (8 days/0 days, 15 days/0 days, and 30 days/0 days) during graft union development in pecan homografts.
Genes 09 00071 g001
Figure 2. GO enrichment of DEGs during the graft process. (A) Significantly enriched GO terms between 8 days and 0 days; (B) significantly enriched GO terms between 15 days and 0 days; (C) significantly enriched GO terms between 30 days and 0 days. Bubbles represent the significant GO terms, and the bubble color gradient represent the magnitude of enrichment corresponding to q-values.
Figure 2. GO enrichment of DEGs during the graft process. (A) Significantly enriched GO terms between 8 days and 0 days; (B) significantly enriched GO terms between 15 days and 0 days; (C) significantly enriched GO terms between 30 days and 0 days. Bubbles represent the significant GO terms, and the bubble color gradient represent the magnitude of enrichment corresponding to q-values.
Genes 09 00071 g002
Figure 3. Significant enriched KEGG pathways during the graft process. (A) Comparison of 8 days/0 days; (B) comparison of 15 days/0 days; (C) comparison of 30 days/0 days.
Figure 3. Significant enriched KEGG pathways during the graft process. (A) Comparison of 8 days/0 days; (B) comparison of 15 days/0 days; (C) comparison of 30 days/0 days.
Genes 09 00071 g003
Figure 4. Determination of the contents of hormone at the graft unions during the pecan grafting process. (A) Indole acetic acid (IAA); (B) zeatin riboside (ZR); and (C) gibberellic acid (GA) at different timepoints. * indicates the significant differences (p < 0.05) between the specific timepoints and the basal level (0 days).
Figure 4. Determination of the contents of hormone at the graft unions during the pecan grafting process. (A) Indole acetic acid (IAA); (B) zeatin riboside (ZR); and (C) gibberellic acid (GA) at different timepoints. * indicates the significant differences (p < 0.05) between the specific timepoints and the basal level (0 days).
Genes 09 00071 g004
Figure 5. Expression patterns of DEGs involved in hormone signaling. The values of log2 fold change are shown in the heat map.
Figure 5. Expression patterns of DEGs involved in hormone signaling. The values of log2 fold change are shown in the heat map.
Genes 09 00071 g005
Figure 6. Expression patterns of DEGs involved in and callus formation. The values of log2 fold change are shown in the heat map.
Figure 6. Expression patterns of DEGs involved in and callus formation. The values of log2 fold change are shown in the heat map.
Genes 09 00071 g006
Figure 7. Expression profiles of DEGs involved in vascular bundle formation. Gene expression values were normalized to z-score. The values of log2 fold change are shown in the heat map.
Figure 7. Expression profiles of DEGs involved in vascular bundle formation. Gene expression values were normalized to z-score. The values of log2 fold change are shown in the heat map.
Genes 09 00071 g007
Figure 8. Expression profiles of DEGs involved in ROS scavenging. The values of log2 fold change are shown in the heat map.
Figure 8. Expression profiles of DEGs involved in ROS scavenging. The values of log2 fold change are shown in the heat map.
Genes 09 00071 g008
Figure 9. Validation RNA-seq data by real-time quantitative RT-PCR. The expression levels at 0 days are considered as 1.
Figure 9. Validation RNA-seq data by real-time quantitative RT-PCR. The expression levels at 0 days are considered as 1.
Genes 09 00071 g009
Figure 10. A putative molecular model of successful grafting in pecan.
Figure 10. A putative molecular model of successful grafting in pecan.
Genes 09 00071 g010
Table 1. Summary for the graft union transcriptome.
Table 1. Summary for the graft union transcriptome.
TranscriptUnigene
Total number140,45583,693
Total length165,440,800 bp74,679,367 np
N50 length1905 bp1350 bp
Mean length1178 bp892 bp
200–300 bp18,665 (13.29%)15,637 (18.68%)
300–500 bp28,705 (20.44%)21,661 (25.88%)
500–1000 bp36,770 (26.18%)24,333 (29.07%)
1000–2000 bp31,127 (22.16%)13,698 (16.37%)
2000+ bp25,188 (17.93%)8364 (9.99%)
Table 2. Summary for the annotation of unigenes.
Table 2. Summary for the annotation of unigenes.
Annotated DatabasesUnigene NumberPercentage (%)300 nt ≤ Length < 1000 ntLength ≥ 1000 nt
COG11,76214.0537476644
GO23,26027.79902810,879
KEGG13,85916.5655676557
KOG21,61225.82818410,759
Pfam25,90930.96883814,658
Swiss-Prot23,24327.77871112,071
NR38,79346.3515,81617,494
Annotated in at least one database40,06947.8816,46317,751

Share and Cite

MDPI and ACS Style

Mo, Z.; Feng, G.; Su, W.; Liu, Z.; Peng, F. Transcriptomic Analysis Provides Insights into Grafting Union Development in Pecan (Carya illinoinensis). Genes 2018, 9, 71. https://0-doi-org.brum.beds.ac.uk/10.3390/genes9020071

AMA Style

Mo Z, Feng G, Su W, Liu Z, Peng F. Transcriptomic Analysis Provides Insights into Grafting Union Development in Pecan (Carya illinoinensis). Genes. 2018; 9(2):71. https://0-doi-org.brum.beds.ac.uk/10.3390/genes9020071

Chicago/Turabian Style

Mo, Zhenghai, Gang Feng, Wenchuan Su, Zhuangzhuang Liu, and Fangren Peng. 2018. "Transcriptomic Analysis Provides Insights into Grafting Union Development in Pecan (Carya illinoinensis)" Genes 9, no. 2: 71. https://0-doi-org.brum.beds.ac.uk/10.3390/genes9020071

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop