Robustness of mathematical models for biological systems

The robustness of mathematical models for biological systems is studied by sensitivity analysis and stochastic simulations. Using a neural network model with three genes as the test problem, we study robustness properties of synthesis and degradation processes. For single parameter robustness, sensitivity analysis techniques are applied for studying parameter variations and stochastic simulations are used for investigating the impact of external noise. Results of sensitivity analysis are consistent with those obtained by stochastic simulations. Stochastic models with external noise can be used for studying the robustness not only to external noise but also to parameter variations. For external noise we also use stochastic models to study the robustness of the function of each gene and that of the system.

A two-phase approach for detecting recombination in nucleotide sequences

Genetic recombination can produce heterogeneous phylogenetic histories within a set of homologous genes. Delineating recombination events is important in the study of molecular evolution, as inference of such events provides a clearer picture of the phylogenetic relationships among different gene sequences or genomes. Nevertheless, detecting recombination events can be a daunting task, as the performance of different recombination-detecting approaches can vary, depending on evolutionary events that take place after recombination. We recently evaluated the effects of post-recombination events on the prediction accuracy of recombination-detecting approaches using simulated nucleotide sequence data. The main conclusion, supported by other studies, is that one should not depend on a single method when searching for recombination events. In this paper, we introduce a two-phase strategy, applying three statistical measures to detect the occurrence of recombination events, and a Bayesian phylogenetic approach in delineating breakpoints of such events in nucleotide sequences. We evaluate the performance of these approaches using simulated data, and demonstrate the applicability of this strategy to empirical data. The two-phase strategy proves to be time-efficient when applied to large datasets, and yields high-confidence results.

Sampling phylogenetic tree space with the generalized Gibbs sampler

The generalized Gibbs sampler (GGS) is a recently developed Markov chain Monte Carlo (MCMC) technique that enables Gibbs-like sampling of state spaces that lack a convenient representation in terms of a fixed coordinate system. This paper describes a new sampler, called the tree sampler, which uses the GGS to sample from a state space consisting of phylogenetic trees. The tree sampler is useful for a wide range of phylogenetic applications, including Bayesian, maximum likelihood, and maximum parsimony methods. A fast new algorithm to search for a maximum parsimony phylogeny is presented, using the tree sampler in the context of simulated annealing. The mathematics underlying the algorithm is explained and its time complexity is analyzed. The method is tested on two large data sets consisting of 123 sequences and 500 sequences, respectively. The new algorithm is shown to compare very favorably in terms of speed and accuracy to the program DNAPARS from the PHYLIP package.

A two-phase strategy for detecting recombination in nucleotide sequences

Genetic recombination can produce heterogeneous phylogenetic histories within a set of homologous genes. Delineating recombination events is important in the study of molecular evolution, as inference of such events provides a clearer picture of the phylogenetic relationships among different gene sequences or genomes. Nevertheless, detecting recombination events can be a daunting task, as the performance of different recombination-detecting approaches can vary, depending on evolutionary events that take place after recombination. We previously evaluated the effects of post-recombination events on the prediction accuracy of recombination-detecting approaches using simulated nucleotide sequence data. The main conclusion, supported by other studies, is that one should not depend on a single method when searching for recombination events. In this paper, we introduce a two-phase strategy, applying three statistical measures to detect the occurrence of recombination events, and a Bayesian phylogenetic approach to delineate breakpoints of such events in nucleotide sequences. We evaluate the performance of these approaches using simulated data, and demonstrate the applicability of this strategy to empirical data. The two-phase strategy proves to be time-efficient when applied to large datasets, and yields high-confidence results.

Detecting environmental impacts on metapopulations of mound spring invertebrates - Assessing an incidence function model

We use a stochastic patch occupancy model of invertebrates in the Mound Springs ecosystem of South Australia to assess the ability of incidence function models to detect environmental impacts on metapopulations. We assume that the probability of colonisation decreases with increasing isolation and the probability of extinction is constant across spring vents. We run the models to quasi-equilibrium, and then impose an impact by increasing the local extinction probability. We sample the output at various times pre- and postimpact, and examine the probability of detecting a significant change in population parameters. The incidence function model approach turns out to have little power to detect environmental impacts on metapopulations with small numbers of patches. (C) 2001 Elsevier Science Ltd. All rights reserved.

Patch dynamics and metapopulation theory: the case of successional species

We present a mathematical framework that combines extinction-colonization dynamics with the dynamics of patch succession. We draw an analogy between the epidemiological categorization of individuals (infected, susceptible, latent and resistant) and the patch structure of a spatially heterogeneous landscape (occupied-suitable, empty-suitable, occupied-unsuitable and empty-unsuitable). This approach allows one to consider life-history attributes that influence persistence in patchy environments (e.g., longevity, colonization ability) in concert with extrinsic processes (e.g., disturbances, succession) that lead to spatial heterogeneity in patch suitability. It also allows the incorporation of seed banks and other dormant life forms, thus broadening patch occupancy dynamics to include sink habitats. We use the model to investigate how equilibrium patch occupancy is influenced by four critical parameters: colonization rate? extinction rate, disturbance frequency and the rate of habitat succession. This analysis leads to general predictions about how the temporal scaling of patch succession and extinction-colonization dynamics influences long-term persistence. We apply the model to herbaceous, early-successional species that inhabit open patches created by periodic disturbances. We predict the minimum disturbance frequency required far viable management of such species in the Florida scrub ecosystem. (C) 2001 Academic Press.

The use and abuse of population viability analysis

A recent study by Brook ef al. empirically tested the performance of population viability analysis (PVA) using data from 21 populations across a wide range of species. The study concluded that PVAs are good at predicting the future dynamics of populations. We suggest that this conclusion is a result of a bias in the studies that Brook et al, included in their analyses, We present arguments that PVAs can only be accurate at predicting extinction probabilities if data are extensive and reliable, and if the distribution of vital rates between individuals and years can be assumed stationary in the future, or if any changes can be accurately predicted. In particular, we note th at although catastrophes are likely to have precipitated many extinctions, estimates of the probability of catastrophes are unreliable.