Comparative phenomics: a new approach to study heterochrony

Understanding the links between development and evolution is one of the major challenges of biology. 'Heterochronies', evolutionary alterations in the timings of development are posited as a key mechanism of evolutionary change, but their quantification requires gross simplification of organismal development. Consequently, how changes in event timings influence development more broadly is poorly understood. Here, we measure organismal development as spectra of energy in pixel values of video, creating high-dimensional landscapes integrating development of all visible form and function. This approach we termed 'Energy proxy traits' (EPTs) is applied alongside previously identified heterochronies in three freshwater pulmonate molluscs (Lymnaea stagnalis, Radix balthica and Physella acuta). EPTs were calculated from time-lapse video of embryonic development to construct a continuous functional time series. High-dimensional transitions in phenotype aligned with major sequence heterochronies between species. Furthermore, differences in event timings between conspecifics were associated with changes in high-dimensional phenotypic space. We reveal EPTs as a powerful approach to considering the evolutionary importance of alterations to developmental event timings. Reimagining the phenotype as energy spectra enabled continuous quantification of developmental changes in high-dimensional phenotypic space, rather than measurement of timings of discrete events. This approach has the possibility to transform how we study heterochrony and development more generally.

A phenomics approach reveals interspecific differences in integrated developmental responses to chronic elevated temperatures

Phenomics, high-dimensional organismal phenotyping, is advanced as a solution to quantifying complex developmental responses to elevated temperatures . 'Energy proxy traits' (EPTs) measure the phenotype as a spectrum of energy values across different temporal frequencies from pixel value fluctuations of video. Whilst they have proven effective in measuring the biology of complex and dynamic developing organisms, their utility in assessing environmental sensitivity of different species is untested. Using EPTs, we assess the relative thermal sensitivities of embryos of three species of freshwater snail with marked differences in their developmental event timings . Embryos of Lymnaea stagnalis, Radix balthica and Physella acuta were videoed hourly for the duration of their embryonic development at two temperatures: 20˚C and 25˚C. The video was used to calculate EPTs for the duration of their embryonic development, and during discrete physiological windows in development. Changes in energy spectra during development identified marked differences in thermal sensitivities between species, and suggest a relatively heightened sensitivity of gross rates of embryonic physiology and behaviour in embryos of R. balthica, developmental window specific thermal responses that reflect ontogenetic differences in observable physiologies, and temperature induced changes in physiological event timing. EPTs enabled comparison of high-dimensional spectral phenotypes, providing a unique capability for assessing sensitivity continuously in developing individuals. Such integrative and scalable phenotyping is prerequisite for improved understanding of the sensitivity of early life stages of different species.

Divergence in Thermal Physiology Could Contribute to Vertical Segregation in Intertidal Ecotypes of Littorina saxatilis

AbstractThermal stress is a potentially important selective agent in intertidal marine habitats, but the role that thermal tolerance might play in local adaptation across shore height has been underexplored. Northwest Spain is home to two morphologically distinct ecotypes of the periwinkle Littorina saxatilis, separated by shore height and subject to substantial differences in thermal stress exposure. However, despite other biotic and abiotic drivers of ecotype segregation being well studied, their thermal tolerance has not been previously characterized. We investigated thermal tolerance across multiple life history stages by employing the thermal death time (TDT) approach to determine (i) whether the two ecotypes differ in thermal tolerance and (ii) how any differences vary with life history stage. Adults of the two ecotypes differed in their thermal tolerance in line with their shore position: the upper-shore ecotype, which experiences more extreme temperatures, exhibited greater endurance of thermal stress compared with the lower-shore ecotype. This difference was most pronounced at the highest temperatures tested. The proximate physiological basis for these differences is unknown but likely due to a multifarious interaction of traits affecting different parts of the TDT curve. Differences in tolerance between ecotypes were less pronounced in early life history stages but increased with ontogeny, suggesting partial divergence of this trait during development. Thermal tolerance could potentially play an important role in maintaining population divergence and genetic segregation between the two ecotypes, since the increased thermal sensitivity of the lower-shore ecotype may limit its dispersal onto the upper shore and so restrict gene flow.

Spectral phenotyping of embryonic development reveals integrative thermodynamic responses

Background

Energy proxy traits (EPTs) are a novel approach to high dimensional organismal phenotyping that quantify the spectrum of energy levels within different temporal frequencies associated with mean pixel value fluctuations from video. They offer significant potential in addressing the phenotyping bottleneck in biology and are effective at identifying lethal endpoints and measuring specific functional traits, but the extent to which they might contribute additional understanding of the phenotype remains unknown. Consequently, here we test the biological significance of EPTs and their responses relative to fundamental thermodynamic principles. We achieve this using the entire embryonic development of Radix balthica, a freshwater pond snail, at different temperatures (20, 25 & 30 °C) and comparing responses against predictions from Arrhenius’ equation (Q₁₀ = 2).

Results

We find that EPTs are thermally sensitive and their spectra of frequency response enable effective high-dimensional treatment clustering throughout organismal development. Temperature-specific deviation in EPTs from thermodynamic predictions were evident and indicative of physiological mitigation, although they differed markedly in their responses from manual measures. The EPT spectrum was effective in capturing aspects of the phenotype predictive of biological outcomes, and suggest that EPTs themselves may reflect levels of energy turnover.

Conclusions

