Applicability of RNA standards for evaluating RT-qPCR assays and platforms
© Devonshire et al; licensee BioMed Central Ltd. 2011
Received: 13 October 2010
Accepted: 18 February 2011
Published: 18 February 2011
The availability of diverse RT-qPCR assay formats and technologies hinder comparability of data between platforms. Reference standards to facilitate platform evaluation and comparability are needed. We have explored using universal RNA standards for comparing the performance of a novel qPCR platform (Fluidigm® BioMark™) against the widely used ABI 7900HT system. Our results show that such standards may form part of a toolkit to evaluate the key performance characteristics of platforms.
Reverse transcription quantitative PCR (RT-qPCR) is rapidly becoming a valuable tool for mRNA biomarker quantification in clinical diagnostics. There has been a proliferation of RT-qPCR assay formats and platforms in recent years due to wider applications of this technology, coupled with improvements in sensitivity, specificity and accuracy of measurements of gene expression. However, there is one intrinsic limitation to the current qPCR platforms, namely lack of controls for cross-platform comparison. Although the manufacturers have developed platform-specific quality controls, they are often not adequate for cross-platform comparisons, particularly for the evaluation and standardization of transcriptomic data due to differences in protocols, data processing and analysis methods. Thus, development of universal RNA standards offers great potential in the validation of data obtained from different RT-qPCR methods. In the present investigation, we have compared the performance of Fluidigm® BioMark™ Integrated microfluidic (henceforth referred to as BioMark) dynamic arrays with the widely used ABI 7900HT real-time PCR platform (henceforth called ABI 7900HT system) using generic RNA standards.
Pre-amplification of RNA or cDNA facilitates the investigation of a large number of genes when the starting material is limiting, such as with tissue biopsies and archival formalin-fixed paraffin-embedded (FFPE) samples [1, 2]. Pre-amplification methods used generally include either linear amplification of RNA or exponential (PCR-based) amplification of cDNA [3–5]. However, concerns have been raised as to whether pre-amplification of samples by exponential amplification introduces bias in expression levels between genes . For the BioMark microfluidic PCR system, each sample in the 48 × 48 dynamic array is distributed amongst 48 different reaction chambers, therefore pre-amplification is recommended for certain applications. However the limit of detection (LOD) of using pre-amplified vs. non-amplified cDNA samples, and its impact on the technical performance of the PCR array have not been fully characterized.
Exogenous RNA controls produced by in vitro transcription are ideal materials for investigating different RT-qPCR kits and methodologies . Recently a panel of RNA controls have been developed for use in gene expression applications by the External RNA Controls Consortium (ERCC), an ad hoc group of 70 members from private, public and academic organizations led by the National Institute of Standards (NIST) [8, 9]. It is hoped that standards developed from these sequences will aid in comparisons of gene expression data generated from various platforms such as microarray, RT-qPCR and next generation sequencing, and also provide quality control of gene expression measurements in the clinical laboratory . Multigene biomarker measurements are at the forefront of a new class of medical devices using in vitro diagnostic multivariate assays, such as MammaPrint and Oncotype Dx in the area of breast cancer prognosis . Since gene expression biomarkers typically encompass a range of transcript abundances and differential expression ratios, it is more appropriate to use multiple RNA standards as quality controls for standardizing such measurements, as opposed to a single transcript at a fixed concentration.
In the current study, we used a sub-set of the 96 ERCC RNA standards (Additional File 1) in order to characterize their performance on a nanofluidic PCR system, the BioMark 48 × 48 dynamic arrays, against a conventional qPCR platform, the ABI 7900HT system. We also investigated the impact of pre-amplification of cDNA samples on the linear range and precision of measurements by nanofluidic qPCR. Two prototype panels were constructed with selected RNA standards containing varying copy number within each panel, and varying ratios between them for mimicking non-differentially and differentially expressed mRNA biomarkers as represented in normal and disease states. The expression profile of the RNA standards was measured using both platforms and the accuracy and precision of their detection were compared.
Linear range of dynamic arrays
One advantage of the BioMark arrays is the capability to analyse a large number of genes in a single sample. In order to facilitate this, up to ~ 2 μL of sample is loaded into each sample inlet of the chip and further distributed in the channels of the microfluidic chip as 48 separate 9 nL reactions using the integrated fluidic circuit (IFC). Thus the original sample is diluted more than 200-fold prior to the PCR reaction. In order to ensure that there are sufficient copies of target molecules in each reaction, Fluidigm® recommends using either RNA samples that do not have a concentration lower than 250 ng total RNA/μL or that a pre-amplification stage is included, whereby the cDNA sample undergoes 14-18 cycles of amplification with a mix of up to 100 different primer pairs (Fluidigm Advanced Development Protocols 3, 5 and 8). In order to further investigate the requirement for pre-amplification, RNA standards were spiked into human total RNA at different concentrations (for sample composition, see Additional File 2) with the aim of mimicking a range of physiological abundances, from highly abundant mRNA transcripts (106 copies/ng total RNA; equivalent to 104 copies per cell) to transcripts only expressed in a sub-population of cells (1 copy/ng total RNA; equivalent to 0.01 copies per cell), based on the RNA content of a cell estimated as 26 pg . A single RT reaction was performed for each RNA sample followed by 3 independent qPCR runs, with replicate assay measurements for each ERCC standard.
LOD of dynamic arrays
qPCR accuracy and precision
Accuracy and precision of linear detection range of PCR platforms
RNA biomarker panels
Concentrations and ratios of ERCC RNA standards in simulated 'normal' and 'disease' panels
Copies/ng total RNA
1 × 105
1 × 105
1 × 102
1.5 × 102
1 × 104
5 × 103
5 × 100
1 × 102
1 × 102
1 × 102
1 × 102
5 × 102
1 × 103
1 × 103
8 × 103
1.2 × 104
1 × 101
1 × 101
1 × 101
1 × 102
Overall, fold change estimation was found to be accurate for both ABI 7900HT system and BioMark arrays (Figures 4A and 4B), with the observed fold change values overlapping with the expected fold change measurements for all standards. For ERCC standards mimicking low abundance transcripts, the technical noise associated with the resulting fold change measurement was considerably greater than higher abundance RNA species. For example, a 20-fold increase in expression level at 5 copies/ng (ERCC-51) is associated with 16-fold and 10-fold difference in the minimum and maximum fold changes detected by the BioMark arrays and ABI 7900HT system respectively. For non-differential expression at low copy numbers (fold change = 1.0; ERCC-113), fold change measurements spanned a range of over 50% and 150% of the expected value on the BioMark and ABI 7900HT system respectively. At levels of abundance exceeding 100 copies/ng, mean fold change measurements were accurate to within 10% of expected values.
In this study we sought to demonstrate the utility of RNA standards for characterisation of a new platform, the BioMark, where PCR reactions are performed in volumes over a 1000-fold lower than on a conventional RT-qPCR instrument (ABI 7900HT system). The requirement for sample pre-amplification for this technology contrasts with standard two-step RT-qPCR approach, therefore the impact of this methodology was also evaluated. Dilutions of RNA standards across a wide physiological range demonstrated that the linear detection range of the BioMark arrays is similar to the ABI 7900HT real-time PCR system, when pre-amplified cDNA is used as the template (Figure 1). The precision of replicate measurements within the array also compared favourably to the intra-run standard deviation of Ct values for the ABI 7900HT system (Figure 3). At copy numbers mimicking medium to high abundance transcripts, the precision of the BioMark arrays is in a similar range to the minimum variation of ~ 0.1 units observed for another nanofluidic PCR array, the OpenArray format .
However, when non-amplified cDNA is quantified using the nanofluidic BioMark arrays, the linear range is severely limited, to only two orders of magnitude (Figure 1). At the lower detection limit of 104 RNA copies per reaction (Ct ~ 27), the variation between measurements increases significantly (Figure 3), whilst below this level of abundance, the rate of PCR failures increases rapidly (Figure 2). Pre-amplification of template cDNA using a preliminary PCR step of 14 cycles improves both the accuracy and precision of the transcript quantification using the dynamic arrays (Table 1). The improved detection of the 10-fold differences in RNA copy numbers between sample series (resulting in an average slope within 6% of expected value) also indicates that the pre-amplification process does not introduce bias into the detection of transcripts which cover a wide dynamic range.
Relative expression measurements are central to gene expression analysis by RT-qPCR and for determining whether a panel of biomarkers has predictive power for disease diagnosis and prognosis . Therefore, we developed two panels of RNA standards in order to further investigate the accuracy of detecting gene expression ratios using the new type of PCR array compared to an established system. The standards were spiked at varying ratios between panels in order to obtain information on how well the methodologies can discriminate between differentially and non-differentially expressed candidate genes at different transcript abundance levels (Figure 4). Our results show good accuracy of observed vs. expected values for both platforms, which is in agreement with previous studies demonstrating good concordance of fold change measurements between the BioMark arrays and the ABI 7900HT system . The precision of the fold change estimation varied according to the abundance of the transcripts, demonstrating increased variation in the observed values for lower concentrations of standards for both nanofluidic and standard real-time PCR approaches. This suggests that the sensitivity of the technique to correctly detect the expected fold change is reduced at low copy numbers (10 RNA copies or less per reaction on the ABI 7900HT system). Dixon et al.  also found that the sensitivity of the OpenArray platform was lower for Ct values corresponding to lower copy numbers, resulting in an increased number of false negative results. The increased variation in fold change detection at low copy numbers is likely to arise due to decreased efficiency at RT stage and increased stochastic variation in the PCR reaction for low target numbers .
We also found that both qPCR platforms were able to accurately detect a 1.5-fold change in mRNA expression, below the 2-fold cut-off which has been cited as a limit to the resolving power of conventional PCR, as it constitutes a difference of less than a single cycle . The BioMark dynamic arrays were recently shown to be able to detect a 1.25-fold difference in DNA copy number by qPCR, with greater levels of precision achievable with the larger number of technical replicates possible with this high-throughput approach . Since assay and sample loadings are in separate inlets on 48 × 48 dynamic arrays, it is possible to increase technical replication by using multiple assay inlets and/or multiple sample inlets. However, it should be noted that replication only at the assay level does not substitute for true sampling variation by the process of taking a sample from a population of molecules.
The use of gene-specific oligonucleotide standards for inter-run and cross-platform calibration has been demonstrated to improve the accuracy of class prediction based on panels of biomarkers . Although ERCC RNA standards do not directly provide information on the performance of biomarker-specific assays, a panel of multiple standards, such as that used here provides a robust means of evaluating platform performance by minimizing confounding effects resulting from differences in assay performance due to individual primer and probe specificity. RNA standards could also serve as calibrator samples between experiments where different sets of potential biomarkers genes are investigated, as well as in the context of a diagnostic assay where the expression of the same panel of genes is quantified. In addition to target gene normalization using a reference gene or panel of reference genes , normalization to an ERCC RNA standard or multiple RNA standards may be a useful control for elucidating technical variation due to RT and qPCR steps .
We conclude that universal RNA standards can provide robust information on the performance characteristics of different RT-qPCR platforms and methodologies. The results obtained using panels of multiple RNA standards indicate that the linear detection range, precision and accuracy of nanofluidic BioMark dynamic arrays are similar to those of an established real-time PCR instrument, the ABI 7900HT system, when pre-amplified cDNA is used as the template. The standards also provide reference values for the range of transcript abundance over which it would be possible to measure non-amplified cDNA on the nanofluidic BioMark high-throughput arrays. Carefully constructed panels of ERCC RNA standards have the potential to act as benchmarks for the calibration and interpretation of biomarker measurements in drug discovery and clinical diagnostics. Further evaluation of these standards is required for potential incorporation into a 'quality metrics' toolkit for assessing their suitability for cross-platform comparisons.
Preparation of in vitro transcribed RNA and samples
In vitro transcribed ERCC RNA standards were produced from ERCC plasmid DNA (courtesy of Dr. Marc Salit, NIST, USA). Plasmid DNA from standards ERCC-13, 25, 42, 51, 81, 84, 95, 99, 113 and 171 was cleaved into a single linear molecule using Bam HI restriction endonuclease (New England Biolabs, UK). 500 ng of plasmid DNA was used for each sample and digested by adding 40 U of Bam HI enzyme in NEB3 buffer provided by the manufacturer. The digestion mixture was incubated at 37°C for 2 hours followed by purification using QiaQuick PCR purification kit with an elution volume of 32 μl. In vitro transcription was carried out with 8 μl digested plasmid DNA using MEGAscript® T7 Kit (Applied Biosystems/Ambion, UK) followed by DNase treatment and clean-up using RNeasy columns (Qiagen, UK). RNA concentration and insert sizes were estimated using the Nanodrop 1000 spectrophotometer (Thermo Scientific, UK) and 2100 Bioanalyzer (Agilent Technologies, USA) respectively. RNA standards were diluted in nuclease free-water and spiked into Universal Human Reference RNA (Stratagene, UK) (final concentration 100 ng/μl). For experiments investigating the linear range of platform detection, standards were spiked at 10-fold intervals between 1 and 106 copies/ng total RNA (Additional File 2). For the simulated 'normal' and 'disease' panels, standards were spiked at various copy numbers and ratios (Table 2).
Reverse transcription and pre-amplification of cDNA
RNA samples were reverse-transcribed using the TaqMan® Reverse Transcription Reagents kit (Applied Biosystems, UK) in 40 μL reactions containing 400 ng total RNA and oligo(dT) primers according to manufacturer's instructions. cDNA samples were diluted to a concentration of 0.5 ng/μL (total RNA equivalent) with nuclease-free water. For experiments investigating the linear range of platform detection (Figures 1, 2, 3), a single RT reaction was performed for each RNA sample whilst for the simulated 'normal' and 'disease' panels (Figure 4), 6 replicate RT reactions were performed. A single aliquot of each cDNA sample, equivalent to 12.5 ng RNA, was pre-amplified with assays corresponding to all 10 standards in a 25 μL volume reaction using TaqMan® PreAmp Mastermix (Applied Biosystems, UK) according to manufacturer's protocol. Following pre-amplification, the samples were diluted 1:5 (v/v) in TE buffer, pH 8.0.
Further information on sample preparation and real-time PCR validation complying with the Minimum Information for Publication of Quantitative Real-Time PCR Experiments (MIQE) guidelines  is available in Additional Files 2 and 3 (MIQE Additional Information and Checklist).
Custom-designed primers and TaqMan® FAM-TAMRA probes for each ERCC standard (Additional File 1) were supplied by Applied Biosystems and a 20 × assay mix was prepared containing 18 μM primer and 5 μM probe (final concentration 900 nM primer and 250 nM probe). qPCR assays were tested initially using a serial dilution of ERCC cDNA and PCR efficiencies calculated (see Additional Data Table 2: MIQE Additional information). All 10 assays were found to have PCR efficiencies of greater than 86%.
BioMark arrays were prepared according to the manufacturer's instructions. TaqMan® assays were diluted 1:1 (v/v) with DA Assay Loading Reagent (Fluidigm®) and 5 μL was added to each assay inlet of the array. Also, 5 μL reaction mix was prepared by mixing 2 × TaqMan® Universal Mastermix (Applied Biosystems), DA Sample Loading Reagent and nuclease-free water containing 2 μL of cDNA or pre-amplified cDNA. The samples were loaded into each sample inlet as per manufacturer's recommendations. Following loading of the assays and samples into the chip by the IFC controller, PCR was performed with the following reactions conditions: 50°C for 2 minutes, 95°C for 10 minutes, followed by 45 cycles of 95°C for 15 seconds and 60°C for 60 seconds. Data was processed by automatic global threshold setting with the same threshold value for all assays and linear baseline correction using BioMark Real-time PCR Analysis software (version 2.1.1). The quality threshold was set at the default setting of 0.65. For experiments investigating the linear range of platform detection (Figures 1, 2, 3), 8 qPCR reactions consisting of 4 assay inlet and 2 sample inlet replicates were performed for each cDNA or pre-amplified cDNA sample. For the simulated 'normal' and 'disease' panels (Figure 4), 12 qPCR reactions consisting of 4 assay inlet and 3 sample inlet replicates were performed for each cDNA sample.
Conventional real-time PCR was performed using ABI 7900HT system in 20 μL reaction volumes containing TaqMan® Universal PCR Master Mix and 2 μL of respective cDNA in optical 96-well plates (Applied Biosystems). Cycling conditions were as those used for the BioMark arrays. Triplicate qPCR reactions were performed for each cDNA sample for all experiments. The threshold fluorescence level was set manually for each plate using SDS software version 2.3 (Applied Biosystems). Following export of Cycle threshold (Ct) data, further data analysis for both platforms was performed in Microsoft® Excel 2003. Comparison of slope and R2 values between pre-amplified and non-amplified cDNA, as a template on the BioMark arrays, was performed as paired t-test in Microsoft® Excel 2003.
The following additional are available with the online version of this paper. Additional data file 1 is a table detailing the primer and probe sequences used for qPCR assays. Additional files 2 and 3 are additional data and a checklist in compliance with the MIQE (Minimum Information for Publication of Quantitative Real-Time PCR Experiments) guidelines.
List of abbreviations
External RNA Controls Consortium
Reverse Transcription Quantitative PCR
limit of detection
integrated fluidic circuit
polymerase chain reaction.
The work described in this paper was funded by the UK National Measurement System. We are grateful to Dr. Marc Salit (NIST, USA) for the provision of ERCC plasmid DNA. We would also like to thank Jesus Minguez (LGC) for statistical analysis and Dr. Bridget Fox (LGC) for assistance with the production of in vitro transcribed RNA standards.
- Vermeulen J, Derveaux S, Lefever S, De Smet E, De Preter K, Yigit N, De Paepe A, Pattyn F, Speleman F, Vandesompele J: RNA pre-amplification enables large-scale RT-qPCR gene-expression studies on limiting sample amounts. BMC Res Notes. 2009, 2: 235-10.1186/1756-0500-2-235.PubMed CentralPubMedView ArticleGoogle Scholar
- Li J, Smyth P, Cahill S, Denning K, Flavin R, Aherne S, Pirotta M, Guenther SM, O'Leary JJ, Sheils O: Improved RNA quality and TaqMan Pre-amplification method (PreAmp) to enhance expression analysis from formalin fixed paraffin embedded (FFPE) materials. BMC Biotechnol. 2008, 8: 10-10.1186/1472-6750-8-10.PubMed CentralPubMedView ArticleGoogle Scholar
- Eberwine J, Yeh H, Miyashiro K, Cao Y, Nair S, Finnell R, Zettel M, Coleman P: Analysis of gene expression in single live neurons. Proc Natl Acad Sci USA. 1992, 89: 3010-3014. 10.1073/pnas.89.7.3010.PubMed CentralPubMedView ArticleGoogle Scholar
- Kurn N, Chen P, Heath JD, Kopf-Sill A, Stephens KM, Wang S: Novel isothermal, linear nucleic acid amplification systems for highly multiplexed applications. Clin Chem. 2005, 51: 1973-1981. 10.1373/clinchem.2005.053694.PubMedView ArticleGoogle Scholar
- Iscove NN, Barbara M, Gu M, Gibson M, Modi C, Winegarden N: Representation is faithfully preserved in global cDNA amplified exponentially from sub-picogram quantities of mRNA. Nat Biotechnol. 2002, 20: 940-943. 10.1038/nbt729.PubMedView ArticleGoogle Scholar
- Ginsberg SD: RNA amplification strategies for small sample populations. Methods. 2005, 37: 229-237. 10.1016/j.ymeth.2005.09.003.PubMedView ArticleGoogle Scholar
- Levesque-Sergerie JP, Duquette M, Thibault C, Delbecchi L, Bissonnette N: Detection limits of several commercial reverse transcriptase enzymes: impact on the low- and high-abundance transcript levels assessed by quantitative RT-PCR. BMC Mol Biol. 2007, 8: 93-10.1186/1471-2199-8-93.PubMed CentralPubMedView ArticleGoogle Scholar
- External RNA Controls Consortium: Proposed methods for testing and selecting the ERCC external RNA controls. BMC Genomics. 2005, 6: 150-10.1186/1471-2164-6-150.PubMed CentralView ArticleGoogle Scholar
- Baker SC, Bauer SR, Beyer RP, Brenton JD, Bromley B, Burrill J, Causton H, Conley MP, Elespuru R, Fero M, Foy C, Fuscoe J, Gao X, Gerhold DL, Gilles P, Goodsaid F, Guo X, Hackett J, Hockett RD, Ikonomi P, Irizarry RA, Kawasaki ES, Kaysser-Kranich T, Kerr K, Kiser G, Koch WH, Lee KY, Liu C, Liu ZL, Lucas A, et al: The External RNA Controls Consortium: a progress report. Nat Methods. 2005, 2: 731-734. 10.1038/nmeth1005-731.PubMedView ArticleGoogle Scholar
- Cronin M, Ghosh K, Sistare F, Quackenbush J, Vilker V, O'Connell C: Universal RNA Reference Materials for Gene Expression. Clin Chem. 2004, 50: 1464-1471. 10.1373/clinchem.2004.035675.PubMedView ArticleGoogle Scholar
- Ross JS: Multigene classifiers, prognostic factors, and predictors of breast cancer clinical outcome. Adv Anat Pathol. 2009, 16: 204-215. 10.1097/PAP.0b013e3181a9d4bf.PubMedView ArticleGoogle Scholar
- Clinical and Laboratory Standards Institute: Use of External RNA Controls in Gene Expression Assays. 2010, Approved Guideline 26 No. 29, MM16-A.Google Scholar
- Stahlberg A, Hakansson J, Xian X, Semb H, Kubista M: Properties of the reverse transcription reaction in mRNA quantification. Clin Chem. 2004, 50: 509-15. 10.1373/clinchem.2003.026161.PubMedView ArticleGoogle Scholar
- Dixon JM, Lubomirski M, Amaratunga D, Morrison TB, Brenan CJ, Ilyin SE: Nanoliter high-throughput RT-qPCR: a statistical analysis and assessment. Biotechniques. 2009, 46: ii-viii. 10.2144/000112746.PubMedView ArticleGoogle Scholar
- Vermeulen J, Pattyn F, De Preter K, Vercruysse L, Derveaux S, Mestdagh P, Lefever S, Hellemans J, Speleman F, Vandesompele J: External oligonucleotide standards enable cross laboratory comparison and exchange of real-time quantitative PCR data. Nucleic Acids Res. 2009, 37: e138-10.1093/nar/gkp721.PubMed CentralPubMedView ArticleGoogle Scholar
- Spurgeon SL, Jones RC, Ramakrishnan R: High throughput gene expression measurement with real time PCR in a microfluidic dynamic array. PLoS One. 2008, 3: e1662-10.1371/journal.pone.0001662.PubMed CentralPubMedView ArticleGoogle Scholar
- Bustin SA, Nolan T: Pitfalls of quantitative real-time reverse-transcription polymerase chain reaction. J Biomol Tech. 2004, 15: 155-166.PubMed CentralPubMedGoogle Scholar
- Lo YM, Lun FM, Chan KC, Tsui NB, Chong KC, Lau TK, Leung TY, Zee BC, Cantor CR, Chiu RW: Digital PCR for the molecular detection of fetal chromosomal aneuploidy. Proc Natl Acad Sci USA. 2007, 104: 13116-13121. 10.1073/pnas.0705765104.PubMed CentralPubMedView ArticleGoogle Scholar
- Weaver S, Dube S, Mir A, Qin J, Sun G, Ramakrishnan R, Jones RC, Livak KJ: Taking qPCR to a higher level: Analysis of CNV reveals the power of high throughput qPCR to enhance quantitative resolution. Methods. 2010, 50: 271-6. 10.1016/j.ymeth.2010.01.003.PubMedView ArticleGoogle 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: RESEARCH0034-10.1186/gb-2002-3-7-research0034.PubMed CentralPubMedView ArticleGoogle Scholar
- Gilsbach R, Kouta M, Bonisch H, Bruss M: Comparison of in vitro and in vivo reference genes for internal standardization of real-time PCR data. Biotechniques. 2006, 40: 173-177. 10.2144/000112052.PubMedView ArticleGoogle Scholar
- Bustin SA, Benes V, Garson JA, Hellemans J, Huggett J, Kubista M, Mueller R, Nolan T, Pfaffl MW, Shipley GL, Vandesompele J, Wittwer CT: The MIQE Guidelines: Minimum Information for Publication of Quantitative Real-Time PCR Experiments. Clin Chem. 2009, 55: 611-622. 10.1373/clinchem.2008.112797.PubMedView ArticleGoogle Scholar
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.