Framing the Salmonidae family phylogenetic portrait: a more complete picture from increased taxon sampling

From sabers to spikes: A newfangled reconstruction of the ancient, giant, sexually dimorphic Pacific salmon, †Oncorhynchus rastrosus (SALMONINAE: SALMONINI)

The impressive †Oncorhynchus rastrosus of the Pacific Northwest's Miocene and Pliocene eras was the largest salmonid ever to live. It sported a hypertrophied premaxilla with a pair of enlarged teeth which the original describers reconstructed as projecting ventrally into the mouth, leading them to assign the species to "Smilodonichthys," a genus now in synonymy. Through CT reconstruction of the holotype and newly collected specimens, we demonstrate that the famed teeth projected laterally like tusks, not ventrally like sabers or fangs. We also expand the original description to characterize sexual dimorphism in mature, breeding individuals. Male and female †Oncorhynchus rastrosus differ in the form of the vomer, rostro-dermethmoid-supraethmoid, and dentary, much as do other extant species of Oncorhynchus. Male specimens possess a more elongate vomer than do females, and female vomers have concave ventral surfaces and prominent median dorsal keels. The dentary of females has no evidence of a kype, though some specimens of †O. rastrosus have a non-uniform density mesial to the tooth bed, which we interpret as a male kype. Unlike extant Oncorhynchus, male and female †O. rastrosus do not differ in premaxilla shape. Because male and females possess hypertrophied premaxillae and lateral premaxillary spikes, the former common name "Sabertoothed Salmon" no longer reflects our understanding of the species' morphology. Accordingly, we redub †O. rastrosus the Spike-Toothed Salmon and postulate that its spikes were multifunctional, serving as defense against predators, in agonism against conspecifics, and as a practical aid to nest construction.

Investigating structural variant, indel and single nucleotide polymorphism differentiation between locally adapted Atlantic salmon populations

Genomic structural variants (SVs) are now recognized as an integral component of intraspecific polymorphism and are known to contribute to evolutionary processes in various organisms. However, they are inherently difficult to detect and genotype from readily available short-read sequencing data, and therefore remain poorly documented in wild populations. Salmonid species displaying strong interpopulation variability in both life history traits and habitat characteristics, such as Atlantic salmon (Salmo salar), offer a prime context for studying adaptive polymorphism, but the contribution of SVs to fine-scale local adaptation has yet to be explored. Here, we performed a comparative analysis of SVs, single nucleotide polymorphisms (SNPs) and small indels (<50 bp) segregating in the Romaine and Puyjalon salmon, two putatively locally adapted populations inhabiting neighboring rivers (Québec, Canada) and showing pronounced variation in life history traits, namely growth, fecundity, and age at maturity and smoltification. We first catalogued polymorphism using a hybrid SV characterization approach pairing both short- (16X) and long-read sequencing (20X) for variant discovery with graph-based genotyping of SVs across 60 salmon genomes, along with characterization of SNPs and small indels from short reads. We thus identified 115,907 SVs, 8,777,832 SNPs and 1,089,321 short indels, with SVs covering 4.8 times more base pairs than SNPs. All three variant types revealed a highly congruent population structure and similar patterns of F ST and density variation along the genome. Finally, we performed outlier detection and redundancy analysis (RDA) to identify variants of interest in the putative local adaptation of Romaine and Puyjalon salmon. Genes located near these variants were enriched for biological processes related to nervous system function, suggesting that observed variation in traits such as age at smoltification could arise from differences in neural development. This study therefore demonstrates the feasibility of large-scale SV characterization and highlights its relevance for salmonid population genomics.

Genetic structure of Ponto-Caspian trout populations shows gene flow among river drainages and supports resident Salmo rizeensis as a genetically distinct taxon

To assess the genetic structure of Ponto-Caspian brown trout (Salmo trutta complex) populations, we analyzed both mitochondrial DNA sequences and genotypes at 10 microsatellite loci of fish caught in the Black Sea and from nine river catchments in Georgia, flowing into either the Black or Caspian seas. The results show that: (1) there is substantial genetic differentiation among Ponto-Caspian trout populations, both among the populations of different nominal species and within those of the same species; (2) the genetic distance between conspecific populations from the Black and Caspian Sea basins exceeds that among the populations within the same basin. Moreover, within drainages, genetic distance correlates with the geographic distance; (3) the Black Sea itself is not a barrier to gene flow among the watersheds draining into the Black Sea; (4) some populations in the headwaters of the rivers draining into the Black Sea Basin fall out of this pattern and likely form a separate, non-anadromous (resident) taxon, previously described from northeastern Turkey as Salmo rizeensis. This hypothesis is supported by mitochondrial DNA phylogeny. The presence of both anadromous and resident populations in a single river basin calls for a substantial re-thinking of speciation patterns and taxonomy of Eurasian brown trout.

Food protein-induced enterocolitis syndrome in response to salmon roe and trout roe

Clinical implications herein, we describe the first case in the medical literature of food protein-induced enterocolitis syndrome in response to specific fish roe.

Ceratonova shasta: a cnidarian parasite of annelids and salmonids

The myxozoan Ceratonova shasta was described from hatchery rainbow trout over 70 years ago. The parasite continues to cause severe disease in salmon and trout, and is recognized as a barrier to salmon recovery in some rivers. This review incorporates changes in our knowledge of the parasite's life cycle, taxonomy and biology and examines how this information has expanded our understanding of the interactions between C. shasta and its salmonid and annelid hosts, and how overarching environmental factors affect this host–parasite system. Development of molecular diagnostic techniques has allowed discrimination of differences in parasite genotypes, which have differing host affinities, and enabled the measurement of the spatio-temporal abundance of these different genotypes. Establishment of the C. shasta life cycle in the laboratory has enabled studies on host–parasite interactions and the availability of transcriptomic data has informed our understanding of parasite virulence factors and host defences. Together, these advances have informed the development of models and management actions to mitigate disease.

Salmonidae Genome: Features, Evolutionary and Phylogenetic Characteristics

The salmon family is one of the most iconic and economically important fish families, primarily possessing meat of excellent taste as well as irreplaceable nutritional and biological value. One of the most common and, therefore, highly significant members of this family, the Atlantic salmon (Salmo salar L.), was not without reason one of the first fish species for which a high-quality reference genome assembly was produced and published. Genomic advancements are becoming increasingly essential in both the genetic enhancement of farmed salmon and the conservation of wild salmon stocks. The salmon genome has also played a significant role in influencing our comprehension of the evolutionary and functional ramifications of the ancestral whole-genome duplication event shared by all Salmonidae species. Here we provide an overview of the current state of research on the genomics and phylogeny of the various most studied subfamilies, genera, and individual salmonid species, focusing on those studies that aim to advance our understanding of salmonid ecology, physiology, and evolution, particularly for the purpose of improving aquaculture production. This review should make potential researchers pay attention to the current state of research on the salmonid genome, which should potentially attract interest in this important problem, and hence the application of new technologies (such as genome editing) in uncovering the genetic and evolutionary features of salmoniforms that underlie functional variation in traits of commercial and scientific importance.

Molecular Phylogeny and Adaptive Mitochondrial DNA Evolution of Salmonids (Pisces: Salmonidae)

