Assessment of Coagulation Homeostasis in Blunt, Penetrating, and Thermal Trauma: Guidance for a Multicenter Systems Biology Approach.

INTRODUCTION: Provisioning care for traumatically injured patients makes conducting research very proximal to injury difficult. These studies also inherently have regulatory barriers to overcome. Here we outline a protocol for acute-phase enrollment of traumatically injured patients into a prospective observational clinical trial with precise and comprehensive sample acquisition in support of a systems biology approach to a research study. METHODS: Experts in trauma, burn, blood coagulation, computational biology, and integrative systems biology developed a prospective study that would capture the natural history of coagulation pathology after traumatic injury. Blood was sampled at admission and serial time points throughout hospitalization. Concurrently, demographic and outcomes data were recorded and on-site point-of-care testing was implemented. Protocols were harmonized across sites and sampling protocols validated through demonstration of feasibility and sample quality assurance testing. A novel data integration platform was developed to store, visualize, and enable large-scale analysis of empirical and clinical data. Regulatory considerations were also addressed in protocol development. RESULTS: A comprehensive Manual of Operations (MOO) was developed and implemented at 3 clinical sites. After regulatory approval, the MOO was followed to collect 5,348 longitudinal samples from 1,547 patients. All samples were collected, processed, and stored per the MOO. Assay results and clinical data were entered into the novel data management platform for analyses. CONCLUSION: We used an iterative, interdisciplinary process to develop a systematic and robust protocol for comprehensive assessment of coagulation in traumatically injured patients. This MOO can be a template for future studies in the acute setting.

Widespread Down-Regulation of Cardiac Mitochondrial and Sarcomeric Genes in Patients With Sepsis.

OBJECTIVES: The mechanism(s) for septic cardiomyopathy in humans is not known. To address this, we measured messenger RNA alterations in hearts from patients who died from systemic sepsis, in comparison to changed messenger RNA expression in nonfailing and failing human hearts. DESIGN: Identification of genes with altered abundance in septic cardiomyopathy, ischemic heart disease, or dilated cardiomyopathy, in comparison to nonfailing hearts. SETTING: ICUs at Barnes-Jewish Hospital, St. Louis, MO. PATIENTS: Twenty sepsis patients, 11 ischemic heart disease, nine dilated cardiomyopathy, and 11 nonfailing donors. INTERVENTIONS: None other than those performed as part of patient care. MEASUREMENTS AND MAIN RESULTS: Messenger RNA expression levels for 198 mitochondrially localized energy production components, including Krebs cycle and electron transport genes, decreased by 43% ± 5% (mean ± SD). Messenger RNAs for nine genes responsible for sarcomere contraction and excitation-contraction coupling decreased by 43% ± 4% in septic hearts. Surprisingly, the alterations in messenger RNA levels in septic cardiomyopathy were both distinct from and more profound than changes in messenger RNA levels in the hearts of patients with end-stage heart failure. CONCLUSIONS: The expression profile of messenger RNAs in the heart of septic patients reveals striking decreases in expression levels of messenger RNAs that encode proteins involved in cardiac energy production and cardiac contractility and is distinct from that observed in patients with heart failure. Although speculative, the global nature of the decreases in messenger RNA expression for genes involved in cardiac energy production and contractility suggests that these changes may represent a short-term adaptive response of the heart in response to acute change in cardiovascular homeostasis.

Evolving beyond the vicious triad: Differential mediation of traumatic coagulopathy by injury, shock, and resuscitation.

