Iodine

by Carlos Gershenson

The performance of public transportation systems affects a large part of the population. Current theory assumes that passengers are served optimally when vehicles arrive at stations with regular intervals. In this paper, it is shown that self-organization can improve the performance of public transportation systems beyond the theoretical optimum by responding adaptively to local conditions. This is possible because of a “slower-is-faster” effect, where passengers wait more time at stations but total travel times are reduced. The proposed self-organizing method uses “antipheromones” to regulate headways, which are inspired by the stigmergy (communication via environment) of some ant colonies.

by Eugeni Belda, Laia Pedrola, Juli Peretó, Juan F. Martínez-Blanch, Arnau Montagud, Emilio Navarro, Javier Urchueguía, Daniel Ramón, Andrés Moya, Manuel Porcar

Background

Insects are associated with microorganisms that contribute to the digestion and processing of nutrients. The European Corn Borer (ECB) is a moth present world-wide, causing severe economical damage as a pest on corn and other crops. In the present work, we give a detailed view of the complexity of the microorganisms forming the ECB midgut microbiota with the objective of comparing the biodiversity of the midgut-associated microbiota and explore their potential as a source of genes and enzymes with biotechnological applications.

Methodological/Principal Findings

A high-throughput sequencing approach has been used to identify bacterial species, genes and metabolic pathways, particularly those involved in plant-matter degradation, in two different ECB populations (field-collected vs. lab-reared population with artificial diet). Analysis of the resulting sequences revealed the massive presence of Staphylococcus warneri and Weissella paramesenteroides in the lab-reared sample. This enabled us to reconstruct both genomes almost completely. Despite the apparently low diversity, 208 different genera were detected in the sample, although most of them at very low frequency. By contrast, the natural population exhibited an even higher taxonomic diversity along with a wider array of cellulolytic enzyme families. However, in spite of the differences in relative abundance of major taxonomic groups, not only did both metagenomes share a similar functional profile but also a similar distribution of non-redundant genes in different functional categories.

Conclusions/Significance

Our results reveal a highly diverse pool of bacterial species in both O. nubilalis populations, with major differences: The lab-reared sample is rich in gram-positive species (two of which have almost fully sequenced genomes) while the field sample harbors mainly gram-negative species and has a larger set of cellulolytic enzymes. We have found a clear relationship between the diet and the midgut microbiota, which reveals the selection pressure of food on the community of intestinal bacteria.

by Armando J. Pinho, Paulo J. S. G. Ferreira, António J. R. Neves, Carlos A. C. Bastos

A finite-context (Markov) model of order yields the probability distribution of the next symbol in a sequence of symbols, given the recent past up to depth . Markov modeling has long been applied to DNA sequences, for example to find gene-coding regions. With the first studies came the discovery that DNA sequences are non-stationary: distinct regions require distinct model orders. Since then, Markov and hidden Markov models have been extensively used to describe the gene structure of prokaryotes and eukaryotes. However, to our knowledge, a comprehensive study about the potential of Markov models to describe complete genomes is still lacking. We address this gap in this paper. Our approach relies on (i) multiple competing Markov models of different orders (ii) careful programming techniques that allow orders as large as sixteen (iii) adequate inverted repeat handling (iv) probability estimates suited to the wide range of context depths used. To measure how well a model fits the data at a particular position in the sequence we use the negative logarithm of the probability estimate at that position. The measure yields information profiles of the sequence, which are of independent interest. The average over the entire sequence, which amounts to the average number of bits per base needed to describe the sequence, is used as a global performance measure. Our main conclusion is that, from the probabilistic or information theoretic point of view and according to this performance measure, multiple competing Markov models explain entire genomes almost as well or even better than state-of-the-art DNA compression methods, such as XM, which rely on very different statistical models. This is surprising, because Markov models are local (short-range), contrasting with the statistical models underlying other methods, where the extensive data repetitions in DNA sequences is explored, and therefore have a non-local character.

by Phillip Senum, Marc Riedel

This paper presents a collection of computational modules implemented with chemical reactions: an inverter, an incrementer, a decrementer, a copier, a comparator, a multiplier, an exponentiator, a raise-to-a-power operation, and a logarithm in base two. Unlike previous schemes for chemical computation, this method produces designs that are dependent only on coarse rate categories for the reactions (“fast” vs. “slow”). Given such categories, the computation is exact and independent of the specific reaction rates. The designs are validated through stochastic simulations of the chemical kinetics.

