Next Article in Journal
Identification of Cumin (Cuminum cyminum) MicroRNAs through Deep Sequencing and Their Impact on Plant Secondary Metabolism
Next Article in Special Issue
Genome-Wide Identification and Expression Analysis of CCO Gene Family in Liriodendron chinense
Previous Article in Journal
Current Progress and Future Prospect of Wheat Genetics Research towards an Enhanced Nitrogen Use Efficiency
Previous Article in Special Issue
Uses, Botanical Characteristics, and Phenological Development of Slender Nightshade (Solanum nigrescens Mart. and Gal.)
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Plants
Volume 12
Issue 9
10.3390/plants12091754
plants-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Effects of Different Donor Ages on the Growth of Cutting Seedlings Propagated from Ancient Platycladus orientalis

by
Yao Dong
1,
Wei Guo
2,
Wenfa Xiao
1,
Jianfeng Liu
3,
Zirui Jia
3,
Xiulian Zhao
3,
Zeping Jiang
1 and
Ermei Chang
3,*
1
Key Laboratory of Forest Ecology and Environment of National Forestry and Grassland Administration, Ecology and Nature Conservation Institute, Chinese Academy of Forestry, Beijing 100091, China
2
Taishan Academy of Forestry Sciences, Taian 271000, China
3
State Key Laboratory of Tree Genetics and Breeding, Key Laboratory Tree Breeding and Cultivation of the National Forestry and Grassland Administration, Research Institute of Forestry, Chinese Academy of Forestry, Beijing 100091, China
*
Author to whom correspondence should be addressed.
Plants 2023, 12(9), 1754; https://doi.org/10.3390/plants12091754
Submission received: 22 March 2023 / Revised: 18 April 2023 / Accepted: 23 April 2023 / Published: 25 April 2023
(This article belongs to the Special Issue Plant Science Research in Nanjing Agricultural University: Commemorating the University’s 120th Anniversary)
Browse Figures
Versions Notes

Abstract

:
The effects of tree age on the growth of cutting seedlings propagated from ancient trees have been an important issue in plant breeding and cultivation. In order to understand seedling growth and stress resistance stability, phenotypic measurements, physiological assays, and high-throughput transcriptome sequencing were performed on sown seedlings propagated from 5-year-old donors and cutting seedlings propagated from 5-, 300-, and 700-year-old Platycladus orientalis donors. In this study, the growth of cutting seedlings propagated from ancient trees was significantly slower; the soluble sugar and chlorophyll contents gradually decreased with the increase in the age of donors, and the flavonoid and total phenolic contents of sown seedlings were higher than those of cutting seedlings. Enrichment analysis of differential genes showed that plant hormone signal transduction, the plant–pathogen interaction, and the flavone and flavonol biosynthesis pathways were significantly up-regulated with the increasing age of cutting seedlings propagated from 300- and 700-year-old donors. A total of 104,764 differentially expressed genes were calculated using weighted gene co-expression network analysis, and 8 gene modules were obtained. Further, 10 hub genes in the blue module were identified, which revealed that the expression levels of JAZ, FLS, RPM1/RPS3, CML, and RPS2 increased with the increase in tree age. The results demonstrated that the age of the donors seriously affected the growth of P. orientalis cutting seedlings and that cutting propagation can preserve the resistance of ancient trees. The results of this study provide important insights into the effects of age on asexually propagated seedlings, reveal potential molecular mechanisms, and contribute to an improvement in the level of breeding and conservation of ancient germplasm resources of P. orientalis trees.

1. Introduction

The conservation of germplasm resources of ancient and notable trees has been paid increasing attention [1]. Ancient P. orientalis has the characteristics of wide adaptability, strong biotic and abiotic resistance, great longevity, and so on [2]. However, due to climate and human factors, a large number of ancient P. orientalis trees have gradually withered, and it is particularly important to preserve the genetic materials. Cutting is now widely applied in the propagation of ancient tree species, and the age of the donors significantly affects the growth of the cutting seedlings [3]. Cutting seedlings of Robinia pseudoacacia and Taxus globosa propagated from young donors grew much better than those from mature trees [4,5]. There are differences between the seedlings propagated using different methods. For example, hybrid poplar seedlings propagated from tissue culture were superior to rooted cuttings for leaf area, leaf length, leaf width, and leaf dry mass [6]. Propagation methods also affect the anti-stress capability of trees. The rooted cuttings of Passiflora edulis were more tolerant to salt stress than sown seedlings and grafted seedlings. Thus, cutting seedlings are also better able to cope with environmental stresses and survive under extreme conditions [7]. However, the stability of the resistance phenotype of cutting seedlings from ancient trees has not been studied so far. Furthermore, the effects of donor age on the growth of cutting seedlings propagated from ancient trees should be studied in depth.
The age of donor trees is known to be related to the growth status and stress resistance of the cutting seedlings [8,9]. In previous studies, our group found that the peroxidase and superoxide dismutase activities were higher in ancient P. orientalis than those in 20-year-old trees, while the contents of soluble protein and chlorophyll were the opposite [10]. It has been demonstrated that the contents of photosynthetic pigments (chlorophyll a, chlorophyll b, and total chlorophyll) and soluble sugars decreased with age [11]. Under abiotic stresses, the over-expressed flavone metabolic and biosynthetic pathways can directly eliminate reactive oxygen species (ROS) in plants, thereby improving the drought tolerance of plants [12,13,14]. Phenolics play an essential role in the defense reactions of plants against pathogens, and the total phenolic content in plants is influenced by age [15,16], indicating that the accumulation of secondary metabolites such as flavonoids and total phenolics is significantly affected by plant age. Different tree ages can lead to different responses of plants to pathogens. For example, adult Dianthus pavonius is much more resistant to infection secondary to inoculation compared to the resistance of young [17]. Comparison experiments between hybrid poplar at different ages showed that young leaves are more susceptible to leaf spot infection than old leaves [18]. Whether the rooted seedlings can effectively preserve the stress resistance of ancient trees needs further research attention.
The regulatory mechanism for the growth and development of ancient trees may affect the growth of their cutting seedlings. For example, the expression patterns and function in some gene families may be affected by tree age, and these genes may be related to plant growth and development [19,20]. Expression of hormone-related genes, ROS scavenger genes, senescence-related genes, and defense response genes were up-regulated with tree age [21]. It has been reported that some genes belonging to the AP2/ERF and bHLH families affect plant growth and development by regulating the biosynthesis and signal transduction pathways of plant hormones. In addition, the expression of genes related to cell division, cell expansion, and differentiation was significantly reduced in ancient trees, especially miR166 and HD-ZIP III involved in cambium activities. Disease resistance-associated genes and defensive secondary metabolites genes were found to be highly expressed in ancient trees [22]. However, the molecular mechanisms for differences in growth between seedlings propagated from donor trees at different ages have not been clarified.
With the advancement in forest tree breeding and development in forest industries, clonal forestry has played an increasingly important role in the cultivation of planted forests throughout the world. In recent years, with the popularization of molecular biology techniques, more and more studies have begun to focus on the mechanisms for the tree age effect on seedling growth. High-throughput sequencing techniques have made it possible to elucidate the effects of donor age on the growth of cutting seedlings. To make a comprehensive evaluation of the stability and reliability of seedlings propagated from ancient trees at different ages, phenotypic measurements, physiological assays, and high-throughput transcriptome sequencing were performed on sown seedlings propagated from 5-year-old donors and cutting seedlings propagated from 5-, 300-, and 700-year-old donors in this study. The aim of this study was to reveal key genes and metabolic pathways involved in the growth and development of P. orientalis seedlings with increasing donor age, which is of important guiding significance for the conservation and management of germplasm resources of excellent ancient trees using clonal propagation.

