Ecological modeling and pest population management: a possible and necessary connection in a changing world

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)

Population dynamics, life stage and ecological modeling in Chrysomya albiceps (Wiedemann) (Diptera: Calliphoridae)

Coordenação de Aperfeiçoamento de Pessoal de Nível Superior (CAPES)

Mathematical modeling and control of population systems: Applications in biological pest control

The aim of this paper is to apply methods from optimal control theory, and from the theory of dynamic systems to the mathematical modeling of biological pest control. The linear feedback control problem for nonlinear systems has been formulated in order to obtain the optimal pest control strategy only through the introduction of natural enemies. Asymptotic stability of the closed-loop nonlinear Kolmogorov system is guaranteed by means of a Lyapunov function which can clearly be seen to be the solution of the Hamilton-Jacobi-Bellman equation, thus guaranteeing both stability and optimality. Numerical simulations for three possible scenarios of biological pest control based on the Lotka-Volterra models are provided to show the effectiveness of this method. (c) 2007 Elsevier B.V. All rights reserved.

Population genetic and dispersal modeling data for Bathymodiolus mussels from the Mid-Atlantic Ridge

The zip folder comprises a text file and a gzipped tar archive. 1) The text file contains individual genotype data for 90 SNPs, 9 microsatellites and the mitochondrial ND4 gene that were determined in deep-sea hydrothermal vent mussels from the Mid-Atlantic Ridge (genus Bathymodiolus). Mussel specimens are grouped according to the population (pop)/location from which they have been sampled (first column). The remaining columns contain the respective allele/haplotype codes for the different genetic loci (names in the header line). The data file is in CONVERT format and can be directly transformed into different input files for population genetic statistics. 2) The tar archive contains NetCDF files with larval dispersal probabilities for simulated annual larval releases between 1998 and 2007. For each simulated vent location (Menez Gwen, Lucky Strike, Rainbow, Vent 1-10) two NetCDF files are given, one for an assumed pelagic larval duration of 1 year and the other one for an assumed pelagic larval duration of 6 months (6m).

Bayesian latent trait modeling of migraine symptom data

Definition of disease phenotype is a necessary preliminary to research into genetic causes of a complex disease. Clinical diagnosis of migraine is currently based on diagnostic criteria developed by the International Headache Society. Previously, we examined the natural clustering of these diagnostic symptoms using latent class analysis (LCA) and found that a four-class model was preferred. However, the classes can be ordered such that all symptoms progressively intensify, suggesting that a single continuous variable representing disease severity may provide a better model. Here, we compare two models: item response theory and LCA, each constructed within a Bayesian context. A deviance information criterion is used to assess model fit. We phenotyped our population sample using these models, estimated heritability and conducted genome-wide linkage analysis using Merlin-qtl. LCA with four classes was again preferred. After transformation, phenotypic trait values derived from both models are highly correlated (correlation = 0.99) and consequently results from subsequent genetic analyses were similar. Heritability was estimated at 0.37, while multipoint linkage analysis produced genome-wide significant linkage to chromosome 7q31-q33 and suggestive linkage to chromosomes 1 and 2. We argue that such continuous measures are a powerful tool for identifying genes contributing to migraine susceptibility.

Estimation of parameters for macroparasite population evolution using approximate Bayesian computation

We estimate the parameters of a stochastic process model for a macroparasite population within a host using approximate Bayesian computation (ABC). The immunity of the host is an unobserved model variable and only mature macroparasites at sacrifice of the host are counted. With very limited data, process rates are inferred reasonably precisely. Modeling involves a three variable Markov process for which the observed data likelihood is computationally intractable. ABC methods are particularly useful when the likelihood is analytically or computationally intractable. The ABC algorithm we present is based on sequential Monte Carlo, is adaptive in nature, and overcomes some drawbacks of previous approaches to ABC. The algorithm is validated on a test example involving simulated data from an autologistic model before being used to infer parameters of the Markov process model for experimental data. The fitted model explains the observed extra-binomial variation in terms of a zero-one immunity variable, which has a short-lived presence in the host.

Birth outcomes after spontaneous or assisted conception among infertile Australian women aged 28 to 36 years: A prospective, population-based study

