Gene Expression in Brain: A Window on Ethanol Dependence, Neuroadaptation, and Preference

This article represents the proceedings of a symposium at the 2002 joint RSA/ISBRA Conference in San Francisco, California. The organizer was Paula L. Hoffman and the co-chairs were Paula L. Hoffman and Michael Miles. The presentations were (1) Introduction and overview of the use of DNA microarrays, by Michael Miles; (2) DNA microarray analysis of gene expression in brains of P and NP rats, by Howard J. Edenberg; (3) Gene expression patterns in brain regions of AA and ANA rats, by Wolfgang Sommer; (4) Patterns of gene expression in brains of selected lines of mice that differ in ethanol tolerance, by Boris Tabakoff; (5) Gene expression profiling related to initial sensitivity and tolerance in gamma-protein kinase C mutants, by Jeanne Wehner; and (6) Gene expression patterns in human alcoholic brain: from microarrays to protein profiles, by Joanne Lewohl.

Model-based clustering in gene expression microarrays: An application to breast cancer data

In microarray studies, the application of clustering techniques is often used to derive meaningful insights into the data. In the past, hierarchical methods have been the primary clustering tool employed to perform this task. The hierarchical algorithms have been mainly applied heuristically to these cluster analysis problems. Further, a major limitation of these methods is their inability to determine the number of clusters. Thus there is a need for a model-based approach to these. clustering problems. To this end, McLachlan et al. [7] developed a mixture model-based algorithm (EMMIX-GENE) for the clustering of tissue samples. To further investigate the EMMIX-GENE procedure as a model-based -approach, we present a case study involving the application of EMMIX-GENE to the breast cancer data as studied recently in van 't Veer et al. [10]. Our analysis considers the problem of clustering the tissue samples on the basis of the genes which is a non-standard problem because the number of genes greatly exceed the number of tissue samples. We demonstrate how EMMIX-GENE can be useful in reducing the initial set of genes down to a more computationally manageable size. The results from this analysis also emphasise the difficulty associated with the task of separating two tissue groups on the basis of a particular subset of genes. These results also shed light on why supervised methods have such a high misallocation error rate for the breast cancer data.

Influence of ageing, heat shock treatment and in vivo total antioxidant status on gene-expression profile and protein synthesis in human peripheral lymphocytes

Ageing results in a progressive, intrinsic and generalised imbalance of the control of regulatory systems. A key manifestation of this complex biological process includes the attenuation of the universal stress response. Here we provide the first global assessment of the ageing process as it affects the heat shock response, utilising human peripheral lymphocytes and cDNA microarray analysis. The genomic approach employed in our preliminary study was supplemented with a proteomic approach. In addition, the current study correlates the in vivo total antioxidant status with the age-related differential gene expression as well as the translational kinetics of heat shock proteins (hsps). Most of the genes encoding stress response proteins on the 4224 element microarray used in this study were significantly elevated after heat shock treatment of lymphocytes obtained from both young and old individuals albeit to a greater extent in the young. Cell signaling and signal transduction genes as well as some oxidoreductases showed varied response. Results from translational kinetics of induction of major hsps, from 0 to 24 It recovery period were broadly consistent with the differential expression of HSC 70 and HSP 40 genes. Total antioxidant levels in plasma from old individuals were found to be significantly lower by comparison with young, in agreement with the widely acknowledged role of oxidant homeostasis in the ageing process. (C) 2002 Elsevier Science Ireland Ltd. All rights reserved.

Expression analysis of four Pinus radiata male cone promoters in the heterologous host Arabidopsis

Four male cone-specific promoters were isolated from the genome of Pinus radiata D. Don, fused to the beta-glucuronidase (GUS) reporter gene and analysed in the heterologous host Arabidopsis thaliana (L.) Heynh. The temporal and spatial activities of the promoters PrCHS1, PrLTP2, PrMC2 and PrMALE1 during seven anther developmental stages are described in detail. The two promoters PrMC2 and PrMALE1 confer an identical GUS expression pattern on Arabidopsis anthers. DNA sequence analysis of the PrMC2 and PrMALE1 promoters revealed an 88% sequence identity over 276 bp and divergence further upstream (

Characterization of a polyubiquitin gene in T. thermophila and of ubiquitin gene expression during sexual reproduction and under stress conditions

A 5-unit polyubiquitin gene, TTU3, was isolated from a T. thermophila genomic library and sequenced. This gene presents an extra triplet coding for Phe, a AGAGA motif and a putative HSE element in its 5'-non-coding region. The ubiquitin gene expression in this ciliate was investigated by Northern blot hybridization in conjugating cells or cells under stress conditions. Exponentially growing cells express two ubiquitin mRNAs of 0.75 and 1.8 kb and a new species of 1.4 kb is induced under hyperthermic stress. During sexual reproduction of the cells (conjugation) the 1.8-kb mRNA is still transcribed whereas the steady-state population of the 0.75 mRNA transcripts is strongly diminished. Southern blot analysis suggests that ubiquitin in T. thermophila constitutes a large family of about ten members.

Studies on BolA and ribonuclease R: two important factors in the control of bacterial gene expression

Dissertation presented to obtain the Ph.D degree in Biology

Gene Expression Profiling and Association Studies Implicate the Neuregulin Signaling Pathway in Behçet's Disease Susceptibility

Behçet's disease (BD) is a complex disease with genetic and environmental risk factors implicated in its etiology; however, its pathophysiology is poorly understood. To decipher BD's genetic underpinnings, we combined gene expression profiling with pathway analysis and association studies. We compared the gene expression profiles in peripheral blood mononuclear cells (PBMCs) of 15 patients and 14 matched controls using Affymetrix microarrays and found that the neuregulin signaling pathway was over-represented among the differentially expressed genes. The Epiregulin (EREG), Amphiregulin (AREG), and Neuregulin-1 (NRG1) genes of this pathway stand out as they are also among the top differentially expressed genes. Twelve haplotype tagging SNPs at the EREG-AREG locus and 15 SNPs in NRG1 found associated in at least one published BD genome-wide association study were tested for association with BD in a dataset of 976 Iranian patients and 839 controls. We found a novel association with BD for the rs6845297 SNP located downstream of EREG, and replicated three associations at NRG1 (rs4489285, rs383632, and rs1462891). Multifactor dimensionality reduction analysis indicated the existence of epistatic interactions between EREG and NRG1 variants. EREG-AREG and NRG1, which are members of the epidermal growth factor (EGF) family, seem to modulate BD susceptibility through main effects and gene–gene interactions. These association findings support a role for the EGF/ErbB signaling pathway inBD pathogenesis that warrants further investigation and highlight the importance of combining genetic and genomic approaches to dissect the genetic architecture of complex diseases.

Controlling gene expression in enterobacteriaceae: studies on sRNAs and strategies for synthetic biology

Dissertation presented to obtain the Master Degree in Molecular, Genetics and Biomedicine

IFN stimulated gene expression in the liver is a better predictor of treatment response in chronic hepatitis c than the IL28b (IFN lambda 3) genotype

Background: Therapy of chronic hepatitis C (CHC) with pegIFNa/ribavirin achieves sustained virologic response (SVR) in ~55%. Pre-activation of the endogenous interferon system in the liver is associated non-response (NR). Recently, genome-wide association studies described associations of allelic variants near the IL28B (IFNλ3) gene with treatment response and with spontaneous clearance of the virus. We investigated if the IL28B genotype determines the constitutive expression of IFN stimulated genes (ISGs) in the liver of patients with CHC. Methods: We genotyped 93 patients with CHC for 3 IL28B single nucleotide polymorphisms (SNPs, rs12979860, rs8099917, rs12980275), extracted RNA from their liver biopsies and quantified the expression of IL28B and of 8 previously identified classifier genes which discriminate between SVR and NR (IFI44L, RSAD2, ISG15, IFI22, LAMP3, OAS3, LGALS3BP and HTATIP2). Decision tree ensembles in the form of a random forest classifier were used to calculate the relative predictive power of these different variables in a multivariate analysis. Results: The minor IL28B allele (bad risk for treatment response) was significantly associated with increased expression of ISGs, and, unexpectedly, with decreased expression of IL28B. Stratification of the patients into SVR and NR revealed that ISG expression was conditionally independent from the IL28B genotype, i.e. there was an increased expression of ISGs in NR compared to SVR irrespective of the IL28B genotype. The random forest feature score (RFFS) identified IFI27 (RFFS = 2.93), RSAD2 (1.88) and HTATIP2 (1.50) expression and the HCV genotype (1.62) as the strongest predictors of treatment response. ROC curves of the IL28B SNPs showed an AUC of 0.66 with an error rate (ERR) of 0.38. A classifier with the 3 best classifying genes showed an excellent test performance with an AUC of 0.94 and ERR of 0.15. The addition of IL28B genotype information did not improve the predictive power of the 3-gene classifier. Conclusions: IL28B genotype and hepatic ISG expression are conditionally independent predictors of treatment response in CHC. There is no direct link between altered IFNλ3 expression and pre-activation of the endogenous system in the liver. Hepatic ISG expression is by far the better predictor for treatment response than IL28B genotype.

A simple genetic basis for complex social behaviour mediates widespread gene expression differences.

A remarkable social polymorphism is controlled by a single Mendelian factor in the fire ant Solenopsis invicta. A genomic element marked by the gene Gp-9 determines whether workers tolerate one or many fertile queens in their colony. Gp-9 was recently shown to be part of a supergene with two nonrecombining variants, SB and Sb. SB/SB and SB/Sb queens differ in how they initiate new colonies, and in many physiological traits, for example odour and maturation rate. To understand how a single genetic element can affect all these traits, we used a microarray to compare gene expression patterns between SB/SB and SB/Sb queens of three different age classes: 1-day-old unmated queens, 11-day-old unmated queens and mated, fully reproductive queens collected from mature field colonies. The number of genes that were differentially expressed between SB/SB and SB/Sb queens of the same age class was smallest in 1-day-old queens, maximal in 11-day-old queens and intermediate in reproductive queens. Gene ontology analysis showed that SB/SB queens upregulate reproductive genes faster than SB/Sb queens. For all age classes, genes inside the supergene were overrepresented among the differentially expressed genes. Consistent with the hypothesized greater number of transposons in the Sb supergene, 13 transposon genes were upregulated in SB/Sb queens. Viral genes were also upregulated in SB/Sb mature queens, consistent with the known greater parasite load in colonies headed by SB/Sb queens compared with colonies headed by SB/SB queens. Eighteen differentially expressed genes between reproductive queens were involved in chemical signalling. Our results suggest that many genes in the supergene are involved in regulating social organization and queen phenotypes in fire ants.

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.

EuroDia: a beta-cell gene expression resource.

Type 2 diabetes mellitus (T2DM) is a major disease affecting nearly 280 million people worldwide. Whilst the pathophysiological mechanisms leading to disease are poorly understood, dysfunction of the insulin-producing pancreatic beta-cells is key event for disease development. Monitoring the gene expression profiles of pancreatic beta-cells under several genetic or chemical perturbations has shed light on genes and pathways involved in T2DM. The EuroDia database has been established to build a unique collection of gene expression measurements performed on beta-cells of three organisms, namely human, mouse and rat. The Gene Expression Data Analysis Interface (GEDAI) has been developed to support this database. The quality of each dataset is assessed by a series of quality control procedures to detect putative hybridization outliers. The system integrates a web interface to several standard analysis functions from R/Bioconductor to identify differentially expressed genes and pathways. It also allows the combination of multiple experiments performed on different array platforms of the same technology. The design of this system enables each user to rapidly design a custom analysis pipeline and thus produce their own list of genes and pathways. Raw and normalized data can be downloaded for each experiment. The flexible engine of this database (GEDAI) is currently used to handle gene expression data from several laboratory-run projects dealing with different organisms and platforms. Database URL: http://eurodia.vital-it.ch.

Transcriptional regulation of CD4 gene expression by T cell factor-1/beta-catenin pathway.

By interacting with MHC class II molecules, CD4 facilitates lineage development as well as activation of Th cells. Expression of physiological levels of CD4 requires a proximal CD4 enhancer to stimulate basic CD4 promoter activity. T cell factor (TCF)-1/beta-catenin pathway has previously been shown to regulate thymocyte survival via up-regulating antiapoptotic molecule Bcl-xL. By both loss and gain of function studies, in this study we show additional function of TCF-1/beta-catenin pathway in the regulation of CD4 expression in vivo. Mice deficient in TCF-1 displayed significantly reduced protein and mRNA levels of CD4 in CD4+ CD8+ double-positive (DP) thymocytes. A transgene encoding Bcl-2 restored survival but not CD4 levels of TCF-1(-/-) DP cells. Thus, TCF-1-regulated survival and CD4 expression are two separate events. In contrast, CD4 levels were restored on DP TCF-1(-/-) cells by transgenic expression of a wild-type TCF-1, but not a truncated TCF-1 that lacks a domain required for interacting with beta-catenin. Furthermore, forced expression of a stabilized beta-catenin, a coactivator of TCF-1, resulted in up-regulation of CD4. TCF-1 or stabilized beta-catenin greatly stimulated activity of a CD4 reporter gene driven by a basic CD4 promoter and the CD4 enhancer. However, mutation of a potential TCF binding site located within the enhancer abrogated TCF-1 and beta-catenin-mediated activation of CD4 reporter. Finally, recruitment of TCF-1 to CD4 enhancer was detected in wild-type but not TCF-1 null mice by chromatin-immunoprecipitation analysis. Thus, our results demonstrated that TCF/beta-catenin pathway enhances CD4 expression in vivo by recruiting TCF-1 to stimulate CD4 enhancer activity.

Pattern and clinical significance of cancer-testis gene expression in head and neck squamous cell carcinoma.

Cancer-testis (CT) antigens comprise families of tumor-associated antigens that are immunogenic in patients with various cancers. Their restricted expression makes them attractive targets for immunotherapy. The aim of this study was to determine the expression of several CT genes and evaluate their prognostic value in head and neck squamous cell carcinoma (HNSCC). The pattern and level of expression of 12 CT genes (MAGE-A1, MAGE-A3, MAGE-A4, MAGE-A10, MAGE-C2, NY-ESO-1, LAGE-1, SSX-2, SSX-4, BAGE, GAGE-1/2, GAGE-3/4) and the tumor-associated antigen encoding genes PRAME, HERV-K-MEL, and NA-17A were evaluated by RT-PCR in a panel of 57 primary HNSCC. Over 80% of the tumors expressed at least 1 CT gene. Coexpression of three or more genes was detected in 59% of the patients. MAGE-A4 (60%), MAGE-A3 (51%), PRAME (49%) and HERV-K-MEL (42%) were the most frequently expressed genes. Overall, the pattern of expression of CT genes indicated a coordinate regulation; however there was no correlation between expression of MAGE-A3/A4 and BORIS, a gene whose product has been implicated in CT gene activation. The presence of MAGE-A and NY-ESO-1 proteins was verified by immunohistochemistry. Analysis of the correlation between mRNA expression of CT genes with clinico-pathological characteristics and clinical outcome revealed that patients with tumors positive for MAGE-A4 or multiple CT gene expression had a poorer overall survival. Furthermore, MAGE-A4 mRNA positivity was prognostic of poor outcome independent of clinical parameters. These findings indicate that expression of CT genes is associated with a more malignant phenotype and suggest their usefulness as prognostic markers in HNSCC.

Liver hyperplasia and paradoxical regulation of glycogen metabolism and glucose-sensitive gene expression in GLUT2-null hepatocytes. Further evidence for the existence of a membrane-based glucose release pathway.

We investigated the impact of GLUT2 gene inactivation on the regulation of hepatic glucose metabolism during the fed to fast transition. In control and GLUT2-null mice, fasting was accompanied by a approximately 10-fold increase in plasma glucagon to insulin ratio, a similar activation of liver glycogen phosphorylase and inhibition of glycogen synthase and the same elevation in phosphoenolpyruvate carboxykinase and glucose-6-phosphatase mRNAs. In GLUT2-null mice, mobilization of glycogen stores was, however, strongly impaired. This was correlated with glucose-6-phosphate (G6P) levels, which remained at the fed values, indicating an important allosteric stimulation of glycogen synthase by G6P. These G6P levels were also accompanied by a paradoxical elevation of the mRNAs for L-pyruvate kinase. Re-expression of GLUT2 in liver corrected the abnormal regulation of glycogen and L-pyruvate kinase gene expression. Interestingly, GLUT2-null livers were hyperplasic, as revealed by a 40% increase in liver mass and 30% increase in liver DNA content. Together, these data indicate that in the absence of GLUT2, the G6P levels cannot decrease during a fasting period. This may be due to neosynthesized glucose entering the cytosol, being unable to diffuse into the extracellular space, and being phosphorylated back to G6P. Because hepatic glucose production is nevertheless quantitatively normal, glucose produced in the endoplasmic reticulum may also be exported out of the cell through an alternative, membrane traffic-based pathway, as previously reported (Guillam, M.-T., Burcelin, R., and Thorens, B. (1998) Proc. Natl. Acad. Sci. U. S. A. 95, 12317-12321). Therefore, in fasting, GLUT2 is not required for quantitative normal glucose output but is necessary to equilibrate cytosolic glucose with the extracellular space. In the absence of this equilibration, the control of hepatic glucose metabolism by G6P is dominant over that by plasma hormone concentrations.