Breaking the constraint on the number of cervical vertebrae in mammals: On homeotic transformations in lorises and pottos

Mammals almost always have seven cervical vertebrae. The strong evolutionary constraint on changes in this number has been broken in sloths and manatees. We have proposed that the extremely low activity and metabolic rates of these species relax the stabilizing selection against changes in the cervical count. Our hypothesis is that strong stabilizing selection in other mammals is largely indirect and due to associated pleiotropic effects, including juvenile cancers. Additional direct selection can occur due to biomechanical problems (thoracic outlet syndrome). Low metabolic and activity rates are thought to diminish these direct and indirect effects. To test this hypothesis within the primates, we have compared the number of cervical vertebrae of three lorisid species with particularly low activity and metabolic rates with those of more active primate species, including with their phylogenetically closest active relatives, the galagids (bushbabies). In support of our hypothesis, we found that 37.6% of the lorisid specimens had an abnormal cervical count, which is a higher percentage than in the other nine primate families, in which the incidence varied from zero to 2.2%. We conclude that our data support the importance of internal selection in constraining evolvability and of a relaxed stabilizing selection for increasing evolvability. Additionally, we discuss that there is no support for a role of the muscularized diaphragm in the evolutionary constraint.

Lifespan prolonging mechanisms and insulin upregulation without fat accumulation in long-lived reproductives of a higher termite

Kings and queens of eusocial termites can live for decades, while queens sustain a nearly maximal fertility. To investigate the molecular mechanisms underlying their long lifespan, we carried out transcriptomics, lipidomics and metabolomics in Macrotermes natalensis on sterile short-lived workers, long-lived kings and five stages spanning twenty years of adult queen maturation. Reproductives share gene expression differences from workers in agreement with a reduction of several aging-related processes, involving upregulation of DNA damage repair and mitochondrial functions. Anti-oxidant gene expression is downregulated, while peroxidability of membranes in queens decreases. Against expectations, we observed an upregulated gene expression in fat bodies of reproductives of several components of the IIS pathway, including an insulin-like peptide, Ilp9. This pattern does not lead to deleterious fat storage in physogastric queens, while simple sugars dominate in their hemolymph and large amounts of resources are allocated towards oogenesis. Our findings support the notion that all processes causing aging need to be addressed simultaneously in order to prevent it.

Exploring copy number variants in deceased fetuses and neonates with abnormal vertebral patterns and cervical ribs

Background

Cervical patterning abnormalities are rare in the general population, but one variant, cervical ribs, is particularly common in deceased fetuses and neonates. The discrepancy between the incidence in the general population and early mortality is likely due to indirect selection against cervical ribs. The cause for the co‐occurrence of cervical ribs and adverse outcome remains unidentified. Copy number variations resulting in gain or loss of specific genes involved in development and patterning could play a causative role.

Methods

Radiographs of 374 deceased fetuses and infants, including terminations of pregnancies, stillbirths and neonatal deaths, were assessed. Copy number profiles of 265 patients were determined using single nucleotide polymorphism array.

Results

274/374 patients (73.3%) had an abnormal vertebral pattern, which was associated with congenital abnormalities. Cervical ribs were present in 188/374 (50.3%) and were more common in stillbirths (69/128 [53.9%]) and terminations of pregnancies (101/188 [53.7%]), compared to live births (18/58, 31.0%). Large (likely) deleterious copy number variants and aneuploidies were prevalent in these patients. None of the rare copy number variants were recurrent or overlapped with candidate genes for vertebral patterning.

Conclusions

The large variety of copy number variants in deceased fetuses and neonates with similar abnormalities of the vertebral pattern probably reflects the etiological heterogeneity of vertebral patterning abnormalities. This genetic heterogeneity corresponds with the hypothesis that cervical ribs can be regarded as a sign of disruption of critical, highly interactive stages of embryogenesis. The vertebral pattern can probably provide valuable information regarding fetal and neonatal outcome.

Mild drought in the vegetative stage induces phenotypic, gene expression, and DNA methylation plasticity in Arabidopsis but no transgenerational effects

