Design, construction, and analysis of cell line arrays and tissue microarrays for gene expression analysis

Cell line array (CMA) and tissue microarray (TMA) technologies are high-throughput methods for analysing both the abundance and distribution of gene expression in a panel of cell lines or multiple tissue specimens in an efficient and cost-effective manner. The process is based on Kononen's method of extracting a cylindrical core of paraffin-embedded donor tissue and inserting it into a recipient paraffin block. Donor tissue from surgically resected paraffin-embedded tissue blocks, frozen needle biopsies or cell line pellets can all be arrayed in the recipient block. The representative area of interest is identified and circled on a haematoxylin and eosin (H&E)-stained section of the donor block. Using a predesigned map showing a precise spacing pattern, a high density array of up to 1,000 cores of cell pellets and/or donor tissue can be embedded into the recipient block using a tissue arrayer from Beecher Instruments. Depending on the depth of the cell line/tissue removed from the donor block 100-300 consecutive sections can be cut from each CMA/TMA block. Sections can be stained for in situ detection of protein, DNA or RNA targets using immunohistochemistry (IHC), fluorescent in situ hybridisation (FISH) or mRNA in situ hybridisation (RNA-ISH), respectively. This chapter provides detailed methods for CMA/TMA design, construction and analysis with in-depth notes on all technical aspects including tips to deal with common pitfalls the user may encounter. © Springer Science+Business Media, LLC 2011.

An evaluation of tumor oxygenation and gene expression in patients with early stage non-small cell lung cancers

Background: To directly assess tumor oxygenation in resectable non - small cell lung cancers (NSCLC) and to correlate tumor pO2 and the selected gene and protein expression to treatment outcomes. Methods: Twenty patients with resectable NSCLC were enrolled. Intraoperative measurements of normal lung and tumor pO2 were done with the Eppendorf polarographic electrode. All patients had plasma osteopontin measurements by ELISA. Carbonic anhydrase-IX (CA IX) staining of tumor sections was done in the majority of patients (n = 16), as was gene expression profiling (n = 12) using cDNA microarrays. Tumor pO2 was correlated with CA IX staining, osteopontin levels, and treatment outcomes. Results: The median tumor pO2 ranged from 0.7 to 46 mm Hg (median, 16.6) and was lower than normal lung pO2 in all but one patient. Because both variables were affected by the completeness of lung deflation during measurement, we used the ratio of tumor/normal lung (T/L) pO2 as a reflection of tumor oxygenation. The median T/L pO 2 was 0.13. T/L pO2 correlated significantly with plasma osteopontin levels (r = 0.53, P = 0.02) and CA IX expression (P = 0.006). Gene expression profiling showed that high CD44 expression was a predictor for relapse, which was confirmed by tissue staining of CD44 variant 6 protein. Other variables associated with the risk of relapse were T stage (P = 0.02), T/L pO2 (P = 0.04), and osteopontin levels (P = 0.001). Conclusions: Tumor hypoxia exists in resectable NSCLC and is associated with elevated expression of osteopontin and CA IX. Tumor hypoxia and elevated osteopontin levels and CD44 expression correlated with poor prognosis. A larger study is needed to confirm the prognostic significance of these factors. © 2006 American Association for Cancer Research.

Transcriptome analyses of amoebic gill disease-affected atlantic salmon (Salmo salar) tissues reveal localized host gene suppression

The transcriptome response of Atlantic salmon (Salmo salar) displaying advanced stages of amoebic gill disease (AGD) was investigated. Naïve smolt were challenged with AGD for 19 days, at which time all fish were euthanized and their severity of infection quantified through histopathological scoring. Gene expression profiles were compared between heavily infected and naïve individuals using a 17 K Atlantic salmon cDNA microarray with real-time quantitative RT-PCR (qPCR) verification. Expression profiles were examined in the gill, anterior kidney, and liver. Twenty-seven transcripts were significantly differentially expressed within the gill; 20 of these transcripts were down-regulated in the AGD-affected individuals compared with naïve individuals. In contrast, only nine transcripts were significantly differentially expressed within the anterior kidney and five within the liver. Again the majority of these transcripts were down-regulated within the diseased individuals. A down-regulation of transcripts involved in apoptosis (procathepsin L, cathepsin H precursor, and cystatin B) was observed in AGD-affected Atlantic salmon. Four transcripts encoding genes with antioxidant properties also were down-regulated in AGD-affected gill tissue according to qPCR analysis. The most up-regulated transcript within the gill was an unknown expressed sequence tag (EST) whose expression was 218-fold (± SE 66) higher within the AGD affected gill tissue. Our results suggest that Atlantic salmon experiencing advanced stages of AGD demonstrate general down-regulation of gene expression, which is most pronounced within the gill. We propose that this general gene suppression is parasite-mediated, thus allowing the parasite to withstand or ameliorate the host response. © 2008 Springer Science+Business Media, LLC.

