Knockout of rapC improves bacillomycin D yield ... - ACS Publications

(SOAP) denovo_v2.04 using the clean data 19, 20. Before assembling, k-mer analysis. 96 was used to estimate the size of genome, the degree of heterozy...
0 downloads 0 Views 1MB Size
Subscriber access provided by UNIV OF DURHAM

Biotechnology and Biological Transformations

Knockout of rapC improves bacillomycin D yield based on de novo genome sequencing of Bacillus amyloliquefaciens fmbJ Jing Sun, Shiquan Qian, Jing Lu, Yanan Liu, Fengxia Lu, Xiaomei Bie, and Zhaoxin Lu J. Agric. Food Chem., Just Accepted Manuscript • DOI: 10.1021/acs.jafc.8b00418 • Publication Date (Web): 12 Apr 2018 Downloaded from http://pubs.acs.org on April 12, 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 32

Journal of Agricultural and Food Chemistry

Knockout of rapC improves bacillomycin D yield based on de novo genome sequencing of Bacillus amyloliquefaciens fmbJ Jing Sun, Shiquan Qian, Jing Lu, Yanan Liu, Fengxia Lu, Xiaomei Bie, Zhaoxin Lu*

College of Food Science and Technology, Nanjing Agricultural University, 1 Weigang, Nanjing, China, 210095

*

Corresponding author, Tel.: +86-25-84396583; Fax: +86-25-84396583. E-mail address: [email protected] (Zhaoxin Lu) Present address: Weigang 1, Nanjing, Jiangsu Province 210095, P. R. China

ACS Paragon Plus Environment

Journal of Agricultural and Food Chemistry

1

Abstract

2

Bacillus amyloliquefaciens, a gram-positive and soil-dwelling bacterium, could

3

produce secondary metabolites that suppress plant pathogens. In this study, we

4

provided the whole genome sequence results of B. amyloliquefaciens fmbJ which had

5

one circular chromosome of 4,193,344 bp with 4,249 genes, 87 transfer RNA genes,

6

and 27 rRNA genes. In addition, fmbJ was found to contain several gene clusters of

7

antimicrobial lipopeptides (bacillomycin D, surfactin, and fengycin) and bacillomycin

8

D homologues were further comprehensively identified. To clarify the influence of

9

rapC regulating the synthesis of lipopeptide on the yield of bacillomycin D, rapC

10

gene in fmbJ was successfully deleted by marker-free method. Finally, it was found

11

that the deletion of rapC gene in fmbJ significantly improved bacillomycin D

12

production from 240.7 ± 18.9 mg/L to 360.8 ± 30.7 mg/L, attributed to the increased

13

the expression of bacillomycin D synthesis-related genes through enhancing the

14

transcriptional level of comA, comP, and phrC. These results showed that the

15

production of bacillomycin D in B. amyloliquefaciens fmbJ might be regulated by the

16

RapC-PhrC system. The findings are expected to advance further agricultural

17

application of bacillus spp. as a promising source of natural bioactive compounds.

18

Keywords:Bacillus amyloliquefaciens, genome sequence, bacillomycin D, rapC,

19

knockout

20

ACS Paragon Plus Environment

Page 2 of 32

Page 3 of 32

Journal of Agricultural and Food Chemistry

21

1. Introduction

22

Many plant diseases caused by pathogenic microorganisms lead to decrease in

23

quality and yield of the important crops. With increasing concern about environmental

24

sustainability and the health of consumers, the application of synthetically chemical

25

fungicides in agriculture is restricted 1. As the plant growth-promoting rhizobacteria

26

(PGPR), Bacillus spp. have shown good prospects for the replacement of chemical

27

fungicides in sustainable agriculture 2. They were well known to secrete compounds

28

promoting plant growth, and produce many secondary metabolites with antimicrobial

29

activity

30

bio-control strain, showed strong antagonistic effect on some fungal pathogens, such

31

as Fusarium graminearum 5, Aspergillus ochraceus 6, A. flavus 7, and Rhizopus

32

stolonifer 8. Its inhibitory effect on these pathogens is mainly due to the produced

33

secondary metabolites. However, their practical application were restricted by poor

34

yield of these secondary metabolites. To improve their yield in B. amyloliquefaciens

35

fmbJ, the sequencing of whole genome will provide a deeper understanding on the

36

biosynthesis of secondary metabolites, and important information for the

37

biotechnological modification of the strain.

3, 4

. Bacillus amyloliquefaciens fmbJ (formerly Bacillus subtilis fmbJ), a

38

B. amyloliquefaciens can be used as a biological control agent to destroy its

39

rivals (such as fungi and bacteria). They not only need complex regulatory networks

40

to respond to the extracellular stimuli variation and control the community

41

differentiation9, but also can secrete antimicrobial active substances (such as

42

lipopeptides). There is the inseparable relationship between many two-component

ACS Paragon Plus Environment

Journal of Agricultural and Food Chemistry

43

systems and other systems and the synthesis of antimicrobial lipopeptides in Bacillus

44

strains. The system of Rap-Phr, as one of the core components of the complicated

45

regulatory network, is composed by response regulator aspartate phosphatase (Rap)

46

and its inhibitory oligopeptide (Phr) 10. Several Rap proteins can regulate the process

47

of community differentiation through reducing phosphorylation modification or

48

binding of the response regulator proteins 11. The literature showed that RapC was

49

able to suppress DNA-binding ability of ComA and the phosphorylation modification

50

of ComA (ComA~P), which combines with target DNA and increases the expression

51

of surfactin biosynthetic genes

52

adjust several subpopulations differentiation in B. subtilis 14. In addition, the study of

53

Yang et al.

54

sporulation, and competent cells in Bacillus sp. However, few studies have explored

55

the regulation of antimicrobial lipopeptides except for surfactin by the Rap-Phr

56

system.

15

12, 13

. Recently, a system of RapP-PhrP was used to

revealed that RapQ-PhrQ system could regulate surfactin production,

57

