Recruitment collapse and population structure of the European eel shaped by local ocean current dynamics

Productive foraging grounds enhance maternal condition and offspring quality in a capital breeding species

Feeding ecology is an essential component of an organism's life, but foraging comes with risks and energetic costs. Species in which populations exhibit more than one feeding strategy, such as sea turtles, are good systems for investigating how feeding ecology impacts life-history traits, reproduction and carried over effects across generations. Here, we investigated how the feeding ecology of loggerhead sea turtles (Caretta caretta) nesting at the Cabo Verde archipelago correlates with reproductive outputs and offspring quality. We determined the feeding ecology of female turtles before and during the breeding season from stable isotope analysis of carbon and nitrogen and correlated isotopic ratio with female and offspring traits. We found that female turtles feeding at higher trophic positions produced larger clutches. We also found that females with higher δ13C values, typical of productive foraging areas, had greater fat reserves, were less likely to be infected by leech parasites and produced heavier offspring. The offspring of infected mothers with higher δ13C values performed best in crawling and self-righting trials than those of non-infected mothers with higher δ13C values. This study shows adult female loggerheads that exploit productive areas build capital reserves that impact their reproductive success and multiple proxies for offspring quality. Overall, our findings provide valuable insights into the complex interplay between feeding ecology and reproductive success, and reveal the transgenerational carry-over effects of both feeding ecology and health on offspring quality in sea turtles.

Large biomass reduction effect on the relative role of climate, fishing, and recruitment on fish population dynamics

Many species around the world have collapsed, yet only some have recovered. A key question is what happens to populations post collapse. Traditionally, marine fish collapses are linked to overfishing, poor climate, and recruitment. We test whether the effect on biomass change from these drivers remains the same after a collapse. We used a regression model to analyse the effect of harvesting, recruitment, and climate variability on biomass change before and after a collapse across 54 marine fish populations around the world. The most salient result was the change in fishing effect that became weaker after a collapse. The change in sea temperature and recruitment effects were more variable across systems. The strongest changes were in the pelagic habitats. The resultant change in the sensitivity to external drivers indicates that whilst biomass may be rebuilt, the responses to variables known to affect stocks may have changed after a collapse. Our results show that a general model applied to many stocks provides useful insights, but that not all stocks respond similarly to a collapse calling for stock-specific models. Stocks respond to environmental drivers differently after a collapse, so caution is needed when using pre-collapse knowledge to advise on population dynamics and management.

Genome-wide methylation in the panmictic European eel (Anguilla anguilla)

The role of methylation in adaptive, developmental and speciation processes has attracted considerable interest, but interpretation of results is complicated by diffuse boundaries between genetic and non-genetic variation. We studied whole genome genetic and methylation variation in the European eel, distributed from subarctic to subtropical environments, but with panmixia precluding genetically based local adaptation beyond single-generation responses. Overall methylation was 70.9%, with hypomethylation predominantly found in promoters and first exons. Redundancy analyses involving juvenile glass eels showed 0.06% and 0.03% of the variance at SNPs to be explained by localities and environmental variables, respectively, with GO terms of genes associated with outliers primarily involving neural system functioning. For CpGs 2.98% and 1.36% of variance was explained by localities and environmental variables. Differentially methylated regions particularly included genes involved in developmental processes, with hox clusters featuring prominently. Life stage (adult versus glass eels) was the most important source of inter-individual variation in methylation, likely reflecting both ageing and developmental processes. Demethylation of transposable elements relative to pure European eel was observed in European X American eel hybrids, possibly representing postzygotic barriers in this system characterized by prolonged speciation and ongoing gene flow. Whereas the genetic data are consistent with a role of single-generation selective responses, the methylation results underpin the importance of epigenetics in the life cycle of eels and suggests interactions between local environments, development and phenotypic variation mediated by methylation variation. Eels are remarkable by having retained eight hox clusters, and the results suggest important roles of methylation at hox genes for adaptive processes.

Potential Effect of the Intrusion of the Kuroshio Current into the South China Sea on Catches of Japanese Eel (Anguilla Japonica) in the South China Sea and Taiwan Strait

This study examined the effect of the intrusion of the Kuroshio Current (KC) into the South China Sea (SCS) and the Taiwan Strait (TS) (SCS–TS region) on changes in catches of larval A. japonica in the traditional fishing ground waters of Gaoping near southwestern Taiwan in the SCS–TS region. First, the oceanic environment and recruitment trends from 1967 to 2019 were investigated based on secondary data. Then, field surveys were conducted to obtain primary data regarding the intrusion of the KC into the SCS, as well as the changes in the fishing sites and catches of A. japonica in the fall and winter of 2014–2015. Hence, the association between oceanic conditions and the number of A. japonica migrating into the SCS–TS region was explored. From 1967 to 2019, the recruitment proportion in the fishing grounds that formed due to the Kuroshio Branch Current (PKSBC) fluctuated significantly. Overall, positive values were observed for the Oceanic Niño Index for each year with a PKSBC > 50%, corresponding to El Niño conditions. In each year with a PKSBC > 70%, a looping path and a warm–core eddy appeared.

New clues on the Atlantic eels spawning behavior and area: the Mid-Atlantic Ridge hypothesis

