|
|
||||||||
1 Unité de Pharmacologie, Institut Curie, 26 rue dUlm, 75005 Paris, France
2 Laboratoire INSERM E0017, Centre René Huguenin, Saint Cloud, France
3 Laboratoire de Génétique Moléculaire, UPRES EA 3618, Faculté des Sciences Pharmaceutiques et Biologiques de Paris, Paris, France
4 Laboratoire dOncobiologie, Centre René Huguenin, Saint Cloud, France
5 Département dAnatomo-Pathologie, Centre René Huguenin, Saint Cloud, France
6 Service dAnatomo-Pathologie, Institut Curie, Paris, France
7 Service de Biostatistiques, Institut Curie, Paris, France
(Requests for offprints should be addressed to P de Cremoux; Email: patricia.de-cremoux{at}curie.net)
| Abstract |
|---|
|
|
|---|
We performed inter-laboratory quality control of quantitative PCR between two independent laboratories, using in-house RT-PCR assays on a series of hormone-related target genes in a retrospective consecutive series of 79 breast tumors. Total RNA was reverse transcribed in a single center. Calibration curves were performed for five target genes (estrogen receptor (ER)
, ERß, progesterone receptor (PR), CYP19 (aromatase) and Ki 67) and for two reference genes (human acidic ribosomal phosphoprotein PO (RPLPO) and TATA box-binding protein (TBP)). Amplification efficiencies of the calibrator were determined for each run and used to calculate mRNA expression. Correlation coefficients were evaluated for each target and each reference gene.
A good correlation was observed for all target and reference genes in both centers using their own protocols and kits (P < 0.0001). The correlation coefficients ranged from 0.90 to 0.98 for the various target genes in the two centers. A good correlation was observed between the level of expression of the ER
and the PR transcripts (P < 0.001). A weak inverse correlation was observed in both centers between ER
and ERß levels, but only when TBP was the reference gene. No other correlation was observed with other parameters.
Real-time PCR assays allow convenient quantification of target mRNA transcripts and quantification of target-derived nucleic acids in clinical specimens. This study addresses the importance of inter-laboratory quality controls for the use of a panel of real-time PCR assays devoted to clinical samples and protocols and to ensure their appropriate accuracy. This can also facilitate exchanges and multicenter comparison of data.
| Introduction |
|---|
|
|
|---|
An analysis of PCR techniques shows that false-positive results due to carry-over and false-negative results due to PCR inhibition are not systematically controlled. These shortcomings of classical PCR should be resolved with the development of standardized real-time PCR assays. Inter-laboratory comparisons should then be possible and regular quality controls will be necessary to ensure the reliability of real-time PCR assays.
Inter-laboratory quality assurance programs for diagnostic tests based on nucleic acid amplification are not widely implemented in clinical laboratories and remain limited to a few tests, mainly in hematology, microbiology and virology. The development of specific inter-laboratory quality assurance trials based on methodologic proficiency testing and directed to the evaluation of any diagnostic molecular target may be useful, particularly in breast cancer patients treated by hormone therapy. Endocrine manipulations are the most effective and least toxic form of systemic therapy currently available for the management of hormone-dependent breast cancer. The panel of available compounds has grown considerably over recent years, with the addition to the older agents, tamoxifen and high-dose progestogens, of type I and II aromatase inhibitors, pure anti-estrogens and, for premenopausal women, luteinizing hormone-releasing hormone analogs. The question of the determinants of resistance or sensitivity to treatment therefore needs to be considered. The estrogen receptor (ER) assay was first proposed by Jensen et al. (1971) to determine which patients are likely to respond to endocrine therapy. It was subsequently established that patients with ER-rich tumors respond to endocrine therapy, whereas patients with ER-negative tumors are unlikely to respond (McGuire et al. 1975). Although response rates of up to 70% are achieved with endocrine therapy in patients with tumors expressing ER and progesterone receptor (PR), most responsive tumors eventually acquire resistance. The mechanisms responsible for the resistance of tumors expressing ER and/or PR may be multifactorial. Molecular biotechnology offers a great potential for developing multiparametric analysis of mRNA expression of clinically relevant biological parameters in estrogen-dependent breast tumors.
We first evaluated endogenous controls for gene expression studies in breast cancer (F Spyratos, unpublished observations). The present study was designed to evaluate the accuracy of quantitative analyses of endocrine-related gene expression in breast cancer using real-time RT-PCR so that they can be used in prospective multicenter clinical trials. Quantitative RT-PCR assays may be particularly useful for quantitative determination of weakly expressed targets such as CYP19 (aromatase) and ERß mRNA for which no consensus immunohistochemistry assay is available.
This quality control was performed between two independent laboratories (Centre René Huguenin (CRH) and Institut Curie (IC)) on a series of 79 breast cancer samples. The endocrine-related target genes were in-house assays designed by each laboratory (ER
, ERß, PR and CYP19). Similar Ki 67 and references gene (human acidic ribosomal phosphoprotein PO (RPLPO) and TATA box-binding protein (TBP)) designs were used in both centers.
| Materials and methods |
|---|
|
|
|---|
A panel of 79 untreated residual breast tissue samples from informed patients undergoing primary surgery for breast cancer were retrospectively analyzed: 39 were provided by the CRH and 40 were provided by the IC. Tumor sampling was performed extemporaneously by a pathologist. Samples were frozen in liquid nitrogen within 20 min after resection using the same procedure in each center. The residual fragment was stored in liquid nitrogen in RNase-free conditions until use. We retrospectively analyzed this residual tissue for ER
, ERß, PR, CYP19 and Ki 67 mRNA expression. Adjacent fragments of frozen material were fixed in alcohol-formalin-acetic acid (AFA) and stained with hematoxylin and eosin in order to determine the proportion of tumor cells in each sample. Clinical staging and histological diagnosis were routinely determined in each case.
Calibration curves
Calibration curves were determined for each target gene using a unique mix of relevant total RNA including: a pool of total RNA from breast tumors known to have high levels of ER
transcripts and proteins, a pool of total RNA from breast tumors known to have high levels of PR transcripts and proteins (Bieche et al. 2001, de Cremoux et al. 2002), an ovarian total RNA, known to express high levels of ERß transcripts (Moore et al. 1998, Pujol et al. 1998) (human ovary total RNA reference 735260 (Stratagene, Amsterdam, The Netherlands)) and a placental total RNA (Moore et al. 1998) (human placenta total RNA reference 64024-1 (Clontech, USA)) known to express high levels of CYP19 gene transcripts. Six-fold diluted cDNA (1:25 to 1:800 for ERß, PR, Ki 67, TBP and RPLPO; 1:25 to 1:3200 for ER
and CYP19) were obtained after reverse transcription of the mix of these relevant total RNA and were amplified to obtain a standard curve and to assess PCR efficiency. The highest level of this calibration curve was used as calibrator in each experiment for each target gene.
RNA extraction
Total RNA was extracted from crushed tumor samples by the RNA plus kit (Bioprobe, Saint Quentin, France) at IC and by RNABle (Eurobio, Les Ulis, France) at CRH, according to the manufacturers instructions. The RNA samples were stored in RNase-free distilled water, at 80 °C. The quality of RNA samples was determined by electrophoresis on 2% agarose gel with ethidium bromide staining. The 18S and 28S RNA bands were visualized under u.v. light. All samples presented a high quality RNA profile. The quantity of RNA was measured by spectrophotometry at 260 nm.
cDNA synthesis
RNA was reverse transcribed in one center (CRH) to avoid possible differences related to variable efficiency of reverse transcriptases. Total RNA was reverse transcribed in a final volume of 20 µl containing 1x reverse transcriptase buffer (500 mM of each deoxynucleotide triphosphate, 3 mM MgCl2, 75 mM KCl and 50 mM TrisHCl (pH 8.3)), 10 units RNase inhibitor (Promega, Madison, WI, USA), 10 mM dithiothreitol, 50 units Superscript II RNase Hreverse transcriptase (Invitrogen), 1.5 mM random hexamers (Pharmacia, Uppsala, Sweden) and 1 µg total RNA (patient samples and calibrator). Samples were incubated at 20 °C for 10 min and 42 °C for 30 min, reverse transcriptase was then inactivated by heating at 99 °C for 5 min and cooling for 5 min.
Real-time PCR amplification
ER
, ERß, PR, CYP19 and Ki 67 transcripts were quantified using real-time quantitative RT-PCR assays. Primers and probes were selected by Oligo 4.0 (National Biosciences, Plymouth, MN, USA) (CRH) or Primer Express (Applied Biosystems, Foster City, CA, USA) (IC) software. The nucleotide sequences were then blasted against dbEST, htgs and nr (the non-redundant set of the GenBank, European Molecular Biology Laboratory (EMBL) and DNA Data Bank of Japan (DDBJ) database sequences) to confirm the total gene specificity of the nucleotide sequences chosen as primers and probes. To avoid amplification of contaminating genomic DNA, one of the two primers was placed in a different exon.
Each laboratory used its own design for ER
, ERß, PR and CYP19 (Bieche et al. 2001, de Cremoux et al. 2002) gene transcript analysis. Common primers and probes designed by CRH were used for Ki 67, TBP and RPLPO (Bieche et al. 2000, Spyratos et al. 2002). The nucleotide and probe sequences are shown in Table 1
.
|
All measurements were performed in duplicate under blind conditions. Two non-template controls were included in each amplification run.
Results are expressed as n-fold differences in target gene expression relative to reference gene (either TBP and RPLPO) and the calibrator was defined as the N target. It was determined as follows: Ntarget = 1+[E(CtcalibratorCtsample)target/E(CtcalibratorCtsample)referencegene], where E is the efficiency of PCR measured using the slope of the calibration curve and Ct is the cycle threshold.
Statistical analysis
Data are expressed as ratios between the quantity of target transcripts and the quantity of the reference genes TBP and RPLPO, in arbitrary units. Non-parametric tests were used. Spearmans rank correlation was used to evaluate the relationship between relative mRNA expression of target genes at IC and CRH, relative to both reference genes.
| Results |
|---|
|
|
|---|
and PR protein content were initially measured by enzyme immunoassay or immunohistochemistry (IHC). The distribution of ER/PR content was as follows: 43 (54%) samples were ER positive and PR positive, 14 (18%) samples were ER negative and PR negative, 17 (22%) samples were ER positive and PR negative and five (6%) samples were ER negative and PR positive.
|
PCR optimization had been performed previously in each laboratory using similar rules as described (Bieche et al. 2001, Spyratos et al. 2002, de Cremoux et al. 2003).
The specificity of all PCR products was also initially determined in both centers by the presence of a single band of the expected size after migration on agarose gels.
The linearity for each target gene and reference gene was observed for more than 5 log dilutions, including the observed samples values.
A calibration curve was generated by analyzing serial dilutions of the pool of RNA, as described in Materials and methods for each target and reference gene, to assess PCR efficiency. The PCR efficiency for ER
, ERß, PR, CYP19, Ki 67, RPLPO and TBP was above 90% for all genes. The intra-assay coefficient of variation (CV) of Ct values was <2% and the interassay CV was <5% as calculated by the standard error.
Correlation of target gene mRNA transcript levels between the two centers for each reference gene
All results are expressed as continuous variables. For each parameter, and for both reference genes, the mean and standard deviation are presented in Table 3
for TBP and for RPLPO as reference genes. The levels of expression were comparable for all target genes in both centers. As shown in Tables 4
and 5
, a statistically significant correlation (P < 0.0001) was demonstrated between the level of gene expression between the two centers for all target genes, considering either TBP or RPLPO as the reference gene.
|
|
|
We observed a positive correlation between ER
and PR transcript levels in both centers and with both reference genes, with correlation coefficients of 0.49 and 0.47 at IC, and 0.41 and 0.42 at CRH, using TBP and RPLPO (P < 0.001) respectively. No correlation was observed in this limited series between ER
and Ki 67 and aromatase transcripts. A weak inverse correlation was observed between ER
and ERß only with TBP as reference gene at IC (r = 0.33) and CRH (r = 0.40) (P < 0.003 and P < 0.0004 respectively).
Interpretation of discrepancies
The expression of all target genes was concordant for all samples but one: for example, samples expressing high levels of ER
and PR transcripts in one center expressed high levels of ER
and PR mRNA transcripts in the other center, except in one case where a 5 Ct difference was observed between the ER
measurement performed at CRH and IC. This sample clearly expressed low levels of ER
mRNA at IC and high levels of ER
mRNA at CRH, using both TBP and RPLPO as reference gene. Sequencing of the two ER
amplicons demonstrated that this sample presented a rare polymorphism in 1559 nt (gtg > ggg) in position 10 of the IC upper primer.
In conclusion, no other major discrepancy was observed in any of the target or reference genes.
| Discussion |
|---|
|
|
|---|
For mRNA detection, we used a quantitative RT-PCR method based on real-time analysis of PCR amplification. This assay is sensitive and allows accurate quantification of mRNA expression, as previously demonstrated for ER
and PR transcripts by both laboratories (CRH and IC) (Bieche et al. 2001, de Cremoux et al. 2002, 2003, Spyratos et al. 2002). The real-time PCR method has several advantages over classical RT-PCR-based methods. It does not require post-PCR sample handling, thereby avoiding problems related to carry-over. It also has a wide dynamic range with very low intra-assay and interassay CV values (<5%).
Real-time RT-PCR specific errors in the quantification of mRNA transcripts are usually correlated to possible variations in starting material which may be a major problem for patient samples for which sampling is uncertain (Freeman et al. 1999). To compensate for sample-to-sample variations of mRNA content, each target gene was normalized with endogenous reference genes (RPLPO, TBP). The RT step also largely contributes to the overall variation of quantitative RT-PCR (Bustin 2000, Ovstebo et al. 2003). However, in this study, we used a common cDNA sample to avoid variations due to this step. Consequently, only the PCR step was evaluated in this quality control.
To compensate for any possible inter-PCR variation between runs, the expression of each target gene was calibrated. In these assays, designed to be used in clinical trials, the standard curve for the calibrator was prepared as a mix of human breast tumor samples and two marketed ovarian and placental RNAs known to present high levels of ERß and CYP19 mRNA. The mix was used in all experiments, allowing us to control the reliability of the method over time. A calibration curve was included in each experiment and analyzed under the same conditions as for unknown samples. PCR efficiency was also calculated from the calibration curve in each experiment and for each transcript. The possible variation from one run to the next was therefore taken into account (as described in Materials and Methods). The use of a calibration curve for quantitative gene expression studies is not widely applied, but has been previously described to be a useful tool (Ovstebo et al. 2003).
Real-time RT-PCR analyses of tumor mRNA in both centers revealed a good correlation for each target gene and each control gene. We used different ER
, ERß, PR and CYP19 specific primer sets and probes and, as expected, we did not observe any differences due to the localization of the amplicon analyzed, except in one case for ER
, corresponding to the presence of a rare polymorphism. The real-time PCR procedure and the PCR kit used did not influence the quantification of mRNA. The amplification conditions of the PCR assay described here were previously optimized and validated by each center (Bieche et al. 2001, de Cremoux et al. 2002). This collaborative study showed a high reproducibility of the PCR assay between two independent centers when identical cDNA were analyzed regardless of the PCR kit and design used.
Furthermore, this good correlation was not only found for genes with high mRNA expression, such as ER
, PR, Ki 67, but also for genes with low mRNA expression, such as CYP19 gene transcripts and ERß. This represents an additional relevant tool for analysis of hormone-dependent gene expression assays in clinical samples, where no other routinely adapted tool is currently available.
Multicenter validation and standardization of PCR have been previously described for the detection of food-borne bacteria transferred from animal to food products, by semiquantitative PCR, using a unique primer design for five pathogens (Hoorfar & Cook 2003, Malorny et al. 2003). The present study represents the first real-time PCR quality control in breast cancer.
In another field, a multicenter European prospective study of the candidate references genes for diagnosis and residual disease detection in leukemic patients by quantitative RT-PCR was recently published (Beillard et al. 2003). This study initially selected three reference genes from a total of 14 (Abelson, ß2 microglobulin and ß glucuronidase) and prospective and retrospective studies on more than 400 samples confirmed that these three genes had a stable mRNA expression. However, the authors of this study propose that only the Abelson gene should be used as a reference gene for quantitative RT-PCR analysis of residual disease, as its expression is not significantly different between normal and leukemic samples at diagnosis, confirming the major importance of these quality control studies before using these tools in prospective clinical trials. The results obtained in this study on snap-frozen material allowed us to use our protocol in a recent multicenter prospective clinical trial of neoadjuvant endocrine therapy whose preliminary results were presented in the last San Antonio Breast Cancer Symposium (Tubiana-Hulin et al. 2003).
In conclusion, we have developed the first independent quality control of mRNA expression of estrogen-dependent target genes in breast carcinoma by quantitative RT-PCR. The good correlation between two independent centers using a different design of PCR analysis allows these assays to be used for prospective clinical trials.
| Acknowledgements |
|---|
| References |
|---|
|
|
|---|
Bieche I, Nogues C, Paradis V, Olivi M, Bedossa P, Lidereau R & Vidaud M 2000 Quantitation of hTERT gene expression in sporadic breast tumors with a real-time reverse transcriptionpolymerase chain reaction assay. Clinical Cancer Research 6 452459.
Bieche I, Parfait BIL, Girault I, Vidaud M & Lidereau R 2001 Quanti.cation of estrogen receptor alpha and beta expression in sporadic breast cancer. Oncogene 20 81098115.[CrossRef][ISI][Medline]
Bustin S 2000 Absolute quanti.cation of mRNA using real-time reverse transcription polymerase chain reaction assays. Journal of Molecular Endocrinology 25 169193.[Abstract]
de Cremoux P, Tran-Perennou C, Elie C, Boudou E, Barbaroux C, Poupon MF, De Rycke Y, Asselain B & Magdelenat H 2002 Quantitation of estradiol receptors alpha and beta and progesterone receptors in human breast tumors by real-time reverse transcription-polymerase chain reaction. Correlation with protein assays. Biochemical Pharmacology 64 507515.[CrossRef][Medline]
de Cremoux P, Tran-Perennou C, Brockdorff B, Boudou E, Brunner N, Magdelenat H & Lykkesfeldt A 2003 Validation of real-time RT-PCR for analysis of human breast cancer cell lines resistant or sensitive to treatment with antiestrogens. Endocrine-Related Cancer 10 409418.[Abstract]
Freeman WM, Walker SJ & Vrana KE 1999 Quantitative RTPCR: pitfalls and potential. BioTechniques 26 112122, 124125.[ISI][Medline]
Hoorfar J & Cook N 2003 Critical aspects of standardization of PCR. Methods in Molecular Biology 216 5164.
Jensen EV, Block GE, Smith S, Kyser K & DeSombre ER 1971 Estrogen receptors and breast cancer response to adrenalectomy. National Cancer Institute Monographs 34 5570.
McGuire WL, Carbone PP, Sears ME & Escher CG 1975 Estrogen Receptors in Human Breast Cancer. An Overview. Eds McGuire WL, Carbone PP & Vollner EP. New York: Raven Press.
Malorny B, Hoorfar J, Bunge C & Helmuth R 2003 Multicenter validation of the analytical accuracy of Salmonella PCR: towards an international standard. Applied Environmental Microbiology 69 290296.
Moore J, McKee D, Slentz-Kesler K, Moore L, Jones S, Horne E, Su J, Kliewer S, Lehmann J & Willson T 1998 Cloning and characterization of human estrogen receptor beta isoforms. Biochemical and Biophysical Research Communications 247 7578.[CrossRef][ISI][Medline]
Ovstebo R, Haug K, Lande K & Kierulf P 2003 PCR-based calibration curves for studies of quantitative gene expression in human monocytes: development and evaluation. Clinical Chemistry 49 425432.
Pujol P, Rey J, Nirde P, Roger P, Gastaldi M, Laffargue F, Rochefort H & Maudelonde T 1998 Differential expression of estrogen receptor-alpha and -beta messenger RNAs as a potential marker of ovarian carcinogenesis. Cancer Research 58 53675373.
Spyratos F, Bouchet C, Tozlu S, Labroquere M, Vignaud S, Becette V, Lidereau R & Bieche I 2002 Prognostic value of uPA, PAI-1 and PAI-2 mRNA expression in primary breast cancer. Anticancer Research 22 29973003.[Medline]
Tubiana-Hulin M, Spyratos F, Becette V, Mauriac L, Romieu G, Bibeau F, Bieche I, Bourgeois H & Chollet P 2003 Phase II study of neo-adjuvant exemestane in postmenopausal patients with operable breast cancer (26th Annual San Antonio Breast Cancer Symposium, December 36). Breast Cancer Research Treatment 82 (Suppl 1) A443.
This article has been cited by other articles:
![]() |
P. de Cremoux, D. Rosenberg, J. Goussard, C. Bremont-Weil, F. Tissier, C. Tran-Perennou, L. Groussin, X. Bertagna, J. Bertherat, and M.-L. Raffin-Sanson Expression of progesterone and estradiol receptors in normal adrenal cortex, adrenocortical tumors, and primary pigmented nodular adrenocortical disease Endocr. Relat. Cancer, June 1, 2008; 15(2): 465 - 474. [Abstract] [Full Text] [PDF] |
||||
![]() |
R. S. Sane, D. J. Buckley, A. R. Buckley, S. C. Nallani, and P. B. Desai Role of Human Pregnane X Receptor in Tamoxifen- and 4-Hydroxytamoxifen-Mediated CYP3A4 Induction in Primary Human Hepatocytes and LS174T Cells Drug Metab. Dispos., May 1, 2008; 36(5): 946 - 954. [Abstract] [Full Text] [PDF] |
||||
![]() |
E. Marangoni, A. Vincent-Salomon, N. Auger, A. Degeorges, F. Assayag, P. de Cremoux, L. de Plater, C. Guyader, G. De Pinieux, J.-G. Judde, et al. A New Model of Patient Tumor-Derived Breast Cancer Xenografts for Preclinical Assays Clin. Cancer Res., July 1, 2007; 13(13): 3989 - 3998. [Abstract] [Full Text] [PDF] |
||||
![]() |
A. B. Dydensborg, E. Herring, J. Auclair, E. Tremblay, and J.-F. Beaulieu Normalizing genes for quantitative RT-PCR in differentiating human intestinal epithelial cells and adenocarcinomas of the colon Am J Physiol Gastrointest Liver Physiol, May 1, 2006; 290(5): G1067 - G1074. [Abstract] [Full Text] [PDF] |
||||
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |