Natural selection on marine carnivores elaborated a diverse family of classical MHC class I genes exhibiting haplotypic gene content variation and allelic polymorphism

Sequence diversity and differences at the highly duplicated MHC-I gene reflect viral susceptibility in sympatric pinniped species

Differences in disease susceptibility among species can result from rapid host-pathogen coevolution and differences in host species ecology that affect the strength and direction of natural selection. Among two sympatric pinniped species that differ in sociality and putative disease exposure, we investigate observed differences in susceptibility through an analysis of a highly variable, duplicated gene family involved in the vertebrate immune response. Using high-throughput amplicon sequencing, we characterize diversity at the two exons that encode the peptide binding region of the major histocompatibility complex class I (MHC-I) gene in harbor (N = 60) and gray (N = 90) seal populations from the Northwest Atlantic. Across species, we identified 106 full-length exon 2 and 103 exon 3 sequence variants and a minimum of 11 duplicated MHC-I loci. The sequence variants clustered in 15 supertypes defined by the physiochemical properties of the peptide binding region, including a putatively novel Northwest Atlantic MHC-I diversity sublineage. Trans-species polymorphisms, dN/dS ratios, and evidence of gene conversion among supertypes are consistent with balancing selection acting on this gene. High functional redundancy suggests particularly strong selection among gray seals at the novel Northwest Atlantic MHC-I diversity sublineage. At exon 2, harbor seals had a significantly greater number of variants per individual than gray seals, but fewer supertypes. Supertype richness and private supertypes are hypothesized to contribute to observed differences in disease resistance between species, as consistently, across the North Atlantic and many disease outbreaks, gray seals appear to be more resistant to respiratory viruses than harbor seals.

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.

Using PacBio SMRT data for identification of class I MHC alleles in a wildlife species, Zalophus californianus (California sea lion)

High allelic polymorphism and association with disease susceptibility has made the genes encoding major histocompatibility complex (MHC) antigen presentation molecules in humans, domesticated animals, and wildlife species of wide interest to ecologists, evolutionary biologists, and health specialists. The often multifaceted polygenism and extreme polymorphism of this immunogenetic system have made it especially difficult to characterize in non-model species. Here we compare and contrast the workflows of traditional Sanger sequencing of plasmid-cloned amplicons to Pacific Biosciences SMRT circular consensus sequencing (CCS) in their ability to capture alleles of MHC class I in a wildlife species where characterization of these genes was absent. We assessed two California sea lions (Zalophus californianus), a species suffering from a high prevalence of an aggressive cancer associated with a sexually transmitted gamma herpesvirus. In this pilot study, SMRT CCS proved superior in identifying more alleles from each animal than the more laborious plasmid cloning/Sanger workflow (12:7, 10:7), and no alleles were identified with the cloning/Sanger approach that were not identified by SMRT CCS. We discuss the advantages and disadvantages of each approach including cost, allele rarefaction, and sequence fidelity.

Evolution of major histocompatibility complex class I genes in the sable Martes zibellina (Carnivora, Mustelidae)

The molecules encoded by major histocompatibility complex (MHC) genes play an essential role in the adaptive immune response among vertebrates. We investigated the molecular evolution of MHC class I genes in the sable Martes zibellina. We isolated 26 MHC class I sequences, including 12 putatively functional sequences and 14 pseudogene sequences, from 24 individuals from two geographic areas of northeast China. The number of putatively functional sequences found in a single individual ranged from one to five, which might be at least 1-3 loci. We found that both balancing selection and recombination contribute to evolution of MHC class I genes in M. zibellina. In addition, we identified a candidate nonclassical MHC class I lineage in Carnivora, which may have preceded the divergence (about 52-57 Mya) of Caniformia and Feliformia. This may contribute to further understanding of the origin and evolution of nonclassical MHC class I genes. Our study provides important immune information of MHC for M. zibellina, as well as other carnivores.

Genomic signatures of population bottleneck and recovery in Northwest Atlantic pinnipeds

