Identification and validation of suitable endogenous reference genes for gene expression studies in human peripheral blood

Arsenic may be involved in fluoride-induced bone toxicity through PTH/PKA/AP1 signaling pathway

Chronic exposure to combined fluoride and arsenic continues to be a major public health problem worldwide, affecting thousands of people. In recent years, more and more researchers began to focus on the interaction between the fluorine and the arsenic. In this study, the selected investigation site was located in China. The study group was selected from people living in fluoride-arsenic polluted areas due to burning coal. The total number of participants was 196; including the fluoride-arsenic anomaly group (130) and the fluoride-arsenic normal group (63). By observing the changes in gene and protein expression of PTH/PKA/AP1 signaling pathway, the results show that fluoride can increase the expression levels of PTH, PKA, and AP1, but arsenic can only affect the expression of AP1; fluoride and arsenic have an interaction on the expression of AP1. Further study found that fluoride and arsenic can affect the mRNA expression level of c-fos gene (AP1 family members), and have an interaction on the expression of c-fos, but not c-jun. The results indicate that PTH/PKA/AP1 signaling pathway may play an important role in bone toxicity of fluoride. Arsenic can affect the expression of c-fos, thereby affecting the expression of transcription factor AP1, indirectly involved in fluoride-induced bone toxicity.

Gene expression profiling of blood in brain arteriovenous malformation patients

Brain arteriovenous malformations (BAVMs) are an important cause of intracranial hemorrhage (ICH) in young adults. Gene expression profiling of blood has led to the identification of stroke biomarkers, and may help identify BAVM biomarkers and illuminate BAVM pathogenesis. It is unknown whether blood gene expression profiles differ between 1) BAVM patients and healthy controls, or 2) unruptured and ruptured BAVM patients at presentation. We characterized blood transcriptional profiles in 60 subjects (20 unruptured BAVM, 20 ruptured BAVM, and 20 healthy controls) using Affymetrix whole genome expression arrays. Expression differences between groups were tested by ANOVA, adjusting for potential confounders. Genes with absolute fold change ≥ 1.2 (false discovery rate corrected p ≤ 0.1) were selected as differentially expressed and evaluated for over-representation in KEGG biological pathways (p ≤ 0.05). Twenty-nine genes were differentially expressed between unruptured BAVM patients and controls, including 13 which may be predictive of BAVM. Patients with ruptured BAVM compared to unruptured BAVM differed in expression of 1490 genes, with over-representation of genes in 8 pathways including MAPK, VEGF, Wnt signaling and several inflammatory pathways. These results suggest clues to the pathogenesis of BAVM and/or BAVM rupture and point to potential biomarkers or new treatment targets.

The expression of inflammatory genes is upregulated in peripheral blood of patients with type 1 diabetes

OBJECTIVE

Our previous gene expression microarray studies identified a number of genes differentially expressed in patients with type 1 diabetes (T1D) and islet autoantibody-positive subjects. This study was designed to validate these gene expression changes in T1D patients and to identify gene expression changes in diabetes complications. RESEARCH DESIGH AND METHODS: We performed high-throughput real-time RT-PCR to validate gene expression changes in peripheral blood mononuclear cells (PBMCs) from a large sample set of 928 T1D patients and 922 control subjects.

RESULTS

Of the 18 genes analyzed here, eight genes (S100A8, S100A9, MNDA, SELL, TGFB1, PSMB3, CD74, and IL12A) had higher expression and three genes (GNLY, PSMA4, and SMAD7) had lower expression in T1D patients compared with control subjects, indicating that genes involved in inflammation, immune regulation, and antigen processing and presentation are significantly altered in PBMCs from T1D patients. Furthermore, one adhesion molecule (SELL) and three inflammatory genes mainly expressed by myeloid cells (S100A8, S100A9, and MNDA) were significantly higher in T1D patients with complications (odds ratio [OR] 1.3-2.6, adjusted P value = 0.005-10(-8)), especially those patients with neuropathy (OR 4.8-7.9, adjusted P value <0.005).

CONCLUSIONS

These findings suggest that inflammatory mediators secreted mainly by myeloid cells are implicated in T1D and its complications.

High-throughput identification of reference genes for research and clinical RT-qPCR analysis of breast cancer samples

BACKGROUND

Quantification and normalization of RT-qPCR data critically depends on the expression of so called reference genes. Our goal was to develop a strategy for the selection of reference genes that utilizes microarray data analysis and combines known approaches for gene stability evaluation and to select a set of appropriate reference genes for research and clinical analysis of breast samples with different receptor and cancer status using this strategy.

METHODS

A preliminary search of reference genes was based on high-throughput analysis of microarray datasets. The final selection and validation of the candidate genes were based on the RT-qPCR data analysis using several known methods for expression stability evaluation: comparative ∆Ct method, geNorm, NormFinder and Haller equivalence test.

RESULTS

