@article{1072, abstract = {Given a finite set of points in Rn and a radius parameter, we study the Čech, Delaunay–Čech, Delaunay (or alpha), and Wrap complexes in the light of generalized discrete Morse theory. Establishing the Čech and Delaunay complexes as sublevel sets of generalized discrete Morse functions, we prove that the four complexes are simple-homotopy equivalent by a sequence of simplicial collapses, which are explicitly described by a single discrete gradient field.}, author = {Bauer, Ulrich and Edelsbrunner, Herbert}, journal = {Transactions of the American Mathematical Society}, number = {5}, pages = {3741 -- 3762}, publisher = {American Mathematical Society}, title = {{The Morse theory of Čech and delaunay complexes}}, doi = {10.1090/tran/6991}, volume = {369}, year = {2017}, } @article{1073, abstract = {Let X and Y be finite simplicial sets (e.g. finite simplicial complexes), both equipped with a free simplicial action of a finite group G. Assuming that Y is d-connected and dimX≤2d, for some d≥1, we provide an algorithm that computes the set of all equivariant homotopy classes of equivariant continuous maps |X|→|Y|; the existence of such a map can be decided even for dimX≤2d+1. This yields the first algorithm for deciding topological embeddability of a k-dimensional finite simplicial complex into Rn under the condition k≤23n−1. More generally, we present an algorithm that, given a lifting-extension problem satisfying an appropriate stability assumption, computes the set of all homotopy classes of solutions. This result is new even in the non-equivariant situation.}, author = {Čadek, Martin and Krcál, Marek and Vokřínek, Lukáš}, issn = {01795376}, journal = {Discrete & Computational Geometry}, number = {4}, pages = {915 -- 965}, publisher = {Springer}, title = {{Algorithmic solvability of the lifting extension problem}}, doi = {10.1007/s00454-016-9855-6}, volume = {54}, year = {2017}, } @article{1074, abstract = {Recently it has become feasible to detect long blocks of nearly identical sequence shared between pairs of genomes. These IBD blocks are direct traces of recent coalescence events and, as such, contain ample signal to infer recent demography. Here, we examine sharing of such blocks in two-dimensional populations with local migration. Using a diffusion approximation to trace genetic ancestry, we derive analytical formulae for patterns of isolation by distance of IBD blocks, which can also incorporate recent population density changes. We introduce an inference scheme that uses a composite likelihood approach to fit these formulae. We then extensively evaluate our theory and inference method on a range of scenarios using simulated data. We first validate the diffusion approximation by showing that the theoretical results closely match the simulated block sharing patterns. We then demonstrate that our inference scheme can accurately and robustly infer dispersal rate and effective density, as well as bounds on recent dynamics of population density. To demonstrate an application, we use our estimation scheme to explore the fit of a diffusion model to Eastern European samples in the POPRES data set. We show that ancestry diffusing with a rate of σ ≈ 50–100 km/√gen during the last centuries, combined with accelerating population growth, can explain the observed exponential decay of block sharing with increasing pairwise sample distance.}, author = {Ringbauer, Harald and Coop, Graham and Barton, Nicholas H}, issn = {00166731}, journal = {Genetics}, number = {3}, pages = {1335 -- 1351}, publisher = {Genetics Society of America}, title = {{Inferring recent demography from isolation by distance of long shared sequence blocks}}, doi = {10.1534/genetics.116.196220}, volume = {205}, year = {2017}, } @inproceedings{10745, abstract = {New ways to investigate and manipulate fluxoid and vortex states of mesoscopic superconducting structures are of great interest. The states with multiple vortices or winding numbers could be useful for the study of vortex interactions and interference effects, the braiding of Majorana bound states by winding vortices, and the development of novel superconducting devices. We demonstrate a methodology based on magnetic force microscopy that allows us to induce, probe and control fluxoid states in thin wall structures comprised of multiple loops. By using micro-magnet as a source of inhomogeneous magnetic field, we can efficiently explore the configuration space of fluxoid states. Scanning over the structure reveals the energy crossing points of the lowest laying fluxoid states. This is due the strong interaction of cantilever with thermally activated fluxoid transitions at points of degeneracy. We show that measured patterns of fluxoid transitions allow to identify the states, investigate their energetics, and manipulate them. Further, we show that the dynamics of driven fluxoid transitions can be described by stochastic resonance model, which provides a unique way of measuring fluxoid transition rate and related energy barrier for chosen transitions even in complicated structures}, author = {Polshyn, Hryhoriy and Naibert, Tyler and Budakian, Raffi}, booktitle = {APS March Meeting 2017}, issn = {0003-0503}, location = {New Orleans, LA, United States}, number = {4}, publisher = {American Physical Society}, title = {{ Probing and controlling fluxoid states in multiply-connected mesoscopic superconducting structures}}, volume = {62}, year = {2017}, } @inbook{1075, author = {Wenzl, Bernhard}, booktitle = {Austria and America: 20th-Century Cross-Cultural Encounters}, editor = {Parker, Joshua and Poole, Ralph}, isbn = {978-3643908124}, pages = {73 -- 80}, publisher = {LIT Verlag Berlin-Münster-Wien-Zürich-London}, title = {{An American in Allied-occupied Austria: John Dos Passos Reports on "The Vienna Frontier"}}, volume = {15}, year = {2017}, } @article{1076, abstract = {Signatures of the Coulomb corrections in the photoelectron momentum distribution during laser-induced ionization of atoms or ions in tunneling and multiphoton regimes are investigated analytically in the case of a one-dimensional problem. A high-order Coulomb-corrected strong-field approximation is applied, where the exact continuum state in the S matrix is approximated by the eikonal Coulomb-Volkov state including the second-order corrections to the eikonal. Although without high-order corrections our theory coincides with the known analytical R-matrix (ARM) theory, we propose a simplified procedure for the matrix element derivation. Rather than matching the eikonal Coulomb-Volkov wave function with the bound state as in the ARM theory to remove the Coulomb singularity, we calculate the matrix element via the saddle-point integration method by time as well as by coordinate, and in this way avoiding the Coulomb singularity. The momentum shift in the photoelectron momentum distribution with respect to the ARM theory due to high-order corrections is analyzed for tunneling and multiphoton regimes. The relation of the quantum corrections to the tunneling delay time is discussed.}, author = {Klaiber, Michael and Daněk, Jiří and Yakaboylu, Enderalp and Hatsagortsyan, Karen and Keitel, Christoph}, issn = {24699926}, journal = { Physical Review A - Atomic, Molecular, and Optical Physics}, number = {2}, publisher = {American Physical Society}, title = {{Strong-field ionization via a high-order Coulomb-corrected strong-field approximation}}, doi = {10.1103/PhysRevA.95.023403}, volume = {95}, year = {2017}, } @article{1077, abstract = {Viral capsids are structurally constrained by interactions among the amino acids (AAs) of their constituent proteins. Therefore, epistasis is expected to evolve among physically interacting sites and to influence the rates of substitution. To study the evolution of epistasis, we focused on the major structural protein of the fX174 phage family by first reconstructing the ancestral protein sequences of 18 species using a Bayesian statistical framework. The inferred ancestral reconstruction differed at eight AAs, for a total of 256 possible ancestral haplotypes. For each ancestral haplotype and the extant species, we estimated, in silico, the distribution of free energies and epistasis of the capsid structure. We found that free energy has not significantly increased but epistasis has. We decomposed epistasis up to fifth order and found that higher-order epistasis sometimes compensates pairwise interactions making the free energy seem additive. The dN/dS ratio is low, suggesting strong purifying selection, and that structure is under stabilizing selection. We synthesized phages carrying ancestral haplotypes of the coat protein gene and measured their fitness experimentally. Our findings indicate that stabilizing mutations can have higher fitness, and that fitness optima do not necessarily coincide with energy minima.}, author = {Fernandes Redondo, Rodrigo A and Vladar, Harold and Włodarski, Tomasz and Bollback, Jonathan P}, issn = {17425689}, journal = {Journal of the Royal Society Interface}, number = {126}, publisher = {Royal Society of London}, title = {{Evolutionary interplay between structure, energy and epistasis in the coat protein of the ϕX174 phage family}}, doi = {10.1098/rsif.2016.0139}, volume = {14}, year = {2017}, } @article{1078, abstract = {One of the key questions in understanding plant development is how single cells behave in a larger context of the tissue. Therefore, it requires the observation of the whole organ with a high spatial- as well as temporal resolution over prolonged periods of time, which may cause photo-toxic effects. This protocol shows a plant sample preparation method for light-sheet microscopy, which is characterized by mounting the plant vertically on the surface of a gel. The plant is mounted in such a way that the roots are submerged in a liquid medium while the leaves remain in the air. In order to ensure photosynthetic activity of the plant, a custom-made lighting system illuminates the leaves. To keep the roots in darkness the water surface is covered with sheets of black plastic foil. This method allows long-term imaging of plant organ development in standardized conditions. }, author = {Von Wangenheim, Daniel and Hauschild, Robert and Friml, Jirí}, journal = {Journal of visualized experiments JoVE}, number = {119}, publisher = {Journal of Visualized Experiments}, title = {{Light sheet fluorescence microscopy of plant roots growing on the surface of a gel}}, doi = {10.3791/55044}, volume = {2017}, year = {2017}, } @article{1079, abstract = {We study the ionization problem in the Thomas-Fermi-Dirac-von Weizsäcker theory for atoms and molecules. We prove the nonexistence of minimizers for the energy functional when the number of electrons is large and the total nuclear charge is small. This nonexistence result also applies to external potentials decaying faster than the Coulomb potential. In the case of arbitrary nuclear charges, we obtain the nonexistence of stable minimizers and radial minimizers.}, author = {Nam, Phan and Van Den Bosch, Hanne}, issn = {13850172}, journal = {Mathematical Physics, Analysis and Geometry}, number = {2}, publisher = {Springer}, title = {{Nonexistence in Thomas Fermi-Dirac-von Weizsäcker theory with small nuclear charges}}, doi = {10.1007/s11040-017-9238-0}, volume = {20}, year = {2017}, } @article{1080, abstract = {Reconstructing the evolutionary history of metastases is critical for understanding their basic biological principles and has profound clinical implications. Genome-wide sequencing data has enabled modern phylogenomic methods to accurately dissect subclones and their phylogenies from noisy and impure bulk tumour samples at unprecedented depth. However, existing methods are not designed to infer metastatic seeding patterns. Here we develop a tool, called Treeomics, to reconstruct the phylogeny of metastases and map subclones to their anatomic locations. Treeomics infers comprehensive seeding patterns for pancreatic, ovarian, and prostate cancers. Moreover, Treeomics correctly disambiguates true seeding patterns from sequencing artifacts; 7% of variants were misclassified by conventional statistical methods. These artifacts can skew phylogenies by creating illusory tumour heterogeneity among distinct samples. In silico benchmarking on simulated tumour phylogenies across a wide range of sample purities (15–95%) and sequencing depths (25-800 × ) demonstrates the accuracy of Treeomics compared with existing methods.}, author = {Reiter, Johannes and Makohon Moore, Alvin and Gerold, Jeffrey and Božić, Ivana and Chatterjee, Krishnendu and Iacobuzio Donahue, Christine and Vogelstein, Bert and Nowak, Martin}, issn = {20411723}, journal = {Nature Communications}, publisher = {Nature Publishing Group}, title = {{Reconstructing metastatic seeding patterns of human cancers}}, doi = {10.1038/ncomms14114}, volume = {8}, year = {2017}, } @article{943, abstract = {Like many developing tissues, the vertebrate neural tube is patterned by antiparallel morphogen gradients. To understand how these inputs are interpreted, we measured morphogen signaling and target gene expression in mouse embryos and chick ex vivo assays. From these data, we derived and validated a characteristic decoding map that relates morphogen input to the positional identity of neural progenitors. Analysis of the observed responses indicates that the underlying interpretation strategy minimizes patterning errors in response to the joint input of noisy opposing gradients. We reverse-engineered a transcriptional network that provides a mechanistic basis for the observed cell fate decisions and accounts for the precision and dynamics of pattern formation. Together, our data link opposing gradient dynamics in a growing tissue to precise pattern formation.}, author = {Zagórski, Marcin P and Tabata, Yoji and Brandenberg, Nathalie and Lutolf, Matthias and Tkacik, Gasper and Bollenbach, Tobias and Briscoe, James and Kicheva, Anna}, issn = {00368075}, journal = {Science}, number = {6345}, pages = {1379 -- 1383}, publisher = {American Association for the Advancement of Science}, title = {{Decoding of position in the developing neural tube from antiparallel morphogen gradients}}, doi = {10.1126/science.aam5887}, volume = {356}, year = {2017}, } @article{9445, abstract = {Cytosine methylation regulates essential genome functions across eukaryotes, but the fundamental question of whether nucleosomal or naked DNA is the preferred substrate of plant and animal methyltransferases remains unresolved. Here, we show that genetic inactivation of a single DDM1/Lsh family nucleosome remodeler biases methylation toward inter-nucleosomal linker DNA in Arabidopsis thaliana and mouse. We find that DDM1 enables methylation of DNA bound to the nucleosome, suggesting that nucleosome-free DNA is the preferred substrate of eukaryotic methyltransferases in vivo. Furthermore, we show that simultaneous mutation of DDM1 and linker histone H1 in Arabidopsis reproduces the strong linker-specific methylation patterns of species that diverged from flowering plants and animals over a billion years ago. Our results indicate that in the absence of remodeling, nucleosomes are strong barriers to DNA methyltransferases. Linker-specific methylation can evolve simply by breaking the connection between nucleosome remodeling and DNA methylation.}, author = {Lyons, David B and Zilberman, Daniel}, issn = {2050-084X}, journal = {eLife}, publisher = {eLife Sciences Publications}, title = {{DDM1 and Lsh remodelers allow methylation of DNA wrapped in nucleosomes}}, doi = {10.7554/elife.30674}, volume = {6}, year = {2017}, } @article{945, abstract = {While chromosome-wide dosage compensation of the X chromosome has been found in many species, studies in ZW clades have indicated that compensation of the Z is more localized and/or incomplete. In the ZW Lepidoptera, some species show complete compensation of the Z chromosome, while others lack full equalization, but what drives these inconsistencies is unclear. Here, we compare patterns of male and female gene expression on the Z chromosome of two closely related butterfly species, Papilio xuthus and Papilio machaon, and in multiple tissues of two moths species, Plodia interpunctella and Bombyx mori, which were previously found to differ in the extent to which they equalize Z-linked gene expression between the sexes. We find that, while some species and tissues seem to have incomplete dosage compensation, this is in fact due to the accumulation of male-biased genes and the depletion of female-biased genes on the Z chromosome. Once this is accounted for, the Z chromosome is fully compensated in all four species, through the up-regulation of Z expression in females and in some cases additional down-regulation in males. We further find that both sex-biased genes and Z-linked genes have increased rates of expression divergence in this clade, and that this can lead to fast shifts in patterns of gene expression even between closely related species. Taken together, these results show that the uneven distribution of sex-biased genes on sex chromosomes can confound conclusions about dosage compensation and that Z chromosome-wide dosage compensation is not only possible but ubiquitous among Lepidoptera.}, author = {Huylmans, Ann K and Macon, Ariana and Vicoso, Beatriz}, issn = {07374038}, journal = {Molecular Biology and Evolution}, number = {10}, pages = {2637 -- 2649}, publisher = {Oxford University Press}, title = {{Global dosage compensation is ubiquitous in Lepidoptera, but counteracted by the masculinization of the Z chromosome}}, doi = {10.1093/molbev/msx190}, volume = {34}, year = {2017}, } @article{946, abstract = {Roots navigate through soil integrating environmental signals to orient their growth. The Arabidopsis root is a widely used model for developmental, physiological and cell biological studies. Live imaging greatly aids these efforts, but the horizontal sample position and continuous root tip displacement present significant difficulties. Here, we develop a confocal microscope setup for vertical sample mounting and integrated directional illumination. We present TipTracker – a custom software for automatic tracking of diverse moving objects usable on various microscope setups. Combined, this enables observation of root tips growing along the natural gravity vector over prolonged periods of time, as well as the ability to induce rapid gravity or light stimulation. We also track migrating cells in the developing zebrafish embryo, demonstrating the utility of this system in the acquisition of high-resolution data sets of dynamic samples. We provide detailed descriptions of the tools enabling the easy implementation on other microscopes.}, author = {Von Wangenheim, Daniel and Hauschild, Robert and Fendrych, Matyas and Barone, Vanessa and Benková, Eva and Friml, Jirí}, journal = {eLife}, publisher = {eLife Sciences Publications}, title = {{Live tracking of moving samples in confocal microscopy for vertically grown roots}}, doi = {10.7554/eLife.26792}, volume = {6}, year = {2017}, } @article{947, abstract = {Viewing the ways a living cell can organize its metabolism as the phase space of a physical system, regulation can be seen as the ability to reduce the entropy of that space by selecting specific cellular configurations that are, in some sense, optimal. Here we quantify the amount of regulation required to control a cell's growth rate by a maximum-entropy approach to the space of underlying metabolic phenotypes, where a configuration corresponds to a metabolic flux pattern as described by genome-scale models. We link the mean growth rate achieved by a population of cells to the minimal amount of metabolic regulation needed to achieve it through a phase diagram that highlights how growth suppression can be as costly (in regulatory terms) as growth enhancement. Moreover, we provide an interpretation of the inverse temperature β controlling maximum-entropy distributions based on the underlying growth dynamics. Specifically, we show that the asymptotic value of β for a cell population can be expected to depend on (i) the carrying capacity of the environment, (ii) the initial size of the colony, and (iii) the probability distribution from which the inoculum was sampled. Results obtained for E. coli and human cells are found to be remarkably consistent with empirical evidence.}, author = {De Martino, Daniele and Capuani, Fabrizio and De Martino, Andrea}, issn = {24700045}, journal = { Physical Review E Statistical Nonlinear and Soft Matter Physics }, number = {1}, publisher = {American Institute of Physics}, title = {{Quantifying the entropic cost of cellular growth control}}, doi = {10.1103/PhysRevE.96.010401}, volume = {96}, year = {2017}, } @inproceedings{949, abstract = {The notion of treewidth of graphs has been exploited for faster algorithms for several problems arising in verification and program analysis. Moreover, various notions of balanced tree decompositions have been used for improved algorithms supporting dynamic updates and analysis of concurrent programs. In this work, we present a tool for constructing tree-decompositions of CFGs obtained from Java methods, which is implemented as an extension to the widely used Soot framework. The experimental results show that our implementation on real-world Java benchmarks is very efficient. Our tool also provides the first implementation for balancing tree-decompositions. In summary, we present the first tool support for exploiting treewidth in the static analysis problems on Java programs.}, author = {Chatterjee, Krishnendu and Goharshady, Amir and Pavlogiannis, Andreas}, editor = {D'Souza, Deepak}, issn = {03029743}, location = {Pune, India}, pages = {59 -- 66}, publisher = {Springer}, title = {{JTDec: A tool for tree decompositions in soot}}, doi = {10.1007/978-3-319-68167-2_4}, volume = {10482}, year = {2017}, } @inproceedings{950, abstract = {Two-player games on graphs are widely studied in formal methods as they model the interaction between a system and its environment. The game is played by moving a token throughout a graph to produce an infinite path. There are several common modes to determine how the players move the token through the graph; e.g., in turn-based games the players alternate turns in moving the token. We study the bidding mode of moving the token, which, to the best of our knowledge, has never been studied in infinite-duration games. Both players have separate budgets, which sum up to $1$. In each turn, a bidding takes place. Both players submit bids simultaneously, and a bid is legal if it does not exceed the available budget. The winner of the bidding pays his bid to the other player and moves the token. For reachability objectives, repeated bidding games have been studied and are called Richman games. There, a central question is the existence and computation of threshold budgets; namely, a value t\in [0,1] such that if\PO's budget exceeds $t$, he can win the game, and if\PT's budget exceeds 1-t, he can win the game. We focus on parity games and mean-payoff games. We show the existence of threshold budgets in these games, and reduce the problem of finding them to Richman games. We also determine the strategy-complexity of an optimal strategy. Our most interesting result shows that memoryless strategies suffice for mean-payoff bidding games. }, author = {Avni, Guy and Henzinger, Thomas A and Chonev, Ventsislav K}, issn = {1868-8969}, location = {Berlin, Germany}, publisher = {Schloss Dagstuhl - Leibniz-Zentrum für Informatik}, title = {{Infinite-duration bidding games}}, doi = {10.4230/LIPIcs.CONCUR.2017.21}, volume = {85}, year = {2017}, } @article{9506, abstract = {Methylation in the bodies of active genes is common in animals and vascular plants. Evolutionary patterns indicate homeostatic functions for this type of methylation.}, author = {Zilberman, Daniel}, issn = {1465-6906}, journal = {Genome Biology}, number = {1}, publisher = {Springer Nature}, title = {{An evolutionary case for functional gene body methylation in plants and animals}}, doi = {10.1186/s13059-017-1230-2}, volume = {18}, year = {2017}, } @article{951, abstract = {Dengue-suppressing Wolbachia strains are promising tools for arbovirus control, particularly as they have the potential to self-spread following local introductions. To test this, we followed the frequency of the transinfected Wolbachia strain wMel through Ae. aegypti in Cairns, Australia, following releases at 3 nonisolated locations within the city in early 2013. Spatial spread was analysed graphically using interpolation and by fitting a statistical model describing the position and width of the wave. For the larger 2 of the 3 releases (covering 0.97 km2 and 0.52 km2), we observed slow but steady spatial spread, at about 100–200 m per year, roughly consistent with theoretical predictions. In contrast, the smallest release (0.11 km2) produced erratic temporal and spatial dynamics, with little evidence of spread after 2 years. This is consistent with the prediction concerning fitness-decreasing Wolbachia transinfections that a minimum release area is needed to achieve stable local establishment and spread in continuous habitats. Our graphical and likelihood analyses produced broadly consistent estimates of wave speed and wave width. Spread at all sites was spatially heterogeneous, suggesting that environmental heterogeneity will affect large-scale Wolbachia transformations of urban mosquito populations. The persistence and spread of Wolbachia in release areas meeting minimum area requirements indicates the promise of successful large-scale population transfo}, author = {Schmidt, Tom and Barton, Nicholas H and Rasic, Gordana and Turley, Andrew and Montgomery, Brian and Iturbe Ormaetxe, Inaki and Cook, Peter and Ryan, Peter and Ritchie, Scott and Hoffmann, Ary and O’Neill, Scott and Turelli, Michael}, issn = {15449173}, journal = {PLoS Biology}, number = {5}, publisher = {Public Library of Science}, title = {{Local introduction and heterogeneous spatial spread of dengue-suppressing Wolbachia through an urban population of Aedes Aegypti}}, doi = {10.1371/journal.pbio.2001894}, volume = {15}, year = {2017}, } @article{952, abstract = {A novel strategy for controlling the spread of arboviral diseases such as dengue, Zika and chikungunya is to transform mosquito populations with virus-suppressing Wolbachia. In general, Wolbachia transinfected into mosquitoes induce fitness costs through lower viability or fecundity. These maternally inherited bacteria also produce a frequency-dependent advantage for infected females by inducing cytoplasmic incompatibility (CI), which kills the embryos produced by uninfected females mated to infected males. These competing effects, a frequency-dependent advantage and frequency-independent costs, produce bistable Wolbachia frequency dynamics. Above a threshold frequency, denoted pˆ, CI drives fitness-decreasing Wolbachia transinfections through local populations; but below pˆ, infection frequencies tend to decline to zero. If pˆ is not too high, CI also drives spatial spread once infections become established over sufficiently large areas. We illustrate how simple models provide testable predictions concerning the spatial and temporal dynamics of Wolbachia introductions, focusing on rate of spatial spread, the shape of spreading waves, and the conditions for initiating spread from local introductions. First, we consider the robustness of diffusion-based predictions to incorporating two important features of wMel-Aedes aegypti biology that may be inconsistent with the diffusion approximations, namely fast local dynamics induced by complete CI (i.e., all embryos produced from incompatible crosses die) and long-tailed, non-Gaussian dispersal. With complete CI, our numerical analyses show that long-tailed dispersal changes wave-width predictions only slightly; but it can significantly reduce wave speed relative to the diffusion prediction; it also allows smaller local introductions to initiate spatial spread. Second, we use approximations for pˆ and dispersal distances to predict the outcome of 2013 releases of wMel-infected Aedes aegypti in Cairns, Australia, Third, we describe new data from Ae. aegypti populations near Cairns, Australia that demonstrate long-distance dispersal and provide an approximate lower bound on pˆ for wMel in northeastern Australia. Finally, we apply our analyses to produce operational guidelines for efficient transformation of vector populations over large areas. We demonstrate that even very slow spatial spread, on the order of 10-20 m/month (as predicted), can produce area-wide population transformation within a few years following initial releases covering about 20-30% of the target area.}, author = {Turelli, Michael and Barton, Nicholas H}, issn = {00405809}, journal = {Theoretical Population Biology}, pages = {45 -- 60}, publisher = {Elsevier}, title = {{Deploying dengue-suppressing Wolbachia: Robust models predict slow but effective spatial spread in Aedes aegypti}}, doi = {10.1016/j.tpb.2017.03.003}, volume = {115}, year = {2017}, }