Comparison of Migratory and Resident Populations of Brown Trout Reveals Candidate Genes for Migration Tendency

Thermotactic behaviour in lacustrine and riverine forms of Salmo trutta and its relevance to an emerging parasitic disease (PKD) in the wake of climate change

The thermotactic response of brown trout (Salmo trutta) was examined with the goal to investigate potential effects of the emerging temperature-dependent fatal trout disease PKD (proliferative kidney disease). First the differences in cold-water preferences of two forms of brown trout, lacustrine (migratory) and riverine, were determined. Second, it was studied whether this preference was changed in fish infected with PKD. The experiment involved a one-week habituation period at 14 °C in a two-chamber runway followed by a week of 3 °C temperature difference between the two runways. The fish could freely move between lanes via an opening at the end where food was provided. The temperature manipulation was repeated twice, and there were 3 trials per experimental group. All fish developed a clear spatial preference in the test. Lacustrine trout demonstrated a preference for warmer water, while riverine trout preferred cooler water. This may increase the risk to PKD in the lacustrine form. Most strikingly, riverine trout experimentally exposed to Tetracapsuloides bryosalmonae, the parasite that causes PKD, demonstrated stronger cold-seeking behaviour than control fish. Cold seeking behaviour suggests the occurrence of a disease-induced behavioural chill response, which may play an important role in disease recovery. This demonstrates the significance of protecting river connectivity and cold-water sanctuaries as management strategies for preserving salmonid populations in a warming climate.

Genomics reveals repeated landlocking of diadromous fish on an isolated island

Landlocking of diadromous fish in freshwater systems can have significant genomic consequences. For instance, the loss of the migratory life stage can dramatically reduce gene flow across populations, leading to increased genetic structuring and stronger effects of local adaptation. These genomic consequences have been well-studied in some mainland systems, but the evolutionary impacts of landlocking in island ecosystems are largely unknown. In this study, we used a genotyping-by-sequencing (GBS) approach to examine the evolutionary history of landlocking in common smelt (Retropinna retropinna) on Chatham Island, a small isolated oceanic island 800 kilometres east of mainland New Zealand. We examined the relationship between Chatham Island and mainland smelt and used coalescent analyses to test the number and timing of landlocking events on Chatham Island. Our genomic analysis, based on 21,135 SNPs across 169 individuals, revealed that the Chatham Island smelt was genomically distinct from the mainland New Zealand fish, consistent with a single ancestral colonisation event of Chatham Island in the Pleistocene. Significant genetic structure was also evident within the Chatham Island smelt, with a diadromous Chatham Island smelt group, along with three geographically structured landlocked groups. Coalescent demographic analysis supported three independent landlocking events, with this loss of diadromy significantly pre-dating human colonisation. Our results illustrate how landlocking of diadromous fish can occur repeatedly across a narrow spatial scale, and highlight a unique system to study the genomic basis of repeated adaptation.

Genetic and environmental drivers of migratory behavior in western burrowing owls and implications for conservation and management

Migration is driven by a combination of environmental and genetic factors, but many questions remain about those drivers. Potential interactions between genetic and environmental variants associated with different migratory phenotypes are rarely the focus of study. We pair low coverage whole genome resequencing with a de novo genome assembly to examine population structure, inbreeding, and the environmental factors associated with genetic differentiation between migratory and resident breeding phenotypes in a species of conservation concern, the western burrowing owl (Athene cunicularia hypugaea). Our analyses reveal a dichotomy in gene flow depending on whether the population is resident or migratory, with the former being genetically structured and the latter exhibiting no signs of structure. Among resident populations, we observed significantly higher genetic differentiation, significant isolation-by-distance, and significantly elevated inbreeding. Among migratory breeding groups, on the other hand, we observed lower genetic differentiation, no isolation-by-distance, and substantially lower inbreeding. Using genotype-environment association analysis, we find significant evidence for relationships between migratory phenotypes (i.e., migrant versus resident) and environmental variation associated with cold temperatures during the winter and barren, open habitats. In the regions of the genome most differentiated between migrants and residents, we find significant enrichment for genes associated with the metabolism of fats. This may be linked to the increased pressure on migrants to process and store fats more efficiently in preparation for and during migration. Our results provide a significant contribution toward understanding the evolution of migratory behavior and vital insight into ongoing conservation and management efforts for the western burrowing owl.

A comparative analysis of multi-biomarker responses to environmental stress: Evaluating differences in landfill leachate and pathogenic oomycete effects between wild and captive Salmo trutta