Objective To examine the extent to which the odds of birth, pregnancy, or adverse birth outcomes are higher among women aged 28 to 36 years who use fertility treatment compared with untreated women. Design Prospective, population-based. Setting Not applicable. Patient(s) Participants in the ALSWH born in 1973 to 1978 who reported on their infertility and use of in vitro fertilization (IVF) or ovulation induction (OI). Intervention(s) Postal survey questionnaires administered as part of ALSWH. Main Outcome Measure(s) Among women treated with IVF or OI and untreated women, the odds of birth outcomes estimated by use of adjusted logistic regression modeling. Result(s) Among 7,280 women, 18.6% (n = 1,376) reported infertility. Half (53.0%) of the treated women gave birth compared with 43.8% of untreated women. Women with prior parity were less likely to use IVF compared with nulliparous women. Women using IVF or OI, respectively, were more likely to have given birth after treatment or be pregnant compared with untreated women. Women using IVF or OI were as likely to have ectopic pregnancies, stillbirths, or premature or low birthweight babies as untreated women. Conclusion(s) More than 40% of women aged 28–36 years reporting a history of infertility can achieve births without using treatment, indicating they are subfertile rather than infertile.

Using a pharmacokinetic model to interpret biomonitoring data of PBDEs in the Australian population

From human biomonitoring data that are increasingly collected in the United States, Australia, and in other countries from large-scale field studies, we obtain snap-shots of concentration levels of various persistent organic pollutants (POPs) within a cross section of the population at different times. Not only can we observe the trends within this population with time, but we can also gain information going beyond the obvious time trends. By combining the biomonitoring data with pharmacokinetic modeling, we can re-construct the time-variant exposure to individual POPs, determine their intrinsic elimination half-lives in the human body, and predict future levels of POPs in the population. Different approaches have been employed to extract information from human biomonitoring data. Pharmacokinetic (PK) models were combined with longitudinal data1, with single2 or multiple3 average concentrations of a cross-sectional data (CSD), or finally with multiple CSD with or without empirical exposure data4. In the latter study, for the first time, the authors based their modeling outputs on two sets of CSD and empirical exposure data, which made it possible that their model outputs were further constrained due to the extensive body of empirical measurements. Here we use a PK model to analyze recent levels of PBDE concentrations measured in the Australian population. In this study, we are able to base our model results on four sets5-7 of CSD; we focus on two PBDE congeners that have been shown3,5,8-9 to differ in intake rates and half-lives with BDE-47 being associated with high intake rates and a short half-life and BDE-153 with lower intake rates and a longer half-life. By fitting the model to PBDE levels measured in different age groups in different years, we determine the level of intake of BDE-47 and BDE-153, as well as the half-lives of these two chemicals in the Australian population.

Understanding mild acid pretreatment of sugarcane bagasse through particle scale modeling

Sugarcane bagasse is an abundant and sustainable resource, generated as a by-product of sugarcane milling. The cellulosic material within bagasse can be broken down into glucose molecules and fermented to produce ethanol, making it a promising feedstock for biofuel production. Mild acid pretreatment hydrolyses the hemicellulosic component of biomass, thus allowing enzymes greater access to the cellulosic substrate during saccharification. A particle-scale mathematical model describing the mild acid pretreatment of sugarcane bagasse has been developed, using a volume averaged framework. Discrete population-balance equations are used to characterise the polymer degradation kinetics, and diffusive effects account for mass transport within the cell wall of the bagasse. As the fibrous material hydrolyses over time, variations in the porosity of the cell wall and the downstream effects on the reaction kinetics are accounted for using conservation of volume arguments. Non-dimensionalization of the model equations reduces the number of parameters in the system to a set of four dimensionless ratios that compare the timescales of different reaction and diffusion events. Theoretical yield curves are compared to macroscopic experimental observations from the literature and inferences are made as to constraints on these “unknown” parameters. These results enable connections to be made between experimental data and the underlying thermodynamics of acid pretreatment. Consequently, the results suggest that data-fitting techniques used to obtain kinetic parameters should be carefully applied, with prudent consideration given to the chemical and physiological processes being modeled.

Modeling predictors of risky drug use behavior among male street laborers in urban Vietnam

Background The application of theoretical frameworks for modeling predictors of drug risk among male street laborers remains limited. The objective of this study was to test a modified version of the IMB (Information-Motivation-Behavioral Skills Model), which includes psychosocial stress, and compare this modified version with the original IMB model in terms of goodness-of-fit to predict risky drug use behavior among this population. Methods In a cross-sectional study, social mapping technique was conducted to recruit 450 male street laborers from 135 street venues across 13 districts of Hanoi city, Vietnam, for face-to-face interviews. Structural equation modeling (SEM) was used to analyze data from interviews. Results Overall measures of fit via SEM indicated that the original IMB model provided a better fit to the data than the modified version. Although the former model was able to predict a lesser variance than the latter (55% vs. 62%), it was of better fit. The findings suggest that men who are better informed and motivated for HIV prevention are more likely to report higher behavioral skills, which, in turn, are less likely to be engaged in risky drug use behavior. Conclusions This was the first application of the modified IMB model for drug use in men who were unskilled, unregistered laborers in urban settings. An AIDS prevention program for these men should not only distribute information and enhance motivations for HIV prevention, but consider interventions that could improve self-efficacy for preventing HIV infection. Future public health research and action may also consider broader factors such as structural social capital and social policy to alter the conditions that drive risky drug use among these men.

