Maize VIVIPAROUS1 Interacts with ABA ... - ACS Publications

Mar 27, 2019 - Yumin Zhang , Qichao Sun , Chunxia Zhang , Guanglong Hao , Chunmei Wang , Lynnette M. A. Dirk , A. Bruce Downie , and Tianyong Zhao...
1 downloads 0 Views 1MB Size
Subscriber access provided by UNIV AUTONOMA DE COAHUILA UADEC

Agricultural and Environmental Chemistry

Maize VIVIPAROUS1 Interacts with ABA INSENSITIVE5 to Regulate GALACTINOL SYNTHASE2 Expression Controlling Seed Raffinose Accumulation Yumin Zhang, Qichao Sun, Chunxia Zhang, Guanglong Hao, Chunmei Wang, Lynnette M. A. Dirk, A. Bruce Downie, and Tianyong Zhao J. Agric. Food Chem., Just Accepted Manuscript • DOI: 10.1021/acs.jafc.9b00322 • Publication Date (Web): 27 Mar 2019 Downloaded from http://pubs.acs.org on March 29, 2019

Just Accepted “Just Accepted” manuscripts have been peer-reviewed and accepted for publication. They are posted online prior to technical editing, formatting for publication and author proofing. The American Chemical Society provides “Just Accepted” as a service to the research community to expedite the dissemination of scientific material as soon as possible after acceptance. “Just Accepted” manuscripts appear in full in PDF format accompanied by an HTML abstract. “Just Accepted” manuscripts have been fully peer reviewed, but should not be considered the official version of record. They are citable by the Digital Object Identifier (DOI®). “Just Accepted” is an optional service offered to authors. Therefore, the “Just Accepted” Web site may not include all articles that will be published in the journal. After a manuscript is technically edited and formatted, it will be removed from the “Just Accepted” Web site and published as an ASAP article. Note that technical editing may introduce minor changes to the manuscript text and/or graphics which could affect content, and all legal disclaimers and ethical guidelines that apply to the journal pertain. ACS cannot be held responsible for errors or consequences arising from the use of information contained in these “Just Accepted” manuscripts.

is published by the American Chemical Society. 1155 Sixteenth Street N.W., Washington, DC 20036 Published by American Chemical Society. Copyright © American Chemical Society. However, no copyright claim is made to original U.S. Government works, or works produced by employees of any Commonwealth realm Crown government in the course of their duties.

Page 1 of 34

Journal of Agricultural and Food Chemistry

1

Maize VIVIPAROUS1 Interacts with ABA INSENSITIVE5 to

2

Regulate GALACTINOL SYNTHASE2 Expression Controlling Seed

3

Raffinose Accumulation

4

Yumin Zhanga,b • Qichao Suna,b • Chunxia Zhang a,b •Guanglong Haoa,b • Chunmei

5

Wangc • Lynnette M.A. Dirk d •A. Bruce Downied • Tianyong Zhaoa,b*

6

Running Title: VP1 and ABI5 Regulate Raffinose Biosynthesis

7

aState

8

Northwest A&F University, Yangling, Shaanxi, 712100, China.

9

bThe

Key Laboratory of Crop Stress Biology for Arid Areas, College of Life Sciences,

Key Laboratory of Biology and Genetics Improvement of Maize in Arid Area of

10

Northwest Region, Ministry of Agriculture, Northwest A&F University, Yangling,

11

Shaanxi, 712100, China.

12

c

13

Northwest A&F University, Yangling, Shaanxi, China.

14

dDepartment

15

Environment, University of Kentucky, Lexington, KY 40546, USA.

16

*Corresponding author

17

Tianyong Zhao, Ph.D

18

Professor

19

Department of Biochemistry and Molecular Biology

20

College of Life Sciences, Northwest A&F University

21

Yangling, Shaanxi 712100, China

22

Email: [email protected]

23

Tel: 86-13324517154

The Biology Teaching and Research Core Facility, College of Life Sciences,

of Horticulture, Seed Biology, College of Agriculture, Food, and

1

ACS Paragon Plus Environment

Journal of Agricultural and Food Chemistry

24

Abstract

25

Raffinose, an oligosaccharide found in many seeds, plays an important role in seed

26

vigor; however, the regulatory mechanism governing raffinose biosynthesis remains

27

unclear. We report here that maize W22 wild type (WT) seeds, but not W22 viviparous1

28

(zmvp1) mutant seeds, start accumulating galactinol and raffinose at 28 days after

29

pollination (DAP). Transcriptome analysis of the zmvp1 embryo showed that the

30

expression of GALACTINOL SYNTHASE2 (GOLS2) was down-regulated relative to

31

WT. Further experiments showed that the expression of ZmGOLS2 was up-regulated

32

by ZmABI5, but not by ZmVP1, and it was further increased by co-expression of

33

ZmABI5 and ZmVP1 in maize protoplasts. ZmABI5 interacted with ZmVP1, while

34

ZmABI5, but not ZmVP1, directly binds to the ZmGOLS2 promoter. Together, all the

35

finding suggest that ZmVP1 interacts with ZmABI5 and regulates ZmGOLS2

36

expression and raffinose accumulation in maize seeds.

37

Key words: maize, seeds, viviparous, ABI5, GALACTINOL SYNTHASE, raffinose

38

biosynthesis.

2

ACS Paragon Plus Environment

Page 2 of 34

Page 3 of 34

Journal of Agricultural and Food Chemistry

39

Introduction

40

Raffinose-family oligosaccharides (RFOs) accumulate during seed maturation1 and

41

play important roles in seed vigor2. There are two commitment steps for raffinose

42

biosynthesis. The first step is to use GALACTINOL SYNTHASE (GolS: EC 2.4.1.123)

43

to catalyze galactose transfer from UDP-galactose to myo-inositol to form galactinol3.

44

The second step is to synthesize raffinose by transfer of the galactosyl unit from

45

galactinol to sucrose, which is catalyzed by RAFFINOSE SYNTHASE (RS, EC

46

2.4.1.82)4. We have previously characterized two maize GALACTINOL SYNTHASE

47

genes5 and the unique maize RAFFINOSE SYNTHASE2. Both ZmGOLS2 and ZmRAFS

48

were expressed in the later stages of seed development2, 5. The seeds of Arabidopsis

49

thaliana atgols1 mutant or AtGOLS1 RNAi lines, with decreased RFOs, completed

50

germination faster than those of wild type6, suggesting that RFOs negatively influence

51

the rapidity of completion of seed germination.

52

Maize VIVIPAROUS1 (ZmVP1) encodes a transcription factor7 orthologous to the

53

Arabidopsis thaliana ABA INSENSITIVE3 (AtABI3)8. ZmVP1 is a key regulator of

54

seed maturation, its expression is induced by abiotic stress and ABA in mature embryos

55

or leaves9. Seeds with null alleles at VP1 result in a loss of sensitivity to ABA, leading

56

to vivipary (precocious germination), and loss of desiccation tolerance10. There are two

57

independent zmvp1 mutant lines, vp1-R and vp1-Mc11. VP1 is not expressed in the vp1-R

58

mutant kernels, and the zmvp1-R seeds accumulate only trace amounts of raffinose in

59

the embryos11. The raffinose content in zmvp1-MC seeds is similar to that of the wild

60

type, and the zmvp1-MC seeds acquire desiccation tolerance, suggesting that raffinose 3

ACS Paragon Plus Environment

Journal of Agricultural and Food Chemistry

61

biosynthesis may occur in the absence of a fully functional vp1 gene product11.

62

Similar to the maize zmvp1 mutant, the RFOs (raffinose and stachyose) accumulate

63

in the Arabidopsis thaliana desiccation-tolerant wild-type and atabi3-1 mutant seeds

64

while the RFOs are not detectable in the more penetrant atabi3-5 mutant lines12. During

65

seed maturation, VP1 activated the anthocyanin pathway and inhibited the germination-

66

specific alpha-amylase genes in aleurone cells in the developing seed13. In Arabidopsis

67

thaliana, ABI3 directly controls the seed chlorophyll breakdown process through

68

regulation of STAY GREEN2 (SGR2) expression during seed maturation14. ABA

69

INSENSITIVE5 (ABI5) is a basic leucine zipper (bZIP) transcription factor which

70

plays a key role in the regulation of seed maturation and seed germination15. Medicago

71

truncatula Mtabi5 seed accumulate fewer RFOs than that of the wild type control

72

during seed maturation15. ABI5 acts downstream of ABI3 and activates late

73

embryogenesis programs16. ABI5 binds to the ABRE motif or G-box in the promoter

74

of some ABA responsive genes, such as AtEm1 and AtEm6, and regulates gene

75

expression17. Despite the extensive studies on the regulatory network of ABI3 and ABI5

76

controlling events during late embryogenesis whether the control of RFO production is

77

directly influenced by these transcription factors is unknown. However, alterations in

78

the abundance of raffinose in maize zmvp1 allelic series; in the Arabidopsis thaliana

79

atabi3 allelic series; in atabi5- or Mtabi5-mutant seed leaves little doubt that both the

80

biochemistry of production and/or genetic mechanism of control of RFO production is

81

linked, directly or indirectly, to this transcription factor cascade in multiple species.

82

In this study, we have demonstrated that the two important transcription factors 4

ACS Paragon Plus Environment

Page 4 of 34

Page 5 of 34

Journal of Agricultural and Food Chemistry

83

ZmVP1 and ZmABI5 in the ABA signaling transduction pathway, interact with each

84

other and synergistically regulate ZmGOLS2 expression and raffinose accumulation in

85

maize seed. These findings provide new insights into the molecular mechanism of

86

raffinose biosynthesis.

87

Materials and methods

88

Plant material

89

