Survey of 548 oncogenic fusion transcripts in thyroid tumors supports the importance of the already established thyroid fusions genes

Exposure-Response Relationships for Pralsetinib in Patients with RET-Altered Thyroid Cancer or RET Fusion-Positive Nonsmall Cell Lung Cancer

Pralsetinib is a highly potent oral kinase inhibitor of oncogenic RET (rearranged during transfection) fusions and mutations. Pralsetinib received approval from the United States Food and Drug Administration for the treatment of patients with metastatic RET fusion-positive non-small cell lung cancer (NSCLC), and received accelerated approval for the treatment of patients with RET fusion-positive thyroid cancer. Exposure-response (ER) analyses of efficacy were performed separately in patients with thyroid cancer and in patients with NSCLC, but data for all patients were pooled for the safety analysis. ER models were developed with time-varying exposure; the effect of covariates was also examined. For patients with NSCLC, a higher starting dose was associated with improved progression-free survival (PFS), but this improvement did not correlate with a higher exposure overall. Significant covariates included sex and baseline Eastern Cooperative Oncology Group (ECOG) score. For patients with thyroid cancer, a higher exposure was associated with improved PFS. Significant covariates included prior systemic cancer therapy and ECOG score. For safety, higher exposure was associated with a greater risk of grade ≥3 anemia, pneumonia, and lymphopenia. Patients with an ECOG score of ≥1 had an increased risk of grade ≥3 pneumonia. Non-White patients had a lower risk of grade ≥3 lymphopenia. ER analysis revealed that higher pralsetinib exposure was associated with improved PFS in thyroid cancer, but not in NSCLC. However, a higher starting dose (ie, 400 vs ≤300 mg daily) was correlated with better PFS for all indications. Higher exposure was also associated with an increased risk of grade ≥3 adverse events (AEs); however, the overall incidence of these events was acceptably low (≤20%). This analysis supports the use of a 400 mg starting dose of pralsetinib, allowing for dose reduction in the event of AEs.

Identification of fusions with potential clinical significance in melanoma

Though uncommon in melanoma, gene fusions may have therapeutic implications. Next generation sequencing-based clinical assays, designed to detect relevant gene fusions, mutations, and copy number changes, were performed on 750 melanomas (375 primary and 375 metastases) at our institution from 2014-2021. These included 599 (80%) cutaneous, 38 (5%) acral, 11 (1.5%) anorectal, 23 (3%) sinonasal, 27 (3.6%) eye (uveal/ conjunctiva), 11 (1.5%) genital (vulva/penile), and 41 (5.5%) melanomas of unknown primary. Sixteen fusions (2%) were detected in samples from 16 patients: 12/599 (2%) cutaneous, 2/38 (5%) acral, 1/9 (11%) vulva, 1/23(4.3%) sinonasal; and 12/16 (75%) fusions were potentially targetable. We identified two novel rearrangements: NAGS::MAST2 and NOTCH1::GNB1; and two fusions that have been reported in other malignancies but not in melanoma: CANT1::ETV4 (prostate cancer) and CCDC6::RET (thyroid cancer). Additional fusions, previously reported in melanoma, included: EML4::ALK, MLPH::ALK, AGAP3::BRAF, AGK::BRAF, CDH3::BRAF, CCT8::BRAF, DIP2B::BRAF, EFNB1::RAF1, LRCH3::RAF1, MAP4::RAF1, RUFY1::RAF1, and ADCY2::TERT. Fusion positive melanomas harbored recurrent alterations in TERT and CDKN2A, among others. Gene fusions were exceedingly rare (0.2%) in BRAF/RAS/NF1-mutant tumors and were detected in 5.6% of triple wild-type melanomas. Interestingly, gene rearrangements were significantly enriched within the subset of triple wild-type melanomas that harbor TERT promoter mutations (18% versus 2%, p < 0.0001). Thirteen (81%) patients were treated with immunotherapy for metastatic disease or in the adjuvant setting. Six of 12 (50%) patients with potentially actionable fusions progressed on immunotherapy, and 3/6 (50%) were treated with targeted agents (ALK and MEK inhibitors), 2 off-label and 1 as part of a clinical trial. One patient with an AGAP3::BRAF fusion positive melanoma experienced a 30-month long response to trametinib. We show that, detecting fusions, especially in triple wild-type melanomas with TERT promoter mutations, may have a clinically significant impact in patients with advanced disease who have failed front-line immunotherapy.

