Ethylene- and shade-induced hypocotyl elongation

0 downloads 0 Views 5MB Size Report
Jun 21, 2016 - ALW Ecogenomics 84410004 to L.A.C.J.V., ALW VENI. 25 ... 41 shade-induced genome-wide gene expression changes in ..... hypocotyl (27 h). ..... addition, under ethylene specifically, Photosystem II and I genes were mostly.
Plant Physiology Preview. Published on June 21, 2016, as DOI:10.1104/pp.16.00725

1

Short title: Ethylene- and shade-induced stress escape

2

Corresponding authors: Rashmi Sasidharan, Ronald Pierik

3

Address: Padualaan 8, 3584 CH, Utrecht, The Netherlands

4

Telephone number: +31 30 2536838

5

Email: [email protected], [email protected]

6 7

Ethylene- and shade-induced hypocotyl elongation share transcriptome

8

patterns and functional regulators

9

Das D.1, St. Onge K.R.1,2, Voesenek L.A.C.J.1, Pierik R.1* and Sasidharan R.1*

10

*shared senior and corresponding authors

11

1

12

Netherlands

13

2

Plant Ecophysiology, Institute of Environmental Biology, Utrecht University, 3584CH Utrecht, The Department of Biological Sciences, CW405 BioSci Building, University of Alberta, T6J2E9, Canada

14 15

Summary: Ethylene and shade share a conserved set of molecular regulators to

16

control plasticity of hypocotyl elongation in Arabidopsis.

17 18

AUTHOR CONTRIBUTIONS

19

R.P., R.S and L.A.C.J.V. conceived the original research plans and project. R.S.,

20

R.P, L.A.C.J.V., K.R.S. supervised the experiments. D.D. performed most of the

21

experiments. D.D., K.R.S., R.P. and R.S. designed the experiments. D.D. and

22

K.R.S. analyzed the data. D.D. wrote the article with contributions of all the authors.

23

FUNDING INFORMATION

24

This research was supported by the Netherlands Organisation for Scientific

25

Research (grant nos. ALW Ecogenomics 84410004 to L.A.C.J.V., ALW VENI

26

86312013 and ALW 82201007 to R.S., ALW VIDI 86412003 to R.P.) and a Utrecht

27

University Scholarship to D.D.

28 29

CORRESPONDING AUTHOR EMAIL: [email protected], [email protected]

30 31 32

1 Downloaded from www.plantphysiol.org on June 22, 2016 - Published by www.plantphysiol.org Copyright © 2016 American Society of Plant Biologists. All rights reserved.

Copyright 2016 by the American Society of Plant Biologists

33

ABSTRACT

34

Plants have evolved shoot elongation mechanisms to escape from diverse

35

environmental stresses such as flooding and vegetative shade. The apparent

36

similarity in growth responses suggests possible convergence of the signalling

37

pathways. Shoot elongation is mediated by passive ethylene accumulating to high

38

concentrations in flooded plant organs and by changes in light quality and quantity

39

under vegetation shade. Here we study hypocotyl elongation as a proxy for shoot

40

elongation and delineated Arabidopsis hypocotyl length kinetics in response to

41

ethylene and shade. Based on these kinetics, we further investigated ethylene and

42

shade-induced genome-wide gene expression changes in hypocotyls and cotyledons

43

separately. Both treatments induced a more extensive transcriptome reconfiguration

44

in the hypocotyls compared to the cotyledons. Bioinformatics analyses suggested

45

contrasting regulation of growth promotion- and photosynthesis-related genes.

46

These analyses also suggested an induction of auxin, brassinosteroid and gibberellin

47

signatures and the involvement of several candidate regulators in the elongating

48

hypocotyls. Pharmacological and mutant analyses confirmed the functional

49

involvement of several of these candidate genes and physiological control points in

50

regulating stress-escape responses to different environmental stimuli. We discuss

51

how these signaling networks might be integrated and conclude that plants, when

52

facing different stresses, utilise a conserved set of transcriptionally regulated genes

53

to modulate and fine tune growth.

54 55

INTRODUCTION

56

All organisms, including plants, assess and respond to both biotic and abiotic factors

57

in their environments (Pierik and De Wit, 2014; Pierik and Testerink, 2014; Franklin

58

et al., 2011; Quint et al., 2016; Osakabe et al., 2014; Voesenek and Bailey-Serres,

59

2015). However, unlike animals, plants, cannot move away from extremes in their

60

surrounding environment but rather rely on various plastic morphological and

61

metabolic responses. Such response traits include changes in plant architecture to

62

escape the stress and optimize resource capture (Pierik and Testerink, 2014;

63

Mickelbart et al., 2015). With energy reserves being invested in escape traits, plants

64

often have lower plant biomass and crop yield (Casal, 2013). Molecular investigation

2 Downloaded from www.plantphysiol.org on June 22, 2016 - Published by www.plantphysiol.org Copyright © 2016 American Society of Plant Biologists. All rights reserved.

65

of the different signalling pathways controlling these traits along with the

66

characterization of underlying molecular components would not only enhance

67

fundamental knowledge of stress-induced plasticity but also benefit crop

68

improvement.

69

Plants are highly sensitive to changes in their light environment. Young plants

70

growing in a canopy experience changes in light quality and quantity due to

71

neighbouring plants and compete to harvest optimum light (Pierik and De Wit, 2014;

72

Casal et al., 2013). When a plant cannot outgrow its neighbours, it experiences

73

complete vegetation shade (hereafter termed as “shade”) which, in addition to low

74

red: far-red (R:FR), is marked by a significant decline in blue light and overall light

75

quantity. These changes initiate so-called Shade Avoidance Syndrome (SAS)

76

responses consisting of petiole, hypocotyl and stem elongation; reduction of

77

cotyledon and leaf expansion; upward movement of leaves (hyponasty), decreased

78

branching and increased apical dominance (Vandenbussche et al., 2005; Franklin,

79

2008; Casal, 2012; Pierik and De Wit, 2014). Shade-induced elongation comprises a

80

complex network of photoreceptor-regulated transcriptional and protein-level

81

regulation involving BASIC HELIX LOOP HELIX (BHLH) and HOMEODOMAIN-

82

LEUCINE ZIPPER, (HD-ZIP) transcription factors and auxin, gibberellin and

83

brassinosteroid hormone genes (Casal, 2012; 2013). Flooding often leads to partial

84

or complete submergence of plants. Water severely restricts gas diffusion and the

85

consequent limited exchange of O2 and CO2 restricts respiration and photosynthesis.

86

Another consequence is the rapid accumulation of the volatile hormone ethylene.

87

Ethylene is considered an important regulator of adaptive responses to flooding,

88

including accelerated shoot elongation responses that bring leaf tips from the water

89

layer into the air (Sasidharan and Voesenek, 2015; Voesenek and Bailey-Serres,

90

2015). In deepwater rice, this flooding-induced elongation response involves

91

ethylene-mediated induction of members of the group VII Ethylene Response

92

Factors (ERF-VII) family, a decline in active abscisic acid (ABA) and consequent

93

increase in gibberellic acid (GA) responsiveness and promotion of GA biosynthesis

94

(Hattori et al., 2009). In submerged Rumex palustris petioles, ethylene also rapidly

95

stimulates cell wall acidification and transcriptional induction of cell wall modification

96

proteins to facilitate rapid elongation (Voesenek and Bailey-Serres, 2015). Shade

97

cues are reported to enhance ethylene production resulting in shade avoidance 3 Downloaded from www.plantphysiol.org on June 22, 2016 - Published by www.plantphysiol.org Copyright © 2016 American Society of Plant Biologists. All rights reserved.

98

phenotypes (Pierik et al., 2004). However, these responses are mediated by

99

ethylene concentrations of a much lower magnitude than that occuring in flooded

100

plant organs (1µl L-1) (Sasidharan and Voesenek, 2015).

101

So far, it is largely unknown to what extent these growth responses to such highly

102

diverse environmental stimuli, share physiological and molecular components

103

through time. A preliminary study in R. palustris showed that GA is a common

104

regulator of responses to both submergence and shade (Pierik et al. 2005). Although

105

submergence is a compound stress, rapid ethylene accumulation is considered an

106

early and reliable flooding signal triggering plant adaptive responses. High ethylene

107

concentrations as occur within submerged plant organs, promote rapid shoot

108

elongation (Voesenek and Bailey-Serres, 2015). This submergence response, which

109

has been extensively characterised in rice and Rumex (Hattori et al. 2009; van Veen

110

et al., 2013), can be almost completely mimicked by the application of saturating (1µl

111

L-1) ethylene concentrations (Sasidharan and Voesenek, 2015). Saturating ethylene

112

concentrations were therefore used here as a submergence mimic. Shade was given

113

as true shade which combines the three known key signals that trigger elongation

114

(red and blue light depletion with relative far-red enrichment).

115

A hypocotyl elongation assay in Arabidopsis thaliana (Col-0) was used as a proxy for

116

shoot elongation under ethylene and shade in order to study to what extent ethylene

117

and shade responses share molecular signalling components. Although ethylene

118

suppresses Arabidopsis hypocotyl elongation in dark, high ethylene concentrations

119

in light (as occurs during submergence) stimulate hypocotyl elongation in

120

Arabidopsis (Smalle et al., 1997; Zhong et al., 2012). Also upon simulated shade,

121

Arabidopsis demonstrates pronounced hypocotyl elongation (Morelli and Ruberti,

122

2000). To capture early physiological responses and gene expression changes in

123

response to ethylene and shade, Arabidopsis seedling hypocotyl elongation and

124

cotyledon expansion were examined over time. The two treatments elicited

125

characteristic hypocotyl growth kinetics. To uncover the transcriptomic changes

126

regulating the elongation response to these signals, an organ-specific genome-wide

127

investigation was carried out on hypocotyls and cotyledons separately at three time-

128

points corresponding to distinct hypocotyl elongation phases. Clustering analyses in

129

combination with biological-enrichment tests allowed identification of gene clusters

130

with expression patterns matching the hypocotyl growth trends across the three time4 Downloaded from www.plantphysiol.org on June 22, 2016 - Published by www.plantphysiol.org Copyright © 2016 American Society of Plant Biologists. All rights reserved.

131

points in both the treatments. Correlation of genome-wide hypocotyl and cotyledon-

132

specific transcriptomic changes to publicly available microarray data on hormone

133

treatments identified enriched hormonal signatures of auxin, brassinosteroid and

134

gibberellin in hypocotyl tissues and several potential growth regulatory candidate

135

genes. Using hormone mutants and chemical inhibitors, we confirmed the combined

136

involvement of these hormones and candidate regulators in the hypocotyl elongation

137

response to ethylene and shade. We suggest that growth responses to diverse

138

environmental stimuli like ethylene and shade converge on a common regulatory

139

module consisting of both positive and negative regulatory proteins that interact with

140

a hormonal triad to achieve a controlled fine-tuned growth response.

141 142 143

RESULTS

144

Delineation of hypocotyl elongation kinetics under ethylene and shade in

145

Arabidopsis seedlings

146

Exogenous application of ethylene (1 µl L-1) in light-grown seedlings resulted in thick,

147

yet elongated hypocotyls and smaller cotyledons as compared to untreated controls.

148

Shade, achieved by the use of a green filter, stimulated strong hypocotyl elongation

149

in seedlings but resulted in mildly smaller cotyledons compared to controls (Fig. 1A).

150

Hypocotyl length increments in the two treatments relative to control were around 2-

151

fold under ethylene and greater than 3-fold under shade (Figs. 1, A and B). For both

152

treatments growth stimulation was strongest in the first two days of treatment (Fig.

153

1C) and faded out in the subsequent days. When combined, shade and ethylene

154

exposure resulted in hypocotyl lengths that were intermediate to the individual

155

treatments (Supplemental Fig. S1).

156

To get a more detailed time-line of the early elongation kinetics, we performed a

157

follow-up experiment with 3 h measurement intervals for the first 33 h (Fig. 1D).

158

Ethylene-mediated stimulation of hypocotyl elongation started only in the middle of

159

the dark period after 15 h. However, under shade, longer hypocotyls were recorded

160

already after 3 h and this rapid stimulation continued until the start of the dark period.

161

Interestingly, accelerated elongation was again observed at around 24 h after the

162

start of the treatments (when the lights were switched on). Based on this time-line,

163

we determined epidermal cell lengths at time-points 0, 3, 7.5, 15 and 27 h (Gendreau 5 Downloaded from www.plantphysiol.org on June 22, 2016 - Published by www.plantphysiol.org Copyright © 2016 American Society of Plant Biologists. All rights reserved.

164

et al., 1997; Fig. 1E), the latter four of which are either prior to start of accelerated

165

growth or during it in response to ethylene and shade. The rapid stimulation of

166

elongation under shade starts at the base of the hypocotyl (3 h) and then progresses

167

all along the hypocotyl with maximum elongation occurring in the middle segment

168

while under ethylene accelerated elongation is observed at the middle-bottom of

169

hypocotyl (27 h).

170 171

Organ-specific transcriptomics in hypocotyl and cotyledon under ethylene and

172

shade

173

The transcriptome response to ethylene and shade in hypocotyl and cotyledon

174

tissues was characterised using Affymetrix Arabidopsis Gene 1.1 ST arrays at three

175

time points of hypocotyl length kinetics (1.5 h, 13.5 h and 25.5 h) (Fig. 1D). Principal

176

component analysis (PCA) (Abdi and Williams, 2010) of all replicate samples for

177

hypocotyl and cotyledon exposed to control, ethylene or shade conditions showed

178

that replicate samples generally clustered together (Fig. 2A). The first principal

179

component (34.2%) separates tissue-specific samples, whereas the second principal

180

component (13.0%) showed separate clustering of the 13.5h samples which falls

181

during the dark period.

182

Hierarchical clustering (HC) (Eisen et al., 1998) of mean absolute expression

183

intensities for the different main samples (combination of 3 replicates) revealed

184

similar trends (Fig. 2B). Fig. 2C shows the distribution of up- and down-regulated

185

differentially expressed genes (DEGs) (genes with adjusted p-value ≤ 0.01) in

186

hypocotyl and cotyledon for ethylene and shade at the three harvest time points 1.5

187

h, 13.5 h, 25.5 h respectively. In both the conditions and tissues, the number of both

188

up- and down-regulated DEGs increased with time. Ethylene regulated substantially

189

more DEGs in the hypocotyl as compared to shade at all harvest time points (Fig.

190

2C). Data analysis identified 6668 and 4741 genes (hereafter termed as “total

191

DEGs”) that were differentially expressed in hypocotyl at one or more of the three

192

tissue harvest time points by ethylene and shade respectively. Interestingly, in the

193

cotyledon, at 1.5 h, the number of significant DEGs under ethylene was higher than

194

under shade but at the subsequent two time points, shade regulated more genes. In

195

the cotyledon, 1197 and 2173 DEGs were identified that were differentially

196

expressed at one or more of the three tissue harvest time points by ethylene and

197

shade respectively. 6 Downloaded from www.plantphysiol.org on June 22, 2016 - Published by www.plantphysiol.org Copyright © 2016 American Society of Plant Biologists. All rights reserved.

198

Interestingly, at 1.5 h, there was more transcriptional regulation in ethylene- than

199

shade-exposed hypocotyls even though subsequent hypocotyl elongation was much

200

more rapid in shade (Figs. 2C and 1D). For ethylene-specific downregulated DEGs

201

at 1.5 h, the topmost enriched GO term was cell wall organization (containing 30

202

genes), which suggested a repression of growth-promoting genes and possible lack

203

of ethylene-mediated elongation at 1.5 h (Supplemental Fig. S2). We also found 6

204

genes (AT1G65310, XYLOGLUCAN ENDOTRANSGLUCOSYLASE / HYDROLASE

205

17 (XTH17); AT5G23870, pectin acetyl esterase; AT3G06770, glycoside hydrolase;

206

AT5G46240, POTASSIUM CHANNEL IN ARABIDOPSIS THALIANA 1 (KAT1);

207

AT1G29460, SMALL AUXIN UPREGULATED RNA 65 (SAUR65) and AT3G02170,

208

LONGIFOLIA 2 (LNG2)) in the shade-specific upregulated 32 genes at 1.5 h, with

209

implied cell expansion

210

elongation response (Sasidharan et al., 2010; Philippar et al., 2004; Nozue et al.,

211

2015; Chae et al., 2012; Lee et al., 2006).

roles and could possibly be associated with the rapid

212 213

Different gene expression clusters contributing to hypocotyl growth in

214

ethylene and shade

215

In order to find specific genes regulating the elongation phenotype under both

216