BACKGROUND: A subset of trauma patients with critical injury present with coagulopathy, portending markedly worse outcomes. Clinical practice is evolving to treat the classical risk factors of hypothermia, hemodilution, and acidosis; however, coagulopathy persists even in the absence of these factors. We sought to determine the relative importance of injury- and shock-specific factors compared with resuscitation-associated factors in coagulopathy after trauma. METHODS: Comprehensive demographic data, laboratory data, and outcomes data were prospectively collected from seven trauma centers over 8 years (November 2003 to August 2011) as part of the Inflammation and the Host Response to Injury Large-Scale Collaborative Program. A total of 1,537 critically injured patients with blunt trauma and hemorrhagic shock were analyzed to evaluate predictors of admission coagulopathy (international normalized ratio [INR] ≥ 1.3), multiorgan failure, and mortality. RESULTS: Of 1,537 patients, 578 (37.6%) had admission INR of 1.3 or greater. Coagulopathic patients had more severe injury, more severe base deficit and lactate levels, as well as lower admission temperature, lower pH, and higher prehospital crystalloid volume (all p < 0.001). Coagulopathic patients required more blood products and mechanical ventilation and had higher rates of nosocomial infection, multiorgan failure, and mortality (all p < 0.02). Injury severity, temperature, and acidosis (all p < 0.02) independently predicted coagulopathy in multivariate analysis, with a significant interaction between lactate and prehospital crystalloid. In Cox regression models, however, coagulopathy itself remained an independent predictor of both multiorgan failure and mortality (p < 0.02) even when adjusted for injury severity, shock, and elements of the vicious triad. CONCLUSION: Most patients with coagulopathy after trauma have mixed risk factors; however, coagulopathy has deleterious effects independent of injury severity, shock, and the vicious triad. Better understanding of the biochemical mechanisms of acute traumatic coagulopathy may facilitate biochemically targeted resuscitation strategies and improve outcomes. LEVEL OF EVIDENCE: Prognostic and epidemiologic study, level II.

Reactivation of multiple viruses in patients with sepsis.

A current controversy is whether patients with sepsis progress to an immunosuppressed state. We hypothesized that reactivation of latent viruses occurred with prolonged sepsis thereby providing evidence of clinically-relevant immunosuppression and potentially providing a means to serially-monitor patients' immune status. Secondly, if viral loads are markedly elevated, they may contribute to morbidity and mortality. This study determined if reactivation of herpesviruses, polyomaviruses, and the anellovirus TTV occurred in sepsis and correlated with severity. Serial whole blood and plasma samples from 560 critically-ill septic, 161 critically-ill non-septic, and 164 healthy age-matched patients were analyzed by quantitative-polymerase-chain-reaction for cytomegalovirus (CMV), Epstein-Barr (EBV), herpes-simplex (HSV), human herpes virus-6 (HHV-6), and TTV. Polyomaviruses BK and JC were quantitated in urine. Detectable virus was analyzed with respect to secondary fungal and opportunistic bacterial infections, ICU duration, severity of illness, and survival. Patients with protracted sepsis had markedly increased frequency of detectable virus. Cumulative viral DNA detection rates in blood were: CMV (24.2%), EBV (53.2%), HSV (14.1%), HHV-6 (10.4%), and TTV (77.5%). 42.7% of septic patients had presence of two or more viruses. The 50% detection rate for herpesviruses was 5-8 days after sepsis onset. A small subgroup of septic patients had markedly elevated viral loads (>104-106 DNA copies/ml blood) for CMV, EBV, and HSV. Excluding TTV, DNAemia was uncommon in critically-ill non-septic patients and in age-matched healthy controls. Compared to septic patients without DNAemia, septic patients with viremia had increased fungal and opportunistic bacterial infections. Patients with detectable CMV in plasma had higher 90-day mortality compared to CMV-negative patients; p<0.05. Reactivation of latent viruses is common with prolonged sepsis, with frequencies similar to those occurring in transplant patients on immunosuppressive therapy and consistent with development of an immunosuppressive state. Whether reactivated latent viruses contribute to morbidity and mortality in sepsis remains unknown.

Smooth muscle cells from abdominal aortic aneurysms are unique and can independently and synergistically degrade insoluble elastin.

