Search

Recent adaptive radiations provide experimental opportunities to parse the relationship between genomic variation and the origins of distinct phenotypes. Sympatric radiations of the charr complex (genus Salvelinus) present a trove for phylogenetic analyses as charrs have repeatedly diversified into multiple morphs with distinct feeding specializations. However, charr species flocks normally comprise only two to three lineages. Dolly Varden charr inhabiting Lake Kronotskoe represent the most extensive radiation described for the genus, containing at least seven lineages, each with defining morphological and ecological traits. Here, we perform the first genome-wide analysis of this species flock to parse the foundations of adaptive change. Our data support distinct, reproductively isolated lineages within the clade. We find that changes in genes associated with thyroid signaling and craniofacial development provided a foundational shift in evolution to the lake. The thyroid axis is further implicated in subsequent lineage partitioning events. These results delineate a genetic scenario for the diversification of specialized lineages and highlight a common axis of change biasing the generation of specific forms during adaptive radiation., Competing Interests: Competing interests The authors declare no competing or financial interests., (© 2024. Published by The Company of Biologists Ltd.)

Meristic characters are often used to differentiate between closely related forms, morphs, and species of fishes, and lend insight into ecology and post-glacial recolonization in taxa with complicated or contentious phylogenies, including the genus Salvelinus. Previous studies of meristics in Salvelinus have focused mostly on individual populations. We collated data from 456 populations/systems across the North American and Russian Arctic and sub-Arctic, and found that counts of pyloric caeca and gill rakers differed consistently between fish visually and/or genetically identified as Arctic char and Dolly Varden across their distributional ranges., (© 2024 His Majesty the King in Right of Canada and The Author(s). Journal of Fish Biology published by John Wiley & Sons Ltd on behalf of Fisheries Society of the British Isles. Reproduced with the permission of the Minister of Fisheries and Oceans Canada.)

Hybridization between divergent lineages can result in losses of distinct evolutionary taxa. Alternatively, hybridization can lead to increased genetic variability that may fuel local adaptation and the generation of novel traits and/or taxa. Here, we examined single-nucleotide polymorphisms generated using genotyping-by-sequencing in a population of Dolly Varden char (Pisces: Salmonidae) that is highly admixed within a contact zone between two subspecies (Salvelinus malma malma, Northern Dolly Varden [NDV] and S. m. lordi, Southern Dolly Varden [SDV]) in southwestern Alaska to assess the spatial distribution of hybrids and to test hypotheses on the origin of the admixed population. Ancestry analysis revealed that this admixed population is composed of advanced generation hybrids between NDV and SDV or advanced backcrosses to SDV; no F 1 hybrids were detected. Coalescent-based demographic modelling supported the origin of this population about 55,000 years ago by secondary contact between NDV and SDV with low levels of contemporary gene flow. Ancestry in NDV and SDV varies within the watershed and ancestry in NDV was positively associated with distance upstream from the sea, contingent on habitat-type sampled, and negatively associated with the number of migrations that individual fish made to the sea. Our results suggest that divergence between subspecies over hundreds of thousands of years may not be associated with significant reproductive isolation, but that elevated diversity owing to hybridization may have contributed to adaptive divergence in habitat use and life history., (© 2024 The Author(s). Molecular Ecology published by John Wiley & Sons Ltd.)

Course-based undergraduate research experiences (CUREs) increase student access to high impact research experiences. CUREs engage students in the scientific process by learning how to pose scientific questions, develop hypotheses, and generate data to test them. Environmental DNA (eDNA) is a growing field of research that is gaining accessibility through decreasing laboratory costs, which can make a foundation for multiple, engaging CUREs. This manuscript describes three case studies that used eDNA in an upper year undergraduate course. The first focusses on a systematic literature review of eDNA metadata reporting. The second describes the biomonitoring of brook trout in southern Ontario using eDNA. The third involves eDNA metabarcoding for freshwater fish detection in southern Ontario. Undergraduates were involved in the development and execution of experiments, scientific communication, the peer review process, and fundraising. Through this manuscript, we show the novel application of eDNA CUREs and provide a roadmap for other instructors interested in implementing similar projects. Interviews with seven students from these courses indicate the benefits experienced from taking these courses. We argue that the use of eDNA in CUREs should be expanded in undergraduate biology programs due to the benefit to students and the increasing accessibility of this technology., Competing Interests: The authors declare no competing interests.

High-temperature stress poses a significant environmental challenge for aquatic organisms, including tsinling lenok trout (Brachymystax lenok tsinlingensis). This study aimed to investigate the role of microRNAs (miRNAs) in inducing liver inflammation in tsinling lenok trout under high-temperature stress. Tsinling lenok trout were exposed to high-temperature conditions (24 °C) for 8 h, and liver samples were collected for analysis. Through small RNA sequencing, we identified differentially expressed miRNAs in the liver of high-temperature-stressed tsinling lenok trout compared to the control group (maintained at 16 °C). Several miRNAs, including novel-m0105-5p and miR-8159-x, showed significant changes in expression levels. Additionally, we conducted bioinformatics analysis to explore the potential target genes of these differentially expressed miRNAs. Our findings revealed that these miRNA target genes are involved in inflammatory response pathways, such as NFKB1 and MAP3K5. The downregulation of novel-m0105-5p and miR-8159-x in the liver of high-temperature-stressed tsinling lenok trout suggests their role in regulating liver inflammatory responses. To validate this, we performed a dual-luciferase reporter assay to confirm the regulatory relationship between miRNAs and target genes. Our results demonstrated that novel-m0105-5p and miR-8159-x enhance the inflammatory response of hepatocytes by promoting the expression of NFKB1 and MAP3K5, respectively. In conclusion, our study provides evidence that high-temperature stress induces liver inflammation in tsinling lenok trout through dysregulation of miRNAs. Understanding the molecular mechanisms underlying the inflammatory response in tsinling lenok trout under high-temperature stress is crucial for developing strategies to mitigate the negative impacts of environmental stressors on fish health and aquaculture production., (© 2024. The Author(s), under exclusive licence to Springer Science+Business Media, LLC, part of Springer Nature.)

The study investigated if gonad maturation in triploid brown trout, Salmo trutta, was entirely suppressed or only delayed, and if triploids could interbreed with diploid counterparts. Ten percent of the total number of 3-year-old triploid S. trutta, 15% of 4-year-old fish, and 17% of 5-year-old fish produced semen. Three and 4 years old triploid fish did not produce eggs, but 15% of the 5-year-old fish did so. The quantity and sperm motility of triploid semen did not differ from diploids, but the sperm concentration was significantly lower. When diploid eggs were fertilized with triploid semen (3n × 2n crosses), the percentage of eyed stage embryos, of hatched larvae, and of normal-shaped larvae did not differ from the diploid controls. Circa 90% of 3n × 2n crosses had a ploidy level of 2.4n. In the remaining percentage of 3n × 2n crosses, the ploidy level was ≥2n and <2.4n. In sperm competition experiments where diploid eggs were fertilized with a mixture of diploid and triploid semen, 52% of the originating larvae had a ploidy level of 2n, 43% of 2.4n, and 5% of the fish were not exactly classified. From the start of feeding to an age of 248 days, the mortality rate of 3n × 2n interploid crosses and of 2n × 2n controls was similar. The growth of interploid crosses was significantly higher than that of controls. In triploid mature females, the egg mass per kilogram of body weight was significantly lower than in diploids. The mass of the non-hardened eggs and the percentile weight increase during hardening did not differ from diploid eggs. When triploid eggs were fertilized with diploid semen (2n × 3n crosses), the development rate to normal hatched larvae was less than 10%. All originating larvae had a ploidy level of 3n. From the start of feeding to an age of 248 days, 2n × 3n crosses had a higher mortality rate (15%) than diploid controls (<5%). Growth of this type of interploid crosses was reduced in comparison to controls. Therefore, triploids introduced into natural waters for recreational fisheries or escaping from farms may interbreed with diploid counterparts. This not only alters the genotypes of local populations but also changes the ploidy levels., (© 2024 Fisheries Society of the British Isles.)