Development of a Low-cost NGS Test for the Evaluation of Thyroid Nodules

Ultrasound-guided fine needle aspiration cytology (FNAC) is the preferred method of identifying malignancy in palpable thyroid nodules using the Bethesda reporting system. However, in around 30-40% of FNACs (Bethesda categories III, IV, and V), the results are indeterminate and surgery is required to confirm malignancy. Out of those who undergo surgery, only 10-40% of patients in these categories are found to have malignancies, thus proving surgery to be unnecessary for some patients or to be incomplete in others. While molecular testing on thyroid FNAC material is part of the American Thyroid Association (ATA) guidelines in evaluating thyroid nodules, it is currently unavailable in India due to cost constraints. In this study, we prospectively collected FNAC samples from sixty-nine patients who presented with palpable thyroid nodules. We designed a cost-effective next-generation sequencing (NGS) test to query multiple variants in the DNA and RNA isolated from the fine needle aspirate. The identification of oncogenic variants was considered to be indicative of malignancy, and confirmed by surgical histopathology. The panel showed an overall sensitivity of 81.25% and a specificity of 100%, while in the case of Bethesda categories III, IV, and V, the sensitivity was higher (87.5%) and the specificity was established at 100%. The panel could thereby serve as a rule-in test for the diagnosis of thyroid cancer and therefore help identify patients who require surgery, especially in the indeterminate Bethesda categories III, IV, and V.

Cerebrospinal Fluid Penetration and Combination Therapy of Entrectinib for Disseminated ROS1/NTRK-Fusion Positive Pediatric High-Grade Glioma

Targeting oncogenic fusion-genes in pediatric high-grade gliomas (pHGG) with entrectinib has emerged as a highly promising therapeutic approach. Despite ongoing clinical studies, to date, no reports on the treatment of cerebrospinal fluid (CSF) disseminated fusion-positive pHGG exist. Moreover, clinically important information of combination with other treatment modalities such as intrathecal therapy, radiotherapy and other targeted agents is missing. We report on our clinical experience of entrectinib therapy in two CSF disseminated ROS1/NTRK-fusion-positive pHGG cases. Combination of entrectinib with radiotherapy or intrathecal chemotherapy appears to be safe and has the potential to act synergistically with entrectinib treatment. In addition, we demonstrate CSF penetrance of entrectinib for the first time in patient samples suggesting target engagement even upon CSF dissemination. Moreover, in vitro analyses of two novel cell models derived from one case with NTRK-fusion revealed that combination therapy with either a MEK (trametinib) or a CDK4/6 (abemaciclib) inhibitor synergistically enhances entrectinib anticancer effects. In summary, our comprehensive study, including clinical experience, CSF penetrance and in vitro data on entrectinib therapy of NTRK/ROS1-fusion-positive pHGG, provides essential clinical and preclinical insights into the multimodal treatment of these highly aggressive tumors. Our data suggest that combined inhibition of NTRK/ROS1 and other therapeutic vulnerabilities enhances the antitumor effect, which should be followed-up in further preclinical and clinical studies.

Genetic differences in follicular thyroid carcinoma between Asian and Western countries: a systematic review

Thyroid cancer is the most common endocrine malignancy, and follicular thyroid carcinoma (FTC) is the second most common thyroid cancer following papillary thyroid carcinoma (PTC). RAS mutation and PAX8/PPARγ rearrangement are the two representative genetic alterations in FTC, and there are studies from various countries on their regional frequencies. In this study, we systematically reviewed all available literature aiming to create a complete global map showing the frequencies of these common oncogenic drivers in FTC and to highlight the trends in Asian and Western countries. We performed a search in two electronic databases and identified 71 studies that fit our criteria from 1,329 studies found with our database search terms. There were 54 articles with 1,143 FTC patients and 39 articles with 764 FTC patients available for calculating the frequency of RAS mutation and PAX8/PPARγ rearrangement, respectively. NRAS mutation was the most frequent RAS mutation in all regions, followed by HRAS and KRAS mutation. The frequency of RAS mutation in Asian countries was higher than Western countries (34% vs. 27%, P=0.006) when the mutation detection method was not taken into account. In contrast, this difference in RAS mutation incidence between Asian and Western countries (28% vs. 25%, P=0.47) did not show up in our subgroup analysis incorporating only studies using direct sequencing method. The reported difference of RAS mutation frequency in the previous literature might not be due to the true prevalence of RAS mutation. They could be attributed to the difference in the detection method. As to PAX8/PPARγ rearrangement, Western countries overall had a much higher prevalence than Asian countries (23% vs. 4%, P<0.001), but some European countries had a low incidence, implying regional heterogeneity of PAX8/PPARγ rearrangement. A substantial lack of mutation data in FTC was found in several regions of the world such as Central Asia, Middle East, Africa, and Central and South America. Our results provide the most comprehensive global status of representative genetic alterations in FTC and highlight the similarities and differences between Asian and Western countries.

