Team:Cambridge-JIC/Informatics

From 2014.igem.org

(Difference between revisions)
 
(23 intermediate revisions not shown)
Line 1: Line 1:
-
{{:Team:Cambridge-JIC/Templates/header}}
+
{{:Team:Cambridge-JIC/Templates/header_prototype3}}
<html>
<html>
Line 12: Line 12:
<body>
<body>
 +
<div align="center"><a href="https://2014.igem.org/wiki/index.php?title=Team:Cambridge-JIC/Informatics&action=edit">Edit this page</a></div>
<h2>Informatics</h2>
<h2>Informatics</h2>
-
To introduce a novel chassis and make full advantage of its properties, we want to computationally analyse its genome in order to:
+
 
 +
To introduce and facilitate future use of the novel chassis <Em>Marchantia polymorpha</Em>, we computationally analysed its genome to:
<ul>
<ul>
-
<li>optimise codon usage in our registry parts and to facilitate future synthetic biology work on Marchantia;</li>
+
<li>find out the most efficient codon usage in order to optimise our Marchantia specific registry parts and facilitate all future synthetic biology work on Marchantia;</li>
-
<li>identify and characterise promoters, in particular looking for strong, inducible, tissue-specific or early development stage promoters.</li>
+
<li>submit a small library of Marchantia promoters to the iGEM registry, in particular looking for those which are strong, inducible, tissue-specific or expressed in an early development stage. </li>
</ul>
</ul>
-
<h3> Codon usage optimisation </h3>
+
<h3 id="Codon-optimisation">Codon usage optimisation </h3>
 +
 