Salmonids are composed of anadromous and freshwater fishes, which is an important model for studying adaptive evolution. Herein, 49 salmonid complete mitochondrial genomes and those of two outgroups were used to infer a robust phylogeny for the family Salmonidae. The BI and RAxML phylogenetic trees based on 13 concatenated mitochondrial protein-coding genes showed well-supported nodes, and topologies were highly congruent. The concatenated 13 mitochondrial protein-coding genes, ND2, ND3, and ND5 genes were shown to have significantly larger dN/dS ratios in anadromous species than in freshwater species of Salmonidae, but the CYTB gene had significantly smaller dN/dS in anadromous species. The FEL analysis identified positively selected sites and negatively selected sites in each mitochondrial protein-coding gene separately. The RELAX program revealed that the ATP8 and CYTB genes supported intensified selection of the anadromous lineages. Our results demonstrated the phylogeny of Salmonidae and explored the mitochondrial DNA evolution pattern between anadromous and freshwater salmonids.

Thermal regime during parental sexual maturation, but not during offspring rearing, modulates DNA methylation in brook charr (Salvelinus fontinalis)

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.

The pink salmon genome: Uncovering the genomic consequences of a two-year life cycle

Pink salmon (Oncorhynchus gorbuscha) adults are the smallest of the five Pacific salmon native to the western Pacific Ocean. Pink salmon are also the most abundant of these species and account for a large proportion of the commercial value of the salmon fishery worldwide. A two-year life history of pink salmon generates temporally isolated populations that spawn either in even-years or odd-years. To uncover the influence of this genetic isolation, reference genome assemblies were generated for each year-class and whole genome re-sequencing data was collected from salmon of both year-classes. The salmon were sampled from six Canadian rivers and one Japanese river. At multiple centromeres we identified peaks of Fst between year-classes that were millions of base-pairs long. The largest Fst peak was also associated with a million base-pair chromosomal polymorphism found in the odd-year genome near a centromere. These Fst peaks may be the result of a centromere drive or a combination of reduced recombination and genetic drift, and they could influence speciation. Other regions of the genome influenced by odd-year and even-year temporal isolation and tentatively under selection were mostly associated with genes related to immune function, organ development/maintenance, and behaviour.

Meta-analysis shows environmental contaminants elevate cortisol levels in teleost fish - Effect sizes depend on contaminant class and duration of experimental exposure

Glucocorticoid hormones (GCs) help vertebrates maintain homeostasis during and following challenging events. Short-term elevations in GC levels are necessary for survival, whereas longer-term changes can lead to reduced reproductive output and immunosuppression. Persistent environmental contaminants (ECs) are widespread globally. Experimental exposure of individuals to ECs is associated with varying GC responses, within, and across, species and contaminants. Individuals exposed to ECs over long durations are expected to have prolonged GC elevations, which likely affect their health. We conducted a meta-analysis to test for a relationship between fish GC levels and experimental exposure to ECs, and to explore potential moderators, including duration of exposure, that could help explain the variation in effect sizes within and between studies. We report almost exclusively on cortisol responses of teleost fish to ECs. Although there was much variation in effect sizes, captive-bred fish exposed to ECs had baseline GC levels 1.5× higher than unexposed fish, and fish exposed to pharmaceuticals (estradiols and stimulants being mainly considered) had baseline GC levels approximately 2.5× higher than unexposed fish. We found that captive-bred and wild-caught fish did not differ in GC levels after exposure to the same classes of ECs - studies on captive bred fish may thus enable inferences about GC responses to ECs for wild species. Furthermore, effect sizes did not differ between baseline and challenge-induced GC measures. In different analyses, duration of exposure was negatively correlated to effect size, suggesting that the GC response may acclimate after chronic exposure to some ECs which could potentially alter the GC response of EC-exposed fish to novel stressors. Future studies should explore the effect of multiple stressors on the fish GC response and perform tests on a broader array of contaminant types and vertebrate classes.

Contribution to the trout of Euphrates River, with description of a new species, and range extension of Salmo munzuricus (Salmoniformes, Salmonidae)

In an effort to reveal the Euphrates trout taxonomy, the Karasu River, which is one of the eastern drainages of the river, was investigated and three independent populations were identified. Result revealed that two populations belonged to Salmo munzuricus, which was known only in Munzur River, while the other population belonged to an unnamed species. Salmo baliki, a new species, is described from the Murat River, a drainage of Euphrates River. It differs from Salmo species in adjacent water by the combination of the following characters: a grayish body; commonly one, rarely two pale black spots behind eye and on cheek; two to seven black spots on opercle; a few black spots on back and upper part of flank, missing on predorsal area; few to numerous large irregular-shaped red spots in median, upper and lower part of flank, surrounded by a large irregular-shaped white ring; the number of black and red spots not increasing in parallel with size; maxilla short and narrow; adipose-fin medium size, no or rarely one or two red spot its posterior edge; 107–118 lateral line scales; 24–28 scales rows between dorsal-in origin and lateral line; 18–22 scale rows between lateral line and anal–fin origin; maxilla length 7.7–9.1% SL in males, 8.2–9.6 in females. Finally, the genetic study of the Cyt b mitochondrial gene confirmed the morphological data, suggesting the separation of S. baliki from other Salmo species.

Hierarchical genetic structure and implications for conservation of the world's largest salmonid, Hucho taimen

Population genetic analyses can evaluate how evolutionary processes shape diversity and inform conservation and management of imperiled species. Taimen (Hucho taimen), the world’s largest freshwater salmonid, is threatened, endangered, or extirpated across much of its range due to anthropogenic activity including overfishing and habitat degradation. We generated genetic data using high throughput sequencing of reduced representation libraries for taimen from multiple drainages in Mongolia and Russia. Nucleotide diversity estimates were within the range documented in other salmonids, suggesting moderate diversity despite widespread population declines. Similar to other recent studies, our analyses revealed pronounced differentiation among the Arctic (Selenge) and Pacific (Amur and Tugur) drainages, suggesting historical isolation among these systems. However, we found evidence for finer-scale structure within the Pacific drainages, including unexpected differentiation between tributaries and the mainstem of the Tugur River. Differentiation across the Amur and Tugur basins together with coalescent-based demographic modeling suggests the ancestors of Tugur tributary taimen likely diverged in the eastern Amur basin, prior to eventual colonization of the Tugur basin. Our results suggest the potential for differentiation of taimen at different geographic scales, and suggest more thorough geographic and genomic sampling may be needed to inform conservation and management of this iconic salmonid.

Transcriptome-Wide Comparisons and Virulence Gene Polymorphisms of Host-Associated Genotypes of the Cnidarian Parasite Ceratonova shasta in Salmonids