A set of five reference genes was identified: ACTB, RPS23, HUWE1, EEF1A1 and SF3A1. The initial selection was based on the analysis of publically available well-annotated microarray datasets containing different breast cancers and normal breast epithelium from breast cancer patients and epithelium from cancer-free patients. The final selection and validation were performed using RT-qPCR data from 39 breast cancer biopsy samples. Three genes from the final set were identified by the means of microarray analysis and were novel in the context of breast cancer assay. We showed that the selected set of reference genes is more stable in comparison not only with individual genes, but also with a system of reference genes used in commercial OncotypeDX test.

CONCLUSION

A selection of reference genes for RT-qPCR can be efficiently performed by combining a preliminary search based on the high-throughput analysis of microarray datasets and final selection and validation based on the analysis of RT-qPCR data with a simultaneous examination of different expression stability measures. The identified set of reference genes proved to be less variable and thus potentially more efficient for research and clinical analysis of breast samples comparing to individual genes and the set of reference genes used in OncotypeDX assay.

Gene expression analysis of peripheral blood cells reveals Toll-like receptor pathway deregulation in colorectal cancer

Colorectal cancer is the leading cause of cancer-related deaths worldwide. The disease is curable when detected at an early stage. However, the compliance rate with current screening recommendations remains poor. An accurate, minimally invasive blood test that has the potential for greater patient compliance would be a welcome addition to the current methods. Recent data have shown that gene expression profile of peripheral blood cells can reflect disease states and thus have diagnostic value. In this study, genome-wide gene expression profiling of peripheral blood cells from 20 healthy controls and 20 colorectal cancer patients were performed using PAXgene™ technology and Affymetrix GeneChip® microarrays. We identified a list of 1,469 genes that were differentially expressed between the healthy controls and cancer patients. Gene annotation and functional enrichment analysis revealed that those genes are mainly related to immune functions. Particularly, a set of genes belonging to the Toll-Like Receptor pathways were up-regulated in the colorectal cancer patients. These findings provide a new understanding of blood gene expression profile in colorectal cancer. Our result may serve as the basis for further development of blood biomarkers for the diagnosis and treatment of colorectal cancer.

Identification and validation of a blood-based 18-gene expression signature in colorectal cancer

PURPOSE

The early detection of colorectal cancer (CRC) is crucial for successful treatment and patient survival. However, compliance with current screening methods remains poor. This study aimed to identify an accurate blood-based gene expression signature for CRC detection.

EXPERIMENTAL DESIGN

Gene expression in peripheral blood samples from 216 patients with CRC tumors and 187 controls was investigated in the study. We first conducted a microarray analysis to select candidate genes that were significantly differentially expressed between patients with cancer and controls. A quantitative reverse transcription PCR assay was then used to evaluate the expression of selected genes. A gene expression signature was identified using a training set (n = 200) and then validated using an independent test set (n = 160).

RESULTS

We identified an 18-gene signature that discriminated the patients with CRC from controls with 92% accuracy, 91% sensitivity, and 92% specificity. The signature performance was further validated in the independent test set with 86% accuracy, 84% sensitivity, and 88% specificity. The area under the receiver operating characteristics curve was 0.94. The signature was shown to be enriched in genes related to immune functions.

CONCLUSIONS

This study identified an 18-gene signature that accurately discriminated patients with CRC from controls in peripheral blood samples. Our results prompt the further development of blood-based gene expression biomarkers for the diagnosis and early detection of CRC.

Identification of Appropriate Reference Genes for qRT-PCR Analysis of Heat-Stressed Mammary Epithelial Cells in Riverine Buffaloes (Bubalus bubalis)

Gene expression studies require appropriate normalization methods for proper evaluation of reference genes. To date, not many studies have been reported on the identification of suitable reference genes in buffaloes. The present study was undertaken to determine the panel of suitable reference genes in heat-stressed buffalo mammary epithelial cells (MECs). Briefly, MEC culture from buffalo mammary gland was exposed to 42 °C for one hour and subsequently allowed to recover at 37 °C for different time intervals (from 30 m to 48 h). Three different algorithms, geNorm, NormFinder, and BestKeeper softwares, were used to evaluate the stability of 16 potential reference genes from different functional classes. Our data identified RPL4, EEF1A1, and RPS23 genes to be the most appropriate reference genes that could be utilized for normalization of qPCR data in heat-stressed buffalo MECs.

Novel biomarkers in type 1 diabetes

