Analysis of gene expression profile during 3T3-L1 preadipocyte differentiation

Cellular differentiation is a process in which a group of differentiation specific genes is programmatically induced. This gene expression program leads to changes in both cellular morphological and physiological phenotypes. Using an 18,376-member cDNA/EST microarray, we analyzed the difference in gene expression profiles between differentiated 3T3-L1 adipocyte and non-differentiated 3T3-L1 preadipocyte. From our study a large number of genes and ESTs were identified as differentially induced or suppressed. In this paper we describe the changes of gene expression profile during 3T3-L1 cell differentiation.     

Gene Expression Analysis

The response of cells to extracellular signals usually requires altered expression of many genes, possibly including several distinct metabolic pathways. In some cases, only a subset of genes involved in such responses are known, which requires techniques to analyze changes in the expression of multiple genes, both known and unknown. Three techniques, two‐dimensional gel electrophoresis, differential display, and gene discovery arrays, provide opportunities for measuring changes in gene expression levels, as well as for identifying novel gene products.     

Detection of Gene Expression in an Individual Cell Type within a Cell Mixture Using Microarray Analysis

Background

A central issue in the design of microarray-based analysis of global gene expression is the choice between using cells of single type and a mixture of cells. This study quantified the proportion of lipopolysaccharide (LPS) induced differentially expressed monocyte genes that could be measured in peripheral blood mononuclear cells (PBMC), and determined the extent to which gene expression in the non-monocyte cell fraction diluted or obscured fold changes that could be detected in the cell mixture.

Methodology/Principal Findings

Human PBMC were stimulated with LPS, and monocytes were then isolated by positive (Mono+) or negative (Mono−) selection. The non-monocyte cell fraction (MonoD) remaining after positive selection of monocytes was used to determine the effect of non-monocyte cells on overall expression. RNA from LPS-stimulated PBMC, Mono+, Mono− and MonoD samples was co-hybridised with unstimulated RNA for each cell type on oligonucleotide microarrays. There was a positive correlation in gene expression between PBMC and both Mono+ (0.77) and Mono− (0.61–0.67) samples. Analysis of individual genes that were differentially expressed in Mono+ and Mono− samples showed that the ability to detect expression of some genes was similar when analysing PBMC, but for others, differential expression was either not detected or changed in the opposite direction. As a result of the dilutional or obscuring effect of gene expression in non-monocyte cells, overall about half of the statistically significant LPS-induced changes in gene expression in monocytes were not detected in PBMC. However, 97% of genes with a four fold or greater change in expression in monocytes after LPS stimulation, and almost all (96–100%) of the top 100 most differentially expressed monocyte genes were detected in PBMC.

Conclusions/Significance

The effect of non-responding cells in a mixture dilutes or obscures the detection of subtle changes in gene expression in an individual cell type. However, for studies in which only the most highly differentially expressed genes are of interest, separating and analysing individual cell types may be unnecessary.     

Influence of TRAIL gene on biomechanical properties of the human leukemic cell line Jurkat

We cloned the cDNA fragment of human TNF-related apoptosis inducing ligand (TRAIL) into RevTet-On, a Tet-regulated and high-level gene expression system. Making use of the TRAIL gene expression system in Jurkat as a cell model, we studied the influence of TRAIL gene on the biomechanics properties of Jurkat through measuring changes of cellular biomechanics properties before and after the TRAIL gene expression, which was induced by adding tetracycline derivative doxycycline (Dox). The results indicated that the TRAIL gene expression led to significant changes in cellular biomechanics properties. The osmotic fragility increased and the cell stiffness increased after the expression of TRAIL gene. Thus, the apoptosis-inducing TRAIL gene caused significant changes in the biomechanics properties of Jurkat cells.     

Analysis of messenger RNA expression by in situ hybridization using RNA probes synthesized via in vitro transcription

