Morphological and molecular characterization of Stomachicola muraenesocis Yamaguti, 1934 (Digenea: Hemiuridae) from the daggertooth pike conger Muraenesox cinereus (Forsskål).

Hemiurid digeneans conspecific with Stomachicola muraenesocis Yamaguti, 1934 (the type species of the genus Stomachicola Yamaguti, 1934) were collected from the stomach of the daggertooth pike conger Muraenesox cinereus (Forsskål) off the Persian Gulf of Iran. This study aimed to provide a detailed characterization of Stom. muraenesocis, including measurements, illustrations and scanning electron microscopy (s.e.m.) representations. Comparisons with the original and previous descriptions revealed morphological and metrical variations in several features (i.e. body size and shape, arrangement of reproductive organs, soma to ecsoma length ratio, position of genital opening, number of vitelline tubules and extension of uterine coils) between Stom. muraenesocis from different hosts and localities. This study presents the first molecular sequence data associated with the small (18S) and large (28S) subunit nuclear ribosomal RNA genes (rDNA) for Stom. muraenesocis. Phylogenetic analyses of the 18S dataset placed Stom. muraenesocis as sister lineage to a clade formed of a group of species of Lecithaster Lühe, 1901 (Lecithasteridae Odhner, 1905). In contrast, phylogenetic analyses based on the 28S consistently recovered a sister relationship between Stom. muraenesocis and representatives of the Hemiuridae Looss, 1899. Further comprehensive phylogenetically based classification in light of morphology and taxonomic history of the Hemiuridae and Lecithasteridae is required to infer phylogenetic affinities and historical biogeography of Stomachicola. A comprehensive list of previously reported species of Stomachicola together with their associated hosts, localities and morphometric data is provided.

Loss of gut microbial diversity in the cultured, agastric fish, Mexican pike silverside (Chirostoma estor: Atherinopsidae).

Teleost fish are the most diverse group of extant vertebrates and have varied digestive anatomical structures and strategies, suggesting they also possess an array of different host-microbiota interactions. Differences in fish gut microbiota have been shown to affect host development, the process of gut colonization, and the outcomes of gene-environment or immune system-microbiota interactions. There is generally a lack of studies on the digestive mechanisms and microbiota of agastric short-intestine fish however, meaning that we do not understand how changes in gut microbial diversity might influence the health of these types of fish. To help fill these gaps in knowledge, we decided to study the Mexican pike silverside (Chirostoma estor) which has a simplified alimentary canal (agastric, short-intestine, 0.7 gut relative length) to observe the diversity and metabolic potential of its intestinal microbiota. We characterized gut microbial populations using high-throughput sequencing of the V3 region in bacterial 16S rRNA genes while searching for population shifts resulting associated with fish development in different environments and cultivation methods. Microbiota samples were taken from the digesta, anterior and posterior intestine (the three different intestinal components) of fish that grew wild in a lake, that were cultivated in indoor tanks, or that were raised in outdoor ponds. Gut microbial diversity was significantly higher in wild fish than in cultivated fish, suggesting a loss of diversity when fish are raised in controlled environments. The most abundant phyla observed in these experiments were Firmicutes and Proteobacteria, particularly of the genera Mycoplasma, Staphylococcus, Spiroplasma, and Aeromonas. Of the 14,161 OTUs observed in this experiment, 133 were found in all groups, and 17 of these, belonging to Acinetobacter, Aeromonas, Pseudomonas, and Spiroplasma genera, were found in all samples suggesting the existence of a core C. estor microbiome. Functional metagenomic prediction of bacterial ecological functions using PICRUSt2 suggested that different intestinal components select for functionally distinct microbial populations with variation in pathways related to the metabolism of amino acids, vitamins, cofactors, and energy. Our results provide, for the first time, information on the bacterial populations present in an agastric, short-gut teleost with commercial potential and show that controlled cultivation of this fish reduces the diversity of its intestinal microbiota.

