Clinical Chemistry Link to Randox Laboratories Web Site
HOME HELP FEEDBACK SUBSCRIPTIONS ARCHIVE SEARCH TABLE OF CONTENTS
QUICK SEARCH:   [advanced]


     

Clinical Chemistry 52: 1988-1996, 2006. First published September 21, 2006; 10.1373/clinchem.2006.073064
This Article
Right arrow Abstract Freely available
Right arrow Full Text (PDF)
Right arrow 073064.Supplemental Data Table 1
Right arrow All Versions of this Article:
clinchem.2006.073064v1
52/11/1988    most recent
Right arrow Submit an electronic Letter to
the Editor about this paper
Right arrow Alert me when this article is cited
Right arrow Alert me when eLetters are posted
Right arrow Alert me if a correction is posted
Right arrow Citation Map
Services
Right arrow Email this article to a friend
Right arrow Similar articles in this journal
Right arrow Similar articles in ISI Web of Science
Right arrow Similar articles in PubMed
Right arrow Alert me to new issues of the journal
Right arrow Download to citation manager
Right arrow reprints & permissions
Citing Articles
Right arrow Citing Articles via ISI Web of Science (1)
Right arrow Citing Articles via Google Scholar
Google Scholar
Right arrow Articles by Stenman, J.
Right arrow Articles by Orpana, A.
Right arrow Search for Related Content
PubMed
Right arrow PubMed Citation
Right arrow Articles by Stenman, J.
Right arrow Articles by Orpana, A.
Related Collections
Right arrow Molecular Diagnostics and Genetics
(Clinical Chemistry. 2006;52:1988-1996.)
© 2006 American Association for Clinical Chemistry, Inc.

Molecular Diagnostics and Genetics

Targeted Gene-Expression Analysis by Genome-Controlled Reverse Transcription-PCR

Jakob Stenman1,2,a,3, Annukka Paju3,3, Oso Rissanen3, Tuomas Tenkanen4,3, Caj Haglund5,6, Jari Räsänen7, Jarmo Salo7, Ulf-Håkan Stenman3,8 and Arto Orpana3,8,9

1 Hospital for Children and Adolescents, University of Helsinki, Helsinki, Finland.
2 Department of Surgery, HUCH, Jorvi Hospital, Espoo, Finland.
3 Departments of Clinical Chemistry,5 Surgery,6 Pathology, and9 Medical Genetics, University of Helsinki, Helsinki, Finland.
4 Finnzymes Oy, Espoo, Finland.
7 Division of General Thoracic and Esophageal Surgery, Department of Cardiothoracic Surgery, HUCH, Helsinki, Finland.
8 HUCH Laboratory Diagnostics, Helsinki, Finland.

aAddress correspondence to this author at: Hospital for Children and Adolescents, University of Helsinki, P.O. Box 281, FIN 00029 HUS, Finland. Fax 358-9-4717-1731; e-mail jakob.stenman{at}helsinki.fi.


   Abstract
Top
 Abstract
Introduction
 Materials and Methods
 Results
 Discussion
 References
 
Background: For gene-expression analysis, which is anticipated to play an important role in classification of tumors and premalignant conditions, PCR-based quantitative assays must have increased diagnostic quantitative accuracy and reproducibility and enable analysis of gene expression in formalin-fixed paraffin-embedded (FFPE) tissue samples.

Methods: We developed a reverse transcription–PCR-based quantitative assay that modifies the cDNA sequence to increase the melting temperature of short (56–64 bp) PCR amplicons, enabling their quantification in-tube by homogeneous melting-curve analysis. We used this method to analyze the expression of 8 genes, 7 potential colon cancer markers, and 1 control in samples obtained from 3 colon carcinoma cell lines, endoscopic biopsy from 8 patients undergoing gastroscopy for Barrett esophagus, and archival FFPE and frozen tissue from 20 patients who underwent surgery for colon carcinoma.

Results: The detection limit of the assay, when optimized for FFPE samples, was 100 copies of cDNA, and the dynamic range was 3 orders of magnitude. A prototype assay containing a panel of 8 genes displayed good reproducibility compared with the commercially available TaqMan® assay (interassay CVs, 5%–20% vs 7%–43%, respectively). Gene-expression analysis was performed successfully in 26 (96%) of 27 endoscopic biopsy specimens, 30 (86%) of 35 archival FFPE samples, and 20 (100%) of 20 archival frozen samples.

Conclusions: This new technology combines the reproducibility of competitive PCR with accurate quantitative detection by in-tube melting-curve analysis, enabling efficient analysis of mRNA profiles in samples with small numbers of cells or small amounts of tissue, as well as in archival FFPE tissues.


   Introduction
Top
 Abstract
 Introduction
Materials and Methods
 Results
 Discussion
 References
 
High-density microarrays are extensively used to define subsets of differentially expressed genes suitable for classification of various forms of cancer (1)(2)(3). Reverse transcription–PCR (RT-PCR)1 has been successfully used for studying the expression of a limited number of genes in a large number of clinical samples to evaluate the clinical relevance of candidate genes (4)(5)(6). Increasing the usefulness of RT-PCR technology in this area, however, requires the achievement of sufficient quantitative accuracy and reproducibility for diagnostic purposes and development of assays for analysis of gene expression in formalin-fixed paraffin-embedded (FFPE) tissue samples (7)(8).

