An evolutionary explanation of female-biased sexual size dimorphism in North Sea plaice, Pleuronectes platessa L.

Sexual size dimorphism (SSD) is caused by differences in selection pressures and life-history trade-offs faced by males and females. Proximate causes of SSD may involve sex-specific mortality, energy acquisition, and energy expenditure for maintenance, reproductive tissues, and reproductive behavior. Using a quantitative, individual-based, eco-genetic model parameterized for North Sea plaice, we explore the importance of these mechanisms for female-biased SSD, under which males are smaller and reach sexual maturity earlier than females (common among fish, but also arising in arthropods and mammals). We consider two mechanisms potentially serving as ultimate causes: (a) Male investments in male reproductive behavior might evolve to detract energy resources that would otherwise be available for somatic growth, and (b) diminishing returns on male reproductive investments might evolve to reduce energy acquisition. In general, both of these can bring about smaller male body sizes. We report the following findings. First, higher investments in male reproductive behavior alone cannot explain the North Sea plaice SSD. This is because such higher reproductive investments require increased energy acquisition, which would cause a delay in maturation, leading to male-biased SSD contrary to observations. When accounting for the observed differential (lower) male mortality, maturation is postponed even further, leading to even larger males. Second, diminishing returns on male reproductive investments alone can qualitatively account for the North Sea plaice SSD, even though the quantitative match is imperfect. Third, both mechanisms can be reconciled with, and thus provide a mechanistic basis for, the previously advanced Ghiselin-Reiss hypothesis, according to which smaller males will evolve if their reproductive success is dominated by scramble competition for fertilizing females, as males would consequently invest more in reproduction than growth, potentially implying lower survival rates, and thus relaxing male-male competition. Fourth, a good quantitative fit with the North Sea plaice SSD is achieved by combining both mechanisms while accounting for sex-specific costs males incur during their spawning season. Fifth, evolution caused by fishing is likely to have modified the North Sea plaice SSD.

Trawl impacts on the relative status of biotic communities of seabed sedimentary habitats in 24 regions worldwide.

Bottom trawling is widespread globally and impacts seabed habitats. However, risks from trawling remain unquantified at large scales in most regions. We address these issues by synthesizing evidence on the impacts of different trawl-gear types, seabed recovery rates, and spatial distributions of trawling intensity in a quantitative indicator of biotic status (relative amount of pretrawling biota) for sedimentary habitats, where most bottom-trawling occurs, in 24 regions worldwide. Regional average status relative to an untrawled state (=1) was high (>0.9) in 15 regions, but <0.7 in three (European) regions and only 0.25 in the Adriatic Sea. Across all regions, 66% of seabed area was not trawled (status = 1), 1.5% was depleted (status = 0), and 93% had status > 0.8. These assessments are first order, based on parameters estimated with uncertainty from meta-analyses; we recommend regional analyses to refine parameters for local specificity. Nevertheless, our results are sufficiently robust to highlight regions needing more effective management to reduce exploitation and improve stock sustainability and seabed environmental status-while also showing seabed status was high (>0.95) in regions where catches of trawled fish stocks meet accepted benchmarks for sustainable exploitation, demonstrating that environmental benefits accrue from effective fisheries management. Furthermore, regional seabed status was related to the proportional area swept by trawling, enabling preliminary predictions of regional status when only the total amount of trawling is known. This research advances seascape-scale understanding of trawl impacts in regions around the world, enables quantitative assessment of sustainability risks, and facilitates implementation of an ecosystem approach to trawl fisheries management globally.

North Sea demersal fisheries prefer specific benthic habitats.

