- Research article
- Open Access
Evolutionary patterns of diadromy in fishes: more than a transitional state between marine and freshwater
BMC Evolutionary Biologyvolume 19, Article number: 168 (2019)
Across the tree of life there are numerous evolutionary transitions between different habitats (i.e., aquatic and terrestrial or marine and freshwater). Many of these dramatic evolutionary shifts parallel developmental shifts that require physiological, anatomical and behavioral changes for survival and reproduction. Diadromy (scheduled movement between marine and freshwater) has been characterized as a behavior that acts as an evolutionary intermediate state between marine and freshwater environments, implying that diadromous lineages are evolutionarily transient. This hypothesis comes with assumptions regarding the rates of evolutionary transitions in and out of diadromy as well as rates of speciation and extinction in diadromous fishes.
Based on a published phylogeny of 7822 species of ray-finned fishes, state speciation and extinction models of evolutionary transition between marine, freshwater, and diadromous species suggest transition rates out of diadromy are 5–100 times higher that transition between marine and freshwater or into diadromy. Additionally, high speciation and low extinction rates separate diadromous fishes from marine and freshwater species. As a result, net diversification (net diversification = speciation – extinction) is about 7–40 times higher in diadromous fishes compared to freshwater and marine respectively. Together the transition, speciation, and extinction rates suggest diadromy is the least stable of the three states.
Evolutionary transitions to diadromy are rare in fishes. However, once established, diversification rates in diadromous lineages are high compared to both marine and freshwater species. Diadromous lineages tend to be more transient than marine or freshwater lineages and are found to give rise to marine and freshwater specialists in addition to diadromous descendants. Although diadromy is not a necessary evolutionary intermediate between marine and freshwater, these results support the interpretation of diadromy as an important, occasionally intermediate state, that contributes to biodiversity in fishes in all environments. This evolutionary instability of diadromous lineages is counteracted by their relatively high diversification rates. These findings highlight the importance of integrating the dynamics of diversification and major evolutionary transitions for understanding macroevolutionary patterns.
Major evolutionary transitions between different habitat types are seen across the tree of life. Transitions between water and land , flightlessness and flight , and marine and freshwater  all require major behavioral, physiological and morphological changes, yet each transition has occurred numerous times [1,2,3,4,5,6,7,8,9]. Explaining the high frequency of these transitions is an important challenge for evolutionary biologists, as these transitions allow species to move into niches not previously occupied. Identifying specific individuals or species that make these transitions as part of their life history could be key in understanding the evolutionary transitions. Life histories in which these major transitions occur ontogenetically are found in many taxonomic groups, including salamanders that reproduce in water and live on land; fishes, crabs and snails that migrate between marine and freshwater; sea turtles that are born on land and hatchlings immediately move to the ocean; and many semi-aquatic insects, such as dragonflies, that have an aquatic larval stage but a terrestrial adulthood with flight. While it is convenient to find examples of species with transitional life histories at the evolutionary transition between lineages occupying different habitats, explicit tests are needed to understand if transitional life histories are acting as an evolutionary intermediate (i.e., a life history that occurs as an intermediate between two other life histories such as semiterrestrial in the case of aquatic ➔ semiterrestrial ➔ terrestrial evolutionary transitions). Here I test the idea that the fishes that migrate between marine and freshwater habitats within their lifetime are also the evolutionary intermediates between completely marine and completely freshwater lineages.
Among fishes, one of the most dramatic changes in habitat type is between marine and freshwater. Diadromy, or the migration between marine and freshwater during a particular life stage [4, 6], requires numerous physiological changes in osmotic and ionic regulation and differs from euryhaline fishes that move across a wide gradient of salinity with regularity [4, 6, 10, 11]. This scheduled movement between environments is of ecological and evolutionary interest, as it is rare (0.8% of the 32,000 species of ray-finned fishes), yet phylogenetically widespread (8% of the nearly 500 families of fishes)  (Fig. 1). In addition, diadromy is observed in other animals including snails , prawns , and crabs [14, 15].
Multiple hypotheses have been proposed to explain the evolutionary significance of diadromy [5, 8, 9, 16]. Many studies focus on the evolutionary transitions in and out of diadromy and propose that diadromy acts as an intermediate stage between marine life histories and freshwater life histories. While the proposed mechanisms range from resource availability for adults based on latitudinal variation in net primary productivity of environments [5, 8] to moving into a safe site for spawning , each of these hypotheses suggests a similar evolutionary trajectory: A marine species has individuals or lineages that venture into freshwater (be it during adulthood or spawning) and transition into a diadromous species. Over time, lineages of the diadromous species stop returning to marine systems, becoming completely freshwater due to some advantage(s) of the new environment (i.e., marine ➔ diadromy ➔ freshwater). Similar, arguments have been constructed in the opposite direction (i.e., freshwater ➔ diadromy ➔ marine) [5, 8]. The natural world provides many examples of evolutionary transitions where diadromy acts as an intermediate, such as sticklebacks (Gasterosteidae), lampreys (Petromyzontiformes), salmonids (Salmonidae), and shads (Clupeidae) [6, 17, 18]. While these examples support the directionality in the previously proposed hypotheses [3, 18], there are also examples of evolutionary transitions involving diadromous species not acting as an intermediate such as freshwater eels (Anguillidae), herrings (Clupeidae), and flagtails (Kuhliidae) [3, 6, 17, 19, 20]. Some of these examples, such as the flagtails, transition from marine to diadromy and back to marine . There are also a number of putative examples of transitions directly between marine and freshwater [21,22,23,24]. These numerous exceptions suggest that diadromy might be more than an evolutionary intermediate.
Modern phylogenetic comparative methods can utilize time calibrated phylogenies to estimate rates of evolutionary transition between life histories (q), speciation (λ) and extinction (μ). These rates can then be incorporated into the hypotheses of diadromy as an evolutionary intermediate. In addition to testing the rates of evolutionary transitions in and out of diadromy, if diadromy is solely a transitional form, it would be expected that diadromous species would show little variation in speciation and extinction rates compared to marine and freshwater species. Drastic differences in speciation and/or extinction rates would suggest that the life history of diadromy affects the net diversification rate (r = λ - μ) of fishes. In fact, diadromous behavior has been proposed to have seemingly conflicting effects on genetic structure and the two components of net diversification [6, 25]. On one hand, diadromous migration between different river catchments is expected to decrease population structure compared to completely freshwater species which cannot move between river systems via the ocean. This should decrease the risk of extinction due to unfavorable habitat changes or stochastic events that would wipe out a completely freshwater species that is fully restricted to that freshwater system. At the same time diadromy could decrease the rate of speciation due to increased gene flow between isolated river systems via migration into and back out of the ocean. On the other hand, diadromous migrations are expected to increase the geographic range of a species compared to completely freshwater taxa, potentially increasing isolation-by-distance and changing the selection pressure put on a population as it expands into new habitats. This might increase rates of ecological speciation via local adaptation or vicariant speciation via landlocking or other barriers to dispersal which can disrupt widely distributed species [6, 13, 25, 26]. There are also conflicting patterns compared to marine species. Diadromous species might have increased population subdivision due to even a small restriction to freshwater habitats for a portion of their lives. This should increase speciation rates in diadromous species relative to marine species (although some marine reef associated fishes have extremely limited dispersal and a high level of population structure [27,28,29]). Comparative analyses on a few select taxa suggest that diadromous fishes have an intermediate level of population structure compared to the more structured freshwater and less structured marine species . Based on the compounding factors including range size shifts, population differentiation, and geneflow in each instance of diadromy, diversification patterns should vary from marine and freshwater suggesting that diadromy may be responsible for shifts in diversification rates, not just a transitional phase. Coupling the assumptions about habitat use associated with transition patterns and the assumptions about speciation and extinction associated with population dynamics may shed light on the true nature of diadromy.
Previous phylogenetic analyses have addressed the evolutionary transitions in and out of diadromy but tend to focus on specific families or orders of fishes [3, 20, 23, 31,32,33,34]. However, these studies can suffer from low statistical power owing to small numbers of evolutionary transitions and little time for speciation and extinction to occur (e.g., Feutry et al.  examine one transition into diadromy and one transition out of diadromy). Some clade-specific studies find no significant difference between marine and freshwater transition rates , while other fish groups show considerably more transitions from marine to freshwater compared to the reverse . Some groups have few transitions into diadromy [6, 20, 34], while others have numerous transitions both in and out of diadromy [3, 6]. These studies also tend to focus on clades with relatively high prevalence of diadromy. Thus, they might provide insight into the timing and mechanisms of particular transitions by examining a few species within an ecological context but have limited scope for generalization to all fishes. To better assess the evolutionary significance of diadromy, this study uses the most complete and time-calibrated tree of 7822 fishes constructed using 13 genes and 60 fossils for time calibration [12, 35] coupled with MuSSE (Multiple State Speciation and Extinction) , a well-established method, to better elucidate the diversification within and transition between three states (i.e., marine, freshwater and diadromous) across ray-finned fishes. The MuSSE model uses a time-calibrated phylogeny to generate rate estimations of speciation (λ) and extinction (μ) within each of the three states, as well as rates of transition between each state (q).
The first aim of this study is to test the hypothesis that diadromy is an intermediate evolutionary stage between marine and freshwater lifestyles (Fig. 2), including a test of the extreme scenario that all transitions between marine and freshwater have to go through diadromy (Fig. 2I). This extreme scenario is contradicted by cases of apparent transition between marine and freshwater without evidence of a diadromous intermediate [22, 24]; however, large scale phylogenetic analysis could support the existence of diadromy as a temporary or “hidden” transitional state. A more realistic null model is one with all transition rates being similar regardless of the direction (Fig. 2II). Additionally, if marine ➔ diadromy ➔ freshwater happens at a different rate from freshwater ➔ diadromy ➔ marine, shifts in rates based on directionality should be observed (Fig. 2III).
The second aim of this study is to test the hypothesis that, even though not all transitions between marine and freshwater go through diadromy, and regardless of the directionality of transition in and out of diadromy, diadromy only acts as a transient stage, resulting in a net diversification rate in diadromous fishes that is lower than either marine or freshwater fishes. This hypothesis suggests that partial movement into new habitats should not result in speciation or extinction, but instead transition out of diadromy (Fig. 2I). An alternative hypothesis predicts an intermediate or higher diversification rate arising from the expected effects of diadromy on population structure, range size, gene flow and extinction risk. The two aims of this study should also be looked at together. The idea that diadromy is primarily a transitional state would be supported if transitions involving diadromy are quantitatively dominant and diversification of diadromous lineages are relatively unimportant. It is important to note that both hypotheses are strictly pattern-oriented assessments of rates and do not address causality in the system. As might be the case in family or order specific studies, no attempt is made here to test why or how diadromy evolves. This study focuses on the quantitative importance of diadromy by assessing transition, speciation, and extinction rates of marine, diadromous and freshwater fishes.
Phylogenetic hypothesis testing
The transition rate from diadromy to freshwater (qDF = 9.37 × 10− 3) is higher than transitions in the reverse direction (qFD = 4.98 × 10− 4) and the transition rate from diadromy to marine (qDM = 4.81 × 10− 1) are also higher than transitions in the reverse direction (qMD = 2.87 × 10− 3). The highest transition rate overall reflected movement from diadromous to marine lifestyles (Table 1A). Transition rates inferred from fitDiscrete analyses corroborated the MuSSE results showing high transition rates out of diadromy. The ARD model was best fit for both the two state and three state models (Table 1B). Both of the best fit models inferred transition rates out of diadromy that were larger than any other transition rate. Regardless of whether the “unknown” species were assigned to marine or freshwater states, the estimated parameter rates from MuSSE remained similar (Additional file 3). Posterior distributions from MuSSE show strong variation between diadromous and both non-diadromous states for both speciation and extinction rates as well as transition out of diadromy as opposed to all other transitions (Fig. 3). BiSSE and HiSSE justified the assumption of “diadromy” as a group for MuSSE analysis. BiSSE and HiSSE analyses resulted in no support for a hidden state implying no sub groupings within diadromous fishes (i.e., each diadromous fish was grouped into the “diadromy-observed” state with > 98% likelihood) (Additional file 4). HiSSE did assign non-diadromous fishes into two rate categories that did not match up with character assignment of marine and freshwater species. Because of the differences in the nature of the three state MuSSE (i.e., diadromous, marine and freshwater) and the two-state BiSSE and HiSSE models (i.e., diadromous vs. non-diadromous), rate comparison cannot be made between the two model types. Additionally, neither BiSSE nor HiSSE can assess diadromy with respect to its role as a transitional state between marine and freshwater because they are combined into one state.
Based on MuSSE, the net diversification rate (r = λ – μ) in diadromous fishes (r = 0.446) is about an order of magnitude higher than in both marine (r = 0.01) and freshwater (r = 0.062) fishes (Table 1, Additional file 2). The observed pattern of freshwater fishes having a higher speciation and extinction rate compared to marine fishes is in line with estimates in New World silversides . State Stability (S) was lowest for diadromy (SD = 0.51 = 1–2.33e− 7-9.37e− 3-4.81e− 1) compared to marine (SM = 0.87 = 1–1.29e− 1-3.36e− 4-2.87e− 3) and freshwater (SF = 0.84 = 1–1.55e− 1-1.77e− 4-4.98e− 4).
Results indicate that diadromy is a rare state (low transition rates into diadromy) that, if obtained, can lead to rapid diversification (high speciation and low extinction rates) and large numbers of descendants specialized in either marine or freshwater (high transition rates out of diadromy) (Fig. 2IV). If the intermediate state hypotheses [5, 8] were true, we would expect the models constraining transition through diadromy (i.e., qFD~qDM, qFM~qDM) or those limiting movement between freshwater and marine (i.e., qFM~ 0 and/or qMF~ 0) to all be better fit models (Additional file 2). In the best-supported model, transitions occur between all states. The top models show highest transition levels while exiting diadromy (Fig. 2IV) and do not show patterns expected from simple interpretations of diadromy as an intermediate state (Fig. 2I, III). Additionally, no model in which the speciation and/or extinction rates in diadromy were constrained to equal marine and/or freshwater rates were well supported. High diversification rate suggest that some population dynamics may be affected by the ecological complexities of diadromy and that diadromy may have some other evolutionary significance. The high diversification rate is greatly impacted by a particularly low extinction rate in diadromy.
Model adequacy testing
Simulations of 500 trees that used the output parameters of the MuSSE model resulted in a distribution of trait frequencies such that the observed values fell within one standard deviation of the mean of the simulated distribution (Fig. 4).
Substantial transition rates between marine and freshwater in the best fit model do not support the hypothesis that all transitions between marine and freshwater move through diadromy. The phylogeny includes many recent transitions, decreasing the likelihood that there is an unobserved diadromous ancestor required in the transition. The patterns observed are neither in support of equal rates across all transitions (Fig. 2II) nor diadromy as an intermediate with one-to-one transition rates from marine to diadromy to freshwater or the other way around (Fig. 2III). The highest transition rates in both the fitDiscrete and MuSSE analyses were those out of diadromy to both marine and freshwater, which is in part responsible for diadromy having the lowest stability (SD = 0.51). Additionally, transitions from marine to diadromy were observed at an intermediate rate and those from freshwater to diadromy and between marine and freshwater showed the lowest rates. While transitions between marine and freshwater may move through diadromy, the evidence does not support labeling diadromy as a necessary transitional state. The hypothesis that even though not all transitions between marine and freshwater go through diadromy, and regardless of the directionality of transition in and out of diadromy, diadromy only acts as a transient stage is also not supported in light of high diversification rate observed in diadromy. This is a function of both high speciation and low extinction rates in diadromy. Instead, an evolutionary significance of diadromy may be that it acts as a potential source of diversity by allowing a lineage to enter a new habitat. This habitat transition may select for behavioral, anatomical or physiological variants that allow lineages to be successful in the new habitat or revert back to the ancestral habitat but with some variation from close relatives.
Many mechanisms could be responsible for the high rates of transition exiting diadromy including landlocked and oceanlocked populations. Freshwater populations can become isolated each time a set of random individuals becomes landlocked due to factors such as: a receding glacier forming isolated lakes, change in flow patterns of rivers, or the formation of natural dams due to tectonic processes. Oceanlocked species can result from a loss of suitable freshwater habitat in a species whose populations are geographically limited to ephemeral island streams. These high transition rates out of diadromy are illustrated by many examples of diadromous species with non-diadromous individuals, populations, or sister species [6, 18, 19, 37]. Species in which ocean- and landlocking has resulted in genetic isolation include alewives, galaxiids, gobies, sticklebacks, and lampreys [6, 18, 38]. With multiple ocean- or landlocking events and a wide enough distribution, multiple completely freshwater species, or completely marine species in the case of eels [19, 39], can emerge from a single diadromous species. Once species have diverged, fine-scale niche preferences can support a diversified clade in their newly occupied habitats . The directionality of observed transitions suggests that diadromy may sometimes be an intermediate state between marine and freshwater life histories and sometimes a temporary state in which lineages persist before returning to their ancestral state. Many papers address the causes and mechanisms of these transitions [5, 8, 9, 20]; however, those are beyond the scope of this study.
The high net-diversification rate of diadromous fishes is somewhat unexpected due to the potential for dispersal (especially in species without a strong tendency to return to their natal sites), since these fishes are capable of movement between distant freshwater sources during their marine phase. A wide range of dispersal distances have been reported in diadromous fishes [26, 41, 42]. Comparisons in some groups, such as minnows, show that marine dispersal is associated with larger range  and in lampreys, dispersal decreases the isolation by distance over large areas . The variation in dispersal time and distance is important in speciation and extinction rates. Short dispersal potential can result in a high level of isolation while still allowing for occasional propagules to venture out and form new populations with low levels of gene flow, a classic scenario for geographic speciation . Long distance dispersal can result in less structured populations compared to poorly dispersing diadromous fishes but could increase the overall range of the species allowing for more opportunities for isolation and transitions out of diadromy, especially in species that exhibit facultative diadromy. These results contradict expectations that diadromy has an intermediate diversification rate associated with an intermediate level of gene flow and range size compared to marine and freshwater fishes. Instead, diadromous fishes have higher diversification rates than either freshwater or marine lineages. The estimated speciation rate is about twice that of freshwater species and four times that of marine species (Table 1). This might arise if diadromous species, rather than having intermediate population structure, experience the best of both worlds for speciation. That is, perhaps diadromous species tend to have large geographic ranges (more like marine species on average), but also high levels of subdivision and local adaptation (more like freshwater species on average). Even more striking, the estimated extinction rate is several orders of magnitude lower for diadromous vs non-diadromous lineages (Table 1). Taken together with the high transition rates, it appears that diadromous ancestors might often escape extinction by leaving relictual freshwater populations in the case of unfavorable conditions in the ocean, and likewise by persisting in the ocean when freshwater populations are extirpated.
Observable patterns of diadromy
While the results presented here neither support nor refute the proposed hypothesis with respect to the causality of the evolution of diadromy, they do not conflict with the possibility of the transition patterns proposed by Gross et al. [5, 8] and Dodson et al. . These results support a hypothesis that a high rate of transitions out of diadromy coupled with the high net diversification rate in diadromy leads to a single species becoming diadromous, diversifying and then many species returning to non-diadromous lives. To date, many existing phylogenetic studies have discussed particular clades that fit the pattern presented by the MuSSE analysis, but not the transitional state hypotheses. Using the flagtails as an example, Feutry et al.  show a marine ancestry that transitioned into diadromy, and once diadromy evolved, multiple speciation events occurred, giving rise to six extant diadromous species of Kuhlia. Then, presumably via the gradual transition of the partially diadromous K. munda, a transition back to a completely marine life history occurred. Of the diadromous species, half of them (i.e., K. malo, K. rupestris and K. munda) exhibit intermediate phenotypic variation and/or partially catadromous life histories. While extant species show only one well defined transition out of diadromy, the loss of total obligate diadromy in some of the diadromous species could be an additional transition out of diadromy. As pointed out by Feutry et al. , the idea proposed by Gross et al. [5, 8] that marine species in tropical regions utilize freshwater habitats because of the increased resource availability does not apply to small island streams that the Kuhlia typically inhabit, as they tend to be nutrient poor.
In salmonids, Alexandrou et al.  show a transition from freshwater to diadromy and back to freshwater within Coregoninae. This transition into diadromy has resulted in eight diadromous species. Nested within the diadromus Coregonus species are three freshwater species (C. nigripinnis, C. zenithicus, and C. hoyi). Nested in those three species is another diadromous species (C. artedi). Either Coregonus artedi was a second transition into diadromy after a transition back to freshwater, or transition from diadromy to freshwater occurred twice in C. nigripinnis and in the C. zenithicus and C. hoyi clade. However, inconsistencies in the placement of various taxa in Coregonus [45, 46] could drastically change the order of events with regards to transitions and any assumptions about the evolution of diadromy could be incorrect. If we are to assume the Crete-Lafreniere et al.  phylogeny is correct, additional transitions out of diadromy are observed: 1) The Ireland endemic C. pollan nested in the widespread C. autumnalis, and 2) A group of multiple species sharing a common ancestor with the polyphyletic C. artedi which inhabit the freshwater of the US Great Lakes species. Some of these transitions (e.g., freshwater to diadromy in temperate regions) do agree with Gross and some do not (e.g. diadromy to freshwater in temperate regions) [5, 8]. In addition to the Coregoninae clade, Alexandrou et al.  show two additional transitions out of diadromy within salmonids in the genera Salmo and Salvelinus. Both Alexandrou et al. and Crete-Lafreniere et al. produced phylogenies with many diadromous species arising after a transition from freshwater, and at least one, if not multiple transitions from diadromy to freshwater. These patterns coincide with the parameters of the best fit MuSSE model.
Within the Anguilliformes, one major transition from marine to diadromy has occurred in the freshwater eels of the genus Anguilla resulting in the 16 species and three subspecies . In agreement with the prediction of Gross et al. [5, 8], we find marine ancestors transitioning to a diadromous life history in a tropical region. There is evidence that in tropical regions individuals of Anguilla marmorata , and in temperate regions A. japonica and A. anguilla , have reverted back to completely marine life histories. Although whether these individuals reproduce or are merely lost at sea remains unknown, the large scale marine ➔ diadromy ➔ marine scenario does not fit the hypotheses proposed by Gross et al. [5, 8] and Dodson et al. .
Many landlocking events have been observed in the threespine stickleback (Gasterosteus aculeatus), for which multiple mechanisms have been proposed [17, 49]. Despite the genetic and life-history variation between diadromous and freshwater G. aculeatus populations, hybridization does occur , though no studies report gene flow between distinct, geographically isolated freshwater populations (although introgressive hybridization between sympatric benthic-limnetic species pairs has occurred ). It is possible that with current population dynamics, each landlocked population could eventually lead to distinct taxa, since adaptive radiation has allowed morphological and behavioral variation to accrue in landlocked populations within a short time. Unlike the sympatric sticklebacks, diadromous and landlocked alewife (Alosa pseudoharengus) do not typically co-occur, and genetic evidence indicates that many of the landlocked populations are the result of separate isolation events such as damming of a river . Additionally, there are morphological, behavioral and genetic differences between landlocked and diadromous populations that have diverged within the last century . Studies have even started to pinpoint specific functional genes that may play a role in the ability to survive landlocking events [52, 53]. Though no speciation has been recognized in A. pseudoharengus, these independent landlocking events (similar to sticklebacks) are the result of multiple transitions out of diadromy .
Complex life histories that require environmental transitions can serve as much more than simple evolutionary stepping stones between different life histories. Diadromy provides species with the ability to enter new environments and diversify, even if the behavior may quickly be lost. While diadromy can act as an evolutionary intermediate between marine and freshwater life histories, it also plays an important role in diversification as well as habitat transitions that are not necessarily between marine and freshwater, but instead marine back to marine or freshwater back to freshwater after diversification from a diadromous life history. With this in mind, continuous studies at the individual, species, and family level are important to understand the mechanism of how and why diadromy is gained and lost. Applying SSE models to other complex life histories may reveal more patterns of biodiversity. Additionally, considering complex life histories that require environmental transitions not only in terms of their costly transitions, but for the advantages of using multiple environments, may allow for new perspectives on multiple topics related to life history evolution.
The 7822 species of ray-finned fishes (Actinopterygii) used in the recently published phylogeny  were each assigned one of three character states: 1) diadromous (n = 180), 2) freshwater (n = 3998), and 3) marine (n = 3600). All species that could not be assigned were classified as “unknown” and dropped from the tree (n = 45) (Additional file 1) [6, 54,55,56]. Diadromous fishes were identified from 28 families (Salmonidae: n = 32, Gobiidae: n = 19, Anguillidae: n = 18, Clupeidae: n = 17, Acipenseridae: n = 16, Eleotridae: n = 14, Osmeridae: n = 9, Galaxiidae: n = 7, Cottidae: n = 6, Achiridae: n = 5, Mugilidae: n = 5, Cyprinidae: n = 4, Retropinnidae: n = 4, Salangidae: n = 4, Ariidae: n = 2, Gasterosteidae: n = 2, Moronidae: n = 2, Pleuronectidae: n = 2, Rhyacichthidae: n = 2, Tetraodontidae: n = 2, Cheimarrichthyidae: n = 1, Gadidae: n = 1, Kuhliidae; n = 1, Lateolabracidae: n = 1, Latidae: n = 1, Plecoglossidae: n = 1, Pseudaphritidae: n = 1, Syngnathidae: n = 1). Character states were determined primarily based on information in: Diadromy in Fishes , Fishbase , Fishes of the World , and Encyclopedia of Life . Species considered as euryhaline, or those that can tolerate a large range of salinity including freshwater, were assigned character states based on their dominant life strategy. Bairdiella chrysoura, for example, is listed as marine, freshwater, and brackish in Fishbase, but since the species is described as moving “to the nursery and feeding areas in estuaries during summer months and sometimes enters freshwaters,” , B. chrysoura was listed as a marine species. When inadequate information was available on the aforementioned sources, other primary literature was used in considering salinity levels during breeding, distribution, and salt tolerance, to better assign species to each of the three categories (additional sources used for character assignment: [23, 24, 57,58,59,60,61,62,63,64,65,66,67,68,69,70,71,72,73,74,75,76]). Additionally, the species listed as both marine and freshwater were confirmed via peer-reviewed evidence, such as in the case of Takifugu poecilonotus, where studies confirmed diadromy . Species such as Alosa pseudoharengus , and many of the freshwater eels  with both diadromous and non-diadromous individuals/populations, were characterized as diadromous since the behavior is dominant within these species. This was preferable to removing or duplicating these species which would alter measures of speciation and transition between the states.
Phylogenetic model selection
Species were assigned character states F (freshwater), M (marine), or D (diadromous) for three-state MuSSE analysis. MuSSE gives estimates for speciation (λ), extinction (μ), and transition between states (qFD; in this example from state F to state D). Net diversification rates (r) were calculated by subtracting extinction rates from speciation rates (r = λ-μ). The MuSSE analysis estimates 12 parameters: λF, λM, λD, μF, μM, μD, qFM, qFD, qMD, qMF, qDF, and qDM. Testing different hypotheses was done by putting constraints on specific parameters. For example, to test the hypothesis that all transitions between marine and freshwater have to go through diadromy as shown in Fig. 2I, constraints were put on transition between marine and freshwater in both directions (qFM = qMF = 0) and on speciation or extinction in diadromy (λD = μD = 0). All remaining parameters (λF, λM, μF, μM, qFD, qMD, qDF, and qDM) were estimated by MuSSE. Looking at Fig. 2III, constraints were put on transitions making estimates of marine to diadromy and diadromy to freshwater equal, as well as freshwater to diadromy and diadromy to marine (qMD = qDF and qFD = qDM) with all remaining state free to vary. State Stability (S) was calculated as one minus the probability of extinction of that state minus the probability of transition out of that state to all other states (Si = 1– μi – qij – qik). Similar to other SSE models (e.g., BiSSE, HiSSE, GeoSSE), rates (λ, μ and q) are each estimated independently using likelihood starting at the tips of the tree and working back in time along each branch in small time increments. Estimates for each parameter are compiled at nodes until likelihoods have been calculated to the root. For a full description of parameter estimation see Maddison et al. , FitzJohn  and Ng and Smith .
MuSSE analyses were run in R  using the “diversitree” package . The correction for incomplete sampling in the tree, the sampling.f = (.26, .21, .46) function, was used to inform the program of the percent of freshwater (3998/15,170) and marine (3600/16,764) species based on published estimates . Because recognized non-diadromous species occasionally exhibit diadromous behavior, the true number of species that are diadromous remains unknown. The percent of diadromous species used (104/223) was calculated based on the inventory provided in “Diadromy in Fishes,”  where only 104 of the 223 species listed were present in the phylogeny. It is important to note that these percentages are representative of the complete phylogeny across each state starting at the root node. There may be some groups that are fully represented and some that are not represented at all. This may be especially pertinent to radiations such as Sicydiinae gobies or cichlids that have lower than average representation compared to Salmonids that have higher than average representation. A total of 27 models were compared using the MuSSE model (Additional file 2). Akaike information criterion (AIC) scores were calculated from the log likelihood value outputs from each model. ∆AIC was calculated by subtracting the lowest AIC value from each AIC score. Akaike weights (wi) for each model was calculated. Bayesian estimates of the parameters using MCMC, as implemented with the mcmc function in DIVERSITREE, were obtained with runs of 10,000 steps. Additionally, to test robustness of character state assignment, models were run with all 45 “unknown” species placed completely in marine or completely in freshwater. Results presented in Table 1 and Additional file 1 do not include these species. All models were run using the single best fit tree reported by Rabosky et al. . No analyses using suboptimal, yet highly supported trees with variation in tree topology or branch lengths were used because no additional topologies were made available for analysis. Additional tree topologies and variation in branch lengths could change the estimated rates produced by all models.
Estimates of extinction rates using phylogenies have been suggested to be inaccurate , however, studies have shown that in certain situations phylogenies can be good tools for extinction estimates [41, 84]. As pointed out by Beaulieu and O’Meara , using a tree with a large number of taxa (> 100) should be adequate for rate estimations involving extinction, as long as the variation in birth rate is not exceptionally high (also see Rabosky ). To see if varying extinction rates affected other parameters, models were compared with all extinction rates constrained to be equal (μF ~ μM ~ μD). Parameter estimates from this model resulted in speciation and transition rates similar to those observed in the best fit models (Table 1).
In some cases, an additional characteristic might be nested within the trait of interest resulting in a skewed rate estimation. For example, if temperate diadromous species have a higher speciation rate then tropical diadromous species, or if particular groups or types of diadromy differ from the rest of the diadromous fishes, then they can skew the rate for the diadromous group as a whole when tested in the traditional BiSSE and MuSSE models. Testing for a hidden state can show that other characteristics of the species are important in understanding the evolutionary rates (for a more detailed discussion see Beaulieu and O’Meara ). To test for a “hidden state” in diadromy, Hidden State Speciation Extinction models (HiSSE)  were also performed. HiSSE assigns a likelihood of each species within a state to a sub state (e.g., “diadromy-observed” and “diadromy-hidden”). HiSSE is a two-state model and was used to test for hidden rates within diadromous and non-diadromous fishes. HiSSE models without hidden states are identical to Binary State Speciation Extinction models (BiSSE) and were also run. While the two state models do not allow for testing the hypotheses of whether diadromy is an intermediate state, identifying a potential hidden state may reveal that a few select taxa are skewing the rates reported in the MuSSE analysis. Using the MarginRecon function with the best fit model parameters, all diadromous species were assigned to the diadromy-observed state with > 98% likelihood, with < 2% likelihood of being in a diadromy-hidden state, suggesting no hidden state in diadromy. As a result, models were run with no hidden state for diadromous fishes.
To support the transition rate estimation reported in MuSSE, fitDiscrete function in geiger 2.0.3  was also used to estimate transition rates in both a two state (diadromous / nondiadromous) and three state (diadromous / marine / freshwater) data set. Models were compared using AIC. For both two and three-state situations, equal-rate (ER), symmetric (SYM), all rates different (ARD), and meristic (meristic) models were tested. The ER model constrains all transitions rates to be equal (qFM = qFD = qMD = qMF = qDF = qDM), SYM constrains rates between two states to be equal (qFM = qMF) but they are free to vary between different states ([qFM = qMF] ≠ [qFD = qDF] ≠ [qDM = qMD]), ARD allows each rate to vary as separate parameter (qFM ≠ qFD ≠ qMD ≠ qMF ≠ qDF ≠ qDM), and meristic requires transitions to occur in a stepwise fashion (i.e., M➔D➔F➔D➔M, but not M➔F).
Model adequacy testing
To test for model adequacy, trees were simulated with the tree.musse function in the diversitree package  incorporating the 12 MuSSE output parameters from the best fit model, with the ancestral state set as marine, 500 trees were simulated to 7778 tips. All trees were simulated with each state being assigned a single estimation of λ and μ that remained consistent throughout the tree. This does not allow the simulation to incorporate mass extinctions or adaptive radiations that have been recorded in the true phylogeny of fishes, and can result in variation of λ and μ at particular locations of the tree. The distribution of the percent of marine, freshwater and diadromous fishes from the 500 simulations were compared to the original percent of taxa in each group. These simulations serve two purposes. First, they validate the MuSSE model by showing the output parameters simulate what is observed in nature. Secondly, they address the potential bias of uneven sampling in the phylogeny used in the study. Showing that results from the tree are representative of what is observed in nature implies that the tree adequately represents the true tree of fishes.
Availability of data and materials
All data generated or analyzed during this study are included in this published article and its supplementary information files.
Akaike information criterion
All rates different
Binary State Speciation and Extinction
Hidden State Speciation and Extinction
Multiple State Speciation and Extinction
Shubin NH, Daeschler EB, Jenkins FA. The pectoral fin of Tiktaalik roseae and the origin of the tetrapod limb. Nature. 2006;440(7085):764–71.
Dudley R, Byrnes G, Yanoviak SP, Borrell B, Brown RM, McGuire JA. Gliding and the functional origins of flight: biomechanical novelty or necessity? Annu Rev Ecol Evol Syst. 2007;38:179–201.
Bloom DD, Lovejoy NR. The evolutionary origins of diadromy inferred from a time-calibrated phylogeny for Clupeiformes (herring and allies). Proc Biol Sci. 2014;281(1778):2013–81.
Myers GS. Usage of anadromous, catadromous and allied terms for migratory fishes. Copeia. 1949;2:89–97.
Gross MR. Evolution of diadromy in fishes. In: Amer Fish Soc Symp, vol. 1; 1987. p. 14–25.
McDowall RM. Diadromy in fishes: migrations between freshwater and marine environments. Portland: Timber Press; 1988.
Hodges MH, Allendorf FW. Population genetics and pattern of larval dispersal of the endemic Hawaiian freshwater amphidromous gastropod Neritina granosa (Prosobranchia: Neritidae). Pac Sci. 1998;52(3):237–49.
Gross MR, Coleman RM, McDowall RM. Aquatic productivity and the evolution of diadromous fish migration. Science. 1988;239(4845):1291–3.
Dodson JJ, Laroche J, Lecomte F. Contrasting evolutionary pathways of anadromy in euteleostean fishes. In: Challenges for diadromous fishes in a dynamic global environment. Am Fish S S, vol. 69; 2009. p. 63–77.
McCormick SD, Anthony PF, Colin JB, editors. Fish physiology: euryhaline, vol. 32. Waltham: Academic Press; 2013.
Black VS. Excretion and osmoregulation. Physiol Fishes. 1957;1:163–205.
Rabosky DL, Santini F, Eastman J, Smith SA, Sidlauskas B, Chang J, et al. Data from: rates of speciation and morphological evolution are correlated across the largest vertebrate radiation: Dryad Digital Repository; 2013b. https://doi.org/10.5061/dryad.j4802.
Castelin M, Feutry P, Hautecoeur M, Marquet G, Wowor D, Zimmermann G, et al. New insight on population genetic connectivity of widespread amphidromous prawn Macrobrachium lar (Fabricius, 1798)(Crustacea: Decapoda: Palaemonidae). Mar Biol. 2013;160(6):1395–406.
Kuribayashi K, Katakura H, Kyono M, Dick MH, Mawatari SF. Round-trip catadromous migration in a Japanese amphipod, Sternomoera rhyaca (Gammaridea: Eusiridae). Zool Sci. 2006;23(9):763–74.
Chen LH, Yang JL, Barnes MD. Effect of natural disturbance on migrations and population dynamics of a catadromous species, Eriocheir japonicus, along the Lo-Mei stream in northern Taiwan. Hydrobiologia. 2014;726(1):179–94.
Parenti LR. Life history patterns and biogeography: an interpretation of diadromy in fishes. Ann Mo Bot Gard. 2008;95(2):232–47.
McPhail JD. Speciation and the evolution of reproductive isolation in the sticklebacks (Gasterosteus) of south-western British Columbia. In: The evolutionary biology of the threespine stickleback; 1994. p. 399–437.
Palkovacs EP, Kirstin BD, Post DM, Caccone A. Independent evolutionary origins of landlocked alewife populations and rapid parallel evolution of phenotypic traits. Mol Ecol. 2008;17(2):582–97.
Tsukamoto K, Izumi N. Do all freshwater eels migrate? Nature. 1998;396(6712):635–6.
Feutry P, Castelin M, Ovenden JR, Dettaï A, Robinet T, Cruaud C, Keith P. Evolution of diadromy in fish: insights from a tropical genus (Kuhlia species). Am Nat. 2013;181(1):52–63.
Lee CE, Bell MA. Causes and consequences of recent freshwater invasions by saltwater animals. Trends Ecol Evol. 1999;14(7):284–8.
Betancur-R R. Molecular phylogenetics supports multiple evolutionary transitions from marine to freshwater habitats in ariid catfishes. Mol Phylogenet Evol. 2010;55:249–58.
Betancur-R R, Ortí G, Pyron RA. Fossil-based comparative analyses reveal ancient marine ancestry erased by extinction in ray-finned fishes. Ecol Lett. 2015;18(5):441–50.
Bloom DD, Lovejoy NR. On the origins of marine-derived freshwater fishes in South America. J Biogeogr. 2017;44(9):1927–38.
McDowall RM. Diadromy, diversity and divergence: implications for speciation processes in fishes. Fish Fish. 2001;2(3):278–85.
Shiao JC, Tzeng CS, Li PC, Bell KN. Upstream migration and marine early life history of amphidromous gobies inferred from otolith increments and microchemistry. Environ Biol Fish. 2015;98(3):933–50.
D'Aloia CC, Bogdanowicz SM, Harrison RG, Buston PM. Seascape continuity plays an important role in determining patterns of spatial genetic structure in a coral reef fish. Mol Ecol. 2014;23(12):2902–13.
Alfaro ME, Santini F, Brock CD. Do reefs drive diversification in marine teleosts? Evidence from the pufferfish and their allies (Order Tetraodontiformes). Evolution. 2007;61(9):2104–26.
Taylor MS, Hellberg ME. Comparative phylogeography in a genus of coral reef fishes: biogeographic and genetic concordance in the Caribbean. Mol Ecol. 2006;15(3):695–707.
Ward RD, Woodwark M, Skibinski DOF. A comparison of genetic diversity levels in marine, freshwater, and anadromous fishes. J Fish Biol. 1994;44(2):213–32.
Bloom DD, Weir JT, Piller KR, Lovejoy NR. Do freshwater fishes diversify faster than marine fishes? A test using state-dependent diversification analyses and molecular phylogenetics of new world silversides (atherinopsidae). Evolution. 2013;67(7):2040–57.
McDowall RM. The evolution of diadromy in fishes (revisited) and its place in phylogenetic analysis. Rev Fish Biol Fisher. 1997;7(4):443–62.
Betancur-R R, Ortí G, Stein AM, Marceniuk AP, Pyron A. R. Apparent signal of competition limiting diversification after ecological transitions from marine to freshwater habitats. Ecol Lett. 2012;15(8):822–30.
Watanabe K, Sakai H, Sanada T, Nishida M. Comparative phylogeography of diadromous and freshwater daces of the genus Tribolodon (Cyprinidae). Ichthyol Res. 2018;1:1–15.
Rabosky DL, Santini F, Eastman J, Smith SA, Sidlauskas B, Chang J, et al. Rates of speciation and morphological evolution are correlated across the largest vertebrate radiation. Nat Commun. 2013a;4:1958.
FitzJohn RG. Diversitree: comparative phylogenetic analyses of diversification in R. Method Ecol Evol. 2012;3(6):1084–92.
Augspurger JM, Warburton M, Closs GP. Life-history plasticity in amphidromous and catadromous fishes: a continuum of strategies. Rev Fish Biol Fisher. 2017;27:1–16.
Bracken FS, Hoelzel AR, Hume JB, Lucas MC. Contrasting population genetic structure among freshwater-resident and anadromous lampreys: the role of demographic history, differential dispersal and anthropogenic barriers to movement. Mol Ecol. 2015;24(6):1188–204.
Daverat F, Limburg KE, Thibault I, Shiao JC, Dodson JJ, Caron F, et al. Phenotypic plasticity of habitat use by three temperate eel species, Anguilla anguilla, A. japonica and A. rostrata. Mar Ecol Prog Ser. 2006;308:231–41.
Glova GJ, Jellyman DJ, Bonnett ML. Factors associated with the distribution and habitat of eels (Anguilla spp.) in three New Zealand lowland streams. NZ J Mar Freshw Res. 1998;32:255–69.
Nee S, Holmes EC, May RM, Harvey PH. Extinction rates can be estimated from molecular phylogenies. Philos T Roy Soc B. 1994;344(1307):77–82.
Feutry P, Valade P, Ovenden JR, Lopez PJ, Keith P. Pelagic larval duration of two diadromous species of Kuhliidae (Teleostei: Percoidei) from Indo-Pacific insular systems. Mar Freshw Res. 2012;63(5):397–402.
Coyne JA, Orr HA. Speciation. Sunderland: Sinauer Association; 2004.
Alexandrou MA, Swartz BA, Matzke NJ, Oakley TH. Genome duplication and multiple evolutionary origins of complex migratory behavior in Salmonidae. Mol Phylogenet Evol. 2013;69(3):514–23.
Crête-Lafrenière A, Weir LK, Bernatchez L. Framing the Salmonidae family phylogenetic portrait: a more complete picture from increased taxon sampling. PLoS One. 2012;7(10):e46662.
Horreo JL. Revisiting the mitogenomic phylogeny of Salmoninae: new insights thanks to recent sequencing advances. PeerJ. 2017;5:e3828.
Inoue JG, Miya M, Miller MJ, Sado T, Hanel R, Hatooka K, et al. Deep-ocean origin of the freshwater eels. Biol Lett. 2010;6(3):363–6.
Arai T, Chino N. Opportunistic migration and habitat use of the giant mottled eel Anguilla marmorata (Teleostei: Elopomorpha). Sci Rep-UK. 2018;8(1):5666.
Bell MA, Andrews CA. Evolutionary consequences of colonization of fresh water by primitively anadromous fishes. In: Streit B, Stadler T, Lively CM, editors. Evolutionary ecology of freshwater animals. Current trends and case studies. Basel: Birkhauser; 1997. p. 323–63.
Velema GJ, Rosenfeld JS, Taylor EB. Effects of invasive American signal crayfish (Pacifastacus leniusculus) on the reproductive behaviour of threespine stickleback (Gasterosteus aculeatus) sympatric species pairs. Can J Zoo. 2012;90(11):1328–38.
Czesny S, Epifanio J, Michalak P. Genetic divergence between freshwater and marine morphs of alewife (Alosa pseudoharengus): a ‘next-generation’sequencing analysis. PLoS One. 2012;7(3):e31803.
Michalak K, Czesny S, Epifanio J, Snyder RJ, Schultz ET, Velotta JP, et al. Beta-thymosin gene polymorphism associated with freshwater invasiveness of alewife (Alosa pseudoharengus). J Exp Zoo Part A. 2014;321(4):233–40.
Lemopoulos A, Uusi-Heikkilä S, Huusko A, Vasemägi A, Vainikka A. Comparison of migratory and resident populations of brown trout reveals candidate genes for migration tendency. Genome Biol Evol. 2018;10(6):1493–503.
Froese R, Pauly D. FishBase. World Wide Web electronic publication. 2014. www.fishbase.org.
Nelson JS. Fishes of the world. New York: Wiley; 2006.
Wilson EO. The encyclopedia of life. Trends Ecol Evol. 2003;18:77–80.
Frost WE. The natural history of the minnow, Phoxinus phoxinus. J Anim Ecol. 1943;12:139–62.
Breder CM. Modes of reproduction in fishes. Garden City: Natural History Press; 1966.
Foster NR. Behavior, development, and early life history of the Asian needlefish, Xenentodon cancila. P Acad Nat Sci Phila. 1973;125:77–88.
Govoni JJ, Merriner JV. The occurrence of ladyfish, Elops saurus, larvae in low salinity waters and another record for Chesapeake Bay. Estuaries. 1978;1(3):205–6.
Linden E, Bengtsson BE, Svanberg O, Sundström G, N. S. E. P. Board. The acute toxicity of 78 chemicals and pesticide formulations against two brackish water organisms, the bleak (Alburnus alburnus) and the harpacticoid Nitocra spinipes. Chemosphere. 1979;8(11–12):843–51.
Allen GR. Snappers of the world: an annotated and illustrated catalogue of lutjanid species known to date. Rome: Food and Agriculture Organization; 1985. FAO Fisheries Synopsis, no. 125, vol. 6
Goto A. Alternative life-history styles of Japanese freshwater sculpins revisited. Environ Bio Fish. 1990;28:101–12.
Marais JFK, Venter DJL. Changes in body composition associated with growth and reproduction in Galeichthys feliceps (Teleostei: Ariidae). S Afr J Marine Sci. 1991;10(1):149–57.
Rosecchi E, Crivelli AJ. Study of a sand smelt (Atherina boyeri Risso 1810) population reproducing in fresh water. Ecol Freshw Fish. 1992;1(2):77–85.
Sakai H. Life-histories and genetic divergence in three species of Tribolodon (Cyprinidae). Mem Fac Fish Hokkaido Univ. 1995;1(2):1–98.
Rainboth WJ. Fishes of the Cambodian Mekong: Food and Agriculture Org; 1996.
Garg SK. Brackishwater carp culture in potentially waterlogged areas using animal wastes as pond fertilizers. Aquacult Int. 1996;4(2):143–55.
Bemis WE, Findeis EK, Grande L. An overview of Acipenseriformes. Environ Biol Fish. 1997;48(1–4):25–71.
Natsumeda T. Year-round local movements of the Japanese fluvial sculpin, Cottus pollux (large egg type), with special reference to the distribution of spawning nests. Ichthyol Res. 1999;46(1):43–8.
Takeshita N, Onikura N, Nagata S, Matsui S, Kimura S. A note on the reproductive ecology of the catadromous fourspine sculpin, Cottus kazika (Scorpaeniformes: Cottidae). Ichthyol Res. 1999;46(3):309–13.
Lovejoy NR, Collette BB. Phylogenetic relationships of new world needlefishes (Teleostei: Belonidae) and the biogeography of transitions between marine and freshwater habitats. Copeia. 2001;2:324–38.
Chen YF. Induced ovulation and embryonic development of ocellated puffer, Takifugu ocellatus. J Appl Ichthyol. 2005;21(2):136–40.
Toyama M, Hironaka M, Yamahama Y, Horiguchi H, Tsukada O, Uto N, et al. Presence of rhodopsin and porphyropsin in the eyes of 164 fishes, representing marine, diadromous, coastal and freshwater species—a qualitative and comparative study. Photochem Photobiol. 2008;84(4):996–1002.
Hua X, Wang W, Yin W, He Q, Jin B, Li J, et al. Phylogeographical analysis of an estuarine fish, Salanx ariakensis (Osmeridae: Salanginae) in the north-western Pacific. J Fish Biol. 2009;75(2):354–67.
Machado I, Vera M, Calliari D, Rodríguez-Graña L. First record of an Elops smithi (Pisces: Elopidae) larva in a south American subtropical-temperate estuary. Mar Biodivers Records. 2012;5:1-5.
Kato A, Doi H, Nakada T, Sakai H, Hirose S. Takifugu obscurus is a euryhaline fugu species very close to Takifugu rubripes and suitable for studying osmoregulation. BMC Physiol. 2005;5(1):18.
Thibault I, Dodson JJ, Caron F, Tzeng WN, Iizuka Y, Shiao JC. Facultative catadromy in American eels: testing the conditional strategy hypothesis. Mar Ecol Prog Ser. 2007;344:219–29.
Maddison WP, Midford PE, Otto SP. Estimating a binary character’s effect on speciation and extinction. Syst Boil. 2007;56(5):701–10.
Ng J, Smith SD. How traits shape trees: new approaches for detecting character state-dependent lineage diversification. J Evolution Biol. 2014;10:2035–45.
R Core Team. R: a language and environment for statistical computing: R Foundation for Statistical Computing; 2013.
Eschmeyer WN, Fricke R, Fong JD, Polack DA. Marine fish diversity: history of knowledge and discovery (Pisces). Zootaxa. 2010;2525(1):19–50.
Rabosky DL. Extinction rates should not be estimated from molecular phylogenies. Evolution. 2010;64(6):1816–24.
Beaulieu JM, O'Meara BC. Extinction can be estimated from moderately sized molecular phylogenies. Evolution. 2015;69(4):1036–43.
Rabosky DL. Challenges in the estimation of extinction from molecular phylogenies: a response to Beaulieu and O’Meara. Evolution. 2016;70(1):218–28.
Beaulieu JM, O'Meara BC. Detecting hidden diversification shifts in models of trait-dependent speciation and extinction. Syst Boil. 2016;65(4):583–601.
Harmon LJ, Weir JT, Brock CD, Glor RE, Challenger W. GEIGER: investigating evolutionary radiations. Bioinformatics. 2008;24(1):129–31.
I would like to thank Benjamin M. Fitzpatrick for his discussion and feedback on all aspects of this project. I would like to thank Jeremy Beaulieu, Brian O’Meara, and Orlando Schwery for their help in SSE model analyses, Darrin Hulsey and Brent Tibbits for their discussions on diadromous fishes during the formation of this project, and Graciela Cabana, Benjamin Keck, and Danial Simberloff for their comments throughout the project as well as Glen Hood and four anonymous reviewers for their comments on the final manuscript. K. Lindsay Hunter copy edited drafts of the manuscript.
Funding for all aspects of this study was provided by the University of Tennessee and the University of Tennessee Department of Ecology and Evolutionary Biology in the form of graduate student support. Funding for open access to this research was provided by University of Tennessee’s Open Publishing Support Fund.
Ethics approval and consent to participate
Consent for publication
The author declares that he has no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Table of Species character state assignment. (XLS 832 kb)
Model constraints for all MuSSE analyses. (DOCX 18 kb)
Output parameters for optimal model in MuSSE analyses with ambiguous species assigned to marine and freshwater. (DOCX 16 kb)
Models constraints for HiSSE and BiSSE analysis. (DOCX 24 kb)