Email updates

Keep up to date with the latest news and content from BMC Evolutionary Biology and BioMed Central.

Open Access Research article

Evolution of dispersal and life history strategies – Tetrahymena ciliates

Else J Fjerdingstad12*, Nicolas Schtickzelle13, Pauline Manhes14, Arnaud Gutierrez15 and Jean Clobert16

Author Affiliations

1 Laboratoire d'Ecologie, CNRS UMR 7625, Université Pierre et Marie Curie, Paris, France

2 Now at the Department of Biology, Queens College, City University of New York, Flushing, NY, USA

3 Biodiversity Research Centre, Université catholique de Louvain, Croix du Sud 4, 1348 Louvain-la-Neuve, Belgium

4 Now at the Department of Biology, Indiana University, Bloomington, IN, USA

5 Laboratoire d'Ecologie, CNRS UMR 7625, Université Pierre et Marie Curie, Paris, France

6 Station d'Ecologie Expérimentale du CNRS à Moulis, Laboratoire Evolution et Diversité Biologique, Moulis, 09200 Saint-Girons, France

For all author emails, please log on.

BMC Evolutionary Biology 2007, 7:133  doi:10.1186/1471-2148-7-133

The electronic version of this article is the complete one and can be found online at: http://www.biomedcentral.com/1471-2148/7/133


Received:15 February 2007
Accepted:6 August 2007
Published:6 August 2007

© 2007 Fjerdingstad et al; licensee BioMed Central Ltd.

This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract

Background

Considerable attention has focused on how selection on dispersal and other core life-history strategies (reproductive effort, survival ability, colonization capacity) may lead to so-called dispersal syndromes. Studies on genetic variation in these syndromes within species could importantly increase our understanding of their evolution, by revealing whether traits co-vary across genetic lineages in the manner predicted by theoretical models, and by stimulating further hypotheses for experimental testing. Yet such studies remain scarce. Here we studied the ciliated protist Tetrahymena thermophila, a particularly interesting organism due to cells being able to transform into morphs differing dramatically in swim-speed. We investigated dispersal, morphological responses, reproductive performance, and survival in ten different clonal strains. Then, we examined whether life history traits co-varied in the manner classically predicted for ruderal species, examined the investment of different strains into short- and putative long-distance dispersal, while considering also the likely impact of semi-sociality (cell aggregation, secretion of 'growth factors') on dispersal strategies.

Results

Very significant among-strain differences were found with regard to dispersal rate, morphological commitment and plasticity, and almost all core life-history traits (e.g. survival, growth performance and strategy), with most of these traits being significantly intercorrelated. Some strains showed high short-distance dispersal rates, high colonization capacity, bigger cell size, elevated growth performance, and good survival abilities. These well performing strains, however, produced fewer fast-swimming dispersal morphs when subjected to environmental degradation than did philopatric strains performing poorly under normal conditions.

Conclusion

Strong evidence was found for a genetic covariation between dispersal strategies and core life history traits in T. thermophila, with a fair fit of observed trait associations with classic colonizer models. However, the well performing strains with high colonization success and short-distance dispersal likely suffered under a long-distance dispersal disadvantage, due to producing fewer fast-swimming dispersal morphs than did philopatric strains. The smaller cell size at carrying capacity of the latter strains and their poor capacity to colonize as individual cells suggest that they may be adapted to greater levels of dependency on clone-mates (stronger sociality). In summary, differential exposure to selection on competitive and cooperative abilities, in conjunction with selective factors targeting specifically dispersal distance, likely contributed importantly to shaping T. thermophila dispersal and life history evolution.

Background

Understanding the selective pressures affecting the evolution of dispersal strategies is of prime importance for a broad range of biological fields, ranging from conservation biology to research on the evolution of species, host-parasite interactions and communities of species [1-7]. Dispersal strategies of living organisms affect the dynamics, and the demographic and genetic structure of their populations [8,9], and impacts often crucially on the survival and reproductive success of individuals [2,10,11]. For example, dispersal has been theoretically demonstrated to profoundly influence the evolution of sociality [12-17] and dispersal strategies affect selection on other core life-history traits, e.g. reproductive strategies and effort, survival capacity, colonization ability, defenses against predators, parasites and diseases [7,18-23].

Co-variation among life-history traits has been the centre of much research in the past decades in the context of adaptive evolution as well as from the point of view of constraints on this evolution due to trade-offs between traits [24-29], and the existence of syndromes encapsulating dispersal strategies has been hypothesized [21,23,30,31]. For example, Baker and Stebbins [30] hypothesized that species living in unstable habitats in metapopulations with a high turn-over develop a set of co-adapted traits where high dispersal and colonization ability are linked to high fecundity and short life span (low survival), allowing excellent exploitation of freed-up patches (successful dispersal and colonization followed by rapid growth of the colonist population) [32]. Such hypotheses have been supported by correlative approaches, especially by comparing species in different types of habitats [23]. Theoretical models equally predict the evolution of syndromes in which dispersal ability influences or is influenced by the evolution of other life history traits [18,33], even if sometimes they come to a different conclusion than verbal or correlative approaches [34]. For example, Ronce and Olivieri [19] predicted a positive association between life span and dispersal while Crowley and McLetchie [35] predicted the reverse, the direction of the trend seemingly depending on details of the functioning of metapopulations [34] and also, we note, on potential constraints (trade-offs) on life-history traits and their capacity to evolve [25-28]. Similar complex theoretical results have also been found for the evolution of social traits and dispersal ability [17].