Phenotypic plasticity is one of the major means by which organisms can manage with environmental factor changes. Captivity-related stress and artificial rearing settings have been shown to dramatically alter fish response plasticity in terms of physiology, behavior, and health, potentially reducing overall fitness and fish survival. Understanding the variations in plasticity between captive-bred (kept in a homogenous environment) and wild fish populations in response to varied environmental pressures is becoming increasingly important, particularly in risk assessment research. In this study, we investigated whether captive-bred trout (Salmo trutta) are more susceptible to stress stimuli than their wild counterparts. In both wild and captive-bred trout, we investigated a battery of biomarkers that depicts the effects at various levels of biological organization in response to landfill leachate as a chemical pollutant, and after exposure to pathogenic oomycetes Saprolegnia parasitica. According to the findings, wild trout were more susceptible to chemical stimuli based on cytogenetic damage and catalase activity changes, whereas captive-bred trout were more sensitive to biological stress as evidenced by changes in overall fish activity and increasing cytogenetic damage in gills erythrocytes. Our findings emphasize the significance of exercising caution when conducting risk assessments of environmental pollutants using captive-bred animals, especially when seeking to extrapolate hazards and better understand the consequences of environmental contamination on wild fish populations. Additional comparative studies are required to investigate the impact of environmental stressors on multi-biomarker responses in both wild and captive fish populations in order to uncover changes in the plasticity of various traits that can result in adaptation or maladaptation to environmental stimuli within these fish populations, affecting data comparability and transferability to wildlife.

Spawning migration behaviour of sea trout (Salmo trutta L.) in a boreal river system: effects of flow conditions and obstacles on migratory activity

In this study, radio telemetry was used to examine the upstream spawning migration behaviour of anadromous brown trout (sea trout), Salmo trutta L., in a boreal river system, the River Isojoki, western Finland. The aim was to study the movement activity and migration characteristics of trout during the upstream spawning migration, as well as to locate the important spawning habitats and study the spawning characteristics. Furthermore, the authors analysed how flow conditions and a hydropower dam, with adjacent fishways, affected the upstream spawning migration. Tagged trout spawned in both the main stem and four tributaries, with spawning taking place from early October to November. The movement activity of radio-tagged trout was influenced by a hydropower dam (Perus dam), with spring migrators spending prolonged periods at the dam area, postponing the migration upstream. Flow conditions affected the total time spent at the dam area, as well as the movement activity in the free-flowing sections above the dam, with increasing flow stimulating activity. In addition, time of river ascent and location of spawning area had a significant effect on the movement activity of tagged trout. These results are further evidence that synergistic effects of flow and migratory obstacles can negatively influence migrations of anadromous fish, regardless of constructed fishways. The management of flow regimes and the efficiency of fishways are vital, as climate change will likely influence the flow and increase the water temperature of boreal river systems, further aggravating issues caused by obstacles.

Whole-genome resequencing reveals recent signatures of selection in five populations of largemouth bass ( Micropterus salmoides)

Largemouth bass ( Micropterus salmoides) is an economically important fish species in North America, Europe, and China. Various genetic improvement programs and domestication processes have modified its genome sequence through selective pressure, leaving nucleotide signals that can be detected at the genomic level. In this study, we sequenced 149 largemouth bass fish, including protospecies (imported from the US) and improved breeds (four domestic breeding populations from China). We detected genomic regions harboring certain genes associated with improved traits, which may be useful molecular markers for practical domestication, breeding, and selection. Subsequent analyses of genetic diversity and population structure revealed that the improved breeds have undergone more rigorous genetic changes. Through selective signal analysis, we identified hundreds of putative selective sweep regions in each largemouth bass line. Interestingly, we predicted 103 putative candidate genes potentially subjected to selection, including several associated with growth (p sst1 and grb10), early development ( klf9, sp4, and sp8), and immune traits ( pkn2, sept2, bcl6, and ripk2). These candidate genes represent potential genomic landmarks that could be used to improve important traits of biological and commercial interest. In summary, this study provides a genome-wide map of genetic variations and selection footprints in largemouth bass, which may benefit genetic studies and accelerate genetic improvement of this economically important fish.

New metabolic signature for Chagas disease reveals sex steroid perturbation in humans and mice

The causative agent of Chagas disease (CD), Trypanosoma cruzi, claims thousands of lives each year. Current diagnostic tools are insufficient to ensure parasitological detection in chronically infected patients has been achieved. A host-derived metabolic signature able to distinguish CD patients from uninfected individuals and assess antiparasitic treatment efficiency is introduced. Serum samples were collected from chronic CD patients, prior to and three years after treatment, and subjected to untargeted metabolomics analysis against demographically matched CD-negative controls. Five metabolites were confirmed by high-resolution tandem mass spectrometry. Several database matches for sex steroids were significantly altered in CD patients. A murine experiment corroborated sex steroid perturbation in T. cruzi-infected mice, particularly in male animals. Proteomics analysis also found increased steroidogenesis in the testes of infected mice. Metabolic alterations identified in this study shed light on the pathogenesis and provide the basis for developing novel assays for the diagnosis and screening of CD patients.

