(ESV1) is required for the control of starch ...

7 downloads 0 Views 4MB Size Report
May 20, 2016 - Mario Coirob, Sebastian Soykb4, Simona Eickeb, Tabea Mettler-Altmannb5, Kuan Jen. 5. Lub, Martin Tricka, Samuel C. Zeemanb, Alison M.
Plant Cell Advance Publication. Published on May 20, 2016, doi:10.1105/tpc.16.00011

1

The starch granule-associated protein EARLY STARVATION1 (ESV1) is required

2

for the control of starch degradation in Arabidopsis thaliana leaves

3 4

Doreen Feikea1, David Seungb1, Alexander Grafa2, Sylvain Bischofb3, Tamaryn Ellickb,

5

Mario Coirob, Sebastian Soykb4, Simona Eickeb, Tabea Mettler-Altmannb5, Kuan Jen

6

Lub, Martin Tricka, Samuel C. Zeemanb, Alison M. Smitha6

7 8 9 10

a

John Innes Centre, Norwich Research Park, Norwich NR4 7UH, United Kingdom

b

Department of Biology, Eidgenössische Technische Hochschule Zürich,

Universitätstrasse 2, CH-8092 Zürich, Switzerland

11 12 13 14 15 16 17 18 19 20 21

1

These authors contributed equally to the work

2

Current address: Department of Molecular Plant Physiology, Max Planck Institute for

Molecular Plant Physiology, D-14476 Potsdam, Germany 3

Current address: Department of Molecular, Cell and Developmental Biology, University

of California at Los Angeles, Los Angeles, CA 90095, USA 4

Current address: Cold Spring Harbor Laboratory, 1 Bungtown Road, Cold Spring

Harbor, NY 11724, USA 5

Current address: Institute of Plant Biochemistry, Cluster of Excellence on Plant

Sciences (CEPLAS), Heinrich-Heine-University, 40225 Düsseldorf, Germany 6

Address correspondence to [email protected].

22 23 24 25

Short title: Control of leaf starch breakdown

26 27 28

The author responsible for distribution of materials integral to the findings presented in

29

this article in accordance with the policy described in the Instructions for Authors

30

(www.plantcell.org) is: Alison Smith ([email protected]).

31 1 ©2016 American Society of Plant Biologists. All Rights Reserved.

32

SYNOPSIS

33

Two proteins present in leaf starch granules are important for the control of starch

34

turnover, allowing plants to match the depletion of starch reserves to the length of the

35

night.

36 37

ABSTRACT

38 39

To uncover components of the mechanism that adjusts the rate of leaf starch

40

degradation to the length of the night, we devised a screen for mutant Arabidopsis

41

thaliana plants in which starch reserves are prematurely exhausted. The mutation in

42

one such mutant, named early starvation 1 or esv1, eliminates a previously-

43

uncharacterized protein. Starch in mutant leaves is degraded rapidly and in a non-linear

44

fashion, so that reserves are exhausted 2 h prior to dawn. The ESV1 protein and a

45

similar uncharacterized A. thaliana protein (named Like ESV1, or LESV) are located in

46

the chloroplast stroma and also bound into starch granules. The region of highest

47

similarity between the two proteins contains a series of near-repeated motifs rich in

48

tryptophan. Both proteins are conserved throughout starch-synthesizing organisms,

49

from angiosperms and monocots to green algae. Analysis of transgenic plants lacking

50

or over-expressing ESV1 or LESV, and of double mutants lacking ESV1 and another

51

protein necessary for starch degradation, leads us to propose that these proteins

52

function in the organization of the starch granule matrix. We argue that their

53

misexpression affects starch degradation indirectly, by altering matrix organization and

54

thus accessibility of starch polymers to starch-degrading enzymes.

55 56

INTRODUCTION

57 58

Normal growth rates in many plants are dependent upon the accumulation of starch as

59

a product of photosynthesis during the day and its controlled utilisation as a carbon

60

source for growth during the night. In Arabidopsis thaliana leaves, up to half of the

61

photosynthetically assimilated carbon accumulates as starch in leaf chloroplasts. At 2

62

night, conversion of starch to sucrose proceeds at a near linear rate, such that about

63

95% of starch is consumed by dawn (Gibon et al., 2004; Gibon et al., 2009; Graf et al.,

64

2010; Sulpice et al., 2014). This pattern of starch biosynthesis and consumption is

65

adjusted in response to changing day length in a manner that optimises the allocation of

66

carbon for growth across a wide range of carbon availability. As day length decreases,

67

more photosynthate is allocated to starch in the day, and its rate of consumption at night

68

decreases so that supplies always last until dawn (e.g. Chatterton and Silvius, 1981;

69

Gibon et al., 2009; Sulpice et al., 2014). Adjustments are also made in response to

70

abrupt, unexpected alterations in day length (Lu et al., 2005). When plants grown in 12

71

h light – 12 h dark cycles are subjected to darkness after only 8 h of light, starch

72

degradation is slower than in a normal night so that reserves again last almost precisely

73

until dawn (Graf et al., 2010). Similarly, when the onset of darkness is delayed by 4 h,

74

starch degradation is faster so that reserves are again exhausted almost precisely at

75

dawn (Scialdone et al., 2013). Variation in the amount of starch present at the end of

76

the day also leads to adjustment of starch degradation rates to permit exhaustion of

77

reserves almost precisely at dawn (Graf et al., 2010, Scialdone et al., 2013).

78 79

We previously determined through experimental and chemical kinetic modelling

80

approaches that the adjustment of starch degradation to match the length of the night is

81

dependent on both the circadian clock and leaf starch content (Graf et al., 2010;

82

Scialdone et al., 2013). It is not known exactly how leaf starch content is perceived such

83

that it influences degradation rate. However, several potentially important features of

84

leaf starch granules are known to be under tight control. These include the timing of

85

granule initiation during leaf development, granule size, shape and number per

86

chloroplast (and therefore surface area), and the degree of crystallinity of the

87

constituent glucan polymers (Roldán et al., 2007; Szydlowski et al., 2009; Crumpton-

88

Taylor et al., 2012; Crumpton-Taylor et al., 2013; Pfister et al., 2014). These features

89

may combine with more specific control mechanisms that modulate the activities of the

90

enzymes of starch breakdown to confer the observed adjustment of the degradation

91

rate to match the length of the night.

3

92

To discover new components that influence starch degradation and mediate adjustment

93

of the rate according to the length of the night, we devised a forward genetic screen to

94

identify mutants defective in this adjustment. The screen made use of the observation

95

that if starch reserves are exhausted prior to dawn, plants rapidly exhibit large

96

transcriptional changes indicative of starvation. This phenomenon is seen in wild-type

97

plants when subjected to an extended night, beyond the point at which starch is

98

exhausted (Gibon et al., 2006), and in the short-period circadian clock mutant cca1 lhy

99

(circadian clock associated1 late elongated hypocotyl), which exhausts its starch

100

reserves about 3 h before dawn when grown in 12 h light – 12 h dark cycles (Graf et al.,

101

2010). Mutant plants unable to synthesize or to degrade starch have a transcriptional

102

signature of starvation throughout most of the night. We previously generated a

103

starvation reporter line of Arabidopsis harbouring a luciferase gene fused to the

104

promoter of a gene that is expressed only after starch reserves are exhausted.

105

Luciferase expression was very low during the normal day-night cycle but increased

106

rapidly after a 2-h extension of the night (Graf et al., 2010). We reasoned that

107

mutagenesis of the reporter line followed by screening for plants with abnormal

108

temporal patterns of luciferase expression should identify mutants defective in the diel

109

control of carbon availability, including aspects of starch granule biosynthesis and

110

degradation, elements of the circadian clock necessary for correct anticipation of dawn,

111

and components that sense or signal carbohydrate status. In the current study, we

112

validated this approach by showing that, as anticipated, the screen identified previously-

113

known classes of mutants defective in starch accumulation or with reduced starch

114

crystallinity.

115 116

The screen also identified a mutant with inappropriately rapid starch degradation at

117

night, such that starch was exhausted well before dawn. We named the mutant early

118

starvation1 (esv1-1). We established that the causal mutation was in an unannotated

119

gene, encoding a protein with no known or predicted role. The protein is highly

120

conserved and tryptophan-rich and is present in the chloroplast stroma and within

121

starch granules. Based on phenotypic characterization of mutant and over-expressing

122

lines, we propose that the protein has an unrecognized and central function in the 4

123

organization of glucan polymers within starch granules. We speculate that abnormal

124

levels of the protein alter the starch granule matrix such that the mechanisms controlling

125

the rate of starch degradation no longer function correctly. We show that another,

