Evolution of T cell receptor beta loci in salmonids

T-cell mediated immunity relies on a vast array of antigen specific T cell receptors (TR). Characterizing the structure of TR loci is essential to study the diversity and composition of T cell responses in vertebrate species. The lack of good-quality genome assemblies, and the difficulty to perform a reliably mapping of multiple highly similar TR sequences, have hindered the study of these loci in non-model organisms. High-quality genome assemblies are now available for the two main genera of Salmonids, Salmo and Oncorhynchus. We present here a full description and annotation of the TRB loci located on chromosomes 19 and 25 of rainbow trout (Oncorhynchus mykiss). To get insight about variations of the structure and composition of TRB locus across salmonids, we compared rainbow trout TRB loci with other salmonid species and confirmed that the basic structure of salmonid TRB locus is a double set of two TRBV-D-J-C loci in opposite orientation on two different chromosomes. Our data shed light on the evolution of TRB loci in Salmonids after their whole genome duplication (WGD). We established a coherent nomenclature of salmonid TRB loci based on comprehensive annotation. Our work provides a fundamental basis for monitoring salmonid T cell responses by TRB repertoire sequencing.

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.

Profiling the T Cell Receptor Alpha/Delta Locus in Salmonids

In jawed vertebrates, two major T cell populations have been characterized. They are defined as α/β or γ/δ T cells, based on the expressed T cell receptor. Salmonids (family Salmonidae) include two key teleost species for aquaculture, rainbow trout (Oncorhynchus mykiss) and Atlantic salmon (Salmo salar) which constitute important models for fish immunology and important targets for vaccine development. The growing interest to decipher the dynamics of adaptive immune responses against pathogens or vaccines has resulted in recent efforts to sequence the immunoglobulin (IG) or antibodies and T cell receptor (TR) repertoire in these species. In this context, establishing a comprehensive and coherent locus annotation is the fundamental basis for the analysis of high-throughput repertoire sequencing data. We therefore decided to revisit the description and annotation of TRA/TRD locus in Atlantic salmon and two strains of rainbow trout (Swanson and Arlee) using the now available high-quality genome assemblies. Phylogenetic analysis of functional TRA/TRD V genes from these three genomes led to the definition of 25 subgroups shared by both species, some with particular feature. A total of 128 TRAJ genes were identified in Salmo, the majority with a close counterpart in Oncorhynchus. Analysis of expressed TRA repertoire indicates that most TRAV gene subgroups are expressed at mucosal and systemic level. The present work on TRA/TRD locus annotation along with the analysis of TRA repertoire sequencing data show the feasibility and advantages of a common salmonid TRA/TRD nomenclature that allows an accurate annotation and analysis of high-throughput sequencing results, across salmonid T cell subsets.

Development of a Multiplex PCR Assay for Genotyping the Fish Pathogen Piscirickettsia salmonis Through Comparative Genomics

Piscirickettsia salmonis is a bacterial pathogen that severely impact the aquaculture in several countries as Canada, Scotland, Ireland, Norway, and Chile. It provokes Piscirickettsiosis outbreaks in the marine phase of salmonid farming, resulting in economic losses. The monophyletic genogroup LF-89 and a divergent genogroup EM-90 are responsible for the most severe Piscirickettsiosis outbreaks in Chile. Therefore, the development of methods for quick genotyping of P. salmonis genogroups in field samples is vital for veterinary diagnoses and understanding the population structure of this pathogen. The present study reports the development of a multiplex PCR for genotyping LF-89 and EM-90 genogroups based on comparative genomics of 73 fully sequenced P. salmonis genomes. The results revealed 2,322 sequences shared between 35 LF-89 genomes, 2,280 sequences in the core-genome of 38 EM-90 genomes, and 331 and 534 accessory coding sequences each genogroup, respectively. A total of 1,801 clusters of coding sequences were shared among all tested genomes of P. salmonis (LF-89 and EM-90), with 253 and 291 unique sequences for LF-89 and EM-90 genogroups, respectively. The Multiplex-1 prototype was chosen for reliable genotyping because of differences in annealing temperatures and respective reaction efficiencies. This method also identified the pathogen in field samples infected with LF-89 or EM-90 strains, which is not possible with other methods currently available. Finally, the genome-based multiplex PCR protocol presented in this study is a rapid and affordable alternative to classical sequencing of PCR products and analyzing the length of restriction fragment polymorphisms.

Evidence for Introgressive Hybridization between Native Dolly Varden (Salvelinus curilus (syn. Salvelinus malma)) and Introduced Brook Trout (Salvelinus fontinalis) in the Nishibetsu River of Hokkaido, Japan

Hybridization is one of the negative outcomes for the introduction of non-native species, which can lead to rapid displacement and genetic extinction of native species. Salmonid fishes have been widely introduced outside of their native ranges for food supply and recreational fishing. Here, we investigate the occurrence of introgressive hybridization among native Dolly Varden (Salvelinus curilus (syn. Salvelinus malma)), white-spotted charr (Salvelinus leucomaenis), and introduced brook trout (Salvelinus fontinalis), in streams of the Nishibetsu River, Hokkaido, Japan. Microsatellite DNA analysis detected five hybrids between native Dolly Varden and introduced brook trout. This is the first evidence for hybridization between native Dolly Varden and introduced brook trout, while the latter has been known to hybridize with many other salmonids. Furthermore, incongruence between mitochondrial DNA and microsatellite DNA analyses suggested introgression among the three Salvelinus species. Further studies to estimate the hybrid fitness are necessary to understand how hybridization among the three species affects the native species.

Know your enemy - transcriptome of myxozoan Tetracapsuloides bryosalmonae reveals potential drug targets against proliferative kidney disease in salmonids

The myxozoan Tetracapsuloides bryosalmonae is a widely spread endoparasite that causes proliferative kidney disease (PKD) in salmonid fish. We developed an in silico pipeline to separate transcripts of T. bryosalmonae from the kidney tissue of its natural vertebrate host, brown trout (Salmo trutta). After stringent filtering, we constructed a partial transcriptome assembly T. bryosalmonae, comprising 3427 transcripts. Based on homology-restricted searches of the assembled parasite transcriptome and Atlantic salmon (Salmo salar) proteome, we identified four protein targets (Endoglycoceramidase, Legumain-like protease, Carbonic anhydrase 2, Pancreatic lipase-related protein 2) for the development of anti-parasitic drugs against T. bryosalmonae. Earlier work of these proteins on parasitic protists and helminths suggests that the identified anti-parasitic drug targets represent promising chemotherapeutic candidates also against T. bryosalmonae, and strengthen the view that the known inhibitors can be effective in evolutionarily distant organisms. In addition, we identified differentially expressed T. bryosalmonae genes between moderately and severely infected fish, indicating an increased abundance of T. bryosalmonae sporogonic stages in fish with low parasite load. In conclusion, this study paves the way for future genomic research in T. bryosalmonae and represents an important step towards the development of effective drugs against PKD.

Standardized IMGT® Nomenclature of Salmonidae IGH Genes, the Paradigm of Atlantic Salmon and Rainbow Trout: From Genomics to Repertoires

In teleost fish as in mammals, humoral adaptive immunity is based on B lymphocytes expressing highly diverse immunoglobulins (IG). During B cell differentiation, IG loci are subjected to genomic rearrangements of V, D, and J genes, producing a unique antigen receptor expressed on the surface of each lymphocyte. During the course of an immune response to infections or immunizations, B cell clones specific of epitopes from the immunogen are expanded and activated, leading to production of specific antibodies. Among teleost fish, salmonids comprise key species for aquaculture. Rainbow trout (Oncorhynchus mykiss) and Atlantic salmon (Salmo salar) are especially important from a commercial point of view and have emerged as critical models for fish immunology. The growing interest to capture accurate and comprehensive antibody responses against common pathogens and vaccines has resulted in recent efforts to sequence the IG repertoire in these species. In this context, a unified and standardized nomenclature of salmonid IG heavy chain (IGH) genes is urgently required, to improve accuracy of annotation of adaptive immune receptor repertoire dataset generated by high-throughput sequencing (AIRRseq) and facilitate comparisons between studies and species. Interestingly, the assembly of salmonids IGH genomic sequences is challenging due to the presence of two large size duplicated IGH loci and high numbers of IG genes and pseudogenes. We used data available for Atlantic salmon to establish an IMGT standardized nomenclature of IGH genes in this species and then applied the IMGT rules to the rainbow trout IGH loci to set up a nomenclature, which takes into account the specificities of Salmonid loci. This unique, consistent nomenclature for Salmonid IGH genes was then used to construct IMGT sequence reference directories allowing accurate annotation of AIRRseq data. The complex issues raised by the genetic diversity of salmon and trout strains are discussed in the context of IG repertoire annotation.