treatments, we used temporal clustering of DEGs based on expression values. Due

217

to distinct hypocotyl length kinetics in response to ethylene and shade (Fig. 1D), we

218

searched for a set of temporally co-expressed genes that could potentially contribute

219

to this treatment-specific kinetics. Time-point based clustering was performed for the

220

6668 ethylene and 4741 shade total DEGs based on the positive or negative

221

magnitude of log2FC for DEGs at the three time points (Fig. 3, A and D). The gene

222

expression patterns in clusters-1 and 5 across the three time points matched the

223

ethylene hypocotyl growth kinetics closely (Fig. 3, B and C). Similarly, gene

224

expression kinetics in clusters-1 and 3 matched the hypocotyl length kinetics in

225

shade (Fig. 3, E and F). These growth pattern matching clusters were termed

226

“positive”. All the clusters with mirror image of gene expression profiles to that of the

227

positive clusters (clusters-8 and 4 in ethylene and clusters-8 and 6 in shade) were

228

termed as “negative” clusters.

7 Downloaded from www.plantphysiol.org on June 22, 2016 - Published by www.plantphysiol.org Copyright © 2016 American Society of Plant Biologists. All rights reserved.

229

Next, a hypergeometric over-representation test for selected MapMan bins (“stress”,

230

“hormone”, “signalling”, “RNA.Regulation of transcription” and “cell wall”) was carried

231

out for the temporal gene clusters (Fig. 3G). Interestingly, “cell wall”, “hormone” and

232

“signalling” were highly co-enriched in positive clusters (cluster-1 and 5 for ethylene

233

and cluster-1 for shade) which hints at co-regulation of genes mapped to these terms

234

during transcriptomic response to ethylene and shade in the hypocotyl.

235

To identify growth promotion related DEGs, we identified DEGs common to both the

236

treatments in the clusters previously designated as ‘positive’ and ‘negative’ clusters

237

(Fig. 3, H and I). 997 DEGs were obtained from a venn diagram between the

238

treatment-specific positive clusters, upregulated at atleast two time-points, and

239

hereafter called “Common Up”. Similarly, 824 DEGs were shared between ethylene

240

and shade negative clusters, were downregulated at atleast two time points and

241

hereafter called “Common Down”.

242

Enriched functional categories in the different gene sets from Venn diagrams of

243

positive and negative clusters, were identified using the GeneCodis tool (Tabas-

244

Madrid et al., 2012) (Fig. 4). In the Common Up set, we found a variety of growth

245

associated GO categories including cell wall modification, hormone (auxin and

246

brassinosteroid) signaling and metabolism, transport processes, tropisms, response

247

to abiotic stimuli and signal transduction. The ethylene-specific set for positive

248

clusters was enriched for ethylene-associated terms as expected, but also for

249

various sugar metabolic, endoplasmic reticulum (ER)-related and protein post-

250

translational modification-related processes. Some of the enriched GO terms in the

251

ethylene-specific set for positive clusters were also found in the Common Up set, but

252

caused by different genes in the same GO category, including those associated with

253

growth, hormones and transport processes. The shade-specific set for positive

254

clusters showed only few clear GO enrichments such as trehalose metabolism,

255

secondary cell wall biogenesis and amino acid metabolism, but also shared some

256

with the Common Up set like shade avoidance and protein phosphorylation and with

257

both Common Up and ethylene-specific set for positive clusters, such as response to

258

auxin and unidimensional growth.

259

In the Common Down set, GO terms associated with photosynthesis, primary and

260

secondary metabolism, response to biotic and abiotic stress, as well as

8 Downloaded from www.plantphysiol.org on June 22, 2016 - Published by www.plantphysiol.org Copyright © 2016 American Society of Plant Biologists. All rights reserved.

261

photomorphogenesis were enriched. The latter is striking given that ethylene does

262

not alter the light environment. The ethylene-specific set for negative clusters

263

included strong enrichment of circadian rhythm and a variety of photosynthesis-and

264

chloroplast-associated GO terms, partially shared with the Common Down set. In the

265

shade-specific set for negative clusters, flavonoid/anthocyanin biosynthesis and

266

response to UV-B and heat were enriched. The shade-specific set for negative

267

clusters shared terms from defence-associated GO categories and cadmium, karrikin

268

response with Common Down set. Terms common to all three sets of negative

269

clusters were related to photomorphogenesis and metabolic process.

270 271

Functional characterization: shared components in ethylene -and shade-

272

mediated regulation of hypocotyl length

273

We classified Common Up and Down set genes into transcriptional regulators,

274

hormone metabolism genes, signalling genes and cell wall genes and also applied a

275

logFC filter (see “Materials and Methods” sub-section “LogFC Filter and Gene

276

Classification”) to obtain a final list of 53 and 8 genes in the two sets respectively

277

(Fig. 5, A and B). We selected a subset of candidate regulators for functional testing

278

and made sure to include two transcription factors since these may be relatively

279

upstream in the convergence of signaling pathways. Obvious targets for these

280

regulators would be different plant hormones, and before zooming in on these, we

281

ran a hormonometer analysis with our transcriptome data to further subset our

282

candidate regulators.

283 284

Transcription factor candidates:

285

ARABIDOPSIS THALIANA HOMEOBOX PROTEIN 28 (ATHB28) is a zinc-finger

286

homeodomain (ZF-HD) transcription factor that showed upto 2-fold induction in the

287

hypocotyls in ethylene and shade respectively. A homozygous null mutant

288

(Supplemental Fig. S3, A-D), “athb28”, showed a significant reduction in hypocotyl

289

lengths under both ethylene and shade (Fig. 6A), consistent with its induction upon

290

both the treatments.

9 Downloaded from www.plantphysiol.org on June 22, 2016 - Published by www.plantphysiol.org Copyright © 2016 American Society of Plant Biologists. All rights reserved.

291

Both shade and ethylene also significantly induced transcript levels of the bHLH

292

transcription factor INCREASED LEAF INCLINATION1 BINDING bHLH1 like-1

293

(IBL1) Arabidopsis hypocotyls. However, the ibl1 (SALK T-DNA insertion line)

294

showed wild-type responses to the treatments (Fig. 6B). The bHLH transcription

295

factors IBL1 and its homolog IBH1 have been implicated in repressing BR-mediated

296

cellular elongation. In an ibl1 mutant, IBH1 is still present and may negatively

297

regulate cell elongation independently of IBL1. The 35S overexpression line, IBL1OE

298

(Zhiponova et al., 2014), had shorter hypocotyls than the wild-type under control

299

conditions. IBL1OE also lacked ethylene and shade-induced hypocotyl elongation

300

implying an inhibitory role for IBL1 (Fig. 6C).

301 302

Hormone candidates: auxin, brassinosteroid and gibberellin

303

To further investigate the significant hormone related changes amongst the growth

304

related DEGs we analysed our data using Hormonometer (Volodarsky et al., 2009).

305

For both treatments, the hormonal signatures across the three time-points for BR

306

and GA most closely matched the hypocotyl elongation kinetics (Figs. 1D and 7A).

307

The analysis also showed significant correlations with auxin responses for all data

308

sets.

309

In the Common Up set 49 genes were present that were all also auxin-regulated,

310

whereas there were 14 that were also BR-regulated. In the Common Down set 16

311

genes were present that can be regulated by auxin, 16 that can also be ABA-

312

regulated and 13 genes that are also JA-regulated. Interestingly, there were no

313

genes for BR in the Common Down set. In addition, genes involved in auxin-

314

conjugation genes (GRETCHEN HAGEN 3 FAMILY PROTEIN 3.17 (GH3.17) and

315

AT5G13370), GA catabolising genes (GIBBERELLIN 2 OXIDASE (GA2OX) 2,

316

GA2OX4 and GA2OX7) and JA augmenting genes (LIPOXYGENASE (LOX) 1,

317

LOX2 and LOX3) were down-regulated.

318

In order to test the possible role of auxin, GA and BR, in mediating shade and

319

ethylene-induced hypocotyl elongation, we first tested the effects of pharmacological

320

inhibitors of these hormones on shade and ethylene-induced hypocotyl elongation.

321

To visualize the auxin effect we treated the pIAA19:GUS auxin response marker line

322

with the auxin transport inhibitor, 1-N-Naphthylphthalamic acid (NPA). As shown in 10 Downloaded from www.plantphysiol.org on June 22, 2016 - Published by www.plantphysiol.org Copyright © 2016 American Society of Plant Biologists. All rights reserved.

323

Fig. 7B, application of NPA (25 µM) inhibited hypocotyl elongation and also strongly

324

reduced staining in the hypocotyl region. This inhibition was rescued by IAA (10 µM).

325

The auxin perception inhibitor, α-(phenylethyl-2-one)-indole-3-acetic acid, PEO-IAA

326

(100 µM), strongly reduced staining in the whole pIAA19:GUS seedling and inhibited

327

elongation in response to both treatments. The significant inhibition of ethylene and

328

shade-induced hypocotyl elongation by the different auxin inhibitors is quantitated in

329

Fig. 7,C-E. The addition of NPA, yucasin (auxin biosynthesis inhibitor) and α-

330

(phenylethyl-2-one)-indole-3-acetic acid (PEO-IAA; auxin antagonist) inhibited

331

hypocotyl elongation under both the treatments confirming that all three aspects of

332

auxin are required for ethylene- and shade-induced hypocotyl elongation. BR

333

biosynthesis inhibitor brassinazole (BRZ) and GA biosynthesis inhibitor paclobutrazol

334

(PBZ) fully inhibited these elongation responses as well (Fig. 7, F and G).

335

To further validate the involvement of auxin, BR and GA in ethylene and shade-

336

induced hypocotyl elongation we tested hypocotyl elongation responses in a variety

337

of hormone mutants, including mutants for candidate genes from Figure 5.

338

Both the auxin receptor (tir1-1) and biosynthesis (wei8-1) mutants showed

339

significantly impaired hypocotyl elongation responses compared to the wild type

340

ethylene and shade response (Fig. 8, A and B). A similar effect was seen in the

341

auxin transport mutant pin3pin4pin7 which had severely reduced hypocotyl

342

elongation in both treatments (Fig. 8C).

343

The GA biosynthesis (ga1-3) and insensitive (gai) mutant both showed a complete

344

lack of hypocotyl elongation in both treatments (Fig. 8, D and E). We also tested the

345

GA biosynthesis mutant ga20ox1-3, since it was identified as a ‘common Up’ gene,

346

induced in response to both treatments (Fig. 5A). The ga20ox1-3 mutant showed a

347

significantly reduced elongation phenotype in both treatments compared to the wild

348

type response (Fig. 8F).

349

The BR receptor (bri1-116) and biosynthesis (dwf4-1) mutants both showed severe

350

hypocotyl elongation phenotypes, and did not respond to either treatment (Fig. 8, G

351

and I). In another biosynthesis mutant, rot3-1, while the ethylene elongation

352

response was absent there was a severely reduced shade response (Fig. 8H). Two

353

BR-metabolism related genes were identified in the common Up set (Fig. 5A):

354

BR6OX1 and BAS1. The bas1-2 mutant showed constitutive elongation in all

355

treatments ((Fig. 8J) confirming a negative role of BR catabolism through BAS1 in 11 Downloaded from www.plantphysiol.org on June 22, 2016 - Published by www.plantphysiol.org Copyright © 2016 American Society of Plant Biologists. All rights reserved.

356

hypocotyl elongation control. Although the BR biosynthetic mutant br6ox1 (cyp85a1-

357

2) mutant did not show any phenotypic alteration (Fig. 8K), a double mutant of

358

BR6OX1 and BR6OX2 (cyp85a1cyp85a2) showed a complete lack of elongation in

359

response to both ethylene and shade (Fig. 8L).

360 361 362

Discussion

363

Accelerated shoot elongation is a common mode of stress escape that allows plants

364

to grow away from stressful conditions (Pierik and Testerink, 2014). Stress escape

365

does, however, come at an energetic cost and is only beneficial if improved

366

conditions are achieved. Here our goal was to establish to what extent shade and

367

ethylene elicit similar responses through shared or distinct molecular pathways. In

368

our study we found distinct elongation kinetics in ethylene and shade for Arabidopsis

369

hypocotyls, differing in both temporal regulation and the degree of response. Shade

370

treatment evoked a rapid, strong and persisting hypocotyl elongation, whereas

371

ethylene initially inhibited elongation and only in the first night period started to

372

promote hypocotyl length (Fig. 1). In both treatments, hypocotyl growth involved

373

enhanced epidermal cell elongation. Previous studies have shown that low R:FR

374

induces ethylene biosynthesis in Arabidopsis (Pierik et al., 2009) and it could be

375

argued that the shade response might act through ethylene. However, ethylene

376

marker genes were not induced in shade, and the ethylene-insensitive ein3eil1

377

mutant retained a full response to shade (Supplemental Fig. S4), ruling out a role for

378

ethylene in the shade response. Interestingly, combining shade and ethylene

379

treatments did not lead to an additive response and instead dampened shade-

380

induced hypocotyl elongation (Supplemental Fig. S1). The growth inhibitory effect of

381

ethylene under shaded conditions could function similar to its effects in limiting

382

hypocotyl elongation in the dark i.e. via induction of negative growth regulators such

383

as ERF1 (Zhong et al., 2012). Transcriptome characterisation of the elongating

384

hypocotyl upon exposure to single shade and ethylene stresses indicated

385

considerable overlap between the two treatments. Thus, a large portion of DEGs

386

under both treatments may contribute to similar processes implying that they target

387

shared genetic components but have treatment-specific upstream regulatory factors.

388

12 Downloaded from www.plantphysiol.org on June 22, 2016 - Published by www.plantphysiol.org Copyright © 2016 American Society of Plant Biologists. All rights reserved.

389

Hypocotyl

growth

390

concurrently under ethylene and shade

391

By identifying gene clusters with expression patterns closely matching the distinct

392

ethylene and shade growth kinetics, we identified positive and negative clusters for

393

the respective treatments. These clusters were also among the bigger clusters that

394

contributed to most of the transcriptomic changes, suggesting that a large part of

395

transcriptomic response are associated with the hypocotyl growth and concurrent

396

biological processes. Functional enrichment analysis for the Common Up set (shared

397

between positive clusters of ethylene and shade) (Fig. 4), suggested an involvement

398

of growth promoting genes. Cell wall genes are all involved in mediating cellular

399

expansion in growing hypocotyls. However, they need to be controlled by either the

400

environmental signal directly or by upstream factors in the signal transduction

401

pathway. The Common Down set (shared between negative clusters of ethylene and

402

shade) (Fig. 4), was highly enriched in photosynthesis-related terms and proteins.

403

The effects of ethylene on photosynthesis can be positive or negative depending on

404

the context (Iqbal et al., 2012; Tholen et al., 2007). Low R:FR treated stems of

405

tomato showed reduced expression of photosynthetic genes (Cagnola et al., 2012).

406

This reduction was mainly due to a decrease in expression of Calvin cycle genes,

407

which we also observed for our Common Down set (Supplemental Table S1). In

408

addition, under ethylene specifically, Photosystem II and I genes were mostly

409

repressed. Thus acceleration of hypocotyl elongation is accompanied by repression

410

of

411

photosynthesis. This was also shown to be true for low R:FR treated elongating

412

stems of tomato (Cagnola et al., 2012). Light capture and carbon fixation are

413

minimized and energy is apparently invested in stimulating growth (Sulpice et al.,

414

2014; Henriques et al., 2014; Lilley et al., 2012). It would be interesting to investigate

415

how photomorphogenic responses are associated with and influence photosynthesis

416

and growth promotion.

genes associated

promotion

with

and

photosynthesis

non-elongation

processes

repression

like

occur

metabolism

and

417 418

Convergence of signaling pathways in response to ethylene and shade in

419

control of hypocotyl elongation

420

In shade avoidance responses, photoreceptors like phyB and cry1/2 would regulate

421

the elongation phenotype via control of Phytochrome Interacting Factor (PIF) levels.

422

However, since these are photoreceptors, it seems unlikely that these proteins 13 Downloaded from www.plantphysiol.org on June 22, 2016 - Published by www.plantphysiol.org Copyright © 2016 American Society of Plant Biologists. All rights reserved.

423

themselves would integrate information from the ethylene pathway as well (Park et

424

al., 2012; Li et al., 2012). It was shown by van Veen et al. (2013), that in the

425

submerged petioles of R. palustris (which displays petiole elongation under complete

426

submergence), early molecular components of light signalling (KIDARI, COP1, PIFs,

427

HD-ZIP IIs) are induced by ethylene independently of any change in light quality.

428

Over expression of PIF5 on the other hand leads to increased ethylene production in

429

etiolated Arabidopsis seedlings, causing inhibition of hypocotyl length (Khanna et al.,

430

2007). Interestingly, the downstream ethylene signal transduction protein EIN3 was

431

shown to physically interact with PIF3 (Zhong et al., 2012). We, therefore, suggest

432

that ethylene and shade might both induce this shared gene pool by, for example,

433

targeting (different) members of the PIF family of transcription factors. Since different

434

PIFs likely regulate the expression of at least partly shared target genes (Leivar and

435

Monte, 2014), this would explain our observed partial overlap in the transcriptional

436

response to shade and ethylene. PIFs are also known to directly bind and regulate

437

expression of other transcription factors like homeodomain (HD) TFs (Capella et al.,

438

2015; Kunihiro et al. 2011) in the control of shade avoidance responses. Indeed, our

439

TAIR motif analysis hinted towards the presence of significantly enriched binding

440

signatures of PIF/MYC proteins (CACATG) as well as HD proteins (TAATTA) in the

441

upstream promoter sequences of the Common Up set genes (Pfeiffer et al., 2014;

442

Kazan and Manners, 2013; Supplemental Fig. S5).

443

Several potential transcriptional regulators were identified in the narrowed down

444

Common Up gene set. A growth-promoting role of KIDARI in regulating elongation in

445

response to shade and ethylene was suggested previously (Hyun and Lee 2006; van

446

Veen et al., 2013). Upregulation of another bHLH encoding gene and a negative

447

regulator of elongation, IBL1 was observed for both treatments (Fig. 6C)). While

