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.

Chicken BAFF--a highly conserved cytokine that mediates B cell survival.

Members of the tumor necrosis factor (TNF) family play key roles in the regulation of inflammation, immune responses and tissue homeostasis. Here we describe the identification of the chicken homologue of mammalian B cell activating factor of the TNF family (BAFF/BLyS). By searching a chicken EST database we identified two overlapping cDNA clones that code for the entire open reading frame of chicken BAFF (chBAFF), which contains a predicted transmembrane domain and a putative furin protease cleavage site like its mammalian counterparts. The amino acid identity between soluble chicken and human BAFF is 76%, considerably higher than for most other known cytokines. The chBAFF gene is most strongly expressed in the bursa of Fabricius. Soluble recombinant chBAFF produced by human 293T cells interacted with the mammalian cell-surface receptors TACI, BCMA and BAFF-R. It bound to chicken B cells, but not to other lymphocytes, and it promoted the survival of splenic chicken B cells in culture. Furthermore, bacterially expressed chBAFF induced the selective expansion of B cells in the spleen and cecal tonsils when administered to young chicks. Our results suggest that like its mammalian counterpart, chBAFF plays an important role in survival and/or proliferation of chicken B cells.

RsmY, a small regulatory RNA, is required in concert with RsmZ for GacA-dependent expression of biocontrol traits in Pseudomonas fluorescens CHA0.

In the plant-beneficial soil bacterium and biocontrol model organism Pseudomonas fluorescens CHA0, the GacS/GacA two-component system upregulates the production of biocontrol factors, i.e. antifungal secondary metabolites and extracellular enzymes, under conditions of slow, non-exponential growth. When activated, the GacS/GacA system promotes the transcription of a small regulatory RNA (RsmZ), which sequesters the small RNA-binding protein RsmA, a translational regulator of genes involved in biocontrol. The gene for a second GacA-regulated small RNA (RsmY) was detected in silico in various pseudomonads, and was cloned from strain CHA0. RsmY, like RsmZ, contains several characteristic GGA motifs. The rsmY gene was expressed in strain CHA0 as a 118 nt transcript which was most abundant in stationary phase, as revealed by Northern blot and transcriptional fusion analysis. Transcription of rsmY was enhanced by the addition of the strain's own supernatant extract containing a quorum-sensing signal and was abolished in gacS or gacA mutants. An rsmA mutation led to reduced rsmY expression, via a gacA-independent mechanism. Overexpression of rsmY restored the expression of target genes (hcnA, aprA) to gacS or gacA mutants. Whereas mutants deleted for either the rsmY or the rsmZ structural gene were not significantly altered in the synthesis of extracellular products (hydrogen cyanide, 2,4-diacetylphloroglucinol, exoprotease), an rsmY rsmZ double mutant was strongly impaired in this production and in its biocontrol properties in a cucumber-Pythium ultimum microcosm. Mobility shift assays demonstrated that multiple molecules of RsmA bound specifically to RsmY and RsmZ RNAs. In conclusion, two small, untranslated RNAs, RsmY and RsmZ, are key factors that relieve RsmA-mediated regulation of secondary metabolism and biocontrol traits in the GacS/GacA cascade of strain CHA0.

Functional interactions between the estrogen receptor and the transcription activator Sp1 regulate the estrogen-dependent transcriptional activity of the vitellogenin A1 io promoter.

Two distinct, TATA box-containing promoters regulate the transcriptional activity of the Xenopus vitellogenin A1 gene. These two promoters are of different strength and are separated by 1.8 kilobase pairs of untranslated sequence. Estrogen receptor (ER) and its ligand, 17beta-estradiol, induce the activity of both promoters. The estrogen response elements (EREs) are located proximal to the downstream i promoter while no ERE-like sequences have been identified in the vicinity of the upstream io promoter. We show here, that transcriptional activity of the upstream io promoter is Sp1-dependent. Moreover, we demonstrate that estrogen inducibility of the io promoter results from functional interactions between the io bound Sp1 and the ER bound at the proximity of i. Functional interactions between Sp1 and ER do not require the presence of a TATA box for transcriptional activation, as is demonstrated using the acyl-CoA oxidase promoter. The relative positions that ER and Sp1 occupy with respect to the initiation site determines whether these two transcription activators can synergize for transcription initiation.

