Molecular systematics of Salmonidae: combined nuclear data yields a robust phylogeny

The genome assembly of the westslope cutthroat trout, Oncorhynchus lewisi, reveals interspecific chromosomal rearrangements with the rainbow trout, Oncorhynchus mykiss

Cutthroat trout (Oncorhynchus clarkii) are popular among anglers throughout their native range along the West Coast and interior of North America. As they colonized the interior of North America, cutthroat trout diverged into several genetically distinct groups. Many of these groups are now threatened by habitat destruction, hybridization with rainbow trout (Oncorhynchus mykiss), and competition from introduced species. These groups were previously classified as subspecies, but recent research suggests that they may represent distinct species. In this study, we produced a chromosomal-level genome assembly and a genetic map for one of the species in the cutthroat trout species complex, the westslope cutthroat trout (Oncorhynchus lewisi-formerly Oncorhynchus clarkii lewisi). We also constructed haplotype-resolved assemblies from a westslope cutthroat-rainbow trout F1 hybrid. We used the new genome assemblies to identify major interspecific chromosomal rearrangements between the 2 sister species, including fusions, fissions, and inversions. These genome assemblies and chromosome data provide valuable insights regarding genetic variation within cutthroat trout and in hybrids between rainbow and cutthroat trout.

Characterization of six new complete mitochondrial genomes of Chiasmodontidae (Scombriformes, Percomorpha) and considerations about the phylogenetic relationships of the family

The fishes of the Chiasmodontidae family, known as swallower fishes, are species adapted to live in deep seas. Several studies have shown the proximity of this family to Tetragonuridae and Amarsipidae. However, the phylogenetic position of this clade related to other Pelagiaria groups remains uncertain even when phylogenomic studies are employed. Since the low number of published mitogenomes, our study aimed to assemble six new mitochondrial genomes of Chiasmodontidae from database libraries to expand the discussion regarding the phylogeny of this group within Scombriformes. As expected, the composition and organization of mitogenomes were stable among the analyzed species, although we detected repetitive sequences in the D-loop of species of the genus Kali not seen in Chiasmodon, Dysalotus, and Pseudoscopelus. Our phylogeny incorporating 51 mitogenomes from several families of Scombriformes, including nine chiasmodontids, recovered interfamilial relationships well established in previous studies, including a clade containing Chiasmodontidae, Amarsipidae, and Tetragonuridae. However, phylogenetic relationships between larger clades remain unclear, with disagreements between different phylogenomic studies. We argue that such inconsistencies are not only due to biases and limitations in the data but mainly to complex biological events in the adaptive irradiation of Scombriformes after the Cretaceous-Paleogene extinction event.

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.

Convergent evolution of polyploid genomes from across the eukaryotic tree of life

By modeling the homoeologous gene losses that occurred in 50 genomes deriving from ten distinct polyploidy events, we show that the evolutionary forces acting on polyploids are remarkably similar, regardless of whether they occur in flowering plants, ciliates, fishes, or yeasts. We show that many of the events show a relative rate of duplicate gene loss before the first postpolyploidy speciation that is significantly higher than in later phases of their evolution. The relatively weak selective constraint experienced by the single-copy genes these losses produced leads us to suggest that most of the purely selectively neutral duplicate gene losses occur in the immediate postpolyploid period. Nearly all of the events show strong evidence of biases in the duplicate losses, consistent with them being allopolyploidies, with 2 distinct progenitors contributing to the modern species. We also find ongoing and extensive reciprocal gene losses (alternative losses of duplicated ancestral genes) between these genomes. With the exception of a handful of closely related taxa, all of these polyploid organisms are separated from each other by tens to thousands of reciprocal gene losses. As a result, it is very unlikely that viable diploid hybrid species could form between these taxa, since matings between such hybrids would tend to produce offspring lacking essential genes. It is, therefore, possible that the relatively high frequency of recurrent polyploidies in some lineages may be due to the ability of new polyploidies to bypass reciprocal gene loss barriers.

Molecular characterization of transient receptor potential vanilloid 4 (TRPV4) gene transcript variant mRNA of chum salmon Oncorhynchus keta in response to salinity or temperature changes

Transient receptor potential vanilloid 4 (TRPV4) is an osmosensory cation channel that respond to an increase in cell volume and participates in various physiological functions. Among organisms in aquatic environments, euryhaline teleost is are suitable experimental models to study ion channel proteins related to physiological functions involving osmosensing. Among the studies of various regulatory molecules that mediate osmotic regulation in fish, however, information is lacking, particularly on the TRP family. This study investigated the structural characteristics of theTRPV4 gene of chum salmon (Oncorhynchus keta) and their responses to changes in salinity and temperature. Interestingly, TRPV4 generates transcript variants of the intron-retention form through alternative splicing, resulting in a frameshift leading to the generation of transcripts of different structures. In particular, TRPV4 x1 and TRPV x2 mRNAs were predominant in the gill and skin including at the lateral line. The expression levels of chum salmon TRPV4 x1 were significantly increased with increase in salinity and temperature, whereas TRPV4 x2 mainly responded to temperature decrease. Overall, these results demonstrate for the first time the effects of salinity and temperature on the expression of two salmonid TRPV4 transcript variants, suggesting their contribution to the regulation of hydromineral balance.

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.

Comparative transcriptome profiling of selected osmotic regulatory proteins in the gill during seawater acclimation of chum salmon (Oncorhynchus keta) fry

Salmonid fishes, chum salmon (Oncorhynchus keta) have the developed adaptive strategy to withstand wide salinity changes from the early life stage. This study investigated gene expression patterns of cell membrane proteins in the gill of chum salmon fry on the transcriptome level by tracking the salinity acclimation of the fish in changing environments ranging from freshwater (0 ppt) to brackish water (17.5 ppt) to seawater (35 ppt). Using GO analysis of DEGs, the known osmoregulatory genes and their functional groups such as ion transport, transmembrane transporter activity and metal ion binding were identified. The expression patterns of membrane protein genes, including pump-mediated protein (NKA, CFTR), carrier-mediated protein (NKCC, NHE3) and channel-mediated protein (AQP) were similar to those of other salmonid fishes in the smolt or adult stages. Based on the protein-protein interaction analysis between transmembrane proteins and other related genes, we identified osmotic-related genes expressed with salinity changes and analyzed their expression patterns. The findings of this study may facilitate the disentangling of the genetic basis of chum salmon and better able an understanding of the osmophysiology of the species.

Toxicity in aquatic model species exposed to a temporal series of three different flowback and produced water samples collected from a horizontal hydraulically fractured well

In the present study, we compared the toxicity and associated chemical characterizations of flowback and produced water (FPW) collected from a single horizontal hydraulically fractured well at different time points during FPW production. Since few studies on whole mixture toxicity related to FPW exist, our aims were to determine both overall toxicity of the FPW mixture in a suite of organisms (Daphnia magna, Lumbriculus variegatus, Danio rerio, and Oncorhynchus mykiss) and also determine if toxicity changes depending on variation in FPW chemical properties as a function of time sampled (1.33, 72, and 228 h FPW samples collected immediately post-well production onset were analyzed in current study). FPW chemical composition was determined via quadra-pole inductively coupled plasma - mass spectrometry/mass spectrometry (ICP-MS/MS), full-scan high performance liquid chromatography/Orbitrap mass spectrometry (HPLC/Orbitrap-MS), and gas chromatography-mass spectrometry (GC-MS). We observed that FPW sampled later in the production process contained higher ion and total dissolved solids concentrations, whereas the highest concentrations of dissolved organic compounds were observed in the earliest FPW sample analyzed. Toxicity associated with FPW exposure was deemed to be species-specific to a certain extent, but general trends revealed the earliest FPW sampled contained highest toxic potential. Accordingly, we theorize that although the saline conditions of FPW are the foremost toxicological drivers to freshwater organisms, dissolved organics associated with FPW significantly contribute to the overall toxicity of exposed organisms.

Iteroparity of stream resident masu salmon Oncorhynchus masou

The degree of iteroparity in stream-resident forms of masu salmon Oncorhynchus masou was examined using mark-recapture studies in natural streams. In a partially migratory population, at least 10% of resident males survived after maturation and repeatedly matured for up to 5 years. In the landlocked amago salmon subspecies, the post-maturation survival rate was at least 7% and repeat maturation was observed for up to 3 years for both males and females.

The genetic status of the Hungarian brown trout populations: exploration of a blind spot on the European map of Salmo trutta studies

BACKGROUND

Analyses of the control region sequences of European brown trout populations' mitrochondrial DNA have revealed five main evolutionary lineages (Atlantic, Danubian, Mediterranean, Adriatic, Marble) mostly relating to the main water basins; however, the hybridization between lineages were increasingly reported. Due to the hydrogeography of Hungary, wild populations should theoretically belong to the Danubian lineage, however, this has not been verified by genetic studies.

METHODS

In our study multiple molecular marker sets (mitochondrial sequence, microsatellites, PCR-RFLP of nuclear markers and sex marker) were used to investigate the genetic composition and population genetics of the brown trout populations in two broodstocks, six wild streams in Hungary and one Serbian population.

RESULTS

The admixture of Atlantic and Danubian lineages in these populations, except the Serbian population with pure Danubian origin, was observed by control region sequences of mitochondrial DNA and PCR-RFLP markers in the nuclear genome, and one unpublished Danubian haplotype was found in Hungarian populations. A sex-specific marker revealed equal gender ratio in broodstocks and Kemence stream, whereas in other wild streams the proportion of female individuals were less than 50%. Structure and principal component analyses based on the alleles of microsatellite loci also revealed overlapping populations, however the populations were still significantly different from each other and were mostly in Hardy-Weinberg equilibrium.

DISCUSSION

Stocking and migration can have a significant genetic impact on trout populations of wild streams, however there are no guidelines or common practices for stocking of small streams in Hungary, thus the genetic background of these populations should be considered when developing conservation actions.

A Narrowing of the Phenotypic Diversity Range after Large Rearrangements of the Karyotype in Salmonidae: The Relationship between Saltational Genome Rearrangements and Gradual Adaptive Evolution

The problem of how a gradual development of ecological and morphological adaptations combines with large genome rearrangements, which have been found to occur in the phylogeny of many groups of organisms, is a matter of discussion in the literature. The objective of this work was to study the problem with the example of salmonids, whose evolution included at least six events of multiple chromosome fusions. Large karyotype rearrangements are associated with a decrease in ecological and morphological diversity in salmonids. In the above example, genome rearrangements seem to distort the function of the genetic systems that are responsible for the occurrence of certain ecological forms in salmonids.

The genetic relationship between extirpated and contemporary Atlantic salmon Salmo salar L. lines from the southern Baltic Sea

Background

The genetic relationship between original Atlantic salmon populations that are now extinct in the southern Baltic Sea and the present-day populations has long been controversial. To investigate and clarify this issue, we successfully genotyped individuals of the historical populations from the Oder and Vistula Rivers using DNA extracted from dried scales with the Atlantic salmon single nucleotide polymorphism array.

Results

Our results showed a global F_ST of 0.2515 for all pairs of loci, which indicates a high level of genetic differentiation among the groups analyzed in this study. Pairwise F_ST values were significant for all comparisons and the highest values were found between present-day reintroduced Slupia River salmon and extinct Vistula River Atlantic salmon. Bayesian analysis of genetic structure revealed the existence of substructures in the extirpated Polish populations and three main clades among studied stocks.

Conclusions

The historical salmon population from the Oder River was genetically closer to present-day salmon from the Neman River than to the historical salmon from the Vistula River. Vistula salmon clearly separated from all other analyzed salmon stocks. It is likely that the origins of the Atlantic salmon population from the Morrum River and the Polish historical native populations are different.

Electronic supplementary material

The online version of this article (doi:10.1186/s12711-016-0208-y) contains supplementary material, which is available to authorized users.

Assessment of yellowtail kingfish (Seriola lalandi) as a surrogate host for the production of southern bluefin tuna (Thunnus maccoyii) seed via spermatogonial germ cell transplantation

Germ cell transplantation is an innovative technology for the production of interspecies surrogates, capable of facilitating easier and more economical management of large-bodied broodstock, such as the bluefin tuna. The present study explored the suitability of yellowtail kingfish (Seriola lalandi) as a surrogate host for transplanted southern bluefin tuna (Thunnus maccoyii) spermatogonial cells to produce tuna donor-derived gametes upon sexual maturity. Germ cell populations in testes of donor T. maccoyii males were described using basic histology and the molecular markers vasa and dead-end genes. The peripheral area of the testis was found to contain the highest proportions of dead-end-expressing transplantable Type A spermatogonia. T. maccoyii Type A spermatogonia-enriched preparations were transplanted into the coelomic cavity of 6–10-day-old post-hatch S. lalandi larvae. Fluorescence microscopy and polymerase chain reaction analysis detected the presence of tuna cells in the gonads of the transplanted kingfish fingerlings at 18, 28, 39 and 75 days after transplantation, indicating that the transplanted cells migrated to the genital ridge and had colonised the developing gonad. T. maccoyii germ cell-derived DNA or RNA was not detected at later stages, suggesting that the donor cells were not maintained in the hosts’ gonads.

Patterns of variation and covariation in the shapes of mandibular bones of juvenile salmonids in the genus Oncorhynchus

What is the nature of evolutionary divergence of the jaw skeleton within the genus Oncorhynchus? How can two associated bones evolve new shapes and still maintain functional integration? Here, we introduce and test a "concordance" hypothesis, in which an extraordinary matching of the evolutionary shape changes of the dentary and angular articular serves to preserve their fitting together. To test this hypothesis, we examined morphologies of the dentary and angular articular at parr (juvenile) stage, and at three levels of biological organization—between salmon and trout, between sister species within both salmon and trout, and among three types differing in life histories within one species, Oncorhynchus mykiss. The comparisons show bone shape divergences among the groups at each level; morphological divergence between salmon and trout is marked even at this relatively early life history stage. We observed substantial matching between the two mandibular bones in both pattern and amount of shape variation, and in shape covariation across species. These findings strongly support the concordance hypothesis, and reflect functional and/or developmental constraint on morphological evolution. We present evidence for developmental modularity within both bones. The locations of module boundaries were predicted from the patterns of evolutionary divergences, and for the dentary, at least, would appear to facilitate its functional association with the angular articular. The modularity results suggest that development has biased the course of evolution.

Patterns of natural selection acting on the mitochondrial genome of a locally adapted fish species

Background

Mitochondrial DNA (mtDNA) is frequently used in population genetic studies and is usually considered as a neutral marker. However, given the functional importance of the proteins encoded by the mitochondrial genome, and the prominent role of mitochondria in cellular energy production, the assumption of neutrality is increasingly being questioned.

Results

We tested for evidence of selection on the mitochondrial genome of the Atlantic salmon, which is a locally adapted and widely farmed species and is distributed across a large latitudinal cline. We analysed 20 independent regions of the salmon mtDNA that represented nine genes (ND1, ND2, ND3, COX1, COX2, ATP6, ND4, ND5, and CYTB). These 20 mtDNA regions were sequenced using a 454 approach from samples collected across the entire European range of this species. We found evidence of positive selection at the ND1, ND3 and ND4 genes, which is supported by at least two different codon-based methods and also by differences in the chemical properties of the amino acids involved. The geographical distribution of some of the mutations indicated to be under selection was not random, and some mutations were private to artic populations. We discuss the possibility that selection acting on the Atlantic salmon mtDNA genome might be related to the need for increased metabolic efficiency at low temperatures.

Conclusions

The analysis of sequences representing nine mitochondrial genes that are involved in the OXPHOS pathway revealed signatures of positive selection in the mitochondrial genome of the Atlantic salmon. The properties of the amino acids involved suggest that some of the mutations that were identified to be under positive selection might have functional implications, possibly in relation to metabolic efficiency. Experimental evidence, and better understanding of regional phylogeographic structuring, are needed to clarify the potential role of selection acting on the mitochondrial genome of Atlantic salmon and other locally adapted fishes.

Electronic supplementary material

The online version of this article (doi:10.1186/s12711-015-0138-0) contains supplementary material, which is available to authorized users.

Genetic analysis of four European huchen (Hucho hucho Linnaeus, 1758) broodstocks from Poland, Germany, Slovakia, and Ukraine: implication for conservation

Four broodstocks of European huchen (Hucho hucho) from: Poland, Germany, Slovakia, and Ukraine were investigated using ten microsatellite DNA loci. Microsatellite DNA analysis was successfully applied for the first time in the Polish broodstock of this fish species. The genetic variation and genetic distance between these broodstocks were evaluated. In addition, we examined the potential effects of a genetic bottleneck on the genetic variation of the broodstocks. The European huchen broodstocks exhibited moderate genetic diversity (PIC = 0.405–0.496 and I = 0.831–1.047) with the exception of German broodstock which presented higher genetic diversity (PIC = 0.590 and I = 1.254). Observed (Ho) and expected (He) heterozygosity across the investigated loci in all broodstocks ranged from 0.434 to 0.686 and from 0.452 to 0.650, respectively. Overall, the studied broodstocks were in Hardy-Weinberg equilibrium (HWE); however, from 8 to 42 % of the loci deviated from HWE in each stock. The Garza-Williamson index (M = 0.146–0.279) and values of the heterozygosity excess revealed a reduction of genetic variation in all studied broodstocks because of the founder or bottleneck effect. The analysis of genetic differentiation (Fst) and Nei’s genetic distance between pairs of broodstocks revealed that Polish and Ukrainian broodstocks of European huchen were characterized by the closest genetic distance. In contrast, the highest genetic divergence parameters (Fst and Nei’s distance) were observed among German, Slovak, and Ukrainian broodstocks.

The complete mitochondrial genome of the endangered Hucho hucho (Salmonidae: Huchen)

Hucho hucho, one of the most endangered members of the family Salmonidae, is endemic to the Danube basin. In this study, the complete mitochondrial genome of H. hucho was sequenced and characterized. The genome is 16,751 bp in length and contains 13 protein-coding genes, 2 ribosomal RNAs, 22 transfer RNAs and a noncoding control region. The gene arrangement and nucleotide composition of the mitochondrial genome are similar to those of H. bleekeri. An 82 bp tandem repeat was identified in the control region, which is variable in length and copy number of repeat between and within species. The complete mitochondrial DNA of H. hucho should be useful to study population genetics, biogeography and adaptive evolution of this lineage.

Evolution of lifespan

Present-day evolutionary theory, modern synthesis and evo-devo, appear to explain evolution. There remain however several points of contention. These include: biological time, direction, macroevolution verses microevolution, ageing and the extent of internal as opposed to external mediation. A new theoretical model for the control of biological time in vertebrates/bilaterians is introduced. Rather than biological time being controlled solely by a molecular cascade domino effect, it is suggested there is also an intracellular oscillatory clock. This clock (life's timekeeper) is synchronised across all cells in an organism and runs at a constant frequency throughout life. Slower frequencies extend lifespan, increase body/brain size and advance behaviour. They also create a time void which could aid additional evolutionary change. Faster frequencies shorten lifespan, reduce body/brain size and diminish behaviour. They are therefore less likely to mediate evolution in vertebrates/mammals. It is concluded that in vertebrates, especially mammals, there is a direction in evolution towards longer lifespan/advanced behaviour. Lifespan extension could equate with macroevolution and subsequent modifications with microevolution. As life's timekeeper controls the rate of ageing it constitutes a new genetic theory of ageing. Finally, as lifespan extension is internally mediated, this suggests a major role for internal mediation in evolution.

Relaxed selection causes microevolution of seawater osmoregulation and gene expression in landlocked Alewives

Ecological transitions from marine to freshwater environments have been important in the creation of diversity among fishes. Evolutionary changes associated with these transitions likely involve modifications of osmoregulatory function. In particular, relaxed selection on hypo-osmoregulation should strongly affect animals that transition into novel freshwater environments. We used populations of the Alewife (Alosa pseudoharengus) to study evolutionary shifts in hypo-osmoregulatory capacity and ion regulation associated with freshwater transitions. Alewives are ancestrally anadromous, but multiple populations in Connecticut have been independently restricted to freshwater lakes; these landlocked populations complete their entire life cycle in freshwater. Juvenile landlocked and anadromous Alewives were exposed to three salinities (1, 20 and 30 ppt) in small enclosures within the lake. We detected strong differentiation between life history forms: landlocked Alewives exhibited reduced seawater tolerance and hypo-osmoregulatory performance compared to anadromous Alewives. Furthermore, gill Na(+)/K(+)-ATPase activity and transcription of genes for seawater osmoregulation (NKCC-Na(+)/K(+)/2Cl(-) cotransporter and CFTR-cystic fibrosis transmembrane conductance regulator) exhibited reduced responsiveness to seawater challenge. Our study demonstrates that adaptations of marine-derived species to completely freshwater life cycles involve partial loss of seawater osmoregulatory performance mediated through changes to ion regulation in the gill.

Comparative genome mapping between Chinook salmon (Oncorhynchus tshawytscha) and rainbow trout (O. mykiss) based on homologous microsatellite loci

Comparative genome mapping can rapidly facilitate the transfer of DNA sequence information from a well-characterized species to one that is less described. Chromosome arm numbers are conserved between members of the teleost family Salmonidae, order Salmoniformes, permitting rapid alignment of large syntenic blocks of DNA between members of the group. However, extensive Robertsonian rearrangements after an ancestral whole-genome duplication event has resulted in different chromosome numbers across Salmonid taxa. In anticipation of the rapid application of genomic data across members of the Pacific salmon genus Oncorhynchus, we mapped the genome of Chinook salmon (O. tshawytscha) by using 361 microsatellite loci and compared linkage groups to those already derived for a well-characterized species rainbow trout (O. mykiss). The Chinook salmon female map length was 1526 cM, the male map 733 cM, and the consensus map between the two sexes was 2206 cM. The average female to male recombination ratio was 5.43 (range 1-42.8 across all pairwise marker comparisons). We detected 34 linkage groups that corresponded with all chromosome arms mapped with homologous loci in rainbow trout and inferred that 16 represented metacentric chromosomes and 18 represented acrocentric chromosomes. Up to 13 chromosomes were conserved between the two species, suggesting that their structure precedes the divergence between Chinook salmon and rainbow trout. However, marker order differed in one of these linkage groups. The remaining linkage group structures reflected independent Robertsonian chromosomal arrangements, possibly after divergence. The putative linkage group homologies presented here are expected to facilitate future DNA sequencing efforts in Chinook salmon.

Genome duplication and multiple evolutionary origins of complex migratory behavior in Salmonidae

Multiple rounds of whole genome duplication have repeatedly marked the evolution of vertebrates, and correlate strongly with morphological innovation. However, less is known about the behavioral, physiological and ecological consequences of genome duplication, and whether these events coincide with major transitions in vertebrate complexity. The complex behavior of anadromy - where adult fishes migrate up rivers from the sea to their natal site to spawn - is well known in salmonid fishes. Some hypotheses suggest that migratory behavior evolved as a consequence of an ancestral genome duplication event, which permitted salinity tolerance and osmoregulatory plasticity. Here we test whether anadromy evolved multiple times within salmonids, and whether genome duplication coincided with the evolution of anadromy. We present a method that uses ancestral character simulation data to plot the frequency of character transitions over a time calibrated phylogenetic tree to provide estimates of the absolute timing of character state transitions. Furthermore, we incorporate extinct and extant taxa to improve on previous estimates of divergence times. We present the first phylogenetic evidence indicating that anadromy evolved at least twice from freshwater salmonid ancestors. Results suggest that genome duplication did not coincide in time with changes in migratory behavior, but preceded a transition to anadromy by 55-50 million years. Our study represents the first attempt to estimate the absolute timing of a complex behavioral trait in relation to a genome duplication event.

Comparative study of genome divergence in salmonids with various rates of genetic isolation

The aim of the study is a comparative investigation of changes that certain genome parts undergo during speciation. The research was focused on divergence of coding and noncoding sequences in different groups of salmonid fishes of the Salmonidae (Salmo, Parasalmo, Oncorhynchus, and Salvelinus genera) and the Coregonidae families under different levels of reproductive isolation. Two basic approaches were used: (1) PCR-RAPD with a 20–22 nt primer design with subsequent cloning and sequencing of the products and (2) a modified endonuclease restriction analysis. The restriction fragments were shown with sequencing to represent satellite DNA. Effects of speciation are found in repetitive sequences. The revelation of expressed sequences in the majority of the employed anonymous loci allows for assuming the adaptive selection during allopatric speciation in isolated char forms.

Molecular cloning and genomic characterization of novel leptin-like genes in salmonids provide new insight into the evolution of the Leptin gene family

In the current study we describe the identification of novel leptin B homologous gene/s in the four salmonid species Atlantic salmon (Salmo salar), rainbow trout (Oncorhynchus mykiss), brown trout (Salmo trutta) and Arctic charr (Salvelinus alpinus). Homology modeling of Salmo salar (Ss) LepB1/B2 suggests that the protein satisfies parameters as long-chain four helical cytokine family and that the basic structural pattern of the protein follows that of human leptin (Zhang et al., 1997). Importantly, the docking studies suggested the SsLepB has binding affinity to the AA residues that identify the leptin binding and FNIII domains of the SsLep receptor (Rønnestad et al., 2010). Phylogenetic analyses support that LepB paralogs have most probably originated by 4R whole genome duplication (WGD) before speciation of the salmonid lineages. LepB1 and LepB2 genes are both present in the two closest relatives, the Atlantic salmon and the brown trout, while rainbow trout and charr have only preserved the long LepB1 variant in their genome. We have defined the sites of SsLepB mRNA expression at key life stages in Atlantic salmon and found that SsLepB1 and SsLepB2, although to different extent, were expressed in redundant and mostly complementary fashion in brain and gills throughout the lifecycle, suggesting that this pair of paralogs is likely undergoing early stages of subfunctionalization. Furthermore, we have quantified the expression profiles of SsLepB genes and of other two recently duplicated salmon leptins (SsLepA1, SsLepA2) during early development and show evidence that in fish, as in mammals and amphibians, leptin could play important roles in growth and development. This study provides an essential groundwork to further elucidate structural and functional evolution of this important hormone in salmonids as well as in other teleosts.

Life's timekeeper

Life's timekeeper is a 'free-running' intracellular oscillator synchronised across all cells. It runs throughout life splitting lifespan into equal length phases. During the maturational period it controls the overall rate of progression whereas in the post-maturational period it controls the overall rate of ageing. This includes the rate of senescence and hence time to death. As such life's timekeeper equates maturational and post-maturational time, hence explains the tight correlation between these time periods that has existed throughout mammalian evolution. Life's timekeeper is proposed to have played an important role in vertebrate evolution. A slower oscillatory frequency results in proportional life phase prolongation. This leads to increased body and brain size, together with extended lifespan. Higher brain centres, neocortex in mammals, are disproportionately enlarged. Hence behavioural capacity is increased. The extended post-maturational period ensures that there is enough time in order that the behavioural advantages can be fully manifest in the environment. A faster oscillatory frequency would result in proportional life phase reduction. This process however would lead to reduced behavioural capacity, and is hence unlikely to be positively selected. Therefore throughout evolution life's timekeeper has operated to extend lifespan. It has hence functioned to promote longevity as opposed to ageing.

Comparative analysis of IgM sub-variants in salmonid fish and identification of a residue in μ3 which is essential for MAb4C10 reactivity

In rainbow trout (Oncorhynchus mykiss) it has been shown that high affinity IgM antibodies have a higher degree of disulfide polymerization and a longer half life time. In the present study, distinct IgM sub-variants related to ancestral tetraploidy in salmonid fish were analyzed to reveal possible characteristic differences between these. A monoclonal antibody (MAb4C10) which distinguishes between IgM-A and IgM-B in Atlantic salmon (Salmo salar) and brown trout (Salmo trutta) was further characterized. It was shown that substitution of a proline located in the loop between the B and C beta strands of the third constant domain (μ3) of salmon μA eliminated MAb4C10 reactivity. Accordingly, the reverse substitution in salmon μB restored MAb4C10 reactivity. Molecular cloning of μ cDNA from arctic char (Salvelinus alpinus) revealed two sub-variants (μA-1 and μA-2), i.e. a similar situation as in Atlantic salmon and brown trout. However, arctic char IgM eluted in one peak by anion exchange chromatography, in contrast to salmon and brown trout IgM that are eluted in two peaks. The only characteristic residue of salmon and brown trout μB is an additional cysteine in the C-terminal part of μ4. Most likely, this cysteine is involved in inter-chain disulfide bonding and influences the elution profiles of IgM-A and IgM-B on anion exchange chromatography. Neither of the μ sub-variants in arctic char have the additional cysteine, and char IgM, as well as salmon and brown trout IgM-A, showed a lower degree of inter-chain disulfide bonding than IgM-B when subjected to denaturation and gel electrophoresis under non-reducing conditions. Hybrids of char/salmon expressed μA-1, μA-2, μA and μB, indicating that there are two paralogous Ig heavy chain gene complexes in the haploid genome of char, like in Atlantic salmon. A comparison of salmonid μ sequences is presented, including representatives of Salmoninae (trout, salmon and char), Thymallinae (grayling) and Coregoninae (whitefish).

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

Considerable research efforts have focused on elucidating the systematic relationships among salmonid fishes; an understanding of these patterns of relatedness will inform conservation- and fisheries-related issues, as well as provide a framework for investigating evolutionary mechanisms in the group. However, uncertainties persist in current Salmonidae phylogenies due to biological and methodological factors, and a comprehensive phylogeny including most representatives of the family could provide insight into the causes of these difficulties. Here we increase taxon sampling by including nearly all described salmonid species (n = 63) to present a time-calibrated and more complete portrait of Salmonidae using a combination of molecular markers and analytical techniques. This strategy improved resolution by increasing the signal-to-noise ratio and helped discriminate methodological and systematic errors from sources of difficulty associated with biological processes. Our results highlight novel aspects of salmonid evolution. First, we call into question the widely-accepted evolutionary relationships among sub-families and suggest that Thymallinae, rather than Coregoninae, is the sister group to the remainder of Salmonidae. Second, we find that some groups in Salmonidae are older than previously thought and that the mitochondrial rate of molecular divergence varies markedly among genes and clades. We estimate the age of the family to be 59.1 MY (CI: 63.2-58.1 MY) old, which likely corresponds to the timing of whole genome duplication in salmonids. The average, albeit highly variable, mitochondrial rate of molecular divergence was estimated as ~0.31%/MY (CI: 0.27-0.36%/MY). Finally, we suggest that some species require taxonomic revision, including two monotypic genera, Stenodus and Salvethymus. In addition, we resolve some relationships that have been notoriously difficult to discern and present a clearer picture of the evolution of the group. Our findings represent an important contribution to the systematics of Salmonidae, and provide a useful tool for addressing questions related to fundamental and applied evolutionary issues.

The fitness consequences of environmental sex reversal in fish: a quantitative review

Environmental sex reversal (ESR) occurs when environmental factors overpower genetic sex-determining factors. The phenomenon of ESR is observed widely in teleost species, where it can be induced by exposing developing fish to endocrine disrupting chemicals (EDCs). EDC-induced ESR has been exploited by the aquaculture industry, while ecological and evolutionary models are also beginning to elucidate the potential roles that sex-reversed individuals play in influencing population dynamics. However, how EDC exposure affects individual fitness remains relatively unknown. To date, many experimental studies have induced sex reversal in fish and measured fitness-as indicated by related traits such as size, survival and gonadal somatic index (GSI), but the reported results vary. Here, we meta-analytically combine the results of 78 studies of induced ESR to gain insight into the fitness of sex-reversed individuals. Overall, our results suggest that the fitness of fish exposed to EDCs is reduced at the time of exposure, with exposed individuals having a smaller size and likely a smaller GSI. Given a period of non-exposure, fish treated with EDCs can regain a size equal to those not exposed, although GSI remains compromised. Interestingly, survival does not appear to be affected by EDC treatment. The published reports that comprise our dataset are, however, based on captive fish and the general small size resulting from exposure is likely to lead to reduced survival in the wild. Additionally, reduced fitness-related parameters are likely to be due to exposure to EDCs rather than ESR itself. We suggest that theoretical models of ESR should account for the fitness-related effects that we report. Whilst we are able to shed light on the physical fitness of EDC-exposed fish, the behaviour of such individuals remains largely untested and should be the focus of future experimental manipulation.

The complete mitochondrial genome of the taimen, Hucho taimen, and its unusual features in the control region

The whole mitochondrial genome of Hucho taimen was firstly sequenced and characterized. The genome is 16,833 bp in length and contains 13 protein-coding genes, 2 ribosomal RNAs, 22 transfer RNAs, and a noncoding control region. Twelve protein-coding genes on the heavy strand showed that the content of A+T was higher than that of G+C, whereas the nd6 protein-coding gene on the light strand displayed an opposite pattern. We described the secondary structure of the origin of light strand (oriL) replication and found that the conserved 5'-GCCGG-3' sequence motif is variable in H. taimen and some other salmonids. We conclude that the control region is variable in length and represents the high A+T content, compared with other mitochondrial control regions available in Salmonidae and other non-salmonids. Additionally, another interesting feature of H. taimen mitogenome is that a T-type mononucleotide microsatellite and an 82 bp tandem repeat were identified in the control region.

The influence of operational sex ratio on the intensity of competition for mates

The evolution and maintenance of secondary sexual characteristics and behavior are heavily influenced by the variance in mating success among individuals in a population. The operational sex ratio (OSR) is often used as a predictor of the intensity of competition for mates, as it describes the relative number of males and females who are ready to mate. We investigate changes in aggression, courtship, mate guarding, and sperm release as a function of changes in the OSR using meta-analytic techniques. As the OSR becomes increasingly biased, aggression increases as competitors attempt to defend mates, but this aggression begins to decrease at an OSR of 1.99, presumably due to the increased costs of competition as rivals become more numerous. Sperm release follows a similar but not significant trend. By contrast, courtship rate decreases as the OSR becomes increasingly biased, whereas mate guarding and copulation duration increase. Overall, predictable behavioral changes occur in response to OSR, although the nature of the change is dependent on the type of mating behavior. These results suggest considerable flexibility of mating system structure within species, which can be predicted by OSR and likely results in variation in the strength of sexual selection.

Detecting signatures of positive selection in partial sequences generated on a large scale: pitfalls, procedures and resources

Studying the actions of selection provides insight into adaptation, population divergence and gene function. Next-generation sequencing produces large amounts of partial sequences, potentially facilitating efforts to detect signatures of selection based on comparisons between synonymous (d(S)) and nonsynonymous (d(N)) substitutions, and single nucleotide polymorphism assays placed in selected genes would improve the ability to study adaptation in population surveys. However, sequences generated by these technologies are typically short. In nonmodel organisms that are a focus of evolutionary studies, the lack of a reference genome that facilitates the assembly of short sequences has limited surveys of positive selection in large numbers of genes. Here, we describe a series of steps to facilitate these surveys. We provide PERL scripts to assist data analysis, and describe the use of commonly available programs. We demonstrate these approaches in six salmon species, which have partially duplicated genomes. We recommend using multiway blast to optimize the number of alignments between partial coding sequences. Reading frames should be manually detected after alignment with sequences in Genbank using the BLASTX program. We encourage the use of a phylogenetic approach to separate orthologs from paralogs in duplicated genomes. Simple simulations on a gene known to have undergone selection in salmon species, transferrin, showed that the ability to detect selection in short sequences (<600 bp) depended on the proportion of codons under selection (1-2%) within that sequence. This relationship was less relevant in longer sequences. In this exploratory study, we detected 11 genes showing evidence of positive selection.