The diversity of functional feeding anatomy is particularly impressive in fishes and correlates with various interspecific ecological specializations. Intraspecific polymorphism can manifest in divergent feeding morphology and ecology, often along a benthic-pelagic axis. Arctic charr (Salvelinus alpinus) is a freshwater salmonid known for morphological variation and sympatric polymorphism and in Lake Þingvallavatn, Iceland, four morphs of charr coexist that differ in preferred prey, behaviour, habitat use, and external feeding morphology. We studied variation in six upper and lower jaw bones in adults of these four morphs using geometric morphometrics and univariate statistics. We tested for allometric differences in bone size and shape among morphs, morph effects on bone size and shape, and divergence along the benthic-pelagic axis. We also examined the degree of integration between bone pairs. We found differences in bone size between pelagic and benthic morphs for two bones (dentary and premaxilla). There was clear bone shape divergence along a benthic-pelagic axis in four bones (dentary, articular-angular, premaxilla and maxilla), as well as allometric shape differences between morphs in the dentary. Notably for the dentary, morph explained more shape variation than bone size. Comparatively, benthic morphs possess a compact and taller dentary, with shorter dentary palate, consistent with visible (but less prominent) differences in external morphology. As these morphs emerged in the last 10,000 years, these results indicate rapid functional evolution of specific feeding structures in arctic charr. This sets the stage for studies of the genetics and development of rapid and parallel craniofacial evolution., Competing Interests: The authors have declared that no competing interests exist., (Copyright: © 2024 Jónsdóttir et al. This is an open access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.)

Background: Mitochondrial DNA (mtDNA) has become a significant tool for exploring genetic diversity and delineating evolutionary links across diverse taxa. Within the group of cold-water fish species that are native to the Indian Himalayan region, Schizothorax esocinus holds particular importance due to its ecological significance and is potentially vulnerable to environmental changes. This research aims to clarify the phylogenetic relationships within the Schizothorax genus by utilizing mitochondrial protein-coding genes., Methods: Standard protocols were followed for the isolation of DNA from S. esocinus. For the amplification of mtDNA, overlapping primers were used, and then subsequent sequencing was performed. The genetic features were investigated by the application of bioinformatic approaches. These approaches covered the evaluation of nucleotide composition, codon usage, selective pressure using nonsynonymous substitution /synonymous substitution (Ka/Ks) ratios, and phylogenetic analysis., Results: The study specifically examined the 13 protein-coding genes of Schizothorax species which belongs to the Schizothoracinae subfamily. Nucleotide composition analysis showed a bias towards A + T content, consistent with other cyprinid fish species, suggesting evolutionary conservation. Relative Synonymous Codon Usage highlighted leucine as the most frequent (5.18%) and cysteine as the least frequent (0.78%) codon. The positive AT-skew and the predominantly negative GC-skew indicated the abundance of A and C. Comparative analysis revealed significant conservation of amino acids in multiple genes. The majority of amino acids were hydrophobic rather than polar. The purifying selection was revealed by the genetic distance and Ka/Ks ratios. Phylogenetic study revealed a significant genetic divergence between S. esocinus and other Schizothorax species with interspecific K2P distances ranging from 0.00 to 8.87%, with an average of 5.76%., Conclusion: The present study provides significant contributions to the understanding of mitochondrial genome diversity and genetic evolution mechanisms in Schizothoracinae, hence offering vital insights for the development of conservation initiatives aimed at protecting freshwater fish species., (© 2024. The Author(s), under exclusive licence to Springer Nature B.V.)

Brook charr is a cold-water species which is highly sensitive to increased water temperatures, such as those associated with climate change. Environmental variation can potentially induce phenotypic changes that are inherited across generations, for instance, via epigenetic mechanisms. Here, we tested whether parental thermal regimes (intergenerational plasticity) and offspring-rearing temperatures (within-generational plasticity) modify the brain transcriptome of Brook charr progeny (fry stage). Parents were exposed to either cold or warm temperatures during final gonad maturation and their progeny were reared at 5 or 8 °C during the first stages of development. Illumina Novaseq6000 was used to sequence the brain transcriptome at the yolk sac resorption stage. The number of differentially expressed genes was very low when comparing fry reared at different temperatures (79 differentially expressed genes). In contrast, 9,050 differentially expressed genes were significantly differentially expressed between fry issued from parents exposed to either cold or warm temperatures. There was a significant downregulation of processes related to neural and synaptic activity in fry originating from the warm parental group vs fry from the cold parental one. We also observed significant upregulation of DNA methylation genes and of the most salient processes associated with compensation to warming, such as metabolism, cellular response to stress, and adaptive immunity., Competing Interests: Conflicts of interest The authors declare no conflict of interest., (© The Author(s) 2024. Published by Oxford University Press on behalf of The Genetics Society of America.)

The 'good genes' hypotheses of sexual selection predict that females prefer males with strong ornaments because they are in good health and vigor and can afford the costs of the ornaments. A key assumption of this concept is that male health and vigor are useful predictors of genetic quality and hence offspring performance. We tested this prediction in wild-caught lake char (Salvelinus umbla) whose breeding coloration is known to reveal aspects of male health. We first reanalyzed results from sperm competition trials in which embryos of known parenthood had been raised singly in either a stress- or non-stress environment. Paternal coloration did not correlate with any measures of offspring performance. However, offspring growth was reduced with higher kinship coefficients between the parents. To test the robustness of these first observations, we collected a new sample of wild males and females, used their gametes in a full-factorial in vitro breeding experiment, and singly raised about 3000 embryos in either a stress- or non-stress environment (stress induced by microbes). Again, paternal coloration did not predict offspring performance, while offspring growth was reduced with higher kinship between the parents. We conclude that, in lake char, the genetic benefits of mate choice would be strongest if females could recognize and avoid genetically related males, while male breeding colors may be more relevant in intra-sexual selection., (© 2024. The Author(s).)

Background: A major goal in evolutionary biology is to understand the processes underlying phenotypic variation in nature. Commonly, studies have focused on large interconnected populations or populations found along strong environmental gradients. However, studies on small fragmented populations can give strong insight into evolutionary processes in relation to discrete ecological factors. Evolution in small populations is believed to be dominated by stochastic processes, but recent work shows that small populations can also display adaptive phenotypic variation, through for example plasticity and rapid adaptive evolution. Such evolution takes place even though there are strong signs of historical bottlenecks and genetic drift. Here we studied 24 small populations of the freshwater fish Arctic charr (Salvelinus alpinus) found in groundwater filled lava caves. Those populations were found within a few km2-area with no apparent water connections between them. We studied the relative contribution of neutral versus non-neutral evolutionary processes in shaping phenotypic divergence, by contrasting patterns of phenotypic and neutral genetic divergence across populations in relation to environmental measurements. This allowed us to model the proportion of phenotypic variance explained by the environment, taking in to account the observed neutral genetic structure., Results: These populations originated from the nearby Lake Mývatn, and showed small population sizes with low genetic diversity. Phenotypic variation was mostly correlated with neutral genetic diversity with only a small environmental effect., Conclusions: Phenotypic diversity in these cave populations appears to be largely the product of neutral processes, fitting the classical evolutionary expectations. However, the fact that neutral processes did not explain fully the phenotypic patterns suggests that further studies can increase our understanding on how neutral evolutionary processes can interact with other forces of selection at early stages of divergence. The accessibility of these populations has provided the opportunity for long-term monitoring of individual fish, allowing tracking how the environment can influence phenotypic and genetic divergence for shaping and maintaining diversity in small populations. Such studies are important, especially in freshwater, as habitat alteration is commonly breaking populations into smaller units, which may or may not be viable., (© 2024. The Author(s).)

Temperature variation is affecting fish biodiversity worldwide, causing changes in geographic distribution, phenotypic structure, and even species extinction. Incubation is a critical stage for stenothermic species, which are vulnerable to large temperature fluctuations, and its effects on the phenotype at later developmental stages are understudied, despite the fact that the phenotype being essential for organism ecology and evolution. In this study, we tested the effects of heat shocks during the embryonic period on the phenotype of Arctic charr (Salvelinus alpinus). We repeatedly quantified multiple phenotypic traits, including morphology, development, and behavior, over a period of 4 months, from hatching to juvenile stage in individuals that had experienced heat shocks (+ 5°C on 24 h, seven times) during their embryonic stage and those that had not. We found that heat shocks led to smaller body size at hatching and a lower sociability. Interestingly, these effects weakened throughout the development of individuals and even reversed in the case of body size. We also found an accelerated growth rate and a higher body condition in the presence of heat shocks. Our study provides evidence that heat shocks experienced during incubation can have long-lasting effects on an individual's phenotype. This highlights the importance of the incubation phase for the development of ectothermic organisms and suggests that temperature fluctuations may have significant ecological and evolutionary implications for Arctic charr. Given the predicted increase in extreme events and the unpredictability of temperature fluctuations, it is critical to further investigate their effects on development by examining fluctuations that vary in frequency and intensity., (© 2024 Fisheries Society of the British Isles.)