Here, the whole genome sequence of B. amyloliquefaciens fmbJ was reported.

58

The strain has remarkable variations in colonial feature with other Bacillus strains

59

isolated in our laboratory. Moreover, we found that there are gene clusters of

60

bacillomycin D in fmbJ genome which are different from B. subtilis 168.

61

Bacillomycin D as a part of the iturin family’s lipopeptide is composed by one

62

β-amino fatty acid and seven α-amino acids

63

antimicrobial bioactivity

64

agriculture. However, the production of bacillomycin D in wild strain is too low to

16

. It was reported to have high

5-7

. Therefore, it has a promising potential application in

ACS Paragon Plus Environment

Page 4 of 32

Page 5 of 32

Journal of Agricultural and Food Chemistry

65

satisfy the production and application. To improve the bacillomycin D yield, we

66

successfully deleted the rapC gene in fmbJ. Furthermore, we present several evidence

67

that the RapC-PhrC can regulate production of bacillomycin D in B.

68

amyloliquefaciens.

69

2. Materials and methods

70

2.1 Strain and cultivation

71

B. amyloliquefaciens fmbJ strain (used name B. subtilis fmbJ) was isolated and

72

characterized in our laboratory. The strain was stored in the China General

73

Microbiological Culture Collection Center as Bacillus sp. (CGMCC 0943). For DNA

74

analysis, the bacteria were cultivated in Luria-Bertani (LB) medium containing 0.5

75

g/L yeast extract, 1 g/L tryptone, and 1 g/L NaCl (pH 7.0) at 37 °C with 180 rpm for

76

24 h. The strain was inoculated in beef extract medium (3% beef extract, 10% peptone,

77

and 5% NaCl (pH 7.2)) and cultivated at 37 °C with 180 rpm as a pre-culture. The

78

Landy medium (20 g/L glucose, 5 g/L L-glutamic acid, 1 g/L yeast extract, 1.0 g/L

79

KH2PO4, 0.5 g/L MgSO4·7H2O; 0.5 g/L KCl, 5.0 mg/L MnSO4, 0.16 mg/L

80

CuSO4·5H2O, 0.15 mg/L FeSO4·7H2O, pH 7.0) was used as a fermentation medium

81

for the preparation of bacillomycin D. The fermentation conditions were 33 °C, 180

82

rpm, and for 72 h. Escherichia coli DH5α was severed as a host for plasmid

83

replication and E. coli JM110 (dam-/dcm-) was utilized for demethylation. They were

84

cultivated in LB medium at 37 °C. If necessary, 100 µg/mL ampicillin and 5 µg/mL

85

erythromycin were added for E. coli and B. amyloliquefaciens, respectively. In the

86

study, all the strains and plasmids are summarized in Table 1. All the chemicals and

ACS Paragon Plus Environment

Journal of Agricultural and Food Chemistry

87

culture medium compositions were bought from Sinopharm Chemical Reagent Co.,

88

Ltd, Nanjing, China.

89

2.2 Whole genome sequencing and assembly

90

The complete genomic sequencing of the fmbJ was carried out by the Beijing

91

Genomic Institute (BGI, Shenzhen, China) by a union of randomly sheared libraries.

92

Illumina Hiseq 4000 sequencing platform (Illumina; CA, USA) and PacBio RSII

93

sequencing platform (Pacific Biosciences, USA) were applied to perform the genomic

94

DNA sequencing, and the evaluations of all generated reads were qualitative 17, 18. De

95

novo assembly was carried out with the Short Oligonucleotides Alignment Program

96

(SOAP) denovo_v2.04 using the clean data

97

was used to estimate the size of genome, the degree of heterozygosis and the degree

98

of duplication. The result exhibited that the genome size of fmbJ was 4.28 Mb (Fig.

99

S1). After genome assembling, the GC distribution of fmbJ was obtained by

100

GC-Depth analysis (Fig. S2).

101

2.3 Genome analysis and annotation

19, 20

. Before assembling, k-mer analysis

102

Putative genes were accomplished through Glimmer_v3.02 to identification 21, 22.

103

All putative genes were checked with databases to acquire their consistent annotations.

104

In order to make sure the biological significance, the best alignment result was

105

selected as annotation. Functional annotation was finished via BLAST with different

106

public databases. Repeating DNA sequences were identified using Tandem Repeat

107

