Search

We consider a homogeneous birth-death process with three different sampling schemes. First, individuals can be sampled through time and included in a reconstructed phylogenetic tree. Second, they can be sampled through time and only recorded as a point 'occurrence' along a timeline. Third, extant individuals can be sampled and included in the reconstructed phylogenetic tree with a fixed probability. We further consider that sampled individuals can be removed or not from the process, upon sampling, with fixed probability. We derive the probability distribution of the population size at any time in the past conditional on the joint observation of a reconstructed phylogenetic tree and a record of occurrences not included in the tree. We also provide an algorithm to simulate ancestral population size trajectories given the observation of a reconstructed phylogenetic tree and occurrences. This distribution can be readily used to draw inferences about the ancestral population size in the field of epidemiology and macroevolution. In epidemiology, these results will allow data from epidemiological case count studies to be used in conjunction with molecular sequencing data (yielding reconstructed phylogenetic trees) to coherently estimate prevalence through time. In macroevolution, it will foster the joint examination of the fossil record and extant taxa to reconstruct past biodiversity., (Copyright © 2020 The Author(s). Published by Elsevier Ltd.. All rights reserved.)

In the wake of the pandemic of coronavirus disease 2019 (COVID-19), contact tracing has become a key element of strategies to control the spread of severe acute respiratory syndrome coronavirus-2 (SARS-CoV-2). Given the rapid and intense spread of SARS-CoV-2, digital contact tracing has emerged as a potential complementary tool to support containment and mitigation efforts. Early modelling studies highlighted the potential of digital contact tracing to break transmission chains, and Google and Apple subsequently developed the Exposure Notification (EN) framework, making it available to the vast majority of smartphones. A growing number of governments have launched or announced EN-based contact tracing apps, but their effectiveness remains unknown. Here, we report early findings of the digital contact tracing app deployment in Switzerland. We demonstrate proof-of-principle that digital contact tracing reaches exposed contacts, who then test positive for SARS-CoV-2. This indicates that digital contact tracing is an effective complementary tool for controlling the spread of SARS-CoV-2. Continued technical improvement and international compatibility can further increase the efficacy, particularly also across country borders.

Emerging infectious diseases pose one of the greatest threats to human health and biodiversity. Phylodynamics is often used to infer epidemiological parameters essential for guiding intervention strategies for human viruses such as severe acute respiratory syndrome coronavirus 2 (SARS-Cov-2). Here, we applied phylodynamics to elucidate the epidemiological dynamics of Tasmanian devil facial tumor disease (DFTD), a fatal, transmissible cancer with a genome thousands of times larger than that of any virus. Despite prior predictions of devil extinction, transmission rates have declined precipitously from ~3.5 secondary infections per infected individual to ~1 at present. Thus, DFTD appears to be transitioning from emergence to endemism, lending hope for the continued survival of the endangered Tasmanian devil. More generally, our study demonstrates a new phylodynamic analytical framework that can be applied to virtually any pathogen., (Copyright © 2020 The Authors, some rights reserved; exclusive licensee American Association for the Advancement of Science. No claim to original U.S. Government Works.)

Estimation of the effective reproductive number Rt is important for detecting changes in disease transmission over time. During the Coronavirus Disease 2019 (COVID-19) pandemic, policy makers and public health officials are using Rt to assess the effectiveness of interventions and to inform policy. However, estimation of Rt from available data presents several challenges, with critical implications for the interpretation of the course of the pandemic. The purpose of this document is to summarize these challenges, illustrate them with examples from synthetic data, and, where possible, make recommendations. For near real-time estimation of Rt, we recommend the approach of Cori and colleagues, which uses data from before time t and empirical estimates of the distribution of time between infections. Methods that require data from after time t, such as Wallinga and Teunis, are conceptually and methodologically less suited for near real-time estimation, but may be appropriate for retrospective analyses of how individuals infected at different time points contributed to the spread. We advise caution when using methods derived from the approach of Bettencourt and Ribeiro, as the resulting Rt estimates may be biased if the underlying structural assumptions are not met. Two key challenges common to all approaches are accurate specification of the generation interval and reconstruction of the time series of new infections from observations occurring long after the moment of transmission. Naive approaches for dealing with observation delays, such as subtracting delays sampled from a distribution, can introduce bias. We provide suggestions for how to mitigate this and other technical challenges and highlight open problems in Rt estimation., Competing Interests: The authors have declared that no competing interests exist.

Infecting large portions of the global population, seasonal influenza is a major burden on societies around the globe. While the global source sink dynamics of the different seasonal influenza viruses have been studied intensively, its local spread remains less clear. In order to improve our understanding of how influenza is transmitted on a city scale, we collected an extremely densely sampled set of influenza sequences alongside patient metadata. To do so, we sequenced influenza viruses isolated from patients of two different hospitals, as well as private practitioners in Basel, Switzerland during the 2016/2017 influenza season. The genetic sequences reveal that repeated introductions into the city drove the influenza season. We then reconstruct how the effective reproduction number changed over the course of the season. While we did not find that transmission dynamics in Basel correlate with humidity or school closures, we did find some evidence that it may positively correlated with temperature. Alongside the genetic sequence data that allows us to see how individual cases are connected, we gathered patient information, such as the age or household status. Zooming into the local transmission outbreaks suggests that the elderly were to a large extent infected within their own transmission network. In the remaining transmission network, our analyses suggest that school-aged children likely play a more central role than pre-school aged children. These patterns will be valuable to plan interventions combating the spread of respiratory diseases within cities given that similar patterns are observed for other influenza seasons and cities., Competing Interests: The authors have declared that no competing interests exist.

Phylogenetic methods can use the sampling times of molecular sequence data to calibrate the molecular clock, enabling the estimation of evolutionary rates and timescales for rapidly evolving pathogens and data sets containing ancient DNA samples. A key aspect of such calibrations is whether a sufficient amount of molecular evolution has occurred over the sampling time window, that is, whether the data can be treated as having come from a measurably evolving population. Here, we investigate the performance of a fully Bayesian evaluation of temporal signal (BETS) in sequence data. The method involves comparing the fit to the data of two models: a model in which the data are accompanied by the actual (heterochronous) sampling times, and a model in which the samples are constrained to be contemporaneous (isochronous). We conducted simulations under a wide range of conditions to demonstrate that BETS accurately classifies data sets according to whether they contain temporal signal or not, even when there is substantial among-lineage rate variation. We explore the behavior of this classification in analyses of five empirical data sets: modern samples of A/H1N1 influenza virus, the bacterium Bordetella pertussis, coronaviruses from mammalian hosts, ancient DNA from Hepatitis B virus, and mitochondrial genomes of dog species. Our results indicate that BETS is an effective alternative to other tests of temporal signal. In particular, this method has the key advantage of allowing a coherent assessment of the entire model, including the molecular clock and tree prior which are essential aspects of Bayesian phylodynamic analyses., (© The Author(s) 2020. Published by Oxford University Press on behalf of the Society for Molecular Biology and Evolution. All rights reserved. For permissions, please e-mail: journals.permissions@oup.com.)