Across its Holarctic range, Arctic charr (Salvelinus alpinus) populations have diverged into distinct trophic specialists across independent replicate lakes. The major aspect of divergence between ecomorphs is in head shape and body shape, which are ecomorphological traits reflecting niche use. However, whether the genomic underpinnings of these parallel divergences are consistent across replicates was unknown but key for resolving the substrate of parallel evolution. We investigated the genomic basis of head shape and body shape morphology across four benthivore-planktivore ecomorph pairs of Arctic charr in Scotland. Through genome-wide association analyses, we found genomic regions associated with head shape (89 SNPs) or body shape (180 SNPs) separately and 50 of these SNPs were strongly associated with both body and head shape morphology. For each trait separately, only a small number of SNPs were shared across all ecomorph pairs (3 SNPs for head shape and 10 SNPs for body shape). Signs of selection on the associated genomic regions varied across pairs, consistent with evolutionary demography differing considerably across lakes. Using a comprehensive database of salmonid QTLs newly augmented and mapped to a charr genome, we found several of the head- and body-shape-associated SNPs were within or near morphology QTLs from other salmonid species, reflecting a shared genetic basis for these phenotypes across species. Overall, our results demonstrate how parallel ecotype divergences can have both population-specific and deeply shared genomic underpinnings across replicates, influenced by differences in their environments and demographic histories., (© 2024 The Authors. Molecular Ecology published by John Wiley & Sons Ltd.)

Teleost gill mucus has a highly diverse microbiota, which plays an essential role in the host's fitness and is greatly influenced by the environment. Arctic char ( Salvelinus alpinus ), a salmonid well adapted to northern conditions, faces multiple stressors in the Arctic, including water chemistry modifications, that could negatively impact the gill microbiota dynamics related to the host's health. In the context of increasing environmental disturbances, we aimed to characterize the taxonomic distribution of transcriptionally active taxa within the bacterial gill microbiota of Arctic char in the Canadian Arctic in order to identify active bacterial composition that correlates with environmental factors. For this purpose, a total of 140 adult anadromous individuals were collected from rivers, lakes, and bays belonging to five Inuit communities located in four distinct hydrologic basins in the Canadian Arctic (Nunavut and Nunavik) during spring (May) and autumn (August). Various environmental factors were collected, including latitudes, water and air temperatures, oxygen concentration, pH, dissolved organic carbon (DOC), salinity, and chlorophyll-a concentration. The taxonomic distribution of transcriptionally active taxa within the gill microbiota was quantified by 16S rRNA gene transcripts sequencing. The results showed differential bacterial activity between the different geographical locations, explained by latitude, salinity, and, to a lesser extent, air temperature. Network analysis allowed the detection of a potential dysbiosis signature (i.e., bacterial imbalance) in fish gill microbiota from Duquet Lake in the Hudson Strait and the system Five Mile Inlet connected to the Hudson Bay, both showing the lowest alpha diversity and connectivity between taxa.IMPORTANCEThis paper aims to decipher the complex relationship between Arctic char ( Salvelinus alpinus ) and its symbiotic microbial consortium in gills. This salmonid is widespread in the Canadian Arctic and is the main protein and polyunsaturated fatty acids source for Inuit people. The influence of environmental parameters on gill microbiota in wild populations remains poorly understood. However, assessing the Arctic char's active gill bacterial community is essential to look for potential pathogens or dysbiosis that could threaten wild populations. Here, we concluded that Arctic char gill microbiota was mainly influenced by latitude and air temperature, the latter being correlated with water temperature. In addition, a dysbiosis signature detected in gill microbiota was potentially associated with poor fish health status recorded in these disturbed environments. With those results, we hypothesized that rapid climate change and increasing anthropic activities in the Arctic might profoundly disturb Arctic char gill microbiota, affecting their survival., Competing Interests: The authors declare no conflict of interest.

Inbreeding depression, that is, the reduction of health and vigour in individuals with high inbreeding coefficients, is expected to increase with environmental, social, or physiological stress. It has therefore been predicted that sexual selection and the associated stress usually lead to higher inbreeding depression in males than in females. However, sex-specific differences in life history may reverse that pattern during certain developmental stages. In some salmonids, for example, female juveniles start developing their gonads earlier than males who instead grow faster. We tested whether the sexes are differently affected by inbreeding during that time. To study the effects of inbreeding coefficients that may be typical for natural populations of brown trout (Salmo trutta), and also to control for potentially confounding maternal or paternal effects, we sampled males and females from the wild, used their gametes in a block-wise full-factorial breeding design to produce 60 full-sib families, released the offspring as yolk-sac larvae into the wild, sampled them 6 months later, identified their genetic sex, and used microsatellites to assign them to their parents. We used whole-genome resequencing to calculate the kinship coefficients for each breeding pair and hence the expected average inbreeding coefficient per family. Juvenile growth could be predicted from these expected inbreeding coefficients and the genetic sex: Females reached lower body sizes with increasing inbreeding coefficient, while no such link could be found in males. This sex-specific inbreeding depression led to the overall pattern that females were on average smaller than males by the end of their first summer., (© 2024 The Authors. Molecular Ecology published by John Wiley & Sons Ltd.)