Anadromous Arctic Char Microbiomes: Bioprospecting in the High Arctic

Northern populations of Arctic char (Salvelinus alpinus) can be anadromous, migrating annually from the ocean to freshwater lakes and rivers in order to escape sub-zero temperatures. Such seasonal behavior demands that these fish and their associated microbiomes adapt to changes in salinity, temperature, and other environmental challenges. We characterized the microbial community composition of anadromous S. alpinus, netted by Inuit fishermen at freshwater and seawater fishing sites in the high Arctic, both under ice and in open water. Bacterial profiles were generated by DNA extraction and high-throughput sequencing of PCR-amplified 16S ribosomal RNA genes. Results showed that microbial communities on the skin and intestine of Arctic char were statistically different when sampled from freshwater or saline water sites. This association was tested using hierarchical Ward's linkage clustering, showing eight distinct clusters in each of the skin and intestinal microbiomes, with the clusters reflecting sampling location between fresh and saline environments, confirming a salinity-linked turnover. This analysis also provided evidence for a core composition of skin and intestinal bacteria, with the phyla Proteobacteria, Firmicutes, and Cyanobacteria presenting as major phyla within the skin-associated microbiomes. The intestine-associated microbiome was characterized by unidentified genera from families Fusobacteriaceae, Comamonadaceae, Pseudomonadaceae, and Vibrionaceae. The salinity-linked turnover was further tested through ordinations that showed samples grouping based on environment for both skin- and intestine-associated microbiomes. This finding implies that core microbiomes between fresh and saline conditions could be used to assist in regulating optimal fish health in aquaculture practices. Furthermore, identified taxa from known psychrophiles and with nitrogen cycling properties suggest that there is additional potential for biotechnological applications for fish farm and waste management practices.

Gradual and Acute Temperature Rise Induces Crossing Endocrine, Metabolic, and Immunological Pathways in Maraena Whitefish (Coregonus maraena)

The complex and still poorly understood nature of thermoregulation in various fish species complicates the determination of the physiological status on the basis of diagnostic marker genes and indicative molecular pathways. The present study aimed to compare the physiological impacts of both gradual and acute temperature rise from 18 to 24°C on maraena whitefish in aquaculture. Microarray-based transcriptome profiles in the liver, spleen and kidney of heat-stressed maraena whitefish revealed the modulation of a significantly higher number of genes in those groups exposed to gradually rising temperatures compared with the acutely stressed groups, which might reflect early adaptation mechanisms. Moreover, we suggest a common set of 11 differentially expressed genes that indicate thermal stress induced by gradual or acute temperature rise in the three selected tissues. Besides the two pathways regulated in both data sets unfolded protein response and aldosterone signaling in epithelial cells, we identified unique tissue- and stress type-specific pathways reflecting the crossroads between signal transduction, metabolic and immunologic pathways to cope with thermal stress. In addition, comparing lists of differentially regulated genes with meta-analyzed published data sets revealed that “acute temperature rise”-responding genes that encode members of the HSP70, HSP90, and HSP40 families; their functional homologs; co-chaperones and stress-signal transducers are well-conserved across different species, tissues and/or cell types and experimental approaches.

Transmission routes maintaining a viral pathogen of steelhead trout within a complex multi-host assemblage

