Characterization of Globin RNA Interference in Gene Expression Profiling of Whole-Blood Samples

How long does the mRNA remains stable in untreated whole bovine blood?

BACKGROUND

High quality and quantity of messenger RNA (mRNA) are required for accuracy of gene expression studies and other RNA-based downstream applications. Since RNA is considered a labile macromolecular prone to degradation, which may result in falsely altered gene expression patterns, several commercial stabilizing reagents have been developed aiming to keep RNA stable for long period. However, for studies involving large number of experimental samples, the high costs related to these specific reagents may constitute a barrier.

METHODS AND RESULTS

In this context the present study was designed aiming to evaluate the stability of mRNA in whole bovine blood collected in EDTA tubes during storage at common fridge (4 °C). Whole blood samples were collected from six Holstein calves and submitted to RNA extraction in each different interval: immediately after blood sampling (< 2 h), at 1-day post-sampling (dps), 2 dps, 3 dps, 7 dps and 14dps intervals. RNA integrity and purity were evaluated, and RT-qPCR assays were run using seven different genes (B2M, ACTB, PPIA, GAPDH, YWHAZ, CD4 and IFN-γ) aiming to evaluate the presence of altered gene transcription during storage. All extracted RNA samples presented high purity, while optimal integrity and unaltered gene expression were observed in whole experimental group up to 3 days of storage.

CONCLUSION

Bovine blood RNA remained stable in K3EDTA tubes for 3 days stored at common fridge and can be successfully and accurately used for gene expression studies.

IMMUNOMEDIATOR GENE TRANSCRIPTION PROFILING IN BELUGA WHALE (DELPHINAPTERUS LEUCAS) CLINICAL CASES

There is an unmet need for specific diagnostics of immune perturbations and inflammation in beluga whale (Delphinapterus leucas) clinical care. Quantitative real-time polymerase chain reaction (qPCR) has been used to measure immunomediator gene transcription in beluga whales. The study hypothesis was that a qPCR-based immunomediator assay would supplement routine clinical data with specific and sensitive information on immune status. Two beluga whale clinical cases provided an opportunity to test this hypothesis: a whale with a skin laceration and a whale with gastrointestinal inflammation. Mitogen-stimulated immunomediator gene transcription (MSIGT) was compared between the cases and healthy contact whales. In both case studies, mitogens increased transcription of IL1B, PTGS2 (Cox-2), TNF, HIF1A, and IL2 but decreased IL10 transcription in peripheral blood mononuclear cells (PBMC) from the abnormal whale over the control. Correlations were identified between most immunomediators tested and one or more standard blood clinical values. Considering all 15 immunomediators tested, the whale with gastrointestinal inflammation had a more unique MSIGT signature than the whale with a laceration. These results support further elucidation of beluga whale PBMC cytokine profiles for use as immune biomarkers.

Preanalytical robustness of blood collection tubes with RNA stabilizers

Background

Efficient blood stabilization is essential to obtaining reliable and comparable RNA analysis data in preclinical operations. PAXgene (Qiagen, Becton Dickinson) and Tempus (Applied Biosystems, Life Technologies) blood collection tubes with RNA stabilizers both avoid preanalytical degradation of mRNA by endogenous nucleases and modifications in specific mRNA concentrations by unintentional up- or down-regulation of gene expression.

Methods

Sixteen different preanalytical conditions were tested in PAXgene and Tempus blood samples from seven donors: different mixing after collection, different fill volumes and different 24-h transport temperature conditions after collection. RNA was extracted by column-based methods. The quality of the extracted RNA was assessed by spectrophotometric quantification, A260/A280 purity ratio, RNA Integrity Number (Agilent Bioanalyzer), miRNA quantative real time polymerase chain reaction (qRT-PCR) on two target miRNAs (RNU-24 and miR-16), mRNA quality index by qRT-PCR on the 3′ and 5′ region of the GAPDH gene, and the PBMC preanalytical score, based on the relative expression levels of the IL8 and EDEM3 coding genes.

Results

When PAXgene RNA and Tempus blood collection tubes were used following the manufacturers’ instructions, there was no statistically or technically significant difference in the output RNA quality attributes. However, the integrity of the RNA extracted from Tempus collection tubes was more sensitive to fill volumes and effective inversion, than to storage temperature, while the integrity of RNA extracted from PAXgene collection tubes was more sensitive to effective inversion and storage temperature than to fill volumes.

Conclusions

Blood collection tubes with different RNA stabilizers present different robustness to common preanalytical variations.

Stabilization of cellular RNA in blood during storage at room temperature: a comparison of cell-free RNA BCT(®) with K3EDTA tubes

Background

Messenger RNA (mRNA) expression levels in blood cells are important in disease diagnosis, prognosis and biomarker discovery research. Accurate measurements of intracellular mRNA levels in blood cells depend upon several pre-analytical factors, including delays in RNA extraction from blood after phlebotomy. Dramatic changes in mRNA expression levels caused by delays in blood sample processing may render such samples unsuitable for gene expression analysis.

Objectives

This study was conducted to evaluate a blood collection tube, cell-free RNA-BCT^® (RNA-BCT), for its ability to stabilize mRNA expression level in blood cells post-phlebotomy using indicator mRNAs in reverse transcription quantitative real-time PCR (RT-qPCR) assays.

Methods

Blood samples from presumed healthy donors were drawn into both RNA-BCT and K₃EDTA tubes and maintained at room temperature (18–22 °C). The samples were processed to obtain white blood cells (WBCs) at days 0, 1, 2 and 3. Total cellular RNA was extracted from WBCs and mRNA concentrations were quantified by RT-qPCR for glyceraldehyde-3-phosphate dehydrogenase (GAPDH), c-fos, and p53 transcripts.