International policy recently adopted commitments to maintain genetic diversity in wild populations to secure their adaptive potential, including metrics to monitor temporal trends in genetic diversity - so-called indicators. A national programme for assessing trends in genetic diversity was recently initiated in Sweden. Relating to this effort, we systematically assess contemporary genome-wide temporal trends (40 years) in wild populations using the newly adopted indicators and whole genome sequencing (WGS). We use pooled and individual WGS data from brown trout (Salmo trutta) in eight alpine lakes in protected areas. Observed temporal trends in diversity metrics (nucleotide diversity, Watterson's ϴ and heterozygosity) lie within proposed acceptable threshold values for six of the lakes, but with consistently low values in lakes above the tree line and declines observed in these northern-most lakes. Local effective population size is low in all lakes, highlighting the importance of continued protection of interconnected systems to allow genetic connectivity for long-term viability of these populations. Inbreeding (F ROH ) spans 10%-30% and is mostly represented by ancient (<1 Mb) runs of homozygosity, with observations of little change in mutational load. We also investigate adaptive dynamics over evolutionarily short time frames (a few generations); identifying putative parallel selection across all lakes within a gene pertaining to skin pigmentation as well as candidates of selection unique to specific lakes and lake systems involved in reproduction and immunity. We demonstrate the utility of WGS for systematic monitoring of natural populations, a priority concern if genetic diversity is to be protected., (© 2023 The Authors. Molecular Ecology published by John Wiley & Sons Ltd.)

Although it is known that the whitefish, an ancient salmonid, expresses three distinct gonadotropin-releasing hormone (GnRH) forms in the brain, it has been thought that the later-evolving salmonids (salmon and trout) had only two types of GnRH: GnRH2 and GnRH3. We now provide evidence for the expression of GnRH1 in the gonads of Atlantic salmon by rapid amplification of cDNA ends, real-time quantitative PCR and immunohistochemistry. We examined six different salmonid genomes and found that each assembly has one gene that likely encodes a viable GnRH1 prepropeptide. In contrast to both functional GnRH2 and GnRH3 paralogs, the GnRH1 homeolog can no longer express the hormone. Furthermore, the viable salmonid GnRH1 mRNA is composed of only three exons, rather than the four exons that build the GnRH2 and GnRH3 mRNAs. Transcribed gnrh1 is broadly expressed (in 17/18 tissues examined), with relative abundance highest in the ovaries. Expression of the gnrh2 and gnrh3 mRNAs is more restricted, primarily to the brain, and not in the gonads. The GnRH1 proximal promoter presents composite binding elements that predict interactions with complexes that contain diverse cell fate and differentiation transcription factors. We provide immunological evidence for GnRH1 peptide in the nucleus of 1-year-old type A spermatogonia and cortical alveoli oocytes. GnRH1 peptide was not detected during other germ cell or reproductive stages. GnRH1 activity in the salmonid gonad may occur only during early stages of development and play a key role in a regulatory network that controls mitotic and/or meiotic processes within the germ cell., (© 2023. The Author(s).)

Most population genomic tools rely on accurate single nucleotide polymorphism (SNP) calling and filtering to meet their underlying assumptions. However, genomic complexity, resulting from structural variants, paralogous sequences, and repetitive elements, presents significant challenges in assembling contiguous reference genomes. Consequently, short-read resequencing studies can encounter mismapping issues, leading to SNPs that deviate from Mendelian expected patterns of heterozygosity and allelic ratio. In this study, we employed the ngsParalog software to identify such deviant SNPs in whole-genome sequencing (WGS) data with low (1.5×) to intermediate (4.8×) coverage for four species: Arctic Char (Salvelinus alpinus), Lake Whitefish (Coregonus clupeaformis), Atlantic Salmon (Salmo salar), and the American Eel (Anguilla rostrata). The analyses revealed that deviant SNPs accounted for 22% to 62% of all SNPs in salmonid datasets and approximately 11% in the American Eel dataset. These deviant SNPs were particularly concentrated within repetitive elements and genomic regions that had recently undergone rediploidization in salmonids. Additionally, narrow peaks of elevated coverage were ubiquitous along all four reference genomes, encompassed most deviant SNPs, and could be partially associated with transposons and tandem repeats. Including these deviant SNPs in genomic analyses led to highly distorted site frequency spectra, underestimated pairwise FST values, and overestimated nucleotide diversity. Considering the widespread occurrence of deviant SNPs arising from a variety of sources, their important impact in estimating population parameters, and the availability of effective tools to identify them, we propose that excluding deviant SNPs from WGS datasets is required to improve genomic inferences for a wide range of taxa and sequencing depths., (© The Author(s) 2023. Published by Oxford University Press on behalf of Society for Molecular Biology and Evolution.)

Aquaculture species are often exposed to acute stressors such as low water levels and handling during routine aquaculture procedures. This might result in oxidative stress by the increased reactive oxygen species (ROS)' production ( e.g. , superoxide anion). The harmful effects of ROS are eliminated by a defense system, referred antioxidant defense system (ADS). sod1 is the first gene involved in the ADS. Therefore, we cloned and characterized the open reading frame of the sod1 in brown trout. Then, we determined the effects of low water level and handling stress on sod1 mRNA expression in the liver and gills at 0 min, 1 and 2 h. The total RNA isolated was used to synthesize cDNA for RT-qPCR analysis. Phylogenetic tree, identity/similarity percentages, genomic organization, and conserved gene synteny analyses were applied to characterize Sod1/ sod1 . While low water level stress upregulated sod1 expression in the liver compared to the control group, no significant differences were observed in the gills between experimental groups. However, brown trout differently responded to handling stress at different time intervals in both tissues. Transcriptional differences were also noted between the sexes. This study contributes to the current understanding of the molecular mechanism between oxidative stress and ADS.

Tsinling lenok trout (Brachymystax lenok tsinlingensis) is a typical cold water fish. High temperature has been shown to damage the liver of fish. However, few studies have investigated the liver apoptosis induced by high temperature stress in fish from the perspective of gene expression and metabolic function. Therefore, we investigated the changes caused by high temperature stress (24 °C) on the liver tissue structure, antioxidant capacity, liver gene expression, and the metabolome of tsinling lenok trout. The transcriptomic results showed that genes associated with apoptosis, such as CASP8, CASP3, PERK, Bcl-6 and TRAIL, were upregulated under high temperature stress. Metabolomic analysis showed that the metabolic pathway of nucleotide synthesis was significantly downregulated, while that of oxygen radical synthesis was significantly upregulated. Integrated analysis showed that after high temperature stress, immune-related signaling pathways in trout were activated and their apoptosis level increased, which might be related to hepatopancreas injury. In addition, abnormalities in the tricarboxylic acid cycle and mitochondrial function were observed, suggesting that functional hypoxia caused by high temperature might be involved fish cell apoptosis. These results provide new insights into the process of cell apoptosis in fish under high temperature stress., Competing Interests: Declaration of competing interest The authors declare that they have no conflicts of interest., (Copyright © 2023 Elsevier B.V. All rights reserved.)

Fish often spawn eggs with ovarian fluids that have been hypothesized to support the sperm of some males over others (cryptic female choice). Alternatively, sperm reactions to ovarian fluids could reveal male strategies. We used wild-caught lake char (Salvelinus umbla) to experimentally test whether sperm react differently to the presence of ovarian fluid, and whether any differential sperm reaction could be predicted by male breeding coloration, male inbreeding coefficients (based of 4150 SNPs) or the kinship coefficients between males and females. Male coloration was positively linked to body size and current health (based on lymphocytosis and thrombocytosis) but was a poor predictor of inbreeding or kinship coefficients. We found that sperm of more colourful males were faster in diluted ovarian fluids than in water only, while sperm of paler males were faster in water than in ovarian fluids. We then let equal numbers of sperm compete for fertilizations in the presence or absence of ovarian fluids and genetically assigned 1464 embryos (from 70 experimental trials) to their fathers. The presence of ovarian fluids significantly increased the success of the more colourful competitors. Sperm of less inbred competitors were more successful when tested in water only than in diluted ovarian fluids. The kinship coefficients had no significant effects on sperm traits or fertilization success in the presence of ovarian fluids, although parallel stress tests on embryos had revealed that females would profit more from mating with least related males rather than most coloured ones. We conclude that sperm of more colourful males are best adapted to ovarian fluids, and that the observed reaction norms suggest male strategies rather than cryptic female choice., (© 2023 The Authors. Molecular Ecology published by John Wiley & Sons Ltd.)

Fish stocking constitutes a common management practice in freshwaters all over the world, to enhance fisheries or to support threatened fish populations. Pervasive detrimental effects may affect the real effectiveness of stocking programs. However, studies assessing the real impacts and relative contribution of stocked trout in wild populations are surprisingly few. The marble trout, Salmo marmoratus (Cuvier 1829), is a critically endangered sub-endemic salmonid in Northern Italy, and an iconic species for recreational fishing and conservation, also representing an emblematic case of species negatively affected by restocking. For instance, marble trout inhabiting the Toce River, the second largest tributary of Lake Maggiore, has been stocked with different hatchery congener trout belonging to the Salmo trutta complex (putative marble trout, Atlantic trout Salmo trutta Linnaeus 1758 and putative Mediterranean trout Salmo ghigii Pomini 1941) over the last decades. Using mitochondrial (D-loop) and nuclear (12 microsatellites and LDH-C1*) markers, we characterised the genetic variability and gene flow among the wild and hatchery individuals of marble trout of this basin, to investigate the effectiveness of stocking activities on the native residual population. Despite extensive hybridization of marble trout with non-native brown trout stocks was shown, the presence of individuals belonging to pure native stock has been detected as well. However, concerns could be advanced regarding its future persistence, due to climatic and hydraulic instabilities or loss of environmental heterogeneity. Moreover, despite ongoing yearly massive stocking activities, a negligible contribution of reared putative marble trout in the wild sample has been documented, suggesting that natural recruitment represents the greatest source of this wild population sustainment. Important adaptive differences between wild and domestic trout are present, likely due to the deleterious long-term effects of the close breeding system of hatcheries. Finally, possible implications for stocking management improvement have been discussed., Competing Interests: Declaration of competing interest The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper., (Copyright © 2023 Elsevier B.V. All rights reserved.)

Salmo caspius Kessler, 1877 is one of the most commercially important species of Salmonidae in the southern basin of the Caspian Sea. The occurrence of its wild populations has undergone sever reduction during the last years. In spite of the yearly restocking activity, still no progress on the recovery of its wild population has been observed. Hence, the present study was done in order to assess the efficiency of the current restocking activity in the southern Caspian basin in term of genetic diversity. DNA extracts of 32 S. caspius from four different groups were screened using 62621 genome-wide single nucleotide polymorphisms (SNP). The overal genetic diversity and Fst values were 0.18 and 0.08, respectively. Considering the observed admixture pattern and the positive values for inbreeding coeficient it seems that S. caspius suffers from its small effective population size. In order to obtain the maximum performance, alonside with expanding the size of brood stocks, rehabilitation of the habitats and spawning rivers of this nationally endangered species is of great importance., Competing Interests: The authors have declared that no competing interests exist., (Copyright: © 2023 Najafikhah et al. This is an open access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.)

Eco-evolutionary interactions following ecosystem change provide critical insight into the ability of organisms to adapt to shifting resource landscapes. Here we explore evidence for the rapid parallel evolution of trout feeding morphology following eco-evolutionary interactions with zooplankton in alpine lakes stocked at different points in time in the Wind River Range (Wyoming, USA). In this system, trout predation has altered the zooplankton species community and driven a decrease in average zooplankton size. In some lakes that were stocked decades ago, we find shifts in gill raker traits consistent with the hypothesis that trout have rapidly adapted to exploit available smaller-bodied zooplankton more effectively. We explore this morphological response in multiple lake populations across two species of trout (cutthroat trout, Oncorhynchus clarkii, and golden trout Oncorhynchus aguabonita) and examine the impact of resource availability on morphological variation in gill raker number among lakes. Furthermore, we present genetic data to provide evidence that historically stocked cutthroat trout populations likely derive from multiple population sources, and incorporate variation from genomic relatedness in our exploration of environmental predictors of feeding morphology. These findings describe rapid adaptation and eco-evolutionary interactions in trout and document an evolutionary response to novel, contemporary ecosystem change., (© The Author(s) 2023. Published by Oxford University Press on behalf of The Society for the Study of Evolution (SSE). All rights reserved. For permissions, please e-mail: journals.permissions@oup.com.)

The longest Baltic River, the Vistula, historically hosted numerous sea trout populations. However, dam construction in the twentieth century drastically reduced the spawning migration rate. Reduced natural reproduction has resulted in a population collapse and, consequentially, a substantial reduction in catches. In response, like other Baltic countries, Poland has initiated an intensive stocking program, mainly involving smolt. Initially, stocking was conducted primarily with offspring of sea trout caught during spawning migration. Currently, due to difficulties in obtaining fish, most stocking involves fish from breeding stocks. Therefore, determining the proportion of fish derived from stocking has become an important issue. Experiments based on traditional tagging did not provide sufficient material for analysis; hence, we decided to use genetic methods based on analysis of relatedness. In this study, we performed parentage-based tagging and an analysis of the origins of parent animals used for artificial spawning in 2013, and offspring returning to the Vistula in subsequent years. We based the analysis on three different algorithms and compared the results, showing that the presented methods were effective for estimating mass stocking success. The study also indicated that a certain level of natural reproduction in the Vistula continues to occur. The proportion of sea trout from spawning in 2013 in Vistula sea trout catches from 2017 to 2018 was approximately 30%., (© 2023. The Author(s).)

Background: The classification of the sub-family Schizothoracinae has been debatable due to the overlap in morphological characters. There are discrepancies between classical taxonomy and molecular taxonomy, as well. In the present study, mitochondrial genes CO-I and Cyt b were sequenced to elucidate the phylogenetic status of three species of the genus Schizothorax., Methods and Results: In total, 29 samples of three species viz., S. plagiostomus, S. progastus, and S. richardsonii, were collected from rivers of Uttarakhand, India. For phylogenetic analyses, 40 sequences of CO-I and 41 sequences of Cyt b of Schizothoracinae species were downloaded from NCBI. The highest genetic divergence based on CO-I (16.08%) is between S. plagiostomus and Ptychobarbus dipogon, while the lowest divergence (0.00%) is between 10 pairs of species. The highest divergence based on Cyt b (19.43%), is between S. niger and Gymnocypris eckloni, while the lowest divergence (0.00%) is between four pairs of species. The divergence (0.00% for CO-I and 2.38% for Cyt b) between S. chongi and S. kozlovi, seems a case of convergent molecular evolution of the CO-I gene and in this case, CO-I alone cannot be used to differentiate these two species., Conclusion: The simultaneous use of two molecular markers along with morphomeristic data is a better strategy for the classification of the sub-family Schizothoracinae. These results will be a resource dataset for determining the taxonomical status of Schizothoracine species and will help in the conservation and commercial production of these commercially important fish species., (© 2023. The Author(s), under exclusive licence to Springer Nature B.V.)

The authors used museomics to reconstruct the mitochondrial genome from two individuals of the Moroccan, endemic and extinct trout, Salmo pallaryi. They further obtained partial data from 21 nuclear genes previously used for trout phylogenetic analyses. Phylogenetic analyses, including publicly available data from the mitochondrial control region and the cytochrome b gene, and the 21 nuclear genes, place S. pallaryi among other North African trouts. mtDNA places S. pallaryi close to Salmo macrostigma within a single North African clade. Although the nuclear coverage of the genome was low, both specimens were independently positioned as sisters to one of two distantly related North African clades, viz. the Atlas clade with the Dades trout, Salmo multipunctatus. Phylogenetic discordance between mtDNA and nuclear DNA phylogenies is briefly discussed. As several specimens that were extracted failed to produce DNA of sufficient quality, the authors discuss potential reasons for the failure. They suggest that museum specimens in poor physical condition may be better for DNA extraction compared to better-preserved ones, possibly related to the innovation of formalin as a fixative before ethanol storage in the early 20th century., (© 2023 The Authors. Journal of Fish Biology published by John Wiley & Sons Ltd on behalf of Fisheries Society of the British Isles.)

The European brown trout, Salmo trutta, is a cold-adapted fish reported as a Least Concern species in the IUCN Red List. This species colonized new territories from southern refuges during the last glacial melting, but during the 20th century suffered from anthropic impacts on its habitats. The long-time survival of the species relies on the genetic diversity within and among populations. Brown trout is among the genetically most diverse vertebrate species; however, native populations in Mediterranean rivers have dramatically suffered of introgressive hybridization from extensive releases of evolutionary distant non-native Atlantic stocks. In addition, in Mediterranean rivers climate change will result in unsuitable conditions for the species during the 21st century. Using brown trout populations at the headstreams of a Pyrenean river as a model, this paper revised how hatchery releases have affected the native gene pools and how environmental and climatic variables controlled the amount of local introgression at intra-basin level. Introgressive hybridization was detected in all studied sites. Ten times larger divergence was observed among populations at tributaries than among populations along the main stem. A highly impacted population distributed in a long transect in the main stem suggested that hatchery fish move towards the main stem wherever released. From already highly impacted populations and despite the cessation of hatchery releases, warmer temperatures and lower precipitation expected from climate change will extend the introgressive hybridization along the basin, contributing to the extinction of the native gene pools. Based on available morphological distinction of native, hatchery and hybrid brown trout, we advocate the involvement of regional social groups (e.g. riverside dwellers, anglers, conservationists, hikers) in citizen science programs to detect the spread of non-native phenotypes along the rivers. These are cheap and fast methods to collaborate with fishery managers in the preservation and recovery of the regional native populations., Competing Interests: Declaration of competing interest The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper., (Copyright © 2022 The Authors. Published by Elsevier B.V. All rights reserved.)

Fish growth can be modulated through genetic selection. However, it is not known whether growth regulatory mechanisms modulated by genetic selection can provide information about phenotypic growth variations among families or populations. Following a five-generation breeding program that selected for the absence of early sexual maturity and increased growth in brook charr we aimed to understand how the genetic selection process modifies the growth regulatory pathway of brook charr at the molecular level. To achieve this, we studied the regulation of growth traits at three different levels: 1) between lines-one under selection, the other not, 2) among-families expressing differences in average growth phenotypes, which we termed family performance, and 3) among individuals within families that expressed extreme growth phenotypes, which we termed slow- and fast-growing. At age 1+, individuals from four of the highest performing and four of the lowest performing families in terms of growth were sampled in both the control and selected lines. The gene expression levels of three reference and ten target genes were analyzed by real-time PCR. Results showed that better growth performance (in terms of weight and length at age) in the selected line was associated with an upregulation in the expression of genes involved in the growth hormone (GH)/insulin growth factor-1 (IGF-1) axis, including the igf-1 receptor in pituitary; the gh-1 receptor and igf-1 in liver; and ghr and igf-1r in white muscle. When looking at gene expression within families, family performance and individual phenotypes were associated with upregulations of the leptin receptor and neuropeptid Y-genes related to appetite regulation-in the slower-growing phenotypes. However, other genes related to appetite (ghrelin, somatostatin) or involved in muscle growth (myosin heavy chain, myogenin) were not differentially expressed. This study highlights how transcriptomics may improve our understanding of the roles of different key endocrine steps that regulate physiological performance. Large variations in growth still exist in the selected line, indicating that the full genetic selection potential has not been reached., (Copyright © 2022 Elsevier Inc. All rights reserved.)

