Team:Valencia Biocampus/Results

From 2014.igem.org

Revision as of 15:05, 17 October 2014 by TonnyESP (Talk | contribs)

Results

Stability

Escherichia coli is the lab rat of the bacterial world. We have performed the most complete characterization of the stability of E. coli in order to both confirm its optimality as a chassis for Synthetic Biology, and to determine the behavior of Biobricks under sub-optimal environmental conditions. We have subjected wild type strains to a range of stressful conditions (see just below "the limits of E. coli") and measured how our Biobricks behaved.

List of protocols here

1 . The Limits of E.coli

Temperature

E.coli is a fully studied and characterized bacterium. The optimal growth of this species is 37-38ºC with continuous agitation (200 r.p.m.) but out of this range the growth slows down. First of all, we measured the growth of E.coli at temperatures ranging from 30 to 50ºC. To achieve our goal, we used a fully original technique: instead of an incubator, we used a thermal cycler programmed to display a gradient of temperatures.

thermal cycler Drawing

This proved a very easy and reproductible way of very precisely incubating at increasing temperatures small (PCR tubes) bacterial cultures.

The results obtained are shown below:

Figure 1
Figure 1: Growth of E. coli strains DH5α (left) and XL1-Blue (right) under different temperatures (30ºC-50ºC) set in a thermal cycler. The two strains display clearly different temperature profiles.

Conclusion

Is E. coli stable at a wide range of temperatures? Yes it is! We have detected some variations between the two strains. For example, DH5α resists low temperatures better than XL1-Blue, but both grow at a very large range of non-optimal temperatures.

Material fatigue (1)

Our modeling sub-team suggested an original idea. Is there something similar to material fatigue in biology? For example, how would an E.coli culture behave if subjected to a temperature fluctuation? We carried out some experiments with our thermal cycler by cycling the temperature of the medium every minute (37-41ºC) for 3h and maintaining controls at fixed temperatures (typically 37 and 41ºC). The results of these experiments are shown in Figure 2.

Figure 2
Figure 2: XL1-Blue (left) and DH5α subjected to cycling (37-41ºC) temperature conditions. Compared to the controls, the former strain did not exhibit any significant change in the growth rate, as deduced by the OD600 after 3 hours (DH5α) and 15 hours (XL1-Blue). By contrast, DH5α significantly increased its growth compared to both controls. The statistical parameters are cited in the text.

The ANOVA analysis performed with the 16 replica of each strain revealed no significant differences in the growth of Xl1-Blue (F=0.870, p=0.433). However, the differences in the growth of DH5α did have statistical significance (F=63.95, p<0.001).

Conclusion

Is E. coli sensitive to fatigue of materials in the form of cycling conditions? Yes! at least DH5α, but in a positive way. Cycling temperatures resulted in an emergent property: increased growth. In order to check whether this effect was due to a mechanical resuspension (convection) linked to the cycling conditions, we performed a specific experiment with dead, antibiotic-treated E. coli cells (Kanamicine and Cloranfenicol) subjected to extreme (30 to 50ºC) thermal fluctuations. The following figure shows that convection forces might not be the reason of such increased growth (no significant differences found in a t-Student test for 2 groups and 16 replicas; t=1.371, p=0.180). A biological explanation may be hypothesized: cycling conditions allow a range of optimal enzymatic activities to be met.

Figure 3
Figure 3: O.D.600 of antibiotic-killed E. coli suspensions of strain XL1-Blue incubated under cycling (30ºC-50ºC) and average (40ºC) conditions. No significant differences suggest that convection forces are not the main key players in the increased growth of E. coli incubated under thermally cycling condition. t-Student test 95% showed no significant differences, as explained in the text.

UV Radiation (1)

UV radiation generally carries out its biological effect by damaging nucleic acids. Depending on the radiation intensity and exposure time the damage caused to the genome is reversible; cells can repair their DNA and continue to grow, although they might carry point mutations. Since UV studies are scarce, our team decided to assay the maximum resistance of E. coli strains to UV radiation of a fixed intensity (340 µW/cm2). The experiment consisted of subjecting spot plated E.coli strains to timed pulses of UV light and comparing the output to a control set of bacteria which received no radiation at all.

Figure 4
Figure 4: Characterization of XL1-Blue (left) and DH5α (right) subjected to UV light (340 µW/cm2). Both strains are rather sensitive, but DH5α resists UV light better than XL1-Blue.

Therefore, E. coli proved very sensitive to UV light We also studied how sensitive these strains are when transformed with our Biobricks. As shown in Fig 5, Biobrick-expressing strains were more sensitive to UV light compared to plasmid-only strains. However, Biobrick7 provided some resistance to radiation in DH5α compared with the controls.

Figure 5
Figure 5: Resistance of DH5α (left) and XL1-Blue (right) to UV pulses (340 µW/cm2). X axis values are given in seconds. Both strains are very sensitive to UV light with a maximum resistance of around 6 seconds.

pH & Salinity (1)

We also characterized E. coli resistance to a range of extreme pH values and salt concentrations. As shown in Table 1, none of the extreme conditions tested resulted in a total inhibition of E. coli’s growth.

Table 1.

pH

5

6

7

8

9

XL1-Blue

+

++

+++

++

+

DH5α

++

+++

+++

++

+

Salt

Without salt

LB

LB+1% salt

LB+2% salt

LB+4%salt

XL1-Blue

+++

+++

++

+

+

DH5α

+++

++

+

+

+

Vacuum (1)

We designed a simple experiment to subject E.coli to Mars-like vacuum. It consisted of streaking Petri dishes and incubating them in deep vacuum during 48 hours. We also maintained control plates, grown at atmospheric pressure. After this incubation period, we measured the diameter of three colonies of each Petri dish. The results, indicating only a moderate decrease in colony size after two days of vacuum incubation are shown below:

Figure 6
Figure 6: Mean Colony Diameter of XL1 (left) and DH5alpha (right) grown under vacuum and atmospheric pressure. Growth differences were statistically significant according to a t-Student test for two groups and three replicas (XL1: t=4.24, p=0.013<0.05 ; DH5: t=9.98, p<0.001).

BL21. Chromosomal T7 polymerase gene is unstable!

BL21 is an Escherichia coli strain that bears the T7 phage polymerase gene into its genome. Therefore, it is used to express proteins that are controlled by T7 phage promoters, as is the case of three of our Biobricks (Biobricks 8, 9 & 10). However, when we transformed the strain with the Biobricks (blue and yellow chromoproteins), we could not measure any signal. After many tests, we hypothesized that the T7 phage polymerase gene was unstable and did not resist stressful situations such as competence or transformation protocols.

We performed a PCR screening to find out whether the T7 gene was present. We used colonies of BL21 transformed with Biobrick 8 (blue chromoprotein) and grew them on LBA medium. We performed a colony-PCR with a couple of T7 specific oligos designed by ourselves targeting a 1.5 kb region of the T7 polymerase coding sequence.

Figure 7
Figure 7: 0.8% agarose gel electrophoresis. We saw that the bacteria forming the colonies did not have the gene that codes for the polymerase.
We then decided to check whether our glycerol stock of bacteria contained the gene. Maybe during the process of competence the culture was contaminated. To our surprise, the gene was present in every single colony we analyzed. Therefore, we decided to do a final colony PCR using the bacteria from the glycerol stock as positive controls, and using colonies from the non-transformed competent stock.

Figure 8
Figure 8: 0,8% agarose gel electrophoresis. In the BL21-glicerol stock the gene was present. Nonetheless, almost all BL21-competent cells lack T7 phage polymerase gene.

Our results reveal that BL21 loses the T7 polymerase gene at some point during the process of making the cells competent, as only three out of the ten colonies analyzed did result in (a rather faint) amplification of the gene.

2 . Behavior of different Biobricks in sub-optimal conditions

We show below a summary of our Biobrick characterization program, in which we have subjected selected Biobrick-expressing strains to a range of environmental stresses. Our results clearly show that whereas UV light, pH variation and salt content do not significantly affect Biobrick expression, sub-optimal temperatures can in some cases decrease the Biobrick output. Anaerobic, vacuum culture resulted in a dramatic decrease of Biobrick 2 (RFP) activity.

Temperature Stress

37ºC and 200 rpm are the standard growth conditions of E.coli. Varying these conditions and measuring the growth after three hours, (O.D.600) as well as fluorescence (Biobrick output), was used to test the effect on Biobrick expression (all data are the average of at least three measurements; fluorescence values are normalized to OD600). We also built a model of the temperature stability of the constructions used in this experiments, as you can see in the graphs.

Figure 9A
Figure 9A: The strain XL1-Blue transformed with a Biobrick which encodes GFP (Biobrick 1, Bb1) was grown at different temperatures in a thermocycler for three hours, the cells were collected and O.D600 was measured. Fluorescence was also determined and compared to the non-transformed/wildtype or control strain. As the plot shows, Biobrick expression is affected by growth at different temperatures, fluorescence (normalized to O.D600) is higher around E.coli optimal growth temperature (37-41 ºC); over or below this rather short temperature range fluorescence decreases. The green line represents the theoretical prediction of our model. Click here to see our model.
Figure 9B
Figure 9B: The strain XL1-Blue transformed with a Biobrick which encodes GFP (Biobrick 3, Bb3) placed under the control of a weak promoter. The micro-cultures were grown at different temperatures in a thermo cycler for three hours, the cells were collected and O.D600 was measured. Fluorescence was also determined and compared to the non-transformed/wildtype or control strain. Fluorescence values in the different temperatures tested are very similar to the control and only slightly higher around optimal temperatures. This indicates a low expression of Bb3 in all cases. The green line represents the theoretical prediction of our model. Click here to see our model.

Material fatigue (2)

As described in the first parts of this section “The limits of E. coli”, the strain DH5α grows better under temperature cycling conditions "thermal fatigue" compared to the optimal growth temperature (37ºC), but the same effect was not observed for the strain XL1-Blue. Hence, we decided to check what was the effect of Biobrick expression in XL1-Blue (Fig 10).

Figure 10
Figure 10: Biobrick 1 expression under different temperatures, including temperature fluctuation (from 37ºC to 41ºC, every minute each) in a thermocycler. Controls were incubated at 37 and 41ºC. After an incubation period of 3 hours, O.D.600 alongside fluorescence was measured. The results show that the expression of the Biobrick is not affected by cycling temperatures. ANOVA test 95% shows no significant differences, as explained in the text.

No significant differences were found among treatments according to an ANOVA test for three groups and 16 replicas (F=1.682, p=0.200) at a 95% level of significance.

pH & Salinity (2)

Large proteins, including enzymes, are affected by pH. Usually, the catalytic properties of the enzymes are lost and metabolism is halted, which eventually leads to death. The pH range determined for E.coli optimal growth is 6 to 7, but does Biobrick expression affect this? The results of the experiment in which the pH of the medium was varied are displayed below. The results show that, as expected, Bb1 displays a higher fluorescence and is stably active regardless the pH value of the medium.

DH5α pH

The following figure shows our experimental results and the prediction of our theoretical model:

Figure 11
Figure 11: The plot displayed above shows the variation in fluorescence of the transformed DH5α cells grown in LB with pH ranging from 5 to 8. Cells were transformed with GFP constructs, one of which is controlled by a strong constitutive promoter (BB1) and the other, BB3, controlled by a weak housekeeping promoter. The output from these is compared to a cell transformed with an empty plasmid (Ø Ampicillin). Click here to see our model. .

