Difference between revisions of "Team:Uppsala/Transcriptomics/Bioinformatics"
| (34 intermediate revisions by 6 users not shown) | |||
| Line 19: | Line 19: | ||
} | } | ||
| + | |||
| + | .side-img a{ | ||
| + | display: inline-block; | ||
| + | color: black; | ||
| + | padding-left: 5px; | ||
| + | text-decoration: none; | ||
| + | } | ||
| + | |||
| + | .inner-card-text a{ | ||
| + | display: inline-block; | ||
| + | color: black; | ||
| + | padding-left: 5px; | ||
| + | text-decoration: none; | ||
| + | } | ||
| + | |||
</style> | </style> | ||
| Line 30: | Line 45: | ||
| − | <div class="svg-wrapper"> | + | <div class="svg-wrapper" id="Project_Description"> |
| Line 167: | Line 182: | ||
<div class="body"> | <div class="body"> | ||
<div class="parallax"></div> | <div class="parallax"></div> | ||
| − | + | <div class="igem-icon"><a href="https://2018.igem.org/Team:Uppsala"><img src="https://static.igem.org/mediawiki/2018/c/cf/T--Uppsala--WormBusterLogo_Black.png"></a></div> | |
<div class= "content blur-box" style="font-size:16px;"> | <div class= "content blur-box" style="font-size:16px;"> | ||
| Line 184: | Line 199: | ||
<div id="toctitle"></div> | <div id="toctitle"></div> | ||
<ul> | <ul> | ||
| − | <li class="toclevel tocsection"><a href="#Project_Description" class="scroll"> <span id="whereYouAre"> | + | <li class="toclevel tocsection"><a href="#Project_Description" class="scroll"> <span id="whereYouAre"> Bioinformatics</span> </a> |
<ul> | <ul> | ||
| − | <li class="toclevel nav-item active"><a href="# | + | <li class="toclevel nav-item active"><a href="#Exp" class="nav-link scroll"> Experiment</a></li> |
| − | <li class="toclevel nav-item"><a href="# | + | <li class="toclevel nav-item"><a href="#Results" class="nav-link scroll"> Results</a></li> |
| − | + | ||
<li class="toclevel nav-item"><a href="#References" class="nav-link scroll"> References </a></li> | <li class="toclevel nav-item"><a href="#References" class="nav-link scroll"> References </a></li> | ||
</ul> | </ul> | ||
| Line 229: | Line 243: | ||
<p>After a succesfull sequencing has been performed and you’re left with raw data containing millions and millions (and millions) of lines of base sequences, all of this needs to be processed and interpreted. This is where the interdisciplinary field of bioinformatics comes in. A vast range of software tools are available, tailored to different kinds of analysis as well as being unique to the different sequencing methods being used.<br><br> | <p>After a succesfull sequencing has been performed and you’re left with raw data containing millions and millions (and millions) of lines of base sequences, all of this needs to be processed and interpreted. This is where the interdisciplinary field of bioinformatics comes in. A vast range of software tools are available, tailored to different kinds of analysis as well as being unique to the different sequencing methods being used.<br><br> | ||
| − | Most of the tools we used were available through the free website Usegalaxy.org which as well let us do the processing on their servers. Because we also made use of nanopore sequencing, tailored tools used for the MinION data were available from their community hub which could be run from a terminal window. </p | + | Most of the tools we used were available through the free website Usegalaxy.org which as well let us do the processing on their servers. Because we also made use of nanopore sequencing, tailored tools used for the MinION data were available from their community hub which could be run from a terminal window. </p> |
| − | <h2>Experiment</h2> | + | <h2 id="Exp">Experiment</h2> |
<p>We decided to create our bioinformatics pipeline from scratch. This was not an easy task however as nanopore technology is novel and many of the available pipelines are tailored to illumina sequencing. Generally though, a basic transcriptomics pipeline looks like the following: Alignment to a reference genome, gene counting and differential gene expression [1]. However a couple of data processing steps were needed for the nanopore data beforehand such as demultiplexing and adapter trimming.</p><br> | <p>We decided to create our bioinformatics pipeline from scratch. This was not an easy task however as nanopore technology is novel and many of the available pipelines are tailored to illumina sequencing. Generally though, a basic transcriptomics pipeline looks like the following: Alignment to a reference genome, gene counting and differential gene expression [1]. However a couple of data processing steps were needed for the nanopore data beforehand such as demultiplexing and adapter trimming.</p><br> | ||
| − | <h3>Demultiplexing and | + | <h3>Demultiplexing and Adapter Trimming</h3> |
</div> | </div> | ||
| Line 246: | Line 260: | ||
<!-- Here you put your paragraphs --> | <!-- Here you put your paragraphs --> | ||
<p>Because the sequencing itself runs pooled samples containing both the barcoded cultured- and control-group samples, the data produced needs to be demultiplexed i.e separated into files containing the reads from respective groups. Because the barcodes used to fingerprint each group is made up of its own base sequence, this also had to be removed or ”trimmed” from the data, leaving us with the pure mRNA sequences. This was achieved using a free nanopore community tool called porechop.</p><br> | <p>Because the sequencing itself runs pooled samples containing both the barcoded cultured- and control-group samples, the data produced needs to be demultiplexed i.e separated into files containing the reads from respective groups. Because the barcodes used to fingerprint each group is made up of its own base sequence, this also had to be removed or ”trimmed” from the data, leaving us with the pure mRNA sequences. This was achieved using a free nanopore community tool called porechop.</p><br> | ||
| − | <h3>Genome | + | <h3>Genome Alignment</h3> |
<p>The base sequences needs to be aligned to the reference genome of the sequenced species in question for the downstream data analysis. This is important because we want to know where each sequence actually lies in the genome and which genes they correspond to. Genome alignment was done using another community tool called minimap2.</p> | <p>The base sequences needs to be aligned to the reference genome of the sequenced species in question for the downstream data analysis. This is important because we want to know where each sequence actually lies in the genome and which genes they correspond to. Genome alignment was done using another community tool called minimap2.</p> | ||
| Line 254: | Line 268: | ||
<div class="side-img" style="background-color:darkolivegreen;"> | <div class="side-img" style="background-color:darkolivegreen;"> | ||
<!-- Here goes the big image to the right --> | <!-- Here goes the big image to the right --> | ||
| − | + | <img src="https://static.igem.org/mediawiki/2018/3/3b/T--Uppsala--Transcriptomics-Demultiplexing.png"> | |
| − | <p><b>Figure 1 | + | <a href="https://static.igem.org/mediawiki/2018/3/3b/T--Uppsala--Transcriptomics-Demultiplexing.png"><p><b>Figure 1.</b> Running demultiplexing and barcode trimming from the terminal. The programme first separates the reads according to barcode and then searches for available possible barcodes to be trimmed off.</p></a> |
</div> | </div> | ||
| Line 264: | Line 278: | ||
<!--End of template with side picture --> | <!--End of template with side picture --> | ||
| − | + | ||
<div class="card-holder"> | <div class="card-holder"> | ||
| − | <h3>Gene | + | <h3>Gene Counting</h3> |
</div> | </div> | ||
| Line 287: | Line 301: | ||
<div class="side-img" style="background-color:darkolivegreen;"> | <div class="side-img" style="background-color:darkolivegreen;"> | ||
<!-- Here goes the big image to the right --> | <!-- Here goes the big image to the right --> | ||
| − | + | ||
| − | <p><b>Figure 2 | + | <img src="https://static.igem.org/mediawiki/2018/a/a9/T--Uppsala--Transcriptomics-Bioinformatics2.png"> |
| + | <a href="https://static.igem.org/mediawiki/2018/a/a9/T--Uppsala--Transcriptomics-Bioinformatics2.png"><p><b>Figure 2.</b> Results of a differential gene expression analysis using Deseq2 on test files. The genes (shown with their gene ID) as well as their mean base length and several statistical results can be seen.</p></a> | ||
</div> | </div> | ||
| Line 298: | Line 313: | ||
| − | <h2>Result</h2> | + | <h2 id="Results">Result</h2> |
| − | <p>The transcriptomics pipeline was tried out and validated using read files available from the internet. The files consisted of two datasets of <i>E. | + | <h3>Validating our Transcriptomics Pipeline</h3> |
| + | <p>The transcriptomics pipeline was tried out and validated using read files available from the internet. The files consisted of two datasets of <i>E. coli</i> (triplicates) cultured in regular LB and a sugar solution respectively.</p><br> | ||
| Line 320: | Line 336: | ||
<div class="inner-card-text"> | <div class="inner-card-text"> | ||
<!-- start of paragraph--> | <!-- start of paragraph--> | ||
| − | <p><b>Figure 3 | + | <a href="https://static.igem.org/mediawiki/2018/a/a9/T--Uppsala--Transcriptomics-Bioinformatics2.png"><p><b>Figure 3.</b> Results of the differential gene expression analysis using Deseq2 on test files. The genes (shown with their gene ID) as well as their mean base length and several statistical results can be seen.</p></a> |
</div> | </div> | ||
<!-- end of paragraph --> | <!-- end of paragraph --> | ||
| Line 327: | Line 343: | ||
<br> | <br> | ||
| − | |||
| − | |||
<img class="content-card-img" src="https://static.igem.org/mediawiki/2018/8/81/T--Uppsala--Transcriptomics-Bioinformatics3.png"> | <img class="content-card-img" src="https://static.igem.org/mediawiki/2018/8/81/T--Uppsala--Transcriptomics-Bioinformatics3.png"> | ||
<div class="inner-card-text"> | <div class="inner-card-text"> | ||
<!-- start of paragraph --> | <!-- start of paragraph --> | ||
| − | + | <a href="https://static.igem.org/mediawiki/2018/8/81/T--Uppsala--Transcriptomics-Bioinformatics3.png"> <p><b>Figure 4.</b> Results of the differential gene expression after filtering for statistical significance and fold change.</p></a> | |
<!-- End of paragraphs --> | <!-- End of paragraphs --> | ||
</div> | </div> | ||
| Line 350: | Line 364: | ||
<!--End of template with side picture --> | <!--End of template with side picture --> | ||
| − | |||
| − | <p>The results after searching for the genes in the NCBI database showed that the most expressed gene from the sugar-cultured <i>E. | + | |
| + | <p>The results after searching for the genes in the NCBI database showed that the most expressed gene from the sugar-cultured <i>E. coli</i> was shown to be involved in a type of sugar system, proving that the pipeline was indeed working.</p><br> | ||
</div> | </div> | ||
| Line 370: | Line 384: | ||
<!--change src to that of the image you want--> | <!--change src to that of the image you want--> | ||
| − | + | <img class="content-card-img" src="https://static.igem.org/mediawiki/2018/4/4c/T--Uppsala--Transcriptomics-Bioinformatics4.png"> | |
<div class="inner-card-text"> | <div class="inner-card-text"> | ||
<!-- start of paragraph--> | <!-- start of paragraph--> | ||
| − | <p><b>Figure 5 | + | <a href="https://static.igem.org/mediawiki/2018/4/4c/T--Uppsala--Transcriptomics-Bioinformatics4.png"><p><b>Figure 5.</b> Highly expressed gene produced from the pipeline matching a glucose specific gene.</p></a> |
</div> | </div> | ||
<!-- end of paragraph --> | <!-- end of paragraph --> | ||
| Line 382: | Line 396: | ||
| − | + | <img class="content-card-img" src="https://static.igem.org/mediawiki/2018/3/3c/T--Uppsala--Transcriptomics-Bioinformatics5.png"> | |
<div class="inner-card-text"> | <div class="inner-card-text"> | ||
<!-- start of paragraph --> | <!-- start of paragraph --> | ||
| − | + | <a href="https://static.igem.org/mediawiki/2018/3/3c/T--Uppsala--Transcriptomics-Bioinformatics5.png"> <p><b>Figure 6.</b> Results of the differential gene expression done on our own data.</p></a> | |
<!-- End of paragraphs --> | <!-- End of paragraphs --> | ||
</div> | </div> | ||
| Line 398: | Line 412: | ||
<div class="card-holder"> | <div class="card-holder"> | ||
| − | + | <br> | |
| − | <p><b>Table 1</b> | + | <h3>Analyzing Our Own Sequencing Data</h3> |
| − | promotor sequence and function in the organism.</p> | + | <p><b>Table 1.</b> The first few genes as a result of the differential gene expression analysis seen in Figure 6 |
| − | <!--Start of template with side picutre --> | + | together with their promotor sequence and function in the organism.</p> |
| + | <!--Start of template with side picutre --> | ||
<!-- Here you put your paragraphs --> | <!-- Here you put your paragraphs --> | ||
| − | <table class= | + | <table class="pgrouptable tablesorter our-table" style="width: 100%;" cellspacing="0" cellpadding="0"> |
| − | + | <thead><tr> | |
| − | <th style= “width: auto”>Gene ID</th> | + | <th style= “width: auto”>Gene ID</th> |
| − | <th style=“width: auto” >Gene name</th> | + | <th style= “width: auto” >Gene name</th> |
| − | <th style=“width: auto” >Promotor sequence</th> | + | <th style= “width: auto” >Promotor sequence</th> |
| − | <th style=“width: auto” >Function</th> | + | <th style= “width: auto” >Function</th> |
| − | </tr></thead> | + | <th style= “width: auto” >Fold change</th> |
| − | <tbody><tr> | + | </tr></thead> |
| + | <tbody><tr> | ||
| − | <td> | + | <td> |
| − | ER3413_45<br> | + | ER3413_45<br> |
| − | ER3413_70<br> | + | ER3413_70<br> |
| − | ER3413_87<br> | + | ER3413_87<br> |
| − | ER3413_126<br> | + | ER3413_126<br> |
| − | ER3413_173 | + | ER3413_173 |
| − | </td> | + | </td> |
| − | <td > | + | <td > |
| − | apaG<br> | + | apaG<br> |
| − | leuA<br> | + | leuA<br> |
| − | murG<br> | + | murG<br> |
| − | panD<br> | + | panD<br> |
| − | frr | + | frr |
| − | </td> | + | </td> |
| − | <td> | + | <td> |
| − | ggcaccatgcagggtcactacgaaatgatcgatgaaa<br> | + | ggcaccatgcagggtcactacgaaatgatcgatgaaa<br> |
| − | ttgacatccgtttttgtatccagtaactctaaaagc<br> | + | ttgacatccgtttttgtatccagtaactctaaaagc<br> |
| − | <br> | + | <p>-</p><br> |
| − | tagacactaaacaaaaatcgggcaatactgcgtga<br> | + | tagacactaaacaaaaatcgggcaatactgcgtga<br> |
| − | ttacccgtaatatgtttaatcagggctatacttagcac | + | ttacccgtaatatgtttaatcagggctatacttagcac |
| − | </td> | + | </td> |
| − | <td> | + | <td> |
| − | protein associated with Co2+ and Mg2+ efflux<br> | + | protein associated with Co2+ and Mg2+ efflux<br> |
| − | 2-isopropylmalate synthase<br> | + | 2-isopropylmalate synthase<br> |
| − | N-acetylglucosaminyl transferase<br> | + | N-acetylglucosaminyl transferase<br> |
| − | putative inner membrane protein<br> | + | putative inner membrane protein<br> |
| − | inner membrane protein, UPF0118 family | + | inner membrane protein, UPF0118 family |
| − | </td> | + | </td> |
| − | </tr><tr> | + | <td> |
| − | </tr></tbody></table> | + | 0.40<br> |
| + | 0.40<br> | ||
| + | 0.40<br> | ||
| + | 0.40<br> | ||
| + | 0.40 | ||
| + | </td> | ||
| + | </tr><tr> | ||
| + | </tr></tbody></table> | ||
<!-- End of Code For TABLE --> | <!-- End of Code For TABLE --> | ||
| − | |||
| − | + | ||
| + | <p>The resuts from our runs unfortunately did not produce as good results as seen above. Due to the major issues with sequencing and actually generating enough data, it can be seen in figure 6 what kind of effects it had. Judging by the adjusted p-values it is clear that even though the genes can indeed be identified as seen in Table 1 the statistical significance is extremely uncertain (the minimal accepted threshold is an adjusted p-value of <= 0.05). Any up-or down regulation of fold-change of interest was not able to be identified either. Looking at these errors it can be assumed that no major change in fold-change as well as low significancy is due to simply not enough data being generated from the prior sequencing step. Because of these facts no gene could be identified as a possible candidate for our reporter system.</p> | ||
| − | + | ||
| − | + | ||
</div> | </div> | ||
| + | <div class="card-holder"> | ||
| + | <h2 id="References">References</h2> | ||
| + | |||
| + | <p><b>[1]</b> Galaxyproject, 2018. Reference-based RNA-Seq data analysis <a href="https://galaxyproject.github.io/training-material/topics/transcriptomics/tutorials/ref-based/tutorial.html">Galaxyproject</a> Date of visit 2018-10-15</p> | ||
| + | </div> | ||
<!-- HERE ENDS THE PORTION WHERE YOU PUT IN YOUR CONTENT--> | <!-- HERE ENDS THE PORTION WHERE YOU PUT IN YOUR CONTENT--> | ||
<div style="height:5em;"></div> | <div style="height:5em;"></div> | ||
</div> | </div> | ||
</div> | </div> | ||
| − | + | </div> | |
| − | + | </body> | |
</html> | </html> | ||
Latest revision as of 23:37, 17 October 2018
Bioinformatics
After a succesfull sequencing has been performed and you’re left with raw data containing millions and millions (and millions) of lines of base sequences, all of this needs to be processed and interpreted. This is where the interdisciplinary field of bioinformatics comes in. A vast range of software tools are available, tailored to different kinds of analysis as well as being unique to the different sequencing methods being used.
Most of the tools we used were available through the free website Usegalaxy.org which as well let us do the processing on their servers. Because we also made use of nanopore sequencing, tailored tools used for the MinION data were available from their community hub which could be run from a terminal window.
Experiment
We decided to create our bioinformatics pipeline from scratch. This was not an easy task however as nanopore technology is novel and many of the available pipelines are tailored to illumina sequencing. Generally though, a basic transcriptomics pipeline looks like the following: Alignment to a reference genome, gene counting and differential gene expression [1]. However a couple of data processing steps were needed for the nanopore data beforehand such as demultiplexing and adapter trimming.
Demultiplexing and Adapter Trimming
Because the sequencing itself runs pooled samples containing both the barcoded cultured- and control-group samples, the data produced needs to be demultiplexed i.e separated into files containing the reads from respective groups. Because the barcodes used to fingerprint each group is made up of its own base sequence, this also had to be removed or ”trimmed” from the data, leaving us with the pure mRNA sequences. This was achieved using a free nanopore community tool called porechop.
Genome Alignment
The base sequences needs to be aligned to the reference genome of the sequenced species in question for the downstream data analysis. This is important because we want to know where each sequence actually lies in the genome and which genes they correspond to. Genome alignment was done using another community tool called minimap2.
Gene Counting
Gene counting basically means that you count how many times each mRNA sequence (aligned over a gene from the previous step) occurs. This in turn directly correlates to the amount of up- or down-regulation of that particular gene. A lot of different tools were available for gene counting but ”featureCounts” was chosen through galaxy.
After the differential gene expression analysis is done the data was filtered twice, one time for the best adjusted P-value and subsequently for the highest (meaning the most significant) fold changes. Left were a couple of candidate genes which could be easily identified by their gene ID through various databases such as NCBI.
Result
Validating our Transcriptomics Pipeline
The transcriptomics pipeline was tried out and validated using read files available from the internet. The files consisted of two datasets of E. coli (triplicates) cultured in regular LB and a sugar solution respectively.
The results after searching for the genes in the NCBI database showed that the most expressed gene from the sugar-cultured E. coli was shown to be involved in a type of sugar system, proving that the pipeline was indeed working.
Analyzing Our Own Sequencing Data
Table 1. The first few genes as a result of the differential gene expression analysis seen in Figure 6 together with their promotor sequence and function in the organism.
| Gene ID | Gene name | Promotor sequence | Function | Fold change |
|---|---|---|---|---|
|
ER3413_45 ER3413_70 ER3413_87 ER3413_126 ER3413_173 |
apaG leuA murG panD frr |
ggcaccatgcagggtcactacgaaatgatcgatgaaa ttgacatccgtttttgtatccagtaactctaaaagc - tagacactaaacaaaaatcgggcaatactgcgtga ttacccgtaatatgtttaatcagggctatacttagcac |
protein associated with Co2+ and Mg2+ efflux 2-isopropylmalate synthase N-acetylglucosaminyl transferase putative inner membrane protein inner membrane protein, UPF0118 family |
0.40 0.40 0.40 0.40 0.40 |
The resuts from our runs unfortunately did not produce as good results as seen above. Due to the major issues with sequencing and actually generating enough data, it can be seen in figure 6 what kind of effects it had. Judging by the adjusted p-values it is clear that even though the genes can indeed be identified as seen in Table 1 the statistical significance is extremely uncertain (the minimal accepted threshold is an adjusted p-value of <= 0.05). Any up-or down regulation of fold-change of interest was not able to be identified either. Looking at these errors it can be assumed that no major change in fold-change as well as low significancy is due to simply not enough data being generated from the prior sequencing step. Because of these facts no gene could be identified as a possible candidate for our reporter system.
References
[1] Galaxyproject, 2018. Reference-based RNA-Seq data analysis Galaxyproject Date of visit 2018-10-15