"
Team:Cambridge-JIC/Marchantia/Codon
From 2014.igem.org
Method
Our dataset was made of large gap read mapping transcripts obtained by mRNA sequencing conducted by the Haseloff Lab on the m. polymorpha Cam strain. Open Reading Frame (ORF) and Coding Sequence (CDS) predictions were made using the CLC bio Transcript Discovery plugin[1].
We used ran blastx on our data [2] to verify the validity of the predicted ORFs and to compare the sequences with the proteins present in Arabidopsis3. A list of candidate genes was compiled from the results and this list was used to calculate a frequency table of codon usage for m.polymorpha.
Results and Discussion
The number of Open Reading Frames (ORFs) resulting from our initial predictions totalled 99 000, which seemed too large to be realistic, and half of these were only 100 amino acids (100 aa) long. By filtering the dataset using a threshold of 300 aa for candidate genes we obtained a normal distribution of lengths, which seemed reasonable. [INSERT FIGURE plot of the distribution of lengths]
In the blastx output, the longest complete sequence match was 40%. This is small as expected given the phylogenetic distance between m.polymorpha and a. thaliana [3][4].
The Codon Usage table we obtained for m. polymorpha is not strikingly similar to that of a. thaliana, as expected from the 400 million years of evolutionary divergence between them [5].
However, there is a similarity in the slight preference for C over other bases at the end of codons and that for G-p-C sites, as can be seen in the codon table.
References
1. Qiagen, CLC bio Transcript Discovery ®, http://www.clcbio.com/clc-plugin/transcript-discovery/#description back to top
