Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Significantly different versions of STAR in star_rsem (2.7.6a) and star (2.6.1d) #498

Closed
chewgl opened this issue Nov 15, 2020 · 3 comments

Comments

@chewgl
Copy link

chewgl commented Nov 15, 2020

Great work on the latest version, the increased modularity is much appreciated.

I did notice that the new star_rsem aligner option uses a newer version of STAR (2.7.6a) that differs from the default (2.6.1d). 2.6.1d was previous specified to maintain iGenomes compatibility, but 2.7.6a would break that. In addition, the same indices cannot be used between these versions.

Are there long term plans for resolving this? Perhaps a note of this could be made in the documentation.

Best regards,
Guo-Liang

@apeltzer
Copy link
Member

Hi Guo-Liang,

yes, there are either plans to update the STAR indices to provide compatible version for STAR v2.7+ or switch to some other solution. Currently very unclear on how we will proceed here, but there are ideas and discussion around this already.

@drpatelh
Copy link
Member

drpatelh commented Nov 16, 2020

Hi @chewgl ! Well spotted! Yes, as @apeltzer mentioned this is a tricky situation because the AWS iGenomes files require an older version of STAR as specified in the module file:

conda (params.enable_conda ? "bioconda::star=2.6.1d" : null)

which is why I chose for the standard --aligner star_salmon route through the pipeline to use that version. Mainly so things don't break for backward compatibility.

However, for the --aligner star_rsem route given that the indices are not on AWS iGenomes and will need to be re-generated anyway with all of the RSEM specific files, I built a multi-tool Biocontainer with the latest version of STAR as specified here

conda (params.enable_conda ? "bioconda::rsem=1.3.3 bioconda::star=2.7.6a" : null)

As you mentioned, this is slightly annoying but it was a toss up between using a much older version of STAR or getting everyone to re-build the indices if they are using --aligner star_rsem.

@drpatelh
Copy link
Member

drpatelh commented Dec 9, 2020

Added a note in the main README summarising this discrepancy @chewgl 58b702f

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Projects
None yet
Development

No branches or pull requests

3 participants