Clusters of time series data may change location and memberships over time; in gene expression data, this occurs as groups of genes or samples respond differently to stimuli or experimental conditions at different times. In order to uncover this underlying temporal structure, we consider dynamic clusters with time-dependent parameters which split and merge over time, enabling cluster memberships to change. These interesting time-dependent structures are useful in understanding the development of organisms or complex organs, and could not be identified using traditional clustering methods. In cell cycle data, these time-dependent structure may provide links between genes and stages of the cell cycle, whilst in developmental data sets they may highlight key developmental transitions.

The manner in which different distributions of synaptic weights onto cortical neurons shape their spiking activity remains open. To characterize a homogeneous neuronal population, we use the master equation for generalized leaky integrate-and-fire neurons with shot-noise synapses. We develop fast semi-analytic numerical methods to solve this equation for either current or conductance synapses, with and without synaptic depression. We show that its solutions match simulations of equivalent neuronal networks better than those of the Fokker-Planck equation and we compute bounds on the network response to non-instantaneous synapses. We apply these methods to study different synaptic weight distributions in feed-forward networks. We characterize the synaptic amplitude distributions using a set of measures, called tail weight numbers, designed to quantify the preponderance of very strong synapses. Even if synaptic amplitude distributions are equated for both the total current and average synaptic weight, distributions with sparse but strong synapses produce higher responses for small inputs, leading to a larger operating range. Furthermore, despite their small number, such synapses enable the network to respond faster and with more stability in the face of external fluctuations.

Background: Segmenting electron microscopy (EM) images of cellular and subcellular processes in the nervous system is a key step in many bioimaging pipelines involving classification and labeling of ultrastructures. However, fully automated techniques to segment images are often susceptible to noise and heterogeneity in EM images (e.g. different histological preparations, different organisms, different brain regions, etc.). Supervised techniques to address this problem are often helpful but require large sets of training data, which are often difficult to obtain in practice, especially across many conditions. Results: We propose a new, principled unsupervised algorithm to segment EM images using a two-step approach: edge detection via salient watersheds following by robust region merging. We performed experiments to gather EM neuroimages of two organisms (mouse and fruit fly) using different histological preparations and generated manually curated ground-truth segmentations. We compared our algorithm against several state-of- the-art unsupervised segmentation algorithms and found superior performance using two standard measures of under-and over-segmentation error. Conclusions: Our algorithm is general and may be applicable to other large-scale segmentation problems for bioimages.

Many animals, including insects, are known to use visual landmarks to orient in their environment. In Drosophila melanogaster, behavioural genetics studies have identified a higher brain structure called the central complex as being required for the fly’s innate responses to vertical visual features and its short- and long-term memory for visual patterns. But whether and how neurons of the fly central complex represent visual features are unknown. Here we use two-photon calcium imaging in head-fixed walking and flying flies to probe visuomotor responses of ring neurons—a class of central complex neurons that have been implicated in landmark-driven spatial memory in walking flies and memory for visual patterns in tethered flying flies. We show that dendrites of ring neurons are visually responsive and arranged retinotopically. Ring neuron receptive fields comprise both excitatory and inhibitory subfields, resembling those of simple cells in the mammalian primary visual cortex. Ring neurons show strong and, in some cases, direction-selective orientation tuning, with a notable preference for vertically oriented features similar to those that evoke innate responses in flies. Visual responses were diminished during flight, but, in contrast with the hypothesized role of the central complex in the control of locomotion, not modulated during walking. Taken together, these results indicate that ring neurons represent behaviourally relevant visual features in the fly’s environment, enabling downstream central complex circuits to produce appropriate motor commands. More broadly, this study opens the door to mechanistic investigations of circuit computations underlying visually guided action selection in the Drosophila central complex.

Abstract We recently described a new form of neural integration and firing in a subset of interneurons, in which evoking hundreds of action potentials over tens of seconds to minutes produces a sudden barrage of action potentials lasting about a minute beyond the inciting stimulation. During this persistent firing, action potentials are generated in the distal axon and propagate retrogradely to the soma. To distinguish this from other forms of persistent firing, we refer to it here as ’retroaxonal barrage firing’, or ’barrage firing’ for short. Its induction is blocked by chemical inhibitors of gap junctions and curiously, stimulation of one interneuron in some cases triggers barrage firing in a nearby, unstimulated interneuron. Beyond these clues, the mechanisms of barrage firing are unknown. Here we report new results related to these mechanisms. Induction of barrage firing was blocked by lowering extracellular calcium, as long as normal action potential threshold was maintained, and it was inhibited by blocking L-type voltage-gated calcium channels. Despite its calcium dependence, barrage firing was not prevented by inhibiting chemical synaptic transmission. Furthermore, loading the stimulated/recorded interneuron with BAPTA did not block barrage firing, suggesting that the required calcium entry occurs in other cells. Finally, barrage firing was normal in mice with deletion of the primary gene for neuronal gap junctions (connexin36), suggesting that non-neuronal gap junctions may be involved. Together, these findings suggest that barrage firing is probably triggered by a multicellular mechanism involving calcium signalling and gap junctions, but operating independently of chemical synaptic transmission.

