Get unbent! R tools for the removal of arching and bending in fish specimens for geometric morphometric shape analysis

Geometric morphometrics is a powerful tool for studying fish body shape; however, body posture can be a hindrance to these analyses. Here I introduce new R language tools for correcting multiple types of bending of 3D data based on the TPS suite (geometric morphometric software) "unbend specimens" methodology. In a sample dataset of darters, these R tools adequately accounted for posture artifacts otherwise evident across multiple principal component axes. I hope these new tools will facilitate the incorporation of 3D landmark data into the comparative analysis of fish body shape.

Is shape in the eye of the beholder? Assessing landmarking error in geometric morphometric analyses on live fish

Geometric morphometrics is widely used to quantify morphological variation between biological specimens, but the fundamental influence of operator bias on data reproducibility is rarely considered, particularly in studies using photographs of live animals taken under field conditions. We examined this using four independent operators that applied an identical landmarking scheme to replicate photographs of 291 live Atlantic salmon (Salmo salar L.) from two rivers. Using repeated measures tests, we found significant inter-operator differences in mean body shape, suggesting that the operators introduced a systematic error despite following the same landmarking scheme. No significant differences were detected when the landmarking process was repeated by the same operator on a random subset of photographs. Importantly, in spite of significant operator bias, small but statistically significant morphological differences between fish from the two rivers were found consistently by all operators. Pairwise tests of angles of vectors of shape change showed that these between-river differences in body shape were analogous across operator datasets, suggesting a general reproducibility of findings obtained by geometric morphometric studies. In contrast, merging landmark data when fish from each river are digitised by different operators had a significant impact on downstream analyses, highlighting an intrinsic risk of bias. Overall, we show that, even when significant inter-operator error is introduced during digitisation, following an identical landmarking scheme can identify morphological differences between populations. This study indicates that operators digitising at least a sub-set of all data groups of interest may be an effective way of mitigating inter-operator error and potentially enabling data sharing.

Evidence for three morphotypes among anadromous Arctic char (Salvelinus alpinus) sampled in the marine environment

Variable resource use and responses to environmental conditions can lead to phenotypic diversity and distinct morphotypes within salmonids, including Arctic char (Salvelinus alpinus). Despite the cultural and economic importance of Arctic char in the Inuvialuit Settlement Region (ISR), limited data exist on the extent and presence of morphological diversity in this region. This is of concern for management given climate change impacts on regional fish populations. The authors investigated morphological diversity in anadromous Arctic char sampled during their summer marine migration-residency period when seasonal harvesting occurs in a coastal mixed-stock fishery. Geometric morphometric analysis was conducted using digital photographs of live Arctic char (n = 103) of which a sub-set was subsequently implanted with acoustic transmitters (n = 90) and released, and their overwintering lakes determined using active acoustic telemetry surveys. Twenty-three morphological landmarks were established and overlaid on digital images, and nine linear measurements of the body and head were recorded. Principle component analysis and K-means clustering based on linear measurements categorised fish into three morphotypes: slender body and slim head (n = 31), small and short head with a small mouth (n = 46) and elongated head shape with large mouth (n = 26). Tagged individuals of the three morphotypes occupied all lakes with no distinction observed. The three Arctic char morphotypes detected in this coastal mixed-stock fishery could represent adaptation to specific feeding-movement behaviours potentially tied to juvenile residency in freshwater systems, efficient exploitation of the marine prey pulse, or are relicts from ancestral types. To the authors' knowledge, this study is the first to identify distinct Arctic char morphotypes occurring in sympatry in the marine environment. Identifying phenotypic diversity will assist management to promote the sustainability of this regional fishery.

Parallel evolution of morphological traits and body shape in littoral and pelagic brook charr, Salvelinus fontinalis, along a gradient of interspecific competition

The parallel evolution of similar ecotypes in response to comparable environmental conditions is believed to reveal the importance of divergent selection in phenotypic diversifying processes. Systems characterized by the presence of multiple replicate populations expressing resource polymorphism thus provide an ideal opportunity to address the occurrence and factors affecting the parallel evolution of ecotypes. Previous studies have shown that brook charr (Salvelinus fontinalis) exhibit resource polymorphism in some Canadian Shield lakes, where a littoral ecotype feeds mainly on zoobenthos and a pelagic ecotype feeds mostly on zooplankton. Using morphological traits and geometric morphometric analyses on 18 native brook charr populations, we explicitly tested (i) whether brook charr ecotypes show parallel evolution across populations (i.e. the same morphological traits discriminate ecotypes among lakes) and (ii) whether interspecific competition decreases the amplitude of morphological differentiation between ecotypes, if any, because brook charr experience some level of competitive exclusion from the littoral habitat in the presence of creek chub or white sucker. We observed a low level of parallel evolution, where the littoral ecotype was overall stouter with longer fins and smaller eyes than the pelagic ecotype. Interspecific competition had no clear impacts on the amplitude of morphological differentiation. We also observed that inter-lake morphological differences are greater than between ecotypes within lakes, suggesting an important effect of local environmental factors on population morphology. Early-stage of diversification as well as phenotypic plasticity and morphological integration could explain why resource polymorphism is still subtle in brook charr populations.

The interaction of resource use and gene flow on the phenotypic divergence of benthic and pelagic morphs of Icelandic Arctic charr (Salvelinus alpinus)

Conceptual models of adaptive divergence and ecological speciation in sympatry predict differential resource use, phenotype-environment correlations, and reduced gene flow among diverging phenotypes. While these predictions have been assessed in past studies, connections among them have rarely been assessed collectively. We examined relationships among phenotypic, ecological, and genetic variation in Arctic charr (Salvelinus alpinus) from six Icelandic localities that have undergone varying degrees of divergence into sympatric benthic and pelagic morphs. We characterized morphological variation with geometric morphometrics, tested for differential resource use between morphs using stable isotopes, and inferred the amount of gene flow from single nucleotide polymorphisms. Analysis of stable isotopic signatures indicated that sympatric morphs showed similar difference in resource use across populations, likely arising from the common utilization of niche space within each population. Carbon isotopic signature was also a significant predictor of individual variation in body shape and size, suggesting that variation in benthic and pelagic resource use is associated with phenotypic variation. The estimated percentage of hybrids between sympatric morphs varied across populations (from 0% to 15.6%) but the majority of fish had genotypes (ancestry coefficients) characteristic of pure morphs. Despite evidence of reduced gene flow between sympatric morphs, we did not detect the expected negative relationship between divergence in resource use and gene flow. Three lakes showed the expected pattern, but morphs in the fourth showed no detectable hybridization and had relatively low differences in resource use between them. This coupled with the finding that resource use and genetic differentiation had differential effects on body shape variation across populations suggests that reproductive isolation maintains phenotypic divergence between benthic and pelagic morphs when the effects of resource use are relatively low. Our ability to assess relationships between phenotype, ecology, and genetics deepens our understanding of the processes underlying adaptive divergence in sympatry.

Investigating the morphological and genetic divergence of arctic char (Salvelinus alpinus) populations in lakes of arctic Alaska

Polymorphism facilitates coexistence of divergent morphs (e.g., phenotypes) of the same species by minimizing intraspecific competition, especially when resources are limiting. Arctic char (Salvelinus sp.) are a Holarctic fish often forming morphologically, and sometimes genetically, divergent morphs. In this study, we assessed the morphological and genetic diversity and divergence of 263 individuals from seven populations of arctic char with varying length-frequency distributions across two distinct groups of lakes in northern Alaska. Despite close geographic proximity, each lake group occurs on landscapes with different glacial ages and surface water connectivity, and thus was likely colonized by fishes at different times. Across lakes, a continuum of physical (e.g., lake area, maximum depth) and biological characteristics (e.g., primary productivity, fish density) exists, likely contributing to characteristics of present-day char populations. Although some lakes exhibit bimodal size distributions, using model-based clustering of morphometric traits corrected for allometry, we did not detect morphological differences within and across char populations. Genomic analyses using 15,934 SNPs obtained from genotyping by sequencing demonstrated differences among lake groups related to historical biogeography, but within lake groups and within individual lakes, genetic differentiation was not related to total body length. We used PERMANOVA to identify environmental and biological factors related to observed char size structure. Significant predictors included water transparency (i.e., a primary productivity proxy), char density (fish·ha-1), and lake group. Larger char occurred in lakes with greater primary production and lower char densities, suggesting less intraspecific competition and resource limitation. Thus, char populations in more productive and connected lakes may prove more stable to environmental changes, relative to food-limited and closed lakes, if lake productivity increases concomitantly. Our findings provide some of the first descriptions of genomic characteristics of char populations in arctic Alaska, and offer important consideration for the persistence of these populations for subsistence and conservation.

Ecological speciation in European whitefish is driven by a large-gaped predator

