Changes in nasal airflow and heat transfer correlate with symptom improvement after surgery for nasal obstruction

Surgeries to correct nasal airway obstruction (NAO) often have less than desirable outcomes, partly due to the absence of an objective tool to select the most appropriate surgical approach for each patient. Computational fluid dynamics (CFD) models can be used to investigate nasal airflow, but variables need to be identified that can detect surgical changes and correlate with patient symptoms. CFD models were constructed from pre- and post-surgery computed tomography scans for 10 NAO patients showing no evidence of nasal cycling. Steady-state inspiratory airflow, nasal resistance, wall shear stress, and heat flux were computed for the main nasal cavity from nostrils to posterior nasal septum both bilaterally and unilaterally. Paired t-tests indicated that all CFD variables were significantly changed by surgery when calculated on the most obstructed side, and that airflow, nasal resistance, and heat flux were significantly changed bilaterally as well. Moderate linear correlations with patient-reported symptoms were found for airflow, heat flux, unilateral allocation of airflow, and unilateral nasal resistance as a fraction of bilateral nasal resistance when calculated on the most obstructed nasal side, suggesting that these variables may be useful for evaluating the efficacy of nasal surgery objectively. Similarity in the strengths of these correlations suggests that patient-reported symptoms may represent a constellation of effects and that these variables should be tracked concurrently during future virtual surgery planning.     

Allelic Diversity and Its Implications for the Rate of Adaptation

Genetic variation is usually estimated empirically from statistics based on population gene frequencies, but alternative statistics based on allelic diversity (number of allelic types) can provide complementary information. There is a lack of knowledge, however, on the evolutionary implications attached to allelic-diversity measures, particularly in structured populations. In this article we simulated multiple scenarios of single and structured populations in which a quantitative trait subject to stabilizing selection is adapted to different fitness optima. By forcing a global change in the optima we evaluated which diversity variables are more strongly correlated with both short- and long-term adaptation to the new optima. We found that quantitative genetic variance components for the trait and gene-frequency-diversity measures are generally more strongly correlated with short-term response to selection, whereas allelic-diversity measures are more correlated with long-term and total response to selection. Thus, allelic-diversity variables are better predictors of long-term adaptation than gene-frequency variables. This observation is also extended to unlinked neutral markers as a result of the information they convey on the demographic population history. Diffusion approximations for the allelic-diversity measures in a finite island model under the infinite-allele neutral mutation model are also provided.     

Nonlinear modelling of renal vasoaction

The control of blood pressure is a complex mixture of neural, hormonal and intrinsic interactions at the level of the heart, kidney and blood vessels. While experimental approaches to understanding these interactions are useful, it remains difficult to conduct experiments to quantify these interactions as the number of parameters increases. Thus, modelling of such physiological systems can offer considerable assistance. Typical mathematical models which describe the ability of the blood vessels to change their diameter (vasoconstriction) assume linearity of operation. However, due to the interaction of multiple vasocontrictive and vasodilative effectors, there is a significant nonlinear response to the influence of neural factors, particularly at higher levels of nerve activity (often seen in subjects with high blood pressure) which leads to low blood flow rates. This paper proposes a number of nonlinear mathematical models for the relationship between neural influences (sympathetic nerve activity (SNA)) and renal blood flow, using a feedback path to model the predominantly nonlinear effect of local vasoactive modulators such as nitric oxide, which oppose the action of SNA. The model structures are motivated by basic physiological principles, while the model parameters are determined using numerical optimisation techniques using open-loop data collected from rabbits. The models were verified by demonstrating correlation between experimental results and model outputs.     

Mathematical modeling of the population dynamics of a distinct interactions type system with local dispersal

