Team:TU-Delft/Sensing

From 2013.igem.org

(Difference between revisions)
 
(33 intermediate revisions not shown)
Line 6: Line 6:
-
<div style="margin-left:30px;margin-right:30px; width:900px;float:left;">   
+
<div style="margin-left:30px;margin-right:30px;margin-top:10px; width:900px;float:left;">   
<h2 align="center">Sensing Device</h2>
<h2 align="center">Sensing Device</h2>
<html>
<html>
<center>
<center>
-
<a href="https://static.igem.org/mediawiki/2013/6/6d/Sensing-AIP.jpg" target="_blank"><img src="https://static.igem.org/mediawiki/2013/6/6d/Sensing-AIP.jpg"width="500px"  height="354px"/></a>
+
<a href="https://static.igem.org/mediawiki/2013/f/fa/PBADRecGFP.jpg" target="_blank"><img src="https://static.igem.org/mediawiki/2013/f/fa/PBADRecGFP.jpg"width="500px"  height="354px"/></a>
</center>
</center>
<center>
<center>
Line 18: Line 18:
<div style="margin-left:30px;margin-right:30px; width:900px;float:left;">   
<div style="margin-left:30px;margin-right:30px; width:900px;float:left;">   
<p align="justify">  
<p align="justify">  
-
In bacteria quorum sensing is an important mechanism for monitoring the state of the population. <i>S. aureus</i> uses its quorum sensing system to regulate a large part of its virulence genes, as having a large population is beneficial for <i>S. aureus</i> during infection. In this project we highjack the native <i>S. aureus</i> system and implement that in the antimicrobial peptide producing host <i>E.coli</i>.  
+
In bacteria quorum sensing is an important mechanism for monitoring the state of the population. <i>S. aureus</i> uses its quorum sensing system to regulate a large part of its virulence genes, as having a large population is beneficial for <i>S. aureus</i> during infection <a href="https://2013.igem.org/Team:TU-Delft/Sensing#references">[1]</a>. In this project we hijack the native <i>S. aureus</i> system and implement that in the antimicrobial peptide expression host <i>E.coli</i>.  
</p>
</p>
<p align="justify">  
<p align="justify">  
Line 24: Line 24:
</p>
</p>
<p align="justify">  
<p align="justify">  
-
The system is composed out of the transmembrane receptor protein AgrC and the cytoplasmic protein AgrA. AgrA is phosphorylated by AgrC upon induction with auto-inducing peptide (AIP), after which it will act as a transcription factor positively acting on the promoter P2. AIP is made from the precursor peptide AgrD, which is circularized and secreted by AgrB (Figure 1).
+
The system is composed out of the transmembrane receptor protein AgrC and the cytoplasmic protein AgrA. AgrA is phosphorylated by AgrC upon induction with auto-inducing peptide (AIP), after which it will act as a transcription factor positively acting on the promoter P2. AIP is made from the precursor peptide AgrD, which is circularized and secreted by AgrB (Figure 1). In the final construct used in Peptidor, pP2 will not be in front of GFP, but in front of the SUMO-peptides as explained <a href="https://2013.igem.org/Team:TU-Delft/Peptides" target="blank">here</a>.
</p>
</p>
-
<br>
+
 