Whole-organismal biology is incredibly complex, and this contributes to the challenge of developing universal phenotyping approaches. Here, we demonstrate the biological relevance of a new holistic approach to phenotyping that is not constrained by preconceived notions of biological importance. Furthermore, we find that EPTs are an effective approach to measuring even the most dynamic life history stages.

Supplementary information

The online version contains supplementary material available at 10.1186/s12859-021-04152-1.

Both maternal and embryonic exposure to mild hypoxia influence embryonic development of the intertidal gastropod Littorina littorea

There is growing evidence that maternal exposure to environmental stressors can alter offspring phenotype and increase fitness. Here, we investigate the relative and combined effects of maternal and developmental exposure to mild hypoxia (65 and 74% air saturation, respectively) on the growth and development of embryos of the marine gastropod Littorina littorea Differences in embryo morphological traits were driven by the developmental environment, whereas the maternal environment and interactive effects of maternal and developmental environment were the main driver of differences in the timing of developmental events. While developmental exposure to mild hypoxia significantly increased the area of an important respiratory organ, the velum, it significantly delayed hatching of veliger larvae and reduced their size at hatching and overall survival. Maternal exposure had a significant effect on these traits, and interacted with developmental exposure to influence the time of appearance of morphological characters, suggesting that both are important in affecting developmental trajectories. A comparison between embryos that successfully hatched and those that died in mild hypoxia revealed that survivors exhibited hypertrophy in the velum and associated pre-oral cilia, suggesting that these traits are linked with survival in low-oxygen environments. We conclude that both maternal and developmental environments shape offspring phenotype in a species with a complex developmental life history, and that plasticity in embryo morphology arising from exposure to even small reductions in oxygen tensions affects the hatching success of these embryos.

A high-throughput and open-source platform for embryo phenomics

Phenomics has the potential to facilitate significant advances in biology but requires the development of high-throughput technologies capable of generating and analysing high-dimensional data. There are significant challenges associated with building such technologies, not least those required for investigating dynamic processes such as embryonic development, during which high rates of temporal, spatial, and functional change are inherently difficult to capture. Here, we present EmbryoPhenomics, an accessible high-throughput platform for phenomics in aquatic embryos comprising an Open-source Video Microscope (OpenVIM) that produces high-resolution videos of multiple embryos under tightly controlled environmental conditions. These videos are then analysed by the Python package Embryo Computer Vision (EmbryoCV), which extracts phenomic data for morphological, physiological, behavioural, and proxy traits during the process of embryonic development. We demonstrate the broad-scale applicability of EmbryoPhenomics in a series of experiments assessing chronic, acute, and multistressor responses to environmental change (temperature and salinity) in >30 million images of >600 embryos of two species with markedly different patterns of development—the pond snail Radix balthica and the marine amphipod Orchestia gammarellus. The challenge of phenomics is significant but so too are the rewards, and it is particularly relevant to the urgent task of assessing complex organismal responses to current rates of environmental change. EmbryoPhenomics can acquire and process data capturing functional, temporal, and spatial responses in the earliest, most dynamic life stages and is potentially game changing for those interested in studying development and phenomics more widely.

EmbryoPhenomics is an open-source technology platform for high-throughput phenome screening of aquatic embryos. This paper demonstrates its application in experiments assessing the sensitivity of aquatic embryos to environmental stress, consisting of more than 600 embryos and more than 30 million images.

Phenomics is the collection of high-dimensional phenotypic data on an organism-wide scale, and it requires high-throughput technologies. However, a lack of technologies for efficiently visualising and measuring whole-organism responses to different environments represents a serious challenge for biologists. This challenge is most apparent when studying complex responses, such as those occurring during the dynamic period of embryonic development, when the phenotype changes markedly through time. Here, we present EmbryoPhenomics (www.embryophenomics.org), a new open-source technological platform comprising high-throughput bioimaging hardware that produces high-resolution video of multiple, developing embryos maintained under controlled environmental conditions and software for automatically quantifying embryo responses from these videos. We demonstrate the broad applicability of EmbryoPhenomics using four experiments assessing responses to global change (elevated temperature and salinity) in which we generate data for more than 600 embryos produced from video comprising more than 30 million images. EmbryoPhenomics was used to capture functional, temporal, and spatial change in morphological, physiological, and behavioural responses in the earliest, most dynamic life stages and addresses a serious bottleneck in biology. Such capabilities are urgently required, particularly within the context of assessing the response of embryos to the current unprecedented rates of global environmental change.

Transcriptomic responses to predator kairomones in embryos of the aquatic snail Radix balthica

The ability of organisms to respond to predation threat by exhibiting induced defenses is well documented, but studies on the potential mechanistic basis for such responses are scarce. Here, we examine the transcriptomic response to predator kairomones of two functionally distinct developmental stages in embryos of the aquatic snail Radix balthica: E8-the stage at which a range-finding trial indicated that kairomone-induced accelerated growth and development first occurred; and E9-the stage at which embryos switched from ciliary- to crawling-driven locomotion. We tested whether expression profiles were influenced by kairomones and whether this influence varied between stages. We also identified potential candidate genes for investigating mechanisms underpinning induced responses. There were 6,741 differentially expressed transcripts between developmental stages, compared to just five in response to predator kairomones. However, on examination of functional enrichment in the transcripts responding to predator kairomones and adopting a less stringent significance threshold, 206 transcripts were identified relating to muscle function, growth, and development, with this response being greater at the later E9 stage. Furthermore, these transcripts included putative annotations for genes identified as responding to predator kairomones in other taxa, including C1q, lectin, and actin domains. Globally, transcript expression appeared reduced in response to predator kairomones and we hypothesize that this might be a result of metabolic suppression, as has been reported in other taxa in response to predation threat.