The differences found among the range of pH values were significant for BB1 (F=2.654, p=0.075) but not significant for BB3 (F=7.6, p=0.001) at a 95% level of confidence.

DH5α salinity

Regarding salinity, a major environmental factor, we also tested the stability of our Biobricks by measuring their output in increasingly salty media (Fig. 12). Again, there were no significant differences associated with salt concentration in the case of BB1 (ANOVA values: F=0.832, p=0.491) but they were significnat for BB3 (ANOVA values: F=3.149, p=0.045).

Figure 12
Figure 12: Using the same approach as the previous experiment, the same strains (DH5 alpha) were grown in mediums with different salt concentrations (LB without salt, regular LB, and LB supplemented with 1% and 2% NaCl). Fluorescence of the GFP Biobricks was measured and results normalized to the O.D.600. As before, Bb1 was significantly expressed compared to the control at values relatively independent of the salt concentration. ANOVA 95% test equal mean values hypothesis rejected, as explained in the text.

UV Radiation (2)

We display below the results of Biobrick behavior upon irradiation by measuring fluorescence of the surviving cells. We found out that the fluorescence was not affected by UV exposure (Figures 13 and 14), since the t-Student test for two groups and three replica showed no significant differences (BB1: t=1.011, p=0.418; BB3: t=1.240, p=0.341)

GFP

Figure 13
Figure 13: Biobrick 1 and 3 behaviors before and after a 2 second UV pulse. Left, DH5α; Right, XL1-Blue. The activity of both Biobricks is not affected by previous exposure to UV. In both cases, t-Student test 95% shows no significant differences, as explained in the text.

RFP

Similar results were obtained with BB2, finding, again, no significant differences under the different treatments (DH5α: t=-0.852, p=0.499; XL1-Blue: t=0.428, p=0.754).

Figure 14
Figure 14: Biobrick 2 behavior before and after a 2 second UV pulse. Left, DH5α; Right, XL1-Blue. The activity this Biobrick is not affected by previous exposure to UV. Unexpectedly Bb2 does not show any significant fluorescence even in the non-UV-treated samples. In both cases, t-Student test 95% shows no significant differences.

Vacuum (2)

Similarly to the experiment carried out in the section “Limits of E.coli” with vacuum, we plated the strains and incubated them in deep vacuum for 48 hours. After this incubation period, cells were recollected and O.D.600 as well as (normalized) fluorescence, determined. Surprisingly, Biobrick 2 proved extremely sensitive to vacuum in both strains, since normalized fluorescence dramatically decreased compared to the normal pressure-incubated cultures (Figure 15).

RFP

Figure 15
Figure 15: Effect of vacuum on the expression of the Biobrick 2 (RFP) in two strains of E.coli. In both cases vacuum-cultured cells display a much lesser normalized fluorescence that the normal pressure-incubated ones (all assays were conducted at 37ºC).

Standardization

Biobrick parts are supposed to be standard, this is, they are expected to behave in a similar and predictable way in different host organisms or “chassis”. In this part of the project we used six different Escherichia coli strains commonly used in molecular biology labs (DH5α, XL1-Blue, JM109, DH10B, HB101, BL21(DE3)) to test whether the behavior of our Biobrick parts was homogeneous among them (Bb1 to Bb6, BBa_K1468000-BBa_K1468005 in the Registry of Standard Biological Parts). Surprisingly, only one of the constructions we tested displayed a relatively standard output among the different strains (BB5).

General protocol

  • Strike the desired strains on LB plates or LB plates supplemented with the appropriate antibiotic (depending on the bacteria we are working with: the wild type, sensitive to the antibiotic, or the transformed ones, resistant to the selective agent). Allow to grow by incubating at 37°C at least 24 hours.
  • Gather a large amount of bacteria with an inoculation loop and transfer to 1 mL of LB or LB supplemented with the antibiotic. Incubate during 20 minutes at 37°C and 200 rpm.
  • Inoculate 100 µL of bacterial suspension in a tube containing 3 mL of LB with antibiotic. Perform four biological replica of each culture. The expression of Bb4 was induced with a heat shock at 42 °C for 2 min. The expression of Bb5 and Bb6 was induced by adding doxycycline and IPTG, respectively, to the medium.
  • Incubate until OD595 is between 0.15-0.4 (1 hour aprox.).
  • Measure OD595.
  • In the case of Biobrick parts 4, 5 and 6, sonicate an aliquot of the culture (it allows us working with the pretein extract, avoiding the interference of cell density in the measurement of the blue, product of the XGal reaction), eliminate cell debris by centrifugation (2 min, max. rpm), add X-Gal, incubate 6 min at RT, and stop the reaction by adding Na2CO3.
  • Measure either GFP fluorescence (Bb1 and Bb3; λex=493nm; λem=505nm), RFP fluorescence (Bb2; λex=576 nm; λem=592 nm), or absorbance at 630 nm (Bb4, Bb5, and Bb6). Fluorescence measurements were performed in a Hitachi F7000 plate reader, whereas a spectrophotometer was used for absorbance measurements.
  • With Excel, normalize fluorescence by cell density, and calculate the average and standard deviation taking into account the four biological replicas. Represent the data.

Results

The behavior of all the Biobrick parts we tested is shown in Figure 1. From our results, it is evident that the expression of all but one of the constructions, Bb5, clearly differed depending on the host strain of choice. Even for Bb5, using an ANOVA test for six groups and 4 replicas, one can reject the null hypothesis of equal mean values (F=4.253, p=0.010) which means that the expresion of the Biobrick is indeed different. Each strain proved to be more “efficient” at hosting a different Biobrick part. For instance, Bb1 and Bb5 are highly expressed in strain JM109, whereas the expression of the other biobricks was low or nearly undetectable. Note that Bb3, containing the weak promoter J23110, could not be detected in four out of the six strains. It is also remarkable the fact that in 1A and in 1B (Bb1 and Bb2 are under the control of the same promoter) the pattern obtained is almost equal among strains except in JM109.