by Georges Natsoulis, John M. Bell, Hua Xu, Jason D. Buenrostro, Heather Ordonez, Susan Grimes, Daniel Newburger, Michael Jensen, Jacob M. Zahn, Nancy Zhang, Hanlee P. Ji

We have developed an integrated strategy for targeted resequencing and analysis of gene subsets from the human exome for variants. Our capture technology is geared towards resequencing gene subsets substantially larger than can be done efficiently with simplex or multiplex PCR but smaller in scale than exome sequencing. We describe all the steps from the initial capture assay to single nucleotide variant (SNV) discovery. The capture methodology uses in-solution 80-mer oligonucleotides. To provide optimal flexibility in choosing human gene targets, we designed an in silico set of oligonucleotides, the Human OligoExome, that covers the gene exons annotated by the Consensus Coding Sequencing Project (CCDS). This resource is openly available as an Internet accessible database where one can download capture oligonucleotides sequences for any CCDS gene and design custom capture assays. Using this resource, we demonstrated the flexibility of this assay by custom designing capture assays ranging from 10 to over 100 gene targets with total capture sizes from over 100 Kilobases to nearly one Megabase. We established a method to reduce capture variability and incorporated indexing schemes to increase sample throughput. Our approach has multiple applications that include but are not limited to population targeted resequencing studies of specific gene subsets, validation of variants discovered in whole genome sequencing surveys and possible diagnostic analysis of disease gene subsets. We also present a cost analysis demonstrating its cost-effectiveness for large population studies.

by Hee-Jung Kim, Hye Joon Joo, Yung Hee Kim, Soyeon Ahn, Jun Chang, Kyu-Baek Hwang, Dong-Hee Lee, Kong-Joo Lee

The molecular basis of heat shock response (HSR), a cellular defense mechanism against various stresses, is not well understood. In this, the first comprehensive analysis of gene expression changes in response to heat shock and MG132 (a proteasome inhibitor), both of which are known to induce heat shock proteins (Hsps), we compared the responses of normal mouse fibrosarcoma cell line, RIF- 1, and its thermotolerant variant cell line, TR-RIF-1 (TR), to the two stresses. The cellular responses we examined included Hsp expressions, cell viability, total protein synthesis patterns, and accumulation of poly-ubiquitinated proteins. We also compared the mRNA expression profiles and kinetics, in the two cell lines exposed to the two stresses, using microarray analysis. In contrast to RIF-1 cells, TR cells resist heat shock caused changes in cell viability and whole-cell protein synthesis. The patterns of total cellular protein synthesis and accumulation of poly-ubiquitinated proteins in the two cell lines were distinct, depending on the stress and the cell line. Microarray analysis revealed that the gene expression pattern of TR cells was faster and more transient than that of RIF-1 cells, in response to heat shock, while both RIF-1 and TR cells showed similar kinetics of mRNA expression in response to MG132. We also found that 2,208 genes were up-regulated more than 2 fold and could sort them into three groups: 1) genes regulated by both heat shock and MG132, (e.g. chaperones); 2) those regulated only by heat shock (e.g. DNA binding proteins including histones); and 3) those regulated only by MG132 (e.g. innate immunity and defense related molecules). This study shows that heat shock and MG132 share some aspects of HSR signaling pathway, at the same time, inducing distinct stress response signaling pathways, triggered by distinct abnormal proteins.

by Benjamin Lehne, Cathryn M. Lewis, Thomas Schlitt

Interpreting Genome-Wide Association Studies (GWAS) at a gene level is an important step towards understanding the molecular processes that lead to disease. In order to incorporate prior biological knowledge such as pathways and protein interactions in the analysis of GWAS data it is necessary to derive one measure of association for each gene. We compare three different methods to obtain gene-wide test statistics from Single Nucleotide Polymorphism (SNP) based association data: choosing the test statistic from the most significant SNP; the mean test statistics of all SNPs; and the mean of the top quartile of all test statistics. We demonstrate that the gene-wide test statistics can be controlled for the number of SNPs within each gene and show that all three methods perform considerably better than expected by chance at identifying genes with confirmed associations. By applying each method to GWAS data for Crohn's Disease and Type 1 Diabetes we identified new potential disease genes.