Maize (Zea mays L.) vp1-R mutant seeds (#326B) in the W22 background, were

90

obtained from the Maize Genetics Cooperation Stock Center and were maintained in

91

Yangling (34°17'47.9"N 108°04'28.2"E). The zmvp1-R homozygous seeds were

92

obtained by self-pollination of the heterozygous stocks. The WT or heterozygous

93

kernels are purple, and the homozygous zmvp1-R mutant seeds are yellow. The WT,

94

zmvp1-R embryos and endosperms were separated from 20, 22, 24, 26, 28, 30, 32, 34

95

and 36 days after pollination (DAP) kernels, then frozen in liquid nitrogen and stored

96

at -75 °Cuntil use.

97

Soluble Carbohydrate Extraction and HPLC–ELSD Analysis of Sugars

98

Soluble sugar extraction followed a published protocol2. Isolated plant tissues (0.5

99

g for embryos, 1.5 g for endosperms) were ground into powder in liquid nitrogen. Four

100

aliquots of 1 mL of 80 % (v/v) ethanol were sequentially added, homogenized to a

101

slurry, and transferred to a 50 mL Falcon tube. Another 3 mL of 80 % ethanol was used

102

to wash the mortar and added to the Falcon tube. The Waters Xbridge amide column

103

was washed using Acetonitrile: H2O (65:35) as a mobile phase at 1.0 mL/min for 5

ACS Paragon Plus Environment

Journal of Agricultural and Food Chemistry

Page 6 of 34

104

separation of soluble carbohydrate components. An evaporative light-scattering

105

detector (ELSD, Waters 2424) was applied to monitor the sugar signal.

106

RNA Seq

107

Maize embryos separated from 26 DAP seeds of W22 or zmvp1-R were sent to

108

BGI-Shenzhen,

China

for

total

RNA

extraction

and

RNA-seq

analysis

109

(https://www.bgi.com/us/sequencing-services/rna-sequencing-solutions/

110

transcriptome-sequencing/ ). RNA library preparation and sequencing were performed

111

with BGISEQ-500 Platform. The data set was analyzed according to the BGI

112

bioinformatics protocols for RNA-seq.

113

Identification of ZmABI5

114

The seed-specific expressed bZIP transcription factors were selected according to

115

the data in the published database18-19. The ZmGOLS2 expression and that of the

116

selected bZIP genes was examined by Pearson correlation analysis using a two-tailed

117

method (IBM SPSS Statistics 20). The expression of 9 bZIP genes were positively

118

correlated with ZmGOLS2 expression. The protein sequences of the 9 selected bZIPs

119

and another 9 known bZIP transcription factors were aligned using CLUSTAL W20 and

120

then input into MEGA721. The neighbor-joining (NJ) method was used to construct a

121

phylogenetic tree22 with the following sequences as given by their accession numbers:

122

GRMZM2G175280 (Zea mays, NP_001147562.1); GRMZM2G060216 (Zea mays,

123

NP_001104893.1);

124

GRMZM2G361611 (Zea mays, NP_001152677.1); GRMZM2G117851 (Zea mays,

GRMZM2G030877

(Zea

6

ACS Paragon Plus Environment

mays,

XP_008677589.1);

Page 7 of 34

Journal of Agricultural and Food Chemistry

125

NP_001105272.2);

GRMZM2G055413

(Zea

mays,

NP_001148839.1);

126

GRMZM2G045236 (Zea mays, NP_001148077.1); GRMZM2G094352 (Zea mays,

127

NP_001131340.1);

128

(Oryza sativa, BAA 83740.1); MtABI5 (Medicago truncatula, XP_003625810.1);

129

AtABF2 (Arabidopsis thaliana, NP_001185157.1); AtABF1 (Arabidopsis thaliana,

130

NP_001189934.1); AtbZIP67 (Arabidopsis thaliana, NP_566870.1); AtbZIP12

131

(Arabidopsis thaliana, NP_565948.1); HvABI5 (Hordeum vulgare, AAO06115.1); and

132

OsABI5 (Oryza sativa, XP_015628684.1).

133

Characterization of zmvp1-R mutant plants

GRMZM2G168079 (Zea mays, AQK 99305.1); OsTRAB1

134

Genomic DNA was isolated from maize embryos using the CTAB method23. PCR

135

was performed to characterize the genotype of the plants using ZmVP1 gene-specific

136

primers (F1 and R1, F2 and R2; as listed in Supplemental Table 1). The primer F2a and

137

R2 were used for RT-PCR analysis of ZmVP1 gene expression (Primers were listed in

138

Supplemental Table 1).

139

RNA extraction and real time RT‑PCR

140

Total RNA was extracted from embryos following a published protocol24. TRIzol

141

Reagent (Takara, Japan) was used for extraction of RNA from protoplasts. The reverse

142

transcription (RT) products were diluted fifty-fold with DEPC-treated water. Different

143

primers (as listed in Supplemental Table 1) were used for real time RT-PCR analysis

144

of gene expression, for ZmGOLS2, primers F3 and R3 were used; for ZmABI5, primers

145

F4 and R4 were used; for ZmGAPDH, primers F5 and R5 were used. The expression of 7

ACS Paragon Plus Environment

Journal of Agricultural and Food Chemistry

146

tested genes was normalized to the ZmGAPDH gene expression. The experiments were

147

repeated at least three times with independent biological samples.

148

Analysis of ZmGOLS2 promoter

149

To further investigate whether ZmVP1 was able to directly bind to the promoter of

150

the ZmGOLS2 gene and up-regulate its transcription, the 1496 bp ZmGOLS2 promoter

151

sequence was analyzed for all putative cis-acting elements using the databases,

152

PlantCARE

153

CARE.html), PLACE (http://www.dna.affrc.go.jp/PLACE/)25 and Neural Network

154

Promoter Prediction (http://www.fruitfly.org/seq_tools/promoter.html).

155

Vector construction

156 157

(http://bioinformatics.psb.ugent.be/webtools/plantcare/html/search_

For all manipulations involving PCR, the primer sequences are listed in Supplemental Table 1.

158

For construction of vectors used for maize protoplasts transformation, the coding

159

sequence of ZmVP1 or ZmABI5 was amplified by PCR from cDNA synthesized from

160

RNA that was isolated from mazie B73 embryos (DAP 20). The PCR products were

161

purified and inserted into pGL3-basic vector between MluI and XbaI as described26.

162

The restriction enzyme sequences were designed in the primers to facilitate cloning (for

163

ZmVP1, primers F6 and R6; for ZmABI5, primers F7 and R7).

164

For construction of vectors used for characterization of the ZmGolS2 promoter

165

activity, overlapping PCR was performed to generate the single mutation or double

166

mutation of the G-box in the ZmGOLS2 promoter in a dual luciferase expression vector, 8

ACS Paragon Plus Environment

Page 8 of 34

Page 9 of 34

Journal of Agricultural and Food Chemistry

167

in which the Rluc expression was controlled by a 1093 bp ZmGOLS2 promoter (named

168

V6)27. For mutation of the ABRE motif at -647 bp, the 5’ part of ZmGOLS2 promoter

169

was amplified by PCR from V6 using primers F8 and R10, the 3’ part of ZmGOLS2

170

promoter was amplified by PCR from V6 using primers F10 and R9. The mutated

171

ZmGOLS2 promoter (M1) was obtained by PCR using a mixture of the 5’ and 3’

172

amplicons as templates with primers F8 and R9. For mutation of the G-box element at

173

-237 bp, the 5’ part of the ZmGOLS2 promoter was amplified by PCR from V6 using

174

primers F8 and R8, and the 3’ part of the ZmGOLS2 promoter was amplified by PCR

175

from V6 using primers F9 and R9. The mutated ZmGOLS2 promoter (M2) was obtained

176

by PCR using the two PCR amplicons as templates with primers F8 and R9. For

177

mutation of both the ABRE motif and the G-box elements at -647 bp and -237 bp, the

178

5’ part of the ZmGOLS2 promoter was amplified by PCR from M1 using primers F8

179

and R8 and the 3’ part of the ZmGOLS2 promoter was amplified by PCR from M1 using

180

primers F9 and R9. The mutated ZmGOLS2 promoter (M3) was obtained by PCR using

181

the two PCR amplicons as templates with primers F8 and R9.

182

For expression of ZmVP1 protein used for antibody preparation, the truncated

183

coding sequence of ZmVP1 with deletion of its acid A domain was amplified by PCR

184

from the cloned vector (used for maize protoplast transformation) using primer F11 and

185

R11. For construction of bacterial expression vectors of ZmABI5, the coding sequence

186

of ZmABI5 was amplified by PCR from the cloned vector (used for maize protoplast

187

transformation) using primer F12 and R12. For both ZmVP1 and ZmABI5, the PCR

188

products were purified and inserted into a pET-28a vector using the EcoRI and HindIII 9

ACS Paragon Plus Environment

Journal of Agricultural and Food Chemistry

189

sites which results in a fusion protein with hexahistidyl tag on both N-terminal and C-

190

terminal. An additional bacterial expression vector of ZmVP1 was also constructed for

191

purification of the protein for EMSA experiments. The DNA fragment of ZmVP1

192

encoding the B2 and B3 domain was amplified by PCR from previously cloned vectors

193

using primer F13 and R13. The PCR products were purified and inserted into the

194

pGEX-4T-1 vector between EcoRI and XhoI.

195

For construction of vectors used for yeast two-hybrid assay, the coding region of

196

ZmVP1 was cloned by PCR from previously described vector using primer F14 and

197

R14. The amplicons were purified and inserted into pGADT7 vector between EcoRI

198

and SacI, which is in frame with the GAL4 DNA-activation domain. The coding region

199

of ZmABI5 was cloned by PCR from a previously described vector using primer F15

200

and R15. The PCR products were purified and inserted into the pGBKT7 vector

201

between EcoRI and SalI, which is in frame with the GAL4 DNA-binding domain.

202

For construction of vectors used for bimolecular fluorescence complementation

203

(BiFC ) assay, the coding region of ZmVP1 was amplified by PCR from previously

204

cloned vector using primer F16 and R16. The PCR products were purified and inserted

205

into PTF-486 NYFP between SpeI and XhoI. The coding region of ZmABI5 was

206

amplified by PCR from previously cloned vectors using primer F17 and R17. The PCR

207

products were purified and inserted into PTF-486 CYFP between SpeI and XhoI.

208

For construction of vectors used in the sub-cellular localization experiment, the

209

coding sequence of ZmVP1 or ZmABI5 without the stop codon was amplified using 10

ACS Paragon Plus Environment

Page 10 of 34

Page 11 of 34

Journal of Agricultural and Food Chemistry

210

PCR from a previously cloned vector (for ZmVP1, primer F18 and R18 were used; for

211

ZmABI5, primer F19 and R19 were used). The PCR products were purified and inserted

212

into pGL3-YFP vector between NheI and MluI28.

213

Antibody preparation, protein extraction and western blot

214

ZmVP1- or ZmABI5-His6 fusion proteins were expressed in BL21(DE3) E. coli

215

strain. The ZmVP1- or ZmABI5-His6 fusion proteins were purified over Nickel

216

columns and were used to immunize the rabbit following a published protocol27.

217

Western Blot analysis of ZmGOLS2 and ZmABI5 protein expression in maize embryos

218

at 20-, 26-, 32- and 36-DAP were performed following a published protocol27.

219

Maize protoplast preparation and transformation

220

The protoplast preparation and culture followed a published protocol26.

221

Preparation of GST‑ZmVP1, His-ZmABI5 fusion protein and electrophoretic

222

mobility shift assay (EMSA)

223

The bacterial expression vector pGEX-4T-1-ZmVP1 was transformed into the E.

224

coli strain BL21 (DE3), and the pET-28a-ZmABI5 was transformed into the E. coli

225

(Rosetta gami2, DE3; EMD Millipore) cells. Bacterial cultures were grown to OD600

226

nm = 0.6 at 37°C, then supplemented with IPTG (the final concentration was 0.42 mM),

227

and further cultured at 25°C for 12 hr. GST-ZmVP1 protein or GST-tag alone (used as

228

a negative control) was purified using a Glutathione Sepharose 4B column (GE

229

Healthcare) according to the manufacturer’s protocol. His-ZmABI5 protein was

230

purified using a GenScript High Affinity Ni-Charged Resin (L00250) according to the 11

ACS Paragon Plus Environment

Journal of Agricultural and Food Chemistry

231

manufacturer’s protocol. For gel mobility shift assays, the doublestranded DNA

232

containing the identified G-box and the mutated G-box were prepared by annealing two

233

complementary oligonucleotides (G-box-F and G-box-R; mG-box-F and mG-box-R;

234

(Supplemental Table 1)) following a published protocol26. The annealed, double-

235

stranded DNA was used directly for DNA EMSA. The binding reaction of DNA and

236

purified protein was conducted using an EMSA kit (Invitrogen, Carlsbad, CA, USA).

237

Mixtures of 10 ng of DNA and different amounts of purified protein were incubated in

238

a binding buffer (25 mM HEPES, 40 mM KCl, 5 mM MgCl2, 1 mM DTT, 1 mM EDTA,

239

and 8% glycerol, pH 8.0) for 30 min at room temperature and then applied to a 6% non-

240

denaturing polyacrylamide gel (29:1 acrylamide/bisacrylamide). Gel electrophoresis

241

and staining followed the kit protocol. The gel image was captured using ChemiDoc

242

MP (Bio-Rad Inc., Hercules, CA, USA).

243

Yeast two-hybrid assay

244

Yeast two-hybrid assay was performed following a published protocol28. The

245

vectors were co-transformed into yeast strain AH109 through the lithium acetate

246

method29. After 3 days culture, the transformed single colonies were picked up from

247

synthetic drop-out (SD)-Leu-Trp (LT) solid medium

248

medium. The exponentially grown yeast cells were centrifuged at 5000 ×g at room

249

temperature and adjusted to OD600 = 0.6 with sterilized double-distilled water. The yeast

250

cells were then diluted 10×, 100× and 1000× with distilled water. Two microliters of

251

the above yeast cells were spotted onto SD-LT and SD-LT-His-Ade (LTHA) medium

252

and grown at 30°C for 2 d (SD-LT plates) or for 5 d (SD-LTHA plates, without or with

and cultured in SD-LT liquid

12

ACS Paragon Plus Environment

Page 12 of 34

Page 13 of 34

Journal of Agricultural and Food Chemistry

253

30 or 60 mM 3-amino-1,2,4-triazole [AT]) before being photographed.

254

Subcellular localization and BiFC assay

255

For subcellular localization of ZmVP1 and ZmABI5 analysis, different vectors

256

and nuclear localization signal (NLS)-mCherry vector were co-transformed into maize

257

protoplasts. For the BiFC assay, vector DNA (10 μg of each) were co-transformed into

258

maize protoplasts. After culture for 12 h, YFP and NLS-mCherry signals were

259

monitored and differential interference contrast (DIC) imaging conducted using a

260

scanning confocal microscope (Andor Revolution WD, UK).

261

Statistical analysis

262

Statistical analysis was performed using SPSS 20.0. Significance (P < 0.05) was

263

assessed by using the Student’s t test or by Duncan's multiple range test.

264

Results

265

Raffinose is not detectable in zmvp1-R mutant seed.

266

The maize viviparous1 (vp1) mutant was characterized by PCR (Supplemental Fig.

267

1A). A 2 kb or a 7 kb DNA fragment was amplified from the wild type or zmvp1-R

268

plants, respectively, using primers which span the insertion site (Supplemental Fig. 1A

269

and B; sequence detail in Supplemental Table 1), while a 2 kb DNA fragment was

270

amplified from both WT and zmvp1-R plants using primers beyond the Mu element-

271

insertion site (Supplemental Fig. 1A and B). PCR and sequencing of the resultant

272

amplicons revealed a typical 5-bp duplication flanking the Mu element (Supplemental

273

Fig. 1A). The transcription of ZmVP1 was detected in WT plants, but not in zmvp1-R 13

ACS Paragon Plus Environment

Journal of Agricultural and Food Chemistry

274

plants, as determined by RT-PCR using F2a and R2 primers (Supplemental Fig. 1C).

275

Consistent with the mRNA expression, ZmVP1 protein accumulation was detected in

276

WT plants, but not in zmvp1-R plants using western blot hybridization (Supplemental

277

Fig. 1D).

278

The zmvp1-R seeds started to complete germination on the cob at 26 DAP and all

279

kernels had completed germination by 36 DAP while WT seed did not complete

280

germination precociously, even at the end of maturation (55 DAP; Fig. 1A). The

281

glucose and sucrose contents in the zmvp1-R embryo were significantly greater than in

282

the WT embryo between 20- and 36-DAP as determined by HPLC (Fig. 1B-C,

283

Supplemental Fig. 2). Myo-inositol accumulation in the zmvp1-R embryo was greater

284

than that of WT at 20- and 22-DAP, and was less than that of WT at 26-, 32-, 34- and

285

36-DAP (Fig. 1D). The raffinose and galactinol amount in WT embryos started to

286

accumulate at 28 DAP and reached to its peak at 32- or 34-DAP and then declined,

287

while they were not detected at any point in the zmvp1-R embryos (Fig. 1E-F). The

288

glucose content in the endosperm of zmvp1-R seeds was higher than that of WT at 34-

289

and 36-DAP (Supplemental Fig. 3A). The sucrose in the endosperm of zmvp1-R seeds

290

was less than that of WT at 32- and 34-DAP, but was higher than that of WT at 36 DAP

291

(Supplemental Fig. 3B). Myo-inositol, galactinol and raffinose content of the

292

endosperm were all below the detection limit in both WT and the zmvp1-R mutant

293

(Supplemental Fig. 4).

294

ZmGOLS2 is down-regulated in the zmvp1-R embryo

295

GALACTINOL SYNTHASE and RAFFINOSE SYNTHASE are two enzymes 14

ACS Paragon Plus Environment

Page 14 of 34

Page 15 of 34

Journal of Agricultural and Food Chemistry

296

that are responsible for raffinose biosynthesis (Fig. 2A). To investigate the molecular

297

mechanism underlying the absence of raffinose in the zmvp1-R embryo, the

298

transcriptome was compared between WT and zmvp1-R developing embryos at 26 DAP,

299

a time point prior to raffinose accumulation but when GALACTINOL SYNTHASE2 and

300

RAFFINOSE SYNTHASE transcript and ZmGOLS2 protein is present in considerable

301

abundance (Fig. 2; Supplemental Fig. 5; Supplemental Table 3). The expression of both

302

ZmGOLS2 and ZmRAFS was significantly decreased in the zmvp1-R embryo as

303

compared with that of WT (Supplemental Table 3). The expression of ZmGOLS2 or

304

ZmRAFS in the zmvp1-R embryo was significantly less than that of WT at 20-, 26-, 32-

305

and 36-DAP as determined by quantitative RT-PCR (Fig. 2B; Supplemental Fig. 5).

306

The ZmGOLS2 protein accumulation in the zmvp1-R embryo at 26-, 32- or 36-DAP

307

was less than that of WT as determined by western blot analysis (Fig. 2C).

308

The predicted cis-elements in the 5’ regulatory region of ZmGOLS2 are depicted

309

in the sequence provided (Supplemental Fig. 6). There was no RY motif in the

310

ZmGOLS2 promoter (Supplemental Fig. 6), suggesting that ZmVP1 may not directly

311

bind to the ZmGOLS2 promoter as this is the known element to which ZmVP1 binds30.

312

ABI5 binds to the G-box or ABRE motif in the promoter of its target genes17. There are

313

3 G-boxes and 2 ABRE elements in the promoter of ZmGOLS2 (Supplemental Fig. 6),

314

suggesting that ABI5 may directly regulate ZmGOLS2 expression.

315

ZmVP1 and ZmABI5 synergistically up-regulate the expression of ZmGOLS2.

316

GFP, ZmVP1 or ZmABI5 expression vectors were constructed (Fig. 3A). The

317

ZmABI5 or ZmVP1/ZmABI5 expression vectors were co-transformed with a variety of 15

ACS Paragon Plus Environment

Journal of Agricultural and Food Chemistry

318

dual luciferase expression vectors into maize protoplasts (Fig. 3A). In these

319

experiments, the Rluc expression in the dual luciferase expression vector was controlled

320

by the ZmGOLS2 promoter (V6), or by a mutated ZmGOLS2 promoter in which one

321

ABRE motif (M1) or one G-box (M2), or both ABRE and G-box (M3) were mutated

322

(Fig. 3A). Rluc expression, regardless of ZmGOLS2 promoter, was not changed by

323

simultaneous ZmVP1 expression (Fig. 3B). However, Rluc expression, when controlled

324

by the V6 or the M1 fragment, was significantly increased by simultaneous ZmABI5

325

over-expression and this was further and dramatically increased by the co-expression

326

of ZmVP1/ZmABI5, as compared with the GFP expression (Fig. 3B). Rluc expression

327

controlled by the M2 or M3 fragment, was not changed by simultaneous expression of

328

ZmVP1, ZmABI5 or ZmVP1/ZmABI5 as compared with the GFP expression (Fig. 3B).

329

These data show that the G-box in the ZmGOLS2 promoter is functional in regulation

330

of its expression in response to ZmABI5. Furthermore, endogenous ZmGOLS2

331

expression and protein abundance, while not enhanced by ZmVP1 expression, was

332

significantly increased in maize protoplasts expressing ZmABI5 and this was further

333

augmented by the co-expression of ZmVP1/ZmABI5 as determined by real time RT-

334

PCR and western blot analysis, as compared with the GFP expression (Fig. 3 C-E).

335

DNA-EMSA results showed that purified ZmABI5 protein was able to bind to the G-

336

box elements of the ZmGOLS2 promoter in vitro, while the purified ZmVP1 protein, or

337

the GST tag failed to bind to the G-box element of the ZmGOLS2 promoter (Fig. 3F).

338

ZmVP1 interacts with ZmABI5

339

There are three G-box motifs in the ZmGOLS2 promoter (Supplemental Fig. 6). It 16

ACS Paragon Plus Environment

Page 16 of 34

Page 17 of 34

Journal of Agricultural and Food Chemistry

340

has been reported that bZIP transcription factors are capable of binding G-box elements

341

due to the core ABRE motif they contain31 and regulate their target genes’ expression.

342

There are 125 bZIP genes predicted in the maize genome18. By searching the published

343

maize embryo and endosperm transcriptome19, we found the expression of 9 bZIP genes

344

that were positively correlated with the expression of ZmGOLS2 in maize embryos

345

(Supplemental Table 4). Phylogenetic analysis of these 9 maize bZIP transcription

346

factors and the known ABI5s from other species showed that GRMZM2G168079 is

347

most closely related to OsABI5 and AtABI5 (Fig. 4A). This GRMZM2G168079

348

protein was named ZmABI5. The mRNA and protein expression of ZmABI5 is

349

significantly decreased in the zmvp1-R embryo as determined by RT-PCR and western

350

blot analysis (Fig. 4B-C).

351

To validate the expected nuclear sub-cellular localization of the two transcription

352

factors, ZmVP1-YFP and ZmABI5-YFP fusion protein expression vectors were

353

constructed and separately transformed into maize protoplasts (Supplemental Fig. 7).

354

The YFP signal was monitored by confocal microscopy to determine protein subcellular

355

residency. The signal of both ZmVP1-YFP and ZmABI5-YFP was detected in the

356

nucleus as anticipated (Supplemental Fig. 7).

357

bZIP transcription factors are also known to interact with ABI3, an Arabidopsis

358

thaliana homolog of VP132. Yeast two-hybrid assay showed that ZmVP1 and ZmABI5

359

interacted with each other (Fig. 5A). The interaction of ZmVP1 and ZmABI5 was

360

further confirmed using BiFC in maize protoplasts (Fig. 5B). The YFP signal was not

361

detected in the protoplasts transformed with control constructs but was detected in 17

ACS Paragon Plus Environment

Journal of Agricultural and Food Chemistry

362 363

protoplasts transformed with ZmVP1-YN and ZmABI5-YC (Fig. 5B). Taken together, our data showed that ZmVP1 interacts with ZmABI5 and

364

regulates ZmGOLS2 expression and raffinose biosynthesis in maize seeds (Fig. 6).

365

Discussion

366

The regulation of raffinose biosynthesis

367

GOLS and RAFS are two key enzymes for raffinose biosynthesis (Fig. 2). The

368

genes encoding GOLS are induced by environmental stresses in plants. Over-

369

expression of GOLS genes isolated from different species increased the amounts of both

370

galactinol and raffinose and improved drought tolerance in transgenic plants27, 33. These

371

results suggest that GALACTINOL SYNTHASE may be a bottleneck for raffinose

372

accumulation. The consistent benefits on plant abiotic stress tolerance of GOLS

373

hyperproduction through transgenic means warrants an examination of how

374

endogenous GOLS is regulated in order to embark on more nuanced manipulations of

375

GOLS transcription, ensuing galactinol and raffinose production, and the level of

376

protection this affords crop productivity in stressful conditions. A better understanding

377

of the various elements acting to boost transcription of GOLS will aid interpretation of

378

intraspecific, natural variability in galactinol and raffinose production upon the

379

induction of stress. To this end, we have found that ZmGOLS2 is induced by

380

dehydration or heat stress5, 26, 34. In addition, we identified that ZmDREB2A regulates

381

the expression of ZmGOLS2 by binding to the DRE elements of the ZmGOLS2

382

promoter in maize leaf protoplasts27, suggesting that dehydration stress or ZmDREB2A

383

up-regulates raffinose biosynthesis in plant leaves. A WRKY transcription factor has 18

ACS Paragon Plus Environment

Page 18 of 34

Page 19 of 34

Journal of Agricultural and Food Chemistry

384

also been characterized to bind to the W-box motif of the Boea hygrometrica

385

GALACTINOL SYNTHASE promoter and regulate its expression35. ZmVP1 is

386

predominantly expressed in maize seed9. ZmVP1 positively regulates ZmGOLS2

387

expression and does so dependent upon, and synergistically with, ZmABI5 (Fig. 3).

388

These data suggest that ZmGOLS2 is regulated by different transcription factors under

389

different spatio-temporal conditions. It now appears as though ZmVP1 and ZmABI5

390

are key regulators for raffinose biosynthesis in the maize embryo.

391

Multiple transcription factors, such as bZIPs, B3s, MYBs, and DOFs, function in

392

the zygotic embryo during seed maturation36. The bZIP53 protein specifically binds to

393

the G-box element in the promoter of ALBUMIN 2S2 gene and the formation of the

394

heterodimer of bZIP53 with bZIP10 or bZIP25, significantly enhances the dimer’s

395

DNA binding activity and synergistically increases its target gene expression over that

396

of bZIP53 alone36. The bZIP53 does not directly interact with ABI3, but the ternary

397

complex formation between the bZIP heterodimers and ABI3 increases target gene

398

expression still further36. Our results showed that ZmABI5 directly interacted with

399

ZmVP1 (Fig. 5) and co-expression of ZmABI5 and ZmVP1 in maize protoplasts

400

synergistically enhanced ZmGOLS2 expression (Fig. 3). Our data cannot exclude the

401

possibility that another bZIP transcription factor, up-regulated by ZmVP1, may

402

participate to form a heterodimer with ZmABI5 and that this is the cause of the observed

403

ZmVP1 stimulation of ZmABI5 transcriptional control of ZmGOLS2. If there was a

404

heterodimer of ZmABI5 and another bZIP factor, it was not required for the interaction

405

of ZmABI5 and ZmVP1 in yeast where the existence of an endogenous bZIP 19

ACS Paragon Plus Environment

Journal of Agricultural and Food Chemistry

406

transcription factor capable of mediating the ZmABI5-ZmVP1 interaction is

407

improbable. The more parsimonious explanation of our results is that ZmABI5 binds

408

the G-box most proximal to the ZmGOLS2 transcriptional start site and that ZmVP1

409

binds ZmABI5 to boost the transcriptional stimulation of ZmGOLS2.

410

RAFFINOSE SYNTHASE (ZmRAFS) is the key enzyme for raffinose

411

biosynthesis. However, to date, the transcriptional regulation of ZmRAFS has not been

412

characterized. Though ZmRAFS mRNA expression was significantly decreased in the

413

zmvp1-R embryo (Supplemental Fig. 5), using the same strategy used to characterize

414

the transcriptional regulation of ZmGOLS2, we found that the ZmRAFS gene was not

415

regulated by ZmVP1, ZmABI5, or their dimer in maize protoplasts (Supplemental Fig.

416

8), suggesting that the regulation of ZmRAFS deviates considerably from that of

417

ZmGOLS2.

418

The dependence of raffinose biosynthesis on the ABA signaling pathway

419

The regulation of raffinose biosynthesis appears to be independent of the abscisic

420

acid (ABA) signaling pathway, at least during dehydration stress in Arabidopsis

421

thaliana37. ZmVP1, like its Arabidopsis thaliana orthologue, ABI3, is responsible for

422

the plant cell to react to ABA9 and ABI3 in Arabidopsis thaliana, regulates the

423

expression of ABI516. Similarly, the expression of ZmABI5 was down-regulated in the

424

zmvp1-R mutant (Fig. 4B, C). VP1/ABI3 controls target gene expression by binding to

425

the RY/Sph element through its B3 domain while ZmABI5 does so by binding to G-

426

box (ABRE) elements. A synopsis of our results follows: in the absence of a functional

427

VP1, neither galactinol nor raffinose accumulate in the maize embryo (Fig. 1), 20

ACS Paragon Plus Environment

Page 20 of 34

Page 21 of 34

Journal of Agricultural and Food Chemistry

428

ZmGOLS2, ZmRAFS, and ZmABI5 transcripts are all repressed below WT in embryos

429

(Fig. 2B; Fig. 4B, C; Supplemental Fig. 5), and ZmGOLS2 and ZmABI5 protein

430

amounts are also much below WT (Fig. 2C, 4C); yet there are no RY/Sph motifs in the

431

ZmGOLS2 promoter to which VP1 could bind while there are three G-box motifs

432

(Supplemental Fig. 6); overexpression of ZmVP1 in maize protoplasts failed to up-

433

regulate the expression of reporter or endogenous ZmGOLS2 or of endogenous

434

ZmABI5 protein (Fig. 3); using the Arabidopsis thaliana ABI3/ABI5 interaction

435

paradigm we identified the maize bZIP most similar to ABI5 and named this putative

436

ortholog ZmABI5 (Fig. 4A); ZmABI5 transcript and protein abundance is severely

437

reduced in the zmvp1 embryo (Fig. 4B, C); ZmVP1 interacts with ZmABI5 in both yeast

438

and in maize protoplasts subjected to BiFC assays (Fig. 5); ZmABI5 overexpression in

439

maize protoplasts results in up-regulation of ZmGOLS2; the hyper-production of

440

ZmVP1, when ZmABI5 is also present in superabundance, synergistically increased

441

ZmGOLS2 expression and protein amounts above that of ZmABI5 alone (Fig. 3). Taken

442

together, these data suggest that ZmVP1 not only regulates the expression of ZmABI5,

443

at least in developing embryos, but also interacted with ZmABI5 to enhance its

444

transcriptional activation of at least ZmGOLS2. There is no difference in AtGOLS1 and

445

AtGOLS2 expression or RFO accumulation between the ABA deficient mutant atnced3

446

and WT under dehydration stress38, suggesting that galactinol biosynthesis, at least

447

from these two AtGOLS family members, was regulated by ABA-independent

448

pathways. Conflicting with this finding, the ABA-deficient and ABA-insensitive

449

aba1:abi3-1 double mutant seeds are viable but with reduced desiccation tolerance; 21

ACS Paragon Plus Environment

Journal of Agricultural and Food Chemistry

450

supplementation of ABA to the aba1:abi3-1 double mutant seed increased RFO

451

biosynthesis and promoted desiccation tolerance39. Raffinose biosynthesis was

452

stimulated by abiotic stress, such as drought-, salt-, or heat-stress, in plant vegetative

453

tissues. This type of regulation is mediated by some transcription factors, such as

454

DREB2A, and can be ABA-independent27. ABA promotes seed maturation and inhibits

455

completion of seed germination in many plant species40. During maize seed maturation,

456

the accumulation of raffinose is regulated by VP1 (Fig. 1), suggesting that the presence

457

of ABA stimulates the biosynthesis of oligosaccharides while suppressing precocious

458

germination.

459

The soluble oligosaccharides in seed and their possible functions

460

The absolute and relative abundance of different soluble oligosaccharides in seeds

461

determines the seed physiological status41. The conversion of the small, soluble

462

carbohydrate pool from predominately reducing monosaccharides to predominately

463

non-reducing sugar oligomers during orthodox seed maturation reduces the absolute

464

number of these dissolved solutes, protects against Maillard reactions by tying up

465

reducing ends in inter-monosaccharide covalent bonds, while creating oligosaccharides

466

capable of forming a glassy state conducive to the intracellular milieu withstanding

467

desiccation42. As well, the relative abundance of sucrose to RFO appears to influence

468

seed vigor2. The zmvp1-R mutant seed exhibits increased glucose and sucrose and

469

decreased galactinol and raffinose, compared with that of WT, and showed a viviparous

470

phenotype (Fig. 1B-D). This phenomenon was also observed in Arabidopsis thaliana

471

abi3-5 mutant seeds43. Raffinose accumulates in the late stages of seed development 22

ACS Paragon Plus Environment

Page 22 of 34

Page 23 of 34

Journal of Agricultural and Food Chemistry

472

and is hydrolyzed during seed germination2, 44. The maize zmvp1-R mutant seed and the

473

Arabidopsis thaliana mutant lines (abi3-5, lec1-1 and fus3-3)45 with reduced raffinose,

474

exhibited a desiccation-intolerant phenotype (a drastically curtailed longevity),

475

supporting the previous findings that the ratio of raffinose to sucrose determines seed

476

longevity in maize seeds2. RFOs are considered to be anti-nutritional46-47 while, on the

477

contrary, RFOs are an important energy resource for beneficial microflora in the small

478

intestine48. Understanding the regulatory mechanism for raffinose biosynthesis will aid

479

interpretation of intraspecific, natural variability in galactinol and raffinose production

480

upon the induction of stress, providing plant breeders with a more sophisticated model

481

of whole-genome alterations directed at using this biochemical pathway to combat

482

productivity loss due to plant stress. At a more fundamental level, an understanding of

483

the precise physical function of RFO-mediated protection of stress susceptible cellular

484

constituents will serve to pinpoint which entities in the cell are most susceptible to stress

485

while demonstrating one means by which plants protect these entities from stress as a

486

paradigm on which to focus synthetic biological research.

487

Supporting Information

488

The information about ZmVP1 mutant characterization, ZmRAFS expression and

489

regulation, Analysis of ZmGOLS2 promoter, HPLC analysis of endosperm sugar profile,

490

subcellular localization of ZmVP1 and ZmABI5, the RNAseq data, the primer list and

491

the abbreviation list were provided as supplemental data. The Supporting Information

492

is available free of charge on the ACS Publications website.

493

Author Information 23

ACS Paragon Plus Environment

Journal of Agricultural and Food Chemistry

494

Corresponding Author Tianyong Zhao; Email: [email protected]; ORCID ID:

495

0000-0002-1278-2842; Co-authors: Lynnette M.A. Dirk ORCID ID: 0000-0002-7564-

496

5095; A. Bruce Downie ORCID ID: 0000-0001-6680-0551.

497

Funding

498

This research was funded by the National Key Research and Development Plan

499

(2018YFD0100901) and the NSFC (31671776) to TZ. For LMAD and ABD, this work

500

was supported by the USDA National Institute of Food and Agriculture, Hatch project

501

1002299.

502

Acknowledgements

503

The HPLC experiments were conducted in The Biology Teaching and Research Core

504

Facility at College of Life Sciences, Northwest A&F University.

505

Notes

506

The authors declare no conflicts of interest. All the abbreviations are listed in

507

Supplemental Table 2.

508

References

509 510 511

1.

512 513 514

2.

515 516

3.

517 518 519

4.

Handley, L. W.; Pharr, D. M.; McFeeters, R. F., Carbohydrate Changes during Maturation of

Cucumber Fruit : Implications for Sugar Metabolism and Transport. Plant Physiol 1983, 72 (2), 498502. Li, T.; Zhang, Y.; Wang, D.; Liu, Y.; Dirk, L. M. A.; Goodman, J.; Downie, A. B.; Wang, J.; Wang, G.;

Zhao, T., Regulation of Seed Vigor by Manipulation of Raffinose Family Oligosaccharides in Maize and Arabidopsis thaliana. Mol Plant 2017, 10 (12), 1540-1555. Saravitz, D. M.; Pharr, D. M.; Carter, T. E., Galactinol synthase activity and soluble sugars in

developing seeds of four soybean genotypes. Plant Physiol 1987, 83 (1), 185-9. Lehle, L.; Tanner, W., The function of myo-inositol in the biosynthesis of raffinose. Purification

and characterization of galactinol:sucrose 6-galactosyltransferase from Vicia faba seeds. Eur J Biochem 1973, 38 (1), 103-10. 24

ACS Paragon Plus Environment

Page 24 of 34

Page 25 of 34

Journal of Agricultural and Food Chemistry

520 521 522

5.

523 524

6.

525 526

7.

527 528

8.

529 530 531

9.

532 533

10. McCarty, D. R.; Carson, C. B.; Stinard, P. S.; Robertson, D. S., Molecular Analysis of viviparous-1:

534 535

11. Brenac, P.; Smith, M. E.; Obendorf, R. L., Raffinose accumulation in maize embryos in the

536 537 538

12. Ooms, J.; Leon-Kloosterziel, K. M.; Bartels, D.; Koornneef, M.; Karssen, C. M., Acquisition of

539 540 541

13. Hoecker, U.; Vasil, I. K.; McCarty, D. R., Integrated control of seed maturation and germination

Zhao, T. Y.; Thacker, R.; Corum, J. W.; Snyder, J. C.; Meeley, R. B.; Obendorf, R. L.; Downie, B.,

Expression of the maize GALACTINOL SYNTHASE gene family: (I) Expression of two different genes during seed development and germination. Physiologia plantarum 2004, 121 (4), 634-646. Jang, J. H.; Shang, Y.; Kang, H. K.; Kim, S. Y.; Kim, B. H.; Nam, K. H., Arabidopsis galactinol

synthases 1 (AtGOLS1) negatively regulates seed germination. Plant Sci 2018, 267, 94-101. McCarty, D. R.; Hattori, T.; Carson, C. B.; Vasil, V.; Lazar, M.; Vasil, I. K., The Viviparous-1

developmental gene of maize encodes a novel transcriptional activator. Cell 1991, 66 (5), 895-905. Suzuki, M.; Kao, C. Y.; Cocciolone, S.; McCarty, D. R., Maize VP1 complements Arabidopsis abi3

and confers a novel ABA/auxin interaction in roots. Plant Journal 2001, 28 (4), 409-418. Cao, X.; Costa, L. M.; Biderre-Petit, C.; Kbhaya, B.; Dey, N.; Perez, P.; McCarty, D. R.; Gutierrez-

Marcos, J. F.; Becraft, P. W., Abscisic acid and stress signals induce Viviparous1 expression in seed and vegetative tissues of maize. Plant Physiol 2007, 143 (2), 720-31.

An Abscisic Acid-Insensitive Mutant of Maize. Plant Cell 1989, 1 (5), 523-532.

absence of a fully functional Vp1 gene product. Planta 1997, 203 (2), 222-228.

Desiccation Tolerance and Longevity in Seeds of Arabidopsis thaliana (A Comparative Study Using Abscisic Acid-Insensitive abi3 Mutants). Plant Physiol 1993, 102 (4), 1185-1191.

programs by activator and repressor functions of Viviparous-1 of maize. Genes Dev 1995, 9 (20), 245969.

542 543 544 545

14. Delmas, F.; Sankaranarayanan, S.; Deb, S.; Widdup, E.; Bournonville, C.; Bollier, N.; Northey, J.

546 547 548

15. Zinsmeister, J.; Lalanne, D.; Terrasson, E.; Chatelain, E.; Vandecasteele, C.; Vu, B. L.;

549 550 551

16. Lopez-Molina, L.; Mongrand, B.; McLachlin, D. T.; Chait, B. T.; Chua, N. H., ABI5 acts downstream

G. B.; McCourt, P.; Samuel, M. A., ABI3 controls embryo degreening through Mendel's I locus. Proceedings of the National Academy of Sciences of the United States of America 2013, 110 (40), E3888-E3894.

Duboislaurent, C.; Geoffriau, E.; Signor, C. L.; Dalmais, M., ABI5 Is a Regulator of Seed Maturation and Longevity in Legumes. Plant Cell 2016, 28 (11), 2735.

of ABI3 to execute an ABA-dependent growth arrest during germination. Plant Journal 2002, 32 (3), 317-328.

552 553 554

17. Carles, C.; Bies-Etheve, N.; Aspart, L.; Leon-Kloosterziel, K. M.; Koornneef, M.; Echeverria, M.;

555 556

18. Wei, K.; Chen, J.; Wang, Y.; Chen, Y.; Chen, S.; Lin, Y.; Pan, S.; Zhong, X.; Xie, D., Genome-wide

Delseny, M., Regulation of Arabidopsis thaliana Em genes: role of ABI5. Plant Journal 2002, 30 (3), 373-383.

analysis of bZIP-encoding genes in maize. DNA Res 2012, 19 (6), 463-76. 25

ACS Paragon Plus Environment

Journal of Agricultural and Food Chemistry

557 558

19. Chen, J.; Zeng, B.; Zhang, M.; Xie, S. J.; Wang, G. K.; Hauck, A.; Lai, J. S., Dynamic Transcriptome

559 560 561

20. Thompson, J. D.; Higgins, D. G.; Gibson, T. J., Clustal-W - Improving the Sensitivity of Progressive

562 563

21. Kumar, S.; Stecher, G.; Tamura, K., MEGA7: Molecular Evolutionary Genetics Analysis Version

564 565

22. Saitou, N.; Nei, M., The neighbor-joining method: a new method for reconstructing phylogenetic

566 567 568

23. Porebski, S.; Bailey, L. G.; Baum, B. R., Modification of a CTAB DNA extraction protocol for plants

569 570

24. Jaakola, L.; Pirttila, A. M.; Halonen, M.; Hohtola, A., Isolation of high quality RNA from bilberry

571 572

25. Higo, K.; Ugawa, Y.; Iwamoto, M.; Korenaga, T., Plant cis-acting regulatory DNA elements

573 574

26. Gu, L.; Han, Z.; Zhang, L.; Downie, B.; Zhao, T., Functional analysis of the 5' regulatory region of

575 576 577

27. Gu, L.; Zhang, Y. M.; Zhang, M. S.; Li, T.; Dirk, L. M. A.; Downie, B.; Zhao, T. Y., ZmGOLS2, a target

578 579 580

28. Zhang, Y. M.; Linghu, J. J.; Wang, D.; Liu, X.; Yu, A. L.; Li, F. T.; Zhao, J. F.; Zhao, T. Y., Foxtail Millet

581 582

29. Gietz, R. D.; Schiestl, R. H.; Willems, A. R.; Woods, R. A., Studies on the transformation of intact

583 584

30. Suzuki, M., .; Kao, C. Y.; Mccarty, D. R., The conserved B3 domain of VIVIPAROUS1 has a

585 586

31. Sarkar, A. K.; Lahiri, A., Specificity determinants for the abscisic acid response element. Febs

Landscape of Maize Embryo and Endosperm Development. Plant physiology 2014, 166 (1), 252-264.

Multiple Sequence Alignment through Sequence Weighting, Position-Specific Gap Penalties and Weight Matrix Choice. Nucleic acids research 1994, 22 (22), 4673-4680.

7.0 for Bigger Datasets. Molecular biology and evolution 2016, 33 (7), 1870-1874.

trees. Mol Biol Evol 1987, 4 (4), 406-25.

containing high polysaccharide and polyphenol components. Plant Molecular Biology Reporter 1997, 15 (1), 8-15.

(Vaccinium myrtillus L.) fruit. Mol Biotechnol 2001, 19 (2), 201-3.

(PLACE) database: 1999. Nucleic Acids Res 1999, 27 (1), 297-300.

the maize GALACTINOL SYNTHASE2 gene. Plant Sci 2013, 213, 38-45.

of transcription factor ZmDREB2A, offers similar protection against abiotic stress as ZmDREB2A. Plant molecular biology 2016, 90 (1-2), 157-170.

CBL4 (SiCBL4) Interacts with SiCIPK24, Modulates Plant Salt Stress Tolerance. Plant Molecular Biology Reporter 2017, 35 (6), 634-646.

yeast cells by the LiAc/SS-DNA/PEG procedure. Yeast 1995, 11 (4), 355-60.

cooperative DNA binding activity. Plant Cell 1997, 9 (5), 799-807.

Open Bio 2013, 3, 101-105.

587 588

32. Nakamura, S.; Lynch, T. J.; Finkelstein, R. R., Physical interactions between ABA response loci of

589 590 591

33. Taji, T.; Ohsumi, C.; Iuchi, S.; Seki, M.; Kasuga, M.; Kobayashi, M.; Yamaguchi-Shinozaki, K.;

592 593

34. Zhao, T. Y.; Martin, D.; Meeley, R. B.; Downie, B., Expression of the maize GALACTINOL

Arabidopsis. Plant Journal 2001, 26 (6), 627-635.

Shinozaki, K., Important roles of drought- and cold-inducible genes for galactinol synthase in stress tolerance in Arabidopsis thaliana. Plant J 2002, 29 (4), 417-26.

SYNTHASE gene family: (II) Kernel abscission, environmental stress and myo-inositol influences 26

ACS Paragon Plus Environment

Page 26 of 34

Page 27 of 34

Journal of Agricultural and Food Chemistry

594 595

accumulation of transcript in developing seeds and callus cells. Physiologia Plantarum 2004, 121 (4),

596 597 598

35. Wang, Z.; Zhu, Y.; Wang, L.; Liu, X.; Liu, Y.; Phillips, J.; Deng, X., A WRKY transcription factor

599 600 601 602

36. Alonso, R.; Onate-Sanchez, L.; Weltmeier, F.; Ehlert, A.; Diaz, I.; Dietrich, K.; Vicente-Carbajosa,

603 604 605 606

37. Urano, K.; Maruyama, K.; Ogata, Y.; Morishita, Y.; Takeda, M.; Sakurai, N.; Suzuki, H.; Saito, K.;

647-655.

participates in dehydration tolerance in Boea hygrometrica by binding to the W-box elements of the galactinol synthase (BhGolS1) promoter. Planta 2009, 230 (6), 1155-66.

J.; Droge-Laser, W., A pivotal role of the basic leucine zipper transcription factor bZIP53 in the regulation of Arabidopsis seed maturation gene expression based on heterodimerization and protein complex formation. Plant Cell 2009, 21 (6), 1747-61.

Shibata, D.; Kobayashi, M.; Yamaguchi-Shinozaki, K.; Shinozaki, K., Characterization of the ABAregulated global responses to dehydration in Arabidopsis by metabolomics. Plant J 2009, 57 (6), 106578.

607 608 609 610

38. Urano, K.; Maruyama, K.; Ogata, Y.; Morishita, Y.; Takeda, M.; Sakurai, N.; Suzuki, H.; Saito, K.;

611 612

39. Ooms, J. J. J.; Wilnier, J. A.; Karssen, C. M., Carbohydrates are not the sole factor determining

613 614

40. Nambara, E.; Okamoto, M.; Tatematsu, K.; Yano, R.; Seo, M.; Kamiya, Y., Abscisic acid and the

615 616 617

41. Bentsink, L.; Alonso-Blanco, C.; Vreugdenhil, D.; Tesnier, K.; Groot, S. P. C.; Koornneef, M.,

618 619 620

42. Santos, J. F.; Dirk, L. M. A.; Downie, A. B.; Sanches, M. F. G.; Vieira, R. D., Reciprocal effect of

621 622 623

43. Ooms, J. J. J.; Leonkloosterziel, K. M.; Bartels, D.; Koornneef, M.; Karssen, C. M., Acquisition of

624 625 626

44. Kaczmarska, K. T.; Chandra-Hioe, M. V.; Zabaras, D.; Frank, D.; Arcot, J., Effect of Germination

627 628 629 630

45. Gonzalez-Morales, S. I.; Chavez-Montes, R. A.; Hayano-Kanashiro, C.; Alejo-Jacuinde, G.; Rico-

41.

631

46. Dierking, E. C.; Bilyeu, K. D., Association of a Soybean Raffinose Synthase Gene with Low

Shibata, D.; Kobayashi, M.; Yamaguchi-Shinozaki, K.; Shinozaki, K., Characterization of the ABAregulated global responses to dehydration in Arabidopsis by metabolomics. Plant Journal 2009, 57 (6), 1065-1078.

desiccation tolerance in seeds of Ambidopsis thaliana. Physiologia plantarum 1994, 90 (3), 431–436.

control of seed dormancy and germination. Seed Science Research 2010, 20 (2), 55-67.

Genetic analysis of seed-soluble oligosaccharides in relation to seed storability of Arabidopsis. Plant physiology 2000, 124 (4), 1595-1604.

parental lines on the physiological potential and seed composition of corn hybrid seeds. Seed Science Research 2017, 27 (3), 206-216.

Desiccation Tolerance and Longevity in Seeds of Arabidopsis-Thaliana - a Comparative-Study Using Abscisic Acid-Insensitive Abi3 Mutants. Plant physiology 1993, 102 (4), 1185-1191.

and Fermentation on Carbohydrate Composition of Australian Sweet Lupin and Soybean Seeds and Flours. J Agr Food Chem 2017, 65 (46), 10064-10073.

Cambron, T. Y.; de Folter, S.; Herrera-Estrella, L., Regulatory network analysis reveals novel regulators of seed desiccation tolerance in Arabidopsis thaliana. Proc Natl Acad Sci U S A 2016, 113 (35), E5232-

27

ACS Paragon Plus Environment

Journal of Agricultural and Food Chemistry

632

Raffinose and Stachyose Seed Phenotype. Plant Genome 2008, 1 (2), 135-145.

633 634 635

47. Bilyeu, K. D.; Wiebold, W. J., Environmental Stability of Seed Carbohydrate Profiles in Soybeans Containing Different Alleles of the Raffinose Synthase 2 (RS2) Gene. J Agr Food Chem 2016, 64 (5), 1071-1078.

636 637 638

2008, 14 (4), 205-8.

639

Figure Captions

640

Figure 1. Raffinose accumulates in wild type W22 maize seed but was not detectable

641

in the zmvp1-R mutant seed during seed maturation.

642

(A) Zmvp1 mutant seed germinated on the cob during seed development stage and its

643

morphology can be distinguished from wild type W22. DAP: days after pollination.

644

Bar = 5 mm. (B-F) Comparison of different sugar contents in embryos between WT

645

and the zmvp1-R mutant. B: Glucose; C: Sucrose; D: myo-inositol; E: Galactinol; F:

646

Raffinose. (µmol•g-1 DW: micromolar per gram dry weight). Values are means ± SE

647

(n=3). **p