Abstract
Taxodium ascendens is a typical tree species with high flood tolerance, and it can generate knee roots in the wetlands. This study investigated the number and size of knee roots and the soil flooding conditions. Furthermore, we also measured physiology, biochemical responses, and the anatomical structure of knee roots and underground roots at different developmental stages. This study aimed to understand the adaptation mechanism of T. ascendens to flooding stress and the formation mechanism of the knee roots. The results showed that the formation of knee roots was significantly affected by the soil water table (P < 0.05). The middle water table was more conducive to the formation of knee roots. In the middle water table, the 1-aminocyclopropane-1-carboxylic acid (ACC) content and ACC synthase activity were significantly lower in the knee roots than in the underground roots. The knee roots at the young-aged stage showed the highest ACC oxidase activity among the development stages of the knee roots. The ethylene release rate was significantly higher in the knee roots than in the underground roots (P < 0.05). Indole-3-acetic acid (IAA) content first increased, then decreased with knee root development. The periderm cells at the apex of the knee roots were dead and had many intercellular spaces, which was beneficial for the growth of T. ascendens. In conclusion, the middle water table induced the ethylene and IAA production, which promoted the formation of knee roots, which improved roots ventilation and flooding tolerance of T. ascendens. The results obtained can provide information about mechanisms of knee roots formation and provide scientific evidence for the afforestation and management under wetland conditions.
Introduction
Nowadays, flood damage is receiving a lot of attention, as flooding phenomena tend to be closely associated with global climate changes (Hirabayashi et al., 2013; Zhou et al., 2020). Flooding is a major abiotic stress that determines species distribution, growth, and yield (Herzog et al., 2016). In China, about 6.6 × 105 km2 land is waterlogged, accounting for 6.6% of the total land area; in 2020, 16 provinces and 30.20 million people were adversely affected by floods. Uneven precipitation and poor drainage result in frequent waterlogging, seriously constrain crop growth and productivity (Kuai et al., 2015; Xiong et al., 2019). Therefore, the selection of a suitable farmland shelter forest is essential. Flooding stress can impede the gas exchange between soil and atmosphere, restricting oxygen diffusion in plant tissues. Since oxygen is critical for mitochondrial respiration, this process cannot be maintained under flood conditions, which seriously constrains plant growth (Voesenek and Baileyserres, 2013).
Flood-tolerant plant species have a series of adaptive mechanisms, such as morphological changes, physiological and biochemical reactions, that can protect the plants from flooding stress (Khan et al., 2020). For example, the structures of some Chenopodiaceae roots change to the herring-bone shape to consume less oxygen and avoid being damaged by toxins (Bouma et al., 2001). The red mangrove (Rhizophora mangle) in coastal intertidal zones formed stilt-like roots for gas exchange and prevention of oxygen loss during flooding stress (Lance and Alison, 2010; Mendez-Alonzo et al., 2015). Avicennia marina can develop finger-like pneumatophores to obtain oxygen for its underground root system under hypoxic conditions (Purnobasuki et al., 2017). Eucalyptus camaldulensis can produce numerous white adventitious roots that float on the water surface to obtain more oxygen (Sena Gomes and Kozlowski, 1980).
Taxodium ascendens is a flood-tolerant species widely distributed in the subtropical riversides, drawdown areas of reservoirs, and ponds in China (Yi et al., 2008; Li et al., 2010). T. ascendens develops knee roots to increase water tolerance under flooding conditions; it is widely cultivated in farmland shelter forests and avenue plantations (Tang et al., 2008; Du et al., 2010). Knee roots are morphological changes of underground roots, which look like a cone bulging from the ground. They result from woody bulges within the secondary root cambium, which increases the thickness of the upper surfaces of the underground roots; the irregular growth of the upper and below parts of roots leads to the special shape (Tang et al., 2008). Sakio and Yamamoto (2002) have reported that only specific flood-tolerant plant species like T. ascendens, Taxodium distichum, and Glyptostrobus pensilis develop knee roots.
Although various studies have reported the effects of waterlogging on the metabolic and physiochemical characteristics of different trees, most of the studies were conducted in greenhouse conditions and focused on small tree seedlings (Conner et al., 1997; Andrson and Pezeshki, 1999; Li et al., 2010; Wang et al., 2016). Few studies were performed using adult trees under field conditions (Sou et al., 2019). Therefore, the adaptation mechanism of T. ascendens roots to waterlogging stress is unclear because the knee roots are only generated in adult trees in the field condition. Tang et al. (2008) found that the middle water table was beneficial for the formation of the knee roots. Yamamoto (1992) reported the flooding stress stimulated ethylene production by stem bark and the apex of knee roots and hypothesized the ethylene and indole-3-acetic acid (IAA) jointly contributed to knee root production. Ethylene can promote plant growth, seed germination, as well as enhance anoxia resistance (Deikman, 1997). Especially, ethylene, as an internal gaseous hormone, is used by plants to sense shifts from aerial to aquatic environments and induce changes in plant morphology and anatomy (Visser and Voesenek, 2004; Jackson, 2007). Sundberg et al. (1991) indicated that endogenous IAA is a critical regulator in the differentiation of the vascular strands, and significantly affects cambial activity and development. 1-Aminocyclopropane-1-carboxylic acid synthetase (ACS) and 1-Aminocyclopropane-1-carboxylic acid oxidase (ACO) are essential enzymes that participate in ethylene production (Dong et al., 1992; Vall-llaura et al., 2020). Application of ACC combined with IAA on the incision of stems of T. distichum seedlings was demonstrated to strongly accelerate the formation of abnormal tracheids (Yamamoto, 1992). Tavares et al. (2018) also suggested that the sensitivity of ethylene and the balance between ethylene and IAA result in forming aerenchyma.
Therefore, we hypothesized that the flooding condition affected the knee roots formation; the synthesis of endogenous hormones can promote the development of knee roots and contribute to the adoption of T. ascendens under flooding conditions. We aimed to explore the suitable water table for the knee roots formation and the growth of T. ascendens; to compare endogenous hormone levels between different development stages of knee roots and conventional underground roots; to understand the adaptation mechanism of T. ascendens under flooding conditions and form a basis for afforestation management under wetland conditions.
Materials and Methods
Research Site
The study was performed in a 28-year-old Taxodium ascendens forest plantation (approximately 500 hectares) with a mean diameter at breast height (DBH) of 25.75 cm and an average tree height of 14.33 m. The plantation was cultivated in a wetland at the Zhaoguan Forest Farm (N32°31′41″, E119°30′35″), Jiangdu County, Jiangsu Province, China. The plant spacing was 1.5 m × 4 m. The soil type was paddy soil with poor air permeability. The plant materials were selected with permission from the Zhaoguan Forest Farm; no ethics approvals were required.
Experimental Design, Analysis of the Soil Water Table and Knee Root Number and Size
A complete-randomized block design was used in this experiment. From October 2013 to October 2016, the soil water table at each site was measured every month using polyethylene pipes (length, 2 m; inner diameter, 11 cm). The pipes were drilled with 4 mm diameter holes (the spacing between each hole was 5 cm), then vertically buried in the soil at each site. According to the underground water table (Figure 1), three experimental sites were chosen based on the soil water table. High water table site (the flooding period was more than 3 months from June to September, maximum submergence depth was about 0.6 m, the annual average water table was about −0.3 m). Middle water table site (the flooding period was 1–2 months from July to August, maximum submergence depth was about 0.2 m, the annual average water table was about −0.7 m). Low water table site (no flooding status, the average yearly water table was about −1.0 m). A total of 21 plots (7 sampling plots at each water-table site) were established. Each plot was 96 m2 (8 × 12 m).
FIGURE 1

