Temperature influences selective mortality during the early life stages of a coral reef fish

Quantifying dispersal variability among nearshore marine populations

Dispersal drives diverse processes from population persistence to community dynamics. However, the amount of temporal variation in dispersal and its consequences for metapopulation dynamics is largely unknown for organisms with environmentally driven dispersal (e.g., many marine larvae, arthropods and plant seeds). Here, we used genetic parentage analysis to detect larval dispersal events in a common coral reef fish, Amphiprion clarkii, along 30 km of coastline consisting of 19 reef patches in Ormoc Bay, Leyte, Philippines. We quantified variation in the dispersal kernel across seven years (2012-2018) and monsoon seasons with 71 parentage assignments from 791 recruits and 1,729 adults. Connectivity patterns differed significantly among years and seasons in the scale and shape but not in the direction of dispersal. This interannual variation in dispersal kernels introduced positive temporal covariance among dispersal routes that theory predicts is likely to reduce stochastic metapopulation growth rates below the growth rates expected from only a single or a time-averaged connectivity estimate. The extent of variation in mean dispersal distance observed here among years is comparable in magnitude to the differences across reef fish species. Considering dispersal variation will be an important avenue for further metapopulation and metacommunity research across diverse taxa.

Reproductive phenology across the lunar cycle: parental decisions, offspring responses, and consequences for reef fish

Most organisms reproduce in a dynamic environment, and life-history theory predicts that this can favor the evolution of strategies that capitalize on good times and avoid bad times. When offspring experience these environmental changes, fitness can depend strongly upon environmental conditions at birth and at later life stages. Consequently, fitness will be influenced by the reproductive decisions of parents (i.e., birth date effects) and developmental decisions (e.g., adaptive plasticity) of their offspring. We explored the consequences of these decisions using a highly iteroparous coral reef fish (the sixbar wrasse, Thalassoma hardwicke) and in a system where both parental and offspring environments vary with the lunar cycle. We tested the hypotheses that (1) reproductive patterns and offspring survival vary across the lunar cycle and (2) offspring exhibit adaptive plasticity in development time. We evaluated temporal variation in egg production from February to June 2017, and corresponding larval developmental histories (inferred from otolith microstructure) of successful settlers and surviving juveniles that were spawned during that same period. We documented lunar-cyclic variation in egg production (most eggs were spawned at the new moon). This pattern was at odds with the distribution of birth dates of settlers and surviving juveniles-most individuals that successfully survived to settlement and older stages were born during the full moon. Consequently, the probability of survival across the larval stage was greatest for offspring born close to the full moon, when egg production was at its lowest. Offspring also exhibited plasticity in developmental duration, adjusting their age at settlement to settle during darker portions of the lunar cycle than expected given their birth date. Offspring born near the new moon tended to be older and larger at settlement, and these traits conveyed a strong fitness advantage (i.e., a carryover effect) through to adulthood. We speculate that these effects (1) are shaped by a dynamic landscape of risk and reward determined by moonlight, which differentially influences adults and offspring, and (2) can explain the evolution of extreme iteroparity in sixbars.

Patterns of selective predation change with ontogeny but not density in a marine fish

Phenotypic variation is prevalent in the early life-history stages of many organisms and provides the basis for selective mortality on size and growth-related traits of older life stages. Densities of organisms can vary widely at important life-history transitions, raising additional questions about the interplay between selection and density-dependent processes. We evaluate density dependence in patterns of selective mortality for a temperate reef fish. Specifically, we exposed pre-settlement and post-settlement stages of the common triplefin (Forsterygion lapillum) to a natural predator and evaluated patterns of selective mortality on early life-history traits as a function of ontogenetic stage and density. We used otoliths to reconstruct the traits of fish that survived versus fish that were consumed (i.e., we recovered otoliths from the guts of predators), and we estimated selection by analysing the relationship between absolute fitness and standardised traits. Absolute fitness was negatively correlated with size and larval growth rate for pre-settlement fish (i.e., larger and faster growing individuals were more likely to be consumed by predators), and this was consistent across the range of densities evaluated. Post-settlement fish experienced no selective mortality. Additionally, absolute fitness was equal across density treatments, suggesting mortality was density-independent. Collectively, these results suggest that patterns of selection change with ontogeny, but may be stable across densities when mortality is density-independent. Shifts in selective mortality for species with distinct life-stages can mask and complicate relationships between traits and fitness, and the importance of such traits may be underappreciated for earlier life stages.

