Next-generation SELEX identifies sequence and structural determinants of splicing factor binding in human pre-mRNA sequence

Many splicing factors interact with both mRNA and pre-mRNA. The identification of these interactions has been greatly improved by the development of in vivo cross-linking immunoprecipitation. However, the output carries a strong sampling bias in favor of RNPs that form on more abundant RNA species like mRNA. We have developed a novel in vitro approach for surveying binding on pre-mRNA, without cross-linking or sampling bias. Briefly, this approach entails specifically designed oligonucleotide pools that tile through a pre-mRNA sequence. The pool is then partitioned into bound and unbound fractions, which are quantified by a two-color microarray. We applied this approach to locating splicing factor binding sites in and around ∼4000 exons. We also quantified the effect of secondary structure on binding. The method is validated by the finding that U1snRNP binds at the 5′ splice site (5′ss) with a specificity that is nearly identical to the splice donor motif. In agreement with prior reports, we also show that U1snRNP appears to have some affinity for intronic G triplets that are proximal to the 5′ss. Both U1snRNP and the polypyrimidine tract binding protein (PTB) avoid exonic binding, and the PTB binding map shows increased enrichment at the polypyrimidine tract. For PTB, we confirm polypyrimidine specificity and are also able to identify structural determinants of PTB binding. We detect multiple binding motifs enriched in the PTB bound fraction of oligonucleotides. These motif combinations augment binding in vitro and are also enriched in the vicinity of exons that have been determined to be in vivo targets of PTB.     

Highly recurring sequence elements identified in eukaryotic DNAs by computer analysis are often homologous to regulatory sequences or protein binding sites.

We have used computer assisted dot matrix and oligonucleotide frequency analyses to identify highly recurring sequence elements of 7-11 base pairs in eukaryotic genes and viral DNAs. Such elements are found much more frequently than expected, often with an average spacing of a few hundred base pairs. Furthermore, the most abundant repetitive elements observed in the ovalbumin locus, the beta-globin gene cluster, the metallothionein gene and the viral genomes of SV40, polyoma, Herpes simplex-1 and Mouse Mammary Tumor Virus were sequences shown previously to be protein binding sites or sequences important for regulating gene expression. These sequences were present in both exons and introns as well as promoter regions. These observations suggest that such sequences are often highly overrepresented within the specific gene segments with which they are associated. Computer analysis of other genetic units, including viral genomes and oncogenes, has identified a number of highly recurring sequence elements that could serve similar regulatory or protein-binding functions. A model for the role of such reiterated sequence elements in DNA organization and function is presented.     

Two-dimensional combinatorial screening identifies specific 6'-acylated kanamycin A- and 6'-acylated neamine-RNA hairpin interactions

Herein, we report the RNA hairpin loops from a six-nucleotide hairpin library that bind 6'-acylated kanamycin A (1) and 6'-acylated neamine (2) identified by two-dimensional combinatorial screening (2DCS). Hairpins selected to bind 1 have K(d)'s ranging from 235 to 1035 nM, with an average K(d) of 618 nM. For 2, the selected hairpins bind with K(d)'s ranging from 135 to 2300 nM, with an average K(d) of 1010 nM. The selected RNA hairpin-ligand interactions are also specific for the ligand that they were selected to bind compared with the other arrayed ligand. For example, the mixture of hairpins selected for 1 on average bind 33-fold more tightly to 1 than to 2, while the mixtures of hairpins selected for 2 on average bind 11-fold more tightly to 2 than to 1. Secondary structure prediction of the selected sequences was completed to determine the motifs that each ligand binds, and the hairpin loop preferences for 1 and 2 were computed. For 1, the preferred hairpin loops contain an adenine separated by at least two nucleotides from a cytosine, for example, ANNCNN (two-tailed p-value = 0.0010) and ANNNCN (two-tailed p-value <0.0001). For 2, the preferred hairpin loops contain both 5'GC and 5'CG steps (two-tailed p-value <0.0001). These results expand the information available on the RNA hairpin loops that bind small molecules and could prove useful for targeting RNA.     

Selective formation of Sed5p-containing SNARE complexes is mediated by combinatorial binding interactions

