Identifying opportunities for high resolution pesticide usage data to improve the efficiency of endangered species pesticide risk assessment.

The US pesticide registration and review process requires regular re-assessment of the risk of pesticide use to species listed under the Endangered Species Act (ESA), yet current assessment methods are inefficient when applied to hundreds of pesticides potentially impacting multiple species across a continent. Thus, many pesticides remain on the market without complete review. We assessed the value of using high resolution pesticide usage data in the risk assessment process to rapidly improve process efficiency. By using data available only in California, we found that high resolution data increased the number of species deemed not likely to be adversely affected by pesticides from <5 % to nearly 50 %. Across the contiguous US, we predicted that 48 % of species would be deemed not likely to be adversely affected using high resolution data, compared to 20 % without. However, if such data were available in just 11 states, 68 % of the available gains in efficiency could be obtained. Overall, using existing high-resolution data in California and a focused collection of such information from 11 other states could reduce risk assessment burden across the contiguous U.S. by one-quarter.

Tracking glioblastoma progression after initial resection with minimal reaction-diffusion models.

We describe a preliminary effort to model the growth and progression of glioblastoma multiforme, an aggressive form of primary brain cancer, in patients undergoing treatment for recurrence of tumor following initial surgery and chemoradiation. Two reaction-diffusion models are used: the Fisher-Kolmogorov equation and a 2-population model, developed by the authors, that divides the tumor into actively proliferating and quiescent (or necrotic) cells. The models are simulated on 3-dimensional brain geometries derived from magnetic resonance imaging (MRI) scans provided by the Barrow Neurological Institute. The study consists of 17 clinical time intervals across 10 patients that have been followed in detail, each of whom shows significant progression of tumor over a period of 1 to 3 months on sequential follow up scans. A Taguchi sampling design is implemented to estimate the variability of the predicted tumors to using 144 different choices of model parameters. In 9 cases, model parameters can be identified such that the simulated tumor, using both models, contains at least 40 percent of the volume of the observed tumor. We discuss some potential improvements that can be made to the parameterizations of the models and their initialization.

Will vaccine-derived protective immunity curtail COVID-19 variants in the US?

Multiple effective vaccines are currently being deployed to combat the COVID-19 pandemic, and are viewed as the major factor in marked reductions of disease burden in regions with moderate to high vaccination coverage. The effectiveness of COVID-19 vaccination programs is, however, significantly threatened by the emergence of new SARS-COV-2 variants that, in addition to being more transmissible than the wild-type (original) strain, may at least partially evade existing vaccines. A two-strain (one wild-type, one variant) and two-group (vaccinated or otherwise) mechanistic mathematical model is designed and used to assess the impact of the vaccine-induced cross-protective efficacy on the spread the COVID-19 pandemic in the United States. Rigorous analysis of the model shows that, in the absence of any co-circulating SARS-CoV-2 variant, the vaccine-derived herd immunity threshold needed to eliminate the wild-type strain can be achieved if 59% of the US population is fully-vaccinated with either the Pfizer or Moderna vaccine. This threshold increases to 76% if the wild-type strain is co-circulating with the Alpha variant (a SARS-CoV-2 variant that is 56% more transmissible than the wild-type strain). If the wild-type strain is co-circulating with the Delta variant (which is estimated to be 100% more transmissible than the wild-type strain), up to 82% of the US population needs to be vaccinated with either of the aforementioned vaccines to achieve the vaccine-derived herd immunity. Global sensitivity analysis of the model reveal the following four parameters as the most influential in driving the value of the reproduction number of the variant strain (hence, COVID-19 dynamics) in the US: (a) the infectiousness of the co-circulating SARS-CoV-2 variant, (b) the proportion of individuals fully vaccinated (using Pfizer or Moderna vaccine) against the wild-type strain, (c) the cross-protective efficacy the vaccines offer against the variant strain and (d) the modification parameter accounting for the reduced infectiousness of fully-vaccinated individuals experiencing breakthrough infection. Specifically, numerical simulations of the model show that future waves or surges of the COVID-19 pandemic can be prevented in the US if the two vaccines offer moderate level of cross-protection against the variant (at least 67%). This study further suggests that a new SARS-CoV-2 variant can cause a significant disease surge in the US if (i) the vaccine coverage against the wild-type strain is low (roughly <66%) (ii) the variant is much more transmissible (e.g., 100% more transmissible), than the wild-type strain, or (iii) the level of cross-protection offered by the vaccine is relatively low (e.g., less than 50%). A new SARS-CoV-2 variant will not cause such surge in the US if it is only moderately more transmissible (e.g., the Alpha variant, which is 56% more transmissible) than the wild-type strain, at least 66% of the population of the US is fully vaccinated, and the three vaccines being deployed in the US (Pfizer, Moderna, and Johnson & Johnson) offer a moderate level of cross-protection against the variant.