Biomarkers are useful tools for research into type 1 diabetes (T1D) for a number of purposes, including elucidation of disease pathogenesis, risk prediction, and therapeutic monitoring. Susceptibility genes and islet autoantibodies are currently the most useful biomarkers for T1D risk prediction. However, these markers do not fully meet the needs of scientists and physicians for several reasons. First, improvement of the specificity and sensitivity is still desirable to achieve better positive predictive values. Second, autoantibodies appear relatively late in the disease process, thus limiting their value in early disease prediction. Third, the currently available biomarkers are not useful for assessing therapeutic outcomes because some are not involved in the disease process (autoantibodies) and others do not change during disease progression (susceptibility genes). Therefore, considerable effort has been devoted to the discovery of novel T1D biomarkers in the last three decades. The advent of high-throughput technologies for genetic, transcriptomic, and proteomic studies has allowed genome-wide examinations of genetic polymorphisms, global gene changes, and protein expression changes in T1D patients and prediabetic subjects. These large-scale studies resulted in the discovery of a large number of susceptibility genes and changes in gene and protein expression. While these studies have provided a number of novel biomarker candidates, their clinical benefits remain to be evaluated in prospective studies, and no new "star biomarker" has been identified until now. Previous studies suggest that significant improvements in study design and analytical methodologies have to be made to identify clinically relevant biomarkers. In this review, we discuss progress, opportunities, challenges, and future directions in the development of T1D biomarkers, mainly by focusing on the genetic, transcriptomic, and proteomic aspects.

The expression of caspases is enhanced in peripheral blood mononuclear cells of autism spectrum disorder patients

Autism and autism spectrum disorders (ASDs) are heterogeneous complex neuro-developmental disorders characterized by dysfunctions in social interaction and communication skills. Their pathogenesis has been linked to interactions between genes and environmental factors. Consistent with the evidence of certain similarities between immune cells and neurons, autistic children also show an altered immune response of peripheral blood mononuclear cells (PBMCs). In this study, we investigated the activation of caspases, cysteinyl aspartate-specific proteases involved in apoptosis and several other cell functions in PBMCs from 15 ASD children compared to age-matched normal healthy developing controls. The mRNA levels for caspase-1, -2, -4, -5 were significantly increased in ASD children as compared to healthy subjects. Protein levels of Caspase-3, -7, -12 were also increased in ASD patients. Our data are suggestive of a possible role of the caspase pathway in ASD clinical outcome and of the use of caspase as potential diagnostic and/or therapeutic tools in ASD management.

Y chromosome gene expression in the blood of male patients with ischemic stroke compared with male controls

BACKGROUND

Sex is suggested to be an important determinant of ischemic stroke risk factors, etiology, and outcome. However, the basis for this remains unclear. The Y chromosome is unique in males. Genes expressed in males on the Y chromosome that are associated with stroke may be important genetic contributors to the unique features of males with ischemic stroke, which would be helpful for explaining sex differences observed between men and women.

OBJECTIVE

We compared Y chromosome gene expression in males with ischemic stroke and male controls.

METHODS

Blood samples were obtained from 40 male patients ≤3, 5, and 24 hours after ischemic stroke and from 41 male controls (July 2003-April 2007). RNA was isolated from blood and was processed using Affymetrix Human U133 Plus 2.0 expression arrays (Affymetrix Inc., Santa Clara, California). Y chromosome genes differentially expressed between male patients with stroke and male control subjects were identified using an ANCOVA adjusted for age and batch. A P < 0.05 and a fold change >1.2 were considered significant.

RESULTS

Seven genes on the Y chromosome were differentially expressed in males with ischemic stroke compared with controls. Five of these genes (VAMP7, CSF2RA, SPRY3, DHRSX, and PLCXD1) are located on pseudoautosomal regions of the human Y chromosome. The other 2 genes (EIF1AY and DDX3Y) are located on the nonrecombining region of the human Y chromosome. The identified genes were associated with immunology, RNA metabolism, vesicle fusion, and angiogenesis.

CONCLUSIONS

Specific genes on the Y chromosome are differentially expressed in blood after ischemic stroke. These genes provide insight into potential molecular contributors to sex differences in ischemic stroke.

Using control genes to correct for unwanted variation in microarray data

Microarray expression studies suffer from the problem of batch effects and other unwanted variation. Many methods have been proposed to adjust microarray data to mitigate the problems of unwanted variation. Several of these methods rely on factor analysis to infer the unwanted variation from the data. A central problem with this approach is the difficulty in discerning the unwanted variation from the biological variation that is of interest to the researcher. We present a new method, intended for use in differential expression studies, that attempts to overcome this problem by restricting the factor analysis to negative control genes. Negative control genes are genes known a priori not to be differentially expressed with respect to the biological factor of interest. Variation in the expression levels of these genes can therefore be assumed to be unwanted variation. We name this method "Remove Unwanted Variation, 2-step" (RUV-2). We discuss various techniques for assessing the performance of an adjustment method and compare the performance of RUV-2 with that of other commonly used adjustment methods such as Combat and Surrogate Variable Analysis (SVA). We present several example studies, each concerning genes differentially expressed with respect to gender in the brain and find that RUV-2 performs as well or better than other methods. Finally, we discuss the possibility of adapting RUV-2 for use in studies not concerned with differential expression and conclude that there may be promise but substantial challenges remain.

Detection of immune cell response to M. tuberculosis-specific antigens by quantitative polymerase chain reaction