The evolution of dispersal syndromes through natural selection requires that the traits in question are under genetic determinism, and we must therefore demonstrate this to corroborate the theoretical view on dispersal syndrome evolution. Detailed studies on genetic variation in life history trait associations may also increase our understanding of the process of dispersal syndrome evolution, by allowing an evaluation of the degree to which the directionality of these patterns follow theoretical predictions or might be limited by trade-offs, and can also suggest additional selective factors that may have shaped syndromes. Nevertheless, only quite few studies test for and explore variation in dispersal syndromes apart from species comparisons, although genetic variance for dispersal itself has been demonstrated for some species (e.g. insects, birds, reptiles, mammals, [36-47]). Some correlative studies have found differences in the life-history of dispersing versus philopatric individuals [10,11] or intra-specific differences in dispersal strategies related with the cause of dispersal and landscape structure [48,49]. Other studies have found clear associations between individual propensity to disperse and some morphological, physiological or behavioral characteristics [50-54]. Such associations may, however, be due to maternal or environmental factors (i.e. forms of phenotypic plasticity, review by [27,47,55,56] as well as to genetic variation. Overall, if we except species with a clear dispersal-dedicated apparatus [24,45,47,58], evidence for genetic variation in dispersal syndromes and close dissections of trait associations are scarce (for a few exceptions [39,41,44], see also [38]). This means that our empirical and experimental insight into the evolution of these syndromes remains relatively poor.

One reason for the scarcity of evidence for genetic variation in dispersal strategies and other life history traits within species is that most studies have been done either on vertebrates, where assessing the genetic determinism of traits associations is difficult, or on invertebrates where an individual following up is often impossible. For this reason, artificial microcosms and clonal organisms with short generation time seem one of the best ways to investigate such syndromes ([23,59], see also [60]). Microorganisms, characterized by the general ease of maintaining large population sizes in the laboratory under controlled environmental conditions, should therefore be very well-suited organisms for studies on the co-variation between dispersal strategies and other core life history strategies. Surprisingly, however, very few studies have addressed the relationships between dispersal rate and core life history traits in microorganisms ([23], but see [60]).

We here present a study on the co-variation of dispersal strategies with other life-history traits in the unicellular, ciliated protozoan Tetrahymena thermophila. This small (60 μm) eukaryote feeds on bacteria and dissolved nutrients in fresh water ponds and streams in America [61,62], and while it is widely used as a model system by molecular and cell biologists [62], it has been surprisingly little studied by evolutionary biologists (but see [63-65]).

Yet the life history characteristics of T. thermophila make it particularly exciting for studies on dispersal. Firstly, this organism lives in habitats likely resembling those of meta-populations with high turn-over of local patches. Studies on life history trait associations in T. thermophila would therefore allow gaining insight into the degree to which central theoretical models can predict dispersal syndrome evolution in ruderal species (e.g. the classic colonizer syndrome model [23,30,31]). Secondly, genetic variation in all core life history traits can be easily estimated in T. thermophila because separate clonal lineages can be kept in the laboratory. Reproduction remains clonal whenever nutrients are present [66] and even under conditions inducing sexual reproduction (such as starvation), conjugation is impossible between clone mates because they carry the same mating type ([67], reviewed by [68]). Thirdly, T. thermophila shows an intriguing co-existence of apparent short- versus long-distance dispersal strategies (by normal cells versus cells that have transformed into elongated morphs with very numerous ciliae and a caudal flagellum [69]). T. thermophila of both morph-types swim about and explore their environment (personal observations), but the much greater swim speed (4–5 times [69,70]) and more directional movements (personal observations, see also [71]) reported for the elongated morph suggest that it is specialized for committed long-distance dispersal. Finally, T. thermophila cells form aggregations and secrete substances favoring the survival of other cells [72-74]; the evolution of dispersal strategies in this organism may therefore also be affected by a balance between kin-benefits and -competition (see also [16]).

Through experiments in the laboratory, we assessed variation in dispersal rate (in a two-patch system) and differences in the colonization capacity of single cells, and studied growth rate, patch carrying capacity, starvation resistance and concomitant changes in cell shape for ten strains of T. thermophila. These experiments and observations allowed us to (1) test for genetic variation and co-variation in core life history traits, (2) examine whether the directionality of life history trait associations across genetic lineages followed predictions for ruderal organisms (in particular as concerns the classic colonizer syndrome and competition-colonization co-existence) and (3) generate testable hypotheses on how the semi-social life style of T. thermophila may affect the evolution of its dispersal strategies.

Results

Growth from low density in presence of nutrients

All strains displayed a similar pattern of logistic growth, first exponential growth then reaching a plateau, where cells became smaller and rounder in a first time due to the rapid cell divisions and returned to their starting size and shape values once again when reaching the plateau. However, strains differed quantitatively. Growth rate and carrying capacity were largely strain-dependent (ΔAIC model with no effect versus model with a strain effect = 222.73), while variation between replicates within strain showed significant variation for some strains only and to a much lesser extent (ΔAIC model with a strain effect versus model with strain and replicate effects = 4.74). At the end of the experiment, where strains were at their carrying capacity, strains also differed in cell density (F9,18 = 10.72, P < 0.0001), size (F9,18= 17.47, P < 0.0001), and shape (F9,18 = 18.95, P < 0.0001).

Strains differed with respect to the co-variation among traits as well. This was shown by a principal component analysis on growth rate, final cell density, final cell size and final cell shape (carrying capacity K was dropped because it correlated tightly with final cell density), followed up by a discriminant analysis on the principal components (PCs). Such discriminant analyses test whether strains can be significantly distinguished from each other on the basis of the variables under study; they do not force artefactual strain separation [94]. All PCA components were significantly implicated in discriminating strains (PC1G: F9,20 = 48.55, P < 0.0001, R2 = 0.96; PC2G: F9,20 = 2.45, P = 0.045, R2 = 0.52; PC3G: F9,20 = 8.40, P < 0.0001, R2 = 0.79) with 80 % of the replicates well classified within their strain. PC1G (explaining 46 % of the variance) represented a factor that was negatively associated with final cell size and positively with final cell shape (Figure 1). This means that at the carrying capacity some strains had small and elongated cells, others rounder and bigger cells. PC2G (explaining 27 % of the variance) represented a contrast between growth rate and final cell density (Figure 1). Thus, some strains grew rapidly in the beginning but reached only a low final cell density (r strategy), whereas others grew slowly in the beginning but reached a high final density (K strategy). PC3G (explaining 20 % of the variance) represented the overall performance of growth (r and K) from low density in the presence of nutrients, contrasting strains performing well (e.g. 7) to strains performing less well (e.g. D3) (Figure 1).

thumbnailFigure 1. Principal component plots expressing correlations between cell morphology and growth variables in presence of nutrients for the ten T. thermophila strains studied. See text for details. PC1G expressed the opposition between strains with big round versus small elongated cells at carrying capacity. PC2G illustrated the r versus K strategy of growth, while PC3G represented overall growth performance.

Starvation in a medium devoid of nutrients

At the start of the experiment, cells were still dividing as if they were in a nutrient-rich environment. Although there was some variation between replicates and strains, a peak was observed at 8 h for density and cell elongation. Then density was steadily decreasing and cell shape rounding. Strains, however, differed significantly in survival rate, as estimated by survival as a density sum over time (F9,20 = 55.70, P < 0.0001). They also differed significantly as regarded cell elongation in response to starvation, as measured by mean (F9,20 = 90.96, P < 0.0001) and variance (F9,20 = 6.11, P = 0.0004) of maximal elongation, elongation persistence (F9,20 = 742.96, P < 0.0001) and frequency of disperser morphs (F9,20 = 15.80, P < 0.0001).

Strains differed regarding the impact of starvation on trait associations, too, as was shown by a principal component analysis on cell survival, elongation, and production of disperser morphs, followed up by a discriminant analysis on the principal components (PCs). The first two PCA components explained 81 % of the variance in the data, were highly significantly implicated in discriminating strains (PC1S: F9,20 = 106.73, P < 0.0001; PC2S: F9,20 = 19.67, P < 0.0001) with all replicates well classified within their strain (P = 1). The first principal component, PC1S, (explaining 57 % of the variance) expressed the overall survival performance, cell elongation and production of disperser morphs in starvation conditions, positively associated with all five included variables (Figure 2). PC2S (explaining 24 % of the variance) represented a contrast of elongation reactions: strains with a larger and persisting mean maximal elongation but a limited variance in elongation (everyone elongates relatively similarly) versus strains with a larger variation between cells, some elongating far more than other, up to becoming real disperser morphs (Figure 2).

thumbnailFigure 2. Principal component plot representing correlations between cell morphology and survival variables under starvation for the ten T. thermophila strains. PC1S is linked to overall survival and elongation capabilities. PC2S represented the cell elongation strategy, opposing strains where all cells elongates similarly for a long time to strains where some cells elongate more than others, up to becoming dispersal morphs.

Dispersal in presence of nutrients

The rates at which T. thermophila cells dispersed from one tube (the start patch) through a connecting tubing to another tube (the target patch) were strongly positively correlated with the degree of cell elongation (Figure 3) and with the initial shape of cells(r = 0.34, n = 55, P = 0.011). A linear model analysis revealed that the differences in dispersal rate were primarily situated between strains (F9,43 = 10.64, P < 0.0001). When controlling for this strain effect, the correlation of dispersal rate with cell elongation was still significant (F1,43 = 21.46, P < 0.0001), but not the correlation with initial shape (F1,43 = 0.15, P = 0.701). Elongation also significantly differed between strains (F9,43 = 3.01, P = 0.007). Our tests hence showed that dispersal rates vary significantly among strains and also supported that elongation of T. thermophila cells is linked with dispersal.

thumbnailFigure 3. Correlation of dispersal rate and cell elongation in the dispersal experiment in presence of nutrients for the ten T. thermophila strains studied.

Single cell colonization capacity in presence of nutrients

Strains differed significantly as concerned the probability of successful colonization in new patches while replicates did not differ significantly (linear model with logit link and binomial error distribution: strain: X29 = 71.42, P < 0.0001; replicate: X22 = 0.000, P = 1.000; replicate*strain: X28 = 5.93, P = 0.655; CV of strain means = 0.595).

Correlation among life-history traits across experiments

Across experiments, strains that showed a great dispersal rate in the two-patch dispersal experiment in presence of nutrients also had a greater probability of successfully colonizing new patches as single cells (Figure 4A). Dispersal and colonization capacities were not associated with the cell sizes and shapes attained at carrying capacity after growth from low density in the presence of nutrients, as shown by their non-significant co-variation with PC1G (dispersal rate: mean Spearman correlation coefficient r = -0.04, n = 30, proportion of significant correlations s = 0%, P = 1; colonization probability: r = -0.06, n = 19, s = 0%, P = 1), nor with the trade-off between growth rate and cell density at carrying capacity (PC2G) (dispersal rate: r = -0.09, n = 30, s = 0%, P = 1; colonization probability: r = 0.20, n = 19, s = 3%, P = 0.998). However, dispersal rate and colonization probability were positively associated with the overall growth performance (PC3G) (dispersal rate: Figure 4B; colonization probability: r = 0.46, n = 19, s = 50%, P < 0.001). Dispersal rate was positively associated with the overall survival capacity and trend to elongate when faced with starvation in a nutrient-free medium (PC1S: Figure 4C), but colonization probability was not (r = 0.03, n = 19, s = 0%, P = 1). Dispersal rate and colonization probability in a nutrient-rich medium were higher for strains where all cells elongated similarly and for a long time but did not produce dispersal morphs, as shown by a negative association with PC2S (dispersal rate: Figure 4D; colonization probability: r = -0.51, n = 19, s = 76% P < 0.001). This elongation strategy (PC2S) was also associated with strains presenting big round (versus small elongated) cells at carrying capacity (PC1G: Figure 4E) and to the overall performance of growth in nutrient-rich medium (PC3G: Figure 4F). Furthermore, strains showing high carrying capacities but low growth rates (K strategy strains) in presence of nutrients also showed superior abilities to survive and elongate under starvation conditions (Figure 4G). Finally, elongation in nutrient-rich and nutrient-free (starvation) conditions were highly correlated (Figure 4H) but elongation was always greatly more pronounced under starvation conditions.

thumbnailFigure 4. Correlation of life-history traits for the ten T. thermophila strains we studied. A) In the experiment on dispersal in presence of nutrients, a higher dispersal rate was associated with larger cell elongation. B) Strains with a high dispersal rate also had a greater probability of successfully colonizing a new patch as single cells. Strains adept at dispersing (C) and colonizing (D) grew faster and reached a higher final cell density in the presence of nutrients (PC3G). Elongation strategy (PC2S) was also associated to cell morphology and performance: strains producing more dispersal morphs and presenting a greater variance in elongation were characterized by small elongated cells (versus big round) at carrying capacity (E) and an overall inferior performance of growth in nutrient-rich medium (F). Strains presenting a K strategy of growth in nutrient-rich medium presented superior abilities to survive and elongate under starvation conditions (G). Note: because the variables whose correlation we studied do not come from the same experiment, each inset indicates the mean Spearman's correlation coefficient computed over 1000 random associations of replicates (r), the proportion of these associations where correlation was significant (s) and the probability to obtain this proportion by chance (P). See text for detail. Graphs, however, display means of five random associations to illustrate within strain variation.