448

PIF4 induces IBH1 and IBL1, IBH1 represses PIF4 targets (Zhiponova et al., 2014).

449

IBH1 and its homolog IBL1 collectively regulate expression of a large number of BR-,

450

GA- and PIF4-regulated genes and this might their mode of action in shade- and

451

ethylene-induced hypocotyl elongation. In addition to these bHLH proteins, we show

452

that the ZF-HD transcription factor ATHB28 is also involved in regulating hypocotyl

453

elongation under ethylene and shade (Fig. 6A). Hong et al. (2011) showed that

454

another ZF-HD protein, MIF1 interacts strongly with four other ZF-HD proteins

455

including ATHB33 and ATHB28. This leads to non-functional MIF1-ATHB

456

heterodimers and inhibition of e.g. ATHB33-regulated expression and growth 14 Downloaded from www.plantphysiol.org on June 22, 2016 - Published by www.plantphysiol.org Copyright © 2016 American Society of Plant Biologists. All rights reserved.

457

promotion.

Transcriptomics

data

for

458

phenotype), shows downregulation of auxin, BR and GA responsive genes and

459

upregulation of ABA genes (Hu and Ma, 2006). We can speculate that MIF1 on one

460

hand and ATHB33 and ATHB28 on the other hand, might target the same set of

461

hormone genes, but in an opposite manner, to control growth. What remains to be

462

studied is how ethylene and shade regulate ZF-HD transcription factors and this will

463

be an important topic for future studies.

464

Well-established targets for above-mentioned PIFs and HD TFs are various aspects

465

of auxin signalling and homeostasis, such as YUCCA biosynthetic enzymes and

466

AUX/IAA proteins for signalling (Kunihiro et al. 2011; Li et al., 2012; Sun et al., 2012;

467

De Smet et al., 2013). Our list of candidate genes (Fig. 5) contained auxin-

468

responsive transcriptional regulator IAA3 and many of the auxin-responsive SAUR

469

genes which have been shown to positively modulate hypocotyl elongation (Kim et

470

al., 1998; Sun et al., 2012; Chae et al., 2012; Spartz et al., 2012) and may act

471

individually or in concert to regulate the phenotype. With reference to elongation

472

responses under shade in Arabidopsis seedlings, auxin seems to play a major role.

473

An increase in free auxin levels and its transport towards epidermal cells in hypocotyl

474

is necessary for low R:FR-mediated hypocotyl elongation (Tao et al., 2008;

475

Keuskamp et al., 2010; Zheng et al., 2016). The importance of YUCCAs and TAA1 in

476

low R:FR responses has been previously demonstrated (Li et al., 2012). It is

477

generally assumed that auxin synthesized in the cotyledons is required to regulate

478

hypocotyl elongation in response to e.g. low R:FR light conditions (Procko et al.

479

2014). Indeed, cotyledons are key regulators of hypocotyl elongation in a

480

phytochrome-dependent way (Endo et al., 2005; Warnasooriya and Montgomery,

481

2009; Estelle, 1998; Tanaka et al., 2002). In our data, hormonometer analysis

482

identified strong induction of auxin-associated genes in the cotyledons in both

483

treatments (Fig. 7). We speculate that the physiological regulation of hypocotyl

484

elongation in our study depends on cotyledons via auxin dynamics. However, we

485

also show that auxin is certainly not the only shared physiological regulator between

486

the ethylene and shade response.

487

GA20OX1 and BR6OX1 expression were up-regulated in patterns that closely

488

matched the hypocotyl elongation profiles (Fig. 5A), and hormonometer analysis also

489

revealed enrichment of GA and BR hormonal signatures in ethylene and shade

35S:MIF1

(displaying

short

hypocotyl

15 Downloaded from www.plantphysiol.org on June 22, 2016 - Published by www.plantphysiol.org Copyright © 2016 American Society of Plant Biologists. All rights reserved.

490

exposed hypocotyls. (Fig. 7). The positive role of GA in flooding-associated shoot

491

elongation (Voesenek and Bailey-Serres, 2015) and shade avoidance (Djakovic-

492

Petrovic et al., 2007) is well established. It is well known that GA20OX1 specifically

493

affects plant height without having any other major phenotypic effects (Barboza et

494

al., 2013; Rieu et al., 2008), and it has been shown to be involved in shade

495

avoidance (Nozue et al. 2015; Hisamatsu et al. 2005). In our data, ga20ox1 knockout

496

showed reduced elongation to shade as well as ethylene, extending its function from

497

controlling shade avoidance to ethylene mediated elongation responses (Fig. 8F).

498

GA biosynthesis via GA20OX1 can be induced by brassinosteroids, suggesting a

499

possible cross-talk between the two growth-promoting hormones (Unterholzner et

500

al., 2015). Although future studies are needed to establish if this crosstalk occurs

501

under the conditions tested here, we do confirm that BR is an important hormone

502

involved in both the responses since several BR mutants showed disturbed

503

elongation responses to ethylene and shade (Fig. 8G-8L). Interestingly, also auxin

504

and brassinosteroids have partially overlapping roles in hypocotyl elongation control

505

(Chapman et al., 2012; Nemhauser et al., 2004), further extending the crosstalk

506

towards a tripartite network. Among the BR mutants tested is the cyp85a1cyp85a2

507

mutant, encoding a double mutant for BR6OX2 and BR6OX1 which showed a

508

complete lack of elongation to both treatments (Fig. 8L) similar to a BRZ treatment

509

(Figs. 7F). BR6OX1 was one of the direct candidate genes identified from the

510

transcriptomics analysis (Figure 5A). A tripartite bHLH transcription factor module

511

consisting of IBH1, PRE and HBI1, has been previously implicated in regulating cell

512

elongation in response to hormonal and environmental signals (Bai et al., 2012).

513

Several BR biosynthesis and signaling genes are direct targets of HBI1, including

514

BR6OX1 (Fan et al. 2014), indicating the possible involvement of this bHLH

515

regulatory module in promoting BR responses during shade and ethylene exposure.

516

Why is there be such an elaborate network of regulators and even hormones

517

involved in controlling unidrectional cell expansion in hypocotyl growth responses?

518

To achieve a controlled growth, feedback loops are likely required, and crosstalk

519

between different routes are probably a necessity to deal with multiple environmental

520

inputs simultaneously. We found BAS1 transcriptional upregulation in hypocotyls in

521

response to ethylene and shade. BAS1 may act to balance the hypocotyl growth

522

promotion mediated by brassinosteroids (castasterone (CS) and brassinolide (BL)) 16 Downloaded from www.plantphysiol.org on June 22, 2016 - Published by www.plantphysiol.org Copyright © 2016 American Society of Plant Biologists. All rights reserved.

523

as it inactivates both CS and BL (Neff et al., 1999; Turk et al., 2005). In shade

524

avoidance research, likewise, HFR1 is induced by PIFs to suppress the growth

525

promotion induced the same PIF proteins, and putative control of DELLAs would

526

modulate GA responses.

527

Conclusion

528

Hypocotyl elongation in response to ethylene and shade treatments is likely

529

regulated at the upstream level by (a) a bHLH module consisting of positive growth

530

regulators, PIFs (PIF3 in ethylene and PIF4 and PIF5 in shade) and inhibitory factors

531

like IBL1 and (b) a homeodomain module, where ZF-HD TFs like ATHB28 may

532

either act in parallel to the bHLH TFs or are regulated by PIFs (similar to induction of

533

HD-ZIP TFs) to transcriptionally target genes related to the growth promoting

534

hormone module (auxin, BR and GA), as hinted by promoter motif analysis. We

535

hypothesize that in Arabidopsis seedlings, shade and ethylene stimulate auxin

536

synthesis in the cotyledons, which is then transported to the hypocotyl to epidermal

537

cell layers where it interacts with both GA and BR to co-ordinately induce hypocotyl

538

elongation. This increased auxin response, indicated by elevated SAUR levels, and

539

likely increased levels of GA and BR as indicated by increased GA20OX1 and

540

BR6OX expression in the hypocotyl, likely act to induce unidirectional epidermal cell

541

wall elongation via upregulation of genes encoding cell-wall modifying proteins,

542

which promote cellular expansion leading to hypocotyl elongation.

543 544

MATERIALS AND METHODS

545

Plant Material and Growth Conditions

546

Around 30 Arabidopsis thaliana (Col-0) and mutant seeds were sown per agar plate

547

containing 1.1 g L-1 Murashige-Skoog (1/4 MS) and 8 g L-1 Plant-agar (0.8%w/v)

548

(both Duchefa Biochemie, The Netherlands). Mutants or overexpression lines used

549

in this work were: pin3pin4pin7 (Blilou et al., 2005); wei8-1 (Stepanova et al., 2008),

550

tir1-1 (Ruegger et al., 1998), dwf4-1 (Azpiroz et al., 1998), rot3-1 (Kim et al., 1998),

551

ga1-3 (Wilson et al., 1992), gai (Talon et al., 1990), ein3eil1 (Binder et al., 2004);

552

bri1-116 (van Esse et al., 2012); ibl1 and IBL1OE (Zhiponova et al., 2014);

553

cyp85a1/cyp85a2 and cyp85a1-2 (Nomura et al., 2005); bas1-2 (Turk et al., 2005);

554

ga20ox1-3 (Hisamatsu et al., 2005). 17 Downloaded from www.plantphysiol.org on June 22, 2016 - Published by www.plantphysiol.org Copyright © 2016 American Society of Plant Biologists. All rights reserved.

555

ibl1 (N657437), cyp85a1 (N681535), wei8-1 (N16407) and ga20ox1-3 (N669422)

556

were obtained from NASC seed-stock centre; bri1-116 was kindly provided by Sacco

557

de Vries Lab, Wageningen University while cyp85a1/cyp85a2 and dwf4-1 were

558

kindly provided by Sunghwa Choe Lab, Seoul National University. IBL1OE

559

overexpression lines were kindly provided by the Jenny Russinova lab (VIB Ghent,

560

Belgium). bas1-2 mutant was kindly provided by Michael Neff lab (Washington State

561

University, USA). Some GA mutants used were in the Ler background. All other

562

mutants were in the Col-0 background. After stratification at 4°C for 3 d in the dark,

563

seeds were transferred for 2 h to control light conditions (see below) and then kept in

564

dark (at 20°C) again for another 15 h. Subsequently seedlings were allowed a period

565

of 24 h growth under control light conditions in Short-Day photoperiod conditions (15

566

h dark/9h light) before being transferred to 22.4 L glass desiccators with air-tight lids

567

for specific treatments. Col-0 genotypes were grown at 21 ± 1 °C. Ler genotypes

568

were grown at 19 ± 1 °C (Supplemental Figure S6).

569

For athb28 (GK-326G12), lines were obtained from NASC (UK). Genotyping was

570

performed

571

(CTAAGTACCGGGAATGTCAGAAG);

572

AGCTA) and LB primer o8474: (ATAATAACGCTG CGGA CATCTACATTTT). For

573

verifying transcript levels, ATHB28_fwd (GGAGAAGATGAAGGAATTTGCA) and

574

ATHB28_rev (TGTTTCTCTTCA TTGCTTGCT) were used.

using

the

following

primers: athb28_rev

for

athb28,

athb28_fwd

