Adaptive Evolution Relieves Nitrogen Catabolite Repression and

Department of Bioengineering, School of Food Sciences and Biotechnology, Zhejiang Gongshang University , Hangzhou 310018 , ... Publication Date (Web):...
0 downloads 0 Views 878KB Size
Subscriber access provided by Kaohsiung Medical University

Food and Beverage Chemistry/Biochemistry

Adaptive Evolution Relieves Nitrogen Catabolite Repression and Decreases Urea Accumulation in Cultures of the Chinese Rice Wine Yeast Strain Saccharomyces cerevisiae XZ-11 Weiping Zhang, Yan Cheng, Yudong Li, Guocheng Du, Guangfa Xie, Huijun Zou, Jingwen Zhou, and Jian Chen J. Agric. Food Chem., Just Accepted Manuscript • DOI: 10.1021/acs.jafc.8b01313 • Publication Date (Web): 08 Jun 2018 Downloaded from http://pubs.acs.org on June 8, 2018

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 36

Journal of Agricultural and Food Chemistry

1

Adaptive Evolution Relieves Nitrogen Catabolite Repression and Decreases Urea

2

Accumulation in Cultures of the Chinese Rice Wine Yeast Strain Saccharomyces

3

cerevisiae XZ-11

4 5

Weiping Zhang†,‡,#, Yan Cheng†, Yudong Li§, Guocheng Du†, Guangfa Xie¶, Huijun

6

Zouǁ, Jingwen Zhou*,†,‡,#, and Jian Chen*,†,#

7 8



9

Biotechnology, Jiangnan University, 1800 Lihu Road, Wuxi, Jiangsu 214122, China;

Key Laboratory of Industrial Biotechnology, Ministry of Education and School of

10



11

University, 1800 Lihu Road, Wuxi, Jiangsu 214122, China

12

§

13

Zhejiang Gongshang University, Hangzhou 310018, China

14



15

China

16

#

17

Jiangnan University, 1800 Lihu Road, Wuxi, Jiangsu 214122, China.

18

ǁ

19

Zhejiang, China.

20

*

21

Jian Chen, Jingwen Zhou

22

Mailing address: School of Biotechnology, Jiangnan University, 1800 Lihu Road,

23

Wuxi, Jiangsu 214122, China.

24

Phone: +86-510-85918312, Fax: +86-510-85918309

25

E-mail: [email protected], [email protected].

National Engineering Laboratory for Cereal Fermentation Technology, Jiangnan

Department of Bioengineering, School of Food Sciences and Biotechnology,

College of Shaoxing Rice Wine, Zhejiang Shuren University, Shaoxing 312028,

Jiangsu Provisional Research Center for Bioactive Product Processing Technology,

Zhejiang Guyuelongshan Shaoxing Wine Company, 13 Yangjiang Road, Shaoxing,

Corresponding authors.

1 ACS Paragon Plus Environment

Journal of Agricultural and Food Chemistry

26

ABSTRACT

27

Urea is the major precursor of ethyl carbamate in Chinese rice wine. Although efforts

28

have been made to decrease urea accumulation, few methods can be applied to the

29

industrial food production due to potential safety concerns. In this study, adaptive

30

laboratory evolution (ALE) followed by high-throughput screening was used to

31

identify low urea-accumulating strains derived from the industrial Chinese rice wine

32

yeast strain Saccharomyces cerevisiae XZ-11. Three evolved strains were obtained

33

that had 47.9%, 16.6%, and 12.4% lower urea concentrations than the wild-type strain.

34

Comparative genomics analysis revealed that genes involved in carbon and nitrogen

35

metabolism evolved quickly. Transcription levels of genes involved in urea

36

metabolism were dramatically upregulated after ALE. This work describes a novel

37

and safe strategy to improve nitrogen utilization of industrial yeast strains involved in

38

food fermentation. The identified genomic variations may also help direct rational

39

genetic engineering of nitrogen metabolism processes to achieve other goals.

40 41

Keywords: Adaptive laboratory evolution; Comparative genomics; High-throughput

42

screening; RT-qPCR; Saccharomyces cerevisiae; Urea

43

2 ACS Paragon Plus Environment

Page 2 of 36

Page 3 of 36

Journal of Agricultural and Food Chemistry

44

INTRODUCTION

45

Chinese rice wine is an important traditional alcoholic beverage in China. It has

46

been widely consumed for the past 5,000 years because of its unique flavor and high

47

nutritional value.1-2 However, harmful components have been detected in Chinese rice

48

wine, of which ethyl carbamate (EC) is the most common.3 EC, also known as

49

urethane, is naturally formed in fermented foods and has been detected in many kinds

50

of alcoholic beverages.3 In 2007, the carcinogen classification of EC was upgraded

51

from Group 2B to Group 2A by the International Agency for Research on Cancer

52

(IARC).4-5 EC is formed by the spontaneous reaction between urea and ethanol and

53

the accumulation of urea is one of the major causes of overaccumulation of EC in

54

Chinese rice wine.3 The average concentration of EC in Chinese rice wine was

55

reported to be 160 µg/kg.6 However, the maximum level of EC in sake, which is

56

another kind of rice wine similar with Chinese rice wine, is limited to 100 µg/kg in

57

America. Although there has yet no limitation about the maximum level of EC in

58

Chinese rice wine, decreasing the concentration of EC is urgent to improve the quality

59

of Chinese rice wine. The key point of reducing EC in Chinese rice wine is how to

60

reduce the accumulation of urea during fermentation process.

61

Urea is considered a non-preferred nitrogen source for yeast.7 When preferred

62

nitrogen sources such as glutamate and glutamine are present, the utilization of

63

non-preferred nitrogen sources is repressed through a process known as nitrogen

64

catabolite repression (NCR).7-10 In Saccharomyces cerevisiae, NCR is regulated by

65

four global transcription factors (TFs), including the two positive TFs, Gln3 and Gat1,

66

and the two negative TFs, Gzf3 and Dal80.11-12 In the presence of preferred nitrogen

67

sources, Gln3 and Gat1 are phosphorylated by the target of rapamycin complex 1

68

(TORC1) and sequestered in the cytoplasm with the help of the repressor, Ure2.13 3 ACS Paragon Plus Environment

Journal of Agricultural and Food Chemistry

69

However, under conditions of nitrogen starvation, Gln3 and Gat1 can be translocated

70

from the cytoplasm into the nucleus, which facilities their interaction with promoter

71

regions of non-preferred nitrogen utilization genes to activate their expression.10, 14-15

72

In addition to NCR, nitrogen metabolism is also regulated by carbon sources.

73

α-ketoglutarate, an intermediate of carbon metabolism, can be converted into

74

