Search PubMed⌕ Search

Biomedical subjects

N J Savill

Publications and source records attributed to N J Savill.

9 recordsLinked to original sources

Mathematical and statistical analysis of the Trypanosoma brucei slender to stumpy transition.

We propose a new model for the Stumpy Induction Factor-induced slender to stumpy transformation of Trypanosoma brucei gambiense cells in immunosuppressed mice. The model is a set of delay differential equations that describe the time-course of the infection. We fit the model, using a maximum-likelihood method, to previously published data on parasitaemia in four mice. The model is shown to be a good fit and parameter estimates and confidence intervals are derived. Our estimated parameter values are consistent with estimates from previous experimental studies. The model predicts the following. Slender cells can be classified as uncommitted, committed and dividing, and committed and non-dividing. A committed slender cell undergoes about 5 divisions before exiting the cell-cycle. Committed slender cells must produce SIF, and stumpy cells must not produce SIF. There are two mechanisms for differentiation, a background differentiation rate, and a SIF-concentration-dependent differentiation rate, which is proportional to SIF concentration. SIF has a half-life of about 1.4 h in mice. We also show, with suitable changes in the parameter values, that the model reflects behaviours seen in other host species and trypanosome strains.

Animals↗

RNA sequence evolution with secondary structure constraints: comparison of substitution rate models using maximum-likelihood methods.

We test models for the evolution of helical regions of RNA sequences, where the base pairing constraint leads to correlated compensatory substitutions occurring on either side of the pair. These models are of three types: 6-state models include only the four Watson-Crick pairs plus GU and UG; 7-state models include a single mismatch state that combines all of the 10 possible mismatches; 16-state models treat all mismatch states separately. We analyzed a set of eubacterial ribosomal RNA sequences with a well-established phylogenetic tree structure. For each model, the maximum-likelihood values of the parameters were obtained. The models were compared using the Akaike information criterion, the likelihood-ratio test, and Cox's test. With a high significance level, models that permit a nonzero rate of double substitutions performed better than those that assume zero double substitution rate. Some models assume symmetry between GC and CG, between AU and UA, and between GU and UG. Models that relaxed this symmetry assumption performed slightly better, but the tests did not all agree on the significance level. The most general time-reversible model significantly outperformed any of the simplifications. We consider the relative merits of all these models for molecular phylogenetics.

Base Pairing↗

Redundant and non-functional guide RNA genes in Trypanosoma brucei are a consequence of multiple genes per minicircle.

The mitochondrial mRNA of the parasitic protozoa Trypanosoma brucei is extensively edited by the insertion, and occasional deletion, of uridine residues. The editing is mediated by over 200 guide RNAs (gRNAs) that are encoded in circular DNA molecules called minicircles. There are some 250 different types of minicircle, called classes, with each encoding several gRNAs. Sequencing of gRNAs and minicircles has revealed a surprising amount of both redundancy, where gRNAs from different minicircle classes edit exactly the same part of an mRNA, and non-functionality, where partial or no complementarity is found between gRNA and mRNA. How does this redundancy and non-functionality arise and persist? We propose the following. Minicircle classes that contain several functional gRNA genes can be lost from the population via drift and replaced by more minicircle classes that contain fewer functional gRNA genes, on the condition that the cells keep a full complement of functional gRNAs. We demonstrate this hypothesis in a computer simulation of a model of minicircle evolution. We show that this process leads to an increasing number of minicircle classes and inevitably to only one functional gRNA per minicircle. Moreover, we show that the genome contains more minicircle classes than is actually necessary for cell survival. We also analyse the available minicircle sequence data and conclude that T. brucei is at a transient stage in this process. In addition, ten new putative gRNAs have been discovered.

Animals↗

Evolution of RNA editing in trypanosome mitochondria.

Two different RNA editing systems have been described in the kinetoplast-mitochondrion of trypanosomatid protists. The first involves the precise insertion and deletion of U residues mostly within the coding regions of maxicircle-encoded mRNAs to produce open reading frames. This editing is mediated by short overlapping complementary guide RNAs encoded in both the maxicircle and the minicircle molecules and involves a series of enzymatic cleavage-ligation steps. The second editing system is a C(34) to U(34) modification in the anticodon of the imported tRNA(Trp), thereby permitting the decoding of the UGA stop codon as tryptophan. U-insertion editing probably originated in an ancestor of the kinetoplastid lineage and appears to have evolved in some cases by the replacement of the original pan-edited cryptogene with a partially edited cDNA. The driving force for the evolutionary fixation of these retroposition events was postulated to be the stochastic loss of entire minicircle sequence classes and their encoded guide RNAs upon segregation of the single kinetoplast DNA network into daughter cells at cell division. A large plasticity in the relative abundance of minicircle sequence classes has been observed during cell culture in the laboratory. Computer simulations provide theoretical evidence for this plasticity if a random distribution and segregation model of minicircles is assumed. The possible evolutionary relationship of the C to U and U-insertion editing systems is discussed.

Animals↗