Thermal strategies vary with life history stage

With both global surface temperatures and the incidence and intensity of extreme temperature events projected to increase, the assessment of species' sensitivity to chronic and acute changes in temperature has become crucial. Sensitivity predictions are based predominantly on adult responses, despite the fact that early life stages may be more vulnerable to thermal challenge. Here, we compared the sensitivity of different life history stages of the intertidal gastropod Littorina obtusata using thermal death time curves, which incorporate the intensity and duration of heat stress, and used these to calculate upper critical thermal limits (CTmax) and sensitivity to temperature change (z). Early (larval) life stages had both a lower CTmax and a lower z than adults, suggesting they are less good at withstanding short-term extreme thermal challenges but better able to survive moderate temperatures in the long term. This result supports the predicted trade-off between acute and chronic tolerance to thermal stress, and is consistent with the different thermal challenges that these stages encounter in the intertidal zone. We conclude that different life history stages employ different thermal strategies that may be adaptive. Our findings caution against the use of predictions of the impact of global warming that are based on only adult responses and, hence, which may underestimate vulnerability.

An annotated draft genome for Radix auricularia (Gastropoda, Mollusca)

Molluscs are the second most species-rich phylum in the animal kingdom, yet only 11 genomes of this group have been published so far. Here, we present the draft genome sequence of the pulmonate freshwater snail Radix auricularia . Six whole genome shotgun libraries with different layouts were sequenced. The resulting assembly comprises 4,823 scaffolds with a cumulative length of 910 Mb and an overall read coverage of 72×. The assembly contains 94.6% of a metazoan core gene collection, indicating an almost complete coverage of the coding fraction. The discrepancy of ∼690 Mb compared with the estimated genome size of R. auricularia (1.6 Gb) results from a high repeat content of 70% mainly comprising DNA transposons. The annotation of 17,338 protein coding genes was supported by the use of publicly available transcriptome data. This draft will serve as starting point for further genomic and population genetic research in this scientifically important phylum.

Regional adaptation defines sensitivity to future ocean acidification

Physiological responses to temperature are known to be a major determinant of species distributions and can dictate the sensitivity of populations to global warming. In contrast, little is known about how other major global change drivers, such as ocean acidification (OA), will shape species distributions in the future. Here, by integrating population genetics with experimental data for growth and mineralization, physiology and metabolomics, we demonstrate that the sensitivity of populations of the gastropod Littorina littorea to future OA is shaped by regional adaptation. Individuals from populations towards the edges of the natural latitudinal range in the Northeast Atlantic exhibit greater shell dissolution and the inability to upregulate their metabolism when exposed to low pH, thus appearing most sensitive to low seawater pH. Our results suggest that future levels of OA could mediate temperature-driven shifts in species distributions, thereby influencing future biogeography and the functioning of marine ecosystems.

Global warming is expected to lead to shifts in species' geographic ranges to track preferred temperatures. Here, the authors show that populations of the common periwinkle vary in their sensitivity to ocean acidification, another major global change driver, which could further restrict range shifts caused by warming.

Heterokairy: a significant form of developmental plasticity?

There is a current surge of research interest in the potential role of developmental plasticity in adaptation and evolution. Here we make a case that some of this research effort should explore the adaptive significance of heterokairy, a specific type of plasticity that describes environmentally driven, altered timing of development within a species. This emphasis seems warranted given the pervasive occurrence of heterochrony, altered developmental timing between species, in evolution. We briefly review studies investigating heterochrony within an adaptive context across animal taxa, including examples that explore links between heterokairy and heterochrony. We then outline how sequence heterokairy could be included within the research agenda for developmental plasticity. We suggest that the study of heterokairy may be particularly pertinent in (i) determining the importance of non-adaptive plasticity, and (ii) embedding concepts from comparative embryology such as developmental modularity and disassociation within a developmental plasticity framework.

Home advantage? Decomposition across the freshwater-estuarine transition zone varies with litter origin and local salinity

Expected increases in the frequency and intensity of storm surges and river flooding may greatly affect the relative salinity of estuarine environments over the coming decades. In this experiment we used detritus from three contrasting environments (marine Fucus vesiculosus; estuarine Spartina anglica; terrestrial Quercus robur) to test the prediction that the decomposition of the different types of litter would be highest in the environment with which they are associated. Patterns of decomposition broadly fitted our prediction: Quercus detritus decomposed more rapidly in freshwater compared with saline conditions while Fucus showed the opposite trend; Spartina showed an intermediate response. Variation in macro-invertebrate assemblages was detected along the salinity gradient but with different patterns between estuaries, suggesting that breakdown rates may be linked in part to local invertebrate assemblages. Nonetheless, our results suggest that perturbation of salinity gradients through climate change could affect the process of litter decomposition and thus alter nutrient cycling in estuarine transition zones. Understanding the vulnerability of estuaries to changes in local abiotic conditions is important given the need to better integrate coastal proceses into a wider management framework at a time when coastlines are increasingly threatened by human activities.

Challenges and opportunities in developmental integrative physiology

