Cell atlas of the Atlantic salmon spleen reveals immune cell heterogeneity and cell-specific responses to bacterial infection.

The spleen is a conserved secondary lymphoid organ that emerged in parallel to adaptive immunity in early jawed vertebrates. Recent studies have applied single cell transcriptomics to reveal the cellular composition of spleen in several species, cataloguing diverse immune cell types and subpopulations. In this study, 51,119 spleen nuclei transcriptomes were comprehensively investigated in the commercially important teleost Atlantic salmon (Salmo salar L.), contrasting control animals with those challenged with the bacterial pathogen Aeromonas salmonicida. We identified clusters of nuclei representing the expected major cell types, namely T cells, B cells, natural killer-like cells, granulocytes, mononuclear phagocytes, endothelial cells, mesenchymal cells, erythrocytes and thrombocytes. We discovered heterogeneity within several immune lineages, providing evidence for resident macrophages and melanomacrophages, infiltrating monocytes, several candidate dendritic cell subpopulations, and B cells at distinct stages of differentiation, including plasma cells and an igt + subset. We provide evidence for twelve candidate T cell subsets, including cd4+ T helper and regulatory T cells, one cd8+ subset, three γδT subsets, and populations double negative for cd4 and cd8. The number of genes showing differential expression during the early stages of Aeromonas infection was highly variable across immune cell types, with the largest changes observed in macrophages and infiltrating monocytes, followed by resting mature B cells. Our analysis provides evidence for a local inflammatory response to infection alongside B cell maturation in the spleen, and upregulation of ccr9 genes in igt + B cells, T helper and cd8+ cells, and monocytes, consistent with the recruitment of immune cell populations to the gut to deal with Aeromonas infection. Overall, this study provides a new cell-resolved perspective of the immune actions of Atlantic salmon spleen, highlighting extensive heterogeneity hidden to bulk transcriptomics. We further provide a large catalogue of cell-specific marker genes that can be leveraged to further explore the function and structural organization of the salmonid immune system.

Mapping the cellular landscape of Atlantic salmon head kidney by single cell and single nucleus transcriptomics.

Single-cell transcriptomics is the current gold standard for global gene expression profiling, not only in mammals and model species, but also in non-model fish species. This is a rapidly expanding field, creating a deeper understanding of tissue heterogeneity and the distinct functions of individual cells, making it possible to explore the complexities of immunology and gene expression on a highly resolved level. In this study, we compared two single cell transcriptomic approaches to investigate cellular heterogeneity within the head kidney of healthy farmed Atlantic salmon (Salmo salar). We compared 14,149 cell transcriptomes assayed by single cell RNA-seq (scRNA-seq) with 18,067 nuclei transcriptomes captured by single nucleus RNA-Seq (snRNA-seq). Both approaches detected eight major cell populations in common: granulocytes, heamatopoietic stem cells, erythrocytes, mononuclear phagocytes, thrombocytes, B cells, NK-like cells, and T cells. Four additional cell types, endothelial, epithelial, interrenal, and mesenchymal cells, were detected in the snRNA-seq dataset, but appeared to be lost during preparation of the single cell suspension submitted for scRNA-seq library generation. We identified additional heterogeneity and subpopulations within the B cells, T cells, and endothelial cells, and revealed developmental trajectories of heamatopoietic stem cells into differentiated granulocyte and mononuclear phagocyte populations. Gene expression profiles of B cell subtypes revealed distinct IgM and IgT-skewed resting B cell lineages and provided insights into the regulation of B cell lymphopoiesis. The analysis revealed eleven T cell sub-populations, displaying a level of T cell heterogeneity in salmon head kidney comparable to that observed in mammals, including distinct subsets of cd4/cd8-negative T cells, such as tcrγ positive, progenitor-like, and cytotoxic cells. Although snRNA-seq and scRNA-seq were both useful to resolve cell type-specific expression in the Atlantic salmon head kidney, the snRNA-seq pipeline was overall more robust in identifying several cell types and subpopulations. While scRNA-seq displayed higher levels of ribosomal and mitochondrial genes, snRNA-seq captured more transcription factor genes. However, only scRNA-seq-generated data was useful for cell trajectory inference within the myeloid lineage. In conclusion, this study systematically outlines the relative merits of scRNA-seq and snRNA-seq in Atlantic salmon, enhances understanding of teleost immune cell lineages, and provides a comprehensive list of markers for identifying major cell populations in the head kidney with significant immune relevance.

