【Why Do We Design This Experiment】

Bacillus subtilis has been widely applied as engineered bacteria, especially in food industry and pharmaceutical industry, for its safety and excellent secretion capacity. Therefore, after comparing characters of distinct mutants we selected Bacillus subtilis WB800N mutant as our engineered bacteria and looked up plenty of papers to select the optimal conditions for our experiment. To our disappointment, very few experiments have been done on WB800N mutant, and most optimization experiments regarding Bacillus subtilis focus solely on the optimization of production of specific proteins produced by Bacillus subtilis. Consider the final goal of our project, it is imperative to design this experiment on our own to find out the best condition for Bacillus subtilis WB800N.

【Methodology】

Any optimization designs will inevitably involve the ideology of Design of Experiment (DOE), which includes several dependent plots. Among them Orthogonal Design and Response Surface Design, RSM for short, are the most common two in biological experiments. Generally, Orthogonal Design consumes less time and has been used more widely, yet it is not logically rigorous in mathematics, and sometimes it overlooks interactions and alias between or among factors. In contrast, RSM is constructed on rigorous mathematical theories and excels in data analysis. Having weighing the features of the two methods carefully, we finally chose RSM.

【Sweeping Factors】

The first step of any methods of DOE is to investigate all variables that affect the results and select controllable factors for the experiment. In terms of this experiment, all factors can be categorized into two kinds: environment factors, like temperature, the rotation speed of the shaker, and the components of the medium. We have looked up several papers about the optimization experiments on Bacillus subtilis, finding the rotation speed of shakers ranging from 100 r/min to 250 r/min, and generally rotation speed only plays a tiny role. Additionally, our lab has only two shakers. While we can place twenty different mediums into one shaker at a time, we must run the shakers every time we alert the speed, which surely consumes longer time. Thus, we fixed the rotation speed of shakers at 200r/min.However, temperature and inoculation time are both vital environment factors whose effects cannot be ignored.
Inoculation amount and pack amount are also two factors that affect results slightly. We fixed them at 5 percent and 30mL/500mL respectively according to earlier authentic experiments.
A typical medium consists of carbon source, nitrogen source and inorganic salt, all of which are essential to ensure the regular metabolism of engineered bacteria. Finally in light of convenience, we infered the components of typical LB medium and determined three independent medium factors: peptone, yeast extract and sodium chloride (NaCl). Peptone provides nitrogen and carbon for the colonies, while yeast extract contains most required inorganic salt, therefore we did not list any inorganic salt except NaCl. We had no idea why NaCl is listed alone, and we suspected the influence of NaCl as yeast extract had already contains sodium.
Thus, we had five independent factors: temperature, inoculation time, peptone, yeast extract and NaCl. We further investigated some papers and defined their ranges. The following table displays their levels, and the unit of peptone, yeast extract and NaCl is g/L:

Factor	Low	High
Temperature	25℃	35℃
25℃	12h	24h
Peptone	5	15
Yeast Extract	2.5	7.5
NaCl	5	15

【Designs&Results】

The methodology of RSM can be divided into two subplots: Central Composite Designs (CCD) and Box-Behnken Designs. Generally the overall runs of Box-Behnken Designs is fewer when the factors are fixed, but Central composite designs are often recommended when the design plan calls for sequential experimentation because these designs can incorporate information from a properly planned factorial experiment. In our experiment, time is more precious than reagents, and as time itself is also an independent factor, Box-Behnken Designs would not have saved any time if adopted. Thus we selected CCD.
CCD itself can also be classified into three subplots: Central Composite Circumscribed Design (CCC), Central Composite Inscribed design(CCI) and Central Composite Face-centered Design(CCF). The alpha value of CCC is related to the number of factors, whereas in CCF α is fixed at 1, and only CCC is rotatable. The rotational invariance empowers CCC to be mathematically preferred, yet the value of alpha in a five-factor-CCC is over 2. In other words, if we adopted CCC, we would get some absurd treatments where the concentration of some specific actual material were negative. If we narrowed down the range to ensure any concentration is positive, the ranges of all three medium factors would be too narrow to yield cogent results. Therefore, we finally selected CCF.
We conducted our experiments according to the following table, which was calculated by Minitab, and the results, which were measure by OD value, were also included:

No.	Temperature	Time	Peptone	Yeast extract	NaCl	OD
1	25	12	5	2.5	5	0.511
2	35	24	5	2.5	5	1.625
3	35	12	15	2.5	5	2.783
4	25	24	15	2.5	5	1.74
5	35	12	5	7.5	5	2.317
6	25	24	5	7.5	5	2.4
7	25	12	15	7.5	5	0.912
8	35	24	15	7.5	5	3
9	35	12	5	2.5	15	2.169
10	25	24	5	2.5	15	1.77
11	25	12	15	2.5	15	0.371
12	35	24	15	2.5	15	2.7
13	25	12	5	7.5	15	0.754
14	35	24	5	7.5	15	2.58
15	35	12	15	7.5	15	3.128
16	25	24	15	7.5	15	2.38
17	30	18	10	5	10	2.908
18	30	18	10	5	10	2.908
19	30	18	10	5	10	1.75
20	30	18	10	5	10	2.908
21	35	12	5	2.5	5	2.082
22	25	24	5	2.5	5	1.75
23	25	12	15	2.5	5	0.508
24	35	24	15	2.5	5	2.6
25	25	12	5	7.5	5	0.989
26	35	24	5	7.5	5	2.8
27	35	12	15	7.5	5	2.782
28	25	24	15	7.5	5	1.7
29	25	12	5	2.5	15	0.508
30	35	24	5	2.5	15	1.338
31	35	12	15	2.5	15	3.061
32	25	24	15	2.5	15	2.2
33	35	12	5	7.5	15	2.167
34	25	24	5	7.5	15	1.53
35	25	12	15	7.5	15	0.555
36	35	24	15	7.5	15	2.9
37	30	18	10	5	10	2.908
38	30	18	10	5	10	2.908
39	30	18	10	5	10	2.908
40	30	18	10	5	10	2.957
41	25	18	10	5	10	1.907
42	35	18	10	5	10
43	30	12	10	5	10	2.652
44	30	24	10	5	10	2.908
45	30	18	5	5	10	2.726
46	30	18	15	5	10	3.042
47	30	18	10	2.5	10	2.598
48	30	18	10	7.5	10	3.124
49	30	18	10	5	5	2.999
50	30	18	10	5	15	2.834
51	30	18	10	5	10	2.908
52	30	18	10	5	10	2.908
53	30	18	10	5	10	2.908
54	30	18	10	5	10	2.908

The result of No.42 medium is destroyed due to some unfortunate reason. Additionally, multiple center points, which means conducting multiple experiments at the center points with identical treatments, is a very common phenomenon in DOE, yet we decided to do only experiment at the center point and reuse its result due to our limited time and reagents.
Estimated Regression Coefficients for OD

Term	Coef	SE Coef	T	P
Constant	2.87625	0.07126	40.361	0.000
Temperature	0.60225	0.05210	11.560	0.000
Time	0.28447	0.05072	5.608	0.000
Peptone	0.18665	0.05072	3.680	0.001
Yeast Extract	0.16776	0.05072	3.308	0.002
NaCl	-0.01626	0.05072	-0.321	0.751
Temperature*Temperature	-0.54900	0.24585	-2.233	0.033
Time*Time	-0.18725	0.19289	-0.971	0.339
Peptone*Peptone	-0.08325	0.19289	-0.432	0.669
Yeast Extract*Yeast Extract	-0.10625	0.19289	-0.551	0.586
NaCl*NaCl	-0.05075	0.19289	-0.263	0.794
Temperature*Time	-0.358338	0.05228	-6.579	0.000
Temperature*Peptone	0.17881	0.05228	3.420	0.002
Temperature*Yeast Extract	0.04544	0.05228	0.869	0.391
Temperature*NaCl	0.01550	0.05228	0.296	0.769
Time*Peptone	0.02575	0.05228	0.493	0.626
Time*Yeast Extract	.06112	0.05228	1.169	0.261
Time*NaCl	-0.00144	0.05228	-0.027	0.978
Peptone*Yeast Extract	-0.07469	0.05228	-1.429	0.163
Peptone*NaCl	0.09150	0.05228	1.750	0.090
Yeast Extract*NaCl	-0.04450	0.05228	-0.851	0.401