glutamate and glutamine. Thus, it is the major amide group donor in amino acid

75

biosynthesis and links carbon and nitrogen metabolisms.12, 16

76

A series of efforts have been made to reduce NCR to decrease the accumulation

77

of urea in Chinese rice wine,3 including remodeling the regulatory and metabolic

78

pathways of nitrogen metabolism.7, 10, 17-20 However, genetic modification methods

79

such as these are not suitable for the optimization of strains used in the food industry

80

due to public safety concerns related to genetically modified organisms (GMOs).

81

Adaptive laboratory evolution (ALE), which does not require genetic modification,

82

has been developed as a powerful tool to domesticate strains and obtain desired

83

characteristics.21 During ALE, cells are cultured under conditions that promote the

84

accumulation of beneficial mutations that allow the cells to adapt to the selection

85

pressure. ALE has been applied to many microbes,21 including Escherichia coli22-26

86

and S. cerevisiae,27-31 to gain desirable characteristics such as increased growth rates

87

and resistance to environmental stresses. Moreover, with the rapid development of

88

massively parallel DNA sequencing technologies, the cost of resequencing microbial

89

genomes has decreased dramatically, facilitating the identification of causal mutations

90

and the investigation of genotype-phenotype relationships.24-26, 30-31

91

In this study, the industrial Chinese rice wine strain S. cerevisiae XZ-11 was

92

propagated for several hundred generations in parallel serial cultures in which urea

93

was the sole nitrogen source to improve its ability to utilize urea. Genome 4 ACS Paragon Plus Environment

Page 4 of 36

Page 5 of 36

Journal of Agricultural and Food Chemistry

94

resequencing and comparative genomics analysis revealed genetic variants to account

95

for the low urea accumulation in cultures of these evolved strains. The strains

96

screened in this study could be used during the fermentation process to improve the

97

quality of Chinese rice wines. Additionally, the results of genomic analysis could

98

serve as a guide for genetic modification of strains in other non-food industries.

99 100

MATERIALS AND METHODS

101

Strains, media, and culture conditions

102

The S. cerevisiae diploid strain XZ-11 used in this study was an industrial

103

Chinese rice wine strain provided by the Guyuelongshan Shaoxing Wine Company

104

(Shaoxing, China).

105

Media for seed culture (YPD) contained 20 g/L glucose, 20 g/L peptone, and 10

106

g/L yeast extract; 20 g/L agar was added when necessary. Nitrogen-free media

107

contained 20 g/L glucose and 1.74 g/L yeast nitrogen base (YNB) without amino acids

108

or ammonium ion. Media for ALE contained nitrogen-free media supplemented with

109

0.03, 0.3, or 1.5 g/L urea as required. Prescreening media contained nitrogen-free

110

media supplemented with 5 g/L urea. First-round screening media contained

111

nitrogen-free media supplemented with 5 g/L ammonium sulfate, 5 g/L glutamate, 5

112

g/L glutamine, and 5 g/L urea. Second-round screening media contained YPD media

113

supplemented with 0.6 g/L urea. Urea detection media contained nitrogen-free media

114

supplemented with 5 g/L urea. Sporulation media contained 20 g/L potassium acetate

115

supplemented with 20 g/L agar. All strains were cultured at 30°C with shaking at 220

116

rpm. Urea concentration measurements were performed after culture for 54 h. For

117

spore dissection, strains were firstly cultured in the YPD medium at 30°C for 24 h.

118

Then cells were collected by centrifugation at 3000 g for 5 min and washed with 5 ACS Paragon Plus Environment

Journal of Agricultural and Food Chemistry

119

sterile water twice. Next, cells were incubated on the sporulation plate at 30°C for 96

120

h.

121 122

Adaptive laboratory evolution

123

Primary ALE experiments were started from wild-type S. cerevisiae XZ-11

124

frozen stocks after overnight activation on YPD plates. A single colony was then

125

inoculated into a 250 mL flask containing 20 mL YPD medium and cultured for 24 h.

126

Three 20 µL aliquots were transferred into three 14 mL test tubes containing 2 mL

127

nitrogen-free media supplemented with 0.03, 0.3, or 1.5 g/L urea for ALE. A total of

128

20 µL of culture was serially passaged after culture for 2 days (1% of the total culture

129

volume was transferred to the subsequent culture).

130 131

Urea detection by high-performance liquid chromatography

132

Urea was detected using a high-performance liquid chromatography32 system

133

(Agilent 1260, Agilent Biotechnologies, Santa Clara, CA, USA) equipped with a

134

ZORBAX Eclipse XDB C18 column (250×4.6 mm, 5 µm; Agilent Biotechnologies).

135

Twenty microliters samples were injected with the mobile phase at 1 mL/min at 35°C.

136

The mobile phase gradient is shown in Table 1. Urea was detected using a

137

fluorescence detector with excitation wavelength of 213 nm and an emission

138

wavelength of 308 nm.

139 140

The high-throughput urea detection method

141

In order to determine urea concentrations using high-throughput screening, a

142

method from previous studies was optimized.33-34 Briefly, 80 µL of culture

6 ACS Paragon Plus Environment

Page 6 of 36

Page 7 of 36

Journal of Agricultural and Food Chemistry

143

supernatant was collected and mixed with 100 µL iron phosphate solution (6 g/L), 20

144

µL mixture of diacetyl monoxime (6 g/L) and thiamine (0.3 g/L) in 96-well plates.

145

Then, the plates were incubated at 100°C for 10 min. The absorbance at 520 nm of

146

each well was immediately measured using a microplate spectrophotometer (BioTek

147

Instruments, Inc., Winooski, VT, USA).

148 149

Single live cell sorting after evolution using flow cytometry

150

Cultures were diluted with saline to 1×105 cells/mL. An equal volume of

151

propidium iodide (PI, 8 µmol/L) was then added to the samples, followed by

152

incubation at 4°C for 20 min. Cells were gated on a forward and side scatter plot and

153

then classified as either PI-stained or non-stained based on their 617 nm fluorescence

154

and side scatter. Non-stained cells were individually sorted into 96-well plates

155

containing 700 µL prescreening media in each well. Finally, the 96-wells plates were

156

incubated at 30°C for 54 h.

157 158

DNA preparation, genome resequencing, and data analysis

159

Genomic DNA was isolated from the adaptive evolved strains 4B, 7H, and 10G

160

using a PureLink™ Genomic Plant DNA Purification Kit (Thermo Fisher Scientific,

161

Waltham, MA, USA) following the manufacturer’s protocol. Paired-end libraries with

162

300 bp inserts were constructed and sequenced using an Illumina Hiseq 2500

163