To mask or not to mask: Modeling the potential for face mask use by the general public to curtail the COVID-19 pandemic.

Face mask use by the general public for limiting the spread of the COVID-19 pandemic is controversial, though increasingly recommended, and the potential of this intervention is not well understood. We develop a compartmental model for assessing the community-wide impact of mask use by the general, asymptomatic public, a portion of which may be asymptomatically infectious. Model simulations, using data relevant to COVID-19 dynamics in the US states of New York and Washington, suggest that broad adoption of even relatively ineffective face masks may meaningfully reduce community transmission of COVID-19 and decrease peak hospitalizations and deaths. Moreover, mask use decreases the effective transmission rate in nearly linear proportion to the product of mask effectiveness (as a fraction of potentially infectious contacts blocked) and coverage rate (as a fraction of the general population), while the impact on epidemiologic outcomes (death, hospitalizations) is highly nonlinear, indicating masks could synergize with other non-pharmaceutical measures. Notably, masks are found to be useful with respect to both preventing illness in healthy persons and preventing asymptomatic transmission. Hypothetical mask adoption scenarios, for Washington and New York state, suggest that immediate near universal (80%) adoption of moderately (50%) effective masks could prevent on the order of 17-45% of projected deaths over two months in New York, while decreasing the peak daily death rate by 34-58%, absent other changes in epidemic dynamics. Even very weak masks (20% effective) can still be useful if the underlying transmission rate is relatively low or decreasing: In Washington, where baseline transmission is much less intense, 80% adoption of such masks could reduce mortality by 24-65% (and peak deaths 15-69%), compared to 2-9% mortality reduction in New York (peak death reduction 9-18%). Our results suggest use of face masks by the general public is potentially of high value in curtailing community transmission and the burden of the pandemic. The community-wide benefits are likely to be greatest when face masks are used in conjunction with other non-pharmaceutical practices (such as social-distancing), and when adoption is nearly universal (nation-wide) and compliance is high.

Weather-driven malaria transmission model with gonotrophic and sporogonic cycles.

Malaria is mainly a tropical disease and its transmission cycle is heavily influenced by environment: The life-cycles of the Anopheles mosquito vector and Plasmodium parasite are both strongly affected by ambient temperature, while suitable aquatic habitat is necessary for immature mosquito development. Therefore, how global warming may affect malaria burden is an active question, and we develop a new ordinary differential equations-based malaria transmission model that explicitly considers the temperature-dependent Anopheles gonotrophic and Plasmodium sporogonic cycles. Mosquito dynamics are coupled to infection among a human population with symptomatic and asymptomatic disease carriers, as well as temporary immunity. We also explore the effect of incorporating diurnal temperature variations upon transmission. Rigorous analysis of the model show that the non-trivial disease-free equilibrium is locally-asymptotically stable when the associated reproduction number is less than unity (this equilibrium is globally-asymptotically for a special case with no density-dependent larval and disease-induced host mortality). Numerical simulations of the model, for the case where the ambient temperature is held constant, suggest a nonlinear, hyperbolic relationship between the reproduction number and clinical malaria burden. Moreover, malaria burden peaks at 29.5 o C when daily ambient temperature is held constant, but this peak decreases with increasing daily temperature variation, to about 23-25 o C. Malaria burden also varies nonlinearly with temperature, such that small temperature changes influent disease mainly at marginal temperatures, suggesting that in areas where malaria is highly endemic, any response to global warming may be highly nonlinear and most typically minimal, while in areas of more marginal malaria potential (such as the East African highlands), increasing temperatures may translate nearly linearly into increased disease potential. Finally, we observe that while explicitly modelling the stages of the Plasmodium sporogonic cycle is essential, explicitly including the stages of the Anopheles gonotrophic cycle is of minimal importance.

Mathematical modeling of climate change and malaria transmission dynamics: a historical review.

Malaria, one of the greatest historical killers of mankind, continues to claim around half a million lives annually, with almost all deaths occurring in children under the age of five living in tropical Africa. The range of this disease is limited by climate to the warmer regions of the globe, and so anthropogenic global warming (and climate change more broadly) now threatens to alter the geographic area for potential malaria transmission, as both the Plasmodium malaria parasite and Anopheles mosquito vector have highly temperature-dependent lifecycles, while the aquatic immature Anopheles habitats are also strongly dependent upon rainfall and local hydrodynamics. A wide variety of process-based (or mechanistic) mathematical models have thus been proposed for the complex, highly nonlinear weather-driven Anopheles lifecycle and malaria transmission dynamics, but have reached somewhat disparate conclusions as to optimum temperatures for transmission, and the possible effect of increasing temperatures upon (potential) malaria distribution, with some projecting a large increase in the area at risk for malaria, but others predicting primarily a shift in the disease's geographic range. More generally, both global and local environmental changes drove the initial emergence of P. falciparum as a major human pathogen in tropical Africa some 10,000 years ago, and the disease has a long and deep history through the present. It is the goal of this paper to review major aspects of malaria biology, methods for formalizing these into mathematical forms, uncertainties and controversies in proper modeling methodology, and to provide a timeline of some major modeling efforts from the classical works of Sir Ronald Ross and George Macdonald through recent climate-focused modeling studies. Finally, we attempt to place such mathematical work within a broader historical context for the "million-murdering Death" of malaria.