Aims: To assess whether the combination of transthoracic echocardiography (TTE) and multidetector computed tomography (MDCT) data affects the grading of aortic stenosis (AS) severity under consideration of the energy loss index (ELI) in patients undergoing transcatheter aortic valve replacement (TAVR)., Methods and Results: Multimodality imaging was performed in 197 patients with symptomatic severe AS undergoing TAVR at the University Hospital Zurich, Switzerland. Fusion aortic valve area index (fusion AVAi) assessed by integrating MDCT derived planimetric left ventricular outflow tract area into the continuity equation was significantly larger as compared to conventional AVAi (0.41 ± 0.1 vs. 0.51 ± 0.1 cm2/m2; P < 0.01). A total of 62 patients (31.4%) were reclassified from severe to moderate AS with fusion AVAi being >0.6 cm2/m2. ELI was obtained for conventional AVAi and fusion AVAi based on sinotubular junction area determined by TTE (ELILTL 0.47 ± 0.1 cm2/m2; fusion ELILTL 0.60 ± 0.1 cm2/m2) and MDCT (ELIMDCT 0.48 ± 0.1 cm2/m2; fusion ELIMDCT 0.61 ± 0.05 cm2/m2). When ELI was calculated with fusion AVAi the effective orifice area was >0.6 cm2/m2 in 85 patients (43.1%). Survival rate 3 years after TAVR was higher in patients reclassified to moderate AS according to multimodality imaging derived ELI (78.8% vs. 67%; P = 0.01)., Conclusion: Multimodality imaging derived ELI reclassifies AS severity in 43% undergoing TAVR and predicts mid-term outcome., (Published on behalf of the European Society of Cardiology. All rights reserved. © The Author(s) 2020. For permissions, please email: journals.permissions@oup.com.)

Next-generation sequencing of pathogen quasispecies within a host yields data sets of tens to hundreds of unique sequences. However, the full data set often contains thousands of sequences, because many of those unique sequences have multiple identical copies. Data sets of this size represent a computational challenge for currently available Bayesian phylogenetic and phylodynamic methods. Through simulations, we explore how large data sets with duplicate sequences affect the speed and accuracy of phylogenetic and phylodynamic analysis within BEAST 2. We show that using unique sequences only leads to biases, and using a random subset of sequences yields imprecise parameter estimates. To overcome these shortcomings, we introduce PIQMEE, a BEAST 2 add-on that produces reliable parameter estimates from full data sets with increased computational efficiency as compared with the currently available methods within BEAST 2. The principle behind PIQMEE is to resolve the tree structure of the unique sequences only, while simultaneously estimating the branching times of the duplicate sequences. Distinguishing between unique and duplicate sequences allows our method to perform well even for very large data sets. Although the classic method converges poorly for data sets of 6,000 sequences when allowed to run for 7 days, our method converges in slightly more than 1 day. In fact, PIQMEE can handle data sets of around 21,000 sequences with 20 unique sequences in 14 days. Finally, we apply the method to a real, within-host HIV sequencing data set with several thousand sequences per patient., (© The Author(s) 2020. Published by Oxford University Press on behalf of the Society for Molecular Biology and Evolution.)

Heterogeneous populations can lead to important differences in birth and death rates across a phylogeny. Taking this heterogeneity into account is necessary to obtain accurate estimates of the underlying population dynamics. We present a new multitype birth-death model (MTBD) that can estimate lineage-specific birth and death rates. This corresponds to estimating lineage-dependent speciation and extinction rates for species phylogenies, and lineage-dependent transmission and recovery rates for pathogen transmission trees. In contrast with previous models, we do not presume to know the trait driving the rate differences, nor do we prohibit the same rates from appearing in different parts of the phylogeny. Using simulated data sets, we show that the MTBD model can reliably infer the presence of multiple evolutionary regimes, their positions in the tree, and the birth and death rates associated with each. We also present a reanalysis of two empirical data sets and compare the results obtained by MTBD and by the existing software BAMM. We compare two implementations of the model, one exact and one approximate (assuming that no rate changes occur in the extinct parts of the tree), and show that the approximation only slightly affects results. The MTBD model is implemented as a package in the Bayesian inference software BEAST 2 and allows joint inference of the phylogeny and the model parameters.[Birth-death; lineage specific rates, multi-type model.]., (© The Author(s) 2020. Published by Oxford University Press, on behalf of the Society of Systematic Biologists.)

Reassortment is an important source of genetic diversity in segmented viruses and is the main source of novel pathogenic influenza viruses. Despite this, studying the reassortment process has been constrained by the lack of a coherent, model-based inference framework. Here, we introduce a coalescent-based model that allows us to explicitly model the joint coalescent and reassortment process. In order to perform inference under this model, we present an efficient Markov chain Monte Carlo algorithm to sample rooted networks and the embedding of phylogenetic trees within networks. This algorithm provides the means to jointly infer coalescent and reassortment rates with the reassortment network and the embedding of segments in that network from full-genome sequence data. Studying reassortment patterns of different human influenza datasets, we find large differences in reassortment rates across different human influenza viruses. Additionally, we find that reassortment events predominantly occur on selectively fitter parts of reassortment networks showing that on a population level, reassortment positively contributes to the fitness of human influenza viruses., Competing Interests: The authors declare no competing interest., (Copyright © 2020 the Author(s). Published by PNAS.)

Protection and neutralization of a vast array of pathogens is accomplished by the tremendous diversity of the B cell receptor (BCR) repertoire. For jawed vertebrates, this diversity is initiated via the somatic recombination of immunoglobulin (Ig) germline elements. While it is clear that the number of these germline segments differs from species to species, the extent of cross-species sequence diversity remains largely uncharacterized. Here we use extensive computational and statistical methods to investigate the sequence diversity and evolutionary relationship between Ig variable (V), diversity (D), and joining (J) germline segments across nine commonly studied species ranging from zebrafish to human. Metrics such as guanine-cytosine (GC) content showed low redundancy across Ig germline genes within a given species. Other comparisons, including amino acid motifs, evolutionary selection, and sequence diversity, revealed species-specific properties. Additionally, we showed that the germline-encoded diversity differs across antibody (recombined V-D-J) repertoires of various B cell subsets. To facilitate future comparative immunogenomics analysis, we created VDJgermlines, an R package that contains the germline sequences from multiple species. Our study informs strategies for the humanization and engineering of therapeutic antibodies.

