The input for the FLUX CAPACITOR is the annotation of a reference transcriptome and reads fromRNAseq technologies aligned to the genome. From the reference annotation, splicing graphs are produced and reads are mapped to corresponding edges in these graphs according to the position where they align in the genomic sequence. The resulting graph with edges labelled by the number of reads can be interpreted as a flow network where each transcript representing a transportation path from its start to its end and consequently each edge a possibly shared segment of transportation along which a certain number of reads per nucleotide -- i.e., a flux -- is observed. Given a density function of reads along a transcript, the expected participation of each transcript in an edge under consideration can be estimated. The basic idea is to cast back from these latter participations and the observed number of reads - allowing for a certain amount of noise - to the original transcript abundancies. To do so, a linear constraint is formalized for each edge, and an optimal solution for the complete set of constraints is found by a standard linear program solver. |
|
Transcriptome genetics using second generation sequencing in a Caucasian population. Montgomery SB, Sammeth M, Gutierrez-Arcelus M, Lach RP, Ingle C, Nisbett J, Guigo R, Dermitzakis ET. Nature. 2010 Apr 1;464(7289):773-7. Epub 2010 Mar 10. PMID: 20220756 |