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

Origins and functional diversification of salinity-responsive Na(+) , K(+) ATPase α1 paralogs in salmonids

The Salmoniform whole-genome duplication is hypothesized to have facilitated the evolution of anadromy, but little is known about the contribution of paralogs from this event to the physiological performance traits required for anadromy, such as salinity tolerance. Here, we determined when two candidate, salinity-responsive paralogs of the Na(+) , K(+) ATPase α subunit (α1a and α1b) evolved and studied their evolutionary trajectories and tissue-specific expression patterns. We found that these paralogs arose during a small-scale duplication event prior to the Salmoniform, but after the teleost, whole-genome duplication. The 'freshwater paralog' (α1a) is primarily expressed in the gills of Salmoniformes and an unduplicated freshwater sister species (Esox lucius) and experienced positive selection in the freshwater ancestor of Salmoniformes and Esociformes. Contrary to our predictions, the 'saltwater paralog' (α1b), which is more widely expressed than α1a, did not experience positive selection during the evolution of anadromy in the Coregoninae and Salmonine. To determine whether parallel mutations in Na(+) , K(+) ATPase α1 may contribute to salinity tolerance in other fishes, we studied independently evolved salinity-responsive Na(+) , K(+) ATPase α1 paralogs in Anabas testudineus and Oreochromis mossambicus. We found that a quarter of the mutations occurring between salmonid α1a and α1b in functionally important sites also evolved in parallel in at least one of these species. Together, these data argue that paralogs contributing to salinity tolerance evolved prior to the Salmoniform whole-genome duplication and that strong selection and/or functional constraints have led to parallel evolution in salinity-responsive Na(+) , K(+) ATPase α1 paralogs in fishes.

Comparative mapping between Coho Salmon (Oncorhynchus kisutch) and three other salmonids suggests a role for chromosomal rearrangements in the retention of duplicated regions following a whole genome duplication event

Whole genome duplication has been implicated in evolutionary innovation and rapid diversification. In salmonid fishes, however, whole genome duplication significantly pre-dates major transitions across the family, and re-diploidization has been a gradual process between genomes that have remained essentially collinear. Nevertheless, pairs of duplicated chromosome arms have diverged at different rates from each other, suggesting that the retention of duplicated regions through occasional pairing between homeologous chromosomes may have played an evolutionary role across species pairs. Extensive chromosomal arm rearrangements have been a key mechanism involved in re-dipliodization of the salmonid genome; therefore, we investigated their influence on degree of differentiation between homeologs across salmon species. We derived a linkage map for coho salmon and performed comparative mapping across syntenic arms within the genus Oncorhynchus, and with the genus Salmo, to determine the phylogenetic relationship between chromosome arrangements and the retention of undifferentiated duplicated regions. A 6596.7 cM female coho salmon map, comprising 30 linkage groups with 7415 and 1266 nonduplicated and duplicated loci, respectively, revealed uneven distribution of duplicated loci along and between chromosome arms. These duplicated regions were conserved across syntenic arms across Oncorhynchus species and were identified in metacentric chromosomes likely formed ancestrally to the divergence of Oncorhynchus from Salmo. These findings support previous studies in which observed pairings involved at least one metacentric chromosome. Re-diploidization in salmon may have been prevented or retarded by the formation of metacentric chromosomes after the whole genome duplication event and may explain lineage-specific innovations in salmon species if functional genes are found in these regions.

Development of a non-lethal method for evaluating transcriptomic endpoints in Arctic grayling (Thymallus arcticus)

With increases in active mining and continued discharge associated with former mine operations, evaluating the health of watersheds in the Canadian Yukon Territory is warranted. Current environmental assessment approaches often employ guidelines established using sentinel species not relevant to Arctic monitoring programs. The present study focused on the successful development of a quantitative real-time polymerase chain reaction (qPCR) assay directed towards the indigenous Arctic grayling (Thymallus arcticus) and examines the feasibility of using non-lethal sampling from the caudal fin as a means for evaluation of mRNA abundance profiles reflective of environmental conditions. In a proof of concept study performed blind, qPCR results from animals in an area with elevated water concentrations of cadmium (Cd) and zinc (Zn) and higher body burdens of Cd, Zn, and lead (Pb) were compared to a reference location in the Yukon Territory. Lower condition factor and a higher abundance of hepatic and caudal fin gene transcripts encoding the metallothionein isoforms (mta/mtb), in addition to elevated heat shock protein 70 (hsp70) and catalase (cat) mRNAs in liver, were observed in fish from the test site. The strong positive correlation between metal body burden and caudal fin mta/mtb mRNA abundance demonstrates a high potential for use of the Arctic grayling assay in non-lethal environmental monitoring programs.