There is renewed interest in whether environmentally induced changes in phenotypes can be heritable. In plants, heritable trait variation can occur without DNA sequence mutations through epigenetic mechanisms involving DNA methylation. However, it remains unknown whether this alternative system of inheritance responds to environmental changes and if it can provide a rapid way for plants to generate adaptive heritable phenotypic variation. To assess potential transgenerational effects induced by the environment, we subjected four natural accessions of Arabidopsis thaliana together with the reference accession Col-0 to mild drought in a multi-generational experiment. As expected, plastic responses to drought were observed in each accession, as well as a number of intergenerational effects of the parental environments. However, after an intervening generation without stress, except for a very few trait-based parental effects, descendants of stressed and non-stressed plants were phenotypically indistinguishable irrespective of whether they were grown in control conditions or under water deficit. In addition, genome-wide analysis of DNA methylation and gene expression in Col-0 demonstrated that, while mild drought induced changes in the DNA methylome of exposed plants, these variants were not inherited. We conclude that mild drought stress does not induce transgenerational epigenetic effects.

Mild drought in the vegetative stage induces phenotypic, gene expression, and DNA methylation plasticity in Arabidopsis but no transgenerational effects

There is renewed interest in whether environmentally induced changes in phenotypes can be heritable. In plants, heritable trait variation can occur without DNA sequence mutations through epigenetic mechanisms involving DNA methylation. However, it remains unknown whether this alternative system of inheritance responds to environmental changes and if it can provide a rapid way for plants to generate adaptive heritable phenotypic variation. To assess potential transgenerational effects induced by the environment, we subjected four natural accessions of Arabidopsis thaliana together with the reference accession Col-0 to mild drought in a multi-generational experiment. As expected, plastic responses to drought were observed in each accession, as well as a number of intergenerational effects of the parental environments. However, after an intervening generation without stress, except for a very few trait-based parental effects, descendants of stressed and non-stressed plants were phenotypically indistinguishable irrespective of whether they were grown in control conditions or under water deficit. In addition, genome-wide analysis of DNA methylation and gene expression in Col-0 demonstrated that, while mild drought induced changes in the DNA methylome of exposed plants, these variants were not inherited. We conclude that mild drought stress does not induce transgenerational epigenetic effects.

Assessing species richness trends: Declines of bees and bumblebees in the Netherlands since 1945

Estimating and predicting temporal trends in species richness is of general importance, but notably difficult because detection probabilities of species are imperfect and many datasets were collected in an opportunistic manner. We need to improve our capabilities to assess richness trends using datasets collected in unstandardized procedures with potential collection bias. Two methods are proposed and applied to estimate richness change, which both incorporate models for sampling effects and detection probability: (a) nonlinear species accumulation curves with an error variance model and (b) Pradel capture-recapture models. The methods are used to assess nationwide temporal trends (1945-2018) in the species richness of wild bees in the Netherlands. Previously, a decelerating decline in wild bee species richness was inferred for part of this dataset. Among the species accumulation curves, those with nonconstant changes in species richness are preferred. However, when analyzing data subsets, constant changes became selected for non-Bombus bees (for samples in collections) and bumblebees (for spatial grid cells sampled in three periods). Smaller richness declines are predicted for non-Bombus bees than bumblebees. However, when relative losses are calculated from confidence intervals limits, they overlap and touch zero loss. Capture-recapture analysis applied to species encounter histories infers a constant colonization rate per year and constant local species survival for bumblebees and other bees. This approach predicts a 6% reduction in non-Bombus species richness from 1945 to 2018 and a significant 19% reduction for bumblebees. Statistical modeling to detect species richness time trends should be systematically complemented with model checking and simulations to interpret the results. Data inspection, assessing model selection bias, and comparisons of trends in data subsets were essential model checking strategies in this analysis. Opportunistic data will not satisfy the assumptions of most models and this should be kept in mind throughout.

Embryonal life histories: Desiccation plasticity and diapause in the Argentinean pearlfish Austrolebias bellottii

Embryos of annual killifish diapause in soil egg banks while ponds are dry. Their rates of development and survival in different developmental stages determine the numbers and stages of embryos at rewetting. In the Argentinean pearlfish Austrolebias bellottii, we investigated plasticity for desiccation in such embryonal life history components across phases of mild desiccation and rewetting and also effects of life history on hatching. In comparison with nonannuals, our data suggest that incidences of diapause have become relatively independent of the occurrence of desiccation, as if they have become genetically assimilated. We found limited survival effects of desiccation, limited developmental delays, and an acceleration of development into the prehatching stage. This response can be adaptive when desiccation informs that an opportunity to hatch approaches. Embryos arrest development in the prehatching stage (diapause DIII) or in the dispersed-cell phase (diapause DI). Parental pair variation in rates of development and survival in the earliest developmental stages affects the fraction of embryos that are in DI at rewetting and the number surviving. Given such effects on life history fitness components, rates during embryonal development seem "visible" to selection and the developmental system can thus adapt when pair variation contains a heritable component. In agreement with expectations for the presence of diversified bet-hedging, some embryos hatched and others not in over half of the clutches with several developed embryos at the moment of rewetting. Hatching probabilities increased for eggs produced later in the experiment, and they increased when embryos were rewetted a second time after two months. This response is opposite of what is expected when age-dependent hatching would be adapted to exploit opportunities for completing another generation before the dry season.

Climate and atmosphere simulator for experiments on ecological systems in changing environments

Grand challenges in global change research and environmental science raise the need for replicated experiments on ecosystems subjected to controlled changes in multiple environmental factors. We designed and developed the Ecolab as a variable climate and atmosphere simulator for multifactor experimentation on natural or artificial ecosystems. The Ecolab integrates atmosphere conditioning technology optimized for accuracy and reliability. The centerpiece is a highly contained, 13-m(3) chamber to host communities of aquatic and terrestrial species and control climate (temperature, humidity, rainfall, irradiance) and atmosphere conditions (O2 and CO2 concentrations). Temperature in the atmosphere and in the water or soil column can be controlled independently of each other. All climatic and atmospheric variables can be programmed to follow dynamical trajectories and simulate gradual as well as step changes. We demonstrate the Ecolab's capacity to simulate a broad range of atmospheric and climatic conditions, their diurnal and seasonal variations, and to support the growth of a model terrestrial plant in two contrasting climate scenarios. The adaptability of the Ecolab design makes it possible to study interactions between variable climate-atmosphere factors and biotic disturbances. Developed as an open-access, multichamber platform, this equipment is available to the international scientific community for exploring interactions and feedbacks between ecological and climate systems.

How to measure maturation: a comparison of probabilistic methods used to test for genotypic variation and plasticity in the decision to mature

Maturation is a developmental trait that plays a key role in shaping organisms' life-history. However, progress in understanding how maturation phenotypes evolve has been held back by confusion over how best to model maturation decisions and a lack of studies comparing genotypic variation in maturation. Here, we fitted probabilistic maturation reaction norms (PMRNs) to data collected from five clones of Daphnia magna and five of Daphnia pulex collected from within and between different populations. We directly compared the utility of modeling approaches that assume maturation to be a process with an instantaneous rate with those that do not by fitting maturation rate and logistic regression models, and emphasize similarities and differences between them. Our results demonstrate that in Daphnia, PMRNs using a logistic regression approach were simpler to use and provided a better fit to the data. The decision to mature was plastic across a range of growth trajectories and dependent upon both body size and age. However, the age effect was stronger in D. magna than D. pulex and varied considerably between clones. Our results support the idea that maturation thresholds can evolve but also suggest that the notion of a threshold based on a single fixed state is an oversimplification that underestimates the adaptability of these important traits.

Phylogenetic reconstruction and shell evolution of the Diplommatinidae (Gastropoda: Caenogastropoda)