Antibody repertoire sequencing provides a molecular fingerprint of current and past pathogens encountered by the immune system. Most repertoire studies in humans require measuring the B cell response in the blood, resulting in a large bias to the IgM isotype. The extent to which the circulating IgM antibody repertoire correlates to lymphoid tissue-resident B cells in the setting of viral infection remains largely uncharacterized. Therefore, we compared the IgM repertoires from both blood and bone marrow (BM) plasma cells (PCs) following acute or chronic lymphocytic choriomeningitis virus (LCMV) infection in mice. Despite previously reported serum alterations between acute and chronic infection, IgM repertoire signatures based on clonal diversity metrics, public clones, network, and phylogenetic analysis were largely unable to distinguish infection cohorts. Our findings, however, revealed mouse-specific repertoire fingerprints between the blood and PC repertoires irrespective of infection status., (Copyright © 2020 Yermanos, Kräutler, Pedrioli, Menzel, Greiff, Stadler, Oxenius and Reddy.)

The reproductive number in Switzerland was between 1.5 and 2 during the first third of March, and has consistently decreased to around 1. After the announcement of the latest strict measure on 20 March 2020, namely that gatherings of more than five people in public spaces are prohibited, the reproductive number dropped significantly below 1; the authors of this study estimate the reproductive number to be between 0.6 and 0.8 in the first third of April.

We study the problem of computing the probability distribution of phylogenetic trees that commonly arise in areas ranging from epidemiology to macroevolution. We focus on homogeneous birth death trees with incomplete sampling and consider observations from three distinct sampling schemes. First, individuals can be sampled and removed, through time, and included in the tree. Second, they can be occurrences which are sampled and removed through time and not included in the tree. Third, extant individuals can be sampled and included in the tree. The outcome of the process is thus composed of the reconstructed phylogenetic tree spanning all individuals sampled and included in the tree, and a timeline of occurrence events which are not placed along the tree. We derive a formula for computing the joint probability density of this outcome, which can readily be used to perform maximum likelihood or Bayesian estimation of the parameters of the model. In the context of epidemiology, our probability density enables the estimation of transmission rates through a joint analysis of epidemiological case count data and phylogenetic trees reconstructed from pathogen sequences. Within macroevolution, our equations form the basis for incorporating fossil occurrences from paleontological databases together with extant species phylogenies for estimating speciation and extinction rates. This work provides the theoretical framework for bridging not only the gap between phylogenetics and epidemiology, but also that between phylogenetics and paleontology., (Copyright © 2019. Published by Elsevier Ltd.)

Phylogenetic comparative methods (PCMs) have been used to study the evolution of quantitative traits in various groups of organisms, ranging from micro-organisms to animal and plant species. A common approach has been to assume a Gaussian phylogenetic model for the trait evolution along the tree, such as a branching Brownian motion (BM) or an Ornstein-Uhlenbeck (OU) process. Then, the parameters of the process have been inferred based on a given tree and trait data for the sampled species. At the heart of this inference lie multiple calculations of the model likelihood, that is, the probability density of the observed trait data, conditional on the model parameters and the tree. With the increasing availability of big phylogenetic trees, spanning hundreds to several thousand sampled species, this approach is facing a two-fold challenge. First, the assumption of a single Gaussian process governing the entire tree is not adequate in the presence of heterogeneous evolutionary forces acting in different parts of the tree. Second, big trees present a computational challenge, due to the time and memory complexity of the model likelihood calculation. Here, we explore a sub-family, denoted G LInv , of the Gaussian phylogenetic models, with the transition density exhibiting the properties that the expectation depends Linearly on the ancestral trait value and the variance is Invariant with respect to the ancestral value. We show that G LInv contains the vast majority of Gaussian models currently used in PCMs, while supporting an efficient (linear in the number of nodes) algorithm for the likelihood calculation. The algorithm supports scenarios with missing data, as well as different types of trees, including trees with polytomies and non-ultrametric trees. To account for the heterogeneity in the evolutionary forces, the algorithm supports models with "shifts" occurring at specific points in the tree. Such shifts can include changes in some or all parameters, as well as the type of the model, provided that the model remains within the G LInv family. This contrasts with most of the current implementations where, due to slow likelihood calculation, the shifts are restricted to specific parameters in a single type of model, such as the long-term selection optima of an OU process, assuming that all of its other parameters, such as evolutionary rate and selection strength, are global for the entire tree. We provide an implementation of this likelihood calculation algorithm in an accompanying R-package called PCMBase. The package has been designed as a generic library that can be integrated with existing or novel maximum likelihood or Bayesian inference tools., (Copyright © 2019 The Author(s). Published by Elsevier Inc. All rights reserved.)

Control of established chronic lymphocytic choriomeningitis virus (LCMV) infection requires the production of neutralizing antibodies, but it remains unknown how the ensemble of antibodies evolves during ongoing infection. Here, we analyze the evolution of antibody responses during acute or chronic LCMV infection, combining quantitative functional assays and time-resolved antibody repertoire sequencing. We establish that antibody responses initially converge in both infection types on a functional and repertoire level, but diverge later during chronic infection, showing increased clonal diversity, the appearance of mouse-specific persistent clones, and distinct phylogenetic signatures. Chronic infection is characterized by a longer-lasting germinal center reaction and a continuous differentiation of plasma cells, resulting in the emergence of higher-affinity plasma cells exhibiting increased antibody secretion rates. Taken together, our findings reveal the emergence of a personalized antibody response in chronic infection and support the concept that maintaining B cell diversity throughout chronic LCMV infection correlates with the development of infection-resolving antibodies., Competing Interests: Declaration of Interests The authors declare no competing interests., (Copyright © 2019 The Author(s). Published by Elsevier Inc. All rights reserved.)

Stochastic birth-death models provide the foundation for studying and simulating evolutionary trees in phylodynamics. A curious feature of such models is that they exhibit fundamental symmetries when the birth and death rates are interchanged. In this article, we first provide intuitive reasons for these known transformational symmetries. We then show that these transformational symmetries (encoded in algebraic identities) are preserved even when individuals at the present are sampled with some probability. However, these extended symmetries require the death rate parameter to sometimes take a negative value. In the last part of this article, we describe the relevance of these transformations and their application to computational phylodynamics, particularly to maximum likelihood and Bayesian inference methods, as well as to model selection., (© The Author(s) 2019. Published by Oxford University Press on behalf of the Society of Systematic Biologists.)

