Une nouvelle approche computationnelle pour la découverte des sites de fixation de facteurs de transcription à l’ADN, adaptée aux données de ChIP-chip et de ChIP-séquençage

Les facteurs de transcription sont des protéines spécialisées qui jouent un rôle important dans différents processus biologiques tel que la différenciation, le cycle cellulaire et la tumorigenèse. Ils régulent la transcription des gènes en se fixant sur des séquences d’ADN spécifiques (éléments cis-régulateurs). L’identification de ces éléments est une étape cruciale dans la compréhension des réseaux de régulation des gènes. Avec l’avènement des technologies de séquençage à haut débit, l’identification de tout les éléments fonctionnels dans les génomes, incluant gènes et éléments cis-régulateurs a connu une avancée considérable. Alors qu’on est arrivé à estimer le nombre de gènes chez différentes espèces, l’information sur les éléments qui contrôlent et orchestrent la régulation de ces gènes est encore mal définie. Grace aux techniques de ChIP-chip et de ChIP-séquençage il est possible d’identifier toutes les régions du génome qui sont liées par un facteur de transcription d’intérêt. Plusieurs approches computationnelles ont été développées pour prédire les sites fixés par les facteurs de transcription. Ces approches sont classées en deux catégories principales: les algorithmes énumératifs et probabilistes. Toutefois, plusieurs études ont montré que ces approches génèrent des taux élevés de faux négatifs et de faux positifs ce qui rend difficile l’interprétation des résultats et par conséquent leur validation expérimentale. Dans cette thèse, nous avons ciblé deux objectifs. Le premier objectif a été de développer une nouvelle approche pour la découverte des sites de fixation des facteurs de transcription à l’ADN (SAMD-ChIP) adaptée aux données de ChIP-chip et de ChIP-séquençage. Notre approche implémente un algorithme hybride qui combine les deux stratégies énumérative et probabiliste, afin d’exploiter les performances de chacune d’entre elles. Notre approche a montré ses performances, comparée aux outils de découvertes de motifs existants sur des jeux de données simulées et des jeux de données de ChIP-chip et de ChIP-séquençage. SAMD-ChIP présente aussi l’avantage d’exploiter les propriétés de distributions des sites liés par les facteurs de transcription autour du centre des régions liées afin de limiter la prédiction aux motifs qui sont enrichis dans une fenêtre de longueur fixe autour du centre de ces régions. Les facteurs de transcription agissent rarement seuls. Ils forment souvent des complexes pour interagir avec l’ADN pour réguler leurs gènes cibles. Ces interactions impliquent des facteurs de transcription dont les sites de fixation à l’ADN sont localisés proches les uns des autres ou bien médier par des boucles de chromatine. Notre deuxième objectif a été d’exploiter la proximité spatiale des sites liés par les facteurs de transcription dans les régions de ChIP-chip et de ChIP-séquençage pour développer une approche pour la prédiction des motifs composites (motifs composés par deux sites et séparés par un espacement de taille fixe). Nous avons testé ce module pour prédire la co-localisation entre les deux demi-sites ERE qui forment le site ERE, lié par le récepteur des œstrogènes ERα. Ce module a été incorporé à notre outil de découverte de motifs SAMD-ChIP.

Evolutionary turnover of mammalian transcription start sites

Alignments of homologous genomic sequences are widely used to identify functional genetic elements and study their evolution. Most studies tacitly equate homology of functional elements with sequence homology. This assumption is violated by the phenomenon of turnover, in which functionally equivalent elements reside at locations that are nonorthologous at the sequence level. Turnover has been demonstrated previously for transcription-factor-binding sites. Here, we show that transcription start sites of equivalent genes do not always reside at equivalent locations in the human and mouse genomes. We also identify two types of partial turnover, illustrating evolutionary pathways that could lead to complete turnover. These findings suggest that the signals encoding transcription start sites are highly flexible and evolvable, and have cautionary implications for the use of sequence-level conservation to detect gene regulatory elements.

