- Research article
- Open Access
Evidence of local adaptation in a waterfall-climbing Hawaiian goby fish derived from coupled biophysical modeling of larval dispersal and post-settlement selection
BMC Evolutionary Biologyvolume 19, Article number: 88 (2019)
Local adaptation of marine and diadromous species is thought to be a product of larval dispersal, settlement mortality, and differential reproductive success, particularly in heterogeneous post-settlement habitats. We evaluated this premise with an oceanographic passive larval dispersal model coupled with individual-based models of post-settlement selection and reproduction to infer conditions that underlie local adaptation in Sicyopterus stimpsoni, an amphidromous Hawaiian goby known for its ability to climb waterfalls.
Our model results demonstrated that larval dispersal is spatio-temporally asymmetric, with more larvae dispersed from the southeast (the Big Island) to northwest (Kaua‘i) along the archipelago, reflecting prevailing conditions such as El Niño/La Niña oscillations. Yet connectivity is nonetheless sufficient to result in homogenous populations across the archipelago. We also found, however, that ontogenetic shifts in habitat can give rise to adaptive morphological divergence when the strength of predation-driven post-settlement selection crosses a critical threshold. Notably, our simulations showed that larval dispersal is not the only factor determining the likelihood of morphological divergence. We found adaptive potential and evolutionary trajectories of S. stimpsoni were greater on islands with stronger environmental gradients and greater variance in larval cohort morphology due to fluctuating immigration.
Contrary to expectation, these findings indicate that immigration can act in concert with selection to favor local adaptation and divergence in species with marine larval dispersal. Further development of model simulations, parameterized to reflect additional empirical estimates of abiotic and biotic factors, will help advance our understanding of the proximate and ultimate mechanisms driving adaptive evolution, population resilience, and speciation in marine-associated species.
The balance of migration and selection can determine the likelihood of local adaptation and evolutionary divergence within species [1,2,3,4]. For instance, migration can result in populations failing to reach local fitness optima for adaptive traits, even when natural selection favors local adaptation, such as when species occupy heterogeneous habitats [5, 6]. The theory underlying outcomes like this has been well studied, often with deterministic models of genetic architecture [7,8,9,10]. However, few natural systems have been examined in which modeled parameters are well constrained by empirical estimates of gene flow and the strength and direction of natural selection [11, 12]. There is a particular deficit in understanding the conditions that favor evolution of local adaptation in oceanic and migratory species with marine dispersing larvae [13,14,15,16,17]. Most prior work has focused on estimates of larval dispersal [1, 18,19,20] rather than integrative analysis of connectivity, post-settlement selection, and survival to reproduction [15, 21, but see 22].
Understanding how interactions between migration and selection influence the evolution of species with marine-dispersing larvae is challenging because the processes that govern population connectivity have not been well quantified. For example, a number of oceanographic transport models [23, 24] suggest that species with a longer pelagic larval duration (PLD) should have more “open” populations, whereas species with a shorter PLD should have more “closed” populations [25, 26]. Yet PLDs are often unknown. Moreover, recent studies using advanced modeling methods have not supported this previously prevailing theory. For example, studies of species with well constrained PLDs (e.g., via otolith microchemistry analyses) indicate that transport models tend to overestimate dispersal distance [27,28,29,30]. Recent comparisons of dispersal models with empirical larval drift surveys also indicate that PLDs above a critical value do not influence population connectivity . In addition to oceanographic features (e.g., eddies, currents, tides), organismal attributes (e.g., larval swimming behavior, vertical migration) can impede dispersal and thus govern connectivity [19, 32,33,34,35,36]. Incongruent estimates of connectivity may also be due to post-settlement selection, which can mediate demographic contributions to larval pools (i.e., individuals that reach settlement sites must survive to reproduction in order to contribute to larval pools). Thus post-settlement selection and differential reproduction arising from habitat heterogeneity can potentially shape connectivity [22, 37,38,39,40] and influence evolutionary trajectories of species with marine larval dispersal .
Spatio-temporal variability is a well-recognized underlying aspect of population connectivity via marine larval dispersal. Not only can larval dispersal vary over space and time [41,42,43,44], so can selective pressures [45,46,47]. For instance, selection can vary ontogenetically if life history stages occur in different habitats (e.g., oceanic versus stream environments). And, though most analyses of connectivity assume that landscape features are deterministically static over ecological timescales [48, 49], shifts in physiographic features can generate further changes in larval dispersal and post-settlement selection [50,51,52,53]. This is well illustrated by predictions that connectivity may be altered as a consequence of climate-driven shifts in ocean currents [54,55,56]. Few attempts have so far been made, however, to determine whether spatio-temporal variability promulgates or constrains the likelihood of local adaptation in species with marine dispersing larvae [57,58,59].
Empirically constrained, coupled modeling of larval transport and post-settlement selection can be an informative approach for assessing the likelihood of local adaptation in species with marine larval dispersal. Advection-diffusion (AD) oceanographic transport models, which allow for directional or stochastic spatial-temporal variation in larval dispersal [19, 60], generate connectivity matrices that can inform spatially-explicit, individual-based models (IBMs) of post-settlement selection and reproduction [61,62,63]. Reciprocally, by following the fate of individuals and entire populations over time [64, 65], IBMs can inform larval transport models by describing variation in larval pool contributions arising from post-settlement selection. Coupled AD-IBM modeling can thus be particularly useful for tracking outcomes of migration and selection across life history stages that occur in different habitats or ecosystems.
Here we present the results of a spatially explicit AD-IBM model developed to examine how migration and selection shape population connectivity and local adaptation in Sicyopterus stimpsoni, an amphidromous fish endemic to the Hawaiian Islands known for its ability to climb waterfalls. We linked a modified Lagrangian transport model of inshore and offshore oceanographic processes for the Hawaiian Islands to IBMs for the islands of Hawai‘i (hereafter referred to as the Big Island), O‘ahu, and Kaua‘i. We utilized the AD-IBM model to assess whether natural selection is sufficient to yield morphological divergence between subpopulations that are connected via marine larval dispersal. The model was run to test the following scenarios: (1) in a closed system (i.e., local reproduction, no immigration), stream topography and discharge yield morphologically divergent subpopulation among the islands; (2) the rate and direction of morphological divergence will change by altering the strength of post-settlement selection of predation evasion and/or climbing performance; (3) alternatively, subpopulations across the archipelago are morphologically homogenous in an open system with no local reproduction (i.e., immigration only); and (4) natural selection will counterbalance high, yet stochastic larval dispersal due to variation in the strength of post-settlement selection or immigration. We then compared our modeling results to empirical estimates of the strength of post-settlement selection [66,67,68] and observed morphological differentiation  to determine the relative influences of larval dispersal and selection shaping empirical patterns of morphological divergence among populations of S. stimpsoni juveniles and adults.
Oceanographic simulations of passive larval dispersal
Passive larval dispersal was not spatially or temporally uniform (Figure 1). In strong El Niño years (mid 2009-mid 2010), dispersal was asymmetric with more larvae dispersed from the southeast to the northwest (i.e., from the Big Island to Kaua‘i). This asymmetry was less pronounced in strong La Niña years (mid 2010-mid 2012). The extent of asymmetry fluctuated during neutral years (mid 2012-early 2014) of the El Niño Southern Oscillation (ENSO) cycle, with an onset of asymmetry (2012-2013) followed by symmetric dispersal (2013-2014), which might have been driven by a brief period of El Niño conditions in mid 2012 (U.S. Department of Commerce, National Oceanographic and Atmospheric Administration, NOAA Research: https://www.esrl.noaa.gov/psd/enso/climaterisks/years/).
The proportion of successful larval settlement and source contributions varied among islands. The Big Island had the highest proportion (2.92%) of successful settlement, whereas Kaua‘i had the lowest proportion (1.4%) compared to all other islands (1.4-2.1%). Despite very low probabilities of successful larval transport and fluctuating patterns of passive larval dispersal, the highest percentage (42%) of local entrainment occurred on the Big Island. The Big Island also served as a source of larvae for all other islands, though with diminishing contributions up the southeast-northwest axis of the archipelago. The second highest level of local entrainment (24%) occurred on Kaua‘i, but unlike the Big Island, Kaua‘i contributed relatively few larvae to other islands. The lowest level of local entrainment (18%) occurred on O‘ahu, which also contributed relatively few larvae to other islands. Slightly higher self-recruitment occurred on Maui and Moloka‘i (23% and 21%, respectively), which contributed a similar number of larvae to other islands (16-23%).
Linking larval dispersal with post-settlement selection
Scenario 1: Isolation without post-settlement selection
Differences in topographic structure and corresponding differences in stream flow did not give rise to divergent morphotypes across the archipelago. Rather, climbing morphotypes evolved and remained prevalent across all islands (Figure 2).
Scenario 2: Isolation with post-settlement selection
When post-settlement selection was allowed to act in conjunction with isolation, morphological divergence arose across islands in accordance with predictions: predation evasion morphologies arose on Kaua‘i, whereas climbing morphologies arose on the Big Island. While morphological divergence was driven by the strength of selection from predation, the amount necessary for the evolution of predation evasion morphotypes varied across islands (Figure 3). Adult morphology does not diverge from larval morphology with increasing selection probabilities, as would be expected in isolation (i.e., larvae are products of only local reproduction and adults are products of that larval composition) (Figure 3 & Additional File 1: Figure S1).
The highest selection differentials occurred on Kaua‘i whereas the lowest occurred on the Big Island (Table 1), with a range of <0.0001 to 0.3909 across all islands. Selection differentials significantly differed between stages (Table 2). With the exception of O‘ahu, where the transition between Stage 9-10 exhibited the highest maximum selection coefficient, the largest maximum and average values of selection differentials on other islands occurred between the transition from Stage 3-4 (i.e., the post-settlement stage; Tables 1 & 2). Stage was the only significant predictor of mean selection differentials on the Big Island, whereas stage and predation probability were significant predictors of mean selection differentials on Kaua‘i. On O‘ahu, all variables were significant predictors of selection differentials.
The RDAs explained 65-89% of the total variance in morphotype evolution across all islands. Morphotype evolution was significantly correlated with year, predation, climbing, and the interaction between predation and climbing (Table 3). However, each parameter explained <1.52% of the variance in morphotype evolution with the exception of predation (Table 3), which explained 16.10% of the observed morphotype variance on Kaua‘i, 19.96% on O‘ahu, and 24.42% on the Big Island.
Scenario 3: Immigration without post-settlement selection
Under conditions of immigration without post-settlement selection, populations on all islands evolved morphotypes that are intermediate between climbing and predation morphologies (Figure 4). Additionally, adult morphotypes were reflective of larval morphotypes.
Scenario 4: Immigration with post-settlement selection
When post-settlement selection was allowed to act in conjunction with immigration, morphological divergence arose across islands in accordance with predictions: predation evasion morphologies arose on Kaua‘i, whereas climbing morphologies arose on the Big Island. Morphological divergence was driven by both the strength of selection from predation and its interaction with immigration (Tables 2 & 3). The amount of predation selection necessary for the evolution of predation evasion morphotypes varied across islands (Figure 5).
Adult morphology diverged from larval morphology with increasing predation probability relative to immigration, as would be expected if post-settlement selection overcomes gene flow from immigration (Figure 5, Additional File 2: Figure S2, & Additional File 3: Figure S3).
The highest selection differentials occurred on Kaua‘i whereas the lowest occurred on the Big Island (Table 1), ranging from <0.0001 to 0.56. Selection differentials differed significantly between stages across all islands (Table 2). And with the exception of O‘ahu, where the transition between Stage 9-10 exhibited the highest maximum selection coefficient, the largest maximum and average values of selection differentials occurred during the transition between Stage 3-4 (i.e., the post-settlement stage; Table 1). Predation probability was the only other predictor of mean selection differentials on O‘ahu and the Big Island. It was not a predictor of mean selection differentials on Kaua‘i, where mean selection differentials were significantly correlated with immigration rate and the interaction between immigration rate and predation probability (Table 2).
The RDAs explained 78-94% of the total observed variance in morphotype evolution across all islands. Morphotype evolution was correlated with year, immigration, predation, and the interaction between immigration and predation (Table 3). Predation explained the highest proportion of variance in morphotype for each island (25.93% - Kaua‘i, 43.23% - O‘ahu, and 33.28% - Big Island; Table 3). The interaction between immigration and predation explained 11.50% of observed variance on O‘ahu, 6.57% on the Big Island, and only 1.93% on Kaua‘i. Immigration and year explained less than 3% and 0.01%, respectively, of morphotype variance on each island.
Coupled AD-IBM modeling demonstrated that the likelihood of local adaptation in a diadromous fish depends on the spatio-temporal balance of pelagic larval dispersal and post-settlement selection . Our results affirm theoretical expectations that passive larval dispersal facilitated by ocean currents can result in homogenization [23, 60, 69], but we also found that dispersal asymmetries can arise due to climate-driven (i.e., ENSO) fluctuations in oceanic conditions. Our model simulations additionally showed that spatio-temporal variation in post-settlement selective pressures can override the homogenizing effects of passive larval dispersal, where differential probabilities of survival and reproduction result in adaptive evolution and morphological divergence among populations. But contrary to the expectation that post-settlement selection pressures from both predation and climbing drive population divergence, we found that predation alone is likely the primary driver of population divergence in S. stimpsoni. Notably, we also found that the amount of larval immigration was not a strong determining factor of morphological evolution. Rather, the strength of predation-driven, post-settlement selection and its interaction with immigration appear to shape morphological divergence across the Hawaiian archipelago, which is consistent with Fisher’s Theorem of Natural Selection [12, 70]. In other words, the strength of selection and rate of change in the mean morphotype increase as trait variance increases. In our simulations, this is driven by the continued influx of maladaptive morphotypes through immigration.
Passive larval dispersal
The complex oceanic hydrodynamics surrounding the Hawaiian Islands makes it difficult to generalize about the dispersal of pelagic propagules across the archipelago . Nonetheless, it has long been hypothesized that mesoscale eddies that form on the leeward side of the Hawaiian Islands prevent larvae from being swept away, and thus promote retention and settlement near islands [72,73,74]. Contrary to this idea, recent studies have found little or no association between mesoscale eddies and local larval retention in the Hawaiian Islands [75, 76], but instead facilitated larval transport to other islands . Although, our AD model simulations do not directly test for the independent influences of these factors, they do indicate that the combination of mesoscale eddies and prevailing current flow patterns (southeast to northwest) contribute to local retention and larval transport among islands [34, 78, 79]. While average estimates of larval settlement were low across the Archipelago (< 3%), we found evidence of local larval retention around each island with the highest proportions of local larval retention occurring around the Big Island (42%). We also found asymmetric larval dispersal from the southeast (i.e., the Big Island) to the northwest (i.e., Kaua‘i). These results suggest that mesoscale eddies could contribute to increased larval transport among islands, which is consistent with the findings of Wren et al. (2016) . However, future work should directly quantify the independent effects of mesoscale eddies from prevailing current flow patterns to understand their respective influences on S. stimpsoni larval dispersal.
While it is evident that larval dispersal and spatial coagulation of larval cohorts  vary because of dynamic nearshore and offshore oceanographic conditions [81,82,83,84,85], our AD modeling also illustrates the importance of considering large-scale, climatic variation in assessments of population connectivity. ENSO and other large-scale climate-driven phenomena can exhibit small temporal-scale variability that enhance or suppress larval dispersal [85,86,87,88,89]. Consistent with evidence of ENSO influencing population connectivity of ocean-dwelling species across the Hawaiian archipelago and elsewhere in the Pacific [19, 90], we found that ENSO likely influences S. stimpsoni larval dispersal (Figure 1). Parallel patterns of temporal and spatial fluctuations in FST estimates from microsatellite markers among subpopulations of S. stimpsoni  serve as further evidence that larval cohort aggregation and composition shift in response to ENSO cycles.
Our AD model was conservative with respect to larval behavior and life history (e.g., the model did not consider larval mortality, vertical migration, or swimming behavior). Incorporating other parameters into the model- including life history variation , larval behavior , and selection ‘at sea’ (i.e., selection against long distance dispersal, which could increase the probability of larval loss)- would likely improve estimates of larval transport, local retention and self-recruitment (i.e., natal homing ). In combination with persistent natural and anthropogenic perturbations in oceanic island streams [92,93,94,95,96,97], theory would predict that such conditions would result in the evolution of dispersal polymorphisms  to protect against extirpation and extinction. Accounting for the possibility of alternative dispersal strategies [99,100,101], which have been found in other diadromous species [102,103,104,105,106,107,108], could shift the balance between larval transport and self-recruitment estimates. Indeed, self-recruitment could be many times greater than predicted by our AD model [1, 109,110,111]. Accordingly, the modeling framework we have so far developed offers a basis for further investigations including, but not limited to, modeling outcomes of alternative dispersal strategies to better understand the implications of life history variation in amphidromous gobies like S. stimpsoni [102, 112, 113].
Settlement site topography
We expected that differences in stream topography (e.g., slope, waterfall locations, and discharge) would result in divergent morphotypes across the archipelago, where steep-sloped streams with fast flows would harbor fish with long, shallow bodies, while shallow-sloped streams with slower flowing water have fish with short, deep bodies. To the contrary, under conditions of no immigration and selection probabilities turned off, we found that streamlined “climbing” morphotypes evolved in all populations, which were initially seeded with a homogenous distribution of morphotypes, regardless of prevailing topography (Figure 2). This finding suggests that, upon initial colonization (likely during the Pleistocene ), S. stimpsoni exhibited a streamlined body shape that better enabled it to penetrate steep-sloped watersheds across the archipelago. At that time, streams on the fully-formed main islands (Maui, Moloka‘i, O‘ahu, and Kaua‘i) were most likely steeper in slope and perhaps more homogenous in topography compared to current conditions. Thus our model results offer support for the idea that initial colonization of the archipelago by amphidromous fauna reflects temporally dynamic, yet time specific opportunities, where habitat suitability and availability shift with island age [114, 115].
Our findings illustrate that strong selection from post-settlement mortality during juvenile recruitment can promote divergence [13, 116, 117] because watersheds across the Hawaiian Islands are heterogeneous. Though environmental gradients on each island are similar in scale, there are stark contrasts between islands across the archipelago - heavily eroded watersheds with high sedimentary loads and waterfalls far inland on older islands, and steep watersheds produced by recent volcanic flows with fast moving clear water and waterfalls close to shore on younger islands. Thus larvae recruiting from non-local sources may encounter highly unsuitable habitat. Prior work has shown, for example, that newly recruiting juveniles are exposed to strong directional selection favoring either climbing performance on the Big Island or predation performance on Kaua‘i [66,67,68, 118]. Our findings also are consistent with inferences that variation in non-linear selection between islands and watersheds can result in complex fitness surfaces that allow for the evolution of locally adapted populations .
IBM simulations intended to assess outcomes of post-settlement selection without immigration (scenario 2) indicated that selection is strongest during recruitment (ontogenetic Stage 3-4). This is consistent with expectations that post-settlement selection ought to be strongest during this transition stage because stream topographic structures (i.e., the proportion of reach types and presence of waterfalls) dictate that selection from predation and climbing occurs soon after stream re-entry. It is worth noting, however, that the influence of associated model parameters (e.g., predation and climbing probabilities, and the interaction of predation and climbing probabilities) on selection differentials varied for each island (Table 2).
Similarly, IBM simulations intended to assess outcomes of post-settlement selection with immigration (scenario 4) indicated that selection is strongest during recruitment (ontogenetic Stage 3-4). But, in contrast to simulations without immigration (scenario 2), allowing for immigration resulted in selection differentials being strongly correlated with either predation probability (O‘ahu and the Big Island) or interactions between immigration and predation probability (Kaua‘i). Notably, because our modeled selection differentials are congruent with empirical estimates of predation selection on S. stimpsoni [66, 68], this finding supports prior interpretations and inferences about the evolution of morphological divergence among populations of S. stimpsoni [53, 66, 68].
Our IBM selection differentials and morphotype results were pooled averages within an island, whereas our empirical estimates of the strength of selection and morphology were measured on a stream basis [66,67,68, 118]. Despite these differences in geographic scale, our simulation results and conclusions are comparable to our empirical findings. In natural populations, the strength of predation selection differentials on fitness-related traits (e.g., mid-body) ranged between 0.021-0.28 [66, 68]. In our simulations, once the mean of the selection differential exceeded 0.03, predation morphologies could evolve on Kaua‘i regardless of immigration rate. However, the degree of predation morphotype (i.e., 6, 7, 8, or 9) did vary with immigration rate (Figure 5). In contrast, even when the selection differentials reached or exceeded this threshold value on the Big Island, predation morphologies never evolved even at low levels of immigration, which is consistent with natural populations of S. stimpsoni residing in Big Island streams .
Theoretical and empirical studies have demonstrated that predators can engender adaptation and persistence of prey in a fluctuating environment [12, 119,120,121]. This may occur through a selective push, in which selection moves the mean of a trait towards a local optimum by predators consuming prey with lower fitness due to the trait in question. It may also occur as a result of the ‘evolutionary hydra effect’, in which predation reduces prey density and then increases prey birthrate, resulting in more selective events per unit time, which effectively reduces generation time . Our results indicate that with post-settlement selection, regardless of immigration (i.e., whether it is included in the model or not, or whether the rate varied), distinct morphologies evolve on each island in accordance with prior predictions (i.e., predation evasion morphologies on Kaua‘i, climbing morphologies on the Big Island). The recovered patterns of divergence were principally driven by the strength of predation-derived selection (Table 3), which is consistent with the pattern of predation-driven divergence found across a broad range of taxa [122,123,124,125,126,127,128,129]. However, the strength of selection required for the evolution of predation evasion morphotypes varied among islands. Much weaker predation-driven selection was necessary on Kaua‘i compared to the Big Island or O‘ahu (Figure 3). This is well illustrated in model runs with immigration, which showed that predation evasion morphotypes do not evolve on the Big Island even when predation-driven selection is strong, whereas predation evasion morphotypes could evolve on Kaua‘i regardless of the immigration rate (Figure 5). This is highlighted by the finding that the interaction of predation-driven selection and immigration, not the strength of predation-driven selection alone, was the driving factor for adaptive evolution on Kaua‘i (Tables 2 & 3). These results suggest that the adaptive potential and adaptive evolutionary trajectory of S. stimpsoni may be greater on islands that have strong environmental gradients and that receive recruits with greater variance in morphology due to immigration (Figure 5).
Contrary to expectations derived from observations of natural populations of S. stimpsoni on the Big Island and Kaua‘i , estimates of the opportunity for selection were greater for the non-primary pressure on each island (i.e., climbing on Kaua‘i, predator evasion on the Big Island) than the primary pressures. Canonical rotations of the nonlinear gamma matrix (i.e., a canonical transformation that identifies the major axes of fitness surfaces and facilitates the detection of stabilizing or disruptive selection) demonstrated that individuals from Kaua‘i and the Big Island occupy regions near their local fitness peaks for some traits . Therefore, selection for predator evasion on Kaua‘i and climbing on the Big Island may be less effective in promoting morphological change in S. stimpsoni than the non-primary pressures because variation of functionally important traits may have been reduced by directional or stabilizing selection. Our model simulations do not recover this relationship for two reasons. First, the modeled morphotypes are on a linear axis of 1-10, with morphotype being coded as a univariate trait. Consequently, the variance for some fitness-related traits could not be reduced by selection and non-selected traits maintain a high level of variance. This increases the potential for greater non-primary selective pressures to operate on each island (i.e., climbing on Kaua‘i and predation on the Big Island) resulting in greater opportunities for selection in those directions; however, this possibility could not be explored in our models. Second, this inconsistency may instead be reflective of much lower connectivity in nature than was accounted for in our oceanographic models, possibly due to a range of factors (e.g., differential larval mortality, larval swimming behavior, vertical migration, or local retention), which would be expected to greatly reduce variance in larval cohort morphology through cohort coagulation [52, 53]. Resolving this incongruence thus warrants continued exploration, parameterization, and quantification of spatio-temporally fluctuating processes that contribute to population connectivity in marine and diadromous species with marine dispersing larvae.
Connectivity and conservation
Further work on population connectivity also stands to advance understanding of population resilience in an ever changing, human-dominated environment. This is particularly pertinent to species endemic to oceanic islands, like the Hawaiian Islands, where escalating growth of the resident human population is imposing ever-greater pressures (e.g., via habitat alteration, pollution, resource divergence) and threats (e.g., increasing temperatures, climate-driven shifts in oceanic currents, invasive species, etc.) [56, 130,131,132]. Many of these concerns are evident across the Hawaiian archipelago, especially on the island of O‘ahu. Home to over 80% of the human population in the Hawaiian Islands, O‘ahu has undergone extensive urbanization over the past century [133, 134]. Populations of S. stimpsoni on O‘ahu are either locally extirpated (i.e., within particular watersheds) or significantly smaller than those found on other islands , which reflects a sensitivity to perturbation and habitat modification [97, 135]. While our model results indicate that S. stimpsoni larvae can recruit and that local adaptation can evolve on O‘ahu, prevailing conditions appear to be overwhelming both processes on the island. It is possible, however, that the extirpation and decline of S. stimpsoni on O‘ahu is a consequence of significantly lower larval dispersal than what we have modeled, exacerbated by small adult populations producing relatively few larvae. It is also possible that larval dispersal is sufficient to ‘seed’ watersheds across O‘ahu, but that in-stream conditions are so unsuitable for S. stimpsoni that post-settlement mortality is too high for local populations to persist. Under either set of conditions, resident populations likely cannot be “rescued” by immigration. Further modeling that accounts for additional abiotic and biotic parameters  could improve understanding of why S. stimpsoni (and other endemic amphidromous species) are nearly extirpated on O‘ahu and thus help identify strategies to best reverse losses and prevent losses on other islands in the archipelago.
Unlike the terrestrial systems in the Hawaiian Islands, which exemplify adaptive radiations [137,138,139], the endemic amphidromous fishes, as well as the marine biodiversity of Hawai‘i, are not the products of evolutionary radiations [140,141,142,143], but see . Even though the Hawaiian gobies reside within the islands like honeycreepers, Drosophila, and silverswords, these gobies are constrained by their phylogenetic history, which includes a life cycle that dictates oceanic larval development resulting in high gene flow amongst populations [53, 114, 135, 145]. It is this contrast between adaptive radiations amongst terrestrial species but not marine species, and the underlying evolutionary mechanisms of adaptation and speciation, for which amphidromous gobies provide distinct insight. Our study demonstrates that despite strong local adaptation resulting from the complex interplay between ocean currents, dispersal strategies, and post-settlement selection, the opportunity for ecological speciation amongst marine-associated organisms is still most likely constrained as a result of gene flow. However, with recent advances in technological and bioinformatic methodologies, phylogeographic and population genetic studies of anadromous and marine species are increasingly finding evidence that barriers to gene flow exist at fine spatial and temporal scales [22, 52, 146,147,148,149,150,151,152,153,154,155]. In combination with further development of oceanographic model simulations, especially ones parameterized to reflect empirical estimates of abiotic and biotic factors (e.g., dispersal potential, larval mortality and swimming behavior, as well as post-settlement selection and ecological conditions), such studies can advance our understanding of adaptive evolution, speciation, and population resilience in an ever-changing aquatic environment.
The amphidromous Hawaiian goby fish, Sicyopterus stimpsoni, is well suited for studying how migration and selection influence the evolution of species with marine-dispersing larvae. While adult subpopulations of S. stimpsoni predominantly reside in upper elevation stream habitat, dispersal among streams occurs via oceanic transport of pelagic larvae [156, 157]. Dispersal distances and among-watershed connectivity are not well known, however, as it is unclear how inshore and offshore oceanic processes affect the spatio-temporal structure of larval dispersal. It is known that recruitment to streams results from pelagic larvae cueing on freshwater plumes , and evidence of life history variation in another goby endemic to the Hawaiian Islands  suggests that the probability of post-larvae returning to natal streams may reflect a combination of active (e.g., habitat use, larval swimming behaviors) and passive (e.g., oceanographic transport) factors. Thus, it is likely that cohorts that recruit back to streams and migrate upstream to adult spawning habitat are composed of individuals originating from a combination of sources (i.e., natal and other streams). Consistent with this, patterns of little to no neutral genetic differentiation in microsatellite and mitochondrial markers, among adult subpopulations across the Hawaiian archipelago, suggest that recruiting cohorts draw from well-mixed larval pools [53, 112, 114, 135, 145, 159, 160].
During upstream migration, both predation and waterfall climbing produce brief, but intensive, episodes of selection on juvenile S. stimpsoni [66,67,68]. Juvenile recruits must escape predation (e.g., by the endemic, piscivorous sleeper, Eleotris sandwicensis), and individuals who survive predation face the additional selective pressure of climbing waterfalls before reaching predator-free adult habitats where S. stimpsoni grow, mature, and reproduce [161, 162]. Both forms of selection may lead to local adaptation in morphology among streams if additive genetic variation exists for shape differences [53, 66,67,68].
The strength of natural selection from predation and waterfall climbing varies according to watershed topography. The geomorphology of Hawaiian watersheds spans a topographic gradient that tracks the progression of erosion with island age . Variation in topography translates to differences in the steepness of stream slope, stream depth and breadth, and surface flow rates [53, 92]. For example, streams on the Big Island are characteristically steep-sloped, and often terminate in waterfalls with little to no estuarine habitat. Individuals recruiting to “Big Island-like” streams must climb waterfalls within a few days of settlement. Therefore, natural selection should favor individuals with streamlined morphologies (i.e., long, shallow bodies) for reduced drag and better waterfall-climbing performance [164, 165]. On the other hand, streams on Kaua‘i are characteristically shallow-sloping with waterfalls that are often located kilometers inland. Individuals recruiting to “Kaua‘i-like” streams thus may have to swim upstream for weeks before escaping predation by climbing waterfalls. Therefore, natural selection should favor individuals with short, deep bodies that facilitate greater thrust production for predator evasion [164,165,166]. Accordingly, differences in stream topography likely constitute a mosaic of selection pressures that underlies observed patterns of divergence in body shape among sub-populations [53, 66,67,68] under conditions of high gene flow.
Modeling passive larval dispersal: oceanographic simulations
Patterns of larval dispersal can change in both space and time [41,42,43,44], therefore we used a modified version of the 2D Lagrangian transport model described in Polovina et al. (1999)  and Wren & Kobayashi (2016)  to simulate larval dispersal and recruitment. As implemented here, the transport model tracks individual larvae (i.e., particles) through physical oceanographic models from spawning to settlement, moving them through ocean velocity fields. We set the diffusivity coefficient to 250 m2/sec (, Jia pers. comm.) and we used a settlement radius of 5 km around each stream mouth. For physical flow fields we used a regional implementation of the daily Hybrid Coordinate Ocean Model (HYCOM) , with a K-profile parameterization (KPP) mixed layer formulation for our current solutions (http://apdrc.soest.hawaii.edu/datadoc/hycom_hawaii_0.04_kpp.php). The domain of the regional HYCOM model spans 194°E-210°E and 16°N-26°N and has a 1/25° horizontal resolution and 28 active vertical layers. This model is eddy resolving, which accurately predicts mesoscale eddies that are often present in the lee of the Hawaiian Islands . For simulation runs, we depth-averaged current velocities in the top 100 m (i.e., the top 7 layers of the HYCOM), excluding the surface layer (0-5 m). We used this approach because it has been shown to best predict the position and settlement of larval reef fishes from empirical studies within the Hawaiian Islands  and there is no evidence from plankton and midwater trawl surveys that larval goby distributions are neustonic or epipelagic (i.e., under the sole influence of the HYCOM surface layer) .
A total of 51 stream mouth locations on Kaua‘i, O‘ahu, Moloka‘i, Maui, and the Big Island were used as release and settlement sites in the model (Figure 6). One hundred virtual larvae were released daily from each stream mouth location from 2 May 2009 until 31 March 2014 (1594 days), totaling over 8 million (8,129,400) virtual larvae released during the simulation. We conducted separate preliminary AD model simulations with PLD set to 55 and 150 days post-release and found slightly but not significantly higher levels of connectivity at 150 PLD compared to 55 PLD. Further, the patterns we observed when running these preliminary simulations were very similar to the overall patterns of connectivity found when using a continuous range of PLDs. Therefore, in an effort to best represent the entire range of possible PLDs of S. stimpsoni, we used a continuous bimodal range of 50 days to 200 days based on estimates of age and PLD from otolith microchemistry studies [102, 172,173,174]. Each larva could settle as early as 50 days after release and up to 200 days after release, after which mortality occurred. In all AD simulations, larvae were assumed to be passive throughout their PLD. Simulation results were averaged across three runs (to account for stochasticity) with daily time steps.
Linking larval dispersal with post-settlement selection: IBM development
Three independent, spatially explicit, individual-based models (IBMs) were developed and parameterized in NetLogo (6.0.3)  for Kaua‘i, O‘ahu, and the Big Island (respectively). We limited our IBMs to these three islands because our prior work quantified selection and morphological divergence on Kaua‘i and the Big Island [53, 66,67,68] allowing for direct comparisons of simulation results with empirical data. In addition, populations of S. stimpsoni on O‘ahu are extremely rare and perhaps locally extirpated, making this island of interest for conservation concerns. Each IBM was divided into habitat cells of 1 km2 and scaled to island size. The habitat cell types represented the primary habitats utilized by S. stimpsoni: offshore ocean and nearshore ocean, which are larvae-only cells, and estuarine and upstream, which are juvenile- and adult-only cells, according to the species life history [81, 176].
We selected eight to ten watersheds on each island representative of habitat type and elevation to include in the island-specific IBMs. Physical characteristics (length, elevation, slope) for each of the selected watersheds were determined from the Atlas of Hawaiian Watersheds and their Aquatic Resources (http://www.hawaiiwatershedatlas.com/index.html). Additional data gathered from Google Earth enabled us to develop metrics for waterfalls in the IBMs. The highest elevation areas within watersheds were excluded from the IBMs because S. stimpsoni do not inhabit headwater reaches [176, 177].
Patch conditions were watershed-specific with weekly discharge estimates based on historical records from United States Geological Survey stream gauges (http://hi.water.usgs.gov). Weekly discharge rates (ft3/sec) were averaged across all years, and an annual sine curve model was calculated for each watershed (Additional File 4: Supplemental Equations). For watersheds without gauge data (22 streams), we used log(stream channel length km), log(watershed area km2), and log(max elevation m) in a principal component analysis to calculate the nearest neighbor distance. Nearest neighbor distance (√((x2-x1)2+(y2-y1)2+(z2-z1)2)) was determined using principal component axes 1-3, and used to assign annual sine curve models of discharge rates to watershed without gauges.
Each IBM followed the survival of individual S. stimpsoni through ten life history stages, from late-stage pelagic larvae to adults on a weekly time step basis. Stage 1 corresponded to the final marine larval stage, with corresponding individuals having reached the nearshore waters surrounding the island but not yet having detected freshwater plumes. Thus, Stage 1 movement is randomly oriented in the offshore marine environment with respect to stream locations. Stage 2 represented the nearshore post-larval stage, during which individuals exhibit strong orientation toward freshwater plumes depending on the width and flow from each watershed. Stage 3 constituted juvenile stream fish with fully developed benthic feeding and climbing structures that orient upstream and continue to migrate until reaching suitable adult habitat and conspecific density. Stages 4 through 10 consisted of reproductively mature adult fishes that retain upstream orientation until reaching suitable adult habitat and conspecific density . Each stage differentiated by growth in approximately 1 cm increments , which influences the probabilities of predation, growth and net reproductive output.
Each individual agent (fish) in the model possessed the following characteristics: stage (1-10), sex (male or female), position (x-y coordinates), direction (degrees), and morphotype (0-1). Morphotype consisted of a single trait on a scale from 0 to 1 that describes a continuum of body shape from long, shallow bodies (climbers) to short, deep bodies (predator evaders). Because S. stimpsoni morphotype heritability is unknown, we treated trait heritability as the average of the two parental morphotypes +/- 5% of the average difference between parental morphotypes, which falls within the range of narrow-sense heritability estimates for morphological traits in fishes [179, 180] and is calculated as:
Each IBM was populated with initial agents of each stage with an equal distribution of all morphotypes. Above a threshold that could cause random extinction, the initial number of agents had little effect on model output, and was therefore set at 3000 gobies to minimize that probability. Upon initiation of the model run, the IBM proceeds as a continuous loop of six subroutines executed on a weekly time step: patch conditions, movement, mortality, reproduction, growth, and immigration. Larvae (Stage 1) begin in the offshore zone of the model and undergo mortality if natural_mortality calculated as:
is greater than a random probability (between 0 and < 1 assigned at each time step), or move into the nearshore zone if they survive to Stage 2. Stage 2 larvae must then find a watershed inlet habitat cell in order to move into a stream or else perish. Upon entering the watershed (Stages 3-10), the probability of moving upstream depends upon the elevation gradient between the current patch and the next patch, the current hydrological discharge conditions of the patch, and the morphotype of the individual. Thus, the probability of upstream movement was calculated as:
where discharge is the logarithmic stream-specific sine curve equations (Supplemental Materials Eq. 1.1) scaled by a factor of 1/3 that converts discharge to a unit of climbing difficulty between 0-1. The parameter depthdiff, a proxy for waterfall height, is the absolute value of the elevation difference between the current patch and the next patch. depthdiff is scaled by a factor of 1/50 that converts this parameter to a unit of climbing difficulty, resulting in a similar distribution of climbing success as observed by Blob et al. 2008, 2010 [67, 68] and a longitudinal distribution of individuals within a stream similar to that observed in nature (Atlas of Hawaiian Watersheds and their Aquatic Resources, http://hawaiiwatershedatlas.com). The climb_threshold parameter is user-defined on a scale from 0 to 1, which is a primary parameter that is changed between our simulated scenarios.
Individuals that fail to climb to the next patch during the current time step remain in the same patch until the next time step or until they experience mortality. Mortality is stage-specific and can result from natural mortality (i.e., aging), predation, or competition. Natural mortality is calculated as in Equation 2. Predation mortality is calculated as:
where predation_risk is user-defined on a scale from 0 to 1 as primary parameter that is changed between our simulated scenarios. An individual experiences predation mortality when predation_mortality is greater than a random probability between 0 and < 1 assigned at each time step. Mortality due to competition is based on carrying-capacity and is calculated as:
where n_individuals_r1 is the number of individuals within a 1 cell radius of an individual and carrying_capacity is user-defined. We set carrying_capacity to 100 gobies per cell, because this value resulted in an asymptotic population growth curve, whereas values below this threshold resulted in population extinction, and above this threshold resulted in exponential population growth. If competition is greater than a random probability between 0 and < 10, then an individual experiences mortality from competition.
Male S. stimpsoni, which are territorial, not only have an elaborate courtship ritual that precedes pair-forming and spawning but also guard fertilized egg clutches [161, 176]. Therefore only individuals of reproductive size (Stages 4-10) were allowed to reproduce if an individual of the opposite sex was within one habitat cell. Reproduction, though year-round, is also linked to stream discharge, which is seasonally variable . Thus we calculated reproduction probability as:
resulting in a normal distribution of probabilities that are scaled to a unit of reproduction between 0 and 1. The number of offspring was determined by:
where births is a user-defined parameter ranging from 0-500, and reproduction is:
Newly-produced larvae were assigned a morphotype as described above.
In nature, both female and male S. stimpsoni have indeterminate growth . In the IBMs, growth rate was expressed as a logistic curve determined by stage dependent size (i.e., faster relative growth of younger, smaller individuals compared to older, larger individuals) and an annual sine curve of temperature (i.e., faster growth in warmer waters) calculated as:
and is scaled by a factor of 1/29 and added to a random number between 0 and < 5 for conversion to a unit of growth. Stage 10 individuals were determined to undergo senescence when the conditions for growth to the next stage were met, which occurred on average at ~165 weeks (i.e., 3.17 years after obtaining the size maximum).
Immigration rules in the IBMs consisted of weekly connectivity matrices from the AD model. Predicted pelagic larval morphotype distributions for these AD connectivity matrices were determined by the slope of the stream from which larvae were released, as stream slope gradient is correlated with morphotype in natural populations (i.e., individuals from steep sloping streams have long, shallow bodies, whereas, shallow sloping streams have individuals with shorter, deeper bodies ). The slopes of the 51 streams in the AD model were determined by regressing the distance from the mouth by the elevation of the site every 1000 meters upstream. The slope values were then used in a linear regression to predict a numerically scaled morphotype for individuals from each of the 51 streams. Predicted morphotype distributions were used to seed immigrant morphotypes (Additional File 5: Figure S4).
Linking larval dispersal with post-settlement selection: IBM simulation scenarios
Scenario 1: Isolation without post-settlement selection
With age, islands progressively erode and eventually subside into the ocean [163, 181]. Thus, topographic differences of Hawaiian watersheds should shape the morphological distributions of S. stimpsoni such that fish on older islands (i.e., Oʻahu and Kauaʻi) have shorter, deeper bodies (e.g., a predator evader morphotype) compared to fish on younger islands (i.e., Molokaʻi, Maui and the Big Island). We therefore ran each island IBM in isolation (i.e., local reproduction only) with climbing and predation selection parameters turned off, in order to assess whether island topographic differences alone give rise to morphological divergence.
Scenario 2: Isolation with post-settlement selection
Similar to patterns of larval dispersal, selective pressures can change over space and time [45,46,47]. Because in-stream selection pressures of predation and climbing are predicted to shape morphological variation, S. stimpsoni in the model should evolve morphotype distributions similar to those observed in nature, consistent with biomechanical predictions of optimal shapes for predator-evasion or waterfall-climbing performance [67, 68]. Accordingly, we ran each island IBM in isolation and with varying climbing and predation post-settlement selection probabilities (incrementally from 0 to 1) to assess whether and how, post-settlement selection gives rise to morphological divergence.
Scenario 3: Immigration without post-settlement selection
Immigration may erode local adaptation because larval recruits from different sources may arrive at settlement sites with suboptimal morphotype distributions . Natal retention might therefore be expected to increase survival and population persistence. Immigration also may be key to population persistence, however, because oceanic island streams are prone to natural and anthropogenic perturbation [92, 93]. We thus ran each island IBM with immigration only (i.e., no local reproduction) and without climbing and predation post-settlement selection to assess whether passive larval dispersal gives rise to morphological homo/heterogeneity and how immigration influences local demography.
Scenario 4: Immigration with post-settlement selection
Post-settlement selection may be strong enough to override the homogenizing effect of immigration via larval dispersal. We thus ran each island IBM with varying (0-1) climbing and predation post-settlement selection probabilities and varying levels of immigration, to assess the potential for morphological divergence of S. stimpsoni under conditions of migration and selection. The extent of immigration was set by weighting the connectivity matrices by 0.25, 0.5, 0.75, and 1, since the connectivity matrices reflect purely passive larval dispersal and thus represent a conservative estimate of dispersal across the archipelago. The weighting parameter therefore estimates the effect of reduced larval dispersal from biotic (i.e., active swimming behavior, vertical migration) or abiotic processes (i.e., temperature, salinity, chlorophyll) demonstrated in other marine and diadromous species  but not taken into account in our larval dispersal model. Because simulations of scenario 2 (see results) showed that climbing selection exerted negligible influence, we left climbing probability set to 1 and did not further analyze climbing for scenario 4. All IBM simulations were run 10 times for 200 generations.
From the oceanographic AD simulations, successful larval settlement was defined as any larvae occurring within a 5 km distance of a stream mouth. To capture inter-annual variation in larval transport success, a rearward probability connectivity matrix for each year was calculated as:
where Sij is the number of larvae release from stream i (source stream) that successfully settled at stream j (receiving stream), and Sj is the total number of larvae that successfully settled at a stream regardless of release stream. Each of the cells of the Pij matrices shows the probability of a particle transported to stream j having originated from stream i each year, averaged across the three simulation replicates, where each row in a matrix sums to 1. The diagonal of the probability matrices shows the amount of local entrainment for each stream, which is defined as the proportion of successfully transported larvae at each stream that originated from that same stream. We also calculated the overall proportion of local entrainment for each island across all years of the model by dividing the number of released larvae from all streams on an island that were successfully transported back to any stream on that same island.
To determine if our modeled selection parameters were within the range of empirical selection differentials [67, 68], we calculated the simulated selection differentials (s) for each island in tests of scenario 2 and scenario 4. Selection differentials were calculated as the difference between the mean morphotype of stage Jn+1 minus the mean morphotype of Jn within each time step of the IBMs . The modeled selection differential encompassed both the effects of climbing and predation selection. Consequently we used generalized linear models (GLMs) to quantify the degree to which age, predation, climbing or the interaction of predation and climbing contributed to our modeled selection differentials from IBM simulations of post-settlement without immigration (scenario 2). For IBM simulations of post-settlement selection with immigration (scenario 4), we used GLMs to quantify the degree to which age, immigration, predation, or the interaction of immigration and predation were contributing to our modeled selection differentials. Because simulations of scenario 2 showed that climbing selection exerted negligible influence on selection differentials (see results), we set the climbing probability in all IBM simulations of scenario 4 to a value of 1. Subsequently we did not use climbing as a parameter in the GLMs for scenario 4.
We used Redundancy Analysis models (RDA), a multiple linear regression ordination method , computed with the vegan package in R , to quantify the contribution of each model parameter to morphotype evolution under conditions of isolation and selection (scenario 2), as well as immigration and selection (scenario 4). We conducted independent RDAs for each island for isolation with selection, with predictor variables of year, predation, climbing, and the interaction of climbing and predation (scenario 2). For models of immigration and selection (scenario 4), we conducted independent RDAs for each island with the predictor variables of year, immigration, predation, and the interaction of immigration and predation. Because simulations of scenario 2 showed that climbing selection exerted negligible influence on morphological change (see results), we set climbing probability in all IBM simulations of scenario 4 to a value of 1. Subsequently we did not include climbing as a parameter in our RDAs for scenario 4.
We estimated the adjusted coefficient of determination (R2adj) for each model, with statistical significance determined using permutation tests to compare observed and randomized model R2adj. Since morphotypes at time t+1 are correlated with morphotypes at time t, we conducted variance partitioning with partial RDAs to estimate the variance in morphotype evolution that was independently explained by each predictor variable [185, 186], thereby controlling for covariance between time points.
Advection-Diffusion-Individual Based Models
El Niño Southern Oscillation
Generalized Linear Models
Hybrid Coordinate Ocean Model
National Oceanographic and Atmospheric Administration
Pelagic Larval Duration
United States Geological Survey
Cowen RK, Paris CB, Srinivasan A. Scaling connectivity in marine populations. Science. 2006;311:522–7.
Hastings A, Botsford L. Persistence of spatial populations depends on returning home. Proc Natl Acad Sci USA. 2006;103:6067–72.
Nosil P, Funk DJ, Ortiz-Barrientos D. Divergent selection and heterogeneous genomic divergence. Mol Ecol. 2009;18:375–402.
Castelin M, Lorion J, Brisset J, Cruaud C, Maestrati P, Utge J, Samadi S. Speciation patterns in gastropods with long-lived larvae from deep-sea seamounts. Mol Ecol. 2012;21:4828–53.
García-Ramos G, Kirkpatrick M. Genetic models of adaptation and gene flow in peripheral populations. Evolution. 1997;51:21–8.
Svenson EI, Calsbeek R. The Adaptive Landscape. Oxford: Oxford University Press; 2012.
Nagylaki T, Lou Y. Evolution under multiallelic migration-selection models. Theor Popul Biol. 2007;72:21–40.
Bürger R. Multilocus selection in subdivided populations I. Convergence properties for weak or strong migration. J Math Biol. 2009;58:939–78.
Yeaman S, Whitlock MC. The genetic architecture of adaptation under migration-selection balance. Evolution. 2011;65:1897–911.
Akerman A, Bürger R. The consequences of gene flow for local adaptation and differentiation: a two-locus two-deme model. J Math Biol. 2014;68:1135–98.
Hendry AP, Taylor EB, McPhail D. Adaptive divergence and the balance between selection and gene flow: lake and stream stickleback in the Misty system. Evolution. 2002;56:1199–216.
Bolnick DI, Nosil P. Natural selection in populations subject to a migration load. Evolution. 2007;61:2229–43.
Hunt HL, Scheibling RE. Roles of early post-settlement mortality in recruitment of benthic marine invertebrates. Mar Ecol-Progr Ser. 1997;155:269–301.
Warner RR. Evolutionary ecology: how to reconcile pelagic dispersal with local adaptation. Coral Reefs. 1997;16:S115–20.
Hamilton SL, Regetz J, Warner RR. Postsettlement survival linked to larval life in marine fish. Proc Natl Acad Sci USA. 2008;105:1561–6.
Sanford E, Kelly MW. Local adaptation in marine invertebrates. Annu Rev Mar Sci. 2011;3:509–35.
Sotka R. Natural selection, larval dispersal, and the geography of the phenotype in the sea. Integr Comp Biol. 2012;52:538–45.
Paris CB, Chérubin LM, Cowen RK. Surfing, spinning, or diving from reef to reef: effects on population connectivity. Mar Ecol-Prog Ser. 2007;347:285–300.
Treml EA, Halpin PN, Urban DL, Pratson LF. Modeling population connectivity by ocean currents, a graph-theoretic approach for marine conservation. Landscape Ecol. 2008;23:19–36.
Foster NL, Paris CB, Kool JT, Baums IB, Stevens JR, Sanchez JA, Bastidas C, Agudelo C, Bush P, Day O, Ferrari R, Gonzalez P, Gore S, Guppy R, McCartney MA, McCoy C, Mendes J, Srinivasan A, Steiner S, Vermeiji MJA, Weil E, Mumby PJ. Connectivity of Caribbean coral populations: complementary insights from empirical and modelled gene flow. Mol Ecol. 2012;21:1143–57.
Conover DO, Clarke LM, Munch SB, Wagner GN. Spatial and temporal scales of adaptive divergence in marine fishes and the implications for conservation. J Fish Biol. 2006;69:21–47.
Toonen RJ, Grosberg RK. Causes of chaos: spatial and temporal genetic heterogeneity in the intertidal anomuran crab Petrolisthes cinctipes. In: Held C, Koenemann S, Schubart CD, editors. Phylogeography and Population Genetics in Crustacea. Boca Raton: CRC Press Crustacean Issues Series; vol. 2011; 2011. p. 73–106.
Cowen RK, Lwisa KMM, Sponaugle S, Paris CB, Olson DB. Connectivity of marine populations: Open or closed? Science. 2000;287:857–9.
Siegel DA, Kinlan BP, Gaylord B, Gaines SD. Lagrangian descriptions of marine larval dispersion. Mar Ecol-Prog Ser. 2003;260:83–96.
Mora C, Sale PF. Are populations of coral reef fish open or closed? Trends Ecol Evol. 2002;17:422–8.
Shanks AL, Grantham BA, Carr MH. Propagule dispersal and the size and spacing of marine reserves. Ecol Appl. 2003;13:S159–69.
Jones GP, Planes S, Thorrold SR. Coral reef fish larvae settle close to home. Curr Biol. 2005;15:314–8.
Taylor MS, Hellberg ME. Marine radiations at small geographic scales: speciation in neotropical reef gobies (Elacatinus). Evolution. 2005;59:374–85.
Weersing K, Toonen RJ. Population genetics, larval dispersal, and connectivity in marine systems. Mar Ecol-Prog Ser. 2009;393:1–12.
Kohn YY, Clements KD. Pelagic larval duration and population connectivity in New Zealand triplefin fishes (Tripterygiidae). Environ Biol Fish. 2011;91:275–86.
Wren JLK, Kobayashi DR. Exploration of the “larval pool”: development and ground-truthing of a larval transport model off leeward Hawaii. PeerJ. 2016;4:e1636.
Gaylord B, Gaines SD. Temperature or transport? Range limits in marine species mediated solely by flow. Am Nat. 2000;155:769–89.
White C, Selkoe KA, Watson J, Siegel DA, Zacherl DC, Toonen RJ. Ocean currents help explain populations genetic structure. Proc R Soc Lond B Biol Sci. 2010;277:1685–94.
Toonen RJ, Andrews KR, Baums IB, Bird CE, Conception GT, Daly-Engel TS, Eble JA, Faucci A, Gaither MR, Iacchei M, Puritz JB, Schultz JK, Skillings DJ, Timmers MA, Bowen BW. Defining boundaries for ecosystem-based management: a multispecies case study of marine connectivity across the Hawaiian archipelago. J Mar Biol. 2011. https://doi.org/10.1155/2011/460173.
Wolanski E, Kingsford MJ. Oceanographic and behavioral assumptions in models of the fate of coral and coral reef fish larvae. J Roy Soc Interface. 2014;11:20140209.
Wong-Ala JATK, Comfort CM, Gove JM, Hixon MA, McManus MA, Powell BS, Whitney JL, Neuheimer AB. How life history characteristics and environmental forcing shape settlement success of coral reef fishes. Front Mar Sci. 2018;5:65.
Hendry AP. Selection against migrants contributes to the rapid evolution of ecologically dependent reproductive isolation. Evol Ecol Res. 2004;6:1219–36.
Nosil P, Vines TH, Funk DJ. Reproductive isolation caused by natural selection against immigrants from divergent habitats. Evolution. 2005;59:705–19.
Pinsky ML, Palumbi SR, Andréfouët S, Purkis SJ. Open and closed seascapes: where does habitat patchiness create populations with high fractions of self-recruitment? Ecol Appl. 2012;22:1257–67.
Ingley SJ, Johnson JB. Divergent natural selection promotes immigrant inviability at early and late stages of evolutionary divergence. Evolution. 2016;70:600–16.
Caley MJ, Carr MH, Hixon MA, Hughes TP, Jones GP, Menge BA. Recruitment and the local dynamics of open marine populations. Annu Rev Ecol Syst. 1996;27:477–500.
Gillanders BM. Temporal and spatial variability in elemental composition of otoliths: implications for determining stock identity and connectivity of populations. Can Jl Fish Aquat Sci. 2002;59:669–79.
Zacherl DC. Spatial and temporal variation in statolith and protoconch trace elements as natural tags to track larval dispersal. Mar Ecol-Prog Ser. 2005;290:145–63.
Pusack TJ, Christie MR, Johnson DW, Stallings CD, Hixon MA. Spatial and temporal patterns of larval dispersal in a coral-reef fish metapopulation: evidence of variable reproductive success. Mol Ecol. 2014;3:3396–408.
Grant PR, Grant BR. Unpredictable evolution in a 30-year study of Darwin’s finches. Science. 2002;296:707–11.
Reimchen TE, Nosil P. Temporal variation in divergent selection on spine number in threespine stickleback. Evolution. 2002;56:2472–83.
Siepielski AM, DiBattista JD, Carlson SM. It’s about time: the temporal dynamics of phenotypic selection in the wild. Ecol Lett. 2009;12:1261–76.
Armsworth PR. Recruitment limitation, population regulation, and larval connectivity in reef fish metapopulations. Ecology. 2002;83:1092–104.
Treml EA, Roberts JJ, Chao Y, Halpin PN, Possingham HP, Riginos C. Reproductive output and duration of the pelagic larval stage determine seascape-wide connectivity of marine populations. Integr Comp Biol. 2012;52:525–37.
Pineda JS, Hare JA, Sponaugle S. Larval transport and dispersal in the coastal ocean and consequences for population connectivity. Oceanography. 2007;20:22–39.
Mitarai S, Siegel D, Watson J, Dong C, McWilliams JC. Quantifying connectivity in the coastal ocean with application to the Southern California Bight. J Geophys Res. 2009;114:C10026. https://doi.org/10.1029/2008JC005166.
Hogan JD, Thiessen RJ, Heath DD. Variability in connectivity indicated by chaotic genetic patchiness within and among populations of a marine fish. Mar Ecol-Prog Ser. 2010;417:263–75.
Moody KN, Hunter SN, Childress MJ, Blob RW, Schoenfuss HL, Blum MJ, Ptacek MB. Local adaptation despite high gene flow in the waterfall-climbing Hawaiian goby Sicyopterus stimpsoni. Mol Ecol. 2015;24:545-563.
Munday PL, Leis LM, Lough JM, Paris CB, Kingsford MJ, Berumen ML, Lambrechts J. Climate change and coral reef connectivity. Coral Reefs. 2009;28:379–95.
Kendall MS, Poti M, Karnauskas KB. Climate change and larval transport in the ocean: fractional effects from physical and physiological factors. Glob Change Biol. 2015;22:1532–47.
Wilson LJ, Fulton CJ, Hogg AMC, Joyce KE, Radford BTM, Fraser CI. Climate-driven changes to ocean circulation and their inferred impacts on marine dispersal patterns. Glob Ecol Biogeogr. 2016;35:923–39.
Galindo HM, Olson DB, Palumbi SR. Seascape genetics: a coupled oceanographic-genetic model predicts population structure of Caribbean corals. Curr Biol. 2006;16:622–6.
Berkley H, Kendall B, Mitarai S, Siegel D. Turbulent dispersal promotes species coexistence. Ecol Lett. 2010;13:360–71.
Watson JR, Kendall BE, Siegel DA, Mitarai S. Changing seascapes, stochastic connectivity, and marine metapopulation dynamics. Am Nat. 2012;180:99–112.
Cowen RK, Sponaugle S. Larval dispersal and marine population connectivity. Annu Rev Mar Sci. 2009;1:443–66.
DeAngelis DL, Rose KA. Individual-based models and approaches in ecology. US: Springer; 1992.
Dunning JB Jr, Stewart DJ, Danielson BJ, Noon BR, Root TL, Lamberson RH, Stevens EE. Spatially explicit population models: current forms and future uses. Ecol Appl. 1995;5:3–11.
Grimm V, Railsback SF. Individual-Based Modeling and Ecology. Princeton: Princeton University Press; 2005.
Huston M, DeAngelis D, Post W. New computer models unify ecological theory. BioScience. 1988;38:682–91.
Judson OP. The rise of the individual-based model in ecology. Trends Ecol Evol. 1994;9:9–14.
Moody KN, Kawano SM, Bridges WC, Blob RW, Schoenfuss HL, Ptacek MB. Contrasting post-settlement selection results in many-to-one mapping of high performance phenotypes in the Hawaiian waterfall-climbing goby Sicyopterus stimpsoni. Evol Ecol. 2017;31:489–516.
Blob RW, Bridges WC, Ptacek MB, Maie T, Cediel RA, Bertolas MM, Julius ML, Schoenfuss HL. Morphological selection in an extreme flow environment: body shape and waterfall-climbing success in the Hawaiian stream fish Sicyopterus stimpsoni. Integr Comp Biol. 2008;48:734–49 (Erratum) 49:732-734 (2009).
Blob RW, Kawano SM, Moody KN, Bridges WC, Maie T, Ptacek MB, Julius ML, Schoenfuss HL. Morphological selection and the evaluation of potential tradeoffs between escape from predators and the climbing of waterfalls in the Hawaiian stream goby Sicyopterus stimpsoni. Integr Comp Biol. 2010;50:1185–99.
Mayr E. Animal Species and Evolution. Cambridge: Harvard University Press; 1963.
Fisher RA. The Genetical Theory of Natural Selection. Oxford: Carendon Press; 1930.
Friedlander A, Kobayashi D, Bowen B, Meyers C, Papastamatiou Y, DeMartini E, Parrish F, Treml E, Currin D, Hilting A, Weiss J, Kelley C, O’Conner R, Parke M, Clark R, Toonen R, Wedding L. Connectivity and integrated ecosystem studies. In: Friedlander A, Keller K, Wedding L, Clarke A, Monaco M, editors. A Marine Biogeographic Assessment of the Northwestern Hawaiian Islands. NOAA Technical Memorandum NOS NCCOS 84. Prepared by NCCOS’s Biogeography Branch in cooperation with the Office of National Marine Sanctuaries Papahanaumokuakea Marine National Monument. MD: Silver Spring; 2009. p. 291–330.
Lobel PS, Robinson AR. Reef fishes at sea: ocean currents and the advection of larvae. In: The Ecology of Deep and Shallow Coral Reefs, Symposium Series for Undersea Research, NOAA’s Undersea Research Program, vol. 1; 1983. p. 29–38.
Lobel PS. Ocean current variability and the spawning season of Hawaiian reef fishes. Environ Biol Fish. 1989;24:161–71.
Lobel PS. Transport of reef lizardfish larvae by an ocean eddy in Hawaiian waters. Dynam Atmos Oceans. 2011;52:119–30.
Fox HE, Haisfield KM, Brown MS, Stevenson TC, Tissot BN, Walsh WJ, Williams ID. Influences of oceanographic and meteorological features on reef fish recruitment in Hawaii. Mar Ecol-Prog Ser. 2012;463:259–72.
Vaz AC, Richards KJ, Jia Y, Paris CB. Mesoscale flow variability and its impact on connectivity for the island of Hawai‘i. Geophys Res Lett. 2013;40:332–7.
Wren JLK, Kobayashi DR, Toonen RJ. Modeled population connectivity across the Hawaiian archipelago. PLoS ONE. 2016;11:e0167626.
Concepcion GT, Baums IB, Toonen RJ. Regional population structure of Montipora capitata across the Hawaiian Archipelago. Bull Mar Sci. 2014;90:257–75.
Conklin EE, Neuheimer AB, Toonen RJ. Modeled larval connectivity of a multi-species reef fish and invertebrate assemblage off the coast of Moloka‘i, Hawai‘i. PeerJ. 2018;6:e5688.
Shima JS, Noonburg EG, Swearer SE. Consequences of variable larval dispersal pathways and resulting phenotypic mixtures to the dynamics of marine metapopulations. Biol Lett. 2015;11:20140778.
Murphy CA, Cowan JH Jr. Production, marine larval retention or dispersal, and recruitment of amphidromous Hawaiian gobies: Issues and implications. Bishop Mus Bull Cult Environ Stud. 2007;3:63–74.
Lobel PS, Robinson AR. Transport and entrapment of fish larvae by ocean mesoscale eddies and currents in Hawaiian water. Deep-Sea Res. 1986;33:483–500.
Price JF, Weller RA, Schudlich RR. Wind-driven ocean currents and Ekman transport. Science. 1987;238:1534–8.
Bryden HL, Roemmich DH, Church JA. Ocean heat transport across 24 N in the Pacific. Deep-Sea Res. 1991;38:297–324.
Chang P, Ji L, Li H. A decadal climate variation in the tropical Atlantic Ocean from thermodynamic air-sea interactions. Nature. 1997;385:516–8.
Stommel H. Varieties of oceanographic experience. Science. 1963;139:572–6.
Pineda J, López M. Temperature, stratification and barnacle larval settlement in two Californian sites. Cont Shelf Res. 2002;22:1183–98.
Edmunds PJ. Unusually high coral recruitment during the 2016 El Nino in Mo’orea, French Polynesia. PLoS ONE. 2017;12:e0185167.
Higgins BA, Pearson D, Mehta RS. El Niño episodes coincide with California moray Gymnothorax mordax settlement around Santa Catalina Island, California. J Fish Biol. 2017;90:1570–83.
Wood S, Baums IB, Paris CB, Ridgwell A, Kessler WS, Hendy EJ. El Niño and coral larval dispersal across the eastern Pacific marine barrier. Nat Commun. 2016;7:12571.
Paris CB, Cowen RK. Direct evidence of a biophysical retention mechanism for coral reef fish larvae. Limnol Oceanogr. 2004;49:964–1979.
Moody KN, Gagne RB, Heim-Ballew H, Alda F, Hain EF, Lisi PJ, Walter RP, Higashi GR, Hogan JD, McIntyre PB, Gilliam JF, Blum MJ. Invasion hotspots and ecological saturation of streams across the Hawaiian archipelago. Cybium. 2017;41:127–56.
McDowall RM. On amphidromy, a distinct form of diadromy in aquatic organisms. Fish Fish. 2007;8:1–13.
Brasher AM. Impacts of human disturbances on biotic communities in Hawaiian streams. BioScience. 2003;53:1052–60.
Resh VH. Stream crossings and the conservation of diadromous invertebrates in South Pacific island streams. Aquat Conserv. 2004;15:313–7.
Ramírez A, Engman A, Rosas KG, Perez-Reyes O, Martinó-Cardona DM. Urban impacts on tropical island streams: Some key aspects influencing ecosystem response. Urban Ecosyst. 2012;15:315–25.
Lisi PJ, Childress ES, Gagne RB, Hain EF, Lamphere BA, Walter RP, Hogan JD, Gilliam JF, Blum MJ, McIntyre PB. Overcoming urban stream syndrome: Trophic flexibility confers resilience in a Hawaiian stream fish. Freshwater Biol. 2018;63:492–502.
Nanninga GB, Berumen ML. The role of individual variation in marine larval dispersal. Front Mar Sci. 2014;1:71.
Heinz S, Mazzucco F, Dieckmann U. Speciation and the evolution of dispersal along environmental gradients. Evol Ecol. 2009;23:53–70.
Nurmi T, Parvinen K. Joint evolution of specialization and dispersal in structure metapopulations. J Theor Biol. 2011;275:78–92.
Berdahl A, Torney CJ, Schertzer E, Levin SA. On the evolutionary interplay between dispersal and local adaptation in heterogeneous environments. Evolution. 2015;69:1361–653.
Hogan JD, McIntyre PB, Blum MJ, Gilliam JF, Bickford N. Consequences of alternative dispersal strategies in a putatively amphidromous fish. Ecology. 2014;95:2397–408.
Sorensen PW, Hobson KA. Stable isotope analysis of amphidromous Hawaiian gobies suggests their larvae spend a substantial period of time in freshwater river plumes. Environ Biol Fish. 2005;74:31–42.
Tsugagawa T, Arai T. Flexible migration of Japanese freshwater gobies Rhinogobius spp. as revealed by otolith Sr:Ca ratios. J Fish Biol. 2008;73:2421–33.
Dodson JJ, Aubin-Horth N, Thériault V, Páez DJ. The evolutionary ecology of alternative migratory tactics in salmonid fishes. Biol Rev. 2013;88:602–25.
Hughes JM, Schmidt DJ, MacDonald JI, Huey JA, Crook DA. Low interbasin connectivity in a facultatively diadromous fish: evidence from genetics and otolith chemistry. Mol Ecol. 2014;23:1000–13.
Augspurger JM, Warburton M, Closs GP. Life-history plasticity in amphidromous and catadromous fishes: a continuum of strategies. Rev Fish Biol Fisher. 2017;27:177–92.
Hicks AS, Jarvis MG, David BO, Waters JM, Norman MC, Closs GP. Lake and species specific patterns of non-diadromous recruitment in amphidromous fish: the importance of local recruitment and habitat requirements. Mar Freshwater Res. 2017;68:2315–23.
Dixson DL, Jones GP, Munday PL, Planes S, Pratchett MS, Srinivasan M, Syms C, Thorrold SR. Coral reef fish smell leaves to find island homes. Proc R Soc Lond B Biol Sci. 2008;275:2831–9.
Berumen ML, Almany GR, Planes S, Jones GP, Saenz-Agudelo P, Thorrold SR. Persistence of self-recruitment and patterns of larval connectivity in a marine protected area network. Ecol Evol. 2012;2:444–52.
Warburton ML, Jarvis MG, Closs GP. Otolith microchemistry indicates regional phylopatry in the larval phase of an amphidromous fish (Gobimorphus hubbsi). New Zeal J Mar Fresh. 2018;52:398–408.
Radtke RL, Kinzie RA III. Evidence of a marine larval stage in endemic Hawaiian stream gobies from isolated high-elevation locations. T Am Fish Soc. 1996;135:613–21.
Bell KNI, Brown JA. Active salinity choice and enhanced swimming endurance in 0 to 8-d-old larvae of diadromous gobies, with emphasis on Sicydium punctatum (Pisces), in Dominica West Indies. Mar Biol. 1995;121:409–17.
Alda F, Gagne RB, Walter RP, Hogan JD, Moody KN, Zink F, McIntyre PB, Gilliam JF, Blum MJ. Colonization and demographic expansion of freshwater fauna across the Hawaiian Archipelago. J Evol Biol. 2016;29:2054–69.
Carine MA. Spatio-temporal relationships of the Macaronesian flora: a relictual series or window of opportunity? Taxon. 2005;54:895–903.
Hedgecock D. Is gene flow from pelagic larval dispersal important in the adaptation and evolution of marine invertebrates? Bull Mar Sci. 1986;39:550–64.
Doherty PJ, Dufour V, Galzin R, Hixon MA, Meekan MG, Planes S. High mortality during settlement is a population bottleneck for a tropical surgeonfish. Ecology. 2004;85:422–8.
Kawano SM, Bridges WC, Schoenfuss HL, Maie T, Blob RW. Differences in locomotor behavior correspond to different patterns of morphological selection in two species of waterfall-climbing gobiid fishes. Evol Ecol. 2013;275:949–69.
Jones AG. A theoretical quantitative genetic study of negative ecological interactions and extinction times in changing environments. BMC Evol Biol. 2008;8:119.
Tseng M, O’Connor MI. Predators modify the evolutionary response of prey to temperature change. Biol Lett. 2015;11:20150798.
Osmond MM, Otto SP, Klausmeier CA. When predators help prey adapt and persist in a changing environment. Am Nat. 2017;190:83–98.
Langerhans RB, Lyman CA, Shokrollahi AM, DeWitt TJ. Predator-driven phenotypic diversification in Gambusia affinis. Evolution. 2004;58:2305–18.
Freeman AS, Byers JE. Divergent induced responses to an invasive predator in marine mussel populations. Science. 2006;313:831–3.
Losos JB, Schoener TW, Langerhans RB, Spiller DA. Rapid temporal reversal in predator-driven natural selection. Science. 2006;314:1111.
Nosil P, Crespi BJ. Experimental evidence that predation promotes divergence in adaptive radiation. Proc Natl Acad Sci USA. 2006;103:909–9095.
Marchinko KB. Predation’s role in repeated phenotypic and genetic divergence of armor in threespine stickelback. Evolution. 2009;63:127–38.
Langerhans RB. Trade-off between steady and unsteady swimming underlies predator-driven divergence in Gambusia affinis. J Evol Biol. 2009;22:1057–75.
Ingley SJ, Billman EJ, Belk MC, Johnson JB. Morphological divergence driven by predation environment within and between species of Brachyrhaphis fishes. PLoS ONE. 2014;9:e90274.
Swaegers J, Strobbe F, McPeek MA, Stoks R. Selection on escape performance during ecological speciation driven by predation. Anim Behav. 2017;124:153–49.
Sax DV, Gaines SD. Species invasions and extinction: The future of native biodiversity on islands. Proc Natl Acad Sci USA. 2008;105:11490–7.
Kingsford RT, Watson EM, Lundquist CJ, Venter O, Hughes L, Johnston EL, Atherton J, Gawel M, Keith DA, Mackey BG, Morley C, Possingham HP, Raynor B, Recher HF, Wilson KA. Major conservation policy issues for biodiversity in Oceania. Conserv Biol. 2009;23:834–40.
Courchamp F, Hoffmann BD, Russell JC, Leclerc C, Bellard C. Climate change, sea-level rise, and conservation: keeping island biodiversity afloat. Trends Ecol Evol. 2014;29:127–30.
Klasner F, Mikami CD. Land use on the island of Oahu, 1998. Honolulu: U.S. Geological Survey, Water Resources Investigations, Report 02-4301; 2003. p. 20.
Oki DS, AMD B. Environmental setting and implications for water quality and aquatic biota, Oahu, Hawaii. Honolulu: U.S. Geological Survey, Water Resources Investigations, Report 03-4156; 2003. p. 98.
Blum MJ, Gilliam JF, McIntyre PB. Development and use of genetic methods for assessing aquatic environmental condition and recruitment dynamics of native stream fishes on Pacific Islands. United States Department of Defense, SERDP RC-1646 2014.
Delevaux JMS, Whittier R, Stamoulis KA, Bremer LL, Jupiter S, Friedlander AM, Poti M, Guannel G, Kurashima N, Friedlander AM, Toonen R, Conklin E, Wiggins C, Knudby A, Goodell W, Burnett K, Yee S, Htun H, Oleson KLL, Wiegner T, Ticktin T. A linked land-sea modeling framework to inform ridge-to-reef management in high oceanic islands. PLoS ONE. 2018;13:e0193230.
Witter MS, Carr GD. Adaptive radiation and genetic differentiation in the Hawaiian silversword alliance (Compositae:Madiinae). Evolution. 1988;42:1278–87.
Lerner HRL, Meyer M, James HF, Hofreiter M, Fleischer RC. Multilocus resolution of phylogeny and timescale in the extant adaptive radiation of Hawaiian honeycreepers. Current Biol. 2011;21:1838–44.
Givish TJ, Millam KC, Mast AR, Paterson TB, Theim TJ, Hip AL, Henss JM, Smith JF, Wood KR, Sytsma KJ. Origin, adaptive radiation and diversification of the Hawaiian lobelaids (Asterales: Campanulaceae). Proc R Soc Lond B Biol Sci. 2009;276:407–16.
Jokiei PL. Ecology, biogeography and evolution in corals in Hawaii. Trends Ecol and Evol. 1987;2:179–82.
Kay EA, Palumbi SR. Endemism and evolution in Hawaiian marine invertebrates. Trends Ecol and Evol. 1987;2:183–6.
Hourigan TF, Reese ES. Mid-ocean isolation and the evolution of Hawaiian reef fishes. Trends Ecol and Evol. 1987;2:187–91.
Thacker CE. Phylogeny of Gobioidei and placement within Acanthomorpha, with a new classification and investigation of diversification and character evolution. Copeia. 2009;2009(1):92–104.
Bird CE. Morphological and behavioral evidence for adaptive diversification of sympatric Hawaiian limpets (Cellana spp.). Integr Comp Biol. 2011;51:466–73.
Fitzsimons JM, Zink RM, Nishimoto RT. Genetic variation in the Hawaiian stream goby Lentipes concolor. Biochem Syst Ecol. 1990;18:81–3.
Rocha LA, Robertson DR, Roman J, Bowen BW. Ecological speciation in tropical reef fishes. Proc R Soc Lond B Biol Sci. 2005;272:573–9.
Bowen BW, Rocha LA, Toonen RJ, Karl SA. The origins of tropical marine biodiversity. Trends Ecol Evol. 2013;28:359–66.
Momigliano P, Jokinen H, Fairmout A, Florin AB, Norkko A, Merilä J. Extraordinary rapid speciation in a marine fish. Proc Natl Acad Sci. 2017;114:6074–9.
Palumbi SR. Marine reserves and ocean neighbourhoods: the spatial scale of marine populations and their management. Annu Rev Env Resour. 2003;29:31–68.
Marko PB. ‘What’s larvae got to do with it?’ Disparate patterns of post-glacial population structure in two benthic marine gastropods with identical dispersal potential. Mol Ecol. 2004;13:597–611.
Siegel DA, Mitarai S, Costello CJ, Gaines SD, Kendall BE, Warner RR, Winters KB. The stochastic nature of larval connectivity among nearshore marine populations. Proc Natl Acad Sci. 2008;105:8974–9.
Teske PR, Papadopoulos I, Newman BK, Dworschak PC, McQuaid CD, Barker NP. Oceanic dispersal barriers, adaptation and larval retention: an interdisciplinary assessment of potential factors maintaining a phylogeographic break between sister lineages of an African prawn. BMC Evol Biol. 2008;8:341–54.
D’Aloia CC, Bogdanowicz SM, Francis RK, Majoris JE, Harrison RG, Buston PM. Patterns, causes, and consequences of marine larval dispersal. Proc Natl Acad Sci. 2015;112:13940–5.
Cuif M, Kaplan DM, Fauvelot C, Lett C, Vigliola L. Monthly variability of self-recruitment for a coral reef damselfish. Coral Reefs. 2015;34:759–70.
Almany GR, Planes S, Thorrold SR, Berumen ML, Bode M, Saenz-Agudelo P, Bonin MC, Frisch AJ, Harrison HB, Messmer V, Nanninga GB. Larval fish dispersal in a coral-reef seascape. Nat Ecol Evol. 2017;1:0148.
Ego K. Life history of freshwater gobies. Freshwater game fish management research. Honolulu, Job Completion Report F-4-R: Territory of Hawaii, Department of Land and Natural Resources; 1956.
Kinzie RA III. Species profiles: life histories and environmental requirements of coastal vertebrates and invertebrates, Pacific Ocean Region; Report 3, Amphidromous macrofauna of Hawaiian island streams. Vicksburg: Report to U.S. Army Engineers Waterways Experiment Station; 1990. 52 p.
Smith RJF, Smith J. Rapid acquisition of directional preferences by migratory juveniles of two amphidromous Hawaiian gobies, Awaous guamensis and Sicyopterus stimpsoni. Environ Biol Fish. 1998;53:275–82.
McDowall RM. Hawaiian biogeography and the islands’ freshwater fish fauna. J Biogeogr. 2003;30:703–10.
Lindstrom DL, Blum MJ, Walter RP, Gagne R, Gilliam JF. Molecular and morphological evidence of distinct evolutionary lineages of Awaous guamensis in Hawai'i and Guam. Copeia. 2012;2:294–301.
Fitzsimons JM, Nishimoto RT, Yuen AR. Courtship and territorial behavior in the native Hawaiian stream goby Sicyopterus stimpsoni. Ichthyol Explor Fres. 1993;4:1–10.
Schoenfuss HL, Blob RW. The importance of functional morphology for fishery conservation and management applications to Hawaiian amphidromous fishes. Bishop Mus Bull Cult Environ Stud. 2007;3:125–41.
Clague DA, Dalrymple GB. Tectonics, geochronology and origin of the Hawaiian-Emperor volcanic chain. In: Decker RW, Wright TL, Stauffer PH, editors. Volcanism in Hawaii. US Geological Survey Professional Paper 1350. Washington DC: US Government Printing Office; 1987. p. 1–54.
Norton SF, Luczkovich JJ, Motta PJ. The roles of ecomorphological studies in comparative biology of fishes. Environ Biol Fish. 1995;44:287–304.
Domenici P. Habitat, body design, and the swimming performance of fish. In: Bels VL, Gasc JP, Casinos A, editors. Vertebrate Biomechanics and Evolution, vol. 2003. Oxford: BIOS Scientific Publishers; 2003. p. 137–60.
Walker JA. Ecological morphology of lacustrine threespine stickleback Gasterosteus aculeatus L. (Gasterosteidae) body shape. Biol J Linn Soc. 1997;61:3–50.
Polovina J, Kleiber P, Kobayashi D. Application of TOPEX-POSEIDON satellite altimetry to simulate transport dynamics of larvae of spiny lobster, Panulirus marginatus, in the Northwestern Hawaiian Islands, 1993-1996. Fish Bull. 1999;97:132–43.
Rivera M, Andrews K, Kobayashi D, Wren J, Kelley D, Roderick G, Toonen R. Genetic analyses and simulations of larval dispersal reveal distinct populations and directional connectivity across the range of the Hawaiian grouper (Epinephelus quernus). J Mar Biol. 2011. https://doi.org/10.1155/2011/765353.
Jai Y, Calil PHR, Chassignet EP, Metzger EJ, Potemra JT, Richards KJ, Wallcraft AJ. Generation of mesoscale eddies in the lee of the Hawaiian Islands. J Geophys Res-Oceans. 2011;116:1978–2012.
Lumpkin CF. Eddies and currents of the Hawaiian Islands. Doctoral dissertation: University of Hawaii; 1998.
Boehlert GW, Mundy BC. Ichthyoplankton vertical distributions near Oahu, Hawaii, 1985-1986: data report. NOAA Technical Memorandum NMFS-TM_NMFS-SWFSC-235; 1996. p. 1–148.
Radke RL, Kinzie RA, Folsom SD. Age at recruitment of Hawaiian freshwater gobies. Environ Biol Fish. 1988;23:205–13.
Radtke RL, Kinzie RA, Shafer DJ. Temporal and spatial variation in length of larval life and size at settlement of the Hawaiian amphidromous goby Lentipes concolor. J Fish Biol. 2001;59:928–38.
Lord C, Brun C, Hautecoeur M, Keith P. Insights on endemism: comparison of the duration of the marine larval phase estimated by otolith microstructural analysis of three amphidromous Sicyopterus species (Gobioidei: Sicydiinae) from Vanuatu and New Caledonia. Ecol Freshw Fish. 2010;19:26–38.
Wilensky U. NetLogo. Center for Connected Learning and Computer-Based Modeling, Northwestern University 1999. http://ccl.northwestern.edu/netlogo/.
Fitzsimons JM, McRae MG. Behavioral ecology of indigenous stream fishes in Hawaii. Bishop Mus Bull Cult Environ Stud. 2007;3:11–21.
Tate DC. The role of behavioral interactions of immature Hawaiian stream fishes (Pisces: Gobioidei) in population dispersal and distribution. Micronesica. 1997;30:51–70.
Nishimoto RT, Kuamo’o DGK. Recruitment of goby postlarvae into Hakalau Stream, Hawai‘i Island. Micronesica. 1997;30:41-49.
Mousseau TA, Roff DA. Natural selection and the heritability of fitness components. Heredity. 1987;59:181–97.
Dingemanse NJ, Van der Plas F, Wright J, Réale D, Schrama M, Roff DA, Van der Zee E, Barber I. Individual experience and evolutionary history of predation affect expression of heritable variation in fish personality and morphology. Proc R Soc Lond B Biol Sci. 2009;276:1285–93.
Price JP, Clague DA. How old is the Hawaiian biota? Geology and phylogeny suggest recent divergence. Proc R Soc Lond B Biol Sci. 2002;269:2429–35.
Lande R, Arnold SJ. The measurement of selection on correlated characters. Evolution. 1983;37:1210–26.
Rao CR. The use and interpretation of principal component analysis in applied research. Sankhya: Indian J Stat Ser A. 1964;26:329–58.
Oksanen J, Blanchet FG, Friendly M, Kindt R, Legendre P, McGlinn D, Minchin PR, O’Hara RB, Simpson GL, Solymos P, Stevens MHH, Szoecs E, Wagner H. vegan: Community ecology package. R package version 2.4-1. 2016. https://CRAN.R-project.org/package=vegan.
Legendre P. Studying beta diversity: ecological variation partitioning by multiple regression and canonical analysis. J Plant Ecol. 2008;1:3–8.
Peres-Neto PR, Legendre P. Estimating and controlling for spatial structure in the study of ecological communities. Glob Ecol Biogeogr. 2010;19:174–84.
We gratefully acknowledge the Hawai‘i Division of Aquatic Resources- especially Bob Nishimoto, Lance Nishiura, Wade Ishikawa, Troy Sakihara, Troy Shimoda, Tim Shindo, and Darrel Kuamo‘o- for logistical support during field work done for this study. We would also like to thank the anonymous reviewers who provided valuable guidance and recommendations.
Funding was provided through National Science Foundation Doctoral Dissertation Improvement Grant (DEB-1310962) to Kristine N. Moody, the Department of Defense Strategic Environmental Research and Development Program awards RC-1646 and RC-2447 to Michael J. Blum, NSF grant IOS-0817794 to Richard W. Blob and Margaret B. Ptacek, NSF grant IOS-0817911 to Heiko L. Schoenfuss and Matthew L. Julius, NSF grant OCE-1260169 to Robert J. Toonen, and National Oceanic and Atmospheric Administration sponsored by the University of Hawaii Sea Grant College Program (UNIHI_SEAGRANT R/EL-43 and R/IR-32 under Project R/SS-13), SOEST, under Institutional Grant No. NA14OAR4170071 from NOAA Office of SeaGrant, Department of Commerce, to Robert J. Toonen. The funding bodies had no role in the design of the study and collection, analysis, and interpretation of data or in writing the manuscript.
Availability of data and materials
Data of model outputs and model codes are available through Open Science Framework (https://doi.org/10.17605/OSF.IO/M5SBH).
Ethics approval and consent to participate
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Figure S1. Simulated counts of larval morphotypes for 200 generations on the islands of Kaua‘i (a), O‘ahu (b), and the Big Island (c) from the individuals-based models of isolation with varying levels (0 to 1) of post-settlement selection of predation (columns) and climbing (rows) (scenario 2). Warm colors represent climbing morphotypes (M1-M4) and cool colors represent predation evasion morphotypes (M7-M10). (PDF 4133 kb)
Figure S2. Simulated counts of immigrant larval morphotypes for 200 generations on the islands of Kaua‘i (a), O‘ahu (b), and the Big Island (c) from the individual-based models of immigration (ranging from 25%-100%) with varying levels (0.25 to 1) of post-settlement predation selection (scenario 4). Warm colors represent climbing morphotypes (M1-M4) and cool colors represent predation evasion morphotypes (M7-M10). (PDF 4221 kb)
Figure S3. Simulated counts of self-recruitment larval morphotypes for 200 generations on the islands of Kaua‘i (a), O‘ahu (b), and the Big Island (c) from the individual-based models of immigration (ranging from 25%-100%) with varying levels (0.25 to 1) of post-settlement predation selection (scenario 4). Warm colors represent climbing morphotypes (M1-M4) and cool colors represent predation evasion morphotypes (M7-M10). (PDF 4769 kb)
Supplemental Equations. Sine curve equations of monthly stream discharge estimates from UGSS stream gauge data. (DOCX 101 kb)
Figure S4. Predicted morphotype frequency distribution and island of origin of settled pelagic larvae from the dispersal. These frequency distributions were used as the inputs for pelagic larval morphology (Stage 1) in the individual-based models that included immigration (scenarios 3, and 4). (PDF 36 kb)