<p>
<p>
-
Our start point was the Marchantia genome and the mRNA transcriptome predicted with the Geneious software (<a href = "http://www.geneious.com/">http://www.geneious.com/</a>). The data was given to us from Jim's lab by Bernardo. 99 000 ORFs were predicted, which seems too large to be realistic. Half of these were only 100 aa long. We set the threshold for candidate genes amongst these at 300 aa, obtaining the expected normal distribution of lengths.
+
Our start point was the Marchantia genome and the mRNA transcriptome predicted with Geneious software (<a href = "http://www.geneious.com/">http://www.geneious.com/</a>). The data was given to us from Jim's lab by Bernardo. 99 000 ORFs were predicted, which seems too large to be a realistic number of genes expressed by this small liverwort. Half of these were only 100 amino acids long. We set the threshold for candidate genes among these at 300 amino acids, obtaining the expected normal distribution of lengths.
</p>
</p>
<p>
<p>
-
Using the BLAST software (<a href="http://blast.ncbi.nlm.nih.gov/Blast.cgi">http://blast.ncbi.nlm.nih.gov/Blast.cgi</a>), we compared the proteins coded by these candidate genes with the proteins present in Arabidopsis, given on Araport (<a href="www.araport.org">www.araport.org</a>). Some of the sequences showed incomplete matches, indicated that our predicted ORFs should be regarded with some vigilance. The longest sequence showed a 40% match, a small number as expected.
+
Using BLAST software (<a href="http://blast.ncbi.nlm.nih.gov/Blast.cgi">http://blast.ncbi.nlm.nih.gov/Blast.cgi</a>), we compared the proteins coded by these candidate genes to the proteins present in <Em>Arabidopsis thaliana</Em>, given on Araport (<a href="www.araport.org">www.araport.org</a>). Some of the sequences showed incomplete matches, indicated that our predicted ORFs should be regarded with some vigilance. The longest sequence showed a 40% match, a small number as expected.
</p>
</p>
Line 32: Line 35:
</p>
</p>
-
<table class="reference" style="width:300px">
+
<img src="https://static.igem.org/mediawiki/2014/a/a9/Cambridge-JIC_Codon_table.png" width = 550> </img>
-
<tr>
+
 
-
<th> codon </th>
+
<p>A .xls version of this table with relevant filters can be found <a href="https://static.igem.org/mediawiki/2014/4/42/Marchantia_codon_table_analysis.xls">here</a>.</p>
-
<th> amino acid </th>
+
 
-
<th> per thousand </th>
+
<h2 id="Promoter">Promoter identification</h2>
-
<th> frequency </th>
+
<p>To identify potential promoter regions in Marchantia's genome, we looked for gene sequences from other species, mainly <i>Arabidopsis thaliana</i> and <i>Physcomitrella patens</i>. We were particularly interested in genes that were:
-
</tr>
+
<ul>
-
<tr>
+
<li>Nitrate inducible</li>
-
<td>aaa</td>
+
<li>Sulphate inducible</li>
-
<td>K</td>
+
<li>Phosphate inducible</li>
-
<td>15.5903</td>
+
<li>Circadian rythm driven</li>
-
<td>0.386139</td>
+
<li>Metabolism related</li>
-
</tr>
+
<li>Development related</li>
-
<tr>
+
<li>Light responsive</li>
-
<td>aag</td>
+
<li>Development</li>
-
<td>K</td>
+
<li>Standard templates for proteins</li>
-
<td>24.7846</td>
+
<li>Standard templates for genes</li>
-
<td>0.613861</td>
+
</ul>
-
</tr>
+
</p>
-
<tr>
+
<p>Reading recent research papers in search for these genes and looking for their sequences on GenBank and ThaleMine, we constructed a list of protein sequences to compare to the Marchantia predicted scaffolds. (For genes where only the nucleotide sequence was available, instead of the protein sequence, a C++ code was written to perform the translation). We ran a tblastn search in BLAST for the best matches between our candidate proteins and the Marchantia scaffolds. Then, for the best matches (~60% and above, with some judging by eye), we annotated the 2kb upstream of the start codon as a potential promoter region.</p>
-
<td>aac</td>
+
<img src="https://static.igem.org/mediawiki/2014/5/5d/Cambridge_JIC_Blast_example.png">
-
<td>N</td>
+
<p>Example of a BLAST hit, matching an inducible nitrate transporter sequence to a Marchantia gene</p>
-
<td>12.1702</td>
+
<p>We identified 30 candidate promoters this way, that we are planning to screen by inserting in a construct driving the yellow fluorescent protein Venus. For each promoter, we will make a construct with and one without amplification by GAL4 and GAL4 UAS, to evaluate the promoter strength and get around any leakages due to the use of GAL4. </p>
-
<td>0.485816</td>
+
 
-
</tr>
+
<!--Edit this page link-->
-
<tr>
+
<div align="center"><a href="https://2014.igem.org/wiki/index.php?title=Team:Cambridge-JIC/Informatics&action=edit">Edit this page</a>
-
<td>aat</td>
+
</div><br>
-
<td>N</td>
+
<!--End of Edit this page link-->
-
<td>12.8809</td>
+
-
<td>0.514184</td>
+
-
</tr>
+
-
<tr>
+
-
<td>aga</td>
+
-
<td>R</td>
+
-
<td>21.4977</td>
+
-
<td>0.204651</td>
+
-
</tr>
+
-
<tr>
+
-
<td>agg</td>
+
-
<td>R</td>
+
-
<td>21.7198</td>
+
-
<td>0.206765</td>
+
-
</tr>
+
-
<tr>
+
-
<td>agc</td>
+
-
<td>S</td>
+
-
<td>22.2528</td>
+
-
<td>0.211839</td>
+
-
</tr>
+
-
<tr>
+
-
<td>agt</td>
+
-
<td>S</td>
+
-
<td>11.6816</td>
+
-
<td>0.111205</td>
+
-
</tr>
+
-
<tr>
+
-
<td>aca</td>
+
-
<td>T</td>
+
-
<td>16.2121</td>
+
-
<td>0.249317</td>
+
-
</tr>
+
-
<tr>
+
-
<td>acg</td>
+
-
<td>T</td>
+
-
<td>17.4114</td>
+
-
<td>0.26776</td>
+
-
</tr>
+
-
<tr>
+
-
<td>acc</td>
+
-
<td>T</td>
+
-
<td>12.4811</td>
+
-
<td>0.19194</td>
+
-
</tr>
+
-
<tr>
+
-
<td>act</td>
+
-
<td>T</td>
+
-
<td>10.5268</td>
+
-
<td>0.161885</td>
+
-
</tr>
+
-
<tr>
+
-
<td>ata</td>
+
-
<td>I</td>
+
-
<td>8.83894</td>
+
-
<td>0.253827</td>
+
-
</tr>
+
-
<tr>
+
-
<td>atg</td>
+
-
<td>M</td>
+
-
<td>20.7426</td>
+
-
<td>1</td>
+
-
</tr>
+
-
<tr>
+
-
<td>atc</td>
+
-
<td>I</td>
+
-
<td>13.858</td>
+
-
<td>0.304094</td>
+
-
</tr>
+
-
<tr>
+
-
<td>att</td>
+
-
<td>I</td>
+
-
<td>12.1258</td>
+
-
<td>0.284969</td>
+
-
</tr>
+
-
<tr>
+
-
<td>gaa</td>
+
-
<td>E</td>
+
-
<td>16.4342</td>
+
-
<td>0.4474</td>
+
-
</tr>
+
-
<tr>
+
-
<td>gag</td>
+
-
<td>E</td>
+
-
<td>20.2985</td>
+
-
<td>0.5526</td>
+
-
</tr>
+
-
<tr>
+
-
<td>gac</td>
+
-
<td>D</td>
+
-
<td>13.4139</td>
+
-
<td>0.521589</td>
+
-
</tr>
+
-
<tr>
+
-
<td>gat</td>
+
-
<td>D</td>
+
-
<td>12.3035</td>
+
-
<td>0.478411</td>
+
-
</tr>
+
-
<tr>
+
-
<td>gga</td>
+
-
<td>G</td>
+
-
<td>15.4126</td>
+
-
<td>0.282343</td>
+
-
</tr>
+
-
<tr>
+
-
<td>ggg</td>
+
-
<td>G</td>
+
-
<td>13.325</td>
+
-
<td>0.244101</td>
+
-
</tr>
+
-
<tr>
+
-
<td>ggc</td>
+
-
<td>G</td>
+
-
<td>15.768</td>
+
-
<td>0.288853</td>
+
-
</tr>
+
-
<tr>
+
-
<td>ggt</td>
+
-
<td>G</td>
+
-
<td>10.0826</td>
+
-
<td>0.184703</td>
+
-
</tr>
+
-
<tr>
+
-
<td>gca</td>
+
-
<td>A</td>
+
-
<td>18.5218</td>
+
-
<td>0.287388</td>
+
-
</tr>
+
-
<tr>
+
-
<td>gcg</td>
+
-
<td>A</td>
+
-
<td>17.8111</td>
+
-
<td>0.276361</td>
+
-
</tr>
+
-
<tr>
+
-
<td>gcc</td>
+
-
<td>A</td>
+
-
<td>13.4139</td>
+
-
<td>0.208132</td>
+
-
</tr>
+
-
<tr>
+
-
<td>gct</td>
+
-
<td>A</td>
+
-
<td>14.702</td>
+
-
<td>0.228119</td>
+
-
</tr>
+
-
<tr>
+
-
<td>gta</td>
+
-
<td>V</td>
+
-
<td>8.17269</td>
+
-
<td>0.173749</td>
+
-
</tr>
+
-
<tr>
+
-
<td>gtg</td>
+
-
<td>V</td>
+
-
<td>15.457</td>
+
-
<td>0.328612</td>
+
-
</tr>
+
-
<tr>
+
-
<td>gtc</td>
+
-
<td>V</td>
+
-
<td>13.3695</td>
+
-
<td>0.28423</td>
+
-
</tr>
+
-
<tr>
+
-
<td>gtt</td>
+
-
<td>V</td>
+
-
<td>10.0382</td>
+
-
<td>0.213409</td>
+
-
</tr>
+
-
<tr>
+
-
<td>caa</td>
+
-
<td>Q</td>
+
-
<td>16.4786</td>
+
-
<td>0.39135</td>
+
-
</tr>
+
-
<tr>
+
-
<td>cag</td>
+
-
<td>Q</td>
+
-
<td>25.6285</td>
+
-
<td>0.60865</td>
+
-
</tr>
+
-
<tr>
+
-
<td>cac</td>
+
-
<td>H</td>
+
-
<td>13.7248</td>
+
-
<td>0.463268</td>
+
-
</tr>
+
-
<tr>
+
-
<td>cat</td>
+
-
<td>H</td>
+
-
<td>15.9012</td>
+
-
<td>0.536732</td>
+
-
</tr>
+
-
<tr>
+
-
<td>cga</td>
+
-
<td>R</td>
+
-
<td>17.7223</td>
+
-
<td>0.16871</td>
+
-
</tr>
+
-
<tr>
+
-
<td>cgg</td>
+
-
<td>R</td>
+
-
<td>15.7235</td>
+
-
<td>0.149683</td>
+
-
</tr>
+
-
<tr>
+
-
<td>cgc</td>
+
-
<td>R</td>
+
-
<td>16.0789</td>
+
-
<td>0.153066</td>
+
-
</tr>
+
-
<tr>
+
-
<td>cgt</td>
+
-
<td>R</td>
+
-
<td>12.3035</td>
+
-
<td>0.117125</td>
+
-
</tr>
+
-
<tr>
+
-
<td>cca</td>
+
-
<td>P</td>
+
-
<td>21.5422</td>
+
-
<td>0.319921</td>
+
-
</tr>
+
-
<tr>
+
-
<td>ccg</td>
+
-
<td>P</td>
+
-
<td>18.9216</td>
+
-
<td>0.281003</td>
+
-
</tr>
+
-
<tr>
+
-
<td>ccc</td>
+
-
<td>P</td>
+
-
<td>10.5712</td>
+
-
<td>0.156992</td>
+
-
</tr>
+
-
<tr>
+
-
<td>cct</td>
+
-
<td>P</td>
+
-
<td>16.301</td>
+
-
<td>0.242084</td>
+
-
</tr>
+
-
<tr>
+
-
<td>cta</td>
+
-
<td>L</td>
+
-
<td>6.88461</td>
+
-
<td>0.0693202</td>
+
-
</tr>
+
-
<tr>
+
-
<td>ctg</td>
+
-
<td>L</td>
+
-
<td>24.4292</td>
+
-
<td>0.245975</td>
+
-
</tr>
+
-
<tr>
+
-
<td>ctc</td>
+
-
<td>L</td>
+
-
<td>20.565</td>
+
-
<td>0.207066</td>
+
-
</tr>
+
-
<tr>
+
-
<td>ctt</td>
+
-
<td>L</td>
+
-
<td>17.3226</td>
+
-
<td>0.174419</td>
+
-
</tr>
+
-
<tr>
+
-
<td>taa</td>
+
-
<td>S</td>
+
-
<td>17.3226</td>
+
-
<td>0.174419</td>
+
-
</tr>
+
-
</table>
+
-
</body>
 
</html>
</html>

Latest revision as of 14:41, 6 October 2014

Cambridge iGEM 2014

  • Keep in contact with us...

Edit this page

Informatics

To introduce and facilitate future use of the novel chassis Marchantia polymorpha, we computationally analysed its genome to:
  • find out the most efficient codon usage in order to optimise our Marchantia specific registry parts and facilitate all future synthetic biology work on Marchantia;
  • submit a small library of Marchantia promoters to the iGEM registry, in particular looking for those which are strong, inducible, tissue-specific or expressed in an early development stage.

Codon usage optimisation

Our start point was the Marchantia genome and the mRNA transcriptome predicted with Geneious software (http://www.geneious.com/). The data was given to us from Jim's lab by Bernardo. 99 000 ORFs were predicted, which seems too large to be a realistic number of genes expressed by this small liverwort. Half of these were only 100 amino acids long. We set the threshold for candidate genes among these at 300 amino acids, obtaining the expected normal distribution of lengths.

Using BLAST software (http://blast.ncbi.nlm.nih.gov/Blast.cgi), we compared the proteins coded by these candidate genes to the proteins present in Arabidopsis thaliana, given on Araport (www.araport.org). Some of the sequences showed incomplete matches, indicated that our predicted ORFs should be regarded with some vigilance. The longest sequence showed a 40% match, a small number as expected.

The relevant DNA sequences from the candidate mRNAs using the BLAST output. Then the different codons were counted in these genes, in order to obtain a codon table for Marchantia. While the table is not strikingly similar to Arabidopsis', we can note a similarity in the slight preference for C over other bases at the end of codons and that for G-p-C sites.

A .xls version of this table with relevant filters can be found here.

Promoter identification

To identify potential promoter regions in Marchantia's genome, we looked for gene sequences from other species, mainly Arabidopsis thaliana and Physcomitrella patens. We were particularly interested in genes that were:

  • Nitrate inducible
  • Sulphate inducible
  • Phosphate inducible
  • Circadian rythm driven
  • Metabolism related
  • Development related
  • Light responsive
  • Development
  • Standard templates for proteins
  • Standard templates for genes

Reading recent research papers in search for these genes and looking for their sequences on GenBank and ThaleMine, we constructed a list of protein sequences to compare to the Marchantia predicted scaffolds. (For genes where only the nucleotide sequence was available, instead of the protein sequence, a C++ code was written to perform the translation). We ran a tblastn search in BLAST for the best matches between our candidate proteins and the Marchantia scaffolds. Then, for the best matches (~60% and above, with some judging by eye), we annotated the 2kb upstream of the start codon as a potential promoter region.

Example of a BLAST hit, matching an inducible nitrate transporter sequence to a Marchantia gene

We identified 30 candidate promoters this way, that we are planning to screen by inserting in a construct driving the yellow fluorescent protein Venus. For each promoter, we will make a construct with and one without amplification by GAL4 and GAL4 UAS, to evaluate the promoter strength and get around any leakages due to the use of GAL4.

Edit this page

Retrieved from "http://2014.igem.org/Team:Cambridge-JIC/Informatics"