The relationship of potential G-quadruplex sequences in cis-upstream regions of the human genome to SP1-binding elements

We have carried out a survey of potential quadruplex structure sequences (PQSS), which occur in the immediate upstream region (500 bp) of human genes. By examining the number and distribution of these we have established that there is a clear link between them and the occurrence of the SP1-binding element ‘GGGCGG’, such that a large number of upstream PQSS incorporate the SP1-binding element.     

Systematic Analysis of Compositional Order of Proteins Reveals New Characteristics of Biological Functions and a Universal Correlate of Macroevolution

We present a novel analysis of compositional order (CO) based on the occurrence of Frequent amino-acid Triplets (FTs) that appear much more than random in protein sequences. The method captures all types of proteomic compositional order including single amino-acid runs, tandem repeats, periodic structure of motifs and otherwise low complexity amino-acid regions. We introduce new order measures, distinguishing between ‘regularity’, ‘periodicity’ and ‘vocabulary’, to quantify these phenomena and to facilitate the identification of evolutionary effects. Detailed analysis of representative species across the tree-of-life demonstrates that CO proteins exhibit numerous functional enrichments, including a wide repertoire of particular patterns of dependencies on regularity and periodicity. Comparison between human and mouse proteomes further reveals the interplay of CO with evolutionary trends, such as faster substitution rate in mouse leading to decrease of periodicity, while innovation along the human lineage leads to larger regularity. Large-scale analysis of 94 proteomes leads to systematic ordering of all major taxonomic groups according to FT-vocabulary size. This is measured by the count of Different Frequent Triplets (DFT) in proteomes. The latter provides a clear hierarchical delineation of vertebrates, invertebrates, plants, fungi and prokaryotes, with thermophiles showing the lowest level of FT-vocabulary. Among eukaryotes, this ordering correlates with phylogenetic proximity. Interestingly, in all kingdoms CO accumulation in the proteome has universal characteristics. We suggest that CO is a genomic-information correlate of both macroevolution and various protein functions. The results indicate a mechanism of genomic ‘innovation’ at the peptide level, involved in protein elongation, shaped in a universal manner by mutational and selective forces.     

Can Natural Proteins Designed with ‘Inverted’ Peptide Sequences Adopt Native-Like Protein Folds?

We have carried out a systematic computational analysis on a representative dataset of proteins of known three-dimensional structure, in order to evaluate whether it would possible to ‘swap’ certain short peptide sequences in naturally occurring proteins with their corresponding ‘inverted’ peptides and generate ‘artificial’ proteins that are predicted to retain native-like protein fold. The analysis of 3,967 representative proteins from the Protein Data Bank revealed 102,677 unique identical inverted peptide sequence pairs that vary in sequence length between 5–12 and 18 amino acid residues. Our analysis illustrates with examples that such ‘artificial’ proteins may be generated by identifying peptides with ‘similar structural environment’ and by using comparative protein modeling and validation studies. Our analysis suggests that natural proteins may be tolerant to accommodating such peptides.     

Duplex Scorpion primers in SNP analysis and FRET applications

Scorpions are fluorogenic PCR primers with a probe element attached at the 5′-end via a PCR stopper. They are used in real-time amplicon-specific detection of PCR products in homogeneous solution. Two different formats are possible, the ‘stem–loop’ format and the ‘duplex’ format. In both cases the probing mechanism is intramolecular. We have shown that duplex Scorpions are efficient probes in real-time PCR. They give a greater fluorescent signal than stem–loop Scorpions due to the vastly increased separation between fluorophore and quencher in the active form. We have demonstrated their use in allelic discrimination at the W1282X locus of the ABCC7 gene and shown that they can be used in assays where fluorescence resonance energy transfer is required.     

Differential Beta-Band Event-Related Desynchronization during Categorical Action Sequence Planning

