White oak (Quercus fabri Hance) regenerated stump sprouts show few senescence symptoms during 40 years’ growth in a natural forest

Background: The relationship between the parental tree’s physiological age and the clonal offspring’s lifespan is unclear. White oak (Quercus fabri Hance) has a high sprouting capability after harvest, with the regenerated sprouts being typical clonal individuals. To determine whether regenerated sprouts undergo rapid senescence compared with the parent, the senescence levels of 5-, 10-, 20- and 40-year-old regenerated stump sprouts in a natural forest were evaluated. The antioxidative abilities and transcriptomes in these leaves and shoots were compared. Results: Older regenerated sprouts still had robust antioxidative systems, with 40-year-old sprouts having lower peroxidation product levels but similar antioxidative enzyme activity levels compared with 5-year-old sprouts. Older leaves had greater transcriptional activities in pathways related to cell growth and division than younger leaves. However, older sprouts had some unhealthy characteristics, such as increased base excision repair levels and upregulated phagosome, proteasome and glycerophospholipid metabolism pathways in 40-year-old leaves, which indicates that DNA damage and tissue remodeling occurred more frequently than in younger leaves. Additionally, plant–pathogen interactions and MAPK signals pathways were upregulated in older shoots, which indicates that older shoots suffered from more pathogen-related biotic stress. Conclusions: The 40-year-old sprouts still had the same vitality level as the 5-year-old sprouts, although the former had some unhealthy characteristics. We conclude that during their rst 40 years of growth, regenerated stump sprouts do not begin to senesce, and that the parental tree’s physiological age does not signicantly affect its clonal offspring’s lifespan.


Background
Sprouts are the main material of vegetative propagation, and many species, including European aspens (Populus tremula L.), Jabuticaba (Plinia jaboticaba Vell) and Populus euphratica, are vegetatively propagated as sprouts (Kramp et al. 2018;Oliveira et al. 2019;Tullus et al. 2020). Previously, it was reported that older parental trees had lower sprouting capabilities (Matula et al. 2012). However, the link between parental tree age and regenerated sprout lifespan is unclear; do regenerated sprouts have shorter lifespans than sexually reproduced plants? Clonal trees are usually found to mature faster, and a trees' growth is limited by the maturation period (Greenwood 1995;Greenwood et al. 2010). This may occur because the clonal materials are mature; therefore, rejuvenation is implemented to delay or reverse maturation (Wendling et al. 2014). Ally et al. (2010) hypothesized that somatic mutations from parents are not eliminated by homologous recombination and natural selection. Consequently, ramets carry somatic mutations, which result in lower fertility rates. However, sprouts grow from calli, in which genes are reprogrammed. The telomere lengths also increase during callus culturing in species like Arabidopsis thaliana, maize, and agave (Rescalvo-Morales et al. 2016;Santos-Serejo and Aguiar-Perecin 2016;Sováková et al. 2018). Thus, the aging status may be reversed in calli. In Scots pine (Pinus Sylvestris Linn.), although DNA methylation levels increase with meristematic age, the growth rates and reproductive capacities of grafted shoots do not differ between old and young donors (Mencuccini et al. 2014). Owing to their long lengths, lifespans of vegetatively propagated perennial trees are hard to measure; consequently, limited studies have been published.
Accelerating senescence and leaf abscission to facilitate the survival of the next generation are strategies of monocarpic plants (mainly annual or biennial) in response to stress, and the senescence of an entire organism occurs readily in monocarpic species (Sade et al. 2018). Unlike monocarpic plants, polycarpic plants (perennial) grow inde nitely, and plants produce new ramets to compensate for lost ones. Therefore, senescence at the organismal level rarely occurs, except when the formation of new ramets cannot offset losses (Munné-Bosch 2015). However, the accelerated senescence of leaves and roots is still harmful to polycarpic species because nutrient acquisition is curtailed. A decline in reproductive capacity is another signal of organismal senescence. White campion (Silene latifolia), a short-lived herb, produces fewer in orescences in old age (Pujol et al. 2014), and the number of viable pollen grains per catkin per ramet decreases with clonal age in male Populus tremuloides (Ally et al. 2010). In Elymus excelsus, a bunchgrass, the reproductive capacity, but not the vegetative growth, decreases with age .
For plant organs, senescence is the nal developmental stage. During this stage, several life-support systems start to decline, tissues stop growing and programmed cell death occurs on a large scale (Thomas 2013). The decline of the antioxidative system is one characteristic of senescence. Decreased antioxidase activity and antioxidant levels lead to reactive oxygen species (ROS) accumulation. ROS react with DNA and metabolites, activate several senescence-regulated NAC transcription factors, induce so-called oxidative damage and ultimately kill the cells (Berni et al. 2019;Kikuchi et al. 2020). In addition, the accumulation of gene mutations and telomere shortening are ongoing during organ life. The DNA repair system continuously upregulates to maintain genetic stability, and apoptosis occurs when DNA damage is irreparable (Abate et al. 2020). Senescence-associated genes, such as those involved in plant hormone signaling (abscisic acidresponsive element-binding transcription factors 2, 3 and 4) (Gao et al. 2016), epigenetic modi cations (Arabidopsis JmjC-domain containing protein 16) (Liu et al. 2019) and pathogen defense (Senescenceassociated genes 13) (Dhar et al. 2020), also play roles in senescence.
White oak (Quercus fabri Hance) is widely distributed in the northern hemisphere, and it takes more than 50 years to grow from a sapling to harvest size (Schweitzer et al. 2019). When white oaks are cut down, the stumps have a high sprouting capability (Weigel and Peng 2002). In this study, to explore whether parental trees' physiological ages in uenced the lifespans of the regenerated sprouts, the senescence levels of 5-, 10-, 20-and 40-year-old white oak regenerated stump sprouts in a natural forest were evaluated and the antioxidative capabilities and transcriptomes in leaves and shoots were compared. We hypothesized that if regenerated sprouts have shorter lifespans, then the leaves and shoots of 40-year-old sprouts should show obvious signs of senescence.