126

related protein is also likely to be involved in organization of glucan polymers within the

127

granule matrix and can indirectly influence the control of degradation. These results are

128

particularly surprising because it has been widely assumed that the formation of the

129

starch granule matrix is a physical process involving self-organization of glucan

130

polymers (Waigh et al., 1998; Waigh et al., 2000). Our discoveries suggest that specific

131

proteins are also required for this process.

132 133 134

RESULTS

135 136

The esv1 Mutant Has a Low Starch Content at the End of the Night

137 138

To identify plants defective in the control of carbon availability at night, we used a

139

forward genetic screen based on an Arabidopsis starvation reporter line (Graf et al.,

140

2010) in which luciferase expression is driven by the promoter of a sugar-repressed

141

gene (At1g10070). Seeds of the reporter line were treated with ethyl methanesulfonate

142

(EMS) to induce point mutations. Then, ten-day-old seedlings of the M2 generation

143

were screened for luciferase-induced bioluminescence at the end of the night, when

144

non-mutagenized seedlings of the reporter line showed no bioluminescence. Isolated

145

mutant plants were allowed to self-pollinate, and bioluminescence measurements were

146

repeated on M3 plants to check for reproducibility.

147 148

The efficacy of this screen was demonstrated by the identification of two anticipated

149

classes of mutants. Mutants that were shown by iodine staining to contain no starch at

150

the end of the day or the night were not examined further. Such mutants are expected

151

to exhibit symptoms of starvation during the night (Gibon et al., 2006). Mutants that

152

stained reddish-brown rather than blue-black with iodine at the end of the day proved to

153

be defective in ISOAMYLASE1 (ISA1) or ISA2 (Supplemental Figure 1). We showed 5

154

previously that isa1 and isa2 mutants accumulate soluble phytoglycogen (which stains

155

reddish-brown with iodine) in place of much of their starch during the day and

156

completely degrade this glucan before the end of the night (Zeeman et al., 1998; Delatte

157

et al., 2005).

158 159

One mutant line exhibited strong bioluminescence at the end of the night, but was

160

shown by iodine staining to have near-normal levels of starch at the end of the day

161

(Figure 1). To test whether the bioluminescence was a consequence of premature

162

exhaustion of starch during the night, rosettes were stained with iodine 2 h before the

163

end of the night. At this point, wild-type rosettes stained a pale color, indicating that

164

some starch remained. Mutant rosettes did not stain, indicating that much less starch

165

remained (Figure 1B). The mutant was named early starvation1, or esv1-1.

166 167

A Splice Site Mutation in At1g42430 Causes the esv1 Phenotype

168 169

The mutation underlying the esv1-1 phenotype was identified by a combination of map-

170

based cloning and whole-genome sequencing. The esv1-1 mutant (Col-0 background)

171

was crossed to a Landsberg erecta (Ler) wild-type plant. F1 plants were allowed to self-

172

pollinate. Seedlings of the F1 generation did not show bioluminescence at the end of

173

the night, but the F2 generation segregated for bioluminescence (Supplemental Figure

174

2A), indicating that the mutation in esv1-1 is recessive. Mapping narrowed the position

175

of the esv1 mutation to a 2.7 Mbp region spanning the centromere of chromosome 1

176

(Supplemental Figure 2, B and C; Supplemental Table 1).

177 178

We used whole genome sequencing to identify the ESV1 locus. DNA from 500 pooled

179

esv1-1 plants in the F2 generation of the mutant x Ler cross was used as a template. In

180

total, 9,000 putative EMS-induced point mutations were identified over the whole

181

genome, two of which were in genes within the target region. One caused an amino

182

acid substitution (glu910lys) in At1g42470, encoding a putative hedgehog receptor

183

located in the plasma membrane. The second was a G-to-A substitution that disrupted

184

the acceptor site at the second intron of At1g42430 (Figure 1C), a gene of unknown 6

185

function. Amplification of full-length cDNA from At1g42430 showed that the mutation

186

results in synthesis of a shorter transcript in esv1-1 than in wild-type plants (Figure 1D).

187

Thus, disruption of the intron splice site may result in splicing at an alternative site and

188

consequent deletions in the transcript. An antiserum was raised to the recombinant

189

protein product of the wild-type At1g42430 cDNA expressed in E. coli. It recognized a

190

protein of the expected mass on immunoblots of extracts of wild-type but not mutant

191

leaves (Figure 1E). The esv1-1 mutant was back-crossed three times to the starvation

192

reporter line before further characterization to minimize the number of unlinked EMS-

193

induced mutations.

194 195

To provide independent evidence that the mutation in At1g42430 was responsible for

196

the early starvation phenotype, we obtained a T-DNA insertion mutant from the GABI-

197

Kat (www.gabi-kat.de) collection (Figure 1C). Genotyping confirmed the presence of the

198

T-DNA insert in the fourth exon and immunoblot analysis confirmed that this insertion

199

mutant (named esv1-2) lacked the protein product (Figure 1E). Rosettes of esv1-2

200

showed no staining with iodine 2 h before the end of the night (Figure 1B), which is

201

consistent with the phenotype brought about by the original esv1-1 allele. To evaluate

202

starvation responses in the esv1-2 mutant, we measured transcript levels of two sugar-

203

repressed genes, At3g59940 and At1g76410. These genes encode an F-box protein

204

KMD4 and a C3HC4 type zinc-finger protein ATL8, respectively, and were used to

205

monitor starvation responses by Graf et al. (2010).

206 207

Like At1g10070, these genes are expressed at very low levels in wild-type plants at the

208

end of both the day and the night. Their expression rises only when sugar levels are

209

abnormally low, for example during an unexpected extension of the night in wild-type

210

plants, and during the normal night in starchless mutants (Bläsing et al., 2005; Osuna

211

et al., 2007; Usadel et al., 2008; Graf et al., 2010). In esv1-2 mutants, levels of

212

transcript of the two genes at the end of the day were similar to those in wild-type

213

plants, but levels at the end of the night were considerably elevated (Supplemental

214

Figure 3).

215 7

216

ESV1 is Similar to a Second Predicted Arabidopsis Protein; Both Are Conserved

217

Throughout the Plant Kingdom

218 219

The predicted ESV1 protein is composed of 426 amino acids, with a predicted mass of

220

about 49 kDa. It contains no previously-annotated domains but has a proline-rich region

221

at the C-terminal end (11 proline residues between amino acids 397 and 425) and is

222

enriched in tryptophan and other aromatic amino acid residues in the C-terminal two-

223

thirds of the protein (~11% of the amino acid residues between amino acids 130 and

224

380 are tryptophans). Within this region are motifs in which single or paired aromatic

225

amino acids are separated by two or three other amino acids, usually including one or

226

more acidic residues (E or D), for example WWETW, WTDKW, WEETWW,

227

WYEKWWEKY, WWEKWGEHY (Figure 2, Supplemental Figure 4).

228

