An Illumina approach to MHC typing of Atlantic salmon

Parallel Selection in Domesticated Atlantic Salmon from Divergent Founders Including on Whole-Genome Duplication-derived Homeologous Regions

Domestication and artificial selection for desirable traits have driven significant phenotypic changes and left detectable genomic footprints in farmed animals. Since the 1960s, intensive breeding has led to the rapid domestication of Atlantic salmon (Salmo salar), with multiple independent events that make it a valuable model for studying early domestication stages and the parallel evolution of populations of different origins subjected to similar selection pressures. Some aquatic species, including Atlantic salmon, have undergone whole-genome duplication (WGD), raising the possibility that genetic redundancy resulting from WGD has contributed to adaptation in captive environments, as seen in plants. Here, we examined the genomic responses to domestication in Atlantic salmon, focusing on potential signatures of parallel selection, including those associated with WGD. Candidate genomic regions under selection were identified by comparing whole-genome sequences from aquaculture and wild populations across 2 independently domesticated lineages (Western Norway and North America) using a genome-wide scan that combined 3 statistical methods: allele frequencies (FST), site frequency (Tajima's D), and haplotype differentiation (XP-EHH). These analyses revealed shared selective sweeps on identical SNPs in major histocompatibility complex (MHC) genes across aquaculture populations. This suggests that a combination of long-term balancing selection and recent human-induced selection has shaped MHC gene evolution in domesticated salmon. Additionally, we observed selective sweeps on a small number of gene pairs in homeologous regions originating from WGD, offering insights into how historical genome duplication events may intersect with recent selection pressures in aquaculture species.

High-Resolution Genotyping of Expressed Equine MHC Reveals a Highly Complex MHC Structure

The Major Histocompatibility Complex (MHC) genes play a key role in a number of biological processes, most notably in immunological responses. The MHCI and MHCII genes incorporate a complex set of highly polymorphic and polygenic series of genes, which, due to the technical limitations of previously available technologies, have only been partially characterized in non-model but economically important species such as the horse. The advent of high-throughput sequencing platforms has provided new opportunities to develop methods to generate high-resolution sequencing data on a large scale and apply them to the analysis of complex gene sets such as the MHC. In this study, we developed and applied a MiSeq-based approach for the combined analysis of the expressed MHCI and MHCII repertoires in cohorts of Thoroughbred, Icelandic, and Norwegian Fjord Horses. The approach enabled us to generate comprehensive MHCI/II data for all of the individuals (n = 168) included in the study, identifying 152 and 117 novel MHCI and MHCII sequences, respectively. There was limited overlap in MHCI and MHCII haplotypes between the Thoroughbred and the Icelandic/Norwegian Fjord horses, showcasing the variation in MHC repertoire between genetically divergent breeds, and it can be inferred that there is much more MHC diversity in the global horse population. This study provided novel insights into the structure of the expressed equine MHC repertoire and highlighted unique features of the MHC in horses.

Functional immune diversity in reindeer reveals a high Arctic population at risk

Climate changes the geographic range of both species as well as pathogens, causing a potential increase in the vulnerability of populations or species with limited genetic diversity. With advances in high throughput sequencing (HTS) technologies, we can now define functional expressed genetic diversity of wild species at a larger scale and identify populations at risk. Previous studies have used genomic DNA to define major histocompatibility complex (MHC) class II diversity in reindeer. Varying numbers of expressed genes found in many ungulates strongly argues for using cDNA in MHC typing strategies to ensure that diversity estimates relate to functional genes. We have used available reindeer genomes to identify candidate genes and established an HTS approach to define expressed MHC class I and class II diversity. To capture a broad diversity we included samples from wild reindeer from Southern Norway, semi-domesticated reindeer from Northern Norway and reindeer from the high Artic archipelago Svalbard. Our data show a medium MHC diversity in semi-domesticated and wild Norwegian mainland reindeer, and low MHC diversity reindeer in Svalbard reindeer. The low immune diversity in Svalbard reindeer provides a potential risk if the pathogenic pressure changes in response to altered environmental conditions due to climate change, or increased human-related activity.

A pan-genome data structure induced by pooled sequencing facilitates variant mining in heterogeneous germplasm

Valuable genetic variation lies unused in gene banks due to the difficulty of exploiting heterogeneous germplasm accessions. Advances in molecular breeding, including transgenics and genome editing, present the opportunity to exploit hidden sequence variation directly. Here we describe the pan-genome data structure induced by whole-genome sequencing of pooled individuals from wild populations of Patellifolia spp., a source of disease resistance genes for the related crop species sugar beet (Beta vulgaris). We represent the pan-genome as a map of reads from pooled sequencing of a heterogeneous population sample to a reference genome, plus a BLAST data base of the mapped reads. We show that this basic data structure can be queried by reference genome position or homology to identify sequence variants present in the wild relative, at genes of agronomic interest in the crop, a process known as allele or variant mining. Further we demonstrate the possibility of cataloging variants in all Patellifolia genomic regions that have corresponding single copy orthologous regions in sugar beet. The data structure, termed a "pooled read archive," can be produced, altered, and queried using standard tools to facilitate discovery of agronomically-important sequence variation.

Supplementary information

The online version contains supplementary material available at 10.1007/s11032-022-01308-6.

Tetraploid Ancestry Provided Atlantic Salmon With Two Paralogue Functional T Cell Receptor Beta Regions Whereof One Is Completely Novel

Protective cellular immune responses have been difficult to study in fish, due to lack of basic understanding of their T cell populations, and tools to study them. Cellular immunity is thus mostly ignored in vaccination and infection studies compared to humoral responses. High throughput sequencing, as well as access to well assembled genomes, now advances studies of cellular responses. Here we have used such resources to describe organization of T cell receptor beta genes in Atlantic salmon. Salmonids experienced a unique whole genome duplication approximately 94 million years ago, which provided these species with many functional duplicate genes, where some duplicates have evolved new functions or sub-functions of the original gene copy. This is also the case for T cell receptor beta, where Atlantic salmon has retained two paralogue T cell receptor beta regions on chromosomes 01 and 09. Compared to catfish and zebrafish, the genomic organization in both regions is unique, each chromosomal region organized with dual variable- diversity- joining- constant genes in a head to head orientation. Sequence identity of the chromosomal constant sequences between TRB01 and TRB09 is suggestive of rapid diversification, with only 67 percent as opposed to the average 82-90 percent for other duplicated genes. Using virus challenged samples we find both regions expressing bona fide functional T cell receptor beta molecules. Adding the 292 variable T cell receptor alpha genes to the 100 variable TRB genes from 14 subgroups, Atlantic salmon has one of the most diverse T cell receptor alpha beta repertoire of any vertebrate studied so far. Perhaps salmonid cellular immunity is more advanced than we have imagined.

DNA metabarcoding across disciplines: sequencing our way to greater understanding across scales of biological organization

DNA metabarcoding describes the use of targeted DNA (i.e., amplicon) sequencing to identify community constituents from a complex sample containing genetic material from multiple organisms, such as water, soil, gut contents, microbiomes, or biofilms. This molecular approach for characterizing mixed DNA samples relies on the development of "universal primers" that allow for effective amplification of target sequences across a broad range of taxa. Armed with optimized lab protocols and rigorous bioinformatics tools, DNA metabarcoding can produce a wealth of information about the hidden biodiversity of various sample types by probing for organisms' molecular footprints. DNA metabarcoding has received considerable popular press over the last few years because of gut microbiome studies in humans and beyond. However, there are many other applications that are continually integrating molecular biology with other fields of study to address questions that have previously been unanswerable, for both prokaryotic and eukaryotic targets. For example, we can now sample mostly-digested gut contents from virtually any organism to learn about ontogeny and foraging ecology. Water samples collected from different locations can be filtered to extract eDNA (i.e., environmental DNA), revealing the biodiversity of fishes and other taxa targeted by carefully selected primer sets. This universal primer metabarcoding approach has even been extended to looking at diverse gene families within single species, which is particularly useful for complex immune system genetics. The purpose of this SICB symposium was to bring together researchers using DNA metabarcoding approaches to (a) showcase the diversity of applications of this technique for addressing questions spanning ecology, evolution, and physiology, and (b) to spark connections among investigators from different fields that are utilizing similar approaches to facilitate optimization and standardization of metabarcoding methods and analyses. The resulting manuscripts from this symposium represent a great diversity of metabarcoding applications and taxonomic groups of interest.

Highly-contiguous bovine genomes underpin accurate functional analyses and updated nomenclature of MHC class I

The major histocompatibility complex (MHC) class I region of cattle is both highly polymorphic and, unlike many species, highly variable in gene content between haplotypes. Cattle MHC class I alleles were historically grouped by sequence similarity in the more conserved 3' end of the coding sequence to form phylogenetic allele groups. This has formed the basis of current cattle MHC class I nomenclature. We presently describe and compare five fully assembled MHC class I haplotypes using the latest cattle and yak genome assemblies. Of the five previously described "pseudogenes" in the cattle MHC class I region, Pseudogene 3 is putatively functional in all haplotypes and Pseudogene 6 and Pseudogene 7 are putatively functional in some haplotypes. This was reinforced by evidence of transcription. Based on full gene sequences as well as 3' coding sequence, we identified distinct subgroups of BoLA-3 and BoLA-6 that represent distinct genetic loci. We further examined allele-specific expression using transcriptomic data revealing that certain alleles are consistently weakly expressed compared to others. These observations will help to inform further studies into how MHC class I region variability influences T cell and natural killer cell functions in cattle.