A promoter-level mammalian expression atlas

Regulated transcription controls the diversity, developmental pathways and spatial organization of the hundreds of cell types that make up a mammal. Using single-molecule cDNA sequencing, we mapped transcription start sites (TSSs) and their usage in human and mouse primary cells, cell lines and tissues to produce a comprehensive overview of mammalian gene expression across the human body. We find that few genes are truly 'housekeeping', whereas many mammalian promoters are composite entities composed of several closely separated TSSs, with independent cell-type-specific expression profiles. TSSs specific to different cell types evolve at different rates, whereas promoters of broadly expressed genes are the most conserved. Promoter-based expression analysis reveals key transcription factors defining cell states and links them to binding-site motifs. The functions of identified novel transcripts can be predicted by coexpression and sample ontology enrichment analyses. The functional annotation of the mammalian genome 5 (FANTOM5) project provides comprehensive expression profiles and functional annotation of mammalian cell-type-specific transcriptomes with wide applications in biomedical research.

Gene expression signatures of radiation response are specific, durable and accurate in mice and humans.

BACKGROUND: Previous work has demonstrated the potential for peripheral blood (PB) gene expression profiling for the detection of disease or environmental exposures. METHODS AND FINDINGS: We have sought to determine the impact of several variables on the PB gene expression profile of an environmental exposure, ionizing radiation, and to determine the specificity of the PB signature of radiation versus other genotoxic stresses. Neither genotype differences nor the time of PB sampling caused any lessening of the accuracy of PB signatures to predict radiation exposure, but sex difference did influence the accuracy of the prediction of radiation exposure at the lowest level (50 cGy). A PB signature of sepsis was also generated and both the PB signature of radiation and the PB signature of sepsis were found to be 100% specific at distinguishing irradiated from septic animals. We also identified human PB signatures of radiation exposure and chemotherapy treatment which distinguished irradiated patients and chemotherapy-treated individuals within a heterogeneous population with accuracies of 90% and 81%, respectively. CONCLUSIONS: We conclude that PB gene expression profiles can be identified in mice and humans that are accurate in predicting medical conditions, are specific to each condition and remain highly accurate over time.

Gene expression disruptions of organism versus organ in Drosophila species hybrids.

Hybrid dysfunctions, such as sterility, may result in part from disruptions in the regulation of gene expression. Studies of hybrids within the Drosophila simulans clade have reported genes expressed above or below the expression observed in their parent species, and such misexpression is associated with male sterility in multigenerational backcross hybrids. However, these studies often examined whole bodies rather than testes or had limited replication using less-sensitive but global techniques. Here, we use a new RNA isolation technique to re-examine hybrid gene expression disruptions in both testes and whole bodies from single Drosophila males by real-time quantitative RT-PCR. We find two early-spermatogenesis transcripts are underexpressed in hybrid whole-bodies but not in assays of testes alone, while two late-spermatogenesis transcripts seem to be underexpressed in both whole-bodies and testes alone. Although the number of transcripts surveyed is limited, these results provide some support for a previous hypothesis that the spermatogenesis pathway in these sterile hybrids may be disrupted sometime after the expression of the early meiotic arrest genes.

Intergenic and genic sequence lengths have opposite relationships with respect to gene expression.