The fascinating and often unlikely shell shapes in the terrestrial micromollusc family Diplommatinidae (Gastropoda: Caenogastropoda) provide a particularly attractive set of multiple morphological traits to investigate evolutionary patterns of shape variation. Here, a molecular phylogenetic reconstruction, based on five genes and 2700 bp, was undertaken for this family, integrated with ancestral state reconstruction and phylogenetic PCA of discrete and quantitative traits, respectively. We found strong support for the Diplommatininae as a monophyletic group, separating the Cochlostomatidae into a separate family. Five main clades appear within the Diplommatininae, corresponding with both coiling direction and biogeographic patterns. A Belau clade (A) with highly diverse (but always sinistral) morphology comprised Hungerfordia, Palaina, and some Diplommatina. Arinia (dextral) and Opisthostoma (sinistroid) are sister groups in clade B. Clade C and D solely contain sinistral Diplommatina that are robust and little ornamented (clade C) or slender and sculptured (clade D). Clade E is dextral but biogeographically diverse with species from all sampled regions save the Caroline Islands. Adelopoma, Diplommatina, Palaina, and Hungerfordia require revision to allow taxonomy to reflect phylogeny, whereas Opisthostoma is clearly monophyletic. Ancestral state reconstruction suggests a sinistral origin for the Diplommatinidae, with three reversals to dextrality.

What life cycle graphs can tell about the evolution of life histories

We analyze long-term evolutionary dynamics in a large class of life history models. The model family is characterized by discrete-time population dynamics and a finite number of individual states such that the life cycle can be described in terms of a population projection matrix. We allow an arbitrary number of demographic parameters to be subject to density-dependent population regulation and two or more demographic parameters to be subject to evolutionary change. Our aim is to identify structural features of life cycles and modes of population regulation that correspond to specific evolutionary dynamics. Our derivations are based on a fitness proxy that is an algebraically simple function of loops within the life cycle. This allows us to phrase the results in terms of properties of such loops which are readily interpreted biologically. The following results could be obtained. First, we give sufficient conditions for the existence of optimisation principles in models with an arbitrary number of evolving traits. These models are then classified with respect to their appropriate optimisation principle. Second, under the assumption of just two evolving traits we identify structural features of the life cycle that determine whether equilibria of the monomorphic adaptive dynamics (evolutionarily singular points) correspond to fitness minima or maxima. Third, for one class of frequency-dependent models, where optimisation is not possible, we present sufficient conditions that allow classifying singular points in terms of the curvature of the trade-off curve. Throughout the article we illustrate the utility of our framework with a variety of examples.

Pollination and protection against herbivory of Nepalese Coelogyninae (Orchidaceae)

PREMISE OF THE STUDY

Although many species in the orchid genus Coelogyne are horticulturally popular, hardly anything is known about their pollination. Pollinators of three species were observed in the field in Nepal. This information is urgently needed because many orchid species in Nepal are endangered. Whether the exudates produced by extrafloral nectaries played a role in protection against herbivory was also investigated.

METHODS

Pollinators of C. flaccida, C. nitida, and Otochilus albus were filmed, captured, and identified. Ant surveys and exclusion experiments were carried out. To investigate whether pollinators are needed for fruit set, plants were wrapped in mesh wire bags. Inflorescence stems were examined with microscopy. Fehling's reagent was used to detect sugars in extrafloral exudates.

KEY RESULTS

Coelogyne flaccida and C. nitida need pollinators to set fruit and are pollinated by wild bees identified as Apis cerana. Otochilus albus was found to be pollinated by Bombus kashmirensis. Extrafloral nectar was found to be exuded by nectary-modified stomata and contained high amounts of sugars. Different species of ants were observed collecting these exudates. A significant difference was found in damage inflicted by flower and leaf-eating beetles between C. nitida plants living in trees with ant nests and those in ant-free trees.

CONCLUSIONS

Floral syndromes include scented and colored trap flowers without reward to their pollinators. All orchids investigated exude extrafloral nectar by nectary-modified stomata. This nectar was found to flow from the phloem to the stomata through intercellular spaces in the outer parenchymatous layer of the inflorescence.

Handedness and asymmetry in scale-eating cichlids: antisymmetries of different strength

