SEQanswers

Go Back   SEQanswers > Bioinformatics > Bioinformatics



Similar Threads
Thread Thread Starter Forum Replies Last Post
Best practice for preprocessing the bam file for copy number estimation ty23991 Bioinformatics 2 05-20-2016 05:48 AM
Normalize deletion number from WGS papori Bioinformatics 0 11-25-2014 12:13 AM
DNA Copy Number arcolombo698 Bioinformatics 0 07-14-2014 02:30 PM
Copy number variation in cancer WGS lethalfang Bioinformatics 31 05-14-2014 11:55 AM
Overlap number discrepancy between VCFTools and BEDTools zxyeo Bioinformatics 2 12-27-2011 05:36 PM

Reply
 
Thread Tools
Old 10-25-2016, 09:45 PM   #1
peerah
Junior Member
 
Location: USA

Join Date: Jan 2015
Posts: 6
Unhappy Discrepancy in plasmid and transposons copy number estimation using WGS vs qPCR

We are trying to estimate plasmid copy number in a bacterium by mapping Illumina reads on corrected PacBio contigs and comparing read coverage between the chromosomal and plasmid contigs. Chromosome is about 5 Mbp and plasmid is 190 kbp.

The depth was good for both contigs (200-1000 reads). The ratio between chromosomal vs plasmid was about 1:2. The coverage of a transposon in the plasmid was even higher (1:4 comparing to chromosomal). I tried to confirm this using qPCR by comparing PCR efficiency-adjusted Ct (not the typical ddCT that assumes all efficiency = 2) between two reference single-copy chromosomal genes (rpoB, gyrA) and three target genes on the plasmid (two inside a transposon and one outside), and everything came down to a tight and solid 1:1. It is not uncommon for a plasmid like this to be single copy. I was wondering if the sequencing depth could be used in this manner without any bias adjustment. I am seriously doubting my ability to design a PCR, but I found nothing wrong with it. Should I trust the simple coverage calculation or the PCR at this point?

The average GC content of chromosome is about 56%. Plasmid GC content is 54%, but the transposon on the plasmid (10 kb) is 62% GC.

Any suggestion would be highly appreciated. Thank you so much!
peerah is offline   Reply With Quote
Reply

Thread Tools

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off




All times are GMT -8. The time now is 01:39 PM.


Powered by vBulletin® Version 3.8.9
Copyright ©2000 - 2018, vBulletin Solutions, Inc.
Single Sign On provided by vBSSO