Dissecting tBHQ induced ARE-driven gene expression through long and short oligonucleotide arrays

FMO1 gene expression independently predicts favorable recurrence-free survival of classical papillary thyroid cancer

Aim: To examine the expression profile of FMO1 in papillary thyroid cancer (PTC) and its prognostic value in recurrence-free survival (RFS). Methods: A retrospective analysis was performed using data from the Cancer Genome Atlas and Human Protein Atlas. Results: The most frequent variants of PTC had decreased FMO1 expression compared with their respective adjacent normal tissues. However, even under the best cut-off model, high FMO1 expression was only significantly associated with better RFS in classical PTC (p < 0.001), but not in other two variants. High FMO1 expression independently predicted favorable RFS (hazard ratio: 0.202; 95% CI: 0.084-0.487; p < 0.001) in classical PTC. Conclusion: High FMO1 expression might serve as a biomarker that independently predicts favorable RFS in classical PTC patients.

Nrf2 and Nrf2-related proteins in development and developmental toxicity: Insights from studies in zebrafish (Danio rerio)

Oxidative stress is an important mechanism of chemical toxicity, contributing to developmental toxicity and teratogenesis as well as to cardiovascular and neurodegenerative diseases and diabetic embryopathy. Developing animals are especially sensitive to effects of chemicals that disrupt the balance of processes generating reactive species and oxidative stress, and those anti-oxidant defenses that protect against oxidative stress. The expression and inducibility of anti-oxidant defenses through activation of NFE2-related factor 2 (Nrf2) and related proteins is an essential process affecting the susceptibility to oxidants, but the complex interactions of Nrf2 in determining embryonic response to oxidants and oxidative stress are only beginning to be understood. The zebrafish (Danio rerio) is an established model in developmental biology and now also in developmental toxicology and redox signaling. Here we review the regulation of genes involved in protection against oxidative stress in developing vertebrates, with a focus on Nrf2 and related cap'n'collar (CNC)-basic-leucine zipper (bZIP) transcription factors. Vertebrate animals including zebrafish share Nfe2, Nrf1, Nrf2, and Nrf3 as well as a core set of genes that respond to oxidative stress, contributing to the value of zebrafish as a model system with which to investigate the mechanisms involved in regulation of redox signaling and the response to oxidative stress during embryolarval development. Moreover, studies in zebrafish have revealed nrf and keap1 gene duplications that provide an opportunity to dissect multiple functions of vertebrate NRF genes, including multiple sensing mechanisms involved in chemical-specific effects.

Comparative analysis of NRF2-responsive gene expression in AcPC-1 pancreatic cancer cell line

NRF2 is a nuclear transcription factor activated in response to oxidative stress and related with metabolizing of xenotoxic materials and ABC transporter mediated drug resistance. We studied the expression of mRNAs under the siRNA-mediated knockdown of NRF2 and tBHQ-treated condition in AsPC-1 metastatic pancreatic cancer cell line to understand the AsPC-1 specific role(s) of NRF2 and further to investigate the relationship between drug resistance and metastatic plasticity and mobility of AsPc1. Here we show that the genes of aldo–keto reductases, cytochrome P450 family, aldehyde dehydrogenase, thioredoxin reductase, ABC transporter and epoxide hydrolase responsible for drug metabolism or oxidative stress concisely responded to NRF2 stabilization and knockdown of NRF2. In addition the expression of PIR, a candidate of oncogene and KISS1, a suppressor of metastasis were affected by NRF2 stabilization and knockdown. Our result provide comprehensive understanding of NRF2 target genes of drug response, oxidative stress response and metastasis in AsPc-1 metastatic pancreatic cancer cell line.

The transcriptional response to oxidative stress during vertebrate development: effects of tert-butylhydroquinone and 2,3,7,8-tetrachlorodibenzo-p-dioxin