Ceratonova shasta is an important myxozoan pathogen affecting the health of salmonid fishes in the Pacific Northwest of North America. Ceratonova shasta exists as a complex of host-specific genotypes, some with low to moderate virulence, and one that causes a profound, lethal infection in susceptible hosts. High throughput sequencing methods are powerful tools for discovering the genetic basis of these host/virulence differences, but deep sequencing of myxozoans has been challenging due to extremely fast molecular evolution of this group, yielding strongly divergent sequences that are difficult to identify, and unavoidable host contamination. We designed and optimized different bioinformatic pipelines to address these challenges. We obtained a unique set of comprehensive, host-free myxozoan RNA-seq data from C. shasta genotypes of varying virulence from different salmonid hosts. Analyses of transcriptome-wide genetic distances and maximum likelihood multigene phylogenies elucidated the evolutionary relationship between lineages and demonstrated the limited resolution of the established Internal Transcribed Spacer marker for C. shasta genotype identification, as this marker fails to differentiate between biologically distinct genotype II lineages from coho salmon and rainbow trout. We further analyzed the data sets based on polymorphisms in two gene groups related to virulence: cell migration and proteolytic enzymes including their inhibitors. The developed single-nucleotide polymorphism-calling pipeline identified polymorphisms between genotypes and demonstrated that variations in both motility and protease genes were associated with different levels of virulence of C. shasta in its salmonid hosts. The prospective use of proteolytic enzymes as promising candidates for targeted interventions against myxozoans in aquaculture is discussed. We developed host-free transcriptomes of a myxozoan model organism from strains that exhibited different degrees of virulence, as a unique source of data that will foster functional gene analyses and serve as a base for the development of potential therapeutics for efficient control of these parasites.

Genomic data support management of anadromous Arctic Char fisheries in Nunavik by highlighting neutral and putatively adaptive genetic variation

Distinguishing neutral and adaptive genetic variation is one of the main challenges in investigating processes shaping population structure in the wild, and landscape genomics can help identify signatures of adaptation to contrasting environments. Arctic Char (Salvelinus alpinus) is an anadromous salmonid and the most harvested fish species by Inuit people, including in Nunavik (Québec, Canada), one of the most recently deglaciated regions in the world. Unlike many other anadromous salmonids, Arctic Char occupy coastal habitats near their natal rivers during their short marine phase restricted to the summer ice-free period. Our main objective was to document putatively neutral and adaptive genomic variation in anadromous Arctic Char populations from Nunavik and bordering regions to inform local fisheries management. We used genotyping by sequencing (GBS) to genotype 18,112 filtered single nucleotide polymorphisms (SNP) in 650 individuals from 23 sampling locations along >2000 km of coastline. Our results reveal a hierarchical genetic structure, whereby neighboring hydrographic systems harbor distinct populations grouped by major oceanographic basins: Hudson Bay, Hudson Strait, Ungava Bay, and Labrador Sea. We found genetic diversity and differentiation to be consistent both with the expected postglacial recolonization history and with patterns of isolation-by-distance reflecting contemporary gene flow. Results from three gene-environment association methods supported the hypothesis of local adaptation to both freshwater and marine environments (strongest associations with sea surface and air temperatures during summer and salinity). Our results support a fisheries management strategy at a regional scale, and other implications for hatchery projects and adaptation to climate change are discussed.

A genomic perspective on an old question: Salmo trouts or Salmo trutta (Teleostei: Salmonidae)?

There are particular challenges in defining the taxonomic status of recently radiated groups due to the low level of phylogenetic signal. Members of the Salmo trutta species-complex, which mostly evolved during and following the Pleistocene, show high morphological and ecological diversity that, along with their very wide geographic distribution, have led to morphological description of 47 extant nominal species. However, many of these species have not been supported by previous phylogenetic studies, which could be partly due to lack of significant genetic differences among them, the limited resolution offered by molecular methods previously used, as well as the often local scale of these studies. The development of next-generation sequencing (NGS) and related analytical tools have enhanced our ability to address such challenging questions. In this study, Genotyping-by-Sequencing (GBS) of 15,169 filtered SNPs and mitochondrial DNA (mtDNA) D-loop sequences were combined to assess the phylogenetic relationships among 166 brown trouts representing 21 described species and three undescribed groups collected from 84 localities throughout their natural distribution in Europe, west Asia, and North Africa. The data were analysed using different clustering algorithms (admixture analysis and discriminant analysis of principal components-DAPC), a Bayes Factor Delimitation (BFD) test, species tree reconstruction, gene flow tests (three- and four-population tests), and Rogue taxa identification tests. Genomic contributions of the Atlantic lineage brown trout were found in all major sea basins excluding the North African and Aral Sea basins, suggesting introgressive hybridization of native brown trouts driven by stocking using strains of the Atlantic lineage. After removing the phylogenetic noise caused by the Atlantic brown trout, admixture clusters and DAPC clustering based on GBS data, respectively, resolved 11 and 13 clusters among the previously described brown trout species, which were also supported by BFD test results. Our results suggest that natural hybridization between different brown trout lineages has probably played an important role in the origin of several of the putative species, including S. marmoratus, S. carpio, S. farioides, S. pellegrini, S. caspius (in the Kura River drainage) and Salmo sp. in the Danube River basin. Overall, our results support a multi-species taxonomy for brown trouts. They also resolve some species in the Adriatic-Mediterranean and Black Sea drainages as members of very closely related genomic clusters that may need taxonomic revision. However, any final conclusions pertaining to the taxonomy of the brown trout complex should be based on an integrative approach combining genomic, morphological, and ecological data. To avoid challenges in taxonomy and conservation of species complexes like brown trouts, it is suggested to describe species based on genomic clusters of populations instead of describing species based only on morphologically differentiated single type populations.

A novel reference dated phylogeny for the genus Spodoptera Guenée (Lepidoptera: Noctuidae: Noctuinae): new insights into the evolution of a pest-rich genus

The noctuid genus Spodoptera currently consists of 31 species with varied host plant breadths, ranging from monophagous and oligophagous non-pest species to polyphagous pests of economic importance. Several of these pest species have become major invaders, colonizing multiple continents outside their native range. Such is the case of the infamous fall armyworm, Spodoptera frugiperda (J.E. Smith), which includes two recognized host strains that have not been treated as separate species. Following its accidental introduction to Africa in 2016, it quickly spread through Africa and Asia to Australia. Given that half the described Spodoptera species cause major crop losses, comparative genomics studies of several Spodoptera species have highlighted major adaptive changes in genetic architecture, possibly relating to their pest status. Several recent population genomics studies conducted on two species enable a more refined understanding of their population structures, migration patterns and invasion processes. Despite growing interest in the genus, the taxonomic status of several Spodoptera species remains unstable and evolutionary studies suffer from the absence of a robust and comprehensive dated phylogenetic framework. We generated mitogenomic data for 14 Spodoptera taxa, which are combined with data from 15 noctuoid outgroups to generate a resolved mitogenomic backbone phylogeny using both concatenation and multi-species coalescent approaches. We combine this backbone with additional mitochondrial and nuclear data to improve our understanding of the evolutionary history of the genus. We also carry out comprehensive dating analyses, which implement three distinct calibration strategies based on either primary or secondary fossil calibrations. Our results provide an updated phylogenetic framework for 28 Spodoptera species, identifying two well-supported ecologically diverse clades that are recovered for the first time. Well-studied larvae in each of these clades are characterized by differences in mandibular shape, with one clade's being more specialized on silica-rich C₄ grasses. Interestingly, the inferred timeframe for the genus suggests an earlier origin than previously thought for the genus: about 17-18 million years ago.

Heterogeneous genetic basis of age at maturity in salmonid fishes