This review explores challenges and opportunities in developmental physiology outlined by a symposium at the 2014 American Physiological Society Intersociety Meeting: Comparative Approaches to Grand Challenges in Physiology. Across animal taxa, adverse embryonic/fetal environmental conditions can alter morphological and physiological phenotypes in juveniles or adults, and capacities for developmental plasticity are common phenomena. Human neonates with body sizes at the extremes of perinatal growth are at an increased risk of adult disease, particularly hypertension and cardiovascular disease. There are many rewarding areas of current and future research in comparative developmental physiology. We present key mechanisms, models, and experimental designs that can be used across taxa to investigate patterns in, and implications of, the development of animal phenotypes. Intraspecific variation in the timing of developmental events can be increased through developmental plasticity (heterokairy), and could provide the raw material for selection to produce heterochrony--an evolutionary change in the timing of developmental events. Epigenetics and critical windows research recognizes that in ovo or fetal development represent a vulnerable period in the life history of an animal, when the developing organism may be unable to actively mitigate environmental perturbations. 'Critical windows' are periods of susceptibility or vulnerability to environmental or maternal challenges, periods when recovery from challenge is possible, and periods when the phenotype or epigenome has been altered. Developmental plasticity may allow survival in an altered environment, but it also has possible long-term consequences for the animal. "Catch-up growth" in humans after the critical perinatal window has closed elicits adult obesity and exacerbates a programmed hypertensive phenotype (one of many examples of "fetal programing"). Grand challenges for developmental physiology include integrating variation in developmental timing within and across generations, applying multiple stressor dosages and stressor exposure at different developmental timepoints, assessment of epigenetic and parental influences, developing new animal models and techniques, and assessing and implementing these designs and models in human health and development.

Combining motion analysis and microfluidics--a novel approach for detecting whole-animal responses to test substances

Small, early life stages, such as zebrafish embryos are increasingly used to assess the biological effects of chemical compounds in vivo. However, behavioural screens of such organisms are challenging in terms of both data collection (culture techniques, drug delivery and imaging) and data evaluation (very large data sets), restricting the use of high throughput systems compared to in vitro assays. Here, we combine the use of a microfluidic flow-through culture system, or BioWell plate, with a novel motion analysis technique, (sparse optic flow - SOF) followed by spectral analysis (discrete Fourier transformation - DFT), as a first step towards automating data extraction and analysis for such screenings. Replicate zebrafish embryos housed in a BioWell plate within a custom-built imaging system were subject to a chemical exposure (1.5% ethanol). Embryo movement was videoed before (30 min), during (60 min) and after (60 min) exposure and SOF was then used to extract data on movement (angles of rotation and angular changes to the centre of mass of embryos). DFT was subsequently used to quantify the movement patterns exhibited during these periods and Multidimensional Scaling and ANOSIM were used to test for differences. Motion analysis revealed that zebrafish had significantly altered movements during both the second half of the alcohol exposure period and also the second half of the recovery period compared to their pre-treatment movements. Manual quantification of tail flicking revealed the same differences between exposure-periods as detected using the automated approach. However, the automated approach also incorporates other movements visible in the organism such as blood flow and heart beat, and has greater power to discern environmentally-driven changes in the behaviour and physiology of organisms. We suggest that combining these technologies could provide a highly efficient, high throughput assay, for assessing whole embryo responses to various drugs and chemicals.

Parent--offspring similarity in the timing of developmental events: an origin of heterochrony?

Understanding the link between ontogeny (development) and phylogeny (evolution) remains a key aim of biology. Heterochrony, the altered timing of developmental events between ancestors and descendants, could be such a link although the processes responsible for producing heterochrony, widely viewed as an interspecific phenomenon, are still unclear. However, intraspecific variation in developmental event timing, if heritable, could provide the raw material from which heterochronies originate. To date, however, heritable developmental event timing has not been demonstrated, although recent work did suggest a genetic basis for intraspecific differences in event timing in the embryonic development of the pond snail, Radix balthica. Consequently, here we used high-resolution (temporal and spatial) imaging of the entire embryonic development of R. balthica to perform a parent-offspring comparison of the timing of twelve, physiological and morphological developmental events. Between-parent differences in the timing of all events were good predictors of such timing differences between their offspring, and heritability was demonstrated for two of these events (foot attachment and crawling). Such heritable intraspecific variation in developmental event timing could be the raw material for speciation events, providing a fundamental link between ontogeny and phylogeny, via heterochrony.

A genetic basis for intraspecific differences in developmental timing?

Heterochrony, altered developmental timing between ancestors and their descendents, has been proposed as a pervasive evolutionary feature and recent analytical approaches have confirmed its existence as an evolutionary pattern. Yet, the mechanistic basis for heterochrony remains unclear and, in particular, whether intraspecific variation in the timing of developmental events generates, or has the potential to generate, future between-species differences. Here we make a key step in linking heterochrony at the inter- and intraspecific level by reporting an association between interindividual variation in both the absolute and relative timing (position within the sequence of developmental events) of key embryonic developmental events and genetic distance for the pond snail, Radix balthica. We report significant differences in the genetic distance of individuals exhibiting different levels of dissimilarity in their absolute and relative timing of developmental events such as spinning activity, eyespot formation, heart ontogeny, and hatching. This relationship between genetic and developmental dissimilarity is consistent with there being a genetic basis for variation in developmental timing and so suggests that intraspecific heterochrony could provide the raw material for natural selection to produce speciation.

Development of cardiovascular function in the marine gastropod Littorina obtusata (Linnaeus)

