ARTICLE   Open Access    

Responses of isolated balsam-fir stem segments to exogenous ACC, IAA, and IBA

More Information
  • Received: 30 April 2024
    Revised: 06 September 2024
    Accepted: 06 September 2024
    Published online: 30 September 2024
    Forestry Research  4 Article number: e033 (2024)  |  Cite this article
  • In this investigation, the effects of exogenous indole-3-acetic acid (IAA), indole-3-butyric acid (IBA), and 1-aminocyclopropane-1-carboxylic acid (ACC) on anatomical development within cultured segments of Abies balsamea (L.) Mill. were compared, using debudded and defoliated leaders produced in the preceding year as bioassay material. In stem apical regions, IAA promoted radial enlargement of pre-existing cortical resin ducts and attending parenchyma enlargement, whereas IBA promoted cell division and expansion of parenchyma on the outer edge of phloem without altering cortical duct shape. Cortical woody ducts, each partially surrounded by cambium, were observed as a novel but infrequent feature. A single cortical woody duct was spatially associated with each mature leaf as its vascular trace, and they were not encountered elsewhere in the cortex, nor were they induced to form in response to any hormone application. An unknown leaf factor induces the development of cortical woody ducts. Both IAA and IBA promoted cell division in the vascular cambium. The common cellular response at the interface between the latewood boundary and cambial zone was the radial expansion of primary-walled fusiform cambial cell derivatives with little if any ensuing tracheary element (TE) differentiation. Enhanced TE production at basal stem positions occurred when ACC was provided with IAA and/or IBA, and an IAA + IBA + ACC combination produced a basal stem response similar to that in untreated segments having intact leaves. The data support the conclusion that IAA, IBA, and ACC have distinct but complementary roles in the overall regulation of the types of cellular differentiation that contribute to cortex histogenesis and diameter growth of balsam-fir leaders.
  • The competition for consumer preference for fresh apples (Malus domestica) from exotic and tropical fruits is intense. Red-fleshed (RF) apple may not only provide a novel point of differentiation and enhanced visual quality, but also a source of increased concentration of potentially health-benefiting compounds within both the fresh fruit and snack/juice markets[1]. Two different types of RF apples have been characterised: Type 1 RF apple has red colouration not only in the fruit core and cortex, but also in vegetative tissues, including stems and leaves; Type 2 RF apples display red pigment only in the fruit cortex[1, 2]. To facilitate trade and lengthen the supply-window, harvested fruit are usually cold stored, which can induce a series of disorders, including physiological breakdown manifesting as a flesh browning disorder (FBD) in RF apples[3, 4]. FBD in RF apples can be caused by senescence, and there is also some evidence to suggest that a large proportion of RF apples are chilling-sensitive (Jason Johnston, Plant & Food Research Hawke's Bay, personal communication).

    Earlier studies suggested that Type 1 RF colour was determined by a promoter mutation of MdMYB10 that has a tandem replication of a myeloblastosis (MYB) binding cis-element (R6) within the promoter, resulting in autoregulation of MdMYB10[5]. However, van Nocker et al.[6] observed a large variation in the degree and pattern of red pigmentation within the cortex among the accessions carrying MdMYB10, and concluded that the presence of this gene alone was not sufficient to ensure the RF colour. A genome-wide association study (GWAS)[3], reported that, in addition to the MdMYB10 gene, other genetic factors (e.g. MdLAR1, a key enzyme in the flavonoid biosynthetic pathway) were associated with RF colour, too. Wang et al.[7] reported that many of the up-regulated genes in RF apples were associated with flavonoid biosynthesis (e.g., chalcone synthase (CHS), chalcone isomerase (CHI), dihydroflavonol 4-reductase (DFR), anthocyanin synthase (ANS), UDP-glucosyltransferase (UGT) and MYB transcription factors). Recently, MdNAC42 was shown to share similar expression patterns in RF fruit with MdMYB10 and MdTTG1, and it interacts with MdMYB10 to participate in the regulation of anthocyanin synthesis in the RF apple Redlove®[8].

    Several transcription factor genes (e.g., MYB, WRKY, bHLH, NAC, ERF, bZIP and HSF) were reported to be differentially expressed during cold-induced morphological and physiological changes in 'Golden Delicious' apples[9]. A study by Zhang et al.[10] showed that ERF1B was capable of interacting with the promotors of anthocyanin and proanthocyanidin (PA) MYB transcription factors, and suggested that ethylene regulation and anthocyanin regulation might be linked in either direction. It was reported that ethylene signal transduction pathway genes or response genes, such as ERS (ethylene response sensor), EIN3 (ethylene-insensitive3) and ERFs (ethylene response factors), may play an important role in the regulatory network of PA biosynthesis[11].

    Espley et al.[12] observed no incidence of FBD in cold-stored fruit of 'Royal Gala', but over-expression of MdMYB10 in 'Royal Gala' resulted in a high rate of FBD in RF fruit, which was hypothesised to be caused by elevated fruit ethylene concentrations before harvest and more anthocyanin, chlorogenic acid (CGA) and pro-cyanidins in RF fruit. In addition, the MYB10 transcription factor was shown to elevate the expression levels of MdACS, MdACO, and MdERF106 ethylene-regulating genes[12]. To elucidate the mechanism regulating the FBD of RF apples, Zuo et al.[13] analysed the transcriptome of tree-ripe apples at 0, 0.5 and 4 h after cutting, and reported that the differentially expressed genes at different sampling points were mainly related to plant–pathogen interactions.

    GWAS is a powerful technique for mining novel functional variants. One of the limitations of GWAS, using SNP arrays, is that they require genotyping of large numbers of individuals, which may be expensive for large populations. DNA pooling-based designs (i.e., bulk segregant analysis) test differential allele distributions in pools of individuals that exhibit extreme phenotypes (XP) in bi-parental populations, large germplasm collections or random mating populations[1416]. In addition to reducing the number of samples to be genotyped, the use of whole genome sequencing (WGS)-based XP-GWAS has the potential to identify small-effect loci and rare alleles via extreme phenotypic selection.

    In this WGS-based XP-GWAS, we investigated the genetic basis of RF and FBD by sequencing the pools of individuals that exhibited extreme phenotypes for these two traits, and analysed the differences in allele frequencies between phenotypic classes. This method combines the simplicity of genotyping pools with superior mapping resolution. We also examined the transcriptome from transgenic apple fruit harbouring the R6:MYB10 promoter as a model for red flesh in apple. Differences in gene expression of a highly pigmented line were compared with expression in control fruit and these genes were then used for comparison with the seedling population. Understanding the genetic basis of the link between RF and FBD will help in design of strategies for selection against FBD in high-quality Type 1 RF apple cultivars.

    A snapshot of visual variation in FBD and WCI is presented in Fig. 1. The average WCI and FBD across all ~900 seedlings ranged from 0 to 7, and from 0% to 58%, respectively. Based on the MLM analysis, the estimated narrow-sense heritability (h2) of WCI and FBD was 0.57 (standard error = 0.18) and 0.09 (standard error = 0.05), respectively. The estimated genetic correlation between WCI and FBD was 0.58, and several fruit quality traits displayed unfavourable correlations with WCI (Supplemental Fig. S1). Seedlings with higher WCI scores were generally characterised by poor firmness and crispness, plus higher astringency and sourness. Estimated phenotypic and genetic correlations between all pairs of traits are listed in Supplemental Table S1.

    Figure 1.  Transverse cross-sections of apple slices showing range in (a) flesh colouration, and (b) flesh browning disorder for Type 1 red-fleshed apple. The weighted cortical intensity (WCI) scores (0−9 scale) and the proportion of the cortex area showing symptoms of flesh browning disorder are also displayed.

    A few seedlings had no red pigment in the cortex, but the average WCI score across all seedlings was 2.25. About two-thirds of the seedlings did not display any FBD symptoms, but among the remaining seedlings, FBD ranged between 1% and 58% (Fig. 2). The average WCI score for the 'low' and 'high' WCI pool was 0.45 and 5.2, respectively, while the average FBD was 0% and 20.6% for the 'low' and 'high' FBD pool, respectively (Supplemental Table S2). The average WCI score of the seedlings in the FBD pools was similar (Low: 5.5; High: 4.6), while the average FBD of the high- and low-WCI pools was 4.5% and 0.1%, respectively.

    Figure 2.  The distribution of weighted cortex intensity (WCI) scores (a) and the internal flesh browning disorder IFBD%; (b) in the population of ~900 apple seedlings. The green and red circles highlight the individuals used to form the 'low' and 'high' pools of samples.

    After filtering, about 204,000 SNPs were used and the average sequencing depth of SNP loci was similar for the two pools (42 vs 44). There was a near-perfect correlation between the Z-test statistics and G-statistics, so only the latter are discussed hereafter. A plot of the G' values, smoothed over 2 Mb windows, is shown for all 17 chromosomes (Chrs) in Supplemental Fig. S2. XP-GWAS identified genomic regions significantly associated with FBD on 12 out of the 17 Chrs (Fig. 3), and putative candidate genes within ±1.0 Mb distanceof the significant G' peaks were identified (Table 1). Additional genomic regions, which did not meet the significance threshold but displayed distinguished G' peaks, were also identified across all chromosomes (Supplemental Table S3).

    Figure 3.  G'-statistics across the linkage groups (LG) showing significant association with the flesh browning disorder (FBD) in apple The horizontal red lines indicate the significance threshold. The putative candidate genes (refer to Table 1) underpinning various G' peaks are also shown.
    Table 1.  A list of the genomic regions associated with internal flesh browning disorder (FBD) in apples. Putative candidate genes residing within these regions are also listed using GDDH13v1.1 reference genome assembly.
    ChrGenomic region (Mb)Putative genes functions
    221.9–23.5Ethylene-responsive element binding factor 13 (MdERF13: MD02G1213600);
    33.1–4.9cinnamate 4-hydroxylase (C4H) enzyme (MD03G1051100, MD03G1050900 and MD03G1051000); MdWRKY2: MD03G1044400; MdWRKY33 (MD03G1057400)
    38.2–9.6ascorbate peroxidase 1 (MdAPX1: MD03G1108200, MD03G1108300)
    314.7–16.6senescence-related MdNAC90 (MD03G1148500)
    336.5–37.5Ethylene response sensor 1 (MdERS1: MD03G1292200); flavonoid biosynthesis protein MdMYB12 (MD03G1297100); Heat shock protein DnaJ (MD03G1296600, MD03G1297000); pectin methylesterase (MdPME) inhibitor protein (MD03G1290800, MD03G1290900, MD03G1291000).
    411.0–13.0phenylalanine and lignin biosynthesis protein MdMYB85 (MD04G1080600)
    422.1–24.4MYB domain protein 1 (MD04G1142200); HSP20-like protein (MD04G1140600); UDP-glucosyltransferase (UGT) proteins UGT85A7 (MD04G1140700, MD04G1140900); UGT85A3 (MD04G1140800); UGT (MD04G1141000, MD04G1141300); UGT85A2 (MD04G1141400); UGT85A4 (MD04G1141500); DNAJ heat shock protein (MD04G1153800, MD04G1153900, MD04G1154100)
    427.6–29.6HCT/HQT regulatory genes MD04G1188000 and MD04G1188400
    629.8–31.7Volz et al. (2013) QTL for IFBD; anthocyanin regulatory proteins MdMYB86 (MD06G1167200); triterpene biosynthesis transcription factor MdMYB66 (MD06G1174200); Cytochrome P450 (MD06G1162600; MD06G1162700, MD06G1162800, MD06G1163100, MD06G1163300, MD06G1163400, MD06G1163500, MD06G1163600, MD06G1163800, MD06G1164000, MD06G1164100, MD06G1164300, MD06G1164400, MD06G1164500, MD06G1164700)
    713.9–15.9heat shock protein 70B (MD07G1116300)
    717.4–19.0Drought-stress WRKY DNA-binding proteins (MdWRKY56: MD07G1131000, MD07G1131400)
    723.6–25.2MdPAL2 (MD07G1172700); drought-stress gene NGA1 (MD07G1162400); DNAJ heat shock family protein (MD07G1162300, MD07G1162200), stress-response protein (MdNAC69: MD07G1163700, MD07G1164000)
    97.9–11.8MD09G1110500 involved in ascorbate oxidase (AO); MdUGT proteins (MD09G1141200, MD09G1141300, MD09G1141500, MD09G1141600, MD09G1141700, MD09G1141800, MD09G1142000, MD09G1142500, MD09G1142600, MD09G1142800, MD09G1142900, MD09G1143000, MD09G1143200, MD09G1143400) involved in flavonoids biosynthesis; heat shock proteins 89.1 (MD09G1122200) and HSP70 (MD09G1137300); Ethylene-forming enzyme MD09G1114800; Anthocyanin regulatory protein MdNAC42 (MD09G1147500, MD09G1147600)
    914.9–16.5triterpene biosynthesis transcription factor protein MdMYB66 (MD09G1183800);
    111.7–3.0ethylene response factor proteins (MdEIN-like 3: MD11G1022400)
    1138.6–40.6Senescence-related gene 1 (MD11G1271400, MD11G1272300, MD11G1272000, MD11G1272100, MD11G1272300, MD11G1272400, and MD11G1272500); Chalcone-flavanone isomerase (CHI) protein (MD11G1273600) and MdbHLH3 (MD11G1286900, MDP0000225680); cytochrome P450 enzyme (MD11G1274000, MD11G1274100, MD11G1274200, MD11G1274300, MD11G1274500, and MD11G1274600); heat shock transcription factor A6B (MdHSFA6B: MD11G1278900) – involved in ABA-mediated heat response and flavonoid biosynthesis.
    122.2–3.5heat shock protein 70-1 ((MD12G1025600, MD12G1025700 and MD12G1026300) and heat shock protein 70 (MD12G1025800 and MD12G1025900 and MD12G1026000); ethylene (MD12G1032000) and auxin-responsive (MD12G1027600) proteins.
    127.2–8.4DNAJ heat shock domain-containing protein (MD12G1065200 and MD12G1067400)
    1327.5–30.5MD13G1257800 involved in polyphenol 4-coumarate:CoA ligase (4CL)
    1337.5–39.5pectin methyl esterase inhibitor superfamily protein MdPMEI (MD13G1278600)
    1425.4–27.5Drought-stress response gene MdWRKY45 (MD14G1154500); chalcone synthase (CHS) family proteins (MD14G1160800 and MD14G1160900); triterpene biosynthesis transcription factor MdMYB66 (MD14G1180700, MD14G1181000, MD14G1180900); anthocyanin biosynthesis protein (MdMYB86: MD14G1172900); cytochrome P450 proteins (MD14G1169000, MD14G1169200, MD14G1169600, MD14G1169700)
    150–1.5Ethylene synthesis proteins (MD15G1020100, MD15G1020300 and MD15G1020500); dihydroflavonol reductase (DFR) gene (MD15G1024100)
    154.6–6.8MdMYB73 (MD15G1076600, MD15G1088000) modulates malate transportation/accumulation via interaction with MdMYB1/10; MdNAC52: MD15G1079400) regulates anthocyanin/PA; heat shock transcription factor B4 (MD15G1080700); stress-response protein MdWRKY7 (MD15G1078200)
    1553.5–54.9MdC3H (MD15G1436500) involved in chlorogenic acid biosynthesis; MdEIN3 (MD15G1441000) involved in regulating ethylene synthesis and anthocyanin accumulation
    168.9–10.9SAUR-like auxin-responsive protein (MD16G1124300) and MdNAC83: (MD16G1125800) associated with fruit ripening; MD16G1140800 regulates proanthocyanidin; MdPAE genes (MD16G1132100, MD16G1140500) regulates ethylene production.
     | Show Table
    DownLoad: CSV

    The ethylene-responsive factor 13 (MdERF13: MD02G1213600) resided within the significant region (21.9–23.5 Mb) on Chr2, while the ascorbate peroxidase 3 (MdAPX3: MD02G1127800, MD02G1132200) resided in the prominent region between 9.45 and 11.28 Mb) (Fig. 3, Table 1). There were several genomic regions showing association with FBD on Chr3. The first region (8.2–9.5 Mb) flanked MdAPX1 (MD03G1108200, MD03G1108300), while MdERF3 (MD03G1194300) and heat shock protein 70 (HSP70: MD03G1201800, D03G1201700) resided within the prominent peak region (25.5–27.5 Mb). Another significant region (36.5–37.5 Mb) at the bottom on Chr3 harboured ethylene response sensor 1 (MdERS1: MD03G1292200), a flavonoid-biosynthesis related protein MdMYB12 (MD03G1297100), HSP DnaJ and pectin methyl esterase (PME) inhibitor proteins (Fig. 3, Table 1).

    The significant G' region (27.6–29.6 Mb) on Chr4 flanked the genes MD04G1188000 and MD04G1188400 involved in the biosynthesis of hydroxycinnamoyl CoA shikimate/quinate hydroxycinnamoyl transferase (HCT/HQT), while the adjacent region (22.1–24.3 Mb) harboured a cluster of UDP-glucosyltransferase (UGT) proteins and HSP (Table 1). Another region between 11.0 and 13.0 Mb encompassed phenylalanine and lignin biosynthesis gene MYB85 (MD04G1080600)[17]. The only significant region associated with FBD on Chr6 spanned between 29.8 and 31.7 Mb, which included a SNP earlier reported associated with FBD[4]. This region also harboured several MYB proteins (MdMYB86: MD06G1167200; MdMYB98: MD06G1172900; MdMYB66: MD06G1174200) and a large cluster of cytochrome P450 proteins (Table 1).

    A significant FBD-associated region (13.9–15.9 Mb) on Chr7 encompassed the HSP 70B (MD07G1116300), while another significant region between 23.6 Mb and 25.2 Mb flanked the gene coding for phenylalanine ammonia lyase 2 (MdPAL2: MD07G1172700), a drought stress gene NGA1 (MD07G1162400), DNAJ HSP, and stress-response protein (MdNAC69: MD07G1163700, MD07G1164000) (Fig. 3, Table 1). A large significant region spanning between 7.9 and 11.8 Mb on Chr9 encompassed the gene MD09G1110500 putatively involved in ascorbate oxidase (AO), HSP (HSP70: MD09G1137300; HSP89.1: MD09G1122200), an ethylene-forming enzyme (MD09G1114800), and MdNAC42 (MD09G1147500, MD09G1147600). Another significant genomic region on Chr9 was between 14.9 and 16.5 Mb, which harbours the MYB domain protein MdMYB66 (MD09G1183800) (Table 1).

    A sharp G' peak region (1.7–3.0 Mb) on Chr11 associated with FBD encompassed ethylene insensitive 3 (MdEIN3: MD11G1022400) along with a cluster of UGT proteins, WD-40, and bHLHL proteins (Table 1). A significant region between 38.6 and 40.6 Mb at the bottom of Chr11 was dominated by clusters of senescence-related genes and cytochrome P450 enzymes. This genomic region also flanked a chalcone-flavanone isomerase (CHI) family protein (MD11G1273600) and a bHLH protein (MdbHLH3: MD11G1286900, MDP0000225680), along with the heat shock transcription factor A6B (HSFA6B: MD11G1278900) (Table 1, Supplemental Table S3).

    The region (2.2–3.3 Mb) associated with FBD on Chr12 flanked genes for the HSP 70 and 70-1, NAC proteins, ethylene and auxin-responsive proteins (Table 1). An adjacent significant region (7.2–8.4 Mb) flanked DNAJ HSP, along with bZIP, bHLH and WD-40 repeat-like proteins (Table 1; Supplemental Table S3). There was a large genomic region on Chr13 showing a significant association with FBD. In this region, the first G' peak (27.5–30.5 Mb) harboured MD13G1257800, which regulates polyphenol 4-coumarate: CoA ligase (4CL) synthesis. The second G' peak region (32.6–34.7 Mb) corresponded to an earlier mapped QTL for flesh browning in white-fleshed apples[18].

    The significant region (25.4–27.5 Mb) on Chr14 harboured various genes for proteins with different putative functions, such as AP2 proteins, WD-40 repeat family proteins, bHLH proteins, MdNAC83 (MD14G1150900), MdWRKY45 (MD14G1154500), chalcone synthase (CHS) family proteins (MD14G1160800 and MD14G1160900), cytochrome P450 proteins, and several MYB domain proteins (MdMYB86: MD14G1172900; MdMYB98: MD14G1179000; MdMYB66: MD14G1180700, MD14G1181000, MD14G1180900) (Fig. 3, Table 1, Supplemental Table S3). A sharp G' peak (15.2–17.2 Mb) on Chr14 did not meet the significance threshold corresponding to the FBD QTL in white-fleshed apples[18].

    The upper 1.5 Mb region on Chr15 associated with FBD encompassed several transcription factor families, including WD-40 repeats, bHLH, bZIP, ethylene synthesis proteins, and dihydroflavonol reductase (DFR) protein (MD15G1024100) (Fig. 3, Table 1; Supplemental Table S3). Another significant region on Chr15 (4.6−6.8 Mb) harboured HSF B4, MdMYB73 (which interacts with MdMYB1/10 to modulate malate transportation) and MdNAC52 (MD15G1079400), which regulates anthocyanin and PA synthesis by directly regulating MdLAR[19]. A significantly associated region at the bottom of Chr15 (53.5–54.9 Mb) harboured MdNAC35 (MD15G1444700), MdC3H (MD15G1436500) involved in the production of p-coumarate 3-hydroxylase (C3H) enzyme, which plays a role in chlorogenic acid biosynthesis, and MdEIN3 (MD15G1441000).

    The significant region between 8.9 and 10.9 Mb on Chr16 harboured a gene for SAUR-like auxin-responsive protein (MD16G1124300) and MdNAC83 (MD16G1125800), both of which have been reportedly associated with apple fruit ripening[20]. This region also encompassed MD16G1140800, which regulates PA[11], and a gene for the pectin acetyl esterase protein MdPAE10 (MD16G1132100) involved in ethylene production and shelf-life[21]. A sharp G' peak region (20.4–22.4 Mb) in the middle of Chr16 flanked MdERF1B (MD16G1216900) and MdMYB15 (MD16G1218000 and MD16G1218900), involved in altering anthocyanin and PA concentrations[10] (Fig. 3, Table 1, Supplemental Table S3).

    The average sequencing depth of the SNP loci (~160,000) retained for marker-trait association was similar for the two WCI pools (41 vs 44). The significant regions were located on Chrs 2, 4, 6, 7, 10, 15 and 16 (Fig. 4). The genomic intervals within ±1.0 Mb of the significant G' peaks, and the putative candidate genes within those intervals, are listed in Table 2. Additional genomic regions, which did not meet the significance threshold but displayed distinct G' peaks, were also identified across most chromosomes ( Supplemental Fig. S2, Supplemental Table S3). A significant region on Chr4 encompassed chalcone synthase (CHS) genes, along with an ERF (MD04G1009000) involved in regulating PA biosynthesis[11].

    Figure 4.  G’-statistics across the linkage groups (LG) showing significant association with the weighted cortex intensity (WCI) in apple. The horizontal red lines indicate the significance threshold. The putative candidate genes (refer to Table 2) underpinning various G' peaks are also shown.
    Table 2.  A list of the genomic regions significantly associated with the weighted cortex intensity (WCI) in apples. Putative candidate genes residing within these regions are also listed using GDDH13v1.1 reference genome assembly.
    ChrGenomic region (Mb)Putative genes functions
    211.2–13.2UDP-glucosyltransferase (UGT) proteins (MD02G1153000, MD02G1153100, MD02G1153200, MD02G1153300; MD02G1153400; MD02G1153500; MD02G1153700; MD02G1153800; MD02G1153900);
    40–1.2Chalcone synthase (CHS) genes (MD04G1003000; MD04G1003300 and MD04G1003400); DNAJ heat shock protein (MD04G1003500); MdERF (MD04G1009000) involved in regulating PA biosynthesis.
    69.5–11.0Ubiquitin protein (MD06G1061100); stress-response WRKY protein MdWRKY21 (MD06G1062800);
    612.1–13.6pectin methylesterase inhibitor superfamily protein (MdPME: MD06G1064700); phenylalanine and lignin biosynthesis gene (MdMYB85)
    616.0–17.6MD06G1071600 (MDP0000360447) involved in leucoanthocyanidin dioxygenase (LDOX) synthesis
    624.0–26.0UDP-glycosyltransferase proteins (MD06G1103300, MD06G1103400, MD06G1103500 and MD06G1103600); auxin response factor 9 (MdARF9: MD06G1111100) and heat shock protein 70 (Hsp 70; MD06G1113000).
    74.1–5.6Ethylene insensitive 3 protein (MdEIN3: MD07G1053500 and MD07G1053800) involved in proanthocyanidins (PA) biosynthesis; ubiquitin-specific protease (MD07G1051000, MD07G1051100, MD07G1051200, MD07G1051500, MD07G1051300, MD07G1051700 and MD07G1051800).
    1015.9–18.3bHLH proteins (MD10G1098900, MD10G1104300, and MD10G1104600); UGT protein (MD10G1101200); UGT 74D1 (MD10G1110800), UGT 74F1 (MD10G1111100), UGT 74F2 (MD10G1111000).
    1035.6–38.6Stress-response WRKY proteins (MdWRKY28: MD10G1266400; MdWRKY65: MD10G1275800). ethylene responsive factors (MdERF2: MD10G1286300; MdERF4: MD10G1290400, and MdERF12: MD10G1290900) and a NAC domain protein (MdNAC73: MD10G1288300); polyphenol oxidase (PPO) genes ((MD10G1298200; MD10G1298300; MD10G1298400; MD10G1298500; MD10G1298700; MD10G1299100; MD10G1299300; MD10G1299400).
    1524.8–26.8heat shock factor 4 (MD15G1283700), drought-stress response WRKY protein 7 (MdWRKY7: MD15G1287300), MdMYB73 (MD15G1288600) involved in ubiquitination and malate synthesis
    1531.7–34.2MYB domain protein 93 (MdMYB93: MD15G1323500) regulates flavonoids and suberin accumulation (Legay et al. 2016); ubiquitin -specific protease 3 (MD15G1318500),
    161.5–3.4MYB domain proteins (MD16G1029400) regulates anthocyanin; senescence-associated gene 12 (MD16G1031600); ethylene response factor proteins (MdERF118: MD16G1043500, and MD16G1047700 (MdRAV1); MdMYB62 ( MD16G1040800) flavonol regulation; malate transporter MdMa2 (MD16G1045000: MDP0000244249), Ubiquitin-like superfamily protein (MD16G1036000),
    165.4–7.4MdMYB88 (MD16G1076100) regulates phenylpropanoid synthesis and ABA-mediated anthocyanin biosynthesis; MdMYB66 (MD16G1093200) regulates triterpene biosynthesis, MdWRKY72 (MD16G1077700) mediates ultraviolet B-induced anthocyanin synthesis.
     | Show Table
    DownLoad: CSV

    There were distinct G' peaks within a large genomic region (spanning between 9.0 and 18.0 Mb) significantly associated with WCI on Chr6 (Fig. 4). The G' peak region 12.1–13.6 Mb encompassed MdMYB85 (MD06G1064300), and a gene for a pectin methyl esterase

    (PME) inhibitor protein (MD06G1064700); while the region between 16.0 and 17.6 Mb flanked the genes involved in leucoanthocyanidin dioxygenase (LDOX) synthesis (Table 2). The auxin response factor 9 (MdARF9: MD06G1111100) and the HSP70 (MD06G1113000) resided in the significant region between 24.0 and 26.0 Mb on Chr6 (Table 2).

    The genomic region (4.1–5.6 Mb) with significant association with WCI on Chr7 flanked MdEIN3 (MD07G1053500, MD07G1053800), which plays an important role in the regulatory network of PA biosynthesis[11]. A significant region (35.6–38.7 Mb) on Chr10 encompassed gene clusters for bHLH and WRKY proteins along with an ethylene repressor factor (MdERF2: MD10G1286300) (Fig. 4, Table 2; Supplemental Table S3). This region also harboured MdERF4 (MD10G1290400), MdERF12 (MD10G1290900) and NAC domain genes (MdNAC73: MD10G1288300), which have been reported to be associated with fruit ripening[20]. In addition, there was a cluster of polyphenol oxidase (PPO) genes residing in this region (Table 2).

    On Chr15, the significant WCI-associated region spanning between 31.7 Mb and 34.2 Mb encompassed several TF families, including MdMYB93 (MD15G1323500) (Fig. 4, Table 2, Supplemental Table S3). A distinguished G' peak region (24.9–26.9 Mb) on Chr15 harboured genes for WD-40 repeat-like proteins, bHLH proteins, redox responsive transcription factor 1 (MD15G1283200), cytochrome P450 proteins, HSF4 (MD15G1283700), and MdWRKY7 (MD15G1287300). The MdMYB73 (MD15G1288600) residing in this region has been shown to interact with MdMYB1/10 and regulates several functions, including cold-stress response, ubiquitination and malate synthesis[22, 23].

    The significant 1.5–3.4 Mb region on Chr16 flanked an anthocyanin repressor MYB protein (MD16G1029400), senescence-associated gene 12 (MD16G1031600), malate transporter MdMa2 (MD16G1045000: MDP0000244249), and two ethylene response factor genes (MdERF118: MD16G1043500; MdRAV1: MD16G1047700), which interact in retaining flesh firmness[24]. The MdMYB62 (MD16G1040800) gene residing in this region is phylogenetically linked to MdMYB8 (MD06G1217200), which plays a major role in flavonoid biosynthesis[25]. Another significant region (5.4–7.4 Mb) on Chr16 harboured several bHLH genes, including MdMYB88 (MD16G1076100) involved in phenylpropanoid synthesis resulting in drought resistance[26] and ABA-mediated anthocyanin production[27]. The MdMYB66 (MD16G1093200) gene, which is involved in triterpene biosynthesis, and the MdWRKY72 (MD16G1077700) gene involved in anthocyanin synthesis, were also present in this region (Table 2, Fig. 4).

    The genomic regions tagged by XP-GWAS were enriched for regulatory functions, and several of these have a paralog (Table 3). For example, paralogues for a trio of genes (MdMYB86, MdMYB98 and Cytochrome 450) resided in the FBD-associated genomic regions on Chr6 and Chr14. Several pairs of genes resided together in the FBD-associated paralog regions; for example, MdNAC90 and germin-like protein 10 on Chr3 and Chr11; WRKY55 and WRKY70 on Chr1 and Chr7; and ERF1 and ERF5 on Chr4 and Chr6. The ethylene response factor MdERF1B had a paralog in the FBD-associated region on Chr13 and Chr16. Interestingly, three copies of MdNAC83 (Chrs 14, 16 17) and four copies of MdMYB66 (Chrs 4, 6, 9 and 14) resided in the FBD-associated regions. A pair of genes (RAV1 and MYB62) resided together in the WCI-associated paralog regions on Chrs13 and 16 (Table 3).

    Table 3.  A list of homologues genes/genomic regions associated with the internal flesh browning disorder (IFBD) and red flesh (WCI) in apples. Putative candidate genes residing within these regions are also listed using GDDH13v1.1 reference genome assembly.
    TraitChr (genomic region: Mb)Gene namePredicted gene IDPutative function
    IFBDChr6 (29.8–31.7 Mb)MdMYB66MD06G1174200Suberin/triterpene deposition
    IFBDChr14 (25.4–27.5 Mb)MdMYB66MD14G1180700, MD14G1181000, MD14G1180900Suberin/triterpene deposition
    IFBDChr4 (7.0–8.1 Mb)MdMYB66MD04G1060200Suberin/triterpene deposition
    IFBDChr9 (14.9–16.5 Mb)MdMYB66MD09G1183800Suberin/triterpene deposition
    IFBDChr6 (29.8–31.7 Mb)MdMYB86MD06G1167200Anthocyanin regulation
    IFBDChr14 (25.4–27.5 Mb)MdMYB86MD14G1172900Anthocyanin regulation
    IFBDChr6 (29.8–31.7 Mb)MdMYB98MD06G1172900Drought stress response
    IFBDChr14 (25.4–27.5 Mb)MdMYB98MD14G1179000Drought stress response
    IFBDChr6 (29.8–31.7 Mb)Cytochrome P450MD06G1162600, MD06G1162700, MD06G1162800, MD06G1163100, MD06G1163300, MD06G1163400, MD06G1163500, MD06G1163600, MD06G1163800, MD06G1164000, MD06G1164100, MD06G1164300, MD06G1164400, MD06G1164500, MD06G1164700Flavonoid and triterpenic metabolism
    IFBDChr14 (25.4–27.5 Mb)Cytochrome P450MD14G1169000, MD14G1169200, MD14G1169600, MD14G1169700Flavonoid and triterpenic metabolism
    IFBDChr3 (14.7–16.6 Mb)MdNAC90MD03G1148500Senescence-related
    IFBDChr11 (16.5–18.5 Mb)MdNAC90MD11G11679000Senescence-related
    IFBDChr14 (25.4–27.5 Mb)MdNAC83MD14G1150900Senescence/ripening-related
    IFBDChr16 (8.9–10.9 Mb)MdNAC83MD16G1125800Senescence/ripening-related
    IFBDChr17 (0–0. 7 Mb)MdNAC83MD17G1010300Senescence/ripening-related
    IFBDChr3 (14.7–16.6 Mb)Germin-like protein 10MD03G1148000Polyphenol oxidase
    IFBDChr11 (16.5–18.5 Mb)Germin-like protein 10MD11G1167000, MD11G1167100, MD11G1167400, MD11G1169200Polyphenol oxidase
    IFBDChr1 (26.9–28.5 Mb)MdWRKY55MD01G1168500Drought stress response
    IFBDChr7 (30.3–31.9 Mb)MdWRKY55MD07G1234600Drought stress response
    IFBDChr1 (26.9–28.5 Mb)MdWRKY70MD01G1168600Drought stress response
    IFBDChr7 (30.3–31.9 Mb)MdWRKY70MD07G1234700Drought stress response
    IFBDChr4 (7.0–8.1 Mb)MdERF1MD04G1058000Ethylene responsive factor
    IFBDChr6 (5.5–7.3 Mb)MdERF1MD06G1051800Ethylene responsive factor
    IFBDChr4 (7.0–8.1 Mb)MdERF5MD04G1058200Ethylene responsive factor
    IFBDChr6 (5.5–7.3 Mb)MdERF5MD06G1051900Ethylene responsive factor
    IFBDChr13 (18.0–20.0 Mb)MdERF1BMD13G1213100Ethylene response factor 1
    IFBDChr16 (20.4–22.4 Mb)MdERF1BMD16G1216900Ethylene response factor 1
    WCIChr13 (2.6–4.4 Mb)MdRAV1MD13G1046100Ethylene responsive factor
    WCIChr16 (1.5–3.4 Mb)MdRAV1MD16G1047700Ethylene responsive factor
    WCIChr13 (2.6–4.4 Mb)MdMYB62MD13G1039900Flavonol biosynthesis
    WCIChr16 (1.5–3.4 Mb)MdMYB62MD16G1040800Flavonol biosynthesis
    IFBDChr9 (7.9–11.8 Mb)MdNAC42MD09G1147500, MD09G1147600Anthocyanin accumulation
    WCIChr17 (11.4–12.4 Mb)MdNAC42MD17G1134400Anthocyanin accumulation
    IFBDChr9 (7.9–11.8 Mb)HSP 70MD09G1137300Heat stress response
    WCIChr17 (11.4–12.4 Mb)HSP 70MD17G1127600Heat stress response
    IFBDChr4 (11.0–13.0 Mb)MdMYB85MD04G1080600Phenylalanine and lignin biosynthesis
    WCIChr6 (12.1–13.6 Mb)MdMYB85MD06G1064300Phenylalanine and lignin biosynthesis
    IFBDChr15 (13.2–14.2 Mb)MdEBF1MD15G1171800Ethylene inhibition
    WCIChr8 (15.2–17.2 Mb)MdEBF1MD08G1150200Ethylene inhibition
    IFBDChr15 (53.5–54.9 Mb)MdEIN3MD15G1441000Ethylene insensitive 3 protein
    WCIChr8 (30.3–31.6 Mb)MdEIN3MD08G1245800Ethylene insensitive 3 protein
    IFBDChr11 (1.7–3.0 Mb)MdEIN3MD11G1022400Ethylene insensitive 3 protein
    WCIChr7 (4.1–5.6 Mb)MdEIN3MD07G1053500, MD07G1053800Ethylene insensitive 3 protein
    IFBDChr15 (4.6–6.8 Mb)MdMYB73MD15G1076600, MD15G1088000Cold-stress response & malate accumulation
    WCIChr15 (24.8–26.8 Mb)MdMYB73MD15G1288600Cold-stress response & malate accumulation
    IFBDChr15 (4.6–6.8M b)MdWRKY7MD15G1078200Anthocyanin accumulation
    WCIChr15 (24.8–26.8 Mb)MdWRKY7MD15G1287300Anthocyanin accumulation
    IFBDChr15 (43.5–45.5 Mb)MdMYB93MD15G1369700Flavonoid & suberin accumulation
    WCIChr15 (31.7–34.2 Mb)MdMYB93MD15G1323500Flavonoid & suberin accumulation
     | Show Table
    DownLoad: CSV

    Paralogs of several regulatory functions were also found in the regions associated with either FBD or WCI; for example, a significant region (7.9–11.8 Mb) harbouring a gene trio (MdNAC42, HSP70 and HSP89.1) on Chr9 was associated with FBD, but the paralogs of this trio also resided in a distinct G' region associated with WCI on Chr17 (Table 3, Supplemental Table S3). Some other examples included MdMYB85 (Chr4 for FBD, and Chr6 for WCI), MdEBF1 (Chr8 for WCI, and Chr15 for FBD), MdEIN3 (Chr8 for WCI, and Chr15 for FBD), and EIN3 (Chr7 for WCI, and Chr11 for FBD). A pair of genes (MdMYB73 and MdWRKY7) resided together in the separate regions associated with FBD (4.6–6.8 Mb) and WCI (24.8–26.8 Mb) within Chr15 (Table 3).

    The red pigmentation in fruit flesh differed amongst transgenic lines, with fruit from line A10 presenting the most deeply pigmented tissues (Supplemental Fig. S3), while those from lines A2 and A4 were similar in having a lower intensity of pigmentation. No pigmentation was observed in the flesh of control 'Royal Gala' (RG) fruit. RNAseq analysis of a representative transgenic line (A2) compared with RG revealed that a total of 1,379 genes were differentially expressed (log 2-fold), with 658 genes upregulated and 721 genes downregulated. This list was then assessed for commonality with the genomic regions and candidate genes from the XP-GWAS.

    Genes that contained mis-sense SNPs (which would result in a change in predicted protein) and that were also differentially expressed between A2 red-fleshed transgenic line and white-fleshed 'Royal Gala' (control) apples included anthocyanin-related flavanone 3-hydroxylase, chalcone synthase and dihydroflavonol 4-reductase (Table 2 & 4). Many more upstream DNA variants (e.g. in potential promoter-controlling elements) were seen in this group of differentially expressed genes (DEGs) that were also in regions underlying WCI or FBD. Genes encoding enzymes that may be involved in FBD, such as a Rho GTPase activating protein, peroxidase, lipoxygenase 1, and ethylene-forming enzyme (ACO4), were DEGs and showed mis-sense SNPs, including a potential stop codon in the Rho GTPase-activating protein (Table 4). This intersection between DNA change and differential expression warrants further research to evaluate the functions of these genes.

    Table 4.  List of candidate genes associated with WCI or FBD in the XP-GWA and R6:MdMYB10 apple datasets.
    Mutations in GWAS apple populationPredicted gene functionExpression in R6 and ‘Royal Gala’ apples
    GeneMissense
    SNPs
    SNP
    Stop
    Upstream
    variants
    DatasetLocusAnnotation and TAIR IDAverage
    RPKM
    R6 flesh
    Average RPKM
    WT flesh
    log2Fold
    Change
    MD02G11322001FBD Supplemental Table S3Chr02:9450615-11289014flavanone 3-hydroxylase
    (F3H, TT6, F3'H) AT3G51240
    387752.47
    MD02G11336003FBD Supplemental Table S3Chr02:9450615-11289014fatty acid desaturase 5
    (FAD5) AT3G15850
    18109.98
    MD02G11537001WCI Table 2Chr02:11278254-13234556UDP-Glycosyltransferase, lignin related AT2G1856010385341.03
    MD03G10592002FBD Supplemental Table S3Chr03:3177994-4910866Peroxidase AT5G053402514.21
    MD03G1143300422FBD Table 1Chr03:14797661-16611554bZIP transcription factor (DPBF2, AtbZIP67) AT3G44460114852-2.80
    MD03G1147700121FBD Table 1Chr03:14797661-16611554Rho GTPase activating protein AT5G22400205891.29
    MD04G1003400331WCI Table 2Chr04:1-1284894Chalcone synthase (CHS, TT4) AT5G1393014433122.34
    MD04G120410063FBD Table 1Chr04:27629536-29612282lipoxygenase 1 (LOX1, ATLOX1) AT1G550205462681.09
    MD06G1160700129FBD Table 1Chr06:29862341-31737341peptide met sulfoxide reductase AT4G251302121254562.04
    MD06G116140014FBD Table 1Chr06:29862341-31737341Pectin lyase-like protein AT5G63180614129530-2.22
    MD07G1240700126FBD Supplemental Table S3Chr07:30357332-31979025Fe superoxide dismutase 2 AT5G511001051826-3.99
    MD07G13069006FBD Supplemental Table S3Chr07:34607570-36531467UDP-glucosyl transferase 78D2 AT5G170506771112.78
    MD08G12491002WCI Supplemental Table S3Chr08:30486569-31607516HSP20-like chaperone (ATHSP22.0) AT4G102503186456.12
    MD09G111400032FBD Table 1Chr09:7999212-11883202fatty acid desaturase 5 (FAD5) AT3G158507608.68
    MD09G11468001FBD Table 1Chr09:7999212-11883202PHYTOENE SYNTHASE (PSY) AT5G1723038456160901.32
    MD10G132810045FBD Supplemental Table S3Chr10:40235258-41736791ethylene-forming enzyme (ACO4) AT1G050108769033910331.21
    MD15G10236001FBD Table 1Chr15:1-1487288jasmonic acid carboxyl methyltransferase (JMT) AT1G19640472111552.07
    MD15G10241008FBD Table 1Chr15:1-1487288dihydroflavonol 4-reductase (DFR, TT3, M318) AT5G428008723091.63
    MD17G11334004WCI Supplemental Table S3Chr17:11422073-12433463PHYTOENE SYNTHASE (PSY) AT5G172306001721.87
    MD17G12606002FBD Supplemental Table S3Chr17:31098955-32776972dehydroascorbate reductase 1 (DHAR3) AT5G167101867−1.78
     | Show Table
    DownLoad: CSV

    Several gene families (e.g. chalcone synthase, UGT, anthocyanin synthase, HSP, PAL, ERF, WRKY proteins, ABA, and bZIP transcription factors) residing in WCI-associated genomic regions have been reported to be associated with RF in apple[7, 12, 28]. Several of these genes are implicated in stress responses, suggesting that flavonoids and anthocyanin biosynthesis could also be associated with stress (e.g., drought, water loss) tolerance of RF apples[7].

    MdMYB73, which may play a role in the cold-stress response[23, 29], resided in the WCI-associated regions. Ethylene is among the various modulators of environmental stresses induced by factors such as drought and cold temperatures[9, 28]. Several ERF proteins (e.g. MdERF1B, MdERF3) interact with the promoters of MYB domain proteins to regulate anthocyanin and proanthocyanidin (PA) accumulation in apple[10, 11, 30]. The WCI-associated genomic regions in our study flanked several ERFs previously reported to have roles in anthocyanin and PA regulation. MdERF4 (MD10G1290400), one of the three ERFs residing in the WCI-associated region on Chr10, has high phylogenetic similarity with MdERF38, which interacts with MdMYB1 to regulate drought-related anthocyanin biosynthesis[31]. Some ERFs (e.g. MdERF2: MD10G1286300; MdERF4: MD10G1290400, and MdERF12: MD10G1290900) and MdNAC73 residing in the WCI-associated regions have been reported to be associated with fruit maturation[20] – suggesting an interaction between ethylene production and RF colour[12].

    The WCI-associated region on Chr15 included a gene MdMYB73 (MD15G1288600) involved in malate acid synthesis. Previous studies have reported that MdMYB1 regulates both anthocyanin and malate accumulation, and perhaps MdBT2 regulates MdMYB73-mediated anthocyanin accumulation[22, 23]. The malate transporter gene MdMa2 and the MdMYB7 gene, involved in the regulation of anthocyanin and flavanols in RF apple[32], resided together in the WCI-associated upper region of LG16. Taken together, these results support the hypothesis that there is interplay between anthocyanin and malate accumulation in the RF apple. There was a cluster of ubiquitin-specific proteases underpinning the WCI-associated regions on Chrs 7 and 15, which is supported by earlier reports suggesting that ubiquitin-specific proteases respond to auxin and might suppress anthocyanin biosynthesis proteins[33]. The auxin response factor 9 (MdARF9: MD06G1111100), which has also been shown to suppress anthocyanin biosynthesis in RF callus samples[33], also resided in the WCI-associated genomic region on LG7.

    Seedlings in both the low- and high-FBD pools carried the MdMYB10 gene, which suggests that MdMYB10 itself is not the causal factor of FBD in RF apples. Long-term cold storage generally results in senescence-related flesh breakdown, and several transcription factor genes (e.g. MYB, WRKY, NAC, ERF, cytochrome P450, and HSP) have been shown to express differentially during long-term cold storage[9]. The FBD-associated genomic regions in our study harboured several ERFs, suggesting that ethylene synthesis proteins may be contributing to cell wall disassembly, allowing PPO enzymes to come into contact with phenolic compounds and potentially leading to FBD symptoms[18].

    Pectin methyl esterase (PME) genes resided in the FBD-associated significant regions on Chrs 3, 13 and 17. Volatile generation and senescence degradation have been suggested to be bio-markers of FBD, and the expression levels of methyl esters were found to be associated with FBD and senescence in 'Fuji' apples after cold storage[34]. The MdPME2 gene has also been reported to be associated with apple flesh firmness and mealiness[21, 35]. The co-location of genes encoding cell wall-degrading enzymes (MdPME) and QTLs for FBD has been reported in apple[18, 36].

    The clusters of cytochrome P450 enzymes and senescence-related genes resided in some of the FBD-associated regions in this study. There are no earlier reports of the involvement of P450 enzymes in apple FBD expression, but some genes related to cytochrome P450 were found to be upregulated during litchi fruit senescence[37]. Pericarp browning in litchi is mainly attributable to the degradation of anthocyanin, and the ABA-initiated oxidation of phenolic compounds by PPO[37]. Flavonoids are among the major polyphenols in RF apples[5], and cytochrome P450 is part of the regulatory mechanism for flavonoid metabolism[38]. Association of FBD with the genomic region harbouring P450 would suggest its role in enhanced polyphenol synthesis causing FBD. The co-occurrence of a senescence-related gene MdNAC90 (MD03G1148500; MD11G11679000) and the PPO regulator germin-like proteins (in the FBD-associated paralog regions on Chrs 3 and 11) lends support to an interplay between senescence and oxidation of phenolics and anthocyanins.

    A significant region on Chr9 encompassed a cluster of UGT proteins, which play a role in the regulation of flavonoids and phenolic compounds, as well as converting phloretin to phloridzin[39]. Cytochrome P450, which resided within several FBD-associated regions, is reportedly involved in flavonoid metabolism, such as chlorogenic acid (CGA) acid and phloridzin, which have also been positively associated with suberin production and cell wall disassembly[40, 41].

    Reactive oxygen species (ROS) play an important role in regulating physiological processes in plants, such as senescence[37]. Legay et al.[42] suggested that MdMYB93 (MD15G1369700), found residing in the FBD-associated regions in this study, plays a critical role in remobilisation of flavonoid/phenolic compounds, which can be utilised for detoxification of ROS in the case of oxidative stress. However, flavonoid biosynthesis has also been linked with suberin production causing cuticle cracks in apples[42, 43], and the development of cuticular cracks could accelerate flesh browning as a result of an enhanced oxidative process[44].

    Several HSP (e.g. HSP70, HSP70-1, HSP60, HSP89.1, and HSP DNAJ) and HSF (e.g. HSF4, HSFB4, HSFA6B) resided in the FBD-associated genomic regions. Ferguson et al.[45] showed that, during summer, apple flesh temperature could reach as high as 43 °C, and that an increase in the expression of HSP in apples was associated with high daily flesh temperatures, suggesting a role of HSP to counter heat stress. HSPs have been reported to interact with AP2/ERFs and to play a role in flavonoid biosynthesis and drought tolerance in apple[46]. Heat stress affects lignin accumulation and its substrate, O-phenols, and has been reported to play role in enzymatic browning[47]. Additionally, HSF that regulate HSP expression have also been reported to be regulated by cold stress to generate heat-induced cold tolerance in banana[48]. HSF1 was shown to transcriptionally regulate the promoters of HSP to enhance chilling tolerance in loquat fruit[49]. Activity of the enzymes (PAL, C4H, 4CL) of the phenylpropanoid pathway was positively correlated with loquat fruit lignification, whilst suppression of their expression by heat shock treatment and low-temperature conditioning significantly reduced fruit lignification[50].

    Wang et al.[7] showed that ascorbate peroxidase (APX) was among the genes that were upregulated in RF apple compared with in white-fleshed apples. Several genes (e.g. MdAPX1, MdAPX3, MdAPX4, and MdDHAR1) involved in ascorbate synthesis resided in the genomic regions associated with FBD. Co-localisation of MdDHAR and ascorbic acid (AsA) synthesis genes in the FBD-linked genomic regions have been reported[51], suggesting that the low AsA content increases fruit susceptibility to FBD[52]. It has been shown that the expression of MdMYB1 and MdDHAR genes was strongly correlated in RF apples, and that AO and APX were upregulated by anthocyanin regulatory genes[31].

    There were several genes associated with CGA biosynthesis residing in the FBD-associated regions on various linkage groups, including the genes MYB19 (MD07G1268000) and MdC3H (MD15G1436500). Higher concentrations of CGA were reported in transgenic apple lines carrying MdMYB10[12], suggesting a role of CGA metabolism in the expression of FBD[53]. Interestingly, some of the FBD-associated regions (e.g. Chrs 9, 11, 13, 14 and 17) reported here in RF apples coincide with those reported earlier for FBD in white-fleshed apples[18, 36, 53], suggesting some common underlying genetic mechanisms.

    As discussed above, ERFs have been reported to be involved in the accumulation of anthocyanin and PA biosynthesis, while ethylene synthesis proteins also contribute to cell membrane breakdown, allowing the PPO enzyme to come into contact with phenolic compounds, potentially leading to FBD symptoms[12, 18]. We observed clusters of anthocyanin biosynthesis proteins (bHLH), ERFs (MdERF2, MdERF4, MdERF12) and PPO genes together in the genomic region associated with WCI on Chr10. The co-occurrence of these gene families perhaps facilitates potential interactions that contribute to the genetic correlation between WCI and FBD. We noted that genes involved in flavonoid regulation and ethylene synthesis occurred together in the FBD-associated regions on several chromosomes (e.g. MD16G1140800 and MdPAE10: MD16G1132100; MdERF1B: MD16G1216900 and MdMYB15: MD16G1218000, MD16G1218900) – suggesting these genes could be in linkage disequilibrium and this would contribute to the expression of WCI and FBD.

    MYB7 (MD16G1029400) resided in the WCI-associated upper region of Chr16, and the expression level of MYB7 was shown to be correlated with that of LAR1 in peach fruit[54]. The MdLAR1 protein (MDP0000376284), which is located about 1.2 Mb upstream of MYB7 (MD16G1029400), was reported to be associated with WCI and FBD in apple[3]. Mellidou et al.[36] reported that 4CL (MD13G1257800 in the FBD-linked region on Chr13), which catalyses the last step of the phenylpropanoid pathway, leading either to lignin or to flavonoids, was upregulated in browning-affected flesh tissues. The gene MdMYB85, involved in the regulation of flavonoid and lignin biosynthesis, resided in the WCI-associated region on Chr6 and FBD-associated paralogous region on Chr4. Metabolic interactions between anthocyanin and lignin biosynthesis have been reported for apple[55] and strawberry[56], while flesh lignification and internal browning during low-temperature storage in a red-fleshed loquat cultivar was shown to be modulated by the interplay between ERF39 and MYB8[57].

    The co-localisation of MdMYB66 and cytochrome P450 proteins, along with the anthocyanin regulatory protein MdMYB86, in paralogous FBD-associated genomic regions on Chr6 and Chr14 suggests that these genes interact as a 'hub' contributing to the WCI-FBD genetic link. The paralogs of some other genes were found to be residing in the regions associated with either WCI or FBD. For example, MdNAC42 and HSP70 co-localised in the FBD-associated region on Chr9, but this same pair of genes also resided in the most prominent region associated with WCI on Chr17. Similarly, paralogs of MdEIN3 resided in the WCI-associated region on Chr7 (MD07G1053500, MD07G1053800) and FBD-associated region on Chr11 (MD11G1022400). Interestingly, the paralogs of MdMYB73 and MdWRKY7 co-localised in the WCI- (24.6–26.8 Mb) and FBD-associated (4.6–6.8 Mb) regions on Chr15. The WCI-associated region at the bottom of Chr11 hosted a cluster of senescence-related genes, along with the anthocyanin biosynthesis gene MdbHLH3 (MD11G1286900), suggesting they might interact in the genetic nexus between FBD and WCI.

    FBD in RF apples can be caused by senescence, injury via extreme temperature exposure (chilling or heat), or enzymatic (cut fruit) reaction. Genes reported to be connected to all three factors were located in various FBD-associated genomic regions in this study. Postharvest strategies that both delay senescence and limit exposure to low temperatures may be needed to manage FBD. We also hypothesise that high ascorbic acid content could help to minimise expression of FBD in Type-1 RF cultivars. The adverse genetic correlation between WCI and FBD appears to arise from dual and/or interactive roles of several transcription factors, which would pose challenges for designing a conventional marker-assisted selection strategy. The use of bivariate genomic BLUP to estimate breeding values to simultaneously improve adversely correlated polygenic traits (e.g. WCI and FBD), could be an alternative approach[58, 59].

    A population of 900 apple seedlings composed of 24 full-sib families was generated in 2011 by selected crossings between six red-leaved pollen parents and six white-fleshed female parents. All six pollen parents inherited their red-leaf phenotype from the same great-grandparent 'Redfield'[1]. Each pollen parent was involved in four crosses, and the female parents were involved in three to six crosses each. Foliage colour of young seedlings is a phenotypic marker for Type 1 RF apple. The main purpose of this trial was to understand the flesh colour variation and FBD in the Type 1 RF seedlings, so only the seedlings with red foliage (i.e. carrying MdMYB10) were kept for this trial. The number of seedlings per family varied from 10 to 95. The seedlings were grafted onto 'M9' rootstock and were planted in duplicate at the Plant & Food Research orchard in Hawke's Bay, New Zealand (39°39′ S, 176°53′ E) in 2015.

    Phenotyping for RF and FBD was conducted over two consecutive fruiting seasons (2017 and 2018). Fruit were harvested once, when judged mature, based on a change in skin background colour from green to yellow, and when the starch pattern index (SPI) was between 1.0 and 2.0 (on a scale of 0 to 7). In each season, six fruit were harvested from each plant and stored for 70 d at 0.5 °C, followed by 7 d at 20 °C before fruit evaluation. Fruit were cut in half across the equator and the proportion of the cortex area (PRA) that was red in colour, and the intensity of the red colour (RI) (= 1 (low) to 9 (high)) was scored. A weighted cortical intensity (WCI) was then calculated (PRA × RI) as an estimation of the amount of red pigment in the fruit. The proportion of the cortex area showing symptoms of FBD was also recorded. WCI and FBD were averaged over all fruit for a particular seedling. Fruit were also assessed for the following eating quality traits on a 1 (lowest) to 9 (highest) scale: firmness, crispness, juiciness, sweetness, sourness and astringency, to understand the genetic correlations of eating-quality traits with WCI and FBD.

    The binary vector pSA277-R6:MYB10 was transferred into Agrobacterium tumefaciens strain LAB4404 by electroporation. Transgenic 'Royal Gala' plants were generated by Agrobacterium-mediated transformation of leaf pieces, using a method previously reported[5]. Wild-type 'Royal Gala' and three independent transgenic lines (A2, A4 and A10) of R6:MdMYB10 were grown under glasshouse conditions in full potting mix with natural light. The resulting fruit were assessed for flesh colour phenotypes at harvest (around 135 d after full bloom). Fruit peel and cortex from three biological replicates were collected and frozen in liquid nitrogen, with each replicate compiled from five pooled mature fruit for each transgenic line or wild-type control.

    Total RNA of 36 samples (3 R6:MYB10 lines and 1 wild type control, 3 time points, 3 biological replicates) was extracted, using Spectrum Plant Total RNA Kit (SIGMA). Removal of genomic DNA contamination and first-strand cDNA synthesis were carried out using the mixture of oligo (dT) and random primers according to the manufacturer's instructions (QuantiTect Reverse Transcription Kit, Qiagen). Real-time qPCR DNA amplification and analysis was carried out using the LightCycler 480 Real-Time PCR System (Roche), with LightCycler 480 software version 1.5. The LightCycler 480 SYBR Green I Master Mix (Roche) was used following the manufacturer's method. The qPCR conditions were 5 min at 95 °C, followed by 45 cycles of 5 s at 95 °C, 5 s at 60°C, and 10 s at 72 °C, followed by 65 °C to 95 °C melting curve detection. The qPCR efficiency of each gene was obtained by analyzing the standard curve of a cDNA serial dilution of that gene. The expression was normalized to Malus × domestica elongation factor 1-alpha MdEF1α (XM_008367439) due to its consistent transcript levels throughout samples, with crossing threshold values changing by less than 2.

    Individual fruit measurements were first averaged for each seedling. As the phenotyping was repeated over two years, we used a mixed linear model (MLM) accounting for this 'permanent environmental effect', as previously described[58]. Pedigree-based additive genetic relationships among seedlings were taken into account for estimation of genetic parameters using ASReml software[60]. Product-moment correlations between best linear unbiased predictions (BLUP) of breeding values of all seedlings for different traits were used as estimates of genetic correlation among traits.

    A selective DNA pooling procedure was adopted to construct DNA pools. A high-pool and a low-pool were constructed separately for the two traits (WCI and FBD). Genomic DNA was extracted from the leaves of selected seedlings, and quantified by fluorimetry using the picogreen reagent (Cat#P11496, Thermo). The low and high pools consisted of 35 seedlings each, and normalised amounts (~300 ng) of DNA from individuals were pooled. The pools were dried down with DNA Stable reagent (Cat#93021001, Biomatrica) in a centrifugal evaporator and shipped for sequencing. Each DNA pool was sequenced using paired-end 125 bp reads on the Illumina HiSeq 2500 platform. The quality of raw sequence reads was checked with FastQC/0.11.2 and MultiQC/1.2. Based on the quality control reports, the reads were aligned to the published apple reference genome GDDH13 v1.1[61] using the program bowtie2/2.3.4.3[62] with trimming from both ends before alignments and aligning in full read length ("-5 6 -3 5 –end-to-end"). The mapping results were marked for duplicate alignments, sorted, compressed and indexed with samtools/1.12[63]. Based on the alignment of binary alignment map (BAM) files of the high and low pools, single nucleotide polymorphism (SNP) identification was performed using samtools/1.12 ('samtools mpileup') and bcftools/1.12 ('bcftools call –mv')[63, 64]. Variant sites with missing genotype in any of the pools, or having the same genotype between the pools, were discarded. To minimise the influence of sequencing quality on association analysis, the identified SNPs were further filtered according to the following criteria: 1) a Phred-scaled quality score > 20; and 2) the read depth in each pool was neither < 35, nor > 500.

    The allele frequencies between each pair of bulk DNAs (low versus high WCI; low versus high FBD) were compared at each SNP locus. Differences in the allele frequencies between the low and high pool were expected to be negligible for unlinked SNP markers, but allele frequency differences would be larger for SNPs closely linked to the underlying quantitative trait loci (QTLs) contributing to the extreme phenotypes. A nonparametric test (G-statistic = 2 × Σni ln(ni/nexp), where ni (i = 1 to 4) represented counts of reference and alternate alleles at a particular SNP generated from sequencing of the low and high pool, and nexp was the expected allele count assuming no allele frequency divergence between the two DNA pools[65].

    We then calculated a modified statistic (G'), which took into account read count variation caused by sampling of segregants as well as variability inherent in short-read sequencing of pooled samples[65]. Using R package QTLseqr[66], firstly a G-statistic was calculated for each SNP marker, and then a weighted average using Nadaraya-Watson kernel was obtained to yield a G' statistic for a sliding genomic window of 2 Mb size. The Nadaraya-Watson method weights neighbouring markers' G-values by their distance from the focal SNP so that closer SNPs receive higher weights. The 95th percentile value of G' was used as a threshold to identify significant hotspots and to identify the putative candidate genes residing within the ±1 Mb region around the G' peak. For comparison purposes, a standard two-sided Z-test[14] was also performed to determine the significance of allele frequency differences at SNP loci between the pools for each trait.

    The GDDH gene models intersecting with the XP-GWA hotspots were pulled out with bedtools/2.30.0 ("bedtools intersect -wo -nonamecheck"). The selected genes were further blasted to TAIR10 ("-evalue 1e-5") and the annotated functions from Arabidopsis genes with the best blast score, the highest % identity, and the longest aligned length, were used. Then the expressions of genes located in the GWA hotspots were extracted from the RNAseq analysis of the R6:MdMYB10 representative transgenic line, and the log 2-fold change between the R6:MdMYB10 and 'Royal Gala' apples were calculated.

    This research was funded in 2017/18 by the Strategic Science Investment Fund of the New Zealand Ministry of Business, Innovation and Employment (MBIE) and from 2019 by the Plant & Food Research Technology Development – Pipfruit programme. We thank our colleague Jason Johnston for providing some pictures of the flesh browning disorder in red-fleshed apples. Richard Volz and Jason Johnston provided constructive comments and suggestions on the manuscript.

  • The authors declare that they have no conflict of interest.

  • Supplemental Table S1 Welch′s t test values p values for stem segments not treated with a lanolin paste.
    Supplemental Fig. S1 Examples of the responses as they appeared on March 19th of stem segments having buds and/or leaves and not treated with a phytohormone. (a)−(c) Apical, mid‐stem and basal regions, respectively, of segments having both buds and leaves present. (d)−(f) Apical, mid‐stem and basal regions, respectively, of segment having only a whorl of terminal buds present.  (g)−(i) Apical, mid‐stem and basal regions of a segment having only leaves present. The inset images in (d), (g) and (i) show the appearance of the cambial zone and adjoining cells at those positions as they appeared in radial section.
  • [1]

    Chaffey N, Cholewa E, Regan S, Sundberg B. 2002. Secondary xylem development in Arabidopsis: a model for wood formation. Physiologia Plantarum 114:594−600

    doi: 10.1034/j.1399-3054.2002.1140413.x

    CrossRef   Google Scholar

    [2]

    Du J, Miura E, Robischon M, Martinez C, Groover A. 2011. The Populus class III HD ZIP transcription factor POPCORONA affects cell differentiation during secondary growth of woody stems. PLoS One 6:e17458

    doi: 10.1371/journal.pone.0017458

    CrossRef   Google Scholar

    [3]

    Savidge RA. 2023. Intrinsic regulation of diameter growth in woody plants. Forests 14:1065

    doi: 10.3390/f14051065

    CrossRef   Google Scholar

    [4]

    Mazur E, Benková E, Friml J. 2016. Vascular cambium regeneration and vessel formation in wounded inflorescence stems of Arabidopsis. Science Reports 6:33754

    doi: 10.1038/srep33754

    CrossRef   Google Scholar

    [5]

    Turley EK, Etchells JP. 2022. Laying it on thick: a study in secondary growth. Journal of Experimental Botany 73:665−79

    doi: 10.1093/jxb/erab455

    CrossRef   Google Scholar

    [6]

    Bueno N, Cuesta C, Centeno ML, Ordás RJ, Alvarez JM. 2021. In vitro plant regeneration in conifers: the role of WOX and KNOX gene families. Genes 12:438

    doi: 10.3390/genes12030438

    CrossRef   Google Scholar

    [7]

    Uddenberg D, Akhter S, Ramachandran P, Sundström JF, Carlsbecker A. 2015. Sequenced genomes and rapidly emerging technologies pave the way for conifer evolutionary developmental biology. Frontiers in Plant Science 6:970

    doi: 10.3389/fpls.2015.00970

    CrossRef   Google Scholar

    [8]

    Zhang J, Eswaran G, Alonso-Serra J, Kucukoglu M, Xiang J, et al. 2019. Transcriptional regulatory framework for vascular cambium development in Arabidopsis roots. Nature Plants 5:1033−42

    doi: 10.1038/s41477-019-0522-9

    CrossRef   Google Scholar

    [9]

    Rashotte AM, Poupart J, Waddell CS, Muday GK. 2003. Transport of the two natural auxins, indole-3-butyric acid and indole-3-acetic acid, in Arabidopsis. Plant Physiology 133:761−72

    doi: 10.1104/pp.103.022582

    CrossRef   Google Scholar

    [10]

    Ludwig-Müller J, Vertocnik A, Town CD. 2005. Analysis of indole-3-butyric acid-induced adventitious root formation on Arabidopsis stem segments. Journal of Experimental Botany 56:2095−105

    doi: 10.1093/jxb/eri208

    CrossRef   Google Scholar

    [11]

    Strader LC, Monroe-Augustus M, Rogers KC, Lin GL, Bartel B. 2008. Arabidopsis iba response5 suppressors separate responses to various hormones. Genetics 180:2019−31

    doi: 10.1534/genetics.108.091512

    CrossRef   Google Scholar

    [12]

    Wang T, Li P, Mu T, Dong G, Zheng C, et al. 2020. Overexpression of UGT74E2, an Arabidopsis IBA glycosyltransferase, enhances seed germination and modulates stress tolerance via ABA signaling in rice. International Journal of Molecular Science 21:7239

    doi: 10.3390/ijms21197239

    CrossRef   Google Scholar

    [13]

    Zolman BK, Bartel B. 2004. An Arabidopsis indole-3-butyric acid-response mutant defective in PEROXIN6, an apparent ATPase implicated in peroxisomal function. Proceedings of the National Academy of Sciences of the United States of America 101:1786−91

    doi: 10.1073/pnas.0304368101

    CrossRef   Google Scholar

    [14]

    Damodaran S, Strader LC. 2019. Indole-3-butyric acid metabolism and transport in Arabidopsis thaliana. Frontiers in Plant Science 10:851

    doi: 10.3389/fpls.2019.00851

    CrossRef   Google Scholar

    [15]

    Yu Y, Meng N, Chen S, Zhang H, Liu Z, et al. 2022. Transcriptomic profiles of poplar (Populus simonii × P. nigra) cuttings during adventitious root formation. Frontiers in Genetics 13:968544

    doi: 10.3389/fgene.2022.968544

    CrossRef   Google Scholar

    [16]

    Van de Poel B, Van Der Straeten D. 2014. 1-aminocyclopropane-1-carboxylic acid (ACC) in plants: more than just the precursor of ethylene! Frontiers in Plant Science 5:640

    doi: 10.3389/fpls.2014.00640

    CrossRef   Google Scholar

    [17]

    Vanderstraeten L, Depaepe T, Bertrand S, Van Der Straeten D. 2019. The ethylene precursor ACC affects early vegetative development independently of ethylene signaling. Frontiers in Plant Science 10:1591

    doi: 10.3389/fpls.2019.01591

    CrossRef   Google Scholar

    [18]

    Vandenbussche F, Vriezen WH, Smalle J, Laarhoven LJJ, Harren FJM, et al. 2003. Ethylene and auxin control the Arabidopsis response to decreased light intensity. Plant Physiology 133:517−27

    doi: 10.1104/pp.103.022665

    CrossRef   Google Scholar

    [19]

    Stepanova AN, Yun J, Likhacheva AV, Alonso JM. 2007. Multilevel interactions between ethylene and auxin in Arabidopsis roots. The Plant Cell 19:2169−85

    doi: 10.1105/tpc.107.052068

    CrossRef   Google Scholar

    [20]

    Poupin MJ, Greve M, Carmona V, Pinedo I. 2016. A complex molecular interplay of auxin and ethylene signaling pathways is involved in Arabidopsis growth promotion by Burkholderia phytofirmans PsJN. Frontiers in Plant Science 7:492

    doi: 10.3389/fpls.2016.00492

    CrossRef   Google Scholar

    [21]

    Zemlyanskaya EV, Omelyanchuk NA, Ubogoeva EV, Mironova VV. 2018. Deciphering auxin-ethylene crosstalk at a systems level. International Journal of Molecular Sciences 19:4060

    doi: 10.3390/ijms19124060

    CrossRef   Google Scholar

    [22]

    Sundberg B, Little CHA, Riding RT, Sandberg G. 1987. Levels of endogenous indole-3-acetic acid in the vascular cambium region of Abies balsamea trees during the activity - rest - quiescence transition. Physiologia Plantarum 71:163−70

    doi: 10.1111/j.1399-3054.1987.tb02862.x

    CrossRef   Google Scholar

    [23]

    Sundberg B, Little CHA. 1987. Effect of defoliation on tracheid production and the level of indole-3-acetic acid in Abies balsamea shoots. Physiologia Plantarum 71:430−35

    doi: 10.1111/j.1399-3054.1987.tb02879.x

    CrossRef   Google Scholar

    [24]

    Little CHA, Sundberg B, Ericsson A. 1990. Induction of acropetal 14C-photosynthate transport and radial growth by indole-3-acetic acid in Pinus sylvestris shoots. Tree Physiology 6:177−89

    doi: 10.1093/treephys/6.2.177

    CrossRef   Google Scholar

    [25]

    Savidge RA. 1988. Auxin and ethylene regulation of diameter growth in trees. Tree Physiology 4:401−14

    doi: 10.1093/treephys/4.4.401

    CrossRef   Google Scholar

    [26]

    Savidge RA, Mutumba GMC, Heald JK, Wareing PF. 1983. Gas chromatography-mass spectroscopy identification of 1-aminocyclopropane-1-carboxylic acid in compressionwood vascular cambium of Pinus contorta Dougl. Plant Physiology 71:434−36

    doi: 10.1104/pp.71.2.434

    CrossRef   Google Scholar

    [27]

    Eklund L, Little CHA. 1995. Interaction between indole-3-acetic acid and ethylene in the control of tracheid production in detached shoots of Abies balsamea. Tree Physiology 15:27−34

    doi: 10.1093/treephys/15.1.27

    CrossRef   Google Scholar

    [28]

    Eklund L, Little CHA. 1996. Laterally applied ethrel causes local increases in radial growth and indole-3-acetic acid concentration in Abies balsamea shoots. Tree Physiology 16:509−13

    doi: 10.1093/treephys/16.5.509

    CrossRef   Google Scholar

    [29]

    Eklund L, Little CHA. 1998. Ethylene evolution, radial growth and carbohydrate concentrations in Abies balsamea shoots ringed with ethrel. Tree Physiology 18:383−91

    doi: 10.1093/treephys/18.6.383

    CrossRef   Google Scholar

    [30]

    Little CHA, Eklund L. 1999. Ethylene in relation to compression wood formation in Abies balsamea shoots. Trees 13:173−77

    doi: 10.1007/PL00009749

    CrossRef   Google Scholar

    [31]

    Little CHA, Eklund L. 2000. Transport of [1-14C]-indole-3-acetic acid in Abies balsamea shoots ringed with ethrel. Trees 15:58−62

    doi: 10.1007/s004680000075

    CrossRef   Google Scholar

    [32]

    Savidge RA. 1982. Regulation of seasonal cambial activity and tracheid differentiation in Pinus contorta Dougl. Thesis. University of Wales, UK.

    [33]

    Savidge RA. 1983. The role of plant hormones in higher plant cellular differentiation. ll. Experiments with the vascular cambium, and sclereid and tracheid differentiation in the pine, Pinus contorta. The Histochemical Journal 15:447−66

    doi: 10.1007/BF01002699

    CrossRef   Google Scholar

    [34]

    Worrall JJ, Little CHA. 1986. An effect of gravity on bud-burst in balsam fir. Tree Physiology 1:47−52

    doi: 10.1093/treephys/1.1.47

    CrossRef   Google Scholar

    [35]

    Little CHA, Lavigne MB. 2002. Gravimorphism in current-year shoots of Abies balsamea: involvement of compensatory growth, indole-3-acetic acid transport and compression wood formation. Tree Physiology 22:311−20

    doi: 10.1093/treephys/22.5.311

    CrossRef   Google Scholar

    [36]

    O'Brien TP, Feder N, McCully ME. 1964. Polychromatic staining of plant cell walls by toluidine blue O. Protoplasma 59:367−73

    doi: 10.1007/BF01248568

    CrossRef   Google Scholar

    [37]

    Savidge RA, Wareing PF. 1981. Plant-growth regulators and the differentiation of vascular elements. In Xylem Cell Development, ed. Barnett JR. Tunbridge Wells, UK: Castle House Publications. pp. 192−235

    [38]

    Savidge RA. 2001. Intrinsic regulation of cambial growth. Journal of Plant Growth Regulation 20:52−77

    doi: 10.1007/s003440010002

    CrossRef   Google Scholar

    [39]

    Rasband WS. 1997−2015. ImageJ. US National Institutes of Health, Bethesda, Maryland, USA. http://imagej.nih.gov/ij

    [40]

    Savidge RA. 2014. Cell biology of bordered-pit formation in balsam-fir trees. Botany 92:495−511

    doi: 10.1139/cjb-2014-0030

    CrossRef   Google Scholar

    [41]

    Welch BL. 1951. On the comparison of several mean values: an alternative approach. Biometrika 38:330−36

    doi: 10.1093/biomet/38.3-4.330

    CrossRef   Google Scholar

    [42]

    Midway S, Robertson M, Flinn S, Kaller M. 2020. Comparing multiple comparisons: practical guidance for choosing the best multiple comparisons test. PeerJ 8:e10387

    doi: 10.7717/peerj.10387

    CrossRef   Google Scholar

    [43]

    Kijidani Y, Wu Z, Savidge RA. 2001. New insight into phytohormone regulation of wood formation in conifers. Phytomorphology 51:185−200

    Google Scholar

    [44]

    Savidge RA, Yuan X, Foerster H. 2023. Gišogenetic variation in white-spruce (Picea glauca (Moench) Voss) trees of Yukon Beringia, Canada. Forests 14:787

    doi: 10.3390/f14040787

    CrossRef   Google Scholar

    [45]

    Savidge RA, Wareing PF. 1982. Apparent auxin production and transport during winter in the non-growing pine tree. Canadian Journal of Botany 60:681−91

    doi: 10.1139/b82-090

    CrossRef   Google Scholar

    [46]

    Savidge RA, Wareing PF. 1981. A tracheid-differentiation factor from pine needles. Planta 153:395−404

    doi: 10.1007/BF00394977

    CrossRef   Google Scholar

    [47]

    Savidge RA. 1994. The tracheid-differentiation factor of conifer needles. International Journal of Plant Science 155:272−90

    doi: 10.1086/297167

    CrossRef   Google Scholar

    [48]

    Savidge RA, Barnett JR. 1993. Protoplasmic changes in cambial cells induced by a tracheid-differentiation factor from pine needles. Journal of Experimental Botany 44:395−405

    doi: 10.1093/jxb/44.2.395

    CrossRef   Google Scholar

    [49]

    Savidge RA. 2007. ACC induction of xylogenesis in horseradish roots. Proceedings 33rd PGRSA Annual Meeting, 2006, Quebec City, Canada. pp. 60−65. www.researchgate.net/publication/311654362_ACC_induction_of_xylogenesis_in_horseradish_roots

    [50]

    Hook DD, Brown CL, Wetmore RH. 1972. Aeration in trees. Botanical Gazette 133:443−54

    doi: 10.1086/336669

    CrossRef   Google Scholar

    [51]

    Laming PB. 1974. On intercellular spaces in the xylem ray parenchyma of Picea abies. Acta Botanica Neerlandica 23:217−23

    doi: 10.1111/j.1438-8677.1974.tb00939.x

    CrossRef   Google Scholar

    [52]

    Kitin P, Fujii T, Abe H, Takata K. 2009. Anatomical features that facilitate radial flow across growth rings and from xylem to cambium in Cryptomeria japonica. Annals of Botany 103:1145−57

    doi: 10.1093/aob/mcp050

    CrossRef   Google Scholar

    [53]

    Shimakura M. 1936. On the expansion of bast cells in conifers. The Botanical Magazine 50:318−23

    doi: 10.15281/jplantres1887.50.318

    CrossRef   Google Scholar

    [54]

    Morgensen HL. 1968. Studies on the bark of the cork bark fir: Abies lasioearpa var. arizonica (merriam) lemmon. I. periderm ontogeny. Journal of the Arizona Academy of Science 5:36−40

    doi: 10.2307/40022822

    CrossRef   Google Scholar

    [55]

    Golonowski WO. 1971. The anatomical structure of the common fir (Abies alba Mill.) bark. I. development of bark tissues. Acta Societatis Botanicorum Poloniae 40:149−81

    doi: 10.5586/asbp.1971.010

    CrossRef   Google Scholar

    [56]

    Grozdits GA, Godkin SE, Keith CT. 1982. The periderms of three North American conifers. Wood Science and Technology 16:305−16

    doi: 10.1007/BF00353159

    CrossRef   Google Scholar

    [57]

    Lev-Yadun S, Liphschitz N. 1989. Sites of first phellogen initiation in conifers. IAWA Journal 10:43−52

    Google Scholar

    [58]

    Xiao W, Molina D, Wunderling A, Ripper D, Vermeer JEM, et al. 2020. Pluripotent pericycle cells trigger different growth outputs by integrating developmental cues into distinct regulatory modules. Current Biology 30:4384−4398.e5

    doi: 10.1016/j.cub.2020.08.053

    CrossRef   Google Scholar

    [59]

    Johns S, Hagihara T, Toyota M, Gilroy S. 2021. The fast and the furious: rapid long-range signaling in plants. Plant Physiology 185:694−706

    doi: 10.1093/plphys/kiaa098

    CrossRef   Google Scholar

    [60]

    Moungsrimuangdee B, Iwanaga F, Yamanaka N, Kodani J, Yamamoto F. 2022. Effects of applied ethrel, jasmonic acid, and salicylic acid on the formation of traumatic resin ducts in the bark of Thujopsis dolabrata cuttings and xylem of Metasequoia glyptostroboides seedlings. Trees 36:793−801

    doi: 10.1007/s00468-021-02250-4

    CrossRef   Google Scholar

    [61]

    Fahn A, Zamski E. 1970. The influence of pressure, wind, wounding and growth substances on the rate of resin duct formation in Pinus halepensis wood. Israel Journal of Botany 19:429−46

    Google Scholar

    [62]

    Evans DE. 2004. Aerenchyma formation. New Phytologist 161:35−49

    doi: 10.1046/j.1469-8137.2003.00907.x

    CrossRef   Google Scholar

    [63]

    Namboodiri KK, Beck CB. 1968. A comparative study of the primary vascular system of conifers. 1. Genera with helical phyllotaxis. American Journal of Botany 55:447−57

    doi: 10.1002/j.1537-2197.1968.tb07398.x

    CrossRef   Google Scholar

    [64]

    Fink S. 1983. The occurrence of adventitious and preventitious buds within the bark of some temperate and tropical trees. American Journal of Botany 70:532−42

    doi: 10.1002/j.1537-2197.1983.tb07880.x

    CrossRef   Google Scholar

  • Cite this article

    Savidge RA. 2024. Responses of isolated balsam-fir stem segments to exogenous ACC, IAA, and IBA. Forestry Research 4: e033 doi: 10.48130/forres-0024-0030
    Savidge RA. 2024. Responses of isolated balsam-fir stem segments to exogenous ACC, IAA, and IBA. Forestry Research 4: e033 doi: 10.48130/forres-0024-0030

Figures(9)  /  Tables(1)

Article Metrics

Article views(1526) PDF downloads(174)

Other Articles By Authors

ARTICLE   Open Access    

Responses of isolated balsam-fir stem segments to exogenous ACC, IAA, and IBA

Forestry Research  4 Article number: e033  (2024)  |  Cite this article

Abstract: In this investigation, the effects of exogenous indole-3-acetic acid (IAA), indole-3-butyric acid (IBA), and 1-aminocyclopropane-1-carboxylic acid (ACC) on anatomical development within cultured segments of Abies balsamea (L.) Mill. were compared, using debudded and defoliated leaders produced in the preceding year as bioassay material. In stem apical regions, IAA promoted radial enlargement of pre-existing cortical resin ducts and attending parenchyma enlargement, whereas IBA promoted cell division and expansion of parenchyma on the outer edge of phloem without altering cortical duct shape. Cortical woody ducts, each partially surrounded by cambium, were observed as a novel but infrequent feature. A single cortical woody duct was spatially associated with each mature leaf as its vascular trace, and they were not encountered elsewhere in the cortex, nor were they induced to form in response to any hormone application. An unknown leaf factor induces the development of cortical woody ducts. Both IAA and IBA promoted cell division in the vascular cambium. The common cellular response at the interface between the latewood boundary and cambial zone was the radial expansion of primary-walled fusiform cambial cell derivatives with little if any ensuing tracheary element (TE) differentiation. Enhanced TE production at basal stem positions occurred when ACC was provided with IAA and/or IBA, and an IAA + IBA + ACC combination produced a basal stem response similar to that in untreated segments having intact leaves. The data support the conclusion that IAA, IBA, and ACC have distinct but complementary roles in the overall regulation of the types of cellular differentiation that contribute to cortex histogenesis and diameter growth of balsam-fir leaders.

    • Arabidopsis thaliana L. ('arabidopsis') and Populus ('poplar' and 'aspen') spp. are evolutionarily distantly removed from conifers but, as annual seed plants, have nevertheless become models for explaining plant growth and development[1,2]. Arabidopsis is likely inadequate for revealing all mechanisms of secondary growth in trees, because it does not produce annual layers of vascular and peridermal tissues, such as those found in woody species of the temperate zones[3]. However, arabidopsis plants can produce some secondary vascular tissue within their flowering stems and roots[1,4,5] and Populus spp. share with arabidopsis many regulatory genes[2,5]. There is growing evidence that organogenesis, histogenesis and cellular differentiation in conifers may involve regulatory genes in common with those eudicots[6], although important differences remain to be resolved[7].

      Arabidopsis investigations revealed that indole-3-acetic acid (IAA) and indole-3-butyric acid (IBA) have distinct regulatory roles[811] and are transmitted in different pathways through plant tissues[1214]. These findings have yet to be confirmed in Populus spp., but IBA is effective at inducing rooting of poplar cuttings[15].

      The ethylene precursor 1-aminocyclopropane-1-carboxylic acid (ACC) has a variety of roles in regulating arabidopsis and poplar development[1618], and many interactions between auxin and ethylene have been described[1921].

      In balsam fir (Abies balsamea (L.) Mill.), research has indicated that IAA[2224], ACC[25,26], and ethylene[2731] have roles as endogenous regulators. IBA as an endogenous hormone has been largely overlooked in conifers, although a full-scan mass spectrum established its presence in the cambial tissue of Pinus contorta Dougl. long ago[32], and exogenous IBA applied to pine stem segments promoted xylogenesis although weakly relative to IAA[33]. Given the above-described advances with eudicots, this investigation aimed to discover what, if any, effects on secondary growth in balsam-fir leaders might occur in response to IAA, IBA, ACC, and combinations thereof.

      Before the balsam-fir segments used in this investigation were treated with hormones, they were debudded and defoliated, so that their only available reserves were those stored within stem tissues. To avoid the well known memory effects induced by gravity on tilted stems and branches[34,35], the stem segments were prepared exclusively from zenith-grown main-stem leaders. After culturing treated stem segments for four weeks under identical conditions, microscopy data were produced and are presented below.

    • On Feb. 16th, dormant main-stem leaders between 24 and 36 cm in length and produced during the previous year's growing season were collected (Fig. 1). The leaders had no obvious defects and displayed perfectly zenith growth (i.e., negative geotropism). The source trees were a wild-type population of saplings, between 1−2 m in height, growing in wild forest at Fredericton, New Brunswick, Canada (45°55'38.56"N, 66°37'54.18"W).

      Figure 1. 

      Flowchart showing the progression of the investigation and examples of stem anatomy at the beginning and end. The stem segment at lower left shows stem diameters below the apical whorl of buds, and the segment at lower right shows diameters at the base of the leader.

      Not all stem segments survive an extended incubation period after having been debudded and defoliated[32]. The explanation for this variability remains uncertain, but with the aim of having at least three replicates for analysis, six stem segments were used for each treatment.

      Using a razor blade, leaders were subdivided under cold tap water into 9-cm-long segments and, except where otherwise noted, all leaves and buds were removed by careful razor-blade trimming. Leaves were cut away through the petiole region just above their circular attachment pads. The terminal bud cluster was fully removed by a transverse stem cut at its base, and axillary buds were sliced off tangentially flush with the epidermis.

      To know how intact buds and leaves influenced diameter growth, three additional treatments were prepared: segments having both buds and leaves; segments having leaves but with buds removed; and segments with buds intact and leaves removed.

      The segments prepared as described above were inserted basal ends down, to a depth of 3.5 cm, into plastic containers (7 cm × 7 cm) containing prewashed Oasis 5,200 horticube foam bricks (Smithers Oasis, Kent, OH, USA 44240). Water to saturate the foam was added to the surface level of the foam and replenished regularly as needed. Six segments were inserted into each container. From the leader's terminal bud whorl to the branch internode below, stems increase slightly in diameter. Pre-examination of eight dormant leaders at several distances below the apical buds revealed no obvious difference in anatomy. Figure 1 illustrates this size difference by showing the stem near the apical whorl of buds on Feb. 16th and a stem segment taken from a more basal leader position on March 19th. Stem segments used in the various treatments were selected randomly from the total pool of stem segments.

    • IAA (Sigma-Aldrich Chemicals I3750) was dissolved in absolute ethanol, and IBA-K+ (Sigma-Aldrich Chemicals I7152) and ACC (Sigma-Aldrich Chemicals A3903) were dissolved in distilled H2O to known concentrations. Appropriate volumetric aliquots of those concentrated solutions, together with additional ethanol and/or H2O were transferred into 50 g of warm (70 °C) liquid lanolin (Fisher L7-500), sealed and shaken thoroughly, to produce lanolin mixtures having 1.0 mL of ethanol and 1.0 mL of H2O. The pastes as so prepared contained 3.0 μmols of each test compound per gram of lanolin. The 'control' paste was 50 g of lanolin that had been supplemented with 1.0 mL of ethanol and 1.0 mL of H2O.

      After thorough mixing, 1.0 ± 0.1 g of warm liquid lanolin was transferred into a pre-weighted gelatin half-capsule, cooled to solidify and the half-capsule upended over the apical end of a stem segment (Fig. 1), repeating the process for each stem segment.

      The experiment was incubated at a daytime temperature of 22 °C, gradually declining after 12 h to 15 °C and climbing back again daily to 22 °C. All stem segments treated with lanolin paste were shielded from light using a loosely fitted cover of aluminum foil.

      Stem segments not provided lanolin were incubated in light. The photoperiod was 16 h using natural daylight illumination supplemented with 60 W incandescent lighting after sunset.

      After a treatment period extending from Feb. 16th to March 19th, the stem segments were transferred into methanol until sectioned for microscopy (Fig. 1).

    • After rinsing off surface methanol, each stem segment was sectioned at three positions (Fig. 1): 1 cm below the apical end, 1 cm above the basal end, and halfway between the two ends. Both transverse and radial sections were cut by hand-held razor blade from the two ends, and transverse sections only were cut from the mid-stem region. Sections were stained for 30 s in 0.05% (w/v) aqueous (unbuffered) toluidine blue O[36], rinsed in distilled water and examined by bright-field microscopy. Photomicroscopy images were obtained at 32× and 500× magnifications.

      Overview (32× magnification) analysis of cross sections was used to assess if there was evidence for circumferential variation in a stem's cambial response, but such variation was mostly lacking. Cell types outward from the latewood boundary were counted in six adjoining radial files, distinguishing tracheids (TEs), radially expanded thin-walled cells centripetal to the cambial zone (RE), fusiform cambial cells of the cambial zone (CZ), and enlarging and enlarged phloem cells (Ph) centrifugal to the CZ[37,38]. The average number per radial file of each cell type was recorded.

      New TEs produced in response to a treatment could be confidently distinguished by reference to the latewood boundary, also by their thinner secondary walls, larger radial diameters, and darker secondary-wall staining than those of preexisting latewood. Cambial derivatives having actively lignifying (blue-green stained) secondary walls were included in the TE count. When uncertain by examination of cross sections if a radially expanded thin-walled cambial derivative was undergoing lignification, radial sections were examined to search for evidence for bordered-pit formation and, when developing bordered pits were found, the cell was considered a TE, whereas if lacking the cell was counted as a radially expanded primary-walled (RE) derivative.

      Radial file counts of phloem cells (Ph) included all cells adjoining and centrifugal to the cambial zone to the cortex boundary that displayed no signs of radial compression. The counts do not necessarily comprise newly produced phloem cells; some probably pre-existed in the dormant state.

      Linear measurements of photomicroscopy images were done using ImageJ software[39] and a 1.00 mm (10 μm divisions) stage micrometer as a calibration standard.

    • SEM preparation and analyses were done using critical point drying (acetone - CO2) and carbon coating using SEM aquisition conditions described previously[40].

    • Calculations of means, standard deviations, and unpaired Welch t tests, assuming normal distributions, were done using Microsoft Excel and an alpha of 0.05 to generate p values on the basis of unequal variance, as was indicated by the data[41,42].

    • In addition to the findings reported herein, Supplemental Table S1 and Supplemental Fig. S1 in support of the methods and results are available online.

    • After 28 d, at the termination of the experimental trial period, visual examination and microscopy data indicated that only four stem segments of some treatments had apical, mid-stem and basal tissues that had all survived. For example, two that had been provided IAA showed no apical swelling, and they displayed cortical shrinkage responses in their lower regions (Fig. 1). Viewed under the microscope, their internal tissues were brown and evidently dead; therefore, they were rejected.

      Segments treated with either IAA or IBA either individually or in combination displayed relatively strong diameter growth to produce swollen stem apices. Stem swelling was most pronounced at approximately 1 cm below the apical end of the stem segment, and though tapering basally was evident for an axial distance of 2−3 cm below the application site (Fig. 1). Lanolin and ACC resulted in little if any apical swelling, but slight basal swelling occurred in response to ACC. Basal swelling was more evident in response to treatments that combined IAA or IBA with ACC, and basal swelling was most obvious on segments having both buds and leaves intact (see Supplemental Fig. S1).

    • At the start of the investigation, the vascular cambium was dormant. First periderm formation was evident at scattered locations around circumferences. Common throughout the cortex were vertical files of small parenchyma cells, some circular and others rectangular in outline, with diameters between 5 and 20 μm. Transverse diameters of cortical resin ducts were found to be between 50 μm and > 200 μm in diameter. Each duct was surrounded by two to three tiers of small-diameter sheath cells. Ducts were conspicuously distributed throughout the cortex, spaced 200 μm or more apart within smaller diameter parenchyma. Cortical ducts having the largest diameters appeared well spaced apart around a circumference occupying a position approximately 1/3 the transverse radial distance from the outermost phloem to the epidermis (Fig. 1, see Feb. 16th). Cortical ducts of smaller diameters were abundantly present both on the inside and outside of that circumferential zone of larger ducts. In the interfacial region between outermost phloem cells and cortical tissue, incompletely developed presumably nascent resin ducts appeared as enlarged thin-walled parenchyma cells having diameters between 50 and 100 μm and lacking one or more of the tiers of circumferentially bounding parenchyma that are associated with mature resin ducts.

      In addition to resin ducts and adjoining parenchyma cells, there were within the dormant cortex aerenchyma surrounded by intercellular spaces. Rarely seen were elongated fibres, some of which were thin-walled though many displayed secondary walls that stained blue-green with toluidine blue. Some but not all parenchyma cells also yielded similar histochemical evidence for lignin or suberin. Widely spaced polyphenolic parenchyma were present but rare in the cortex; they were more commonly associated with the phloem.

    • Low magnification cross sections of the apical regions of the stem segments show the three tissue systems (cortical, vascular, pith) present (Figs 2, 3). Higher magnification images are provided in Figs 47.

      Figure 2. 

      Overview cross sections of stem segments at an apical position approximately 1 cm below the hormone application site. The text at upper right in each photomicrograph indicates the treatment. (a) Black arrow points to evidence for the first periderm enveloping the cortex; the white arrow points to an enlarged cortical resin duct; numerous small ones are also evident. (b) Labels indicate cortex, phloem (Ph) and secondary xylem (xylem) locations. (c) Thinner section than those shown in Fig. 2a & b but otherwise similar. (d) Parenchyma have proliferated and preexisting resin ducts are enlarged or stretched radially. (e) Note the diffuse porosity throughout the cortex and the band of large diameter resin ducts immediately centrifugal to the phloem. (f) Radial enlargement of preexisting resin ducts and proliferation of cortical parenchyma, but lacking is an external band of new resin ducts just beyond the phloem. (g) Preexisting resin ducts are radially enlarged (arrow). (h) Circumferential band of resin ducts has developed external to the phloem. (i) A similar response to that shown in Fig. 2e.

      Figure 3. 

      Cross sectional anatomy of the cortex surrounding secondary xylem (staining blue green) at approximately 1 cm below the apical ends of stem segments. (a) Lanolin control, showing cortical resin ducts (d), mature phloem (ph) and secondary xylem (x). The arrow indicates dormant vascular cambium. (b) ACC treatment, showing evidence for browning (arrow) of secondary cell walls in the cortex and phloem. (c) IAA treatment, radially elongated ducts (white arrow) and greatly enlarged parenchyma (black arrow). (d) IBA treatment, numerous enlarged parenchyma (white arrows); ducts retained their circular appearance.

      Figure 4. 

      Higher magnification images of cortical tissue. (a) Early stage of cortical duct formation near the outer phloem interface; the wall of a plasmolysing cell is arrowed. (b) Primary wall (arrowed) fully collapsed. (c) Tonoplast (arrowed) of an expanded vacuole in a greatly enlarged cell. (d) Walls of several collapsed cells in the vicinity of a partially collapsed are arrowed; note wall thickness and the staining reaction possibly indicative of lignin or suberin. (e) An early duct (d) surrounded by two tiers of sheath cells with collapsed and collapsing cells (arrowed) nearby. (f) Active production of cells (one arrowed) to produce the inner sheath tier surrounding a duct (d). (g) A duct (d) surrounded by three tiers of sheath cells having contents probably indicative of resin production.

      Figure 5. 

      Woody duct formation as observed in cross sections (a)−(e) and (i), and radial longitudinal sections (f)−(h) of balsam fir cortex. (a) Early stage of a developing woody duct (d) encased by variable number of parenchyma tiers (arrow) and most of the duct opening filled with intrusive parenchyma. (b) A slightly more advanced stage with the former duct opening completely filled with parenchyma, some differentiating woody elements, and nascent cambium (c) developing on its outer periphery. The arrow indicates a ray-like string of enlarged parenchyma cells bisecting the woody element population. (c) A more advanced stage of cambium (c) formation in a developing woody duct. (d) The cambium (c) in this woody duct is fully developed but the woody elements appear to be at different stages of secondary wall formation (arrow). (e) A fully mature woody duct bisected by a radial file of parenchyma cells (p). (f) An intrusive tip (arrow) of a parenchyma cell elongating within a cortical duct. (g) Low magnification, showing a longitudinal strand (arrow) containing woody elements and non-woody parenchyma and running axially through the cortex. (h) A longitudinal section through a cortical woody duct showing its cambium (c) with elongated nuclei (black arrow) and internalized parenchyma (p). The duct's woody elements all appear to have annular ribs of the primary xylem type (white arrow). (i) This low magnification cross section through the circular attachment pad of a mature leaf shows the spatial association of the leaf base to a woody duct (arrow).

      Figure 6. 

      Developmental responses, as seen in cross sections, of stem segments the apical ends of which were treated with lanolin (only), IAA in lanolin, IBA in lanolin, and a combination of IAA and IBA in lanolin. In each of the four columns of photomicrographs, the apical end is shown at the top, followed by the mid-stem region, and the basal end at the bottom. Lanolin column: (a) Non-dividing CZ, several phloem cells per radial file; (b) non-dividing CZ; (c) same as Fig. 6a but with 2−3 RE cells and, as shown in the inset, in scattered locations around the circumference a single TE per radial file. IAA column: (d) new TEs; (e) a single RE and a single thin-walled TE per radial file; (f) RE cells without TEs. IBA column: (g) new TEs; (h) RE only; (i) no RE or TE cells. IAA + IBA column: (j) new TEs; (k) RE cells only; (l) RE cells only.

      Figure 7. 

      Examples of xylogenic responses, as viewed in cross sections, of stem segments the apical ends of which were treated with a combination of IAA + ACC in lanolin, IBA + ACC in lanolin, and IAA + IBA + ACC in lanolin. In each of the three columns of photomicrographs, the apical end is shown at the top, followed by the mid-stem region, and the basal end at the bottom. IAA + ACC column: (a) Several new TEs per radial file; (b) 1−2 TEs; (c) 1−2 RE cells only per radial file. IBA + ACC column: (d) several new TEs; (e) 1−2 RE cells and no TEs; (f) a single RE cell, 1−2 new TEs and a very narrow CZ. IAA + IBA + ACC column: (g) a single RE cell and no new TEs; (h) a narrow CZ and a single TE; (i) several new TE cells per radial file.

      Compared to the strong cortical growth that resulted in swelling at the apical ends of stem segments in response to all treatments involving IAA or IBA (Fig. 1), the vascular cambium's contribution to radial growth in those same stem segments was weak (Table 1), not apparent in Fig. 2 and barely visible in Fig. 3. The pith showed no evidence of any response to any treatment.

      Table 1.  Microscopy data and P values (95% confidence) summarized for the eight treatments.

      Treatment Mean number of cells per radial file
      Apical end Midstem Basal end
      lanolin
      Replicate Ph CZ RE SL Total Ph CZ RE SL Total Ph CZ RE SL Total
      1 3 4 1 0 8 4 4 1 0 9 5 5 1 0 11
      2 4 5 1 1 11 5 4 1 0 10 4 6 1 2 13
      3 4 6 0 0 10 3 5 1 0 9 3 4 2 0 9
      4 2 4 0 0 6 4 6 0 0 10 3 4 1 0 8
      Mean 3.3 4.8 0.5 0.3 8.8 4.0 4.8 0.8 0.0 9.5 3.8 4.8 1.3 0.5 10.3
      Std. Dev. 1.0 1.0 0.6 0.5 2.2 0.8 1.0 0.5 0.0 0.6 1.0 1.0 0.5 1.0 2.2
      ACC
      Replicate Ph CZ RE SL Total Ph CZ RE SL Total Ph CZ RE SL Total
      1 4 4 2 0 10 2 6 0 0 8 4 6 1 2 13
      2 3 4 2 0 9 3 5 0 0 8 3 5 2 1 11
      3 4 4 1 1 10 4 4 1 0 9 3 4 2 0 9
      4 3 5 0 0 8 5 5 1 0 11 4 4 2 0 10
      Mean 3.5 4.3 1.3 0.3 9.3 3.5 5.0 0.5 0.0 9.0 3.5 4.8 1.8 0.8 10.8
      Std. Dev. 0.6 0.5 1.0 0.5 1.0 1.3 0.8 0.6 0.0 1.4 0.6 1.0 0.5 1.0 1.7
      P: ACC vs lanolin 0.337 0.201 0.119 0.500 0.350 0.271 0.353 0.269 1.000 0.274 0.337 0.500 0.104 0.365 0.367
      IAA
      Replicate Ph CZ RE SL Total Ph CZ RE SL Total Ph CZ RE SL Total
      1 5 6 5 3 19 5 4 1 2 12 3 4 2 1 10
      2 4 4 1 3 12 2 2 2 0 6 3 0 4 0 7
      3 2 6 3 2 13 3 6 1 1 11 4 6 3 0 13
      4 2 4 2 3 11 3 4 1 2 10 2 2 4 0 8
      Mean 3.3 5.0 2.8 2.8 13.8 3.3 4.0 1.3 1.3 9.8 3.0 3.0 3.3 0.3 9.5
      Std. Dev. 1.5 1.2 1.7 0.5 3.6 1.3 1.6 0.5 1.0 2.6 0.8 2.6 1.0 0.5 2.6
      P: IAA vs lanolin 0.500 0.375 0.036 0.000 0.032 0.181 0.233 0.104 0.040 0.432 0.140 0.138 0.008 0.338 0.340
      P: IAA vs ACC 0.386 0.149 0.095 0.000 0.042 0.395 0.165 0.049 0.040 0.319 0.180 0.138 0.022 0.201 0.231
      IBA
      Replicate Ph CZ RE SL Total Ph CZ RE SL Total Ph CZ RE SL Total
      1 4 8 4 2 18 6 4 2 0 12 3 6 0 0 9
      2 5 6 2 1 14 5 2 2 1 10 3 2 2 0 7
      3 2 4 1 2 9 2 2 3 0 7 2 2 2 0 6
      4 2 4 1 1 8 2 6 2 0 10 2 2 2 0 6
      Mean 3.3 5.5 2.0 1.5 12.3 3.8 3.5 2.3 0.3 9.8 2.5 3.0 1.5 0.0 7.0
      Std. Dev. 1.5 1.9 1.4 0.6 4.6 2.1 1.9 0.5 0.5 2.1 0.6 2.0 1.0 0.0 1.4
      P: IBA vs lanolin 0.500 0.259 0.061 0.009 0.120 0.416 0.151 0.003 0.196 0.414 0.038 0.092 0.338 0.196 0.028
      P: IBA vs ACC 0.386 0.143 0.209 0.009 0.144 0.423 0.111 0.002 0.196 0.287 0.025 0.092 0.338 0.108 0.008
      P:IBA vs IAA 0.500 0.337 0.262 0.009 0.314 0.348 0.353 0.015 0.065 0.500 0.180 0.500 0.022 0.196 0.081
      IAA + ACC
      Replicate Ph CZ RE SL Total Ph CZ RE SL Total Ph CZ RE SL Total
      1 3 7 3 1 14 2 5 3 0 10 3 7 3 2 15
      2 4 8 2 1 15 3 6 2 2 13 3 8 3 1 15
      3 3 3 4 2 12 2 4 2 0 8 2 6 1 3 12
      4 2 6 2 3 13 1 4 2 2 9 2 5 0 0 7
      Mean 3.0 6.0 2.8 1.8 13.5 2.0 4.8 2.3 1.0 10.0 2.5 6.5 1.8 1.5 12.3
      Std. Dev. 0.8 2.2 1.0 1.0 1.3 0.8 1.0 0.5 1.2 2.2 0.6 1.3 1.5 1.3 3.8
      P: IAA + ACC vs lanolin 0.353 0.174 0.005 0.022 0.007 0.007 0.500 0.003 0.091 0.341 0.038 0.038 0.282 0.135 0.202
      P: IAA + ACC vs ACC 0.180 0.102 0.034 0.022 0.001 0.053 0.353 0.002 0.091 0.236 0.025 0.038 0.500 0.195 0.254
      P: IAA + ACC vs IAA 0.391 0.227 0.500 0.065 0.451 0.077 0.233 0.015 0.375 0.444 0.180 0.033 0.076 0.074 0.141
      P: IAA + ACC vs IBA 0.391 0.371 0.209 0.337 0.318 0.096 0.151 0.500 0.149 0.436 0.500 0.016 0.396 0.051 0.031
      IBA + ACC
      Replicate Ph CZ RE SL Total Ph CZ RE SL Total Ph CZ RE SL Total
      1 5 6 1 5 17 3 5 1 0 9 5 4 1 2 12
      2 6 4 2 2 14 3 4 1 2 10 5 4 2 1 12
      3 3 3 1 2 9 3 4 1 2 10 2 2 1 2 7
      4 6 8 3 5 22 4 6 1 0 11 4 2 3 0 9
      Mean 5.0 5.3 1.8 3.5 15.5 3.3 4.8 1.0 1.0 10.0 4.0 3.0 1.8 1.3 10.0
      Std. Dev. 1.4 2.2 1.0 1.7 5.4 0.5 1.0 0.0 1.2 0.8 1.4 1.2 1.0 1.0 2.4
      P: IBA + ACC vs lanolin 0.046 0.350 0.038 0.014 0.042 0.089 0.500 0.196 0.091 0.180 0.390 0.030 0.201 0.160 0.442
      P: IBA + ACC vs ACC 0.061 0.219 0.244 0.014 0.052 0.368 0.353 0.091 0.091 0.139 0.274 0.030 0.500 0.244 0.318
      P: IBA + ACC vs IAA 0.070 0.425 0.178 0.229 0.307 0.500 0.233 0.196 0.375 0.433 0.139 0.500 0.034 0.065 0.395
      P: IBA + ACC vs IBA 0.070 0.435 0.390 0.050 0.200 0.333 0.151 0.008 0.149 0.416 0.061 0.500 0.365 0.040 0.045
      P: IBA + ACC vs IAA + ACC 0.030 0.323 0.095 0.071 0.261 0.024 0.500 0.008 0.500 0.500 0.061 0.003 0.500 0.384 0.181
      IAA + IBA
      Replicate Ph CZ RE SL Total Ph CZ RE SL Total Ph CZ RE SL Total
      1 2 2 4 2 10 4 4 1 0 9 5 6 4 0 15
      2 2 2 2 2 8 3 6 1 2 12 1 4 4 0 9
      3 2 2 2 2 8 3 8 1 3 15 3 8 1 3 15
      4 3 2 2 1 8 2 5 3 2 12 2 6 1 1 10
      Mean 2.3 2.0 2.5 1.8 8.5 3.0 5.8 1.5 1.8 12.0 2.8 6.0 2.5 1.0 12.3
      Std. Dev. 0.5 0.0 1.0 0.5 1.0 0.8 1.7 1.0 1.3 2.4 1.7 1.6 1.7 1.4 3.2
      P: IAA + IBA vs lanolin 0.065 0.005 0.010 0.003 0.423 0.067 0.178 0.122 0.034 0.066 0.178 0.123 0.124 0.293 0.174
      P: IAA + IBA vs ACC 0.009 0.001 0.061 0.003 0.160 0.271 0.235 0.073 0.034 0.045 0.228 0.123 0.229 0.390 0.225
      P: IAA + IBA vs IAA 0.140 0.007 0.405 0.015 0.028 0.376 0.095 0.338 0.276 0.129 0.402 0.053 0.242 0.189 0.118
      P: IAA + IBA vs IBA 0.140 0.018 0.293 0.269 0.102 0.268 0.065 0.122 0.046 0.105 0.398 0.030 0.183 0.126 0.019
      P: IAA + IBA vs IAA + ACC 0.089 0.017 0.365 0.500 0.001 0.067 0.178 0.122 0.207 0.134 0.398 0.324 0.269 0.310 0.500
      P: IAA + IBA vs IBA + ACC 0.012 0.030 0.160 0.067 0.040 0.312 0.178 0.196 0.207 0.101 0.152 0.014 0.242 0.390 0.155
      IAA + IBA + ACC
      Replicate Ph CZ RE SL Total Ph CZ RE SL Total Ph CZ RE SL Total
      1 2 5 2 0 9 5 6 1 1 13 3 6 1 5 15
      2 5 8 1 3 17 4 2 2 0 8 4 4 1 2 11
      3 3 6 3 0 12 5 3 3 0 11 4 8 1 3 16
      4 4 2 2 3 11 3 4 3 0 10 4 2 2 1 9
      Mean 3.5 5.3 2.0 1.5 12.3 4.3 3.8 2.3 0.3 10.5 3.8 5.0 1.3 2.8 12.8
      Std. Dev. 1.3 2.5 0.8 1.7 3.4 1.0 1.7 1.0 0.5 2.1 0.5 2.6 0.5 1.7 3.3
      P:IAA + IBA + ACC vs lanolin 0.384 0.364 0.014 0.124 0.072 0.353 0.178 0.022 0.196 0.207 0.500 0.433 0.500 0.037 0.131
      P: IAA + IBA + ACC vs ACC 0.500 0.243 0.140 0.124 0.088 0.195 0.126 0.013 0.196 0.142 0.269 0.433 0.104 0.050 0.168
      P: IAA + IBA + ACC vs IAA 0.405 0.432 0.235 0.124 0.283 0.128 0.420 0.065 0.065 0.336 0.089 0.158 0.008 0.028 0.089
      P: IAA + IBA + ACC vs IBA 0.405 0.440 0.500 0.500 0.500 0.341 0.426 0.500 0.500 0.313 0.009 0.135 0.338 0.024 0.016
      P: IAA + IBA + ACC vs IAA + ACC 0.271 0.333 0.140 0.406 0.266 0.006 0.178 0.500 0.149 0.375 0.009 0.176 0.282 0.145 0.424
      P: IAA + IBA + ACC vs IBA + ACC 0.084 0.500 0.353 0.077 0.179 0.065 0.178 0.040 0.149 0.339 0.378 0.114 0.201 0.095 0.117
      P: IAA + IBA + ACC vs IAA + IBA 0.074 0.040 0.235 0.399 0.056 0.048 0.074 0.160 0.046 0.194 0.166 0.271 0.124 0.084 0.418
      Legend (t tests): values in green, p ≤ 0.05; values in yellow, p = 0.05−0.10; values without any color, p > 0.10.

      In stem segments of all treatments, axial resin ducts circumscribed by tiers of sheath cells were present (Fig. 3eg) throughout the cortex. As already noted for dormant stems, the duct opening appeared to be at several stages of ontogeny from phloem to epidermis.

      The cortex in stem segments treated with lanolin (Fig. 3a) or ACC (Fig. 3b) appeared similar and unchanged from cortex anatomy as viewed at the start of the investigation. It appears possible that ACC somewhat enhanced cortical aerenchyma development and altered secondary wall thickness and chemistry of pre-existing parenchyma (Fig. 3b); however, those aspects remain for future investigation.

      In swollen apical regions of stem segments treated with IAA, cortical ducts became greatly stretched radially, possibly in response to radial enlargement of adjoining parenchyma connected through cell-wall bonds to the ducts. Many cortical parenchyma cells underwent expansive primary-wall growth in response to IAA treatment, but no evidence was found for development of new resin ducts (Figs 2d, 3c).

      In swollen apical regions of stem segments treated with IBA, preexisting cortical tissue became altered into a more diffuse tissue comprising exceptionally large primary-walled cells with numerous intercellular gaps (Figs 2e, 3d). A pinkish-colored reaction of primary cell walls to toluidine blue (Fig. 3d) is evidence of recently produced cell walls, probably due to strong radial expansion in pre-existing cell walls that enveloped small diameter parenchyma cells. In two of the four examined stem segments treated with IBA, the cortex displayed a phloem-enveloping zone of enlarged cortical resin ducts (see also Fig. 2e, h), but as can be seen in Fig. 3d this zone was not invariably present, and therefore it could not be decided if they were newly formed in response to the IBA treatment.

      Figure 4 shows several aspects of cortical resin duct formation as observed in stem segments that served as controls to the phytohormone treatments. There was abundant evidence for cell division, expansion, making and breaking of intercellular bonds, and cellular differentiation (Fig. 4ag). Collapsed cells and intercellular spaces were present throughout (Fig. 4ae). Plasmolysis occurred in thin-walled parenchyma (Fig. 4a, b), but thereafter those collapsed cells developed thickened walls (Fig. 4d). Non-collapsed primary-walled parenchyma in the vicinity of collapsed cells became greatly enlarged, then separated schizogeneously (Fig. 4ac), followed within developing ducts by what appeared to be lysigeneous bursting (Fig. 4cg).

      Circumferential tiers of parenchyma cells were produced around the duct by control of the cell-division plane (Fig. 4f). Cells at the exterior surface of the duct sheath enlarged and produced thickened secondary cell walls in the tiers of the sheath border (Fig. 4eg). With advanced duct development, those tiers of circumferential sheath cells changed biochemically (compare Fig. 3a & g), presumably in support of resin formation and secretion.

      In addition to resin ducts, cortical woody ducts were present as a novel anatomical feature (Fig. 5ae). Examination of cross sections of 26 randomly selected distantly spaced stem positions from eight stem segments fixed on the starting day of the experiment yielded only nine woody ducts in a total cortical resin duct population estimated at > 1,500. After the March 19th conclusion of the investigation, woody ducts were found present in lanolin control and all hormone treatments, but no convincing evidence was observed for any experimental treatment having either increased their frequency or altered their diameter.

      Additional investigation of dormant stem segments collected and fixed in winter revealed an invariable spatial association of cortical woody ducts with mature leaf attachment points to the stem (Fig. 5i). Depending on their stage of development, woody duct transverse diameters were determined to be between 80 μm and 200 μm. Investigation of dormant stem segments collected and fixed in winter revealed an invariable and obligate spatial association of cortical woody ducts with mature leaf attachment points to the stem (Fig. 5i). Woody ducts extended axially through the cortex for at least 1 mm (Fig. 5g, h). Over that longitudinal distance, some portions of the woody duct displayed non-woody parenchyma while other axial positions in the same duct displayed woody elements (Fig. 5g), possibly explainable in terms of varied extents of woody duct maturation at different axial positions (Fig. 5ae).

      Cortical woody ducts basal to mature leaves were initiated through similar processes described above for cortical resin ducts (Fig. 4), but as the woody duct channel opened up it became occupied by extending parenchyma cells (Fig. 5a, b, f) that subsequently divided to create ordered radial files of elongated small-diameter cells within the duct followed by the production of annular ribs of lignified secondary walls similar to those of primary xylem (Fig. 5be, h). No bordered pits were seen in the woody elements. Ray-like nucleated parenchyma cells, uniseriate and rarely bi- or tri-seriate, sub-divided the radial files of woody elements of each cortical woody duct; these duct-bisecting parenchyma evidently emerged early in woody duct formation (Fig. 5ac). As with vascular tissues, files of woody elements displayed radial polarity.

      Also associated with cortical woody ducts and not with other resin ducts of the cortex was an arc of cambium-like cells having elongated nuclei. These cambial cells arose within or near the radially external border of the outermost tier of the sheath cells surrounding the developing woody duct (Fig. 5be, h). Formation of the woody-duct cambium was not essential for the cells that differentiated into woody elements within the duct to be produced, as the cambium formed after, rather than before, the appearance of the duct's woody cells (Fig. 5a, b).

    • Supplemental Table S1 provides radial file cell number data for dormant stem segments at the start of the experiment.

      After the 28-d experimental treatment, no evidence was found for traumatic resin canals (TRC) on the inner (centripetal) side of the vascular cambium having been produced as a response to any hormone treatment. Debudded stem segments having intact leaves also displayed no TRCs. In contrast, segments having only intact buds and those having both buds and leaves produced readily detectable TRCs in the first formed earlywood.

      Table 1 presents each of the four analyzed stem segments at its three examined positions cell counts per radial file (averages based on six radial files) of the several cell types observed. Means, standard deviations, and p values are provided. Green boxes indicate results whereby the Welch t test it is 95% or more probable that the null hypothesis should be rejected. Yellow boxes in Table 1 have p values between 0.05 and 0.10, too high to reject the null hypothesis but possibly an indication that if a larger sample size had been tested, rejection may have been indicated.

      Tests comparing stem-segment positions and cell types support the interpretation of qualitative differences seen during microscopy examinations. IAA and IBA promoted some xylogenesis in the apical portions of the stem segments, but the response at lower positions was limited to radial expansion of primary-walled cambial derivatives; they did not become TEs despite the 28-d incubation period. TEs differentiated at all three stem-segment positions in only one of the four segments provided IAA, and in none of those provided IBA.

      IAA + IBA in combination yielded increased radial file cell numbers in both the CZ and the zone of earlywood TEs, particularly in the apical region of stem segments; however, only two of the four segments produced new TEs in the basal region. A similar result was produced by the IAA + ACC combination. In contrast, the IAA + IBA + ACC treatment favored xylogenesis primarily in the basal region of stem segments, with very little TE production occurring above (Table 1).

      Examples of stem-segment responses to lanolin, IAA, IBA and a combination of IAA + IBA are provided in Fig. 6. Protoplasm of vascular cambium in segments receiving lanolin alone (Fig. 6ac) cleared from its previously condensed dormant state, but otherwise there was no change in three of the four stem segments. The fourth lanolin control segment produced one new TE per radial file in two locations, only, around the entire circumference of the apical region, and two TEs per radial file were present but only at one circumferential location in the basal region (see Fig. 6c inset). In the mid-stem of that same segment, the cambial zone remained unproductive and adjoining the latewood.

      IAA had its greater TE-inducing effect at apical ends (Fig. 6d), weaker at mid-stem regions (Fig. 6e) and absent in three of the four segments at basal ends (Fig. 6f, Table 1). IBA also promoted TE differentiation at apical ends (Fig. 4g) but not at lower stem-segment positions (Fig. 6h, i, Table 1). The response to an IAA + IBA combination in two of the four stem segments was similar to that of IBA alone (Fig. 6jl); however, the other two stem segments displayed new xylem at each of the three examined positions. IAA + IBA clearly had a positive effect on the number of cells per radial file in the cambial zone (Table 1).

      Figure 7 provides examples of how stem segments responded to auxin + ACC mixtures. IAA + ACC elicited a relatively strong xylogenic response at apical ends (Fig. 7a), weaker at mid-stems (Fig. 7b) and strongest at stem bases in three of the segments (Table 1), but the fourth failed to respond in a like manner (Fig. 7c). IBA + ACC produced xylogenic results similar to those of IAA + ACC (Fig. 7df; Table 1). The triplet combination of IAA + IBA + ACC promoted xylogenesis in basal locations but in only one of the four mid-stem positions and in only two of the four apical regions (Fig. 7gi; Table 1).

      It was observed in radial sections of the vascular cambium region that ACC, both by itself and in combination with an auxin, enhanced the dark coloration of coarse cell-wall-adjoining lines visible between axial and radial elements. These mostly horizontal radial 'lines' are intercellular spaces between procumbent ray cell walls (Fig. 8a, b). Although narrow, the spaces are evident when viewed in the SEM (Fig. 8c).

      Figure 8. 

      Sections of balsam-fir cambium and xylem that display variable darkening of intercellular spaces. (a) Brightfield radial section of cambial zone (cz) bordering latewood (LW) and stained with toluidine blue. The white arrow indicates a darkened compound middle lamella region and the black arrow an absence of similar darkening. (b) The white arrows indicate intercellular spaces between procumbent xylem ray cells that stained less intensively to toluidine blue. (c) SEM of a xylem ray in radial section; the arrows point to intercellar spaces between the procumbent ray cells. (d) Tangential section showing two rays on the centripetal periphery of the cambial zone. The arrow points to an intercellular space between the radial wall of a fusiform cell and that of a ray cell. (e) Tangential section of a 2-celled ray in the cambial zone, the small arrow pointing to evidence for an intercellular space between fusiform and ray cell walls. An axially oriented 'spear tip' (large arrow) appears to contain particulate matter and to be intruding between what had been adjoining walls of two fusiform cells. (f) Radial section (interference contrast optics) of phloem (Ph), cambial zone (CZ), radially expanding cambial derivative (RE) and differentiating TEs in proximity to latewood (LW). The arrows point to accumulations of insoluble matter within or paralleling the axial walls.

      The dark coloration within intercellular spaces when viewed by brightfield microscopy extended radially along procumbent rays from mature xylem through cambial zone into mature phloem, although in the cambial zone the staining reaction to toluidine blue, a metachromatic dye, was different from that in phloem and xylem, indicating chemical differences (Fig. 8a). Figure 8d & e are tangential sections viewed by differential interference contrast optics and showing the triangular shapes of the intercellular spaces in the cambial zone. The large arrow in Fig. 8e points to a spearhead-shaped intrusion between walls of two fusiform cells, perhaps related to the deposits of material associated with the axial walls of fusiform cells at an early stage of TE differentiation (Fig. 8f).

      Data for those stem segments cultured in this experiment but not treated with a lanolin paste are provided in Supplemental Table S1 and Supplemental Fig. S1.

    • The primary focus of this investigation was to discover if qualitative differences in tissue anatomy might emerge, following applications of ACC, IAA, IBA, and their combinations to the apical ends of balsam-fir stem segments. The working hypothesis based on earlier research was that varied expressions of xylogenesis, e.g., formation of primary vs secondary xylem TEs should be visible after those hormones were applied at identical concentrations dispersed in lanolin. Unexpectedly, differences in cortical development in response to hormone treatment were also found. In addition, several cellular phenotypes, possibly entirely novel, also were observed as natural phenomena of non-treated stems.

      The concept of vascular development being regulated principally by auxin originated in herbaceous species[37] and, when applied to conifers, has limitations depending on the age and stage of development of the stem segments investigated[3]. Previous investigations with several conifer species revealed that IAA was entirely ineffective in promoting xylogenesis when stem segments older than two years were subjected to the same treatment that stimulated xylogenesis in young stems[38,43]. For this reason, only stem segments prepared from leaders that had grown in the preceding year were investigated.

      At the start, the cambium in the investigated stem segments was fully dormant; the stems lacked any evidence of cell division or any stage of earlywood formation. Under the bioassay's culture conditions, a four-week incubation period was assumed to be more than ample for cambial reactivation, cell production and completion of cellular differentiation to occur. After stopping the experiment, that assumption was confirmed by observations that earlywood tracheids had matured.

      Because all segments receiving hormone treatment were first debudded and defoliated, lacked roots, and were provided only water, the metabolism underlying the observed growth and development responses can be assumed to have drawn upon storage reserves within the cortex, ray, and pith parenchyma cells. Remarkably, almost all segments appeared still healthy after 28 d, and the few having any evidence of dead cambium were rejected.

      Although the experimental timeframe proved adequate for newly produced cambial derivatives to become fully differentiated, evidence for TE differentiation was present at only some, not all, of the three examined stem-segment positions, even in treatments where xylogenesis appeared to be more strongly induced. Within stem segments displaying earlywood TEs, examination of three positions (apical, basal, and mid-stem) revealed other positions where no TEs had been produced anywhere around the stem circumference. At those positions, centripetal cambial derivatives were enlarged but remained primary walled and lacked evidence for initiation of bordered-pit development, the earliest stage indicating cellular commitment to TE differentiation[40].

      In relation to phloegenesis, similar radially enlarged cambial derivatives accumulated centrifugal to the cambial zone, presumably as nascent phloem cells. Many of those derivatives on the cambium's centrifugal side were primary-walled, not fully differentiated as either sieve cells or other phloem elements. This incomplete development of phloem was particularly apparent in response to IBA treatments.

      Xylogenesis occurred in response to IAA, IBA, and IAA + IBA, but compared to cortical growth responses at the apical ends of those stem segments, all xylogenic responses were relatively weak. Thus, no conclusions are offered about which hormonal treatment was more effective based on radial file cell number counts, mainly because cortical development appeared to be a more favored response but also given the heterogeneous anatomical responses associated with the vascular cambium and its derivative cells. The p values (Table 1, Supplemental Table S1) are provided merely as indications of possibly productive avenues of future research.

      Observations on the presence/absence of cell types, in particular of TEs, appear to be more important than any deductions that might be inferred based on quantitative data. For example, in two of the four stem segments treated with IAA + ACC, no TEs were produced at their mid-stem region despite apical and basal regions having produced new TEs (Table 1). Earlier papers concerned with hormonal effects in the balsam-fir stem segment bioassays have overlooked this natural variation in responsiveness, perhaps because only single rather than multiple positions in stem segments were observed.

      The leaders used to produce stem segments in this investigation were from a single population of sapling trees growing nearby one another on the same site, and all leaders were selected based on being overtly healthy, straight, and of zenith orientation. Based on their past growth performance, it can be imagined that all would have grown similarly strongly in both height and diameter had they remained intact on their source trees and received contributions from roots, leaves, and buds. But again, the three examined positions within each stem segment produced different responses, even those within the lanolin control, and variation among the four replicates further corroborated this intrinsic variability.

      Longitudinal variation over lengths of tree trunks in springtime cambial reactivation, the onset of cambial dormancy, and rates of xylem development are established natural processes[3,37,38]. With current knowledge, somatic genomic variation at different positions within the stem of the tree cannot be confidently excluded as a plausible explanation for such variation[3,44]. However, it seems more probable that the unequal variance revealed by microscopy in this study could be explainable in terms of intrinsic metabolic differences between stem positions and between genotypes, for example in their storage and endogenous hormone reserves, the status of dormancy release or other physiogenetic considerations, in particular, mechanisms of uptake and internal transmission of exogenous hormones.

      Differences between stem positions in their xylogenic responsiveness to auxin were earlier noted in other conifer species[32,33]. However, those qualitatively different developmental responses had no explanation and were relegated to the realm of artifacts. The observations made in this investigation negate an artifactual interpretation. The relatively small diameters of balsam-fir stem segments enabled complete cross-sectional examination of the cambial region around entire stem circumferences in the microscope. Moreover, the bioassayed stem segments when first collected as leaders were all of the zenith orientation, and the stem segments were subsequently maintained in their vertical orientation throughout the investigation. Therefore, it is assumed that the well-known effect of gravity on inducing circumferentially unequal cambial growth responses was not a complicating factor in this investigation.

      The inconsistency of the cambial response over the length of the stem segment seems to indicate that the cambium competence for xylogenesis varies, in terms of the intrinsic ability of cambial cells to receive, transmit, or respond to hormones. The observed heterogeneity indicates the individuality and complexity of the numerous underlying physiological processes that function during cambial growth and xylogenesis[3]. Before initiating further experiments, there is a need to know how to better characterize, screen, and standardize the intrinsic variation existing within the specimens selected for investigation. When debudded-defoliated stem segments from leaders are intended as the bioassay specimens, it can be recommended that measures be taken to ensure that all initially are investigated and confirmed to have identical apical and basal diameters and equal numbers and distributions of leaves and buds. In reality, achieving such standardization will not be easy, because casual examination reveals that each individual leader on a tree grows somewhat uniquely from those of neighboring trees. As already noted, there is also need during microscopy to observe full stem circumferences at apical, mid-stem and basal stem-segment positions.

      Earlier investigations into balsam-fir stem-segment bioassay responses to hormones, auxin, in particular, have retained mature leaves intact on the stem segments. Young developing leaves have been described as the primary source of IAA in balsam-fir trees[23], but investigations with other conifer species indicated that mature leaves and dormant buds both contributed IAA to cambium[3,23,32,45]. In pine species, IAA was exported from mature leaves during winter dormancy as well as during the growing season[45]. Past research has also revealed that mature conifer leaves when intact on stem segments effectively promoted xylogenesis[32,38,43,4648], and the same was again observed in this investigation (Supplemental Table S1).

      In the present investigation, buds left on stem segments whether or not defoliated remained dormant throughout the trial period, and cambium nevertheless reactivated and produced a small amount of new xylem. None of the budbreak, shoot elongation, or new leaf development was necessary for the resumption of diameter growth. The regulation of dormancy release in apical and lateral meristems must differ, as the vascular cambium in stem segments reactivated while the shoot apical meristems remained dormant. Similar findings have been noted in other conifer species[37,38].

      It is noteworthy that traumatic resin canals were produced in spatiotemporal association with xylogenesis in non-experimental stem segments having buds, and TRCs were not produced in debudded-defoliated stem segments subjected to any of the hormone treatments, despite the razor-blade wounding that was done to remove leaves and buds. Two plausible explanations can be suggested, that buds export a TRC-inducing factor (such as methyljasmonic acid) or, alternatively, that the lanolin carrier used to provide hormones to segments either blocked or absorbed the TRC-inducing factor.

      In addition to IAA, mature leaves export many other compounds, possibly including ACC for which there is growing evidence for a role in xylogenesis[3]. Exogenous ACC applied by itself to needle stumps on stem segments resulted in xylogenesis independently of preceding cambial cell division or radial expansion of cambial derivatives[43]. The inductive effect of ACC was observed to be similar to that of mature leaves, resulting in both cambial fusiform and cambial ray cells transforming into tracheary elements, mostly of primary-xylem cell types having altered fine structure[48] and lacking the bordered pits characteristic of secondary-xylem tracheids[40].

      In another investigation, an IAA + ACC combination was more effective than IAA alone at inducing cambial cell division, radial enlargement and differentiation of cambial derivatives into new tracheids[43]. Similar complementary and seemingly synergistic effects of auxin and ACC were noted in relation to xylogenesis in Armoracia rusticana roots cultivated in vitro[49].

      In the present investigation, no convincing evidence was found for applied ACC having induced xylogenesis independently, though it appeared to enhance the promotion induced by both IAA and IBA. Plausibly, the present observations can be reconciled with earlier findings on the basis that they showed that little if any ACC oxidation to ethylene occurred when ACC was supplied to cambium cells cultured in vitro except when IAA was also provided[25]. Because the present investigation began with debudded and defoliated dormant stem segments, it may be that insufficient endogenous IAA was present to promote ethylene production.

      Segments having mature leaves, only, yielded the stronger xylogenic response, but only at the basal stem-segment location. Of the several hormone treatments investigated, debudded-defoliated stem segments treated with an IAA + IBA + ACC combination responded most similarly to segments having intact mature leaves (compare Table 1 with Supplemental Table S1). Investigations to quantify how much IAA, IBA, and ACC are present in mature balsam-fir leaves and allocated to vascular development remains to be done.

      Small intercellular spaces that are oriented radially between tiers of procumbent ray cells have frequently been considered as transverse pathways for gas exchange[5052]. This investigation has provided evidence that those pathways have continuity from xylem to non-collapsed phloem, possibly extending even into the cortex. It could not be conclusively demonstrated but seems reasonable to suspect that aqueous fluids - those continually present in the secondary xylem, vascular cambium and secondary phloem - would diffuse into those intercellular spaces and traverse those tissues. Some intercellular spaces stained with toluidine blue while others displayed darkened contents. Based on this and observations made in an earlier investigation[43], it can be suggested that these intercellular spaces between procumbent ray cell walls do more than transfer gases radially.

      Excepting reactivation of dormant cambium and subsequent development of vascular tissues[23,24,2731], little knowledge yet exists about how overall stem development occurs during the balsam-fir leader's second year of life. Advances have been made in understanding post-cortical stages of periderm formation in fir trees[5357]. For example, large-diameter parenchyma cells in the bark of Abies firma and A. homolepis, were observed long ago[53,54]. However, descriptions of the cortex associated with the first two years of a fir tree stem's life generally have been brief, portraying a relatively simple tissue system. A prevalent but assumptive concept is that epidermal, cortical, vascular, and pith development of stems during their first year arise as largely predetermined structures within the dormant shoot apical meristem.

      This study revealed the balsam-fir leader cortex to be a complex milieu of different cell types, shapes, and evident specializations, an impressive diversity of mostly living cells capable of growing and differentiating. Such complex histogenesis requires regulation; it cannot simply arise through predetermined developments in the apical meristem.

      The promotion by auxin treatments of cortical diameter growth, viewed macroscopically as stem swelling, was particularly visible for approximately 3 cm below the apical ends of stem segments (Fig. 1). The cellular expansion activity underlying swelling of the cortex indicates that both IAA and IBA can move through cortical tissues, in the process causing small-diameter primary-walled parenchyma to enlarge into quite large cells. This is an important observation worthy of further research, as there has been a tendency to assume that auxin transmission in stems is restricted to cambium and phloem. The axial continuity and zonation of enlarged cortical parenchyma suggests the possibility that a signal translocated intercellularly through the tissue while still at only the brick stage is the explanation for other cells extending or enlarging radially[3,38,48].

      Previous investigations into young arabidopsis seedlings revealed IAA and IBA to move preferentially in the basipetal direction through the hypocotyl, both at similar transport rates (8−10 mm·h−1)[9]. Nevertheless, the transmission of IAA evidently was facilitated by a different protein complex than that of IBA[9]. The arabidopsis hypocotyl is anatomically, morphologically and evolutionarily distant from balsam-fir stem segments, but given the observed differences in cortical responses of balsam-fir stem segments to IAA and IBA, it seems possible that different mechanisms of transmission of those two auxins through cortex might explain the observed differences in cortical development.

      On the other hand, IBA has been described as merely a precursor of IAA in arabidopsis and other eudicots[11,13,14]. That concept is difficult to reconcile with the distinguishably different growth responses of the cortex of balsam fir in response to IAA and IBA. This seems to be evidence that those two endogenous auxins may not only have different transport requirements but also different regulatory effects in balsam-fir stems.

      Some cortical parenchyma responded to IAA and IBA by dividing and enlarging to increase the girth of the cortex. In the case of IAA, radial expansion of parenchyma cells resulted in distortion of preexisting resin-duct shape, as viewed in transverse sections, from rounded to radially extended, but this shape change was not attended by evidence for enlargement of the ducts. The change from a circular to radially extended shape indicates cell-wall bonding between the outer tiers of the resin duct and the adjoining parenchyma, although all around were aerenchyma and intercellular spaces. In the case of IBA, enlarged parenchyma arose proximal to the outermost phloem and gave the impression of being nascent ducts. Other primary-walled parenchyma in the same region and further out in the cortex displayed radial extension growth.

      ACC did not appear to affect the anatomy of the cortex, but it did appear to alter parenchyma cell wall thickness and staining reactions to toluidine blue; however, more intensive research is needed to corroborate these provisional interpretations.

      Sites of cortical resin duct formation appeared to be where numbers of collapsed cells among enlarged parenchyma were greater. Concomitant plasmolysis and turgid expansion in neighboring cells of the cortical tissue system occurs. If this is the explanation, the aqueous environment is not hypertonic for the general cellular population. Plasmolysis-susceptible cells either have inadequate osmotica or are compromised in their ability to regulate osmotic pressure, relative to their neighbors which, by their swelling, evidence a hypotonic environment. Cellular collapse augments intercellular space formation, but how such changes in the cortex may trigger non-burst living parenchyma cells to produce tiers of sheath cells surrounding enlarged parenchyma cells remain unclear.

      IAA has been reported to be a signal for periderm formation in arabidopsis roots[58], but no evidence was seen for that in this investigation. More research is needed, but based on observed cortical duct development, it could be that IBA is part of the explanation. Many additional signaling mechanisms, including wounding and pressure, deserve consideration[5961].

      Sites displaying cortical resin ducts were mainly within longitudinal zones of enlarged parenchyma. However, ducts were absent from similar zones where enlarged parenchyma appeared as schizogeneous aerenchyma[62]. As thin-walled parenchyma cells enlarge, some burst thereby producing intercellular spaces in support of aerenchyma. The explanation for bursting could not be resolved; it may have occurred through catalyzed lysigenous development, or it may have resulted from simple physical rupturing of primary cell walls enriched in pectin but deficient in constraining cellulose microfibrils.

      Excepting a study of the helical distribution of procambial leaf primordia in the developing cortex a few mm basal to the apical meristem[63], no description of how leaf traces develop and persist in Abies spp. could be found in the botanical literature (an exhaustive search was not possible). The development as observed in this investigation indicates that the leaf trace begins unusually, as a cortical resin duct lacking procambium or other vascular tissue. The cortical duct channel is invaded by intrusively growing axial parenchyma cells that subsequently subdivide, by repeated anticlinal and periclinal divisions, to populate the duct with well-ordered radial files of small diameter cells, most of which differentiate into woody elements having annular ribs, similar to those of the primary xylem. Before the onset of this xylogenesis, it may be that phloem translocatory cells differentiate around the perimeter of the duct; living cells are abundantly present around the developing woody duct; however, none could be definitively identified as sieve cells (see Fig. 5). Viewed in cross sections, each cortical woody duct is bisected radially by one or two radial files of enlarged thin-walled parenchyma and, as previously noted in relation to general gišogenesis, those ray cells may fulfill an essential role in signal transduction in support of vascular development[3].

      Earlier findings indicate that cambium formation is induced in response to the polar transport of auxin and intercellular compression[37], and after a cambium forms its maintenance in terms of fusiform cell shape and length require an ongoing supply of auxin[33]. However, auxin also is assumed to initiate xylogenesis, but this investigation revealed that xylogenesis within cortical woody ducts occurred before, evidently independently of, subsequent formation of the woody duct's cambium. Despite the other observed effects on cellular differentiation within the cortex, no evidence was found for any of exogenous ACC, IAA, and/or IBA inducing cambium to form or altering any aspect of cortical woody duct anatomy or formation.

      Examination of 3-year-old balsam-fir stem regions revealed that the overall anatomical structure of woody ducts in the periderm remained unchanged from that of the cortical woody duct originally produced during leader growth (Fig. 9).

      Figure 9. 

      Hand-cut sections of untreated balsam-fir dormant stems: (a), (c) and (e) are unstained sections; (b), (d), (f) and (g) are stained with toluidine blue. (a) Radial section of dormant cortex showing the abscission zone (a) and axially oriented leaf trace (t) of the mature first-year leaf (l), also showing a resin duct (d), latewood (LW) and vascular cambium (c). (b) Cross section basal to the leaf abscission zone (a) of the leader showing the cortical leaf trace (t); the inset at higher magnification shows the trace as a woody duct. (c) An unstained cross section at midway along the length of a 3-year-old balsam fir needle, showing its central vascular strand (s). (d) Tangential section of a 3-year-old balsam fir needle base at its stem attachment point showing the singular leaf trace (t) entering the stem. (e) Radial section of a 3-year-old balsam fir stem showing the vascular strand in the leaf (s) passing through the abscission zone (a) into the woody duct leaf trace (t) in the periderm. The trace traverses the cambium (c) and runs radially through three annual layers of secondary xylem (arrow) to the pith (p). (f) Cross section through the trace-cambium junction of a 3-year-old balsam fir stem, with vascular cambium (c) and latewood (LW) also indicated. Note that most leaf-trace cells are radially elongated but thin walled. (g) Cross section through latewood (LW) of the second year and earlywood (EW) of the third year in a 3-year-old balsam fir stem showing the size and anatomy of the leaf trace. (h) Cross section through second year earlywood in a 3-year-old balsam fir stem; the trace comprises mostly tracheids but also has living cells (arrow) that appear to be sieve elements.

      While in peridermal tissue, the woody duct's cambium remains dormant, unchanged from the structure produced in the cortex. Assuming that the woody duct cambium remains capable of reactivation, its persistence as a dormant peridermal meristem may explain how preventitious buds sometimes form in Abies bark[64].

      In contrast, at the junction of the leaf trace with the vascular cambium, the evidence indicates that the woody duct cambium becomes meristematically active, followed by radial cell elongation and differentiation of its derivative cells at the trace's interface with the secondary xylem to produce xylem, phloem and parenchyma elements of the trace (Fig. 9). The continuing radial growth of the trace in pace with that of secondary xylem presumably serves to provide the mature leaf with water and nutriment, and on this basis it can be suggested that activity of the woody duct cambium may be the primary basis for the well-known multi-year longevity of Abies leaves. However, this interpretation of meristematic activity by the woody duct cambium at the trace-cambium junctions seems to be logically flawed on the basis that IAA when transported basipetally via the leaf trace is nevertheless without effect on living cells in the cortex or periderm but then promotes its growth at the trace-cambium junction. On the other hand, serial transverse sections prepared over the axial length of the woody duct leaf trace in the cortex revealed that some places were entirely lacking of tracheary elements within the duct channel, which seems to indicate variable transmission or reception of the xylogenesis signal. Interruptions in primary xylem continuity over the length of the cortical trace also indicates that conduction of water and nutriment into the leaf from the secondary xylem of the stem must be principally through living cells circumscribing the duct.

    • This investigation revealed that isolated dormant stem segments lacking buds, leaves, and roots are competent to produce a variety of cell types in response to applications of micromolar concentrations of indole-3-acetic acid (IAA), indole-3-butyric acid (IBA), 1-aminocyclopropane-1-carboxylic acid (ACC) and combinations of those three. Those three hormones in combination appeared to mimic the overall effect of mature leaves left on phytohormone-untreated stem segments in their regulation of cambial growth, secondary phloem, and xylem development. IAA, IBA, and ACC individually yielded some xylogenesis, and ACC with IAA or IBA appeared to enhance the response.

      This study provides new information about the stem cortex, its diverse cell types and their varied sizes, shapes, wall thicknesses, polarities, and probable specializations. Cortical growth responses below hormone application sites were stronger than those of vascular development. IBA promoted the early stages of cortical resin duct formation, whereas IAA promoted radial expansion of existing resin ducts, and both promoted radial elongation of parenchyma. The cortical woody duct, formation of which precedes leaf trace development and is dependent on a factor from mature leaves, is an entirely novel anatomical observation deserving further research.

    • The author confirms sole responsibility for the following: study conception and design, data collection, analysis and interpretation of results, and manuscript preparation.

    • All data generated or analyzed during this study are included in this published article and its supplementary information files. References 32, 37, 43 and 49 lack a doi but are available at www.researchgate.net/profile/Rodney-Savidge.

    • This research was independently financed by the author. University of New Brunswick Libraries made literature access possible.

      • The author declares that there is no conflict of interest.

      • Supplemental Table S1 Welch′s t test values p values for stem segments not treated with a lanolin paste.
      • Supplemental Fig. S1 Examples of the responses as they appeared on March 19th of stem segments having buds and/or leaves and not treated with a phytohormone. (a)−(c) Apical, mid‐stem and basal regions, respectively, of segments having both buds and leaves present. (d)−(f) Apical, mid‐stem and basal regions, respectively, of segment having only a whorl of terminal buds present.  (g)−(i) Apical, mid‐stem and basal regions of a segment having only leaves present. The inset images in (d), (g) and (i) show the appearance of the cambial zone and adjoining cells at those positions as they appeared in radial section.
      • Copyright: © 2024 by the author(s). Published by Maximum Academic Press, Fayetteville, GA. This article is an open access article distributed under Creative Commons Attribution License (CC BY 4.0), visit https://creativecommons.org/licenses/by/4.0/.
    Figure (9)  Table (1) References (64)
  • About this article
    Cite this article
    Savidge RA. 2024. Responses of isolated balsam-fir stem segments to exogenous ACC, IAA, and IBA. Forestry Research 4: e033 doi: 10.48130/forres-0024-0030
    Savidge RA. 2024. Responses of isolated balsam-fir stem segments to exogenous ACC, IAA, and IBA. Forestry Research 4: e033 doi: 10.48130/forres-0024-0030

Catalog

  • About this article

/

DownLoad:  Full-Size Img  PowerPoint
Return
Return