Term	Mark
OD	F
Temperature	T
Time	T
Peptone	P
Yeast Extract	Y
NaCl	C

StdOrder	OD	Fit	SE	Fit	Residual	StResid
1	1	1.740	1.757	0.207	-0.017	-0.08
2	2	0.912	0.612	0.207	0.300	1.42
3	3	1.770	1.563	0.207	0.207	0.98
4	4	2.317	2.593	0.208	-0.276	-1.31
5	5	2.169	1.958	0.208	0.211	1.00
6	6	3.000	2.979	0.208	0.021	0.10
7	7	1.750	2.876	0.071	-1.126	-3.92R
8	8	0.511	0.619	0.207	-0.108	-0.51
9	9	2.580	2.253	0.208	0.327	1.55
10	10	3.128	3.036	0.208	0.092	0.44
11	11	2.380	2.091	0.207	0.289	1.37
12	12	2.783	2.697	0.208	0.086	0.41
13	13	2.700	2.758	0.208	-0.058	-0.28
14	14	2.908	2.876	0.071	0.032	0.11
15	15	2.908	2.876	0.071	0.032	0.11
16	16	0.371	0.761	0.207	-0.390	-1.85
17	17	2.400	2.329	0.207	0.071	0.34
18	18	1.625	1.742	0.208	-0.117	-0.56
19	19	0.754	0.647	0.207	0.107	0.51
20	20	2.908	2.876	0.071	0.032	0.11
21	21	3.124	2.938	0.203	0.186	0.87
22	22	2.652	2.405	0.203	0.247	1.15
23	23	2.908	2.876	0.071	0.032	0.11
25	25	2.908	2.876	0.071	0.032	0.11
26	26	2.598	2.602	0.203	-0.004	-0.02
27	27	2.908	2.876	0.071	0.032	0.11
28	28	1.907	1.725	0.260	0.182	1.29
29	29	2.908	2.876	0.071	0.032	0.11
30	30	2.999	2.842	0.203	0.157	0.73
31	31	2.726	2.606	0.203	0.120	0.56
32	32	3.042	2.980	0.203	0.062	0.29
33	33	2.834	2.809	0.203	0.025	0.12
34	34	2.908	2.973	0.203	-0.065	-0.30
35	35	0.508	0.550	0.207	-0.042	-0.20
36	36	0.989	0.980	0.207	0.009	0.04
37	37	1.338	1.644	0.208	-0.306	-1.45
38	38	2.908	2.876	0.071	0.032	0.11
39	39	2.908	2.876	0.071	0.032	0.11
40	40	2.082	2.051	0.208	0.031	0.15
41	41	1.750	1.724	0.207	0.026	0.12
42	42	2.167	2.323	0.208	-0.156	-0.74
43	43	2.957	2.876	0.071	0.081	0.28
44	44	2.800	2.529	0.208	0.271	1.29
45	45	0.555	0.645	0.207	-0.090	-0.43
46	46	2.900	3.069	0.208	-0.169	-0.80
47	47	3.061	2.970	0.208	0.091	0.43
48	48	2.600	2.491	0.208	0.109	0.52
49	49	2.782	2.940	0.208	-0.158	-0.75
50	50	2.200	1.963	0.207	0.237	1.12
51	51	0.508	0.464	0.207	0.044	0.21
52	52	1.530	1.991	0.207	-0.461	-2.18R
53	53	1.700	2.064	0.207	-0.364	-1.73
54	54	2.908	2.876	0.071	0.032	0.11