Individual symmetry is believed to be advantageous and reflecting developmental stability, but frequency-dependent selection can also maintain polymorphisms of asymmetric phenotypes. There are many examples of so-called antisymmetry, where mirror image morphs occur at equal frequencies. With very few exceptions, these are caused by nongenetic variation. One notable exception is handedness and mouth bending variation in the scale-eating cichlid Perissodus microlepis, which has been suggested to be an example of antisymmetry determined by a single genetic locus of large effect. Here, we report that this handedness and mouth bending asymmetry are not jointly and exclusively determined by a single major locus. We found no evidence of a major locus for asymmetry and some support for a major handedness locus. Also, asymmetry is plastic in this species: it can change in adults. We suggest that behavioral handedness in this system precedes and guides morphological asymmetry.

Female mating preference functions predict sexual selection against hybrids between sibling species of cichlid fish

The evolutionary outcome of interspecific hybridization, i.e. collapse of species into a hybrid swarm, persistence or even divergence with reinforcement, depends on the balance between gene flow and selection against hybrids. If female mating preferences are open-ended but sign-inversed between species, they can theoretically be a source of such selection. Cichlid fish in African lakes have sustained high rates of speciation despite evidence for widespread hybridization, and sexual selection by female choice has been proposed as important in the origin and maintenance of species boundaries. However, it had never been tested whether hybridizing species have open-ended preference rules. Here we report the first experimental test using Pundamilia pundamilia, Pundamilia nyererei and their hybrids in three-way choice experiments. Hybrid males are phenotypically intermediate. Wild-caught females of both species have strong preferences for conspecific over heterospecific males. Their responses to F1 hybrid males are intermediate, but more similar to responses to conspecifics in one species and more similar to responses to heterospecifics in the other. We suggest that their mate choice mechanism may predispose haplochromine cichlids to maintain and perhaps undergo phenotypic diversification despite hybridization, and that species differences in female preference functions may predict the potential for adaptive trait transfer between hybridizing species.

Do large dogs die young?

In most animal taxa, longevity increases with body size across species, as predicted by the oxidative stress theory of aging. In contrast, in within-species comparisons of mammals and especially domestic dogs (e.g. Patronek et al., '97; Michell, '99; Egenvall et al., 2000; Speakman et al., 2003), longevity decreases with body size. We explore two datasets for dogs and find support for a negative relationship between size and longevity if we consider variation across breeds. Within breeds, however, the relationship is not negative and is slightly, but significantly, positive in the larger of the two datasets. The negative across-breed relationship is probably the consequence of short life spans in large breeds. Artificial selection for extremely high growth rates in large breeds appears to have led to developmental diseases that seriously diminish longevity.

The Interplay between Behavior and Morphology in the Evolutionary Dynamics of Resource Specialization

We analyze the consequences of diet choice behavior for the evolutionary dynamics of foraging traits by means of a mathematical model. The model is characterized by the following features. Consumers feed on two different substitutable resources that are distributed in a fine-grained manner. On encounter with a resource item, consumers decide whether to attack it so as to maximize their energy intake. Simultaneously, evolutionary change occurs in morphological traits involved in the foraging process. The assumption here is that evolution is constrained by a trade-off in the consumer's ability to forage on the alternative resources. The model predicts that flexible diet choice behavior can guide the direction of evolutionary change and mediate coexistence of different consumer types. Such polymorphisms can evolve from a monomorphic population at evolutionary branching points and also at points where a small genetic change in a trait can provoke a sharp instantaneous and nongenetic change in choice behavior. In the case of weak trade-offs, the evolutionary dynamics of a dimorphic consumer population can lead to alternative evolutionarily stable communities. The robustness of these predictions is checked with individual-based simulations and by relaxing the assumption of optimally foraging consumers.

Extreme selection in humans against homeotic transformations of cervical vertebrae

Why do all mammals, except for sloths and manatees, have exactly seven cervical vertebrae? In other vertebrates and other regions, the vertebral number varies considerably. We investigated whether natural selection constrains the number of cervical vertebrae in humans. To this end, we determined the incidence of cervical ribs and other homeotic vertebral changes in radiographs of deceased human fetuses and infants, and analyzed several existing datasets on the incidence in infants and adults. Our data show that homeotic transformations that change the number of cervical vertebrae are extremely common in humans, but are strongly selected against: almost all individuals die before reproduction. Selection is most probably indirect, caused by a strong coupling of such changes with major congenital abnormalities. Changes in the number of thoracic vertebrae appear to be subject to weaker selection, in good correspondence with the weaker evolutionary constraint on these numbers. Our analysis highlights the role of prenatal selection in the conservation of our common body plan.