Applying quantile regression for modeling equivalent property damage only crashes to identify accident blackspots

Hot spot identification (HSID) aims to identify potential sites—roadway segments, intersections, crosswalks, interchanges, ramps, etc.—with disproportionately high crash risk relative to similar sites. An inefficient HSID methodology might result in either identifying a safe site as high risk (false positive) or a high risk site as safe (false negative), and consequently lead to the misuse the available public funds, to poor investment decisions, and to inefficient risk management practice. Current HSID methods suffer from issues like underreporting of minor injury and property damage only (PDO) crashes, challenges of accounting for crash severity into the methodology, and selection of a proper safety performance function to model crash data that is often heavily skewed by a preponderance of zeros. Addressing these challenges, this paper proposes a combination of a PDO equivalency calculation and quantile regression technique to identify hot spots in a transportation network. In particular, issues related to underreporting and crash severity are tackled by incorporating equivalent PDO crashes, whilst the concerns related to the non-count nature of equivalent PDO crashes and the skewness of crash data are addressed by the non-parametric quantile regression technique. The proposed method identifies covariate effects on various quantiles of a population, rather than the population mean like most methods in practice, which more closely corresponds with how black spots are identified in practice. The proposed methodology is illustrated using rural road segment data from Korea and compared against the traditional EB method with negative binomial regression. Application of a quantile regression model on equivalent PDO crashes enables identification of a set of high-risk sites that reflect the true safety costs to the society, simultaneously reduces the influence of under-reported PDO and minor injury crashes, and overcomes the limitation of traditional NB model in dealing with preponderance of zeros problem or right skewed dataset.

Use of simple pharmacokinetic modeling to characterize exposure of Australians to perfluorooctanoic acid and perfluorooctane sulfonic acid

Perflurooctanoic acid (PFOA) and perfluorooctane sulfonic acid (PFOS) have been used for a variety of applications including fluoropolymer processing, fire-fighting foams and surface treatments since the 1950s. Both PFOS and PFOA are polyfluoroalkyl chemicals (PFCs), man-made compounds that are persistent in the environment and humans; some PFCs have shown adverse effects in laboratory animals. Here we describe the application of a simple one compartment pharmacokinetic model to estimate total intakes of PFOA and PFOS for the general population of urban areas on the east coast of Australia. Key parameters for this model include the elimination rate constants and the volume of distribution within the body. A volume of distribution was calibrated for PFOA to a value of 170ml/kgbw using data from two communities in the United States where the residents' serum concentrations could be assumed to result primarily from a known and characterized source, drinking water contaminated with PFOA by a single fluoropolymer manufacturing facility. For PFOS, a value of 230ml/kgbw was used, based on adjustment of the PFOA value. Applying measured Australian serum data to the model gave mean+/-standard deviation intake estimates of PFOA of 1.6+/-0.3ng/kgbw/day for males and females >12years of age combined based on samples collected in 2002-2003 and 1.3+/-0.2ng/kg bw/day based on samples collected in 2006-2007. Mean intakes of PFOS were 2.7+/-0.5ng/kgbw/day for males and females >12years of age combined based on samples collected in 2002-2003, and 2.4+/-0.5ng/kgbw/day for the 2006-2007 samples. ANOVA analysis was run for PFOA intake and demonstrated significant differences by age group (p=0.03), sex (p=0.001) and date of collection (p<0.001). Estimated intake rates were highest in those aged >60years, higher in males compared to females, and higher in 2002-2003 compared to 2006-2007. The same results were seen for PFOS intake with significant differences by age group (p<0.001), sex (p=0.001) and date of collection (p=0.016).

Maternal psychological distress mediates the relationship between asthma and physician visits in a population-based sample of adolescents

Objective This study examined whether maternal psychological distress mediates the relationship between presence of adolescent asthma and number of physician visits, and whether the association between maternal psychological distress and physician visits is moderated by adolescent general health. Methods Data were obtained from the Mater University Study of Pregnancy and included 4025 adolescents. Path analysis was used to examine mediating and moderating effects. Results Maternal psychological distress was found to partially mediate the relationship between adolescent asthma and number of physician visits, accounting for 25% of the effect of adolescent asthma on physician visits (p = .046). There was no evidence to suggest that adolescent general health moderated the association between maternal psychological distress and physician visits (p = .093). Conclusions The findings suggest that maternal psychological distress is associated with increased physician visits, regardless of adolescents' general health. Lowering maternal psychological distress may serve to reduce health care utilization and costs among adolescents with asthma.