2. Results

2.1. Growth and Morphological Characteristics of Sown and Cutting Seedlings of P. orientalis Propagated from Donors at Different Ages

The growth and morphological characteristics of cutting seedlings of P. orientalis propagated from donors at different ages were different (Figure 1). The results showed that the basal stem diameters and plant heights of cutting seedlings propagated from 5-year-old donors were significantly greater than those from 300- and 700-year-old donors (p < 0.05). No significant differences were found between sown and cutting seedlings, both of which were propagated from 5-year-old donors, and neither were there differences between cutting seedlings propagated from 300- and 700-year-old donors, respectively (Figure 2). Therefore, the age of donors significantly affected the growth of the rooted cuttings, and the growth of rooted cuttings propagated from ancient trees was significantly slower.

2.2. Analysis of Soluble Sugar, Chlorophyll, Flavonoid, and Total Phenol Content

The contents of soluble sugar and chlorophyll (Chl. a, Chl. b, and Chl. a + b) gradually decreased with the increase in the age of donors (Figure 3). The contents of Chl. a and Chl. (a + b) in sown seedlings propagated from 5-year-old donors were significantly greater than that of rooted cuttings propagated from 5-year-old donors. The contents of soluble sugar and chlorophyll in cutting seedlings propagated from 5-year-old donors were significantly greater than those from 300- and 700-year-old donors. This showed that the sown seedlings grew more vigorously than the cuttings. The contents of flavonoid and total phenol in sown seedlings propagated from 5-year-old donors were higher than those in cutting seedlings propagated from 5-year-old donors (Figure 3). The contents of flavonoid and total phenol in the cutting seedlings increased with the increasing age of donors.

2.3. Transcriptome Sequencing and Reassembly of P. orientalis Unigenes

A total of 105,530,209 reads were assembled into 104,764 unigenes, and the average length of the assembled unigenes was 1007.31 bp. The N50 value was 1858 bp, and the GC content was 38.75%. Assessment of the assembly integrity was performed using BUSCO scored C: 73.9%. There were 23,652 unigenes of 501–1000 bp, 9559 unigenes of 501–1000 bp, and 20,641 unigenes of >1500 bp. Most unigenes were in the range of 0–500 bp, with an abundance of 49% (Table 1), indicating that the sequencing data met the strict quality requirements for subsequent transcriptome analysis.

2.4. Functional Classification of P. orientalis Unigenes

To assign gene ontology annotation for the unigenes, the BLAST algorithm was used to search seven databases: NR, NT, Swiss-Prot, Pfam, COG, GO, and KEGG. The results showed 26,226 (54.58%) unigenes with homologous sequences in at least one database and 3105 (6.46%) unigenes in all seven databases. Among them, 24,080 unigenes (50.12% of all unigenes) were annotated with a significant BLAST result in the Nr database, 8380 unigenes (17.44%) were annotated in the Nt database, and 18,483 unigenes (38.47%) were annotated in the Swiss-Prot database. Annotation results for unigenes are shown in Table 2.
In this study, a total of 104,764 unigenes were annotated in the GO database. The GO terms were placed in three categories including cellular component (6258 unigenes), biological process (5524 unigenes), and molecular function (5702 unigenes). The GO terms enriched in cellular components were cell part (1766 unigenes) and membrane part (1526 unigenes). The GO terms enriched in the biological process were cellular process and metabolic process, with 1856 and 1541 unigenes, respectively. There were only 109 and 207 unigenes in molecular function with transcription regulator activity and structural molecule activity, respectively (Figure 4).
Unigenes were mainly enriched in nine KEGG pathways. The results indicate that the plant hormone signal transduction pathway has the most differentially expressed genes (DEGs) in B5 vs. C300, while the DEGs in A5 vs. B5 and A5 vs. D700 were mostly enriched in the flavone and flavonol biosynthesis pathways. In addition, the plant–pathogen interaction pathway was significantly enriched with DEGs in A5 vs. D700 (Table 3, Figure 5). Therefore, gene expression may be affected by the propagation method and age of the donor. With increasing tree age, the cuttings of P. orientalis had greater stress resistance and hormone signal transduction capabilities. The biosynthesis of flavone and flavonol in cutting seedlings propagated from 700-year-old donors was enhanced, and secondary metabolism was improved in cutting seedlings propagated from 5-year-old donors, suggesting that they grew vigorously, while the cutting seedlings propagated from 300- and 700-year-old donors had a higher degree of lignification.

2.5. Functional Classification of DEGs

To further investigate the functional classifications of DEGs, we performed a K-means cluster analysis, and the 104,764 unigenes were classified into 20 functional cluster profiles (Figure 6). The genes in profile 19 were mainly involved in the response to bacterium. These genes, which had higher expression levels in cutting seedlings than in sown seedlings, were continuously up-regulated with increasing age, indicating that the cutting seedlings of ancient trees have greater resistance to pathogens. Profile 3 and profile 10 were mainly enriched in response to stimulus and signaling pathways, and they were up-regulated in the cutting seedlings propagated from 700-year-old donors, suggesting that with the influence of donor age, the cutting seedlings propagated from 700-year-old donors grew with modulation of complex defensive response. The genes involved in profile 17, enriched in reproduction, had higher expression levels in sown seedlings than in cutting seedlings, both of which were propagated from 5-year-old donors, and then their expression levels tended to be stable.
Genes involved in profile 0 were enriched in the defense response and were continuously down-regulated with increasing tree age, indicating that the resistance ability of these cutting seedlings was affected by the age of the donors. Genes involved in profiles 1 and 18 exhibited a trend in decreased expression levels except for seedlings propagated from 300-year-old donors. GO functional enrichment demonstrated that these genes were enriched in response to stimulus and cellular component organization. In the three groups of cutting seedlings propagated from donors at different ages, expression levels of the genes in profile 2 were stable, suggesting that the expression of these genes was related to the propagation method.

2.6. Changes in Age-Related Changes in Three Enrichment Pathways

The enriched KEGG pathways included plant hormone signal transduction, the plant–pathogen interaction, and flavone and flavonol biosynthesis, the genes of which were significantly up-regulated in cutting seedlings propagated from ancient donors (Figure 7). GID2/SLY1 (Gibberellin-insensitive dwarf 2/Slender Rice 1), JAR1 (Jasmonate Resistant 1), EFR (Elongation factor Tu receptor), FLS2 (Flavonol synthase 2) and CYP75A (Flavonoid 3′,5′-hydroxylase) were up-regulated with increasing tree age, indicating that with the increase in tree age, cutting seedlings have stronger resistance to stresses. The expression levels of genes simultaneously involved in plant hormone signal transduction and the plant–pathogen interaction in sown seedlings propagated from 5-year-old donors were higher than those in cutting seedlings propagated from 5-, 300-, and 700-year-old donors. These pathways included SNRK2 (Sucrose non-fermenting-1-related protein kinase 2s) and JAR1_4_6, as well as the anti-disease genes RPM1/RPS3 (Resistance to Pseudomonas syringae pv. maculicola 1), RPS2 (Resistance to P. syringae 2) and CML (Calmodulin-like protein), suggesting that cutting seedlings were less resistant to stress than sown seedlings.

2.7. Weighted Gene Co-Expression Network Analysis (WGCNA) Identifying Conserved Differentially Expressed Genes