Oxidative stress is an important mechanism of chemical toxicity, contributing to teratogenesis and to cardiovascular and neurodegenerative diseases. Developing animals may be especially sensitive to chemicals causing oxidative stress. The developmental expression and inducibility of anti-oxidant defenses through activation of NF-E2-related factor 2 (NRF2) affect susceptibility to oxidants, but the embryonic response to oxidants is not well understood. To assess the response to chemically mediated oxidative stress and how it may vary during development, zebrafish embryos, eleutheroembryos, or larvae at 1, 2, 3, 4, 5, and 6 days post fertilization (dpf) were exposed to DMSO (0.1%), tert-butylhydroquinone (tBHQ; 10 µM) or 2,3,7,8-tetrachlorodibenzo-p-dioxin (TCDD; 2 nM) for 6 hr. Transcript abundance was assessed by real-time qRT-PCR and microarray. qRT-PCR showed strong (4- to 5-fold) induction of gstp1 by tBHQ as early as 1 dpf. tBHQ also induced gclc (2 dpf), but not sod1, nqo1, or cyp1a. TCDD induced cyp1a but none of the other genes. Microarray analysis showed that 1477 probes were significantly different among the DMSO-, tBHQ-, and TCDD-treated eleutheroembryos at 4 dpf. There was substantial overlap between genes induced in developing zebrafish and a set of marker genes induced by oxidative stress in mammals. Genes induced by tBHQ in 4-dpf zebrafish included those involved in glutathione synthesis and utilization, signal transduction, and DNA damage/stress response. The strong induction of hsp70 determined by microarray was confirmed by qRT-PCR and by use of transgenic zebrafish expressing enhanced green fluorescent protein (EGFP) under control of the hsp70 promoter. Genes strongly down-regulated by tBHQ included mitfa, providing a molecular explanation for the loss of pigmentation in tBHQ-exposed embryos. These data show that zebrafish embryos are responsive to oxidative stress as early as 1 dpf, that responsiveness varies with development in a gene-specific manner, and that the oxidative stress response is substantially conserved in vertebrate animals.

Resveratrol inhibits estrogen-induced breast carcinogenesis through induction of NRF2-mediated protective pathways

The importance of estrogens in the etiology of breast cancer is widely recognized. Estrogen-induced oxidative stress has been implicated in this carcinogenic process. Resveratrol (Res), a natural antioxidant phytoestrogen has chemopreventive effects against a variety of illnesses including cancer. The objective of the present study was to characterize the mechanism(s) of Res-mediated protection against estrogen-induced breast carcinogenesis. Female August Copenhagen Irish rats were treated with 17β-estradiol (E2), Res and Res + E2 for 8 months. Cotreatment of rats with Res and E2 inhibited E2-mediated proliferative changes in mammary tissues and significantly increased tumor latency and reduced E2-induced breast tumor development. Resveratrol treatment alone or in combination with E2 significantly upregulated expression of nuclear factor erythroid 2-related factor 2 (NRF2) in mammary tissues. Expression of NRF2-regulated antioxidant genes NQO1, SOD3 and OGG1 that are involved in protection against oxidative DNA damage was increased in Res- and Res + E2-treated mammary tissues. Resveratrol also prevented E2-mediated inhibition of detoxification genes AOX1 and FMO1. Inhibition of E2-mediated alterations in NRF2 promoter methylation and expression of NRF2 targeting miR-93 after Res treatment indicated Res-mediated epigenetic regulation of NRF2 during E2-induced breast carcinogenesis. Resveratrol treatment also induced apoptosis and inhibited E2-mediated increase in DNA damage in mammary tissues. Increased apoptosis and decreased DNA damage, cell migration, colony and mammosphere formation in Res- and Res + E2-treated MCF-10A cells suggested a protective role of Res against E2-induced mammary carcinogenesis. Small-interfering RNA-mediated silencing of NRF2 inhibited Res-mediated preventive effects on the colony and mammosphere formation. Taken together, these results suggest that Res inhibits E2-induced breast carcinogenesis via induction of NRF2-mediated protective pathways.