The above associations were well represented by the first two axes of a final principal component analysis (Figure 5) performed on the seven variables summarizing the four experiments (dispersal rate, colonization probability, survival and elongation in presence of nutrients PC1G, contrast between growth rate and final cell density PC2G, growth performance PC3G, survival and elongation under starvation PC1S, elongation strategy PC2S; Figure 4). The first axis of this comprehensive PCA explained on average 39% of the variance (SD: 1.6%; range: 34% to 44%, over 1000 random associations of replicates across experiments) and the second axis 23% of the variance (SD: 1.4%; range: 19% to 28%). These two axes allowed discriminating strains very efficiently, with only on average 16% (SD: 7.7%; minimum, median, maximum: 0%, 16%, 37%, respectively) of the replicates not correctly classified. Some strains (e.g. 7 and D2) were always perfectly discriminated, and misclassifications only concerned a single replicate for the majority of the strains; only B strain was more frequently misclassified for all replicates.

thumbnailFigure 5. Summary principal component plot representing the correlations between the seven variables summarizing the four experiments. Graph displays means of five random associations to illustrate within strain variation; see text for detail.

Discussion

Genetic variation in dispersal and life history strategies

Our study provided strong evidence for genetic variation in dispersal and life history strategies in Tetrahymena thermophila protozoans. Strains varied considerably and significantly in nearly all aspects of life history (growth performance, survival, cell morphology changes, and degree of phenotypic plasticity of cell shape) under different environmental conditions, and showed significant differences in the associations of these traits under our experimental conditions (axenic food resource or food free medium). Life history trait differences were in many cases significantly related with dispersal strategies, which also varied significantly between strains. Among-strain differences in life-history traits and trait associations were not just due to variation in strain quality, although some T. thermophila strains did perform better than others as regarded a series of core fitness traits (survival, dispersal, growth performance, and colonization ability). This is because strains also differed with respect to various trait associations that cannot as such be considered to represent differential quality but rather different strategies. Examples are the observed among-strain differences in growth rates r versus final population density in presence of nutrients (PC2G), cell shape versus size in presence of nutrients (PC1G), and the opposition between a high degree of phenotypic plasticity of cell shape (high variance in cell elongation) with some cells turning into dispersal morphs versus a more durable and greater, but less plastic, mean cell elongation under starvation (PC2S). The T. thermophila ten strains hence differed also with respect to the patterns of their investment along different trait axes. Hence, our study adds to the limited evidence (e.g. [24,38,39,41-47,58]) for genetic variance in dispersal syndromes within a species. Also, the among-strain variance in the degree of plasticity of cell shape found in T. thermophila is similar to recent findings of genetic variance in plasticity for core life-history traits in other animals, e.g. Caenorhabditis elegans nematodes [55] and in dispersal related morphological traits in pea aphids (a winged-non-winged polyphenism) [47].

The life-history trait association differences found among T. thermophila strains should be attributed to among-strain differences in genes of the macronucleus (somatic nucleus characteristic for ciliated protozoans; reviews by [75,76]), not directly the micronucleus (germline nucleus, not transcribed [77]). This is because macronucleus genes represent a modified, often variably amplified and rearranged subset of micronucleus genes [75,76,78]. The genetic trait variance expressed among clonal lines is therefore not necessarily the same that would be found in crossing experiments (sexual reproduction; involving micronuclei). Another curiosity of ciliates, the amitotic divisions of macronuclei making daughter cells from a clonal division potentially have slightly different alleles or amounts of alleles ([78], reviews by [75,76]) did not invalidate our study. Such differences among cells within a clone line would only make our tests for genetic differences among strains more conservative. Equally, non-genetic 'maternal' effects are not likely to have contributed importantly to among-strain differences and so confounded our study, because all strains were maintained under the same conditions for some 700 generations before the start of experiments (> one year). Finally, the observed strain differences in morphological responses to starvation were very unlikely to be artefactual. We used standard techniques (centrifugation, decanting of super-natant and addition of water) to eliminate nutrients from the medium [e.g. [62,69]], and any damage during centrifugation would rather lead to ciliae loss and so to a decreased, not a greatly enhanced swimming speed. Cells, moreover, started swimming within seconds after centrifugations indicating perfectly preserved cell integrity, verified also under the microscope (E.J.F and N.S. personal observations). Also, all strains were subjected to precisely the same manipulations (same centrifugal forces and length of centrifugation).

Covariation of dispersal strategies with life history traits

Dispersal strategies, growth performance, and survival

