Question

Why are there nucleotides in the reference genome which are also minor alleles in the 1000 genomes project

0

Entering edit mode

2.8 years ago

nhaus ▴ 300

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 • 855 views

ADD COMMENT • link updated 2.8 years ago by WouterDeCoster 47k • written 2.8 years ago by nhaus ▴ 300

score 2 · Accepted Answer · 2021-06-30

2

Entering edit mode

2.8 years ago

WouterDeCoster 47k

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.

ADD COMMENT • link 2.8 years ago by WouterDeCoster 47k

0

Entering edit mode

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.

ADD REPLY • link 2.8 years ago by i.sudbery 19k

0

Entering edit mode

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

ADD REPLY • link 2.8 years ago by WouterDeCoster 47k

0

Entering edit mode

Thanks! Big misconception on my side

ADD REPLY • link 2.8 years ago by nhaus ▴ 300