Results

While blood cells isolated from K₃EDTA tubes showed significant changes in cellular mRNA concentrations for GAPDH, c-fos, and p53, these mRNAs concentrations were stable in blood drawn into RNA-BCT.

Conclusion

The reagent in the RNA-BCT device stabilizes cellular mRNA concentrations for GAPDH, c-fos and p53 for at least three days at room temperature.

Increasing gene discovery and coverage using RNA-seq of globin RNA reduced porcine blood samples

Background

Transcriptome analysis of porcine whole blood has several applications, which include deciphering genetic mechanisms for host responses to viral infection and vaccination. The abundance of alpha- and beta-globin transcripts in blood, however, impedes the ability to cost-effectively detect transcripts of low abundance. Although protocols exist for reduction of globin transcripts from human and mouse/rat blood, preliminary work demonstrated these are not useful for porcine blood Globin Reduction (GR). Our objectives were to develop a porcine specific GR protocol and to evaluate the GR effects on gene discovery and sequence read coverage in RNA-sequencing (RNA-seq) experiments.

Results

A GR protocol for porcine blood samples was developed using RNase H with antisense oligonucleotides specifically targeting porcine hemoglobin alpha (HBA) and beta (HBB) mRNAs. Whole blood samples (n = 12) collected in Tempus tubes were used for evaluating the efficacy and effects of GR on RNA-seq. The HBA and HBB mRNA transcripts comprised an average of 46.1% of the mapped reads in pre-GR samples, but those reads reduced to an average of 8.9% in post-GR samples. Differential gene expression analysis showed that the expression level of 11,046 genes were increased, whereas 34 genes, excluding HBA and HBB, showed decreased expression after GR (FDR <0.05). An additional 815 genes were detected only in post-GR samples.

Conclusions

Our porcine specific GR primers and protocol minimize the number of reads of globin transcripts in whole blood samples and provides increased coverage as well as accuracy and reproducibility of transcriptome analysis. Increased detection of low abundance mRNAs will ensure that studies relying on transcriptome analyses do not miss information that may be vital to the success of the study.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-954) contains supplementary material, which is available to authorized users.

Reduction of systematic bias in transcriptome data from human peripheral blood mononuclear cells for transportation and biobanking

Transportation of samples is essential for large-scale biobank projects. However, RNA degradation during pre-analytical operations prior to transportation can cause systematic bias in transcriptome data, which may prevent subsequent biomarker identification. Therefore, to collect high-quality biobank samples for expression analysis, specimens must be transported under stable conditions. In this study, we examined the effectiveness of RNA-stabilizing reagents to prevent RNA degradation during pre-analytical operations with an emphasis on RNA from peripheral blood mononuclear cells (PBMCs) to establish a protocol for reducing systematic bias. To this end, we obtained PBMCs from 11 healthy volunteers and analyzed the purity, yield, and integrity of extracted RNA after performing pre-analytical operations for freezing PBMCs at −80°C. We randomly chose 7 samples from 11 samples individually, and systematic bias in expression levels was examined by real-time quantitative reverse transcription polymerase chain reaction (qRT-PCR), RNA sequencing (RNA-Seq) experiments and data analysis. Our data demonstrated that omission of stabilizing reagents significantly lowered RNA integrity, suggesting substantial degradation of RNA molecules due to pre-analytical freezing. qRT-PCR experiments for 19 selected transcripts revealed systematic bias in the expression levels of five transcripts. RNA-Seq for 25,223 transcripts also suggested that about 40% of transcripts were systematically biased. These results indicated that appropriate reduction in systematic bias is essential in protocols for collection of RNA from PBMCs for large-scale biobank projects. Among the seven commercially available stabilizing reagents examined in this study, qRT-PCR and RNA-Seq experiments consistently suggested that RNALock, RNA/DNA Stabilization Reagent for Blood and Bone Marrow, and 1-Thioglycerol/Homogenization solution could reduce systematic bias. On the basis of the results of this study, we established a protocol to reduce systematic bias in the expression levels of RNA transcripts isolated from PBMCs. We believe that these data provide a novel methodology for collection of high-quality RNA from PBMCs for biobank researchers.

Comparison of blood RNA extraction methods used for gene expression profiling in amyotrophic lateral sclerosis

Amyotrophic lateral sclerosis (ALS) is a neurodegenerative disease that causes death within a mean of 2–3 years from symptom onset. There is no diagnostic test and the delay from symptom onset to diagnosis averages 12 months. The identification of prognostic and diagnostic biomarkers in ALS would facilitate earlier diagnosis and faster monitoring of treatments. Gene expression profiling (GEP) can help to identify these markers as well as therapeutic targets in neurological diseases. One source of genetic material for GEP in ALS is peripheral blood, which is routinely accessed from patients. However, a high proportion of globin mRNA in blood can mask important genetic information. A number of methods allow safe collection, storage and transport of blood as well as RNA stabilisation, including the PAXGENE and TEMPUS systems for the collection of whole blood and LEUKOLOCK which enriches for the leukocyte population. Here we compared these three systems and assess their suitability for GEP in ALS. We collected blood from 8 sporadic ALS patients and 7 controls. PAXGENE and TEMPUS RNA extracted samples additionally underwent globin depletion using GlobinClear. RNA was amplified and hybridised onto Affymetrix U133 Plus 2.0 arrays. Lists of genes differentially regulated in ALS patients and controls were created for each method using the R package PUMA, and RT-PCR validation was carried out on selected genes. TEMPUS/GlobinClear, and LEUKOLOCK produced high quality RNA with sufficient yield, and consistent array expression profiles. PAXGENE/GlobinClear yield and quality were lower. Globin depletion for PAXGENE and TEMPUS uncovered the presence of over 60% more transcripts than when samples were not depleted. TEMPUS/GlobinClear and LEUKOLOCK gene lists respectively contained 3619 and 3047 genes differentially expressed between patients and controls. Real-time PCR validation revealed similar reliability between these two methods and gene ontology analyses revealed similar pathways differentially regulated in disease compared to controls.