Finder (version 4.04, http://tandem.bu.edu/trf/trf.html)

108

minisatellite DNA and microsatellite DNA were determined by the length and number

ACS Paragon Plus Environment

23

. The selection of the

Page 6 of 32

Page 7 of 32

Journal of Agricultural and Food Chemistry

109

of repeat units. The tRNA genes and the tRNA secondary structure were predicted

110

using tRNAscan-SE 24, and the genes of rRNA were identified via BLAST according

111

to the current rRNA database. The result of rRNA database blasting is precise but not

112

complete, and rRNAmmer was used to predict rRNA when there was no homology

113

reference 25. The sRNA genes were identified using Rfam database 26. The structural

114

variation in fmbJ genome was determined according to that of the reference bacterium

115

based on Mummer

116

information is summarized in Fig. S3, S4. Meanwhile, the Core-Pan genes in fmbJ

117

and four other Bacillus. strains (B. amylpliquefaciens FZB42, B. amyloliquefaciens

118

NAU-B3, B. amyloliquefaciens Y2, and B. subtilis 168) were confirmed using

119

BLAST based on the method of Qin et al. 28.

120

2.4 Lipopeptide extraction and HPLC/ESI/CID-MS analysis

27

, including amino acid level and nucleic acid level. Detailed

121

After centrifugation, the fermentation broth supernatant was collected, and its pH

122

was adjusted to 2.0 by 4 M HCl. Then, the solution was stored at 4 °C until further

123

treatment. Subsequently, the precipitation was collected using centrifugation and

124

discarding the supernatant. The right amount of methanol was added to the

125

precipitation to extraction of the substance, and the pH was adjusted to 7.0. After

126

10,000 g × 10 min centrifugation, the crude products of bacillomycin D were obtained.

127

Bacillomycin D crude products were separated by Sephadex LH-20 column.

128

Preparation HPLC (Waters 600, USA) with a C18 column (Eclipse XDB, 5µm

129

4.6×250 mm, Agilent, USA) was used for further purification. Water with 0.1%

130

trifluoroacetic acid (TFA) and acetonitrile with 0.1% TFA were utilized for moving

ACS Paragon Plus Environment

Journal of Agricultural and Food Chemistry

Page 8 of 32

131

phase. The structure of bacillomycin D was identified by high performance liquid

132

chromatography/electrospray

133

spectrometry (HPLC/ESI/CID-MS) and Surveyor-LCQ DECA XP Plus of Thermo

134

Finnegan (Thermo Electron Corporation, San Jose, CA, USA). Specific HPLC and

135

MS conditions were carried out according to Gong et al. 7 and Qian et al. 29.

136

2.5 Strain construction

ionization/collision-induced

dissociation-mass

137

The temperature-sensitive vector pCBS was used to delete the rapC gene of

138

strain fmbJ. The upstream/ downstream regions of rapC were amplified with two

139

pairs of primers rapC-P1: 5′-gtcgacGAAGAAACGAAGCGGATG (SalI restriction

140

site

141

5′-CAAATAACAAACCATTCCTTCACCCTCCCCATCCA,

142

GGGGAGGGTGAAGGAATGGTTTGTTATTTGTTTAG,

143

agatctGCAGGAACTTCAAGCAGA (BglII restriction site underlined). The two

144

fragments were spliced via splicing by overlap extension (SOE) PCR, followed by

145

insertion of the spliced fragments (1407 bp) into pCBS using the SalI/ BamHI

146

restriction sites (because BamHI and BglII are the same tailed enzymes), resulting in

147

pCBS△rapC. The recombined plasmid pCBS△rapC was transferred into fmbJ

148

successfully through electro transformation technology

149

the rapC gene in B. amyloliquefaciens fmbJ was conducted according to the

150

previously procedure reported

151

selected. Then the rapC-deleted mutants were validated by PCR amplification and

152

sorting by the primer pair rapC-P1 and rapC-P4.

underlined),

rapC-P2: rapC-P3: and

rapC-P4:

5′5′-

30

. A marker-free deletion of

30

. Colonies without erythromycin resistance were

ACS Paragon Plus Environment

Page 9 of 32

Journal of Agricultural and Food Chemistry

153

2.6 Reverse transcription quantitative real-time PCR (RT-qPCR)

154

Total RNA was extracted with Trizol Reagent (TranGen Biotech, Beijing, China)

155

according to the protocol of manufacturer. RNA quality was examined by

156

electrophoresis on a 2.0% agarose gel, and its amount was analyzed by a

157

spectrophotometer (NanoDrop 2000, Thermo Scientific, USA). For RT-qPCR analysis,

158

1 µg RNA sample was used for cDNA synthesis with 5X All-In-One RT MasterMix

159

(AccuRT Genomic DNA Removal Kit; Applied Biological Materials Inc., Canada)

160

according to the protocol of manufacturer. Real-time PCR was carried out with a

161

mixture containing 1 µL cDNA, 0.2 µM forward primer, 0.2 µM reverse primer, 10 µL

162

Hieff™ qPCR SYBR Green Master Mix (High Rox Plus) (Yeasen, Shanghai, China),

163

and ddH2O in 20 µL of total volume. DNA was amplified with Real-Time PCR

164

System (StepOnePlus™,Applied Biosystems, USA) to analyze the expression of target

165

genes (bymA, bymB, bymC, bymD, TE, comA, comP, and phrC) under the PCR

166

procedure below: denaturation 5 min at 95 °C and 40 cycles of 95 °C for 10 s, 60 °C

167

for 30 s. The primers used for amplification of the reference gene 16sRNA and target

168

genes were presented in Table S1. The relative fold change of the target genes

169

expression was evaluated by the calculation of the 2-△△Ct. The threshold cycle (Ct)

170

values were obtained by the Real-Time PCR System software (StepOnePlusTM,

171

Applied Biosystems,USA)) 31.

172

2.7 Statistical analysis

173

All statistical analyses were carried out with one-way analysis of variance

174

(ANOVA) by SPSS (SPSS version 17.0, IBM, USA). After checking analysis results

ACS Paragon Plus Environment

Journal of Agricultural and Food Chemistry

175

of ANOVA, the p-value was given. Duncan's test was used to examine significance

176

below 0.05. Means ± standard deviation (SD) and triplicate of assay were used to

177

express all results.

178

3. Results

179

3.1. General genomic features of B. amyloliquefaciens fmbJ

180

The complete genomic sequencing of B. amyloliquefaciens fmbJ was carried out

181

by Illumina Hiseq 4000 sequencing platform to produce 605 Mb clean data and

182

Pacbio RSII sequencing platform to produce 594 Mb clean data. Based on the

183

assemble result of sample fmbJ, we found that the genome contains a 4,193,344 bp

184

circular chromosome, with 45.98% G+C (Fig. 1). With 134 tandem repeats (including

185

103 minisatellite DNA and 10 microsatellite DNA), 11 small RNA (sRNA), 87

186

transfer RNA genes, and 27 rRNA genes (Table S2), the whole genome of fmbJ was

187

predicted to contain 4,249 genes. Compared with B. amylpliquefaciens FZB42, fmbJ

188

may have a much more complicated gene regulation, because it has more number of

189

genes than FZB42 (3892 genes) (Table 2). Although there was no obvious difference

190

in the size of genome, to find the functional differences and similarities among five

191

strains, the core and pan genes were analyzed by taking the genome of B.

192

amyloliquefaciens FZB42 as a reference, B. amyloliquefaciens NAU-B3, B.

193

amyloliquefaciens Y2, B. subtilis 168, and fmbJ as query

194

fmbJ genome, the core genes organized by the five Bacillus strains consisted of 3,342

195