The 104,764 DEGs were calculated using WGCNA, and 8 modules were obtained. The trend in gene expression levels in the blue module was significantly correlated with increasing tree age (R = −0.585, P = 0.0457; Figure 8a), suggesting that genes in this module were positively associated with age in P. orientalis. The blue module contained 3321 DEGs, most of which were highly expressed in seedlings propagated from 300- and 700-year-old donors. A total of 17 genes related to the plant–pathogen interaction were involved in the MAPK signaling pathway, including 7 members of the plant receptor-like kinase FLS2 gene family. FLS2 is involved in recognizing plant flagellin as a signal of bacterial existence and activating the defense response [23], and its expression levels gradually increased during the growth of Arabidopsis thaliana seedlings [24]. In the blue module, ten genes considered hub genes were involved in key biosynthetic pathways. The expression levels of four plant–pathogen interaction-related genes (EFR/FLS2, RPM1/RPS3, CML, and RPS2) in the module increased with increasing tree age (Figure 8b), indicating that cutting seedlings retain the strong stress resistance of ancient trees. These results indicated that the pathogen defense capability of P. orientalis cutting seedlings was gradually enhanced with the increasing age of the donors.

2.8. Identification of Hub Genes in the Blue Module

Ten genes in the blue module were involved in key biosynthetic pathways, and they were considered hub genes in the module. The expression levels of four plant–pathogen interaction-related genes (EFR/FLS2, RPM1/RPS3, CML, and RPS2) in the module increased with the increase in tree age (Figure 8b), indicating that cutting seedlings retain the strong stress resistance of ancient trees. Among the two genes (CCR (Cinnamoyl-CoA reductase) and F3H (Flavanone 3-hydroxylase)) involved in flavonoid biosynthesis, the expression trend in F3H increased with tree age. F3H plays an important role in plant fiber development and the anthocyanin synthesis pathway [25]. The differential expressions of these genes affect the accumulation of secondary metabolites and leaf color in cutting seedlings. JAZ (Jasmonate ZIM-domain protein) and JAR1_4_6, and K14486/ARF (Auxin response factor) participate in plant hormone signal transduction, among which JAZ and JAR1_4_6 are members of the JAZ family and are highly expressed in ancient tree cutting seedlings. JAZ is a key regulator of the JA signaling pathway, which can regulate the biosynthesis of secondary metabolites such as terpenes, alkaloids, and flavonoids. This suggests that with the increase in the age of cutting seedlings, their secondary metabolism increases, and they become more sensitive to stress.

2.9. qRT-PCR Validation of Selected Genes from Blue Module

The expression levels of eight hub genes associated with plant hormone signal transduction, the plant–pathogen interaction, and phenylpropanoid biosynthesis in the blue module in the WGCNA analysis were determined using qRT-PCR. The primers for the eight selected unigenes are listed in Supplementary Table S1. The expression patterns in the eight candidate unigenes were similar to those in the transcriptome data (Figure 9). Therefore, the qRT-PCR analysis verified the data from the RNA-seq analysis, showing that our results are high reliability.

3. Discussion

3.1. Increasing Age of the Donors Affected the Basal Stem Diameter and Plant Height of Cutting Seedlings

Cuttage propagation has a positive significance for preserving excellent germplasm resources with stress resistance from ancient P. orientalis, and tree age significantly affects the growth and stress resistance of cutting seedlings [26]. In this study, the basal stem diameter and plant height of cutting seedlings decreased significantly with the increase in the age of donors, and sown seedlings grew more vigorously than cutting seedlings, both of which were propagated from 5-year-old donors. This is consistent with the research results reported indicating that the height, diameter, and rooting ability of Larix laricina decreased with aging [27]. The height of Picea abies cutting seedlings was significantly greater than that of sown seedlings [28]. Changes in the contents of soluble sugars and chlorophyll can be used as important physiological indicators for evaluating plant growth. The soluble sugar and chlorophyll contents of Populus tremuloides decreased with aging, respectively [29]. The soluble sugar and chlorophyll contents in cutting seedlings propagated from 5-year-old P. orientalis donors were greater than those from 300- and 700-year-old trees. Collectively, donor age has effects on the growth of cutting seedlings, and there is a decreasing tendency as the age of the tree increases. This may be due to the fact that ancient trees have more complex branching systems and higher requirements for environmental conditions, such as temperature, water, and nutrients, and are more susceptible to external environmental influences [30,31], which may lead to limitations in the growth of cutting seedlings propagated from ancient trees as well.

3.2. Increasing Age of the Donors Affected the Changes in Physiological Indexes of Cutting Seedlings

Flavonols participating in cell wall hardening are important in the plant stress response [32]. Flavonoid and flavonol biosynthetic pathway-related genes are up-regulated in old Ginkgo Biloba leaves [22], which is consistent with the up-regulation of FLS2 and CYP75A in 300- and 700-year-old P. orientalis seedlings in this study. The transcription levels of the FLS2 and CYP75 family genes increased with increasing tree age. Research on Taxus chinensis came to a similar conclusion [33]. CYP75 belongs to the largest cytochrome P450 (CYP) superfamily of the plant enzyme family for plant metabolism, which participates in the biosynthesis of most plant secondary metabolites to adapt to biotic and abiotic stresses [34,35]. Both CCR and F3H are involved in flavonoid biosynthesis, and the expression level of F3H in P. orientalis cutting seedlings increased with donor age. F3H plays an important role in plant fiber development and the anthocyanin synthesis pathway [25]. Secondary metabolism in plants gradually increases with the age of the tree, which may be due to the fact that old seedlings are affected by various internal and external factors during their long-term growth and gradually develop their own defense mechanisms to cope with environmental challenges, and these defense mechanisms can be realized through the synthesis of secondary metabolites [36,37]. The secondary metabolite content and related gene expression in cutting seedlings are influenced by the donor trees. Thus, the age of the donors also has effects on the expression of secondary metabolite-related genes in P. orientalis cutting seedlings.

3.3. Regulation of Hub Gene in the Stress Resistance of Ancient Cutting Seedlings

Trees at different ages respond differently to pathogens. The age of the donor trees affects the function of the immune system in the cutting seedlings and the functioning of the resistance mechanisms to pathogens. This study showed that plant–pathogen interaction pathway-related disease-resistance proteins RPM1/RPS3 and CML were up-regulated in cutting seedlings propagated from 300- and 700-year-old donors. RPS3 can endow P. syringae strains with specific disease resistance [38]. CML is a calmodulin-like protein that promotes plant innate immunity through a flagellin-dependent signaling pathway [39]. CML in P. abies promotes the root growth of seedlings by sensing the fluctuation in cytoplasmic Ca2+ during the stress response [40]. The regulation of these genes involved in stress resistance in cutting seedlings propagated from ancient P. orientalis needs further study in the future. Seedlings propagated from cuttings are genetically identical to the donor trees, and their genetic background is similar to that of the donor trees, so their disease resistance may also be similar to that of the donor tree [41]. These results indicate that cuttage can retain the growth properties of ancient P. orientalis, especially stress resistance.
GA regulates plant growth and defense in response to various stimuli [42]. Our data showed that GA and JA synthesis-related genes such as GID2/SLY1 and JAR1 were up-regulated in cutting seedlings propagated from 300- and 700-year-old P. orientalis donors. GID2, a receptor similar to hormone-sensitive lipase in gibberellin, plays a role in the response to stresses [43,44]. JAR1 up-regulation was consistent with the research results indicating that JA-related genes were up-regulated in aged A. thaliana [45]. However, the expression levels of genes involved in the JA synthesis pathway decreased in aged G. biloba branches [46]. Therefore, the regulatory roles of JA in cutting seedlings propagated from 300- and 700-year-old ancient P. orientalis need further study. These results demonstrate that cutting propagation well preserves the ability of aged trees to enhance their resistance to biotic and abiotic stresses, as well as excellent genetic resources.
Based on the collective results of this study, we constructed framework plots showing phenotypes, physiological indicators, and molecular expression levels of sown seedlings propagated from 5-year-old donors and cutting seedlings propagated from 5-, 300-, and 700-year-old donors (Figure 10).