The association of T. thermophila dispersal rates in presence of nutrients with other life history traits fits partly with the classic 'colonizer syndrome' [30] envisioned typical for species living in variable patches in metapopulations (see also [21,23,31]). Strains with great short-distance dispersal rates in presence of nutrients had a better growth performance, and high colonization abilities as predicted ([30,31,79] but see [80]), but contrary to classical predictions they also had a better mean strain survival rate than did less dispersive strains. An actual negative co-variance (trade-off) between two life-history traits could be masked by differences among strains in overall condition, because some strains would have more resources to invest both in growth and in survival [28,29]. Our experiments were not set up to explicitly explore this issue. All our strains were, however, kept in the same conditions of nutrient availability in the growth experiments, and so should have had the same amount of resources available. Also we note that our findings could be explained by some alternative dispersal evolution models that do predict positive associations of survival with dispersal rates and reproductive effort for ruderal species, the precise expectations depending on landscape characteristics (productivity, demographics, spatio-temporal structure [34], review by [23]), that are only little known for T. thermophila.

The opposite end of the spectrum, namely T. thermophila strains dispersing little, growing less, and colonizing less well as single cells in presence of food, and surviving food-stress poorly, represented relative philopatry as long as the environment remained good, and likely long-distance dispersal (via fast-swimming dispersal morphs) when the environment turned bad. Such a long-distance colonizer strategy could allow coexisting with better competitor strains (see above). Alternatively, 'philopatric' strains may come from habitats with less local spatio-temporal variability (see e.g. [23,32]) but more catastrophic patch degradations (favouring the production of the dispersal morph). Finally, the variation among strains in dispersal strategies and life history could represent non-adaptive variance around an adaptive mean. In all cases, relative philopatry in presence of food did not signify an adaptation to a more intensive exploitation of already colonized patches (see e.g. [81,83,84,86]), because the growth performance of philopatric strains was inferior and also philopatry was not related significantly with r versus K strategies. Nevertheless, the lower dispersal rate of philopatric strains suggests that they freely tolerate higher densities without dispersing massively. Density-dependent dispersal strategies are indeed predicted by several theoretical models and have been reported in a number of species (e.g. [4,85-88]).

Potential density-dependence of dispersal strategies is particularly interesting given the primitive form of cooperation exhibited by T. thermophila. Cells secrete certain chemical compounds ('growth substances') that increase the survival of other cells [72,73] and also actively form aggregations [74], which should enhance the effective concentration of these substances. Those T. thermophila strains that disperse little even when at high density and are less good at colonizing as single cells, could be strains adapted to and dependent on greater cooperation than the more dispersive strains. Consistent with this, strains producing many dispersal morphs (generally philopatric strains) have smaller cell sizes (at carrying capacity; Figure 4E), and so each cell likely has fewer resources than is the case in more dispersive strains. The inferior rates of survival under starvation for philopatric strains corroborate this idea. If our hypothesis is correct, the philopatric strains should show a stronger degree of aggregation and only disperse much when densities are very high and benefits of 'cooperation' outweighed by costs of conspecific competition. We are currently testing this.

Dispersal morphs

Morphological differences between dispersers and non-dispersers were found in our study, as in a series of other species ([41,58,79], see also [10,11]). Morph diversity was, however, greater than suggested by previous morphological studies on T. thermophila [69]. Dispersing cells were more elongated than non-dispersers in our two-patch dispersal experiment, as expected, and strains that elongated much in presence of nutrients also showed greater mean cell elongation under starvation (but a lower variance and hence lower plasticity). Elongated disperser cells, however, differed in shape from the more elongated cells produced under starvation, and strains dispersing much in presence of nutrients surprisingly produced only few of the fast-swimming dispersal morphs described by Nelsen and Debault [69]. Importantly, the differential production of dispersive morph types by different strains was linked to variation in core life history trait associations (growth, survival, colonization, short-distance dispersal), with strains performing poorly under normal conditions producing more fast-swimming dispersal morphs under starvation.

That the fast-swimming very elongated dispersal morphs likely constitute a commitment to long-distance dispersal in degraded environments was supported by our study. Firstly, our dispersal experiment showed that normal cell morphs or partially elongated cells were very capable of short-distance dispersal [69,70]. Secondly, our starvation experiment confirmed that considerable time is required for a cell to transform itself into the very elongated dispersal morph (4–8 hours in this study, which corresponds to half a generation; see temperature dependent data in [69,70]). Hence transformation, likely also energetically costly, should only occur when cells need to surpass the mobility capacities of normal or partially elongated morphs. Strains producing many fast-swimming morphs should therefore benefit from a dispersal-distance advantage, which may help explain the coexistence of low performance strains with strains doing well (see-colonization-competition models for coexistence [81], see also [82] and references in these).

Other selective factors could also have contributed to the observed variation among strains in dispersal distance strategies. For example different strains may come from habitats that vary with respect to the selective pressures of kin competition versus spatio-temporal variation in habitat quality, distributions of resources and kin in space, or the shapes of distance-dependent dispersal cost functions (e.g. [17,21,81,82,86,87,89]; see also [85]). Finally, we cannot reject that strain variation in life history traits could represent non-adaptive variation around an adaptive mean strategy. Intimate knowledge on the habitats of origin of different strains or, better still, experimental evolution studies will be required to test this.

Conclusion

We found strong evidence for genetic variation in dispersal syndromes in T. thermophila protozoans, because strains differed in overall life-history associations, with dispersal strategies, colonization capacity, survival, reproductive performance, and cell shape plasticity showing complex patterns of variation. While some strains fit rather well with the classic colonizer syndrome (high short-distance dispersal rates, high colonization capacity, and elevated growth performance), they were also characterized by good survival abilities and produced few of the putative long-distance dispersal morphs when subjected to environmental degradation. Poorly performing, locally philopatric strains, by contrast, produced relatively many of these fast-swimming dispersal morphs, and so likely benefit from a dispersal-distance advantage which may facilitate their persistence. Finally, the smaller cell size of these latter strains at carrying capacity and their poorer skills at colonizing as individual cells, suggest that they may be adapted to greater levels of dependency on clone-mate cells (stronger sociality). Overall, differential exposure to selection on competitive and cooperative abilities, in conjunction with selective factors targeting specifically dispersal distance, likely contributed importantly to shaping T. thermophila dispersal and life history evolution. Ongoing studies on aggregation behavior and density-dependence of dispersal strategies will explore this further.

Methods

Culture conditions, and strain conditioning

Stocks of 10 T. thermophila strains (Table 1) were maintained in an axenic (bacteria-free) culture medium (2% proteose peptone, Difco no. 3; 0.2 % yeast Difco; demineralized water; all autoclaved after mixing) (see also [62]), and new cultures were made every three to five weeks by transferring 500 μl of an old stock to a 50 ml tube filled with 40 ml of new culture medium. All work was done under sterile conditions under a flow-hood. Culture tubes were kept at 27°C at a 12 h light/12 h dark cycle in an incubator with sun-glow bulbs emitting light corresponding precisely to natural sunlight.

Table 1. Tetrahymena thermophila strains used

Data collection using digital imaging

A procedure based on image analysis was developed to collect data from digital pictures taken under the microscope. For each sample to be measured, aliquots of 10 μl were taken by pipette after vortexing for 10 seconds at medium speed to homogenize the cell suspension. Each aliquot was then placed in a chamber of a Plexiglas cell count slide (KOVA Glasstic Slide 10 no grids from Hycor Biomedical Inc., California, USA, reference 87146). One grey level picture was taken per chamber using a digital camera (Nikon Coolpix 4500; settings: focal length at 8 mm, fixed exposure at 1/60 s, f/6, manual focus fixed to infinity with focus obtained manually using the microscope) mounted via an eyepiece (Olympus Coolpix digital coupler) on a microscope (Olympus CX-41; settings: 10 × oculars, 4x/013PhL planar lens, condensator placed on phase 2, maximum light); an example of such a picture is shown in Figure 6. Conversion of measurement units on pictures (pixels) into real units (μm and μm2) were obtained from a picture made with the same setup of a similar cell count slide with a calibrated grid of known size and volume (KOVA Glasstic Slide 10 with counting grids, reference 87144).