ddRAD-seq reveals the genetic structure and detects signals of selection in Italian brown trout

Background

Brown trout is one of the most widespread fresh-water fish species in Europe. The evolutionary history of and phylogenetic relationships between brown trout populations are complex, and this is especially true for Italian populations, which are heavily influenced in different ways by stocking practices. The characterization of the genetic structure of Italian brown trout populations may give information on the risk of losing endemic Italian populations due to lack of genetic diversity or to admixture with stocking populations. The identification of signatures of selection, and the information deriving from dense genotyping data will help genotype-informed breeding programs. We used a ddRAD-seq approach to obtain more than 100,000 single nucleotide polymorphisms (SNPs), and to characterize the population structure and signatures of selection in 90 brown trout samples.

Results

Italian brown trout populations are genetically differentiated, although the stocking practices have introduced strong admixture in endemic Italian trout, especially with the Atlantic lineage. Most of the analysed populations showed high levels of kinship and inbreeding. We detected putative signatures of selection using different approaches, and investigated if the regions were enriched for functional categories. Several regions putatively under selection and characterized by a reduction in heterozygosity across all the studied populations are enriched for genes involved in the response to viral infections.

Conclusions

Our results, which show evidence of admixture with the Atlantic lineage (commonly used for stocking), confirm the need for controlling stocking practices, in order to avoid the erosion of the endemic gene pool; given the apparently high levels of kinship and inbreeding in local populations, our results also show the need to take action for increasing gene diversity. In addition, we used the genetically-distinct lineages to detect signatures of selection and we identified putative signatures of selection in several regions associated with resistance to infectious diseases. These constitute candidate regions for the study of resistance to infections in wild and farmed trout.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12711-022-00698-7.

Drivers of neutral and adaptive differentiation in pike (Esox lucius) populations from contrasting environments

Understanding how eco-evolutionary processes and environmental factors drive population differentiation and adaptation are key challenges in evolutionary biology of relevance for biodiversity protection. Differentiation requires at least partial reproductive separation, which may result from different modes of isolation such as geographic isolation (allopatry) or isolation by distance (IBD), resistance (IBR), and environment (IBE). Despite that multiple modes might jointly influence differentiation, studies that compare the relative contributions are scarce. Using RADseq, we analyse neutral and adaptive genetic diversity and structure in 11 pike (Esox lucius) populations from contrasting environments along a latitudinal gradient (54.9-63.6°N), to investigate the relative effects of IBD, IBE and IBR, and to assess whether the effects differ between neutral and adaptive variation, or across structural levels. Patterns of neutral and adaptive variation differed, probably reflecting that they have been differently affected by stochastic and deterministic processes. The importance of the different modes of isolation differed between neutral and adaptive diversity, yet were consistent across structural levels. Neutral variation was influenced by interactions among all three modes of isolation, with IBR (seascape features) playing a central role, wheares adaptive variation was mainly influenced by IBE (environmental conditions). Taken together, this and previous studies suggest that it is common that multiple modes of isolation interactively shape patterns of genetic variation, and that their relative contributions differ among systems. To enable identification of general patterns and understand how various factors influence the relative contributions, it is important that several modes are simultaneously investigated in additional populations, species and environmental settings.

Genetic Causes and Consequences of Sympatric Morph Divergence in Salmonidae: A Search for Mechanisms

Repeatedly and recently evolved sympatric morphs exhibiting consistent phenotypic differences provide natural experimental replicates of speciation. Because such morphs are observed frequently in Salmonidae, this clade provides a rare opportunity to uncover the genomic mechanisms underpinning speciation. Such insight is also critical for conserving salmonid diversity, the loss of which could have significant ecological and economic consequences. Our review suggests that genetic differentiation among sympatric morphs is largely nonparallel apart from a few key genes that may be critical for consistently driving morph differentiation. We discuss alternative levels of parallelism likely underlying consistent morph differentiation and identify several factors that may temper this incipient speciation between sympatric morphs, including glacial history and contemporary selective pressures. Our synthesis demonstrates that salmonids are useful for studying speciation and poses additional research questions to be answered by future study of this family. Expected final online publication date for the Annual Review of Animal Biosciences, Volume 10 is February 2022. Please see http://www.annualreviews.org/page/journal/pubdates for revised estimates.

