Paucity of class I MHC gene heterogeneity between individuals in the endangered Hawaiian monk seal population

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.

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.

Effect of the enzyme and PCR conditions on the quality of high-throughput DNA sequencing results

Library preparation protocols for high-throughput DNA sequencing (HTS) include amplification steps in which errors can build up. In order to have confidence in the sequencing data, it is important to understand the effects of different Taq polymerases and PCR amplification protocols on the DNA molecules sequenced. We compared thirteen enzymes in three different marker systems: simple, single copy nuclear gene and complex multi-gene family. We also tested a modified PCR protocol, which has been suggested to reduce errors associated with amplification steps. We find that enzyme choice has a large impact on the proportion of correct sequences recovered. The most complex marker systems yielded fewer correct reads, and the proportion of correct reads was greatly affected by the enzyme used. Modified cycling conditions did reduce the number of incorrect sequences obtained in some cases, but enzyme had a much greater impact on the number of correct reads. Thus, the coverage required for the safe identification of genotypes using one of the low quality enzymes could be seven times larger than with more efficient enzymes in a biallelic system with equal amplification of the two alleles. Consequently, enzyme selection for downstream HTS has important consequences, especially in complex genetic systems.

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

Pinnipeds, marine carnivores, diverged from terrestrial carnivores ~45 million years ago, before their adaptation to marine environments. This lifestyle change exposed pinnipeds to different microbiota and pathogens, with probable impact on their MHC class I genes. Investigating this question, genomic sequences were determined for 71 MHC class I variants: 27 from harbor seal and 44 from gray seal. These variants form three MHC class I gene lineages, one comprising a pseudogene. The second, a candidate nonclassical MHC class I gene, comprises a nonpolymorphic transcribed gene related to dog DLA-79 and giant panda Aime-1906. The third is the diversity lineage, which includes 62 of the 71 seal MHC class I variants. All are transcribed, and they minimally represent six harbor and 12 gray seal MHC class I genes. Besides species-specific differences in gene number, seal MHC class I haplotypes exhibit gene content variation and allelic polymorphism. Patterns of sequence variation, and of positions for positively selected sites, indicate the diversity lineage genes are the seals' classical MHC class I genes. Evidence that expansion of diversity lineage genes began before gray and harbor seals diverged is the presence in both species of two distinctive sublineages of diversity lineage genes. Pointing to further expansion following the divergence are the presence of species-specific genes and greater MHC class I diversity in gray seals than harbor seals. The elaboration of a complex variable family of classical MHC class I genes in pinnipeds contrasts with the single, highly polymorphic classical MHC class I gene of dog and giant panda, terrestrial carnivores.

Evolution of MHC class I genes in the European badger (Meles meles)

The major histocompatibility complex (MHC) plays a central role in the adaptive immune system and provides a good model with which to understand the evolutionary processes underlying functional genes. Trans-species polymorphism and orthology are both commonly found in MHC genes; however, mammalian MHC class I genes tend to cluster by species. Concerted evolution has the potential to homogenize different loci, whereas birth-and-death evolution can lead to the loss of orthologs; both processes result in monophyletic groups within species. Studies investigating the evolution of MHC class I genes have been biased toward a few particular taxa and model species. We present the first study of MHC class I genes in a species from the superfamily Musteloidea. The European badger (Meles meles) exhibits moderate variation in MHC class I sequences when compared to other carnivores. We identified seven putatively functional sequences and nine pseudogenes from genomic (gDNA) and complementary (cDNA) DNA, signifying at least two functional class I loci. We found evidence for separate evolutionary histories of the α1 and α2/α3 domains. In the α1 domain, several sequences from different species were more closely related to each other than to sequences from the same species, resembling orthology or trans-species polymorphism. Balancing selection and probable recombination maintain genetic diversity in the α1 domain, evidenced by the detection of positive selection and a recombination event. By comparison, two recombination breakpoints indicate that the α2/α3 domains have most likely undergone concerted evolution, where recombination has homogenized the α2/α3 domains between genes, leading to species-specific clusters of sequences. Our findings highlight the importance of analyzing MHC domains separately.

Major Histocompatibility Complex (MHC) markers in conservation biology

Human impacts through habitat destruction, introduction of invasive species and climate change are increasing the number of species threatened with extinction. Decreases in population size simultaneously lead to reductions in genetic diversity, ultimately reducing the ability of populations to adapt to a changing environment. In this way, loss of genetic polymorphism is linked with extinction risk. Recent advances in sequencing technologies mean that obtaining measures of genetic diversity at functionally important genes is within reach for conservation programs. A key region of the genome that should be targeted for population genetic studies is the Major Histocompatibility Complex (MHC). MHC genes, found in all jawed vertebrates, are the most polymorphic genes in vertebrate genomes. They play key roles in immune function via immune-recognition and -surveillance and host-parasite interaction. Therefore, measuring levels of polymorphism at these genes can provide indirect measures of the immunological fitness of populations. The MHC has also been linked with mate-choice and pregnancy outcomes and has application for improving mating success in captive breeding programs. The recent discovery that genetic diversity at MHC genes may protect against the spread of contagious cancers provides an added impetus for managing and protecting MHC diversity in wild populations. Here we review the field and focus on the successful applications of MHC-typing for conservation management. We emphasize the importance of using MHC markers when planning and executing wildlife rescue and conservation programs but stress that this should not be done to the detriment of genome-wide diversity.

Range-wide genetic connectivity of the Hawaiian monk seal and implications for translocation

The Hawaiian monk seal (Monachus schauinslandi) is one of the most critically endangered marine mammals. Less than 1200 individuals remain, and the species is declining at a rate of approximately 4% per year as a result of juvenile starvation, shark predation, and entanglement in marine debris. Some of these problems may be alleviated by translocation; however, if island breeding aggregates are effectively isolated subpopulations, moving individuals may disrupt local adaptations. In these circumstances, managers must balance the pragmatic need of increasing survival with theoretical concerns about genetic viability. To assess range-wide population structure of the Hawaiian monk seal, we examined an unprecedented, near-complete genetic inventory of the species (n =1897 seals, sampled over 14 years) at 18 microsatellite loci. Genetic variation was not spatially partitioned ((w) =-0.03, p = 1.0), and a Bayesian clustering method provided evidence of one panmictic population (K =1). Pairwise F(ST) comparisons (among 7 island aggregates over 14 annual cohorts) did not reveal temporally stable, spatial reproductive isolation. Our results coupled with long-term tag-resight data confirm seal movement and gene flow throughout the Hawaiian Archipelago. Thus, human-mediated translocation of seals among locations is not likely to result in genetic incompatibilities.

Evolution and survival of marine carnivores did not require a diversity of killer cell Ig-like receptors or Ly49 NK cell receptors

Ly49 lectin-like receptors and killer cell Ig-like receptors (KIR) are structurally unrelated cell surface glycoproteins that evolved independently to function as diverse NK cell receptors for MHC class I molecules. Comparison of primates and various domesticated animals has shown that species have either a diverse Ly49 or KIR gene family, but not both. In four pinniped species of wild marine carnivore, three seals and one sea lion, we find that Ly49 and KIR are each represented by single, orthologous genes that exhibit little polymorphism and are transcribed to express cell surface protein. Pinnipeds are therefore species in which neither Ly49 nor KIR are polygenic, but retain the ancestral single-copy state. Whereas pinniped Ly49 has been subject to purifying selection, we find evidence for positive selection on KIR3DL during pinniped evolution. This selection, which focused on the D0 domain and the stem, points to the functionality of the KIR and most likely led to the sea lion's loss of D0. In contrast to the dynamic and rapid evolution of the KIR and Ly49 genes in other species, the pinniped KIR and Ly49 have been remarkably stable during the >33 million years since the last common ancestor of seals and sea lions. These results demonstrate that long-term survival of placental mammal species need not require a diverse system of either Ly49 or KIR NK cell receptors.

Extremely low genetic diversity in the endangered Hawaiian monk seal (Monachus schauinslandi)

Hunted to near extinction in the late 19th century, the endangered and endemic Hawaiian monk seal (Monachus schauinslandi) exhibits low variation at all molecular markers tested to date. Here we confirm extreme paucity of genetic diversity, finding polymorphisms at only 8 of 154 microsatellite loci tested (143 novel species-specific loci, 10 loci from Antarctic seals, and 1 previously characterized locus). This screening revealed unprecedentedly low levels of allelic diversity and heterozygosity (A = 1.1, H(e) = 0.026). Subsequent analyses of 2409 Hawaiian monk seals at the 8 polymorphic loci provide evidence for a bottleneck (P = 0.002), but simulations indicate low genetic diversity (H(e) < 0.09) prior to recorded human influence. There is little indication of contemporary inbreeding (F(IS) = 0.018) or population structure (K = 1 population). Minimal genetic variation did not prevent partial recovery by the late 1950s and may not be driving the current population decline to approximately 1200 seals. Nonetheless, genotyping nearly every individual living during the past 25 years sets a new benchmark for low genetic diversity in an endangered species.

Transmission of a fatal clonal tumor by biting occurs due to depleted MHC diversity in a threatened carnivorous marsupial

A fatal transmissible tumor spread between individuals by biting has emerged in the Tasmanian devil (Sarcophilus harrisii), a carnivorous marsupial. Here we provide genetic evidence establishing that the tumor is clonal and therefore foreign to host devils. Thus, the disease is highly unusual because it is not just a tumor but also a tissue graft, passed between individuals without invoking an immune response. The MHC plays a key role in immune responses to both tumors and grafts. The most common mechanism of immune evasion by tumors is down-regulation of classical cell surface MHC molecules. Here we show that this mode of immune escape does not occur. However, because the tumor is a graft, it should still be recognized and rejected by the host's immune system due to foreign cell surface antigens. Mixed lymphocyte responses showed a lack of alloreactivity between lymphocytes of different individuals in the affected population, indicating a paucity of MHC diversity. This result was verified by genotyping, providing a conclusive link between a loss of MHC diversity and spread of a disease through a wild population. This novel disease arose as a direct result of loss of genetic diversity and the aggressive behavior of the host species. The neoplastic clone continues to spread although the population, and, without active disease control by removal of affected animals and the isolation of disease-free animals, the Tasmanian devil faces extinction.