Friday, May 29, 2015

5 29 2015 Designing Primers for Sequencing Target areas

Today after some conversations between Brent and Steven based on the previous primer checks (here and here) the idea arose that not all populations may share the same isoforms of the target genes. To determine if the genes are the same in all populations, its necessary to produce a sequence containing the primer region plus 100 bp on either side for sequencing in sub samples from each population. This put me in a unique situation. I needed to not only have the primer for the target gene, but also a primer that encompasses the target gene + original primer + 100 bp on either side of the sequence. I decided to use the NCBI Primer Blast, that I used to generate the original primers, to create the new primers.

First I need to figure out where in the sequence does the original primer amplify. To do this I have to upload the original sequence and the primer sequences into NCBI. For this I'm looking at the TLR2.1 primer since it has the most interesting results from the primer check.

The original primer sequences for TLR2.1:
FWD ACAAAGATTCCACCCGGCAA
REV  ACACCAACGACAGGAAGTGG
Product Size  109 bp

The original sequence for target gene:

tcaaaacagattttggagtgtaacgtctttaaaaattactaatacagcaacaatgaaatcattttgccatgaattgcttcacgatttgccagaacagcttctggccatttctattcattgtccactctacatgagatgctgacttgagaagcaccttcagagattgattaatatgttcagaacgaatatcttgtaatatgatcagcaatagcttgtcaccgcctccgtcagccaaagtgtggttggcaatggcggcttcgtacttacaccattgatcgtccaagaaattattggacagcaccaaaatgaattttctgctgacttcaatgctttccagaaagtcatccacaaagattccacccggcaagatatctctctcgtgcaaacaaagcttgtacttgtttctctctaccatctgaacaagttcagacattatccacttcctgtcgttggtgttgtaagcaacaaaaccatcatagaggaattcatcatctgtaaatcttttataccccatgggttttctgtttttgcatgtgtatatataatacttaatgtcccaacgaaactttctcaacagaacaatggtcaaaatcaccaaggacactaatgcaaacgaagcgcatactagaatgacataggggcttattttgtgacattcggacgcaggatcaaagtgtgttatactttttcctttcaagtctgacggggaagcacagatgtattgatgaggatatccaacagtcctattttggttgttttcaacccagagtcggaaccattccagttggcatccacagtccaatggattcaaggacacatcaattcgaaaatcttttttcttccaaaggaaagcaggcagagaagattcattgatactgctcaaacgattgccacgaaaaatcaattcattcaatctagtcaaattttgaacaccgctctttcgaatttctgttattcggttaacactgtagtcaacggacgttagctgtgaaaaatcagtcagtgtagtcatttgatcgatttcgccattaacaatagccaaaactgacagatcctttagaggtgatagcatctctgcaacatcagtgttggatagatcggtgttcattattttcaatgcttttatatttctacattcagaaaagatgtatttgaagcctcctacgcttttttccaccatacgtccgatagagagtgttcggagggaattggatgataatgacctaccatctaaagcgactccatgcaagttatcaatgatcaaagtttgcaaacgattcatattcaaaaaggaatacaaggaaatgtcatgaattgatgtgtgtgaaacattaaaaaattcaagggatttcaggcagtagccctcttgtatgaaatttgttgtttcccgaataccattcactgccaggtgcaattgttgcaagtttgggaaggcggctatacgtttctcattggataaacaaaactcgggaatcttctcaaaagaattgttagtaaaatacagttcttttaaagaaggaactgagctatttggcaatttataccttgagattaggttgcttcctaaatcgattttgtggatattgattaatttagcgaaggatttaaaag


First I run the sequence and the primers through NCBI Primer Blast:

Once the Primer Blast returns the primer target. I can then find where the primer begins and ends in the original sequence.

The forward primer begins at roughly 340 bp and the reverse primer begins at roughly 460 bp. Using these numbers I subtract 100 bp from the start of the forward primer and add 100 bp to the beginning of the reverse primer. This extends the range to 240 bp - 560 bp which encompasses the original target and primer area as well as 100 bp on either side. To give the Primer Blast some leeway to produce the new primers, I give a 100 bp range for the forward and reverse primers. The range for the forward primer becomes 140-240 bp. The range for the reverse primer becomes 560-660 bp. I also required the product size to be between 325 bp and 1000 bp long. This is so that the 100 bp on either side plus the 109 bp original target are amplified.