Complete mitochondrial genome of Sakhalin taimen Parahucho perryi (Salmoniformes, Salmonidae) without two frame-disrupting indels in the ND4 gene

The complete mitochondrial genome was sequenced in two individuals of Sakhalin taimen Parahucho perryi. The genome sizes were 16,652 bp and 16,653 bp in the two isolates. Gene arrangement, base composition, and size of the two sequenced haplotypes are very similar to the P. perryi genome previously published (KC897653), but do not contain two frame-disrupting indels in the ND4 gene.

The rainbow trout genome provides novel insights into evolution after whole-genome duplication in vertebrates

Vertebrate evolution has been shaped by several rounds of whole-genome duplications (WGDs) that are often suggested to be associated with adaptive radiations and evolutionary innovations. Due to an additional round of WGD, the rainbow trout genome offers a unique opportunity to investigate the early evolutionary fate of a duplicated vertebrate genome. Here we show that after 100 million years of evolution the two ancestral subgenomes have remained extremely collinear, despite the loss of half of the duplicated protein-coding genes, mostly through pseudogenization. In striking contrast is the fate of miRNA genes that have almost all been retained as duplicated copies. The slow and stepwise rediploidization process characterized here challenges the current hypothesis that WGD is followed by massive and rapid genomic reorganizations and gene deletions.

Although whole-genome duplications (WGDs) are rare events, they have an important role in shaping vertebrate evolution. Here, the authors sequence the rainbow trout genome and show that rediploidization after WGD occurs in a slow and stepwise manner.

A dense linkage map for Chinook salmon (Oncorhynchus tshawytscha) reveals variable chromosomal divergence after an ancestral whole genome duplication event

Comparisons between the genomes of salmon species reveal that they underwent extensive chromosomal rearrangements following whole genome duplication that occurred in their lineage 58−63 million years ago. Extant salmonids are diploid, but occasional pairing between homeologous chromosomes exists in males. The consequences of re-diploidization can be characterized by mapping the position of duplicated loci in such species. Linkage maps are also a valuable tool for genome-wide applications such as genome-wide association studies, quantitative trait loci mapping or genome scans. Here, we investigated chromosomal evolution in Chinook salmon (Oncorhynchus tshawytscha) after genome duplication by mapping 7146 restriction-site associated DNA loci in gynogenetic haploid, gynogenetic diploid, and diploid crosses. In the process, we developed a reference database of restriction-site associated DNA loci for Chinook salmon comprising 48528 non-duplicated loci and 6409 known duplicated loci, which will facilitate locus identification and data sharing. We created a very dense linkage map anchored to all 34 chromosomes for the species, and all arms were identified through centromere mapping. The map positions of 799 duplicated loci revealed that homeologous pairs have diverged at different rates following whole genome duplication, and that degree of differentiation along arms was variable. Many of the homeologous pairs with high numbers of duplicated markers appear conserved with other salmon species, suggesting that retention of conserved homeologous pairing in some arms preceded species divergence. As chromosome arms are highly conserved across species, the major resources developed for Chinook salmon in this study are also relevant for other related species.

Linkage maps of the Atlantic salmon (Salmo salar) genome derived from RAD sequencing

BACKGROUND

Genetic linkage maps are useful tools for mapping quantitative trait loci (QTL) influencing variation in traits of interest in a population. Genotyping-by-sequencing approaches such as Restriction-site Associated DNA sequencing (RAD-Seq) now enable the rapid discovery and genotyping of genome-wide SNP markers suitable for the development of dense SNP linkage maps, including in non-model organisms such as Atlantic salmon (Salmo salar). This paper describes the development and characterisation of a high density SNP linkage map based on SbfI RAD-Seq SNP markers from two Atlantic salmon reference families.

RESULTS

Approximately 6,000 SNPs were assigned to 29 linkage groups, utilising markers from known genomic locations as anchors. Linkage maps were then constructed for the four mapping parents separately. Overall map lengths were comparable between male and female parents, but the distribution of the SNPs showed sex-specific patterns with a greater degree of clustering of sire-segregating SNPs to single chromosome regions. The maps were integrated with the Atlantic salmon draft reference genome contigs, allowing the unique assignment of ~4,000 contigs to a linkage group. 112 genome contigs mapped to two or more linkage groups, highlighting regions of putative homeology within the salmon genome. A comparative genomics analysis with the stickleback reference genome identified putative genes closely linked to approximately half of the ordered SNPs and demonstrated blocks of orthology between the Atlantic salmon and stickleback genomes. A subset of 47 RAD-Seq SNPs were successfully validated using a high-throughput genotyping assay, with a correspondence of 97% between the two assays.

CONCLUSIONS

This Atlantic salmon RAD-Seq linkage map is a resource for salmonid genomics research as genotyping-by-sequencing becomes increasingly common. This is aided by the integration of the SbfI RAD-Seq SNPs with existing reference maps and the draft reference genome, as well as the identification of putative genes proximal to the SNPs. Differences in the distribution of recombination events between the sexes is evident, and regions of homeology have been identified which are reflective of the recent salmonid whole genome duplication.

A well-constrained estimate for the timing of the salmonid whole genome duplication reveals major decoupling from species diversification

Whole genome duplication (WGD) is often considered to be mechanistically associated with species diversification. Such ideas have been anecdotally attached to a WGD at the stem of the salmonid fish family, but remain untested. Here, we characterized an extensive set of gene paralogues retained from the salmonid WGD, in species covering the major lineages (subfamilies Salmoninae, Thymallinae and Coregoninae). By combining the data in calibrated relaxed molecular clock analyses, we provide the first well-constrained and direct estimate for the timing of the salmonid WGD. Our results suggest that the event occurred no later in time than 88 Ma and that 40-50 Myr passed subsequently until the subfamilies diverged. We also recovered a Thymallinae-Coregoninae sister relationship with maximal support. Comparative phylogenetic tests demonstrated that salmonid diversification patterns are closely allied in time with the continuous climatic cooling that followed the Eocene-Oligocene transition, with the highest diversification rates coinciding with recent ice ages. Further tests revealed considerably higher speciation rates in lineages that evolved anadromy--the physiological capacity to migrate between fresh and seawater--than in sister groups that retained the ancestral state of freshwater residency. Anadromy, which probably evolved in response to climatic cooling, is an established catalyst of genetic isolation, particularly during environmental perturbations (for example, glaciation cycles). We thus conclude that climate-linked ecophysiological factors, rather than WGD, were the primary drivers of salmonid diversification.

The covert world of fish biofluorescence: a phylogenetically widespread and phenotypically variable phenomenon

The discovery of fluorescent proteins has revolutionized experimental biology. Whereas the majority of fluorescent proteins have been identified from cnidarians, recently several fluorescent proteins have been isolated across the animal tree of life. Here we show that biofluorescence is not only phylogenetically widespread, but is also phenotypically variable across both cartilaginous and bony fishes, highlighting its evolutionary history and the possibility for discovery of numerous novel fluorescent proteins. Fish biofluorescence is especially common and morphologically variable in cryptically patterned coral-reef lineages. We identified 16 orders, 50 families, 105 genera, and more than 180 species of biofluorescent fishes. We have also reconstructed our current understanding of the phylogenetic distribution of biofluorescence for ray-finned fishes. The presence of yellow long-pass intraocular filters in many biofluorescent fish lineages and the substantive color vision capabilities of coral-reef fishes suggest that they are capable of detecting fluoresced light. We present species-specific emission patterns among closely related species, indicating that biofluorescence potentially functions in intraspecific communication and evidence that fluorescence can be used for camouflage. This research provides insight into the distribution, evolution, and phenotypic variability of biofluorescence in marine lineages and examines the role this variation may play.

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.

Targeted sequence capture and resequencing implies a predominant role of regulatory regions in the divergence of a sympatric lake whitefish species pair (Coregonus clupeaformis)

