Name: merging absolute and relative quantitative pcr data to quantify stat3 splice variant transcripts
Uploaded: 2016-10-09T13:00:00
Duration: 11 min 19 s
Description: Watch this Scientific Journal Video about Merging Absolute and Relative Quantitative PCR Data to Quantify STAT3 Splice Variant Transcripts at JoVE.com

The authors would like to acknowledge the National Institutes of Health-NHLBI for the Program Project Grant on the Role of Eosinophils in Airway Inflammation and Remodeling: P01HL088584 (PI: N. Jarjour), and the University of Wisconsin Carbone Cancer Center and Department of Medicine for intramural funding. We thank Douglas Annis for cloning the four STAT3 variants.

NOTE: Peripheral blood eosinophils were received without identifying information in accord with a protocol approved (#2013-1570) by the University of Wisconsin-Madison Center for Health Sciences Institutional Review Board. Signed informed consent from the donor was obtained for the use of each sample in research.
1. Creating Plasmids as Template Standards
<ol>
	<li>Create primers for an amplicon spanning both STAT3 splice sites, with 5&#39;- KpnI and NheI restriction sites (and 5&#39;-extensions) as per Table 1a. 
	NOTE: KpnI and NheI sites were chosen to enable inserts to be ligated into a modified pET-28a vector with kanamycin resistance25,26. KpnI site had been added into the multiple cloning site.</li>
	<li>Using freshly donated eosinophils, prepare duplicate samples of 2.5 x 106 cells/ml in cell culture medium (Roswell Park Memorial Institute (RPMI)1640 medium). Incubate for 1 hr at 37 &#176;C under 5% CO2 and 10% humidity.
	<ol>
		<li>Prepare complementary (c)DNA from samples (at least 1 x 106 cells per preparation) as per manufacturer&#39;s instructions using the RNA extraction and cDNA synthesis kits.</li>
	</ol>
	</li>
	<li>From template cDNA prepared in step 1.2.1, amplify STAT3 splice variants using amplification PCR as per Table 2a.</li>
	<li>Resolve amplicons in 2% agarose Tris-acetate-ethylenediaminetetraacetate (TAE) gel27. Excise bands from the gel and purify according to gel excision kit instructions.</li>
	<li>Cut amplicons and plasmid with appropriate restriction enzymes (overview of restriction in reference27), using volumes, incubation times and temperatures recommended by enzyme provider; and purify as in step 1.4. Ligate restricted amplicons into plasmid under provider-recommended conditions.</li>
	<li>Prepare a negative control (restricted plasmid DNA without insert but with ligase) and a positive control (unrestricted plasmid with kanamycin resistance). Perform standard plasmid transformations of competent E. coli DH5&#945;28 using ligation mixtures27.</li>
	<li>Incubate transformed E. coli in 1 ml Luria-Broth (LB) medium (prepared as instructed27) shaking for 1 hr at 37 &#176;C. Spread transformants on LB plates containing 50 &#181;g/ml kanamycin and incubate at 37 &#176;C overnight.</li>
	<li>Pick several colonies from the STAT3 &#945; and STAT3 &#946; plates using sterile toothpicks27. Transfer each to 2 ml LB. Grow cultures overnight at 37 &#176;C in a shaking incubator.</li>
	<li>Purify plasmids from bacterial cultures by following instructions provided in a DNA purification kit. Store plasmids at -20 or -80 &#176;C.</li>
	<li>Sequence plasmids from several colonies by preparing 20 &#181;l sequencing reactions. Pipette 3 &#181;l sequencing reaction buffer, 2 &#181;l sequencing reaction mix, 12 &#181;l ultrapure water, 1 &#181;l plasmid as template and 2 &#181;l sequencing primer. 
	NOTE: If using pET-28a vector, use sequencing primer 5&#39;-TAA TAC GAC TCA CTA TAG GGG-3&#39;.
	<ol>
		<li>Assess electrophoretograms of plasmids encoding each variant: S&#945;, S&#946;, &#916;S&#945; and &#916;S&#946;25.</li>
	</ol>
	</li>
	<li>Measure concentration of pure plasmid DNA in ng/&#181;l. If reading absorbance at 260 nm using a spectrophotometer, calculate DNA concentration using the Beer-Lambert law29: 
	<img alt="Equation 1" src="/files/ftp_upload/54473/54473eq1.jpg" /> 
	Where c = concentration, A = absorbance, &#949; (extinction coefficient) = 0.02 (&#181;g/ml)-1 cm-1 for double-stranded DNA and L = path length of light (usually 1 cm).</li>
	<li>Using the molecular weights of STAT3 splice variant cDNA amplicons and the vector, calculate the copy number per &#181;l. 
	NOTE: Molecular weight per base pair (bp) is estimated as 650 Da. 
	<img alt="Equation 2" src="/files/ftp_upload/54473/54473eq2.jpg" /></li>
</ol>
2. Analyzing Primer Specificity for Absolute qPCR
<ol>
	<li>Prepare 1 in 10 dilution series of STAT3 plasmids, with ~108 to 103 copies per &#181;l (20 &#181;l) using ultrapure water. Measure concentration of most concentrated (~108 copies per &#181;l, see step 1.11).</li>
	<li>Use diluted plasmids to prepare four &#34;non-target&#34; mixes (i.e., STAT3&#916;S&#946; negative control containing equal concentrations of STAT3 S&#945;, &#916;S&#945;, S&#946; but no &#916;S&#946;) at 106 copies per &#181;l each. Prepare a &#34;target&#34; mix with equal concentrations of all four templates each at 106 copies per &#181;l.</li>
	<li>Prepare primer pair solutions for each splice variant (see Table 1b and Figure 1) by making ~60 &#181;l of primers each at a 7 &#181;M concentration (final concentration in sample will be 560 nM).</li>
	<li>Dilute DNA polymerase/dsDNA-binding dye mix 7:5 with pure H2O.</li>
	<li>Set up Assay 1 in a 96-well PCR plate as shown in template (Table 3).
	<ol>
		<li>Add 21 &#181;l diluted DNA polymerase/dsDNA-binding dye mix, then 2 &#181;l of primer mix and 2 &#181;l of template (as per Table 2b). Instead of template, add 2 &#181;l filtered H2O to the no template control (NTC) well. Set up reactions in duplicate to assess repeatability30.</li>
	</ol>
	</li>
	<li>Seal qPCR plate with adhesive cover and centrifuge for 5 min at 1,200 x g at 12 &#176;C.</li>
	<li>If using software listed on Materials, set up experiment as described.
	<ol>
		<li>Turn on qPCR machine and insert sealed qPCR plate. Click File &#62; New &#62; Next &#62; Select &#34;New detector&#34;, choose reporter, quencher and provide name. Set up a &#34;new detector&#34; for each splice variant.</li>
		<li>Select Next and set up standards as prompted on the Set Up Sample Plate page, ensuring quantities are entered into the table and appropriate detectors are selected. Click finish.</li>
		<li>Set up cycling as per Table 2b. Ensure that fluorescence reading (to determine Ct) occurs during the 72 &#176;C step of each cycle. Select Run.</li>
		<li>When run is complete, click the Results tab &#62; Amplification plot tab. Evaluate output amplification plot to assess data quality (look for an exponential plot, as in Figure 2). 
		NOTE: Ensure amplification plots are mostly exponential and that the cycle threshold can be set at the exponential part of the plot. Ensure NTCs are not amplified and other standard points demonstrate a low Ct value with a high &#916;Rn.</li>
		<li>To export results to spreadsheet software, click File &#62; Export &#62; Results.</li>
	</ol>
	</li>
</ol>
3. Assessing Absolute qPCR Assay Specificity and Repeatability
<ol>
	<li>Repeatability
	<ol>
		<li>Use the data from Step 2.7.5 to calculate standard deviation of Ct (threshold cycle for fluorescence) values. 
		<img alt="Equation 3" src="/files/ftp_upload/54473/54473eq3.jpg" /> 
		Where N = number of samples (2 if done in duplicate), <img alt="Equation 4" src="/files/ftp_upload/54473/54473eq4.jpg" /> = sample&#39;s Ct and <img alt="Equation 5" src="/files/ftp_upload/54473/54473eq5.jpg" /> = sample Ct mean. Check that the standard deviation of Ct values of duplicates is &#8804;0.2. Check that Ct values in all but NTC wells are smaller than 38.</li>
		<li>If repeatability is poor optimize cycling parameters by either increasing or decreasing annealing time.</li>
	</ol>
	</li>
	<li>Plot log copy number (as determined in Step 1.12 and prepared in Step 2.1) vs Ct value to create a standard curve; yielding the following equation: 
	<img alt="Equation 6" src="/files/ftp_upload/54473/54473eq6.jpg" /> 
	Where y is Ct, m is the gradient, x is log(copy number) and b is the intercept value. 
	NOTE: The R2 value of the linear regression will be &#8805;0.95 for good data.</li>
	<li>Calculate amplification efficiency using standard curve and the equation: 
	<img alt="Equation 7" src="/files/ftp_upload/54473/54473eq7.jpg" /> 
	NOTE: Aim for efficiency &#8805;75%.</li>
	<li>Calculate specificity from qPCR assay 1 using Too&#39;s formula: 
	<img alt="Equation 8" src="/files/ftp_upload/54473/54473eq8.jpg" /> 
	Where &#916;CtToo = Cttarget - Ctnon-target , describing the difference in threshold cycle number for specific and non-specific amplification 
	NOTE: Specificity should exceed four orders of magnitude (i.e., specificity factor &#8805;4).
	<ol>
		<li>If specificity is poor, optimize by lowering primer concentration. Set up assays (as described in steps 2.5-2.7) with final primer concentrations ranging from 100 nM to 500 nM.</li>
	</ol>
	</li>
</ol>
4. Performing Relative qPCR Assays
<ol>
	<li>Treat Cells with Cytokines to Promote Transcription.
	<ol>
		<li>For freshly donated eosinophils, prepare duplicate samples of 2.5 x 106 cells/ml in cell culture medium (RPMI 1640 medium) and treat with 50 ng/ml interleukin-3 (IL3) and/or 50 ng/ml tumor necrosis factor &#945; (TNF&#945;) for periods of 3, 6, 9 and 20 hr at 37 &#176;C under 5% CO2 and 10% humidity.</li>
	</ol>
	</li>
	<li>Prepare cDNA from eosinophil samples (at least 1 x 106 cells per preparation) as per manufacturer&#39;s instructions using the RNA extraction and cDNA synthesis kits.</li>
	<li>Prepare serial dilutions of two sample cDNA aliquots (control or resting samples), with a dilution range from &#34;1&#34; to &#34;1 in 1024&#34; dilution.
	<ol>
		<li>For the 1 in 4 dilution, pipette 1 &#181;l cDNA and 3 &#181;l ultrapure water.</li>
		<li>For the 1 in 16 dilution, pipette 1 &#181;l of the 1 in 4 dilution and 3 &#181;l ultrapure water, etc.</li>
	</ol>
	</li>
	<li>Set up PCR in 96-well plate with 25 &#181;l reactions as described steps 2.3-2.5 but with several primer concentrations for pan-STAT3 and GUSB (see template in Table 4).</li>
	<li>Add &#34;new detector&#34; for GUSB and follow steps 2.6-2.7, using the cycling parameters described in Table 2c.</li>
	<li>Generate standard curves (see steps 3.2 and 3.3) to determine which primer concentrations result in 95-100% amplification efficiency.</li>
	<li>Set up plate with prepared cDNA samples (from Step 4.2) and optimized primer concentrations (see Table 2c for cycling parameters, reagents and Table 5 96-well plate template).</li>
	<li>Perform steps 2.6-2.7. Repeat assay at least once and check data quality.</li>
	<li>Calculate the average Ct value of duplicate sample Cts for both STAT3 and housekeeping gene GUSB.</li>
	<li>Obtain &#916;Ct values for each sample. 
	<img alt="Equation 9" src="/files/ftp_upload/54473/54473eq9.jpg" /></li>
	<li>Designate resting sample (usually has lowest concentration of transcript of interest) as Sample 0. Calculate &#916;&#916;Ct31 for each sample (here designated as sample X): 
	<img alt="Equation 10" src="/files/ftp_upload/54473/54473eq10.jpg" /></li>
	<li>Calculate fold-increase for each sample: 
	<img alt="Equation 11" src="/files/ftp_upload/54473/54473eq11.jpg" /></li>
	<li>Reproducibility
	<ol>
		<li>Calculate the coefficient of variation (CV) of copy number for pan-STAT3 and GUSB. 
		<img alt="Equation 12" src="/files/ftp_upload/54473/54473eq12.jpg" /> 
		Where N = number of samples, = sample&#39;s calculated copy number and = sample mean. Check that CV of each sample (multiple assays) is &#8804;10%.</li>
	</ol>
	</li>
</ol>
5. Analyzing Absolute qPCR Data for Unknown Samples
<ol>
	<li>Compare Ct values obtained in relative qPCR (as described in step 4.9) for housekeeping gene GUSB to ensure samples&#39; cDNA concentrations are within one order of magnitude.
	<ol>
		<li>Dilute cDNA accordingly (e.g., if sample 1 has a Ct value of ~17.5, and other samples&#39; Ct values are ~21, sample 1 is roughly 2(21-17.5) = 11 times more concentrated. Perform a 1:1 dilution with ultrapure water within the same range without diluting more than necessary). 
		NOTE: Vastly different starting concentrations will bias the linear regression analysis (step 4.12).</li>
	</ol>
	</li>
	<li>Set up absolute qPCR assay of samples. Prepare assay 2a template plate for S&#945; and S&#946; as per Table 6; and prepare assay 2b for &#916;S&#945; and &#916;S&#946; as per Table 7. Carry out assays as described in steps 2.6-2.7.3 (and Table 2b). Repeat assays at least once.</li>
	<li>Check data quality and export as described in steps 2.7.4-2.7.5.</li>
	<li>Set up absolute qPCR 96 well plate with 25 &#181;l reactions using pan-STAT3 primers (primers in Table 1c, template in Table 8) at 400 &#181;M and a mix of all four plasmids or a single plasmid at listed total concentrations. 
	NOTE: Efficiency does not matter for absolute qPCR, but optimizing efficiency of these primers is important for relative qPCR (Step 4).</li>
	<li>Seal qPCR plate with adhesive cover and centrifuge for 5 min at 1,200 x g at 12 &#176;C.</li>
	<li>Set up &#34;new detector&#34; for pan-STAT3 as in steps 2.7.1-2.7.3.</li>
	<li>Set up cycling for pan-STAT3 as per Table 2c (same conditions as relative qPCR). Ensure that fluorescence reading (to determine Ct) occurs during the 72 &#176;C step of each cycle. Repeat assay at least once to ensure reproducibility.</li>
	<li>Check data quality and export (see steps 2.7.4-2.7.5).
	<ol>
		<li>If amplification plots are not exponential (see Figure 2), dilute sample cDNA (1:10) and repeat assay as described (step 5.7).</li>
	</ol>
	</li>
	<li>Using the equation of the standard curve (see 3.2, Figure 3) and the Ct values of samples, calculate the copy number of each splice variant and of pan-STAT3 in each sample. 
	<img alt="Equation 13" src="/files/ftp_upload/54473/54473eq13.jpg" /></li>
	<li>Plot log (absolute qPCR copy number values obtained for pan-STAT3) vs log (cumulative absolute qPCR copy number values for STAT3 splice variants). 
	NOTE: Ideal linear regression and slope values would both be ~1.</li>
	<li>Calculate repeatability (Step 3.1) and reproducibility (Step 4.13).</li>
</ol>
6. Merging Absolute and Relative qPCR Data
<ol>
	<li>Multiply the proportion of variant with the total STAT3 fold-increase to calculate fold-increase of each splice variant.</li>
	<li>Calculate the standard deviations (see Step 3.1.1) of the absolute and relative values to account for propagation of error. Determine standard error of measurement (SEM) as shown: 
	<img alt="Equation 14" src="/files/ftp_upload/54473/54473eq14.jpg" /></li>
</ol>

<ol>
	<li>Hiller, M., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=Phylogenetically+widespread+alternative+splicing+at+unusual+GYNGYN+donors.">Phylogenetically widespread alternative splicing at unusual GYNGYN donors.</a> Genome Biol. 7 (7), R65 (2006).</li><li>Treacy, M. N., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=Twin+of+I-POU:+a+two+amino+acid+difference+in+the+I-POU+homeodomain+distinguishes+an+activator+from+an+inhibitor+of+transcription.">Twin of I-POU: a two amino acid difference in the I-POU homeodomain distinguishes an activator from an inhibitor of transcription.</a> Cell. 68 (3), 491-505 (1992).</li><li>Schindler, S., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=Alternative+splicing+at+NAGNAG+acceptors+in+Arabidopsis+thaliana+SR+and+SR-related+protein-coding+genes.">Alternative splicing at NAGNAG acceptors in Arabidopsis thaliana SR and SR-related protein-coding genes.</a> BMC Genomics. 9, 159 (2008).</li><li>Szafranski, K., Kramer, M. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=It's+a+bit+over,+is+that+ok?+The+subtle+surplus+from+tandem+alternative+splicing.">It's a bit over, is that ok? The subtle surplus from tandem alternative splicing.</a> RNA Biol. 12 (2), 115-122 (2015).</li><li>Wang, M., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=Alternative+splicing+at+GYNNGY+5'+splice+sites:+more+noise,+less+regulation.">Alternative splicing at GYNNGY 5' splice sites: more noise, less regulation.</a> Nucleic Acids Res. 42 (22), 13969-13980 (2014).</li><li>Hiller, M., Platzer, M. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=Widespread+and+subtle:+alternative+splicing+at+short-distance+tandem+sites.">Widespread and subtle: alternative splicing at short-distance tandem sites.</a> Trends Genet. 24 (5), 246-255 (2008).</li><li>Tsai, K. W., Lin, W. C. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=Quantitative+analysis+of+wobble+splicing+indicates+that+it+is+not+tissue+specific.">Quantitative analysis of wobble splicing indicates that it is not tissue specific.</a> Genomics. 88 (6), 855-864 (2006).</li><li>Tadokoro, K., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=Frequent+occurrence+of+protein+isoforms+with+or+without+a+single+amino+acid+residue+by+subtle+alternative+splicing:+the+case+of+Gln+in+DRPLA+affects+subcellular+localization+of+the+products.">Frequent occurrence of protein isoforms with or without a single amino acid residue by subtle alternative splicing: the case of Gln in DRPLA affects subcellular localization of the products.</a> J Hum Genet. 50 (8), 382-394 (2005).</li><li>Bradley, R. K., Merkin, J., Lambert, N. J., Burge, C. B. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=Alternative+splicing+of+RNA+triplets+is+often+regulated+and+accelerates+proteome+evolution.">Alternative splicing of RNA triplets is often regulated and accelerates proteome evolution.</a> PLoS Biol. 10 (1), e1001229 (2012).</li><li>Zheng, C. L., Fu, X. D., Gribskov, M. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=Characteristics+and+regulatory+elements+defining+constitutive+splicing+and+different+modes+of+alternative+splicing+in+human+and+mouse.">Characteristics and regulatory elements defining constitutive splicing and different modes of alternative splicing in human and mouse.</a> RNA. 11 (12), 1777-1787 (2005).</li><li>Schaefer, T. S., Sanders, L. K., Nathans, D. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=Cooperative+transcriptional+activity+of+Jun+and+Stat3+beta,+a+short+form+of+Stat3.">Cooperative transcriptional activity of Jun and Stat3 beta, a short form of Stat3.</a> Proc Natl Acad Sci U S A. 92 (20), 9097-9101 (1995).</li><li>Waitkus, M. S., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=Signal+integration+and+gene+induction+by+a+functionally+distinct+STAT3+phosphoform.">Signal integration and gene induction by a functionally distinct STAT3 phosphoform.</a> Mol Cell Biol. 34 (10), 1800-1811 (2014).</li><li>Srivastava, J., DiGiovanni, J. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=Non-canonical+Stat3+signaling+in+cancer.">Non-canonical Stat3 signaling in cancer.</a> Mol Carcinog. , (2015).</li><li>Holland, S. M., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=STAT3+mutations+in+the+hyper-IgE+syndrome.">STAT3 mutations in the hyper-IgE syndrome.</a> N Engl J Med. 357 (16), 1608-1619 (2007).</li><li>Maritano, D., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=The+STAT3+isoforms+alpha+and+beta+have+unique+and+specific+functions.">The STAT3 isoforms alpha and beta have unique and specific functions.</a> Nat Immunol. 5 (4), 401-409 (2004).</li><li>Hevehan, D. L., Miller, W. M., Papoutsakis, E. T. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=Differential+expression+and+phosphorylation+of+distinct+STAT3+proteins+during+granulocytic+differentiation.">Differential expression and phosphorylation of distinct STAT3 proteins during granulocytic differentiation.</a> Blood. 99 (5), 1627-1637 (2002).</li><li>Yoo, J. Y., Huso, D. L., Nathans, D., Desiderio, S. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=Specific+ablation+of+Stat3beta+distorts+the+pattern+of+Stat3-responsive+gene+expression+and+impairs+recovery+from+endotoxic+shock.">Specific ablation of Stat3beta distorts the pattern of Stat3-responsive gene expression and impairs recovery from endotoxic shock.</a> Cell. 108 (3), 331-344 (2002).</li><li>Stahl, N., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=Choice+of+STATs+and+other+substrates+specified+by+modular+tyrosine-based+motifs+in+cytokine+receptors.">Choice of STATs and other substrates specified by modular tyrosine-based motifs in cytokine receptors.</a> Science. 267 (5202), 1349-1353 (1995).</li><li>Zheng, M., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=A+mix+of+S+and+DeltaS+variants+of+STAT3+enable+survival+of+activated+B-cell-like+diffuse+large+B-cell+lymphoma+cells+in+culture.">A mix of S and DeltaS variants of STAT3 enable survival of activated B-cell-like diffuse large B-cell lymphoma cells in culture.</a> Oncogenesis. 4, 184 (2016).</li><li>Livak, K. J., Schmittgen, T. D. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=Analysis+of+relative+gene+expression+data+using+real-time+quantitative+PCR+and+the+2(-Delta+Delta+C(T))+Method.">Analysis of relative gene expression data using real-time quantitative PCR and the 2(-Delta Delta C(T)) Method.</a> Methods. 25 (4), 402-408 (2001).</li><li>Zipper, H., Brunner, H., Bernhagen, J., Vitzthum, F. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=Investigations+on+DNA+intercalation+and+surface+binding+by+SYBR+Green+I,+its+structure+determination+and+methodological+implications.">Investigations on DNA intercalation and surface binding by SYBR Green I, its structure determination and methodological implications.</a> Nucleic Acids Res. 32 (12), e103 (2004).</li><li>Kelly, E. A., Liu, L. Y., Esnault, S., Quinchia Johnson, B. H., Jarjour, N. N. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=Potent+synergistic+effect+of+IL-3+and+TNF+on+matrix+metalloproteinase+9+generation+by+human+eosinophils.">Potent synergistic effect of IL-3 and TNF on matrix metalloproteinase 9 generation by human eosinophils.</a> Cytokine. 58 (2), 199-206 (2012).</li><li>Whelan, J. A., Russell, N. B., Whelan, M. A. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=A+method+for+the+absolute+quantification+of+cDNA+using+real-time+PCR.">A method for the absolute quantification of cDNA using real-time PCR.</a> J Immunol Methods. 278 (1-2), 261-269 (2003).</li><li>Too, H. P. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=Real+time+PCR+quantification+of+GFRalpha-2+alternatively+spliced+isoforms+in+murine+brain+and+peripheral+tissues.">Real time PCR quantification of GFRalpha-2 alternatively spliced isoforms in murine brain and peripheral tissues.</a> Brain Res Mol Brain Res. 114 (2), 146-153 (2003).</li><li>Turton, K. B., Annis, D. S., Rui, L., Esnault, S., Mosher, D. F. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=Ratios+of+Four+STAT3+Splice+Variants+in+Human+Eosinophils+and+Diffuse+Large+B+Cell+Lymphoma+Cells.">Ratios of Four STAT3 Splice Variants in Human Eosinophils and Diffuse Large B Cell Lymphoma Cells.</a> PloS One. 10 (5), e0127243 (2015).</li><li>Maurer, L. M., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=Extended+binding+site+on+fibronectin+for+the+functional+upstream+domain+of+protein+F1+of+Streptococcus+pyogenes.">Extended binding site on fibronectin for the functional upstream domain of protein F1 of Streptococcus pyogenes.</a> J Biol Chem. 285 (52), 41087-41099 (2010).</li><li>Ausubel, F. M., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=.">.</a> Current protocols in molecular biology. , (1987).</li><li>Grant, S. G., Jessee, J., Bloom, F. R., Hanahan, D. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=Differential+plasmid+rescue+from+transgenic+mouse+DNAs+into+Escherichia+coli+methylation-restriction+mutants.">Differential plasmid rescue from transgenic mouse DNAs into Escherichia coli methylation-restriction mutants.</a> Proc Natl Acad Sci U S A. 87 (12), 4645-4649 (1990).</li><li>Kocsis, L., Herman, P., Eke, A. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=The+modified+Beer-Lambert+law+revisited.">The modified Beer-Lambert law revisited.</a> Phys Med Biol. 51 (5), N91-N98 (2006).</li><li>Bustin, S. A., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=The+MIQE+guidelines:+minimum+information+for+publication+of+quantitative+real-time+PCR+experiments.">The MIQE guidelines: minimum information for publication of quantitative real-time PCR experiments.</a> Clin Chem. 55 (4), 611-622 (2009).</li><li>Pfaffl, M. W. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=A+new+mathematical+model+for+relative+quantification+in+real-time+RT-PCR.">A new mathematical model for relative quantification in real-time RT-PCR.</a> Nucleic Acids Res. 29 (9), e45 (2001).</li><li>Bustin, S. A. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=.">.</a> A-Z of quantitative PCR. , (2004).</li><li>Zhang, H. F., Lai, R. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=STAT3+in+cancer-friend+or+foe.">STAT3 in cancer-friend or foe.</a> Cancers (Basel). 6 (3), 1408-1440 (2014).</li><li>Caldenhoven, E., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=STAT3beta,+a+splice+variant+of+transcription+factor+STAT3,+is+a+dominant+negative+regulator+of+transcription.">STAT3beta, a splice variant of transcription factor STAT3, is a dominant negative regulator of transcription.</a> J Biol Chem. 271 (22), 13221-13227 (1996).</li><li>Zammarchi, F., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=Antitumorigenic+potential+of+STAT3+alternative+splicing+modulation.">Antitumorigenic potential of STAT3 alternative splicing modulation.</a> Proc Natl Acad Sci U S A. 108 (43), 17779-17784 (2011).</li><li>Liu, Y., Beyer, A., Aebersold, R. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=On+the+Dependency+of+Cellular+Protein+Levels+on+mRNA+Abundance.">On the Dependency of Cellular Protein Levels on mRNA Abundance.</a> Cell. 165 (3), 535-550 (2016).</li><li>Narimatsu, M., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=Tissue-specific+autoregulation+of+the+stat3+gene+and+its+role+in+interleukin-6-induced+survival+signals+in+T+cells.">Tissue-specific autoregulation of the stat3 gene and its role in interleukin-6-induced survival signals in T cells.</a> Mol Cell Biol. 21 (19), 6615-6625 (2001).</li><li>Szafranski, K., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=Physiological+state+co-regulates+thousands+of+mammalian+mRNA+splicing+events+at+tandem+splice+sites+and+alternative+exons.">Physiological state co-regulates thousands of mammalian mRNA splicing events at tandem splice sites and alternative exons.</a> Nucleic Acids Res. 42 (14), 8895-8904 (2014).</li><li>Sundin, M., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=Novel+STAT3+mutation+causing+hyper-IgE+syndrome:+studies+of+the+clinical+course+and+immunopathology.">Novel STAT3 mutation causing hyper-IgE syndrome: studies of the clinical course and immunopathology.</a> J Clin Immunol. 34 (4), 469-477 (2014).</li><li>Oltean, S., Bates, D. O. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=Hallmarks+of+alternative+splicing+in+cancer.">Hallmarks of alternative splicing in cancer.</a> Oncogene. 33 (46), 5311-5318 (2014).</li><li>Boultwood, J., Dolatshad, H., Varanasi, S. S., Yip, B. H., Pellagatti, A. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=The+role+of+splicing+factor+mutations+in+the+pathogenesis+of+the+myelodysplastic+syndromes.">The role of splicing factor mutations in the pathogenesis of the myelodysplastic syndromes.</a> Adv Biol Regul. 54, 153-161 (2014).</li><li>Tilgner, H., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=Comprehensive+transcriptome+analysis+using+synthetic+long-read+sequencing+reveals+molecular+co-association+of+distant+splicing+events.">Comprehensive transcriptome analysis using synthetic long-read sequencing reveals molecular co-association of distant splicing events.</a> Nat Biotechnol. 33 (7), 736-742 (2015).</li><li>Hiller, M., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=Widespread+occurrence+of+alternative+splicing+at+NAGNAG+acceptors+contributes+to+proteome+plasticity.">Widespread occurrence of alternative splicing at NAGNAG acceptors contributes to proteome plasticity.</a> Nat Genet. 36 (12), 1255-1257 (2004).</li><li>Rosenfeld, J. A., Malhotra, A. K., Lencz, T. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?db=PubMed&cmd=Search&doptcmdl=Citation&defaultField=Title+Word&term=Novel+multi-nucleotide+polymorphisms+in+the+human+genome+characterized+by+whole+genome+and+exome+sequencing.">Novel multi-nucleotide polymorphisms in the human genome characterized by whole genome and exome sequencing.</a> Nucleic Acids Res. 38 (18), 6102-6111 (2010).</li></ol>

The authors have nothing to disclose.

We developed this protocol in order to assess the levels and proportions of STAT3 splice variant transcripts in eosinophils and lymphoma cells and learn whether cytokine stimulation affected the levels and proportions. STAT3 is of particular interest because of its pleiotropic and uncertain functionality, with conflicting reports on whether it acts as an oncoprotein or tumor suppressor in cancer (reviewed in reference33). Differences in STAT3 &#945; and &#946; splice variant function had been characterized previously34,35, and our protocol facilitated a knock-down/re-expression analysis that suggests a need for an optimal ratio of S and &#916;S transcripts19.
Accurate quantification of distinct splice variants will facilitate further investigations of relating heterogeneous STAT3 function to splice variant composition. The protocol integrates absolute and relative qPCR data, combining the ability of absolute qPCR to calculate splice variant proportions, and relative qPCR to measure changes in overall STAT3 expression. This approach allows one to distinguish subtle differences in sequence and simultaneously measure splicing ratios at two alternative splice sites more than 50 nucleotides apart. Determining the ratios of the splicing events individually would not have yielded the remarkable finding that a co-splicing bias existed such that &#916;S&#946; levels are higher than anticipated if uses of the two sites are randomly spliced25.
Critically, absolute qPCR with the use of plasmid calibration curves enables quantification (at sub-optimal efficiency) of splice variations that result in highly similar sequences. We anticipate a novel subtle splicing qPCR assay should take roughly two months to optimize. Key steps in assay development are the creation of STAT3 plasmids used in generating standard curves for absolute qPCR; experimentally determining optimal primer sequences and cycling parameters to ensure specificity and reproducibility; and the integration of relative qPCR data derived from quantifying pan-STAT3 expression relative to GAPDH expression. The correlation of copy number quantified by pan-STAT3 versus cumulative quantification (S&#945;+&#916;S&#945;+S&#946;+&#916;S&#946;) shows that the protocol produces reliable results.
A caveat of the technique is the extensive validation process. It is necessary to assess intra-assay variability (repeatability), interassay variability (reproducibility) and specificity. The protocol outlines ways to get numeric outputs for these parameters. We deemed efficiency &#8805;75%, specificity factor &#8805;4, coefficient of variation (reproducibility) &#8804;10% and Ct standard deviation (repeatability) &#8804;0.2 as suitable thresholds30. Mutations or deletions in the STAT3 sequence amino acids 1-690 will not be discovered by this protocol, although they may influence splicing ratios. Transcript proportions might not be proportional to proteoform proportions36.
Since samples have differing starting amounts of total cDNA, absolute qPCR is suitable for comparing copy numbers of splice variants within a sample but not for inter-sample comparison unless coupled with relative qPCR using an established housekeeping gene. The method described conforms to MIQE qPCR guidelines for reproducibility30. PCR cycling parameters and primer concentrations may need to be modified to obtain reproducible data if other equipment is used. Perfect specificity is not possible without drastically compromising efficiency, but the target was amplified more efficiently by greater than four of orders of magnitude.
Linear DNA is more easily amplified than circular. If an alternate plasmid does not provide satisfactory standard curves (R2 &#60; 0.95), consider linearizing the plasmid by single site restriction prior to quantification. Optimizing qPCR is crucial for obtaining good quality data (Figure 1). Most qPCR protocols rely on two-step cycling, and machines are optimized accordingly. Non-uniform heating of the heating block may be exacerbated in three-step cycling, contributing to poor repeatability. Assays must be set up under sterile conditions with filter pipette tips and ultrapure water, ideally in a dedicated laminar flow hood. Because contaminants can lead to inconsistent results, assays should be set up under sterile conditions with filter pipette tips and ultrapure water, ideally in a dedicated laminar flow hood. For more information about qPCR optimization, refer to Bustin et al.32
Quantifying STAT3 may lead to greater insight in a number of contexts. STAT3 auto-regulates its own expression37, and the protocol described above may help to elucidate whether ratios of STAT3 splice variants contribute to regulating this positive feedback loop. The protocol could be used to study shifts in splice variant ratios as observed in cells at differing densities38 or over the course of development: it is known that the STAT3 &#945;/&#946; ratio changes at the protein level during hematopoiesis16. Sundin et al. found that an intronic single nucleotide polymorphism biased splicing of exon 12 in STAT3 of a patient with Job&#39;s syndrome39. It is conceivable that one of the many SNPs present in the introns between exon 21 and 22, or exon 22 and 23 may contribute to splicing ratios of &#916;S/S and &#945;/&#946; respectively. The assay could be used to quantify STAT3 transcripts in cancerous cells, where mutations or changes in splicing regulation may introduce bias to the splicing process40. Mutations in splicing factors (like SF3B1), as observed in myelodysplastic syndromes41 may also lead to changes that can be measured by this protocol.
More broadly, this approach specifically detects co-association in splicing, which is not feasible with conventional RNA-Seq, nor standard qPCR. While the phenomenon of mutually exclusive exon splicing demonstrates coordination of splicing decisions, the co-association of other splicing events has not been well-researched. A recently described alternative method, in which RNA-Seq was modified so as to interrogate full-length cDNA, suggests distant splicing events are more co-dependent than previously thought42.
STAT3 contains a donor tandem splice site. Acceptor tandem splice sites are more frequent43 and the principles of the outlined protocol could serve as a starting point for developing assays for coincidence detection of NAGNAG splicing and other splicing events within 200 nucleotides. Other potential applications include quantification of coincidence of other subtle sequence differences, like indels or double/triple nucleotide polymorphisms44.

Short-range (tandem) alternative splicing, where alternate acceptor or donor sites are in close proximity to one another, is common in mammals1, invertebrates2 and plants3. It is estimated that 20% of mammalian genes contain alternative splice sites 2-12 nucleotides apart4. Many of these sites are 3 nucleotides apart and result in inclusion or exclusion of a single codon. There is debate about the nature of splicing regulation at these sites5,6 with some arguing that the splicing motif differences are so subtle that selection is stochastic7, while others infer regulation based on tissue specificity8.
Tandem splice site selection has been analyzed semi-quantitatively using modified capillary electrophoresis7, and high-resolution gel electrophoresis8. RNA-Seq (RNA sequencing) reads can be used to quantify the splicing ratios at each splice site. In this way, RNA-Seq data has provided insight into the regulation of tandem splice sites9. It has also enabled prediction of expected splice variant ratios based on nucleotide motif10. Most of the emphasis on splicing that includes or excludes a single codon has been on the more commonly occurring tandem acceptor splice sites, known as NAGNAGs (where N = any nucleotide).
Tandem donor alternative splice sites including or excluding a single codon (GYNGYN recognition motif, where Y = pyrimidine) are less common than tandem acceptors. Signal transducer and activator of transcription 3 (STAT3) is a key gene undergoing tandem donor alternative splicing1,11. The tandem donor splice sites join exons 21 and 22 and result in the inclusion or exclusion of the codon for Serine-701 (S or &#916;S respectively)1,11. Downstream alternative acceptor sites (40 nucleotides apart from each other) joining exons 22 and 23a/b result in the inclusion of either the 55 terminal residues of the transactivation domain (&#945;) or a truncated transactivation domain with 7 unique C-terminal residues (&#946;)11. Therefore, four splice variants are possible.
STAT3 protein is a transcription factor and major signal integrator in numerous cell types12 and when mutated its constitutive activation contributes to several cancer phenotypes (reviewed in reference13). Job&#39;s Syndrome, an immunodeficiency disorder characterized by high levels of IgE, is also caused by mutations in STAT3 (reviewed in reference14). Distinct roles for STAT3 &#945; and &#946; splice variant proteins have been previously described15. Initially, STAT3 &#946; was thought to act in a dominant-negative manner16, antagonizing STAT3 &#945;&#39;s transcriptional activity, but subsequent work suggested it has independent target genes17. Despite the subtlety of tandem splicing, there is reason to believe the absence or presence of Serine-701 (Ser701) influences function. Not only is Ser701 in close proximity to Tyrosine-705 (the residue phosphorylated in STAT3 activation18), but a recent study suggests that STAT3 S and &#916;S splice variants are both necessary for viability of STAT3-addicted Diffuse Large B Cell Lymphoma (DLBLCL) cells19. The biological relevance remains to be explored. Given that splice variant composition could influence function, we endeavored to discover whether the ratio was perturbed by cytokine stimulation in eosinophils.
We initially attempted to explore the linkage between the two splicing events by using PCR specific for STAT3 &#945; and &#946; splice variants, followed by cleavage of products with a restriction enzyme specific for the S splice variants, AfeI. Densitometry of products indicated inclusion of Ser701 was roughly ten times more common than its omission (&#916;S) in both STAT3 &#945; and &#946; (data not shown). However, this semi-quantitative approach was not sufficiently reproducible, and could not be used effectively to measure all four splice variants simultaneously. To analyze proportions of each of the four splice variants, it was necessary to establish a quantitative PCR (qPCR) protocol that yielded tight technical (several assays of a given sample) replicates.
Relative qPCR relies on comparison of a gene of interest to a standard or housekeeping gene known to be expressed at a particular level20 and is appropriate when the gene of interest and housekeeping gene are amplified with similar efficiency. A double stranded (ds) DNA-binding fluorescent (cyanine) dye binds to PCR amplicons21, and after a certain number of cycles, sufficient amplification has occurred for fluorescence to be detectable. The higher the initial level of the transcript, the lower the threshold cycle (Ct) value. Since the concentration of cDNA preparations differs, one needs to compare the transcript&#39;s concentration with the concentration of a transcript known to be expressed at a consistent level in all samples, like glucuronidase-&#946; (GUSB) in eosinophils22.
Relative qPCR is not feasible for highly similar sequences, as seen in splice variants resulting from tandem splicing. The stringent conditions required to specifically amplify the splice variants result in decreased efficiency. Instead, absolute quantification must be used23. This entails preparing a standard curve with known concentrations of the spliced transcript of interest, and ensuring PCR conditions optimize specificity24. As described, absolute and relative qPCR data for a particular gene can be merged to inform understanding of the gene&#39;s expression in a particular cell type, in this case STAT3 in variously stimulated eosinophils25.
Herein, STAT3 splice variant quantification is described with the expectation that the method can be adapted to targeted studies of other tandem splicing events. Optimization was a lengthy process, where several primer pairs at various concentrations and numerous iterations of cycling parameters were tested over the course of a few months. The key features of the protocol are primer specificity validation and quantification based on standard curves with known concentrations of the splice variants. Relative qPCR in conjunction proved helpful for our application but is not necessary.

<table border="0" cellpadding="0" cellspacing="0" width="1102">
	<tbody>
		<tr height="21">
			<td height="21" style="height:21px;width:351px;">MJ Research PTC-200 Thermal Cycler</td>
			<td style="width:259px;">GMI</td>
			<td style="width:175px;">N/A</td>
			<td style="width:320px;">Used for standard PCR</td>
		</tr>
		<tr height="21">
			<td height="21" style="height:21px;">7500 Real-Time PCR System</td>
			<td style="width:259px;">Applied Biosystems</td>
			<td>N/A</td>
			<td>qPCR machine</td>
		</tr>
		<tr height="21">
			<td height="21" style="height:21px;">GS-6R</td>
			<td style="width:259px;">Beckman Coulter</td>
			<td>N/A</td>
			<td>centrifuge for 96-well plates</td>
		</tr>
		<tr height="21">
			<td height="21" style="height:21px;">Nanodrop 2000 sprectrophotometer</td>
			<td style="width:259px;">ThermoFisher Scientific</td>
			<td>N/A</td>
			<td></td>
		</tr>
		<tr height="21">
			<td height="21" style="height:21px;">RPMI-1640 medium</td>
			<td style="width:259px;">Sigma Aldrich</td>
			<td>R8758</td>
			<td>cell culture medium</td>
		</tr>
		<tr height="21">
			<td height="21" style="height:21px;">PfuTurbo DNA Polymerase</td>
			<td style="width:259px;">Agilent Technologies</td>
			<td>600410</td>
			<td>DNA polymerase for standard PCR</td>
		</tr>
		<tr height="21">
			<td height="21" style="height:21px;width:351px;">KpnI</td>
			<td>New England Biosciences</td>
			<td>R0142S</td>
			<td></td>
		</tr>
		<tr height="21">
			<td height="21" style="height:21px;width:351px;">NheI</td>
			<td>New England Biosciences</td>
			<td>R0131S</td>
			<td></td>
		</tr>
		<tr height="21">
			<td height="21" style="height:21px;">SYBR Green PCR Master Mix</td>
			<td style="width:259px;">Qiagen</td>
			<td>330523</td>
			<td>qPCR, DNA polymerase/dsDNA-binding dye mix</td>
		</tr>
		<tr height="21">
			<td height="21" style="height:21px;width:351px;">Rneasy Mini Kit</td>
			<td style="width:259px;">Qiagen</td>
			<td>74204</td>
			<td>RNA extraction kit</td>
		</tr>
		<tr height="21">
			<td height="21" style="height:21px;">SuperScript III First-Strand Synthesis System</td>
			<td style="width:259px;">Invitrogen (ThermoFisher Scientific)</td>
			<td>18080-044</td>
			<td>cDNA synthesis kit</td>
		</tr>
		<tr height="21">
			<td height="21" style="height:21px;width:351px;">Primers</td>
			<td style="width:259px;">Integrated DNA Technology</td>
			<td style="width:175px;">N/A</td>
			<td></td>
		</tr>
		<tr height="21">
			<td height="21" style="height:21px;width:351px;">NEBuffer 1.1</td>
			<td>New England Biosciences</td>
			<td>B7201S</td>
			<td></td>
		</tr>
		<tr height="21">
			<td height="21" style="height:21px;width:351px;">GenePure LE Agarose</td>
			<td>ISC BioExpress</td>
			<td>E-3120-500</td>
			<td>component of TAE gel</td>
		</tr>
		<tr height="21">
			<td height="21" style="height:21px;width:351px;">Pipettors</td>
			<td>Major lab suppliers (MLS)</td>
			<td>N/A</td>
			<td></td>
		</tr>
		<tr height="21">
			<td height="21" style="height:21px;width:351px;">Filter pipette tips</td>
			<td>Neptune Scientific</td>
			<td>BT10XL, BT20, BT200</td>
			<td></td>
		</tr>
		<tr height="21">
			<td height="21" style="height:21px;">EU One Piece Thin Wall Plate</td>
			<td style="width:259px;">MidSci</td>
			<td>ABI7501</td>
			<td></td>
		</tr>
		<tr height="21">
			<td height="21" style="height:21px;">ThermalSeal A Sealing Film</td>
			<td style="width:259px;">MidSci</td>
			<td style="width:175px;">TSA-100</td>
			<td style="width:320px;">96 well plate seal</td>
		</tr>
		<tr height="21">
			<td height="21" style="height:21px;">pET-Elmer (variant of pET-28a)</td>
			<td style="width:259px;">Novagen; modified in Mosher lab</td>
			<td style="width:175px;">N/A</td>
			<td>Details in PMID: 20947497</td>
		</tr>
		<tr height="21">
			<td height="21" style="height:21px;width:351px;">Wizard Plus SV Minipreps DNA purification system</td>
			<td style="width:259px;">Promega</td>
			<td style="width:175px;">A1460</td>
			<td>Plasmid purification</td>
		</tr>
		<tr height="21">
			<td height="21" style="height:21px;width:351px;">BigDye Terminator v3.1 Cycle Sequencing Kit</td>
			<td style="width:259px;">ThermoFisher Scientific</td>
			<td style="width:175px;">4337455</td>
			<td>Sequencing kit</td>
		</tr>
		<tr height="21">
			<td height="21" style="height:21px;width:351px;">QIAEX II Gel Extraction kit</td>
			<td style="width:259px;">Qiagen</td>
			<td style="width:175px;">20021</td>
			<td>Amplicon purification</td>
		</tr>
		<tr height="21">
			<td height="21" style="height:21px;width:351px;">DH5&#945; competent cells</td>
			<td style="width:259px;">ThermoFisher Scientific</td>
			<td style="width:175px;">18265-017</td>
			<td>available from several providers, see PMID: 2162051</td>
		</tr>
		<tr height="21">
			<td height="21" style="height:21px;width:351px;">kanamycin</td>
			<td style="width:259px;">Research Products International Corp.</td>
			<td style="width:175px;">K22000-5.0</td>
			<td></td>
		</tr>
		<tr height="21">
			<td height="21" style="height:21px;width:351px;">Tris base</td>
			<td style="width:259px;">ThermoFisher Scientific</td>
			<td style="width:175px;">BP152-5</td>
			<td>component of TAE buffer</td>
		</tr>
		<tr height="21">
			<td height="21" style="height:21px;width:351px;">Acetic acid, glacial</td>
			<td style="width:259px;">ThermoFisher Scientific</td>
			<td style="width:175px;">A38C-212</td>
			<td>component of TAE buffer</td>
		</tr>
		<tr height="42">
			<td height="42" style="height:42px;">EDTA (Ethylenediaminetetraacetic acid)</td>
			<td style="width:259px;">Sigma Chemical Company (Sigma Aldrich)</td>
			<td style="width:175px;">E-5134</td>
			<td>component of TAE buffer</td>
		</tr>
		<tr height="21">
			<td height="21" style="height:21px;width:351px;">Bacto Tryptone</td>
			<td style="width:259px;">BD Biosciences</td>
			<td style="width:175px;">211705</td>
			<td>component of Luria Broth</td>
		</tr>
		<tr height="21">
			<td height="21" style="height:21px;width:351px;">Bacto Yeast extract, technical</td>
			<td style="width:259px;">BD Biosciences</td>
			<td style="width:175px;">288620</td>
			<td>component of Luria Broth</td>
		</tr>
		<tr height="21">
			<td height="21" style="height:21px;width:351px;">Sodium chloride</td>
			<td style="width:259px;">ThermoFisher Scientific</td>
			<td style="width:175px;">S271-10</td>
			<td>component of Luria Broth</td>
		</tr>
		<tr height="21">
			<td height="21" style="height:21px;width:351px;">Sodium hydroxide</td>
			<td style="width:259px;">ThermoFisher Scientific</td>
			<td style="width:175px;">SS255-1</td>
			<td>component of Luria Broth</td>
		</tr>
		<tr height="21">
			<td height="21" style="height:21px;width:351px;">Bacto Agar</td>
			<td style="width:259px;">BD Biosciences</td>
			<td>214010</td>
			<td>component of Luria Broth plate</td>
		</tr>
		<tr height="42">
			<td height="42" style="height:42px;">Lasergene SeqBuilder</td>
			<td style="width:259px;">DNASTAR</td>
			<td style="width:175px;"></td>
			<td style="width:320px;">Figure 1 generated using Lasergene SeqBuilder software version 12.2.0 (DNASTAR)</td>
		</tr>
	</tbody>
</table>

merging absolute and relative quantitative pcr data to quantify stat3 splice variant transcripts

Human signal transducer and activator of transcription 3 (STAT3) is one of many genes containing a tandem splicing site. Alternative donor splice sites 3 nucleotides apart result in either the inclusion (S) or exclusion (&#916;S) of a single residue, Serine-701. Further downstream, splicing at a pair of alternative acceptor splice sites result in transcripts encoding either the 55 terminal residues of the transactivation domain (&#945;) or a truncated transactivation domain with 7 unique residues (&#946;). As outlined in this manuscript, measuring the proportions of STAT3's four spliced transcripts (S&#945;, S&#946;, &#916;S&#945; and &#916;S&#946;) was possible using absolute qPCR (quantitative polymerase chain reaction). The protocol therefore distinguishes and measures highly similar splice variants. Absolute qPCR makes use of calibrator plasmids and thus specificity of detection is not compromised for the sake of efficiency. The protocol necessitates primer validation and optimization of cycling parameters. A combination of absolute qPCR and efficiency-dependent relative qPCR of total STAT3 transcripts allowed a description of the fluctuations of STAT3 splice variants' levels in eosinophils treated with cytokines. The protocol also provided evidence of a co-splicing interdependence between the two STAT3 splicing events. The strategy based on a combination of the two qPCR techniques should be readily adaptable to investigation of co-splicing at other tandem splicing sites.

Good quality qPCR data will generate a sigmoidal amplification plot (Figure 2a), signifying exponential increase in transcripts over the course of cycling. The presence of too much template can result in a high fluorescence background, meaning an inappropriate baseline is established in the first few cycles. If the data do not provide an exponential curve (Figure 2b), further optimization is necessary (outlined in steps 3.1 and 3.4). For further information about troubleshooting qPCR results, refer to reference32. The standard curves generated for the template calibrator plasmids will indicate the efficiency of amplification (curve for STAT3 S&#945; shown in Figure 3). Efficiencies between 83 and 95% were observed under the conditions described. The equation for specificity (Step 3.4) assumes equal efficiency, which is unlikely25, so the specificity is likely to be greater than specificity factor suggests.
In order to evaluate the congruency of STAT3 levels, the absolute values of each of the four splice variants were measured as well as the level of total STAT3, the latter using primers amplifying a region common to all four splice variants (Figure 4). Ideally the linear regression (indication of correlation) and slope (ratio of pan-STAT3 to summed splice variants) should both be close to 1.
The absolute qPCR data are presented as pie charts to show the proportions of the four splice variants over time post-stimulation with cytokines (Figure 5a). Resting eosinophils (0 hr) had the smallest proportion of STAT3 S&#945;, although this variant was always the most abundant. Multiplying the fraction of STAT3 &#946; splice variants (S&#946;+&#916;S&#946;) by the fraction of STAT3 &#916;S splice variants (&#916;S&#945;+&#916;S&#946;) consistently gave a value lower than the experimentally-recorded value for &#916;S&#946;. If the splicing events were independent, one would expect that multiplying fraction of &#916;S variants by the fraction of &#946; variants would give a value that agrees with the experimentally-determined value. Higher levels of &#916;S&#946; than expected from independent splicing events suggests a co-splicing bias exists.
Merging absolute and relative qPCR data demonstrated that levels of all STAT3 splice variants increased post-stimulation with cytokines IL3 and TNF&#945;, with levels peaking 6 hr post-stimulation (Figure 5b-e). For three of the four splice variants, transcript levels were roughly 3 times higher in IL3+TNF&#945; treated eosinophils (6 hr) compared to eosinophils in media at the same time point. STAT3 S&#945; levels were 3.5 times higher in IL3+TNF&#945; treated eosinophils compared to eosinophils in media at this time point. The greatest uncertainty (largest standard error of measurement) was seen in &#916;S&#946; (Figure 5e), which comprises the smallest fraction of total STAT3 in all samples. This was not surprising, as lower levels are associated with higher Ct values. Requiring more cycles to reach the threshold cycle will compound uncertainty due to cycle-to-cycle variation in efficiency of amplification.
<img alt="Figure 1" src="/files/ftp_upload/54473/54473fig1.jpg" /> 
Figure 1: Schematic of primer pairs used to perform qPCR of STAT3 splice variants and pan-STAT3. Primers used to specifically amplify each of the STAT3 splice variants (S&#945;, S&#946;, &#916;S&#945; and &#916;S&#946; respectively) are shown. Forward primers (STAT3 &#34;S&#34; and &#34;&#916;S&#34;) span the junction between exons 21 and 22. <a href="https://www.jove.com/files/ftp_upload/54473/54473fig1large.jpg" target="_blank">Please click here to view a larger version of this figure.</a>
<img alt="Figure 2" src="/files/ftp_upload/54473/54473fig2.jpg" /> 
Figure 2: Amplification plots of qPCR data. (a) Sigmoidal amplification plots means reliable amplification. These data were obtained from qPCR of two serial dilutions of plasmid containing STAT3 S&#945;, with each pair of colored lines representing the fluorescence levels of duplicate diluted samples over the course of 40 cycles (x-axis). The most concentrated sample (green-grey) was sufficiently amplified by cycle 17 (with dsDNA-binding dye proportional to fluorescence, shown on the y-axis) to exceed the threshold fluorescence value (baseline shown as green arrow). Its Ct value would be 17. (b) Non-exponential plots suggest that the background-fluorescence threshold was not correctly established in the first few cycles. This could be due to the presence of an inhibitor, or highly-concentrated template or primers. <a href="https://www.jove.com/files/ftp_upload/54473/54473fig2large.jpg" target="_blank">Please click here to view a larger version of this figure.</a>
<img alt="Figure 3" src="/files/ftp_upload/54473/54473fig3.jpg" /> 
Figure 3: Standard curve of log (copy number of STAT3 S&#945;) vs Ct. There is a linear relationship between the log of each STAT3 splice variant&#39;s copy number and threshold cycle (Ct). Creating a standard curve from plasmid DNA imitates the cDNA of samples and thus provides a better measurement than a curve created from diluted PCR amplicons. The data presented are Ct values obtained from qPCR of two serial dilutions of plasmid containing STAT3 S&#945;. From this curve, the copy number present in each sample can be interpolated, and amplification efficiency calculated (83.9%). Although y-intercepts are less reproducible than slope, the intercept suggests 42.2 cycles would be necessary to be certain no target DNA is present. Error bars indicate SEM, n = 3. Comparable curves were constructed for STAT3 S&#946;, &#916;S&#945; and &#916;S&#946; (not shown). <a href="https://www.jove.com/files/ftp_upload/54473/54473fig3large.jpg" target="_blank">Please click here to view a larger version of this figure.</a>
<img alt="Figure 4" src="/files/ftp_upload/54473/54473fig4.jpg" /> 
Figure 4: Comparison of quantified pan-STAT3 vs cumulative STAT3 splice variants. The regression of added STAT3 splice variants vs total STAT3 should have slope (ratio of pan- to cumulative) and R2 value (correlation) close to 1. Values from 17 samples (eosinophils and DLBCL) included. Error bars indicate SEM of x-to-y determinations, n &#8805; 2 for each. Figure adapted from reference20. <a href="https://www.jove.com/files/ftp_upload/54473/54473fig4large.jpg" target="_blank">Please click here to view a larger version of this figure.</a>
<img alt="Figure 5" src="/files/ftp_upload/54473/54473fig5.jpg" /> 
Figure 5: STAT3 splice variant levels fluctuated over the course of cytokine treatment. (a) Pie charts indicating percentages of each STAT3 splice variant in eosinophils during treatment with IL3 and TNF&#945;. (b-e) Changes in STAT3 splice variants in eosinophils treated with various combinations of cytokines, measured by combining relative and absolute qPCR data. STAT3 S&#945; (b), S&#946; (c), &#916;S&#945; (d), and &#916;S&#946; (e) levels fluctuated over time post-stimulation. Levels initially increased, peaking 6 hr post-stimulation. The IL3+TNF&#945; combination elicited higher expression of all four STAT3 splice variants than IL3 alone. SEM calculated for each data point accounting for propagation of error. <a href="https://www.jove.com/files/ftp_upload/54473/54473fig5large.jpg" target="_blank">Please click here to view a larger version of this figure.</a>
<img alt="Table 1" src="/files/ftp_upload/54473/54473tbl1.jpg" /> 
Table 1: Primers used for amplification (a), absolute (b) and relative (c) quantitative PCR. Cloning primers have restriction sequences and 5&#39;-extensions for efficient cutting. KpnI and NheI restriction sites are indicated in bold. <a href="https://www.jove.com/files/ftp_upload/54473/54473tbl1large.jpg" target="_blank">Please click here to view a larger version of this table.</a>
<img alt="Table 2" src="/files/ftp_upload/54473/54473tbl2.jpg" /> 
Table 2. PCR cycling parameters (left) and reagent volumes (right) for amplification (a), relative (b), absolute (c) quantitative PCR. <a href="https://www.jove.com/files/ftp_upload/54473/54473tbl2large.jpg" target="_blank">Please click here to view a larger version of this table.</a>
<img alt="Table 3" src="/files/ftp_upload/54473/54473tbl3.jpg" /> 
Table 3: Template for absolute qPCR plasmid calibration assay. This assay is necessary to assess reproducibility and efficiency, as well as generating standard curves from which to interpolate data. The &#34;non-target&#34; mixes will give an estimate of specificity. Optimization may be necessary to achieve consistency. <a href="https://www.jove.com/files/ftp_upload/54473/54473tbl3large.jpg" target="_blank">Please click here to view a larger version of this table.</a>
<img alt="Table 4" src="/files/ftp_upload/54473/54473tbl4.jpg" /> 
Table 4: Template for relative qPCR (pan-STAT3 and housekeeping gene GUSB) calibration assay. Unlike absolute qPCR, the point of this assay is to determine conditions at which amplification efficiency is ~100%. <a href="https://www.jove.com/files/ftp_upload/54473/54473tbl4large.jpg" target="_blank">Please click here to view a larger version of this table.</a>
<img alt="Table 5" src="/files/ftp_upload/54473/54473tbl5.jpg" /> 
Table 5: Template for relative qPCR sample assay for measuring pan-STAT3 and housekeeping gene GUSB. Standard curves are repeated together with samples to ensure comparable efficiency of the assays. <a href="https://www.jove.com/files/ftp_upload/54473/54473tbl5large.jpg" target="_blank">Please click here to view a larger version of this table.</a>
<img alt="Table 6" src="/files/ftp_upload/54473/54473tbl6.jpg" /> 
Table 6: Template for absolute qPCR sample assay for measuring S variants. Standard curves are repeated together with samples to ensure comparable efficiency of the assays. <a href="https://www.jove.com/files/ftp_upload/54473/54473tbl6large.jpg" target="_blank">Please click here to view a larger version of this table.</a>
<img alt="Table 7" src="/files/ftp_upload/54473/54473tbl7.jpg" /> 
Table 7: Template for absolute qPCR sample assay for measuring &#916;S variants. Standard curves are repeated together with samples to ensure comparable efficiency of the assays. <a href="https://www.jove.com/files/ftp_upload/54473/54473tbl7large.jpg" target="_blank">Please click here to view a larger version of this table.</a>
<img alt="Table 8" src="/files/ftp_upload/54473/54473tbl8.jpg" /> 
Table 8: Template for absolute qPCR sample assay for measuring pan-STAT3. <a href="https://www.jove.com/files/ftp_upload/54473/54473tbl8large.jpg" target="_blank">Please click here to view a larger version of this table.</a>

Watch this Scientific Journal Video about Merging Absolute and Relative Quantitative PCR Data to Quantify STAT3 Splice Variant Transcripts at JoVE.com

Merging Absolute and Relative Quantitative PCR Data to Quantify STAT3 Splice Variant Transcripts

Tandem splicing events occur at sites less than 12 nucleotides apart. Quantifying ratios of such splice variants is feasible using an absolute quantitative PCR approach. This manuscript describes how splice variants of the gene STAT3, in which two splicing events results in Serine-701 inclusion/exclusion and &#945;/&#946; C-termini, can be quantified.

Human signal transducer and activator of transcription 3 (STAT3) is one of many genes containing a tandem splicing site. Alternative donor splice sites 3 ...

Research

Education

JoVE Journal

JoVE Core

Molecular Biology

Genetics

Unit 1

Additional Roles of RNA

SCIENTISTS IN ACTION

Detection of microRNA Expression in Peritoneal Membrane of Rats Using Quantitative Real-time PCR

MicroRNAs (miRNAs) are small noncoding RNAs that regulate messenger RNA expression post-transcriptionally. The miRNA expression profile has been ...

Bidirectional Retroviral Integration Site PCR Methodology and Quantitative Data Analysis Workflow

Integration Site (IS) assays are a critical component of the study of retroviral integration sites and their biological significance. In recent retroviral ...

Adapting 3' Rapid Amplification of CDNA Ends to Map Transcripts in Cancer

Maturation of eukaryotic mRNAs involves 3&#39; end formation, which involves the addition of a poly(A) tail. In order to map the 3&#39; end of a gene, the ...

A Fast and Quantitative Method for Post-translational Modification and Variant Enabled Mapping of Peptides to Genomes

Cross-talk between genes, transcripts, and proteins is the key to cellular responses; hence, analysis of molecular levels as distinct entities is slowly ...

Determining the Likelihood of Variant Pathogenicity Using Amino Acid-level Signal-to-Noise Analysis of Genetic Variation

Advancements in the cost and speed of next generation genetic sequencing have generated an explosion of clinical whole exome and whole genome testing. ...

Use of Single Molecule Fluorescent In Situ Hybridization (SM-FISH) to Quantify and Localize mRNAs in Murine Oocytes

Use of Single Molecule Fluorescent In Situ Hybridization SM-FISH to Quantify and Localize mRNAs in Murine Oocytes

Current methods routinely used to quantify mRNA in oocytes and embryos include digital reverse-transcription polymerase chain reaction (dPCR), ...

Discrimintion and Mapping of the Primary and Processed Transcripts in Maize Mitochondrion Using a Circular RT-PCR-based Strategy

In plant mitochondria, some steady-state transcripts have 5&#39; triphosphate derived from transcription initiation (primary transcripts), while the ...

Detection of Human Immunodeficiency Virus Type 1 (HIV-1) Antisense Protein (ASP) RNA Transcripts in Patients by Strand-Specific RT-PCR

Detection of Human Immunodeficiency Virus Type 1 HIV-1 Antisense Protein ASP RNA Transcripts in Patients by Strand-Specific RT-PCR

In retroviruses, antisense transcription has been described in both human immunodeficiency virus type 1 (HIV-1) and human T-lymphotropic virus 1 (HTLV-1). ...

High-Throughput Quantitative RT-PCR in Single and Bulk C. elegans Samples Using Nanofluidic Technology

High-Throughput Quantitative RT-PCR in Single and Bulk C. elegans Samples Using Nanofluidic Technology

This paper presents a high-throughput reverse transcription quantitative PCR (RT-qPCR) assay for Caenorhabditis elegans that is fast, robust, and highly ...

Quantitative Real-Time PCR Evaluation of microRNA Expressions in Mouse Kidney with Unilateral Ureteral Obstruction

MicroRNAs (miRNAs) are single stranded, non-coding RNA molecules that typically regulate gene expression at the post-transcriptional level by binding to ...