Instruct: Given a scientific claim, retrieve relevant passages that support or refute the claim.
Query: Sequence conservation in gene regulatory regions between species is a poor predictor of conserved functionality.
The mammalian radiation has corresponded with rapid changes in noncoding regions of the genome, but we lack a comprehensive understanding of regulatory evolution in mammals. Here, we track the evolution of promoters and enhancers active in liver across 20 mammalian species from six diverse orders by profiling genomic enrichment of H3K27 acetylation and H3K4 trimethylation. We report that rapid evolution of enhancers is a universal feature of mammalian genomes. Most of the recently evolved enhancers arise from ancestral DNA exaptation, rather than lineage-specific expansions of repeat elements. In contrast, almost all liver promoters are partially or fully conserved across these species. Our data further reveal that recently evolved enhancers can be associated with genes under positive selection, demonstrating the power of this approach for annotating regulatory adaptations in genomic sequences. These results provide important insight into the functional genetics underpinning mammalian regulatory evolution.
Mammalian genomes are extensively transcribed outside the borders of protein-coding genes. Genome-wide studies recently demonstrated that cis-regulatory genomic elements implicated in transcriptional control, such as enhancers and locus-control regions, represent major sites of extragenic noncoding transcription. Enhancer-templated transcripts provide a quantitatively small contribution to the total amount of cellular nonribosomal RNA; nevertheless, the possibility that enhancer transcription and the resulting enhancer RNAs may, in some cases, have functional roles, rather than represent mere transcriptional noise at accessible genomic regions, is supported by an increasing amount of experimental data. In this article we review the current knowledge on enhancer transcription and its functional implications.
The prediction of novel pre-microRNA (miRNA) from genomic sequence has received considerable attention recently. However, the majority of studies have focused on the human genome. Previous studies have demonstrated that sensitivity (correctly detecting true miRNA) is sustained when human-trained methods are applied to other species, however they have failed to report the dramatic drop in specificity (the ability to correctly reject non-miRNA sequences) in non-human genomes. Considering the ratio of true miRNA sequences to pseudo-miRNA sequences is on the order of 1:1000, such low specificity prevents the application of most existing tools to non-human genomes, as the number of false positives overwhelms the true predictions. We here introduce a framework (SMIRP) for creating species-specific miRNA prediction systems, leveraging sequence conservation and phylogenetic distance information. Substantial improvements in specificity and precision are obtained for four non-human test species when our framework is applied to three different prediction systems representing two types of classifiers (support vector machine and Random Forest), based on three different feature sets, with both human-specific and taxon-wide training data. The SMIRP framework is potentially applicable to all miRNA prediction systems and we expect substantial improvement in precision and specificity, while sustaining sensitivity, independent of the machine learning technique chosen.
Genome function, replication, integrity, and propagation rely on the dynamic structural organization of chromosomes during the cell cycle. Genome folding in interphase provides regulatory segmentation for appropriate transcriptional control, facilitates ordered genome replication, and contributes to genome integrity by limiting illegitimate recombination. Here, we review recent high-resolution chromosome conformation capture and functional studies that have informed models of the spatial and regulatory compartmentalization of mammalian genomes, and discuss mechanistic models for how CTCF and cohesin control the functional architecture of mammalian chromosomes.
Detection of new genomic control elements is critical in understanding transcriptional regulatory networks in their entirety. We studied the genome-wide binding locations of three key regulatory proteins (POU5F1, also known as OCT4; NANOG; and CTCF) in human and mouse embryonic stem cells. In contrast to CTCF, we found that the binding profiles of OCT4 and NANOG are markedly different, with only ∼5% of the regions being homologously occupied. We show that transposable elements contributed up to 25% of the bound sites in humans and mice and have wired new genes into the core regulatory network of embryonic stem cells. These data indicate that species-specific transposable elements have substantially altered the transcriptional circuitry of pluripotent stem cells.
Enhancers control the correct temporal and cell-type-specific activation of gene expression in multicellular eukaryotes. Knowing their properties, regulatory activity and targets is crucial to understand the regulation of differentiation and homeostasis. Here we use the FANTOM5 panel of samples, covering the majority of human tissues and cell types, to produce an atlas of active, in vivo-transcribed enhancers. We show that enhancers share properties with CpG-poor messenger RNA promoters but produce bidirectional, exosome-sensitive, relatively short unspliced RNAs, the generation of which is strongly related to enhancer activity. The atlas is used to compare regulatory programs between different cells at unprecedented depth, to identify disease-associated regulatory single nucleotide polymorphisms, and to classify cell-type-specific and ubiquitous enhancers. We further explore the utility of enhancer redundancy, which explains gene expression strength rather than expression patterns. The online FANTOM5 enhancer atlas represents a unique resource for studies on cell-type-specific enhancers and gene regulation.
Synthetic gene expression is highly sensitive to intragenic compositional context (promoter structure, spacing regions between promoter and coding sequences, and ribosome binding sites). However, much less is known about the effects of intergenic compositional context (spatial arrangement and orientation of entire genes on DNA) on expression levels in synthetic gene networks. We compare expression of induced genes arranged in convergent, divergent, or tandem orientations. Induction of convergent genes yielded up to 400% higher expression, greater ultrasensitivity, and dynamic range than divergent- or tandem-oriented genes. Orientation affects gene expression whether one or both genes are induced. We postulate that transcriptional interference in divergent and tandem genes, mediated by supercoiling, can explain differences in expression and validate this hypothesis through modeling and in vitro supercoiling relaxation experiments. Treatment with gyrase abrogated intergenic context effects, bringing expression levels within 30% of each other. We rebuilt the toggle switch with convergent genes, taking advantage of supercoiling effects to improve threshold detection and switch stability.
Chromatin insulators are DNA elements that regulate the level of gene expression either by preventing gene silencing through the maintenance of heterochromatin boundaries or by preventing gene activation by blocking interactions between enhancers and promoters. CCCTC-binding factor (CTCF), a ubiquitously expressed 11-zinc-finger DNA-binding protein, is the only protein implicated in the establishment of insulators in vertebrates. While CTCF has been implicated in diverse regulatory functions, CTCF has only been studied in a limited number of cell types across human genome. Thus, it is not clear whether the identified cell type-specific differences in CTCF-binding sites are functionally significant. Here, we identify and characterize cell type-specific and ubiquitous CTCF-binding sites in the human genome across 38 cell types designated by the Encyclopedia of DNA Elements (ENCODE) consortium. These cell type-specific and ubiquitous CTCF-binding sites show uniquely versatile transcriptional functions and characteristic chromatin features. In addition, we confirm the insulator barrier function of CTCF-binding and explore the novel function of CTCF in DNA replication. These results represent a critical step toward the comprehensive and systematic understanding of CTCF-dependent insulators and their versatile roles in the human genome.
Gene expression is controlled by the complex interaction of transcription factors binding to promoters and other regulatory DNA elements. One common characteristic of the genomic regions associated with regulatory proteins is a pronounced sensitivity to DNase I digestion. We generated genome-wide high-resolution maps of DNase I hypersensitive (DH) sites from both seedling and callus tissues of rice (Oryza sativa). Approximately 25% of the DH sites from both tissues were found in putative promoters, indicating that the vast majority of the gene regulatory elements in rice are not located in promoter regions. We found 58% more DH sites in the callus than in the seedling. For DH sites detected in both the seedling and callus, 31% displayed significantly different levels of DNase I sensitivity within the two tissues. Genes that are differentially expressed in the seedling and callus were frequently associated with DH sites in both tissues. The DNA sequences contained within the DH sites were hypomethylated, consistent with what is known about active gene regulatory elements. Interestingly, tissue-specific DH sites located in the promoters showed a higher level of DNA methylation than the average DNA methylation level of all the DH sites located in the promoters. A distinct elevation of H3K27me3 was associated with intergenic DH sites. These results suggest that epigenetic modifications play a role in the dynamic changes of the numbers and DNase I sensitivity of DH sites during development.
UNLABELLED Artificially synthesized short interfering RNAs (siRNAs) are widely used in functional genomics to knock down specific target genes. One ongoing challenge is to guarantee that the siRNA does not elicit off-target effects. Initial reports suggested that siRNAs were highly sequence-specific; however, subsequent data indicates that this is not necessarily the case. It is still uncertain what level of similarity and other rules are required for an off-target effect to be observed, and scoring schemes have not been developed to look beyond simple measures such as the number of mismatches or the number of consecutive matching bases present. We created design rules for predicting the likelihood of a non-specific effect and present a web server that allows the user to check the specificity of a given siRNA in a flexible manner using a combination of methods. The server finds potential off-target matches in the corresponding RefSeq database and ranks them according to a scoring system based on experimental studies of specificity. AVAILABILITY The server is available at http://informatics-eskitis.griffith.edu.au/SpecificityServer.
The most highly conserved noncoding elements (HCNEs) in mammalian genomes cluster within regions enriched for genes encoding developmentally important transcription factors (TFs). This suggests that HCNE-rich regions may contain key regulatory controls involved in development. We explored this by examining histone methylation in mouse embryonic stem (ES) cells across 56 large HCNE-rich loci. We identified a specific modification pattern, termed "bivalent domains," consisting of large regions of H3 lysine 27 methylation harboring smaller regions of H3 lysine 4 methylation. Bivalent domains tend to coincide with TF genes expressed at low levels. We propose that bivalent domains silence developmental genes in ES cells while keeping them poised for activation. We also found striking correspondences between genome sequence and histone methylation in ES cells, which become notably weaker in differentiated cells. These results highlight the importance of DNA sequence in defining the initial epigenetic landscape and suggest a novel chromatin-based mechanism for maintaining pluripotency.
Transcription regulation in higher eukaryotes is controlled by regulatory elements such as enhancers that are recognized by transcription factors. In many cases regulatory elements can be located at distances up to several megabases from their target genes. Recent evidence shows that long-range control of gene expression can be mediated through direct physical interactions between genes and these regulatory elements. Such looping interactions can be detected using the chromosome conformation capture (3C) methodology. Although 3C is experimentally straightforward, to draw meaningful conclusions one must carefully design 3C experiments and implement the conscientious use of controls. The general guidelines presented here should help experimental design and minimize misinterpretation of 3C experiments.
Despite being one of the first eukaryotic transcriptional regulatory elements identified, the sequence of a native TATA box and its significance remain elusive. Applying criteria associated with TATA boxes we queried several Saccharomyces genomes and arrived at the consensus TATA(A/T)A(A/T)(A/G). Approximately 20% of yeast genes contain a TATA box. Strikingly, TATA box-containing genes are associated with responses to stress, are highly regulated, and preferentially utilize SAGA rather than TFIID when compared to TATA-less promoters. Transcriptional regulation in yeast appears to be mechanistically bipolar, possibly reflecting a need to balance inducible stress-related responses with constitutive housekeeping functions.
Studies of replication timing provide a handle into previously impenetrable higher-order levels of chromosome organization and their plasticity during development. Although mechanisms regulating replication timing are not clear, novel genome-wide studies provide a thorough survey of the extent to which replication timing is regulated during most of the early cell fate transitions in mammals, revealing coordinated changes of a defined set of 400–800 kb chromosomal segments that involve at least half the genome. Furthermore, changes in replication time are linked to changes in sub-nuclear organization and domain-wide transcriptional potential, and tissue-specific replication timing profiles are conserved from mouse to human, suggesting that the program has developmental significance. Hence, these studies have provided a solid foundation for linking megabase level chromosome structure to function, and suggest a central role for replication in domain-level genome organization.
We report a database of genes responsive to the Myc oncogenic transcription factor. The database Myc Target Gene prioritizes candidate target genes according to experimental evidence and clusters responsive genes into functional groups. We coupled the prioritization of target genes with phylogenetic sequence comparisons to predict c-Myc target binding sites, which are in turn validated by chromatin immunoprecipitation assays. This database is essential for the understanding of the genetic regulatory networks underlying the genesis of cancers.
Archaea contain a variety of sequence-independent DNA binding proteins consistent with the evolution of several different, sometimes overlapping and exchangeable solutions to the problem of genome compaction. Some of these proteins undergo residue-specific post-translational lysine acetylation or methylation, hinting at analogues of the histone modifications that regulate eukaryotic chromatin structure and transcription. Archaeal transcription initiation most closely resembles the eukaryotic RNA polymerase II (RNAPII) system, but Archaea do not appear to have homologues of the multisubunit complexes that remodel eukaryotic chromatin and activate RNAPII initiation. In contrast, they have sequence-specific regulators that repress and perhaps activate archaeal transcription by mechanisms superficially similar to the bacterial paradigm of regulating promoter binding by RNAP. Repressors compete with archaeal TATA-box binding protein (TBP) and TFB for the TATA-box and TFB-recognition elements (BRE) of the archaeal promoter, or with archaeal RNAP for the site of transcription initiation. Transcript-specific regulation by repressors binding to sites of transcript initiation is consistent with such sites having very little sequence conservation. However, most Archaea have only one TBP and/or TFB that presumably must therefore bind to similar TATA-box and BRE sequences upstream of most genes. Repressors that function by competing with TBP and/or TFB binding must therefore also make additional contacts with transcript-specific regulatory sites adjacent or remote from the TATA-box/BRE region. The fate of the archaeal TBP and TFB following transcription initiation remains to be determined. Based on functional homology with their eukaryotic RNAPII-system counterparts, archaeal TBP and possibly also TFB should remain bound to the TATA-box/BRE region after transcription initiation. However, this seems unlikely as it might limit repressor competition at this site to only the first round of transcription initiation.
Evolutionary related multisubunit RNA polymerases (RNAPs) transcribe the genomes of all living organisms. Whereas the core subunits of RNAPs are universally conserved in all three domains of life-indicative of a common evolutionary descent-this only applies to one RNAP-associated transcription factor-Spt5, also known as NusG in bacteria. All other factors that aid RNAP during the transcription cycle are specific for the individual domain or only conserved between archaea and eukaryotes. Spt5 and its bacterial homologue NusG regulate gene expression in several ways by (i) modulating transcription processivity and promoter proximal pausing, (ii) coupling transcription and RNA processing or translation, and (iii) recruiting termination factors and thereby silencing laterally transferred DNA and protecting the genome against double-stranded DNA breaks. This review discusses recent discoveries that identify Spt5-like factors as evolutionary conserved nexus for the regulation and coordination of the machineries responsible for information processing in the cell.
The yeast ribosomal protein gene RPL32 of Saccharomyces cerevisiae is of particular interest for two reasons: 1) it is adjacent to another ribosomal protein gene, RP29, whose divergent transcription may be driven from the same control sequences, and 2) it appears that the splicing of its transcript is regulated by the product of the gene, ribosomal protein in L32. RPL32 has been analyzed in detail. It is essential for cell growth. Its sequence predicts L32 to be a protein of 105 amino acids, somewhat basic near the NH2 terminus, rather acidic near the COOH terminus, and homologous to ribosomal protein L30 of mammals. The reading frame has been confirmed by partial NH2-terminal analysis of L32. The nucleotide sequence also predicts an intron of 230 nucleotides, which begins with the unusual sequence GTCAGT and ends 40 nucleotides downstream of the consensus sequence TAC-TAAC. The intron has been confirmed by determination of the sequence of a cDNA clone. Transcription initiates 58 nucleotides upstream of the AUG initiation codon, and the polyadenylation site occurs 100 nucleotides downstream of the termination codon. Regulation of the transcription of ribosomal protein genes has been linked to two related consensus sequences. Analysis of the intergenic region between RP29 and RPL32 reveals three copies of these sequences. A deletion removing all three sequences reduces synthesis of a L32-LacZ fusion protein by more than 90%. Some residual activity, however, remains.
Polycomb Group (PcG) proteins mediate heritable gene silencing by modifying chromatin structure. An essential PcG complex, PRC1, compacts chromatin and inhibits chromatin remodeling. In Drosophila melanogaster, the intrinsically disordered C-terminal region of PSC (PSC-CTR) mediates these noncovalent effects on chromatin, and is essential for viability. Because the PSC-CTR sequence is poorly conserved, the significance of its effects on chromatin outside of Drosophila was unclear. The absence of folded domains also made it difficult to understand how the sequence of PSC-CTR encodes its function. To determine the mechanistic basis and extent of conservation of PSC-CTR activity, we identified 17 metazoan PSC-CTRs spanning chordates to arthropods, and examined their sequence features and biochemical properties. PSC-CTR sequences are poorly conserved, but are all highly charged and structurally disordered. We show that active PSC-CTRs--which bind DNA tightly and inhibit chromatin remodeling efficiently--are distinguished from less active ones by the absence of extended negatively charged stretches. PSC-CTR activity can be increased by dispersing its contiguous negative charge, confirming the importance of this property. Using the sequence properties defined as important for PSC-CTR activity, we predicted the presence of active PSC-CTRs in additional diverse genomes. Our analysis reveals broad conservation of PSC-CTR activity across metazoans. This conclusion could not have been determined from sequence alignments. We further find that plants that lack active PSC-CTRs instead possess a functionally analogous PcG protein, EMF1. Thus, our study suggests that a disordered domain with dispersed negative charges underlies PRC1 activity, and is conserved across metazoans and plants.
The transcriptional regulator Spx plays a key role in maintaining the redox homeostasis of Bacillus subtilis cells exposed to disulfide stress. Defects in Spx were previously shown to lead to differential expression of numerous genes but direct and indirect regulatory effects could not be distinguished. Here we identified 283 discrete chromosomal sites potentially bound by the Spx-RNA polymerase (Spx-RNAP) complex using chromatin immunoprecipitation of Spx. Three quarters of these sites were located near Sigma(A)-dependent promoters, and upon diamide treatment, the fraction of the Spx-RNAP complex increased in parallel with the number and occupancy of DNA sites. Correlation of Spx-RNAP-binding sites with gene differential expression in wild-type and Δspx strains exposed or not to diamide revealed that 144 transcription units comprising 275 genes were potentially under direct Spx regulation. Spx-controlled promoters exhibited an extended -35 box in which nucleotide composition at the -43/-44 positions strongly correlated with observed activation. In vitro transcription confirmed activation by oxidized Spx of seven newly identified promoters, of which one was also activated by reduced Spx. Our study globally characterized the Spx regulatory network, revealing its role in the basal expression of some genes and its complex interplay with other stress responses.
To characterize the sequence features surrounding the translation initiation sites on the genome of Synechocystis sp. strain 6803, the total proteins extracted from the cell were resolved by two-dimensional electrophoresis, and the amino-terminal sequences of the relatively abundant protein spots were determined. By comparison of the determined amino-terminal sequences with the nucleotide sequence of the entire genome, the translation initiation sites of a total of 72 proteins were successfully assigned on the genome. The sequence features emerged from the nucleotide sequences at and surrounding the translation initiation sites were as follows: (1) In addition to the three initiation codons, ATG, GTG, and TTG, evidence was obtained that ATT was also used as a rare initiation codon; (2) the core sequences (GAGG, GGAG and AGGA) of the Shine-Dalgarno sequence were identified in the appropriate position preceding the 35 initiation sites (48.6%); and (3) the preferential sequence surrounding the initiation codons was formulated as 5'-YY[...]R-3' where Y and R denote pyrimidine and purine nucleotides, respectively, and three dots represent the initiation codons. The result obtained would provide valuable information for improvement of the gene-finding software, and the approach used in this study should be applicable for comprehensive analysis of the expression profiles of cellular proteins.
The reconstruction of gene regulatory networks underlying cell differentiation from high-throughput gene expression and chromatin data remains a challenge. Here, we derive dynamic gene regulatory networks for human myeloid differentiation using a 5-day time series of RNA-seq and ATAC-seq data. We profile HL-60 promyelocytes differentiating into macrophages, neutrophils, monocytes, and monocyte-derived macrophages. We find a rapid response in the expression of key transcription factors and lineage markers that only regulate a subset of their targets at a given time, which is followed by chromatin accessibility changes that occur later along with further gene expression changes. We observe differences between promyelocyte- and monocyte-derived macrophages at both the transcriptional and chromatin landscape level, despite using the same differentiation stimulus, which suggest that the path taken by cells in the differentiation landscape defines their end cell state. More generally, our approach of combining neighboring time points and replicates to achieve greater sequencing depth can efficiently infer footprint-based regulatory networks from long series data.
SR proteins are well-characterized RNA binding proteins that promote exon inclusion by binding to exonic splicing enhancers (ESEs). However, it has been unclear whether regulatory rules deduced on model genes apply generally to activities of SR proteins in the cell. Here, we report global analyses of two prototypical SR proteins, SRSF1 (SF2/ASF) and SRSF2 (SC35), using splicing-sensitive arrays and CLIP-seq on mouse embryo fibroblasts (MEFs). Unexpectedly, we find that these SR proteins promote both inclusion and skipping of exons in vivo, but their binding patterns do not explain such opposite responses. Further analyses reveal that loss of one SR protein is accompanied by coordinated loss or compensatory gain in the interaction of other SR proteins at the affected exons. Therefore, specific effects on regulated splicing by one SR protein actually depend on a complex set of relationships with multiple other SR proteins in mammalian genomes.
RNases P and MRP are ribonucleoprotein complexes involved in tRNA and rRNA processing, respectively. The RNA subunits of these two enzymes are structurally related to each other and play an essential role in the enzymatic reaction. Both of the RNAs have a highly conserved helical region, P4, which is important in the catalytic reaction. We have used a bioinformatics approach based on conserved elements to computationally analyze available genomic sequences of eukaryotic organisms and have identified a large number of novel nuclear RNase P and MRP RNA genes. For MRP RNA for instance, this investigation increases the number of known sequences by a factor of three. We present secondary structure models of many of the predicted RNAs. Although all sequences are able to fold into the consensus secondary structure of P and MRP RNAs, a striking variation in size is observed, ranging from a Nosema locustae MRP RNA of 160 nt to much larger RNAs, e.g. a Plasmodium knowlesi P RNA of 696 nt. The P and MRP RNA genes appear in tandem in some protists, further emphasizing the close evolutionary relationship of these RNAs.
In the present study, the relationship between short interfering RNA (siRNA) sequence and RNA interference (RNAi) effect was extensively analyzed using 62 targets of four exogenous and two endogenous genes and three mammalian and Drosophila cells. We present the rules that may govern siRNA sequence preference and in accordance with which highly effective siRNAs essential for systematic mammalian functional genomics can be readily designed. These rules indicate that siRNAs which simultaneously satisfy all four of the following sequence conditions are capable of inducing highly effective gene silencing in mammalian cells: (i) A/U at the 5' end of the antisense strand; (ii) G/C at the 5' end of the sense strand; (iii) at least five A/U residues in the 5' terminal one-third of the antisense strand; and (iv) the absence of any GC stretch of more than 9 nt in length. siRNAs opposite in features with respect to the first three conditions give rise to little or no gene silencing in mammalian cells. Essentially the same rules for siRNA sequence preference were found applicable to DNA-based RNAi in mammalian cells and in ovo RNAi using chick embryos. In contrast to mammalian and chick cells, little siRNA sequence preference could be detected in Drosophila in vivo RNAi.