Package: arv
arv — a fast 23andMe parser for PythonArv (Norwegian; “heritage” or “inheritance”) is a Python module for parsing raw
23andMe genome files. It lets you lookup SNPs from RSIDs.fromarvimportload,unphased_matchasmatchgenome=load("genome.txt")print("You are a{gender}with{color}eyes and{complexion}skin.".format(gender="man"ifgenome.y_chromosomeelse"woman",complexion="light"ifgenome["rs1426654"]=="AA"else"dark",color=match(genome["rs12913832"],{"AA":"brown","AG":"brown or green","GG":"blue"})))For my genome, this little program produces:You are a man with blue eyes and light skin.The parser is insanely fast, having been written in finely tuned C++, exposed
via Cython. A 2013 Xeon machine I’ve tested on parses a 24 Mb file into a hash
table in about 78 ms. The newer 23andMe files are smaller, and parses in a mere
62 ms!Works with Python 2.7+ and 3+. Installable with pip!$pipinstall--upgradearvSee below for software requirements.Important disclaimerIt’s very important to tell you that I, the author of arv, am merely ahobbyist! Iama professional software developer, butnota geneticist,
biologist, medical doctor or anything like that.Because of that, this software may not only look weird to people in the field,
it may also contain serious errors. If you find any problem whatsoever, please
submit a GitHub issue.This a slightly modified version of what I wrote for the original software
called “dna-traits”, and the same goes for this software:In addition to the GPL v3 licensing terms, and given that this code deals with
health-related issues, I want to stress that the provided code most likely
contains errors, or invalid genome reports. Results from this code must be
interpreted as HIGHLY SPECULATIVE and may even be downright INCORRECT. Always
consult an expert (medical doctor, geneticist, etc.) for guidance. I take NO
RESPONSIBILITY whatsoever for any consequences of using this code, including
but not limited to loss of life, money, spouses, self-esteem and so on. Use at
YOUR OWN RISK.The indended use is for casual, educational purposes. If this code is used for
research purposes, please cross-check key results with other software: The
parser code may contain serious errors, for example.An interesting story about the research part: I once released a pretty good
Mersenne Twister PRNG for C++ that ended up being used in research. Turned out
the engine had bugs, and by the time I had fixed them, a poor researcher had
already produced results with it (hopefully not published; I don’t know). The
guy had to go back and fix his stuff, and I felt terribly bad about it.So beware!InstallationThe recommended way is to install from PyPi.$pipinstallarvThis will most likely build Arv from source. The package will automatically
install Cython, but it doesn’t check if you have a C++11 compiler. Furthermore,
it passes some additional compilation flags that are specific to clang/gcc.If you have problems runningpip install arv, please open an issue on
GitHub with as much detail as possible (g++/clang++--version,uname-a,python--versionand so on).If you set the environment variableARV_DEBUG, it will build with full
warnings and debug symbols.You can also install it locally throughsetup.py. The following builds and
tests, but does not install, arv:$pythonsetup.pytestIf you set the environment variableARV_BENCHMARKto a genome filename and
run the tests, it will perform a short benchmark, reporting the best parsing
time on it. You can also setARV_BENCHMARK_COUNT=<number>to change how
many times it should parse the given file.UsageFirst you need to dump the raw genome file from 23andMe. You’ll find it under
the raw genome browser, and download the file. You may have to unzip it first:
The parser works on the pure text files.Then you load the genome in Python with>>>genome=arv.load("filename.txt")>>>genome<Genome:SNPs=960613,name='filename.txt'>To see if there are any Y-chromosomes present in the genome,>>>genome.y_chromosomeTrueThe genome provides adict-like interface. To get a given SNP, just enter the RSID.>>>genome["rs123"]>>>snp<SNP:chromosome=7position=24966446genotype='AA'>>>>snp.chromosome7>>>snp.position24966446>>>snp.genotype<Genotype'AA'>TheGenotypeobject can be converted to a string withstr, but it also
allows rich comparisons with strings directly:>>>snp.genotype=="AA"Trueyou can get its complement with the~-operator.>>>type(snp.genotype)<class'_arv.Genotype'>>>>~snp.genotype<Genotype'TT'>The complement is important due to eah SNPs orientation. All of 23andMe SNPs
are oriented towards the positive (“plus”) strand, based on theGRCh37reference human genome assembly
build. But some SNPs on SNPedia are given with theminus orientation.For example, to determine if the human in question is likely lactose tolerant
or not, we can look atrs4988235.
SNPedia reports itsStabilizedorientation to be minus, so we need to use the
complement:>>>genome["rs4988235"].genotype<Genotype'AA'>>>>~genome["rs4988235"].genotype<Genotype'TT'>By reading a fewGWASresearch
papers, we can build a rule to determine a human’s likelihood for lactose
tolerance:>>>arv.unphased_match(~genome["rs4988235"].genotype,{"TT":"Likely lactose tolerant","TC":"Likely lactose tolerant","CC":"Likely lactose intolerant",None:"Unable to determine (genotype not present)"})'Likely lactose tolerant'Note that reading GWAS papers for hobbyists can be a bit tricky. If you are a
hobbyist, be sure to spend some time reading the paper closely, checking up
SNPs on places likeSNPedia,dnSNPandOpenSNP. Finally, have fun, but be extremely careful
about drawing conclusions from your results.Command line interfaceYou can also invokearvfrom the command line:$python-marv--helpFor example, you can drop into a Python REPL like so:$python-marv--replgenome.txtgenome.txt...960614SNPs,maleType`genome`toseetheparsed23andMerawgenomefile>>>genome<Genome:SNPs=960614,name='genome.txt'>>>>genome["rs123"]<SNP:chromosome=7position=24966446genotype=<Genotype'AA'>>If you specify several files, you can access them through the variablegenomes.The example at the top of this document can be run with--example:$python-marv--examplegenome.txtgenome.txt...960614SNPs,malegenome.txt...AmanwithblueeyesandlightskinLicenseCopyright 2017 Christian Stigen LarsenDistributed under the GNU GPL v3 or later. See the file COPYING for the full
license text. This software makes use of open source software; see LICENSES for
details.