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

Origin of the Laurentian Great Lakes fish fauna through upward adaptive radiation cascade prior to the Last Glacial Maximum

The evolutionary histories of adaptive radiations can be marked by dramatic demographic fluctuations. However, the demographic histories of ecologically-linked co-diversifying lineages remain understudied. The Laurentian Great Lakes provide a unique system of two such radiations that are dispersed across depth gradients with a predator-prey relationship. We show that the North American Coregonus species complex ("ciscoes") radiated rapidly prior to the Last Glacial Maximum (80-90 kya), a globally warm period, followed by rapid expansion in population size. Similar patterns of demographic expansion were observed in the predator species, Lake Charr (Salvelinus namaycush), following a brief time lag, which we hypothesize was driven by predator-prey dynamics. Diversification of prey into deep water created ecological opportunities for the predators, facilitating their demographic expansion, which is consistent with an upward adaptive radiation cascade. This study provides a new timeline and environmental context for the origin of the Laurentian Great Lakes fish fauna, and firmly establishes this system as drivers of ecological diversification and rapid speciation through cyclical glaciation.

Genomic structural variation in Barramundi Perch Lates calcarifer and potential roles in speciation and adaptation

Advancements in genome sequencing and assembly techniques have increased the documentation of structural variants in wild organisms. Of these variants, chromosomal inversions are especially prominent due to their large size and active recombination suppression between alternative homokaryotypes. This suppression enables the 2 forms of the inversion to be maintained and allows the preservation of locally adapted alleles. The Barramundi Perch (BP; Lates calcarifer) is a widespread species complex with 3 main genetic lineages located in the biogeographic regions of Australia and New Guinea (AUS + NG), Southeast Asia (SEA), and the Indian Subcontinent (IND). BP are typically considered to be a protandrous sequential hermaphrodite species that exhibits catadromy. Freshwater occupancy and intraspecific variation in life history (e.g. partially migratory populations) exist and provide opportunities for strongly divergent selection associated with, for example, salinity tolerance, swimming ability, and marine dispersal. Herein, we utilize genomic data generated from all 3 genetic lineages to identify and describe 3 polymorphic candidate chromosomal inversions. These candidate chromosomal inversions appear to be fixed for ancestral variants in the IND lineage and for inverted versions in the AUS + NG lineage and exhibit variation in all 3 inversions in the SEA lineage. BP have a diverse portfolio of life history options that includes migratory strategy as well as sexual system (i.e. hermaphroditism and gonochorism). We propose that the some of the life history variabilities observed in BP may be linked to inversions and, in doing so, we present genetic data that might be useful in enhancing aquaculture production and population management.

The development of the digestive system and the fate of the yolk syncytial layer in postembryogenesis of Stenodus leucichthys nelma (Teleostei)

Stenodus leucichthys nelma is an economically important species for cold-water aquaculture. Unlike other Coregoninae, S. leucichthys nelma is a piscivore. Here, we describe in detail the development of the digestive system and the yolk syncytial layer from hatching to early juvenile stage using histological and histochemical methods to determine their common and specific characteristics and to test the hypothesis that the digestive system of S. leucichthys nelma rapidly acquires adult features. The digestive tract differentiates at hatching and starts to function before the transition to mixed feeding. The mouth and anus are open, mucous cells and taste buds are present in the buccopharyngeal cavity and esophagus, pharyngeal teeth have erupted, the stomach primordium is seen, the intestinal epithelium with mucous cells is folded and the intestinal valve is observed; the epithelial cells of the postvalvular intestine contain supranuclear vacuoles. The liver blood vessels are filled with blood. The cells of exocrine pancreas are loaded with zymogen granules, and at least two islets of Langerhans are present. However, the larvae remain dependent on maternal yolk and lipids for a long time. The adult features of the digestive system develop gradually, the most significant changes take place approximately from 31 to 42 days posthatching. Then, the gastric glands and pyloric caeca buds appear, the U-shaped stomach with glandular and aglandular regions develops, the swim bladder inflates, the number of islets of Langerhans increases, the pancreas becomes scattered, and the yolk syncytial layer undergoes programmed death during the larval-to-juvenile transition. During postembryonic development, the mucous cells of the digestive system contain neutral mucosubstances.

