<div dir="ltr"><div><div>Thank you guys for your suggestions. So you do not suggest to use RNA-seq data from another study, even I assemble them separately and then provide both assemblies into Maker as a comma separated list. The issues you mentioned do exist, but some people did collect RNA-seq data from different individuals and used them for gene annotation (e.g., doi:10.1038/ng.3198). But thank you for your suggestions, I will think about it.<br><br></div>Best<br></div>Quanwei <br></div><div class="gmail_extra"><br><div class="gmail_quote">2017-01-31 16:05 GMT-05:00 Fields, Christopher J <span dir="ltr"><<a href="mailto:cjfields@illinois.edu" target="_blank">cjfields@illinois.edu</a>></span>:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">I agree with Mike.  I also suggest not combining RNA-Seqs from different runs (e.g. different studies) even if they are from the same tissue, development stage etc. There are many other factors (biological variation, sample quality, sequencing chemistry or technology differences, etc) that can significantly and negatively impact trx assembly quality.<br>
<br>
chris<br>
<div class="HOEnZb"><div class="h5"><br>
On 1/31/17, 1:26 PM, "maker-devel on behalf of Michael Campbell" <<a href="mailto:maker-devel-bounces@yandell-lab.org">maker-devel-bounces@yandell-<wbr>lab.org</a> on behalf of <a href="mailto:michael.s.campbell1@gmail.com">michael.s.campbell1@gmail.com</a>> wrote:<br>
<br>
    I would probably try merging the replicates but not the tissues. You can then pass the output files to MAKER in a comma separated list in the opts file.<br>
<br>
    Example:<br>
    est=/PATH/TO/file1.fsata,/<wbr>PATH/TO/file2.fasta<br>
<br>
    Good luck,<br>
    Mike<br>
<br>
    > On Jan 31, 2017, at 2:08 PM, Quanwei Zhang <<a href="mailto:qwzhang0601@gmail.com">qwzhang0601@gmail.com</a>> wrote:<br>
    ><br>
    > Hello:<br>
    ><br>
    > I am trying to assemble transcripts using RNA-seq data by the tool Trinity, which will be used for gene annotation for Maker. Now I have data from two tissues with two replicates each. Should I merge all four samples to get one assembly file? Or should I merge replicates of each tissue separately and use the two assembly files as input of Maker. Merging all samples into one, we will have much higher coverage level, but I think there may be some genes expressed by tissue-specific isoforms. So I not sure whether I should merge RNA-seq from different tissues.<br>
    > What's more, I find some published RNA-seq data from another individual (and also for different tissue from us) for the same species. Should I merge all RNA-seq together (across individuals and tissues)? Or should I generate different transcript assembly and use all those assemblies as input to Maker?<br>
    ><br>
    > Thanks<br>
    > Best<br>
    > Quanwei<br>
    > ______________________________<wbr>_________________<br>
    > maker-devel mailing list<br>
    > <a href="mailto:maker-devel@box290.bluehost.com">maker-devel@box290.bluehost.<wbr>com</a><br>
    > <a href="http://box290.bluehost.com/mailman/listinfo/maker-devel_yandell-lab.org" rel="noreferrer" target="_blank">http://box290.bluehost.com/<wbr>mailman/listinfo/maker-devel_<wbr>yandell-lab.org</a><br>
<br>
<br>
    ______________________________<wbr>_________________<br>
    maker-devel mailing list<br>
    <a href="mailto:maker-devel@box290.bluehost.com">maker-devel@box290.bluehost.<wbr>com</a><br>
    <a href="http://box290.bluehost.com/mailman/listinfo/maker-devel_yandell-lab.org" rel="noreferrer" target="_blank">http://box290.bluehost.com/<wbr>mailman/listinfo/maker-devel_<wbr>yandell-lab.org</a><br>
<br>
<br>
</div></div></blockquote></div><br></div>