Distinct biotic interactions in multi-species communities are a ubiquitous force in the natural ecosystem, and this force is an essential determinant of community stability and species coexistence outcomes. We conduct numerical simulations and bifurcation analysis of partial differential equations to gain better understanding and ecological insights into how predation (a), predator handling time (h), and local dispersal affect multi-species community dynamics. This system consists of resource-mutualist-exploiter-competitor interactions and local dispersal. From the inspection of our numerical simulations and co-dimension one bifurcation analysis findings, we discover several critical values that correspond to transcritical bifurcation, subcritical and supercritical Hopf bifurcations. This occurs as we vary the bifurcation parameters a and h in this complex ecological system under symmetric and asymmetric dispersal scenarios. Furthermore, the interplay between these local bifurcation points results in an exciting co-dimension two bifurcations, i.e., Bogdanov-Takens and cusp bifurcation points, respectively, which act as the synchronization points in this complex ecological system. From an ecological viewpoint, we find that (i) the effect of the no-dispersal scenario supports the maintenance of species biodiversity when the predation strength is moderate; (ii) symmetric dispersal induces both subcritical and supercritical Hopf bifurcation and support species diversity for moderate predation strength; and (iii) asymmetric dispersal promotes species diversity as it simplifies the bifurcation changes in dynamics by eliminating the subcritical bifurcations that trigger uncertainty, and this dispersal mechanism mediates species coexistence outcomes. Fundamentally, stable limit cycles have been reported as predator handling time varies in some ecological models; however, we observed in our bifurcation analysis the emergence of the unstable limit cycle as predator handling time changes. We discover that intense predator handling time destabilizes this complex ecological community. In general, our results demonstrate the influential roles of predation, predator handling time, and local dispersal in determining this system’s coexistence dynamics. This knowledge provides a better understanding of species conservation and biological control management.     

Modeling heterogeneity in nest survival data

Current statistical methods for estimating nest survival rates assume that nests are identical in their propensity to succeed. However, there are several biological reasons to question this assumption. For example, experience of the nest builder, number of nest helpers, genetic fitness of individuals, and site effects may contribute to an inherent disparity between nests with respect to their daily mortality rates. Ignoring such heterogeneity can lead to incorrect survival estimates. Our results show that constant survival models can seriously underestimate overall survival in the presence of heterogeneity. This paper presents a flexible random-effects approach to model heterogeneous nest survival data. We illustrate our methods through data on redwing blackbirds.     

Small-sample confidence regions in exponential families

This article presents an algorithm for small-sample conditional confidence regions for two or more parameters for any discrete regression model in the generalized linear interactive model family. Regions are constructed by careful inversion of conditional hypothesis tests. This method presupposes the use of approximate or exact techniques for enumerating the sample space for some components of the vector of sufficient statistics conditional on other components. Such enumeration may be performed exactly or by exact or approximate Monte Carlo, including the algorithms of Kolassa and Tanner (1994, Journal of the American Statistical Association 89, 697-702; 1999, Biometrics 55, 246-251). This method also assumes that one can compute certain conditional probabilities for a fixed value of the parameter vector. Because of a property of exponential families, one can use this set of conditional probabilities to directly compute the conditional probabilities associated with any other value of the vector of the parameters of interest. This observation dramatically reduces the computational effort required to invert the hypothesis test to obtain the confidence region. To construct a region with confidence level 1 - alpha, the algorithm begins with a grid of values for the parameters of interest. For each parameter vector on the grid (corresponding to the current null hypothesis), one transforms the initial set of conditional probabilities using exponential tilting and then calculates the p value for this current null hypothesis. The confidence region is the set of parameter values for which the p value is at least alpha.     

Numerical simulations of surface plasmon resonance system for monitoring DNA hybridization and detecting protein-lipid film interactions