A versatile nuclei extraction protocol for single nucleus sequencing in non-model species-Optimization in various Atlantic salmon tissues.

The use of single cell sequencing technologies has exploded over recent years, and is now commonly used in many non-model species. Sequencing nuclei instead of whole cells has become increasingly popular, as it does not require the processing of samples immediately after collection. Here we present a highly effective nucleus isolation protocol that outperforms previously available method in challenging samples in a non-model specie. This protocol can be successfully applied to extract nuclei from a variety of tissues and species.

Conserved and divergent arms of the antiviral response in the duplicated genomes of salmonid fishes.

Antiviral innate immunity is orchestrated by the interferon system, which appeared in ancestors of jawed vertebrates. Interferon upregulation induces hundreds of interferon-stimulated-genes (ISGs) with effector or regulatory functions. Here we investigated the evolutionary diversification of ISG responses through comparison of two salmonid fishes, accounting for the impact of sequential whole genome duplications ancestral to teleosts and salmonids. We analysed the transcriptomic response of the IFN pathway in the head kidney of rainbow trout and Atlantic salmon, which separated 25-30 Mya. We identified a large set of ISGs conserved in both species and cross-referenced them with zebrafish and human ISGs. In contrast, around one-third of salmonid ISG lacked orthologs in human, mouse, chicken or frog, and often between rainbow trout and Atlantic salmon, revealing a fast-evolving, lineage-specific arm of the antiviral response. This study also provides a key resource for in-depth functional analysis of ISGs in salmonids of commercial significance.

Organized B cell sites in cartilaginous fishes reveal the evolutionary foundation of germinal centers.

The absence of germinal centers (GCs) in cartilaginous fishes lies at odds with data showing that nurse sharks can produce robust antigen-specific responses and affinity mature their B cell repertoires. To investigate this apparent incongruity, we performed RNA sequencing on single nuclei, allowing us to characterize the cell types present in the nurse shark spleen, and RNAscope to provide in situ cellular resolution of key marker gene expression following immunization with R-phycoerythrin (PE). We tracked PE to the splenic follicles where it co-localizes with CXCR5high centrocyte-like B cells and a population of putative T follicular helper (Tfh) cells, surrounded by a peripheral ring of Ki67+ AID+ CXCR4+ centroblast-like B cells. Further, we reveal selection of mutations in B cell clones dissected from these follicles. We propose that the B cell sites identified here represent the evolutionary foundation of GCs, dating back to the jawed vertebrate ancestor.

Independent rediploidization masks shared whole genome duplication in the sturgeon-paddlefish ancestor.

Whole genome duplication (WGD) is a dramatic evolutionary event generating many new genes and which may play a role in survival through mass extinctions. Paddlefish and sturgeon are sister lineages that both show genomic evidence for ancient WGD. Until now this has been interpreted as two independent WGD events due to a preponderance of duplicate genes with independent histories. Here we show that although there is indeed a plurality of apparently independent gene duplications, these derive from a shared genome duplication event occurring well over 200 million years ago, likely close to the Permian-Triassic mass extinction period. This was followed by a prolonged process of reversion to stable diploid inheritance (rediploidization), that may have promoted survival during the Triassic-Jurassic mass extinction. We show that the sharing of this WGD is masked by the fact that paddlefish and sturgeon lineage divergence occurred before rediploidization had proceeded even half-way. Thus, for most genes the resolution to diploidy was lineage-specific. Because genes are only truly duplicated once diploid inheritance is established, the paddlefish and sturgeon genomes are thus a mosaic of shared and non-shared gene duplications resulting from a shared genome duplication event.

Single cell genomics as a transformative approach for aquaculture research and innovation.