A theoretical study of random segregation of minicircles in trypanosomatids.

The kinetoplast (k) DNA network of trypanosomatids is made up of approximately 50 maxicircles and the order of 10(4) minicircles. It has been proposed, based on various observations and experiments, that the minicircles are randomly segregated between daughter cells when the parent cell divides. In this paper, this random segregation hypothesis is theoretically tested in a population dynamics model to see if it can account for the observed phenomena. The hypothesis is shown to successfully explain, in Leishmania tarentolae, the observation that there are a few major and many minor minicircle classes, the fluctuations of minicircle class copy numbers over time, the loss of non-essential minicircle classes, the long survival times of a few of these classes and that these classes are likely to be the major classes within the population. Implications of the model are examined for trypanosomatids in general, leading to several predictions. The model predicts variation in network size within a population, variation in the average network size and large-scale changes in class copy number over long time-scales, an evolutionary pressure towards larger network sizes, the selective advantage of non-random over random segregation, very strong selection for the amplified class in Crithidia fasciculata if its minicircles undergo random segregation and that Trypanosoma brucei may use sexual reproduction to maintain its viability.

Animals↗

Competition and dispersal in predator-prey waves.

Dispersing predators and prey can exhibit complex spatio-temporal wave-like patterns if the interactions between them cause oscillatory dynamics. We study the effect of these predator-prey density waves on the competition between prey populations and between predator populations with different dispersal strategies. We first describe 1- and 2-dimensional simulations of both discrete and continuous predator-prey models. The results suggest that any population that diffuses faster, disperses farther, or is more likely to disperse will exclude slower diffusing, shorter dispersing, or less likely dispersing populations, everything else being equal. It also appears that it does not matter whether time, space, or state are discrete or continuous, nor what the exact interactions between the predators and prey are. So long as waves exist the competition between populations occurs in a similar fashion. We derive a theory that qualitatively explains the observed behaviour and calculate approximate analytical solutions that describe, to a reasonable extent, these behaviours. Predictions about the cost of dispersal are tested. If strong enough, cost can reverse the populations' relative competitive strengths or lead to coexistence because of the effect of spiral wave cores. The theory is also able to explain previous results of simulations of coexistence in host-parasitoid models (Comins, H. N., and Massell, M. P., 1996, J. Theor. Biol. 183, 19-28).

Animals↗

Spatially induced speciation prevents extinction: the evolution of dispersal distance in oscillatory predator-prey models.

In a discrete-generation, individual-oriented model of predator-prey interactions that exhibits oscillations, we show that the self-structuring of the populations into spiral waves induces a selection pressure for ever-increasing dispersal distances in both populations. As the dispersal distances increase, the sizes of the spatial patterns increase, until they are too large to fit into the limited space. The patterns are then lost and the predators go extinct. This scenario, is, however, not the only outcome. A second selection pressure induced by the spatial boundary can cause reduction of the dispersal distances. Depending on the relative strengths of the two selection pressures, the predators and prey may speciate to give coexistence between short-dispersing boundary quasi-species and far-dispersing spiral quasi-species. Now, when pattern loss occurs, the predators switch to predating on the boundary prey quasi-species and do not go extinct. Also, if the populations reproduce sexually, local gene flow can inhibit the evolution of increasing dispersal distances, and hence the spatial patterns are not lost. Speciation and coexistence can also occur in the sexually reproducing species.

Biological Evolution↗

Self-reinforcing spatial patterns enslave evolution in a host-parasitoid system.

Spatially structured models of host-parasitoid interactions exhibit self-structuring into spatial patterns such as spiral waves and turbulence. We discuss the consequences of these patterns in an eco-evolutionary model of host-parasitoid interactions with evolution of the parasitoids' ability to disperse towards dense populations of hosts (termed the aggregation strength). It turns out that the direction of, and the time-scale over which the evolutionary selection pressure acts depends on the type of spatial pattern a parasitoid finds itself in. Evolution tends to reinforce the existence of the prevalent local pattern. Moreover, there is also competition between the patterns that ultimately determines the eco-evolutionary attractor. It is the interaction between multiple processes across spatial and temporal scales that leads to the rich meso-scale behaviour. Predicting the evolutionary outcome from statistical measures and subprocesses is shown to give incorrect and conflicting answers. Comparison with the behaviours of the complex Ginzburg-Landau equation shows striking similarities on which we comment.

Animals↗

Evolutionary stagnation due to pattern-pattern interactions in a coevolutionary predator-prey model.

We consider a spatially structured model of a coevolutionary predator-prey system with interactions in a one-dimensional phenotype space. We show that in phenotype space predators and prey organize themselves into distinct clusters of phenotypes called quasi-species. The prey quasi-species also cluster in patches in real space. As the prey quasi-species evolve away from the predator quasi-species (in phenotype space) the prey patch size reduces and the single predator quasi-species is inhibited from evolving toward either of the two prey species. We show that it is the interaction between the phenotype space patterns (quasi-species) and the real space patterns (patches) that inhibit the predators from evolving.

Biological Evolution↗