sequencer (Illumina, San Diego, CA, USA) at Shanghai Biochip Co., Ltd. (Shanghai,

164

China). The quality of raw data was monitored using FastQC35 and filtered using

165

PRINSEQ.36 The resulting clean data sets were mapped to the genome sequence of

166

the wild-type strain XZ-11 using Bowtie2.37 Genomic variants were analyzed using 7 ACS Paragon Plus Environment

Journal of Agricultural and Food Chemistry

167

samtools38 and VCFtools.39 HMMcopy was used to analyze chromosome copy

168

numbers.40 KaKs_Caculator was applied to measure the gene evolutionary rates of the

169

evolved strains based on genetic variation.41

170 171

RNA preparation, cDNA synthesis, and RT-qPCR assays

172

After overnight incubation on YPD plates, single colonies were cultured to the

173

early log phase in liquid YPD medium for 12 h. Total RNA was extracted from each

174

sample using a RiboPure RNA Purification Kit (Thermo Fisher Scientific) following

175

the manufacturer’s instructions. The cDNA was synthesized using the PrimeScript™

176

RT Reagent Kit with gDNA Eraser (Takara-Bio, Dalian, China). RT-qPCR

177

experiments were performed on a LightCycler 480 system (Roche Applied Science,

178

Mannheim, Germany) using a SYBR® Premix Ex Taq™ Kit (Takara-Bio). Primers

179

used in this assay are listed in Table S1. ACT1 was selected as the housekeeping gene

180

and fold-change was calculated using the 2-∆∆Ct method. Biological triplicates for each

181

sample were analyzed.

182 183

Simulation of Chinese rice wine fermentation

184

Simulations of Chinese rice wine fermentation performed in this study were

185

similar to those described in previous research.42 First, 0.75 kg glutinous rice was

186

soaked for 48 h in 0.9 L water at 28°C. Next, the glutinous rice was washed and

187

steamed for 20 min. Then, 1.5 mL yeast cultured overnight was mixed with the

188

steamed rice and incubated at 28°C. Forty-eight hours later, 75 g Chinese koji

189

(provided by the Guyuelongshan Shaoxing Wine Company) and 0.9 L water were

190

added and incubated at 28°C for 2–3 days to complete the mashing process. The

8 ACS Paragon Plus Environment

Page 8 of 36

Page 9 of 36

Journal of Agricultural and Food Chemistry

191

temperature was then lowered to 18°C for post-fermentation for 16–20 days. The

192

concentration of EC in the stimulated Chinese rice wine was determined by gas

193

chromatography with mass spectrometry (GS/MS) as described by Ryu et al..43

194 195

RESULTS

196

Acquisition of low urea-accumulating strains by adaptive laboratory evolution

197

S. cerevisiae XZ-11 was subjected to continuous cultivation under selective

198

culture conditions in which urea was the sole nitrogen source presented at one of three

199

different concentrations. During the ALE process, a constant passage size (1% of the

200

total culture volume) was used. Log phase cultures were transferred to tubes

201

containing fresh media every 48 h. Cultures were sampled and subjected to batch

202

fermentation every 60 transfers. As mutations accrued and gained dominance within

203

the separate populations, the amount of residual urea decreased (Figure 1). Compared

204

to the initial culture, the percentage of urea residue in cultures containing 0.03, 0.3,

205

and 1.5 g/L urea as sole nitrogen source decreased by 30.1%, 34.0%, and 34.7%,

206

respectively, after 360 rounds of transfer. These results show that stronger selective

207

pressure led to faster adaptation to culture condition.

208 209

High-throughput screening and characteristics of low urea-accumulating strains

210

Using flow cytometry sorting, a total of 4,721 single colonies were isolated from

211

the three end-point evolved populations. Results of high-throughput screening after

212

culture in first-round screening media showed that the amount of residual urea in the

213

cultures of most evolved strains was significantly lower than that in cultures of the

214

wild-type strain (Figure 2a). Preliminary screening revealed nine strains with low urea

215

concentrations that were then subjected to further fermentation. After 54 h of batch 9 ACS Paragon Plus Environment

Journal of Agricultural and Food Chemistry

216

fermentation using YPD medium, the urea concentrations in cultures of four of the

217

nine selected strains were reduced by 5.1%−53.9% (Figure 2b). Further fermentation

218

analysis of the three strains with the lowest urea concentration, named 10G, 7H, and

219

4B, in medium containing 0.6 g/L urea as the sole nitrogen source revealed that their

220

urea consumption capabilities were almost doubled and their growth rates were much

221

faster than the wild-type strain XZ-11 (Figure 2c and d). Compared to the wild-type

222

strain, the residual urea concentration in cultures of 10G, 7H, and 4B were 37.4%,

223

39.9%, and 43.9% lower, respectively.

224

Beneficial mutations gained from ALE were likely located on only one copy of

225

the chromosome and could be lost during budding dissection. Therefore, the three

226

evolved strains were subjected to serial spore dissection. Urea consumption of the

227

homozygotes of each budding generation was found to be stable after six generations

228

(Table 2), suggesting that beneficial mutations may spread evenly across different

229

chromosomes. In addition, in simulated Chinese rice wine fermentation, the final urea

230

concentrations in cultures of the evolved strains 4B, 7H, and 10G were 47.9%, 16.6%,

231

and 12.4% lower than that of the wild-type strain, respectively (Figure 2e). Moreover,

232

the volume percent of ethanol, one of the most important indicators of the quality of

233

Chinese rice wine, was not significantly changed (Figure 2f). The EC concentrations

234

in the cultures of the evolved strains 4B, 7H, and 10G were 40.5%, 28.7%, and 18.6%

235

lower than that of the wild-type strain, respectively (Figure 2g).

236 237

Adaptive laboratory evolution triggers multiple genomic variants

238

Genome resequencing was performed to analyze the genome of strains as they

239

adapted to the urea media. After mapping sequence reads of evolved strains to the

240

genome of wild-type strain XZ-11, many genetic variants were detected (Table 3 and 10 ACS Paragon Plus Environment

Page 10 of 36

Page 11 of 36

Journal of Agricultural and Food Chemistry

241

Table S2). The number of detected genetic variants in the genome of strain 4B (388)

242

was much lower than that in the genomes of strains 7H (19,169) and 10G (19,312).

243

Fifty-five unique variants corresponding to changes in the coding sequences of 19

244

genes were identified in the genome of strain 4B (Figure S1, Table S2). Functional

245

annotation of these 19 genes revealed that most were involved in transportation of

246

hexoses, including fructose, glucose, and mannose. The percentages of non-exonic

247

variants in the genomes of strains 4B, 7H, and 10G were approximately 76.8%, 40.9%,

248

and 41.1%, respectively. Considering that non-coding sequences account for only

