Modular cis-regulatory organization of developmentally expressed genes: two genes transcribed territorially in the sea urchin embryo, and additional examples.

The cis-regulatory systems that control developmental expression of two sea urchin genes have been subjected to detailed functional analysis. Both systems are modular in organization: specific, separable fragments of the cis-regulatory DNA each containing multiple transcription factor target sites execute particular regulatory subfunctions when associated with reporter genes and introduced into the embryo. The studies summarized here were carried out on the CyIIIa gene, expressed in the embryonic aboral ectoderm and on the Endo16 gene, expressed in the embryonic vegetal plate, archenteron, and then midgut. The regulatory systems of both genes include modules that control particular aspects of temporal and spatial expression, and in both the territorial boundaries of expression depend on a combination of negative and positive functions. In both genes different regulatory modules control early and late embryonic expression. Modular cis-regulatory organization is widespread in developmentally regulated genes, and we present a tabular summary that includes many examples from mouse and Drosophila. We regard cis-regulatory modules as units of developmental transcription control, and also of evolution, in the assembly of transcription control systems.

Sequence analysis of the cis-regulatory regions of the bithorax complex of Drosophila.

The bithorax complex (BX-C) of Drosophila, one of two complexes that act as master regulators of the body plan of the fly, has now been entirely sequenced and comprises approximately 315,000 bp, only 1.4% of which codes for protein. Analysis of this sequence reveals significantly overrepresented DNA motifs of unknown, as well as known, functions in the non-protein-coding portion of the sequence. The following types of motifs in that portion are analyzed: (i) concatamers of mono-, di-, and trinucleotides; (ii) tightly clustered hexanucleotides (spaced < or = 5 bases apart); (iii) direct and reverse repeats longer than 20 bp; and (iv) a number of motifs known from biochemical studies to play a role in the regulation of the BX-C. The hexanucleotide AGATAC is remarkably overrepresented and is surmised to play a role in chromosome pairing. The positions of sites of highly overrepresented motifs are plotted for those that occur at more than five sites in the sequence, when < 0.5 case is expected. Expected values are based on a third-order Markov chain, which is the optimal order for representing the BXCALL sequence.

Elucidating the cis-regulatory landscape of retinal development and regeneration

Thesis (Ph.D.)--University of Washington, 2016-04

Functional characterisation of putative Cis-Regulatory Risk Loci for breast cancer

Dissertação de Mestrado, Qualidade em Análises - Erasmus Mundus, Faculdade de Ciências e Tecnologia, Universidade do Algarve, 2015

Analysis of opo cis-regulatory landscape uncovers Vsx2 requirement in early eye morphogenesis

The self-organized morphogenesis of the vertebrate optic cup entails coupling the activation of the retinal gene regulatory network to the constriction-driven infolding of the retinal epithelium. Yet the genetic mechanisms underlying this coordination remain largely unexplored. Through phylogenetic footprinting and transgenesis in zebrafish, here we examine the cis-regulatory landscape of opo, an endocytosis regulator essential for eye morphogenesis. Among the different conserved enhancers identified, we isolate a single retina-specific element (H6_10137) and show that its activity depends on binding sites for the retinal determinant Vsx2. Gain- and loss-of-function experiments and ChIP analyses reveal that Vsx2 regulates opo expression through direct binding to this retinal enhancer. Furthermore, we show that vsx2 knockdown impairs the primary optic cup folding. These data support a model by which vsx2, operating through the effector gene opo, acts as a central transcriptional node that coordinates neural retina patterning and optic cup invagination in zebrafish.

Site-directed mutagenesis of the ATM promoter: Consequences for response to proliferation and ionizing radiation

Although ATM, the protein defective in ataxia-telangiectasia (A-T), is activated primarily by radiation, there is also evidence that expression of the protein can be regulated by both radiation and growth factors. Computer analysis of the ATM promoter proximal 700-bp sequence reveals a number of potentially important cis-regulatory sequences. Using nucleotide substitutions to delete putative functional elements in the promoter of ATM, we examined the importance of some of these sites for both the basal and the radiation-induced activity of the promoter. In lymphoblastoid cells, most of the mutations in transcription factor consensus sequences [Sp1(1), Sp1(2), Cre, Ets, Xre, gammaIre(2), a modified AP1 site (Fse), and GCF] reduced basal activity to various extents, whereas others [gammaIre(1), NF1, Myb] left basal activity unaffected. In human skin fibroblasts, results were generally the same, but the basal activity varied up to 8-fold in these and other cell lines. Radiation activated the promoter approximately 2.5-fold in serum-starved lymphoblastoid cells, reaching a maximum by 3 hr, and all mutated elements equally blocked this activation. Reduction in Sp1 and AP1 DNA binding activity by serum starvation was rapidly reversed by exposure of cells to radiation. This reduction was not evident in A-T cells, and the response to radiation was less marked. Data provided for interaction between ATM and Sp1 by protein binding and co-immunoprecipitation could explain the altered regulation of Sp1 in A-T cells. The data described here provide additional evidence that basal and radiation-induced regulation of the ATM promoter is under multifactorial control. (C) 2003 Wiley-Liss, Inc.