This is the first comprehensive region wide, spatially explicit epidemiologic analysis of surveillance data of the aquatic viral pathogen infectious hematopoietic necrosis virus (IHNV) infecting native salmonid fish. The pathogen has been documented in the freshwater ecosystem of the Pacific Northwest of North America since the 1950s, and the current report describes the disease ecology of IHNV during 2000-2012. Prevalence of IHNV infection in monitored salmonid host cohorts ranged from 8% to 30%, with the highest levels observed in juvenile steelhead trout. The spatial distribution of all IHNV-infected cohorts was concentrated in two sub-regions of the study area, where historic burden of the viral disease has been high. During the study period, prevalence levels fluctuated with a temporal peak in 2002. Virologic and genetic surveillance data were analyzed for evidence of three separate but not mutually exclusive transmission routes hypothesized to be maintaining IHNV in the freshwater ecosystem. Transmission between year classes of juvenile fish at individual sites (route 1) was supported at varying levels of certainty in 10%-55% of candidate cases, transmission between neighboring juvenile cohorts (route 2) was supported in 31%-78% of candidate cases, and transmission from adult fish returning to the same site as an infected juvenile cohort was supported in 26%-74% of candidate cases. The results of this study indicate that multiple specific transmission routes are acting to maintain IHNV in juvenile fish, providing concrete evidence that can be used to improve resource management. Furthermore, these results demonstrate that more sophisticated analysis of available spatio-temporal and genetic data is likely to yield greater insight in future studies.

Adverse Husbandry of Maraena Whitefish Directs the Immune System to Increase Mobilization of Myeloid Cells and Proinflammatory Responses

Adverse life circumstances evoke a common “conserved transcriptional response to adversity” (CTRA) in mammalian leukocytes. To investigate whether this pattern is preserved in lower vertebrates, maraena whitefish (Coregonus maraena) were exposed for 9 days to different stocking densities: ~10 kg/m³ (low density), ~33 kg/m³ (moderate), ~60 kg/m³ (elevated), and ~100 kg/m³ (high). Transcriptome profiling in the liver and kidney of individuals from each group suggested that crowding conditions activate stress-related signaling and effector pathways. Remarkably, about one-quarter of the genes differentially expressed under crowding conditions were involved in the activation of immune pathways such as acute-phase response and interleukin/TNF signaling attended by the simultaneous reduction of antiviral potency. Network analysis confirmed the complex interdigitation of immune- and stress-relevant pathways with interleukin-1 playing a central role. Antibody-based techniques revealed remarkable changes in the blood composition of whitefish and demonstrated the correlation between increasing stocking densities and elevated number of myeloid cells together with the increased phagocytic activity of peripheral blood leukocytes. In line with current studies in mammals, we conclude that crowding stress triggers in whitefish hallmarks of a CTRA, indicating that the stress-induced molecular mechanisms regulating the immune responses not only are conserved within mammals but were established earlier in evolution.

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.

Migration costs drive convergence of threshold traits for migratory tactics

Partial migration of some, but not all, members of a population is a common form of migration. We evaluated how migration costs influence which members migrate in 10 populations of two salmonid species. The migratory patterns of both species were evaluated based on the size at maturity for resident males, which is the threshold trait that determines the migratory tactics used within a population. In both species, this size was smaller in males located further from the sea, where migration costs are presumably higher. Moreover, the threshold sizes at maturity in males were correlated between both species. Our results suggest that migration costs are a significant convergent selective force on migratory tactics and life-history traits in nature.

L-Kynurenine, an amino acid identified as a sex pheromone in the urine of ovulated female masu salmon

Many animals employ sex pheromones to find mating partners during their reproductive seasons. However, most sex pheromones of vertebrates remain to be identified. Over the past 20 years, steroids and prostaglandins have been identified as sex pheromones in several fishes. These pheromones are broadly termed "hormonal pheromones" because they or their precursors act as hormones in these fishes. Hitherto, no other type of sex pheromone has been unambiguously identified in teleost fish. Here we report the identification of a "nonhormonal pheromone" in teleost fish. The urine of the reproductively mature female masu salmon (Oncorhynchus masou) contains a male-attracting pheromone. Bioassay-guided fractionation yielded an active compound that was identical to L-kynurenine in spectral and chromatographic properties. L-Kynurenine is a major metabolite of L-tryptophan in vertebrates. This pheromone elicits a male-specific behavior at even picomolar concentrations; its electrophysiological threshold is 10(-14) M. L-Kynurenine is a reasonable substance for female masu salmon to advertise their readiness for mating.