Tandem repeat sequence variation as causative cis-eQTLs for protein-coding gene expression variation: the case of CSTB.

Association studies have revealed expression quantitative trait loci (eQTLs) for a large number of genes. However, the causative variants that regulate gene expression levels are generally unknown. We hypothesized that copy-number variation of sequence repeats contribute to the expression variation of some genes. Our laboratory has previously identified that the rare expansion of a repeat c.-174CGGGGCGGGGCG in the promoter region of the CSTB gene causes a silencing of the gene, resulting in progressive myoclonus epilepsy. Here, we genotyped the repeat length and quantified CSTB expression by quantitative real-time polymerase chain reaction in 173 lymphoblastoid cell lines (LCLs) and fibroblast samples from the GenCord collection. The majority of alleles contain either two or three copies of this repeat. Independent analysis revealed that the c.-174CGGGGCGGGGCG repeat length is strongly associated with CSTB expression (P = 3.14 × 10(-11)) in LCLs only. Examination of both genotyped and imputed single-nucleotide polymorphisms (SNPs) within 2 Mb of CSTB revealed that the dodecamer repeat represents the strongest cis-eQTL for CSTB in LCLs. We conclude that the common two or three copy variation is likely the causative cis-eQTL for CSTB expression variation. More broadly, we propose that polymorphic tandem repeats may represent the causative variation of a fraction of cis-eQTLs in the genome.

MFG-E8/lactadherin regulates cyclins D1/D3 expression and enhances the tumorigenic potential of mammary epithelial cells.

Milk fat globule-EGF factor 8 (MFG-E8) is a glycoprotein highly expressed in breast cancer that contributes to tumor progression through largely undefined mechanisms. By analyzing publicly available gene expression profiles of breast carcinomas, we found that MFG-E8 is highly expressed in primary and metastatic breast carcinomas, associated with absent estrogen receptor expression. Immunohistochemistry analysis of breast cancer biopsies revealed that MFG-E8 is expressed on the cell membrane as well as in the cytoplasm and nucleus. We also show that increased expression of MFG-E8 in mammary carcinoma cells increases their tumorigenicity in immunodeficient mice, and conversely, its downregulation reduces their in vivo growth. Moreover, expression of MFG-E8 in immortalized mammary epithelial cells promotes their growth and branching in three-dimensional collagen matrices and induces the expression of cyclins D1/D3 and N-cadherin. A mutant protein unable to bind integrins can in part exert these effects, indicating that MFG-E8 function is only partially dependent on integrin activation. We conclude that MFG-E8-dependent signaling stimulates cell proliferation and the acquisition of mesenchymal properties and contributes to mammary carcinoma development.

Cardiac raptor ablation impairs adaptive hypertrophy, alters metabolic gene expression, and causes heart failure in mice.

Background- Cardiac hypertrophy involves growth responses to a variety of stimuli triggered by increased workload. It is an independent risk factor for heart failure and sudden death. Mammalian target of rapamycin (mTOR) plays a key role in cellular growth responses by integrating growth factor and energy status signals. It is found in 2 structurally and functionally distinct multiprotein complexes called mTOR complex (mTORC) 1 and mTORC2. The role of each of these branches of mTOR signaling in the adult heart is currently unknown. Methods and Results- We generated mice with deficient myocardial mTORC1 activity by targeted ablation of raptor, which encodes an essential component of mTORC1, during adulthood. At 3 weeks after the deletion, atrial and brain natriuretic peptides and β-myosin heavy chain were strongly induced, multiple genes involved in the regulation of energy metabolism were altered, but cardiac function was normal. Function deteriorated rapidly afterward, resulting in dilated cardiomyopathy and high mortality within 6 weeks. Aortic banding-induced pathological overload resulted in severe dilated cardiomyopathy already at 1 week without a prior phase of adaptive hypertrophy. The mechanism involved a lack of adaptive cardiomyocyte growth via blunted protein synthesis capacity, as supported by reduced phosphorylation of ribosomal S6 kinase 1 and 4E-binding protein 1. In addition, reduced mitochondrial content, a shift in metabolic substrate use, and increased apoptosis and autophagy were observed. Conclusions- Our results demonstrate an essential function for mTORC1 in the heart under physiological and pathological conditions and are relevant for the understanding of disease states in which the insulin/insulin-like growth factor signaling axis is affected such as diabetes mellitus and heart failure or after cancer therapy.