Study sites and plants
The study sites were located at Fengshushan Forestry Farm (29°11′N, 117°32′E) in Jiangxi Province, Southern China, which experiences a subtropical monsoon climate. The mean annual , monthly minimum (January) and monthly maximum (July) temperatures are 17.8°C, 10.9°C and 41.8°C, respectively. The average length of the growing season is 263 days, and the mean annual precipitation is 1,805 mm. All the white oaks form one population, with sprouts being regenerated from stumps and allowed to grow for 5, 10, were the sprouts' main-stem diameters in the same age group. Fresh leaves and shoots attached to leaves on regenerated sprouts were collected in mid-July. We used pruning shears to collect leaf and shoot samples.
Because the samples were collected from one population, there were only four samples in both the 5a and 40a groups (n = 4) and six samples in both the 10a and 20a groups (n = 6). All the leaves and shoots were stored on dry ice immediately after collection, transported to the laboratory and then stored at −80°C.

Estimates of antioxidative capabilities
A certain weight of leaf or shoot sample was crushed in a prechilled pestle and crushed with a mortar in an appropriate volume of prechilled extraction buffer. After uniform homogenization and centrifugation at 4,000 g for 15 min at 4°C, the supernatant was collected and used for biochemical analyses.
Antioxidative capabilities were estimated by quantifying peroxidation products and antioxidase activity levels.

RNA sequencing and analysis
Total RNA was extracted using TRIzol reagent (Invitrogen, California, USA) following the manufacturer's instructions. Quality was checked using gel electrophoresis and an Agilent 2100 Bioanalyzer (Agilent Technologies, CA, USA). Pooled libraries were sequenced on an Illumina HiSeq platform with 2× 150-bp paired-end reads. After removing the adaptor and low-quality reads, the remaining clean reads were de novo assembled using Trinity software (Haas et al. 2013). Differentially expressed genes (DEGs) were determined using DESeq (R package) with the criteria fold change > 2 and p-value < 0.05 (Anders and Huber 2010). All the genes were annotated in Kyoto Encyclopedia of Genes and Genomes (KEGG) databases for functional classi cations. The KEGG enrichment was performed using TBtools (Chen et al. 2020).
The gene expression pattern analysis was perform using Short Time-series Expression Miner software (Ernst and Bar-Joseph 2006) on the OmicShare tools platform (www.omicshare.com/tools).

Statistical analysis
Data were evaluated using a one-way ANOVA followed by Tukey's multiple comparisons tests. Values are expressed as means ± SEMs.

