• Welcome to Phoenix Rising!

    Created in 2008, Phoenix Rising is the largest and oldest forum dedicated to furthering the understanding of and finding treatments for complex chronic illnesses such as chronic fatigue syndrome (ME/CFS), fibromyalgia (FM), long COVID, postural orthostatic tachycardia syndrome (POTS), mast cell activation syndrome (MCAS), and allied diseases.

    To become a member, simply click the Register button at the top right.

Big Data App to Explore Genomes for Clinical Relevance, Rare Variants, Drug Response, etc (Free)

Moof

Senior Member
Messages
778
Location
UK
Would anyone be interested in doing a patient collaborative study looking at whole genomes or am I a lone nutter? I have a small list candidate genes/risks that I want to verify prevalence of and also want to see if the algorithms can find other risks. I have also developed algorithms that can look at all copy number variants and structural variants in a genome as well. The algorithm could automatically assess loss of function, missense, CADD scores, Polyphen/SIFT deletriousness scores, etc.

Well, you have access to my info...you're welcome to experiment with it! I can't remember whether I uploaded the CNV and SF variant files, but I can do so if needed.
 

kday

Senior Member
Messages
369
Update rolled out with support for rare and uncommon SNPs that match ref and general fixes of bugs.

Made indel calling a little more liberal with a little magic so if you don't completely match either the ref or alt allele, the indel may show anyway! This increased sensitivity has seemed ok with my testing, but this increased sensitivity can be at the expense of specificity.

This indel algorithm is important, because quite a few Indels can be missed with typical techniques since the alleles can vary and aren't always consistent for the same indel!
 

Rufous McKinney

Senior Member
Messages
13,249
@kday

I re-entered my DNA and was able to save the reports this time...thank you.

So mine says: hetero for Ehlers-Danlos syndrome 7A

If that is autosomal dominant, then: in theory I have this.

I've got "issues"..but I don't "have" full bore EDS.

so: I still don;t really understand what this stuff means. Or why would it say Likely benign- when it also says: you have it.

Since I am collagen obsessed at the moment, I'd like to understand this: should I start a thread? Does somebody else understand this?
 

kday

Senior Member
Messages
369
@kday

I re-entered my DNA and was able to save the reports this time...thank you.

So mine says: hetero for Ehlers-Danlos syndrome 7A

If that is autosomal dominant, then: in theory I have this.

I've got "issues"..but I don't "have" full bore EDS.

so: I still don;t really understand what this stuff means. Or why would it say Likely benign- when it also says: you have it.

Since I am collagen obsessed at the moment, I'd like to understand this: should I start a thread? Does somebody else understand this?
If it says it's likely benign, then it's likely benign. However, the submitters that classified it could have been wrong. If there is only 1 submission and no research, you can't know for sure. But if the CADD score is over 20 without 1 submssion and no research, I would treat it as a suspect and bookmark it to see what happens with further research. If it's over 15, it's a potential suspect in my eyes, but variants with a score over 20 usually bear more weight from my observations. But that's not always the case as CADD scores can be off.

Also look at frequency. If it's not very rare, it probably isn't a causative variant. I personally haven't seen a true case of genetic EDS yet in the genomes I've looked at.
 
Last edited:

Judee

Psalm 46:1-3
Messages
4,461
Location
Great Lakes
Ah, nevermind my previous post. Perhaps my mind is getting too creative. I guess I am just sick of things moving so slow. It's been a decade for me and I feel like I nearly gave up completely not too long ago.

Please don't give up. I'm not techy enough to understand all of what you're doing but it seems like some very good and tremendously important work. Dream big if you can. My faith, family and dreams (that I can't live out right now) are all I have left.

It's been almost 40 years for me and I keep hitting new baseline lows. Just when I think it can't get any worse, it does. :(

So for yourself and please, for us...Dream Big!!! (just pace yourself a bit so you don't completely wear out.)
 

Moof

Senior Member
Messages
778
Location
UK
I re-uploaded my indel file, which worked fine. None showed up at all, so that's one less thing to worry about!
 

nanonug

Senior Member
Messages
1,709
Location
Virginia, USA
Pretty cool tool, thanks!

Interesting that all of the identified "clinically significant" variants are mitochondrial.
 

Attachments

  • Nanonug Variant Report.pdf
    662.8 KB · Views: 39

Mel9

Senior Member
Messages
995
Location
NSW Australia
I have tried to get whole genome done but the lab in Italy caused some problems and after three months the packaged returned to me ‘unclaimed’

They sent me another kit
 

kday

Senior Member
Messages
369
@nanonug. That report is completely wrong. I believe the issue is that you uploaded a file generated by 23andMe many years ago using a different genome assembly.

Can you confirm the date in your file by unzipping it and sending me the header of the file. Because if so, you need to download an updated file at https://you.23andme.com/tools/data/download/

I also need to make an error for very old files with a link to download an updated version.
 

wigglethemouse

Senior Member
Messages
776
Because if so, you need to download an updated file at https://you.23andme.com/tools/data/download/
Even with more recent files, some MT locations were offset by 1 and have have since been corrected. Well worth everyone with 23andMe data downloading a recent file before trying your app. It's from the same sample as before, but just cleaned up data as 23andMe learns more about mistakes in their chips or rsID's change.
 

Seven7

Seven
Messages
3,444
Location
USA
It's designed to work very well with such files.
THank you! bizare my report from them, said I had a lot of drug interactions and here I do not or any genetic condition either even though I think the EDS gene is mutated
 

kday

Senior Member
Messages
369
Thank you for sharing your app @kday

Is there any web you can get the HLA variants with your VCF file ?
HLA is a complex topic, and is a bit over my head, but the short answer is yes.

My app does not calculate HLA types, but they can be calculated from WGS data. There are tools that do this.

edit: I think specifically, you need to FASTQ files to calculate HLA genotype. At this point in time, an application on your local machine is better suited than a web application because of the amount of data.
 
Last edited: