Question: assembly of mouse transcriptome from BAM file
1
gravatar for regulatinggeneexpression
18 months ago by

hello, We work with a mouse strain different from the strain that was used in generating the standard mouse reference genome (mm9 or mm10). The Sanger institute has done NGS on the strain we work with and has a BAM file available on its website.

How can I use that BAM file to assemble a transcriptome that I can use as a reference for analyzing RNA-Seq data from this particular strain? I only care about protein-coding orf so I do not need to do de novo genome assembly.

Thanks.

rna-seq • 898 views
ADD COMMENTlink modified 17 months ago • written 18 months ago by regulatinggeneexpression10
1

What kind of BAM file is it? If it is for WGS then you can't use it directly to assemble a transcriptome. If it is from RNAseq data then you could use one of the options mentioned below by @grant after extracting the reads from that BAM file.

ADD REPLYlink written 18 months ago by genomax70k

hello mmfansier, Thanks for suggesting mmseq. I found a genome for my mouse strain. I tried opening the genome in IGV. I had to gunzip, and make a .genome file as described in IGV website (https://software.broadinstitute.org/software/igv/LoadGenome). IGV does not read the fasta genome file. Do you have any advice? Thanks.

ADD REPLYlink written 17 months ago by regulatinggeneexpression10

Please use ADD COMMENT/ADD REPLY when responding to existing posts to keep threads logically organized. This comment belongs under @mmfansler's answer.

IGV does not read the fasta genome file.

Are you getting an error? Are there just fasta sequences of transcripts in the file?

ADD REPLYlink modified 17 months ago • written 17 months ago by genomax70k

Just to test, I downloaded a transcriptome from there, unzipped, then indexed (igvtools index transcriptome.fa). This loaded fine in IGV. As mentioned by @genomax, we'll need more details to help further.

ADD REPLYlink written 17 months ago by mmfansler320
0
gravatar for grant.hovhannisyan
18 months ago by
grant.hovhannisyan1.7k wrote:

Depends on what kind of analysis you want to do further with your data. But in general case you have two main options:

  1. De novo transcriptome assembly, using for example Trinity software.
  2. Reference guided transcriptome assembly, using for example StringTie

In general, reference based assembly is more accurate than de novo assembly, and since most probably genomes of mice strains are not very different (this is my guess), I think reference based assembly is more suitable.

ADD COMMENTlink modified 18 months ago • written 18 months ago by grant.hovhannisyan1.7k

Should note, Trinity also has a reference-guided mode. I agree re: reference-guided assemblies, with one caveat: this depends on the quality of the reference genome assembly; I wouldn't recommend this with a poor quality genome assembly, where genes may be fragmented. With mouse this isn't an issue.

ADD REPLYlink written 18 months ago by Chris Fields2.1k

Thank you, I moved my question to a more relevant threat although not recent one

ADD REPLYlink modified 17 months ago • written 17 months ago by F3.4k

Did you post this in the wrong thread by mistake?

ADD REPLYlink modified 17 months ago • written 17 months ago by genomax70k

Sorry, no, I noticed about transcriptome and trinity as my read counts come from RNA-seq and RSEM method I thought to ask here. Sorry if was irrelevant. I seek recent threats but nothing about transcriptome, thus I thought not to create a new post for this question

ADD REPLYlink modified 17 months ago • written 17 months ago by F3.4k
0
gravatar for mmfansler
18 months ago by
mmfansler320
MSKCC | New York, NY
mmfansler320 wrote:

mmseq is a project that has transcriptomes assembled for the Sanger mouse strains as well as tools to do your own assembly if needed (e.g.,mouse_strain_transcriptome.sh).

ADD COMMENTlink modified 18 months ago • written 18 months ago by mmfansler320
Please log in to add an answer.

Help
Access

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 2183 users visited in the last hour