Latest technological developments in evolutionary biology bring new challenges in documenting the intricate genetic architecture of species in the process of divergence. Sympatric populations of lake whitefish represent one of the key systems to investigate this issue. Despite the value of random genotype-by-sequencing methods and decreasing cost of sequencing technologies, it remains challenging to investigate variation in coding regions, especially in the case of recently duplicated genomes as in salmonids, as this greatly complicates whole genome resequencing. We thus designed a sequence capture array targeting 2773 annotated genes to document the nature and the extent of genomic divergence between sympatric dwarf and normal whitefish. Among the 2728 genes successfully captured, a total of 2182 coding and 10,415 noncoding putative single-nucleotide polymorphisms (SNPs) were identified after applying a first set of basic filters. A genome scan with a quality-refined selection of 2203 SNPs identified 267 outlier SNPs in 210 candidate genes located in genomic regions potentially involved in whitefish divergence and reproductive isolation. We found highly heterogeneous FST estimates among SNP loci. There was an overall low level of coding polymorphism, with a predominance of noncoding mutations among outliers. The heterogeneous patterns of divergence among loci confirm the porous nature of genomes during speciation with gene flow. Considering that few protein-coding mutations were identified as highly divergent, our results, along with previous transcriptomic studies, imply that changes in regulatory regions most likely had a greater role in the process of whitefish population divergence than protein-coding mutations. This study is the first to demonstrate the efficiency of large-scale targeted resequencing for a nonmodel species with such a large and unsequenced genome.

Pike and salmon as sister taxa: detailed intraclade resolution and divergence time estimation of Esociformes + Salmoniformes based on whole mitochondrial genome sequences

The increasing number of taxa and loci in molecular phylogenetic studies of basal euteleosts has brought stability in a controversial area. A key emerging aspect to these studies is a sister Esociformes (pike) and Salmoniformes (salmon) relationship. We evaluate mitochondrial genome support for a sister Esociformes and Salmoniformes hypothesis by surveying many potential outgroups for these taxa, employing multiple phylogenetic approaches, and utilizing a thorough sampling scheme. Secondly, we conduct a simultaneous divergence time estimation and phylogenetic inference in a Bayesian framework with fossil calibrations focusing on relationships within Esociformes+Salmoniformes. Our dataset supports a sister relationship between Esociformes and Salmoniformes; however the nearest relatives of Esociformes+Salmoniformes are inconsistent among analyses. Within the order Esociformes, we advocate for a single family, Esocidae. Subfamily relationships within Salmonidae are poorly supported as Salmoninae sister to Thymallinae+Coregoninae.

Mitochondrial DNA variation and introgression in Siberian taimen Hucho taimen

Siberian taimen Hucho taimen is the largest representative of the family Salmonidae inhabiting rivers of northern Eurasia. The species is under intensive aquaculture activity. To monitor natural taimen populations we have sequenced a portion (8,141 bp) of the mitochondrial (mt) genome in 28 specimens of H. taimen from six localities in the Amur River basin. Nucleotide variability is low (π = 0.0010), but structured in two divergent haplotype groups. A comparison of the data with the GenBank H. taimen mt genome (HQ897271) reveals significant differences between them in spite of the fact that the fish specimens come from neighboring geographical areas. The distribution of divergence is non-uniform with two highly pronounced divergent regions centered on two genes, ND3 and ND6. To clarify the pattern of divergence we sequenced the corresponding portion of the mt genome of lenok Brachymystax tumensis and analyzed the GenBank complete mt genomes of related species. We have found that the first and second divergent regions are identical between the GenBank H. taimen and two lenok subspecies, B. lenok and B. lenok tsinlingensis, respectively. Consequently, both divergent regions represent introgressed mtDNA resulting from intergeneric hybridization between the two lenok subspecies and H. taimen. Introgression is, however, not detected in our specimens. This plus the precise identity of the introgressed fragments between the donor and the recipient GenBank sequence suggests that the introgression is local and very recent, probably due to artificial manipulations involving taimen – lenok intergeneric hybridization. Human-mediated hybridization may become a major threat to aquatic biodiversity. Consequently we suggest that due attention needs to be given to this threat by means of responsible breeding program management, so as to prevent a potential spread of hybrid fishes that could jeopardize the resilience of locally adapted gene pools of the native H. taimen populations.

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.