A primate study reported the existence of neurons from the dorso-lateral prefrontal cortex which fired prior to executing categorical action sequences. The authors suggested these activities may represent abstract level information. Here, we aimed to find the neurophysiological representation of planning categorical action sequences at the population level in healthy humans. Previous human studies have shown beta-band event-related desynchronization (ERD) during action planning in humans. Some of these studies showed different levels of ERD according to different types of action preparation. Especially, the literature suggests that variations in cognitive factors rather than physical factors (force, direction, etc) modulate the level of beta-ERD. We hypothesized that the level of beta-band power will differ according to planning of different categorical sequences. We measured magnetoencephalography (MEG) from 22 subjects performing 11 four-sequence actions - each consisting of one or two of three simple actions - in 3 categories; ‘Paired (ooxx)’, ‘Alternative (oxox)’ and ‘Repetitive (oooo)’ (‘o’ and ‘x’ each denoting one of three simple actions). Time-frequency representations were calculated for each category during the planning period, and the corresponding beta-power time-courses were compared. We found beta-ERD during the planning period for all subjects, mostly in the contralateral fronto-parietal areas shortly after visual cue onset. Power increase (transient rebound) followed ERD in 20 out of 22 subjects. Amplitudes differed among categories in 20 subjects for both ERD and transient rebound. In 18 out of 20 subjects ‘Repetitive’ category showed the largest ERD and rebound. The current result suggests that beta-ERD in the contralateral frontal/motor/parietal areas during planning is differentiated by the category of action sequences.     

Rational design of DNA sequences for nanotechnology, microarrays and molecular computers using Eulerian graphs

Nucleic acids are molecules of choice for both established and emerging nanoscale technologies. These technologies benefit from large functional densities of ‘DNA processing elements’ that can be readily manufactured. To achieve the desired functionality, polynucleotide sequences are currently designed by a process that involves tedious and laborious filtering of potential candidates against a series of requirements and parameters. Here, we present a complete novel methodology for the rapid rational design of large sets of DNA sequences. This method allows for the direct implementation of very complex and detailed requirements for the generated sequences, thus avoiding ‘brute force’ filtering. At the same time, these sequences have narrow distributions of melting temperatures. The molecular part of the design process can be done without computer assistance, using an efficient ‘human engineering’ approach by drawing a single blueprint graph that represents all generated sequences. Moreover, the method eliminates the necessity for extensive thermodynamic calculations. Melting temperature can be calculated only once (or not at all). In addition, the isostability of the sequences is independent of the selection of a particular set of thermodynamic parameters. Applications are presented for DNA sequence designs for microarrays, universal microarray zip sequences and electron transfer experiments.     

Polyethylene glycol binding alters human telomere G-quadruplex structure by conformational selection

Polyethylene glycols (PEGs) are widely used to perturb the conformations of nucleic acids, including G-quadruplexes. The mechanism by which PEG alters G-quadruplex conformation is poorly understood. We describe here studies designed to determine how PEG and other co-solutes affect the conformation of the human telomeric quadruplex. Osmotic stress studies using acetonitrile and ethylene glycol show that conversion of the ‘hybrid’ conformation to an all-parallel ‘propeller’ conformation is accompanied by the release of about 17 water molecules per quadruplex and is energetically unfavorable in pure aqueous solutions. Sedimentation velocity experiments show that the propeller form is hydrodynamically larger than hybrid forms, ruling out a crowding mechanism for the conversion by PEG. PEGs do not alter water activity sufficiently to perturb quadruplex hydration by osmotic stress. PEG titration experiments are most consistent with a conformational selection mechanism in which PEG binds more strongly to the propeller conformation, and binding is coupled to the conformational transition between forms. Molecular dynamics simulations show that PEG binding to the propeller form is sterically feasible and energetically favorable. We conclude that PEG does not act by crowding and is a poor mimic of the intranuclear environment, keeping open the question of the physiologically relevant quadruplex conformation.     

Kissing complex-mediated dimerisation of HIV-1 RNA: coupling extended duplex formation to ribozyme cleavage

Initiation of retroviral genomic RNA dimerisation is mediated by the mutual interaction of the dimerisation initiation site (DIS) stem–loops near to the 5′ end of the RNA. This process is thought to involve formation of a transient ‘kissing’ complex over the self-complementary loop bases, which then refolds into a more stable extended interaction. We have developed a novel experimental system that allows us to clearly detect the extended duplex in vitro. Ribozyme sequences were incorporated into or adjacent to the type 1 human immunodeficiency virus DIS stem, leading to the formation of a functional ribozyme only in the extended duplex conformer. Here we show that extended duplex formation results in ribozyme cleavage, thus demonstrating the double-stranded nature of the extended complex and confirming that refolding occurs via melting of the DIS stems. Loop complementarity is essential for extended duplex formation but no sequence requirements for the loops were observed. Efficiency of extended duplex formation is dependent on the strength of the loop–loop interaction, temperature, the magnesium concentration and is strongly accelerated by the viral nucleocapsid protein NCp7. Our ribozyme-coupled approach should be applicable to the analyses of other refolding processes involving RNA loop–loop interactions.     