Single cell genomics encompasses a suite of rapidly maturing technologies that measure the molecular profiles of individual cells within target samples. These approaches provide a large up-step in biological information compared to long-established 'bulk' methods that profile the average molecular profiles of all cells in a sample, and have led to transformative advances in understanding of cellular biology, particularly in humans and model organisms. The application of single cell genomics is fast expanding to non-model taxa, including aquaculture species, where numerous research applications are underway with many more envisaged. In this review, we highlight the potential transformative applications of single cell genomics in aquaculture research, considering barriers and potential solutions to the broad uptake of these technologies. Focusing on single cell transcriptomics, we outline considerations for experimental design, including the essential requirement to obtain high quality cells/nuclei for sequencing in ectothermic aquatic species. We further outline data analysis and bioinformatics considerations, tailored to studies with the under-characterized genomes of aquaculture species, where our knowledge of cellular heterogeneity and cell marker genes is immature. Overall, this review offers a useful source of knowledge for researchers aiming to apply single cell genomics to address biological challenges faced by the global aquaculture sector though an improved understanding of cell biology.

Two parallel chromosome-level reference genomes to support restoration and aquaculture of European flat oyster Ostrea edulis.

This volume of Evolutionary Applications sees the publication of two genomes for the European native flat oyster Ostrea edulis, a species of significant evolutionary, ecological and commercial value. Each is a highly contiguous chromosome-level assembly from individuals of different genetic backgrounds, which have been benchmarked against one another. This situation has resulted from the serendipitous discovery that two independent research groups were both deep into the process of building, annotating and investigating separately produced assemblies. Due to constraints with funder requirements and the need to recognize early career researchers for their work, alongside the technical challenge of integrating assemblies from two very different genomes, there was limited capacity to merge the sequences into one publication at the stage of discovery. This issue is likely to become very common over the next few years until the technologies for working with multiple genomes at once, for example, graph genomes, become commonplace in nonmodel species. Consequently, both of our teams have decided to collaborate rather than compete, recognizing the benefit to copublishing two separate genome resources for the research community, each with distinct scientific investigations, and working collaboratively to benchmark the assemblies.

Chromosome-level reference genome for European flat oyster (Ostrea edulis L.).

The European flat oyster (Ostrea edulis L.) is a bivalve naturally distributed across Europe, which was an integral part of human diets for centuries, until anthropogenic activities and disease outbreaks severely reduced wild populations. Despite a growing interest in genetic applications to support population management and aquaculture, a reference genome for this species is lacking to date. Here, we report a chromosome-level assembly and annotation for the European Flat oyster genome, generated using Oxford Nanopore, Illumina, Dovetail OmniC™ proximity ligation and RNA sequencing. A contig assembly (N50: 2.38 Mb) was scaffolded into the expected karyotype of 10 pseudochromosomes. The final assembly is 935.13 Mb, with a scaffold-N50 of 95.56 Mb, with a predicted repeat landscape dominated by unclassified elements specific to O. edulis. The assembly was verified for accuracy and completeness using multiple approaches, including a novel linkage map built with ddRAD-Seq technology, comprising 4016 SNPs from four full-sib families (eight parents and 163 F1 offspring). Annotation of the genome integrating multitissue transcriptome data, comparative protein evidence and ab-initio gene prediction identified 35,699 protein-coding genes. Chromosome-level synteny was demonstrated against multiple high-quality bivalve genome assemblies, including an O. edulis genome generated independently for a French O. edulis individual. Comparative genomics was used to characterize gene family expansions during Ostrea evolution that potentially facilitated adaptation. This new reference genome for European flat oyster will enable high-resolution genomics in support of conservation and aquaculture initiatives, and improves our understanding of bivalve genome evolution.

A single genomic region involving a putative chromosome rearrangement in flat oyster (Ostrea edulis) is associated with differential host resilience to the parasite Bonamia ostreae.