Eukaryotic genomes are mostly composed of noncoding DNA whose role is still poorly understood. Studies in several organisms have shown correlations between the length of the intergenic and genic sequences of a gene and the expression of its corresponding mRNA transcript. Some studies have found a positive relationship between intergenic sequence length and expression diversity between tissues, and concluded that genes under greater regulatory control require more regulatory information in their intergenic sequences. Other reports found a negative relationship between expression level and gene length and the interpretation was that there is selection pressure for highly expressed genes to remain small. However, a correlation between gene sequence length and expression diversity, opposite to that observed for intergenic sequences, has also been reported, and to date there is no testable explanation for this observation. To shed light on these varied and sometimes conflicting results, we performed a thorough study of the relationships between sequence length and gene expression using cell-type (tissue) specific microarray data in Arabidopsis thaliana. We measured median gene expression across tissues (expression level), expression variability between tissues (expression pattern uniformity), and expression variability between replicates (expression noise). We found that intergenic (upstream and downstream) and genic (coding and noncoding) sequences have generally opposite relationships with respect to expression, whether it is tissue variability, median, or expression noise. To explain these results we propose a model, in which the lengths of the intergenic and genic sequences have opposite effects on the ability of the transcribed region of the gene to be epigenetically regulated for differential expression. These findings could shed light on the role and influence of noncoding sequences on gene expression.

A noisy linear map underlies oscillations in cell size and gene expression in bacteria.

During bacterial growth, a cell approximately doubles in size before division, after which it splits into two daughter cells. This process is subjected to the inherent perturbations of cellular noise and thus requires regulation for cell-size homeostasis. The mechanisms underlying the control and dynamics of cell size remain poorly understood owing to the difficulty in sizing individual bacteria over long periods of time in a high-throughput manner. Here we measure and analyse long-term, single-cell growth and division across different Escherichia coli strains and growth conditions. We show that a subset of cells in a population exhibit transient oscillations in cell size with periods that stretch across several (more than ten) generations. Our analysis reveals that a simple law governing cell-size control-a noisy linear map-explains the origins of these cell-size oscillations across all strains. This noisy linear map implements a negative feedback on cell-size control: a cell with a larger initial size tends to divide earlier, whereas one with a smaller initial size tends to divide later. Combining simulations of cell growth and division with experimental data, we demonstrate that this noisy linear map generates transient oscillations, not just in cell size, but also in constitutive gene expression. Our work provides new insights into the dynamics of bacterial cell-size regulation with implications for the physiological processes involved.

The Role of Thymidylate Synthase Induction in Modulating p53-regulated Gene Expression in Response to 5-Fluorouracil and Antifolates

Thymidylate synthase (TS) is a critical target for chemotherapeutic agents such as 5-fluorouracil (5-FU) and antifolates such as tomudex (TDX),multitargeted antifolate, and ZD9331. Using the MCF-7 breast cancer line, we have developed p53 wild-type (M7TS90) and null (M7TS90-E6) isogenic lines with inducible TS expression (approximately 6-fold induction compared with control after 48 h). In the M7TS90 line, inducible TS expression resulted in a moderate approximately 3-fold increase in 5-FU IC-50(72 h) dose and a dramatic >20-fold increase in the IC-50(72 h) doses of TDX, multitargeted antifolate, and ZD9331. S-phase cell cycle arrest and apoptosis induced by the antifolates were abrogated by TS induction. In contrast, cell cycle arrest and apoptosis induced by 5-FU was unaffected by TS expression levels. Inactivation of p53 significantly increased resistance to 5-FU and the antifolates with IC-50(72 h) doses for 5-FU and TDX of >100 and >10 microM, respectively, in the M7TS90-E6 cell line. Furthermore, p53 inactivation completely abrogated the cell cycle arrest and apoptosis induced by 5-FU. The antifolates induced S-phase arrest in the p53 null cell line; however, the induction of apoptosis by these agents was significantly reduced compared with p53 wild-type cells. Both inducible TS expression and the addition of exogenous thymidine (10 microM) blocked p53 and p21 induction by the antifolates but not by 5-FU in the M7TS90 cell line. Similarly, inducible TS expression and exogenous thymidine abrogated antifolate but not 5-FU-mediated up-regulation of Fas/CD95 in M7TS90 cells. Our results indicate that in M7TS90 cells, inducible TS expression modulates p53 and p53 target gene expression in response to TS-targeted antifolate therapies but not to 5-FU.

Cytosolic phospholipase A2 gene expression in rat mesangial cells is regulated post-transcriptionally