Complete mitochondrial genome of the endangered Sakhalin taimen Parahucho perryi (Salmoniformes, Salmonidae)

The complete mitochondrial genome sequence of Parahucho perryi has been obtained by the next generation sequencing, which contained 22 tRNA genes, 13 protein-coding genes, 2 rRNA genes and non-coding control region with the total length of 16,651 bp. The gene content, arrangement, codon usage and base composition of P. perryi mitogenome are identical to those observed in salmonids and most other teleost fishes. The sequence data could provide useful information for the studies on molecular systematics and conservation genetics.

The sexually dimorphic on the Y-chromosome gene (sdY) is a conserved male-specific Y-chromosome sequence in many salmonids

All salmonid species investigated to date have been characterized with a male heterogametic sex-determination system. However, as these species do not share any Y-chromosome conserved synteny, there remains a debate on whether they share a common master sex-determining gene. In this study, we investigated the extent of conservation and evolution of the rainbow trout (Oncorhynchus mykiss) master sex-determining gene, sdY (sexually dimorphic on the Y-chromosome), in 15 different species of salmonids. We found that the sdY sequence is highly conserved in all salmonids and that sdY is a male-specific Y-chromosome gene in the majority of these species. These findings demonstrate that most salmonids share a conserved sex-determining locus and also strongly suggest that sdY may be this conserved master sex-determining gene. However, in two whitefish species (subfamily Coregoninae), sdY was found both in males and females, suggesting that alternative sex-determination systems may have also evolved in this family. Based on the wide conservation of sdY as a male-specific Y-chromosome gene, efficient and easy molecular sexing techniques can now be developed that will be of great interest for studying these economically and environmentally important species.

Robust detection of rare species using environmental DNA: the importance of primer specificity

Environmental DNA (eDNA) is being rapidly adopted as a tool to detect rare animals. Quantitative PCR (qPCR) using probe-based chemistries may represent a particularly powerful tool because of the method’s sensitivity, specificity, and potential to quantify target DNA. However, there has been little work understanding the performance of these assays in the presence of closely related, sympatric taxa. If related species cause any cross-amplification or interference, false positives and negatives may be generated. These errors can be disastrous if false positives lead to overestimate the abundance of an endangered species or if false negatives prevent detection of an invasive species. In this study we test factors that influence the specificity and sensitivity of TaqMan MGB assays using co-occurring, closely related brook trout (Salvelinus fontinalis) and bull trout (S. confluentus) as a case study. We found qPCR to be substantially more sensitive than traditional PCR, with a high probability of detection at concentrations as low as 0.5 target copies/µl. We also found that number and placement of base pair mismatches between the Taqman MGB assay and non-target templates was important to target specificity, and that specificity was most influenced by base pair mismatches in the primers, rather than in the probe. We found that insufficient specificity can result in both false positive and false negative results, particularly in the presence of abundant related species. Our results highlight the utility of qPCR as a highly sensitive eDNA tool, and underscore the importance of careful assay design.

Use of multiple markers demonstrates a cryptic western refugium and postglacial colonisation routes of Atlantic salmon (Salmo salar L.) in Northwest Europe

Glacial and postglacial processes are known to be important determinants of contemporary population structuring for many species. In Europe, refugia in the Italian, Balkan and Iberian peninsulas are believed to be the main sources of species colonising northern Europe after the glacial retreat; however, there is increasing evidence of small, cryptic refugia existing north of these for many cold-tolerant species. This study examined the glacial history of Atlantic salmon in western Europe using two independent classes of molecular markers, microsatellites (nuclear) and mitochondrial DNA variation. Alongside the well-documented refuge in the Iberian Peninsula, evidence for a cryptic refuge in northwest France is also presented. Critically, methods utilised to estimate divergence times between the refugia indicated that salmon in these two regions had diverged a long time before the last glacial maximum; coalescence analysis (as implemented in the program IMa2) estimated divergence times at around 60 000 years before present. Through the examination of haplotype frequencies, previously glaciated areas of northwest Europe, that is, Britain and Ireland, appear to have been colonised from salmon expanding out of both refugia, with the southwest of England being the primary contact zone and exhibiting the highest genetic diversity.