Excess information at bacteriophage T7 genomic promoters detected by a random cloning technique.

In our previous analysis of the information at binding sites on nucleic acids, we found that most of the sites examined contain the amount of information expected from their frequency in the genome. The sequences at bacteriophage T7 promoters are an exception, because they are far more conserved (35 bits of information content) than should be necessary to distinguish them from the background of the Escherichia coli genome (17 bits). To determine the information actually used by the T7 RNA polymerase, promoters were chemically synthesized with many variations and those that function well in an in vivo assay were sequenced. Our analysis shows that the polymerase uses 18 bits of information, so the sequences at phage genomic promoters have significantly more information than the polymerase needs. The excess may represent the binding site of another protein.     

U2AF binding selects for the high conservation of the C. elegans 3' splice site

Caenorhabditis elegans is unusual among animals in having a highly conserved octamer sequence at the 3' splice site: UUUU CAG/R. This sequence can bind to the essential heterodimeric splicing factor U2AF, with U2AF65 contacting the U tract and U2AF35 contacting the splice site itself (AG/R). Here we demonstrate a strong correspondence between binding to U2AF of RNA oligonucleotides with variant octamer sequences and the frequency with which such variations occur in splice sites. C. elegans U2AF has a strong preference for the octamer sequence and exerts much of the pressure for 3' splice sites to have the precise UUUUCAG/R sequence. At two positions the splice site has a very strong preference for U even though alternative bases can also bind tightly to U2AF, suggesting that evolution can select against sequences that may have a relatively modest reduction in binding. Although pyrimidines are frequently present at the first base in the exon, U2AF has a very strong bias against them, arguing there is a mechanism to compensate for weakened U2AF binding at this position. Finally, the C in the consensus sequence must remain adjacent to the AG/R rather than to the stretch of U's, suggesting this C is recognized by U2AF35.     

A genome-wide analysis of Cas9 binding specificity using ChIP-seq and targeted sequence capture

Clustered regularly interspaced short palindromic repeat (CRISPR) RNA-guided nucleases have gathered considerable excitement as a tool for genome engineering. However, questions remain about the specificity of target site recognition. Cleavage specificity is typically evaluated by low throughput assays (T7 endonuclease I assay, target amplification followed by high-throughput sequencing), which are limited to a subset of potential off-target sites. Here, we used ChIP-seq to examine genome-wide CRISPR binding specificity at gRNA-specific and gRNA-independent sites for two guide RNAs. RNA-guided Cas9 binding was highly specific to the target site while off-target binding occurred at much lower intensities. Cas9-bound regions were highly enriched in NGG sites, a sequence required for target site recognition by Streptococcus pyogenes Cas9. To determine the relationship between Cas9 binding and endonuclease activity, we applied targeted sequence capture, which allowed us to survey 1200 genomic loci simultaneously including potential off-target sites identified by ChIP-seq and by computational prediction. A high frequency of indels was observed at both target sites and one off-target site, while no cleavage activity could be detected at other ChIP-bound regions. Our results confirm the high-specificity of CRISPR endonucleases and demonstrate that sequence capture can be used as a high-throughput genome-wide approach to identify off-target activity.     

In vivo footprinting studies suggest a role for chromatin in transcription of the human 7SK gene

Mutation and deletion analyses of mammalian class III small nuclear RNA genes transcribed by RNA polymerase (pol) III have defined three functional promoter elements: a distal sequence element (DSE) at around -220, a proximal sequence element (PSE) at around -60 and a TATA box at around -30. Although binding studies have identified factors that bind to these sites in vitro, it is not known exactly how proteins interact with the promoters of these genes in vivo. In this study, we have used dimethyl sulphate and DNase I treatment of HeLa cells and nuclei, respectively, followed by linker-mediated polymerase chain reaction, to obtain in vivo footprints of proteins binding to the promoter of the Pol III-transcribed 7SK gene. Our results show that most of the characterised promoter elements of this gene are protected in vivo in these cells, and the pattern of DNase I protection suggests that a nucleosome lies between the DSE and the PSE. Methylation protection was also seen upstream of the DSE over a sequence corresponding to the binding site of a POZ domain-containing protein, ZID, which interacts with components of histone deacetylase complexes. These findings suggest that chromatin structure plays a role in the cascade of protein-DNA interactions that regulate expression of this pol III-transcribed gene.     

