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.

Common microgeographical selection patterns revealed in four European conifers.

Microgeographical adaptation occurs when the effects of directional selection persist despite gene flow. Traits and genetic loci under selection can then show adaptive divergence, against the backdrop of little differentiation at other traits or loci. How common such events are and how strong the selection is that underlies them remain open questions. Here, we discovered and analysed microgeographical patterns of genomic divergence in four European and Mediterranean conifers with widely differing life-history traits and ecological requirements (Abies alba MIll., Cedrus atlantica [Endl.] Manetti, Pinus halepensis Mill. and Pinus pinaster Aiton) by screening pairs from geographically close forest stands sampled along steep ecological gradients. We inferred patterns of genomic divergence by applying a combination of divergence outlier detection methods, demographic modelling, Approximate Bayesian Computation inferences and genomic annotation to genomic data. Surprisingly for such small geographical scales, we showed that selection is strong in all species but generally affects different loci in each. A clear signature of selection was systematically detected on a fraction of the genome, of the order of 0.1%-1% of the loci depending on the species. The novel modelling method we designed for estimating selection coefficients showed that the microgeographical selection coefficient scaled by population size (Ns) was 2-30. Our results convincingly suggest that selection maintains within-population diversity at microgeographical scales in spatially heterogeneous environments. Such genetic diversity is likely to be a major reservoir of adaptive potential, helping populations to adapt under fluctuating environmental conditions.

Inferring selection in instances of long-range colonization: The Aleppo pine (Pinus halepensis) in the Mediterranean Basin.

Teasing apart the effects of natural selection and demography on current allele frequencies is challenging, due to both processes leaving a similar molecular footprint. In particular, when attempting to identify selection in species that have undergone a recent range expansion, the increase in genetic drift at the edges of range expansions ("allele surfing") can be a confounding factor. To address this potential issue, we first assess the long-range colonization history of the Aleppo pine across the Mediterranean Basin, using molecular markers. We then look for single nucleotide polymorphisms (SNPs) involved in local adaptation using: (a) environmental correlation methods (bayenv2), focusing on bioclimatic variables important for the species' adaptation (i.e., temperature, precipitation and water availability); and (b) FST -related methods (pcadapt). To assess the rate of false positives caused by the allele surfing effect, these results are compared with results from simulated SNP data that mimics the species' past range expansions and the effect of genetic drift, but with no selection. We find that the Aleppo pine shows a previously unsuspected complex genetic structure across its range, as well as evidence of selection acting on SNPs involved with the response to bioclimatic variables such as drought. This study uses an original approach to disentangle the confounding effects of drift and selection in range margin populations. It also contributes to the increased evidence that plant populations are able to adapt to new environments despite the expected accumulation of deleterious mutations that takes place during long-range colonizations.

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.

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.

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.

Looking for Local Adaptation: Convergent Microevolution in Aleppo Pine (Pinus halepensis).

Finding outlier loci underlying local adaptation is challenging and is best approached by suitable sampling design and rigorous method selection. In this study, we aimed to detect outlier loci (single nucleotide polymorphisms, SNPs) at the local scale by using Aleppo pine (Pinus halepensis), a drought resistant conifer that has colonized many habitats in the Mediterranean Basin, as the model species. We used a nested sampling approach that considered replicated altitudinal gradients for three contrasting sites. We genotyped samples at 294 SNPs located in genomic regions selected to maximize outlier detection. We then applied three different statistical methodologies-Two Bayesian outlier methods and one latent factor principal component method-To identify outlier loci. No SNP was an outlier for all three methods, while eight SNPs were detected by at least two methods and 17 were detected only by one method. From the intersection of outlier SNPs, only one presented an allelic frequency pattern associated with the elevational gradient across the three sites. In a context of multiple populations under similar selective pressures, our results underline the need for careful examination of outliers detected in genomic scans before considering them as candidates for convergent adaptation.