The rivers of southern England and northern France which drain into the English Channel contain several genetically unique groups of trout (Salmo trutta L.) that have suffered dramatic declines in numbers over the past 40 years. Knowledge of levels and patterns of genetic diversity is essential for effective management of these vulnerable populations. Using restriction site-associated DNA sequencing (RADseq) data, we describe the development and characterisation of a panel of 95 single nucleotide polymorphism (SNP) loci for trout from this region and investigate their applicability and variability in both target (i.e., southern English) and non-target trout populations from northern Britain and Ireland. In addition, we present three case studies which demonstrate the utility and resolution of these genetic markers at three levels of spatial separation:(a) between closely related populations in nearby rivers, (b) within a catchment and (c) when determining parentage and familial relationships between fish sampled from a single site, using both empirical and simulated data. The SNP loci will be useful for population genetic and assignment studies on brown trout within the UK and beyond., (© 2022 The Authors. Journal of Fish Biology published by John Wiley & Sons Ltd on behalf of Fisheries Society of the British Isles.)

Genetic diversity is the basis for population adaptation and long-term survival, yet rarely considered in biodiversity monitoring. One key issue is the need for useful and straightforward indicators of genetic diversity. We monitored genetic diversity over 40 years (1970-2010) in metapopulations of brown trout (Salmo trutta) inhabiting 27 small mountain lakes representing 10 lake systems in central Sweden using >1200 fish per time point. We tested six newly proposed indicators; three were designed for broad, international use in the UN Convention on Biological Diversity (CBD) and are currently applied in several countries. The other three were recently elaborated for national use by a Swedish science-management effort and applied for the first time here. The Swedish indicators use molecular genetic data to monitor genetic diversity within and between populations (indicators ΔH and ΔF ST , respectively) and assess the effective population size (N e -indicator). We identified 29 genetically distinct populations, all retained over time. Twelve of the 27 lakes harboured more than one population indicating that brown trout biodiversity hidden as cryptic, sympatric populations are more common than recognized. The N e indicator showed values below the threshold (N e  ≤ 500) in 20 populations with five showing N e  < 100. Statistically significant genetic diversity reductions occurred in several populations. Metapopulation structure appears to buffer against diversity loss; applying the indicators to metapopulations suggest mostly acceptable rates of change in all but one system. The CBD indicators agreed with the Swedish ones but provided less detail. All these indicators are appropriate for managers to initiate monitoring of genetic biodiversity., (© 2022 The Authors. Molecular Ecology published by John Wiley & Sons Ltd.)