thumbnailFigure 6. A) Typical digital picture of Tetrahymena thermophila cells used to extract quantitative variables. B) Cells (D2 strain) in normal growth condition (with nutrients) showing large size and a largely rounded shape. C) Cells (20 strain) eight hours after nutrient removal, showing reduced size and (middle cell) the typical elongated shape of the fast-swimming dispersal morph described by Nelsen and Debault [69].

Pictures were analyzed using the public domain image-analysis program Scion Image [92] with a macro that we developed to automate the treatment of large numbers of pictures. The macro first converted each grey level pixel into black or white, with the threshold set at 200 as it proved to optimize discrimination of cells from dust and scratches. Then the "Analyze Particles" function was used to export to a text file a list of parameters for each particle present on the picture. Particles smaller than 35 or bigger than 500 pixels were discarded because these sizes are far outside normal sizes for T. thermophila cells [62] and hence represented dust. For each cell, the following variables were measured: size, color density (mean, mode and std; using the grey level picture), position (X and Y coordinates), perimeter, and three variables from the best fitting ellipse: major axis, minor axis and angle (angle between the major axis and a line parallel to the x-axis of the image). A shape variable was computed as the ratio major/minor axis for each cell; a perfect circle had shape = 1 and shape increased for more elongated cells. The density of cells per picture was obtained from the number of particles identified.

To ensure that results were comparable between all the pictures, the method was designed using fixed parameters only, both for taking the pictures (microscope and camera parameters) and analyzing them (Scion Image macro parameters) with settings described above. Extensive fine-tuning of each step was performed to define an optimal method, and results of picture analyses were manually checked on several hundreds of cells to verify the validity of measures recorded. Artefacts (missed cells, or dust considered as cells) were less than 5% of the results and distributed randomly among pictures; therefore they did not bias the results even if they did increase the random noise in the data.

Experimental design

For each of our 10 T. thermophila strains (Table 1), we studied (1) growth from low density in presence of nutrients, (2) survival in a medium devoid of nutrients (starvation), (3) dispersal in presence of nutrients, and (4) single cell colonization capacity in presence of nutrients. Experiments were started from cell cultures inoculated one week before. 10 pictures of the initial tube were taken and density of cells was measured so that experiments could be started with a precise cell density by diluting the solution accordingly. All experimental setups were autoclaved before experiments and conserved in the incubator under the same conditions as strain stocks (see above).

Growth from low density in presence of nutrients

Growth rate and carrying capacity in nutrient-rich conditions were examined by placing approximately 50 000 cells in a 12 ml polypropylene culture tube (Greiner Bio-One), adding nutrient medium till a total volume of 5 ml to reach a starting cell density of 10 000 cells/ml. Three replicate tubes were made per strain. Pictures were taken at six different times (from 0 to 192 h).

From the raw cell measurements at the replicate level, we estimated five major variables describing cell population growth: carrying capacity K and growth rate r, and final (i.e. at carrying capacity) cell density, size and shape (Table 2). Growth was estimated by fitting the following logistic growth function to successive densities Dt over time t (t > 0):

Table 2. List of variables used in statistical analyses

<a onClick="popup('http://www.biomedcentral.com/1471-2148/7/133/mathml/M1','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1471-2148/7/133/mathml/M1">View MathML</a>

with D0 being the initial density at time 0, and two parameters to be estimated: growth rate r and carrying capacity K. The fitting was achieved through least-squares non linear programming using PROC NLP in SAS [93]. Statistical difference between replicates and strains were tested by comparing a model with equality constraints of K and r parameters between replicates of each strain, a model with equality constrains of K and r parameters between strains, and a model where parameters K and r were free to vary between strains and replicates; model selection was achieved using the AIC criterion [94-96].

Survival in a medium devoid of nutrients (starvation)

To remove all traces of nutrients from cell cultures, the contents of five tubes per strain were concentrated and washed with water (similar to [70]). This was done by centrifugating tubes at 2000 r.p.m. for three minutes, after which the supernatant (approximately 35 ml) was aspirated using a vacuum pump. The pellet, containing the cells, was then resuspended by adding 35 ml of autoclaved demineralized water and then tubes were vortexed gently for 5 seconds to detach any cells sticking to tube walls. This process was repeated four times for each tube. After the fifth centrifugation, the cells were left in a volume of 5 ml (i.e. no water was added) but still vortexed gently. Hence, this washing process had the effect of diluting the nutrient solution originally present in the tube 4096-fold, effectively leaving the cells nutrient-deprived. The five tubes of concentrated cells per strain were then pooled into a single initial tube. A volume containing approximately 200 000 cells was placed in a 12 ml polystyrene culture tube and autoclaved demineralized water was added to increase the total volume to 2 ml. Three replicate tubes were set up for each strain. For each replicate, pictures were taken at 14 time steps up to 408 h (every four hours during the first 24 h, this delay increasing to 24 h and 48 h in the course of the experiment).

From the raw cell measurements at the replicate level, we estimated five major variables describing the changes in cell density and shape: survival as a density sum over time, mean and variance of maximal cell elongation, elongation persistence as a sum of cell elongation over time, and the frequency of putative disperser morphs after 8 h, time when the maximal elongation was observed (Table 2; see also Figure 6 and Additional file 1).

Additional file 1. Fast, directional swimming behavior of dispersal morphs compared to other cells. This video clip shows the rapid, directional swimming behavior of T. thermophila dispersal morphs compared to other cells.

Format: WMV Size: 2.5MB Download fileOpen Data

Dispersal in presence of nutrients

Dispersal rate was measured using a two patch system setup consisting of two tubes (1.5 ml polypropylene test tubes, Greiner Bio-One) connected by a small horizontal pipe (silicone Tygon tubing, 6 mm diameter, 17 mm long) inserted through a hole drilled on the side of each tube. For each setup, a cell suspension volume corresponding to 300 000 cells of a given strain was taken from a culture tube (the initial tube) and introduced in one patch (the start tube), by the use of a pipette. By pipette we then filled up the target tube with cell-free nutrient medium, filling up in this way also the connecting tube and the remaining free space in the start tube. Both start and target tube contained 1 ml in the filled-up state. The setup was then incubated at 24°C. Pilot experiments using an ink solution instead of a cell suspension had verified that filling up the setup and placing it in the incubator did not induce perturbations that could displace cells from the start patch.

After 17 h we estimated the number of cells that had swum from the start to the target patch. This was done by closing both ends of the connecting tube with clamps and then pipetting off the solution in both start and target tubes and taking pictures for the two tubes. The experiment was performed over several days in a semi-block design due to constraints on manipulation time. A total of six replicates were done for each strain, in the manner of three replicates per strain for two out of four blocks (dates).

From the raw cell measurements, we estimated three major variables at the replicate level: initial cell shape, dispersal rate and cell elongation (Table 2).

Single cell colonization capacity in presence of nutrients

To assess colonization capacity, individual cells were isolated from a one-week old cell culture tube using a handheld 10 μl pipetteman and a binocular dissection microscope. Each cell was placed separately in a 12 ml test tube with nutrient medium (the new patch to be colonized), and 10 replicates (10 different tubes) were created per strain. The entire set-up was replicated one week later, leading to a total of 20 single-cell tube experiments per strain. In the second experiment, however, some tubes became infected and were discarded. All tubes were incubated at standard temperatures in the incubator for eight days, after which the success of single cells to colonize the new patches by surviving and undergoing multiple cell divisions was determined.