Results
Oxidative stress does not increase in older regenerated sprouts Oxidative stress-related indices were different for leaves and shoots, and they were not linearly correlated with sprout ages (Figure 1). In leaves, SOD and POD activities showed no signi cant changes among the four tested ages, while the CAT activities in the 10a and 40a groups were lower than in the 5a and 20a groups. The MDA, H 2 O 2 and O 2 − levels were lowest in the 40a group. In shoots, SOD and CAT activity levels in the 10a group were lower than in the other groups, while the POD activity in the 20a group was greater than in the 5a and 40a groups. However, the MDA, H 2 O 2 and O 2 − levels in shoots were not signi cantly different among the four tested ages. Although the enzyme activities differed among the four age groups, the transcriptional levels of the SOD, POD and CAT families were not signi cantly changed ( Figure 2).

Enrichment analysis reveals some unhealthy characteristics
After removing pathways with no signi cance (p > 0.05) and those that were only distinguished in a single comparison, a KEGG enrichment analysis of DEGs indicated that several pathways changed during sprout growth (Tables 1 and 2). Programmed cell death-related pathways, including apoptosis, autophagy, and DNA mismatch repair, were not markedly upregulated in older organs (Signorelli et al. 2019).
In leaves, cyanoamino acid metabolism, avonoid biosynthesis, phenylpropanoid biosynthesis, brassinosteroid biosynthesis, starch and sucrose metabolism, cutin, suberine and wax biosynthesis, and fatty acid elongation pathways varied across age, showing changes in four of six comparisons (Table 1 and S1).
The majority of DEGs enriched in pathways related to growth (starch and sucrose metabolism, cutin, suberine and wax biosynthesis, fatty acid elongation and RNA transport), cell division (DNA replication and both porphyrin and chlorophyll metabolism), and death (base excision repair, phagosome, proteasome and glycerophospholipid metabolism) were upregulated at greater ages (Figure 4a).
The number of DEGs in shoots was an order of magnitude lower than in leaves, with only nine pathways varying across age, and no altered pathways were found in the S5:S40 comparison (Table 2 and   S2). Half the DEGs in the starch and sucrose metabolism pathway were downregulated in older shoots, while the other half were upregulated. This trend was unlike that in leaves (Figure 4b). The plant-pathogen interaction pathway changed with sprout growth. The majority of the DEGs in this pathway were downregulated at the 5a-10a stage but upregulated at the 10a-20a and 20a-40a stages. Accordingly, to respond to pathogens, the phagosome pathway was upregulated at 20a. The MAPK signaling pathway, a senescence-related pathway , was upregulated at the 10a-20a and 20a-40a stages. The phagosome pathway was upregulated in the S20 group, and the plant hormone signal transduction pathway was upregulated in older shoots (Figure 4b).
Genes were separated on the basis of expression pro les In leaves, six gene expression pro les were identi ed during regenerated sprout growth ( Figure 5). The down-trending expression pro les (0/2 and 7) contained genes related to energy, sugar, lipid, amino acid and nitrogen metabolism, which decreased with sprout age. The up-trending expression pro les (10 and 17/19) contained genes related to DNA replication, DNA repair and signal transduction, which increased with sprout age. Expression pro le 18 contained genes related to MAPK signal and chlorophyll metabolism that were highly expressed in the 20a group.
In shoots, ve gene expression pro les were identi ed during regenerated sprout growth ( Figure 6). The down-trending expression pro les (0 and 2) contained genes related to transporter, amino acid metabolism and nitrogen metabolism, which decreased with sprout age. The up-trending expression pro les (12 and 17/19) contained genes related to signal transduction, plant-pathogen interactions and translation, which increased with sprout age. Expression pro les 11/18 contained genes related to the ubiquitin system, cytochrome P450 and proteasome that were highly expressed in the 20a group.

