Results of System 1

In system one, we aim at boosting NAD+ levels for its critical role in mitochondrial function. We choose Lactobacillus Plantarum L168 as the chassis and pLDHLH673 as the vector. The gene circuit begins with the promoter ldhl, followed by the PncA gene, whose overexpression enables our bacteria to produce additional NAD+ for the host.

The DNA Level

$Lactobacillus Plantarum$ L168, a strain we obtained from our PI, Xingyin Liu's lab, has been proven to be able to alleviate social behavior deficits in animal models of autism. It carries the PncA gene itself and our aim is to amplify the gene, insert it into the vector, and transform the vector back to the L168 which enables the engineering bacteria to overexpress the PncA encoding protein. Thus, we designed primers, performed PCR, and ran gel to verify the successful amplification of the PncA gene.

Figure 1. Proof of the successful amplification of the PncA gene

As shown in the figure above, after PCR, the bands we obtained are near 600 bp, which was consistent with the target band (621 bp). This further indicated that we have successfully amplified the PncA gene.

We then cut bands out of the agarose gel and purify the DNA samples. Using pLDHLH673 as the vector, we constructed the pLDHLH673-PncA plasmid expression system, the main components of which is the PncA gene.

To prove our construction is successful, we performed enzyme digestion verification on the plasmid, adding restriction endonuclease NdeI and XbaI for double enzyme digestion, and ran the agarose gel to verify.

Figure 2. Proof of the successful construction of the pLDHLH673-PncA plasmid

As shown in the figure above, after enzyme digestion, the lower bands we obtained are near 600 bp, which was consistent with the target band (632 bp). The upper bands we obtained were higher than the marker and were near 5500bp, which was as we expected. This further indicated that have successfully constructed the pLDHLH673-PncA plasmid.

Although the construction was quite successful, we encountered difficulty in the electric transformation step. In fact, our original design used pMG36e as the vector. To amplify the plasmid, we need to transform it into Escherichia coli DH5a. However, every time we performed the transformation, extracted the plasmid, and ran the gel, we found no expected band on the gel. Considering that pMG36e was a low-copy plasmid, we switched our plasmid extraction kit from FastPure Plasmid Mini Kit (Vazyme) to TIANprep Midi Plasmid Kit (TIANGEN). However, the result was the same. Then, we changed the condition during the heat shock of the plasmid but the result showed no change. In this case, we had to change our vector to pLDHLMCS then to pLDHLH673, and finally constructed the target plasmid successfully. By this time, however, with few days left to finish the project, we were trapped in the electric transformation. After careful thinking and consulting our PI, we came to several explanations. First, the plasmid extraction kit we used was ineffective. In the lysing bacteria step, the color of the solution should have become clear after adding the reagent, but it remained unchanged. Second, the chassis bacteria we chose, Lactobacillus plantarum, are gram-positive bacteria that have different physiological characteristics from the gram-negative bacteria commonly used in plasmid extraction.

In fact, in the design stage, we discussed about using the chassis bacteria that are more common in synthetic biology. Due to its capacity to ameliorate social behavior deficiencies in animal models of autism, however, we elected to continue using it. Although we have not fully verified the functionality of this system, we believe that in the future, when given enough time, we can demonstrate its functionality and viability.

Function Verification

In our pre-experiments, we knocked out the pncA gene in Lactobacillus plantarum L168 to determine whether it would produce fewer NAD+.

Figure 3. NAD+ level and pncA knocked out

As shown in the figure above, the bacteria with pncA knocked out produce fewer NAD+ compared with that of Wildtype of LP-L168 supernatant through mass spectrometry assay, which implicated that our experimental design is reasonable. Despite the fact that this result does not directly demonstrate the NAD+-boosting ability of our engineering bacteria, it does provide indirect evidence.

Results of System 2

In system two, we aim at eliminating heavy metal ions such as Hg2+ and Pb2+ in the intestine of autistic children. We choose Escherichia coli 1917 as the chassis and pET-28a(+) as the vector. The gene circuit begins with two sensors which can be induced by Hg2+ and Pb2+ respectively and follows the Lpp-OmpA-MT gene, which encodes the protein that can anchor to the bacterial outer membrane and chelate the heavy metal ions. For biosafety, a kill switch that can be induced by arabinose is assembled downstream.

The DNA Level