The first truly quantitative gene-expression assays were based on competitive coamplification by RT-PCR of the wild-type cDNA template with an in vitro-generated RNA (9)(10) or DNA (11) internal standard and, later, universal internal standards containing concatenated primer sequences for multiple genes (12). The use of competitive RT-PCR decreased in favor of real-time PCR (13)(14) until Sequenom® Corporation launched its quantitative gene-expression application for the MassARRAY® platform (15). We have shown that highly accurate end-point detection can be achieved in a competitive RT-PCR system with highly homologous competing templates (16), provided that the amplification reaction originates from a sufficient amount of cDNA molecules (17).

Although shorter amplicons are essential for analysis of the fragmented mRNA in FFPE samples (4), RT-PCR amplicons are generally designed to span at least 1 exon–exon boundary to enable separation of the mRNA-derived amplicons from genomic DNA (gDNA), complicating the design of PCR amplicons shorter than ~100 nucleotides. Different approaches have been used to enable differential competitive amplification of templates derived from cDNA and gDNA within the boundaries of a single exon, including a technique for measuring the amount of a specific mRNA transcript in relation to the amount of gDNA (18) and a similar approach to enable differential amplification of viral RNA and DNA lacking introns (19). Taking this concept further, we developed an RT-PCR technique for quantitative measurement of the relative amounts of multiple mRNA transcripts in a cell or in a tissue sample containing as few as 1000 cells (20). We analyzed the expression of 8 genes. Seven of these genes, trypsinogen (Try or PRSS1, protease serine 1)2 , matrix-metalloproteinase-2 (MMP2), MMP7, MMP9, urokinase-type plasminogen activator (uPA or PLAU, plasminogen activator urokinase), uPA receptor (uPAR or PLAUR, plasminogen activator urokinase receptor), and mutS homolog 2 (MSH2) are potential markers for colon cancer (21)(22)(23)(24). Glyceraldehyde-3-phosphate dehydrogenase (GAPD) was included in the assay as a positive control.


   Materials and Methods
Top
 Abstract
 Introduction
 Materials and Methods
Results
 Discussion
 References
 
cell culture
COLO 205, DLD-1, and HT-29 cells were obtained from American Type Culture Collection (Rockville, MD). The cells were grown at 37 °C in a humidified atmosphere with 5% CO2. COLO 205 cells were maintained in RPMI 1640 medium supplemented with 10% heat-inactivated fetal bovine serum, 2 mmol/L L-glutamine, 100 000 units/L penicillin, and 100 mg/L streptomycin. DLD-1 cells were maintained in RPMI 1640 medium supplemented with 10% heat-inactivated fetal bovine serum, 2 mmol/L L-glutamine, 100 000 units/L penicillin, 100 mg/L streptomycin, 14 mmol/L sodium bicarbonate, 25 mmol/L glucose, 10 mmol/L HEPES, and 1 mmol/L sodium pyruvate. HT-29 cells were maintained in McCoy 5a medium supplemented with 10% heat-inactivated fetal bovine serum, 1.5 mmol/L L-glutamine, 100 000 units/L penicillin, and 100 mg/L streptomycin.

patient samples
We obtained endoscopic biopsy samples from 8 patients undergoing gastroscopy for intestinal metaplasia of the lower esophagus (Barrett esophagus) at Helsinki University Central Hospital. Biopsy specimens ~1–5 mm3 in size were taken from the unaffected upper part of the esophagus, the area of intestinal metaplasia in the lower part of the esophagus, and from the corpus of the stomach. We also used archival FFPE and frozen tissue samples obtained during 2000–2004 from 20 patients who underwent surgery for colon carcinoma. The study was carried out with the approval of the ethics committee of the Helsinki University Central Hospital. Informed consent was obtained from all study patients.

Tissue specimens were macroscopically divided in the operating room into parallel samples for freezing and formalin fixation. The frozen tissue samples were primarily snap frozen in liquid nitrogen and stored at –80 °C for 5 years. The corresponding FFPE samples were primarily fixed in formalin, embedded in paraffin, and stored at room temperature for the same time period. RNA was extracted from 10 8-µm FFPE whole sections and 10 10-µm frozen sections, and tissue morphology was evaluated from an adjacent hematoxylin and eosin–stained section (FFPE samples) or toluidine blue–stained section (frozen samples).

in vitro synthesis of CRNA
We obtained full-length Try-2 and GAPD cDNA clones from American Type Culture Collection and MMP2, MMP7, MMP9, uPA, uPAR, and MSH2 clones from the Incyte Genomics IRAL cDNA clone library (Incyte Genomics). We confirmed sequences of the cDNA clones by sequencing from both ends on an ABI Prism 310 genetic analyzer (Applied Biosystems) with the ABI Prism dye terminator cycle sequencing core method and AmpliTaq DNA polymerase. With the cDNA clones as templates, we generated cRNA by in vitro transcription with a MEGAscript High Yield Transcription Kit and purification with MEGAclear (Ambion).