SUMMARY: Sequence database searches are an essential part of molecular biology, providing information about the function and evolutionary history of proteins, RNA molecules and DNA sequence elements. We present a tool for DNA/DNA sequence comparison that is built on the HMMER framework, which applies probabilistic inference methods based on hidden Markov models to the problem of homology search. This tool, called nhmmer, enables improved detection of remote DNA homologs, and has been used in combination with Dfam and RepeatMasker to improve annotation of transposable elements in the human genome. AVAILABILITY: nhmmer is a part of the new HMMER3.1 release. Source code and documentation can be downloaded from http://hmmer.org. HMMER3.1 is freely licensed under the GNU GPLv3 and should be portable to any POSIX-compliant operating system, including Linux and Mac OS/X. CONTACT: wheelert@janelia.hhmi.org.

We demonstrate a meaningful prospective power analysis for an (admittedly idealized) illustrative connectome inference task. Modeling neurons as vertices and synapses as edges in a simple random graph model, we optimize the trade-off between the number of (putative) edges identified and the accuracy of the edge identification procedure. We conclude that explicit analysis of the quantity/quality trade-off is imperative for optimal neuroscientific experimental design. In particular, identifying edges faster/more cheaply, but with more error, can yield superior inferential performance.

A variant Hb zeta(2)beta(s)(2) that is formed from sickle hemoglobin (Hb S; alpha(2)beta(s)(2)) by exchanging adult alpha-globin with embryonic zeta-globin subunits shows promise as a therapeutic agent for sickle-cell disease (SCD). Hb zeta(2)beta(s)(2) inhibits the polymerization of deoxygenated Hb S in vitro and reverses characteristic features of SCD in vivo in mouse models of the disorder. When compared with either Hb S or with normal human adult Hb A (alpha(2)beta(2)), Hb zeta(2)beta(s)(2) exhibits atypical properties that include a high oxygen affinity, reduced cooperativity, a weak Bohr effect and blunted 2,3-diphosphoglycerate allostery. Here, the 1.95 angstrom resolution crystal structure of human Hb zeta(2)beta(s)(2) that was expressed in complex transgenic knockout mice and purified from their erythrocytes is presented. When fully liganded with carbon monoxide, Hb zeta(2)beta(s)(2) displays a central water cavity, a zeta 1-beta(s)2 (or zeta 2-beta(s)1) interface, intersubunit salt-bridge/hydrogen-bond interactions, C-terminal beta His146 salt-bridge interactions, and a beta-cleft, that are highly unusual for a relaxed hemoglobin structure and are more typical of a tense conformation. These quaternary tense-like features contrast with the tertiary relaxed-like conformations of the zeta 1-beta(s1) dimer and the CD and FG corners, as well as the overall structures of the heme cavities. This crystallographic study provides insights into the altered oxygen-transport properties of Hb zeta(2)beta(s)(2) and, moreover, decouples tertiary- and quaternary-structural events that are critical to Hb ligand binding and allosteric function.

The evolution of phenotypic similarities between species, known as convergence, illustrates that populations can respond predictably to ecological challenges. Convergence often results from similar genetic changes, which can emerge in two ways: the evolution of similar or identical mutations in independent lineages, which is termed parallel evolution; and the evolution in independent lineages of alleles that are shared among populations, which I call collateral genetic evolution. Evidence for parallel and collateral evolution has been found in many taxa, and an emerging hypothesis is that they result from the fact that mutations in some genetic targets minimize pleiotropic effects while simultaneously maximizing adaptation. If this proves correct, then the molecular changes underlying adaptation might be more predictable than has been appreciated previously.

A quantitative description of animal social behaviour is informative for behavioural biologists and clinicians developing drugs to treat social disorders. Social interaction in a group of animals has been difficult to measure because behaviour develops over long periods of time and requires tedious manual scoring, which is subjective and often non-reproducible. Computer-vision systems with the ability to measure complex social behaviour automatically would have a transformative impact on biology. Here, we present a method for tracking group-housed mice individually as they freely interact over multiple days. Each mouse is bleach-marked with a unique fur pattern. The patterns are automatically learned by the tracking software and used to infer identities. Trajectories are analysed to measure behaviour as it develops over days, beyond the range of acute experiments. We demonstrate how our system may be used to study the development of place preferences, associations and social relationships by tracking four mice continuously for five days. Our system enables accurate and reproducible characterisation of wild-type mouse social behaviour and paves the way for high-throughput long-term observation of the effects of genetic, pharmacological and environmental manipulations.