The primers that NCBI now covers the entire area that needs to be sequenced. All I have to do is choose primers that have low self complimentarity.


After looking at the offered primers I chose one that should give a 425 bp product with low self complimentarity.

This is the quickest way I've found to take the primers I produced originally and create primers to sequence the target area of interest. If this works I'll do this with all the primers of interest that worked in the previous primer runs.

Thursday, May 28, 2015

5 28 2015 New Primer Check (13-23)

Today I completed the primer checks for the new primers. This run is less promising than yesterdays as it doesn't have any primers that show any real difference in them. There were more than a couple that failed in some manner.

To make the working stock of the primers.

  1. 90 ul of Nuclease free (or Nanopure) H20 to 1.5 ml tube
  2. 10 ul of Stock primers
  3. Vortex briefly
I made these carefully in order to produce the following 12 pairs of primers. 

1626HSP70c_FWDAGGAAAGGTCGGGAGAGGAAJH5/21/20152055O.luridaHeat shock 70 kDa protein 12A
1625HSP70c_REVACCTCGGACTTTGGACGAACJH5/21/20152055O.luridaHeat shock 70 kDa protein 12A
1624p29ING4_FWDTACCTTTGGGCTTCACCGTCJH5/21/20152055O.luridaInhibitor of growth protein 4 (p29ING4)
1623p29ING4_REVGTCCATCACACACCCCTCAGJH5/21/20152055O.luridaInhibitor of growth protein 4 (p29ING4)
1622CerS2_FWDTTGTCGGTCTCCTCCTGCTAJH5/21/20152055O.luridaCeramide synthase 2 (CerS2) (LAG1 longevity assurance homolog 2)
1621CerS2_REVCCGTCTTCTGAGCCATCGTTJH5/21/20152055O.luridaCeramide synthase 2 (CerS2) (LAG1 longevity assurance homolog 2)
1620GABABR1_FWDCCGAGGAGGACACGAAACTCJH5/21/20152055O.luridaGamma-aminobutyric acid type B receptor subunit 1 (GABA-B receptor 1) (GABA-B-R1) (GABA-BR1) (GABABR1) (Gb1)
1619GABABR1_REVCGGACAGGTTCTGGATTCCGJH5/21/20152055O.luridaGamma-aminobutyric acid type B receptor subunit 1 (GABA-B receptor 1) (GABA-B-R1) (GABA-BR1) (GABABR1) (Gb1)
1618HSP70d_FWDTTTGTCTCACCGGCTTTGTGJH5/21/20152055O.luridaHeat shock 70 kDa protein 6 (Heat shock 70 kDa protein B')
1617HSP70d_REVGACATGAGACCAAAGACGCCJH5/21/20152055O.luridaHeat shock 70 kDa protein 6 (Heat shock 70 kDa protein B')
1616THRa_FWDGACACTATCCTCACTCGGCGJH5/21/20152055O.luridaThyroid hormone receptor alpha (Nuclear receptor subfamily 1 group A member 1)
1615THRa_REVGGGTGCCGAGTAAACAAGGAJH5/21/20152055O.luridaThyroid hormone receptor alpha (Nuclear receptor subfamily 1 group A member 1)
1614Defensin_FWDTCTAGCGGAGTTTGTTGGGGJH5/21/20152055O.luridaBig defensin
1613Defensin_REVATGGCTGTCGGAGGAGGATTJH5/21/20152055O.luridaBig defensin
1612GRB2_FWDAACTTTGTCCACCCAGACGGJH5/21/20152055O.luridaGrowth factor receptor-bound protein 2 (Adapter protein GRB2) (Protein Ash) (SH2/SH3 adapter GRB2)
1611GRB2_REVCCAGTTGCAGTCCACTTCCTJH5/21/20152055O.luridaGrowth factor receptor-bound protein 2 (Adapter protein GRB2) (Protein Ash) (SH2/SH3 adapter GRB2)
1610H3.3_FWDCACGCTCTCCTCGAATCCTCJH5/21/20152055O.luridaHistone H3.3
1609H3.3_REVAAGTTGCCTTTCCAGCGTCTJH5/21/20152055O.luridaHistone H3.3
1608H2A.V_FWDTGCTTTCTGTGTGCCCTTCTJH5/21/20152055O.luridaHistone H2A.V (H2A.F/Z) (Fragment)
1607H2A.V_REVTATCACACCCCGTCACTTGCJH5/21/20152055O.luridaHistone H2A.V (H2A.F/Z) (Fragment)
1606H2A_FWDGCTGGGGTTTTTCTGGGTCTJH5/21/20152055O.luridaHistone H2A
1605H2A_REVGGAACTACGCCGAGAGAGTGJH5/21/20152055O.luridaHistone H2A
After producing the working stocks I then made a 10 reaction master mix. 

Master Mix reagent table
VolumeReactions X10
Ssofast Evagreen MM10100
FWD Primer0.55
REV Primer0.55
Nanopure H2O880
cDNA1
Protocol:
  1. Added Ssofast to each of 12 tubes
  2. Added Nanopure water to each of 12 tubes
  3. Carefully added FWD primer, then Reverse primer to the appropriate tube
  4. Vortex briefly
  5. One Master Mix used per column in plate
  6. Pipette 19 ul of master mix to each well for the appropriate column
  7. Either No Template Control (NTC) or sample for each Row in the plate
  8. 1 ul Sample/NTC per well in the appropriate row
Table Layout:
HSP70cP29INGCerS2GABABRHSP70dTHRaDefensinGRB2H3.3H2A.VH2A
123456789101112
NTCNTCNTCNTCNTCNTCNTCNTCNTCNTCNTC
NT1NT1NT1NT1NT1NT1NT1NT1NT1NT1NT1
HT1HT1HT1HT1HT1HT1HT1HT1HT1HT1HT1
ST1ST1ST1ST1ST1ST1ST1ST1ST1ST1ST1
NC1NC1NC1NC1NC1NC1NC1NC1NC1NC1NC1
HC1HC1HC1HC1HC1HC1HC1HC1HC1HC1HC1
SC1SC1SC1SC1SC1SC1SC1SC1SC1SC1SC1
NTCNTCNTCNTCNTCNTCNTCNTCNTCNTCNTC

qPCR Program:
Sybr New Plate+Sybr cDNA 60 melt 2 Read
StepTemperatureTime
Initiation95 C10 min
Elongation95 C30 sec
55 C1 min
Read
72 C30 sec
Read
Repeat Elongation 39 times
Termination95 C1 min
55 C1 sec
Melt Curve Manual ramp 0.2C per sec Read 0.5 C55 - 95 C30 sec
21 C10 min
End

Results:
HSP70c Amplification
 HSP70c Melt Curve

p29ING Amplification
 p29ING Melt Curve

CerS2 Amplification
 CerS2 Melt Curve

GABABR Amplification
 GABABR Melt Curve

HSP70d Amplification
 HSP70d Melt Curve

THRa Amplification
 THRa Melt Curve

Defensin Amplification
 Defensin Melt Curve

GRB2 Amplification
 GRB2 Melt Curve

H3.3 Amplification
 H3.3 Melt Curve

H2A.V Amplification
 H2A.V Melt Curve

H2A Amplification
 H2A Melt Curve

Posted below are assumptions based on the very limited data from this check. It is only being produced to eliminate useless primers from future tests and highlight possibly interesting primers. Only with further testing will we be able to determine true trends. 

CerS2, HSP70d, THRa, and Defensin failed due to no amplification or poor melt curve conditions. 

HSP70c produced a weird bump inline with the other melt curves. This primer should probably not be used due to the inability to eliminate this issue.

H3.3, H2A.V, and H2A all showed similar results with H2A being very similar in all samples. These could be candidates for normalizing genes to use with Actin. 

p29ING also appeared to have uniform expression in all samples. This could also be another normalizing gene.

GABABR appeared to have down regulation in Heat Stress for Fidalgo and Dabob but there's no way to tell unless we run more replicates. 

GRB2 had some variable expression between heat shock and control samples but we'll need to run further tests to ensure this. 

Tomorrow I should finalize a list of the primers of interest.  Hopefully next week I'll be able to run them with a full complement of samples and normalizing genes. 

You can find the raw qPCR data here.

5 27 2015 New Primer Check (1-12)

Yesterday I received the new primers that I designed using the oly transcriptome. You can read about that here. Sam rehydrated the original stocks and I produced work stocks from them to use in a plate.

To make the working stock of the primers.

  1. 90 ul of Nuclease free (or Nanopure) H20 to 1.5 ml tube
  2. 10 ul of Stock primers
  3. Vortex briefly
I made these carefully in order to produce the following 12 pairs of primers. 


Hspb11_FWDATGTTTCCTGGTCTCCGTCAJH5/21/20152055O.luridaHeat shock protein beta-11 (Hspb11) (Placental protein 25) (PP25)
Hspb11_REVCATCAACGCCAGGGGAACTTJH5/21/20152055O.luridaHeat shock protein beta-11 (Hspb11) (Placental protein 25) (PP25)
GDF-8_FWDCCGTGGATGTCGCAGAAAGAJH5/21/20152055O.luridaGrowth/differentiation factor 8 (GDF-8) (Myostatin) (Myostatin-1) (zfMSTN-1) (Myostatin-B)
GDF-8_REVCTGCTTTCTCCGTCCCCTTTJH5/21/20152055O.luridaGrowth/differentiation factor 8 (GDF-8) (Myostatin) (Myostatin-1) (zfMSTN-1) (Myostatin-B)
HSP70b_FWDAAGTACCTTGGGGAGCTTGCJH5/21/20152055O.luridaHeat shock 70 kDa protein 12B
HSP70b_REVTCCACAGACTTTCCTCCCCAJH5/21/20152055O.luridaHeat shock 70 kDa protein 12B
GRP-78_FWDGAGAAACCACGCAGGGAGAAJH5/21/20152055O.lurida78 kDa glucose-regulated protein (GRP-78) (Heat shock 70 kDa protein 5) (Immunoglobulin heavy chain-binding protein) (BiP)
GRP-78_REVCATCAGCATCGAAGGCAACGJH5/21/20152055O.lurida78 kDa glucose-regulated protein (GRP-78) (Heat shock 70 kDa protein 5) (Immunoglobulin heavy chain-binding protein) (BiP)
CARM1_FWDTGGTTATCAACAGCCCCGACJH5/21/20152055O.luridaHistone-arginine methyltransferase CARM1 (EC 2.1.1.-) (EC 2.1.1.125) (Coactivator-associated arginine methyltransferase 1) (Protein arginine N-methyltransferase 4)
CARM1_REVGTTGTTGACCCCAGGAGGAGJH5/21/20152055O.luridaHistone-arginine methyltransferase CARM1 (EC 2.1.1.-) (EC 2.1.1.125) (Coactivator-associated arginine methyltransferase 1) (Protein arginine N-methyltransferase 4)
BMP-2_FWDTGAAGGAACGACCAAAGCCAJH5/21/20152055O.luridaBone morphogenetic protein 2 (BMP-2) (Bone morphogenetic protein 2A) (BMP-2A)
BMP-2_REVTCCGGTTGAAGAACCTCGTGJH5/21/20152055O.luridaBone morphogenetic protein 2 (BMP-2) (Bone morphogenetic protein 2A) (BMP-2A)
PGE/EP4_FWDACAGCGACGGACGATTTTCTJH5/21/20152055O.luridaProstaglandin E2 receptor EP4 subtype (PGE receptor EP4 subtype) (PGE2 receptor EP4 subtype) (Prostanoid EP4 receptor)
PGE/EP4_REVATGGCAGACGTTACCCAACAJH5/21/20152055O.luridaProstaglandin E2 receptor EP4 subtype (PGE receptor EP4 subtype) (PGE2 receptor EP4 subtype) (Prostanoid EP4 receptor)
CRAF1_FWDAGCAGGGCATCAAACTCTCCJH5/21/20152055O.luridaTNF receptor-associated factor 3 (EC 6.3.2.-) (CD40 receptor-associated factor 1) (CRAF1) (TRAFAMN)
CRAF1_REVACAAGTCGCACTGGCTACAAJH5/21/20152055O.luridaTNF receptor-associated factor 3 (EC 6.3.2.-) (CD40 receptor-associated factor 1) (CRAF1) (TRAFAMN)
NFKBina_FWDGATGGCGGTGCATGTGTTAGJH5/21/20152055O.luridaNF-kappa-B inhibitor alpha (I-kappa-B-alpha) (IkB-alpha) (IkappaBalpha) (REL-associated protein pp40)
NFKBina_REVCGAGGAGAACCTTGTGCAGTJH5/21/20152055O.luridaNF-kappa-B inhibitor alpha (I-kappa-B-alpha) (IkB-alpha) (IkappaBalpha) (REL-associated protein pp40)
PGRP-S_FWDGAGACTTCACCTCGCACCAAJH5/21/20152055O.luridaPeptidoglycan recognition protein 1 (Peptidoglycan recognition protein short) (PGRP-S)
PGRP-S_REVAACTGGTTTGCCCGACATCAJH5/21/20152055O.luridaPeptidoglycan recognition protein 1 (Peptidoglycan recognition protein short) (PGRP-S)
TLR2.1_FWDACAAAGATTCCACCCGGCAAJH5/21/20152055O.luridaToll-like receptor 2 type-1
TLR2.1_REVACACCAACGACAGGAAGTGGJH5/21/20152055O.luridaToll-like receptor 2 type-1
GDF-8b_FWDAACTGATTCTGCTCGTCGCAJH5/21/20152055O.luridaGrowth/differentiation factor 8 (GDF-8) (Myostatin)
GDF-8b_REVTGTTCTTCCACCCACCACTGJH5/21/20152055O.luridaGrowth/differentiation factor 8 (GDF-8) (Myostatin)

After producing the working stocks I then made a 10 reaction master mix. 

Master Mix reagent table
VolumeReactions X10
Ssofast Evagreen MM 10100
FWD Primer0.55
REV Primer0.55
Nanopure H2O880
cDNA1
Protocol:
  1. Added Ssofast to each of 12 tubes
  2. Added Nanopure water to each of 12 tubes
  3. Carefully added FWD primer, then Reverse primer to the appropriate tube
  4. Vortex briefly
  5. One Master Mix used per column in plate
  6. Pipette 19 ul of master mix to each well for the appropriate column
  7. Either No Template Control (NTC) or sample for each Row in the plate
  8. 1 ul Sample/NTC per well in the appropriate row
Table Layout:
Hspb11GDF-8HSP70bGRF-78CARM1BMP2PGE/EP4CRAF1NFKBinaPGRP-STLR2.1GDF-8b
123456789101112
NTCNTCNTCNTCNTCNTCNTCNTCNTCNTCNTCNTC
NT1NT1NT1NT1NT1NT1NT1NT1NT1NT1NT1NT1
HT1HT1HT1HT1HT1HT1HT1HT1HT1HT1HT1HT1
ST1ST1ST1ST1ST1ST1ST1ST1ST1ST1ST1ST1
NC1NC1NC1NC1NC1NC1NC1NC1NC1NC1NC1NC1
HC1HC1HC1HC1HC1HC1HC1HC1HC1HC1HC1HC1
SC1SC1SC1SC1SC1SC1SC1SC1SC1SC1SC1SC1
NTCNTCNTCNTCNTCNTCNTCNTCNTCNTCNTCNTC
qPCR Program:
Sybr New Plate+Sybr cDNA 60 melt 2 Read
StepTemperatureTime
Initiation95 C10 min
Elongation95 C30 sec
55 C1 min
Read
72 C30 sec
Read
Repeat Elongation 39 times
Termination95 C1 min
55 C1 sec
Melt Curve Manual ramp 0.2C per sec Read 0.5 C55 - 95 C30 sec
21 C10 min
End

Results:
HSPb11 Amplification

HSPb11 Melt Curve

GDF-8 Amplification
 GDF-8 Melt Curve
HSP70b Amplification

 HSP70b Melt Curve
GRF-78 Amplification


GRF-78 Melt Curve
CARM1 Amplification

 CARM1 Melt Curve
BMP2 Amplification


BMP2 Melt Curve

PGE/EP4 Amplification

PGE/EP4 Melt Curve
CRAF1 Amplification


CRAF1 Melt Curve
NFKBina Amplification

 NFKBina Melt Curve

PGRP-S Amplification
 PGRP-S Melt Curve
TLR2.1 Amplification

 TLR2.1 Melt Curve

GDF-8b Amplification
 GDF-8b Melt Curve


There was some very interesting results. 

CARM1, BMP2, HSP11b, and PGE/EP4 All showed down regulation from heat shock in all populations.

GDF-8 Had High expression in Fidalgo

CRAF1 had Up regulation in the Oyster Bay population and down regulation in others

PGRP-S had Up regulation in the Dabob population and down regulation in the others

TLR2.1 Had the most interesting response. It was down regulated in the Fidalgo bay population, Equal Expression in the Oyster Bay population, and did not amplify in either Dabob sample. 

HSP70b, GRF-78, NFKBina and GDF-8b had severe issues with the melt curve and odd expression in the samples. These four primers are probably not useful. 

I'm testing the next batch of primers now and will report the results when the come in. 

You can see the raw qPCR file here.