MADAP, a flexible clustering tool for the interpretation of one-dimensional genome annotation data.

A recurring task in the analysis of mass genome annotation data from high-throughput technologies is the identification of peaks or clusters in a noisy signal profile. Examples of such applications are the definition of promoters on the basis of transcription start site profiles, the mapping of transcription factor binding sites based on ChIP-chip data and the identification of quantitative trait loci (QTL) from whole genome SNP profiles. Input to such an analysis is a set of genome coordinates associated with counts or intensities. The output consists of a discrete number of peaks with respective volumes, extensions and center positions. We have developed for this purpose a flexible one-dimensional clustering tool, called MADAP, which we make available as a web server and as standalone program. A set of parameters enables the user to customize the procedure to a specific problem. The web server, which returns results in textual and graphical form, is useful for small to medium-scale applications, as well as for evaluation and parameter tuning in view of large-scale applications, requiring a local installation. The program written in C++ can be freely downloaded from ftp://ftp.epd.unil.ch/pub/software/unix/madap. The MADAP web server can be accessed at http://www.isrec.isb-sib.ch/madap/.

The alarmin HMGB1 enhances CXCR4-induced activities by binding CXCL12

A variety of chemokines and inflammatory molecules are concomitantly produced at target sites of leukocyte trafficking and homing, accounting for the complex cellular responses occurring in homeostasis and inflammation. The chemokine CXCL12 plays an essential and unique role in homeostatic regulation of leukocyte traffic and tissue regeneration. The chromatin protein HMGB1 is released by dying and distressed cells, and acts as a Damage Associated Molecular Pattern or alarmin, promoting cell migration towards the site of tissue damage. We show here that HMGB1 synergises with CXCL12 by forming a heterocomplex that we characterized by NMR chemical shift mapping. The heterocomplex enhances CXCR4-induced responses on cells of the immune system, acting exclusively through the CXCL12 receptor CXCR4, and not through the HMGB1 receptors RAGE, TLR2 and TLR4. FRET analysis show that CXCL12 and CXCL12+HMGB1 promote a different conformational change in the homodimer CXCR4. The enhancement induced by HMGB1 on CXCL12-induced migration is selective, since little changes in migration of neutrophils and PreB 300.19-CCR2+ or -CCR7+ are observed towards CXCL8 and CCR2 or CCR7 agonists. HMGB1 also promotes CXCL 12 release, which is ultimately responsible for the chemoattractant activities of HMGB1. This study highlights the role of HMGB1 in promoting CXCL12-dependent cell migration, and suggests a cooperative role of these two molecules in tissue repair as well as in pathological conditions, such as rheumatoid arthritis.

Functional interaction between the estrogen receptor and CTF1: analysis of the vitellogenin gene B1 promoter in yeast.

Eukaryotic gene expression depends on a complex interplay between the transcriptional apparatus and chromatin structure. We report here a yeast model system for investigating the functional interaction between the human estrogen receptor (hER) and CTF1, a member of the CTF/NFI transcription factor family. We show that a CTF1-fusion protein and the hER transactivate a synthetic promoter in yeast in a synergistic manner. This interaction requires the proline-rich transactivation domain of CTF1. When the natural estrogen-dependent vitellogenin B1 promoter is tested in yeast, CTF1 and CTF1-fusion proteins are unable to activate transcription, and no synergy is observed between hER, which activates the B1 promoter, and these factors. Chromatin structure analysis on this promoter reveals positioned nucleosomes at -430 to -270 (+/-20 bp) and at -270 to - 100 (+/-20 bp) relative to the start site of transcription. The positions of the nucleosomes remain unchanged upon hormone-dependent transcriptional activation of the promoter, and the more proximal nucleosome appears to mask the CTF/NFI site located at - 101 to -114. We conclude that a functional interaction of hER with the estrogen response element located upstream of a basal promoter occurs in yeast despite the nucleosomal organization of this promoter, whereas the interaction of CTF1 with its target site is apparently precluded by a nucleosome.

Seventy-five genetic loci influencing the human red blood cell.