Pan-Cancer Analysis Reveals the Diverse Landscape of Novel Sense and Antisense Fusion Transcripts

Gene fusions that contribute to oncogenicity can be explored for identifying cancer biomarkers and potential drug targets. To investigate the nature and distribution of fusion transcripts in cancer, we examined the transcriptome data of about 9,000 primary tumors from 33 different cancers in TCGA (The Cancer Genome Atlas) along with cell line data from CCLE (Cancer Cell Line Encyclopedia) using ChimeRScope, a novel fusion detection algorithm. We identified several fusions with sense (canonical, 39%) or antisense (non-canonical, 61%) transcripts recurrent across cancers. The majority of the recurrent non-canonical fusions found in our study are novel, unexplored, and exhibited highly variable profiles across cancers, with breast cancer and glioblastoma having the highest and lowest rates, respectively. Overall, 4,344 recurrent fusions were identified from TCGA in this study, of which 70% were novel. Additional analysis of 802 tumor-derived cell line transcriptome data across 20 cancers revealed significant variability in recurrent fusion profiles between primary tumors and corresponding cell lines. A subset of canonical and non-canonical fusions was validated by examining the structural variation evidence in whole-genome sequencing (WGS) data or by Sanger sequencing of fusion junctions. Several recurrent fusion genes identified in our study show promise for drug repurposing in basket trials and present opportunities for mechanistic studies.

Targetable Gene Fusions Associate With the IDH Wild-Type Astrocytic Lineage in Adult Gliomas

Gene fusions involving oncogenes have been reported in gliomas and may serve as novel therapeutic targets. Using RNA-sequencing, we interrogated a large cohort of gliomas to assess for the incidence of targetable genetic fusions. Gliomas (n = 390) were profiled using the ArcherDx FusionPlex Assay. Fifty-two gene targets were analyzed and fusions with preserved kinase domains were investigated. Overall, 36 gliomas (9%) harbored a total of 37 potentially targetable fusions, the majority of which were found in astrocytomas (n = 34). Within this lineage 11% (25/235) of glioblastomas, 12% (5/42) of anaplastic astrocytomas, 8% (2/25) of grade II astrocytomas, and 33% (2/6) of pilocytic astrocytoma harbored targetable fusions. Fusions were significantly more frequent in IDH wild-type tumors (12%, n = 31/261) relative to IDH mutants (4%; n = 4/109) (p = 0.011). No fusions were seen in oligodendrogliomas. The most frequently observed therapeutically targetable fusions were in FGFR (n = 12), MET (n = 11), and NTRK (n = 8). Several additional novel fusions that have not been previously described in gliomas were identified including EGFR:VWC2 and FGFR3:NBR1. In summary, targetable gene fusions are enriched in IDH wild-type high-grade astrocytic tumors, which will influence enrollment in and interpretation of clinical trials of glioma patients.

CRABP1, C1QL1 and LCN2 are biomarkers of differentiated thyroid carcinoma, and predict extrathyroidal extension

Background

The prognostic variability of thyroid carcinomas has led to the search for accurate biomarkers at the molecular level. Follicular thyroid carcinoma (FTC) is a typical example of differentiated thyroid carcinomas (DTC) in which challenges are faced in the differential diagnosis.

Methods

We used high-throughput paired-end RNA sequencing technology to study four cases of FTC with different degree of capsular invasion: two minimally invasive (mFTC) and two widely invasive FTC (wFTC). We searched by genes differentially expressed between mFTC and wFTC, in an attempt to find biomarkers of thyroid cancer diagnosis and/or progression. Selected biomarkers were validated by real-time quantitative PCR in 137 frozen thyroid samples and in an independent dataset (TCGA), evaluating the diagnostic and the prognostic performance of the candidate biomarkers.

Results

We identified 17 genes significantly differentially expressed between mFTC and wFTC. C1QL1, LCN2, CRABP1 and CILP were differentially expressed in DTC in comparison with normal thyroid tissues. LCN2 and CRABP1 were also differentially expressed in DTC when compared with follicular thyroid adenoma. Additionally, overexpression of LCN2 and C1QL1 were found to be independent predictors of extrathyroidal extension in DTC.

Conclusions

We conclude that the underexpression of CRABP1 and the overexpression of LCN2 may be useful diagnostic biomarkers in thyroid tumours with questionable malignity, and the overexpression of LCN2 and C1QL1 may be useful for prognostic purposes.

Electronic supplementary material

The online version of this article (10.1186/s12885-017-3948-3) contains supplementary material, which is available to authorized users.

RAS Mutations, and RET/PTC and PAX8/PPAR-gamma Chromosomal Rearrangements Are Also Prevalent in Benign Thyroid Lesions: Implications Thereof and A Systematic Review

BACKGROUND

Molecular markers associated with thyroid malignancy are increasingly being used as differential diagnostic tools for thyroid nodules. However, little has been reported recently regarding the prevalence of these markers in benign lesions. The literature was systematically reviewed to examine studies that reported on the prevalence of these markers in benign thyroid lesions.

METHODS

Appropriate studies published between January 1, 2000, and April 30, 2015, and cataloged in PubMed, Embase, Cochrane, Scopus, and Web of Science databases were searched for by combining different keywords for "thyroid tumor" with both general and specific keywords for "molecular marker" by using "AND" as the Boolean operator. All studies meeting criteria that reported the prevalence of RAS mutations, and RET/PTC and PAX8/PPAR-gamma chromosomal rearrangements in benign thyroid lesions were included for study.

RESULTS

A total of 64 articles (including 8162 patients, of whom 42.5% had benign lesions) that met all the study criteria were systematically reviewed and abstracted. Among 35 studies examining RAS mutations, the reported prevalence of RAS mutation in benign lesions ranged from 0% to 48%. In 38 studies examining RET/PTC rearrangements, the prevalence in benign lesions ranged from 0% to 68%. PAX8/PPAR-gamma rearrangements were examined in 27 studies, with the prevalence in benign lesions ranging from 0% to 55%.

CONCLUSION

The presence of these biomarkers and the tremendous variation in reports of their prevalence in benign lesions suggests the need for caution when including these markers in diagnostic decisions. Further understanding of the importance of these markers, as well as newly discovered markers of thyroid malignancy, may be required in order to avoid overtreatment of patients with benign thyroid tumors.

Global analysis of somatic structural genomic alterations and their impact on gene expression in diverse human cancers

Tumor genomes are mosaics of somatic structural variants (SVs) that may contribute to the activation of oncogenes or inactivation of tumor suppressors, for example, by altering gene copy number amplitude. However, there are multiple other ways in which SVs can modulate transcription, but the general impact of such events on tumor transcriptional output has not been systematically determined. Here we use whole-genome sequencing data to map SVs across 600 tumors and 18 cancers, and investigate the relationship between SVs, copy number alterations (CNAs), and mRNA expression. We find that 34% of CNA breakpoints can be clarified structurally and that most amplifications are due to tandem duplications. We observe frequent swapping of strong and weak promoters in the context of gene fusions, and find that this has a measurable global impact on mRNA levels. Interestingly, several long noncoding RNAs were strongly activated by this mechanism. Additionally, SVs were confirmed in telomere reverse transcriptase (TERT) upstream regions in several cancers, associated with elevated TERT mRNA levels. We also highlight high-confidence gene fusions supported by both genomic and transcriptomic evidence, including a previously undescribed paired box 8 (PAX8)-nuclear factor, erythroid 2 like 2 (NFE2L2) fusion in thyroid carcinoma. In summary, we combine SV, CNA, and expression data to provide insights into the structural basis of CNAs as well as the impact of SVs on gene expression in tumors.

The landscape of kinase fusions in cancer