Selective inhibition of the human tie-1 promoter with triplex-forming oligonucleotides targeted to ets binding sites

The Tie receptors (Tie-1 and Tie-2/Tek) are essential for angiogenesis and vascular remodeling/integrity. Tie receptors are up-regulated in tumor-associated endothelium, and their inhibition disrupts angiogenesis and can prevent tumor growth as a consequence. To investigate the potential of anti-gene approaches to inhibit tie gene expression for anti-angiogenic therapy, we have examined triple-helical (triplex) DNA formation at 2 tandem Ets transcription factor binding motifs (designated E-1 and E-2) in the human tie-1 promoter. Various tie-1 promoter deletion/mutation luciferase reporter constructs were generated and transfected into endothelial cells to examine the relative activities of E-1 and E-2. The binding of antiparallel and parallel (control) purine motif oligonucleotides (21-22 bp) targeted to E-1 and E-2 was assessed by plasmid DNA fragment binding and electrophoretic mobility shift assays. Triplex-forming oligonucleotides were incubated with tie-1 reporter constructs and transfected into endothelial cells to determine their activity. The Ets binding motifs in the E-1 sequence were essential for human tie-1 promoter activity in endothelial cells, whereas the deletion of E-2 had no effect. Antiparallel purine motif oligonucleotides targeted at E-1 or E-2 selectively formed strong triplex DNA (K(d) approximately 10(-7) M) at 37 degrees C. Transfection of tie-1 reporter constructs with triplex DNA at E-1, but not E-2, specifically inhibited tie-1 promoter activity by up to 75% compared with control oligonucleotides in endothelial cells. As similar multiple Ets binding sites are important for the regulation of several endothelial-restricted genes, this approach may have broad therapeutic potential for cancer and other pathologies involving endothelial proliferation/dysfunction.

The human sulfotransferase SULT1A1 gene is regulated in a synergistic manner by Sp1 and GA binding protein

Human sulfotransferase SULT1A1 is an important phase II xenobiotic metabolizing enzyme that is highly expressed in the liver and mediates the sulfonation of drugs, carcinogens, and steroids. Until this study, the transcriptional regulation of the SULT1A subfamily had been largely unexplored. Preliminary experiments in primary human hepatocytes showed that SULT1A mRNA levels were not changed in response to nuclear receptor activators, such as dexamethasone and 3-methylcolanthrene, unlike other metabolizing enzymes. Using HepG2 cells, the high activity of the TATA-less SULT1A1 promoter was shown to be dependent on the presence of Sp1 and Ets transcription factor binding sites (EBS), located within - 112 nucleotides from the transcriptional start site. The homologous promoter of the closely related SULT1A3 catecholamine sulfotransferase, which is expressed at negligible levels in the adult liver, displayed 70% less activity than SULT1A1. This was shown to be caused by a two-base pair difference in the EBS. The Ets transcription factor GA binding protein (GABP) was shown to bind the SULT1A1 EBS and could transactivate the SULT1A1 promoter in Drosophila melanogaster S2 cells. Cotransfection of Sp1 could synergistically enhance GABP-mediated activation by 10-fold. Although Sp1 and GABP alone could induce SULT1A3 promoter activity, the lack of the EBS on this promoter prevented a synergistic interaction between the two factors. This study reports the first insight into the transcriptional regulation of the SULT1A1 gene and identifies a crucial difference in regulation of the closely related SULT1A3 gene, which accounts for the two enzymes' differential expression patterns observed in the adult liver.

Genomic organization of the Neurospora crassa gsn gene: possible involvement of the STRE and HSE elements in the modulation of transcription during heat shock