Anaemia is a chief determinant of global ill health, contributing to cognitive impairment, growth retardation and impaired physical capacity. To understand further the genetic factors influencing red blood cells, we carried out a genome-wide association study of haemoglobin concentration and related parameters in up to 135,367 individuals. Here we identify 75 independent genetic loci associated with one or more red blood cell phenotypes at P < 10(-8), which together explain 4-9% of the phenotypic variance per trait. Using expression quantitative trait loci and bioinformatic strategies, we identify 121 candidate genes enriched in functions relevant to red blood cell biology. The candidate genes are expressed preferentially in red blood cell precursors, and 43 have haematopoietic phenotypes in Mus musculus or Drosophila melanogaster. Through open-chromatin and coding-variant analyses we identify potential causal genetic variants at 41 loci. Our findings provide extensive new insights into genetic mechanisms and biological pathways controlling red blood cell formation and function.

PPARalpha governs glycerol metabolism.

Glycerol, a product of adipose tissue lipolysis, is an important substrate for hepatic glucose synthesis. However, little is known about the regulation of hepatic glycerol metabolism. Here we show that several genes involved in the hepatic metabolism of glycerol, i.e., cytosolic and mitochondrial glycerol 3-phosphate dehydrogenase (GPDH), glycerol kinase, and glycerol transporters aquaporin 3 and 9, are upregulated by fasting in wild-type mice but not in mice lacking PPARalpha. Furthermore, expression of these genes was induced by the PPARalpha agonist Wy14643 in wild-type but not PPARalpha-null mice. In adipocytes, which express high levels of PPARgamma, expression of cytosolic GPDH was enhanced by PPARgamma and beta/delta agonists, while expression was decreased in PPARgamma(+/-) and PPARbeta/delta(-/-) mice. Transactivation, gel shift, and chromatin immunoprecipitation experiments demonstrated that cytosolic GPDH is a direct PPAR target gene. In line with a stimulating role of PPARalpha in hepatic glycerol utilization, administration of synthetic PPARalpha agonists in mice and humans decreased plasma glycerol. Finally, hepatic glucose production was decreased in PPARalpha-null mice simultaneously fasted and exposed to Wy14643, suggesting that the stimulatory effect of PPARalpha on gluconeogenic gene expression was translated at the functional level. Overall, these data indicate that PPARalpha directly governs glycerol metabolism in liver, whereas PPARgamma regulates glycerol metabolism in adipose tissue.

In vivo activation of PPAR target genes by RXR homodimers.

The ability of a retinoid X receptor (RXR) to heterodimerize with many nuclear receptors, including LXR, PPAR, NGF1B and RAR, underscores its pivotal role within the nuclear receptor superfamily. Among these heterodimers, PPAR:RXR is considered an important signalling mediator of both PPAR ligands, such as fatty acids, and 9-cis retinoic acid (9-cis RA), an RXR ligand. In contrast, the existence of an RXR/9-cis RA signalling pathway independent of PPAR or any other dimerization partner remains disputed. Using in vivo chromatin immunoprecipitation, we now show that RXR homodimers can selectively bind to functional PPREs and induce transactivation. At the molecular level, this pathway requires stabilization of the homodimer-DNA complexes through ligand-dependent interaction with the coactivator SRC1 or TIF2. This pathway operates both in the absence and in the presence of PPAR, as assessed in cells carrying inactivating mutations in PPAR genes and in wild-type cells. In addition, this signalling pathway via PPREs is fully functional and can rescue the severe hypothermia phenotype observed in fasted PPARalpha-/- mice. These observations have important pharmacological implications for the development of new rexinoid-based treatments.

Models that include supercoiling of topological domains reproduce several known features of interphase chromosomes.

Understanding the structure of interphase chromosomes is essential to elucidate regulatory mechanisms of gene expression. During recent years, high-throughput DNA sequencing expanded the power of chromosome conformation capture (3C) methods that provide information about reciprocal spatial proximity of chromosomal loci. Since 2012, it is known that entire chromatin in interphase chromosomes is organized into regions with strongly increased frequency of internal contacts. These regions, with the average size of ∼1 Mb, were named topological domains. More recent studies demonstrated presence of unconstrained supercoiling in interphase chromosomes. Using Brownian dynamics simulations, we show here that by including supercoiling into models of topological domains one can reproduce and thus provide possible explanations of several experimentally observed characteristics of interphase chromosomes, such as their complex contact maps.

Scaffold attachment factor B1 directly interacts with nuclear receptors in living cells and represses transcriptional activity.

