Question: RNA-STAR Failing using hg38 but works for hg19
0
gravatar for aklie
3 months ago by
aklie0
aklie0 wrote:

Hi,

I am trying to align two collections of 4 single-end RNA-Seq fastq files (8 fastq files total) from human skin samples to a reference genome using RNA STAR. I have no problem aligning to the hg19 reference genome but run into memory issues when aligning to the hg38. I've tried twice for each reference genome and have been successful twice with hg19 and failed twice with hg38. The error message I receive when using hg38 as the reference is:

"This job was terminated because it used more memory than it was allocated."

I'm hoping someone can help me understand what's going on here and what I might do to fix it!

Thanks!

star memory usage rna-seq hisat2 • 162 views
ADD COMMENTlink modified 3 months ago by Jennifer Hillman Jackson25k • written 3 months ago by aklie0
2
gravatar for Jennifer Hillman Jackson
3 months ago by
United States
Jennifer Hillman Jackson25k wrote:

Hello,

The job is probably right at the edge of the memory limits available on our clusters. I would suggest trying a few more test runs versus hg38. Nodes can process more that one job at a time, and by chance sometimes two large jobs end up on the same node, triggering this type of error.

You might also consider doing some more QA/QC on your input data, it might help. Reviewing tool options and testing out alternatives that are not as compute-intensive can also lead to successful jobs.

If that doesn't work out, it is good to know that RNA STAR uses a lot of memory. The target genome size is one important factor, although I know that many users map against hg38 successfully with STAR, so it is not the only factor.

When there are failures with STAR that cannot be resolved with input data/parameter changes, HISAT2 is the recommended alternative.

The settings for HISAT2 when matter and need to be adjusted for the downstream tools you plan to use. Please see the Galaxy RNA-seq tutorials here for example usage, including QA/QC: https://galaxyproject.org/learn/

Support FAQs: https://galaxyproject.org/support

Thanks! Jen, Galaxy team

ADD COMMENTlink modified 3 months ago • written 3 months ago by Jennifer Hillman Jackson25k

Thank you for such a detailed answer! That makes some sense given STAR memory requirements and the fact that the hg38 is a larger reference for STAR to build a suffix array. I wanted to run through pipelines with multiple different aligners anyway so I'll see if I can get it working and if not I just won't include it. Thanks for your help!

ADD REPLYlink written 3 months ago by aklie0
Please log in to add an answer.

Help
Access

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