Mapping of the genes coding for the two major vaccinia virus core polypeptides.

We have mapped the genes coding for two major structural polypeptides of the vaccinia virus core by hybrid selection and transcriptional mapping. First, RNA was selected by hybridization to restriction fragments of the vaccinia virus genome, translated in vitro and the products were immunoprecipitated with antibodies against the two polypeptides. This approach allowed us to map the genes to the left hand end of the largest Hind III restriction fragment of 50 kilobase pairs. Second, transcriptional mapping of this region of the genome revealed the presence of the two expected RNAs. Both RNAs are transcribed from the leftward reading strand and the 5'-ends of the genes are separated by about 7.5 kilobase pairs of DNA. Thus, two genes encoding structural polypeptides with a similar location in the vaccinia virus particle are clustered at approximately 105 kilobase pairs from the left hand end of the 180 kilobase pair vaccinia virus genome.

Landscape of transcription in human cells.

Eukaryotic cells make many types of primary and processed RNAs that are found either in specific subcellular compartments or throughout the cells. A complete catalogue of these RNAs is not yet available and their characteristic subcellular localizations are also poorly understood. Because RNA represents the direct output of the genetic information encoded by genomes and a significant proportion of a cell's regulatory capabilities are focused on its synthesis, processing, transport, modification and translation, the generation of such a catalogue is crucial for understanding genome function. Here we report evidence that three-quarters of the human genome is capable of being transcribed, as well as observations about the range and levels of expression, localization, processing fates, regulatory regions and modifications of almost all currently annotated and thousands of previously unannotated RNAs. These observations, taken together, prompt a redefinition of the concept of a gene.

Signal search analysis server.

Signal search analysis is a general method to discover and characterize sequence motifs that are positionally correlated with a functional site (e.g. a transcription or translation start site). The method has played an instrumental role in the analysis of eukaryotic promoter elements. The signal search analysis server provides access to four different computer programs as well as to a large number of precompiled functional site collections. The programs offered allow: (i) the identification of non-random sequence regions under evolutionary constraint; (ii) the detection of consensus sequence-based motifs that are over- or under-represented at a particular distance from a functional site; (iii) the analysis of the positional distribution of a consensus sequence- or weight matrix-based sequence motif around a functional site; and (iv) the optimization of a weight matrix description of a locally over-represented sequence motif. These programs can be accessed at: http://www.isrec.isb-sib.ch/ssa/.

HTPSELEX--a database of high-throughput SELEX libraries for transcription factor binding sites.

HTPSELEX is a public database providing access to primary and derived data from high-throughput SELEX experiments aimed at characterizing the binding specificity of transcription factors. The resource is primarily intended to serve computational biologists interested in building models of transcription factor binding sites from large sets of binding sequences. The guiding principle is to make available all information that is relevant for this purpose. For each experiment, we try to provide accurate information about the protein material used, details of the wet lab protocol, an archive of sequencing trace files, assembled clone sequences (concatemers) and complete sets of in vitro selected protein-binding tags. In addition, we offer in-house derived binding sites models. HTPSELEX also offers reasonably large SELEX libraries obtained with conventional low-throughput protocols. The FTP site contains the trace archives and database flatfiles. The web server offers user-friendly interfaces for viewing individual entries and quality-controlled download of SELEX sequence libraries according to a user-defined sequencing quality threshold. HTPSELEX is available from ftp://ftp.isrec.isb-sib.ch/pub/databases/htpselex/ and http://www.isrec.isb-sib.ch/htpselex.

