search

ncborcherding / scRepertoire

A toolkit for single-cell immune profiling
https://www.borch.dev/uploads/screpertoire/
MIT License
301 stars 49 forks source link

New cloneCall option to more closely match `immunarch` options? #396

Open jeremymsimon opened 1 month ago

jeremymsimon commented 1 month ago

Hi @ncborcherding I'm wondering if it's at all possible to include another option in cloneCall to be CDR3 AA + V + J, which would more closely match immunarch's TrackClonotypes(aa+v+j) function.

strict in your options includes the C gene, which is one notch more stringent, and aa of course only includes the CDR3 amino acid sequence alone.

Since we're using both it would be amazing if we could be sure all groupings are done consistently

Thanks!

ncborcherding commented 1 month ago

Hey @jeremymsimon

I think you make a great point - that ideally there would be a matching option between the 2 workflows - immunarch and scRepertoire. I think the fundamental problem with making the change would be interupting/changing clonal definitions for users, which makes me hesitate.

I think a work around would be to modify the CTstrict column using something like the stringr package. As along as the header is the same, there should be no break in any of the downstream functions.

I will think more on this though - I have been really interested in making scRepertoire as compatible with other workflows.

Thanks, Nick

jeremymsimon commented 1 month ago

That would be great! Keep me posted, happy to help with testing as well

jeremymsimon commented 1 week ago

Hi @ncborcherding any updates on this?