A fast detection of fusion genes from paired-end RNA-seq data

LMNA::NTRK1 Fusion-positive Leiomyosarcoma: Discrepancy between DNA-based Comprehensive Genomic Profiling and RNA Sequencing

A 26-year-old man presented with a tumor in the left soleus muscle. The tumor was diagnosed as a locally advanced leiomyosarcoma. The patient was treated with irradiation followed by wide resection. One year after surgery, the patient presented with multiple lung metastases. Despite aggressive sequential chemotherapy, systemic metastatic tumors continued to develop. To explore therapeutic options for the patient, we performed DNA-based CGP with FoundationOne® CDx (F1). F1 identified an out-of-strand rearrangement of the NOS1AP::NTRK1 gene, which has not been previously reported. In contrast, RNA sequencing revealed an in-frame LMNA::NTRK1 gene, which is an oncogenic fusion gene.

RAF1 gene fusions are recurrent driver events in infantile fibrosarcoma-like mesenchymal tumors

Infantile fibrosarcomas (IFS) and congenital mesoblastic nephroma (CMN) are rare myofibroblastic tumors of infancy and early childhood commonly harboring the ETV6::NTRK3 gene fusion. IFS/CMN are considered as tumors with an 'intermediate prognosis' as they are locally aggressive, but rarely metastasize, and generally have a favorable outcome. A fraction of IFS/CMN-related neoplasms are negative for the ETV6::NTRK3 gene rearrangement and are characterized by other chimeric proteins promoting MAPK signaling upregulation. In a large proportion of these tumors, which are classified as IFS-like mesenchymal neoplasms, the contributing molecular events remain to be identified. Here, we report three distinct rearrangements involving RAF1 among eight ETV6::NTRK3 gene fusion-negative tumors with an original histological diagnosis of IFS/CMN. The three fusion proteins retain the entire catalytic domain of the kinase. Two chimeric products, GOLGA4::RAF1 and LRRFIP2::RAF1, had previously been reported as driver events in different cancers, whereas the third, CLIP1::RAF1, represents a novel fusion protein. We demonstrate that CLIP1::RAF1 acts as a bona fide oncoprotein promoting cell proliferation and migration through constitutive upregulation of MAPK signaling. We show that the CLIP1::RAF1 hyperactive behavior does not require RAS activation and is mediated by constitutive 14-3-3 protein-independent dimerization of the chimeric protein. As previously reported for the ETV6::NTRK3 fusion protein, CLIP1::RAF1 similarly upregulates PI3K-AKT signaling. Our findings document that RAF1 gene rearrangements represent a recurrent event in ETV6::NTRK3-negative IFS/CMN and provide a rationale for the use of inhibitors directed to suppress MAPK and PI3K-AKT signaling in these cancers. © 2024 The Pathological Society of Great Britain and Ireland.

RTCpredictor: identification of read-through chimeric RNAs from RNA sequencing data

Read-through chimeric RNAs are being recognized as a means to expand the functional transcriptome and contribute to cancer tumorigenesis when mis-regulated. However, current software tools often fail to predict them. We have developed RTCpredictor, utilizing a fast ripgrep tool to search for all possible exon-exon combinations of parental gene pairs. We also added exonic variants allowing searches containing common SNPs. To our knowledge, it is the first read-through chimeric RNA specific prediction method that also provides breakpoint coordinates. Compared with 10 other popular tools, RTCpredictor achieved high sensitivity on a simulated and three real datasets. In addition, RTCpredictor has less memory requirements and faster execution time, making it ideal for applying on large datasets.

Chimeric RNAs reveal putative neoantigen peptides for developing tumor vaccines for breast cancer

Introduction

We present here a strategy to identify immunogenic neoantigen candidates from unique amino acid sequences at the junctions of fusion proteins which can serve as targets in the development of tumor vaccines for the treatment of breastcancer.

Method

We mined the sequence reads of breast tumor tissue that are usually discarded as discordant paired-end reads and discovered cancer specific fusion transcripts using tissue from cancer free controls as reference. Binding affinity predictions of novel peptide sequences crossing the fusion junction were analyzed by the MHC Class I binding predictor, MHCnuggets. CD8+ T cell responses against the 15 peptides were assessed through in vitro Enzyme Linked Immunospot (ELISpot).

Results

