Line 1: | Line 1: | ||
− | <html | + | <html> |
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
+ | <head> | ||
+ | <link href="https://2021.igem.org/wiki/index.php?title=Template:Vilnius-Lithuania/styles/main&action=raw&ctype=text/css" rel="stylesheet" type="text/css" /> | ||
+ | <link href="https://2021.igem.org/wiki/index.php?title=Template:Vilnius-Lithuania/styles/fonts&action=raw&ctype=text/css" rel="stylesheet" type="text/css" /> | ||
+ | <link href="https://2021.igem.org/wiki/index.php?title=Template:Vilnius-Lithuania/styles/background&action=raw&ctype=text/css" rel="stylesheet" type="text/css" /> | ||
+ | <link href="https://2021.igem.org/wiki/index.php?title=Template:Vilnius-Lithuania/styles/contentpage-desktop&action=raw&ctype=text/css" rel="stylesheet" type="text/css" /> | ||
+ | <link href="https://2021.igem.org/wiki/index.php?title=Template:Vilnius-Lithuania/styles/contentpage-mobile&action=raw&ctype=text/css" rel="stylesheet" type="text/css" /> | ||
+ | <link href="https://2021.igem.org/wiki/index.php?title=Template:Vilnius-Lithuania/styles/sideindex-desktop&action=raw&ctype=text/css" rel="stylesheet" type="text/css" /> | ||
+ | <link href="https://2021.igem.org/wiki/index.php?title=Template:Vilnius-Lithuania/styles/sideindex-mobile&action=raw&ctype=text/css" rel="stylesheet" type="text/css" /> | ||
+ | <link href="https://2021.igem.org/wiki/index.php?title=Template:Vilnius-Lithuania/styles/navbar&action=raw&ctype=text/css" rel="stylesheet" type="text/css" /> | ||
+ | <link href="https://2021.igem.org/wiki/index.php?title=Template:Vilnius-Lithuania/styles/navmenu-desktop&action=raw&ctype=text/css" rel="stylesheet" type="text/css" /> | ||
+ | <link href="https://2021.igem.org/wiki/index.php?title=Template:Vilnius-Lithuania/styles/navmenu-mobile&action=raw&ctype=text/css" rel="stylesheet" type="text/css" /> | ||
+ | <link href="https://2021.igem.org/wiki/index.php?title=Template:Vilnius-Lithuania/styles/accesibility-menu-desktop&action=raw&ctype=text/css" rel="stylesheet" type="text/css" /> | ||
+ | <link href="https://2021.igem.org/wiki/index.php?title=Template:Vilnius-Lithuania/styles/accesibility-menu-mobile&action=raw&ctype=text/css" rel="stylesheet" type="text/css" /> | ||
+ | <link href="https://2021.igem.org/wiki/index.php?title=Template:Vilnius-Lithuania/styles/button-pill&action=raw&ctype=text/css" rel="stylesheet" type="text/css" /> | ||
+ | <link href="https://2021.igem.org/wiki/index.php?title=Template:Vilnius-Lithuania/styles/footer&action=raw&ctype=text/css" rel="stylesheet" type="text/css" /> | ||
+ | <link href="https://2021.igem.org/wiki/index.php?title=Template:Vilnius-Lithuania/styles/table&action=raw&ctype=text/css" rel="stylesheet" type="text/css" /> | ||
+ | <link href="https://2021.igem.org/wiki/index.php?title=Template:Vilnius-Lithuania/styles/bypass-block&action=raw&ctype=text/css" rel="stylesheet" type="text/css" /> | ||
+ | <script type="text/x-mathjax-config"> MathJax.Hub.Config({ TeX: { equationNumbers: { autoNumber: "AMS" } } }); </script> | ||
+ | <script src="https://2021.igem.org/common/MathJax-2.5-latest/MathJax.js?config=TeX-AMS-MML_HTMLorMML"></script> | ||
+ | </head> | ||
− | + | <body> | |
− | + | <div class="navbar-container"> | |
− | <div class=" | + | <nav class="navbar"> <a class="bypass-block-link visually-hidden visible-when-focused" href="#main-content">Skip to main content</a> |
− | <div> | + | <a class="navbar-logo" href="/Team:Vilnius-Lithuania"> <img alt="AmebyeLogo" src="https://static.igem.org/mediawiki/2021/b/b0/T--Vilnius-Lithuania--amebyeLogo.svg" /> AmeBye </a> |
− | + | <ul class="nav-menu"></ul> | |
− | + | </nav> | |
+ | <div class="progress-container"> | ||
+ | <div class="progress-bar"></div> | ||
+ | </div> | ||
+ | <div class="animation-container"> | ||
+ | <script> | ||
+ | document.querySelector('body').classList.add("content-hidden"); window.addEventListener("load", function() { setTimeout(function() { document.querySelector('body').classList.remove("content-hidden"); document.querySelector('.animation-container').classList.add("loaded"); document.querySelector('.animation-container').innerHTML=""; }, 200); }); | ||
+ | </script> | ||
+ | <div id="loading-animation"></div> | ||
+ | <div>Loading...</div> | ||
+ | <script src="https://2021.igem.org/wiki/index.php?title=Template:Vilnius-Lithuania/scripts/lottie&action=raw&ctype=text/javascript" type="text/javascript"></script> | ||
+ | <script src="https://2021.igem.org/wiki/index.php?title=Template:Vilnius-Lithuania/scripts/loading_animation&action=raw&ctype=text/javascript" type="text/javascript"></script> | ||
+ | <script> | ||
+ | var animation = bodymovin.loadAnimation({ container: document.getElementById('loading-animation'), animationData: loading_animation, renderer: 'svg', loop: true, autoplay: true, }); | ||
+ | </script> | ||
+ | </div> | ||
+ | </div> | ||
+ | <div class="background"> | ||
+ | <canvas id="background-canvas"></canvas> | ||
+ | <canvas id="canvas-transition"></canvas> | ||
+ | <div class="app-header"> | ||
+ | <div class="app-header-image-wrapper" id="img"> | ||
+ | <h1 id="title"> MEASUREMENT </h1> <img alt="Header" class="cover" src="https://static.igem.org/mediawiki/2021/1/1b/T--Vilnius-Lithuania--Measurement.JPG" /> </div> | ||
+ | </div> | ||
+ | <div class="app-content"> | ||
+ | <div class="app-content-text"> | ||
+ | <div class="content-page-container"> | ||
+ | <h3 class="index-headline"> Introduction </h3> | ||
+ | <p> To construct an efficient naringenin biosynthesis pathway in probiotic bacteria, selected promoters strength in our chassis had to be evaluated. This was done by measuring the fluorescence of bacteria containing different promoters | ||
+ | and sfGFP protein (<a href="http://parts.igem.org/Part:BBa_K1365020" target="”_blank”">BBa_K1365020</a>). What is more, during the same set of experiments, mRNA cyclization (<a href="http://parts.igem.org/Part:BBa_K3904217" target="”_blank”">BBa_K3904217</a>) | ||
+ | performance was tested in our system. </p> | ||
+ | <h3 class="index-headline"> Calibration </h3> | ||
+ | <p> Before the actual experiments, we calibrated our measurements with an iGEM measurement kit according to the <a href="https://www.protocols.io/view/calibration-protocol-plate-reader-fluorescence-cal-6zrhf56" target="”_blank”">standard protocol</a>. | ||
+ | This calibration allowed choosing the most appropriate gain settings in the plate reader and standardizing our results. </p> | ||
+ | <h3 class="index-headline">Experiment</h3> | ||
+ | <p>The protocol used for measurement of fluorescence and OD600.</p> | ||
+ | <p>Preparation of bacteria:</p> | ||
+ | <ol> | ||
+ | <li>Transform bacteria with plasmid containing your target genes.</li> | ||
+ | <li>Pick three colonies from each plate and inoculate in 5 mL of LB medium with the right antibiotic.</li> | ||
+ | <li>Make a 1:10 dilution of each overnight culture in LB medium with the right antibiotic.</li> | ||
+ | <li>Measure OD600 and write down the values of OD600.</li> | ||
+ | <li>Dilute the culture to 0.05 OD600 in an eppendorf tube. The sample is placed on ice.</li> | ||
+ | <li>The samples should be laid out according to the plate down below:</li> | ||
+ | </ol> | ||
+ | <div class="figure-container"> <img src="https://static.igem.org/mediawiki/2021/a/a6/T--Vilnius-Lithuania--fluorescence_protocol.png" /> | ||
+ | <div> <b>Fig. 1.</b> Fluorescence measurement plate </div> | ||
+ | </div> | ||
+ | <h3 class="index-headline"> Measurement </h3> | ||
+ | <p> At first, the fluorescence of <i>E. coli</i> Nissle 1917 bacteria containing pTRKH2+sfGFP and pTRKH2+loop+sfGFP plasmids with different promoters (Table 1) was measured. </p> | ||
+ | <div class="table-container"> | ||
+ | <div class="table-headline"> <b> Table. 1. </b> Characterized promoters list </div> | ||
+ | <table class="table table-bordered table-hover table-condensed"> | ||
+ | <thead> | ||
+ | <tr> | ||
+ | <th title="Field #1">Part</th> | ||
+ | <th title="Field #2">Name</th> | ||
+ | </tr> | ||
+ | </thead> | ||
+ | <tbody> | ||
+ | <tr> | ||
+ | <td><a href="http://parts.igem.org/Part:BBa_K1033225" target="”_blank”">BBa_K1033225</a></td> | ||
+ | <td>CP44 promoter</td> | ||
+ | </tr> | ||
+ | <tr> | ||
+ | <td><a href="http://parts.igem.org/Part:BBa_K1033222" target="”_blank”">BBa_K1033222</a></td> | ||
+ | <td>CP29 promoter</td> | ||
+ | </tr> | ||
+ | <tr> | ||
+ | <td><a href="http://parts.igem.org/Part:BBa_K1033220" target="”_blank”">BBa_K1033220</a></td> | ||
+ | <td>CP8 promoter</td> | ||
+ | </tr> | ||
+ | <tr> | ||
+ | <td><a href="http://http://parts.igem.org/Part:BBa_J23118" target="”_blank”">BBa_J23118</a></td> | ||
+ | <td>Anderson promoter</td> | ||
+ | </tr> | ||
+ | <tr> | ||
+ | <td><a href="http://http://parts.igem.org/Part:BBa_J23117" target="”_blank”">BBa_J23117</a></td> | ||
+ | <td>Anderson promoter</td> | ||
+ | </tr> | ||
+ | <tr> | ||
+ | <td><a href="http://http://parts.igem.org/Part:BBa_J23115" target="”_blank”">BBa_J23115</a></td> | ||
+ | <td>Anderson promoter</td> | ||
+ | </tr> | ||
+ | <tr> | ||
+ | <td><a href="http://http://parts.igem.org/Part:BBa_J23114" target="”_blank”">BBa_J23114</a></td> | ||
+ | <td>Anderson promoter</td> | ||
+ | </tr> | ||
+ | <tr> | ||
+ | <td><a href="http://http://parts.igem.org/Part:BBa_J23113" target="”_blank”">BBa_J23113</a></td> | ||
+ | <td>Anderson promoter</td> | ||
+ | </tr> | ||
+ | <tr> | ||
+ | <td><a href="http://http://parts.igem.org/Part:BBa_J23107" target="”_blank”">BBa_J23107</a></td> | ||
+ | <td>Anderson promoter</td> | ||
+ | </tr> | ||
+ | <tr> | ||
+ | <td><a href="http://http://parts.igem.org/Part:BBa_J23106" target="”_blank”">BBa_J23106</a></td> | ||
+ | <td>Anderson promoter</td> | ||
+ | </tr> | ||
+ | <tr> | ||
+ | <td><a href="http://http://parts.igem.org/Part:BBa_J23103" target="”_blank”">BBa_J23103</a></td> | ||
+ | <td>Anderson promoter</td> | ||
+ | </tr> | ||
+ | </tbody> | ||
+ | </table> | ||
+ | </div> | ||
+ | <p> Since mRNA cyclization system performance was not as expected, we did a deeper analysis of the reference article [1] and decided to repeat the experiment while using longer sequence protein TAL fused with sfGFP. What is more, these | ||
+ | measurements were conducted in <i>E. coli</i> DH5Α and at different temperatures of 37°C and 24°C, since during the construction of mRNA cyclization system, we noticed that it "activates" and begins to fluorescent after incubation | ||
+ | at the room temperature. </p> | ||
+ | <p> Each experiment had three controls: positive – fluorescent <i>E. coli</i> DH5Α with J23101 Anderson promoter, negative – non-fluorescent bacteria, contamination – media with an antibiotic. </p> | ||
+ | <p> Overall, the detailed descriptions of the procedure can be found on <a href="https://2021.igem.org/Team:Vilnius-Lithuania/Experiments" target="”_blank”">our Experiments page.</a> </p> | ||
+ | <h3 class="index-headline">Methodology</h3> | ||
+ | <h4>Functional data analysis</h4> | ||
+ | <p> Functional data analysis (FDA) can be applied more and more broadly every year. At a high level, FDA arises when at least one of the variables is viewed as a function. When discrete data is made continuous, we can employ FDA for multivariate | ||
+ | statistical analysis. </p> | ||
+ | <h4>Smoothing</h4> | ||
+ | <p> First, we define a set of functions - <b>basis functions</b> - that we will use to approximate discrete data. Secondly, we set up a matrix of coefficients to define the function as a linear combination of the basis functions. In our | ||
+ | case, we have employed <i>splines</i> or <i>monomial</i> basis which is $$x(t) = \sum_{k=1}^{K} c_k \phi_k(t) = c'\phi(t)$$ written as </p> | ||
+ | <div class="copyable-code-container"> | ||
+ | <div> databasis = create.monomial.basis(c(measure_beginning, measure_ending),number_of_basis_functions) | ||
+ | <br/> data.fd = smooth.basis(y=data.matrix, fdParobj= databasis) | ||
+ | <br/> </div> | ||
+ | <div class="btn" role="button" tabindex="0"> | ||
+ | <svg fill="none" viewbox="0 0 24 34" xmlns="http://www.w3.org/2000/svg"> | ||
+ | <path d="M20.8001 11.9722H10.9002C9.68516 11.9722 8.7002 13.3757 8.7002 15.1071V29.2145C8.7002 30.9459 9.68516 32.3495 10.9002 32.3495H20.8001C22.0151 32.3495 23 30.9459 23 29.2145V15.1071C23 13.3757 22.0151 11.9722 20.8001 11.9722Z" | ||
+ | stroke="#0D4B5E" stroke-linecap="round" stroke-linejoin="round" stroke-width="2"></path> | ||
+ | <path d="M4.29997 21.3773H3.19998C2.61651 21.3773 2.05693 21.047 1.64436 20.4591C1.23178 19.8712 1 19.0738 1 18.2423V4.13497C1 3.30352 1.23178 2.50613 1.64436 1.91821C2.05693 1.33029 2.61651 1 3.19998 1H13.0999C13.6833 1 14.2429 1.33029 14.6555 1.91821C15.0681 2.50613 15.2999 3.30352 15.2999 4.13497V5.70245" | ||
+ | stroke="#0D4B5E" stroke-linecap="round" stroke-linejoin="round" stroke-width="2"></path> | ||
+ | </svg> | ||
+ | </div> | ||
+ | </div> | ||
+ | <h4>Statistical analysis</h4> | ||
+ | <p> Next descriptive statistical tools can be used, for instance, we used mean and standard deviation to analyse the variability of the experiment. $$\overline{x}(t) = N^{-1}\sum_{i}x_i(t) \quad \text{and} \quad s(t) = (N-1)^{-1} \sum_{i} | ||
+ | \left[ x_i(t) - \overline{x}(t) \right]^2.$$ </p> | ||
+ | <p> <b>FDA</b> package has an immense number of built-in functions which can be found in the <a href="”https://cran.r-project.org/web/packages/fda/fda.pdf”" target="”_blank”">documentation.</a> </p> | ||
+ | <div class="copyable-code-container"> | ||
+ | <div> mean = mean.fd(data.fd$fd) | ||
+ | <br/> std = std.fd(data.fd$fd) | ||
+ | <br/> | ||
+ | <br/> # You can use basic mathematical operations to work with fd objects | ||
+ | <br/> two_std = mean + std * 2 | ||
+ | <br/> </div> | ||
+ | <div class="btn" role="button" tabindex="0"> | ||
+ | <svg fill="none" viewbox="0 0 24 34" xmlns="http://www.w3.org/2000/svg"> | ||
+ | <path d="M20.8001 11.9722H10.9002C9.68516 11.9722 8.7002 13.3757 8.7002 15.1071V29.2145C8.7002 30.9459 9.68516 32.3495 10.9002 32.3495H20.8001C22.0151 32.3495 23 30.9459 23 29.2145V15.1071C23 13.3757 22.0151 11.9722 20.8001 11.9722Z" | ||
+ | stroke="#0D4B5E" stroke-linecap="round" stroke-linejoin="round" stroke-width="2"></path> | ||
+ | <path d="M4.29997 21.3773H3.19998C2.61651 21.3773 2.05693 21.047 1.64436 20.4591C1.23178 19.8712 1 19.0738 1 18.2423V4.13497C1 3.30352 1.23178 2.50613 1.64436 1.91821C2.05693 1.33029 2.61651 1 3.19998 1H13.0999C13.6833 1 14.2429 1.33029 14.6555 1.91821C15.0681 2.50613 15.2999 3.30352 15.2999 4.13497V5.70245" | ||
+ | stroke="#0D4B5E" stroke-linecap="round" stroke-linejoin="round" stroke-width="2"></path> | ||
+ | </svg> | ||
+ | </div> | ||
+ | </div> | ||
+ | <h4>Landmark registration and phase-plane plot</h4> | ||
+ | <p> There are two kinds of variations: amplitude and phase. Amplitude variation is a variation when intensity of growth spurt varies and phase variation is a variation when growth spurt varies from process to process. To avoid amplitude | ||
+ | variation we use landmark registration because we are interested in the shift that appears in multiple observations. </p> | ||
+ | <p> Landmark registration aligns a moment of time that indicates when the common trend started by estimating a strictly increasing nonlinear transformation of time. It takes all the moments of time into a common time value. </p> | ||
+ | <p> Phase-plane plot is a simple acceleration against velocity plot which can be used to determine maximum potential or kinetic energy to analyse similarity between groups in the experiment. However, first of all we have to find the first | ||
+ | and second derivatives of the functional data. These derivatives are described as: </p> | ||
+ | <p> However, as we have seen in VapXD analysis, it is more useful for harmonic processes. </p> | ||
+ | <div class="copyable-code-container"> | ||
+ | <div> # Second parameter indicates a number of derivatives taken | ||
+ | <br/> Accel.fd = deriv.fd(data.fd$fd,2) | ||
+ | <br/> Accel_mean.fd= mean.fd(accel.fd) | ||
+ | <br/> | ||
+ | <br/> # Vector to register the same point of every function | ||
+ | <br/> registration_vector = rep(0,number_of_processes) | ||
+ | <br/> time = seq(measure_beginning, measure_ending, length = measure_ending) | ||
+ | <br/> | ||
+ | <br/> # Create a window to see the warping mutual point visually | ||
+ | <br/> par(mfrow = c(1,1)) | ||
+ | <br/> for (icase in 1:number_of_processes) { | ||
+ | <br/> accveci = predict(accel.fd[icase], time) | ||
+ | <br/> windows() | ||
+ | <br/> plot(time, accveci, "l", ylim = c(measure_ylim_lower,measure_ylim_upper), | ||
+ | <br/> xlab="Time, h", ylab = "OD 600", | ||
+ | <br/> main=paste("Case", icase)) | ||
+ | <br/> lines(c(measure_beginning, measure_ending), c(0,0), lty = 2) | ||
+ | <br/> | ||
+ | <br/> # How many locators you want, in other word how many mutual points observed processes have in common | ||
+ | <br/> registration_vector [icase] = locator(1) | ||
+ | <br/> } | ||
+ | <br/> | ||
+ | <br/> # Averaging registered values | ||
+ | <br/> registration_vector_means= mean(sapply(registration_vector , mean)) | ||
+ | <br/> | ||
+ | <br/> # Creating a basis for time transforming | ||
+ | <br/> wbasisLM = create.bspline.basis(c(measure_beginning, measure_ending), number_of_knots, number_of_functions, c(measure_beginning, registration_vector_means, measure_ending)) | ||
+ | <br/> Wfd = fd(matrix(c(0:3)), wbasis) | ||
+ | <br/> WfdPar = fdPar(WfdLM,1,10^(-12)) | ||
+ | <br/> | ||
+ | <br/> # Putting everything to landmarker function | ||
+ | <br/> regList = landmarkreg(accel.fd, sapply(registration_vector,mean), registration_vector_means, WfdPar, TRUE) | ||
+ | <br/> reg_accel.fd = regList$regfd | ||
+ | <br/> | ||
+ | <br/> # Find mean of already registered observations | ||
+ | <br/> accelmean.fd = mean.fd(reg_accel.fd) | ||
+ | <br/> warp.fd = regList$warpfd | ||
+ | <br/> W.fd = regList$Wfd | ||
+ | <br/> | ||
+ | <br/> </div> | ||
+ | <div class="btn" role="button" tabindex="0"> | ||
+ | <svg fill="none" viewbox="0 0 24 34" xmlns="http://www.w3.org/2000/svg"> | ||
+ | <path d="M20.8001 11.9722H10.9002C9.68516 11.9722 8.7002 13.3757 8.7002 15.1071V29.2145C8.7002 30.9459 9.68516 32.3495 10.9002 32.3495H20.8001C22.0151 32.3495 23 30.9459 23 29.2145V15.1071C23 13.3757 22.0151 11.9722 20.8001 11.9722Z" | ||
+ | stroke="#0D4B5E" stroke-linecap="round" stroke-linejoin="round" stroke-width="2"></path> | ||
+ | <path d="M4.29997 21.3773H3.19998C2.61651 21.3773 2.05693 21.047 1.64436 20.4591C1.23178 19.8712 1 19.0738 1 18.2423V4.13497C1 3.30352 1.23178 2.50613 1.64436 1.91821C2.05693 1.33029 2.61651 1 3.19998 1H13.0999C13.6833 1 14.2429 1.33029 14.6555 1.91821C15.0681 2.50613 15.2999 3.30352 15.2999 4.13497V5.70245" | ||
+ | stroke="#0D4B5E" stroke-linecap="round" stroke-linejoin="round" stroke-width="2"></path> | ||
+ | </svg> | ||
+ | </div> | ||
+ | </div> | ||
+ | <h4>Linear regression to compare genome and plasmid cases</h4> | ||
+ | <p> We have used a simple linear regression method to predict OD 600 level for plasmid and genome cases given an explanatory variable x in logarithmic base (log(x)) since the independent variable varies in a wide range from 2E-1 to 2E-6. | ||
+ | Plasmid and genome models achieved high R scores 91.1 percent, 91.8 percent respectively which indicates that models explain more than 91 precent of variability. </p> | ||
+ | <div class="figure-container center"> <img alt="" src="https://static.igem.org/mediawiki/2021/b/bb/T--Vilnius-Lithuania--linear-regression.png" /> | ||
+ | <div> <b> Fig. 2. </b> Linear regression plot</div> | ||
+ | </div> | ||
+ | <div class="copyable-code-container"> | ||
+ | <div> # Create linear regression for y-plasmid, genome data, x - initial OD 600 | ||
+ | <br/> model_plasmid = lm(y[,1]~I(log(x))) | ||
+ | <br/> model_genoma = lm(y[,2]~I(log(x))) | ||
+ | <br/> | ||
+ | <br/> plot(y = y[,1], x = log(x), xlab = "log(x)", ylab = "OD 600", main = "Linear regression on plasmid data", lwd = 4, col = cl[3]) | ||
+ | <br/> abline(model_plasmid$coefficients, lwd = 3) # Call model coefficients to draw the model linear regression | ||
+ | <br/> | ||
+ | <br/> plot(y = y[,2], x = log(x), xlab = "log(x)", ylab = "OD 600", main = "Linear regression on genome data", lwd = 4, col = cl[3]) | ||
+ | <br/> abline(model_genoma$coefficients, lwd = 3) | ||
+ | <br/> | ||
+ | <br/> # Genoma linear regression oefficient R = 91.1 | ||
+ | <br/> summary(model_genoma) | ||
+ | <br/> | ||
+ | <br/> # Plasmid linear regression coefficient R = 91.9 % | ||
+ | <br/> summary(model_plasmid) | ||
+ | <br/> | ||
+ | <br/> # Creating x values vector for model to make predictions for | ||
+ | <br/> genoma_value <- data.frame(x = seq(0.01, 0.06, by=0.0025)) | ||
+ | <br/> plasmid_value <- data.frame(x = seq(0.01, 0.06, by=0.0025)) | ||
+ | <br/> | ||
+ | <br/> # Predicts post-experiment OD 600 value for a given x value pre-experiment OD 600 | ||
+ | <br/> genoma_pred = predict(model_genoma,genoma_value) | ||
+ | <br/> plasmid_pred = predict(model_plasmid,plasmid_value) | ||
+ | <br/> | ||
+ | <br/> print(genoma_pred) # Output is y prediction for given x, find a closest to t value to investigate x that satisfy your investigation | ||
+ | <br/> print(plasmid_pred) | ||
+ | <br/> | ||
+ | <br/> </div> | ||
+ | <div class="btn" role="button" tabindex="0"> | ||
+ | <svg fill="none" viewbox="0 0 24 34" xmlns="http://www.w3.org/2000/svg"> | ||
+ | <path d="M20.8001 11.9722H10.9002C9.68516 11.9722 8.7002 13.3757 8.7002 15.1071V29.2145C8.7002 30.9459 9.68516 32.3495 10.9002 32.3495H20.8001C22.0151 32.3495 23 30.9459 23 29.2145V15.1071C23 13.3757 22.0151 11.9722 20.8001 11.9722Z" | ||
+ | stroke="#0D4B5E" stroke-linecap="round" stroke-linejoin="round" stroke-width="2"></path> | ||
+ | <path d="M4.29997 21.3773H3.19998C2.61651 21.3773 2.05693 21.047 1.64436 20.4591C1.23178 19.8712 1 19.0738 1 18.2423V4.13497C1 3.30352 1.23178 2.50613 1.64436 1.91821C2.05693 1.33029 2.61651 1 3.19998 1H13.0999C13.6833 1 14.2429 1.33029 14.6555 1.91821C15.0681 2.50613 15.2999 3.30352 15.2999 4.13497V5.70245" | ||
+ | stroke="#0D4B5E" stroke-linecap="round" stroke-linejoin="round" stroke-width="2"></path> | ||
+ | </svg> | ||
+ | </div> | ||
+ | </div> | ||
+ | <h3 class="index-headline"> Results </h3> | ||
+ | <h4>pTRKH2+sfGFP and pTRKH2+loop+sfGFP</h4> | ||
+ | <p> After the first set of measurements with different promoters, it was seen that the mRNA cyclization system was not working as it was expected. For this reason, only results obtained from pTRKH2+sfGFP with different promoters were analyzed. | ||
+ | They allowed us to conclude promoters’ strength. In the Fig. 3, the fluorescence and OD600 ratio were compared. The results showed that p-slpA (<a href="http://parts.igem.org/Part:BBa_K3904228" target="”_blank”">BBa_K3904712</a>) | ||
+ | was the strongest, and further, it was selected as a reference point for other promoters evaluation. The trio of the most effective promoters was p-slpA, CP44 (<a href="http://parts.igem.org/Part:BBa_K1033225" target="”_blank”">BBa_K1033225</a>) | ||
+ | and CP8 (<a href="http://parts.igem.org/Part:BBa_K1033220" target="”_blank”">BBa_K1033220</a>). </p> | ||
+ | <div class="figure-container center"> <img alt="" src="https://static.igem.org/mediawiki/2021/e/e7/T--Vilnius-Lithuania--promoter-old-1.png" /> | ||
+ | <div> <b> Fig. 3. </b> Fluorescence and OD600 ratio comparison of different promoters</div> | ||
+ | </div> | ||
+ | <h4>pTRKH2+TAL-sfGFP and pTRKH2+loop+TAL-sfGFP</h4> | ||
+ | <p> In the Fig. 4, one can see that there is no significant difference in the mRNA cyclization performance in 37°C and 24°C because the shift of the curves is only affected by the lower OD600 values in 24°C. The bottom graph illustrates | ||
+ | averaged data at 37°C. Interestingly, it was noticed that mRNA cyclization does not allow the accumulation of the sfGFP protein. After some time, the system reaches equilibrium, and the fluorescence/OD600 ratio stabilizes. </p> | ||
+ | <div class="figure-container center"> <img alt="" src="https://static.igem.org/mediawiki/2021/c/cf/T--Vilnius-Lithuania--promoter-new-1.png" /> | ||
+ | <div> <b> Fig. 4. </b> mRNA cyclization system performance</div> | ||
+ | </div> | ||
+ | <h4>Analysis of the cyclization system’s stability</h4> | ||
+ | <p> Fig. 5 represents stability of the processes without and with a cyclization system at 24°C </p> | ||
+ | <p> In the right graph we observe: relying on FDA methodology we were able to obtain a confidence interval (CI) of 95 percent, which is superbly narrow; it indicates the non-intuitive stability of the loop system. </p> | ||
+ | <p> On the other hand, the process in the left graph is observed as barely stable - it obtains a surprisingly wide CI, which has got two groups of the processes that localize around different CI boundaries. The first group (2 standard | ||
+ | deviations) provides an output of bacteria colony growth and the second group (-2 standard deviations) gives an opposite outcome. </p> | ||
+ | <div class="figure-container center"> <img alt="" src="https://static.igem.org/mediawiki/2021/c/c1/T--Vilnius-Lithuania--promoter-new-2.png" /> | ||
+ | <div> <b> Fig. 5. </b> Stability comparison graphs of system without (left) and with (right) loop</div> | ||
+ | </div> | ||
+ | <div class="figure-container center"> <img alt="" src="https://static.igem.org/mediawiki/2021/2/22/T--Vilnius-Lithuania--promoters_no_loop_registration.png" /> | ||
+ | <div> <b> Fig. 6. </b> Landmark registration of promoters (BBa1033225, BBa1033220, BBa1033222) </div> | ||
+ | </div> | ||
+ | <h3 class="index-headline"> Conclusions </h3> | ||
+ | <p> All in all, promoters' strength was successfully evaluated, and following these data, promoters were selected for the naringenin metabolic pathway. On the other hand, one of our chosen protein synthesis enhancing mechanisms - mRNA | ||
+ | cyclization - appeared not to work as expected and described in the Yang et al. article [1]. In our case, it did not increase protein production. Interestingly, protein quantity stabilization was noticed as a novel property of | ||
+ | this system and may be further tested. </p> | ||
+ | </div> | ||
+ | <div class="references-wrapper"> | ||
+ | <div class="breaker"> </div> | ||
+ | <h2> References </h2> | ||
+ | <div class="references-container"> | ||
+ | <div class="number">1.</div> | ||
+ | <div>Yang, J., Han, Y. H., Im, J., & Seo, S. W. (2021). Synthetic protein quality control to enhance full-length translation in bacteria. Nature Chemical Biology, 17(4), 421–427. <a href="https://doi.org/10.1038/s41589-021-00736-3" | ||
+ | target="_blank">To the article.</a></div> | ||
+ | </div> | ||
+ | </div> | ||
</div> | </div> | ||
− | <div | + | <div class="index-container"> |
− | + | <div class="index-header"> </div> | |
− | + | <div class="index-content"> </div> | |
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | < | + | |
− | + | ||
</div> | </div> | ||
− | + | </div> | |
− | + | <footer> | |
− | + | <div class="logo-igem"> | |
− | + | <object data="https://static.igem.org/mediawiki/2021/f/ff/T--Vilnius-Lithuania--iGEM-2021.svg"> </object> | |
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
</div> | </div> | ||
− | <div | + | <div class="social-container"> |
− | + | <div>FOLLOW US</div> | |
− | + | <div> | |
− | + | <a class="placeholder-social-icon" href="https://www.facebook.com/VilniusiGEM"> <img src="https://static.igem.org/mediawiki/2021/3/36/T--Vilnius-Lithuania--facebook.svg" /> </a> | |
− | + | <a class="placeholder-social-icon" href="https://www.instagram.com/igem_vilnius/"> <img src="https://static.igem.org/mediawiki/2021/6/64/T--Vilnius-Lithuania--instagram.svg" /> </a> | |
− | + | <a class="placeholder-social-icon" href="https://www.linkedin.com/company/vilnius-igem/"> <img src="https://static.igem.org/mediawiki/2021/e/e7/T--Vilnius-Lithuania--linkedin.svg" /> </a> | |
− | + | </div> | |
− | + | ||
− | + | ||
</div> | </div> | ||
− | + | <div class="mail-container"> | |
− | + | <div>CONTACT US</div> <a href="mailto:info@vilniusigem.lt">info@vilniusigem.lt</a> </div> | |
− | + | <div class="grid-sponsors"> | |
− | + | <div> | |
− | + | <div> | |
− | + | <object data="https://static.igem.org/mediawiki/2021/d/dc/T--Vilnius-Lithuania--VU.svg"> </object> | |
− | + | </div> | |
− | + | <div> | |
− | + | <object data="https://static.igem.org/mediawiki/2021/b/bf/T--Vilnius-Lithuania--Termofisher.svg"> </object> | |
− | + | </div> | |
− | + | <div> | |
− | + | <object data="https://static.igem.org/mediawiki/2021/1/10/T--Vilnius-Lithuania--CityOfVilnius.svg"> </object> | |
− | + | </div> | |
− | + | </div> | |
− | + | <div> | |
− | + | <div> | |
− | + | <object data="https://static.igem.org/mediawiki/2021/b/bb/T--Vilnius-Lithuania--GMC.svg"> </object> | |
− | + | </div> | |
− | + | <div> | |
− | + | <object data="https://static.igem.org/mediawiki/2021/9/98/T--Vilnius-Lithuania--Nanodiagnostika.svg"> </object> | |
− | + | </div> | |
− | + | <div> | |
− | + | <object data="https://static.igem.org/mediawiki/2021/1/16/T--Vilnius-Lithuania--Telesoftas.svg"> </object> | |
− | + | </div> | |
− | + | <div> | |
− | + | <object data="https://static.igem.org/mediawiki/2021/7/70/T--Vilnius-Lithuania--Kopicentras.svg"> </object> | |
− | + | </div> | |
− | + | </div> | |
− | + | <div> | |
− | + | <div> | |
− | + | <object data="https://static.igem.org/mediawiki/2021/d/df/T--Vilnius-Lithuania--SnapGene.svg"> </object> | |
− | + | </div> | |
− | + | <div> | |
− | + | <object data="https://static.igem.org/mediawiki/2021/3/3e/T--Vilnius-Lithuania--Laborama.svg"> </object> | |
− | + | </div> | |
− | + | <div> | |
− | + | <object data="https://static.igem.org/mediawiki/2021/c/c1/T--Vilnius-Lithuania--Biotecha.svg"> </object> | |
− | + | </div> | |
− | + | <div> | |
− | + | <object data="https://static.igem.org/mediawiki/2021/4/4f/T--Vilnius-Lithuania--Grida.svg"> </object> | |
− | + | </div> | |
− | + | </div> | |
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
</div> | </div> | ||
− | + | </footer> | |
− | + | </div> | |
− | + | <script src="https://2021.igem.org/wiki/index.php?title=Template:Vilnius-Lithuania/scripts/navigationTabs&action=raw&ctype=text/javascript" type="text/javascript"></script> | |
− | + | <script src="https://2021.igem.org/wiki/index.php?title=Template:Vilnius-Lithuania/scripts/background&action=raw&ctype=text/javascript" type="text/javascript"></script> | |
− | + | <script src="https://2021.igem.org/wiki/index.php?title=Template:Vilnius-Lithuania/scripts/contentpage&action=raw&ctype=text/javascript" type="text/javascript"></script> | |
− | + | <script src="https://2021.igem.org/wiki/index.php?title=Template:Vilnius-Lithuania/scripts/navbar&action=raw&ctype=text/javascript" type="text/javascript"></script> | |
− | + | <script> | |
− | + | contentPage( "Sections", true, 300, ) | |
− | + | </script> | |
− | + | <script src="https://2021.igem.org/wiki/index.php?title=Template:Vilnius-Lithuania/scripts/GlslCanvas&action=raw&ctype=text/javascript" type="text/javascript"></script> | |
− | + | <script src="https://2021.igem.org/wiki/index.php?title=Template:Vilnius-Lithuania/scripts/backgroundTransition&action=raw&ctype=text/javascript" type="text/javascript"></script> | |
− | + | </body> | |
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | + | ||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
− | |||
</html> | </html> |
Revision as of 00:05, 22 October 2021
MEASUREMENT
Introduction
To construct an efficient naringenin biosynthesis pathway in probiotic bacteria, selected promoters strength in our chassis had to be evaluated. This was done by measuring the fluorescence of bacteria containing different promoters and sfGFP protein (BBa_K1365020). What is more, during the same set of experiments, mRNA cyclization (BBa_K3904217) performance was tested in our system.
Calibration
Before the actual experiments, we calibrated our measurements with an iGEM measurement kit according to the standard protocol. This calibration allowed choosing the most appropriate gain settings in the plate reader and standardizing our results.
Experiment
The protocol used for measurement of fluorescence and OD600.
Preparation of bacteria:
- Transform bacteria with plasmid containing your target genes.
- Pick three colonies from each plate and inoculate in 5 mL of LB medium with the right antibiotic.
- Make a 1:10 dilution of each overnight culture in LB medium with the right antibiotic.
- Measure OD600 and write down the values of OD600.
- Dilute the culture to 0.05 OD600 in an eppendorf tube. The sample is placed on ice.
- The samples should be laid out according to the plate down below:
Measurement
At first, the fluorescence of E. coli Nissle 1917 bacteria containing pTRKH2+sfGFP and pTRKH2+loop+sfGFP plasmids with different promoters (Table 1) was measured.
Part | Name |
---|---|
BBa_K1033225 | CP44 promoter |
BBa_K1033222 | CP29 promoter |
BBa_K1033220 | CP8 promoter |
BBa_J23118 | Anderson promoter |
BBa_J23117 | Anderson promoter |
BBa_J23115 | Anderson promoter |
BBa_J23114 | Anderson promoter |
BBa_J23113 | Anderson promoter |
BBa_J23107 | Anderson promoter |
BBa_J23106 | Anderson promoter |
BBa_J23103 | Anderson promoter |
Since mRNA cyclization system performance was not as expected, we did a deeper analysis of the reference article [1] and decided to repeat the experiment while using longer sequence protein TAL fused with sfGFP. What is more, these measurements were conducted in E. coli DH5Α and at different temperatures of 37°C and 24°C, since during the construction of mRNA cyclization system, we noticed that it "activates" and begins to fluorescent after incubation at the room temperature.
Each experiment had three controls: positive – fluorescent E. coli DH5Α with J23101 Anderson promoter, negative – non-fluorescent bacteria, contamination – media with an antibiotic.
Overall, the detailed descriptions of the procedure can be found on our Experiments page.
Methodology
Functional data analysis
Functional data analysis (FDA) can be applied more and more broadly every year. At a high level, FDA arises when at least one of the variables is viewed as a function. When discrete data is made continuous, we can employ FDA for multivariate statistical analysis.
Smoothing
First, we define a set of functions - basis functions - that we will use to approximate discrete data. Secondly, we set up a matrix of coefficients to define the function as a linear combination of the basis functions. In our case, we have employed splines or monomial basis which is $$x(t) = \sum_{k=1}^{K} c_k \phi_k(t) = c'\phi(t)$$ written as
data.fd = smooth.basis(y=data.matrix, fdParobj= databasis)
Statistical analysis
Next descriptive statistical tools can be used, for instance, we used mean and standard deviation to analyse the variability of the experiment. $$\overline{x}(t) = N^{-1}\sum_{i}x_i(t) \quad \text{and} \quad s(t) = (N-1)^{-1} \sum_{i} \left[ x_i(t) - \overline{x}(t) \right]^2.$$
FDA package has an immense number of built-in functions which can be found in the documentation.
std = std.fd(data.fd$fd)
# You can use basic mathematical operations to work with fd objects
two_std = mean + std * 2
Landmark registration and phase-plane plot
There are two kinds of variations: amplitude and phase. Amplitude variation is a variation when intensity of growth spurt varies and phase variation is a variation when growth spurt varies from process to process. To avoid amplitude variation we use landmark registration because we are interested in the shift that appears in multiple observations.
Landmark registration aligns a moment of time that indicates when the common trend started by estimating a strictly increasing nonlinear transformation of time. It takes all the moments of time into a common time value.
Phase-plane plot is a simple acceleration against velocity plot which can be used to determine maximum potential or kinetic energy to analyse similarity between groups in the experiment. However, first of all we have to find the first and second derivatives of the functional data. These derivatives are described as:
However, as we have seen in VapXD analysis, it is more useful for harmonic processes.
Accel.fd = deriv.fd(data.fd$fd,2)
Accel_mean.fd= mean.fd(accel.fd)
# Vector to register the same point of every function
registration_vector = rep(0,number_of_processes)
time = seq(measure_beginning, measure_ending, length = measure_ending)
# Create a window to see the warping mutual point visually
par(mfrow = c(1,1))
for (icase in 1:number_of_processes) {
accveci = predict(accel.fd[icase], time)
windows()
plot(time, accveci, "l", ylim = c(measure_ylim_lower,measure_ylim_upper),
xlab="Time, h", ylab = "OD 600",
main=paste("Case", icase))
lines(c(measure_beginning, measure_ending), c(0,0), lty = 2)
# How many locators you want, in other word how many mutual points observed processes have in common
registration_vector [icase] = locator(1)
}
# Averaging registered values
registration_vector_means= mean(sapply(registration_vector , mean))
# Creating a basis for time transforming
wbasisLM = create.bspline.basis(c(measure_beginning, measure_ending), number_of_knots, number_of_functions, c(measure_beginning, registration_vector_means, measure_ending))
Wfd = fd(matrix(c(0:3)), wbasis)
WfdPar = fdPar(WfdLM,1,10^(-12))
# Putting everything to landmarker function
regList = landmarkreg(accel.fd, sapply(registration_vector,mean), registration_vector_means, WfdPar, TRUE)
reg_accel.fd = regList$regfd
# Find mean of already registered observations
accelmean.fd = mean.fd(reg_accel.fd)
warp.fd = regList$warpfd
W.fd = regList$Wfd
Linear regression to compare genome and plasmid cases
We have used a simple linear regression method to predict OD 600 level for plasmid and genome cases given an explanatory variable x in logarithmic base (log(x)) since the independent variable varies in a wide range from 2E-1 to 2E-6. Plasmid and genome models achieved high R scores 91.1 percent, 91.8 percent respectively which indicates that models explain more than 91 precent of variability.
model_plasmid = lm(y[,1]~I(log(x)))
model_genoma = lm(y[,2]~I(log(x)))
plot(y = y[,1], x = log(x), xlab = "log(x)", ylab = "OD 600", main = "Linear regression on plasmid data", lwd = 4, col = cl[3])
abline(model_plasmid$coefficients, lwd = 3) # Call model coefficients to draw the model linear regression
plot(y = y[,2], x = log(x), xlab = "log(x)", ylab = "OD 600", main = "Linear regression on genome data", lwd = 4, col = cl[3])
abline(model_genoma$coefficients, lwd = 3)
# Genoma linear regression oefficient R = 91.1
summary(model_genoma)
# Plasmid linear regression coefficient R = 91.9 %
summary(model_plasmid)
# Creating x values vector for model to make predictions for
genoma_value <- data.frame(x = seq(0.01, 0.06, by=0.0025))
plasmid_value <- data.frame(x = seq(0.01, 0.06, by=0.0025))
# Predicts post-experiment OD 600 value for a given x value pre-experiment OD 600
genoma_pred = predict(model_genoma,genoma_value)
plasmid_pred = predict(model_plasmid,plasmid_value)
print(genoma_pred) # Output is y prediction for given x, find a closest to t value to investigate x that satisfy your investigation
print(plasmid_pred)
Results
pTRKH2+sfGFP and pTRKH2+loop+sfGFP
After the first set of measurements with different promoters, it was seen that the mRNA cyclization system was not working as it was expected. For this reason, only results obtained from pTRKH2+sfGFP with different promoters were analyzed. They allowed us to conclude promoters’ strength. In the Fig. 3, the fluorescence and OD600 ratio were compared. The results showed that p-slpA (BBa_K3904712) was the strongest, and further, it was selected as a reference point for other promoters evaluation. The trio of the most effective promoters was p-slpA, CP44 (BBa_K1033225) and CP8 (BBa_K1033220).
pTRKH2+TAL-sfGFP and pTRKH2+loop+TAL-sfGFP
In the Fig. 4, one can see that there is no significant difference in the mRNA cyclization performance in 37°C and 24°C because the shift of the curves is only affected by the lower OD600 values in 24°C. The bottom graph illustrates averaged data at 37°C. Interestingly, it was noticed that mRNA cyclization does not allow the accumulation of the sfGFP protein. After some time, the system reaches equilibrium, and the fluorescence/OD600 ratio stabilizes.
Analysis of the cyclization system’s stability
Fig. 5 represents stability of the processes without and with a cyclization system at 24°C
In the right graph we observe: relying on FDA methodology we were able to obtain a confidence interval (CI) of 95 percent, which is superbly narrow; it indicates the non-intuitive stability of the loop system.
On the other hand, the process in the left graph is observed as barely stable - it obtains a surprisingly wide CI, which has got two groups of the processes that localize around different CI boundaries. The first group (2 standard deviations) provides an output of bacteria colony growth and the second group (-2 standard deviations) gives an opposite outcome.
Conclusions
All in all, promoters' strength was successfully evaluated, and following these data, promoters were selected for the naringenin metabolic pathway. On the other hand, one of our chosen protein synthesis enhancing mechanisms - mRNA cyclization - appeared not to work as expected and described in the Yang et al. article [1]. In our case, it did not increase protein production. Interestingly, protein quantity stabilization was noticed as a novel property of this system and may be further tested.