Figure 1
Figure 1. Behavior of Biobrick parts in six different E. coli strains. Fluorescence intensity of Bb1, Bb2, and Bb3 (A, B, and C, respectively) was corrected by the fluorescence of cells containing an empty plasmid. Absorbance values of Bb4, Bb5, and Bb6 (D, E, and F, respectively) was corrected by the absorbance displayed by non-transformed cells (bacteria containing an empty plasmid where not used because the white-blue screening sytem present in the vector).

Standarization modelling, click here to see

Discussion and conclusions

We demonstrate here a difference in the output of the Biobricks depending on the host strain. If we assume that they are in the same physiological state this results could be related to epistasis, since epistasis describes how gene interactions can affect phenotypes (Miko, I. Epistasis: Gene interaction and phenotype effects. Nature Education (2008) 1(1):197). In other words, the behaviour of a particular gene strongly depends on the presence of others, on the genetic background of each Escherichia coli strain.

The different behavior can be a consequence of the synthesis rate, which also depends on the genetic background. The measures were made in the exponential phase of growing, where the bacteria are constantly dividing, in such a way that the protein synthesized is also being divided. Maybe, if the same experiment was carried out when the bacterial cells had reached the stationary phase, the output obtained could be more homogeneous. We are currently testing that.

In any case, with the results we obtained, we have to conclude that Biobrick parts do not behave in a standard way even whitin strains of the same bacterial species.

Orthogonality

In Synthetic Biology, two constructions can be considered orthogonal when they only interact at specific and predictable interphases, and do not disturb each other. We have studied this desirable feature of Biobrick parts by combining two constructions in the same cell and comparing this output with the one produced by simple transformations. To do this, we used both standard fluorimetry assays and flow cytometry, which allowed to get information at the cell level.Then, we wanted to go one step further: what is the effect of a simple transformation (a plasmid with a Biobrick part) into the cell architecture? To get some insights about this, we performed a proteomic analysis in which the whole proteome of an E. coli strain transformed with a Biobrick part was compared to that of the non-transformed, control strain.

1. Standard fluorimetry assays

We performed experiments with two different combinations of Biobrick parts containing a fluorescent protein under the control of a promoter sequence from Anderson’s promoter collection: Bb1 (GFP with strong promoter J23104) + Bb2 (RFP with strong promoter J23104) and Bb2 + Bb3 (GFP with less strong promoter J23110). In all cases, several controls were used: wild type (non-transformed) cells, simple transformations with either the biobrick or an empty plasmid, and co-transformations of the biobricks and the empty plasmids. It has to be noted that there was a high background emission of red fluorescence by the cells in all cases. Our results are shown in Figures 1 and 2.

General protocol

  • Strike the desired strains on LB plates with the appropriate antibiotic. Allow to grow by incubating at 37°C at least 24 hours.
  • Gather a large amount of bacteria with an inoculation loop and transfer to 1 mL of LB. Incubate during 20 minutes at room temperature.
  • Inoculate 100 µL of bacterial suspension in a tube containing 3 mL of LB with antibiotic. Perform three biological replicas of each combination of Biobrick parts.
  • Incubate until OD600 is between 0.1-0.3.
  • Measure OD600.
  • Measure GFP fluorescence (exc.= 493nm; emis.=505nm) and RFP fluorescence (exc.= 576 nm; emis.= 592 nm) in a FP6200 spectrofluorimeter (Jasco, Easton, MD) fluorometer using standard plastic cuvettes.
  • With excel, normalize fluorescence by cell density, calculate the average and standard deviation taking into account the three biological replicas, and represent the data.

Results

Figure 1
Figure 1. Fluorescence emission in XL1-Blue Escherichia coli strains co-expressing Biobrick parts 1, 2, and 3 reveals asymmetry of behaviors (lack of orthogonality). Left, cells containing both Biobrick parts 1 and 2, Biobrick part 1 and an empty kanamycin plasmid (e.k.p), Biobrick part 2 and an empty ampicillin plasmid (e.a.p.), Biobrick part 1, and Biobrick part 2. Right: Cells containing both Biobrick parts 2 and 3, Biobrick part 2 and an empty kanamycin plasmid (e.k.p), Biobrick part 3 and an empty ampicillin plasmid (e.a.p.), Biobrick part 2, and Biobrick part 3. XL1-Blue E. coli cells containing no plasmids or empty plasmids were used as controls. Both green (505 nm) and red (592 nm) fluorescence intensity (FI) were measured for all samples and normalized by the optical density at 600 nm (OD600).

Figure 2
Figure 2. Fluorescence emission by DH5α E. coli strains co-expressing Biobrick parts 1 and 2 also reveals non orthogonal expression of Biobrick parts. Cells containing both Biobrick parts 1 and 2, Biobrick part 1 and an empty kanamycin plasmid (e.k.p), Biobrick part 2 and an empty ampicillin plasmid (e.a.p.), Biobrick part 1, and Biobrick part 2 were analyzed. XL1-Blue E. coli cells containing no plasmids or empty plasmids were used as controls. Both green (505 nm) and red (592 nm) fluorescence intensity (FI) were measured for all samples and normalized by the optical density at 600 nm (OD600).

When in co-transformation, Bb1 and Bb2 showed a clearly asymmetric expression . Bb1 was highly expressed in comparison to Bb2 (2-fold in strain XL1-Blue and 3-fold in DH5αα). A similar phenomenon was observed when either Bb1 or Bb2 were cotransformed with empty plasmids. Again, less fluorescence was detected in comparison to the simple transformants. Surprisingly, this fluorescence decay was strikingly different between both E. coli strains: Bb1 strongly lowered its expression in DH5α, whereas Bb2 showed a higher decay in XL1-Blue cells.