Monthly variation of underground water table during the experimental period from October 2013 to October 2016 in different study sites. The error bars represent standard deviation (SD) of the underground water table.
The tree height, DBH, and number of knee roots were measured in October 2015. The knee roots’ height and average diameter were also measured at half the height of the knee roots. Due to the irregular surface of the knee roots, we carefully wrapped the knee roots with tape (width, 1.62 cm). The surface area was calculated according to the length and width of the tape (with no overlap in the wrapping process). The survey of underground roots with knee roots and without knee roots was randomly conducted in 15 plots (1 m × 1 m) within a range of 2 m from the trunk of the T. ascendens. The underground roots (depth, 0–50 cm) were excavated, cleaned, dried, and weighed.
Sampling of Knee Roots and Assay of Physiological Indicators
The sampling was also conducted in August 2016. Since the best growth condition and the highest knee roots numbers were observed in the middle water table, the middle water table sites were chosen to investigate the development stages of the knee roots. The knee roots were divided into three stages based on size and age: young-aged stage (KY, Figure 2A), middle-aged stage (KM, Figure 2B), and old-aged stage (KO, Figure 2C). The knee roots at the young-aged stage were less than 5 cm in height and less than 5 years of age. The knee roots at the middle-aged stage were 5–10 cm in height and 5–10 years of age. The knee roots at the old-aged stage were more than 10 cm in height and more than 10 years of age. The age of the knee roots was determined using the annual rings in the transverse section. The knee root samples were divided into the swollen (upside) and the non-swollen part (underside). Underground roots approximately 10 mm in diameter were collected as the control. After removing the bark, the fresh root tissues, including cambium and phloem, were collected for physiological and biochemical analysis.
FIGURE 2