European flat oyster (Ostrea edulis) is an ecologically and economically important marine bivalve, that has been severely affected by the intracellular parasite Bonamia ostreae. In this study, a flat oyster SNP array (~14,000 SNPs) was used to validate previously reported outlier loci for divergent selection associated with B. ostreae exposure in the Northeast Atlantic Area. A total of 134 wild and hatchery individuals from the North Sea, collected in naïve (NV) and long-term affected (LTA) areas, were analysed. Genetic diversity and differentiation were related to the sampling origin (wild vs. hatchery) when using neutral markers, and to bonamiosis status (NV vs. LTA) when using outlier loci for divergent selection. Two genetic clusters appeared intermingled in all sampling locations when using outlier loci, and their frequency was associated with their bonamiosis status. When both clusters were compared, outlier data sets showed high genetic divergence (F ST > 0.25) unlike neutral loci (F ST not ≠ 0). Moreover, the cluster associated with LTA samples showed much higher genetic diversity and significant heterozygote excess with outlier loci, but not with neutral data. Most outliers mapped on chromosome 8 (OE-C8) of the flat oyster genome, supporting a main genomic region underlying resilience to bonamiosis. Furthermore, differentially expressed genes previously reported between NV and LTA strains showed higher mapping density on OE-C8. A range of relevant immune functions were specifically enriched among genes annotated on OE-C8, providing hypotheses for resilience mechanisms to an intracellular parasite. The results suggest that marker-assisted selection could be applied to breed resilient strains of O. edulis to bonamiosis, if lower parasite load and/or higher viability of the LTA genetic cluster following B. ostreae infection is demonstrated.

Single cell transcriptomics of Atlantic salmon (Salmo salar L.) liver reveals cellular heterogeneity and immunological responses to challenge by Aeromonas salmonicida.

The liver is a multitasking organ with essential functions for vertebrate health spanning metabolism and immunity. In contrast to mammals, our understanding of liver cellular heterogeneity and its role in regulating immunological status remains poorly defined in fishes. Addressing this knowledge gap, we generated a transcriptomic atlas of 47,432 nuclei isolated from the liver of Atlantic salmon (Salmo salar L.) contrasting control fish with those challenged with a pathogenic strain of Aeromonas salmonicida, a problematic bacterial pathogen in global aquaculture. We identified the major liver cell types and their sub-populations, revealing poor conservation of many hepatic cell marker genes utilized in mammals, while identifying novel heterogeneity within the hepatocyte, lymphoid, and myeloid lineages. This included polyploid hepatocytes, multiple T cell populations including γδ T cells, and candidate populations of monocytes/macrophages and dendritic cells. A dominant hepatocyte population radically remodeled its transcriptome following infection to activate the acute phase response and other defense functions, while repressing routine functions such as metabolism. These defense-specialized hepatocytes showed strong activation of genes controlling protein synthesis and secretion, presumably to support the release of acute phase proteins into circulation. The infection response further involved up-regulation of numerous genes in an immune-cell specific manner, reflecting functions in pathogen recognition and killing, antigen presentation, phagocytosis, regulation of inflammation, B cell differentiation and T cell activation. Overall, this study greatly enhances our understanding of the multifaceted role played by liver immune and non-immune cells in host defense and metabolic remodeling following infection and provides many novel cell-specific marker genes to empower future studies of this organ in fishes.

Genome-Wide Association and Genomic Prediction of Growth Traits in the European Flat Oyster (Ostrea edulis).

The European flat oyster (Ostrea edulis) is a bivalve mollusc that was once widely distributed across Europe and represented an important food resource for humans for centuries. Populations of O. edulis experienced a severe decline across their biogeographic range mainly due to overexploitation and disease outbreaks. To restore the economic and ecological benefits of European flat oyster populations, extensive protection and restoration efforts are in place within Europe. In line with the increasing interest in supporting restoration and oyster farming through the breeding of stocks with enhanced performance, the present study aimed to evaluate the potential of genomic selection for improving growth traits in a European flat oyster population obtained from successive mass-spawning events. Four growth-related traits were evaluated: total weight (TW), shell height (SH), shell width (SW) and shell length (SL). The heritability of the growth traits was in the low-moderate range, with estimates of 0.45, 0.37, 0.22, and 0.32 for TW, SH, SW and SL, respectively. A genome-wide association analysis revealed a largely polygenic architecture for the four growth traits, with two distinct QTLs detected on chromosome 4. To investigate whether genomic selection can be implemented in flat oyster breeding at a reduced cost, the utility of low-density SNP panels was assessed. Genomic prediction accuracies using the full density panel were high (> 0.83 for all traits). The evaluation of the effect of reducing the number of markers used to predict genomic breeding values revealed that similar selection accuracies could be achieved for all traits with 2K SNPs as for a full panel containing 4,577 SNPs. Only slight reductions in accuracies were observed at the lowest SNP density tested (i.e., 100 SNPs), likely due to a high relatedness between individuals being included in the training and validation sets during cross-validation. Overall, our results suggest that the genetic improvement of growth traits in oysters is feasible. Nevertheless, and although low-density SNP panels appear as a promising strategy for applying GS at a reduced cost, additional populations with different degrees of genetic relatedness should be assessed to derive estimates of prediction accuracies to be expected in practical breeding programmes.

Exploration of the Nurse Shark (Ginglymostoma cirratum) Plasma Immunoproteome Using High-Resolution LC-MS/MS.

Many animals of scientific importance lack species-specific reagents (e.g., monoclonal antibodies) for in-depth studies of immune proteins. Mass spectrometry (MS)-based proteomics has emerged as a useful method for monitoring changes in protein abundance and modifications in non-model species. It can be used to quantify hundreds of candidate immune molecules simultaneously without the generation of new reagents. Here, we used MS-based proteomics to identify and quantify candidate immune proteins in the plasma of the nurse shark (Ginglymostoma cirratum), a cartilaginous fish and representative of the most basal extant vertebrate lineage with an immunoglobulin-based immune system. Mass spectrometry-based LC-MS/MS was performed on the blood plasma of nurse sharks immunized with human serum albumin (n=4) or sham immunized (n=1), and sampled at days 0 (baseline control), 1, 2, 3, 5, 7, 14, 21, 28, 25, 42 and 49. An antigen-specific antibody response was experimentally confirmed post-immunization. To provide a high-quality reference to identify proteins, we assembled and annotated a multi-tissue de novo transcriptome integrating long- and short-read sequence data. This comprised 62,682 contigs containing open reading frames (ORFs) with a length >80 amino acids. Using this transcriptome, we reliably identified 626 plasma proteins which were broadly categorized into coagulation, immune, and metabolic functional groups. To assess the feasibility of performing LC-MS/MS proteomics in nurse shark in the absence of species-specific protein annotations, we compared the results to an alternative strategy, mapping peptides to proteins predicted in the genome assembly of a related species, the whale shark (Rhincodon typus). This approach reliably identified 297 proteins, indicating that useful data on the plasma proteome may be obtained in many instances despite the absence of a species-specific reference protein database. Among the plasma proteins defined against the nurse shark transcriptome, fifteen showed consistent changes in abundance across the immunized shark individuals, indicating a role in the immune response. These included alpha-2-macroglobulin (A2M) and a novel protein yet to be characterized in diverse vertebrate lineages. Overall, this study enhances genetic and protein-level resources for nurse shark research and vastly improves our understanding of the elasmobranch plasma proteome, including its remodelling following immune stimulation.

Genome-Wide Reconstruction of Rediploidization Following Autopolyploidization across One Hundred Million Years of Salmonid Evolution.