Isolation of RNA from equine peripheral blood cells: comparison of methods

Gene expression studies in equine research involve the use of whole blood samples as a vital source of RNA. To determine the optimal method for RNA isolation from equine whole blood, we compared three RNA isolation strategies using different commercially available kits to evaluate the yield and quality of equine RNA. All 3 methods produced RNA with high quality. Though it did not produce the highest yield, combining the quality, yield and the need for the downstream application in our project, LeukoLOCK™ total RNA isolation system was the best RNA extraction method.

Genome-wide comparison of two RNA-stabilizing reagents for transcriptional profiling of peripheral blood

Peripheral whole blood is relatively easily obtained for monitoring gene expression for biomarker discovery using transcriptomic platforms such as genome-wide microarrays. However, whole blood provides challenges caused by sensitivity for ex vivo incubation and overrepresentation of globin mRNAs. We compared the performance of 2 commercial whole blood preservation methods, TEMPUS (Applied Biosystems, Foster City, CA) and PAXgene (PreAnalytiX, Qiagen BD, Valencia, CA), using 2 RNA amplification protocols and high-density microarrays. Performance of commercial globin mRNA reduction protocol also was studied. Human peripheral blood samples collected with TEMPUS and PAXgene Blood RNA tubes were amplified with the RiboAmp OA 1 Round RNA Amplification Kit (Arcturus; Applied Biosystems) and the Affymetrix (Santa Clara, CA) small sample protocol. Affymetrix globin reduction protocol was applied for total RNA samples. Samples amplified with RiboAmp were hybridized on Illumina Sentrix HumanRef-8 Expression BeadChips (Illumina Inc, San Diego, CA) and subjected to statistical analyses. RiboAmp mRNA amplification did not notably amplify globin mRNA that is overrepresented in RNA isolated by both TEMPUS and PAXgene preservation. Enzymatic depletion of globin transcript reduced the quality of total RNA and is thus not recommendable. Microarray analysis showed acceptable correlation within and between the RNA preservation methods, but altogether 443 transcripts were differentially expressed between RNA samples preserved in TEMPUS and PAXgene tubes. We demonstrated that the 2 tested blood RNA-preservation methods combined with RiboAmp mRNA amplification may be used for microarray experiments without the need for a prior globin RNA reduction. However, because genes involved in immune cell functions and gene regulatory pathways were differentially expressed as a result of the technical bias between the preservation methods, they should not be used in the same analytic setting.

Whole blood genome-wide gene expression profile in males after prolonged wakefulness and sleep recovery

Although the specific functions of sleep have not been completely elucidated, the literature has suggested that sleep is essential for proper homeostasis. Sleep loss is associated with changes in behavioral, neurochemical, cellular, and metabolic function as well as impaired immune response. Using high-resolution microarrays we evaluated the gene expression profiles of healthy male volunteers who underwent 60 h of prolonged wakefulness (PW) followed by 12 h of sleep recovery (SR). Peripheral whole blood was collected at 8 am in the morning before the initiation of PW (Baseline), after the second night of PW, and one night after SR. We identified over 500 genes that were differentially expressed. Notably, these genes were related to DNA damage and repair and stress response, as well as diverse immune system responses, such as natural killer pathways including killer cell lectin-like receptors family, as well as granzymes and T-cell receptors, which play important roles in host defense. These results support the idea that sleep loss can lead to alterations in molecular processes that result in perturbation of cellular immunity, induction of inflammatory responses, and homeostatic imbalance. Moreover, expression of multiple genes was downregulated following PW and upregulated after SR compared with PW, suggesting an attempt of the body to re-establish internal homeostasis. In silico validation of alterations in the expression of CETN3, DNAJC, and CEACAM genes confirmed previous findings related to the molecular effects of sleep deprivation. Thus, the present findings confirm that the effects of sleep loss are not restricted to the brain and can occur intensely in peripheral tissues.

Efficient recovery of whole blood RNA--a comparison of commercial RNA extraction protocols for high-throughput applications in wildlife species

Background

Since the emergence of next generation sequencing platforms, unprecedented opportunities have arisen in the study of natural vertebrate populations. In particular, insights into the genetic and epigenetic mechanisms of adaptation can be revealed through study of the expression profiles of genes. However, as a pre-requisite to expression profiling, care must be taken in RNA preparation as factors like DNA contamination, RNA integrity or transcript abundance can affect downstream applications. Here, we evaluated five commonly used RNA extraction methods using whole blood sampled under varying conditions from 20 wild carnivores.

Results

Despite the use of minute starting volumes, all methods produced quantifiable RNA extracts (1.4 – 18.4 μg) with varying integrity (RIN 4.6 - 7.7), the latter being significantly affected by the storage and extraction method used. We observed a significant overall effect of the extraction method on DNA contamination. One particular extraction method, the LeukoLOCK™ filter system, yielded high RNA integrity along with low DNA contamination and efficient depletion of hemoglobin transcripts highly abundant in whole blood. In a proof of concept sequencing experiment, we found globin RNA transcripts to occupy up to ¼ of all sequencing reads if libraries were not depleted of hemoglobin prior to sequencing.