Photographs of knee roots of young-aged stage (A), middle-aged stage (B), and old-aged stage (C). The knee roots at the young-aged stage were less than 5 cm in height and less than 5 years of age. The knee roots at the middle-aged stage were 5–10 cm in height and 5–10 years of age. The knee roots at the old-aged stage were more than 10 cm in height and more than 10 years of age.
Ethylene Release Rate of the Roots
According to Kurepin et al. (2010), the ethylene release rate was determined with minor revisions. Fresh root tissues (0.6 g) were sealed in a closed bottle, incubated for 4 h at 30°C. Then, an airtight syringe injected 0.5 mL of gas samples into a gas chromatograph (Agilent 7890A; Agilent Technologies Inc. United States). The gas chromatograph was equipped with a flame ionization detector, electron capture detector and HP-PLOT Al2O3 S column (50.0 m × 0.53 mm × 0.15 μm, Agilent Technologies Inc. United States). The initial temperature of the column was set as 100°C for 10 min. Then the column was heated to 180°C for 90 min. The carrier gas (nitrogen) flow rate was 30 mL min–1. The flow rate of hydrogen was 35 mL min–1. The airflow rate was 450 mL min–1. The external standard method was used to calculate the ethylene concentration of the samples.
1-Aminocyclopropane-1-Carboxylic Acid Content
The ACC content was determined according to Hoffman et al. (1983). The root tissue samples were cut into small pieces and mixed completely. Then, the samples (1.0 g) were homogenized in 8 mL of ethanol (95%) and centrifuged (8,000 × g) at 4°C for 15 min. The supernatant was transferred to a plastic bottle, 6 mL of ethanol (80%) was added to the residue and shaken for 30 min. After centrifugation, the supernatant was transferred to the previous bottle and dried in a water bath at 40°C. The dried residue was dissolved with 5 mL of distilled water, then centrifuged (8,000 × g) at 4°C for 10 min. The supernatant was used as the ACC fluid; 1 mL of the fluid and 40 μL of HgCl2 (25 mmol⋅L–1) were added into test tubes (20 mL) closed with rubber stoppers. Then, a syringe injected 1 mL of NaOCI–NaOH (v:v = 2:1) into the test tubes. After mixing up, gas samples (0.5 mL) were collected from the tubes to determine the ethylene concentration via gas chromatography. The ACC concentration of fresh root samples was calculated using the following formula:
where C is the ethylene concentration measured using gas chromatography (nL⋅L–1); VL is the volume of sample bottles without solution (mL); V is the volume of the extracting solution (mL); R is the transfer coefficient from ACC to ethylene; V1 is the volume of the extracting solution used for measurement (mL); V2 is the gas sample volume for gas chromatography (mL); W is the fresh weight of the root sample (g); 22.4 is 1 molar gas constant under normal atmospheric conditions (L⋅mol–1).
1-Aminocyclopropane-1-Carboxylic Acid Synthetase Activity
The ACS activity was determined according to Mehta et al. (1988). Briefly, 1 g of root tissues were ground in a mortar with the extraction buffer (2 mL) and centrifuged at 10,000 × g for 20 min at 4°C. The extraction buffer solution contained 400 mmol⋅L–1 potassium phosphate buffer solution (pH 8.5), 1 mmol⋅L–1 ethylene diamine tetraacetic acid, 0.5% (v:v) β-mercaptoethanol, and 10 μmol⋅L–1 pyridoxal phosphate (PLP). The supernatant was used to determine ACC enzyme activity; 0.4 mL of the enzyme extract and 1.6 mL of the buffer solution (containing 50 μmol⋅L–1 S-adenosylmethionine (SAM), 10 μmol⋅L–1 PLP, and 50 mmol⋅L–1 Hepes-KOH, pH 8.5) were added to the test tubes. After incubation (32°C, 1 h), 0.1 mL of mercuric chloride (500 mmol⋅L–1) stopped the reaction. The test tubes were closed, incubated in ice water for 5 min. 0.2 mL of 5% NaOCI-NaOH (v:v = 2:1) was injected into the test tubes. After shaking, a 0.5 mL gas sample was collected to determine the ethylene content using gas chromatography.
1-Aminocyclopropane-1-Carboxylic Acid Oxidase Activity
According to Dong et al. (1992), the ACO activity was determined. 0.5 g of root tissues were ground in a mortar with the extraction buffer (2 mL) and centrifuged at 12,000 × g for 10 min at 4°C. The extraction buffer solution contained 100 mmol⋅L–1 Tris-HCl (pH 7.5), 10% (v:v) glycerin, 30 mmol⋅L–1 sodium ascorbate, 5% (v:v) polyvinylpyrrolidone, 0.1 mmol⋅L–1 FeSO4, and 5 mmol⋅L–1 DTT. The supernatant was used to determine ACO activity; 0.2 mL of the enzyme extract and 1.8 mL of the buffer solution (containing 100 mmol⋅L–1 Tris-HCl (pH 7.5), 10% (v:v) glycerin, 30 mmol⋅L–1 sodium ascorbate, 30 mmol⋅L–1 NaHCO3, 1.0 mmol⋅L–1 ACC, and 0.1 mmol⋅L–1 FeSO4) were added to test tubes (20 mL) closed with rubber stoppers. After incubation (35°C, 20 min), a 0.5 mL gas sample was collected to determine the ethylene content using gas chromatography.
Indole-3-Acetic Acid Content
The IAA content was determined according to Chen and Zhao (2008) and Yuan et al. (2008). After homogenization, 1 g of root tissues were extracted in precooled 80% methanol and butylated hydroxytoluene (1 mmol⋅L–1) overnight at 4°C. After centrifugation, the supernatant was passed through a C18 Sep-Pak Cartridge (Waters Corp., Milford, United States), and dried in a rotary evaporator (RE-2000A, China). The residue was collected in 0.8 ml mobile phase (consisting of 23% (v:v) methanol and 23% (v:v) acetonitrile in double distilled water supplemented with 0.1% (v:v) phosphoric acid), filtered through a 0.25 mm filter for High Performance Liquid Chromatography (HPLC) analysis. The IAA content was determined using the external standard method with a Waters 2695 Alliance HPLC (Waters Corp.) equipped with a Symmery C18 column (Waters Corp.) (4.6 × 250 mm, 5 μm).
Anatomical Analysis of Knee Roots and Underground Roots
The thickness of the rhytidome, size of the secondary phloem, phloem parenchyma and earlywood tracheids from the apex of underground roots and knee roots at different development stages were determined according to Yamamoto (1992). The anatomical analysis was performed according to Yamamoto (1992) and Zhang et al. (2013). Three samples containing bark and currently produced xylem from the apex of knee roots at the middle-aged stage and mid-sized underground roots were obtained (Figure 3). Small pieces (10 × 10 × 10 mm) of these root materials were fixed in FAA solution (formalin: acetic acid: ethanol: water, 5:5:60:30, v:v) for 24 h, rinsed in water, dehydrated in ethanol, and sealed in a paraffin block. The samples were transversely sectioned (10 μm) with a slide microtome (Leica RM2125RT), dewaxed, stained with safranin-fast green solution, and oven-dried at 40°C. The anatomical structure of the samples was observed under a light microscope (Olympus, Japan), the pictures were taken and processed using an image analysis software (DT2000, China).
FIGURE 3

Photographs of a root segment with a knee root (left) and transverse section of a knee root (right). Vertical bars indicate the portion with the transverse section of the knee root. Letters indicate positions of samples for microscopic observations. The horizontal bar equals 2 cm. A, apex; B, base; S, side.
Statistical Analysis
All data were subjected to a one-way analysis of variance by using IBM SPSS Statistics 19.0 (IBM Corp., Armonk, NY, United States). Pearson’s correlation analysis was performed using IBM SPSS Statistics 19.0 (IBM Corp). The data were presented as mean ± standard deviation (M ± SD) values, and the means of 3 replicates were evaluated using Duncan’s test at a significance level of 0.05.
Results
Relationships Between the Morphological Characteristics of Knee Roots and Underground Water Table
The results showed that the water table significantly affected the formation and distribution of knee roots (P < 0.05). The knee root density in the middle water table was significantly higher than that in the high water table and low water table (P < 0.05). The height and diameter of the knee roots were also observed to be significantly higher in the middle water table site (P < 0.05) (Table 1).
TABLE 1
| Site | Annual water table (m) | Flooding period (month year–1) | Knee root density (root m–2) | Knee root height (cm) | Knee root diameter (cm) |
| High water table | −0.3 | > 3 | 0.66 ± 0.02 b | 6.65 ± 1.36 b | 3.73 ± 0.12 b |
| Middle water table | −0.7 | 1–2 | 1.61 ± 0.64 a | 7.64 ± 1.42 a | 4.11 ± 0.51 a |
| Low water table | −1.0 | 0 | 0.65 ± 0.02 b | 6.72 ± 1.60 b | 3.70 ± 0.25 b |
Number and size of knee roots of Taxodium ascendens at different water table sites.
Data represent mean ± SD values of seven replication sites.
Values followed by the same letter are not significantly different at P < 0.05, according to Duncan’s multiple range tests. The height of the knee roots refers to the height from the ground to the apex of the knee roots. The knee root diameter refers to the average diameter at half the height of the knee roots.
Effects of the Knee Roots on the Growth of Taxodium ascendens
As shown in Figure 4, the weight of underground roots with knee roots was significantly higher than that without knee roots (P < 0.05). Furthermore, the correlation between tree height and DBH and knee roots in the middle water table was analyzed. Our results showed that the diameter at breast height and tree height of T. ascendens were positively and significantly correlated with the number and size of knee roots (Table 2). The number and the weight of underground roots were significantly correlated with the mean height and surface area of knee roots.
FIGURE 4

