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.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-14-818) contains supplementary material, which is available to authorized users.     

Identification of novel deregulated RNA metabolism-related genes in non-small cell lung cancer

Lung cancer is a leading cause of cancer death worldwide. Several alterations in RNA metabolism have been found in lung cancer cells; this suggests that RNA metabolism-related molecules are involved in the development of this pathology. In this study, we searched for RNA metabolism-related genes that exhibit different expression levels between normal and tumor lung tissues. We identified eight genes differentially expressed in lung adenocarcinoma microarray datasets. Of these, seven were up-regulated whereas one was down-regulated. Interestingly, most of these genes had not previously been associated with lung cancer. These genes play diverse roles in mRNA metabolism: three are associated with the spliceosome (ASCL3L1, SNRPB and SNRPE), whereas others participate in RNA-related processes such as translation (MARS and MRPL3), mRNA stability (PCBPC1), mRNA transport (RAE), or mRNA editing (ADAR2, also known as ADARB1). Moreover, we found a high incidence of loss of heterozygosity at chromosome 21q22.3, where the ADAR2 locus is located, in NSCLC cell lines and primary tissues, suggesting that the downregulation of ADAR2 in lung cancer is associated with specific genetic losses. Finally, in a series of adenocarcinoma patients, the expression of five of the deregulated genes (ADAR2, MARS, RAE, SNRPB and SNRPE) correlated with prognosis. Taken together, these results support the hypothesis that changes in RNA metabolism are involved in the pathogenesis of lung cancer, and identify new potential targets for the treatment of this disease.     

Identification of novel candidate genes by exome sequencing in Tunisian familial male breast cancer patients

Molecular Biology Reports - Male Breast Cancer (MBC) is a rare and aggressive disease that is associated with genetic factors. Mutations in BRCA1 and BRCA2 account for...     

Identification of novel mutations in FFPE lung adenocarcinomas using DEPArray sorting technology and next-generation sequencing

Formalin-fixed paraffin-embedded (FFPE) tissues are utilized as the standard diagnostic method in pathology laboratories. However, admixture of unwanted tissues and shortage of normal samples, which can be used to detect somatic mutation, are considered critical factors to accurately diagnose cancer. To explore these challenges, we sorted the pure tumor cells from 22 FFPE lung adenocarcinoma tissues via Di-Electro-Phoretic Array (DEPArray) technology, a new cell sorting technology, and analyzed the variants with next-generation sequencing (NGS) for the most accurate analysis. The allele frequencies of the all gene mutations were improved by 1.2 times in cells sorted via DEPArray (tumor suppressor genes, 1.3–10.1 times; oncogenes, 1.3–2.6 times). We identified 16 novel mutations using the sequencing from sorted cells via DEPArray technology, compared to detecting 4 novel mutation by the sequencing from unsorted cells. Using this analysis, we also revealed that five genes (TP53, EGFR, PTEN, RB1, KRAS, and CTNNB1) were somatically mutated in multiple homogeneous lung adenocarcinomas. Together, we sorted pure tumor cells from 22 FFPE lung adenocarcinomas by DEPArray technology and identified 16 novel somatic mutations. We also established the precise genomic landscape for more accurate diagnosis in 22 lung adenocarcinomas with mutations detected in pure tumor cells. The results obtained in this study could offer new avenues for the treatment and the diagnosis of squamous cell lung cancers.     

Assessment of metagenomic assembly using simulated next generation sequencing data

Due to the complexity of the protocols and a limited knowledge of the nature of microbial communities, simulating metagenomic sequences plays an important role in testing the performance of existing tools and data analysis methods with metagenomic data. We developed metagenomic read simulators with platform-specific (Sanger, pyrosequencing, Illumina) base-error models, and simulated metagenomes of differing community complexities. We first evaluated the effect of rigorous quality control on Illumina data. Although quality filtering removed a large proportion of the data, it greatly improved the accuracy and contig lengths of resulting assemblies. We then compared the quality-trimmed Illumina assemblies to those from Sanger and pyrosequencing. For the simple community (10 genomes) all sequencing technologies assembled a similar amount and accurately represented the expected functional composition. For the more complex community (100 genomes) Illumina produced the best assemblies and more correctly resembled the expected functional composition. For the most complex community (400 genomes) there was very little assembly of reads from any sequencing technology. However, due to the longer read length the Sanger reads still represented the overall functional composition reasonably well. We further examined the effect of scaffolding of contigs using paired-end Illumina reads. It dramatically increased contig lengths of the simple community and yielded minor improvements to the more complex communities. Although the increase in contig length was accompanied by increased chimericity, it resulted in more complete genes and a better characterization of the functional repertoire. The metagenomic simulators developed for this research are freely available.     

基于PacBio测序数据的蜜蜂球囊菌转录因子、融合基因及RNA编辑事件的鉴定

【目的】本研究旨在利用已获得的PacBio单分子实时(single molecule real-time, SMRT)测序数据对蜜蜂球囊菌Ascosphaera apis菌丝(AaM)和孢子(AaS)中的转录因子(TF)、融合基因和RNA编辑事件进行鉴定和分析,以期丰富蜜蜂球囊菌的相关信息,并为进一步探究它们的功能提供理论依据。【方法】利用BLASTx工具将AaM和AaS的全长转录本序列比对到Nr, Swiss-Prot和KEGG数据库以获得一致性最高的蛋白序列,再利用hmmscan软件将上述蛋白序列比对到Plant TFdb数据库以获得TF的分类及注释信息。采用TOFU软件中的fusion_finder.py程序进行融合基因的预测,进而分析融合基因的序列和位置信息。使用SAMtools预测AaM和AaS中的RNA编辑事件,再利用ANNOVAR软件对RNA编辑事件进行注释,进而采用相关生物信息学软件对RNA编辑位点基因进行GO功能和KEGG通路注释。【结果】在AaS中共鉴定到17个TF家族的213个TF,其中C2H2家族包含的TF成员最多。在AaM和AaS中分别鉴定到921和510个融合基因,二者共有的融合基因为510个,特有的融合基因分别为411和0个。在AaM和AaS中分别鉴定到547和191次RNA编辑事件,其中AaM中同义单核苷酸突变的数量最多,AaS中非同义单核苷酸突变的数量最多。此外,在AaM中鉴定到12种碱基替换类型,其中发生C->T的RNA编辑事件数量最多,达到158次;在AaS中鉴定到9种碱基替换类型,其中发生C->T和G->T的RNA编辑事件数量最多,均有42次。AaM和AaS中RNA编辑位点基因分别涉及19和24个GO功能条目;此外还能注释到11和20条KEGG通路。【结论】蜜蜂球囊菌的菌丝和孢子中含有丰富的TF、融合基因和RNA编辑位点;转录因子C2H2家族与蜜蜂球囊菌菌丝和孢子的生长发育和细胞活动具有潜在关联;RNA编辑事件的碱基替换类型在蜜蜂球囊菌和其他物种中具有物种特异性;RNA编辑可能在蜜蜂球囊菌菌丝和孢子的生长和代谢中发挥作用。