Alteration of TAC1 expression in Prunus species leads to pleiotropic shoot phenotypes

Prunus persica (peach) trees carrying the “Pillar” or “Broomy” trait (br) have vertically oriented branches caused by loss-of-function mutations in a gene called TILLER ANGLE CONTROL 1 (TAC1). TAC1 encodes a protein in the IGT gene family that includes LAZY1 and DEEPER ROOTING 1 (DRO1), which regulate lateral branch and root orientations, respectively. Here we found that some of the native TAC1 alleles in the hexaploid plum species Prunus domestica, which has a naturally more upright stature, contained a variable length trinucleotide repeat within the same exon 3 region previously found to be disrupted in pillar peach trees. RNAi silencing of TAC1 in plum resulted in trees with severely vertical branch orientations similar to those in pillar peaches but with an even narrower profile. In contrast, PpeTAC1 overexpression in plum led to trees with wider branch angles and more horizontal branch orientations. Pillar peach trees and transgenic plum lines exhibited pleiotropic phenotypes, including differences in trunk and branch diameter, stem growth, and twisting branch phenotypes. Expression profiling of pillar peach trees revealed differential expression of numerous genes associated with biotic and abiotic stress, hormone responses, plastids, reactive oxygen, secondary, and cell wall metabolism. Collectively, the data provide important clues for understanding TAC1 function and show that alteration of TAC1 expression may have broad applicability to agricultural and ornamental tree industries.


Introduction
Given the limited availability of agricultural land and water, increases in crop productivity will require higherdensity agricultural production to keep pace with the growing worldwide demand for food. During key stages of the domestication of cereal crops such as rice and maize, considerable gains in productivity were achieved through modifications in plant architecture. This included reduced tillering and breeding for upright tiller and leaf angles, which minimizes competition and increases the efficiency of light capture under crowded conditions 1,2 . There is a growing body of evidence that pillar or columnar shapes have the potential to increase crop productivity in fruit trees as they have for cereal crops. The "Broomy" (br) trait in peach was previously bred to produce commercial quality peach and nectarine varieties with a pillar or, in the case of heterozygous individuals, an upright stature [3][4][5] . The narrower canopy of these trees was shown to improve dry matter partitioning due to more effective light interception when the angle of the sun is less than vertical 6,7 . In addition to branch angle, pleotropic phenotypes have also been described in pillar peach trees, including fewer sylleptic branches, shorter branches, longer internodes, and increased auxin and auxin-tocytokinin ratios, features possibly associated with increased apical dominance 8 . A similar trait in apple called columnar or Co is also being bred for integration into high-density orchard systems 9 . Co trees also have narrow branch angles, fewer sylleptic shoots, and higher apical dominance 10 . But unlike pillar peaches, columnar apples produce numerous short spurs and have thickened stems 10 . The Co locus has been shown to contain a Ty3type transposable element that is not inserted within a known gene but may act, in part, through altered expression of a neighboring 2-oxoglutarate (2OG)-Fe(II) oxygenase 11,12 . The ultimate utility of pillar or columnar traits to improve productivity in orchards will require a greater understanding of how the underlying genes function and their roles in overall plant development.
To date, there is relatively little known about the genes and molecular mechanisms that contribute to changes in tree branch orientation and much of what is known comes from work in monocots. Three genes were previously identified in cereal crops that are associated with upright tiller and/or leaf angle, TILLER ANGLE CON-TROL 1 (TAC1) and two zinc finger proteins LOOSE PLANT ARCHITECTURE 1 (LPA1) and PROSTRATE GROWTH 1 (PROG1) 2,13,14 . LPA1 and its Arabidopsis ortholog SHOOT GRAVITROPISM 5 (SGR5) regulate shoot gravitropism upstream of amyloplast sedimentation 14,15 . No orthologs of rice PROG1 have yet been identified in dicots.
We previously showed that pillar traits in peach are caused by mutations within a TAC1 ortholog 16 . Changes in TAC1 have since been linked to upright tiller or branch angles in other plant species including Arabidopsis, Miscanthus sinensis, Brassica napus L. (rapeseed) and Populus species (poplar) [17][18][19] . TAC1 was found to be expressed to varying degrees in shoot tips, branch nodes, leaf axils, and flower buds in Arabidopsis, peach, and poplar 13,16,19,20 . In field-grown peach trees, TAC1 expression was seasonal, being highest in early spring and coinciding with lower auxin levels 21 .
In most plant genomes, TAC1 typically occurs as a single gene and putative TAC1 proteins contain no known domains or motifs nor have any molecular functions been ascribed to date 16 . However, TAC1 shares distant sequence similarity with other known plant architecture regulators including LAZY1 and DEEPER ROOTING 1 (DRO1), which together comprise a larger gene family dubbed IGT for a conserved amino acid triplet found within domain II of the protein 16,22,23 . In rice and Arabidopsis, lazy1 mutants display the opposite phenotype as tac1, having wider lateral tillers or branches, while dro1 mutants have wider lateral root angles [23][24][25][26][27][28] . In Arabidopsis, DRO family members were recently shown to also play a minor role in contributing to shoot architecture 28,29 . LAZY1 mutant phenotypes are associated with reduced gravitropic responses. These proteins appear to function downstream of statolith sedimentation but upstream of auxin translocation [24][25][26][28][29][30] , though no specific molecular or cellular functions are known at this time. The primary difference between TAC1 and LAZY/ DRO proteins is the presence of a C-terminal domain (domain V) containing a putative Ethylene-responsive element binding factor-associated Amphiphilic Repression motif, which is not present in TAC1 proteins 16,23,28 . A role for TAC1 in regulating gravitropism has not been established; however, based on the phenotypic, phylogenetic, and expression patterns of TAC1 and LAZY1, we hypothesized that TAC1 may alter shoot growth angle by negatively regulating LAZY1 16 .
Here we further characterized TAC1 in Prunus species to better understand both the functional role of this gene and the potential applicability to producing upright fruit trees for high-density orchard systems.