Weight of underground roots with or without knee roots. Values followed by the same letter(s) are not significantly different at P < 0.05, according to Duncan’s multiple range tests. Error bars are SD of the mean; n = 3.
TABLE 2
| Number of knee roots | Mean height of knee roots | Mean diameter of knee roots | Surface area of knee roots | |
| Tree height | 0.45* | 0.46* | 0.52* | 0.52* |
| Tree DBH | 0.66** | 0.61** | 0.65** | 0.59** |
| Number of underground roots | 0.87** | 0.48 | 0.86** | |
| Weight of underground roots | 0.93** | 0.53 | 0.95** |
Correlation coefficient between number and size of knee roots and DBH and tree height of Taxodium ascendens at Zhaoguan Forest Farm.
DBH, diameter at breast height.
*P = 0.05 significance; **P = 0.01 significance.
1-Aminocyclopropane-1-Carboxylic Acid Content and 1-Aminocyclopropane-1-Carboxylic Acid Synthetase and 1-Aminocyclopropane-1-Carboxylic Acid Oxidase Enzyme Activities of the Roots
The ACC content ranged from 0.08 to 0.29 nmol g–1FW among the four stages. There was no significant difference in the ACC content among the knee roots at different development stages (Figure 5A). However, the ACC content was significantly lower in the knee roots than underground roots (P < 0.05). There were no significant differences in the ACS activity among the development stages of the knee roots (Figure 5B); underground roots exhibited significantly higher ACS activity as compared to knee roots (P < 0.05). KY showed significantly higher ACO activity among the development stages of the knee roots (P < 0.05); The ACO activity decreased with the growth of knee roots (Figure 5C). There was no significant difference between UR and KO.
FIGURE 5

ACC content (A) and ACS (B) and ACO (C) activities in different roots of Taxodium ascendens. KO, KM, and KY refer to the knee roots at different development stages: old-aged stage, middle-aged stage, and young-aged stage, respectively. UR means underground roots. ACS, 1-Aminocyclopropane-1-carboxylic acid synthetase; ACO, 1-Aminocyclopropane-1-carboxylic acid oxidase. Values followed by the same letter(s) are not significantly different at P < 0.05, according to Duncan’s multiple range tests. Error bars are SD of the mean; n = 3.
Endogenous Hormone Release Rates
Different ethylene release rates were observed in the underground roots and knee roots (Figure 6A). The ethylene release rate was significantly higher in the knee roots than in the underground roots (P < 0.05). In addition, significant differences were observed among the different development stages of the knee roots (P < 0.05). The maximum ethylene release rate from the knee roots was observed in KO, and the minimum, in KM.
FIGURE 6

Ethylene release rate (A) and IAA content (B) from different root types of Taxodium ascendens. Values followed by the same letter(s) are not significantly different at P < 0.05, according to Duncan’s multiple range tests. Error bars are SD of the mean; n = 3.
The IAA content was in the order of UR > KM > KO > KY (Figure 6B), with no statistically significant differences between UR and KM (P > 0.05). There was no significant difference between KO and KY (P > 0.05).
Anatomical Structure of the Knee Roots and Underground Roots
As shown in Table 3, the thickness of the rhytidome increased with the development of the knee roots. The thickness of the rhytidome in the knee roots was significantly greater than that in the underground roots (P < 0.05). The radial width of the secondary phloem followed the order of UR < KY < KM < KO. The size of phloem parenchyma increased with the development of the knee roots, the phloem parenchyma of KY was significantly smaller than that of UR (P < 0.05). The tangential width of the tracheids also increased with the development of the knee roots. However, the radial width of the earlywood tracheids decreased with the development of the knee roots.
TABLE 3
| KY | KM | KO | UR | ||
| Thickness of the rhytidome (mm) | 1.53 ± 0.23c | 2.51 ± 0.41b | 4.07 ± 0.58a | 0.24 ± 0.08d | |
| Radial width of the secondary phloem (mm) | 0.31 ± 0.03c | 0.45 ± 0.03b | 1.25 ± 0.03a | 0.22 ± 0.02d | |
| Size of the phloem parenchyma | Tangential width (μm) | 25.07 ± 2.82b | 26.49 ± 5.73ab | 29.78 ± 3.54a | 28.49 ± 2.93a |
| Radial width (μm) | 10.72 ± 0.87a | 11.08 ± 1.22a | 11.14 ± 1.47a | 12.75 ± 2.10a | |
| Size of the tracheids | Tangential width (μm) | 26.99 ± 2.96b | 29.47 ± 7.28b | 36.86 ± 5.49a | 29.61 ± 5.31b |
| Radial width (μm) | 37.83 ± 7.19ab | 35.39 ± 5.77b | 33.34 ± 8.39b | 40.37 ± 4.89a |
The thickness of the rhytidome, size of the secondary phloem, phloem parenchyma, and earlywood tracheids from the apex of underground roots and knee roots at different development stages.
KO, KM, and KY refer to the knee roots at different development stages: old-aged stage, middle-aged stage, and young-aged stage, respectively. UR means underground roots. Values followed by the same letter(s) are not significantly different at P < 0.05, according to Duncan’s multiple range tests. Data represent mean ± SD values of 9 replicates.
The anatomical structures of middle-aged stage knee roots and underground roots were shown in Figure 7. The knee roots had 3–4 layers of rhytidome, which were formed by integrating the periderm and phloem, and some parts of the periderm had branches. The periderm was composed of many layers of cork cells, and the cork layer, cork-forming layer, and inner cork layer were closely overlapped. The phloem, isolated from the periderm on the apex side of the knee roots, was dead, the cells were broken, and the arrangement of cells was loose and porous. The phloem parenchyma cells were close to the cambium and rectangular. The cell wall of the phloem fiber was thickened and showed an increase in lignification. The phloem ray of knee roots was more expanded than that of underground roots. The knee roots were mainly composed of the secondary xylem; in the cross-section, xylem tracheids were arranged in order; the early tracheids were rectangular, square, or polygonal; the late tracheids were smaller than the early tracheids; the cell wall was thicker.
FIGURE 7