INTRODUCTION: The future protection of marine biodiversity through good conservation planning requires both the identification of key habitats with unique ecological characteristics and detailed knowledge of their human utilization through fisheries. Demersal fisheries are important disturbers of benthic habitats. They often have a heterogeneous spatial distribution, pressurizing particular habitats with high abundances of target species. For the North Sea, we quantified the commonness/rarity of habitats in relation to the environmental determinants of so-called fishing hotspots, to support better-informed conservation planning of benthic habitats in this intensively used continental shelf. METHODS: We first distinguished 9 main seascapes in the study area based on seabed morphology. Secondly, we determined average fishing intensity and fishing hotspots using VMS-data for the three dominant Dutch fisheries from 2008 to 2015: beam-trawlers targeting sole Solea solea (Beam-Sole), beam-trawlers targeting plaice Pleuronectes platessa (Beam-Plaice), and otter-trawlers targeting Norway lobster Nephrops norvegicus and demersal fish (Otter-Mix). Within the seascapes subjected to >80% of the fishing activity, nineteen environmental factors (summarized by PCA) were used to ecologically characterize fishing hotspot locations using MaxEnt response modelling. RESULTS: We found that all three fisheries target highly specific, uncommon habitats. Beam-Sole fishers targeted warmer, shallow, dynamic, nearshore habitats, and within these specifically the depressions between sand ridges. Beam-Plaice fishers mainly targeted the exposed, non-muddy flanks of the Dogger Bank and similar large-scale elevations (50-75 km) where especially the ridges of smaller sand banks are used. Otter-Mix fisheries concentrated in areas with low bed shear stress, located in muddy, relatively deeper areas. IMPLICATIONS: This study is the first to provide insight in benthic habitat types that are frequently targeted by fishers in the North Sea. We demonstrated unequal exploitation pressure between seabed habitats, with the majority of hotspots in the less common habitats. Our results hence contribute to a more effective, ecologically informed planning for the protection and monitoring of all seabed habitats and biodiversity of the North Sea.

Bottom trawl fishing footprints on the world's continental shelves.

Bottom trawlers land around 19 million tons of fish and invertebrates annually, almost one-quarter of wild marine landings. The extent of bottom trawling footprint (seabed area trawled at least once in a specified region and time period) is often contested but poorly described. We quantify footprints using high-resolution satellite vessel monitoring system (VMS) and logbook data on 24 continental shelves and slopes to 1,000-m depth over at least 2 years. Trawling footprint varied markedly among regions: from <10% of seabed area in Australian and New Zealand waters, the Aleutian Islands, East Bering Sea, South Chile, and Gulf of Alaska to >50% in some European seas. Overall, 14% of the 7.8 million-km2 study area was trawled, and 86% was not trawled. Trawling activity was aggregated; the most intensively trawled areas accounting for 90% of activity comprised 77% of footprint on average. Regional swept area ratio (SAR; ratio of total swept area trawled annually to total area of region, a metric of trawling intensity) and footprint area were related, providing an approach to estimate regional trawling footprints when high-resolution spatial data are unavailable. If SAR was ≤0.1, as in 8 of 24 regions, there was >95% probability that >90% of seabed was not trawled. If SAR was 7.9, equal to the highest SAR recorded, there was >95% probability that >70% of seabed was trawled. Footprints were smaller and SAR was ≤0.25 in regions where fishing rates consistently met international sustainability benchmarks for fish stocks, implying collateral environmental benefits from sustainable fishing.

Estimating sensitivity of seabed habitats to disturbance by bottom trawling based on the longevity of benthic fauna.

Bottom fishing such as trawling and dredging may pose serious risks to the seabed and benthic habitats, calling for a quantitative assessment method to evaluate the impact and guide management to develop mitigation measures. We provide a method to estimate the sensitivity of benthic habitats based on the longevity composition of the invertebrate community. We hypothesize that long-lived species are more sensitive to trawling mortality due to their lower pace of life (i.e., slower growth, late maturation). We analyze data from box-core and grab samples taken from 401 stations in the English Channel and southern North Sea to estimate the habitat-specific longevity composition of the benthic invertebrate community and of specific functional groups (i.e., suspension feeders and bioturbators), and examine how bottom trawling affects the longevity biomass composition. The longevity biomass composition differed between habitats governed by differences in sediment composition (gravel and mud content) and tidal bed-shear stress. The biomass proportion of long-lived species increased with gravel content and decreased with mud content and shear stress. Bioturbators had a higher median longevity than suspension feeders. Trawling, in particular by gears that penetrate the seabed >2 cm, shifted the community toward shorter-lived species. Changes from bottom trawling were highest in habitats with many long-lived species (hence increasing with gravel content, decreasing with mud content). Benthic communities in high shear stress habitats were less affected by bottom trawling. Using these relationships, we predicted the sensitivity of the benthic community from bottom trawling impact at large spatial scale (the North Sea). We derived different benthic sensitivity metrics that provide a basis to estimate indicators of trawling impact on a continuous scale for the total community and specific functional groups. In combination with high resolution data of trawling pressure, our approach can be used to monitor and assess trawling impact and seabed status at the scale of the region or broadscale habitat and to compare the environmental impact of bottom-contacting fishing gears across fisheries.

Global analysis of depletion and recovery of seabed biota after bottom trawling disturbance.

Bottom trawling is the most widespread human activity affecting seabed habitats. Here, we collate all available data for experimental and comparative studies of trawling impacts on whole communities of seabed macroinvertebrates on sedimentary habitats and develop widely applicable methods to estimate depletion and recovery rates of biota after trawling. Depletion of biota and trawl penetration into the seabed are highly correlated. Otter trawls caused the least depletion, removing 6% of biota per pass and penetrating the seabed on average down to 2.4 cm, whereas hydraulic dredges caused the most depletion, removing 41% of biota and penetrating the seabed on average 16.1 cm. Median recovery times posttrawling (from 50 to 95% of unimpacted biomass) ranged between 1.9 and 6.4 y. By accounting for the effects of penetration depth, environmental variation, and uncertainty, the models explained much of the variability of depletion and recovery estimates from single studies. Coupled with large-scale, high-resolution maps of trawling frequency and habitat, our estimates of depletion and recovery rates enable the assessment of trawling impacts on unprecedented spatial scales.

Using marine reserves to manage impact of bottom trawl fisheries requires consideration of benthic food-web interactions.

Marine protected areas (MPAs) are widely used to protect exploited fish species as well as to conserve marine habitats and their biodiversity. They have also become a popular management tool for bottom trawl fisheries, a common fishing technique on continental shelves worldwide. The effects of bottom trawling go far beyond the impact on target species, as trawls also affect other components of the benthic ecosystem and the seabed itself. This means that for bottom trawl fisheries, MPAs can potentially be used not only to conserve target species but also to reduce impact of these side effects of the fishery. However, predicting the protective effects of MPAs is complicated because the side effects of trawling potentially alter the food-web interactions between target and non-target species. These changes in predatory and competitive interactions among fish and benthic invertebrates may have important ramifications for MPAs as tools to manage or mitigate the effects of bottom trawling. Yet, in current theory regarding the functioning of MPAs in relation to bottom trawl fisheries, such predatory and competitive interactions between species are generally not taken into account. In this study, we discuss how food-web interactions that are potentially affected by bottom trawling may alter the effectiveness of MPAs to protect (1) biodiversity and marine habitats, (2) fish populations, (3) fisheries yield, and (4) trophic structure of the community. We make the case that in order to be applicable for bottom trawl fisheries, guidelines for the implementation of MPAs must consider their potential food-web effects, at the risk of failing management.

Conservation physiology of marine fishes: state of the art and prospects for policy.

The state of the art of research on the environmental physiology of marine fishes is reviewed from the perspective of how it can contribute to conservation of biodiversity and fishery resources. A major constraint to application of physiological knowledge for conservation of marine fishes is the limited knowledge base; international collaboration is needed to study the environmental physiology of a wider range of species. Multifactorial field and laboratory studies on biomarkers hold promise to relate ecophysiology directly to habitat quality and population status. The 'Fry paradigm' could have broad applications for conservation physiology research if it provides a universal mechanism to link physiological function with ecological performance and population dynamics of fishes, through effects of abiotic conditions on aerobic metabolic scope. The available data indicate, however, that the paradigm is not universal, so further research is required on a wide diversity of species. Fish physiologists should interact closely with researchers developing ecological models, in order to investigate how integrating physiological information improves confidence in projecting effects of global change; for example, with mechanistic models that define habitat suitability based upon potential for aerobic scope or outputs of a dynamic energy budget. One major challenge to upscaling from physiology of individuals to the level of species and communities is incorporating intraspecific variation, which could be a crucial component of species' resilience to global change. Understanding what fishes do in the wild is also a challenge, but techniques of biotelemetry and biologging are providing novel information towards effective conservation. Overall, fish physiologists must strive to render research outputs more applicable to management and decision-making. There are various potential avenues for information flow, in the shorter term directly through biomarker studies and in the longer term by collaborating with modellers and fishery biologists.

Warming temperatures and smaller body sizes: synchronous changes in growth of North Sea fishes.

Decreasing body size has been proposed as a universal response to increasing temperatures. The physiology behind the response is well established for ectotherms inhabiting aquatic environments: as higher temperatures decrease the aerobic capacity, individuals with smaller body sizes have a reduced risk of oxygen deprivation. However, empirical evidence of this response at the scale of communities and ecosystems is lacking for marine fish species. Here, we show that over a 40-year period six of eight commercial fish species in the North Sea examined underwent concomitant reductions in asymptotic body size with the synchronous component of the total variability coinciding with a 1-2 °C increase in water temperature. Smaller body sizes decreased the yield-per-recruit of these stocks by an average of 23%. Although it is not possible to ascribe these phenotypic changes unequivocally to temperature, four aspects support this interpretation: (i) the synchronous trend was detected across species varying in their life history and life style; (ii) the decrease coincided with the period of increasing temperature; (iii) the direction of the phenotypic change is consistent with physiological knowledge; and (iv) no cross-species synchrony was detected in other species-specific factors potentially impacting growth. Our findings support a recent model-derived prediction that fish size will shrink in response to climate-induced changes in temperature and oxygen. The smaller body sizes being projected for the future are already detectable in the North Sea.

Evolutionary impact assessment: accounting for evolutionary consequences of fishing in an ecosystem approach to fisheries management.

Managing fisheries resources to maintain healthy ecosystems is one of the main goals of the ecosystem approach to fisheries (EAF). While a number of international treaties call for the implementation of EAF, there are still gaps in the underlying methodology. One aspect that has received substantial scientific attention recently is fisheries-induced evolution (FIE). Increasing evidence indicates that intensive fishing has the potential to exert strong directional selection on life-history traits, behaviour, physiology, and morphology of exploited fish. Of particular concern is that reversing evolutionary responses to fishing can be much more difficult than reversing demographic or phenotypically plastic responses. Furthermore, like climate change, multiple agents cause FIE, with effects accumulating over time. Consequently, FIE may alter the utility derived from fish stocks, which in turn can modify the monetary value living aquatic resources provide to society. Quantifying and predicting the evolutionary effects of fishing is therefore important for both ecological and economic reasons. An important reason this is not happening is the lack of an appropriate assessment framework. We therefore describe the evolutionary impact assessment (EvoIA) as a structured approach for assessing the evolutionary consequences of fishing and evaluating the predicted evolutionary outcomes of alternative management options. EvoIA can contribute to EAF by clarifying how evolution may alter stock properties and ecological relations, support the precautionary approach to fisheries management by addressing a previously overlooked source of uncertainty and risk, and thus contribute to sustainable fisheries.

When does fishing lead to more fish? Community consequences of bottom trawl fisheries in demersal food webs.

Bottom trawls are a globally used fishing gear that physically disturb the seabed and kill non-target organisms, including those that are food for the targeted fish species. There are indications that ensuing changes to the benthic invertebrate community may increase the availability of food and promote growth and even fisheries yield of target fish species. If and how this occurs is the subject of ongoing debate, with evidence both in favour and against. We model the effects of trawling on a simple ecosystem of benthivorous fish and two food populations (benthos), susceptible and resistant to trawling. We show that the ecosystem response to trawling depends on whether the abundance of benthos is top-down or bottom-up controlled. Fishing may result in higher fish abundance, higher (maximum sustainable) yield and increased persistence of fish when the benthos which is the best-quality fish food is also more resistant to trawling. These positive effects occur in bottom-up controlled systems and systems with limited impact of fish feeding on benthos, resembling bottom-up control. Fishing leads to lower yields and fish persistence in all configurations where susceptible benthos are more profitable prey. Our results highlight the importance of mechanistic ecosystem knowledge as a requirement for successful management.

Thermal preference of juvenile Dover sole (Solea solea) in relation to thermal acclimation and optimal growth temperature.

Dover sole (Solea solea) is an obligate ectotherm with a natural thermal habitat ranging from approximately 5 to 27°C. Thermal optima for growth lie in the range of 20 to 25°C. More precise information on thermal optima for growth is needed for cost-effective Dover sole aquaculture. The main objective of this study was to determine the optimal growth temperature of juvenile Dover sole (Solea solea) and in addition to test the hypothesis that the final preferendum equals the optimal growth temperature. Temperature preference was measured in a circular preference chamber for Dover sole acclimated to 18, 22 and 28°C. Optimal growth temperature was measured by rearing Dover sole at 19, 22, 25 and 28°C. The optimal growth temperature resulting from this growth experiment was 22.7°C for Dover sole with a size between 30 to 50 g. The temperature preferred by juvenile Dover sole increases with acclimation temperature and exceeds the optimal temperature for growth. A final preferendum could not be detected. Although a confounding effect of behavioural fever on temperature preference could not be entirely excluded, thermal preference and thermal optima for physiological processes seem to be unrelated in Dover sole.

Estimating age at maturation and energy-based life-history traits from individual growth trajectories with nonlinear mixed-effects models.