Human cancer genomes harbour a variety of alterations leading to the deregulation of key pathways in tumour cells. The genomic characterization of tumours has uncovered numerous genes recurrently mutated, deleted or amplified, but gene fusions have not been characterized as extensively. Here we develop heuristics for reliably detecting gene fusion events in RNA-seq data and apply them to nearly 7,000 samples from The Cancer Genome Atlas. We thereby are able to discover several novel and recurrent fusions involving kinases. These findings have immediate clinical implications and expand the therapeutic options for cancer patients, as approved or exploratory drugs exist for many of these kinases.

Kinases activated by gene fusions represent potentially important targets for the development of cancer drugs. Here, the authors develop a method for detecting gene fusion events in RNA sequencing data from The Cancer Genome Atlas and identify several novel recurrent fusions involving kinases.

RNA sequencing identifies multiple fusion transcripts, differentially expressed genes, and reduced expression of immune function genes in BRAF (V600E) mutant vs BRAF wild-type papillary thyroid carcinoma

CONTEXT

The BRAF V600E mutation (BRAF-MUT) confers an aggressive phenotype in papillary thyroid carcinoma, but unidentified additional genomic abnormalities may be required for full phenotypic expression.

OBJECTIVE

RNA sequencing (RNA-Seq) was performed to identify genes differentially expressed between BRAF-MUT and BRAF wild-type (BRAF-WT) tumors and to correlate changes to patient clinical status.

DESIGN

BRAF-MUT and BRAF-WT tumors were identified in patients with T1N0 and T2-3N1 tumors evaluated in a referral medical center. Gene expression levels were determined (RNA-Seq) and fusion transcripts were detected. Multiplexed capture/detection and digital counting of mRNA transcripts (nCounter, NanoString Technologies) validated RNA-Seq data for immune system-related genes.

PATIENTS

BRAF-MUT patients included nine women, three men; nine were TNM stage I and three were stage III. Three (25%) had tumor infiltrating lymphocytes. BRAF-WT included five women, three men; all were stage I, and five (62.5%) had tumor infiltrating lymphocytes.

RESULTS

RNA-Seq identified 560 of 13 085 genes differentially expressed between BRAF-MUT and BRAF-WT tumors. Approximately 10% of these genes were related to MetaCore immune function pathways; 51 were underexpressed in BRAF-MUT tumors, whereas 4 (HLAG, CXCL14, TIMP1, IL1RAP) were overexpressed. The four most differentially overexpressed immune genes in BRAF-WT tumors (IL1B; CCL19; CCL21; CXCR4) correlated with lymphocyte infiltration. nCounter confirmed the RNA-Seq expression level data. Eleven different high-confidence fusion transcripts were detected (four interchromosomal; seven intrachromosomal) in 13 of 20 tumors. All in-frame fusions were validated by RT-PCR.

CONCLUSION

BRAF-MUT papillary thyroid cancers have reduced expression of immune/inflammatory response genes compared with BRAF-WT tumors and correlate with lymphocyte infiltration. In contrast, HLA-G and CXCL14 are overexpressed in BRAF-MUT tumors. Sixty-five percent of tumors had between one and three fusion transcripts. Functional studies will be required to determine the potential role of these newly identified genomic abnormalities in contributing to the aggressiveness of BRAF-MUT and BRAF-WT tumors.

Exploration of the gene fusion landscape of glioblastoma using transcriptome sequencing and copy number data

BACKGROUND

RNA-seq has spurred important gene fusion discoveries in a number of different cancers, including lung, prostate, breast, brain, thyroid and bladder carcinomas. Gene fusion discovery can potentially lead to the development of novel treatments that target the underlying genetic abnormalities.

RESULTS

In this study, we provide comprehensive view of gene fusion landscape in 185 glioblastoma multiforme patients from two independent cohorts. Fusions occur in approximately 30-50% of GBM patient samples. In the Ivy Center cohort of 24 patients, 33% of samples harbored fusions that were validated by qPCR and Sanger sequencing. We were able to identify high-confidence gene fusions from RNA-seq data in 53% of the samples in a TCGA cohort of 161 patients. We identified 13 cases (8%) with fusions retaining a tyrosine kinase domain in the TCGA cohort and one case in the Ivy Center cohort. Ours is the first study to describe recurrent fusions involving non-coding genes. Genomic locations 7p11 and 12q14-15 harbor majority of the fusions. Fusions on 7p11 are formed in focally amplified EGFR locus whereas 12q14-15 fusions are formed by complex genomic rearrangements. All the fusions detected in this study can be further visualized and analyzed using our website: http://ivygap.swedish.org/fusions.

