Development and benchmarking of improved computational methods for transcript-level expression analysis using RNA-seq data

Abstracts are not currently available in GtR for all funded research. This is normally because the abstract was not required at the time of proposal submission, but may be because it included sensitive information such as personal details.

RNA-seq technology enables the discovery and quantification of multiple transcripts for each gene, including different gene isoforms and different allelic forms. We propose the development of a Bayesian inference approach for inferring the concentration of different transcripts present in a sample by using a probabilistic model of mapped reads. By using a Bayesian inference approach we will capture the level of inherent uncertainty in our estimates of transcript expression levels due to mapping ambiguity, technical noise, read depth limitations and biological noise. We will include the possibility of discovering unannotated isoforms. The use of a read-level probabilistic model will allow us to incorporate information about read density biases and read mapping quality scores. We will apply the model to quantify allele-specific isoform expression which is particularly challenging in complex genomes such as the hexaploid wheat that can express genes from a set of three diploid genomes. We will develop a transcript-level benchmark dataset for method evaluation in which different gene isoforms are spiked in at known concentrations against a natural background. R-code implementing our methods for transcript-level inference and benchmarking will be disseminated through the Bioconductor project. We will extend the existing puma Bioconductor package for noise propagation in microarray analysis so that the methods there can be applied to transcript-level expression data with an associated multivariate uncertainty distribution.

Communication and Engagement: We will publish papers in open access peer-reviewed journals so that the academic community are made aware of developments. Software will be implemented as open source Bioconductor packages. A public benchmark will lead to better practice by allowing a publically available comparison of competing methods. We have close links to TGAC and the other MRC hubs and we will ensure that all of these groups are made aware of the tools developed and their application. The CGR, as a NERC and MRC hub, also works with a large bioinformatics community and will train new users in working with this software.

Collaboration and Co-production: The investigators are also engaged in many other BBSRC projects which can adopt the methodology developed here to add value to those projects. These projects also provide excellent application data for this proposal. Many of these projects involve short read sequencing of economically important species and comparative analysis to model species and we will identify other projects where the software will be deployed and ensure that their feedback is reflected in the development of the software.

Exploitation and Application: As this tool will be deployed primarily for academic research we do not intend to protect its application. It will be made freely available to the user community through a suitable open source license.

Capacity and Involvement: We are involved in supervising BBSRC and MRC funded Ph.D. students who will benefit from this research as they will be directly using the software developed and we regularly employed sixth form students to undertake research activities in the lab. Both SITRAN and the CGR undertake a wide range of outreach activities to industry, the academic community and the general public and actively engage with the media at local, national and international level

Impact Activity Deliverables and Milestones: Computational Biology developments will be presented at international conferences. Four key papers and associated software will be published along with a benchmarking website.