-
<h4>AIP Sensing Protocol: </h4>
+
-
<ol>
+
-
<li> Grow the construct    <a href="http://parts.igem.org/Part:BBa_K1022100" style="text-decoration: none"" target="_blank">BBa_K1022100</a> (pBAD AIP Receiver GFP) overnight. </li>
+
-
<li> Make 1/50 dilutions of the overnight grown culture. Check for OD at 600nm till 0.1.</li>
+
-
<li> Then induce with 0.1% Arabinose. Keep aside samples with No Arabinose, No AIPs induction to use as control. BL21 cells and Const GFP are also used as negative and positive controls.</li>
+
-
<li> Track the OD till 0.5 and then induce with AIP 1 µM and 10 µM. </li>
+
-
<li> After 3 hours of incubation time, check for GFP signals on the FACS (Fluorescence-activated cell sorting). </li>
+
-
</ol>
+
<br>
<br>
<h4>AIP Sensing Experiment: </h4>
<h4>AIP Sensing Experiment: </h4>
<p align="justify">
<p align="justify">
-
The construct <a href="http://parts.igem.org/Part:BBa_K1022100" style="text-decoration: none"" target="_blank">BBa_K1022100</a>  (pBAD AIP Receiver GFP) is the test construct for the first module of the project in order to sense the auto-inducing peptides (AIPs) produced by <i>Staphylococcus</i> aureus.  
+
The construct <a href="http://parts.igem.org/Part:BBa_K1022100" style="text-decoration: none"" target="_blank">BBa_K1022100</a>  ,pBAD AIP Receiver GFP, (fig. 2) is the test construct for the first module of the project in order to sense the auto-inducing peptides (AIPs) produced by <i>S. aureus</i>.  
</p>
</p>
<h4>Diagram of the Part</h4>
<h4>Diagram of the Part</h4>
-
In this sensor test-construct AgrC and AgrA are expressed, after which AgrC can act as a receptor for AIP. After binding AIP, AgrC will phosphorylate AgrA which then starts acting as a postitive transcription factor for pP2 <a href="https://2013.igem.org/Team:TU-Delft/Sensing#references">[1]</a>. Therefore, upon induction by AIP, the construct will give a GFP signal. This is tested using the following experiment.  
+
In this sensor test-construct AgrC and AgrA are expressed, after which AgrC can act as a receptor for AIP. After binding AIP, AgrC will phosphorylate AgrA which then starts acting as a postitive transcription factor for pP2 <a href="https://2013.igem.org/Team:TU-Delft/Sensing#references">[2]</a>. Therefore, upon induction by AIP, the construct will give a GFP signal. This is tested using the following experiment.  
 +
 
 +
<center>
 +
<img src="https://static.igem.org/mediawiki/2013/3/38/Part100.jpg">
 +
<p>Figure 2: Part:BBa_K1022100 </p>
 +
</center>
 +
 
 +
 
<br><br>
<br><br>
<h4>Experimental Set Up</h4>
<h4>Experimental Set Up</h4>
<p align="justify">
<p align="justify">
-
Cells are induced with 0.1 % Arabinose (to induce pBAD) and AIP’s, and with only AIP. Induction with AIP is done at 1µM and 10 µM. 3 hours after induction fluorescence-activated cell sorting (FACS) is used to check for GFP signals. The controls are BL21(DE3) cells, uninduced plasmid containing BL21(DE3) and a constitutively expressed GFP construct containing strain. The AIP used to induce the cells was synthesized by Eurogentec.
+
Cells are induced with 0.1 % Arabinose (to induce pBAD) and AIP, and with only AIP. AIP Induction is done at O.D 0.5 with 10 µM AIP. 3 hours after induction fluorescence-activated cell sorting (FACS) is used to check for GFP signal. The controls are BL21(DE3) cells, AIP receiver GFP containing BL21(DE3)without any induction and a strain constitutively expressing GFP. The AIP used to induce the cells was synthesized by Eurogentec.
</p>
</p>
 +
 +