BACKGROUND: The purpose of this study was to further elucidate the role of the vascular smooth muscle cells (SMCs) in abdominal aortic aneurysm (AAA) disease. We hypothesized that that AAA SMCs are unique and actively participate in the process of degrading the aortic matrix. METHODS: Whole-genome expression profiles of SMCs from AAAs, nondilated abdominal aorta (NAA), and carotid endarterectomy (CEA) were compared. We quantified elastolytic activity by culturing SMCs in [(3)H]elastin-coated plates and measuring solubilized tritium in the media after 7 days. Matrix metalloproteinase (MMP)-2 and MMP-9 production was assessed using real-time polymerase chain reaction, zymography, and Western blotting. RESULTS: Each SMC type exhibited a unique gene expression pattern. AAA SMCs had greater elastolytic activity than NAA-SMCs (+68%; P < .001) and CEA-SMCs (+45%; P < .001). Zymography showed an increase of active MMP-2 (62 kD) in media from AAA SMCs. AAA SMCs demonstrated twofold greater expression of MMP-2 messenger (m)RNA (P < .05) and 7.3-fold greater MMP-9 expression (P < .01) than NAA-SMCs. Culture with U937 monocytes caused a synergistic increase of elastolysis by AAA SMCs (41%; P < .001) but not NAA-SMCs or CEA-SMCs (P = .99). Coculture with U937 caused a large increase in MMP-9 mRNA in AAA-SMCs and NAA-SMCs (P < .001). MMP-2 mRNA expression was not affected. Western blots of culture media showed a fourfold increase of MMP-9 (92 kD) protein only in AAA-SMCs/U937 but not in NAA-SMCs/U937 (P < .001) and a large increase in active-MMP2 (62 kD), which was less apparent in NAA-SMCs/U937 media (P < .01). CONCLUSIONS: AAA-SMCs have a unique gene expression profile and a proelastolytic phenotype that is augmented by macrophages. This may occur by a failure of post-transcriptional control of MMP-9 synthesis.

Genomic responses in mouse models poorly mimic human inflammatory diseases.

A cornerstone of modern biomedical research is the use of mouse models to explore basic pathophysiological mechanisms, evaluate new therapeutic approaches, and make go or no-go decisions to carry new drug candidates forward into clinical trials. Systematic studies evaluating how well murine models mimic human inflammatory diseases are nonexistent. Here, we show that, although acute inflammatory stresses from different etiologies result in highly similar genomic responses in humans, the responses in corresponding mouse models correlate poorly with the human conditions and also, one another. Among genes changed significantly in humans, the murine orthologs are close to random in matching their human counterparts (e.g., R(2) between 0.0 and 0.1). In addition to improvements in the current animal model systems, our study supports higher priority for translational medical research to focus on the more complex human conditions rather than relying on mouse models to study human inflammatory diseases.

A genomic storm in critically injured humans.

Human survival from injury requires an appropriate inflammatory and immune response. We describe the circulating leukocyte transcriptome after severe trauma and burn injury, as well as in healthy subjects receiving low-dose bacterial endotoxin, and show that these severe stresses produce a global reprioritization affecting >80% of the cellular functions and pathways, a truly unexpected "genomic storm." In severe blunt trauma, the early leukocyte genomic response is consistent with simultaneously increased expression of genes involved in the systemic inflammatory, innate immune, and compensatory antiinflammatory responses, as well as in the suppression of genes involved in adaptive immunity. Furthermore, complications like nosocomial infections and organ failure are not associated with any genomic evidence of a second hit and differ only in the magnitude and duration of this genomic reprioritization. The similarities in gene expression patterns between different injuries reveal an apparently fundamental human response to severe inflammatory stress, with genomic signatures that are surprisingly far more common than different. Based on these transcriptional data, we propose a new paradigm for the human immunological response to severe injury.

Clinical microfluidics for neutrophil genomics and proteomics.

Neutrophils have key roles in modulating the immune response. We present a robust methodology for rapidly isolating neutrophils directly from whole blood with 'on-chip' processing for mRNA and protein isolation for genomics and proteomics. We validate this device with an ex vivo stimulation experiment and by comparison with standard bulk isolation methodologies. Last, we implement this tool as part of a near-patient blood processing system within a multi-center clinical study of the immune response to severe trauma and burn injury. The preliminary results from a small cohort of subjects in our study and healthy controls show a unique time-dependent gene expression pattern clearly demonstrating the ability of this tool to discriminate temporal transcriptional events of neutrophils within a clinical setting.