Introduction: Urban transmission patterns of influenza viruses are complex and poorly understood, and multiple factors may play a critical role in modifying transmission. Whole genome sequencing (WGS) allows the description of patient-to-patient transmissions at highest resolution. The aim of this study is to explore urban transmission patterns of influenza viruses in high detail by combining geographical, epidemiological and immunological data with WGS data., Methods and Analysis: The study is performed at the University Hospital Basel, University Children's Hospital Basel and a network of paediatricians and family doctors in the Canton of Basel-City, Switzerland. The retrospective study part includes an analysis of PCR-confirmed influenza cases from 2013 to 2018. The prospective study parts include (1) a household survey regarding influenza-like illness (ILI) and vaccination against influenza during the 2015/2016 season; (2) an analysis of influenza viruses collected during the 2016/2017 season using WGS-viral genomic sequences are compared with determine genetic relatedness and transmissions; and (3) measurement of influenza-specific antibody titres against all vaccinated and circulated strains during the 2016/2017 season from healthy individuals, allowing to monitor herd immunity across urban quarters. Survey data and PCR-confirmed cases are linked to data from the Statistics Office of the Canton Basel-City and visualised using geo-information system mapping. WGS data will be analysed in the context of patient epidemiological data using phylodynamic analyses, and the obtained herd immunity for each quarter. Profound knowledge on the key geographical, epidemiological and immunological factors influencing urban influenza transmission will help to develop effective counter measurements., Ethics and Dissemination: The study is registered and approved by the regional ethics committee as an observational study (EKNZ project ID 2015-363 and 2016-01735). It is planned to present the results at conferences and publish the data in scientific journals., Trial Registration Number: NCT03010007., Competing Interests: Competing interests: None declared., (© Author(s) (or their employer(s)) 2019. Re-use permitted under CC BY-NC. No commercial re-use. See rights and permissions. Published by BMJ.)

Phylogenetic comparative methods are widely used to understand and quantify the evolution of phenotypic traits, based on phylogenetic trees and trait measurements of extant species. Such analyses depend crucially on the underlying model. Gaussian phylogenetic models like Brownian motion and Ornstein-Uhlenbeck processes are the workhorses of modeling continuous-trait evolution. However, these models fit poorly to big trees, because they neglect the heterogeneity of the evolutionary process in different lineages of the tree. Previous works have addressed this issue by introducing shifts in the evolutionary model occurring at inferred points in the tree. However, for computational reasons, in all current implementations, these shifts are "intramodel," meaning that they allow jumps in 1 or 2 model parameters, keeping all other parameters "global" for the entire tree. There is no biological reason to restrict a shift to a single model parameter or, even, to a single type of model. Mixed Gaussian phylogenetic models (MGPMs) incorporate the idea of jointly inferring different types of Gaussian models associated with different parts of the tree. Here, we propose an approximate maximum-likelihood method for fitting MGPMs to comparative data comprising possibly incomplete measurements for several traits from extant and extinct phylogenetically linked species. We applied the method to the largest published tree of mammal species with body- and brain-mass measurements, showing strong statistical support for an MGPM with 12 distinct evolutionary regimes. Based on this result, we state a hypothesis for the evolution of the brain-body-mass allometry over the past 160 million y., Competing Interests: The authors declare no conflict of interest., (Copyright © 2019 the Author(s). Published by PNAS.)

Beneficial and deleterious mutations cause the fitness of lineages to vary across a phylogeny and thereby shape its branching structure. While standard phylogenetic models do not allow mutations to feedback and shape trees, birth-death models can account for this feedback by letting the fitness of lineages depend on their type. To date, these multi-type birth-death models have only been applied to cases where a lineage's fitness is determined by a single character state. We extend these models to track sequence evolution at multiple sites. This approach remains computationally tractable by tracking the genotype and fitness of lineages probabilistically in an approximate manner. Although approximate, we show that we can accurately estimate the fitness of lineages and site-specific mutational fitness effects from phylogenies. We apply this approach to estimate the population-level fitness effects of mutations in Ebola and influenza virus, and compare our estimates with in vitro fitness measurements for these mutations., Competing Interests: DR, TS No competing interests declared, (© 2019, Rasmussen and Stadler.)

Population dynamics can be inferred from genetic sequence data by using phylodynamic methods. These methods typically quantify the dynamics in unstructured populations or assume migration rates and effective population sizes to be constant through time in structured populations. When considering rates to vary through time in structured populations, the number of parameters to infer increases rapidly and the available data might not be sufficient to inform these. Additionally, it is often of interest to know what predicts these parameters rather than knowing the parameters themselves. Here, we introduce a method to  infer the predictors for time-varying migration rates and effective population sizes by using a generalized linear model (GLM) approach under the marginal approximation of the structured coalescent. Using simulations, we show that our approach is able to reliably infer the model parameters and its predictors from phylogenetic trees. Furthermore, when simulating trees under the structured coalescent, we show that our new approach outperforms the discrete trait GLM model. We then apply our framework to a previously described Ebola virus dataset, where we infer the parameters and its predictors from genome sequences while accounting for phylogenetic uncertainty. We infer weekly cases to be the strongest predictor for effective population size and geographic distance the strongest predictor for migration. This approach is implemented as part of the BEAST2 package MASCOT, which allows us to jointly infer population dynamics, i.e. the parameters and predictors, within structured populations, the phylogenetic tree, and evolutionary parameters.

Modern phylodynamic methods interpret an inferred phylogenetic tree as a partial transmission chain providing information about the dynamic process of transmission and removal (where removal may be due to recovery, death, or behavior change). Birth-death and coalescent processes have been introduced to model the stochastic dynamics of epidemic spread under common epidemiological models such as the SIS and SIR models and are successfully used to infer phylogenetic trees together with transmission (birth) and removal (death) rates. These methods either integrate analytically over past incidence and prevalence to infer rate parameters, and thus cannot explicitly infer past incidence or prevalence, or allow such inference only in the coalescent limit of large population size. Here, we introduce a particle filtering framework to explicitly infer prevalence and incidence trajectories along with phylogenies and epidemiological model parameters from genomic sequences and case count data in a manner consistent with the underlying birth-death model. After demonstrating the accuracy of this method on simulated data, we use it to assess the prevalence through time of the early 2014 Ebola outbreak in Sierra Leone., (© The Author(s) 2019. Published by Oxford University Press on behalf of the Society for Molecular Biology and Evolution.)

Mass extinction events (MEEs), defined as significant losses of species diversity in significantly short time periods, have attracted the attention of biologists because of their link to major environmental change. MEEs have traditionally been studied through the fossil record, but the development of birth-death models has made it possible to detect their signature based on extant-taxa phylogenies. Most birth-death models consider MEEs as instantaneous events where a high proportion of species are simultaneously removed from the tree ("single pulse" approach), in contrast to the paleontological record, where MEEs have a time duration. Here, we explore the power of a Bayesian Birth-Death Skyline (BDSKY) model to detect the signature of MEEs through changes in extinction rates under a "time-slice" approach. In this approach, MEEs are time intervals where the extinction rate is greater than the speciation rate. Results showed BDSKY can detect and locate MEEs but that precision and accuracy depend on the phylogeny's size and MEE intensity. Comparisons of BDSKY with the single-pulse Bayesian model, CoMET, showed a similar frequency of Type II error and neither model exhibited Type I error. However, while CoMET performed better in detecting and locating MEEs for smaller phylogenies, BDSKY showed higher accuracy in estimating extinction and speciation rates., (© 2019 The Author(s). Evolution published by Wiley Periodicals, Inc. on behalf of The Society for the Study of Evolution.)