Identification of alternative splicing events regulated by the oncogenic factor SRSF1 in lung cancer

Abnormal alternative splicing has been associated with cancer. Genome-wide microarrays can be used to detect differential splicing events. In this study, we have developed ExonPointer, an algorithm that uses data from exon and junction probes to identify annotated cassette exons. We used the algorithm to profile differential splicing events in lung adenocarcinoma A549 cells after downregulation of the oncogenic serine/arginine-rich splicing factor 1 (SRSF1). Data were generated using two different microarray platforms. The PCR-based validation rate of the top 20 ranked genes was 60% and 100%. Functional enrichment analyses found a substantial number of splicing events in genes related to RNA metabolism. These analyses also identified genes associated with cancer and developmental and hereditary disorders, as well as biologic processes such as cell division, apoptosis, and proliferation. Most of the top 20 ranked genes were validated in other adenocarcinoma and squamous cell lung cancer cells, with validation rates of 80% to 95% and 70% to 75%, respectively. Moreover, the analysis allowed us to identify four genes, ATP11C, IQCB1, TUBD1, and proline-rich coiled-coil 2C (PRRC2C), with a significantly different pattern of alternative splicing in primary non-small cell lung tumors compared with normal lung tissue. In the case of PRRC2C, SRSF1 downregulation led to the skipping of an exon overexpressed in primary lung tumors. Specific siRNA downregulation of the exon-containing variant significantly reduced cell growth. In conclusion, using a novel analytical tool, we have identified new splicing events regulated by the oncogenic splicing factor SRSF1 in lung cancer.

Alzheimer's disease: redox dysregulation as a common denominator for diverse pathogenic mechanisms

Alzheimer's disease (AD) is the most common cause of dementia and a progressive neurodegeneration that appears to result from multiple pathogenic mechanisms (including protein misfolding/aggregation, involved in both amyloid β-dependent senile plaques and tau-dependent neurofibrillary tangles), metabolic and mitochondrial dysfunction, excitoxicity, calcium handling impairment, glial cell dysfunction, neuroinflammation, and oxidative stress. Oxidative stress, which could be secondary to several of the other pathophysiological mechanisms, appears to be a major determinant of the pathogenesis and progression of AD. The identification of oxidized proteins common for mild cognitive impairment and AD suggests that key oxidation pathways are triggered early and are involved in the initial progression of the neurodegenerative process. Abundant data support that oxidative stress, also considered as a main factor for aging, the major risk factor for AD, can be a common key element capable of articulating the divergent nature of the proposed pathogenic factors. Pathogenic mechanisms influence each other at different levels. Evidence suggests that it will be difficult to define a single-target therapy resulting in the arrest of progression or the improvement of AD deterioration. Since oxidative stress is present from early stages of disease, it appears as one of the main targets to be included in a clinical trial. Exploring the articulation of AD pathogenic mechanisms by oxidative stress will provide clues for better understanding the pathogenesis and progression of this dementing disorder and for the development of effective therapies to treat this disease.

Treatment-induced oxidative stress and cellular antioxidant capacity determine response to bortezomib in mantle cell lymphoma

PURPOSE

Proteasome inhibition disrupts protein homeostasis and induces apoptosis. Up to 50% of patients with relapsed mantle cell lymphoma (MCL) respond to bortezomib. We used gene expression profiling to investigate the connection between proteasome inhibition, cellular response, and clinical efficacy.

EXPERIMENTAL DESIGN

We assessed transcriptional changes in primary tumor cells from five patients during treatment with bortezomib in vivo, and in 10 MCL cell lines exposed to bortezomib in vitro, on Affymetrix microarrays. Key findings were confirmed by western blotting.

RESULTS