Genome assembly, structural variants, and genetic differentiation between lake whitefish young species pairs (Coregonus sp.) with long and short reads

Nascent pairs of ecologically differentiated species offer an opportunity to get a better glimpse at the genetic architecture of speciation. Of particular interest is our recent ability to consider a wider range of genomic variants, not only single-nucleotide polymorphisms (SNPs), thanks to long-read sequencing technology. We can now identify structural variants (SVs) such as insertions, deletions and other rearrangements, allowing further insights into the genetic architecture of speciation and how different types of variants are involved in species differentiation. Here, we investigated genomic patterns of differentiation between sympatric species pairs (Dwarf and Normal) belonging to the lake whitefish (Coregonus clupeaformis) species complex. We assembled the first reference genomes for both C. clupeaformis sp. Normal and C. clupeaformis sp. Dwarf, annotated the transposable elements and analysed the genomes in the light of related coregonid species. Next, we used a combination of long- and short-read sequencing to characterize SVs and genotype them at the population scale using genome-graph approaches, showing that SVs cover five times more of the genome than SNPs. We then integrated both SNPs and SVs to investigate the genetic architecture of species differentiation in two different lakes and highlighted an excess of shared outliers of differentiation. In particular, a large fraction of SVs differentiating the two species correspond to insertions or deletions of transposable elements (TEs), suggesting that TE accumulation may represent a key component of genetic divergence between the Dwarf and Normal species. Together, our results suggest that SVs may play an important role in speciation and that, by combining second- and third-generation sequencing, we now have the ability to integrate SVs into speciation genomics.

Genomic and environmental influences on resilience in a cold-water fish near the edge of its range

Small, isolated populations present a challenge for conservation. The dueling effects of selection and drift in a limited pool of genetic diversity make the responses of small populations to environmental perturbations erratic and difficult to predict. This is particularly true at the edge of a species range, where populations often persist at the limits of their environmental tolerances. Populations of cisco, Coregonus artedi, in inland lakes have experienced numerous extirpations along the southern edge of their range in recent decades, which are thought to result from environmental degradation and loss of cold, well-oxygenated habitat as lakes warm. Yet, cisco extirpations do not show a clear latitudinal pattern, suggesting that local environmental factors and potentially local adaptation may influence resilience. Here, we used genomic tools to investigate the nature of this pattern of resilience. We used restriction site-associated DNA capture (Rapture) sequencing to survey genomic diversity and differentiation in southern inland lake cisco populations and compared the frequency of deleterious mutations that potentially influence fitness across lakes. We also examined haplotype diversity in a region of the major histocompatibility complex involved in stress and immune system response. We correlated these metrics to spatial and environmental factors including latitude, lake size, and measures of oxythermal habitat and found significant relationships between genetic metrics and broad and local factors. High levels of genetic differentiation among populations were punctuated by a phylogeographic break and residual patterns of isolation-by-distance. Although the prevalence of deleterious mutations and inbreeding coefficients was significantly correlated with latitude, neutral and non-neutral genetic diversity were most strongly correlated with lake surface area. Notably, differences among lakes in the availability of estimated oxythermal habitat left no clear population genomic signature. Our results shed light on the complex dynamics influencing these isolated populations and provide valuable information for their conservation.

Deciphering reticulate evolution of the largest group of polyploid vertebrates, the subfamily cyprininae (Teleostei: Cypriniformes)