Protected polymorphism and evolutionary stability in pleiotropic models with trait-specific dominance

When alleles have pleiotropic effects on a number of quantitative traits, the degree of dominance between a pair of alleles can be different for each trait. Such trait-specific dominance has been studied previously in models for the maintenance of genetic variation by antagonistic effects of an allele on two fitness components. By generalizing these models to an arbitrary number of fitness components or other phenotypic traits with different degrees of dominance, I show that genetic polymorphism is generally impossible without antagonistic fitness effects of different traits and without trait-specific dominance. I also investigate dominance and pleiotropy from a more long-term evolutionary perspective, allowing for the study of general ecological scenarios, and I discuss the effects of trait-specific dominance on evolutionary stability criteria. When selection is mainly directional and only trait-specific dominance and antagonism cause the emergence of polymorphism, then these polymorphisms can be overtaken by single mutants again, such that they are probably short-lived on an evolutionary time scale. Near evolutionarily singular points where directional selection is absent, trait-specific dominance and overdominance facilitate the emergence of polymorphism and cause evolutionary divergence in some cases. An important outcome of these models is that trait-specific dominance allows for the emergence of genetic polymorphisms without a selective disadvantage for heterozygotes. This removes the scope for the evolution of assortative mate choice and affects dominance modification. Sympatric speciation by disruptive ecological selection requires this heterozygote disadvantage in order to evolve, and therefore it becomes less plausible if the emergence of genetic polymorphism usually occurs via trait-specific dominance and antagonistic effects.

Disruptive selection and then what?

Disruptive selection occurs when extreme phenotypes have a fitness advantage over more intermediate phenotypes. The phenomenon is particularly interesting when selection keeps a population in a disruptive regime. This can lead to increased phenotypic variation while disruptive selection itself is diminished or eliminated. Here, we review processes that increase phenotypic variation in response to disruptive selection and discuss some of the possible outcomes, such as sympatric species pairs, sexual dimorphisms, phenotypic plasticity and altered community assemblages. We also identify factors influencing the likelihoods of these different outcomes.

A new perspective on developmental plasticity and the principles of adaptive morph determination

Organisms can have divergent paths of development leading to alternative phenotypes, or morphs. The choice of developmental path may be set by environmental cues, the individual's genotype, or a combination of the two. Using individual-based simulation and analytical investigation, we explore the idea that from the viewpoint of a developmental switch, genetic morph determination can sometimes be regarded as adaptive developmental plasticity. We compare the possibilities for the evolution of environmental and genetic morph determination and combinations of the two in situations with spatial variation in conditions. We find that the accuracy of environmental cues in predicting coming selective conditions is important for environmental morph determination, in accordance with previous results, and that genetic morph determination is favored in a similar way by the accuracy of genetic cues, in the form of selectively maintained gene frequency differences between local populations. Restricted gene flow and strong selection acting on the phenotypic alternatives produce clearer gene frequency differences and lead to greater accuracy of genetic cues. For combined environmental and genetic morph determination, we show that the developmental machinery can evolve toward efficiently combining information in environmental and genetic cues for the purpose of predicting coming selective conditions.

The evolution of resource specialization through frequency-dependent and frequency-independent mechanisms

Levins's fitness set approach has shaped the intuition of many evolutionary ecologists about resource specialization: if the set of possible phenotypes is convex, a generalist is favored, while either of the two specialists is predicted for concave phenotype sets. An important aspect of Levins's approach is that it explicitly excludes frequency-dependent selection. Frequency dependence emerged in a series of models that studied the degree of character displacement of two consumers coexisting on two resources. Surprisingly, the evolutionary dynamics of a single consumer type under frequency dependence has not been studied in detail. We analyze a model of one evolving consumer feeding on two resources and show that, depending on the trait considered to be subject to evolutionary change, selection is either frequency independent or frequency dependent. This difference is explained by the effects different foraging traits have on the consumer-resource interactions. If selection is frequency dependent, then the population can become dimorphic through evolutionary branching at the trait value of the generalist. Those traits with frequency-independent selection, however, do indeed follow the predictions based on Levins's fitness set approach. This dichotomy in the evolutionary dynamics of traits involved in the same foraging process was not previously recognized.