Glycogen synthase, an enzyme involved in glycogen biosynthesis, is regulated by phosphorylation and by the allosteric ligand glucose-6-phosphate (G6P). In addition, enzyme levels can be regulated by changes in gene expression. We recently cloned a cDNA for glycogen synthase (gsn) from Neurospora crassa, and showed that gsn transcription decreased when cells were exposed to heat shock (shifted from 30degreesC to 45degreesC). In order to understand the mechanisms that control gsn expression, we isolated the gene, including its 5' and 3' flanking regions, from the genome of N. crassa. An ORF of approximately 2.4 kb was identified, which is interrupted by four small introns (II-V). Intron I (482 bp) is located in the 5'UTR region. Three putative Transcription Initiation Sites (TISs) were mapped, one of which lies downstream of a canonical TATA-box sequence (5'-TGTATAAA-3'). Analysis of the 5'-flanking region revealed the presence of putative transcription factor-binding sites, including Heat Shock Elements (HSEs) and STress Responsive Elements (STREs). The possible involvement of these motifs in the negative regulation of gsn transcription was investigated using Electrophoretic Mobility Shift Assays (EMSA) with nuclear extracts of N. crassa mycelium obtained before and after heat shock, and DNA fragments encompassing HSE and STRE elements from the 5'-flanking region. While elements within the promoter region are involved in transcription under heat shock, elements in the 5'UTR intron may participate in transcription during vegetative growth. The results thus suggest that N. crassa possesses trans-acting elements that interact with the 5'-flanking region to regulate gsn transcription during heat shock and vegetative growth.

Positive and negative regulation of endogenous genes by designed transcription factors

Gene regulation by imposed localization was studied by using designed zinc finger proteins that bind 18-bp DNA sequences in the 5′ untranslated regions of the protooncogenes erbB-2 and erbB-3. Transcription factors were generated by fusion of the DNA-binding proteins to repression or activation domains. When introduced into cells these transcription factors acted as dominant repressors or activators of, respectively, endogenous erbB-2 or erbB-3 gene expression. Significantly, imposed regulation of the two genes was highly specific, despite the fact that the transcription factor binding sites targeted in erbB-2 and erbB-3 share 15 of 18 nucleotides. Regulation of erbB-2 gene expression was observed in cells derived from several species that conserve the DNA target sequence. Repression of erbB-2 in SKBR3 breast cancer cells inhibited cell-cycle progression by inducing a G1 accumulation, suggesting the potential of designed transcription factors for cancer gene therapy. These results demonstrate the willful up- and down-regulation of endogenous genes, and provide an additional means to alter biological systems.

Exploring the DNA-binding specificities of zinc fingers with DNA microarrays

A key step in the regulation of networks that control gene expression is the sequence-specific binding of transcription factors to their DNA recognition sites. A more complete understanding of these DNA–protein interactions will permit a more comprehensive and quantitative mapping of the regulatory pathways within cells, as well as a deeper understanding of the potential functions of individual genes regulated by newly identified DNA-binding sites. Here we describe a DNA microarray-based method to characterize sequence-specific DNA recognition by zinc-finger proteins. A phage display library, prepared by randomizing critical amino acid residues in the second of three fingers of the mouse Zif268 domain, provided a rich source of zinc-finger proteins with variant DNA-binding specificities. Microarrays containing all possible 3-bp binding sites for the variable zinc fingers permitted the quantitation of the binding site preferences of the entire library, pools of zinc fingers corresponding to different rounds of selection from this library, as well as individual Zif268 variants that were isolated from the library by using specific DNA sequences. The results demonstrate the feasibility of using DNA microarrays for genome-wide identification of putative transcription factor-binding sites.

Sequence-selective carbohydrate-DNA interaction: dimeric and monomeric forms of the calicheamicin oligosaccharide interfere with transcription factor function.

The synthetic oligosaccharide moiety of the antibiotic calicheamicin and the head-to-head dimer of this oligosaccharide are known to bind to the minor groove of DNA in a sequence-selective manner preferring distinct target sequences. We tested these carbohydrates for their ability to interfere with transcription factor function. The oligosaccharides inhibit binding of transcription factors to DNA in a sequence-selective manner, probably by inducing a conformational change in DNA structure. They also interfere with transcription by polymerase II in vitro. The effective concentrations of the oligosaccharides for inhibition of transcription factor binding and for transcriptional inhibition are in the micromolar range. The dimer is a significantly more active inhibitor than is the monomer.