However, the result was completely different when Bb2 and Bb3 were present in the same cell. In the XL1-Blue strain, the fluorescence intensity of Bb2 was as high as that of the single transformant, whereas the fluorescence of Bb3 –containing a less strong promoter- dropped dramatically. Fluorescence intensity of Bb3 was not detectable in strain DH5α under our experimental conditions.

With these results with co-transformants, we have to conclude that Biobrick parts, even if controlled by the same promoter and cloned into the same plasmid, are not expressed at similar rates. But, is this lack of orthogonality stable and reproducible? In other words, is the asymmetric output always the same or does it fluctuate in time or among cells? We were able to answer this question by using flow cytometry. Check our results below!

2. Flow cytometry

To study cell-level fluorescence we used flow cytometry. To prepare the samples, XL1-Blue cultures (wild type, simple transformants for Bb1 and Bb2, and cotransformants with Bb1+Bb2) were grown in LB with the appropriate antibiotic, harvested and resuspended in 1mL PBS. From now on, we will refer to Bb1 as GFP, and Bb2 as RFP.

Orthogonality at the cell level

Our first result, shown in Figure 3, is the confirmation of non-orthogonal behavior as deduced by a higher expression of GFP in comparison to RFP in cotransformed XL1-Blue cultures. There is four to five times more expression of GFP in comparison to RFP (Fig 3D). These results are similar to those we previously obtained by fluorimetry (with younger cultures, 5 h), in which an imbalanced GFP:RFP rate of about 2:1 was detected. It has to be stressed that, unexpectedly, only a fraction (around 10%) of the cells exhibited activity of both reporter proteins, indicating that the fluorescence of a culture as it is usually measured (fluorimetry) represents in fact and average of a very diverse pool of individual outputs, including mono-active cells. We plan to repeat the assay using the reverse conditions: GFP in a Kan resistant plasmid, and RFP in an Amp resistant plasmid.

Figure 3
Figure 3. Dot plots of XL1-Blue reveal a non-orthogonal behavior. Dot plots obtained by flow cytometry. A, B, C & D correspond to wild type, GFP transformant, RFP transformant and cotransformed cells, respectively. FITC-A, or X axis, indicates green fluorescence (505 nm), whereas PE-A, or Y axis, indicates red fluorescence (592 nm).

Orthogonality changes in time

Figure 4 shows a clear correlation of fluorescense and time in mono-transformants. Fluorescence is very low for all samples at 5h (approximately 30% of cells transformed with Bb1 express GFP but only 1% of cells transformed with Bb2 express RFP); moderate at 10h for GFP mono-transformants (41% and 1.5% respectively); and maximum after 20h, when GFP is expressed by 80% of the cells, and RFP is expressed by approximately 45% of the cells. Regarding co-transformants, fluorescence also increased with time but even after 20h only a fraction (10%) of doubly fluorescent cells were detected.
Figure 4
Figure 4. XL1-Blue dot plots indicating an increase of fluorescent protein expression with time. Samples of the four variations (wild type, GFP transformants, RFP transformants, and cotransformed cells) are taken at 5, 10 and 20 hours. FITC-A, or X axis, represents green fluorescence (505 nm). PE-A, or Y axis, represents red fluorescence (592 nm).

These results indicate that higher fluorescence levels are reached when cells are reaching stationary phase and that even long incubation times fail to yield a majoritary population of co-transformed cells with dual GFP and RFP activity.

Fluorescense rates within cultures: a very stable imbalance

Our last assay aimed at answering this paradoxical question: is the lack of orthogonality a stable trait? In other words, is the GFP:RFP imbalance always the same? The answer, according to Figure 5 is clear: Yes, independent cultures exhibit a striking similarity in their flow cytometry plots. Almost half of the cell population only displays GFP activity, whereas only a relatively small fraction is either displaying RFP fluorescence or both.
Figure 5
Figure 5. Independent cultures of cotransformed cells indicate the lack of orthogonality is stable. Analysis by flow cytometry of XL1-Blue cotransformed cells. One dot plot is shown for each independent culture. FITC-A indicates green fluorescence (505 nm). PE-A indicates red fluorescence (592 nm). Below: diagrams representing the percentage of cells that express no fluorescence, RFP, GFP, or both.

3. Proteomics

Our experiments involving E. coli cells co-expressing two different Biobrick parts revealed a lack of orthogonality between them. But what is going on with all the other parts naturally present in the chassis? Does the expression of a Biobrick part interact or interfere with any of the genes expressed by the host? In order to test this hypothesis, we performed a proteomic analysis on the simplest scenario: we compared the proteome of the DH5α E. coli strain expressing Biobrick part 1 (consisting of a green fluorescent protein under the control of a constitutive promoter) with that of the same strain carrying the empty cloning plasmid and the non-transformed wild-type strain.

General protocol

  • Set up three independent 5-mL cultures of the strains named before as we explained above for the standard fluorimetry assays.
  • When an OD value between 0.2 and 0.4 is reached, pellet the cells at 4ºC and then wash them with sterile ice-cooled PBS buffer.
  • Resuspend the cells in 1 mL of ice-cooled sterile PBS buffer.
After this, we kept the tubes on ice and brought them to the Proteomics lab of the University of Valencia. The proteomic analysis was carried out in the SCSIE_university of Valencia Proteomics Unit , a member of ISCIII ProteoRed Proteomics Platform. Once there, proteins were isolated and quantified with the supervision of the technicians. Finally, proteins were digested and labeled by using the iTRAQ technology, which allowed us to work with all the samples in the mass spectrometer at a time and then assign each identification to a particular sample.