249

25% of the S. cerevisiae genome, genomic variants appeared to be enriched in

250

non-coding regions, suggesting they were not randomly distributed.

251

The ratio of the number of non-synonymous substitutions per non-synonymous

252

site (Ka) to the number of synonymous substitutions per synonymous site (Ks), Ka/Ks,

253

is used to estimate the evolutionary rate of genes under selection.44 We identified 15,

254

718, and 831 fast-evolving genes (Ka/Ks value > 1) in strains 4B, 7H, and 10G,

255

respectively. Some of these genes are involved in the amino acid metabolism pathway,

256

including AAT2, which is involved in converting α-ketoglutarate to glutamate,45 and

257

PRO1, which is involved in converting glutamate to ornithine and entering into the

258

urea cycle. ARG2, another gene involved in the generation of ornithine from

259

glutamate,46 was identified as a fast-evolving gene in strains 4B and 7H but not in

260

strain 10G. The duplicated chromosomal regions in the genomes of the evolved strains

261

are shown in Figure 3. In general, these were evenly distributed throughout the

262

genomes of the three evolved strains. Remarkably, a large region (approximately 38.5

263

kb) in chromosome XII was found to undergo massive amplification with a final copy

264

number of approximately 10 (Figure 3). Twenty-two annotated genes are located in

265

this region (Table S3), including IDP2, which encodes the isocitrate dehydrogenase 11 ACS Paragon Plus Environment

Journal of Agricultural and Food Chemistry

266

responsible for converting oxaloacetate into α-ketoglutarate.

267 268

Adaptive laboratory evolution leads to upregulation of genes related to urea

269

utilization

270

The transcriptional level of genes related to urea utilization, including genes

271

involved in urea metabolism and regulatory pathways (Figure 4), were examined after

272

culture in YPD media. In the three evolved strains, these genes were significantly

273

upregulated compared to the wild-type strain (Figure 5). Genes directly involved in

274

the metabolism of urea, such as DUR1,2, DUR3, and CAR1, were upregulated by

275

81.2-, 8.1-, and 24.8-fold, respectively, in strain 4B compared to the wild-type strain.

276

Furthermore, genes encoding the positive regulators such as GLN3 and GAT1 were

277

upregulated to a greater extent than those encoding the negative regulators such as

278

GZF3, DAL80, and URE2. These results may suggest that beneficial mutations in

279

GLN3 and GAT1 coding regions contributed to the derepression of NCR. In addition,

280

an obvious correlation between the expression levels of urea metabolism genes and

281

the urea consumption ability of strains was observed (Figure 5). Generally,

282

upregulation of genes was highest in strain 4B and lowest in strain 10G (Figure 5).

283 284

DISCUSSION

285

In this study, a novel strategy integrating ALE and high-throughput screening

286

was established to identify Chinese rice yeast strains that accumulate low urea

287

concentrations during fermentation. Three adaptive evolved strains, named 4B, 7H,

288

and 10G, were selected using this strategy. Compared to the wild-type strain XZ-11,

289

the residual concentrations of urea in cultures of 4B, 7H, and 10G were reduced by

290

47.9%, 16.6%, and 12.4%, respectively. Meanwhile, the fermentation performances of 12 ACS Paragon Plus Environment

Page 12 of 36

Page 13 of 36

Journal of Agricultural and Food Chemistry

291

these strains remained similar to that of strain XZ-11, while the EC concentrations in

292

cultures of evolved strains 4B, 7H, and 10G were reduced by 40.5%, 28.7%, and

293

18.6% than that in culture of wild-type strain. NCR regulators and genes involved in

294

urea catabolism were significantly upregulated in low urea-accumulating strains.

295

Genome resequencing and comparative genomic analysis revealed that many variants

296

accumulated in the genomes of evolved strains during ALE. Genes involved in linking

297

carbon metabolism and the urea cycle pathway, including ATT1, PRO1, and ARG2,

298

evolved more quickly.

299

Urea is considered to be the major precursor of EC in Chinese rice wine.32

300

Previously, the reduction of urea in alcoholic beverages was achieved mainly by

301

genetically modifying fermentation strains. For example, promoting urea

302

degradation by enhancing the expression of DUR1,2 and blocking urea formation by

303

deleting CAR1 decreased urea accumulation by 75.6% and 86.9%, respectively.17, 19

304

Similarly, stabilizing the urea permease Dur3 by mutating its potential ubiquitination

305

sites decreased urea accumulation by 29.7%.47 Furthermore, eliminating the

306

repression of urea utilization by mutating potential phosphorylation sites and

307

truncating the localization regulation signal of the NCR regulator Gln3 reduced urea

308

concentrations by 63%.10 Although these methods greatly reduce urea concentrations,

309

the potential safety risks associated with GMOs limit their application in industrial

310

processes. In contrast, strains screened after ALE exhibit much lower urea

311

accumulation than the wild-type strain, are free from genetic modification, and are

312

safe for use in fermented beverage production. Moreover, adaptive evolved strains

313

exhibit similar fermentation properties to the wild-type strain in simulations of

314

Chinese rice wine fermentation.

315

In the past decades, ALE has been applied to many microbes to obtain desired 13 ACS Paragon Plus Environment

Journal of Agricultural and Food Chemistry

316

characteristics without genetic modification.21 During ALE, cells are cultured under

317

significant selection pressure, which promotes the enrichment of beneficial

318

mutations and allows for the cells to adapt to the selection pressure. However, a

319

simple high-throughput screening method is essential to isolate adaptive evolved

320

strains with desired characteristics. Based on our previous work, we implemented a

321

flow cytometry-based sorting method to isolate single colonies from large evolved

322

populations.34 This is a much simpler and more convenient method than the

323

traditional dilution separation method, which requires single colonies to be picked up,

324

which is laborious. Moreover, urea concentrations can be rapidly and conveniently

325

determined by measuring the absorbance of each well at 520 nm. By combining flow

326

cytometry sorting with a convenient urea detection method, nine evolved strains with

327

low urea accumulation levels were easily selected from 4,721 candidate colonies in

328

this study.

329

The utilization of urea is controlled by the systematic regulation of nitrogen

330

metabolism. Many regulatory pathways have been found to be involved in nitrogen

331

metabolism, including the SPS-sensor system,48 the TOR pathway,49 the general

332

amino acid pathway,50-51 and NCR.7, 10, 52 In this study, the transcription levels of urea

333

cycle genes and regulators of the NCR and TOR pathways, including CAR1, were

334

found to be upregulated, suggesting that, overall, transcriptional activity was activated

335

after ALE. However, greater activation of DUR1,2 compared to CAR1 may lead to

336