<a href="https://2013.igem.org/Team:TU-Delft/Protocol_13#protocol_13" target="_blank"><i>Full protocol</i></a>
<p align="justify">
<p align="justify">
-
The cells induced with arabinose all died, which could be seen by the fact the signal given by the FACS, this could be due the fact AgrA and/or AgrC are toxic in high concentrations. The fact <i>E. coli</i> is still able to detect AIPs without induction of the Agr genes is explainable due the fact pBAD is a bit leaky, making low transcription possible, apparently only a small amount of activated AgrA is required to give a measurable level of GFP.  
+
The cells induced with arabinose all died, which could be seen by the fact the signal given by the FACS. FACS gives a fluorescent signal for every cell passing the reader, live cells give a higher signal than cell debris. From the fact the arabinose induced cells showed a lower signal than the BL21(DE3) control we concluded these cells were dead. This cell lysis could be due the fact AgrA and/or AgrC is toxic in high concentrations. The detection of AIP by <i>E. coli</i> without induction of the Agr genes is explainable due the fact pBAD is a bit leaky, making low levels of transcription possible. Apparently only a small amount of activated AgrA is required to give a measurable level of GFP.  
</p>
</p>
<p align="justify">
<p align="justify">
-
Earlier work on the detection of AIPs by Gram-negatives has been done, but with unclear or non-convincing results <a href="https://2013.igem.org/Team:TU-Delft/Sensing#references">[1,2]</a>. Often the outer-membrane of E. coli was seen as an obstacle for the implementation of auto-inducing peptide sensing systems in Gram negatives. Although, after intensive literature research we found a molecule very similar to the AIP molecule used in these experiments, that is known to pass the outer membrane. The molecule bacitracin [4] has been shown to interfere with cell wall synthesis in Gram negative species, therefore being able to pass the outer membrane as cell wall synthesis occurs in the periplasm. Comparison of the molecular structures of bacitracin and AIP regarding size, molecular weight, polarity and charge gave us an indication AIP should be possible to pass through the outer membrane of <i>E. coli</i> (fig 1).
+
Earlier work on the detection of AIPs by Gram-negatives has been done, but with unclear or non-convincing results <a href="https://2013.igem.org/Team:TU-Delft/Sensing#references">[3,4]</a>. Often the outer-membrane of E. coli was seen as an obstacle for the implementation of auto-inducing peptide sensing systems in Gram negatives. In the AIP sensing experiments preformed by the iGEM Cambridge 2007 team a FepA construct was used as this is known to make the outer membrane more permeable<a href="https://2013.igem.org/Team:TU-Delft/Sensing#references">[8]</a>. Although, after intensive literature research we found a molecule very similar to the AIP molecule used in these experiments, that is known to pass the outer membrane. The molecule bacitracin [5] has been shown to interfere with cell wall synthesis in Gram negative species, therefore being able to pass the outer membrane as cell wall synthesis occurs in the periplasm. Comparison of the molecular structures of bacitracin and AIP regarding size, molecular weight, polarity and charge (fig. 3) gave us an indication AIP should be possible to pass through the outer membrane of <i>E. coli</i> without any modifications done on the outer membrane.
</p>
</p>
 +
 +
<center>
 +
<img src="https://static.igem.org/mediawiki/2013/1/13/BacitracinAip.png">
 +
<p >Figure 3: Molecule structures of bacitracin and autoinducing peptides  </p>
 +
</center>
 +
<h4>Results</h4>
<h4>Results</h4>
<p align="justify">
<p align="justify">
-
Using Cyflogic 1.2.1 software <a href="https://2013.igem.org/Team:TU-Delft/Sensing#references">[5]</a>, FACS data was analysed, showing a clear induction of the AIP receiver GFP construct by AIP. The difference between the blue and the red line (218) indicates the difference in fluorescent signal of uninduced and induced cells. Untransformed Bl21(DE3) and constitutively expressed GFP were used as respectively negative and positive control (fig. 2). The significant difference between the uninduced and the induced samples is determined with a P-value >10-4 as described in the next paragraph.  
+
Using Cyflogic 1.2.1 software <a href="https://2013.igem.org/Team:TU-Delft/Sensing#references">[6]</a>, FACS data was analysed, showing a clear induction of the AIP receiver GFP construct by AIP. The difference between the blue and the red line (218) indicates the difference in fluorescent signal of uninduced and induced cells. Untransformed Bl21(DE3) and constitutively expressed GFP were used as respectively negative and positive control (fig. 4). The significant difference between the uninduced and the induced samples is determined with a P-value >10<sup>-4</sup> as described in the next paragraph.  
</p>
</p>
Line 65: Line 72:
<img src="https://static.igem.org/mediawiki/2013/f/f0/FinalSensing.jpg">
<img src="https://static.igem.org/mediawiki/2013/f/f0/FinalSensing.jpg">
<p align="justify">
<p align="justify">
-
Fig. 2. Fluorescent signal measured with the FACS. The white line refers to the control BL21(DE3) cells. In blue the uninduced pBAD AIP receiver GFP is shown. The red line represents the AIP induced receiver cells. The fluorescence intensity scale is logarithmic and the total number of cells counted per sample is 10,000.
+
Figure 4: Fluorescent signal measured with the FACS. The white line refers to the control BL21(DE3) cells. In blue the uninduced pBAD AIP receiver GFP is shown. The red line represents the AIP induced receiver cells. The green line shows cells constitutively expressing GFP. The fluorescence intensity scale is logarithmic and the total number of cells counted per sample is 10,000.
</p>
</p>
</center>
</center>
Line 71: Line 78:
<h4>Significance test</h4>
<h4>Significance test</h4>
<p align="justify">
<p align="justify">
-
To test if the difference between the uninduced and induced result is significant a statistic test is performed as in <a href="https://2013.igem.org/Team:TU-Delft/Sensing#references">[6]</a>. Here it is assumed that both results are normal distributions, as the number of datasets is large (central limit theorem) <a href="https://2013.igem.org/Team:TU-Delft/Sensing#references">[6]</a>. Then the critical value equation becomes a simple one:
+
To test if the difference between the uninduced and induced result is significant a statistic test is performed as in <a href="https://2013.igem.org/Team:TU-Delft/Sensing#references">[7]</a>. Here it is assumed that both results are normal distributions, as the number of datasets is large (central limit theorem) <a href="https://2013.igem.org/Team:TU-Delft/Sensing#references">[7]</a>. Then the critical value equation becomes a simple one:
</p>
</p>
<center>
<center>
-
<img src="2013.igem.org/wiki/images/7/7f/EquationDelft.png">
+
<img src="https://static.igem.org/mediawiki/2013/a/a2/Eq_stat.jpg" width="250">
<p align="justify">
<p align="justify">
-
Based on the induced and uninduced data as in Table 1, this t-value is calculated as 52.7. This is a very high value and corresponds to a one-sided p value of less than 10-4. This means that there is a significant difference between the induced and uninduced curves.  
+
Based on the induced and uninduced data as in Table 1, this t-value is calculated as 52.7. This is a very high value and corresponds to a one-sided p value of less than 10<sup>-4</sup>. This means that there is a significant difference between the induced and uninduced curves.  
</p>
</p>
</center>
</center>
Line 184: Line 191:
                 <table >
                 <table >
                     <tr>
                     <tr>