The molluscan cardiovascular system typically incorporates a transient extracardiac structure, the larval heart, early in development, but the functional importance of this structure is unclear. We documented the ontogeny and regulatory ability of the larval heart in relation to two other circulatory structures, the true heart and the velum, in the intertidal gastropod Littorina obtusata. There was a mismatch between the appearance of the larval heart and the velum. Velar lobes appeared early in development (day 4), but the larval heart did not begin beating until day 13. The beating of the larval heart reached a maximum on day 17 and then decreased until the structure itself disappeared (day 24). The true heart began to beat on day 17. Its rate of beating increased as that of the larval heart decreased, possibly suggesting a gradual shift from a larval heart-driven to a true heart-driven circulation. The true heart was not sensitive to acutely declining PO2 shortly after it began to beat, but increased in activity in response to acutely declining PO2 by day 21. Larval heart responses were similar to those of the true heart, with early insensitivity to declining PO2 (day 13) followed by a response by day 15. Increased velum-driven rotational activity under acutely declining PO2 was greatest in early developmental stages. Together, these findings point to cardiovascular function in L. obtusata larvae being the result of a complex interaction between velum, larval and true heart activities, with the functions of the three structures coinciding but their relative importance changing throughout larval development.

Exposure to elevated temperature and Pco(2) reduces respiration rate and energy status in the periwinkle Littorina littorea

In the future, marine organisms will face the challenge of coping with multiple environmental changes associated with increased levels of atmospheric Pco(2), such as ocean warming and acidification. To predict how organisms may or may not meet these challenges, an in-depth understanding of the physiological and biochemical mechanisms underpinning organismal responses to climate change is needed. Here, we investigate the effects of elevated Pco(2) and temperature on the whole-organism and cellular physiology of the periwinkle Littorina littorea. Metabolic rates (measured as respiration rates), adenylate energy nucleotide concentrations and indexes, and end-product metabolite concentrations were measured. Compared with values for control conditions, snails decreased their respiration rate by 31% in response to elevated Pco(2) and by 15% in response to a combination of increased Pco(2) and temperature. Decreased respiration rates were associated with metabolic reduction and an increase in end-product metabolites in acidified treatments, indicating an increased reliance on anaerobic metabolism. There was also an interactive effect of elevated Pco(2) and temperature on total adenylate nucleotides, which was apparently compensated for by the maintenance of adenylate energy charge via AMP deaminase activity. Our findings suggest that marine intertidal organisms are likely to exhibit complex physiological responses to future environmental drivers, with likely negative effects on growth, population dynamics, and, ultimately, ecosystem processes.

Predator cues alter the timing of developmental events in gastropod embryos

Heterochrony, differences in the timing of developmental events between descendent species and their ancestors, is a pervasive evolutionary pattern. However, the origins of such timing changes are still not resolved. Here we show, using sequence analysis, that exposure to predator cues altered the timing of onset of several developmental events in embryos of two closely related gastropod species: Radix balthica and Radix auricularia. These timing alterations were limited to certain events and were species-specific. Compared with controls, over half (62%) of exposed R. auricularia embryos had a later onset of body flexing and an earlier occurrence of the eyes and the heart; in R. balthica, 67 per cent of exposed embryos showed a later occurrence of mantle muscle flexing and an earlier attachment to, and crawling on, the egg capsule wall. The resultant developmental sequences in treated embryos converged, and were more similar to one another than were the sequences of the controls for both species. We conclude that biotic agents can elicit altered event timing in developing gastropod embryos. These changes were species-specific, but did not occur in all individuals. Such developmental plasticity in the timing of developmental events could be an important step in generating interspecific heterochrony.

Alarm substance from adult zebrafish alters early embryonic development in offspring

Alarm substances elicit behavioural responses in a wide range of animals but effects on early embryonic development are virtually unknown. Here we investigated whether skin injury-induced alarm substances caused physiological responses in embryos produced by two Danio species (Danio rerio and Danio albolineatus). Both species showed more rapid physiological development in the presence of alarm substance, although there were subtle differences between them: D. rerio had advanced muscle contraction and heart function, whereas D. albolineatus had advanced heart function only. Hence, alarm cues from injured or dying fish may be of benefit to their offspring, inducing physiological responses and potentially increasing their inclusive fitness.

Experimental determination of a Viviparus contectus thermometry equation

Experimental measurements of the (18)O/(16)O isotope fractionation between the biogenic aragonite of Viviparus contectus (Gastropoda) and its host freshwater were undertaken to generate a species-specific thermometry equation. The temperature dependence of the fractionation factor and the relationship between Deltadelta(18)O (delta(18)O(carb.) - delta(18)O(water)) and temperature were calculated from specimens maintained under laboratory and field (collection and cage) conditions. The field specimens were grown (Somerset, UK) between August 2007 and August 2008, with water samples and temperature measurements taken monthly. Specimens grown in the laboratory experiment were maintained under constant temperatures (15 degrees C, 20 degrees C and 25 degrees C) with water samples collected weekly. Application of a linear regression to the datasets indicated that the gradients of all three experiments were within experimental error of each other (+/-2 times the standard error); therefore, a combined (laboratory and field data) correlation could be applied. The relationship between Deltadelta(18)O (delta(18)O(carb.) - delta(18)O(water)) and temperature (T) for this combined dataset is given by: T = - 7.43( + 0.87, - 1.13)*Deltadelta18O + 22.89(+/- 2.09) (T is in degrees C, delta(18)O(carb.) is with respect to Vienna Pee Dee Belemnite (VPDB) and delta(18)O(water) is with respect to Vienna Standard Mean Ocean Water (VSMOW). Quoted errors are 2 times standard error).Comparisons made with existing aragonitic thermometry equations reveal that the linear regression for the combined Viviparus contectus equation is within 2 times the standard error of previously reported aragonitic thermometry equations. This suggests there are no species-specific vital effects for Viviparus contectus. Seasonal delta(18)O(carb.) profiles from specimens retrieved from the field cage experiment indicate that during shell secretion the delta(18)O(carb.) of the shell carbonate is not influenced by size, sex or whether females contained eggs or juveniles.