The Sargasso Sea has long been considered as the only spawning area for Atlantic eels, despite the absence of direct observations. The present study raises a novel scenario, deviating from Schmidt’s dogma, begins with a review of historical and recent observations that were combined to build up a global theory on spawning ecology and migration behavior of Atlantic eels. From this, it is argued that a favorable spawning area could be located eastward of Sargasso Sea at the intersection between the Mid-Atlantic Ridge and the oceanic fronts. Ocean circulation models combined with 3D particle-tracking method confirmed that spawning at this specific area would result in larval distribution fitting the field observation. This study explores the hypothesis that leptocephali are able to swim and orientate to reach their specific growth areas. It proposes a novel framework about spawning ecology, based on orientation, navigation and meeting cues of silver eels to the spawning area. Together this framework may serve as a stepping-stone for solving the long-lasting mystery of eel reproduction which first came out 2,400 years ago and promotes the understanding of oceanic migration and reproduction of marine organisms.

Microbial tropicalization driven by a strengthening western ocean boundary current

Western boundary currents (WBCs) redistribute heat and oligotrophic seawater from the tropics to temperate latitudes, with several displaying substantial climate change-driven intensification over the last century. Strengthening WBCs have been implicated in the poleward range expansion of marine macroflora and fauna, however, the impacts on the structure and function of temperate microbial communities are largely unknown. Here we show that the major subtropical WBC of the South Pacific Ocean, the East Australian Current (EAC), transports microbial assemblages that maintain tropical and oligotrophic (k-strategist) signatures, to seasonally displace more copiotrophic (r-strategist) temperate microbial populations within temperate latitudes of the Tasman Sea. We identified specific characteristics of EAC microbial assemblages compared with non-EAC assemblages, including strain transitions within the SAR11 clade, enrichment of Prochlorococcus, predicted smaller genome sizes and shifts in the importance of several functional genes, including those associated with cyanobacterial photosynthesis, secondary metabolism and fatty acid and lipid transport. At a temperate time-series site in the Tasman Sea, we observed significant reductions in standing stocks of total carbon and chlorophyll a, and a shift towards smaller phytoplankton and carnivorous copepods, associated with the seasonal impact of the EAC microbial assemblage. In light of the substantial shifts in microbial assemblage structure and function associated with the EAC, we conclude that climate-driven expansions of WBCs will expand the range of tropical oligotrophic microbes, and potentially profoundly impact the trophic status of temperate waters.

A mechanical approach to understanding the impact of the nematode Anguillicoloides crassus on the European eel swimbladder

One of the most detrimental factors in the drastic decline of the critically endangered European eel (Anguilla anguilla) was the inadvertent introduction of the invasive nematode Anguillicoloides crassus Infection primarily affects the swimbladder, a gas-filled organ that enables the eel to control its depth in the water. A reduction in swimbladder function may be fatal for eel undergoing their spawning migration to the Sargasso Sea, a journey of over 5000 km. Although the physiological damage caused by this invasive parasite is well studied through the use of quantifiable gross pathological indices, providing a good measure of the swimbladder health status, they cannot separate the role of mechanical and morphological damage. Our study examined the appropriateness of three commonly used indices as a measure of mechanical damage by performing uniaxial tensile tests on swimbladder specimens obtained from an infected eel population. When the test results were compared with the gross pathological indices it was found that thickness correlated most strongly with mechanical damage, both confirming and, more importantly, explaining the counterintuitive findings of earlier work. In a damaged swimbladder, the immune response leads to a trade-off; increasing wall thickness raises the pressure required for organ rupture but decreases strength. The results indicate that for moderate infection the mechanical integrity of the swimbladder can be maintained. For severe infection, however, a reduction in mechanical integrity may reach a tipping point, thereby affecting the successful completion of their oceanic migration.

Larval behaviour, dispersal and population connectivity in the deep sea

Ecosystem connectivity is an essential consideration for marine spatial planning of competing interests in the deep sea. Immobile, adult communities are connected through freely floating larvae, depending on new recruits for their health and to adapt to external pressures. We hypothesize that the vertical swimming ability of deep-sea larvae, before they permanently settle at the bottom, is one way larvae can control dispersal. We test this hypothesis with more than \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$3\times 10^{8}$$\end{document}3×108 simulated particles with a range of active swimming behaviours embedded within the currents of a high-resolution ocean model. Despite much stronger horizontal ocean currents, vertical swimming of simulated larvae can have an order of magnitude impact on dispersal. These strong relationships between larval dispersal, pathways, and active swimming demonstrate that lack of data on larval behaviour traits is a serious impediment to modelling deep-sea ecosystem connectivity; this uncertainty greatly limits our ability to develop ecologically coherent marine protected area networks.

Spatio-temporal genetic tagging of a cosmopolitan planktivorous shark provides insight to gene flow, temporal variation and site-specific re-encounters

Migratory movements in response to seasonal resources often influence population structure and dynamics. Yet in mobile marine predators, population genetic consequences of such repetitious behaviour remain inaccessible without comprehensive sampling strategies. Temporal genetic sampling of seasonally recurring aggregations of planktivorous basking sharks, Cetorhinus maximus, in the Northeast Atlantic (NEA) affords an opportunity to resolve individual re-encounters at key sites with population connectivity and patterns of relatedness. Genetic tagging (19 microsatellites) revealed 18% of re-sampled individuals in the NEA demonstrated inter/multi-annual site-specific re-encounters. High genetic connectivity and migration between aggregation sites indicate the Irish Sea as an important movement corridor, with a contemporary effective population estimate (Ne) of 382 (CI = 241-830). We contrast the prevailing view of high gene flow across oceanic regions with evidence of population structure within the NEA, with early-season sharks off southwest Ireland possibly representing genetically distinct migrants. Finally, we found basking sharks surfacing together in the NEA are on average more related than expected by chance, suggesting a genetic consequence of, or a potential mechanism maintaining, site-specific re-encounters. Long-term temporal genetic monitoring is paramount in determining future viability of cosmopolitan marine species, identifying genetic units for conservation management, and for understanding aggregation structure and dynamics.

