Mathematical models of calcium and tight junctions in normal and reconstructed epidermis

This project investigated the calcium distributions of the skin, and the growth patterns of skin substitutes grown in the laboratory, using mathematical models. The research found that the calcium distribution in the upper layer of the skin is controlled by three different mechanisms, not one as previously thought. The research also suggests that tight junctions, which are adhesions between neighbouring skin cells, cannot be solely responsible for the differences in the growth patterns of skin substitutes and normal skin.

Crinkles in connectivity: combining genetics and other types of biological data to estimate movement and interbreeding between populations

Marine species generally have large population sizes, continuous distributions and high dispersal capacity. Despite this, they are often subdivided into separate populations, which are the basic units of fisheries management. For example, populations of some fisheries species across the deep water of the Timor Trench are genetically different, inferring minimal movement and interbreeding. When connectivity is higher than the Timor Trench example, but not so high that the populations become one, connectivity between populations is crinkled. Crinkled connectivity occurs when migration is above the threshold required to link populations genetically, but below the threshold for demographic links. In future, genetic estimates of connectivity over crinkled links could be uniquely combined with other data, such as estimates of population size and tagging and tracking data, to quantify demographic connectedness between these types of populations. Elasmobranch species may be ideal targets for this research because connectivity between populations is more likely to be crinkled than for finfish species. Fisheries stock-assessment models could be strengthened with estimates of connectivity to improve the strategic and sustainable harvesting of biological resources.

Bayesian Inference for Retrospective Population Genetics Models Using Markov Chain Monte Carlo Methods

Genetics, the science of heredity and variation in living organisms, has a central role in medicine, in breeding crops and livestock, and in studying fundamental topics of biological sciences such as evolution and cell functioning. Currently the field of genetics is under a rapid development because of the recent advances in technologies by which molecular data can be obtained from living organisms. In order that most information from such data can be extracted, the analyses need to be carried out using statistical models that are tailored to take account of the particular genetic processes. In this thesis we formulate and analyze Bayesian models for genetic marker data of contemporary individuals. The major focus is on the modeling of the unobserved recent ancestry of the sampled individuals (say, for tens of generations or so), which is carried out by using explicit probabilistic reconstructions of the pedigree structures accompanied by the gene flows at the marker loci. For such a recent history, the recombination process is the major genetic force that shapes the genomes of the individuals, and it is included in the model by assuming that the recombination fractions between the adjacent markers are known. The posterior distribution of the unobserved history of the individuals is studied conditionally on the observed marker data by using a Markov chain Monte Carlo algorithm (MCMC). The example analyses consider estimation of the population structure, relatedness structure (both at the level of whole genomes as well as at each marker separately), and haplotype configurations. For situations where the pedigree structure is partially known, an algorithm to create an initial state for the MCMC algorithm is given. Furthermore, the thesis includes an extension of the model for the recent genetic history to situations where also a quantitative phenotype has been measured from the contemporary individuals. In that case the goal is to identify positions on the genome that affect the observed phenotypic values. This task is carried out within the Bayesian framework, where the number and the relative effects of the quantitative trait loci are treated as random variables whose posterior distribution is studied conditionally on the observed genetic and phenotypic data. In addition, the thesis contains an extension of a widely-used haplotyping method, the PHASE algorithm, to settings where genetic material from several individuals has been pooled together, and the allele frequencies of each pool are determined in a single genotyping.

Hazard Function Models to Estimate Mortality Rates Affecting Fish Populations with Application to the Sea Mullet (Mugil cephalus) Fishery on the Queensland Coast (Australia)

Fisheries management agencies around the world collect age data for the purpose of assessing the status of natural resources in their jurisdiction. Estimates of mortality rates represent a key information to assess the sustainability of fish stocks exploitation. Contrary to medical research or manufacturing where survival analysis is routinely applied to estimate failure rates, survival analysis has seldom been applied in fisheries stock assessment despite similar purposes between these fields of applied statistics. In this paper, we developed hazard functions to model the dynamic of an exploited fish population. These functions were used to estimate all parameters necessary for stock assessment (including natural and fishing mortality rates as well as gear selectivity) by maximum likelihood using age data from a sample of catch. This novel application of survival analysis to fisheries stock assessment was tested by Monte Carlo simulations to assert that it provided unbiased estimations of relevant quantities. The method was applied to the data from the Queensland (Australia) sea mullet (Mugil cephalus) commercial fishery collected between 2007 and 2014. It provided, for the first time, an estimate of natural mortality affecting this stock: 0.22±0.08 year −1 .

Role of Temperature in the Biological Activity of a Boreal Forest