Drivers of neutral and adaptive differentiation in pike (Esox lucius) populations from contrasting environments.

Understanding how eco-evolutionary processes and environmental factors drive population differentiation and adaptation are key challenges in evolutionary biology of relevance for biodiversity protection. Differentiation requires at least partial reproductive separation, which may result from different modes of isolation such as geographic isolation (allopatry) or isolation by distance (IBD), resistance (IBR), and environment (IBE). Despite that multiple modes might jointly influence differentiation, studies that compare the relative contributions are scarce. Using RADseq, we analyse neutral and adaptive genetic diversity and structure in 11 pike (Esox lucius) populations from contrasting environments along a latitudinal gradient (54.9-63.6°N), to investigate the relative effects of IBD, IBE and IBR, and to assess whether the effects differ between neutral and adaptive variation, or across structural levels. Patterns of neutral and adaptive variation differed, probably reflecting that they have been differently affected by stochastic and deterministic processes. The importance of the different modes of isolation differed between neutral and adaptive diversity, yet were consistent across structural levels. Neutral variation was influenced by interactions among all three modes of isolation, with IBR (seascape features) playing a central role, wheares adaptive variation was mainly influenced by IBE (environmental conditions). Taken together, this and previous studies suggest that it is common that multiple modes of isolation interactively shape patterns of genetic variation, and that their relative contributions differ among systems. To enable identification of general patterns and understand how various factors influence the relative contributions, it is important that several modes are simultaneously investigated in additional populations, species and environmental settings.

The rise and fall of the ancient northern pike master sex-determining gene.

The understanding of the evolution of variable sex determination mechanisms across taxa requires comparative studies among closely related species. Following the fate of a known master sex-determining gene, we traced the evolution of sex determination in an entire teleost order (Esociformes). We discovered that the northern pike (Esox lucius) master sex-determining gene originated from a 65 to 90 million-year-old gene duplication event and that it remained sex linked on undifferentiated sex chromosomes for at least 56 million years in multiple species. We identified several independent species- or population-specific sex determination transitions, including a recent loss of a Y chromosome. These findings highlight the diversity of evolutionary fates of master sex-determining genes and the importance of population demographic history in sex determination studies. We hypothesize that occasional sex reversals and genetic bottlenecks provide a non-adaptive explanation for sex determination transitions.

Environmental DNA detection of endangered and invasive species in Kejimkujik National Park and Historic Site.

The use of environmental DNA (eDNA) allows the early detection of aquatic species at low densities (e.g., elusive and invasive species), which otherwise could be challenging to monitor using conventional techniques. Here, we assess the ability of eDNA sampling to detect the presence or absence of one species at risk (Blanding's turtle) and two invasive species (chain pickerel and smallmouth bass) in Kejimkujik National Park and National Historic Site, Nova Scotia, where the aquatic system is highly acidic and rich in organic compounds. Five replicates of 1 L water samples were taken per sampling site. Water filtration and eDNA extractions were performed on-site, while qPCR reactions were performed in the laboratory using species-specific assays. Samples were treated with an inhibition removal kit and analyzed pre- and post-inhibition removal. Despite the low pH and PCR inhibitors in water samples, our results showed positive eDNA detections in almost all expected positive sites (except in one site for Blanding's turtle). Detections of the target species were also observed at sites where their presence was previously unknown. Our study supports the advantage of eDNA to monitor species at low densities, revealing new distributions or recently invaded areas. We also demonstrate how eDNA can directly instruct management strategies in Kejimkujik.

Fate of MHCII in salmonids following 4WGD.

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

Environmental exposure of northern pike to a primary wastewater effluent: Impact on the lipidomic profile and lipid metabolism.