Ocean currents and the population genetic signature of fish migrations

Animal migrations are a fascinating and global phenomenon, yet they are often difficult to study and sometimes poorly understood. Here, we build on classic ecological theory by hypothesizing that some enigmatic spawning migrations across coastal marine habitats can be inferred from the population genetic signature of larval dispersal by ocean currents. We test this assumption by integrating spatially realistic simulations of alternative spawning migration routes, associated patterns of larval dispersal, and associated variation in the population genetic structure of eastern Australian sea mullet (Mugil cephalus). We then use simulation results to assess the implications of alternative spawning destinations for larval replenishment, and we contrast simulated against measured population genetic variation. Both analyses suggest that the spawning migrations of M. cephalus in eastern Australia are likely to be localized (approximately 100 km along the shore), and that spawning is likely to occur in inshore waters. Our conclusions are supported by multiple lines of evidence available through independent studies, but they challenge the more traditional assumption of a single, long-distance migration event with subsequent offshore spawning in the East Australian Current. More generally, our study operationalizes classic theory on the relationship between fish migrations, ocean currents, and reproductive success. However, rather than confirming the traditionally assumed adaptation of migratory behavior to dominant ocean current flow, our findings support the concept of a genetically measurable link between fish migrations and local oceanographic conditions, specifically water temperature and coastal retention of larvae. We believe that future studies using similar approaches for high resolution and spatially realistic ecological-genetic scenario testing can help rapidly advance our understanding of key ecological processes in many other marine species.

Influence of ocean circulation and the Kuroshio large meander on the 2018 Japanese eel recruitment season

Japanese eel (Anguilla japonica) recruitment to Japan was very low during the early 2017−2018 recruitment season when most glass eels are usually caught, but catches increased in the late recruitment season when recruitment usually decreases. Concurrently, the Kuroshio meander south of Japan had formed again after the previous event ended in 2005. The role of the large meander and ocean circulation features such as the North Equatorial Current (NEC) in the unusual 2017−2018 Japanese eel recruitment timing-pattern was investigated using a three-dimensional particle tracking model that simulated swimming behaviors of virtual larvae (v-larvae) in addition to their drift in ocean currents. Four recruitment seasons were selected for when the Kuroshio large meander was present (2004−2005, 2017−2018) or absent (2009−2010, 2015−2016), and when NEC was shifted north (2004−2005, 2015−2016) or south (2009−2010, 2017−2018). The simulated recruitment timing-patterns were similar to the actual recent-year recruitment, with no early recruitment period v-larvae arrival to southern Japan and increased late period recruitment occurring. Rather than being related to the presence of the Kuroshio large meander, the late arrival appeared to be caused by a southward shifted, weak North Equatorial Current near the spawning area, a longer Subtropical Countercurrent eddy region retention time, and a weak Kuroshio during the early migration and recruitment period of those years. In the late recruitment period, the Kuroshio was stronger than other selected years near the East China Sea and south of Japan and v-larvae were transported more rapidly. The Kuroshio large meander may influence local eel recruitment in Japan, and the recirculation formed by the large meander could potentially enhance recruitment to the Tokai region. Oriented (northwestward) swimming v-larvae were less affected by the large meander, and showed higher recruitment success than those using along-current swimming. Although the Kuroshio large meander did not seem to be responsible for the unusual recruitment pattern in 2018, how Japanese eel larvae and glass eels actually cross out of the Kuroshio and reach coastal waters in Japan remains to be explored.

Animal Navigation: The Eel's Magnetic Guide to the Gulf Stream

The geographic distribution of migratory species can span thousands of kilometers. Yet, traits that enable large-scale migrations are poorly understood. A recent study demonstrates that juvenile eels use the Earth's magnetism for their dispersal, with possible implications for their evolution.

Effect of larval swimming in the western North Pacific subtropical gyre on the recruitment success of the Japanese eel

The possible effect of directional larval swimming on the recruitment success of the Japanese eel, Anguilla japonica, was examined with a three-dimensional particle-tracking ocean circulation model using horizontal northwestward swimming and diel vertical migration (DVM). Four separate experiments included virtual larvae (v-larvae) movement from the spawning area over 290 days (total migration) and 160 days (stage A), from the STCC eddy region in 70 days (stage B), and from the origin of the Kuroshio in 60 days (stage C) to evaluate the effect of directional swimming and DVM compared to simple drifting. Passive or random swimming were not the most effective strategies for larvae dispersing from the spawning area because most v-larvae remained south of 20°N without entering the Kuroshio. Northwestward swimming resulted in wider dispersion and a better chance of successful recruitment, with v-larvae becoming widely distributed in the STCC eddy zone, arriving at the east coast of the Philippines (stage A), escaping the STCC eddy area and reaching the Kuroshio (stage B), and crossing the Kuroshio into the East China Sea shelf (stage C). DVM slightly shortened the migration period due to faster shallow layer ocean currents during nighttime. The NEC transported non-swimming v-larvae westward to the Kuroshio and occasionally northward into the Subtropical Countercurrent (STCC) area where eddies transported v-larvae westward into the Kuroshio, but less than with swimming. Directional swimming increased recruitment success, northwestward swimming was more effective than other directions, and a slower swimming speed was still better than no/random swimming in sensitivity tests. The present study demonstrated a first view of the possibility that Japanese eel larvae might be able to use a strategy of single-direction swimming to increase arrival at their recruitment areas.