overall reductions in urea concentrations after ALE. In addition, the upregulation of

337

IDP2 and the fast-evolving genes ATT2, PRO1, and ARG2 were confirmed by

338

comparative genomics analysis. These genes are involved in the metabolic pathway

339

by which α-ketoglutarate is converted to glutamate and finally directed to ornithine.

340

Importantly,

α-ketoglutarate

is essential for linking nitrogen and carbon 14

ACS Paragon Plus Environment

Page 14 of 36

Page 15 of 36

Journal of Agricultural and Food Chemistry

341

metabolism.53-54 These results suggested that the remodeling of carbon metabolism

342

plays an important role in reducing urea accumulation during ALE.

343

Three Chinese rice wine strains with low urea accumulation properties were

344

selected using ALE, a genetic modification-free approach. However, the level of urea

345

reduction was still relatively low compared to that achieved by genetic modification

346

or urease addition.17-19,

347

concentrations could be performed. With the help of the high-throughput screening

348

method established here, large mutation libraries can be constructed using other

349

efficient approaches in addition to ALE, such as random mutagenesis. Identification

350

of carbon catabolism-related genes by comparative genomic analysis of the adaptive

351

evolved and wild-type strains may provide new approaches for reducing urea

352

accumulation and could provide new clues for coordinating the metabolic flux of

353

carbon and nitrogen sources.

55

Further evolution in the presence of higher urea

354

15 ACS Paragon Plus Environment

Journal of Agricultural and Food Chemistry

Page 16 of 36

355 356

ASSOCIATED CONTENT

357

Supporting Information

358

Table S1, primers used in this study; Table S2, genetic variants in the genome of

359

adaptive evolution improved strains 4B, 7H, and 10G; Table S3, functional annotation

360

of the genes in the massive amplified region in the chromosome XII of evolved strains;

361

and Figure S1, comparative analysis of the variants among three evolved strains.

362 363 364

AUTHOR INFORMATION

365

Corresponding Authors

366

*Phone:

367

[email protected] (Jingwen Zhou), [email protected] (Jian Chen).

+86-510-85914317,

Fax:

+86-510-85914317.

E-mail:

368 369

Funding

370

This work was supported by The National Key Research and Development Program

371

of China (2017YFC1600403), the National Natural Science Foundation of China

372

(31670095, 31770097, 31671836), the Key Research and Development Program of

373

Jiangsu Province (BE2016689), the Fundamental Research Funds for the Central

374

Universities (JUSRP51701A), the Six Talent Peaks Project in Jiangsu Province

375

(2015-JY-005), the Distinguished Professor Project of Jiangsu Province.

376 377

Notes

378

The authors declare no conflict of competing financial interests.

16 ACS Paragon Plus Environment

Page 17 of 36

Journal of Agricultural and Food Chemistry

379

REFERNCES

380

1. McGovern, P. E.; Zhang, J.; Tang, J.; Zhang, Z.; Hall, G. R.; Moreau, R. A.;

381

Nunez, A.; Butrym, E. D.; Richards, M. P.; Wang, C. S.; Cheng, G.; Zhao, Z.; Wang,

382

C., Fermented beverages of pre- and proto-historic China. Proc. Natl. Acad. Sci. U. S.

383

A. 2004, 101 (51), 17593-17598.

384 385

2. Guangfa, X., Functional constituents and health function of Chinese rice wine. Liquor Making 2008, 35, 14-16.

386

3. Zhao, X.; Du, G.; Zou, H.; Fu, J.; Zhou, J.; Chen, J., Progress in preventing

387

the accumulation of ethyl carbamate in alcoholic beverages. Trends Food Sci. Tech.

388

2013, 32 (2), 97-107.

389 390 391

4. IARC, Alcohol Consumption and Ethyl Carbamate. Volume 96, IARC. 2011; Vol. 31, p 2402. 5. Baan, R.; Straif, K.; Grosse, Y.; Secretan, B.; El Ghissassi, F.; Bouvard, V.;

392

Altieri, A.; Cogliano, V., Carcinogenicity of alcoholic beverages. Lancet Oncol. 2007,

393

8 (4), 292-3.

394

6. Wu, P.; Pan, X.; Wang, L.; Shen, X.; Yang, D., A survey of ethyl carbamate in

395

fermented foods and beverages from Zhejiang, China. Food Control 2012, 23 (1),

396

286-288.

397

7. Zhao, X.; Zou, H.; Chen, J.; Du, G.; Zhou, J., The modification of Gat1p in

398

nitrogen catabolite repression to enhance non-preferred nitrogen utilization in

399

Saccharomyces cerevisiae. Sci. Rep. 2016, 6, 21603.

400

8. Zhao, S.; Zhao, X.; Zou, H.; Fu, J.; Du, G.; Zhou, J.; Chen, J., Comparative

401

proteomic analysis of Saccharomyces cerevisiae under different nitrogen sources. J.

402

Proteomics 2014, 101, 102-112.

403

9. Zhao, X.; Zou, H.; Du, G.; Chen, J.; Zhou, J., Effects of nitrogen catabolite 17 ACS Paragon Plus Environment

Journal of Agricultural and Food Chemistry

404

repression-related amino acids on the flavour of rice wine. J. I. Brewing 2015, 121 (4),

405

581-588.

406

10. Zhao, X.; Zou, H.; Fu, J.; Zhou, J.; Du, G.; Chen, J., Metabolic engineering

407

of the regulators in nitrogen catabolite repression to reduce the production of ethyl

408

carbamate in a model rice wine system. App. Environ. Microb. 2014, 80 (1), 392-398.

409

11. Cooper, T. G., Transmitting the signal of excess nitrogen in Saccharomyces

410

cerevisiae from the Tor proteins to the GATA factors: connecting the dots. FEMS.

411

Microbiol. Rev. 2002, 26 (3), 223-238.

412 413

12. Magasanik, B.; Kaiser, C. A., Nitrogen regulation in Saccharomyces cerevisiae. Gene 2002, 290 (1-2), 1-18.

414

13. Bertram, P. G.; Choi, J. H.; Carvalho, J.; Ai, W.; Zeng, C.; Chan, T. F.; Zheng,

415

X. F., Tripartite regulation of Gln3p by TOR, Ure2p, and phosphatases. J. Biol. Chem.

416

2000, 275 (46), 35727-35733.

417

14. Cox, K. H.; Rai, R.; Distler, M.; Daugherty, J. R.; Coffman, J. A.; Cooper, T.

418

G., Saccharomyces cerevisiae GATA sequences function as TATA elements during

419

nitrogen catabolite repression and when Gln3p is excluded from the nucleus by

420

overproduction of Ure2p. J. Biol. Chem. 2000, 275 (23), 17611-17618.