genes were sharing more than 50% identity to each other (Fig. 2a). They are necessary

196

for growth. Some genes are special genes when they are contained only by one of the

ACS Paragon Plus Environment

28

. Out of 4,249 genes in

Page 10 of 32

Page 11 of 32

Journal of Agricultural and Food Chemistry

197

bacteria. 5,143 pan genes were detected in five strains (Fig. 2b). Furthermore, we

198

analyzed the differences in five strains after removing the core genes (Fig. 3). The

199

results indicated that the difference between fmbJ and B. subtilis 168 was the largest,

200

with little difference among fmbJ, B. amyloliquefaciens NAU-B3 and B.

201

amyloliquefaciens Y2. The same results were found in fmbJ genomic structural

202

variation including nucleic acid level (Fig. S3) and amino acid level (Fig. S4).

203

After annotation, 2,956 CDSs (coding sequences) were assigned to putative

204

biological function, while 1,293 CDSs were considered as unknown function proteins.

205

For all CDSs without allocation functions, 1,238 CDSs were consistent with

206

conserved hypothetical proteins, whereas 55 CDSs were not homologous with any

207

previously reported sequences. By the analysis for the Cluster of Orthologous Groups

208

(COG), 3,424 CDSs were allocated to one or more COG functional groups. In

209

addition, 64 genes were related to controlling the bacterial mobility, 213 genes were

210

connected with the biosynthesis of membrane and cell wall, and 121 genes were

211

associated with transport, catabolism and the biosynthesis of secondary metabolites

212

(Fig. 4). The genome encodes a large number of pathways, including the production

213

of large amounts of antimicrobial substances. These genes may be favorable to

214

promote the growth of fmbJ strain and to protect its plant host from pathogens.

215

3.2. Cluster of non-ribosomal biosynthesized bacillomycin D in fmbJ

216

Three gene clusters (including bacillomycin D, fengycin and surfactin) for the

217

non-ribosomal biosynthesis of lipopeptides were found in fmbJ genome. The 37,250

218

bp bmy gene cluster was inserted in the fmbJ genome, which comprised four genes

ACS Paragon Plus Environment

Journal of Agricultural and Food Chemistry

219

(bmyD, bmyA, bmyB, bmyC). By compared the genome sequence of B.

220

amyloliquefaciens FZB42 (a strain of bacillomycin D production) with fmbJ, it was

221

found that bmyD (002202), bmyA (002201), bmyB (002200), and bmyC (002199) are

222

very similar to those of FZB42. BmyD, bmyA, bmyB, and bmyC are the synthetase

223

genes encoding bacillomycin D. The first gene bmyD codes a putative malonyl

224

coenzyme A transacylase which is nearly identical to the function of FabD in fatty

225

acid synthesis

226

synthetase, which are constituted as their relatives in the operons of mycosubtilin and

227

iturin A 32. A thioesterase (TE) domain is at the C-terminal end of bmyC, which has

228

the function of bacillomycin D cyclization and release.

229

3.3. HPLC-ESI mass spectrometry identification of bacillomycin D

16

. BmyA-C code the peptide forming subunits of bacillomycin D

230

Bacillomycin D was prepared by Sephadex LH-20 chromatography combined

231

with HPLC. Its liquid chromatographic analysis showed that 11 active peaks were

232

obtained in Fig. 5. These active peaks were further analyzed by MS/MS. The

233

corresponding compounds were confirmed by match up with the MS data of

234

bacillomycin D

235

contained seven peptide structures (-Asn-Tyr-Asn-Pro-Glu-Ser-Thr-). The MS of

236

peaks 1-3 displayed their [M + H]+ ion peaks at m/z 1030.60, 1031.30, and 1032.30,

237

respectively. Through analysis and comparison, we found that these peaks could be

238

identified to bacillomycin D of C14 fatty acid chain. The MS of C15 bacillomycin D

239

gave two [M+H]+ ion peaks at m/z 1045.30 (peaks 4) and 1046.31(peaks 5). The MS

240

of peaks 6-9 were confirmed as C16 bacillomycin D, which showed [M+H]+ ion peaks

7, 33, 34

. The results (Fig. S5) showed that the isolated bacillomycin D

ACS Paragon Plus Environment

Page 12 of 32

Page 13 of 32

Journal of Agricultural and Food Chemistry

241

at m/z 1058.76, 1059.18, 1060.20, and 1061.50. In addition, the MS of C17

242

bacillomycin D gave two [M+H]+ ion peaks at m/z 1072.91 (peak 10) and

243

1073.61(peak 11). These molecular weight intervals of 14 (-CH2-) were homologues

244

of bacillomycin D 7.

245

3.4. The influence of rapC on bacillomycin D biosynthesis

246

In Bacillus, the phosphatase RapC was found to involve in the biosynthesis of

247

lipopeptide surfactin 35. Bacillomycin D and surfactin are lipopeptides, and does rapC

248

gene also regulate bacillomycin D biosynthesis? Furthermore, we found that there is

249

rapC (003893) gene in fmbJ from the results of genome sequencing. Subsequently,

250

the rapC gene in fmbJ was successfully knocked out using marker-free biological

251

technique. In this study, the lipopeptides collected from fmbJ strain and its mutant

252

fmbJ△rapC strain were detected by HPLC. Fig. 6 showed the production of

253

bacillomycin D between in fmbJ strain and its mutant fmbJ△rapC in 24-120 h. As

254

shown in Fig. 6, there was no obvious difference between the bacillomycin D

255

production of the deletion of rapC mutant and that of wild strain during 48 h.

256

However, the bacillomycin D production of the mutant was remarkably higher than

257

that of fmbJ after 72 h (Fig. 6). Its maximum yield was 360.8 mg/L, which was 1.5

258

times higher than the original yield of 240.7 mg/L. These results indicated that the

259

rapC gene could negatively regulate the synthesis of bacillomycin D.

260

3.5. The relationship between rapC and bacillomycin D genes and related signal

261

factors

262

To further validate the regulation of bacillomycin D production by rapC,