Detection of selection signatures in the genome of a farmed population of anadromous rainbow trout (Oncorhynchus mykiss)

Domestication processes and artificial selection are likely to leave signatures that can be detected at a molecular level in farmed rainbow trout (Oncorhynchus mykiss). These signatures of selection are genomic regions that contain functional genetic variants conferring a higher fitness to their bearers. We genotyped 749 rainbow trout from a commercial population using a rainbow trout Axiom 57 K SNP array panel and identified putative genomic regions under selection using the pcadapt, Composite Likelihood Ratio (CLR) and Integrated Haplotype Score (iHS) methods. After applying quality-control pipelines and statistical analyses, we detected 12, 96 and 16 SNPs putatively under selection, associated with 96, 781 and 115 candidate genes, respectively. Several of these candidate genes were associated with growth, early development, reproduction, behavior and immune system traits. In addition, some of the SNPs were found in interesting regions located in autosomal inversions on Omy05 and Omy20. These findings could represent a genome-wide map of selection signatures in farmed rainbow trout and could be important in explaining domestication and selection for genetic traits of commercial interest.

Long-term monitoring of a brown trout (Salmo trutta) population reveals kin-associated migration patterns and contributions by resident trout to the anadromous run

BACKGROUND

In species showing partial migration, as is the case for many salmonid fishes, it is important to assess how anthropogenic pressure experienced by migrating individuals affects the total population. We focused on brown trout (Salmo trutta) from the Guddal River in the Norwegian Hardanger Fjord system, which encompasses both resident and anadromous individuals. Aquaculture has led to increased anthropogenic pressure on brown trout during the marine phase in this region. Fish traps in the Guddal River allow for sampling all ascending anadromous spawners and descending smolts. We analyzed microsatellite DNA markers from all individuals ascending in 2006-2016, along with all emigrating smolts in 2017. We investigated (1) if there was evidence for declines in census numbers and effective population size during that period, (2) if there was association between kinship and migration timing in smolts and anadromous adults, and (3) to what extent resident trout were parents of outmigrating smolts.

RESULTS

Census counts of anadromous spawners showed no evidence for a decline from 2006 to 2016, but were lower than in 2000-2005. Estimates of effective population size also showed no trends of declines during the study period. Sibship reconstruction of the 2017 smolt run showed significant association between kinship and migration timing, and a similar association was indicated in anadromous spawners. Parentage assignment of 2017 smolts with ascending anadromous trout as candidate parents, and assuming that unknown parents represented resident trout, showed that 70% of smolts had at least one resident parent and 24% had two resident parents.

CONCLUSIONS

The results bear evidence of a population that after an initial decline has stabilized at a lower number of anadromous spawners. The significant association between kinship and migration timing in smolts suggests that specific episodes of elevated mortality in the sea could disproportionally affect some families and reduce overall effective population size. Finally, the results based on parentage assignment demonstrate a strong buffering effect of resident trout in case of elevated marine mortality affecting anadromous trout, but also highlight that increased mortality of anadromous trout, most of which are females, may lower overall production in the system.

Differences in growth between offspring of anadromous and freshwater brown trout Salmo trutta

In this study, individual growth of juvenile offspring of anadromous and freshwater resident brown trout Salmo trutta and crosses between the two from the River Imsa, Norway, was estimated. The juveniles were incubated until hatching at two temperatures (±S.D.), either 4.4 ± 1.5°C or 7.1 ± 0.6°C. Growth rate was estimated for 22 days in August-September when the fish on average were c. 8 g in wet mass, and the estimates were standardized to 1 g fish dry mass. Offspring of anadromous S. trutta grew better at both 15 and 18°C than offspring of freshwater resident S. trutta or offspring of crosses between the two S. trutta types. This difference appears not to result from a maternal effect because anadromous S. trutta grew better than the hybrids with anadromous mothers. Instead, this appears to be an inherited difference between the anadromous and the freshwater resident fish lending support to the hypothesis that anadromous and freshwater resident S. trutta in this river differ in genetic expression. Egg incubation temperature of S. trutta appeared not to influence the later growth as reported earlier from the studies of Atlantic salmon Salmo salar.

Alternative migratory tactics in brown trout (Salmo trutta) are underpinned by divergent regulation of metabolic but not neurological genes

