@article{157, abstract = {Social dilemmas occur when incentives for individuals are misaligned with group interests 1-7 . According to the 'tragedy of the commons', these misalignments can lead to overexploitation and collapse of public resources. The resulting behaviours can be analysed with the tools of game theory 8 . The theory of direct reciprocity 9-15 suggests that repeated interactions can alleviate such dilemmas, but previous work has assumed that the public resource remains constant over time. Here we introduce the idea that the public resource is instead changeable and depends on the strategic choices of individuals. An intuitive scenario is that cooperation increases the public resource, whereas defection decreases it. Thus, cooperation allows the possibility of playing a more valuable game with higher payoffs, whereas defection leads to a less valuable game. We analyse this idea using the theory of stochastic games 16-19 and evolutionary game theory. We find that the dependence of the public resource on previous interactions can greatly enhance the propensity for cooperation. For these results, the interaction between reciprocity and payoff feedback is crucial: neither repeated interactions in a constant environment nor single interactions in a changing environment yield similar cooperation rates. Our framework shows which feedbacks between exploitation and environment - either naturally occurring or designed - help to overcome social dilemmas.}, author = {Hilbe, Christian and Šimsa, Štepán and Chatterjee, Krishnendu and Nowak, Martin}, journal = {Nature}, number = {7713}, pages = {246 -- 249}, publisher = {Nature Publishing Group}, title = {{Evolution of cooperation in stochastic games}}, doi = {10.1038/s41586-018-0277-x}, volume = {559}, year = {2018}, } @article{384, abstract = {Can orthologous proteins differ in terms of their ability to be secreted? To answer this question, we investigated the distribution of signal peptides within the orthologous groups of Enterobacterales. Parsimony analysis and sequence comparisons revealed a large number of signal peptide gain and loss events, in which signal peptides emerge or disappear in the course of evolution. Signal peptide losses prevail over gains, an effect which is especially pronounced in the transition from the free-living or commensal to the endosymbiotic lifestyle. The disproportionate decline in the number of signal peptide-containing proteins in endosymbionts cannot be explained by the overall reduction of their genomes. Signal peptides can be gained and lost either by acquisition/elimination of the corresponding N-terminal regions or by gradual accumulation of mutations. The evolutionary dynamics of signal peptides in bacterial proteins represents a powerful mechanism of functional diversification.}, author = {Hönigschmid, Peter and Bykova, Nadya and Schneider, René and Ivankov, Dmitry and Frishman, Dmitrij}, journal = {Genome Biology and Evolution}, number = {3}, pages = {928 -- 938}, publisher = {Oxford University Press}, title = {{Evolutionary interplay between symbiotic relationships and patterns of signal peptide gain and loss}}, doi = {10.1093/gbe/evy049}, volume = {10}, year = {2018}, } @article{563, abstract = {In continuous populations with local migration, nearby pairs of individuals have on average more similar genotypes than geographically well separated pairs. A barrier to gene flow distorts this classical pattern of isolation by distance. Genetic similarity is decreased for sample pairs on different sides of the barrier and increased for pairs on the same side near the barrier. Here, we introduce an inference scheme that utilizes this signal to detect and estimate the strength of a linear barrier to gene flow in two-dimensions. We use a diffusion approximation to model the effects of a barrier on the geographical spread of ancestry backwards in time. This approach allows us to calculate the chance of recent coalescence and probability of identity by descent. We introduce an inference scheme that fits these theoretical results to the geographical covariance structure of bialleleic genetic markers. It can estimate the strength of the barrier as well as several demographic parameters. We investigate the power of our inference scheme to detect barriers by applying it to a wide range of simulated data. We also showcase an example application to a Antirrhinum majus (snapdragon) flower color hybrid zone, where we do not detect any signal of a strong genome wide barrier to gene flow.}, author = {Ringbauer, Harald and Kolesnikov, Alexander and Field, David and Barton, Nicholas H}, journal = {Genetics}, number = {3}, pages = {1231--1245}, publisher = {Genetics Society of America}, title = {{Estimating barriers to gene flow from distorted isolation-by-distance patterns}}, doi = {10.1534/genetics.117.300638}, volume = {208}, year = {2018}, } @article{135, abstract = {The Fluid Implicit Particle method (FLIP) reduces numerical dissipation by combining particles with grids. To improve performance, the subsequent narrow band FLIP method (NB‐FLIP) uses a FLIP‐based fluid simulation only near the liquid surface and a traditional grid‐based fluid simulation away from the surface. This spatially‐limited FLIP simulation significantly reduces the number of particles and alleviates a computational bottleneck. In this paper, we extend the NB‐FLIP idea even further, by allowing a simulation to transition between a FLIP‐like fluid simulation and a grid‐based simulation in arbitrary locations, not just near the surface. This approach leads to even more savings in memory and computation, because we can concentrate the particles only in areas where they are needed. More importantly, this new method allows us to seamlessly transition to smooth implicit surface geometry wherever the particle‐based simulation is unnecessary. Consequently, our method leads to a practical algorithm for avoiding the noisy surface artifacts associated with particle‐based liquid simulations, while simultaneously maintaining the benefits of a FLIP simulation in regions of dynamic motion.}, author = {Sato, Takahiro and Wojtan, Christopher J and Thuerey, Nils and Igarashi, Takeo and Ando, Ryoichi}, issn = {0167-7055}, journal = {Computer Graphics Forum}, number = {2}, pages = {169 -- 177}, publisher = {Wiley}, title = {{Extended narrow band FLIP for liquid simulations}}, doi = {10.1111/cgf.13351}, volume = {37}, year = {2018}, } @article{316, abstract = {Self-incompatibility (SI) is a genetically based recognition system that functions to prevent self-fertilization and mating among related plants. An enduring puzzle in SI is how the high diversity observed in nature arises and is maintained. Based on the underlying recognition mechanism, SI can be classified into two main groups: self- and non-self recognition. Most work has focused on diversification within self-recognition systems despite expected differences between the two groups in the evolutionary pathways and outcomes of diversification. Here, we use a deterministic population genetic model and stochastic simulations to investigate how novel S-haplotypes evolve in a gametophytic non-self recognition (SRNase/S Locus F-box (SLF)) SI system. For this model the pathways for diversification involve either the maintenance or breakdown of SI and can vary in the order of mutations of the female (SRNase) and male (SLF) components. We show analytically that diversification can occur with high inbreeding depression and self-pollination, but this varies with evolutionary pathway and level of completeness (which determines the number of potential mating partners in the population), and in general is more likely for lower haplotype number. The conditions for diversification are broader in stochastic simulations of finite population size. However, the number of haplotypes observed under high inbreeding and moderate to high self-pollination is less than that commonly observed in nature. Diversification was observed through pathways that maintain SI as well as through self-compatible intermediates. Yet the lifespan of diversified haplotypes was sensitive to their level of completeness. By examining diversification in a non-self recognition SI system, this model extends our understanding of the evolution and maintenance of haplotype diversity observed in a self recognition system common in flowering plants.}, author = {Bodova, Katarina and Priklopil, Tadeas and Field, David and Barton, Nicholas H and Pickup, Melinda}, journal = {Genetics}, number = {3}, pages = {861--883}, publisher = {Genetics Society of America}, title = {{Evolutionary pathways for the generation of new self-incompatibility haplotypes in a non-self recognition system}}, doi = {10.1534/genetics.118.300748}, volume = {209}, year = {2018}, } @article{190, abstract = {The German cockroach, Blattella germanica, is a worldwide pest that infests buildings, including homes, restaurants, and hospitals, often living in unsanitary conditions. As a disease vector and producer of allergens, this species has major health and economic impacts on humans. Factors contributing to the success of the German cockroach include its resistance to a broad range of insecticides, immunity to many pathogens, and its ability, as an extreme generalist omnivore, to survive on most food sources. The recently published genome shows that B. germanica has an exceptionally high number of protein coding genes. In this study, we investigate the functions of the 93 significantly expanded gene families with the aim to better understand the success of B. germanica as a major pest despite such inhospitable conditions. We find major expansions in gene families with functions related to the detoxification of insecticides and allelochemicals, defense against pathogens, digestion, sensory perception, and gene regulation. These expansions might have allowed B. germanica to develop multiple resistance mechanisms to insecticides and pathogens, and enabled a broad, flexible diet, thus explaining its success in unsanitary conditions and under recurrent chemical control. The findings and resources presented here provide insights for better understanding molecular mechanisms that will facilitate more effective cockroach control.}, author = {Harrison, Mark and Arning, Nicolas and Kremer, Lucas and Ylla, Guillem and Belles, Xavier and Bornberg Bauer, Erich and Huylmans, Ann K and Jongepier, Evelien and Puilachs, Maria and Richards, Stephen and Schal, Coby}, journal = {Journal of Experimental Zoology Part B: Molecular and Developmental Evolution}, pages = {254--264}, publisher = {Wiley}, title = {{Expansions of key protein families in the German cockroach highlight the molecular basis of its remarkable success as a global indoor pest}}, doi = {10.1002/jez.b.22824}, volume = {330}, year = {2018}, } @article{404, abstract = {We construct martingale solutions to stochastic thin-film equations by introducing a (spatial) semidiscretization and establishing convergence. The discrete scheme allows for variants of the energy and entropy estimates in the continuous setting as long as the discrete energy does not exceed certain threshold values depending on the spatial grid size $h$. Using a stopping time argument to prolongate high-energy paths constant in time, arbitrary moments of coupled energy/entropy functionals can be controlled. Having established Hölder regularity of approximate solutions, the convergence proof is then based on compactness arguments---in particular on Jakubowski's generalization of Skorokhod's theorem---weak convergence methods, and recent tools on martingale convergence. }, author = {Fischer, Julian L and Grün, Günther}, journal = {SIAM Journal on Mathematical Analysis}, number = {1}, pages = {411 -- 455}, publisher = {Society for Industrial and Applied Mathematics }, title = {{Existence of positive solutions to stochastic thin-film equations}}, doi = {10.1137/16M1098796}, volume = {50}, year = {2018}, } @misc{9813, abstract = {File S1 contains figures that clarify the following features: (i) effect of population size on the average number/frequency of SI classes, (ii) changes in the minimal completeness deficit in time for a single class, and (iii) diversification diagrams for all studied pathways, including the summary figure for k = 8. File S2 contains the code required for a stochastic simulation of the SLF system with an example. This file also includes the output in the form of figures and tables.}, author = {Bod'ová, Katarína and Priklopil, Tadeas and Field, David and Barton, Nicholas H and Pickup, Melinda}, publisher = {Genetics Society of America}, title = {{Supplemental material for Bodova et al., 2018}}, doi = {10.25386/genetics.6148304.v1}, year = {2018}, } @article{5780, abstract = {Bioluminescence is found across the entire tree of life, conferring a spectacular set of visually oriented functions from attracting mates to scaring off predators. Half a dozen different luciferins, molecules that emit light when enzymatically oxidized, are known. However, just one biochemical pathway for luciferin biosynthesis has been described in full, which is found only in bacteria. Here, we report identification of the fungal luciferase and three other key enzymes that together form the biosynthetic cycle of the fungal luciferin from caffeic acid, a simple and widespread metabolite. Introduction of the identified genes into the genome of the yeast Pichia pastoris along with caffeic acid biosynthesis genes resulted in a strain that is autoluminescent in standard media. We analyzed evolution of the enzymes of the luciferin biosynthesis cycle and found that fungal bioluminescence emerged through a series of events that included two independent gene duplications. The retention of the duplicated enzymes of the luciferin pathway in nonluminescent fungi shows that the gene duplication was followed by functional sequence divergence of enzymes of at least one gene in the biosynthetic pathway and suggests that the evolution of fungal bioluminescence proceeded through several closely related stepping stone nonluminescent biochemical reactions with adaptive roles. The availability of a complete eukaryotic luciferin biosynthesis pathway provides several applications in biomedicine and bioengineering.}, author = {Kotlobay, Alexey A. and Sarkisyan, Karen and Mokrushina, Yuliana A. and Marcet-Houben, Marina and Serebrovskaya, Ekaterina O. and Markina, Nadezhda M. and Gonzalez Somermeyer, Louisa and Gorokhovatsky, Andrey Y. and Vvedensky, Andrey and Purtov, Konstantin V. and Petushkov, Valentin N. and Rodionova, Natalja S. and Chepurnyh, Tatiana V. and Fakhranurova, Liliia and Guglya, Elena B. and Ziganshin, Rustam and Tsarkova, Aleksandra S. and Kaskova, Zinaida M. and Shender, Victoria and Abakumov, Maxim and Abakumova, Tatiana O. and Povolotskaya, Inna S. and Eroshkin, Fedor M. and Zaraisky, Andrey G. and Mishin, Alexander S. and Dolgov, Sergey V. and Mitiouchkina, Tatiana Y. and Kopantzev, Eugene P. and Waldenmaier, Hans E. and Oliveira, Anderson G. and Oba, Yuichi and Barsova, Ekaterina and Bogdanova, Ekaterina A. and Gabaldón, Toni and Stevani, Cassius V. and Lukyanov, Sergey and Smirnov, Ivan V. and Gitelson, Josef I. and Kondrashov, Fyodor and Yampolsky, Ilia V.}, issn = {00278424}, journal = {Proceedings of the National Academy of Sciences of the United States of America}, number = {50}, pages = {12728--12732}, publisher = {National Academy of Sciences}, title = {{Genetically encodable bioluminescent system from fungi}}, doi = {10.1073/pnas.1803615115}, volume = {115}, year = {2018}, } @article{428, abstract = {The plant hormone gibberellic acid (GA) is a crucial regulator of growth and development. The main paradigm of GA signaling puts forward transcriptional regulation via the degradation of DELLA transcriptional repressors. GA has also been shown to regulate tropic responses by modulation of the plasma membrane incidence of PIN auxin transporters by an unclear mechanism. Here we uncovered the cellular and molecular mechanisms by which GA redirects protein trafficking and thus regulates cell surface functionality. Photoconvertible reporters revealed that GA balances the protein traffic between the vacuole degradation route and recycling back to the cell surface. Low GA levels promote vacuolar delivery and degradation of multiple cargos, including PIN proteins, whereas high GA levels promote their recycling to the plasma membrane. This GA effect requires components of the retromer complex, such as Sorting Nexin 1 (SNX1) and its interacting, microtubule (MT)-associated protein, the Cytoplasmic Linker-Associated Protein (CLASP1). Accordingly, GA regulates the subcellular distribution of SNX1 and CLASP1, and the intact MT cytoskeleton is essential for the GA effect on trafficking. This GA cellular action occurs through DELLA proteins that regulate the MT and retromer presumably via their interaction partners Prefoldins (PFDs). Our study identified a branching of the GA signaling pathway at the level of DELLA proteins, which, in parallel to regulating transcription, also target by a nontranscriptional mechanism the retromer complex acting at the intersection of the degradation and recycling trafficking routes. By this mechanism, GA can redirect receptors and transporters to the cell surface, thus coregulating multiple processes, including PIN-dependent auxin fluxes during tropic responses.}, author = {Salanenka, Yuliya and Verstraeten, Inge and Löfke, Christian and Tabata, Kaori and Naramoto, Satoshi and Glanc, Matous and Friml, Jirí}, journal = {PNAS}, number = {14}, pages = { 3716 -- 3721}, publisher = {National Academy of Sciences}, title = {{Gibberellin DELLA signaling targets the retromer complex to redirect protein trafficking to the plasma membrane}}, doi = {10.1073/pnas.1721760115}, volume = {115}, year = {2018}, } @article{62, abstract = {Imaging is a dominant strategy for data collection in neuroscience, yielding stacks of images that often scale to gigabytes of data for a single experiment. Machine learning algorithms from computer vision can serve as a pair of virtual eyes that tirelessly processes these images, automatically detecting and identifying microstructures. Unlike learning methods, our Flexible Learning-free Reconstruction of Imaged Neural volumes (FLoRIN) pipeline exploits structure-specific contextual clues and requires no training. This approach generalizes across different modalities, including serially-sectioned scanning electron microscopy (sSEM) of genetically labeled and contrast enhanced processes, spectral confocal reflectance (SCoRe) microscopy, and high-energy synchrotron X-ray microtomography (μCT) of large tissue volumes. We deploy the FLoRIN pipeline on newly published and novel mouse datasets, demonstrating the high biological fidelity of the pipeline’s reconstructions. FLoRIN reconstructions are of sufficient quality for preliminary biological study, for example examining the distribution and morphology of cells or extracting single axons from functional data. Compared to existing supervised learning methods, FLoRIN is one to two orders of magnitude faster and produces high-quality reconstructions that are tolerant to noise and artifacts, as is shown qualitatively and quantitatively.}, author = {Shabazi, Ali and Kinnison, Jeffery and Vescovi, Rafael and Du, Ming and Hill, Robert and Jösch, Maximilian A and Takeno, Marc and Zeng, Hongkui and Da Costa, Nuno and Grutzendler, Jaime and Kasthuri, Narayanan and Scheirer, Walter}, journal = {Scientific Reports}, number = {1}, publisher = {Nature Publishing Group}, title = {{Flexible learning-free segmentation and reconstruction of neural volumes}}, doi = {10.1038/s41598-018-32628-3}, volume = {8}, year = {2018}, } @article{437, abstract = {Dendritic cells (DCs) are sentinels of the adaptive immune system that reside in peripheral organs of mammals. Upon pathogen encounter, they undergo maturation and up-regulate the chemokine receptor CCR7 that guides them along gradients of its chemokine ligands CCL19 and 21 to the next draining lymph node. There, DCs present peripherally acquired antigen to naïve T cells, thereby triggering adaptive immunity.}, author = {Leithner, Alexander F and Renkawitz, Jörg and De Vries, Ingrid and Hauschild, Robert and Haecker, Hans and Sixt, Michael K}, journal = {European Journal of Immunology}, number = {6}, pages = {1074 -- 1077}, publisher = {Wiley-Blackwell}, title = {{Fast and efficient genetic engineering of hematopoietic precursor cells for the study of dendritic cell migration}}, doi = {10.1002/eji.201747358}, volume = {48}, year = {2018}, } @article{617, abstract = {Insects are exposed to a variety of potential pathogens in their environment, many of which can severely impact fitness and health. Consequently, hosts have evolved resistance and tolerance strategies to suppress or cope with infections. Hosts utilizing resistance improve fitness by clearing or reducing pathogen loads, and hosts utilizing tolerance reduce harmful fitness effects per pathogen load. To understand variation in, and selective pressures on, resistance and tolerance, we asked to what degree they are shaped by host genetic background, whether plasticity in these responses depends upon dietary environment, and whether there are interactions between these two factors. Females from ten wild-type Drosophila melanogaster genotypes were kept on high- or low-protein (yeast) diets and infected with one of two opportunistic bacterial pathogens, Lactococcus lactis or Pseudomonas entomophila. We measured host resistance as the inverse of bacterial load in the early infection phase. The relationship (slope) between fly fecundity and individual-level bacteria load provided our fecundity tolerance measure. Genotype and dietary yeast determined host fecundity and strongly affected survival after infection with pathogenic P. entomophila. There was considerable genetic variation in host resistance, a commonly found phenomenon resulting from for example varying resistance costs or frequency-dependent selection. Despite this variation and the reproductive cost of higher P. entomophila loads, fecundity tolerance did not vary across genotypes. The absence of genetic variation in tolerance may suggest that at this early infection stage, fecundity tolerance is fixed or that any evolved tolerance mechanisms are not expressed under these infection conditions.}, author = {Kutzer, Megan and Kurtz, Joachim and Armitage, Sophie}, issn = {1420-9101}, journal = {Journal of Evolutionary Biology}, number = {1}, pages = {159 -- 171}, publisher = {Wiley}, title = {{Genotype and diet affect resistance, survival, and fecundity but not fecundity tolerance}}, doi = {10.1111/jeb.13211}, volume = {31}, year = {2018}, } @article{5888, abstract = {Despite the remarkable number of scientific breakthroughs of the last 100 years, the treatment of neurodevelopmental disorders (e.g., autism spectrum disorder, intellectual disability) remains a great challenge. Recent advancements in genomics, such as whole-exome or whole-genome sequencing, have enabled scientists to identify numerous mutations underlying neurodevelopmental disorders. Given the few hundred risk genes that have been discovered, the etiological variability and the heterogeneous clinical presentation, the need for genotype — along with phenotype- based diagnosis of individual patients has become a requisite. In this review we look at recent advancements in genomic analysis and their translation into clinical practice.}, author = {Tarlungeanu, Dora-Clara and Novarino, Gaia}, issn = {2092-6413}, journal = {Experimental & Molecular Medicine}, number = {8}, publisher = {Springer Nature}, title = {{Genomics in neurodevelopmental disorders: an avenue to personalized medicine}}, doi = {10.1038/s12276-018-0129-7}, volume = {50}, year = {2018}, } @article{295, abstract = {We prove upper and lower bounds on the ground-state energy of the ideal two-dimensional anyon gas. Our bounds are extensive in the particle number, as for fermions, and linear in the statistics parameter (Formula presented.). The lower bounds extend to Lieb–Thirring inequalities for all anyons except bosons.}, author = {Lundholm, Douglas and Seiringer, Robert}, journal = {Letters in Mathematical Physics}, number = {11}, pages = {2523--2541}, publisher = {Springer}, title = {{Fermionic behavior of ideal anyons}}, doi = {10.1007/s11005-018-1091-y}, volume = {108}, year = {2018}, } @article{555, abstract = {Conventional wisdom has it that proteins fold and assemble into definite structures, and that this defines their function. Glycosaminoglycans (GAGs) are different. In most cases the structures they form have a low degree of order, even when interacting with proteins. Here, we discuss how physical features common to all GAGs — hydrophilicity, charge, linearity and semi-flexibility — underpin the overall properties of GAG-rich matrices. By integrating soft matter physics concepts (e.g. polymer brushes and phase separation) with our molecular understanding of GAG–protein interactions, we can better comprehend how GAG-rich matrices assemble, what their properties are, and how they function. Taking perineuronal nets (PNNs) — a GAG-rich matrix enveloping neurons — as a relevant example, we propose that microphase separation determines the holey PNN anatomy that is pivotal to PNN functions.}, author = {Richter, Ralf and Baranova, Natalia and Day, Anthony and Kwok, Jessica}, journal = {Current Opinion in Structural Biology}, pages = {65 -- 74}, publisher = {Elsevier}, title = {{Glycosaminoglycans in extracellular matrix organisation: Are concepts from soft matter physics key to understanding the formation of perineuronal nets?}}, doi = {10.1016/j.sbi.2017.12.002}, volume = {50}, year = {2018}, } @article{448, abstract = {Around 150 million years ago, eusocial termites evolved from within the cockroaches, 50 million years before eusocial Hymenoptera, such as bees and ants, appeared. Here, we report the 2-Gb genome of the German cockroach, Blattella germanica, and the 1.3-Gb genome of the drywood termite Cryptotermes secundus. We show evolutionary signatures of termite eusociality by comparing the genomes and transcriptomes of three termites and the cockroach against the background of 16 other eusocial and non-eusocial insects. Dramatic adaptive changes in genes underlying the production and perception of pheromones confirm the importance of chemical communication in the termites. These are accompanied by major changes in gene regulation and the molecular evolution of caste determination. Many of these results parallel molecular mechanisms of eusocial evolution in Hymenoptera. However, the specific solutions are remarkably different, thus revealing a striking case of convergence in one of the major evolutionary transitions in biological complexity.}, author = {Harrison, Mark and Jongepier, Evelien and Robertson, Hugh and Arning, Nicolas and Bitard Feildel, Tristan and Chao, Hsu and Childers, Christopher and Dinh, Huyen and Doddapaneni, Harshavardhan and Dugan, Shannon and Gowin, Johannes and Greiner, Carolin and Han, Yi and Hu, Haofu and Hughes, Daniel and Huylmans, Ann K and Kemena, Karsten and Kremer, Lukas and Lee, Sandra and López Ezquerra, Alberto and Mallet, Ludovic and Monroy Kuhn, Jose and Moser, Annabell and Murali, Shwetha and Muzny, Donna and Otani, Saria and Piulachs, Maria and Poelchau, Monica and Qu, Jiaxin and Schaub, Florentine and Wada Katsumata, Ayako and Worley, Kim and Xie, Qiaolin and Ylla, Guillem and Poulsen, Michael and Gibbs, Richard and Schal, Coby and Richards, Stephen and Belles, Xavier and Korb, Judith and Bornberg Bauer, Erich}, journal = {Nature Ecology and Evolution}, number = {3}, pages = {557--566}, publisher = {Springer Nature}, title = {{Hemimetabolous genomes reveal molecular basis of termite eusociality}}, doi = {10.1038/s41559-017-0459-1}, volume = {2}, year = {2018}, } @article{723, abstract = {Escaping local optima is one of the major obstacles to function optimisation. Using the metaphor of a fitness landscape, local optima correspond to hills separated by fitness valleys that have to be overcome. We define a class of fitness valleys of tunable difficulty by considering their length, representing the Hamming path between the two optima and their depth, the drop in fitness. For this function class we present a runtime comparison between stochastic search algorithms using different search strategies. The (1+1) EA is a simple and well-studied evolutionary algorithm that has to jump across the valley to a point of higher fitness because it does not accept worsening moves (elitism). In contrast, the Metropolis algorithm and the Strong Selection Weak Mutation (SSWM) algorithm, a famous process in population genetics, are both able to cross the fitness valley by accepting worsening moves. We show that the runtime of the (1+1) EA depends critically on the length of the valley while the runtimes of the non-elitist algorithms depend crucially on the depth of the valley. Moreover, we show that both SSWM and Metropolis can also efficiently optimise a rugged function consisting of consecutive valleys.}, author = {Oliveto, Pietro and Paixao, Tiago and Pérez Heredia, Jorge and Sudholt, Dirk and Trubenova, Barbora}, journal = {Algorithmica}, number = {5}, pages = {1604 -- 1633}, publisher = {Springer}, title = {{How to escape local optima in black box optimisation when non elitism outperforms elitism}}, doi = {10.1007/s00453-017-0369-2}, volume = {80}, year = {2018}, } @article{321, abstract = {The twelve papers in this special section focus on learning systems with shared information for computer vision and multimedia communication analysis. In the real world, a realistic setting for computer vision or multimedia recognition problems is that we have some classes containing lots of training data and many classes containing a small amount of training data. Therefore, how to use frequent classes to help learning rare classes for which it is harder to collect the training data is an open question. Learning with shared information is an emerging topic in machine learning, computer vision and multimedia analysis. There are different levels of components that can be shared during concept modeling and machine learning stages, such as sharing generic object parts, sharing attributes, sharing transformations, sharing regularization parameters and sharing training examples, etc. Regarding the specific methods, multi-task learning, transfer learning and deep learning can be seen as using different strategies to share information. These learning with shared information methods are very effective in solving real-world large-scale problems.}, author = {Darrell, Trevor and Lampert, Christoph and Sebe, Nico and Wu, Ying and Yan, Yan}, journal = {IEEE Transactions on Pattern Analysis and Machine Intelligence}, number = {5}, pages = {1029 -- 1031}, publisher = {IEEE}, title = {{Guest editors' introduction to the special section on learning with Shared information for computer vision and multimedia analysis}}, doi = {10.1109/TPAMI.2018.2804998}, volume = {40}, year = {2018}, } @misc{9841, abstract = {Around 150 million years ago, eusocial termites evolved from within the cockroaches, 50 million years before eusocial Hymenoptera, such as bees and ants, appeared. Here, we report the 2-Gb genome of the German cockroach, Blattella germanica, and the 1.3-Gb genome of the drywood termite Cryptotermes secundus. We show evolutionary signatures of termite eusociality by comparing the genomes and transcriptomes of three termites and the cockroach against the background of 16 other eusocial and non-eusocial insects. Dramatic adaptive changes in genes underlying the production and perception of pheromones confirm the importance of chemical communication in the termites. These are accompanied by major changes in gene regulation and the molecular evolution of caste determination. Many of these results parallel molecular mechanisms of eusocial evolution in Hymenoptera. However, the specific solutions are remarkably different, thus revealing a striking case of convergence in one of the major evolutionary transitions in biological complexity.}, author = {Harrison, Mark C. and Jongepier, Evelien and Robertson, Hugh M. and Arning, Nicolas and Bitard-Feildel, Tristan and Chao, Hsu and Childers, Christopher P. and Dinh, Huyen and Doddapaneni, Harshavardhan and Dugan, Shannon and Gowin, Johannes and Greiner, Carolin and Han, Yi and Hu, Haofu and Hughes, Daniel S. T. and Huylmans, Ann K and Kemena, Carsten and Kremer, Lukas P. M. and Lee, Sandra L. and Lopez-Ezquerra, Alberto and Mallet, Ludovic and Monroy-Kuhn, Jose M. and Moser, Annabell and Murali, Shwetha C. and Muzny, Donna M. and Otani, Saria and Piulachs, Maria-Dolors and Poelchau, Monica and Qu, Jiaxin and Schaub, Florentine and Wada-Katsumata, Ayako and Worley, Kim C. and Xie, Qiaolin and Ylla, Guillem and Poulsen, Michael and Gibbs, Richard A. and Schal, Coby and Richards, Stephen and Belles, Xavier and Korb, Judith and Bornberg-Bauer, Erich}, publisher = {Dryad}, title = {{Data from: Hemimetabolous genomes reveal molecular basis of termite eusociality}}, doi = {10.5061/dryad.51d4r}, year = {2018}, } @article{32, abstract = {The functional role of AMPA receptor (AMPAR)-mediated synaptic signaling between neurons and oligodendrocyte precursor cells (OPCs) remains enigmatic. We modified the properties of AMPARs at axon-OPC synapses in the mouse corpus callosum in vivo during the peak of myelination by targeting the GluA2 subunit. Expression of the unedited (Ca2+ permeable) or the pore-dead GluA2 subunit of AMPARs triggered proliferation of OPCs and reduced their differentiation into oligodendrocytes. Expression of the cytoplasmic C-terminal (GluA2(813-862)) of the GluA2 subunit (C-tail), a modification designed to affect the interaction between GluA2 and AMPAR-binding proteins and to perturb trafficking of GluA2-containing AMPARs, decreased the differentiation of OPCs without affecting their proliferation. These findings suggest that ionotropic and non-ionotropic properties of AMPARs in OPCs, as well as specific aspects of AMPAR-mediated signaling at axon-OPC synapses in the mouse corpus callosum, are important for balancing the response of OPCs to proliferation and differentiation cues. In the brain, oligodendrocyte precursor cells (OPCs) receive glutamatergic AMPA-receptor-mediated synaptic input from neurons. Chen et al. show that modifying AMPA-receptor properties at axon-OPC synapses alters proliferation and differentiation of OPCs. This expands the traditional view of synaptic transmission by suggesting neurons also use synapses to modulate behavior of glia.}, author = {Chen, Ting and Kula, Bartosz and Nagy, Balint and Barzan, Ruxandra and Gall, Andrea and Ehrlich, Ingrid and Kukley, Maria}, journal = {Cell Reports}, number = {4}, pages = {852 -- 861.e7}, publisher = {Elsevier}, title = {{In Vivo regulation of Oligodendrocyte processor cell proliferation and differentiation by the AMPA-receptor Subunit GluA2}}, doi = {10.1016/j.celrep.2018.09.066}, volume = {25}, year = {2018}, } @article{5672, abstract = {The release of IgM is the first line of an antibody response and precedes the generation of high affinity IgG in germinal centers. Once secreted by freshly activated plasmablasts, IgM is released into the efferent lymph of reactive lymph nodes as early as 3 d after immunization. As pentameric IgM has an enormous size of 1,000 kD, its diffusibility is low, and one might wonder how it can pass through the densely lymphocyte-packed environment of a lymph node parenchyma in order to reach its exit. In this issue of JEM, Thierry et al. show that, in order to reach the blood stream, IgM molecules take a specific micro-anatomical route via lymph node conduits.}, author = {Reversat, Anne and Sixt, Michael K}, issn = {00221007}, journal = {Journal of Experimental Medicine}, number = {12}, pages = {2959--2961}, publisher = {Rockefeller University Press}, title = {{IgM's exit route}}, doi = {10.1084/jem.20181934}, volume = {215}, year = {2018}, } @article{458, abstract = {We consider congruences of straight lines in a plane with the combinatorics of the square grid, with all elementary quadrilaterals possessing an incircle. It is shown that all the vertices of such nets (we call them incircular or IC-nets) lie on confocal conics. Our main new results are on checkerboard IC-nets in the plane. These are congruences of straight lines in the plane with the combinatorics of the square grid, combinatorially colored as a checkerboard, such that all black coordinate quadrilaterals possess inscribed circles. We show how this larger class of IC-nets appears quite naturally in Laguerre geometry of oriented planes and spheres and leads to new remarkable incidence theorems. Most of our results are valid in hyperbolic and spherical geometries as well. We present also generalizations in spaces of higher dimension, called checkerboard IS-nets. The construction of these nets is based on a new 9 inspheres incidence theorem.}, author = {Akopyan, Arseniy and Bobenko, Alexander}, journal = {Transactions of the American Mathematical Society}, number = {4}, pages = {2825 -- 2854}, publisher = {American Mathematical Society}, title = {{Incircular nets and confocal conics}}, doi = {10.1090/tran/7292}, volume = {370}, year = {2018}, } @inproceedings{5788, abstract = {In two-player games on graphs, the players move a token through a graph to produce an infinite path, which determines the winner or payoff of the game. Such games are central in formal verification since they model the interaction between a non-terminating system and its environment. We study bidding games in which the players bid for the right to move the token. Two bidding rules have been defined. In Richman bidding, in each round, the players simultaneously submit bids, and the higher bidder moves the token and pays the other player. Poorman bidding is similar except that the winner of the bidding pays the “bank” rather than the other player. While poorman reachability games have been studied before, we present, for the first time, results on infinite-duration poorman games. A central quantity in these games is the ratio between the two players’ initial budgets. The questions we study concern a necessary and sufficient ratio with which a player can achieve a goal. For reachability objectives, such threshold ratios are known to exist for both bidding rules. We show that the properties of poorman reachability games extend to complex qualitative objectives such as parity, similarly to the Richman case. Our most interesting results concern quantitative poorman games, namely poorman mean-payoff games, where we construct optimal strategies depending on the initial ratio, by showing a connection with random-turn based games. The connection in itself is interesting, because it does not hold for reachability poorman games. We also solve the complexity problems that arise in poorman bidding games.}, author = {Avni, Guy and Henzinger, Thomas A and Ibsen-Jensen, Rasmus}, isbn = {9783030046118}, issn = {03029743}, location = {Oxford, UK}, pages = {21--36}, publisher = {Springer}, title = {{Infinite-duration poorman-bidding games}}, doi = {10.1007/978-3-030-04612-5_2}, volume = {11316}, year = {2018}, } @article{150, abstract = {A short, 14-amino-acid segment called SP1, located in the Gag structural protein1, has a critical role during the formation of the HIV-1 virus particle. During virus assembly, the SP1 peptide and seven preceding residues fold into a six-helix bundle, which holds together the Gag hexamer and facilitates the formation of a curved immature hexagonal lattice underneath the viral membrane2,3. Upon completion of assembly and budding, proteolytic cleavage of Gag leads to virus maturation, in which the immature lattice is broken down; the liberated CA domain of Gag then re-assembles into the mature conical capsid that encloses the viral genome and associated enzymes. Folding and proteolysis of the six-helix bundle are crucial rate-limiting steps of both Gag assembly and disassembly, and the six-helix bundle is an established target of HIV-1 inhibitors4,5. Here, using a combination of structural and functional analyses, we show that inositol hexakisphosphate (InsP6, also known as IP6) facilitates the formation of the six-helix bundle and assembly of the immature HIV-1 Gag lattice. IP6 makes ionic contacts with two rings of lysine residues at the centre of the Gag hexamer. Proteolytic cleavage then unmasks an alternative binding site, where IP6 interaction promotes the assembly of the mature capsid lattice. These studies identify IP6 as a naturally occurring small molecule that promotes both assembly and maturation of HIV-1.}, author = {Dick, Robert and Zadrozny, Kaneil K and Xu, Chaoyi and Schur, Florian and Lyddon, Terri D and Ricana, Clifton L and Wagner, Jonathan M and Perilla, Juan R and Ganser, Pornillos Barbie K and Johnson, Marc C and Pornillos, Owen and Vogt, Volker}, issn = {1476-4687}, journal = {Nature}, number = {7719}, pages = {509–512}, publisher = {Nature Publishing Group}, title = {{Inositol phosphates are assembly co-factors for HIV-1}}, doi = {10.1038/s41586-018-0396-4}, volume = {560}, year = {2018}, } @article{303, abstract = {The theory of tropical series, that we develop here, firstly appeared in the study of the growth of pluriharmonic functions. Motivated by waves in sandpile models we introduce a dynamic on the set of tropical series, and it is experimentally observed that this dynamic obeys a power law. So, this paper serves as a compilation of results we need for other articles and also introduces several objects interesting by themselves.}, author = {Kalinin, Nikita and Shkolnikov, Mikhail}, journal = {Discrete and Continuous Dynamical Systems- Series A}, number = {6}, pages = {2827 -- 2849}, publisher = {AIMS}, title = {{Introduction to tropical series and wave dynamic on them}}, doi = {10.3934/dcds.2018120}, volume = {38}, year = {2018}, } @article{282, abstract = {Adaptive introgression is common in nature and can be driven by selection acting on multiple, linked genes. We explore the effects of polygenic selection on introgression under the infinitesimal model with linkage. This model assumes that the introgressing block has an effectively infinite number of genes, each with an infinitesimal effect on the trait under selection. The block is assumed to introgress under directional selection within a native population that is genetically homogeneous. We use individual-based simulations and a branching process approximation to compute various statistics of the introgressing block, and explore how these depend on parameters such as the map length and initial trait value associated with the introgressing block, the genetic variability along the block, and the strength of selection. Our results show that the introgression dynamics of a block under infinitesimal selection is qualitatively different from the dynamics of neutral introgression. We also find that in the long run, surviving descendant blocks are likely to have intermediate lengths, and clarify how the length is shaped by the interplay between linkage and infinitesimal selection. Our results suggest that it may be difficult to distinguish introgression of single loci from that of genomic blocks with multiple, tightly linked and weakly selected loci.}, author = {Sachdeva, Himani and Barton, Nicholas H}, journal = {Genetics}, number = {4}, pages = {1279 -- 1303}, publisher = {Genetics Society of America}, title = {{Introgression of a block of genome under infinitesimal selection}}, doi = {10.1534/genetics.118.301018}, volume = {209}, year = {2018}, } @inproceedings{108, abstract = {Universal hashing found a lot of applications in computer science. In cryptography the most important fact about universal families is the so called Leftover Hash Lemma, proved by Impagliazzo, Levin and Luby. In the language of modern cryptography it states that almost universal families are good extractors. In this work we provide a somewhat surprising characterization in the opposite direction. Namely, every extractor with sufficiently good parameters yields a universal family on a noticeable fraction of its inputs. Our proof technique is based on tools from extremal graph theory applied to the \'collision graph\' induced by the extractor, and may be of independent interest. We discuss possible applications to the theory of randomness extractors and non-malleable codes.}, author = {Obremski, Marciej and Skorski, Maciej}, location = {Vail, CO, USA}, publisher = {IEEE}, title = {{Inverted leftover hash lemma}}, doi = {10.1109/ISIT.2018.8437654}, volume = {2018}, year = {2018}, } @inproceedings{160, abstract = {We present layered concurrent programs, a compact and expressive notation for specifying refinement proofs of concurrent programs. A layered concurrent program specifies a sequence of connected concurrent programs, from most concrete to most abstract, such that common parts of different programs are written exactly once. These programs are expressed in the ordinary syntax of imperative concurrent programs using gated atomic actions, sequencing, choice, and (recursive) procedure calls. Each concurrent program is automatically extracted from the layered program. We reduce refinement to the safety of a sequence of concurrent checker programs, one each to justify the connection between every two consecutive concurrent programs. These checker programs are also automatically extracted from the layered program. Layered concurrent programs have been implemented in the CIVL verifier which has been successfully used for the verification of several complex concurrent programs.}, author = {Kragl, Bernhard and Qadeer, Shaz}, location = {Oxford, UK}, pages = {79 -- 102}, publisher = {Springer}, title = {{Layered Concurrent Programs}}, doi = {10.1007/978-3-319-96145-3_5}, volume = {10981}, year = {2018}, } @article{82, abstract = {In experimental cultures, when bacteria are mixed with lytic (virulent) bacteriophage, bacterial cells resistant to the phage commonly emerge and become the dominant population of bacteria. Following the ascent of resistant mutants, the densities of bacteria in these simple communities become limited by resources rather than the phage. Despite the evolution of resistant hosts, upon which the phage cannot replicate, the lytic phage population is most commonly maintained in an apparently stable state with the resistant bacteria. Several mechanisms have been put forward to account for this result. Here we report the results of population dynamic/evolution experiments with a virulent mutant of phage Lambda, λVIR, and Escherichia coli in serial transfer cultures. We show that, following the ascent of λVIR-resistant bacteria, λVIRis maintained in the majority of cases in maltose-limited minimal media and in all cases in nutrient-rich broth. Using mathematical models and experiments, we show that the dominant mechanism responsible for maintenance of λVIRin these resource-limited populations dominated by resistant E. coli is a high rate of either phenotypic or genetic transition from resistance to susceptibility—a hitherto undemonstrated mechanism we term "leaky resistance." We discuss the implications of leaky resistance to our understanding of the conditions for the maintenance of phage in populations of bacteria—their “existence conditions.”.}, author = {Chaudhry, Waqas and Pleska, Maros and Shah, Nilang and Weiss, Howard and Mccall, Ingrid and Meyer, Justin and Gupta, Animesh and Guet, Calin C and Levin, Bruce}, journal = {PLoS Biology}, number = {8}, publisher = {Public Library of Science}, title = {{Leaky resistance and the conditions for the existence of lytic bacteriophage}}, doi = {10.1371/journal.pbio.2005971}, volume = {16}, year = {2018}, } @article{4, abstract = {We present a data-driven technique to instantly predict how fluid flows around various three-dimensional objects. Such simulation is useful for computational fabrication and engineering, but is usually computationally expensive since it requires solving the Navier-Stokes equation for many time steps. To accelerate the process, we propose a machine learning framework which predicts aerodynamic forces and velocity and pressure fields given a threedimensional shape input. Handling detailed free-form three-dimensional shapes in a data-driven framework is challenging because machine learning approaches usually require a consistent parametrization of input and output. We present a novel PolyCube maps-based parametrization that can be computed for three-dimensional shapes at interactive rates. This allows us to efficiently learn the nonlinear response of the flow using a Gaussian process regression. We demonstrate the effectiveness of our approach for the interactive design and optimization of a car body.}, author = {Umetani, Nobuyuki and Bickel, Bernd}, journal = {ACM Trans. Graph.}, number = {4}, publisher = {ACM}, title = {{Learning three-dimensional flow for interactive aerodynamic design}}, doi = {10.1145/3197517.3201325}, volume = {37}, year = {2018}, } @article{566, abstract = {We consider large random matrices X with centered, independent entries which have comparable but not necessarily identical variances. Girko's circular law asserts that the spectrum is supported in a disk and in case of identical variances, the limiting density is uniform. In this special case, the local circular law by Bourgade et. al. [11,12] shows that the empirical density converges even locally on scales slightly above the typical eigenvalue spacing. In the general case, the limiting density is typically inhomogeneous and it is obtained via solving a system of deterministic equations. Our main result is the local inhomogeneous circular law in the bulk spectrum on the optimal scale for a general variance profile of the entries of X. }, author = {Alt, Johannes and Erdös, László and Krüger, Torben H}, journal = {Annals Applied Probability }, number = {1}, pages = {148--203}, publisher = {Institute of Mathematical Statistics}, title = {{Local inhomogeneous circular law}}, doi = {10.1214/17-AAP1302}, volume = {28}, year = {2018}, } @article{106, abstract = {The goal of this article is to introduce the reader to the theory of intrinsic geometry of convex surfaces. We illustrate the power of the tools by proving a theorem on convex surfaces containing an arbitrarily long closed simple geodesic. Let us remind ourselves that a curve in a surface is called geodesic if every sufficiently short arc of the curve is length minimizing; if, in addition, it has no self-intersections, we call it simple geodesic. A tetrahedron with equal opposite edges is called isosceles. The axiomatic method of Alexandrov geometry allows us to work with the metrics of convex surfaces directly, without approximating it first by a smooth or polyhedral metric. Such approximations destroy the closed geodesics on the surface; therefore it is difficult (if at all possible) to apply approximations in the proof of our theorem. On the other hand, a proof in the smooth or polyhedral case usually admits a translation into Alexandrov’s language; such translation makes the result more general. In fact, our proof resembles a translation of the proof given by Protasov. Note that the main theorem implies in particular that a smooth convex surface does not have arbitrarily long simple closed geodesics. However we do not know a proof of this corollary that is essentially simpler than the one presented below.}, author = {Akopyan, Arseniy and Petrunin, Anton}, journal = {Mathematical Intelligencer}, number = {3}, pages = {26 -- 31}, publisher = {Springer}, title = {{Long geodesics on convex surfaces}}, doi = {10.1007/s00283-018-9795-5}, volume = {40}, year = {2018}, } @article{275, abstract = {Lymphatic endothelial cells (LECs) release extracellular chemokines to guide the migration of dendritic cells. In this study, we report that LECs also release basolateral exosome-rich endothelial vesicles (EEVs) that are secreted in greater numbers in the presence of inflammatory cytokines and accumulate in the perivascular stroma of small lymphatic vessels in human chronic inflammatory diseases. Proteomic analyses of EEV fractions identified > 1,700 cargo proteins and revealed a dominant motility-promoting protein signature. In vitro and ex vivo EEV fractions augmented cellular protrusion formation in a CX3CL1/fractalkine-dependent fashion and enhanced the directional migratory response of human dendritic cells along guidance cues. We conclude that perilymphatic LEC exosomes enhance exploratory behavior and thus promote directional migration of CX3CR1-expressing cells in complex tissue environments.}, author = {Brown, Markus and Johnson, Louise and Leone, Dario and Májek, Peter and Vaahtomeri, Kari and Senfter, Daniel and Bukosza, Nora and Schachner, Helga and Asfour, Gabriele and Langer, Brigitte and Hauschild, Robert and Parapatics, Katja and Hong, Young and Bennett, Keiryn and Kain, Renate and Detmar, Michael and Sixt, Michael K and Jackson, David and Kerjaschki, Dontscho}, journal = {Journal of Cell Biology}, number = {6}, pages = {2205 -- 2221}, publisher = {Rockefeller University Press}, title = {{Lymphatic exosomes promote dendritic cell migration along guidance cues}}, doi = {10.1083/jcb.201612051}, volume = {217}, year = {2018}, } @article{158, abstract = {The angiosperm seed is composed of three genetically distinct tissues: the diploid embryo that originates from the fertilized egg cell, the triploid endosperm that is produced from the fertilized central cell, and the maternal sporophytic integuments that develop into the seed coat1. At the onset of embryo development in Arabidopsis thaliana, the zygote divides asymmetrically, producing a small apical embryonic cell and a larger basal cell that connects the embryo to the maternal tissue2. The coordinated and synchronous development of the embryo and the surrounding integuments, and the alignment of their growth axes, suggest communication between maternal tissues and the embryo. In contrast to animals, however, where a network of maternal factors that direct embryo patterning have been identified3,4, only a few maternal mutations have been described to affect embryo development in plants5–7. Early embryo patterning in Arabidopsis requires accumulation of the phytohormone auxin in the apical cell by directed transport from the suspensor8–10. However, the origin of this auxin has remained obscure. Here we investigate the source of auxin for early embryogenesis and provide evidence that the mother plant coordinates seed development by supplying auxin to the early embryo from the integuments of the ovule. We show that auxin response increases in ovules after fertilization, due to upregulated auxin biosynthesis in the integuments, and this maternally produced auxin is required for correct embryo development.}, author = {Robert, Hélène and Park, Chulmin and Gutièrrez, Carla and Wójcikowska, Barbara and Pěnčík, Aleš and Novák, Ondřej and Chen, Junyi and Grunewald, Wim and Dresselhaus, Thomas and Friml, Jirí and Laux, Thomas}, journal = {Nature Plants}, number = {8}, pages = {548 -- 553}, publisher = {Nature Publishing Group}, title = {{Maternal auxin supply contributes to early embryo patterning in Arabidopsis}}, doi = {10.1038/s41477-018-0204-z}, volume = {4}, year = {2018}, } @article{152, abstract = {Complex I has an essential role in ATP production by coupling electron transfer from NADH to quinone with translocation of protons across the inner mitochondrial membrane. Isolated complex I deficiency is a frequent cause of mitochondrial inherited diseases. Complex I has also been implicated in cancer, ageing, and neurodegenerative conditions. Until recently, the understanding of complex I deficiency on the molecular level was limited due to the lack of high-resolution structures of the enzyme. However, due to developments in single particle cryo-electron microscopy (cryo-EM), recent studies have reported nearly atomic resolution maps and models of mitochondrial complex I. These structures significantly add to our understanding of complex I mechanism and assembly. The disease-causing mutations are discussed here in their structural context.}, author = {Fiedorczuk, Karol and Sazanov, Leonid A}, journal = {Trends in Cell Biology}, number = {10}, pages = {835 -- 867}, publisher = {Elsevier}, title = {{Mammalian mitochondrial complex I structure and disease causing mutations}}, doi = {10.1016/j.tcb.2018.06.006}, volume = {28}, year = {2018}, } @inproceedings{310, abstract = {A model of computation that is widely used in the formal analysis of reactive systems is symbolic algorithms. In this model the access to the input graph is restricted to consist of symbolic operations, which are expensive in comparison to the standard RAM operations. We give lower bounds on the number of symbolic operations for basic graph problems such as the computation of the strongly connected components and of the approximate diameter as well as for fundamental problems in model checking such as safety, liveness, and coliveness. Our lower bounds are linear in the number of vertices of the graph, even for constant-diameter graphs. For none of these problems lower bounds on the number of symbolic operations were known before. The lower bounds show an interesting separation of these problems from the reachability problem, which can be solved with O(D) symbolic operations, where D is the diameter of the graph. Additionally we present an approximation algorithm for the graph diameter which requires Õ(n/D) symbolic steps to achieve a (1 +ϵ)-approximation for any constant > 0. This compares to O(n/D) symbolic steps for the (naive) exact algorithm and O(D) symbolic steps for a 2-approximation. Finally we also give a refined analysis of the strongly connected components algorithms of [15], showing that it uses an optimal number of symbolic steps that is proportional to the sum of the diameters of the strongly connected components.}, author = {Chatterjee, Krishnendu and Dvorák, Wolfgang and Henzinger, Monika H and Loitzenbauer, Veronika}, location = {New Orleans, Louisiana, United States}, pages = {2341 -- 2356}, publisher = {ACM}, title = {{Lower bounds for symbolic computation on graphs: Strongly connected components, liveness, safety, and diameter}}, doi = {10.1137/1.9781611975031.151}, year = {2018}, } @article{436, abstract = {There has been significant interest recently in using complex quantum systems to create effective nonreciprocal dynamics. Proposals have been put forward for the realization of artificial magnetic fields for photons and phonons; experimental progress is fast making these proposals a reality. Much work has concentrated on the use of such systems for controlling the flow of signals, e.g., to create isolators or directional amplifiers for optical signals. In this Letter, we build on this work but move in a different direction. We develop the theory of and discuss a potential realization for the controllable flow of thermal noise in quantum systems. We demonstrate theoretically that the unidirectional flow of thermal noise is possible within quantum cascaded systems. Viewing an optomechanical platform as a cascaded system we show here that one can ultimately control the direction of the flow of thermal noise. By appropriately engineering the mechanical resonator, which acts as an artificial reservoir, the flow of thermal noise can be constrained to a desired direction, yielding a thermal rectifier. The proposed quantum thermal noise rectifier could potentially be used to develop devices such as a thermal modulator, a thermal router, and a thermal amplifier for nanoelectronic devices and superconducting circuits.}, author = {Barzanjeh, Shabir and Aquilina, Matteo and Xuereb, André}, journal = {Physical Review Letters}, number = {6}, publisher = {American Physical Society}, title = {{Manipulating the flow of thermal noise in quantum devices}}, doi = {10.1103/PhysRevLett.120.060601}, volume = {120}, year = {2018}, } @article{5858, abstract = {Spatial patterns are ubiquitous on the subcellular, cellular and tissue level, and can be studied using imaging techniques such as light and fluorescence microscopy. Imaging data provide quantitative information about biological systems; however, mechanisms causing spatial patterning often remain elusive. In recent years, spatio-temporal mathematical modelling has helped to overcome this problem. Yet, outliers and structured noise limit modelling of whole imaging data, and models often consider spatial summary statistics. Here, we introduce an integrated data-driven modelling approach that can cope with measurement artefacts and whole imaging data. Our approach combines mechanistic models of the biological processes with robust statistical models of the measurement process. The parameters of the integrated model are calibrated using a maximum-likelihood approach. We used this integrated modelling approach to study in vivo gradients of the chemokine (C-C motif) ligand 21 (CCL21). CCL21 gradients guide dendritic cells and are important in the adaptive immune response. Using artificial data, we verified that the integrated modelling approach provides reliable parameter estimates in the presence of measurement noise and that bias and variance of these estimates are reduced compared to conventional approaches. The application to experimental data allowed the parametrization and subsequent refinement of the model using additional mechanisms. Among other results, model-based hypothesis testing predicted lymphatic vessel-dependent concentration of heparan sulfate, the binding partner of CCL21. The selected model provided an accurate description of the experimental data and was partially validated using published data. Our findings demonstrate that integrated statistical modelling of whole imaging data is computationally feasible and can provide novel biological insights.}, author = {Hross, Sabrina and Theis, Fabian J. and Sixt, Michael K and Hasenauer, Jan}, issn = {17425689}, journal = {Journal of the Royal Society Interface}, number = {149}, publisher = {Royal Society Publishing}, title = {{Mechanistic description of spatial processes using integrative modelling of noise-corrupted imaging data}}, doi = {10.1098/rsif.2018.0600}, volume = {15}, year = {2018}, } @article{16, abstract = {We report quantitative evidence of mixing-layer elastic instability in a viscoelastic fluid flow between two widely spaced obstacles hindering a channel flow at Re 1 and Wi 1. Two mixing layers with nonuniform shear velocity profiles are formed in the region between the obstacles. The mixing-layer instability arises in the vicinity of an inflection point on the shear velocity profile with a steep variation in the elastic stress. The instability results in an intermittent appearance of small vortices in the mixing layers and an amplification of spatiotemporal averaged vorticity in the elastic turbulence regime. The latter is characterized through scaling of friction factor with Wi and both pressure and velocity spectra. Furthermore, the observations reported provide improved understanding of the stability of the mixing layer in a viscoelastic fluid at large elasticity, i.e., Wi 1 and Re 1 and oppose the current view of suppression of vorticity solely by polymer additives.}, author = {Varshney, Atul and Steinberg, Victor}, journal = {Physical Review Fluids}, number = {10}, publisher = {American Physical Society}, title = {{Mixing layer instability and vorticity amplification in a creeping viscoelastic flow}}, doi = {10.1103/PhysRevFluids.3.103303}, volume = {3}, year = {2018}, } @article{43, abstract = {The initial amount of pathogens required to start an infection within a susceptible host is called the infective dose and is known to vary to a large extent between different pathogen species. We investigate the hypothesis that the differences in infective doses are explained by the mode of action in the underlying mechanism of pathogenesis: Pathogens with locally acting mechanisms tend to have smaller infective doses than pathogens with distantly acting mechanisms. While empirical evidence tends to support the hypothesis, a formal theoretical explanation has been lacking. We give simple analytical models to gain insight into this phenomenon and also investigate a stochastic, spatially explicit, mechanistic within-host model for toxin-dependent bacterial infections. The model shows that pathogens secreting locally acting toxins have smaller infective doses than pathogens secreting diffusive toxins, as hypothesized. While local pathogenetic mechanisms require smaller infective doses, pathogens with distantly acting toxins tend to spread faster and may cause more damage to the host. The proposed model can serve as a basis for the spatially explicit analysis of various virulence factors also in the context of other problems in infection dynamics.}, author = {Rybicki, Joel and Kisdi, Eva and Anttila, Jani}, journal = {PNAS}, number = {42}, pages = {10690 -- 10695}, publisher = {National Academy of Sciences}, title = {{Model of bacterial toxin-dependent pathogenesis explains infective dose}}, doi = {10.1073/pnas.1721061115}, volume = {115}, year = {2018}, } @article{13, abstract = {We propose a new method for fabricating digital objects through reusable silicone molds. Molds are generated by casting liquid silicone into custom 3D printed containers called metamolds. Metamolds automatically define the cuts that are needed to extract the cast object from the silicone mold. The shape of metamolds is designed through a novel segmentation technique, which takes into account both geometric and topological constraints involved in the process of mold casting. Our technique is simple, does not require changing the shape or topology of the input objects, and only requires off-the- shelf materials and technologies. We successfully tested our method on a set of challenging examples with complex shapes and rich geometric detail. © 2018 Association for Computing Machinery.}, author = {Alderighi, Thomas and Malomo, Luigi and Giorgi, Daniela and Pietroni, Nico and Bickel, Bernd and Cignoni, Paolo}, journal = {ACM Trans. Graph.}, number = {4}, publisher = {ACM}, title = {{Metamolds: Computational design of silicone molds}}, doi = {10.1145/3197517.3201381}, volume = {37}, year = {2018}, } @article{137, abstract = {Fluorescent sensors are an essential part of the experimental toolbox of the life sciences, where they are used ubiquitously to visualize intra- and extracellular signaling. In the brain, optical neurotransmitter sensors can shed light on temporal and spatial aspects of signal transmission by directly observing, for instance, neurotransmitter release and spread. Here we report the development and application of the first optical sensor for the amino acid glycine, which is both an inhibitory neurotransmitter and a co-agonist of the N-methyl-d-aspartate receptors (NMDARs) involved in synaptic plasticity. Computational design of a glycine-specific binding protein allowed us to produce the optical glycine FRET sensor (GlyFS), which can be used with single and two-photon excitation fluorescence microscopy. We took advantage of this newly developed sensor to test predictions about the uneven spatial distribution of glycine in extracellular space and to demonstrate that extracellular glycine levels are controlled by plasticity-inducing stimuli.}, author = {Zhang, William and Herde, Michel and Mitchell, Joshua and Whitfield, Jason and Wulff, Andreas and Vongsouthi, Vanessa and Sanchez Romero, Inmaculada and Gulakova, Polina and Minge, Daniel and Breithausen, Björn and Schoch, Susanne and Janovjak, Harald L and Jackson, Colin and Henneberger, Christian}, journal = {Nature Chemical Biology}, number = {9}, pages = {861 -- 869}, publisher = {Nature Publishing Group}, title = {{Monitoring hippocampal glycine with the computationally designed optical sensor GlyFS}}, doi = {10.1038/s41589-018-0108-2}, volume = {14}, year = {2018}, } @article{276, abstract = {Directed migration of cells relies on their ability to sense directional guidance cues and to interact with pericellular structures in order to transduce contractile cytoskeletal- into mechanical forces. These biomechanical processes depend highly on microenvironmental factors such as exposure to 2D surfaces or 3D matrices. In vivo, the majority of cells are exposed to 3D environments. Data on 3D cell migration are mostly derived from intravital microscopy or collagen-based in vitro assays. Both approaches offer only limited controlla-bility of experimental conditions. Here, we developed an automated microfluidic system that allows positioning of cells in 3D microenvironments containing highly controlled diffusion-based chemokine gradients. Tracking migration in such gradients was feasible in real time at the single cell level. Moreover, the setup allowed on-chip immunocytochemistry and thus linking of functional with phenotypical properties in individual cells. Spatially defined retrieval of cells from the device allows down-stream off-chip analysis. Using dendritic cells as a model, our setup specifically allowed us for the first time to quantitate key migration characteristics of cells exposed to identical gradients of the chemokine CCL19 yet placed on 2D vs in 3D environments. Migration properties between 2D and 3D migration were distinct. Morphological features of cells migrating in an in vitro 3D environment were similar to those of cells migrating in animal tissues, but different from cells migrating on a surface. Our system thus offers a highly controllable in vitro-mimic of a 3D environment that cells traffic in vivo.}, author = {Frick, Corina and Dettinger, Philip and Renkawitz, Jörg and Jauch, Annaïse and Berger, Christoph and Recher, Mike and Schroeder, Timm and Mehling, Matthias}, journal = {PLoS One}, number = {6}, publisher = {Public Library of Science}, title = {{Nano-scale microfluidics to study 3D chemotaxis at the single cell level}}, doi = {10.1371/journal.pone.0198330}, volume = {13}, year = {2018}, } @article{283, abstract = {Light represents the principal signal driving circadian clock entrainment. However, how light influences the evolution of the clock remains poorly understood. The cavefish Phreatichthys andruzzii represents a fascinating model to explore how evolution under extreme aphotic conditions shapes the circadian clock, since in this species the clock is unresponsive to light. We have previously demonstrated that loss-of-function mutations targeting non-visual opsins contribute in part to this blind clock phenotype. Here, we have compared orthologs of two core clock genes that play a key role in photic entrainment, cry1a and per2, in both zebrafish and P. andruzzii. We encountered aberrantly spliced variants for the P. andruzzii per2 transcript. The most abundant transcript encodes a truncated protein lacking the C-terminal Cry binding domain and incorporating an intronic, transposon-derived coding sequence. We demonstrate that the transposon insertion leads to a predominantly cytoplasmic localization of the cavefish Per2 protein in contrast to the zebrafish ortholog which is distributed in both the nucleus and cytoplasm. Thus, it seems that during evolution in complete darkness, the photic entrainment pathway of the circadian clock has been subject to mutation at multiple levels, extending from opsin photoreceptors to nuclear effectors.}, author = {Ceinos, Rosa Maria and Frigato, Elena and Pagano, Cristina and Frohlich, Nadine and Negrini, Pietro and Cavallari, Nicola and Vallone, Daniela and Fuselli, Silvia and Bertolucci, Cristiano and Foulkes, Nicholas S}, journal = {Scientific Reports}, number = {1}, publisher = {Nature Publishing Group}, title = {{Mutations in blind cavefish target the light regulated circadian clock gene period 2}}, doi = {10.1038/s41598-018-27080-2}, volume = {8}, year = {2018}, } @inproceedings{81, abstract = {We solve the offline monitoring problem for timed propositional temporal logic (TPTL), interpreted over dense-time Boolean signals. The variant of TPTL we consider extends linear temporal logic (LTL) with clock variables and reset quantifiers, providing a mechanism to specify real-time constraints. We first describe a general monitoring algorithm based on an exhaustive computation of the set of satisfying clock assignments as a finite union of zones. We then propose a specialized monitoring algorithm for the one-variable case using a partition of the time domain based on the notion of region equivalence, whose complexity is linear in the length of the signal, thereby generalizing a known result regarding the monitoring of metric temporal logic (MTL). The region and zone representations of time constraints are known from timed automata verification and can also be used in the discrete-time case. Our prototype implementation appears to outperform previous discrete-time implementations of TPTL monitoring,}, author = {Elgyütt, Adrian and Ferrere, Thomas and Henzinger, Thomas A}, location = {Beijing, China}, pages = {53 -- 70}, publisher = {Springer}, title = {{Monitoring temporal logic with clock variables}}, doi = {10.1007/978-3-030-00151-3_4}, volume = {11022}, year = {2018}, } @article{76, abstract = {Consider a fully-connected synchronous distributed system consisting of n nodes, where up to f nodes may be faulty and every node starts in an arbitrary initial state. In the synchronous C-counting problem, all nodes need to eventually agree on a counter that is increased by one modulo C in each round for given C>1. In the self-stabilising firing squad problem, the task is to eventually guarantee that all non-faulty nodes have simultaneous responses to external inputs: if a subset of the correct nodes receive an external “go” signal as input, then all correct nodes should agree on a round (in the not-too-distant future) in which to jointly output a “fire” signal. Moreover, no node should generate a “fire” signal without some correct node having previously received a “go” signal as input. We present a framework reducing both tasks to binary consensus at very small cost. For example, we obtain a deterministic algorithm for self-stabilising Byzantine firing squads with optimal resilience f<n/3, asymptotically optimal stabilisation and response time O(f), and message size O(log f). As our framework does not restrict the type of consensus routines used, we also obtain efficient randomised solutions.}, author = {Lenzen, Christoph and Rybicki, Joel}, journal = {Distributed Computing}, publisher = {Springer}, title = {{Near-optimal self-stabilising counting and firing squads}}, doi = {10.1007/s00446-018-0342-6}, year = {2018}, } @article{530, abstract = {Inclusion–exclusion is an effective method for computing the volume of a union of measurable sets. We extend it to multiple coverings, proving short inclusion–exclusion formulas for the subset of Rn covered by at least k balls in a finite set. We implement two of the formulas in dimension n=3 and report on results obtained with our software.}, author = {Edelsbrunner, Herbert and Iglesias Ham, Mabel}, journal = {Computational Geometry: Theory and Applications}, pages = {119 -- 133}, publisher = {Elsevier}, title = {{Multiple covers with balls I: Inclusion–exclusion}}, doi = {10.1016/j.comgeo.2017.06.014}, volume = {68}, year = {2018}, } @article{307, abstract = {Spontaneous emission spectra of two initially excited closely spaced identical atoms are very sensitive to the strength and the direction of the applied magnetic field. We consider the relevant schemes that ensure the determination of the mutual spatial orientation of the atoms and the distance between them by entirely optical means. A corresponding theoretical description is given accounting for the dipole-dipole interaction between the two atoms in the presence of a magnetic field and for polarizations of the quantum field interacting with magnetic sublevels of the two-atom system. }, author = {Redchenko, Elena and Makarov, Alexander and Yudson, Vladimir}, journal = { Physical Review A - Atomic, Molecular, and Optical Physics}, number = {4}, publisher = {American Physical Society}, title = {{Nanoscopy of pairs of atoms by fluorescence in a magnetic field}}, doi = {10.1103/PhysRevA.97.043812}, volume = {97}, year = {2018}, } @article{279, abstract = {Background: Natural selection shapes cancer genomes. Previous studies used signatures of positive selection to identify genes driving malignant transformation. However, the contribution of negative selection against somatic mutations that affect essential tumor functions or specific domains remains a controversial topic. Results: Here, we analyze 7546 individual exomes from 26 tumor types from TCGA data to explore the portion of the cancer exome under negative selection. Although we find most of the genes neutrally evolving in a pan-cancer framework, we identify essential cancer genes and immune-exposed protein regions under significant negative selection. Moreover, our simulations suggest that the amount of negative selection is underestimated. We therefore choose an empirical approach to identify genes, functions, and protein regions under negative selection. We find that expression and mutation status of negatively selected genes is indicative of patient survival. Processes that are most strongly conserved are those that play fundamental cellular roles such as protein synthesis, glucose metabolism, and molecular transport. Intriguingly, we observe strong signals of selection in the immunopeptidome and proteins controlling peptide exposition, highlighting the importance of immune surveillance evasion. Additionally, tumor type-specific immune activity correlates with the strength of negative selection on human epitopes. Conclusions: In summary, our results show that negative selection is a hallmark of cell essentiality and immune response in cancer. The functional domains identified could be exploited therapeutically, ultimately allowing for the development of novel cancer treatments.}, author = {Zapata, Luis and Pich, Oriol and Serrano, Luis and Kondrashov, Fyodor and Ossowski, Stephan and Schaefer, Martin}, journal = {Genome Biology}, publisher = {BioMed Central}, title = {{Negative selection in tumor genome evolution acts on essential cellular functions and the immunopeptidome}}, doi = {10.1186/s13059-018-1434-0}, volume = {19}, year = {2018}, }