Population increases over the past several decades provide natural settings in which to study the evolutionary processes that occur during bottleneck, growth, and spatial expansion. We used parallel natural experiments of historical decline and subsequent recovery in two sympatric pinniped species in the Northwest Atlantic, the gray seal (Halichoerus grypus atlantica) and harbor seal (Phoca vitulina vitulina), to study the impact of recent demographic change in genomic diversity. Using restriction site-associated DNA sequencing, we assessed genomic diversity at over 8,700 polymorphic gray seal loci and 3,700 polymorphic harbor seal loci in samples from multiple cohorts collected throughout recovery over the past half-century. Despite significant differences in the degree of genetic diversity assessed in the two species, we found signatures of historical bottlenecks in the contemporary genomes of both gray and harbor seals. We evaluated temporal trends in diversity across cohorts, as well as compared samples from sites at both the center and edge of a recent gray seal range expansion, but found no significant change in genomewide diversity following recovery. We did, however, find that the variance and degree of allele frequency change measured over the past several decades were significantly different from neutral expectations of drift under population growth. These two cases of well-described demographic history provide opportunities for critical evaluation of current approaches to simulating and understanding the genetic effects of historical demographic change in natural populations.

The unique evolution of the pig LRC, a single KIR but expansion of LILR and a novel Ig receptor family

The leukocyte receptor complex (LRC) encodes numerous immunoglobulin (Ig)-like receptors involved in innate immunity. These include the killer-cell Ig-like receptors (KIR) and the leukocyte Ig-like receptors (LILR) which can be polymorphic and vary greatly in number between species. Using the recent long-read genome assembly, Sscrofa11.1, we have characterized the porcine LRC on chromosome 6. We identified a ~ 197-kb region containing numerous LILR genes that were missing in previous assemblies. Out of 17 such LILR genes and fragments, six encode functional proteins, of which three are inhibitory and three are activating, while the majority of pseudogenes had the potential to encode activating receptors. Elsewhere in the LRC, between FCAR and GP6, we identified a novel gene that encodes two Ig-like domains and a long inhibitory intracellular tail. Comparison with two other porcine assemblies revealed a second, nearly identical, non-functional gene encoding a short intracellular tail with ambiguous function. These novel genes were found in a diverse range of mammalian species, including a pseudogene in humans, and typically consist of a single long-tailed receptor and a variable number of short-tailed receptors. Using porcine transcriptome data, both the novel inhibitory gene and the LILR were highly expressed in peripheral blood, while the single KIR gene, KIR2DL1, was either very poorly expressed or not at all. These observations are a prerequisite for improved understanding of immune cell functions in the pig and other species.

Co-evolution of MHC class I and variable NK cell receptors in placental mammals

Shaping natural killer (NK) cell functions in human immunity and reproduction are diverse killer cell immunoglobulin-like receptors (KIRs) that recognize polymorphic MHC class I determinants. A survey of placental mammals suggests that KIRs serve as variable NK cell receptors only in certain primates and artiodactyls. Divergence of the functional and variable KIRs in primates and artiodactyls predates placental reproduction. Among artiodactyls, cattle but not pigs have diverse KIRs. Catarrhine (humans, apes, and Old World monkeys) and platyrrhine (New World monkeys) primates, but not prosimians, have diverse KIRs. Platyrrhine and catarrhine systems of KIR and MHC class I are highly diverged, but within the catarrhines, a stepwise co-evolution of MHC class I and KIR is discerned. In Old World monkeys, diversification focuses on MHC-A and MHC-B and their cognate lineage II KIR. With evolution of C1-bearing MHC-C from MHC-B, as informed by orangutan, the focus changes to MHC-C and its cognate lineage III KIR. Evolution of C2 from C1 and fixation of MHC-C drove further elaboration of MHC-C-specific KIR, as exemplified by chimpanzee. In humans, the evolutionary trajectory changes again. Emerging from reorganization of the KIR locus and selective attenuation of KIR avidity for MHC class I are the functionally distinctive KIR A and KIR B haplotypes.

Patterns of MHC-G-Like and MHC-B Diversification in New World Monkeys