4. Materials and Methods

4.1. Plant Materials and Growth Conditions

At a point in the growing season (June 2021), healthy growing sown seedlings propagated from 5-year-old P. orientalis donors were selected as controls from the experimental nursery of the Chinese Academy of Forestry located in Beijing. The experimental materials were 5-year-old P. orientalis cutting seedlings propagated from donors at different ages (5-, 300-, and 700-year-old). Phenotype characterization, physiological tests, and transcriptomic analysis were carried out with three biological replicates per experiment. The freshly collected leaf material was immediately frozen in liquid nitrogen and stored at −80 °C until processing for total RNA extraction.

4.2. Phenotypic Measurements

A completely randomized design was performed throughout the experiments. Basal stem diameters of seedlings were measured at 2 cm above the ground using an INSIZE 1108-150 electronic caliper, and the heights were measured from the base of the main stem to the tip of the highest leaf with a tape measure.

4.3. Determination of Soluble Sugar and Chlorophyll Contents

Approximately 0.2 g of healthy seedling leaves were harvested for the measurement of the soluble sugar content. Soluble sugar was analyzed using an anthrone colorimetric method [47].
Chlorophyll was extracted with approximately 0.1 g ground leaf material using acetone and 90% ethanol (v/v) at 4 °C under low light intensity for 3 h until the leaf material turned blanch. Chlorophyll a, chlorophyll b, and total chlorophyll contents were determined using absorbance readings.

4.4. Determination of Flavonoids and Total Phenolics

Total flavonoid content of the prepared extracts was determined using the aluminum chloride method [48]. First, 0.02 g of finely ground leaf material was extracted with 60% ethanol (v/v) at 60 °C. Then, 108 µL of supernatant was taken from the top of the tube, and 6 µL of NaNO2 (1:20 w/v) and 6 µL of AlCl3 (1:10 w/v) was added and then incubated for 6 min at 25 °C. Then, 80 µL sodium hydroxide (1 N) was added to the mixture. Finally, after mixing thoroughly, the absorbance was measured at 510 nm with a UV visible spectrophotometer using distilled water as the blank. The standard calibration curve was prepared using Rutin ranging from 0 to 1200 µg/mL. The content of flavonoids was determined using the absorbance.
The total phenolic content was determined using the Folin–Ciocalteu colorimetric method [49]. Briefly, 0.1 g of finely ground leaf material was extracted with 60% ethanol (v/v) for 2 h with shaking at 60 °C. Then, 10 µL of the extracted sample was mixed with 50 µL of Folin–Ciocalteu reagent at room temperature for 2 min. Next, 50 µL of 5% (w/v) sodium carbonate solution was added to stop the reaction, and then 90 µL of distilled water was added to make 200 µL. The mixture was incubated at 25 °C for 10 min, and the absorbance was measured at 760 nm. Gallic acid content ranging from 0 to 300 µg/mL was prepared, and the calibration curve was obtained using a linear fit (r2 = 0.9961).

4.5. Transcriptome Sequencing Methods

Total RNA was extracted using an RNAprep Pure Plant Kit (Tiangen, Beijing, China). The quality of the RNA was assessed using gel electrophoresis, and the concentration was measured using a NanoDrop 2000 spectrophotometer (Thermo Scientific, Wilmington, DE, USA). Then, cDNA libraries were constructed with RNA and sequenced using the Illumina NovaSeq 6000 platform. Raw reads were cleaned by removing adaptors and low-quality sequences. Clean reads were then aligned to the reference genome using the spliced read mapper TopHat2 (version 2.0.7) [50] and also mapped to the de novo P. orientalis transcriptome assembly using Cufflinks (version 2.2.1) [51]. Protein-coding genes were predicted using the NCBI NR, SWISS-PROT, GO, COG, KOG, eggNOG, and KEGG databases.

4.6. Quantitative Real-Time PCR Analysis

First-strand cDNA was synthesized from 1 µg total RNA using a Prime Script™ Reagent Kit (Takara, Dalian, China). The quantitative real-time PCR was performed using KAPA SYBR FAST qPCR Master Mix (Kapa Biosystems, Boston, MA, USA). Primers were designed using Premier-BLAST [52], and the sequences are depicted in Supplementary Table S1. Thermal cycling conditions were as follows: 5 s at 95 °C and 30 s at 60 °C for 40 cycles. Actin was selected as a reference gene [53], and the relative expression data were calculated using the 2−ΔΔct method [54].

4.7. Statistical Analysis of Data

The significance of differences among the mean values for growth and physiological characteristics were calculated using one-way analysis of variance (ANOVA) and expressed as the mean ± standard error (SE), followed with a Duncan’s test at p < 0.05. All statistical analyses were performed with SPSS 22.0 (SPSS Inc., Chicago, IL, USA). Transcript levels were quantified as fragments per kilobase of transcript per million mapped reads (FPKM). DEGs were identified using the R package DESeq2 (V 1.32.0) [55]. Important DEGs were identified according to the following criteria: each paired comparison |fold-change| ≥ 2 and FDR ≤ 0.001. The number of DEGs in different comparisons was visualized using the R package UpSetR (version 1.4.0) [56]. The weighted gene co-expression network was constructed using the WGCNA package in R software. Soft thresholds were calculated using the pick soft threshold function provided by WGCNA, and the gene cluster heat map was drawn using Tbtools [57].

5. Conclusions

This study combined growth, physiological indicators, and transcriptome sequencing to comprehensively analyze the effects of P. orientalis donors at different ages on the growth of cutting seedlings. Donor age affected the basal stem diameter and plant height of cutting seedlings, and the sown seedlings grew markedly better than the cutting seedlings. With increasing donor age, the contents of soluble sugars and chlorophyll decreased, while the content of secondary metabolites increased, and the expression levels of plant hormone signal transduction, the plant–pathogen interaction, and flavone and flavonol biosynthesis pathway-related genes (JAZ, FLS, RPM1/RPS3, CML, and RPS2) were significantly up-regulated. The results demonstrated that the growth and stress resistance of P. orientalis cutting seedlings were greatly affected by the age of donors and that cutting seedlings can maintain the stress resistance of ancient trees, also indicating the reliability and stability of clonal propagation. These results provide a theoretical basis for revealing the molecular mechanisms underlying the effects of donor age on the growth of cutting seedlings and are important guidelines for the conservation and management of excellent ancient tree germplasm resources using asexual propagation. Future research needs to combine the photosynthetic characteristics of seedlings, other physiological indicators, and multi-omics in a comprehensive analysis to investigate the effects of ancient tree age on cutting seedlings in depth.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/plants12091754/s1, Table S1: Primers for the candidate genes designed for qRT-PCR.

Author Contributions

Conceptualization, Y.D., W.G. and E.C.; data curation, Y.D., W.G. and E.C.; formal analysis, Y.D. and W.G.; funding acquisition, E.C.; investigation, Y.D., W.G, W.X., J.L. and Z.J. (Zirui Jia); project administration, E.C.; resources, Y.D., W.G., X.Z., Z.J. (Zeping Jiang) and W.X.; software, Y.D., W.G. and E.C.; supervision, W.G., E.C. and W.X.; writing—original draft, Y.D. and W.G.; writing—review and editing, Y.D., W.G. and E.C. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported by the State Natural Science Fund projects (31870664), and the open fund of the State Key Laboratory of Subtropical Silviculture (SKLSS-KF2022-09). The National Non-Profit Research Institutions of the Chinese Academy of Forestry [Grant numbers: CAFYBB2020ZD002-4].

Data Availability Statement