Whole-Genome Sequencing of 84 Japanese Eels Reveals Evidence against Panmixia and Support for Sympatric Speciation

The Japanese eel (Anguilla japonica), European eel (Anguilla anguilla), and American eel (Anguilla rostrata) are migratory, catadromous, temperate zone fish sharing several common life cycle features. The population genetics of panmixia in these eel species has already been investigated. Our extensive population genetics analysis was based on 1400 Gb of whole-genome sequence (WGS) data from 84 eels. It demonstrated that a Japanese eel group from the Kuma River differed from other populations of the same species. Even after removing the potential adapted/selected single nucleotide polymorphism (SNP) data, and with very small differences (fixation index [Fst] = 0.01), we obtained results consistently indicating that panmixia does not occur in Japanese eels. The life cycle of the Japanese eel is well-established and the Kuma River is in the center of its habitat. Nevertheless, simple reproductive isolation is not the probable cause of non-panmixia in this species. We propose that the combination of spawning area subdivision, philopatry, and habitat preference/avoidance accounts for the non-panmixia in the Japanese eel population. We named this hypothesis the “reproductive isolation like subset mapping” (RISM) model. This finding may be indicative of the initial stages of sympatric speciation in these eels.

Head shape disparity impacts pollutant accumulation in European eel

Several aspects of the life cycle of the critically endangered European eel (Anguilla anguilla) remain poorly understood. One such aspect is the broad-versus narrow-head dimorphism, and how this impacts their overall performance at different stages of their life cycle. At the yellow eel stage, the phenotypes show a trophic divergence. We investigated whether pollutant accumulation is affected by this disparity. We show that broad-headed eels contained higher concentrations of mercury and several lipophilic organic pollutants, compared to narrow-headed ones, irrespective of their fat content. The hereby confirmed link between the phenotypic disparity, its associated feeding ecology and its impact on pollutant accumulation thus raises further concerns about their migratory and reproductive success. Considering that pollution is an important contributor to the European eel's decline, our results demonstrate that broad-headed eels are more vulnerable to detrimental pollutant accumulation. This compromises their successful contribution to their population's reproduction and its restoration.

Potential impact of ocean circulation on the declining Japanese eel catches

Recruitment of Japanese eels, Anguilla japonica, has declined in recent decades possibly due to both anthropogenic and ocean-atmosphere factors. The potential impact of ocean circulation on the decreasing Japanese eel catches in the western North Pacific was examined based on a three-dimensional particle-tracking method, in which virtual larvae (v-larvae) were programmed to swim horizontally and vertically, in addition to being transported by ocean currents after being released in their North Equatorial Current (NEC) spawning area. Transport patterns varied among years between 1993 and 2013, and dispersion of v-larvae towards East Asia decreased in the last two decades, especially for the western Taiwan and Japan regions. In recent years, instead of entering the Kuroshio and moving towards East Asia as in the 1990s’, more v-larvae tended to enter the southern areas due to the weakening of the NEC and strengthening of subsurface southward flow near the spawning area. Changes in ocean circulation in the western Pacific appear to be caused by the weakening of subtropical and tropical wind stress curl in the past two decades. This suggests that decadal changes in ocean circulation have occurred that affect the larval migration success of the Japanese eel to their recruitment areas.

Anguillicola crassus infection affects mRNA expression levels in gas gland tissue of European yellow and silver eel

Using Illumina sequencing, we investigated transcriptional changes caused by the nematode Anguillicola crassus within yellow and silver eels by comparing swimbladder samples of uninfected yellow with infected yellow eels, and uninfected silver with infected silver eels, respectively. In yellow eel gas gland, the infection caused a modification of steady state mRNA levels of 1675 genes, most of them being upregulated. Functional annotation analysis based on GO terms was used to categorize identified genes with regard to swimbladder metabolism or response to the infection. In yellow eels, the most prominent category was 'immune response', including various inflammatory components, complement proteins, and immunoglobulins. The elevated expression of several glucose and monocarboxylate transporters indicated an attempt to maintain the level of glucose metabolism, even in due to the infection thickened swimbladder tissue. In silver eel swimbladder tissue, on the contrary, the mRNA levels of only 291 genes were affected. Genes in the categories 'glucose metabolism' and 'ROS metabolism' barely responded to the infection and even the reaction of the immune system was much less pronounced compared to infected yellow eels. However, in the category 'extracellular matrix', the mRNA levels of several mucin genes were strongly elevated, suggesting increased mucus production as a defense reaction against the parasite. The present study revealed a strong reaction to an Anguillicola crassus infection on mRNA expression levels in swimbladder tissue of yellow eels, whereas in silver eels the changes ware almost negligible. A possible explanation for this difference is that the silvering process requires so much energy that there is not much scope to cope with the additional challenge of a nematode infection. Another possible explanation could be that gas-secreting activity of the silver eel swimbladder was largely reduced, which could coincide with a reduced responsiveness to other challenges, like a nematode infection.

Anguillicola crassus impairs the silvering-related enhancements of the ROS defense capacity in swimbladder tissue of the European eel (Anguilla anguilla)