ACS Paragon Plus Environment

Journal of Agricultural and Food Chemistry

Page 14 of 32

263

RT-qPCR was used to detect the expression of bacillomycin D synthetase genes

264

(bmyA, bmyB, bmyC, bmyD and TE) and relative regulatory signal genes (comA,

265

comP and phrC) in fmbJ and its mutant. The results were displayed in Fig. 7A-D,

266

which presented that genes expression had a certain time-space properties. As shown

267

in Fig. 7A, compared with fmbJ strain, the expression of bacillomycin D synthetase

268

genes was decreased in fmbJ△rapC, except for bmyC and TE at 12h, whereas the

269

signal gene comP increased significantly at this time, 1.8 times as much as fmbJ strain.

270

As can be seen from Fig. 7B, TE and signal genes have increased significantly,

271

especially that the expression of the comA gene was 18.9 times as high as that of fmbJ

272

at 24 h. All the bacillomycin D synthetase genes and signal genes were significantly

273

increased at 36 h and 48 h (Fig. 7 C, D). Bacillomycin D was a kind of secondary

274

metabolites which was produced after the bacterial growth to a certain stage. It was

275

also explained that rapC knockout improved bacillomycin D synthase genes

276

expression in the late stage of strain growth, thus increasing bacillomycin D yield.

277

This confirmed that the rapC gene negatively regulated the expression of

278

bacillomycin D synthase genes in strain fmbJ.

279

4. Discussion

280

The genomic sequencing of fmbJ provides the theoretical basis for

281

biotechnological means to improve the bacillomycin D production. Bacillomycin D

282

that composed of 14–17 carbon atoms of β-amino fatty acid chain and a seven amino

283

acids of cyclic peptide, had high antifungal and antitumor activities

284

the signal proteins (such as ComA, ComP, PhrC and RapC) were also closely related

ACS Paragon Plus Environment

5, 7, 29, 36

. Some of

Page 15 of 32

Journal of Agricultural and Food Chemistry

16, 29, 37

285

to the biosynthesis of lipopeptides

. The fmbJ genome showed that as a major

286

transcriptional regulator, comA (001020) was probably to be related to histidine

287

kinases comP (001019) through its phosphorylated form comA~P 9. PhrC (003892),

288

as a regulator of phosphatase rapC (003893), plays an inhibitory role 11. As shown in

289

Fig. 7, the deletion of rapC gene led to an increase in the expression level of three

290

genes (comA, comP and phrC) in varying degrees, compared to the original strain.

291

The systems of Rap-Phr have been found to be take part in biosynthesis of

292

lipopeptides through the mutual effect with major regulatory factor (ComA)15. When

293

the density of cell population reaches a very high level, Phr peptides are re-entered

294

inside the cell by the oligopeptide permease (Opp), then they inhibit cognate Rap

295

proteins activities

296

involved in regulation of bacillomycin D production in B. amyloliquefaciens fmbJ.

297

The finding was in close agreement to Yang’s researchresult that RapQ-PhrQ system

298

could regulate the production of surfactin in B. subtilis OKB105

299

transcriptional regulatory factor was related to the regulation of lipopeptide

300

production and sporulation in Bacillus

301

ComA enable to initiate the srfA operon expression by binding with promotor, so the

302

synthesis of surfactin was promoted

303

with ComA restrained the ability of response regulator to bind its target DNA

304

promoter. It was found in the work that the knockout of rapC gene increased the

305

expression of comA and comP genes, which would be beneficial for starting the

306

transcription and for improving the transcriptional ability of bacillomycin D

38

. Our results revealed that the RapC-PhrC system might be

13

15

. ComA, a

15, 37

. ComA~P, the product phosphorylated

. Core et al.

12

reported that RapC interaction

ACS Paragon Plus Environment

Journal of Agricultural and Food Chemistry

307

synthetase genes, indicated the bacillomycin D synthesis was of similar regulation

308

mode to surfactin production.

309

In summary, whole genome sequencing shows that B. amyloliquefaciens fmbJ

310

has a good prospect of bio-control and could produce many kinds of bioactive

311

secondary metabolites, especially bacillomycin D. To improve bacillomycin D yield,

312

the rapC regulation gene in the wild strain fmbJ was successfully knocked out by

313

marker-free method. It was found by deletion of RapC gene that RapC-PhrC system

314

in B. amyloliquefaciens fmbJ are involved in the regulation of bacillomycin D and

315

regulate the expression of comA, comP, and phrC genes, as well bacillomycin D

316

production of rapC-null strain was significantly increased. Further studies using the

317

multi-omics and biotransformation techniques of B. amyloliquefaciens fmbJ on the

318

basis of genome sequencing will contribute to clarify the function of the putative

319

genes and the synthetic pathways of bioactive substances.

320 321

Acknowledgements

322

The work was supported by the National Natural Science Foundation of China

323

(No. 31571887), Agricultural Innovation Foundation of Jiangsu Province (CX

324

16-1058), Jiangsu Collaborative Innovation Center (2011 plan) of Meat Production

325

and Processing, Quality and Safety Control.

326 327 328

ACS Paragon Plus Environment

Page 16 of 32

Page 17 of 32

Journal of Agricultural and Food Chemistry

329

References

330

1.

331

the Council of 21 October 2009 establishing a framework for Community action to

332

achieve the sustainable use of pesticides. In 2009; Vol. 309, pp 71-86.

333

2.

334

root-associated Bacillus amyloliquefaciens FZB42 – a review. Frontiers in

335

Microbiology 2015, 6, 780.

336

3. Wang, B. B.; Shen, Z. Z.; Zhang, F. G.; Waseem, R. Z.; Yuan, J.; Huang, R.; Ruan,

337

Y. Z.; Li, R.; Shen, Q. R., Bacillus amyloliquefaciens strain W19 can promote growth

338

and yield and suppress fusarium wilt in banana under greenhouse and field conditions.

339

Pedosphere 2016, 26, 733-744.

340

4.