Fossil information is essential for estimating species divergence times, and can be integrated into Bayesian phylogenetic inference using the fossilized birth-death (FBD) process. An important aspect of palaeontological data is the uncertainty surrounding specimen ages, which can be handled in different ways during inference. The most common approach is to fix fossil ages to a point estimate within the known age interval. Alternatively, age uncertainty can be incorporated by using priors, and fossil ages are then directly sampled as part of the inference. This study presents a comparison of alternative approaches for handling fossil age uncertainty in analysis using the FBD process. Based on simulations, we find that fixing fossil ages to the midpoint or a random point drawn from within the stratigraphic age range leads to biases in divergence time estimates, while sampling fossil ages leads to estimates that are similar to inferences that employ the correct ages of fossils. Second, we show a comparison using an empirical dataset of extant and fossil cetaceans, which confirms that different methods of handling fossil age uncertainty lead to large differences in estimated node ages. Stratigraphic age uncertainty should thus not be ignored in divergence time estimation and instead should be incorporated explicitly.

Asymmetric cell division is a major mechanism generating cell diversity. As cell cycle duration varies among cells in mammalian tissue culture cells, we asked whether their division asymmetry contributes to this variability. We identify among sibling cells an outlier using hierarchical clustering on cell cycle durations of granddaughter cells obtained by lineage tracking of single histone2B-labelled MDCKs. Remarkably, divisions involving outlier cells are not uniformly distributed in lineages, as shown by permutation tests, but appear to emerge from asymmetric divisions taking place at non-stochastic levels: a parent cell influences with 95% confidence and 0.5% error the unequal partitioning of the cell cycle duration in its two progenies. Upon ninein downregulation, this variability propagation is lost, and outlier frequency and variability in cell cycle durations in lineages is reduced. As external influences are not detectable, we propose that a cell-autonomous process, possibly involved in cell specialisation, determines cell cycle duration variability.

Elaboration of Bayesian phylogenetic inference methods has continued at pace in recent years with major new advances in nearly all aspects of the joint modelling of evolutionary data. It is increasingly appreciated that some evolutionary questions can only be adequately answered by combining evidence from multiple independent sources of data, including genome sequences, sampling dates, phenotypic data, radiocarbon dates, fossil occurrences, and biogeographic range information among others. Including all relevant data into a single joint model is very challenging both conceptually and computationally. Advanced computational software packages that allow robust development of compatible (sub-)models which can be composed into a full model hierarchy have played a key role in these developments. Developing such software frameworks is increasingly a major scientific activity in its own right, and comes with specific challenges, from practical software design, development and engineering challenges to statistical and conceptual modelling challenges. BEAST 2 is one such computational software platform, and was first announced over 4 years ago. Here we describe a series of major new developments in the BEAST 2 core platform and model hierarchy that have occurred since the first release of the software, culminating in the recent 2.5 release., Competing Interests: The authors have declared that no competing interests exist.

We compared two sample preparation protocols for whole genome sequencing of influenza A viruses. Each protocol was assessed using cDNA quantity and quality and the resulting mean genome coverage after sequencing. Both protocols produced acceptable result for samples with high viral load, whereas one protocol performed slightly better with limited virus count., (Copyright © 2019 The Authors. Published by Elsevier B.V. All rights reserved.)

Human immunodeficiency virus type 1 (HIV-1) non-B subtype infections occurred in Belgium since the 1980s, mainly amongst migrants and heterosexuals, whereas subtype B predominated in men-having-sex-with-men (MSM). In the last decade, the diagnosis of F1 sub-subtype in particular has increased substantially, which prompted us to perform a detailed reconstruction of its epidemiological history. To this purpose, the Belgian AIDS Reference Laboratories collected HIV-1 pol sequences from all sub-subtype F1-infected patients for whom genotypic drug resistance testing was requested as part of routine clinical follow-up. This data was complemented with HIV-1 pol sequences from countries with a high burden of F1 infections or a potential role in the global origin of sub-subtype F1. The molecular epidemiology of the Belgian subtype F1 epidemic was investigated using Bayesian phylogenetic inference and transmission dynamics were characterized based on birth-death models. F1 sequences were retained from 297 patients diagnosed and linked to care in Belgium between 1988 and 2015. Phylogenetic inference indicated that among the 297 Belgian F1 sequences, 191 belonged to a monophyletic group that mainly contained sequences from people likely infected in Belgium (OR 26.67, 95% CI 9.59-74.15), diagnosed in Flanders (OR 7.28, 95% CI 4.23-12.53), diagnosed at a recent stage of infection (OR 7.19, 95% CI 2.88-17.95) or declared to be MSM (OR 34.8, 95% CI 16.0-75.6). Together with a Spanish clade, this Belgian clade was embedded in the genetic diversity of Brazilian subtype F1 strains and most probably emerged after one or only a few migration events from Brazil to the European continent before 2002. The origin of the Belgian outbreak was dated back to 2002 (95% higher posterior density 2000-2004) and birth-death models suggested that its extensive growth had been controlled ( R e < 1) by 2012, coinciding with a time period where delay in antiretroviral treatment initiation substantially declined. In conclusion, phylogenetic reconstruction of the Belgian HIV-1 sub-subtype F1 epidemic illustrates the introduction and substantial dissemination of viral strains in a geographically restricted risk group that was most likely controlled by effective treatment as prevention.

Rapidly evolving pathogens, such as viruses and bacteria, accumulate genetic change at a similar timescale over which their epidemiological processes occur, such that, it is possible to make inferences about their infectious spread using phylogenetic time-trees. For this purpose it is necessary to choose a phylodynamic model. However, the resulting inferences are contingent on whether the model adequately describes key features of the data. Model adequacy methods allow formal rejection of a model if it cannot generate the main features of the data. We present TreeModelAdequacy, a package for the popular BEAST2 software that allows assessing the adequacy of phylodynamic models. We illustrate its utility by analyzing phylogenetic trees from two viral outbreaks of Ebola and H1N1 influenza. The main features of the Ebola data were adequately described by the coalescent exponential-growth model, whereas the H1N1 influenza data were best described by the birth-death susceptible-infected-recovered model.

Despite increasing access to antiretrovirals, HIV incidence in rural KwaZulu-Natal remains among the highest ever reported in Africa. While many epidemiological factors have been invoked to explain such high incidence, widespread human mobility and viral movement suggest that transmission between communities may be a major source of new infections. High cross-community transmission rates call into question how effective increasing the coverage of antiretroviral therapy locally will be at preventing new infections, especially if many new cases arise from external introductions. To help address this question, we use a phylodynamic model to reconstruct epidemic dynamics and estimate the relative contribution of local transmission versus external introductions to overall incidence in KwaZulu-Natal from HIV-1 phylogenies. By comparing our results with population-based surveillance data, we show that we can reliably estimate incidence from viral phylogenies once viral movement in and out of the local population is accounted for. Our analysis reveals that early epidemic dynamics were largely driven by external introductions. More recently, we estimate that 35 per cent (95% confidence interval: 20-60%) of new infections arise from external introductions. These results highlight the growing need to consider larger-scale regional transmission dynamics when designing and testing prevention strategies.