A 20-nucleotide element in the intestinal fatty acid binding protein gene modulates its cell lineage-specific, differentiation-dependent, and cephalocaudal patterns of expression in transgenic mice.

A sequence of epithelial cell proliferation, allocation to four principal lineages, migration-associated differentiation, and cell loss occurs along the crypt-villus axis of the mouse intestine. The sequence is completed in a few days and is recapitulated throughout the life-span of the animal. We have used an intestine-specific fatty acid binding protein gene, Fabpi, as a model for studying regulation of gene expression in this unique developmental system. Promoter mapping studies in transgenic mice identified a 20-bp cis-acting element (5'-AGGTGGAAGCCATCACACTT-3') that binds small intestinal nuclear proteins and participates in the control of Fabpi's cephalocaudal, differentiation-dependent, and cell lineage-specific patterns of expression. Immunocytochemical studies using confocal and electron microscopy indicate that it does so by acting as a suppressor of gene expression in the distal small intestine/colon, as a suppressor of gene activation in proliferating and nonproliferating cells located in the crypts of Lieberkühn, and as a suppressor of expression in the growth factor and defensin-producing Paneth cell lineage. The 20-bp domain has no obvious sequence similarities to known transcription factor binding sites. The three functions modulated by this compact element represent the types of functions required to establish and maintain the intestine's remarkably complex spatial patterns of gene expression. The transgenes described in this report also appear to be useful in characterizing the crypt's stem cell hierarchy.

Mating system and avpr1a promoter variation in primates.

It has been suggested that primate mating and social behaviours may be influenced by variation in promoter region repetitive DNA of the vasopressin receptor 1a gene (avpr1a). We show that male mating behaviour does not covary in a simple way with promoter repetitive DNA in 12 Old World primates. We found that one microsatellite (-553 bp upstream) was present in all species, irrespective of their behaviour. By contrast, two microsatellites (-3956 and -3625 bp upstream) were present only in some species, yet this variation did not correlate with behaviour. These findings agree with a recent comparative analysis of voles and show that the variation in repetitive DNA in the avpr1a promoter region does not generally explain variation in male mating behaviour. Phylogenetic analysis revealed a GAGTA motif that has been independently deleted three times and involved in another larger deletion. Importantly, the presence/absence of this GAGTA motif leads to changes in predicted transcription factor-binding sites. Given the repeated loss of this motif, we speculate that it might be of functional relevance. We suggest that such non-repetitive variation, either in indels or in sequence variation, are likely to be important in explaining interspecific variation in avpr1a expression.

Social insect genomes exhibit dramatic evolution in gene composition and regulation while preserving regulatory features linked to sociality.

Genomes of eusocial insects code for dramatic examples of phenotypic plasticity and social organization. We compared the genomes of seven ants, the honeybee, and various solitary insects to examine whether eusocial lineages share distinct features of genomic organization. Each ant lineage contains ∼4000 novel genes, but only 64 of these genes are conserved among all seven ants. Many gene families have been expanded in ants, notably those involved in chemical communication (e.g., desaturases and odorant receptors). Alignment of the ant genomes revealed reduced purifying selection compared with Drosophila without significantly reduced synteny. Correspondingly, ant genomes exhibit dramatic divergence of noncoding regulatory elements; however, extant conserved regions are enriched for novel noncoding RNAs and transcription factor-binding sites. Comparison of orthologous gene promoters between eusocial and solitary species revealed significant regulatory evolution in both cis (e.g., Creb) and trans (e.g., fork head) for nearly 2000 genes, many of which exhibit phenotypic plasticity. Our results emphasize that genomic changes can occur remarkably fast in ants, because two recently diverged leaf-cutter ant species exhibit faster accumulation of species-specific genes and greater divergence in regulatory elements compared with other ants or Drosophila. Thus, while the "socio-genomes" of ants and the honeybee are broadly characterized by a pervasive pattern of divergence in gene composition and regulation, they preserve lineage-specific regulatory features linked to eusociality. We propose that changes in gene regulation played a key role in the origins of insect eusociality, whereas changes in gene composition were more relevant for lineage-specific eusocial adaptations.

