This is a test version of Biostars. For the public version, visit https://www.biostars.org.
[ASK| About CDR3 gene region extraction and clustering

Dear all Biostar community members,

First of all, I apologize for asking about beginner things. This is my first time posting here. I have VH gene FASTQ data of bio-panning products to be analyzed. I usually use USEARCH in LINUX System in order to merge, pipelining and cluster. But, the result is always the full length of the VH gene. Recently, I would extract HCDR3 of the whole population of each bio-panning for better understanding. But, I am really confused about how to extract just the HCDR3 region. Would anyone mind helping me with my problem? I am a beginner in the python language programming.

Thank you

cdr3 antibody
• 1,098 views • link
updated by Jesse • Jesse • written by muhamamadfeisaljatnika • muhamamadfeisaljatnika

1 answer

If you mean you have full variable region sequences (V+D+J segments included) and you just want to extract the CDRH3 region from those, I'd recommend IgBLAST. It has a command-line version you can run with multiple output formats including standard AIRR TSV that makes it easy to parse out the cdr3 column. The command-line version takes a little bit of setup because you have to provide it with the databases of V+D+J germline sequences plus a table of J gene attributes (the "auxiliary data") if you want that junction info, though. There's a guide for the command-line version to get that set up.

• 0 views • link
written by Jesse • Jesse

Log in to answer this question.

More posts like this