Detecting haplotype-specific transcript variation in long reads with FLAIR2

Abstract

BackgroundRNA-Seq has brought forth significant discoveries regarding aberrations in RNA processing, implicating these RNA variants in a variety of diseases. Aberrant splicing and single nucleotide variants in RNA have been demonstrated to alter transcript stability, localization, and function. In particular, the upregulation of ADAR, an enzyme which mediates adenosine-to-inosine editing, has been previously linked to an increase in the invasiveness of lung ADC cells and associated with splicing regulation. Despite the functional importance of studying splicing and SNVs, short read RNA-Seq has limited the communitys ability to interrogate both forms of RNA variation simultaneously. ResultsWe employed long-read technology to obtain full-length transcript sequences, elucidating cis-effects of variants on splicing changes at a single molecule level. We have developed a computational workflow that augments FLAIR, a tool that calls isoform models expressed in long-read data, to integrate RNA variant calls with the associated isoforms that bear them. We generated nanopore data with high sequence accuracy of H1975 lung adenocarcinoma cells with and without knockdown of ADAR. We applied our workflow to identify key inosine-isoform associations to help clarify the prominence of ADAR in tumorigenesis. ConclusionsUltimately, we find that a long-read approach provides valuable insight toward characterizing the relationship between RNA variants and splicing patterns. HighlightsO_LIFLAIR2 has improved transcript isoform detection and incorporates sequence variants for haplotype-specific transcript detection. C_LIO_LIIn addition to haplotype-specific variant detection, it identifies transcript-specific RNA editing C_LIO_LIAble to identify haplotype-specific transcript isoform bias in expression C_LIO_LILong-read sequencing identifies hyperedited transcripts that are missed from short-read sequencing methods for a more comprehensive identification of ADAR targets C_LI