<html>
<head>
<meta http-equiv="Content-Type" content="text/html; charset=Windows-1252">
<style type="text/css" style="display:none"><!--P{margin-top:0;margin-bottom:0;} --></style>
</head>
<body dir="ltr" style="font-size:12pt;color:#000000;background-color:#FFFFFF;font-family:Calibri,Arial,Helvetica,sans-serif;">
So I have assembled my transcriptome with Trinity using the jaccard clip option and I have run maker with and without corrected_est_fusion.
<p><br>
</p>
<p>I have then use SNAP to train/filter it with:</p>
<p><br>
</p>
<p>maker2zff specie.all.gff</p>
<p><br>
</p>
<p>Here are my results: </p>
<p><br>
</p>
<p>Number of gene after maker -> Number of gene after maker2zff</p>
<p><br>
</p>
<p>- Without corrected_est_fusion: 21621 -> 13875</p>
<p>- With corrected_est_fusion: 16850 -> 9098</p>
<p><br>
</p>
<p>1 )If I understand well how works corrected_est_fusion, because it prevents gene merging, shouldn't be the invert ?
</p>
<p>Normally I should find more genes with corrected_est_fusion right ?</p>
<p><br>
</p>
<p>2) I think I should find something like 13000-14000 genes for my specie. SHould I go with the "Without corrected_est_fusion" for the 2nd iteration of maker ?</p>
<p><br>
</p>
<p> Thanks for your help <br>
</p>
<p><br>
</p>
<p><br>
</p>
<div id="Signature">
<div name="divtagdefaultwrapper" style="font-family:Calibri,Arial,Helvetica,sans-serif; font-size:; margin:0">
<br>
Patrick Tran Van<br>
<br>
Groups Chapuisat, Robinson-Rechavi & Schwander<br>
Department of Ecology and Evolution<br>
University of Lausanne<br>
Le Biophore<br>
CH-1015 Lausanne<br>
Switzerland<br>
Office 3206<br>
<br>
</div>
</div>
<div style="word-wrap:break-word">
<hr tabindex="-1" style="display:inline-block; width:98%">
<div id="divRplyFwdMsg" dir="ltr"><font style="font-size:11pt" face="Calibri, sans-serif" color="#000000"><b>From:</b> Carson Holt <carsonhh@gmail.com><br>
<b>Sent:</b> Monday, June 26, 2017 11:38 PM<br>
<b>To:</b> Patrick Tran Van<br>
<b>Cc:</b> maker-devel@yandell-lab.org<br>
<b>Subject:</b> Re: [maker-devel] Advice on my pipeline</font>
<div> </div>
</div>
<div>Sorry the option is —> correct_est_fusion
<div class=""><br class="">
</div>
<div class="">It is in the maker_opts.ctl file.</div>
<div class=""><br class="">
</div>
<div class="">I would use both SNAP and Augustus on a few large contigs then review the results manually. If one of them is not behaving well, then drop it. If both behave well (i.e. correlate well with evidence alignemnts) then keep them both.</div>
<div class=""><br class="">
</div>
<div class="">—Carson</div>
<div class=""><br class="">
</div>
<div class=""><br class="">
</div>
<div class=""><br class="">
<div>
<blockquote type="cite" class="">
<div class="">On Jun 26, 2017, at 3:48 AM, Patrick Tran Van <<a href="mailto:Patrick.TranVan@unil.ch" class="">Patrick.TranVan@unil.ch</a>> wrote:</div>
<br class="Apple-interchange-newline">
<div class="">
<div class="" style="margin-top:0px; margin-bottom:0px; font-family:Calibri,Arial,Helvetica,sans-serif; font-size:16px; font-style:normal; font-weight:normal; letter-spacing:normal; text-align:start; text-indent:0px; text-transform:none; white-space:normal; word-spacing:0px; background-color:rgb(255,255,255)">
Thanks for your answer.</div>
<div class="" style="margin-top:0px; margin-bottom:0px; font-family:Calibri,Arial,Helvetica,sans-serif; font-size:16px; font-style:normal; font-weight:normal; letter-spacing:normal; text-align:start; text-indent:0px; text-transform:none; white-space:normal; word-spacing:0px; background-color:rgb(255,255,255)">
<br class="">
</div>
<div class="" style="margin-top:0px; margin-bottom:0px; font-family:Calibri,Arial,Helvetica,sans-serif; font-size:16px; font-style:normal; font-weight:normal; letter-spacing:normal; text-align:start; text-indent:0px; text-transform:none; white-space:normal; word-spacing:0px; background-color:rgb(255,255,255)">
1) Do you think that adding a Augustus training in addition to SNAP at the step 3 and 5 will add more confidence (instead of adding Augustus only for the final round) ?</div>
<div class="" style="margin-top:0px; margin-bottom:0px; font-family:Calibri,Arial,Helvetica,sans-serif; font-size:16px; font-style:normal; font-weight:normal; letter-spacing:normal; text-align:start; text-indent:0px; text-transform:none; white-space:normal; word-spacing:0px; background-color:rgb(255,255,255)">
Because I am using autoAug for this and it tooks a while to compute ..</div>
<div class="" style="margin-top:0px; margin-bottom:0px; font-family:Calibri,Arial,Helvetica,sans-serif; font-size:16px; font-style:normal; font-weight:normal; letter-spacing:normal; text-align:start; text-indent:0px; text-transform:none; white-space:normal; word-spacing:0px; background-color:rgb(255,255,255)">
<br class="">
</div>
<div class="" style="margin-top:0px; margin-bottom:0px; font-family:Calibri,Arial,Helvetica,sans-serif; font-size:16px; font-style:normal; font-weight:normal; letter-spacing:normal; text-align:start; text-indent:0px; text-transform:none; white-space:normal; word-spacing:0px; background-color:rgb(255,255,255)">
2) I don't see this option : 'avoid_est_fusion=1' . I have tried to add it but I got this error:</div>
<div class="" style="margin-top:0px; margin-bottom:0px; font-family:Calibri,Arial,Helvetica,sans-serif; font-size:16px; font-style:normal; font-weight:normal; letter-spacing:normal; text-align:start; text-indent:0px; text-transform:none; white-space:normal; word-spacing:0px; background-color:rgb(255,255,255)">
<br class="">
</div>
<div class="" style="margin-top:0px; margin-bottom:0px; font-family:Calibri,Arial,Helvetica,sans-serif; font-size:16px; font-style:normal; font-weight:normal; letter-spacing:normal; text-align:start; text-indent:0px; text-transform:none; white-space:normal; word-spacing:0px; background-color:rgb(255,255,255)">
WARNING: Invalid option 'avoid_est_fusion' in control file maker_opts.ctl</div>
<div class="" style="margin-top:0px; margin-bottom:0px; font-family:Calibri,Arial,Helvetica,sans-serif; font-size:16px; font-style:normal; font-weight:normal; letter-spacing:normal; text-align:start; text-indent:0px; text-transform:none; white-space:normal; word-spacing:0px; background-color:rgb(255,255,255)">
<br class="">
</div>
<div class="" style="margin-top:0px; margin-bottom:0px; font-family:Calibri,Arial,Helvetica,sans-serif; font-size:16px; font-style:normal; font-weight:normal; letter-spacing:normal; text-align:start; text-indent:0px; text-transform:none; white-space:normal; word-spacing:0px; background-color:rgb(255,255,255)">
(I am using v 2.31.8 )<br class="">
</div>
<div class="" style="margin-top:0px; margin-bottom:0px; font-family:Calibri,Arial,Helvetica,sans-serif; font-size:16px; font-style:normal; font-weight:normal; letter-spacing:normal; text-align:start; text-indent:0px; text-transform:none; white-space:normal; word-spacing:0px; background-color:rgb(255,255,255)">
<br class="">
</div>
<div id="Signature" class="" style="font-family:Calibri,Arial,Helvetica,sans-serif; font-size:16px; font-style:normal; font-weight:normal; letter-spacing:normal; text-align:start; text-indent:0px; text-transform:none; white-space:normal; word-spacing:0px; background-color:rgb(255,255,255)">
<div name="divtagdefaultwrapper" class="" style="font-family:Calibri,Arial,Helvetica,sans-serif; margin:0px">
<br class="">
Patrick Tran Van<br class="">
<br class="">
Groups Chapuisat, Robinson-Rechavi & Schwander<br class="">
Department of Ecology and Evolution<br class="">
University of Lausanne<br class="">
Le Biophore<br class="">
CH-1015 Lausanne<br class="">
Switzerland<br class="">
Office 3206<br class="">
<br class="">
</div>
</div>
<div class="" style="font-family:Calibri,Arial,Helvetica,sans-serif; font-size:16px; font-style:normal; font-weight:normal; letter-spacing:normal; text-align:start; text-indent:0px; text-transform:none; white-space:normal; word-spacing:0px; background-color:rgb(255,255,255); word-wrap:break-word">
<hr tabindex="-1" class="" style="display:inline-block; width:1067.21875px">
<div id="divRplyFwdMsg" dir="ltr" class=""><font class="" style="font-size:11pt" face="Calibri, sans-serif"><b class="">From:</b><span class="Apple-converted-space"> </span>Carson Holt <<a href="mailto:carsonhh@gmail.com" class="">carsonhh@gmail.com</a>><br class="">
<b class="">Sent:</b><span class="Apple-converted-space"> </span>Monday, June 5, 2017 8:29 PM<br class="">
<b class="">To:</b><span class="Apple-converted-space"> </span>Patrick Tran Van<br class="">
<b class="">Cc:</b><span class="Apple-converted-space"> </span><a href="mailto:maker-devel@yandell-lab.org" class="">maker-devel@yandell-lab.org</a><br class="">
<b class="">Subject:</b><span class="Apple-converted-space"> </span>Re: [maker-devel] Advice on my pipeline</font>
<div class=""> </div>
</div>
<div class="">Your plan sounds good. A couple of related notes.
<div class=""><br class="">
</div>
<div class="">Insect genomes tend to have high gene density, so gene merging will be the primary difficulty. You can avoid merging of mRNA-seq evidence by using options like jaccard_clip in Trinity. Then use avoid_est_fusion=1 inside of MAKER.</div>
<div class=""><br class="">
</div>
<div class="">Also it is more convenient to do each run in the same directory rather than supplying the previous run as GFF3 input. MAKER will automatically recycle previous results archived in the run directory when you do this. Using the maker_gff option
is really more for getting data into the run from jobs performed a long time ago (so they can’t be run in the same directory).</div>
<div class=""><br class="">
</div>
<div class="">—Carson</div>
<div class=""><br class="">
</div>
<div class=""><br class="">
<div class="">
<blockquote type="cite" class="">
<div class="">On Jun 2, 2017, at 3:56 AM, Patrick Tran Van <<a href="mailto:Patrick.TranVan@unil.ch" class="">Patrick.TranVan@unil.ch</a>> wrote:</div>
<br class="Apple-interchange-newline">
<div class="">
<div class="" style="margin-top:0px; margin-bottom:0px; font-family:Calibri,Arial,Helvetica,sans-serif; font-size:16px; font-style:normal; font-weight:normal; letter-spacing:normal; text-align:start; text-indent:0px; text-transform:none; white-space:normal; word-spacing:0px; background-color:rgb(255,255,255)">
Hello,<br class="">
<br class="">
This is my first time running Maker for an insect genome annotation.<span class="Apple-converted-space"> </span><br class="">
<br class="">
I have found various resources and tried to make a consensus, I am looking for your thoughts and advices about my pipeline, if I can improve something or doing useless things:<br class="">
<br class="">
<br class="">
What I have:<br class="">
- RNA evidence: transcriptome<br class="">
- Proteine evidence: swissprot/uniprot + busco protein set of insect<br class="">
- Cegma and busco results of my genome<br class="">
<br class="">
<br class="">
1) Train SNAP with CEGMA<br class="">
<br class="">
2) Run (run A) maker with repeat masking with transcript, protein, the new SNAP file (from step 1) and augustus file (from busco).<br class="">
<br class="">
3) Create SNAP model from run A.<br class="">
<br class="">
4) Run (run B ) with the new SNAP (done at step 3) with options turned off (est2genome=0) and (protein2genome=0) data, provide gff file (maker_gff=run_A.gff), turn off repeat masking (rm_pass=1), and use previous mapping results (altest_pass=1 and protein_pass=1).<br class="">
<br class="">
5) Create SNAP model from run B.<br class="">
<br class="">
6) Run (run C) with the new SNAP (done at step 5) with options turned off (est2genome=0) and (protein2genome=0) data, provide gff file (maker_gff=run_B.gff), turn off repeat masking (rm_pass=1), and use previous mapping results (altest_pass=1 and protein_pass=1).<br class="">
<br class="">
7) Create SNAP model from run C AND Create Augustus gene model from run C<br class="">
<br class="">
8) Run (run D) with the new SNAP (done at step 7) + AUGUSTUS file (step 7) with options turned off (est2genome=0) and (protein2genome=0) data, provide gff file (maker_gff=run_C.gff), turn off repeat masking (rm_pass=1), and use previous mapping results (altest_pass=1
and protein_pass=1). + Use keep_preds=1<br class="">
<br class="">
<br class="">
<br class="">
Does it seems coherent ?</div>
<div class="" style="margin-top:0px; margin-bottom:0px; font-family:Calibri,Arial,Helvetica,sans-serif; font-size:16px; font-style:normal; font-weight:normal; letter-spacing:normal; text-align:start; text-indent:0px; text-transform:none; white-space:normal; word-spacing:0px; background-color:rgb(255,255,255)">
<br class="">
</div>
<div class="" style="margin-top:0px; margin-bottom:0px; font-family:Calibri,Arial,Helvetica,sans-serif; font-size:16px; font-style:normal; font-weight:normal; letter-spacing:normal; text-align:start; text-indent:0px; text-transform:none; white-space:normal; word-spacing:0px; background-color:rgb(255,255,255)">
Cheers,<br class="">
</div>
<div id="Signature" class="" style="font-family:Calibri,Arial,Helvetica,sans-serif; font-size:16px; font-style:normal; font-weight:normal; letter-spacing:normal; text-align:start; text-indent:0px; text-transform:none; white-space:normal; word-spacing:0px; background-color:rgb(255,255,255)">
<div name="divtagdefaultwrapper" class="" style="font-family:Calibri,Arial,Helvetica,sans-serif; margin:0px">
<br class="">
Patrick Tran Van<br class="">
<br class="">
Groups Chapuisat, Robinson-Rechavi & Schwander<br class="">
Department of Ecology and Evolution<br class="">
University of Lausanne<br class="">
Le Biophore<br class="">
CH-1015 Lausanne<br class="">
Switzerland<br class="">
Office 3206<br class="">
<br class="">
</div>
</div>
<span class="" style="font-family:Calibri,Arial,Helvetica,sans-serif; font-size:16px; font-style:normal; font-weight:normal; letter-spacing:normal; text-align:start; text-indent:0px; text-transform:none; white-space:normal; word-spacing:0px; background-color:rgb(255,255,255); float:none; display:inline!important">_______________________________________________</span><br class="" style="font-family:Calibri,Arial,Helvetica,sans-serif; font-size:16px; font-style:normal; font-weight:normal; letter-spacing:normal; text-align:start; text-indent:0px; text-transform:none; white-space:normal; word-spacing:0px; background-color:rgb(255,255,255)">
<span class="" style="font-family:Calibri,Arial,Helvetica,sans-serif; font-size:16px; font-style:normal; font-weight:normal; letter-spacing:normal; text-align:start; text-indent:0px; text-transform:none; white-space:normal; word-spacing:0px; background-color:rgb(255,255,255); float:none; display:inline!important">maker-devel
mailing list</span><br class="" style="font-family:Calibri,Arial,Helvetica,sans-serif; font-size:16px; font-style:normal; font-weight:normal; letter-spacing:normal; text-align:start; text-indent:0px; text-transform:none; white-space:normal; word-spacing:0px; background-color:rgb(255,255,255)">
<a href="mailto:maker-devel@box290.bluehost.com" class="" style="font-family:Calibri,Arial,Helvetica,sans-serif; font-size:16px; font-style:normal; font-weight:normal; letter-spacing:normal; text-align:start; text-indent:0px; text-transform:none; white-space:normal; word-spacing:0px; background-color:rgb(255,255,255)">maker-devel@box290.bluehost.com</a><br class="" style="font-family:Calibri,Arial,Helvetica,sans-serif; font-size:16px; font-style:normal; font-weight:normal; letter-spacing:normal; text-align:start; text-indent:0px; text-transform:none; white-space:normal; word-spacing:0px; background-color:rgb(255,255,255)">
<a href="http://box290.bluehost.com/mailman/listinfo/maker-devel_yandell-lab.org" class="" style="font-family:Calibri,Arial,Helvetica,sans-serif; font-size:16px; font-style:normal; font-weight:normal; letter-spacing:normal; text-align:start; text-indent:0px; text-transform:none; white-space:normal; word-spacing:0px; background-color:rgb(255,255,255)">http://box290.bluehost.com/mailman/listinfo/maker-devel_yandell-lab.org</a><br class="" style="font-family:Calibri,Arial,Helvetica,sans-serif; font-size:16px; font-style:normal; font-weight:normal; letter-spacing:normal; text-align:start; text-indent:0px; text-transform:none; white-space:normal; word-spacing:0px; background-color:rgb(255,255,255)">
</div>
</blockquote>
</div>
<br class="">
</div>
</div>
</div>
<br class="Apple-interchange-newline">
</div>
</blockquote>
</div>
<br class="">
</div>
</div>
</div>
</body>
</html>