421

15. Kulkarni, A.; Buford, T. D.; Rai, R.; Cooper, T. G., Differing responses of

422

Gat1 and Gln3 phosphorylation and localization to rapamycin and methionine

423

sulfoximine treatment in Saccharomyces cerevisiae. FEMS. Yeast. Res. 2006, 6 (2),

424

218-229.

425

16. Chisholm, G.; Cooper, T. G., Isolation and characterization of mutants that

426

produce the allantoin-degrading enzymes constitutively in Saccharomyces cerevisiae.

427

Mol. Cell. Biol. 1982, 2 (9), 1088-1095.

428

17. Wu, D.; Li, X.; Lu, J.; Chen, J.; Xie, G.; Zhang, L., The overexpression of 18 ACS Paragon Plus Environment

Page 18 of 36

Page 19 of 36

Journal of Agricultural and Food Chemistry

429

DUR1,2 and deletion of CAR1 in an industrial Saccharomyces cerevisiae strain and

430

effects on nitrogen catabolite repression in Chinese rice wine production. J. I.

431

Brewing 2016, 122 (3), 480-485.

432

18. Wu, D.; Li, X.; Lu, J.; Chen, J.; Zhang, L.; Xie, G., Constitutive expression of

433

the DUR1,2 gene in an industrial yeast strain to minimize ethyl carbamate production

434

during Chinese rice wine fermentation. FEMS Microbiol. Lett. 2016, 363 (1), fnv214.

435

19. Wu, D.; Li, X.; Shen, C.; Lu, J.; Chen, J.; Xie, G., Decreased ethyl carbamate

436

generation during Chinese rice wine fermentation by disruption of CAR1 in an

437

industrial yeast strain. Int. J. Food Microbiol. 2014, 180, 19-23.

438

20. Zhang, P.; Du, G.; Zou, H.; Chen, J.; Xie, G.; Shi, Z.; Zhou, J., Effects of

439

three permeases on arginine utilization in Saccharomyces cerevisiae. Sci. Rep. 2016, 6,

440

20910.

441 442

21. Dragosits, M.; Mattanovich, D., Adaptive laboratory evolution-principles and applications for biotechnology. Microb. Cell. Fact. 2013, 12, 64.

443

22. Fu, W.; Guomundsson, O.; Paglia, G.; Herjolfsson, G.; Andresson, O. S.;

444

Palsson, B. O.; Brynjolfsson, S., Enhancement of carotenoid biosynthesis in the green

445

microalga Dunaliella salina with light-emitting diodes and adaptive laboratory

446

evolution. Appl. Microbiol. Biot. 2013, 97 (6), 2395-2403.

447

23. LaCroix, R. A.; Palsson, B. O.; Feist, A. M., A model for designing adaptive

448

laboratory evolution experiments. Appl. Environ. Microbiol. 2017, 83 (8),

449

e03115-03116.

450

24. LaCroix, R. A.; Sandberg, T. E.; O'Brien, E. J.; Utrilla, J.; Ebrahim, A.;

451

Guzman, G. I.; Szubin, R.; Palsson, B. O.; Feist, A. M., Use of adaptive laboratory

452

evolution to discover key mutations enabling rapid growth of Escherichia coli K-12

453

MG1655 on glucose minimal medium. Appl. Environ. Microbiol. 2015, 81 (1), 17-30. 19 ACS Paragon Plus Environment

Journal of Agricultural and Food Chemistry

454

25. Luan, G.; Cai, Z.; Li, Y.; Ma, Y., Genome replication engineering assisted

455

continuous evolution (GREACE) to improve microbial tolerance for biofuels

456

production. Biotechnol. Biofuels 2013, 6 (1), 137.

457

26. Sandberg, T. E.; Pedersen, M.; LaCroix, R. A.; Ebrahim, A.; Bonde, M.;

458

Herrgard, M. J.; Palsson, B. O.; Sommer, M.; Feist, A. M., Evolution of Escherichia

459

coli to 42 degrees C and subsequent genetic engineering reveals adaptive mechanisms

460

and novel mutations. Mol. Biol. Evol. 2014, 31 (10), 2647-2662.

461

27.Gonzalez-Ramos, D.; Gorter de Vries, A. R.; Grijseels, S. S.; van Berkum, M.

462

C.; Swinnen, S.; van den Broek, M.; Nevoigt, E.; Daran, J. M.; Pronk, J. T.; van Maris,

463

A. J., A new laboratory evolution approach to select for constitutive acetic acid

464

tolerance in Saccharomyces cerevisiae and identification of causal mutations.

465

Biotechnol. Biofuels 2016, 9, 173.

466

28. Peris, D.; Moriarty, R. V.; Alexander, W. G.; Baker, E.; Sylvester, K.; Sardi,

467

M.; Langdon, Q. K.; Libkind, D.; Wang, Q. M.; Bai, F. Y.; Leducq, J. B.; Charron, G.;

468

Landry, C. R.; Sampaio, J. P.; Goncalves, P.; Hyma, K. E.; Fay, J. C.; Sato, T. K.;

469

Hittinger, C. T., Hybridization and adaptive evolution of diverse Saccharomyces

470

species for cellulosic biofuel production. Biotechnol. Biofuels 2017, 10, 78.

471

29. Sato, T. K.; Tremaine, M.; Parreiras, L. S.; Hebert, A. S.; Myers, K. S.;

472

Higbee, A. J.; Sardi, M.; McIlwain, S. J.; Ong, I. M.; Breuer, R. J.; Avanasi

473

Narasimhan, R.; McGee, M. A.; Dickinson, Q.; La Reau, A.; Xie, D.; Tian, M.; Reed,

474

J. L.; Zhang, Y.; Coon, J. J.; Hittinger, C. T.; Gasch, A. P.; Landick, R., Directed

475

evolution reveals unexpected epistatic interactions that alter metabolic regulation and

476

enable anaerobic xylose use by Saccharomyces cerevisiae. PLoS. Genet. 2016, 12

477

(10), e1006372.

478

30. Satomura, A.; Miura, N.; Kuroda, K.; Ueda, M., Molecular analysis of 20 ACS Paragon Plus Environment

Page 20 of 36

Page 21 of 36

Journal of Agricultural and Food Chemistry

479

adaptational strategies under heat stress by whole-genome sequencing of

480

laboratory-evolved Saccharomyces cerevisiae. FASEB J. 2016, 30 (1 Supplement),

481

585.2-585.2.

482

31. Zeng, W. Y.; Tang, Y. Q.; Gou, M.; Sun, Z. Y.; Xia, Z. Y.; Kida, K.,

483

Comparative transcriptomes reveal novel evolutionary strategies adopted by