rna extraction from cultured cells, endoscopic biopsy specimens, and frozen tissue sections
With the RNeasy Mini Kit (Qiagen), used according to the manufacturer’s instructions, we extracted total RNA from the cultured cells, the 1–5 mm3 endoscopic biopsies, and 10 10-µm sections of each frozen sample. Total RNA extracted from cultured cells was DNase-treated by adding to 10 µg of RNA, 5 µL 10x DNase buffer (100 mmol/L Tris, pH 7.5, 25 mmol/L MgCl2, 5 mmol/L CaCl2), and 10 units of DNase I (Roche Applied Science) and incubating for 30 min at 37 °C, after which the DNase was inactivated for 5 min at 70 °C. RNA was then precipitated with 5 µL of 2 mol/L sodium acetate (pH 4.0) and 140 µL of 100% ethanol overnight at –20 °C, after which the sample was centrifuged at 5200g for 30 min at 4 °C. The supernatant was removed carefully and the RNA pellet was washed with 200 µL of 70% ethanol, after which the sample was centrifuged at 5200g for 20 min at 4 °C. The supernatant was removed, and the RNA was dissolved in 20 µL of RNase-free water. RNA extracted from endoscopic biopsies and frozen sections was DNase-treated during the extraction on the membrane of the RNeasy column, according to the manufacturer’s instructions.

rna extraction from ffpe tissue samples
RNA was extracted from 10 8-µm FFPE tissue sections of each sample. Paraffin was dissolved in octane (Sigma) for 2 x 10 min at 37 °C and then washed for 2 min with absolute ethanol. After centrifugation at 2000g for 1.5 min, we added 300 µL of 4 mol/L guanidium isothiocyanate buffer and 1 µg of yeast tRNA (Sigma). The samples were incubated at 95 °C for 30 min, followed by centrifugation at 5200g for 2 min (25). Total RNA was extracted by use of a modification of the protocol of Chomczynski and Sacchi as follows (26). Thirty µL of 2 mol/L sodium acetate and 300 µL of water-saturated phenol were added to the sample, after which the sample was mixed thoroughly by repeated inversion and incubated on ice for 30 min. We transferred the sample to a prespun Phase Lock Gel 2 mL Heavy tube (Eppendorf), added 90 µL of chloroform–isoamyl alcohol (49:1), and mixed by repeated gentle inversion. After centrifugation at 5200g for 5 min to separate the phases, we added 300 µL of phenol–chloroform–isoamyl alcohol (50:49:1) to the aqueous phase and mixed thoroughly. After centrifuging at 5200g for 5 min, the resulting aqueous phase was collected and the RNA was precipitated by adding a double volume of absolute ethanol and incubating at –70 °C overnight. The precipitate was washed twice with 500 µL of 70% ethanol and centrifuged at 5200g for 10 min at 4 °C to repellet the sample. After the final wash, the pellet was dried at room temperature and dissolved in 20 µL of RNase-free water. The extracted total RNA was then DNase-treated as described for cultured cells.