A new method is presented to estimate individuals' (1) age at maturation, (2) energy acquisition rate, (3) energy expenditure for body maintenance, and (4) reproductive investment, and the multivariate distribution of these traits in a population. The method relies on adjusting a conceptual energy allocation model to individual growth curves using nonlinear mixed-effects modelling. The method's performance was tested using simulated growth curves for a range of life-history types. Individual age at maturation, energy acquisition rate and the sum of maintenance and reproductive investment rates, and their multivariate distribution, were accurately estimated. For the estimation of maintenance and reproductive investment rates separately, biases were observed for life-histories with a large imbalance between these traits. For low reproductive investment rates and high maintenance rates, reproductive investment rate estimates were strongly biased whereas maintenance rate estimates were not, the reverse holding in the opposite situation. The method was applied to individual growth curves back-calculated from otoliths of North Sea plaice (Pleuronectes platessa) and from scales of Norwegian spring spawning herring (Clupea harengus). For plaice, maturity ogives derived from our individual estimates of age at maturation were almost identical to the maturity ogives based on gonad observation in catch samples. For herring, we observed 51.5% of agreement between our individual estimates and those directly obtained from scale reading, with a difference lower than 1 year in 97% of cases. We conclude that the method is a powerful tool to estimate the distribution of correlated life-history traits for any species for which individual growth curves are available.

Multiple paternity analysis in the thornback ray Raja clavata L.

Skates (Rajidae) are characterized by slow growth rate, low fecundity, and late maturity and are thus considered to be vulnerable to exploitation. Although understanding mating systems and behavior are important for long-term conservation and fisheries management, this aspect of life history is poorly understood in skates. Using 5 highly polymorphic microsatellite loci, we analyzed egg clutches collected from 4 female Raja clavata captured in the wild to test for multiple paternity. Using the reconstructed multilocus genotypes method to explain the progeny genotype array, we showed that all 4 clutches were sired by a minimum of 4-6 fathers and, thus, female thornback rays are polyandrous. Whether polyandry in R. clavata is natural or a consequence of overexploitation remains uncertain. This is the first report of multiple paternity in a rajiform species and any oviparous elasmobranch.

Phylogeography and population structure of thornback rays (Raja clavata L., Rajidae).

The phylogeography of thornback rays (Raja clavata) was assessed from European waters, using five nuclear microsatellite loci and mitochondrial cytochome b sequences. Strong regional differentiation was found between the Mediterranean basin, the Azores and the European continental shelf. Allelic and haplotype diversities were high in Portuguese populations, consistent with the existence of a refugium along the Iberian Peninsula. Unexpectedly, high diversity was also found in the English Channel/North Sea area. The lowest genetic diversity was found in the Black Sea. Populations sampled from the Mediterranean, Adriatic and Black Seas were characterized by a single mitochondrial haplotype. This haplotype was also the most ancestral and widespread outside of the Mediterranean basin except for the Azores. Populations from the Azores were dominated by a second ancestral haplotype which was shared with British populations. Results from multidimensional scaling, amova and nested clade analysis indicate that British waters are a secondary contact zone recolonized from at least two refugia--one around the Iberian Peninsula and one possibly in the Azores. Links to a potential refugium known as the Hurd Deep, between Cornwall and Brittany, are discussed. Finally, a historical demographic analysis indicates that thornback ray populations started to expand between 580,000 and 362,000 years ago, which suggests that the Last Glacial Maximum (20,000 years ago) had mainly affected the distribution of populations rather than population size.

Low effective population size and evidence for inbreeding in an overexploited flatfish, plaice (Pleuronectes platessa L.).

Overexploitation and subsequent collapses of major worldwide fisheries has made it clear that marine stocks are no inexhaustible. Unfortunately, the perception remains that marine fished are resilient to large population reductions, as even a commercially 'collapsed' stock will still consist of millions of individuals. Coupled with this notion is the idea that fisheries can, therefore, have little effect on the genetic diversity of stocks. We used DNA from archived otoliths collected between 1924 and 1972 together with 2002 juvenile;s tissue to estimate effective population size (Ne) in plaice (Pleuronrctes platessa). Ne was estimated at 20,000 in the North Sea and 2000 in Iceland. These values are five orders of magnitude smaller than the estimated census size foe the two locations. Populations examined between 1924 and 1960 were in Hardy-Weinberg equilibrium, whereas populations examined after 1970 were not. Extensive testing was performed to rule out genotyping artefacts and Wahlund effects. The significant heterozygote deficiencies found from 1970 onward were attributed to inbreeding. The emergence of inbreeding between 1905 and 19070 coincides with the increase in fishing mortality after World War II. Although the biological mechanisms remain speculative, our demonstration of inbreeding signals the need for understanding the social and mating behaviour in commercially important fishes.