Realistic artificial DNA sequences as negative controls for computational genomics

A common practice in computational genomic analysis is to use a set of ‘background’ sequences as negative controls for evaluating the false-positive rates of prediction tools, such as gene identification programs and algorithms for detection of cis-regulatory elements. Such ‘background’ sequences are generally taken from regions of the genome presumed to be intergenic, or generated synthetically by ‘shuffling’ real sequences. This last method can lead to underestimation of false-positive rates. We developed a new method for generating artificial sequences that are modeled after real intergenic sequences in terms of composition, complexity and interspersed repeat content. These artificial sequences can serve as an inexhaustible source of high-quality negative controls. We used artificial sequences to evaluate the false-positive rates of a set of programs for detecting interspersed repeats, ab initio prediction of coding genes, transcribed regions and non-coding genes. We found that RepeatMasker is more accurate than PClouds, Augustus has the lowest false-positive rate of the coding gene prediction programs tested, and Infernal has a low false-positive rate for non-coding gene detection. A web service, source code and the models for human and many other species are freely available at http://repeatmasker.org/garlic/.     

Conceptual learning by miniature brains

Concepts act as a cornerstone of human cognition. Humans and non-human primates learn conceptual relationships such as ‘same’, ‘different’, ‘larger than’, ‘better than’, among others. In all cases, the relationships have to be encoded by the brain independently of the physical nature of objects linked by the relation. Consequently, concepts are associated with high levels of cognitive sophistication and are not expected in an insect brain. Yet, various works have shown that the miniature brain of honeybees rapidly learns conceptual relationships involving visual stimuli. Concepts such as ‘same’, ‘different’, ‘above/below of’ or ‘left/right are well mastered by bees. We review here evidence about concept learning in honeybees and discuss both its potential adaptive advantage and its possible neural substrates. The results reviewed here challenge the traditional view attributing supremacy to larger brains when it comes to the elaboration of concepts and have wide implications for understanding how brains can form conceptual relations.     

Identifying the Computational Requirements of an Integrated Top-Down-Bottom-Up Model for Overt Visual Attention within an Active Vision System

Computational visual attention systems have been constructed in order for robots and other devices to detect and locate regions of interest in their visual world. Such systems often attempt to take account of what is known of the human visual system and employ concepts, such as ‘active vision’, to gain various perceived advantages. However, despite the potential for gaining insights from such experiments, the computational requirements for visual attention processing are often not clearly presented from a biological perspective. This was the primary objective of this study, attained through two specific phases of investigation: 1) conceptual modeling of a top-down-bottom-up framework through critical analysis of the psychophysical and neurophysiological literature, 2) implementation and validation of the model into robotic hardware (as a representative of an active vision system). Seven computational requirements were identified: 1) transformation of retinotopic to egocentric mappings, 2) spatial memory for the purposes of medium-term inhibition of return, 3) synchronization of ‘where’ and ‘what’ information from the two visual streams, 4) convergence of top-down and bottom-up information to a centralized point of information processing, 5) a threshold function to elicit saccade action, 6) a function to represent task relevance as a ratio of excitation and inhibition, and 7) derivation of excitation and inhibition values from object-associated feature classes. The model provides further insight into the nature of data representation and transfer between brain regions associated with the vertebrate ‘active’ visual attention system. In particular, the model lends strong support to the functional role of the lateral intraparietal region of the brain as a primary area of information consolidation that directs putative action through the use of a ‘priority map’.     

Fundamental structural units of the Escherichia coli nucleoid revealed by atomic force microscopy

A small container of several to a few hundred µm³ (i.e. bacterial cells and eukaryotic nuclei) contains extremely long genomic DNA (i.e. mm and m long, respectively) in a highly organized fashion. To understand how such genomic architecture could be achieved, Escherichia coli nucleoids were subjected to structural analyses under atomic force microscopy, and found to change their structure dynamically during cell growth, i.e. the nucleoid structure in the stationary phase was more tightly compacted than in the log phase. However, in both log and stationary phases, a fundamental fibrous structure with a diameter of ~80 nm was found. In addition to this ‘80 nm fiber’, a thinner ‘40 nm fiber’ and a higher order ‘loop’ structure were identified in the log phase nucleoid. In the later growth phases, the nucleoid turned into a ‘coral reef structure’ that also possessed the 80 nm fiber units, and, finally, into a ‘tightly compacted nucleoid’ that was stable in a mild lysis buffer. Mutant analysis demonstrated that these tight compactions of the nucleoid required a protein, Dps. From these results and previously available information, we propose a structural model of the E.coli nucleoid.     