-
                        <td></td>
 
                         <td>
                         <td>
                       Dataset
                       Dataset
Line 203: Line 209:
                         <td>
                         <td>
                           242.6
                           242.6
 +
                        </td>
 +
                        <td>10.000</td>
 +
                   
 +
                    </tr>
 +
    <tr>
 +
                        <td>
 +
                        Induced
 +
                        </td>
 +
                        <td >
 +
                          683
 +
                        </td>
 +
                        <td>
 +
                          334.5
                         </td>
                         </td>
                         <td>10.000</td>  
                         <td>10.000</td>  
                        
                        
                     </tr>
                     </tr>
-
 
</table>
</table>
Line 214: Line 232:
<h4>Discussion</h4>
<h4>Discussion</h4>
<p align="justify">
<p align="justify">
-
Taking the results described above into account, we are able to sense the presence of a Gram positive species with the use of E. coli transformed with our AIP receiver GFP construct. To our knowledge, we are the first to show the interspecies communication between Gram negative and Gram positive bacterial species. In this case the sensing of the AIP of <i>S. aureus</i> by <i>E. coli</i>. For the final application described in the introduction (link) page though, the level of activation of the pP2 promoter by AgrA might be too low. In order to make the sensing system able to activate a peptide production system capable of producing a concentration high enough to kill or inhibit the growth of S. aureus, an amplification circuit could be used.  
+
Taking the results described above into account, we are able to sense the presence of a Gram positive species with the use of <i>E. coli</i> transformed with our AIP receiver GFP construct. To our knowledge, we are the first to show the interspecies communication between Gram negative and Gram positive bacterial species. In this case the sensing of the AIP of <i>S. aureus</i> by <i>E. coli</i>. For the final application described in the introduction (link) page though, the level of activation of the pP2 promoter by AgrA might be too low. In order to make the sensing system able to activate a peptide production system capable of producing a concentration high enough to kill or inhibit the growth of S. aureus, an amplification circuit could be used.  
</p>
</p>
Line 220: Line 238:
<h2 align="center">References</h2>
<h2 align="center">References</h2>
<ol>
<ol>
 +
