This is a test version of Biostars. For the public version, visit https://www.biostars.org.
Matching complementary genotypes from 2 different files

I have 2 files tab separated as follows:

file1:

rs8192678 GG VO2max PPARGC1A Higher
rs2282679 AC Vitamin D deficiency

file2:

    rs8192678 GG 
    rs2282679 AC 
    rs8192678 CC
    rs2282679 AG

result: where even its complementary should also match with file1

rs8192678 GG GG VO2max PPARGC1A Higher 
rs8192678 GG CC VO2max PPARGC1A Higher
rs2282679 AC AC Vitamin D deficiency
rs2282679 AC AG Vitamin D deficiency

Thank you

genomics
• 2,140 views • link
updated by Kevin Blighe Kevin Blighe • written by mittu1602 mittu1602

1 answer

This will match just on the rs ID (the piped sort command at the end is not necessary):

awk 'BEGIN {FS="\t"} FNR==NR {key=$1; arrayVars[key]=$0; next} {key=$1; if (arrayVars[key]) split(arrayVars[key], var, "\t"); if (arrayVars[key]) print var[1]"\t"var[2]"\t"$2"\t"var[3]}' file1 file2 | sort -r -k1
rs8192678   GG  GG  VO2max PPARGC1A Higher
rs8192678   GG  CC  VO2max PPARGC1A Higher
rs2282679   AC  AG  Vitamin D deficiency
rs2282679   AC  AC  Vitamin D deficiency

-----------------------------------

This will match on rs ID and genotype, should you require that in the future:

awk 'BEGIN {FS="\t"} FNR==NR {key=$1$2; arrayVars[key]=$0; next} {key=$1$2; if (arrayVars[key]) split(arrayVars[key], var, "\t"); if (arrayVars[key]) print var[1]"\t"var[2]"\t"$2"\t"var[3]}' file1 file2 | sort -r -k1
rs8192678   GG  GG  VO2max PPARGC1A Higher
rs2282679   AC  AC  Vitamin D deficiency

If you want an explanation, then let me know. There are undoubtedly many other ways to do this.

Kevin

• 0 views • link
written by Kevin Blighe Kevin Blighe

Hi Kevin, Thanks for the reply. But my aim is to look for complementary genotypes and there are high chances that file2 would also contain genotypes that are not complementary to file1 genotypes. if something can be looked in that way. but none the less, thank you so much.

• 0 views • link
written by mittu1602 mittu1602

Please explain better with some sample input and output

• 0 views • link
written by Kevin Blighe Kevin Blighe

File1:

 rs492602   **CC**  Vitamin B12 deficiency  FUT2    Higher levels of vitamin B12        
    rs492602    CT  Vitamin B12 deficiency  FUT2    Normal levels of vitamin B12        
    rs492602    TT     Vitamin B12 deficiency   FUT2    Normal levels of vitamin B12

File2:

    rs492602    **GG**  exm-rs492602
    rs492602    **CC**  exm-rs492602

Result:

rs492602**CC** **GG** Vitamin B12 deficiency    FUT2    Higher levels of vitamin B12    
rs492602**CC** **CC** Vitamin B12 deficiency    FUT2    Higher levels of vitamin B12
• 0 views • link
written by mittu1602 mittu1602

Are those asterisks in the actual data or did you just put them there for this example?

• 0 views • link
written by Kevin Blighe Kevin Blighe

for this example only, wanted to bold them to stand out!

• 0 views • link
written by mittu1602 mittu1602

I see. Yes, the formatting doesn't always come out right on these forums.

You can try this:

cat file1 
rs492602    CC  Vitamin B12 deficiency  FUT2    Higher levels of vitamin B12
rs492602    CT  Vitamin B12 deficiency  FUT2    Normal levels of vitamin B12
rs492602    TT  Vitamin B12 deficiency  FUT2    Normal levels of vitamin B12

cat file2
rs492602    GG  exm-rs492602
rs492602    CC  exm-rs492602

paste <(cut -f2 file1 | tr [ATGC] [TACG]) file1 | awk 'BEGIN {FS="\t"} FNR==NR {key1=$2$3; key2=$2$1; arrayVars[key1]=$0; arrayVars[key2]=$0; next} {key=$1$2; if (arrayVars[key]) split(arrayVars[key], var, "\t"); if (arrayVars[key]) print var[2]"\t"var[3]"\t"$2"\t"var[4]"\t"var[5]"\t"var[6]}' - file2
rs492602    CC  GG  Vitamin B12 deficiency  FUT2    Higher levels of vitamin B12
rs492602    CC  CC  Vitamin B12 deficiency  FUT2    Higher levels of vitamin B12

Here I find the complementary base with the first part of the pipe with paste <(cut -f2 file1 | tr [ATGC] [TACG]) file1, and this is then piped into awk as stdin. If there are also more columns to print out, just add more var[2]"\t" to the final print command.

Hope that this works!

• 0 views • link
written by Kevin Blighe Kevin Blighe

Thank you so much @Kevin it really worked and helped me. Much appreciated for your time over this! Can u also suggest me that if I can use it as a shell script? Happy New Year!

• 0 views • link
written by mittu1602 mittu1602

Yes this should be fine within a shell script (bash)

• 0 views • link
written by Kevin Blighe Kevin Blighe

Log in to answer this question.

More posts like this