To determine colonization success we took out a 10 μl aliquot from each experimental tube, after gentle vortexing, placed it on a slide and counted the number of cells directly under binoculars. We then assigned a 1 or a 0 to each experimental tube according to whether any cells were found or not (presence of at least one cell versus none), and calculated (over replicate tubes) the probability of successful colonization for each strain.

Statistical analyses

The basic data for studying life-history trait covariation between strains were the estimates at the replicate level of 14 parameters (Table 2). Covariation between these parameters was studied through Spearman's correlation [97] and principal component analysis [98], and differences between strains by generalized linear models [99] and discriminant analyses [98], all implemented using SAS software.

Because replicates of a given strain were not linked between experiments, we used a permutation at the replicate level procedure to correlate parameters from different experiments. This procedure prevented discarding the information on variation between replicates of a given strain. The replicates of a given strain were randomly associated across experiments 1000 times, and a correlation was computed for each random association. The mean Spearman's correlation (r) was reported as covariation measure between the two parameters studied. Statistical significance of this correlation (showing that it differed significantly from zero), however, was based on the probability of obtaining the observed proportion of significant correlations by chance (in the 1000 simulations). The distribution of this proportion of significant correlations under the null hypothesis of no correlation was also obtained by resampling, with observed values randomized for replicate and strain, breaking any existing correlation. In total, 1000 such random associations were done and the proportion of significant correlations computed. This procedure was repeated another 1000 times, and the P-value for the test was computed as the proportion of random associations with a proportion of significant correlations greater than or equal to the one observed.

Authors' contributions

The work presented here was carried out in collaboration between all authors. EJF, NS and JC defined the research theme. EJF and NS designed methods and experiments, carried out the laboratory experiments, analyzed the data, interpreted the results and wrote the paper. PM and AG co-designed the dispersal and colonization experiments, and co-worked on associated data collection and their interpretation. JC co-designed experiments, discussed analyses, interpretation, and presentation. All authors have contributed to, seen and approved the manuscript.

Acknowledgements

We thank P. Doerder for very kindly donating strains, M. van Baalen for endorsing E.J.F.'s extended stay at Université Pierre et Marie Curie, and two anonymous reviewers for constructive comments. E.J.F. was supported by the European Research Training Network ModLife (Modern Life-History Theory and Its Application to the Management of Natural Resources), funded through the Human Potential Programme of the European Commission (Contract HPRN-CT-2000-00051). N.S. is a postdoctoral researcher of the Belgian Fund for Scientific Research (Chargé de Recherches F.R.S.-FNRS) and acknowledges a travel grant from the same institution (FRFC 2.4587.02F). P.M. was supported by a M.Sc. stipend from the Centre Regional des Oeuvres Universitaires et Scolaires de Paris. J.C. acknowledges support from Agence Nationale de Recherche (ANR-05-BLAN-0265-02).