Sleep loss reduces the DNA-binding of BMAL1, CLOCK, and NPAS2 to specific clock genes in the mouse cerebral cortex.

We have previously demonstrated that clock genes contribute to the homeostatic aspect of sleep regulation. Indeed, mutations in some clock genes modify the markers of sleep homeostasis and an increase in homeostatic sleep drive alters clock gene expression in the forebrain. Here, we investigate a possible mechanism by which sleep deprivation (SD) could alter clock gene expression by quantifying DNA-binding of the core-clock transcription factors CLOCK, NPAS2, and BMAL1 to the cis-regulatory sequences of target clock genes in mice. Using chromatin immunoprecipitation (ChIP), we first showed that, as reported for the liver, DNA-binding of CLOCK and BMAL1 to target clock genes changes in function of time-of-day in the cerebral cortex. Tissue extracts were collected at ZT0 (light onset), -6, -12, and -18, and DNA enrichment of E-box or E'-box containing sequences was measured by qPCR. CLOCK and BMAL1 binding to Cry1, Dbp, Per1, and Per2 depended on time-of-day, with maximum values reached at around ZT6. We then observed that SD, performed between ZT0 and -6, significantly decreased DNA-binding of CLOCK and BMAL1 to Dbp, consistent with the observed decrease in Dbp mRNA levels after SD. The DNA-binding of NPAS2 and BMAL1 to Per2 was also decreased by SD, although SD is known to increase Per2 expression in the cortex. DNA-binding to Per1 and Cry1 was not affected by SD. Our results show that the sleep-wake history can affect the clock molecular machinery directly at the level of chromatin binding thereby altering the cortical expression of Dbp and Per2 and likely other targets. Although the precise dynamics of the relationship between DNA-binding and mRNA expression, especially for Per2, remains elusive, the results also suggest that part of the reported circadian changes in DNA-binding of core clock components in tissues peripheral to the suprachiasmatic nuclei could, in fact, be sleep-wake driven.

Unique Alterations of an Ultraconserved Non-Coding Element in the 3 ' UTR of ZIC2 in Holoprosencephaly

Coding region alterations of ZIC2 are the second most common type of mutation in holoprosencephaly (HPE). Here we use several complementary bioinformatic approaches to identify ultraconserved cis-regulatory sequences potentially driving the expression of human ZIC2. We demonstrate that an 804 bp element in the 3' untranslated region (3'UTR) is highly conserved across the evolutionary history of vertebrates from fish to humans. Furthermore, we show that while genetic variation of this element is unexpectedly common among holoprosencephaly subjects (6/528 or >1%), it is not present in control individuals. Two of six proband-unique variants are de novo, supporting their pathogenic involvement in HPE outcomes. These findings support a general recommendation that the identification and analysis of key ultraconserved elements should be incorporated into the genetic risk assessment of holoprosencephaly cases.

The effects of DNA cytosine methylation on H -ras promoter structure and function

The effect of DNA cytosine methylation on H-ras promoter activity was assessed using a transient expression system employing the plasmid H-rasCAT (NaeI H-ras promoter linked to the chloramphenicol acetyltransferase (CAT) gene). This 551 bp promoter is 80% GC rich, enriched with 168 CpG dinucleotides, and contains six functional GC box elements which represent major DNA methylation target sites. Prokaryotic methyltransferases HhaI (CGm$\sp5$CG) and HpaII (Cm$\sp5$CGG) alone or in combination with a human placental methyltransferase (HP MTase) were used to introduce methyl groups at different CpG sites within the promoter. To test for functional promoter activity, the methylated plasmids were introduced into CV-1 cells and CAT activity assessed 48 h post-transfection. Methylation at specific HhaI and HpaII sites reduced CAT expression by 70%, whereas more extensive methylation at generalized CpG sites with HP MTase inactivated the promoter $>$95%. The inhibition of H-ras promoter activity was not attributable to methylation-induced differences in DNA uptake or stability in the cell, topological form of the plasmid, or methylation effects in nonpromoter regions. We also observed that DNA cytosine methylation of a 360 bp promoter fragment by HP MTase induced a local change in DNA conformation. Using three independent methodologies (nitrocellulose filter binding assays, gel mobility shifts, and Southwestern blots), we determined that this change in promoter conformation affected the interaction of nuclear proteins with cis-regulatory sequences residing in the promoter region. The results provide evidence to suggest that DNA methylation may regulate gene expression by inducing changes in local promoter conformation which in turn alters the interactions between DNA and protein factors required for transcription. The results provide supportive evidence for the hypothesis of Cedar and Riggs, who postulated that DNA methylation may regulate gene expression by altering the binding affinities of proteins for DNA. ^

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.

Identification of Conserved Splicing Motifs in Mutually Exclusive Exons of 15 Insect Species

