| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |
|
|
|
|||||||
|
|||||||||
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
INNOVATIVE METHODOLOGY
Canadian Institutes of Health Research Group in Functional Development and Physiopathology of the Digestive Tract, Département d'Anatomie et de Biologie Cellulaire, Faculté de Médecine, Université de Sherbrooke, Sherbrooke, Québec, Canada
Submitted 23 May 2005 ; accepted in final form 26 December 2005
 |
ABSTRACT |
|---|
 TOP ABSTRACT MATERIALS AND METHODSRESULTSDISCUSSIONGRANTSREFERENCES |
|---|
-actin and glyceraldehyde-3-phosphate dehydrogenase. It has, however, become clear that these and other normalizing genes can display modulated patterns of expression across tissue types and during complex cellular processes such as cell differentiation and cancer progression. Our objective was to set the basis for identifying normalizing genes that displayed stable expression during enterocytic differentiation and between healthy tissue and adenocarcinomas of the human colon. We thus identified novel potential normalizing genes using previously generated cDNA microarray data and examined the alterations of expression of two of these genes as well as seven commonly used normalizing genes during the enterocytic differentiation process and between matched pairs of resection margins and primary carcinomas of the human colon using real-time RT-PCR. We found that ribosomal phosphoprotein P0 was particularly stable in all intestinal epithelial cell extracts, thereby representing a particularly robust housekeeping reference gene for the assessment of gene expression during the human enterocytic differentiation process. On the other hand, -2-microglobulin generated the best score as a normalizing gene for comparing human colon primary carcinomas with their corresponding normal mucosa of the resection margin, although others were found to represent acceptable alternatives. In conclusion, we identified and characterized specific normalizing genes that should significantly improve quantitative mRNA studies related to both the differentiation process of the human intestinal epithelium and adenocarcinomas of the human colon. This approach should also be useful to validate normalizing genes in other intestinal contexts.
mRNA analysis; intestinal cell differentiation; colon cancer cells
The differentiation of a tissue (33) or the formation of carcinomas (7) are complex cellular processes that contribute to the difficulty of assessing gene expression changes using RT-PCR. Indeed, inherent difficulties in identifying a "gold standard" normalizing gene for cDNA input in RT-PCR analysis experiments have become apparent from reports of variation of expression across tissue types (26), developmental stages (9), and cancer types (1, 11). One proposed alternative for addressing this problem is to use a weighted average of several normalizing genes (24, 32). However, this option suffers from the drawback that it only normalizes a gene of interest against the average variation of many normalizing genes. While the ubiquitous normalizing gene may not exist, another more promising strategy is to identify a single normalizing gene for a single tissue type (27), pathology (12), or experimental design (5, 16, 18, 28).
The aim of the present study was to identify and validate such potential normalizing genes to be used for studies of proliferating and differentiating human intestinal epithelial cells as well as human colon adenocarcinomas. In the present study, we identified novel candidate normalizing genes using cDNA microarray data from human intestinal epithelial cell models reflecting different stages of enterocytic differentiation and selected two of these for validation along with seven commonly used normalizing genes. We compared the variation of expression of these genes by validating their stability during enterocytic differentiation using real-time qRT-PCR on stringently purified and quantified RNA samples. We found that the ribosomal phosphoprotein P0 (RPLP0) gene (10) is the most suited for normalization purposes during the enterocytic differentiation process. Among the nine chosen candidate genes further tested for variation of expression between adenocarcinomas and their corresponding histologically normal resection margins, we found that the -2-microgobulin (B2M) gene is best for normalizing gene expression, although other genes were found to represent acceptable alternatives.
|
MATERIALS AND METHODS |
|---|
|
TOPABSTRACTMATERIALS AND METHODSRESULTSDISCUSSIONGRANTSREFERENCES |
|---|
Tissues. Specimens of the ileum from five fetuses ranging from 16 to 20 wk of age (postfertilization) were obtained after legal therapeutic abortion. Only specimens obtained rapidly (within 60 min) were used. The project was in accordance with the protocol approved by the Institutional Human Research Review Committee at the Université de Sherbrooke for the use of human material. Fetal small intestinal epithelia were separated from the mesenchyme using Matrisperse (Collaborative Biomedical Products, Becton-Dickenson Labware; Mississauga, ON, Canada), and fraction purity was confirmed by RT-PCR as described previously (21, 23).
Cell culture. Two epithelial cell lines, which together with the samples obtained from tissues recapitulate the crypt-villus axis of the fetal intestine (20), were used in this study. The HIEC-6 cell line was generated from the normal fetal human intestine (22). These cells express a number of crypt cell markers but no villus cell markers, appear to be unable to differentiate, and are thus considered to be intestinal stemlike cells (20, 22). HIEC cells were used at passages 14–17 and grown as described (22). The Caco-2/15 cell line, a stable clone of the parent Caco-2 cell line (25), has been characterized elsewhere (4, 31). Although of cancerous origin, these cells are unique in that upon confluence, they spontaneously undergo a gradual villus-like enterocytic differentiation process, similar to that observed in the epithelium of the intact fetal small and large intestine (19, 20). Cells between passages 14 and 22, used for the generation of cDNA microarrays, were cultured in plastic dishes as previously described (31) and analyzed at 0–25 days postconfluence.
Three other human cell lines were used to construct a reference sample for the cDNA microarrays: 1) the HT-29 cell line, a cell line derived from a carcinoma of the colon (3); 2) the A549 cell line, a pulmonary adenocarcinoma cell line (30) obtained from Dr. Jacques Bérard; and 3) the SKOV3 cell line, a cell line derived from an ovarian carcinoma obtained from Dr. Claudine Rancourt (17). These cells were grown as previously described (3, 17, 30).
qRT-PCR
Tissues, RNA purification, quantification, and RT of samples used for real-time qRT-PCR. Total RNA from our RNA bank was originally purified using TriPure (Roche; Indianapolis, IN) and stored at –80°C. The integrity of triplicate total RNA samples of proliferatively active HIEC-6 cells, differentiated fetal epithelial cells (see Tissues), and Caco-2/15 cells at several differentiation stages [subconfluent (SC), 0 days postconfluent (PC), 2–5 days PC, 10–15 days PC, and 25–35 days PC] was verified by visual analysis of 18S and 28S bands on an ethidium bromide-stained 1.5% agarose gel. The same procedure was followed for the seven patient-matched pairs of resection margins (RM) and primary adenocarcinoma samples. These samples were provided by the Cooperative Human Tissue Network, which is funded by the National Cancer Institute, and consisted of one well-differentiated, five moderately differentiated, and one poorly differentiated colon adenocarcinoma. It is noteworthy that in contrast to the cell lines and villus cells representing pure epithelial cell preparations, the colonic tissues used herein were heterogeneous, containing both epithelial and stromal cells. Because the RT reaction is extremely sensitive to variations in salt concentration and contamination by carry-over reagents from previous steps (15), the samples were repurified using RNeasy columns (Qiagen; Mississauga, ON, Canada) to obtain equivalent salt concentrations and purity of the samples. The RiboGreen RNA Quantitation Kit (Molecular Probes; Hornby, ON, Canada) was used according to the manufacturer's instructions to accurately quantify the samples. Fluorescence emission was read by an Mx3000P (Stratagene; La Jolla, CA) functioning as a plate reader. First-strand cDNA synthesis using Omniscript (Qiagen) was performed on 2 µg total RNA using a combination of oligo(dT)12–18 and random hexamers (Amersham Pharmacia; Baie d'Urfé, QC, Canada) as primers. To ensure an RT performance of equal quality, all samples were reverse transcribed simultaneously using a single mastermix.
Real-time PCR.
All reactions were performed in an Mx3000P (Stratagene) starting with 10 min of Taq activation at 95°C, followed by 40 cycles of melting (95°C, 30 s), primer annealing at the temperature appropriate for each primer (55–60°C, 45 s), and extension (72°C, 45 s) ending with a melting curve analysis to validate the specificity of the PCR products. Fluorescence data were acquired after each annealing step. Amplification efficiencies ranged from 93% to 104%. The genes investigated were 5-methyltetrahydrofolate-homocystine methyltransferase (MTR; Accession No. NM_000254), RNA polymerase II (DNA directed) polypeptide A (POLR2A; Accession No. NM_000937), B2M (Accession No. NM_004048), -actin (ACTB; Accession No. NM_001101), GAPD (Accession No. NM_002046), ribosomal protein S14 (RPS14; Accession No. NM_005617), mannosidase-
, class 1B, member 1 (MAN1B1; Accession No. NM_016219), RPLP0 (Accession No. BC019014), hypoxanthine phosphoribosyltransferase 1 (HPRT1; Accession No. NM_000194), sucrase-isomaltase (-glucosidase) (SI; Accession No. NM_001041), and c-Myc (MYC; Accession No. NM_002467). Table 1 summarizes the primers used in this study. Primers were generated using the primer formation software Primer3 (http://frodo.wi.mit.edu/cgi-bin/primer3/primer3_www.cgi) with attention given to avoiding primer-dimer formation by stringent use of the maximum 3' self-complimentarity function of the Primer3 program. The absence of primer-dimers was verified postamplification by melting curve analysis.
|
All statistical analyses were performed using the statistical software Graph Pad Prism 3.02 (Graph Pad Software; San Diego, CA). For linear regression analysis, we applied the following numeric denotation to the different cellular samples representing different stages of differentiation: HIEC: 1, Caco-2/15 SC: 2, Caco-2/15 0 PC: 3, Caco-2/15 2–5 PC: 4, Caco-2/15 10–15PC: 5, Caco-2/15 25–35 PC: 6, and villus epithelial: 7. These numeric values were then plotted against the average cycle threshold (Ct) values obtained for the individual gene at the corresponding differentiation stage. Goodness of fit analysis (r2) was performed on the obtained linear regression lines to evaluate the strength of the regression and to assess the usability of the individual genes as normalizing genes during differentiation. ANOVA was performed as a one-way analysis of variance without post hoc testing, and paired t-tests were performed as two-tailed tests. A P value below <0.05 was considered significant in all analyses. The Normfinder applet (2) was courtesy of Dr. Ørntoft (Aarhus University Hospital; Aarhus, Denmark). Briefly, this applet describes expression values and performs both a separate analysis of the expression values in sample subgroups as well as an estimation of the inter- and intragroup variations providing a "stability factor." The geNorm applet (32) was obtained from the Gene Quantification homepage (www.gene-quantification.org). Briefly, this applet performs an analysis of expression stability leading to a stability index value. Furthermore, the applet provides a normalization factor based on the geometric average of all input normalization genes.
|
RESULTS |
|---|
|
TOPABSTRACTMATERIALS AND METHODSRESULTSDISCUSSIONGRANTSREFERENCES |
|---|
Samples represent different stages of differentiation.
We initially verified that our stringently purified and reverse transcribed RNA samples truly represented differentiating cells by analyzing the expression level of the brush-border enzyme SI, a well-known and accepted marker of differentiation. Bescause our samples were stringently purified, quantified, and reverse transcribed, we compared the obtained Ct values directly. The Ct value represents the number of thermal cycles needed to reach a certain threshold value of fluorescence and is inversely correlated to the amount of input cDNA. Under perfect conditions, a 
Ct of 1 equals a doubling of input material. The HIEC-6 and SC Caco-2/15 samples did not yield detectable amplicons and were thus attributed the nomination of nondetectable (ND). As expected, Fig. 1A demonstrates a dramatic increase in transcript levels as a function of differentiation. The decline shown in the Ct value corresponds to an 
325-fold increase in gene expression from newly confluent Caco-2/15 cells to Matrisperse-extracted villus epithelial cells, demonstrating that the samples indeed represent intestinal epithelial cells at different stages of differentiation.
|
|
Analysis of intra- and intergroup variance. m And r2 values are obtained using the full data set and are thus an "average" of the entire data set. Thus a certain tolerance of outliers that can "disappear" in the full data set is a precondition of these methods. To evaluate if any such outliers were present at any stage(s) of differentiation, we performed one-way ANOVA on the data from the nine candidate genes and found that five of the genes had statistically significantly (P < 0.05) different levels of expression from one differentiation stage to another. The four genes that did not have statistically relevant alterations of expression from differentiation stage to differentiation stage were RPLP0, B2M, RPS14, and MTR (Table 2).
We next evaluated the SD of the Ct values of the nine candidate genes. The four genes with the lowest SDs were RPLP0, B2M, RPS14, and ACTB (0.43, 0.58, 0.61, and 0.61, respectively), whereas POLR2A, MTR, and HPRT1 were the only genes scoring SD values higher than 1.000 (Table 2), representing an input cDNA variation of more than twofold (24). These three genes all showed high variation in the nondifferentiated stages, whereas the differentiated stages displayed considerably less variation (Fig. 1B). To exclude the possibility of experimental error as the cause of the high variation, we repeated the PCR run, resulting in analogous results (data not shown). Finally, we analyzed the data sets using the Excel applets Normfinder and geNorm, which take into account both the intra- and the inter-group variation when calculating the stability factor of a given set of data (2, 32). This approach identified RPLP0 and RPS14 as the two genes with the least variability (Table 2). Overall, we thus found that RPLP0 was the gene that most frequently was identified as having the least variation using all the different statistical approaches.
Analysis of normalizing gene modulation in differentiating Caco-2/15 cells only. To further characterize the data sets for enterocyte differentiation, we investigated the same parameters as above on the nine candidate genes for the Caco-2/15 cell line only. The results are summarized in Table 2. POLR2, HPRT1, RPS14, and MTR displayed a strong tendency toward modulation during the course of Caco-2/15 differentiation (m = –0.82, –0.35, –0.31, and –0.52, respectively), whereas GAPD, MAN1B1, and, to a certain extent, ACTB showed less modulation (m = –0.04, 0.08, and 0.12, respectively). RPLP0 retained a low m value (–0.07), whereas B2M showed less stable expression during Caco-2/15 differentiation (–0.13; Table 2) than with all three cell lines (–0.03; Table 2). Likewise, it was MAN1B1, GAPD, and RPLP0 that scored the lowest values for the r2 test (0.05, 0.01, and 0.04, respectively) with Caco-2/15 cells only. Using ANOVA, we also obtained statistically significant differences in expression between differentiation stages for only POLR2A and HPRT1 (0.01 and 0.03, respectively). The SD data were not significantly modulated by the exclusion of the data from the HIEC and villus cells except for HPRT1 (0.76 vs. 1.47). The three genes with the lowest SDs were RPLP0, MAN1B1, and ACTB (0.43, 0.4, and 0.41, respectively). Finally, analysis of stability values showed RPLP0, RPS14, and MAN1B1 to be the three genes with the best scores (0.31, 0.28, and 0.29, respectively).
Healthy Tissue Versus Adenocarcinomas
Cancer samples were first tested for the expression of the oncogene MYC in the tumor samples compared with the corresponding RMs. As expected (13), we found a significant upregulation of the expression of MYC in the cancer samples compared with the RMs (Fig. 2A). The Ct value of 1.69 between the RMs and adenocarcinoma samples represented a 3.22-fold (21.69, considering 100% PCR efficiency) upregulation (P < 0.01).
|
Ct values from the patient-paired samples to absolute (positive) numbers to analyze the average gene modulation between the RMs and adenocarcinomas for the nine candidate genes. As shown in Fig. 2B, B2M scored the lowest overall absolute Ct value at 0.32. Three other genes stood out with low differences in expression between the tumors and RMs: ACTB, MAN1B1, and MTR (Ct: 0.43, 0.49, and 0.48, respectively). In particular, POLR2A, which displayed hardly any modulation when the crude average was observed (Fig. 2A), showed a high modulation of gene expression using absolute values (Fig. 2B). Expression of genes of interest: averages of normalizing genes versus selected normalizing genes. The expression levels of genes of interest for both intestinal cell differentiation (SI) and colon adenocarcinomas (MYC) were analyzed relative to the geometric average (23, 31) of either all nine analyzed normalizing genes or the best three or four normalizing genes as determined previously (see above) versus the single best gene for each experimental setting. As shown in Table 3, all three approaches led to similar conclusions, although some differences in the levels of induction were observed.
|
|
DISCUSSION |
|---|
|
TOPABSTRACTMATERIALS AND METHODSRESULTSDISCUSSIONGRANTSREFERENCES |
|---|
We employed several different statistical analyses to analyze the levels of gene expression in the intestinal cell samples representing distinct differentiation stages. Distinct statistical analyses were performed to consider complementary parameters and information relative to each gene. Individually, this allowed us to address specific parameters such as trend of expression, internal variation, and relative variation as well as inter- and intragroup variation. Collectively, these analyses led to the identification of optimal housekeeping genes. Indeed, some of the genes, such as POLR2A, MTR, and HPRT1, were rapidly discarded based on their high SD values, which could be ascribed to biological rather than input cDNA variation because it was not seen with the other genes. In fact, considering all the various analyses performed, only three genes stood up to our assessment: RPLP0, RPS14, and B2M. Both RPLP0 and B2M had near-zero slope values; however, B2M was, in our opinion, difficult to accept as the "gold standard" for normalization during differentiation of enterocytes due to the relative high variation observed in the undifferentiated samples, which was further reflected by the poor stability value obtained using the Normfinder and geNorm applets (2, 32). These analytic tools scored RPS14 and RPLP0 as the most stable normalizing genes, although RPS14 displayed a tendency toward a negative slope. Therefore, taking into consideration all the parameters, RPLP0 was selected on the basis that it clearly displayed the best overall stability. Although other experimental conditions may require further validation for the identification of housekeeping genes, it is noteworthy that these analyses were performed across two very different cell lines (HIEC-6 and Caco-2/15) as well as in primary tissue (epithelial villus fractions of the small intestine) strengthening the potency of RPLP0 as the ideal normalizing gene for mRNA quantification in human intestinal epithelial cells. Interestingly, performing the same analyses on Caco-2/15 cells only led to the identification of additional genes, such as GAPD and MAN1B1, that appeared to perform nearly as well as RPLP0. This finding of stable expression in a larger subset of candidate genes in a single cell line is consistent with a more homogeneous cellular background. This further underscores both the remarkably stable expression of RPLP0 across all the cell lines as well as the necessity to carefully validate potential normalizing genes under particular experimental settings.
For the identification of normalizing genes to be used for studies related to adenocarcinomas in the colon, we mainly considered the average modulation obtained using the absolute Ct values rather than the crude average, which does not always reflect individual pair variation. With this approach, we found that the classic normalizing gene B2M performed slightly better than others such as ACTB, MTR, and MAN1B1, whereas other classic normalizing genes such as GAPD and POLR2A appeared to significantly vary between control and tumor tissue. Minute alterations in B2M expression as a consequence of tumor development is in accordance with reported results from the colon (6) and the kidney (14). The strength of these studies, including ours, is the use of matched primary tumors and corresponding histologically normal RMs, enabling an accurate assessment of the difference in gene expression levels between healthy and adenocarcinoma tissue. Our finding that ACTB is little affected by tumor formation is also in agreement with previous studies (6, 14). Therefore, based on the current analyses, it appears that B2M represents the best normalizing gene for studies related to gene expression in adenocarcinomas of the human colon, although ACTB, MTR, and MAN1B1 represent acceptable alternatives. It cannot, however, be excluded that these normalizing genes, as well as others not tested herein, might also be acceptable controls for distinct subtypes of colonic malignancy.
In conclusion, we herein identified RPLP0 as a very well suited, if not ideal, normalizing gene for mRNA quantification in human intestinal epithelial cells at various stages of differentiation, whereas B2M was confirmed to be the best housekeeping gene for studying gene expression in human colon primary carcinomas. Interestingly, a comparison of the expression levels of genes of interest for intestinal cell differentiation (SI) and colon carcinomas (MYC) using these single optimal housekeeping genes versus an average of several normalizing genes (23, 31) led to similar results. The identification of different optimal normalizing genes for studies of differentiation versus adenocarcinoma formation is not surprising considering the fact that the differentiation studies were performed on pure epithelial cell cultures and extracts, whereas the mRNA for cancers and resections margins were extracted from mucosa. In support of this, differences in gene expression profiles between primary cancers and the surrounding stroma have been recently demonstrated (29).
|
GRANTS |
|---|
|
TOPABSTRACTMATERIALS AND METHODSRESULTSDISCUSSIONGRANTSREFERENCES |
|---|
|
ACKNOWLEDGMENTS |
|---|
|
FOOTNOTES |
|---|
The costs of publication of this article were defrayed in part by the payment of page charges. The article must therefore be hereby marked "advertisement" in accordance with 18 U.S.C. Section 1734 solely to indicate this fact.
|
REFERENCES |
|---|
|
TOPABSTRACTMATERIALS AND METHODSRESULTSDISCUSSIONGRANTSREFERENCES |
|---|
This article has been cited by other articles:
|
|
 |
|
  L A Simms, J D Doecke, M D Walsh, N Huang, E V Fowler, and G L Radford-Smith Reduced {alpha}-defensin expression is associated with inflammation and not NOD2 mutation status in ileal Crohn's disease Gut, July 1, 2008; 57(7): 903 - 910. [Abstract] [Full Text] [PDF]
|
|
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |
| Visit Other APS Journals Online |