We constructed the pET28a-MT plasmid, transformed it into Escherichia coli Nissle 1917, and extracted the plasmid to prove the success of our transformation. First, we performed enzyme digestion verification by adding restriction endonucleases SphI and Xhol for double enzyme digestion and running an agarose gel. Second, we performed the plasmid PCR verification, by designing primers to amplify the target segment and running an agarose gel.

Figure 4. Proof of the successful transformation (A) Enzyme digestion verification. Lanes 1: digestion sample of the extracted plasmid. (B) plasmid PCR verification. Lanes 1: PCR sample from the extracted plasmid.

As shown in the figure above, after enzyme digestion, the sample from the extracted plasmid ran out bands near 2000 bp, which was consistent with the target band (2399 bp). After plasmid PCR, the sample from the extracted plasmid also ran out bands near 2000 bp, which was consistent with the target band (2399 bp). This further indicated that our plasmid transformation was successful.

The Protein Level

To verify the expression of our functional genes Lpp-OmpA-MT, we induced the bacteria at different concentrations of lead ion salt solution. After lysing the bacteria and performing the centrifugation, we took the supernatant or precipitation at the first or third hour. Then we ran the SDS-PAGE to check MT protein level for all samples

Figure 5. SDS-PAGE of Lpp-OmpA-MT expression in different induction conditions Lane 1-6: 1h. Lane 7-12: 3h. Lane 1,3,5,7,11: Supernatant sample. Lane 2,4,6,8,10,12: Precipitation sample. Lane 13: negative control of the supernatant sample. Lane 14: negative control of the precipitation sample.

The pre-stained protein ladder is from 10 to 180 kDa. We induced the protein expression at different concentrations of lead ion salt solution (112 mg/L, 61mol/L, and 30.5 mol/L) and used double distilled water as the negative control. As shown in the figure above, the bands we obtained were nearly 40 kDa, which was almost twice the molecular weight we expected. After careful consultation with our PI and searching literature, we found it has been reported that MT easily form a dimer in vivo. Hence, we believed this indicated that Lpp-OmpA-MT proteins would combine and form a dimer in the bacteria, which confirmed the expression of our target protein.

Then, using ImageJ, we quantitatively analyzed the band brightness and obtained the following results: (ImageJ analysis: measure the brightness of corresponding lanes, minus the background brightness, and compare it with the marker to get the protein expression of each lane)

Figure 6. Pb$^{2+}$ concentration and the corresponding protein expression. (A) Supernatant. (B) Precipitation.

As shown in the figure above, the protein expression was greater in the precipitation than in the supernatant and was greater after three hours than after one, which was consistent with our expectations. Using the above data, we solved the model of inducer-protein expression. More information can be found on the Model page.

Function Verification

To determine whether our engineering bacteria can grow in the presence of heavy metals, we measured the growth curve. The most frequent model for describing the growth kinetics of bacteria is the logistic equation, which is frequently used to fit growth curves based on optical density (OD) in order to comprehend the dynamics of biological systems. Its differential expression is $\dfrac{1}{N}\dfrac{\partial N_t}{\partial t} = r(1-\dfrac{N_t}{K})$, by variation, we can obtain$N_t = \dfrac{K}{1 + \frac{K - N_0}{N_0} \times e^{-rt}}$

We use MATLAB to fit the data and obtain the figure and formulas shown below.

Figure 7. growth curve
Table 1. fit results of the growth curve
Variable Formula R$^2$

Wild Type in LB

$N_t = \dfrac{1.4240}{1 + 28.0422 \times e^{1.5709t}}$

0.9979

Wild Type in Pb(NO$_3$)$_2$ solution

$N_t = \dfrac{1.3153}{1 + 18.8537 \times e^{1.7823t}}$

0.9964

pET28a-MT in LB

$N_t = \dfrac{1.5713}{1 + 11.5395 \times e^{1.1985t}}$

0.9798

pET28a-MT in Pb(NO$_3$)$_2$ solution

$N_t = \dfrac{1.4928}{1 + 16.9268 \times e^{1.5685t}}$

0.9956

We used untransformed Escherichia coli Nissle 1917 (wild type) as the control. The concentration of Pb(NO$_3$)$_2$ solution is 112 mg/L. As shown in the figure above, the growth of our engineered bacteria will not be affected by heavy metals. Also, at the logarithmic phase, when the bacteria are most sensitive to changes in the environment, our engineering bacteria grew faster than the wild type, which indirectly proved that our engineering bacteria have the ability to eliminate heavy metals.

As to prove that our engineering bacteria can eliminate heavy metals directly, we put the engineering bacteria into the metal salt solution Pb(NO$_3$)$_2$ and cultured them for 2 hours. Then, we took the supernatant and measured the concentration of the metal salt solution using atomic absorption photometer.

Figure 8. proof of the heavy metal elimination ability

We used untransformed Escherichia coli Nissle 1917 (wild type) as the negative control. As shown in the figure above, the transformed bacteria eliminated more heavy metals than the wild type, and the difference was significant($P$$<$0.05). This indicated that our engineered bacteria were able to eliminate heavy metals. In the future, we will adjust the reaction time with the metal salt solution to further optimize the experiment condition. Kill Switch Verification

To validate the function of our arabinose-induced kill switch, we inoculated 3% fresh bacterial solution in the culture medium and added 0 mM and 10 mM arabinose respectively as the control and experimental groups. Every 1 hour, 2 hours, and 4 hours, we draw 100 ul of bacteria to the smear plate and cultured at 37 ℃ overnight. The colony counts are shown below.

Figure 9. colony count in the kill switch verification after dilution (the dilution multiple is 10$^9$) (A) 1h. (B) 2h. (C).3h
Figure 10. kill switch verification

As shown in the figure above, in one, two, and four hours, the number of colonies with arabinose induction are all lower than that without arabinose induction. This indicated the effectiveness of our kill switch.

Results of System 3

In system three, we aim at increasing the bioavailability of our engineering probiotics. We wrap them in alginate, coat them with whey protein and perform freeze-drying. Microencapsulation enables our engineering bacteria to go through the human digestive tract environment with less damage and fatalities. The experiments are conducted using the engineered probiotics in systems one and two, respectively. For convenience, our results are referred to by their chassis names.

Morphology

To observe the morphology of the microencapsulated bacteria, we took electron microscope photos. Since the morphology observation was considered to be the same for different chassis, we used Escherichia coli as an example.

Figure 11. Electron microscopy $(Escherichia coli)$

As shown in the figure above, the black structure is the bacteria and the greyish mist surrounding the bacteria is the coating of sodium alginate.

Then, we investigated the microencapsulation yield and effect of Lactobacillus Plantarum and Escherichia coli, respectively.

  • $Lactobacillus$ $Plantarum$
  • $Escherichia$ $coli$

Microencapsulation Yield

To determine whether sufficient bacteria remain after microencapsulation, we counted the number of our engineering bacteria before and after microencapsulation and calculated the microencapsulation yield.

Figure 12. Microencapsulation Yield of $Lactobacillus Plantarum.$ (A) Bacterial count before microencapsulation. (B) Bacterial count after microencapsulation (Dissolved in Sodium Citrate).
Table 2 Microencapsulation yield of $Lactobacillus plantarum&. (with dilution multiples calculated)
Num Bacterial count before microencapsulation $N_{L0}$ Bacterial count after microencapsulation (Dissolved in Sodium Citrate) $N_{L1}$

1

$7.2 \times 10^9$

$6.1 \times 10^{10}$

2

$4.7 \times 10^9$

$7.1 \times 10^{10}$

3

$3.4 \times 10^9$

$7.1 \times 10^{10}$

Mean

$7.1 \times 10^{10}$

$6.5 \times 10^{10}$

Microencapsulation yield of $Lactobacillus Plantarum$ $=\dfrac{lgN_{L1}}{lgN_{L0}} \times 100\% = 111.38\%$

The yield of microencapsulation was as expected, which indicated that sufficient engineering bacteria remained after microencapsulation. During freeze-drying, the engineering bacteria may continue to grow, which is normal and may lead to the yield exceeding 100%.

Microencapsulation Effect

To verify the effect of microencapsulation, we prepared simulated gastric juice and intestinal juice. Then we used them to process our engineering bacteria. The procedure was called simu-digest for convenience.

We defined survival rate as the ratio of the bacteria count after digestion to the bacteria count before digestion, which is

Survival rate = $\dfrac{\text{bacteria count after digestion}} {\text{ bacteria count before digestion}} \times 100\%$