Comparative modular analysis of gene expression in vertebrate organs.

ABSTRACT: BACKGROUND: The degree of conservation of gene expression between homologous organs largely remains an open question. Several recent studies reported some evidence in favor of such conservation. Most studies compute organs' similarity across all orthologous genes, whereas the expression level of many genes are not informative about organ specificity. RESULTS: Here, we use a modularization algorithm to overcome this limitation through the identification of inter-species co-modules of organs and genes. We identify such co-modules using mouse and human microarray expression data. They are functionally coherent both in terms of genes and of organs from both organisms. We show that a large proportion of genes belonging to the same co-module are orthologous between mouse and human. Moreover, their zebrafish orthologs also tend to be expressed in the corresponding homologous organs. Notable exceptions to the general pattern of conservation are the testis and the olfactory bulb. Interestingly, some co-modules consist of single organs, while others combine several functionally related organs. For instance, amygdala, cerebral cortex, hypothalamus and spinal cord form a clearly discernible unit of expression, both in mouse and human. CONCLUSIONS: Our study provides a new framework for comparative analysis which will be applicable also to other sets of large-scale phenotypic data collected across different species.

MAR elements regulate the probability of epigenetic switching between active and inactive gene expression.

Gene expression often cycles between active and inactive states in eukaryotes, yielding variable or noisy gene expression in the short-term, while slow epigenetic changes may lead to silencing or variegated expression. Understanding how cells control these effects will be of paramount importance to construct biological systems with predictable behaviours. Here we find that a human matrix attachment region (MAR) genetic element controls the stability and heritability of gene expression in cell populations. Mathematical modeling indicated that the MAR controls the probability of long-term transitions between active and inactive expression, thus reducing silencing effects and increasing the reactivation of silent genes. Single-cell short-terms assays revealed persistent expression and reduced expression noise in MAR-driven genes, while stochastic burst of expression occurred without this genetic element. The MAR thus confers a more deterministic behavior to an otherwise stochastic process, providing a means towards more reliable expression of engineered genetic systems.

Hepatitis B Virus e Antigen Physically Associates With Receptor-Interacting Serine/Threonine Protein Kinase 2 and Regulates IL-6 Gene Expression.

We previously reported that hepatitis B virus (HBV) e antigen (HBeAg) inhibits production of interleukin 6 by suppressing NF-κB activation. NF-κB is known to be activated through receptor-interacting serine/threonine protein kinase 2 (RIPK2), and we examined the mechanisms of interleukin 6 regulation by HBeAg. HBeAg inhibits RIPK2 expression and interacts with RIPK2, which may represent 2 mechanisms through which HBeAg blocks nucleotide-binding oligomerization domain-containing protein 1 ligand-induced NF-κB activation in HepG2 cells. Our findings identified novel molecular mechanisms whereby HBeAg modulates intracellular signaling pathways by targeting RIPK2, supporting the concept that HBeAg could impair both innate and adaptive immune responses to promote chronic HBV infection.

ExpressionView--an interactive viewer for modules identified in gene expression data.

SUMMARY: ExpressionView is an R package that provides an interactive graphical environment to explore transcription modules identified in gene expression data. A sophisticated ordering algorithm is used to present the modules with the expression in a visually appealing layout that provides an intuitive summary of the results. From this overview, the user can select individual modules and access biologically relevant metadata associated with them. AVAILABILITY: http://www.unil.ch/cbg/ExpressionView. Screenshots, tutorials and sample data sets can be found on the ExpressionView web site.

MicroRNAs shape circadian hepatic gene expression on a transcriptome-wide scale.

A considerable proportion of mammalian gene expression undergoes circadian oscillations. Post-transcriptional mechanisms likely make important contributions to mRNA abundance rhythms. We have investigated how microRNAs (miRNAs) contribute to core clock and clock-controlled gene expression using mice in which miRNA biogenesis can be inactivated in the liver. While the hepatic core clock was surprisingly resilient to miRNA loss, whole transcriptome sequencing uncovered widespread effects on clock output gene expression. Cyclic transcription paired with miRNA-mediated regulation was thus identified as a frequent phenomenon that affected up to 30% of the rhythmic transcriptome and served to post-transcriptionally adjust the phases and amplitudes of rhythmic mRNA accumulation. However, only few mRNA rhythms were actually generated by miRNAs. Overall, our study suggests that miRNAs function to adapt clock-driven gene expression to tissue-specific requirements. Finally, we pinpoint several miRNAs predicted to act as modulators of rhythmic transcripts, and identify rhythmic pathways particularly prone to miRNA regulation.DOI: http://dx.doi.org/10.7554/eLife.02510.001.

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

Les cancer-testis antigènes appartiennent à la famille des antigènes tumoraux spécifiques. Ils ont montré un pouvoir immunogène chez les patients porteurs de différents cancers. En effet, ils stimulent sélectivement les lymphocytes cytotoxiques, et leur expression spécifique dans les tissus tumoraux en fait une cible idéale pour une vaccination antitumorale. Le but de cette étude est d'identifier l'expression de certains de ces antigènes, d'analyser leur valeur pronostique et de déterminer la meilleure cible antigénique pour permettre une immunothérapie spécifique dans les carcinomes épidermoïdes des voies aérodigestives supérieures. Le profil et le taux d'expression de 12 cancer-testis antigènes (MAGE-A1, MAGE-A3, MAGE-A4, MAGEA10, MAGE-C2, NY-ESO-1, LAGE-1, SSX-2, SSX-4, BAGE, GAGE-1/2, GAGE-3/4) et de 3 autres antigènes tumoraux spécifiques (PRAME, HERV-K-MEL, NA-17A) ont été évalués par RT-PCR sur 57 échantillons de cancers ORL primaires. Les paramètres tumoraux et cliniques ont été prospectivement collectés afin de corréler ces données avec le résultat de nos investigations immunobiologiques. Quatre-vingt-huit pour cent des tumeurs expriment au moins 1 antigène. Une co-expression de 3 gènes ou plus est détectée chez 59% des patients. MAGE-A4 (60%), MAGE-A3 (51%), PRAME (49%) et HERV-K-MEL (42%) sont les gènes le plus fréquemment exprimés. Ils sont totalement absents des muqueuses saines avoisinantes. La présence de MAGE-A et NY-ESO-1 à la surface des cellules a été vérifiée par immunohistochimie. Nos analyses statistiques ont permis d'identifier une diminution de la survie liée au cancer chez les patients porteurs d'une tumeur exprimant de multiples cancer-testis antigènes et notamment MAGE-A4 dont l'expression indépendante d'autres éléments cliniques s'associe statistiquement à un taux de survie diminué. Nos résultats ont permis d'identifier un rôle pronostique de l'expression des gènes associés aux tumeurs dont l'expression est apparemment liée à un phénotype de malignité plus élevé. Cette constatation, corroborée par l'identification parallèle d'un infiltrat lymphocytaire spécifique confirme l'utilité potentielle de certains cancer-testis antigènes comme cible pour une immunothérapie ciblée dans les carcinomes des voies aérodigestives supérieures

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

The CD8 T cell response generatedby gene-based vaccines is importantfor protective immunity againstmany infectious diseases but its complexityis incompletely understood.Here, we report that different vaccinesencoding HIV Env elicit qualitativelydistinct CD8 T cells that wereidentified by patterns of gene expressionin individual cells. Three alternativeprime-boost vector combinationsstimulated antigen-specific CD8 Tcell populations of similar magnitudeand function by intracellular cytokinestaining; however, single cell geneexpression profiling enabled the discriminationof distinct CM and EMCD8 cells elicited by the three vaccines.Two previously unrecognizedCD8 T cell subsets have been definedby their coexpression of Eomes,Cxcr3 and Ccr7; or Klrk1, Klrg1 andCcr5 in CM and EM cells respectively.

eVOC: a controlled vocabulary for unifying gene expression data.

Expression data contribute significantly to the biological value of the sequenced human genome, providing extensive information about gene structure and the pattern of gene expression. ESTs, together with SAGE libraries and microarray experiment information, provide a broad and rich view of the transcriptome. However, it is difficult to perform large-scale expression mining of the data generated by these diverse experimental approaches. Not only is the data stored in disparate locations, but there is frequent ambiguity in the meaning of terms used to describe the source of the material used in the experiment. Untangling semantic differences between the data provided by different resources is therefore largely reliant on the domain knowledge of a human expert. We present here eVOC, a system which associates labelled target cDNAs for microarray experiments, or cDNA libraries and their associated transcripts with controlled terms in a set of hierarchical vocabularies. eVOC consists of four orthogonal controlled vocabularies suitable for describing the domains of human gene expression data including Anatomical System, Cell Type, Pathology and Developmental Stage. We have curated and annotated 7016 cDNA libraries represented in dbEST, as well as 104 SAGE libraries,with expression information,and provide this as an integrated, public resource that allows the linking of transcripts and libraries with expression terms. Both the vocabularies and the vocabulary-annotated libraries can be retrieved from http://www.sanbi.ac.za/evoc/. Several groups are involved in developing this resource with the aim of unifying transcript expression information.

Evolution under monogamy feminizes gene expression in Drosophila melanogaster.

Many genes have evolved sexually dimorphic expression as a consequence of divergent selection on males and females. However, because the sexes share a genome, the extent to which evolution can shape gene expression independently in each sex is controversial. Here, we use experimental evolution to reveal suboptimal sex-specific expression for much of the genome. By enforcing a monogamous mating system in populations of Drosophila melanogaster for over 100 generations, we eliminated major components of selection on males: female choice and male-male competition. If gene expression is subject to sexually antagonistic selection, relaxed selection on males should cause evolution towards female optima. Monogamous males and females show this pattern of feminization in both the whole-body and head transcriptomes. Genes with male-biased expression patterns evolved decreased expression under monogamy, while genes with female-biased expression evolved increased expression, relative to polygamous populations. Our results demonstrate persistent and widespread evolutionary tension between male and female adaptation.

Stem cell-related "self-renewal" signature and high epidermal growth factor receptor expression associated with resistance to concomitant chemoradiotherapy in glioblastoma.

PURPOSE: Glioblastomas are notorious for resistance to therapy, which has been attributed to DNA-repair proficiency, a multitude of deregulated molecular pathways, and, more recently, to the particular biologic behavior of tumor stem-like cells. Here, we aimed to identify molecular profiles specific for treatment resistance to the current standard of care of concomitant chemoradiotherapy with the alkylating agent temozolomide. PATIENTS AND METHODS: Gene expression profiles of 80 glioblastomas were interrogated for associations with resistance to therapy. Patients were treated within clinical trials testing the addition of concomitant and adjuvant temozolomide to radiotherapy. RESULTS: An expression signature dominated by HOX genes, which comprises Prominin-1 (CD133), emerged as a predictor for poor survival in patients treated with concomitant chemoradiotherapy (n = 42; hazard ratio = 2.69; 95% CI, 1.38 to 5.26; P = .004). This association could be validated in an independent data set. Provocatively, the HOX cluster was reminiscent of a "self-renewal" signature (P = .008; Gene Set Enrichment Analysis) recently characterized in a mouse leukemia model. The HOX signature and EGFR expression were independent prognostic factors in multivariate analysis, adjusted for the O-6-methylguanine-DNA methyltransferase (MGMT) methylation status, a known predictive factor for benefit from temozolomide, and age. Better outcome was associated with gene clusters characterizing features of tumor-host interaction including tumor vascularization and cell adhesion, and innate immune response. CONCLUSION: This study provides first clinical evidence for the implication of a "glioma stem cell" or "self-renewal" phenotype in treatment resistance of glioblastoma. Biologic mechanisms identified here to be relevant for resistance will guide future targeted therapies and respective marker development for individualized treatment and patient selection.

Analysis of gene expression patterns in animals

During my PhD, my aim was to provide new tools to increase our capacity to analyse gene expression patterns, and to study on a large-scale basis the evolution of gene expression in animals. Gene expression patterns (when and where a gene is expressed) are a key feature in understanding gene function, notably in development. It appears clear now that the evolution of developmental processes and of phenotypes is shaped both by evolution at the coding sequence level, and at the gene expression level.Studying gene expression evolution in animals, with complex expression patterns over tissues and developmental time, is still challenging. No tools are available to routinely compare expression patterns between different species, with precision, and on a large-scale basis. Studies on gene expression evolution are therefore performed only on small genes datasets, or using imprecise descriptions of expression patterns.The aim of my PhD was thus to develop and use novel bioinformatics resources, to study the evolution of gene expression. To this end, I developed the database Bgee (Base for Gene Expression Evolution). The approach of Bgee is to transform heterogeneous expression data (ESTs, microarrays, and in-situ hybridizations) into present/absent calls, and to annotate them to standard representations of anatomy and development of different species (anatomical ontologies). An extensive mapping between anatomies of species is then developed based on hypothesis of homology. These precise annotations to anatomies, and this extensive mapping between species, are the major assets of Bgee, and have required the involvement of many co-workers over the years. My main personal contribution is the development and the management of both the Bgee database and the web-application.Bgee is now on its ninth release, and includes an important gene expression dataset for 5 species (human, mouse, drosophila, zebrafish, Xenopus), with the most data from mouse, human and zebrafish. Using these three species, I have conducted an analysis of gene expression evolution after duplication in vertebrates.Gene duplication is thought to be a major source of novelty in evolution, and to participate to speciation. It has been suggested that the evolution of gene expression patterns might participate in the retention of duplicate genes. I performed a large-scale comparison of expression patterns of hundreds of duplicated genes to their singleton ortholog in an outgroup, including both small and large-scale duplicates, in three vertebrate species (human, mouse and zebrafish), and using highly accurate descriptions of expression patterns. My results showed unexpectedly high rates of de novo acquisition of expression domains after duplication (neofunctionalization), at least as high or higher than rates of partitioning of expression domains (subfunctionalization). I found differences in the evolution of expression of small- and large-scale duplicates, with small-scale duplicates more prone to neofunctionalization. Duplicates with neofunctionalization seemed to evolve under more relaxed selective pressure on the coding sequence. Finally, even with abundant and precise expression data, the majority fate I recovered was neither neo- nor subfunctionalization of expression domains, suggesting a major role for other mechanisms in duplicate gene retention.

Plants respond to pathogen infection by enhancing the antifungal gene expression of root-associated bacteria.

Plant health and fitness widely depend on interactions with soil microorganisms. Some bacteria such as pseudomonads can inhibit pathogens by producing antibiotics, and controlling these bacteria could help improve plant fitness. In the present study, we tested whether plants induce changes in the antifungal activity of root-associated bacteria as a response to root pathogens. We grew barley plants in a split-root system with one side of the root system challenged by the pathogen Pythium ultimum and the other side inoculated with the biocontrol strain Pseudomonas fluorescens CHA0. We used reporter genes to follow the expression of ribosomal RNA indicative of the metabolic state and of the gene phlA, required for production of 2,4-diacetylphloroglucinol, a key component of antifungal activity. Infection increased the expression of the antifungal gene phlA. No contact with the pathogen was required, indicating that barley influenced gene expression by the bacteria in a systemic way. This effect relied on increased exudation of diffusible molecules increasing phlA expression, suggesting that communication with rhizosphere bacteria is part of the pathogen response of plants. Tripartite interactions among plants, pathogens, and bacteria appear as a novel determinant of plant response to root pathogens.