In northern latitudes, temperature is the key factor driving the temporal scales of biological activity, namely the length of the growing season and the seasonal efficiency of photosynthesis. The formation of atmospheric concentrations of biogenic volatile organic compounds (BVOCs) are linked to the intensity of biological activity. However, interdisciplinary knowledge of the role of temperature in the biological processes related to the annual cycle and photosynthesis and atmospheric chemistry is not fully understood. The aim of this study was to improve understanding of the role of temperature in these three interlinked areas: 1) onset of growing season, 2) photosynthetic efficiency and 3) BVOC air concentrations in a boreal forest. The results present a cross-section of the role of temperature on different spatial (southern northern boreal), structural (tree forest stand - forest) and temporal (day-season- year) scales. The fundamental status of the Thermal Time model in predicting the onset of spring recovery was confirmed. However, it was recommended that sequential models would be more appropriate tools when the onset of the growing season is estimated under a warmer climate. A similar type of relationship between photosynthetic efficiency and temperature history was found in both southern and northern boreal forest stands. This result draws attention to the critical question of the seasonal efficiency of coniferous species to emit organic compounds under a warmer climate. New knowledge about the temperature dependence of the concentrations of biogenic volatile organic compounds in a boreal forest stand was obtained. The seasonal progress and the inter-correlation of BVOC concentrations in ambient air indicated a link to biological activity. Temperature was found to be the main driving factor for the concentrations. However, in addition to temperature, other factors may play a significant role here, especially when the peak concentrations are studied. There is strong evidence that the spring recovery and phenological events of many plant species have already advanced in Europe. This study does not fully support this observation. In a boreal forest, changes in the annual cycle, especially the temperature requirement in winter, would have an impact on the atmospheric BVOC composition. According to this study, more joint phenological and BVOC field observations and laboratory experiments are still needed to improve these scenarios.

Diaryl dihydropyrazole-3-carboxamides with significant in vivo antiobesity activity related to CB1 receptor antagonism: Synthesis, biological evaluation, and molecular modeling in the homology model

A number of analogues of diaryl dihydropyrazole-3-carboxamides have been synthesized. Their activities were evaluated for appetite suppression and body weight reduction in animal models. Depending on the chemical modification of the selected dihydropyrazole scaffold, the lead compoundsthe bisulfate salt of (±)-5-(4-chlorophenyl)-1-(2,4-dichlorophenyl)-4,5-dihydro-1H-pyrazole-3-carboxylic acid morpholin-4-ylamide 26 and the bisulfate salt of (−)-5-(4-chlorophenyl)-1-(2,4-dichlorophenyl)-4,5-dihydro-1H-pyrazole-3-carboxylic acid morpholin-4-ylamide 30showed significant body weight reduction in vivo, which is attributed to their CB1 antagonistic activity and exhibited a favorable pharmacokinetic profile. The molecular modeling studies also showed interactions of two isomers of (±)-5-(4-chlorophenyl)-1-(2,4-dichlorophenyl)-4,5-dihydro-1H-pyrazole-3-carboxylic acid morpholin-4-ylamide 9 with CB1 receptor in the homology model similar to those of N-piperidino-5-(4-chlorophenyl)-1-(2,4-dichlorophenyl)-4-methyl-3-pyrazole-carboxamide (rimonabant) 1 and 4S-(−)-3-(4-chlorophenyl)-N-methyl-N‘-[(4-chlorophenyl)-sulfonyl]-4-phenyl-4,5-dihydro-1H-pyrazole-1-carboxamidine (SLV-319) 2.

NUVIEW:software for display and interactive manipulation of nucleic acid models

The NUVIEW software package allows skeletal models of any double helical nucleic acid molecule to be displayed out a graphics monitor and to apply various rotations, translations and scaling transformations interactively, through the keyboard. The skeletal model is generated by connecting any pair of representative points, one from each of the bases in the basepair. In addition to the above mentioned manipulations, the base residues can be identified by using a locator and the distance between any pair of residues can be obtained. A sequence based color coded display allows easy identification of sequence repeats, such as runs of Adenines. The real time interactive manipulation of such skeletal models for large DNA/RNA double helices, can be used to trace the path of the nucleic acid chain in three dimensions and hence get a better idea of its topology, location of linear or curved regions, distances between far off regions in the sequence etc. A physical picture of these features will assist in understanding the relationship between base sequence, structure and biological function in nucleic acids.

Models of Dispersal and Diversification

