Abstract
Accurate detection of oncogenic gene fusions is becoming increasingly important given the availability of highly effective targeted therapies. However, their identification in clinical practice remains challenging due to the rarity of individual events, diversity of partner genes, and variability of breakpoint locations. Conventional approaches such as immunohistochemistry (IHC) and fluorescence in situ hybridization (FISH) lack multiplexing capacity and demonstrate variable sensitivity and specificity, while direct identification of fusion transcripts in whole-transcriptome sequencing (RNA-seq) profiles provides broader applicability but limited sensitivity, as fusion junctions are frequently supported by a minimal number of reads or even no reads at all. In this study, a novel approach was employed to accurately detect clinically actionable RET (REarranged during Transfection) fusions. This approach entailed the measurement of the imbalance in RNA-seq read coverage of potential fusion oncogenes at their 3' and 5' exons. A total of 1327 experimental solid tumor RNA-seq profiles were screened, including 154 non-small cell lung cancer and 221 thyroid cancer samples. The RET status was validated in 78 selected cases by targeted NGS and Sanger sequencing. An analysis of the coverage imbalance was conducted, which enabled the accurate discrimination between true and false positive RET fusions. This approach outperformed other methods and yielded 100% sensitivity and specificity with optimized thresholds. The findings were validated using an independent cohort of 79 thyroid cancer cases, confirming the reliability of the results. Among the 18 RET fusion-positive samples, one was identified as an extremely rare case (RUFY3::RET), and two were determined to be novel fusions (FN1::RET, PPP1R21::RET). The findings of this study demonstrate that exon coverage imbalance analysis serves as a robust complement to computational RNA-seq analysis pipelines for the detection of clinically relevant RET fusions.