The long-term evolutionary impacts of whole-genome duplication (WGD) are strongly influenced by the ensuing rediploidization process. Following autopolyploidization, rediploidization involves a transition from tetraploid to diploid meiotic pairing, allowing duplicated genes (ohnologs) to diverge genetically and functionally. Our understanding of autopolyploid rediploidization has been informed by a WGD event ancestral to salmonid fishes, where large genomic regions are characterized by temporally delayed rediploidization, allowing lineage-specific ohnolog sequence divergence in the major salmonid clades. Here, we investigate the long-term outcomes of autopolyploid rediploidization at genome-wide resolution, exploiting a recent "explosion" of salmonid genome assemblies, including a new genome sequence for the huchen (Hucho hucho). We developed a genome alignment approach to capture duplicated regions across multiple species, allowing us to create 121,864 phylogenetic trees describing genome-wide ohnolog divergence across salmonid evolution. Using molecular clock analysis, we show that 61% of the ancestral salmonid genome experienced an initial "wave" of rediploidization in the late Cretaceous (85-106 Ma). This was followed by a period of relative genomic stasis lasting 17-39 My, where much of the genome remained tetraploid. A second rediploidization wave began in the early Eocene and proceeded alongside species diversification, generating predictable patterns of lineage-specific ohnolog divergence, scaling in complexity with the number of speciation events. Using gene set enrichment, gene expression, and codon-based selection analyses, we provide insights into potential functional outcomes of delayed rediploidization. This study enhances our understanding of delayed autopolyploid rediploidization and has broad implications for future studies of WGD events.

Genomic Epidemiology of Salmonid Alphavirus in Norwegian Aquaculture Reveals Recent Subtype-2 Transmission Dynamics and Novel Subtype-3 Lineages.

Viral disease poses a major barrier to sustainable aquaculture, with outbreaks causing large economic losses and growing concerns for fish welfare. Genomic epidemiology can support disease control by providing rapid inferences on viral evolution and disease transmission. In this study, genomic epidemiology was used to investigate salmonid alphavirus (SAV), the causative agent of pancreas disease (PD) in Atlantic salmon. Our aim was to reconstruct SAV subtype-2 (SAV2) diversity and transmission dynamics in recent Norwegian aquaculture, including the origin of SAV2 in regions where this subtype is not tolerated under current legislation. Using nanopore sequencing, we captured ~90% of the SAV2 genome for n = 68 field isolates from 10 aquaculture production regions sampled between 2018 and 2020. Using time-calibrated phylogenetics, we infer that, following its introduction to Norway around 2010, SAV2 split into two clades (SAV2a and 2b) around 2013. While co-present at the same sites near the boundary of Møre og Romsdal and Trøndelag, SAV2a and 2b were generally detected in non-overlapping locations at more Southern and Northern latitudes, respectively. We provide evidence for recent SAV2 transmission over large distances, revealing a strong connection between Møre og Romsdal and SAV2 detected in 2019/20 in Rogaland. We also demonstrate separate introductions of SAV2a and 2b outside the SAV2 zone in Sognefjorden (Vestland), connected to samples from Møre og Romsdal and Trøndelag, respectively, and a likely 100 km Northward transmission of SAV2b within Trøndelag. Finally, we recovered genomes of SAV2a and SAV3 co-infecting single fish in Rogaland, involving novel SAV3 lineages that diverged from previously characterized strains >25 years ago. Overall, this study demonstrates useful applications of genomic epidemiology for tracking viral disease spread in aquaculture.

Comparative regulomics supports pervasive selection on gene dosage following whole genome duplication.

BACKGROUND: Whole genome duplication (WGD) events have played a major role in eukaryotic genome evolution, but the consequence of these extreme events in adaptive genome evolution is still not well understood. To address this knowledge gap, we used a comparative phylogenetic model and transcriptomic data from seven species to infer selection on gene expression in duplicated genes (ohnologs) following the salmonid WGD 80-100 million years ago. RESULTS: We find rare cases of tissue-specific expression evolution but pervasive expression evolution affecting many tissues, reflecting strong selection on maintenance of genome stability following genome doubling. Ohnolog expression levels have evolved mostly asymmetrically, by diverting one ohnolog copy down a path towards lower expression and possible pseudogenization. Loss of expression in one ohnolog is significantly associated with transposable element insertions in promoters and likely driven by selection on gene dosage including selection on stoichiometric balance. We also find symmetric expression shifts, and these are associated with genes under strong evolutionary constraints such as ribosome subunit genes. This possibly reflects selection operating to achieve a gene dose reduction while avoiding accumulation of "toxic mutations". Mechanistically, ohnolog regulatory divergence is dictated by the number of bound transcription factors in promoters, with transposable elements being one likely source of novel binding sites driving tissue-specific gains in expression. CONCLUSIONS: Our results imply pervasive adaptive expression evolution following WGD to overcome the immediate challenges posed by genome doubling and to exploit the long-term genetic opportunities for novel phenotype evolution.