Salmincola markewitschi Shedko et Shedko, 2002 (Copepoda: Lernaeopodidae) is an ectoparasitic copepod mainly infecting the buccal cavities of white-spotted charr Salvelinus leucomaenis (Pallas) (Salmonidae). This species has only been recorded from Northeast Asia, where a morphologically similar congener Salmincola carpionis (Krøyer, 1837) is also distributed, using the same host species. These copepods are hard to distinguish from each other because of their similarities. We thus examined the newly collected specimens morphologically and genetically from five populations of white-spotted charr in Japan. Most of the specimens were morphologically consistent with S. markewitschi but showed great variations in the numbers of spines on the exopods of the antennae, shape of the maxilliped myxal palps, and the bulla diameter. Consequently, some specimens shared characteristics with S. carpionis. In addition to the mophological continuities, genetic analyses of 28S rDNA and COI mitochondrial DNA confirmed that all specimens belong to a single species. Further taxonomic revisions are required to draw conclusions of whether S. markewitschi is a valid species different from S. carpionis, by collecting samples from across their wide distributional ranges, such as Europe, North America, and Northeast Asia. A key to identification of species of Salmincola Wilson, 1915 occurring in Japan is also provided.

AbstractThe rise in temperature associated with climate change may threaten the persistence of stenothermal organisms with limited capacities for beneficial thermal acclimation. We investigated the capacity for within-generation and transgenerational thermal responses in brook trout ( Salvelinus fontinalis ), a cold-adapted salmonid. Adult fish were acclimated to temperatures within (10°C) and above (21°C) their thermal optimum for 6 mo before spawning, then mated in a full factorial breeding design to produce offspring from cold- and warm-acclimated parents and bidirectional crosses between parents from both temperature treatments. Offspring from families were subdivided and reared at two acclimation temperatures representing their current (15°C) and anticipated future (19°C) habitat temperatures. Offspring thermal physiology was measured as the rate of oxygen consumption (Mo 2 ) during an acute change in temperature (increase of 2°C h -1 ) to observe their Mo 2 -temperature relationship. We recorded resting Mo 2 , peak (highest achieved, thermally induced) Mo 2 , and critical thermal maximum (CTM) as performance metrics. Although limited, within-generation plasticity was greater than transgenerational plasticity, with offspring warm acclimation elevating CTM by 0.5°C but slightly lowering peak thermally induced Mo 2 . Transgenerational plasticity was evident as a slightly elevated resting Mo 2 and a shift of the Mo 2 -temperature relationship to higher rates overall in offspring from warm-acclimated parents. Furthermore, offspring whose parents were warm acclimated were in worse condition than those whose parents were cold acclimated. Both parents contributed to offspring thermal responses; however, the paternal effect was stronger. Despite the existence of within-generation and transgenerational plasticity in brook trout, it is unlikely that these will be sufficient for coping with long-term changes to environmental temperatures.

Effective conservation actions to counteract the current decline of populations and species require a deep knowledge on their genetic structure. We used Single Nucleotide Polymorphisms (SNPs) to infer the population structure of the highly threatened freshwater pearl mussel Margaritifera margaritifera in the Iberian Peninsula. A total of 130 individuals were collected from 26 locations belonging to 16 basins. We obtained 31,692 SNPs through Genotyping by Sequencing (GBS) and used this dataset to infer population structure. Genetic diversity given as observed heterozygosity was low. Pairwise F ST comparisons revealed low levels of genetic differentiation among geographically close populations. Up to 3 major genetic lineages were determined: Atlantic, Cantabrian and Douro. This structure suggests a close co-evolutionary process with brown trout (Salmo trutta), the primordial fish host of this mussel in the studied area. Some sub-basins showed some genetic structuring, whereas in others no intrapopulation differentiation was found. Our results confirm that genetic conservation units do not match individual basins, and that knowledge about the genetic structure is necessary before planning recovery plans that may involve relocation or restocking. The same reasoning should be applied to strictly freshwater species that are sessile or have restricted dispersal abilities and are currently imperiled worldwide., (© 2022. The Author(s).)

Saprolegnia parasitica causes saprolegniosis, a disease responsible for significant economic losses in aquaculture and declines of fish populations in the wild, but the knowledge of its distribution and prevalence in the environment is limited. We developed a fast, sensitive and specific S. parasitica droplet digital PCR (ddPCR) assay and demonstrated its applicability for the detection and quantification of the pathogen in environmental samples: swab DNA collected from the host (trout skin, surface of eggs) and environmental DNA extracted from water. The developed assay was used to assess how abiotic (i.e. physico-chemical parameters of the water) and biotic (health status of the host) factors influence the S. parasitica load in the environment. The pathogen load in water samples was positively correlated with some site-specific abiotic parameters such as electrical conductivity (EC) and calcium, while fluorides were negatively correlated, suggesting that physico-chemical parameters are important for determining S. parasitica load in natural waters. Furthermore, skin swabs of injured trout had significantly higher pathogen load than swabs collected from healthy fish, confirming that S. parasitica is a widespread opportunistic pathogen. Our results provide new insights into various environmental factors that influence the distribution and abundance of S. parasitica., (© 2022. The Author(s).)

Changes in DNA methylation in specific coding or non-coding regions can influence development and potentially divergence in traits within species and groups. While the impact of epigenetic variation on developmental pathways associated with evolutionary divergence is the focus of intense investigation, few studies have looked at recently diverged systems. Phenotypic diversity between closely related populations of Arctic charr (Salvelinus alpinus), which diverged within the last 10,000 years, offers an interesting ecological model to address such effects. Using bisulphite sequencing, we studied general DNA methylation patterns during development in the four sympatric morphs of Arctic charr from Lake Thingvallavatn. The data revealed strong differences between developmental timepoints and between morphs (mainly along the benthic-limnetic axis), both at single CpG sites and in 1000 bp-regions. Genes located close to differentially methylated CpG sites were involved in nucleosome assembly, regulation of osteoclast differentiation, and cell-matrix adhesion. Differentially methylated regions were enriched in tRNA and rRNA sequences, and half of them were located close to transcription start sites. The expression of 14 genes showing methylation differences over time or between morphs was further investigated by qPCR and nine of these were found to be differentially expressed between morphs. Four genes (ARHGEF37-like, H3-like, MPP3 and MEGF9) showed a correlation between methylation and expression. Lastly, histone gene clusters displayed interesting methylation differences between timepoints and morphs, as well as intragenic methylation variation. The results presented here provide a motivation for further studies on the contribution of epigenetic traits, such as DNA methylation, to phenotypic diversity and developmental mechanisms., (© 2022 John Wiley & Sons Ltd.)

Sympatric adaptive phenotypic divergence should be underlain by genomic differentiation between subpopulations. When divergence drives similar patterns of phenotypic and ecological variation within species we expect evolution to draw on common allelic variation. We investigated divergence histories and genomic signatures of adaptive divergence between benthic and pelagic morphs of Icelandic Arctic charr. Divergence histories for each of four populations were reconstructed using coalescent modelling and 14,187 single nucleotide polymorphisms. Sympatric divergence with continuous gene flow was supported in two populations while allopatric divergence with secondary contact was supported in one population; we could not differentiate between demographic models in the fourth population. We detected parallel patterns of phenotypic divergence along benthic-pelagic evolutionary trajectories among populations. Patterns of genomic differentiation between benthic and pelagic morphs were characterized by outlier loci in many narrow peaks of differentiation throughout the genome, which may reflect the eroding effects of gene flow on nearby neutral loci. We then used genome-wide association analyses to relate both phenotypic (body shape and size) and ecological (carbon and nitrogen stable isotopes) variation to patterns of genomic differentiation. Many peaks of genomic differentiation were associated with phenotypic and ecological variation in the three highly divergent populations, suggesting a genomic basis for adaptive divergence. We detected little evidence for a parallel genomic basis of differentiation as most regions and outlier loci were not shared among populations. Our results show that adaptive divergence can have varied genomic consequences in populations with relatively recent common origins, similar divergence histories, and parallel phenotypic divergence., (© 2022 John Wiley & Sons Ltd.)