The raw sequencing data generated from this study have been deposited in NCBI SRA (https://www.ncbi.nlm.nih.gov/sra, accessed on 15 April 2023) under the accession number PRJNA947522.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Nolan, V.; Reader, T.; Gilbert, F.; Atkinson, N. The ancient tree inventory: A summary of the results of a 15 year citizen science project recording ancient, veteran and notable trees across the UK. Biodivers. Conserv. 2020, 29, 3103–3129. [Google Scholar] [CrossRef]
  2. Zhang, S.; Zhang, L.; Chai, Y.; Wang, F.; Li, Y.; Su, L.; Zhao, Z. Physiology and proteomics research on the leaves of ancient Platycladus orientalis (L.) during winter. J. Proteomics 2015, 126, 263–278. [Google Scholar] [CrossRef] [PubMed]
  3. Pizarro, A.; Díaz-Sala, C. Cellular dynamics during maturation-related decline of adventitious root formation in forest tree species. Physiol. Plant. 2019, 165, 73–80. [Google Scholar] [CrossRef] [PubMed]
  4. Muñoz-Gutiérrez, L.; Vargas-Hernández, J.J.; López-Upton, J.; Soto-Hernández, M. Effect of cutting age and substrate temperature on rooting of Taxus globosa. New Forest. 2009, 38, 187–196. [Google Scholar] [CrossRef]
  5. Swamy, S.L.; Puri, S.; Singh, A.K. Effect of auxins (IBA and NAA) and season on rooting of juvenile and mature hardwood cuttings of Robinia pseudoacacia and Grewia optiva. New Forest. 2002, 23, 143–157. [Google Scholar] [CrossRef]
  6. Ďurkovič, J.; Husárová, H.; Javoříková, L.; Čaňová, I.; Šuleková, M.; Kardošová, M.; Lukáčik, I.; Mamoňová, M.; Lagaňa, R. Physiological, vascular and nanomechanical assessment of hybrid poplar leaf traits in micropropagated plants and plants propagated from root cuttings: A contribution to breeding programs. Plant Physiol. Bioch. 2017, 118, 449–459. [Google Scholar] [CrossRef]
  7. Moura, R.D.S.; Soares, T.L.; Lima, L.K.S.; Gheyi, H.R.; Jesus, O.N.; Coelho Filho, M.A. Salinity-induced changes in biometric, physiological and anatomical parameters of Passiflora edulis Sims plants propagated by different methods. Arch. Agron. Soil Sci. 2020, 66, 1692–1706. [Google Scholar] [CrossRef]
  8. Rteiu, A.M.V.; Ştef, R.; Grozea, I.; Butnariu, M. Allelopathy potential of aesculus hippocastanum extracts assessed by phytobiologicaltest method using zea mays. Environ. Eng. Manag. J. 2015, 14, 1313–1321. [Google Scholar]
  9. Butnariu, M.; Samfira, I.; Kasapsaraçoğlu, I.; Adina, N.; Negrea, P. Allelopathic effects of Pteridium aquilinum alcoholic extract on seed germination and seedling growth of Poa pratensis. Allelopath. J. 2015, 35, 227–236. [Google Scholar]
  10. Chang, E.; Yao, X.; Zhang, J.; Deng, N.; Jiang, Z.; Shi, S. De novo characterization of Platycladus orientalis transcriptome 2 and analysis of its gene expression during aging. PeerJ Prepr. 2017, 5, e2866v1. [Google Scholar] [CrossRef]
  11. Turfan, N.; Alay, M.; Sariyildiz, T. Effect of tree age on chemical compounds of ancient Anatolian black pine (Pinus nigra subsp. pallasiana) needles in Northwest Turkey. Iforest Biogeosciences For. 2018, 11, 406. [Google Scholar] [CrossRef]
  12. Bostan, C.; Butnariu, M.; Butu, M.; Ortan, A.; Butu, A.; Rodino, S.; Parvu, C. Allelopathic effect of Festuca rubra on perennial grasses. Rom. Biotech. Lett. 2013, 18, 8190–8196. [Google Scholar]
  13. Barbat, A.; Rodino, S.; Petrache, P.; Butu, M.; BUTNARIU, M. Microencapsulation of the allelochemical compounds and study of their release from different products. Dig. J. Nanomater. Bios. 2013, 8, 945–953. [Google Scholar]
  14. Samfira, I.; BUTNARIU, M.; Rodino, S.; Butu, M. Structural investigation of mistletoe plants from various hosts exhibiting diverse lignin phenotypes. Dig. J. Nanomater. Bios. 2013, 8, 1679–1686. [Google Scholar]
  15. Alba, T.M.; Tessaro, E.; Sobottka, A.M. Seasonal effect on phenolic content and antioxidant activity of young, mature and senescent leaves from Anredera cordifolia (Ten.) Steenis (Basellaceae). Braz. J. Biol. 2022, 84, e254174. [Google Scholar] [CrossRef]
  16. Dong, Y.; Xiao, W.; Guo, W.; Liu, Y.; Nie, W.; Huang, R.; Tan, C.; Jia, Z.; Liu, J.; Jiang, Z.; et al. Effects of donor ages and propagation methods on seedling growth of Platycladus orientalis (L.) Franco in winter. Int. J. Mol. Sci. 2023, 24, 7170. [Google Scholar] [CrossRef]
  17. Bruns, E.L.; Antonovics, J.; Carasso, V.; Hood, M. Transmission and temporal dynamics of anther-smut disease (Microbotryum) on alpine carnation (Dianthus pavonius). J. Ecol. 2017, 105, 1413–1424. [Google Scholar] [CrossRef]
  18. Dunnell, K.L.; LeBoldus, J.M. The correlation between septoria leaf spot and stem canker resistance in hybrid poplar. Plant Dis. 2017, 101, 464–469. [Google Scholar] [CrossRef]
  19. Huang, L.; Zeng, Y.; Yang, S.; Zhou, H.; Xu, J.; Zhou, Y.; Wang, G. Transcriptome analysis of gene expression profiles reveals wood formation mechanisms in Chinese fir at different stand ages. Heliyon 2023, 9, e14861. [Google Scholar] [CrossRef]
  20. Guo, H.; Zhong, Q.; Tian, F.; Zhou, X.; Tan, X.; Luo, Z. Transcriptome analysis reveals putative induction of floral initiation by old leaves in tea-oil tree (Camellia oleifera ‘changlin53’). Int. J. Mol. Sci. 2022, 23, 13021. [Google Scholar] [CrossRef]
  21. Chang, E.; Zhang, J.; Deng, N.; Yao, X.; Liu, J.; Zhao, X.; Jiang, Z.; Shi, S. Transcriptome differences between 20- and 3,000-year-old Platycladus orientalis reveal that ROS are involved in senescence regulation. Electron. J. Biotechn. 2017, 29, 68–77. [Google Scholar] [CrossRef]
  22. Wang, L.; Cui, J.; Jin, B.; Zhao, J.; Xu, H.; Lu, Z.; Li, W.; Li, X.; Li, L.; Liang, E.; et al. Multifeature analyses of vascular cambial cells reveal longevity mechanisms in old Ginkgo biloba trees. Proc. Natl. Acad. Sci. USA 2020, 117, 2201–2210. [Google Scholar] [CrossRef] [PubMed]
  23. Fonseca, J.P.; Lakshmanan, V.; Boschiero, C.; Mysore, K.S. The pattern recognition receptor FLS2 Can Shape the Arabidopsis rhizosphere microbiome β-diversity but not EFR1 and CERK1. Plants 2022, 11, 1323. [Google Scholar] [CrossRef] [PubMed]
  24. Zou, Y.; Wang, S.; Zhou, Y.; Bai, J.; Huang, G.; Liu, X.; Zhang, Y.; Tang, D.; Lu, D. Transcriptional regulation of the immune receptor FLS2 controls the ontogeny of plant innate immunity. Plant Cell 2018, 30, 2779–2794. [Google Scholar] [CrossRef]
  25. Dai, M.; Kang, X.; Wang, Y.; Huang, S.; Guo, Y.; Wang, R.; Chao, N.; Liu, L. Functional characterization of Flavanone 3-Hydroxylase (F3H) and its role in anthocyanin and flavonoid biosynthesis in mulberry. Molecules 2022, 27, 3341. [Google Scholar] [CrossRef]
  26. Singh, B.; Yadav, R.; Bhatt, B.P. Effects of mother tree ages, different rooting mediums, light conditions and auxin treatments on rooting behaviour of Dalbergia sissoo branch cuttings. J. Forestry Res. 2011, 22, 53–57. [Google Scholar] [CrossRef]
  27. Greenwood, M.S.; Hopper, C.A.; Hutchison, K.W. Maturation in Larch 1: I. Effect of age on shoot growth, foliar characteristics, and DNA methylation. Plant Physiol. 1989, 90, 406–412. [Google Scholar] [CrossRef]
  28. Gemmel, P.; Orlander, G.; Hogberg, K.A. Norway spruce cuttings perform better than seedlings of the same genetic origin. Silvae Genet. 1991, 40, 198–202. [Google Scholar] [CrossRef]
  29. Donaldson, J.R.; Stevens, M.T.; Barnhill, H.R.; Lindroth, R.L. Age-related shifts in leaf chemistry of clonal aspen (Populus tremuloides). J. Chem. Ecol. 2006, 32, 1415–1429. [Google Scholar] [CrossRef]
  30. Serkova, A.A.; Tarelkina, T.V.; Galibina, N.A.; Nikerova, K.M.; Moshchenskaya, Y.L.; Sofronova, I.N.; Nikolaeva, N.N.; Ivanova, D.S.; Semenova, L.I.; Novitskaya, L.L. Changes in the differentiation program of birch cambial derivatives following trunk girdling. Forests 2022, 13, 1171. [Google Scholar] [CrossRef]
  31. Gilhen-Baker, M.; Roviello, V.; Beresford-Kroeger, D.; Roviello, G.N. Old growth forests and large old trees as critical organisms connecting ecosystems and human health. A review. Environ. Chem. Lett. 2022, 20, 1529–1538. [Google Scholar] [CrossRef]
  32. Daryanavard, H.; Postiglione, A.E.; Mühlemann, J.K.; Muday, G.K. Flavonols modulate plant development, signaling, and stress responses. Curr. Opin. Plant Biol. 2023, 72, 102350. [Google Scholar] [CrossRef]
  33. Zhang, L.; Sun, X.; Wilson, I.W.; Shao, F.; Qiu, D. Identification of the genes involved in anthocyanin biosynthesis and accumulation in Taxus chinensis. Genes 2019, 10, 982. [Google Scholar] [CrossRef]
  34. Singh, A.; Panwar, R.; Mittal, P.; Hassan, M.I.; Singh, I.K. Plant cytochrome P450s: Role in stress tolerance and potential applications for human welfare. Int. J. Biol. Macromol. 2021, 184, 874–886. [Google Scholar] [CrossRef]
  35. Das, A.; Begum, K.; Akhtar, S.; Ahmed, R.; Tamuli, P.; Kulkarni, R.; Banu, S. Genome-wide investigation of Cytochrome P450 superfamily of Aquilaria agallocha: Association with terpenoids and phenylpropanoids biosynthesis. Int. J. Biol. Macromol. 2023, 234, 123758. [Google Scholar] [CrossRef]
  36. Liu, Z.; Carpenter, S.B.; Bourgeois, W.J.; Yu, Y.; Constantin, R.J.; Falcon, M.J.; Adams, J.C. Variations in the secondary metabolite camptothecin in relation to tissue age and season in Camptotheca acuminata. Tree Physiol. 1998, 18, 265–270. [Google Scholar] [CrossRef]
  37. Rao, M.J.; Ahmed, U.; Ahmed, M.H.; Duan, M.; Wang, J.; Wang, Y.; Wang, L. Comparison and quantification of metabolites and their antioxidant activities in young and mature leaves of sugarcane. ACS Food Sci. Technol. 2021, 1, 362–373. [Google Scholar] [CrossRef]
  38. Wang, X.; Feng, H.; Chang, Y.; Ma, C.; Wang, L.; Hao, X.; Li, A.; Cheng, H.; Wang, L.; Cui, P.; et al. Population sequencing enhances understanding of tea plant evolution. Nat. Commun. 2020, 11, 4447. [Google Scholar] [CrossRef]
  39. Zhang, J.; Zou, A.; Wen, Y.; Wei, X.; Liu, C.; Lv, X.; Ma, X.; Fan, G.; Sun, X. SlCML55, a novel Solanum lycopersicum calmodulin-like gene, negatively regulates plant immunity to Phytophthora pathogens. Sci. Hortic.-Amst. 2022, 299, 111049. [Google Scholar] [CrossRef]
  40. Wen, Z.; Terhonen, E.; Asiegbu, F.O. The dark septate endophyte Phialocephala sphaeroides confers growth fitness benefits and mitigates pathogenic effects of Heterobasidion on Norway spruce. Tree Physiol. 2022, 42, 891–906. [Google Scholar] [CrossRef]
  41. Kumar, P.; Patel, P.K.; Sonkar, M.K. Propagation through juvenile shoot cuttings in difficult-to-root Dalbergia latifolia—examining role of endogenous IAA in adventitious rooting. Plant Physiol. Rep. 2022, 27, 242–249. [Google Scholar] [CrossRef]
  42. Khan, N.; Bano, A.; Ali, S.; Babar, M.A. Crosstalk amongst phytohormones from planta and PGPR under biotic and abiotic stresses. Plant Growth Regul. 2020, 90, 189–203. [Google Scholar] [CrossRef]
  43. Schwechheimer, C. Understanding gibberellic acid signaling—Are we there yet? Curr. Opin. Plant Biol. 2008, 11, 9–15. [Google Scholar] [CrossRef] [PubMed]
  44. Hirano, K.; Ueguchi-Tanaka, M.; Matsuoka, M. GID1-mediated gibberellin signaling in plants. Trends Plant Sci. 2008, 13, 192–199. [Google Scholar] [CrossRef] [PubMed]
  45. Wasternack, C. 9—Jasmonates—Biosynthesis and role in stress responses and developmental processes. In Plant Cell Death Processes; Noodén, L.D., Ed.; Academic Press: San Diego, CA, USA, 2004; pp. 143–155. ISBN 978-0-12-520915-1. [Google Scholar] [CrossRef]
  46. Yan, J.; Zhang, S.; Tong, M.; Lu, J.; Wang, T.; Xu, Y.; Li, W.; Wang, L. Physiological and genetic analysis of leaves from the resprouters of an old Ginkgo biloba tree. Forests 2021, 12, 1255. [Google Scholar] [CrossRef]
  47. Wu, L.; Deng, Z.; Cao, L.; Meng, L. Effect of plant density on yield and quality of perilla sprouts. Sci. Rep. 2020, 10, 9937. [Google Scholar] [CrossRef]
  48. Shraim, A.M.; Ahmed, T.A.; Rahman, M.M.; Hijji, Y.M. Determination of total flavonoid content by aluminum chloride assay: A critical evaluation. LWT 2021, 150, 111932. [Google Scholar] [CrossRef]
  49. Cicco, N.; Lanorte, M.T.; Paraggio, M.; Viggiano, M.; Lattanzio, V. A reproducible, rapid and inexpensive Folin–Ciocalteu micro-method in determining phenolics of plant methanol extracts. Microchem. J. 2009, 91, 107–110. [Google Scholar] [CrossRef]
  50. Kim, D.; Pertea, G.; Trapnell, C.; Pimentel, H.; Kelley, R.; Salzberg, S.L. TopHat2: Accurate alignment of transcriptomes in the presence of insertions, deletions and gene fusions. Genome Biol. 2013, 14, R36. [Google Scholar] [CrossRef]
  51. Trapnell, C.; Roberts, A.; Goff, L.; Pertea, G.; Kim, D.; Kelley, D.R.; Pimentel, H.; Salzberg, S.L.; Rinn, J.L.; Pachter, L. Differential gene and transcript expression analysis of RNA-seq experiments with TopHat and Cufflinks. Nat. Protoc. 2012, 7, 562–578. [Google Scholar] [CrossRef]
  52. Goh, H.-H.; Abu Bakar, S.; Kamal Azlan, N.D.; Zainal, Z.; Mohd Noor, N. Transcriptional reprogramming during Garcinia-type recalcitrant seed germination of Garcinia mangostana. Sci. Hortic.-Amst. 2019, 257, 108727. [Google Scholar] [CrossRef]
  53. Marum, L.; Miguel, A.; Ricardo, C.P.; Miguel, C. Reference gene selection for quantitative real-time PCR normalization in Quercus suber. PLoS ONE 2012, 7, e35113. [Google Scholar] [CrossRef]
  54. Livak, K.J.; Schmittgen, T.D. Analysis of relative gene expression data using real-time quantitative PCR and the 2−ΔΔCT method. Methods 2001, 25, 402–408. [Google Scholar] [CrossRef]
  55. Wang, L.; Feng, Z.; Wang, X.; Wang, X.; Zhang, X. DEGseq: An R package for identifying differentially expressed genes from RNA-seq data. Bioinformatics 2010, 26, 136–138. [Google Scholar] [CrossRef]
  56. Conway, J.R.; Lex, A.; Gehlenborg, N. UpSetR: An R package for the visualization of intersecting sets and their properties. Bioinformatics 2017, 33, 2938–2940. [Google Scholar] [CrossRef]
  57. Chen, C.; Chen, H.; Zhang, Y.; Thomas, H.R.; Frank, M.H.; He, Y.; Xia, R. TBtools: An integrative toolkit developed for interactive analyses of big biological data. Mol. Plant 2020, 13, 1194–1202. [Google Scholar] [CrossRef]
Plants 12 01754 g001 550
Figure 1. The 5-year-old sown seedlings and cutting seedlings propagated from 5-, 300-, and 700-year-old P. orientalis donors used as research material.
Figure 1. The 5-year-old sown seedlings and cutting seedlings propagated from 5-, 300-, and 700-year-old P. orientalis donors used as research material.
Plants 12 01754 g001
Plants 12 01754 g002 550
Figure 2. Variation in basal stem diameters and heights of sown seedlings and cutting seedlings of P. orientalis propagated from differently aged donors. (a) Basal stem diameters. (b) Heights. Note: A5: sown seedlings propagated from 5-year-old donors, B5: cutting seedlings propagated from 5-year-old donors, C300: cutting seedlings propagated from 300-year-old donors, and D700: cutting seedlings propagated from 700-year-old donors. Values are the mean of three replicates of samples, and bars represent standard errors. Means with the different letters are significantly different at p < 0.05 calculated using Duncan’s multiple range test.
Figure 2. Variation in basal stem diameters and heights of sown seedlings and cutting seedlings of P. orientalis propagated from differently aged donors. (a) Basal stem diameters. (b) Heights. Note: A5: sown seedlings propagated from 5-year-old donors, B5: cutting seedlings propagated from 5-year-old donors, C300: cutting seedlings propagated from 300-year-old donors, and D700: cutting seedlings propagated from 700-year-old donors. Values are the mean of three replicates of samples, and bars represent standard errors. Means with the different letters are significantly different at p < 0.05 calculated using Duncan’s multiple range test.
Plants 12 01754 g002
Plants 12 01754 g003 550
Figure 3. Changes in the contents of soluble sugar, chlorophyll a, chlorophyll b, chlorophyll a+b, flavonoid, and total phenol of sown and cutting seedlings propagated from differently aged donors of P. orientalis. Note: A5: sown seedlings propagated from 5-year-old donors, B5: cutting seedlings propagated from 5-year-old donors, C300: cutting seedlings propagated from 300-year-old donors, and D700: cutting seedlings propagated from 700-year-old donors. Values are the mean of three replicates of samples, and bars represent standard errors. Means with different letters are significantly different at p < 0.05 calculated using Duncan’s multiple range test.
Figure 3. Changes in the contents of soluble sugar, chlorophyll a, chlorophyll b, chlorophyll a+b, flavonoid, and total phenol of sown and cutting seedlings propagated from differently aged donors of P. orientalis. Note: A5: sown seedlings propagated from 5-year-old donors, B5: cutting seedlings propagated from 5-year-old donors, C300: cutting seedlings propagated from 300-year-old donors, and D700: cutting seedlings propagated from 700-year-old donors. Values are the mean of three replicates of samples, and bars represent standard errors. Means with different letters are significantly different at p < 0.05 calculated using Duncan’s multiple range test.
Plants 12 01754 g003
Plants 12 01754 g004 550
Figure 4. Gene ontology (GO) classifications of the assembled unigenes. The main functional annotations in the biological process, cellular component, and molecular function categories are presented. Bars represent the number of P. orientalis proteins with BLASTX matches to each GO term. Note: A5: sown seedlings propagated from 5-year-old donors, B5: cutting seedlings propagated from 5-year-old donors, C300: cutting seedlings propagated from 300-year-old donors, and D700: cutting seedlings propagated from 700-year-old donors.
Figure 4. Gene ontology (GO) classifications of the assembled unigenes. The main functional annotations in the biological process, cellular component, and molecular function categories are presented. Bars represent the number of P. orientalis proteins with BLASTX matches to each GO term. Note: A5: sown seedlings propagated from 5-year-old donors, B5: cutting seedlings propagated from 5-year-old donors, C300: cutting seedlings propagated from 300-year-old donors, and D700: cutting seedlings propagated from 700-year-old donors.
Plants 12 01754 g004
Plants 12 01754 g005 550
Figure 5. Heatmap showing unigenes based on the Kyoto Encyclopedia of Genes and Genomes (KEGG) categorizations in pairwise comparisons. Significantly enriched KEGG pathways (corrected p-value < 0.05) in five comparisons (A5-/B5-year-old, C300-/D700-year-old, A5-/D700-year-old, B5-/D700-year-old, and B5/C300-year-old). Note: A5: sown seedlings propagated from 5-year-old donors, B5: cutting seedlings propagated from 5-year-old donors, C300: cutting seedlings propagated from 300-year-old donors, and D700: cutting seedlings propagated from 700-year-old donors.
Figure 5. Heatmap showing unigenes based on the Kyoto Encyclopedia of Genes and Genomes (KEGG) categorizations in pairwise comparisons. Significantly enriched KEGG pathways (corrected p-value < 0.05) in five comparisons (A5-/B5-year-old, C300-/D700-year-old, A5-/D700-year-old, B5-/D700-year-old, and B5/C300-year-old). Note: A5: sown seedlings propagated from 5-year-old donors, B5: cutting seedlings propagated from 5-year-old donors, C300: cutting seedlings propagated from 300-year-old donors, and D700: cutting seedlings propagated from 700-year-old donors.
Plants 12 01754 g005
Plants 12 01754 g006 550
Figure 6. K-means clustering of gene expression profiles. Genes with different expression patterns were grouped into 20 clusters. The number of genes in each cluster and the functional classification of each cluster are shown. Note: A5: sown seedlings propagated from 5-year-old donors, B5: cutting seedlings propagated from 5-year-old donors, C300: cutting seedlings propagated from 300-year-old donors, and D700: cutting seedlings propagated from 700-year-old donors.
Figure 6. K-means clustering of gene expression profiles. Genes with different expression patterns were grouped into 20 clusters. The number of genes in each cluster and the functional classification of each cluster are shown. Note: A5: sown seedlings propagated from 5-year-old donors, B5: cutting seedlings propagated from 5-year-old donors, C300: cutting seedlings propagated from 300-year-old donors, and D700: cutting seedlings propagated from 700-year-old donors.
Plants 12 01754 g006
Plants 12 01754 g007 550
Figure 7. Genes representing different functional classes and showing age-dependent changes in transcript levels: (a) plant hormone signal transduction, (b) plant–pathogen interaction, and (c) flavone and flavonol biosynthesis. Note: A5: sown seedlings propagated from 5-year-old donors, B5: cutting seedlings propagated from 5-year-old donors, C300: cutting seedlings propagated from 300-year-old donors, and D700: cutting seedlings propagated from 700-year-old donors.
Figure 7. Genes representing different functional classes and showing age-dependent changes in transcript levels: (a) plant hormone signal transduction, (b) plant–pathogen interaction, and (c) flavone and flavonol biosynthesis. Note: A5: sown seedlings propagated from 5-year-old donors, B5: cutting seedlings propagated from 5-year-old donors, C300: cutting seedlings propagated from 300-year-old donors, and D700: cutting seedlings propagated from 700-year-old donors.
Plants 12 01754 g007
Plants 12 01754 g008 550
Figure 8. (a) Correlations between modules and traits. The horizontal axis represents different samples, and the vertical axis represents different modules. The column of numbers on the left represents the number of genes in the module, and each group of data on the right represents the correlation coefficients and significant p-values between the modules and the phenotypes. Red and blue represent large and small correlations between modules and phenotypes, respectively. (b) Genes representing different functional classes and showing age-dependent changes in transcript levels in blue modules. Note: A5: sown seedlings propagated from 5-year-old donors, B5: cutting seedlings propagated from 5-year-old donors, C300: cutting seedlings propagated from 300-year-old donors, and D700: cutting seedlings propagated from 700-year-old donors.
Figure 8. (a) Correlations between modules and traits. The horizontal axis represents different samples, and the vertical axis represents different modules. The column of numbers on the left represents the number of genes in the module, and each group of data on the right represents the correlation coefficients and significant p-values between the modules and the phenotypes. Red and blue represent large and small correlations between modules and phenotypes, respectively. (b) Genes representing different functional classes and showing age-dependent changes in transcript levels in blue modules. Note: A5: sown seedlings propagated from 5-year-old donors, B5: cutting seedlings propagated from 5-year-old donors, C300: cutting seedlings propagated from 300-year-old donors, and D700: cutting seedlings propagated from 700-year-old donors.
Plants 12 01754 g008
Plants 12 01754 g009 550
Figure 9. Comparison of the expression profiles for eight representative genes measured using RNA-seq and qRT-PCR. Columns represent expression levels determined using RNA-seq in FPKM values, while lines represent expression levels determined using qRT-PCR. Note: A5: sown seedlings propagated from 5-year-old donors, B5: cutting seedlings propagated from 5-year-old donors, C300: cutting seedlings propagated from 300-year-old donors, and D700: cutting seedlings propagated from 700-year-old donors. The error bars represent the SD from three replicates. The differential expression analysis was conducted using the 2−ΔΔct method.
Figure 9. Comparison of the expression profiles for eight representative genes measured using RNA-seq and qRT-PCR. Columns represent expression levels determined using RNA-seq in FPKM values, while lines represent expression levels determined using qRT-PCR. Note: A5: sown seedlings propagated from 5-year-old donors, B5: cutting seedlings propagated from 5-year-old donors, C300: cutting seedlings propagated from 300-year-old donors, and D700: cutting seedlings propagated from 700-year-old donors. The error bars represent the SD from three replicates. The differential expression analysis was conducted using the 2−ΔΔct method.
Plants 12 01754 g009
Plants 12 01754 g010 550
Figure 10. Framework plots showing phenotypes, physiological indicators, and molecular expression levels of sown seedlings propagated from 5-year-old donors and cutting seedlings propagated from 5-, 300-, and 700-year-old donors. Red represents high values, yellow represents medium values, and blue represents low values.
Figure 10. Framework plots showing phenotypes, physiological indicators, and molecular expression levels of sown seedlings propagated from 5-year-old donors and cutting seedlings propagated from 5-, 300-, and 700-year-old donors. Red represents high values, yellow represents medium values, and blue represents low values.
Plants 12 01754 g010
Table
Table 1. De novo assembly of the P. orientalis transcriptome.
Table 1. De novo assembly of the P. orientalis transcriptome.
ItemsNumber
Total number104,764
Total base105,530,209
Average length (bp)1007.31
N501858
GC percent38.75%
BUSCO scoreC: 73.9%
Number of transcripts between 0 and 50050,912 (49%)
Number of transcripts between 501 and 100023,652 (23%)
Number of transcripts between 1001 and 15009559 (9%)
Number of transcripts >150020,641 (19%)
Table
Table 2. Functional annotation of the P. orientalis transcriptome.
Table 2. Functional annotation of the P. orientalis transcriptome.
Number of UnigenesPercentage (%)
Annotated in NR24,08050.12
Annotated in NT838017.44
Annotated in KO778316.19
Annotated in SwissProt18,48338.47
Annotated in PFAM18,54638.6
Annotated in GO20,02941.68
Annotated in KOG987720.55
Annotated in all Databases31056.46
Annotated in at least one Database26,22654.58
Total unigenes48,044100
Table
Table 3. Number of DEGs in pairwise comparisons.
Table 3. Number of DEGs in pairwise comparisons.
Specific
Up-Regulated		Specific
Down-Regulated		PPI
(Nodes)		Percentage
A5 vs. B512253847321.60%
A5 vs. C300331968312353.80%
A5 vs. D7009325564520%
B5 vs. C3005569731.70%
B5 vs. D7002718140.60%
C300 vs. D7005012182.30%
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Dong, Y.; Guo, W.; Xiao, W.; Liu, J.; Jia, Z.; Zhao, X.; Jiang, Z.; Chang, E. Effects of Different Donor Ages on the Growth of Cutting Seedlings Propagated from Ancient Platycladus orientalis. Plants 2023, 12, 1754. https://doi.org/10.3390/plants12091754

AMA Style

Dong Y, Guo W, Xiao W, Liu J, Jia Z, Zhao X, Jiang Z, Chang E. Effects of Different Donor Ages on the Growth of Cutting Seedlings Propagated from Ancient Platycladus orientalis. Plants. 2023; 12(9):1754. https://doi.org/10.3390/plants12091754

Chicago/Turabian Style

Dong, Yao, Wei Guo, Wenfa Xiao, Jianfeng Liu, Zirui Jia, Xiulian Zhao, Zeping Jiang, and Ermei Chang. 2023. "Effects of Different Donor Ages on the Growth of Cutting Seedlings Propagated from Ancient Platycladus orientalis" Plants 12, no. 9: 1754. https://doi.org/10.3390/plants12091754

APA Style

Dong, Y., Guo, W., Xiao, W., Liu, J., Jia, Z., Zhao, X., Jiang, Z., & Chang, E. (2023). Effects of Different Donor Ages on the Growth of Cutting Seedlings Propagated from Ancient Platycladus orientalis. Plants, 12(9), 1754. https://doi.org/10.3390/plants12091754

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