3
$\begingroup$

I bought a "Whole Genome Sequence" kit, which gave me the following VCF files (this service did not include any analysis, just the raw data):

cnv.vcf indel.vcf snp.vcf sv.vcf

My understanding is that these files show the differences of my DNA relative to a reference human genome (hg19 in my case)

What tools and resources can I use to find out which known alleles/variants of genes I have?

An example: Blood type

According to Wikipedia, there are 6 common alleles of the ABO gene, and which alleles I have determine my blood type.

Looking at my VCF results, I have these SNPs and INDELs within that gene:

<removed>

Given this, how would I figure out my blood type?

I understand that my data may not be 100% accurate, and that many phenotypes are the result of more than one gene, environmental factors, or are not well studied at all, but this should work for some simple things ie blood type, COMT, BHLHE41, CYP1A2, etc

Improve this question
$\endgroup$
7
  • $\begingroup$ For the blood locus. Try and find online the genomic position (i.e. the number like 133273540 in the POS field of the .vcf) of the ABO blood locus and look it up on your dataset. Remember to make sure you find it on the gr37/hg19 build. $\endgroup$
    – user438383
    Jan 27, 2021 at 10:07
  • $\begingroup$ Wait, did you got yourself yourself sequenced?? I thought "my DNA" means "my DNA sample" not "my own DNA"! $\endgroup$
    – Kamil S Jaron
    Jan 27, 2021 at 12:43
  • 3
    $\begingroup$ Please don't just share your own genetic information online, especially if you are not fully understanding what you are looking at. It is not something to be posted everywhere, although some loci might be relatively harmless you shouldn't take this lightly. $\endgroup$
    – Wouter De Coster
    Jan 27, 2021 at 13:26
  • 3
    $\begingroup$ Please ping me if you decide you want this censored. You should be very aware of the long-term consequences for you and your current and future relatives of sharing your genetic information on the internet in a post. $\endgroup$
    – Devon Ryan
    Jan 27, 2021 at 13:35
  • 1
    $\begingroup$ @WouterDeCoster also thanks for the warning $\endgroup$
    – Knio
    Jan 27, 2021 at 19:08

1 Answer 1

Reset to default
1
$\begingroup$

If this is your own DNA and you are trying to analyse yourself you should talk to a genetic counselor. I personally think it's cool, but I deeply understand the potential consequences.

What you want is a variant annotation tool. There are many, but I normally write my own. Which is simply a case of pulling down a database and using pandas to match variants with annotation.

There are many databases. The biggest popular database is gnomad and the biggest disease database is clinvar. Others include omim, snpdb and 1k genomes. I'm not sure about blood sorry but I'm sure you can find one.

Improve this answer
$\endgroup$

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.