Transcriptional activity relies on coregulators that modify the chromatin structure and serve as bridging factors between transcription factors and the basal transcription machinery. Using the DE domain of human peroxisome proliferator-activated receptor gamma (PPARgamma) as bait in a yeast two-hybrid screen of a human adipose tissue library, we isolated the scaffold attachment factor B1 (SAFB1/HET/HAP), which was previously shown to be a corepressor of estrogen receptor alpha. We show here that SAFB1 has a very broad tissue expression profile in human and is also expressed all along mouse embryogenesis. SAFB1 interacts in pull-down assays not only with PPARgamma but also with all nuclear receptors tested so far, albeit with different affinities. The association of SAFB1 and PPARgamma in vivo is further demonstrated by fluorescence resonance energy transfer (FRET) experiments in living cells. We finally show that SAFB1 is a rather general corepressor for nuclear receptors. Its change in expression during the early phases of adipocyte and enterocyte differentiation suggests that SAFB1 potentially influences cell proliferation and differentiation decisions.

Genome sequence of the pea aphid Acyrthosiphon pisum

Aphids are important agricultural pests and also biological models for studies of insect-plant interactions, symbiosis, virus vectoring, and the developmental causes of extreme phenotypic plasticity. Here we present the 464 Mb draft genome assembly of the pea aphid Acyrthosiphon pisum. This first published whole genome sequence of a basal hemimetabolous insect provides an outgroup to the multiple published genomes of holometabolous insects. Pea aphids are host-plant specialists, they can reproduce both sexually and asexually, and they have coevolved with an obligate bacterial symbiont. Here we highlight findings from whole genome analysis that may be related to these unusual biological features. These findings include discovery of extensive gene duplication in more than 2000 gene families as well as loss of evolutionarily conserved genes. Gene family expansions relative to other published genomes include genes involved in chromatin modification, miRNA synthesis, and sugar transport. Gene losses include genes central to the IMD immune pathway, selenoprotein utilization, purine salvage, and the entire urea cycle. The pea aphid genome reveals that only a limited number of genes have been acquired from bacteria; thus the reduced gene count of Buchnera does not reflect gene transfer to the host genome. The inventory of metabolic genes in the pea aphid genome suggests that there is extensive metabolite exchange between the aphid and Buchnera, including sharing of amino acid biosynthesis between the aphid and Buchnera. The pea aphid genome provides a foundation for post-genomic studies of fundamental biological questions and applied agricultural problems.

dKDM5/LID regulates H3K4me3 dynamics at the transcription-start site (TSS) of actively transcribed developmental genes

H3K4me3 is a histone modification that accumulates at the transcription-start site (TSS) of active genes and is known to be important for transcription activation. The way in which H3K4me3 is regulated at TSS and the actual molecular basis of its contribution to transcription remain largely unanswered. To address these questions, we have analyzed the contribution of dKDM5/LID, the main H3K4me3 demethylase in Drosophila, to the regulation of the pattern of H3K4me3. ChIP-seq results show that, at developmental genes, dKDM5/LID localizes at TSS and regulates H3K4me3. dKDM5/LID target genes are highly transcribed and enriched in active RNApol II and H3K36me3, suggesting a positive contribution to transcription. Expression-profiling show that, though weakly, dKDM5/LID target genes are significantly downregulated upon dKDM5/LID depletion. Furthermore, dKDM5/LID depletion results in decreased RNApol II occupancy, particularly by the promoter-proximal Pol lloser5 form. Our results also show that ASH2, an evolutionarily conserved factor that locates at TSS and is required for H3K4me3, binds and positively regulates dKDM5/LID target genes. However, dKDM5/LID and ASH2 do not bind simultaneously and recognize different chromatin states, enriched in H3K4me3 and not, respectively. These results indicate that, at developmental genes, dKDM5/LID and ASH2 coordinately regulate H3K4me3 at TSS and that this dynamic regulation contributes to transcription.

Integrated methodologies to study human transcriptional regulation from functional genomics data