Measuring the pace at which speciation and extinction occur is fundamental to understanding the origin and evolution of biodiversity. Both the fossil record and molecular phylogenies of living species can provide independent estimates of speciation and extinction rates, but often produce strikingly divergent results. Despite its implications, the theoretical reasons for this discrepancy remain unknown. Here, we reveal a conceptual and methodological basis able to reconcile palaeontological and molecular evidence: discrepancies are driven by different implicit assumptions about the processes of speciation and species evolution in palaeontological and neontological analyses. We present the "birth-death chronospecies" model that clarifies the definition of speciation and extinction processes allowing for a coherent joint analysis of fossil and phylogenetic data. Using simulations and empirical analyses we demonstrate not only that this model explains much of the apparent incongruence between fossils and phylogenies, but that differences in rate estimates are actually informative about the prevalence of different speciation modes.

The fast evolution of pathogenic viruses has allowed for the development of phylodynamic approaches that extract information about the epidemiological characteristics of viral genomes. Thanks to advances in whole genome sequencing, they can be applied to slowly evolving bacterial pathogens like Mycobacterium tuberculosis. In this study, we investigate and compare the epidemiological dynamics underlying two M. tuberculosis outbreaks using phylodynamic methods. Specifically, we (i) test if the outbreak data sets contain enough genetic variation to estimate short-term evolutionary rates and (ii) reconstruct epidemiological parameters such as the effective reproduction number. The first outbreak occurred in the Swiss city of Bern (1987-2012) and was caused by a drug-susceptible strain belonging to the phylogenetic M. tuberculosis Lineage 4. The second outbreak was caused by a multidrug-resistant (MDR) strain of Lineage 2, imported from the Wat Tham Krabok (WTK) refugee camp in Thailand into California. There is little temporal signal in the Bern data set and moderate temporal signal in the WTK data set. Thanks to its high sampling proportion (90%) the Bern outbreak allows robust estimation of epidemiological parameters despite the poor temporal signal. Conversely, there is much uncertainty in the epidemiological estimates concerning the sparsely sampled (9%) WTK outbreak. Our results suggest that both outbreaks peaked around 1990, although they were only recognized as outbreaks in 1993 (Bern) and 2004 (WTK). Furthermore, individuals were infected for a significantly longer period (around 9 years) in the WTK outbreak than in the Bern outbreak (4-5 years). Our work highlights both the limitations and opportunities of phylodynamic analysis of outbreaks involving slowly evolving pathogens: (i) estimation of the evolutionary rate is difficult on outbreak time scales and (ii) a high sampling proportion allows quantification of the age of the outbreak based on the sampling times, and thus allows for robust estimation of epidemiological parameters., (Copyright © 2018 The Authors. Published by Elsevier B.V. All rights reserved.)

Motivation: The structured coalescent is widely applied to study demography within and migration between sub-populations from genetic sequence data. Current methods are either exact but too computationally inefficient to analyse large datasets with many sub-populations, or make strong approximations leading to severe biases in inference. We recently introduced an approximation based on weaker assumptions to the structured coalescent enabling the analysis of larger datasets with many different states. We showed that our approximation provides unbiased migration rate and population size estimates across a wide parameter range., Results: We extend this approach by providing a new algorithm to calculate the probability of the state of internal nodes that includes the information from the full phylogenetic tree. We show that this algorithm is able to increase the probability attributed to the true sub-population of a node. Furthermore we use improved integration techniques, such that our method is now able to analyse larger datasets, including a H3N2 dataset with 433 sequences sampled from five different locations., Availability and Implementation: The presented methods are part of the BEAST2 package MASCOT, the Marginal Approximation of the Structured COalescenT. This package can be downloaded via the BEAUti package manager. The source code is available at https://github.com/nicfel/Mascot.git., Supplementary Information: Supplementary data are available at Bioinformatics online.

Background: Tracking recent transmission is a vital part of controlling widespread pathogens such as Mycobacterium tuberculosis. Multiple methods with specific performance characteristics exist for detecting recent transmission chains, usually by clustering strains based on genotype similarities. With such a large variety of methods available, informed selection of an appropriate approach for determining transmissions within a given setting/time period is difficult., Methods: This study combines whole genome sequence (WGS) data derived from 324 isolates collected 2005-2010 in Kinshasa, Democratic Republic of Congo (DRC), a high endemic setting, with phylodynamics to unveil the timing of transmission events posited by a variety of standard genotyping methods. Clustering data based on Spoligotyping, 24-loci MIRU-VNTR typing, WGS based SNP (Single Nucleotide Polymorphism) and core genome multi locus sequence typing (cgMLST) typing were evaluated., Findings: Our results suggest that clusters based on Spoligotyping could encompass transmission events that occurred almost 200 years prior to sampling while 24-loci-MIRU-VNTR often represented three decades of transmission. Instead, WGS based genotyping applying low SNP or cgMLST allele thresholds allows for determination of recent transmission events, e.g. in timespans of up to 10 years for a 5 SNP/allele cut-off., Interpretation: With the rapid uptake of WGS methods in surveillance and outbreak tracking, the findings obtained in this study can guide the selection of appropriate clustering methods for uncovering relevant transmission chains within a given time-period. For high resolution cluster analyses, WGS-SNP and cgMLST based analyses have similar clustering/timing characteristics even for data obtained from a high incidence setting., (Copyright © 2018 The Authors. Published by Elsevier B.V. All rights reserved.)

Antibody evolution studies have been traditionally limited to either tracing a single clonal lineage (B cells derived from a single V-(D)-J recombination) over time or examining bulk functionality changes (e.g., tracing serum polyclonal antibody proteins). Studying a single B cell disregards the majority of the humoral immune response, whereas bulk functional studies lack the necessary resolution to analyze the co-existing clonal diversity. Recent advances in high-throughput sequencing (HTS) technologies and bioinformatics have made it possible to examine multiple co-evolving antibody monoclonal lineages within the context of a single repertoire. A plethora of accompanying methods and tools have been introduced in hopes of better understanding how pathogen presence dictates the global evolution of the antibody repertoire. Here, we provide a comprehensive summary of the tremendous progress of this newly emerging field of systems phylogeny of antibody responses. We present an overview encompassing the historical developments of repertoire phylogenetics, state-of-the-art tools, and an outlook on the future directions of this fast-advancing and promising field.