The occurrence of alternative morphs within populations is common, but the underlying molecular mechanisms remain poorly understood. Many animals, for example, exhibit facultative migration, where two or more alternative migratory tactics (AMTs) coexist within populations. In certain salmonid species, some individuals remain in natal rivers all their lives, while others (in particular, females) migrate to sea for a period of marine growth. Here, we performed transcriptional profiling ("RNA-seq") of the brain and liver of male and female brown trout to understand the genes and processes that differentiate between migratory and residency morphs (AMT-associated genes) and how they may differ in expression between the sexes. We found tissue-specific differences with a greater number of genes expressed differentially in the liver (n = 867 genes) compared with the brain (n = 10) between the morphs. Genes with increased expression in resident livers were enriched for Gene Ontology terms associated with metabolic processes, highlighting key molecular-genetic pathways underlying the energetic requirements associated with divergent migratory tactics. In contrast, smolt-biased genes were enriched for biological processes such as response to cytokines, suggestive of possible immune function differences between smolts and residents. Finally, we identified evidence of sex-biased gene expression for AMT-associated genes in the liver (n = 12) but not the brain. Collectively, our results provide insights into tissue-specific gene expression underlying the production of alternative life histories within and between the sexes, and point toward a key role for metabolic processes in the liver in mediating divergent physiological trajectories of migrants versus residents.

Inversion breakpoints and the evolution of supergenes

The coexistence of discrete morphs that differ in multiple traits is common within natural populations of many taxa. Such morphs are often associated with chromosomal inversions, presumably because the recombination suppressing effects of inversions help maintain alternate adaptive combinations of alleles across the multiple loci affecting these traits. However, inversions can also harbour selected mutations at their breakpoints, leading to their rise in frequency in addition to (or independent from) their role in recombination suppression. In this review, we first describe the different ways that breakpoints can create mutations. We then critically examine the evidence for the breakpoint-mutation and recombination suppression hypotheses for explaining the existence of discrete morphs associated with chromosomal inversions. We find that the evidence that inversions are favoured due to recombination suppression is often indirect. The evidence that breakpoints harbour mutations that are adaptive is also largely indirect, with the characterization of inversion breakpoints at the sequence level being incomplete in most systems. Direct tests of the role of suppressed recombination and breakpoint mutations in inversion evolution are thus needed. Finally, we emphasize how the two hypotheses of recombination suppression and breakpoint mutation can act in conjunction, with implications for understanding the emergence of supergenes and their evolutionary dynamics. We conclude by discussing how breakpoint characterization could improve our understanding of complex, discrete phenotypic forms in nature.

Know your enemy - transcriptome of myxozoan Tetracapsuloides bryosalmonae reveals potential drug targets against proliferative kidney disease in salmonids

The myxozoan Tetracapsuloides bryosalmonae is a widely spread endoparasite that causes proliferative kidney disease (PKD) in salmonid fish. We developed an in silico pipeline to separate transcripts of T. bryosalmonae from the kidney tissue of its natural vertebrate host, brown trout (Salmo trutta). After stringent filtering, we constructed a partial transcriptome assembly T. bryosalmonae, comprising 3427 transcripts. Based on homology-restricted searches of the assembled parasite transcriptome and Atlantic salmon (Salmo salar) proteome, we identified four protein targets (Endoglycoceramidase, Legumain-like protease, Carbonic anhydrase 2, Pancreatic lipase-related protein 2) for the development of anti-parasitic drugs against T. bryosalmonae. Earlier work of these proteins on parasitic protists and helminths suggests that the identified anti-parasitic drug targets represent promising chemotherapeutic candidates also against T. bryosalmonae, and strengthen the view that the known inhibitors can be effective in evolutionarily distant organisms. In addition, we identified differentially expressed T. bryosalmonae genes between moderately and severely infected fish, indicating an increased abundance of T. bryosalmonae sporogonic stages in fish with low parasite load. In conclusion, this study paves the way for future genomic research in T. bryosalmonae and represents an important step towards the development of effective drugs against PKD.

Investigating Diadromy in Fishes and Its Loss in an -Omics Era

Diadromy, the predictable movements of individuals between marine and freshwater environments, is biogeographically and phylogenetically widespread across fishes. Thus, despite the high energetic and potential fitness costs involved in moving between distinct environments, diadromy appears to be an effective life history strategy. Yet, the origin and molecular mechanisms that underpin this migratory behavior are not fully understood. In this review, we aim first to summarize what is known about diadromy in fishes; this includes the phylogenetic relationship among diadromous species, a description of the main hypotheses regarding its origin, and a discussion of the presence of non-migratory populations within diadromous species. Second, we discuss how recent research based on -omics approaches (chiefly genomics, transcriptomics, and epigenomics) is beginning to provide answers to questions on the genetic bases and origin(s) of diadromy. Finally, we suggest future directions for -omics research that can help tackle questions on the evolution of diadromy.

