Matt:LabNotes/2014-11-19: Difference between revisions

From ZhangLabWiki
Jump to navigation Jump to search
>Mzcai
>Mzcai
Line 58: Line 58:
*390 genes
*390 genes


[[]]
[[File:20141119_Constitutive8reads_ProbesperGene.JPG]]


*Remove probes belonging to genes with <8 probes per gene
*Remove probes belonging to genes with <8 probes per gene

Revision as of 19:06, 19 November 2014

Design FISSEQ Padlock Probeset for in situ cDNA Capture (newGenes, constitutive exons & contigs, variable capture arm length)

  • Previous designs (Matt:LabNotes/2014-11-1 & Matt:LabNotes/2014-11-10) had exactly 25bp capture arms due to some different code in ppDesigner
    • In jobFiles, H1_plus_H2_Len was the exact length the two arms had to equal
    • Changed ppDesigner so H1_plus_H2_Len is the maximum length the two arms can equal
      • Each arm can now vary between 20-25bp as was originally intended

Changes to ppDesigner Code

get_probes.pl

  • Lines 563 and 677
    • Commented out: next if(! (length($a_oligo)+length($b_oligo) == $H1_plus_H2_Len) );
    • Added: next if( length($a_oligo)+length($b_oligo) > $H1_plus_H2_Len );

output_text.pl

  • Error of dividing by 0 caused by $searchEnd being <= $searchStart (due to sortTargetFilesByChr.pl removing 25bp from each end even if total target < 50bp)
  • Inserted in lines 26-28
 if($searchEnd <= $searchStart) {
      $searchEnd = $searchStart + 1;
 }

Constitutive Exon Targets

 run ppDesigner.sh
 cat outputFile_chr*.txt > outputFile_1gap.txt
 perl ConvertToZeroGapProbe.pl
  • 7,027 probes successfully designed
  • 5,845 exons
  • 407 genes
 perl Probes2fasta.pl < outputFile_0gap.txt > outputFile_0gap.fa
 /home/kunzhang/softwares/Novocraft/novocraft/novoalign -d /home/kunzhang/RNAseq/Data/CommonFiles/refMrna.ndx -f outputFile_0gap.fa -F FA -r ALL > outputFile_0gap_novoalign_refMrna.out &
 /home/kunzhang/softwares/Novocraft/novocraft/novoalign -d /home/kunzhang/softwares/Novocraft/novocraft/human_g1k_v37 -f outputFile_0gap.fa -F FA -r ALL > outputFile_0gap_novoalign_human_g1k_v37.out &
 perl CleanupProbelist.pl
  • 5,629 probes
  • 4,735 exons
  • 406 genes
 perl  ConvertProbesHg38_2_BedHg19.pl
 coverageBed -abam BrainBulkTissueRNA_mappedhg19.sorted.bam -b outputFile_0gap_hg19.bed > outputFile_BrainBulkTissueRNA.cov

File:20141119 Constitutive RankOrderProbesOverlapReads.JPG

 perl FilterProbelist_LowReadCounts.pl > outputFile_0gap_10reads.txt
 perl RevComp_Probes.pl
  • 3,126 probes
  • 2,655 exons
  • 378 genes
 perl FilterProbelist_LowReadCounts.pl > outputFile_0gap_8reads.txt
 perl RevComp_Probes.pl
  • 3,403 probes
  • 2,879 exons
  • 390 genes

File:20141119 Constitutive8reads ProbesperGene.JPG

  • Remove probes belonging to genes with <8 probes per gene
    • Outputs: outputFile_0gap_8reads_8pergene.txt & outputFile_0gap_8reads_8pergene_RevComp.txt
 perl FilterProbelist_LowProbesPerGene.pl
 mv outputFile* newGenes_variableArms_ppDesignerFiles

Contiguous Exon Targets

 run ppDesigner.sh
 cat outputFile_chr*.txt > outputFile_1gap.txt
 perl ConvertToZeroGapProbe.pl
  • 10,429 probes successfully designed
  • 6,435 exons
  • 406 genes
 perl Probes2fasta.pl < outputFile_0gap.txt > outputFile_0gap.fa
 /home/kunzhang/softwares/Novocraft/novocraft/novoalign -d /home/kunzhang/RNAseq/Data/CommonFiles/refMrna.ndx -f outputFile_0gap.fa -F FA -r ALL > outputFile_0gap_novoalign_refMrna.out &
 /home/kunzhang/softwares/Novocraft/novocraft/novoalign -d /home/kunzhang/softwares/Novocraft/novocraft/human_g1k_v37 -f outputFile_0gap.fa -F FA -r ALL > outputFile_0gap_novoalign_human_g1k_v37.out &
 perl CleanupProbelist.pl
  • 7,390 probes
  • 4,969 exons
  • 406 genes
 perl  ConvertProbesHg38_2_BedHg19.pl
 coverageBed -abam BrainBulkTissueRNA_mappedhg19.sorted.bam -b outputFile_0gap_hg19.bed > outputFile_BrainBulkTissueRNA.cov

File:20141119 Contig RankOrderProbesOverlapReads.JPG

 perl FilterProbelist_LowReadCounts.pl > outputFile_0gap_10reads.txt
 perl RevComp_Probes.pl
  • 4,250 probes
  • 2,830 exons
  • 390 genes
 perl FilterProbelist_LowReadCounts.pl > outputFile_0gap_8reads.txt
 perl RevComp_Probes.pl
  • 4,598 probes
  • 3,073 exons
  • 394 genes
  • Remove probes belonging to genes with <8 probes per gene
    • Outputs: outputFile_0gap_8reads_8pergene.txt & outputFile_0gap_8reads_8pergene_RevComp.txt
 perl FilterProbelist_LowProbesPerGene.pl
 mv outputFile* newGenes_variableArms_contigs_ppDesignerFiles