SEQanswers

Go Back   SEQanswers > Bioinformatics > Bioinformatics



Similar Threads
Thread Thread Starter Forum Replies Last Post
Why repeat masking is important prior to gene model construction during de novo genom zxyeo General 0 02-26-2014 10:39 PM
Repeat masker - how to use it for estimating repeat expression levels? Alex234 RNA Sequencing 0 11-28-2013 04:44 AM
pre-filtering before running alignment: help needed angerusso Bioinformatics 2 11-15-2011 07:49 PM
Repeat masking Amy S Bioinformatics 2 06-02-2011 04:32 PM
Trimming SOLiD data pre-alignment Michael.James.Clark SOLiD 0 03-04-2011 01:47 PM

Reply
 
Thread Tools
Old 08-20-2014, 01:19 AM   #1
WhatsOEver
Senior Member
 
Location: Germany

Join Date: Apr 2012
Posts: 215
Default Pre-alignment repeat masking

Hi @all,

I'm working currently on human whole genome NGS data and have some questions about repeat masking:

1) Would you say it's advisable to do masking before mapping to the genome? Why - Why not?

2) What tools could you suggest to run masking on reads? I am familiar with - and tried already - repeatmasker/repeatmodeler but it's working damn slow on reads. I thought of doing a mapping against the repbase?! Alternatively, I could use a repeat masked genome (where repeats are removed, not only in lower case, of course) and hope that repeats are not mapped anymore.

3) Through different filtering criteria, I have to discard already ~8% of my reads before doing any alignments against the reference. I think including a repeat-filtering will increase this drastically. Is there any possibility that this will have an effect on downstream analysis? Do you include such things in your calculation when estimating required read numbers to get a certain coverage?

4) (Not directly related to the question) I have recently tried a mapping against the "new" hg38 and saw that I just got more multi-reads. Having a brief look on the data, these were mainly repetitive sequences... Has anyone observed similar things?

Thanks for reading and any suggestions
WhatsOEver is offline   Reply With Quote
Old 08-20-2014, 01:49 AM   #2
Cytosine
Member
 
Location: Belgium

Join Date: Mar 2014
Posts: 22
Default

Hi,

1) I guess it depends on what you want to achieve. For resequencing runs I would use repeat masking on the genome, for CNV analyses I would not mask the genome.

2) Why mask reads? That doesn't really make sense... Mask the genome you are trying to map to. Hard-masking will prevent the reads from mapping to those regions.

3) Masking repeat regions will probably lead to more unmapped reads, but that should not affect any downstream analyses with great effect. I don't remember ever doing read number correction because of repetitive regions, I guess you could do it.

4) Did not have a look at the hg38 yet.
Cytosine is offline   Reply With Quote
Old 08-20-2014, 01:53 AM   #3
dpryan
Devon Ryan
 
Location: Freiburg, Germany

Join Date: Jul 2011
Posts: 3,480
Default

1) Are we talking about masking the genome or the reads? You can just download the repeatmasked genome, though at least if it's only soft-masked most aligners will just ignore that. In general, I'd say it's not advisable to ignore repeat regions without good reason (it's a pretty big chunk of the genome, afterall).

2) Don't mask the reads, there's no point to that. If you want to filter out reads with long homopolymer runs or that are pure di/tri-nucleotide repeats then that'll save a little time with mapping, but you'd just be spending it with filtering.

3) Certainly if you filter out reads that would otherwise align then you're biasing against some reasons. Aligning doesn't take that long, so just attempt to align them and then filter out reads with crappy MAPQs.

4) Depending on what you included in the reference that might not be very surprising. I know the new release contains a LOT of haplotypes and also centromeric sequence, all of which will increase the mapping rate and the numbers of multimappers. That doesn't surprise me really.
dpryan is offline   Reply With Quote
Old 08-20-2014, 10:32 AM   #4
Brian Bushnell
Super Moderator
 
Location: Walnut Creek, CA

Join Date: Jan 2014
Posts: 2,707
Default

As Devon said, don't map to a masked genome unless you have a specific reason to do so. Often genomes are masked to increase mapping speed, as things that map to repeats are often uninformative, yet can take N times longer to map (where N is the number of repeats) compared to unique sequences; therefore, masking repeats can in some cases make things orders-of-magnitude faster without losing information. But that's mainly an issue with very large conglomerations of assorted sequences. When I was working with human NGS data, neither I nor anyone that I knew used a masked human reference for anything. That would make a few of the reads that should have mapped to repeat areas instead map to unique areas, leading to false variation calls.

Last edited by Brian Bushnell; 08-20-2014 at 10:36 AM.
Brian Bushnell is offline   Reply With Quote
Old 08-20-2014, 10:42 AM   #5
swbarnes2
Senior Member
 
Location: San Diego

Join Date: May 2008
Posts: 912
Default

Don't mask repeated regions before mapping. If a read is repetitive, you want it to map to where it belongs. You don't want it forced to map to the wrong place, because you masked away the right place.

If it really bothers you, filter away repetitive regions from your .bam after mapping.
swbarnes2 is offline   Reply With Quote
Old 08-21-2014, 04:44 AM   #6
WhatsOEver
Senior Member
 
Location: Germany

Join Date: Apr 2012
Posts: 215
Default

Thanks for the answers!
I'm convinced it was a bad idea
I'll do repeat-filtering after mapping then.
WhatsOEver is offline   Reply With Quote
Reply

Thread Tools

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off




All times are GMT -8. The time now is 12:14 AM.


Powered by vBulletin® Version 3.8.9
Copyright ©2000 - 2020, vBulletin Solutions, Inc.
Single Sign On provided by vBSSO