The analysis of the spatial patterning of mRNA expression is critically important for assigning functional and physiological significance to a given gene product. Given the tens of thousands of mRNAs in the mammalian genome, a full assessment of individual gene functions would ideally be overlaid upon knowledge of the specific cell types expressing each mRNA. In situ hybridization approaches represent a molecular biological/histological method that can reveal cellular patterns of mRNA expression. Here, we present detailed procedures for the detection of specific mRNAs using radioactive RNA probes in tissue sections followed by autoradiographic detection. These methods allow for the specific and sensitive detection of spatial patterns of mRNA expression, thereby linking mRNA expression with cell type and function. Radioactive detection methods also facilitate semi-quantitative analyses of changes in mRNA gene expression.     

Life Cycle Analysis of Kidney Gene Expression in Male F344 Rats

Age is a predisposing condition for susceptibility to chronic kidney disease and progression as well as acute kidney injury that may arise due to the adverse effects of some drugs. Age-related differences in kidney biology, therefore, are a key concern in understanding drug safety and disease progression. We hypothesize that the underlying suite of genes expressed in the kidney at various life cycle stages will impact susceptibility to adverse drug reactions. Therefore, establishing changes in baseline expression data between these life stages is the first and necessary step in evaluating this hypothesis. Untreated male F344 rats were sacrificed at 2, 5, 6, 8, 15, 21, 78, and 104 weeks of age. Kidneys were collected for histology and gene expression analysis. Agilent whole-genome rat microarrays were used to query global expression profiles. An ANOVA (p<0.01) coupled with a fold-change>1.5 in relative mRNA expression, was used to identify 3,724 unique differentially expressed genes (DEGs). Principal component analyses of these DEGs revealed three major divisions in life-cycle renal gene expression. K-means cluster analysis identified several groups of genes that shared age-specific patterns of expression. Pathway analysis of these gene groups revealed age-specific gene networks and functions related to renal function and aging, including extracellular matrix turnover, immune cell response, and renal tubular injury. Large age-related changes in expression were also demonstrated for the genes that code for qualified renal injury biomarkers KIM-1, Clu, and Tff3. These results suggest specific groups of genes that may underlie age-specific susceptibilities to adverse drug reactions and disease. This analysis of the basal gene expression patterns of renal genes throughout the life cycle of the rat will improve the use of current and future renal biomarkers and inform our assessments of kidney injury and disease.     

Cellular and transcriptomic analysis of NS0 cell response during exposure to hypoxia

Mammalian cells have the ability to alter their gene expression in order to survive or adapt to a variety of environment stresses including hypoxic stress. Maintaining oxygen supply has been accepted as essential for cell survival and growth. To determine the cellular and molecular changes which take place under oxygen deprivation, an NS0 cell line producing a human-mouse chimeric antibody was cultured under hypoxic conditions (<1%). Various cellular parameters such as viability, productivity, metabolism, apoptosis and cell cycle were studied and notable changes were shown to be accompanied by changes in metabolic rates. When the cells where exposed to hypoxia for 48 h, cell growth was suppressed and cell death was detected. To better understand and explore the mechanisms underpinning these biological alterations and to identify the genes involved in the genetic reprogramming, genome-wide analyses were performed using GeneChip Mouse Genome arrays. The gene expression profiling generated by the microarray technique revealed that hypoxia, even in the early stages (12h), induces significant changes in gene expression in NS0 cells. The primary responses to hypoxia within the cells were: (1) the up-regulation of pathways such as glycolysis that ultimately lead to alternative routes of ATP generation and increased oxygen availability; and (2) the down-regulation of genes involved in purine/pyrimidine and one carbon pool metabolisms required for DNA and RNA synthesis. By combining gene expression and physiological changes under hypoxia, it was possible to explore the mechanisms of hypoxia-induced alterations in more depth.     

Comparative Analysis of Hox Gene Expression in the Polychaete Chaetopterus: Implications for the Evolution of Body Plan Regionalization

