- Methodology article
- Open Access
Identification of suitable endogenous control genes for microRNA gene expression analysis in human breast cancer
BMC Molecular Biology volume 9, Article number: 76 (2008)
The discovery of microRNAs (miRNAs) added an extra level of intricacy to the already complex system regulating gene expression. These single-stranded RNA molecules, 18–25 nucleotides in length, negatively regulate gene expression through translational inhibition or mRNA cleavage. The discovery that aberrant expression of specific miRNAs contributes to human disease has fueled much interest in profiling the expression of these molecules. Real-time quantitative PCR (RQ-PCR) is a sensitive and reproducible gene expression quantitation technique which is now being used to profile miRNA expression in cells and tissues. To correct for systematic variables such as amount of starting template, RNA quality and enzymatic efficiencies, RQ-PCR data is commonly normalised to an endogenous control (EC) gene, which ideally, is stably-expressed across the test sample set. A universal endogenous control suitable for every tissue type, treatment and disease stage has not been identified and is unlikely to exist, so, to avoid introducing further error in the quantification of expression data it is necessary that candidate ECs be validated in the samples of interest. While ECs have been validated for quantification of mRNA expression in various experimental settings, to date there is no report of the validation of miRNA ECs for expression profiling in breast tissue. In this study, the expression of five miRNA genes (let-7a, miR-10b, miR-16, miR-21 and miR-26b) and three small nucleolar RNA genes (RNU19, RNU48 and Z30) was examined across malignant, benign and normal breast tissues to determine the most appropriate normalisation strategy. This is the first study to identify reliable ECs for analysis of miRNA by RQ-PCR in human breast tissue.
Approximately 98% of the human transcriptome is non-protein-coding RNA (ncRNA) [1, 2]. The fraction of ncRNA believed to be functional in cells was once limited to the well-characterised transfer and ribosomal RNAs. However, this fraction has recently been expanded to include microRNAs (miRNAs), a class of short, single-stranded RNAs that target, through nucleotide complementarity, specific messenger RNAs, enabling them to negatively modulate gene expression. First characterised in 1993, miRNAs were initially shown to be involved in the control of developmental timing in Caenhorhabditis elegans . Now, 15 years later, 541 human miRNAs have been submitted to the most recent edition of the online miRNA sequence repository, miRBase [4, 5]. A very small proportion of identified miRNAs have verified roles, in processes such as cell proliferation  and apoptosis . It may be some time before the full catalogue of biologically-functional miRNAs is compiled but the involvement of miRNAs in the regulation of cancer-related genes at the post-transcriptional level has already led to miRNAs being hailed as a novel class of tumour suppressor genes and oncogenes and the coining of the term "oncomiR" [reviewed in  and ]. Further elucidating our limited understanding of the mechanisms of metastasis [10, 11]; the spread of cancer cells from the primary neoplasm to distant organs, has become a major focus in miRNA studies [12–14]. Gene signatures from these studies will contribute to our understanding of the multi-step processes of metastasis and may also enable advanced indication of the likelihood of tumour invasion and metastasis based on the characteristics of the primary tumour.
The expression of miRNAs has been studied using traditional, semi-quantitative methods such as northern blotting , bead-based flow-cytometry  and microarray technology . However, by far the method of choice for expression quantitation is real-time quantitative PCR (RQ-PCR) due to its sensitivity, wide dynamic range and low template requirements. The technique has itself been revolutionised in recent years with the development of stem-loop primers that specifically convert the mature, functional miRNA into its DNA complement . Furthermore, a multiplex stem-loop RQ-PCR format is currently being refined to allow multiple miRNAs to be transcribed simultaneously [19, 20].
To produce reliable relative RQ-PCR data, corrections must be made for variation between reactions introduced during the steps from sample preparation to amplification. Incorporating an endogenous control (EC) gene into the experimental design is an effective method of normalising the data but candidate ECs must be tested on a representative number of the sample population if not the entire sample set [21–23]. Use of an unreliable EC may lead to inaccurate, unreliable results and previous studies show that mRNA expression can be made to appear up- or down-regulated by one order of magnitude based solely on the choice of EC .
Using the Medical Subject Heading (MeSH) terms microRNAs, neoplasm and reverse transcriptase polymerase chain reaction, a recent PubMed search returned 42 articles, 5 of which detailed miRNA RQ-PCR expression profiling studies using human neoplastic breast samples. The EC(s) used in these studies were let-7a and miR-16 , U6 small nuclear RNA (snRNA) and tRNA for initiator methionine , 18S rRNA , glyceraldehyde-3-phosphate dehydrogenase (GAPDH)  and in one article, the EC used for RQ-PCR analysis was not given . There is currently no consensus on suitable ECs for quantitative analysis of miRNA by RQ-PCR in human breast tissue. Concern has been expressed regarding the use of ribosomal RNAs in normalisation strategies as they can be expressed at much greater levels than the target RNA making it very difficult to quantify an rRNA and a rare transcript in the same RNA dilution [22, 29, 30]. Moreover, there is evidence of rRNA deregulation in apoptosis . It is also clear that the traditionally-used but seldom validated GAPDH and β-actin (ABTB) genes are not suitable endogenous controls for some studies [32–34]. Our aim was to identify from a panel of RNA species similar to miRNAs, suitable EC candidates for miRNA analysis by RQ-PCR.
The expression of eight small RNAs was determined in 36 fresh-frozen breast tissues; three small nucleolar RNAs (snoRNAs, RNU19, RNU48 and Z30) and five miRNAs (let-7a, miR-10b, miR-16, miR-21 and miR-26b). The five miRNAs were chosen as candidate ECs for this study based on their known expression in human breast tissue and/or their previous use as an EC gene for miRNA RQ-PCR analysis . The three snoRNA genes were selected from a panel of ten commercially available TaqMan control assays. SnoRNAs, found within the nucleolus, range from 60–300 nt in length and chemically modify rRNA , small nuclear RNA (snRNA)  and mRNA  through their recognition of short sequences in the target molecule and recruitment of associated proteins to the site. MiR-30*, previously referred to as miR-30a-3p, targets RNA involved in several cancer-related biological processes  and was chosen as a target gene to investigate the effect of EC gene selection on relative quantitation. Downregulation of miR-30* has previously been shown to increase transcription of mRNAs involved in processes such as angiogenesis (thrombospondin I, cysteine-rich, angiogenic inducer) and cell cycle transition (cyclin-dependent kinase 6) . Samples consisted of malignant (n = 21), benign (n = 5) and normal (n = 5) breast tissues. Malignant tumour tissues were representative of all tumour grades and hormone receptor states. The malignant breast tumour tissues were divided into three groups depending on the patient's disease progression in the five years following removal of the primary tumour; tumours from patients who did not develop metastases (metastases-free, MF, n = 13), those from patients who developed bone metastases (BM, n = 7) and those from patients who developed visceral and bone metastases (VBM, n = 6).
Concentrations of small RNA ranged from 45–431 ng μL-1. The distribution of miRNA yields (RNA in the 10–35 nt range) across the various malignant groups, benign and normal tissues was determined using the Agilent Small RNA Assay (Table 1). Percentage of miRNA in the small RNA fraction ranged from 12%–98%. For the majority of samples, miRNA comprised 26–75% of total small RNA. High miRNA yields (> 75% miRNA) were seen in samples from the BM and MF groups. Low yields (< 25% miRNA) were seen in all groups apart from the normal tissue group. The large RNA samples, extracted separately but at the same time as the small RNA used for this study, all had an RNA integrity number (RIN) ≥ 7.
Relative Expression Quantitation
The threshold cycle (Ct) is the amplification cycle number at which the fluorescence generated within a reaction rises above a defined threshold fluorescence . The eight candidate ECs displayed a wide expression range with Ct values between 18 and 36. MiR-16 and miR-21 showed relatively high expression with median Cts of 21, while let-7a, miR-10b, miR-26b and RNU48 were moderately abundant with median Cts of between 23 and 27. Z30 had lower abundance with a median Ct of 29. RNU19 expression was very low in these samples with Cts ranging from 26.2 to 38.9. It was decided to exclude RNU19 from further analysis due to its low expression. Ct values were converted to relative quantities (Q.Rel) using the formula: E-ΔCt, where E = PCR amplification efficiency and ΔCt = average Ct, test sample-average Ct, calibrator sample. There was no significant difference in variance between genes (P > 0.05, Fig. 1B). The relative quantities did not differ significantly between the MF, BM, VBM and BEN groups for any of the candidate ECs (P > 0.05; Fig. 1A). As NormFinder and geNorm assume candidates are not differentially expressed between groups, this analysis is necessary to validate use of these methodologies . Equivalent expression between the tumour tissues (benign and malignant) and the normal breast tissues (used as the calibrator) was confirmed for each of the seven candidate ECs using the equivalence test and a fold change cut-off of 3. All genes, with the exception of Z30, were also equivalently expressed between the malignant and benign tumour groups using a fold change cut-off of 3.
Stability Candidate EC Expression
Stability of the candidate ECs was examined using geNorm and NormFinder. The ranking of the candidates, as determined by these programmes, is shown in Table 2. The lower the stability value, the higher the gene stability. With a stability value of 0.312, NormFinder selected let-7a as the most stably expressed single gene. The best combination of two genes, let-7a and RNU48, further reduced the NormFinder stability value to 0.221. Two out of the seven genes showed geNorm stability measures (M) below the default limit of 1.5 and the programme identified let-7a and miR-16 as the most stable pair of ECs (M = 0.978).
Effect of EC on Relative Quantity of miR-30*
To assess the effect of EC on relative quantitation of the target gene, miR-30*, this miRNA was normalised using each of the candidate EC genes in turn. Depending on the normaliser, miR-30* expression was either significantly different between tissue groups (P < 0.05) or no differences were detected. When normalised using miR-26b, ranked in the top four mosts Table 3 andidates by both geNorm and NormFinder (Table 2), no differences were detected between tissue groups. Normalisation to all other individual ECs detected significant differences between the BM and VBM tissue groups. Only normalisation with RNU48 detected a significant difference in miR-30* expression between the MF and VBM tissue groups. GeNorm selected let-7a and miR-16 as the most stable EC pair and let-7a was selected as the single most stable EC gene using NormFinder. Thus the effect of using either let-7a as a single gene or using the recommended EC pair, let-7a and miR-16, on miR-30* expression was assessed. Significant differences in miR-30* expression were detected between tissue groups using either the one EC (P = 0.007) or the two EC (P = 0.01) approach, however the BM and MF tissue groups were found to be significantly different using the EC pair, let-7a and miR-16 but this was not detected when let-7a was used as the sole EC gene. The lowest pairwise variation value (V) was generated using the top five genes from the panel, indicating that this would be the most stable EC gene set to use. MiR-30* was normalised using the top two EC genes and using the top five EC genes to assess what effect this would have on miR-30* relative quantification. Significant differences in miR-30* expression were detected between the tissue groups using either the top two ECs (P = 0.01) or the top five ECs (P = 0.002, Fig. 2B), however the post-hoc analyses varied slightly in that the two gene normalisation detected a difference between the BEN and MF tissue groups not detected by the five EC gene approach. Conversely, the 5 gene approach identified a significant difference in miR-30* expression between the MF and VBM tissues, not detected when using the most stable pair. Both normalisation strategies did detect significant differences between the BM vs MF groups, the BM vs VBM groups and the BEN vs VBM groups.
In recent years, miRNAs have emerged as key players in tightly-controlled biological processes such as proliferation , apoptosis [7, 42] and tumour invasion . MiRNAs, first implicated in malignancy in 2002 , are known to be deregulated and/or mutated in numerous cancers including breast cancer  and there is evidence to suggest that miRNA expression profiles may be more accurate in classifying breast cancer subtypes than mRNA expression profiles .
Adaptation of traditional RNA isolation and reverse transcription protocols has facilitated the application of RQ-PCR to the study of miRNA expression. Mature miRNAs were amplified and quantified using PCR for the first time in 2005 [18, 45] and recent developments include a 220-plex RQ-PCR allowing the analysis of multiple miRNAs from single cells . The high sensitivity of RQ-PCR demands appropriate normalisation to correct for non-biological variation and the use of EC genes remains the most commonly used method. The issue of carefully selecting and validating EC genes has already been discussed for a number of experimental systems in the context of RQ-PCR for mRNA [21, 29, 46] however, this issue has not yet been addressed in relation to the relative quantitation of miRNA in breast tissue. This issue is particularly pertinent to the area of miRNA studies utilising RQ-PCR since it is still common practice to synthesise a gene-specific cDNA for each sample, using miRNA-specific primers, thereby introducing additional non-biological variation not accrued during the synthesis of cDNA from mRNA when using random or oligo-dT primers.
This paper describes the first systematic assessment of candidate ECs for the normalisation of miRNA RQ-PCR data in breast cancer. In the rapidly evolving field of miRNA reasearch, consensus has not yet been reached on how best to tackle this issue. Numerous RNA species, including rRNAs, tRNAs, snRNAs and miRNAs have previously been used as ECs in miRNA RQ-PCR studies of the breast. Concern has been expressed regarding the use of rRNAs in normalisation strategies as they can be expressed at much greater levels than the target RNA making it very difficult to quantify a rRNA and a rare transcript in the same RNA dilution [22, 29, 30]. Moreover, a role for rRNA in apoptosis  and cancer  has been reported. A proportion of snRNAs and snoRNAs may exhibit tissue-specific and developmental regulation  emphasising the need for validation of commercially-available control assays. U6 snRNA (RNU6B), commonly used to normalise miRNA RQ-PCR data [49, 50] was found to be less stably expressed than let-7a and miR-16, the EC pair proposed by this study .
This is the first report detailing the percentage of miRNA retrieved in small RNA-enriched fractions of primary breast tissue. RNAs detected using the Agilent small RNA assay include miRNAs, smaller ribosomal RNAs such as the 5.85S (154 nt) and 5S (117 nt) subunits, transfer RNAs (73–93 nt) and snoRNAs (60–300 nt). Perhaps unsurprisingly, we found the proportion of miRNA in the small RNA sample ranged from 12–98%. The varying miRNA yields were well distributed amongst the tissue groups. The variation in ratio was not dependent on the type of tissue, on the RNA extraction or on the total yield of the RNA. In this laboratory we also found a much lower proportion of miRNA, ranging from 1–5%, in small RNA extracted from commonly used breast cancer cell lines such as MCF-7, SK-BR-3, T47D and ZR-75-1 (data not shown). A suitable EC gene will have to reflect such changes in the global miRNA population. The variation in miRNA yields and ratios may reflect genomic alterations, common in cancer (reviewed in . It has been shown previously that miRNA frequently map to such regions of instability . This finding raises concerns over how much of the small RNA used for cDNA reactions and other applications is actually the RNA species of interest, and this is especially relevant in studies employing non-miRNA ECs.
Normfinder and geNorm were used to identify suitable ECs for the relative quantitation of miRNA in fresh-frozen primary breast tissue. There was no effect of tissue group on scaled EC expression (P > 0.05, Fig. 1). As previously stated  this is an important validation prior to use of geNorm and NormFinder as these models assume candidates are not differentially expressed between experimental groups. The absence of a significant difference in EC expression between groups does not necessarily equate to equivalent expression. Equivalence of expression was assessed using the equivalence test . Equivalent expression between the independent tumour and normal breast tissues was confirmed for all ECs using a fold change cut-off of 3. Equivalent expression between the malignant and benign tumour groups was also assessed and was confirmed for all ECs with the exception of Z30 using the same cut-off. Using the benign, MF, BM and VBM subgroups, NormFinder calculated the intra- and intergroup variations and identified let-7a as the single most stable EC with a stability value of 0.312. However, the use of more than one EC is believed to increase the accuracy of quantitation compared to the use of a single EC [30, 54, 55] and use of let-7a alone would therefore not be recommended. The EC gene pair, let-7a and RNU48 had an improved NormFinder stability value of 0.221.
GeNorm generates a gene-stability measure (M) which may be defined as the average pairwise variation (V) for one candidate EC gene compared to all other candidate EC genes. Stepwise exclusion of the gene with the highest M value results in recalculation of M values for the remaining genes and ultimately, the identification of the most stable pair . The wide range in M values depicts the high variability detected in candidate EC gene stability (Table 2). The differences in miR-30* expression detected between the tissue groups varied greatly depending on which single EC was used for normalisation. For example, in the BEN and BM breast tissues, the expression of miR-30* could be made to appear up- or down-regulated relative to normal breast tissue depending on the EC gene used (Fig. 2A). These results draw particular attention to the potential effect of EC choice on the outcome of a study and demonstrates the need for validation of candidate ECs to produce reliable expression data. In geNorm, a normalisation factor (NF) is generated for each sample using the geometric average of the expression of the most stable EC genes. The pairwise variation value, V is the variation between two sequential NFs (Vn/n+1, where n = the number of ECs used). The recommended pairwise variation of 0.15 is a guideline value and is not intended as an absolute cut-off. This guideline value may not always be achievable  but should be considered, particularly if small expression differences are to be measured. The lowest pairwise variation value was achieved when the top 5 candidate genes were used as ECs (0.220, Fig. 3B). Let-7a and miR-16 were identified as the most stable pair of EC genes using geNorm. Significant differences in miR-30* expression were detected between the tissue groups using either the top two ECs (P = 0.01) or the top five ECs (P = 0.002, Fig. 2B). The post-hoc analyses revealed both normalisation approaches detected significant differences between the BM vs MF groups, the BM vs VBM groups and the BEN vs VBM groups but each approach detected an additional intergroup difference not detected by the other approach. The number of genes to use in a normalisation strategy is in most cases, a trade off between required resolution and practicality and for most purposes the EC gene combination let-7a and miR-16 should suffice. Both genes had the lowest stability values, as determined by geNorm and NormFinder (Table 2).
A tumour suppressor role for let-7a in lung tissues seems likely due to its widespread downregulation in tumour versus normal lung tissues as well as the identification of an oncogenic target, RAS, in this tissue . However, it is unclear whether let-7a is implicated in breast cancer since the results of recent studies have been equivocal . Whilst deletion of the miR-16 gene has been implicated in the development of chronic lymphocytic leukemia , a specific role for this miRNA in breast cancer has not been identified. From a panel of 345 miRNAs, miR-16 was selected in the top 15 most stably-expressed miRNAs across 40 normal human tissue types . A microarray study  which looked at the expression of 287 miRNAs in various normal and tumour tissues, not including breast tissue, selected a panel of suitable EC genes based on a number of criteria including high and consistent expression of the miRNA across the tissues. Depending on the tissue sample set, both let-7a and miR-16 were ranked in the top 10–15 most stably-expressed miRNAs, supporting the findings of the present study.
The tissues used in this study are clinically and pathologically diverse (see Table 3) making this study of interest to a broad spectrum of the breast cancer research community. Recent findings would suggest that, unlike mRNAs, the miRNA fraction present in FFPE tissues is relatively unaffected by the fixation process and that miRNAs extracted from these tissues may be accurately profiled using RQ-PCR . Thus, the ECs identified in this study may also prove useful for miRNA RQ-PCR analysis of FFPE breast tissues.
MiRNA expression studies utilising RQ-PCR should begin with the careful selection of appropriate ECs for normalisation to ensure accurate quantitation of this very exciting class of molecules. This study indicates an appropriate strategy to validate ECs for any miRNA RQ-PCR study and has identified a reliable two-gene normaliser for use in breast cancer studies. We recommend the combined use of Let-7a and miR-16 in this context.
Primary breast tumour tissues (n = 31) were obtained from patients during primary curative resection, at Galway University Hospital, Galway, Ireland. Samples were categorised into benign (n = 5) or malignant groups (n = 26) according to standard histopathological parameters. The malignant breast tumour tissues were divided into three groups depending on the patient's disease progression in the five years following removal of the primary tumour; tumours from patients who did not develop metastases (metastases-free, MF, n = 13), those from patients who developed bone metastases (BM, n = 7) and those from patients who developed visceral and bone metastases (VBM, n = 6).
Clinical data relating to the malignant tumour tissues used in this study are shown in Table 3. RNA from normal tissues (n = 5), recovered from patients undergoing reduction mastopexy surgery were used as calibrator samples for RQ-PCR analysis. Tissues were immediately snap-frozen in liquid nitrogen and stored at -80°C until RNA extraction. Prior written and informed consent was obtained from each patient and the study was approved by the ethics review board of Galway University Hospital. Clinical data were obtained from the Breast Cancer Database at the Department of Surgery, Galway University Hospital.
Approximately 100 mg of tissue was homogenised in 1–2 mL QIAzol (Qiagen, Crawley, UK) using a bench-top homogeniser (Polytron PT1600E, Kinematica AG, Littau-Luzem, Switzerland. Large (> 200 nt) and small RNA (< 200 nt) fractions were isolated separately using the RNeasy® Plus Mini Kit and RNeasy MinElute® Cleanup Kit (Qiagen, West Sussex, UK) according to the Supplementary Protocol: Purification of miRNA from animal cells. A portion of the purified large and small RNA was aliquotted for quantitative and qualitative analysis using NanoDrop1000® spectrophotometry and the Agilent 2100 Bioanalyzer respectively. The remaining RNA was stored at -80°C until further use.
MicroRNA concentration and purity were assessed using the NanoDrop1000® spectrophotometer (NanoDrop Technologies Inc, Wilmington, DE, USA). The small-RNA enriched fraction was analysed using the Small RNA Assay with the Agilent 2100 Bioanalyzer (Fig. 4) (Agilent Technologies, Palo Alto, CA, USA). For this assay, samples were diluted to 1 ng/μL, within the quantitative and qualitative range of the assay. Integrity of the large RNA fraction (> 200 nt) was assessed using the RNA 6000 Nano LabChip Series II Assay (Agilent Technologies). An RNA integrity number (RIN) is generated for each sample based on the ratio of ribosomal bands and also the presence or absence of degradation products on the electrophoretic image. A threshold value of RIN ≥ 7 was applied.
Candidate Endogenous Control Genes
The small nuclear and small nucleolar genes were chosen from the ten Human TaqMan MicroRNA Assay Controls available from Applied Biosystems (Foster city, CA, USA) at the time of study. The five miRNA genes were selected based on their known expression in breast cancer tissues and/or based on their previous use as an EC gene in a breast cancer study . Known functions of the candidates are listed in Table 4.
cDNA Synthesis and RQ-PCR
Each reaction was primed using a gene-specific stem-loop primer. The primer and probe sequences for let-7a, miR-10b, miR-16, miR-26b and miR-30* were as previously published . Where sequences were available, primers were obtained from MWG Biotech (Ebersberg, Germany). Otherwise, assays containing the RT stem-loop primer and the PCR primers and probes were used (Applied Biosystems, Foster City, CA, USA). Small RNA (5 ng) was transcribed using MultiScribe Reverse Transcriptase (Applied Biosystems). The reaction was performed using a GeneAmp PCR system 9700 thermal cycler (Applied Biosystems). An RT-negative control was included in each batch of reactions. The PCR reactions were carried out in a final volume of 20 μL using an ABI Prism 7000 Sequence Detection System (Applied Biosystems). Briefly, reactions consisted of 1.33 μL cDNA, 1× TaqMan Universal PCR Master mix, No AmpErase UNG, 0.2 μM TaqMan® Probe (Applied Biosystems), 1.5 μM forward primer and 0.7 μM reverse primer. The PCR reactions were initiated with a 10 min incubation at 95°C followed by 40 cycles of 95°C for 15 sec and 60°C for 60 sec. cDNA, synthesised from pooled normal breast tissue, was included on each 96-well plate as an interassay control and calibrator. All reactions were performed in triplicate. The threshold standard deviation for intra- and inter-assay replicates was 0.3. PCR amplification efficiencies were calculated for each candidate EC RQ-PCR assay using the formula E = (10-1/slope-1) × 100, using the slope of the plot, Ct versus log input of cDNA. PCR amplification efficiencies for each EC candidate are shown in Table 4.
Relative quantities (Q.rel) for each candidate EC gene were calculated from cycle threshold (Ct) values scaled to a calibrator sample (pool of 5 normal tissues) and corrected for efficiency of amplification (E) according to the formula: Q.rel. = E-ΔCt, with ΔCt = average Ct test sample-average Ct calibrator sample. To calculate the expression of the target gene miR-30* relative to each of the EC candidates, the ΔΔCt method was used, with ΔΔCt = (Ct target gene, test sample-Ct endogenous control gene, test sample)-(Ct target gene, calibrator sample-Ct endogenous control gene, calibrator sample). Again, quantities were corrected for efficiency of amplification and the errors were calculated as previously described .
Stability of EC expression was analysed using two freely-available programmes; geNorm and NormFinder. GeNorm (vs. 3.4) is a Visual Basic Application for Excel that calculates a gene-stability measure (M) for all candidate EC genes in a given set of samples and determines the most reliable pair of ECs, showing greatest stability of expression ratio across samples. It is based on a pair-wise comparison model. NormFinder , an excel-add-in, uses an ANOVA-based model to estimate intra- and inter-group variation. It combines these estimates to produce a stability value for each candidate. NormFinder indicates the single most stable EC and EC pair where the stability of the latter is greater than that of the single EC. Prior to geNorm and NormFinder analysis, Ct values were converted into Q.rel values (E = -ΔCt), as detailed above. For Normfinder analysis samples were grouped into metastases free (MF, n = 13), bone metastases (BM, n = 7), visceral and bone metastases (VBM, n = 6) and benign (BEN, n = 5) as described above. Statistical analyses were performed using Minitab (vs 15; Minitab Ltd., Coventry, UK). Distribution of data was determined using the Anderson-Darling normality test and parametric tests were used where appropriate. Levene's statistic was used to assess if there was a significant difference in variance between genes. The equivalence test was used to assess whether genes were equivalently expressed between tumour (benign and malignant) and normal breast tissues and between malignant and benign breast tissues . ANOVA, Fisher's least significant difference tests and Kruskal Wallis tests were applied to determine the effect of EC on target gene expression. P values < 0.05 were considered statistically significant.
Mattick JS: Non-coding RNAs: the architects of eukaryotic complexity. EMBO Rep 2001, 2(11):986-991. 10.1093/embo-reports/kve230.
Mattick JS: Challenging the dogma: the hidden layer of non-protein-coding RNAs in complex organisms. Bioessays 2003, 930-935. 10.1002/bies.10332
Lee RC, Feinbaum RL, Ambros V: The C. elegans heterochronic gene lin-4 encodes small RNAs with antisense complementarity to lin-14. Cell 1993, 75(5):843-854. 10.1016/0092-8674(93)90529-Y.
Griffiths-Jones S, Grocock RJ, van Dongen S, Bateman A, Enright AJ: miRBase: microRNA sequences, targets and gene nomenclature. Nucleic Acids Res 2006, 34(Database issue):D140-4. 10.1093/nar/gkj112.
Griffiths-Jones S: The microRNA Registry. Nucleic Acids Res 2004, 32(Database issue):D109-11. 10.1093/nar/gkh023.
O'Donnell KA, Wentzel EA, Zeller KI, Dang CV, Mendell JT: c-Myc-regulated microRNAs modulate E2F1 expression. Nature 2005, 435(7043):839-843. 10.1038/nature03677.
Cimmino A, Calin GA, Fabbri M, Iorio MV, Ferracin M, Shimizu M, Wojcik SE, Aqeilan RI, Zupo S, Dono M, Rassenti L, Alder H, Volinia S, Liu CG, Kipps TJ, Negrini M, Croce CM: miR-15 and miR-16 induce apoptosis by targeting BCL2. Proc Natl Acad Sci U S A 2005, 102(39):13944-13949. 10.1073/pnas.0506654102.
Kent OA, Mendell JT: A small piece in the cancer puzzle: microRNAs as tumor suppressors and oncogenes. Oncogene 2006, 25(46):6188-6196. 10.1038/sj.onc.1209913.
Esquela-Kerscher A, Slack FJ: Oncomirs - microRNAs with a role in cancer. Nat Rev Cancer 2006, 6(4):259-269. 10.1038/nrc1840.
Fidler IJ: The pathogenesis of cancer metastasis: the 'seed and soil' hypothesis revisited. Nat Rev Cancer 2003, 3(6):453-458. 10.1038/nrc1098.
Gupta GP, Massague J: Cancer metastasis: building a framework. Cell 2006, 127(4):679-695. 10.1016/j.cell.2006.11.001.
Tavazoie SF, Alarcon C, Oskarsson T, Padua D, Wang Q, Bos PD, Gerald WL, Massague J: Endogenous human microRNAs that suppress breast cancer metastasis. Nature 2008, 451(7175):147-152. 10.1038/nature06487.
Ma L, Teruya-Feldstein J, Weinberg RA: Tumour invasion and metastasis initiated by microRNA-10b in breast cancer. Nature 2007, 449(7163):682-688. 10.1038/nature06174.
Budhu A, Jia HL, Forgues M, Liu CG, Goldstein D, Lam A, Zanetti KA, Ye QH, Qin LX, Croce CM, Tang ZY, Wang XW: Identification of metastasis-related microRNAs in hepatocellular carcinoma. Hepatology 2008, 47(3):897-907. 10.1002/hep.22160.
Sempere LF, Freemantle S, Pitha-Rowe I, Moss E, Dmitrovsky E, Ambros V: Expression profiling of mammalian microRNAs uncovers a subset of brain-expressed microRNAs with possible roles in murine and human neuronal differentiation. Genome Biol 2004, 5(3):R13. 10.1186/gb-2004-5-3-r13.
Lu J, Getz G, Miska EA, Alvarez-Saavedra E, Lamb J, Peck D, Sweet-Cordero A, Ebert BL, Mak RH, Ferrando AA, Downing JR, Jacks T, Horvitz HR, Golub TR: MicroRNA expression profiles classify human cancers. Nature 2005, 435(7043):834-838. 10.1038/nature03702.
Miska EA, Alvarez-Saavedra E, Townsend M, Yoshii A, Sestan N, Rakic P, Constantine-Paton M, Horvitz HR: Microarray analysis of microRNA expression in the developing mammalian brain. Genome Biol 2004, 5(9):R68. 10.1186/gb-2004-5-9-r68.
Chen C, Ridzon DA, Broomer AJ, Zhou Z, Lee DH, Nguyen JT, Barbisin M, Xu NL, Mahuvakar VR, Andersen MR, Lao KQ, Livak KJ, Guegler KJ: Real-time quantification of microRNAs by stem-loop RT-PCR. Nucleic Acids Res 2005, 33(20):e179. 10.1093/nar/gni178.
Tang F, Hajkova P, Barton SC, Lao K, Surani MA: MicroRNA expression profiling of single whole embryonic stem cells. Nucleic Acids Res 2006, 34(2):e9. 10.1093/nar/gnj009.
Lao K, Xu NL, Yeung V, Chen C, Livak KJ, Straus NA: Multiplexing RT-PCR for the detection of multiple miRNA species in small samples. Biochem Biophys Res Commun 2006, 343(1):85-89. 10.1016/j.bbrc.2006.02.106.
McNeill RE, Miller N, Kerin MJ: Evaluation and validation of candidate endogenous control genes for real-time quantitative PCR studies of breast cancer. BMC Mol Biol 2007, 8: 107. 10.1186/1471-2199-8-107.
Wong ML, Medrano JF: Real-time PCR for mRNA quantitation. Biotechniques 2005, 39(1):75-85.
Schmittgen TD, Zakrajsek BA: Effect of experimental treatment on housekeeping gene expression: validation by real-time, quantitative RT-PCR. J Biochem Biophys Methods 2000, 46(1-2):69-81. 10.1016/S0165-022X(00)00129-9.
Mattie MD, Benz CC, Bowers J, Sensinger K, Wong L, Scott GK, Fedele V, Ginzinger D, Getts R, Haqq C: Optimized high-throughput microRNA expression profiling provides novel biomarker assessment of clinical prostate and breast cancer biopsies. Mol Cancer 2006, 5: 24. 10.1186/1476-4598-5-24.
Liu CG, Calin GA, Meloon B, Gamliel N, Sevignani C, Ferracin M, Dumitru CD, Shimizu M, Zupo S, Dono M, Alder H, Bullrich F, Negrini M, Croce CM: An oligonucleotide microchip for genome-wide microRNA profiling in human and mouse tissues. Proc Natl Acad Sci U S A 2004, 101(26):9740-9744. 10.1073/pnas.0403293101.
Kulshreshtha R, Ferracin M, Wojcik SE, Garzon R, Alder H, Agosto-Perez FJ, Davuluri R, Liu CG, Croce CM, Negrini M, Calin GA, Ivan M: A microRNA signature of hypoxia. Mol Cell Biol 2007, 27(5):1859-1867. 10.1128/MCB.01395-06.
Chang XZ, Li DQ, Hou YF, Wu J, Lu JS, Di GH, Jin W, Ou ZL, Shen ZZ, Shao ZM: Identification of the functional role of peroxiredoxin 6 in the progression of breast cancer. Breast Cancer Res 2007, 9(6):R76. 10.1186/bcr1789.
Liang Z, Wu H, Reddy S, Zhu A, Wang S, Blevins D, Yoon Y, Zhang Y, Shim H: Blockade of invasion and metastasis of breast cancer cells via targeting CXCR4 with an artificial microRNA. Biochem Biophys Res Commun 2007, 363(3):542-546. 10.1016/j.bbrc.2007.09.007.
Bustin SA: Absolute quantification of mRNA using real-time reverse transcription polymerase chain reaction assays. J Mol Endocrinol 2000, 25(2):169-193. 10.1677/jme.0.0250169.
Vandesompele J, De Preter K, Pattyn F, Poppe B, Van Roy N, De Paepe A, Speleman F: Accurate normalization of real-time quantitative RT-PCR data by geometric averaging of multiple internal control genes. Genome Biol 2002, 3(7):RESEARCH0034. 10.1186/gb-2002-3-7-research0034.
Nadano D, Sato TA: Caspase-3-dependent and -independent degradation of 28 S ribosomal RNA may be involved in the inhibition of protein synthesis during apoptosis initiated by death receptor engagement. J Biol Chem 2000, 275(18):13967-13973. 10.1074/jbc.275.18.13967.
Rondinelli RH, Epner DE, Tricoli JV: Increased glyceraldehyde-3-phosphate dehydrogenase gene expression in late pathological stage human prostate cancer. Prostate Cancer Prostatic Dis 1997, 1(2):66-72. 10.1038/sj.pcan.4500208.
Tokunaga K, Nakamura Y, Sakata K, Fujimori K, Ohkubo M, Sawada K, Sakiyama S: Enhanced expression of a glyceraldehyde-3-phosphate dehydrogenase gene in human lung cancers. Cancer Res 1987, 47(21):5616-5619.
Weisman AS, Tixier-Vidal A, Gourdji D: Thyrotropin-releasing hormone increases the levels of c-fos and beta-actin mRNA in GH3/B6 pituitary tumor cells. In Vitro Cell Dev Biol 1987, 23(8):585-590. 10.1007/BF02620977.
Cavaille J, Bachellerie JP: SnoRNA-guided ribose methylation of rRNA: structural features of the guide RNA duplex influencing the extent of the reaction. Nucleic Acids Res 1998, 26(7):1576-1587. 10.1093/nar/26.7.1576.
Ganot P, Jady BE, Bortolin ML, Darzacq X, Kiss T: Nucleolar factors direct the 2'-O-ribose methylation and pseudouridylation of U6 spliceosomal RNA. Mol Cell Biol 1999, 19(10):6906-6917.
Cavaille J, Buiting K, Kiefmann M, Lalande M, Brannan CI, Horsthemke B, Bachellerie JP, Brosius J, Huttenhofer A: Identification of brain-specific and imprinted small nucleolar RNA genes exhibiting an unusual genomic organization. Proc Natl Acad Sci U S A 2000, 97(26):14311-14316. 10.1073/pnas.250426397.
Nakamoto M, Jin P, O'Donnell WT, Warren ST: Physiological identification of human transcripts translationally regulated by a specific microRNA. Hum Mol Genet 2005, 14(24):3813-3821. 10.1093/hmg/ddi397.
Nolan T, Hands RE, Bustin SA: Quantification of mRNA using real-time RT-PCR. Nat Protoc 2006, 1(3):1559-1582. 10.1038/nprot.2006.236.
Ohl F, Jung M, Radonic A, Sachs M, Loening SA, Jung K: Identification and validation of suitable endogenous reference genes for gene expression studies of human bladder cancer. J Urol 2006, 175(5):1915-1920. 10.1016/S0022-5347(05)00919-5.
Haller F, Kulle B, Schwager S, Gunawan B, von Heydebreck A, Sultmann H, Fuzesi L: Equivalence test in quantitative reverse transcription polymerase chain reaction: confirmation of reference genes suitable for normalization. Anal Biochem 2004, 335(1):1-9. 10.1016/j.ab.2004.08.024.
Chan JA, Krichevsky AM, Kosik KS: MicroRNA-21 is an antiapoptotic factor in human glioblastoma cells. Cancer Res 2005, 65(14):6029-6033. 10.1158/0008-5472.CAN-05-0137.
Calin GA, Dumitru CD, Shimizu M, Bichi R, Zupo S, Noch E, Aldler H, Rattan S, Keating M, Rai K, Rassenti L, Kipps T, Negrini M, Bullrich F, Croce CM: Frequent deletions and down-regulation of micro- RNA genes miR15 and miR16 at 13q14 in chronic lymphocytic leukemia. Proc Natl Acad Sci U S A 2002, 99(24):15524-15529. 10.1073/pnas.242606799.
Huang Q, Gumireddy K, Schrier M, le Sage C, Nagel R, Nair S, Egan DA, Li A, Huang G, Klein-Szanto AJ, Gimotty PA, Katsaros D, Coukos G, Zhang L, Pure E, Agami R: The microRNAs miR-373 and miR-520c promote tumour invasion and metastasis. Nat Cell Biol 2008, 10(2):202-210. 10.1038/ncb1681.
Raymond CK, Roberts BS, Garrett-Engele P, Lim LP, Johnson JM: Simple, quantitative primer-extension PCR assay for direct monitoring of microRNAs and short-interfering RNAs. Rna 2005, 11(11):1737-1744. 10.1261/rna.2148705.
Tricarico C, Pinzani P, Bianchi S, Paglierani M, Distante V, Pazzagli M, Bustin SA, Orlando C: Quantitative real-time reverse transcription polymerase chain reaction: normalization to rRNA or single housekeeping genes is inappropriate for human tissue biopsies. Anal Biochem 2002, 309(2):293-300. 10.1016/S0003-2697(02)00311-1.
Chan MW, Wei SH, Wen P, Wang Z, Matei DE, Liu JC, Liyanarachchi S, Brown R, Nephew KP, Yan PS, Huang TH: Hypermethylation of 18S and 28S ribosomal DNAs predicts progression-free survival in patients with ovarian cancer. Clin Cancer Res 2005, 11(20):7376-7383. 10.1158/1078-0432.CCR-05-1100.
Runte M, Huttenhofer A, Gross S, Kiefmann M, Horsthemke B, Buiting K: The IC-SNURF-SNRPN transcript serves as a host for multiple small nucleolar RNA species and as an antisense RNA for UBE3A. Hum Mol Genet 2001, 10(23):2687-2700. 10.1093/hmg/10.23.2687.
Takamizawa J, Konishi H, Yanagisawa K, Tomida S, Osada H, Endoh H, Harano T, Yatabe Y, Nagino M, Nimura Y, Mitsudomi T, Takahashi T: Reduced expression of the let-7 microRNAs in human lung cancers in association with shortened postoperative survival. Cancer Res 2004, 64(11):3753-3756. 10.1158/0008-5472.CAN-04-0637.
Pineles BL, Romero R, Montenegro D, Tarca AL, Han YM, Kim YM, Draghici S, Espinoza J, Kusanovic JP, Mittal P, Hassan SS, Kim CJ: Distinct subsets of microRNAs are expressed differentially in the human placentas of patients with preeclampsia. Am J Obstet Gynecol 2007, 196(3):261 e1-6. 10.1016/j.ajog.2007.01.008.
Peltier HJ, Latham GJ: Normalization of microRNA expression levels in quantitative RT-PCR assays: Identification of suitable reference RNA targets in normal and cancerous human solid tissues. Rna 2008.
Popescu NC: Genetic alterations in cancer as a result of breakage at fragile sites. Cancer Lett 2003, 192(1):1-17. 10.1016/S0304-3835(02)00596-7.
Calin GA, Sevignani C, Dumitru CD, Hyslop T, Noch E, Yendamuri S, Shimizu M, Rattan S, Bullrich F, Negrini M, Croce CM: Human microRNA genes are frequently located at fragile sites and genomic regions involved in cancers. Proc Natl Acad Sci U S A 2004, 101(9):2999-3004. 10.1073/pnas.0307323101.
Andersen CL, Jensen JL, Orntoft TF: Normalization of real-time quantitative reverse transcription-PCR data: a model-based variance estimation approach to identify genes suited for normalization, applied to bladder and colon cancer data sets. Cancer Res 2004, 64(15):5245-5250. 10.1158/0008-5472.CAN-04-0496.
Bustin SA, Nolan T: Pitfalls of quantitative real-time reverse-transcription polymerase chain reaction. J Biomol Tech 2004, 15(3):155-166.
Goossens K, Van Poucke M, Van Soom A, Vandesompele J, Van Zeveren A, Peelman LJ: Selection of reference genes for quantitative real-time PCR in bovine preimplantation embryos. BMC Dev Biol 2005, 5: 27. 10.1186/1471-213X-5-27.
Johnson SM, Grosshans H, Shingara J, Byrom M, Jarvis R, Cheng A, Labourier E, Reinert KL, Brown D, Slack FJ: RAS is regulated by the let-7 microRNA family. Cell 2005, 120(5):635-647. 10.1016/j.cell.2005.01.014.
Liang Y, Ridzon D, Wong L, Chen C: Characterization of microRNA expression profiles in normal human tissues. BMC Genomics 2007, 8: 166. 10.1186/1471-2164-8-166.
Bortolin ML, Kiss T: Human U19 intron-encoded snoRNA is processed from a long primary transcript that possesses little potential for protein coding. Rna 1998, 4(4):445-454.
Ganot P, Bortolin ML, Kiss T: Site-specific pseudouridine formation in preribosomal RNA is guided by small nucleolar RNAs. Cell 1997, 89(5):799-809. 10.1016/S0092-8674(00)80263-9.
Kiss T, Bortolin ML, Filipowicz W: Characterization of the intron-encoded U19 RNA, a new mammalian small nucleolar RNA that is not associated with fibrillarin. Mol Cell Biol 1996, 16(4):1391-1400.
Galardi S, Fatica A, Bachi A, Scaloni A, Presutti C, Bozzoni I: Purified box C/D snoRNPs are able to reproduce site-specific 2'-O-methylation of target RNA in vitro. Mol Cell Biol 2002, 22(19):6663-6668. 10.1128/MCB.22.19.6663-6668.2002.
Zhou H, Chen YQ, Du YP, Qu LH: The Schizosaccharomyces pombe mgU6-47 gene is required for 2'-O-methylation of U6 snRNA at A41. Nucleic Acids Res 2002, 30(4):894-902. 10.1093/nar/30.4.894.
The authors would like to acknowledge the National Breast Cancer Research Institute (NBCRI) for their continued financial support. AJL is supported by a Clinician Scientist fellowship award from Molecular Medicine Ireland. We gratefully acknowledge Ms. Emer Hennessy for continued technical assistance and for curation of the Department of Surgery, BioBank, NUIG. We also wish to thank Ms. Catherine Curran for collation of clinical and histological data.
PAD performed the experiments, was responsible for data analyses and drafted the manuscript. REM contributed throughout the experiment, critically reviewed the manuscript and participated in data analysis. AL contributed to RQ-PCR analysis and preliminary data analysis. MJK participated clinically in sample provision and in critical examination of the manuscript. NM conceived, designed and supervised experimental work and manuscript editing. All authors read and approved the final manuscript.
About this article
Cite this article
Davoren, P.A., McNeill, R.E., Lowery, A.J. et al. Identification of suitable endogenous control genes for microRNA gene expression analysis in human breast cancer. BMC Molecular Biol 9, 76 (2008). https://doi.org/10.1186/1471-2199-9-76
- Endogenous Control
- Normal Breast Tissue
- Stable Pair
- Endogenous Control Gene
- Tissue Group