# Early transcriptional response to gravistimulation in poplar without phototropic confounding factors

David Lopez, Jérôme Franchel, Jean-Stéphane Venisse, Joël J. Drevet, Philippe Label, Catherine Coutand, Patricia Roeckel-Drevet

► To cite this version:

David Lopez, Jérôme Franchel, Jean-Stéphane Venisse, Joël J. Drevet, Philippe Label, et al.. Early transcriptional response to gravistimulation in poplar without phototropic confounding factors. AoB Plants, 2021, 13 (1), 10.1093/aobpla/plaa071. hal-03149688

HAL Id: hal-03149688

https://hal.inrae.fr/hal-03149688

Submitted on 23 Feb 2021

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
OXFORD

<figure><img src="image_2.png" /></figure>

AoBP
AoB PLANTS

AoB PLANTS 2021, Vol. 13, No. 1

doi:10.1093/aobpla/plaa071

Advance Access Publication December 31, 2020
Studies

STUDIES

# Early transcriptional response to gravistimulation in poplar without phototropic confounding factors

David Lopez¹,², Jérôme Franchel³, Jean-Stéphane Venisse³, Joël R. Drevet⁴, Philippe Label³, Catherine Coutand⁵ and Patricia Roeckel-Drevet³*

¹CIRAD, UMR AGAP, F-34398 Montpellier, France, ²AGAP, Univ Montpellier, CIRAD, INRAE, Institut Agro, Montpellier, France, ³Université Clermont Auvergne, INRAE, PIAF, Campus Universitaire des Cézeaux, 1 Impasse Amélie Murat, TSA 60026, 63178 Aubière Cedex, France, ⁴Université Clermont Auvergne, GReD INSERM U1103-CNRS UMR 6293, Faculté de Médecine, CRBC (Centre de Recherche Bio-Clinique), 28 Place Henri Dunant, 63000 Clermont-Ferrand, France, ⁵INRAE, UR 115 PSH, Centre de recherche PACA, 228, route de l'aérodrome, CS 40509, 84914 Avignon Cedex 9, France

*Corresponding author's e-mail address: patricia.drevet@uca.fr

Form & Function. Chief Editor: Kate McCulloh

Associate Editor: Brian Atwell

## Abstract

In response to gravistimulation under anisotropic light, tree stems showing an active cambium produce reaction wood that redirects the axis of the trees. Several studies have described transcriptomic or proteomic models of reaction wood relative to the opposite wood. However, the mechanisms leading to the formation of reaction wood are difficult to decipher because so many environmental factors can induce various signalling pathways leading to this developmental reprogramming. Using an innovative isotropic device where the phototropic response does not interfere with gravistimulation we characterized the early molecular responses occurring in the stem of poplar after gravistimulation in an isotropic environment, and without deformation of the stem. After 30 min tilting at 35° under anisotropic light, we collected the upper and lower xylems from the inclined stems. Controls were collected from vertical stems. We used a microarray approach to identify differentially expressed transcripts. High-throughput real-time PCR allowed a kinetic experiment at 0, 30, 120 and 180 min after tilting at 35°, with candidate genes. We identified 668 differentially expressed transcripts, from which we selected 153 candidates for additional Fluidigm qPCR assessment. Five candidate co-expression gene clusters have been identified after the kinetic monitoring of the expression of candidate genes. Gene ontology analyses indicate that molecular reprogramming of processes such as 'wood cell expansion', 'cell wall reorganization' and 'programmed cell death' occur as early as 30 min after gravistimulation. Of note is that the change in the expression of different genes involves a fine regulation of gibberellin and brassinosteroid pathways as well as flavonoid and phosphoinositide pathways. Our experimental set-up allowed the identification of genes regulated in early gravitropic response without the bias introduced by phototropic and stem bending responses.

**Keywords:** Cell wall; Fluidigm; gravitropism; gravity; isotropic light; phototropism; Populus alba; Populus tremula; qPCR; tension wood; transcriptome; xylem poplar.

## Introduction

The perception of Earth's gravity leads to a gravity-driven growth process called gravitropism. In higher plants, shoots grow upwards (negative gravitropism), while roots grow downwards

(positive gravitropism). Gravitropism includes the perception of gravity, signal transmission and growth response. Many studies were conducted on roots in which the key gravity detection site

Received: 2 July 2020; Editorial decision: 1 December 2020; Accepted: 30 December 2020

© The Author(s) 2020. Published by Oxford University Press on behalf of the Annals of Botany Company.
This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted reuse, distribution, and reproduction in any medium, provided the original work is properly cited.
(columella cells at the root tip) is separated from the bending zone (extension zone) (for a review of the topic, see Baldwin *et al.* 2013). Several plant species were also used to study the response to gravistimulation of aerial parts including maize coleoptiles, pulvini barley, Pisum hypocotyls, Arabidopsis inflorescence stems and hypocotyls (for a review, see Hashiguchi *et al.* 2013). In these herbaceous species, vertical root or shoot movements are caused by differential elongation between the upper and lower parts of the gravistimulated organs. In trees, two different motors allow the curvature of the stem as reorientation is obtained by differential growth in areas of primary elongation, while it is obtained by differential cellular differentiation between the two sides of the stem in areas with active cambium (for a review, see Groover 2016). Consequently, the cambium produces an asymmetric 'reaction wood' (Côté and Day 1964). In gymnosperms, 'reaction wood' is called 'compression wood' and develops on the underside of the inclined rods where it generates a compressive force to push the stem up (Timell 1986). In arboreal dicotyledonous angiosperms, the 'reaction wood' is called 'tension wood' (TW) and forms on the upper face of the inclined stem when the tree begins the reorientation process (for a review, see Du and Yamamamoto 2007). The 'opposite wood' forms on the underside and is anatomically similar to the 'normal wood' formed by straight stems. Tension wood is a particular type of secondary xylem that generates high internal tension forces, mainly due to the longitudinal shrinkage of maturing wood cells (Kubler 1959). Ruelle (2014) reviewed the variations in the anatomy and ultrastructure of TW compared to the 'opposite wood' (opposite side of the stem) and 'normal wood' (in a non-tilted tree). In many angiosperms, TW is characterized by the presence of gelatinous fibres (G-layer) of almost pure cellulose.

The molecular basis of the early events that followed gravistimulation and led to the formation of TW is still fragmentary. The literature review showed that most studies using comprehensive approaches such as proteomics and transcriptomics have examined the formation of 'reaction wood' at developmental stages when it is already histologically observable and/or when the tree undergoes recovery movement (for a review, see Tocquard et al. 2014). For example, Andersson-Gunneras et al. (2006) pioneering work reported a complete quantification of gene expression and related metabolites occurring during TW formation in poplar. In this study, xylem samples from plants that had been tilted for 3 weeks were used, focussing the analysis on the stage of the secondary wall formation. Chen et al. (2015) characterized gene expression patterns, in TW, opposite and normal wood in 30-year-old Populus tomentosa. Gerttula et al. (2015) produced genomic transcriptomes for 3-month-old Populus tremula × alba, stimulated either by gravistimulation (placed horizontally for 2 days), gibberellic acid treatment or ARBORKNOX2 expression.

In the present work, we have focussed on deciphering transcriptomic events caused by stem tilting before the beginning of the righting movement. For that reason, xylem samples were collected after 30 min of tilting, before the appearance of gravitational motion and TW formation. In addition, to our knowledge, previous experimental designs did not distinguish gravisensing from mechanosensing signalling pathways (Lopez et al. 2014; Groover 2016). For example, although Azri et al. (2009, 2013) have characterized the proteome of young poplars inclined at 35° for 45 min. In that particular set-up, the inclined organ induces flexion under its own weight, which can be considered as a thigmomorphogenetic stimulus (Coutand 2010). To avoid such situation, we limited our analyses to the part of

the stem showing secondary growth. The trees were planted a week before use, so there was no bending due to the weight of the tree after tipping. Finally, and most importantly, phototropic and gravitropic responses as well as the autotropic response, which is essentially the straightening of a curved organ, are processes in constant interaction under natural conditions that were not dissociated in the experimental situations reported so far (Bastien et al. 2015).

In a recent study, we succeeded in dissociating these complex stimuli thanks to an original device that allowed us to combine isotropic lighting and tilting of young poplars (Coutand et al. 2019), thus avoiding the phototropic response. In such a device, we showed that the gravitropic response was partly distinct when comparing an anisotropic light condition to an isotropic light condition (Coutand et al. 2019). To better understand the molecular basis of these responses, we used here this device that allows us to distinguish between gravitropic, mechanosensitive and phototropic responses in plants in order to study the molecular adjustments due to gravistimulation on stems in the early phases of 'reaction wood' formation.

# Materials and Methods

## Plant material

The study was conducted on young poplars (Populus tremula × Populus alba cv 717 1B4) produced by in vitro propagation (Azri et al. 2009). They were placed in growth chambers after acclimation (photoperiod: 16 h/8 h; day: 24 °C, night: 18 °C). For transcriptome studies, plants with straight stem, aged 3 months, with a height of between 40 and 50 cm and an average diameter at half height of 10 mm were carefully and loosely staked avoiding any mechanical stress 1 week before use. The side of the stem to the stake was marked with dots using a black marker.

## Isotropic light device and tilting of the plants

The spheres used for gravistimulation were 150 cm wide and uniformly bright (i.e. isotropic light; Coutand et al. 2019) (Fig. 1). Growing conditions were 16 h/8 h photoperiod; 24 °C (day)/18 °C (night). The temperature was maintained by a controlled cooling system (blyss, WAP 267EC).

One week after staking (staking took place in the growth chamber), the plants reached 55 cm (SD = 3.7 cm, n = 76). Staked trees were placed in the sphere (one plant per sphere). First of all, the plants were left in an upright position for 24 h. Then, 7 h after the beginning of the photoperiod, the plants were tilted 35° from the vertical. The control plants were kept straight in the sphere. When tilted, the stake was oriented beside the upper face of the stem to prevent the stem from bending due to its own weight (Fig. 2A).

## Tissue sampling

After treatment, for each plant, we selected a stem segment (15 cm long, 5 cm above ground level, i.e. internodes 4–11) that did not contain any TW according to Astra-blue/Safranine staining (Azri et al. 2009). In this part of the stem, the cambium is active and a secondary xylem is formed. Stem segments 0.5 mm long were taken from internodes 4, 7 and 11 for further microscopic checks. The rest of the stem was divided longitudinally into two parts using a scalpel (Fig. 2B). The lower (opposite the stake) and upper (towards the stake) sides of the stems were collected; the bark was peeled and discarded (Fig. 2C). Lower xylem (XL) and upper xylem (XU) were instantly frozen in liquid nitrogen and
<figure><img src="image_3.png" /><figcaption>**Figure 1.** Experimental device of isotropic light to disentangle gravitropic and phototropic plant responses. The isotropic light is obtained from circular neon tubes arranged in a turtle pattern at the periphery of the white plastic hemisphere. The light is concentrated towards the inside of the hemisphere by metallic plates covering the device. The poplar plant is placed in the centre of the device, and the hemispheres are tightly connected (Coutand et al. 2019 for the full description of the isotropic device). (A) Straight young poplar plant. (B) Plant subjected to tilting.</figcaption></figure>

stored. We used the same protocol for straight stems (controls). The xylem to the stake was used for the comparison with XU, while the xylem opposite the stake was used for the comparison with XL.

## RNA extraction

When no TW was detected in the stem by microscopic controls, the xylem samples were crushed in liquid nitrogen. The total RNA was extracted from XU and XL samples and controls (also arbitrarily referred to as XU and XL) using the Qiagen RNaeasy Plant Kit according to the supplier's instructions (Qiagen, Courtaboeuf, France).

## Transcriptome studies using DNA microarray technology