Lipids play important roles in growth, reproduction, locomotion, and migration of fish. Municipal effluents, which are complex mixtures of biological and chemical compounds including flame retardants, have been shown to alter lipid metabletabolism in environmentally and experimentally exposed fish. Down-regulation of several genes coding for fatty acid metabolism enzymes has previously been reported in male northern pike (Esox lucius) collected in the St. Lawrence River (QC, Canada) downstream of a major primary wastewater treatment plant (WWTP) point of discharge. The main objective of this study was to investigate the effects of exposure to the Montreal's WWTP effluent on the lipidomic profile (i.e., fatty acids, acylcarnitines, and phospholipids) as well as the transcription of genes related to lipid metabolism in the liver of northern pike collected upstream and downstream of this WWTP effluent. Halogenated flame retardant concentrations were also determined in pike liver and used as markers of exposure to this effluent. Greater concentrations of saturated and monounsaturated lysophosphatidylcholines (LPCs) and lower concentrations of polyunsaturated LPCs were determined in the liver of pike collected downstream of the WWTP compared to those collected upstream. Lower mRNA levels of peroxisome proliferator-activated receptor alpha (pparα), a major regulator of lipid metabolism, were also measured in pike exposed to Montreal's WWTP effluent. In addition, the relative contributions (%) of LPC 18:2 and LPC14:0, compounds used as markers of obesity and inflammation, were significantly correlated with halogenated flame retardant concentrations and fish girth. Results of the present study suggest that chronic environmental exposure to a primary WWTP effluent can modulate the transcription of genes related to lipid metabolism, and hence affect the hepatic phospholipid composition of pike from the St. Lawrence River.

Catch-related and genetic outcome of adult northern pike Esox lucius stocking in a large river system.

Genetic introgression from stocked adult northern pike Esox lucius to a wild self-recruiting population was detected in a large river system and some stocked E. lucius survived up to two spawning seasons and dispersed over several kilometres in the river. Moreover, the catch rate of stocked E. lucius by anglers was low (9.6%), hence suggesting that the efficiency of stocking activity is questionable.

Species richness in the African pike genus Hepsetus: a perfect match between genetics and morphology.

In this study, morphological differences were found that corroborate earlier results that showed the existence of six species within the previously monospecific African pike genus Hepsetus. Additional genetic data (coI, mtDNA and rag1, nDNA) confirm the morphology-based species delineations. Deep genetic divergences imply a relatively old age for diversification within the genus. An identification key for the six species is provided in the present study.

Higher-order organisation of extremely amplified, potentially functional and massively methylated 5S rDNA in European pikes (Esox sp.).

BACKGROUND: Pikes represent an important genus (Esox) harbouring a pre-duplication karyotype (2n = 2x = 50) of economically important salmonid pseudopolyploids. Here, we have characterized the 5S ribosomal RNA genes (rDNA) in Esox lucius and its closely related E. cisalpinus using cytogenetic, molecular and genomic approaches. Intragenomic homogeneity and copy number estimation was carried out using Illumina reads. The higher-order structure of rDNA arrays was investigated by the analysis of long PacBio reads. Position of loci on chromosomes was determined by FISH. DNA methylation was analysed by methylation-sensitive restriction enzymes. RESULTS: The 5S rDNA loci occupy exclusively (peri)centromeric regions on 30-38 acrocentric chromosomes in both E. lucius and E. cisalpinus. The large number of loci is accompanied by extreme amplification of genes (>20,000 copies), which is to the best of our knowledge one of the highest copy number of rRNA genes in animals ever reported. Conserved secondary structures of predicted 5S rRNAs indicate that most of the amplified genes are potentially functional. Only few SNPs were found in genic regions indicating their high homogeneity while intergenic spacers were more heterogeneous and several families were identified. Analysis of 10-30 kb-long molecules sequenced by the PacBio technology (containing about 40% of total 5S rDNA) revealed that the vast majority (96%) of genes are organised in large several kilobase-long blocks. Dispersed genes or short tandems were less common (4%). The adjacent 5S blocks were directly linked, separated by intervening DNA and even inverted. The 5S units differing in the intergenic spacers formed both homogeneous and heterogeneous (mixed) blocks indicating variable degree of homogenisation between the loci. Both E. lucius and E. cisalpinus 5S rDNA was heavily methylated at CG dinucleotides. CONCLUSIONS: Extreme amplification of 5S rRNA genes in the Esox genome occurred in the absence of significant pseudogenisation suggesting its recent origin and/or intensive homogenisation processes. The dense methylation of units indicates that powerful epigenetic mechanisms have evolved in this group of fish to silence amplified genes. We discuss how the higher-order repeat structures impact on homogenisation of 5S rDNA in the genome.