The Hox genes are widely regarded as candidates for involvementin major evolutionary changes in body plan organization. Weexamine Hox gene expression data for several taxa, in relationto recent work on the polychaete annelid Chaetopterus. The workin Chaetopterus shows the basic conservation of colinearityof anterior expression boundaries seen in other groups. It alsoreveals novel patterns including early expression in the larvalgrowth zone and later formation of posterior boundaries thatcorrelate with morphological transitions in the polychaete bodyplan. The polychaete gene expression pattern is compared withthose of Hox gene homologs in other taxa to reveal differencesthat represent evolutionary changes in Hox gene regulation betweenlineages. Correlations between Hox gene expression differencesand morphological differences are examined, focussing on a numberof cases in which posterior Hox gene expression boundaries correlatewith morphological transitions. Differential regulation of theseposterior expression boundaries is proposed as a possible mechanismfor changes body plan regionalization.     

Accurate detection of aneuploidies in array CGH and gene expression microarray data

MOTIVATION: Chromosomal copy number changes (aneuploidies) are common in cell populations that undergo multiple cell divisions including yeast strains, cell lines and tumor cells. Identification of aneuploidies is critical in evolutionary studies, where changes in copy number serve an adaptive purpose, as well as in cancer studies, where amplifications and deletions of chromosomal regions have been identified as a major pathogenetic mechanism. Aneuploidies can be studied on whole-genome level using array CGH (a microarray-based method that measures the DNA content), but their presence also affects gene expression. In gene expression microarray analysis, identification of copy number changes is especially important in preventing aberrant biological conclusions based on spurious gene expression correlation or masked phenotypes that arise due to aneuploidies. Previously suggested approaches for aneuploidy detection from microarray data mostly focus on array CGH, address only whole-chromosome or whole-arm copy number changes, and rely on thresholds or other heuristics, making them unsuitable for fully automated general application to gene expression datasets. There is a need for a general and robust method for identification of aneuploidies of any size from both array CGH and gene expression microarray data. RESULTS: We present ChARM (Chromosomal Aberration Region Miner), a robust and accurate expectation-maximization based method for identification of segmental aneuploidies (partial chromosome changes) from gene expression and array CGH microarray data. Systematic evaluation of the algorithm on synthetic and biological data shows that the method is robust to noise, aneuploidal segment size and P-value cutoff. Using our approach, we identify known chromosomal changes and predict novel potential segmental aneuploidies in commonly used yeast deletion strains and in breast cancer. ChARM can be routinely used to identify aneuploidies in array CGH datasets and to screen gene expression data for aneuploidies or array biases. Our methodology is sensitive enough to detect statistically significant and biologically relevant aneuploidies even when expression or DNA content changes are subtle as in mixed populations of cells. AVAILABILITY: Code available by request from the authors and on Web supplement at http://function.cs.princeton.edu/ChARM/     

Meta-analysis of differentiating mouse embryonic stem cell gene expression kinetics reveals early change of a small gene set

Stem cell differentiation involves critical changes in gene expression. Identification of these should provide endpoints useful for optimizing stem cell propagation as well as potential clues about mechanisms governing stem cell maintenance. Here we describe the results of a new meta-analysis methodology applied to multiple gene expression datasets from three mouse embryonic stem cell (ESC) lines obtained at specific time points during the course of their differentiation into various lineages. We developed methods to identify genes with expression changes that correlated with the altered frequency of functionally defined, undifferentiated ESC in culture. In each dataset, we computed a novel statistical confidence measure for every gene which captured the certainty that a particular gene exhibited an expression pattern of interest within that dataset. This permitted a joint analysis of the datasets, despite the different experimental designs. Using a ranking scheme that favored genes exhibiting patterns of interest, we focused on the top 88 genes whose expression was consistently changed when ESC were induced to differentiate. Seven of these (103728_at, 8430410A17Rik, Klf2, Nr0b1, Sox2, Tcl1, and Zfp42) showed a rapid decrease in expression concurrent with a decrease in frequency of undifferentiated cells and remained predictive when evaluated in additional maintenance and differentiating protocols. Through a novel meta-analysis, this study identifies a small set of genes whose expression is useful for identifying changes in stem cell frequencies in cultures of mouse ESC. The methods and findings have broader applicability to understanding the regulation of self-renewal of other stem cell types.     