Lake-dwelling fish that form species pairs/flocks characterized by body size divergence are important model systems for speciation research. Although several sources of divergent selection have been identified in these systems, their importance for driving the speciation process remains elusive. A major problem is that in retrospect, we cannot distinguish selection pressures that initiated divergence from those acting later in the process. To address this issue, we studied the initial stages of speciation in European whitefish (Coregonus lavaretus) using data from 358 populations of varying age (26-10,000 years). We find that whitefish speciation is driven by a large-growing predator, the northern pike (Esox lucius). Pike initiates divergence by causing a largely plastic differentiation into benthic giants and pelagic dwarfs: ecotypes that will subsequently develop partial reproductive isolation and heritable differences in gill raker number. Using an eco-evolutionary model, we demonstrate how pike's habitat specificity and large gape size are critical for imposing a between-habitat trade-off, causing prey to mature in a safer place or at a safer size. Thereby, we propose a novel mechanism for how predators may cause dwarf/giant speciation in lake-dwelling fish species.

Heritable morphological differentiation in salmonids from two distinct stream types

We tested for phenotype-to-habitat associations in brown trout Salmo trutta populations from two ecologically different habitat types; i.e., groundwater and surface-water-fed streams. Additionally, we raised captive offspring from two such populations under standardised conditions to test whether potential phenotypic differentiation would be passed on to offspring. We found analogous differentiation by habitat in multiple wild populations. Some of these morphological differences were at least partially inherited by offspring. We suggest that this could have implications for both scientists and fisheries authorities studying or managing trout populations.

Changing environmental gradients over forty years alter ecomorphological variation in Guadalupe Bass Micropterus treculii throughout a river basin

Understanding the degree of intraspecific variation within and among populations is a key aspect of predicting the capacity of a species to respond to anthropogenic disturbances. However, intraspecific variation is usually assessed at either limited temporal, but broad spatial scales or vice versa, which can make assessing changes in response to long-term disturbances challenging. We evaluated the relationship between the longitudinal gradient of changing flow regimes and land use/land cover patterns since 1980 and morphological variation of Guadalupe Bass Micropterus treculii throughout the Colorado River Basin of central Texas. The Colorado River Basin in Texas has experienced major alterations to the hydrologic regime due to changing land- and water-use patterns. Historical collections of Guadalupe Bass prior to rapid human-induced change present the unique opportunity to study the response of populations to varying environmental conditions through space and time. Morphological differentiation of Guadalupe Bass associated with temporal changes in flow regimes and land use/land cover patterns suggests that they are exhibiting intraspecific trait variability, with contemporary individuals showing increased body depth, in response to environmental alteration through time (specifically related to an increase in herbaceous land cover, maximum flows, and the number of low pulses and high pulses). Additionally, individuals from tributaries with increased hydrologic alteration associated with urbanization or agricultural withdrawals tended to have a greater distance between the anal and caudal fin. These results reveal trait variation that may help to buffer populations under conditions of increased urbanization and sprawl, human population growth, and climate risk, all of which impose novel selective pressures, especially on endemic species like Guadalupe Bass. Our results contribute an understanding of the adaptability and capacity of an endemic population to respond to expected future changes based on demographic or climatic projection.

Morphological divergence between three Arctic charr morphs - the significance of the deep-water environment

Morphological divergence was evident among three sympatric morphs of Arctic charr (Salvelinus alpinus (L.)) that are ecologically diverged along the shallow-, deep-water resource axis in a subarctic postglacial lake (Norway). The two deep-water (profundal) spawning morphs, a benthivore (PB-morph) and a piscivore (PP-morph), have evolved under identical abiotic conditions with constant low light and temperature levels in their deep-water habitat, and were morphologically most similar. However, they differed in important head traits (e.g., eye and mouth size) related to their different diet specializations. The small-sized PB-morph had a paedomorphic appearance with a blunt head shape, large eyes, and a deep body shape adapted to their profundal lifestyle feeding on submerged benthos from soft, deep-water sediments. The PP-morph had a robust head, large mouth with numerous teeth, and an elongated body shape strongly related to their piscivorous behavior. The littoral spawning omnivore morph (LO-morph) predominantly utilizes the shallow benthic–pelagic habitat and food resources. Compared to the deep-water morphs, the LO-morph had smaller head relative to body size. The LO-morph exhibited traits typical for both shallow-water benthic feeding (e.g., large body depths and small eyes) and planktivorous feeding in the pelagic habitat (e.g., streamlined body shape and small mouth). The development of morphological differences within the same deep-water habitat for the PB- and PP-morphs highlights the potential of biotic factors and ecological interactions to promote further divergence in the evolution of polymorphism in a tentative incipient speciation process. The diversity of deep-water charr in this study represents a novelty in the Arctic charr polymorphism as a truly deep-water piscivore morph has to our knowledge not been described elsewhere.