Sed5p is the only syntaxin family member required for protein transport through the yeast Golgi and it is known to bind up to nine other soluble N-ethylmaleimide-sensitive factor attachment receptor (SNARE) proteins in vivo. We describe in vitro binding experiments in which we identify ternary and quaternary Sed5p-containing SNARE complexes. The formation of SNARE complexes among these endoplasmic reticulum- and Golgi-localized proteins requires Sed5p and is syntaxin-selective. In addition, Sed5p-containing SNARE complexes form selectively and this selectivity is mediated by Sed5p-containing intermediates that discriminate among subsequent binding partners. Although many of these SNAREs have overlapping distributions in vivo, the SNAREs that form complexes with Sed5p in vitro reflect their functionally distinct locales. Although SNARE-SNARE interactions are promiscuous and a single SNARE protein is often found in more than one complex, both the biochemical as well as genetic analyses reported here suggest that this is not a result of nonselective direct substitution of one SNARE for another. Rather our data are consistent with the existence of multiple (perhaps parallel) trafficking pathways where Sed5p-containing SNARE complexes play overlapping and/or distinct functional roles.     

Differential specificities and simultaneous occupancy of human MutSalpha nucleotide binding sites

We have examined the permissible nucleotide occupancy states of human MutSalpha. The MSH2.MSH6 heterodimer binds 1 mol of ADP and 1 mol of adenosine 5'-O-(thiotriphosphate) (ATPgammaS), with a K(d) for each nucleotide of about 1 microm. Anisotropy measurements using BODIPY TR and BODIPY FL fluorescent derivatives of ADP and 5'-adenylyl-beta,gamma-imidodiphosphate (AMPPNP) also indicate an interaction stoichiometry of 1 mol of ADP and 1 mol of triphosphate analogue per MutSalpha heterodimer. Di- and triphosphate sites can be simultaneously occupied as judged by sequential filling of the two binding site classes with differentially radiolabeled ADP and ATPgammaS and by fluorescence resonance energy transfer between BODIPY TR- and BODIPY FL-labeled ADP and AMPPNP. ATP hydrolysis by MutSalpha is accompanied by a pre-steady-state burst of ADP formation, and analysis of MutSalpha-bound nucleotide during the first turnover has demonstrated the presence of both ADP and ATP. Simultaneous presence of ADP and a nonhydrolyzable ATP analogue modulates MutSalpha.heteroduplex interaction in a manner that is distinct from that observed in the presence of ADP or nonhydrolyzable triphosphate alone, and it is unlikely that this effect is due to the presence of a mixed population of binary complexes between MutSalpha and ADP or a triphosphate analogue. These findings imply that MutSalpha has two nucleotide binding sites with differential specificities for ADP and ATP and suggest that the ADP.MutSalpha.ATP ternary complex has an important role in mismatch repair.     

Combinatorial interactions between positive and negative cis-acting elements control spatial patterns of 4CL-1 expression in transgenic tobacco

The phenylpropanoid enzyme 4-coumarate:coenzyme A ligase (4CL) plays a key role in linking general phenylpropanoid metabolism to end-product specific biosynthetic pathways. During vascular system and floral organ differentiation, the parsley 4CL-1 gene is expressed in a restricted set of tissues and cell types where 4CL activity is required to supply precursors for the synthesis of diverse phenylpropanoid-derived products such as lignin and flavonoids. In order to localize cis -acting elements which specify complex patterns of 4CL-1 expression, we analyzed the expression of internally deleted promoter fragment— GUS fusions in tobacco plants and parsley protoplasts. Elements located between −244 and −78 were required for most aspects of developmentally regulated expression. Within this region, three separate promoter domains containing partially redundant cis -elements directed vascular-specific expression when combined with a TATA-proximal domain. A negative cis -acting element which represses phloem expression was revealed in one of the domains and appears to be responsible for restricting vascular expression to the xylem. Distinct but overlapping promoter domain combinations were required for expression in floral organs, suggesting that different combinations of cis -acting elements may direct expression in different organs. Gel retardation assays were used to demonstrate the formation of DNA-protein complexes between factors present in nuclear extracts of parsley tissue culture cells and various tobacco organs and a 4CL-1 promoter fragment. Competition experiments showed that complex formation required the presence of a 42 bp promoter domain shown to be critical for 4CL-1 expression in vascular and floral tissues. The results are discussed in light of the coordinate expression of 4CL and other phenylpropanoid genes.     

Fugu and human sequence comparison identifies novel human genes and conserved non-coding sequences