Is adult stem cell aging driven by conflicting modes of chromatin remodeling?

Epigenetic control of gene expression by chromatin remodeling is critical for adult stem cell function. A decline in stem cell function is observed during aging, which is accompanied by changes in the chromatin structure that are currently unexplained. Here, we hypothesize that these epigenetic changes originate from the limited cellular capability to inherit epigenetic information. We suggest that spontaneous loss of histone modification, due to fluctuations over short time scales, gives rise to long-term changes in DNA methylation and, accordingly, in gene expression. These changes are assumed to impair stem cell function and, thus, to contribute to aging. We discuss cell replication as a major source of fluctuations in histone modification patterns. Gene silencing by our proposed mechanism can be interpreted as a manifestation of the conflict between the stem cell plasticity required for tissue regeneration and the permanent silencing of potentially deleterious genomic sequences.     

Tumorigenicity of simian virus 40-hepatocyte cell lines: effect of in vitro and in vivo passage on expression of liver-specific genes and oncogenes.

Five simian virus 40 (SV40)-hepatocyte cell lines were examined for tumorigenicity and the effect of in vitro passage on the expression of four liver-specific genes (albumin, transferrin, alpha 1-antitrypsin, and phosphoenolpyruvate carboxykinase), two oncogenes (c-Ha-ras and c-raf), and two genes associated with hepatocarcinogenesis (alpha-fetoprotein and placental-type glutathione-S-transferase). At low passage (12 to 22), all five cell lines expressed the four liver-specific genes at levels similar to those in the liver and were not tumorigenic or were weakly tumorigenic. At high passage (33 to 61), the cell lines formed carcinomas, and four out of five cell lines produced primary tumors that metastasized. At least two cell lines produced well-differentiated hepatocellular carcinomas that expressed liver-specific RNAs. Levels of expression of liver-specific genes changed with time in culture. Some of the changes in liver-specific gene expression in the tumor tissue (such as for the phosphoenolpyruvate carboxykinase gene) paralleled those that occurred with in vitro passage, while other changes (such as for the albumin gene) did not parallel those that occurred with in vitro passage. Correlations between enhanced expression of c-Ha-ras and tumorigenic potential and between the process of SV40 immortalization and induced expression of c-raf and glutathione-S-transferase-P were observed. Induction of alpha-fetoprotein was detected with in vitro and in vivo passage only in the CWSV14 cell line and was paralleled by diminished albumin expression. In conclusion, we developed a model system with five SV40-hepatocyte cell lines, tumors induced by them, and tumor cell lines to examine changes in gene expression that accompany the progression from a normal cell to a hepatocellular carcinoma. Because the SV40-hepatocyte cell lines and tumor cell lines remain highly differentiated and vary in the magnitude of expression of specific genes, they can be used to study the molecular mechanisms regulating gene expression, in particular those regulating specific genes associated with differentiation.     

Influence of expressed TRAIL on biophysical properties of the human leukemic cell line Jurkat

The cDNA fragment of human TRAIL (TNF-related apoptosis inducing ligand) was cloned into RevTet-On, a Tetregulated and high-level gene expression system. The gene expression system was constructed in a human leukemic cell line: Jurkat. By using RevTet-On TRAIL gene expression system in Jurkat as a cell model, we studied the influence of TRAIL gene on the changes of cellular apoptosis before and after the TRAIL gene expression, which was induced by adding tetracycline derivative doxycycline (Dox). The results indicated that the cellular apoptosis ratio was largely dependent on the TRAIL gene expression level. Moreover, it was found that the apoptosis-inducing TRAIL could cause significant changes in the biophysical properties of Jurkat cells. The cell surface charge density decreased, the membrane fluidity declined, the elastic coefficients K1 increased, and the proportion of α-helix in membrane protein secondary structure decreased. Thus, the apoptosis-inducing TRAIL gene caused significant changes on the biomechanic properties of Jurkat cells.