This article has Open Peer Review reports available.
Improved Methodology for Assessment of mRNA Levels in Blood of Patients with FMR1 Related Disorders
© Godler et al; licensee BioMed Central Ltd. 2009
Received: 16 December 2008
Accepted: 09 June 2009
Published: 09 June 2009
Elevated levels of FMR1 mRNA in blood have been implicated in RNA toxicity associated with a number of clinical conditions. Due to the extensive inter-sample variation in the time lapse between the blood collection and RNA extraction in clinical practice, the resulting variation in mRNA quality significantly confounds mRNA analysis by real-time PCR.
Here, we developed an improved method to normalize for mRNA degradation in a sample set with large variation in rRNA quality, without sample omission. Initially, RNA samples were artificially degraded, and analyzed using capillary electrophoresis and real-time PCR standard curve method, with the aim of defining the best predictors of total RNA and mRNA degradation.
We found that: (i) the 28S:18S ratio and RNA quality indicator (RQI) were good predictors of severe total RNA degradation, however, the greatest changes in the quantity of different mRNAs (FMR1, DNMT1, GUS, B2M and GAPDH) occurred during the early to moderate stages of degradation; (ii) chromatographic features for the 18S, 28S and the inter-peak region were the most reliable predictors of total RNA degradation, however their use for target gene normalization was inferior to internal control genes, of which GUS was the most appropriate. Using GUS for normalization, we examined in the whole blood the relationship between the FMR1 mRNA and CGG expansion in a non-coding portion of this gene, in a sample set (n = 30) with the large variation in rRNA quality. By combining FMR1 3' and 5' mRNA analyses the confounding impact of mRNA degradation on the correlation between FMR1 expression and CGG size was minimized, and the biological significance increased from p = 0.046 for the 5' FMR1 assay, to p = 0.018 for the combined FMR1 3' and 5' mRNA analysis.
Our observations demonstrate that, through the use of an appropriate internal control and the direct analysis of multiple sites of target mRNA, samples that do not conform to the conventional rRNA criteria can still be utilized to obtain biologically/clinically relevant data. Although, this strategy clearly has application for improved assessment of FMR1 mRNA toxicity in blood, it may also have more general implications for gene expression studies in fresh and archival tissues.
The mutations in an X linked FMR1 gene  have some unique qualities. While large CGG expansions (>200) 3' of its promoter give rise to fragile X syndrome – a neurodevelopmental disorder caused by silencing of the gene, smaller expansions (permutations – PM, 55–200 repeats) lead to late onset pathologies through entirely different mechanisms, involving a toxic gain of function of this gene's transcripts [2–4]. Fragile X Tremour Ataxia Syndrome (FXTAS) is the most prevalent disorder associated with PM alleles occurring in 50% of all older carrier males . It is progressive neurodegenerative disorder manifesting with ataxia and/or intension or other tremors, cognitive decline, psychiatric involvement and characteristic MRI and histopathological features [5–8]. In females, these alleles are associated with a significant increase in premature ovarian failure . A toxicity of an elevated FMR1 mRNA which is believed to be pathogenic in these conditions, has been supported in animal models where it leads to increased cell death in Drosophila  and mice . Although, in humans the level of mRNA is usually assessed in whole blood, there is a significant correlation between FXTAS associated with these levels and the typical neurological , and psychiatric  manifestations. This implies that the assessments of FMR1 mRNA levels in carriers of small CGG expansions, apart from being essential in the mechanisms involved in RNA toxicity, may also have a diagnostic and prognostic significance.
FMR1 expression is primarily regulated by the state of methylation of its promoter. DNMT1 is an enzyme pivotal in this process [14–16]. This implies that DNMT1 may be indirectly involved in the mechanisms leading to FMR1 related disorders. Thus, expression of both genes, FMR1 and DNMT1, is of interest in fresh and archival samples from individuals with diseases characterized by abnormal expansions within the CGG region. Unfortunately, archival samples and samples extracted from blood have variable total RNA and mRNA degradation, which we have found to be a confounding factor in the gene expression analysis, unless a method was established to normalize for FMR1 and DNMT1 mRNA degradation.
Previous real-time PCR studies normalized gene expression against total cell number, total RNA or mRNA of internal control gene/s, reviewed in . The main constraints associated with using total cell number for normalization is that the RNA has to be of a very high quality and the cell counts have to be accurate, which is not possible for solid tissues. Total RNA has been particularly useful for normalization using absolute and relative standard curve methods for real-time PCR, and has been shown to produce biologically relevant and highly reproducible data . However, the most frequently applied normalization method involves the use of internal control gene/s based on the rational that the mRNA levels of these control genes in the tested samples should reflect the effect of variables on the quantity and quality of the target gene mRNA, and at the same time be minimally effected by the tested variables [19, 20]. Previous real-time PCR expression assays for FMR1 and DNMT1 have utilized GUS and HPRT as internal control genes [2, 21, 22]. However, these assays did not take into consideration the appropriateness of these genes as controls in relation to the quantity and quality of the starting material. A number of studies have suggested that confirmation of total RNA quality prior to real-time PCR quantification via capillary electrophoresis is essential in obtaining meaningful gene expression data .
In this study the Experion system (Bio-rad) was used to assess the quality of rRNA by subjective/visual interpretation of the chromatographs for absence of peaks other than those of 18S, 28S and 5S, and/or by examining the 28S:18S ratio and the RNA Qiality Indicator (RQI) values. However, these approaches have their limitations which we addressed in the study. We also evaluated in depth the use of internal control genes and the capillary electrophoresis features for the normalization of mRNA degradation in partially degraded and intact RNA from patient lymphoblasts and blood. We then defined the best predictor of target gene mRNA degradation, utilizing it for normalization in the whole blood RNA with a high variability of rRNA quality. In this sample set we examined the relationship between the FMR1 mRNA and the CGG expansion in a non-coding portion of this gene, and found that by incorporating the 3' and 5' mRNA analysis data from patients with small to intermediate CGG expansions, we have significantly improved upon the current approach used to examine the FMR1 mRNA toxicity, demonstrating the clinical relevance of our strategy.
Cell Line Details
MCRI Cell Line ID:
CGG repeat number
Healthy control male
Healthy control male
Healthy control male
Healthy control male
Healthy control male
Grey zone male
Grey zone male
Premutation female carrier
FRAX female carrier
Degradation of total RNA
Total RNA was extracted and purified using the RNeasy extraction kit (Qiagen Inc., Hilden, Germany). RNA concentrations were measured in triplicate using a NanoDrop ND-1000 Spectrophotometer, with purity being determined by the A260/A280 ratio using the expected values between 1.8 and 2. Each RNA sample was then diluted to 2 ng/ul. Sixty ul of each sample was left at room temperature for 0, 18, 24 and 96 hours. At each time point an aliquot of 15 ul was collected from each sample, 5 ul of which were used for capillary electrophoresis (Bio-rad) and 2 ul for cDNA synthesis.
RNA samples from blood of CGG expansion carriers
The original study that investigated the relationship between FMR1 expression and CGG expansion was approved by the La Trobe University Human Ethics Committee, the Southern Tasmanian Health and Medical Human Research Ethics Committees, and by the Institutional Review Board of the University of California at Davis. The methods used for CGG repeat sizing and FMR1 mRNA assessments, as well as the results of correlation between these two measures in a larger sample, which included the present participants, were reported in the earlier study . An aliquot of total RNA was originally isolated from 3 ml of peripheral blood using Tempus Blood RNA tubes as previously described in: Loesch, et al. (2007) . After more than 1 year of storage at -80°C, the total RNA quality was assessed using capillary electrophoresis. All samples were diluted to 5 ng/ul, prior to reverse transcription. Here we utilized these earlier data on CGG repeat size, but validated and conducted new assays on the expression of FMR1ex3.4, FMR1ex13.14 and GUS.
Reverse Transcription Real-time PCR
Reverse transcription was performed one reaction per sample using the Multiscribe Reverse Transcription System, 50 units/ul (Applied Biosystems). The 7900HT Fast Real Time PCR (Applied Biosystems) was used to quantify FMR1ex3.4, FMR1ex13.14, DNMT1, GAPDH, B2M, and GUS, using the relative standard curve method. The target gene and the internal control gene dynamic linear ranges were performed on a series of doubling dilutions of RNA (160-0.5 ng/ul) of a selected peripheral bood mononuclear cell sample. Previously published sequences were used for RT-PCR primers and TaqMan probe for: FMR1exon3/4 and GUS ; FMR1exon13/exon14 , and DNMT1 . FMR1exon3/4, FMR1exon13/14 and DNMT1 primers and probes were used at concentrations of 18 uM and 2 uM, respectively. GAPDH and B2M primer/probe mixes were obtained from PrimerDesign (PerfectProbe ge-PP-12-hu kit) and used at concentration of 2 uM. Each sample was assayed in duplicate 10 μl PCR reactions, consisting of 5.8 mM MgCl2, 1 μl Buffer A (Applied Biosystems), 3.35 μl RNase-free water, 1.2 mM dNTPs, 0.01 units/μl of AmpliTaq Gold, 0.5 μl of TaqMan probe and 0.5 μl forward and 0.5 μl reverse primers, and 1 μl of the reverse transcription (cDNA) reaction. The annealing temperature for thermal cycling protocol was 60°C for 40 cycles. Samples were quantified in arbitrary units (au) in relation to the standard curves performed on each plate.
Generalized estimating equations (GEE) were used for the objective assessment of the chromatographs and mRNA integrity. Normalization of FMR1 and DNMT1 mRNA to 18S, 28S, GAPDH and GUS, expressed as a function of the degradation time, was assessed by fitting the values to a simple linear regression. The sign test was then used to examine the null hypothesis that the median slope was zero against the two-sided alternative of it not being equal to 0. The relationship between the qPCR results for FMR1ex3.4 and FMR1ex13.14 assays with the CGG expansion size was assessed using a significance test in a linear regression. The analysis was conducted using the publicly available R statistical computing package, version 1.191 (R development Core team, 2004).
Subjective assessment of chromatographs and mRNA integrity of FMR1 and DNMT1
Sample 496 (Figure 2; Panel I) displayed minor differences in the capillary electrophoresis profiles between 0 and 24 hours, with the most prominent change being observed between 24 and 96 hours (a large increase in the inter-peak area and the fast region, and a marked decrease in the 28S peak). This was also reflected in the 28S:18S ratios and RQI values, where there was no change between 0 and 24 hours, but a large decrease between 24 and 96 hours. In contrast, FMR1ex13.14 mRNA level decreased by more than half within the first 18 hours and remained at this level for the following time points. A similar trend was observed for FMR1ex3.4 mRNA, which decreased between 18 and 24 hours of incubation, and remained at the same level through to 96 hours. In contrast, there were no major changes in DNMT1 mRNA quantities for sample 496 throughout the time course.
Sample 497 was different to all the other RNA examined, as it displayed a second, smaller 28S peak at 0 hours (Figure 2; Pannel II E). The origin of this peak, which has been previously reported using the analogous system (Schroeder et al.; United states patent publication 2006/0246577 A1) is unknown; the absence of any other RNA degradation markers such as an increase in the inter-peak area and/or the fast region argue against a plausible explanation that it could represent 28S rRNA degradation. It is also unclear why a significant change in the capillary elecrophoresis profile for this sample was observed at 18 hours, when the larger 28S peak(s) almost completely disappeared (Figure 1, Panel II F), while at 24 hours there was an increase in the 28S:18S ratio and RQI (Figure 1, Panel II F and G). This may be linked to the structure of the anomolous 28S peak, the unusual integrity of which is beyond the scope of this study. Interestingly, the mRNA levels for FMR1 and DNMT1 did not change between 0 and 18 hours. The greatest increase in total RNA degradation of sample 497 was observed between 24 and 96 hours, as indicated by a significant increase in the inter-peak area and the fast region and a marked decrease in the 28S:18S ratio and the RQI value (Figure 2, Panel II G and H). However, again this trend was not mirrored by the degradation of FMR1 and DNMT1 mRNA, which showed no major differences between 24 and 96 hours.
For the sample 488, within the first 18 hours there was no difference in RQI, however the 28S:18S ratio decreased from 2.5 to 1.57, with a slight increase in the inter-peak area within the first 18 hours (Figure 2, Panel III, I and J). At 24 and 96 hours we observed a further decrease in the 28S:18S ratio which was mirrored by a moderate decrease in RQI, and a prominent increase in the inter peak area and the fast region (Figure 2, Panel III, K and L). In contrast, there was only a slight decrease in DNMT1 mRNA level within the first 18 hours, which remained at the same level for the following time points.
Samples 584 and 585 demonstrated similar kinetics of total RNA degradation. At 0 hours these samples had similar RQI values, but different 28S:18S ratios (Figure 2, Panels IV and V). A striking increase in total RNA degradation was observed between 24 and 96 hours. During this period for sample 584 the 28S:18S ratio dropped from 1.2 to 0, and RQI from 9.3 to 2.9 (Figure 2, Panels IV; O and P), while for sample 585 the 28S:18S ratio dropped from 1.64 to 0.77, and RQI from 9.7 to 7.7 (Figure 2, Panel V; S and T). This was related to a significant increase in the inter peak area and the fast region. In contrast, mRNA quantities of FMR1ex3.4 and FMR1ex13.14 were decreased by approximately half between 0 and 18 hours in both samples 584 and 585, and for DNMT1 in sample 584 about 3 fold between 0 and 24 hours.
The similarities in FMR1 mRNA and total RNA degradation kinetics in both 584 and 585 samples may be related to these cell lines harbouring grey zone alleles, whereas samples 496 and 497 that showed different degradation kinetics were from premutation carriers (Table 1). Although, an in-depth investigation of this relationship is beyond the scope of this manuscript, the differences in FMR1ex3.4 mRNA degradation kinetics between the premutation and grey zone cell lines (Figure 2) were consistent with previous studies showing that the increased length of the CGG tract correlates with increased mRNA stability through hairpin formation within the 5'UTR repeat region .
Together these data suggested that the total RNA degradation rate moderately varies between the samples from different cell lines with the most prominent changes being observed between 24 and 96 hours. This was poorly correlated with the profile of FMR1 and DNMT1 mRNA degradation that predominantly occurred within the first 18 to 24 hours, indicating that there was no clear correlation between the rate of total RNA degradation from the subjective assessment of the chromatographs and mRNA degradation as determined by real-time PCR. Because FMR1 mRNA stability may be related to the size of the CGG repeat within its UTR and the pathology of FMR1 related disorders , in the following sections we have established a method to normalize for mRNA degradation independent of the CGG expansion size, so that the clinical relevance of the CGG related FMR1 mRNA toxicity can be identified in samples with variable rRNA quality.
Objective assessment of the chromatographs and mRNA integrity
Objective assessment of total RNA and mRNA degradation was performed in RNA samples from 14 cell lines (Table 1). Ten features were measured from each chromatograph and 6 variables measured using real-time PCR from the corresponding cDNA samples at 4 paired time points. The RNA from three fragile X cell lines were excluded from the FMR1 real-time PCR analysis for the RNA degradation study, as they had no FMR1 expression. The first aim of this approach was to objectively delineate which features of the chromatographs and gene expression profiles (GAPDH, B2M, GUS) could be used as predictors of the total RNA degradation as reflected by the degradation time. The second aim was to objectively delineate whether FMR1 and DNMT1 mRNA degradation correlated with the degradation time, and if not, which features could be used as predictors of the target gene mRNA integrity.
Based on the subjective assessment, we arbitrarily divided the degree of total RNA degradation into four categories: 0 hours – intact; 18 hours – early degradation; 24 hours – moderate to severe degradation; 96 hours – severe degradation, under the assumption that RNAs from different cell lines follow this progressive trend of degradation at the four time points. GEE were then utilised to provide an estimation of which parameters most closely reflected FMR1 and DNMT1 mRNA integrity in early versus moderate versus late degradation stages, and through the time course as a whole.
Correlation of RNA degradation time with capillary electrophoresis features and mRNA quantities
Predominant size of
RNA detected and
(for mRNA only)
0 to 18 hours comparison
0 to 24 hours comparison
0 to 96 hours comparison
18S % of total area (n = 53)
28S % of total area (n = 53)
5S % of total area (n = 53)
Fast region % of total area (n = 53)
Inter-peak region % of total area (n = 53)
Post region % of total area (n = 53)
Lower Marker % of total area (n = 53)
Small fragmentation region % of total area (n = 53)
28S:18S ratio (n = 53)
DNMT1 mRNA (au) (n = 41)
5084 bp – 142 bp
FMR1ex.3.4 mRNA (au) (n = 43)
4350 bp – 122 bp
FMR1ex13.14 mRNA (au) (n = 43)
4350 bp – 89 bp
GUS mRNA (au) (n = 52)
2200 bp – 80 bp
B2M mRNA (au) (n = 35)
976 bp – 114 bpo
GAPDH mRNA (au) (n = 52)
1310 bp – 110 bp
GUS and GAPDH mRNA degradation closely reflected moderate to severe degradation of total RNA (p < 0.001), best represented by the small fragmentation region % area. In contrast, B2M mRNA could only be used to predict early total RNA degradation (p = 0.093). Furthermore, it was not associated with any other chromatographic and gene expression parameters examined. For the target transcripts, DNMT1 and FMR1ex3/4, mRNA degradation was closely associated with moderate to severe degradation of total RNA (p < 0.001), and was best predicted by the small fragmentation region % area, and by GUS and GAPDH mRNA, all of which had similar degradation kinetics. In contrast, FMR1ex13.14 mRNA degradation appeared to be closely associated with early, moderate and severe total RNA degradation (p < 0.001), and was best predicted by % areas of 18S, 28S and the inter-peak region (p < 0.001). This analysis demonstrated that the degradation of mRNA was different between most internal controls and target genes examined, and that the degradation kinetics of specific mRNAs were not necessarily the same as those for total RNA and rRNA.
Assessment of FMR1 and DNMT1 mRNA normalized to capillary electrophoresis and real-time PCR predictors of mRNA degradation
Combining the analysis of 3' and 5' mRNA sites and GUS normalization minimizes the confounding impact of mRNA degradation in RNA samples from whole blood with highly variable rRNA quality
We have previously demonstrated using freshly extracted RNA, that FMR1 expression was significantly elevated in carriers of CGG expansion, compared with normal controls of a similar age, and that the expression was proportional to the size of CGG expansions within the grey zone and lower premutation range . Subsequent analysis of these stored samples revealed a high variability in rRNA quality, which posed major confounder concerns. We have therefore determined whether these samples could provide clinically relevant data using the normalization criteria tested in the study.
A number of samples (samples 339, 354, 358 and 360 – colour coded), did not follow the common pattern of correlation between the FMR1ex3.4/GUS and FMR1ex13.14/GUS assays (Figure 4A), and the FMR1ex3.4/GUS and FMR1ex13.14/GUS assays with the CGG expansion size (Figure 4B). Two of these samples, 339 and 360 were also outliers for rRNA quality assessment with the 28S:18S ratios of 3.3 and 0.24, and RQI values of 10 and 3, respectively (Figure 4D). Furthermore, there was no uniform correlation between rRNA integrity and FMR1 mRNA quality for both the FMR1ex3.4/GUS and FMR1ex13.14/GUS assays in these samples. Each of these samples appeared to be significantly affected at either the 5' or 3' sites (Figure 4A, B and 4C). Only one sample (360), with poor rRNA profile, appeared to have FMR1 mRNA integrity compromised at both the 5' and 3' sites (Figure 4B and 4D). The confounding impact of these outliers was minimized when the data for the FMR1 mRNA 3' and 5' end analyses were combined, as the significance of the FMR1 correlation with CGG expansion size increased to p = 0.018.
The specific aim of this study was to establish an optimal method to normalize for the degradation of target gene FMR1 mRNA in a sample set showing large variability in rRNA quality, and demonstrate it's clinical/biological relevance. Previous studies have normalized FMR1 expression to GUS [2, 21]. However, we questioned the appropriateness of GUS as an internal control for our sample set, where we have observed a large variation in rRNA degradation, particularly as there are no previous studies that examined the rates of mRNA degradation for both FMR1 and GUS. Since different mRNA species degrade depending on their length and secondary structure [18, 23, 29], we assessed if GUS was a suitable control for FMR1 mRNA degradation, and if not, which of the capillary electrophoresis and real-time PCR parameters would provide a better normalization method. Because of the potential relevance of DNMT1 to FMR1 gene regulation [30–33], we have included a parallel analysis of DNMT1 as as target gene in some of the studies.
We initially demonstrated using the Experion system that for the artificially degraded RNA samples, both the 28S:18S ratio and the RQI were most useful as predictors of severe RNA degradation, whereas the greatest changes in stability of different transcripts examined occurred during early to moderate stages of RNA degradation. Thus, the 28S:18S ratio and the RQI were not suitable predictors of mRNA stability, at least in our settings. Since, the RQI is closely related to a more widely used RNA Integrity Number (RIN) from an analogous Agilent system (Bio-Rad electrophoresis technical note 5761), our findings suggest that RIN may be also inappropriate as a normalization tool in our settings. Furthermore, the subjective assessment of general chromatographic features did not provide a useful estimation of mRNA degradation, as real-time PCR could still be used to obtain biologically relevant mRNA data in samples with chromatographs indicating severe rRNA degradation. In another approach, we determined that of the 10 selected chromatographic features, 18S, 28S and the inter-peak region % areas were the most reliable predictors of total RNA degradation when examined as a function of the degradation time. However, the normalization of the target genes to 18S and 28S chromatographic features was found to be inferior to the use of the internal control genes.
Primarily these observations indicate that the degradation kinetics of rRNA may be heavily size dependent, as the small 5S rRNA subunit, 160 nucleotides, was found to be a good predictor of only early to moderate RNA degradation, whereas the 18S and 28S % areas, 1770 and 3770 nucleotides respectively, could be used to predict RNA degradation at early, moderate and late stages of degradation. For mRNA quantitation, by qPCR, however, it has been previously suggested that the length of the amplicon, rather than that of the whole mRNA molecule, may be a more important indicator of degradation kinetics, particularly because fragmentation of a long mRNA may only result in a loss of the molecule for qPCR detection if the RNA break occurs within the product sequence .
We've found that the location of the amplicon, may be just as important as the size in determining the effect of mRNA degradation on qPCR performance. This was initially observed by examining the differences in FMR1ex3.4 (5') and FMR1ex13.14 (3') qPCR data throughout the RNA degradation time course. Since most FMR1 transcripts contain the 195-bp exon 14 , and all contain exons3/4, both assays target FMR1 transcripts of a similar size and abundance, which is experimentally reflected by the highly significant correlation between the two assays for high quality total RNA samples from the peripheral blood of patients with small to intermediate size expansions (Additional file 1). Since the amplicons for both assays are also of a similar size, if the location of the amplicon is unimportant in determining RNA degradation kinetics, its effects should be similar on the qPCR results for both assays. However there were clear differences between the 5' and 3' FMR1 results for the sample set with the high variability in total RNA quality, which may be explained by the 3' alternative splice site being more susceptible to RNA breaks than the 5' common region. Thus, thermal stress at the 3' site may result in the loss of more molecules for qPCR detection, than that at the 5' site, which is consistent with our findings in the RNA samples extracted from whole blood where RNA degradation had far less impact on the biological significance of the FMR1ex3.4 data than on the FMR1ex13.14 results (Figure 4B).
We have also found that due to the differences between degradation kinetics of mRNA and rRNA, normalization of target mRNA levels using rRNA profiles, as indicated by capillary electrophoresis, was inferior to the use of the internal control genes. Therefore, we proposed an alternative approach, independent of capillary electrophoresis, where by targeting each transcript at multiple sites for qPCR analysis, the confounding impact of mRNA breaks at any specific location on qPCR performance can be minimized. Specifically, in a previous study we've found that the FMR1 mRNA levels were proportional to the size of CGG expansions within the grey zone and lower PM range . In 30 of the original samples we've found the rRNA quality to vary significantly. Half of the samples would have had to be omitted from further analysis unless an alternative method to capillary electrophoresis was found to normalize for mRNA degradation.
We examined FMR1 expression in these samples without omission by targeting FMR1 transcript at 5' and 3' ends standardized to GUS, as we have shown that this strategy is superior to the use of rRNA profiles from the capillary electrophoresis analyses of the artificially degraded samples. Of the 30 samples we found four outliers that had a significant impact on qPCR performance particularly for the 3' FMR1ex13.14 assay. In 3 out of the 4 outliers this was not directly related to rRNA quality. By combining FMR1 mRNA 3' and 5' end analyses, the confounding impact of these outliers on the correlation between FMR1 expression and CGG size was minimized, and the statistical significance of the data doubled.
In summary, we have demonstrated that in artificially degraded RNA samples a number of chromatographic features including 18S, 28S, the inter-peak region, the 28S:18S ratio and RQI can be used as predictors of different stages of total RNA degradation. However, their use for normalization of target gene mRNA degradation was inferior to the use of internal control genes, of which GUS was the most appropriate as it closely reflected the target gene mRNA degradation kinetics. For the target FMR1 mRNA we've shown that the location of the amplicon, may be just as important as the size of its transcript in determining the effect of mRNA degradation on qPCR performance. Furthermore, we've found that by targeting the FMR1 transcript at multiple sites, the confounding impact of mRNA breaks due to fragmentation within any specific qPCR product sequence, and the subsequent loss of mRNA molecules for qPCR detection, was minimized. In clinical practice, this approach may be extremely useful, as there is extensive variation in the time lapse between the blood collection and RNA extraction for different samples, which contributes variation in rRNA and mRNA quality. Our strategy allows for a more accurate comparison between samples without sample loss.
Whilst we have shown that this approach is applicable for the FMR1 mRNA in RNA extracted from whole blood in our patient subset, there are also potential implications on the development of diagnostic tests for the levels of FMR1 mRNA toxicity associated with a number of clinical conditions [2–9, 12, 13]. These findings are also likely to have a broader application to expression studies of other genes using precious archival or fresh blood samples with a large variation in rRNA quality where sample omission is not an option. More importantly, the approach presented in this manuscript may be useful to any diagnostic application where mRNA integrity may be compromised , including future developments of methods for detection of placental mRNA in maternal plasma, foetal RNA markers for non-invasive prenatal diagnosis of pregnancy associated diseases and foetal chromosomal aneuploidies, as well as clinical uses of plasma or whole blood RNA for detection of tumour RNA in cancer diagnostics.
This study was supported by the March of Dimes Foundation, USA [6-FY06-318 to KHAC]; National Institute of Child Health and Human Development grant, USA [HD36071 to DZL]; Senior Principal Research Fellowship of the National Health and Medical Research Council of Australia [334300 to KHAC].
- Verkerk AJ, Pieretti M, Sutcliffe JS, Fu YH, Kuhl DP, Pizzuti A, Reiner O, Richards S, Victoria MF, Zhang FP: Identification of a gene (FMR-1) containing a CGG repeat coincident with a breakpoint cluster region exhibiting length variation in fragile X syndrome. Cell. 1991, 65 (5): 905-914. 10.1016/0092-8674(91)90397-H.View ArticlePubMedGoogle Scholar
- Tassone F, Hagerman RJ, Taylor AK, Gane LW, Godfrey TE, Hagerman PJ: Elevated levels of FMR1 mRNA in carrier males: a new mechanism of involvement in the fragile-X syndrome. Am J Hum Genet. 2000, 66 (1): 6-15. 10.1086/302720.View ArticlePubMedPubMed CentralGoogle Scholar
- Allen EG, He W, Yadav-Shah M, Sherman SL: A study of the distributional characteristics of FMR1 transcript levels in 238 individuals. Hum Genet. 2004, 114 (5): 439-447. 10.1007/s00439-004-1086-x.View ArticlePubMedGoogle Scholar
- Loesch DZ, Bui QM, Huggins RM, Mitchell RJ, Hagerman RJ, Tassone F: Transcript levels of the intermediate size or grey zone fragile X mental retardation 1 alleles are raised, and correlate with the number of CGG repeats. J Med Genet. 2007, 44 (3): 200-204. 10.1136/jmg.2006.043950.View ArticlePubMedGoogle Scholar
- Hagerman RJ, Leehey M, Heinrichs W, Tassone F, Wilson R, Hills J, Grigsby J, Gage B, Hagerman PJ: Intention tremor, parkinsonism, and generalized brain atrophy in male carriers of fragile X. Neurology. 2001, 57 (1): 127-130.View ArticlePubMedGoogle Scholar
- Jacquemont S, Orrico A, Galli L, Sahota PK, Brunberg JA, Anichini C, Leehey M, Schaeffer S, Hagerman RJ, Hagerman PJ: Spastic paraparesis, cerebellar ataxia, and intention tremor: a severe variant of FXTAS?. J Med Genet. 2005, 42 (2): e14-10.1136/jmg.2004.024190.View ArticlePubMedPubMed CentralGoogle Scholar
- Loesch DZ, Churchyard A, Brotchie P, Marot M, Tassone F: Evidence for, and a spectrum of, neurological involvement in carriers of the fragile X pre-mutation: FXTAS and beyond. Clin Genet. 2005, 67 (5): 412-417. 10.1111/j.1399-0004.2005.00425.x.View ArticlePubMedGoogle Scholar
- Greco CM, Hagerman RJ, Tassone F, Chudley AE, Del Bigio MR, Jacquemont S, Leehey M, Hagerman PJ: Neuronal intranuclear inclusions in a new cerebellar tremor/ataxia syndrome among fragile X carriers. Brain. 2002, 125 (Pt 8): 1760-1771. 10.1093/brain/awf184.View ArticlePubMedGoogle Scholar
- Tejada MI, Garcia-Alegria E, Bilbao A, Martinez-Bouzas C, Beristain E, Poch M, Ramos-Arroyo MA, Lopez B, Fernandez Carvajal I, Ribate MP: Analysis of the molecular parameters that could predict the risk of manifesting premature ovarian failure in female premutation carriers of fragile X syndrome. Menopause. 2008, 15 (5): 945-949. 10.1097/gme.0b013e3181647762.View ArticlePubMedGoogle Scholar
- Sofola OA, Jin P, Qin Y, Duan R, Liu H, de Haro M, Nelson DL, Botas J: RNA-binding proteins hnRNP A2/B1 and CUGBP1 suppress fragile X CGG premutation repeat-induced neurodegeneration in a Drosophila model of FXTAS. Neuron. 2007, 55 (4): 565-571. 10.1016/j.neuron.2007.07.021.View ArticlePubMedPubMed CentralGoogle Scholar
- Van Dam D, Errijgers V, Kooy RF, Willemsen R, Mientjes E, Oostra BA, De Deyn PP: Cognitive decline, neuromotor and behavioural disturbances in a mouse model for fragile-X-associated tremor/ataxia syndrome (FXTAS). Behav Brain Res. 2005, 162 (2): 233-239. 10.1016/j.bbr.2005.03.007.View ArticlePubMedGoogle Scholar
- Cohen S, Masyn K, Adams J, Hessl D, Rivera S, Tassone F, Brunberg J, DeCarli C, Zhang L, Cogswell J: Molecular and imaging correlates of the fragile X-associated tremor/ataxia syndrome. Neurology. 2006, 67 (8): 1426-1431. 10.1212/01.wnl.0000239837.57475.3a.View ArticlePubMedGoogle Scholar
- Bacalman S, Farzin F, Bourgeois JA, Cogswell J, Goodlin-Jones BL, Gane LW, Grigsby J, Leehey MA, Tassone F, Hagerman RJ: Psychiatric phenotype of the fragile X-associated tremor/ataxia syndrome (FXTAS) in males: newly described fronto-subcortical dementia. J Clin Psychiatry. 2006, 67 (1): 87-94.View ArticlePubMedGoogle Scholar
- Schermelleh L, Haemmer A, Spada F, Rosing N, Meilinger D, Rothbauer U, Cardoso MC, Leonhardt H: Dynamics of Dnmt1 interaction with the replication machinery and its role in postreplicative maintenance of DNA methylation. Nucleic Acids Res. 2007, 35 (13): 4301-4312. 10.1093/nar/gkm432.View ArticlePubMedPubMed CentralGoogle Scholar
- Takebayashi S, Tamura T, Matsuoka C, Okano M: Major and essential role for the DNA methylation mark in mouse embryogenesis and stable association of DNMT1 with newly replicated regions. Mol Cell Biol. 2007, 27 (23): 8243-8258. 10.1128/MCB.00899-07.View ArticlePubMedPubMed CentralGoogle Scholar
- Brylawski BP, Chastain PD, Cohen SM, Cordeiro-Stone M, Kaufman DG: Mapping of an origin of DNA replication in the promoter of fragile X gene FMR1. Exp Mol Pathol. 2007, 82 (2): 190-196. 10.1016/j.yexmp.2006.10.004.View ArticlePubMedGoogle Scholar
- Bustin SA, Benes V, Nolan T, Pfaffl MW: Quantitative real-time RT-PCR – a perspective. J Mol Endocrinol. 2005, 34 (3): 597-601. 10.1677/jme.1.01755.View ArticlePubMedGoogle Scholar
- Bustin SA: Quantification of mRNA using real-time reverse transcription PCR (RT-PCR): trends and problems. J Mol Endocrinol. 2002, 29 (1): 23-39. 10.1677/jme.0.0290023.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedPubMed CentralGoogle Scholar
- Pietrobono R, Tabolacci E, Zalfa F, Zito I, Terracciano A, Moscato U, Bagni C, Oostra B, Chiurazzi P, Neri G: Molecular dissection of the events leading to inactivation of the FMR1 gene. Hum Mol Genet. 2005, 14 (2): 267-277. 10.1093/hmg/ddi024.View ArticlePubMedGoogle Scholar
- Langer F, Dingemann J, Kreipe H, Lehmann U: Up-regulation of DNA methyltransferases DNMT1, 3A, and 3B in myelodysplastic syndrome. Leuk Res. 2005, 29 (3): 325-329. 10.1016/j.leukres.2004.08.004.View ArticlePubMedGoogle Scholar
- Fleige S, Walf V, Huch S, Prgomet C, Sehm J, Pfaffl MW: Comparison of relative mRNA quantification models and the impact of RNA integrity in quantitative real-time RT-PCR. Biotechnol Lett. 2006, 28 (19): 1601-1613. 10.1007/s10529-006-9127-2.View ArticlePubMedGoogle Scholar
- Maddalena A, Richards CS, McGinniss MJ, Brothman A, Desnick RJ, Grier RE, Hirsch B, Jacky P, McDowell GA, Popovich B: Technical standards and guidelines for fragile X: the first of a series of disease-specific supplements to the Standards and Guidelines for Clinical Genetics Laboratories of the American College of Medical Genetics. Quality Assurance Subcommittee of the Laboratory Practice Committee. Genet Med. 2001, 3 (3): 200-205.View ArticlePubMedPubMed CentralGoogle Scholar
- Fu YH, Kuhl DP, Pizzuti A, Pieretti M, Sutcliffe JS, Richards S, Verkerk AJ, Holden JJ, Fenwick RG, Warren ST: Variation of the CGG repeat at the fragile X site results in genetic instability: resolution of the Sherman paradox. Cell. 1991, 67 (6): 1047-1058. 10.1016/0092-8674(91)90283-5.View ArticlePubMedGoogle Scholar
- Heitz D, Devys D, Imbert G, Kretz C, Mandel JL: Inheritance of the fragile X syndrome: size of the fragile X premutation is a major determinant of the transition to full mutation. J Med Genet. 1992, 29 (11): 794-801. 10.1136/jmg.29.11.794.View ArticlePubMedPubMed CentralGoogle Scholar
- Fleige S, Pfaffl MW: RNA integrity and the effect on the real-time qRT-PCR performance. Mol Aspects Med. 2006, 27 (2–3): 126-139. 10.1016/j.mam.2005.12.003.View ArticlePubMedGoogle Scholar
- Napierala M, Michalowski D, de Mezer M, Krzyzosiak WJ: Facile FMR1 mRNA structure regulation by interruptions in CGG repeats. Nucleic Acids Res. 2005, 33 (2): 451-463. 10.1093/nar/gki186.View ArticlePubMedPubMed CentralGoogle Scholar
- Galzitskaya O, Finkelstein AV: Folding rate dependence on the chain length of RNA-like heteropolymers. Fold Des. 1998, 3 (2): 69-78. 10.1016/S1359-0278(98)00012-1.View ArticlePubMedGoogle Scholar
- Chiurazzi P, Pomponi MG, Pietrobono R, Bakker CE, Neri G, Oostra BA: Synergistic effect of histone hyperacetylation and DNA demethylation in the reactivation of the FMR1 gene. Hum Mol Genet. 1999, 8 (12): 2317-2323. 10.1093/hmg/8.12.2317.View ArticlePubMedGoogle Scholar
- Wohrle D, Salat U, Hameister H, Vogel W, Steinbach P: Demethylation, reactivation, and destabilization of human fragile X full-mutation alleles in mouse embryocarcinoma cells. Am J Hum Genet. 2001, 69 (3): 504-515. 10.1086/322739.View ArticlePubMedPubMed CentralGoogle Scholar
- Pietrobono R, Pomponi MG, Tabolacci E, Oostra B, Chiurazzi P, Neri G: Quantitative analysis of DNA demethylation and transcriptional reactivation of the FMR1 gene in fragile X cells treated with 5-azadeoxycytidine. Nucleic Acids Res. 2002, 30 (14): 3278-3285. 10.1093/nar/gkf434.View ArticlePubMedPubMed CentralGoogle Scholar
- Tabolacci E, Pietrobono R, Moscato U, Oostra BA, Chiurazzi P, Neri G: Differential epigenetic modifications in the FMR1 gene of the fragile X syndrome after reactivating pharmacological treatments. Eur J Hum Genet. 2005, 13 (5): 641-648. 10.1038/sj.ejhg.5201393.View ArticlePubMedGoogle Scholar
- Huang T, Li LY, Shen Y, Qin XB, Pang ZL, Wu GY: Alternative splicing of the FMR1 gene in human fetal brain neurons. Am J Med Genet. 1996, 64 (2): 252-255. 10.1002/(SICI)1096-8628(19960809)64:2<252::AID-AJMG3>3.0.CO;2-W.View ArticlePubMedGoogle Scholar
- Wong BC, Chiu RW, Tsui NB, Chan KC, Chan LW, Lau TK, Leung TN, Lo YM: Circulating placental RNA in maternal plasma is associated with a preponderance of 5' mRNA fragments: implications for noninvasive prenatal diagnosis and monitoring. Clin Chem. 2005, 51 (10): 1786-1795. 10.1373/clinchem.2005.052340.View ArticlePubMedGoogle Scholar
- The pre-publication history for this paper can be accessed here:http://www.biomedcentral.com/1472-6890/9/5/prepub
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.