Phylogeny and evolution of the aspartyl protease family from clinically relevant Candida species

Aspartyl proteases are a class of enzymes that include the yeast aspartyl proteases and secreted aspartyl protease (Sap) superfamilies. Several Sap superfamily members have been demonstrated or suggested as virulence factors in opportunistic pathogens of the genus Candida. Candida albicans, Candida tropicalis, Candida dubliniensis and Candida parapsilosis harbour 10, four, eight and three SAP genes, respectively. In this work, genome mining and phylogenetic analyses revealed the presence of new members of the Sap superfamily in C. tropicalis (8), Candida guilliermondii (8), C. parapsilosis(11) and Candida lusitaniae (3). A total of 12 Sap families, containing proteins with at least 50% similarity, were discovered in opportunistic, pathogenic Candida spp. In several Sap families, at least two subfamilies or orthologous groups were identified, each defined by > 90% sequence similitude, functional similarity and synteny among its members. No new members of previously described Sap families were found in a Candida spp. clinical strain collection; however, the universality of SAPT gene distribution among C. tropicalis strains was demonstrated. In addition, several features of opportunistic pathogenic Candida species, such as gene duplications and inversions, similitude, synteny, putative transcription factor binding sites and genome traits of SAP gene superfamily were described in a molecular evolutionary context.

Multiple non-collinear TF-map alignments of promoter regions

Background: The analysis of the promoter sequence of genes with similar expression patterns isa basic tool to annotate common regulatory elements. Multiple sequence alignments are on thebasis of most comparative approaches. The characterization of regulatory regions from coexpressedgenes at the sequence level, however, does not yield satisfactory results in manyoccasions as promoter regions of genes sharing similar expression programs often do not shownucleotide sequence conservation.Results: In a recent approach to circumvent this limitation, we proposed to align the maps ofpredicted transcription factors (referred as TF-maps) instead of the nucleotide sequence of tworelated promoters, taking into account the label of the corresponding factor and the position in theprimary sequence. We have now extended the basic algorithm to permit multiple promotercomparisons using the progressive alignment paradigm. In addition, non-collinear conservationblocks might now be identified in the resulting alignments. We have optimized the parameters ofthe algorithm in a small, but well-characterized collection of human-mouse-chicken-zebrafishorthologous gene promoters.Conclusion: Results in this dataset indicate that TF-map alignments are able to detect high-levelregulatory conservation at the promoter and the 3'UTR gene regions, which cannot be detectedby the typical sequence alignments. Three particular examples are introduced here to illustrate thepower of the multiple TF-map alignments to characterize conserved regulatory elements inabsence of sequence similarity. We consider this kind of approach can be extremely useful in thefuture to annotate potential transcription factor binding sites on sets of co-regulated genes fromhigh-throughput expression experiments.

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/.

Integrative analysis of the Regulatory Region of the FGFR3 Oncogene

The study of transcriptional regulation often needs the integration of diverse yet independent data. In the present work, sequence conservation, predic-tion of transcription factor binding sites (TFBS) and gene expression analysis have been applied to the detection of putative transcription factor (TF) modules in the regulatory region of the FGFR3 oncogene. Several TFs with conserved binding sites in the FGFR3 regulatory region have shown high positive or negative corre-lation with FGFR3 expression both in urothelial carcinoma and in benign nevi. By means of conserved TF cluster analysis, two different TF modules have been iden-tified in the promoter and first intron of FGFR3 gene. These modules contain acti-vating AP2, E2F, E47 and SP1 binding sites plus motifs for EGR with possible repressor function.