aneuploidy detection in DNA sequencing
1
0
Entering edit mode
7.8 years ago

I am planning to perform DNA seq for detection of aneluploidy condition in circulatory DNA. Is there a set rule how much sequencing depth should be used. Secondly I read few papers but could not find a set software to detect aneluploidy from DNA sequencing. Is there a set software or tool that can be useful. Any suggestion will be helpful. I am working with human genome and using Illumina nex seq platform for whole genome DNA seq.

Thanks

sequencing • 2.4k views
ADD COMMENT
0
Entering edit mode

Please add information with regard to the organism you are working on and the sequencing technology you are planning to use.

ADD REPLY
0
Entering edit mode

My Bad. I have edited information in original post.

ADD REPLY
0
Entering edit mode
7.8 years ago
willgilks ▴ 360

Hi joerodger2017

I'm not sure what the minimum depth you use is. Maybe you could sequence the same individual at different depths as a pilot study.

To count the number of reads mapping to each chromosome, you can use the command: samtools idxstats alignments.bam (http://samtools.sourceforge.net/). Divide this number by the length of each chromosome to get a reads-per-base for each chromosome. This value is the same across all autosomes of an individual (low variance) ... unless there is a gain or loss of material. Sex chromosome dosage could be incorporated.

To detect sub-chromosomal structural variation, you could use the CNV-detection program, Genomestrip, except adjust the settings so that only really large events are screened-for, http://software.broadinstitute.org/software/genomestrip/ . I used it more generally with open-acces code, here https://f1000research.com/articles/5-2644/v3.

ADD COMMENT
0
Entering edit mode

I'm not sure what the minimum depth you use is. Maybe you could sequence the same individual at different depths as a pilot study.

Less expensive solution: sequence an individual once and use in silico downsampling to estimate required/sufficient depths.

ADD REPLY

Login before adding your answer.

Traffic: 2404 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6