The MHC class I (MHC-I) region in New World monkeys (Platyrrhini) has remained relatively understudied. To evaluate the diversification patterns and transcription behavior of MHC-I in Platyrrhini, we first analyzed public genomic sequences from the MHC-G-like subregion in Saimiri boliviensis, Ateles geoffroyi and Callicebus moloch, and from the MHC-B subregion in Saimiri boliviensis. While S. boliviensis showed multiple copies of both MHC-G-like (10) and –B (15) loci, A. geoffroyi and C. moloch had only three and four MHC-G-like genes, respectively, indicating that not all Platyrrhini species have expanded their MHC-I loci. We then sequenced MHC-G-like and -B cDNAs from nine Platyrrhini species, recovering two to five unique cDNAs per individual for both loci classes. In two Saguinus species, however, no MHC-B cDNAs were found. In phylogenetic trees, MHC-G-like cDNAs formed genus-specific clusters whereas the MHC-B cDNAs grouped by Platyrrhini families, suggesting a more rapid diversification of the former. Furthermore, cDNA sequencing in 12 capuchin monkeys showed that they transcribe at least four MHC-G-like and five MHC-B polymorphic genes, showing haplotypic diversity for gene copy number and signatures of positive natural selection at the peptide binding region. Finally, a quantitative index for MHC:KIR affinity was proposed and tested to predict putative interacting pairs. Altogether, our data indicate that i) MHC-I genes has expanded differentially among Platyrrhini species, ii) Callitrichinae (tamarins and marmosets) MHC-B loci have limited or tissue-specific expression, iii) MHC-G-like genes have diversified more rapidly than MHC-B genes, and iv) the MHC-I diversity is generated mainly by genetic polymorphism and gene copy number variation, likely promoted by natural selection for ligand binding.

Phocine distemper virus: current knowledge and future directions

Phocine distemper virus (PDV) was first recognized in 1988 following a massive epidemic in harbor and grey seals in north-western Europe. Since then, the epidemiology of infection in North Atlantic and Arctic pinnipeds has been investigated. In the western North Atlantic endemic infection in harp and grey seals predates the European epidemic, with relatively small, localized mortality events occurring primarily in harbor seals. By contrast, PDV seems not to have become established in European harbor seals following the 1988 epidemic and a second event of similar magnitude and extent occurred in 2002. PDV is a distinct species within the Morbillivirus genus with minor sequence variation between outbreaks over time. There is now mounting evidence of PDV-like viruses in the North Pacific/Western Arctic with serological and molecular evidence of infection in pinnipeds and sea otters. However, despite the absence of associated mortality in the region, there is concern that the virus may infect the large Pacific harbor seal and northern elephant seal populations or the endangered Hawaiian monk seals. Here, we review the current state of knowledge on PDV with particular focus on developments in diagnostics, pathogenesis, immune response, vaccine development, phylogenetics and modeling over the past 20 years.

Canine DLA-79 gene: an improved typing method, identification of new alleles and its role in graft rejection and graft-versus-host disease

Developing a preclinical canine model that predicts outcomes for hematopoietic cell transplantation in humans requires a model that mimics the degree of matching between human donor and recipient major histocompatibility complex (MHC) genes. The polymorphic class I and class II genes in mammals are typically located in a single chromosome as part of the MHC complex. However, a divergent class I gene in dogs, designated dog leukocyte antigen-79 (DLA-79), is located on chromosome 18 while other MHC genes are on chromosome 12. This gene is not taken into account while DLA matching for transplantation. Though divergent, this gene shares significant similarity in sequence and exon-intron architecture with other class I genes, and is transcribed. Little is known about the polymorphisms of DLA-79 and their potential role in transplantation. This study was aimed at exploring the reason for high rate of rejection seen in DLA-matched dogs given reduced intensity conditioning, in particular, the possibility that DLA-79 allele mismatches may be the cause. We found that about 82% of 407 dogs typed were homozygous for a single, reference allele. Owing to the high prevalence of a single allele, 87 of the 108 dogs (∼80%) transplanted were matched for DLA-79 with their donor. In conclusion, we have developed an efficient method to type alleles of a divergent MHC gene in dogs and identified two new alleles. We did not find any statistical correlation between DLA-79 allele disparity and graft rejection or graft-versus-host disease, among our transplant dogs.