By comparing the mean survival rates of Microencapsulated and Unmicroencapsulation bacteria, we can verify the effect of microencapsulation. Noted, the count of bacteria unmicroencapsulated and microencapsulated before digestion is identical to the count of bacteria before and after microencapsulation in the Microencapsulation Yield section, so we reused the data of them.

Figure 13. Microencapsulation Effect of $Lactobacillus plantarum$. (A) Unmicroencapsulated bacteria. (B) Microencapsulated bacteria.
Table 3 Microencapsulation Effect of $Lactobacillus Plantarum$.
Num. Unmicroencapsulated bacteria Microencapsulated bacteria

Before digestion

After digestion

Before digestion

After digestion

1

$7.2 \times 10^9$

$9.2 \times 10^8$

$6.1 \times 10^{10}$

$7.0 \times 10^{10}$

2

$4.7 \times 10^9$

$8.6 \times 10^8$

$6.3 \times 10^{10}$

$5.6 \times 10^{10}$

3

$3.4 \times 10^9$

$6.2 \times 10^8$

$7.1 \times 10^{10}$

$4.7 \times 10^{10}$

Mean

$5.1 \times 10^9$

$8.1 \times 10^8$

$6.5 \times 10^{10}$

$5.8 \times 10^{10}$

According to the data in the table above, we can know that the mean survival rate of the unmicroencapsulated bacteria was 15.69% whereas the mean survival rate of the microencapsulated bacteria was 89.23%. This indicated that the microencapsulation of Lactobacillus Plantarum is effective.

Microencapsulation Yield

To determine whether sufficient bacteria remain after microencapsulation, we counted the number of our engineering bacteria before and after microencapsulation and calculated the microencapsulation yield.

Figure 14. Microencapsulation Yield of $Escherichia coli$. (A) Bacterial count before microencapsulation. (B) Bacterial count after microencapsulation (Dissolved in Sodium Citrate).
Num Bacterial count before microencapsulation &N_{E0}& Bacterial count after microencapsulation (Dissolved in Sodium Citrate) &N_{E1}$

1

$5.7 \times 10^{11}$

$5.9 \times 10^{10}$

2

$1.1 \times 10^{12}$

$4.7 \times 10^{10}$

3

$1.0 \times 10^{12}$

$5.2 \times 10^{10}$

Mean

$8.9 \times 10^{11}$

$5.3 \times 10^{10}$

Microencapsulation yield of $Lactobacillus Plantarum$ $=\dfrac{lgN_{E1}}{lgN_{E0}} \times 100\% = 89.73\%$

The yield of microencapsulation was as expected, which indicated that sufficient engineering bacteria remained after microencapsulation.

Microencapsulation Effect

To verify the effect of microencapsulation, we prepared simulated gastric juice and intestinal juice. Then we used them to process our engineering bacteria. The procedure was called simu-digest for convenience.

We defined survival rate as the ratio of the bacteria count after digestion to the bacteria count before digestion, which is

Survival rate = $\dfrac{\text{bacteria count after digestion}} {\text{ bacteria count before digestion}} \times 100\%$

By comparing the mean survival rates of Microencapsulated and Unmicroencapsulation bacteria, we can verify the effect of microencapsulation. Noted, the count of bacteria unmicroencapsulated and microencapsulated before digestion is identical to the count of bacteria before and after microencapsulation in the Microencapsulation Yield section, so we reused the data of them.

Figure 15. Microencapsulation Effect of $Escherichia coli.$. (A) Unmicroencapsulated bacteria. (B) Microencapsulated bacteria.
Table 2 Microencapsulation yield of $Lactobacillus plantarum&. (with dilution multiples calculated)
Num. Unmicroencapsulated bacteria Microencapsulated bacteria

Before digestion

After digestion

Before digestion

After digestion

1

$5.7 \times 10^{11}$

$6.0 \times 10^{10}$

$5.2 \times 10^{10}$

$3.1 \times 10^{10}$

2

$1.1 \times 10^{12}$

$4.2 \times 10^{10}$

$5.9 \times 10^{10}$

$3.2 \times 10^{10}$

3

$1.0 \times 10^{12}$

$1.5 \times 10^{11}$

$4.7 \times 10^{10}$

$3.2 \times 10^{10}$

Mean

$8.9 \times 10^{11}$

$8.4 \times 10^{10}$

$5.27 \times 10^{10}$

$3.17 \times 10^{10}$