Transverse anatomical structure of knee roots and underground roots. Horizontal line = 100 μm. (A) Rhytidome at the apex of the knee roots (middle stage); (B) rhytidome of mid-sized underground roots; (C) phloem at the apex of the knee roots (middle stage); (D) phloem of mid-sized underground roots; (E) xylem at the apex of the knee roots (middle stage); (F) xylem of mid-sized underground roots; rh, rhytidome; pe, periderm; ph, phloem; pr, phloem ray; pf, phloem fiber; pc, phloem parenchyma cell; wr, wood ray.
The rhytidome layers of the underground roots were lesser than those of the knee roots (1–2 layers). The phloem cells isolated by periderm were loosely arranged and more orderly than the knee root. The composition and arrangement of secondary phloem were the same as the knee root, but the width of the phloem was smaller than that of the knee roots. The xylem tracheids were arranged in order, and the shape was similar to the knee roots.
Discussion
In plants, flood tolerance is related to changes in anatomical and morphological characteristics (Jackson and Colmer, 2005; Hua et al., 2017). The formation of knee roots is a morphological adaptation of T. ascendens to environmental stress (Figure 8). In the present study, T. ascendens with knee roots had more underground roots and showed better tree growth, suggesting that the knee roots are beneficial for T. ascendens under flooding conditions. Our results also indicate that the middle water table is more suitable for the formation and development of knee roots, which is similar to the study of Tang et al. (2008), who reported that T. ascendens formed adventitious roots in the high water table and knee roots in the middle water table. The adventitious roots were usually attached to the tree trunk or bark gap, while knee roots were mainly vertical woody outgrowths rising from lateral roots (Rogers, 2021). The low water table was also not conducive to the formation of knee roots. The reason might be the low water table was easier for T. ascendens to obtain nutrients and oxygen; there was no need to form the knee roots to resist the flooding stress (Chen et al., 2021; Middleton et al., 2021).
FIGURE 8

Knee roots of Taxodium ascendens in the wetland.
Ethylene can trigger programmed cell death during aerenchyma formation, which is vital in modifying the response of plants to flooding stress (Mignolli et al., 2020). Ethylene synthesis includes two essential processes: catalysis of S-adenosylmethionine to ACC by ACS and catalysis of ACC to ethylene by ACO (Vanderstraeten et al., 2019). ACC can be synthesized under hypoxia conditions in the underground roots; ACS activity is increased in plant roots under flood stress (Dong et al., 1992; Zhou et al., 2001; Williams and Golden, 2002). In our study, the submerged underground roots of T. ascendens were exposed to anoxic conditions, which may stimulate the expression of the ACS gene (Van der Straeten et al., 2001; Rieu et al., 2005). Thus the significantly higher ACS activity was observed in the underground roots. However, the formation of knee roots alleviated the anoxic stress of T. ascendens, possibly resulting in significantly lower ACS activity in the knee roots. Oxygen is necessary for the process of ACC to ethylene. ACO activity can also be induced by anoxic stress (Vriezen et al., 1999; Bailey-Serres and Voesenek, 2008; Buttò et al., 2020). Thus, in the present study, ACO activity decreased with the development of knee roots, which suggested that the oxygen status was improved by the formation and development of knee roots. When the underground roots of T. ascendens were exposed to flooding in the growing season, the activity of ACS was enhanced by anaerobic stress, which led to the highest accumulation of ACC. But the lowest ethylene was observed in underground roots. We believe it was due to the middle water table and the formation of the knee roots. When the water table receded, the upper surface of the underground roots can receive oxygen easily, which contributes to the ethylene synthesis by ACO enzyme; ethylene maybe promotes the formation of knee roots, which resulted in significantly higher ethylene content in the knee roots than the underground roots (Figure 6A). The highest ethylene content was observed in the old-aged stage of the knee roots, which probably contributing to the cell death and lignification in the knee roots formation process (Pesquet and Tuominen, 2011). Besides, ethylene can indirectly regulate gene expression through the signal transduction pathway in plants tissues (Sairam et al., 2008; Zhang et al., 2014). Fan et al. (2018) reported ThRAP2.3 is one of the critical downstream-response ethylene response factor genes, which can respond to the ethylene signal under flooding stress in a Taxodium hybrid “Zhongshanshan 406.” The molecular regulation of ethylene in T. ascendens is still unknown, which merits further study.
The phytohormone IAA promotes root development and adventitious root formation (Visser and Voesenek, 2004; Kitomiy et al., 2008). IAA also participates in the development and activities of the cambium (Uggla et al., 1998; Bhalerao and Fischer, 2014). Our study suggests that, with the development of the knee roots, the IAA content first increased, then decreased; IAA content was lower in the knee roots at the young-aged and old-aged stages and higher in those at the middle-aged stage. High IAA content is beneficial for cambium cell enlargement, expansion, and division (Bhalerao and Fischer, 2016), whereas low IAA content is beneficial for secondary cell-wall deposition and lignification (Fajstavr et al., 2018). Our results were similar to the conclusion of Cui et al. (1999, 2000), who reported that endogenous IAA concentration in the cambium of Broussonetia papyrifera increased when the cambium formed the immature phloem and xylem; IAA concentration decreased when the immature vascular cells differentiated toward maturation. This result indicates IAA concentrations are different in the division stage of cambium cells and the differentiation stage of cambium derivative cells.
Under flood conditions, anoxia stimulates ethylene production and accumulation in roots, which contributed to the programmed cell death in the cortex tissue (He et al., 1996; Drew et al., 2000). In our study, the periderm cells at the apex of the knee roots were dead, arranged loosely, and had a large number of intercellular spaces, which is conducive to gas exchange between the knee roots and air. Because the knee roots were exposed to air to resist the adverse effects of external environmental factors, the periderm of the knee roots was thicker than that of the underground roots (Table 3). The phloem rays at the apex of the knee roots expanded when compared with the underground roots, possibly due to the ethylene accumulation (Pallardy, 2011).
In conclusion, consistent with our hypothesis, the flooding resistance and growth condition of T. ascendens are related to the formation of knee roots and enhancement of air permeability. Our study also suggested middle water table was conducive to the formation of the knee roots compared with the high water table and low water table. The ethylene and IAA may affect the formation of knee roots (Figure 9). In this study, we realized that the adaptation mechanism of T. ascendens to flooding stress is preliminary. The underlying formation mechanism of knee roots is still unclear. We recommend further studies of the formation mechanism of knee roots should be explored to understand the mechanism of flooding resistance.
FIGURE 9

