The evolution of the natural killer complex; a comparison between mammals using new high-quality genome assemblies and targeted annotation.

Natural killer (NK) cells are a diverse population of lymphocytes with a range of biological roles including essential immune functions. NK cell diversity is in part created by the differential expression of cell surface receptors which modulate activation and function, including multiple subfamilies of C-type lectin receptors encoded within the NK complex (NKC). Little is known about the gene content of the NKC beyond rodent and primate lineages, other than it appears to be extremely variable between mammalian groups. We compared the NKC structure between mammalian species using new high-quality draft genome assemblies for cattle and goat; re-annotated sheep, pig, and horse genome assemblies; and the published human, rat, and mouse lemur NKC. The major NKC genes are largely in the equivalent positions in all eight species, with significant independent expansions and deletions between species, allowing us to propose a model for NKC evolution during mammalian radiation. The ruminant species, cattle and goats, have independently evolved a second KLRC locus flanked by KLRA and KLRJ, and a novel KLRH-like gene has acquired an activating tail. This novel gene has duplicated several times within cattle, while other activating receptor genes have been selectively disrupted. Targeted genome enrichment in cattle identified varying levels of allelic polymorphism between the NKC genes concentrated in the predicted extracellular ligand-binding domains. This novel recombination and allelic polymorphism is consistent with NKC evolution under balancing selection, suggesting that this diversity influences individual immune responses and may impact on differential outcomes of pathogen infection and vaccination

A pan-cetacean MHC amplicon sequencing panel developed and evaluated in combination with genome assemblies

This study was funded by a Royal Society Research Grants for Research Fellows (RGF\R1\181014) to E.C.G. E.C.G. is funded by a Royal Society University Research Fellowship (UF160081 & URF\R\221020). F.E. is supported by a University of St Andrews School of Biology Ph.D. scholarship and a Royal Society Research Fellows Enhancement Award (RGF\EA\180213 to E.C.G). E.L.C. is funded by a Rutherford Discovery Fellowship from the Royal Society of New Zealand Te Apārangi.The major histocompatibility complex (MHC) is a highly polymorphic gene family that is crucial in immunity, and its diversity can be effectively used as a fitness marker for populations. Despite this, MHC remains poorly characterised in non-model species (e.g., cetaceans: whales, dolphins and porpoises) as high gene copy number variation, especially in the fast-evolving class I region, makes analyses of genomic sequences difficult. To date, only small sections of class I and IIa genes have been used to assess functional diversity in cetacean populations. Here, we undertook a systematic characterisation of the MHC class I and IIa regions in available cetacean genomes. We extracted full-length gene sequences to design pan-cetacean primers that amplified the complete exon2 from MHC class I and IIa genes in one combined sequencing panel. We validated this panel in 19 cetacean species and described 354 alleles for both classes. Furthermore, we identified likely assembly artefacts for many MHC class I assemblies based on the presence of class I genes in the amplicon data compared to missing genes from genomes. Finally, we investigated MHC diversity using the panel in 25 humpback and 30 southern right whales, including four paternity trios for humpback whales. This revealed copy-number variable class I haplotypes in humpback whales, which is likely a common phenomenon across cetaceans. These MHC alleles will form the basis for a cetacean branch of the Immuno-Polymorphism Database (IPD-MHC), a curated resource intended to aid in the systematic compilation of MHC alleles across several species, to support conservation initiatives.Peer reviewe

Hector's and Maui's dolphin MHC and microsatellite genotyping data

The file contains MHC genotyping data from SSCP and microsatellite data for 162 individual Hector's and Maui's dolphins from New Zealand. The individuals are from two subspecies and two populations of the Hector's dolphins

Data from: The influence of selection on MHC DQA and DQB haplotypes in the endemic New Zealand Hector’s and Māui dolphins