Results

A general view on the protein profile of each sample was firstly observed on an SDS-PAGE gel (Figure 6). Strong differences were detected, as expected, in the bands corresponding to the GFP and the β-lactamase encoded by the ampicillin-resistance gene of the pUC57 plasmid used for cloning.

Figure 6
Figure 6. SDS-PAGE gel showing the protein profile of the samples used for proteomics analysis. WT: wild-type E. coli strain DH5α (2 replicates); Bb1: Biobrick part 1-expressing strain (3 replicates); ᴓ: strain carrying the empty cloning plasmid. Red and green arrows indicate β-lactamase and GFP, respectively.

An average of 1500 proteins per sample were identified with the mass spectrometer. Among them, 472 proteins could be properly quantified. Figure 7 shows the variation of a subset of proteins among samples. Notice that chaperone proteins are specially over-represented in the strain expressing Bb1, whereas proteins involved in response to oxidative stress are particularly abundant in the strains carrying the empty plasmid.

Figure 7
Figure 7. Representation of a subset of differentially expressed proteins among samples. Red and green colors indicate strong over-representation and under-representation, respectively. Blue arrows indicating β-lactamase and GFP; orange arrows indicating chaperone proteins; purple arrows indicating oxidative stress-related proteins.

But this apparent variation in the proteome changes when applying some statistics. A relative comparison of the average expression level of each one of the quantified proteins in the three groups of samples (WT strain, Bb1-expressing strain, and empty plasmid-carrying strain) is shown in Figure 8. Those proteins displaying the highest alterations in their expression levels are marked with numbers and identified in the table on the right.

Figure 8
Figure 8. Differential expression level of the proteins quantified by iTRAQ mass spectrometry. Each of the 472 quantified proteins has been represented for each group of samples: WT (blue line), Bb1-expressing strain (green line), and e.a.p. (empty Amp plasmid)-carrying strain (red line). Proteins on the X axis ordered in decreasing confidence level.

As deduced in Figure 7, a very low number of proteins strongly changed their expression in the strains carrying Bb1 or the empty plasmid. Several statistical tests (FDR-corrected ANOVA, T-test, p-value-weighted analysis) were performed with these proteins, and only three of them proved to be significant in all the tests. These proteins were the GFP (only present in Bb1 strain), the β-lactamase (present in both Bb1 and ᴓ strains), and the chaperone protein DnaJ, which was particularly abundant in the strain expressing Bb1. The high expression of chaperone proteins was expected, since many molecules of GFP need to be properly folded in order to avoid the formation of inclusion bodies and cellular toxicity. From this experiment, we have to conclude that the expression of a simple Biobrick part such as a fluorescent protein in an E. coli cell does not alter the proteomic architecture of the chassis in a significant way. Therefore, Biobrick part 1 could be considered orthogonal with respect to chassis' parts.

4. Conclusions

  • Surprisingly, standard fluorimetry assays show an assymetric output of red and green fluorescent proteins in co-transformed cells.
  • Flow cytometry confirms this asymmetric activity and strongly suggests the co-existence of sub-populations of cells with different outputs. This asymmetry proved very reproducible among cultures.
  • The activity of fluorescent proteins dramatically changes over time, reaching a peak when cells are in the stationary phase. Despite being under the control of the same promoter, GFP is expressed earlier than RFP. To the question “are two Biobricks orthogonal?” we thus suggest to add this other question: “when are they orthogonal?”
  • Proteomic analysis reveals the robustness of E. coli confronted to heterologous gene expression: only minor proteomic changes were found compared to the controls. An orthogonal relation between Biobrick part 1 and the proteomic architecture (chassis) is thus confirmed.

Open License (and Responsible Research and Innovation)

Not unlike scientific research, RRI is about common sense, about combining self-profit with others’ profit, about taking care of our planet. We believe that being responsible is not the opposite of being selfish; responsibility is just a broad-minded way of selfishness that points at the long-term benefit. And there is no long-term benefit without a sustainable society in a sustainable environment. What is good for our neighbors and for our home must be good for us as well.

Our Proposal: Responsible Research and Innovation as a tool to choose within the IP ecosystem.

Since we, our Human Practices sub-team, are a biotechnology and a law student, we are also glad to see that another key concept in IP is a beautiful science-law metaphor: Jane Calvert & Drew Endy's “diverse ecology” scenario. This image fits well with our interdisciplinary spirit and we found it very useful because it evokes very elegantly the range of different solutions out there to choose from in order to use a scientific invention in a given way.

Having reached this point, it is obvious which the next step is: there is a whole ecosystem of IP figures (see the poetic cover illustration by Paula, another team member), and one has to be chosen. How? Our bet is simple: by using RRI principles (mutual learning, human resources management, fundamental rights respect, research formation and transparency) as the main guide in the decision process.

The following figure illustrates what we propose:

Figure 1
Figure 1. RRI-based decision process on a diverse ecology scenario. From a diverse array of IP protection possibilities (diverse ecology scenario), we propose RRI as a framework to choose the most suitable one. In order to assist in the decision process, experts’ opinion is the major force, and we propose here two tools to contribute to this process: a common language exemplified by our dictionaries (Annex I), and a formula-based assistance procedure (Annex II).

In our view, the decision process has to fall upon the responsibility of the experts or examiners, as it is the case today. But we propose that a formation on RRI for patent examiners and other IP actors is introduced. In the same way that scientists need a basic law background for transferring their knowledge into applied solutions, we envisage a future in which everybody involved in the production of a given IP figure (a patent, for example, from scientists and engineers to ethics specialists, lawyers and other experts) should share a common interdisciplinary formation on RRI. A common language is the first step towards a common goal.