A user's guide to the encyclopedia of DNA elements (ENCODE).

The mission of the Encyclopedia of DNA Elements (ENCODE) Project is to enable the scientific and medical communities to interpret the human genome sequence and apply it to understand human biology and improve health. The ENCODE Consortium is integrating multiple technologies and approaches in a collective effort to discover and define the functional elements encoded in the human genome, including genes, transcripts, and transcriptional regulatory regions, together with their attendant chromatin states and DNA methylation patterns. In the process, standards to ensure high-quality data have been implemented, and novel algorithms have been developed to facilitate analysis. Data and derived results are made available through a freely accessible database. Here we provide an overview of the project and the resources it is generating and illustrate the application of ENCODE data to interpret the human genome.

RecA-mediated strand exchange traverses substitutional heterologies more easily than deletions or insertions.

RecA protein in bacteria and its eukaryotic homolog Rad51 protein are responsible for initiation of strand exchange between homologous DNA molecules. This process is crucial for homologous recombination, the repair of certain types of DNA damage and for the reinitiation of DNA replication on collapsed replication forks. We show here, using two different types of in vitro assays, that in the absence of ATP hydrolysis RecA-mediated strand exchange traverses small substitutional heterologies between the interacting DNAs, whereas small deletions or insertions block the ongoing strand exchange. We discuss evolutionary implications of RecA selectivity against insertions and deletions and propose a molecular mechanism by which RecA can exert this selectivity.

Vaccinia virus produces late mRNAs by discontinuous synthesis.

We describe the unusual structure of a vaccinia virus late mRNA. In these molecules, the protein-coding sequences of a major late structural polypeptide are preceded by long leader RNAs, which in some cases are thousands of nucleotides long. These sequences map to different regions of the viral genome and in one instance are separated from the late gene by more than 100 kb of DNA. Moreover, the leader sequences map either upstream or downstream of the late gene, are transcribed from either DNA strand, and are fused to the late gene coding sequence via a poly(A) stretch. This demonstrates that vaccinia virus produces late mRNAs by tagging the protein-coding sequences onto the 3' end of other RNAs.

Peroxisome proliferator-activated receptors: finding the orphan a home.

The peroxisome proliferator-activated receptor (PPAR) is a member of the steroid hormone receptor superfamily and is activated by a variety of fibrate hypolipidaemic drugs and non-genotoxic rodent hepatocarcinogens that are collectively termed peroxisome proliferators. A key marker of peroxisome proliferator action is the peroxisomal enzyme acyl CoA oxidase, which is elevated about ten fold in the livers of treated rodents. Additional peroxisome proliferator responsive genes include other peroxisomal beta-oxidation enzymes and members of the cytochrome P450 IVA family. A peroxisome proliferator response element (PPRE), consisting of an almost perfect direct repeat of the sequence TGACCT spaced by a single base pair, has been identified in the upstream regulatory sequences of each of these genes. The retinoid X receptor (RXR) forms a heterodimer with PPAR and binds to the PPRE. Furthermore, the RXR ligand, 9-cis retinoic acid, enhances PPAR action. Retinoids may therefore modulate the action of peroxisome proliferators and PPAR may interfere with retinoid action, perhaps providing one mechanism to explain the toxicity of peroxisome proliferators. Interestingly, a variety of fatty acids can activate PPAR supporting the suggestion that fatty acids, or their acyl CoA derivatives, may be the natural ligands of PPAR and that the physiological role of PPAR is to regulate fatty acid homeostasis. Taken together, the discovery of PPAR has opened up new opportunities in understanding how lipid homeostasis is regulated, how the fibrate hypolipidaemic drugs may act and should lead to improvements in the assessment of human risk from peroxisome proliferators based upon a better understanding of their mechanism of action.

DNA binding properties of peroxisome proliferator-activated receptor subtypes on various natural peroxisome proliferator response elements. Importance of the 5'-flanking region.