484

Saccharomyces cerevisiae with improved xylose utilization capability. Appl.

485

Microbiol. Biotechnol. 2017, 101 (4), 1753-1767.

486

32. Xing, J.; Zhong, Q.; Xiong, Z.; Zou, H.; Zhou, J.; Xie, G., Determination of

487

urea content in yellow rice wine by high-performance liquid chromatography

488

combined with fluorescence detection after pre-column derivatisation. Liquor-Making

489

Sci. Tech. 2011, 104-107.

490

33. Alam, M. S.; Casareto, B. E.; Suzuki, Y.; Sultana, R.; Suzuki, T.,

491

Optimization of dissolved urea measurements in coastal waters with the combination

492

of a single reagent and high temperature. J. Oceanogr. 2017, 73 (2), 249-258.

493

34. Cheng, Y.; Du, G.; Zhou, J.; Chen, J., High-throughput screening of mutant

494

strain to reduce the accumulation of ethyl carbamate precursor in rice wine

495

fermentation. Acta. Microbiol. Sin. 2017, 57 (10), 1517-1526.

496 497 498 499 500 501 502 503

35. Andrews, S., FastQC: a quality control tool for high throughput sequence data. 2010. 36. Schmieder, R.; Edwards, R., Quality control and preprocessing of metagenomic datasets. Bioinformatics 2011, 27 (6), 863-864. 37. Langmead, B.; Salzberg, S. L., Fast gapped-read alignment with Bowtie 2. Nat. Methods. 2012, 9 (4), 357-359. 38. Li, H.; Handsaker, B.; Wysoker, A.; Fennell, T.; Ruan, J.; Homer, N.; Marth, G.; Abecasis, G.; Durbin, R.; Genome Project Data Processing, S., The Sequence 21 ACS Paragon Plus Environment

Journal of Agricultural and Food Chemistry

504

Alignment/Map format and SAMtools. Bioinformatics 2009, 25 (16), 2078-9.

505

39. Danecek, P.; Auton, A.; Abecasis, G.; Albers, C. A.; Banks, E.; DePristo, M.

506

A.; Handsaker, R. E.; Lunter, G.; Marth, G. T.; Sherry, S. T.; McVean, G.; Durbin, R.;

507

Genomes Project Analysis, G., The variant call format and VCFtools. Bioinformatics

508

2011, 27 (15), 2156-8.

509 510

40. Lai, D.; Ha, G.; Shah, S., HMMcopy: copy number prediction with correction for GC and mappability bias for HTS data. R package version 2016, 1 (0).

511

41. Wang, D.; Zhang, Y.; Zhang, Z.; Zhu, J.; Yu, J., KaKs_Calculator 2.0: a

512

toolkit incorporating gamma-series methods and sliding window strategies. Genomics

513

Proteomics Bioinformatics 2010, 8 (1), 77-80.

514

42. Zhou, W.; Fang, R.; Chen, Q., Effect of gallic and protocatechuic acids on the

515

metabolism of ethyl carbamate in Chinese yellow rice wine brewing. Food Chem.

516

2017, 233, 174-181.

517

43. Ryu, D.; Choi, B.; Kim, E.; Park, S.; Paeng, H.; Kim, C. I.; Lee, J. Y.; Yoon,

518

H. J.; Koh, E., Determination of ethyl carbamate in alcoholic beverages and

519

fermented foods sold in Korea. Toxicol. Res. 2015, 31 (3), 289-297.

520 521 522

44. Hurst, L. D., The Ka/Ks ratio: diagnosing the form of sequence evolution. Trends. Genet. 2002, 18 (9), 486. 45.Verleur, N.; Elgersma, Y.; Van Roermund, C. W.; Tabak, H. F.; Wanders, R. J.,

523

Cytosolic aspartate aminotransferase encoded by the AAT2 gene is targeted to the

524

peroxisomes in oleate-grown Saccharomyces cerevisiae. Eur. J. Biochem. 1997, 247

525

(3), 972-980.

526

46. Abadjieva, A.; Pauwels, K.; Hilven, P.; Crabeel, M., A new yeast metabolon

527

involving at least the two first enzymes of arginine biosynthesis: acetylglutamate

528

synthase activity requires complex formation with acetylglutamate kinase. J. Biol. 22 ACS Paragon Plus Environment

Page 22 of 36

Page 23 of 36

Journal of Agricultural and Food Chemistry

529 530

Chem. 2001, 276 (46), 42869-42880. 47. Zhang, P.; Du, G.; Zou, H.; Xie, G.; Chen, J.; Shi, Z.; Zhou, J., Mutant

531

potential ubiquitination sites in Dur3p enhance the urea and ethyl carbamate reduction

532

in a model rice wine system. J.Agr. Food Chem. 2017, 65 (8), 1641-1648.

533

48. Forsberg, H.; Ljungdahl, P. O., Genetic and biochemical analysis of the yeast

534

plasma membrane Ssy1p-Ptr3p-Ssy5p sensor of extracellular amino acids. Mol. Cell.

535

Biol. 2001, 21 (3), 814-826.

536

49. Conrad, M.; Schothorst, J.; Kankipati, H. N.; Van Zeebroeck, G.;

537

Rubio-Texeira, M.; Thevelein, J. M., Nutrient sensing and signaling in the yeast

538

Saccharomyces cerevisiae. FEMS. Microbiol. Rev. 2014, 38 (2), 254-299.

539 540 541

50. Hinnebusch, A. G., Translational regulation of GCN4 and the general amino acid control of yeast. Annu. Rev. Microbiol. 2005, 59, 407-450. 51. Hinnebusch, A. G.; Natarajan, K., Gcn4p, a master regulator of gene

542

expression, is controlled at multiple levels by diverse signals of starvation and stress.

543

Eukaryot. Cell 2002, 1 (1), 22-32.

544 545

52. Hofman-Bang, J., Nitrogen catabolite repression in Saccharomyces cerevisiae. Mol. Biotechnol. 1999, 12 (1), 35-73.

546

53. Cogoni, C.; Valenzuela, L.; Gonzalez-Halphen, D.; Olivera, H.; Macino, G.;

547

Ballario, P.; Gonzalez, A., Saccharomyces cerevisiae has a single glutamate synthase

548

gene coding for a plant-like high-molecular-weight polypeptide. J. Bacteriol. 1995,

549

177 (3), 792-798.

550 551 552 553

54. Liu, Z.; Butow, R. A., Mitochondrial retrograde signaling. Annu. Rev. Genet. 2006, 40, 159-185. 55. Yang, L. Q.; Wang, S. H.; Tian, Y. P., Purification, properties, and application of a novel acid urease from Enterobacter sp. Appl. Biochem. Biotechnol. 2010, 160 23 ACS Paragon Plus Environment