Role of alpha-transinducing factor (VP16) in the induction of alpha genes within the context of viral genomes.

In herpes simplex virus 1, the five alpha genes are induced by alpha-transinducing factor (alpha TIF; VP16), a virion protein, acting in concert with Oct-1 and other cellular proteins on a cis-acting site in the promoter domain of alpha genes. Because alpha TIF is an essential virion protein, its function as an inducer can best be evaluated only by mutating the cis-acting site. Earlier we reported on a series of 17 mutations in and around the cis-acting site of a 275-bp alpha 27 promoter fused to a reporter gene and recombined into the viral genome. These recombinant viruses were tested in Vero cells in the presence of cycloheximide, and we demonstrated that mutations in the sequence required for Oct-1 binding abolished transactivation whereas mutations in the alpha TIF-dependent GARAT sequence decreased but did not abolish transactivation. We now report that (i) in limited-passage human embryonic lung cells, alpha gene expression from promoters mutated in the GARAT sequences is often higher and more variable than in Vero cells, (ii) in the absence of cycloheximide, the mutant viruses show less significant impairment of reporter gene expression, (iii) Oct-1 can bind either to the overlapping octamer element or to various TAATGARAT sequences with differing degrees of binding strength and these relative binding levels correlate well with levels of gene expression observed in infected cells, (iv) in the cis-acting site upstream of the alpha 4 gene, no degenerate overlapping Oct-1 sequence exists, and therefore in this instance Oct-1 must be binding directly to the TAATGARAT sequence, (v) extension of the alpha 27 promoter by an additional 1,334 bp results in much higher expression of the reporter gene as a result of additional upstream cis-acting sites, and (vi) obliteration of the most proximal Oct-1 binding element within the 275-bp promoter dramatically reduces gene expression even in the presence of the additional upstream cis-acting sites.     

Development and refinement of pregnane X receptor (PXR) DNA binding site model using information theory: insights into PXR-mediated gene regulation

The pregnane X receptor (PXR) acts as a receptor to induce gene expression in response to structurally diverse xenobiotics through binding as a heterodimer with the 9-cis retinoic acid receptor (RXR) to enhancers in target gene promoters. We identified and estimated the affinities of novel PXR/RXR binding sites in regulated genes and additional genomic targets of PXR with an information theory-based model of the PXR/RXR binding site. Our initial PXR/RXR model, the result of the alignment of 15 previously characterized binding sites, was used to scan the promoters of known PXR target genes. Sites from these genes, with information contents of >8 bits bound by PXR/RXR in vitro, were used to revise the information weight matrix; this procedure was repeated by screening for progressively weaker binding sites. After three iterations of refinement, the model was based on 48 validated PXR/RXR binding sites and has an average information content (Rsequence) of 14.43 +/- 3.21 bits. A scan of the human genome predicted novel PXR/RXR binding sites in the promoters of UGT1A3 (19.78 bits at -8040 and 16.37 bits at -6930) and UGT1A6 (12.74 bits at -9216), both of which were identified previously as targets for PXR. These sites were subsequently demonstrated to specifically bind PXR/RXR in competition electrophoretic mobility shift assays. A strong PXR site was also predicted upstream of the CASP10 gene (18.69 bits at -7872) and was validated by binding studies and reporter assays as a PXR responsive element. This suggests that the PXR-mediated response extends beyond genes involved in drug biotransformation and transport.     

Nucleotide sequence from the genetic left end of bacteriophage T7 DNA to the beginning of gene 4

