A transcription factor map as revealed by a genome-wide gene expression analysis of whole-blood mRNA transcriptome in multiple sclerosis

Background Several lines of evidence suggests that transcription factors are involved in the pathogenesis of Multiple Sclerosis (MS) but a complete mapping the whole network has been elusive. One of the reasons is that there are several clinical subtypes of MS and transcription factors which may be involved in one subtype may not be in others. We investigated the possibility that this network could be mapped using microarray technologies and modern bioinformatics methods on a dataset from whole blood in 99 untreated MS patients (36 Relapse Remitting MS, 43 Primary Progressive MS, and 20 Secondary Progressive MS) and 45 age-matched healthy controls, Methodology/Principal Findings We have used two different analytical methodologies: a differential expression analysis and a differential co-expression analysis, which have converged on a significant number of regulatory motifs that seem to be statistically overrepresented in genes which are either differentially expressed (or differentially co-expressed) in cases and controls (e.g. V$KROX_Q6, p-value < 3.31E-6; V$CREBP1_Q2, p-value < 9.93E-6, V$YY1_02, p-value < 1.65E-5). Conclusions/significance: Our analysis uncovered a network of transcription factors that potentially dysregulate several genes in MS or one or more of its disease subtypes. Analysing the published literature we have found that these transcription factors are involved in the early T-lymphocyte specification and commitment as well as in oligodendrocytes dedifferentiation and development. The most significant transcription factors motifs were for the Early Growth response EGR/KROX family, ATF2, YY1 (Yin and Yang 1), E2F-1/DP-1 and E2F-4/DP-2 heterodimers, SOX5, and CREB and ATF families.

Ethanol exposure alters monoamine oxidase gene-expression in neuronal SH-SY5Y cells

Abstract: Monoamine Oxidase (MAO) enzymes catabolise, and thus modulate abundance of, neurotransmitters in the brain. Variation in MAO enzyme activity has been linked to alcohol abuse behaviour, although the molecular mechanisms underlying this association are not understood. The present study evaluated relative gene-transcript abundance of MAO-A and MAO-B in the SH-SY5Y human neuroblastoma cell-line in response to ethanol exposure and following ethanol withdrawal. We found that each isoform of MAO was significantly transcriptionally up-regulated 55-80% in response to 100mM ethanol exposure. This trend was maintained following prolonged exposures (24 h-72 h) and with short exposures (24 h) followed by a period of ethanol withdrawal, suggesting that the transcriptional regulation is the result of a cellular change occurring within the first 24 hours of ethanol exposure. These results suggest a role for MAO transcriptional regulation in the complex neurobiochemical changes underlying alcohol addiction.

Global gene expression in skeletal muscle from well-trained strength and endurance athletes

