Daniel:Notebook/ComboLock/CProbeDesign: Difference between revisions
>Djacobse |
>Djacobse |
||
Line 9: | Line 9: | ||
[[Image:C_Probes.png|600px]] | [[Image:C_Probes.png|600px]] | ||
==Designing Capture Arm== | |||
The capture arm design aims to create 2 sequences that are orthogonal to the human genome/transcriptome, so that these sequences will not have non-specific binding. The two flanking sequences (one for padlock and one for latch) will have an 8bp region between them that will be a barcode unique to each C-probe. | |||
Each C probe of a pair will have: | |||
#Padlock matching region, this will be a different sequence for each of the pair of C probes, but will be the same across pairs | |||
#Latch matching region, this will be a different sequence for each of the pair of C probes, but will be the same across pairs | |||
#8 bp barcode region, this will be unique for every C probe | |||
<ol> | |||
<li>Randomly generate a large list of 18bp sequences</li> | |||
python primergenerator.py -a -n $nprimers -l $plength -o primers.fa | |||
<li>BLAST sequences against genome</li> | |||
blastn -db hg19.fa -query primers.fa -num_descriptions 5 -task blastn-short -out primers.blastout.txt -outfmt 7 | |||
<li>Reorder the sequences based on BLAST results (oligos with fewest hits are best)</li> | |||
python primerblast.analyzer.py -n $nprimers -l $plength -o $OUTFILE | |||
<li>Generate random 8bp barcodes and make sure they are not too similar to one another (I use >=2bp of difference) </li> | |||
python barcodegenerator.py -n $nbarcodes -l $blength | |||
<li>Select two oligo sequences and link them together with barcodes</li> | |||
[[Media:Makecprobearm.sh.txt|bash script]] | |||
<li>BLAST sequences against genome and transcriptome</li> | |||
</ol> | |||
==Designing Transcript Matching Region== | ==Designing Transcript Matching Region== | ||
Revision as of 16:56, 21 April 2016
C Probe Design
- 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.
Designing Capture Arm
The capture arm design aims to create 2 sequences that are orthogonal to the human genome/transcriptome, so that these sequences will not have non-specific binding. The two flanking sequences (one for padlock and one for latch) will have an 8bp region between them that will be a barcode unique to each C-probe.
Each C probe of a pair will have:
- Padlock matching region, this will be a different sequence for each of the pair of C probes, but will be the same across pairs
- Latch matching region, this will be a different sequence for each of the pair of C probes, but will be the same across pairs
- 8 bp barcode region, this will be unique for every C probe
- Randomly generate a large list of 18bp sequences python primergenerator.py -a -n $nprimers -l $plength -o primers.fa
- BLAST sequences against genome blastn -db hg19.fa -query primers.fa -num_descriptions 5 -task blastn-short -out primers.blastout.txt -outfmt 7
- Reorder the sequences based on BLAST results (oligos with fewest hits are best) python primerblast.analyzer.py -n $nprimers -l $plength -o $OUTFILE
- Generate random 8bp barcodes and make sure they are not too similar to one another (I use >=2bp of difference) python barcodegenerator.py -n $nbarcodes -l $blength
- Select two oligo sequences and link them together with barcodes bash script
- BLAST sequences against genome and transcriptome
Designing Transcript Matching Region
- Taken from Frei et al paper
- README for probe design software
- Must install R, Primer3, and Blast+ before using
- Make blast databases
- make one for repetitive regions (from Repbase) makeblastdb -in repbase.fa -dbtype nucl
- Make one for transcriptome; use only high quality bases (from NCBI ftp) makeblastdb -in rna_human_high_qual.fa -parse_seqids -dbtype nucl
- In R: library(devtools) install_github("nolanlab/PLAYRDesign") library(PLAYRDesign)
- Next make a library from UCSC in R: library(GenomicFeatures) txdb <- makeTxDbFromUCSC(genome = "hg19", tablename = "refGene") saveDb(txdb, "PUT THE PATH TO THE OUTPUT FILE HERE")
- Finally, run PLAYRDesign
- Select a file (any file) in the working directory you wish to use
- You must have a configuartion file in the directory you wish to use (see below)
- You can hit escape any time in the R window to stop PLAYRDesign
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.