Explore the vast range of titles available.

RNA-sequencing (RNA-seq) has a wide variety of applications, but no single analysis pipeline can be used in all cases. We review all of the major steps in RNA-seq data analysis, including experimental design, quality control, read alignment, quantification of gene and transcript levels, visualization, differential gene expression, alternative splicing, functional analysis, gene fusion detection and eQTL mapping. We highlight the challenges associated with each step. We discuss the analysis of small RNAs and the integration of RNA-seq with other functional genomics techniques. Finally, we discuss the outlook for novel technologies that are changing the state of the art in transcriptomics.

Glioblastoma (GBM) is the most common type of glioma and is uniformly fatal. Currently, tumour heterogeneity and mutation acquisition are major impedances for tailoring personalized therapy. We collected blood and tumour tissue samples from 25 GBM patients and 25 blood samples from healthy controls. Cell‐free DNA (cfDNA) was extracted from the plasma of GBM patients and from healthy controls. Tumour DNA was extracted from fresh tumour samples. Extracted DNA was sequenced using a whole‐genome sequencing procedure. We also collected 180 tumour DNA datasets from GBM patients publicly available at the TCGA/PANCANCER project. These data were analysed for mutations and gene–gene fusions that could be potential druggable targets. We found that plasma cfDNA concentrations in GBM patients were significantly elevated (22.6 ± 5 ng·mL−1), as compared to healthy controls (1.4 ± 0.4 ng·mL−1) of the same average age. We identified unique mutations in the cfDNA and tumour DNA of each GBM patient, including some of the most frequently mutated genes in GBM according to the COSMIC database (TP53, 18.75%; EGFR, 37.5%; NF1, 12.5%; LRP1B, 25%; IRS4, 25%). Using our gene–gene fusion database, ChiTaRS 5.0, we identified gene–gene fusions in cfDNA and tumour DNA, such as KDR–PDGFRA and NCDN–PDGFRA, which correspond to previously reported alterations of PDGFRA in GBM (44% of all samples). Interestingly, the PDGFRA protein fusions can be targeted by tyrosine kinase inhibitors such as imatinib, sunitinib, and sorafenib. Moreover, we identified BCR–ABL1 (in 8% of patients), COL1A1–PDGFB (8%), NIN–PDGFRB (8%), and FGFR1–BCR (4%) in cfDNA of patients, which can be targeted by analogues of imatinib. ROS1 fusions (CEP85L–ROS1 and GOPC–ROS1), identified in 8% of patient cfDNA, might be targeted by crizotinib, entrectinib, or larotrectinib. Thus, our study suggests that integrated analysis of cfDNA plasma concentration, gene mutations, and gene–gene fusions can serve as a diagnostic modality for distinguishing GBM patients who may benefit from targeted therapy. These results open new avenues for precision medicine in GBM, using noninvasive liquid biopsy diagnostics to assess personalized patient profiles. Moreover, repeated detection of druggable targets over the course of the disease may provide real‐time information on the evolving molecular landscape of the tumour. Personalized therapy of patients with glioblastoma (GBM) is challenging owing to tumour heterogeneity. Here, we extracted and sequenced cell‐free DNA (cfDNA) from the plasma of 25 GBM patients and tumour DNA from fresh tumour samples. We found that cfDNA concentrations in the plasma of GBM patients were significantly elevated, as compared to healthy controls. Moreover, we identified unique mutations and gene–gene fusions in the cfDNA and tumour DNA of GBM patients, some of which could be therapeutically targeted by tyrosine kinase inhibitors.

Matthew Meyerson and colleagues report exome sequencing of solitary fibrous tumor (SFT), a rare mesenchymal tumor. They identify a NAB2 - STAT6 fusion in 55% of cases. Solitary fibrous tumors (SFTs) are rare mesenchymal tumors. Here, we describe the identification of a NAB2-STAT6 fusion from whole-exome sequencing of 17 SFTs. Analysis in 53 tumors confirmed the presence of 7 variants of this fusion transcript in 29 tumors (55%), representing a lower bound for fusion frequency at this locus and suggesting that the NAB2-STAT6 fusion is a distinct molecular feature of SFTs.

Arul Chinnaiyan and colleagues identify NAB2 - STAT6 fusions in 52 of 52 solitary fibrous tumor cases. Overexpression of this fusion induced cell proliferation, which could be suppressed by knockdown of EGR1 . A 44-year old woman with recurrent solitary fibrous tumor (SFT)/hemangiopericytoma was enrolled in a clinical sequencing program including whole-exome and transcriptome sequencing. A gene fusion of the transcriptional repressor NAB2 with the transcriptional activator STAT6 was detected. Transcriptome sequencing of 27 additional SFTs identified the presence of a NAB2 - STAT6 gene fusion in all tumors. Using RT-PCR and sequencing, we detected this fusion in all 51 SFTs, indicating high levels of recurrence. Expression of NAB2-STAT6 fusion proteins was confirmed in SFT, and the predicted fusion products harbor the early growth response (EGR)-binding domain of NAB2 fused to the activation domain of STAT6. Overexpression of the NAB2 - STAT6 gene fusion induced proliferation in cultured cells and activated the expression of EGR-responsive genes. These studies establish NAB2 - STAT6 as the defining driver mutation of SFT and provide an example of how neoplasia can be initiated by converting a transcriptional repressor of mitogenic pathways into a transcriptional activator.