Ecology and evolutionary biology is the study of life on this planet. One of the many methods applied to answering the great diversity of questions regarding the lives and characteristics of individual organisms, is the utilization of mathematical models. Such models are used in a wide variety of ways. Some help us to reason, functioning as aids to, or substitutes for, our own fallible logic, thus making argumentation and thinking clearer. Models which help our reasoning can lead to conceptual clarification; by expressing ideas in algebraic terms, the relationship between different concepts become clearer. Other mathematical models are used to better understand yet more complicated models, or to develop mathematical tools for their analysis. Though helping us to reason and being used as tools in the craftmanship of science, many models do not tell us much about the real biological phenomena we are, at least initially, interested in. The main reason for this is that any mathematical model is a simplification of the real world, reducing the complexity and variety of interactions and idiosynchracies of individual organisms. What such models can tell us, however, both is and has been very valuable throughout the history of ecology and evolution. Minimally, a model simplifying the complex world can tell us that in principle, the patterns produced in a model could also be produced in the real world. We can never know how different a simplified mathematical representation is from the real world, but the similarity models do strive for, gives us confidence that their results could apply. This thesis deals with a variety of different models, used for different purposes. One model deals with how one can measure and analyse invasions; the expanding phase of invasive species. Earlier analyses claims to have shown that such invasions can be a regulated phenomena, that higher invasion speeds at a given point in time will lead to a reduction in speed. Two simple mathematical models show that analysis on this particular measure of invasion speed need not be evidence of regulation. In the context of dispersal evolution, two models acting as proof-of-principle are presented. Parent-offspring conflict emerges when there are different evolutionary optima for adaptive behavior for parents and offspring. We show that the evolution of dispersal distances can entail such a conflict, and that under parental control of dispersal (as, for example, in higher plants) wider dispersal kernels are optimal. We also show that dispersal homeostasis can be optimal; in a setting where dispersal decisions (to leave or stay in a natal patch) are made, strategies that divide their seeds or eggs into fractions that disperse or not, as opposed to randomized for each seed, can prevail. We also present a model of the evolution of bet-hedging strategies; evolutionary adaptations that occur despite their fitness, on average, being lower than a competing strategy. Such strategies can win in the long run because they have a reduced variance in fitness coupled with a reduction in mean fitness, and fitness is of a multiplicative nature across generations, and therefore sensitive to variability. This model is used for conceptual clarification; by developing a population genetical model with uncertain fitness and expressing genotypic variance in fitness as a product between individual level variance and correlations between individuals of a genotype. We arrive at expressions that intuitively reflect two of the main categorizations of bet-hedging strategies; conservative vs diversifying and within- vs between-generation bet hedging. In addition, this model shows that these divisions in fact are false dichotomies.

Some aspects of information processing in biological vision

An attempt is made to present some challenging problems (mainly to the technically minded researchers) in the development of computational models for certain (visual) processes which are executed with, apparently, deceptive ease by the human visual system. However, in the interest of simplicity (and with a nonmathematical audience in mind), the presentation is almost completely devoid of mathematical formalism. Some of the findings in biological vision are presented in order to provoke some approaches to their computational models, The development of ideas is not complete, and the vast literature on biological and computational vision cannot be reviewed here. A related but rather specific aspect of computational vision (namely, detection of edges) has been discussed by Zucker, who brings out some of the difficulties experienced in the classical approaches.Space limitations here preclude any detailed analysis of even the elementary aspects of information processing in biological vision, However, the main purpose of the present paper is to highlight some of the fascinating problems in the frontier area of modelling mathematically the human vision system.

Influence of Spatial Scales on the Performance of Saturated Subsurface Flow and Transport Models

Numerical modeling of saturated subsurface flow and transport has been widely used in the past using different numerical schemes such as finite difference and finite element methods. Such modeling often involves discretization of the problem in spatial and temporal scales. The choice of the spatial and temporal scales for a modeling scenario is often not straightforward. For example, a basin-scale saturated flow and transport analysis demands larger spatial and temporal scales than a meso-scale study, which in turn has larger scales compared to a pore-scale study. The choice of spatial-scale is often dictated by the computational capabilities of the modeler as well as the availability of fine-scale data. In this study, we analyze the impact of different spatial scales and scaling procedures on saturated subsurface flow and transport simulations.

Analytical solutions for diffuse fluorescence spectroscopy/imaging in biological tissues. Part II: comparison and validation

The analytical solutions for the coupled diffusion equations that are encountered in diffuse fluorescence spectroscopy/ imaging for regular geometries were compared with the well-established numerical models, which are based on the finite element method. Comparison among the analytical solutions obtained using zero boundary conditions and extrapolated boundary conditions (EBCs) was also performed. The results reveal that the analytical solutions are in close agreement with the numerical solutions, and solutions obtained using EBCs are more accurate in obtaining the mean time of flight data compared to their counterpart. The analytical solutions were also shown to be capable of providing bulk optical properties through a numerical experiment using a realistic breast model. (C) 2013 Optical Society of America

