Abstract:BackgroundThe 15 sibling species of the Paramecium aurelia cryptic species complex emerged after a whole genome duplication that occurred tens of millions of years ago. Given extensive knowledge of the genetics and epigenetics of Paramecium acquired over the last century, this species complex offers a uniquely powerful system to investigate the consequences of whole genome duplication in a unicellular eukaryote as well as the genetic and epigenetic mechanisms that drive speciation. High quality Paramecium gene models are important for research using this system. The major aim of the work reported here was to build an improved gene annotation pipeline for the Paramecium lineage.ResultsWe generated oriented RNA-Seq transcriptome data across the sexual process of autogamy for the model species Paramecium tetraurelia. We determined, for the first time in a ciliate, candidate P. tetraurelia transcription start sites using an adapted Cap-Seq protocol. We developed TrUC, multi-threaded Perl software that in conjunction with TopHat mapping of RNA-Seq data to a reference genome, predicts transcription units for the annotation pipeline. We used EuGene software to combine annotation evidence. The high quality gene structural annotations obtained for P. tetraurelia were used as evidence to improve published annotations for 3 other Paramecium species. The RNA-Seq data were also used for differential gene expression analysis, providing a gene expression atlas that is more sensitive than the previously established microarray resource.ConclusionsWe have developed a gene annotation pipeline tailored for the compact genomes and tiny introns of Paramecium species. A novel component of this pipeline, TrUC, predicts transcription units using Cap-Seq and oriented RNA-Seq data. TrUC could prove useful beyond Paramecium, especially in the case of high gene density. Accurate predictions of 3′ and 5′ UTR will be particularly valuable for studies of gene expression (e.g. nucleosome positioning, identification of cis regulatory motifs). The P. tetraurelia improved transcriptome resource, gene annotations for P. tetraurelia, P. biaurelia, P. sexaurelia and P. caudatum, and Paramecium-trained EuGene configuration are available through ParameciumDB (http://paramecium.i2bc.paris-saclay.fr). TrUC software is freely distributed under a GNU GPL v3 licence (https://github.com/oarnaiz/TrUC).

Isolation and expression of two genes encoding eukaryotic release factor 1 from Paramecium tetraurelia

Different modes of stop codon restriction by the Stylonychia and Paramecium eRF1 translation termination factors

Class I release factors in ciliates with variant genetic codes

Functional Analysis of E. Coli Class II Release Factor RF3 Mutants

Crystal Structure and Functional Analysis of the Eukaryotic Class II Release Factor Erf3 from S. Pombe

Stop codon recognition in ciliates: Euplotes release factor does not respond to reassigned UGA codon

Stop codon recognition in the early-diverged protozoans Giardia lamblia and Trichomonas vaginalis

The Crystal Structure of Human Eukaryotic Release Factor Erf1—mechanism of Stop Codon Recognition and Peptidyl-tRNA Hydrolysis

Omnipotent decoding potential resides in eukaryotic translation termination factor eRF1 of variant-code organisms and is modulated by the interactions of amino acid sequences within domain 1

Genetic interaction between yeast Saccharomyces cerevisiae release factors and the decoding region of 18 S rRNA

A precedented nuclear genetic code with all three termination codons reassigned as sense codons in the syndinean Amoebophrya sp. ex Karlodinium veneficum

Structural Insights into Erf3 and Stop Codon Recognition by Erf1.

Functional Activity of Isoform 2 of Human eRF1

Structural basis for stop codon recognition in eukaryotes

Genome-wide Evolutionary Analysis of the Noncoding RNA Genes and Noncoding DNA of Paramecium Tetraurelia

Decoding accuracy in eRF1 mutants and its correlation with pleiotropic quantitative traits in yeast

The coding/non-coding overlapping architecture of the gene encoding the Drosophila pseudouridine synthase

The Sub-Chromosomic Macronuclear Pheromone Genes of the Ciliate Euplotes raikovi: Comparative Structural Analysis and Insights into the Mechanism of Expression

Analysis of micronuclear, macronuclear and cDNA sequences encoding the regulatory subunit of cAMP-dependent protein kinase of Euplotes octocarinatus: evidence for a ribosomal frameshift

Improved methods and resources for paramecium genomics: transcription units, gene annotation and gene expression