S = 0.295758 PRESS = 7.78904
R-Sq = 92.25% R-Sq(pred) = 78.45% R-Sq(adj) = 87.41%
Suppose we redefine the factors accoding to the following table:
According to the ANOVA calculated by minitab, we got the expression of OD:
f(T,t,p,y,c)=2.87625+0.60225T+0.28447t+0.18665p+0.16776y-0.01626c-0.549T^2-0.18725t^2-0.08325p^2-0.10625y^2-0.05075c^2-0.35558Tt+0.17881Tp+0.04544Ty+0.0155Tc+0.02575tp+0.06112ty-0.00144tc-0.07469py+0.09150pc-0.04450yc
P represents confidence coefficient, which is a key judgment to check the reliability of the fitting function. In other words, if P=0.05, the probability that this term is wrong is 5%. The coefficient of determination (R) was calculated to be 0.9225, indicating that the model could explain 92% of the variability .From the above table we can identify eight statistically significant and reliable terms:
Constant;
Temperature;
Time;
Yeast Extract;
Peptone;
Temperature*Temperature;
Temperature*Time;
Temperature*Yeast Extract;
The influences of linear terms predominated, except NaCl, which substantiated our suspicion whereas most square terms and interaction terms were ignorable and statistically unreliable. Temperature and time and two most influential factor.
As our world is three-dimensional but the intact response surface is six-dimensional, it is impossible to draw the intact surface. Yet we could fix some factors to lower the dimensional, which empowers us to imagine the full surface. Here are some surfaces and contours of our fitting surface, we can extrapolate this super surface by combining these pictures:

【Designs&Results】

In our experiment , we discovered that the promoter P43 is too weak to express amilCP for a blue-color display. So we redesigned the circuits , aiming to find a more feasible plan, and test the part amilCP simultaneously. The new disign contains two circuits:

1 Pgrac-amilCP-Terminator

We use the Promoter grac , a promotor with lac promotor ,on the PHT vector to express amilCP, instead of Promoter 43. The Promoter grac is strong enough ,so we can easily see the results theoretically.

pgrac-amilCP overlap PCR（Right three tracks） PSBC3 Double Digestion（Left four tracks）

2 P43-amilcp-SigB-Terminator

We plan to use a positive feedback to magnify the expressing of amilCP. Because the P43 is a sigma B factor binding promotor, we designed a circuit, that the P43 is fused to the sigma B factor. We hope this could increase the response of P43.
Course of time limitation, we haven’t put this design to experiment.

【designing of the suicide system】

We design a circuit of killing switch based on its endogenous genetic system.
In B.subtilis, when it comes to the stationary phase, the environmental pressure increases and nutrition becomes limited, so B begin to produce spores. Now the community will be divided into two different parts. One of them are trying to kill others to get enough nutrient , delaying the production of spores and achieving a competitive advantage. Killing is mediated by the exported toxic protein SdpC. SdpI will appear on the membrane surface to avoid itself from being damaged. SdpI could bind free SdpC and autopressor SdpR, to remove SdpR’s inhibition against I and R, to produce more SdpI to offset SdpC, finally guaranteeing the subgroup alive, thereby delaying the spores production.
https://static.igem.org/mediawiki/2013/2/2b/Reporter_3.png We transfer SdpC which is fused by promoter SdpI/R into high copy plasmids in order to damage the balance of the system, thereby killing whole colony. When SdpC appears, SdpI on the membrane will bind free SdpC and adsorb SdpR to cease its inhibition against SdpI P/R, trying to produce more SdpI. At the same time, it will activate the promoter SdpR/I in our circuits and generate more SdpC.The system would fall into an infinite loop, and according to our modeling ,the amount of SdpC increases beyond the ability of SdpI.Thus,the cells with protection mechanism will crack and die because of too much SdpC. All above forms the killing device. We Also designed a test circuit,which contains promotor grac and sdpABC only,aiming to determine the ability of SdpC.

Results

https://static.igem.org/mediawiki/2013/e/ec/Reporter_4.tif
Colony PCR E.coli PHT43 + Promotor grac + SdpABC
https://static.igem.org/mediawiki/2013/9/92/Reporter_5.tif
PHT43 + Promotor SdpRI + SdpABC Enzyme digestion

【References】

Parallel pathways of repression and antirepression governing the transition to stationary phase in Bacillus subtilis AV Banse, A Chastanet, L Rahn-Lee…,PNAS ,2008

Reporter System

B.Subtilis Culture

Mice Modeling

Team:USTC CHINA/Modeling/B.SubtilisCulture

From 2013.igem.org

【Why Do We Design This Experiment】

【Methodology】

【Sweeping Factors】

【Designs&Results】

【Designs&Results】

1 Pgrac-amilCP-Terminator

2 P43-amilcp-SigB-Terminator

【designing of the suicide system】

Results

【References】