References

  1. Boulinier T, McCoy KD, Sorci G: Dispersal and parasitism. In Dispersal. Edited by Clobert J, Danchin E, Dhondt AA, Nichols D. Oxford: Oxford University Press; 2001:169-179. OpenURL

  2. Clobert J, Danchin E, Dhondt AA, Nichols D, Editors: Dispersal. Oxford: Oxford University Press; 2001.

  3. Bullock JM, Kenward RE, Hails RS, Editors: Dispersal ecology. Oxford: Blackwell;; 2002.

  4. Clobert J, Ims RA, Rousset F: Causes, mechanisms and consequences of dispersal. In Ecology, genetics and evolution of metapopulation. Edited by Hanski I, Gaggiotti OE. Amsterdam: Academic Press; 2004:307-335. OpenURL

  5. Bowler DE, Benton TG: Causes and consequences of animal dispersal strategies: relating individual behaviour to spatial dynamics.

    Biol Rev 2005, 80:205-225. PubMed Abstract | Publisher Full Text OpenURL

  6. Telschow A, Hammerstein A, Werren JH: The effect of Wolbachia versus genetic incompatibilities on reinforcement and speciation.

    Evolution 2005, 59:1607-1619. PubMed Abstract OpenURL

  7. Yukilevich R, True JR: Divergent outcomes of reinforcement speciation: the relative importance of assortative mating and migration modification.

    Am Nat 2006, 167:638-654. PubMed Abstract | Publisher Full Text OpenURL

  8. Hanski I, Gilpin ME, Editors: Metapopulation Biology, Ecology, Genetics, and Evolution. San Diego: Academic Press; 1997.

  9. Hanski I, Gaggiotti OE, Editors: Ecology, genetics, and evolution of metapopulations. Amsterdam: Academic press; 2004.

  10. Bélichon SJ, Clobert J, Massot M: Are there differences in fitness components between philopatric and dispersing individuals?

    Acta Oecol 1996, 17:503-517. OpenURL

  11. Murren CJ, Julliard R, Schlichting CD, Clobert J: Dispersal, individual phenotype, and phenotypic plasticity. In Dispersal. Edited by Clobert J, Danchin E, Dhondt AA, Nichols D. Oxford: Oxford University Press; 2001:261-282. OpenURL

  12. Keller L, (ed): Queen number and sociality in the insects. Oxford Science Publications, Oxford; 1993.

  13. Van Baalen M, Rand DA: The unit of selection in viscuous populations and the evolution of altruism.

    J Theor Biol 1998, 193:631-648. PubMed Abstract | Publisher Full Text OpenURL

  14. Koella J: The spatial spread of altruism versus the evolutionary response of egoists.

    Proc R Soc Lond B 2000, 267:1979-1985. Publisher Full Text OpenURL

  15. Perrin N, Lehmann L: Is sociality driven by the costs of dispersal or the benefits of philopatry ? A role for kin-discrimination mechanisms.

    Am Nat 2001, 158:471-487. Publisher Full Text OpenURL

  16. West SA, Pen I, Griffin AS: Cooperation and competition between relatives.

    Science 2002, 296:72-75. PubMed Abstract | Publisher Full Text OpenURL

  17. Le Galliard J-F, Ferriere R, Dieckmann U: Adaptive evolution of social traits: origin, trajectories, and correlations of altruism and mobility.

    Am Nat 2005, 165:206-224. PubMed Abstract | Publisher Full Text OpenURL

  18. Olivieri I, Gouyon P-H: Evolution of migration rate and other traits: the metapopulation effect. In Metapopulation Biology, Ecology, Genetics and Evolution. Edited by Hanski I, Gilpin ME. San Diego: Academic Press; 1997:293-324. OpenURL

  19. Ronce O, Olivieri I: Evolution of reproductive effort in a metapopulation with local extinctions and ecological succession.

    Am Nat 1997, 150:220-249. Publisher Full Text OpenURL

  20. Gandon S, Michalakis Y: Multiple causes of the evolution of dispersal. In Dispersal. Edited by Clobert J, Danchin E, Dhondt AA, Nichols D. Oxford: Oxford University Press; 2001:155-167. OpenURL

  21. Ronce O, Olivieri I, Clobert J, Danchin E: Perspectives on the study of dispersal evolution. In Dispersal. Edited by Clobert J, Danchin E, Dhondt AA, Nichols D. Oxford: Oxford University Press; 2001:341-357. OpenURL

  22. Weitkamp LA: The effects of predation on dispersal. In Dispersal. Edited by Clobert J, Danchin E, Dhondt AA, Nichols D. Oxford: Oxford University Press; 2001:180-188. OpenURL

  23. Ronce O, Olivieri I: Life history evolution in metapopulations. In Ecology, genetics and evolution of metapopulation. Edited by Hanski I, Gaggiotti OE. Amsterdam: Academic Press; 2004:227-257. OpenURL

  24. Roff DA, Stirling G, Fairbairn DJ: The evolution of threshold traits: A quantitative genetic analysis of the physiological and life-history correlates of wing dimorphism in the sand cricket.

    Evolution 1997, 51:1910-1919. Publisher Full Text OpenURL

  25. Stearns SC: The Evolution of Life Histories. Oxford University Press, Oxford; 1992.

  26. Roff DA: Life history evolution. Sunderland, MA; Sinauer Associates; 2002.

  27. Malausa T, Guillemaud T, Lapchin L: Combining genetic variation and phenotypic plasticity in tradeoff modeling.

    Oikos 2005, 110:130-338. Publisher Full Text OpenURL

  28. Roff DA, Fairbairn DJ: The evolution of trade-offs: where are we?

    J Evol Biol 2007, 20:433-447. PubMed Abstract | Publisher Full Text OpenURL

  29. Van Noordwijk AJ, deJong G: Acquisition and allocation of resources: their influence on variation in life history tactics.

    Am Nat 1986, 128:137-142. Publisher Full Text OpenURL

  30. Baker HG, Stebbins GL: The genetics of colonizing species. Academic Press, New York; 1965.

  31. Grime JP: Evidence for the existence of three primary strategies in plants and its relevance to ecological and evolutionary theory.

    Am Nat 1977, 111:1194. Publisher Full Text OpenURL

  32. McPeek MA, Holt RD: The evolution of dispersal in spatially and temporally varying environments.

    Am Nat 1992, 140:1010-1027. Publisher Full Text OpenURL

  33. Jansen VAA, Mulder GSEE: Evolving biodiversity.

    Ecol Letters 1999, 2:379-386. Publisher Full Text OpenURL

  34. Ronce O, Perret F, Olivieri I: Landscape dynamics and evolution of colonizers syndromes: Interactions between reproductive effort and dispersal in a metapopulation.

    Evol Res 2000, 14:233-260. OpenURL

  35. Crowley PH, McLetchie DN: Trade-offs and spatial life-history strategies in classical metapopulations.

    Am Nat 2002, 159:190-208. Publisher Full Text OpenURL

  36. Gu H, Danthanarayana W: Quantitative genetic analysis of dispersal in Epiphyas postvittana I. Genetic variation in flight capacity.

    Heredity 1992, 68:53-60. OpenURL

  37. Li J, Margolies DC: Quantitative genetics of aerial dispersal behaviour and life-history traits in Tetranychus urticae .

    Heredity 1993, 70:544-552. OpenURL

  38. Trefilov A, Berard J, Krawczak M, Schmidtke J: Natal dispersal in rhesus macaques is related to serotonin transporter gene promoter variation.

    Behav Genet 2000, 30:295-301. PubMed Abstract | Publisher Full Text OpenURL

  39. Krackow S: Motivational and heritable determinants of dispersal latency in wild male house mice (Mus musculus musculus).

    Ethology 2003, 109:671-689. Publisher Full Text OpenURL

  40. Massot M, Huey RB, Tsuji J, Van-Berkum FH: Genetic, prenatal and postnatal correlates of dispersal in hatchling fence lizards (Sceloporus occidentalis).

    Behav Ecol 2003, 14:650-655. Publisher Full Text OpenURL

  41. Sinervo B, Clobert J: Morphs, dispersal behavior, genetic similarity and the evolution of cooperation.

    Science 2003, 300:1949-1951. PubMed Abstract | Publisher Full Text OpenURL

  42. Bargum K, Boomsma JJ, Sundström L: A genetic component to size in queens of the ant, Formica truncorum .

    Behav Ecol Sociobiol 2004, 57:9-16. Publisher Full Text OpenURL

  43. Pasinelli G, Schiegg K, Walters JR: Genetic and environmental influences on natal dispersal distance in a resident bird species.

    Am Nat 2004, 164:660-669. PubMed Abstract | Publisher Full Text OpenURL

  44. Sinervo B, Calsbeek R, Comendant T, Both C, Adamopoulou C, Clobert J: Genetic and maternal determinants of effective dispersal: the effect of sire genotype and size at birth in side-blotched lizards.

    Am Nat 2006, 168:88-99. PubMed Abstract | Publisher Full Text OpenURL

  45. Lushai G, Loxdale HD, Brookes CP, von Mende N, Harrington R, Hardie J: Genotypic variation among different phenotypes within aphid clones.

    Proc R Soc Lond B 1997, 1382:725-730. Publisher Full Text OpenURL

  46. Osborne KA, Robichon A, Burgess E, Butland S, Shaw RA, Coulthard A, Pereira HS, Greenspan RJ, Sokolowski MB: Natural behavior polymorphism due to cGMP-dependent protein kinase of Drosophila .

    Science 1997, 277:874-876. Publisher Full Text OpenURL

  47. Braendle C, Friebe I, Caillaud MC, Stern DL: Genetic variation for an aphid wing polyphenism is genetically linked to a naturally occurring wing polymorphism.

    Proc R Soc Lond B 2005, 1563:657-664. Publisher Full Text OpenURL

  48. Léna J-P, Clobert J, de Fraipont M, Lecomte J, Guyot G: The relative influence of density and kin competition on dispersal in the common lizard.

    Behav Ecol 1998, 9:500-507. Publisher Full Text OpenURL

  49. Schtickzelle N, Mennechez G, Baguette M: Dispersal depression with habitat fragmentation in the bog fritillary butterfly.

    Ecology 2006, 87:1057-1065. PubMed Abstract | Publisher Full Text OpenURL

  50. O'Riain MJ, Jarvis JUM, Faulkes GG: A dispersive morph in the naked mole-rat.

    Nature 1996, 380:619-621. PubMed Abstract | Publisher Full Text OpenURL

  51. Sundström L: Dispersal polymorphism and physiological condition of males and females in the ant Formica truncorum .

    Behav Ecol 1995, 6:132-139. Publisher Full Text OpenURL

  52. de Fraipont M, Clobert J, John-Alder H, Meylan S: Pre-natal stress increases offspring philopatry.

    J Anim Ecol 2000, 69:404-413. Publisher Full Text OpenURL

  53. Dufty AM, Clobert J, Moller A-P: Hormones, developmental plasticity and adaptation.

    Trends Ecol Evol 2002, 17:190-196. Publisher Full Text OpenURL

  54. Aragon P, Clobert J, Massot M: Individual dispersal status influences space use of conspecific residents in the common lizard, Lacerta vivipara.

    Behav Ecol Sociobiol, in press. OpenURL

  55. Gutteling EW, Riksen JAG, Bakker J, Kammenga JE: Mapping phenotypic plasticity and genotype-environment interactions affecting life-history traits in Caenorhabditis elegans.

    Heredity 2007, 98:28-37. PubMed Abstract | Publisher Full Text OpenURL

  56. DeWitt TJ, Sih A, Wilson DS: Costs and limits of phenotypic plasticity.

    Trends Ecol Evol 1998, 13:77-81. Publisher Full Text OpenURL

  57. Via S, Gomulklewicz R, De Jong G, Scheiner SM, Schlichting CD, Van Tienderen PH: Adaptive phenotypic plasticity: consensus and controversy.

    Trends Ecol Evol 1995, 10:212-217. Publisher Full Text OpenURL

  58. Zera AJ, Denno RF: Physiology and ecology of dispersal polymorphism in insects.

    Ann Rev Entomol 1997, 42:207-230.

    Harrison RG: Dispersal polymorphisms in insects. Ann Rev Ecol System 1980, 11: 95–118.

    Publisher Full Text OpenURL

  59. Friedenberg NA: Experimental evolution of dispersal in spatiotemporally variable microcosms.

    Ecol Letters 2003, 6:953-959. Publisher Full Text OpenURL

  60. Jessup CM, Kassen R, Forde SE, Kerr B, Buckling A, Rainey PB, Bohannan BJM: Big questions, small worlds:microbial model systems in ecology.

    Trends Ecol Evol 2004, 19:189-197. PubMed Abstract | Publisher Full Text OpenURL

  61. Nanney DL, McCoy JW: Characterization of the species of the Tetrahymena pyriformis complex.

    Transac Am Microscop Soc 1976, 95:664-682. Publisher Full Text OpenURL

  62. Asai DJ, Forney JD: Tetrahymena thermophila. In Methods Cell Biol. Volume 62. Academic Press, San Diego; 1999. OpenURL

  63. Doerder FP, Gates L, Eberhardt F, Arslanyolu M: High frequency of sex and equal frequency of mating types in natural populations of the ciliate Tetrahymena thermophila .

    Proc Natl Acad Sci USA 1995, 92:8715-8718. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  64. Laakso J, Loytynoja K, Kaitala V: Environmental noise and population dynamics of the ciliated protozoa Tetrahymena thermophila in aquatic microcosms.

    Oikos 2003, 102:663-671. Publisher Full Text OpenURL

  65. Fryxell JM, Smith IM, Lynn DH: Evaluation of alternate harvesting strategies using experimental microcosms.

    Oikos 2005, 111:143-149. Publisher Full Text OpenURL

  66. Frankel J: Cell biology of Tetrahymena thermophila 1999.

    In Tetrahymena thermophila Methods Cell Biol Edited by Asai DJ, Forney JD. 1999, 62:27-125. OpenURL

  67. Doerder FP, Arslanyolu M, Saad Y, Kaczmarek M, Mendoza M, Mita F: Ecological genetics of Tetrahymena thermophila : Mating types, i-Antigens, Multiple Alleles, Epistasis.

    J Eukar Microbiol 1996, 43:95-100. Publisher Full Text OpenURL

  68. Nanney DL, Simon EH: Laboratory and evolutionary history of Tetrahymena thermophila.

    In Tetrahymena thermophila. Methods Cell Biol Edited by Asai DJ, Forney JD . 1999, 62:3-25. OpenURL

  69. Nelsen EA, Debault LE: Transformation in Tetrahymena pyriformis : description of an inducible phenotype.

    J Protozool 1978, 25:113-119. PubMed Abstract OpenURL

  70. Nelsen EA: Transformation in Tetrahymena thermophila : development of an inducible phenotype.

    Dev Biol 1978, 66:17-31. PubMed Abstract | Publisher Full Text OpenURL

  71. Schtickzelle N, Joiris A, Van Dyck H, M Baguette: Quantitative analysis of changes in movement behaviour within and outside habitat in a specialist butterfly.

    BMC Evolutionary Biology 2007, 7:4.

    doi:10.1186/1471-2148-7-4

    PubMed Abstract | BioMed Central Full Text | PubMed Central Full Text OpenURL

  72. Wheatley DN, Christensen S, Schousboe P, Rasmussen L: Signaling in growth and death: adequate nutrition alone may not be sufficient for ciliates. A minireview.

    Cell Biol Internatl 1993, 17:817-823. Publisher Full Text OpenURL

  73. Christensen ST, Sorensen H, Beyer NH, Kristiansen K, Rasmussen L, Rasmussen MI: Cell death in Tetrahymena thermophila : New observations on culture conditions.

    Cell Biol Internatl 2001, 25:509-519. Publisher Full Text OpenURL

  74. Wille JR, Ehret C: Circadian rythm of pattern formation in populations of a free-swimming organism, Tetrahymena .

    J Protozool 1968, 15:792-795. PubMed Abstract OpenURL

  75. Prescott DN: The DNA of ciliated protozoa.

    Microbiol Rev 1994, 58:223-267. OpenURL

  76. Karrer KM: Tetrahymena genetics: two nuclei are better than one.

    In Tetrahymena thermophila. Methods Cell Biol Edited by Asai DJ, Forney JD. 1999, 62:127-186. OpenURL

  77. Mayo KA, Orias E: Further evidence for the lack of gene expression in the Tetrahymena micronucleus.

    Genetics 1981, 98:747-762. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  78. Sonneborn TM: Tetrahymena pyriformis . In Handbook of genetics, plants, plant viruses and protists. Edited by King RC. Plenum, New York; 1974:433-467. OpenURL

  79. Hanski I, Saastamoinen M, Ovaskainen O: Dispersal-related life-history trade-offs in a butterfly metapopulation.

    J Anim Ecol 2006, 75:91-100. PubMed Abstract | Publisher Full Text OpenURL

  80. Baguette M, Schtickzelle N: Negative relationship between dispersal distance and demography in butterfly metapopulations.

    Ecology 2006, 87:648-654. PubMed Abstract | Publisher Full Text OpenURL

  81. Holmes EE, Wilson HB: Running from trouble: long-distance dispersal and the competitive coexistence of inferior species.

    Am Nat 1998, 151:578-586. Publisher Full Text OpenURL

  82. Johst K, Brandl R, Eber S: Metapopulation persistence in dynamic landscapes: the role of dispersal distance.

    Oikos 2002, 98:263-270. Publisher Full Text OpenURL

  83. Lambin X, Aars J, Piertney SB: Dispersal, intraspecific competition, kin competition and kin facilitation: a review of the empirical evidence. In Dispersal. Edited by Clobert J, Danchin E, Dhondt AA, Nichols JD. Oxford, Oxford University Press; 2001:110-122. OpenURL

  84. Poethke HJ, Hovestadt H: Evolution of density- and patch-size dependent dispersal rates.

    Proc R Soc Lond B 2002, 269:637-645. Publisher Full Text OpenURL

  85. Mathias A, Kisdi E, Olivieri I: Divergent evolution of dispersal in a heterogeneous landscape.

    Evolution 2001, 55:246-259. PubMed Abstract OpenURL

  86. Johst K, Doebeli M, Brandl R: Evolution of complex dynamics in spatially structured populations.

    Proc R Soc Lond B 1999, 266:1147-1154. Publisher Full Text OpenURL

  87. Rousset F, Gandon S: Evolution of the distribution of dispersal distance under distance-dependent cost of dispersal.

    J Evol Biol 2002, 15:515-523. Publisher Full Text OpenURL

  88. Travis JMJ, Murrel DJ, Dytham C: The evolution of density-dependent dispersal.

    Proc R Soc Lond B 1999, 266:1471-2954. Publisher Full Text OpenURL

  89. Müller-Landau HC, Levin SA, Keymer JE: Theoretical perspectives on evolution of long-distance dispersal and the example of specialized pests.

    Ecology 2003, 84:1957-1967. Publisher Full Text OpenURL

  90. Culture Collection of Algae and Protozoa [http://www.ccap.ac.uk] webcite

  91. American Type Culture Collection, ATCC The global biosource center [http://www.atcc.org] webcite

  92. National Institutes of Health, Scion Corporation home page [http://www.scioncorp.com] webcite

  93. SAS Institute, documentation [http://support.sas.com/91doc/docMainpage.jsp] webcite

  94. Johnson JB, Omland KS: Model selection in ecology and evolution.

    Trends Ecol Evol 2004, 19:101-108. PubMed Abstract | Publisher Full Text OpenURL

  95. Richards SA: Testing ecological theory using the information-theoretic approach: Examples and cautionary results.

    Ecology 2005, 86:2805-2814. Publisher Full Text OpenURL

  96. Burnham KP, Anderson DR: Model Selection and Multimodel Inference: A Practical Information-Theoretic Approach. 2nd edition. Springer, New York; 2002.

  97. Zar JH: Biostatistical Analysis. Fourth edition. Prentice Hall, Upper Saddle River, New Jersey; 1999.

  98. McGarigal K, Cushman S, Stafford S: Multivariate statistics for wildlife and ecology research. Springer, New York; 2000.

  99. McCullagh P, Nelder JA: Generalized Linear Models. Chapman and Hall, London; 1989.