The compact genome of the pufferfish, Fugu rubripes, has been proposed as a 'reference' genome to aid in annotating and analysing the human genome. We have annotated and compared 85 kb of Fugu sequence containing 17 genes with its homologous loci in the human draft genome and identified three 'novel' human genes that were missed or incompletely predicted by the previous gene prediction methods. Two of the novel genes contain zinc finger domains and are designated ZNF366 and ZNF367. They map to human chromosomes 5q13.2 and 9q22.32, respectively. The third novel gene, designated C9orf21, maps to chromosome 9q22.32. This gene is unique to vertebrates, and the protein encoded by it does not contain any known domains. We could not find human homologs for two Fugu genes, a novel chemokine gene and a kinase gene. These genes are either specific to teleosts or lost in the human lineage. The Fugu-human comparison identified several conserved non-coding sequences in the promoter and intronic regions. These sequences, conserved during 450 million years of vertebrate evolution, are likely to be involved in gene regulation. The 85 kb Fugu locus is dispersed over four human loci, occupying about 1.5 Mb. Contiguity is conserved in the human genome between six out of 16 Fugu gene pairs. These contiguous chromosomal segments should share a common evolutionary history dating back to the common ancestor of mammals and teleosts. We propose contiguity as strong evidence to identify orthologous genes in distant organisms. This study confirms the utility of the Fugu as a supplementary tool to uncover and confirm novel genes and putative gene regulatory regions in the human genome.     

Assessing the biological significance of primary structure consensus patterns using sequence databanks. I. Heat-shock and glucocorticoid control elements in eukaryotic promoters

We describe FORTRAN 77 software allowing for convenient searchingof any segmented and ambiguous pattern in the currently availableprotein or nucleotide sequence databanks. For proteins, thissoftware can be instrumental in defining conserved functionaldomains among non–homologous overall primary structures.For nucleic acids, it is used in detecting complex and/or lowconsensus structural or regulatory patterns. As first applicationswe have studied the distribution of short consensus sequencesbelieved to characterize heat–shock and glucocorticoidregulated promoters. This analysis allowed an evaluation ofthe specificity, probable role and thus biological significanceof various regions of these consensus. In addition, the expressionof several known genes are predicted to be heat–shockor glucocorticoid sensitive. Received on March 25, 1985; accepted on March 28, 1985     

DNA sequence patterns in human,mouse, and rabbit immunoglobulin kappa-genes

Summary DNA sequences of the human, mouse, and rabbit immunoglobulin kappa-gene (J-C regions) are compared with respect to various DNA patterns, including dyad symmetry pairings, runs of nucleotides, repeat clusters, and repeats that occur with unusually high frequency. The significant dyad symmetry pairings within each of the sequences emphasize the two control-enhancer elements of the J₅-C intron. Dyad symmetry pairs between the J-C region and a number of kappa variable (V)-gene domains suggest differences in the affinities between the V and J segments. It is the consensus heptamer rather than the consensus nonamer that embodies the longest V-J dyad symmetry combinations. In the rabbit there are long runs and repeat clusters of the sequences that identify regions of high duplication; these regions are absent in the human and mouse sequences. High-frequency oligonucleotides feature the consensus nonamer 5 to the J segments, especially in the mouse sequence.     

Motion-based analysis of spatial patterns by the human visual system

BACKGROUND: It is known that the visibility of patterns presented through stationary multiple slits is significantly improved by pattern movements. This study investigated whether this spatiotemporal pattern interpolation is supported by motion mechanisms, as opposed to the general belief that the human visual cortex initially analyses spatial patterns independent of their movements. RESULTS: Psychophysical experiments showed that multislit viewing could not be ascribed to such motion-irrelevant factors as retinal painting by tracking eye movements or an increase in the number of views by pattern movements. Pattern perception was more strongly impaired by the masking noise moving in the same direction than by the noise moving in the opposite direction, which indicates the direction selectivity of the pattern interpolation mechanism. A direction-selective impairment of pattern perception by motion adaptation also indicates the direction selectivity of the interpolation mechanism. Finally, the map of effective spatial frequencies, estimated by a reverse-correlation technique, indicates observers' perception of higher spatial frequencies, the recovery of which is theoretically impossible without the aid of motion information. CONCLUSIONS: These results provide clear evidence against the notion of separate analysis of pattern and motion. The visual system uses motion mechanisms to integrate spatial pattern information along the trajectory of pattern movement in order to obtain clear perception of moving patterns. The pattern integration mechanism is likely to be direction-selective filtering by V1 simple cells, but the integration of the local pattern information into a global figure should be guided by a higher-order motion mechanism such as MT pattern cells.