"
Team:XMU-China/Content4
From 2013.igem.org
| Line 185: | Line 185: | ||
<div class="clear"></div><center> | <div class="clear"></div><center> | ||
| - | <img src="https://static.igem.org/mediawiki/2013/a/a7/XMU-China_s14.png" width = | + | <img src="https://static.igem.org/mediawiki/2013/a/a7/XMU-China_s14.png" width =900px height=150px class="border alignleft" alt="Figure 3" style="margin-left:50px"/></center> |
<div class="clear"></div> | <div class="clear"></div> | ||
| - | <p><small style="font-size:16px"></small></p | + | <p><small style="font-size:16px"> Figure 14</small></p> |
| - | + | ||
<img src="https://static.igem.org/mediawiki/2013/6/60/XMU-China_s13.png" width =350px height=300px class="border alignleft" alt="Figure 2" style="margin-left:140px"/> | <img src="https://static.igem.org/mediawiki/2013/6/60/XMU-China_s13.png" width =350px height=300px class="border alignleft" alt="Figure 2" style="margin-left:140px"/> | ||
<img src="https://static.igem.org/mediawiki/2013/f/fe/XMU-China_s15.png" width =350px height=300px class="border alignleft" alt="Figure 3" /></center> | <img src="https://static.igem.org/mediawiki/2013/f/fe/XMU-China_s15.png" width =350px height=300px class="border alignleft" alt="Figure 3" /></center> | ||
<div class="clear"></div> | <div class="clear"></div> | ||
| + | <p><small style="font-size:16px"> Figure 14 & 15 show the results given by Gene OS.</small></p> | ||
| + | <div class="clear"></div><center> | ||
| - | |||
| - | |||
<br /><br /> | <br /><br /> | ||
<h4>ImageMe: Comparing the results with manual calculation</h4> | <h4>ImageMe: Comparing the results with manual calculation</h4> | ||
Revision as of 05:25, 27 September 2013
Due to the complexity of our circuit, it is easier to understand the mechanism through modeling. In addition, the suitable modeling gives the results more quickly and correctly. A good mathematic model can also help us compare the affecting strength of different factors.
Software & Function Introduction
Based on the DDEs, “Gene oscillate simulator” has two main modes of simulations. The basic simulation is designed for the oscillate circuit with LuxI, aiiA and LuxR, while the advanced simulation adds the effect of ndh gene.
Figure 1 shows the first
interface of the Gene OS.
In each of the modes, users can make a series of graphs though designing the different parameters. The numbers of graphs mainly depend on the different parameters and the modes.
Figure 2 shows the basic simulation user interface. Figure 3 shows the choose dialog box.
The meaning of each parameter can be displayed when click each parameter or right click each input box. For the people first use this software, we designed the “One example” function which can set the parameters suitably in a single time.
To get an idea about the different conditions in which oscillations could occur, we created a function which allows users to enter the range in which the variables should be varied and compare their results in the results. Due to the arithmetic speed, we finally apply this function in fluid flow speed and cell density, which are the most important parameters in the model. For example, people can set the cell density range from 0.1 to 0.9 with the step 0.2 by not selecting the “Fix the cell density” item.
A white noise system was coded to simulate the noisy environment in the medium. This function makes the model closer to the actual environment. Users can set the strength by clicking “add white noise” box and inputting the decibels.
A script calculates amplify and period of the oscillation is also added as an additional function. This allows our software detect and compare the characteristic of certain oscillation.
After setting all the parameters, users can save and output these settings into a standard 2003 excel file through the file menu. Loading is also supported which makes more convince for transmitting the parameters.
Before draw the graphs users need to choose the plot mode when range parameters exist. For example, when set the cell density the software will ask user to draw separately or put curves with different cell density together.
Figure 4 & 5 show some example output of Gene OS.
Then clicking the “Draw” button, the software iterates over the values and plots graphs of all combinations possible for that range of values. After calculating, at most 3 kinds of graphs can be given, (1) concentration alteration of four main protein (2) fluorescence alteration under different cell density and fluid flow speed (3) amplify and period alteration under different cell density and fluid flow speed.
To protect the software and computer, we designed an error-stop system which could stop the simulation as soon as one error was detected. This system also is called when users click “Cancel” button to stop the calculation. But different with error happened; this stop process takes about 1 min. The arithmetic speed of this software depends on the complexity of the parameters. Choosing both cell density range and fluid flow range to calculate can sometimes take more than 10 minutes. It is mainly because the large numbers of nested for-loops in the script and the slowly calculate speed about MATLAB core.
You can download the MATLAB code here.
Figure 6 shows some example output of Gene OS.
Background
We usually got a great number of fluorescent images every time when we finished the 8-hours-long microfluidic experiment. However, to get the data we have to calculate each tray in each image though using the software. What’s worse, the software could only process one image each time. Obviously, a more humanized software is needed if we do not want to spend much time on the manual work.
Figure 7 shows the first interface of the ImageMe.
Coding in MATLAB
As we have not found such software on the net, we decided to design our own software. Thus, we used MATLAB coding the script and designed the user interface again for easy use.
We called this one “ImageMe” which was related with our team name. The most advantage of ImageMe is that it can analyze all the trays in the all images once.
Arithmetic
Thanks to the powerful image processing ability of MATLAB, the main arithmetic is concise and efficient. First ImageMe changes the 256-color image to gray scale image and then uses edge processing to make the edge of each tray easier to detect. This process also reduces the negative effect of the background. The final step is to calculate the total gray level using numerical integration.
Figure 8 shows the main interface. Figure 9 shows the input dialog box.
Software & function introduction
Comparing with the Gene OS, ImageMe has more concise interface arrangement and easier to use. Click “Start” to enter the main interface.
This is the main interface of ImageMe. Here, users can load the images by clicking the “Load” bottom and then choose the images. Then, ImageMe will pre-read all the images with “.jpg” suffix in this root directory. Make sure that no other .jpg images exist in the same folder. Click “Remove” bottom before changing the images to be analyzed.
After clicking the “Begin” bottom, ImageMe will ask users to input two parameters (which are also the only two need to input through the whole process) to assist the analysis more correctly. The two parameters are very easy to get comparing with those in Gene OS. First one is the lines of the trays and the second one is the rows. To get correct and precise results, please do not change the view or magnification times while taking images.
Figure 10 to 12 show the graphs given by ImageMe.
It is common that some images are indistinct or the edges are difficult to detect. When these situations occur to one image, ImageMe will use fuzzy calculation to analyze this one and make sure the result has the most correctness.
The real time message will be written on the screen when some indistinct images are detected or some errors happened.
If the analysis finished successfully, the graphs will show on the screen automatically. At most two kinds of graphs can be given, (1) fluorescence strength change in each tray (2) fluorescence strength change in each line.
You can download the MATLAB code here. A package of test images is included.
Gene OS: Comparing various results with different cell density
This example illustrates the importance of the cell density by comparing various results with different cell density.
Cell density is been normalized so that it is greater than 0 and less than or equal to 1. Thus we vary the cell density from 0.1 to 0.9 with the step-size 0.25 and keep other parameters the same.
Then, choose “Cell density” and start calculation. From the results we can see that under this specific condition too small or too big cell density all cannot lead to a steady oscillation. Only the cell density near 0.35 gives the expected curve. This is mainly because that cell in small density cannot product enough LuxR, AHL, aiiA & GFP to start the oscillation. This can be confirmed by the little fluorescence strength. While too many cells surely lead to high fluorescence strength and slow mass transfer speed which cover up the oscillation.
Figure 13 shows the whole parameters.
Figure 14
Figure 14 & 15 show the results given by Gene OS.
ImageMe: Comparing the results with manual calculation
To confirm the accuracy of ImageMe, we compared the two group results with manual calculated data. Due to the different arithmetic, the fluorescence strength is not equal, but the tendency indicated in the each graph is comparable.
First we compare them under the weak fluorescence situation.
From above, we can see that the general trend of each graph is both decreasing. However, ImageMe shows that the fluorescence strength of the last image is almost half of the first one while the manual calculation results only indicate a 15 percent decrease. The results are more amazing when we compare the 44th image (middle of the group) and the last one. ImageMe shows a 40 percent decrease while manual calculation results even implies a small increase in the fluorescence strength.
Now we compare these results with some of the images.
Thus, from these images we can conclude that the fluorescence strength is decreasing significantly from the first one to the last. In addition, it is obviously that fluorescence strength of the last image is less than the middle one which means a great computation error exists in the manual calculation.
Second, we compare them under the strong fluorescence situation.
In the second group, ImageMe shows a total 13 percent fluorescence strength increase while manual calculation indicates 10 percent. Comparing the 58th image (minimum of fluorescence strength) with the last one, ImageMe shows a 29 percent increase while manual calculation shows 32 percent. From the analysis, we can see the general trends of the two graphs match better than the first group.
Now we compare these results with some of the images again.
We choose the 54th and 58th (time 270 & 290) image because ImageMe shows a strange peak at 54th and a following valley at 58th.
From these images we can observe that the fluorescence strength is increasing a little from the first one to the last. In addition, the fluorescence of 58th seems the weakest in the four and the 54th obviously stronger than 58th. This result confirms ImageMe is sbetter again.
Therefore, ImageMe adapts the both situation while manual calculation has significant computation error when the total fluorescence strength is weak which means only gives reliable results under strong fluorescence.