Background: During alternative splicing, the inclusion of an exon in the final mRNA molecule is determined by nuclear proteins that bind cis-regulatory sequences in a target pre-mRNA molecule. A recent study suggested that the regulatory codes of individual RNA-binding proteins may be nearly immutable between very diverse species such as mammals and insects. The model system Drosophila melanogaster therefore presents an excellent opportunity for the study of alternative splicing due to the availability of quality EST annotations in FlyBase. Methods: In this paper, we describe an in silico analysis pipeline to extract putative exonic splicing regulatory sequences from a multiple alignment of 15 species of insects. Our method, ESTs-to-ESRs (E2E), uses graph analysis of EST splicing graphs to identify mutually exclusive (ME) exons and combines phylogenetic measures, a sliding window approach along the multiple alignment and the Welch’s t statistic to extract conserved ESR motifs. Results: The most frequent 100% conserved word of length 5 bp in different insect exons was “ATGGA”. We identified 799 statistically significant “spike” hexamers, 218 motifs with either a left or right FDR corrected spike magnitude p-value < 0.05 and 83 with both left and right uncorrected p < 0.01. 11 genes were identified with highly significant motifs in one ME exon but not in the other, suggesting regulation of ME exon splicing through these highly conserved hexamers. The majority of these genes have been shown to have regulated spatiotemporal expression. 10 elements were found to match three mammalian splicing regulator databases. A putative ESR motif, GATGCAG, was identified in the ME-13b but not in the ME-13a of Drosophila N-Cadherin, a gene that has been shown to have a distinct spatiotemporal expression pattern of spliced isoforms in a recent study. Conclusions: Analysis of phylogenetic relationships and variability of sequence conservation as implemented in the E2E spikes method may lead to improved identification of ESRs. We found that approximately half of the putative ESRs in common between insects and mammals have a high statistical support (p < 0.01). Several Drosophila genes with spatiotemporal expression patterns were identified to contain putative ESRs located in one exon of the ME exon pairs but not in the other.

Assaying the regulatory potential of mammalian conserved non-coding sequences in human cells.

BACKGROUND: Conserved non-coding sequences in the human genome are approximately tenfold more abundant than known genes, and have been hypothesized to mark the locations of cis-regulatory elements. However, the global contribution of conserved non-coding sequences to the transcriptional regulation of human genes is currently unknown. Deeply conserved elements shared between humans and teleost fish predominantly flank genes active during morphogenesis and are enriched for positive transcriptional regulatory elements. However, such deeply conserved elements account for <1% of the conserved non-coding sequences in the human genome, which are predominantly mammalian. RESULTS: We explored the regulatory potential of a large sample of these 'common' conserved non-coding sequences using a variety of classic assays, including chromatin remodeling, and enhancer/repressor and promoter activity. When tested across diverse human model cell types, we find that the fraction of experimentally active conserved non-coding sequences within any given cell type is low (approximately 5%), and that this proportion increases only modestly when considered collectively across cell types. CONCLUSIONS: The results suggest that classic assays of cis-regulatory potential are unlikely to expose the functional potential of the substantial majority of mammalian conserved non-coding sequences in the human genome.

Pancreatic-specific expression of the glucose transporter type 2 gene: identification of cis-elements and islet-specific trans-acting factors.

A defect in glucose sensing of the pancreatic beta-cells has been observed in several animal models of type II diabetes and has been correlated with a reduced gene expression of the glucose transporter type 2 (Glut2). In a transgenic mouse model, expression of Glut2 antisense RNA in pancreatic beta-cells has recently been shown to be associated with an impaired glucose-induced insulin secretion and the development of diabetes. To identify factors that may be involved in the specific decrease of Glut2 in the beta-cells of the diabetic animal, an attempt was made to localize the cis-elements and trans-acting factors involved in the control of Glut2 expression in the endocrine pancreas. It was demonstrated by transient transfection studies that only 338 base pairs (bp) of the murine Glut2 proximal promoter are needed for reporter gene expression in pancreatic islet-derived cell lines, whereas no activity was detected in nonpancreatic cells. Three cis-elements, GTI, GTII, and GTIII, have been identified by DNAse I footprinting and gel retardation experiments within these 338 bp. GTI and GTIII bind distinct but ubiquitously expressed trans-acting factors. On the other hand, nuclear proteins specifically expressed in pancreatic cell lines interact with GTII, and their relative abundance correlates with endogenous Glut2 expression. These GTII-binding factors correspond to nuclear proteins of 180 and 90 kilodaltons as defined by Southwestern analysis. The 180-kilodalton factor is present in pancreatic beta-cell lines but not in an alpha-cell line. Mutation of the GTI or GTIII cis-elements decreases transcriptional activity directed by the 338-bp promoter, whereas mutation of GTII increases gene transcription. Thus negative and positive regulatory sequences are identified within the proximal 338 bp of the GLUT2 promoter and may participate in the islet-specific expression of the gene by binding beta-cell specific trans-acting factors.