Journal of Agricultural and Food Chemistry

554

(2), 303-313.

555

24 ACS Paragon Plus Environment

Page 24 of 36

Page 25 of 36

Journal of Agricultural and Food Chemistry

556

FIGURE CAPTIONS

557

Figure 1. Urea residue percentages of evolved populations. The percentage of urea

558

in cultures of evolved populations tested every 60 transfers compared to the original

559

strain. Square, circle, and triangle symbols indicate populations that evolved with 0.03,

560

0.3, and 1.5 g/L urea as the sole nitrogen source, respectively. Error bars represent

561

standard deviations from three biological replicates.

562 563

Figure 2. Fermentation profiles of evolved strains. a, The percentage of urea

564

relative to the wild-type (WT) strain XZ-11 of 4,721 single colonies cultured for 54 h

565

after adaptive laboratory evolution in nitrogen-free media containing 5 g/L each

566

glutamine, glutamate, (NH4)2SO4, and urea. b, the concentration of urea in the media

567

of nine evolved strains and the WT strain XZ-11 after culture for 54 h in YPD

568

containing 5 g/L each glutamine, glutamate, (NH4)2SO4, and urea. c, d, The growth

569

and urea concentration curves of the three selected evolved strains and the WT strain

570

after culture in nitrogen-free media containing 5 g/L urea. e, f, g, Final concentrations

571

of urea, ethanol and EC of the three selected evolved strains and the WT strain after

572

Chinese rice wine fermentation simulations. Error bars represent standard deviations

573

from three biological replicates.

574 575

25 ACS Paragon Plus Environment

Journal of Agricultural and Food Chemistry

576

Figure 3. Chromosomal copy number variants in the genomes of evolved strains.

577

Three adaptive evolved strains, 4B, 7H, and 10G, were selected after adaptive

578

laboratory evolution. Chromosome copy numbers were analyzed based on sequencing

579

data using HMMcopy. Individual chromosomes, indicated by Roman numerals, are

580

separated by dashed lines.

581 582

Figure 4. The metabolic and regulatory pathways of urea. Urea is mainly

583

generated from the catabolism of arginine in the urea cycle. The product of DUR1,2

584

mediates the complete degradation of urea to CO2. DUR3 encodes the permease

585

responsible for urea import and export. The activation of these urea metabolism genes

586

is dependent on the positive regulators of nitrogen catabolite repression (NCR), Gln3

587

and Gat1. However, when nitrogen sources are replete, Gln3 and Gat1 are sequestered

588

in the cytoplasm, which represses expression of genes involved in urea metabolism.

589

The TOR pathway triggers phosphorylation of Gln3 and Gat1, blocking their

590

translocation from the cytoplasm into the nucleus.

591 592 593

Figure 5. Transcriptional analysis of genes involved in urea metabolism and

594

regulatory pathways. The three evolved strains, 4B, 7H, and 10G, were cultured in

595

YPD media supplemented with 5 g/L urea for 12 h. The expression levels of the genes

596

involved in urea metabolism and regulatory pathways were measured by RT-qPCR

597

using the wild-type strain XZ-11 as the control. Data were normalized to the

598

expression level of the ACT1 gene. Error bars represent standard deviations of three

599

biological replicates. Black, dark gray, and light gray bars indicate the 7H, 4B, and

600

10G strains, respectively. 26 ACS Paragon Plus Environment

Page 26 of 36

Page 27 of 36

Journal of Agricultural and Food Chemistry

TABLES Table 1. Mobile phase gradients for urea detection

a

Time (min)

A%a

B%b

C%c

0.00

80

20

0

0.06

80

20

0

12.06

50

50

0

12.80

0

50

50

15.60

0

100

0

21.60

0

100

0

23.80

0

20

80

24.50

80

20

0

30.00

80

20

0

31.00

80

20

0

20 mM sodium acetate; bacetonitrile; cH2O.

27 ACS Paragon Plus Environment

Journal of Agricultural and Food Chemistry

Page 28 of 36

Table 2. Genetic stability of adaptive evolved strains

Strains

Residual urea (mg/L) 1st

2nd

3rd

4th

5th

6th

Generation

Generation

Generation

Generation

Generation

Generation

WT

15.48 ± 0.16

15.2 ± 0.12

15.6 ± 0.11

15.3 ± 0.15

15.5 ± 0.08

15.4 ± 0.08

4B

9.5 ± 0.15

9.6 ± 0.15

9.5 ± 0.17

9.6 ± 0.24

9.7 ± 0.16

9.6 ± 0.19

7H

10.7 ± 0.06

10.6 ± 0.10

10.6 ± 0.11

10.5 ± 0.07

10.5 ± 0.10

10.6 ± 0.07

10G

10.4 ± 0.07

10.6 ± 0.10

10.7 ± 0.15

10.7 ± 0.16

10.6 ± 0.10

10.8 ± 0.09

28 ACS Paragon Plus Environment

Page 29 of 36

Journal of Agricultural and Food Chemistry

Table 3. Genetic variations identified in adaptive evolved strains

4B

7H

10G

Total

Hom a

Het b

Total

Hom

Het

Total

Hom

Het

90

39

51

11324

10628

696

11373

10683

690

Deletion

3

3

0

110

98

12

115

107

8

Insertion

20

18

2

170

138

32

175

143

32

MNV c

0

0

0

198

131

67

196

138

58

Replacement

1

1

0

9

6

3

11

8

3

SNV d

66

17

49

10837

10255

582

10876

10287

589

298

195

103

7845

7809

36

7939

7912

27

Exonic

Non-Exonic a

Hom, homozygous; bHet, heterozygous; cMNV, multi-nucleotide variant; dSNV,

single-nucleotide variant.

29 ACS Paragon Plus Environment

Journal of Agricultural and Food Chemistry

FIGURE GRAPHICS Figure 1

30 ACS Paragon Plus Environment

Page 30 of 36

Page 31 of 36

Journal of Agricultural and Food Chemistry

Figure 2

31 ACS Paragon Plus Environment

Journal of Agricultural and Food Chemistry

32 ACS Paragon Plus Environment

Page 32 of 36

Page 33 of 36

Journal of Agricultural and Food Chemistry

Figure 3

33 ACS Paragon Plus Environment

Journal of Agricultural and Food Chemistry

Figure 4

34 ACS Paragon Plus Environment

Page 34 of 36

Page 35 of 36

Journal of Agricultural and Food Chemistry

Figure 5

35 ACS Paragon Plus Environment

Journal of Agricultural and Food Chemistry

The TOC graphic

36 ACS Paragon Plus Environment

Page 36 of 36