The nucleotide sequence running from the genetic left end of bacteriophage T7 DNA to within the coding sequence of gene 4 is given, except for the internal coding sequence for the gene 1 protein, which has been determined elsewhere. The sequence presented contains nucleotides 1 to 3342 and 5654 to 12,100 of the approximately 40,000 base-pairs of T7 DNA. This sequence includes: the three strong early promoters and the termination site for Escherichia coli RNA polymerase: eight promoter sites for T7 RNA polymerase; six RNAase III cleavage sites; the primary origin of replication of T7 DNA; the complete coding sequences for 13 previously known T7 proteins, including the anti-restriction protein, protein kinase, DNA ligase, the gene 2 inhibitor of E. coli RNA polymerase, single-strand DNA binding protein, the gene 3 endonuclease, and lysozyme (which is actually an N-acetylmuramyl-l-alanine amidase); the complete coding sequences for eight potential new T7-coded proteins; and two apparently independent initiation sites that produce overlapping polypeptide chains of gene 4 primase. More than 86% of the first 12,100 base-pairs of T7 DNA appear to be devoted to specifying amino acid sequences for T7 proteins, and the arrangement of coding sequences and other genetic elements is very efficient. There is little overlap between coding sequences for different proteins, but junctions between adjacent coding sequences are typically close, the termination codon for one protein often overlapping the initiation codon for the next. For almost half of the potential T7 proteins, the sequence in the messenger RNA that can interact with 16 S ribosomal RNA in initiation of protein synthesis is part of the coding sequence for the preceding protein. The longest non-coding region, about 900 base-pairs, is at the left end of the DNA. The right half of this region contains the strong early promoters for E. coli RNA polymerase and the first RNAase III cleavage site. The left end contains the terminal repetition (nucleotides 1 to 160), followed by a striking array of repeated sequences (nucleotides 175 to 340) that might have some role in packaging the DNA into phage particles, and an A · T-rich region (nucleotides 356 to 492) that contains a promoter for T7 RNA polymerase, and which might function as a replication origin.     

Compilation and analysis of Bacillus subtilis sigma A-dependent promoter sequences: evidence for extended contact between RNA polymerase and upstream promoter DNA.

Sequence analysis of 236 promoters recognized by the Bacillus subtilis sigma A-RNA polymerase reveals an extended promoter structure. The most highly conserved bases include the -35 and -10 hexanucleotide core elements and a TG dinucleotide at position -15, -14. In addition, several weakly conserved A and T residues are present upstream of the -35 region. Analysis of dinucleotide composition reveals A2- and T2-rich sequences in the upstream promoter region (-36 to -70) which are phased with the DNA helix: An tracts are common near -43, -54 and -65; Tn tracts predominate at the intervening positions. When compared with larger regions of the genome, upstream promoter regions have an excess of An and Tn sequences for n > 4. These data indicate that an RNA polymerase binding site affects DNA sequence as far upstream as -70. This sequence conservation is discussed in light of recent evidence that the alpha subunits of the polymerase core bind DNA and that the promoter may wrap around RNA polymerase.     

A nuclear protein is required for thyroid hormone receptor binding to an inhibitory half-site in the epidermal growth factor receptor promoter.

The epidermal growth factor (EGF) receptor (EGFR) promoter is negatively regulated by thyroid hormone and retinoic acid. This regulation can be mapped to a 36-basepair GC-rich region of the promoter (EGFR P/E) that functions autonomously as a promoter and an enhancer when placed in front of the thymidine kinase gene TATA element. Direct high affinity binding of the thyroid hormone receptor (T3R) to this element requires a nuclear protein. Through ion exchange chromatography and gel filtration of HeLa nuclear extract, this activity was identified as a protein of approximately 67 kilodaltons. This protein did not bind to DNA alone, but greatly augmented T3R binding to the EGFR P/E sequence in gel mobility shift and DNA precipitation assays. When combined with the T3R auxillary protein (TRAP), the T3R migrated as a larger complex on the DNA. Chemical cross-linking identified this complex as a heterodimer between T3R and TRAP. T3R-TRAP binds to a 7-basepair site in the EGFR P/E (GGGACTC) that has weak homology to a consensus thyroid response element half-site. Thus, on this element, T3R-TRAP heterodimers contact the DNA primarily on a single site that comprises an inhibitory thyroid response element.