TXGP ens63 reference

From Marcotte Lab
Jump to navigationJump to search

Overview

One of the most interesting questions we can ask with X. laevis genome would be how many genes it has. To construct gene models, we are mainly focusing on de novo transcriptome assembly approach with our RNA-seq data. However, de novo transcriptome assembly programs generate many 'false positive' transcripts. Also, because of allotetraploidy in X. laevis, transcriptome data may contain many transcript variants for each gene. So, to estimate the gene model from transcriptome data precisely, we would like to combine all transcripts candidates for each gene together, and analyze them separately. Sequence-based clustering is a natural way to do this, but we need to optimize parameters, such as %identity to define a cluster. To get some ideas for this, we have looked at genes and transcripts of several well-studied organisms.

Genes & Transcripts

This figure shows total number of genes and transcripts in each organisms. The number on top of green bar means total number of transcripts, and the number on top of blue bar means total number of genes (based on EnsEMBL v.63 annotation). The number on top of cyan bar means the number of genes that contain only one transcript.

ens63_gene_tx.small.png

Clustering of transcripts

We clustered all transcripts for each organism with usearch program with different %id cutoff. The number on top of red bar means the ratio of 'the number of clusters' to 'the number of genes'. The number on top of pink bar means the ratio of 'the number of clusters having more than one gene' to 'the number of clusters'. Although they may be very closed paralogous genes, we considered this number as 'clustering error'.

  • Ultimate goal is to match the number of genes to the number of clusters in all organisms.
  • Human and mouse have too much transcripts compared to other organisms, so we would allow 1.5x more clusters than total number of genes in these organisms.
  • Although X. tropicalis is the closest model organism, it does not have many transcripts yet. So we use D. rerio(zebrafish) to estimate 'optimal number of clusters'.
  • We would like to control 'clustering error' less than 0.10 (it may be little bit higher than conventional cutoff, i.e. 0.05. But, as mentioned earlier, it may also contain many paralogous genes, so it is unlikely that all of them are clustering errors.)

With these criteria, we determined '%id>0.80' as an optimal cutoff for transcripts clustering for gene model estimation. It does not mean that we discarded all other sequences in clusters except one representative sequence. We just grouped them for further analysis.

ens63.gene_vs_clusters.uc090.small.png ens63.gene_vs_clusters.uc085.small.png ens63.gene_vs_clusters.uc080.small.png ens63.gene_vs_clusters.uc075.small.png ens63.gene_vs_clusters.uc070.small.png ens63.gene_vs_clusters.uc060.small.png

Length of transcripts

We also looked at the distribution of transcript length with same data. Except few genes in C. elegans, and D. melanogaster, there is no gene that can make transcripts longer than 10kbp. So we could consider a transcript longer than this would be a false positive. Also, more than 90% of transcripts are longer than 500 bp except human('bottom10pct' below;maybe because a lot of short transcripts are annotated on human genome), so we may discard assembled transcripts shorter than 500bp if they are much more than 10% of total transcripts.

ens63.len_dist.small.png

HUMAN (cDNA): mean=1600.3, median=936, bottom10pct=442
HUMAN (seed): mean=1713.9, median=890, bottom10pct=432
MOUSE (cDNA): mean=1881.5, median=1234, bottom10pct=427
MOUSE (seed): mean=2112.6, median=1556, bottom10pct=454
XENTR (cDNA): mean=1893.8, median=1547, bottom10pct=645
XENTR (seed): mean=1987.2, median=1638, bottom10pct=699
DANRE (cDNA): mean=1805.1, median=1363, bottom10pct=540
DANRE (seed): mean=2072.5, median=1609, bottom10pct=594
CHICK (cDNA): mean=1719.9, median=1338, bottom10pct=442
CHICK (seed): mean=1795.2, median=1432, bottom10pct=516
DROME (cDNA): mean=2359.7, median=1870, bottom10pct=575
DROME (seed): mean=2107.6, median=1615, bottom10pct=557
CAEEL (cDNA): mean=1415.2, median=1124, bottom10pct=406
CAEEL (seed): mean=1399.6, median=1092, bottom10pct=447