The discovery of genes whose expression has been modulated by inclination has been carried out using a microarray approach. RNA samples were obtained from the upper or lower xylem of four inclined plants (30 min, 35°) and four straight plants. For each comparison, a technical repeat with fluorochrome inversion including four biological repetitions was performed (i.e. 4 dye-switch hybridizations per comparison).

The hybridization, scanning and analysis of the biochips were performed by the Plant Genomics Research Unit (Evry, France) according to their standard procedures. From 200 ng total RNA, amplification and labelling of cDNA with Cy3-dUTP or Cy5-dUTP (Perkin-Elmer-NEN Life Science Products) was performed with the TransPlex Complete Whole Transcriptome Amplification Kit (WT2A; Sigma-Aldrich, Inc.) as recommended by the manufacturer. The Populus microarray based on Roche-NimbleGen technology (Roche, Boulogne-Billancourt, France) was used. A single high-density Populus microarray slide contains 12 chambers, each containing 135 029 primers representing 49 474 gene models based on the annotation of the Populus trichocarpa v1.1 genome. Hybridization and washing were carried out in accordance with the recommendations of NimbleGen (Roche NimbleGen Technologies, Inc.). A two-micron scan was performed with the InnoScan900 scanner (InnopsysR, Carbonne, France) and the raw data were extracted with the MapixR software (InnopsysR, Carbonne, France). The measured fold variation was corrected according to Bonferroni, followed by a Benjamini-Hochberg P-value adjustment (Benjamini and Hochberg 2000). Biochip data from this study have been deposited at the Gene Expression Omnibus (http://www.ncbi.nlm.nih.gov/geo/, access number GSE86951) and CATdb (http://urgv.evry.inra.fr/CATdb/; Project: 12plex-poplar-2012_01) according to MIAME standards.

## Statistical analysis of microarray data

For each array, the raw data included the logarithm of the median pixel intensity of the characteristics at 635 nm (red) and 532 nm (green) wavelengths. For each array, global normalization according to intensity was performed using the LOESS procedure (Yang et al. 2002) to correct for dye bias. The differential analysis was based on the average of log-ratios on duplicate probes and technical repetitions. Therefore, the number of data available for each gene is equal to the number of biological replicates and was used to calculate the moderate t-test (Smyth 2004). Under the null hypothesis, Limma (version 3.34.8) found no evidence that specific variances varied from one probe to another and, therefore, the moderate t-statistics was assumed to follow a standard normal distribution. To control the rate of false discoveries, the corrected P-values found were calculated using the optimized FDR approach (Storey 2003). We considered the probes with an adjusted P-value ≤ 0.05 to be differentially expressed.

The analysis was done with the R software (version 3.4.3, R Core Team (2017). R: A language and environment for statistical calculation. R Foundation for Statistical Computing, Vienna, Austria. https://www.R-project.org/). The SqueezeVar function of the Limma library was used to smooth specific variances by calculating empirical Bayesian posterior averages. The library kerfdr (version 2.0.1) was used to calculate the adjusted P-values.

Data management, statistical analysis and hierarchical grouping were performed using R. Hierarchical grouping was calculated using the hclust function (set of statistics) using Euclidean distances and the Ward method. The false discovery rate was calculated using the p adjust function (set of statistics) using the Benjamini and Yekutieli (2001) method.

## Validation of microarray results

For the qPCR validation, we selected 153 genes amongst the most differentially expressed genes from biochip analyses. Before proceeding with the qPCR analysis, the cDNA was synthesized from 0.5 µg total RNA using SuperScript III (Invitrogen, USA) with
<figure><img src="image_4.png" /><figcaption>Figure 2. Scheme of the procedures. (A) Isotropic device containing a stalked plant. (B) Longitudinal splitting of the stem. (C) Isolation of XU samples used for RNA extraction and microscopical controls. (D) Microarray experiment. (E) Candidate gene exploration.</figcaption></figure>

a 1:1 mixture of oligo-dT and random hexamers. The primer pairs were designed using the Quantprime online service (www.quantprime.de) [see Supporting Information—Table S1]. The specificity of the primers was ensured by checking the melting curve and sequencing the PCR products on the cDNA. The primers were validated when the efficiency was close to 100 % on the slopes of the standard curve. Relative quantification with respect to vertical controls was performed using the 2⁻ΔΔCq method developed by Pfaffl (2001). The expression levels of the 153 genes obtained by qPCR and microarray analysis were compared using an analysis of variance followed by a Tukey Honest Significance Difference (HSD) post hoc test.

## Kinetics study of candidate gene expression

The plants were either straight or inclined at 35° for 30, 60, 120 or 180 min in the isotropic device. Eight independent biological replicates (individual plants) were produced for each condition. XU and XL cDNA samples were used to monitor the expression of the 153 candidate genes by real-time qPCR at high-throughput. The experiments were performed according to the Fluidigm gene expression protocol without modification using the Biomark platform (Fluidigm, USA) and the 96x96 IFC gene expression (integrated fluidic circuit). Relative quantification was performed using the 2⁻ΔΔCq method developed by Pfaffl (2001). Expression levels were related to the geometric mean of the 153 expression levels of the candidate genes proposed by the Bestkeeper reference research software (Vandesompele

et al. 2002; Pfaffl et al. 2004), but using a larger subset of genes. No significant Tukey post hoc HSD was found in the change in gene expression (P < 0.05, data not shown) between qPCR and microarray, which validated the microarray results.

## Bioinformatics analyses

The functional annotation was retrieved using a standalone version of the BioMart software (http://biomart.org) and the genomic annotation P. trichocarpa v3.0 (http://phytozome.jgi.doe.gov). The singular AgriGO enrichment analysis (http://bioinfo.cau.edu.cn/agriGO/, Zhou et al. 2010) was used for the analysis of gene set enrichment (GSEA) using the genetic ontology terms (GO) annotated for genes identified on the DNA chip and the default AgriGO parameters, adapted to large data sets (Fisher test, Yekutieli multiple test adjustment method and significance level at 0.05). We evaluated the possible bias of the GO term enrichment procedure when using the Nimblegene poplar network by comparing the whole gene models, based on P. trichocarpa V1.1, present on the Nimblegene array with the P. trichocarpa V3 genome gene models. No significant representativeness bias was detected. Thus, we considered that the enrichment of the term GO based on Nimblegene data was valid when using GO annotations of the poplar genome V3. Further analyses were carried out using PopGenIE (http://popgenie.org, Sjödin et al. 2009; Sundell et al. 2015) for the analysis of qPCR clusters and the discovery of miRNA targets.
# Results

## Neutral discovery of early modulated genes

In the isotropic device, the phototropic response cannot occur. This device was used to identify genes regulated by gravistimulation in forming 'reaction wood' using DNA chip technology (Nimblegene microarray). Analyses of xylem samples were carried out by comparing the expression levels of vertical (control) and inclined plants in the lower part of the xylem (XL, i.e. the face facing the ground when it is inclined) and in the upper part of the xylem (XU). Array hybridizations resulted in a subset of 2355 gene models (poplar genome version 1.1) expressed differentially after Benjamini and Hochberg as well as Bonferroni multiple test corrections were made (P-value < 0.05; see Supporting Information—Table S2). Our experimental set-up provided well-controlled conditions, as evidenced by the reproducibility of the microarray experiments [see Supporting Information—Fig. S1].

From the list of 2355 genes, 668 showed significant differential expression between the two conditions (annotation was performed using the poplar genome version 3; see Supporting

Information—Table S3). Of these, 214 and 250 genes were modulated exclusively in XU or XL, respectively, while 204 were modulated on both sides of the xylem in response to gravity (Fig. 3A). For this last set of genes, there is no tendency for 'side-specific' expression, as confirmed by the hierarchical grouping of biological repetitions showing that the XL and XU samples are mixed. A high-resolution version of the heatmaps containing the names of the gene models is available in Supporting Information—Fig. S2. As shown in the expression maps and density histogram (Figs 3A and B), the most differentially expressed genes are within a range of 2-fold variation between inclined and control trees. Figure 3C shows, a total of 379 genes were always regulated upward (133 + 121 + 125) and 262 genes (81 + 56 + 125) were always regulated downward.

## Functional exploration

In order to clarify the role of modulated genes in response to gravistimulation, the GSEA was performed using GO term enrichment with the three different ontologies available (Biological Process (BP); Molecular Function (MF); Cell Component (CC)) and calculated with AgriGO. Considering the

<figure><img src="image_5.png" /><figcaption>**Figure 3.** Identification of poplar genes with expression significantly regulated in response to a 30min, 35° gravi-stimulation in xylem tissues. (A) Heatmaps of differentially expressed genes found in Upper Xylem (XU, 214), Lower Xylem (XL) and in both XU and XL. Gene models are hierarchically clustered on the Y-axis and the four biological repeats are hierarchically clustered on the X-axis. (B) Key of the heatmaps, fold changes are represented on the X-axis and frequencies on the Y-axis. (C) Detail by tissue and trend of the gene expression patterns. Note that 27 genes showed a significant regulation in both XU and XL but with opposite trends. However, none of these "opposite" trends were statistically significant hence, they are not accounted in the "XU and XL" column.</figcaption></figure>
<table><caption>Table 1. GO term enrichment analysis of the 668 differentially expressed xylem genes (35° tilted vs. control).</caption><thead><tr><th>GO term</th><th>Ontology</th><th>Description</th><th>Number in input list</th><th>Number in BG/Ref</th><th>P-value</th><th>FDR</th><th>GO term</th><th>Ontology</th><th>Description</th><th>Number in input list</th><th>Number in BG/Ref</th><th>P-value</th><th>FDR</th></tr></thead><tbody><tr><td>GO:0005975</td><td>P</td><td>Carbohydrate metabolic process</td><td>71</td><td>2018</td><td>5.10E-10</td><td>1.50E-06</td><td>GO:0046271</td><td>P</td><td>Phenylpropanoid catabolic process</td><td>5</td><td>37</td><td>0.00046</td><td>0.025</td></tr><tr><td>GO:0034641</td><td>P</td><td>Cellular nitrogen compound metabolic process</td><td>45</td><td>1031</td><td>2.40E-09</td><td>2.30E-06</td><td>GO:0042180</td><td>P</td><td>Cellular ketone metabolic process</td><td>49</td><td>1886</td><td>0.00056</td><td>0.03</td></tr><tr><td>GO:0005976</td><td>P</td><td>Polysaccharide metabolic process</td><td>34</td><td>639</td><td>2.40E-09</td><td>2.30E-06</td><td>GO:0070297</td><td>P</td><td>Regulation of two-component signal transduction system (phosphorelay)</td><td>5</td><td>40</td><td>0.00064</td><td>0.032</td></tr><tr><td>GO:0044264</td><td>P</td><td>Cellular polysaccharide metabolic process</td><td>26</td><td>457</td><td>5.30E-08</td><td>3.80E-05</td><td>GO:0010104</td><td>P</td><td>Regulation of ethylene-mediated signalling pathway</td><td>5</td><td>40</td><td>0.00064</td><td>0.032</td></tr><tr><td>GO:0016051</td><td>P</td><td>Carbohydrate biosynthetic process</td><td>31</td><td>641</td><td>9.40E-08</td><td>5.40E-05</td><td>GO:0010105</td><td>P</td><td>Negative regulation of ethylene-mediated signalling pathway</td><td>5</td><td>40</td><td>0.00064</td><td>0.032</td></tr><tr><td>GO:0006519</td><td>P</td><td>Cellular amino acid and derivative metabolic process</td><td>57</td><td>1688</td><td>1.10E-07</td><td>5.40E-05</td><td>GO:0070298</td><td>P</td><td>Negative regulation of two-component signal transduction system (phosphorelay)</td><td>5</td><td>40</td><td>0.00064</td><td>0.032</td></tr><tr><td>GO:0006073</td><td>P</td><td>Cellular glucan metabolic process</td><td>20</td><td>298</td><td>1.60E-07</td><td>6.50E-05</td><td>GO:0009063</td><td>P</td><td>Cellular amino acid catabolic process</td><td>9</td><td>144</td><td>0.00067</td><td>0.033</td></tr><tr><td>GO:0044262</td><td>P</td><td>Cellular carbohydrate metabolic process</td><td>41</td><td>1093</td><td>6.10E-07</td><td>0.00022</td><td>GO:0009266</td><td>P</td><td>Response to temperature stimulus</td><td>35</td><td>1229</td><td>0.00077</td><td>0.037</td></tr><tr><td>GO:0044042</td><td>P</td><td>Glucan metabolic process</td><td>20</td><td>329</td><td>7.00E-07</td><td>0.00022</td><td>GO:0009310</td><td>P</td><td>Amine catabolic process</td><td>9</td><td>149</td><td>0.00084</td><td>0.04</td></tr><tr><td>GO:0019748</td><td>P</td><td>Secondary metabolic process</td><td>44</td><td>1235</td><td>8.90E-07</td><td>0.00026</td><td>GO:0022611</td><td>P</td><td>Dormancy process</td><td>5</td><td>43</td><td>0.00087</td><td>0.04</td></tr><tr><td>GO:0009066</td><td>P</td><td>Aspartate family amino acid metabolic process</td><td>13</td><td>145</td><td>1.20E-06</td><td>0.00032</td><td>GO:0046395</td><td>P</td><td>Carboxylic acid catabolic process</td><td>12</td><td>253</td><td>0.00097</td><td>0.043</td></tr><tr><td>GO:0034637</td><td>P</td><td>Cellular carbohydrate biosynthetic process</td><td>24</td><td>483</td><td>1.70E-06</td><td>0.00037</td><td>GO:0000097</td><td>P</td><td>Sulfur amino acid biosynthetic process</td><td>6</td><td>67</td><td>0.00095</td><td>0.043</td></tr><tr><td>GO:0007018</td><td>P</td><td>Microtubule-based movement</td><td>14</td><td>174</td><td>1.60E-06</td><td>0.00037</td><td>GO:0016054</td><td>P</td><td>Organic acid catabolic process</td><td>12</td><td>253</td><td>0.00097</td><td>0.043</td></tr><tr><td>GO:0009308</td><td>P</td><td>Amine metabolic process</td><td>41</td><td>1171</td><td>3.20E-06</td><td>0.00061</td><td>GO:0009873</td><td>P</td><td>Ethylene-mediated signalling pathway</td><td>10</td><td>185</td><td>0.00099</td><td>0.043</td></tr><tr><td>GO:0006555</td><td>P</td><td>Methionine metabolic process</td><td>10</td><td>88</td><td>3.00E-06</td><td>0.00061</td><td>GO:0009312</td><td>P</td><td>Oligosaccharide biosynthetic process</td><td>7</td><td>95</td><td>0.0011</td><td>0.046</td></tr><tr><td>GO:0033692</td><td>P</td><td>Cellular polysaccharide biosynthetic process</td><td>18</td><td>311</td><td>4.70E-06</td><td>0.00085</td><td>GO:0016491</td><td>F</td><td>Oxidoreductase activity</td><td>97</td><td>3249</td><td>1.30E-09</td><td>1.30E-06</td></tr></tbody></table>
<table><caption>Table 1. Continued</caption><thead><tr><th>GO term</th><th>Ontology</th><th>Description</th><th>Number in input list</th><th>Number in BG/Ref</th><th>P-value</th><th>FDR</th><th>GO term</th><th>Ontology</th><th>Description</th><th>Number in input list</th><th>Number in BG/Ref</th><th>P-value</th><th>FDR</th></tr></thead><tbody><tr><td>GO:0005978</td><td>P</td><td>Glycogen biosynthetic process</td><td>6</td><td>24</td><td>5.80E-06</td><td>0.00098</td><td>GO:0005506</td><td>F</td><td>Iron ion binding</td><td>49</td><td>1268</td><td>1.90E-08</td><td>9.30E-06</td></tr><tr><td>GO:0000271</td><td>P</td><td>Polysaccharide biosynthetic process</td><td>19</td><td>354</td><td>7.20E-06</td><td>0.0012</td><td>GO:0046906</td><td>F</td><td>Tetrapyrrole binding</td><td>34</td><td>819</td><td>6.80E-07</td><td>0.00022</td></tr><tr><td>GO:0044106</td><td>P</td><td>Cellular amine metabolic process</td><td>33</td><td>890</td><td>9.70E-06</td><td>0.0015</td><td>GO:0020037</td><td>F</td><td>Heme binding</td><td>31</td><td>744</td><td>2.00E-06</td><td>0.00047</td></tr><tr><td>GO:0006575</td><td>P</td><td>Cellular amino acid derivative metabolic process</td><td>33</td><td>908</td><td>1.40E-05</td><td>0.0021</td><td>GO:0016638</td><td>F</td><td>Oxidoreductase activity, acting on the CH-NH<sub>2</sub> group of donors</td><td>10</td><td>94</td><td>5.20E-06</td><td>0.001</td></tr><tr><td>GO:0042221</td><td>P</td><td>Response to chemical stimulus</td><td>114</td><td>4909</td><td>1.50E-05</td><td>0.0021</td><td>GO:0009055</td><td>F</td><td>Electron carrier activity</td><td>43</td><td>1346</td><td>1.60E-05</td><td>0.0026</td></tr><tr><td>GO:0042546</td><td>P</td><td>Cell wall biogenesis</td><td>17</td><td>313</td><td>1.90E-05</td><td>0.0025</td><td>GO:0004091</td><td>F</td><td>Carboxylesterase activity</td><td>24</td><td>564</td><td>2.00E-05</td><td>0.0028</td></tr><tr><td>GO:0009409</td><td>P</td><td>Response to cold</td><td>31</td><td>855</td><td>2.70E-05</td><td>0.0034</td><td>GO:0004553</td><td>F</td><td>Hydrolase activity, hydrolyzing O-glycosyl compounds</td><td>29</td><td>773</td><td>2.70E-05</td><td>0.0033</td></tr><tr><td>GO:0009415</td><td>P</td><td>Response to water</td><td>24</td><td>584</td><td>3.40E-05</td><td>0.004</td><td>GO:0016758</td><td>F</td><td>Transferase activity, transferring hexosyl groups</td><td>27</td><td>707</td><td>3.80E-05</td><td>0.0037</td></tr><tr><td>GO:0005977</td><td>P</td><td>Glycogen metabolic process</td><td>7</td><td>52</td><td>3.50E-05</td><td>0.004</td><td>GO:0046914</td><td>F</td><td>Transition metal ion binding</td><td>101</td><td>4311</td><td>3.80E-05</td><td>0.0037</td></tr><tr><td>GO:0006112</td><td>P</td><td>Energy reserve metabolic process</td><td>7</td><td>54</td><td>4.40E-05</td><td>0.0046</td><td>GO:0004497</td><td>F</td><td>Monooxygenase activity</td><td>25</td><td>638</td><td>5.00E-05</td><td>0.004</td></tr><tr><td>GO:0000096</td><td>P</td><td>Sulfur amino acid metabolic process</td><td>10</td><td>122</td><td>4.20E-05</td><td>0.0046</td><td>GO:0019825</td><td>F</td><td>Oxygen binding</td><td>19</td><td>409</td><td>4.80E-05</td><td>0.004</td></tr><tr><td>GO:0009067</td><td>P</td><td>Aspartate family amino acid biosynthetic process</td><td>8</td><td>75</td><td>4.50E-05</td><td>0.0046</td><td>GO:0008509</td><td>F</td><td>Anion transmembrane transporter activity</td><td>14</td><td>251</td><td>7.70E-05</td><td>0.0057</td></tr><tr><td>GO:0009698</td><td>P</td><td>Phenylpropanoid metabolic process</td><td>25</td><td>639</td><td>5.10E-05</td><td>0.005</td><td>GO:0016762</td><td>F</td><td>Xyloglucan:xyloglucosyl transferase activity</td><td>6</td><td>44</td><td>0.00012</td><td>0.0082</td></tr><tr><td>GO:0006520</td><td>P</td><td>Cellular amino acid metabolic process</td><td>31</td><td>890</td><td>5.60E-05</td><td>0.0053</td><td>GO:0016757</td><td>F</td><td>Transferase activity, transferring glycosyl groups</td><td>31</td><td>954</td><td>0.00019</td><td>0.011</td></tr><tr><td>GO:0009250</td><td>P</td><td>Glucan biosynthetic process</td><td>12</td><td>184</td><td>6.20E-05</td><td>0.0057</td><td>GO:0016798</td><td>F</td><td>Hydrolase activity, acting on glycosyl bonds</td><td>29</td><td>864</td><td>0.00018</td><td>0.011</td></tr><tr><td>GO:0010817</td><td>P</td><td>Regulation of hormone levels</td><td>18</td><td>387</td><td>7.40E-05</td><td>0.0065</td><td>GO:0003824</td><td>F</td><td>Catalytic activity</td><td>344</td><td>19164</td><td>0.00023</td><td>0.013</td></tr><tr><td>GO:0019439</td><td>P</td><td>Aromatic compound catabolic process</td><td>7</td><td>59</td><td>7.30E-05</td><td>0.0065</td><td>GO:0005200</td><td>F</td><td>Structural constituent of cytoskeleton</td><td>7</td><td>73</td><td>0.00025</td><td>0.013</td></tr><tr><td>GO:0009832</td><td>P</td><td>Plant-type cell wall biogenesis</td><td>14</td><td>252</td><td>8.00E-05</td><td>0.0068</td><td>GO:0015149</td><td>F</td><td>Hexose transmembrane transporter activity</td><td>6</td><td>59</td><td>0.00051</td><td>0.026</td></tr></tbody></table>
<table><caption>Table 1. Continued</caption><thead><tr><th>GO term</th><th>Ontology</th><th>Description</th><th>Number in input list</th><th>Number in BG/Ref</th><th>P-value</th><th>FDR</th><th>GO term</th><th>Ontology</th><th>Description</th><th>Number in input list</th><th>Number in BG/Ref</th><th>P-value</th><th>FDR</th></tr></thead><tbody><tr><td>GO:0009719</td><td>P</td><td>Response to endogenous stimulus</td><td>69</td><td>2734</td><td>9.60E-05</td><td>0.0079</td><td>GO:0022891</td><td>F</td><td>Substrate-specific transmembrane transporter activity</td><td>48</td><td>1842</td><td>0.0006</td><td>0.027</td></tr><tr><td>GO:0009834</td><td>P</td><td>Secondary cell wall biogenesis</td><td>9</td><td>112</td><td>0.00012</td><td>0.0093</td><td>GO:0008471</td><td>F</td><td>Laccase activity</td><td>5</td><td>39</td><td>0.00058</td><td>0.027</td></tr><tr><td>GO:0009414</td><td>P</td><td>Response to water deprivation</td><td>22</td><td>560</td><td>0.00013</td><td>0.01</td><td>GO:0016641</td><td>F</td><td>Oxidoreductase activity, acting on the CH-NH<sub>2</sub>, group of donors, oxygen as acceptor</td><td>6</td><td>63</td><td>0.00071</td><td>0.031</td></tr><tr><td>GO:0042445</td><td>P</td><td>Hormone metabolic process</td><td>13</td><td>232</td><td>0.00013</td><td>0.01</td><td>GO:0043169</td><td>F</td><td>Cation binding</td><td>132</td><td>6449</td><td>0.0008</td><td>0.033</td></tr><tr><td>GO:0051258</td><td>P</td><td>Protein polymerization</td><td>10</td><td>147</td><td>0.00018</td><td>0.013</td><td>GO:0043167</td><td>F</td><td>Ion binding</td><td>132</td><td>6472</td><td>0.00091</td><td>0.037</td></tr><tr><td>GO:0009086</td><td>P</td><td>Methionine biosynthetic process</td><td>5</td><td>30</td><td>0.00019</td><td>0.014</td><td>GO:0015145</td><td>F</td><td>Monosaccharide transmembrane transporter activity</td><td>6</td><td>67</td><td>0.00095</td><td>0.037</td></tr><tr><td>GO:0009628</td><td>P</td><td>Response to abiotic stimulus</td><td>82</td><td>3504</td><td>0.00023</td><td>0.016</td><td>GO:0016705</td><td>F</td><td>Oxidoreductase activity, acting on paired donors, with incorporation or reduction of molecular oxygen</td><td>16</td><td>415</td><td>0.0013</td><td>0.047</td></tr><tr><td>GO:0010114</td><td>P</td><td>Response to red light</td><td>10</td><td>154</td><td>0.00026</td><td>0.017</td><td>GO:0012505</td><td>C</td><td>Endomembrane system</td><td>98</td><td>3339</td><td>2.50E-09</td><td>1.10E-06</td></tr><tr><td>GO:0010876</td><td>P</td><td>Lipid localization</td><td>5</td><td>32</td><td>0.00025</td><td>0.017</td><td>GO:0009501</td><td>C</td><td>Amyloplast</td><td>7</td><td>25</td><td>5.00E-07</td><td>0.00011</td></tr><tr><td>GO:0019252</td><td>P</td><td>Starch biosynthetic process</td><td>5</td><td>35</td><td>0.00037</td><td>0.022</td><td>GO:0045298</td><td>C</td><td>Tubulin complex</td><td>7</td><td>33</td><td>2.50E-06</td><td>0.00035</td></tr><tr><td>GO:0009741</td><td>P</td><td>Response to brassinosteroid stimulus</td><td>13</td><td>258</td><td>0.00035</td><td>0.022</td><td>GO:0005618</td><td>C</td><td>Cell wall</td><td>37</td><td>1051</td><td>8.60E-06</td><td>0.00082</td></tr><tr><td>GO:0043436</td><td>P</td><td>Oxoacid metabolic process</td><td>49</td><td>1849</td><td>0.00037</td><td>0.022</td><td>GO:0031225</td><td>C</td><td>Anchored to membrane</td><td>23</td><td>501</td><td>9.50E-06</td><td>0.00082</td></tr><tr><td>GO:0019752</td><td>P</td><td>Carboxylic acid metabolic process</td><td>49</td><td>1849</td><td>0.00037</td><td>0.022</td><td>GO:0048046</td><td>C</td><td>Apoplast</td><td>15</td><td>285</td><td>8.00E-05</td><td>0.0057</td></tr><tr><td>GO:0006082</td><td>P</td><td>Organic acid metabolic process</td><td>49</td><td>1852</td><td>0.00038</td><td>0.023</td><td>GO:0030312</td><td>C</td><td>External encapsulating structure</td><td>37</td><td>1207</td><td>0.00014</td><td>0.0076</td></tr><tr><td>GO:0046274</td><td>P</td><td>Lignin catabolic process</td><td>5</td><td>36</td><td>0.00041</td><td>0.024</td><td>GO:0009505</td><td>C</td><td>Plant-type cell wall</td><td>30</td><td>896</td><td>0.00014</td><td>0.0076</td></tr><tr><td>GO:0010033</td><td>P</td><td>Response to organic substance</td><td>76</td><td>3252</td><td>0.00042</td><td>0.024</td><td>GO:0031224</td><td>C</td><td>Intrinsic to membrane</td><td>96</td><td>4289</td><td>0.0003</td><td>0.014</td></tr><tr><td>GO:0009725</td><td>P</td><td>Response to hormone stimulus</td><td>62</td><td>2522</td><td>0.00043</td><td>0.024</td><td>GO:0009705</td><td>C</td><td>Plant-type vacuole membrane</td><td>8</td><td>113</td><td>0.00062</td><td>0.024</td></tr><tr><td>GO:0050896</td><td>P</td><td>Response to stimulus</td><td>192</td><td>9863</td><td>0.00045</td><td>0.025</td><td>GO:0005576</td><td>C</td><td>Extracellular region</td><td>32</td><td>1072</td><td>0.00061</td><td>0.024</td></tr></tbody></table>

Ontology "P" indicates Biological Process, ontology "F" indicates Molecular Function, and ontology "C" indicates Cellular Component. "Number in input list" is the number of occurrence of the Gene Ontology term in the input list (668 genes). "Number in BG/Ref" is the observed occurrence of the Gene Ontology term in the Populus trichocarpa genome. P-value represents the statistical significance of the enrichment test and FDR is the false discovery rate set to 0.05.
668 genes significantly differentially expressed for GSEA, 104 GO terms had a significantly different frequency of occurrence and a wide range of BPs were affected (Table 1; see Supporting Information—Fig. S3). Among the 67 significant BP terms brought forward, 'hormone (ethylene and brassinosteroids) level response and regulation', 'cell wall biogenesis', 'amino acid metabolism', 'phenylpropanoid metabolism', 'lignin catabolic process', 'carbohydrate metabolism' and 'microtubule-based movements' were the most prominent. For the MF terms,

'ion binding' was the most over-represented term, followed by 'oxidoreductase activity' and 'transmembrane transport'. Several CC terms of interest were affected including 'membrane and endomembrane systems', 'cell wall', 'apoplast', 'amyloplast', 'vacuole' and 'tubulin complex' [see Supporting Information—Fig. S3].

Tables 2 and 3 list the significant GO terms for the 379 genes regulated upwards and the 262 genes regulated downwards, respectively. It is interesting to note that most

<table><caption>**Table 2.** GO term enrichment analysis of the 379 upregulated xylem genes (35° tilted vs. control).</caption><thead><tr><th>GO term</th><th>Ontology</th><th>Description</th><th>Number in input list</th><th>Number in BG/Ref</th><th>P-value</th><th>FDR</th></tr></thead><tbody><tr><td>GO:0005977</td><td>P</td><td>Glycogen metabolic process</td><td>5</td><td>52</td><td>0.00016</td><td>0.014</td></tr><tr><td>GO:0007017</td><td>P</td><td>Microtubule-based process</td><td>16</td><td>605</td><td>0.00018</td><td>0.014</td></tr><tr><td>GO:0006112</td><td>P</td><td>Energy reserve metabolic process</td><td>5</td><td>54</td><td>0.00019</td><td>0.015</td></tr><tr><td>GO:0009250</td><td>P</td><td>Glucan biosynthetic process</td><td>8</td><td>184</td><td>0.00036</td><td>0.026</td></tr><tr><td>GO:0009873</td><td>P</td><td>Ethylene-mediated signalling pathway</td><td>8</td><td>185</td><td>0.00037</td><td>0.026</td></tr><tr><td>GO:0009415</td><td>P</td><td>Response to water</td><td>15</td><td>584</td><td>0.00038</td><td>0.026</td></tr><tr><td>GO:0009741</td><td>P</td><td>Response to brassinosteroid stimulus</td><td>9</td><td>258</td><td>0.00074</td><td>0.049</td></tr><tr><td>GO:0009834</td><td>P</td><td>Secondary cell wall biogenesis</td><td>8</td><td>112</td><td>1.30E-05</td><td>0.0016</td></tr><tr><td>GO:0044262</td><td>P</td><td>Cellular carbohydrate metabolic process</td><td>28</td><td>1093</td><td>1.30E-06</td><td>0.00033</td></tr><tr><td>GO:0009409</td><td>P</td><td>Response to cold</td><td>26</td><td>855</td><td>1.40E-07</td><td>8.30E-05</td></tr><tr><td>GO:0044042</td><td>P</td><td>Glucan metabolic process</td><td>13</td><td>329</td><td>1.50E-05</td><td>0.0016</td></tr><tr><td>GO:0005975</td><td>P</td><td>Carbohydrate metabolic process</td><td>48</td><td>2018</td><td>1.50E-09</td><td>2.80E-06</td></tr><tr><td>GO:0006519</td><td>P</td><td>Cellular amino acid and derivative metabolic process</td><td>34</td><td>1688</td><td>1.60E-05</td><td>0.0016</td></tr><tr><td>GO:0051258</td><td>P</td><td>Protein polymerization</td><td>10</td><td>147</td><td>1.70E-06</td><td>0.00038</td></tr><tr><td>GO:0007018</td><td>P</td><td>Microtubule-based movement</td><td>13</td><td>174</td><td>1.70E-08</td><td>1.50E-05</td></tr><tr><td>GO:0009266</td><td>P</td><td>Response to temperature stimulus</td><td>27</td><td>1229</td><td>2.90E-05</td><td>0.0029</td></tr><tr><td>GO:0034637</td><td>P</td><td>Cellular carbohydrate biosynthetic process</td><td>17</td><td>483</td><td>3.40E-06</td><td>0.00068</td></tr><tr><td>GO:0044264</td><td>P</td><td>Cellular polysaccharide metabolic process</td><td>18</td><td>457</td><td>3.70E-07</td><td>0.00017</td></tr><tr><td>GO:0006575</td><td>P</td><td>Cellular amino acid derivative metabolic process</td><td>22</td><td>908</td><td>4.10E-05</td><td>0.0038</td></tr><tr><td>GO:0034641</td><td>P</td><td>Cellular nitrogen compound metabolic process</td><td>26</td><td>1031</td><td>4.10E-06</td><td>0.00073</td></tr><tr><td>GO:0042546</td><td>P</td><td>Cell wall biogenesis</td><td>12</td><td>313</td><td>4.30E-05</td><td>0.0038</td></tr><tr><td>GO:0009832</td><td>P</td><td>Plant-type cell wall biogenesis</td><td>12</td><td>252</td><td>5.40E-06</td><td>0.00083</td></tr><tr><td>GO:0006073</td><td>P</td><td>Cellular glucan metabolic process</td><td>13</td><td>298</td><td>5.50E-06</td><td>0.00083</td></tr><tr><td>GO:0005976</td><td>P</td><td>Polysaccharide metabolic process</td><td>21</td><td>639</td><td>6.90E-07</td><td>0.00022</td></tr><tr><td>GO:0000271</td><td>P</td><td>Polysaccharide biosynthetic process</td><td>14</td><td>354</td><td>7.10E-06</td><td>0.00098</td></tr><tr><td>GO:0016051</td><td>P</td><td>Carbohydrate biosynthetic process</td><td>21</td><td>641</td><td>7.30E-07</td><td>0.00022</td></tr><tr><td>GO:0033692</td><td>P</td><td>Cellular polysaccharide biosynthetic process</td><td>13</td><td>311</td><td>8.60E-06</td><td>0.0011</td></tr><tr><td>GO:0016758</td><td>F</td><td>Transferase activity, transferring hexosyl groups</td><td>18</td><td>707</td><td>0.00011</td><td>0.022</td></tr><tr><td>GO:0004553</td><td>F</td><td>Hydrolase activity, hydrolyzing O-glycosyl compounds</td><td>18</td><td>773</td><td>0.00032</td><td>0.043</td></tr><tr><td>GO:0016641</td><td>F</td><td>Oxidoreductase activity, acting on the CH-NH<sub>2</sub> group of donors, oxygen as acceptor</td><td>5</td><td>63</td><td>0.00037</td><td>0.043</td></tr><tr><td>GO:0005507</td><td>F</td><td>Copper ion binding</td><td>10</td><td>301</td><td>0.00055</td><td>0.048</td></tr><tr><td>GO:0016491</td><td>F</td><td>Oxidoreductase activity</td><td>48</td><td>3249</td><td>0.00057</td><td>0.048</td></tr><tr><td>GO:0005506</td><td>F</td><td>Iron ion binding</td><td>24</td><td>1268</td><td>0.00068</td><td>0.05</td></tr><tr><td>GO:0016638</td><td>F</td><td>Oxidoreductase activity, acting on the CH-NH<sub>2</sub> group of donors</td><td>7</td><td>94</td><td>3.60E-05</td><td>0.011</td></tr><tr><td>GO:0005200</td><td>F</td><td>Structural constituent of cytoskeleton</td><td>7</td><td>73</td><td>7.90E-06</td><td>0.0046</td></tr><tr><td>GO:0030312</td><td>C</td><td>External encapsulating structure</td><td>24</td><td>1207</td><td>0.00034</td><td>0.014</td></tr><tr><td>GO:0016020</td><td>C</td><td>Membrane</td><td>111</td><td>9319</td><td>0.00065</td><td>0.024</td></tr><tr><td>GO:0031224</td><td>C</td><td>Intrinsic to membrane</td><td>59</td><td>4289</td><td>0.00076</td><td>0.025</td></tr><tr><td>GO:0044425</td><td>C</td><td>Membrane part</td><td>71</td><td>5441</td><td>0.00087</td><td>0.025</td></tr><tr><td>GO:0044430</td><td>C</td><td>Cytoskeletal part</td><td>17</td><td>790</td><td>0.0011</td><td>0.029</td></tr><tr><td>GO:0031225</td><td>C</td><td>Anchored to membrane</td><td>21</td><td>501</td><td>1.40E-08</td><td>3.40E-06</td></tr><tr><td>GO:0012505</td><td>C</td><td>Endomembrane system</td><td>63</td><td>3339</td><td>2.30E-08</td><td>3.40E-06</td></tr><tr><td>GO:0005618</td><td>C</td><td>Cell wall</td><td>24</td><td>1051</td><td>4.50E-05</td><td>0.0033</td></tr><tr><td>GO:0045298</td><td>C</td><td>Tubulin complex</td><td>7</td><td>33</td><td>6.30E-08</td><td>6.10E-06</td></tr><tr><td>GO:0048046</td><td>C</td><td>Apoplast</td><td>11</td><td>285</td><td>8.40E-05</td><td>0.0047</td></tr><tr><td>GO:0009505</td><td>C</td><td>Plant-type cell wall</td><td>21</td><td>896</td><td>9.60E-05</td><td>0.0047</td></tr></tbody></table>

Ontology "P" indicates Biological Process, ontology "F" indicates Molecular Function, and ontology "C" indicates Cellular Component. "Number in input list" is the number of occurrence of the Gene Ontology term in the input list (379 genes). "Number in BG/Ref" is the observed occurrence of the Gene Ontology term in the *Populus trichocarpa* genome. P-value represent the statistical significance of the enrichment test and FDR is the false discovery rate set to 0.05.
<table><caption>Table 3. GO term enrichment analysis of the 262 downregulated xylem genes (35° tilted vs. control).</caption><thead><tr><th>GO term</th><th>Ontology</th><th>Description</th><th>Number in input list</th><th>Number in BG/Ref</th><th>P-value</th><th>FDR</th></tr></thead><tbody><tr><td>GO:0010876</td><td>P</td><td>Lipid localization</td><td>5</td><td>32</td><td>2.80E-06</td><td>0.0037</td></tr><tr><td>GO:0009067</td><td>P</td><td>Aspartate family amino acid biosynthetic process</td><td>6</td><td>75</td><td>9.70E-06</td><td>0.0065</td></tr><tr><td>GO:0034641</td><td>P</td><td>Cellular nitrogen compound metabolic process</td><td>19</td><td>1031</td><td>2.20E-05</td><td>0.0088</td></tr><tr><td>GO:0042221</td><td>P</td><td>Response to chemical stimulus</td><td>52</td><td>4909</td><td>3.40E-05</td><td>0.0088</td></tr><tr><td>GO:0009066</td><td>P</td><td>Aspartate family amino acid metabolic process</td><td>7</td><td>145</td><td>3.90E-05</td><td>0.0088</td></tr><tr><td>GO:0015698</td><td>P</td><td>Inorganic anion transport</td><td>7</td><td>144</td><td>3.80E-05</td><td>0.0088</td></tr><tr><td>GO:0000097</td><td>P</td><td>Sulfur amino acid biosynthetic process</td><td>5</td><td>67</td><td>7.50E-05</td><td>0.013</td></tr><tr><td>GO:0019748</td><td>P</td><td>Secondary metabolic process</td><td>20</td><td>1235</td><td>7.60E-05</td><td>0.013</td></tr><tr><td>GO:0009308</td><td>P</td><td>Amine metabolic process</td><td>19</td><td>1171</td><td>0.00011</td><td>0.017</td></tr><tr><td>GO:0050896</td><td>P</td><td>Response to stimulus</td><td>85</td><td>9863</td><td>0.00013</td><td>0.017</td></tr><tr><td>GO:0006555</td><td>P</td><td>Methionine metabolic process</td><td>5</td><td>88</td><td>0.00025</td><td>0.031</td></tr><tr><td>GO:0044106</td><td>P</td><td>Cellular amine metabolic process</td><td>15</td><td>890</td><td>0.00041</td><td>0.042</td></tr><tr><td>GO:0006520</td><td>P</td><td>Cellular amino acid metabolic process</td><td>15</td><td>890</td><td>0.00041</td><td>0.042</td></tr><tr><td>GO:0009755</td><td>P</td><td>Hormone-mediated signalling pathway</td><td>16</td><td>1003</td><td>0.00047</td><td>0.046</td></tr><tr><td>GO:0006820</td><td>P</td><td>Anion transport</td><td>8</td><td>297</td><td>0.00055</td><td>0.049</td></tr><tr><td>GO:0016491</td><td>F</td><td>Oxidoreductase activity</td><td>43</td><td>3249</td><td>9.10E-07</td><td>0.00041</td></tr><tr><td>GO:0009055</td><td>F</td><td>Electron carrier activity</td><td>24</td><td>1346</td><td>2.90E-06</td><td>0.00066</td></tr><tr><td>GO:0046906</td><td>F</td><td>Tetrapyrrole binding</td><td>17</td><td>819</td><td>1.40E-05</td><td>0.0018</td></tr><tr><td>GO:0020037</td><td>F</td><td>Heme binding</td><td>16</td><td>744</td><td>1.60E-05</td><td>0.0018</td></tr><tr><td>GO:0015103</td><td>F</td><td>Inorganic anion transmembrane transporter activity</td><td>7</td><td>154</td><td>5.70E-05</td><td>0.0051</td></tr><tr><td>GO:0005506</td><td>F</td><td>Iron ion binding</td><td>20</td><td>1268</td><td>0.00011</td><td>0.008</td></tr><tr><td>GO:0008509</td><td>F</td><td>Anion transmembrane transporter activity</td><td>8</td><td>251</td><td>0.00018</td><td>0.012</td></tr><tr><td>GO:0016209</td><td>F</td><td>Antioxidant activity</td><td>8</td><td>297</td><td>0.00055</td><td>0.03</td></tr><tr><td>GO:0022891</td><td>F</td><td>Substrate-specific transmembrane transporter activity</td><td>23</td><td>1842</td><td>0.00089</td><td>0.044</td></tr><tr><td>GO:0019825</td><td>F</td><td>Oxygen binding</td><td>9</td><td>409</td><td>0.001</td><td>0.045</td></tr><tr><td>GO:0016684</td><td>F</td><td>Oxidoreductase activity, acting on peroxide as acceptor</td><td>7</td><td>262</td><td>0.0013</td><td>0.046</td></tr><tr><td>GO:0042803</td><td>F</td><td>Protein homodimerization activity</td><td>7</td><td>268</td><td>0.0014</td><td>0.046</td></tr><tr><td>GO:0008289</td><td>F</td><td>Lipid binding</td><td>11</td><td>606</td><td>0.0014</td><td>0.046</td></tr><tr><td>GO:0004601</td><td>F</td><td>Peroxidase activity</td><td>7</td><td>262</td><td>0.0013</td><td>0.046</td></tr></tbody></table>

Ontology "P" indicates Biological Process, ontology "F" indicates Molecular Function, and ontology "C" indicates Cellular Component. "Number in input list" is the number of occurrence of the Gene Ontology term in the input list (262 genes). "Number in BG/Ref" is the observed occurrence of the Gene Ontology term in the *Populus trichocarpa* genome. P-value represent the statistical significance of the enrichment test and FDR is the false discovery rate set to 0.05.

enriched GO terms are exclusive to genes regulated upwards or downwards. Since 27 genes with opposite tendencies were included in either group, they were analysed separately (Table 4).

By looking more precisely at CC terms impacted by the gravitational stimulation treatment, we discovered how the cells integrate the stimulus. Figure 4 is a graphical analysis of the CC term enrichment related to modulated genes. XU-specific genes appeared to be significantly linked to the 'cell wall' while the 'plasma membrane' was the only significant XL-specific CC. XU- and XL-modulated genes covered the 'apoplast' and the 'endomembrane system'. When considering genes according to their expression tendency, no significant CC was identified for the 262 genes regulated downstream. The greatest variety of CC was found in upregulated genes including: 'apoplast', 'cell wall', 'tubulin complex', 'amyloplast', 'endomembrane and membrane anchored'. For the 27 genes with upwards or downwards regulation on the xylem side, the 'cell wall' compartment was significantly impacted, in the same way as the specific XU genes (Fig. 4). Finally, it is interesting to note that among the genes expressed differently after gravistimulation, 36 were annotated as putative transcription factor (TF) [see Supporting Information—Table S4].

Exploration of candidate genes

The experiments on DNA chips identified 668 genes with significant differential expression after 30 min of inclination at 35°. A subset (153 genes) was used to further investigate their involvement in gravitropic response by high-throughput quantitative PCR in a kinetic study (0, 30, 60, 120 and 180 min at 35°). As a prerequisite, we ensured the results of the biochips by reanalysing the expression of candidate genes found in the biochips experiment on five independent biological replicates. We found significant pairwise correlations between the two approaches (i.e. 0.67, 0.50, 0.71, 0.72, 0.67) allowing for a reliable interpretation of the data produced subsequently. We then performed a hierarchical grouping within the inclination expression data that led to the identification of five co-expression groups identified thereafter as Cluster #1 to Cluster #5 (Fig. 5).

Cluster #1 identified 19 genes mainly downregulated in XL after gravistimulation. A large bulk (8/19) concerned the metabolism of sugars (glycosyl groups and/or carbohydrates), whereas another group of five genes had something to do with auxin and/or ethylene responses and/or metabolism. Lastly, four genes were involved in transmembrane transport. Cluster #2 comprised the highest number of genes (81) and a wide panel of metabolic processes. Using a KEGG enrichment analysis, three
<table><caption>**Table 4.** GO term enrichment analysis of the 27 differentially expressed genes with opposite trends in XU and XL (35° tilted vs. control).</caption><thead><tr><th>GO term</th><th>Ontology</th><th>Description</th><th>Number in input list</th><th>Number in BG/Ref</th><th>P-value</th><th>FDR</th></tr></thead><tbody><tr><td>GO:0044042</td><td>P</td><td>Glucan metabolic process</td><td>6</td><td>329</td><td>6.90E-08</td><td>1.60E-06</td></tr><tr><td>GO:0006073</td><td>P</td><td>Cellular glucan metabolic process</td><td>6</td><td>298</td><td>3.90E-08</td><td>1.60E-06</td></tr><tr><td>GO:0044264</td><td>P</td><td>Cellular polysaccharide metabolic process</td><td>6</td><td>457</td><td>4.60E-07</td><td>7.00E-06</td></tr><tr><td>GO:0005976</td><td>P</td><td>Polysaccharide metabolic process</td><td>6</td><td>639</td><td>3.10E-06</td><td>3.60E-05</td></tr><tr><td>GO:0044262</td><td>P</td><td>Cellular carbohydrate metabolic process</td><td>7</td><td>1093</td><td>5.10E-06</td><td>4.60E-05</td></tr><tr><td>GO:0005975</td><td>P</td><td>Carbohydrate metabolic process</td><td>7</td><td>2018</td><td>0.00025</td><td>0.0019</td></tr><tr><td>GO:0016758</td><td>F</td><td>Transferase activity, transferring hexosyl groups</td><td>5</td><td>707</td><td>8.90E-05</td><td>0.0033</td></tr><tr><td>GO:0016757</td><td>F</td><td>Transferase activity, transferring glycosyl groups</td><td>5</td><td>954</td><td>0.00036</td><td>0.0066</td></tr><tr><td>GO:0005506</td><td>F</td><td>Iron ion binding</td><td>5</td><td>1268</td><td>0.0013</td><td>0.016</td></tr><tr><td>GO:0005618</td><td>C</td><td>Cell wall</td><td>5</td><td>1051</td><td>0.00055</td><td>0.019</td></tr><tr><td>GO:0030312</td><td>C</td><td>External encapsulating structure</td><td>5</td><td>1207</td><td>0.001</td><td>0.019</td></tr></tbody></table>

Ontology "P" indicates Biological Process, ontology "F" indicates Molecular Function, and ontology "C" indicates Cellular Component. "Number in input list" is the number of occurrence of the Gene Ontology term in the input list (27 genes). "Number in BG/Ref" is the observed occurrence of the Gene Ontology term in the *Populus trichocarpa* genome. P-value represent the statistical significance of the enrichment test and FDR is the false discovery rate set to 0.05.

major pathways were brought forward including photosystem II oxygen-evolving enhancer protein 1 (i.e. OEE), the SAUR protein family (i.e. auxin responsive genes) and β-tubulin (i.e. cytoskeleton). Briefly, cell redox players, putative mechanical sensors, various kinases, several brassinoid, ethylene- and auxin-related genes and TFs were well-represented in that cluster. It is also in Cluster #2 that we found genes already known to be involved in lignin synthesis and cell wall modification. Cluster #3 identified three Fasciclin-like arabinoglycan genes, expression of which was bimodal in XU and XL reaching two maxima at 60 and 180 min after gravistimulation. Cluster #4 was the second largest cluster with 42 genes over-expressed in XL and whose expression was not modified in XU. Functional annotations highlighted phenylpropanoid, flavonoid and terpenoid biosynthesis. Of note within that cluster was the significant BP enrichment in 'transmembrane transport' especially concerning sugars. 'Lipid metabolic process' was the second biological process pinpointed. Finally, Cluster #5 was the smallest cluster identified with only three co-expressed genes upregulated in both XL and XU. Functional annotations suggested that oxidative stress and oligopeptide mobilization were the processes concerned.

## Discussion

### The experimental device allowed to identify genes regulated by inclination

This study aimed to better understand the molecular processes that occur in the secondary xylem of the woody stem after short-term tilting while separating the phototropic effect from gravitrophic stimulation. To do so, we used an innovative experimental design (Coutand et al. 2019) to identify a list of differentially expressed genes specifically responding to gravitropism without interference with the phototropic response. Stem staking in our condition was shown to be effective in preventing stem deformation as assessed by the absence of induction of PtaZFP2, a well-characterized reporter gene for stem bending in poplar (Martin et al. 2009). It is interesting to note that we were able to record very early molecular responses, whereas previous studies, which did not use our innovative device and are therefore difficult to compare, found significant differential gene expression only after 8 h of gravistimulation (Zinkgraft et al. 2018). It is difficult to say whether this difference is due to the distinctive statistical power of the studies or whether it is due

to other factors. The idea that an early molecular adjustment could occur in the gravitropic response is however supported by our earlier observation that changes in the poplar stem proteome were recorded as early as 45 min after tilting (Azri et al. 2009).

### Genes supporting cell wall reorganization are modulated in 'opposite and tension wood' forming zones

Among the many gene classes identified in this study, we have chosen to first focus on the genes involved in cell wall formation and modification. The cell component 'cell wall' was indeed one of the CCs of interest and this choice was strongly supported by the biological process terms that were significantly brought forward in our analysis including: 'cell wall biogenesis', 'carbohydrate metabolic process' and 'catabolic lignin process'. Within that compartment, our experimental set-up allowed the identification of genes involved in wood formation specifically regulated by inclination independently of the phototropic stimulation response. It was interesting to note that as early as after 30 min of gravistimulation modulation of gene expression was recorded on both sides of the inclined stem (Fig. 3), thus concerning the whole secondary xylem. From the analysis of the GO enrichment (Fig. 4), it appears that the changes mainly affect the cell wall compartment on the upper surface, which corresponds to the installation of G-layer fibres (Jourez et al. 2001). For the lower xylem, the plasma membrane is the significant cellular component affected, which could indicate the importance of control between adjacent compartments. Some genes, encoding mainly TFs involved in the formation of the secondary cell wall, appear to be regulated in both XU and XL. This is the case with Potri.001G112200, a KNAT7 homologue (class II KNOX gene) that has been shown elsewhere to be a key negative regulator of secondary cell wall formation (Li et al. 2012). Bhargava et al. (2013) suggested that KNAT7 could be involved in the suppression of lignin biosynthesis in stems. The analysis of promoters of genes encoding such TF sensitive to gravitational stimulation could provide a good opportunity to decipher molecular events upstream.

Rose et al. (2004) suggested that cell wall remodelling involves a complex interplay of enzymatic and non-enzymatic oxidative and hydrolytic changes. In agreement with this we found in Cluster #2 several redox-regulating players including: four oxygen-evolving enhancer proteins, a glutathione S-transferase, thioredoxin h, three 2OG-Fe(II) oxygenases and a peroxydase.
<figure><img src="image_6.png" /><figcaption>**Figure 4.** Graphical results of AgriGO cellular component terms analysis for differentially expressed genes put in different subgroups. Upper xylem specific (XU: 214 genes, 184 annotated), lower xylem specific (XL: 250 genes, 206 annotated), Genes modulated in both sides of xylem (204 genes, 181 annotated), Up-regulated genes (379 genes, 352 annotated), Up- and down-regulated (27 genes, 24 annotated). No significant cellular component term was found in the down-regulated subset. Relationships between terms are materialised by arrows described in the key. Coloured terms indicate a statistical signification: the more significant a term is, the redder the colour is. Term frequency found in *Populus trichocarpa* genome is indicated at the lower right corner and is used for statistical significance test.</figcaption></figure>

We also found within that same cluster eight carbohydrate active enzymes (CAZymes), four glycosyltransferases (GTs) and four glycoside hydrolases (GHs) that could participate in

modulating the composition and structure of gravistimulated xylem cell wall. The fact that they were all clustered together supports a synergistic effect on cell wall properties as suggested
<figure><img src="image_7.png" /><figcaption>Figure 5. qPCR heatmap + clusters.</figcaption></figure>

earlier by Tenhaken (2015). In addition, the observation that these genes also clustered with genes involved in regulation of cell expansion (namely Potri.009G067100, Potri.009G067200 and Potri.008G16280) suggests that they all may be mobilized during cell wall modifications accompanying cell expansion.

A well-established process is that of the reprogrammation of sugar metabolism when normal wood is compared to TW (Mizrachi et al. 2014). Several reports (Lafarguette et al. 2004;

Paux et al. 2005; Andersson-Gunneras et al. 2006) showed that in TW (upper side), the upregulation of cellulose synthase activity was responsible for the production of the so-called gelatinous layer (G-layer). In our experimental setting, reprogrammation of sugar metabolism appears to occur on both sides of the stem shortly after gravistimulation as evidenced by the increased expression of two genes encoding cellulose synthases in XL and XU (CESA4 Potri.004G059600, Potri.002G257900). In addition,
the expression of two cellulose synthase-like genes (CSL, β-1,4-mannan synthases) was downregulated in the two xylem sides. In poplar it was shown that Fasciclin-like arabinogalactan proteins (FLA) were induced after gravistimulation and thought to participate in TW G-layer synthesis (Lafarguette et al. 2004; Fagerstedt et al. 2014; MacMillan et al. 2015; Wang et al. 2015). We show here that a bulk of 10 FLA were upregulated after 30 min of gravistimulation. Cluster #1 also brought forward genes involved in sugar transport/removal/hydrolysis (Potri.005G260300, and the CAZymes: Potri.001G028100, Potri.004G044100, Potri.004G113500, Potri.002G080700, Potri.003G112600) suggesting that sugar repartition is modified in xylem after inclination.

In agreement with the observation made by Villalobos et al. (2012) that 'opposite wood' (lower side) has a greater lignin content than TW we found four upregulated LIM TFs (Potri.002G157300, Potri.002G118000, Potri.009G087200, Potri.014G080900) that were suspected to play a role in lignin biosynthesis (Chen et al. 2015). Our analysis also brought forward genes involved in monolignols production such as phenylalanine ammonia-lyase (PtPAL1.2, Potri.006G126800) expression of which was repressed in XL and, 4-coumarate:CoA ligase (4CL, Potri.003G188500) repressed in XU along with caffeoyl-CoA O-methyltransferase 1 (CCoAOMT1, Potri.009G099800). In addition, Potri.016G112100, a laccase 1 homolog was found upregulated in XU although it is not yet consensually admitted that it participates in the lignification process (Gavnholt and Larsen 2002; Berthet et al. 2012). Consistent with the idea that TW has a lower lignin content we found that PtSAMS1 (Potri.013G004100), a S-adenosylmethionine synthetase (SAMS) suspected to play a central role in lignin biosynthesis (Vander Mijnsbrugge et al. 2000; Shen et al. 2002), was inhibited in XU.

Finally, we also found that two sequences (Potri.011G006800 and Potri.004G009600) encoding two T-proteins of the glycine decarboxylase complex (GDC) participating in serine biosynthesis were downregulated in XU. This is in line with the suggestion that GDC may allow woody perennials to cope with C1 demands of lignification (Rajinikanth et al. 2007). It is also consistent with the report of GDC involvement in the context of compression wood formation (Villalobos et al. 2012).

## Impact of gravistimulation on wood cell expansion

Some of the significant BP, MF and CC terms brought forward by the present study strongly suggest that wood cell expansion is a process impacted by gravistimulation. This is in agreement with a growing number of studies linking gibberellins and TW (Eriksson et al. 2000; Mauriat and Moritz 2009; Nugroho et al. 2012, 2013; for a review, see Felten and Sundberg 2013). The prevailing view, as reviewed recently by Ruelle (2014), is that in poplar TW fibres are longer and thinner when compared to that in 'opposite wood', suggesting complex regulation of cell expansion processes. In agreement with these observations, we found that Potri.012G132400 and Potri.015G134600, two active gibberellins (GA20OX) coding genes that have been linked elsewhere with cell expansion were modulated (Israelsson et al. 2005). Besides gibberellins (GA), other phytohormones are known to play a role in cell expansion. In particular, brassinosteroids (BR) were shown to be involved in stem reorientation and cell elongation (Bienert et al. 2012; Vandenbussche et al. 2013). Interestingly, we found that a homolog to BRI1 (brassinosteroid-insensitive1, an A. thaliana cell surface LRR-S/T kinase receptor for BR) (here identified as Potri.007G078100), and a serine carboxypeptidase (SCP) homolog to brassinosteroid insensitive1 suppressor (BRS1) (here identified as Potri.011G046600) were found modulated within Cluster #2 in response to gravistimulation. Of course, we

cannot rule out the hypothesis that the listed genes involved in detection/response to GA and BR have a role in the regulation of cell differentiation and secondary cell wall biogenesis as suggested by Du et al. (2019). In addition, GA may regulate key genes involved in BR response (Gerttula et al. 2015).

As cell expansion is regulated by the plasticity of the primary cell wall (Cosgrove 2005) it was consistent to see that after our stimuli several genes encoding proteins involved in loosening tension between primary cell wall constituents were found modulated, including expansin, expansin-like and pectin xyloglucans modifier genes (Potri.007G083400, Potri.010G167200, Potri.013G134300, Potri.002G060500, Potri.013G005700, Potri.014G115000, Potri.018G084300).

At last, supporting the involvement of cell expansion in response to gravistimulation we recorded the upregulation of Potri.004G117200 in XL. This sequence is a homolog of the A. thaliana COBRA protein, an extracellular glycosyl-phosphatidyl inositol-anchored protein required for polarized cell expansion through cellulose microfibril orientation (Schindelman et al. 2001; Roudier et al. 2005).

## Could concerted regulations of flavonoid and phosphoinositides pathways occur after gravistimulation?

Following gravistimulation, signal transduction leads to complex biochemical and genomic responses at the subcellular level. In Cluster #4, we observed the co-expression of genes involved in phenylpropanoid biosynthesis (PAL1—Potri.006G126800) and flavonoid (Potri.002G033600, Potri.003G100200) pathways. The TF LIM1 (Potri.002G118000) is also found in this cluster which agrees with Chen et al. (2015) who have found co-expression in 'opposite wood' of several LIM TFs with PAL genes. Several studies suggest that the synthesis or deposition of flavonoids regulates auxin transport during gravitropic responses (Geisler et al. 2005; Buer et al. 2007, 2013; Santelia et al. 2008; Lewis et al. 2011). Gayomba et al. (2016) demonstrated their role in auxin transport inhibition and antioxidant activity in Arabidopsis.

Furthermore, it is known that phosphoinositides regulate vesicle trafficking (for a review, see Krishnamoorthy et al. 2014) and auxin-dependent development (Tejos et al. 2014). Phosphatidyl inositol transfer protein SEC14 (Potri.004G157600) is co-expressed in Cluster #4 with CLATHRIN assembly protein (Potri.006G066900) indicating the regulation of vesicle trafficking. Lipid metabolism is also regulated (Potri.006G166800, Potri.012G12390).

In this context we hypothesize that concerted regulations of the flavonoid and phosphoinositides pathways led to auxin modulation after gravistimulation, probably through vesicle trafficking.

In the research areas of nutrition and dietetics, it has been shown that flavonoids exert modulating actions on the cellular system through direct action on various signalling pathways such as those of phosphoinositide 3-kinase, and protein kinase C (Mansuri et al. 2014).

## What about cytoskeleton involvement in xylem following gravistimulation?

On the understanding that the involvement of the cytoskeleton following shoot gravistimulation is still an unsolved question (Song et al. 2019), we noticed that CC terms such as 'membrane and endomembrane systems', 'cell wall', 'apoplast', 'amyloplast', 'vacuole' and 'tubulin complex' were affected in our gravistimulation experiment. The observation that the expression
of tubulin genes, as well as that of an actin-depolymerizing factor (Potri.001G106200) and a F-actin capping protein alpha subunit (Potri.013G017000) is modulated suggests that modification of the cell cytoskeleton is partly a response to gravistimulation as suggested earlier by Nick (2011). It was shown already in *Arabidopsis* stem that the actin microfilaments may be considered as negative regulators of the gravitropic response considering their interaction with the amyloplasts through the link to SGR9 (Kolesnikov *et al.* 2016). Song *et al.* (2019) highlighted the role of the actin-bonding protein Rice Morphology Determinant (RMD) in the regulation of the gravitropic response of the light-grown shoots. In different studies, microtubules were also shown to participate in gravitropism mainly via cellulose deposition in xylem (Fagerstedt *et al.* 2014; Araniti *et al.* 2016).

To conclude, our experimental isotropic light experimental set-up allowing the strict analysis of gravitropic stimulation of trees not polluted by the phototropic response highlighted a list of 668 xylem-regulated genes involved in early gravitropic response. Some of the genes pinpointed here were already reported to participate in gravity response confirming the pertinence of our approach and attesting that they were indeed related to the gravity stimuli and not to another confounding factor (phototropic response, stem bending response). Interestingly, more than a third of these genes are still not annotated for gene function. We consider that this particular gene subset, especially those still not annotated, constitutes a good starting point to expand our current knowledge on gravitropism. This work opens new research prospects towards the understanding of events triggering 'reaction wood' formation.

## Supporting Information

The following additional information is available in the online version of this article—

Table S1. List of qPCR primers.

Table S2. List of the 2355 differentially expressed genes based on poplar annotation V1.

Table S3. List of the 668 differentially expressed genes based on poplar annotation V3.

Table S4. List of transcription factors.

Figure S1. Correlation analysis of microarray data.

Figure S2. High-resolution version of Fig. 3 with gene model names.

Figure S3. GO term enrichment of Biological Processes (BP), Molecular Functions (MF) and Cellular Components (CC) of the 668 differentially expressed genes in XU and XL.

## Sources of Funding

This work was supported by the 'Agence Nationale de la Recherche' (ANR-11-BSV7-0012).

## Conflict of Interest

None declared.

## Contributions by the Authors

D.L. and J.F. performed the experiments and collected data. DL and PRD designed the study, analysed the data and wrote the first draft. J.S.V., P.L., C.C. and J.R.D. critically reviewed the manuscript before final approval.

## Acknowledgements

We thank Brigitte Girard, Christelle Boisselet and Sylvaine Blateyron for the production and care of the plants, RNA extraction and their help during the measurements. We would like to thank Stéphane Ploquin for his technical skills in improving the spheres. We thank the Plant Genomics Research Unit (URGV) platform for their advice.

## Literature Cited

Andersson-Gunnerås S, Mellerowicz EJ, Love J, Segerman B, Ohmiya Y, Coutinho PM, Nilsson P, Henrissat B, Moritz T, Sundberg B. 2006. Biosynthesis of cellulose-enriched tension wood in *Populus*: global analysis of transcripts and metabolites identifies biochemical and developmental regulators in secondary wall biosynthesis. *The Plant Journal* 45:144–165.

Araniti F, Graña E, Krasuska U, Bogatek R, Reigosa MJ, Abenavoli MR, Sánchez-Moreiras AM. 2016. Loss of gravitropism in farnesene-treated *Arabidopsis* is due to microtubule malformations related to hormonal and ROS unbalance. *PLoS One* 11:e0160202.

Azri W, Brunel N, Franchel J, Ben Rejeb I, Jacquot JP, Julien JL, Herbette S, Roeckel-Drevet P. 2013. Putative involvement of Thioredoxin h in early response to gravitropic stimulation of poplar stems. *Journal of Plant Physiology* 170:707–711.

Azri W, Chambon C, Herbette S, Brunel N, Coutand C, Leplé JC, Ben Rejeb I, Ammar S, Julien JL, Roeckel-Drevet P. 2009. Proteome analysis of apical and basal regions of poplar stems under gravitropic stimulation. *Physiologia Plantarum* 136:193–208.

Baldwin KL, Strohm AK, Masson PH. 2013. Gravity sensing and signal transduction in vascular plant primary roots. *American Journal of Botany* 100:126–142.

Bastien R, Douady S, Moulia B. 2015. A unified model of shoot tropism in plants: photo-, gravi- and propio-ception. *PLoS Computational Biology* 11:e1004037.

Benjamini Y, Hochberg Y. 2000. On the adaptative control of false discovery rate in multiple testing with independent statistics. *Journal of Educational and Behavioral Statistics* 25:60–83.

Benjamini Y, Yekutieli D. 2001. The control of the false discovery rate in multiple testing under dependency. *Annals of Statistics* 29:1165–1188.

Berthet S, Thévenin J, Baratiny D, Demont-Caulet N, Debeaujon I, Bidzinski P, Leplé J-C, Huis R, Hawkins S, Gomez LD, Lapierre C, Jouanin L. 2012. Role of plant laccases in lignin polymerization. Chapter 5. In: Jacquot J-P, Gadal P, Jouanin L, Lapierre C, eds. *Lignins—biosynthesis, biodegradation and bioengineering. advances in botanical research*. London: Academic Press, Elsevier, 145–172.

Bhargava A, Ahad A, Wang S, Mansfield SD, Haughn GW, Douglas CJ, Ellis BE. 2013. The interacting MYB75 and KNAT7 transcription factors modulate secondary cell wall deposition both in stems and seed coat in *Arabidopsis*. *Planta* 237:1199–1211.

Bienert MD, Delannoy M, Navarre C, Boutry M. 2012. NtSCP1 from tobacco is an extracellular serine carboxypeptidase III that has an impact on cell elongation. *Plant Physiology* 158:1220–1229.

Buer CS, Kordbacheh F, Truong TT, Hocart CH, Djordjevic MA. 2013. Alteration of flavonoid accumulation patterns in transparent testa mutants disturbs auxin transport, gravity responses, and imparts long-term effects on root and shoot architecture. *Planta* 238:171–189.

Buer CS, Muday GK, Djordjevic MA. 2007. Flavonoids are differentially taken up and transported long distances in *Arabidopsis*. *Plant Physiology* 145:478–490.

Chen J, Chen B, Zhang D. 2015. Transcript profiling of *Populus tomentosa* genes in normal, tension, and opposite wood by RNA-seq. *BMC Genomics* 16:164.

Cosgrove DJ. 2005. Growth of the plant cell wall. *Nature Reviews Molecular Cell Biology* 6:850–861.

Côté WAJ, Day AC. 1964. *Anatomy and ultrastructure of reaction wood*. New York: Syracuse University Press.

Coutand C. 2010. Mechanosensing and thigmomorphogenesis, a physiological and biomechanical point of view. *Plant Science* 179:168–182.
Coutand C, Adam B, Ploquin S, Moulia B. 2019. A method for the quantification of phototropic and gravitropic sensitivities of plants combining an original experimental device with model-assisted phenotyping: exploratory test of the method on three hardwood tree species. *PLoS One* 14:e0209973.

Du J, Gerttula S, Li Z, Zhao ST, Liu YL, Liu Y, Lu MZ, Groover AT. 2019. Brassinosteroid regulation of wood formation in poplar. *The New Phytologist* 225:1516–1530.

Du S, Yamamoto F. 2007. An overview of the biology of reaction wood formation. *Journal of Integrative Plant Biology* 49:131–143.

Eriksson ME, Israelsson M, Olsson O, Moritz T. 2000. Increased gibberellin biosynthesis in transgenic trees promotes growth, biomass production and xylem fiber length. *Nature Biotechnology* 18:784–788.

Fagerstedt KV, Mellerowicz E, Gorshkova T, Riel K, Joseleau J-P. 2014. Cell wall polymers in reaction wood. In: Gardiner B, Barnett J, Saranpää P, Joseph G, eds. *The biology of reaction wood*. Montpellier: Springer (Springer series in wood science), 37–106.

Felten J, Sundberg B. 2013. Biology, chemistry and structure of tension wood. In: Fromm J, ed. *Cellular aspects of wood formation*. Plant Cell Monographs 20. Berlin, Heidelberg: Springer-Verlag, 203–224.

Gavnholt B, Larsen K. 2002. Molecular biology of plant laccases in relation to lignin formation. *Physiologia Plantarum* 116:273–280.

Gayomba SR, Watkins JM, Muday GK. 2016. Flavonols regulate plant growth and development through regulation of auxin transport and cellular redox status. *Recent Advances in Polyphenol Research* 5:143–170.

Geisler M, Blakeslee JJ, Bouchard R, Lee OR, Vincenzetti V, Bandyopadhyay A, Titapiwatanakun B, Peer WA, Bailly A, Richards EL, Ejendal KF, Smith AP, Baroux C, Grossniklaus U, Müller A, Hrycyna CA, Dudler R, Murphy AS, Martinoia E. 2005. Cellular efflux of auxin catalyzed by the *Arabidopsis* MDR/PGP transporter AtPGP1. *The Plant Journal* 44:179–194.

Gerttula S, Zinkgraf M, Muday GK, Lewis DR, Ibatullin FM, Brumer H, Hart F, Mansfield SD, Filkov V, Groover A. 2015. Transcriptional and hormonal regulation of gravitropism of woody stems in *Populus*. *The Plant Cell* 27:2800–2813.

Groover A. 2016. Gravitropisms and reaction woods of forest trees - evolution, functions and mechanisms. *The New Phytologist* 211:790–802.

Hashiguchi Y, Tasaka M, Morita MT. 2013. Mechanism of higher plant gravity sensing. *American Journal of Botany* 100:91–100.

Israelsson M, Sundberg B, Moritz T. 2005. Tissue-specific localization of gibberellins and expression of gibberellin-biosynthetic and signaling genes in wood-forming tissues in aspen. *The Plant Journal* 44:494–504.

Jourez B, Riboux A, Leclercq A. 2001. Anatomical characteristics of tension wood and opposite wood in young inclined stems of poplar (*Populus Euramerica* cv "Ghoy"). *IAWA Journal* 22:133–157.

Kolesnikov YS, Kretynin SV, Volotovsky ID, Kordyum EL, Ruelland E, Kravets VS. 2016. Molecular mechanisms of gravity perception and signal transduction in plants. *Protoplasma* 253:987–1004.

Krishnamoorthy P, Sanchez-Rodriguez C, Heilmann I, Persson S. 2014. Regulatory roles of phosphoinositides in membrane trafficking and their potential impact on cell-wall synthesis and re-modelling. *Annals of Botany* 114:1049–1057.

Kubler H. 1959. Studies on growth stresses in trees. 1. The origin of growth stresses and the stresses in transverse direction. *Holz als Roh- und Werkstoff* 17:1–9.

Lafarguette F, Leple JC, Dejardin A, Laurans F, Costa G, Lesage-Descauses MC, Pilate G. 2004. Poplar genes encoding fasciclin-like arabinogalactan proteins are highly expressed in tension wood. *The New Phytologist* 164:107–121.

Lewis DR, Ramirez MV, Miller ND, Vallabhaneni P, Ray WK, Helm RF, Winkel BS, Muday GK. 2011. Auxin and ethylene induce flavonol accumulation through distinct transcriptional networks. *Plant Physiology* 156:144–164.

Li E, Bhargava A, Qiang W, Friedmann MC, Forneris N, Savidge RA, Johnson LA, Mansfield SD, Ellis BE, Douglas CJ. 2012. The Class II KNOX gene KNAT7 negatively regulates secondary wall formation in *Arabidopsis* and is functionally conserved in *Populus*. *The New Phytologist* 194:102–115.

Lopez D, Tocquard K, Venisse JS, Legué V, Roeckel-Drevet P. 2014. Gravity sensing, a largely misunderstood trigger of plant orientated growth. *Frontiers in Plant Science* 5:610.

MacMillan CP, Taylor L, Bi Y, Southerton SG, Evans R, Spokevicius A. 2015. The fasciclin-like arabinogalactan protein family of *Eucalyptus grandis* contains members that impact wood biology and biomechanics. *The New Phytologist* 206:1314–1327.

Mansuri ML, Parihar P, Solanki I, Parihar MS. 2014. Flavonoids in modulation of cell survival signalling pathways. *Genes & Nutrition* 9:400.

Martin L, Leblanc-Fournier N, Azri W, Lenne C, Henry C, Coutand C, Julien JL. 2009. Characterization and expression analysis under bending and other abiotic factors of PtaZFP2, a poplar gene encoding a Cys2/His2 zinc finger protein. *Tree Physiology* 29:125–136.

Mauriat M, Moritz T. 2009. Analyses of GA20ox- and GID1-over-expressing aspen suggest that gibberellins play two distinct roles in wood formation. *The Plant Journal* 58:989–1003.

Mizrachi E, Maloney VJ, Silberbauer J, Hefer CA, Berger DK, Mansfield SD. Myburg AA. 2014. Investigating the molecular underpinnings underlying morphology and changes in carbon partitioning during tension wood formation in *Eucalyptus*. *The New Phytologist* 206:1351–1363.

Nick P. 2011. Mechanics of the cytoskeleton. In: Wojtaszek P, ed. *Signaling and communication in plants*. *mechanical integration of plant cells and plants*. Heidelberg: Springer, 53–90.

Nugroho WD, Nakaba S, Yamagishi Y, Begum S, Marsoem SN, Ko JH, Jin HO, Funada R. 2013. Gibberellin mediates the development of gelatinous fibres in the tension wood of inclined *Acacia mangium* seedlings. *Annals of Botany* 112:1321–1329.

Nugroho WD, Yamagishi Y, Nakaba S, Fukuhara S, Begum S, Marsoem SN, Ko JH, Jin HO, Funada R. 2012. Gibberellin is required for the formation of tension wood and stem gravitropism in *Acacia mangium* seedlings. *Annals of Botany* 110:887–895.

Paux E, Carocha V, Marques C, Mendes de Sousa A, Borralho N, Sivadon P, Grima-Pettenati J. 2005. Transcript profiling of *Eucalyptus* xylem genes during tension wood formation. *The New Phytologist* 167:89–100.

Pfaffl MW. 2001. A new mathematical model for relative quantification in real-time RT-PCR. *Nucleic Acids Research* 29:e45.

Pfaffl MW, Tichopad A, Prgomet C, Neuvians TP. 2004. Determination of stable housekeeping genes, differentially regulated target genes and sample integrity: BestKeeper–Excel-based tool using pair-wise correlations. *Biotechnology Letters* 26:509–515.

Rajinikanth M, Harding SA, Tsai CJ. 2007. The glycine decarboxylase complex multienzyme family in *Populus*. *Journal of Experimental Botany* 58:1761–1770.

Rose JK, Saladié M, Catalá C. 2004. The plot thickens: new perspectives of primary cell wall modification. *Current Opinion in Plant Biology* 7:296–301.

Roudier F, Fernandez AG, Fujita M, Himmelspach R, Borner GH, Schindelman G, Song S, Baskin TI, Dupree P, Wasteneys GO, Benfey PN. 2005. COBRA, an *Arabidopsis* extracellular glycosyl-phosphatidyl inositol-anchored protein, specifically controls highly anisotropic expansion through its involvement in cellulose microfibril orientation. *The Plant Cell* 17:1749–1763.

Ruelle J. 2014. Morphology, anatomy and ultrastructure of reaction wood. In: Gardiner B, Barnett J, Saranpää P, Gril J, eds. *The biology of reaction wood*. Berlin, Heidelberg: Springer-Verlag, 13–35.

Santelia D, Henrichs S, Vincenzetti V, Sauer M, Bigler L, Klein M, Bailly A, Lee Y, Friml J, Geisler M, Martinoia E. 2008. Flavonoids redirect PIN-mediated polar auxin fluxes during root gravitropic responses. *The Journal of Biological Chemistry* 283:31218–31226.

Schindelman G, Morikami A, Jung J, Baskin TI, Carpita NC, Derbyshire P, McCann MC, Benfey PN. 2001. COBRA encodes a putative GPI-anchored protein, which is polarly localized and necessary for oriented cell expansion in *Arabidopsis*. *Genes & Development* 15:1115–1127.

Shen B, Li C, Tarczynski MC. 2002. High free-methionine and decreased lignin content result from a mutation in the *Arabidopsis* S-adenosyl-L-methionine synthetase 3 gene. *The Plant Journal* 29:371–380.

Sjödin A, Street NR, Sandberg G, Gustafsson P, Jansson S. 2009. The *Populus* Genome Integrative Explorer (PopGenIE): a new resource for exploring the *Populus* genome. *The New Phytologist* 182:1013–1025.

Smyth GK. 2004. Linear models and empirical Bayes methods for assessing differential expression in microarray experiments. *Statistical Applications in Genetics and Molecular Biology* 3:3.
Song Y, Li G, Nowak J, Zhang X, Xu D, Yang X, Huang G, Liang W, Yang L, Wang C, Bulone V, Nikoloski Z, Hu J, Persson S, Zhang D. 2019. The rice actin-binding protein RMD regulates light-dependent shoot gravitropism. *Plant Physiology* **181**:630–644.

Storey JD. 2003. The positive false discovery rate: a Bayesian interpretation and the q-value. *Annals of Statistics* **31**:2013–2015.

Sundell D, Mannapperuma C, Netotea S, Delhomme N, Lin YC, Sjödin A, Van de Peer Y, Jansson S, Hvidsten TR, Street NR. 2015. The plant genome integrative explorer resource: PlantGenIE.org. *The New Phytologist* **208**:1149–1156.

Tejos R, Sauer M, Vanneste S, Palacios-Gomez M, Li H, Heilmann M, van Wijk R, Vermeer JE, Heilmann I, Munnik T, Friml J. 2014. Bipolar plasma membrane distribution of phosphoinositides and their requirement for auxin-mediated cell polarity and patterning in *Arabidopsis*. *The Plant Cell* **26**:2114–2128.

Tenhaken R. 2015. Cell wall remodeling under abiotic stress. *Frontiers in Plant Science* **5**:771.

Timell TE. 1986. *Compression wood in gymnosperms*. Berlin: Springer.

Tocquard K, Lopez D, Decourteix M, Thibaut B, Julien JL, Label P, Leblanc-Fournier N, Roeckel-Drevet P. 2014. The molecular mechanisms of reaction wood induction. In: Gardiner B, Barnett J, Saranpää P, Gril J, eds. *The biology of reaction wood*. Berlin, Heidelberg: Springer-Verlag, 107–138.

Vandenbussche F, Callebert P, Zadnikova P, Benkova E, Van Der Straeten D. 2013. Brassinosteroid control of shoot gravitropism interacts with ethylene and depends on auxin signaling components. *American Journal of Botany* **100**:215–225.

Vander Mijnsbrugge K, Meyermans H, Van Montagu M, Bauw G, Boerjan W. 2000. Wood formation in poplar: identification, characterization, and seasonal variation of xylem proteins. *Planta* **210**:589–598.

Vandesompele J, De Preter K, Pattyn F, Poppe B, Van Roy N, De Paepe A, Speleman F. 2002. Accurate normalization of real-time quantitative RT-PCR data by geometric averaging of multiple internal control genes. *Genome Biology* **3**:RESEARCH0034.

Villalobos DP, Díaz-Moreno SM, Said el-SS, Cañas RA, Osuna D, Van Kerckhoven SH, Bautista R, Claros MG, Cánovas FM, Cantón FR. 2012. Reprogramming of gene expression during compression wood formation in pine: coordinated modulation of S-adenosylmethionine, lignin and lignan related genes. *BMC Plant Biology* **12**:100.

Wang H, Jiang C, Wang C, Yang Y, Yang L, Gao X, Zhang H. 2015. Antisense expression of the fasciclin-like arabinogalactan protein FLA6 gene in *Populus* inhibits expression of its homologous genes and alters stem biomechanics and cell wall composition in transgenic trees. *Journal of Experimental Botany* **66**:1291–1302.

Yang YH, Dudoit S, Luu P, Lin DM, Peng V, Ngai J, Speed TP. 2002. Normalization for cDNA microarray data: a robust composite method addressing single and multiple slide systematic variation. *Nucleic Acids Research* **30**:e15.

Zhou D, Zhou X, Ling Y, Zhang Z, Su Z. 2010. AgriGO: a GO analysis toolkit for the agricultural community. *Nucleic Acids Research* **38**: doi:10.1093/nar/gkq310.

Zinkgraf M, Gerttula S, Zhao S, Filkov V, Groover A. 2018. Transcriptional and temporal response of *Populus* stems to gravi-stimulation. *Journal of Integrative Plant Biology* **60**:578–590.