CONCLUSIONS

Our study highlights the prevalence of gene fusions as one of the major genomic abnormalities in GBM. The majority of the fusions are private fusions, and a minority of these recur with low frequency. A small subset of patients with fusions of receptor tyrosine kinases can benefit from existing FDA approved drugs and drugs available in various clinical trials. Due to the low frequency and rarity of clinically relevant fusions, RNA-seq of GBM patient samples will be a vital tool for the identification of patient-specific fusions that can drive personalized therapy.

Possible formation of mitochondrial-RNA containing chimeric or trimeric RNA implies a post-transcriptional and post-splicing mechanism for RNA fusion

Human cells are known to express many chimeric RNAs, i.e. RNAs containing two genes' sequences. Wondering whether there also is trimeric RNA, i.e. an RNA containing three genes' sequences, we wrote simple computer code to screen human expression sequence tags (ESTs) deposited in different public databases, and obtained hundreds of putative trimeric ESTs. We then used NCBI Blast and UCSC Blat browsers to further analyze their sequences, and identified 61 trimeric and two tetrameric ESTs (one EST containing four different sequences). We also identified 57 chimeric, trimeric or teterameric ESTs that contained both mitochondrial (mt) RNA and nuclear RNA (nRNA), i.e. were mtRNA-nRNA fusions. In some trimeric ESTs, the downstream partner was fused to the poly-A tail of the upstream partner, which, together with the mtRNA-nRNA fusions, suggests a possible new mechanism for RNA fusion that occurs after both transcription and splicing have been terminated, and possibly outside the nucleus, in contrast to the two current hypothetical mechanisms, trans-splicing and transcriptional-slippage, that occur in the nucleus. The mt-sequences in the mtRNA-nRNA fusions had pseudogenes in the nucleus but, surprisingly, localized mainly in chromosomes 1 and 5. In some mtRNA-nRNA fusions, as well as in some ESTs that were derived only from mtRNA, the mt-sequences might be cis- or trans-spliced. Actually, we cloned a new cis-spliced mtRNA, coined as 16SrRNA-s. Hence, mtDNA may not always be intron-less. Fusion of three or more RNAs to one, fusion of nRNA to mtRNA, and cis- or trans-splicing of mtRNA should all enlarge the cellular RNA repertoire, in turn enlarging the cellular functions. Therefore, future experimental verification of the existence of these novel classes of fusion RNAs and spliced mtRNAs in human cells should significantly advance our understanding of biology and medicine.

New methods as alternative or corrective measures for the pitfalls and artifacts of reverse transcription and polymerase chain reactions (RT-PCR) in cloning chimeric or antisense-accompanied RNA

We established new methods for cloning cDNA ends that start with reverse transcription (RT) and soon proceed with the synthesis of the second cDNA strand, avoiding manipulations of fragile RNA. Our 3′-end cloning method does not involve poly-dT primers and polymerase chain reactions (PCR), is low in efficiency but high in fidelity and can clone those RNAs without a poly-A tail. We also established a cDNA protection assay to supersede RNA protection assay. The protected cDNA can be amplified, cloned and sequenced, enhancing sensitivity and fidelity. We report that RT product using gene-specific primer (GSP) cannot be gene- or strand-specific because RNA sample contains endogenous random primers (ERP). The gene-specificity may be improved by adding a linker sequence at the 5′-end of the GSP to prime RT and using the linker as a primer in the ensuing PCR. The strand-specificity may be improved by using strand-specific DNA oligos in our protection assay. The CDK4 mRNA and TSPAN31 mRNA are transcribed from the opposite DNA strands and overlap at their 3′ ends. Using this relationship as a model, we found that the overlapped sequence might serve as a primer with its antisense as the template to create a wrong-template extension in RT or PCR. We infer that two unrelated RNAs or cDNAs overlapping at the 5′- or 3′-end might create a spurious chimera in this way, and many chimeras with a homologous sequence may be such artifacts. The ERP and overlapping antisense together set complex pitfalls, which one should be aware of.