Anadromous char as an alternate food choice to marine animals: a synthesis of Hg concentrations, population features and other influencing factors

This study was conducted to confirm sporadic measurements made over the late 1970s to the early 1990 s which determined that mercury (Hg) concentrations were low in anadromous char across Arctic and subarctic Canada including northern Québec and Labrador. Over 2004-2013, anadromous char populations across northern Canada were investigated at 20 sites for Hg concentrations and life history characteristics. Hg concentrations were extremely low in anadromous char muscle, typically <0.05 μg/g (wet weight) and, at each location, generally increased with fish length, age and nitrogen isotope (δ(15)N) ratio and decreased with condition factor and %lipid; correlations with carbon isotope (δ(13)C) ratio were inconsistent. Location and year were significant variables influencing Hg concentrations over the study area; longitude and latitude also were significant influencing variables. Char length, weight, age, condition factor and lipid content explained additional variance. A tendency towards higher Hg concentrations with increasing latitude may be partially related to decreasing growth of char towards the north. However, Hg concentrations in char were positively correlated with growth rates suggesting that Hg concentrations in char also were higher in the more productive study areas, including to the west where mainland riverine inputs of terrestrial carbon, nutrients, and Hg were greater. The data base for assessing time trends in char was limited by the small number of years investigated at most locations, variable fish size across years, small sample size, etc. Where temporal trends were detected, they were of increase on the long term (1970s, 1980s or early 1990 s to the present) but of decrease on the short term (early 2000s to present) with Nain (Labrador) showing the converse pattern. Higher Hg concentrations were also related to lower condition factor and cooler springs. Hg concentrations in anadromous char are compared with other terrestrial, aquatic and marine vertebrates in traditional diets. The known information on anadromous char is reviewed including population features, habitat, and harvests. Future Hg trend monitoring should focus on specific locations and harvest areas within these areas to better assess trends and influencing factors.

Morphological change and phenotypic plasticity in native and non-native pumpkinseed sunfish in response to competition

Non-indigenous species are oftentimes exposed to ecosystems with unfamiliar species, and organisms that exhibit a high degree of phenotypic plasticity may be better able to contend with the novel competitors that they may encounter during range expansion. In this study, differences in morphological plasticity were investigated using young-of-year pumpkinseed sunfish (Lepomis gibbosus) from native North American and non-native European populations. Two Canadian populations, isolated from bluegill sunfish (L. macrochirus) since the last glaciation, and two Spanish populations, isolated from bluegill since their introduction in Europe, were reared in a common environment using artificial enclosures. Fish were subjected to allopatric (without bluegill) or sympatric (with bluegill) conditions, and differences in plasticity were tested through a MANOVA of discriminant function scores. All pumpkinseed populations exhibited dietary shifts towards more benthivorous prey when held with bluegill. Differences between North American and European populations were observed in body dimensions, gill raker length and pelvic fin position. Sympatric treatments induced an increase in body width and a decrease in caudal peduncle length in native fish; non-native fish exhibited longer caudal peduncle lengths when held in sympatry with bluegill. Overall, phenotypic plasticity influenced morphological divergence less than genetic factors, regardless of population. Contrary to predictions, pumpkinseeds from Europe exhibited lower levels of phenotypic plasticity than Canadian populations, suggesting that European pumpkinseeds are more canalized than their North American counterparts.

Resource polymorphism and diversity of Arctic charr Salvelinus alpinus in a series of isolated lakes

Morphological, dietary and life-history variation in Arctic charr Salvelinus alpinus were characterized from three geographically proximate, but isolated lakes and one large lake into which they drain in south-western Alaska. Polymorphism was predicted to occur in the first three lakes because S. alpinus tend to become polymorphic in deep, isolated lakes with few co-occurring species. Only one morph was evident in the large lake and two of the three isolated lakes. In the third isolated lake, Lower Tazimina Lake, small and large morphs were found, the latter including two forms differing in growth rate. The small morph additionally differed from the two large forms by having more gill rakers and a deeper body than same-sized individuals of the large morph, consuming more limnetic and fewer benthic resources, having a greater gonado-somatic index and maturing at a smaller size. The two large forms consumed only slightly different foods (more terrestrial insects were consumed by the medium-growth form; more snails by the high-growth form). Trends in consumption of resources with body shape also differed between lakes. Variability in life history of S. alpinus in these Alaskan lakes was as broad as that found elsewhere. This variability is important for understanding lake ecosystems of remote regions where this species is commonly dominant.