Despite the rarity of polyploidy in animals, some groups with polyploid species exhibit complicated and interesting patterns of reticulate evolution. Here we focus on fishes in the subfamily Cyprininae, the largest polyploid group of vertebrates. The large number of polyploid taxa poses significant challenges for phylogenetic and evolutionary studies on this subfamily. In this study, we cloned and sequenced three single-copy nuclear loci to investigate the evolution of polyploidy in the Cyprininae. Topologies of nuclear gene trees were compared with a newly reconstructed mitochondrial tree. The data provided herein corroborate the hybrid origins of the tribes Torini, Cyprinini, Spinibarbini, Barbini, and also Probarbini. Based on results from this study and previous studies, we hypothesize that at least 13 independent polyploidization events have occurred during the evolution of the Cyprininae. We offer hypotheses on the origin of each polyploid group and show that a diploid group or the diploid ancestor of a polyploid group might have served as progenitor of one or two other polyploid groups. The evolutionary history of Cyprinine (since its first divergence) can be divided into three stages: the "Diploid stage" (69.2-43.4 Ma or million years ago), the "Tetraploidization stage" (43.4-18.9 Ma), and the "Hexaploidization stage" (18.9 Ma to present). The second stage is when all tetraploidization events happened, while the last stage is when all hexaploidization events and most genus- or species-specific polyploidization events occurred. The post-polyploidization dynamics of polyploid groups are complicated and warrant more genomic level studies. We showed that the subfamily Cyprininae may represent a more complicated polyploid system than most, if not all, other vertebrates and some plants, if one or more of the following factors are considered: number of polyploid species, number of different ploidy levels, and number and type of independent polyploidization events.

A candidate chromosome inversion in Arctic charr (Salvelinus alpinus) identified by population genetic analysis techniques

The "genomics era" has allowed questions to be asked about genome organization and genome architecture of non-model species at a rate not previously seen. Analyses of these genome-wide datasets have documented many examples of novel structural variants (SVs) such as chromosomal inversions, copy number variants, and chromosomal translocations, many of which have been linked to adaptation. The salmonids are a taxonomic group with abundant genome-wide datasets due to their importance in aquaculture and fisheries. However, the number of documented SVs in salmonids is surprisingly low and is most likely due to removing loci in high linkage disequilibrium when analyzing structure and gene flow. Here we re-analyze RAD-seq data from several populations of Arctic charr (Salvelinus alpinus) and document a novel ∼1.2 MB SV at the distal end of LG12. This variant contains 15 protein-coding genes connected to a wide-range of functions including cell adhesion and signal transduction. Interestingly, we studied the frequency of this polymorphism in four disjointed populations of charr-one each from Nunavut, Newfoundland, Eastern Russia, and Scotland-and found evidence of the variant only in Nunavut, Canada, suggesting the polymorphism is novel and recently evolved.

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.

A de novo chromosome-level genome assembly of Coregonus sp. "Balchen": One representative of the Swiss Alpine whitefish radiation

Salmonids are of particular interest to evolutionary biologists due to their incredible diversity of life-history strategies and the speed at which many salmonid species have diversified. In Switzerland alone, over 30 species of Alpine whitefish from the subfamily Coregoninae have evolved since the last glacial maximum, with species exhibiting a diverse range of morphological and behavioural phenotypes. This, combined with the whole genome duplication which occurred in the ancestor of all salmonids, makes the Alpine whitefish radiation a particularly interesting system in which to study the genetic basis of adaptation and speciation and the impacts of ploidy changes and subsequent rediploidization on genome evolution. Although well-curated genome assemblies exist for many species within Salmonidae, genomic resources for the subfamily Coregoninae are lacking. To assemble a whitefish reference genome, we carried out PacBio sequencing from one wild-caught Coregonus sp. "Balchen" from Lake Thun to ~90× coverage. PacBio reads were assembled independently using three different assemblers, falcon, canu and wtdbg2 and subsequently scaffolded with additional Hi-C data. All three assemblies were highly contiguous, had strong synteny to a previously published Coregonus linkage map, and when mapping additional short-read data to each of the assemblies, coverage was fairly even across most chromosome-scale scaffolds. Here, we present the first de novo genome assembly for the Salmonid subfamily Coregoninae. The final 2.2-Gb wtdbg2 assembly included 40 scaffolds, an N50 of 51.9 Mb and was 93.3% complete for BUSCOs. The assembly consisted of ~52% transposable elements and contained 44,525 genes.