The effect of CO2 acidified sea water and reduced salinity on aspects of the embryonic development of the amphipod Echinogammarus marinus (Leach)

We investigated the effect of CO(2) acidified sea water (S=35, 22 and 10(PSU)) on embryonic development of the intertidal amphipod Echinogammarus marinus (Leach). Low pH, but not low salinity (22(PSU)), resulted in a more protracted embryonic development in situ although the effect was only evident at low salinity. However reduced salinity, not pH, exerted a strong significant effect, on numbers and calcium content of hatchlings. Females exposed to low salinity (10(PSU)) did not carry eggs through to hatching. There was no significant difference in the number of viable hatchlings between females cultured in 22 and 35(PSU) but the exoskeleton of the juveniles at 22(PSU) contained significantly less calcium. Ocean acidification may affect aspects of E. marinus development but exposure to realistic low salinities appear, in the short term, to be more important in impacting development than exposure to CO(2) acidified sea water at levels predicted for 300 years time.

Embryonic rotational behaviour in the pond snail Lymnaea stagnalis: influences of environmental oxygen and development stage

Responses of freshwater organisms to environmental oxygen tensions (PO(2)) have focused on adult (i.e. late developmental) stages, yet responses of embryonic stages to changes in environmental PO(2) must also have implications for organismal biology. Here we assess how the rotational behaviour of the freshwater snail Lymnaea stagnalis changes during development in response to conditions of hypoxia and hyperoxia. As rotation rate is linked to gas mixing in the fluid surrounding the embryo, we predicted that it would increase under hypoxic conditions but decrease under hyperoxia. Contrary to predictions, early, veliger stage embryos showed no change in their rotation rate under hyperoxia, and later, hippo stage embryos showed only a marginally significant increase in rotation under these conditions. Predictions for hypoxia were broadly supported, however, with both veliger and hippo stages showing a marked hypoxia-related increase in their rotation rates. There were also subtle differences between developmental stages, with hippos responding at PO(2)s (50% air saturation) greater than those required to elicit a similar response in veligers (20% air saturation). Differences between developmental stages also occurred on return to normoxic conditions following hypoxia: rotation in veligers returned to pre-exposure levels, whereas there was a virtual cessation in embryos at the hippo stage, likely the result of overstimulation of oxygen sensors driving ciliary movement in later, more developed embryos. Together, these findings suggest that the spinning activity of L. stagnalis embryos varies depending on environmental PO(2)s and developmental stage, increasing during hypoxia to mix capsular contents and maintain a diffusive gradient for oxygen entry into the capsule from the external environment ("stir-bar" theory of embryonic rotational behaviour).

Comparing the strength of behavioural plasticity and consistency across situations: animal personalities in the hermit crab Pagurus bernhardus

Many phenotypic traits show plasticity but behaviour is often considered the 'most plastic' aspect of phenotype as it is likely to show the quickest response to temporal changes in conditions or 'situation'. However, it has also been noted that constraints on sensory acuity, cognitive structure and physiological capacities place limits on behavioural plasticity. Such limits to plasticity may generate consistent differences in behaviour between individuals from the same population. It has recently been suggested that these consistent differences in individual behaviour may be adaptive and the term 'animal personalities' has been used to describe them. In many cases, however, a degree of both behavioural plasticity and relative consistency is probable. To understand the possible functions of animal personalities, it is necessary to determine the relative strength of each tendency and this may be achieved by comparison of statistical effect sizes for tests of difference and concordance. Here, we describe a new statistical framework for making such comparisons and investigate cross-situational plasticity and consistency in the duration of startle responses in the European hermit crab Pagurus bernhardus, in the field and the laboratory. The effect sizes of tests for behavioural consistency were greater than for tests of behavioural plasticity, indicating for the first time the presence of animal personalities in a crustacean model.

Phylogenetic relatedness and ecological interactions determine antipredator behavior

Interspecific recognition of alarm cues among guild members through "eavesdropping" may allow prey to fine-tune antipredator responses. This process may be linked to taxonomic relatedness but might also be influenced by local adaptation to recognize alarm cues from sympatric species. We tested this hypothesis using antipredator responses of a freshwater gastropod Lymnaea stagnalis (L.) to alarm cues from damaged conspecific and 10 heterospecific gastropod species. As predicted, the magnitude of antipredator response decreased significantly with increasing phylogenetic distance, but increased when species were naturally sympatric (defined as species cohabiting in the same water body) with the source population of L. stagnalis. The responses to sympatric species were higher overall, and the relationship between genetic distance and alarm cue response was stronger when tested with sympatric species. This is the first study to demonstrate that population sympatry influences innate antipredator responses to alarm cues from intraguild members and suggests that responses based on phylogenetic relationships can be modified through local adaptation. Such adaptation to heterospecific alarm cues suggests that species could be at a disadvantage when they encounter novel intraguild members resulting from species invasion or range expansion due to a reduction in the presence of reliable information about predation risk.