reverse transcription
We performed reverse transcription with the DyNAmo SYBR Green 2-step qRT-PCR Kit (Finnzymes) according to the manufacturer’s instructions. Briefly, 1 µg of total RNA from cultured cells, 0.5–1 µg of total RNA from endoscopic biopsies, 60–200 ng of total RNA from frozen sections, or 5 µL of FFPE-derived RNA (the concentration of which could not be directly measured) was transcribed in a 20-µL reaction volume containing 1x reverse transcription buffer, 10 pmol of each of the sequence-modifying antisense primers (see Table 1 in the Data Supplement that accompanies the online version of this article at http://www.clinchem.org/content/vol52/issue11), and 2 µL of M-MuLV RNase H+ reverse transcriptase. The reverse transcription reaction was incubated at 42 °C for 50 min, and then the reverse transcriptase was inactivated for 5 min at 85 °C.

exonuclease i treatment
After reverse transcription, the unbound sequence-modifying reverse transcription primers were degraded by exonuclease I (New England Biolabs) in a 25-µL reaction volume in 1x restriction buffer (67 mmol/L glycineKOH, 6.7 mmol/L MgCl2, 10 mmol/L ß-mercaptoethanol, pH 9.5) and with 2 units of exonuclease I at 37 °C for 60 min, followed by denaturation at 80 °C for 20 min.

polymerase chain reaction
We amplified 2.5 µL of the exonuclease I-treated reverse transcription products together with 2 ng (FFPE samples) or 20 ng (cultured cells and endoscopic biopsies) of human gDNA (Roche Applied Science) in gene-specific amplification reactions with the DyNAmo SYBR Green 2-step qRT-PCR Kit (Finnzymes) according to the manufacturer’s instructions. For PCR primers used, see Table 1 in the online Data Supplement. We used locked nucleic acids to optimize the melting temperature of some primers (27). We tested primer pairs by analyzing the amplification products by standard agarose gel electrophoresis and by melting-curve analysis, as described. The accepted primer pairs showed no nonspecific or dimer amplification on an agarose gel and no nonspecific peaks in the temperature range of the specific cDNA and gDNA peaks in the melting-curve analysis. PCR was performed on an ABI 7500 Fast (Applied Biosystems) instrument for 30 cycles (cultured cells and endoscopic biopsies) or 35 cycles (FFPE samples), with denaturation at 95 °C (15-s hold) and annealing at 57 °C (1-min hold). Fluorescence was measured at the end of annealing.

melting curve acquisition and analysis
Immediately after amplification, we detected and quantified amplification products by melting-curve analysis on the ABI 7500 Fast instrument, with an additional denaturation at 95 °C and continuous melting curve acquisition from 57 °C to 98 °C with a 0.1 °C/s ramp rate. We obtained a derivative melting curve plot by use of the default settings of the instrument. We determined background fluorescence for the 1st peak by forward baseline extension and for the 2nd peak by backward extension from a 5 °C temperature range adjacent to the peak. The areas of the cDNA- and gDNA-derived melting curve peaks were calculated after subtraction of the background, by integrating the negative derivative raw data measuring points at 0.3 °C intervals over a temperature range of 5 °C. The cDNA:gDNA ratio was calculated from these peak areas.

determination of the dynamic range of the assay
The measuring range of the assay was determined by analyzing a 10-fold dilution series of sequence-modified cDNA reverse-transcribed from 101 copies of cRNA of single genes in relation to the reference DNA template containing 20 ng of human gDNA corresponding to ~104 copies of each gene.

determination of the imprecision of the genome-controlled rt-pcr and taqman® assays
To determine the intra- and interassay CVs of the genome-controlled RT-PCR assay, we used 20 ng of gDNA as an internal calibrator and performed 8 parallel multiplexed cDNA synthesis reactions from total RNA extracted from COLO 205, HT-29, and DLD-1 cells. To determine intraassay variation, we divided pooled cDNA synthesis products treated with exonuclease I into 8 parallel sets of gene-specific amplification reactions [2.5 µL (104 cells) per reaction]. Each set contained gene-specific amplification reactions for each of the analyzed genes. For determining interassay variation, 8 separate exonuclease I–treated cDNA synthesis products were transferred to 8 parallel sets of amplification reactions. With 2 ng of gDNA as an internal calibrator, we determined the interassay variation of the genome-controlled RT-PCR assay for FFPE samples by performing 8 separate multiplexed cDNA synthesis reactions from total RNA extracted from 8 FFPE samples. Each sample consisted of 10 8-µm sections from the same archival colon cancer paraffin specimen.

To determine imprecision of the TaqMan assay, we used the High-Capacity cDNA Archive Kit (Applied Biosystems) according to the manufacturer’s instructions to perform single (intraassay) or 8 parallel (interassay) cDNA synthesis reactions with random hexamers. To amplify cDNA reverse transcribed from 104 cells, we used TaqMan Universal PCR Master Mix (Applied Biosystems) according to the manufacturer’s instruction and inventoried TaqMan assays Hs00605631_g1 for Try-1, Hs00234422_m1 for MMP2, Hs00159163_m1 for MMP7, Hs00234579_m1 for MMP9, Hs00170182_m1 for uPA, Hs00182181_m1 for uPAR, Hs00179887_m1 for MSH2, Hs99999905_m1 for GAPD (as an external standard), and GAPD endogenous control (VIC/TAMRA Probe).

statistical analysis
We used the Mann–Whitney U-test (2-sided) to estimate the statistical significance of the differences in the expression of the analyzed genes.


   Results
Top
 Abstract
 Introduction
 Materials and Methods
 Results
Discussion
 References
 
genome-controlled rt-pcr assay procedure
The 3-step analytical procedure is shown in Fig. 1 (20). The structure of the sequence-modifying reverse transcription primer is shown in Fig. 2 . Detection of the coamplified sequence-modified cDNA and gDNA templates by melting-curve analysis is shown in Fig. 3 .


Figure 1
View larger version (17K):
[in this window]
[in a new window]
  		Figure 1. Principle of the genome-controlled RT-PCR assay.

Step 1: modification of the cDNA sequence during reverse transcription. Reverse transcription of the mRNA of the sample was carried out with gene-specific sequence-modifying primers. This is performed for all genes of interest in a single reaction. A sequence modification of 5 or 6 A or T to G or C nucleotide substitutions was generated in each of the cDNAs during reverse transcription. The nucleotide sequence modifications were located in such a manner that they would be incorporated in the amplification products but do not affect the nucleotide sequence of the binding sites of the PCR primers. These sequence modifications were incorporated in the cDNA-derived amplicon and led to a 4–5 °C increase in melting temperature in relation to the corresponding gDNA-derived amplicon. The reaction was carried out under nondenaturing conditions under which any DNA contained in the sample remained double-stranded and did not function as a template for cDNA synthesis. Step 2: inactivation of redundant primers after reverse transcription. After completion of the reverse transcription, remaining sequence-modifying primers were removed by treatment with single-stranded DNA-specific exonuclease, which left double-stranded mRNA-cDNA hybrids intact. This step is essential to prevent the primers from being transferred to the PCR step, as this would cause modification of the gDNA internal standard and nonspecific amplification because of primer-dimers. Step 3: coamplification of the sequence-modified cDNA and gDNA templates and melting-curve analysis of the PCR amplicons. The gDNA was added to the multiplexed reverse transcription reaction before it was divided into individual gene-specific amplification reactions. In each of these amplification reactions, the sequence-modified cDNA was coamplified in a competitive manner with a universal internal standard consisting of gDNA with a single pair of gene-specific PCR primers. After amplification, the ratio of the cDNA and gDNA-derived amplification products was measured by melting-curve analysis (28) (Figs. 3Up and 4Up ). The gene-specific cDNA/gDNA ratios determined in each amplification reaction were combined to generate a sample-specific profile, which reflected the relative amounts of mRNA transcripts present in the sample (Fig. 5Up ).


Figure 2
View larger version (7K):
[in this window]
[in a new window]
  		Figure 2. Structure of the sequence-modifying reverse transcription primers.

The sequence-modifying primers comprise 3 functional segments; a 5'-terminal segment containing the entire nucleotide sequence of the antisense primers used in subsequent amplification reactions, a central segment containing 5 or 6 A or T to C or G substitutions, giving rise to a 4–5 °C increase in melting temperature of the subsequent cDNA-derived PCR amplicon, and a 3'-terminal segment complementary to the mRNA sequence of the corresponding gene to ensure specific annealing to the target mRNA.


Figure 3
View larger version (31K):
[in this window]
[in a new window]
  		Figure 3. Amplification and melting-curve analysis of a 10:1 dilution series of cDNA reverse transcribed from 1010 copies of Try-1 cRNA.

(A), real time amplification curves from 8 separate gene-specific amplification reactions containing 10-fold dilutions, corresponding to 10–108 copies of Try cRNA, and 20 ng of gDNA. (B), melting curves obtained from the same reactions after amplification. A temperature shift of 4–5 °C allowed for sufficient separation of the gDNA- and cDNA-derived PCR amplicons for all of the analyzed genes.

measuring range and detection limit
For analysis of fresh tissue samples or cultured cells with 20 ng of gDNA (containing ~104 copies of each gene) as an internal standard, the linear range of the assay was ~3 orders of magnitude (interval, 103 to ~106 copies of target cDNA) (Fig. 4 ). The detection limit of the assay was 1000 copies of target cDNA, corresponding to a cDNA:gDNA ratio of 0.1. Above this limit, the CVs of the assay were ≤20% for all the analyzed genes (Figs. 4 and 5 ).


Figure 4
View larger version (15K):
[in this window]
[in a new window]
  		Figure 4. Dynamic range and detection limit of the assay.

Tenfold dilutions of cDNA reverse transcribed from 1010 copies of cRNA from each of the analyzed genes were separately coamplified with 20 ng of human gDNA, corresponding to ~104 copies of each gene. (A), melting curve peak areas for the 10-fold dilutions of MMP9 cDNA and the gDNA internal standard. (B), standard curve for MMP9 displaying the linear interval of 103–106 copies of MMP9 cDNA. The linear range of the assay was ~3 orders of magnitude for all the analyzed genes. The detection limit of the assay was 1000 copies of target cDNA, corresponding to a cDNA/gDNA ratio of 0.1 (red line), with 20 ng of gDNA as an internal standard. Above this amount, the CV was <20% for all of the analyzed genes.


Figure 5
View larger version (31K):
[in this window]
[in a new window]
  		Figure 5. Imprecision of the genome-controlled RT-PCR assay and the commercially available TaqMan gene-expression assay (Applied Biosystems).

A panel of 7 genes was analyzed in 8 repetitive experiments on 3 different colon carcinoma cell lines (COLO 205, DLD-1, and HT-29), and the imprecision of the genome-controlled RT-PCR assay (left panels) and the TaqMan assay (right panels) was compared. A cutoff of 0.1, corresponding to 1000 copies of cDNA, was used to define the detection limit of the genome-controlled RT-PCR assay. The TaqMan assay was multiplexed with GAPD to obtain optimal reproducibility. The correlation between the 2 different techniques was good, showing a discrepant result only in the expression levels of the MSH2 gene. The mean (error bars, SD) values of the cDNA/gDNA ratios are shown. The interassay CVs are displayed above the diagram for each of the analyzed genes. Values for genes expressed above the detection limit (genome-controlled RT-PCR) or above the lowest standard (TaqMan) are shown in boldface.

When we analyzed FFPE samples, we decreased the amount of internal standard to 2 ng of gDNA (containing ~103 copies of each gene). This decrease was necessary because the RNA yield from FFPE samples was lower than that from fresh biopsy samples. Thus, the linear interval of the assay was decreased to 102 to 105 copies of target cDNA. The detection limit was 100 copies of target cDNA, corresponding to a cDNA:gDNA ratio of 0.1. As a trade-off, the imprecision of the assay was slightly increased, with CVs ≤34% for all of the analyzed genes.

imprecision
For the genome-controlled RT-PCR technique, the intra- and interassay CVs were 5%–15% and 5%–20%, respectively. The intraassay CV of the TaqMan assay was 15%–45% with GAPD as an external standard and 10%–42% with GAPD multiplexing. The corresponding interassay CVs were 14%–43% and 7%–43%, respectively. The correlation between the different techniques was good; discrepant results occurred only for the expression amounts of the MSH2 gene. MSH2 expression was considerably higher when analyzed by the TaqMan technique, but the CV was high (43%), suggesting a technical problem in the TaqMan assay (Fig. 5Up ).

When we analyzed 8 replicate FFPE samples with 2 ng of gDNA as an internal standard, we observed expression above the detection limit (cDNA:gDNA ratio, >0.1) for GAPD, Try, MMP7, MMP9, uPA, and uPAR, whereas MMP2 and MSH2 were expressed below the detection limit. For these genes, the interassay variation of the genome-controlled RT-PCR assay, including deparaffinization and RNA extraction steps, was 11%–34%. When we analyzed the same samples with the TaqMan assay, expression of MMP9, uPAR, MMP2, and MSH2 was detected, with interassay CVs of 63%, 44%, 130%, and 110%, respectively. Expression amounts of MMP2 and MSH2 were below the lowest standard used in TaqMan assay, i.e., 100 copies of mRNA. For FFPE sample analysis, however, Try, MMP7, uPA, and GAPD did not amplify in the TaqMan assay, possibly because of the relatively long PCR amplicons for these genes (interval, 101–122 bp). In the genome-controlled RT-PCR assay the amplicon length of the expressed genes was 56–64 bp, whereas in the TaqMan PCR assay it was 54–87 bp.

analysis of biopsy samples and archival ffpe and frozen specimens
To demonstrate genome-controlled RT-PCR expression profiling, we analyzed fresh endoscopic biopsy samples from 8 patients undergoing follow-up for intestinal metaplasia of the lower esophagus (Barrett esophagus) (n = 27), and 1–5-year-old archival FFPE surgical specimens of colon carcinoma (n = 20) and healthy colonic mucosa (n = 15). The correlation between FFPE and frozen archival samples was tested in samples (n = 10) from 5 patients undergoing surgery for colon carcinoma.

We obtained sufficient amounts of mRNA from 26 of 27 single endoscopic biopsies. Distinct patterns of gene expression were observed in different anatomic regions: GAPD was consistently highly expressed in all regions. Samples of esophageal tissue displayed low expression of the rest of the analyzed genes; intestinal metaplasia tissue displayed consistent expression of uPAR and occasional expression of MMP-2, MMP-7, and uPA; and stomach tissue displayed relatively high expression of Try and uPAR (data not shown).

We obtained sufficient amounts of mRNA from 18 of 20 of the carcinoma FFPE tissue samples and from 12 of 15 healthy FFPE tissue samples. Carcinoma samples displayed variable expression of the analyzed genes, with significantly increased expression of MMP9 (P = 0.003) and uPA (P = 0.02) compared with healthy samples. Healthy colonic mucosa displayed high expression of GAPD and low expression of the rest of the analyzed genes.

We obtained sufficient amounts of mRNA from all 20 frozen carcinoma samples. We analyzed 15 colon cancer tissue specimens with similar morphology to test correlations between archival frozen tissue and FFPE specimens. The correlation between frozen and FFPE tumor tissue samples was good (Fig. 6 ).


Figure 6
View larger version (24K):
[in this window]
[in a new window]
  		Figure 6. Correlation between FFPE and frozen samples in 3 representative colon cancer tissue specimens.

The specimens were primarily divided for formalin fixation vs freezing macroscopically in the operating room. Specimens were stored for 5 years at room temperature (FFPE) and –80 °C (frozen tissue), respectively. RNA was extracted from 10 8-µm FFPE whole sections and 10 10-µm frozen sections, and tissue morphology was evaluated from an adjacent hematoxylin and eosin–stained section (FFPE samples) or toluidine blue–stained section (frozen samples) (not shown).


   Discussion
Top
 Abstract
 Introduction
 Materials and Methods
 Results
 Discussion
References
 
Genome-controlled RT-PCR allowed accurate measurement of the relative expression amounts of any mRNA species in small cell and tissue samples. This technique has several advantages over conventional RT-PCR. The assay is robust and reproducible because competitive coamplification of the gDNA internal standard occurs in every amplification reaction for any amount of target cDNA template. The detection limit and imprecision of the assay can be controlled by adjustment of the amount of internal standard. Differentiation between expressed and genomic sequences is achieved by modification of the cDNA transcript and does not rely on intronic sequences. This characteristic facilitates use of PCR amplicons short enough for analysis of FFPE samples and allows use of the untranslated regions for differential analysis of closely related genes. Carryover of gDNA from the RNA extraction step remains unmodified during reverse transcription and thus does not lead to false-positive results. Quantitative detection of the PCR products is done in-tube by melting-curve analysis of the PCR amplicons with a saturating DNA dye, such as SYBR Green or LC Green, lowering the costs of high-throughput applications.

The primer pairs for genome-controlled RT-PCR, in contrast to those of conventional RT-PCR (including real-time and competitive PCR), are designed to generate PCR amplicons within the boundaries of a single exon. The cDNA- and gDNA-derived amplicons are equal in length and differ only by a few nucleotide substitutions that produce a slight but sufficient difference in melting temperatures of the respective amplicons, ensuring truly competitive amplification of cDNA and gDNA templates with a single pair of primers. Short PCR amplicons are a prerequisite for melting-curve analysis of the PCR amplicons (28) and are also essential for analyzing mRNA in archival FFPE samples, in which most of the RNA is degraded (4). We used PCR amplicons of only 56–64 bp, a length that was efficient for amplifying the mRNA in FFPE samples. When we analyzed these samples with the TaqMan assay, we detected expression only for genes with <100 bp amplicons. With short PCR amplicons, melting temperature differences of ~4–5 °C yielded excellent discrimination power in the melting-curve analysis, with a dynamic range of ~3 orders of magnitude. The melting curve analyses were done on an ABI 7500 instrument, but we obtained comparable results with a Roche LightCycler. The 96-well format of these instruments allows PCRs to be performed on affordable benchtop thermal cyclers and uses the capacity of the detection instrument for melting curve analyses in high-throughput applications.

In the genome-controlled RT-PCR technique, the sequence, and thus the melting temperature, of the cDNA-derived amplicons is modified in the reverse transcription step. When the reactions are performed under nondenaturing conditions, any gDNA present in the reaction remains unmodified because the reverse transcriptase enzyme cannot use double-stranded DNA as a template (18)(19). In malignant tissues, the gDNA is often aneuploid, and the genetic aberrations vary from sample to sample. Therefore, we used as internal standard a fixed amount of diploid gDNA, in which the relative copy number of the analyzed genes was constant. In theory, carryover of aneuploid gDNA from the RNA isolation step could lead to analytical inaccuracy by skewing the relative copy number of the genes in the gDNA internal standard. This problem was eliminated by DNase digestion of the extracted RNA, which decreased DNA carry-over to a degreethat did not affect the quantitative analysis. For analysis of nonmalignant cells, RNA and DNA can be coisolated, and the gDNA contained in the analyzed cells can be used as an internal standard (18). In such a setup, the gene-specific cDNA:gDNA ratios will directly reflect the number of mRNA copies per genome, thus allowing for absolute quantification of mRNA production in the analyzed cells. This approach can also be used for analyzing changes in RNA production in cultured cells, e.g., after knockdown with short inhibitory RNAs, when possible aberrations in the genome of the analyzed cells remain constant.

In conclusion, genome-controlled RT-PCR enables highly reproducible and cost-efficient analysis of mRNA production in small samples, such as needle and endoscopic biopsy specimens and FFPE tissue specimens. This technique provides an efficient research tool for large-scale analysis of mRNA profiles in archival FFPE samples, as well as for monitoring changes in entire gene-expression pathways and networks in cell cultures. Because assay imprecision can be controlled by adjusting the amount of internal standard, in the future this technique may be used for diagnostic gene-expression analysis.


   Acknowledgments
 
This study was supported by grants from Finska Läkaresällskapet.


   Footnotes
 
3 These authors contributed equally to this work.

1 Nonstandard abbreviations: RT-PCR, reverse transcription–PCR; FFPE, formalin-fixed paraffin-embedded; gDNA, genomic DNA.

2 Human genes: Try, trypsinogen; MMP2/7/9, matrix-metalloproteinase-2/-7/-9; uPA, urokinase-type plasminogen activator; uPAR, urokinase-type plasminogen activator receptor; MSH2, Mut S homolog 2; GAPD, glyceraldehyde-3-phosphate dehydrogenase.


   References
Top
 Abstract
 Introduction
 Materials and Methods
 Results
 Discussion
 References
 

  1. Golub TR, Slonim DK, Tamayo P, Huard C, Gaasenbeek M, Mesirov JP, et al. Molecular classification of cancer: class discovery and class prediction by gene expression monitoring. Science 1999;286:531-537.[Abstract/Free Full Text]
  2. Emmert-Buck MR, Strausberg RL, Krizman DB, Bonaldo MF, Bonner RF, Bostwick DG, et al. Molecular profiling of clinical tissues specimens: feasibility and applications. J Mol Diagn 2000;2:60-66.[Free Full Text]
  3. Alizadeh AA, Ross DT, Perou CM, van de Rijn M. Towards a novel classification of human malignancies based on gene expression patterns. J Pathol 2001;195:41-52.[CrossRef][ISI][Medline] [Order article via Infotrieve]
  4. Cronin M, Pho M, Dutta D, Stephans JC, Shak S, Kiefer MC, et al. Measurement of gene expression in archival paraffin-embedded tissues: development and performance of a 92-gene reverse transcriptase-polymerase chain reaction assay. Am J Pathol 2004;164:35-42.[Abstract/Free Full Text]
  5. Bertucci F, Salas S, Eysteries S, Nasser V, Finetti P, Ginestier C, et al. Gene expression profiling of colon cancer by DNA microarrays and correlation with histoclinical parameters. Oncogene 2004;23:1377-1391.[CrossRef][ISI][Medline] [Order article via Infotrieve]
  6. Lossos IS, Czerwinski DK, Alizadeh AA, Wechser MA, Tibshirani R, Botstein D, et al. Prediction of survival in diffuse large-B-cell lymphoma based on the expression of six genes. N Engl J Med 2004;350:1828-1837.[Abstract/Free Full Text]
  7. Bustin SA. Absolute quantification of mRNA using real-time reverse transcription polymerase chain reaction assays. J Mol Endocrinol 2000;25:169-193.[Abstract]
  8. Bernard PS, Wittwer CT. Real-time PCR technology for cancer diagnostics. Clin Chem 2002;48:1178-1185.[Abstract/Free Full Text]
  9. Wang AM, Doyle MV, Mark DF. Quantitation of mRNA by the polymerase chain reaction. Proc Natl Acad Sci U S A 1989;86:9717-9721.[Abstract/Free Full Text]
  10. Becker-Andre M, Hahlbrock K. Absolute mRNA quantification using the polymerase chain reaction (PCR). A novel approach by a PCR aided transcript titration assay (PATTY). Nucleic Acids Res 1989;17:9437-9446.[Abstract/Free Full Text]
  11. Gilliland G, Perrin S, Blanchard K, Bunn HF. Analysis of cytokine mRNA and DNA: detection and quantitation by competitive polymerase chain reaction. Proc Natl Acad Sci U S A 1990;87:2725-2729.[Abstract/Free Full Text]
  12. Siebert PD, Larrick JW. PCR MIMICS: competitive DNA fragments for use as internal standards in quantitative PCR. Biotechniques 1993;14:244-249.[ISI][Medline] [Order article via Infotrieve]
  13. Holland PM, Abramson RD, Watson R, Gelfand DH. Detection of specific polymerase chain reaction product by utilizing the 5'—-3' exonuclease activity of Thermus aquaticus DNA polymerase. Proc Natl Acad Sci U S A 1991;88:7276-7280.[Abstract/Free Full Text]
  14. Gibson UE, Heid CA, Williams PM. A novel method for real time quantitative RT-PCR. Genome Res 1996;6:995-1001.[Abstract/Free Full Text]
  15. Ding C, Cantor CR. A high-throughput gene expression analysis technique using competitive PCR and matrix-assisted laser desorption ionization time-of-flight MS. Proc Natl Acad Sci U S A 2003;100:3059-3064.[Abstract/Free Full Text]
  16. Stenman J, Finne P, Stahls A, Grenman R, Stenman UH, Palotie A, et al. Accurate determination of relative messenger RNA levels by RT-PCR. Nat Biotechnol 1999;17:720-722.[CrossRef][ISI][Medline] [Order article via Infotrieve]
  17. Stenman J, Orpana A. Accuracy in amplification. Nat Biotechnol 2001;19:1011-1012.[ISI][Medline] [Order article via Infotrieve]
  18. Zenilman ME, Graham W, Tanner K, Shuldiner AR. Competitive reverse-transcriptase polymerase chain reaction without an artificial internal standard. Anal Biochem 1995;224:339-346.[CrossRef][ISI][Medline] [Order article via Infotrieve]
  19. Joo CH, Lee H, Kim E, Lee B, Cho YK, Kim YK. Differential amplifying RT-PCR: a novel RT-PCR method to differentiate mRNA from its DNA lacking intron. J Virol Methods 2002;100:71-81.[CrossRef][ISI][Medline] [Order article via Infotrieve]
  20. Stenman J, Paju A, Rissanen O, Orpana A. Method for quantitative and/or comparative measurement of mRNA expression levels in small biological samples. European Patent Register WO2005116248; May 25, 2005..
  21. Zeng ZS, Huang Y, Cohen AM, Guillem JG. Prediction of colorectal cancer relapse and survival via tissue RNA levels of matrix metalloproteinase-9. J Clin Oncol 1996;14:3133-3140.[Abstract]
  22. Yang JL, Seetoo D, Wang Y, Ranson M, Berney CR, Ham JM, et al. Urokinase-type plasminogen activator and its receptor in colorectal cancer: independent prognostic factors of metastasis and cancer-specific survival and potential therapeutic targets. Int J Cancer 2000;89:431-439.[CrossRef][ISI][Medline] [Order article via Infotrieve]
  23. Yamamoto H, Iku S, Adachi Y, Imsumran A, Taniguchi H, Nosho K, et al. Association of trypsin expression with tumour progression and matrilysin expression in human colorectal cancer. J Pathol 2003;199:176-184.[CrossRef][ISI][Medline] [Order article via Infotrieve]
  24. Parc Y, Gueroult S, Mourra N, Serfaty L, Flejou JF, Tiret E, et al. Prognostic significance of microsatellite instability determined by immunohistochemical staining of MSH2 and MLH1 in sporadic T3N0M0 colon cancer. Gut 2004;53:371-375.[Abstract/Free Full Text]
  25. Danenberg K, Danenberg P, Swenson S. Isolation of RNA, DNA, and proteins from formalin-fixed paraffin-embedded tissue specimens. US Patent US2006199197; May 7, 2006..
  26. Chomczynski P, Sacchi N. Single-step method of RNA isolation by acid guanidinium thiocyanate-phenol-chloroform extraction. Anal Biochem 1987;162:156-159.[ISI][Medline] [Order article via Infotrieve]
  27. Koshkin AA, Singh SK, Nielsen P, Rajwanshi VK, Kumar R, Meldgaard M, et al. LNA (locked nucleic acids): synthesis of the adenine, cytosine, guanine, 5-methylcytosine, thymine and uracil bicyclonucleoside monomers, oligomerisation, and unprecedented nucleic acid recognition. Tetrahedron 1998;54:3607-3630.[CrossRef][ISI]
  28. Wittwer CT, Reed GH, Gundry CN, Vandersteen JG, Pryor RJ. High-resolution genotyping by amplicon melting analysis using LCGreen. Clin Chem 2003;49:853-860.[Abstract/Free Full Text]




This Article
Right arrow Abstract Freely available
Right arrow Full Text (PDF)
Right arrow 073064.Supplemental Data Table 1
Right arrow All Versions of this Article:
clinchem.2006.073064v1
52/11/1988    most recent
Right arrow Submit an electronic Letter to
the Editor about this paper
Right arrow Alert me when this article is cited
Right arrow Alert me when eLetters are posted
Right arrow Alert me if a correction is posted
Right arrow Citation Map
Services
Right arrow Email this article to a friend
Right arrow Similar articles in this journal
Right arrow Similar articles in ISI Web of Science
Right arrow Similar articles in PubMed
Right arrow Alert me to new issues of the journal
Right arrow Download to citation manager
Right arrow reprints & permissions
Citing Articles
Right arrow Citing Articles via ISI Web of Science (1)
Right arrow Citing Articles via Google Scholar
Google Scholar
Right arrow Articles by Stenman, J.
Right arrow Articles by Orpana, A.
Right arrow Search for Related Content
PubMed
Right arrow PubMed Citation
Right arrow Articles by Stenman, J.
Right arrow Articles by Orpana, A.
Related Collections
Right arrow Molecular Diagnostics and Genetics

HOME HELP FEEDBACK SUBSCRIPTIONS ARCHIVE SEARCH TABLE OF CONTENTS