HIV patients form clusters in HIV transmission networks. Accurate identification of these transmission clusters is essential to effectively target public health interventions. One reason for clustering is that the underlying contact network contains many local communities. We present a new maximum-likelihood method for identifying transmission clusters caused by community structure, based on phylogenetic trees. The method employs a multi-state birth-death (MSBD) model which detects changes in transmission rate, which are interpreted as the introduction of the epidemic into a new susceptible community, i.e. the formation of a new cluster. We show that the MSBD method is able to reliably infer the clusters and the transmission parameters from a pathogen phylogeny based on our simulations. In contrast to existing cutpoint-based methods for cluster identification, our method does not require that clusters be monophyletic nor is it dependent on the selection of a difficult-to-interpret cutpoint parameter. We present an application of our method to data from the Swiss HIV Cohort Study. The method is available as an easy-to-use R package., (© 2018 The Authors.)

Most stored-grain pest insects increase their population within a relatively short time, causing serious damage to stored products. Sitophilus oryzae (L.) is one of the world's major stored-grain pest insects and was chosen as the model insect for our studies. This study compared the efficacy of three different dusts under laboratory conditions: aluminum dust (nanostructured alumina), DiatomiD ® , and Protect-It ® (commercial diatomaceous earth). Parental survival, grain damage, and progeny production were measured at 250 and 500 ppm in treated wheat. The tests were conducted in 400 mL galvanized steel jars, an experimental model used for the first time to measure the effectiveness of nanostructured alumina, since most studies have been typically performed in small petri dishes. Parental survival obtained was highest in the untreated controls, followed in decreasing order by DiatomiD ® , Protect-It ® , and nanostructured alumina (NSA). NSA caused the greatest mortality. All treatments significantly reduced grain weight loss and frass production in wheat infested by S. oryzae . The degree of progeny (F1) suppression was directly related to the product and treatment rate, progeny being significantly suppressed by NSA in wheat followed by Protect-It ® and DiatomiD ® . Therefore, NSA had a greater impact on insect population dynamics.

A birth-death-sampling model gives rise to phylogenetic trees with samples from the past and the present. Interpreting "birth" as branching speciation, "death" as extinction, and "sampling" as fossil preservation and recovery, this model - also referred to as the fossilized birth-death (FBD) model - gives rise to phylogenetic trees on extant and fossil samples. The model has been mathematically analyzed and successfully applied to a range of datasets on different taxonomic levels, such as penguins, plants, and insects. However, the current mathematical treatment of this model does not allow for a group of temporally distinct fossil specimens to be assigned to the same species. In this paper, we provide a general mathematical FBD modeling framework that explicitly takes "stratigraphic ranges" into account, with a stratigraphic range being defined as the lineage interval associated with a single species, ranging through time from the first to the last fossil appearance of the species. To assign a sequence of fossil samples in the phylogenetic tree to the same species, i.e., to specify a stratigraphic range, we need to define the mode of speciation. We provide expressions to account for three common speciation modes: budding (or asymmetric) speciation, bifurcating (or symmetric) speciation, and anagenetic speciation. Our equations allow for flexible joint Bayesian analysis of paleontological and neontological data. Furthermore, our framework is directly applicable to epidemiology, where a stratigraphic range is the observed duration of infection of a single patient, "birth" via budding is transmission, "death" is recovery, and "sampling" is sequencing the pathogen of a patient. Thus, we present a model that allows for incorporation of multiple observations through time from a single patient., (Copyright © 2018 The Authors. Published by Elsevier Ltd.. All rights reserved.)

Bayesian phylogenetics aims at estimating phylogenetic trees together with evolutionary and population dynamic parameters based on genetic sequences. It has been noted that the clock rate, one of the evolutionary parameters, decreases with an increase in the sampling period of sequences. In particular, clock rates of epidemic outbreaks are often estimated to be higher compared with the long-term clock rate. Purifying selection has been suggested as a biological factor that contributes to this phenomenon, since it purges slightly deleterious mutations from a population over time. However, other factors such as methodological biases may also play a role and make a biological interpretation of results difficult. In this paper, we identify methodological biases originating from the choice of tree prior, that is, the model specifying epidemiological dynamics. With a simulation study we demonstrate that a misspecification of the tree prior can upwardly bias the inferred clock rate and that the interplay of the different models involved in the inference can be complex and nonintuitive. We also show that the choice of tree prior can influence the inference of clock rate on real-world Ebola virus (EBOV) datasets. While commonly used tree priors result in very high clock-rate estimates for sequences from the initial phase of the epidemic in Sierra Leone, tree priors allowing for population structure lead to estimates agreeing with the long-term rate for EBOV., Competing Interests: The authors declare no conflict of interest., (Copyright © 2018 the Author(s). Published by PNAS.)

Pathogen traits, such as the virulence of an infection, can vary significantly between patients. A major challenge is to measure the extent to which genetic differences between infecting strains explain the observed variation of the trait. This is quantified by the trait's broad-sense heritability, H2. A recent discrepancy between estimates of the heritability of HIV-virulence has opened a debate on the estimators' accuracy. Here, we show that the discrepancy originates from model limitations and important lifecycle differences between sexually reproducing organisms and transmittable pathogens. In particular, current quantitative genetics methods, such as donor-recipient regression of surveyed serodiscordant couples and the phylogenetic mixed model (PMM), are prone to underestimate H2, because they neglect or do not fit to the loss of resemblance between transmission partners caused by within-host evolution. In a phylogenetic analysis of 8,483 HIV patients from the United Kingdom, we show that the phenotypic correlation between transmission partners decays with the amount of within-host evolution of the virus. We reproduce this pattern in toy-model simulations and show that a phylogenetic Ornstein-Uhlenbeck model (POUMM) outperforms the PMM in capturing this correlation pattern and in quantifying H2. In particular, we show that POUMM outperforms PMM even in simulations without selection-as it captures the mentioned correlation pattern-which has not been appreciated until now. By cross-validating the POUMM estimates with ANOVA on closest phylogenetic pairs, we obtain H2 ≈ 0.2, meaning ∼20% of the variation in HIV-virulence is explained by the virus genome both for European and African data., (© The Author(s) 2018. Published by Oxford University Press on behalf of the Society for Molecular Biology and Evolution.)

Understanding macroevolutionary processes using phylogenetic trees is a challenging and complex process that draws on mathematics, computer science and biology. Given the development of complex mathematical models and the growing computational processing power, simulation tools are becoming increasingly popular.In order to simulate phylogenetic trees, most evolutionary biologists are forced to build their own algorithms or use existing tools built on different platforms and/or as standalone programmes. The absence of a simulation tool accommodating for user-chosen model specifications limits, amongst others, model testing and pipelining with approximate Bayesian computation methods or other subsequent statistical analysis.We introduce "TreeSimGM ," an r-package simulation tool for phylogenetic trees under a general Bellman and Harris model. This package allows the user to specify any desired probability distribution for the waiting times until speciation and extinction (e.g. age-dependent speciation/extinction). Upon speciation, the user can specify whether one descendant species corresponds to the ancestor species inheriting its age or whether both descendant species are new species of age 0. Moreover, it is possible to scale the waiting time to speciation/extinction for newly formed species. Thus, "TreeSimGM" not only allows the user to simulate stochastic phylogenetic trees assuming several popular existing models, such as the Yule model, the constant-rate birth-death model, and proportional to distinguishable arrangement models, but it also allows the user to formulate new models for exploration. A short explanation of the supported models and a few examples of how to use our package are presented here.As an r-package, "TreeSimGM" allows flexible and powerful stochastic phylogenetic tree simulations. Moreover, it facilitates the pipelining of outputs or inputs with other functions in r. "TreeSimGM" contributes to the tools available to the r community in the fields of ecology and evolution, is freely available under the GPL-2 licence and can be downloaded at https://cran.r-project.org/web/packages/TreeSimGM.

