Identification of novel universal housekeeping genes by statistical analysis of microarray data

Deciphering the sex bias in housekeeping gene expression in adipose tissue: a comprehensive meta-analysis of transcriptomic studies

BACKGROUND

As the housekeeping genes (HKG) generally involved in maintaining essential cell functions are typically assumed to exhibit constant expression levels across cell types, they are commonly employed as internal controls in gene expression studies. Nevertheless, HKG may vary gene expression profile according to different variables introducing systematic errors into experimental results. Sex bias can indeed affect expression display, however, up to date, sex has not been typically considered as a biological variable.

METHODS

In this study, we evaluate the expression profiles of six classical housekeeping genes (four metabolic: GAPDH, HPRT, PPIA, and UBC, and two ribosomal: 18S and RPL19) to determine expression stability in adipose tissues (AT) of Homo sapiens and Mus musculus and check sex bias and their overall suitability as internal controls. We also assess the expression stability of all genes included in distinct whole-transcriptome microarrays available from the Gene Expression Omnibus database to identify sex-unbiased housekeeping genes (suHKG) suitable for use as internal controls. We perform a novel computational strategy based on meta-analysis techniques to identify any sexual dimorphisms in mRNA expression stability in AT and to properly validate potential candidates.

RESULTS

Just above half of the considered studies informed properly about the sex of the human samples, however, not enough female mouse samples were found to be included in this analysis. We found differences in the HKG expression stability in humans between female and male samples, with females presenting greater instability. We propose a suHKG signature including experimentally validated classical HKG like PPIA and RPL19 and novel potential markers for human AT and discarding others like the extensively used 18S gene due to a sex-based variability display in adipose tissue. Orthologs have also been assayed and proposed for mouse WAT suHKG signature. All results generated during this study are readily available by accessing an open web resource ( https://bioinfo.cipf.es/metafun-HKG ) for consultation and reuse in further studies.

CONCLUSIONS

This sex-based research proves that certain classical housekeeping genes fail to function adequately as controls when analyzing human adipose tissue considering sex as a variable. We confirm RPL19 and PPIA suitability as sex-unbiased human and mouse housekeeping genes derived from sex-specific expression profiles, and propose new ones such as RPS8 and UBB.

Identification and Validation of Reference Genes for Expression Analysis in Nitrogen-Fixing Bacteria under Environmental Stress

Reference genes, also referred to as housekeeping genes (HKGs), play an important role in gene expression analysis by serving as an internal control. These HKGs are usually involved in basic cellular functions and their expression should remain at relatively constant levels. Quantitative reverse transcription-polymerase chain reaction (qRT-PCR) has been used to measure gene expression. Since the normalization of gene expression data depends on baseline expression of HKGs, it is important to identify and verify true HKGs for the qRT-PCR analysis. The goal of this study is to identify and confirm HKGs in Bradyrhizobium diazoefficiens, a nitrogen fixing bacterium which forms a symbiotic relationship with soybean. By revealing such HKGs, the normalization of gene expression would be more robust, reliable, and consistent. Here, we analyzed previous gene expression data for B. diazoefficiens under multiple environmental conditions. As a result, we identified seven constitutively expressed genes among 8453 genes across all conditions. Their fold-change values were within a range of −1.25-fold < x < 1.25-fold. We adopted GeNorm, NormFinder, and comparative ∆Ct methods to rank the seven candidate genes based on their expression stability. To validate these potential HKGs, we measured their expression in various experimental conditions, such as heat, pH, and heavy metal stress. The HKGs that were found in B. diazoefficiens were also applied in closely related species by identifying their homologs.

Conservation and architecture of housekeeping genes in the model marine diatom Thalassiosira pseudonana

Housekeeping genes (HKGs) are constitutively expressed with low variation across tissues/conditions. They are thought to be highly conserved and fundamental to cellular maintenance, with distinctive genomic features. Here, we identify 1505 HKGs in the unicellular marine diatom Thalassiosira pseudonana based on an RNA-seq analysis of 232 samples taken under 12 experimental conditions over 0-72 h. We identify promising internal reference genes (IRGs) for T. pseudonana from the most stably expressed HKGs. A comparative analysis indicates < 18% of HKGs in T. pseudonana have orthologs in other eukaryotes, including other diatom species. Contrary to work on human tissues, T. pseudonana HKGs are longer than non-HKGs, due to elongated introns. More ancient HKGs tend to be shorter than more recent HKGs, and expression levels of HKGs decrease more rapidly with gene length relative to non-HKGs. Our results indicate that HKGs are highly variable across the tree of life and thus unlikely to be universally fundamental for cellular maintenance. We hypothesize that the distinct genomic features of HKGs of T. pseudonana may be a consequence of selection pressures associated with high expression and low variance across conditions.

The role and robustness of the Gini coefficient as an unbiased tool for the selection of Gini genes for normalising expression profiling data

We recently introduced the Gini coefficient (GC) for assessing the expression variation of a particular gene in a dataset, as a means of selecting improved reference genes over the cohort ('housekeeping genes') typically used for normalisation in expression profiling studies. Those genes (transcripts) that we determined to be useable as reference genes differed greatly from previous suggestions based on hypothesis-driven approaches. A limitation of this initial study is that a single (albeit large) dataset was employed for both tissues and cell lines. We here extend this analysis to encompass seven other large datasets. Although their absolute values differ a little, the Gini values and median expression levels of the various genes are well correlated with each other between the various cell line datasets, implying that our original choice of the more ubiquitously expressed low-Gini-coefficient genes was indeed sound. In tissues, the Gini values and median expression levels of genes showed a greater variation, with the GC of genes changing with the number and types of tissues in the data sets. In all data sets, regardless of whether this was derived from tissues or cell lines, we also show that the GC is a robust measure of gene expression stability. Using the GC as a measure of expression stability we illustrate its utility to find tissue- and cell line-optimised housekeeping genes without any prior bias, that again include only a small number of previously reported housekeeping genes. We also independently confirmed this experimentally using RT-qPCR with 40 candidate GC genes in a panel of 10 cell lines. These were termed the Gini Genes. In many cases, the variation in the expression levels of classical reference genes is really quite huge (e.g. 44 fold for GAPDH in one data set), suggesting that the cure (of using them as normalising genes) may in some cases be worse than the disease (of not doing so). We recommend the present data-driven approach for the selection of reference genes by using the easy-to-calculate and robust GC.

Stimulation of cell invasion by the Golgi Ion Channel GAAP/TMBIM4 via an H2O2-Dependent Mechanism

The mechanisms by which the Golgi apparatus (GA) impacts on cell invasion are poorly understood. The human Golgi Anti-Apoptotic Protein (hGAAP, also known as TMBIM4) is a highly conserved Golgi cation channel that modulates intracellular Ca²⁺ fluxes. Human GAAP is expressed in all human tissues, is essential for cell viability and provides resistance against a range of apoptotic stresses. Furthermore, hGAAP enhances adhesion and cell migration by increasing the turnover of focal adhesions due to activation of store-operated Ca²⁺ entry. Here, we describe a GA-derived mechanism that controls cell invasion. The overexpression of hGAAP stimulates 3-dimensional proteolytic cell invasion by a mechanism that is dependent on the accumulation of intracellular hydrogen peroxide_, which might be produced by the hGAAP-dependent stimulation of mitochondrial respiration. These findings provide new insight into the complex mechanisms by which Ca²⁺ and reactive oxygen species signaling contribute to cell invasion and to the role of the GA in these processes.

Selection of best reference genes for qRT-PCR analysis of human neural stem cells preconditioned with hypoxia or baicalein-enriched fraction extracted from Oroxylum indicum medicinal plant

Whilst the potential of neural stem cell (NSC)-based treatment is recognized worldwide and seems to offer a promising therapeutic option for stroke treatments, there is currently no full understanding regarding the effects of hypoxic and baicalein-enriched fraction (BEF) preconditioning approaches on the therapeutic potential of these cells for stroke. The potential of preconditioned NSC can be determined based on the expression of several key neuroprotective genes using qRT-PCR technique. However, prior to that, it is imperative and extremely important to carefully select reference gene(s) for accurate qRT-PCR data normalization to avoid error in data interpretation. This study aimed to evaluate the stability of ten candidate reference genes via comprehensive analysis using three algorithms software: geNorm, NormFinder and BestKeeper. Our results revealed that HPRT1 and RPL13A were the most reliable reference genes for BEF-preconditioned NSCs, but ironically, HPRT1 was ranked as the least stable reference gene for hypoxic-preconditioned NSCs. On the other hand, RPLP1 and RPL13A were selected as the most stably expressed pair of reference genes for hypoxic-preconditioned NSCs. In conclusion, this study has pointed out the importance of identifying valid reference genes and has presented the first significant validation on best reference genes recommended for qRT-PCR study involves NSC preconditioned with hypoxia or with BEF extracted from Oroxylum indicum medicinal plant.

Golgi anti-apoptotic protein: a tale of camels, calcium, channels and cancer

Golgi anti-apoptotic protein (GAAP), also known as transmembrane Bax inhibitor-1 motif-containing 4 (TMBIM4) or Lifeguard 4 (Lfg4), shares remarkable amino acid conservation with orthologues throughout eukaryotes, prokaryotes and some orthopoxviruses, suggesting a highly conserved function. GAAPs regulate Ca²⁺ levels and fluxes from the Golgi and endoplasmic reticulum, confer resistance to a broad range of apoptotic stimuli, promote cell adhesion and migration via the activation of store-operated Ca²⁺ entry, are essential for the viability of human cells, and affect orthopoxvirus virulence. GAAPs are oligomeric, multi-transmembrane proteins that are resident in Golgi membranes and form cation-selective ion channels that may explain the multiple functions of these proteins. Residues contributing to the ion-conducting pore have been defined and provide the first clues about the mechanistic link between these very different functions of GAAP. Although GAAPs are naturally oligomeric, they can also function as monomers, a feature that distinguishes them from other virus-encoded ion channels that must oligomerize for function. This review summarizes the known functions of GAAPs and discusses their potential importance in disease.

GeneGini: Assessment via the Gini Coefficient of Reference "Housekeeping" Genes and Diverse Human Transporter Expression Profiles

The expression levels of SLC or ABC membrane transporter transcripts typically differ 100- to 10,000-fold between different tissues. The Gini coefficient characterizes such inequalities and here is used to describe the distribution of the expression of each transporter among different human tissues and cell lines. Many transporters exhibit extremely high Gini coefficients even for common substrates, indicating considerable specialization consistent with divergent evolution. The expression profiles of SLC transporters in different cell lines behave similarly, although Gini coefficients for ABC transporters tend to be larger in cell lines than in tissues, implying selection. Transporter genes are significantly more heterogeneously expressed than the members of most non-transporter gene classes. Transcripts with the stablest expression have a low Gini index and often differ significantly from the "housekeeping" genes commonly used for normalization in transcriptomics/qPCR studies. PCBP1 has a low Gini coefficient, is reasonably expressed, and is an excellent novel reference gene. The approach, referred to as GeneGini, provides rapid and simple characterization of expression-profile distributions and improved normalization of genome-wide expression-profiling data.

Identification of genes associated with histologic tumor grade of esophageal squamous cell carcinoma

The present study aimed to identify the genes associated with the histologic tumor grade of patients with esophageal squamous cell carcinoma (ESCC) and to provide valuable information for the identification of potential diagnostic biomarkers in ESCC. Tumor samples of ESCC patients retrieved from The Cancer Genome Atlas were divided into Grade 1 (well-differentiated; G1), Grade 2 (moderately-differentiated; G2) and Grade 3 (poorly-differentiated; G3) groups in accordance with the clinical record of the tumor grade of ESCC patients. The genes associated with tumor grade were identified. The signaling pathways of identified genes were enriched from the Kyoto Encyclopedia of Genes and Genomes (KEGG). The diagnostic value of candidate genes was assessed by receiver operating characteristic analysis. We used the GSE23400 dataset generated from the Gene Expression Omnibus to examine the expression levels of candidate genes in ESCC tissues compared to matched mucosa tissues. In total, 440 genes positively correlated with tumor grade and 882 genes negatively correlated with tumor grade were identified. There were 310 differentially expressed genes (DEGs) between G1 and G2, 184 DEGs between G2 and G3, and 710 DEGs between G1 and G3. There were 1322 genes associated with tumor grade that were significantly enriched in pathways in cancer and the phospholipase D signaling pathway. Cyclin-dependent kinase inhibitor 1A, golgin A7 family member B and transforming growth factor B1-induced anti-apoptotic factor 1 (TIAF1) had potential diagnostic value for discriminating ESCC patients with G1 from those with G3. TIAF1 was significantly down-regulated in ESCC. The results of the present study comprise useful groundwork with respect to determining the tumorigenesis mechanism in ESCC and discovering potential diagnostic biomarkers for ESCC.

Systematic identification of human housekeeping genes possibly useful as references in gene expression studies

The ideal reference, or control, gene for the study of gene expression in a given organism should be expressed at a medium-high level for easy detection, should be expressed at a constant/stable level throughout different cell types and within the same cell type undergoing different treatments, and should maintain these features through as many different tissues of the organism. From a biological point of view, these theoretical requirements of an ideal reference gene appear to be best suited to housekeeping (HK) genes. Recent advancements in the quality and completeness of human expression microarray data and in their statistical analysis may provide new clues toward the quantitative standardization of human gene expression studies in biology and medicine, both cross- and within-tissue. The systematic approach used by the present study is based on the Transcriptome Mapper tool and exploits the automated reassignment of probes to corresponding genes, intra- and inter-sample normalization, elaboration and representation of gene expression values in linear form within an indexed and searchable database with a graphical interface recording quantitative levels of expression, expression variability and cross-tissue width of expression for more than 31,000 transcripts. The present study conducted a meta-analysis of a pool of 646 expression profile data sets from 54 different human tissues and identified actin γ 1 as the HK gene that best fits the combination of all the traditional criteria to be used as a reference gene for general use; two ribosomal protein genes, RPS18 and RPS27, and one aquaporin gene, POM121 transmembrane nucleporin C, were also identified. The present study provided a list of tissue- and organ-specific genes that may be most suited for the following individual tissues/organs: Adipose tissue, bone marrow, brain, heart, kidney, liver, lung, ovary, skeletal muscle and testis; and also provides in these cases a representative, quantitative portrait of the relative, typical gene-expression profile in the form of searchable database tables.

Influence of a pre-stimulation with chronic low-dose UVB on stress response mechanisms in human skin fibroblasts

Exposure to solar ultraviolet type B (UVB), through the induction of cyclobutane pyrimidine dimer (CPD), is the major risk factor for cutaneous cancer. Cells respond to UV-induced CPD by triggering the DNA damage response (DDR) responsible for signaling DNA repair, programmed cell death and cell cycle arrest. Underlying mechanisms implicated in the DDR have been extensively studied using single acute UVB irradiation. However, little is known concerning the consequences of chronic low-dose of UVB (CLUV) on the DDR. Thus, we have investigated the effect of a CLUV pre-stimulation on the different stress response pathways. We found that CLUV pre-stimulation enhances CPD repair capacity and leads to a cell cycle delay but leave residual unrepaired CPD. We further analyzed the consequence of the CLUV regimen on general gene and protein expression. We found that CLUV treatment influences biological processes related to the response to stress at the transcriptomic and proteomic levels. This overview study represents the first demonstration that human cells respond to chronic UV irradiation by modulating their genotoxic stress response mechanisms.

Identification of reliable reference genes for qRT-PCR studies of the developing mouse mammary gland

Cell growth and differentiation are often driven by subtle changes in gene expression. Many challenges still exist in detecting these changes, particularly in the context of a complex, developing tissue. Quantitative reverse transcription polymerase chain reaction (qRT-PCR) allows relatively high-throughput evaluation of multiple genes and developmental time points. Proper quantification of gene expression levels by qRT-PCR requires normalization to one or more reference genes. Traditionally, these genes have been selected based on their presumed “housekeeping” function, with the implicit assumption that they are stably expressed over the entire experimental set. However, this is rarely tested empirically. Here we describe the identification of novel reference genes for the mouse mammary gland based on their stable expression in published microarray datasets. We compared eight novel candidate reference genes (Arpc3, Clock, Ctbp1, Phf7, Prdx1, Sugp2, Taf11 and Usp7) to eight traditional ones (18S, Actb, Gapdh, Hmbs, Hprt, Rpl13a, Sdha and Tbp) and analysed all genes for stable expression in the mouse mammary gland from pre-puberty to adulthood using four different algorithms (GeNorm, DeltaCt, BestKeeper and NormFinder). Prdx1, Phf7 and Ctbp1 were validated as novel and reliable, tissue-specific reference genes that outperform traditional reference genes in qRT-PCR studies of postnatal mammary gland development.

Identification of reference genes for quantitative real-time PCR studies in human cell lines under copper and zinc exposure

Accurate quantification depends on normalization of the measured gene expression data. In particular, gene expression studies with exposure to metals are challenging due their toxicity and redox-active properties. Here, we assessed the stability of potential reference genes in three cell lines commonly used to study metal cell metabolism: Caco-2 (colon), HepG2 (liver) and THP-1 (peripheral blood) under copper (Cu) or zinc (Zn) exposure. We used combined statistical tools to identify the best reference genes from a set of eleven candidates, which included traditional "housekeeping" genes such as GAPDH and B-ACTIN, in cell lines exposed to high and low, Zn and Cu concentrations. The expression stabilities of ATP5B (ATP synthase) and CYC1 (subunits of the cytochrome) were the highest considering the effect of Zn and Cu treatments whereas SDHA (succinate dehydrogenase) was found to be the most unstable gene. Even though the transcriptional effect of Zn and Cu is very different in term of redox properties, the same best reference genes were identified when Zn or Cu treatments were analyzed together. Our results indicate that ATP5B/CYC1 are the best candidates for reference genes after metal exposure, which can be used as a suitable starting point to evaluate gene expression with other metals or in different cell types in human models.

A strategy to identify housekeeping genes suitable for analysis in breast cancer diseases

BACKGROUND

The selection of suitable internal control genes is crucial for proper interpretation of real-time PCR data. Here we outline a strategy to identify housekeeping genes that could serve as suitable internal control for comparative analyses of gene expression data in breast cancer cell lines and tissues obtained by high throughput sequencing and quantitative real-time PCR (qRT-PCR).

METHODS

The strategy proposed includes the large-scale screening of potential candidate reference genes from RNA-seq data as well as their validation by qRT-PCR, and careful examination of reference data from the International Cancer Genome Consortium, The Cancer Genome Atlas and Gene Expression Omnibus repositories.

RESULTS

The identified set of reference genes, also called novel housekeeping genes that includes CCSER2, SYMPK, ANKRD17 and PUM1, proved to be less variable and thus potentially more accurate for research and clinical analyses of breast cell lines and tissue samples compared to the traditional housekeeping genes used to this end.

DISCUSSION

These results highlight the importance of a massive evaluation of housekeeping genes for their relevance as internal control for optimized intra- and inter-assay comparison of gene expression.

CONCLUSION

We developed a strategy to identify and evaluate the significance of housekeeping genes as internal control for the intra- and inter-assay comparison of gene expression in breast cancer that could be applied to other tumor types and diseases.

Beta-2-microglobulin (B2M) expression in the urinary sediment correlates with clinical markers of kidney disease in patients with type 1 diabetes

PURPOSE

After observing variation in the expression of the housekeeping gene B2M in cells of the urinary sediment during a study of candidate genes potentially involved in diabetic kidney disease (DKD), we hypothesized that B2M mRNA expression in the urinary sediment could reflect the presence of DKD.

METHODS

qPCR was used to quantify B2M mRNA expression in cells of the urinary sediment of 51 type 1 diabetes (T1D) patients (61% women, 33.5 [27.0-39.7] years old, with diabetes duration of 21.0 [15.0-28.0] years and HbA1c of 8.2% [7.3-8.9]; median [interquartile interval]) sorted according to the diabetic nephropathy (DN) stages; 8 focal segmental glomerulosclerosis (FSGS) patients and 10 healthy controls. B2M mRNA expression was also evaluated in human embryonic kidney epithelium-like (HEK-293) cells exposed to 25mM glucose and to albumin in order to mimic, respectively, a diabetic and a proteinuric milieu.

RESULTS

No differences were found in B2M mRNA expression among healthy controls, FSGS and T1D patients. Nonetheless B2M mRNA expression was higher in the group composed by T1D patients with incipient or overt DN combined with FSGS patients versus T1D patients without DN combined with healthy controls (P=0.0007). B2M mRNA expression was higher in T1D patients with incipient or overt DN versus without DN (P=0.03). B2M mRNA expression positively correlated with albuminuria in the overall T1D population (r=0.43; P=0.01) and negatively correlated with estimated glomerular filtration rate in male T1D patients (r=- 0.57; P=0.01). Increased B2M expression was observed in HEK-293 cells exposed to 25mM glucose and to albumin.

CONCLUSIONS

Β2M mRNA expression in cells of the urinary sediment is higher in T1D patients with DKD and in patients with FSGS in comparison to healthy subjects, maybe reflecting a tubulointerstitial injury promoted by albumin. Given the proinflammatory nature of B2M, we suggest that this protein contributes to diabetic (and possibly, to non-diabetic) tubulopathy.

Identification and Validation of Reference Genes for RT-qPCR Studies of Hypoxia in Squamous Cervical Cancer Patients

Hypoxia is an adverse factor in cervical cancer, and hypoxia-related gene expression could be a powerful biomarker for identifying the aggressive hypoxic tumors. Reverse transcription quantitative PCR (RT-qPCR) is a valuable method for gene expression studies, but suitable reference genes for data normalization that are independent of hypoxia status and clinical parameters of cervical tumors are lacking. In the present work, we aimed to identify reference genes for RT-qPCR studies of hypoxia in squamous cervical cancer. From 422 candidate reference genes selected from the literature, we used Illumina array-based expression profiles to identify 182 genes not affected by hypoxia in cervical cancer, i.e. genes regulated by hypoxia in eight cervical cancer cell lines or correlating with the hypoxia-associated dynamic contrast-enhanced magnetic resonance imaging parameter A_Brix in 42 patients, were excluded. Among the 182 genes, nine candidates (CHCHD1, GNB2L1, IPO8, LASP1, RPL27A, RPS12, SOD1, SRSF9, TMBIM6) that were not associated with tumor volume, stage, lymph node involvement or disease progression in array data of 150 patients, were selected for further testing by RT-qPCR. geNorm and NormFinder analyses of RT-qPCR data of 74 patients identified CHCHD1, SRSF9 and TMBIM6 as the optimal set of reference genes, with stable expression both overall and across patient subgroups with different hypoxia status (A_Brix) and clinical parameters. The suitability of the three reference genes were validated in studies of the hypoxia-induced genes DDIT3, ERO1A, and STC2. After normalization, the RT-qPCR data of these genes showed a significant correlation with Illumina expression (P<0.001, n = 74) and A_Brix (P<0.05, n = 32), and the STC2 data were associated with clinical outcome, in accordance with the Illumina data. Thus, CHCHD1, SRSF9 and TMBIM6 seem to be suitable reference genes for studying hypoxia-related gene expression in squamous cervical cancer samples by RT-qPCR. Moreover, STC2 is a promising prognostic hypoxia biomarker in cervical cancer.

A Simple Method for Optimization of Reference Gene Identification and Normalization in DNA Microarray Analysis

Background

Comparative DNA microarray analyses typically yield very large gene expression data sets that reflect complex patterns of change. Despite the wealth of information that is obtained, the identification of stable reference genes is required for normalization of disease- or drug-induced changes across tested groups. This is a prerequisite in quantitative real-time reverse transcription-PCR (qRT-PCR) and relative RT-PCR but rare in gene microarray analysis. The goal of the present study was to outline a simple method for identification of reliable reference genes derived from DNA microarray data sets by comparative statistical analysis of software-generated and manually calculated candidate genes.

Material/Methods

DNA microarray data sets derived from whole-blood samples obtained from 14 Zucker diabetic fatty (ZDF) rats (7 lean and 7 diabetic obese) were used for the method development. This involved the use of software-generated filtering parameters to accomplish the desired signal-to-noise ratios, 75^th percentile signal manual normalizations, and the selection of reference genes as endogenous controls for target gene expression normalization.

Results

The combination of software-generated and manual normalization methods yielded a group of 5 stably expressed, suitable endogenous control genes which can be used in further target gene expression determinations in whole blood of ZDF rats.

Conclusions

This method can be used to correct for potentially false results and aid in the selection of suitable endogenous control genes. It is especially useful when aimed to aid the software in cases of borderline results, where the expression and/or the fold change values are just beyond the pre-established set of acceptable parameters.

Consensus reference gene(s) for gene expression studies in human cancers: end of the tunnel visible?

BACKGROUND

Gene expression studies are increasingly used to provide valuable information on the diagnosis and prognosis of human cancers. Also, for in vitro and in vivo experimental cancer models gene expression studies are widely used. The complex algorithms of differential gene expression analyses require normalization of data against a reference or normalizer gene, or a set of such genes. For this purpose, mostly invariant housekeeping genes are used. Unfortunately, however, there are no consensus (housekeeping) genes that serve as reference or normalizer for different human cancers. In fact, scientists have employed a wide range of reference genes across different types of cancer for normalization of gene expression data. As a consequence, comparisons of these data and/or data harmonizations are difficult to perform and challenging. In addition, an inadequate choice for a reference gene may obscure genuine changes and/or result in erroneous gene expression data comparisons.

METHODS

In our effort to highlight the importance of selecting the most appropriate reference gene(s), we have screened the literature for gene expression studies published since the turn of the century on thirteen of the most prevalent human cancers worldwide.

CONCLUSIONS

Based on the analysis of the data at hand, we firstly recommend that in each study the suitability of candidate reference gene(s) should carefully be evaluated in order to yield reliable differential gene expression data. Secondly, we recommend that a combination of PPIA and either GAPDH, ACTB, HPRT and TBP, or appropriate combinations of two or three of these genes, should be employed in future studies, to ensure that results from different studies on different human cancers can be harmonized. This approach will ultimately increase the depth of our understanding of gene expression signatures across human cancers.

Modifications in stromal extracellular matrix of aged corneas can be induced by ultraviolet A irradiation

With age, structural and functional changes can be observed in human cornea. Some studies have shown a loss of corneal transparency and an increase in turbidity associated with aging. These changes are caused by modifications in the composition and arrangement of extracellular matrix in the corneal stroma. In human skin, it is well documented that exposure to solar radiation, and mainly to the UVA wavelengths, leads to phenotypes of photoaging characterized by alteration in extracellular matrix of the dermis. Although the cornea is also exposed to solar radiation, the extracellular matrix modifications observed in aging corneas have been mainly attributed to chronological aging and not to solar exposure. To ascertain the real implication of UVA exposure in extracellular matrix changes observed with age in human cornea, we have developed a model of photoaging by chronically exposing corneal stroma keratocytes with a precise UVA irradiation protocol. Using this model, we have analyzed UVA-induced transcriptomic and proteomic changes in corneal stroma. Our results show that cumulative UVA exposure causes changes in extracellular matrix that are found in corneal stromas of aged individuals, suggesting that solar exposure catalyzes corneal aging. Indeed, we observe a downregulation of collagen and proteoglycan gene expression and a reduction in proteoglycan production and secretion in response to cumulative UVA exposure. This study provides the first evidence that chronic ocular exposure to sunlight affects extracellular matrix composition and thus plays a role in corneal changes observed with age.

Do housekeeping genes exist?

The searching of human housekeeping (HK) genes has been a long quest since the emergence of transcriptomics, and is instrumental for us to understand the structure of genome and the fundamentals of biological processes. The resolved genes are frequently used in evolution studies and as normalization standards in quantitative gene-expression analysis. Within the past 20 years, more than a dozen HK-gene studies have been conducted, yet none of them sampled human tissues completely. We believe an integration of these results will help remove false positive genes owing to the inadequate sampling. Surprisingly, we only find one common gene across 15 examined HK-gene datasets comprising 187 different tissue and cell types. Our subsequent analyses suggest that it might not be appropriate to rigidly define HK genes as expressed in all tissue types that have diverse developmental, physiological, and pathological states. It might be beneficial to use more robustly identified HK functions for filtering criteria, in which the representing genes can be a subset of genome. These genes are not necessarily the same, and perhaps need not to be the same, everywhere in our body.

Golgi anti-apoptotic proteins are highly conserved ion channels that affect apoptosis and cell migration

Golgi anti-apoptotic proteins (GAAPs) are multitransmembrane proteins that are expressed in the Golgi apparatus and are able to homo-oligomerize. They are highly conserved throughout eukaryotes and are present in some prokaryotes and orthopoxviruses. Within eukaryotes, GAAPs regulate the Ca²⁺ content of intracellular stores, inhibit apoptosis, and promote cell adhesion and migration. Data presented here demonstrate that purified viral GAAPs (vGAAPs) and human Bax inhibitor 1 form ion channels and that vGAAP from camelpox virus is selective for cations. Mutagenesis of vGAAP, including some residues conserved in the recently solved structure of a related bacterial protein, BsYetJ, altered the conductance (E207Q and D219N) and ion selectivity (E207Q) of the channel. Mutation of residue Glu-207 or -178 reduced the effects of GAAP on cell migration and adhesion without affecting protection from apoptosis. In contrast, mutation of Asp-219 abrogated the anti-apoptotic activity of GAAP but not its effects on cell migration and adhesion. These results demonstrate that GAAPs are ion channels and define residues that contribute to the ion-conducting pore and affect apoptosis, cell adhesion, and migration independently.

Selection of proper reference genes for the cyanobacterium Synechococcus PCC 7002 using real-time quantitative PCR

Synechococcus sp. PCC 7002 is known to be tolerant to most of the environmental factors in natural habitats of Cyanobacteria. Gene expression can be easily studied in this cyanobacterium, as its complete genome sequence is available. These properties make Synechococcus sp. PCC 7002 an appropriate model organism for biotechnological applications. To study the gene expression in Cyanobacteria, real-time quantitative PCR (qPCR) can be used, but as this is a highly sensitive method, data standardization is indicated between samples. The most commonly used strategy is normalization against internal reference genes. Synechococcus sp. PCC 7002 has not yet been evaluated for the best reference genes. In this work, six candidate genes were analyzed for this purpose. Cyanobacterial cultures were exposed to several stress conditions, and three different algorithms were used for ranking the reference genes: geNorm, NormFinder, and BestKeeper. Moreover, gene expression stability value M and single-control normalization error E were calculated. Our data provided a list of reference genes that can be used in qPCR experiments in Synechococcus sp. PCC 7002.

Partitioning the human transcriptome using HKera, a novel classifier of housekeeping and tissue-specific genes

High-throughput transcriptomic experiments have made it possible to classify genes that are ubiquitously expressed as housekeeping (HK) genes and those expressed only in selective tissues as tissue-specific (TS) genes. Although partitioning a transcriptome into HK and TS genes is conceptually problematic owing to the lack of precise definitions and gene expression profile criteria for the two, information whether a gene is an HK or a TS gene can provide an initial clue to its cellular and/or functional role. Consequently, the development of new and novel HK (TS) classification methods has been a topic of considerable interest in post-genomics research. Here, we report such a development. Our method, called HKera, differs from the others by utilizing a novel property of HK genes that we have previously uncovered, namely that the ranking order of their expression levels, as opposed to the expression levels themselves, tends to be preserved from one tissue to another. Evaluated against multiple benchmark sets of human HK genes, including one recently derived from second generation sequencing data, HKera was shown to perform significantly better than five other classifiers that use different methodologies. An enrichment analysis of pathway and gene ontology annotations showed that HKera-predicted HK and TS genes have distinct functional roles and, together, cover most of the ontology categories. These results show that HKera is a good transcriptome partitioner that can be used to search for, and obtain useful expression and functional information for, novel HK (TS) genes.

PaGenBase: a pattern gene database for the global and dynamic understanding of gene function

Pattern genes are a group of genes that have a modularized expression behavior under serial physiological conditions. The identification of pattern genes will provide a path toward a global and dynamic understanding of gene functions and their roles in particular biological processes or events, such as development and pathogenesis. In this study, we present PaGenBase, a novel repository for the collection of tissue- and time-specific pattern genes, including specific genes, selective genes, housekeeping genes and repressed genes. The PaGenBase database is now freely accessible at http://bioinf.xmu.edu.cn/PaGenBase/. In the current version (PaGenBase 1.0), the database contains 906,599 pattern genes derived from the literature or from data mining of more than 1,145,277 gene expression profiles in 1,062 distinct samples collected from 11 model organisms. Four statistical parameters were used to quantitatively evaluate the pattern genes. Moreover, three methods (quick search, advanced search and browse) were designed for rapid and customized data retrieval. The potential applications of PaGenBase are also briefly described. In summary, PaGenBase will serve as a resource for the global and dynamic understanding of gene function and will facilitate high-level investigations in a variety of fields, including the study of development, pathogenesis and novel drug discovery.

hGAAP promotes cell adhesion and migration via the stimulation of store-operated Ca2+ entry and calpain 2

Golgi antiapoptotic proteins (GAAPs) are highly conserved Golgi membrane proteins that inhibit apoptosis and promote Ca(2+) release from intracellular stores. Given the role of Ca(2+) in controlling cell adhesion and motility, we hypothesized that human GAAP (hGAAP) might influence these events. In this paper, we present evidence that hGAAP increased cell adhesion, spreading, and migration in a manner that depended on the C-terminal domain of hGAAP. We show that hGAAP increased store-operated Ca(2+) entry and thereby the activity of calpain at newly forming protrusions. These hGAAP-dependent effects regulated focal adhesion dynamics and cell migration. Indeed, inhibition or knockdown of calpain 2 abrogated the effects of hGAAP on cell spreading and migration. Our data reveal that hGAAP is a novel regulator of focal adhesion dynamics, cell adhesion, and migration by controlling localized Ca(2+)-dependent activation of calpain.

Human housekeeping genes, revisited

Housekeeping genes are involved in basic cell maintenance and, therefore, are expected to maintain constant expression levels in all cells and conditions. Identification of these genes facilitates exposure of the underlying cellular infrastructure and increases understanding of various structural genomic features. In addition, housekeeping genes are instrumental for calibration in many biotechnological applications and genomic studies. Advances in our ability to measure RNA expression have resulted in a gradual increase in the number of identified housekeeping genes. Here, we describe housekeeping gene detection in the era of massive parallel sequencing and RNA-seq. We emphasize the importance of expression at a constant level and provide a list of 3804 human genes that are expressed uniformly across a panel of tissues. Several exceptionally uniform genes are singled out for future experimental use, such as RT-PCR control genes. Finally, we discuss both ways in which current technology can meet some of past obstacles encountered, and several as yet unmet challenges.

The role of vascular endothelial growth factor and matrix metalloproteinases in canine lymphoma: in vivo and in vitro study

BACKGROUND

Canine lymphoma represents the most frequent haematopoietic cancer and it shares some similarities with human non-Hodgkin lymphoma. Matrix metalloproteinases (MMPs) and vascular endothelial growth factor (VEGF) play a coordinated role during invasion and proliferation of malignant cells; however, little is known about their role in canine haematologic malignancies. The aim of this study was to investigate the mRNA and protein expression of VEGF and the most relevant MMPs in canine lymphoma. Lymph node aspirates from 26 B-cell and 21 T-cell lymphomas were collected. The protein expression levels of MMP-9, MMP-2 and VEGF-A were evaluated by immunocytochemistry, and the mRNA levels of MMP-2, MMP-9, MT1-MMP, TIMP-1, TIMP-2, RECK, VEGF-A and VEGF-164 were measured using quantitative RT-PCR.

RESULTS

MT1-MMP, TIMP-1 and RECK mRNA levels were significantly higher in T-cell lymphomas than in B-cell lymphomas. Higher mRNA and protein levels of MMP-9 and VEGF-A were observed in T-cell lymphomas than in B-cell lymphomas and healthy control lymph nodes. A positive correlation was found between MMP-9 and VEGF-A in T-cell lymphomas. Moreover, MMP-9, MT1-MMP, TIMP-1 and VEGF-A were expressed at the highest levels in high-grade T-cell lymphomas.

CONCLUSIONS

This study provides new information on the expression of different MMPs and VEGF in canine lymphoma, suggesting a possible correlation between different MMPs and VEGF, immunophenotype and prognosis.

Human and viral Golgi anti-apoptotic proteins (GAAPs) oligomerize via different mechanisms and monomeric GAAP inhibits apoptosis and modulates calcium

Golgi anti-apoptotic proteins (GAAPs) are hydrophobic proteins resident in membranes of the Golgi complex. They protect cells from a range of apoptotic stimuli, reduce the Ca²⁺ content of intracellular stores, and regulate Ca²⁺ fluxes. GAAP was discovered in camelpox virus, but it is highly conserved throughout evolution and encoded by all eukaryote genomes examined. GAAPs are part of the transmembrane Bax inhibitor-containing motif (TMBIM) family that also includes other anti-apoptotic and Ca²⁺-modulating membrane proteins. Most TMBIM members show multiple bands when analyzed by SDS-PAGE, suggesting that they may be oligomeric. However, the molecular mechanisms of oligomerization, the native state of GAAPs in living cells and the functional significance of oligomerization have not been addressed. TMBIM members are thought to have evolved from an ancestral GAAP. Two different GAAPs, human (h) and viral (v)GAAP were therefore selected as models to examine oligomerization of TMBIM family members. We show that both hGAAP and vGAAP in their native states form oligomers and that oligomerization is pH-dependent. Surprisingly, hGAAP and vGAAP do not share the same oligomerization mechanism. Oligomerization of hGAAP is independent of cysteines, but oligomerization of vGAAP depends on cysteines 9 and 60. A mutant vGAAP that is unable to oligomerize revealed that monomeric vGAAP retains both its anti-apoptotic function and its effect on intracellular Ca²⁺ stores. In conclusion, GAAP can oligomerize in a pH-regulated manner, and monomeric GAAP is functional.

Transcriptome Analysis of Spermophilus lateralis and Spermophilus tridecemlineatus Liver Does Not Suggest the Presence of Spermophilus-Liver-Specific Reference Genes

The expressions of reference genes used in gene expression studies are assumed to be stable under most circumstances. However, studies had demonstrated that genes assumed to be stably expressed in a species are not necessarily stably expressed in other organisms. This study aims to evaluate the likelihood of genus-specific reference genes for liver using comparable microarray datasets from Spermophilus lateralis and Spermophilus tridecemlineatus. The coefficient of variance (CV) of each probe was calculated and there were 178 probes common between the lowest 10% CV of both datasets (n = 1258). All 3 lists were analysed by NormFinder. Our results suggest that the most invariant probe for S. tridecemlineatus was 02n12, while that for S. lateralis was 24j21. However, our results showed that Probes 02n12 and 24j21 are ranked 8644 and 926 in terms of invariancy for S. lateralis and S. tridecemlineatus respectively. This suggests the lack of common liver-specific reference probes for both S. lateralis and S. tridecemlineatus. Given that S. lateralis and S. tridecemlineatus are closely related species and the datasets are comparable, our results do not support the presence of genus-specific reference genes.

Selection of novel reference genes for use in the human central nervous system: a BrainNet Europe Study

The use of an appropriate reference gene to ensure accurate normalisation is crucial for the correct quantification of gene expression using qPCR assays and RNA arrays. The main criterion for a gene to qualify as a reference gene is a stable expression across various cell types and experimental settings. Several reference genes are commonly in use but more and more evidence reveals variations in their expression due to the presence of on-going neuropathological disease processes, raising doubts concerning their use. We conducted an analysis of genome-wide changes of gene expression in the human central nervous system (CNS) covering several neurological disorders and regions, including the spinal cord, and were able to identify a number of novel stable reference genes. We tested the stability of expression of eight novel (ATP5E, AARS, GAPVD1, CSNK2B, XPNPEP1, OSBP, NAT5 and DCTN2) and four more commonly used (BECN1, GAPDH, QARS and TUBB) reference genes in a smaller cohort using RT-qPCR. The most stable genes out of the 12 reference genes were tested as normaliser to validate increased levels of a target gene in CNS disease. We found that in human post-mortem tissue the novel reference genes, XPNPEP1 and AARS, were efficient in replicating microarray target gene expression levels and that XPNPEP1 was more efficient as a normaliser than BECN1, which has been shown to change in expression as a consequence of neuronal cell loss. We provide herein one more suitable novel reference gene, XPNPEP1, with no current neuroinflammatory or neurodegenerative associations that can be used for gene quantitative gene expression studies with human CNS post-mortem tissue and also suggest a list of potential other candidates. These data also emphasise the importance of organ/tissue-specific stably expressed genes as reference genes for RNA studies.

Expression of the aryl hydrocarbon receptor pathway and cyclooxygenase-2 in dog tumors

In humans, the aryl hydrocarbon receptor (AHR) gene battery constitutes a set of contaminant-responsive genes, which have been recently shown to be involved in the regulation of several patho-physiological conditions, including tumorigenesis. As the domestic dog represents a valuable animal model in comparative oncology, mRNA levels of cytochromes P450 1A1, 1A2 and 1B1 (CYP1A1, 1A2 and 1B1), AHR, AHR nuclear translocator (ARNT), AHR repressor (AHRR, whose partial sequence was here obtained) and cyclooxygenase-2 (COX2) were measured in dog control tissues (liver, skin, mammary gland and bone), in 47 mast cell tumors (MCTs), 32 mammary tumors (MTs), 5 osteosarcoma (OSA) and related surgical margins. Target genes were constitutively expressed in the dog, confirming the available human data. Furthermore, their pattern of expression in tumor biopsies was comparable to that already described in a variety of human cancers; in particular, both AHR and COX2 genes were up-regulated and positively correlated, while CYP1A1 and CYP1A2 mRNAs were generally poorly expressed. This work demonstrated for the first time that target mRNAs are expressed in neoplastic tissues of dogs, thereby increasing the knowledge about dog cancer biology and confirming this species as an useful animal model for comparative studies on human oncology.

Six-transmembrane topology for Golgi anti-apoptotic protein (GAAP) and Bax inhibitor 1 (BI-1) provides model for the transmembrane Bax inhibitor-containing motif (TMBIM) family

The Golgi anti-apoptotic protein (GAAP) is a hydrophobic Golgi protein that regulates intracellular calcium fluxes and apoptosis. GAAP is highly conserved throughout eukaryotes and some strains of vaccinia virus (VACV) and camelpox virus. Based on sequence, phylogeny, and hydrophobicity, GAAPs were classified within the transmembrane Bax inhibitor-containing motif (TMBIM) family. TMBIM members are anti-apoptotic and were predicted to have seven-transmembrane domains (TMDs). However, topology prediction programs are inconsistent and predicted that GAAP and other TMBIM members have six or seven TMDs. To address this discrepancy, we mapped the transmembrane topology of viral (vGAAP) and human (hGAAP), as well as Bax inhibitor (BI-1). Data presented show a six-, not seven-, transmembrane topology for vGAAP with a putative reentrant loop at the C terminus and both termini located in the cytosol. We find that this topology is also conserved in hGAAP and BI-1. This places the charged C terminus in the cytosol, and mutation of these charged residues in hGAAP ablated its anti-apoptotic function. Given the highly conserved hydrophobicity profile within the TMBIM family and recent phylogenetic data indicating that a GAAP-like protein may have been the ancestral progenitor of a subset of the TMBIM family, we propose that this vGAAP topology may be used as a model for the remainder of the TMBIM family of proteins. The topology described provides valuable information on the structure and function of an important but poorly understood family of proteins.

Matrix metalloproteinases and their inhibitors in canine mammary tumors

Background

Malignant canine mammary tumors represent 50% of all neoplasms in female dogs. Matrix metalloproteinases (MMPs) and tissue inhibitors of metalloproteinases (TIMPs) are thought to be involved in tumor progression, and they are also associated with the reactive stroma, which provides structural and vascular support for tumor growth.

Results

MMP-2, MMP-9 and MT1-MMP were expressed at both the mRNA and protein levels in tumor samples. MMP-2 and MMP-9 immunohistochemical reactions were evident both in the epithelial tumor cells and in the stromal compartment to varying degrees; in particular, the intensity of the MMP-2 staining was stronger in the stromal fibroblasts close to epithelial tumor cells in simple carcinomas than in adenomas. These data were supported by gelatin-zymography; bands for the active form of MMP-2 were found in 94% of carcinoma samples, compared with 17% of benign tumor samples. The gene expression and immunohistochemical results for MT1-MMP were comparable to those for MMP-2. The immunoreactivity for MMP-13 and TIMP-2 was lower in carcinomas than in adenomas, confirming the mRNA data for MMP-13 and the other MMP inhibitors that were evaluated. The active form of MMP-9, but not the active form of MMP-2, was identified in the plasma of all of the tested dogs.

Conclusions

Our findings suggest that MMP-9, MMP-2 and MT1-MMP, which are synthesized by epithelial cancer cells and cancer-associated fibroblasts, play an important role in malignant canine mammary tumors. The reduction of MMP-13 and TIMP-2 could also be a significant step in malignant transformation. MMP-2 and MT1-MMP could be further evaluated as future biomarkers for predicting the progression and prognosis of canine mammary tumors.

Ensuring the safety of vaccine cell substrates by massively parallel sequencing of the transcriptome

Massively parallel, deep, sequencing of the transcriptome coupled with algorithmic analysis to identify adventitious agents (MP-Seq™) is an important adjunct in ensuring the safety of cells used in vaccine production. Such cells may harbour novel viruses whose sequences are unknown or latent viruses that are only expressed following stress to the cells. MP-Seq is an unbiased and comprehensive method to identify such viruses and other adventitious agents without prior knowledge of the nature of those agents. Here we demonstrate its utility as part of an integrated approach to identify and characterise potential contaminants within commonly used virus and vaccine production cell lines. Through this analysis, in combination with more traditional approaches, we have excluded the presence of porcine circoviruses in the ATCC Vero cell bank (CCL-81), however, we found that a full length betaretrovirus related to SRV can be expressed in these cells, a factor that may be of importance in the production of certain vaccines. Similarly, insect cells are proving to be valuable for the production of virus like particles and sub-unit vaccines, but they can harbour a range of latent viruses. We show that following MP-Seq of the Trichoplusia ni (High Five cell line) transcriptome we were able to detect a contaminating, latent nodavirus and identify an expressed errantivirus genome. Collectively, these studies have reinforced the role of MP-Seq as an integral tool for the identification of contaminating agents in vaccine cell substrates.

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.

Identification of reference genes across physiological states for qRT-PCR through microarray meta-analysis

Background

The accuracy of quantitative real-time PCR (qRT-PCR) is highly dependent on reliable reference gene(s). Some housekeeping genes which are commonly used for normalization are widely recognized as inappropriate in many experimental conditions. This study aimed to identify reference genes for clinical studies through microarray meta-analysis of human clinical samples.

Methodology/Principal Findings

After uniform data preprocessing and data quality control, 4,804 Affymetrix HU-133A arrays performed by clinical samples were classified into four physiological states with 13 organ/tissue types. We identified a list of reference genes for each organ/tissue types which exhibited stable expression across physiological states. Furthermore, 102 genes identified as reference gene candidates in multiple organ/tissue types were selected for further analysis. These genes have been frequently identified as housekeeping genes in previous studies, and approximately 71% of them fall into Gene Expression (GO:0010467) category in Gene Ontology.

Conclusions/Significance

Based on microarray meta-analysis of human clinical sample arrays, we identified sets of reference gene candidates for various organ/tissue types and then examined the functions of these genes. Additionally, we found that many of the reference genes are functionally related to transcription, RNA processing and translation. According to our results, researchers could select single or multiple reference gene(s) for normalization of qRT-PCR in clinical studies.

Measuring 'normalcy' in plant gene expression after herbivore attack

Plants make drastic changes to their transcriptome to appropriately respond to environmental change, and the regulation of genes that are specific to abiotic and biotic stresses is a key to plant survival. The coordination of defence gene transcription is often coupled with significant adjustments in the levels of expression of primary metabolic and structural genes to relocate resources, repair damage and/or induce senescence. This complicates the process of finding suitable 'housekeeping' or reference genes to use in measurements of gene expression by real-time reverse transcription (RT-PCR) in response to herbivore attack. Several software programs have been developed to identify candidate reference genes, but measurement of their expression may still not yield an appropriate gene or suite of genes for normalization. This is especially true in plant-herbivore interactions where tissue damage is immediate and continuous. Here, we show that 12 traditional reference genes customarily used in RT-PCR analysis are not stably expressed after insect attack. We describe the pitfalls of using traditional reference genes and why insect attack may be affecting whole cell metabolism. We propose a method using RNA quantification in combination with an external spike of commercially available mRNA as normalization factors in studies involving herbivory, multiple stress treatments or species where stable reference genes are unknown.

Microarray-driven validation of reference genes for quantitative real-time polymerase chain reaction in a rat vocal fold model of mucosal injury

Relative quantification by normalization against a stably expressed reference gene is a widely used data analysis method in microarray and quantitative real-time polymerase chain reaction (qRT-PCR) platforms; however, recent evidence suggests that many commonly utilized reference genes are unstable in certain experimental systems and situations. The primary aim of this study, therefore, was to screen and identify stably expressed reference genes in a well-established rat model of vocal fold mucosal injury. We selected and evaluated the expression stability of nine candidate reference genes. Ablim1, Sptbn1, and Wrnip1 were identified as stably expressed in a model-specific microarray dataset and were further validated as suitable reference genes in an independent qRT-PCR experiment using 2(-DeltaCT) and pairwise comparison-based (geNorm) analyses. Parallel analysis of six commonly used reference genes identified Sdha as the only stably expressed candidate in this group. Sdha, Sptbn1, and the geometric mean of Sdha and Sptbn1 each provided accurate normalization of target gene Tgfb1; Gapdh, the least stable candidate gene in our dataset, provided inaccurate normalization and an invalid experimental result. The stable reference genes identified here are suitable for accurate normalization of target gene expression in vocal fold mucosal injury experiments.

Tissue banking of diagnostic lung cancer biopsies for extraction of high quality RNA

INTRODUCTION

There is a clear need to develop a practical approach to obtain high quality RNA for gene expression analysis from lung cancer patients. Current approaches are restricted to using material from surgical resection specimens. We systematically investigated whether high quality RNA could be obtained from routine lung cancer diagnostic biopsies to determine the optimum method.

METHODS

Extra biopsies were taken at diagnosis from patients later confirmed to have lung cancer. Comparisons were made between RNA extracted from samples snap frozen in liquid nitrogen and those treated with an RNA preservative before freezing. Further comparisons were made between biopsies taken by different methods.

RESULTS

Acceptable RNA for gene expression analysis was extracted from 72% of lung cancer biopsies. Use of an RNA preservative for storage allowed the extraction of higher quality, more intact RNA from biopsies gathered by both endobronchial forceps and transbronchial needle aspiration. High quality RNA could also be extracted from computed tomography-guided needle core biopsies.

CONCLUSION

Banking lung cancer biopsy specimens by storage in an RNA preservative solution will allow use of a broader spectrum of lung cancers for gene expression analysis. We describe a model that makes personalized medicine for lung cancer patients a more practical proposition.

A sequence-based approach to identify reference genes for gene expression analysis

BACKGROUND

An important consideration when analyzing both microarray and quantitative PCR expression data is the selection of appropriate genes as endogenous controls or reference genes. This step is especially critical when identifying genes differentially expressed between datasets. Moreover, reference genes suitable in one context (e.g. lung cancer) may not be suitable in another (e.g. breast cancer). Currently, the main approach to identify reference genes involves the mining of expression microarray data for highly expressed and relatively constant transcripts across a sample set. A caveat here is the requirement for transcript normalization prior to analysis, and measurements obtained are relative, not absolute. Alternatively, as sequencing-based technologies provide digital quantitative output, absolute quantification ensues, and reference gene identification becomes more accurate.

METHODS

Serial analysis of gene expression (SAGE) profiles of non-malignant and malignant lung samples were compared using a permutation test to identify the most stably expressed genes across all samples. Subsequently, the specificity of the reference genes was evaluated across multiple tissue types, their constancy of expression was assessed using quantitative RT-PCR (qPCR), and their impact on differential expression analysis of microarray data was evaluated.

RESULTS

We show that (i) conventional references genes such as ACTB and GAPDH are highly variable between cancerous and non-cancerous samples, (ii) reference genes identified for lung cancer do not perform well for other cancer types (breast and brain), (iii) reference genes identified through SAGE show low variability using qPCR in a different cohort of samples, and (iv) normalization of a lung cancer gene expression microarray dataset with or without our reference genes, yields different results for differential gene expression and subsequent analyses. Specifically, key established pathways in lung cancer exhibit higher statistical significance using a dataset normalized with our reference genes relative to normalization without using our reference genes.

CONCLUSIONS

Our analyses found NDUFA1, RPL19, RAB5C, and RPS18 to occupy the top ranking positions among 15 suitable reference genes optimal for normalization of lung tissue expression data. Significantly, the approach used in this study can be applied to data generated using new generation sequencing platforms for the identification of reference genes optimal within diverse contexts.

Gene expression results in lipopolysaccharide-stimulated monocytes depend significantly on the choice of reference genes

Background

Gene expression in lipopolysaccharide (LPS)-stimulated monocytes is mainly studied by quantitative real-time reverse transcription PCR (RT-qPCR) using GAPDH (glyceraldehyde 3-phosphate dehydrogenase) or ACTB (beta-actin) as reference gene for normalization. Expression of traditional reference genes has been shown to vary substantially under certain conditions leading to invalid results. To investigate whether traditional reference genes are stably expressed in LPS-stimulated monocytes or if RT-qPCR results are dependent on the choice of reference genes, we have assessed and evaluated gene expression stability of twelve candidate reference genes in this model system.

Results

Twelve candidate reference genes were quantified by RT-qPCR in LPS-stimulated, human monocytes and evaluated using the programs geNorm, Normfinder and BestKeeper. geNorm ranked PPIB (cyclophilin B), B2M (beta-2-microglobulin) and PPIA (cyclophilin A) as the best combination for gene expression normalization in LPS-stimulated monocytes. Normfinder suggested TBP (TATA-box binding protein) and B2M as the best combination. Compared to these combinations, normalization using GAPDH alone resulted in significantly higher changes of TNF-α (tumor necrosis factor-alpha) and IL10 (interleukin 10) expression. Moreover, a significant difference in TNF-α expression between monocytes stimulated with equimolar concentrations of LPS from N. meningitides and E. coli, respectively, was identified when using the suggested combinations of reference genes for normalization, but stayed unrecognized when employing a single reference gene, ACTB or GAPDH.

Conclusions

Gene expression levels in LPS-stimulated monocytes based on RT-qPCR results differ significantly when normalized to a single gene or a combination of stably expressed reference genes. Proper evaluation of reference gene stabiliy is therefore mandatory before reporting RT-qPCR results in LPS-stimulated monocytes.

Normalization with genes encoding ribosomal proteins but not GAPDH provides an accurate quantification of gene expressions in neuronal differentiation of PC12 cells

Background

Gene regulation at transcript level can provide a good indication of the complex signaling mechanisms underlying physiological and pathological processes. Transcriptomic methods such as microarray and quantitative real-time PCR require stable reference genes for accurate normalization of gene expression. Some but not all studies have shown that housekeeping genes (HGKs), β-actin (ACTB) and glyceraldehyde-3-phosphate dehydrogenase (GAPDH), which are routinely used for normalization, may vary significantly depending on the cell/tissue type and experimental conditions. It is currently unclear if these genes are stably expressed in cells undergoing drastic morphological changes during neuronal differentiation. Recent meta-analysis of microarray datasets showed that some but not all of the ribosomal protein genes are stably expressed. To test the hypothesis that some ribosomal protein genes can serve as reference genes for neuronal differentiation, a genome-wide analysis was performed and putative reference genes were identified based on stability of expressions. The stabilities of these potential reference genes were then analyzed by reverse transcription quantitative real-time PCR in six differentiation conditions.

Results

Twenty stably expressed genes, including thirteen ribosomal protein genes, were selected from microarray analysis of the gene expression profiles of GDNF and NGF induced differentiation of PC12 cells. The expression levels of these candidate genes as well as ACTB and GAPDH were further analyzed by reverse transcription quantitative real-time PCR in PC12 cells differentiated with a variety of stimuli including NGF, GDNF, Forskolin, KCl and ROCK inhibitor, Y27632. The performances of these candidate genes as stable reference genes were evaluated with two independent statistical approaches, geNorm and NormFinder.

Conclusions

The ribosomal protein genes, RPL19 and RPL29, were identified as suitable reference genes during neuronal differentiation of PC12 cells, regardless of the type of differentiation conditions. The combination of these two novel reference genes, but not the commonly used HKG, GAPDH, allows robust and accurate normalization of differentially expressed genes during PC12 differentiation.

Using microarray technology to select housekeeping genes in Chinese hamster ovary cells

In the present study, we have identified species-specific housekeeping genes (HKGs) for Chinese Hamster Ovary (CHO) cells using data from microarray gene expression profiling. HKGs suitable for quantitative RT-PCR normalization should display relatively stable expression levels across experimental conditions. We analyzed transcription profiles of several IgG-producing recombinant CHO cell lines under numerous culture conditions using a custom CHO DNA microarray platform and calculated relative expression variability from 124 microarrays. We selected a novel panel of candidate HKGs based on their low variability in expression from the microarray data. Compared to three traditional HKGs (Gapdh, Actb, and B2m), the majority of genes on this panel demonstrated lower or equal variability. Each candidate HKG was then validated using qRT-PCR. Final selection of CHO-specific HKGs include Actr5, Eif3i, Hirip3, Pabpn1, Vezt, Cog1, and Yaf2. The results reported here provide a useful tool for gene expression analyses in CHO cells, a critical expression platform used in biotherapeutics.

Dynamism in gene expression across multiple studies

In this study we develop methods of examining gene expression dynamics, how and when genes change expression, and demonstrate their application in a meta-analysis involving over 29,000 microarrays. By defining measures across many experimental conditions, we have a new way of characterizing dynamics, complementary to measures looking at changes in absolute variation or breadth of tissues showing expression. We show conservation in overall patterns of dynamism across three species (human, mouse, and rat) and show associations with known disease-related genes. We discuss the enriched functional properties of the sets of genes showing different patterns of dynamics and show that the differences in expression dynamics is associated with the variety of different transcription factor regulatory sites. These results can influence thinking about the selection of genes for microarray design and the analysis of measurements of mRNA expression variation in a global context of expression dynamics across many conditions, as genes that are rarely differentially expressed between experimental conditions may be the subject of increased scrutiny when they significantly vary in expression between experimental subsets.

Validation of reference genes for quantitative expression analysis by real-time RT-PCR in Saccharomyces cerevisiae

Background

Real-time RT-PCR is the recommended method for quantitative gene expression analysis. A compulsory step is the selection of good reference genes for normalization. A few genes often referred to as HouseKeeping Genes (HSK), such as ACT1, RDN18 or PDA1 are among the most commonly used, as their expression is assumed to remain unchanged over a wide range of conditions. Since this assumption is very unlikely, a geometric averaging of multiple, carefully selected internal control genes is now strongly recommended for normalization to avoid this problem of expression variation of single reference genes. The aim of this work was to search for a set of reference genes for reliable gene expression analysis in Saccharomyces cerevisiae.

Results

From public microarray datasets, we selected potential reference genes whose expression remained apparently invariable during long-term growth on glucose. Using the algorithm geNorm, ALG9, TAF10, TFC1 and UBC6 turned out to be genes whose expression remained stable, independent of the growth conditions and the strain backgrounds tested in this study. We then showed that the geometric averaging of any subset of three genes among the six most stable genes resulted in very similar normalized data, which contrasted with inconsistent results among various biological samples when the normalization was performed with ACT1. Normalization with multiple selected genes was therefore applied to transcriptional analysis of genes involved in glycogen metabolism. We determined an induction ratio of 100-fold for GPH1 and 20-fold for GSY2 between the exponential phase and the diauxic shift on glucose. There was no induction of these two genes at this transition phase on galactose, although in both cases, the kinetics of glycogen accumulation was similar. In contrast, SGA1 expression was independent of the carbon source and increased by 3-fold in stationary phase.

Conclusion

In this work, we provided a set of genes that are suitable reference genes for quantitative gene expression analysis by real-time RT-PCR in yeast biological samples covering a large panel of physiological states. In contrast, we invalidated and discourage the use of ACT1 as well as other commonly used reference genes (PDA1, TDH3, RDN18, etc) as internal controls for quantitative gene expression analysis in yeast.