Discussion
Perennial plants grow continuously throughout their lives, with the meristem producing and rejuvenating new organs to compensate for those that are lost. The organism prevents senescence through continuous growth, which make it di cult to observe whole perennial plant senescence (Munné-Bosch 2015). The early aging of new meristematic products, that is, leaves and shoots, can indirectly prove plant senescence (Ryser et al. 2020). To ascertain whether regenerated sprouts succeed the parent's age, which indicates that regenerated sprouts have shorter lifespans and age faster than sexually reproduced sprouts, the senescence levels of leaves and shoots over 40 years of growth were estimated. We did not nd enough physiological or transcriptional evidence to support that regenerated sprouts senesce in 40 years.
Senescence is a programmed cell death phenomenon that involves oxidative stress, which is indicated by a decrease in antioxidative enzyme activity levels and peroxidation product accumulation. The peroxidation products lead to further DNA damage and cell membrane disruption (Juvany et al. 2013;Sade et al. 2018). In plants, respiration and photosynthetic activities are vigorous in summer. High levels of oxygen free radicals are produced in electron transport chains, and the antioxidative defense system is important in preventing oxidative damage (Juvany et al. 2013). Here, several antioxidative enzymes showed lower activities in older sprouts, but peroxidation products did not accumulate (Figures 1 and 2). This indicated that the 35-year period only slightly in uenced the antioxidative capabilities both at the transcriptional and physiological levels. The 40-year-old regenerated sprouts still maintained the capacity to eliminate peroxidation products.
This trend was different from in senescent plant organs. During Tulipa gesneriana petal senescence, antioxidative enzyme activity and ROS levels increase , while, ROS levels increased but antioxidative enzyme activity levels decreased during rice leaf senescence .
During cell senescence, macromolecules tend to be catabolized rather than synthesized; consequently, nutrients are transferred to healthy cells (Avila-Ospina et al. 2014;Qi et al. 2020;Shinozaki et al. 2020). For example, during the bamboo shoot postharvest stage, pathways, including autophagy, fatty acid degradation and amino acid metabolism, are upregulated and accompanied by organism senescence . In this study, the main differences between ages was related to metabolism (Figures 4-6), with older regenerated sprouts having greater synthetic reactions in biosynthetic pathways, including avonoid, brassinosteroid, wax, fatty acids and amino acids, compared with younger ones, and even the gene expression and cell division (RNA transport and DNA replication) levels increased with age (Table 1, Figure  4a). These transcriptional changes indicate that 40-year-old sprouts still produce enough materials to maintain a large organism. However, the base excision repair pathway was upregulated in the L40 group, which indicates that there was more DNA damage in 40-year-old leaves . Furthermore, the upregulation of phagosome, proteasome and glycerophospholipid metabolism suggested that tissue remodeling was accompanied by DNA damage in the L40 group (Figure 4a). Although the autophagy pathway was not activated in the 40a group, the accumulation of DNA damage was an indicator of declining cell health, because it induces terminal programed cell death when the DNA damage exceeds the repair system's threshold (Martin and Wood 2019;Yamada and Coffman 2005).
Unlike leaves, shoots act as transporters rather than producers, and they have longer lifespans in deciduous plants. Thus, strong metabolic activities are not required for shoot cell production and growth; consequently, the transcriptional differences in shoots were much smaller than in leaves. The main transcriptional variation was in the plant-pathogen interaction pathway, which changed with age (Table 2), which indicates that shoots faced age-related biotic stresses. Both effector-triggered immunity and pathogen associated molecular patterns (PMAP)-triggered immunity (PTI) were in uenced by sprout age (Table S2). The PTI responses produce ROS to inhibit pathogen growth (Boller and Felix 2009), and the trends in shoot SOD and CAT activity levels were similar to the trend in the plant-pathogen interaction pathway. PTI responses activate the MAPK pathway (Meng and Zhang 2013), and in this study, the MAPK pathway was simultaneously upregulated in older shoots (Figure 4b). MAPK signals induce events, including plant senescence, ROS generation and hypersensitive response, in Arabidopsis thaliana (Kawasaki et al. 2017;Matsuoka et al. 2015;Zhou et al. 2009).
In natural forests, sprouts are not arti cially rejuvenated. Therefore, if regenerated sprouts inherit the parent's age-related characteristics, then senescence symptoms should be more conspicuous. However, we did not nd evidence that sprouts had inherited their parent's age-related characteristics. This has been observed in other species. A 10-year phenological study in Brazil found that Diospyros lasiocalyx trees in three forests reproduced mainly through vegetative propagation, and widespread senescence in the D. lasiocalyx populations was not observed (Aguiar et al. 2020). Pyrenean saxifrage is a monocarpic perennial species that escapes monocarpic senescence at high altitudes by clonally producing lateral rosettes that share the same root, which results in the plant becoming polycarpic (Cotado and Munné-Bosch 2020). Lomatia tasmanica is an endangered species with only one population, and vegetative propagation is the only L. tasmanica reproduction method. The clonal age has been determined to be 43,600 years (Lynch et al. 1998).
Nevertheless, advanced clonal ages diminish or even delete the sexual propagation capacity, as seen in species such as L. tasmanica (Lynch et al. 1998), P. tremuloides (Ally et al. 2010) and Lemna minor (Barks and Laird 2015). Although these species still propagate vegetatively and are not senescent, the lack of sexual propagation decreases gene functional redundancy and population stress resistance.
Many forest tree species lose their regenerative capacities with age or maturation. Trees cut at more advanced ages have lower organ regenerative capabilities or regenerate organs having the parental characteristics resulting in low growth rates (Diaz-Sala 2014). However, a study reported that meristems from old-growth trees grafted onto juvenile rootstock are just as vigorous as those from mid-aged trees, although the former exhibited the old-tree parental characteristics (Greenwood et al. 2010). This indicated that meristem vigor does not decline with age and that the excepted lifespans of clonal trees are restrained by natural selection rather than genetic factors. Another study found that trees with a habitat-speci c optimum body size have minimal mortality probability and maximal fertility in populations (Dani and Kodandaramaiah 2019). In this study, we hypothesized that two factors protected white oak sprouts from senescence, as follows: (1) sprouts were regenerated from stumps and their body size were smaller than parental, which indicated that the habitat was suitable for sprout growth and allowed the sprouts to absorb nutrients through the parental roots; and (2) callus cells reversed age-related effects during sprout formation, and consequently, the regenerated sprouts did not receive the parent's age-related characteristics.
Undeniably, this study has some defects, as follows.
(1) Technique-related limitations. The ages of sprout parents may not be consistent and errors in predicted ages may in uence the sprout phenotypes. All the parental stumps had similar diameters, and the trees are from one population. The similar genetic background ensures that any age-related errors are not large.
(2) The white oak lifespan is more than 400 years; therefore, a long time period is required for senescence symptoms to become obvious (Johnson and Abrams 2009). The results in this study indicate that 40-year-old sprouts are very healthy and are unlikely to senesce soon. (3) The experimental forest in the present study is located in a relatively comfortable environment and faces less abiotic and biotic stresses than other forests, such as windbreak forests and commercial plantations. The effects of stress resistance on regenerated sprouts have not been estimated; therefore, climate, interspeci c competition and nutrient limitations may accelerate the aging of sprouts, especially clonal trees of the same genotype (Sloan and Sayer 2015;Zhang et al. 2018).