Understanding the genetic basis of repeated evolution of the same phenotype across taxa is a fundamental aim in evolutionary biology and has applications in conservation and management. However, the extent to which interspecific life-history trait polymorphisms share evolutionary pathways remains underexplored. Here, we address this gap by studying the genetic basis of a key life-history trait, age at maturity, in four species of Pacific salmonids (genus Oncorhynchus) that exhibit intra- and interspecific variation in this trait-Chinook Salmon, Coho Salmon, Sockeye Salmon, and Steelhead Trout. We tested for associations in all four species between age at maturity and two genome regions, six6 and vgll3, that are strongly associated with the same trait in Atlantic Salmon (Salmo salar). We also conducted a genome-wide association analysis in Steelhead to assess whether additional regions were associated with this trait. We found the genetic basis of age at maturity to be heterogeneous across salmonid species. Significant associations between six6 and age at maturity were observed in two of the four species, Sockeye and Steelhead, with the association in Steelhead being particularly strong in both sexes (p = 4.46 × 10^-9 after adjusting for genomic inflation). However, no significant associations were detected between age at maturity and the vgll3 genome region in any of the species, despite its strong association with the same trait in Atlantic Salmon. We discuss possible explanations for the heterogeneous nature of the genetic architecture of this key life-history trait, as well as the implications of our findings for conservation and management.

Effects of Maternally Transferred Egg Selenium on Embryo-Larval Survival, Growth, and Development in Arctic Grayling (Thymallus arcticus)

Selenium (Se) toxicity to fish is primarily manifested via maternal transfer to the eggs, which may result in adverse effects on larval survival and development. The present study assessed the effects of egg Se concentrations derived via maternal transfer on early life-stage development, survival, and growth of Arctic grayling (Thymallus arcticus), a salmonid species not previously assessed for Se sensitivity. Fish gametes were collected from 4 streams in Alaska known to exhibit a range of egg Se concentrations. Eggs were fertilized and reared in the laboratory from hatch through post-swim-up. Larvae were assessed for survival, length, and weight, as well as deformities (skeletal, craniofacial, fin-fold) and edema based on a graduated severity index. Eggs from a total of 47 females were collected, with egg Se concentrations ranging from 3.3 to 33.9 mg kg^-1  dry weight. No relationships were observed between larval endpoints evaluated and parent females' egg, muscle, or whole-body Se concentrations. Therefore, Se 10% effective concentrations (EC10s) were defined as the maximum measured Se concentrations: >33.9, >17.6, and >19.7 mg kg^-1  dry weight for eggs, muscle, and whole-body tissue, respectively. Collectively, these data indicate that Arctic grayling are relatively insensitive to maternally transferred Se compared to other fish species. Environ Toxicol Chem 2021;40:380-389. © 2020 SETAC.

MHC class I evolution; from Northern pike to salmonids

Background

Salmonids are of major importance both as farmed and wild animals. With the changing environment comes changes in pathogenic pressures so understanding the immune system of all salmonid species is of essence. Major histocompatibility complex (MHC) genes are key players in the adaptive immune system signalling infection to responding T-cells populations. Classical MHC class I (MHCI) genes, defined by high polymorphism, broad expression patterns and peptide binding ability, have a key role in inducing immunity. In salmonids, the fourth whole genome duplication that occurred 94 million years ago has provided salmonids with duplicate MHCI regions, while Northern Pike, a basal sister clade to salmonids, represent a species which has not experienced this whole genome duplication.

Results

Comparing the gene organization and evolution of MHC class I gene sequences in Northern pike versus salmonids displays a complex picture of how many of these genes evolved. Regional salmonid Ia and Ib Z lineage gene duplicates are not orthologs to the Northern pike Z lineage sequences. Instead, salmonids have experienced unique gene duplications in both duplicate regions as well as in the Salmo and Oncorhynchus branch. Species-specific gene duplications are even more pronounced for some L lineage genes.

Conclusions

Although both Northern pike as well as salmonids have expanded their U and Z lineage genes, these gene duplications occurred separately in pike and in salmonids. However, the similarity between these duplications suggest the transposable machinery was present in a common ancestor. The salmonid MHCIa and MHCIb regions were formed during the 94 MYA since the split from pike and before the Oncorhynchus and Salmo branch separated. As seen in tetrapods, the non-classical U lineage genes are diversified duplicates of their classical counterpart. One MHCI lineage, the L lineage, experienced massive species-specific gene duplications after Oncorhynchus and Salmo split approximately 25 MYA. Based on what we currently know about L lineage genes, this large variation in number of L lineage genes also signals a large functional diversity in salmonids.

Differential migration in Pacific salmon and trout: Patterns and hypotheses

Migrations affect the population dynamics, life history, evolution, and connections of animals to natural ecosystems and humans. Many species and populations display partial migration (some individuals migrate and some do not), and differential migration (migration distance varies). Partial migration is widely distributed in fishes but the term differential migration is much less commonly applied, despite the occurrence of this phenomenon. This paper briefly reviews the extent of differential migration in Pacific salmon and trout (genus Oncorhynchus), a very extensively studied group. Three hypotheses are presented to explain the patterns among species: 1) phylogenetic relationships, 2) the prevalence of partial migration (i.e., variation in anadromy), and 3) life history patterns (iteroparous or semelparous, and duration spent feeding at sea prior to maturation). Each hypothesis has some support but none is consistent with all patterns. The prevalence of differential migration, ranging from essentially non-existent to common within a species, reflects phylogeny and life history, interacting with the geographic features of the region where juvenile salmon enter the ocean. Notwithstanding the uncertain evolution of this behavior, it has very clear implications for salmon conservation, as it strongly affects exposure to predators, patterns of fishery exploitation and also uptake of toxic contaminants.

Global systematic diversity, range distributions, conservation and taxonomic assessments of graylings (Teleostei: Salmonidae; Thymallus spp.)

Graylings (Thymallus) are among the less well-studied groups of salmonid fishes, especially across their Asian distribution range. Here we perform a comprehensive global review of their phylogeography, systematic diversity and range distributions, including biogeographic reconstruction and assessment of both conservation and taxonomic status of each species. Based on a mitogenomic phylogenetic analysis, three approaches to the delineation of molecular operational units, and evaluation of 15 a-priori defined species, we provide biological support for the recognition of 13 grayling species, plus two additional species tentatively. Several instances of paraphyly and its potential effect on systematic inferences are discussed. Overall, the genus displays increasing species diversity and decreasing range size from higher to lower latitudes and ancestral trait reconstruction supports an East Asian origin for extant diversity, most likely centred in the Amur River drainage. Europe’s colonization by Thymallus took place as early as the late Miocene, at least two colonisations of North America are supported, and multiple dispersal events likely took place into Western Siberia. The conservation status for the 15 taxa was estimated to be: 6 least concern, 1 near-threatened, 2 vulnerable, 3 endangered and 3 data deficient.

Fate of MHCII in salmonids following 4WGD

Major histocompatibility complex (MHC) genes are key players in the adaptive immunity providing a defense against invading pathogens. Although the basic structures are similar when comparing mammalian and teleost MHC class II (MHCII) molecules, there are also clear-cut differences. Based on structural requirements, the teleosts non-classical MHCII molecules do not comply with a function similar to the human HLA-DM and HLA-DO, i.e., assisting in peptide loading and editing of classical MHCII molecules. We have previously studied the evolution of teleost class II genes identifying various lineages and tracing their phylogenetic occurrence back to ancient ray-finned fishes. We found no syntenic MHCII regions shared between cyprinids, salmonids, and neoteleosts, suggesting regional instabilities. Salmonids have experienced a unique whole genome duplication 94 million years ago, providing them with the opportunity to experiment with gene duplicates. Many salmonid genomes have recently become available, and here we set out to investigate how MHCII has evolved in salmonids using Northern pike as a diploid sister phyla, that split from the salmonid lineage prior to the fourth whole genome duplication (4WGD) event. We identified 120 MHCII genes in pike and salmonids, ranging from 11 to 20 genes per species analyzed where DB-group genes had the most expansions. Comparing the MHC of Northern pike with that of Atlantic salmon and other salmonids species provides a tale of gene loss, translocations, and genome rearrangements.