Nanopore whole genome sequencing and partitioned phylogenetic analysis supports a new salmonid alphavirus genotype (SAV7).

Salmon pancreas disease virus, more commonly known as salmonid alphavirus (SAV), is a single-stranded positive sense RNA virus and the causative agent of pancreas disease and sleeping disease in salmonids. In this study, a unique strain of SAV previously isolated from ballan wrasse was subjected to whole genome sequencing using nanopore sequencing. In order to accurately examine the evolutionary history of this strain in comparison to other SAV strains, a partitioned phylogenetic analysis was performed to account for variation in the rate of evolution for both individual genes and codon positions. Partitioning the genome alignments almost doubled the observed branch lengths in the phylogenetic tree when compared to the more common approach of applying one model of substitution across the genome and significantly increased the statistical fit of the best-fitting models of nucleotide substitution. Based on the genomic data, a valid case can be made for the viral strain examined in this study to be considered a new SAV genotype. In addition, this study adds to a growing number of studies in which SAV has been found to infect non-salmonid fish, and as such we have suggested that the viral species name be amended to the more inclusive 'piscine alphavirus'.

Plasma Proteome Responses in Salmonid Fish Following Immunization.

Vaccination plays a critical role in the protection of humans and other animals from infectious diseases. However, the same vaccine often confers different protection levels among individuals due to variation in genetics and/or immunological histories. While this represents a well-recognized issue in humans, it has received little attention in fish. Here we address this knowledge gap in a proteomic study of rainbow trout (Oncorhynchus mykiss, Walbaum), using non-lethal repeated blood sampling to establish the plasma protein response of individual fish following immunization. Six trout were immunized with adjuvanted hen egg-white lysozyme (HEL) and peripheral blood sampled at ten time points from day 0 to day 84 post-injection. We confirm that an antigen-specific antibody response to HEL was raised, showing differences in timing and magnitude among individuals. Using label-free liquid chromatography-mass spectrometry, we quantified the abundance of 278 plasma proteins across the timecourse. As part of the analysis, we show that this approach can distinguish many (but not all) duplicated plasma proteins encoded by paralogous genes retained from the salmonid-specific whole genome duplication event. Global variation in the plasma proteome was predominantly explained by individual differences among fish. However, sampling day explained a major component of variation in abundance for a statistically defined subset of 41 proteins, representing 15% of those detected. These proteins clustered into five groups showing distinct temporal responses to HEL immunization at the population level, and include classical immune (e.g. complement system members) and acute phase molecules (e.g. apolipoproteins, haptoglobins), several enzymes and other proteins supporting the immune response, in addition to evolutionarily conserved molecules that are as yet uncharacterized. Overall, this study improves our understanding of the fish plasma proteome, provides valuable marker proteins for different phases of the immune response, and has implications for vaccine development and the design of immune challenge experiments.

The structural variation landscape in 492 Atlantic salmon genomes.

Structural variants (SVs) are a major source of genetic and phenotypic variation, but remain challenging to accurately type and are hence poorly characterized in most species. We present an approach for reliable SV discovery in non-model species using whole genome sequencing and report 15,483 high-confidence SVs in 492 Atlantic salmon (Salmo salar L.) sampled from a broad phylogeographic distribution. These SVs recover population genetic structure with high resolution, include an active DNA transposon, widely affect functional features, and overlap more duplicated genes retained from an ancestral salmonid autotetraploidization event than expected. Changes in SV allele frequency between wild and farmed fish indicate polygenic selection on behavioural traits during domestication, targeting brain-expressed synaptic networks linked to neurological disorders in humans. This study offers novel insights into the role of SVs in genome evolution and the genetic architecture of domestication traits, along with resources supporting reliable SV discovery in non-model species.