This paper presents a simple method to extract information about thin organic films from surface plasmon resonance (SPR) spectra. From numerical simulations it was found that a shift (Δθ _SPR) of an absorption peak in the SPR spectrum was directly proportional to the product of the thin organic film thickness and the refractive index difference between the thin organic film and a buffer soaking the sample. It was also found that Δθ _SPR was not sensitive to the thin organic film support of a gold film and a glass cover slip. Relationships between Δθ _SPR and distributions of macromolecule structures, in the thin organic films were theoretically established. Formulae were derived for a homemade SPR system to calculate length, transverse area, density and surface concentration of macromolecules in the thin organic film. The validity of these treatments was checked by precisely measuring the size of a single distearoylphosphatidylcholine molecule on a gold-supported phospholipid film; by quantitatively monitoring hybridization of synthesized oligonucleotides strands based on a biotin/avidin system; and by quantitatively detecting the steric hindrance of rabbit C-reactive protein specifically bound to phospholipid monolayers composed of synthesized lipids. Received: 4 May 1998 / Revised version: 27 July 1998 / Accepted: 27 August 1998     

Optimal HIV treatment by maximising immune response

We present an optimal control model of drug treatment of the human immunodeficiency virus (HIV). Our model is based upon ordinary differential equations that describe the interaction between HIV and the specific immune response as measured by levels of natural killer cells. We establish stability results for the model. We approach the treatment problem by posing it as an optimal control problem in which we maximise the benefit based on levels of healthy CD4+ T cells and immune response cells, less the systemic cost of chemotherapy. We completely characterise the optimal control and compute a numerical solution of the optimality system via analytic continuation.Research supported by the Natural Science and Engineering Research Council (NSERC) and the Mathematics of Information Technology and Complex Systems (MITACS) of Canada     

Competition between near and far dispersers in spatially structured habitats

Competitive interactions and invasibility between short- and long-distance dispersal was investigated in a population on a heterogeneous landscape with spatial correlations in habitat types, and where the driving interaction between individuals is competition for space. Stochastic spatially explicit simulations were used, along with differential equation models based on pair approximations. Conditions under which either dispersal strategy can successfully invade the other were determined, as a function of the amount and clustering of suitable habitat and the relative costs involved in the two dispersal strategies. Long-distance dispersal, which reduces intraspecific competition, is sometimes advantageous even where aggregation of suitable habitat would otherwise favor short-distance dispersal, although certain habitat distributions can lead to either strategy being dominant. Coexistence is also possible on some landscapes, where the spatial structure of the populations partitions suitable sites according to the number of suitable neighboring sites. Mutual competitive exclusion, where whichever strategy is established first cannot be invaded, is also possible. All of these results are observed even when there is no intrinsic difference in the two strategies' costs, such as mortality or competitive abilities.     

Short-term changes in consumption and oviposition rates of Neoseiulus californicus strains (Acari: Phytoseiidae) after a diet shift

Short-term effects on consumption and oviposition rates of four strainsof Neoseiulus californicus (McGregor) after a diet shiftwere evaluated. The new feeding conditions experienced by the predators weresixfixed densities of eggs or protonymphs of Tetranychusurticae Koch placed on excised strawberry leaflet discs andmaintained under laboratory conditions (25±1°C,75–85% RH, 16L: 8D). The observations were made on the first and thefifthday of the experiment. The phytoseiids came from three long-term mass-rearedstrains fed on T. urticae, Dermatophagoides farinaeHughes,or Quercus spp. pollen, respectively. The fourth strainwascollected directly in a strawberry field. Time since diet transfer can be addedto the factors (i.e. feeding history and prey density) already known to affectthe functional and numerical responses of N. californicus,both when it feeds on prey eggs and protonymphs. If consumption rates wereaveraged over all strains and densities, 9.04 and 11.41 eggs, and 6.97 and 6.48protonymphs were consumed on the first and the fifth day, respectively. If thesame was done for oviposition rates, predators feeding on eggs produced 1.46and2.36 eggs/female/day, whereas predators feeding on protonymphs produced 1.35and2.29 eggs/female/day. Time had the greatest impact on the functional responseofthe strain that had previously fed on tetranychids, while an effect of time onthe numerical response was detectable in all strains.