Genome sequencing and analysis of two early-flowering cherry (Cerasus × kanzakura) varieties, ‘Kawazu-zakura’ and ‘Atami-zakura’

To gain genetic insights into the early-flowering phenotype of ornamental cherry, also known as sakura, we determined the genome sequences of two early-flowering cherry (Cerasus × kanzakura) varieties, ‘Kawazu-zakura’ and ‘Atami-zakura’. Because the two varieties are interspecific hybrids, likely derived from crosses between Cerasus campanulata (early-flowering species) and Cerasus speciosa, we employed the haplotype-resolved sequence assembly strategy. Genome sequence reads obtained from each variety by single-molecule real-time sequencing (SMRT) were split into two subsets, based on the genome sequence information of the two probable ancestors, and assembled to obtain haplotype-phased genome sequences. The resultant genome assembly of ‘Kawazu-zakura’ spanned 519.8 Mb with 1,544 contigs and an N50 value of 1,220.5 kb, while that of ‘Atami-zakura’ totalled 509.6 Mb with 2,180 contigs and an N50 value of 709.1 kb. A total of 72,702 and 69,528 potential protein-coding genes were predicted in the genome assemblies of ‘Kawazu-zakura’ and ‘Atami-zakura’, respectively. Gene clustering analysis identified 2,634 clusters uniquely presented in the C. campanulata haplotype sequences, which might contribute to its early-flowering phenotype. Genome sequences determined in this study provide fundamental information for elucidating the molecular and genetic mechanisms underlying the early-flowering phenotype of ornamental cherry tree varieties and their relatives.     

Breaking an Epigenetic Chromatin Switch: Curious Features of Hysteresis in Saccharomyces cerevisiae Telomeric Silencing

In addition to gene network switches, local epigenetic modifications to DNA and histones play an important role in all-or-none cellular decision-making. Here, we study the dynamical design of a well-characterized epigenetic chromatin switch: the yeast SIR system, in order to understand the origin of the stability of epigenetic states. We study hysteresis in this system by perturbing it with a histone deacetylase inhibitor. We find that SIR silencing has many characteristics of a non-linear bistable system, as observed in conventional genetic switches, which are based on activities of a few promoters affecting each other through the abundance of their gene products. Quite remarkably, our experiments in yeast telomeric silencing show a very distinctive pattern when it comes to the transition from bistability to monostability. In particular, the loss of the stable silenced state, upon increasing the inhibitor concentration, does not seem to show the expected saddle node behavior, instead looking like a supercritical pitchfork bifurcation. In other words, the ‘off’ state merges with the ‘on’ state at a threshold concentration leading to a single state, as opposed to the two states remaining distinct up to the threshold and exhibiting a discontinuous jump from the ‘off’ to the ‘on’ state. We argue that this is an inevitable consequence of silenced and active regions coexisting with dynamic domain boundaries. The experimental observations in our study therefore have broad implications for the understanding of chromatin silencing in yeast and beyond.     

An in vitro selection scheme for oligonucleotide probes to discriminate between closely related DNA sequences

Using an in vitro selection, we have obtained oligonucleotide probes with high discriminatory power against multiple, similar nucleic acid sequences, which is often required in diagnostic applications for simultaneous testing of such sequences. We have tested this approach, referred to as iterative hybridizations, by selecting probes against six 22-nt-long sequence variants representing human papillomavirus, (HPV). We have obtained probes that efficiently discriminate between HPV types that differ by 3–7nt. The probes were found effective to recognize HPV sequences of the type 6, 11, 16, 18 and a pair of type 31 and 33, either when immobilized on a solid support or in a reverse configuration, as well to discriminate HPV types from the clinical samples. This methodology can be extended to generate diagnostic kits that rely on nucleic acid hybridization between closely related sequences. In this approach, instead of adjusting hybridization conditions to the intended set of probe–target pairs, we ‘adjust’, through in vitro selection, the probes to the conditions we have chosen. Importantly, these conditions have to be ‘relaxed’, allowing the formation of a variety of not fully complementary complexes from which those that efficiently recognize and discriminate intended from non-intended targets can be readily selected.