Plasticity in the timing of physiological development: Physiological heterokairy — What is it, how frequent is it, and does it matter?

The study of developmental sequences of physiological traits could be an important way of placing comparative developmental physiology (CDP) within the research agenda being forged by work on developmental plasticity. Here we focus on the concept of heterokairy defined by Spicer & Burggren in 2003 as changes in the timing of physiological development in an individual. The role of this concept in the future of the CDP is discussed. First we provide an historical perspective of the ideas that have led to the investigation of sequences in CDP. This is followed by a re-examination and clarification of the definition of physiological heterokairy before empirical case studies that (explicitly or implicitly) demonstrate physiological heterokairy are reviewed. We suggest that physiological heterokairy can be demonstrated through a wide range of invertebrate and vertebrate examples. However, care must be taken when inferring that heterokairy as a pattern is always the result of heterokairic processes as there is evidence that physiological heterokairy could result from the altered timing of both homologous or analogous physiological mechanisms. We conclude by discussing the potential link between heterokairy and heterochrony and suggest that the investigation of this link should be a major goal for workers in both CDP and developmental plasticity.

An integrative approach identifies developmental sequence heterochronies in freshwater basommatophoran snails

Adopting an integrative approach to the study of sequence heterochrony, we compared the timing of developmental events encompassing a mixture of developmental stages and functional traits in the embryos of 12 species of basommatophoran snails in an explicit phylogenetic framework. PARSIMOV analysis demonstrated clear functional heterochronies associated both with basal branches within the phylogeny and with terminal speciation events. A consensus of changes inferred under both accelerated transformation and delayed transformation optimizations identified four heterochronies where the direction of movement was known plus six twin heterochronies where the relative movements of the two events could not be assigned. On average, 0.5 and 0.58 events were inferred to have changed their position in the developmental sequence on internal and terminal branches of the phylogeny, respectively; these values are comparable with frequencies of sequence heterochrony reported in mammals. Directional heterochronies such as the early occurrence of body flexing in relation to the ontogeny of the eye spots, heart beat, and free swimming events occurred convergently and/or at different levels (i.e., familial, generic, and species) within the phylogeny. Such a functional approach to the study of developmental sequences has highlighted the possibility that heterochrony may have played a prominent role in the evolution of this group of invertebrates.

Do developmental mode and dispersal shape abundance-occupancy relationships in marine macroinvertebrates?

1. Dispersal is a crucial process in maintaining population structures in many organisms, and is hypothesized as a process underlying the interspecific relationship between abundance and distribution. Here we examined whether there was a link between the dispersal and developmental modes of marine macroinvertebrates and the slopes and elevations of interspecific abundance-occupancy relationships. We predicted that if within-site retention of larvae ranks in the order brooders > lecithotrophs > planktotrophs, for any given level of mean abundance, occupancy should increase in the order brooders < lecithotrophs < planktotrophs. We also predicted that propensity to form metapopulations should be greater for planktonic dispersers (i.e. lecithotrophs and planktotrophs combined) than for non-planktonic (i.e. brooders), resulting in steeper abundance-occupancy relationships for the former. 2. Predictions were tested using a data set for 362 subtidal marine macroinvertebrates occurring across 446 1-km(2) grid squares around the British Isles; analyses were performed on the data set as a whole and for separate phyla. 3. The total data set had a Z-transformed effect size of 0.79, within the confidence intervals described by Blackburn et al. (2006; Journal of Animal Ecology, 75, 1426-1439), and was consistently present with relatively homogeneous effect size in separate analyses of polychaetes, crustaceans, molluscs and echinoderms. 4. In all cases, planktonic dispersing organisms showed an abundance-occupancy relationship with greater elevation than that for non-planktonic organisms; in polychaetes the elevation of slopes was in the rank order planktotrophs > lecithotrophs > brooders. No differences between the slopes of the abundance-occupancy relationship were apparent for different dispersal modes either within, or across phyla. 5. We conclude that dispersal capacity may play an important part in determining the elevation of the abundance-occupancy relationship, the corollary of low dispersal in the marine realm being greater local retention of larvae and greater local population abundance at low extents of geographical distribution.

A microplate freshwater copepod bioassay for evaluating acute and chronic effects of chemicals

Toxicity test protocols for establishing the effect of zinc on the freshwater copepod Bryocamptus zschokkei are reported. In the absence of food, larval life-history stages were more sensitive than adult females to acute zinc exposure with a 96-h median lethal concentration of 0.62 mg Zn/L (0.52-0.73 mg Zn/L, 95% confidence intervals) for copepodids. The acute toxicity of zinc to adult females was also significantly reduced in the presence of food (a leaf disc). The main toxic effect of zinc in a life-cycle test was a reduction in the number of offspring per brood at 0.48 mg Zn/L. As this endpoint corresponded with zinc concentrations causing reduced juvenile survival, prolonged embryonic development times, and mortality during hatching, it appears that zinc had a direct toxic effect on the current brood rather than an indirect effect on egg production via maternal reallocation of resources. The lowest-observed-effect concentration for reproduction of measured zinc concentrations of 0.48 mg Zn/L corresponded with environmental concentrations of zinc causing reduced population densities of this copepod in the field. It is proposed that B. zschokkei is an ecologically relevant test species for evaluating the potential effects of contaminants on freshwater invertebrates and for giving insight into the mode of action of a test chemical.