One third of the world's population is latently infected with Mycobacterium tuberculosis (Mtb) and up to 10% of infected individuals develop active tuberculosis (TB) in their lifetime. Among the major challenges in the control of TB is the implementation of sensitive methods for detection of latent tuberculosis infection (LTBI). Currently, in vitro interferon gamma release assays, yielding single value readout, are used as an alternative to the traditional tuberculin skin test for the diagnosis of LTBI. More complex characterization of immune status of LTBI individuals, however, is desirable for indication of LTBI subjects for preventative chemotherapy. Here we describe a quantitative polymerase chain reaction (qPCR) for determination of expression levels of 14 genes, additional to interferon gamma, which was applied for comparison of the specific Mtb-antigen immune response of blood cells from healthy, latently infected, and TB individuals. With the use of principal component analysis and discriminant analysis, a pattern of mRNA levels of 6 genes was identified, allowing discrimination of healthy individuals from active TB and LTBI subjects. These results open the way to development of multimarker qPCR for the detection of LTBI.

Evaluation of internal reference genes for quantitative expression analysis by real-time PCR in ovine whole blood

The use of reference genes is commonly accepted as the most reliable approach to normalize qRT-PCR and to reduce possible errors in the quantification of gene expression. The most suitable reference genes in sheep have been identified for a restricted range of tissues, but no specific data on whole blood are available. The aim of this study was to identify a set of reference genes for normalizing qRT-PCR from ovine whole blood. We designed 11 PCR assays for commonly employed reference genes belonging to various functional classes and then determined their expression stability in whole blood samples from control and disease-stressed sheep. SDHA and YWHAZ were considered the most suitable internal controls as they were stably expressed regardless of disease status according to both geNorm and NormFinder software; furthermore, geNorm indicated SDHA/HPRT, YWHAZ/GAPDH and SDHA/YWHAZ as the best reference gene combinations in control, disease-stressed and combined sheep groups, respectively. Our study provides a validated panel of optimal control genes which may be useful for the identification of genes differentially expressed by qRT-PCR in a readily accessible tissue, with potential for discovering new physiological and disease markers and as a tool to improve production traits (e.g., by identifying expression Quantitative Trait Loci). An additional outcome of the study is a set of intron-spanning primer sequences suitable for gene expression experiments employing SYBR Green chemistry on other ovine tissues and cells.

The X-chromosome has a different pattern of gene expression in women compared with men with ischemic stroke

BACKGROUND AND PURPOSE

Differences in ischemic stroke between men and women have been mainly attributed to hormonal effects. However, sex differences in immune response to ischemia may exist. We hypothesized that differential expression of X-chromosome genes in blood immune cells contribute to differences between men and women with ischemic stroke.

METHODS

RNA levels of 683 X-chromosome genes were measured on Affymetrix U133 Plus2.0 microarrays. Blood samples from patients with ischemic stroke were obtained at ≤ 3 hours, 5 hours, and 24 hours (n=61; 183 samples) after onset and compared with control subjects without symptomatic vascular diseases (n=109). Sex difference in X-chromosome gene expression was determined using analysis of covariance (false discovery rate ≤ 0.05, fold change ≥ 1.2).

RESULTS

At ≤ 3, 5, and 24 hours after stroke, there were 37, 140, and 61 X-chromosome genes, respectively, that changed in women; and 23, 18, and 31 X-chromosome genes that changed in men. Female-specific genes were associated with post-translational modification, small-molecule biochemistry, and cell-cell signaling. Male-specific genes were associated with cellular movement, development, cell-trafficking, and cell death. Altered sex specific X-chromosome gene expression occurred in 2 genes known to be associated with human stroke, including galactosidase A and IDS, mutations of which result in Fabry disease and Hunter syndrome, respectively.

CONCLUSIONS

There are differences in X-chromosome gene expression between men and women with ischemic stroke. Future studies are needed to decipher whether these differences are associated with sexually dimorphic immune response, repair or other mechanisms after stroke, or whether some of them represent risk determinants.

Evaluation of Normalization Methods for RNA-Seq Gene Expression Estimation

Statistical inferences on RNA-Seq data, e.g., detecting differential gene expression, are meaningful only after proper normalization. However, there is no consensus for choosing a normalization procedure from among the many existing procedures. We evaluated several RNA-Seq normalization procedures by (1) correlating estimated RNA-Seq expression values to those of microarrays, (2) examining the concordance of stable and differential gene detection between the platforms, and (3) applying the procedures to simulated RNA-Seq data. Results suggested that RNA-Seq normalization procedures have little effect on both inter-platform gene expression correlation as well as inter-platform concordance of genes detected as stably or differentially expressed. However, the results of simulated analysis suggested that some normalization procedures are more robust to changes in distribution of differentially expressed genes. These results may provide guidance for selecting RNA-Seq normalization procedures.

The effects of hypertonic fluid administration on the gene expression of inflammatory mediators in circulating leucocytes in patients with septic shock: a preliminary study

Objective

This study was designed to investigate the effect of hypertonic fluid administration on inflammatory mediator gene expression in patients with septic shock.