A schematic summary figure showing the adaption mechanism of Taxodium ascendens to flooding stress.
Conclusion
In the present study, the middle water table significantly enhanced the formation and distribution of knee roots in T. ascendens. Furthermore, the knee roots were beneficial for the development of T. ascendens. The underground roots showed significantly higher ACC content and ACS activities than knee roots. Besides, the ethylene release rate was significantly higher in the knee roots than in the underground roots. The highest IAA content was observed in the middle-aged knee roots. The thickness of the rhytidome increased with the development of the knee roots. Cells of the periderm at the apex of the knee roots were dead, arranged loosely, and had a large number of intercellular spaces that improved internal gas diffusion. In conclusion, the middle water table induced IAA and ethylene production, which changed the morphology and anatomy of T. ascendens roots. The formation of knee roots can improve roots ventilation and growth of T. ascendens.
Publisher’s Note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Statements
Data availability statement
The original contributions presented in the study are included in the article/supplementary material, further inquiries can be directed to the corresponding author.
Author contributions
LW and LT designed the experiment. ZQ wrote the draft. LT and ZQ revised the manuscript. LW and ZQ carried out the field and laboratory experiments. All authors have contributed to the corrections of the manuscript and approve it for submission.
Funding
The funders had no role in the experiment design, data analysis, decision to publish or preparation of the manuscript. This study was funded by the National Natural Science Foundation of China (No. 31170566) and the Priority Academic Program Development of Jiangsu Higher Education Institutions (PAPD).
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
References
1
AndrsonP. H.PezeshkiS. R. (1999). The effects of intermittent flooding on seedlings of three forest species.Photosynthetica37543–552. 10.1023/A:1007163206642
2
Bailey-SerresJ.VoesenekL. A. C. J. (2008). Flooding stress: acclimations and genetic diversity.Annu. Rev. Plant Biol.59313–339. 10.1146/annurev.arplant.59.032607.092752
3
BhaleraoR. P.FischerU. (2014). Auxin gradients across wood—instructive or incidental?Physiol. Plant15143–51. 10.1111/ppl.12134
4
BhaleraoR. P.FischerU. (2016). Environmental and hormonal control of cambial stem cell dynamics.J. Exp. Bot.6879–87. 10.1093/jxb/erw466
5
BoumaT. J.NielsenK. L.Van HalJ.KoutstaalB. (2001). Root system topology and diameter distribution of species from habitats differing in inundation frequency.Funct. Ecol.15360–369. 10.1046/j.1365-2435.2001.00523.x
6
ButtòV.DeslauriersA.RossiS.RozenbergP.ShishovV.MorinH. (2020). The role of plant hormones in tree-ring formation.Trees34315–335. 10.1007/s00468-019-01940-4
7
ChenD.ZhaoJ. (2008). Free IAA in stigmas and styles during pollen germination and pollen tube growth of Nicotiana tabacum.Physiol. Plantarum.134202–215. 10.1111/j.1399-3054.2008.01125.x
8
ChenZ.ArifM.WangC.ChenX.LiC. (2021). Effects of hydrological regime on foliar decomposition and nutrient release in the riparian zone of the three gorges reservoir, China.Front. Plant Sci.12:661865. 10.3389/fpls.2021.661865
9
ConnerW. H.McLeodK. W.McCarronJ. K. (1997). Flooding and salinity effects on growth and survival of four common forested wetland species.Wetl. Ecol. Manag.599–109. 10.1023/A:1008251127131
10
CuiK.WangX.DuanJ. (2000). Change in the concentration and tissue - localization of endogenous IAA during the regeneration after girdling in Broussonetia papyrifera (L.) Vent.Acta Scientiarum. Natural. Univ. Pekinensis36495–502. 10.13209/j.0479-8023.2000.071(in Chinese),
11
CuiK.WangX.DuanJ.WangS. Y. (1999). Change of endogenous IAA during cambial activity and IAA binding protein in cambial cells in Broussonetia papyrifera.Acta Bot. Sin.411082–1085. 10.3321/j.issn:1672-9072.1999.10.011(in Chinese),
12
DeikmanJ. (1997). Molecular mechanisms of ethylene regulation of gene transcription.Physiol. Plant100561–566. 10.1034/j.1399-3054.1997.1000318.x
13
DongJ. G.Fernández-MaculetJ. C.YangS. F. (1992). Purification and characterization of 1-aminocyclopropane-1-carboxylate oxidase from apple fruit.Proc. Natl. Acad. Sci. U.S.A.899789–9793. 10.1073/pnas.89.20.9789
14
DrewM. C.HeC. J.MorganP. W. (2000). Programmed cell death and aerenchyma formation in roots.Trends Plant Sci.5123–127. 10.1016/s1360-1385(00)01570-3
15
DuQ.ZhongQ. C.WangK. Y. (2010). Root effect of three vegetation types on shoreline stabilization of Chongming Island, Shanghai.Pedosphere20692–701. 10.1016/S1002-0160(10)60059-8
16
FajstavrM.PaschováZ.GiagliK.VavrèíkH.GrycV.UrbanJ. (2018). Auxin (IAA) and soluble carbohydrate seasonal dynamics monitored during xylogenesis and phloemogenesis in Scots pine.iForest11553–562. 10.3832/ifor2734-011
17
FanW.YangY.WangZ.YinY.YuC.ShiQ.et al (2018). Molecular cloning and expression analysis of three ThERFs involved in the response to waterlogging stress of Taxodium ‘Zhongshanshan406’, and subcellular localization of the gene products.Peer J.6:e4434. 10.7717/peerj.4434
18
HeC. J.FinlaysonS. A.DrewM. C.JordanW. R.MorganP. W. (1996). Ethylene biosynthesis during aerenchyma formation in roots of maize subjected to mechanical impedance and hypoxia.Plant Physiol.1121679–1685. 10.1104/pp.112.4.1679
19
HerzogM.StrikerG. G.ColmerT. D.PedersenO. (2016). Mechanisms of waterlogging tolerance in wheat - a review of root and shoot physiology.Plant Cell Environ.391068–1086. 10.1111/pce.12676
20
HirabayashiY.MahendranR.KoiralaS.KonoshimaL.YamazakiD.WatanabeS.et al (2013). Global flood risk under climate change.Nat. Clim. Change3816–821. 10.1038/NCLIMATE1911
21
HoffmanN. E.LiuY.YangS. F. (1983). Changes in 1- (malonylamino) cyclopropane-1-carboxylic acid content in wilted wheat leaves in relation to their ethylene production rates and 1-aminocyclopropane-1-carboxylic acid content.Planta157518–523. 10.1007/BF00396882
22
HuaJ.HanL.WangZ.GuC.YinY. (2017). Morpho-anatomical and photosynthetic responses of Taxodium hybrid ‘Zhongshanshan’ 406 to prolonged flooding.Flora23129–37.
23
JacksonM. B. (2007). Ethylene-promoted elongation: an adaptation to submergence stress.Ann. Bot. Lond.101229–248. 10.1093/aob/mcm237
24
JacksonM. B.ColmerT. D. (2005). Response and adaptation by plants to flooding stress. Ann. Bot. Lond. 96, 501–505. 10.1093/aob/mci205
25
KhanM.IqbalR.TrivelliniA.ChhillarH.ChopraP.FerranteA.et al (2020). The significance and functions of ethylene in flooding stress tolerance in plants.Environ. Exp. Bot.179:104188. 10.1016/j.envexpbot.2020.104188
26
KitomiyY.OgawaA.KitanohI.InukaiY. (2008). CRL4 regulates crown root formation through auxin transport in rice.Plant Root.219–28. 10.3117/plantroot.2.19
27
KuaiJ.ZhouZ.WangY.MengY.ChenB.ZhaoW. (2015). The effects of short-term waterlogging on the lint yield and yield components of cotton with respect to boll position.Eur. J. Agron.6761–74. 10.1016/j.eja.2015.03.005
28
KurepinL. V.WaltonL. J.ReidD. M.PharisR. P.ChinnappaC. C. (2010). Growth and ethylene evolution by shade and sun ecotypes of Stellaria longipes in response to varied light quality and irradiance.Plant Cell Environ.29647–652. 10.1111/j.1365-3040.2005.01443.x
29
LanceS. E.AlisonB. (2010). Characterization of cork warts and aerenchyma in leaves of Rhizophova mangle and Rhizophora racemosa.J. Torrey Bot. Soc.13730–38. 10.2307/40864968
30
LiC.ZhongZ.GengY.SchneiderR. (2010). Comparative studies on physiological and biochemical adaptation of Taxodium distichum and Taxodium ascendens seedlings to different soil water regimes.Plant Soil.329481–494. 10.1007/s11104-009-0174-z
31
MehtaA. M.JordanR. L.AndersonJ. D.MattooA. K. (1988). Identification of a unique isoform of 1-aminocyclopropane -1-carboxylic acid synthase by monoclonal antibody.Proc. Natl. Acad. Sci. U.S.A.858810–8814. 10.1073/pnas.85.23.8810
32
Mendez-AlonzoR.MoctezumaC.OrdonezV. R.AngelesG.MartinezA. J.Lopez-PortilloJ. (2015). Root biomechanics in Rhizophora mangle: anatomy, morphology and ecology of mangrove’s flying buttresses.Ann. Bot Lond.115833–840. 10.1093/aob/mcv002
33
MiddletonB. A.LeiT.VillegasO.LiuX. (2021). Regeneration trends along climate gradients in Taxodium distichum forests of the southeastern United States.For. Ecol. Manag.497:119485. 10.1016/j.foreco.2021.119485
34
MignolliF.TodaroJ. S.VidozM. L. (2020). Internal aeration and respiration of submerged tomato hypocotyls are enhanced by ethylene-mediated aerenchyma formation and hypertrophy.Physiol. Plantarum.16949–63. 10.1111/ppl.13044
35
PallardyS. G. (2011). Physiology of Woody Plants.Beijing: Science Press.
36
PesquetE.TuominenH. (2011). Ethylene stimulates tracheary element differentiation in Zinnia elegans cell cultures.New Phytol.190138–149. 10.1111/j.1469-8137.2010.03600.x
37
PurnobasukiH.PurnamaP. R.KobayashiK. (2017). Morphology of four root types and anatomy of root-root junction in relation gas pathway of Avicennia marina (Forsk) Vierh roots.Int. J. Plant Res.30:100. 10.5958/2229-4473.2017.00143.4
38
RieuI.CristescuS. M.HarrenF. J. M.HuibersW.VoesenekL. A. C. J.MarianiC.et al (2005). RP-ACS1, a flooding-induced 1-aminocyclopropane-1-carboxylate synthase gene of Rumex palustris, is involved in rhythmic ethylene production.J. Exp. Bot.56841–849. 10.1093/jxb/eri078
39
RogersG. K. (2021). Bald Cypress knees, Taxodium distichum (Cupressaceae): an anatomical study, with functional implications.Flora278:151788. 10.1016/j.flora.2021.151788
40
SairamR. K.KumuthaD.EzhilmathiK.DeshmukhP. S.SrivastavaG. C. (2008). Physiology and biochemistry of waterlogging tolerance in plants.Bio Plantarum52401–412. 10.1007/s10535-008-0084-6
41
SakioH.YamamotoF. (2002). Ecology of Riparian Forests.Tokyo: University of Tokyo Press.
42
Sena GomesA. R.KozlowskiT. T. (1980). Effects of flooding on Eucalyptus camaldulensis and Eucalyptus globulus seedlings.Oecologia46139–142. 10.1007/BF00540117
43
SouH. D.MasumoriM.KurokochiH.TangeT. (2019). Histological observation of primary and secondary aerenchyma formation in adventitious roots of Syzygium kunstleri (King) Bahadur and R.C.Gaur grown in hypoxic medium.Forests10:137. 10.3390/f10020137
44
SundbergB.LittleC. H.CuiK.SandbergG. (1991). Level of endogenous indole-3-acetic acid in the stem of Pinus sylvestris in relation to the seasonal variation of cambial activity.Plant Cell Environ.14241–246. 10.1111/j.1365-3040.1991.tb01342.x
45
TangL. Z.HuangB. L.HaibaraK.TodaH. (2008). Ecological adaption mechanisms of roots to flooded soil and respiration characteristics of knee roots of Taxodium ascendens.J. Plant Ecol.321258–1267. 10.3773/j.issn.1005-264x.2008.06.006 (in Chinese).
46
TavaresE. Q. P.GrandisA.LembkeC. G.SouzaG. M.PurgattoE.De SouzaA. P.et al (2018). Roles of auxin and ethylene in aerenchyma formation in sugarcane roots.Plant Signal. Behav.13:e1422464. 10.1080/15592324.2017.1422464
47
UgglaC.MellerowiczE. J.SundbergB. (1998). Indole-3-acetic acid controls cambial growth in Scots pine by positional signaling.Plant Physiol.117113–121. 10.1104/pp.117.1.113
48
Vall-llauraN.Giné-BordonabaJ.UsallJ.LarrigaudièreC.TeixidóN.TorresR. (2020). Ethylene biosynthesis and response factors are differentially modulated during the interaction of peach petals with Monilinia laxa or Monilinia fructicola.Plant Sci.299:110599. 10.1016/j.plantsci.2020.110599
49
Van der StraetenD.ZhouZ.PrinsenE.Van OnckelenH. A.Van MontaguM. C. (2001). A comparative molecular physiological study of submergence response in lowland and deepwater rice.Plant Physiol.125955–968. 10.2307/4279720
50
VanderstraetenL.DepaepeT.BertrandS.Van Der StraetenD. (2019). The ethylene precursor ACC affects early vegetative development independently of ethylene signaling.Front. Plant Sci.10:1591. 10.3389/fpls.2019.01591
51
VisserE. J. W.VoesenekL. A. C. J. (2004). Acclimation to soil flooding – sensing and signal-transduction.Plant Soil.254197–214. 10.1007/s11104-004-1650-0
52
VoesenekL. A.BaileyserresJ. (2013). Flooding tolerance: O2 sensing and survival strategies.Curr. Opin. Plant Biol.16647–653. 10.1016/j.pbi.2013.06.008
53
VriezenW. H.HulzinkR.MarianiC.VoesenekL. A. C. J. (1999). 1-Aminocyclopropane-1-carboxylate oxidase activity limits ethylene biosynthesis in Rumex palustris during submergence.Plant Physiol.121189–195. 10.1104/pp.121.1.189
54
WangC.LiC.HongW.WeiX.HanW. (2016). Effects of long-term periodic submergence on photosynthesis and growth of Taxodium distichum and Taxodium ascendens saplings in the hydro-fluctuation zone of the three gorges reservoir of China.PLoS One11:e0162867. 10.1371/journal.pone.0162867
55
WilliamsO. J.GoldenK. D. (2002). Purification and characterization of ACC oxidase from Artocarpus altilis.Plant Physiol. Biochem.40273–279. 10.1016/s0981-9428(02)01379-7
56
XiongQ.CaoC.ShenT.ZhongL.HeH.ChenX. (2019). Comprehensive metabolomic and proteomic analysis in biochemical metabolic pathways of rice spikes under drought and submergence stress.BBA Proteins Proteomics1867237–247. 10.1016/j.bbapap.2019.01.001
57
YamamotoF. (1992). Effects of depth of flooding on growth and anatomy of stems and knee roots of Taxodium distichum.IAWA J.1393–104. 10.1163/22941932-90000560
58
YiY. H.FanD. Y.XieZ. Q.ChenF. Q. (2008). The effects of waterlogging on photosynthesis-related eco-physiological processes in the seedlings of Quercus variabilis and Taxodium ascendens.Acta Ecol. Sin.286025–6033. 10.3321/j.issn:1000-0933.2008.12.030
59
YuanJ.ChenD.RenY. J.ZhangX. L.ZhaoJ. (2008). Characteristic and expression analysis of a metallothionein gene, OsMT2b, down-regulated by cytokinin suggest functions in root development and seed embryo germination of rice (Oryza sativa L.).Plant Physiol.1461637–1650. 10.1104/pp.107.110304
60
ZhangJ.GouX.ZhaoZ.LiuW.ZhangF.CaoZ.et al (2013). Improved method of obtaining micro-core paraffin sections in dendroecological research.J. Plant Ecol.37972–977. 10.3724/SP.J.1258.2013.00100(in Chinese),
61
ZhangJ.YuJ.WenC. K. (2014). An alternate route of ethylene receptor signaling.Front. Plant Sci.5:648. 10.3389/fpls.2014.00648
62
ZhouW.ChenF.MengY.ChandrasekaranU.LuoX.YangW.et al (2020). Plant waterlogging/flooding stress responses: from seed germination to maturation.Plant Physiol. Biochem.148228–236. 10.1016/j.plaphy.2020.01.020
63
ZhouZ.VriezenW.Van CaeneghemW.MontaguM. V.StraetenD. V. D. (2001). Rapid induction of a novel ACC synthase gene in deepwater rice seedlings upon complete submergence.Euphytica121137–143. 10.1023/A:1012059425624
Summary
Keywords
anatomic structure, ethylene, flooding, formation mechanism, knee root, Taxodium ascendens
Citation
Qian Z, Wu L and Tang L (2022) Effects of Flooding and Endogenous Hormone on the Formation of Knee Roots in Taxodium ascendens. Front. Plant Sci. 13:803619. doi: 10.3389/fpls.2022.803619
Received
28 October 2021
Accepted
11 January 2022
Published
03 February 2022
Volume
13 - 2022
Edited by
Bhumi Nath Tripathi, Indira Gandhi National Tribal University, India
Reviewed by
Ghazala Mustafa, Quaid-i-Azam University, Pakistan; Guzel Kudoyarova, Institute of Biology of Karelian Research Centre, Russian Academy of Sciences (RAS), Russia
Updates
Copyright
© 2022 Qian, Wu and Tang.
This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Luozhong Tang, luozhongtang@njfu.edu.cn
This article was submitted to Plant Abiotic Stress, a section of the journal Frontiers in Plant Science
Disclaimer
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.