A Monte Carlo method to estimate cell population heterogeneity from cell snapshot data.

Variation is characteristic of all living systems. Laboratory techniques such as flow cytometry can probe individual cells, and, after decades of experimentation, it is clear that even members of genetically identical cell populations can exhibit differences. To understand whether variation is biologically meaningful, it is essential to discern its source. Mathematical models of biological systems are tools that can be used to investigate causes of cell-to-cell variation. From mathematical analysis and simulation of these models, biological hypotheses can be posed and investigated, then parameter inference can determine which of these is compatible with experimental data. Data from laboratory experiments often consist of "snapshots" representing distributions of cellular properties at different points in time, rather than individual cell trajectories. These data are not straightforward to fit using hierarchical Bayesian methods, which require the number of cell population clusters to be chosen a priori. Nor are they amenable to standard nonlinear mixed effect methods, since a single observation per cell is typically too few to estimate parameter variability. Here, we introduce a computational sampling method named "Contour Monte Carlo" (CMC) for estimating mathematical model parameters from snapshot distributions, which is straightforward to implement and does not require that cells be assigned to predefined categories. The CMC algorithm fits to snapshot probability distributions rather than raw data, which means its computational burden does not, like existing approaches, increase with the number of cells observed. Our method is appropriate for underdetermined systems, where there are fewer distinct types of observations than parameters to be determined, and where observed variation is mostly due to variability in cellular processes rather than experimental measurement error. This may be the case for many systems due to continued improvements in resolution of laboratory techniques. In this paper, we apply our method to quantify cellular variation for three biological systems of interest and provide Julia code enabling others to use this method.

The power of evolutionary rescue is constrained by genetic load.

The risk of extinction faced by small isolated populations in changing environments can be reduced by rapid adaptation and subsequent growth to larger, less vulnerable sizes. Whether this process, called evolutionary rescue, is able to reduce extinction risk and sustain population growth over multiple generations is largely unknown. To understand the consequences of adaptive evolution as well as maladaptive processes in small isolated populations, we subjected experimental Tribolium castaneum populations founded with 10 or 40 individuals to novel environments, one more favorable, and one resource poor, and either allowed evolution, or constrained it by replacing individuals one-for-one each generation with those from a large population maintained in the natal environment. Replacement individuals spent one generation in the target novel environment before use to standardize effects due to the parental environment. After eight generations we mixed a subset of surviving populations to facilitate admixture, allowing us to estimate drift load by comparing performance of mixed to unmixed groups. Evolving populations had reduced extinction rates, and increased population sizes in the first four to five generations compared to populations where evolution was constrained. Performance of evolving populations subsequently declined. Admixture restored their performance, indicating high drift load that may have overwhelmed the beneficial effects of adaptation in evolving populations. Our results indicate that evolution may quickly reduce extinction risk and increase population sizes, but suggest that relying solely on adaptation from standing genetic variation may not provide long-term benefits to small isolated populations of diploid sexual species, and that active management facilitating gene flow may be necessary for longer term persistence.

Dynamics of prion disease transmission in mule deer.

Chronic wasting disease (CWD), a contagious prion disease of the deer family, has the potential to severely harm deer populations and disrupt ecosystems where deer occur in abundance. Consequently, understanding the dynamics of this emerging infectious disease, and particularly the dynamics of its transmission, has emerged as an important challenge for contemporary ecologists and wildlife managers. Although CWD is contagious among deer, the relative importance of pathways for its transmission remains unclear. We developed seven competing models, and then used data from two CWD outbreaks in captive mule deer and model selection to compare them. We found that models portraying indirect transmission through the environment had 3.8 times more support in the data than models representing transmission by direct contact between infected and susceptible deer. Model-averaged estimates of the basic reproductive number (R0) were 1.3 or greater, indicating likely local persistence of CWD in natural populations under conditions resembling those we studied. Our findings demonstrate the apparent importance of indirect, environmental transmission in CWD and the challenges this presents for controlling the disease.