Reconstructing larval growth and habitat use in an amphidromous goby using otolith increments and microchemistry

High-resolution analysis of growth increments, trace element chemistry and oxygen isotope ratios (δ¹⁸ O) in otoliths were combined to assess larval and post-larval habitat use and growth of Awaous stamineus, an amphidromous goby native to Hawai'i. Otolith increment widths indicate that all individuals experience a brief period of rapid growth during early life as larvae and that the duration of this growth anomaly is negatively correlated with larval duration. A protracted high-growth period early in larval life is associated with a lower ratio of Sr:Ca, which may reflect low salinity conditions in nearshore habitats. A distinct shift in δ¹⁸ O (range: 4-5‰) is closely associated with the metamorphic mark in otoliths, indicating that larval metamorphosis occurs promptly upon return to fresh water. Strontium and other trace elements are not as tightly coupled to the metamorphosis mark, but confirm the marine-to-freshwater transition. Integration of microstructural and microchemical approaches reveals that larvae vary substantially in growth rate, possibly in association with habitat differences. Although time and financial costs make it difficult to achieve large sample sizes, present results show that examining even a small number of individuals can lead to novel inferences about early life history in diadromous fishes and illustrates the value of integrating analyses.

Demographic and reproductive plasticity across the depth distribution of a coral reef fish

As humans expand into natural environments, populations of wild organisms may become relegated to marginal habitats at the boundaries of their distributions. In the ocean, mesophotic coral ecosystems (30–150 m) at the depth limit of photosynthetic reefs are hypothesized to act as refuges that are buffered from anthropogenic and natural disturbances, yet the viability and persistence of subpopulations in these peripheral habitats remains poorly understood. To assess the potential for mesophotic reefs to support robust coral reef fish populations, we compared population density and structure, growth, size, and reproductive output of the bicolor damselfish (Stegastes partitus) from shallow (<10 m), deep shelf (20–30 m), and mesophotic reefs (60–70 m) across the Florida Platform. Population densities decreased and size and age distributions shifted toward older and larger individuals in deeper habitats. Otolith-derived ages indicated that S. partitus found on mesophotic reefs reach larger asymptotic sizes and have longer lifespans than fish in shallower habitats. Based on measurements of oocyte area and batch fecundity, mesophotic fish also have higher reproductive investment. These demographic patterns indicate that mesophotic fish populations composed of large, fecund individuals produce high condition larvae and rely on longevity of individuals for population persistence and viability.

Productivity and sea surface temperature are correlated with the pelagic larval duration of damselfishes in the Red Sea

We examined the variation of pelagic larval durations (PLDs) among three damselfishes, Dascyllus aruanus, D. marginatus, and D. trimaculatus, which live under the influence of an environmental gradient in the Red Sea. PLDs were significantly correlated with latitude, sea surface temperature (SST), and primary production (CHLA; chlorophyll a concentrations). We find a consistent decrease in PLDs with increasing SST and primary production (CHLA) towards the southern Red Sea among all species. This trend is likely related to higher food availability and increased metabolic rates in that region. We suggest that food availability is a potentially stronger driver of variation in PLD than temperature, especially in highly oligotrophic regions. Additionally, variations in PLDs were particularly high among specimens of D. marginatus, suggesting a stronger response to local environmental differences for endemic species. We also report the first average PLD for this species over a broad geographic range (19.82 ± 2.92 days).

Climate change and larval transport in the ocean: fractional effects from physical and physiological factors

Changes in larval import, export, and self-seeding will affect the resilience of coral reef ecosystems. Climate change will alter the ocean currents that transport larvae and also increase sea surface temperatures (SST), hastening development, and shortening larval durations. Here, we use transport simulations to estimate future larval connectivity due to: (1) physical transport of larvae from altered circulation alone, and (2) the combined effects of altered currents plus physiological response to warming. Virtual larvae from islands throughout Micronesia were moved according to present-day and future ocean circulation models. The Hybrid Coordinate Ocean Model (HYCOM) spanning 2004-2012 represented present-day currents. For future currents, we altered HYCOM using analysis from the National Center for Atmospheric Research Community Earth System Model, version 1-Biogeochemistry, Representative Concentration Pathway 8.5 experiment. Based on the NCAR model, regional SST is estimated to rise 2.74 °C which corresponds to a ~17% decline in larval duration for some taxa. This reduction was the basis for a separate set of simulations. Results predict an increase in self-seeding in 100 years such that 62-76% of islands experienced increased self-seeding, there was an average domainwide increase of ~1-3% points in self-seeding, and increases of up to 25% points for several individual islands. When changed currents alone were considered, approximately half (i.e., random) of all island pairs experienced decreased connectivity but when reduced PLD was added as an effect, ~65% of connections were weakened. Orientation of archipelagos relative to currents determined the directional bias in connectivity changes. There was no universal relationship between climate change and connectivity applicable to all taxa and settings. Islands that presently export large numbers of larvae but that also maintain or enhance this role into the future should be the focus of conservation measures that promote long-term resilience of larval supply.

Born small, die young: Intrinsic, size-selective mortality in marine larval fish

Mortality during the early stages is a major cause of the natural variations in the size and recruitment strength of marine fish populations. In this study, the relation between the size-at-hatch and early survival was assessed using laboratory experiments and on field-caught larvae of the European sardine (Sardina pilchardus). Larval size-at-hatch was not related to the egg size but was significantly, positively related to the diameter of the otolith-at-hatch. Otolith diameter-at-hatch was also significantly correlated with survival-at-age in fed and unfed larvae in the laboratory. For sardine larvae collected in the Bay of Biscay during the spring of 2008, otolith radius-at-hatch was also significantly related to viability. Larval mortality has frequently been related to adverse environmental conditions and intrinsic factors affecting feeding ability and vulnerability to predators. Our study offers evidence indicating that a significant portion of fish mortality occurs during the endogenous (yolk) and mixed (yolk /prey) feeding period in the absence of predators, revealing that marine fish with high fecundity, such as small pelagics, can spawn a relatively large amount of eggs resulting in small larvae with no chances to survive. Our findings help to better understand the mass mortalities occurring at early stages of marine fish.

Characteristics of settling coral reef fish are related to recruitment timing and success

Many marine populations exhibit high variability in the recruitment of young into the population. While environmental cycles and oceanography explain some patterns of replenishment, the role of other growth-related processes in influencing settlement and recruitment is less clear. Examination of a 65-mo. time series of recruitment of a common coral reef fish, Stegastes partitus, to the reefs of the upper Florida Keys revealed that during peak recruitment months, settlement stage larvae arriving during dark lunar phases grew faster as larvae and were larger at settlement compared to those settling during the light lunar phases. However, the strength and direction of early trait-mediated selective mortality also varied by settlement lunar phase such that the early life history traits of 2–4 week old recruit survivors that settled across the lunar cycle converged to more similar values. Similarly, within peak settlement periods, early life history traits of settling larvae and selective mortality of recruits varied by the magnitude of the settlement event: larvae settling in larger events had longer PLDs and consequently were larger at settlement than those settling in smaller pulses. Traits also varied by recruitment habitat: recruits surviving in live coral habitat (vs rubble) or areas with higher densities of adult conspecifics were those that were larger at settlement. Reef habitats, especially those with high densities of territorial conspecifics, are more challenging habitats for young fish to occupy and small settlers (due to lower larval growth and/or shorter PLDs) to these habitats have a lower chance of survival than they do in rubble habitats. Settling reef fish are not all equal and the time and location of settlement influences the likelihood that individuals will survive to contribute to the population.

Phenotypic variation and selective mortality as major drivers of recruitment variability in fishes

An individual's phenotype will usually influence its probability of survival. However, when evaluating the dynamics of populations, the role of selective mortality is not always clear. Not all mortality is selective, patterns of selective mortality may vary, and it is often unknown how selective mortality compares or interacts with other sources of mortality. As a result, there is seldom a clear expectation for how changes in the phenotypic composition of populations will translate into differences in average survival. We address these issues by evaluating how selective mortality affects recruitment of fish populations. First, we provide a quantitative review of selective mortality. Our results show that most of the mortality during early life is selective, and that variation in phenotypes can have large effects on survival. Next, we describe an analytical framework that accounts for variation in selection, while also describing the amount of selective mortality experienced by different cohorts recruiting to a single population. This framework is based on reconstructing fitness surfaces from phenotypic selection measurements, and can be employed for either single or multiple traits. Finally, we show how this framework can be integrated with models of density-dependent survival to improve our understanding of recruitment variability and population dynamics.

The relationship between offspring size and fitness: integrating theory and empiricism

How parents divide the energy available for reproduction between size and number of offspring has a profound effect on parental reproductive success. Theory indicates that the relationship between offspring size and offspring fitness is of fundamental importance to the evolution of parental reproductive strategies: this relationship predicts the optimal division of resources between size and number of offspring, it describes the fitness consequences for parents that deviate from optimality, and its shape can predict the most viable type of investment strategy in a given environment (e.g., conservative vs. diversified bet-hedging). Many previous attempts to estimate this relationship and the corresponding value of optimal offspring size have been frustrated by a lack of integration between theory and empiricism. In the present study, we draw from C. Smith and S. Fretwell's classic model to explain how a sound estimate of the offspring size--fitness relationship can be derived with empirical data. We evaluate what measures of fitness can be used to model the offspring size--fitness curve and optimal size, as well as which statistical models should and should not be used to estimate offspring size--fitness relationships. To construct the fitness curve, we recommend that offspring fitness be measured as survival up to the age at which the instantaneous rate of offspring mortality becomes random with respect to initial investment. Parental fitness is then expressed in ecologically meaningful, theoretically defensible, and broadly comparable units: the number of offspring surviving to independence. Although logistic and asymptotic regression have been widely used to estimate offspring size-fitness relationships, the former provides relatively unreliable estimates of optimal size when offspring survival and sample sizes are low, and the latter is unreliable under all conditions. We recommend that the Weibull-1 model be used to estimate this curve because it provides modest improvements in prediction accuracy under experimentally relevant conditions.

Response to ocean acidification in larvae of a large tropical marine fish, Rachycentron canadum

Currently, ocean acidification is occurring at a faster rate than at any time in the last 300 million years, posing an ecological challenge to marine organisms globally. There is a critical need to understand the effects of acidification on the vulnerable larval stages of marine fishes, as there is potential for large ecological and economic impacts on fish populations and the human economies that rely on them. We expand upon the narrow taxonomic scope found in the literature today, which overlooks many life history characteristics of harvested species, by reporting on the larvae of Rachycentron canadum (cobia), a large, highly mobile, pelagic-spawning, widely distributed species with a life history and fishery value contrasting other species studied to date. We raised larval cobia through the first 3 weeks of ontogeny under conditions of predicted future ocean acidification to determine effects on somatic growth, development, otolith formation, swimming ability, and swimming activity. Cobia exhibited resistance to treatment effects on growth, development, swimming ability, and swimming activity at 800 and 2100 μatm pCO2 . However, these scenarios resulted in a significant increase in otolith size (up to 25% larger area) at the lowest pCO2 levels reported to date, as well as the first report of significantly wider daily otolith growth increments. When raised under more extreme scenarios of 3500 and 5400 μatm pCO2 , cobia exhibited significantly reduced size-at-age (up to 25% smaller) and a 2-3 days developmental delay. The robust nature of cobia may be due to the naturally variable environmental conditions this species currently encounters throughout ontogeny in coastal environments, which may lead to an increased acclimatization ability even during long-term exposure to stressors.

Seasonal and spatial heterogeneity of recent sea surface temperature trends in the Caribbean Sea and southeast Gulf of Mexico

Recent changes in ocean temperature have impacted marine ecosystem function globally. Nevertheless, the responses have depended upon the rate of change of temperature and the season when the changes occur, which are spatially variable. A rigorous statistical analysis of sea surface temperature observations over 25 years was used to examine spatial variability in overall and seasonal temperature trends within the wider Caribbean. The basin has experienced high spatial variability in rates of change of temperature. Most of the warming has been due to increases in summer rather than winter temperatures. However, warming was faster in winter in the Loop Current area and the south-eastern Caribbean, where the annual temperature ranges have contracted. Waters off Florida, Cuba and the Bahamas had a tendency towards cooling in winter, increasing the amplitude of annual temperature ranges. These detailed patterns can be used to elucidate ecological responses to climatic change in the region.