Conclusion

By carefully choosing the appropriate RNA extraction method, whole blood can become a valuable source for high-throughput applications like expression arrays or transcriptome sequencing from natural populations. Additionally, candidate genes showing signs of selection could subsequently be genotyped in large population samples using whole blood as a source for RNA without harming individuals from rare or endangered species.

Increased sensitivity of next generation sequencing-based expression profiling after globin reduction in human blood RNA

Background

Transcriptome analysis is of great interest in clinical research, where significant differences between individuals can be translated into biomarkers of disease. Although next generation sequencing provides robust, comparable and highly informative expression profiling data, with several million of tags per blood sample, reticulocyte globin transcripts can constitute up to 76% of total mRNA compromising the detection of low abundant transcripts. We have removed globin transcripts from 6 human whole blood RNA samples with a human globin reduction kit and compared them with the same non-reduced samples using deep Serial Analysis of Gene Expression.

Results

Globin tags comprised 52-76% of total tags in our samples. Out of 21,633 genes only 87 genes were detected at significantly lower levels in the globin reduced samples. In contrast, 11,338 genes were detected at significantly higher levels in the globin reduced samples. Removing globin transcripts allowed us to also identify 2112 genes that could not be detected in the non-globin reduced samples, with roles in cell surface receptor signal transduction, G-protein coupled receptor protein signalling pathways and neurological processes.

Conclusions

The reduction of globin transcripts in whole blood samples constitutes a reproducible and reliable method that can enrich data obtained from next generation sequencing-based expression profiling.

Optimal techniques for mRNA extraction from neonatal salivary supernatant

BACKGROUND

Gene expression profiling of the salivary supernatant is emerging as a new and important source of real-time, systemic, biological information. However, existing technologies prevent RNA extraction of small quantities found in neonatal salivary supernatant.

OBJECTIVE

The aim of this study was to develop techniques to enhance extraction of cell-free RNA from neonatal salivary supernatant.

METHODS