Abstract : The human body is composed of a huge number of cells acting together in a concerted manner. The current understanding is that proteins perform most of the necessary activities in keeping a cell alive. The DNA, on the other hand, stores the information on how to produce the different proteins in the genome. Regulating gene transcription is the first important step that can thus affect the life of a cell, modify its functions and its responses to the environment. Regulation is a complex operation that involves specialized proteins, the transcription factors. Transcription factors (TFs) can bind to DNA and activate the processes leading to the expression of genes into new proteins. Errors in this process may lead to diseases. In particular, some transcription factors have been associated with a lethal pathological state, commonly known as cancer, associated with uncontrolled cellular proliferation, invasiveness of healthy tissues and abnormal responses to stimuli. Understanding cancer-related regulatory programs is a difficult task, often involving several TFs interacting together and influencing each other's activity. This Thesis presents new computational methodologies to study gene regulation. In addition we present applications of our methods to the understanding of cancer-related regulatory programs. The understanding of transcriptional regulation is a major challenge. We address this difficult question combining computational approaches with large collections of heterogeneous experimental data. In detail, we design signal processing tools to recover transcription factors binding sites on the DNA from genome-wide surveys like chromatin immunoprecipitation assays on tiling arrays (ChIP-chip). We then use the localization about the binding of TFs to explain expression levels of regulated genes. In this way we identify a regulatory synergy between two TFs, the oncogene C-MYC and SP1. C-MYC and SP1 bind preferentially at promoters and when SP1 binds next to C-NIYC on the DNA, the nearby gene is strongly expressed. The association between the two TFs at promoters is reflected by the binding sites conservation across mammals, by the permissive underlying chromatin states 'it represents an important control mechanism involved in cellular proliferation, thereby involved in cancer. Secondly, we identify the characteristics of TF estrogen receptor alpha (hERa) target genes and we study the influence of hERa in regulating transcription. hERa, upon hormone estrogen signaling, binds to DNA to regulate transcription of its targets in concert with its co-factors. To overcome the scarce experimental data about the binding sites of other TFs that may interact with hERa, we conduct in silico analysis of the sequences underlying the ChIP sites using the collection of position weight matrices (PWMs) of hERa partners, TFs FOXA1 and SP1. We combine ChIP-chip and ChIP-paired-end-diTags (ChIP-pet) data about hERa binding on DNA with the sequence information to explain gene expression levels in a large collection of cancer tissue samples and also on studies about the response of cells to estrogen. We confirm that hERa binding sites are distributed anywhere on the genome. However, we distinguish between binding sites near promoters and binding sites along the transcripts. The first group shows weak binding of hERa and high occurrence of SP1 motifs, in particular near estrogen responsive genes. The second group shows strong binding of hERa and significant correlation between the number of binding sites along a gene and the strength of gene induction in presence of estrogen. Some binding sites of the second group also show presence of FOXA1, but the role of this TF still needs to be investigated. Different mechanisms have been proposed to explain hERa-mediated induction of gene expression. Our work supports the model of hERa activating gene expression from distal binding sites by interacting with promoter bound TFs, like SP1. hERa has been associated with survival rates of breast cancer patients, though explanatory models are still incomplete: this result is important to better understand how hERa can control gene expression. Thirdly, we address the difficult question of regulatory network inference. We tackle this problem analyzing time-series of biological measurements such as quantification of mRNA levels or protein concentrations. Our approach uses the well-established penalized linear regression models where we impose sparseness on the connectivity of the regulatory network. We extend this method enforcing the coherence of the regulatory dependencies: a TF must coherently behave as an activator, or a repressor on all its targets. This requirement is implemented as constraints on the signs of the regressed coefficients in the penalized linear regression model. Our approach is better at reconstructing meaningful biological networks than previous methods based on penalized regression. The method is tested on the DREAM2 challenge of reconstructing a five-genes/TFs regulatory network obtaining the best performance in the "undirected signed excitatory" category. Thus, these bioinformatics methods, which are reliable, interpretable and fast enough to cover large biological dataset, have enabled us to better understand gene regulation in humans.

A multiplicity of factors contributes to selective RNA polymerase III occupancy of a subset of RNA polymerase III genes in mouse liver.