Validation of the riboleukogram to detect ventilator-associated pneumonia after severe injury.

OBJECTIVE: We hypothesized that circulating leukocyte RNA profiles or "riboleukograms" detect ventilator-associated pneumonia after blunt trauma. SUMMARY BACKGROUND DATA: A pilot microarray study of 11 ventilator-associated pneumonia (VAP) patients suggested that 85 leukocyte genes can be used to diagnose VAP. Validation of this gene set to detect VAP was tested using data from an independent patient cohort. METHODS: A total of 158 intubated blunt trauma patients were enrolled at 5 centers, where 57 (36%) developed VAP. Patient age was 34.2 ± 11.1 years; 65% were male. Circulating leukocyte GeneChip U133 2.0 expression values were measured at time 0.5, 1, 4, 7, 14, 21, and 28 days after injury. DChip normalized leukocyte transcriptional profiles were analyzed using repeated measures logistic regression. A compound covariate model based on leukocyte gene transcriptional profiles in a training subset of patients was tested to determine predictive accuracy for VAP 4 days prior to clinical diagnosis in the test subset. RESULTS: Using gene expression values measured on each study day at an FDR <0.05, 27 (32%) of the 85 genes were associated with the diagnosis of VAP 1 to 4 days before diagnosis. However, the compound covariate model based on these 85-genes did not predict VAP in the test cohort better than chance (P = 0.27). In contrast, a compound covariate model based upon de novo transcriptional analysis of the 158 patients predicted VAP better than chance 4 days before diagnosis with a sensitivity of 57% and a specificity of 69%. CONCLUSION: Our results validate those described in a pilot study, confirming that riboleukograms are associated with the development of VAP days prior to clinical diagnosis. Similarly, a riboleukogram predictive model tested on a larger cohort of 158 patients was better than chance at predicting VAP days prior to clinical diagnosis.

A genomic score prognostic of outcome in trauma patients.

Traumatic injuries frequently lead to infection, organ failure, and death. Health care providers rely on several injury scoring systems to quantify the extent of injury and to help predict clinical outcome. Physiological, anatomical, and clinical laboratory analytic scoring systems (Acute Physiology and Chronic Health Evaluation [APACHE], Injury Severity Score [ISS]) are utilized, with limited success, to predict outcome following injury. The recent development of techniques for measuring the expression level of all of a person's genes simultaneously may make it possible to develop an injury scoring system based on the degree of gene activation. We hypothesized that a peripheral blood leukocyte gene expression score could predict outcome, including multiple organ failure, following severe blunt trauma. To test such a scoring system, we measured gene expression of peripheral blood leukocytes from patients within 12 h of traumatic injury. cRNA derived from whole blood leukocytes obtained within 12 h of injury provided gene expression data for the entire genome that were used to create a composite gene expression score for each patient. Total blood leukocytes were chosen because they are active during inflammation, which is reflective of poor outcome. The gene expression score combines the activation levels of all the genes into a single number which compares the patient's gene expression to the average gene expression in uninjured volunteers. Expression profiles from healthy volunteers were averaged to create a reference gene expression profile which was used to compute a difference from reference (DFR) score for each patient. This score described the overall genomic response of patients within the first 12 h following severe blunt trauma. Regression models were used to compare the association of the DFR, APACHE, and ISS scores with outcome. We hypothesized that patients with a total gene response more different from uninjured volunteers would tend to have poorer outcome than those more similar. Our data show that for measures of poor outcome, such as infections, organ failures, and length of hospital stay, this is correct. DFR scores were associated significantly with adverse outcome, including multiple organ failure, duration of ventilation, length of hospital stay, and infection rate. The association remained significant after adjustment for injury severity as measured by APACHE or ISS. A single score representing changes in gene expression in peripheral blood leukocytes within hours of severe blunt injury is associated with adverse clinical outcomes that develop later in the hospital course. Assessment of genome-wide gene expression provides useful clinical information that is different from that provided by currently utilized anatomic or physiologic scores.

Comparison of longitudinal leukocyte gene expression after burn injury or trauma-hemorrhage in mice.