In a process called silvering, European eels prepare for their long-distance migration from European freshwater systems to the Sargasso Sea for reproduction. During this journey, eels perform extended diel vertical migrations, and the concomitant changes in hydrostatic pressure significantly affect the swimbladder, functioning as a buoyancy organ. As the swimbladder is primarily filled with oxygen, the tissue has to cope with extreme hyperoxic conditions, which typically are accompanied by the generation of reactive oxygen species (ROS) and oxidative stress. In addition, since the introduction of the parasitic nematode Anguillicola crassus in the early 1980s, swimbladder function of most of the European eels is impaired by the infection with this parasite. However, the exact pathways to detoxify ROS and how these pathways are affected by silvering or the infection are still unknown. In swimbladder and muscle tissue from uninfected and infected yellow, and from uninfected and infected silver eels, we measured the level of lipid peroxidation, which increases with ROS stress. To assess the capacity of the ROS defense systems, we analyzed the activities of superoxide dismutase (SOD), catalase, glutathione peroxidase (GPx) and glutathione reductase (GR), and determined the concentration of the antioxidant glutathione (GSH + GSSG). In swimbladder tissue, we found increased concentrations of GSH + GSSG as well as higher activities of SOD, GPx and GR, suggesting that SOD and the glutathione cycle are important for ROS detoxification. Comparing swimbladder tissue of uninfected yellow with uninfected silver eels, the concentration of GSH + GSSG and the activity of SOD were higher after silvering, corresponding with lower levels of lipid peroxidation. Whereas in yellow eels the infection with A. crassus had no effect, in silver eels the capacity to cope with ROS was significantly impaired. In muscle tissue, silvering or the infection only affected the activity of SOD but in exactly the same way as in swimbladder tissue.

Evolution: Ocean Models Reveal Life in Deep Seas

Even though the deep sea represents the largest area in the world, evolution of species from those environments remains largely unstudied. A series of recent papers indicate that combining molecular tools with biophysical models can help us resolve some of these deep mysteries.

Biophysical and Population Genetic Models Predict the Presence of "Phantom" Stepping Stones Connecting Mid-Atlantic Ridge Vent Ecosystems

Deep-sea hydrothermal vents are patchily distributed ecosystems inhabited by specialized animal populations that are textbook meta-populations. Many vent-associated species have free-swimming, dispersive larvae that can establish connections between remote populations. However, connectivity patterns among hydrothermal vents are still poorly understood because the deep sea is undersampled, the molecular tools used to date are of limited resolution, and larval dispersal is difficult to measure directly. A better knowledge of connectivity is urgently needed to develop sound environmental management plans for deep-sea mining. Here, we investigated larval dispersal and contemporary connectivity of ecologically important vent mussels (Bathymodiolus spp.) from the Mid-Atlantic Ridge by using high-resolution ocean modeling and population genetic methods. Even when assuming a long pelagic larval duration, our physical model of larval drift suggested that arrival at localities more than 150 km from the source site is unlikely and that dispersal between populations requires intermediate habitats ("phantom" stepping stones). Dispersal patterns showed strong spatiotemporal variability, making predictions of population connectivity challenging. The assumption that mussel populations are only connected via additional stepping stones was supported by contemporary migration rates based on neutral genetic markers. Analyses of population structure confirmed the presence of two southern and two hybridizing northern mussel lineages that exhibited a substantial, though incomplete, genetic differentiation. Our study provides insights into how vent animals can disperse between widely separated vent habitats and shows that recolonization of perturbed vent sites will be subject to chance events, unless connectivity is explicitly considered in the selection of conservation areas.

Does asymmetric gene flow among matrilines maintain the evolutionary potential of the European eel?

Using evolutionary theory to predict the dynamics of populations is one of the aims of evolutionary conservation. In endangered species, with geographic range extending over continuous areas, the predictive capacity of evolutionary-based conservation measures greatly depends on the accurate identification of reproductive units. The endangered European eel (Anguilla anguilla) is a highly migratory fish species with declining population due to a steep recruitment collapse in the beginning of the 1980s. Despite punctual observations of genetic structure, the population is viewed as a single panmictic reproductive unit. To understand the possible origin of the detected structure in this species, we used a combination of mitochondrial and nuclear loci to indirectly evaluate the possible existence of cryptic demes. For that, 403 glass eels from three successive cohorts arriving at a single location were screened for phenotypic and genetic diversity, while controlling for possible geographic variation. Over the 3 years of sampling, we consistently identified three major matrilines which we hypothesized to represent demes. Interestingly, not only we found that population genetic models support the existence of those matriline-driven demes over a completely panmictic mode of reproduction, but also we found evidence for asymmetric gene flow amongst those demes. We uphold the suggestion that the detection of demes related to those matrilines reflect a fragmented spawning ground, a conceptually plausible consequence of the low abundance that the European eel has been experiencing for three decades. Furthermore, we suggest that this cryptic organization may contribute to the maintenance of the adaptive potential of the species.

Evaluating the adaptive potential of the European eel: is the immunogenetic status recovering?

The recent increased integration of evolutionary theory into conservation programs has greatly improved our ability to protect endangered species. A common application of such theory links population dynamics and indices of genetic diversity, usually estimated from neutrally evolving markers. However, some studies have suggested that highly polymorphic adaptive genes, such as the immune genes of the Major Histocompatibility Complex (MHC), might be more sensitive to fluctuations in population dynamics. As such, the combination of neutrally- and adaptively-evolving genes may be informative in populations where reductions in abundance have been documented. The European eel (Anguilla anguilla) underwent a drastic and well-reported decline in abundance in the late 20th century and still displays low recruitment. Here we compared genetic diversity indices estimated from neutral (mitochondrial DNA and microsatellites) and adaptive markers (MHC) between two distinct generations of European eels. Our results revealed a clear discrepancy between signatures obtained for each class of markers. Although mtDNA and microsatellites showed no changes in diversity between the older and the younger generations, MHC diversity revealed a contemporary drop followed by a recent increase. Our results suggest ongoing gain of MHC genetic diversity resulting from the interplay between drift and selection and ultimately increasing the adaptive potential of the species.