Design and setting

Prospective, randomized, controlled, double-blind clinical study in a 15-bed mixed intensive care unit in a tertiary referral teaching hospital.

Interventions

Twenty-four patients, who met standard criteria for septic shock, were randomized to receive a bolus of hypertonic fluid (HT, 250 ml 6% HES/7.2% NaCl) or isotonic fluid (IT, 500 ml 6% HES/0.9% NaCl) administered over 15 minutes. Randomization and study fluid administration was within 24 hours of ICU admission for all patients. This trial is registered with ANZCTR.org.au as ACTRN12607000259448.

Results

Blood samples were taken immediately before and 4, 8, 12, and 24 hours after fluid administration. Real-time reverse transcriptase polymerase chain reaction (RT rtPCR) was used to quantify mRNA expression of different inflammatory mediators in peripheral leukocytes. In the HT group, compared with the IT group, levels of gene expression of MMP9 and L-selectin were significantly suppressed (p = 0.0002 and p = 0.007, respectively), and CD11b gene expression tended to be elevated (p = NS). No differences were found in the other mediators examined.

Conclusions

In septic shock patients, hypertonic fluid administration compared with isotonic fluid may modulate expression of genes that are implicated in leukocyte-endothelial interaction and capillary leakage.

The study was performed at the Intensive Care Department, Waikato Hospital, and at the Molecular Genetics Laboratory, University of Waikato, Hamilton, New Zealand.

Trial registration

Australia and New Zealand Clinical Trials Register (ANZCTR): ACTRN12607000259448

Identification of Phosphoglycerate Kinase 1 (PGK1) as a reference gene for quantitative gene expression measurements in human blood RNA

Background

Blood is a convenient sample and increasingly used for quantitative gene expression measurements with a variety of diseases including chronic fatigue syndrome (CFS). Quantitative gene expression measurements require normalization of target genes to reference genes that are stable and independent from variables being tested in the experiment. Because there are no genes that are useful for all situations, reference gene selection is an essential step to any quantitative reverse transcription-PCR protocol. Many publications have described appropriate genes for a wide variety of tissues and experimental conditions, however, reference genes that may be suitable for the analysis of CFS, or human blood RNA derived from whole blood as well as isolated peripheral blood mononuclear cells (PBMCs), have not been described.

Findings

Literature review and analyses of our unpublished microarray data were used to narrow down the pool of candidate reference genes to six. We assayed whole blood RNA from Tempus tubes and cell preparation tube (CPT)-collected PBMC RNA from 46 subjects, and used the geNorm and NormFinder algorithms to select the most stable reference genes. Phosphoglycerate kinase 1 (PGK1) was one of the optimal normalization genes for both whole blood and PBMC RNA, however, additional genes differed for the two sample types; Ribosomal protein large, P0 (RPLP0) for PBMC RNA and Peptidylprolyl isomerase B (PPIB) for whole blood RNA. We also show that the use of a single reference gene is sufficient for normalization when the most stable candidates are used.

Conclusions

We have identified PGK1 as a stable reference gene for use with whole blood RNA and RNA derived from PBMC. When stable genes are selected it is possible to use a single gene for normalization rather than two or three. Optimal normalization will improve the ability of results from PBMC RNA to be compared with those from whole blood RNA and potentially allows comparison of gene expression results from blood RNA collected and processed by different methods with the intention of biomarker discovery. Results of this study should facilitate large-scale molecular epidemiologic studies using blood RNA as the target of quantitative gene expression measurements.

Molecular markers and mechanisms of stroke: RNA studies of blood in animals and humans

Whole genome expression microarrays can be used to study gene expression in blood, which comes in part from leukocytes, immature platelets, and red blood cells. Since these cells are important in the pathogenesis of stroke, RNA provides an index of these cellular responses to stroke. Our studies in rats have shown specific gene expression changes 24  hours after ischemic stroke, hemorrhage, status epilepticus, hypoxia, hypoglycemia, global ischemia, and following brief focal ischemia that simulated transient ischemic attacks in humans. Human studies show gene expression changes following ischemic stroke. These gene profiles predict a second cohort with >90% sensitivity and specificity. Gene profiles for ischemic stroke caused by large-vessel atherosclerosis and cardioembolism have been described that predict a second cohort with >85% sensitivity and specificity. Atherosclerotic genes were associated with clotting, platelets, and monocytes, and cardioembolic genes were associated with inflammation, infection, and neutrophils. These gene profiles predicted the cause of stroke in 58% of cryptogenic patients. These studies will provide diagnostic, prognostic, and therapeutic markers, and will advance our understanding of stroke in humans. New techniques to measure all coding and noncoding RNAs along with alternatively spliced transcripts will markedly advance molecular studies of human stroke.

Transcriptional profiling to identify biomarkers of disease and drug response