<li>Textbook of bacteriology,available on:<a href="http://textbookofbacteriology.net/staph_2.html" target="blank"> www.textbookofbacteriology</a></li>
<li>Philip Mdowell, Zina Affas, et al.,Structure, activity and evolution of the group I thiolactone peptide quorum-sensing system of Staphylococcus aureus, Molecular Microbiology, Volume 41, Issue 2, pages 503–512, Jul 2001.</li>
<li>Philip Mdowell, Zina Affas, et al.,Structure, activity and evolution of the group I thiolactone peptide quorum-sensing system of Staphylococcus aureus, Molecular Microbiology, Volume 41, Issue 2, pages 503–512, Jul 2001.</li>
<li><a href="https://2007.igem.org/wiki/index.php/Cambridge" target="blank">https://2007.igem.org/wiki/index.php/Cambridge</a></li>
<li><a href="https://2007.igem.org/wiki/index.php/Cambridge" target="blank">https://2007.igem.org/wiki/index.php/Cambridge</a></li>
<li><a href="https://2012.igem.org/Team:HIT-Harbin" target="blank">https://2012.igem.org/Team:HIT-Harbin</a></li>
<li><a href="https://2012.igem.org/Team:HIT-Harbin" target="blank">https://2012.igem.org/Team:HIT-Harbin</a></li>
<li>T.J. Pollock, L. Thorne, Mechanism of bacitracin resistance in gram-negative bacteria that synthesize exopolysaccharides, J Bacteriol ; 176(20): 6229–6237, Oct 1994 </li>
<li>T.J. Pollock, L. Thorne, Mechanism of bacitracin resistance in gram-negative bacteria that synthesize exopolysaccharides, J Bacteriol ; 176(20): 6229–6237, Oct 1994 </li>
-
<li><a href="http://www.cyflogic.com/" target="blank">http://www.cyflogic.com/</a>
+
<li>Cyflogic 1.2.1., available on:<a href="http://www.cyflogic.com/" target="blank">http://www.cyflogic.com/</a>.
</li>
</li>
-
<li>F.M. Dekker, C. Kraaikaamp et al., ‘A Modern Introduction to Probability and Statistics’, January 2005, Springer.</li>
+
<li>F.M. Dekker, C. Kraaikaamp et al., ‘A Modern Introduction to Probability and Statistics’, Springer,January 2005.</li>
 +
<li><a href="https://2007.igem.org/wiki/index.php/Cambridge/Signalling_project#FepA_permeability_device" target="blank">www.2007.igem.org/wiki/cambridge</a></li>
</ol>
</ol>

Latest revision as of 14:18, 4 October 2013


Sensing Device

Figure 1: Schematic diagram of the sensing part


In bacteria quorum sensing is an important mechanism for monitoring the state of the population. S. aureus uses its quorum sensing system to regulate a large part of its virulence genes, as having a large population is beneficial for S. aureus during infection [1]. In this project we hijack the native S. aureus system and implement that in the antimicrobial peptide expression host E.coli.

We used the quorum sensing system of S. aureus in our peptidor in order to make sure the antimicrobial peptides will only be produced when S. aureus is present. This way there is less chance of other microbes becoming resistant against the peptides and therefore indirectly, through horizontal gene transfer, preventing S. aureus becoming resistant.

The system is composed out of the transmembrane receptor protein AgrC and the cytoplasmic protein AgrA. AgrA is phosphorylated by AgrC upon induction with auto-inducing peptide (AIP), after which it will act as a transcription factor positively acting on the promoter P2. AIP is made from the precursor peptide AgrD, which is circularized and secreted by AgrB (Figure 1). In the final construct used in Peptidor, pP2 will not be in front of GFP, but in front of the SUMO-peptides as explained here.


AIP Sensing Experiment:

The construct BBa_K1022100 ,pBAD AIP Receiver GFP, (fig. 2) is the test construct for the first module of the project in order to sense the auto-inducing peptides (AIPs) produced by S. aureus.

Diagram of the Part

In this sensor test-construct AgrC and AgrA are expressed, after which AgrC can act as a receptor for AIP. After binding AIP, AgrC will phosphorylate AgrA which then starts acting as a postitive transcription factor for pP2 [2]. Therefore, upon induction by AIP, the construct will give a GFP signal. This is tested using the following experiment.

Figure 2: Part:BBa_K1022100



Experimental Set Up

Cells are induced with 0.1 % Arabinose (to induce pBAD) and AIP, and with only AIP. AIP Induction is done at O.D 0.5 with 10 µM AIP. 3 hours after induction fluorescence-activated cell sorting (FACS) is used to check for GFP signal. The controls are BL21(DE3) cells, AIP receiver GFP containing BL21(DE3)without any induction and a strain constitutively expressing GFP. The AIP used to induce the cells was synthesized by Eurogentec.

Full protocol