The axon degeneration gene SARM1 is evolutionarily distinct from other TIR domain-containing proteins.

Many forms of neurodegenerative disease are characterized by Wallerian degeneration, an active program of axonal destruction. Recently, the important player which enacts Wallerian degeneration was discovered, the multidomain protein SARM1. Since the SARM1 protein has classically been thought of as an innate immune molecule, its role in Wallerian degeneration has raised questions on the evolutionary forces acting on it. Here, we synthesize a picture of SARM1's evolution through various organisms by examining the molecular and genetic changes of SARM1 and the genes around it. Using proteins that possess domains homologous to SARM1, we established distances and Ka/Ks values through 5671 pairwise species-species comparisons. We demonstrate that SARM1 diverged across species in a pattern similar to other SAM domain-containing proteins. This is surprising, because it was expected that SARM1 would behave more like its TIR domain relatives. Going along with this divorce from TIR, we also noted that SARM1's TIR is under stronger purifying selection than the rest of the TIR domain-containing proteins (remaining highly conserved). In addition, SARM1's synteny analysis reveals that the surrounding gene cluster is highly conserved, functioning as a potential nexus of gene functionality across species. Taken together, SARM1 demonstrates a unique evolutionary pattern, separate from the TIR domain protein family.

Potential of Environmental DNA to Evaluate Northern Pike (Esox lucius) Eradication Efforts: An Experimental Test and Case Study.

Determining the success of invasive species eradication efforts is challenging because populations at very low abundance are difficult to detect. Environmental DNA (eDNA) sampling has recently emerged as a powerful tool for detecting rare aquatic animals; however, detectable fragments of DNA can persist over time despite absence of the targeted taxa and can therefore complicate eDNA sampling after an eradication event. This complication is a large concern for fish eradication efforts in lakes since killed fish can sink to the bottom and slowly decay. DNA released from these carcasses may remain detectable for long periods. Here, we evaluated the efficacy of eDNA sampling to detect invasive Northern pike (Esox lucius) following piscicide eradication efforts in southcentral Alaskan lakes. We used field observations and experiments to test the sensitivity of our Northern pike eDNA assay and to evaluate the persistence of detectable DNA emitted from Northern pike carcasses. We then used eDNA sampling and traditional sampling (i.e., gillnets) to test for presence of Northern pike in four lakes subjected to a piscicide-treatment designed to eradicate this species. We found that our assay could detect an abundant, free-roaming population of Northern pike and could also detect low-densities of Northern pike held in cages. For these caged Northern pike, probability of detection decreased with distance from the cage. We then stocked three lakes with Northern pike carcasses and collected eDNA samples 7, 35 and 70 days post-stocking. We detected DNA at 7 and 35 days, but not at 70 days. Finally, we collected eDNA samples ~ 230 days after four lakes were subjected to piscicide-treatments and detected Northern pike DNA in 3 of 179 samples, with a single detection at each of three lakes, though we did not catch any Northern pike in gillnets. Taken together, we found that eDNA can help to inform eradication efforts if used in conjunction with multiple lines of inquiry and sampling is delayed long enough to allow full degradation of DNA in the water.

Testing for Local Adaptation to Spawning Habitat in Sympatric Subpopulations of Pike by Reciprocal Translocation of Embryos.