ALK, ROS1, and RET kinase fusions are important predictive biomarkers of tyrosine kinase inhibitors (TKIs) in non‐small‐cell lung cancer (NSCLC). Analysis of cell‐free DNA (cfDNA) provides a noninvasive method to identify gene changes in tumor cells. The present study sought to use cfRNA and cfDNA for identifying fusion genes. A reliable protocol was established to detect fusion genes using cfRNA and assessed the analytical validity and clinical usefulness in 30 samples from 20 cases of fusion‐positive NSCLC. The results of cfRNA‐based assays were compared with tissue biopsy and cfDNA‐based liquid biopsy (Guardant360 plasma next‐generation sequencing [NGS] assay). The overall sensitivity of the cfRNA‐based assay was 26.7% (8/30) and that of cfDNA‐based assay was 16.7% (3/18). When analysis was limited to the samples collected at chemo‐naïve or progressive disease status and available for both assays, the sensitivity of the cfRNA‐based assay was 77.8% (7/9) and that of cfDNA‐based assay was 33.3% (3/9). Fusion gene identification in cfRNA was correlated with treatment response. These results suggest that the proposed cfRNA assay is a useful diagnostic test for patients with insufficient tissues to facilitate effective administration of first‐line treatment and is a useful tool to monitor the progression of NSCLC for consideration of second‐line treatments. cfRNA‐ and cfDNA‐based assays are evaluated in 20 cases of fusion‐positive NSCLC. cfRNA assay was superior to cfDNA assay for the detection of gene fusions. The results of the cfRNA assay were consistent with the therapeutic effect.

André Oliveira and colleagues identify a recurrent translocation in biphenotypic sinonasal sarcomas generating a PAX3 - MAML3 fusion gene. The resulting protein is a potent transcriptional activator of PAX3 response elements and is associated with aberrant expression of genes involved in neuroectodermal and myogenic differentiation. Biphenotypic sinonasal sarcoma (SNS) is a newly described tumor of the nasal and paranasal areas. Here we report a recurrent chromosomal translocation in SNS, t(2;4)(q35;q31.1), resulting in a PAX3-MAML3 fusion protein that is a potent transcriptional activator of PAX3 response elements. The SNS phenotype is characterized by aberrant expression of genes involved in neuroectodermal and myogenic differentiation, closely simulating the developmental roles of PAX3.

Abstract 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.

Thyroid cancer is a common endocrine malignancy that encompasses well-differentiated as well as dedifferentiated cancer types. The latter tumors have high mortality and lack effective therapies. Using a paired-end RNA-sequencing approach, we report the discovery of rearrangements involving the anaplastic lymphoma kinase (ALK) gene in thyroid cancer. The most common of these involves a fusion between ALK and the striatin (STRN) gene, which is the result of a complex rearrangement involving the short arm of chromosome 2. STRN-ALK leads to constitutive activation of ALK kinase via dimerization mediated by the coiled-coil domain of STRN and to a kinase-dependent, thyroid-stimulating hormone-independent proliferation of thyroid cells. Moreover, expression of STRN-ALK transforms cells in vitro and induces tumor formation in nude mice. The kinase activity of STRN-ALK and the ALK-induced cell growth can be blocked by the ALK inhibitors crizotinib and TAE684. In addition to well-differentiated papillary cancer, STRN-ALK was found with a higher prevalence in poorly differentiated and anaplastic thyroid cancers, and it did not overlap with other known driver mutations in these tumors. Our data demonstrate that STRN-ALK fusion occurs in a subset of patients with highly aggressive types of thyroid cancer and provide initial evidence suggesting that it may represent a therapeutic target for these patients.

Three subfamilies of grasses, the Ehrhartoideae, Panicoideae and Pooideae, provide the bulk of human nutrition and are poised to become major sources of renewable energy. Here we describe the genome sequence of the wild grass Brachypodium distachyon (Brachypodium), which is, to our knowledge, the first member of the Pooideae subfamily to be sequenced. Comparison of the Brachypodium, rice and sorghum genomes shows a precise history of genome evolution across a broad diversity of the grasses, and establishes a template for analysis of the large genomes of economically important pooid grasses such as wheat. The high-quality genome sequence, coupled with ease of cultivation and transformation, small size and rapid life cycle, will help Brachypodium reach its potential as an important model system for developing new energy and food crops.

Background Events of gene fusion have been reported in several organisms. However, the general role of gene fusion as part of new gene origination remains unknown. Results We conduct genome-wide interrogations of four Oryza genomes by designing and implementing novel pipelines to detect fusion genes. Based on the phylogeny of ten plant species, we detect 310 fusion genes across four Oryza species. The estimated rate of origination of fusion genes in the Oryza genus is as high as 63 fusion genes per species per million years, which is fixed at 16 fusion genes per species per million years and much higher than that in flies. By RNA sequencing analysis, we find more than 44% of the fusion genes are expressed and 90% of gene pairs show strong signals of purifying selection. Further analysis of CRISPR/Cas9 knockout lines indicates that newly formed fusion genes regulate phenotype traits including seed germination, shoot length and root length, suggesting the functional significance of these genes. Conclusions We detect new fusion genes that may drive phenotype evolution in Oryza. This study provides novel insights into the genome evolution of Oryza.