Question: Uses for low-coverage human genome sequence
1
gravatar for Will
4.2 years ago by
Will4.5k
United States
Will4.5k wrote:

I have a cohort of patients that we are doing HIV deep sequencing of integrated provirus and because of the technical methods we are using we also get some of host (human) genomc sequence ... I'm estimating 2 billion bp, so ~1X coverage. Its not going to be enough coverage to capture any clinically relevant human SNPs (unless we get extremely lucky).

The data is a combination of Illumina NextSeq and PacBio reads (mostly Illumina) on what will ultimately be >500 patients.

I already know what I'm going to do with the HIV sequence but I'm looking to squeeze something interesting out of this "leftover" data.

Any ideas? 

genome • 1.3k views
ADD COMMENTlink modified 4.2 years ago by mikhail.shugay3.3k • written 4.2 years ago by Will4.5k

Check out this: http://www.ncbi.nlm.nih.gov/pmc/articles/PMC3400344/

ADD REPLYlink written 4.2 years ago by lh331k

Hmm, that seemed to be exome sequencing (or at least from my 60 second scan). We are getting whole genome sequence, there's no way to enrich for exomes. But maybe using this and pooling would work.

ADD REPLYlink written 4.2 years ago by Will4.5k

The main point of that paper is not about exome sequencing. It is saying that you can call pretty good variants from off-target reads with imputation. The off-target reads are similar to your "undesired" human reads at ~1X coverage. There is another post today about using off-target reads for CNV calling.

ADD REPLYlink modified 4.2 years ago • written 4.2 years ago by lh331k
0
gravatar for mikhail.shugay
4.2 years ago by
mikhail.shugay3.3k
Czech Republic, Brno, CEITEC
mikhail.shugay3.3k wrote:

If you have a large cohort you can do pooling and compare minor alleles for patient groups. Calling SNPs in pooled samples would require to alter variant frequency model assumptions, please consider this and similar papers.

ADD COMMENTlink written 4.2 years ago by mikhail.shugay3.3k

That's pretty cool. We could almost do a mini case/control GWAS study ... not sure if it would work but its worth a try.

ADD REPLYlink written 4.2 years ago by Will4.5k
Please log in to add an answer.

Help
Access

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 924 users visited in the last hour