Transcriptome and proteome profiling of adventitious root development in hybrid larch (Larix kaempferi × Larix olgensis)
BMC Plant Biology volume 14, Article number: 305 (2014)
Hybrids of larch (Larix kaempferi × Larix olgensis) are important afforestation species in northeastern China. They are routinely propagated via rooted stem cuttings. Despite the importance of rooting, little is known about the regulation of adventitious root development in larch hybrids. 454 GS FLX Titanium technology represents a new method for characterizing the transcriptomes of non-model species. This method can be used to identify differentially expressed genes, and then two-dimensional difference gel electrophoresis (2D-DIGE) and matrix-assisted laser desorption-ionization time-of-flight mass spectrometry (MALDI-TOF/TOF MS) analyses can be used to analyze their corresponding proteins. In this study, we analyzed semi-lignified cuttings of two clones of L. kaempferi × L. olgensis with different rooting capacities to study the molecular basis of adventitious root development.
We analyzed two clones; clone 25-5, with strong rooting capacity, and clone 23-12, with weak rooting capacity. We constructed four cDNA libraries from 25-5 and 23-12 at two development stages. Sequencing was conducted using the 454 pyrosequencing platform. A total of 957832 raw reads was produced; 95.07% were high-quality reads, and were assembled into 45137 contigs and 61647 singletons. The functions of the unigenes, as indicated by their Gene Ontology annotation, included diverse roles in the molecular functions, biological processes, and cellular component categories. We analyzed 75 protein spots (-fold change ≥2, P ≤0.05) by 2D-DIGE, and identified the differentially expressed proteins using MALDI-TOF/TOF MS. A joint analysis of transcriptome and proteome showed genes related to two pathways, polyamine synthesis and stress response, might play an important role on adventitious root development.
These results provide fundamental and important information for research on the molecular mechanism of adventitious root development. We also demonstrated for the first time the combined use of two important technologies as a powerful approach to advance research on non-model, but otherwise important, larch species.
Larches are very useful for afforestation because of their fast juvenile growth . Among the various larches, hybrids of Larix kaempferi × Larix olgensis are widely planted in the mountainous regions of northeastern China for their high-quality timber, their pest resistance, and their amenity value . The use of rooted stem cuttings is becoming the most popular method to propagate these hybrids. Adventitious root formation is a key step for stem cuttings, and it is a complex process affected by many factors -. Despite the ecological and economical importance of larch, little genomic research has been conducted on this species. It would be very useful to understand the biological process of adventitious root formation and development in larch.
In our previous study, two clones of L. kaempferi × L. olgensis, 25-5 and 23-12, are obviously different on the rooting rates, which were 100.0%, 94.7%, 93.3%, 96.3% and 25.1%, 19.6%, 21.5%, 22.9%, in 2001, 2002, 2003, and 2007, respectively. So, the two clones are good materials for studying the biological process of adventitious root formation. Additionally, we had examined the development process of adventitious roots in semi-lignified cuttings of hybrid larch , which consisted of three key stages: 14-18 days after cutting (DAC) for cell dedifferentiation and division; 25-35 DAC for meristem formation and development; 50-55 DAC for root formation and elongation. The early stages, 14 and 25 DAC, were crucial for the expression and regulation of rooting-related genes. These anatomical and physiological researches on adventitious roots development in hybrid larch make the further study on the molecular mechanisms governing adventitious root development feasible.
Although adventitious root formation is the key developmental process for asexual propagation of most plants, the molecular mechanisms of adventitious root formation are still poorly understood , because adventitious root formation is a complex quantitative genetic trait regulated by both environmental and endogenous factors. Only a limited number of molecular studies of adventitious root formation have been performed. Brinker et al.  identified 220 differentially expressed genes during different developmental stages of adventitious root formation in Pinus contorta hypocotyls which were treated with the auxin indole-3-butyric acid. Butler et al.  focused on gene expression patterns during adventitious root formation, and a number of mRNAs during adventitious root formation in apple were identified. Sorin et al.  identified some molecular markers which positively or negatively correlated with adventitious root formation by a proteomic analysis of different mutant genotypes of Arabidopsis. Chen et al.  provided a method suitable to study de novo root organogenesis. Ahkami et al.  identify specific genes determining the initiation and formation of adventitious roots by a microarray-based transcriptome analysis in the stem base of the cuttings of Petunia hybrida (line W115). Subramaniyam et al.  obtained a comprehensive transcript expression profiling for adventitious roots of Panax ginseng Meyer. However, very few genes were characterized for adventitious root formation of larch, due to the limited genome information.
RNA sequencing represents a powerful and rapid method for obtaining functional information on a genome-wide scale, especially for the large genomes and non-model species. Wu et al.  studied gene expression in North American ginseng root samples at seven developmental stages by 454 sequencing. Their results suggested that ginsenoside biosynthesis occurs at distinct stages of development. Firon et al.  identified the molecular mechanisms involved in the initiation of storage root formation in sweet potato by RNA sequencing. Alagna et al.  enriched the very small amount of sequence data currently available for olive, and identified genes involved in fruit quality via 454 pyrosequencing. Extensive analyses of transcriptome profiles - have identified functional genes with roles in specific biological processes. However, changes in mRNA transcript levels do not always reflect changes in protein levels, and there are many ways in which cellular events are regulated at the protein level. Therefore, a complementary protein analysis is necessary to clarify the molecular events at particular stages of development. Two-dimensional difference gel electrophoresis (2D-DIGE) combined with matrix-assisted laser desorption-ionization time-of-flight mass spectrometry (MALDI-TOF/TOF MS) has many advantages for protein analysis. Li et al.  detected significant changes in the abundance of certain proteins in cucumber roots in response to hypoxia using 2D-DIGE. Yan et al.  used a proteomic analysis to study the responses to growth restriction in maize plants supplied with sufficient nitrogen. Wang et al.  identified proteins involved in brassinosteroid regulation of rice growth using 2D-DIGE followed by MS. The results of those studies illustrated that transcriptomic or proteomic analyses are effective methods to find genes and proteins involved in regulating biological processes. Yang et al.  combined use of proteomic and transcriptomic analyses of M14 leaves and roots under NaCl treatment in sugar beet, these analyses revealed candidate genes and proteins for detailed functional characterization.
In this study, to gain a better understanding of adventitious root development in hybrid larch, we evaluated the transcriptomic profiles of two clones with different rooting capacity at two important stages of adventitious root development. Some of the genes identified as being differentially expressed between the two clones were further analyzed at the protein level. These data can remarkably expand the transcriptome and proteome database of larch, and provides a valuable platform to understand the molecular basis of regulation of adventitious root development in larch and other conifers.
Results and discussion
Compared with other techniques, the 454 pyrosequencing technique is a better alternative for de novo sequencing and analysis of the transcriptomes of non-model plants, because of the high number of reads generated per run and the low sequencing error rate in the obtained contigs . This technology has been used to construct a number of expressed sequence tag (EST) libraries from plants including including maize , chestnut , olive , North American ginseng , and sweet potato . Zhang et al.  used 454 sequencing to analyze Japanese larch during somatic embryogenesis. In this study, we constructed four cDNA libraries from two clones at two important stages of adventitious root development. Each of the cDNA libraries was sequenced using the 454 GS FLX Titanium platform. A total of 957,832 raw reads were produced from the four cDNA libraries. Most of the raw reads were between 450 and 500 bp long. After trimming adapter sequences and removing short sequences of less than 50 bases, a total of 910,607 high-quality reads were generated from the four cDNA libraries. Therefore, 95.07% of raw reads were useful. The length distribution of these high-quality reads is shown in Figure 1. Approximately 90% of the reads were between 300 and 500 bp in length.
These 910,607 high-quality reads were assembled into 45137 contigs (a contig is a set of overlapping RNA segments derived from a single genetic source) and 61647 singletons (A single gene that identifies the comparable region in a target genome, but does not cluster with other genes to form a segment is called a singleton). The distribution of contig lengths is shown in Figure 2; 68.79% (31048) of contigs were 300–800 bp long, and 24.06% (10859) of contigs were 400–500 bp long. The median contig length was 739 bp, while the average contig length was 663 bp. The lengths of unigenes were similar to those obtained in previous transcriptome analyses of Artemisia , chestnut , American ginseng , and Cucurbita pepo . The data shown here represent the first large-scale effort to generate cDNA resources and analyze the transcriptomes of hybrid larch during adventitious root development. These resources have been made publically available; the raw data have been submitted to the Short Read Archive (SRA) division of the Genbank repository (accession number SRP015266). The transcriptomic sequences should be useful for gene discovery in larch, and will be helpful for elucidating the genomic basis of adventitious root development.
Functional annotation and classification
All unique sequences were first compared with sequences in the non-redundant database (nt) of the NCBI using the BLASTN algorithm. Then, they were compared with sequences in the two major public protein databases (Methods) using the BLASTX algorithm. When the E-value cutoff was set at 10-5, a total of 25965 contig sequences were annotated, which accounted for 57.5% of all contig sequences.
In principle, the higher the number of reads assembled in a contig, the higher the number of mRNA molecules encoding that gene in a given tissue sample. The read numbers of more than 700 are shown in Additional file 1. Among of the 58 contigs, 4 were involved in cell wall remodeling during plant growth and development, 10 were involved in the adaptive stress response and regulation of expression of many stress-responsive genes, and 2 were involved in polyamine synthesis. The others were involved in protein translation, synthesis, and degradation. Plant cell walls play many critical roles during plant growth, including regulation of cell differentiation, intercellular adhesion and communication, and defense against invasions by pests and pathogens -, knowing which genes are involved in the formation and remodeling of plant cell walls is of great importance. Cell wall remodeling-related genes detected in this study should be related to the wound of cuttings. Stress-related genes could enhance survival of various cell types against stress. Polyamines are key regulators in cell growth and differentiation, the detection of polyamine synthesis-related genes should be related to the different development stages of adventitious root. The detailed discussion of stress-related genes and polyamine synthesis-related genes could be found in the following sections.
We used the Gene Ontology (GO) classification system to determine the possible functions of tagged genes. The GO ontology analysis showed that the distributions of gene functions for cDNA sequences from 25-5 and 23-12 on 14 days after cutting (DAC) and 25 DAC were similar (Figures 3 and 4). The functions of the identified genes covered various biological processes. Among 106784 unigenes (45137 contigs +61647 singletons), 13316 unigenes were classified into the biological process (BP) category, 5590 unigenes were classified into cellular component (CC), and 10885 unigenes were classified into molecular function (MF). The BP, CC, and MF categories in the GO classification system contained 34 subgroups (Figures 3 and 4), which covered the whole profile of gene expression in hybrid larch. The numbers of genes involved in some categories were significant difference between 14 DAC and 25 DAC (Figures 3 and 4). For 25-5 (Figures 3), the “catalytic activity” and “binding” genes were detected increased at 25 DAC than 14 DAC (p value were 0.0307 and 0.0484, Additional file 2), which was almost no change for 23-12 (Figures 4). “catalytic activity” and “binding” genes, involved in many biological processes, were detected increased at adventitious root meristem formation stage (25DAC), which would play an important effect on root formation. For 23-12 (Figures 4), the “reproductive process” genes were detected significant decrease at 25 DAC (p value was 0.0335), which was almost no change for 25-5 (Figures 3). Although the “immune system process” genes have no significant changes between 14 DAC and 25 DAC for 25-5 and 23-12 separately, the genes increased 29.8% at 25 DAC for 25-5, decreased 18.6% at 25 DAC for 23-12, which would be better for 25-5 to form adventitious roots. Adventitious root formation is a complex process regulated by many factors, the “signaling process” genes were significant increase at 25 DAC for 25-5 (p value was 0.0241), but almost no change for 23-12, which indicated signaling process involved in regulation of adventitious root meristem formation (25DAC). The detailed information of the number of detected genes between 14 DAC and 25 DAC and the p values of each GO category were shown in Additional file 2.
Differentially expressed proteins
The profiles of proteins extracted from cutting tissue on 14 DAC, 25 DAC, and 35 DAC were investigated using a proteomics approach. We detected 86 differentially expressed proteins (fold changes ≥2, P ≤0.05) by 2D-DIGE. These proteins were subsequently analyzed using DeCyder™ software (Figure 5, spots 1–75). Approximately one-third of the differentially expressed proteins (32.4%) were present at 25 DAC (the initial cell formation stage in adventitious roots). We successfully identified 75 proteins using MALDI-TOF/TOF-MS. Detailed information for these proteins are shown in Additional file 3. For most of the proteins, the theoretical molecular mass (MM) and calculated isoelectric point (pI) was consistent with experimental data, as judged from the location of spots on 2D-DIGE gels (Figure 5; Additional file 3). However, there were some exceptions, such as spots 7, 8, 9, 15, 29, 41, and 66, which had theoretical MMs greater than the corresponding identified proteins, and spots 32, 37, 48, 49, 54, and 71, which had theoretical MMs smaller than the corresponding identified proteins (see Additional file 3). Several proteins, had the same functional annotation, showed in several different spots (spots 10, 11, 25, 26, and 27) on the gel (Figure 5). This phenomenon has also been reported in previous studies . There are several possible explanations for this, including post-translational modification, protein degradation, different encoding genes, and protein translation from alternatively spliced mRNAs. Another possibility is that these spots represented protein homologs unique to larch.
The 75 identified proteins could be divided into eight categories according to their functional annotation (see Additional file 3). Of these proteins, 14.7% were involved in carbohydrate and energy metabolism, 6.7% in protein translation and degradation, 34.7% in cell defense, 9.3% in hormone-related activities, 5.3% in signal transduction, 5.3% in RNA transcription and processing, and a few were involved in other metabolic pathways. Proteins of unknown function made up 24% of the 75 identified proteins. Figure 6A,B shows the abundance of some of the differentially expressed proteins at different stages of adventitious root development.
One notable result was that the expression patterns of S-adenosylmethionine synthetase 2 (SAMS2, spot 21, Figure 6A,B) and S-adenosyl methionine decarboxylase (SAMDC, spot 22), which are involved in polyamine synthesis, and a late-embryogenesis-abundant protein (spot 33) and an intracellular pathogenesis-related protein (spot 45), which are involved in stress responses, were the same as those detected for their encoding genes in the transcriptomic analysis. The detailed discussions of the effects of these genes on adventitious root development are showed in “elucidation of important signaling pathways” part.
Elucidation of important signaling pathways
Polyamine synthesis pathway
Polyamines are small polycationic molecules; the main polyamines in living organisms are putrescine (Put), spermidine (Spd), and spermine (Spm) . They are found ubiquitously in living organisms and function in a wide variety of biological processes. Polyamine biosynthesis involves a series of steps catalyzed by specific enzymes ,, including SAMDC and SAMS. In this study, the contigs UN_lyscdhit99_32363 and UN_lyscdhit99_25323 were annotated as SAMS2 and SAMDC proenzymes, respectively, in the transcriptome data (Table 1). The read numbers assembled in these two genes in 25-5 (strong rooting capacity) were higher than those in 23-12 (weak rooting capacity) at the root cell dedifferentiation and division stage (14 DAC) and at the root initial cell formation stage (25 DAC). In the proteomic profile, two proteins were identified as SAMS2 (spot 21) and SAMDC (spot 22); importantly, these two proteins were expressed at higher levels in 25-5 than in 23-12 on both 14 DAC and 25 DAC (Table 1). The SAMDC protein was also up-regulated at 35 DAC (see Additional file 3). This result that polyamines play a role in regulating root development is consistent with the results of other studies. For example, Tisi et al.  showed that polyamine catabolism strongly affected root development and xylem differentiation in maize (Zea mays). Niemi et al.  showed that polyamines accelerated adventitious root formation and increased subsequent root growth of Scots pine (Pinus sylvestris L.) hypocotyl cuttings in vitro. Tang and Newton  demonstrated that polyamines promoted root elongation and growth in regenerated Virginia pine plantlets. In other studies, polyamines were shown to improve the rooting rate , regulate vascular development , and enhance tolerance of the root to salinity-alkalinity stress . In this study, the expression of SAMDC and SAMS2 at the transcriptional and translational levels demonstrated that polyamines play an important role in regulating adventitious root development in hybrid larch, possibly via increasing root cell division and enhancing tolerance of the roots to stresses.
Stress response pathway
Cuttings excised from plants are subjected to several stresses, such as wounding and pathogen attack. The responses of cuttings to stresses are complex and involve a number of metabolic changes to protect cells and macromolecules. In this study, Contigs of UN_lyscdhit99_17285 and UN_lyscdhit99_23395 were annotated as an embryo-abundant protein and a putative intracellular pathogenesis-related protein, respectively (Table 1). Both of these proteins are related to various stress responses. We identified corresponding proteins in the proteomic profile; the embryo-abundant protein corresponded to spot 33, and the intracellular pathogenesis-related protein corresponded to spot 45.
The late embryogenesis abundant (LEA) protein family is a large family that includes proteins that accumulate at late stages of seed development or in vegetative tissues in response to some stresses. Many studies have shown that LEA proteins have important effects on salt tolerance in plants. Aghaei et al.  showed that a LEA protein was up-regulated in the hypocotyl and root of soybean treated with 100 mM NaCl. Zhang et al.  identified a putative salt-tolerance gene LEA1 in Thellungiella salsuginea, and confirmed that TsLEA1 conferred salt tolerance when expressed in yeast and transgenic Arabidopsis. Bai et al.  demonstrated that transgenic tobacco expressing the LEA3-1 protein encoded by a gene from Medicago sativa showed enhanced salt tolerance. Park et al.  showed that LEA14 might positively regulate the response to dehydration by enhancing root cell lignification in sweet potato. Overexpression of SmLEA showed faster root elongation and a higher salt and drought tolerance in Salvia miltiorrhiza . In this study, a LEA protein was annotated in the transcriptome data. The read number assembled in this gene was higher in clone 25-5 (strong rooting capacity) than in clone 23-12 (weak rooting capacity) at the root cell dedifferentiation and division stage (14 DAC). Consistent with this, the corresponding protein was expressed at higher levels in 25-5 than in 23-12 at 14 DAC (Table 1). However, the read number assembled in LEA was lower in 25-5 than in 23-12 at the root initial cell formation stage (25 DAC), and there was no significant difference in the read number assembled in LEA between 25-5 and 23-12 at 35 DAC (see Additional file 3). Together, these findings conjectured that this LEA might positively regulate the response to various stresses by enhancing lignification of cuttings, mainly at the early stage of adventitious root development.
Pathogenesis-related (PR) proteins are plant proteins induced by abiotic and biotic stresses . They accumulate around damaged cell walls to protect plants against infection by fungi, bacteria, or viruses . Many studies have shown that PR proteins have important effects on root development. Bantignies et al.  showed that a PR-10-like protein was constitutively expressed at all stages of root development in roots of white lupin. Borghi et al.  showed that PR genes affected potassium homeostasis in Arabidopsis thaliana. Takeuchi et al.  identified a root-specific PR protein induced by drought and salt treatments in rice. Immunohistochemical analyses showed that this protein strongly accumulated in cortex cells around the vascular system of roots. Koehler et al.  identified PR proteins associated with cold tolerance in strawberry. The abundance of PRprotein1.2 significantly increased in wheat seminal roots after 7 days of waterlogging . In this study, the expression of a PR protein at transcriptional and translational levels implied that it has a role in adventitious root development. It was significantly up-regulated in 25-5 (strong rooting capacity) at the adventitious root cell dedifferentiation and division stage (14 DAC) (Table 1). Therefore, it may play a crucial role in protecting the cuttings from damage caused by bacteria and fungi, mainly at the early stage after cutting, thus allowing normal adventitious root development.
In the present study, two clones of hybrid larch, different on rooting capacity, were analyzed from RNA and protein level. A transcriptome database of adventitious root development was generated included a total of 957832 raw reads, which would be useful to future molecular studies of larch. The protein profile of adventitious root development was also analyzed, 75 proteins were identified, which enriched the protein database. A joint analysis of transcriptome and proteome showed.that genes related to polyamine synthesis and the stress response might play an important role on adventitious root development. This study provides fundamental and important information for subsequent studies on the molecular mechanism of adventitious roots development. Furthermore, we demonstrated for the first time the combined use of two important technologies as a powerful approach for studying non-model, but otherwise important, larch species.
In this study, two clones of L. kaempferi × L. olgensis, 25-5 (strong rooting capacity) and 23-12 (weak rooting capacity), and two early stages of adventitious root development, the root cell dedifferentiation and division stage (14 DAC) and the root initial cell formation stage (25 DAC) , are selected to study the molecular basis of developmental regulation of adventitious root formation. All plant material used in this study was propagated at Dagujia forest farm in Liaoning province in northeastern China (longitude of 42°2'48''N, latitude of 124°47'48''E; altitude of 225–394 m; and annual average temperature of 5.6°C) . A completely randomized block design with three replicates was used. Four to five cuttings were randomly taken from each replicate, and cuttings from each replicate were separated according to clones and development stages. A 5-mm-thick portion from the bottom of each cutting was taken for RNA and protein extraction. Images of the stem cuttings at different stages of adventitious root development are shown in Figure 7A. The 5-mm-thick portion used for analyses is illustrated in Figure 7B.
Our experimental research complied with institutional, national, or international guidelines. Field studies were conducted in accordance with local legislation, and with the permissions of forestry administrative department. We did not use any endangered species and complied with the Convention on the Trade in Endangered Species of Wild Fauna and Flora. Voucher specimens (LH 2013001, LH 2013002, LH 2013003, LH 2013004), were identified by Prof. Xiaoshan Wang of the research institute of forestry, Bejing, and have been deposited in the herbarium of Chinese Academy of Forestry.
RNA extraction and cDNA library construction
Total RNA was prepared by the method of Sánchez et al. . The concentration and quality of total RNA were determined using a NanoDrop 1000 spectrophotometer (Thermo Scientific, Wilmington, DE, USA) and checked on 1% agarose gels before proceeding (see Additional file 4). Approximately 1 μg total RNA was converted into cDNA using a SMART cDNA synthesis kit (Clontech, Palo Alto, CA, USA), optimizing the conditions to obtain a large quantity of clean cDNA in a small volume. For second-strand synthesis, PCR was carried out using a small aliquot (1/10th volume) of the primary template and Advantage 2 Polymerase Mix (Clontech). The thermal cycling program was as follows: initial denaturation at 95°C for 60 s, followed by 18 cycles of denaturation at 95°C for 15 s and annealing at 65°C for 30 s, followed by extension at 68°C for 3 min. For each sample and all PCR reactions, the PCR products were purified using a PureLink PCRTM purification kit (Invitrogen, Carlsbad, CA, USA). Double-stranded cDNA was quantified with a spectrophotometer (NanoDrop 1000, Thermo Scientific). The products were checked on a 2% agarose gel to verify cDNA quality and fragment length.
sequencing and assembly
Digested cDNA was recovered with a QIAquick PCR Purification kit (Qiagen, Hilden, Germany). Approximately 7 μg ds cDNA was sheared via nebulization into small fragments, and sequenced using the GS FLX Titanium kit. Raw unprocessed EST sequences generated from this study have been submitted to the Short Read Archive (SRA) division of Genbank. The 454 SFF file containing raw sequences and sequence quality information can be accessed through the SRA web site under the accession number SRP015266. Using the GS FLX pyrosequencing software, we selected high-quality sequences (>99.5% accuracy on single base reads) for further processing and assembly. Adapter trimming and poly (A/T) and short sequence (<50 bp) removal were performed by in-house Perl scripts to obtain clean ESTs. We used Newbler software (provided with the Roche GS FLX sequencer) for sequence assembly. The quality score threshold was set at 40.
Functional annotation and classification
The assembled unique transcripts were compared with the sequences in the non-redundant database of GenBank using the BLASTN algorithm to find and remove ribosomal RNA sequences . The remaining sequences that putatively encoded proteins were searched against the NCBI non-redundant protein (Nr) database (http://www.ncbi.nlm.nih.gov/) using the BLASTX algorithm, and against the UniProt protein database (http://www.uniprot.org/help/uniprotkb) and the KEGG protein database (http://www.genome.jp/kegg/). A typical cut-off value of E <1.0-5 was used. The GO  system was used for functional classification of sequences. GO provides a structured and controlled vocabulary to describe gene products according to three categories: molecular function, biological process, and cellular component (http://geneontology.org/).
Gene expression analysis
Normalization between the libraries
We did normalization between the libraries, including the comparison of the number of unigenes detected in each library of GO classification, and the total number of reads in each library when found the difference between two libraries. The detailed information of normalization was shown in Additional file 5.
Protein extraction and preparation
For each sample, 2 g tissue was ground into a fine powder in liquid nitrogen for protein extraction. The lyophilized powder was homogenized in lysis buffer (7 M urea, 2 M thiourea, 4% 3-[(3-cholamidopropyl) dimethylammonio]-1-propanesulfonate (CHAPS), 10 mM Tris) for 30 min at room temperature with repeated shaking. The undissolved powder was removed from the homogenate by centrifugation at 40 000 × g for 60 min at 4°C. The supernatant was stored in aliquots at –80°C. Protein concentration was determined with a Bradford assay kit (BioRad, Hercules, CA, USA) using albumin diluted in lysis buffer as the internal standard.
Protein labeling with CyDye DIGE fluor
Sample lysates were labeled with Cy2, Cy3, and Cy5 following the protocols described in the Ettan DIGE User Manual (18-1164-40 Edition AA, GE Healthcare, Buckinghamshire, UK). Typically, 50 μg lysate was labeled with 400 pmol Cy3 or Cy5, while the same volume of a pooled standard containing equal quantities of all samples was labeled with Cy2. Labeling reactions were carried out in the dark on ice for 30 min before quenching with 1 mL 10 mM lysine for 10 min on ice. These labeled samples were then combined for 2D-DIGE analysis.
2D-DIGE and image analysis
We conducted 2D-DIGE as follows: The IPG strips (24 cm, pH 3–10, and NL) were rehydrated with labeled samples in the dark overnight with rehydration buffer (7 M urea, 4% w/v CHAPS, 20 mM dithiothreitol (DTT), and 1% v/v IPG buffer with a trace amount of bromophenol blue). First-dimension IEF was performed using an Ettan IPGphor System (GE Healthcare) for a total of 67 kV h at 20°C. The strips were then treated with a two-step reduction and an alkylation step prior to the second-dimension SDS-PAGE. After equilibration with a solution containing 6 M urea, 30% glycerol, 2% sodium dodecyl sulfate (SDS), 50 mM Tris-Cl, pH 8.8, and 0.5% w/v DTT, the strips were treated with the same solution except that it contained 4.5% w/v iodoacetamide instead of DTT. The strips were overlaid onto 12% polyacrylamide gels (20 × 24 cm), immobilized to a low-fluorescence glass plate, and electrophoresed for 12–18 h at 30 mA per gel using an Ettan DALT Twelve System (GE Healthcare). The Cy2, Cy3, and Cy5-labeled images were acquired using a Typhoon 9410 scanner (GE Healthcare) at excitation/emission wavelengths of 488/520 nm, 532/580 nm, and 633/670 nm, respectively. The DIGE images were analyzed by DeCyder 6.5 following the Ettan DIGE User Manual (GE Healthcare). Spot detection was performed using the DIA (differential in-gel analysis) module with the estimated number of spots set to 2500. Only spots present in all three replicate gels and qualitatively consistent in size and shape were considered. After removing artifact spots by manual editing, the DIGE images were further analyzed using the DeCyder BVA (biological variation analysis) module. A 0.05 significance level was used to define differences between groups when analyzing parallel spots. One-way ANOVA and Student-Newman-Keuls tests were conducted using the SAS software package version 8.2 (SAS Institute).
In-gel tryptic digestion and MALDI-TOF/TOF-MS analysis
Separate preparative gels were run to obtain sufficient amounts of protein for MS analysis. These gels were fixed and stained with colloidal Coomassie brilliant blue (CBB). Proteins of interest, as defined by the 2D- DIGE/DeCyder analysis, were excised from the CBB-stained gels for a modified in-gel tryptic digestion procedure. Gel pieces were first discolored in 50% acetonitrile and 25 mM ammonium bicarbonate, and then subjected to reduction and alkylation in 10 mM DTT and 55 mM iodoacetic acid, respectively. Following vacuum drying, the gel pieces were incubated with sequencing-grade modified trypsin (Promega, Madison, WI, USA) at a final concentration of 0.01 mg/mL in 25 mM ammonium bicarbonate for 16 h at 37°C. Supernatants were collected, vacuum dried, and then redissolved in 50% acetonitrile and 0.1% trifluoroacetic acid (TFA) for MS analysis using an ABI 4800 Proteomics Analyzer MALDI-TOF/TOF-MS. The TOF spectra were recorded in the positive ion reflector mode with a mass range from 700 to 4000 Da, and 10 of the strongest peaks in each sample were chosen for MS/MS analysis. The spectra were corrected by an external standard method using trypsin-treated myoglobin peptides. The MS/MS results were searched using GPS (Applied Biosystems, Foster City, CA, USA) - MASCOT (Matrix Science, London, UK) with the following criteria: NCBInr database; species restriction, green plants; MS tolerance was set at 6100 ppm and MS/MS at 60.6 Da; at most one missed cleavage site; fixed modification was carbamidomethyl (Cys) and variable modification was oxidation (Met); and cleavage by trypsin was on the C-terminal side of Lys and Arg unless the next residue was Pro. If peptides matched to multiple members of a protein family, or if a protein appeared under different names and accession numbers, the entry with the highest score was selected. In addition, the theoretical molecular weights and pI of the identified proteins were calculated using the Peptide Mass program (http://web.expasy.org/peptide_mass/).
Availability of supporting data
The raw sequences data supporting the results of this article are available in the Short Read Archive (SRA) (accession number SRP015266), http://www.ncbi.nlm.nih.gov/sra.
Pâques LE: Roles of European and Japanese larch in the genetic control of growth, architecture and wood quality traits in interspecific hybrids (Larix × eurolepis Henry). Ann For Sci. 2004, 61: 25-33. 10.1051/forest:2003081.
Han H, Sun XM, Xie YH, Feng J, Zhang SG: Anatomical and physiological effects of phytohormones on adventitious roots development in Larix kaempferi × L. olgensis . Silvae Genetica. 2013, 62 (3): 96-104.
Leyser O: Auxin distribution and plant pattern formation: how many angels can dance on the point of PIN. Cell. 2005, 121: 819-822. 10.1016/j.cell.2005.06.005.
Liu H, Wang S, Yu X, Yu J, He X, Zhang S, Shou H, Wu P: ARL1, a LOB-domain protein required for adventitious root formation in rice. Plant J. 2005, 43: 47-56. 10.1111/j.1365-313X.2005.02434.x.
Cheng Y, Dai X, Zhao Y: Auxin biosynthesis by the YUCCA flavin monooxygenases control the formation of floral organs and vascular tissues in Arabidopsis . Genes Dev. 2006, 20: 1790-1799. 10.1101/gad.1415106.
Han H, Zhang SG, Sun XM: A review on the molecular mechanism of plants rooting modulated by auxin. Afr J Biotechnol. 2009, 8 (3): 348-353.
Sorin C, Bussell JD, Camus I, Ljung K, Kowalczyk M: Auxin and light control of adventitious rooting in Arabidopsis require ARGONAUTE1. Plant Cell. 2005, 17: 1343-1359. 10.1105/tpc.105.031625.
Brinker M, van Zyl L, Liu WB, Craig D, Sederoff RR: Microarray analyses of gene expression during adventitious root development in Pinus contorta . Plant Physiol. 2004, 135: 1526-1539. 10.1104/pp.103.032235.
Butler E, Gallagher TF: Gene expression during adventitious root formation in apple. Symp Soc Exp Biol. 1998, 51: 79-84.
Sorin C, Negroni L, Balliau T, Corti H, Jacquemot MP: Proteomic analysis of different mutant genotypes of Arabidopsis led to the identification of 11 proteins correlating with adventitious root development. Plant Physiol. 2006, 140: 349-364. 10.1104/pp.105.067868.
Chen XD, Qu YL, Sheng LH, Liu JC, Huang H, Xu L: A simple method suitable to study de novo root organogenesis. Front Plant Sci. 2014, 5: 208-10.3389/fpls.2014.00208.
Ahkami A, Scholz U, Steuernagel B, Strickert M, Haensch KT, Druege U, Reinhardt D, Nouri E, von Wirén N, Franken P, Hajirezaei MR: Comprehensive transcriptome analysis unravels the existence of crucial genes regulating primary metabolism during adventitious root formation in Petunia hybrida. PLoS One. 2014, 9 (6): e100997-10.1371/journal.pone.0100997.
Subramaniyam S, Mathiyalagan R, Natarajan S, Kim YJ, Jang MG, Park JH, Yang DC: Transcript expression profiling for adventitious roots of Panax ginseng Meyer. Gene. 2014, 546 (1): 89-96. 10.1016/j.gene.2014.05.024.
Wu D, Austin RS, Zhou SJ, Brown D: The root transcriptome for North American ginseng assembled and profiled across seasonal development. BMC Genomics. 2013, 14: 564-10.1186/1471-2164-14-564.
Firon N, LaBonte D, Villordon A, Kfir Y, Solis J, Lapis E, Perlman TS, Doron-Faigenboim A, Hetzroni A, Althan L, Nadir LA: Transcriptional profiling of sweetpotato (Ipomoea batatas) roots indicates down-regulation of lignin biosynthesis and up-regulation of starch biosynthesis at an early stage of storage root formation. BMC Genomics. 2013, 14: 460-10.1186/1471-2164-14-460.
Alagna F, D’Agostino N, Torchia L, Servili M, Rao R, Pietrella M, Giuliano G, Chiusano LM, Baldoni L, Perrotta G: Comparative 454 pyrosequencing of transcripts from two olive genotypes during fruit development. BMC Genomics. 2009, 10: 399-10.1186/1471-2164-10-399.
Yang HY, Hu LW, Hurek T, Hurek BR: Global characterization of the root transcriptome of a wild species of rice, Oryza longistaminata, bydeep sequencing. BMC Genomics. 2010, 11: 705-10.1186/1471-2164-11-705.
Lin X, Zhang J, Li Y, Luo H, Wu Q, Sun C, Song J, Li X, Wei J, Lu A, Qian Z, Khan IA, Chen S: Functional genomics of a living fossil tree Ginkgo based on next generation sequencing technology. Physiol Plant. 2011.
Wu Q, Sun C, Luo H, Li Y, Niu Y, Sun Y, Lu A, Chen S: Transcriptome analysis of Taxus cuspidata needles based on 454 pyrosequencing. Planta Med. 2011, 77 (4): 394-400. 10.1055/s-0030-1250331.
Li J, McConkey BJ, Cheng Z, Guo S, Glick BR: Identification of plant growth-promoting bacteria-responsive proteins in cucumber roots under hypoxic stress using a proteomic approach. J Proteomics. 2013, 84: 119-131. 10.1016/j.jprot.2013.03.011.
Yan H, Li K, Ding H, Liao C, Li X, Yuan L, Li C: Root morphological and proteomic responses to growth restriction in maize plants supplied with sufficient N. J Plant Physiol. 2011, 168 (10): 1067-1075. 10.1016/j.jplph.2010.12.018.
Wang F, Bai MY, Deng Z, Oses-Prieto JA, Burlingame AL, Lu T, Chong K, Wang ZY: Proteomic study identifies proteins involved in brassinosteroid regulation of rice growth. J Integr Plant Biol. 2010, 52 (12): 1075-1085. 10.1111/j.1744-7909.2010.00992.x.
Yang L, Ma C, Wang L, Chen S, Li H: Salt stress induced proteome and transcriptome changes in sugar beet monosomic addition line M14. J Plant Physiol. 2012, 169 (9): 839-850. 10.1016/j.jplph.2012.01.023.
Barakat A, DiLoreto DS, Zhang Y, Smith C, Baier K, Powell WA, Wheeler N, Sederoff R, Carlson JE: Comparison of the transcriptomes of American chestnut (Castanea dentata) and Chinese chestnut (Castanea mollissima) in response to the chestnut blight infection. BMC Plant Biol. 2009, 9: 51-10.1186/1471-2229-9-51.
Vega-Arreguín JC, Ibarra-Laclette E, Jiménez-Moraila B, Martínez O, Vielle-Calzada JP, Herrera-Estrella L, Herrera-Estrella A: Deep sampling of the Palomero maize transcriptome by a high throughput strategy of pyrosequencing. BMC Genomics. 2009, 10: 299-10.1186/1471-2164-10-299.
Zhang Y, Zhang SG, Han SY, Li XM, Qi LW: Transcriptome profiling and in silico analysis of somatic embryos in Japanese larch (Larix leptolepis). Plant Cell Rep. 2012, 31: 1637-1657. 10.1007/s00299-012-1277-1.
Wang W, Wang YJ, Zhang Q, Qi Y, Guo DJ: Global characterization of Artemisia annua glandular trichome transcriptome using 454 pyrosequencing. BMC Genomics. 2009, 10: 465-10.1186/1471-2164-10-465.
Sun C, Li Y, Wu Q, Luo HM, Sun YZ, Song JY, Lui EMK, Chen SL: De novo sequencing and analysis of the American ginseng root transcriptome using a GS FLX Titanium platform to discover putative genes involved in ginsenoside biosynthesis. BMC Genomics. 2010, 11: 262-10.1186/1471-2164-11-262.
Blanca J, Cañizares J, Roig C, Ziarsolo P, Nuez F, Picó B: Transcriptome characterization and high throughput SSRs and SNPs discovery in Cucurbita pepo (Cucurbitaceae). BMC Genomics. 2011, 12: 104-10.1186/1471-2164-12-104.
Jones DA, Takemoto D: Plant innate immunity - direct and indirect recognition of general and specific pathoge-associated molecules. Curr Opin Immunol. 2004, 16 (1): 48-62. 10.1016/j.coi.2003.11.016.
Scheible WR, Pauly M: Glycosyltransferases and cell wall biosynthesis: novel players and insights. Curr Opin Plant Biol. 2004, 7 (3): 285-295. 10.1016/j.pbi.2004.03.006.
Vorwerk S, Somerville S, Somerville C: The role of plant cell wall polysaccharide composition in disease resistance. Trends Plant Sci. 2004, 9 (4): 203-209. 10.1016/j.tplants.2004.02.005.
Wang S, Yin YB, Ma Q, Tang XJ, Hao DY, Xu Y: Genome-scale identification of cell-wall related genes in Arabidopsis based on co-expression network analysis. BMC Plant Biol. 2012, 12: 138-10.1186/1471-2229-12-138.
Ndimba BK, Chivasa S, Simon WJ, Slabas AR: Identification of Arabidopsis salt and osmotic stress responsive proteins using two-dimensional difference gel electrophoresis and mass spectrometry. Proteomics. 2005, 5 (16): 4185-4196. 10.1002/pmic.200401282.
Liu JH, Kitashiba H, Wang J, Ban Y, Moriguchi T: Polyamines and their ability to provide environmental stress tolerance to plants. Plant Biotechnol. 2007, 24: 117-126. 10.5511/plantbiotechnology.24.117.
Bagni N, Tassoni A: Biosynthesis, oxidation and conjugation of aliphatic polyamines in higher plants. Amino Acids. 2001, 20: 301-317. 10.1007/s007260170046.
Alcázar R, Altabella T, Marco F, Bortolotti C, Reymond M, Koncz C, Carrasco P, Tiburcio AF: Polyamines: molecules with regulatory functions in plant abiotic stress tolerance. Planta. 2010, 231: 1237-1249. 10.1007/s00425-010-1130-0.
Tisi A, Federico R, Moreno S, Lucretti S, Moschou NP, Roubelakis-Angelakis AK, Angelini R, Cona A: Perturbation of polyamine catabolism can strongly affect root development and xylem differentiation. Plant Physiol. 2011, 157: 200-215. 10.1104/pp.111.173153.
Niemi K, Häggman H, Sarjala T: Effects of exogenous diamines on the interaction between ectomycorrhizal fungi and adventitious root formation in Scots pine in vitro. Tree Physiol. 2002, 22 (6): 373-381. 10.1093/treephys/22.6.373.
Tang W, Newton RJ: Polyamines promote root elongation and growth by increasing root cell division in regenerated Virginia pine (Pinus virginiana Mill.) plantlets. Plant Cell Rep. 2005, 24 (10): 581-589. 10.1007/s00299-005-0021-5.
Le Guen-Le SF, Hourmant A: Stimulation of putrescine biosynthesis via the ornithine decarboxylase pathway by gibberellic acid in the in vitro rooting of globe artichoke (Cynara scolymus). Plant Growth Regul. 2001, 35: 277-284. 10.1023/A:1014494013416.
Vera-Sirera F, Minguet EG, Singh SK, Ljung K, Tuominen H, Blázquez MA, Carbonell J: Role of polyamines in plant vascular development. Plant Physiol Biochem. 2010, 48 (7): 534-539. 10.1016/j.plaphy.2010.01.011.
Hu X, Zhang Y, Shi Y, Zhang Z, Zou Z, Zhang H, Zhao J: Effect of exogenous spermidine on polyamine content and metabolism in tomato exposed to salinity-alkalinity mixed stress. Plant Physiol Biochem. 2012, 57C: 200-209. 10.1016/j.plaphy.2012.05.015.
Aghaei K, Ehsanpour AA, Shah AH, Komatsu S: Proteome analysis of soybean hypocotyl and root under salt stress. Amino Acids. 2009, 36 (1): 91-98. 10.1007/s00726-008-0036-7.
Zhang Y, Li Y, Lai J, Zhang H, Liu Y, Liang L, Xie Q: Ectopic expression of a LEA protein gene TsLEA1 from Thellungiella salsuginea confers salt-tolerance in yeast and Arabidopsis . Mol Biol Rep. 2012, 39 (4): 4627-4633. 10.1007/s11033-011-1254-8.
Bai YQ, Yang QC, Kang JM, Sun Y, Gruber M, Chao YH: Isolation and functional characterization of a Medicago sativa L. gene, MsLEA3-1 . Mol Biol Rep. 2012, 39: 2883-2892. 10.1007/s11033-011-1048-z.
Park SC, Kim YH, Jeong JC, Kim CY, Lee HS, Bang JW, Kwak SS: Sweetpotato late embryogenesis abundant 14 (IbLEA14) gene influences lignification and increases osmotic- and salt stress-tolerance of transgenic calli. Planta. 2011, 233 (3): 621-634. 10.1007/s00425-010-1326-3.
Wu Y, Liu C, Kuang J, Ge Q, Zhang Y, Wang Z: Overexpression of SmLEA enhances salt and drought tolerance in Escherichia coli and Salvia miltiorrhiza. Protoplasma. 2014, 251 (5): 1191-1199. 10.1007/s00709-014-0626-z.
Van Loon Rep LC, Pieterse CMJ: Significance of inducible defense-related proteins in infected plants. Annu Rev Phytopathol. 2006, 44: 135-162. 10.1146/annurev.phyto.44.070505.143425.
Liu JJ, Ekramoddoullah AKM: The family 10 of plant pathogenesis-related proteins: their structure, regulation, and function in response to biotic and abiotic stresses. Physiol Mol Plant Pathol. 2006, 68: 3-13. 10.1016/j.pmpp.2006.06.004.
Bantignies B, Séguin J, Muzac I, Dédaldéchamp F, Gulick P, Ibrahim R: Direct evidence for ribonucleolytic activity of a PR-10-like protein from white lupin roots. Plant Mol Biol. 2000, 42 (6): 871-881. 10.1023/A:1006475303115.
Borghi M, Rus A, Salt DE: Loss-of-function of constitutive expresser of pathogenesis related genes5 affects potassium homeostasis in Arabidopsis thaliana . PLoS One. 2011, 6 (10): e26360-10.1371/journal.pone.0026360. doi:10.1371/journal.pone.0026360,
Takeuchi K, Gyohda A, Tominaga M, Kawakatsu M, Hatakeyama A, Ishii N, Shimaya K, Nishimura T, Riemann M, Nick P, Hashimoto M, Komano T, Endo A, Okamoto T, Jikumaru Y, Kamiya Y, Terakawa T, Koshiba T: RSOsPR10 expression in response to environmental stresses is regulated antagonistically by jasmonate/ethylene and salicylic acid signaling pathways in rice roots. Plant Cell Physiol. 2011, 52 (9): 1686-1696. 10.1093/pcp/pcr105.
Koehler G, Wilson RC, Goodpaster JV, Sønsteby A, Lai XY, Witzmann FA, You JS, Rohloff J, Randall SK, Alsheikh M: Proteomic study of low temperature responses in strawberry cultivars (Fragaria x ananassa Duchesne) that differ in cold tolerance. Plant Physiol. 2012.
Haque ME, Abe F, Mori M, Oyanagi A, Komatsu S, Kawaguchi K: Characterization of a wheat pathogenesis-related protein, TaBWPR-1.2, in seminal roots in response to waterlogging stress. J Plant Physiol. 2014, 171 (8): 602-609. 10.1016/j.jplph.2013.12.003.
Sánchez C, Vielba JM, Ferro E, Covelo G, Sole A, Abarca D, de Mier BS, Díaz-Sala C: Two SCARECROW-LIKE genes are induced in response to exogenous auxin in rooting-competent cutting of distantly related forest species. Tree Physiol. 2007, 27: 1459-1470. 10.1093/treephys/27.10.1459.
Altschul SF, Gish W, Miller W, Myers EW, Lipman DJ: Basic local alignment search tool. J Mol Biol. 1990, 215 (3): 403-410. 10.1016/S0022-2836(05)80360-2.
Consortium GO: The Gene Ontology project in 2008. Nucleic Acids Res. 2008, 36: D440-D444. 10.1093/nar/gkm883.
Audic S, Claverie JM: The significance of digital gene expression profiles. Genome Res. 1997, 7 (10): 986-995.
This study was financially supported by the National Science and Technology Support Program (2012BAD01B01) and the National Science Foundation (30972393). The authors thank Dr. Liang Zhang for the service of cDNA library construction and 454 sequencing. The authors are also grateful to Edanz (http://www.liwenbianji.cn/) for their helpful suggestions on improving the manuscript.
The authors declare that they have no competing interests.
HH performed sample collection, RNA extraction, protein extraction, 2D-DIGE, MALDI-TOF/TOF-MS, analyzed the data, and drafted the manuscript; XS supervised HH in bioinformatic work, performed bioinformatic analysis of the data, and helped to write and edit the manuscript; YX helped with MALDI-TOF/TOF-MS analysis; JF helped with plant preparation and sample collection; SZ coordinated the project design, and edited the manuscript. All authors read and approved the final manuscript.
Electronic supplementary material
Additional file 2: The detailed information of the number of detected genes between 14 DAC and 25 DAC in clones 25-5 and 23-12, and the p values of each GO category. (XLS 38 KB)
Additional file 3: Identification of differentially expressed proteins between clones 23-12 and 25-5 of L. kaempferi × L. olgensis using MALDI-TOF/TOF-MS. (PDF 21 KB)
Authors’ original submitted files for images
Below are the links to the authors’ original submitted files for images.
About this article
Cite this article
Han, H., Sun, X., Xie, Y. et al. Transcriptome and proteome profiling of adventitious root development in hybrid larch (Larix kaempferi × Larix olgensis). BMC Plant Biol 14, 305 (2014). https://doi.org/10.1186/s12870-014-0305-4