We uncovered 20 novel fusion transcripts from 75 breast tumors of 3 subtypes: TNBC, HER2+, and HR+. Of these, the NSFP1-LRRC37A2 fusion transcript was selected for further study. The 3833 bp chimeric RNA predicted by the consensus fusion junction sequence is consistent with a read-through transcription of the 5'-gene NSFP1-Pseudo gene NSFP1 (NSFtruncation at exon 12/13) followed by trans-splicing to connect withLRRC37A2 located immediately 3' through exon 1/2. A total of 15 different 8-mer neoantigen peptides discovered from the NSFP1 and LRRC37A2 truncations were predicted to bind to a total of 35 unique MHC class I alleles with a binding affinity of IC50<500nM.); 1 of which elicited a robust immune response.

Conclusion

Our data provides a framework to identify immunogenic neoantigen candidates from fusion transcripts and suggests a potential vaccine strategy to target the immunogenic neopeptides in patients with tumors carrying the NSFP1-LRRC37A2 fusion.

On de novo Bridging Paired-end RNA-seq Data

The high-throughput short-reads RNA-seq protocols often produce paired-end reads, with the middle portion of the fragments being unsequenced. We explore if the full-length fragments can be computationally reconstructed from the sequenced two ends in the absence of the reference genome-a problem here we refer to as de novo bridging. Solving this problem provides longer, more informative RNA-seq reads, and benefits downstream RNA-seq analysis such as transcript assembly, expression quantification, and splicing differential analysis. However, de novo bridging is a challenging and complicated task owing to alternative splicing, transcript noises, and sequencing errors. It remains unclear if the data provides sufficient information for accurate bridging, let alone efficient algorithms that determine the true bridges. Methods have been proposed to bridge paired-end reads in the presence of reference genome (called reference-based bridging), but the algorithms are far away from scaling for de novo bridging as the underlying compacted de Bruijn graph (cdBG) used in the latter task often contains millions of vertices and edges. We designed a new truncated Dijkstra's algorithm for this problem, and proposed a novel algorithm that reuses the shortest path tree to avoid running the truncated Dijkstra's algorithm from scratch for all vertices for further speeding up. These innovative techniques result in scalable algorithms that can bridge all paired-end reads in a cdBG with millions of vertices. Our experiments showed that paired-end RNA-seq reads can be accurately bridged to a large extent. The resulting tool is freely available at https://github.com/Shao-Group/rnabridge-denovo.

RTCpredictor: Identification of Read-Through Chimeric RNAs from RNA Sequencing Data

