@article{3735, abstract = {The visual system is challenged with extracting and representing behaviorally relevant information contained in natural inputs of great complexity and detail. This task begins in the sensory periphery: retinal receptive fields and circuits are matched to the first and second-order statistical structure of natural inputs. This matching enables the retina to remove stimulus components that are predictable (and therefore uninformative), and primarily transmit what is unpredictable (and therefore informative). Here we show that this design principle applies to more complex aspects of natural scenes, and to central visual processing. We do this by classifying high-order statistics of natural scenes according to whether they are uninformative vs. informative. We find that the uninformative ones are perceptually nonsalient, while the informative ones are highly salient, and correspond to previously identified perceptual mechanisms whose neural basis is likely central. Our results suggest that the principle of efficient coding not only accounts for filtering operations in the sensory periphery, but also shapes subsequent stages of sensory processing that are sensitive to high-order image statistics.}, author = {Gasper Tkacik and Prentice, Jason S and Victor,Jonathan D and Balasubramanian, Vijay}, journal = {PNAS}, number = {42}, pages = {18149 -- 18154}, publisher = {National Academy of Sciences}, title = {{Local statistics in natural scenes predict the saliency of synthetic textures}}, doi = {10.1073/pnas.0914916107}, volume = {107}, year = {2010}, } @article{3736, abstract = {In retina and in cortical slice the collective response of spiking neural populations is well described by "maximum-entropy" models in which only pairs of neurons interact. We asked, how should such interactions be organized to maximize the amount of information represented in population responses? To this end, we extended the linear-nonlinear-Poisson model of single neural response to include pairwise interactions, yielding a stimulus-dependent, pairwise maximum-entropy model. We found that as we varied the noise level in single neurons and the distribution of network inputs, the optimal pairwise interactions smoothly interpolated to achieve network functions that are usually regarded as discrete–stimulus decorrelation, error correction, and independent encoding. These functions reflected a trade-off between efficient consumption of finite neural bandwidth and the use of redundancy to mitigate noise. Spontaneous activity in the optimal network reflected stimulus-induced activity patterns, and single-neuron response variability overestimated network noise. Our analysis suggests that rather than having a single coding principle hardwired in their architecture, networks in the brain should adapt their function to changing noise and stimulus correlations.}, author = {Gasper Tkacik and Prentice, Jason S and Balasubramanian, Vijay and Schneidman, Elad}, journal = {PNAS}, number = {32}, pages = {14419 -- 14424}, publisher = {National Academy of Sciences}, title = {{Optimal population coding by noisy spiking neurons}}, doi = {10.1073/pnas.1004906107}, volume = {107}, year = {2010}, } @article{3738, abstract = {Central to the functioning of a living cell is its ability to control the readout or expression of information encoded in the genome. In many cases, a single transcription factor protein activates or represses the expression of many genes. As the concentration of the transcription factor varies, the target genes thus undergo correlated changes, and this redundancy limits the ability of the cell to transmit information about input signals. We explore how interactions among the target genes can reduce this redundancy and optimize information transmission. Our discussion builds on recent work [Tkacik, Phys. Rev. E 80, 031920 (2009)], and there are connections to much earlier work on the role of lateral inhibition in enhancing the efficiency of information transmission in neural circuits; for simplicity we consider here the case where the interactions have a feed forward structure, with no loops. Even with this limitation, the networks that optimize information transmission have a structure reminiscent of the networks found in real biological systems.}, author = {Walczak, Aleksandra M and Gasper Tkacik and Bialek, William S}, journal = {Physical Review E Statistical Nonlinear and Soft Matter Physics}, number = {4}, publisher = {American Institute of Physics}, title = {{Optimizing information flow in small genetic networks. II. Feed-forward interactions}}, doi = {10.1103/PhysRevE.81.041905}, volume = {81}, year = {2010}, } @unpublished{3743, abstract = {These are notes for a set of 7 two-hour lectures given at the 2010 Summer School on Quantitative Evolutionary and Comparative Genomics at OIST, Okinawa, Japan. The emphasis is on understanding how biological systems process information. We take a physicist's approach of looking for simple phenomenological descriptions that can address the questions of biological function without necessarily modeling all (mostly unknown) microscopic details; the example that is developed throughout the notes is transcriptional regulation in genetic regulatory networks. We present tools from information theory and statistical physics that can be used to analyze noisy nonlinear biological networks, and build generative and predictive models of regulatory processes.}, author = {Gasper Tkacik}, booktitle = {ArXiv}, pages = {1 -- 52}, publisher = {ArXiv}, title = {{From statistical mechanics to information theory: understanding biophysical information-processing systems}}, volume = {q-bio.MN}, year = {2010}, } @article{3748, abstract = {The chemotaxis signalling network in Escherichia coli that controls the locomotion of bacteria is a classic model system for signal transduction1, 2. This pathway modulates the behaviour of flagellar motors to propel bacteria towards sources of chemical attractants. Although this system relaxes to a steady state in response to environmental changes, the signalling events within the chemotaxis network are noisy and cause large temporal variations of the motor behaviour even in the absence of stimulus3. That the same signalling network governs both behavioural variability and cellular response raises the question of whether these two traits are independent. Here, we experimentally establish a fluctuation–response relationship in the chemotaxis system of living bacteria. Using this relationship, we demonstrate the possibility of inferring the cellular response from the behavioural variability measured before stimulus. In monitoring the pre- and post-stimulus switching behaviour of individual bacterial motors, we found that variability scales linearly with the response time for different functioning states of the cell. This study highlights that the fundamental relationship between fluctuation and response is not constrained to physical systems at thermodynamic equilibrium4 but is extensible to living cells5. Such a relationship not only implies that behavioural variability and cellular response can be coupled traits, but it also provides a general framework within which we can examine how the selection of a network design shapes this interdependence}, author = {Park, Heungwon and Pontius, William and Calin Guet and Marko, John F and Emonet,Thierry and Cluzel,Philippe}, journal = {Nature}, pages = {819 -- 823}, publisher = {Nature Publishing Group}, title = {{Interdependence of behavioural variability and response to small stimuli in bacteria}}, doi = {10.1038/nature09551}, volume = {468}, year = {2010}, } @article{3749, abstract = {In E. coli, chemotactic behavior exhibits perfect adaptation that is robust to changes in the intracellular concentration of the chemotactic proteins, such as CheR and CheB. However, the robustness of the perfect adaptation does not explicitly imply a robust chemotactic response. Previous studies on the robustness of the chemotactic response relied on swarming assays, which can be confounded by processes besides chemotaxis, such as cellular growth and depletion of nutrients. Here, using a high-throughput capillary assay that eliminates the effects of growth, we experimentally studied how the chemotactic response depends on the relative concentration of the chemotactic proteins. We simultaneously measured both the chemotactic response of E. coli cells to L: -aspartate and the concentrations of YFP-CheR and CheB-CFP fusion proteins. We found that the chemotactic response is fine-tuned to a specific ratio of [CheR]/[CheB] with a maximum response comparable to the chemotactic response of wild-type behavior. In contrast to adaptation in chemotaxis, that is robust and exact, capillary assays revealed that the chemotactic response in swimming bacteria is fined-tuned to wild-type level of the [CheR]/[CheB] ratio.}, author = {Park, Heungwon and Calin Guet and Emonet,Thierry and Cluzel,Philippe}, journal = {Current Microbiology}, number = {3}, pages = {764 -- 769}, publisher = {Springer}, title = {{Fine-tuning of chemotactic response in E. coli determined by high-throughput capillary assay}}, doi = {10.1007/s00284-010-9778-z}, volume = {62}, year = {2010}, } @article{3759, abstract = {We propose a mesh-based surface tracking method for fluid animation that both preserves fine surface details and robustly adjusts the topology of the surface in the presence of arbitrarily thin features like sheets and strands. We replace traditional re-sampling methods with a convex hull method for connecting surface features during topological changes. This technique permits arbitrarily thin fluid features with minimal re-sampling errors by reusing points from the original surface. We further reduce re-sampling artifacts with a subdivision-based mesh-stitching algorithm, and we use a higher order interpolating subdivision scheme to determine the location of any newly-created vertices. The resulting algorithm efficiently produces detailed fluid surfaces with arbitrarily thin features while maintaining a consistent topology with the underlying fluid simulation.}, author = {Wojtan, Christopher J and Thürey, Nils and Gross, Markus and Turk, Greg}, journal = {ACM Transactions on Graphics}, number = {4}, publisher = {ACM}, title = {{Physics-inspired topology changes for thin fluid features}}, doi = {10.1145/1778765.1778787}, volume = {29}, year = {2010}, } @article{3761, abstract = {We present an algorithm for creating realistic animations of characters that are swimming through fluids. Our approach combines dynamic simulation with data-driven kinematic motions (motion capture data) to produce realistic animation in a fluid. The interaction of the articulated body with the fluid is performed by incorporating joint constraints with rigid animation and by extending a solid/fluid coupling method to handle articulated chains. Our solver takes as input the current state of the simulation and calculates the angular and linear accelerations of the connected bodies needed to match a particular motion sequence for the articulated body. These accelerations are used to estimate the forces and torques that are then applied to each joint. Based on this approach, we demonstrate simulated swimming results for a variety of different strokes, including crawl, backstroke, breaststroke, and butterfly. The ability to have articulated bodies interact with fluids also allows us to generate simulations of simple water creatures that are driven by simple controllers.}, author = {Kwatra, Nipun and Wojtan, Christopher J and Carlson, Mark and Essa, Irfan and Mucha, Peter and Turk, Greg}, journal = {IEEE Transactions on Visualization and Computer Graphics}, number = {1}, pages = {70 -- 80}, publisher = {IEEE}, title = {{Fluid simulation with articulated bodies}}, doi = {10.1109/TVCG.2009.66}, volume = {16}, year = {2010}, } @article{3766, abstract = {We present an approach to simulate flows driven by surface tension based on triangle meshes. Our method consists of two simulation layers: the first layer is an Eulerian method for simulating surface tension forces that is free from typical strict time step constraints. The second simulation layer is a Lagrangian finite element method that simulates sub-grid scale wave details on the fluid surface. The surface wave simulation employs an unconditionally stable, symplectic time integration method that allows for a high propagation speed due to strong surface tension. Our approach can naturally separate the grid-and sub-grid scales based on a volumepreserving mean curvature flow. As our model for the sub-grid dynamics enforces a local conservation of mass, it leads to realistic pinch off and merging effects. In addition to this method for simulating dynamic surface tension effects, we also present an efficient non-oscillatory approximation for capturing damped surface tension behavior. These approaches allow us to efficiently simulate complex phenomena associated with strong surface tension, such as Rayleigh-Plateau instabilities and crown splashes, in a short amount of time.}, author = {Thürey, Nils and Wojtan, Christopher J and Gross, Markus and Turk, Greg}, journal = {ACM Transactions on Graphics}, number = {4}, publisher = {ACM}, title = {{A multiscale approach to mesh-based surface tension flows}}, doi = {10.1145/1778765.1778785}, volume = {29}, year = {2010}, } @article{3772, author = {Barton, Nicholas H}, journal = {PLoS Genetics}, number = {6}, publisher = {Public Library of Science}, title = {{Understanding adaptation in large populations}}, doi = {10.1371/journal.pgen.1000987}, volume = {6}, year = {2010}, } @article{3773, abstract = {If distinct biological species are to coexist in sympatry, they must be reproductively isolated and must exploit different limiting resources. A two-niche Levene model is analysed, in which habitat preference and survival depend on underlying additive traits. The population genetics of preference and viability are equivalent. However, there is a linear trade-off between the chances of settling in either niche, whereas viabilities may be constrained arbitrarily. With a convex trade-off, a sexual population evolves a single generalist genotype, whereas with a concave trade-off, disruptive selection favours maximal variance. A pure habitat preference evolves to global linkage equilibrium if mating occurs in a single pool, but remarkably, evolves to pairwise linkage equilibrium within niches if mating is within those niches--independent of the genetics. With a concave trade-off, the population shifts sharply between a unimodal distribution with high gene flow and a bimodal distribution with strong isolation, as the underlying genetic variance increases. However, these alternative states are only simultaneously stable for a narrow parameter range. A sharp threshold is only seen if survival in the 'wrong' niche is low; otherwise, strong isolation is impossible. Gene flow from divergent demes makes speciation much easier in parapatry than in sympatry.}, author = {Barton, Nicholas H}, journal = {Philosophical Transactions of the Royal Society of London. Series B, Biological Sciences}, number = {1547}, pages = {1825 -- 1840}, publisher = {Royal Society}, title = {{What role does natural selection play in speciation?}}, doi = {10.1098/rstb.2010.0001}, volume = {365}, year = {2010}, } @article{3774, abstract = {1. Hybridisation with an invasive species has the potential to alter the phenotype and hence the ecology of a native counterpart. 2. Here data from populations of native red deer Cervus elaphus and invasive sika deer Cervus nippon in Scotland is used to assess the extent to which hybridisation between them is causing phenotypic change. This is done by regression of phenotypic traits against genetic hybrid scores. 3. Hybridisation is causing increases in the body weight of sika-like deer and decreases in the body weight of red-like females. Hybridisation is causing increases in jaw length and increases in incisor arcade breadth in sika-like females. Hybridisation is also causing decreases in incisor arcade breadth in red-like females. 4. There is currently no evidence that hybridisation is causing changes in the kidney fat weight or pregnancy rates of either population. 5. Increased phenotypic similarity between the two species is likely to lead to further hybridisation. The ecological consequences of this are difficult to predict.}, author = {Senn, Helen and Swanson, Graeme and Goodman, Simon and Barton, Nicholas H and Pemberton, Josephine}, journal = {Journal of Animal Ecology}, number = {2}, pages = {414 -- 425}, publisher = {Wiley-Blackwell}, title = {{Phenotypic correlates of hybridisation between red and sika deer (genus Cervus)}}, doi = {10.1111/j.1365-2656.2009.01633.x}, volume = {79}, year = {2010}, } @article{3776, abstract = {The prevalence of recombination in eukaryotes poses one of the most puzzling questions in biology. The most compelling general explanation is that recombination facilitates selection by breaking down the negative associations generated by random drift (i.e. Hill-Robertson interference, HRI). I classify the effects of HRI owing to: deleterious mutation, balancing selection and selective sweeps on: neutral diversity, rates of adaptation and the mutation load. These effects are mediated primarily by the density of deleterious mutations and of selective sweeps. Sequence polymorphism and divergence suggest that these rates may be high enough to cause significant interference even in genomic regions of high recombination. However, neither seems able to generate enough variance in fitness to select strongly for high rates of recombination. It is plausible that spatial and temporal fluctuations in selection generate much more fitness variance, and hence selection for recombination, than can be explained by uniformly deleterious mutations or species-wide selective sweeps.}, author = {Barton, Nicholas H}, journal = {Philosophical Transactions of the Royal Society of London. Series B, Biological Sciences}, number = {1552}, pages = {2559 -- 2569}, publisher = {Royal Society}, title = {{Genetic linkage and natural selection}}, doi = {10.1098/rstb.2010.0106}, volume = {365}, year = {2010}, } @article{3777, abstract = {Under the classical view, selection depends more or less directly on mutation: standing genetic variance is maintained by a balance between selection and mutation, and adaptation is fuelled by new favourable mutations. Recombination is favoured if it breaks negative associations among selected alleles, which interfere with adaptation. Such associations may be generated by negative epistasis, or by random drift (leading to the Hill-Robertson effect). Both deterministic and stochastic explanations depend primarily on the genomic mutation rate, U. This may be large enough to explain high recombination rates in some organisms, but seems unlikely to be so in general. Random drift is a more general source of negative linkage disequilibria, and can cause selection for recombination even in large populations, through the chance loss of new favourable mutations. The rate of species-wide substitutions is much too low to drive this mechanism, but local fluctuations in selection, combined with gene flow, may suffice. These arguments are illustrated by comparing the interaction between good and bad mutations at unlinked loci under the infinitesimal model.}, author = {Barton, Nicholas H}, journal = {Philosophical Transactions of the Royal Society of London. Series B, Biological Sciences}, number = {1544}, pages = {1281 -- 1294}, publisher = {Royal Society}, title = {{Mutation and the evolution of recombination}}, doi = {10.1098/rstb.2009.0320}, volume = {365}, year = {2010}, } @article{3779, abstract = {Crosses between closely related species give two contrasting results. One result is that species hybrids may be inferior to their parents, for example, being less fertile [1]. The other is that F1 hybrids may display superior performance (heterosis), for example with increased vigour [2]. Although various hypotheses have been proposed to account for these two aspects of hybridisation, their biological basis is still poorly understood [3]. To gain further insights into this issue, we analysed the role that variation in gene expression may play. We took a conserved trait, flower asymmetry in Antirrhinum, and determined the extent to which the underlying regulatory genes varied in expression among closely related species. We show that expression of both genes analysed, CYC and RAD, varies significantly between species because of cis-acting differences. By making a quantitative genotype-phenotype map, using a range of mutant alleles, we demonstrate that the species lie on a plateau in gene expression-morphology space, so that the variation has no detectable phenotypic effect. However, phenotypic differences can be revealed by shifting genotypes off the plateau through genetic crosses. Our results can be readily explained if genomes are free to evolve within an effectively neutral zone in gene expression space. The consequences of this drift will be negligible for individual loci, but when multiple loci across the genome are considered, we show that the variation may have significant effects on phenotype and fitness, causing a significant drift load. By considering these consequences for various gene-expression-fitness landscapes, we conclude that F1 hybrids might be expected to show increased performance with regard to conserved traits, such as basic physiology, but reduced performance with regard to others. Thus, our study provides a new way of explaining how various aspects of hybrid performance may arise through natural variation in gene activity.}, author = {Rosas, Ulises and Barton, Nicholas H and Copsey, Lucy and Barbier De Reuille, Pierre and Coen, Enrico}, journal = {PLoS Biology}, number = {7}, publisher = {Public Library of Science}, title = {{Cryptic variation between species and the basis of hybrid performance}}, doi = {10.1371/journal.pbio.1000429}, volume = {8}, year = {2010}, } @inproceedings{3782, abstract = {In cortex surface segmentation, the extracted surface is required to have a particular topology, namely, a two-sphere. We present a new method for removing topology noise of a curve or surface within the level set framework, and thus produce a cortical surface with correct topology. We define a new energy term which quantifies topology noise. We then show how to minimize this term by computing its functional derivative with respect to the level set function. This method differs from existing methods in that it is inherently continuous and not digital; and in the way that our energy directly relates to the topology of the underlying curve or surface, versus existing knot-based measures which are related in a more indirect fashion. The proposed flow is validated empirically.}, author = {Chen, Chao and Freedman, Daniel}, booktitle = { Conference proceedings MCV 2010}, location = {Beijing, China}, pages = {31 -- 42}, publisher = {Springer}, title = {{Topology noise removal for curve and surface evolution}}, doi = {10.1007/978-3-642-18421-5_4}, volume = {6533}, year = {2010}, } @article{3783, abstract = {MICROSATELIGHT is a Perl/Tk pipeline with a graphical user interface that facilitates several tasks when scoring microsatellites. It implements new subroutines in R and PERL and takes advantage of features provided by previously developed freeware. MICROSATELIGHT takes raw genotype data and automates the peak identification through PeakScanner. The PeakSelect subroutine assigns peaks to different microsatellite markers according to their multiplex group, fluorochrome type, and size range. After peak selection, binning of alleles can be carried out 1) automatically through AlleloBin or 2) by manual bin definition through Binator. In both cases, several features for quality checking and further binning improvement are provided. The genotype table can then be converted into input files for several population genetics programs through CREATE. Finally, Hardy–Weinberg equilibrium tests and confidence intervals for null allele frequency can be obtained through GENEPOP. MICROSATELIGHT is the only freely available public-domain software that facilitates full multiplex microsatellite scoring, from electropherogram files to user-defined text files to be used with population genetics software. MICROSATELIGHT has been created for the Windows XP operating system and has been successfully tested under Windows 7. It is available at http://sourceforge.net/projects/microsatelight/.}, author = {Palero, Ferran and González Candelas, Fernando and Pascual, Marta}, journal = {Journal of Heredity}, number = {2}, pages = {247 -- 249}, publisher = {Oxford University Press}, title = {{Microsatelight – Pipeline to expedite microsatellite analysis}}, doi = {10.1093/jhered/esq111}, volume = {102}, year = {2010}, } @article{3785, abstract = {Most fisheries involving spiny lobsters of the genus Palinurus have been over exploited during the last decades, so there is a raising concern about management decisions for these valuable resources. A total of 13 microsatellite DNA loci recently developed in Palinurus elephas were assayed in order to assess genetic diversity levels in every known species of the genus. Microsatellite markers gave amplifications and showed polymorphism in all species, with gene diversity values varying from 0.65060.077 SD (Palinurus barbarae) to 0.79260.051 SD (Palinurus elephas). Most importantly, when depth distribution was taken into account, shallower waters pecies consistently showed larger historical effective population sizes than their deeper-water counterparts. This could explain why deeper-water species are more sensitive to overfishing, and would indicate that overexploitation may have a larger impact on their long-term genetic diversity.}, author = {Palero, Ferran and Abello, Pere and Macpherson, E. and Matthee, C. and Pascual, Marta}, issn = {1937-240X}, journal = {Journal of Crustacean Biology}, number = {4}, pages = {658 -- 663}, publisher = {Oxford University Press}, title = {{Genetic diversity levels in fishery-exploited spiny lobsters of the Genus Palinurus (Decapoda: Achelata)}}, doi = {10.1651/09-3192.1}, volume = {30}, year = {2010}, } @article{3786, abstract = {Four rare palinurid phyllosoma larvae, one mid-stage and three final stage, were found among the unclassified collections in the Crustacea Section, Natural History Museum, London. Detailed morphological analysis of the larvae indicated that they belong to several Palinustus species given the presence of incipient blunt-horns, length of antennula, length ratio of segments of antennular peduncle, distribution of pereiopod spines, and shape of uropods and telson. Moreover, the size of the final-stage larvae agrees with that expected given the size of the recently described puerulus stage of Palinustus mossambicus. This constitutes the first description of a complete phyllosoma assigned to Palinustus species. The phyllosoma described in the present study include the largest Palinuridae larva ever found.}, author = {Palero, Ferran and Guerao, Guillermo and Clark, Paul and Abello, Pere}, journal = {Zootaxa}, number = {1}, pages = {42 -- 58}, publisher = {Magnolia Press}, title = {{Final-stage phyllosoma of Palinustus A. Milne-Edwards, 1880 (Crustacea: Decapoda: Achelata: Palinuridae)-The first complete description}}, doi = {10.11646/zootaxa.2403.1.4}, volume = {2403}, year = {2010}, } @article{3787, abstract = {DNA samples were extracted from ethanol and formalin-fixed decapod crustacean tissue using a new method based on Tetramethylsilane (TMS)-Chelex. It is shown that neither an indigestible matrix of cross-linked protein nor soluble PCR inhibitors impede PCR success when dealing with formalin-fixed material. Instead, amplification success from formalin-fixed tissue appears to depend on the presence of unmodified DNA in the extracted sample. A staining method that facilitates the targeting of samples with a high content of unmodified DNA is provided.}, author = {Palero, Ferran and Hall, Sally and Clark, Paul and Johnston, David and Mackenzie Dodds, Jackie and Thatje, Sven}, journal = {Scientia Marina}, number = {3}, pages = {465 -- 470}, publisher = {Consejo Superior de Investigaciones Científicas}, title = {{DNA extraction from formalin-fixed tissue: new light from the deep sea}}, doi = {10.3989/scimar.2010.74n3465}, volume = {74}, year = {2010}, }