PURPOSE: We used gene microarray analysis to compare the global expression profile of genes involved in adaptation to training in skeletal muscle from chronically strength-trained (ST), endurance-trained (ET), and untrained control subjects (Con). METHODS: Resting skeletal muscle samples were obtained from the vastus lateralis of 20 subjects (Con n = 7, ET n = 7, ST n = 6; trained [TR] groups >8 yr specific training). Total RNA was extracted from tissue for two color microarray analysis and quantative (Q)-PCR. Trained subjects were characterized by performance measures of peak oxygen uptake V?O 2peak) on a cycle ergometer and maximal concentric and eccentric leg strength on an isokinetic dynamometer. RESULTS: Two hundred and sixty-three genes were differentially expressed in trained subjects (ET + ST) compared with Con (P < 0.05), whereas 21 genes were different between ST and ET (P < 0.05). These results were validated by reverse transcriptase polymerase chain reaction for six differentially regulated genes (EIFSJ, LDHB, LMO4, MDH1, SLC16A7, and UTRN. Manual cluster analyses revealed significant regulation of genes involved in muscle structure and development in TR subjects compared with Con (P < 0.05) and expression correlated with measures of performance (P < 0.05). ET had increased whereas ST had decreased expression of gene clusters related to mitochondrial/oxidative capacity (P ?‰Currency sign 0.05). These mitochondrial gene clusters correlated with V?O2peak (P < 0.05). V?O2peak also correlated with expression of gene clusters that regulate fat and carbohydrate oxidation (P < 0.05). CONCLUSION: We demonstrate that chronic training subtly coregulates numerous genes from important functional groups that may be part of the long-term adaptive process to adapt to repeated training stimuli.

Analysis of kinase gene expression in the frontal cortex of suicide victims: Implications of fear and stress

Suicide is a serious public health issue that results from an interaction between multiple risk factors including individual vulnerabilities to complex feelings of hopelessness, fear, and stress. Although kinase genes have been implicated in fear and stress, including the consolidation and extinction of fearful memories, expression profiles of those genes in the brain of suicide victims are less clear. Using gene expression microarray data from the Online Stanley Genomics Database 1 and a quantitative PCR, we investigated the expression profiles of multiple kinase genes including the calcium calmodulin-dependent kinase (CAMK), the cyclin-dependent kinase, the mitogen-activated protein kinase (MAPK), and the protein kinase C (PKC) in the prefrontal cortex (PFC) of mood disorder patients died with suicide (N = 45) and without suicide (N = 38). We also investigated the expression pattern of the same genes in the PFC of developing humans ranging in age from birth to 49 year (N = 46). The expression levels of CAMK2B, CDK5, MAPK9, and PRKCI were increased in the PFC of suicide victims as compared to non-suicide controls (false discovery rate, FDR-adjusted p < 0.05, fold change >1.1). Those genes also showed changes in expression pattern during the postnatal development (FDR-adjusted p < 0.05). These results suggest that multiple kinase genes undergo age-dependent changes in normal brains as well as pathological changes in suicide brains. These findings may provide an important link to protein kinases known to be important for the development of fear memory, stress associated neural plasticity, and up-regulation in the PFC of suicide victims. More research is needed to better understand the functional role of these kinase genes that may be associated with the pathophysiology of suicide

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.

Gene expression analysis of diagnostic biopsies predicts pathological response to neoadjuvant chemoradiotherapy of esophageal cancer

OBJECTIVE: This study explored gene expression differences in predicting response to chemoradiotherapy in esophageal cancer. PURPOSE:: A major pathological response to neoadjuvant chemoradiation is observed in about 40% of esophageal cancer patients and is associated with favorable outcomes. However, patients with tumors of similar histology, differentiation, and stage can have vastly different responses to the same neoadjuvant therapy. This dichotomy may be due to differences in the molecular genetic environment of the tumor cells. BACKGROUND DATA: Diagnostic biopsies were obtained from a training cohort of esophageal cancer patients (13), and extracted RNA was hybridized to genome expression microarrays. The resulting gene expression data was verified by qRT-PCR. In a larger, independent validation cohort (27), we examined differential gene expression by qRT-PCR. The ability of differentially-regulated genes to predict response to therapy was assessed in a multivariate leave-one-out cross-validation model. RESULTS: Although 411 genes were differentially expressed between normal and tumor tissue, only 103 genes were altered between responder and non-responder tumor; and 67 genes differentially expressed >2-fold. These included genes previously reported in esophageal cancer and a number of novel genes. In the validation cohort, 8 of 12 selected genes were significantly different between the response groups. In the predictive model, 5 of 8 genes could predict response to therapy with 95% accuracy in a subset (74%) of patients. CONCLUSIONS: This study has identified a gene microarray pattern and a set of genes associated with response to neoadjuvant chemoradiation in esophageal cancer. The potential of these genes as biomarkers of response to treatment warrants further investigation. Copyright © 2009 by Lippincott Williams & Wilkins.

The use of artificial MicroRNA technology to control gene expression in Arabidopsis thaliana

In plants, double-stranded RNA (dsRNA) is an effective trigger of RNA silencing, and several classes of endogenous small RNA (sRNA), processed from dsRNA substrates by DICER-like (DCL) endonucleases, are essential in controlling gene expression. One such sRNA class, the microRNAs (miRNAs) control the expression of closely related genes to regulate all aspects of plant development, including the determination of leaf shape, leaf polarity, flowering time, and floral identity. A single miRNA sRNA silencing signal is processed from a long precursor transcript of nonprotein-coding RNA, termed the primary miRNA (pri-miRNA). A region of the pri-miRNA is partially self-complementary allowing the transcript to fold back onto itself to form a stem-loop structure of imperfectly dsRNA. Artificial miRNA (amiRNA) technology uses endogenous pri-miRNAs, in which the miRNA and miRNA*(passenger strand of the miRNA duplex) sequences have been replaced with corresponding amiRNA/ amiRNA*sequences that direct highly efficient RNA silencing of the targeted gene. Here, we describe the rules for amiRNA design, as well as outline the PCR and bacterial cloning procedures involved in the construction of an amiRNA plant expression vector to control target gene expression in Arabidopsis thaliana. © 2014 Springer Science+Business Media New York.

Validation of β-actin used as endogenous control for gene expression analysis in mechanobiology studies : amendments

We write in response to the letter by Liu et al. [1] commenting on our article, ‘‘Mesenchymal Stem Cells Regulate Angiogenesis According to Their Mechanical Environment’’ [2]. The study by Liu et al. demonstrates that the commonly used endogeneous reference gene (ERG), b-actin, is upregulated by mechanical loading, indicating a potential bias in the determined target gene expression when normalizing to b-actin, such as in our report on unchanged vascular endothelial growth factor (VEGF) and hypoxia-inducible factors (HIF)-1a mRNA levels in mechanically loaded mesenchymal stem cells (MSCs).

β-Arrestin2 regulates RANKL and Ephrins gene expression in response to bone remodeling in mice

PTH-stimulated intracellular signaling is regulated by the cytoplasmic adaptor molecule barrestin. We reported that the response of cancellous bone to intermittent PTH is reduced in b-arrestin22/2 mice and suggested that b-arrestins could influence the bone mineral balance by controlling RANKL and osteoprotegerin (OPG) gene expression. Here, we study the role of b-arrestin2 on the in vitro development and activity of bone marrow (BM) osteoclasts (OCs) and Ephrins ligand (Efn), and receptor (Eph) mRNA levels in bone in response to PTH and the changes of bone microarchitecture in wildtype (WT) and barrestin2 2/2 mice in models of bone remodeling: a low calcium diet (LoCa) and ovariectomy (OVX). The number of PTH-stimulated OCs was higher in BM cultures from b-arrestin22/2 compared with WT, because of a higher RANKL/OPG mRNA and protein ratio, without directly influencing osteoclast activity. In vivo, high PTH levels induced by LoCa led to greater changes in TRACP5b levels in b-arrestin22/2 compared with WT. LoCa caused a loss of BMD and bone microarchitecture, which was most prominent in b-arrestin22/2. PTH downregulated Efn and Eph genes in b-arrestin22/2, but not WT. After OVX, vertebral trabecular bone volume fraction and trabecular number were lower in b-arrestin22/2 compared with WT. Histomorphometry showed that OC number was higher in OVX-b-arrestin22/2 compared with WT. These results indicate that b-arrestin2 inhibits osteoclastogenesis in vitro, which resulted in decreased bone resorption in vivo by regulating RANKL/OPG production and ephrins mRNAs. As such, b-arrestins should be considered an important mechanism for the control of bone remodeling in response to PTH and estrogen deprivation.

Direct repression of MYB by ZEB1 suppresses proliferation and epithelial gene expression during epithelial-to-mesenchymal transition of breast cancer cells

Introduction Epithelial-to-mesenchymal transition (EMT) promotes cell migration and is important in metastasis. Cellular proliferation is often downregulated during EMT, and the reverse transition (MET) in metastases appears to be required for restoration of proliferation in secondary tumors. We studied the interplay between EMT and proliferation control by MYB in breast cancer cells. Methods MYB, ZEB1, and CDH1 expression levels were manipulated by lentiviral small-hairpin RNA (shRNA)-mediated knockdown/overexpression, and verified with Western blotting, immunocytochemistry, and qRT-PCR. Proliferation was assessed with bromodeoxyuridine pulse labeling and flow cytometry, and sulforhodamine B assays. EMT was induced with epidermal growth factor for 9 days or by exposure to hypoxia (1% oxygen) for up to 5 days, and assessed with qRT-PCR, cell morphology, and colony morphology. Protein expression in human breast cancers was assessed with immunohistochemistry. ZEB1-MYB promoter binding and repression were determined with Chromatin Immunoprecipitation Assay and a luciferase reporter assay, respectively. Student paired t tests, Mann–Whitney, and repeated measures two-way ANOVA tests determined statistical significance (P < 0.05). Results Parental PMC42-ET cells displayed higher expression of ZEB1 and lower expression of MYB than did the PMC42-LA epithelial variant. Knockdown of ZEB1 in PMC42-ET and MDA-MB-231 cells caused increased expression of MYB and a transition to a more epithelial phenotype, which in PMC42-ET cells was coupled with increased proliferation. Indeed, we observed an inverse relation between MYB and ZEB1 expression in two in vitro EMT cell models, in matched human breast tumors and lymph node metastases, and in human breast cancer cell lines. Knockdown of MYB in PMC42-LA cells (MYBsh-LA) led to morphologic changes and protein expression consistent with an EMT. ZEB1 expression was raised in MYBsh-LA cells and significantly repressed in MYB-overexpressing MDA-MB-231 cells, which also showed reduced random migration and a shift from mesenchymal to epithelial colony morphology in two dimensional monolayer cultures. Finally, we detected binding of ZEB1 to MYB promoter in PMC42-ET cells, and ZEB1 overexpression repressed MYB promoter activity. Conclusions This work identifies ZEB1 as a transcriptional repressor of MYB and suggests a reciprocal MYB-ZEB1 repressive relation, providing a mechanism through which proliferation and the epithelial phenotype may be coordinately modulated in breast cancer cells.

Assessment of gene expression of intracellular calcium channels, pumps and exchangers with epidermal growth factor-induced epithelial-mesenchymal transition in a breast cancer cell line

Background Epithelial-mesenchymal transition (EMT) is a process implicated in cancer metastasis that involves the conversion of epithelial cells to a more mesenchymal and invasive cell phenotype. In breast cancer cells EMT is associated with altered store-operated calcium influx and changes in calcium signalling mediated by activation of cell surface purinergic receptors. In this study, we investigated whether MDA-MB-468 breast cancer cells induced to undergo EMT exhibit changes in mRNA levels of calcium channels, pumps and exchangers located on intracellular calcium storing organelles, including the Golgi, mitochondria and endoplasmic reticulum (ER). Methods Epidermal growth factor (EGF) was used to induce EMT in MDA-MB-468 breast cancer cells. Serum-deprived cells were treated with EGF (50 ng/mL) for 12 h and gene expression was assessed using quantitative RT-PCR. Results and conclusions These data reveal no significant alterations in mRNA levels of the Golgi calcium pump secretory pathway calcium ATPases (SPCA1 and SPCA2), or the mitochondrial calcium uniporter (MCU) or Na+/Ca2+ exchanger (NCLX). However, EGF-induced EMT was associated with significant alterations in mRNA levels of specific ER calcium channels and pumps, including (sarco)-endoplasmic reticulum calcium ATPases (SERCAs), and inositol 1,4,5-trisphosphate receptor (IP3R) and ryanodine receptor (RYR) calcium channel isoforms. The most prominent change in gene expression between the epithelial and mesenchymal-like states was RYR2, which was enriched 45-fold in EGF-treated MDA-MB-468 cells. These findings indicate that EGF-induced EMT in breast cancer cells may be associated with major alterations in ER calcium homeostasis.

Transient gene expression in Medicago truncatula leaves via agroinfiltration

Transient expression is a powerful method for the functional characterization of genes. In this chapter, we outline a protocol for the transient expression of constructs in Medicago truncatula leaves using Agrobacterium tumefaciens infiltration. Using quantitative real-time PCR we demonstrate that the infiltration of a construct containing the LEGUME ANTHOCYANIN PRODUCTION 1 (LAP1) transcription factor results in the strong upregulation of key biosynthetic genes and the accumulation of anthocyanin pigment in the leaves after just 3 days. Thus, this method provides a rapid and powerful way to the discovery of downstream targets of M. truncatula transcription factors.

Effect of 5′UTR introns on gene expression in Arabidopsis thaliana

Background The majority of introns in gene transcripts are found within the coding sequences (CDSs). A small but significant fraction of introns are also found to reside within the untranslated regions (5′UTRs and 3′UTRs) of expressed sequences. Alignment of the whole genome and expressed sequence tags (ESTs) of the model plant Arabidopsis thaliana has identified introns residing in both coding and non-coding regions of the genome. Results A bioinformatic analysis revealed some interesting observations: (1) the density of introns in 5′UTRs is similar to that in CDSs but much higher than that in 3′UTRs; (2) the 5′UTR introns are preferentially located close to the initiating ATG codon; (3) introns in the 5′UTRs are, on average, longer than introns in the CDSs and 3′UTRs; and (4) 5′UTR introns have a different nucleotide composition to that of CDs and 3′UTR introns. Furthermore, we show that the 5′UTR intron of the A. thaliana EFIα-A3 gene affects the gene expression and the size of the 5′UTR intron influences the level of gene expression. Conclusion Introns within the 5′UTR show specific features that distinguish them from introns that reside within the coding sequence and the 3′UTR. In the EFIα-A3 gene, the presence of a long intron in the 5′UTR is sufficient to enhance gene expression in plants in a size dependent manner.

Gene electrotransfer into murine skeletal muscle: A systematic analysis of parameters for long-term gene expression

Skeletal muscle is an attractive target tissue for delivery of therapeutic genes, since it is well vascularized, easily accessible, and has a high capacity for protein synthesis. For efficient transfection in skeletal muscle, several protocols have been described, including delivery of low voltage electric pulses and a combination of high and low voltage electric pulses. The aim of this study was to determine the influence of different parameters of electrotransfection on short-term and long-term transfection efficiency in murine skeletal muscle, and to evaluate histological changes in the treated tissue. Different parameters of electric pulses, different time lags between plasmid DNA injection and application of electric pulses, and different doses of plasmid DNA were tested for electrotransfection of tibialis cranialis muscle of C57BI/6 mice using DNA plasmid encoding green fluorescent protein (GFP). Transfection efficiency was assessed on frozen tissue sections one week after electrotransfection using a fluorescence microscope and also noninvasively, followed by an in vivo imaging system using a fluorescence stereo microscope over a period of several months. Histological changes in muscle were evaluated immediately or several months after electrotransfection by determining infiltration of inflammatory mononuclear cells and presence of necrotic muscle fibers. The most efficient electrotransfection into skeletal muscle of C57BI/6 mice in our experiments was achieved when one high voltage (HV) and four low voltage (LV) electric pulses were applied 5 seconds after the injection of 30 μg of plasmid DNA. This protocol resulted in the highest short-term as well as long-term transfection. The fluorescence intensity of the transfected area declined after 2-3 weeks, but GFP fluorescence was still detectable 18 months after electrotransfection. Extensive inflammatory mononuclear cell infiltration was observed immediately after the electrotransfection procedure using the described parameters, but no necrosis or late tissue damage was observed. This study showed that electric pulse parameters, time lag between the injection of DNA and application of electric pulses, and dose of plasmid DNA affected the duration of transgene expression in murine skeletal muscle. Therefore, transgene expression in muscle can be controlled by appropriate selection of electrotransfection protocol.