The discovery, biological qualification and analytical validation of genomic biomarkers requires extensive collaborations between individuals with expertise in biology, statistics, bioinformatics, chemistry, clinical medicine, regulatory science and so on. For clinical utility, blood-borne biomarkers (e.g., mRNA and miRNA) of organ damage, drug toxicity and/or response would be preferred to those that are tissue based. Currently used biomarkers such as serum creatinine (indicating renal dysfunction) denote organ damage whether caused by disease, physical injury or drugs. Therefore, it is anticipated that studies of disease will discover biomarkers that can also be used to identify drug-induced injury and vice versa. This article describes transcriptomic blood-borne biomarkers that have been reported to be connected with disease and drug toxicity. Much more qualification and validation needs to be carried out before many of these biomarkers can prove useful. Discussed here are some of the lessons learned and roadblocks to success.

Development of commercial gene-expression-based signatures: review of the scientific strategies

Many scientific articles have been published that use gene-expression-based technologies to discriminate a trait of interest, typically a disease subgroup, within a patient population. However, few gene-expression-based signatures have at present reached the market and become a financially and clinically successful product. The technological, scientific and medical challenges, the regulatory environment and the financial considerations are all essential parts of the development process. Here we discuss the scientific aspects of successfully developing a gene-expression-based signature and review the global strategy of six products that made it to the market. We also present a point-to-point guide that should help researchers to successfully develop genomic signatures, thus paving the way towards personalized medicine.

Evaluation of endogenous control genes for gene expression studies across multiple tissues and in the specific sets of fat- and muscle-type samples of the pig

To normalize a set of quantitative real-time PCR (q-PCR) data, it is essential to determine an optimal number/set of housekeeping genes, as the abundance of housekeeping genes can vary across tissues or cells during different developmental stages, or even under certain environmental conditions. In this study, of the 20 commonly used endogenous control genes, 13, 18 and 17 genes exhibited credible stability in 56 different tissues, 10 types of adipose tissue and five types of muscle tissue, respectively. Our analysis clearly showed that three optimal housekeeping genes are adequate for an accurate normalization, which correlated well with the theoretical optimal number (r ≥ 0.94). In terms of economical and experimental feasibility, we recommend the use of the three most stable housekeeping genes for calculating the normalization factor. Based on our results, the three most stable housekeeping genes in all analysed samples (TOP2B, HSPCB and YWHAZ) are recommended for accurate normalization of q-PCR data. We also suggest that two different sets of housekeeping genes are appropriate for 10 types of adipose tissue (the HSPCB, ALDOA and GAPDH genes) and five types of muscle tissue (the TOP2B, HSPCB and YWHAZ genes), respectively. Our report will serve as a valuable reference for other studies aimed at measuring tissue-specific mRNA abundance in porcine samples.

RefGenes: identification of reliable and condition specific reference genes for RT-qPCR data normalization

Background

RT-qPCR is a sensitive and increasingly used method for gene expression quantification. To normalize RT-qPCR measurements between samples, most laboratories use endogenous reference genes as internal controls. There is increasing evidence, however, that the expression of commonly used reference genes can vary significantly in certain contexts.

Results

Using the Genevestigator database of normalized and well-annotated microarray experiments, we describe the expression stability characteristics of the transciptomes of several organisms. The results show that a) no genes are universally stable, b) most commonly used reference genes yield very high transcript abundances as compared to the entire transcriptome, and c) for each biological context a subset of stable genes exists that has smaller variance than commonly used reference genes or genes that were selected for their stability across all conditions.

Conclusion

We therefore propose the normalization of RT-qPCR data using reference genes that are specifically chosen for the conditions under study. RefGenes is a community tool developed for that purpose. Validation RT-qPCR experiments across several organisms showed that the candidates proposed by RefGenes generally outperformed commonly used reference genes. RefGenes is available within Genevestigator at http://www.genevestigator.com.

Selection of reference genes for gene expression studies in ultraviolet B-irradiated human skin fibroblasts using quantitative real-time PCR

Background

Reference genes are frequently used to normalise mRNA levels between different samples. The expression level of these genes, however, may vary between tissues or cells and may change under certain circumstances. Cytoskeleton genes have served as multifunctional tools for experimental studies as reference genes. Our previous studies have demonstrated that the expression of vimentin, one cytoskeletal protein, was increased in ultraviolet B (UVB)-irradiated fibroblasts. Thus, we examined the expression of other cytoskeleton protein genes, ACTB (actin, beta), TUBA1A (tubulin, alpha 1a), and TUBB1 (tubulin, beta 1), in human dermal fibroblasts irradiated by UVB to determine which of these candidates were the most appropriate reference genes.

Results

Quantitative real-time PCR followed by analysis with the NormFinder and geNorm software programmes was performed. The initial screening of the expression patterns demonstrated that the expression of VIM was suppressed after UVB irradiation at doses ≥25 mJ/cm² and that the expression of TUBA1A was significantly reduced by UVB doses ≥75 mJ/cm² in cultured human dermal fibroblasts. The analysis of the experimental data revealed ACTB to be the most stably expressed gene, followed by GAPDH (aglyceraldehyde-3-phosphate dehydrogenase), under these experimental conditions. By contrast, VIM was found to be the least stable gene. The combination of ACTB and TUBB1 was revealed to be the gene pair that introduced the least systematic error into the data normalisation.