by Katja Kornysheva, Ricarda I. Schubotz

Integrating auditory and motor information often requires precise timing as in speech and music. In humans, the position of the ventral premotor cortex (PMv) in the dorsal auditory stream renders this area a node for auditory-motor integration. Yet, it remains unknown whether the PMv is critical for auditory-motor timing and which activity increases help to preserve task performance following its disruption. 16 healthy volunteers participated in two sessions with fMRI measured at baseline and following rTMS (rTMS) of either the left PMv or a control region. Subjects synchronized left or right finger tapping to sub-second beat rates of auditory rhythms in the experimental task, and produced self-paced tapping during spectrally matched auditory stimuli in the control task. Left PMv rTMS impaired auditory-motor synchronization accuracy in the first sub-block following stimulation (p<0.01, Bonferroni corrected), but spared motor timing and attention to task. Task-related activity increased in the homologue right PMv, but did not predict the behavioral effect of rTMS. In contrast, anterior midline cerebellum revealed most pronounced activity increase in less impaired subjects. The present findings suggest a critical role of the left PMv in feed-forward computations enabling accurate auditory-motor timing, which can be compensated by activity modulations in the cerebellum, but not in the homologue region contralateral to stimulation.

by Juan Rodriguez-Paris, Marta L. Tamayo, Nancy Gelvez, Iris Schrijver

Mutations in the GJB2 gene, which encodes connexin 26, are a frequent cause of congenital non-syndromic sensorineural hearing loss. Two large deletions, del(GJB6-D13S1830) and del(GJB6-D13S1854), which truncate GJB6 (connexin 30), cause hearing loss in individuals homozygous, or compound heterozygous for these deletions or one such deletion and a mutation in GJB2. Recently, we have demonstrated that the del(GJB6-D13S1830) deletion contributes to hearing loss due to an allele-specific lack of GJB2 mRNA expression and not as a result of digenic inheritance, as was postulated earlier. In the current study we investigated the smaller del(GJB6-D13S1854) deletion, which disrupts the expression of GJB2 at the transcriptional level in a manner similar to the more common del(GJB6-D13S1830) deletion. Interestingly, in the presence of this deletion, GJB2 expression remains minimally but reproducibly present. The relative allele-specific expression of GJB2 was assessed by reverse-transcriptase PCR and restriction digestions in three probands who were compound heterozygous for a GJB2 mutation and del(GJB6-D13S1854). Each individual carried a different sequence variant in GJB2. All three individuals expressed the mutated GJB2 allele in trans with del(GJB6-D13S1854), but expression of the GJB2 allele in cis with the deletion was almost absent. Our study clearly corroborates the hypothesis that the del(GJB6-D13S1854), similar to the larger and more common del(GJB6-D13S1830), removes (a) putative cis-regulatory element(s) upstream of GJB6 and narrows down the region of location.

by Seif Eldawlatly, Karim G. Oweiss

Correlation among neocortical neurons is thought to play an indispensable role in mediating sensory processing of external stimuli. The role of temporal precision in this correlation has been hypothesized to enhance information flow along sensory pathways. Its role in mediating the integration of information at the output of these pathways, however, remains poorly understood. Here, we examined spike timing correlation between simultaneously recorded layer V neurons within and across columns of the primary somatosensory cortex of anesthetized rats during unilateral whisker stimulation. We used Bayesian statistics and information theory to quantify the causal influence between the recorded cells with millisecond precision. For each stimulated whisker, we inferred stable, whisker-specific, dynamic Bayesian networks over many repeated trials, with network similarity of 83.3±6% within whisker, compared to only 50.3±18% across whiskers. These networks further provided information about whisker identity that was approximately 6 times higher than what was provided by the latency to first spike and 13 times higher than what was provided by the spike count of individual neurons examined separately. Furthermore, prediction of individual neurons' precise firing conditioned on knowledge of putative pre-synaptic cell firing was 3 times higher than predictions conditioned on stimulus onset alone. Taken together, these results suggest the presence of a temporally precise network coding mechanism that integrates information across neighboring columns within layer V about vibrissa position and whisking kinetics to mediate whisker movement by motor areas innervated by layer V.