Unconventional natural gas development (fracking) has been a rapidly expanding technique used for the extraction of natural gas from the Marcellus Shale formation in Pennsylvania. There remains a knowledge gap regarding the ecological impacts of fracking, especially regarding the long-term health of native Brook trout (Salvelinus fontinalis) populations. During the summer of 2015, Brook trout were sampled from twelve streams located in forested, northwestern Pennsylvania in order to evaluate the impacts of fracking on Brook trout. Four stream sites were undisturbed (no fracking activity), three had a developed well pad without fracking activity, and five had active fracking with natural gas production. Liver tissue was isolated from two to five fish per stream and underwent RNA-Seq analysis to identify differentially expressed genes between ecosystems with differing fracking status. Data were analyzed individually and with samples pooled within-stream to account for hierarchical data structure and variation in sample coverage within streams. Differentially expressed and differentially alternatively spliced genes had functions related to lipid and steroid metabolism, mRNA processing, RNA polymerase and protein regulation. Unique to our study, genes related to xenobiotic and stress responses were found as well as potential markers for endocrine disruption and saline adaptation that were identified in watersheds with active fracking activity. These results support the utility of RNA-Seq to assess trout health and suggest detrimental impacts of fracking on sensitive trout populations., (© 2022. The Author(s), under exclusive licence to Springer Science+Business Media, LLC, part of Springer Nature.)

Brown trout ( Salmo trutta ), like many other freshwater species, is threated by the release in its natural environment of alien species and the restocking with allochthonous conspecific stocks. Many conservation projects are ongoing and several morphological and genetic tools have been proposed to support activities aimed to restore genetic integrity status of native populations. Nevertheless, due to the complexity of degree of introgression reached up after many generations of crossing, the use of dichotomous key and molecular markers, such as mtDNA, LDH-C1 * and microsatellites, are often not sufficient to discriminate native and admixed specimens at individual level. Here we propose a reduced panel of ancestry-informative SNP markers (AIMs) to support on field activities for Mediterranean trout management and conservation purpose. Starting from the genotypes data obtained on specimens sampled in the main two Molise's rivers (Central-Southern Italy), a 47 AIMs panel was identified and validated on simulated and real hybrid population datasets, mainly through a Machine Learning approach based on Random Forest classifier. The AIMs panel proposed may represent an interesting and cost-effective tool for monitoring the level of introgression between native and allochthonous trout population for conservation purpose and this methodology could be also applied in other species.

Background: The cnidarian myxozoan parasite Tetracapsuloides bryosalmonae causes chronic proliferative kidney disease (PKD) in salmonids. This parasite is a serious threat to wild and cultured salmonids. T. bryosalmonae undergoes intra-luminal sporogonic development in the kidney of brown trout (Salmo trutta) and the viable spores are released via urine. We investigated the alternative splicing pattern in the posterior kidney of brown trout during PKD., Results: RNA-seq data were generated from the posterior kidney of brown trout collected at 12 weeks post-exposure to T. bryosalmonae. Subsequently, this data was mapped to the brown trout genome. About 153 significant differently expressed alternatively spliced (DEAS) genes, (delta PSI = 5%, FDR P-value < 0.05) were identified from 19,722 alternatively spliced events. Among the DEAS genes, the least and most abundant alternative splicing types were alternative 5' splice site (5.23%) and exon skipping (70.59%), respectively. The DEAS genes were significantly enriched for sodium-potassium transporter activity and ion homeostasis (ahcyl1, atp1a3a, atp1a1a.1, and atp1a1a.5). The protein-protein interaction network analysis enriched two local network clusters namely cation transporting ATPase C-terminus and Sodium/potassium ATPase beta chain cluster, and mixed inclusion of Ion homeostasis and EF-hand domain cluster. Furthermore, the human disease-related salmonella infection pathway was significantly enriched in the protein-protein interaction network., Conclusion: This study provides the first baseline information about alternative splicing in brown trout during PKD. The generated data lay a foundation for further functional molecular studies in PKD - brown trout infection model. The information generated from the present study can help to develop therapeutic strategies for PKD in the future., (© 2022. The Author(s).)

Two related salmonids, Arctic char (Salvelinus alpinus) and lake whitefish (Coregonus clupeaformis) sampled from the high Arctic region of Nunavut, Canada are anadromous fish, migrating annually from the same ice-covered freshwater waterbodies to spend summers in the marine waters of the Arctic Ocean. Microbiota associated with the skin-associated mucus undergo community change coincident with migration, and irrespective of this turnover, antibiotic resistance was detected in mixed bacterial cultures initiated with mucus samples. Although as expected most bacteria were unculturable, however, 5/7 isolates showed susceptibility to a panel of five common antibiotics. The fish were sampled under severe conditions and at remote locations far from human habitation. Regardless, two isolates, 'Carnobacterium maltaromaticum sm-2' and 'Arthrobacter citreus sm', showed multi-resistance to two or more antibiotics including ampicillin and streptomycin indicating multiple resistance genes. It is unknown if these fish bacteria have 'natural' resistance phenotypes or if resistance has been acquired. As result of these observations, we urge long-term monitoring of drug-resistant bacteria in the region and caution the assumption of a lack of drug-resistant organisms even in such extreme environments., (© 2021 The Authors. Environmental Microbiology Reports published by Society for Applied Microbiology and John Wiley & Sons Ltd.)

Stream-dwelling white-spotted charr, Salvelinus leucomaenis , populations tend to have unique color and spot patterns in different regions and may even display stream-specific patterns. An extreme edge of such diversity, found in individuals with atypical body color patterns (so-called nagaremon-type charr, a rare morphotype of Salvelinus leucomaenis [hereafter termed nagaremon-charr]), sympatrically occurring with normal-charr, has been reported from only six small isolated populations in Japan. Based on morphological and ecological perspectives, nagaremon-charr has been considered as an intraspecific color variant of white-spotted charr, although the genetic status of nagaremon-charr has not been determined. In this study, genetic diversity and population structure of the nagaremon-charr in a tributary of the Ane River (Lake Biwa system) were investigated through microsatellite and mtDNA analyses. Nagaremon-charr and sympatric normal-charr in the tributary shared the mtDNA haplotypes and were assigned to the same cluster in the STRUCTURE analysis and discriminant analysis of principal components (DAPC). These results suggested that nagaremon-charr in the Ane River is an intra-populational specific color variant of white-spotted charr. Above a waterfall, nagaremon-charr specimens exhibited extremely reduced genetic diversity, indicating that genetic drift may account for the fixation of the nagaremon-morphotype. Normal-charr below the waterfall clustered separately from hatchery-reared charr, indicative of native status of the former. Thus, both nagaremon-charr and normal-charr in the entire Ane River tributary should be conserved.

Genetic evidence of selection for complex and polygenically regulated phenotypes can easily become masked by neutral population genetic structure and phenotypic plasticity. Without direct evidence of genotype-phenotype associations it can be difficult to conclude to what degree a phenotype is heritable or a product of environment. Common garden laboratory studies control for environmental stochasticity and help to determine the mechanism that regulate traits. Here we assess lipid content, growth, weight, and length variation in full and hybrid F 1 crosses of deep and shallow water sympatric lake charr ecotypes reared for nine years in a common garden experiment. Redundancy analysis (RDA) and quantitative-trait-loci (QTL) genomic scans are used to identify associations between genotypes at 19,714 single nucleotide polymorphisms (SNPs) aligned to the lake charr genome and individual phenotypes to determine the role that genetic inheritance plays in ecotype phenotypic diversity. Lipid content, growth, length, and weight differed significantly among lake charr crosses throughout the experiment suggesting that pedigree plays a large roll in lake charr development. Polygenic scores of 15 SNPs putatively associated with lipid content and/or condition factor indicated that ecotype distinguishing traits are polygenically regulated and additive. A QTL identified on chromosome 38 contained >200 genes, some of which were associated with lipid metabolism and growth, demonstrating the complex nature of ecotype diversity. The results of our common garden study further indicate that lake charr ecotypes observed in nature are predetermined at birth and that ecotypes differ fundamentally in lipid metabolism and growth., (© 2022 The Authors. Molecular Ecology published by John Wiley & Sons Ltd.)