A Study of Speed of the Boundary Element Method as applied to the Realtime Computational Simulation of Biological Organs

In this work, possibility of simulating biological organs in realtime using the Boundary Element Method (BEM) is investigated. Biological organs are assumed to follow linear elastostatic material behavior, and constant boundary element is the element type used. First, a Graphics Processing Unit (GPU) is used to speed up the BEM computations to achieve the realtime performance. Next, instead of the GPU, a computer cluster is used. Results indicate that BEM is fast enough to provide for realtime graphics if biological organs are assumed to follow linear elastostatic material behavior. Although the present work does not conduct any simulation using nonlinear material models, results from using the linear elastostatic material model imply that it would be difficult to obtain realtime performance if highly nonlinear material models that properly characterize biological organs are used. Although the use of BEM for the simulation of biological organs is not new, the results presented in the present study are not found elsewhere in the literature.

Correct biological timing in Arabidopsis requires multiple light-signaling pathways.

Circadian oscillators provide rhythmic temporal cues for a range of biological processes in plants and animals, enabling anticipation of the day/night cycle and enhancing fitness-associated traits. We have used engineering models to understand the control principles of a plant's response to seasonal variation. We show that the seasonal changes in the timing of circadian outputs require light regulation via feed-forward loops, combining rapid light-signaling pathways with entrained circadian oscillators. Linear time-invariant models of circadian rhythms were computed for 3,503 circadian-regulated genes and for the concentration of cytosolic-free calcium to quantify the magnitude and timing of regulation by circadian oscillators and light-signaling pathways. Bioinformatic and experimental analysis show that rapid light-induced regulation of circadian outputs is associated with seasonal rephasing of the output rhythm. We identify that external coincidence is required for rephasing of multiple output rhythms, and is therefore important in general phase control in addition to specific photoperiod-dependent processes such as flowering and hypocotyl elongation. Our findings uncover a fundamental design principle of circadian regulation, and identify the importance of rapid light-signaling pathways in temporal control.

Factors to consider when using native biological control organisms to manage exotic plants

Biological control of exotic plant populations with native organisms appears to be increasing, even though its success to date has been limited. Although many researchers and managers feel that native organisms are easier to use and present less risk to the environment this may not be true. Developing a successful management program with a native insect is dependent on a number of critical factors that need to be considered. Information is needed on the feeding preference of the agent, agent effectiveness, environmental regulation of the agent, unique requirements of the agent, population maintenance of the agent, and time to desired impact. By understanding these factors, researchers and managers can develop a detailed protocol for using the native biological control agent for a specific target plant. . We found E. lecontei in 14 waterbodies, most of which were in eastern Washington. Only one lake with weevils was located in western Washington. Weevils were associated with both Eurasian ( Myriophyllum spicatum L.) and northern watermilfoil ( M. sibiricum K.). Waterbodies with E. lecontei had significantly higher ( P < 0.05) pH (8.7 ± 0.2) (mean ± 2SE), specific conductance (0.3 ± 0.08 mS cm -1 ) and total alkalinity (132.4 ± 30.8 mg CaCO 3 L -1 ). We also found that weevil presence was related to surface water temperature and waterbody location ( = 24.3, P ≤ 0.001) and of all the models tested, this model provided the best fit (Hosmer- Lemeshow goodness-of-fit = 4.0, P = 0.9). Our results suggest that in Washington State E. lecontei occurs primarily in eastern Washington in waterbodies with pH ≥ 8.2 and specific conductance ≥ 0.2 mS cm -1 . Furthermore, weevil distribution appears to be correlated with waterbody location (eastern versus western Washington) and surface water temperature.

Stock assessment of protogynous fish: evaluating measures of spawning biomass used to estimate biological reference points

In stock assessments, recruitment is typically modeled as a function of females only. For protogynous stocks, however, disproportionate fishing on males increases the possibility of reduced fertilization rates. To incorporate the importance of males in protogynous stocks, assessment models have been used to predict recruitment not just from female spawning biomass (Sf), but also from that of males (Sm) or both sexes (Sb). We conducted a simulation study to evaluate the ability of these three measures to estimate biological reference points used in fishery management. Of the three, Sf provides best estimates if the potential for decreased fertilization is weak, whereas Sm is best only if the potential is very strong. In general, Sb estimates the true reference points most closely, which indicates that if the potential for decreased fertilization is moderate or unknown, Sb should be used in assessments of protogynous stocks. Moreover, for a broad range of scenarios, relative errors from Sf and Sb occur in opposite directions, indicating that estimates from these measures could be used to bound uncertainty.