Unconfigured Ad

Collapse
X
 
  • Time
  • Show
Clear All
new posts
  • Newsbot!
    Banned
    • Feb 2008
    • 1331

    #1

    PubMed: Integrative analysis of the melanoma transcriptome.

    Syndicated from PubMed RSS Feeds

    Integrative analysis of the melanoma transcriptome.

    Genome Res. 2010 Feb 23;

    Authors: Berger MF, Levin JZ, Vijayendran K, Sivachenko A, Adiconis X, Maguire J, Johnson LA, Robinson J, Verhaak RG, Sougnez C, Onofrio RC, Ziaugra L, Cibulskis K, Laine E, Barretina J, Winckler W, Fisher DE, Getz G, Meyerson M, Jaffe DB, Gabriel SB, Lander ES, Dummer R, Gnirke A, Nusbaum C, Garraway LA

    Global studies of transcript structure and abundance in cancer cells enable the systematic discovery of aberrations that contribute to carcinogenesis, including gene fusions, alternative splice isoforms, and somatic mutations. We developed a systematic approach to characterize the spectrum of cancer-associated mRNA alterations through integration of transcriptomic and structural genomic data, and we applied this approach to generate new insights into melanoma biology. Using paired-end massively parallel sequencing of cDNA (RNA-seq) together with analyses of high-resolution chromosomal copy number data, we identified 11 novel melanoma gene fusions produced by underlying genomic rearrangements, as well as 12 novel readthrough transcripts. We mapped these chimeric transcripts to base-pair resolution and traced them to their genomic origins using matched chromosomal copy number information. We also used these data to discover and validate base-pair mutations that accumulated in these melanomas, revealing a surprisingly high rate of somatic mutation and lending support to the notion that point mutations constitute the major driver of melanoma progression. Taken together, these results may indicate new avenues for target discovery in melanoma, while also providing a template for large-scale transcriptome studies across many tumor types.

    PMID: 20179022 [PubMed - as supplied by publisher]



    More...
  • Jon_Keats
    Senior Member
    • Mar 2010
    • 279

    #2
    Anyone know how (they/you can) mapped to the reference "transcriptome" and then cross-mapped those alignments to the genome to capture exon-exon junctions? I really like this analysis concept. Seems simple enough to align a dataset to a transcriptome reference but then how to you assign these alignments to genome coordinates?

    Comment

    • newbee
      Junior Member
      • Apr 2010
      • 1

      #3
      BAM file flagstats

      I have never used RNA seq data -- so I decided to repeat some analysis from above paper:

      I download supplemental data from

      http://www.ncbi.nlm.nih.gov/sites/en...&term=GSE17593[Accession]&cmd=search



      Run samtools flagstat on bam file :

      I have hard time understanding the flagstats output

      For e.g.
      samtools flagstats GSM506401_MEWO.paired.bwa.remapped.bam > GSM506401_MEWO.paired.bwa.remapped.bam.flagstats

      GSM506401_MEWO.paired.bwa.remapped.bam.flagstats
      10361094 in total
      0 QC failure
      3784092 duplicates
      10361094 mapped (100.00%)
      10361094 paired in sequencing
      5180547 read1
      5180547 read2
      0 properly paired (0.00%)
      10361094 with itself and mate mapped
      0 singletons (0.00%)
      0 with mate mapped to a different chr
      0 with mate mapped to a different chr (mapQ>=5)

      Sample :
      SL-XAN_2_FC30BV1AAXX:1:47:1108:971 161 chrM 2 37 51M = 416 0 ATCACAGGTCTATCACCCTATTAACCACTC
      ACGGGAGCTCTCCATGCATTT 77777777777777777777777777777766644122122222222.222 XT:A:U NM:i:0 X0:i:1 X1:i:0 XM:i:0 XO:i:0 XG:i:0
      MD:Z:51
      SL-XAN_2_FC30BV1AAXX:1:17:1367:1052 97 chrM 3 37 51M = 411 0 TCACAGGTCTATCACCCTATTAACCACTCA
      CGGGAGCTCTCCATGCATTTT 77777777777777777777777777777774+44-222222/2222222" XT:A:U NM:i:1 X0:i:1 X1:i:0 XM:i:1 XO:i:0 XG:i:0
      MD:Z:50G0
      SL-XAN_2_FC30BV1AAXX:1:48:1493:1675 161 chrM 3 37 51M = 396 0 TCACAGGGCTATCACCCTATTAACCACTCA
      CGGGAGCGCTCCATGCATTTG 66,66666666666666666(66666666$66'6+,2'22!)'2,22,222 XT:A:U NM:i:2 X0:i:1 X1:i:0 XM:i:2 XO:i:0 XG:i:0
      MD:Z:7T29T13
      SL-XAN_2_FC30BV1AAXX:1:2:1699:495 1121 chrM 3 37 51M = 411 0 TCACAGGTCTATCACCCTATTAACCACTCA
      CGGGAGCTCTCCATGCATTTT 7753577770757777777777777-777-71+1'/22221/%2/2*222" XT:A:U NM:i:1 X0:i:1 X1:i:0 XM:i:1 XO:i:0 XG:i:0
      MD:Z:50G0
      SL-XAN_2_FC30BV1AAXX:1:94:521:1276 97 chrM 4 37 51M = 326 0 CACAGGTCTATCACCCTATTAACCACTCAC
      GGGAGCTCTCCATGCATTTGG 77777777777777777777777777777707777222222/2/2*222/% XT:A:U NM:i:0 X0:i:1 X1:i:0 XM:i:0 XO:i:0 XG:i:0
      MD:Z:51



      a) I am assuming these are mapped bam files (From paper : "Alignments were performed using the Burrows-Wheeler Alignment Tool (BWA), allowing up to four mismatches with the reference")
      b) Please can somebody help me understand flagstats -- it is 100% mapped but there are "0 properly paired" reads
      c) Please also let me know if you have any suggestions for software package -- which can be used explore this data.

      Thanks
      Last edited by newbee; 04-20-2010, 03:01 PM.

      Comment

      Latest Articles

      Collapse

      • SEQadmin2
        Multiomic Liquid Biopsy: What DNA/RNA Panels Reveal That DNA Alone Can’t
        by SEQadmin2



        Tumors contain complex mixtures of cell populations that evolve in response to treatment and host-specific factors. Tissue biopsies provide only static snapshots of a tumor’s state. But what if we could follow that evolution over time using a minimally invasive technique?

        Liquid biopsies can detect tumor-derived material circulating in the bloodstream, providing a way to study treatment response, detect emerging resistance, and discover new biomarkers. Circulating
        ...
        10-14-2026, 09:55 AM
      • SEQadmin2
        New Genomics Technologies Take Aim at Long-Standing Limits
        by SEQadmin2


        Researchers using sequencing and genomics tools often have to make trade-offs. They can choose between speed or scale, short reads or long-range information, or targeted panels or a view of the whole transcriptome. New technologies that have been released this year are built to address those tough choices.

        We asked six companies the same four questions to learn about their latest products. The new technologies bring a lot to the table, including rethinking sequencing
        ...
        09-28-2026, 10:25 AM

      ad_right_rmr

      Collapse

      News

      Collapse

      Topics Statistics Last Post
      Started by SEQadmin2, 10-15-2026, 10:52 AM
      0 responses
      26 views
      0 reactions
      Last Post SEQadmin2  
      Started by SEQadmin2, 10-14-2026, 05:39 AM
      0 responses
      33 views
      0 reactions
      Last Post SEQadmin2  
      Started by SEQadmin2, 10-13-2026, 09:56 AM
      0 responses
      31 views
      0 reactions
      Last Post SEQadmin2  
      Started by SEQadmin2, 10-09-2026, 11:14 AM
      0 responses
      28 views
      0 reactions
      Last Post SEQadmin2  
      Working...