Drug resistant HIV is a major threat to the long-term efficacy of antiretroviral treatment. Around 10% of ART-naïve patients in Europe are infected with drug-resistant HIV type 1. Hence it is important to understand the dynamics of transmitted drug resistance evolution. Thanks to routinely performed drug resistance tests, HIV sequence data is increasingly available and can be used to reconstruct the phylogenetic relationship among viral lineages. In this study we employ a phylodynamic approach to quantify the fitness costs of major resistance mutations in the Swiss HIV cohort. The viral phylogeny reflects the transmission tree, which we model using stochastic birth-death-sampling processes with two types: hosts infected by a sensitive or resistant strain. This allows quantification of fitness cost as the ratio between transmission rates of hosts infected by drug resistant strains and transmission rates of hosts infected by drug sensitive strains. The resistance mutations 41L, 67N, 70R, 184V, 210W, 215D, 215S and 219Q (nRTI-related) and 103N, 108I, 138A, 181C, 190A (NNRTI-related) in the reverse trancriptase and the 90M mutation in the protease gene are included in this study. Among the considered resistance mutations, only the 90M mutation in the protease gene was found to have significantly higher fitness than the drug sensitive strains. The following mutations associated with resistance to reverse transcriptase inhibitors were found to be less fit than the sensitive strains: 67N, 70R, 184V, 219Q. The highest posterior density intervals of the transmission ratios for the remaining resistance mutations included in this study all included 1, suggesting that these mutations do not have a significant effect on viral transmissibility within the Swiss HIV cohort. These patterns are consistent with alternative measures of the fitness cost of resistance mutations. Overall, we have developed and validated a novel phylodynamic approach to estimate the transmission fitness cost of drug resistance mutations.

Introduction: Extended-spectrum beta-lactamases (ESBL)-producing Enterobacteriaceae were first described in relation with hospital-acquired infections. In the 2000s, the epidemiology of ESBL-producing organisms changed as especially ESBL-producing Escherichia coli was increasingly described as an important cause of community-acquired infections, supporting the hypothesis that in more recent years ESBL-producing Enterobacteriaceae have probably been imported into hospitals rather than vice versa. Transmission of ESBL-producing Enterobacteriaceae is complicated by ESBL genes being encoded on self-transmissible plasmids, which can be exchanged among the same and different bacterial species. The aim of this research project is to quantify hospital-wide transmission of ESBL-producing Enterobacteriaceae on both the level of bacterial species and the mobile genetic elements and to determine if hospital-acquired infections caused by ESBL producers are related to strains and mobile genetic elements predominantly circulating in the community or in the healthcare setting. This distinction is critical in prevention since the former emphasises the urgent need to establish or reinforce antibiotic stewardship programmes, and the latter would call for more rigorous infection control., Methods and Analysis: This protocol presents an observational study that will be performed at the University Hospital Basel and in the city of Basel, Switzerland. ESBL-producing Enterobacteriaceae will be collected from any specimens obtained by routine clinical practice or by active screening in both inpatient and outpatient settings, as well as from wastewater samples and foodstuffs, both collected monthly over a 12-month period for analyses by whole genome sequencing. Bacterial chromosomal, plasmid and ESBL-gene sequences will be compared within the cohort to determine genetic relatedness and migration between humans and their environment., Ethics and Dissemination: This study has been approved by the local ethics committee (Ethikkommission Nordwest-und Zentralschweiz) as a quality control project (Project-ID 2017-00100). The results of this study will be published in peer-reviewed medical journals, communicated to participants, the general public and all relevant stakeholders., Competing Interests: Competing interests: ST-S is a member of the Astellas Fidaxomicin Advisory Board and a member of the MSD C. difficile Advisory Board. She has received unrestricted research grants from Astellas., (© Article author(s) (or their employer(s) unless otherwise stated in the text of the article) 2018. All rights reserved. No commercial use is permitted unless otherwise expressly granted.)

Stem cells play a central role in the regeneration and repair of multicellular organisms. However, it remains far from trivial to reliably identify them. Despite decades of work, current techniques to isolate hematopoietic stem cells (HSCs) based on cell-surface markers only result in 50% purity, i.e. half of the sorted cells are not stem cells when functionally tested. Modern microscopy techniques allow us to follow single cells and their progeny for up to weeks in vitro, while recording the cell fates and lifetime of each individual cell. This cell tracking generates so-called lineage trees. Here, we propose statistical techniques to determine if the initial cell in a lineage tree was a HSC. We apply these techniques to murine hematopoietic lineage trees, revealing that 18% of the trees in our HSC dataset display a unique signature, and this signature is compatible with these trees having started from a true stem cell. Assuming 50% purity of HSC empirical datasets, this corresponds to a 0.35 power of the test, and the type-1-error is estimated to be 0.047. In summary, this study shows that statistical analysis of lineage trees could improve the classification of cells, which is currently done based on bio-markers only. Our statistical techniques are not limited to mammalian stem cell biology. Any type of single cell lineage trees, be it from bacteria, single cell eukaryotes, or single cells in a multicellular organism can be investigated. We expect this to contribute to a better understanding of the molecules influencing cellular dynamics at the single cell level., (Copyright © 2017 The Authors. Published by Elsevier Ltd.. All rights reserved.)

Reticulate species evolution, such as hybridization or introgression, is relatively common in nature. In the presence of reticulation, species relationships can be captured by a rooted phylogenetic network, and orthologous gene evolution can be modeled as bifurcating gene trees embedded in the species network. We present a Bayesian approach to jointly infer species networks and gene trees from multilocus sequence data. A novel birth-hybridization process is used as the prior for the species network, and we assume a multispecies network coalescent prior for the embedded gene trees. We verify the ability of our method to correctly sample from the posterior distribution, and thus to infer a species network, through simulations. To quantify the power of our method, we reanalyze two large data sets of genes from spruces and yeasts. For the three closely related spruces, we verify the previously suggested homoploid hybridization event in this clade; for the yeast data, we find extensive hybridization events. Our method is available within the BEAST 2 add-on SpeciesNetwork, and thus provides an extensible framework for Bayesian inference of reticulate evolution., (© The Author(s) 2017. Published by Oxford University Press on behalf of the Society for Molecular Biology and Evolution.)