Strong balancing selection on the Major Histocompatibility Complex (MHC) can lead to different patterns in gene frequencies and neutral genomic variation within species. We investigated diversity and geographic structure of MHC genes DQA and DQB, as well as their inferred functional haplotypes, from two regional populations (East and West Coast) of the endangered Hector’s dolphin (Cephalorhynchus hectori hectori) and the critically endangered Māui dolphin (Cephalorhynchus hectori maui) (West Coast, North Island), and contrasted these results with patterns from neutral microsatellites. The Māui had the lowest number of alleles for DQA (2) and DQB (3), consistent with strong genetic drift acting on this remnant population. However, the two retained DQA alleles are amongst the most divergent combinations of all four alleles found across the Hector’s metapopulation, potentially reflecting the retention of divergent alleles due to balancing selection. The high frequency of the divergent DQB04 allele also gave this population the highest nucleotide diversity for DQB. Strong differentiation was evident for DQA, DQB, and DQA-DQB haplotypes between the regional populations of Hector’s dolphins (FST>0.213) and both subspecies (FST>0.311). Differentiation was generally greater than observed at neutral microsatellite loci, suggesting the influence of selection between geographically proximate East and West Coast populations. This might be the result of spatial differences in directional selection on those opposite coastlines. Additionally, measures of the ratio of non-synonymous to synonymous substitutions (dN/dS) were consistent with balancing selection over evolutionary time. Together, these results suggest a complex interplay of balancing selection, directional selection, local fidelity, and genetic drift

A pan-cetacean MHC amplicon sequencing panel developed and evaluated in combination with genome assemblies

The major histocompatibility complex (MHC) is a highly polymorphic gene family that is crucial in immunity, and its diversity can be effectively used as a fitness marker for populations. Despite this, MHC remains poorly characterised in non-model species (e.g., cetaceans: whales, dolphins and porpoises) as high gene copy number variation, especially in the fast-evolving class I region, makes analyses of genomic sequences difficult. To date, only small sections of class I and IIa genes have been used to assess functional diversity in cetacean populations. Here, we undertook a systematic characterisation of the MHC class I and IIa regions in available cetacean genomes. We extracted full-length gene sequences to design pan-cetacean primers that amplified the complete exon2 from MHC class I and IIa genes in one combined sequencing panel. We validated this panel in 19 cetacean species and described 354 alleles for both classes. Furthermore, we identified likely assembly artefacts for many MHC class I assemblies based on the presence of class I genes in the amplicon data compared to missing genes from genomes. Finally, we investigated MHC diversity using the panel in 25 humpback and 30 southern right whales, including four paternity trios for humpback whales. This revealed copy-number variable class I haplotypes in humpback whales, which is likely a common phenomenon across cetaceans. These MHC alleles will form the basis for a cetacean branch of the Immuno-Polymorphism Database (IPD-MHC), a curated resource intended to aid in the systematic compilation of MHC alleles across several species, to support conservation initiatives

The evolution of the natural killer complex; a comparison between mammals using new high-quality genome assemblies and targeted annotation

Natural killer (NK) cells are a diverse population of lymphocytes with a range of biological roles including essential immune functions. NK cell diversity is in part created by the differential expression of cell surface receptors which modulate activation and function, including multiple subfamilies of C-type lectin receptors encoded within the NK complex (NKC). Little is known about the gene content of the NKC beyond rodent and primate lineages, other than it appears to be extremely variable between mammalian groups. We compared the NKC structure between mammalian species using new high-quality draft genome assemblies for cattle and goat; re-annotated sheep, pig, and horse genome assemblies; and the published human, rat, and mouse lemur NKC. The major NKC genes are largely in the equivalent positions in all eight species, with significant independent expansions and deletions between species, allowing us to propose a model for NKC evolution during mammalian radiation. The ruminant species, cattle and goats, have independently evolved a second KLRC locus flanked by KLRA and KLRJ, and a novel KLRH-like gene has acquired an activating tail. This novel gene has duplicated several times within cattle, while other activating receptor genes have been selectively disrupted. Targeted genome enrichment in cattle identified varying levels of allelic polymorphism between the NKC genes concentrated in the predicted extracellular ligand-binding domains. This novel recombination and allelic polymorphism is consistent with NKC evolution under balancing selection, suggesting that this diversity influences individual immune responses and may impact on differential outcomes of pathogen infection and vaccination. ELECTRONIC SUPPLEMENTARY MATERIAL: The online version of this article (doi:10.1007/s00251-017-0973-y) contains supplementary material, which is available to authorized users