A primary objective of the large collaborative project entitled "Inflammation and the Host Response to Injury" was to identify leukocyte genes that are differentially expressed after two different types of injury in mouse models and to test the hypothesis that both forms of injury would induce similar changes in gene expression. We report here the genes that are expressed in white blood cells (WBCs) and in splenocytes at 2 h, 1 day, 3 days, and 7 days after burn and sham injury or trauma-hemorrhage (T-H) and sham T-H. Affymetrix Mouse Genome 430 2.0 GeneChips were used to profile gene expression, and the results were analyzed by dCHIP, BRB Array Tools, and Ingenuity Pathway Analysis (IPA) software. We found that the highest number of genes differentially expressed following burn injury were at day 1 for both WBCs (4,989) and for splenocytes (4,715) and at day 1 for WBCs (1,167) and at day 3 for splenocytes (1,117) following T-H. The maximum overlap of genes that were expressed after both forms of injury were at day 1 in WBCs (136 genes) and at day 7 in splenocytes (433 genes). IPA revealed that the cell-to-cell signaling, cell death, immune response, antiapoptosis, and cell cycle control pathways were affected most significantly. In summary, this report provides a database of genes that are modulated in WBCs and splenocytes at sequential time points after burn or T-H in mice and reveals that relatively few leukocyte genes are expressed in common after these two forms of injury.

Cell-specific expression and pathway analyses reveal alterations in trauma-related human T cell and monocyte pathways.

Monitoring genome-wide, cell-specific responses to human disease, although challenging, holds great promise for the future of medicine. Patients with injuries severe enough to develop multiple organ dysfunction syndrome have multiple immune derangements, including T cell apoptosis and anergy combined with depressed monocyte antigen presentation. Genome-wide expression analysis of highly enriched circulating leukocyte subpopulations, combined with cell-specific pathway analyses, offers an opportunity to discover leukocyte regulatory networks in critically injured patients. Severe injury induced significant changes in T cell (5,693 genes), monocyte (2,801 genes), and total leukocyte (3,437 genes) transcriptomes, with only 911 of these genes common to all three cell populations (12%). T cell-specific pathway analyses identified increased gene expression of several inhibitory receptors (PD-1, CD152, NRP-1, and Lag3) and concomitant decreases in stimulatory receptors (CD28, CD4, and IL-2Ralpha). Functional analysis of T cells and monocytes confirmed reduced T cell proliferation and increased cell surface expression of negative signaling receptors paired with decreased monocyte costimulation ligands. Thus, genome-wide expression from highly enriched cell populations combined with knowledge-based pathway analyses leads to the identification of regulatory networks differentially expressed in injured patients. Importantly, application of cell separation, genome-wide expression, and cell-specific pathway analyses can be used to discover pathway alterations in human disease.

Microfluidic isolation of leukocytes from whole blood for phenotype and gene expression analysis.

Technologies that enable the isolation of cell subtypes from small samples of complex populations will greatly facilitate the implementation of proteomics and genomics to human diseases. Transcriptome analysis of blood requires the depletion of contaminating erythrocytes. We report an automated microfluidic device to rapidly deplete erythrocytes from whole blood via deionized water lysis and to collect enriched leukocytes for phenotype and genomic analyses. Starting with blood from healthy subjects, we demonstrate the utility of this microfluidic cassette and lysis protocol to prepare unstimulated leukocytes, and leukocytes stimulated ex vivo with Staphylococcal enterotoxin B, which mimics some of the cellular effects seen in patients with severe bacterial infections. Microarrays are used to assess the global gene expression response to enterotoxin B. The results demonstrate that this system can isolate unactivated leukocytes from small blood samples without any significant loss, which permits more information to be obtained from subsequent analysis, and will be readily applicable to clinical settings.

The number of genes changing expression after chronic exposure to code division multiple access or frequency DMA radiofrequency radiation does not exceed the false-positive rate.