Genomic basis of the loss of diadromy in Galaxias maculatus: Insights from reciprocal transplant experiments

Diadromy is known for having major effects on the distribution and richness of aquatic species, and so does its loss. The loss of diadromy has led to the diversification of many species, yet research focusing on understanding its molecular basis and consequences are limited. This is particularly true for amphidromous species despite being the most abundant group of diadromous species. Galaxias maculatus, an amphidromous species and one of the most widely distributed fishes in the Southern Hemisphere, exhibits many instances of nonmigratory or resident populations. The existence of naturally replicated resident populations in Patagonia can serve as an ideal system for the study of the mechanisms that lead to the loss of the diadromy and its ecological and evolutionary consequences. Here, we studied two adjacent river systems in which resident populations are genetically differentiated yet derived from the same diadromous population. By combining a reciprocal transplant experiment with genomic data, we showed that the two resident populations followed different evolutionary pathways by exhibiting a differential response in their capacity to survive in salt water. While one resident population was able to survive salt water, the other was not. Genomic analyses provided insights into the genes that distinguished (a) migratory from nonmigratory populations; (b) populations that can vs those that cannot survive a saltwater environment; and (c) between these resident populations. This study demonstrates that the loss of diadromy can be achieved by different pathways and that environmental (selection) and random (genetic drift) forces shape this dynamic evolutionary process.

Indirect genetic control of migration in a salmonid fish

Migration is a complex trait that often has genetic underpinnings. However, it is unclear if migratory behaviour itself is inherited (direct genetic control), or if the decision to migrate is instead the outcome of a set of physiological traits (indirect genetic control). For steelhead/rainbow trout (Oncorhynchus mykiss), migration is strongly linked to a large genomic region across their range. Here, we demonstrate a shared allelic basis between early life growth rate and migratory behaviour. Next, we demonstrate that early life growth differs among resident/migratory genotypes in wild juveniles several months prior to migration, with resident genotypes achieving a larger size in their first few months of life than migratory genotypes. We suggest that the genetic basis of migration is likely indirect and mediated by physiological traits such as growth rate. Evolutionary benefits of this indirect genetic mechanism likely include flexibility among individuals and persistence of life-history diversity within and among populations.

Whole genome re-sequencing reveals recent signatures of selection in three strains of farmed Nile tilapia (Oreochromis niloticus)

Nile tilapia belongs to the second most cultivated group of fish in the world, mainly because of its favorable characteristics for production. Genetic improvement programs and domestication process of Nile tilapia may have modified the genome through selective pressure, leaving signals that can be detected at the molecular level. In this work, signatures of selection were identified using genome-wide SNP data, by two haplotype-based (iHS and Rsb) and one F_ST based method. Whole-genome re-sequencing of 326 individuals from three strains (A, B and C) of farmed tilapia maintained in Brazil and Costa Rica was carried out using Illumina HiSeq 2500 technology. After applying conventional SNP-calling and quality-control filters, ~ 1.3 M high-quality SNPs were inferred and used as input for the iHS, Rsb and F_ST based methods. We detected several candidate genes putatively subjected to selection in each strain. A considerable number of these genes are associated with growth (e.g. NCAPG, KLF3, TBC1D1, TTN), early development (e.g. FGFR3, PFKFB3), and immunity traits (e.g. NLRC3, PIGR, MAP1S). These candidate genes represent putative genomic landmarks that could be associated to traits of biological and commercial interest in farmed Nile tilapia.

Northern European Salmo trutta (L.) populations are genetically divergent across geographical regions and environmental gradients

The salmonid fish Brown trout is iconic as a model for the application of conservation genetics to understand and manage local interspecific variation. However, there is still scant information about relationships between local and large-scale population structure, and to what extent geographical and environmental variables are associated with barriers to gene flow. We used information from 3,782 mapped SNPs developed for the present study and conducted outlier tests and gene-environment association (GEA) analyses in order to examine drivers of population structure. Analyses comprised >2,600 fish from 72 riverine populations spanning a central part of the species' distribution in northern Europe. We report hitherto unidentified genetic breaks in population structure, indicating strong barriers to gene flow. GEA loci were widely spread across genomic regions and showed correlations with climatic, abiotic and geographical parameters. In some cases, individual loci showed consistent GEA across the geographical regions Britain, Europe and Scandinavia. In other cases, correlations were observed only within a sub-set of regions, suggesting that locus-specific variation was associated with local processes. A paired-population sampling design allowed us to evaluate sampling effects on detection of outlier loci and GEA. Two widely applied methods for outlier detection (pcadapt and bayescan) showed low overlap in loci identified as statistical outliers across sub-sets of data. Two GEA analytical approaches (LFMM and RDA) showed good correspondence concerning loci associated with specific variables, but LFMM identified five times more statistically significant associations than RDA. Our results emphasize the importance of carefully considering the statistical methods applied for the hypotheses being tested in outlier analysis. Sampling design may have lower impact on results if the objective is to identify GEA loci and their population distribution. Our study provides new insights into trout populations, and results have direct management implications in serving as a tool for identification of conservation units.