We tested for local adaption in early life-history traits by performing a reciprocal translocation experiment with approximately 2,500 embryos of pike (Esox lucius) divided in paired split-family batches. The experiment indicated local adaptation in one of the two subpopulations manifested as enhanced hatching success of eggs in the native habitat, both when compared to siblings transferred to a non-native habitat, and when compared to immigrant genotypes from the other subpopulation. Gene-by-environment effects on viability of eggs and larvae were evident in both subpopulations, showing that there existed genetic variation allowing for evolutionary responses to divergent selection, and indicating a capacity for plastic responses to environmental change. Next, we tested for differences in female life-history traits. Results uncovered that females from one population invested more resources into reproduction and also produced more (but smaller) eggs in relation to their body size compared to females from the other population. We suggest that these females have adjusted their reproductive strategies as a counter-adaptation because a high amount of sedimentation on the eggs in that subpopulations spawning habitat might benefit smaller eggs. Collectively, our findings point to adaptive divergence among sympatric subpopulations that are physically separated only for a short period during reproduction and early development-which is rare. These results illustrate how combinations of translocation experiments and field studies of life-history traits might infer about local adaptation and evolutionary divergence among populations. Local adaptations in subdivided populations are important to consider in management and conservation of biodiversity, because they may otherwise be negatively affected by harvesting, supplementation, and reintroduction efforts targeted at endangered populations.

Causes and consequences of intra-specific variation in vertebral number.

Intraspecific variation in vertebral number is taxonomically widespread. Much scientific attention has been directed towards understanding patterns of variation in vertebral number among individuals and between populations, particularly across large spatial scales and in structured environments. However, the relative role of genes, plasticity, selection, and drift as drivers of individual variation and population differentiation remains unknown for most systems. Here, we report on patterns, causes and consequences of variation in vertebral number among and within sympatric subpopulations of pike (Esox lucius). Vertebral number differed among subpopulations, and common garden experiments indicated that this reflected genetic differences. A QST-FST comparison suggested that population differences represented local adaptations driven by divergent selection. Associations with fitness traits further indicated that vertebral counts were influenced both by stabilizing and directional selection within populations. Overall, our study enhances the understanding of adaptive variation, which is critical for the maintenance of intraspecific diversity and species conservation.

Mitochondrial genome of Esox flaviae (Southern pike): announcement and comparison with other Esocidae.

Pikes are fish species belonging to order Esociformes, family Esocidae, genus Esox. Species of the genus Esox are common, large, and economically important for food and fishing activities. Recently, a new species, southern pike E. flaviae, was described for a well-studied area such as Italy, using also two mtDNA markers: cox1 and cytb. A scant number of remnant populations of the species persist in Italy, threatened by habitat loss and degradation and by competition and possible hybridization with E. lucius, massively and recurrently stocked to sustain angling pressure. The availability of new mtDNA markers will possibly contribute to the conservation of the species. Currently, whole mitogenome information for the genus is available only for E. lucius and for E. reichertii. The aim of the present paper is to report novel mitogenomic information for southern pike.

Evolutionary divergence of adult body size and juvenile growth in sympatric subpopulations of a top predator in aquatic ecosystems.

Evolutionary theory predicts that different selective regimes may contribute to divergent evolution of body size and growth rate among populations, but most studies have focused on allopatric populations. Here, we studied five sympatric subpopulations of anadromous northern pike (Esox lucius) in the Baltic Sea subjected to allopatric habitats for a short period of their life cycle due to homing behavior. We report differences in adult body size among subpopulations that were in part due to variation in growth rate. Body size of emigrating juveniles also differed among subpopulations, and differences remained when individuals were reared in a common environment, thus indicating evolutionary divergence among subpopulations. Furthermore, a QST-FST comparison indicated that differences had evolved due to divergent selection rather than genetic drift, possibly in response to differences in selective mortality among spawning habitats during the allopatric life stage. Adult and juvenile size were negatively correlated across subpopulations, and reconstruction of growth trajectories of adult fishes suggested that body size differences developed gradually and became accentuated throughout the first years of life. These results represent rare evidence that sympatric subpopulations can evolve differences in key life-history traits despite being subjected to allopatric habitats during only a very short fraction of their life.

