@article{8586, abstract = {Cryo-electron microscopy (cryo-EM) of cellular specimens provides insights into biological processes and structures within a native context. However, a major challenge still lies in the efficient and reproducible preparation of adherent cells for subsequent cryo-EM analysis. This is due to the sensitivity of many cellular specimens to the varying seeding and culturing conditions required for EM experiments, the often limited amount of cellular material and also the fragility of EM grids and their substrate. Here, we present low-cost and reusable 3D printed grid holders, designed to improve specimen preparation when culturing challenging cellular samples directly on grids. The described grid holders increase cell culture reproducibility and throughput, and reduce the resources required for cell culturing. We show that grid holders can be integrated into various cryo-EM workflows, including micro-patterning approaches to control cell seeding on grids, and for generating samples for cryo-focused ion beam milling and cryo-electron tomography experiments. Their adaptable design allows for the generation of specialized grid holders customized to a large variety of applications.}, author = {Fäßler, Florian and Zens, Bettina and Hauschild, Robert and Schur, Florian KM}, issn = {1047-8477}, journal = {Journal of Structural Biology}, keywords = {electron microscopy, cryo-EM, EM sample preparation, 3D printing, cell culture}, number = {3}, publisher = {Elsevier}, title = {{3D printed cell culture grid holders for improved cellular specimen preparation in cryo-electron microscopy}}, doi = {10.1016/j.jsb.2020.107633}, volume = {212}, year = {2020}, } @article{6418, abstract = {Males and females of Artemia franciscana, a crustacean commonly used in the aquarium trade, are highly dimorphic. Sex is determined by a pair of ZW chromosomes, but the nature and extent of differentiation of these chromosomes is unknown. Here, we characterize the Z chromosome by detecting genomic regions that show lower genomic coverage in female than in male samples, and regions that harbor an excess of female-specific SNPs. We detect many Z-specific genes, which no longer have homologs on the W, but also Z-linked genes that appear to have diverged very recently from their existing W-linked homolog. We assess patterns of male and female expression in two tissues with extensive morphological dimorphism, gonads, and heads. In agreement with their morphology, sex-biased expression is common in both tissues. Interestingly, the Z chromosome is not enriched for sex-biased genes, and seems to in fact have a mechanism of dosage compensation that leads to equal expression in males and in females. Both of these patterns are contrary to most ZW systems studied so far, making A. franciscana an excellent model for investigating the interplay between the evolution of sexual dimorphism and dosage compensation, as well as Z chromosome evolution in general.}, author = {Huylmans, Ann K and Toups, Melissa A and Macon, Ariana and Gammerdinger, William J and Vicoso, Beatriz}, issn = {1759-6653}, journal = {Genome biology and evolution}, number = {4}, pages = {1033--1044}, publisher = {Oxford University Press}, title = {{Sex-biased gene expression and dosage compensation on the Artemia franciscana Z-chromosome}}, doi = {10.1093/gbe/evz053}, volume = {11}, year = {2019}, } @article{6442, abstract = {This paper investigates the use of fundamental solutions for animating detailed linear water surface waves. We first propose an analytical solution for efficiently animating circular ripples in closed form. We then show how to adapt the method of fundamental solutions (MFS) to create ambient waves interacting with complex obstacles. Subsequently, we present a novel wavelet-based discretization which outperforms the state of the art MFS approach for simulating time-varying water surface waves with moving obstacles. Our results feature high-resolution spatial details, interactions with complex boundaries, and large open ocean domains. Our method compares favorably with previous work as well as known analytical solutions. We also present comparisons between our method and real world examples.}, author = {Schreck, Camille and Hafner, Christian and Wojtan, Christopher J}, journal = {ACM Transactions on Graphics}, number = {4}, publisher = {ACM}, title = {{Fundamental solutions for water wave animation}}, doi = {10.1145/3306346.3323002}, volume = {38}, year = {2019}, } @article{6859, abstract = {V (vacuolar)/A (archaeal)-type adenosine triphosphatases (ATPases), found in archaeaand eubacteria, couple ATP hydrolysis or synthesis to proton translocation across theplasma membrane using the rotary-catalysis mechanism. They belong to the V-typeATPase family, which differs from the mitochondrial/chloroplast F-type ATP synthasesin overall architecture. We solved cryo–electron microscopy structures of the intactThermus thermophilusV/A-ATPase, reconstituted into lipid nanodiscs, in three rotationalstates and two substates. These structures indicate substantial flexibility betweenV1and Voin a working enzyme, which results from mechanical competition between centralshaft rotation and resistance from the peripheral stalks. We also describedetails of adenosine diphosphate inhibition release, V1-Votorque transmission, andproton translocation, which are relevant for the entire V-type ATPase family.}, author = {Zhou, Long and Sazanov, Leonid A}, issn = {1095-9203}, journal = {Science}, number = {6455}, publisher = {AAAS}, title = {{Structure and conformational plasticity of the intact Thermus thermophilus V/A-type ATPase}}, doi = {10.1126/science.aaw9144}, volume = {365}, year = {2019}, } @article{6978, abstract = {In pipes and channels, the onset of turbulence is initially dominated by localizedtransients, which lead to sustained turbulence through their collective dynamics. In thepresent work, we study numerically the localized turbulence in pipe flow and elucidate astate space structure that gives rise to transient chaos. Starting from the basin boundaryseparating laminar and turbulent flow, we identify transverse homoclinic orbits, thepresence of which necessitates a homoclinic tangle and chaos. A direct consequence ofthe homoclinic tangle is the fractal nature of the laminar-turbulent boundary, which wasconjectured in various earlier studies. By mapping the transverse intersections between thestable and unstable manifold of a periodic orbit, we identify the gateways that promote anescape from turbulence.}, author = {Budanur, Nazmi B and Dogra, Akshunna and Hof, Björn}, journal = {Physical Review Fluids}, number = {10}, pages = {102401}, publisher = {American Physical Society}, title = {{Geometry of transient chaos in streamwise-localized pipe flow turbulence}}, doi = {10.1103/PhysRevFluids.4.102401}, volume = {4}, year = {2019}, } @article{6848, abstract = {Proton-translocating transhydrogenase (also known as nicotinamide nucleotide transhydrogenase (NNT)) is found in the plasma membranes of bacteria and the inner mitochondrial membranes of eukaryotes. NNT catalyses the transfer of a hydride between NADH and NADP+, coupled to the translocation of one proton across the membrane. Its main physiological function is the generation of NADPH, which is a substrate in anabolic reactions and a regulator of oxidative status; however, NNT may also fine-tune the Krebs cycle1,2. NNT deficiency causes familial glucocorticoid deficiency in humans and metabolic abnormalities in mice, similar to those observed in type II diabetes3,4. The catalytic mechanism of NNT has been proposed to involve a rotation of around 180° of the entire NADP(H)-binding domain that alternately participates in hydride transfer and proton-channel gating. However, owing to the lack of high-resolution structures of intact NNT, the details of this process remain unclear5,6. Here we present the cryo-electron microscopy structure of intact mammalian NNT in different conformational states. We show how the NADP(H)-binding domain opens the proton channel to the opposite sides of the membrane, and we provide structures of these two states. We also describe the catalytically important interfaces and linkers between the membrane and the soluble domains and their roles in nucleotide exchange. These structures enable us to propose a revised mechanism for a coupling process in NNT that is consistent with a large body of previous biochemical work. Our results are relevant to the development of currently unavailable NNT inhibitors, which may have therapeutic potential in ischaemia reperfusion injury, metabolic syndrome and some cancers7,8,9.}, author = {Kampjut, Domen and Sazanov, Leonid A}, issn = {1476-4687}, journal = {Nature}, number = {7773}, pages = {291–295}, publisher = {Springer Nature}, title = {{Structure and mechanism of mitochondrial proton-translocating transhydrogenase}}, doi = {10.1038/s41586-019-1519-2}, volume = {573}, year = {2019}, } @phdthesis{6435, abstract = {Social insect colonies tend to have numerous members which function together like a single organism in such harmony that the term ``super-organism'' is often used. In this analogy the reproductive caste is analogous to the primordial germ cells of a metazoan, while the sterile worker caste corresponds to somatic cells. The worker castes, like tissues, are in charge of all functions of a living being, besides reproduction. The establishment of new super-organismal units (i.e. new colonies) is accomplished by the co-dependent castes. The term oftentimes goes beyond a metaphor. We invoke it when we speak about the metabolic rate, thermoregulation, nutrient regulation and gas exchange of a social insect colony. Furthermore, we assert that the super-organism has an immune system, and benefits from ``social immunity''. Social immunity was first summoned by evolutionary biologists to resolve the apparent discrepancy between the expected high frequency of disease outbreak amongst numerous, closely related tightly-interacting hosts, living in stable and microbially-rich environments, against the exceptionally scarce epidemic accounts in natural populations. Social immunity comprises a multi-layer assembly of behaviours which have evolved to effectively keep the pathogenic enemies of a colony at bay. The field of social immunity has drawn interest, as it becomes increasingly urgent to stop the collapse of pollinator species and curb the growth of invasive pests. In the past decade, several mechanisms of social immune responses have been dissected, but many more questions remain open. I present my work in two experimental chapters. In the first, I use invasive garden ants (*Lasius neglectus*) to study how pathogen load and its distribution among nestmates affect the grooming response of the group. Any given group of ants will carry out the same total grooming work, but will direct their grooming effort towards individuals carrying a relatively higher spore load. Contrary to expectation, the highest risk of transmission does not stem from grooming highly contaminated ants, but instead, we suggest that the grooming response likely minimizes spore loss to the environment, reducing contamination from inadvertent pickup from the substrate. The second is a comparative developmental approach. I follow black garden ant queens (*Lasius niger*) and their colonies from mating flight, through hibernation for a year. Colonies which grow fast from the start, have a lower chance of survival through hibernation, and those which survive grow at a lower pace later. This is true for colonies of naive and challenged queens. Early pathogen exposure of the queens changes colony dynamics in an unexpected way: colonies from exposed queens are more likely to grow slowly and recover in numbers only after they survive hibernation. In addition to the two experimental chapters, this thesis includes a co-authored published review on organisational immunity, where we enlist the experimental evidence and theoretical framework on which this hypothesis is built, identify the caveats and underline how the field is ripe to overcome them. In a final chapter, I describe my part in two collaborative efforts, one to develop an image-based tracker, and the second to develop a classifier for ant behaviour.}, author = {Casillas Perez, Barbara E}, issn = {2663-337X}, keywords = {Social Immunity, Sanitary care, Social Insects, Organisational Immunity, Colony development, Multi-target tracking}, pages = {183}, publisher = {Institute of Science and Technology Austria}, title = {{Collective defenses of garden ants against a fungal pathogen}}, doi = {10.15479/AT:ISTA:6435}, year = {2019}, } @article{134, abstract = {The current state of the art in real-time two-dimensional water wave simulation requires developers to choose between efficient Fourier-based methods, which lack interactions with moving obstacles, and finite-difference or finite element methods, which handle environmental interactions but are significantly more expensive. This paper attempts to bridge this long-standing gap between complexity and performance, by proposing a new wave simulation method that can faithfully simulate wave interactions with moving obstacles in real time while simultaneously preserving minute details and accommodating very large simulation domains. Previous methods for simulating 2D water waves directly compute the change in height of the water surface, a strategy which imposes limitations based on the CFL condition (fast moving waves require small time steps) and Nyquist's limit (small wave details require closely-spaced simulation variables). This paper proposes a novel wavelet transformation that discretizes the liquid motion in terms of amplitude-like functions that vary over space, frequency, and direction, effectively generalizing Fourier-based methods to handle local interactions. Because these new variables change much more slowly over space than the original water height function, our change of variables drastically reduces the limitations of the CFL condition and Nyquist limit, allowing us to simulate highly detailed water waves at very large visual resolutions. Our discretization is amenable to fast summation and easy to parallelize. We also present basic extensions like pre-computed wave paths and two-way solid fluid coupling. Finally, we argue that our discretization provides a convenient set of variables for artistic manipulation, which we illustrate with a novel wave-painting interface.}, author = {Jeschke, Stefan and Skrivan, Tomas and Mueller Fischer, Matthias and Chentanez, Nuttapong and Macklin, Miles and Wojtan, Christopher J}, journal = {ACM Transactions on Graphics}, number = {4}, publisher = {ACM}, title = {{Water surface wavelets}}, doi = {10.1145/3197517.3201336}, volume = {37}, year = {2018}, } @article{44, abstract = {Recent realization of a kinetically constrained chain of Rydberg atoms by Bernien et al., [Nature (London) 551, 579 (2017)] resulted in the observation of unusual revivals in the many-body quantum dynamics. In our previous work [C. J. Turner et al., Nat. Phys. 14, 745 (2018)], such dynamics was attributed to the existence of “quantum scarred” eigenstates in the many-body spectrum of the experimentally realized model. Here, we present a detailed study of the eigenstate properties of the same model. We find that the majority of the eigenstates exhibit anomalous thermalization: the observable expectation values converge to their Gibbs ensemble values, but parametrically slower compared to the predictions of the eigenstate thermalization hypothesis (ETH). Amidst the thermalizing spectrum, we identify nonergodic eigenstates that strongly violate the ETH, whose number grows polynomially with system size. Previously, the same eigenstates were identified via large overlaps with certain product states, and were used to explain the revivals observed in experiment. Here, we find that these eigenstates, in addition to highly atypical expectation values of local observables, also exhibit subthermal entanglement entropy that scales logarithmically with the system size. Moreover, we identify an additional class of quantum scarred eigenstates, and discuss their manifestations in the dynamics starting from initial product states. We use forward scattering approximation to describe the structure and physical properties of quantum scarred eigenstates. Finally, we discuss the stability of quantum scars to various perturbations. We observe that quantum scars remain robust when the introduced perturbation is compatible with the forward scattering approximation. In contrast, the perturbations which most efficiently destroy quantum scars also lead to the restoration of “canonical” thermalization.}, author = {Turner, C J and Michailidis, Alexios and Abanin, D A and Serbyn, Maksym and Papić, Z}, journal = {Physical Review B}, number = {15}, publisher = {American Physical Society}, title = {{Quantum scarred eigenstates in a Rydberg atom chain: Entanglement, breakdown of thermalization, and stability to perturbations}}, doi = {10.1103/PhysRevB.98.155134}, volume = {98}, year = {2018}, } @article{1087, abstract = {Using extensive direct numerical simulations, the dynamics of laminar-turbulent fronts in pipe flow is investigated for Reynolds numbers between and 5500. We here investigate the physical distinction between the fronts of weak and strong slugs both by analysing the turbulent kinetic energy budget and by comparing the downstream front motion to the advection speed of bulk turbulent structures. Our study shows that weak downstream fronts travel slower than turbulent structures in the bulk and correspond to decaying turbulence at the front. At the downstream front speed becomes faster than the advection speed, marking the onset of strong fronts. In contrast to weak fronts, turbulent eddies are generated at strong fronts by feeding on the downstream laminar flow. Our study also suggests that temporal fluctuations of production and dissipation at the downstream laminar-turbulent front drive the dynamical switches between the two types of front observed up to.}, author = {Song, Baofang and Barkley, Dwight and Hof, Björn and Avila, Marc}, issn = {0022-1120}, journal = {Journal of Fluid Mechanics}, pages = {1045 -- 1059}, publisher = {Cambridge University Press}, title = {{Speed and structure of turbulent fronts in pipe flow}}, doi = {10.1017/jfm.2017.14}, volume = {813}, year = {2017}, } @article{1118, abstract = {Sharp wave-ripple (SWR) oscillations play a key role in memory consolidation during non-rapid eye movement sleep, immobility, and consummatory behavior. However, whether temporally modulated synaptic excitation or inhibition underlies the ripples is controversial. To address this question, we performed simultaneous recordings of excitatory and inhibitory postsynaptic currents (EPSCs and IPSCs) and local field potentials (LFPs) in the CA1 region of awake mice in vivo. During SWRs, inhibition dominated over excitation, with a peak conductance ratio of 4.1 ± 0.5. Furthermore, the amplitude of SWR-associated IPSCs was positively correlated with SWR magnitude, whereas that of EPSCs was not. Finally, phase analysis indicated that IPSCs were phase-locked to individual ripple cycles, whereas EPSCs were uniformly distributed in phase space. Optogenetic inhibition indicated that PV+ interneurons provided a major contribution to SWR-associated IPSCs. Thus, phasic inhibition, but not excitation, shapes SWR oscillations in the hippocampal CA1 region in vivo.}, author = {Gan, Jian and Weng, Shih-Ming and Pernia-Andrade, Alejandro and Csicsvari, Jozsef L and Jonas, Peter M}, journal = {Neuron}, number = {2}, pages = {308 -- 314}, publisher = {Elsevier}, title = {{Phase-locked inhibition, but not excitation, underlies hippocampal ripple oscillations in awake mice in vivo}}, doi = {10.1016/j.neuron.2016.12.018}, volume = {93}, year = {2017}, } @article{470, abstract = {This paper presents a method for simulating water surface waves as a displacement field on a 2D domain. Our method relies on Lagrangian particles that carry packets of water wave energy; each packet carries information about an entire group of wave trains, as opposed to only a single wave crest. Our approach is unconditionally stable and can simulate high resolution geometric details. This approach also presents a straightforward interface for artistic control, because it is essentially a particle system with intuitive parameters like wavelength and amplitude. Our implementation parallelizes well and runs in real time for moderately challenging scenarios.}, author = {Jeschke, Stefan and Wojtan, Christopher J}, issn = {0730-0301}, journal = {ACM Transactions on Graphics}, number = {4}, publisher = {ACM}, title = {{Water wave packets}}, doi = {10.1145/3072959.3073678}, volume = {36}, year = {2017}, } @article{804, abstract = {Polysaccharides (carbohydrates) are key regulators of a large number of cell biological processes. However, precise biochemical or genetic manipulation of these often complex structures is laborious and hampers experimental structure–function studies. Molecular Dynamics (MD) simulations provide a valuable alternative tool to generate and test hypotheses on saccharide function. Yet, currently used MD force fields often overestimate the aggregation propensity of polysaccharides, affecting the usability of those simulations. Here we tested MARTINI, a popular coarse-grained (CG) force field for biological macromolecules, for its ability to accurately represent molecular forces between saccharides. To this end, we calculated a thermodynamic solution property, the second virial coefficient of the osmotic pressure (B22). Comparison with light scattering experiments revealed a nonphysical aggregation of a prototypical polysaccharide in MARTINI, pointing at an imbalance of the nonbonded solute–solute, solute–water, and water–water interactions. This finding also applies to smaller oligosaccharides which were all found to aggregate in simulations even at moderate concentrations, well below their solubility limit. Finally, we explored the influence of the Lennard-Jones (LJ) interaction between saccharide molecules and propose a simple scaling of the LJ interaction strength that makes MARTINI more reliable for the simulation of saccharides.}, author = {Schmalhorst, Philipp S and Deluweit, Felix and Scherrers, Roger and Heisenberg, Carl-Philipp J and Sikora, Mateusz K}, issn = {1549-9618}, journal = {Journal of Chemical Theory and Computation}, number = {10}, pages = {5039 -- 5053}, publisher = {American Chemical Society}, title = {{Overcoming the limitations of the MARTINI force field in simulations of polysaccharides}}, doi = {10.1021/acs.jctc.7b00374}, volume = {13}, year = {2017}, } @article{727, abstract = {Actin filaments polymerizing against membranes power endocytosis, vesicular traffic, and cell motility. In vitro reconstitution studies suggest that the structure and the dynamics of actin networks respond to mechanical forces. We demonstrate that lamellipodial actin of migrating cells responds to mechanical load when membrane tension is modulated. In a steady state, migrating cell filaments assume the canonical dendritic geometry, defined by Arp2/3-generated 70° branch points. Increased tension triggers a dense network with a broadened range of angles, whereas decreased tension causes a shift to a sparse configuration dominated by filaments growing perpendicularly to the plasma membrane. We show that these responses emerge from the geometry of branched actin: when load per filament decreases, elongation speed increases and perpendicular filaments gradually outcompete others because they polymerize the shortest distance to the membrane, where they are protected from capping. This network-intrinsic geometrical adaptation mechanism tunes protrusive force in response to mechanical load.}, author = {Mueller, Jan and Szep, Gregory and Nemethova, Maria and De Vries, Ingrid and Lieber, Arnon and Winkler, Christoph and Kruse, Karsten and Small, John and Schmeiser, Christian and Keren, Kinneret and Hauschild, Robert and Sixt, Michael K}, issn = {0092-8674}, journal = {Cell}, number = {1}, pages = {188 -- 200}, publisher = {Cell Press}, title = {{Load adaptation of lamellipodial actin networks}}, doi = {10.1016/j.cell.2017.07.051}, volume = {171}, year = {2017}, } @article{1350, abstract = {The hippocampal CA3 region plays a key role in learning and memory. Recurrent CA3–CA3 synapses are thought to be the subcellular substrate of pattern completion. However, the synaptic mechanisms of this network computation remain enigmatic. To investigate these mechanisms, we combined functional connectivity analysis with network modeling. Simultaneous recording fromup to eight CA3 pyramidal neurons revealed that connectivity was sparse, spatially uniform, and highly enriched in disynaptic motifs (reciprocal, convergence,divergence, and chain motifs). Unitary connections were composed of one or two synaptic contacts, suggesting efficient use of postsynaptic space. Real-size modeling indicated that CA3 networks with sparse connectivity, disynaptic motifs, and single-contact connections robustly generated pattern completion.Thus, macro- and microconnectivity contribute to efficient memory storage and retrieval in hippocampal networks.}, author = {Guzmán, José and Schlögl, Alois and Frotscher, Michael and Jonas, Peter M}, journal = {Science}, number = {6304}, pages = {1117 -- 1123}, publisher = {American Association for the Advancement of Science}, title = {{Synaptic mechanisms of pattern completion in the hippocampal CA3 network}}, doi = {10.1126/science.aaf1836}, volume = {353}, year = {2016}, } @inproceedings{1363, abstract = {When aiming to seamlessly integrate a fluid simulation into a larger scenario (like an open ocean), careful attention must be paid to boundary conditions. In particular, one must implement special "non-reflecting" boundary conditions, which dissipate out-going waves as they exit the simulation. Unfortunately, the state of the art in non-reflecting boundary conditions (perfectly-matched layers, or PMLs) only permits trivially simple inflow/outflow conditions, so there is no reliable way to integrate a fluid simulation into a more complicated environment like a stormy ocean or a turbulent river. This paper introduces the first method for combining nonreflecting boundary conditions based on PMLs with inflow/outflow boundary conditions that vary arbitrarily throughout space and time. Our algorithm is a generalization of stateof- the-art mean-flow boundary conditions in the computational fluid dynamics literature, and it allows for seamless integration of a fluid simulation into much more complicated environments. Our method also opens the door for previously-unseen postprocess effects like retroactively changing the location of solid obstacles, and locally increasing the visual detail of a pre-existing simulation.}, author = {Bojsen-Hansen, Morten and Wojtan, Christopher J}, location = {Anaheim, CA, USA}, number = {4}, publisher = {ACM}, title = {{Generalized non-reflecting boundaries for fluid re-simulation}}, doi = {10.1145/2897824.2925963}, volume = {35}, year = {2016}, } @article{2944, abstract = {We propose a two-step procedure for estimating multiple migration rates in an approximate Bayesian computation (ABC) framework, accounting for global nuisance parameters. The approach is not limited to migration, but generally of interest for inference problems with multiple parameters and a modular structure (e.g. independent sets of demes or loci). We condition on a known, but complex demographic model of a spatially subdivided population, motivated by the reintroduction of Alpine ibex (Capra ibex) into Switzerland. In the first step, the global parameters ancestral mutation rate and male mating skew have been estimated for the whole population in Aeschbacher et al. (Genetics 2012; 192: 1027). In the second step, we estimate in this study the migration rates independently for clusters of demes putatively connected by migration. For large clusters (many migration rates), ABC faces the problem of too many summary statistics. We therefore assess by simulation if estimation per pair of demes is a valid alternative. We find that the trade-off between reduced dimensionality for the pairwise estimation on the one hand and lower accuracy due to the assumption of pairwise independence on the other depends on the number of migration rates to be inferred: the accuracy of the pairwise approach increases with the number of parameters, relative to the joint estimation approach. To distinguish between low and zero migration, we perform ABC-type model comparison between a model with migration and one without. Applying the approach to microsatellite data from Alpine ibex, we find no evidence for substantial gene flow via migration, except for one pair of demes in one direction.}, author = {Aeschbacher, Simon and Futschik, Andreas and Beaumont, Mark}, journal = {Molecular Ecology}, number = {4}, pages = {987 -- 1002}, publisher = {Wiley-Blackwell}, title = {{Approximate Bayesian computation for modular inference problems with many parameters: the example of migration rates. }}, doi = {10.1111/mec.12165}, volume = {22}, year = {2013}, } @article{2962, abstract = {The choice of summary statistics is a crucial step in approximate Bayesian computation (ABC). Since statistics are often not sufficient, this choice involves a trade-off between loss of information and reduction of dimensionality. The latter may increase the efficiency of ABC. Here, we propose an approach for choosing summary statistics based on boosting, a technique from the machine learning literature. We consider different types of boosting and compare them to partial least squares regression as an alternative. To mitigate the lack of sufficiency, we also propose an approach for choosing summary statistics locally, in the putative neighborhood of the true parameter value. We study a demographic model motivated by the re-introduction of Alpine ibex (Capra ibex) into the Swiss Alps. The parameters of interest are the mean and standard deviation across microsatellites of the scaled ancestral mutation rate (θanc = 4 Ne u), and the proportion of males obtaining access to matings per breeding season (ω). By simulation, we assess the properties of the posterior distribution obtained with the various methods. According to our criteria, ABC with summary statistics chosen locally via boosting with the L2-loss performs best. Applying that method to the ibex data, we estimate θanc ≈ 1.288, and find that most of the variation across loci of the ancestral mutation rate u is between 7.7×10−4 and 3.5×10−3 per locus per generation. The proportion of males with access to matings is estimated to ω ≈ 0.21, which is in good agreement with recent independent estimates.}, author = {Aeschbacher, Simon and Beaumont, Mark and Futschik, Andreas}, journal = {Genetics}, number = {3}, pages = {1027 -- 1047}, publisher = {Genetics Society of America}, title = {{A novel approach for choosing summary statistics in approximate Bayesian computation}}, doi = {10.1534/genetics.112.143164}, volume = {192}, year = {2012}, }