According to the data in the table above, we can know that the mean survival rate of the unmicroencapsulated bacteria was 9.44% whereas the mean survival rate of the microencapsulated bacteria was 60.12%. This indicated that the microencapsulation of Lactobacillus Plantarum is effective.

Results of System 4

In system four, we aim at sensing the elevated concentration of lactate to detect the potential mitochondrial dysfunction in autistic children. We choose Escherichia coli DH5α as the chassis and pSB4K5 as the vector. The gene circuit begins with a sensor that can be induced by lactate and follows the lacZ gene, which encodes the protein β-Galactosidases that can have a color reaction with X-gal.

The DNA Level

We constructed the ALPaGA-lacZ plasmid and transformed it into Escherichia coli DH5α. To prove the success of our transformation, we designed primers, performed plasmid PCR and colony PCR to amplify the LldR gene, and ran the agarose gal.

Figure 16. Proof of the successful transformation. (A) plasmid PCR. Lanes 1-2: PCR sample form transformed DH5α. Lanes 3-4: PCR sample form untransformed DH5α. (B) colony PCR. Lanes 1-4: PCR sample form transformed DH5α. Lanes 5-8: PCR sample form untransformed DH5α.

As shown in the figure above, after plasmid PCR and colony PCR, the sample from the transformed DH5α both ran out bands near 750 bp, which was consistent with the target band (777 bp). This further indicated that our plasmid transformation was successful.

For our reporter gene lacZ, in order to determine whether our chassis bacteria Escherichia coli DH5α themselves carry this gene, which may interfere with the results of our functional analysis. We designed primers, performed colony PCR to amplify the lacZ gene, and ran agarose gel.

Figure 17. PCR amplification of lacZ reporter gene. Lanes 1-2: PCR sample form transformed DH5α. Lanes 3-4: PCR sample form untransformed DH5α.

As shown in the figure above, after colony PCR, the sample from the transformed DH5α ran out bands near 3000 bp, which was consistent with the target band (3075bp). The sample from the untransformed DH5α can also run out bands near 3000 bp but with much weaker light. This indicated that our chassis bacteria Escherichia coli DH5α themselves carried the lacZ gene and could express them in micro amounts.

The protein Level

To verify the expression of our reporter genes lacZ, we induced the bacteria at different lactate concentrations. After lysing the bacteria and performing the centrifugation, we took the supernatant and ran SDS-PAGE.

Figure 18. SDS-PAGE of β-galactosidases expression in different induction conditions.

The pre-stained protein ladder is from 10 to 180 kDa. We induced the protein expression at different lactate concentrations (1 mol/L, 10-3 mol/L, and 10-6 mol/L) and used double distilled water as the negative control. As shown in the figure above, the bands we obtained were close to the target band (113 kDa). This confirmed the expression of our target protein.

Then, using ImageJ, we quantitatively analyzed the band brightness and obtained the following results: (ImageJ analysis: measure the brightness of corresponding lanes, minus the background brightness, and compare it with the marker to get the protein expression of each lane)

Figure 19. Lactate concentration and the corresponding protein expression

As shown in the figure above, as the lactate concentration increased, the protein expression also increased gradually, which was consistent with our expectations. Using the above data, we solved the model of inducer-protein expression. More information can be found on the Model page.

Function Verification

The function of this system is to sense the elevated concentration of lactate and then detect potential mitochondrial dysfunction in autistic children. To verify this, we used filter paper, a common supply in the laboratory, to perform the experiment. First, we dropped the bacteria solution onto the filter paper. After 30 minutes, our engineering bacteria were firmly attached to the filter paper. Then we added lactate and x-gal onto the filter paper and observed the color changes.

Figure 20. color of the filter paper under different lactate concentrations

We performed the function analysis at different lactate concentrations (1 mol/L, 10$^{-3}$ mol/L, and 10$^{-6}$ mol/L) and used double distilled water as the negative control.

As shown in the figure above, without the induction of lactate, the lacZ protein produced by the chassis bacteria was insufficient to color the filter paper. When adding 10$^{-6}$-1 mol/L lactate, the filter paper turned blue, and as the concentration of lactate increased, the blue became darker. This confirmed that our lactate test strip was functional.

However, the color change was not very significant to the naked eye, which may be related to our plasmid being a low copy, the material of filter paper, the lower activity of bacteria after fixing on the test paper, and the unsuitable reaction environment.