Principal dynamic mode analysis of the Hodgkin-Huxley equations.

We develop an autoregressive model framework based on the concept of Principal Dynamic Modes (PDMs) for the process of action potential (AP) generation in the excitable neuronal membrane described by the Hodgkin-Huxley (H-H) equations. The model's exogenous input is injected current, and whenever the membrane potential output exceeds a specified threshold, it is fed back as a second input. The PDMs are estimated from the previously developed Nonlinear Autoregressive Volterra (NARV) model, and represent an efficient functional basis for Volterra kernel expansion. The PDM-based model admits a modular representation, consisting of the forward and feedback PDM bases as linear filterbanks for the exogenous and autoregressive inputs, respectively, whose outputs are then fed to a static nonlinearity composed of polynomials operating on the PDM outputs and cross-terms of pair-products of PDM outputs. A two-step procedure for model reduction is performed: first, influential subsets of the forward and feedback PDM bases are identified and selected as the reduced PDM bases. Second, the terms of the static nonlinearity are pruned. The first step reduces model complexity from a total of 65 coefficients to 27, while the second further reduces the model coefficients to only eight. It is demonstrated that the performance cost of model reduction in terms of out-of-sample prediction accuracy is minimal. Unlike the full model, the eight coefficient pruned model can be easily visualized to reveal the essential system components, and thus the data-derived PDM model can yield insight into the underlying system structure and function.

A nonlinear autoregressive Volterra model of the Hodgkin-Huxley equations.

We propose a new variant of Volterra-type model with a nonlinear auto-regressive (NAR) component that is a suitable framework for describing the process of AP generation by the neuron membrane potential, and we apply it to input-output data generated by the Hodgkin-Huxley (H-H) equations. Volterra models use a functional series expansion to describe the input-output relation for most nonlinear dynamic systems, and are applicable to a wide range of physiologic systems. It is difficult, however, to apply the Volterra methodology to the H-H model because is characterized by distinct subthreshold and suprathreshold dynamics. When threshold is crossed, an autonomous action potential (AP) is generated, the output becomes temporarily decoupled from the input, and the standard Volterra model fails. Therefore, in our framework, whenever membrane potential exceeds some threshold, it is taken as a second input to a dual-input Volterra model. This model correctly predicts membrane voltage deflection both within the subthreshold region and during APs. Moreover, the model naturally generates a post-AP afterpotential and refractory period. It is known that the H-H model converges to a limit cycle in response to a constant current injection. This behavior is correctly predicted by the proposed model, while the standard Volterra model is incapable of generating such limit cycle behavior. The inclusion of cross-kernels, which describe the nonlinear interactions between the exogenous and autoregressive inputs, is found to be absolutely necessary. The proposed model is general, non-parametric, and data-derived.

The evolutionary impact of androgen levels on prostate cancer in a multi-scale mathematical model.

BACKGROUND: Androgens bind to the androgen receptor (AR) in prostate cells and are essential survival factors for healthy prostate epithelium. Most untreated prostate cancers retain some dependence upon the AR and respond, at least transiently, to androgen ablation therapy. However, the relationship between endogenous androgen levels and cancer etiology is unclear. High levels of androgens have traditionally been viewed as driving abnormal proliferation leading to cancer, but it has also been suggested that low levels of androgen could induce selective pressure for abnormal cells. We formulate a mathematical model of androgen regulated prostate growth to study the effects of abnormal androgen levels on selection for pre-malignant phenotypes in early prostate cancer development. RESULTS: We find that cell turnover rate increases with decreasing androgen levels, which may increase the rate of mutation and malignant evolution. We model the evolution of a heterogeneous prostate cell population using a continuous state-transition model. Using this model we study selection for AR expression under different androgen levels and find that low androgen environments, caused either by low serum testosterone or by reduced 5alpha-reductase activity, select more strongly for elevated AR expression than do normal environments. High androgen actually slightly reduces selective pressure for AR upregulation. Moreover, our results suggest that an aberrant androgen environment may delay progression to a malignant phenotype, but result in a more dangerous cancer should one arise. CONCLUSIONS: The model represents a useful initial framework for understanding the role of androgens in prostate cancer etiology, and it suggests that low androgen levels can increase selection for phenotypes resistant to hormonal therapy that may also be more aggressive. Moreover, clinical treatment with 5alpha-reductase inhibitors such as finasteride may increase the incidence of therapy resistant cancers.