Environmental calcium modifies induced defences in snails

Inducible defences are adaptive phenotypes that arise in response to predation threats. Such plasticity incurs costs to individuals, but there has been little interest in how such induced traits in animals may be constrained by environmental factors. Here, we demonstrate that calcium availability interacts with predation cues to modify snail shell growth and form. Small snails increased their growth and were heavier when exposed to fish chemical cues, but this response was calcium limited. There was also an interactive effect of fish cues and calcium on the shell growth of larger snails, but shell strength and aperture narrowness were affected by calcium alone. For small snails, behavioural avoidance was greatest for snails exhibiting least morphological plasticity, suggesting a trade-off. There was no trade-off of somatic growth with plasticity. We suggest that the expression of defensive traits in molluscs can be constrained by calcium availability, which has implications for molluscan ecology and evolution.

A copepod life-cycle test and growth model for interpreting the effects of lindane

A full life-cycle test was performed to measure the effects of lindane (3.2-3,200 microg l(-1)) on the survival, development and reproduction of the freshwater copepod Bryocamptus zschokkei. This copepod survived at relatively high concentrations of lindane compared with other freshwater crustaceans with a 10 day LC50 of 241 microg l(-1) (95% CL of 141-440). 'Equiproportional development', which assumes that each moult stage represents a specific proportion of the total development time, and is not affected by processes that influence metabolism such as temperature and food quality, was used to determine the mode of action of lindane on development in B. zschokkei. Development to adult was significantly longer at 100 microg l(-1) lindane compared with the controls, however, development remained equiproportional regardless of lindane exposure. Increased development times, therefore, are not due to a direct effect of lindane on the moulting process but are due probably to reduced food intake or increased metabolism through the stress imposed by toxicant exposure. Although the survival data suggest that B. zschokkei is relatively tolerant of lindane exposure, reproduction was affected at low lindane concentrations. At 32 microg l(-1) lindane, significantly fewer eggs and viable offspring were produced per female compared with the solvent control. At very low lindane concentrations (3.2 and 10 microg l(-1)), there was a significant increase in the numbers of offspring produced per female compared with the controls and this is interpreted as a hormesis effect. In conclusion, a full life-cycle test demonstrated B. zschokkei is relatively sensitive to lindane compared with other freshwater crustaceans. Incorporating a copepod growth model (equiproportional development) into the life-cycle test design, provided information on the dominant mode of action of the toxicant.

Inter- and intraspecific trait compensation of defence mechanisms in freshwater snails

Trait compensation occurs when mechanically independent adaptations are negatively correlated. Here, we report the first study to demonstrate trait compensation in predator-defence adaptations across several species. Freshwater pulmonate snails exposed experimentally to predation chemical cues from fishes and crushed conspecifics showed clear interspecific differences in their behavioural avoidance responses, which were negatively correlated with shell crush resistance. The type of avoidance response varied between species: thin-shelled species (Lymnaea stagnalis and Physa fontinalis) moved to the water-line or out of the water, while those with thick shells moved under cover or showed a mixed response. There were also intraspecific size-linked differences, with an ontogenetic increase in shell strength accompanied by a decrease in behavioural avoidance. Such trait compensation in response to predation has important implications for interspecific interactions and food-web dynamics.

Population structure and dispersal in the Canary Island caddisfly Mesophylax aspersus (Trichoptera, Limnephilidae)

Population genetic structure of the circum-Mediterranean caddisfly Mesophylax aspersus (Trichoptera, Limnephilidae) on the Canary Islands was investigated by studying allozyme variation at nine putative loci in five populations. Genetic variability, population structure and gene flow were compared with data in the literature for continental taxa to assess the effect of isolation of island populations on the genetic structure. Larvae were collected from streams on the islands of Tenerife (one population), La Gomera (two populations in the same catchment) and La Palma (two populations in different catchments). Genetic variability within populations was high relative to that recorded previously for continental Trichoptera, e.g. mean heterozygosity was 0.119--0.336 (0.035--0.15 in continental taxa). Highly significant population structuring was observed (mean F(ST)=0.250), and there was significant within-population structuring (mean F(IS)=0.098). Populations from the same catchment or island were no more similar than populations from different islands, which suggests that occasional long-distance dispersal, both between and within islands, is the predominant influence on the population structure. This dispersal ability has contributed to the colonization of most permanent streams on the Canary Islands by M. aspersus.

The relationship between trace metal contamination and stream meiofauna

The effect of chemical contaminants on freshwater meiofaunal communities is poorly understood and meiofauna rarely form part of environmental impact assessments in fresh waters. The community composition and diversity of meiofauna in streams of southwest England (Cornwall) representing a gradient in trace metal contamination were investigated. Multivariate and univariate statistical techniques were used to correlate community composition, diversity and environmental variables. Meiofaunal communities were very different at sites with high compared with low metal concentrations. Copper, either singly, or in combination with other environmental variables (aluminium, zinc or dissolved organic carbon), was the most important correlate with community composition. Not all meiofaunal species were adversely affected by metals, however, and some taxa, most notably certain cyclopoid copepods, were abundant at high metal concentrations. Moreover, sites with high metal concentrations were not significantly less diverse than sites of low metal concentrations. These data suggest that metal contamination significantly alters the composition of stream meiofaunal assemblages and highlight the importance of including meiofauna in impact studies of fresh waters.