Association Mapping Based on a Common-Garden Migration Experiment Reveals Candidate Genes for Migration Tendency in Brown Trout

A better understanding of the environmental and genetic contribution to migratory behavior and the evolution of traits linked to migration is crucial for fish conservation and fisheries management. Up to date, a few genes with unequivocal influence on the adoption of alternative migration strategies have been identified in salmonids. Here, we used a common garden set-up to measure individual migration distances of generally highly polymorphic brown trout Salmo trutta from two populations. Fish from the assumedly resident population showed clearly shorter migration distances than the fish from the assumed migratory population at the ages of 2 and 3 years. By using two alternative analytical pipelines with 22186 and 18264 SNPs obtained through RAD-sequencing, we searched for associations between individual migration distance, and both called genotypes and genotype probabilities. None of the SNPs showed statistically significant individual effects on migration after correction for multiple testing. By choosing a less stringent threshold, defined as an overlap of the top 0.1% SNPs identified by the analytical pipelines, GAPIT and Angsd, we identified eight candidate genes that are potentially linked to individual migration distance. While our results demonstrate large individual and population level differences in migration distances, the detected genetic associations were weak suggesting that migration traits likely have multigenic control.

Anadromy, potamodromy and residency in brown trout Salmo trutta: the role of genes and the environment

Brown trout Salmo trutta is endemic to Europe, western Asia and north-western Africa; it is a prominent member of freshwater and coastal marine fish faunas. The species shows two resident (river-resident, lake-resident) and three main facultative migratory life histories (downstream-upstream within a river system, fluvial-adfluvial potamodromous; to and from a lake, lacustrine-adfluvial (inlet) or allacustrine (outlet) potamodromous; to and from the sea, anadromous). River-residency v. migration is a balance between enhanced feeding and thus growth advantages of migration to a particular habitat v. the costs of potentially greater mortality and energy expenditure. Fluvial-adfluvial migration usually has less feeding improvement, but less mortality risk, than lacustrine-adfluvial or allacustrine and anadromous, but the latter vary among catchments as to which is favoured. Indirect evidence suggests that around 50% of the variability in S. trutta migration v. residency, among individuals within a population, is due to genetic variance. This dichotomous decision can best be explained by the threshold-trait model of quantitative genetics. Thus, an individual's physiological condition (e.g., energy status) as regulated by environmental factors, genes and non-genetic parental effects, acts as the cue. The magnitude of this cue relative to a genetically predetermined individual threshold, governs whether it will migrate or sexually mature as a river-resident. This decision threshold occurs early in life and, if the choice is to migrate, a second threshold probably follows determining the age and timing of migration. Migration destination (mainstem river, lake, or sea) also appears to be genetically programmed. Decisions to migrate and ultimate destination result in a number of subsequent consequential changes such as parr-smolt transformation, sexual maturity and return migration. Strong associations with one or a few genes have been found for most aspects of the migratory syndrome and indirect evidence supports genetic involvement in all parts. Thus, migratory and resident life histories potentially evolve as a result of natural and anthropogenic environmental changes, which alter relative survival and reproduction. Knowledge of genetic determinants of the various components of migration in S. trutta lags substantially behind that of Oncorhynchus mykiss and other salmonines. Identification of genetic markers linked to migration components and especially to the migration-residency decision, is a prerequisite for facilitating detailed empirical studies. In order to predict effectively, through modelling, the effects of environmental changes, quantification of the relative fitness of different migratory traits and of their heritabilities, across a range of environmental conditions, is also urgently required in the face of the increasing pace of such changes.

Evolutionary patterns of diadromy in fishes: more than a transitional state between marine and freshwater

Background

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.

Results

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.

Conclusion

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.

Electronic supplementary material

The online version of this article (10.1186/s12862-019-1492-2) contains supplementary material, which is available to authorized users.

Length polymorphisms at two candidate genes explain variation of migratory behaviors in blackpoll warblers (Setophaga striata)