341

Richter, T.; Borriss, R., Extracellular phytase activity of Bacillus amyloliquefaciens

342

FZB45 contributes to its plant-growth-promoting effect. Microbiology 2002, 148,

343

2097-2109.

344

5.

345

inhibition of Fusarium graminearum and reduction of deoxynivalenol production in

346

wheat grain by bacillomycin D. Journal of Stored Products Research 2018, 75, 21-28.

347

6.

348

activity mode of Aspergillus ochraceus by bacillomycin D and its inhibition of

349

ochratoxin A (OTA) production in food samples. Food Control 2016, 60, 281-288.

350

7.

European Parliament, Directive 2009/128/EC of the European Parliament and of

Chowdhury, S. P.; Hartmann, A.; Gao, X.; Borriss, R., Biocontrol mechanism by

Idriss, E. E.; Makarewicz, O.; Farouk, A.; Rosner, K.; Greiner, R.; Bochow, H.;

Sun, J.; Li, W.; Liu, Y.; Lin, F.; Huang, Z.; Lu, F.; Bie, X.; Lu, Z., Growth

Qian, S.; Lu, H.; Sun, J.; Zhang, C.; Zhao, H.; Lu, F.; Bie, X.; Lu, Z., Antifungal

Gong, Q.; Zhang, C.; Lu, F.; Zhao, H.; Bie, X.; Lu, Z., Identification of

ACS Paragon Plus Environment

Journal of Agricultural and Food Chemistry

351

bacillomycin D from Bacillus subtilis fmbJ and its inhibition effects against

352

Aspergillus flavus. Food Control 2014, 36, 8-14.

353

8.

354

Bacillus subtilis fmbJ on apoptosis and necrosis in Rhizopus stolonifer. J Microbiol.

355

2014, 52, 675-680.

356

9.

357

fates in Bacillus subtilis. FEMS Microbiology Reviews 2010, 34, 134-149.

358

10. Auchtung, J. M.; Lee, C. A.; Grossman, A. D., Modulation of the

359

ComA-dependent quorum response in Bacillus subtilis by multiple Rap proteins and

360

Phr peptides. Journal of Bacteriology 2006, 188, 5273-5285.

361

11. Bongiorni, C.; Ishikawa, S.; Stephenson, S.; Ogasawara, N.; Perego, M.,

362

Synergistic regulation of competence development in Bacillus subtilis by two Rap-Phr

363

systems. Journal of Bacteriology 2005, 187, 4353-4361.

364

12. Core, L.; Perego, M., TPR-mediated interaction of RapC with ComA inhibits

365

response regulator-DNA binding for competence development in Bacillus subtilis.

366

Molecular Microbiology 2003, 49, 1509-1522.

367

13. Nakano, M. M.; Xia, L. A.; Zuber, P., Transcription initiation region of the srfA

368

operon, which is controlled by the comP-comA signal transduction system in Bacillus

369

subtilis. Journal of Bacteriology 1991, 173, 5487-5493.

370

14. Parashar, V.; Konkol, M. A.; Kearns, D. B.; Neiditch, M. B., A plasmid-encoded

371

phosphatase regulates Bacillus subtilis biofilm architecture, sporulation, and genetic

372

competence. Journal of Bacteriology 2013, 195, 2437-2448.

Tang, Q.; Bie, X.; Lu, Z.; Lv, F.; Tao, Y.; Qu, X., Effects of fengycin from

López, D.; Kolter, R., Extracellular signals that define distinct and coexisting cell

ACS Paragon Plus Environment

Page 18 of 32

Page 19 of 32

Journal of Agricultural and Food Chemistry

373

15. Yang, Y.; Wu, H.-J.; Lin, L.; Zhu, Q.-q.; Borriss, R.; Gao, X.-W., A plasmid-born

374

Rap-Phr system regulates surfactin production, sporulation and genetic competence in

375

the heterologous host, Bacillus subtilis OKB105. Applied Microbiology and

376

Biotechnology 2015, 99, 7241-7252.

377

16. Chen, X. H.; Koumoutsi, A.; Scholz, R.; Schneider, K.; Vater, J.; Süssmuth, R.;

378

Piel, J.; Borriss, R., Genome analysis of Bacillus amyloliquefaciens FZB42 reveals its

379

potential for biocontrol of plant pathogens. Journal of Biotechnology 2009, 140,

380

27-37.

381

17. Li, X.; Kui, L.; Zhang, J.; Xie, Y.; Wang, L.; Yan, Y.; Wang, N.; Xu, J.; Li, C.;

382

Wang, W.; van Nocker, S.; Dong, Y.; Ma, F.; Guan, Q., Improved hybrid de novo

383

genome assembly of domesticated apple (Malus x domestica). GigaScience 2016, 5,

384

35.

385

18. Quail, M. A.; Smith, M.; Coupland, P.; Otto, T. D.; Harris, S. R.; Connor, T. R.;

386

Bertoni, A.; Swerdlow, H. P.; Gu, Y., A tale of three next generation sequencing

387

platforms: comparison of Ion Torrent, Pacific Biosciences and Illumina MiSeq

388

sequencers. BMC Genomics 2012, 13, 341.

389

19. Li, R.; Zhu, H.; Ruan, J.; Qian, W.; Fang, X.; Shi, Z.; Li, Y.; Li, S.; Shan, G.;

390

Kristiansen, K., De novo assembly of human genomes with massively parallel short

391

read sequencing. Genome research 2010, 20, 265-272.

392

20. Li, R.; Li, Y.; Kristiansen, K.; Wang, J., SOAP: short oligonucleotide alignment

393

program. Bioinformatics 2008, 24, 713-714.

394

21. Delcher, A. L.; Bratke, K. A.; Powers, E. C.; Salzberg, S. L., Identifying bacterial

ACS Paragon Plus Environment

Journal of Agricultural and Food Chemistry

395

genes and endosymbiont DNA with Glimmer. Bioinformatics 2007, 23, 673-679.

396