Experiments with cultured C3H 10T 1/2 cells were performed to determine if exposure to cell phone radiofrequency (RF) radiations induce changes in gene expression. Following a 24 h exposure of 5 W/kg specific adsorption rate, RNA was extracted from the exposed and sham control cells for microarray analysis on Affymetrix U74Av2 Genechips. Cells exposed to 0.68 Gy of X-rays with a 4-h recovery were used as positive controls. The number of gene expression changes induced by RF radiation was not greater than the number of false positives expected based on a sham versus sham comparison. In contrast, the X-irradiated samples showed higher numbers of probe sets changing expression level than in the sham versus sham comparison.

Commonality and differences in leukocyte gene expression patterns among three models of inflammation and injury.

The aim of this study was to compare gene expression profiles of leukocytes from blood (white blood cells; WBCs) and spleen harvested at an early time point after injury or sham injury in mice subjected to trauma/hemorrhage, burn injury, or lipopolysaccharide (LPS) infusion at three experimental sites. Groups of injured or LPS-infused animals and sham controls were killed at 2 h after injury and resuscitation, blood and spleen were harvested, and leukocyte populations were recovered after erythrocyte lysis. RNA was extracted from postlysis leukocyte populations. Complementary RNA was synthesized from each RNA sample and hybridized to microarrays. A large number (500-1,400) of genes were differentially expressed at the 2-h time point in injured or LPS-infused vs. sham animals. Thirteen of the differentially expressed genes in blood, and 46 in the spleen, were upregulated or downregulated in common among all three animal models and may represent a common, early transcriptional response to systemic inflammation from a variety of causes. The majority of these genes could be assigned to pathways involved in the immune response and cell death. The up- or downregulation of a cohort of 23 of these genes was validated by RT-PCR. This large-scale microarray analysis shows that, at the 2-h time point, there is marked alteration in leukocyte gene expression in three animal models of injury and inflammation. Although there is some commonality among the models, the majority of the differentially expressed genes appear to be uniquely associated with the type of injury and/or the inflammatory stimulus.

Expression of the proto-oncogene Fos after exposure to radiofrequency radiation relevant to wireless communications.

In this study the expression levels of the proto-oncogene Fos were measured after exposure to radiofrequency (RF) radiation at two relatively high specific absorption rates (SARs) of 5 and 10 W/kg for three types of modulated signals: 847.74 MHz code division multiple access (CDMA), 835.62 MHz frequency division multiple access (FDMA), and 836.55 MHz time division multiple access (TDMA). This work was undertaken to confirm a previous report by Goswami et al. (Radiat. Res. 151, 300-309, 1999) that CDMA and FDMA radiation caused small but statistically significant increases in Fos levels as cells entered plateau phase during exposure. No effects on Myc or Jun levels were observed in that study. Therefore, in the present study, analyses were restricted to Fos expression during the transition from exponential growth to plateau phase. Fos expression was measured using the real-time polymerase chain reaction (RT-PCR) technique. Serum-stimulated C3H 10T(1/2) cells were used as a positive control for Fos expression. Possible influences of final cell number or pH variability on Fos expression were evaluated. Expression of Fos mRNA in C3H 10T(1/2) cells was not significantly different from that found after sham exposure at either SAR level for any signal modulation. Therefore, the results of Goswami et al. could not be confirmed.

A network-based analysis of systemic inflammation in humans.

Oligonucleotide and complementary DNA microarrays are being used to subclassify histologically similar tumours, monitor disease progress, and individualize treatment regimens. However, extracting new biological insight from high-throughput genomic studies of human diseases is a challenge, limited by difficulties in recognizing and evaluating relevant biological processes from huge quantities of experimental data. Here we present a structured network knowledge-base approach to analyse genome-wide transcriptional responses in the context of known functional interrelationships among proteins, small molecules and phenotypes. This approach was used to analyse changes in blood leukocyte gene expression patterns in human subjects receiving an inflammatory stimulus (bacterial endotoxin). We explore the known genome-wide interaction network to identify significant functional modules perturbed in response to this stimulus. Our analysis reveals that the human blood leukocyte response to acute systemic inflammation includes the transient dysregulation of leukocyte bioenergetics and modulation of translational machinery. These findings provide insight into the regulation of global leukocyte activities as they relate to innate immune system tolerance and increased susceptibility to infection in humans.