Epigenetic inheritance can result in plastic responses to changing environments being faithfully transmitted to offspring. However, it remains unclear how epigenetic mechanisms such as DNA methylation can contribute to multigenerational acclimation and adaptation to environmental stressors. Brook charr ( Salvelinus fontinalis ), an economically important salmonid, is highly sensitive to thermal stress and is of conservation concern in the context of climate change. We studied the effects of temperature during parental sexual maturation and offspring rearing on whole-genome DNA methylation in brook charr juveniles (fry). Parents were split between warm and cold temperatures during sexual maturation, mated in controlled breeding designs, then offspring from each family were split between warm (8°C) and cold (5°C) rearing environments. Using whole-genome bisulfite sequencing, we found 188 differentially methylated regions (DMRs) due to parental maturation temperature after controlling for family structure. By contrast, offspring rearing temperature had a negligible effect on offspring methylation. Stable intergenerational inheritance of DNA methylation and minimal plasticity in progeny could result in the transmission of acclimatory epigenetic states to offspring, priming them for a warming environment. Our findings have implications pertaining to the role of intergenerational epigenetic inheritance in response to ongoing climate change.

AbstractThe potentially significant genetic consequences associated with the loss of migratory capacity of diadromous fishes that have become landlocked in freshwater are poorly understood. Consistent selective pressures associated with freshwater residency may drive repeated differentiation both between allopatric landlocked and anadromous populations and within landlocked populations (resulting in sympatric morphs). Alternatively, the strong genetic drift anticipated in isolated landlocked populations could hinder consistent adaptation, limiting genetic parallelism. Understanding the degree of genetic parallelism underlying differentiation has implications for both the predictability of evolution and management practices. We employed an 87k single-nucleotide polymorphism (SNP) array to examine the genetic characteristics of landlocked and anadromous Arctic char ( Salvelinus alpinus ) populations from five drainages within Labrador, Canada. One gene was detected as an outlier between sympatric, size-differentiated morphs in each of two landlocked lakes. While no single locus differentiated all replicate pairs of landlocked and anadromous populations, several SNPs, genes, and paralogs were consistently detected as outliers in at least 70% of these pairwise comparisons. A significant C-score suggested that the amount of shared outlier SNPs across all paired landlocked and anadromous populations was greater than expected by chance. Our results indicate that despite their isolation, selection due to the loss of diadromy may drive consistent genetic responses in landlocked populations.

The present experiment tested if temperature during embryogenesis and parental heritage affected the migratory behaviour of young brown trout Salmo trutta. Two parental forms were used, a freshwater resident form and an anadromous form, both from the same river system but geographically isolated since 1993-95. Four groups of young S. trutta were produced and reared from (a) freshwater resident parents spawning in a tributary to the River Imsa, Norway, (b) anadromous parents spawning in the main stem of the same river system, (c) resident male × anadromous female parents and (d) resident female × anadromous male parents. The eggs were incubated until first exogenous feeding in River Imsa water, either unheated or heated c. 2.7°C above ambient temperature. Thereafter, all fish experienced the same ambient river temperature until release. Groups were released below an impassable waterfall 900 m upstream of the mouth of the River Imsa, either as age-0 in October 2019 or as age-1 in May 2020. About 7.5% of the released fish moved downstream and were captured in a trap at the outlet. For any given body size, the proportion of warm incubated trout that moved downstream was greater than the proportion of cold incubated trout. It was also found that most emigrants of the October-released S. trutta were caught within a month of release. Also, most May-released S. trutta emigrated in October. The offspring of the freshwater resident parents emigrated to a larger extent than offspring of anadromous parents. Thus, the difference in emigration with regard to embryonic temperature was phenotypically plastic and may be associated with an epigenetic effect of the thermal conditions during early development. The effect of parental origin suggests there may be genetic divergence between the geographically isolated populations., (© 2022 The Authors. Journal of Fish Biology published by John Wiley & Sons Ltd on behalf of Fisheries Society of the British Isles.)

Animal behaviour is increasingly recognised as critical to the prediction of non-native species success and impacts. Rainbow trout and brown trout have been introduced globally, but there appear to be differences in their patterns of invasiveness and ecological impact. Here, we investigated whether diploid rainbow trout and diploid and triploid brown trout differ among several key behavioural measures linked to invasiveness and impact. We assessed activity, boldness, aggression, and feeding, using open field, novel object, shelter, mirror, feeding, and functional response experiments. We also tested within each fish type for behavioural syndromes comprising correlations among activity, boldness and aggression. Rainbow trout were more active and aggressive but less bold than diploid and triploid brown trout. In small groups, however, rainbow trout were bolder than both types of brown trout. Diploid brown trout were more active and bolder than triploids when tested individually, and had a higher functional response than both rainbow trout and triploid brown trout. In terms of behavioural syndromes, there was no association between activity and boldness in rainbow trout, however, there was in both brown trout types. The increased activity and aggression of rainbow trout may reflect an increased stress response to novel situations, with this response reduced in a group. These results suggest that rainbow trout do not manage their energy budgets effectively, and may explain why they have limited survival as invaders. In addition, the lower functional response of rainbow trout may explain why they are implicated in fewer ecological impacts, and the triploidy treatment also appears to lower the potential impact of brown trout. Comparative analyses of multiple behaviours of invasive species and genetic variants may thus be key to understanding and predicting invader success and ecological impacts., (© 2022. The Author(s).)

Here, we present an annotated, chromosome-anchored, genome assembly for Lake Trout (Salvelinus namaycush) - a highly diverse salmonid species of notable conservation concern and an excellent model for research on adaptation and speciation. We leveraged Pacific Biosciences long-read sequencing, paired-end Illumina sequencing, proximity ligation (Hi-C) sequencing, and a previously published linkage map to produce a highly contiguous assembly composed of 7378 contigs (contig N50 = 1.8 Mb) assigned to 4120 scaffolds (scaffold N50 = 44.975 Mb). Long read sequencing data were generated using DNA from a female double haploid individual. 84.7% of the genome was assigned to 42 chromosome-sized scaffolds and 93.2% of Benchmarking Universal Single Copy Orthologues were recovered, putting this assembly on par with the best currently available salmonid genomes. Estimates of genome size based on k-mer frequency analysis were highly similar to the total size of the finished genome, suggesting that the entirety of the genome was recovered. A mitochondrial genome assembly was also produced. Self-versus-self synteny analysis allowed us to identify homeologs resulting from the salmonid specific autotetraploid event (Ss4R) as well as regions exhibiting delayed rediploidization. Alignment with three other salmonid genomes and the Northern Pike (Esox lucius) genome also allowed us to identify homologous chromosomes in related taxa. We also generated multiple resources useful for future genomic research on Lake Trout, including a repeat library and a sex-averaged recombination map. A novel RNA sequencing data set for liver tissue was also generated in order to produce a publicly available set of annotations for 49,668 genes and pseudogenes. Potential applications of these resources to population genetics and the conservation of native populations are discussed., (© 2021 The Authors. Molecular Ecology Resources published by John Wiley & Sons Ltd.)

Five ecologically and phenotypically divergent ecomorphs of the genus Salmo are known from a landlocked alpine lake in the Caucasus, Lake Sevan. It is an example of sympatric diversification within a species-rich lineage with predominate mode of speciation being allopatric. The diversification of Sevan trouts was accompanied by spawning resource partitioning. Four lacustrine ecomorphs with different temporal-spatial spawning strategies and divergent morphology and coloration evolved along with a fifth ecomorph, brook trout, inhabiting the tributaries. Unfortunately, the Sevan trout diversity was almost destroyed by human activity, with two ecomorphs becoming extinct in the 1980s. We performed reconstruction of the evolutionary history of Sevan trouts based on high-throughput sequencing of both contemporary and historical DNA (∼ 50 y.o.) of all Sevan trout ecomorphs. Our study of complete mitogenomes along with genome-wide SNP data revealed the monophyly of four lacustrine ecomorphs and local brook trout, all derived from the anadromous form Caspian salmon, S. caspius. The species tree suggests a scenario of stepwise evolution from riverine to lacustrine spawning. Three genomic clusters were revealed, of which two refer to the riverine and lacustrine spawners within the flock of Sevan trouts (with F ST value = 0.069). A few SNP outliers under selection were discovered that could be responsible for assortative mating based on visual recognition. The Holocene climatic oscillations and the desiccation of tributaries could have played an important role in the origin of lacustrine spawning. The relationships between lacustrine ecomorphs were not yet fully resolved. This radiation warrants further investigation., (Copyright © 2021 Elsevier Inc. All rights reserved.)