The analysis of reaction norms for age and size at maturity using maturation rate models

Reaction norms for age and size at maturity are being analyzed to answer important questions about the evolution of life histories. A new statistical method is developed in the framework of time-to-event data analysis, which circumvents shortcomings in currently available approaches. The method emphasizes the estimation of age- and size-dependent maturation rates. Individual probabilities of maturation during any given time interval follow by integrating maturation rate along the growth curve. The integration may be performed in different ways, over ages or sizes or both, corresponding to different assumptions on how individuals store the operational history of the maturation process. Data analysis amounts to fitting generalized nonlinear regression models to a maturation status variable. This technique has three main advantages over existing methods: (1) treating maturation as a stochastic process enables one to specify a rate of maturation; (2) age and size at which maturation occurs do not have to be observed exactly, and bias arising from approximations and interpolations is avoided; (3) ages at which sizes are measured and maturation status are observed can differ between individuals. An application to data on the springtail Folsomia candida is presented. Models with age-dependent integration of maturation rates were preferred. The analysis demonstrates a significant size dependence of the maturation rate but no age dependence.

Adaptation and constraint in the evolution of environmental sex determination

When environments differentially influence male and female performance, environmental sex determination (ESD) might evolve. The conclusion from several previous theoretical models was that reaction norms for sex determination should have a single, sharp threshold, with only females being produced in some environments and only males in others. These reaction norms can be disadvantageous in fluctuating environments, however, because they lead to sex-ratio fluctuations. We analysed the evolution of ESD, looking for equilibrium strategies in unconstrained as well as constrained strategy spaces. We identified situations where a single-threshold reaction norm is not evolutionarily stable. In these cases, we found stable strategies in the form of complex reaction norms, showing an oscillatory pattern of sex determination with respect to variation in an environmental variable. Considering that constraints could prevent such phenotypes from being realized, we found that certain randomized reaction norms, with probabilistic sex determination for a range of environments, would achieve nearly the same fitness. We also investigated reaction norms constrained to have a single threshold and found that genetic polymorphism in the environmental threshold value could evolve, producing a similar effect as a randomized reaction norm. We argue that the appearance of genetic variation can be regarded as an alternative outcome when constraints prevent the evolution of a more complex or a randomized strategy.

The evolution of environmental and genetic sex determination in fluctuating environments

Twenty years ago, Bulmer and Bull suggested that disruptive selection, produced by environmental fluctuations, can result in an evolutionary transition from environmental sex determination (ESD) to genetic sex determination (GSD). We investigated the feasibility of such a process, using mutation-limited adaptive dynamics and individual-based computer simulations. Our model describes the evolution of a reaction norm for sex determination in a metapopulation setting with partial migration and variation in an environmental variable both within and between local patches. The reaction norm represents the probability of becoming a female as a function of environmental state and was modeled as a sigmoid function with two parameters, one giving the location (i.e., the value of the environmental variable for which an individual has equal chance of becoming either sex) and the other giving the slope of the reaction norm for that environment. The slope can be interpreted as being set by the level of developmental noise in morph determination, with less noise giving a steeper slope and a more switchlike reaction norm. We found convergence stable reaction norms with intermediate to large amounts of developmental noise for conditions characterized by low migration rates, small differential competitive advantages between the sexes over environments, and little variation between individual environments within patches compared to variation between patches. We also considered reaction norms with the slope parameter constrained to a high value, corresponding to little developmental noise. For these we found evolutionary branching in the location parameter and a transition from ESD toward GSD, analogous to the original analysis by Bulmer and Bull. Further evolutionary change, including dominance evolution, produced a polymorphism acting as a GSD system with heterogamety. Our results point to the role of developmental noise in the evolution of sex determination.