Cargando…

The “Grep” Command But Not FusionMap, FusionFinder or ChimeraScan Captures the CIC-DUX4 Fusion Gene from Whole Transcriptome Sequencing Data on a Small Round Cell Tumor with t(4;19)(q35;q13)

Whole transcriptome sequencing was used to study a small round cell tumor in which a t(4;19)(q35;q13) was part of the complex karyotype but where the initial reverse transcriptase PCR (RT-PCR) examination did not detect a CIC-DUX4 fusion transcript previously described as the crucial gene-level outc...

Descripción completa

Detalles Bibliográficos
Autores principales: Panagopoulos, Ioannis, Gorunova, Ludmila, Bjerkehagen, Bodil, Heim, Sverre
Formato: Online Artículo Texto
Lenguaje:English
Publicado: Public Library of Science 2014
Materias:
Acceso en línea:https://www.ncbi.nlm.nih.gov/pmc/articles/PMC4064965/
https://www.ncbi.nlm.nih.gov/pubmed/24950227
http://dx.doi.org/10.1371/journal.pone.0099439
_version_ 1782322003377651712
author Panagopoulos, Ioannis
Gorunova, Ludmila
Bjerkehagen, Bodil
Heim, Sverre
author_facet Panagopoulos, Ioannis
Gorunova, Ludmila
Bjerkehagen, Bodil
Heim, Sverre
author_sort Panagopoulos, Ioannis
collection PubMed
description Whole transcriptome sequencing was used to study a small round cell tumor in which a t(4;19)(q35;q13) was part of the complex karyotype but where the initial reverse transcriptase PCR (RT-PCR) examination did not detect a CIC-DUX4 fusion transcript previously described as the crucial gene-level outcome of this specific translocation. The RNA sequencing data were analysed using the FusionMap, FusionFinder, and ChimeraScan programs which are specifically designed to identify fusion genes. FusionMap, FusionFinder, and ChimeraScan identified 1017, 102, and 101 fusion transcripts, respectively, but CIC-DUX4 was not among them. Since the RNA sequencing data are in the fastq text-based format, we searched the files using the “grep” command-line utility. The “grep” command searches the text for specific expressions and displays, by default, the lines where matches occur. The “specific expression” was a sequence of 20 nucleotides from the coding part of the last exon 20 of CIC (Reference Sequence: NM_015125.3) chosen since all the so far reported CIC breakpoints have occurred here. Fifteen chimeric CIC-DUX4 cDNA sequences were captured and the fusion between the CIC and DUX4 genes was mapped precisely. New primer combinations were constructed based on these findings and were used together with a polymerase suitable for amplification of GC-rich DNA templates to amplify CIC-DUX4 cDNA fragments which had the same fusion point found with “grep”. In conclusion, FusionMap, FusionFinder, and ChimeraScan generated a plethora of fusion transcripts but did not detect the biologically important CIC-DUX4 chimeric transcript; they are generally useful but evidently suffer from imperfect both sensitivity and specificity. The “grep” command is an excellent tool to capture chimeric transcripts from RNA sequencing data when the pathological and/or cytogenetic information strongly indicates the presence of a specific fusion gene.
format Online
Article
Text
id pubmed-4064965
institution National Center for Biotechnology Information
language English
publishDate 2014
publisher Public Library of Science
record_format MEDLINE/PubMed
spelling pubmed-40649652014-06-25 The “Grep” Command But Not FusionMap, FusionFinder or ChimeraScan Captures the CIC-DUX4 Fusion Gene from Whole Transcriptome Sequencing Data on a Small Round Cell Tumor with t(4;19)(q35;q13) Panagopoulos, Ioannis Gorunova, Ludmila Bjerkehagen, Bodil Heim, Sverre PLoS One Research Article Whole transcriptome sequencing was used to study a small round cell tumor in which a t(4;19)(q35;q13) was part of the complex karyotype but where the initial reverse transcriptase PCR (RT-PCR) examination did not detect a CIC-DUX4 fusion transcript previously described as the crucial gene-level outcome of this specific translocation. The RNA sequencing data were analysed using the FusionMap, FusionFinder, and ChimeraScan programs which are specifically designed to identify fusion genes. FusionMap, FusionFinder, and ChimeraScan identified 1017, 102, and 101 fusion transcripts, respectively, but CIC-DUX4 was not among them. Since the RNA sequencing data are in the fastq text-based format, we searched the files using the “grep” command-line utility. The “grep” command searches the text for specific expressions and displays, by default, the lines where matches occur. The “specific expression” was a sequence of 20 nucleotides from the coding part of the last exon 20 of CIC (Reference Sequence: NM_015125.3) chosen since all the so far reported CIC breakpoints have occurred here. Fifteen chimeric CIC-DUX4 cDNA sequences were captured and the fusion between the CIC and DUX4 genes was mapped precisely. New primer combinations were constructed based on these findings and were used together with a polymerase suitable for amplification of GC-rich DNA templates to amplify CIC-DUX4 cDNA fragments which had the same fusion point found with “grep”. In conclusion, FusionMap, FusionFinder, and ChimeraScan generated a plethora of fusion transcripts but did not detect the biologically important CIC-DUX4 chimeric transcript; they are generally useful but evidently suffer from imperfect both sensitivity and specificity. The “grep” command is an excellent tool to capture chimeric transcripts from RNA sequencing data when the pathological and/or cytogenetic information strongly indicates the presence of a specific fusion gene. Public Library of Science 2014-06-20 /pmc/articles/PMC4064965/ /pubmed/24950227 http://dx.doi.org/10.1371/journal.pone.0099439 Text en © 2014 Panagopoulos et al http://creativecommons.org/licenses/by/4.0/ This is an open-access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are properly credited.
spellingShingle Research Article
Panagopoulos, Ioannis
Gorunova, Ludmila
Bjerkehagen, Bodil
Heim, Sverre
The “Grep” Command But Not FusionMap, FusionFinder or ChimeraScan Captures the CIC-DUX4 Fusion Gene from Whole Transcriptome Sequencing Data on a Small Round Cell Tumor with t(4;19)(q35;q13)
title The “Grep” Command But Not FusionMap, FusionFinder or ChimeraScan Captures the CIC-DUX4 Fusion Gene from Whole Transcriptome Sequencing Data on a Small Round Cell Tumor with t(4;19)(q35;q13)
title_full The “Grep” Command But Not FusionMap, FusionFinder or ChimeraScan Captures the CIC-DUX4 Fusion Gene from Whole Transcriptome Sequencing Data on a Small Round Cell Tumor with t(4;19)(q35;q13)
title_fullStr The “Grep” Command But Not FusionMap, FusionFinder or ChimeraScan Captures the CIC-DUX4 Fusion Gene from Whole Transcriptome Sequencing Data on a Small Round Cell Tumor with t(4;19)(q35;q13)
title_full_unstemmed The “Grep” Command But Not FusionMap, FusionFinder or ChimeraScan Captures the CIC-DUX4 Fusion Gene from Whole Transcriptome Sequencing Data on a Small Round Cell Tumor with t(4;19)(q35;q13)
title_short The “Grep” Command But Not FusionMap, FusionFinder or ChimeraScan Captures the CIC-DUX4 Fusion Gene from Whole Transcriptome Sequencing Data on a Small Round Cell Tumor with t(4;19)(q35;q13)
title_sort “grep” command but not fusionmap, fusionfinder or chimerascan captures the cic-dux4 fusion gene from whole transcriptome sequencing data on a small round cell tumor with t(4;19)(q35;q13)
topic Research Article
url https://www.ncbi.nlm.nih.gov/pmc/articles/PMC4064965/
https://www.ncbi.nlm.nih.gov/pubmed/24950227
http://dx.doi.org/10.1371/journal.pone.0099439
work_keys_str_mv AT panagopoulosioannis thegrepcommandbutnotfusionmapfusionfinderorchimerascancapturesthecicdux4fusiongenefromwholetranscriptomesequencingdataonasmallroundcelltumorwitht419q35q13
AT gorunovaludmila thegrepcommandbutnotfusionmapfusionfinderorchimerascancapturesthecicdux4fusiongenefromwholetranscriptomesequencingdataonasmallroundcelltumorwitht419q35q13
AT bjerkehagenbodil thegrepcommandbutnotfusionmapfusionfinderorchimerascancapturesthecicdux4fusiongenefromwholetranscriptomesequencingdataonasmallroundcelltumorwitht419q35q13
AT heimsverre thegrepcommandbutnotfusionmapfusionfinderorchimerascancapturesthecicdux4fusiongenefromwholetranscriptomesequencingdataonasmallroundcelltumorwitht419q35q13
AT panagopoulosioannis grepcommandbutnotfusionmapfusionfinderorchimerascancapturesthecicdux4fusiongenefromwholetranscriptomesequencingdataonasmallroundcelltumorwitht419q35q13
AT gorunovaludmila grepcommandbutnotfusionmapfusionfinderorchimerascancapturesthecicdux4fusiongenefromwholetranscriptomesequencingdataonasmallroundcelltumorwitht419q35q13
AT bjerkehagenbodil grepcommandbutnotfusionmapfusionfinderorchimerascancapturesthecicdux4fusiongenefromwholetranscriptomesequencingdataonasmallroundcelltumorwitht419q35q13
AT heimsverre grepcommandbutnotfusionmapfusionfinderorchimerascancapturesthecicdux4fusiongenefromwholetranscriptomesequencingdataonasmallroundcelltumorwitht419q35q13