6
$\begingroup$

I have always thought that 1/0 and 0/1 in VCF genotype fields are equivalent. And yet, GATK uses both. For example, these are two variants called in the same sample and the same run of GATK 4.1.4.0:

chr7    117120317   .   ATTCATTGTTTTGAAAGAAAGATGGAAGAATGAACTGAAG    A   748.97  .   AC=1;AF=0.5;AN=2;DP=64;ExcessHet=3.0103;FS=0;MLEAC=1;MLEAF=0.5;MQ=60;QD=11.89;SOR=7.223 GT:AD:DP:GQ:PL:SB   1/0:0,36:63:99:2294,1042,933:0,0,0,36
chr7    117120306   .   GACGCGCT    G   516.6   .   AC=1;AF=0.5;AN=2;BaseQRankSum=0.826;DP=70;ExcessHet=3.0103;FS=0;MLEAC=1;MLEAF=0.5;MQ=60;MQRankSum=0;QD=8.2;ReadPosRankSum=-5.786;SOR=0.976  GT:AD:DP:GQ:PL:SB   0/1:37,26:63:99:524,0,1203:0,20,0,26

And here are the same lines redacted for clarity:

chr7    117120317   .   AT..AG  A   [...]   GT:AD:DP:GQ:PL:SB   1/0:0,36:63:99:2294,1042,933:0,0,0,36
chr7    117120306   .   GA..CT  G   [...]   GT:AD:DP:GQ:PL:SB   0/1:37,26:63:99:524,0,1203:0,20,0,26

As you can see, one is genotyped as 1/0 and the other as 0/1. The VCF specs say:

GT (String): Genotype, encoded as allele values separated by either of / or |. The allele values are 0 for the reference allele (what is in the REF field), 1 for the first allele listed in ALT, 2 for the second allele list in ALT and so on. For diploid calls examples could be 0/1, 1 | 0, or 1/2, etc. Haploid calls, e.g. on Y, male non-pseudoautosomal X, or mitochondrion, are indicated by having only one allele value. A triploid call might look like 0/0/1. If a call cannot be made for a sample at a given locus, ‘.’ must be specified for each missing allele in the GT field (for example ‘./.’ for a diploid genotype and ‘.’ for haploid genotype). The meanings of the separators are as follows (see the PS field below for more details on incorporating phasing information into the genotypes):

  • / : genotype unphased
  • | : genotype phased

This suggests that although 1|0 is not the same as 0|1, the unphased 0/1 and 1/0 genotypes are completely equivalent. Which makes perfect sense given that the N/N genotypes are unphased, so n/N must be equivalent to N/n.

So, my questions are: why does gatk produce some lines with 0/1 and others with 1/0? Is it a bug? Am I missing something? I can't imagine a developer writing code to pick either 0/1 or 1/0 randomly, so presumably they indicate some difference in what gatk sees about a variant. So are 0/1 and 1/0 not 100% equivalent somehow?


I have seen How is the GT field in a VCF file defined? which suggests the two are equivalent, but this seems to be contradicted by GATK's own output which is why I am asking this again to get a clear and preferably referenced answer.

$\endgroup$
3
  • $\begingroup$ @M__ this was from a human sample. But you can't know dominant or recessive unless you also know what gene, if any, is affected by the variant and what that gene's mode of inheritance is with respect to a specific disease, so I don't think it can be about AD vs AR. $\endgroup$
    – terdon
    Mar 10, 2022 at 16:20
  • $\begingroup$ Okay finally got it. The convention is it should be in ascending ordering, i.e. 0/1. The question is, why the convention broken? $\endgroup$
    – M__
    Mar 23, 2022 at 15:48
  • 1
    $\begingroup$ I just realized that's almost certainly on me. I just took a couple of variants from a file randomly, and I hadn't really noticed the positions. I'm pretty sure they were in the right (ascending) order in the original file. Which lends more credence to your suggestion, @M__. $\endgroup$
    – terdon
    Mar 23, 2022 at 21:23

1 Answer 1

2
$\begingroup$

Convention The convention is it should be in ascending ordering, i.e. 0/1. The question is, why the convention broken?

Phasing I think, at a guess, the answer is it depends what reference genome you are using. I would suggest that this result would not occur if the reference genome is mum (or dad) and the alleles are occurring in the child. The output supports this idea because it is marked with '/', i.e. they are not phased, thus the reference genome is neither mum nor dad.

Example provided The thing I would like to point out is in the example given 117120306 and 117120317 are not randomly scattered across the genome but appear to be 11 bp apart. For a human genome those are very strongly linked.

Linkage hypothesis Thus I would suggest the answer is linkage in this example, given the reference genome isn't a parent. The mutations are occurring on different strands therefore better represented as 0/1 and then 1/0, i.e. one strand comprises 'A' allele and the opposite strand is the 'G' allele.

Why is unphased linkage important In this scenario the interpretation could be hugely important because it would be stating that each allele has come from different parents. This would be more important than maintaining convention. In an inherited disorder there is a risk the child has a dominant phenotype, when each parent could have been heterozygous. Clearly it depends whether the mutation is forcing an amino acid mutation, at 11bp it is certainly possible, i.e. mutation in the first codon upstream and mutation in the second codon downstream.

More generally in any population analysis, linkage is an enormously important issue and central to any calculation. The point is simply that most analysis infers independent segregation of each allele and in this example that remains a good assumption because the allele (I am suggesting) are on opposite strands. If they were on the same strand it would create a bias away from panmixia.

The hypothesis is formed on this example alone and if the alleles were 10 megabp apart it would be impossible to conclude they were linked, therefore the convention of 0/1 should hold. That's my guess, its not a bad guess.

$\endgroup$
5
  • 1
    $\begingroup$ This is an interesting suggestion. I just picked two random variants from the file, so I will have to investigate a bit to see if what you suggest is actually the case more generally. I'd love to see this backed up by a gatk source (post, or dev, or source file or whatever). Certainly plausible, anyway, thanks! $\endgroup$
    – terdon
    Mar 23, 2022 at 16:29
  • $\begingroup$ Hope it helps @terdon. Its a hypothesis to work on. Geraldine@Broad will know of course. $\endgroup$
    – M__
    Mar 23, 2022 at 16:41
  • 1
    $\begingroup$ I'm honored by your confidence in me, @M__ — But in this case I don't know for sure 😅 I don't recall seeing something like this before. And if it was a phasing thing I would expect there to be phase related annotations. Are you seeing this in more than one record? And is this single or multi sample, have you done any subsetting, etc? FWIW you might get a definitive answer faster on the GATK forum. $\endgroup$ Mar 28, 2022 at 17:35
  • $\begingroup$ Alas, but regardless my confidence remains @Geraldine_VdAuwera 🙂 Thanks. $\endgroup$
    – M__
    Mar 29, 2022 at 19:57
  • $\begingroup$ Thanks, @Geraldine_VdAuwera, I'll try there. $\endgroup$
    – terdon
    Mar 30, 2022 at 11:55

Your Answer

By clicking “Post Your Answer”, you agree to our terms of service and acknowledge that you have read and understand our privacy policy and code of conduct.

Not the answer you're looking for? Browse other questions tagged or ask your own question.