Application of genome-wide expression analysis to human health and disease.

The application of genome-wide expression analysis to a large-scale, multicentered program in critically ill patients poses a number of theoretical and technical challenges. We describe here an analytical and organizational approach to a systematic evaluation of the variance associated with genome-wide expression analysis specifically tailored to study human disease. We analyzed sources of variance in genome-wide expression analyses performed with commercial oligonucleotide arrays. In addition, variance in gene expression in human blood leukocytes caused by repeated sampling in the same subject, among different healthy subjects, among different leukocyte subpopulations, and the effect of traumatic injury, were also explored. We report that analytical variance caused by sample processing was acceptably small. Blood leukocyte gene expression in the same individual over a 24-h period was remarkably constant. In contrast, genome-wide expression varied significantly among different subjects and leukocyte subpopulations. Expectedly, traumatic injury induced dramatic changes in apparent gene expression that were greater in magnitude than the analytical noise and interindividual variance. We demonstrate that the development of a nation-wide program for gene expression analysis with careful attention to analytical details can reduce the variance in the clinical setting to a level where patterns of gene expression are informative among different healthy human subjects, and can be studied with confidence in human disease.

Whole blood and leukocyte RNA isolation for gene expression analyses.

The analysis of gene expression data in clinical medicine has been plagued by the lack of a critical evaluation of accepted methodologies for the collection, processing, and labeling of RNA. In the present report, the reliability of two commonly used techniques to isolate RNA from whole blood or its leukocyte compartment was compared by examining their reproducibility, variance, and signal-to-noise ratios. Whole blood was obtained from healthy subjects and was either untreated or stimulated ex vivo with Staphylococcus enterotoxin B (SEB). Blood samples were also obtained from trauma patients but were not stimulated with SEB ex vivo. Total RNA was isolated from whole blood with the PAXgene proprietary blood collection system or from isolated leukocytes. Biotin-labeled cRNA was hybridized to Affymetrix GeneChips. The Pearson correlation coefficient for gene expression measurements in replicates from healthy subjects with both techniques was excellent, exceeding 0.985. Unsupervised analyses, including hierarchical cluster analysis, however, revealed that the RNA isolation method resulted in greater differences in gene expression than stimulation with SEB or among different trauma patients. The intraclass correlation, a measure of signal-to-noise ratio, of the difference between SEB-stimulated and unstimulated blood from healthy subjects was significantly higher in leukocyte-derived samples than in whole blood: 0.75 vs. 0.46 (P = 0.002). At the P < 0.001 level of significance, twice as many probe sets discriminated between SEB-stimulated and unstimulated blood with leukocyte isolation than with PAXgene. The findings suggest that the method of RNA isolation from whole blood is a critical variable in the design of clinical studies using microarray analyses.

Expression profiling using random genomic DNA microarrays identifies differentially expressed genes associated with three major developmental stages of the protozoan parasite Leishmania major.

To complete its life cycle, protozoan parasites of the genus Leishmania undergo at least three major developmental transitions. However, previous efforts to identify genes showing stage regulated changes in transcript abundance have yielded relatively few. Here we used expression profiling to assess changes in transcript abundance in three stages: replicating promastigotes and infective non-replicating metacyclics, which occur in the sand fly vector, and in the amastigote stage residing with macrophage phagolysosomes in mammals. Microarrays were developed containing 11,484 PCR products that included a number of known genes and 10,464 random 1 kb genomic DNA fragments. Arrays were hybridized in triplicate and genes showing two-fold or greater changes in 2/3 experiments were scored as differentially expressed. Remarkably, only about one percent of the DNAs expression varied by this criteria, in either stage comparison. Northern blot analysis confirmed the predicted change in mRNA abundance for most of these (68%). This set of genes included most of those previously identified in the literature as differentially regulated as well as a number of novel genes. Notably, Leishmania maxicircle transcripts showed strong up-regulation in metacyclic and amastigote parasites, probably associated with changes in parasite energy metabolism. However, current data suggest that expression profiling using shotgun DNA libraries significantly underestimates the extent of regulated transcripts.