Migratory behaviors such as the timing and duration of migration are genetically inherited and can be under strong natural selection, yet we still know very little about the specific genes or molecular pathways that control these behaviors. Studies in candidate genes Clock and Adcyap1 have revealed that both of these loci can be significantly correlated with migratory behaviors in birds, though observed relationships appear to vary across species. We investigated geographic genetic structure of Clock and Adcyap1 in four populations of blackpoll warblers (Setophaga striata), a Neotropical-Nearctic migrant that exhibits geographic variation in migratory timing and duration across its boreal breeding distribution. Further, we used data on migratory timing and duration, obtained from light-level geolocator trackers to investigate candidate genotype-phenotype relationships at the individual level. While we found no geographic structure in either candidate gene, we did find evidence that candidate gene lengths are correlated with five of the six migratory traits. Maximum Clock allele length was significantly and negatively associated with spring arrival date. Minimum Adcyap1 allele length was significantly and negatively associated with spring departure date and positively associated with fall arrival date at the wintering grounds. Additionally, we found a significant interaction between Clock and Adcyap1 allele lengths on both spring and fall migratory duration. Adcyap1 heterozygotes also had significantly shorter migration duration in both spring and fall compared to homozygotes. Our results support the growing body of evidence that Clock and Adcyap1 allele lengths are correlated with migratory behaviors in birds.

Genomic signatures and correlates of widespread population declines in salmon

Global losses of biodiversity are occurring at an unprecedented rate, but causes are often unidentified. Genomic data provide an opportunity to isolate drivers of change and even predict future vulnerabilities. Atlantic salmon (Salmo salar) populations have declined range-wide, but factors responsible are poorly understood. Here, we reconstruct changes in effective population size (Ne) in recent decades for 172 range-wide populations using a linkage-based method. Across the North Atlantic, Ne has significantly declined in >60% of populations and declines are consistently temperature-associated. We identify significant polygenic associations with decline, involving genomic regions related to metabolic, developmental, and physiological processes. These regions exhibit changes in presumably adaptive diversity in declining populations consistent with contemporary shifts in body size and phenology. Genomic signatures of widespread population decline and associated risk scores allow direct and potentially predictive links between population fitness and genotype, highlighting the power of genomic resources to assess population vulnerability.

Experimental crossbreeding reveals strain-specific variation in mortality, growth and personality in the brown trout (Salmo trutta)

Although hybridization between populations with low genetic diversity may induce heterosis, it can also lead to reduced fitness of hybrid offspring through outbreeding depression and loss of local adaptations. Using a half-sib mating design, we studied on brown trout (Salmo trutta) how hybridization of migratory hatchery-strain females with males from various strains would affect early mortality, growth and personality in F₁ offspring. No differences in mortality or alevin body length were found between the crossing groups by the end of the yolk-sac stage. At later developmental stages, higher mortality and slower growth in one of the geographically distant hybrid groups indicated potential outbreeding depression. The personality component indicating boldness and exploration tendency showed fairly low genetic variation and no phenotypic differences among the crossing groups while the personality component related to freezing behavior indicated stronger freezing responses in the purebred and local cross strain when compared to the two other strains. However, the purebred hatchery strain possessed stronger additive genetic tendency for boldness and explorative behavior, and weaker genetic tendency for freezing behavior, when compared to the wild × hatchery hybrid group. Our results add to the cumulating evidence of risks related to the stocking of fish strains from non-native origins.

Comparing RADseq and microsatellites for estimating genetic diversity and relatedness - Implications for brown trout conservation

The conservation and management of endangered species requires information on their genetic diversity, relatedness and population structure. The main genetic markers applied for these questions are microsatellites and single nucleotide polymorphisms (SNPs), the latter of which remain the more resource demanding approach in most cases. Here, we compare the performance of two approaches, SNPs obtained by restriction-site-associated DNA sequencing (RADseq) and 16 DNA microsatellite loci, for estimating genetic diversity, relatedness and genetic differentiation of three, small, geographically close wild brown trout (Salmo trutta) populations and a regionally used hatchery strain. The genetic differentiation, quantified as F ST, was similar when measured using 16 microsatellites and 4,876 SNPs. Based on both marker types, each brown trout population represented a distinct gene pool with a low level of interbreeding. Analysis of SNPs identified half- and full-siblings with a higher probability than the analysis based on microsatellites, and SNPs outperformed microsatellites in estimating individual-level multilocus heterozygosity. Overall, the results indicated that moderately polymorphic microsatellites and SNPs from RADseq agreed on estimates of population genetic structure in moderately diverged, small populations, but RADseq outperformed microsatellites for applications that required individual-level genotype information, such as quantifying relatedness and individual-level heterozygosity. The results can be applied to other small populations with low or moderate levels of genetic diversity.