The three subtypes of the peroxisome proliferator-activated receptors (PPARalpha, beta/delta, and gamma) form heterodimers with the 9-cis-retinoic acid receptor (RXR) and bind to a common consensus response element, which consists of a direct repeat of two hexanucleotides spaced by one nucleotide (DR1). As a first step toward understanding the molecular mechanisms determining PPAR subtype specificity, we evaluated by electrophoretic mobility shift assays the binding properties of the three PPAR subtypes, in association with either RXRalpha or RXRgamma, on 16 natural PPAR response elements (PPREs). The main results are as follows. (i) PPARgamma in combination with either RXRalpha or RXRgamma binds more strongly than PPARalpha or PPARbeta to all natural PPREs tested. (ii) The binding of PPAR to strong elements is reinforced if the heterodimerization partner is RXRgamma. In contrast, weak elements favor RXRalpha as heterodimerization partner. (iii) The ordering of the 16 natural PPREs from strong to weak elements does not depend on the core DR1 sequence, which has a relatively uniform degree of conservation, but correlates with the number of identities of the 5'-flanking nucleotides with respect to a consensus element. This 5'-flanking sequence is essential for PPARalpha binding and thus contributes to subtype specificity. As a demonstration of this, the PPARgamma-specific element ARE6 PPRE is able to bind PPARalpha only if its 5'-flanking region is exchanged with that of the more promiscuous HMG PPRE.

Transcript profiling suggests that differential metabolic adaptation of mice to a high fat diet is associated with changes in liver to muscle lipid fluxes.

Genetically homogenous C57Bl/6 mice display differential metabolic adaptation when fed a high fat diet for 9 months. Most become obese and diabetic, but a significant fraction remains lean and diabetic or lean and non-diabetic. Here, we performed microarray analysis of "metabolic" transcripts expressed in liver and hindlimb muscles to evaluate: (i) whether expressed transcript patterns could indicate changes in metabolic pathways associated with the different phenotypes, (ii) how these changes differed from the early metabolic adaptation to short term high fat feeding, and (iii) whether gene classifiers could be established that were characteristic of each metabolic phenotype. Our data indicate that obesity/diabetes was associated with preserved hepatic lipogenic gene expression and increased plasma levels of very low density lipoprotein and, in muscle, with an increase in lipoprotein lipase gene expression. This suggests increased muscle fatty acid uptake, which may favor insulin resistance. In contrast, the lean mice showed a strong reduction in the expression of hepatic lipogenic genes, in particular of Scd-1, a gene linked to sensitivity to diet-induced obesity; the lean and non-diabetic mice presented an additional increased expression of eNos in liver. After 1 week of high fat feeding the liver gene expression pattern was distinct from that seen at 9 months in any of the three mouse groups, thus indicating progressive establishment of the different phenotypes. Strikingly, development of the obese phenotype involved re-expression of Scd-1 and other lipogenic genes. Finally, gene classifiers could be established that were characteristic of each metabolic phenotype. Together, these data suggest that epigenetic mechanisms influence gene expression patterns and metabolic fates.

How do thyroid hormone receptors bind to structurally diverse response elements?

Three classes of thyroid hormone response elements have been described. They are composed of two half-sites arranged either as a palindromic, a direct repeat or as an inverted palindromic array. Receptor homodimers as well as heterodimers can bind to all three types of response element. While the ligand binding domain of the receptors provides the major dimerization surface, asymmetric contacts between the DNA binding domains are necessary for binding to a direct repeat. Moreover, some recent findings suggest that in TR, compared to RXR, the ligand binding domain has a 180 degrees rotation with respect to the DNA binding domain. This feature could explain the preferential binding of the RXR-TR heterodimer to the direct repeat response element, in which RXR exclusively binds the 5' half-site, and of the TR homodimer to the inverted palindrome response element.