The cells induced with arabinose all died, which could be seen by the fact the signal given by the FACS. FACS gives a fluorescent signal for every cell passing the reader, live cells give a higher signal than cell debris. From the fact the arabinose induced cells showed a lower signal than the BL21(DE3) control we concluded these cells were dead. This cell lysis could be due the fact AgrA and/or AgrC is toxic in high concentrations. The detection of AIP by E. coli without induction of the Agr genes is explainable due the fact pBAD is a bit leaky, making low levels of transcription possible. Apparently only a small amount of activated AgrA is required to give a measurable level of GFP.

Earlier work on the detection of AIPs by Gram-negatives has been done, but with unclear or non-convincing results [3,4]. Often the outer-membrane of E. coli was seen as an obstacle for the implementation of auto-inducing peptide sensing systems in Gram negatives. In the AIP sensing experiments preformed by the iGEM Cambridge 2007 team a FepA construct was used as this is known to make the outer membrane more permeable[8]. Although, after intensive literature research we found a molecule very similar to the AIP molecule used in these experiments, that is known to pass the outer membrane. The molecule bacitracin [5] has been shown to interfere with cell wall synthesis in Gram negative species, therefore being able to pass the outer membrane as cell wall synthesis occurs in the periplasm. Comparison of the molecular structures of bacitracin and AIP regarding size, molecular weight, polarity and charge (fig. 3) gave us an indication AIP should be possible to pass through the outer membrane of E. coli without any modifications done on the outer membrane.

Figure 3: Molecule structures of bacitracin and autoinducing peptides

Results

Using Cyflogic 1.2.1 software [6], FACS data was analysed, showing a clear induction of the AIP receiver GFP construct by AIP. The difference between the blue and the red line (218) indicates the difference in fluorescent signal of uninduced and induced cells. Untransformed Bl21(DE3) and constitutively expressed GFP were used as respectively negative and positive control (fig. 4). The significant difference between the uninduced and the induced samples is determined with a P-value >10-4 as described in the next paragraph.

Figure 4: Fluorescent signal measured with the FACS. The white line refers to the control BL21(DE3) cells. In blue the uninduced pBAD AIP receiver GFP is shown. The red line represents the AIP induced receiver cells. The green line shows cells constitutively expressing GFP. The fluorescence intensity scale is logarithmic and the total number of cells counted per sample is 10,000.


Significance test

To test if the difference between the uninduced and induced result is significant a statistic test is performed as in [7]. Here it is assumed that both results are normal distributions, as the number of datasets is large (central limit theorem) [7]. Then the critical value equation becomes a simple one:

Based on the induced and uninduced data as in Table 1, this t-value is calculated as 52.7. This is a very high value and corresponds to a one-sided p value of less than 10-4. This means that there is a significant difference between the induced and uninduced curves.

Table 1: Data of the induced and non-induced FACS
Dataset Mean (μ) Standard deviation (σ) Number of samples (N)
Non-induced 465 242.6 10.000
Induced 683 334.5 10.000

Discussion

Taking the results described above into account, we are able to sense the presence of a Gram positive species with the use of E. coli transformed with our AIP receiver GFP construct. To our knowledge, we are the first to show the interspecies communication between Gram negative and Gram positive bacterial species. In this case the sensing of the AIP of S. aureus by E. coli. For the final application described in the introduction (link) page though, the level of activation of the pP2 promoter by AgrA might be too low. In order to make the sensing system able to activate a peptide production system capable of producing a concentration high enough to kill or inhibit the growth of S. aureus, an amplification circuit could be used.

References

  1. Textbook of bacteriology,available on: www.textbookofbacteriology
  2. Philip Mdowell, Zina Affas, et al.,Structure, activity and evolution of the group I thiolactone peptide quorum-sensing system of Staphylococcus aureus, Molecular Microbiology, Volume 41, Issue 2, pages 503–512, Jul 2001.
  3. https://2007.igem.org/wiki/index.php/Cambridge
  4. https://2012.igem.org/Team:HIT-Harbin
  5. T.J. Pollock, L. Thorne, Mechanism of bacitracin resistance in gram-negative bacteria that synthesize exopolysaccharides, J Bacteriol ; 176(20): 6229–6237, Oct 1994
  6. Cyflogic 1.2.1., available on:http://www.cyflogic.com/.
  7. F.M. Dekker, C. Kraaikaamp et al., ‘A Modern Introduction to Probability and Statistics’, Springer,January 2005.
  8. www.2007.igem.org/wiki/cambridge