Rapture facilitates inexpensive and high-throughput parent-based tagging in salmonids

Accurate methods for tracking individuals are crucial to the success of fisheries and aquaculture management. Management of migratory salmonid populations, which are important for the health of many economies, ecosystems, and indigenous cultures, is particularly dependent on data gathered from tagged fish. However, the physical tagging methods currently used have many challenges including cost, variable marker retention, and information limited to tagged individuals. Genetic tracking methods combat many of the problems associated with physical tags, but have their own challenges including high cost, potentially difficult marker design, and incompatibility of markers across species. Here we show the feasibility of a new genotyping method for parent-based tagging (PBT), where individuals are tracked through the inherent genetic relationships with their parents. We found that Rapture sequencing, a combination of restriction-site associated DNA and capture sequencing, provides sufficient data for parentage assignment. Additionally, the same capture bait set, which targets specific restriction-site associated DNA loci, can be used for both Rainbow Trout Oncorhynchus mykiss and Chinook Salmon Oncorhynchus tshawytscha. We input 248 single nucleotide polymorphisms from 1,121 samples to parentage assignment software and compared parent-offspring relationships of the spawning pairs recorded in a hatchery. Interestingly, our results suggest sperm contamination during hatchery spawning occurred in the production of 14% of offspring, further confirming the need for genetic tagging in accurately tracking individuals. PBT with Rapture successfully assigned progeny to parents with a 98.86% accuracy with sufficient genetic data. Cost for this pilot study was approximately $3 USD per sample. As costs vary based on the number of markers used and individuals sequenced, we expect that when implemented at a large-scale, per sample costs could be further decreased. We conclude that Rapture PBT provides a cost-effective and accurate alternative to the physical coded wire tags, and other genetic-based methods.

The sockeye salmon genome, transcriptome, and analyses identifying population defining regions of the genome

Sockeye salmon (Oncorhynchus nerka) is a commercially and culturally important species to the people that live along the northern Pacific Ocean coast. There are two main sockeye salmon ecotypes—the ocean-going (anadromous) ecotype and the fresh-water ecotype known as kokanee. The goal of this study was to better understand the population structure of sockeye salmon and identify possible genomic differences among populations and between the two ecotypes. In pursuit of this goal, we generated the first reference sockeye salmon genome assembly and an RNA-seq transcriptome data set to better annotate features of the assembly. Resequenced whole-genomes of 140 sockeye salmon and kokanee were analyzed to understand population structure and identify genomic differences between ecotypes. Three distinct geographic and genetic groups were identified from analyses of the resequencing data. Nucleotide variants in an immunoglobulin heavy chain variable gene cluster on chromosome 26 were found to differentiate the northwestern group from the southern and upper Columbia River groups. Several candidate genes were found to be associated with the kokanee ecotype. Many of these genes were related to ammonia tolerance or vision. Finally, the sex chromosomes of this species were better characterized, and an alternative sex-determination mechanism was identified in a subset of upper Columbia River kokanee.

Adaptive and maladaptive genetic diversity in small populations: Insights from the Brook Charr (Salvelinus fontinalis) case study

Investigating the relative importance of neutral versus selective processes governing the accumulation of genetic variants is a key goal in both evolutionary and conservation biology. This is particularly true in the context of small populations, where genetic drift can counteract the effect of selection. Using Brook Charr (Salvelinus fontinalis) from Québec, Canada, as a case study, we investigated the importance of demographic versus selective processes governing the accumulation of both adaptive and maladaptive mutations in closed versus open and connected populations to assess gene flow effect. This was achieved by using 14,779 high-quality filtered SNPs genotyped among 1,416 fish representing 50 populations from three life history types: lacustrine (closed populations), riverine and anadromous (connected populations). Using the PROVEAN algorithm, we observed a considerable accumulation of putative deleterious mutations across populations. The absence of correlation between the occurrence of putatively beneficial or deleterious mutations and local recombination rate supports the hypothesis that genetic drift might be the main driver of the accumulation of such variants. However, despite a lower genetic diversity observed in lacustrine than in riverine or anadromous populations, lacustrine populations do not exhibit more deleterious mutations than the two other history types, suggesting that the negative effect of genetic drift in lacustrine populations may be mitigated by that of relaxed purifying selection. Moreover, we also identified genomic regions associated with anadromy, as well as an overrepresentation of transposable elements associated with variation in environmental variables, thus supporting the importance of transposable elements in adaptation.

Comparative Genomic Analyses and a Novel Linkage Map for Cisco (Coregonus artedi) Provide Insights into Chromosomal Evolution and Rediploidization Across Salmonids

Whole-genome duplication (WGD) is hypothesized to be an important evolutionary mechanism that can facilitate adaptation and speciation. Genomes that exist in states of both diploidy and residual tetraploidy are of particular interest, as mechanisms that maintain the ploidy mosaic after WGD may provide important insights into evolutionary processes. The Salmonidae family exhibits residual tetraploidy, and this, combined with the evolutionary diversity formed after an ancestral autotetraploidization event, makes this group a useful study system. In this study, we generate a novel linkage map for cisco (Coregonus artedi), an economically and culturally important fish in North America and a member of the subfamily Coregoninae, which previously lacked a high-density haploid linkage map. We also conduct comparative genomic analyses to refine our understanding of chromosomal fusion/fission history across salmonids. To facilitate this comparative approach, we use the naming strategy of protokaryotype identifiers (PKs) to associate duplicated chromosomes to their putative ancestral state. The female linkage map for cisco contains 20,292 loci, 3,225 of which are likely within residually tetraploid regions. Comparative genomic analyses revealed that patterns of residual tetrasomy are generally conserved across species, although interspecific variation persists. To determine the broad-scale retention of residual tetrasomy across the salmonids, we analyze sequence similarity of currently available genomes and find evidence of residual tetrasomy in seven of the eight chromosomes that have been previously hypothesized to show this pattern. This interspecific variation in extent of rediploidization may have important implications for understanding salmonid evolutionary histories and informing future conservation efforts.

Addressing incomplete lineage sorting and paralogy in the inference of uncertain salmonid phylogenetic relationships

Recent and continued progress in the scale and sophistication of phylogenetic research has yielded substantial advances in knowledge of the tree of life; however, segments of that tree remain unresolved and continue to produce contradicting or unstable results. These poorly resolved relationships may be the product of methodological shortcomings or of an evolutionary history that did not generate the signal traits needed for its eventual reconstruction. Relationships within the euteleost fish family Salmonidae have proven challenging to resolve in molecular phylogenetics studies in part due to ancestral autopolyploidy contributing to conflicting gene trees. We examine a sequence capture dataset from salmonids and use alternative strategies to accommodate the effects of gene tree conflict based on aspects of salmonid genome history and the multispecies coalescent. We investigate in detail three uncertain relationships: (1) subfamily branching, (2) monophyly of Coregonus and (3) placement of Parahucho. Coregoninae and Thymallinae are resolved as sister taxa, although conflicting topologies are found across analytical strategies. We find inconsistent and generally low support for the monophyly of Coregonus, including in results of analyses with the most extensive dataset and complex model. The most consistent placement of Parahucho is as sister lineage of Salmo.

Evaluating the effect of reference genome divergence on the analysis of empirical RADseq datasets

The advent of high-throughput sequencing (HTS) has made genomic-level analyses feasible for nonmodel organisms. A critical step of many HTS pipelines involves aligning reads to a reference genome to identify variants. Despite recent initiatives, only a fraction of species has publically available reference genomes. Therefore, a common practice is to align reads to the genome of an organism related to the target species; however, this could affect read alignment and bias genotyping. In this study, I conducted an experiment using empirical RADseq datasets generated for two species of salmonids (Actinopterygii; Teleostei; Salmonidae) to address these questions. There are currently reference genomes for six salmonids of varying phylogenetic distance. I aligned the RADseq data to all six genomes and identified variants with several different genotypers, which were then fed into population genetic analyses. Increasing phylogenetic distance between target species and reference genome reduced the proportion of reads that successfully aligned and mapping quality. Reference genome also influenced the number of SNPs that were generated and depth at those SNPs, although the affect varied by genotyper. Inferences of population structure were mixed: increasing reference genome divergence reduced estimates of differentiation but similar patterns of population relationships were found across scenarios. These findings reveal how the choice of reference genome can influence the output of bioinformatic pipelines. It also emphasizes the need to identify best practices and guidelines for the burgeoning field of biodiversity genomics.

Mapping of Adaptive Traits Enabled by a High-Density Linkage Map for Lake Trout

Understanding the genomic basis of adaptative intraspecific phenotypic variation is a central goal in conservation genetics and evolutionary biology. Lake trout (Salvelinus namaycush) are an excellent species for addressing the genetic basis for adaptive variation because they express a striking degree of ecophenotypic variation across their range; however, necessary genomic resources are lacking. Here we utilize recently-developed analytical methods and sequencing technologies to (1) construct a high-density linkage and centromere map for lake trout, (2) identify loci underlying variation in traits that differentiate lake trout ecophenotypes and populations, (3) determine the location of the lake trout sex determination locus, and (4) identify chromosomal homologies between lake trout and other salmonids of varying divergence. The resulting linkage map contains 15,740 single nucleotide polymorphisms (SNPs) mapped to 42 linkage groups, likely representing the 42 lake trout chromosomes. Female and male linkage group lengths ranged from 43.07 to 134.64 centimorgans, and 1.97 to 92.87 centimorgans, respectively. We improved the map by determining coordinates for 41 of 42 centromeres, resulting in a map with 8 metacentric chromosomes and 34 acrocentric or telocentric chromosomes. We use the map to localize the sex determination locus and multiple quantitative trait loci (QTL) associated with intraspecific phenotypic divergence including traits related to growth and body condition, patterns of skin pigmentation, and two composite geomorphometric variables quantifying body shape. Two QTL for the presence of vermiculations and spots mapped with high certainty to an arm of linkage group Sna3, growth related traits mapped to two QTL on linkage groups Sna1 and Sna12, and putative body shape QTL were detected on six separate linkage groups. The sex determination locus was mapped to Sna4 with high confidence. Synteny analysis revealed that lake trout and congener Arctic char (Salvelinus alpinus) are likely differentiated by three or four chromosomal fissions, possibly one chromosomal fusion, and 6 or more large inversions. Combining centromere mapping information with putative inversion coordinates revealed that the majority of detected inversions differentiating lake trout from other salmonids are pericentric and located on acrocentric and telocentric linkage groups. Our results suggest that speciation and adaptive divergence within the genus Salvelinus may have been associated with multiple pericentric inversions occurring primarily on acrocentric and telocentric chromosomes. The linkage map presented here will be a critical resource for advancing conservation oriented genomic research on lake trout and exploring chromosomal evolution within and between salmonid species.

Parallel selection on ecologically relevant gene functions in the transcriptomes of highly diversifying salmonids

Background

Salmonid fishes are characterised by a very high level of variation in trophic, ecological, physiological, and life history adaptations. Some salmonid taxa show exceptional potential for fast, within-lake diversification into morphologically and ecologically distinct variants, often in parallel; these are the lake-resident charr and whitefish (several species in the genera Salvelinus and Coregonus). To identify selection on genes and gene categories associated with such predictable diversifications, we analysed 2702 orthogroups (4.82 Mbp total; average 4.77 genes/orthogroup; average 1783 bp/orthogroup). We did so in two charr and two whitefish species and compared to five other salmonid lineages, which do not evolve in such ecologically predictable ways, and one non-salmonid outgroup.

Results

All selection analyses are based on Coregonus and Salvelinus compared to non-diversifying taxa. We found more orthogroups were affected by relaxed selection than intensified selection. Of those, 122 were under significant relaxed selection, with trends of an overrepresentation of serine family amino acid metabolism and transcriptional regulation, and significant enrichment of behaviour-associated gene functions. Seventy-eight orthogroups were under significant intensified selection and were enriched for signalling process and transcriptional regulation gene ontology terms and actin filament and lipid metabolism gene sets. Ninety-two orthogroups were under diversifying/positive selection. These were enriched for signal transduction, transmembrane transport, and pyruvate metabolism gene ontology terms and often contained genes involved in transcriptional regulation and development. Several orthogroups showed signs of multiple types of selection. For example, orthogroups under relaxed and diversifying selection contained genes such as ap1m2, involved in immunity and development, and slc6a8, playing an important role in muscle and brain creatine uptake. Orthogroups under intensified and diversifying selection were also found, such as genes syn3, with a role in neural processes, and ctsk, involved in bone remodelling.

Conclusions

Our approach pinpointed relevant genomic targets by distinguishing among different kinds of selection. We found that relaxed, intensified, and diversifying selection affect orthogroups and gene functions of ecological relevance in salmonids. Because they were found consistently and robustly across charr and whitefish and not other salmonid lineages, we propose these genes have a potential role in the replicated ecological diversifications.

Enhanced osmoregulatory ability marks the smoltification period in developing chum salmon (Oncorhynchus keta)

The freshwater (FW) life of chum salmon is short, as they migrate to the ocean soon after emergence from the substrate gravel of natal waters. The alevins achieve seawater (SW) acclimating ability at an early developmental stage and the details of smoltification are not clear. We examined the stage-dependent SW acclimating ability in chum salmon alevins and found a sharp increase in SW tolerance during development that resembles the physiological parr-smolt transformation seen in other salmonids. Perturbation of plasma Na⁺ after SW exposure was prominent from the hatched embryo stage to emerged alevins, but the plasma Na⁺ became highly stable and more resistant to perturbation soon after complete absorption of yolk. Marker gene expression for SW-ionocytes including Na/K-ATPase (NKA α1b), Na-K-Cl cotransporter 1a (NKCC1a), Na/H exchanger 3a (NHE3a), cystic fibrosis transmembrane conductance regulators (CFTR I and CFTR II) were all upregulated profoundly at the same stage when the alevins were challenged by SW, suggesting that the stability of plasma Na⁺ concentration was partly a result of elevated osmoregulatory capability. FW-ionocyte markers including NKA α1a and NHE3b were consistently downregulated independent of stage by SW exposure, suggesting that embryos at all stages respond to salinity challenge, but the increase in SW osmoregulatory capability is restricted to the developmental stage after emergence. We propose that the "smoltification period" is condensed and integrated into the early development of chum salmon, and our results can be extrapolated to the future studies on hormonal controls and developmental triggers for smoltification in salmonids.

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.