Two saliva samples (10-100 μl) were serially collected from newborns (36-41 weeks' gestation) (n = 13) and stabilized. Total RNA was extracted from salivary supernatant with the use of two modified extraction techniques: Qiagen RNAprotect® Saliva Mini Kit (method 1) and the QIAamp Viral RNA Mini Kit (method 2). Quantitative RT-PCR amplification for GAPDH was performed on extracted salivary samples. Statistical analyses were performed on mean threshold cycle (Ct) levels to compare RNA yield from each protocol. Paired microarray analyses were made between neonatal whole saliva and supernatant (n = 3) to discern gene expression differences between these biolayers.

RESULTS

mRNA was successfully extracted and amplified from all salivary supernatant samples. Extraction with method 2 yielded more RNA than with method 1 (p = 0.008). There was a 7.5% discordance between paired gene expression analyses for whole saliva and supernatant. Genes that were statistically significantly upregulated in supernatant highlighted 16 distinct biological functions not seen in whole saliva. Conversely, only two biological functions were unique to whole saliva.

CONCLUSION

Neonatal cell-free salivary supernatant mRNA may be readily extracted and utilized on downstream applications. These technical enhancements allow for further exploration of the diagnostic potential of the neonatal salivary transcriptome.

Chronic myelogenous leukemia: monitoring response to therapy

Molecular monitoring is a key component of the management of patients with chronic myeloid leukemia. The current recommendation is that molecular monitoring be performed in place of cytogenetic assessment when a major molecular response (MMR) is achieved. With the more potent kinase inhibitors nilotinib and dasatinib now approved as front-line therapy, more patients will achieve an MMR and will benefit from molecular monitoring. There is a strong correlation between certain BCR-ABL1 levels and the cytogenetic response, which means that molecular monitoring may act as a surrogate for cytogenetic response, but only if the BCR-ABL1 values are converted to the international reporting scale. Furthermore, improvements in the limit of BCR-ABL1 detection and reduction of intra-assay variability are ongoing issues of importance for molecular monitoring. Standardization of molecular methods to accurately assess the patient response also remains a challenge, despite the recent certification of international scale reference materials.

Gene expression profiling of human whole blood samples with the Illumina WG-DASL assay

Background

Microarray-based gene expression analysis of peripheral whole blood is a common strategy in the development of clinically relevant biomarker panels for a variety of human diseases. However, the results of such an analysis are often plagued by decreased sensitivity and reliability due to the effects of relatively high levels of globin mRNA in whole blood. Globin reduction assays have been shown to overcome such effects, but they require large amounts of total RNA and may induce distinct gene expression profiles. The Illumina whole genome DASL assay can detect gene expression levels using partially degraded RNA samples and has the potential to detect rare transcripts present in highly heterogeneous whole blood samples without the need for globin reduction. We assessed the utility of the whole genome DASL assay in an analysis of peripheral whole blood gene expression profiles.

Results

We find that gene expression detection is significantly increased with the use of whole genome DASL compared to the standard IVT-based direct hybridization. Additionally, globin-probe negative whole genome DASL did not exhibit significant improvements over globin-probe positive whole genome DASL. Globin reduction further increases the detection sensitivity and reliability of both whole genome DASL and IVT-based direct hybridization with little effect on raw intensity correlations. Raw intensity correlations between total RNA and globin reduced RNA were 0.955 for IVT-based direct hybridization and 0.979 for whole genome DASL.

Conclusions

Overall, the detection sensitivity of the whole genome DASL assay is higher than the IVT-based direct hybridization assay, with or without globin reduction, and should be considered in conjunction with globin reduction methods for future blood-based gene expression studies.

Balanced caloric macronutrient composition downregulates immunological gene expression in human blood cells-adipose tissue diverges

Cardiovascular disease, obesity, and type 2 diabetes are conditions characterized by low-grade systemic inflammation, strongly influenced by lifestyle, but the mechanisms that link these characteristics are poorly understood. Our first objective was to investigate if a normocaloric diet with a calorically balanced macronutrient composition influenced immunological gene expression. Findings regarding the suitability of blood as biological material in nutrigenomics and gene expression profiling have been inconclusive. Our second objective was to compare blood and adipose tissue sample quality in terms of adequacy for DNA-microarray analyses, and to determine tissue-specific gene expression patterns. Blood and adipose tissue samples were collected for gene expression profiling from three obese men before, during, and after a 28-day normocaloric diet intervention where each meal contained an approximately equal caloric load of macronutrients. Time series analyses of blood gene expression revealed a cluster of downregulated genes involved in immunological processes. Blood RNA quality and yield were satisfactory, and DNA-microarray analysis reproducibility was similar in blood and adipose tissue. Gene expression correlation between blood and adipose tissue varied according to gene function, and was especially low for genes involved in immunological and metabolic processes. This suggests that diet composition is of importance in inflammatory processes in blood cells. The findings also suggest that a systems biology approach, in which tissues are studied in parallel, should be employed to fully understand the impact of dietary challenges on the human body.

Transcriptome profiling of whole blood cells identifies PLEK2 and C1QB in human melanoma

Background

Developing analytical methodologies to identify biomarkers in easily accessible body fluids is highly valuable for the early diagnosis and management of cancer patients. Peripheral whole blood is a “nucleic acid-rich” and “inflammatory cell-rich” information reservoir and represents systemic processes altered by the presence of cancer cells.

Methodology/Principal Findings

We conducted transcriptome profiling of whole blood cells from melanoma patients. To overcome challenges associated with blood-based transcriptome analysis, we used a PAXgene™ tube and NuGEN Ovation™ globin reduction system. The combined use of these systems in microarray resulted in the identification of 78 unique genes differentially expressed in the blood of melanoma patients. Of these, 68 genes were further analyzed by quantitative reverse transcriptase PCR using blood samples from 45 newly diagnosed melanoma patients (stage I to IV) and 50 healthy control individuals. Thirty-nine genes were verified to be differentially expressed in blood samples from melanoma patients. A stepwise logit analysis selected eighteen 2-gene signatures that distinguish melanoma from healthy controls. Of these, a 2-gene signature consisting of PLEK2 and C1QB led to the best result that correctly classified 93.3% melanoma patients and 90% healthy controls. Both genes were upregulated in blood samples of melanoma patients from all stages. Further analysis using blood fractionation showed that CD45⁻ and CD45⁺ populations were responsible for the altered expression levels of PLEK2 and C1QB, respectively.

Conclusions/Significance

The current study provides the first analysis of whole blood-based transcriptome biomarkers for malignant melanoma. The expression of PLEK2, the strongest gene to classify melanoma patients, in CD45⁻ subsets illustrates the importance of analyzing whole blood cells for biomarker studies. The study suggests that transcriptome profiling of blood cells could be used for both early detection of melanoma and monitoring of patients for residual disease.

Global transcriptomic profiling using small volumes of whole blood: a cost-effective method for translational genomic biomarker identification in small animals

Blood is an ideal tissue for the identification of novel genomic biomarkers for toxicity or efficacy. However, using blood for transcriptomic profiling presents significant technical challenges due to the transcriptomic changes induced by ex vivo handling and the interference of highly abundant globin mRNA. Most whole blood RNA stabilization and isolation methods also require significant volumes of blood, limiting their effective use in small animal species, such as rodents. To overcome these challenges, a QIAzol-based RNA stabilization and isolation method (QSI) was developed to isolate sufficient amounts of high quality total RNA from 25 to 500 μL of rat whole blood. The method was compared to the standard PAXgene Blood RNA System using blood collected from rats exposed to saline or lipopolysaccharide (LPS). The QSI method yielded an average of 54 ng total RNA per μL of rat whole blood with an average RNA Integrity Number (RIN) of 9, a performance comparable with the standard PAXgene method. Total RNA samples were further processed using the NuGEN Ovation Whole Blood Solution system and cDNA was hybridized to Affymetrix Rat Genome 230 2.0 Arrays. The microarray QC parameters using RNA isolated with the QSI method were within the acceptable range for microarray analysis. The transcriptomic profiles were highly correlated with those using RNA isolated with the PAXgene method and were consistent with expected LPS-induced inflammatory responses. The present study demonstrated that the QSI method coupled with NuGEN Ovation Whole Blood Solution system is cost-effective and particularly suitable for transcriptomic profiling of minimal volumes of whole blood, typical of those obtained with small animal species.

Gene expression profiling in whole blood identifies distinct biological pathways associated with obesity

Background

Obesity is reaching epidemic proportions and represents a significant risk factor for cardiovascular disease, diabetes, and cancer.

Methods

To explore the relationship between increased body mass and gene expression in blood, we conducted whole-genome expression profiling of whole blood from seventeen obese and seventeen well matched lean subjects. Gene expression data was analyzed at the individual gene and pathway level and a preliminary assessment of the predictive value of blood gene expression profiles in obesity was carried out.

Results

Principal components analysis of whole-blood gene expression data from obese and lean subjects led to efficient separation of the two cohorts. Pathway analysis by gene-set enrichment demonstrated increased transcript levels for genes belonging to the "ribosome", "apoptosis" and "oxidative phosphorylation" pathways in the obese cohort, consistent with an altered metabolic state including increased protein synthesis, enhanced cell death from proinflammatory or lipotoxic stimuli, and increased energy demands. A subset of pathway-specific genes acted as efficient predictors of obese or lean class membership when used in Naive Bayes or logistic regression based classifiers.

Conclusion

This study provides a comprehensive characterization of the whole blood transcriptome in obesity and demonstrates that the investigation of gene expression profiles from whole blood can inform and illustrate the biological processes related to regulation of body mass. Additionally, the ability of pathway-related gene expression to predict class membership suggests the feasibility of a similar approach for identifying clinically useful blood-based predictors of weight loss success following dietary or surgical interventions.

cDNA targets improve whole blood gene expression profiling and enhance detection of pharmocodynamic biomarkers: a quantitative platform analysis

Background

Genome-wide gene expression profiling of whole blood is an attractive method for discovery of biomarkers due to its non-invasiveness, simple clinical site processing and rich biological content. Except for a few successes, this technology has not yet matured enough to reach its full potential of identifying biomarkers useful for clinical prognostic and diagnostic applications or in monitoring patient response to therapeutic intervention. A variety of technical problems have hampered efforts to utilize this technology for identification of biomarkers. One significant hurdle has been the high and variable concentrations of globin transcripts in whole blood total RNA potentially resulting in non-specific probe binding and high background. In this study, we investigated and quantified the power of three whole blood profiling approaches to detect meaningful biological expression patterns.

Methods

To compare and quantify the impact of different mitigation technologies, we used a globin transcript spike-in strategy to synthetically generate a globin-induced signature and then mitigate it with the three different technologies. Biological differences, in globin transcript spiked samples, were modeled by supplementing with either 1% of liver or 1% brain total RNA. In order to demonstrate the biological utility of a robust globin artifact mitigation strategy in biomarker discovery, we treated whole blood ex vivo with suberoylanilide hydroxamic acid (SAHA) and compared the overlap between the obtained signatures and signatures of a known biomarker derived from SAHA-treated cell lines and PBMCs of SAHA-treated patients.

Results

We found cDNA hybridization targets detect at least 20 times more specific differentially expressed signatures (2597) between 1% liver and 1% brain in globin-supplemented samples than the PNA (117) or no treatment (97) method at FDR = 10% and p-value < 3x10-3. In addition, we found that the ex vivo derived gene expression profile was highly concordant with that of the previously identified SAHA pharmacodynamic biomarkers.

Conclusions

We conclude that an amplification method for gene expression profiling employing cDNA targets effectively mitigates the negative impact on data of abundant globin transcripts and greatly improves the ability to identify relevant gene expression based pharmacodynamic biomarkers from whole blood.

Gene expression patterns of dengue virus-infected children from nicaragua reveal a distinct signature of increased metabolism

Background

Infection with dengue viruses (DENV) leads to a spectrum of disease outcomes. The pathophysiology of severe versus non-severe manifestations of DENV infection may be driven by host responses, which could be reflected in the transcriptional profiles of peripheral blood immune cells.

Methodology/Principal Findings

We conducted genome-wide microarray analysis of whole blood RNA from 34 DENV-infected children in Nicaragua collected on days 3–6 of illness, with different disease manifestations. Gene expression analysis identified genes that are differentially regulated between clinical subgroups. The most striking transcriptional differences were observed between dengue patients with and without shock, especially in the expression of mitochondrial ribosomal proteins associated with protein biosynthesis. In the dengue hemorrhagic fever patients, one subset of differentially expressed genes encode neutrophil-derived anti-microbial peptides associated with innate immunity. By performing a meta-analysis of our dataset in conjunction with previously published datasets, we confirmed that DENV infection in vivo is associated with large changes to protein and nucleic acid metabolism. Additionally, whereas in vitro infection leads to an increased interferon signature, this was not consistently observed from in vivo patient samples, suggesting that the interferon response in vivo is relatively transient and was no longer observed by days 3–6 of illness.

Conclusions/Significance

These data highlight important differences between different manifestations of severity during DENV infection as well as identify some commonalities. Compilation of larger datasets in the future across multiple studies, as we have initiated in this report, may well lead to better prediction of disease manifestation via a systems biology approach.

Dengue is a widespread viral disease for which over 3 billion people are at risk. There are no drug treatments or vaccines available for this disease. It is also difficult for physicians to predict which patients are at highest risk for the severe manifestations known as dengue hemorrhagic fever (DHF) and dengue shock syndrome (DSS). We used genome-wide transcriptional profiling analysis to study peripheral blood responses to dengue among patients from Nicaragua. We found that patients with severe manifestations involving shock had very different transcriptional profiles from dengue patients with mild and moderate illness. We then compared our results with other microarray experiments on dengue patients available from public databases and confirmed that dengue is often associated with large changes to the metabolic processes within cells. This approach could identify prognostic markers for severe dengue as well as provide a better understanding of the pathophysiology associated with different grades of disease severity.

The effects of globin on microarray-based gene expression analysis of mouse blood

The use of mouse blood as a model for human blood is often considered in the development of clinically relevant, gene expression-based disease biomarkers. However, the ability to derive biologically meaningful insights from microarray-based gene expression patterns in mouse whole blood, as in human whole blood, is hindered by high levels of globin mRNA. In order to characterize the effects of globin reduction on gene expression of peripheral mouse blood, we performed gene set enrichment analysis on genes identified as expressed in blood via microarray-based genome-wide transcriptome analysis. Depletion of globin mRNA enhanced the quality of microarray data as shown by improved gene expression detection and increased sensitivity. Compared to genes expressed in whole blood, genes detected as expressed in blood following globin reduction were enriched for low abundance transcripts implicated in many biological pathways, including development, g-protein signaling, and immune response. Broadly, globin reduction resulted in improved detection of expressed genes that serve as molecular binding proteins and enzymes in cellular metabolism, intracellular transport/localization, transcription, and translation, as well as genes that potentially could act as biomarkers for diseases such as schizophrenia. These significantly enriched pathways overlap considerably with those identified in globin-reduced human blood suggesting that globin-reduced mouse blood gene expression studies may be useful for identifying genes relevant to human disease. Overall, the results of this investigation provide a better understanding of the impact of reducing globin transcripts in mouse blood and highlight the potential of microarray-based, globin-reduced, mouse blood gene expression studies in biomarker development.

Variability of gene expression profiles in human blood and lymphoblastoid cell lines

BACKGROUND

Readily accessible samples such as peripheral blood or cell lines are increasingly being used in large cohorts to characterise gene expression differences between a patient group and healthy controls. However, cell and RNA isolation procedures and the variety of cell types that make up whole blood can affect gene expression measurements. We therefore systematically investigated global gene expression profiles in peripheral blood from six individuals collected during two visits by comparing five of the following cell and RNA isolation methods: whole blood (PAXgene), peripheral blood mononuclear cells (PBMCs), lymphoblastoid cell lines (LCLs), CD19 and CD20 specific B-cell subsets.

RESULTS

Gene expression measurements were clearly discriminated by isolation method although the reproducibility was high for all methods (range rho = 0.90-1.00). The PAXgene samples showed a decrease in the number of expressed genes (P < 1*10(-16)) with higher variability (P < 1*10(-16)) compared to the other methods. Differentially expressed probes between PAXgene and PBMCs were correlated with the number of monocytes, lymphocytes, neutrophils or erythrocytes. The correlations (rho = 0.83; rho = 0.79) of the expression levels of detected probes between LCLs and B-cell subsets were much lower compared to the two B-cell isolation methods (rho = 0.98). Gene ontology analysis of detected genes showed that genes involved in inflammatory responses are enriched in B-cells CD19 and CD20 whereas genes involved in alcohol metabolic process and the cell cycle were enriched in LCLs.

CONCLUSION

Gene expression profiles in blood-based samples are strongly dependent on the predominant constituent cell type(s) and RNA isolation method. It is crucial to understand the differences and variability of gene expression measurements between cell and RNA isolation procedures, and their relevance to disease processes, before application in large clinical studies.

Validation of extraction methods for total RNA and miRNA from bovine blood prior to quantitative gene expression analyses

The benefit and precision of blood diagnosis by quantitative real-time PCR (qPCR) is limited by sampling procedures and RNA extraction methods. We have compared five different RNA extraction protocols from bovine blood regarding RNA and miRNA yield, quality, and most reproducible data in the qRT-PCR with the lowest point of quantification. Convincing results in terms of highest quantity, quality, and best performance for mRNA qPCR were obtained by leukocyte extraction following blood lysis as well as extraction of PAXgene stabilized blood. The best microRNA qPCR results were obtained for samples extracted by the leukocyte extraction method.

A practical platform for blood biomarker study by using global gene expression profiling of peripheral whole blood

Background

Although microarray technology has become the most common method for studying global gene expression, a plethora of technical factors across the experiment contribute to the variable of genome gene expression profiling using peripheral whole blood. A practical platform needs to be established in order to obtain reliable and reproducible data to meet clinical requirements for biomarker study.

Methods and Findings

We applied peripheral whole blood samples with globin reduction and performed genome-wide transcriptome analysis using Illumina BeadChips. Real-time PCR was subsequently used to evaluate the quality of array data and elucidate the mode in which hemoglobin interferes in gene expression profiling.

We demonstrated that, when applied in the context of standard microarray processing procedures, globin reduction results in a consistent and significant increase in the quality of beadarray data. When compared to their pre-globin reduction counterparts, post-globin reduction samples show improved detection statistics, lowered variance and increased sensitivity. More importantly, gender gene separation is remarkably clearer in post-globin reduction samples than in pre-globin reduction samples. Our study suggests that the poor data obtained from pre-globin reduction samples is the result of the high concentration of hemoglobin derived from red blood cells either interfering with target mRNA binding or giving the pseudo binding background signal.

Conclusion

We therefore recommend the combination of performing globin mRNA reduction in peripheral whole blood samples and hybridizing on Illumina BeadChips as the practical approach for biomarker study.

Interindividual and interethnic variation in genomewide gene expression: insights into the biological variation of gene expression and clinical implications

BACKGROUND

Analysis of gene expression in peripheral blood samples is increasingly being applied in biomarker studies of disease diagnosis and prognosis. Although knowledge of interindividual and interethnic variation in gene expression is required to set ethnicity-specific reference intervals and to select reference genes and preferred markers from a list of candidate genes, few studies have attempted to characterize such biological variation on a genomewide scale.

METHODS

The genomewide expression profiles of 11 355 transcripts expressed among 210 multiethnic individuals of the HapMap project were obtained and analyzed; 4 replicates were included for each sample. The total biological CV in gene expression (CV(b)) was partitioned into interindividual (CV(g)), inter-ethnic group (CV(e)), and residual components by random-effects mixed models.

RESULTS

CV(g) was the major component of CV(b), and the differences among transcripts were large (up to 38%). Distinct groups of genes were characterized by CV values and expression levels. Of the genes with lowest biological variation (CV(b) < 1.5%), 35 genes were highly expressed, whereas 32 had intermediate or low expression. Although CV(g) was almost always greater than CV(e), we identified 10 genes in which ethnic variation predominated (range, 8%-18%). On the other hand, 17 annotated genes were highly variable with CV(g) values ranging between 15% and 38%.

CONCLUSIONS

Genomewide analysis of gene expression variation demonstrated biological differences among transcripts. Transcripts with the least biological variation are better candidates for reference genes, whereas those with low interindividual variation may be good disease markers. The presence of interethnic variation suggests that ethnicity-specific reference intervals may be necessary.

Gene expression profiling of whole blood: comparison of target preparation methods for accurate and reproducible microarray analysis

Background

Peripheral blood is an accessible and informative source of transcriptomal information for many human disease and pharmacogenomic studies. While there can be significant advantages to analyzing RNA isolated from whole blood, particularly in clinical studies, the preparation of samples for microarray analysis is complicated by the need to minimize artifacts associated with highly abundant globin RNA transcripts. The impact of globin RNA transcripts on expression profiling data can potentially be reduced by using RNA preparation and labeling methods that remove or block globin RNA during the microarray assay. We compared four different methods for preparing microarray hybridization targets from human whole blood collected in PAXGene tubes. Three of the methods utilized the Affymetrix one-cycle cDNA synthesis/in vitro transcription protocol but varied treatment of input RNA as follows: i. no treatment; ii. treatment with GLOBINclear; or iii. treatment with globin PNA oligos. In the fourth method cDNA targets were prepared with the Ovation amplification and labeling system.

Results

We find that microarray targets generated with labeling methods that reduce globin mRNA levels or minimize the impact of globin transcripts during hybridization detect more transcripts in the microarray assay compared with the standard Affymetrix method. Comparison of microarray results with quantitative PCR analysis of a panel of genes from the NF-kappa B pathway shows good correlation of transcript measurements produced with all four target preparation methods, although method-specific differences in overall correlation were observed. The impact of freezing blood collected in PAXGene tubes on data reproducibility was also examined. Expression profiles show little or no difference when RNA is extracted from either fresh or frozen blood samples.

Conclusion

RNA preparation and labeling methods designed to reduce the impact of globin mRNA transcripts can significantly improve the sensitivity of the DNA microarray expression profiling assay for whole blood samples. While blockage of globin transcripts during first strand cDNA synthesis with globin PNAs resulted in the best overall performance in this study, we conclude that selection of a protocol for expression profiling studies in blood should depend on several factors, including implementation requirements of the method and study design. RNA isolated from either freshly collected or frozen blood samples stored in PAXGene tubes can be used without altering gene expression profiles.

Intraplatform reproducibility and technical precision of gene expression profiling in 4 laboratories investigating 160 leukemia samples: the DACH study

BACKGROUND

Gene expression profiling has the potential to offer consistent, objective diagnostic test results once a standardized protocol has been established. We investigated the robustness, precision, and reproducibility of microarray technology.

METHODS

One hundred sixty individual patient samples representing 11 subtypes of acute and chronic leukemias, myelodysplastic syndromes, and nonleukemia as a control group were centrally collected and diagnosed as part of the daily routine in the Munich Leukemia Laboratory. The custom AmpliChip Leukemia research microarray was used for technical analyses of quadruplicate mononuclear cell lysates in 4 different laboratories in Germany (D), Austria (A), and Switzerland (CH) (the DACH study).

RESULTS

Total-RNA preparations were successfully performed in 637 (99.5%) of 640 cases. Mean differences between pairs of laboratories in the total-RNA yield from the same sample ranged from 0.02 mug to 1.03 mug. Further processing produced 622 successful in vitro transcription reactions (97.6%); the mean differences between laboratories in the cRNA yield from the same sample ranged from 0.40 mug to 6.18 mug. After hybridization to microarrays, a mean of 47.6%, 46.5%, 46.2%, and 46.4% of probe sets were detected as present for the 4 laboratories, with mean signal-intensity scaling factors of 3.1, 3.7, 4.0, and 4.2, respectively. In unsupervised hierarchical cluster and principal component analyses, replicates from the same patient always clustered closely together, with no indications of any association between gene expression profiles due to different operators or laboratories.

CONCLUSIONS

Microarray analysis can be performed with high interlaboratory reproducibility and with comparable quality and high technical precision across laboratories.

Improved microarray gene expression profiling of virus-infected cells after removal of viral RNA

Background

Sensitivity and accuracy are key points when using microarrays to detect alterations in gene expression under different conditions. Critical to the acquisition of reliable results is the preparation of the RNA. In the field of virology, when analyzing the host cell's reaction to infection, the often high representation of viral RNA (vRNA) within total RNA preparations from infected cells is likely to interfere with microarray analysis. Yet, this effect has not been investigated despite the many reports that describe gene expression profiling of virus-infected cells using microarrays.

Results

In this study we used coronaviruses as a model to show that vRNA indeed interferes with microarray analysis, decreasing both sensitivity and accuracy. We also demonstrate that the removal of vRNA from total RNA samples, by means of virus-specific oligonucleotide capturing, significantly reduced the number of false-positive hits and increased the sensitivity of the method as tested on different array platforms.

Conclusion

We therefore recommend the specific removal of vRNA, or of any other abundant 'contaminating' RNAs, from total RNA samples to improve the quality and reliability of microarray analyses.