This is a test version of Biostars. For the public version, visit https://www.biostars.org.
Why are there nucleotides in the reference genome which are also minor alleles in the 1000 genomes project

Hello,

I was wondering if you could answer a conceptional question for me.

I have noticed, that there are quite some SNPs, where the reference nucleotide is not the most common nucleotide in the majority of populations.

One such example would be the SNP rs13303010. If I take a look at this SNP in ensembl, it says that "G" is the reference nucleotide and "A" is the alternative nucleotide. This is the case for both GRCh37 and GRCh38. When looking at the population genetics, it is clear that for all populations except the African, the A nucleotide is much more common. The "G" reference allele only has a allele frequency of 37% across all populations.

So my question is, why the reference assembly (GRCh38) was not changed at that position when moving from GRCh37 to GRCh38.

I thought the reference assembly is supposed to be somewhat like a consensus sequence of all humans and therefore "A" should be put as the reference because "A" has the highest allele frequency.

I would be very thankful if somebody could explain to me why the reference is not changed in such a case.

Cheers.

reference genotype minor allele
• 1,787 views • link
updated by WouterDeCoster WouterDeCoster • written by nhaus • nhaus

1 answer

I thought the reference assembly is supposed to be somewhat like a consensus sequence of all humans and therefore "A" should be put as the reference because "A" has the highest allele frequency

Well, no, this is totally not the case. The reference genome is a mosaic, collapsed haploid, derived from multiple anonymous donors (some more than others), but at every location, just one contig is chosen. There is no consensus made.

That is also neither straightforward nor correct. You would for example create haplotypes that do not exist.

• 0 views • link
written by WouterDeCoster WouterDeCoster

While this is true, there was an attempt to make more of the reference alleles the major allele in GRCh38 than in previous genome versions.

• 0 views • link
written by i.sudbery i.sudbery

Oh yes, some have been corrected, but you absolutely shouldn't rely on that.

• 0 views • link
written by WouterDeCoster WouterDeCoster

Thanks! Big misconception on my side

• 0 views • link
written by nhaus • nhaus

Log in to answer this question.

More posts like this