@inproceedings{397, abstract = {Concurrent sets with range query operations are highly desirable in applications such as in-memory databases. However, few set implementations offer range queries. Known techniques for augmenting data structures with range queries (or operations that can be used to build range queries) have numerous problems that limit their usefulness. For example, they impose high overhead or rely heavily on garbage collection. In this work, we show how to augment data structures with highly efficient range queries, without relying on garbage collection. We identify a property of epoch-based memory reclamation algorithms that makes them ideal for implementing range queries, and produce three algorithms, which use locks, transactional memory and lock-free techniques, respectively. Our algorithms are applicable to more data structures than previous work, and are shown to be highly efficient on a large scale Intel system. }, author = {Arbel Raviv, Maya and Brown, Trevor A}, isbn = {978-1-4503-4982-6}, location = {Vienna, Austria}, number = {1}, pages = {14 -- 27}, publisher = {ACM}, title = {{Harnessing epoch-based reclamation for efficient range queries}}, doi = {10.1145/3178487.3178489}, volume = {53}, 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{398, abstract = {Objective: To report long-term results after Pipeline Embolization Device (PED) implantation, characterize complex and standard aneurysms comprehensively, and introduce a modified flow disruption scale. Methods: We retrospectively reviewed a consecutive series of 40 patients harboring 59 aneurysms treated with 54 PEDs. Aneurysm complexity was assessed using our proposed classification. Immediate angiographic results were analyzed using previously published grading scales and our novel flow disruption scale. Results: According to our new definition, 46 (78%) aneurysms were classified as complex. Most PED interventions were performed in the paraophthalmic and cavernous internal carotid artery segments. Excellent neurologic outcome (modified Rankin Scale 0 and 1) was observed in 94% of patients. Our data showed low permanent procedure-related mortality (0%) and morbidity (3%) rates. Long-term angiographic follow-up showed complete occlusion in 81% and near-total obliteration in a further 14%. Complete obliteration after deployment of a single PED was achieved in all standard aneurysms with 1-year follow-up. Our new scale was an independent predictor of aneurysm occlusion in a multivariable analysis. All aneurysms with a high flow disruption grade showed complete occlusion at follow-up regardless of PED number or aneurysm complexity. Conclusions: Treatment with the PED should be recognized as a primary management strategy for a highly selected cohort with predominantly complex intracranial aneurysms. We further show that a priori assessment of aneurysm complexity and our new postinterventional angiographic flow disruption scale predict occlusion probability and may help to determine the adequate number of per-aneurysm devices.}, author = {Dodier, Philippe and Frischer, Josa and Wang, Wei and Auzinger, Thomas and Mallouhi, Ammar and Serles, Wolfgang and Gruber, Andreas and Knosp, Engelbert and Bavinzski, Gerhard}, journal = {World Neurosurgery}, pages = {e568--e578}, publisher = {Elsevier}, title = {{Immediate flow disruption as a prognostic factor after flow diverter treatment long term experience with the pipeline embolization device}}, doi = {10.1016/j.wneu.2018.02.096}, volume = {13}, 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}, } @article{426, abstract = {Sperm cells are the most morphologically diverse cells across animal taxa. Within species, sperm and ejaculate traits have been suggested to vary with the male's competitive environment, e.g., level of sperm competition, female mating status and quality, and also with male age, body mass, physiological condition, and resource availability. Most previous studies have based their conclusions on the analysis of only one or a few ejaculates per male without investigating differences among the ejaculates of the same individual. This masks potential ejaculate-specific traits. Here, we provide data on the length, quantity, and viability of sperm ejaculated by wingless males of the ant Cardiocondyla obscurior. Males of this ant species are relatively long-lived and can mate with large numbers of female sexuals throughout their lives. We analyzed all ejaculates across the individuals' lifespan and manipulated the availability of mating partners. Our study shows that both the number and size of sperm cells transferred during copulations differ among individuals and also among ejaculates of the same male. Sperm quality does not decrease with male age, but the variation in sperm number between ejaculates indicates that males need considerable time to replenish their sperm supplies. Producing many ejaculates in a short time appears to be traded-off against male longevity rather than sperm quality.}, author = {Metzler, Sina and Schrempf, Alexandra and Heinze, Jürgen}, journal = {Journal of Insect Physiology}, pages = {284--290}, publisher = {Elsevier}, title = {{Individual- and ejaculate-specific sperm traits in ant males}}, doi = {10.1016/j.jinsphys.2017.12.003}, volume = {107}, 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{280, abstract = {Flowers have a species-specific functional life span that determines the time window in which pollination, fertilization and seed set can occur. The stigma tissue plays a key role in flower receptivity by intercepting pollen and initiating pollen tube growth toward the ovary. In this article, we show that a developmentally controlled cell death programme terminates the functional life span of stigma cells in Arabidopsis. We identified the leaf senescence regulator ORESARA1 (also known as ANAC092) and the previously uncharacterized KIRA1 (also known as ANAC074) as partially redundant transcription factors that modulate stigma longevity by controlling the expression of programmed cell death-associated genes. KIRA1 expression is sufficient to induce cell death and terminate floral receptivity, whereas lack of both KIRA1 and ORESARA1 substantially increases stigma life span. Surprisingly, the extension of stigma longevity is accompanied by only a moderate extension of flower receptivity, suggesting that additional processes participate in the control of the flower's receptive life span.}, author = {Gao, Zhen and Daneva, Anna and Salanenka, Yuliya and Van Durme, Matthias and Huysmans, Marlies and Lin, Zongcheng and De Winter, Freya and Vanneste, Steffen and Karimi, Mansour and Van De Velde, Jan and Vandepoele, Klaas and Van De Walle, Davy and Dewettinck, Koen and Lambrecht, Bart and Nowack, Moritz}, journal = {Nature Plants}, number = {6}, pages = {365 -- 375}, publisher = {Nature Publishing Group}, title = {{KIRA1 and ORESARA1 terminate flower receptivity by promoting cell death in the stigma of Arabidopsis}}, doi = {10.1038/s41477-018-0160-7}, volume = {4}, year = {2018}, } @article{503, abstract = {Buffers are essential for diluting bacterial cultures for flow cytometry analysis in order to study bacterial physiology and gene expression parameters based on fluorescence signals. Using a variety of constitutively expressed fluorescent proteins in Escherichia coli K-12 strain MG1655, we found strong artifactual changes in fluorescence levels after dilution into the commonly used flow cytometry buffer phosphate-buffered saline (PBS) and two other buffer solutions, Tris-HCl and M9 salts. These changes appeared very rapidly after dilution, and were linked to increased membrane permeability and loss in cell viability. We observed buffer-related effects in several different E. coli strains, K-12, C and W, but not E. coli B, which can be partially explained by differences in lipopolysaccharide (LPS) and outer membrane composition. Supplementing the buffers with divalent cations responsible for outer membrane stability, Mg2+ and Ca2+, preserved fluorescence signals, membrane integrity and viability of E. coli. Thus, stabilizing the bacterial outer membrane is essential for precise and unbiased measurements of fluorescence parameters using flow cytometry.}, author = {Tomasek, Kathrin and Bergmiller, Tobias and Guet, Calin C}, journal = {Journal of Biotechnology}, pages = {40 -- 52}, publisher = {Elsevier}, title = {{Lack of cations in flow cytometry buffers affect fluorescence signals by reducing membrane stability and viability of Escherichia coli strains}}, doi = {10.1016/j.jbiotec.2018.01.008}, volume = {268}, 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}, } @inproceedings{183, abstract = {Fault-localization is considered to be a very tedious and time-consuming activity in the design of complex Cyber-Physical Systems (CPS). This laborious task essentially requires expert knowledge of the system in order to discover the cause of the fault. In this context, we propose a new procedure that AIDS designers in debugging Simulink/Stateflow hybrid system models, guided by Signal Temporal Logic (STL) specifications. The proposed method relies on three main ingredients: (1) a monitoring and a trace diagnostics procedure that checks whether a tested behavior satisfies or violates an STL specification, localizes time segments and interfaces variables contributing to the property violations; (2) a slicing procedure that maps these observable behavior segments to the internal states and transitions of the Simulink model; and (3) a spectrum-based fault-localization method that combines the previous analysis from multiple tests to identify the internal states and/or transitions that are the most likely to explain the fault. We demonstrate the applicability of our approach on two Simulink models from the automotive and the avionics domain.}, author = {Bartocci, Ezio and Ferrere, Thomas and Manjunath, Niveditha and Nickovic, Dejan}, location = {Porto, Portugal}, pages = {197 -- 206}, publisher = {Association for Computing Machinery, Inc}, title = {{Localizing faults in simulink/stateflow models with STL}}, doi = {10.1145/3178126.3178131}, 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}, } @misc{9810, 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}, publisher = {Public Library of Science}, title = {{Numerical data used in figures}}, doi = {10.1371/journal.pbio.2005971.s008}, year = {2018}, }