MCL cell lines exposed to bortezomib in vitro showed upregulation of endoplasmic reticulum and oxidative stress response pathways. Gene expression changes were strongest in bortezomib-sensitive cells and these cells were also more sensitive to oxidative stress induced by H2O2. Purified tumor cells obtained at several timepoints during bortezomib treatment in 5 previously untreated patients with leukemic MCL showed strong activation of the antioxidant response controlled by NRF2. Unexpectedly, activation of this homeostatic program was significantly stronger in tumors with the best clinical response. Consistent with its proapoptotic function, we found upregulation of NOXA in circulating tumor cells of responding patients. In resistant cells, gene expression changes in response to bortezomib were limited and upregulation of NOXA was absent. Interestingly, at baseline, bortezomib-resistant cells displayed a relatively higher expression of the NRF2 gene-expression signature than sensitive cells (P < 0.001).

CONCLUSION

Bortezomib triggers an oxidative stress response in vitro and in vivo. High cellular antioxidant capacity contributes to bortezomib resistance.

Gene expression analysis of RNA purified from embryonic stem cells and embryoid body-derived cells using a high-throughput microarray platform

In this unit, starting with purified RNA, experimental protocols for performing microarray expression analysis of embryonic stem cell lines compared to their corresponding differentiated embryocidal bodies are described. Methods for data analysis are suggested, with the goal of determining which genes are differentially expressed between the preparations. As an example, the use of the Affymetrix microarray expression platform is described, but alternative experimental options for analysis of RNA transcript levels are also summarized. This unit suggests quality control metrics, summarizes the critical parameters necessary for obtaining reproducible experimental results, and outlines quantitative PCR methods for validating microarray results.

Identification of a peripheral blood transcriptional biomarker panel associated with operational renal allograft tolerance

Long-term allograft survival generally requires lifelong immunosuppression (IS). Rarely, recipients display spontaneous "operational tolerance" with stable graft function in the absence of IS. The lack of biological markers of this phenomenon precludes identification of potentially tolerant patients in which IS could be tapered and hinders the development of new tolerance-inducing strategies. The objective of this study was to identify minimally invasive blood biomarkers for operational tolerance and use these biomarkers to determine the frequency of this state in immunosuppressed patients with stable graft function. Blood gene expression profiles from 75 renal-transplant patient cohorts (operational tolerance/acute and chronic rejection/stable graft function on IS) and 16 healthy individuals were analyzed. A subset of samples was used for microarray analysis where three-class comparison of the different groups of patients identified a "tolerant footprint" of 49 genes. These biomarkers were applied for prediction of operational tolerance by microarray and real-time PCR in independent test groups. Thirty-three of 49 genes correctly segregated tolerance and chronic rejection phenotypes with 99% and 86% specificity. The signature is shared with 1 of 12 and 5 of 10 stable patients on triple IS and low-dose steroid monotherapy, respectively. The gene signature suggests a pattern of reduced costimulatory signaling, immune quiescence, apoptosis, and memory T cell responses. This study identifies in the blood of kidney recipients a set of genes associated with operational tolerance that may have utility as a minimally invasive monitoring tool for guiding IS titration. Further validation of this tool for safe IS minimization in prospective clinical trials is warranted.

Three methods for optimization of cross-laboratory and cross-platform microarray expression data

Microarray gene expression data becomes more valuable as our confidence in the results grows. Guaranteeing data quality becomes increasingly important as microarrays are being used to diagnose and treat patients (1–4). The MAQC Quality Control Consortium, the FDA's Critical Path Initiative, NCI's caBIG and others are implementing procedures that will broadly enhance data quality. As GEO continues to grow, its usefulness is constrained by the level of correlation across experiments and general applicability. Although RNA preparation and array platform play important roles in data accuracy, pre-processing is a user-selected factor that has an enormous effect. Normalization of expression data is necessary, but the methods have specific and pronounced effects on precision, accuracy and historical correlation. As a case study, we present a microarray calibration process using normalization as the adjustable parameter. We examine the impact of eight normalizations across both Agilent and Affymetrix expression platforms on three expression readouts: (1) sensitivity and power, (2) functional/biological interpretation and (3) feature selection and classification error. The reader is encouraged to measure their own discordant data, whether cross-laboratory, cross-platform or across any other variance source, and to use their results to tune the adjustable parameters of their laboratory to ensure increased correlation.