Additionally, we propose in this report two tools that might contribute to avoid misunderstandings and to make a simplified decision process available for a wide audience. In order to accomplish the former objective, we have prepared the two common-languages short dictionaries listed in Annex II. To achieve the latter, we propose, for the first time for the best of our knowledge, a math-based equation that could contribute to make IP issues available to a wide audience (for a complete description see Annex I). We have prepared the equation in two steps. First, we combined the patentability requirements (novelty, inventive step, industrial application) in a mathematically logical way to reach our draft equation (F0). Then, we submitted the formula to several experts to have their feedback, and thus modify F0 accordingly to yield F1 (this reflexive process, combining design with improvement rounds, elegantly evokes a biotechnology strategy called directed evolution).

Figure 2
Figure 2. Human decisions are only acceptable if no help is used? Left, chess player. Right, graphic representation of the equation F0 we initially proposed for patentability determination, as a function of the Industrial Application (I) and the Inventive step (A). See Annex II for details.

The very idea of using maths to decide whether something is patentable may seem odd. The opinion of the law experts was useful; but we wanted to have a broader feedback from other social actors. We thus carried out a survey to probe the opinion of iGEMers, scientists/engineers and law students, but also of general population to compare their view on our proposal. The C section of Annex II describes the results of this survey that can be summarize as follows: the vast majority of the 526 respondents (70 to around 80%) found that a math-based approach could be a useful tool, but should not substitute humans as the final decision makers. We could not agree more on that.

Report

Original watercolor by Paula Villaescusa.

PDF File Annex I. Patentability index PDF File Annex 2. Brief dictionary of IP and SB

The seat of the ST2OOL

One of the goals of our project was to select standard, stable, and orthogonal parts naturally occurring in nature. To do that, we carried out a functional metagenomics strategy aiming at selecting promoters from a library of metagenomic DNA. As a first step, we wanted to isolate metagenomic DNA from environmental samples, but instead of using traditional manual kits, we decided to build a robot able to automatically extract metagenomic DNA for us… It is a pleasure for our team to introduce you our TOOL robot:

“The TOOL” robot

How does it work?

The main goal of this machine was the automation of the DNA isolation process. This robot can be used for both E. coli (mini/maxipreps) and for DNA extraction from metagenomic samples.

The robotic system was developed in a continuous philosophy and all the components were assembled on a piece of wood. In short, it consisted of the following modules:

  • The control and powering system: the system was controlled by an Arduino Mega 2560 microchip, powered by a computer via USB. The other components of the machine were powered by an ATX power supply.
  • The injection system, which consisted of several screw stepper motors controlled by an A4988 driver. We chose this kind of motors because they allowed a quite accurate control of the position and the speed of the injections. A torque of 35N*m was also included. The system had three-way stopcocks in order to refill the syringes easily.
  • The temperature modules: high temperatures were achieved with a system that combined hot water and Peltier cells, whereas low temperatures were achieved by ice cooling. Temperature was measured with a dbs1820 temperature sensor.
  • The mixing: a homogeneous mixture of all the reagents was obtained thanks to the use of tubes allowing a turbulent regime.

pictures robot

The complete process was as follows:

  1. Manually inject a suspension of the cells into 3 mL of ice-cold solution I.
  2. From this step until step 8, the procedure is automatic (indicated in italics).
  3. 200 µL of solution II are injected and mixed during the flow into the silicone tubes.
  4. The mixture is moved through the silicone tubes to a water bath at 65 ºC and incubated for 30 minutes.
  5. 60 µL of ice-cold Potassium Acetate 3M (pH 5.0) are injected and mixed during the flow into the silicone tubes.
  6. The mixture is moved through the silicone tubes to an ice bath and incubated 20 minutes at -20ºC.
  7. The resulting suspension is then filtered through a membrane and insoluble particles are removed.
  8. An equal volume of isopropanol is injected and mixed with to the filtered solution during the flow into the silicone tubes, and the mixture is incubated at room temperature for 5 minutes.
  9. At this point, the solution is collected in Eppendorf tubes and these tubes are centrifuged at maximum speed for 10 minutes.
  10. The supernatant is discarded and the pellet washed with 500 µL of 70% (v/v) ethanol.
  11. The tubes are centrifuged again at maximum speed for 3 minutes and the supernatant is discarded.
  12. Finally, the pellet is dried and resuspended in 50 µL of water or elution buffer.

Results

As a proof of concept, we used a suspension of soil bacteria (obtained after several mild centrifugations of a suspension of compost in PBS sterile buffer) to isolate metagenomic DNA. Up to now, we have been able to perform steps 1-8 with our robot, and we are still working hard on the implementation of the final steps of the protocol. Figure 1A shows the Nanodrop analysis on the metagenomic DNA we obtained from the suspension. We isolated pure (A260/A280 and A260/A230 ratios close to 2) and quite concentrated (90.5 ng/µL) metagenomic DNA, the integrity of which was also checked on a 0.8% agarose gel (Figure 1B). Interestingly, we also observed the bands corresponding to the 16S and 23S ribosomal RNA on the gel.

Figure 1
Figure 1. A) Screenshot of the Nanodrop output of compost metagenomic DNA isolated with the robot as described in the main text. B) TBE 0.8% agarose gel showing the nucleic acids isolated with our robot. Lane 1: M.W. marker; lanes 2-5: metagenomic DNA aliquots. Arrows indicate metagenomic DNA (upper part) and 16S/23S ribosomal RNA.

Small budget, big ideas

Table 1 shows a list of all the components and materials used to build the robot. We decided to buy cheap and easy-to-find materials, so anyone can re-build or improve the robot! We only spent around 500€:

Table 1.Materials and components used during the construction of “The TOOL” robot. The price of each component is indicated on the right.
MATERIAL NUMBER Price (€)
ARDUINO MEGA 2560 1 12
WIRES 100 10
PROTOBOARD 1 3
DRIVER A4988 5 15
SCREW STEPPER MOTOR 4 130
SYRINGE 4 1
SCREW S 50 5
ALUMINUM BARS 1 3
WOOD TABLE 3 3
CLAMP 8 3
DRAWER GUIDE 4 3
SILICONE TUBES 5 1
PERISTALTIC PUMP 2 10
BEAKERS 3 7
SAW 12
FISH 1 2
DRILL 1 20
STOPCOCK 3 WAYS 6 12
SCROLL SAW 1 8
FILTER 1 per sample -
FILTER HOLDER (own design) 1 4
AID SILICONE TUBES 4 5
RESISTOR (400 h) 4 3
CAPACITOR(0.5F) 1 2
LIMIT SWITCH 2 3
WELDER 1 15
TIN COIL 1 11
RULE 1 3
LUB 1 6
SILICONE PISTOL 1 6
SILICONE 4 1
SCREW DRIVER 2 5
ATX POWER SUPPLY 1 10
dbs1820 SENSOR 2 6
PELTIER CELL 3 12
RELAY 1 6
PERISTALTIC PUMP 2 16
TRANSISTOR 2 1
Hours of work… A lot! A lot!

Selecting natural standard promoters

Once we had the metagenomic DNA extraction performed by our robot, we wanted to clone a library of fragments coming from this DNA in the pAcGFP1-1 plasmid (Clontech Laboratories, Inc.). This plasmid contains a GFP protein placed downstream the MCS, so if a promoter sequence is inserted in the polylinker, the expression of the fluorescent protein is triggered and thus fluorescent transformants can be easily detected.

We digested both the plasmid and the metagenomic DNA with the restriction enzymes EcoRI and BamHI and then purified the resulting fragments (Figure 2).

Figure 2
Figure 2. TBE 0.8% agarose gel showing the plasmid and the metagenomic DNA fragments obtained. Lane 1: MW marker; lane 2: undigested mgDNA; lane 3: digested mgDNA; lane 4: digested pAcGFP1-1 plasmid.

We set up a ligation reaction, incubated it overnight, and then transformed E. coli DH5α competent cells. We obtained around 300 colonies, but none was fluorescent… We performed a colony PCR with oligonucleotides flanking the MCS of the pAcGFP1-1 plasmid in order to check whether the transformant colonies were carrying inserts of metagenomic DNA, and confirmed their presence as you can see in Figure 3.

Figure 3
Figure 3. TBE 0.8% agarose gel showing the results of the colony-PCRs. Lanes 1 and 10: MW marker; lanes 2-9 and 11-16: independent colonies; lane 17: negative control.

We are working now on improving the efficiency of the transformation, so we can increase the probability of "catching" a good promoter sequence… we hope to succeed soon!

We summarize here our main contributions modeling ST$^2$OOL. Our hypothesis, formulae and analysis can be found in the modeling section.

Stability We have modeled the variation of fluorescence per cell of Biobricks with temperature, pH and salt concentration by considering changes in their transcription and synthesis rates. Our model captures the measured shape of the fluorescence per cell as a function of temperature. We have found that transcription rate goes up with temperature by the increase of free σ factors and synthesis rate goes down above the optimal temperature. Unfortunately, this is not all the story. Our model dramatically fails to fit the amount of variation, which is much larger experimentally than the predicted one. Let us compare theoretical predictions and experimental results. In both studies, strain XL1-Blue transformed with Biobricks Bb1 or Bb3, we find that the ratio between theoretical prediction and experimental result is the same at temperatures lower than the optimal (in fact, accidentally too close given the error bars). This fact suggests that the extra temperature dependence not included in our model should come from cell mechanisms independent on parameters of the Biobricks (promoter strength, ...). Quite on the contrary, the comparison of theoretical predictions and experimental results at larger temperatures show that our model has not yet capture the dependence of the magnitude of the decrease in activity with the Biobrick parameters.

Standardization We have built an app that computes the standardization-stability index and generates a promoter sequence for a given expression rate. Our wet-lab results show that de Biobrick 5 is the more standard one. Our model predicts that systems with less modules or factors will be the more stable and standard. Biobrick 5 has a tetracycline depending promoter, whose limitating step is the binding of tetracycline to TetR. In this scenario, the intracellular concentration of aTc and the binding to the TetR repressor do not depend on the strain, so there is only one factor in which the strain could have an important effect: the TetR intracellular concentration, which has a constitutive expression. Analyzing the other Biobricks, we found that the strains with less expression in the constitutive promoter are the ones with higher expression when we use the Ptet promoter (but JM109 and BB1), which agree with our prediction.

Orthogonality We have studied deterministic and stochastic models of the plasmid growth to learn about the unspecific interaction between two plasmids in a cell. We found that this system is unstable and the only way to keep the plasmid levels is introducing a resistance. However, it is not a good way to maintain those levels: bacteria which have non-desired levels of plasmids are killed. This leads to a lower yield. Also, there is an important variability, which can drive the system away from the expected behaviour. In conclusion, plasmids with identical origin of replication are not orthogonal per se. Although we can introduce use different resistances to tune the system, it could not behave as expected due to a great variability. As a solution, we could use plasmids with different mechanisms of replication, but they could show different dynamics, making the system even mor difficult to control. Therefore, new plasmids must be designed to make engineered cells more reliable and predictable to make progress in synthetic biology.

Open License We have modeled the susceptibility to patent a human creation and define the factors that influence it. Inventive is the most relevant variable but is regulated by the industrial application. The highest inventive step can be reduced by almost a factor 2 with the lowest industrial application. On the contrary, the lowest inventive step is not affected by the industrial application. After computing all the discussed factors, the patentability index can run from 0 to 10. The reader can find a detailed discussion of the questionnaires, tests, analysis of the index and references in the area Human Practices.