First report of three complete mitochondrial genomes of the long-finned charr Salvethymus svetovidovi Chereshnev et Skopetz, 1990 (Salmoniformes: Salmonidae) with phylogenetic consideration

The complete mitochondrial genome was sequenced in three individuals of long-finned charr Salvethymus svetovidovi from Lake El’gygytgyn (Chukotka Peninsula, Russia). The genome sequences are 16,655 bp in size and the gene arrangement, composition, and size are similar to the charr genomes published previously. The difference between the three genomes studied is low, 0.07%. Our results support the phylogenetic closeness of Sl. svetovidovi with representatives of the genus Salvelinus and their origin from a common ancestor. A placement of Sl. svetovidovi in the phylogenetic tree is strictly defined and this taxon should be included in the genus Salvelinus.

Cloning, characterization and tissue-specific expression of the antimicrobial peptide hepcidin from caspian trout (Salmo caspius) and the antibacterial activity of the synthetic peptide

Antimicrobial peptides have a wide range of antimicrobial activity and widely occur in different organisms including mollusks, crustaceans and vertebrates. Hepcidins are a group of cysteine-rich antimicrobial peptides that are active against a variety of pathogens including gram-positive and gram-negative bacteria, as well as viruses. In this study, the hepcidin gene of Caspian trout (CtHep) was identified and characterized. Our results showed that CtHep cDNA has a 267-bp Open Reading Frame (ORF), which is translated to 88 amino acids. The CtHep was classified in the HAMP1 class of hepcidins. Comparison of DNA and cDNA sequences showed that CtHep has 3 exons and 2 introns. The signal, prodomain and mature part of CtHep have 24, 39 and 25 amino acids, respectively. The mature peptide has a molecular weight of 2881.43 Da and a theoretical isoelectric point of 8.53. The expression of CtHep mRNA was detected in different tissues of healthy and infected fish. CtHep expression in the liver, head kidney, spleen and skin was significantly enhanced after bacterial challenge. Expression of CtHep in different embryonic development stages was also substantial. Antibacterial activity of synthetic CtHep peptides was investigated against a number of Gram-positive and Gram-negative bacteria. CtHep inhibited some pathogenic bacteria such as Streptococcus iniae and Aeromonas hydrophila. In the in vivo experiment, CtHep upregulated the cytokines IL-6 and TNF-α in both kidney and spleen tissues after 24 h of the peptide injection. In conclusion, our study showed that CtHep plays an important role in the immune system of Caspian trout and also in the embryonic stages. Moreover, CtHep peptide has a potential to be used as an antimicrobial therapeutic agent as well as an immunostimulant in aquaculture.

Expression and activity of lipid and oxidative metabolism enzymes following elevated temperature exposure and thyroid hormone manipulation in juvenile lake whitefish (Coregonus clupeaformis)

Temperature has unequivocal effects on several aspects of fish physiology, but the full extent of its interaction with key endocrine signaling systems to influence metabolic function remains unknown. The aim of the current study was to assess the individual and combined effects of elevated temperature and hyperthyroidism on hepatic metabolism in juvenile lake whitefish by quantifying mRNA abundance and activity of key metabolic enzymes. Fish were exposed to 13 (control), 17 or 21 °C for 0, 4, 8 or 24 days in the presence or absence of low-T₄ (1 µg × g body weight^-1) or high-T₄ (10 µg × g body weight^-1) treatment. Our results demonstrate moderate sensitivity to elevated temperature in this species, characterized by short-term changes in mRNA abundance of several metabolic enzymes and long-term declines in citrate synthase (CS) and cytochrome c oxidase (COX) activities. T₄-induced hyperthyroidism also had several short-term effects on mRNA abundance of metabolic transcripts, including depressions in acetyl-coA carboxylase β (accβ) and carnitine palmitoyltransferase 1β (cpt1β), and stabilization of cs mRNA levels; however, these effects were primarily limited to elevated temperature groups, indicating temperature-dependent effects of exogenous T₄ treatment in this species. In contrast, maximal CS and COX activities were not altered by hyperthyroidism at any temperature. Collectively, our data suggest that temperature has the potential to manipulate thyroid hormone physiology in juvenile lake whitefish and, under warm-conditions, hyperthyroidism may suppress certain elements of the β-oxidation pathway without substantial impacts on overall cellular oxidative capacity.

Sex-Specific Co-expression Networks and Sex-Biased Gene Expression in the Salmonid Brook Charr Salvelinus fontinalis

Networks of co-expressed genes produce complex phenotypes associated with functional novelty. Sex differences in gene expression levels or in the structure of gene co-expression networks can cause sexual dimorphism and may resolve sexually antagonistic selection. Here we used RNA-sequencing in the salmonid Brook Charr Salvelinus fontinalis to characterize sex-specific co-expression networks in the liver of 47 female and 53 male offspring. In both networks, modules were characterized for functional enrichment, hub gene identification, and associations with 15 growth, reproduction, and stress-related phenotypes. Modules were then evaluated for preservation in the opposite sex, and in the congener Arctic Charr Salvelinus alpinus Overall, more transcripts were assigned to a module in the female network than in the male network, which coincided with higher inter-individual gene expression and phenotype variation in the females. Most modules were preserved between sexes and species, including those involved in conserved cellular processes (e.g., translation, immune pathways). However, two sex-specific male modules were identified, and these may contribute to sexual dimorphism. To compare with the network analysis, differentially expressed transcripts were identified between the sexes, revealing a total of 16% of expressed transcripts as sex-biased. For both sexes, there was no overrepresentation of sex-biased genes or sex-specific modules on the putative sex chromosome. Sex-biased transcripts were also not overrepresented in sex-specific modules, and in fact highly male-biased transcripts were enriched in preserved modules. Comparative network analysis and differential expression analyses identified different aspects of sex differences in gene expression, and both provided new insights on the genes underlying sexual dimorphism in the salmonid Brook Charr.

Integration and modularity of teleostean pectoral fin shape and its role in the diversification of acanthomorph fishes

Phenotypic integration and modularity describe the strength and pattern of interdependencies between traits. Integration and modularity have been proposed to influence the trajectory of evolution, either acting as constraints or facilitators. Here, we examine trends in the integration and modularity of pectoral fin morphology in teleost fishes using geometric morphometrics. We compare the fin shapes of the highly diverse radiation of acanthomorph fishes to lower teleosts. Integration and modularity are measured using two-block partial least squares analysis and the covariance ratio coefficient between the radial bones and lepidotrichia of the pectoral fins. We show that the fins of acanthomorph fishes are more tightly integrated but also more morphologically diverse and faster evolving compared to nonacanthomorph fishes. The main pattern of shape covariation in nonacanthomorphs is concordant with the main trajectory of evolution between nonacanthomorphs and acanthomorphs. Our findings support a facilitating role for integration during the acanthomorph diversification. Potential functional consequences and developmental mechanisms of fin integration are discussed.

Interspecific germ cell transplantation: a new light in the conservation of valuable Balkan trout genetic resources?

Interspecific transplantation of germ cells from the brown trout Salmo trutta m. fario and the European grayling Thymallus thymallus into rainbow trout Oncorhynchus mykiss recipients was carried out in order to improve current practices in conservation of genetic resources of endangered salmonid species in the Balkan Peninsula. Current conservation methods mainly include in situ efforts such as the maintenance of purebred individuals in isolated streams and restocking with purebred fingerlings; however, additional ex situ strategies such as surrogate production are needed. Steps required for transplantation such as isolation of high number of viable germ cells and fluorescent labeling of germ cells which are to be transplanted have been optimized. Isolated and labeled brown trout and grayling germ cells were intraperitoneally transplanted into 3 to 5 days post hatch rainbow trout larvae. Survival of the injected larvae was comparable to the controls. Sixty days after transplantation, fluorescently labeled donor cells were detected within the recipient gonads indicating successful incorporation of germ cells (brown trout spermatogonia and oogonia-27%; grayling spermatogonia-28%; grayling oogonia-23%). PCR amplification of donor mtDNA CR fragments within the recipient gonads additionally corroborated the success of incorporation. Overall, the transplantation method demonstrated in this study presents the first step and a possible onset of the application of the germ cell transplantation technology in conservation and revitalization of genetic resources of endangered and endemic species or populations of salmonid fish and thus give rise to new or improved management strategies for such species.

Assessing the Performance of Ks Plots for Detecting Ancient Whole Genome Duplications

Genomic data have provided evidence of previously unknown ancient whole genome duplications (WGDs) and highlighted the role of WGDs in the evolution of many eukaryotic lineages. Ancient WGDs often are detected by examining distributions of synonymous substitutions per site (Ks) within a genome, or “Ks plots.” For example, WGDs can be detected from Ks plots by using univariate mixture models to identify peaks in Ks distributions. We performed gene family simulation experiments to evaluate the effects of different Ks estimation methods and mixture models on our ability to detect ancient WGDs from Ks plots. The simulation experiments, which accounted for variation in substitution rates and gene duplication and loss rates across gene families, tested the effects of WGD age and gene retention rates following WGD on inferring WGDs from Ks plots. Our simulations reveal limitations of Ks plot analyses. Strict interpretations of mixture model analyses often overestimate the number of WGD events, and Ks plot analyses typically fail to detect WGDs when ≤10% of the duplicated genes are retained following the WGD. However, WGDs can accurately be characterized over an intermediate range of Ks. The simulation results are supported by empirical analyses of transcriptomic data, which also suggest that biases in gene retention likely affect our ability to detect ancient WGDs. Although our results indicate mixture model results should be interpreted with great caution, using node-averaged Ks estimates and applying more appropriate mixture models can improve the accuracy of detecting WGDs.

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

Candidate genes associated with migration have been identified in multiple taxa: including salmonids, many of whom perform migrations requiring a series of physiological changes associated with the freshwater–saltwater transition. We screened over 5,500 SNPs for signatures of selection related to migratory behavior of brown trout Salmo trutta by focusing on ten differentially migrating freshwater populations from two watersheds (the Koutajoki and the Oulujoki). We found eight outlier SNPs potentially associated with migratory versus resident life history using multiple (≥3) outlier detection approaches. Comparison of three migratory versus resident population pairs in the Koutajoki watershed revealed seven outlier SNPs, of which three mapped close to genes ZNF665-like, GRM4-like, and PCDH8-like that have been previously associated with migration and smoltification in salmonids. Two outlier SNPs mapped to genes involved in mucus secretion (ST3GAL1-like) and osmoregulation (C14orf37-like). The last two strongly supported outlier SNPs mapped to thermally induced genes (FNTA1-like, FAM134C-like). Within the Oulujoki, the only consistent outlier SNP mapped close to a gene (EZH2) that is associated with compensatory growth in fasted trout. Our results suggest that a relatively small yet common set of genes responsible for physiological functions associated with resident and migratory life histories is evolutionarily conserved.

Hybridization between two bitterling fish species in their sympatric range and a river where one species is native and the other is introduced

The distributions of two bitterling fish (subfamily: Acheilognathinae), Tanakia lanceolata and T. limbata, overlap in western Japan. Acheilognathinae fish lay their eggs in the gills of freshwater bivalves, and the early juvenile stage develops in the gills. Populations of freshwater bivalves are declining worldwide, which has limited the number of spawning substrate for bitterlings. T. limbata has been artificially introduced to some rivers in Ehime, Japan, where it coexists with native T. lanceolata, and some hybrids have been observed. We collected both species from several sites in western Japan, and from the Kunichi River system in Ehime, and analyzed genetic population structure based on six microsatellite loci and sequences of the mitochondrial cytochrome b gene. Structure analysis identified three genetically distinct populations: T. lanceolata, T. limbata “West Kyushu”, and T. limbata “Setouchi”. Two clades of T. limbata were also supported by molecular phylogenetic analyses based on cytochrome b. Hybrids in Ehime originated mostly from interbreeding between male T. lanceolata and female T. limbata “West Kyushu”, and made up 10.2% of all collected fish, suggesting that hybrids occurred frequently between females of colonizing species and males of native species. On the other hand, interspecific hybrids were detected at rates of 40.0%, 20.0%, and 17.6% in the Ima River (Fukuoka), Midori River (Kumamoto), and Kase River (Saga), respectively, which are naturally sympatric regions. We found a few T. limbata “Setouchi” in the Midori and Kase Rivers, which were supposed to be introduced from other regions, coexisting with native T. limbata “West Kyushu”, and this cryptic invasion may have triggered the interspecific hybridization. These results suggest that artificial introduction of a fish species, a decline in the unionid population, and degradation of habitat have caused broad hybridization of bitterlings in western Japan.

Population genetic analysis of aquaculture salmonid populations in China using a 57K rainbow trout SNP array

Various salmonid species are cultivated in cold water aquaculture. However, due to limited genomic data resources, specific high-throughput genotyping tools are not available to many of the salmonid species. In this study, a 57K single nucleotide polymorphism (SNP) array for rainbow trout (Oncorhynchus mykiss) was utilized to detect polymorphisms in seven salmonid species, including Hucho taimen, Oncorhynchus masou, Salvelinus fontinalis, Brachymystax lenok, Salvelinus leucomaenis, O. kisutch, and O. mykiss. The number of polymorphic markers per population ranged from 3,844 (O. kisutch) to 53,734 (O. mykiss), indicating that the rainbow trout SNP array was applicable as a universal genotyping tool for other salmonid species. Among the six other salmonid populations from four genera, 28,882 SNPs were shared, whereas 525 SNPs were polymorphic in all four genera. The genetic diversity and population relationships of the seven salmonid species were studied by principal component analysis (PCA). The phylogenetic relationships among populations were analyzed using the maximum likelihood method, which indicated that the shared SNP markers provide reliable genomic information for population genetic analyses in common aquaculture salmonid fishes. Furthermore, this obtained genomic information may be applicable for population genetic evaluation, marker-assisted breeding, and propagative parent selection in fry production.