Conclusion

The data herein provide evidence that ACTB and TUBB1 are suitable reference genes in human skin fibroblasts irradiated by UVB, whereas VIM and TUBA1A are not and should therefore be excluded as reference genes in any gene expression studies involving UVB-irradiated human skin fibroblasts.

Leukocyte count affects expression of reference genes in canine whole blood samples

Background

The dog is frequently used as a model for hematologic human diseases. In this study the suitability of nine potential reference genes for quantitative RT-PCR studies in canine whole blood was investigated.

Findings

The expression of these genes was measured in whole blood samples of 263 individual dogs, representing 73 different breeds and a group of 40 mixed breed dogs, categorized into healthy dogs and dogs with internal and hematological diseases, and dogs that underwent a surgical procedure. GeNorm analysis revealed that a combination of 5 to 6 of the most stably expressed genes constituted a stable normalizing factor. Evaluation of the expression revealed different ranking of reference genes in Normfinder and GeNorm. The disease category and the white blood cell count significantly affected reference gene expression.

Conclusions

The discrepancy between the ranking of reference genes in this study by Normfinder and Genorm can be explained by differences between the experimental groups such as "disease category" and "WBC count". This stresses the importance of assessing the expression stability of potential reference genes for gene experiments in canine whole blood anew for each specific experimental condition.

mRNA-based skin identification for forensic applications

Although the identification of human skin cells is of important relevance in many forensic cases, there is currently no reliable method available. Here, we present a highly specific and sensitive messenger RNA (mRNA) approach for skin identification, meeting the key requirements in forensic analyses. We examined 11 candidate genes with skin-specific expression, as ascertained from expression databases and the literature, as well as five candidate reference genes ascertained from previous studies, in skin samples and in other forensically relevant tissues. We identified mRNA transcripts from three genes CDSN, LOR and KRT9, showing strong over-expression in skin samples relative to samples from forensic body fluids, making them suitable markers for skin identification. Out of the candidate reference genes tested, only ACTB showed similarly high expression in skin and body-fluid samples, providing a suitable reference marker for quantitative real-time PCR (qPCR) analysis of skin. Analyses of palmar and thumbprint skin samples indicate that our qPCR approach for the three skin-targeted mRNA markers, as well as the reference mRNA marker ACTB, is highly sensitive, allowing successful detection of minute amounts of skin material including full, half and quarter thumbprints, albeit with decreased success in decreasing print material. Furthermore, thumbprints stored for 6.5 months provided similar results relative to freshly analysed samples implying reasonable time-wise stability of the three skin-targeted mRNAs as well as the ACTB reference mRNA. Our study represents the first attempt towards reliable mRNA-based skin identification in forensic applications with particular relevance for future trace/touched object analyses in forensic case work. Although the approach for skin identification introduced here can be informative when applied on its own, we recommend for increased reliability the integration of (one or more of) the skin-targeted mRNA markers presented here into multiplex assays additionally including mRNA markers targeting alternative cell types expected in forensic samples.

Increased peripheral blood expression of electron transport chain genes in bipolar depression

OBJECTIVE

  To identify specific genetic pathways showing altered expression in peripheral blood of depressed subjects with bipolar disorder (BPD).

METHODS

  Illumina Sentrix BeadChip (Human-6v2) microarrays containing >48,000 transcript probes were used to measure levels of gene expression in peripheral blood from 20 depressed subjects with BPD and in 15 healthy control subjects. Quantitative real-time reverse transcription polymerase chain reaction (qRT-PCR) was used to confirm a subset of these differences.

RESULTS

  A total of 1,180 genes were differentially expressed between subjects with BPD and healthy controls (fold-change >1.3, false discovery rate-corrected p < 0.05, covaried for age and sex). Of these, 559 genes were up-regulated in BPD subjects and 621 were down-regulated. Surprisingly, there was no difference between medicated (n = 11) and unmedicated (n = 9) subjects with BPD for any of these genes. Pathway analysis using GeneGo MetaCore software showed that the most significantly affected pathway was the mitochondrial electron transport chain (ETC). Of the 85 objects (genes or proteins) in this pathway, 22 were up-regulated and 2 down-regulated in subjects with BPD. qRT-PCR confirmed up-regulation of nuclear encoded ETC genes in complexes I, III, IV, and V and, in addition, demonstrated up-regulation of mitochondrially encoded genes in each of these complexes.

CONCLUSION

  These results suggest that increased expression of multiple components of the mitochondrial ETC may be a primary deficit in bipolar depression, rather than an effect of medication.

Signatures of cardioembolic and large-vessel ischemic stroke

OBJECTIVE

