Daniel:Notebook/ComboLock/CProbeDesign: Difference between revisions

From ZhangLabWiki
Jump to navigation Jump to search
>Djacobse
>Djacobse
Line 2: Line 2:


[[Daniel:Notebook/ComboLock|Back to Notebook]]
[[Daniel:Notebook/ComboLock|Back to Notebook]]
*[http://www.nature.com/nmeth/journal/v13/n3/abs/nmeth.3742.html Frei et al paper] C probe desing is based on


This page describes the methods used for designing C-Probes for Combo Lock.  The illustration below shows what a C-probe is expected to look like.  The descriptions given include the original design sizes.
This page describes the methods used for designing C-Probes for Combo Lock.  The illustration below shows what a C-probe is expected to look like.  The descriptions given include the original design sizes.

Revision as of 16:40, 21 April 2016

C Probe Design

Back to Notebook

This page describes the methods used for designing C-Probes for Combo Lock. The illustration below shows what a C-probe is expected to look like. The descriptions given include the original design sizes.

File:C Probes.png

Designing Transcript Matching Region

    Must install R, Primer3, and Blast+ before using
  1. Make blast databases
  2. make one for repetitive regions (from Repbase)
  3. makeblastdb -in repbase.fa -dbtype nucl
  4. Make one for transcriptome; use only high quality bases (from NCBI ftp)
  5. makeblastdb -in rna_human_high_qual.fa -parse_seqids -dbtype nucl
  6. In R:
  7. library(devtools) install_github("nolanlab/PLAYRDesign") library(PLAYRDesign)
  8. Next make a library from UCSC in R:
  9. library(GenomicFeatures) txdb <- makeTxDbFromUCSC(genome = "hg19", tablename = "refGene") saveDb(txdb, "PUT THE PATH TO THE OUTPUT FILE HERE")
  10. Finally, run PLAYRDesign
    1. Select a file (any file) in the working directory you wish to use
    2. You must have a configuartion file in the directory you wish to use (see below)
    3. You can hit escape any time in the R window to stop PLAYRDesign
  11. PLAYRDesign.run()

PLAYRDesign Configuration File

Needed to run PLAYRDesign. The name of the file must be playrdesign_conf.txt and must be in your working directory. An example is below. Note that if there are any spaces in your file (including after the path name) PLAYRDesign will not work properly

BLASTN_EXEC=/usr/bin/blastn           (The full path to the blastn executable)
BLASTN_DB=/opt/BLAST/                 (The directory containing your BLAST database files)
PRIMER3_EXEC=/usr/bin/primer3_core    (The full path to the primer3_core executable)
PRIMER3_CONFIG=/opt/primer3_config/   (The primer3_config directory that is found in the primer3 distribution, see above)
PLAYRDESIGN_DATA=/opt/PLAYRDesign_data   (The directory containing the EST and exon data, see above)

PLAYRDesign Output

I have attached a sample of the output probes from the C-probe designer: C Probe Analysis SOD1/VIM (pdf)

C Probe Analysis SOD1/VIM (docx)


Looking at the file, probes, I have highlighted the sequences that match the mRNA (mRNA sequence used is below the probe). From looking at two different pairs of probes for SOD1 and VIM (8 probes total), it is clear that the first 24bp are genomic matching. There then appears to be a 10 base poly-A spacer region, followed by the relevant barcodes. I have used the 24bp transcript matching region and 10bp spacer in my own designs.