22. Delcher, A. L.; Harmon, D.; Kasif, S.; White, O.; Salzberg, S. L., Improved

397

microbial gene identification with Glimmer. Nucleic Acids Research 1999, 27,

398

4636-4641.

399

23. Benson, G., Tandem repeats finder: a program to analyze DNA sequences.

400

Nucleic Acids Research 1999, 27, 573-80.

401

24. Lowe, T. M.; Eddy, S. R., tRNAscan-SE: a program for improved detection of

402

transfer RNA genes in genomic sequence. Nucleic Acids Research 1997, 25, 955-64.

403

25. Lagesen, K.; Hallin, P.; Rodland, E. A.; Staerfeldt, H. H.; Rognes, T.; Ussery, D.

404

W., RNAmmer: consistent and rapid annotation of ribosomal RNA genes. Nucleic

405

Acids Research 2007, 35, 3100-8.

406

26. Gardner, P. P.; Daub, J.; Tate, J. G.; Nawrocki, E. P.; Kolbe, D. L.; Lindgreen, S.;

407

Wilkinson, A. C.; Finn, R. D.; Griffiths-Jones, S.; Eddy, S. R.; Bateman, A., Rfam:

408

updates to the RNA families database. Nucleic Acids Research 2009, 37, D136-D140.

409

27. Kurtz, S.; Phillippy, A.; Delcher, A. L.; Smoot, M.; Shumway, M.; Antonescu, C.;

410

Salzberg, S. L., Versatile and open software for comparing large genomes. Genome

411

Biology 2004, 5, R12.

412

28. Qin, J.; Li, R.; Raes, J.; Arumugam, M.; Burgdorf, K. S.; Manichanh, C.; Nielsen,

413

T.; Pons, N.; Levenez, F.; Yamada, T.; Mende, D. R.; Li, J.; Xu, J.; Li, S.; Li, D.; Cao,

414

J.; Wang, B.; Liang, H.; Zheng, H.; Xie, Y.; Tap, J.; Lepage, P.; Bertalan, M.; Batto,

415

J.-M.; Hansen, T.; Le Paslier, D.; Linneberg, A.; Nielsen, H. B.; Pelletier, E.; Renault,

416

P.; Sicheritz-Ponten, T.; Turner, K.; Zhu, H.; Yu, C.; Li, S.; Jian, M.; Zhou, Y.; Li, Y.;

ACS Paragon Plus Environment

Page 20 of 32

Page 21 of 32

Journal of Agricultural and Food Chemistry

417

Zhang, X.; Li, S.; Qin, N.; Yang, H.; Wang, J.; Brunak, S.; Dore, J.; Guarner, F.;

418

Kristiansen, K.; Pedersen, O.; Parkhill, J.; Weissenbach, J.; Bork, P.; Ehrlich, S. D.;

419

Wang, J., A human gut microbial gene catalogue established by metagenomic

420

sequencing. Nature 2010, 464, 59-65.

421

29. Qian, S.; Lu, H.; Meng, P.; Zhang, C.; Lv, F.; Bie, X.; Lu, Z., Effect of inulin on

422

efficient production and regulatory biosynthesis of bacillomycin D in Bacillus subtilis

423

fmbJ. Bioresource Technology 2015, 179, 260-267.

424

30. Zakataeva, N.; Nikitina, O.; Gronskiy, S.; Romanenkov, D.; Livshits, V., A simple

425

method to introduce marker-free genetic modifications into the chromosome of

426

naturally nontransformable Bacillus amyloliquefaciens strains. Applied Microbiology

427

and Biotechnology 2010, 85, 1201-1209.

428

31. Wang, P.; Guo, Q.; Ma, Y.; Li, S.; Lu, X.; Zhang, X.; Ma, P., DegQ regulates the

429

production of fengycins and biofilm formation of the biocontrol agent Bacillus

430

subtilis NCD-2. Microbiological Research 2015, 178, 42-50.

431

32. Moyne, A. L.; Cleveland, T. E.; Tuzun, S., Molecular characterization and

432

analysis of the operon encoding the antifungal lipopeptide bacillomycin D. FEMS

433

Microbiology Letters 2004, 234, 43-49.

434

33. Tanaka, K.; Ishihara, A.; Nakajima, H., Isolation of anteiso-C-17, iso-C-17,

435

iso-C-16, and iso-C-15 bacillomycin D from Bacillus amyloliquefaciens SD-32 and

436

their antifungal activities against plant pathogens. Journal of Agricultural and Food

437

Chemistry 2014, 62, 1469-1476.

438

34. Xu, Z.; Shao, J.; Li, B.; Yan, X.; Shen, Q.; Zhang, R., Contribution of

ACS Paragon Plus Environment

Journal of Agricultural and Food Chemistry

439

bacillomycin D in Bacillus amyloliquefaciens SQR9 to antifungal activity and biofilm

440

formation. Applied and Environmental Microbiology 2013, 79, 808-815.

441

35. Duitman, E. H.; Wyczawski, D.; Boven, L. G.; Venema, G.; Kuipers, O. P.;

442

Hamoen, L. W., Novel methods for genetic transformation of natural Bacillus subtilis

443

isolates used to study the regulation of the mycosubtilin and surfactin synthetases.

444

Applied and Environmental Microbiology 2007, 73, 3490-3496.

445

36. Hajare, S. N.; Subramanian, M.; Gautam, S.; Sharma, A., Induction of apoptosis

446

in human cancer cells by a Bacillus lipopeptide bacillomycin D. Biochimie 2013, 95,

447

1722-1731.

448

37. Koumoutsi, A.; Chen, X.-H.; Vater, J.; Borriss, R., DegU and YczE positively

449

regulate the synthesis of bacillomycin D by Bacillus amyloliquefaciens strain FZB42.

450

Applied and Environmental Microbiology 2007, 73, 6953-6964.

451

38. Tjalsma, H.; Bolhuis, A.; Jongbloed, J. D. H.; Bron, S.; van Dijl, J. M., Signal

452

peptide-dependent protein transport in Bacillus subtilis: a genome-based survey of the

453

secretome. Microbiology and Molecular Biology Reviews 2000, 64, 515-547.

454

ACS Paragon Plus Environment

Page 22 of 32

Page 23 of 32

Journal of Agricultural and Food Chemistry

Figure Captions and Tables Fig.1 The circular map of the Bacillus amyloliquefaciens fmbJ genome for several specific genomic features. Innermost circle (1st): scale (bps); 2nd circle: GC skew-(violet); 3rd circle: GC skew+ (green); 4th circle: GC content (black); 5th circle: all genes are labeled color according to their function. Fig. 2 Dilution curve of bacterial genes. A, core genes; B, pan genes. Fig. 3 Heat-map after core gene deletion. Fig. 4 COG annotation of sample bacillus. fmbJ. Fig. 5 Analysis of bacillomycin D by HPLC. Peaks 1-11: bacillomycin D. Fig. 6 The influence of rapC gene knockout in fmbJ on bacillomycin D production. The strains fmbJ and fmbJ△rapC were cultured in100 mL Landy at 33°C for 24-120 h. * and ** were significantly different from controls at 0.05 and 0.01. Fig. 7 Effect of rapC gene knockout on relative expression of bacillomycin D synthase genes and signal genes. The strains fmbJ and fmbJ△rapC were cultured in 100 mL Landy at 33°C for 12 h, 24 h, 36 h, and 48 h (A-D). * and ** were significantly different from controls at 0.05 and 0.01.

Tables Table 1 Microorganisms and plasmids presented in the study. Table 2 General features of the genomes of the fmbJ and other Bacillus species

ACS Paragon Plus Environment

Journal of Agricultural and Food Chemistry

Page 24 of 32

Tables Table 1 Microorganisms and plasmids presented in the study. Strains or plasmids

Characteristics

Source

Wild-type strain, producer of bacillomycin

Laboratory stock

B. amyloliquefaciens strain fmbJ

D fmbJ∆rapC

rapC deletion mutant, derivative of strain

Current study

fmbJ Trans 5α (E. coli DH 5α)

F-φ80 lac Z∆M15 ∆(lacZYA-arg F) U169

TransGen Biotech

endA1 recA1 hsdR17(rk-,mk+) supE44λ- thi -1 gyrA96 relA1 phoA Trans 110 (E. coli JM110)

rpsL (Str R) thr leu thi-1 lacY galK galT ara

TransGen Biotech

tonA tsx dam dcm supE44 Δ (lac-proAB) /F′[traD36 proAB lacIqlacZΔM15] Plasmids pMAD

E.

coli

and

B.

subtilis

shuttle,

Laboratory stock

temperature-sensitive vector. Apr Emr (9666 bp) pCBS

pMAD

with

minor

modification.

The

Current study

3928-6049 bases in pMAD were removed, and Pamy, SamyE and lacZ were added in that location. Apr Emr (8102 bp) pCBS∆rapC

pCBS with rapC deletion box. Apr Emr

Apr, Emr indicate resistance to ampicillin and erythromycin, respectively.

ACS Paragon Plus Environment

Current study

Page 25 of 32

Journal of Agricultural and Food Chemistry

Table 2 General features of the genomes of the fmbJ and other Bacillus species Species Genome size (bp)

B. amyloliquefaciens

B. amyloliquefaciens

B. amyloliquefaciens

B. amyloliquefaciens

B. subtilis

fmbJ

FZB42

NAU-B3

Y2

168

4,193,344

3,918,589

4,204,608

4,238,624

4,215,606

gene

4,249

3,892

4,194

4,246

4,420

tRNA

87

88

91

86

86

rRNA

27

29

30

29

30

GC%

45.98

46.50

45.99

45.90

43.50

ACS Paragon Plus Environment

Journal of Agricultural and Food Chemistry

Fig. 1 The circular map of the Bacillus amyloliquefaciens fmbJ genome for several specific genomic features. Innermost circle (1st): scale (bps); 2nd circle: GC skew-(violet); 3rd circle: GC skew+ (green); 4th circle: GC content (black); 5th circle: all genes are labeled color according to their function. 101x76mm (300 x 300 DPI)

ACS Paragon Plus Environment

Page 26 of 32

Page 27 of 32

Journal of Agricultural and Food Chemistry

Fig. 2 Dilution curve of bacterial genes. A, core genes; B, pan genes. 168x84mm (300 x 300 DPI)

ACS Paragon Plus Environment

Journal of Agricultural and Food Chemistry

Fig. 3 Heat-map after core gene deletion. 82x77mm (300 x 300 DPI)

ACS Paragon Plus Environment

Page 28 of 32

Page 29 of 32

Journal of Agricultural and Food Chemistry

Fig. 4 COG annotation of sample bacillus. fmbJ. 101x72mm (300 x 300 DPI)

ACS Paragon Plus Environment

Journal of Agricultural and Food Chemistry

Fig. 5 Analysis of bacillomycin D by HPLC. Peaks 1-11: bacillomycin D. 101x65mm (300 x 300 DPI)

ACS Paragon Plus Environment

Page 30 of 32

Page 31 of 32

Journal of Agricultural and Food Chemistry

Fig. 6 The influence of rapC gene knockout in fmbJ on bacillomycin D production. The strains fmbJ and fmbJ△rapC were cultured in100 mL Landy at 33°C for 24-120 h. * and ** were significantly different from controls at 0.05 and 0.01. 212x174mm (300 x 300 DPI)

ACS Paragon Plus Environment

Journal of Agricultural and Food Chemistry

Fig. 7 Effect of rapC gene knockout on relative expression of bacillomycin D synthase genes and signal genes. The strains fmbJ and fmbJ△rapC were cultured in100 mL Landy at 33°C for 12 h, 24 h, 36 h, and 48 h (A-D). * and ** were significantly different from controls at 0.05 and 0.01. 229x170mm (300 x 300 DPI)

ACS Paragon Plus Environment

Page 32 of 32