Identification of a predictive gene expression signature of cervical lymph node metastasis in oral squamous cell carcinoma

An accurate assessment of the cervical lymph node metastasis status in oral cavity cancer not only helps predict the prognosis of patients, but also helps surgeons to perform the appropriate treatment. We investigated the utilization of microarray technology focusing on the differences in gene expression profiles between primary tumors of oral squamous cell carcinoma that had metastasized to cervical lymph nodes and those that had not metastasized in the hope of finding new biomarkers to serve for diagnosis and treatment of oral cavity cancer. To design this experiment, we prepared two groups: the learning case group with 30 patients and the test case group with 13 patients. All tissue samples were performed using laser captured microdissection to yield cancer cells, and RNA was isolated from purified cancer cells. To identify a predictive gene expression signature, the different gene expressions between the two groups with and without metastasis in the learning case (n = 30) were analyzed, and the 85 genes expressed differentially were selected. Subsequently, to construct a more accurate prediction model, we further selected the genes with a high power for prediction from the 85 genes using the AdaBoost algorithm. The eight candidate genes, DCTD, IL-15, THBD, GSDML, SH3GL3, PTHLH, RP5-1022P6 and C9orf46, were selected to achieve the minimum error rate. Quantitative reverse transcription-polymerase chain reaction was carried out to validate the selected genes. From these statistical methods, the prediction model was constructed including the eight genes and this model was evaluated by using the test case group. The results in 12 of 13 cases ( approximately 92.3%) were predicted correctly.

Experimental comparison and cross-validation of the Affymetrix and Illumina gene expression analysis platforms

The growth in popularity of RNA expression microarrays has been accompanied by concerns about the reliability of the data especially when comparing between different platforms. Here, we present an evaluation of the reproducibility of microarray results using two platforms, Affymetrix GeneChips and Illumina BeadArrays. The study design is based on a dilution series of two human tissues (blood and placenta), tested in duplicate on each platform. The results of a comparison between the platforms indicate very high agreement, particularly for genes which are predicted to be differentially expressed between the two tissues. Agreement was strongly correlated with the level of expression of a gene. Concordance was also improved when probes on the two platforms could be identified as being likely to target the same set of transcripts of a given gene. These results shed light on the causes or failures of agreement across microarray platforms. The set of probes we found to be most highly reproducible can be used by others to help increase confidence in analyses of other data sets using these platforms.

Gene expression evidence for remodeling of lateral hypothalamic circuitry in cocaine addiction

By using high-density oligonucleotide arrays, we profiled gene expression in reward-related brain regions of rats that developed escalated cocaine intake after extended access to cocaine (6 h per day). Rats allowed restricted daily access to cocaine (only 1 h) that displayed a stable level of cocaine intake and cocaine naive rats were used for controls. Four analysis methods were compared: Affymetrix microarray suite 4 and microarray suite 5, which use perfect-match-minus-mismatch models, and dchip and rma, which use perfect-match-only models to generate expression values. Results were validated by RT-PCR in individual animals from an independent replication of the experiment. A small number of genes was associated with escalated cocaine intake (ESC genes). Unexpectedly, of the brain regions examined [prefrontal cortex, nucleus accumbens, septum, lateral hypothalamus (LH), amygdala, and ventral tegmental area], the LH was the most transcriptionally responsive in escalation of cocaine intake. Most of the ESC genes identified are also expressed during synaptogenesis and synaptic plasticity and include genes that code for several presynaptic and postsynaptic proteins involved in neurotransmission. These results suggest that LH intrinsic circuitry undergoes a structural reorganization during escalation of cocaine use. This remodeling of LH circuitry could contribute to the chronic deficit in reward function that has been hypothesized to drive the transition to drug addiction. Results also support the value of using multiple analysis strategies to identify the most robust changes in gene expression and to compensate for the biases that affect each strategy.