A BLAST search (http://www.arabidopsis.org/Blast/index.jsp) showed that an

229

unannotated protein encoded by At3g55760 has the highest sequence similarity to

230

ESv1: the predicted amino acid sequences are 38% identical. We refer to this second

231

protein as LIKE ESV1, or LESV. It is composed of 578 amino acids with a predicted

232

mass of 66 kDa. LESV1 lacks the proline-rich C-terminal region of ESV1 but shares the

233

tryptophan-rich region with aromatic motifs similar to those in ESV1. Within this region,

234

the two proteins have similar numbers of tryptophan residues and the highest level of

235

overall identity (Figure 2A, Supplemental Figure 4A). Unlike ESV1, LESV is represented

236

on commonly-used Arabidopsis microarrays. Examination of publicly available

237

microarray data revealed that LESV is expressed in all organs of the plant. In leaves,

238

transcript levels are high at the end of the night and low during much of the day

239

(Supplemental Figure 5). LESV is strongly co-expressed with several genes encoding

240

enzymes of starch metabolism: the top twenty co-expressed genes from Atted-II

241

(www.atted.jp: Obayashi et al., 2009) encode enzymes including ISOAMYLASE3

242

(ISA3), α-GLUCAN PHOSPHORYLASE2 (PHS2), GLUCAN, WATER DIKINASE (GWD,

243

also known as STARCH EXCESS1, SEX1), STARCH-EXCESS 4 (SEX4),

244

DISPROPORTIONATING ENZYME1 (DPE1) and DPE2, STARCH BRANCHING

245

ENZYME3 (BE3), α-AMYLASE3 (AMY3) and ADP-GLUCOSE

8

246

PYROPHOSPHORYLASE LARGE SUBUNIT4 (APL4) [Supplemental Table 2; see

247

Smith (2012) for a discussion of the roles of these enzymes].

248 249

Genes homologous to ESV1 and LESV are found in land plants and green algae,

250

including microalgae in the genus Ostreococcus which have the smallest known

251

genomes among eukaryotes (Supplemental Figure 4B, Supplemental Table 3,

252

Supplemental Data Set 1). The genes are apparently absent from prokaryotes, red

253

algae, and other eukaryotic life-forms. ESV1 was included in a phylogenomic inventory

254

of proteins specific to the plant lineage (GreenCut2: Karpowicz et al., 2011). The high

255

level of evolutionary conservation of ESV1 and LESV protein sequences suggests that

256

both proteins have important functions that are specific to the Viridiplantae.

257 258

ESV1 and LESV are Chloroplastic and Bind to Starch Granules

259 260

Chloroplastic locations for ESV1 and LESV were suggested by previous publications on

261

the chloroplast proteome [e.g. Bayer et al. (2011) for ESV1; Kleffmann et al. (2004),

262

Peltier et al. (2006) for LESV]. The ChloroP program predicts a 56-amino acid N-

263

terminal chloroplast transit peptide (cTP) for LESV but no cTP for ESV1. However, the

264

predicted proteins most similar to ESV1 from Glycine max, Prunus persica, Vitis

265

vinifera, Manihot esculenta and Populus trichocarpa have putative transit peptides. We

266

sought to confirm the plastidial localisation of the Arabidopsis proteins by transiently

267

expressing them in leaves of woodland tobacco (Nicotiana sylvestris) as C-terminal

268

fusions to Yellow Fluorescent Protein (YFP). For both proteins, YFP fluorescence was

269

exclusively located in chloroplasts and was associated with discrete bodies likely to be

270

starch granules (Figure 2B). Stable expression of these constructs in Arabidopsis gave

271

the same result (Supplemental Figure 6A). To provide more information about the

272

location of the proteins within chloroplasts, we also expressed the fusion proteins

273

transiently in the pgm mutant of N. sylvestris, which lacks chloroplastic

274

phosphoglucomutase and therefore cannot synthesize starch (Hanson and McHale,

275

1988). In the pgm mutant, YFP appeared as a diffuse signal in the chloroplast stroma

276

for both fusion proteins rather than being associated with discrete structures. This 9

277

suggests that the fluorescent structures in chloroplasts of wild-type N. sylvestris are

278

indeed starch granules.

279 280

ESV1 and LESV were both partitioned between the soluble and insoluble (starch-

281

containing) fraction of leaves. Immunoblotting of extracts of wild-type N. sylvestris

282

leaves transiently expressing ESV1 or LESV as a YFP fusion revealed that the fusion

283

proteins were more abundant in the insoluble fraction – the expected location for starch-

284

bound proteins - than in the soluble fraction. Both proteins were mostly soluble when

285

transiently expressed in pgm leaves, which lack starch (Supplemental Figure 6B). In

286

wild-type Arabidopsis leaves, ESV1 and LESV proteins were also present in both

287

insoluble and soluble fractions (Figure 2C). The proteins were largely soluble at the end

288

of the night when starch levels were low, and they were largely insoluble at the end of

289

the day when starch level were maximal. ESV1 and LESV were present in proteins

290

extracted from purified starch from wild-type Arabidopsis leaves and were among the

291

most abundant granule-bound proteins in these leaves (Figure 2D). As the ESV1 and

292

LESV genes are conserved throughout the plant kingdom, we examined whether the

293

proteins are present in starches from species of economic significance. Both proteins

294

were present in starch purified from roots of cassava, tubers of potato, and mature

295

grains of rice and maize (Supplemental Table 4).

296 297

esv1 Mutants Cannot Adjust the Rate of Starch Degradation to the Length of the

298

Night

299 300

To characterize the defect in starch degradation in the esv1 mutants, we analyzed

301

starch contents in rosettes during a normal 12-h night and in rosettes subjected to

302

darkness after only 8 h of light (an early night) (Figure 3). Under both conditions, as

303

expected, wild-type plants degraded starch at a near linear rate and did not exhaust

304

their reserves before the end of the night. The rates of starch degradation were faster in

305

esv1 mutant plants than in wild-type plants, such that starch reserves were completely

306

exhausted 2 h before the end of the night in normal nights and 4 h before the end of the

307

night in early nights (Figure 3, A to C). Extrapolation of the near-linear phase of starch 10

308

degradation for esv1 plants, over the first 8 h of the night, revealed that these rates

309

(estimated from linear regression analysis) would exhaust starch reserves after

310

approximately 9 h in both normal (12-h) nights and early nights. Thus, wild-type plants

311

adjusted the rate of starch degradation to the length of the night, but esv1 mutant plants

312

failed to make this adjustment correctly.

313 314

To provide more information about the alterations in the control of starch degradation in

315

esv1 plants, we subjected batches of esv1 and wild-type plants to different light levels

316

over a single day to generate plants with very different starch levels at the start of the

317

12-h night. As expected, the rate of starch degradation in wild-type plants was adjusted

318

such that the starch reserves lasted for 12 h over a wide range of end-of-day starch

319

contents. In esv1 plants, the time at which starch reserves were exhausted depended

320

upon the initial starch content. Reserves were exhausted after 9 h of darkness in plants

321

with 7.5 mg starch g-1 fresh weight at the end of the day, and after 6 h of darkness in

322

plants with 2 mg starch g-1 fresh weight at the end of the day (Supplemental Figure 7, A

323

and B). Whereas in wild-type plants, the relative rate of starch degradation (fraction of

324

end-of-day starch degraded per unit time) was independent of starch content, in esv1

325

plants, the relative rate decreased as starch content decreased (Supplemental Figure

326

7C). Thus, in esv1, the dependence of starch degradation on time until dawn has

327

decreased or been lost.

328

In most experiments, esv1 mutants had lower rates of net starch accumulation than

329

wild-type plants during the day and hence, lower starch contents at the onset of

330

darkness (Figure 3). Interestingly, the levels of both sucrose and maltose in esv1-2

331

plants were higher than in wild-type plants during the day but lower than in wild-type

332

plants towards the end of the night (Figure 3D and E). Elevation of sucrose levels during

333

the day is commonly observed in Arabidopsis mutants with defects in pathways of

334

starch biosynthesis or degradation (e.g. Chia et al., 2004; Comparot-Moss et al., 2010;

335

Mugford et al., 2014). Elevated levels of maltose are generally indicative of starch

336

degradation (Weise et al., 2004 ). Taken together, these results indicate that esv1

337

mutants usually accumulate less starch than wild-type plants because they degrade

338

starch during the day as well as during the night. The low levels of sucrose and maltose 11

339

towards the end of night in esv1 are consistent with the observed premature starch

340

depletion in these plants.

341 342

We investigated whether loss of LESV also affects starch turnover by examining two T-

343

DNA insertion mutants (Figure 1C). The absence of LESV protein from these mutants

344

was confirmed by immunoblotting with a specific antiserum (Figure 1E). Diel patterns of

345

starch accumulation and loss in the lesv mutants were generally similar to those of wild-

346

type plants (Figure 3, A, F and G). Sucrose and maltose levels in the lesv-1 mutant

347

were also similar to those of wild-type plants over the day-night cycle (Figure 3, D and

348

E). Starch turnover and sucrose and maltose levels throughout the day-night cycle in an

349

esv1-2 lesv-1 double mutant resembled those of the esv1 parent (Figure 3, A, D and E).

350 351

Lack of ESV1 Affects Starch Contents in Non-Photosynthetic Tissues

352 353

We investigated whether ESV1 is important for the control of starch content in non-

354

photosynthetic parts of the plant (Figure 4). Detailed examination of iodine-stained

355

leaves showed that cells immediately adjacent to veins accumulate less starch during

356

the day in esv1 mutants than in wild-type plants. This netted pattern of starch

357

accumulation in esv1 leaves was particularly pronounced after an extended light period

358

of 24 h (Figure 4A). Staining for starch also revealed lower starch contents in esv1 than

359

in wild-type plants in columella cells of root caps, stems, flowers and siliques (Figure 4,

360

B to D).

361

Plants lacking ESV1 had several developmental and morphological phenotypes that

362

may be attributable to defective starch metabolism. After 24 d of growth, rosettes of

363

esv1 plants were smaller than those of wild-type plants (Figure 4E), and their fresh

364

weights were about 40% lower. Mutant plants also flowered later than wild-type plants.

365

Other mutants with defects in diel starch turnover in leaves are slow-growing in 12 h

366

light, 12 h dark cycles, e.g., the essentially starchless mutants pgm and adg1 (lacking

367

the small subunit of ADPglucose pyrophosphorylase) and the starch-degradation

368

mutants sex1 (lacking glucan, water dikinase) as well as bam1 bam3 (lacking β-

369

amylases 1 and 3) (Caspar et al., 1985; Lin et al., 1988; Yu et al., 2001; Fulton et al., 12

370

2008), and some exhibit delayed flowering, e.g., sex1 and pgi (lacking plastidial

371

phosphoglucomutase) (Eimert et al., 1995; Corbesier et al., 1998; Yu et al., 2000). The

372

axillary shoots of esv1 flower stems grew with a wider angle from the main stem than

373

those of wild-type plants (Figure 4F). Wide stem angles are characteristic of mutants

374

defective in the stem gravitropic response, which is dependent on the presence of

375

starch granules in endodermal cells (Weise and Kiss, 1999; Fujihira et al., 2000;

376

Tanimoto et al., 2008). Whereas the endodermis of wild-type and lesv stems contained

377

large starch granules, no starch was visible in this tissue in esv1 stems (Figure 4D).

378 379

Over-Expression of ESV1 and LESV Proteins Alters Leaf Starch Content

380 381

To provide more insight into the importance of ESV1 and LESV proteins in controlling

382

starch content in leaves, we generated plants with elevated amounts of these proteins

383

by constitutively expressing them as YFP fusions (Supplemental Fig 8A). Expression of

384

ESV1-YFP in either wild-type or esv1 mutant plants resulted in elevated starch

385

contents. In some lines, starch content was elevated at the end of the night but not at

386

the end of the day, whereas in the line with the highest level of expression of ESV1-YFP

387

(line 3-2), starch content was strongly elevated throughout the day-night cycle (Figure

388

3G).

389 390

We also generated plants in which ESV1 protein without a tag was expressed in an

391

esv1 background. These plants had smaller increases in ESV1 protein than those

392

expressing ESV1-YFP (Supplemental Figure 8B). Nonetheless, as for plants expressing

393

ESV1-YFP, starch contents of leaves were elevated relative to those of wild-type plants

394

at both the end of the day and the end of the night (Supplemental Figure 8C), showing

395

that this phenotype was not an artefact arising from the use of a YFP fusion protein.

396

Given this striking effect of ESV1 over-expression, we examined whether LESV over-

397

expression affected starch content. Constitutive expression of LESV-YFP in a wild-type

398

background had little impact on starch content at the end of the day, but starch content

399

10 h into the night and at the end of the night was about three-fold lower than in wild-

400

type plants (Figure 3G, Supplemental Figure 8A). These data suggested that LESV13

401

overexpressing plants might be starving at the end of the night. Consistent with this

402

expectation, the transcript abundance of two starvation genes was higher in LESV-

403

overexpressing plants than in wild-type plants at the end of the night, and was similar to

404

that in esv1 plants at this point (Supplemental Figure 3).

405 406

Altered Levels of ESV1 and LESV Affect Starch Granule Morphology, Number and

407

Composition but Have Only Minor Effects on Amylopectin Structure

408

Starch granules from leaves of esv1 and lesv mutants and the double mutant esv1 lesv

409

were roughly discoid in shape, like those of wild-type leaves. However, granules of the

410

mutants were less regular in shape than wild-type granules: they consistently had lobed

411

or wavy outlines (Figure 5). ESV1 over-expressing plants had thicker, larger granules

412

than wild-type plants. LESV over-expressing plants had many more granules per

413

chloroplast than wild-type plants. These granules were highly variable in size and

414

shape, and many were much smaller than granules of wild-type plants.

415 416

The starch of both esv1 and lesv mutants consistently contained about 60% more

417

amylose than that of wild-type plants (Figure 6). It seemed likely that this change in

418

starch composition was a secondary rather than the primary consequence of the

419

mutations, since the esv1 and lesv mutations have the same effect on amylose content

420

but radically different effects on the rate of starch degradation. To investigate this

421

further, we examined whether loss of amylose affects starch degradation in esv1

422

mutants by generating a double mutant lacking both ESV1 and granule-bound starch

423

synthase (GBSS), the enzyme responsible for amylose biosynthesis (Seung et al.,

424

2015). As expected, the esv1 gbss double mutant lacked amylose. Its starch content

425

was lower than that of wild-type or gbss plants and comparable to that of esv1 at both

426

the end of the day and two hours before the end of the night (Supplemental Figure 9A).

427

Thus, loss of ESV1 results in faster starch degradation regardless of whether amylose

428

is present in the starch.

429 430

Over-expression of ESV1 and LESV had very different effects on starch composition. In

431

ESV1-over-expressing plants, amylose content was similar to that of esv1 mutants – 14

432

about 60% higher than in wild-type plants. Over-expression of LESV drastically reduced

433

amylose content (Figure 6A; Supplemental Figure 9B). The variation in amylose content

434

caused by changes in amounts of ESV1 or LESV was not attributable to variation in the

435

amount of GBSS in the starch. There was no obvious relationship between GBSS

436

protein levels and amylose levels in plants with different levels of ESV1 or LESV

437

proteins (Supplemental Figure 9C). Loss or over-expression of ESV1 caused only minor

438

alterations in the levels of starch-bound phosphate. Changes in LESV content had

439

greater effects: loss of LESV reduced starch phosphate content by about 25%, whereas

440

over-expression increased it by about 75% (Supplemental Figure 9D).

441 442

We investigated whether altered levels of ESV1 or LESV proteins affect the structure of

443

the major starch polymer amylopectin. The chain length distribution of amylopectin

444

differed slightly between esv1, lesv, ESV1 and LESV over-expressing lines and wild-

445

type starch (Figure 6B), but these differences were small in comparison to those seen in

446

Arabidopsis mutants deficient in enzymes of starch biosynthesis (e.g. Pfister et al.,

447

2014). The chain length distributions of the β-limit dextrin of amylopectin from esv1 and

448

lesv mutants were also similar to those of wild-type plants (Supplemental Figure 9F).

449

These results indicate that neither chain elongation nor branching/debranching during

450

amylopectin biosynthesis are strongly affected by the loss of ESV proteins.

451 452

Loss of ESV1 Modifies but Does Not Abolish the Effects of Mutations That

453

Reduce Starch Degradation

454 455

We considered the possibility that ESV1 interacts with a specific enzyme of starch

456

degradation and inhibits its activity. To test this idea, we created a series of double

457

mutant plants that lacked both ESV1 and a specific enzyme shown previously to be

458

involved in starch degradation in the chloroplast at night (Figure 7). Interpretation of the

459

phenotypes of mutants lacking enzymes of starch degradation is complicated by

460

redundancy within the pathway: for most enzymes, loss does not prevent starch

461

degradation but rather reduces the proportion of end-of-day starch that is consumed

462

during the night (Scialdone et al., 2013). Nonetheless, if ESV1 inhibits one specific 15

463

enzyme of starch degradation, the loss of both this target enzyme and ESV1 should

464

give a starch turnover phenotype identical to that given by the loss of the target enzyme

465

alone. The loss of both ESV and any enzyme other than its direct target might be

466

expected to give a phenotype different from that of either parent mutant.

467 468

We first investigated the effect of losing ESV1 on the phenotype of the sex1 mutant,

469

which lacks GWD . This enzyme catalyzes the first step in degradation: the

470

phosphorylation of the C6-position of glucose residues in amylopectin (Yu et al., 2001;

471

Ritte et al., 2002). Starch degradation in sex1 mutants is severely impaired and starch

472

accumulates to high levels. It is thought that phosphorylation at the granule surface by

473

GWD disrupts the organization of the starch polymers and facilitates degradation by

474

starch-hydrolyzing enzymes (Blennow and Engelsen, 2010; Hejazi et al., 2008; Hejazi

475

et al., 2009; Hejazi et al., 2010). The esv1 sex1 mutants had leaf starch contents that

476

were different from those of either parent. Starch contents were significantly lower than

477

those of the sex1 mutant, but much higher than those of esv1 and wild-type plants.

478

(Figure 7A, Supplemental Figure 9F). However, the starch phenotype of the double

479

mutant was very different in other parts of the plant. Whereas sex1 mutants accumulate

480

very high levels of starch throughout the plant, including root caps (Caspar et al., 1991;

481

Yu et al., 2001), the esv1 sex1 mutant had very low levels of starch in root cap cells,

482

similar to the esv1 parent (Figure 4B). The esv1 sex1 mutant also had little starch in the

483

major veins and petioles of mature leaves and around the meristem. These regions had

484

very high starch contents in the sex1 parent (Figure 7B).

485 486

As for the esv1 sex1 mutant, leaves of other double mutants lacking both ESV1 and a

487

protein necessary for starch degradation had starch turnover phenotypes different from

488

those of their parents (Figure 7, C and D). Loss of ESV1 reduced the severity of the

489

starch-excess phenotype in all mutant backgrounds examined, but double mutants

490

retained more starch in the rosette at the end of the night than esv1 and wild-type

491

plants. The starch contents of leaves of the esv1 sex4 and esv1 pwd double mutants

492

were reduced relative to that of sex4 (lacking a glucan phosphate phosphatase: Kötting

493

et al., 2009) and pwd (lacking phosphoglucan water dikinase: Baunsgaard et al., 2005; 16

494

Kötting et al., 2005) by about one third or less and were always greatly in excess of the

495

starch content of wild-type plants. Mutants lacking ESV1 and either the major

496

chloroplastic β-amylase BAM3, the β-amylase-like protein BAM4 (Fulton et al., 2008) or

497

the glucan phosphate phosphatase-like protein LSF1 (Comparot-Moss et al., 2010) had

498

starch contents that were reduced by 50% or more relative to their starch excess parent

499

at the end of the day and the end of the night. Because loss of ESV1 reduced the

500

severity of the starch excess phenotype in all of the starch degradation mutants, we

501

suggest that its role in starch degradation is not primarily as a negative regulator of any

502

one of the five starch-degradation proteins affected in the mutants. However, these

503

results must be interpreted with caution because it remains possible that – as for sex1 -

504

the impact of loss of ESV1 on the starch-excess phenotypes is different in

505

photosynthetic and non-photosynthetic cells of the rosette.

506 507

We also investigated whether the high rate of starch degradation in the esV1 mutant

508

results from activation of a degradative enzyme that is redundant with respect to starch

509

degradation in wild-type plants. Accordingly, we crossed esV1-1 with mutants lacking

510

either the β-amylase BAM1 or the α-amylase AMY3. Loss of either of these

511

chloroplastic proteins has no effect on starch turnover in a wild-type background (Yu et

512

al., 2005; Fulton et al., 2008). We found that the esV1 bam1 and esV1 amy3 double

513

mutants had the same very low starch contents at the end of the night as the esV1

514

parent (Figure 7, C and E); thus, the high rate of starch degradation in esv1 does not

515

require the activity of either BAM1 or AMY3.

516 517

DISCUSSION

518 519

This work describes two proteins with central and previously unsuspected roles in plant

520

carbohydrate metabolism. We show that the presence and/or levels of the ESV1 and

521

LESV proteins are crucial for normal patterns of starch biosynthesis and degradation in

522

the Arabidopsis plant, and hence for normal growth and productivity. Our results lead us

523

to propose that both proteins are involved in determining the conformation of the starch

524

granule matrix, and that the perturbations of starch turnover in plants with altered levels 17

525

of the proteins are indirect consequences of abnormal starch granule structures. Below

526

we present the evidence for this proposal and discuss the wider significance of our

527

findings for the understanding of carbon storage and allocation in plants.

528 529

ESV1 and LESV Are Important for Normal Starch Turnover in Leaves

530 531

Our starvation reporter screen identified ESV1 as a protein required for normal rates of

532

starch degradation in Arabidopsis leaves at night. Whereas in wild-type plants, starch

533

reserves are used at an essentially linear rate at night such that they are exhausted

534

almost exactly at dawn, in the absence of ESV1, starch reserves are consumed in a

535

non-linear manner and are exhausted prior to dawn. Accelerated starch degradation at

536

night was previously reported in the double mutant cca1 lhy, which lacks two MYB

537

transcription factors that are central components of the circadian clock and thus

538

anticipates dawn several hours too early (Graf et al., 2010). It seems highly unlikely that

539

the rapid starch degradation in esv1 mutants is brought about by a defective circadian

540

clock. The mutant does not exhibit short-period phenotypes such as early flowering, and

541

the ESV1 protein is located in the chloroplast whereas known clock components are

542

nuclear and cytosolic.

543 544

The esv1 phenotype is similar in several respects to that of isa mutants. The isa1 and

545

isa2 mutants also exhaust their glucan reserves prior to dawn and appear to degrade

546

some storage glucan during the day as well as at night (Delatte et al., 2005). The

547

reserve glucans in these mutants are largely soluble rather than in granular form

548

(Zeeman et al., 1998; Delatte et al., 2005). Loss of normal control of degradation is

549

expected in this case. Our chemical kinetic models capable of explaining the

550

dependency of starch degradation on starch content and time until dawn require the

551

interaction of controlling molecules with the solid surface of the starch granule, and they

552

predict that control will be lost if glucans are soluble (Scialdone et al., 2013). esv1

553

mutants accumulate glucans as starch granules, so loss of control of glucan

554

degradation at night is not due to a radical change in glucan solubility. However, the

555

loss of control may be generally analogous to that in isa mutants in that it may stem 18

556

from increased accessibility of glucans to starch-degrading enzymes. We discuss this

557

possibility below.

558 559

LESV is less important for normal starch turnover than ESV1. The pattern of starch

560

turnover in the lesv mutant is similar to that of wild-type plants. The involvement of

561

LEST in processes underpinning normal starch metabolism is apparent from the

562

significantly altered levels of amylose and of granule-bound phosphate in the lesv

563

mutant and from the effects of LESV over-expression in Arabidopsis leaves. As with

564

esv1 mutants, plants with elevated LESV consume starch too rapidly during the night so

565

that reserves are exhausted before dawn. By contrast, over-expression of ESV1 results

566

in very high starch levels throughout the day-night cycle. Overall, our results show that

567

control of ESV1 and LESV protein levels is essential for normal starch turnover in

568

Arabidopsis leaves.

569 570

The ESV1 and LESV Proteins Are Directly Associated with Starch Granules

571 572

Experiments employing fluorescently-tagged proteins confirmed findings from previous

573

proteomic studies that both ESV1 and LESV proteins are located inside the chloroplast.

574

Together with cell fractionations, these experiments revealed that both proteins are at

575

least in part associated with starch granules. Changes during the day-night cycle in the

576

fraction of the proteins that is granule bound suggest that both are incorporated into

577

granules as they grow during the day and are then released as they are degraded at

578

night. Although LESV transcript levels vary several-fold over the day-night cycle, being

579

low at the end of the day and high at the end of the night, the total amount of LESV

580

protein is similar at these two time points. Thus, it appears that large daily changes in

581

the transcript level of LESV are not reflected at the protein level, a feature that LESV

582

shares with several other proteins of starch metabolism (Skeffington et al., 2014).

583

Equivalent information is not available for ESV1 because it is not represented on

584

commonly-used Arabidopsis microarrays. Most of the proteins previously reported to be

585

present in both starch granules and the stroma are enzymes of starch biosynthesis and

586

degradation or proteins that interact directly with these enzymes. They include various 19

587

isoforms of starch synthase and starch-branching enzyme (e.g. Denyer et al., 1993; Mu-

588

Forster et al., 1995; Grimaud et al., 2008) and the enzymes responsible for

589

phosphorylation and dephosphorylation of the starch granule surface (e.g. Yu et al.,

590

2001; Kötting et al., 2009; Santelia et al., 2011). The dual location of ESV1 and LESV

591

proteins in the stroma and the starch granule thus indicates a direct role in the

592

biosynthesis, assembly and/or degradation of starch granules.

593 594

Proteins that bind to starch usually possess either Carbohydrate Binding Modules

595

(CBMs) or Surface binding Sites (SBSs) or both. Arabidopsis proteins with one or more

596

of these domains include starch synthases, starch branching and debranching

597

enzymes, amylases, enzymes that phosphorylate or dephosphorylate the starch granule

598

surface, and the recently-characterized PTST (PROTEIN TARGETING TO STARCH)

599

protein that facilitates the binding to starch of granule-bound starch synthase (e.g.

600

Palopoli et al., 2006; Glaring et al., 2011; Meekins et al., 2013; Seung et al., 2015).

601

Starch-binding CBMs fall into distinct classes, for example, GWD and AMY3 both

602

possess two CBM45 domains in tandem and isoamylases, branching enzymes and

603

PTST possess CBM48 domains. SBSs are generally highly variable and require

604

experimental definition (e.g. Cockburn et al., 2013; Meekins et al., 2014). Neither ESV

605

nor LESV possesses a recognizable CBM, and further work will be required to discover

606

whether they possess SBSs. However, the observation that both proteins contain a

607

large domain strongly enriched in conserved tryptophan residues, and with conserved

608

phenylalanine and tyrosine resides, leads us to speculate that they may be capable of

609

binding to specific glucan structures. Known starch-binding CBMs are 90-130 amino

610

acids long and typically have four conserved aromatic amino acid residues that align the

611

protein with the non-polar faces of glucose residues in starch chains (Christiansen et al.,

612

2009). Although they are often found singly, CBMs also occur in tandem – for example

613

STARCH SYNTHASE3 from Arabidopsis has three CBMs (Palopoli et al., 2006). The

614

tryptophan-rich regions of ESV1 and LESV contain over 35 conserved aromatic amino

615

acid residues arranged in motifs. It seems possible that the repeated motifs may

616

mediate binding to numerous glucans or facilitate interaction with long glucans at

617

numerous contact points. 20

618 619

ESV1 Probably Acts Upstream of Enzymes of Starch Degradation

620 621

To discover whether ESV1 might inhibit the action of a protein required for starch

622

degradation, we examined the effect of loss of ESV1 on the phenotypes of starch-

623

excess mutants lacking individual proteins involved in starch degradation. Loss of ESV1

624

reduced but did not abolish the starch excess phenotypes of the sex1, sex4, lsf1, bam3

625

and bam4 mutants. Loss of ESV1 from plants lacking AMY3 or BAM1 - starch-

626

degrading enzymes that are largely redundant in a wild-type background (Yu et al.,

627

2005; Fulton et al., 2008) - resulted in esv1-like phenotypes. The pathway of starch

628

degradation is not amenable to straightforward genetic analysis because of multiple

629

redundancies. Nonetheless, these results allow us to make three proposals about the

630

role of ESV1 as follows.

631 632

First, ESV1 is unlikely to act through direct inhibition of SEX1, PWD, SEX4, LSF1,

633

BAM3 or BAM4. Second, the accelerated rate of starch degradation in the esv1 mutant

634

is largely catalyzed by the same enzymes and proteins responsible for starch

635

degradation in wild-type plants. SEX1, SEX4, LSF1, BAM3 and BAM4 are all required

636

for the accelerated rate of degradation in the esv1 mutant, and we found no evidence

637

that this accelerated rate requires either of the normally “redundant” enzymes, AMY3

638

and BAM1. Third, in mesophyll cells, ESV1 may act at or prior to the start of the

639

pathway of starch degradation. Loss of ESV1 affected the starch excess phenotypes of

640

mutants lacking both downstream and initial enzymes of the pathway. The observation

641

that it had smaller effects on the starch excess phenotypes of sex1, pwd and sex4 than

642

on those of mutants lacking β-amylases could indicate that ESV1 function is linked more

643

closely to the phosphorylation/dephosphorylation of the granule surface than to the

644

actions of downstream hydrolytic enzymes. However, as discussed above, there is no

645

evidence that ESV1 directly inhibits the enzymes of phosphorylation/dephosphorylation,

646

and its loss and overexpression have only small effects on levels of starch-bound

647

phosphate.

648 21

649

ESV1 May Play a Role in Determining the Organization of Glucans in the Granule

650

Matrix

651 652

Given that ESV1 binds to starch granules and appears to act upstream of the pathway

653

of starch degradation, it seems possible that it influences the accessibility of glucans

654

within the matrix to proteins involved in starch degradation. Any such influence is

655

unlikely to be due to a direct effect of ESV1 on starch polymer composition or structure.

656

Loss or over-expression of either ESV1 or LESV affected the amylose content of starch,

657

and had minor effects on amylopectin chain length distribution. However, these changes

658

are highly unlikely to account for the accelerated rate of starch degradation in the esv1

659

mutant. First, both esv1 and lesv mutants have elevated amylose contents, but only the

660

esv1 mutant has accelerated starch degradation. Second, abolition of amylose from the

661

esv1 mutant (in the esv1 gbss mutant) does not alter its starch degradation phenotype.

662

Third, the alterations in chain length distribution of amylopectin are much smaller than in

663

many mutants that do not exhibit accelerated starch degradation (e.g. starch synthase

664

mutants: Szydlowski et al., 2009; Pfister et al., 2014).

665 666

Our data indicate that ESV1 and LESV directly influence the organization of the starch

667

granule matrix rather than the biosynthesis of the starch polymers that comprise it. The

668

esv1 and lesv mutants have irregular granule shapes, and over-expression of ESV1 and

669

LESV dramatically alters granule size, morphology and - in the case of LESV –

670

abundance. ESV1 over-expressers have large, thick granules, whereas LESV over-

671

expressers have increased numbers of granules of very variable size and shape. These

672

effects are unlikely to be due to altered amylopectin structure in the over-expressing

673

lines: the alterations in amylopectin structure are minor compared with those of other

674

starch mutants that retain normal granule shapes and numbers.

675 676

Little is known about how starch polymers are organized to form the starch granule

677

matrix. Waigh and colleagues (Waigh et al., 1998; Waigh et al., 2000) showed that

678

many physical properties of starch – including gelatinization, freezing/thawing and

679

hydration/dehydration - can be explained by modelling it as a chiral side-chain liquid 22

680

crystalline polymer, a self-organizing structure. By extrapolation, and in the absence of

681

a biological explanation, it has been assumed that starch chains synthesized at the

682

granule surface in vivo assemble into double helices that self-organize through physical

683

processes, forming the semi-crystalline lamellar structure of the matrix. The discovery of

684

the ESV1 and LESV proteins raises the possibility that correct assembly of the matrix

685

involves proteins that bind to starch polymers, in addition to physical processes.

686 687

A direct role for ESV1 in the correct organization of the granule matrix could explain why

688

its loss perturbs the rate of starch degradation. We suggest that the degree to which

689

hydrolytic enzymes can access starch polymers is in part determined by ESV1. By

690

helping to confer a high level of matrix organization, ESV1 could restrict access by

691

hydrolytic enzymes. If, in the absence of ESV1, the matrix is less organized, it could be

692

much more susceptible to degradation by hydrolytic enzymes than in wild-type plants.

693

In wild-type plants, control over the rate of starch degradation may be exerted via

694

modulation of the level of starch phosphorylation. Phosphorylation of glucose residues

695

at the granule surface by GWD and PWD is thought to interfere with the regular packing

696

of double helices, thus rendering the starch polymers more accessible to hydrolytic

697

enzymes (Blennow and Engelsen, 2010; Hejazi et al., 2008; Hejazi et al., 2009; Hejazi

698

et al., 2010). We recently presented evidence that phosphorylation may be the point at

699

which circadian clock control over starch degradation is exerted. In the absence of

700

PWD, the rate of starch degradation could no longer be correctly adjusted to the length

701

of the night (Scialdone et al., 2013). It seems plausible that an altered level of starch

702

matrix organization in esv1 mutants may render the modulation of phosphorylating

703

activities less effective as a means of controlling the rate of degradation.

704 705

Although there is no clear starch turnover phenotype in the lesv mutant, the altered

706

amylose and starch phosphate contents of the mutant starch and the large effects of

707

LESV over-expression on starch degradation and granule size, shape and number

708

indicate that LESV too may play a role in matrix organization and in modulating the

709

effects of phosphorylation on the accessibility of starch polymers to hydrolytic enzymes.

710

However ESV1 and LESV appear to have opposite roles in these respects, ESV1 23

711

promoting a high level of organization of the granule matrix and LESV potentially

712

reducing the level of organization. These differences in function may stem from

713

divergent features of the protein sequences. The two proteins share a tryptophan-rich

714

domain, in which many of the aromatic amino acid residues are contained within distinct

715

motifs. In ESV1 sequences, there are 39 conserved aromatic residues within this

716

domain. LESV sequences contain 37 conserved aromatic residues in the same region,

717

30 of which are in the same amino acid position as in ESV1 (Supplemental Figure 4A).

718

However, motif sequences differ between the two proteins. Comparisons of ESV1 and

719

LESV proteins from a wide range of starch-synthesizing organisms shows that many

720

motifs in the domain are conserved within either ESV1 or LESV but not shared between

721

them. ESV1 also possesses a proline-rich terminal domain that is absent from LESV,

722

and LESV has a unique >100-amino-acid N-terminal domain. Further information about

723

the roles of these proteins will come from investigation of the starch- and glucan-binding

724

properties of wild-type and mutant forms combined with detailed physical analyses of

725

starches from plants expressing wild-type and mutant forms at different levels.

726 727

The Importance of ESV1 and LESV for Starch Degradation Varies between Organs

728 729

Examination of the esv1 mutant indicates that the importance of the ESV1 protein for

730

normal starch turnover varies from one organ and tissue to another. In leaves, loss of

731

ESV1 accelerates starch degradation and also reduces the extent of starch

732

accumulation during the day. The reduction in starch accumulation may be due to the

733

occurrence of starch degradation in the light as well as in the dark: the accumulation of

734

maltose – the main product of starch degradation in the chloroplast (Weise et al., 2004)

735

- in esv1 leaves during the day supports this idea. Thus, loss of ESV1 may render leaf

736

starch granules partially accessible to hydrolytic enzymes during the day, as well as

737

increasing access for these enzymes during the night.

738 739

In other organs of the plant, loss of ESV1 has a more radical effect on starch turnover.

740

The starch content of root columella cells and of stems in esv1 plants is very much

741

lower in mutant than in wild-type plants. In sex1 mutant plants, loss of ESV1 has little 24

742

effect on the starch excess phenotype in mesophyll cells of the leaf but abolishes starch

743

accumulation in root cap, vein and petiole cells. These observations indicate that the

744

importance of ESV1 for normal starch metabolism is dependent on the dynamics of

745

starch turnover. In wild-type leaves, little or no starch degradation occurs during most of

746

the light period: the rate of starch biosynthesis and starch accumulation are the same.

747

In the absence of ESV1, some degradation occurs during the day, but this rate is far

748

lower than the rate of biosynthesis. In other organs of the plant, starch may be subject

749

to simultaneous biosynthesis and degradation at all times. In the embryo, for example,

750

the effects of loss of starch degrading enzymes on starch content in early development

751

reveal that starch biosynthesis and degradation occur simultaneously throughout the

752

period of accumulation and loss of starch (Andriotis et al., 2010). Thus, we anticipate

753

that loss of ESV1 effectively prevents starch storage in organs other than the leaf

754

because it permits faster degradation without a change in the rate of biosynthesis.

755 756

Are ESV1 and LESV part of the control mechanism that adjusts the rate of leaf starch

757

degradation according to time of dawn and starch content? Such a role cannot be ruled

758

out, but at present there is no evidence that the actions of these proteins are subject to

759

either transcriptional or post-translational control on a short-term basis. The wide

760

occurrence and conservation of both proteins, and the importance of ESV1 in many

761

organs of the plant, indicate that they are fundamental components of the starch

762

biosynthesis and turnover machinery in general.

763 764

METHODS

765 766

Plant Materials and Growth Conditions

767 768

All Arabidopsis thaliana mutants were in the Col-0 background. Starch-excess mutants

769

were described previously (amy3, Yu et al., 2005; bam mutants, Fulton et al., 2008;

770

gbss, Seung et al., 2015; isa1, Delatte et al., 2005; lsf1, Comparot-Moss et al., 2010;

771

sex1, Yu et al., 2001; sex4, Kötting et al., 2009; pwd, Kötting et al., 2005). The lines 25

772

were as follows: amy3-2 (SAIL_613 D12); bam1-1 (SALK_039895); bam3-1 (CS92461);

773

bam4-1 (SALK_037355); gbss (GABI_914G01), isa1 (SALK_042704); lsf1-1

774

(SALK_100036); sex1-3 (Yu et al., 2001); sex4-3 (SALK_102567); pwd

775

(SALK_110814). Double mutants were selected from the F2 of crosses between esv1

776

and a second mutant line. Retention of both T-DNA insertions was checked for all

777

double mutant lines. Primers are listed in Supplemental Table 5.

778 779

Unless otherwise stated, Arabidopsis thaliana plants were grown on soil in growth

780

cabinets or controlled environment rooms with a 12 h light/12 h dark cycle, at 20°C and

781

at 150-180 µmol photons m-2 s-1 (metal halide lamps).

782 783

Generation and Mutagenesis of the Starvation Reporter Line

784 785

An A. thaliana starvation reporter line was established by transforming wild-type plants

786

with a construct consisting of a fusion of the starvation responsive promoter from

787

At1g10070 and the luciferase (LUC) gene (Graf et al., 2010). To create M0 seeds, 400

788

mg of homozygous T4 seeds of the reporter line were incubated in a 50-mL tube in 10

789

mL ethyl-methanesulphonate (EMS) solution [0.15% (v/v) EMS, 0.02% (v/v) Tween20]

790

on a rotating shaker at room temperature for 20 h. Seeds were washed with 10 mL

791

0.02% (v/v) Tween20 solution for 10 min on a rotating shaker. This washing step was

792

repeated 12 more times, and the seeds were then combined with 330 mL 0.1% (w/v)

793

aq. agar and pipetted on soil in 121 cm2 pots (10 mL per pot). The pots were held at 4°C

794

for 3 d and transferred to a controlled environment room. At 12 d old, five seedlings

795

were transferred to each of 1500 pots (121 cm2) in a greenhouse and allowed to set M2

796

seed. Seed from each pot was collected as a single pool.

797 798

Imaging of Bioluminescence

799 800

Seedlings or plants grown on soil were sprayed with luciferin solution [0.8 mM luciferin,

801

0.01% (v/v) Triton X-100] 24 h before imaging, then sprayed again 60 min before

802

imaging. Plants were transferred to the NightOwl CCD camera system (Berthold 26

803

Technologies; www.berthold.com) and bioluminescence was assayed using either

804

Indigo software (imaging settings: 1 min exposure time, 2x2 binning) or WinLight

805

software (imaging settings: 1 min exposure time, medium resolution, pixel binning 4x4,

806

single frame accumulation according to the manufacturer’s instructions).

807 808

Mapping of the ESV1 Gene

809 810

The esv1 mutant was out-crossed to the A. thaliana accession Landsberg erecta (Ler).

811

About 20 F1 plants were allowed to self-pollinate and set seed. F2 plants were grown for

812

10 d before identification of individuals showing bioluminescence by the end of the

813

normal night.

814 815

DNA was extracted from individual homozygous F2 plants. To obtain a rough map

816

position for the mutation, each plant was genotyped using genetic markers that were

817

distributed over the five chromosomes (Supplemental Table 1). Recombination

818

frequencies for each marker were calculated as the percentage of Ler polymorphisms

819

detected at that locus. To obtain a smaller interval, new polymorphic markers upstream

820

and downstream of the marker indicating the lowest recombination frequency were

821

analyzed. New markers were obtained using the A. thaliana mapping platform (Hou et

822

al., 2010; http://amp.genomics.org.cn/).

823 824

Identification of the ESV1 Gene

825 826

The mutation underlying the esv1 phenotype was identified by genome resequencing.

827

About 500 homozygous mutant plants were selected from the F2 of the cross between

828

the mutant and Ler. DNA was prepared from nuclei extracted from leaves. Library

829

construction from 5 µg of RNA-free gDNA, cluster generation, and sequencing on one

830

lane on the Illumina GAIIx platform were carried out by The Genome Analysis Centre

831

(Norwich, UK). The raw data were analysed as follows. Maq v0.71 (Li et al., 2008) was

832

used to align the 35.6M 100 base, paired-end reads against the TAIR8 Col-0 reference

833

sequence, producing about 28X coverage, and thus to generate a list of raw SNPs. The 27

834

maq.pl Perl script was employed to filter the SNPs on quality criteria, and the survivors

835

were used as input to a post-processing script that first eliminated SNPs that

836

corresponded to known Col-0/Ler polymorphisms (http://signal.salk.edu/atg1001/data/)

837

and then retained only EMS candidates from the remainder. The output of this script

838

was a GFF file that was loaded into a local instance of the GBrowse genome browser

839

(Stein et al., 2002), together with the TAIR8 pseudochromosome sequences and gene

840

model annotations, allowing visual inspection through a web browser. By interrogating

841

the GBrowse MySQL database with a Perl script using Bio::DB::GFF methods, a

842

genome-wide list of EMS candidates (G/C -> A/T) within annotated gene sequences

843

and inferred to induce either non-synonymous codon or donor/acceptor splice site

844

mutations was produced, which was then further refined based on chromosomal

845

location.

846 847

Quantification of Transcripts for Starvation Marker Genes

848

Quantitative reverse-transcription PCR (RT-qPCR) was used to measure transcript levels of

849

starvation marker genes, using the primers and methods described in Graf et al. (2010). Briefly,

850

total RNA was extracted from entire Arabidopsis rosettes harvested at the end of night, using an

851

RNeasy Plant RNA purification kit (www.qiagen.com). Following DNase treatment, 2 μg total

852

RNA was used for reverse transcription using RevertAid Reverse Transcriptase

853

(www.thermofisher.com), and qPCR analysis was performed using the Fast SYBR Green

854

master mix together with a 7500 Fast Real-Time PCR system (Applied Biosystems,

855

www.thermofisher.com). Transcript levels were calculated relative to the YLS8 housekeeping

856

gene. The primer pairs are listed in Supplemental Table 5.

857 858

Measurement of Starch and Sugar Contents

859 860

Starch and soluble sugars were extracted and quantified as previously described

861

(Critchley et al., 2001; Delatte et al., 2005; Martinis et al., 2014). Briefly, rosettes were

862

harvested into liquid nitrogen and ground to a fine powder in a ball mill. The powder was

863

suspended and agitated in ice-cold 0.7 M perchloric acid. Following centrifugation, the

864

pellet was washed three times in 80% (v/v) ethanol, resuspended in water, heated to

865

gelatinize the starch, digested with α-amylase and amyloglucosidase, and assayed 28

866

enzymatically for glucose. For sugars, the supernatant was neutralized, passed through

867

sequential cation- and anion-exchange columns (Dowex 50 and Dowex 1), and

868

analysed by High Performance Anion Exchange Chromatography with Pulsed

869

Amperometric Detection (HPAEC-PAD; www.dionex.com).

870 871

Visualization of Starch in Tissues

872 873

For iodine staining, tissues were decolorized in hot 80% (v/v) ethanol, rinsed in water,

874

stained in Lugol’s iodine solution and rinsed again. For visualization of starch in the

875

stem, tissue sections were stained using the modified Pseudo-Schiff Propidium Iodide

876

(mPS-PI) staining method described by Truernit et al. (2008). Stained sections were

877

imaged using confocal laser scanning microscopy, as described below.

878 879

Analysis of Starch Structure and Composition

880 881

The chain length distribution of amylopectin was profiled as described by Streb et al.

882

(2008). Briefly, starch in the pellet from the perchloric acid extraction (see above) was

883

debranched with isoamylase from Pseudomonas sp. (www.sigmaaldrich.com) and

884

pullulanase M1 from Klebsiella planticola (www.megazyme.com). The resulting glucan

885

chains were purified by passage through sequential cation- and anion-exchange

886

columns and analyzed by HPAEC-PAD.

887 888

Analysis of Granule Morphology, Amylose Content and Starch-Bound Phosphate

889 890

For the determination of amylose content, starch-bound phosphate, and granule

891

morphology by Scanning Electron Microscopy (SEM), starch granules were purified

892

from four-week–old Arabidopsis rosettes as described by Seung et al. (2015). The

893

apparent amylose content of the starch was determined using the iodine colorimetry-

894

based method described by Zeeman et al. (2002). Starch-bound phosphate was

895

quantified as described in Santelia et al. (2011). Hydrolyzed starch was

896

dephosphorylated with Antarctic phosphatase (New England Biolabs; www.neb.com), 29

897

and the phosphate released was quantified using malachite green. Granule morphology

898

was examined using a Merlin Field Emission Scanning Electron Microscope

899

(www.zeiss.com).

900 901

To visualize granule morphology within chloroplasts, segments from young leaves of

902

three-week-old plants were fixed in glutaraldehyde followed by osmium tetroxide and

903

embedded in Epon resin as described in Seung et al. (2015). Light microscopy images

904

of toluidine blue-stained sections were acquired on an AxioImager Z2 microscope fitted

905

with a 100× oil-immersion lens with 1.4 numerical aperture and an AxioCam

906

monochrome camera (www.zeiss.com).

907 908

For transmission electron microscopy, ultrathin (70 µm) sections were cut with a

909

diamond knife and placed on formvar carbon-coated copper grids, stained with 2% (w/v)

910

uranyl acetate and Reynold’s lead citrate and imaged with a FEI Morgagni 268 electron

911

microscope (www.fei.com). Pictures are representative of sections from two individual

912

plants per genotype.

913 914

Expression Vectors for YFP-Fusion Proteins in Planta and Plant Transformation

915 916

The coding sequences for Arabidopsis ESV1 and LESV were amplified from the full

917

length cDNA clones RAFL09-78-O20 and RAFL16-10-H06, respectively (RIKEN

918

Bioresource Centre; epd.brc.riken.jp), using primers flanked with attB recombination

919

sites. Primers are listed in Supplemental Table 5. The amplified inserts were then

920

recombined into the Gateway-compatible entry vector pDONR221 (Invitrogen;

921

www.thermofisher.com) and then recombined into the expression vector, pB7YWG2

922

(Karimi et al., 2002), downstream of the CaMV 35S promoter and in frame with the C-

923

terminal YFP tag. Expression constructs were transformed into Agrobacterium

924

tumefaciens strain GV3101.

925 926

Transient expression of YFP-tagged proteins in Nicotiana sylvestris (wild-type and pgm)

927

leaves was achieved by infiltrating A. tumefaciens cells into the abaxial epidermis 30

928

(Seung et al., 2015). Stable transformation of Arabidopsis was also carried out using the

929

floral-dipping method, as described by Zhang et al. (2006). Transformants were

930

identified in the T1 generation based on their resistance to the herbicide Basta.

931

Homozygous plants were identified in the T2 generation based on the segregation rates

932

of the Basta-resistance gene.

933 934

Detection of ESV1 and LESV Proteins by Silver Staining and Immunoblotting

935 936

For extraction of total (soluble and insoluble) proteins, two young leaves from individual

937

four-week-old rosettes were harvested and homogenized using a pestle in

938

microcentrifuge tubes in 300 µL extraction medium [40 mM Tris-HCl, pH 6.8, 5 mM

939

MgCl2, 2% (w/v) SDS, Complete Protease Inhibitor (www.roche.com)]. Insoluble debris

940

was pelleted at 20,000g. The protein concentration of the supernatant was determined

941

using the Pierce BCA Protein Assay kit (www.thermofisher.com), and the indicated

942

amounts of protein were loaded onto SDS-PAGE gels.

943 944

For the fractionation of soluble and insoluble proteins from N. sylvestris leaves, 7 mm

945

leaf discs were collected from transformed leaves at the end of the photoperiod three

946

days after infiltration, and snap frozen in liquid N2. Discs were homogenized using a

947

pestle in microcentrifuge tubes in 100 µL extraction medium [40 mM Tris-HCl, pH 6.8, 5

948

mM MgCl2, Complete Protease Inhibitor (www.roche.com)]. Insoluble debris was

949

pelleted at 20,000g. The pellet was washed once in extraction medium, then

950

resuspended in 100 µL SDS-PAGE loading medium [50 mM Tris-HCl, pH 6.8, 3% (w/v)

951

glycerol, 2% (w/v) SDS, 100 mM DTT and 0.005% (w/v) bromophenol blue) .The

952

suspension was heated at 95°C for 5 min, and insoluble debris were removed by

953

centrifugation. The supernatant was diluted with 10X SDS-PAGE loading medium.

954 955

Granule-bound proteins were extracted from purified starch granules (prepared as

956

described above) using the method described by Seung et al. (2015).

957

31

958

Silver staining was performed with the Silver Stain Plus kit (www.biorad.com). For

959

immunoblotting, proteins were transferred onto a PVDF membrane following SDS-

960

PAGE and probed with antisera specific to ESV1 or LESV. Antisera were raised in

961

rabbits against recombinant ESV1 or LESV proteins expressed in and purified from E.

962

coli. YFP-tagged ESV proteins were detected with an anti-GFP antiserum

963

(www.clontech.com). Plant actin was detected with a commercial monoclonal antibody

964

(Sigma A0480). Dilutions of antisera were as follows: Anti-ESV1, 1:1000; anti-LESV,

965

1:3000; anti-GFP, 1:5000; anti-actin,1:10000.

966 967

Confocal Laser Scanning Microscopy

968 969

Confocal laser scanning microscopy was carried out on an LSM 780 confocal

970

microscope (Carl Zeiss), with a 40X water-immersion lens (1.1 numerical aperture). For

971

the acquisition of YFP signal, the excitation beam was produced with an argon laser set

972

at 514 nm, and emitted light was captured between 518 to 557 nm. The

973

autofluorescence of chlorophyll was captured between 662 to 721 nm. Images were

974

processed with ImageJ software (http://rsbweb.nih.gov/ij/).

975 976

Phylogenetic Analysis

977 978

To build the phylogenetic tree, ESV1 and LESV sequences were retrieved from the

979

NCBI and 1000 plants (1KP; Johnson et al., 2012; http://www.onekp.com) databases

980

using BLASTp. The alignment was constructed using the MAFFT server (Katoh and

981

Standley, 2013) with the "Auto" alignment strategy. The tree was built using MEGA

982

software version 6 (Tamura et al., 2013), using an LG model, four gamma categories for

983

rate variation, a SPR level 5 method for heuristic search, and a neighbor-joining tree as

984

starting tree. 1,000 bootstrap replicates were used to assess branch support (branch

985

support