The cause of stroke remains unknown or cryptogenic in many patients. We sought to determine whether gene expression signatures in blood can distinguish between cardioembolic and large-vessel causes of stroke, and whether these profiles can predict stroke etiology in the cryptogenic group.

METHODS

A total of 194 samples from 76 acute ischemic stroke patients were analyzed. RNA was isolated from blood and run on Affymetrix U133 Plus2.0 microarrays. Genes that distinguish large-vessel from cardioembolic stroke were determined at 3, 5, and 24 hours following stroke onset. Predictors were evaluated using cross-validation and a separate set of patients with known stroke subtype. The cause of cryptogenic stroke was predicted based on a model developed from strokes of known cause and identified predictors.

RESULTS

A 40-gene profile differentiated cardioembolic stroke from large-vessel stroke with >95% sensitivity and specificity. A separate 37-gene profile differentiated cardioembolic stroke due to atrial fibrillation from nonatrial fibrillation causes with >90% sensitivity and specificity. The identified genes elucidate differences in inflammation between stroke subtypes. When applied to patients with cryptogenic stroke, 17% are predicted to be large-vessel and 41% to be cardioembolic stroke. Of the cryptogenic strokes predicted to be cardioembolic, 27% were predicted to have atrial fibrillation.

INTERPRETATION

Gene expression signatures distinguish cardioembolic from large-vessel causes of ischemic stroke. These gene profiles may add valuable diagnostic information in the management of patients with stroke of unknown etiology though they need to be validated in future independent, large studies.

Gene expression profiling of blood for the prediction of ischemic stroke

BACKGROUND AND PURPOSE

A blood-based biomarker of acute ischemic stroke would be of significant value in clinical practice. This study aimed to (1) replicate in a larger cohort our previous study using gene expression profiling to predict ischemic stroke; and (2) refine prediction of ischemic stroke by including control groups relevant to ischemic stroke.

METHODS

Patients with ischemic stroke (n=70, 199 samples) were compared with control subjects who were healthy (n=38), had vascular risk factors (n=52), and who had myocardial infarction (n=17). Whole blood was drawn ≤3 hours, 5 hours, and 24 hours after stroke onset and from control subjects. RNA was processed on whole genome microarrays. Genes differentially expressed in ischemic stroke were identified and analyzed for predictive ability to discriminate stroke from control subjects.

RESULTS

The 29 probe sets previously reported predicted a new set of ischemic strokes with 93.5% sensitivity and 89.5% specificity. Sixty- and 46-probe sets differentiated control groups from 3-hour and 24-hour ischemic stroke samples, respectively. A 97-probe set correctly classified 86% of ischemic strokes (3 hour+24 hour), 84% of healthy subjects, 96% of vascular risk factor subjects, and 75% with myocardial infarction.

CONCLUSIONS

This study replicated our previously reported gene expression profile in a larger cohort and identified additional genes that discriminate ischemic stroke from relevant control groups. This multigene approach shows potential for a point-of-care test in acute ischemic stroke.

Microarray analysis of relative gene expression stability for selection of internal reference genes in the rhesus macaque brain

Background

Normalization of gene expression data refers to the comparison of expression values using reference standards that are consistent across all conditions of an experiment. In PCR studies, genes designated as "housekeeping genes" have been used as internal reference genes under the assumption that their expression is stable and independent of experimental conditions. However, verification of this assumption is rarely performed. Here we assess the use of gene microarray analysis to facilitate selection of internal reference sequences with higher expression stability across experimental conditions than can be expected using traditional selection methods.

We recently demonstrated that relative gene expression from qRT-PCR data normalized using GAPDH, ALG9 and RPL13A expression values mirrored relative expression using quantile normalization in Robust Multichip Analysis (RMA) on the Affymetrix^® GeneChip^® rhesus Macaque Genome Array.

Having shown that qRT-PCR and Affymetrix^® GeneChip^® data from the same hormone replacement therapy (HRT) study yielded concordant results, we used quantile-normalized gene microarray data to identify the most stably expressed among probe sets for prospective internal reference genes across three brain regions from the HRT study and an additional study of normally menstruating rhesus macaques (cycle study). Gene selection was limited to 575 previously published human "housekeeping" genes. Twelve animals were used per study, and three brain regions were analyzed from each animal. Gene expression stabilities were determined using geNorm, NormFinder and BestKeeper software packages.

Results

Sequences co-annotated for ribosomal protein S27a (RPS27A), and ubiquitin were among the most stably expressed under all conditions and selection criteria used for both studies. Higher annotation quality on the human GeneChip^® facilitated more targeted analysis than could be accomplished using the rhesus GeneChip^®. In the cycle study, multiple probe sets annotated for actin, gamma 1 (ACTG1) showed high signal intensity and were among the most stably expressed.

Conclusions

Using gene microarray analysis, we identified genes showing high expression stability under various sex-steroid environments in different regions of the rhesus macaque brain. Use of quantile-normalized microarray gene expression values represents an improvement over traditional methods of selecting internal reference genes for PCR analysis.