Optimal immune specificity at the intersection of host life history and parasite epidemiology

Hosts diverge widely in how, and how well, they defend themselves against infection and immunopathology. Why are hosts so heterogeneous? Both epidemiology and life history are commonly hypothesized to influence host immune strategy, but the relationship between immune strategy and each factor has commonly been investigated in isolation. Here, we show that interactions between life history and epidemiology are crucial for determining optimal immune specificity and sensitivity. We propose a demographically-structured population dynamics model, in which we explore sensitivity and specificity of immune responses when epidemiological risks vary with age. We find that variation in life history traits associated with both reproduction and longevity alters optimal immune strategies–but the magnitude and sometimes even direction of these effects depends on how epidemiological risks vary across life. An especially compelling example that explains previously-puzzling empirical observations is that depending on whether infection risk declines or rises at reproductive maturity, later reproductive maturity can select for either greater or lower immune specificity, potentially illustrating why studies of lifespan and immune variation across taxa have been inconclusive. Thus, the sign of selection on the life history-immune specificity relationship can be reversed in different epidemiological contexts. Drawing on published life history data from a variety of chordate taxa, we generate testable predictions for this facet of the optimal immune strategy. Our results shed light on the causes of the heterogeneity found in immune defenses both within and among species and the ultimate variability of the relationship between life history and immune specificity.

Characterizing the Observation Bias in Gravitational-wave Detections and Finding Structured Population Properties

The observed distributions of the source properties from gravitational-wave (GW) detections are biased due to the selection effects and detection criteria in the detections, analogous to the Malmquist bias. In this work, this observation bias is investigated through its fundamental statistical and physical origins. An efficient semi-analytical formulation for its estimation is derived, which is as accurate as the standard method of numerical simulations, with only a millionth of the computational cost. Then, the estimated bias is used for unmodeled inferences on the binary black hole population. These inferences show additional structures, specifically two peaks in the joint mass distribution around binary masses ∼10 M ⊙ and ∼30 M ⊙. Example ready-to-use scripts and some produced data sets for this method are shared in an online repository.

Dynamic maximum entropy provides accurate approximation of structured population dynamics

Realistic models of biological processes typically involve interacting components on multiple scales, driven by changing environment and inherent stochasticity. Such models are often analytically and numerically intractable. We revisit a dynamic maximum entropy method that combines a static maximum entropy with a quasi-stationary approximation. This allows us to reduce stochastic non-equilibrium dynamics expressed by the Fokker-Planck equation to a simpler low-dimensional deterministic dynamics, without the need to track microscopic details. Although the method has been previously applied to a few (rather complicated) applications in population genetics, our main goal here is to explain and to better understand how the method works. We demonstrate the usefulness of the method for two widely studied stochastic problems, highlighting its accuracy in capturing important macroscopic quantities even in rapidly changing non-stationary conditions. For the Ornstein-Uhlenbeck process, the method recovers the exact dynamics whilst for a stochastic island model with migration from other habitats, the approximation retains high macroscopic accuracy under a wide range of scenarios in a dynamic environment.

Dynamics of a new stage-structured population model with transient and nontransient impulsive effects in a polluted environment

In this paper, we consider a new stage-structured population model with transient and nontransient impulsive effects in a polluted environment. By using the theories of impulsive differential equations, we obtain the globally asymptotically stable condition of a population-extinction solution; we also present the permanent condition for the investigated system. The results indicate that the nontransient and transient impulsive harvesting rate play important roles in system permanence. Finally, numerical analyses are carried out to illustrate the results. Our results provide effective methods for biological resource management in a polluted environment.

Employing phylogenetic tree shape statistics to resolve the underlying host population structure

Background Host population structure is a key determinant of pathogen and infectious disease transmission patterns. Pathogen phylogenetic trees are useful tools to reveal the population structure underlying an epidemic. Determining whether a population is structured or not is useful in informing the type of phylogenetic methods to be used in a given study. We employ tree statistics derived from phylogenetic trees and machine learning classification techniques to reveal an underlying population structure. Results In this paper, we simulate phylogenetic trees from both structured and non-structured host populations. We compute eight statistics for the simulated trees, which are: the number of cherries; Sackin, Colless and total cophenetic indices; ladder length; maximum depth; maximum width, and width-to-depth ratio. Based on the estimated tree statistics, we classify the simulated trees as from either a non-structured or a structured population using the decision tree (DT), K-nearest neighbor (KNN) and support vector machine (SVM). We incorporate the basic reproductive number ($$R_0$$ R 0 ) in our tree simulation procedure. Sensitivity analysis is done to investigate whether the classifiers are robust to different choice of model parameters and to size of trees. Cross-validated results for area under the curve (AUC) for receiver operating characteristic (ROC) curves yield mean values of over 0.9 for most of the classification models. Conclusions Our classification procedure distinguishes well between trees from structured and non-structured populations using the classifiers, the two-sample Kolmogorov-Smirnov, Cucconi and Podgor-Gastwirth tests and the box plots. SVM models were more robust to changes in model parameters and tree size compared to KNN and DT classifiers. Our classification procedure was applied to real -world data and the structured population was revealed with high accuracy of $$92.3\%$$ 92.3 % using SVM-polynomial classifier.