Validate PCR DNA primers against whole de novel genome
1
0
Entering edit mode
2.3 years ago
qwzhang0601 ▴ 70

I want to test gene copy number of some genes. I designed the primer using primer3 and what to test whether the primers are specific to my genes. I know some tools can do this. For example, I tried primer-blast, but it can not work when I test it against the whole genome (about 3G base pairs). I also know others (e.g., UCSC In-Silico PCR), but since our genome is new, I can not provide my genome. Do anyone knows some tools can handle large genome to check specificity of primers?

Thanks

primers in sillico • 566 views
ADD COMMENT
1
Entering edit mode
2.3 years ago
trausch ★ 1.6k

Dicey can do this. You first need to build an FM-Index and the standard FASTA index of your genome

dicey index -o new_genome.fa.fm9 new_genome.fa.gz

samtools faidx new_genome.fa.gz

Then you can search primer pairs across the entire genome.

echo -e ">Primer_forward\nGCCCCATAGGTTTTGAACTCA\n>Primer_reverse\nTGATTTGTCTGTAATTGCCAGC" > primers.fa

dicey search -c 45 -g new_genome.fa.gz primers.fa

The default output is a JSON file that you can convert to text using the provided script (scripts/json2txt.py). Dicey is available as a static binary or on Bioconda.

ADD COMMENT
0
Entering edit mode

Many thanks. Is there a publication and manual for dicey? I am looking at MFEprimer-3.0 (https://www.mfeprimer.com/), have you tried this before?

Thank you

ADD REPLY
0
Entering edit mode

No publication yet, there is only the GitHub README. But I am happy to help if you encounter problems with dicey. Just open an issue in the repo.

ADD REPLY
0
Entering edit mode

Ok. Thank you very much.

ADD REPLY

Login before adding your answer.

Traffic: 1138 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