Individual heterogeneity in life histories and eco-evolutionary dynamics.

Individual heterogeneity in life history shapes eco-evolutionary processes, and unobserved heterogeneity can affect demographic outputs characterising life history and population dynamical properties. Demographic frameworks like matrix models or integral projection models represent powerful approaches to disentangle mechanisms linking individual life histories and population-level processes. Recent developments have provided important steps towards their application to study eco-evolutionary dynamics, but so far individual heterogeneity has largely been ignored. Here, we present a general demographic framework that incorporates individual heterogeneity in a flexible way, by separating static and dynamic traits (discrete or continuous). First, we apply the framework to derive the consequences of ignoring heterogeneity for a range of widely used demographic outputs. A general conclusion is that besides the long-term growth rate lambda, all parameters can be affected. Second, we discuss how the framework can help advance current demographic models of eco-evolutionary dynamics, by incorporating individual heterogeneity. For both applications numerical examples are provided, including an empirical example for pike. For instance, we demonstrate that predicted demographic responses to climate warming can be reversed by increased heritability. We discuss how applications of this demographic framework incorporating individual heterogeneity can help answer key biological questions that require a detailed understanding of eco-evolutionary dynamics.

Complete mitochondrial genome of Esox reichertii (Amur pike).

The whole mitochondrial genome of Esox reichertii (fish) was first sequenced and characterized. It was determined to be 16,909 bp long, which contains the control region (CR), the origin of light-strand replication (OL), 22 transfer RNA genes, 2 ribosomal genes and 13 protein-coding genes. Overall base composition of the complete mitochondrial DNA was 28.65% A, 28.67% T, 27.28% C and 15.41% G, with 57.32%AT.

The genome and linkage map of the northern pike (Esox lucius): conserved synteny revealed between the salmonid sister group and the Neoteleostei.

The northern pike is the most frequently studied member of the Esociformes, the closest order to the diverse and economically important Salmoniformes. The ancestor of all salmonids purportedly experienced a whole-genome duplication (WGD) event, making salmonid species ideal for studying the early impacts of genome duplication while complicating their use in wider analyses of teleost evolution. Studies suggest that the Esociformes diverged from the salmonid lineage prior to the WGD, supporting the use of northern pike as a pre-duplication outgroup. Here we present the first genome assembly, reference transcriptome and linkage map for northern pike, and evaluate the suitability of this species to provide a representative pre-duplication genome for future studies of salmonid and teleost evolution. The northern pike genome sequence is composed of 94,267 contigs (N50 = 16,909 bp) contained in 5,688 scaffolds (N50 = 700,535 bp); the total scaffolded genome size is 878 million bases. Multiple lines of evidence suggest that over 96% of the protein-coding genome is present in the genome assembly. The reference transcriptome was constructed from 13 tissues and contains 38,696 transcripts, which are accompanied by normalized expression data in all tissues. Gene-prediction analysis produced a total of 19,601 northern pike-specific gene models. The first-generation linkage map identifies 25 linkage groups, in agreement with northern pike's diploid karyotype of 2N = 50, and facilitates the placement of 46% of assembled bases onto linkage groups. Analyses reveal a high degree of conserved synteny between northern pike and other model teleost genomes. While conservation of gene order is limited to smaller syntenic blocks, the wider conservation of genome organization implies the northern pike exhibits a suitable approximation of a non-duplicated Protacanthopterygiian genome. This dataset will facilitate future studies of esocid biology and empower ongoing examinations of the Atlantic salmon and rainbow trout genomes by facilitating their comparison with other major teleost groups.

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.