Simulating the Oceanic Migration of Silver Japanese Eels

The oceanic migration of silver Japanese eels starts from their continental growth habitats in East Asia and ends at the spawning area near the West Mariana Ridge seamount chain. However, the actual migration routes remain unknown. In this study, we examined the possible oceanic migration routes and strategies of silver Japanese eels using a particle tracking method in which virtual eels (v-eels) were programmed to move vertically and horizontally in an ocean circulation model (Japan Coastal Ocean Predictability Experiment 2, JCOPE2). Four horizontal swimming strategies were tested: random heading, true navigation (readjusted heading), orientation toward the spawning area (fixed heading), and swimming against the Kuroshio. We found that all strategies, except random swimming, allowed v-eels swimming at 0.65 m s-1 to reach the spawning area within eight months after their departure from the south coast of Japan (end of the spawning season). The estimated minimum swimming speed required to reach the area spawning within eight months was 0.1 m s-1 for true navigation, 0.12 m s-1 for constant compass heading, and 0.35 m s-1 for swimming against the Kuroshio. The lowest swimming speed estimated from tracked Japanese eels at sea was 0.03 m.s-1, which would not allow them to reach the spawning area within eight months, through any of the tested orientation strategies. Our numerical experiments also showed that ocean circulation significantly affected the migration of Japanese v-eels. A strong Kuroshio could advect v-eels further eastward. In addition, western Pacific ocean currents accelerated the migration of navigating v-eels. The migration duration was shortened in years with a stronger southward flow, contributed by a stronger recirculation south of Japan, an enhanced subtropical gyre, or a higher southward Kuroshio velocity.

Impacts of Interannual Ocean Circulation Variability on Japanese Eel Larval Migration in the Western North Pacific Ocean

The Japanese eel larvae hatch near the West Mariana Ridge seamount chain and travel through the North Equatorial Current (NEC), the Kuroshio, and the Subtropical Countercurrent (STCC) region during their shoreward migration toward East Asia. The interannual variability of circulation over the subtropical and tropical regions of the western North Pacific Ocean is affected by the Philippines–Taiwan Oscillation (PTO). This study examines the effect of the PTO on the Japanese eel larval migration routes using a three-dimensional (3D) particle tracking method, including vertical and horizontal swimming behavior. The 3D circulation and hydrography used for particle tracking are from the ocean circulation reanalysis produced by the Japan Coastal Ocean Predictability Experiment 2 (JCOPE2). Our results demonstrate that bifurcation of the NEC and the strength and spatial variation of the Kuroshio affect the distribution and migration of eel larvae. During the positive phase of PTO, more virtual eels (“v-eels”) can enter the Kuroshio to reach the south coast of Japan and more v-eels reach the South China Sea through the Luzon Strait; the stronger and more offshore swing of the Kuroshio in the East China Sea leads to fewer eels entering the East China Sea and the onshore movement of the Kuroshio to the south of Japan brings the eels closer to the Japanese coast. Significant differences in eel migration routes and distributions regulated by ocean circulation in different PTO phases can also affect the otolith increment. The estimated otolith increment suggests that eel age tends to be underestimated after six months of simulation due to the cooler lower layer temperature. Underestimation is more significant in the positive PTO years due to the wide distribution in higher latitudes than in the negative PTO years.

A global viability assessment of the European eel

The global European eel (Anguilla anguilla) stock is critically endangered according to the IUCN, and the European Commission has urged the development of conservation plans aimed to ensure its viability. However, the complex life cycle of this panmictic species, which reproduces in the open ocean but spends most of its prereproductive life in continental waters (thus embracing a huge geographic range and a variety of habitat types), makes it difficult to assess the long-term effectiveness of conservation measures. The interplay between local and global stressors raises intriguing cross-scale conservation challenges that require a comprehensive modelling approach to be addressed. We developed a full life cycle model of the global European eel stock, encompassing both the oceanic and the continental phases of eel's life, and explicitly allowing for spatial heterogeneity in vital rates, availability of suitable habitat and settlement potential via a metapopulation approach. We calibrated the model against a long-term time series of global European eel catches and used it to hindcast the dynamics of the stock in the past and project it over the 21st century under different management scenarios. Although our analysis relies on a number of inevitable simplifying assumptions and on data that may not embrace the whole range of variation in population dynamics at the small spatiotemporal scale, our hindcast is consistent with the general pattern of decline of the stock over recent decades. The results of our projections suggest that (i) habitat loss played a major role in the European eel decline; (ii) the viability of the global stock is at risk if appropriate protection measures are not implemented; (iii) the recovery of spawner escapement requires that fishing mortality is significantly reduced; and (iv) the recovery of recruitment might not be feasible if reproductive output is not enhanced.

Nano-tags for neonates and ocean-mediated swimming behaviours linked to rapid dispersal of hatchling sea turtles

Dispersal during juvenile life stages drives the life-history evolution and dynamics of many marine vertebrate populations. However, the movements of juvenile organisms, too small to track using conventional satellite telemetry devices, remain enigmatic. For sea turtles, this led to the paradigm of the 'lost years' since hatchlings disperse widely with ocean currents. Recently, advances in the miniaturization of tracking technology have permitted the application of nano-tags to track cryptic organisms. Here, the novel use of acoustic nano-tags on neonate loggerhead turtle hatchlings enabled us to witness first-hand their dispersal and behaviour during their first day at sea. We tracked hatchlings distances of up to 15 km and documented their rapid transport (up to 60 m min(-1)) with surface current flows passing their natal areas. Tracking was complemented with laboratory observations to monitor swimming behaviours over longer periods which highlighted (i) a positive correlation between swimming activity levels and body size and (ii) population-specific swimming behaviours (e.g. nocturnal inactivity) suggesting local oceanic conditions drive the evolution of innate swimming behaviours. Knowledge of the swimming behaviours of small organisms is crucial to improve the accuracy of ocean model simulations used to predict the fate of these organisms and determine resultant population-level implications into adulthood.