Conclusions
We found that in a natural forest, regenerated stump sprouts did not show negative signs of aging in the rst 40 years of growth. Older sprouts still had robust antioxidative systems, while endocrine disorders, strong catabolism and massive autophagy were undetected, compared with younger sprouts. However, older sprouts had a few unhealthy characteristics, such as increased levels of DNA damage in leaf cells and more frequent tissue remodeling. Additionally, shoots suffered more biotic stress from pathogens, and MAPK signals were activated. The data indicated that although invisible injuries were accumulating, 40-year-old sprouts were still in the prime of their lives. The parental tree's physiological age did not signi cantly shorten the lifespans of the clonal offspring. Thus, although regenerated sprouts were not as strong as parental trees, they maintain the parental genotypes and help to sustain genetic diversity in harvested forests. PTI: pathogen associated molecular patterns (PMAP)-triggered immunity.

Declarations
Ethics approval and consent to participate Not applicable.

Consent for publication
Not applicable.

Availability of data and materials
Data sharing is not applicable to this article as no datasets were generated or analyzed during the current study.

Competing interests
The authors declare that they have no competing interests.

Author contributions
Honggang Sun and Liwen Wu contributed to the conception and design of the study; Sisheng Wu organized the data and performed the statistical analysis; Honggang Sun and Sisheng Wu wrote the manuscript. All authors approved the submitted version.

Figure 6
Gene expression patterns in white oak shoots. Each box represents an expression pro le. The upper and lower numbers in the boxes are the pro le serial number and the p-value, respectively. Boxes of the same color represent similar expression pro les having high correlation rates. Bar plots represent the KEGG enrichment results for the genes in the pro le, with the 20 most signi cant pathways being represented.

Supplementary Files
This is a list of supplementary les associated with this preprint. Click to download. TableS2.xlsx