SEQanswers

Go Back   SEQanswers > Bioinformatics > Bioinformatics



Similar Threads
Thread Thread Starter Forum Replies Last Post
Zorro maximum size flobpf Bioinformatics 1 11-13-2011 09:39 PM
Splitting fasta AND quality on masked sequence mucku Bioinformatics 3 12-15-2010 12:53 PM
bwa masked sequence? rand Bioinformatics 2 11-09-2009 07:02 AM
why consed finds discrepancies with masked n sequence alig Illumina/Solexa 3 06-12-2009 01:32 PM
Masked/Unmasked Reference Genome ytmnd85 General 5 05-31-2009 04:52 PM

Reply
 
Thread Tools
Old 04-13-2011, 09:13 AM   #1
glacerda
Member
 
Location: Brazil

Join Date: Aug 2008
Posts: 27
Smile Zorro: The Masked Assembler (first public release)

ZORRO is an hybrid sequencing technology assembler. It takes 2 sets of pre-assembled contigs and merge them into a more contiguous and consistent assembly. We have already tested Zorro with Illumina Solexa and 454 from some of organisms varying from 3Mb to 100Mb. The main caracteristic of Zorro is the treatment before and after assembly to avoid errors.

The ZORRO project is maintained by Gustavo Lacerda, Ramon Vidal and Marcelo Carazzole and were first used in this Yeast assembly: Genome structure of a Saccharomyces cerevisiae strain widely used in bioethanol production

ZORRO needs to be better documented and has not undergone enough testing. If you want to discuss the pipeline you can join the mailing list: zorro-google group

ZORRO PIPELINE
Zorro is based on the minimus2 pipeline (AMOS package) and uses MuMMer,
AMOS and bowtie in its internals. Zorro takes 2 contigs fasta files as
input (representing assembled contigs from a whole genome assembly)
and one fasta file containing some of the reads used for assembly
(only 10X coverage is enough, more will slow down the pipeline and
consume more resources).

Zorro initial phase detect inconsistencies in the assemblies and split
the contigs where they occur. Next, zorro counts k-mers (default k=22)
in the reads and use the k-mer count table to detect and mask repeats
in both assembly1 and assembly2. After repeat masking, zorro uses nucmer
to detect overlaps between assembly1 and assembly2 (no overlaps between
contigs from the same assembly are allowed). All overlaps found in this
phase are expected to be between unique regions (because repeats are
masked). The overlaps are used to layout and generate consensus for the
merged contigs, using AMOS tools. The merged contigs are built using the
unmasked contigs, so the final merged assembly should include the repeat
regions.

Another round of assembly, less stringent, tries to merge contigs that
were not included in the first Zorro phase. All the contigs are outputted
to <prefix>.ZORRO.fasta. We recommend the use of SSPACE to scaffold the
ZORRO contigs.


Zorro Website: www.lge.ibi.unicamp.br/zorro
glacerda is offline   Reply With Quote
Old 04-16-2011, 04:01 AM   #2
RLB_84
Junior Member
 
Location: Washington, DC

Join Date: Nov 2010
Posts: 4
Default

hi, I read about your pipeline and it seems quite interesting, as I'm a newbie user of MUMmer. So I hope I could try your assembler as soon as possible. By the way, just one question about "masking": when I consider a contig dataset returned by Newbler, I should keep in mind that many of them can be repeated and I can't figure it from the contig fasta file (I need the other Newbler outputs). So, how does your assembler deal with the probable occurrence of a contig in more than one genome locus?
RLB_84 is offline   Reply With Quote
Old 04-16-2011, 05:00 AM   #3
glacerda
Member
 
Location: Brazil

Join Date: Aug 2008
Posts: 27
Default

Hi RLB_84, that's a very good question.

In addition to the contig fasta files, Zorro takes as input the reads file (a subsample of WGS reads). The reads are used only to allow us to identify repeats in the contigs. I will explain technically

1-Zorro counts the occurences of 22-mers in the reads file supplied by the user
2-22-mer words that are unique in the genome should occur proportionally to the genome coverage
3-22-mer words that represent repeats should occur at least twice the peak coverage
4-we select the 22-mer words that occur at least twice the mode of the distribution. These 22-mer words are used to mask the contig files using bowtie.

This technique is used by many ab initio repeat detection software. We do not need to screen repeat libraries and, even if newbler (or other software) has collapsed the repeats, we coulod still detect them.
glacerda is offline   Reply With Quote
Old 04-16-2011, 06:11 AM   #4
RLB_84
Junior Member
 
Location: Washington, DC

Join Date: Nov 2010
Posts: 4
Default

Thanks for the clarification, that sounds good! So, considering I'm trying to merge two sets of contigs returned by Newbler and AbySS, this approach can make feasible to compare directly the coverage of the datasets, as Newbler comes out with a bunch of files (very useful about coverage and so on), while AbySS output needs more processing.

Does the "less stringent" assembly phase take into account the repeat prediction and try to infer this information in the assembly itself?
RLB_84 is offline   Reply With Quote
Reply

Tags
assembly, repeats, zorro

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 03:07 AM.


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