Cytosolic phospholipase A2 (cPLA2) is thought to be the rate-limiting enzyme in the arachidonic acid/eicosanoid cascade. The ability of various agonists to increase steady-state cPLA2 mRNA levels has previously been reported. The current study delineates the contributions of transcriptional and post-transcriptional processes to the regulation of cPLA2 gene expression in response to a variety of agonists in cultured rat glomerular mesangial cells. Epidermal growth factor, platelet-derived growth factor, serum and phorbol myristate acetate all increase the half-life of cPLA2 mRNA transcripts, indicating a role for post-transcriptional modulation of gene expression. The presence of three ATTTA motifs in the 3' untranslated region (3'UTR) of the rat cPLA2 cDNA is ascertained. Heterologous expression of chimeric constructs with different 3'UTRs ligated into the 3' end of the luciferase coding region reveals that the presence of the cPLA2 3'UTR results in reduced luciferase activity compared with constructs without the cPLA2 3'UTR. Furthermore, the luciferase activity in the constructs with the cPLA2 3'UTR is increased in response to the same agonists which stabilize endogenous cPLA2 mRNA. A negligible effect of these agonists on transcriptional control of cPLA2 is evident using promoter-reporter constructs expressed in transient and stable transfectants. Taken together, these results indicate predominant post-transcriptional regulation of cPLA2 mRNA levels.

IGFBP7 promoter methylation and gene expression analysis in prostate cancer

PURPOSE: IGFBP7 belongs to a family of insulin-like growth factor-1 regulatory binding proteins. IGFBP7 hypermethylation is associated with its down-regulation in various carcinomas. In prostate cancer IGFBP7 down-regulation has been widely reported but to our knowledge the mechanisms behind this event are unknown. We performed a denaturing high performance liquid chromatography screening and validation strategy to profile the methylation status of IGFBP7 in prostate cancer.

MATERIALS AND METHODS: We combined denaturing high performance liquid chromatography and bisulfite sequencing to examine IGFBP7 methylation in a panel of prostate cancer cell lines. Quantitative methylation specific polymerase chain reaction was used to determine methylation levels in prostate tissue specimens of primary prostate cancer, histologically benign prostate adjacent to tumor, high grade prostatic intraepithelial neoplasia and benign prostatic hyperplasia. IGFBP7 gene expression was measured by quantitative methylation specific polymerase chain reaction in cell lines and tissue specimens.

RESULTS: IGFBP7 was methylated in the 4 prostate cancer cell lines DU145, LNCaP, PC-3 and 22RV1. Quantitative methylation specific polymerase chain reaction analysis revealed that promoter methylation was associated with decreased IGFBP7 expression. Quantitative methylation specific polymerase chain reaction showed that IGFBP7 methylation was more frequently detected in prostate cancer (60% (31/52)) and high grade prostatic intraepithelial neoplasia (40% (6/15)) samples compared to histologically benign prostate adjacent to tumor (10%) and benign prostatic hyperplasia (0%) samples.

CONCLUSIONS: To our knowledge this is the first report of aberrant IGFBP7 promoter hypermethylation and concurrent IGFBP7 gene silencing in prostate cancer cell lines. Results demonstrate that CpG methylation of IGFBP7 may represent a novel biomarker of prostate cancer and pre-invasive neoplasms. Thus, future examination of IGFBP7 methylation and expression in a larger patient cohort, including bodily fluids, is justified to further evaluate its role in a diagnostic and prognostic setting.

Achieving hypoxia-inducible gene expression in tumors

Hypoxia is an inevitable feature of solid tumors and a common cause of treatment failure. Hypoxia acts as a trigger to genetic instability, apoptosis and possibly metastases. The adaptive response to cellular hypoxia involves the modulation of the synthesis of multiple proteins controlling processes such as glucose homeostasis, angiogenesis, vascular permeability and inflammation. The hypoxia responsive element (HRE) sequences isolated from oxygen-responsive genes have been shown to selectively induce gene expression in response to hypoxia when placed upstream of a promoter. The levels of induced gene expression were dependent on the number of HRE copies and the oxygen tension. Hypoxia-mediated cancer gene therapy strategies may represent a promising mean to significantly improve the efficacy of standard radiation therapy and chemotherapy approaches.

Histone deacetylase inhibitors suppress thymidylate synthase gene expression and synergize with the fluoropyrimidines in colon cancer cells