Local scale, coastal currents influence recruitment to freshwater populations in the European eel Anguilla anguilla: a case study from the Isle of Man

This study examines juvenile Anguilla anguilla (<30 cm) abundance in five study catchments on the Isle of Man. Preliminary results suggest that juvenile abundance is negatively correlated with increasing coastal current speed at river mouth entry (P < 0·05). These findings indicate that at least under some circumstances, tidally driven coastal currents may influence recruitment to freshwater habitats; therefore, it is presumed that high coastal current speed at the entry to river mouths may reduce the likelihood of freshwater entry.

RAD sequencing highlights polygenic discrimination of habitat ecotypes in the panmictic American eel

The two primary ways that species respond to heterogeneous environments is through local adaptation and phenotypic plasticity. The American eel (Anguilla rostrata) presents a paradox; despite inhabiting drastically different environments [1], the species is panmictic [2, 3]. Spawning takes place only in the southern Sargasso Sea in the Atlantic Ocean [1]. Then, the planktonic larvae (leptocephali) disperse to rearing locations from Cuba to Greenland, and juveniles colonize either freshwater or brackish/saltwater habitats, where they spend 3-25 years before returning to the Sargasso Sea to spawn as a panmictic species. Depending on rearing habitat, individuals exhibit drastically different ecotypes [4-6]. In particular, individuals rearing in freshwater tend to grow slowly and mature older and are more likely to be female in comparison to individuals that rear in brackish/saltwater [4, 6]. The hypothesis that phenotypic plasticity alone can account for all of the differences was not supported by three independent controlled experiments [7-10]. Here, we present a genome-wide association study that demonstrates a polygenic basis that discriminates these habitat-specific ecotypes belonging to the same panmictic population. We found that 331 co-varying loci out of 42,424 initially considered were associated with the divergent ecotypes, allowing a reclassification of 89.6%. These 331 SNPs are associated with 101 genes that represent vascular and morphological development, calcium ion regulation, growth and transcription factors, and olfactory receptors. Our results are consistent with divergent natural selection of phenotypes and/or genotype-dependent habitat choice by individuals that results in these genetic differences between habitats, occurring every generation anew in this panmictic species.

Relationship between European eel Anguilla anguilla infection with non-native parasites and swimming behaviour on encountering accelerating flow

The effect of Anguillicola crassus, Pseudodactylogyrus bini and Pseudodactylogyrus anguillae infection on the behaviour of downstream migrating adult European eels Anguilla anguilla as they encountered accelerating water velocity, common at engineered structures where flow is constricted (e.g. weirs and bypass systems), was evaluated in an experimental flume. The probability of reacting to, and rejecting, the velocity gradient was positively related to A. crassus larval, adult and total abundance. High abundance of Pseudodactylogyrus spp. reduced this effect, but A. crassus was the strongest parasitic factor associated with fish behaviour, and abundance was positively related to delay in downstream passage. Delayed downstream migration at hydraulic gradients associated with riverine anthropogenic structures could result in additional energetic expenditure for migrating A. anguilla already challenged by A. crassus infection.

Swimbladder function and the spawning migration of the European eel Anguilla anguilla

The spawning migration of the European eel is an extensive journey over 5000 to 7000 km from the European coast to the Sargasso Sea. Eels do not feed during this journey and on-board fuels must be sufficient to support the journey of 3.5 to 6 month, as well as sexual maturation and the spawning activity. Swimming of eels appears to be quite energy efficient compared to other fish species, and elevated hydrostatic pressure has been shown to even reduce the costs of transport. Recent studies revealed, however, that during traveling eels perform extensive diurnal migrations and swim at a depth of about 100-300 m at night time, but go down to 600-1000 m at day time. At a depth of 200 m eels are exposed to a hydrostatic pressure of 21 atmospheres (2.13 MPa), while at 800 m hydrostatic pressure increases to 81 atmospheres (8.21 MPa). Accordingly, without any compensation at a depth of 800 m swimbladder volume will be reduced to about 25% of the volume established with neutral buoyancy at 200 m. Consequently, these diurnal changes in depth must be taken into consideration for a calculation of the energy requirements of the spawning migration. Without compensation a compression of the swimbladder will result in a status of negative buoyancy, which makes swimming more costly. Trying to keep the status of neutral buoyancy during descent by gas secretion into the swimbladder in turn requires metabolic activity to enhance swimbladder perfusion and for acid production of the gas gland cells to stimulate gas secretion. During ascent gas is passively removed from the swimbladder in the resorbing section and in the blood transported to the gills, where it is lost into the water. Accordingly, the swimbladder appears to be a crucial organ for the spawning migration. It can be assumed that an impairment of swimbladder function for example due to an infection with the nematode Anguillicola crassus significantly threatens the success of the spawning migration.

PCB and organochlorine pesticide burden in eels in the lower Thames River (UK)