Peach and plum germplasm
The peach population used for structural measurements and expression profiling was derived from a selfpollinated tree heterozygous for the pillar trait (KV991636) in 2002. KV991636 was a seedling from a cross between Kv93065 (a pillar tree) as the female parent and pollen from the peach variety "Weeping White". Leaf samples from the plum variety "Improved French" were used for Illumina genome sequencing. For plum transformation and associated controls, seedlings derived from open-pollinated "Bluebyrd" or "President" trees were used.

Plum Illumina sequencing
Genomic DNA was extracted from leaves of "Improved French" plum trees as previously described 16 . Briefly, leaves were pulverized in liquid nitrogen and genomic DNA was extracted using the EZNA TM High Performance (HP) DNA Kit (Omega Bio-Tek Inc., http://www. omegabiotek.com) with the addition of 2% polyvinylpyrrolidone-40 (PVP-40) (w/v) to CPL buffer and 2-mercaptoethanol. Genomic DNA quantity was assessed using the Quant-iT PicoGreen Kit (Invitrogen, Carlsbad, CA).
Two micrograms of purified DNA was provided to David H Murdock Research Institute, Kannapolis, NC for library construction and sequencing. Both a paired-end and a mate-pair library were constructed with average insert sizes of 375 bp and 2950, respectively, and sequenced using an Illumina HiSeq 2000. A total of 194,856,870 100 bp paired-end reads and 158,319,386 mate-pair reads were obtained. Reads were assembled against a 5 kb TAC1 genomic region (Chr#2: 19,657,662,152) (peach genome version 1) using the CLC Genomics Workbench (Qiagen, USA) reference assembly tool with 2 modifications to the default settings: length fraction = 0.7, similarity fraction = 0.9 31 . The TAC1 genic region had an average coverage of~100 reads per base representing 14× coverage of the hexaploid plum genome. Variant detection was performed using CLC Genomics Workbench Basic Variant Detection tool with the following modifications to the default parameters: Ploidy = 6, minimum coverage = 20, minimum count = 4, minimum frequency = 8.0. Alignments of reads spanning the GAT repetitive element within exon 3 were extracted and manually refined using BioEdit (http://www.mbio. ncsu.edu/BioEdit/bioedit) to correct for errors and identify the six variant alleles.

Tree measurements
Peach tree branch measurements were taken from 10 standard trees (TAC/TAC1) and 10 pillar (tac1/tac1) at the end of the 2012 growing season. Fifty measurements in total were used for both 2012 branch growth and diameter, with five branches measured per tree. The 2012 shoot growth (length) measurements taken from June 11 to August 14 were all from the same branches. Trunk diameters and heights were taken from 50 individual standard and pillar trees. All peach trees were from 4-year-old trees in the field that were part of a pillar/ upright/standard (1:2:1) segregating population. Trunk and stem diameters were measured 4 cm from the base using calipers. Field-grown RNA interference (RNAi) and overexpression (OE) TAC1 transgenic plums were measured in the spring of 2017 as follows: trunk and branch diameters were measured using digital calipers, with trunk diameter taken at 4 inches from the base, and branch diameters taken at the base of the branch; branch diameters, lengths, and angles were taken for the first four major branches of each tree; branch angles were measured using a protractor aligned to the main trunk and reported as degrees from vertical. Transgenic plums were compared to untransformed seedlings from their respective standard plum varieties, "Bluebyrd" or "President". Measurements were taken for three of the RNAi lines (Lines 1, 3, and 6), all three OE lines. Between four and six replicate trees were measured per line and four branches per tree were measured for angles and lengths.

Tree models
Tree models were acquired using a robotic vision system called Robotic System for Tree Shape Estimation, or RoTSE, which is described in detail in Tabb and Medeiros (2017a) 32 . RoTSE consists of a small truck, robot arm, cameras, generator and other incidentals. The truck was parked near the tree, and a background unit consisting of blue background material was parked on the other side of the tree. The robot arm outfitted with cameras acquired still frame images. These images were processed via algorithms described in Tabb and Medeiros (2017a), Tabb (2013), and Tabb (2014) to generate estimates of the tree's shape, which was subsequently skeletonized with the method of Tabb and Medeiros (2017b) [32][33][34][35] .

Plasmid construction
To generate the TAC1-silencing vector, PpeTAC1-HG, PpeTAC1-specific primer sequences Fwd 5′-TGG GTTTGCTGGGAATGTGA-3′ and Rev 5′-CAGCTG GTTTCTGAACAATGGC-3′ were used to PCR amplify a 380-base-pair cDNA fragment from peach genomic DNA. The resulting fragment was cloned into the pENTR-D TOPO (Invitrogen) vector per the manufacturer's specifications and sequenced for verification. PpeTAC1_300 pENTR-D TOPO was recombined with the RNAisilencing vector pHellsgate 8 (Commonwealth Scientific and Industrial Research Organisation (CSIRO), Australia) using gateway technology to create PpeTAC1-HG (Invitrogen, Carlsbad, CA). To generate the 35S::TAC1 OE vector, PpeTAC1 primer sequences Fwd 5′-GAATT-CAATTCGCTCACAAAATATGAAG-3′ and Rev 5′-CC TTGTGTGCACTGAATTAAGGATCC-3′ were used to amplify the full-length PpeTAC1-coding sequence from peach RNA purified from apical shoots. The resulting fragment was digested and cloned into a modified pBIN-ARS vector (called pBIN-AFRS) using the EcoRI and BamHI restriction sites behind the Cauliflower Mosaic Virus 35S promoter (vector sequence and map provided in Datafile S1).

Plum transformation
PpeTAC1-HG and 35S::PpeTAC1 vectors were transformed into the Agrobacterium tumefaciens strain GV3101. Agrobacterium-mediated transformation of plum hypocotyl slices was performed as previously described in Petri et al. (2012). Owing to the limited availability of seeds, the RNAi-silencing construct was transformed into the cultivar "President" and the OE construct was transformed into the cultivar "Bluebyrd". Four RNAi lines and three OE lines were obtained under Kanamycin selection, propagated via tissue culture, rooted, and transplanted to soil during 2013.
Plum RNA extraction and real-time quantitative polymerase chain reaction (qPCR) To confirm silencing or OE of PpeTAC1 in transgenic plums, qPCR was performed on plum transgenic and control lines. Total RNA was extracted from apical shoots of 1-year-old, greenhouse-grown plants using the SQ Total RNA Extraction Kit (Omega Biotech, Norcross, GA) per the manufacturer's instructions. qPCR from total RNA was performed using the SYBR Green One-Step qPCR Kit (Invitrogen, Carlsbad, CA) and run in an ABI 7900HT Sequence Detection System. Each reaction was run in triplicate using 50 ng of RNA. Transgene expression was determined from 4 to 5 independent vegetatively propagated biological replicates from each transgenic line and two control plants. Expression for the controls was the average of two biological replicates. Three technical replicates were performed for each biological replicate. Quantification was performed using a standard curve derived from a serially diluted standard RNA run in parallel. A dissociation curve was run to verify that a single desired amplified product was obtained from each reaction. The PpeTAC1 primers used [For 5′-TTTG CCAAGAAACTCATCCCTCGC and Rev 5′-GCTGC TTCTGGCCATCTGATTTGT] were designed flanking an intron to amplify a 145 bp fragment from both the PpeTAC1 transgene and all native PpoTAC1 alleles.
RNAseq experiments[experimental data is being deposited into NCBI GEO and accession number will be provided] RNA extraction from the tips of peach tree branches (~4-6 cm in length) was performed as previously described 16,36 . Briefly, frozen tissue was ground in liquid N 2 using a mortar and pestle and subjected to RNA extraction using the E.Z.N.A. SQ Total RNA Kit (Omega Biotek Inc.). RNA was prepared according to the manufacturer's instructions with the exception that 2% PVP was added to the RCL buffer. For expression profiling studies,~6-10 cm of each of the three shoot tips (with leaves removed) were cut from four standard and four pillar trees from the KV991636 population (described above). Approximately 4 μg of total RNA for each tree was provided to the Cornell Weill Genomics center (Ithaca, NY) for single-read 50 bp sequencing. Raw reads were trimmed to remove reads <50 bp in length. Remaining reads were aligned to the peach genome version 2 using CLC "RNAseq" function under the following settings: maximum number of hits for a read = 10, count paired reads as two = No, expression value = RPKM, no global alignment, similarity fraction = 0.95, length fraction = 0.8, mismatch cost = 2, insertion cost = 3, deletion cost = 3 37 . The "Differential Gene Expression" function was used to identify differentially expressed genes (DEGs) in a pairwise fashion. Expression means were calculated from the four biological replicates of pillar and standard, respectively. DEGs with a P-value ≤0.05 and having a minimum average of 5 reads in both pillar and standard genotypes were used for further analysis. The DEG list was manually categorized based on published data in order to identify coordinated changes within known biological, signaling, and/or metabolic pathways 38 .

Pillar peach trees display pleotropic phenotypes
Previous studies found a number of morphological differences between pillar and standard peach growth habits in addition to branch angle including stem diameter, branch length, and internode distance 8 . In these reports, comparisons were performed between different cultivars and could not definitively account for genetic background differences that are not directly associated with the pillar trait. To directly compare between trees of the same genetic background, we measured tree height, trunk diameter, and branch length using sibling trees from a peach population segregating for the pillar trait ( Figure 1 and Table 1). Pillar peaches displayed pleotropic phenotypes and had reduced trunk diameter and stem diameters and reduced branch lengths but showed no differences in overall tree height or fruit count (Table 1 and Fig. 1b, c).

TAC1allelic variation in Prunus domestica
To better understand the role of TAC1 in regulating tree architecture, we evaluated TAC1 gene function in the closely related species P. domestica (plum). Unlike peach, plum can be readily transformed via Agrobacteriummediated strategies 39,40 . However, the use of plum is complicated by the fact that it is a hexaploid species thought to have arisen as an interspecific cross between Prunus cerasifera and Prunus spinosa 40,41 . It is also important to note that plum trees tend to have a naturally more upright architecture than peach 42,43 . To address this, we first characterized the native TAC1 alleles present in plum. The genome of the commercial plum variety  Note that apparent decreases in branch length throughout the season are due to small (<1 cm) differences in subsequent measurements of branches after they ceased growing. c Fruit set in pillar and standard trees. Note that decreases in count occur normally as a consequence of June drop "Improved French" was sequenced to an average haploid coverage of 102× via a combination of Illumina 100 bp paired-end reads and 100 bp mate-pair reads and aligned to a 5 kb fragment spanning the TAC1 gene from peach genome version 1 31 . Next, single-nucleotide polymorphisms (SNPs) and insertion/deletion (indel) variants of the native PdoTAC1 alleles were identified relative to the peach TAC1 sequence. Overall, there was a high degree of sequence identity between peach and plum. A total of 16 variants were identified within coding sequences, 5 of which were tightly clustered (Table S1). Manual inspection of this region revealed that the five clustered SNPs were errantly called due to misalignment caused by the presence of a variable simple sequence repeat (SSR) located within exon 3 that comprised multiple TAC1 alleles (Fig. 2a). Incidentally, this SSR occurred in the same location (between conserved domains II and III) as the previously identified insertion event in pillar peach (Fig. 2b) 16 . The SSR sequence consists of a repeated Asp codon (GAT 7 ). Plum alleles contained repeats ranging from GAT 6 to GAT 15 and the longest alleles included two additional 3′ Asp codons (GAC GAT). While the impact of the additional Asp residues on TAC1 protein function is unknown, the results suggest that this region may be unstable in Prunus species and could contribute to phenotypic variability.

TAC1 silenced and OX plums
Given the high level of sequence identity between the PpeTAC1 and PdoTAC1 sequences, we used a previously cloned peach TAC1 cDNA sequence to create silencing (via RNAi) and OE (under the 35S promoter) constructs for transformation into plum. Four putative RNAi transgenic lines and three putative OE lines were obtained. Owing to limited seed for the hypocotyl-based transformation method, the RNAi plums were generated using the plum cultivar President, while the OE lines were generated using seed from the cultivar Bluebyrd. These two varieties differed slightly in their growth habit (Table 2). Thus the RNAi and OE lines were only evaluated with respect to their parental cultivar background.
Using primers designed to detect all PdoTAC1 alleles, qPCR confirmed that all four RNAi lines showed significantly reduced TAC1 expression (Fig. 3a). The 3 plum OE lines showed 2-20-fold increases in TAC1 expression   (Fig. 3a). Replicate trees rooted from tissue culture explants for each line were planted to the field under APHIS permit and monitored for phenotypes. During their second year of growth, select trees were photographed in the winter and the resulting images were used to generate three-dimensional reconstructions and skeleton images. Detailed phenotypic measurements were taken from these trees during the third growing season. All four PpeTAC1 RNAi lines showed a clear pillar growth habit marked by vertically oriented secondary branches and significantly narrower branch angles when compared to untransformed trees of the same cultivar (Figs. 3b, c). In contrast to peach pillar trees, plum RNAi lines had significantly longer lateral branches and thicker trunks ( Table 2). The RNAi lines also displayed a waving or twisting of the upper branches, a phenotype that has been occasionally observed in pillar peach (Fig. 3b). When their leaves were present, the PpeTAC1 OE lines displayed few visibly obvious architectural phenotypes. However, skeletonized images revealed that their overall architecture was different (Fig. 3c). OE lines had, on average, significantly wider branch angles, thicker stems, and were overall larger trees relative to controls (Table 3; Fig. 3b, c; Figure S1).

Expression profiling of pillar trees
Currently, little is known about the mechanism of action for TAC1 in regulating branch angle. To begin to assess this, we performed an RNAseq expression profiling study using shoot tips from 4 standard and 4 pillar individuals of a segregating peach population (raw data is available for download at the NCBI Gene Expression Omnibus (GEO) accession # GSE112649). A total of 641 DEGs (P-value > 0.05) showing fold changes >1.5 were identified, with more genes being upregulated in pillar (453) than downregulated (188) (Dataset S1). To evaluate potential changes in regulatory and/or metabolic pathways, all DEGs were examined using MapMan 38 . Gene enrichment analysis revealed significant differential expression of genes involved in biotic and abiotic stress responses along with several other categories that are commonly associated with stress/defense pathways, including reactive oxygen, cell wall, and secondary   (Fig. 4, Figure S2). Manual categorization of all DEGs revealed that nearly one third of the genes in the biotic stress category 33 were comprised by NBS-LRR disease resistance genes (R genes) (Fig. 4). Among genes associated with hormone functions, the most were associated with jasmonic acid (JA) signaling including JA metabolism and terpene biosynthesis genes, which were universally upregulated in pillar trees. Genes associated with auxin and ethylene signaling were also largely upregulated while those related to brassinosteroid signaling were repressed (Fig. 4, Dataset S1). Other notable genes included MORE AXILLARY GROWTH 1 (MAX1) and MAX2, which regulate branch outgrowth and were previously shown to be differentially expressed in pillar peach trees (Dataset S1) 21 . Three DEGs encoding 2OGand Fe(II)-dependent oxygenases were identified; however, none of them appeared to be orthologous to the putative apple columnar gene MdCo31 (Dataset S1) 11 . Also of potential functional significance was a gene we recently found to be associated with the peach weeping growth habit (WEEP), which was slightly downregulated in pillar trees (Hollender et al., in press).

Discussion
Over the next half century, the Food and Agriculture Organization of the United Nations predicts that meeting the growing demand for food will require substantial increases to crop productivity without concomitant expansion of total farm acreage 44 . Crops with reduced levels of TAC1 expression could help meet this challenge having demonstrated potential for narrowing plant architectural profiles and increasing productivity via highdensity planting. Here we evaluated the effect of altered TAC1 expression in both peach and plum trees.
Segregating pillar peach siblings carrying a knockout mutation of PpeTAC1 had consistently shorter branches and thinner trunks and lateral branches while tree height was unchanged. Fruit yield showed no significant differences between pillar and standard trees. These data confirm a previous report that the pillar cultivar "Crimson Rocket" produced shorter and fewer sylleptic branches compared to a standard peach cultivar and that differences in branch growth rate occurred primarily in spring 8 .
We previously discovered that "Italian" pillar peach trees have a repetitive element inserted within a repetitive region of PpeTAC1 Exon 3 16 . Surprisingly, Illumina DNA sequencing of the plum (P. domestica) cultivar "Improved French" revealed the presence of a variable length GAT repeat in this same location resulting in six distinct alleles within the hexaploid genome. These changes resulted in some alleles encoding long runs of aspartic acid residues 13-17 amino acids in length. A single shorter allele (GAT 6 ) was similar in length to that found in peach (GAT 7 ). Coincidentally, many plum cultivars including "Improved French", "President", and "Bluebyrd" plum (controls depicted in Fig. 3b), have a more upright stature than commonly observed in other Prunus species 42,43 . We speculate that the additional Asp residues in P. domestica alleles may result in TAC1 proteins that have reduced function or are less stable leading to a more upright tree shape. In animals, these genetic "stutters" within coding sequences have been shown to underlie genetic disorders such as Huntington's disease where CAG repeat expansion produces longer tracts of glutamines within the huntingtin protein leading to protein instability 45 .
Silencing of PdoTAC1 via RNAi resulted in plum trees having architectural profiles that were significantly narrower than those observed in peach. Branch growth in these trees was extremely vertical producing slender but dense canopies. Mean branch angles across all RNAi plum lines were 17 o versus pillar peaches, which have been reported to be as wide as 40 o46 . The RNAi plum lines also exhibited additional differences from the pillar peach including longer branch lengths and thicker trunk diameters. In contrast, transgenic plums overexpressing PpeTAC1 displayed wider branch angles, thicker branches, and increased overall tree height relative to nontransgenic controls. Given the lack of knowledge about the cellular function of TAC1, it is difficult to conceive of a model by which these diverse phenotypic effects could be explained aside from generalities regarding potential changes in levels of hormones or their translocation, perception, and/or signaling. A prior study attributed some of the phenotypic variation within pillar peaches to ET (6) BR (3) CK GA (5) SL (5) ABA (6) IAA (10) JA (19) Hormone PR genes (11) R genes (33) All other (74) BioƟc/ abioƟc stress PlasƟd (21) ROS (27) Secondary metabolism (20) Cell wall (26) General Categories higher auxin levels and higher auxin-cytokinin ratios found throughout the canopy 8 . The contrasting results in plum suggest that TAC1 RNAi plums may have hormone profiles that are either different from peach or that the phenotypic effects result from other unknown genetic or physiological variables. Collectively, these data show that biotechnology strategies can be used to engineer tree shapes by altering TAC1 levels; however, these may be accompanied by additional phenotypic effects that are less predictable.
Expression profiling data support a hypothesis that loss of TAC1 leads to changes in key plant hormones as responses to auxin and jasmonic acid were consistently upregulated. Surprisingly, pillar trees exhibited pronounced changes in biotic/abiotic stress signaling including numerous R genes and pathogenesis-related genes. This effect may be caused by or associated with the upregulation of JA metabolism and/or signaling, which is a known defense hormone. Categories related to defense/ stress responses were also differentially expressed including reactive oxygen, secondary metabolism, and cell wall biosynthesis. Overall, these expression changes are similar to what was previously reported for apple columnar trees, which exhibited changes in genes associated with cell wall metabolism, defense, and JA responses 47,48 . At this time, we cannot rule out the possibility that the TAC1 locus is tightly linked to genetic variation within an unrelated gene that affects biotic stress signaling pathways. While the data support the possibility that TAC1 functions via altered defense hormone profiles and/or translocation, additional in-depth molecular studies will be needed to better understand these relationships.