Despite recent therapeutic advances, the response rates to chemotherapy for patients with metastatic colon cancer remain at approximately 50% with the fluoropyrimidine, 5-fluorouracil (5-FU), continuing to serve as the foundation chemotherapeutic agent for the treatment of this disease. Previous studies have demonstrated that overexpression of thymidylate synthase (TS) is a key determinant of resistance to 5-FU-based chemotherapy. Therefore, there is a significant need to develop alternative therapeutic strategies to overcome TS-mediated resistance. In this study, we demonstrate that the histone deacetylase inhibitors (HDACi) vorinostat and LBH589 significantly downregulate TS gene expression in a panel of colon cancer cell lines. Downregulation of TS was independent of p53, p21 and HDAC2 expression and was achievable in vivo as demonstrated by mouse xenograft models. We provide evidence that HDACi treatment leads to a potent transcriptional repression of the TS gene. Combination of the fluoropyrimidines 5-FU or FUdR with both vorinostat and LBH589 enhanced cell cycle arrest and growth inhibition. Importantly, the downstream effects of TS inhibition were significantly enhanced by this combination including the inhibition of acute TS induction and the enhanced accumulation of the cytotoxic nucleotide intermediate dUTP. These data demonstrate that HDACi repress TS expression at the level of transcription and provides the first evidence suggesting a direct mechanistic link between TS downregulation and the synergistic interaction observed between HDACi and 5-FU. This study provides rationale for the continued clinical evaluation of HDACi in combination with 5-FU-based therapies as a strategy to overcome TS-mediated resistance.

A coactivator role of CARM1 in the dysregulation of β-catenin activity in colorectal cancer cell growth and gene expression

Aberrant activation of Wnt/β-catenin signaling, resulting in the expression of Wnt-regulated oncogenes, is recognized as a critical factor in the etiology of colorectal cancer. Occupancy of β-catenin at promoters of Wnt target genes drives transcription, but the mechanism of β-catenin action remains poorly understood. Here, we show that CARM1 (coactivator-associated arginine methyltransferase 1) interacts with β-catenin and positively modulates β-catenin-mediated gene expression. In colorectal cancer cells with constitutively high Wnt/β-catenin activity, depletion of CARM1 inhibits expression of endogenous Wnt/β-catenin target genes and suppresses clonal survival and anchorage-independent growth. We also identified a colorectal cancer cell line (RKO) with a low basal level of β-catenin, which is dramatically elevated by treatment with Wnt3a. Wnt3a also increased the expression of a subset of endogenous Wnt target genes, and CARM1 was required for the Wnt-induced expression of these target genes and the accompanying dimethylation of arginine 17 of histone H3. Depletion of β-catenin from RKO cells diminished the Wnt-induced occupancy of CARM1 on a Wnt target gene, indicating that CARM1 is recruited to Wnt target genes through its interaction with β-catenin and contributes to transcriptional activation by mediating events (including histone H3 methylation) that are downstream from the actions of β-catenin. Therefore, CARM1 is an important positive modulator of Wnt/β-catenin transcription and neoplastic transformation, and may thereby represent a novel target for therapeutic intervention in cancers involving aberrantly activated Wnt/β-catenin signaling.

Single-cell gene-expression profiling reveals qualitatively distinct CD8 T cells elicited by different gene-based vaccines.

CD8 T cells play a key role in mediating protective immunity against selected pathogens after vaccination. Understanding the mechanism of this protection is dependent upon definition of the heterogeneity and complexity of cellular immune responses generated by different vaccines. Here, we identify previously unrecognized subsets of CD8 T cells based upon analysis of gene-expression patterns within single cells and show that they are differentially induced by different vaccines. Three prime-boost vector combinations encoding HIV Env stimulated antigen-specific CD8 T-cell populations of similar magnitude, phenotype, and functionality. Remarkably, however, analysis of single-cell gene-expression profiles enabled discrimination of a majority of central memory (CM) and effector memory (EM) CD8 T cells elicited by the three vaccines. Subsets of T cells could be defined based on their expression of Eomes, Cxcr3, and Ccr7, or Klrk1, Klrg1, and Ccr5 in CM and EM cells, respectively. Of CM cells elicited by DNA prime-recombinant adenoviral (rAd) boost vectors, 67% were Eomes(-) Ccr7(+) Cxcr3(-), in contrast to only 7% and 2% stimulated by rAd5-rAd5 or rAd-LCMV, respectively. Of EM cells elicited by DNA-rAd, 74% were Klrk1(-) Klrg1(-)Ccr5(-) compared with only 26% and 20% for rAd5-rAd5 or rAd5-LCMV. Definition by single-cell gene profiling of specific CM and EM CD8 T-cell subsets that are differentially induced by different gene-based vaccines will facilitate the design and evaluation of vaccines, as well as enable our understanding of mechanisms of protective immunity.