- Research article
- Open Access
- Open Peer Review
A longitudinal study of gene expression in healthy individuals
BMC Medical Genomicsvolume 2, Article number: 33 (2009)
The use of gene expression in venous blood either as a pharmacodynamic marker in clinical trials of drugs or as a diagnostic test requires knowledge of the variability in expression over time in healthy volunteers. Here we defined a normal range of gene expression over 6 months in the blood of four cohorts of healthy men and women who were stratified by age (22–55 years and > 55 years) and gender.
Eleven immunomodulatory genes likely to play important roles in inflammatory conditions such as rheumatoid arthritis and infection in addition to four genes typically used as reference genes were examined by quantitative reverse transcription-polymerase chain reaction (qRT-PCR), as well as the full genome as represented by Affymetrix HG U133 Plus 2.0 microarrays.
Gene expression levels as assessed by qRT-PCR and microarray were relatively stable over time with ~2% of genes as measured by microarray showing intra-subject differences over time periods longer than one month. Fifteen genes varied by gender. The eleven genes examined by qRT-PCR remained within a limited dynamic range for all individuals. Specifically, for the seven most stably expressed genes (CXCL1, HMOX1, IL1RN, IL1B, IL6R, PTGS2, and TNF), 95% of all samples profiled fell within 1.5–2.5 Ct, the equivalent of a 4- to 6-fold dynamic range. Two subjects who experienced severe adverse events of cancer and anemia, had microarray gene expression profiles that were distinct from normal while subjects who experienced an infection had only slightly elevated levels of inflammatory markers.
This study defines the range and variability of gene expression in healthy men and women over a six-month period. These parameters can be used to estimate the number of subjects needed to observe significant differences from normal gene expression in clinical studies. A set of genes that varied by gender was also identified as were a set of genes with elevated expression in a subject with iron deficiency anemia and another subject being treated for lung cancer.
Gene expression profiling studies in venous blood are used to explore transcriptional differences between diseased and healthy individuals, identify biomarkers that may identify the appropriate therapy for individual patients, monitor a therapy and pharmacodynamic responses to new drugs, or define prognosis. For example, tumor molecular profiles have been identified that select preferred therapy and predict recurrence [1–3]. For autoimmune diseases, where the disease tissue is not readily accessible, peripheral blood has become a surrogate tissue by default . Venous blood contains cells that will migrate into and/or have migrated from the inflamed tissue, potentially allowing the identification of patients with disease by expression profiling of whole blood. To date, studies of rheumatoid arthritis, systemic lupus erythematosus (SLE), Crohn's disease, and multiple sclerosis (MS) patients have identified expression profiles that may be useful in diagnosis [5, 6].
Knowledge of the variability of gene expression over time in healthy men and women will help in the design of studies where drugs are expected to modulate gene expression in the diseased population. It will be particularly important to clearly establish the extent of normal variation for cases where the difference between normal and diseased gene expression profiles is subtle.
Clinical studies need to control for multiple sources of variability since gene expression can be influenced by psychological stress [7, 8], exercise [9, 10] and eating . In vitro, inflammation markers are rapidly induced in peripheral blood mononuclear cells (PBMC) by infectious agents [12–14]. In vivo gene expression profiles of individuals exposed to a bacterial toxin appeared to be different from those of healthy individuals not exposed to the toxin . Gene expression is also sensitive to the method by which the RNA from blood is isolated [16, 17].
Several previous studies have examined gene expression in the peripheral blood of normal individuals [15, 18–21]. These investigators found inter-individual variation as well as genes that uniquely identified the subjects being studied. The results were obtained using either peripheral blood mononuclear cells (PBMC) or whole blood. Gene expression profiles differ greatly among the different blood cell types [22, 23]. However, the isolation of the different cell types requires capable laboratories that are not always available at all centers involved in large clinical trials. Whole blood can be sampled reliably with the PAXgene Blood RNA System [24, 25]. Thus, the studies here were carried out using this collection method.
The gene expression studies mentioned previously were largely limited to a single platform for assays (microarray or qRT-PCR), small numbers of subjects with analyses often restricted to a short sampling time frame, or even a single time point [15, 18–21]. This study represents the largest reported study of longitudinal genomic sampling utilizing both microarray and qRT-PCR in an age and gender stratified healthy patient population, to the best of our knowledge. Moreover, the study was performed under controlled conditions with all samples drawn in the morning from fasted individuals at fixed times relative to the start of the study.
The present study extends the information base for healthy individuals by measuring the stability and dynamic range of gene expression in healthy subjects over an extended period of time, six months. qRT-PCR was used to define the range of gene expression in 28 individuals. Eleven immunomodulatory genes likely to play important roles in inflammatory conditions and four reference genes were examined. In addition, 22 subjects were assessed by high-density oligonucleotide microarray, two of whom experienced a severe adverse event during the six-month study period.
The study was conducted in accordance with good clinical practice, the Declaration of Helsinki, and appropriate regulatory guidelines. An independent ethics committee approved the protocol. Subjects gave written informed consent to participate in the study.
Blood samples were obtained from 80 healthy volunteers enrolled at a single site in Strasbourg, France. All subjects were enrolled in January and early February, 2003 and had their last visit in August. Subjects were enrolled into one of four cohorts:
Cohort 1: males 20–55 years of age
Cohort 2: males > 55 years of age
Cohort 3: females 20–55 years of age
Cohort 4: females > 55 years of age
Subjects were enrolled in the study if they were at least 20 years of age and were healthy, as determined by medical history, physical examination, and standard laboratory test values. Exclusion criteria included abnormal laboratory test results, BMI > 32, hypertension, and a positive urine pregnancy test, or positive urine HCV, HBV or HIV test. Any subject who had a concomitant disease or condition that could interfere with clinical evaluations, was an active smoker, or was taking medication on a regular basis was also excluded. Former smoking status was not captured.
A subset of 20 subjects was randomly selected from this study population for microarray analysis. The subset consisted of five subjects from each cohort. Two additional subjects who suffered from serious adverse events during the study were added subsequently. Samples from those 22 subjects were taken at five time points (Baseline, Day 14, Day 28, Day 90 and Day 180).
Quantitative RT-PCR was originally performed on all 80 subjects at each of the five time points for all analytes described in the study. An analysis of the data revealed that mean gene expression in all analytes varied significantly by time point. Due to flaws in the initial experimental design, we could not discern the degree to which the differences were attributable to a true biological effect (e.g. a seasonal effect of gene expression) or experimental bias (e.g. a "batch" effect).
The qPCR study was therefore repeated with 28 subjects chosen at random from the larger set of 80 where sufficient RNA remained. The 28-subject study was designed to distinguish between biological and experimental variability, and to minimize experimental variability. The 28 subjects were assayed by qRT-PCR at four time points (Baseline, Day 28, Day 90 and Day 180). Samples from eight subjects were analyzed on both qRT-PCR and microarray platforms.
RNA processing and quantitation
Whole blood (2.5 mL) was collected into PAXgene tubes (Becton-Dickinson Diagnostics; Hombrechtikon, Switzerland) and frozen immediately at -20°C. The tubes were shipped frozen and stored at -80°C for a period ranging from one week to six months before RNA was extracted. RNA extractions were performed in batches in the order in which the PAXgene tubes were received. Total RNA was extracted using the PAXgene 96 Blood RNA kit (Qiagen; Hilden, Germany) according to the manufacturer's instructions. RNA quality was assessed on an Agilent Bioanalyzer 2100 (Agilent Technologies; Palo Alto, CA). Any RNA sample with an RNA integrity number (RIN) < 6.5 was rejected and a replacement sample was obtained from a replicate PAXgene tube. The mean yield of RNA from all extractions with RIN ≥ 6.5 was 7.0 μg per tube (1 SD = 2.6 μg).
Quantitative RT-PCR assays used the fluorescent dye SYBR Green to monitor amplicon formation in a one step format in which the RT step was performed in the same tube as the PCR. All 384-well plates were profiled on an ABI Prism 7900 HT Sequence Detection System (Applied Biosystems Inc.; Foster City, CA). All reactions were performed with a single reagent lot. The cycling parameters were as follows: Step 1: 50°C 2 min for Uracyl-N-Glycosylase (UNG) digest followed by 95°C 1 min to destroy UNG activity; Step 2: 60°C 30 minutes for reverse transcription; Step 3: 50 Cycles of two-step PCR. Each cycle included 95°C 30 sec (denaturation) followed by 60°C 30 sec (re-annealing and extension) Step 4: Slow ramp from 60°C to 95°C to collect data for dissociation curves. UNG was used to eliminate the possibility of carryover contamination.
One well of RNA diluent was run as a negative control for each RNA specimen. RNA diluent consisted of 10 mM TRIS pH 8.0, 0.1 mM EDTA, 0.020 mg/mL poly r(A) RNA (GE Healthcare), and 0.09% sodium azide (w/v). "High input" and "low input" positive controls assessing expression of B2M were run on each plate. The target for the positive controls was human blood peripheral leukocyte total RNA (Clontech) (100 ng for high input and 1 ng for low input). After a run, the plate was inspected to see that Ct values for high input and low input controls fell within a specified range.
Samples were quantified using ribogreen, a nucleic acid stain. After the concentration was determined, the samples were diluted in RNA diluent. 2 ng of RNA were used in a 10 μl single-tube RT-PCR reaction for all assays except IL-6, where 20 ng were used. All reactions were set up in triplicate using a Biomek FX Laboratory Automation Workstation (Beckman Coulter; Fullerton, CA). The primers used for each assay and gene accession numbers are shown in Additional File 1.
Whole genome expression profiles were generated for all samples at scheduled time points (Days 1, 14, 28, 90, and 180) in the subset of 22 subjects with the exception of the Day 180 time point for subject 174 who died prior to the end of the study. A total number of 109 microarrays were hybridized. The Day 1, 90, and 180 samples were obtained in January/February, April/May, and July/August, respectively, and the mRNA was extracted in the order received (see "RNA processing and quantitation"). Additionally, the samples were processed for microarrays in two batches. The first batch was with samples up to Day 28, and the second with samples taken at Day 90 and Day 180. One microgram of total RNA from each blood sample was used to generate biotinylated cRNA, using the Agilent Low RNA Input Linear Amplification procedure (Version 2), with minor modifications (Agilent Technologies GmbH; Waldbronn, Germany); modification: biotinylated UTP and CTP (Enzo Life Sciences; Farmingdale, NY) were used for the in vitro transcription reaction, each at a final concentration of 1 mM. cRNA was purified using Promega SV 96 DNA isolation kit (automated method), as per the manufacturer's instructions (Promega Corporation, Madison, WI) or Qiagen's RNeasy purification kit (manual method) (Qiagen GmbH; Hilden, Germany); performance equivalency of the manual versus automated methods were assessed in terms of yield, quality and reproducibility with control universal human reference RNA and control whole blood RNA (data not shown). The cRNA samples were hybridized overnight to Affymetrix U133 Plus 2.0 full genome oligonucleotide arrays and then stained with Streptavidin-Phycoerythrin according to the manufacturer's instructions (Affymetrix Inc, CA, USA). Arrays were scanned using a GeneChip Scanner 3000 (Affymetrix) and signal intensities were calculated automatically by GeneChip Operating Software (GCOS, version 1.0; Affymetrix). Gene signal intensities were computed using the MAS 5.0 algorithm (component of GCOS 1.0 software). Signal intensities were normalized using a quantile-quantile method . All normalized data were log2-transformed prior to analysis to down-weight the influence of high expression values. The data have been deposited in NCBI's Gene Expression Omnibus (Edgar et al., 2002) and are accessible through GEO Series accession number GSE16028 http://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE16028.
Units of cycle threshold (Ct), were used for the expression data for all analyses. Outliers were defined as data points falling ≥ 3 standard deviations (SD) outside the mean Ct value, averaged over all time points and individuals, within a gene. A graphical assessment was performed to look for possible seasonal effects.
In addition to the normalization to total RNA input that was performed, a strategy of normalization to four reference genes as a function of average gene expression weighted by their variances was considered. Several aspects of our study led us to conclude that reference gene normalization was introducing several sources of systematic variability into the data. This variability was estimated to be of appreciable size (analyses not shown) and was not present in the total RNA normalized data because it had been controlled for through study design. This variability was due to differences among reactions and differences across plates. The differences among reactions were due to the necessity of having reference genes measured in external reactions, because SYBR green assays cannot be multiplexed. The differences among plates were due to our 28-subject substudy design where a single gene was assayed per plate, which allowed us to make all analyses on within-plate comparisons only. It is known that the precision that can be achieved with reference gene normalization is limited by the variability present in the reference genes . Since we could not increase the precision of the immunomodulatory gene expression data with the reference genes, and we had alternatively controlled for appreciable amounts of systematic variability through study design, the strategy of normalization to reference genes was abandoned.
A linear mixed model with fixed effects for gender, the interaction of gender and time, and age, and a random effect for subject was applied to the qRT-PCR data. From this model, effects were tested and within-individual variation in gene expression over time was estimated. The analysis was implemented using SAS Proc Mixed.
An exploratory analysis was first performed to assess outliers and to determine if any clusters of subjects could be observed. For this purpose, a correspondence analysis was used. This analysis reduces the complexity of the data and facilitates their interpretation by finding combinations of genes that best explain the variability in the entire data set . The analysis was performed using the statistical package XlStat 6.0 (Addinsoft; New York, NY).
Identification of differentially expressed genes
Of the potential 54,675 probe sets on the U133 Plus 2.0 microarray, 34,573 probe sets were present in at least one sample (among the 109 microarrays hybridized) and were used for the analysis. A statistical analysis was performed to determine if expression profiles could indicate significant gender, age, and time effects on gene expression. A linear mixed-effect model  was built independently for each of the 34,573 selected probe sets using log2-transformed expression values and the analysis was implemented using the function lme in SPLUS.
The model used can be expressed as:
where y is the normalized signal intensity of probe set p in individual i at time t, β gender , β age , and β time are respectively the gender, age and time effects and α is the individual random effect. The variances were denoted (between-subjects variance) for and (within-subject variance) for ε itp .
The significance of the effects was tested with t-tests. Because of the numerous t-tests performed, we needed to use a correction procedure, called the False Discovery Rate (FDR) , which controls the expected proportion of genes erroneously identified as differentially expressed genes. For this study, an FDR of 5% was chosen. For each probe set, a Shapiro test was performed to check whether the linear mixed-effect models correctly fit the data. Only 17,329 probe sets for which the Shapiro P-value > 0.1% were kept for further analysis.
Using the variance components estimated with the linear mixed-effect model built for each of the 17, 329 probe sets, we determined for each probe set p the intraclass correlation coefficient among expression measures within a subject over time by calculating the term , with the variance computed between subjects, and the within-subject variance. Values of the correlation range between 1 (i.e. no variation within subject) and 0 (i.e. high variation within subject).
We used the graphical representation proposed by Bland and Altman  to determine the level of agreement between gene expression values measured at different time points. This type of graph plots the difference between two measures as a function of their means. A good agreement between different measures is observed when the cloud of data points is located around 0 and from the absence of a trend (i.e., slope of a robust linear fit of the scatter plot close to 0). In the graphs shown the differences on the vertical axis and the means on the horizontal axis were not computed from the normalized expression values but from the effects determined by the linear mixed-effect model using time as a factor.
Subjects who experienced an adverse event had an additional sample collected following recovery from their infections. Eighteen of the eighty subjects in the study had a total of 24 infections. Infections in 7 patients were ongoing at the time of a scheduled sample collection, and infections in 11 patients ended within fourteen days before a collection. The other seven infections in these patients occurred outside these time frames and were not included in this analysis. There was a wide range of times between the scheduled blood draw and the recovery blood draw (9–51 days), with most recovery period blood draws occurring between 9–22 days. The time between the end of the infection and the recovery period blood draw ranged from 14–97 days, with most occurring between 14 to 23 days after the end of the event.
Results and Discussion
An overview of the study population and hematology measures is presented first, followed by gene expression data from both individuals who remained healthy throughout the study, and those who experienced an adverse event.
Study population characteristics
The healthy volunteers, 79 Caucasians and 1 Oriental, were from eastern France and stratified into four cohorts of 20 individuals each based on age (22 – 55 years and > 55 years) and gender (Table 1). The average age of males versus females in each of the two age groups was similar: 34.5 ± 7.8 (mean ± SD) years for males versus 32.3 ± 9.5 years for females in the age 20 – 55 years group, and 59.7 ± 3.1 versus 59.5 ± 4.5, respectively, in the over 55 group. Patients were excluded from the trial if they took medication on a regular basis including birth control, hormone replacement therapy, or herbal medicines or used tobacco. The most frequently used medication for adverse events that occurred during the study was paracetamol (30 subjects), followed by hexetedine (4), phloroglucinol (4), and amoxicylline (4). One subject died from lung cancer prior to the end of the study, while all other subjects completed the five scheduled visits on Days 1, 14, 28, 90, and 180 per protocol.
A subset of 28 individuals was chosen from the larger study population and analyzed for the expression of selected genes by qRT-PCR at four time points (Baseline, Day 28, Day 90 and Day 180) (Table 1). Gene expression of twenty-two subjects was assessed by high-density microarray. Samples from eight subjects were analyzed on both qRT-PCR and microarray platforms. All other samples in the two subsets were assayed on only one platform.
Complete blood cell counts were obtained from all 80 subjects at all five time points. Neutrophils and lymphocytes, which together account for ~90% of white blood cells, had similar distributions (Table 2 and Additional File 2). Most subjects had cell counts that were in the lower half of the normal laboratory range for each cell type. The within subject estimate of variability for neutrophils was 21% of the mean, and was 16% for lymphocytes. The ratio of the top to the bottom of the range of values measured in the study was 6.3 for neutrophils and 5.7 for lymphocytes.
Range of expression in selected immunomodulatory and reference genes by qRT-PCR
The RNA expression of eleven inflammatory markers (CXCL1, HMOX1, ICAM1, IL1B, IL1RN, IL6R, MMP9, PTGS2, SERPINE1, TGFB1, TNF), and four genes commonly used as reference genes in blood [32, 33] (B2M, 18S rRNA, GAPDH and PPP1CA), was evaluated by qRT-PCR (Figure 1). The dynamic range of expression for all the genes, both inflammatory and reference, was relatively narrow. More specifically, 95% of all samples profiled for CXCL1, HMOX1, IL1B, IL1RN, IL6R, PTGS2 and TNF had a range of gene expression that fell within 1.5–2.5 Ct, an approximately 4–6 fold range in transcript number, similar to the range for neutrophil and lymphocyte cell numbers. The gene that exhibited the highest variability was TGFB1, in which 95% of all samples spanned 4.8 Ct, or a 28-fold range. Shapiro-Wilkes tests and Quantile-Quantile analysis demonstrated that expression was approximately log-normally distributed for most targets, with slight deviations in the tails (data not shown). Median expression values for each immunomodulatory gene evaluated remained stable across each of the four time points (Figure 2). Each qRT-PCR assay was performed using an equivalent amount of RNA from all samples as input. Normalizing against any of the four reference genes was not effective at increasing precision of the observed inflammatory gene expression (see Additional File 3).
McLoughlin et al. evaluated expression of 48 inflammation and immune-related transcripts by qRT-PCR in healthy individuals at a single time point, including 10 of the genes in the present study. They found a similarly narrow range of expression between subjects, including standard deviation values that were within 0.2 Ct of that observed here for most targets examined. Two of the targets with the highest variability in our study, MMP9 and SERPINE1, showed similarly higher variability in that study.
While PPP1CA and B2M had a relatively narrow range of expression, GAPDH and 18S rRNA were more variable among individuals than several of the immunomodulatory genes and therefore were not suitable for normalizing the data set. GAPDH is elevated in T cells upon activation [32, 34] and can fluctuate in response to the changing energy demands of cells . The results presented here support the notion that normalizing qRT-PCR data by input RNA amount is a reasonable approach and may introduce less variability than the use of a single reference gene.
Variation in gene expression over time in samples hybridized on microarrays
Pairwise comparisons were performed to determine whether the mean expression of each probe set on the array varied between any two time points. None of the differences seen at Days 14 and 28 compared to Day 1 remained significant after correcting for multiple comparisons (Figure 3). The Bland-Altman plot  (Figure 4) shows there is no difference for samples separated by two weeks (e.g. when Day 14 is compared to Day 1 or when Day 28 is compared to Day 14). However, when Day 90 was compared to Day 180, 248 probe sets were found to be differentially expressed, a higher number than would be expected by random chance. Those 248 probe sets correspond to 157 unique genes (excluding probe sets which are not correctly mapped to the genome). Among those 157 genes, 66 genes were found associated with apoptotic activities with log2-fold change ranging from -0.9 to 0.6. It was not possible to draw a conclusion from the comparison of the gene signal intensities measured at Day 1 and at Day 90 (and at Day 1 and Day 180) as samples were processed separately in two batches. The differences we found may be the result of a seasonal effect, a true biological effect, or experimental bias. The observed stability of the gene expression profiles over a one-month period is consistent with results from other studies [19, 21].
Within-individual variability in gene expression over time
The estimate of total within-individual variance of the eleven inflammatory genes shown in Table 3 reflects the amount of deviation around the mean trend of gene expression for a target gene over time. Total within-individual variance estimates for the eleven inflammatory genes ranged from 0.220 to 0.838 Ct 2, which was 1–3% of the mean gene expression for a given gene. The two genes with the highest within-individual variance estimates were MMP9 and SERPINE1, while TNF had the smallest. From the mixed effects models for each mRNA measured by qRT-PCR, it was estimated that the correlation of gene expression among time points in a single individual is less than 50% for all genes assayed.
One measure of the temporal variability of a gene is how well its expression is correlated in any individual from one point in time to the next. Using the variance components estimated with the linear mixed effect models, we determined the correlation coefficient among expression measures within a subject over time for each of the 17,329 probe sets satisfying the Shapiro test (see Methods). Only 10% of the probe sets exhibited a coefficient > 0.4 (Figure 5). However, when only probe sets demonstrating robust transcriptional expression with a log2-transformed signal intensity ≥ 7 (raw signal intensity > 100) were considered, the percentage with an intraclass correlation coefficient > 0.4 increased to 23%. The observation that those genes expressed at low levels had the lowest intra-individual correlations may be partly due to the known characteristic of probe sets at the lower end of the quantitative range of a microarray to have poorer reproducibility . For example, Dobbin et al.  found that RNA profiles from tumor tissues had weaker correlations within an individual for those genes with low levels of expression. Inspection of the probe sets with the lowest correlation coefficient (< 0.1) but still with a log2 expression signal intensity ≥ 10 showed that ~90% either do not map to a unique gene locus or do not map to the current version of the human genome, as already observed by Zhang .
Age and gender effects
No genes analyzed by either qRT-PCR or microarray were observed to be differentially regulated between the two age groups. An age effect has been observed previously with RNA profiling studies [19, 21, 39], and might have been expected here, given the large number of genes measured. One explanation is that the broad age range within a given group in this study may have prevented the detection of such an effect. Another is that blood might be less susceptible to senescence effects when compared to other tissues such as brain or muscle [40, 41].
Females showed a small but statistically significant 1.3–1.5-fold increase in gene expression relative to males for the qRT-PCR targets CXCL1, HMOX1 and ICAM1 (see Additional File 4. On the microarray, a total of 78 unique genes were found to be differentially regulated by gender. As expected, the XIST gene was found to be highly down-regulated in males compared to females, with an average log2-transformed signal intensity equal to 50 in females and roughly 20 times lower expression in males. Fifteen of the 78 genes had a log2-transformed signal intensity ≥ 7.0 (Table 4). For all of these genes, the mean difference between sexes was less than 20%. Four X chromosome genes were identified as differentially expressed. One of these, RPS4X, is known to escape X inactivation  while another, E1F1AX, has family members known to escape X inactivation. Only 23 out of the 78 probe sets mapping on the Y chromosome were present in at least one sample and their log2-transformed signal intensities were always below 7. The majority of the highly expressed probe sets that were differentially regulated by gender are located on autosomal chromosomes (Table 4). Among those genes, flotillin1 (FLOT1), which was down-regulated in males compared to females, was recently identified as an estrogen-responsive gene .
Effect of two severe adverse events on gene expression
A correspondence analysis (see Methods) was performed on a set of 34,573 probe sets for all the samples from the 22 individuals studied with microarrays. This analysis looked for combinations of genes that would explain the overall variability in the data set and identify possible outliers. Almost all individuals across all time points were in close proximity to one another (Figure 6). However, data points from two subjects, numbered 121 and 174, were outside of the main cluster.
Subject 121 had anemia secondary to rectal bleeding with a progressive drop in red cell count from baseline to Day 90. The subject's bleeding was stopped surgically on Day 100; on Day 180 red cell parameters had returned to the normal range. The correspondence analysis data points on Days 1, 14, and 28 are outside of the main cluster with the Day 90 profile reaching a maximum distance from the others. The Day 180 point has returned to normal. RBC-associated genes with elevated expression levels included ferrochetalase, carbonic anhydrase, ALAS2, erythrocyte membrane protein band 4.2, glycophorin A and B, and 2,3-bisphosphoglycerate mutase. Some of the genes over-expressed in subject 121, were correlated previously to reticulocyte expressed genes: SLC4A1, EPP42, BCL2L1 and BNIp3L . ALAS2 and carbonic anhydrase have also been identified as genes that are down-regulated in anemia during acute renal allograft rejection .
Subject 174 developed lung cancer and began chemotherapy on approximately Day 70. She was a strong outlier in the correspondence analysis on Day 90. Her monocyte count on that day was severely depressed to 6.5% of the baseline value (from 0.46 × 109/L to 0.03 × 109/L). All other white cell types were within normal limits on Days 1, 14, 28, and 90. For subject 174, 31 genes were found to be up-regulated. The gene ACRBP was elevated four- fold above the other subjects on Day 90. ACRBP is a member of the cancer/testis family of antigens, is immunogenic, and has been detected in different tumor types . The growth arrest-specific (GAS)2-like 1 gene, an actin-associated protein expressed at high levels in growth-arrested cells , was found to be 1.5-fold higher in subject 174, compared to other individuals at Day 90. The changes in gene expression observed in this subject correlated with the onset of chemotherapy and may be largely attributable to treatment. However, for both these cases data is required in larger numbers of patients to draw firm conclusions.
Effect of infection on gene expression
Eighteen of the eighty subjects in the study had a total of 24 infections. The infections were predominantly upper respiratory infections (sinusitis, rhinitis, nasopharyngitis, or pharyngitis), in addition to gastroenteritis, bronchitis, and skin infection and oral herpes simplex. Immunomodulatory gene expression in subjects with active infection at the time of sampling or within two weeks were slightly elevated relative to the expression levels observed two weeks after recovery and at all other time points (Figure 7). However, the expression of TGFB1 was still increased in samples taken after recovery.
In this study with 80 healthy male and female subjects ≥ 20 years of age, a normal range of gene expression was established for 11 inflammation-related genes and four housekeeping genes. Each gene was stably expressed, independent of age and gender. There was no apparent correlation between values more than three standard deviations from the mean with adverse events or hematology values. Nineteen subjects who had infections at the time of or within two weeks before their scheduled blood draws gave an additional sample two weeks after the resolution of their infections. The samples collected during the active infection showed slight elevations from the expression levels seen in uninfected subjects. However, the samples taken two weeks later showed increased expression of TGFβ1.
Gene expression levels, as measured by microarray analysis, appeared to be constant over one month; however, over three months, a small percentage of genes appeared to vary. We observed that intra-individual gene correlations differ greatly depending on signal intensity. A small proportion of genes were found to be differentially regulated according to gender. Differential gene regulation by age (in subjects 25–55 years of age versus subjects > 55 years of age) was not observed. Elevated expression levels of red blood cell-associated genes were observed in one subject who experienced progressive anemia secondary to blood loss. Included among those genes were ferrochetalase, carbonic anhydrase, ALAS2, erythrocyte membrane protein band 4.2, glycophorin A and B, and 2,3-bisphosphoglycerate mutase.
Use of RNA expression studies in clinical trials
Several lessons can be drawn from this study regarding the design and conduct of future clinical trials that test new pharmacological entities and employ gene expression as an assessment:
Estimates of sample sizes required to achieve a certain power to detect changes in gene expression of a particular magnitude can be calculated from the estimates of variance components given in the results section, results from other studies, and the estimated effect size from in vitro work.
Samples should be analyzed at the end of the experiment rather than on an ongoing basis. A bias was introduced when samples from the present study were processed and analyzed with different reagent lots over time in sequential batches This type of bias, also recently described by Yang et al. , could lead to the identification of many artifacts among differentially expressed gene sets.
Small elevations in inflammatory gene expression produced during upper respiratory infections generally return to normal levels of expression within two weeks after the infection has resolved.
Time and/or seasonal effects may be a factor in trials lasting longer than one month.
Gender and age effects are not likely to be problematic in populations from 20 to 65 years of age. The assessment of gene expression profiles in children or very old individuals might lead to different conclusions.
Perou CM, Sorlie T, Eisen MB, Rijn van de M, Jeffrey SS, Rees CA, Pollack JR, Ross DT, Johnsen H, Akslen LA, et al: Molecular portraits of human breast tumours. Nature. 2000, 406: 747-752. 10.1038/35021093.
van't Veer LJ, Dai H, Vijver van de MJ, He YD, Hart AA, Mao M, Peterse HL, Kooy van der K, Marton MJ, Witteveen AT, et al: Gene expression profiling predicts clinical outcome of breast cancer. Nature. 2002, 415: 530-536. 10.1038/415530a.
Paik S, Shak S, Tang G, Kim C, Baker J, Cronin M, Baehner FL, Walker MG, Watson D, Park T, et al: A multigene assay to predict recurrence of tamoxifen-treated, node-negative breast cancer. N Engl J Med. 2004, 351: 2817-2826. 10.1056/NEJMoa041588.
Burczynski ME, Dorner AJ: Transcriptional profiling of peripheral blood cells in clinical pharmacogenomic studies. Pharmacogenomics. 2006, 7: 187-202. 10.2217/146224188.8.131.52.
Aune TM, Maas K, Parker J, Moore JH, Olsen NJ: Profiles of gene expression in human autoimmune disease. Cell Biochem Biophys. 2004, 40: 81-96. 10.1385/CBB:40:2:081.
Baechler EC, Batliwalla FM, Reed AM, Peterson EJ, Gaffney PM, Moser KL, Gregersen PK, Behrens TW: Gene expression profiling in human autoimmunity. Immunol Rev. 2006, 210: 120-137. 10.1111/j.0105-2896.2006.00367.x.
Bierhaus A, Wolf J, Andrassy M, Rohleder N, Humpert PM, Petrov D, Ferstl R, von Eynatten M, Wendt T, Rudofsky G, et al: A mechanism converting psychosocial stress into mononuclear cell activation. Proc Natl Acad Sci USA. 2003, 100: 1920-1925. 10.1073/pnas.0438019100.
Nagabhushan M, Mathews HL, Witek-Janusek L: Aberrant nuclear expression of AP-1 and NFkappaB in lymphocytes of women stressed by the experience of breast biopsy. Brain Behav Immun. 2001, 15: 78-84. 10.1006/brbi.2000.0589.
Petersen AM, Pedersen BK: The anti-inflammatory effect of exercise. J Appl Physiol. 2005, 98: 1154-1162. 10.1152/japplphysiol.00164.2004.
Sureda A, Ferrer MD, Tauler P, Maestre I, Aguilo A, Cordova A, Tur JA, Roche E, Pons A: Intense physical activity enhances neutrophil antioxidant enzyme gene expression. Immunocytochemistry evidence for catalase secretion. Free Radic Res. 2007, 41: 874-883. 10.1080/10715760701416459.
van Erk MJ, Blom WA, van Ommen B, Hendriks HF: High-protein and high-carbohydrate breakfasts differentially change the transcriptome of human blood cells. Am J Clin Nutr. 2006, 84: 1233-1241.
Boldrick JC, Alizadeh AA, Diehn M, Dudoit S, Liu CL, Belcher CE, Botstein D, Staudt LM, Brown PO, Relman DA: Stereotyped and specific gene expression programs in human innate immune responses to bacteria. Proc Natl Acad Sci USA. 2002, 99: 972-977. 10.1073/pnas.231625398.
Spurrell JC, Wiehler S, Zaheer RS, Sanders SP, Proud D: Human airway epithelial cells produce IP-10 (CXCL10) in vitro and in vivo upon rhinovirus infection. Am J Physiol Lung Cell Mol Physiol. 2005, 289: L85-95. 10.1152/ajplung.00397.2004.
Stockl J, Vetr H, Majdic O, Zlabinger G, Kuechler E, Knapp W: Human major group rhinoviruses downmodulate the accessory function of monocytes by inducing IL-10. J Clin Invest. 1999, 104: 957-965. 10.1172/JCI7255.
McLoughlin K, Turteltaub K, Bankaitis-Davis D, Gerren R, Siconolfi L, Storm K, Cheronis J, Trollinger D, Macejak D, Tryon V, Bevilacqua M: Limited dynamic range of immune response gene expression observed in healthy blood donors using RT-PCR. Mol Med. 2006, 12: 185-195. 10.2119/2006-00018.McLoughlin.
Baechler EC, Batliwalla FM, Karypis G, Gaffney PM, Moser K, Ortmann WA, Espe KJ, Balasubramanian S, Hughes KM, Chan JP, et al: Expression levels for many genes in human peripheral blood cells are highly sensitive to ex vivo incubation. Genes Immun. 2004, 5: 347-353. 10.1038/sj.gene.6364098.
Hartel C, Bein G, Muller-Steinhardt M, Kluter H: Ex vivo induction of cytokine mRNA expression in human blood samples. J Immunol Methods. 2001, 249: 63-71. 10.1016/S0022-1759(00)00334-3.
Campbell C, Vernon SD, Karem KL, Nisenbaum R, Unger ER: Assessment of normal variability in peripheral blood gene expression. Dis Markers. 2002, 18: 201-206.
Eady JJ, Wortley GM, Wormstone YM, Hughes JC, Astley SB, Foxall RJ, Doleman JF, Elliott RM: Variation in gene expression profiles of peripheral blood mononuclear cells from healthy volunteers. Physiol Genomics. 2005, 22: 402-411. 10.1152/physiolgenomics.00080.2005.
Radich JP, Mao M, Stepaniants S, Biery M, Castle J, Ward T, Schimmack G, Kobayashi S, Carleton M, Lampe J, Linsley PS: Individual-specific variation of gene expression in peripheral blood leukocytes. Genomics. 2004, 83: 980-988. 10.1016/j.ygeno.2003.12.013.
Whitney AR, Diehn M, Popper SJ, Alizadeh AA, Boldrick JC, Relman DA, Brown PO: Individuality and variation in gene expression patterns in human blood. Proc Natl Acad Sci USA. 2003, 100: 1896-1901. 10.1073/pnas.252784499.
Cobb JP, Mindrinos MN, Miller-Graziano C, Calvano SE, Baker HV, Xiao W, Laudanski K, Brownstein BH, Elson CM, Hayden DL, et al: Application of genome-wide expression analysis to human health and disease. Proc Natl Acad Sci USA. 2005, 102: 4801-4806. 10.1073/pnas.0409768102.
Palmer C, Diehn M, Alizadeh AA, Brown PO: Cell-type specific gene expression profiles of leukocytes in human peripheral blood. BMC Genomics. 2006, 7: 115-10.1186/1471-2164-7-115.
Fan H, Hegde PS: The transcriptome in blood: challenges and solutions for robust expression profiling. Curr Mol Med. 2005, 5: 3-10. 10.2174/1566524053152861.
Kim SJ, Dix DJ, Thompson KE, Murrell RN, Schmid JE, Gallagher JE, Rockett JC: Effects of storage, RNA extraction, genechip type, and donor sex on gene expression profiling of human whole blood. Clin Chem. 2007, 53: 1038-1045. 10.1373/clinchem.2006.078436.
Bolstad BM, Irizarry RA, Astrand M, Speed TP: A comparison of normalization methods for high density oligonucleotide array data based on variance and bias. Bioinformatics. 2003, 19: 185-193. 10.1093/bioinformatics/19.2.185.
Huggett J, Dheda K, Bustin S, Zumla A: Real-time RT-PCR normalisation; strategies and considerations. Genes Immun. 2005, 6: 279-284. 10.1038/sj.gene.6364190.
Fellenberg K, Hauser NC, Brors B, Neutzner A, Hoheisel JD, Vingron M: Correspondence analysis applied to microarray data. Proc Natl Acad Sci USA. 2001, 98: 10781-10786. 10.1073/pnas.181597298.
Pinheiro JC, Bates DM: Mixed-effects models in S and S-PLUS. 2000, New York: Springer
Benjamini Y, Hochberg Y: Controlling the False Discovery Rate – a Practical and Powerful Approach to Multiple Testing. Journal of the Royal Statistical Society Series B-Methodological. 1995, 57: 289-300.
Bland JM, Altman DG: Statistical methods for assessing agreement between two methods of clinical measurement. Lancet. 1986, 1: 307-310.
Bas A, Forsberg G, Hammarstrom S, Hammarstrom ML: Utility of the housekeeping genes 18S rRNA, beta-actin and glyceraldehyde-3-phosphate-dehydrogenase for normalization in real-time quantitative reverse transcriptase-polymerase chain reaction analysis of gene expression in human T lymphocytes. Scand J Immunol. 2004, 59: 566-573. 10.1111/j.0300-9475.2004.01440.x.
Silver N, Best S, Jiang J, Thein SL: Selection of housekeeping genes for gene expression studies in human reticulocytes using real-time PCR. BMC Mol Biol. 2006, 7: 33-10.1186/1471-2199-7-33.
Radonic A, Thulke S, Mackay IM, Landt O, Siegert W, Nitsche A: Guideline to reference gene selection for quantitative real-time PCR. Biochem Biophys Res Commun. 2004, 313: 856-862. 10.1016/j.bbrc.2003.11.177.
Barber RD, Harmer DW, Coleman RA, Clark BJ: GAPDH as a housekeeping gene: analysis of GAPDH mRNA expression in a panel of 72 human tissues. Physiol Genomics. 2005, 21: 389-395. 10.1152/physiolgenomics.00025.2005.
Li X, Kim J, Zhou J, Gu W, Quigg R: Use of signal thresholds to determine significant changes in microarray data analyses. Genet Mol Biol. 2005, 28: 191-200.
Dobbin KK, Beer DG, Meyerson M, Yeatman TJ, Gerald WL, Jacobson JW, Conley B, Buetow KH, Heiskanen M, Simon RM, et al: Interlaboratory comparability study of cancer gene expression analysis using oligonucleotide microarrays. Clin Cancer Res. 2005, 11: 565-572.
Zhang J, Finney RP, Clifford RJ, Derr LK, Buetow KH: Detecting false expression signals in high-density oligonucleotide arrays by an in silico approach. Genomics. 2005, 85: 297-308. 10.1016/j.ygeno.2004.11.004.
Tang Y, Lu A, Ran R, Aronow BJ, Schorry EK, Hopkin RJ, Gilbert DL, Glauser TA, Hershey AD, Richtand NW, et al: Human blood genomics: distinct profiles for gender, age and neurofibromatosis type 1. Brain Res Mol Brain Res. 2004, 132: 155-167. 10.1016/j.molbrainres.2003.10.014.
Welle S, Brooks AI, Delehanty JM, Needler N, Thornton CA: Gene expression profile of aging in human muscle. Physiol Genomics. 2003, 14: 149-159.
Galvin JE, Ginsberg SD: Expression profiling in the aging brain: a perspective. Ageing Res Rev. 2005, 4: 529-547. 10.1016/j.arr.2005.06.009.
Fisher EM, Beer-Romero P, Brown LG, Ridley A, McNeil JA, Lawrence JB, Willard HF, Bieber FR, Page DC: Homologous ribosomal protein genes on the human X and Y chromosomes: escape from X inactivation and possible implications for Turner syndrome. Cell. 1990, 63: 1205-1218. 10.1016/0092-8674(90)90416-C.
Punyadeera C, Dassen H, Klomp J, Dunselman G, Kamps R, Dijcks F, Ederveen A, de Goeij A, Groothuis P: Oestrogen-modulated gene expression in the human endometrium. Cell Mol Life Sci. 2005, 62: 239-250. 10.1007/s00018-004-4435-y.
Chua MS, Barry C, Chen X, Salvatierra O, Sarwal MM: Molecular profiling of anemia in acute renal allograft rejection using DNA microarrays. Am J Transplant. 2003, 3: 17-22. 10.1034/j.1600-6143.2003.30104.x.
Ono T, Kurashige T, Harada N, Noguchi Y, Saika T, Niikawa N, Aoe M, Nakamura S, Higashi T, Hiraki A, et al: Identification of proacrosin binding protein sp32 precursor as a human cancer/testis antigen. Proc Natl Acad Sci USA. 2001, 98: 3282-3287. 10.1073/pnas.041625098.
Schneider C, King RM, Philipson L: Genes specifically expressed at growth arrest of mammalian cells. Cell. 1988, 54: 787-793. 10.1016/S0092-8674(88)91065-3.
Yang H, Harrington CA, Vartanian K, Coldren CD, Hall R, Churchill GA: Randomization in laboratory procedure is key to obtaining reproducible microarray results. PLoS ONE. 2008, 3: e3724-10.1371/journal.pone.0003724.
The pre-publication history for this paper can be accessed here:http://www.biomedcentral.com/1755-8794/2/33/prepub
We thank Emanuelle Erhart-Gaudel and the staff of the Roche Clinical Pharmacology Unit in Strasbourg, France, for conducting the clinical trial, David Trollinger of Source Precision Medicine for discussions regarding the genes chosen for qRT-PCR, Xiaoning Wu and Sean Chien for providing technical support for the qRT-PCR experiments, and Laurent Essioux and Brian Rhees for comments regarding the statistical analyses. The study was funded by the Pharmaceutical and Diagnostic Divisions of Hoffman-LaRoche. All authors read and approved the final manuscript.
This study was funded by Roche and sponsored by Roche Palo Alto.
CK processed and quantified RNA from PAXgene tubes; developed assays, performed experiments and analyzed data for qRT-PCR, and drafted part of the manuscript. PM performed the microarray experiments, analyzed data and drafted part of the manuscript. CF, MT and MWS performed the microarray experiments. AJ performed statistical analyses on the qRT-PCR data and drafted part of the manuscript. GD performed statistical analyses on the microarray data and drafted part of the manuscript. MN performed statistical analyses on the adverse events and hematology data. LS and TN performed experiments for qRT-PCR. JC conceived of the study and its experimental design, participated in the interpretation of results, and helped prepare the manuscript.
Chris Karlovich, Guillemette Duchateau-Nguyen, Andrea Johnson, Patricia McLoughlin contributed equally to this work.
Electronic supplementary material
Authors’ original submitted files for images
Below are the links to the authors’ original submitted files for images.