The genomic loci occupied by RNA polymerase (RNAP) III have been characterized in human culture cells by genome-wide chromatin immunoprecipitations, followed by deep sequencing (ChIP-seq). These studies have shown that only ∼40% of the annotated 622 human tRNA genes and pseudogenes are occupied by RNAP-III, and that these genes are often in open chromatin regions rich in active RNAP-II transcription units. We have used ChIP-seq to characterize RNAP-III-occupied loci in a differentiated tissue, the mouse liver. Our studies define the mouse liver RNAP-III-occupied loci including a conserved mammalian interspersed repeat (MIR) as a potential regulator of an RNAP-III subunit-encoding gene. They reveal that synteny relationships can be established between a number of human and mouse RNAP-III genes, and that the expression levels of these genes are significantly linked. They establish that variations within the A and B promoter boxes, as well as the strength of the terminator sequence, can strongly affect RNAP-III occupancy of tRNA genes. They reveal correlations with various genomic features that explain the observed variation of 81% of tRNA scores. In mouse liver, loci represented in the NCBI37/mm9 genome assembly that are clearly occupied by RNAP-III comprise 50 Rn5s (5S RNA) genes, 14 known non-tRNA RNAP-III genes, nine Rn4.5s (4.5S RNA) genes, and 29 SINEs. Moreover, out of the 433 annotated tRNA genes, half are occupied by RNAP-III. Transfer RNA gene expression levels reflect both an underlying genomic organization conserved in dividing human culture cells and resting mouse liver cells, and the particular promoter and terminator strengths of individual genes.

The role of PKCzeta in amikacin-induced apoptosis in the cochlea: prevention by aspirin.

Aminoglycoside antibiotics are ototoxic, inducing irreversible sensorineural hearing loss mediated by oxidative and excitotoxic stresses. The NF-kappaB pathway is involved in the response to aminoglycoside damage in the cochlea. However, the molecular mechanisms of this ototoxicity remain unclear. We investigated the expression of PKCzeta, a key regulator of NF-kappaB activation, in response to aminoglycoside treatment. Amikacin induced PKCzeta cleavage and nuclear translocation. These events were concomitant with chromatin condensation and paralleled the decrease in NF-kappaB (p65) levels in the nucleus. Amikacin also induced the nuclear translocation of apoptotic inducing factor (AIF). Prior treatment with aspirin prevented PKCzeta cleavage and nuclear translocation. Thus, aspirin counteracts the early effects of amikacin, thereby protecting hair cells and spiral ganglion neurons. These results demonstrate that PKCzeta acts as sentinel connecting specific survival pathways to mediate cellular responses to amikacin ototoxicity.

Myc's other life: stem cells and beyond.

Over the last three decades genetic and biochemical studies have revealed the pleiotropic effects of the Myc oncoprotein. While cell line studies have defined the intracellular processes regulated by Myc such as proliferation, differentiation, and metabolic growth, in vivo studies have confirmed these functions, and revealed roles in acquisition and maintenance of stem cell properties. These roles may be partially mediated by Myc's capacity to modify the chromatin landscape on a global scale. Myc also regulates numerous protein-coding transcripts, and many noncoding RNAs (rRNAs, tRNAs, and miRNAs). As Myc activity directly correlates with protein expression, further complexity is provided by post-translational modifications that regulate Myc in normal stem cells or deregulate it in malignant stem cells.

Transcriptional potentiation of the vitellogenin B1 promoter by a combination of both nucleosome assembly and transcription factors: an in vitro dissection.

The Xenopus laevis vitellogenin B1 promoter was assembled into nucleosomes in an oocyte extract. Subsequent RNA polymerase II-dependent transcription from these DNA templates fully reconstituted in chromatin in a HeLa nuclear extract was increased 50-fold compared with naked DNA. Remarkably, under specific conditions, production of a high level of transcripts occurred at very low DNA (1 ng/microliter) and HeLa nuclear protein (1.6 micrograms/microliters) concentrations. When partially reconstituted templates were used, transcription efficiency was intermediate between that of fully reconstituted and naked DNA. These results implicate chromatin in the process of the transcriptional activation observed. Depletion from the oocyte assembly extract of an NF-I-like factor which binds in the promoter region upstream of the TATA box (-114 to -101) or deletion from the promoter of the region interacting with this factor reduced the transcriptional efficiency of the assembled templates by a factor of 5, but transcription of these templates was still 10 times higher than that of naked DNA. Together, these results indicate that the NF-I-like factor participates in the very efficient transcriptional potentiation of the vitellogenin B1 promoter which occurs during nucleosome assembly.