The black piranha (Serrasalmus rhombeus), a widely spread species in the rivers of the Amazon basin, plays a vital role as both key predator and important prey. Despite its essential contribution to ecosystem stability, there is a lack of information regarding its genetic diversity and population dynamics in the central Amazon region. As the Amazon continues to undergo environmental changes in the context of growing anthropogenic threats, such knowledge is fundamental for assist in the conservation of this species. This study is the first to analyze the genetic diversity and population structure of S. rhombeus in the central Amazon region using high-resolution genomic data. We employed a Genotyping-by-Sequencing approach with 248 samples across 14 study sites from various tributaries, encompassing diverse water types (black, white and clear water) and characterized by 34 physiochemical parameters. The data reveals low diversity accompanied by pronounced signs of inbreeding in half of the sites and robust genetic differentiation and variation among sites and within-sites. Surprisingly, we also found evidence of higher dispersal capacity than previously recognized. Our analysis exposed a complex and high population structure with genetic groups exclusive to some sites. Gene flow was low and some groups presented ambiguous genealogical divergence index (gdi) signals, suggesting the occurrence of potential cryptic species. Moreover, our results suggest that the the population structure of black piranha appears more influenced by historical events than contemporary factors. These results underscore the need to give greater attention to this keystone species, for which no regulatory framework or conservation strategies is presently in effect.

Single Nucleotide Polymorphism (SNP) panels are powerful tools for assessing genetic population structure and dispersal of fish that can enhance management practices for commercial, recreational, and subsistence mixed-stock fisheries. Arctic Char(Salvelinus alpinus), Brook Trout (Salvelinus fontinalis), and Lake Whitefish(Coregonus clupeaformis) are amongst the top harvested and consumed fish species in northern Indigenous communities in Canada, contributing significantly to food security, culture, and economy. Genomic resources like SNP panels, however, have not been widely accessible to Indigenous fisheries managers. We developed Genotyping-in-Thousands by sequencing (GT-seq) panels for population assignment and mixed-stock analyses for three salmonids to support fisheries co-management in northern Canada. Using low-coverage Whole-Genome Sequencing data from 418 individuals across source populations in Cambridge Bay (Nunavut), Great Slave Lake (Northwest Territories), James Bay (Québec), and Mistassini Lake (Québec), we developed a bioinformatic SNP filtering workflow to select informative SNP markers from genotype likelihoods. These markers were used to design GT-seq panels enabling high-throughput genotyping. Three GT-seq panels yielded an average of 413 autosomal loci and were validated with an average assignment accuracy of 98.5%. These GT-seq panels emerge as powerful tools for assessing population structure and quantifying the relative contributions of populations/stocks in mixed stock fisheries across multiple regions. Interweaving these novel genomic-derived tools with Traditional Ecological Knowledge will ensure the sustainable harvest of three culturally important salmonids in northern Indigenous communities, contributing to food security programs and the economy in northern Canada.

Amazonia is characterized by very heterogeneous riverscapes dominated by two drastically divergent water types: black (ion-poor, dissolved organic carbonate rich and acidic) and white (nutrient rich and turbid) waters. Recent phylogeographic and genomic studies have associated the ecotone formed by these environments to ecologically driven speciation in fishes. With the objective of better understanding the evolutionary forces behind the Amazonian Teleostean diversification, we sampled 240 Mesonauta festivus from 12 sites on a wide area of the Amazonian basin. These sites included three confluences of black and white water environments to seek for repeated evidences of ecological speciation at these ecotones. Results obtained through our genetic assessment based on 41,268 SNPs contrast with previous findings and supports a low structuring power of water types. Conversely, we detected a strong pattern of isolation by unidirectional downstream water current and evidence of past events of vicariance potentially linked to the Amazon River formation. Using a combination of population genetic, phylogeographic analysis and environmental association models, we decomposed the spatial variance from the environmental genetic variance specifically to assess which evolutive forces have shaped inter-population differences in M. festivus’ genome. Our sampling design comprising four major Amazonian rivers and three confluences of black and white water rivers supports the possibility that past studies potentially confounded ecological speciation with a site effect unrepresentative of the full Amazonian watershed. While ecological speciation admittedly played a role in Amazonian fish species diversification, we argue that neutral evolutionary processes explain most of the divergence between M. festivus populations.

Associations between host genotype and the microbiome of holobionts have been shown in a variety of animal clades, but studies on teleosts mostly show weak associations. Our study aimed to explore these relationships in four sympatric Serrasalmidae (i.e. piranha) teleosts from an Amazonian lake, using datasets from the hosts genomes (SNPs from GBS), skin and gut microbiomes (16S rRNA metataxonomics), and diets (COI metabarcoding) from the same fish individuals. Firstly, we investigated whether there were significant covariations of microbiome and fish genotypes at the inter and intraspecific scales. We also assessed the extent of co-variation between Serrasalmidae diet and microbiome, to isolate genotypic differences from dietary effects on community structure. We observed a significant covariation of skin microbiomes and host genotypes at interspecific (R2=24.4%) and intraspecific (R2=6.2%) scales, whereas gut microbiomes correlated poorly with host genotypes. Serrasalmidae diet composition was significantly correlated to fish genotype only at the interspecific scale (R2=5.4%), but did not covary with gut microbiome composition (Mantel R=-0.04; only 6 microbiome taxa involved). Secondly, we tested whether microbial taxa represent reliable host traits to complement host genotypic variations in these species. By using an NMDS ordination-based approach, we observed that subsets of the skin and gut microbiomes selected by a machine-learning Random Forest algorithm can complement host genotypic variations by increasing significantly the average interspecific differentiation. The complementarity of genome and microbiome variations suggests that combining both markers could potentially benefit our understanding of the evolution of Serrasalmidae in future studies.

Here we present an annotated, chromosome-anchored, genome assembly for Lake Trout (Salvelinus namaycush) – a highly diverse salmonid species of notable conservation concern and an excellent model for research on adaptation and speciation. We leveraged Pacific Biosciences long-read sequencing, paired-end Illumina sequencing, proximity ligation (Hi-C), and a previously published linkage map to produce a highly contiguous assembly composed of 7,378 contigs (contig N50 = 1.8 mb) assigned to 4,120 scaffolds (scaffold N50 = 44.975 mb). 84.7% of the genome was assigned to 42 chromosome-sized scaffolds and 93.2% of Benchmarking Universal Single Copy Orthologs were recovered, putting this assembly on par with the best currently available salmonid genomes. Estimates of genome size based on k-mer frequency analysis were highly similar to the total size of the finished genome, suggesting that the entirety of the genome was recovered. A mitome assembly was also produced. Self-vs-self synteny analysis allowed us to identify homeologs resulting from the Salmonid specific autotetraploid event (Ss4R) and alignment with three other salmonid species allowed us to identify homologous chromosomes in other species. We also generated multiple resources useful for future genomic research on Lake Trout including a repeat library and a sex averaged recombination map. A novel RNA sequencing dataset was also used to produce a publicly available set of gene annotations using the National Center for Biotechnology Information Eukaryotic Genome Annotation Pipeline. Potential applications of these resources to population genetics and the conservation of native populations are discussed.