(TAACCAACTGAGCTATTCC

575 576

Treatments

577

Control

578

(Photosynthetically Active Radiation, PAR = 140 μmol m-2 s-1, blue light (400 nm-500

579

nm) = 29 μmol m-2 s-1 and R:FR = 2.1). Ethylene treatments were started by injecting

580

ethylene into the desiccators (with 1 µl L-1 final concentration in the dessicator) and

581

levels were verified with a Gas Chromatograph (GC955; Synspec, The Netherlands).

582

Shade treatment was started by putting desiccators under a single layer of Lee Fern

583

Green Filter (Lee Hampshire, UK) (PAR = 40 μmol m-2s-1, blue light (400 nm-500

584

nm) = 3 μmol m-2 s-1 and R:FR = 0.45). For growth curve experiments, two plates

585

with 15 seedlings per treatment per genotype distributed over two desiccators were

586

used. For mutant analyses, one plate with 15 seedlings per treatment per genotype

587

was used.

and

ethylene

desiccators

were

kept

in

control

light

conditions

588 18 Downloaded from www.plantphysiol.org on June 22, 2016 - Published by www.plantphysiol.org Copyright © 2016 American Society of Plant Biologists. All rights reserved.

589

Imaging and Hypocotyl Length Measurements

590

For hypocotyl elongation assays, experiments were replicated twice. Seedling plates

591

were collected from the dessicators. Seedlings were flattened on the agar plates to

592

reveal the full extent of their hypocotyl, and images of the seedlings were obtained

593

by scanning the plates using EPSON Perfection v370 Photo scanner (Epson Europe

594

BV, The Netherlands). Hypocotyl lengths were measured from these images using

595

ImageJ (http://rsbweb.nih.gov/ij/) and values (per data point) were obtained for n ≥

596

30-60 seedlings. Final values for the data was obtained by taking mean ± S.E. for

597

values from the two independent experiments.

598 599

Epidermal Cell Length Measurements

600

Seedlings were mounted on microscopic slides and covered with a cover slip.

601

Hypocotyl epidermal cells were imaged using Olympus AX70 (20x objective, Nikon

602

DXM1200 camera), after which cell lengths were measured using ImageJ software

603

tool (http://rsbweb.nih.gov/ij/).

604 605

Microarray Tissue Harvest, RNA Isolation and Array Hybridization

606

Seedlings were dissected using BD PrecisionGlide Hypodermic 27 Gauge 1 1/4"

607

Grey Needle with outer diameter = 0.41 mm (Becton Dickinson B.V., The

608

Netherlands) to separately harvest the hypocotyl and cotyledon + shoot apical

609

meristem (hereafter termed “cotyledon”). The roots were discarded. Samples were

610

harvested at 1.5 h, 13.5 h and 25.5 h after the start of the treatments. For 13.5 h time

611

point, which occurs during the dark period, dissection was carried out under low

612

intensity green safelight (≈ 5 μmol m-2 s-1). For minimizing the effects of green light

613

on gene expression, seedlings were kept in the dark until dissected and dissection

614

was carried out in several rounds with each round involving maximum 2-3 seedlings.

615

In total, 3 replicate experiments were carried out. In each replicate experiment,

616

tissues were harvested from two independent technical replicates (each with 25

617

seedlings from the two technical replicate plates as mentioned in “Treatments”

618

section above). Harvested material was immediately frozen in liquid nitrogen and

619

stored at -80°C until further use.

620

Frozen tissue was ground using tissue lyser and total RNA was isolated using

621

RNeasy Mini Kit (QIAGEN, Germany). QIAGEN RNase-Free DNase set was used to

622

eliminate Genomic DNA contamination by performing on-column DNase digestion. 19 Downloaded from www.plantphysiol.org on June 22, 2016 - Published by www.plantphysiol.org Copyright © 2016 American Society of Plant Biologists. All rights reserved.

623

Extracted RNA was verified (for quality check) and quantified using NanoDropTM ND-

624

1000 Spectrophotometer (Isogen Life Science, De Meern, The Netherlands).

625

RNA samples were sent to AROS (AROS Applied Biotechnology A/S, Aarhus,

626

Denmark). RNA was repurified on low-elution QIAGEN RNeasy columns, re-

627

quantified with NanoDropTM 8000 UV-Vis Spectrophotometer and checked for quality

628

with Agilent 2100 Bioanalyzer (Agilent Technologies, Palo Alto, CA). RNA samples

629

with RNA Integrity Number (RIN) value of > 7.5 were considered for further use. 50

630

ng RNA from each of the two independent technical replicates was mixed in the ratio

631

1:1 and the pooled sample was considered as one biological replicate for

632

hybridization experiments. Thus, 3 biological replicates were obtained for the 3

633

replicate experiments. 100 ng of RNA sample was processed for cDNA synthesis,

634

fragmentation and labelling was carried out for the RNA samples. The samples were

635

hybridized to the Affymetrix Gene 1.1 ST Arabidopsis Array Plate and washed on an

636

Affymetrix GeneAtlas system followed by scanning of arrays at AROS Applied

637

Biotechnology (http://arosab.com/services/microarrays/gene-expression/).

638 639

Microarray Data Analysis

640

Scanned arrays in the form of .cel files (provided by AROS) were checked for quality

641

control using Affymetrix Expression Console Software and an in-house script in R

642

and

643

(Bioconductor “oligo” and “pd.aragene.1.1.st” ). Bioconductor was used for Robust

644

Multi-array Average (RMA) normalization of raw data at Gene level to obtain

645

summarized signal intensity values for all genes present on the array (log2 format).

646

Principal component analysis was carried out using Affymetrix Expression Console

647

Software (http://www.affymetrix.com/) and dendrogram of all microarray samples

648

according to the mean signal intensity values was generated using R (“plot”

649

package). Bioconductor (“Limma” package) was used for carrying out differential

650

expression analysis.

Bioconductor

(http://www.r-project.org./;

http://www.bioconductor.org/)

651 652

Temporal Clustering And Bioinformatics Analysis

653

We clustered the list of total DEGs (defined as number of DEGs that were regulated

654

at atleast one of the three time points) under ethylene and shade based on positive

655

or negative regulation at each of the three time-points. With three time-points and

656

two directions of expression (positive or negative), 23 = 8 possible trends can occur 20 Downloaded from www.plantphysiol.org on June 22, 2016 - Published by www.plantphysiol.org Copyright © 2016 American Society of Plant Biologists. All rights reserved.

657

and accordingly as many clusters were obtained. DEGs were clustered temporally

658

based on log2FC.

659 660

MapMan bin overrepresentation using hypergeometric test was done using R (“stats”

661

package) and adjusted p-values for the statistical significance of enrichment were

662

converted into negative logarithm score and plotted as a heatmap. A score above

663

1.3 for this score was considered significant.

664

GeneCodis (http://genecodis.cnb.csic.es/compareanalysis) webtool was used for GO

665

analysis of different sets obtained from venn of Positive clusters and of Negative

666

Clusters. A score above 1.3 for negative log of adjusted p-values was considered

667

significant.

668 669

LogFC Filter and Gene Classification

670

In order to narrow down the genes for functional characterization (from the list of

671

classified genes), we utilised a logFC filter. To narrow down the Common Up set, a

672

filter of log2FC < 0.5 at 1.5 h and log2FC ≥ 0.5 at both 13.5 h and 25.5 h for ethylene

673

and log2FC ≥ 0.5 at both 1.5 h and 25.5 h for shade was applied. To narrow down

674

the Common Down set, we applied a filter of log2FC > - 0.5 at 1.5 h and log2FC ≤ -

675

0.5 at both 13.5 h and 25.5 h under ethylene and of log2FC ≤ - 0.5 at both 1.5 h and

676

25.5 h under shade. These resulting group of genes were then classified based on

677

MapMan based classification for the terms: “RNA.regulation of transcription”, “Cell

678

wall”,

679

(http://planttfdb.cbi.pku.edu.cn/index.php?sp=Ath)

680

(http://plntfdb.bio.uni-potsdam.de/v3.0/index.php?sp_id=ATH) were also included as

681

a source for gene-classification to select for additional transcriptional regulators.

“Signalling”

and

“Hormone

metabolism”.

Genes

from

Plant

and

Potsdam

TFDB TFDB

682 683

Hormone Correlational Analysis

684

Hormonometer software (http://genome.weizmann.ac.il/hormonometer/) was used to

685

evaluate transcriptional similarities between the transcriptome data obtained here

686

and the published, indexed list of those elicited by exogenous application of plant

687

hormones. Arabidopsis gene locus IDs were converted to Affymetrix GeneChip

688

identifiers using the “at to AGI converter” tool (The Bio-Analytic Resource for Plant

689

Biology, http://bar.utoronto.ca/). We used the new Affymetrix aragen1.1st arrays (28k

690

genes) for transcriptomics but the hormonometer data is based on 3’ ATH1 arrays 21 Downloaded from www.plantphysiol.org on June 22, 2016 - Published by www.plantphysiol.org Copyright © 2016 American Society of Plant Biologists. All rights reserved.

691

(22k genes). Accordingly, many locus IDs could not be included (those which were

692

newly incorporated in the aragene1.1st arrays) in this analysis. In few other cases

693

where multiple ATH1 GeneChip IDs for one locus ID was obtained, all IDs were

694

retained.

695 696

Pharmacological Treatments

697

Auxin transport was inhibited by use of 25 µM NPA (Duchefa Biochemie, The

698

Netherlands) (Petrásek et al., 2003). Auxin perception was blocked by use of 100 µM

699

PEO-IAA (Hayashi et al., 2008). Auxin biosynthesis was blocked by use of 50 µM

700

yucasin (Nishimura et al., 2014). 2 µM brassinazole (TCI Europe, Japan) was used

701

to inhibit BR biosynthesis (Asami et al., 2000). 2 µM paclobutrazol (Duchefa

702

Biochemie, The Netherlands) was used to inhibit GA biosynthesis (Rademacher,

703

2000). All chemicals were dissolved in respective solvents (DMSO or ETOH) with

704

final solvent concentration in media < 0.1% to prevent toxicity due to solvents. All

705

chemicals were applied by pipetting 150 µl of chemical solutions or mock solvents as

706

a thin film over the MS-agar media in the petri plates and then allowing the solution

707

to diffuse through the medium before starting the treatments.

708 709

GUS Staining and Imaging

710

For GUS assays, seedlings were transferred immediately from treatments to a GUS

711

staining solution (1 mM X-Gluc (5-bromo-4-chloro-3-indolyl-D-glucoronide) (Duchefa

712

Biochem, The Netherlands) in 100 mM Sodium phosphate buffer (pH=7.0) along with

713

0.1 % Triton X-100, 0.5 mM each of Potassium ferrocyanide (K4Fe(CN)6 and

714

Potassium ferricyanide (K3Fe(CN)6 and 10 mM of EDTA (Merck Darmstadt,

715

Germany)) and kept at 37°C overnight. Seedlings were bleached in 70% ethanol for

716

1 day before capturing images.

717 718

Statistical Analysis and Graphing

719

1-way ANOVAs followed by Tukey's HSD Post-Hoc tests were performed on the

720

measurements obtained in hypocotyl length/cotyledon area kinetics to assess

721

statistically significant differences between mean hypocotyl length/cotyledon area

722

under ethylene or shade relative to control at the same time point.

723

A two-way ANOVA followed by a Tukey's HSD Post-Hoc test was used for pairwise

724

multiple comparison. For hypocotyl elongation assays, statistical significance was 22 Downloaded from www.plantphysiol.org on June 22, 2016 - Published by www.plantphysiol.org Copyright © 2016 American Society of Plant Biologists. All rights reserved.

725

indicated by the use of different letters. All statistical analyses were done in the R

726

software environment. Graphs were plotted using Prism 6 software (GraphPad

727

Software, USA).

728 729

Accession numbers

730

.CEL files utilised in the organ-specific transcriptomics for hypocotyl and cotyledon

731

tissues are posted with the GEO accession series GSE83212.

732 733

Supplemental data

734

The following materials are available in the online version of this article.

735

Supplemental Figure S1. Hypocotyl elongation response in wild-type Col-0 under

736

control, ethylene, shade and combination (ethylene + shade) treatments.

737

Supplemental Figure S2. Venn diagram of gene intersection between up- and

738

down- regulated DEGs of ethylene and shade separately at 1.5 h.

739

Supplemental Figure S3. Genotyping and transcript level verification for (A)-(D)

740

athb28.

741

Supplemental Figure S4. Hypocotyl elongation response in ethylene signalling

742

mutant, ein3eil1 under ethylene and shade.

743

Supplemental Figure S5. TAIR motif analysis for Common Up and Down genesets.

744

Supplemental Figure S6. Hypocotyl length of Ler under control and ethylene

745

conditions when grown at 220 C day / 200 C night and 200 C day / 180 C night

746

temperature regime.

747

Supplemental Table S1. Photosynthesis gene proportion in genesets from Venn of

748

negative clusters.

749

ACKNOWLEDGMENTS

750

We would like to thank all the group members of Plant Ecophysiology, Utrecht

751

University, The Netherlands for their help in the harvest for the transcriptomics

752

experiment.

23 Downloaded from www.plantphysiol.org on June 22, 2016 - Published by www.plantphysiol.org Copyright © 2016 American Society of Plant Biologists. All rights reserved.

753 754

FIGURE LEGENDS

755

Figure 1. Physiological responses, hypocotyl length and epidermal cell length

756

kinetics under ethylene and shade in Arabidopsis (Col-0) seedlings. 1-day old

757

seedlings were exposed to control conditions (PAR = 140 μmol m-2 s-1, blue light =

758

29 μmol m-2 s-1and R:FR = 2.1), elevated ethylene (1 µl L-1, PAR = 140 μmol m-2 s-1,

759

blue light = 29 μmol m-2 s-1 and R:FR = 2.1) or shade (PAR = 40 μmol m-2 s-1, blue

760

light = 3 μmol m-2 s-1 and R:FR = 0.45). A, Figure shows representative seedlings

761

displaying typical phenotypes after 96 h of exposure to control, ethylene or shade

762

conditions. B, Mean hypocotyl lengths at 0 h, 24 h, 48 h, 72 h and 96 h under control

763

(open circle), ethylene (closed circle) and shade (open triangle) are shown. C, Rate

764

of increase in hypocotyl length. Differences between mean hypocotyl lengths of

765

subsequent time points averaged over 1 d time interval for control, ethylene and

766

shade are shown. D, Detailed hypocotyl length kinetics. 1-day old seedlings were

767

exposed to control (open circle), ethylene (closed circle) and shade (open triangle)

768

and measured at 3 h time intervals (s, e: first time point at which shade (s) or

769

ethylene (e) treatment lead to significantly longer hypocotyls). Data are mean ± S.E.

770

(n=60) for (A)-(D). Shaded area denotes the 15 h dark period in the 15 h dark/9 h

771

light photoperiodic growth condition. e and s denote the first point of statistically

772

significant differences in hypocotyl length or cotyledon area relative to control for

773

ethylene and shade respectively. E, Epidermal cell length kinetics. 1-day old

774

seedlings were exposed to control, ethylene or shade. Mean cell length ± S.E. (n ≥

775

10) for epidermal cells of the Arabidopsis hypocotyl at 0 h control (black line) and at

776

3 h, 7.5 h, 15 h and 27 h control (orange line), ethylene (blue line) and shade (green

777

line) is shown. Apex denotes the hypocotyl-cotyledon junction and base denotes the

778

hypocotyl-root junction.

779

Figure 2. Overall description of microarray data. A, Principal component analysis

780

(PCA) and hierarchical clustering (HC) was used to describe the structure in the

781

microarray data. Expression intensities for all genes on the array for all 54 hypocotyl

782

and cotyledon samples (3 time points, 3 treatments and 3 replicates) were projected

783

onto the first three principal components. B, Hierarchical clustering was used to

784

group 18 main samples (according to mean expression intensity of 3 replicates for

24 Downloaded from www.plantphysiol.org on June 22, 2016 - Published by www.plantphysiol.org Copyright © 2016 American Society of Plant Biologists. All rights reserved.

785

each main sample) into a dendrogram. C, Distribution of differentially expressed

786

genes (DEGs) in hypocotyl and cotyledon samples at three time-points in response

787

to ethylene and shade. DEGs obtained for each time-point were plotted separately

788

as up-regulated, down-regulated (adjusted p-value ≤ 0.01 and log2FC > or < 0) or

789

non-significant (adjusted p-value > 0.01). Bar length denotes total DEGs obtained

790

after combining DEGs from all 3 time points.

791

Figure 3. Temporal gene expression clusters for ethylene (A) and shade (D);

792

hypocotyl growth curve for ethylene (B) and shade (E); clusters with gene expression

793

matching hypocotyl growth kinetics in ethylene (C) and shade (F). Heatmap for

794

temporal clusters (A, D) based on the log2FC at the three microarray time points

795

(grey box represents dark; arrows indicate heatmap time-points and treatments

796

include control (open circle), ethylene (closed circle), shade (open triangle)). Yellow

797

denotes up-regulation and blue denotes down-regulation. Two gene expression

798

clusters (C, F) with mean log2FC temporal pattern resembling the hypocotyl length

799

kinetics (B, E) were named positive clusters. G, Heatmap for hypergeometric

800

enrichment of selected MapMan bins for temporal clusters under ethylene and

801

shade. Horizontal-axis denotes the cluster number. More intense colours indicate

802

higher statistical significance. Grey color indicates non-significant score or absence

803

of genes in the bin. H, Venn intersection for positive clusters (clusters with gene

804

expression pattern matching the hypocotyl length kinetics) from ethylene and shade

805

to obtain “Common Up” genes. I, Venn intersection for negative clusters (mirror

806

image to Positive clusters) from ethylene and shade to obtain “Common Down”

807

genes.

808

Figure 4. Gene Ontology (GO) enrichment analysis using GeneCodis for positive

809

and negative clusters. Adjusted p-values for statistical significance of GO enrichment

810

were converted into negative logarithm score (base10) (> 1.3 are considered

811

significant). Heatmap colours denote this score and more intense colours indicate

812

higher statistical significance. White color indicates non-significant score or absence

813

of genes in the GO terms.

814 815

Figure 5. Heatmap of (A) Common Up and (B) Common Down set of genes

816

classified into the categories: transcriptional regulator, hormone metabolism,

25 Downloaded from www.plantphysiol.org on June 22, 2016 - Published by www.plantphysiol.org Copyright © 2016 American Society of Plant Biologists. All rights reserved.

817

signalling and cell wall genes after application of a log2fold-change (log2FC) filter.

818

Log2FC at the three time-points 1.5 h, 13.5 h, 25.5 h for ethylene and shade

819

microarray dataset identified in the color scheme of the heatmap. Genes were

820

categorized according to Mapman bin gene classification (shown on left side).

821 822

Figure 6. Hypocotyl length measurements for (A) athb28 (B) ibl1 (C) IBL1-OE

823

following 96 h of control, ethylene and shade treatments. Data represents mean ±

824

S.E. (n=30 seedlings). Different letters above the bars indicate significant differences

825

(2-way ANOVA followed by Tukey’s HSD Post-Hoc pairwise comparison).

826 827

Figure 7. A, Identification of enriched hormonal signatures in the ethylene- and

828

shade-induced Arabidopsis transcriptome. Ethylene and shade-induced hypocotyl

829

and cotyledon transcriptomes were analyzed for hormonal signatures using the tool

830

Hormonometer (Valdorsky et al., 2009) to establish correlations with expression data

831

in an established hormonal transcriptome database. Positive correlations were

832

colored yellow and negative correlations blue. Significant correlations were identified

833

with absolute correlation values of 0.3 and higher. Numbers in the cells represent the

834

exact correlation values. Rows denote hormone treatments that are indicated by the

835

name of the hormone and the duration of hormone treatment. Columns denote

836

ethylene and shade transcriptome in the hypocotyl and cotyledon at the three time-

837

points of tissue harvest. The magnitude of correlation in gene expression is indicated

838

by the color scale at top right side. B, Effect of auxin transport inhibitor NPA (25 µM),

839

IAA (10 µM), NPA (25 µM) + IAA (10 µM) and auxin perception inhibitor PEO-IAA

840

(100 µM) on GUS staining of the pIAA19:GUS lines. For NPA and PEO-IAA effect in

841

GUS assay, seedlings were exposed to 2 days of treatment conditions. Arabidopsis

842

(Col-0) seedlings were treated with chemical inhibitors for (C) auxin transport (NPA),

843

(D)

844

brassinosteroid biosynthesis (BRZ) and (G) gibberellin biosynthesis (PBZ) at the

845

indicated concentrations in the legend and length was measured following 96h of

846

ethylene and shade. Hypocotyl length was measured following 96 h of ethylene and

847

shade. Mean ± S.E. was calculated for 30 seedlings. Different letters above the bars

848

indicate significant differences from a 2-way ANOVA followed by Tukey’s HSD Post-

849

Hoc pairwise comparison. Abbreviations: PEO-IAA is α-(phenylethyl-2-one)-indole-3-

auxin

biosynthesis

(Yucasin),

(E)

auxin

perception

(PEO-IAA),

(F)

26 Downloaded from www.plantphysiol.org on June 22, 2016 - Published by www.plantphysiol.org Copyright © 2016 American Society of Plant Biologists. All rights reserved.

850

acetic acid; NPA is 1-N-Naphthylphthalamic acid; BRZ is Brassinazole; PBZ is

851

Paclobutrazol.

852 853

Figure 8. Hypocotyl elongation response in response to ethylene and shade in (B)–

854

(D) auxin, (E)–(G) gibberellin (GA) and (H)–(J) brassinosteroid (BR) mutants. Mean

855

± S.E. was calculated for 30 seedlings. Different letters above the bars indicate

856

significant differences from a 2-way ANOVA followed by Tukey’s HSD Post-Hoc

857

pairwise comparison.

858 859 860 861 862 863 864 865 866 867 868 869 870 871 872 873 874 875 876 877 878 879 880 881 882 883 884 885 886 887 888 889 890 891 892 893

LITERATURE CITED Abdi, H., and Williams, L. J. (2010). Principal component analysis. Wiley Interdisciplinary Reviews: Computational Statistics, 2(4), 433-459. Achard, P., Cheng, H., Grauwe, L. De, Decat, J., Schoutteten, H., Moritz, T., Straeten, D. Van Der, Peng, J., and Harberd, N.P. (2006). Integration of Plant Responses to Environmentally Activated Phytohormonal Signals. Science 311: 91–94. Asami, T., Min, Y.K., Nagata, N., Yamagishi, K., Takatsuto, S., Fujioka, S., Murofushi, N., Yamaguchi, I., and Yoshida, S. (2000). Characterization of brassinazole, a triazole-type brassinosteroid biosynthesis inhibitor. Plant Physiol. 123: 93–100. Azpiroz, R., Wu, Y., LoCascio, J.C., and Feldmann, K.A. (1998). An Arabidopsis brassinosteroid-dependent mutant is blocked in cell elongation. The Plant Cell 10: 219–230. Bai, M.Y., Fan, M., Oh, E., and Wang, Z.Y. (2012) A triple helix-loop-helix/basic helix-loop-helix cascade controls cell elongation downstream of multiple hormonal and environmental signaling pathways in Arabidopsis. Plant Cell. 24: 4917-4929 Barboza, L., Effgen, S., Alonso-blanco, C., Kooke, R., Keurentjes, J.J.B., and Koornneef, M. (2013). Arabidopsis semidwarfs evolved from independent mutations in GA20ox1, ortholog to green revolution dwarf alleles in rice and barley. Proc. Natl. Acad. Sci. U. S. A. 110: 15818–15823. Binder, B.M., Mortimore, L.A., Stepanova, A.N., Ecker, J.R., and Bleecker, A.B. (2004). Short-Term Growth Responses to Ethylene in Arabidopsis. Plant Physiol. 136: 2921–2927. Blilou, I., Xu, J., Wildwater, M., Willemsen, V., Paponov, I., Friml J, Heidstra R., Aida M., Palme K. & Scheres, B. (2005). The PIN auxin efflux facilitator network controls growth and patterning in Arabidopsis roots. Nature, 433(7021), 39-44. Cagnola, J. I., Ploschuk, E., Benech-Arnold, T., Finlayson, S. A., and Casal, J. J. (2012). Stem transcriptome reveals mechanisms to reduce the energetic cost of shade-avoidance responses in tomato. Plant Physiol., 160(2), 1110-1119. Capella, M., Ribone, P. A., Arce, A. L., and Chan, R. L. (2015). Arabidopsis thaliana HomeoBox 1 (AtHB1), a Homedomain‐Leucine Zipper I (HD‐Zip I) 27 Downloaded from www.plantphysiol.org on June 22, 2016 - Published by www.plantphysiol.org Copyright © 2016 American Society of Plant Biologists. All rights reserved.

894 895 896 897 898 899 900 901 902 903 904 905 906 907 908 909 910 911 912 913 914 915 916 917 918 919 920 921 922 923 924 925 926 927 928 929 930 931 932 933 934 935 936 937 938 939 940 941

transcription factor, is regulated by PHYTOCHROME‐INTERACTING FACTOR 1 to promote hypocotyl elongation. New Phytologist, 207, 669-682. Casal, J.J. (2013). Photoreceptor signaling networks in plant responses to shade. Annu. Rev. Plant Biol. 64: 403–27. Casal, J.J. (2012). Shade avoidance. Arabidopsis Book 10: e0157. Chae, K., Isaacs, C.G., Reeves, P.H., Maloney, G.S., Muday, G.K., Nagpal, P., and Reed, J.W. (2012). Arabidopsis SMALL AUXIN UP RNA63 promotes hypocotyl and stamen filament elongation. Plant J. 71: 684–697. Chapman, E.J., Greenham, K., Castillejo, C., Sartor, R., Bialy, A., Sun, T.P., and Estelle, M. (2012). Hypocotyl transcriptome reveals auxin regulation of growthpromoting genes through GA-dependent and -independent pathways. PLoS One 7: e36210. De Smet, I., Lau, S., Ehrismann, J.S., Axiotis, I., Kolb, M., Kientz, M., Weijers, D., and Jurgens, G. (2013). Transcriptional repression of BODENLOS by HDZIP transcription factor HB5 in Arabidopsis thaliana. J. Exp. Bot. 64: 3009–3019. Djakovic-Petrovic, T., Wit, M. De, Voesenek, L.A.C.J., and Pierik, R. (2007). DELLA protein function in growth responses to canopy signals. Plant J. 51: 117– 126. Eisen, M.B., Spellman, P.T., Brown, P.O., and Botstein, D. (1998). Cluster analysis and display of genome-wide expression patterns. Proc. Natl. Acad. Sci. U. S. A. 95: 14863–14868. Endo, M., Nakamura, S., Araki, T., Mochizuki, N., and Nagatani, A. (2005). Phytochrome B in the Mesophyll Delays Flowering by Suppressing FLOWERING LOCUS T Expression in Arabidopsis Vascular Bundles. The Plant Cell 17: 1941–1952. Estelle, M. (1998). Polar Auxin Transport: New Support for an Old Model. The Plant Cell 10: 1775–1778. Fan, M., Bai, M. Y., Kim, J. G., Wang, T., Oh, E., Chen, L., Park, C.H., Son S., Kim S., Mudgett M.B. and Wang, Z. Y. (2014). The bHLH transcription factor HBI1 mediates the trade-off between growth and pathogen-associated molecular pattern–triggered immunity in Arabidopsis. The Plant Cell Online, 26(2), 828-841. Franklin, K. A., Lee, S.H., Patel, D., Kumar, S.V., Spartz, A.K., Gu, C., Ye, S., Yu, P., Breen, G., Cohen, J.D., Wigge, P.A., and Gray, W.M. (2011). Phytochrome-interacting factor 4 (PIF4) regulates auxin biosynthesis at high temperature. Proc. Natl. Acad. Sci. U. S. A. 108: 20231–5. Franklin, K. A. (2008). Shade avoidance. New Phytol. 179: 930–944. Gendreau, E., Traas, J., Desnos, T., Grandjean, O., Caboche, M., and Höfte, H. (1997). Cellular basis of hypocotyl growth in Arabidopsis thaliana. Plant Physiol. 114: 295–305. Hattori, Y., Nagai, K., Furukawa, S., Song, X., Kawano, R., Sakakibara, H., Wu, J., Matsumoto, T., Yoshimura, A., Kitano, H., Matsuoka, M., Mori, H., Ashikari, M. (2009). The ethylene response factors SNORKEL1 and SNORKEL2 allow rice to adapt to deep water. Nature 460: 1026–1030. Hayashi, K., Tan, X., Zheng, N., Hatate, T., Kimura, Y., Kepinski, S., and Nozaki, H. (2008). Small-molecule agonists and antagonists of F-box protein-substrate interactions in auxin perception and signaling. Proc. Natl. Acad. Sci. U. S. A. 105: 5632–5637.

28 Downloaded from www.plantphysiol.org on June 22, 2016 - Published by www.plantphysiol.org Copyright © 2016 American Society of Plant Biologists. All rights reserved.

942 943 944 945 946 947 948 949 950 951 952 953 954 955 956 957 958 959 960 961 962 963 964 965 966 967 968 969 970 971 972 973 974 975 976 977 978 979 980 981 982 983 984 985 986 987 988 989 990 991

Henriques, R., Bögre, L., Horváth, B., and Magyar, Z. (2014). Balancing act: Matching growth with environment by the TOR signalling pathway. J. Exp. Bot. 65: 2691–2701. Hisamatsu, T., King, R.W., Helliwell, C. a, and Koshioka, M. (2005). The involvement of Gibberellin 20-Oxidase genes in phytochrome-regulated petiole elongation of Arabidopsis. Plant Physiol. 138: 1106–1116. Hong, S.Y., Kim, O.K., Kim, S.G., Yang, M.S., and Park, C.M. (2011). Nuclear import and DNA binding of the ZHD5 transcription factor is modulated by a competitive peptide inhibitor in Arabidopsis. J. Biol. Chem. 286: 1659–1668. Hu, W. and Ma, H. (2006). Characterization of a novel putative zinc finger gene MIF1: Involvement in multiple hormonal regulation of Arabidopsis development. Plant J. 45: 399–422. Hyun, Y. and Lee, I. (2006). KIDARI, encoding a non-DNA Binding bHLH protein, represses light signal transduction in Arabidopsis thaliana. Plant Mol. Biol. 61: 283–96. Iqbal, N., Khan, N. A., Nazar, R., and da Silva, J. A. T. (2012). Ethylene-stimulated photosynthesis results from increased nitrogen and sulfur assimilation in mustard types that differ in photosynthetic capacity. Environ. Exp. Bot. 78: 84– 90. Kazan, K. and Manners, J.M. (2013). MYC2: The master in action. Mol. Plant 6: 686–703. Keuskamp, D.H., Pollmann, S., Voesenek, L.A.C.J., Peeters, A.J.M., and Pierik, R. (2010). Auxin transport through PIN-FORMED 3 (PIN3) controls shade avoidance and fitness during competition. Proc. Natl. Acad. Sci. U. S. A. 107: 22740–22744. Keuskamp, D.H., Sasidharan, R., Vos, I., Peeters, A.J.M., Voesenek, L.A.C.J., and Pierik, R. (2011). Blue-light-mediated shade avoidance requires combined auxin and brassinosteroid action in Arabidopsis seedlings. Plant J. 67: 208–217. Khanna, R., Shen, Y., Marion, C.M., Tsuchisaka, A., Theologis, A., Schafer, E., and Quail, P.H. (2007). The basic helix-loop-helix transcription factor PIF5 acts on ethylene biosynthesis and phytochrome signaling by distinct mechanisms. Plant Cell. 19: 3915-3929. Kim, B.C., Soh, M.S., Hong, S.H., Furuya, M., and Nam, H.G. (1998). Photomorphogenic development of the Arabidopsis shy2-1D mutation and its interaction with phytochromes in darkness. Plant J. 15: 61–68. Kunihiro, A., Yamashino, T., Nakamichi, N., Niwa, Y., Nakanishi, H., and Mizuno, T. (2011). PHYTOCHROME-INTERACTING FACTOR 4 and 5 (PIF4 and PIF5) activate the Homeobox ATHB2 and auxin-inducible IAA29 genes in the coincidence mechanism underlying photoperiodic control of plant growth of Arabidopsis thaliana. The Plant Cell Physiol. 52: 1315–1329. Lee, Y.K., Kim, G., Kim, I., Park, J., Kwak, S., Choi, G., and Chung, W. (2006). LONGIFOLIA1 and LONGIFOLIA2, two homologous genes, regulate longitudinal cell elongation in Arabidopsis. Development 133: 4305–4314. Leivar, P. and Monte, E. (2014). PIFs: systems integrators in plant development. The Plant Cell 26: 56–78. Li, L. et al. (2012). Linking photoreceptor excitation to changes in plant architecture. Genes Dev. 26: 785–790. Lilley, J.L.S., Gee, C.W., Sairanen, I., Ljung, K., and Nemhauser, J.L. (2012). An Endogenous Carbon-Sensing Pathway Triggers Increased Auxin Flux and Hypocotyl Elongation. Plant Physiol. 160: 2261–2270. 29 Downloaded from www.plantphysiol.org on June 22, 2016 - Published by www.plantphysiol.org Copyright © 2016 American Society of Plant Biologists. All rights reserved.

992 993 994 995 996 997 998 999 1000 1001 1002 1003 1004 1005 1006 1007 1008 1009 1010 1011 1012 1013 1014 1015 1016 1017 1018 1019 1020 1021 1022 1023 1024 1025 1026 1027 1028 1029 1030 1031 1032 1033 1034 1035 1036 1037 1038 1039 1040 1041

Mickelbart, M. V., Hasegawa, P.M., and Bailey-Serres, J. (2015). Genetic mechanisms of abiotic stress tolerance that translate to crop yield stability. Nat. Rev. Genet. 16: 237–251. Morelli, G., Ruberti, I., Molecolare, B., Sapienza, L., and Moro, P.A. (2000). Update on Light Signaling Shade Avoidance Responses. Driving Auxin along Lateral Routes. Plant Physiol. 122: 621–626. Neff, M.M., Nguyen, S.M., Malancharuvil, E.J., Fujioka, S., Noguchi, T., Seto, H., Tsubuki, M., Honda, T., Takatsuto, S., Yoshida, S., and Chory, J. (1999). BAS1: A gene regulating brassinosteroid levels and light responsiveness in Arabidopsis. Proc. Natl. Acad. Sci. 96: 15316–15323. Nemhauser, J.L., Mockler, T.C., and Chory, J. (2004). Interdependency of brassinosteroid and auxin signaling in Arabidopsis. PLoS Biol. 2: E258. Nishimura, T., Hayashi, K.I., Suzuki, H., Gyohda, A., Takaoka, C., Sakaguchi, Y., Matsumoto, S., Kasahara, H., Sakai, T., Kato, J.I., Kamiya, Y., and Koshiba, T. (2014). Yucasin is a potent inhibitor of YUCCA, a key enzyme in auxin biosynthesis. Plant J. 77: 352–366. Nomura, T., Kushiro, T., Yokota, T., Kamiya, Y., Bishop, G.J., and Yamaguchi, S. (2005). The last reaction producing brassinolide is catalyzed by cytochrome P-450s, CYP85A3 in tomato and CYP85A2 in Arabidopsis. J. Biol. Chem. 280: 17873–17879. Nozue, K., Tat, A. V., Kumar Devisetty, U., Robinson, M., Mumbach, M.R., Ichihashi, Y., Lekkala, S., and Maloof, J.N. (2015). Shade Avoidance Components and Pathways in Adult Plants Revealed by Phenotypic Profiling. PLOS Genet. 11: e1004953. Osakabe, Y., Osakabe, K., Shinozaki, K., and Tran, L. S. P. (2014). Response of plants to water stress. Front. Plant Sci. 5(86): 10-3389. Park, E., Park, J., Kim, J., Nagatani, A., Lagarias, J.C., and Choi, G. (2012). Phytochrome B inhibits binding of phytochrome-interacting factors to their target promoters. Plant J. 72: 537–546. Petrásek, J., Cerná, A., Schwarzerová, K., Elckner, M., Morris, D. A., and Zazímalová, E. (2003). Do phytotropins inhibit auxin efflux by impairing vesicle traffic? Plant Physiol. 131: 254–263. Pfeiffer, A., Shi, H., Tepperman, J.M., Zhang, Y., and Quail, P.H. (2014). Combinatorial Complexity in a Transcriptionally Centered Signaling Hub in Arabidopsis. Mol. Plant 7: 1598–1618. Philippar, K., Ivashikina, N., Ache, P., Christian, M., Lüthen, H., Palme, K., and Hedrich, R. (2004). Auxin activates KAT1 and KAT2, two K+-channel genes expressed in seedlings of Arabidopsis thaliana. Plant J. 37: 815–827. Pierik, R., Cuppens, M.L., Voesenek, L.A.C.J. and Visser, E.J. (2004) Interactions between ethylene and gibberellins in phytochrome-mediated shade avoidance responses in tobacco. Plant Physiol. 136: 2928-2936. Pierik, R., Djakovic-Petrovic, T., Keuskamp, D.H., de Wit, M., and Voesenek, L.A.C.J. (2009). Auxin and ethylene regulate elongation responses to neighbor proximity signals independent of gibberellin and DELLA proteins in Arabidopsis. Plant Physiol. 149: 1701–1712. Pierik, R., Millenaar, F.F., Peeters, A. J.M., and Voesenek, L.A.C.J. (2005). New perspectives in flooding research: The use of shade avoidance and Arabidopsis thaliana. Ann. Bot. 96: 533–540. Pierik, R. and Testerink, C. (2014). The art of being flexible: how to escape from shade, salt and drought. Plant Physiol. 166: 5–22. 30 Downloaded from www.plantphysiol.org on June 22, 2016 - Published by www.plantphysiol.org Copyright © 2016 American Society of Plant Biologists. All rights reserved.

1042 1043 1044 1045 1046 1047 1048 1049 1050 1051 1052 1053 1054 1055 1056 1057 1058 1059 1060 1061 1062 1063 1064 1065 1066 1067 1068 1069 1070 1071 1072 1073 1074 1075 1076 1077 1078 1079 1080 1081 1082 1083 1084 1085 1086 1087 1088 1089 1090 1091

Pierik, R. and De Wit, M. (2014). Shade avoidance: Phytochrome signalling and other aboveground neighbour detection cues. J. Exp. Bot. 65: 2815–2824. Procko, C., Crenshaw, C.M., Ljung, K., Noel, J.P., and Chory, J. (2014). Cotyledon-Generated Auxin Is Required for Shade-Induced Hypocotyl Growth in Brassica rapa. Plant Physiol. 165: 1285–1301. Quint, M., Delker, C., Franklin, K. A., Wigge, P. A., Halliday, K. J., and van Zanten, M. (2016). Molecular and genetic control of plant thermomorphogenesis. Nature Plants. 2: 15190. Rademacher, W. (2000). GROWTH RETARDANTS: Effects on Gibberellin Biosynthesis and Other Metabolic Pathways. Annu. Rev. Plant Physiol. Plant Mol. Biol. 51: 501–531. Rieu, I., Ruiz-Rivero, O., Fernandez-Garcia, N., Griffiths, J., Powers, S.J., Gong, F., Linhartova, T., Eriksson, S., Nilsson, O., Thomas, S.G., Phillips, A.L., and Hedden, P. (2008). The gibberellin biosynthetic genes AtGA20ox1 and AtGA20ox2 act, partially redundantly, to promote growth and development throughout the Arabidopsis life cycle. Plant J. 53: 488–504. Ruegger, M., Dewey, E., Gray, W.M., Hobbie, L., Turner, J., and Estelle, M. (1998). The TIR1 protein of Arabidopsis functions in auxin response and is related to human SKP2 and yeast Grr1p. Genes Dev. 12: 198–207. Sasidharan, R., Chinnappa, C.C., Staal, M., Elzenga, J.T.M., Yokoyama, R., Nishitani, K., Voesenek, L.A.C.J., and Pierik, R. (2010). Light qualitymediated petiole elongation in Arabidopsis during shade avoidance involves cell wall modification by xyloglucan endotransglucosylase/hydrolases. Plant Physiol. 154: 978–990. Sasidharan, R. and Voesenek, L.A.C.J. (2015). Ethylene-mediated acclimations to flooding stress. Plant Physiol.: pp.00387. Smalle, J., Haegman, M., Kurepa, J., Van Montagu M, and Straeten, D. V (1997). Ethylene can stimulate Arabidopsis hypocotyl elongation in the light. Proc. Natl. Acad. Sci. U. S. A. 94: 2756–2761. Spartz, A.K., Lee, S.H., Wenger, J.P., Gonzalez, N., Itoh, H., Inzé, D., Peer, W.A., Murphy, A.S., Overvoorde, P.J., and Gray, W.M. (2012). The SAUR19 subfamily of SMALL AUXIN UP RNA genes promote cell expansion. Plant J. 70: 978–990. Stepanova, A.N., Robertson-Hoyt, J., Yun, J., Benavente, L.M., Xie, D.Y., Doležal, K., Schlereth, A., Jürgens, G., and Alonso, J.M. (2008). TAA1Mediated Auxin Biosynthesis Is Essential for Hormone Crosstalk and Plant Development. Cell 133: 177–191. Sulpice, R., Flis, A., Ivakov, A. A., Apelt, F., Krohn, N., Encke, B., Abel, C., Feil, R., Lunn, J.E., and Stitt, M. (2014). Arabidopsis coordinates the diurnal regulation of carbon allocation and growth across a wide range of Photoperiods. Mol. Plant 7: 137–155. Sun, J., Qi, L., Li, Y., Chu, J., and Li, C. (2012). PIF4-mediated activation of YUCCA8 expression integrates temperature into the auxin pathway in regulating Arabidopsis hypocotyl growth. PLoS Genet. 8: e1002594. Tabas-Madrid, D., Nogales-Cadenas, R., and Pascual-Montano, A. (2012). GeneCodis 3: a non-redundant and modular enrichment analysis tool for functional genomics. Nucleic Acids Res. 40: W478–W483. Talon, M., Koornneef, M., and Zeevaart, J. A. D. (1990). Accumulation of C19gibberellins in the gibberellin-insensitive dwarf mutant gai of Arabidopsis thaliana (L.) Heynh. Planta 182: 501–505. 31 Downloaded from www.plantphysiol.org on June 22, 2016 - Published by www.plantphysiol.org Copyright © 2016 American Society of Plant Biologists. All rights reserved.

1092 1093 1094 1095 1096 1097 1098 1099 1100 1101 1102 1103 1104 1105 1106 1107 1108 1109 1110 1111 1112 1113 1114 1115 1116 1117 1118 1119 1120 1121 1122 1123 1124 1125 1126 1127 1128 1129 1130 1131 1132 1133 1134 1135 1136 1137 1138 1139 1140

Tanaka, S., Nakamura, S., Mochizuki, N., and Nagatani, A. (2002). Phytochrome in cotyledons regulates the expression of genes in the hypocotyl through auxindependent and -independent pathways. The Plant Cell Physiol. 43: 1171–1181. Tao, Y., Ferrer, J.L., Ljung, K., Pojer, F., Hong, F., Long, J.A., Li, L., Moreno, J.E., Bowman, M.E., Ivans, L.J., Cheng, Y., Lim, J., Zhao, Y., Ballaré, C.L., Sandberg, G., Noel, J.P. and Chory, J. (2008). Rapid Synthesis of Auxin via a New Tryptophan-Dependent Pathway Is Required for Shade Avoidance in Plants. Cell 133: 164–176. Tholen, D., Pons, T.L., Voesenek, L.A.C.J., and Poorter, H. (2007). Ethylene insensitivity results in down-regulation of rubisco expression and photosynthetic capacity in tobacco. Plant Physiol. 144: 1305–1315. Turk, E.M. et al. (2005). BAS1 and SOB7 act redundantly to modulate Arabidopsis photomorphogenesis via unique brassinosteroid inactivation mechanisms. Plant J. 42: 23–34. Unterholzner, S.J., Rozhon, W., Papacek, M., Ciomas, J., Lange, T., Kugler, K.G., Mayer, K.F., Sieberer, T., and Poppenberger, B. (2015). Brassinosteroids Are Master Regulators of Gibberellin Biosynthesis in Arabidopsis. The Plant Cell 27: 2261–2272. Van Esse, G.W., van Mourik, S., Stigter, H., ten Hove, C. a., Molenaar, J., and de Vries, S.C. (2012). A Mathematical Model for BRASSINOSTEROID INSENSITIVE1-Mediated Signaling in Root Growth and Hypocotyl Elongation. Plant Physiol. 160: 523–532. Vandenbussche, F., Pierik, R., Millenaar, F.F., Voesenek, L.A.C.J., and Van Der Straeten, D. (2005). Reaching out of the shade. Curr. Opin. Plant Biol. 8: 462– 468. van Veen, H., Mustroph, A., Barding, G.A., Vergeer-van Eijk, M., WelschenEvertman, R.A.M., Pedersen, O., Visser, E.J.W., Larive, C.K., Pierik, R., Bailey-Serres, J., Voesenek, L.A.C.J., and Sasidharan, R. (2013). Two Rumex species from contrasting hydrological niches regulate flooding tolerance through distinct mechanisms. The Plant Cell 25: 4691–707. Voesenek L.A.C.J.; and Bailey-serres, J. (2015). Flood adaptive traits and processes : an overview. New Phytol. 206: 57–73. Volodarsky, D., Leviatan, N., Otcheretianski, A., and Fluhr, R. (2009). HORMONOMETER: a tool for discerning transcript signatures of hormone action in the Arabidopsis transcriptome. Plant Physiol. 150: 1796–1805. Warnasooriya, S.N. and Montgomery, B.L. (2009). Detection of spatial-specific phytochrome responses using targeted expression of biliverdin reductase in Arabidopsis. Plant Physiol. 149: 424–33. Wilson, R.N., Heckman, J.W., and Somerville, C.R. (1992). Gibberellin Is Required for Flowering in Arabidopsis thaliana under Short Days. Plant Physiol. 100: 403– 408. Zhiponova, M.K., Morohashi, K., Vanhoutte, I., Machemer-Noonan, K., Revalska, M., Van Montagu, M., Grotewold, E., and Russinova, E. (2014). Helix-loop-helix/basic helix-loop-helix transcription factor network represses cell elongation in Arabidopsis through an apparent incoherent feed-forward loop. Proc. Natl. Acad. Sci. U. S. A. 111: 2824–9. Zheng, Z., Guo, Y., Novak, O., Chen, W., Ljung, K., Noel, J.P., and Chory, J. (2016). Local auxin metabolism regulates environment-induced hypocotyl elongation. Nat. Plants. 21: 16025. doi: 10.1038/nplants.2016.25

32 Downloaded from www.plantphysiol.org on June 22, 2016 - Published by www.plantphysiol.org Copyright © 2016 American Society of Plant Biologists. All rights reserved.

1141 1142 1143 1144 1145 1146 1147 1148 1149

Zhong, S., Shi, H., Xue, C., Wang, L., Xi, Y., Li, J., Quail, P.H., Deng, X.W., and Guo, H. (2012). A molecular framework of light-controlled phytohormone action in Arabidopsis. Curr. Biol. 22: 1530–1535.

33 Downloaded from www.plantphysiol.org on June 22, 2016 - Published by www.plantphysiol.org Copyright © 2016 American Society of Plant Biologists. All rights reserved.

Downloaded from www.plantphysiol.org on June 22, 2016 - Published by www.plantphysiol.org Copyright © 2016 American Society of Plant Biologists. All rights reserved.

Downloaded from www.plantphysiol.org on June 22, 2016 - Published by www.plantphysiol.org Copyright © 2016 American Society of Plant Biologists. All rights reserved.

Downloaded from www.plantphysiol.org on June 22, 2016 - Published by www.plantphysiol.org Copyright © 2016 American Society of Plant Biologists. All rights reserved.

Downloaded from www.plantphysiol.org on June 22, 2016 - Published by www.plantphysiol.org Copyright © 2016 American Society of Plant Biologists. All rights reserved.

Downloaded from www.plantphysiol.org on June 22, 2016 - Published by www.plantphysiol.org Copyright © 2016 American Society of Plant Biologists. All rights reserved.

Downloaded from www.plantphysiol.org on June 22, 2016 - Published by www.plantphysiol.org Copyright © 2016 American Society of Plant Biologists. All rights reserved.

Downloaded from www.plantphysiol.org on June 22, 2016 - Published by www.plantphysiol.org Copyright © 2016 American Society of Plant Biologists. All rights reserved.

Downloaded from www.plantphysiol.org on June 22, 2016 - Published by www.plantphysiol.org Copyright © 2016 American Society of Plant Biologists. All rights reserved.

Parsed Citations Abdi, H., and Williams, L. J. (2010). Principal component analysis. Wiley Interdisciplinary Reviews: Computational Statistics, 2(4), 433-459. Pubmed: Author and Title CrossRef: Author and Title Google Scholar: Author Only Title Only Author and Title

Achard, P., Cheng, H., Grauwe, L. De, Decat, J., Schoutteten, H., Moritz, T., Straeten, D. Van Der, Peng, J., and Harberd, N.P. (2006). Integration of Plant Responses to Environmentally Activated Phytohormonal Signals. Science 311: 91-94. Pubmed: Author and Title CrossRef: Author and Title Google Scholar: Author Only Title Only Author and Title

Asami, T., Min, Y.K., Nagata, N., Yamagishi, K., Takatsuto, S., Fujioka, S., Murofushi, N., Yamaguchi, I., and Yoshida, S. (2000). Characterization of brassinazole, a triazole-type brassinosteroid biosynthesis inhibitor. Plant Physiol. 123: 93-100. Pubmed: Author and Title CrossRef: Author and Title Google Scholar: Author Only Title Only Author and Title

Azpiroz, R., Wu, Y., LoCascio, J.C., and Feldmann, K.A. (1998). An Arabidopsis brassinosteroid-dependent mutant is blocked in cell elongation. The Plant Cell 10: 219-230. Pubmed: Author and Title CrossRef: Author and Title Google Scholar: Author Only Title Only Author and Title

Bai, M.Y., Fan, M., Oh, E., and Wang, Z.Y. (2012) A triple helix-loop-helix/basic helix-loop-helix cascade controls cell elongation downstream of multiple hormonal and environmental signaling pathways in Arabidopsis. Plant Cell. 24: 4917-4929 Pubmed: Author and Title CrossRef: Author and Title Google Scholar: Author Only Title Only Author and Title

Barboza, L., Effgen, S., Alonso-blanco, C., Kooke, R., Keurentjes, J.J.B., and Koornneef, M. (2013). Arabidopsis semidwarfs evolved from independent mutations in GA20ox1, ortholog to green revolution dwarf alleles in rice and barley. Proc. Natl. Acad. Sci. U. S. A. 110: 15818-15823. Pubmed: Author and Title CrossRef: Author and Title Google Scholar: Author Only Title Only Author and Title

Binder, B.M., Mortimore, L.A., Stepanova, A.N., Ecker, J.R., and Bleecker, A.B. (2004). Short-Term Growth Responses to Ethylene in Arabidopsis. Plant Physiol. 136: 2921-2927. Pubmed: Author and Title CrossRef: Author and Title Google Scholar: Author Only Title Only Author and Title

Blilou, I., Xu, J., Wildwater, M., Willemsen, V., Paponov, I., Friml J, Heidstra R., Aida M., Palme K. & Scheres, B. (2005). The PIN auxin efflux facilitator network controls growth and patterning in Arabidopsis roots. Nature, 433(7021), 39-44. Pubmed: Author and Title CrossRef: Author and Title Google Scholar: Author Only Title Only Author and Title

Cagnola, J. I., Ploschuk, E., Benech-Arnold, T., Finlayson, S. A., and Casal, J. J. (2012). Stem transcriptome reveals mechanisms to reduce the energetic cost of shade-avoidance responses in tomato. Plant Physiol., 160(2), 1110-1119. Pubmed: Author and Title CrossRef: Author and Title Google Scholar: Author Only Title Only Author and Title

Capella, M., Ribone, P. A., Arce, A. L., and Chan, R. L. (2015). Arabidopsis thaliana HomeoBox 1 (AtHB1), a Homedomain-Leucine Zipper I (HD-Zip I) transcription factor, is regulated by PHYTOCHROME-INTERACTING FACTOR 1 to promote hypocotyl elongation. New Phytologist, 207, 669-682. Pubmed: Author and Title CrossRef: Author and Title Google Scholar: Author Only Title Only Author and Title

Casal, J.J. (2013). Photoreceptor signaling networks in plant responses to shade. Annu. Rev. Plant Biol. 64: 403-27. Pubmed: Author and Title CrossRef: Author and Title Google Scholar: Author Only Title Only Author and Title

Casal, J.J. (2012). Shade avoidance. Arabidopsis Book 10: e0157. Pubmed: Author and Title CrossRef: Author and Title Google Scholar: Author Only Title Only Author and Title

Chae, K., Isaacs, C.G., Reeves, P.H., Maloney, G.S., Muday, G.K., Nagpal, P., and Reed, J.W. (2012). Arabidopsis SMALL AUXIN UP RNA63 promotes hypocotyl and stamen filament elongation. Plant J. 71: 684-697. Pubmed: Author and Title CrossRef: Author and Title Google Scholar: Author Only Title Only Author and Title Downloaded from www.plantphysiol.org on June 22, 2016 - Published by www.plantphysiol.org Copyright © 2016 American Society of Plant Biologists. All rights reserved.

Chapman, E.J., Greenham, K., Castillejo, C., Sartor, R., Bialy, A., Sun, T.P., and Estelle, M. (2012). Hypocotyl transcriptome reveals auxin regulation of growth-promoting genes through GA-dependent and -independent pathways. PLoS One 7: e36210. Pubmed: Author and Title CrossRef: Author and Title Google Scholar: Author Only Title Only Author and Title

De Smet, I., Lau, S., Ehrismann, J.S., Axiotis, I., Kolb, M., Kientz, M., Weijers, D., and Jurgens, G. (2013). Transcriptional repression of BODENLOS by HD-ZIP transcription factor HB5 in Arabidopsis thaliana. J. Exp. Bot. 64: 3009-3019. Pubmed: Author and Title CrossRef: Author and Title Google Scholar: Author Only Title Only Author and Title

Djakovic-Petrovic, T., Wit, M. De, Voesenek, L.A.C.J., and Pierik, R. (2007). DELLA protein function in growth responses to canopy signals. Plant J. 51: 117-126. Pubmed: Author and Title CrossRef: Author and Title Google Scholar: Author Only Title Only Author and Title

Eisen, M.B., Spellman, P.T., Brown, P.O., and Botstein, D. (1998). Cluster analysis and display of genome-wide expression patterns. Proc. Natl. Acad. Sci. U. S. A. 95: 14863-14868. Pubmed: Author and Title CrossRef: Author and Title Google Scholar: Author Only Title Only Author and Title

Endo, M., Nakamura, S., Araki, T., Mochizuki, N., and Nagatani, A. (2005). Phytochrome B in the Mesophyll Delays Flowering by Suppressing FLOWERING LOCUS T Expression in Arabidopsis Vascular Bundles. The Plant Cell 17: 1941-1952. Pubmed: Author and Title CrossRef: Author and Title Google Scholar: Author Only Title Only Author and Title

Estelle, M. (1998). Polar Auxin Transport: New Support for an Old Model. The Plant Cell 10: 1775-1778. Pubmed: Author and Title CrossRef: Author and Title Google Scholar: Author Only Title Only Author and Title

Fan, M., Bai, M. Y., Kim, J. G., Wang, T., Oh, E., Chen, L., Park, C.H., Son S., Kim S., Mudgett M.B. and Wang, Z. Y. (2014). The bHLH transcription factor HBI1 mediates the trade-off between growth and pathogen-associated molecular pattern-triggered immunity in Arabidopsis. The Plant Cell Online, 26(2), 828-841. Pubmed: Author and Title CrossRef: Author and Title Google Scholar: Author Only Title Only Author and Title

Franklin, K. A., Lee, S.H., Patel, D., Kumar, S.V., Spartz, A.K., Gu, C., Ye, S., Yu, P., Breen, G., Cohen, J.D., Wigge, P.A., and Gray, W.M. (2011). Phytochrome-interacting factor 4 (PIF4) regulates auxin biosynthesis at high temperature. Proc. Natl. Acad. Sci. U. S. A. 108: 20231-5. Pubmed: Author and Title CrossRef: Author and Title Google Scholar: Author Only Title Only Author and Title

Franklin, K. A. (2008). Shade avoidance. New Phytol. 179: 930-944. Pubmed: Author and Title CrossRef: Author and Title Google Scholar: Author Only Title Only Author and Title

Gendreau, E., Traas, J., Desnos, T., Grandjean, O., Caboche, M., and Höfte, H. (1997). Cellular basis of hypocotyl growth in Arabidopsis thaliana. Plant Physiol. 114: 295-305. Pubmed: Author and Title CrossRef: Author and Title Google Scholar: Author Only Title Only Author and Title

Hattori, Y., Nagai, K., Furukawa, S., Song, X., Kawano, R., Sakakibara, H., Wu, J., Matsumoto, T., Yoshimura, A., Kitano, H., Matsuoka, M., Mori, H., Ashikari, M. (2009). The ethylene response factors SNORKEL1 and SNORKEL2 allow rice to adapt to deep water. Nature 460: 1026-1030. Pubmed: Author and Title CrossRef: Author and Title Google Scholar: Author Only Title Only Author and Title

Hayashi, K., Tan, X., Zheng, N., Hatate, T., Kimura, Y., Kepinski, S., and Nozaki, H. (2008). Small-molecule agonists and antagonists of F-box protein-substrate interactions in auxin perception and signaling. Proc. Natl. Acad. Sci. U. S. A. 105: 5632-5637. Pubmed: Author and Title CrossRef: Author and Title Google Scholar: Author Only Title Only Author and Title

Henriques, R., Bögre, L., Horváth, B., and Magyar, Z. (2014). Balancing act: Matching growth with environment by the TOR signalling pathway. J. Exp. Bot. 65: 2691-2701. Pubmed: Author and Title CrossRef: Author and Title Google Scholar: Author Only Title Only Author and Title

Hisamatsu, T., King, R.W., Helliwell, C. a, and Koshioka, M. (2005). The involvement of Gibberellin 20-Oxidase genes in Downloaded from www.plantphysiol.org on June 22, 2016 - Published by www.plantphysiol.org Copyright © 2016 American Society of Plant Biologists. All rights reserved.

phytochrome-regulated petiole elongation of Arabidopsis. Plant Physiol. 138: 1106-1116. Pubmed: Author and Title CrossRef: Author and Title Google Scholar: Author Only Title Only Author and Title

Hong, S.Y., Kim, O.K., Kim, S.G., Yang, M.S., and Park, C.M. (2011). Nuclear import and DNA binding of the ZHD5 transcription factor is modulated by a competitive peptide inhibitor in Arabidopsis. J. Biol. Chem. 286: 1659-1668. Pubmed: Author and Title CrossRef: Author and Title Google Scholar: Author Only Title Only Author and Title

Hu, W. and Ma, H. (2006). Characterization of a novel putative zinc finger gene MIF1: Involvement in multiple hormonal regulation of Arabidopsis development. Plant J. 45: 399-422. Pubmed: Author and Title CrossRef: Author and Title Google Scholar: Author Only Title Only Author and Title

Hyun, Y. and Lee, I. (2006). KIDARI, encoding a non-DNA Binding bHLH protein, represses light signal transduction in Arabidopsis thaliana. Plant Mol. Biol. 61: 283-96. Pubmed: Author and Title CrossRef: Author and Title Google Scholar: Author Only Title Only Author and Title

Iqbal, N., Khan, N. A., Nazar, R., and da Silva, J. A. T. (2012). Ethylene-stimulated photosynthesis results from increased nitrogen and sulfur assimilation in mustard types that differ in photosynthetic capacity. Environ. Exp. Bot. 78: 84-90. Pubmed: Author and Title CrossRef: Author and Title Google Scholar: Author Only Title Only Author and Title

Kazan, K. and Manners, J.M. (2013). MYC2: The master in action. Mol. Plant 6: 686-703. Pubmed: Author and Title CrossRef: Author and Title Google Scholar: Author Only Title Only Author and Title

Keuskamp, D.H., Pollmann, S., Voesenek, L.A.C.J., Peeters, A.J.M., and Pierik, R. (2010). Auxin transport through PIN-FORMED 3 (PIN3) controls shade avoidance and fitness during competition. Proc. Natl. Acad. Sci. U. S. A. 107: 22740-22744. Pubmed: Author and Title CrossRef: Author and Title Google Scholar: Author Only Title Only Author and Title

Keuskamp, D.H., Sasidharan, R., Vos, I., Peeters, A.J.M., Voesenek, L.A.C.J., and Pierik, R. (2011). Blue-light-mediated shade avoidance requires combined auxin and brassinosteroid action in Arabidopsis seedlings. Plant J. 67: 208-217. Pubmed: Author and Title CrossRef: Author and Title Google Scholar: Author Only Title Only Author and Title

Khanna, R., Shen, Y., Marion, C.M., Tsuchisaka, A., Theologis, A., Schafer, E., and Quail, P.H. (2007). The basic helix-loop-helix transcription factor PIF5 acts on ethylene biosynthesis and phytochrome signaling by distinct mechanisms. Plant Cell. 19: 39153929. Pubmed: Author and Title CrossRef: Author and Title Google Scholar: Author Only Title Only Author and Title

Kim, B.C., Soh, M.S., Hong, S.H., Furuya, M., and Nam, H.G. (1998). Photomorphogenic development of the Arabidopsis shy2-1D mutation and its interaction with phytochromes in darkness. Plant J. 15: 61-68. Pubmed: Author and Title CrossRef: Author and Title Google Scholar: Author Only Title Only Author and Title

Kunihiro, A., Yamashino, T., Nakamichi, N., Niwa, Y., Nakanishi, H., and Mizuno, T. (2011). PHYTOCHROME-INTERACTING FACTOR 4 and 5 (PIF4 and PIF5) activate the Homeobox ATHB2 and auxin-inducible IAA29 genes in the coincidence mechanism underlying photoperiodic control of plant growth of Arabidopsis thaliana. The Plant Cell Physiol. 52: 1315-1329. Pubmed: Author and Title CrossRef: Author and Title Google Scholar: Author Only Title Only Author and Title

Lee, Y.K., Kim, G., Kim, I., Park, J., Kwak, S., Choi, G., and Chung, W. (2006). LONGIFOLIA1 and LONGIFOLIA2, two homologous genes, regulate longitudinal cell elongation in Arabidopsis. Development 133: 4305-4314. Pubmed: Author and Title CrossRef: Author and Title Google Scholar: Author Only Title Only Author and Title

Leivar, P. and Monte, E. (2014). PIFs: systems integrators in plant development. The Plant Cell 26: 56-78. Pubmed: Author and Title CrossRef: Author and Title Google Scholar: Author Only Title Only Author and Title

Li, L. et al. (2012). Linking photoreceptor excitation to changes in plant architecture. Genes Dev. 26: 785-790. Pubmed: Author and Title CrossRef: Author and Title Downloaded from on June 22, 2016 - Published by www.plantphysiol.org Google Scholar: Author Only Title Only Author andwww.plantphysiol.org Title Copyright © 2016 American Society of Plant Biologists. All rights reserved.

Lilley, J.L.S., Gee, C.W., Sairanen, I., Ljung, K., and Nemhauser, J.L. (2012). An Endogenous Carbon-Sensing Pathway Triggers Increased Auxin Flux and Hypocotyl Elongation. Plant Physiol. 160: 2261-2270. Pubmed: Author and Title CrossRef: Author and Title Google Scholar: Author Only Title Only Author and Title

Mickelbart, M. V., Hasegawa, P.M., and Bailey-Serres, J. (2015). Genetic mechanisms of abiotic stress tolerance that translate to crop yield stability. Nat. Rev. Genet. 16: 237-251. Pubmed: Author and Title CrossRef: Author and Title Google Scholar: Author Only Title Only Author and Title

Morelli, G., Ruberti, I., Molecolare, B., Sapienza, L., and Moro, P.A. (2000). Update on Light Signaling Shade Avoidance Responses. Driving Auxin along Lateral Routes. Plant Physiol. 122: 621-626. Pubmed: Author and Title CrossRef: Author and Title Google Scholar: Author Only Title Only Author and Title

Neff, M.M., Nguyen, S.M., Malancharuvil, E.J., Fujioka, S., Noguchi, T., Seto, H., Tsubuki, M., Honda, T., Takatsuto, S., Yoshida, S., and Chory, J. (1999). BAS1: A gene regulating brassinosteroid levels and light responsiveness in Arabidopsis. Proc. Natl. Acad. Sci. 96: 15316-15323. Pubmed: Author and Title CrossRef: Author and Title Google Scholar: Author Only Title Only Author and Title

Nemhauser, J.L., Mockler, T.C., and Chory, J. (2004). Interdependency of brassinosteroid and auxin signaling in Arabidopsis. PLoS Biol. 2: E258. Pubmed: Author and Title CrossRef: Author and Title Google Scholar: Author Only Title Only Author and Title

Nishimura, T., Hayashi, K.I., Suzuki, H., Gyohda, A., Takaoka, C., Sakaguchi, Y., Matsumoto, S., Kasahara, H., Sakai, T., Kato, J.I., Kamiya, Y., and Koshiba, T. (2014). Yucasin is a potent inhibitor of YUCCA, a key enzyme in auxin biosynthesis. Plant J. 77: 352-366. Pubmed: Author and Title CrossRef: Author and Title Google Scholar: Author Only Title Only Author and Title

Nomura, T., Kushiro, T., Yokota, T., Kamiya, Y., Bishop, G.J., and Yamaguchi, S. (2005). The last reaction producing brassinolide is catalyzed by cytochrome P-450s, CYP85A3 in tomato and CYP85A2 in Arabidopsis. J. Biol. Chem. 280: 17873-17879. Pubmed: Author and Title CrossRef: Author and Title Google Scholar: Author Only Title Only Author and Title

Nozue, K., Tat, A. V., Kumar Devisetty, U., Robinson, M., Mumbach, M.R., Ichihashi, Y., Lekkala, S., and Maloof, J.N. (2015). Shade Avoidance Components and Pathways in Adult Plants Revealed by Phenotypic Profiling. PLOS Genet. 11: e1004953. Pubmed: Author and Title CrossRef: Author and Title Google Scholar: Author Only Title Only Author and Title

Osakabe, Y., Osakabe, K., Shinozaki, K., and Tran, L. S. P. (2014). Response of plants to water stress. Front. Plant Sci. 5(86): 103389. Pubmed: Author and Title CrossRef: Author and Title Google Scholar: Author Only Title Only Author and Title

Park, E., Park, J., Kim, J., Nagatani, A., Lagarias, J.C., and Choi, G. (2012). Phytochrome B inhibits binding of phytochromeinteracting factors to their target promoters. Plant J. 72: 537-546. Pubmed: Author and Title CrossRef: Author and Title Google Scholar: Author Only Title Only Author and Title

Petrásek, J., Cerná, A., Schwarzerová, K., Elckner, M., Morris, D. A., and Zazímalová, E. (2003). Do phytotropins inhibit auxin efflux by impairing vesicle traffic? Plant Physiol. 131: 254-263. Pubmed: Author and Title CrossRef: Author and Title Google Scholar: Author Only Title Only Author and Title

Pfeiffer, A., Shi, H., Tepperman, J.M., Zhang, Y., and Quail, P.H. (2014). Combinatorial Complexity in a Transcriptionally Centered Signaling Hub in Arabidopsis. Mol. Plant 7: 1598-1618. Pubmed: Author and Title CrossRef: Author and Title Google Scholar: Author Only Title Only Author and Title

Philippar, K., Ivashikina, N., Ache, P., Christian, M., Lüthen, H., Palme, K., and Hedrich, R. (2004). Auxin activates KAT1 and KAT2, two K+-channel genes expressed in seedlings of Arabidopsis thaliana. Plant J. 37: 815-827. Pubmed: Author and Title CrossRef: Author and Title Google Scholar: Author Only Title Only Author and Title Downloaded from www.plantphysiol.org on June 22, 2016 - Published by www.plantphysiol.org Copyright © 2016 American Society of Plant Biologists. All rights reserved.

Pierik, R., Cuppens, M.L., Voesenek, L.A.C.J. and Visser, E.J. (2004) Interactions between ethylene and gibberellins in phytochrome-mediated shade avoidance responses in tobacco. Plant Physiol. 136: 2928-2936. Pubmed: Author and Title CrossRef: Author and Title Google Scholar: Author Only Title Only Author and Title

Pierik, R., Djakovic-Petrovic, T., Keuskamp, D.H., de Wit, M., and Voesenek, L.A.C.J. (2009). Auxin and ethylene regulate elongation responses to neighbor proximity signals independent of gibberellin and DELLA proteins in Arabidopsis. Plant Physiol. 149: 17011712. Pubmed: Author and Title CrossRef: Author and Title Google Scholar: Author Only Title Only Author and Title

Pierik, R., Millenaar, F.F., Peeters, A. J.M., and Voesenek, L.A.C.J. (2005). New perspectives in flooding research: The use of shade avoidance and Arabidopsis thaliana. Ann. Bot. 96: 533-540. Pubmed: Author and Title CrossRef: Author and Title Google Scholar: Author Only Title Only Author and Title

Pierik, R. and Testerink, C. (2014). The art of being flexible: how to escape from shade, salt and drought. Plant Physiol. 166: 5-22. Pubmed: Author and Title CrossRef: Author and Title Google Scholar: Author Only Title Only Author and Title

Pierik, R. and De Wit, M. (2014). Shade avoidance: Phytochrome signalling and other aboveground neighbour detection cues. J. Exp. Bot. 65: 2815-2824. Pubmed: Author and Title CrossRef: Author and Title Google Scholar: Author Only Title Only Author and Title

Procko, C., Crenshaw, C.M., Ljung, K., Noel, J.P., and Chory, J. (2014). Cotyledon-Generated Auxin Is Required for Shade-Induced Hypocotyl Growth in Brassica rapa. Plant Physiol. 165: 1285-1301. Pubmed: Author and Title CrossRef: Author and Title Google Scholar: Author Only Title Only Author and Title

Quint, M., Delker, C., Franklin, K. A., Wigge, P. A., Halliday, K. J., and van Zanten, M. (2016). Molecular and genetic control of plant thermomorphogenesis. Nature Plants. 2: 15190. Pubmed: Author and Title CrossRef: Author and Title Google Scholar: Author Only Title Only Author and Title

Rademacher, W. (2000). GROWTH RETARDANTS: Effects on Gibberellin Biosynthesis and Other Metabolic Pathways. Annu. Rev. Plant Physiol. Plant Mol. Biol. 51: 501-531. Pubmed: Author and Title CrossRef: Author and Title Google Scholar: Author Only Title Only Author and Title

Rieu, I., Ruiz-Rivero, O., Fernandez-Garcia, N., Griffiths, J., Powers, S.J., Gong, F., Linhartova, T., Eriksson, S., Nilsson, O., Thomas, S.G., Phillips, A.L., and Hedden, P. (2008). The gibberellin biosynthetic genes AtGA20ox1 and AtGA20ox2 act, partially redundantly, to promote growth and development throughout the Arabidopsis life cycle. Plant J. 53: 488-504. Pubmed: Author and Title CrossRef: Author and Title Google Scholar: Author Only Title Only Author and Title

Ruegger, M., Dewey, E., Gray, W.M., Hobbie, L., Turner, J., and Estelle, M. (1998). The TIR1 protein of Arabidopsis functions in auxin response and is related to human SKP2 and yeast Grr1p. Genes Dev. 12: 198-207. Pubmed: Author and Title CrossRef: Author and Title Google Scholar: Author Only Title Only Author and Title

Sasidharan, R., Chinnappa, C.C., Staal, M., Elzenga, J.T.M., Yokoyama, R., Nishitani, K., Voesenek, L.A.C.J., and Pierik, R. (2010). Light quality-mediated petiole elongation in Arabidopsis during shade avoidance involves cell wall modification by xyloglucan endotransglucosylase/hydrolases. Plant Physiol. 154: 978-990. Pubmed: Author and Title CrossRef: Author and Title Google Scholar: Author Only Title Only Author and Title

Sasidharan, R. and Voesenek, L.A.C.J. (2015). Ethylene-mediated acclimations to flooding stress. Plant Physiol.: pp.00387. Pubmed: Author and Title CrossRef: Author and Title Google Scholar: Author Only Title Only Author and Title

Smalle, J., Haegman, M., Kurepa, J., Van Montagu M, and Straeten, D. V (1997). Ethylene can stimulate Arabidopsis hypocotyl elongation in the light. Proc. Natl. Acad. Sci. U. S. A. 94: 2756-2761. Pubmed: Author and Title CrossRef: Author and Title Google Scholar: Author Only Title Only Author and Title

Spartz, A.K., Lee, S.H., Wenger, J.P., Gonzalez, N., Itoh, H., Inzé, D., Peer, W.A., Murphy, A.S., Overvoorde, P.J., and Gray, W.M. Downloaded from www.plantphysiol.org on June 22, 2016 - Published by www.plantphysiol.org Copyright © 2016 American Society of Plant Biologists. All rights reserved.

(2012). The SAUR19 subfamily of SMALL AUXIN UP RNA genes promote cell expansion. Plant J. 70: 978-990. Pubmed: Author and Title CrossRef: Author and Title Google Scholar: Author Only Title Only Author and Title

Stepanova, A.N., Robertson-Hoyt, J., Yun, J., Benavente, L.M., Xie, D.Y., Doležal, K., Schlereth, A., Jürgens, G., and Alonso, J.M. (2008). TAA1-Mediated Auxin Biosynthesis Is Essential for Hormone Crosstalk and Plant Development. Cell 133: 177-191. Pubmed: Author and Title CrossRef: Author and Title Google Scholar: Author Only Title Only Author and Title

Sulpice, R., Flis, A., Ivakov, A. A., Apelt, F., Krohn, N., Encke, B., Abel, C., Feil, R., Lunn, J.E., and Stitt, M. (2014). Arabidopsis coordinates the diurnal regulation of carbon allocation and growth across a wide range of Photoperiods. Mol. Plant 7: 137-155. Pubmed: Author and Title CrossRef: Author and Title Google Scholar: Author Only Title Only Author and Title

Sun, J., Qi, L., Li, Y., Chu, J., and Li, C. (2012). PIF4-mediated activation of YUCCA8 expression integrates temperature into the auxin pathway in regulating Arabidopsis hypocotyl growth. PLoS Genet. 8: e1002594. Pubmed: Author and Title CrossRef: Author and Title Google Scholar: Author Only Title Only Author and Title

Tabas-Madrid, D., Nogales-Cadenas, R., and Pascual-Montano, A. (2012). GeneCodis 3: a non-redundant and modular enrichment analysis tool for functional genomics. Nucleic Acids Res. 40: W478-W483. Pubmed: Author and Title CrossRef: Author and Title Google Scholar: Author Only Title Only Author and Title

Talon, M., Koornneef, M., and Zeevaart, J. A. D. (1990). Accumulation of C19-gibberellins in the gibberellin-insensitive dwarf mutant gai of Arabidopsis thaliana (L.) Heynh. Planta 182: 501-505. Pubmed: Author and Title CrossRef: Author and Title Google Scholar: Author Only Title Only Author and Title

Tanaka, S., Nakamura, S., Mochizuki, N., and Nagatani, A. (2002). Phytochrome in cotyledons regulates the expression of genes in the hypocotyl through auxin-dependent and -independent pathways. The Plant Cell Physiol. 43: 1171-1181. Pubmed: Author and Title CrossRef: Author and Title Google Scholar: Author Only Title Only Author and Title

Tao, Y., Ferrer, J.L., Ljung, K., Pojer, F., Hong, F., Long, J.A., Li, L., Moreno, J.E., Bowman, M.E., Ivans, L.J., Cheng, Y., Lim, J., Zhao, Y., Ballaré, C.L., Sandberg, G., Noel, J.P. and Chory, J. (2008). Rapid Synthesis of Auxin via a New Tryptophan-Dependent Pathway Is Required for Shade Avoidance in Plants. Cell 133: 164-176. Pubmed: Author and Title CrossRef: Author and Title Google Scholar: Author Only Title Only Author and Title

Tholen, D., Pons, T.L., Voesenek, L.A.C.J., and Poorter, H. (2007). Ethylene insensitivity results in down-regulation of rubisco expression and photosynthetic capacity in tobacco. Plant Physiol. 144: 1305-1315. Pubmed: Author and Title CrossRef: Author and Title Google Scholar: Author Only Title Only Author and Title

Turk, E.M. et al. (2005). BAS1 and SOB7 act redundantly to modulate Arabidopsis photomorphogenesis via unique brassinosteroid inactivation mechanisms. Plant J. 42: 23-34. Pubmed: Author and Title CrossRef: Author and Title Google Scholar: Author Only Title Only Author and Title

Unterholzner, S.J., Rozhon, W., Papacek, M., Ciomas, J., Lange, T., Kugler, K.G., Mayer, K.F., Sieberer, T., and Poppenberger, B. (2015). Brassinosteroids Are Master Regulators of Gibberellin Biosynthesis in Arabidopsis. The Plant Cell 27: 2261-2272. Pubmed: Author and Title CrossRef: Author and Title Google Scholar: Author Only Title Only Author and Title

Van Esse, G.W., van Mourik, S., Stigter, H., ten Hove, C. a., Molenaar, J., and de Vries, S.C. (2012). A Mathematical Model for BRASSINOSTEROID INSENSITIVE1-Mediated Signaling in Root Growth and Hypocotyl Elongation. Plant Physiol. 160: 523-532. Pubmed: Author and Title CrossRef: Author and Title Google Scholar: Author Only Title Only Author and Title

Vandenbussche, F., Pierik, R., Millenaar, F.F., Voesenek, L.A.C.J., and Van Der Straeten, D. (2005). Reaching out of the shade. Curr. Opin. Plant Biol. 8: 462-468. Pubmed: Author and Title CrossRef: Author and Title Google Scholar: Author Only Title Only Author and Title

van Veen, H., Mustroph, A., Barding, G.A., Vergeer-van Eijk, M., Welschen-Evertman, R.A.M., Pedersen, O., Visser, E.J.W., Larive, C.K., Pierik, R., Bailey-Serres, J., Voesenek, L.A.C.J., and Sasidharan, R. (2013). Two Rumex species from contrasting hydrological Downloaded from www.plantphysiol.org on June 22, 2016 - Published by www.plantphysiol.org Copyright © 2016 American Society of Plant Biologists. All rights reserved.

niches regulate flooding tolerance through distinct mechanisms. The Plant Cell 25: 4691-707. Pubmed: Author and Title CrossRef: Author and Title Google Scholar: Author Only Title Only Author and Title

Voesenek L.A.C.J.; and Bailey-serres, J. (2015). Flood adaptive traits and processes?: an overview. New Phytol. 206: 57-73. Pubmed: Author and Title CrossRef: Author and Title Google Scholar: Author Only Title Only Author and Title

Volodarsky, D., Leviatan, N., Otcheretianski, A., and Fluhr, R. (2009). HORMONOMETER: a tool for discerning transcript signatures of hormone action in the Arabidopsis transcriptome. Plant Physiol. 150: 1796-1805. Pubmed: Author and Title CrossRef: Author and Title Google Scholar: Author Only Title Only Author and Title

Warnasooriya, S.N. and Montgomery, B.L. (2009). Detection of spatial-specific phytochrome responses using targeted expression of biliverdin reductase in Arabidopsis. Plant Physiol. 149: 424-33. Pubmed: Author and Title CrossRef: Author and Title Google Scholar: Author Only Title Only Author and Title

Wilson, R.N., Heckman, J.W., and Somerville, C.R. (1992). Gibberellin Is Required for Flowering in Arabidopsis thaliana under Short Days. Plant Physiol. 100: 403-408. Pubmed: Author and Title CrossRef: Author and Title Google Scholar: Author Only Title Only Author and Title

Zhiponova, M.K., Morohashi, K., Vanhoutte, I., Machemer-Noonan, K., Revalska, M., Van Montagu, M., Grotewold, E., and Russinova, E. (2014). Helix-loop-helix/basic helix-loop-helix transcription factor network represses cell elongation in Arabidopsis through an apparent incoherent feed-forward loop. Proc. Natl. Acad. Sci. U. S. A. 111: 2824-9. Pubmed: Author and Title CrossRef: Author and Title Google Scholar: Author Only Title Only Author and Title

Zheng, Z., Guo, Y., Novak, O., Chen, W., Ljung, K., Noel, J.P., and Chory, J. (2016). Local auxin metabolism regulates environmentinduced hypocotyl elongation. Nat. Plants. 21: 16025. doi: 10.1038/nplants.2016.25 Pubmed: Author and Title CrossRef: Author and Title Google Scholar: Author Only Title Only Author and Title

Zhong, S., Shi, H., Xue, C., Wang, L., Xi, Y., Li, J., Quail, P.H., Deng, X.W., and Guo, H. (2012). A molecular framework of lightcontrolled phytohormone action in Arabidopsis. Curr. Biol. 22: 1530-1535. Pubmed: Author and Title CrossRef: Author and Title Google Scholar: Author Only Title Only Author and Title

Downloaded from www.plantphysiol.org on June 22, 2016 - Published by www.plantphysiol.org Copyright © 2016 American Society of Plant Biologists. All rights reserved.

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19

SUPPLEMENTAL METHODS Functional Enrichment Analysis g:Profiler (http://biit.cs.ut.ee/gprofiler/) webtool was used for finding enriched GO terms in the DEGs at 1.5 hrs with Best per parent (moderate filtering) setting from g:Profiler for “Hierarchical filtering” option in order to select only the most statistically significant term from a particular hierarchical group. A score above 1.3 for negative log of adjusted p-values was considered significant. Motif Enrichment Analysis Webtool (https://www.arabidopsis.org/tools/bulk/motiffinder/index.jsp) was used for identifying motif overrepresentation (using hypergeometric test) within the 1000 bp upstream sequences (relative to Transcription Start Site) of genes. MapMan Photosynthesis Genes Mappings excel sheet “Ath_AFFY_ATH1_TAIR10” at (http://mapman.gabipd.org/ web/guest/mapmanstore) was used to assign photosynthesis genes to the different categories in Supplemental Table S1.

20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 1

39 40 41 42

SUPPLEMENTAL FIGURES AND TABLES

43

44 45 46 47 48

Supplemental Figure S1. Effect of 7 d control, ethylene, shade and combination (ethylene + shade) treatments on hypocotyl length of Arabidopsis seedlings in short day photoperiod (9 h light / 15 h dark). Data represents mean ± S.E. (n=30-50). Different letters above the bars indicate significant differences (1way ANOVA followed by Post-Hoc pairwise comparison).

49 50 51 52

Supplemental Figure S2. g:Profiler analysis for enriched Gene ontology (GO) IDs (adj. p-value ≤ 0.05) or genes for DEGs (A) up- or (B) down- regulated at 1.5 hrs in the hypocotyl in response to ethylene or shade. Biological process (BP) terms corresponding to enriched GO IDs have been tabulated on the right side.

2

53 54 55 56 57

Supplemental Figure S3. Genotyping and transcript level verification. A, T-DNA insertion position in GABIKAT line GK-326G12 (athb28) B, Genotyping of athb28, (C) and (D) Null transcript and downregulation level of ATHB28 in athb28 (vs. Col-0). Abbreviations: In (B) LP: left primer, RP: right primer, LBo8474: GABI-KAT sequencing primer.

58

59 60 61 62 63

Supplemental Figure S4. Ethylene signalling mutant ein3eil1 shows hypocotyl elongation response similar to wild-type in response to shade but totally lacks elongation in response to ethylene. Hypocotyl length was measured following 96 h of ethylene and shade. Mean ± S.E. was calculated for 60 seedlings. Different letters indicate significant differences from a 2-way ANOVA followed by Tukey’s HSD Post-Hoc pairwise comparison.

64 65 66 67 68 69 70

Supplemental Figure S5. TAIR motif analysis and functional validation of PIFs and MYCs in hypocotyl elongation. TAIR motif overrepresentation (Lamesch et al., 2012). Top 2 hexamer motifs are shown for each set along with the name of motifs (relevance). % Microarray and % Genome indicate percentage of genes with the motif relative to the total number of genes in the gene set and in the genome respectively.

3

71 72 73 74

Supplemental Figure S6. Hypocotyl length of Landsberg erecta ecotype of Arabidopsis thaliana under control and 1 µlL-1 ethylene conditions when grown at 22⁰C day / 20⁰C night and 20⁰C day / 18⁰C night temperature regime.

75 76

77 78 79 80

Supplemental Table S1. Proportion of photosynthesis genes in different gene sets from the Venn diagram of negative clusters. Numbers in the table denote the % of total Mapman-assigned genes for the different photosynthetic processes. Bold numbers indicate top 3 photosynthetic categories.

81 82

LITERATURE CITED

83 84 85

J. Reimand, M. Kull, H. Peterson, J. Hansen, J. Vilo (2007). g:Profiler - a web-based toolset for functional profiling of gene lists from large-scale experiments NAR. 35: W193-W200

86 87 88 89 90

Nagel, A., Thimm, O., Redestig, H., Blaesing, O.E., Steinhauser, D., Gibon, Y., Morcuende, R., Weicht, D., Meyer, S., and Stitt, M. (2005). Extension of the Visualization Tool MapMan to Allow Statistical Analysis of Arrays, Display of Corresponding Genes, and Comparison with Known Responses. Plant Physiol. 138: 1195–1204.

91 92 93

4