Read-through chimeric RNAs are gaining attention in cancer and other research fields, yet current tools often fail in predicting them. We have thus developed the first read-through chimeric RNA specific prediction method, RTCpredictor, utilizing a fast ripgrep algorithm to search for all possible exon-exon combinations of parental gene pairs. Compared with other ten popular tools, RTCpredictor achieved top performance on both simulated and real datasets. We randomly selected up to 30 candidate read-through chimeras predicted from each software method and experimentally validated a total of 109 read-throughs and on this set, RTCpredictor outperformed all the other methods. In addition, RTCpredictor ( https://github.com/sandybioteck/RTCpredictor ) has less memory requirements and faster execution time.

Exon skipping in genes encoding lineage-defining myogenic transcription factors in rhabdomyosarcoma

Rhabdomyosarcoma (RMS) is a childhood sarcoma composed of myoblast-like cells, which suggests a defect in terminal skeletal muscle differentiation. To explore potential defects in the differentiation program, we searched for mRNA splicing variants in genes encoding transcription factors driving skeletal muscle lineage commitment and differentiation. We studied two RMS cases and identified altered splicing resulting in "skipping" the second of three exons in MYOD1. RNA-Seq data from 42 tumors and additional RMS cell lines revealed exon 2 skipping in both MYOD1 and MYF5 but not in MYF6 or MYOG. Complementary molecular analysis of MYOD1 mRNA found evidence for exon skipping in 5 additional RMS cases. Functional studies showed that so-called MYODΔEx2 protein failed to robustly induce muscle-specific genes, and its ectopic expression conferred a selective advantage in cultured fibroblasts and an RMS xenograft. In summary, we present previously unrecognized exon skipping within MYOD1 and MYF5 in RMS, and we propose that alternative splicing can represent a mechanism to alter the function of these two transcription factors in RMS.

Feasibility of whole genome and transcriptome profiling in pediatric and young adult cancers

The utility of cancer whole genome and transcriptome sequencing (cWGTS) in oncology is increasingly recognized. However, implementation of cWGTS is challenged by the need to deliver results within clinically relevant timeframes, concerns about assay sensitivity, reporting and prioritization of findings. In a prospective research study we develop a workflow that reports comprehensive cWGTS results in 9 days. Comparison of cWGTS to diagnostic panel assays demonstrates the potential of cWGTS to capture all clinically reported mutations with comparable sensitivity in a single workflow. Benchmarking identifies a minimum of 80× as optimal depth for clinical WGS sequencing. Integration of germline, somatic DNA and RNA-seq data enable data-driven variant prioritization and reporting, with oncogenic findings reported in 54% more patients than standard of care. These results establish key technical considerations for the implementation of cWGTS as an integrated test in clinical oncology.

Cancer whole-genome and transcriptome sequencing (cWGTS) has been challenging to implement in clinical settings. Here, the authors develop a workflow to deliver robust cWGTS analyses and reports within clinically-relevant timeframes for paediatric, adolescent and young adult solid tumour patients.

Utility of Multimodality Molecular Profiling for Pediatric Patients with Central Nervous System Tumors

Background

As our molecular understanding of pediatric central nervous system (CNS) tumors evolves, so too do diagnostic criteria, prognostic biomarkers, and clinical management decision-making algorithms. Here, we explore the clinical utility of wide-breadth assays including whole-exome sequencing (WES), RNA sequencing (RNAseq), and methylation array profiling as an addition to more conventional diagnostic tools for pediatric CNS tumors.

Methods

This study comprises an observational, prospective cohort followed at a single academic medical center over three years. Paired tumor and normal control specimens from 53 enrolled pediatric patients with CNS tumors underwent WES. A subset of cases also underwent RNAseq (n=28) and/or methylation array analysis (n=27).

Results

RNAseq identified driver and/or targetable fusions in 7/28 cases, including potentially targetable NTRK fusions, and uncovered possible rationalized treatment options based on outlier gene expression in 23/28 cases. Methylation profiling added diagnostic confidence (8/27 cases) or diagnostic subclassification endorsed by the WHO (10/27 cases). WES detected clinically pertinent Tier 1 or Tier 2 variants in 36/53 patients. Of these, 16/17 SNVs/indels and 10/19 copy number alterations would have been detected by current in-house conventional tests including targeted sequencing panels.

Conclusions

Over a heterogeneous set of pediatric tumors, RNAseq and methylation profiling frequently yielded clinically relevant information orthogonal to conventional methods while WES demonstrated clinically-relevant added-value primarily via copy number assessment. Longitudinal cohorts comparing targeted molecular pathology workup versus broader genomic approaches including therapeutic selection based on RNA-expression data will be necessary to further evaluate the clinical benefits of these modalities in practice.

Fusion Gene Detection Using Whole-Exome Sequencing Data in Cancer Patients

Several fusion genes are directly involved in the initiation and progression of cancers. Numerous bioinformatics tools have been developed to detect fusion events, but they are mainly based on RNA-seq data. The whole-exome sequencing (WES) represents a powerful technology that is widely used for disease-related DNA variant detection. In this study, we build a novel analysis pipeline called Fuseq-WES to detect fusion genes at DNA level based on the WES data. The same method applies also for targeted panel sequencing data. We assess the method to real datasets of acute myeloid leukemia (AML) and prostate cancer patients. The result shows that two of the main AML fusion genes discovered in RNA-seq data, PML-RARA and CBFB-MYH11, are detected in the WES data in 36 and 63% of the available samples, respectively. For the targeted deep-sequencing of prostate cancer patients, detection of the TMPRSS2-ERG fusion, which is the most frequent chimeric alteration in prostate cancer, is 91% concordant with a manually curated procedure based on four other methods. In summary, the overall results indicate that it is challenging to detect fusion genes in WES data with a standard coverage of ∼ 15–30x, where fusion candidates discovered in the RNA-seq data are often not detected in the WES data and vice versa. A subsampling study of the prostate data suggests that a coverage of at least 75x is necessary to achieve high accuracy.

ChiTaH: a fast and accurate tool for identifying known human chimeric sequences from high-throughput sequencing data

Fusion genes or chimeras typically comprise sequences from two different genes. The chimeric RNAs of such joined sequences often serve as cancer drivers. Identifying such driver fusions in a given cancer or complex disease is important for diagnosis and treatment. The advent of next-generation sequencing technologies, such as DNA-Seq or RNA-Seq, together with the development of suitable computational tools, has made the global identification of chimeras in tumors possible. However, the testing of over 20 computational methods showed these to be limited in terms of chimera prediction sensitivity, specificity, and accurate quantification of junction reads. These shortcomings motivated us to develop the first ‘reference-based’ approach termed ChiTaH (Chimeric Transcripts from High–throughput sequencing data). ChiTaH uses 43,466 non–redundant known human chimeras as a reference database to map sequencing reads and to accurately identify chimeric reads. We benchmarked ChiTaH and four other methods to identify human chimeras, leveraging both simulated and real sequencing datasets. ChiTaH was found to be the most accurate and fastest method for identifying known human chimeras from simulated and sequencing datasets. Moreover, especially ChiTaH uncovered heterogeneity of the BCR-ABL1 chimera in both bulk and single-cells of the K-562 cell line, which was confirmed experimentally.

SeekFusion - A Clinically Validated Fusion Transcript Detection Pipeline for PCR-Based Next-Generation Sequencing of RNA

Detecting gene fusions involving driver oncogenes is pivotal in clinical diagnosis and treatment of cancer patients. Recent developments in next-generation sequencing (NGS) technologies have enabled improved assays for bioinformatics-based gene fusions detection. In clinical applications, where a small number of fusions are clinically actionable, targeted polymerase chain reaction (PCR)-based NGS chemistries, such as the QIAseq RNAscan assay, aim to improve accuracy compared to standard RNA sequencing. Existing informatics methods for gene fusion detection in NGS-based RNA sequencing assays traditionally use a transcriptome-based spliced alignment approach or a de-novo assembly approach. Transcriptome-based spliced alignment methods face challenges with short read mapping yielding low quality alignments. De-novo assembly-based methods yield longer contigs from short reads that can be more sensitive for genomic rearrangements, but face performance and scalability challenges. Consequently, there exists a need for a method to efficiently and accurately detect fusions in targeted PCR-based NGS chemistries. We describe SeekFusion, a highly accurate and computationally efficient pipeline enabling identification of gene fusions from PCR-based NGS chemistries. Utilizing biological samples processed with the QIAseq RNAscan assay and in-silico simulated data we demonstrate that SeekFusion gene fusion detection accuracy outperforms popular existing methods such as STAR-Fusion, TOPHAT-Fusion and JAFFA-hybrid. We also present results from 4,484 patient samples tested for neurological tumors and sarcoma, encompassing details on some novel fusions identified.

Comparative study of bioinformatic tools for the identification of chimeric RNAs from RNA Sequencing

Chimeric RNAs are gaining more and more attention as they have broad implications in both cancer and normal physiology. To date, over 40 chimeric RNA prediction methods have been developed to facilitate their identification from RNA sequencing data. However, a limited number of studies have been conducted to compare the performance of these tools; additionally, previous studies have become outdated as more software tools have been developed within the last three years. In this study, we benchmarked 16 chimeric RNA prediction software, including seven top performers in previous benchmarking studies, and nine that were recently developed. We used two simulated and two real RNA-Seq datasets, compared the 16 tools for their sensitivity, positive prediction value (PPV), F-measure, and also documented the computational requirements (time and memory). We noticed that none of the tools are inclusive, and their performance varies depending on the dataset and objects. To increase the detection of true positive events, we also evaluated the pair-wise combination of these methods to suggest the best combination for sensitivity and F-measure. In addition, we compared the performance of the tools for the identification of three classes (read-through, inter-chromosomal and intra-others) of chimeric RNAs. Finally, we performed TOPSIS analyses and ranked the weighted performance of the 16 tools.

LongGF: computational algorithm and software tool for fast and accurate detection of gene fusions by long-read transcriptome sequencing

BACKGROUND

Long-read RNA-Seq techniques can generate reads that encompass a large proportion or the entire mRNA/cDNA molecules, so they are expected to address inherited limitations of short-read RNA-Seq techniques that typically generate < 150 bp reads. However, there is a general lack of software tools for gene fusion detection from long-read RNA-seq data, which takes into account the high basecalling error rates and the presence of alignment errors.

RESULTS

In this study, we developed a fast computational tool, LongGF, to efficiently detect candidate gene fusions from long-read RNA-seq data, including cDNA sequencing data and direct mRNA sequencing data. We evaluated LongGF on tens of simulated long-read RNA-seq datasets, and demonstrated its superior performance in gene fusion detection. We also tested LongGF on a Nanopore direct mRNA sequencing dataset and a PacBio sequencing dataset generated on a mixture of 10 cancer cell lines, and found that LongGF achieved better performance to detect known gene fusions over existing computational tools. Furthermore, we tested LongGF on a Nanopore cDNA sequencing dataset on acute myeloid leukemia, and pinpointed the exact location of a translocation (previously known in cytogenetic resolution) in base resolution, which was further validated by Sanger sequencing.

CONCLUSIONS

In summary, LongGF will greatly facilitate the discovery of candidate gene fusion events from long-read RNA-Seq data, especially in cancer samples. LongGF is implemented in C++ and is available at https://github.com/WGLab/LongGF .

NDRG1-PLAG1 and TRPS1-PLAG1 Fusion Genes in Chondroid Syringoma

BACKGROUND/AIM

Chondroid syringoma is a rare benign tumor emanating from sweat glands. Although rearrangements of the pleomorphic adenoma gene 1 (PLAG1) have been reported in such tumors, information on PLAG1 fusion genes is very limited.

MATERIALS AND METHODS

Cytogenetic, fluorescence in situ hybridization, RNA sequencing, array comparative genomic hybridization, reverse transcription polymerase chain reaction, and Sanger sequencing analyses were performed on two chondroid syringoma cases.

RESULTS

Both tumors had structural rearrangements of chromosome 8. An NDRG1-PLAG1 transcript was found in the first tumor in which exon 3 of PLAG1 was fused with exon 1 of NDRG1. A TRPS1-PLAG1 chimeric transcript was detected in the second chondroid syringoma in which exon 2 or exon 3 of PLAG1 was fused with exon 1 of TRPS1.

CONCLUSION

The NDRG1-PLAG1 and TRPS1-PLAG1 resemble other PLAG1 fusion genes inasmuch as the expression of PLAG1 comes under the control of the NDRG1 or TRPS1 promoter.

Cytogenetic and molecular genetic methods for chromosomal translocations detection with reference to the KMT2A/MLL gene

Acute leukemias (ALs) are often associated with chromosomal translocations, in particular, KMT2A/MLL gene rearrangements. Identification or confirmation of these translocations is carried out by a number of genetic and molecular methods, some of which are routinely used in clinical practice, while others are primarily used for research purposes. In the clinic, these methods serve to clarify diagnoses and monitor the course of disease and therapy. On the other hand, the identification of new translocations and the confirmation of known translocations are of key importance in the study of disease mechanisms and further molecular classification. There are multiple methods for the detection of rearrangements that differ in their principle of operation, the type of problem being solved, and the cost-result ratio. This review is intended to help researchers and clinicians studying AL and related chromosomal translocations to navigate this variety of methods. All methods considered in the review are grouped by their principle of action and include karyotyping, fluorescence in situ hybridization (FISH) with probes for whole chromosomes or individual loci, PCR and reverse transcription-based methods, and high-throughput sequencing. Another characteristic of the described methods is the type of problem being solved. This can be the discovery of new rearrangements, the determination of unknown partner genes participating in the rearrangement, or the confirmation of the proposed rearrangement between the two genes. We consider the specifics of the application, the basic principle of each method, and its pros and cons. To illustrate the application, examples of studying the rearrangements of the KMT2A/MLL gene, one of the genes that are often rearranged in AL, are mentioned.

Next-generation diagnostics for precision oncology: Preanalytical considerations, technical challenges, and available technologies

Molecular diagnostics as the centrepiece of precision oncology has gone through revolutionary developments over the last decade, becoming tremendously broad, deep and precise with still ongoing advancements. In the majority of scenarios, treatment selection for cancer patients without any type of molecular characterization is no longer conceivable. Considering the impact of sample quality on the reliability of molecular analyses and the importance of the results for the fate of an individual patient, it is surprising how sparsely preanalytical and analytical requirements are addressed scientifically. Standardization and rigorous quality assessment continue to play only a marginal role in the field. Within this review, we will systematically discuss influencing preanalytic parameters and technology setups affecting molecular test results. We will shed light on the specifics of different analytes, technical modalities, and analysis pipelines. The review will have a certain focus on broad molecular genetic tumour testing with next generation sequencing but will go beyond that including other molecular diagnostic modalities and will give a glimpse into the future of molecular testing.

Transcriptomics and solid tumors: The next frontier in precision cancer medicine

Transcriptomics, which encompasses assessments of alternative splicing and alternative polyadenylation, identification of fusion transcripts, explorations of noncoding RNAs, transcript annotation, and discovery of novel transcripts, is a valuable tool for understanding cancer mechanisms and identifying biomarkers. Recent advances in high-throughput technologies have enabled large-scale gene expression profiling. Importantly, RNA expression profiling of tumor tissue has been successfully used to determine clinically actionable molecular alterations. The WINTHER precision medicine clinical trial was the first prospective trial in diverse solid malignancies that assessed both genomics and transcriptomics to match treatments to specific molecular alterations. The use of transcriptome analysis in WINTHER and other trials increased the number of targetable -omic changes compared to genomic profiling alone. Other applications of transcriptomics involve the evaluation of tumor and circulating noncoding RNAs as predictive and prognostic biomarkers, the improvement of risk stratification by the use of prognostic and predictive multigene assays, the identification of fusion transcripts that drive tumors, and an improved understanding of the impact of DNA changes as some genomic alterations are silenced at the RNA level. Finally, RNA sequencing and gene expression analysis have been incorporated into clinical trials to identify markers predicting response to immunotherapy. Many issues regarding the complexity of the analysis, its reproducibility and variability, and the interpretation of the results still need to be addressed. The integration of transcriptomics with genomics, proteomics, epigenetics, and tumor immune profiling will improve biomarker discovery and our understanding of disease mechanisms and, thereby, accelerate the implementation of precision oncology.

Using equivalence class counts for fast and accurate testing of differential transcript usage

Background: RNA sequencing has enabled high-throughput and fine-grained quantitative analyses of the transcriptome. While differential gene expression is the most widely used application of this technology, RNA-seq data also has the resolution to infer differential transcript usage (DTU), which can elucidate the role of different transcript isoforms between experimental conditions, cell types or tissues. DTU has typically been inferred from exon-count data, which has issues with assigning reads unambiguously to counting bins, and requires alignment of reads to the genome. Recently, approaches have emerged that use transcript quantification estimates directly for DTU. Transcript counts can be inferred from 'pseudo' or lightweight aligners, which are significantly faster than traditional genome alignment. However, recent evaluations show lower sensitivity in DTU analysis compared to exon-level analysis. Transcript abundances are estimated from equivalence classes (ECs), which determine the transcripts that any given read is compatible with. Recent work has proposed performing a variety of RNA-seq analysis directly on equivalence class counts (ECCs). Methods: Here we demonstrate that ECCs can be used effectively with existing count-based methods for detecting DTU. We evaluate this approach on simulated human and drosophila data, as well as on a real dataset through subset testing. Results: We find that ECCs have similar sensitivity and false discovery rates as exon-level counts but can be generated in a fraction of the time through the use of pseudo-aligners. Conclusions: We posit that equivalence class read counts are a natural unit on which to perform differential transcript usage analysis.

Fast and accurate differential transcript usage by testing equivalence class counts

Background: RNA sequencing has enabled high-throughput and fine-grained quantitative analyses of the transcriptome. While differential gene expression is the most widely used application of this technology, RNA-seq data also has the resolution to infer differential transcript usage (DTU), which can elucidate the role of different transcript isoforms between experimental conditions, cell types or tissues. DTU has typically been inferred from exon-count data, which has issues with assigning reads unambiguously to counting bins, and requires alignment of reads to the genome. Recently, approaches have emerged that use transcript quantifications estimates directly for DTU. Transcript counts can be inferred from 'pseudo' or lightweight aligners, which are significantly faster than traditional genome alignment. However, recent evaluations show lower sensitivity in DTU analysis. Transcript abundances are estimated from equivalence classes (ECs), which determine the transcripts that any given read is compatible with. Recent work has proposed performing differential expression testing directly on equivalence class read counts (ECs). Methods: Here we demonstrate that ECs can be used effectively with existing count-based methods for detecting DTU. We evaluate this approach on simulated human and drosophila data, as well as on a real dataset through subset testing. Results: We find that ECs counts have similar sensitivity and false discovery rates as exon-level counts but can be generated in a fraction of the time through the use of pseudo-aligners. Conclusions: We posit that equivalence class read counts are a natural unit on which to perform many types of analysis.