rfishbase: exploring, manipulating and visualizing FishBase data from R

Creating functional groups of marine fish from categorical traits

Background

Functional groups serve two important functions in ecology: they allow for simplification of ecosystem models and can aid in understanding diversity. Despite their important applications, there has not been a universally accepted method of how to define them. A common approach is to cluster species on a set of traits, validated through visual confirmation of resulting groups based primarily on expert opinion. The goal of this research is to determine a suitable procedure for creating and evaluating functional groups that arise from clustering nominal traits.

Methods

To do so, we produced a species by trait matrix of 22 traits from 116 fish species from Tasman Bay and Golden Bay, New Zealand. Data collected from photographs and published literature were predominantly nominal, and a small number of continuous traits were discretized. Some data were missing, so the benefit of imputing data was assessed using four approaches on data with known missing values. Hierarchical clustering is utilised to search for underlying data structure in the data that may represent functional groups. Within this clustering paradigm there are a number of distance matrices and linkage methods available, several combinations of which we test. The resulting clusters are evaluated using internal metrics developed specifically for nominal clustering. This revealed the choice of number of clusters, distance matrix and linkage method greatly affected the overall within- and between- cluster variability. We visualise the clustering in two dimensions and the stability of clusters is assessed through bootstrapping.

Results

Missing data imputation showed up to 90% accuracy using polytomous imputation, so was used to impute the real missing data. A division of the species information into three functional groups was the most separated, compact and stable result. Increasing the number of clusters increased the inconsistency of group membership, and selection of the appropriate distance matrix and linkage method improved the fit.

Discussion

We show that the commonly used methodologies used for the creation of functional groups are fraught with subjectivity, ultimately causing significant variation in the composition of resulting groups. Depending on the research goal dictates the appropriate strategy for selecting number of groups, distance matrix and clustering algorithm combination.

Drivers and rates of stock assessments in the United States

Fisheries management is most effective when based on scientific estimates of sustainable fishing rates. While some simple approaches allow estimation of harvest limits, more data-intensive stock assessments are generally required to evaluate the stock's biomass and fishing rates relative to sustainable levels. Here we evaluate how stock characteristics relate to the rate of new assessments in the United States. Using a statistical model based on time-to-event analysis and 569 coastal marine fish and invertebrate stocks landed in commercial fisheries, we quantify the impact of region, habitat, life-history, and economic factors on the annual probability of being assessed. Although the majority of landings come from assessed stocks in all regions, less than half of the regionally-landed species currently have been assessed. As expected, our time-to-event model identified landed tonnage and ex-vessel price as the dominant factors determining increased rates of new assessments. However, we also found that after controlling for landings and price, there has been a consistent bias towards assessing larger-bodied species. A number of vulnerable groups such as rockfishes (Scorpaeniformes) and groundsharks (Carcharhiniformes) have a relatively high annual probability of being assessed after controlling for their relatively small tonnage and low price. Due to relatively low landed tonnage and price of species that are currently unassessed, our model suggests that the number of assessed stocks will increase more slowly in future decades.

Pulsed evolution shaped modern vertebrate body sizes

The diversity of forms found among animals on Earth is striking. Despite decades of study, it has been difficult to reconcile the patterns of diversity seen between closely related species with those observed when studying single species on ecological timescales. We propose a set of models, called Lévy processes, to attempt to reconcile rapid evolution between species with the relatively stable distributions of phenotypes seen within species. These models, which have been successfully used to model stock market data, allow for long periods of stasis followed by bursts of rapid change. We find that many vertebrate groups are well fitted by Lévy models compared with models for which traits evolve toward a stationary optimum or evolve in an incremental and wandering manner.

The relative importance of different modes of evolution in shaping phenotypic diversity remains a hotly debated question. Fossil data suggest that stasis may be a common mode of evolution, while modern data suggest some lineages experience very fast rates of evolution. One way to reconcile these observations is to imagine that evolution proceeds in pulses, rather than in increments, on geological timescales. To test this hypothesis, we developed a maximum-likelihood framework for fitting Lévy processes to comparative morphological data. This class of stochastic processes includes both an incremental and a pulsed component. We found that a plurality of modern vertebrate clades examined are best fitted by pulsed processes over models of incremental change, stationarity, and adaptive radiation. When we compare our results to theoretical expectations of the rate and speed of regime shifts for models that detail fitness landscape dynamics, we find that our quantitative results are broadly compatible with both microevolutionary models and observations from the fossil record.

A global database on freshwater fish species occurrence in drainage basins

A growing interest is devoted to global-scale approaches in ecology and evolution that examine patterns and determinants of species diversity and the threats resulting from global change. These analyses obviously require global datasets of species distribution. Freshwater systems house a disproportionately high fraction of the global fish diversity considering the small proportion of the earth's surface that they occupy, and are one of the most threatened habitats on Earth. Here we provide complete species lists for 3119 drainage basins covering more than 80% of the Earth surface using 14953 fish species inhabiting permanently or occasionally freshwater systems. The database results from an extensive survey of native and non-native freshwater fish species distribution based on 1436 published papers, books, grey literature and web-based sources. Alone or in combination with further datasets on species biological and ecological characteristics and their evolutionary history, this database represents a highly valuable source of information for further studies on freshwater macroecology, macroevolution, biogeography and conservation.

Improving understanding of the functional diversity of fisheries by exploring the influence of global catch reconstruction

Functional diversity is thought to enhance ecosystem resilience, driving research focused on trends in the functional composition of fisheries, most recently with new reconstructions of global catch data. However, there is currently little understanding of how accounting for unreported catches (e.g. small-scale and illegal fisheries, bycatch and discards) influences functional diversity trends in global fisheries. We explored how diversity estimates varied among reported and unreported components of catch in 2010, and found these components had distinct functional fingerprints. Incorporating unreported catches had little impact on global-scale functional diversity patterns. However, at smaller, management-relevant scales, the effects of incorporating unreported catches were large (changes in functional diversity of up to 46%). Our results suggest there is greater uncertainty about the risks to ecosystem integrity and resilience from current fishing patterns than previously recognized. We provide recommendations and suggest a research agenda to improve future assessments of functional diversity of global fisheries.

Integrating lipid storage into general representations of fish energetics

Fish, even of the same species, can exhibit substantial variation in energy density (energy per unit wet weight). Most of this variation is due to differences in the amount of storage lipids. In addition to their importance as energy reserves for reproduction and for survival during unfavourable conditions, the accumulation of lipids represents a large energetic flux for many species, so figuring out how this energy flux is integrated with other major energy fluxes (growth, reproduction) is critical for any general theory of organismal energetics. Here, we synthesize data from a wide range of fish species and identify patterns of intraspecific variation in energy storage, and use these patterns to formulate a general model of energy allocation between growth, lipid storage and reproduction in fishes. From the compiled data we identified two patterns: (1) energy density increases with body size during the juvenile period, but is invariant with body size within the adult size range for most species, and (2) energy density changes across seasons, with depletion over winter, but increases fastest in periods of transition between favourable and unfavourable conditions for growth (i.e. fall). Based on these patterns we propose DEBlipid, a simple, general model of energy allocation that is closely related to a simplified version of Dynamic Energy Budget theory, DEBkiss. The crux of the model is that assimilated energy is partitioned, with κ fraction of energy allocated to pay maintenance costs first, and the surplus allocated to growth, and 1 - κ fraction of assimilated energy is allocated to accumulating storage lipids during the juvenile phase, and later to reproduction as adults. This mechanism, in addition to capturing the two patterns that motivated the model, was able to predict lipid dynamics in a novel context, the migration of anadromous fish from low-food freshwater to high-food marine environments. Furthermore, the model was used to explain intra and interspecific variation in reproductive output based on patterns of lipid accumulation as juveniles. Our results suggest that many seemingly complex, adaptive energy allocation strategies in response to ontogeny, seasonality and habitat quality can emerge from a simple physiological heuristic.

Sex-related mercury bioaccumulation in fish from the Madeira River, Amazon

Sex plays an important role in the kinetics and dynamics of methylmercury in some animals. Although fish is the main source of mercury exposure to consumers, the role of sex in fish-Hg bioaccumulation is less known. We studied total Hg (THg) concentrations in 2538 samples (males=1052, females=1486) of fish from different trophic levels (herbivorous, planctivorous, detritivorous, omnivorous, carnivorous, piscivorous); for each species we made a post hoc estimation of the minimum number of samples required to detect variance-based differences between sexes. Only five of the 41 studied species showed significant difference between sexes; but, no consistent dominant pattern of THg concentrations favored either sex. When grouped by trophic levels, overall mean difference in THg concentrations between males and females were not statistically significant. Correlation analysis showed sex-dependent THg bio-accumulation as a function of condition factor was statistically significant and negative for all trophic levels (detritivorous, herviborous, omnivorous, planctivorous, carnivorous, and piscivorous).

CONCLUSIONS

Sex is not the main driver of Hg bioaccumulation in most Amazonian fish species; however, studies have to consider the minimum number of samples required to ascertain sex effects on THg bioaccumulation. Therefore, neither the surveillance of environmental pollution nor the current food advisories based on muscle THg need to change because of fish sex.

A pharyngeal jaw evolutionary innovation facilitated extinction in Lake Victoria cichlids

Evolutionary innovations, traits that give species access to previously unoccupied niches, may promote speciation and adaptive radiation. Here, we show that such innovations can also result in competitive inferiority and extinction. We present evidence that the modified pharyngeal jaws of cichlid fishes and several marine fish lineages, a classic example of evolutionary innovation, are not universally beneficial. A large-scale analysis of dietary evolution across marine fish lineages reveals that the innovation compromises access to energy-rich predator niches. We show that this competitive inferiority shaped the adaptive radiation of cichlids in Lake Tanganyika and played a pivotal and previously unrecognized role in the mass extinction of cichlid fishes in Lake Victoria after Nile perch invasion.

Co-extinction in a host-parasite network: identifying key hosts for network stability

Parasites comprise a substantial portion of total biodiversity. Ultimately, this means that host extinction could result in many secondary extinctions of obligate parasites and potentially alter host-parasite network structure. Here, we examined a highly resolved fish-parasite network to determine key hosts responsible for maintaining parasite diversity and network structure (quantified here as nestedness and modularity). We evaluated four possible host extinction orders and compared the resulting co-extinction dynamics to random extinction simulations; including host removal based on estimated extinction risk, parasite species richness and host level contributions to nestedness and modularity. We found that all extinction orders, except the one based on realistic extinction risk, resulted in faster declines in parasite diversity and network structure relative to random biodiversity loss. Further, we determined species-level contributions to network structure were best predicted by parasite species richness and host family. Taken together, we demonstrate that a small proportion of hosts contribute substantially to network structure and that removal of these hosts results in rapid declines in parasite diversity and network structure. As network stability can potentially be inferred through measures of network structure, our findings may provide insight into species traits that confer stability.

Mercury in fish of the Madeira river (temporal and spatial assessment), Brazilian Amazon

The Madeira River is the largest tributary of the Amazon River Basin and one of the most impacted by artisanal gold-mining activities, deforestation for agricultural projects, and recent hydroelectric reservoirs. Total Hg (and methylmercury-MeHg) concentrations was determined in 3182 fish samples of 84 species from different trophic levels as a function of standard size. Species at the top of the trophic level (Piscivorous, Carnivorous) showed the highest mean total Hg concentrations (51-1242 µg/kg), Planctivorous and Omnivorous species showed intermediate total Hg concentrations (26-494 µg/kg), while Detritivorous and Herbivorous species showed the lowest range of mean total Hg concentrations (9-275 µg/kg). Significant correlations between fish size (standard length) and total Hg concentrations were seen for Planctivorous (r=0.474, p=0.0001), Piscivorous (r=0.459, p=0.0001), Detritivorous (r=0.227, p=0.0001), Carnivorous (r=0.212, p=0.0001), and Herbivorous (r=0.156, p=0.01), but not for the Omnivorous species (r=-0.064, p=0.0685). Moreover, fish trophic levels influenced the ratio of MeHg to total Hg (ranged from 70% to 92%). When adjusted for standard body length, significant increases in Hg concentrations in the last 10 years were species specific. Spatial differences, albeit significant for some species, were not consistent with time trends for environmental contamination from past alluvial gold mining activities. Fish-Hg bioaccumulation is species specific but fish feeding strategies are the predominant influence in the fish-Hg bioaccumulation pattern.

A preliminary inventory of the catfishes of the lower Rio Nhamundá, Brazil (Ostariophysi, Siluriformes)

The Rio Nhamundá is a poorly-known clearwater river draining the southern Guiana Shield of Brazil. In this study we report the findings of a preliminary ichthyological survey, focusing on catfishes (Siluriformes). We identify a total of 36 species (31 genera, seven families) from the Nhamundá, including 11 species already recorded from the river. Overall, our survey results show that even rapid surveys can provide important information on Amazon fish biodiversity, suggesting potential new species, providing range extensions for nominal species, and additionally highlighting taxa in need of taxonomic revision and genetic study. As well as the traditional forms of data collected on biodiversity surveys (i.e. preserved specimen vouchers), our study also provides "new" types of data in the form of DNA barcodes and images of fishes exhibiting colouration in life, information that will be invaluable in future work addressing difficult groups. O Rio Nhamundá é um rio de água clara, pouco conhecido, que drena parte do Escudo das Guianas em território brasileiro. Nesse estudo, nós reportamos os resultados de um levantamento ictiofaunístico preliminar dessa área, tendo como foco os bagres (Siluriformes). Nós identificamos um total de 36 espécies (31 gêneros, sete famílias) provenientes de nossa coleta, e adicionamos 11 espécies já conhecidas para o rio. De maneira geral, os resultados de nossa pesquisa mostram que mesmo levantamentos rápidos podem gerar informações importantes sobre a biodiversidade de peixes amazônicos, sugerindo potenciais espécies novas, ampliando a área de distribuição de espécies, além de apontar a necessidade de revisões taxonômicas e estudos genéticos para alguns taxa. Para além das formas tradicionais de dados coletados em pesquisas de biodiversidade (i.e. espécimes preservados), nosso estudo fornece "novas" formas de dados, como DNA barcodes e imagens com o padrão de coloração dos espécimes vivos, informações essas que serão de valor inestimável para futuros estudos que abordem grupos taxonômicos difíceis.

Two waves of colonization straddling the K-Pg boundary formed the modern reef fish fauna

Living reef fishes are one of the most diverse vertebrate assemblages on Earth. Despite its prominence and ecological importance, the origins and assembly of the reef fish fauna is poorly described. A patchy fossil record suggests that the major colonization of reef habitats must have occurred in the Late Cretaceous and early Palaeogene, with the earliest known modern fossil coral reef fish assemblage dated to 50 Ma. Using a phylogenetic approach, we analysed the early evolutionary dynamics of modern reef fishes. We find that reef lineages successively colonized reef habitats throughout the Late Cretaceous and early Palaeogene. Two waves of invasion were accompanied by increasing morphological convergence: one in the Late Cretaceous from 90 to 72 Ma and the other immediately following the end-Cretaceous mass extinction. The surge in reef invasions after the Cretaceous–Palaeogene boundary continued for 10 Myr, after which the pace of transitions to reef habitats slowed. Combined, these patterns match a classic niche-filling scenario: early transitions to reefs were made rapidly by morphologically distinct lineages and were followed by a decrease in the rate of invasions and eventual saturation of morphospace. Major alterations in reef composition, distribution and abundance, along with shifts in climate and oceanic currents, occurred during the Late Cretaceous and early Palaeogene interval. A causal mechanism between these changes and concurrent episodes of reef invasion remains obscure, but what is clear is that the broad framework of the modern reef fish fauna was in place within 10 Myr of the end-Cretaceous extinction.

Evolutionary origin of the Scombridae (tunas and mackerels): members of a paleogene adaptive radiation with 14 other pelagic fish families

Uncertainties surrounding the evolutionary origin of the epipelagic fish family Scombridae (tunas and mackerels) are symptomatic of the difficulties in resolving suprafamilial relationships within Percomorpha, a hyperdiverse teleost radiation that contains approximately 17,000 species placed in 13 ill-defined orders and 269 families. Here we find that scombrids share a common ancestry with 14 families based on (i) bioinformatic analyses using partial mitochondrial and nuclear gene sequences from all percomorphs deposited in GenBank (10,733 sequences) and (ii) subsequent mitogenomic analysis based on 57 species from those targeted 15 families and 67 outgroup taxa. Morphological heterogeneity among these 15 families is so extraordinary that they have been placed in six different perciform suborders. However, members of the 15 families are either coastal or oceanic pelagic in their ecology with diverse modes of life, suggesting that they represent a previously undetected adaptive radiation in the pelagic realm. Time-calibrated phylogenies imply that scombrids originated from a deep-ocean ancestor and began to radiate after the end-Cretaceous when large predatory epipelagic fishes were selective victims of the Cretaceous-Paleogene mass extinction. We name this clade of open-ocean fishes containing Scombridae “Pelagia” in reference to the common habitat preference that links the 15 families.