Question: Rearrangement of genome - can I try a local denovo assembly
gravatar for philippe.mcqt
4.4 years ago by
philippe.mcqt30 wrote:


I have a bunch of WGS sequences for a variety of cultivars at a decent coverage (20X). I am interested in a specific region of the genome that I suspect as undergo rearrangement in resistant cultivars. I have alignment data to a reference that does not have the rearrangement. I would like to reconstruct that region and compare it to non resistant cultivars. How would you do that ?

I can perform a de novo assembly on the whole datasets for every cultivars and look for the region of interest (via blast I guess) to see if any contigs match the region. This seems like a bit overkill.

I would like to select only the reads of the region of interest but not sure how to since the reference does not have this rearrangement. Have the reads describing it been thrown away during the mapping, not sure... Could I use the region that is flanking and conserved as a core and then extend this region by capturing overlapping reads ? Not sure how to do it, any suggestions ?

sequencing assembly genome • 982 views
ADD COMMENTlink modified 4.4 years ago by bdelolmo10 • written 4.4 years ago by philippe.mcqt30
gravatar for bdelolmo
4.4 years ago by
bdelolmo10 wrote:

If the reference genome is well stablished I would suggest you to try some SV algorithms that are available. Paired End Mapping (PEM) algorithms need explicitly paired-end reads. If you have single-end you can try read depth (RD) or split-read methods (SR).

There is a wide variety, I give you some examples: - Breakdancer (PEM). - Pindel (SR). - Delly (PEM+SR). - CREST (assembly of split reads). - Lumpy (PEM + SR and can integrate RD methods). - CNVnator for analyzing RD.

It is better to combine some of them, or to use algorithms that combine different methodologies, such as PEM and SR.

ADD COMMENTlink written 4.4 years ago by bdelolmo10

I have PEM reads. Thanks for suggestion I am looking into it.

ADD REPLYlink written 4.4 years ago by philippe.mcqt30
Please log in to add an answer.


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