Thirty-five European eels (Anguilla anguilla), caught in 2007 in the river Thames upstream and downstream of both London and the tidal limit, were analysed for PCBs and organochlorine pesticides. Most chemicals were detectable in every fish, although they have been banned or severely restricted for many years. In general, the tidal eels were more contaminated than upstream ones, which was related to their higher lipid contents. The ICES7 indicator PCB concentrations ranged overall from 4.2 to 124μgkg(-1) fresh weight with averages of 33 and 56μgkg(-1) for the upstream and tidal eels; 3.5-104μgkg(-1), average 26 and 48μgkg(-1) of that were ICES6 PCBs. Total DDT was on average 16μgkg(-1) (1.7-38μgkg(-1)) upstream and 18μgkg(-1) (8.6-35μgkg(-1)) downstream with about half of that provided by pp'DDE. Lindane (γ-HCH) was found at up to 2.8μgkg(-1) (averages 0.58 and 1.1μgkg(-1) upstream and downstream) and hexachlorobenzene (HCB) was on average 1.9 and 2.5μgkg(-1) in the two groups with a maximum of 6.4μgkg(-1) in each. Therefore all individuals passed the European Environmental Quality Standard (EQS) of 10μgkg(-1) for HCB. PCB contamination was fairly typical for recent UK eel data, whilst DDE and lindane concentrations were lower than most previous UK eel studies, perhaps reflecting a downward trend. Although not as highly contaminated as some eels from previous UK and European studies, the presence of so many of these chemicals, with their known health effects may represent a stress for the fish or higher predators, such as birds.

A century of research on the larval distributions of the Atlantic eels: a re-examination of the data

The spawning areas of the Atlantic freshwater eels were discovered about a century ago by the Danish scientist Johannes Schmidt who after years of searching found newly hatched larvae of the European eel, Anguilla anguilla, and the American eel, Anguilla rostrata, in the southern Sargasso Sea. The discovery showed that anguillid eels migrate thousands of kilometers to offshore spawning areas for reproduction, and that their larvae, called leptocephali, are transported equally long distances by ocean currents to their continental recruitment areas. The spawning sites were found to be related to oceanographic conditions several decades later by German and American surveys from 1979 to 1989 and by a Danish survey in 2007 and a German survey in 2011. All these later surveys showed that spawning occurred within a restricted latitudinal range, between temperature fronts within the Subtropical Convergence Zone of the Sargasso Sea. New data and re-examinations of Schmidt's data confirmed his original conclusions about the two species having some overlap in spawning areas. Although there have been additional collections of leptocephali in various parts of the North Atlantic, and both otolith research and transport modelling studies have subsequently been carried out, there is still a range of unresolved questions about the routes of larval transport and durations of migration. This paper reviews the history and basic findings of surveys for anguillid leptocephali in the North Atlantic and analyses a new comprehensive database that includes 22612 A. anguilla and 9634 A. rostrata leptocephali, which provides a detailed view of the spatial and temporal distributions and size of the larvae across the Atlantic basin and in the Mediterranean Sea. The differences in distributions, maximum sizes, and growth rates of the two species of larvae are likely linked to the contrasting migration distances to their recruitment areas on each side of the basin. Anguilla rostrata leptocephali originate from a more western spawning area, grow faster, and metamorphose at smaller sizes of <70 mm than the larvae of A. anguilla, which mostly are spawned further east and can reach sizes of almost 90 mm. The larvae of A. rostrata spread west and northwest from the spawning area as they grow larger, with some being present in the western Caribbean and eastern Gulf of Mexico. Larvae of A. anguilla appear to be able to reach Europe by entering the Gulf Stream system or by being entrained into frontal countercurrents that transport them directly northeastward. The larval duration of A. anguilla is suggested to be quite variable, but gaps in sampling effort prevent firm conclusions. Although knowledge about larval behaviour is lacking, some influences of directional swimming are implicated by the temporal distributions of the largest larvae. Ocean-atmosphere changes have been hypothesized to affect the survival of the larvae and cause reduced recruitment, so even after about a century following the discovery of their spawning areas, mysteries still remain about the marine life histories of the Atlantic eels.

Genome-wide single-generation signatures of local selection in the panmictic European eel

Next-generation sequencing and the collection of genome-wide data allow identifying adaptive variation and footprints of directional selection. Using a large SNP data set from 259 RAD-sequenced European eel individuals (glass eels) from eight locations between 34 and 64(o) N, we examined the patterns of genome-wide genetic diversity across locations. We tested for local selection by searching for increased population differentiation using F(ST) -based outlier tests and by testing for significant associations between allele frequencies and environmental variables. The overall low genetic differentiation found (F(ST) = 0.0007) indicates that most of the genome is homogenized by gene flow, providing further evidence for genomic panmixia in the European eel. The lack of genetic substructuring was consistent at both nuclear and mitochondrial SNPs. Using an extensive number of diagnostic SNPs, results showed a low occurrence of hybrids between European and American eel, mainly limited to Iceland (5.9%), although individuals with signatures of introgression several generations back in time were found in mainland Europe. Despite panmixia, a small set of SNPs showed high genetic differentiation consistent with single-generation signatures of spatially varying selection acting on glass eels. After screening 50 354 SNPs, a total of 754 potentially locally selected SNPs were identified. Candidate genes for local selection constituted a wide array of functions, including calcium signalling, neuroactive ligand-receptor interaction and circadian rhythm. Remarkably, one of the candidate genes identified is PERIOD, possibly related to differences in local photoperiod associated with the >30° difference in latitude between locations. Genes under selection were spread across the genome, and there were no large regions of increased differentiation as expected when selection occurs within just a single generation due to panmixia. This supports the conclusion that most of the genome is homogenized by gene flow that removes any effects of diversifying selection from each new generation.